(* 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[ 18064, 574] NotebookOptionsPosition[ 16340, 514] NotebookOutlinePosition[ 16853, 533] CellTagsIndexPosition[ 16810, 530] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{"LiH", " ", "molecule", " ", "with", " ", RowBox[{"Li", "[", RowBox[{"1", "s2", " ", "2", "s"}], "]"}]}], ",", " ", "overlap", ",", " ", RowBox[{"parametrized", " ", "ES", " ", "31.12", ".2012"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"Variables", ":", " ", "P"}], ",", "Q", ",", "R", ",", "S"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"with", " ", StyleBox["exchange", FontColor->RGBColor[1, 0, 0]]}], " ", "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"Clear", "[", RowBox[{ "k1", ",", "k2", ",", "s1", ",", "s2", ",", " ", "Z", ",", " ", "P", ",", "Q", ",", "R", ",", "S"}], "]"}], ";"}]}]], "Input", CellChangeTimes->{{3.5659569125308685`*^9, 3.565956917678878*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"Z", "=", "3.0"}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"T", "=", RowBox[{ RowBox[{"2.25", "*", RowBox[{"k1", "/", RowBox[{"P", "^", "2"}]}]}], "+", RowBox[{"9.0", "*", RowBox[{"k2", "/", RowBox[{"Q", "^", "2"}]}]}]}]}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"Vne", "=", RowBox[{ RowBox[{ RowBox[{"-", "3.0"}], "*", RowBox[{"Z", "/", "P"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"2.0", "*", RowBox[{"Z", "/", "Q"}]}], ")"}], "*", RowBox[{"(", RowBox[{"1.5", "-", RowBox[{"0.5", "*", RowBox[{ RowBox[{"(", RowBox[{"R", "/", "Q"}], ")"}], "^", "2"}]}]}], ")"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"2", "/", "Q"}], ")"}], "*", RowBox[{"(", RowBox[{"1.5", "-", RowBox[{"0.5", "*", RowBox[{ RowBox[{"(", RowBox[{"S", "/", "Q"}], ")"}], "^", "2"}]}]}], ")"}]}], "-", RowBox[{"2.0", "/", RowBox[{"(", RowBox[{"R", "+", "S"}], ")"}]}]}]}], ";"}]], "Input"], Cell[BoxData[{ RowBox[{ RowBox[{"Vee", "=", RowBox[{ RowBox[{"3.0", "*", RowBox[{"s1", "/", "P"}]}], "+", RowBox[{"3.0", "*", RowBox[{"s2", "/", "Q"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "Q"}], ")"}], "*", RowBox[{"(", RowBox[{"1.5", "-", RowBox[{"0.5", "*", RowBox[{ RowBox[{"(", RowBox[{"R", "/", "Q"}], ")"}], "^", "2"}]}], "-", RowBox[{"0.3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", "/", "Q"}], ")"}], "^", "2"}]}]}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vex", " ", "=", " ", RowBox[{ RowBox[{"-", StyleBox["1.2", FontColor->RGBColor[1, 0, 0]]}], StyleBox["*", FontColor->RGBColor[1, 0, 0]], StyleBox[ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"P", "/", "Q"}], ")"}], "^", "3"}], "/", "Q"}], FontColor->RGBColor[1, 0, 0]]}]}], StyleBox[";", FontColor->RGBColor[1, 0, 0]]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vee", " ", "=", " ", RowBox[{"Vee", "+", "Vex"}]}], ";"}]}], "Input", CellChangeTimes->{{3.5659557885800943`*^9, 3.5659558088601303`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"Vnn", "=", RowBox[{"Z", "/", RowBox[{"(", RowBox[{"R", "+", "S"}], ")"}]}]}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ StyleBox["(*", FontSize->12], StyleBox[" ", FontSize->12], RowBox[{"Parameters", " ", "in", " ", "list", " ", "c"}], StyleBox[" ", FontSize->12], StyleBox["*)", FontSize->12]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{"c", "=", RowBox[{"{", RowBox[{ RowBox[{"k1", " ", "->", " ", "1.0"}], ",", " ", RowBox[{"k2", "->", " ", "1.0"}], ",", " ", RowBox[{"s1", " ", "->", " ", "0.3"}], ",", " ", RowBox[{"s2", " ", "->", " ", "0.3"}]}], "}"}]}], ";"}], " ", "*)"}], " ", "\[IndentingNewLine]", RowBox[{ RowBox[{"c", "=", RowBox[{"{", RowBox[{ RowBox[{"k1", " ", "\[Rule]", " ", "1.01613950"}], ",", " ", RowBox[{"k2", " ", "->", " ", "0.65858068"}], ",", RowBox[{"s1", "->", " ", "0.25742249"}], ",", RowBox[{"s2", "->", " ", "0.22958899"}]}], "}"}]}], ";"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"func", "=", RowBox[{ RowBox[{"T", "+", "Vne", "+", "Vee", "+", "Vnn"}], " ", "/.", " ", "c"}]}], ";"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"t", "=", RowBox[{"FindMinimum", "[", RowBox[{"func", ",", RowBox[{"{", RowBox[{"P", ",", "0.557"}], "}"}], ",", RowBox[{"{", RowBox[{"Q", ",", "2.99"}], "}"}], ",", RowBox[{"{", RowBox[{"R", ",", "1.49"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"S", ",", " ", "1.49"}], "}"}]}], "]"}], " "}]], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "8.03041514112117`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.5571059532347319`"}], ",", RowBox[{"Q", "\[Rule]", "2.9878874996438496`"}], ",", RowBox[{"R", "\[Rule]", "1.4939437503693909`"}], ",", RowBox[{"S", "\[Rule]", "1.4939437503693915`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.5659559086691055`*^9, 3.5659559824572353`*^9, 3.5659561686435623`*^9, 3.5659562232592583`*^9, 3.565956264755331*^9, 3.5659563493074794`*^9, 3.565956388775549*^9, {3.565956422471608*^9, 3.565956450879258*^9}, { 3.5659565650714583`*^9, 3.5659565811082864`*^9}, 3.5659567740962257`*^9, { 3.5659568051090803`*^9, 3.5659568167623005`*^9}}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"u", "=", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}], ";"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Vne", " ", "/.", "u"}]], "Input"], Cell[BoxData[ RowBox[{"-", "20.50581631217891`"}]], "Output", CellChangeTimes->{ 3.5659559086847057`*^9, 3.565955982472835*^9, 3.565956168659162*^9, 3.565956223274858*^9, 3.565956264755331*^9, 3.565956349323079*^9, 3.565956388791149*^9, {3.565956422471608*^9, 3.565956450894858*^9}, { 3.5659565650714583`*^9, 3.565956581123886*^9}, 3.565956774111826*^9, { 3.56595680512468*^9, 3.5659568167623005`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"Vee", " ", "/.", "u"}], " ", "/.", " ", "c"}]], "Input"], Cell[BoxData["3.4409321461347626`"], "Output", CellChangeTimes->{ 3.565955908700306*^9, 3.5659559824884353`*^9, 3.565956168674762*^9, 3.5659562232904577`*^9, 3.565956264770931*^9, 3.565956349323079*^9, 3.565956388806749*^9, {3.565956422487208*^9, 3.565956450894858*^9}, { 3.565956565087058*^9, 3.5659565811394863`*^9}, 3.5659567741274257`*^9, { 3.56595680512468*^9, 3.5659568167779007`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Vex", " ", "/.", " ", "u"}]], "Input"], Cell[BoxData[ RowBox[{"-", "0.0026033806871996043`"}]], "Output", CellChangeTimes->{ 3.565955908700306*^9, 3.565955982504035*^9, 3.565956168674762*^9, 3.5659562232904577`*^9, 3.5659562647865305`*^9, 3.5659563493386793`*^9, 3.565956388806749*^9, {3.565956422487208*^9, 3.5659564509104576`*^9}, { 3.565956565087058*^9, 3.5659565811394863`*^9}, 3.5659567741274257`*^9, { 3.56595680514028*^9, 3.5659568167779007`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Vnn", " ", "/.", " ", "u"}]], "Input"], Cell[BoxData["1.0040538672417294`"], "Output", CellChangeTimes->{ 3.5659559087159057`*^9, 3.565955982519635*^9, 3.565956168690362*^9, 3.565956223306058*^9, 3.5659562647865305`*^9, 3.5659563493542795`*^9, 3.5659563888223486`*^9, {3.565956422502808*^9, 3.5659564509104576`*^9}, { 3.5659565651026583`*^9, 3.5659565811550865`*^9}, 3.5659567741430254`*^9, { 3.56595680514028*^9, 3.5659568167935004`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"Vne", "+", "Vee", "+", "Vnn"}], ")"}]}], "/", "T"}], " ", "/.", "u"}], "/.", "c"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"R", "+", "S"}], ")"}], "*", "0.529177"}], " ", "/.", "u"}]}], "Input"], Cell[BoxData["1.999999997937829`"], "Output", CellChangeTimes->{ 3.5659559087315054`*^9, 3.565955982519635*^9, 3.565956168690362*^9, 3.565956223306058*^9, 3.5659562648021307`*^9, 3.5659563493542795`*^9, 3.5659563888223486`*^9, {3.565956422502808*^9, 3.565956450926058*^9}, { 3.5659565651026583`*^9, 3.5659565811550865`*^9}, 3.5659567741430254`*^9, { 3.56595680515588*^9, 3.5659568168091*^9}}], Cell[BoxData["1.5811213439784468`"], "Output", CellChangeTimes->{ 3.5659559087315054`*^9, 3.565955982519635*^9, 3.565956168690362*^9, 3.565956223306058*^9, 3.5659562648021307`*^9, 3.5659563493542795`*^9, 3.5659563888223486`*^9, {3.565956422502808*^9, 3.565956450926058*^9}, { 3.5659565651026583`*^9, 3.5659565811550865`*^9}, 3.5659567741430254`*^9, { 3.56595680515588*^9, 3.5659568168091*^9}}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"plot1", "=", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "R"}], ",", "0"}], "}"}], ",", "P"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "R"}], ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", "\n", "\t\t\t", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "Q"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{"S", ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", "\n", "\t\t\t", RowBox[{"{", RowBox[{ RowBox[{"Thickness", "[", "0.012", "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "R"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"S", ",", "0"}], "}"}]}], "}"}], "]"}]}], "}"}]}], "}"}], "]"}], "/.", " ", "u"}]}], ";"}], " "}]], "Input", CellChangeTimes->{ 3.565955928761941*^9, {3.565956646144801*^9, 3.5659566483444047`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ StyleBox[ RowBox[{"Show", "[", RowBox[{"plot1", ",", " ", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", "->", "Automatic"}], ",", " ", RowBox[{"PlotRange", "->", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "3.5"}], ",", "3.5"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "3.5"}], ",", "3.5"}], "}"}]}], "}"}]}], ",", RowBox[{"Axes", "->", "True"}], ",", RowBox[{"GridLines", "->", "Automatic"}], ",", RowBox[{"Frame", "->", "True"}]}], "}"}]}], "]"}], FontSize->12]], "Input"], Cell[BoxData[ GraphicsBox[{CircleBox[{-1.4939437503693909`, 0}, 0.5571059532347319], DiskBox[{-1.4939437503693909`, 0}, 0.1], CircleBox[{0, 0}, 2.9878874996438496`], DiskBox[{1.4939437503693915`, 0}, 0.1], {Thickness[0.012], LineBox[{{-1.4939437503693909`, 0}, {1.4939437503693915`, 0}}]}}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, PlotRange->{{-3.5, 3.5}, {-3.5, 3.5}}]], "Output", CellChangeTimes->{ 3.565955908762706*^9, 3.565955982550835*^9, 3.565956168705962*^9, 3.565956223321658*^9, 3.5659562648021307`*^9, 3.5659563493698797`*^9, 3.565956388837949*^9, {3.565956422518408*^9, 3.565956450941658*^9}, { 3.5659565651182585`*^9, 3.5659565811706867`*^9}, 3.5659567741586256`*^9, { 3.56595680517148*^9, 3.5659568168091*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"plot2", "=", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.6", "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "R"}], ",", "0", ",", "0"}], "}"}], ",", "P"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Darker", "[", RowBox[{"Red", ",", "1"}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "R"}], ",", "0", ",", "0"}], "}"}], ",", "0.1"}], "]"}]}], "}"}], ",", "\n", "\t ", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.4", "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", "Q"}], "]"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"Darker", "[", "Red", "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", "0.04"}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"S", ",", "0", ",", "0"}], "}"}], ",", "0.1"}], "]"}]}], "}"}]}], "}"}], "]"}], "/.", " ", "u"}]}], ";"}], "\n", RowBox[{" ", RowBox[{ StyleBox["Show", FontSize->12], StyleBox["[", FontSize->12], RowBox[{ StyleBox["plot2", FontSize->12], StyleBox[",", FontSize->12], StyleBox[" ", FontSize->12], RowBox[{ StyleBox["{", FontSize->12], RowBox[{ StyleBox[ RowBox[{"AspectRatio", "->", "Automatic"}], FontSize->12], StyleBox[",", FontSize->12], StyleBox[" ", FontSize->12], RowBox[{ StyleBox["Boxed", FontSize->12], StyleBox[" ", FontSize->12], StyleBox["\[Rule]", FontSize->12], " ", "False"}]}], StyleBox["}", FontSize->12]}]}], StyleBox["]", FontSize->12]}], " "}]}], "Input", CellChangeTimes->{{3.5659560360121293`*^9, 3.5659561637451534`*^9}, { 3.565956202838822*^9, 3.565956217518448*^9}, {3.565956255629315*^9, 3.565956259997322*^9}, {3.5659563161886215`*^9, 3.5659563818803368`*^9}, { 3.565956411894789*^9, 3.5659564417532415`*^9}, {3.5659564820793123`*^9, 3.565956572902672*^9}, 3.565956629218771*^9, {3.5659567525369873`*^9, 3.5659567615226035`*^9}, {3.5659568000546713`*^9, 3.565956812035492*^9}}], Cell[BoxData[ Graphics3DBox[{ {Opacity[0.6], SphereBox[{-1.4939437503693909`, 0, 0}, 0.5571059532347319]}, {RGBColor[0, 0, 0], SphereBox[{-1.4939437503693909`, 0, 0}, 0.1]}, {Opacity[0.4], SphereBox[{0, 0, 0}, 2.9878874996438496`]}, {RGBColor[ NCache[ Rational[2, 3], 0.6666666666666666], 0, 0], SphereBox[{0, 0, 0}, 0.04], SphereBox[{1.4939437503693915`, 0, 0}, 0.1]}}, AspectRatio->Automatic, AutomaticImageSize->True, Boxed->False, ImageSize->{346.6398050027404, 394.85696195823203`}, ViewPoint->{0.009656280498412662, -2.2041657841858426`, 2.5674033481460117`}, ViewVertical->{0.02448601044034235, 0.22681482941647008`, 0.9736300470145183}]], "Output", CellChangeTimes->{ 3.5659561687215624`*^9, 3.5659562233372583`*^9, 3.565956264817731*^9, 3.56595634938548*^9, 3.565956388853549*^9, {3.565956422534008*^9, 3.565956450957258*^9}, {3.5659565651338587`*^9, 3.565956581186287*^9}, 3.565956774174226*^9, {3.56595680518708*^9, 3.5659568168247004`*^9}}] }, Open ]] }, WindowSize->{776, 740}, WindowMargins->{{356, Automatic}, {Automatic, 12}}, PrintingCopies->1, PrintingPageRange->{Automatic, Automatic}, PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}}, CellLabelAutoDelete->True, 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[557, 20, 908, 24, 112, "Input"], Cell[1468, 46, 69, 2, 31, "Input"], Cell[1540, 50, 254, 9, 31, "Input"], Cell[1797, 61, 803, 28, 52, "Input"], Cell[2603, 91, 1212, 41, 72, "Input"], Cell[3818, 134, 144, 5, 31, "Input"], Cell[3965, 141, 237, 10, 31, "Input"], Cell[4205, 153, 710, 19, 52, "Input"], Cell[4918, 174, 161, 5, 31, "Input"], Cell[CellGroupData[{ Cell[5104, 183, 376, 11, 31, "Input"], Cell[5483, 196, 753, 16, 31, "Output"] }, Open ]], Cell[6251, 215, 121, 4, 31, "Input"], Cell[CellGroupData[{ Cell[6397, 223, 57, 1, 31, "Input"], Cell[6457, 226, 420, 7, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6914, 238, 92, 2, 31, "Input"], Cell[7009, 242, 408, 6, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7454, 253, 62, 1, 31, "Input"], Cell[7519, 256, 432, 7, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7988, 268, 62, 1, 31, "Input"], Cell[8053, 271, 416, 6, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8506, 282, 336, 12, 52, "Input"], Cell[8845, 296, 409, 6, 31, "Output"], Cell[9257, 304, 410, 6, 31, "Output"] }, Open ]], Cell[9682, 313, 1364, 41, 72, "Input"], Cell[CellGroupData[{ Cell[11071, 358, 641, 19, 52, "Input"], Cell[11715, 379, 804, 17, 375, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[12556, 401, 2745, 87, 132, "Input"], Cell[15304, 490, 1020, 21, 410, "Output"] }, Open ]] } ] *) (* End of internal cache information *)