c = {k0 -> 0.94630662, 
        k1 ->  1.42576921, k2 ->  1.21505765, k3 ->  1.01808921, 
        sig0 ->  0.30001954, sig1 ->  0.40033732, sig2 ->  0.30034772, 
        z1 -> 9.0, z2 -> 3.0}
ekin = 
        2.25*k0/P^2 + 2.25*k1/Q^2 + 6.75*k2/R^2 + 
          2.25*k3/X^2;
vne = -3*z1/P - 3*z2/X;
vee = 3*sig0/P + 3*sig2/X;
asin = Sin[Y];
sin2 = asin^2;
acos = Sqrt[1 - sin2];
vee = vee + 3*sig1/Q + 9*sig1/R;
r = P + Q;
q = P + R;
o = Q + X;
p = r + o;
q2 = q*q;
r2 = r*r;
qr = q*r;
p2 = p*p;
vee = vee + 4/r + 4/p + 4/o;
vne = vne - 2*z1/r - 2*z1/p;
vne = vne - 2*z2/o - 2*z2/p;
vnn = z1*z2/p;
vee = vee + 12/q;
vne = vne - 6*z1/q;
vee = 
        vee + 12/Sqrt[2*qr*acos + q2 + r2];
vne = 
        vne - 6*z2/Sqrt[2*q*p*acos + p2 + q2];
vee = 
        vee + 12/Sqrt[2*q*p*acos + p2 + q2];
vee = vee + 12/(Sqrt[3]*q*asin);
vir = -(vne + vnn + vee)/ekin  /. 
          c;
etot = ekin + vne + vnn + vee  /.  c
eel = vne + ekin + vee;
N[t = 
        FindMinimum[
          etot, {P, 0.16}, {Q,  1.08}, {R,  0.85}, {X,  0.56}, {Y,  
            1.8}], 50]
u = t[[2]]
N[vir  /.  c  /.  u, 10]
N[(P + 2*Q + X)*0.529177  /.  u, 10]
N[(-t[[1]] - 99.684205 - 7.465127)*27.2114 ,  
      10]
plot1 = 
      Graphics[{Circle[{0, 0}, P], Circle[{P + Q, 0}, Q], 
            Circle[{P + Q + Q + X, 0}, X], Disk[{0, 0}, 0.08], 
            Disk[{P + Q + Q + X, 0}, 0.08],  
            Line[{{0, 0}, {P + Q + Q + X, 0}}]}]  /.  
        u
Show[
      plot1, {AspectRatio -> Automatic}]