(* 
      generalized H2 with free start of nuclei *) 
 (* n   nuclear charges 0.3  ..  2.2,  
      0.1 *) 
T  =  
        2.25/R^2;
Vee  =  1.2/R;
d12  =  
        Sqrt[(x1 - x2)^2 + (y1 - y2)^2 + (z1 - 
                  z2)^2];
Vnn  =  n*n/d12;
r1sq  =  x1^2 + y1^2 + z1^2;
r2sq  =  
        x2^2 + y2^2 + z2^2;

Vne  =  -(3*(n + n)  -  (n*r1sq + n*r2sq)/R^2)/
          R;

func  =  T  +  Vee  +  Vnn  +  Vne;
   
t  =  
        Table[FindMinimum[
              func,  {R, {0.5, 
                  1.5}}, 
     {x1, {0.4, 
                  2}}, {y1, {0.1, 2}}, {z1, {0.2, 
                  2}}, {x2, {-0.3, -2}}, {y2, {0.3, 
                  2}}, {z2, {-0.3, -2}}, {MaxIterations  ->  
                  100}] + n^2, {n, 0.3, 2.2, 
            0.1}] ; 
 u  =  
        t[[{1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 
              19, 20}, 1]];

v  =  
        Table[n, {n, 0.3, 2.2, 0.1}];

ListPlot[{{v[[1]], 
          u[[1]]}, {v[[2]], u[[2]]}, {v[[3]], 
          u[[3]]}, {v[[4]], u[[4]]}, {v[[5]], 
          u[[5]]}, {v[[6]], u[[6]]}, {v[[7]], 
          u[[7]]}, {v[[8]], u[[8]]}, {v[[9]], 
          u[[9]]}, {v[[10]], u[[10]]}, {v[[11]], 
          u[[11]]}, {v[[12]], u[[12]]}, {v[[13]], 
          u[[13]]}, {v[[14]], u[[14]]}, {v[[15]], 
          u[[15]]}, {v[[16]], u[[16]]}, {v[[17]], 
          u[[17]]}, {v[[18]], u[[18]]}, {v[[19]], 
          u[[19]]}, {v[[20]], u[[20]]}}]