(* Ne like ions F- to Ar(+8) nach Kimball 30.04.2011 *)

Clear[k0,k2,sig0,sig2,c]; n = 1.;
 c = {k0 -> 0.989, k2 -> 1.1,
     sig0 -> 0.3, sig2 -> 0.3};

T = 2.25*n*k0/P^2+4.5*(n+1.0)*k2/R^2;

Vee=3.0*n*sig0/P+(n+1.0)*6.0*sig2/R+16/(P+R)+24/Sqrt[8/3]/(P+R);

Vne=-3.0*n*z/P-8*z/(P+R);
      
func = T + Vee + Vne /. c;

t=Table[FindMinimum[func, {P, 0.09},{R, 0.5}], {z,9,18}]
                        

u = t[[{1,2,3,4,5,6,7,8,9,10},2]];
v = u[[{1,2,3,4,5,6,7,8,9,10},1]];
w = u[[{1,2,3,4,5,6,7,8,9,10},2]];
x = Table[N[P+R /. v[[i]] /. w[[i]]],{i,1,10}]
(* Kriterium für Overlap der Aussenwolken: muss <= 0 sein *)
y = Table[N[(Sqrt[3/2]-1)*R-P /. v[[i]] /. w[[i]]],{i,1,10}]