r**2*(-4*Qa/d**4 + 2*U_md/d**2) + 4*Qa/d**2 - U_md