(* 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[     13029,        399]
NotebookOptionsPosition[     12467,        378]
NotebookOutlinePosition[     12891,        394]
CellTagsIndexPosition[     12848,        391]
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[{"6.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", ",", "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[{
    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[{"Table", "[", 
    RowBox[{
     RowBox[{"N", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"t", "[", 
         RowBox[{"[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"z", "-", "9"}], "}"}], ",", "1"}], "]"}], "]"}], "*", 
        "27.21138"}], ",", "10"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"z", ",", "10", ",", "18"}], "}"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.566231343910093*^9, 3.5662313445964937`*^9}, 
   3.566231395717784*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "130.472923715349`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.1539519882592071`"}], ",", 
       RowBox[{"Q", "\[Rule]", "1.5021631064446892`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.6091890012608211`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "163.60496596248268`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.1396446133015732`"}], ",", 
       RowBox[{"Q", "\[Rule]", "1.1694915878086667`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.5323126162033147`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "200.81036198325643`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.1278166417476058`"}], ",", 
       RowBox[{"Q", "\[Rule]", "1.0213424140990899`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.46822262918739344`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "242.0497567467124`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.11787158702926995`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.7827227493919515`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.4294169460993777`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "287.35683133122245`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.10914150179699829`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.6808021945935068`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.3904917456762293`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "336.69007146478504`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.10178517434979706`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.5984136395572477`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.35723022655716885`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "390.0343836911412`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.09538550355656433`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.5321361898385282`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.3286833448331945`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "447.3753527312848`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.08971757005194607`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.48673224658307723`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.3046823485335471`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "508.715840930971`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.08467803833964592`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.45077832728826267`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.2832470060766872`"}]}], "}"}]}], "}"}]}], 
  "}"}]], "Output"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"-", "1.994036628913992`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.99570461713758`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.9986212368073328`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.0014046362798914`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.0012364441674286`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.0010303240522145`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.9996858059150486`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.001320220051573`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.0020906183255605`"}], "}"}]}], "}"}]], "Output"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"-", "313.3287481367081`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "388.5188458572932`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "470.67419663434964`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "562.3362565771624`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "661.6268425755761`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "769.2447963197454`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "885.4980213583713`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1008.0058192368742`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1138.7950730581877`"}], "}"}]}], "}"}]], "Output"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", "51.6001745587442`", "}"}], ",", 
   RowBox[{"{", "60.60313638052499`", "}"}], ",", 
   RowBox[{"{", "68.77622046689837`", "}"}], ",", 
   RowBox[{"{", "78.57625805855152`", "}"}], ",", 
   RowBox[{"{", "87.26804182426098`", "}"}], ",", 
   RowBox[{"{", "96.2111962175376`", "}"}], ",", 
   RowBox[{"{", "105.30666896431984`", "}"}], ",", 
   RowBox[{"{", "113.8449689469338`", "}"}], ",", 
   RowBox[{"{", "122.42470305777194`", "}"}]}], "}"}]], "Output"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"-", "3550.3483069293734`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "4451.916898692181`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "5464.327067863944`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "6586.507909742355`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "7819.3759329498`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "9161.801476855422`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "10613.373827685446`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "12173.700725805029`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "13842.860059592205`"}], "}"}]}], "}"}]], "Output"]
}, Open  ]]
},
WindowSize->{663, 815},
WindowMargins->{{Automatic, 540}, {Automatic, 10}},
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, 6339, 195, 392, "Input"],
Cell[6921, 219, 2920, 77, 192, "Output"],
Cell[9844, 298, 697, 20, 52, "Output"],
Cell[10544, 320, 695, 20, 52, "Output"],
Cell[11242, 342, 512, 11, 52, "Output"],
Cell[11757, 355, 694, 20, 52, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)