(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 8.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[      4763,        163]
NotebookOptionsPosition[      3753,        130]
NotebookOutlinePosition[      4449,        153]
CellTagsIndexPosition[      4406,        150]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
(* Cs atom Kimball, Ne,Ar,Ni,Kr,Pd-centers, 4 5sp-spheres tetrahedral \
28.12.2011 
   vne contribution of each shell roughly parametrized to fda *)
Clear[k1,k2,k3,k4,k5,k6,k7,k8sig1,sig2,sig3,sig4,sig5,sig6,sig7,c,z,R1,R2,R3,\
R4,R5,R6,R7,R8S2,S3,S4,S5,S6,S7];
c = {k1 -> 1.0, k2 -> 1.0, k3 -> 1.0, k4 -> 1.0, k5 -> 1.0, k6 -> 1.0, k7 -> \
1.0, k8 -> 1.0, sig1 -> 0.3, sig2 -> 0.3, sig3 -> 0.3, sig4 -> 0.3, sig5 -> \
0.3, sig6 -> 0.3, sig7 -> 0.3}; z=55.;

(* He+Ne shell *)
T = (2.*9./8.)*k1/R1^2+(8.*9./8.)*k2/R2^2 /. c;
ad = Sqrt[3./8.]; 
Vee=3.0*sig1/R1+12.*sig2/R2+16/(R1+R2)+24*ad/(R1+R2) /. c;
Vne=-3.0*z/R1-8.0*z/(R1+R2);
S2 = R2*4^(1/3);

(* Ar shell *)
T = T + (8.*9./8.)*k3/R3^2 /. c;
Vee = Vee+12.*sig3/R3+80./(S2+R3)+24.*ad/(S2+R3) /. c;
Vne = Vne-8.35*z/(S2+R3);
S3 = R3*4^(1/3);

(* Ni shell trigbipyr, d10 with 5*2 charges; *)
T = T + (10.*9./8.)*k4/R4^2 /. c;
Vee = Vee+5.*3.*sig4/R4+(180.+25.898766)/(S3+R4) /. c; (* \
2+24/Sqrt[2]+12/Sqrt[3] *)
Vne = Vne - 10.5*z/(S3+R4);
S4 = R4*5^(1/3);

(* Kr shell *)
T = T + (8.*9./8.)*k5/R5^2 /. c;
Vee = Vee+12.*sig5/R5+224./(S4+R5)+24.*ad/(S4+R5) /. c;
Vne = Vne-9.7*z/(S4+R5);
S5 = R5*4^(1/3);

(* Pd shell trigbipyr, d10 with 5*2 charges *)
T = T + (10.*9./8)*k6/R6^2 /. c;
Vee = Vee+5.*3.*sig6/R6+(360.+25.898766)/(S5+R6) /. c;
Vne = Vne - 10.95*z/(S5+R6);
S6 = R6*5^(1/3);

(* Xe shell *)
T = T + (8.*9./8.)*k7/R7^2 /. c;
Vee = Vee+12.*sig7/R7+(368+24.*ad)/(S6+R7) /. c;
Vne = Vne-8.3*z/(S6+R7);
S7 = R7*4^(1/3);

(* Cs outer sphere *)
T = T + 1.125*k8/R8^2 /. c;
Vee = Vee+54/(S7+R8);
Vne = Vne - z/(S7+R8);


func = T + Vee + Vne;

t = FindMinimum[func, {R1,0.0234209}, {R2,0.0660675}, {R3,0.1135267}, \
{R4,0.1653891},{R5,0.2159044},{R6,0.3677186},{R7,0.8571924},{R8,4.0272848},{\
Method -> Automatic}, {MaxIterations -> 500}]

N[Vne /. c /. t[[2]],10]
N[Vee /. c /. t[[2]],10]
N[-(Vee+Vne)/T /. c /. t[[2]],10]
N[(S2) /. c /. t[[2]],10]
N[(S3) /. c /. t[[2]],10]
N[(S4) /. c /. t[[2]],10]
N[(S5) /. c /. t[[2]],10]
N[(S6) /. c /. t[[2]],10]
N[S7 /. c /. t[[2]],10]
\
\>", "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "7561.888135605378`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"R1", "\[Rule]", "0.023420888686931596`"}], ",", 
     RowBox[{"R2", "\[Rule]", "0.0660674573993974`"}], ",", 
     RowBox[{"R3", "\[Rule]", "0.1135267355417589`"}], ",", 
     RowBox[{"R4", "\[Rule]", "0.1653891429512644`"}], ",", 
     RowBox[{"R5", "\[Rule]", "0.21590443220613187`"}], ",", 
     RowBox[{"R6", "\[Rule]", "0.36771858876367874`"}], ",", 
     RowBox[{"R7", "\[Rule]", "0.8571924003590957`"}], ",", 
     RowBox[{"R8", "\[Rule]", "4.027284799996753`"}]}], "}"}]}], 
  "}"}]], "Output"],

Cell[BoxData[
 RowBox[{"-", "17970.471636108236`"}]], "Output"],

Cell[BoxData["2846.6953548597903`"], "Output"],

Cell[BoxData["1.9999999986725943`"], "Output"],

Cell[BoxData["0.10487555137666761`"], "Output"],

Cell[BoxData["0.18021245942550362`"], "Output"],

Cell[BoxData["0.28281145628813587`"], "Output"],

Cell[BoxData["0.34272692280861045`"], "Output"],

Cell[BoxData["0.6287899419317905`"], "Output"],

Cell[BoxData["1.3607081180691742`"], "Output"]
}, Open  ]]
},
PrintingStyleEnvironment->"Printout",
WindowSize->{874, 854},
WindowMargins->{{2, Automatic}, {Automatic, 10}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
PrintingOptions->{"PrintCellBrackets"->False,
"PrintMultipleHorizontalPages"->False,
"PrintRegistrationMarks"->False,
"PrintingMargins"->{{34, 14.125}, {56.6875, 56.6875}}},
PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}},
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (November 20, 2012)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 2061, 70, 1068, "Input"],
Cell[2643, 94, 632, 14, 52, "Output"],
Cell[3278, 110, 63, 1, 31, "Output"],
Cell[3344, 113, 46, 0, 31, "Output"],
Cell[3393, 115, 46, 0, 31, "Output"],
Cell[3442, 117, 47, 0, 31, "Output"],
Cell[3492, 119, 47, 0, 31, "Output"],
Cell[3542, 121, 47, 0, 31, "Output"],
Cell[3592, 123, 47, 0, 31, "Output"],
Cell[3642, 125, 46, 0, 31, "Output"],
Cell[3691, 127, 46, 0, 31, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)