(* 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[     16225,        483]
NotebookOptionsPosition[     15784,        465]
NotebookOutlinePosition[     16207,        481]
CellTagsIndexPosition[     16164,        478]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"C", "-", 
     RowBox[{"Atom", " ", "and", " ", "homologues", " ", "1", "s", 
      RowBox[{"(", "2", ")"}], "2", "s", 
      RowBox[{"(", "2", ")"}], " ", "2", "p", 
      RowBox[{"(", "2", ")"}], " ", "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[{"4", "*", 
          RowBox[{"Q", "^", "2"}]}], ")"}]}]}], "+", 
      RowBox[{"4", "/", 
       RowBox[{"R", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Vne", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"-", "3"}], "*", 
       RowBox[{"z", "/", "P"}]}], "-", 
      RowBox[{"3", "*", 
       RowBox[{"z", "/", "Q"}]}], "-", 
      RowBox[{"2", "*", 
       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[{"4", "*", 
       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", "*", 
       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[{"0.548528", "/", "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.27"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"Q", ",", "1.6"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"R", ",", "2.0"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"MaxIterations", "\[Rule]", "200"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"Method", "\[Rule]", "Newton"}], "}"}]}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"z", ",", "6", ",", "18"}], "}"}]}], "]"}]}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"Virial", " ", "ratio"}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"N", "[", 
       RowBox[{
        RowBox[{"Vir", " ", "/.", " ", 
         RowBox[{"t", "[", 
          RowBox[{"[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"z", "-", "5"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", 
        "10"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"z", " ", ",", "6", ",", "18"}], "}"}]}], "]"}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"N", "[", 
       RowBox[{
        RowBox[{"Vne", " ", "/.", " ", 
         RowBox[{"t", "[", 
          RowBox[{"[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"z", "-", "5"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", 
        "10"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"z", " ", ",", "6", ",", "18"}], "}"}]}], "]"}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"N", "[", 
       RowBox[{
        RowBox[{"Vee", " ", "/.", " ", 
         RowBox[{"t", "[", 
          RowBox[{"[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"z", "-", "5"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", 
        "10"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"z", " ", ",", "6", ",", "18"}], "}"}]}], "]"}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"Etot", " ", "in", " ", "eV"}], " ", "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"N", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"t", "[", 
          RowBox[{"[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"z", "-", "5"}], "}"}], ",", "1"}], "]"}], "]"}], "*", 
         "27.21138"}], ",", "10"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"z", ",", "6", ",", "18"}], "}"}]}], "]"}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"PA", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"t", "[", 
        RowBox[{"[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"z", "-", "5"}], "}"}], ",", "2", ",", "1"}], "]"}], "]"}],
        ",", 
       RowBox[{"{", 
        RowBox[{"z", ",", "6", ",", "18"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"QA", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"t", "[", 
        RowBox[{"[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"z", "-", "5"}], "}"}], ",", "2", ",", "2"}], "]"}], "]"}],
        ",", 
       RowBox[{"{", 
        RowBox[{"z", ",", "6", ",", "18"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"RA", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"t", "[", 
        RowBox[{"[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"z", "-", "5"}], "}"}], ",", "2", ",", "3"}], "]"}], "]"}],
        ",", 
       RowBox[{"{", 
        RowBox[{"z", ",", "6", ",", "18"}], "}"}]}], "]"}]}], 
    ";"}]}]}]], "Input",
 CellChangeTimes->{{3.5648433523441534`*^9, 3.564843391640623*^9}, {
  3.5648434938676023`*^9, 3.5648435191396465`*^9}, {3.564843588232168*^9, 
  3.5648435965313826`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "37.54762137638885`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.263204963166588`"}], ",", 
       RowBox[{"Q", "\[Rule]", "1.5414013734294478`"}], ",", 
       RowBox[{"R", "\[Rule]", "1.7669172600613985`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "53.552327924755865`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.22400317124036692`"}], ",", 
       RowBox[{"Q", "\[Rule]", "1.2707498945369435`"}], ",", 
       RowBox[{"R", "\[Rule]", "1.1290222850583227`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "72.58727951672662`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.19496664985718593`"}], ",", 
       RowBox[{"Q", "\[Rule]", "1.0674794873826865`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.8530105404224914`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "94.63366763304896`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.17259028185749714`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.9137104906388269`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.6934491837746999`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "119.68548274326108`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.15481709896906395`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.7962091107435088`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.5868435018278179`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "147.74042292026166`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.14035921537294688`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.7045048134963547`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.5096781015363541`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "178.79739644991557`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.1283684565090239`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.6313080973311753`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.45090754663420296`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "212.85579837636027`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.11826327169286124`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.5716779463726016`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.40451949311168683`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "249.91525616140422`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.10963155236196531`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.5222289576778898`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.36691159096276427`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "289.9755231034424`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.10217310010946866`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.4805911196355063`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.33577546978536915`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "333.03642711894736`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.09566405451640556`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.4450659780475779`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.30955616463884156`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "379.09784343264187`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.08993408093699519`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.4144088478503389`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.287164631228629`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "428.1596788168225`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.08485127778914286`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.38768868409557133`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.2678138723489608`"}]}], "}"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.564843601788592*^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[{"PA", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}], ",", 
         RowBox[{"Circle", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"0", ",", "0"}], "}"}], ",", 
           RowBox[{"Q", " ", "/.", " ", 
            RowBox[{"QA", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"Dashing", "[", 
            RowBox[{"{", 
             RowBox[{"0.02", ",", "0.02"}], "}"}], "]"}], ",", 
           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[{"{", 
       RowBox[{
        RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", 
        RowBox[{"Axes", " ", "->", " ", "True"}], ",", 
        RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", 
        RowBox[{"PlotRange", " ", "\[Rule]", " ", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"-", "3.6"}], ",", "3.6"}], "}"}], ",", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"-", "3.6"}], ",", "3.6"}], "}"}]}], "}"}]}], ",", " ", 
        RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "13"}], "}"}]}], "]"}], 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.5648433523441534`*^9, 3.564843391640623*^9}, {
   3.5648434938676023`*^9, 3.5648435191396465`*^9}, 3.564843588232168*^9}]
},
WindowSize->{734, 815},
WindowMargins->{{214, Automatic}, {Automatic, 7}},
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, 7489, 233, 432, "Input"],
Cell[8071, 257, 4249, 113, 272, "Output"]
}, Open  ]],
Cell[12335, 373, 3445, 90, 192, "Input"]
}
]
*)

(* End of internal cache information *)