(* K atom nach Kimball, Aussen-sp-Kugeln 19.12.2011 *)
Clear[k1,k2,k3,k4,sig1,sig2,sig3,c,z,R1,R2,R3,R4,a,b,k];
k = {k1 -> 1.0, k2 -> 1.0, k3->1.0, k4->1.0, sig1 -> 0.4, sig2 -> 0.4, sig3 -> 0.4};
z=21.0;
T = 2.25*k1/R1^2+9.0*k2/R2^2+9.0*k3/R3^2+1.125*k4/R4^2 /. k;
a = (R1+R2)/3+Sqrt[(R2+R3)^2-8*(R1+R2)^2/9];
b = a/3+Sqrt[(R3+R4)^2-8*a*a/9];
c = b+Sqrt[((R1+R2)/3)^2+8*(R1+R2)^2/9];
Vee = 3.0*sig1/R1+4.0*3.0*sig2/R2+4.0*3.0*sig3/R3+16/(R1+R2)+24/((R1+R2)*Sqrt[8/3])+
      16/a+48/(R2+R3)+16/(R1+R2+a)+24*Sqrt[3/8]/a+2/b+6/(R3+R4)+2/(a+b)+2/(b-R1-R2)+6/c /. k;
Vne = -3.0*z/R1-8*z/(R1+R2)-8*z/a-z/b;
Vir = (Vee+Vne)/T;
func = T+Vee+Vne;

t = FindMinimum[func, {R1,{0.064,0.066}}, {R2,{0.218,0.221}}, {R3,{0.54,0.55}}, {R4,{0.96,1.1}}, {Method -> Automatic},{MaxIterations -> 500}]


N[-Vir /. t[[2]],10]
Show[Graphics[{Circle[{0,0},R1],
      Circle[{0,0},R2],Circle[{0,0},R3],Circle[{0,R1+R2},R2],
        Circle[{-0.866*(R1+R2),-(R1+R2)/2},R2],Circle[{0.866*(R1+R2),-(R1+R2)/2},R2],
        Circle[{-0.866*(R2+R3),(R2+R3)/2},R3],Circle[{0.866*(R2+R3),(R2+R3)/2},R3],
        Circle[{0,-(R2+R3)},R3],Circle[{0,0},R3],
        Circle[{-0.866*(R3+R4),-(R3+R4)/2},R4]}] /. t[[2]],{AspectRatio -> Automatic,
    Axes -> True,GridLines -> Automatic, PlotRange -> {{-2.5,1.5},{-2.0,1.0}}, Frame -> True}]