(* Li2 2 Li +  in Li - Li sphere;  
      Kimb_Loebl ext sphere not realizable! *) 
 (* no exchange! for molecule,  
      nor atoms *) 
Clear[c, T, Vne, Vee, Vnn, 
      u]
T = 4*2.25*k2/P^2 + 2*2.25*k1/Q^2;
Vne = -2*3*Z/Q - 
          2*Z/R - (2*
                Z/P)*(3 - (R/P)^2);
Vee = 
        4*(3 - (R/P)^2 - 0.6*(Q/P)^2)/P + 3*s2/P + 2*3*s1/Q + 
          2/R;
Vnn = Z*Z/(2*R);
func = T + Vne + Vee + Vnn;
  
c = {k1  ->  1.004,  k2 ->  1.0,  s1 -> 0.3, 
          s2 ->  0.3,  
          Z  ->  3.0};
func  =  func  /.  c
N[t = FindMinimum[func,  {P, 4.0}, {Q, 0.5  }, {R, 2.0}], 
      30]
u  =  t[[2]];
  
N[Vne  /.  u  /.  c,  10]
  
N[Vee  /.  u  /. c,  10]
  
N[Vnn  /.  u  /. c, 10]
  
N[-(Vne + Vee + Vnn)/T  /.  u  /.  c,  
      10]
  
N[2*R*0.529177  /.  u, 10]
  
N[(t[[1]] + 2*7.41587)*27.2114,  
      10]
plot1 = 
      Graphics[{Circle[{-R, 0}, Q], Circle[{R, 0}, Q], Circle[{0, 0}, P], 
            Disk[{-R, 0}, 0.08], Disk[{R, 0}, 0.08],  
            Line[{{-R, 0}, {R, 0}}]}]  /.  u
Show[plot1, {AspectRatio -> Automatic}]