T  =  
        2.25/R^2;
Vee  =  1.2/R;
d12  =  
        Sqrt[(x1 - x2)^2 + (y1 - y2)^2 + (z1 - 
                  z2)^2];
Vnn  =  1/d12;
r1sq  =  x1^2 + y1^2 + z1^2;
r2sq  =  
        x2^2 + y2^2 + z2^2;

Vne  =  -(6  -  (r1sq + r2sq)/R^2)/
          R;

func  =  T  +  Vee  +  Vnn  +  Vne;
   
N[t  =  
        FindMinimum[
          func,  {R, {1.0, 1.5}}, 
     {x1, {0.4, 
              2}}, {y1, {0.1, 2}}, {z1, {0.2, 
              2}}, {x2, {-0.3, -2}}, {y2, {0.8, 
              2}}, {z2, {-0.3, -2}}, {MaxIterations  ->  
              100}], 20]  
u  =  t[[2]];
  
Sqrt[ x1^2 + y1^2 + z1^2]  /.  u
Sqrt[x2^2 + y2^2 + z2^2]  /.  u
  
(x1 + x2)/2  /.  u
(y1 + y2)/2 /.  u
(z1 + z2)/
        2 /. u