expand Q*(delta(r+a)-delta(r-a)) at r = a