(* 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[     10952,        252]
NotebookOptionsPosition[      9834,        215]
NotebookOutlinePosition[     10665,        243]
CellTagsIndexPosition[     10622,        240]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
(* [NaCl6/6]infG Kimball, Ne/Ar-centers,  17.2.2012 *)
Clear[k1,k2,k3,k4,k5,sig1,sig2,sig3,sig4,sig5,c,z,R1,R2,R3,R4,R5,S2,S3,S4];
c = {k1 -> 0.99000000, k2 -> 2.65000000, k3 -> 0.99000000, k4 -> 1.22000000, \
k5 -> 3.42000000, 
sig1 -> 0.3000000, sig2 -> 0.3000000, sig3 -> 0.3000000, sig4 -> 0.3000000, \
sig5 -> 0.3000000};
z=11.0000000000; z1=17.0000000000; sq23=1.5874010519682;
ad = Sqrt[3./8.]; M = 1.74756000000000000; (* Madelung constant for rocksalt \
structure *)

(* He and Neon shell of Na+ *)
T = 2.25*k1/R1^2.+9.0*k2/R2^2. /. c;
Vee = 3.0*sig1/R1+12.0*sig2/R2+16.0/(R1+R2)+24.0*ad/(R1+R2) /. c;
Vne = -3.0*z/R1-8.0*z/(R1+R2);
S2 = R2*sq23;

(* He and Neon shell of Cl- *)
T = T + 2.25*k3/R3^2.+9.0*k4/R4^2. /. c;
Vee = Vee + 3.0*sig3/R3+12.0*sig4/R4+16.0/(R3+R4)+24.0*ad/(R3+R4) /. c;
Vne = Vne -3.0*z1/R3-8.0*z1/(R3+R4);
S3 = R4*sq23;

(* Argon shell, full for Cl- *)
T = T + 9.0*k5/R5^2. /. c;
Vee = Vee+12.0*sig5/R5+80.0/(S3+R5)+24.0*ad/(S3+R5) /. c;
Vne = Vne-8.35*z1/(S3+R5);
S4 = R5*sq23;

(* Na+Cl- ionpair; these terms sum up to -1/d *)
d=S2+S4 /. c;
(* Vnn=z*z1/d;
Veeg=(z-1.0)*(z1+1.0)/d;
Vneg=-(z*(z1+1.0)+(z-1.0)*z1)/d;
Vee=Vee+Veeg;
Vne=Vne+Vneg; *)
func = T + Vee + Vne - M/d;

t = FindMinimum[func, {R1,0.136}, {R2,1.131},{R3,0.082228}, {R4,0.33557}, \
{R5,2.2276785},
{Method -> Automatic}, {MaxIterations -> 500}, {WorkingPrecision -> 14}]

N[Vne /. c /. t[[2]],10]
N[Vee /. c /. t[[2]],10]
N[Vnn /. c /. t[[2]],10]
N[-(Vee + Vne - M/d)/T /. c /. t[[2]],10]
N[d*0.529177 /. c /. t[[2]],10]
N[S2*0.529177 /. c /. t[[2]],10]
N[S4*0.529177 /. c /. t[[2]],10]\
\>", "Input",
 CellChangeTimes->{
  3.541223915163168*^9, {3.541224038621784*^9, 3.5412240886666727`*^9}, {
   3.541224392227606*^9, 3.5412245953711624`*^9}, {3.5412249175273285`*^9, 
   3.5412249232525387`*^9}, {3.541225255439522*^9, 3.5412253027700047`*^9}, {
   3.5412253428932753`*^9, 3.5412253501472883`*^9}, {3.54122542501182*^9, 
   3.5412254389582443`*^9}, 3.541225488619931*^9, {3.541225734242363*^9, 
   3.541225763898015*^9}, {3.5412301812838545`*^9, 3.5412302729340153`*^9}, {
   3.5412303362077265`*^9, 3.541230373975393*^9}, {3.5412304249406824`*^9, 
   3.5412304431927147`*^9}, {3.541230484064786*^9, 3.5412304847979875`*^9}, {
   3.5412305371048794`*^9, 3.541230552658107*^9}, {3.5412306241842327`*^9, 
   3.5412306496590767`*^9}, {3.541230699953565*^9, 3.5412307299524183`*^9}, {
   3.541230930974371*^9, 3.5412309378227835`*^9}, {3.5413114875226707`*^9, 
   3.54131154978238*^9}, {3.5413115830640383`*^9, 3.541311584421241*^9}, {
   3.541311679316208*^9, 3.5413117218886824`*^9}, 3.5413118175168505`*^9, 
   3.541311854925716*^9, {3.5413142573436174`*^9, 3.5413143147673182`*^9}, {
   3.5413147386512628`*^9, 3.54131475415769*^9}, {3.5413148884115257`*^9, 
   3.541314894058736*^9}, {3.5413149865824986`*^9, 3.5413150341469817`*^9}, {
   3.541318510208914*^9, 3.5413186224355106`*^9}, {3.541318683166417*^9, 
   3.5413187100452642`*^9}, {3.5413188984467955`*^9, 3.541318899055197*^9}, 
   3.5413189349508595`*^9, {3.5413190146825995`*^9, 3.541319070577498*^9}, {
   3.541319380986843*^9, 3.5413193818760448`*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"FindMinimum", "::", "precw"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"The precision of the argument function \
(\[NoBreak]\\!\\(2.2275`\\/R1\\^2.` - 32.1`\\/R1 + \
23.849999999999998`\\/R2\\^2.` + 3.5999999999999996`\\/R2 - \
57.303061543300934`\\/\\(R1 + R2\\) + 2.2275`\\/R3\\^2.` - 50.1`\\/R3 + \
10.98`\\/R4\\^2.` + 3.5999999999999996`\\/R4 - 105.30306154330094`\\/\\(R3 + \
R4\\) + 30.78`\\/R5\\^2.` + 3.5999999999999996`\\/R5 - 47.253061543300916`\\/\
\\(\\(\\(1.5874010519682`\\\\ R4\\)\\) + R5\\) - \
1.74756`15.954589770191005\\/\\(\\(\\(1.5874010519682`\\\\ R2\\)\\) + \
\\(\\(1.5874010519682`\\\\ R5\\)\\)\\)\\)\[NoBreak]) is less than \
WorkingPrecision (\[NoBreak]\\!\\(14.`\\)\[NoBreak]). \\!\\(\\*ButtonBox[\\\"\
\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \
\\\"FindMinimum::precw\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{
  3.5412249399913673`*^9, 3.541225280415166*^9, 3.541225318994033*^9, 
   3.5412253572765007`*^9, 3.5412254518750668`*^9, 3.541225496170345*^9, 
   3.5412257714952283`*^9, 3.541230279127226*^9, 3.5412303954254303`*^9, {
   3.5412304308842926`*^9, 3.5412304482939234`*^9}, 3.5412304920052*^9, 
   3.5412305657777295`*^9, 3.5412307372220306`*^9, {3.5413115673948107`*^9, 
   3.5413115904584513`*^9}, 3.541311727364292*^9, 3.5413118225712595`*^9, 
   3.541314274737648*^9, 3.5413143196657267`*^9, 3.54131457207417*^9, 
   3.541314764391308*^9, 3.5413149004547467`*^9, {3.5413149957553144`*^9, 
   3.5413150414009943`*^9}, 3.541318631218326*^9, {3.5413186884080267`*^9, 
   3.5413187171588774`*^9}, 3.541318904936407*^9, 3.5413191158331776`*^9, 
   3.5413193876168547`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "572.27295403097806048458839361316175123732`14."}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
     "R1", "\[Rule]", "0.13598731791930248700744855920351950655`14."}], ",", 
     RowBox[{
     "R2", "\[Rule]", "1.13073435077069399046036639270977099125`14."}], ",", 
     RowBox[{
     "R3", "\[Rule]", "0.08222747482613548011195255623362290357`14."}], ",", 
     RowBox[{
     "R4", "\[Rule]", "0.33556645565914517456825198173509848336`14."}], ",", 
     RowBox[{
     "R5", "\[Rule]", "2.22624337734693024232620703957525264659`14."}]}], 
    "}"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.5413115904740515`*^9, 3.5413117273798923`*^9, 3.5413118225556593`*^9, 
   3.541314274737648*^9, 3.5413143196657267`*^9, 3.54131457207417*^9, 
   3.541314764391308*^9, 3.5413149004547467`*^9, {3.5413149957553144`*^9, 
   3.5413150414009943`*^9}, 3.541318631218326*^9, {3.541318688423627*^9, 
   3.5413187171588774`*^9}, 3.5413189049520073`*^9, 3.5413191158487773`*^9, 
   3.5413193876168547`*^9}],

Cell[BoxData[
 RowBox[{"-", "1309.3415894900456`"}]], "Output",
 CellChangeTimes->{
  3.5413115904740515`*^9, 3.5413117273798923`*^9, 3.5413118225556593`*^9, 
   3.541314274737648*^9, 3.5413143196657267`*^9, 3.54131457207417*^9, 
   3.541314764391308*^9, 3.5413149004547467`*^9, {3.5413149957553144`*^9, 
   3.5413150414009943`*^9}, 3.541318631218326*^9, {3.541318688423627*^9, 
   3.5413187171588774`*^9}, 3.5413189049520073`*^9, 3.5413191158487773`*^9, 
   3.5413193876168547`*^9}],

Cell[BoxData["165.12362342221388`"], "Output",
 CellChangeTimes->{
  3.5413115904740515`*^9, 3.5413117273798923`*^9, 3.5413118225556593`*^9, 
   3.541314274737648*^9, 3.5413143196657267`*^9, 3.54131457207417*^9, 
   3.541314764391308*^9, 3.5413149004547467`*^9, {3.5413149957553144`*^9, 
   3.5413150414009943`*^9}, 3.541318631218326*^9, {3.541318688423627*^9, 
   3.5413187171588774`*^9}, 3.5413189049520073`*^9, 3.5413191158487773`*^9, 
   3.5413193876168547`*^9}],

Cell[BoxData["35.09186765776569`"], "Output",
 CellChangeTimes->{
  3.5413115904740515`*^9, 3.5413117273798923`*^9, 3.5413118225556593`*^9, 
   3.541314274737648*^9, 3.5413143196657267`*^9, 3.54131457207417*^9, 
   3.541314764391308*^9, 3.5413149004547467`*^9, {3.5413149957553144`*^9, 
   3.5413150414009943`*^9}, 3.541318631218326*^9, {3.541318688423627*^9, 
   3.5413187171588774`*^9}, 3.5413189049520073`*^9, 3.5413191158487773`*^9, 
   3.5413193876168547`*^9}],

Cell[BoxData["2.000000000080898`"], "Output",
 CellChangeTimes->{
  3.5413115904740515`*^9, 3.5413117273798923`*^9, 3.5413118225556593`*^9, 
   3.541314274737648*^9, 3.5413143196657267`*^9, 3.54131457207417*^9, 
   3.541314764391308*^9, 3.5413149004547467`*^9, {3.5413149957553144`*^9, 
   3.5413150414009943`*^9}, 3.541318631218326*^9, {3.541318688423627*^9, 
   3.5413187171588774`*^9}, 3.5413189049520073`*^9, 3.5413191158487773`*^9, 
   3.541319387632455*^9}],

Cell[BoxData["2.8199154278441894`"], "Output",
 CellChangeTimes->{
  3.5413115904740515`*^9, 3.5413117273798923`*^9, 3.5413118225556593`*^9, 
   3.541314274737648*^9, 3.5413143196657267`*^9, 3.54131457207417*^9, 
   3.541314764391308*^9, 3.5413149004547467`*^9, {3.5413149957553144`*^9, 
   3.5413150414009943`*^9}, 3.541318631218326*^9, {3.541318688423627*^9, 
   3.5413187171588774`*^9}, 3.5413189049520073`*^9, 3.5413191158487773`*^9, 
   3.541319387632455*^9}],

Cell[BoxData["0.9498350894093092`"], "Output",
 CellChangeTimes->{
  3.5413115904740515`*^9, 3.5413117273798923`*^9, 3.5413118225556593`*^9, 
   3.541314274737648*^9, 3.5413143196657267`*^9, 3.54131457207417*^9, 
   3.541314764391308*^9, 3.5413149004547467`*^9, {3.5413149957553144`*^9, 
   3.5413150414009943`*^9}, 3.541318631218326*^9, {3.541318688423627*^9, 
   3.5413187171588774`*^9}, 3.5413189049520073`*^9, 3.5413191158487773`*^9, 
   3.541319387632455*^9}],

Cell[BoxData["1.8700803384348803`"], "Output",
 CellChangeTimes->{
  3.5413115904740515`*^9, 3.5413117273798923`*^9, 3.5413118225556593`*^9, 
   3.541314274737648*^9, 3.5413143196657267`*^9, 3.54131457207417*^9, 
   3.541314764391308*^9, 3.5413149004547467`*^9, {3.5413149957553144`*^9, 
   3.5413150414009943`*^9}, 3.541318631218326*^9, {3.541318688423627*^9, 
   3.5413187171588774`*^9}, 3.5413189049520073`*^9, 3.5413191158487773`*^9, 
   3.541319387632455*^9}]
}, Open  ]]
},
PrintingStyleEnvironment->"Printout",
WindowSize->{1382, 940},
WindowMargins->{{185, Automatic}, {-4, Automatic}},
DockedCells->FEPrivate`FrontEndResource[
 "FEExpressions", "CompatibilityToolbar"],
PrintingCopies->1,
PrintingPageRange->{32000, 32000},
PrintingOptions->{"Magnification"->1.,
"PaperOrientation"->"Portrait",
"PaperSize"->{580.92, 829.0799999999999},
"PrintCellBrackets"->False,
"PrintMultipleHorizontalPages"->False,
"PrintRegistrationMarks"->False,
"PrintingMargins"->{{33.99987438408, 14.125039523099998`}, {56.68752817317, \
56.68752817317}}},
FrontEndVersion->"8.0 for Microsoft Windows (64-bit) (October 6, 2011)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 3131, 70, 762, "Input"],
Cell[3713, 94, 1761, 29, 59, "Message"],
Cell[5477, 125, 1051, 23, 30, "Output"],
Cell[6531, 150, 483, 8, 30, "Output"],
Cell[7017, 160, 466, 7, 30, "Output"],
Cell[7486, 169, 465, 7, 30, "Output"],
Cell[7954, 178, 463, 7, 30, "Output"],
Cell[8420, 187, 464, 7, 30, "Output"],
Cell[8887, 196, 464, 7, 30, "Output"],
Cell[9354, 205, 464, 7, 30, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)