Clear[d12, d13, d14, d23, d24, 
        d34];
c2 = {k1 ->  1.02480344,  k2 ->  1.232443643, 
          s1 ->  0.27877144, 
          s2 ->  0.28887416};
T  =  k2*2.25/R^2  + k1*2.25/P^2;
d12  = P + Q + R - x;
d13  =  x;
d14  =  P + R - x;
d23  =  P + R + Q;
d24  =  Q;
d34  = P + R;

Vee  =  3*s2/R + 3*s1/P + 4/d34;
Vnn  =  3/d12;
Vne  =  -3*(3 - (d13/P)^2)/P - 
          6/d14 - (3 - (d24/R)^2)/R - 
          2/d23;

func  =  
      T  +  Vee  +  Vnn  +  Vne  /.  c2
   
N[t  =  
        FindMinimum[
          func,  {R, {1.8, 2.2}}, {P, {0.5, 0.6}}, {Q, {0.4, 
              0.5}}, {x, {0.01, 0.02}}, {MaxIterations  ->  100}], 
      20]  
u  =  t[[2]];
  
d12   /. u
  
d13  /.  u
d14  /.  u
d23  /.  u
d24  /.  u
d34  /. u
N[(Vnn + Vee + Vne)/T  /.  u  /.  c2,  
      10]
  
N[-6*x/P^3 - 3/d12^2 + 6/d14^2  /.  u ,  
      10] 
  
N[2*Q/R^3 + 2/d23^2 - 3/d12^2  /.  u,  10]