(* 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[     52401,       1181]
NotebookOptionsPosition[     50403,       1110]
NotebookOutlinePosition[     50888,       1128]
CellTagsIndexPosition[     50845,       1125]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
C2H6 Molecule after Kimball, non overlapping, parameterized with MP2/6-31G(d) \
ES 04 June 2013\
\>", "Subsubtitle",
 CellChangeTimes->{
  3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, {
   3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 
   3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, {
   3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 
   3.578896407337781*^9}, 3.5789031386419754`*^9, {3.5789304190931296`*^9, 
   3.5789304212303333`*^9}, {3.5789317063449907`*^9, 
   3.5789317142698045`*^9}, {3.5792337731768165`*^9, 3.579233777451224*^9}, 
   3.5792341967955604`*^9, 3.579234656987569*^9, {3.579237662522448*^9, 
   3.579237690805298*^9}, {3.57924170619635*^9, 3.579241706555151*^9}, {
   3.579280081796549*^9, 3.57928008223335*^9}, {3.5793187338499007`*^9, 
   3.5793187469539237`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"Clear", "[", 
   RowBox[{
   "z", ",", "\[Sigma]", ",", "k", ",", "nc", ",", "R1", ",", "R2", ",", "R3",
     ",", "p", ",", "Ekin", ",", "Epot", ",", "vne", ",", "vnn", ",", "vee", 
    ",", "vir", ",", "ow", ",", "cn", ",", "wi", ",", "c2", ",", "ec", ",", 
    "\[IndentingNewLine]", "dwnr", ",", "dwr", ",", "dnr", ",", "d", ",", 
    "d1", ",", "d2", ",", "d3", ",", "d4", ",", "d5", ",", "d6", ",", "d7", 
    ",", "d8", ",", "d9", ",", "nn", ",", "n", ",", "pi", ",", "i", ",", "j", 
    ",", "t", ",", "\[Lambda]", ",", "cs", ",", "ss", ",", "wi"}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"bohr", "=", "0.529177"}], ";", "  ", 
   RowBox[{"rad", "=", "57.2957795"}], ";"}], "\[IndentingNewLine]", " ", 
  RowBox[{"(*", "Parameterlists", " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"k", "=", 
    RowBox[{"{", 
     RowBox[{
     "1.02669238", ",", "1.02669238", ",", "1.36467050", ",", "1.20088548", 
      ",", "1.20088548", ",", "1.20088548", ",", "1.20088548", ",", 
      "1.20088548", ",", "1.20088548"}], "}"}]}], ";"}], " "}], "\n", 
 RowBox[{
  RowBox[{"\[Sigma]", "=", 
   RowBox[{"{", 
    RowBox[{
    "0.30102681", ",", "0.30102681", ",", "0.34863843", ",", "0.34858036", 
     ",", "0.34858036", ",", "0.34858036", ",", "0.34858036", ",", 
     "0.34858036", ",", "0.34858036"}], "}"}]}], ";"}]}], "Input",
 CellChangeTimes->{
  3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, {
   3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 
   3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, {
   3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 
   3.578896407337781*^9}, {3.578897895970396*^9, 3.5788979299784555`*^9}, {
   3.578901683206219*^9, 3.5789017130022717`*^9}, 3.5789017480087333`*^9, 
   3.5789031559112062`*^9, {3.578930508668487*^9, 3.5789305512877617`*^9}, {
   3.578930588213027*^9, 3.5789305892270284`*^9}, {3.578930623203888*^9, 
   3.5789307850073724`*^9}, {3.5789308197642336`*^9, 
   3.5789308212306356`*^9}, {3.5789308916959596`*^9, 
   3.5789309168276043`*^9}, {3.578930968635295*^9, 3.578931053967445*^9}, {
   3.5789310902531085`*^9, 3.5789312758466344`*^9}, {3.57898059849006*^9, 
   3.57898062724091*^9}, {3.5789833698569264`*^9, 3.578983385831354*^9}, {
   3.578983821063319*^9, 3.5789838692362037`*^9}, {3.578984093892198*^9, 
   3.5789842193788185`*^9}, {3.578984259408489*^9, 3.5789842993913593`*^9}, 
   3.57898438517591*^9, {3.578984422584776*^9, 3.5789845278225603`*^9}, 
   3.578994102766264*^9, {3.578994136384323*^9, 3.5789944171536164`*^9}, {
   3.5789944872445393`*^9, 3.5789944880089407`*^9}, {3.578994526821809*^9, 
   3.5789946056019473`*^9}, {3.578994644555216*^9, 3.5789949234213057`*^9}, {
   3.5790013546913004`*^9, 3.5790013995725794`*^9}, {3.579001432987838*^9, 
   3.5790014462946615`*^9}, {3.5790031307327833`*^9, 3.579003231402541*^9}, {
   3.5790863790628405`*^9, 3.5790864565325766`*^9}, {3.579087975881403*^9, 
   3.5790880551589375`*^9}, {3.5790943224634066`*^9, 3.579094402850005*^9}, {
   3.5792337813044305`*^9, 3.5792337913664484`*^9}, {3.5792339896271963`*^9, 
   3.5792339924820013`*^9}, 3.5792355216502876`*^9, {3.579240085619104*^9, 
   3.579240093216317*^9}, {3.579240186317281*^9, 3.579240221245742*^9}, {
   3.5792416408634357`*^9, 3.5792416835139103`*^9}, 3.5792418843954635`*^9, {
   3.579254526481946*^9, 3.579254528197949*^9}, {3.579256065565049*^9, 
   3.5792562478981695`*^9}, {3.579256291921447*^9, 3.579256298769859*^9}, {
   3.5792570670400085`*^9, 3.5792571152284927`*^9}, {3.579257266377158*^9, 
   3.5792572671259594`*^9}, {3.5792573858077683`*^9, 
   3.5792574224522324`*^9}, {3.579258705944487*^9, 3.579258744117754*^9}, {
   3.5792587970174465`*^9, 3.5792588039126587`*^9}, 3.5792608573214664`*^9, {
   3.579260900939143*^9, 3.579260918489174*^9}, {3.579273252603754*^9, 
   3.5792733376239038`*^9}, 3.579277735474428*^9, {3.579278201057246*^9, 
   3.579278202508048*^9}, {3.5792790593863535`*^9, 3.5792791008356266`*^9}, {
   3.579279199115799*^9, 3.5792792082574153`*^9}, {3.5792797138855033`*^9, 
   3.5792797526515713`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Terms", " ", "analogous", " ", "to", " ", "CH4_mat"}], " ", 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"z", "=", "6.0"}], ";", " ", 
    RowBox[{"nc", "=", "9"}], ";", " ", 
    RowBox[{"nn", "=", "2"}], ";", 
    RowBox[{"(*", 
     RowBox[{
      RowBox[{"wi", "=", 
       RowBox[{"\[Pi]", "-", 
        RowBox[{"ArcCos", "[", 
         RowBox[{
          RowBox[{"-", "1"}], "/", "3"}], "]"}]}]}], ";"}], "*)"}], " ", 
    RowBox[{"cs", "=", 
     RowBox[{"Cos", "[", "wi", "]"}]}], ";", " ", 
    RowBox[{"ss", "=", 
     RowBox[{"Sin", "[", "wi", "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"dwr", "=", 
     RowBox[{"Table", "[", 
      RowBox[{"0", ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "nc"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"dnr", "=", "dwr"}], ";", 
    RowBox[{"dwnr", "=", "dwr"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"c2", "=", 
     RowBox[{"Table", "[", 
      RowBox[{"0", ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ow", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "2"}], ",", 
       RowBox[{"-", "2"}], ",", 
       RowBox[{"-", "2"}], ",", 
       RowBox[{"-", "2"}], ",", 
       RowBox[{"-", "2"}], ",", 
       RowBox[{"-", "2"}], ",", 
       RowBox[{"-", "2"}], ",", 
       RowBox[{"-", "2"}], ",", 
       RowBox[{"-", "2"}]}], "}"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"cn", "=", 
     RowBox[{"{", 
      RowBox[{
      "z", ",", "z", ",", "0", ",", "1", ",", "1", ",", "1", ",", "1", ",", 
       "1", ",", "1"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"c", "=", 
     RowBox[{"{", 
      RowBox[{
      "R1", ",", "R1", ",", "R2", ",", "R3", ",", "R3", ",", "R3", ",", "R3", 
       ",", "R3", ",", "R3"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ec", "=", 
     RowBox[{"{", 
      RowBox[{
      "0", ",", "0", ",", "0", ",", "p", ",", "p", ",", "p", ",", "p", ",", 
       "p", ",", "p"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"d1", "=", 
     RowBox[{"R1", "+", "R2"}]}], ";"}], " ", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"d2", "=", 
     RowBox[{"d1", "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"R1", "+", "R3"}], ")"}], "*", "cs"}]}]}], ";", "         ", 
    RowBox[{"d3", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"R1", "+", "R3"}], ")"}], "*", "ss"}]}], ";", "        ", 
    RowBox[{"d4", "=", " ", 
     RowBox[{"d3", "/", "2"}]}], ";", " ", 
    RowBox[{"d5", "=", " ", 
     RowBox[{"d4", "*", 
      RowBox[{"Sqrt", "[", "3", "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"d6", "=", 
     RowBox[{"d1", "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"R1", "+", "R3", "+", "p"}], ")"}], "*", "cs"}]}]}], ";", " ", 
    RowBox[{"d7", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"R1", "+", "R3", "+", "p"}], ")"}], "*", "ss"}]}], ";", " ", 
    RowBox[{"d8", "=", 
     RowBox[{"d7", "/", "2"}]}], ";", " ", 
    RowBox[{"d9", "=", 
     RowBox[{"d8", "*", 
      RowBox[{"Sqrt", "[", "3", "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"w", " ", "=", " ", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"d1", ",", " ", "0", ",", " ", "0"}], "}"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "d1"}], ",", " ", "0", ",", "0"}], "}"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", "\[IndentingNewLine]", 
       "   ", 
       RowBox[{"{", 
        RowBox[{"d2", ",", " ", "d3", ",", " ", "0"}], "}"}], ",", "    ", 
       RowBox[{"{", 
        RowBox[{"d2", ",", " ", 
         RowBox[{"-", "d4"}], ",", " ", "d5"}], "}"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{"d2", ",", " ", 
         RowBox[{"-", "d4"}], ",", " ", 
         RowBox[{"-", "d5"}]}], "}"}], ",", "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "d2"}], ",", 
         RowBox[{"-", "d3"}], ",", "0"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "d2"}], ",", "d4", ",", 
         RowBox[{"-", "d5"}]}], "}"}], ",", "   ", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "d2"}], ",", "d4", ",", "d5"}], "}"}]}], "}"}]}], ";"}],
    " ", "\n", 
   RowBox[{
    RowBox[{"n", " ", "=", " ", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"d1", ",", " ", "0", ",", " ", "0"}], "}"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "d1"}], ",", " ", "0", ",", "0"}], "}"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", "\[IndentingNewLine]", 
       "   ", 
       RowBox[{"{", 
        RowBox[{"d6", ",", " ", "d7", ",", " ", "0"}], "}"}], ",", "   ", 
       RowBox[{"{", 
        RowBox[{"d6", ",", " ", 
         RowBox[{"-", "d8"}], ",", " ", "d9"}], "}"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{"d6", ",", " ", 
         RowBox[{"-", "d8"}], ",", " ", 
         RowBox[{"-", "d9"}]}], "}"}], ",", "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "d6"}], ",", 
         RowBox[{"-", "d7"}], ",", "0"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "d6"}], ",", "d8", ",", 
         RowBox[{"-", "d9"}]}], "}"}], ",", "  ", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "d6"}], ",", "d8", ",", "d9"}], "}"}]}], "}"}]}], ";"}],
    " ", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"i", " ", "=", " ", "1"}], ",", " ", 
      RowBox[{"i", " ", "\[LessEqual]", " ", "nc"}], ",", " ", 
      RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{
        RowBox[{"c2", "[", 
         RowBox[{"[", "i", "]"}], "]"}], " ", "=", " ", 
        RowBox[{
         RowBox[{"k", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "/", 
         RowBox[{
          RowBox[{"c", "[", 
           RowBox[{"[", "i", "]"}], "]"}], "^", "2"}]}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"For", "[", 
        RowBox[{
         RowBox[{"j", " ", "=", " ", "1"}], ",", " ", 
         RowBox[{"j", " ", "\[LessEqual]", " ", "nc"}], ",", " ", 
         RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"d", " ", "=", " ", 
           RowBox[{"Norm", "[", 
            RowBox[{
             RowBox[{"w", "[", 
              RowBox[{"[", "i", "]"}], "]"}], " ", "-", " ", 
             RowBox[{"w", "[", 
              RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"dwr", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", 
           RowBox[{"If", "[", 
            RowBox[{
             RowBox[{"i", " ", "\[NotEqual]", " ", "j"}], ",", " ", 
             RowBox[{"1", "/", "d"}], ",", 
             RowBox[{"1.5", "*", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"\[Sigma]", "[", 
                 RowBox[{"[", "i", "]"}], "]"}], "/", 
                RowBox[{"c", "[", 
                 RowBox[{"[", "i", "]"}], "]"}]}], ")"}]}]}], "]"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"d", "=", " ", 
           RowBox[{"Norm", "[", 
            RowBox[{
             RowBox[{"n", "[", 
              RowBox[{"[", "i", "]"}], "]"}], " ", "-", " ", 
             RowBox[{"n", "[", 
              RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"dnr", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", 
           RowBox[{"If", "[", 
            RowBox[{
             RowBox[{"i", " ", "\[NotEqual]", " ", "j"}], ",", " ", 
             RowBox[{"1", "/", "d"}], ",", " ", "0"}], "]"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"d", "=", " ", 
           RowBox[{"Norm", "[", 
            RowBox[{
             RowBox[{"w", "[", 
              RowBox[{"[", "i", "]"}], "]"}], " ", "-", " ", 
             RowBox[{"n", "[", 
              RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"dwnr", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", 
           RowBox[{"If", "[", 
            RowBox[{
             RowBox[{"i", " ", "\[NotEqual]", " ", "j"}], ",", " ", 
             RowBox[{"1", "/", "d"}], ",", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"1.5", " ", "-", " ", 
                RowBox[{"0.5", "*", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"ec", "[", 
                    RowBox[{"[", "i", "]"}], "]"}], "/", 
                    RowBox[{"c", "[", 
                    RowBox[{"[", "i", "]"}], "]"}]}], ")"}], "^", "2"}]}]}], 
               ")"}], "/", 
              RowBox[{"c", "[", 
               RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}]}]}]}], "]"}]}]}], 
     "]"}], ";"}], " ", "\n", 
   RowBox[{
    RowBox[{"Ekin", " ", "=", " ", 
     RowBox[{
      RowBox[{"-", 
       RowBox[{"(", 
        RowBox[{"9", "/", "8"}], ")"}]}], "*", 
      RowBox[{"ow", " ", ".", " ", "c2"}]}]}], ";"}], " ", "\n", 
   RowBox[{
    RowBox[{"vee", " ", "=", " ", 
     RowBox[{"0.5", "*", 
      RowBox[{"ow", " ", ".", " ", 
       RowBox[{"(", 
        RowBox[{"ow", " ", ".", " ", "dwr"}], ")"}]}]}]}], ";"}], " ", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"vnn", " ", "=", " ", 
     RowBox[{"0.5", "*", 
      RowBox[{"cn", " ", ".", " ", 
       RowBox[{"(", 
        RowBox[{"cn", " ", ".", " ", "dnr"}], ")"}]}]}]}], ";"}], " ", "\n", 
   RowBox[{
    RowBox[{"vne", " ", "=", " ", 
     RowBox[{"cn", " ", ".", " ", 
      RowBox[{"(", 
       RowBox[{"ow", " ", ".", " ", "dwnr"}], ")"}]}]}], ";"}], " ", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Epot", " ", "=", " ", 
     RowBox[{"vne", " ", "+", " ", "vee", " ", "+", " ", "vnn"}]}], ";"}], 
   " ", "\n", 
   RowBox[{
    RowBox[{"func", " ", "=", " ", 
     RowBox[{"Ekin", " ", "+", " ", "Epot"}]}], ";"}], " "}]}]], "Text", \
"Input",
 CellChangeTimes->{{3.5792338519413548`*^9, 3.5792338963078327`*^9}, {
   3.579233998924813*^9, 3.579234002029218*^9}, 3.579235530823104*^9, {
   3.579254579459639*^9, 3.579254604497683*^9}, {3.5792546701893983`*^9, 
   3.5792549883675575`*^9}, {3.579255040643249*^9, 3.5792550505492663`*^9}, {
   3.579255231415984*^9, 3.5792553234093456`*^9}, {3.579255363563816*^9, 
   3.5792558209098196`*^9}, {3.5792558541534777`*^9, 3.579255854403078*^9}, {
   3.5792559106099772`*^9, 3.57925596911008*^9}, 3.5792560056141443`*^9, {
   3.5792566762281218`*^9, 3.579256679098527*^9}, 3.5792567540254583`*^9, {
   3.5792579024650755`*^9, 3.5792579089702873`*^9}, {3.579258007921261*^9, 
   3.57925801311607*^9}, {3.5792585291025763`*^9, 3.579258546933408*^9}, {
   3.5792589863081794`*^9, 3.579258990551387*^9}, {3.57925925426985*^9, 
   3.579259268918276*^9}, {3.5792608668842835`*^9, 3.579260889363923*^9}, {
   3.579260931390397*^9, 3.579260934323202*^9}, {3.579265466119326*^9, 
   3.5792654959777784`*^9}, {3.5792725895401897`*^9, 3.57927260670022*^9}, {
   3.579273860661622*^9, 3.579273922609331*^9}, {3.579274366742111*^9, 
   3.579274384229742*^9}, {3.5792744185342026`*^9, 3.579274459203474*^9}, {
   3.579274876145406*^9, 3.5792748782202096`*^9}, {3.579274985329998*^9, 
   3.579274996889618*^9}, {3.5792755145297275`*^9, 3.579275563435813*^9}, {
   3.579275828589479*^9, 3.5792758995072036`*^9}, {3.5792759677573233`*^9, 
   3.5792759978809767`*^9}, {3.5792761061763663`*^9, 
   3.5792761209339924`*^9}, {3.5792761561432543`*^9, 3.579276190619315*^9}, {
   3.5792762354381933`*^9, 3.5792762371853967`*^9}, {3.5792762674650497`*^9, 
   3.579276268151451*^9}, {3.579277193841477*^9, 3.5792772111887074`*^9}, {
   3.579277403396645*^9, 3.579277409106255*^9}, {3.5792775332200727`*^9, 
   3.579277607897404*^9}, 3.579277692855153*^9, 3.5792779685076375`*^9, 
   3.579279041992323*^9, {3.5792792325310574`*^9, 3.5792792443246784`*^9}, {
   3.579279484814701*^9, 3.579279489479109*^9}, 3.5792796121265244`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"vir", "=", 
   RowBox[{
    RowBox[{"-", "Epot"}], "/", "Ekin"}]}], ";"}]], "Input",
 CellChangeTimes->{
  3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, {
   3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 
   3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, {
   3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 
   3.5788964212998056`*^9}, {3.57889769338844*^9, 3.5788976974288473`*^9}, {
   3.5788979843601513`*^9, 3.578897986466155*^9}, {3.578902015393203*^9, 
   3.578902039620045*^9}, {3.5789024960924473`*^9, 3.578902497106449*^9}, {
   3.579239787767781*^9, 3.579239823413843*^9}, {3.579256422977277*^9, 
   3.5792564375789022`*^9}, 3.579256757613465*^9, {3.5792569515842056`*^9, 
   3.5792569528634076`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"func", " ", "=", 
   RowBox[{"Ekin", "+", "Epot"}]}], ";"}]], "Input",
 CellChangeTimes->{
  3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, {
   3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 
   3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, {
   3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 
   3.57889639532576*^9}, {3.5788978015434303`*^9, 3.5788978862671785`*^9}, {
   3.5788979169212327`*^9, 3.5788979187152357`*^9}, {3.5789023036037087`*^9, 
   3.578902308065317*^9}, {3.578902360762209*^9, 3.57890236119901*^9}, {
   3.578902644932308*^9, 3.578902662295139*^9}, {3.578903281959427*^9, 
   3.578903341645132*^9}, {3.5789033823612037`*^9, 3.5789034163380632`*^9}, {
   3.579238017719472*^9, 3.5792380215726786`*^9}, {3.579238063037551*^9, 
   3.5792381365136805`*^9}, {3.5792383627608776`*^9, 
   3.5792384035705495`*^9}, {3.579238444427021*^9, 3.579238476875078*^9}, {
   3.5792385863872705`*^9, 3.579238636338558*^9}, {3.579238680533436*^9, 
   3.5792387403283405`*^9}, {3.579238925812667*^9, 3.5792389587287245`*^9}, {
   3.5792390123772187`*^9, 3.5792390173224277`*^9}, {3.579239059130501*^9, 
   3.579239067835316*^9}, {3.579239175990306*^9, 3.579239176988708*^9}, {
   3.5792396087730665`*^9, 3.5792396113470707`*^9}, {3.5792396550583477`*^9, 
   3.579239678083988*^9}, {3.579239718129258*^9, 3.5792397203444624`*^9}, {
   3.57923983269586*^9, 3.579239865268717*^9}, {3.579240600888009*^9, 
   3.579240673256536*^9}, {3.579240710556202*^9, 3.579240711445403*^9}, 
   3.579241293248025*^9, {3.5792413542909327`*^9, 3.579241380077778*^9}, {
   3.579241480526354*^9, 3.5792415281532373`*^9}, {3.579256407065249*^9, 
   3.5792564083132515`*^9}, {3.5792568138671637`*^9, 3.579256890338498*^9}, {
   3.57925694283259*^9, 3.5792569442365923`*^9}, {3.5792572545055375`*^9, 
   3.579257255191939*^9}, 3.579257291852003*^9, {3.579260811769387*^9, 
   3.579260834217826*^9}, {3.5792609483476267`*^9, 3.5792609527312346`*^9}, {
   3.5792609956469097`*^9, 3.5792610021521215`*^9}, 3.5792610533046107`*^9, 
   3.579265453093303*^9, 3.579273355392335*^9, {3.579277746144847*^9, 
   3.579277749530053*^9}, {3.5792779342655773`*^9, 3.579277940255988*^9}, 
   3.579278128049118*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"t", " ", "=", " ", 
   RowBox[{"NMinimize", "[", 
    RowBox[{
     RowBox[{"{", "func", "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"R1", ",", " ", "0.25", ",", " ", "0.27"}], "}"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{"R2", ",", " ", "1.19", ",", " ", "1.22"}], "}"}], ",", " ", 
       "\n", "   ", 
       RowBox[{"{", 
        RowBox[{"R3", ",", " ", "1.24", ",", " ", "1.28"}], "}"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{"p", ",", " ", "0.53", ",", "0.56"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"wi", ",", "1.2", ",", "1.25"}], "}"}]}], "}"}]}], "]"}]}], 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{
  3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, {
   3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 
   3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, {
   3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 
   3.57889639532576*^9}, {3.5788978015434303`*^9, 3.5788978862671785`*^9}, {
   3.5788979169212327`*^9, 3.5788979187152357`*^9}, {3.5789023036037087`*^9, 
   3.578902308065317*^9}, {3.578902360762209*^9, 3.57890236119901*^9}, {
   3.578902644932308*^9, 3.578902662295139*^9}, {3.578903281959427*^9, 
   3.578903341645132*^9}, {3.5789033823612037`*^9, 3.5789034163380632`*^9}, {
   3.579238017719472*^9, 3.5792380215726786`*^9}, {3.579238063037551*^9, 
   3.5792381365136805`*^9}, {3.5792383627608776`*^9, 
   3.5792384035705495`*^9}, {3.579238444427021*^9, 3.579238476875078*^9}, {
   3.5792385863872705`*^9, 3.579238636338558*^9}, {3.579238680533436*^9, 
   3.5792387403283405`*^9}, {3.579238925812667*^9, 3.5792389587287245`*^9}, {
   3.5792390123772187`*^9, 3.5792390173224277`*^9}, {3.579239059130501*^9, 
   3.579239067835316*^9}, {3.579239175990306*^9, 3.579239176988708*^9}, {
   3.5792396087730665`*^9, 3.5792396113470707`*^9}, {3.5792396550583477`*^9, 
   3.579239678083988*^9}, {3.579239718129258*^9, 3.5792397203444624`*^9}, {
   3.57923983269586*^9, 3.579239865268717*^9}, {3.579240600888009*^9, 
   3.579240673256536*^9}, {3.579240710556202*^9, 3.579240711445403*^9}, 
   3.579241293248025*^9, {3.5792413542909327`*^9, 3.579241380077778*^9}, {
   3.579241480526354*^9, 3.5792415281532373`*^9}, {3.579256407065249*^9, 
   3.5792564083132515`*^9}, {3.5792568138671637`*^9, 3.579256890338498*^9}, {
   3.57925694283259*^9, 3.5792569442365923`*^9}, {3.5792572545055375`*^9, 
   3.579257255191939*^9}, 3.579257291852003*^9, {3.579260811769387*^9, 
   3.579260834217826*^9}, {3.5792609483476267`*^9, 3.5792609527312346`*^9}, {
   3.5792609956469097`*^9, 3.5792610021521215`*^9}, 3.5792610533046107`*^9, 
   3.579265453093303*^9, 3.579273355392335*^9, {3.579277746144847*^9, 
   3.579277749530053*^9}, {3.5792779342655773`*^9, 3.579277940255988*^9}, 
   3.579278128049118*^9, 3.579279475142684*^9, {3.579279593234891*^9, 
   3.5792795983985004`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "79.4941838100811`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"R1", "\[Rule]", "0.2622399840848178`"}], ",", 
     RowBox[{"R2", "\[Rule]", "1.2000052975561055`"}], ",", 
     RowBox[{"R3", "\[Rule]", "1.261980036069545`"}], ",", 
     RowBox[{"p", "\[Rule]", "0.5440916298444888`"}], ",", 
     RowBox[{"wi", "\[Rule]", "1.2270391929835192`"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.5793183918037*^9, 3.579318436872179*^9, 
  3.5793185152623167`*^9, 3.579318922735032*^9, 3.5793189866951447`*^9, 
  3.5793191723666706`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"u", " ", "=", " ", 
   RowBox[{"t", "[", 
    RowBox[{"[", "2", "]"}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{
  3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, {
   3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 
   3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, {
   3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 
   3.57889639532576*^9}, 3.578896439614238*^9, {3.578902116762181*^9, 
   3.5789021273389993`*^9}, 3.5789025446709323`*^9, 3.5789025890062103`*^9, {
   3.5789028686835012`*^9, 3.578902880882723*^9}, {3.5789029284316063`*^9, 
   3.578902967977676*^9}, {3.5789030373665977`*^9, 3.5789031226987476`*^9}, {
   3.578903474307765*^9, 3.5789034887065907`*^9}, {3.578930301359723*^9, 
   3.5789303077401342`*^9}, {3.578930371263446*^9, 3.5789303949442873`*^9}, 
   3.578930448764382*^9, 3.579101184396015*^9, {3.5792316082810135`*^9, 
   3.57923163459826*^9}, {3.5792316901987576`*^9, 3.579231853687045*^9}, 
   3.579231923980768*^9, {3.579232395803597*^9, 3.5792324440544815`*^9}, 
   3.579273624695608*^9, {3.5792777675480843`*^9, 3.5792777701220894`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Total", " ", "Energy"}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"Etot", "=", 
   RowBox[{"t", "[", 
    RowBox[{"[", "1", "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{
  3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, {
   3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 
   3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, {
   3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 
   3.57889639532576*^9}, 3.578896439614238*^9, {3.578902116762181*^9, 
   3.5789021273389993`*^9}, 3.5789025446709323`*^9, 3.5789025890062103`*^9, {
   3.5789028686835012`*^9, 3.578902880882723*^9}, {3.5789029284316063`*^9, 
   3.578902967977676*^9}, {3.5789030373665977`*^9, 3.5789031226987476`*^9}, {
   3.578903474307765*^9, 3.5789034887065907`*^9}, {3.578930301359723*^9, 
   3.5789303077401342`*^9}, {3.578930371263446*^9, 3.5789303949442873`*^9}, 
   3.578930448764382*^9, 3.579101184396015*^9, {3.5792316082810135`*^9, 
   3.57923163459826*^9}, {3.5792316901987576`*^9, 3.579231853687045*^9}, 
   3.579231923980768*^9, {3.579232395803597*^9, 3.5792324440544815`*^9}, {
   3.579233677423848*^9, 3.5792336815734553`*^9}, 3.579235541368722*^9, {
   3.5792734364032774`*^9, 3.5792734535945077`*^9}, {3.579277779684906*^9, 
   3.5792777872665195`*^9}, 3.579278239121313*^9}],

Cell[BoxData[
 RowBox[{"-", "79.4941838100811`"}]], "Output",
 CellChangeTimes->{3.5793183918037*^9, 3.579318436887779*^9, 
  3.5793185152623167`*^9, 3.579318922735032*^9, 3.579318986710745*^9, 
  3.579319172382271*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Kinetic", " ", "energy"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"Ekin", " ", "/.", "u"}]}]], "Input",
 CellChangeTimes->{
  3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, {
   3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 
   3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, {
   3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 
   3.57889639532576*^9}, 3.578896439614238*^9, {3.578902116762181*^9, 
   3.5789021273389993`*^9}, 3.5789025446709323`*^9, 3.5789025890062103`*^9, {
   3.5789028686835012`*^9, 3.578902880882723*^9}, {3.5789029284316063`*^9, 
   3.578902967977676*^9}, {3.5789030373665977`*^9, 3.5789031226987476`*^9}, {
   3.578903474307765*^9, 3.5789034887065907`*^9}, {3.578930301359723*^9, 
   3.5789303077401342`*^9}, {3.578930371263446*^9, 3.5789303949442873`*^9}, 
   3.578930448764382*^9, 3.579101184396015*^9, {3.5792316082810135`*^9, 
   3.57923163459826*^9}, {3.5792316901987576`*^9, 3.579231853687045*^9}, 
   3.579231923980768*^9, 3.5792319734952555`*^9, {3.579233685333062*^9, 
   3.579233688671468*^9}, 3.5792355483575344`*^9, {3.579256527825061*^9, 
   3.5792565289950633`*^9}, 3.5792735029841948`*^9, {3.579277801290944*^9, 
   3.579277804130149*^9}}],

Cell[BoxData["79.49419683114962`"], "Output",
 CellChangeTimes->{3.5793183918192997`*^9, 3.579318436887779*^9, 
  3.579318515277917*^9, 3.5793189227506323`*^9, 3.579318986726345*^9, 
  3.579319172382271*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{" ", 
  RowBox[{"(*", 
   RowBox[{"NUCLEUS", "-", 
    RowBox[{"ELECTRON", " ", "POTENTIAL", " ", "ENERGY"}]}], "  ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"vne", " ", "/.", "u"}]}]], "Input",
 CellChangeTimes->{
  3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, {
   3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 
   3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, {
   3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 
   3.57889639532576*^9}, 3.578896439614238*^9, {3.578902116762181*^9, 
   3.5789021273389993`*^9}, 3.5789025446709323`*^9, 3.5789025890062103`*^9, {
   3.5789028686835012`*^9, 3.578902880882723*^9}, {3.5789029284316063`*^9, 
   3.578902967977676*^9}, {3.5789030373665977`*^9, 3.5789031226987476`*^9}, {
   3.578903474307765*^9, 3.5789034887065907`*^9}, {3.578930301359723*^9, 
   3.5789303077401342`*^9}, {3.578930371263446*^9, 3.5789303949442873`*^9}, 
   3.578930448764382*^9, 3.579101184396015*^9, {3.5792316082810135`*^9, 
   3.57923163459826*^9}, {3.5792316901987576`*^9, 3.579231853687045*^9}, {
   3.579231923980768*^9, 3.579231930080379*^9}, 3.5792319760224595`*^9, {
   3.5792336918850737`*^9, 3.5792336945838785`*^9}, 3.5792355580451517`*^9, {
   3.579256537949479*^9, 3.5792565381210794`*^9}, 3.5792735029997945`*^9, {
   3.5792778082173557`*^9, 3.5792778109629607`*^9}}],

Cell[BoxData[
 RowBox[{"-", "270.426103019044`"}]], "Output",
 CellChangeTimes->{3.5793183918192997`*^9, 3.5793184369033794`*^9, 
  3.579318515277917*^9, 3.5793189227662325`*^9, 3.579318986726345*^9, 
  3.579319172397871*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"ELECTRON", "-", 
    RowBox[{"ELECTRON", " ", "POTENTIAL", " ", "ENERGY"}]}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"vee", " ", "/.", "u"}]}]], "Input",
 CellChangeTimes->{
  3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, {
   3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 
   3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, {
   3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 
   3.57889639532576*^9}, 3.578896439614238*^9, {3.578902116762181*^9, 
   3.5789021273389993`*^9}, 3.5789025446709323`*^9, 3.5789025890062103`*^9, {
   3.5789028686835012`*^9, 3.578902880882723*^9}, {3.5789029284316063`*^9, 
   3.578902967977676*^9}, {3.5789030373665977`*^9, 3.5789031226987476`*^9}, {
   3.578903474307765*^9, 3.5789034887065907`*^9}, {3.578930301359723*^9, 
   3.5789303077401342`*^9}, {3.578930371263446*^9, 3.5789303949442873`*^9}, 
   3.578930448764382*^9, 3.579101184396015*^9, {3.5792316082810135`*^9, 
   3.57923163459826*^9}, {3.5792316901987576`*^9, 3.579231853687045*^9}, {
   3.579231923980768*^9, 3.5792319350255876`*^9}, 3.5792319781128635`*^9, {
   3.579233697766284*^9, 3.5792337003714886`*^9}, 3.579235564737563*^9, {
   3.5792565420990863`*^9, 3.579256542270686*^9}, 3.5792566391000566`*^9, 
   3.5792735029997945`*^9, {3.5792778131625648`*^9, 3.5792778153309684`*^9}}],

Cell[BoxData["69.40599530071324`"], "Output",
 CellChangeTimes->{3.5793183918192997`*^9, 3.5793184369033794`*^9, 
  3.579318515293517*^9, 3.5793189227662325`*^9, 3.579318986741945*^9, 
  3.579319172397871*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"NUCLEUS", "-", 
    RowBox[{"NUCLEUS", " ", "POTENTIAL", " ", "ENERGY"}]}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"vnn", " ", "/.", "u"}]}]], "Input",
 CellChangeTimes->{
  3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, {
   3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 
   3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, {
   3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 
   3.57889639532576*^9}, 3.578896439614238*^9, {3.578902116762181*^9, 
   3.5789021273389993`*^9}, 3.5789025446709323`*^9, 3.5789025890062103`*^9, {
   3.5789028686835012`*^9, 3.578902880882723*^9}, {3.5789029284316063`*^9, 
   3.578902967977676*^9}, {3.5789030373665977`*^9, 3.5789031226987476`*^9}, {
   3.578903474307765*^9, 3.5789034887065907`*^9}, {3.578930301359723*^9, 
   3.5789303077401342`*^9}, {3.578930371263446*^9, 3.5789303949442873`*^9}, 
   3.578930448764382*^9, 3.579101184396015*^9, {3.5792316082810135`*^9, 
   3.57923163459826*^9}, {3.5792316901987576`*^9, 3.579231853687045*^9}, {
   3.579231923980768*^9, 3.579231939689996*^9}, 3.5792319799224663`*^9, {
   3.579233703163893*^9, 3.5792337068142996`*^9}, 3.5792355762971835`*^9, {
   3.579256545499892*^9, 3.5792565456870923`*^9}, 3.5792735029997945`*^9, {
   3.579277817780173*^9, 3.579277822959382*^9}}],

Cell[BoxData["42.031727077100044`"], "Output",
 CellChangeTimes->{3.5793183918349*^9, 3.579318436918979*^9, 
  3.579318515293517*^9, 3.5793189227818327`*^9, 3.5793189867575445`*^9, 
  3.579319172413471*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{" ", 
  RowBox[{"(*", 
   RowBox[{"C", "-", 
    RowBox[{"H", " ", "distance", "  ", "\[CapitalARing]"}]}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"ch", " ", "=", 
     RowBox[{
      RowBox[{"R1", "+", "R3", "+", "p"}], "/.", "u"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"dch", " ", "=", " ", 
    RowBox[{
    "ch", "*", "bohr", 
     "                                             "}]}]}]}]], "Input",
 CellChangeTimes->{
  3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, {
   3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 
   3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, {
   3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 
   3.57889639532576*^9}, 3.578896439614238*^9, {3.578902116762181*^9, 
   3.5789021273389993`*^9}, 3.5789025446709323`*^9, 3.5789025890062103`*^9, {
   3.5789028686835012`*^9, 3.578902880882723*^9}, {3.5789029284316063`*^9, 
   3.578902967977676*^9}, {3.5789030373665977`*^9, 3.5789031226987476`*^9}, {
   3.578903474307765*^9, 3.5789034887065907`*^9}, {3.578930301359723*^9, 
   3.5789303077401342`*^9}, {3.578930371263446*^9, 3.5789303949442873`*^9}, 
   3.578930448764382*^9, 3.579101184396015*^9, {3.5792316082810135`*^9, 
   3.57923163459826*^9}, {3.5792316901987576`*^9, 3.579231853687045*^9}, {
   3.579231923980768*^9, 3.5792319444324045`*^9}, 3.57923198196607*^9, {
   3.579232330751483*^9, 3.579232335509491*^9}, {3.5792337100279055`*^9, 
   3.579233713132311*^9}, 3.579235583473196*^9, {3.5792565803815536`*^9, 
   3.5792566233284287`*^9}, {3.579257470812318*^9, 3.5792574914355536`*^9}, 
   3.579273503015394*^9, {3.5792778318357973`*^9, 3.5792778352678037`*^9}}],

Cell[BoxData["1.0945029540114424`"], "Output",
 CellChangeTimes->{3.5793183918349*^9, 3.579318436918979*^9, 
  3.579318515309117*^9, 3.5793189227818327`*^9, 3.5793189867731447`*^9, 
  3.579319172413471*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{" ", 
  RowBox[{"(*", 
   RowBox[{"C", "-", 
    RowBox[{"C", " ", "distance", "  ", "\[CapitalARing]"}]}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"cc", " ", "=", "  ", 
     RowBox[{
      RowBox[{"2", "*", "d1"}], "  ", "/.", "u"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"dcc", "=", 
    RowBox[{
    "cc", "*", "bohr", 
     "                                              "}]}]}]}]], "Input",
 CellChangeTimes->{
  3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, {
   3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 
   3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, {
   3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 
   3.57889639532576*^9}, 3.578896439614238*^9, {3.578902116762181*^9, 
   3.5789021273389993`*^9}, 3.5789025446709323`*^9, 3.5789025890062103`*^9, {
   3.5789028686835012`*^9, 3.578902880882723*^9}, {3.5789029284316063`*^9, 
   3.578902967977676*^9}, {3.5789030373665977`*^9, 3.5789031226987476`*^9}, {
   3.578903474307765*^9, 3.5789034887065907`*^9}, {3.578930301359723*^9, 
   3.5789303077401342`*^9}, {3.578930371263446*^9, 3.5789303949442873`*^9}, 
   3.578930448764382*^9, 3.579101184396015*^9, {3.5792316082810135`*^9, 
   3.57923163459826*^9}, {3.5792316901987576`*^9, 3.579231853687045*^9}, {
   3.579231923980768*^9, 3.5792319537300205`*^9}, 3.579231983994074*^9, 
   3.5792322853398027`*^9, 3.579232323949871*^9, {3.579233721228725*^9, 
   3.57923372403673*^9}, 3.579235589853607*^9, {3.5792566088672037`*^9, 
   3.5792566162616167`*^9}, 3.579273503015394*^9, {3.5792778392614107`*^9, 
   3.579277842319016*^9}}],

Cell[BoxData["1.5475731428057977`"], "Output",
 CellChangeTimes->{3.5793183918505*^9, 3.579318436934579*^9, 
  3.579318515309117*^9, 3.579318922797433*^9, 3.5793189867731447`*^9, 
  3.579319172429071*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"angle", " ", "CCH"}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{"\[Pi]", "-", "wi"}], ")"}], "*", "rad"}], " ", "/.", " ", 
   "u"}]}]], "Input",
 CellChangeTimes->{{3.5792797775492153`*^9, 3.579279833818514*^9}, {
  3.579279914813856*^9, 3.579279930460684*^9}}],

Cell[BoxData["109.69583291985901`"], "Output",
 CellChangeTimes->{3.5793183918505*^9, 3.579318436934579*^9, 
  3.5793185153247166`*^9, 3.5793189228130326`*^9, 3.579318986788745*^9, 
  3.579319172444671*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Virial", " ", "Theorem"}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"vir", "/.", "u", "      "}]}]], "Input",
 CellChangeTimes->{
  3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, {
   3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 
   3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, {
   3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 
   3.57889639532576*^9}, 3.578896439614238*^9, {3.578902116762181*^9, 
   3.5789021273389993`*^9}, 3.5789025446709323`*^9, 3.5789025890062103`*^9, {
   3.5789028686835012`*^9, 3.578902880882723*^9}, {3.5789029284316063`*^9, 
   3.578902967977676*^9}, {3.5789030373665977`*^9, 3.5789031226987476`*^9}, {
   3.578903474307765*^9, 3.5789034887065907`*^9}, {3.578930301359723*^9, 
   3.5789303077401342`*^9}, {3.578930371263446*^9, 3.5789303949442873`*^9}, 
   3.578930448764382*^9, 3.579101184396015*^9, {3.5792316082810135`*^9, 
   3.57923163459826*^9}, {3.5792316901987576`*^9, 3.579231853687045*^9}, {
   3.579231923980768*^9, 3.579231986490078*^9}, {3.579233726953935*^9, 
   3.5792337305575414`*^9}, 3.579235598917223*^9, {3.579240004389761*^9, 
   3.579240005544163*^9}, {3.5792566498484755`*^9, 3.579256654731284*^9}, 
   3.579273503015394*^9, {3.579277847763426*^9, 3.57927785046223*^9}}],

Cell[BoxData["1.9999998362010167`"], "Output",
 CellChangeTimes->{3.5793183918505*^9, 3.579318436950179*^9, 
  3.5793185153247166`*^9, 3.5793189228286324`*^9, 3.579318986788745*^9, 
  3.579319172444671*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "Atomization", " ", "Energy", " ", "at", " ", "0", " ", "K", " ", "in", 
    " ", 
    RowBox[{"kcal", "/", "mol"}]}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "Etot"}], " ", "-", " ", 
      RowBox[{"nn", "*", "37.78899"}], "  ", "-", 
      RowBox[{
       RowBox[{"(", " ", 
        RowBox[{
         RowBox[{"2", "*", "nn"}], " ", "+", " ", "2"}], ")"}], "*", 
       "0.50177"}]}], ")"}], "*", "627.5096"}], "  ", "/.", "u"}]}]], "Input",
 CellChangeTimes->{
  3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, {
   3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 
   3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, {
   3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 
   3.57889639532576*^9}, 3.578896439614238*^9, {3.578902116762181*^9, 
   3.5789021273389993`*^9}, 3.5789025446709323`*^9, 3.5789025890062103`*^9, {
   3.5789028686835012`*^9, 3.578902880882723*^9}, {3.5789029284316063`*^9, 
   3.578902967977676*^9}, {3.5789030373665977`*^9, 3.5789031226987476`*^9}, {
   3.578903474307765*^9, 3.5789034887065907`*^9}, {3.578930301359723*^9, 
   3.5789303077401342`*^9}, {3.578930371263446*^9, 3.5789303949442873`*^9}, 
   3.578930448764382*^9, 3.579101184396015*^9, {3.5792316082810135`*^9, 
   3.57923163459826*^9}, {3.5792316901987576`*^9, 3.579231853687045*^9}, {
   3.579231923980768*^9, 3.5792319607344327`*^9}, {3.5792337327883453`*^9, 
   3.579233738778756*^9}, {3.579234893452784*^9, 3.579234901361998*^9}, 
   3.5792356093536415`*^9, {3.5792569813646574`*^9, 3.579256983189861*^9}, 
   3.5792575835225153`*^9, {3.579259029379855*^9, 3.5792590354014654`*^9}, 
   3.579259348790416*^9, {3.5792789528537664`*^9, 3.579278982649819*^9}, {
   3.579279122972065*^9, 3.5792791254368696`*^9}, {3.5793189037497993`*^9, 
   3.5793189165262213`*^9}, {3.5793189559318905`*^9, 3.5793189714695177`*^9}, 
   3.579319165705459*^9}],

Cell[BoxData["568.262534430467`"], "Output",
 CellChangeTimes->{3.579319172460271*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"plot0", "=", 
   RowBox[{
    RowBox[{"Graphics3D", "[", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Opacity", "[", "0.6", "]"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Darker", "[", 
          RowBox[{"Orange", ",", "0.4"}], "]"}], ",", 
         RowBox[{"Sphere", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"w", "[", 
             RowBox[{"[", "1", "]"}], "]"}], "}"}], ",", 
           RowBox[{"c", "[", 
            RowBox[{"[", "1", "]"}], "]"}]}], "]"}], ",", 
         RowBox[{"Sphere", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"w", "[", 
             RowBox[{"[", "2", "]"}], "]"}], "}"}], ",", 
           RowBox[{"c", "[", 
            RowBox[{"[", "2", "]"}], "]"}]}], "]"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Lighter", "[", 
          RowBox[{"Orange", ",", "0.6"}], "]"}], ",", 
         RowBox[{"Sphere", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"w", "[", 
             RowBox[{"[", "3", "]"}], "]"}], "}"}], ",", 
           RowBox[{"c", "[", 
            RowBox[{"[", "3", "]"}], "]"}]}], "]"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Darker", "[", 
          RowBox[{"Red", ",", "1.0"}], "]"}], ",", 
         RowBox[{"Sphere", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"n", "[", 
             RowBox[{"[", "1", "]"}], "]"}], "}"}], ",", "0.08"}], "]"}], ",", 
         RowBox[{"Sphere", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"n", "[", 
             RowBox[{"[", "2", "]"}], "]"}], "}"}], ",", "0.08"}], "]"}]}], 
        "}"}]}], "}"}], "]"}], " ", "/.", "u"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"plot1", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics3D", "[", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Opacity", "[", "0.6", "]"}], ",", 
         RowBox[{"Sphere", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"w", "[", 
             RowBox[{"[", "i", "]"}], "]"}], "}"}], ",", 
           RowBox[{"c", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"Darker", "[", 
            RowBox[{"Red", ",", "1.0"}], "]"}], ",", 
           RowBox[{"Sphere", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{"n", "[", 
               RowBox[{"[", "i", "]"}], "]"}], "}"}], ",", "0.08"}], "]"}]}], 
          "}"}]}], "}"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "4", ",", "nc"}], "}"}]}], "]"}], "   ", "/.", " ", 
    "u"}]}], "  ", ";"}], "\n", 
 RowBox[{"Show", "[", 
  RowBox[{"plot0", ",", "plot1", ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", 
     RowBox[{"Axes", " ", "->", " ", "True"}], ",", 
     RowBox[{"SphericalRegion", "\[Rule]", "True"}]}], "}"}]}], 
  "]"}]}], "Input",
 CellChangeTimes->{
  3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, {
   3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 
   3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, {
   3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 
   3.57889639532576*^9}, 3.578896439614238*^9, {3.578902116762181*^9, 
   3.5789021273389993`*^9}, 3.5789025446709323`*^9, 3.5789025890062103`*^9, {
   3.5789028686835012`*^9, 3.578902880882723*^9}, {3.5789029284316063`*^9, 
   3.578902967977676*^9}, {3.5789030373665977`*^9, 3.5789031226987476`*^9}, {
   3.578903474307765*^9, 3.5789034887065907`*^9}, {3.578930301359723*^9, 
   3.5789303077401342`*^9}, {3.578930371263446*^9, 3.5789303949442873`*^9}, 
   3.578930448764382*^9, 3.579101184396015*^9, {3.5792316082810135`*^9, 
   3.57923163459826*^9}, {3.5792316901987576`*^9, 3.579231853687045*^9}, {
   3.579231923980768*^9, 3.5792319607344327`*^9}, {3.5792337327883453`*^9, 
   3.579233738778756*^9}, {3.579234893452784*^9, 3.579234901361998*^9}, 
   3.5792356093536415`*^9, {3.5792569813646574`*^9, 3.579256983189861*^9}, 
   3.5792575835225153`*^9, 3.579257631055799*^9, 3.5792736529160576`*^9, {
   3.579277858823845*^9, 3.5792778625678515`*^9}, {3.579278807960712*^9, 
   3.5792788207215347`*^9}, {3.5793171650019455`*^9, 3.5793171770919666`*^9}, 
   3.579317208448022*^9, {3.5793172568705063`*^9, 3.579317269631329*^9}, {
   3.579317305105791*^9, 3.579317313077405*^9}, {3.579317359003886*^9, 
   3.5793175029765387`*^9}, {3.579317557857435*^9, 3.579317599291108*^9}, {
   3.5793176332991676`*^9, 3.579317667198027*^9}, {3.5793177909218445`*^9, 
   3.5793178397499304`*^9}, {3.5793180136434355`*^9, 3.579318057136312*^9}, {
   3.579318124949631*^9, 3.5793181894869447`*^9}, {3.5793182499838505`*^9, 
   3.579318295489131*^9}, {3.579318365767254*^9, 3.5793183737544684`*^9}, {
   3.579318429025365*^9, 3.579318430616568*^9}, {3.579318492673477*^9, 
   3.579318505325099*^9}}],

Cell[BoxData[
 Graphics3DBox[{
   {Opacity[0.6], 
    {RGBColor[0.6, 0.3, 0.], 
     SphereBox[{{1.4622452816409233`, 0, 0}}, 0.2622399840848178], 
     SphereBox[{{-1.4622452816409233`, 0, 0}}, 0.2622399840848178]}, 
    {RGBColor[1., 0.8, 0.6], SphereBox[{{0, 0, 0}}, 1.2000052975561055`]}, 
    {RGBColor[0., 0., 0.], SphereBox[{{1.4622452816409233`, 0, 0}}, 0.08], 
     SphereBox[{{-1.4622452816409233`, 0, 0}}, 0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{{1.975948255396272, 1.4350456176004571`, 0}}, 
     1.261980036069545], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{{2.1593217084087852`, 1.9473051986695669`, 0}}, 0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{{1.975948255396272, -0.7175228088002286, 1.242785960431525}}, 
     1.261980036069545], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{{2.1593217084087852`, -0.9736525993347834, 
      1.6864157709693481`}}, 0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{{1.975948255396272, -0.7175228088002286, -1.242785960431525}}, 
     1.261980036069545], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{{
      2.1593217084087852`, -0.9736525993347834, -1.6864157709693481`}}, 
      0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{{-1.975948255396272, -1.4350456176004571`, 0}}, 
     1.261980036069545], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{{-2.1593217084087852`, -1.9473051986695669`, 0}}, 0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{{-1.975948255396272, 0.7175228088002286, -1.242785960431525}}, 
     1.261980036069545], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{{-2.1593217084087852`, 
      0.9736525993347834, -1.6864157709693481`}}, 0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{{-1.975948255396272, 0.7175228088002286, 1.242785960431525}}, 
     1.261980036069545], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{{-2.1593217084087852`, 0.9736525993347834, 
      1.6864157709693481`}}, 0.08]}}},
  AspectRatio->Automatic,
  Axes->True,
  ImageSize->{407.0116941416884, 407.},
  SphericalRegion->True,
  ViewPoint->{1.5866929698981913`, -2.636731848607672, 1.4071426998757615`},
  ViewVertical->{0.13189268980890168`, -0.18824632091226143`, 
   1.2652809336016606`}]], "Output",
 CellChangeTimes->{3.5793183918661003`*^9, 3.5793184369657793`*^9, 
  3.579318515355917*^9, 3.5793189228442326`*^9, 3.5793189868199453`*^9, 
  3.579319172475871*^9}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{990, 810},
WindowMargins->{{452, Automatic}, {Automatic, 25}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}},
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 901, 15, 33, "Subsubtitle"],
Cell[1483, 39, 4208, 70, 152, "Input"],
Cell[5694, 111, 12629, 325, 692, "Text"],
Cell[18326, 438, 843, 15, 31, "Input"],
Cell[19172, 455, 2308, 34, 31, "Input"],
Cell[CellGroupData[{
Cell[21505, 493, 2988, 51, 72, "Input"],
Cell[24496, 546, 605, 14, 31, "Output"]
}, Open  ]],
Cell[25116, 563, 1202, 19, 31, "Input"],
Cell[CellGroupData[{
Cell[26343, 586, 1407, 23, 52, "Input"],
Cell[27753, 611, 219, 4, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[28009, 620, 1330, 21, 52, "Input"],
Cell[29342, 643, 207, 3, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[29586, 651, 1425, 23, 52, "Input"],
Cell[31014, 676, 225, 4, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[31276, 685, 1442, 23, 52, "Input"],
Cell[32721, 710, 209, 3, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[32967, 718, 1412, 23, 52, "Input"],
Cell[34382, 743, 206, 3, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[34625, 751, 1763, 33, 72, "Input"],
Cell[36391, 786, 206, 3, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[36634, 794, 1708, 33, 72, "Input"],
Cell[38345, 829, 204, 3, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[38586, 837, 369, 10, 52, "Input"],
Cell[38958, 849, 206, 3, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[39201, 857, 1382, 21, 52, "Input"],
Cell[40586, 880, 206, 3, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[40829, 888, 2048, 38, 52, "Input"],
Cell[42880, 928, 86, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[43003, 934, 5062, 117, 152, "Input"],
Cell[48068, 1053, 2307, 53, 454, "Output"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)