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 - 12*s1*Q^2/P^3;
Vnn = Z*Z/(2*R);
func = T + Vne + Vee + Vnn;
  
c = {k1  ->  0.9852,  k2 ->  1.0,  s1 -> 0.32, 
          s2 ->  0.39,  
          Z  ->  3.0};
func  =  func  /.  c
N[t = FindMinimum[func,  {P, 5.04764}, {Q, 0.56255}, {R, 2.5262}], 
      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, Axes -> True, 
        GridLines -> Automatic, 
        PlotRange -> {{-6, 6}, {-6, 6}}, 
        Frame -> True}]