(* 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[     13652,        403]
NotebookOptionsPosition[     13253,        386]
NotebookOutlinePosition[     13677,        402]
CellTagsIndexPosition[     13634,        399]
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"}], ";", " ", 
    RowBox[{"n", "=", "9"}], ";"}], "\[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[{"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.17"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"Q", ",", "1.78"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"R", ",", "0.7"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"MaxIterations", "\[Rule]", "800"}], ",", " ", 
          RowBox[{"Method", "\[Rule]", " ", "\"\<Gradient\>\""}]}], "}"}]}], 
       "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"z", ",", "9", ",", "n"}], "}"}]}], "]"}]}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"Virial", " ", "ratio"}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"N", "[", 
       RowBox[{
        RowBox[{"Vir", " ", "/.", " ", 
         RowBox[{"t", "[", 
          RowBox[{"[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"z", "-", "8"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", 
        "10"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"z", " ", ",", "9", ",", "n"}], "}"}]}], "]"}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"N", "[", 
       RowBox[{
        RowBox[{"Vne", " ", "/.", " ", 
         RowBox[{"t", "[", 
          RowBox[{"[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"z", "-", "8"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", 
        "10"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"z", " ", ",", "9", ",", "n"}], "}"}]}], "]"}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"N", "[", 
       RowBox[{
        RowBox[{"Vee", " ", "/.", " ", 
         RowBox[{"t", "[", 
          RowBox[{"[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"z", "-", "8"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", 
        "10"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"z", " ", ",", "9", ",", "n"}], "}"}]}], "]"}], ";"}], 
   "\[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", ",", "n"}], "}"}]}], "]"}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"PA", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"t", "[", 
        RowBox[{"[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"z", "-", "8"}], "}"}], ",", "2", ",", "1"}], "]"}], "]"}],
        ",", 
       RowBox[{"{", 
        RowBox[{"z", ",", "9", ",", "n"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"QA", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"t", "[", 
        RowBox[{"[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"z", "-", "8"}], "}"}], ",", "2", ",", "2"}], "]"}], "]"}],
        ",", 
       RowBox[{"{", 
        RowBox[{"z", ",", "9", ",", "n"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"RA", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"t", "[", 
        RowBox[{"[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"z", "-", "8"}], "}"}], ",", "2", ",", "3"}], "]"}], "]"}],
        ",", 
       RowBox[{"{", 
        RowBox[{"z", ",", "9", ",", "n"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   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", ",", 
       RowBox[{"n", "-", "9"}]}], "}"}]}], "]"}], 
   "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.564848410426808*^9, 3.5648484320172453`*^9}, {
  3.5648484779281263`*^9, 3.564848590903525*^9}, {3.5648487364985805`*^9, 
  3.564848751459007*^9}, {3.5648489799838085`*^9, 3.5648489835094147`*^9}, {
  3.56484904355392*^9, 3.564849046689525*^9}, {3.5648491026780233`*^9, 
  3.564849109323635*^9}, {3.564849151506109*^9, 3.5648492156066217`*^9}, {
  3.564849352840063*^9, 3.564849393743335*^9}, {3.564849429170997*^9, 
  3.564849460464652*^9}, {3.564849522427961*^9, 3.5648495884784765`*^9}, {
  3.56484964700978*^9, 3.564849673732627*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"-", "99.1184135543036`"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"P", "\[Rule]", "0.1715184559693351`"}], ",", 
      RowBox[{"Q", "\[Rule]", "1.780188025548774`"}], ",", 
      RowBox[{"R", "\[Rule]", "0.7107674371199213`"}]}], "}"}]}], "}"}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.5648484490212755`*^9, 3.564848486086941*^9, 3.5648485311866198`*^9, 
   3.5648485985319386`*^9, 3.564848765889032*^9, 3.5648489906854267`*^9, 
   3.5648490552383404`*^9, 3.564849116655648*^9, 3.5648491633777304`*^9, {
   3.564849193922584*^9, 3.5648492221742334`*^9}, 3.5648493660376863`*^9, 
   3.5648494009661474`*^9, 3.564849438047413*^9, 3.564849473475075*^9, 
   3.5648496058725076`*^9, {3.5648496532029905`*^9, 3.5648496805810385`*^9}}]
}, Open  ]]
},
WindowSize->{869, 815},
WindowMargins->{{Automatic, 509}, {Automatic, 44}},
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, 11839, 342, 612, "Input"],
Cell[12421, 366, 816, 17, 31, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)