(* 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[     16078,        470]
NotebookOptionsPosition[     15508,        448]
NotebookOutlinePosition[     15935,        465]
CellTagsIndexPosition[     15892,        462]
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.694975", "/", "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[{"t", " ", "=", " ", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"FindMinimum", "[", 
       RowBox[{"func", ",", 
        RowBox[{"{", 
         RowBox[{"P", ",", "0.22"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"Q", ",", "1.0"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"R", ",", "1.1"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"MaxIterations", "\[Rule]", "120"}], "}"}]}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"z", ",", "6", ",", "18"}], "}"}]}], "]"}]}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"Virial", " ", "ratio"}], " ", "*)"}], "\[IndentingNewLine]", 
   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[{"Table", "[", 
    RowBox[{
     RowBox[{"N", "[", 
      RowBox[{
       RowBox[{"Vne", " ", "/.", " ", 
        RowBox[{"t", "[", 
         RowBox[{"[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"z", "-", "5"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", 
       "10"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"z", " ", ",", "6", ",", "18"}], "}"}]}], "]"}], 
   "\[IndentingNewLine]", 
   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[{"Table", "[", 
    RowBox[{
     RowBox[{"N", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"t", "[", 
         RowBox[{"[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"z", "-", "5"}], "}"}], ",", "1"}], "]"}], "]"}], "*", 
        "27.21138"}], ",", "10"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"z", ",", "6", ",", "18"}], "}"}]}], "]"}], 
   "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.5410885842887406`*^9, 3.541088596550362*^9}, {
  3.541088652819661*^9, 3.541088678575306*^9}, {3.5410887162649727`*^9, 
  3.5410887354998064`*^9}, {3.5410887682754636`*^9, 3.5410887727370715`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "37.46718700172893`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.2632072695483324`"}], ",", 
       RowBox[{"Q", "\[Rule]", "1.5281837081918972`"}], ",", 
       RowBox[{"R", "\[Rule]", "1.8768334635575599`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "53.42500872397746`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.22400254144931234`"}], ",", 
       RowBox[{"Q", "\[Rule]", "1.2596701700962223`"}], ",", 
       RowBox[{"R", "\[Rule]", "1.1720246678940636`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "72.4177535111813`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.1949644343280144`"}], ",", 
       RowBox[{"Q", "\[Rule]", "1.060700380386825`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.8749428714237419`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "94.42445125996912`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.17258773846426234`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.9096787937186024`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.7066204596175716`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "119.43777604227729`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.15481471128444707`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.7936979900474184`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.5956412851727786`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "147.45484735250648`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.14035711694396563`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.7028549888984699`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.515986657018611`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "178.4743081729923`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.12836665599167651`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.6301721321612119`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.45566189669147067`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "212.49541935025294`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.11826173796095947`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.5708646603872615`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.40823618500009545`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "249.5177340374536`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.10963024722321929`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.5216275680289877`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.3698997836273564`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "289.5409612160212`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.10217198605687451`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.48013428024993227`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.33823188027508655`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "332.56490074461027`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.09566309935622641`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.44471099512958023`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.3116121198355315`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "378.58940920396304`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.0899332576525295`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.41412761976258966`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.2889112818334707`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "427.6143804774629`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.0848505645126518`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.3874621538068578`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.2693165198542921`"}]}], "}"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.541088565303507*^9, 3.54108860116797*^9, {3.5410886616336765`*^9, 
   3.541088684269316*^9}, 3.541088739930214*^9, 3.541088778696282*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"-", "1.9999999980833432`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.000000001177959`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.0000000002916307`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.0000000052606897`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.000000001183142`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.9999999999850904`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.000000001637582`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.9999999948019578`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.00000000068918`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.000000002605458`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.000000004956512`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.9999999997460671`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.000000000020479`"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.541088565303507*^9, 3.54108860116797*^9, {3.5410886616336765`*^9, 
   3.541088684269316*^9}, 3.541088739930214*^9, 3.541088778696282*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"-", "86.55960538310434`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "122.36510805972092`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "164.01284241447644`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "211.59633811142822`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "265.15504229808056`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "324.7026720596104`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "390.24454090304613`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "461.7830306724663`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "539.3193450976307`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "622.8541665598946`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "712.3879069786078`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "807.9208371259058`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "909.4531341939205`"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.541088565303507*^9, 3.54108860116797*^9, {3.5410886616336765`*^9, 
   3.541088684269316*^9}, 3.541088739930214*^9, 3.541088778696282*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", "11.62523130783475`", "}"}], ",", 
   RowBox[{"{", "15.515090674698447`", "}"}], ",", 
   RowBox[{"{", "19.17733541323309`", "}"}], ",", 
   RowBox[{"{", "22.747436088227776`", "}"}], ",", 
   RowBox[{"{", "26.279490354837865`", "}"}], ",", 
   RowBox[{"{", "29.792977352399053`", "}"}], ",", 
   RowBox[{"{", "33.29592484932799`", "}"}], ",", 
   RowBox[{"{", "36.792190867400194`", "}"}], ",", 
   RowBox[{"{", "40.283877194685935`", "}"}], ",", 
   RowBox[{"{", "43.77224488223907`", "}"}], ",", 
   RowBox[{"{", "47.25810713774912`", "}"}], ",", 
   RowBox[{"{", "50.74201862184337`", "}"}], ",", 
   RowBox[{"{", "54.22437324775207`", "}"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.541088565303507*^9, 3.54108860116797*^9, {3.5410886616336765`*^9, 
   3.541088684269316*^9}, 3.541088739930214*^9, 3.541088778696282*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"-", "1019.5338630351065`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1453.7682138914658`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1970.5870095390885`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2569.4196245264984`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "3250.066710241303`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "4012.4498841510476`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "4856.532219932399`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "5782.293604199086`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "6789.721877632084`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "7878.809121214414`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "9049.549888823873`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "10301.940277824535`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "11635.977400636822`"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.541088565303507*^9, 3.54108860116797*^9, {3.5410886616336765`*^9, 
   3.541088684269316*^9}, 3.541088739930214*^9, 3.541088778696282*^9}]
}, Open  ]]
},
WindowSize->{604, 815},
WindowMargins->{{405, Automatic}, {Automatic, 3}},
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, 6228, 189, 392, "Input"],
Cell[6810, 213, 4355, 113, 259, "Output"],
Cell[11168, 328, 1143, 31, 50, "Output"],
Cell[12314, 361, 1143, 31, 69, "Output"],
Cell[13460, 394, 884, 18, 69, "Output"],
Cell[14347, 414, 1145, 31, 69, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)