(* 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[ 18239, 538] NotebookOptionsPosition[ 17675, 517] NotebookOutlinePosition[ 18098, 533] CellTagsIndexPosition[ 18055, 530] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"Ne", "-", RowBox[{"Atom", " ", "and", " ", "homologues", " ", "1", "s", RowBox[{"(", "2", ")"}], "2", "s", RowBox[{"(", "2", ")"}], " ", "2", "p", RowBox[{"(", "6", ")"}]}]}], " ", ",", " ", 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[{"12", "/", RowBox[{"R", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vne", "=", RowBox[{ RowBox[{ RowBox[{"-", "3"}], "*", RowBox[{"z", "/", "P"}]}], "-", RowBox[{"3", "*", RowBox[{"z", "/", "Q"}]}], "-", RowBox[{"6", "*", 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[{"12", "*", 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[{"4.5", "*", RowBox[{"sig", "/", "R"}]}], "+", RowBox[{"12", "*", 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[{"6.582336", "/", "R"}]}]}], ";"}], " ", RowBox[{"(*", " ", RowBox[{"8.48528", "/", "R"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"func", "=", RowBox[{"T", "+", "Vne", "+", "Vee"}]}], ";", " ", RowBox[{"Vir", " ", "=", " ", RowBox[{ RowBox[{"(", RowBox[{"Vne", "+", "Vee"}], ")"}], "/", "T"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"Etot", ",", " ", "P", ",", " ", "Q"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"t", " ", "=", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{"FindMinimum", "[", RowBox[{"func", ",", RowBox[{"{", RowBox[{"P", ",", "0.15"}], "}"}], ",", RowBox[{"{", RowBox[{"Q", ",", "1.6"}], "}"}], ",", RowBox[{"{", RowBox[{"R", ",", "0.6"}], "}"}], ",", RowBox[{"{", RowBox[{"MaxIterations", "\[Rule]", "800"}], "}"}], ",", RowBox[{"{", RowBox[{"Method", "\[Rule]", "\"\\""}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"z", ",", "10", ",", "18"}], "}"}]}], "]"}]}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"Virial", " ", "ratio"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"N", "[", RowBox[{ RowBox[{"Vir", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "9"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", "10"}], "]"}], ",", RowBox[{"{", RowBox[{"z", " ", ",", "10", ",", "18"}], "}"}]}], "]"}], "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"N", "[", RowBox[{ RowBox[{"Vne", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "9"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", "10"}], "]"}], ",", RowBox[{"{", RowBox[{"z", " ", ",", "10", ",", "18"}], "}"}]}], "]"}], "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"N", "[", RowBox[{ RowBox[{"Vee", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "9"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", "10"}], "]"}], ",", RowBox[{"{", RowBox[{"z", " ", ",", "10", ",", "18"}], "}"}]}], "]"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"Etot", " ", "in", " ", "eV"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"N", "[", RowBox[{ RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "9"}], "}"}], ",", "1"}], "]"}], "]"}], "*", "27.21138"}], ",", "10"}], "]"}], ",", RowBox[{"{", RowBox[{"z", ",", "10", ",", "18"}], "}"}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"PA", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "9"}], "}"}], ",", "2", ",", "1"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"z", ",", "10", ",", "18"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"QA", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "9"}], "}"}], ",", "2", ",", "2"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"z", ",", "10", ",", "18"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"RA", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "9"}], "}"}], ",", "2", ",", "3"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"z", ",", "10", ",", "18"}], "}"}]}], "]"}]}], ";"}]}]}]], "Input", CellChangeTimes->{{3.565285708771845*^9, 3.565285709801447*^9}, { 3.5652857568979297`*^9, 3.5652857982848024`*^9}, {3.565285856628905*^9, 3.5652858574401064`*^9}, {3.565285894817772*^9, 3.56528589956018*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "129.90549560209547`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.15398421744914842`"}], ",", RowBox[{"Q", "\[Rule]", "1.5885139078383728`"}], ",", RowBox[{"R", "\[Rule]", "0.6177358807009392`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "162.97478956201732`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.13969192675004327`"}], ",", RowBox[{"Q", "\[Rule]", "1.2913825954113067`"}], ",", RowBox[{"R", "\[Rule]", "0.5327797060592275`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "200.04569312208884`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.12782468311414308`"}], ",", RowBox[{"Q", "\[Rule]", "1.02443385856661`"}], ",", RowBox[{"R", "\[Rule]", "0.4724387728144514`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "241.2194433693671`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.11781805218407913`"}], ",", RowBox[{"Q", "\[Rule]", "0.74977754076882`"}], ",", RowBox[{"R", "\[Rule]", "0.4360242734454343`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "286.45459716903`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.1092617149883094`"}], ",", RowBox[{"Q", "\[Rule]", "0.6466370326233891`"}], ",", RowBox[{"R", "\[Rule]", "0.39616563882769434`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "335.6994761497142`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.10186175056008497`"}], ",", RowBox[{"Q", "\[Rule]", "0.5772654286442558`"}], ",", RowBox[{"R", "\[Rule]", "0.3614697967847427`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "388.9487782394322`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.0953990268490955`"}], ",", RowBox[{"Q", "\[Rule]", "0.5235313645874632`"}], ",", RowBox[{"R", "\[Rule]", "0.33197443067159177`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "446.20084439126464`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.0897062911586054`"}], ",", RowBox[{"Q", "\[Rule]", "0.47979765093440596`"}], ",", RowBox[{"R", "\[Rule]", "0.3067791621916184`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "507.45488567684595`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.08465380331289858`"}], ",", RowBox[{"Q", "\[Rule]", "0.44320969672752136`"}], ",", RowBox[{"R", "\[Rule]", "0.28506814451379575`"}]}], "}"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.5652859287946315`*^9, 3.56532867906854*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"-", "1.9999999999999964`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999999630884`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999999999936`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999954166932`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999999806777`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999999904157`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999999998836`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999999964242`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999999988225`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.5652859287946315`*^9, 3.56532867906854*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"-", "310.83962654506786`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "385.6667237935207`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "469.17780773429934`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "561.9234344188847`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "661.3820931128294`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "768.7124464586709`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "884.0137357751692`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1007.303537477298`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1138.5874621119622`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.5652859287946315`*^9, 3.56532867906854*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", "51.02863534087648`", "}"}], ",", RowBox[{"{", "59.71714466347031`", "}"}], ",", RowBox[{"{", "69.08642149012044`", "}"}], ",", RowBox[{"{", "79.48454657456786`", "}"}], ",", RowBox[{"{", "88.4728987692344`", "}"}], ",", RowBox[{"{", "97.31349415602482`", "}"}], ",", RowBox[{"{", "106.11617929625945`", "}"}], ",", RowBox[{"{", "114.90184869317312`", "}"}], ",", RowBox[{"{", "123.67769075767262`", "}"}]}], "}"}]], "Output", CellChangeTimes->{3.5652859287946315`*^9, 3.56532867906854*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"Do", "[", RowBox[{ RowBox[{"Show", "[", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"P", "/.", RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}]}]}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"Q", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}]}]}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Dashing", "[", RowBox[{"{", RowBox[{"0.02", ",", "0.02"}], "}"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{ RowBox[{"-", "R"}], "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}]}]}], "}"}], ",", RowBox[{"R", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}]}]}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "R"}], "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}]}], ",", "0"}], "}"}], ",", RowBox[{"R", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}]}]}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"R", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}]}], ",", "0"}], "}"}], ",", RowBox[{"R", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}]}]}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{"R", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}]}]}], "}"}], ",", RowBox[{"R", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}]}]}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"R", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}]}]}], "]"}]}], "}"}]}], "}"}], "]"}], ",", 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", ",", "9"}], "}"}]}], "]"}], "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{{3.565285708771845*^9, 3.565285709801447*^9}, { 3.5652857568979297`*^9, 3.5652857982848024`*^9}, {3.565285856628905*^9, 3.5652858574401064`*^9}, 3.565285894817772*^9, {3.5653285198234606`*^9, 3.5653286510352907`*^9}}] }, WindowSize->{663, 815}, WindowMargins->{{57, Automatic}, {Automatic, 40}}, 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, 7737, 238, 472, "Input"], Cell[8319, 262, 2983, 78, 192, "Output"], Cell[11305, 342, 766, 21, 31, "Output"], Cell[12074, 365, 760, 21, 52, "Output"], Cell[12837, 388, 579, 12, 52, "Output"] }, Open ]], Cell[13431, 403, 4240, 112, 232, "Input"] } ] *) (* End of internal cache information *)