(* P atom nach Kimball, Aussen-sp-Kugeln 18.12.2011/10.04.2013 *)
Clear[k1,k2,k3,sig1,sig2,sig3,c,z,R1,R2,R3,a];
c = {k1 -> 0.9905, k2 -> 0.9907, k3->0.99108, sig1 -> 0.3, sig2 -> 0.32, sig3 -> 0.32};

T = 2.25*k1/R1^2+9.0*k2/R2^2+5.625*k3/R3^2 /. c;
a = (R1+R2)/3+Sqrt[(R2+R3)^2-8*(R1+R2)^2/9];
Vee = 3.0*sig1/R1+4.0*3.0*sig2/R2+3.0*sig3/R3+16/(R1+R2)+24/((R1+R2)*Sqrt[8/3])+
    10/a+30/(R2+R3)+10/(R1+R2+a)+9*Sqrt[3/8]/a /. c;
Vne=-3.0*z/R1-8*z/(R1+R2)-5*z/a;

func = T + Vee + Vne;

t = Table[FindMinimum[func, {R1,0.04}, {R2,0.25}, {R3,0.6}, {MaxIterations -> 500, Method -> "Newton"}], {z,15,26}]

Table[N[Vne /. c /. t[[{z-14},2]],10], {z,15,26}]
Table[N[Vee /. c /. t[[{z-14},2]],10], {z,15,26}]
Table[N[-(Vee+Vne)/T /. c /. t[[{z-14},2]],10], {z,15,26}]
Table[N[a /. c /. t[[{z-14},2]],10], {z,15,26}]