(* 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[ 17367, 512] NotebookOptionsPosition[ 16756, 489] NotebookOutlinePosition[ 17184, 506] CellTagsIndexPosition[ 17141, 503] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{"Li", "-", RowBox[{"Atom", " ", "and", " ", "homologues", " ", "1", "s", RowBox[{"(", "2", ")"}], "2", "s", " ", "concentric", " ", "spherical"}]}], ",", RowBox[{ RowBox[{ "with", " ", "exchange", " ", "correction", " ", "page", " ", "40"}], ";", " ", RowBox[{"4.12", ".2011"}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Clear", "[", RowBox[{"z", ",", "sig"}], "]"}], ";", " ", RowBox[{"sig", "=", "0.29"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"T", "=", RowBox[{ RowBox[{"9", "/", RowBox[{"(", RowBox[{"4", "*", RowBox[{"P", "^", "2"}]}], ")"}]}], "+", RowBox[{"9", "*", RowBox[{"4", "/", RowBox[{"(", RowBox[{"8", "*", RowBox[{"Q", "^", "2"}]}], ")"}]}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vne", "=", RowBox[{ RowBox[{ RowBox[{"-", "3"}], "*", RowBox[{"z", "/", "P"}]}], "-", RowBox[{"3", "*", RowBox[{"z", "/", RowBox[{"(", RowBox[{"2", "*", "Q"}], ")"}]}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vee", "=", RowBox[{ RowBox[{ RowBox[{"+", RowBox[{"(", RowBox[{"2", "/", "Q"}], ")"}]}], "*", RowBox[{"(", RowBox[{ RowBox[{"3", "/", "2"}], "-", RowBox[{ RowBox[{"(", RowBox[{"3", "/", "10"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "/", "Q"}], ")"}], "^", "2"}]}]}], ")"}]}], "-", RowBox[{"3", "*", "sig", "*", RowBox[{ RowBox[{"P", "^", "2"}], "/", RowBox[{"Q", "^", "3"}]}]}], " ", "+", RowBox[{"3", "*", RowBox[{"sig", "/", "P"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"func", "=", RowBox[{"T", "+", "Vne", "+", "Vee"}]}], ";", " ", RowBox[{"Vir", " ", "=", " ", RowBox[{ RowBox[{"(", RowBox[{"Vne", "+", "Vee"}], ")"}], "/", "T"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]"}], RowBox[{"(*", " ", RowBox[{"Etot", ",", " ", "P", ",", " ", "Q"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"t", " ", "=", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{"FindMinimum", "[", RowBox[{"func", ",", RowBox[{"{", RowBox[{"P", ",", "0.6"}], "}"}], ",", RowBox[{"{", RowBox[{"Q", ",", "4"}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"z", ",", "3", ",", "18"}], "}"}]}], "]"}]}], "\[IndentingNewLine]", "\[IndentingNewLine]"}], RowBox[{"(*", " ", RowBox[{"Virial", " ", "ratio"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"N", "[", RowBox[{ RowBox[{"Vir", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "2"}], " ", "}"}], ",", "2"}], "]"}], "]"}]}], ",", "10"}], "]"}], ",", RowBox[{"{", RowBox[{"z", " ", ",", "3", ",", "18"}], "}"}]}], "]"}], "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"N", "[", RowBox[{ RowBox[{"Vne", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "2"}], " ", "}"}], ",", "2"}], "]"}], "]"}]}], ",", "10"}], "]"}], ",", RowBox[{"{", RowBox[{"z", " ", ",", "3", ",", "18"}], "}"}]}], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"N", "[", RowBox[{ RowBox[{"Vee", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "2"}], " ", "}"}], ",", "2"}], "]"}], "]"}]}], ",", "10"}], "]"}], ",", RowBox[{"{", RowBox[{"z", " ", ",", "3", ",", "18"}], "}"}]}], "]"}], "\[IndentingNewLine]", "\[IndentingNewLine]"}], RowBox[{"(*", " ", RowBox[{"Etot", " ", "in", " ", "eV"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"N", "[", RowBox[{ RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "2"}], "}"}], ",", "1"}], "]"}], "]"}], "*", "27.21138"}], ",", "10"}], "]"}], ",", RowBox[{"{", RowBox[{"z", ",", "3", ",", "18"}], "}"}]}], "]"}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"PA", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "2"}], "}"}], ",", "2", ",", "1"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"z", ",", "3", ",", "18"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"QA", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "2"}], "}"}], ",", "2", ",", "2"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"z", ",", "3", ",", "18"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Do", "[", RowBox[{ RowBox[{"Show", "[", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"P", "/.", " ", RowBox[{"PA", "[", RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"Q", " ", "/.", " ", RowBox[{"QA", "[", RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", RowBox[{"Axes", " ", "->", " ", "True"}], ",", RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "6"}], ",", "6"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "6"}], ",", "6"}], "}"}]}], "}"}]}], ",", " ", RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", "1", ",", "16"}], "}"}]}], "]"}], "\[IndentingNewLine]"}]}]}]], "Input"], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lstol"}], "MessageName"], RowBox[{ ":", " "}], "\<\"The line search decreased the step size to within the \ tolerance specified by AccuracyGoal and PrecisionGoal but was unable to find \ a sufficient decrease in the function. You may need more than \ \[NoBreak]\\!\\(MachinePrecision\\)\[NoBreak] digits of working precision to \ meet these tolerances. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lstol\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.541087277536845*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lstol"}], "MessageName"], RowBox[{ ":", " "}], "\<\"The line search decreased the step size to within the \ tolerance specified by AccuracyGoal and PrecisionGoal but was unable to find \ a sufficient decrease in the function. You may need more than \ \[NoBreak]\\!\\(MachinePrecision\\)\[NoBreak] digits of working precision to \ meet these tolerances. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lstol\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.541087277536845*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "7.471268110264797`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.5536756460973181`"}], ",", RowBox[{"Q", "\[Rule]", "5.845825636762729`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "14.273383530314938`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.4045977073635724`"}], ",", RowBox[{"Q", "\[Rule]", "2.9175196667627397`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "23.328486304007562`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.31876153657133016`"}], ",", RowBox[{"Q", "\[Rule]", "1.9489063499094903`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "34.63525431988847`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.26295784257489596`"}], ",", RowBox[{"Q", "\[Rule]", "1.465316072098628`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "48.193014651325264`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.22377272454754635`"}], ",", RowBox[{"Q", "\[Rule]", "1.174940305409747`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "64.00140705094233`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.19474615550003446`"}], ",", RowBox[{"Q", "\[Rule]", "0.9810574039382058`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "82.06022473952459`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.17238208425255291`"}], ",", RowBox[{"Q", "\[Rule]", "0.8423260944198183`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "102.36934149755291`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.15462339002076406`"}], ",", RowBox[{"Q", "\[Rule]", "0.73809594754815`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "124.92867623286747`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.140180675200997`"}], ",", RowBox[{"Q", "\[Rule]", "0.6568946228012585`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "149.73817459025213`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.12820470917136068`"}], ",", RowBox[{"Q", "\[Rule]", "0.5918350303370835`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "176.7977988241952`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.11811339197844363`"}], ",", RowBox[{"Q", "\[Rule]", "0.5385307814601937`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "206.10752192957304`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.10949438873548693`"}], ",", RowBox[{"Q", "\[Rule]", "0.49405468548800996`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "237.6673240876434`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.10204745398342989`"}], ",", RowBox[{"Q", "\[Rule]", "0.45637805914196744`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "271.4771904313095`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.09554877121594765`"}], ",", RowBox[{"Q", "\[Rule]", "0.4240502515136475`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "307.5371095932887`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.0898280949349076`"}], ",", RowBox[{"Q", "\[Rule]", "0.3960062729398866`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "345.84707273564646`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.08475362327918382`"}], ",", RowBox[{"Q", "\[Rule]", "0.37144657156327526`"}]}], "}"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.5410872775524454`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"-", "1.9999999886482893`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999988427315`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999999770892`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "2.000000000177633`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "2.0000000008552425`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "2.0000000012599557`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "2.000000003456254`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999974676255`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999997676778`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "2.000000000293964`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "2.0000000008338223`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "2.000000000579578`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999999605091`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.999999986079369`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999970867284`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "2.000000002763393`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.5410872775680456`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"-", "17.02478434372861`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "31.715631860413094`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "50.90543254069606`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "74.59405669255355`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "102.78184170257293`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "135.4690495055627`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "172.65585292607298`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "214.342366238848`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "260.5286637383898`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "311.21479847610857`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "366.4008070723234`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "426.0867163734708`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "490.2725459837656`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "558.9583143224937`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "632.1440225875147`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "709.8296872459407`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.5410872775680456`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", "2.0822480383873425`", "}"}], ",", RowBox[{"{", "3.1688647832650805`", "}"}], ",", RowBox[{"{", "4.248459932146458`", "}"}], ",", RowBox[{"{", "5.323548058928986`", "}"}], ",", RowBox[{"{", "6.395812441139123`", "}"}], ",", RowBox[{"{", "7.4662354843169485`", "}"}], ",", RowBox[{"{", "8.535403730644765`", "}"}], ",", RowBox[{"{", "9.603682984504697`", "}"}], ",", RowBox[{"{", "10.671311243631223`", "}"}], ",", RowBox[{"{", "11.738449339621937`", "}"}], ",", RowBox[{"{", "12.805209571351014`", "}"}], ",", RowBox[{"{", "13.871672633780172`", "}"}], ",", RowBox[{"{", "14.937897799093115`", "}"}], ",", RowBox[{"{", "16.003929680740946`", "}"}], ",", RowBox[{"{", "17.069802504998094`", "}"}], ",", RowBox[{"{", "18.135542730359262`", "}"}]}], "}"}]], "Output", CellChangeTimes->{3.5410872775680456`*^9}] }, Open ]] }, WindowSize->{690, 815}, WindowMargins->{{341, Automatic}, {Automatic, 26}}, DockedCells->FEPrivate`FrontEndResource[ "FEExpressions", "CompatibilityToolbar"], 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, 7189, 219, 532, "Input"], Cell[7771, 243, 681, 13, 73, "Message"], Cell[8455, 258, 681, 13, 73, "Message"], Cell[9139, 273, 4189, 118, 164, "Output"], Cell[13331, 393, 1243, 35, 50, "Output"], Cell[14577, 430, 1236, 35, 69, "Output"], Cell[15816, 467, 924, 19, 69, "Output"] }, Open ]] } ] *) (* End of internal cache information *)