(* 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[ 18599, 544] NotebookOptionsPosition[ 18029, 522] NotebookOutlinePosition[ 18457, 539] CellTagsIndexPosition[ 18414, 536] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"F", "-", RowBox[{"Atom", " ", "and", " ", "homologues", " ", "1", "s", RowBox[{"(", "2", ")"}], "2", "s", RowBox[{"(", "2", ")"}], " ", "2", "p", RowBox[{"(", "5", ")"}]}]}], ",", " ", RowBox[{"all", " ", "p", " ", "equal", " ", "size", " ", RowBox[{"(", "?", ")"}]}], ",", RowBox[{ RowBox[{ "with", " ", "exchange", " ", "correction", " ", "page", " ", "40"}], ";", " ", RowBox[{"8.12", ".2011"}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Clear", "[", RowBox[{"z", ",", "sig"}], "]"}], ";", " ", RowBox[{"sig", "=", "0.24"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"T", "=", RowBox[{ RowBox[{"9", "/", RowBox[{"(", RowBox[{"4", "*", RowBox[{"P", "^", "2"}]}], ")"}]}], "+", RowBox[{"9", "*", RowBox[{"4", "/", RowBox[{"(", RowBox[{"4", "*", RowBox[{"Q", "^", "2"}]}], ")"}]}]}], "+", RowBox[{"10", "/", RowBox[{"R", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vne", "=", RowBox[{ RowBox[{ RowBox[{"-", "3"}], "*", RowBox[{"z", "/", "P"}]}], "-", RowBox[{"3", "*", RowBox[{"z", "/", "Q"}]}], "-", RowBox[{"5", "*", RowBox[{"z", "/", "R"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vee", "=", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"2", "/", "Q"}], ")"}], "*", RowBox[{"(", RowBox[{"3", "-", RowBox[{ RowBox[{"(", RowBox[{"6", "/", "10"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "/", "Q"}], ")"}], "^", "2"}]}]}], ")"}]}], "-", RowBox[{"6", "*", "sig", "*", RowBox[{ RowBox[{"P", "^", "2"}], "/", RowBox[{"Q", "^", "3"}]}]}], " ", "+", " ", RowBox[{"3", "*", RowBox[{"sig", "/", "P"}]}], "+", RowBox[{"3", "*", RowBox[{"sig", "/", "Q"}]}], "+", RowBox[{"10", "*", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{ RowBox[{"(", RowBox[{"3", "/", "20"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "/", "R"}], ")"}], "^", "2"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "32"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "/", "R"}], ")"}], "^", "3"}]}]}], ")"}], "/", "R"}]}], "+", RowBox[{"3", "*", RowBox[{"sig", "/", "R"}]}], "+", RowBox[{"10", "*", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{ RowBox[{"(", RowBox[{"3", "/", "20"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"Q", "/", "R"}], ")"}], "^", "2"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "32"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"Q", "/", "R"}], ")"}], "^", "3"}]}]}], ")"}], "/", "R"}]}], "+", RowBox[{"4.974012", "/", "R"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"func", "=", RowBox[{"T", "+", "Vne", "+", "Vee"}]}], ";", " ", RowBox[{"Vir", " ", "=", " ", RowBox[{ RowBox[{"(", RowBox[{"Vne", "+", "Vee"}], ")"}], "/", "T"}]}], ";"}], "\[IndentingNewLine]"}], RowBox[{"(*", " ", RowBox[{"Etot", ",", " ", "P", ",", " ", "Q"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"t", " ", "=", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{"FindMinimum", "[", RowBox[{"func", ",", RowBox[{"{", RowBox[{"P", ",", "0.23"}], "}"}], ",", RowBox[{"{", RowBox[{"Q", ",", "1.4"}], "}"}], ",", RowBox[{"{", RowBox[{"R", ",", "1.4"}], "}"}], ",", RowBox[{"{", RowBox[{"MaxIterations", "\[Rule]", "800"}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"z", ",", "9", ",", "18"}], "}"}]}], "]"}]}], "\[IndentingNewLine]"}], RowBox[{"(*", " ", RowBox[{"Virial", " ", "ratio"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"N", "[", RowBox[{ RowBox[{"Vir", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "8"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", "10"}], "]"}], ",", RowBox[{"{", RowBox[{"z", " ", ",", "9", ",", "18"}], "}"}]}], "]"}], "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"N", "[", RowBox[{ RowBox[{"Vne", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "8"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", "10"}], "]"}], ",", RowBox[{"{", RowBox[{"z", " ", ",", "9", ",", "18"}], "}"}]}], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"N", "[", RowBox[{ RowBox[{"Vee", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "8"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", "10"}], "]"}], ",", RowBox[{"{", RowBox[{"z", " ", ",", "9", ",", "18"}], "}"}]}], "]"}], "\[IndentingNewLine]"}], RowBox[{"(*", " ", RowBox[{"Etot", " ", "in", " ", "eV"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"N", "[", RowBox[{ RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "8"}], "}"}], ",", "1"}], "]"}], "]"}], "*", "27.21138"}], ",", "10"}], "]"}], ",", RowBox[{"{", RowBox[{"z", ",", "9", ",", "18"}], "}"}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"PA", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "8"}], "}"}], ",", "2", ",", "1"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"z", ",", "9", ",", "18"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"QA", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "8"}], "}"}], ",", "2", ",", "2"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"z", ",", "9", ",", "18"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"RA", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "8"}], "}"}], ",", "2", ",", "3"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"z", ",", "9", ",", "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[{"Dashing", "[", RowBox[{"{", RowBox[{"0.01", ",", "0.01"}], "}"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{ RowBox[{"-", "R"}], "/.", " ", RowBox[{"RA", "[", RowBox[{"[", "i", "]"}], "]"}]}]}], "}"}], ",", RowBox[{"R", " ", "/.", " ", RowBox[{"RA", "[", RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "R"}], "/.", " ", RowBox[{"RA", "[", RowBox[{"[", "i", "]"}], "]"}]}], ",", "0"}], "}"}], ",", RowBox[{"R", " ", "/.", " ", RowBox[{"RA", "[", RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"R", "/.", " ", RowBox[{"RA", "[", RowBox[{"[", "i", "]"}], "]"}]}], ",", "0"}], "}"}], ",", RowBox[{"R", " ", "/.", " ", RowBox[{"RA", "[", RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{"R", "/.", " ", RowBox[{"RA", "[", RowBox[{"[", "i", "]"}], "]"}]}]}], "}"}], ",", RowBox[{"R", " ", "/.", " ", RowBox[{"RA", "[", RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"R", " ", "/.", " ", RowBox[{"RA", "[", RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}]}], "}"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", RowBox[{"Axes", " ", "->", " ", "True"}], ",", RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "2"}], ",", "2"}], "}"}]}], "}"}]}], ",", " ", RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", "1", ",", "10"}], "}"}]}], "]"}], "\[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.5410877543984833`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "99.118413591089`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.17151811516117538`"}], ",", RowBox[{"Q", "\[Rule]", "1.7802313774807526`"}], ",", RowBox[{"R", "\[Rule]", "0.7107404555821463`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "127.20192653643646`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.15396951052221758`"}], ",", RowBox[{"Q", "\[Rule]", "1.3611975054186312`"}], ",", RowBox[{"R", "\[Rule]", "0.6037439931858775`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "159.14654426899915`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.13968594764946102`"}], ",", RowBox[{"Q", "\[Rule]", "0.8707916972219262`"}], ",", RowBox[{"R", "\[Rule]", "0.5593241365504921`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "194.884959735859`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.12781953650869107`"}], ",", RowBox[{"Q", "\[Rule]", "0.7535184461457581`"}], ",", RowBox[{"R", "\[Rule]", "0.4919673296981902`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "234.37651027426395`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.11780849308316321`"}], ",", RowBox[{"Q", "\[Rule]", "0.6670897079094822`"}], ",", RowBox[{"R", "\[Rule]", "0.4383819123368835`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "277.6203315993553`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.10924960630320682`"}], ",", RowBox[{"Q", "\[Rule]", "0.5991994427155761`"}], ",", RowBox[{"R", "\[Rule]", "0.3951616872247946`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "324.61594455821455`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.10184849895806411`"}], ",", RowBox[{"Q", "\[Rule]", "0.5441381818608664`"}], ",", RowBox[{"R", "\[Rule]", "0.3596457726765442`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "375.36301460264855`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.09538529534169142`"}], ",", RowBox[{"Q", "\[Rule]", "0.4984799920331864`"}], ",", RowBox[{"R", "\[Rule]", "0.3299671152448272`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "429.86129067160755`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.08969246447886695`"}], ",", RowBox[{"Q", "\[Rule]", "0.4599638627322984`"}], ",", RowBox[{"R", "\[Rule]", "0.3048049535184025`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "488.11057786826547`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.08464011821830424`"}], ",", RowBox[{"Q", "\[Rule]", "0.4270159570053794`"}], ",", RowBox[{"R", "\[Rule]", "0.2832050084743732`"}]}], "}"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.541087754414083*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"-", "1.9999999995279782`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "2.0000000012833032`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999987714148`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999997376123`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "2.0000000000735496`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999990948696`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "2.0000000011994308`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999996916868`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "2.0000000006419065`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999999835707`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.541087754414083*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"-", "235.89859407100434`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "299.6997470371396`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "372.47374854200865`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "451.38227494426445`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "537.7811844548413`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "631.6769470046786`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "733.0708362036398`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "841.9633532235322`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "958.3547754407872`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1082.2452925825676`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.541087754414083*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", "37.661766842040315`", "}"}], ",", RowBox[{"{", "45.29589412750536`", "}"}], ",", RowBox[{"{", "54.18065980848532`", "}"}], ",", RowBox[{"{", "61.612355421410996`", "}"}], ",", RowBox[{"{", "69.02816392355183`", "}"}], ",", RowBox[{"{", "76.43628355468532`", "}"}], ",", RowBox[{"{", "83.83894747656528`", "}"}], ",", RowBox[{"{", "91.2373239025058`", "}"}], ",", RowBox[{"{", "98.63219437350267`", "}"}], ",", RowBox[{"{", "106.02413683801737`", "}"}]}], "}"}]], "Output", CellChangeTimes->{3.541087754414083*^9}] }, Open ]] }, WindowSize->{663, 815}, WindowMargins->{{451, Automatic}, {Automatic, 20}}, 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, 11233, 331, 632, "Input"], Cell[11815, 355, 683, 13, 73, "Message"], Cell[12501, 370, 3272, 84, 202, "Output"], Cell[15776, 456, 815, 23, 30, "Output"], Cell[16594, 481, 809, 23, 50, "Output"], Cell[17406, 506, 607, 13, 50, "Output"] }, Open ]] } ] *) (* End of internal cache information *)