(* H2+sp-Molecule: H2+ 0.5*(1s,1p) *)
func=0.5*9/(8*P^2)+0.5*9/(8*Q^2)-(1/(2*P))*(3-(R+Q)^2/P^2)-
2*0.25/(R+2*Q)-(1/(2*Q))*(3-R^2/Q^2)+
(0.5/Q)*(1-(3/20)*(P/Q)^2+(1/32)*(P/Q)^3)+
1/(2*(R+Q))

t=FindMinimum[func,{P,1.5},{Q,1},{R,0.2}]
u=t[[2]]
d=2*(Q+R)*0.529177 /. u
plot1=Graphics[{
Circle[{0,0},P],Circle[{-Q,0},Q],
Circle[{Q,0},Q],Disk[{Q+R,0},0.08],
Disk[{-(Q+R),0},0.08],
Line[{{-(Q+R),0},{Q+R,0}}]}] /. u

Show[plot1,{AspectRatio->Automatic}]