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

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

(* CreatedBy='Mathematica 9.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[      3824,        122]
NotebookOptionsPosition[      2967,         93]
NotebookOutlinePosition[      3663,        116]
CellTagsIndexPosition[      3620,        113]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
(* Fe++ ion Kimball, Ne,Ar,Ni-centers, 6 4s3d-spheres(1) octahedral \
21.12.2011/15.07.2013 *)
Clear[k1,k2,k3,k4,k5,sig1,sig2,sig3,sig4,sig5,c,z,R1,R2,R3,R4,R5,S2,S3,S4];
c = {k1 -> 1.0, k2 -> 0.99, k3 -> 0.97, k4 -> 0.95, k5 -> 0.8,sig1 -> 0.3, \
sig2 -> 0.302, sig3 -> 0.302, sig4 -> 0.302, sig5 -> 0.302}; z=26.;

(* Ne shell *)
T = 2.25*k1/R1^2+9.*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.*z/(R1+R2);
S2 = R2*4^(1/3);

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

(* Fe+++ shell octahedron, full with 6/6 charge each*)
T = T + (6.*9./8)*k4/R4^2 /. c;
Vee = Vee+(90.+15/Sqrt[2])/(S3+R4);
Vne = Vne - 5.6*z/(S3+R4);
func = T + Vee + Vne;

t = FindMinimum[func, {R1,0.051424}, {R2,0.16087}, {R3,0.38299}, \
{R4,0.8377},{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+R3) /. c /. t[[2]],10]
N[(S3+R4) /. c /. t[[2]],10]
\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5828667242007036`*^9, 3.582866793246425*^9}, {
   3.5828671304722176`*^9, 3.582867131174219*^9}, {3.5828672073491526`*^9, 
   3.582867222605979*^9}, 3.582867338623383*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "1265.6477691977166`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"R1", "\[Rule]", "0.05141901843901389`"}], ",", 
     RowBox[{"R2", "\[Rule]", "0.16072269421396745`"}], ",", 
     RowBox[{"R3", "\[Rule]", "0.3794572633858755`"}], ",", 
     RowBox[{"R4", "\[Rule]", "0.8402240494108533`"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.5828668212328744`*^9, 3.582867246926422*^9}],

Cell[BoxData[
 RowBox[{"-", "2944.5642788863856`"}]], "Output",
 CellChangeTimes->{3.5828668212328744`*^9, 3.5828672469420223`*^9}],

Cell[BoxData["413.2687386224652`"], "Output",
 CellChangeTimes->{3.5828668212328744`*^9, 3.5828672469420223`*^9}],

Cell[BoxData["1.9999999985236918`"], "Output",
 CellChangeTimes->{3.5828668212328744`*^9, 3.5828672469420223`*^9}],

Cell[BoxData["0.6345886372562907`"], "Output",
 CellChangeTimes->{3.5828668212328744`*^9, 3.5828672469420223`*^9}],

Cell[BoxData["1.442574908486566`"], "Output",
 CellChangeTimes->{3.5828668212328744`*^9, 3.5828672469420223`*^9}]
}, Open  ]]
},
PrintingStyleEnvironment->"Printout",
WindowSize->{746, 763},
WindowMargins->{{57, Automatic}, {Automatic, 13}},
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) (January 25, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 1315, 39, 592, "Input"],
Cell[1897, 63, 454, 11, 31, "Output"],
Cell[2354, 76, 131, 2, 31, "Output"],
Cell[2488, 80, 113, 1, 31, "Output"],
Cell[2604, 83, 114, 1, 31, "Output"],
Cell[2721, 86, 114, 1, 31, "Output"],
Cell[2838, 89, 113, 1, 31, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)