(* 
      Li2 g2at! no overlap! from Paramki 1998 *) 
c = {k0 -> 1.09926089, k1 -> 0.64854319, 
        sig0 -> 0.20631886, sig1 -> 0.18863924, 
        z -> 3.0}
ekin = 
        2.0*2.25*k0/P^2 + 
          2.25*k1/Q^2;              
r = P + Q;   p = 2.0*r;
vne = -6.0*z/P - 18.0/r;
vee = 
        2.0*3.0*sig0/P + 3.0*sig1/Q + 10.0/r;
vnn = z*z/p;
vir = -(vne + vnn + vee)/ekin  /.  
          c;             
etot = ekin + vne + vnn + vee  /.  c
N[t = FindMinimum[etot, {P, 0.58}, {Q, 1.94}], 
      20]
u = t[[2]];
N[vir  /.  u, 10]
N[p*0.529177  /.  u,  10]
N[(-t[[1]] - 2*7.4334)*27.2114, 
      10]
plot1 = 
      Graphics[{Circle[{-r, 0}, Q], Circle[{r, 0}, Q], Circle[{0, 0}, Q], 
            Disk[{-r, 0}, 0.08], Disk[{r, 0}, 0.08],  
            Line[{{-r, 0}, {r, 0}}]}]  /.  u
Show[
      plot1, {AspectRatio -> Automatic}]

(* 
      Li2 oo g2 no overlap fromParamki *) 
c1 = {k0 -> 1.00199347, k1 -> 1.07858954, sig0 -> 0.29788384, 
        sig1 -> 0.30644253, z -> 3.0}
ekin = 2.0*2.25*k0/P^2 + 9.0*k1/Q^2;
p = 2.0*R;
vne = -6.0*z/P - 6.0*(3.0 - (R/Q)^2)/Q - 
          6.0/R;
vee = 
        2.0*3.0*sig0/P + 3.0*sig1/Q + 
          4.0*(3.0 - (R/Q)^2 - 1.8*(P/Q)^2)/Q + 
          2.0/R;
vnn = z*z/p;
vir = -(vne + vnn + vee)/ekin  /.  
          c1;              
etot = ekin + vne + vnn + vee  /.  c1
N[t = FindMinimum[etot, {P, 0.5}, {Q, 5.0}, {R, 2.5}], 
      20]
u = t[[2]];
N[vir  /.  u, 10]
N[2*R*0.529177  /.  u,  10]
N[(-t[[1]] - 2*7.44047)*27.2114, 
      10]
plot1 = 
      Graphics[{Circle[{-R, 0}, P], Circle[{R, 0}, P], Circle[{0, 0}, Q], 
            Disk[{-R, 0}, 0.08], Disk[{R, 0}, 0.08],  
            Line[{{-R, 0}, {R, 0}}]}]  /.  u
Show[
      plot1, {AspectRatio -> 
          Automatic}]