(* 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[ 56524, 1604] NotebookOptionsPosition[ 55001, 1554] NotebookOutlinePosition[ 55346, 1569] CellTagsIndexPosition[ 55303, 1566] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ StyleBox[ RowBox[{"(*", " ", RowBox[{ RowBox[{ StyleBox[ RowBox[{ StyleBox["H", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "2O"}]], " ", "with", " ", RowBox[{ StyleBox[ StyleBox[ StyleBox[ StyleBox["O", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "[", RowBox[{ RowBox[{"1", "s2", " ", StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["2", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "sp"}], ",", " ", RowBox[{ StyleBox[ StyleBox[ StyleBox[ StyleBox["2", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "sp", " ", "LP"}]}], "]"}]}], ",", RowBox[{ RowBox[{"~", "no"}], " ", "overlap"}], ",", " ", RowBox[{"parametrized", " ", "ES", " ", StyleBox[ RowBox[{"0", StyleBox[ StyleBox[ StyleBox["4", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], ".01"}]], ".2013", "\[IndentingNewLine]", RowBox[{"Variables", ":", "\[IndentingNewLine]", StyleBox[ RowBox[{" ", StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]], StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], RowBox[{ RowBox[{ StyleBox[ RowBox[{" ", " "}], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], StyleBox[ RowBox[{ StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], " "}]], RowBox[{"P", StyleBox[ RowBox[{" ", StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]], StyleBox[ StyleBox[ StyleBox[ StyleBox["O", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]}], "-", RowBox[{"core", " ", "radius", "\[IndentingNewLine]", "Q", StyleBox[ RowBox[{" ", StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]], StyleBox[ StyleBox[ StyleBox[ StyleBox["O", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}], "-", RowBox[{"H", " ", "cloud", " ", "radius", "\[IndentingNewLine]", "R", StyleBox[ RowBox[{" ", StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]], "lone", " ", "pair", " ", "radius", "\[IndentingNewLine]", "S", StyleBox[ RowBox[{" ", StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]], "distance", " ", "proton", " ", "to", " ", "center", " ", "of", " ", StyleBox[ RowBox[{ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["O", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "H"}]], " ", "cloud", "\[IndentingNewLine]", "\[Alpha]", StyleBox[ RowBox[{" ", StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]], "angle", " ", "lonePair"}], "-", StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["O", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "-", RowBox[{ StyleBox[ RowBox[{ StyleBox[ StyleBox[ StyleBox[ StyleBox["l", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["o", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "n", StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["e", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Pa", StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["i", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["r", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]], "\[IndentingNewLine]", "\[Beta]", StyleBox[ RowBox[{" ", StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]], "angle", " ", StyleBox[ StyleBox[ StyleBox[ StyleBox["H", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}], "-", StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["O", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "-", "H"}]}]}]}], " ", "*)"}], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, 3.566143344240477*^9, 3.566143438043442*^9, 3.566143491551536*^9, {3.5661922351554785`*^9, 3.5661923427020674`*^9}, {3.566201924306572*^9, 3.566201971387455*^9}, 3.566202322310071*^9, {3.5663083624775496`*^9, 3.5663083893563967`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Clear", "[", RowBox[{ "wc", ",", "k1", ",", "k2", ",", "k3", ",", "s1", ",", "s2", ",", "s3", ",", "par"}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Z", "=", "8.0"}], ";", " ", RowBox[{"w3", "=", RowBox[{"Sqrt", "[", "3", "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", "core", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Tc", "=", RowBox[{"2.25", "*", RowBox[{"k1", "/", RowBox[{"P", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vnec", "=", RowBox[{ RowBox[{"-", "3.0"}], "*", RowBox[{"Z", "/", "P"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veec", "=", RowBox[{"3.0", "*", RowBox[{"s1", "/", "P"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"XH", " ", "bond"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Tb", "=", RowBox[{"2.25", "*", RowBox[{"k2", "/", RowBox[{"Q", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vneb", "=", RowBox[{ RowBox[{"-", "2.0"}], "*", RowBox[{ RowBox[{"(", RowBox[{"1.5", "-", RowBox[{"0.5", "*", RowBox[{ RowBox[{"(", RowBox[{"S", "/", "Q"}], ")"}], "^", "2"}]}]}], ")"}], "/", "Q"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veeb", "=", RowBox[{"3.0", "*", RowBox[{"s2", "/", "Q"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"lone", " ", "pair"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Tp", "=", RowBox[{"2.25", "*", RowBox[{"k3", "/", RowBox[{"R", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veep", "=", RowBox[{"3.0", "*", RowBox[{"s3", "/", "R"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"core", "-", "bond"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vnecb", "=", RowBox[{ RowBox[{ RowBox[{"-", "2.0"}], "*", RowBox[{"Z", "/", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}]}], "-", RowBox[{"2", "/", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Veecb", "=", RowBox[{"4.0", "/", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Vnncb", "=", RowBox[{"Z", "/", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"core", "-", RowBox[{"lone", " ", "pair"}]}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vnecp", "=", RowBox[{ RowBox[{"-", "2.0"}], "*", RowBox[{"Z", "/", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veecp", "=", RowBox[{"4.0", "/", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"bond", "-", "bond"}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"d1", " ", "w", RowBox[{"(", "NH", ")"}]}], "-", RowBox[{"w", RowBox[{"(", "NH", ")"}]}]}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d1", "=", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "*", RowBox[{"Sqrt", "[", RowBox[{"2", "*", RowBox[{"(", RowBox[{"1", "-", RowBox[{"Cos", "[", "\[Beta]", "]"}]}], ")"}]}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veebb", "=", RowBox[{"4.0", "/", "d1"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"d3", " ", "H"}], "-", "H"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d3", "=", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "*", RowBox[{"Sqrt", "[", RowBox[{"2", "*", RowBox[{"(", RowBox[{"1", "-", RowBox[{"Cos", "[", "\[Beta]", "]"}]}], ")"}]}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Vnnbb", "=", RowBox[{"1.0", "/", "d3"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"d5", " ", "w", RowBox[{"(", "OH", ")"}]}], "-", "H"}], ",", " ", RowBox[{"\[Beta]", " ", "is", " ", "angle", " ", "HOH"}]}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d5", "=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "^", "2"}], "-", RowBox[{"2", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "*", RowBox[{"Cos", "[", "\[Beta]", "]"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Vnebb", "=", RowBox[{ RowBox[{"-", "2.0"}], "/", "d5"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"bond", "-", RowBox[{"lone", " ", "pair"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"d2", " ", "w", RowBox[{"(", "OH", ")"}]}], "-", RowBox[{"w", RowBox[{"(", "LP", ")"}]}]}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d2", "=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "^", "2"}], "+", RowBox[{"2", "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "*", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "*", RowBox[{"Cos", "[", RowBox[{"\[Alpha]", "/", "2"}], "]"}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veebp", " ", "=", " ", RowBox[{"4", "/", "d2"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"d4", " ", "w", RowBox[{"(", "LP", ")"}]}], "-", "H"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d4", "=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "^", "2"}], "+", RowBox[{"2", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "*", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "*", RowBox[{"Cos", "[", RowBox[{"\[Alpha]", "/", "2"}], "]"}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Vnebp", " ", "=", RowBox[{ RowBox[{"-", "2.0"}], "/", "d4"}]}], " ", ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", StyleBox[ RowBox[{"(*", RowBox[{ RowBox[{"lone", " ", "pair"}], "-", RowBox[{"lone", " ", "pair"}]}], "*)"}], "Text", FontFamily->"Arial Narrow"], "\[IndentingNewLine]", StyleBox[ RowBox[{"(*", " ", RowBox[{ RowBox[{"none", " ", "for", " ", "NH3"}], ",", " ", RowBox[{"but", " ", "for", " ", "H2O", " ", "and", " ", "HF"}]}], " ", "*)"}], "Text", FontFamily->"Arial Narrow"], "\[IndentingNewLine]", StyleBox[ RowBox[{"(*", " ", RowBox[{ RowBox[{ "H2O", " ", "has", " ", "two", " ", "independent", " ", "angles", " ", "as", " ", RowBox[{"variables", ":", " ", RowBox[{"H", "-", "O", "-", RowBox[{"H", " ", "and", " ", "LP"}], "-", "O", "-", StyleBox[ RowBox[{ StyleBox["L", "Text", FontFamily->"Arial Narrow"], StyleBox[ StyleBox["P", "Text", FontFamily->"Arial Narrow"], "Text", FontFamily->"Arial Narrow"]}]]}]}]}], ","}], " ", "*)"}], "Text", FontFamily->"Arial Narrow"]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ StyleBox["d6", "Input"], StyleBox[" ", "Input"], StyleBox["=", "Input"], StyleBox[" ", "Input"], RowBox[{ StyleBox[ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "Input"], "*", RowBox[{"Sqrt", "[", RowBox[{"2", "*", RowBox[{"(", RowBox[{"1", "-", RowBox[{"Cos", "[", "\[Alpha]", "]"}]}], ")"}]}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veepp", "=", RowBox[{"4", "/", "d6"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", RowBox[{"sum", " ", "terms"}], "*)"}]}]}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, 3.566143344240477*^9, 3.566143438043442*^9, {3.5661923476784763`*^9, 3.566192348645678*^9}, { 3.566192420296604*^9, 3.566192420858205*^9}, {3.5661924684694886`*^9, 3.566192482369113*^9}, {3.5661925519608355`*^9, 3.5661925980901165`*^9}, { 3.566192662409029*^9, 3.5661926728298473`*^9}, {3.566194733031866*^9, 3.5661947591463118`*^9}, {3.566201172322851*^9, 3.5662012049425087`*^9}, { 3.566201275782233*^9, 3.5662012764842343`*^9}, {3.566201607563616*^9, 3.5662016091548185`*^9}, {3.5662016860473537`*^9, 3.5662017455770583`*^9}, {3.566201840690425*^9, 3.5662018510176435`*^9}, { 3.5662021072324934`*^9, 3.56620213362774*^9}, {3.566202199116655*^9, 3.566202199381855*^9}, {3.5663049674119854`*^9, 3.5663050328385005`*^9}, { 3.566305081167385*^9, 3.566305085270192*^9}, {3.56630512948067*^9, 3.566305132819076*^9}, {3.5663056216303344`*^9, 3.56630564170757*^9}, 3.566305710550491*^9, {3.5663057513601627`*^9, 3.566305753419366*^9}, { 3.566307371173808*^9, 3.5663074141830835`*^9}, 3.566307457207959*^9, 3.5663085639675035`*^9, {3.5663156573399625`*^9, 3.5663156653739767`*^9}, { 3.5663157659941535`*^9, 3.5663157671173553`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"T", "=", RowBox[{"Tc", "+", " ", RowBox[{"2", "*", "Tb"}], "+", " ", RowBox[{"2", "*", "Tp"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vne", "=", RowBox[{"Vnec", "+", RowBox[{"2", "*", "Vneb"}], "+", " ", RowBox[{"2", "*", "Vnecb"}], "+", RowBox[{"2", "*", "Vnecp"}], "+", " ", RowBox[{"2", "*", "Vnebb"}], "+", RowBox[{"4", "*", "Vnebp"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vee", "=", RowBox[{"Veec", "+", RowBox[{"2", "*", "Veeb"}], "+", " ", RowBox[{"2", "*", "Veep"}], "+", RowBox[{"2", "*", "Veecb"}], "+", RowBox[{"2", "*", "Veecp"}], "+", " ", "Veebb", "+", RowBox[{"4", "*", "Veebp"}], " ", "+", " ", "Veepp"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ "Vnn", "=", " ", RowBox[{ RowBox[{"2", "*", "Vnncb"}], "+", " ", "Vnnbb"}]}], ";"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, 3.566143344240477*^9, { 3.5661435118159714`*^9, 3.5661435187891836`*^9}, {3.5661470921646585`*^9, 3.5661471085134873`*^9}, {3.566194770659132*^9, 3.566194777023943*^9}, { 3.566200993811738*^9, 3.566201133478783*^9}, {3.5662012401361704`*^9, 3.5662012446913786`*^9}, {3.5662014429677267`*^9, 3.566201443654128*^9}, { 3.566201792735941*^9, 3.5662017938591433`*^9}, 3.5662022754943886`*^9, { 3.5662025519112744`*^9, 3.566202565436498*^9}, {3.566305456238844*^9, 3.566305473289674*^9}, {3.5663055060497313`*^9, 3.5663055070949335`*^9}, { 3.566308222654504*^9, 3.5663082634017754`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"Parameters", ":", " ", "G3"}], "//", RowBox[{ RowBox[{ RowBox[{ RowBox[{"6", "-", RowBox[{"311", "G", RowBox[{"(", RowBox[{"d", ",", "p"}], ")"}], " ", "can", " ", "be", " ", "read", " ", "for", " ", "LiH"}]}], " ", "..."}], " ", "HF", "\[IndentingNewLine]", "in", " ", "table", " ", "from", " ", RowBox[{"www", ".", "kimball"}]}], "-", RowBox[{ RowBox[{ RowBox[{"model", ".", "org"}], "/", "Kimball"}], "/", RowBox[{"params", ".", "xls"}]}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"par", "=", RowBox[{"{", RowBox[{ RowBox[{"k1", "\[Rule]", "0.92438558"}], ",", RowBox[{"k2", "\[Rule]", "1.42313501"}], ",", RowBox[{"k3", "\[Rule]", "1.32139847"}], ",", RowBox[{"s1", "\[Rule]", "0.43732345"}], ",", RowBox[{"s2", "\[Rule]", "0.31779336"}], ",", RowBox[{"s3", "\[Rule]", "0.44741336"}]}], "}"}]}], ";"}]}]], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.566143371150524*^9}, {3.5663152612084665`*^9, 3.5663152707400837`*^9}, { 3.566315729474489*^9, 3.566315733312096*^9}, {3.566315805368623*^9, 3.5663158109066324`*^9}, {3.566316282745061*^9, 3.566316287222269*^9}, { 3.5663163274391394`*^9, 3.5663163954864593`*^9}}], Cell[BoxData[ FormBox[ StyleBox[ RowBox[{ StyleBox[ RowBox[{"(", "*"}]], RowBox[{ RowBox[{"par", " ", "=", " ", RowBox[{"{", RowBox[{ RowBox[{"k1", "\[Rule]", "0.92438558"}], ",", RowBox[{"k2", "\[Rule]", "1.42313501"}], ",", RowBox[{"k3", "\[Rule]", "1.32139847"}], ",", RowBox[{"s1", "\[Rule]", "0.43732345"}], ",", RowBox[{"s2", "\[Rule]", "0.31779336"}], ",", RowBox[{"s3", "\[Rule]", "0.44741336"}]}], "}"}]}], ";", "\[IndentingNewLine]", "par"}], StyleBox[ RowBox[{"*", ")"}]]}], "Input"], InputForm]], "Text", CellChangeTimes->{{3.566315254204055*^9, 3.566315321642973*^9}, { 3.566315583583033*^9, 3.566315632286319*^9}, {3.566315713468861*^9, 3.566315719209671*^9}, {3.566315796461007*^9, 3.5663158011722155`*^9}, { 3.566315868345933*^9, 3.566315958139691*^9}, {3.5663160990079384`*^9, 3.5663161075099535`*^9}, {3.566316271918642*^9, 3.5663162785954537`*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"func", "=", RowBox[{ RowBox[{"T", "+", "Vne", "+", "Vee", "+", "Vnn"}], "/.", "par"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"t", "=", RowBox[{"FindMinimum", "[", RowBox[{"func", ",", RowBox[{"{", RowBox[{"P", ",", "0.18"}], "}"}], ",", RowBox[{"{", RowBox[{"Q", ",", "1.0"}], "}"}], ",", RowBox[{"{", RowBox[{"R", ",", " ", "1.1"}], "}"}], ",", RowBox[{"{", RowBox[{"S", ",", "0.6"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"\[Alpha]", ",", "1.97"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Beta]", ",", "1.97"}], "}"}], ",", RowBox[{"MaxIterations", "\[Rule]", "500"}], ",", " ", RowBox[{"Method", " ", "\[Rule]", " ", "\"\\""}]}], "]"}]}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.5661433819769435`*^9}, {3.5662024426798825`*^9, 3.5662024734587364`*^9}, { 3.566305156858718*^9, 3.566305166265535*^9}, {3.5663079842392845`*^9, 3.566307984972486*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "76.3817728692884`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.17819997281079103`"}], ",", RowBox[{"Q", "\[Rule]", "1.0173065142262625`"}], ",", RowBox[{"R", "\[Rule]", "1.12517354005205`"}], ",", RowBox[{"S", "\[Rule]", "0.5897462840993366`"}], ",", RowBox[{"\[Alpha]", "\[Rule]", "2.1112470512532235`"}], ",", RowBox[{"\[Beta]", "\[Rule]", "1.8036924790864421`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.5663163013246937`*^9, 3.566316402428471*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"u", "=", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Vne", "/.", "u"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vee", "/.", "u"}], "/.", "par"}], "\[IndentingNewLine]", RowBox[{"Vnn", "/.", "u"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"Vne", "+", "Vee", "+", "Vnn"}], ")"}], "/", "T"}], "/.", "u"}], "/.", "par"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "*", "0.529177"}], "/.", "u"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"angle", " ", "HNH"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"180", "/", "\[Pi]"}], "*", "\[Beta]"}], "/.", "u"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"angle", " ", "LP"}], "-", "N", "-", "LP"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"180", "/", "\[Pi]"}], "*", "\[Alpha]"}], " ", "/.", " ", "u"}], "\[IndentingNewLine]", RowBox[{"d1", " ", "/.", "u"}], "\[IndentingNewLine]", RowBox[{"d3", " ", "/.", "u"}], "\[IndentingNewLine]", RowBox[{"d5", " ", "/.", "u"}], "\[IndentingNewLine]", RowBox[{"d2", " ", "/.", "u"}], "\[IndentingNewLine]", RowBox[{"d4", " ", "/.", "u"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.5661433961729684`*^9}, {3.5663078436050377`*^9, 3.5663078735570903`*^9}}], Cell[BoxData[ RowBox[{"-", "198.35190730103852`"}]], "Output", CellChangeTimes->{3.566316301340294*^9, 3.5663164024440713`*^9}], Cell[BoxData["36.269025578002974`"], "Output", CellChangeTimes->{3.566316301340294*^9, 3.5663164024440713`*^9}], Cell[BoxData["9.319335984458716`"], "Output", CellChangeTimes->{3.566316301340294*^9, 3.5663164024440713`*^9}], Cell[BoxData[ RowBox[{"-", "2.`"}]], "Output", CellChangeTimes->{3.566316301340294*^9, 3.5663164024440713`*^9}], Cell[BoxData["0.9447147056716415`"], "Output", CellChangeTimes->{3.566316301340294*^9, 3.5663164024440713`*^9}], Cell[BoxData["103.34396659114164`"], "Output", CellChangeTimes->{3.566316301340294*^9, 3.5663164024440713`*^9}], Cell[BoxData["120.96554554624991`"], "Output", CellChangeTimes->{3.566316301340294*^9, 3.5663164024596715`*^9}], Cell[BoxData["1.8756856345332227`"], "Output", CellChangeTimes->{3.566316301340294*^9, 3.5663164024596715`*^9}], Cell[BoxData["2.8009659614062543`"], "Output", CellChangeTimes->{3.566316301340294*^9, 3.5663164024596715`*^9}], Cell[BoxData["2.366755647766626`"], "Output", CellChangeTimes->{3.566316301340294*^9, 3.5663164024596715`*^9}], Cell[BoxData["2.019953840865406`"], "Output", CellChangeTimes->{3.566316301340294*^9, 3.5663164024596715`*^9}], Cell[BoxData["2.5115384660834215`"], "Output", CellChangeTimes->{3.566316301340294*^9, 3.5663164024596715`*^9}] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{"wd", " ", "=", " ", RowBox[{ RowBox[{"Sin", "[", "\[Alpha]", "]"}], " ", "/.", " ", "u"}]}], ";", " ", RowBox[{"wc", "=", RowBox[{ RowBox[{"Cos", "[", "\[Alpha]", "]"}], " ", "/.", " ", "u"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{"Projection", " ", "onto", " ", "HHH"}], "-", "plane"}], " ", "*)"}]}]}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.566143410119393*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"plot1", "=", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "P"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "R"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "/", "2"}]}], ",", RowBox[{"wd", "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "/", "2"}]}]}], "}"}], ",", "R"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "/", "2"}]}], ",", RowBox[{ RowBox[{"-", "wd"}], "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "/", "2"}]}]}], "}"}], ",", "Q"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"wd", "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}], ",", "0"}], "}"}], ",", "Q"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "/", "2"}]}], ",", RowBox[{"wd", "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "/", "2"}]}]}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "/", "2"}]}], ",", RowBox[{ RowBox[{"-", "wd"}], "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "/", "2"}]}]}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"wd", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}], ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Thickness", "[", "0.012", "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"wd", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}], ",", "0"}], "}"}]}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "/", "2"}]}], ",", RowBox[{"wd", "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "/", "2"}]}]}], "}"}]}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "/", "2"}]}], ",", RowBox[{ RowBox[{"-", "wd"}], "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "/", "2"}]}]}], "}"}]}], "}"}], "]"}]}], "}"}]}], "}"}], "]"}], "/.", "u"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plot1", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "3"}], ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "}"}]}], ",", RowBox[{"Axes", "\[Rule]", "True"}], ",", RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", RowBox[{"Frame", "\[Rule]", "True"}]}], "}"}]}], "]"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.566143410119393*^9}, 3.566308742431816*^9}], Cell[BoxData[ GraphicsBox[{CircleBox[{0, 0}, 0.17819997281079103`], DiskBox[{0, 0}, 0.1], CircleBox[{0, 0}, 1.12517354005205], CircleBox[{-0.5125595746552611, 0.8877792252088051}, 1.12517354005205], CircleBox[{-0.5125595746552611, -0.8877792252088051}, 1.0173065142262625`], CircleBox[{1.0251191493105223`, 0}, 1.0173065142262625`], DiskBox[{-0.7654064707701024, 1.3257228958158}, 0.1], DiskBox[{-0.7654064707701024, -1.3257228958158}, 0.1], DiskBox[{1.5308129415402048`, 0}, 0.1], {Thickness[0.012], LineBox[{{0, 0}, {1.5308129415402048`, 0}}], LineBox[{{0, 0}, {-0.7654064707701024, 1.3257228958158}}], LineBox[{{0, 0}, {-0.7654064707701024, -1.3257228958158}}]}}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, PlotRange->{{-3, 3}, {-3, 3}}]], "Output", CellChangeTimes->{3.566316301371494*^9, 3.566316402475271*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"plot3", " ", "=", " ", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", " ", "P"}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", " ", "0.1"}], "]"}], ",", " ", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", RowBox[{"P", "+", "R"}]}], "}"}], ",", " ", "R"}], "]"}], ",", " ", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}], "/", "2"}]}], ",", " ", RowBox[{"wd", "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}], "/", "2"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}]}], "}"}], ",", " ", "R"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.5", "]"}], ",", " ", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}], "/", "2"}]}], ",", " ", RowBox[{ RowBox[{"-", "wd"}], "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}], "/", "2"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}]}], "}"}], ",", " ", "Q"}], "]"}], ",", " ", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"wd", "*", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], ",", " ", "0", ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}]}], "}"}], ",", " ", "Q"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Darker", "[", "Red", "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q", " ", "+", " ", "S"}], ")"}], "/", "2"}]}], ",", " ", RowBox[{"wd", "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q", " ", "+", " ", "S"}], ")"}], "/", "2"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}], "}"}], ",", " ", "0.06"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q", " ", "+", " ", "S"}], ")"}], "/", "2"}]}], ",", " ", RowBox[{ RowBox[{"-", "wd"}], "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q", " ", "+", " ", "S"}], ")"}], "/", "2"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}], "}"}], ",", " ", "0.06"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"wd", "*", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q", " ", "+", " ", "S"}], ")"}]}], ",", " ", "0", ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}], "}"}], ",", " ", "0.06"}], "]"}]}], "}"}]}], "}"}], "]"}], " ", "/.", " ", "u"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plot3", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "->", " ", "Automatic"}], ",", " ", RowBox[{"Boxed", " ", "\[Rule]", " ", "False"}]}], "}"}]}], "]"}]}], "Input", CellChangeTimes->{{3.5661464147646685`*^9, 3.566146439287912*^9}, { 3.566146505400828*^9, 3.566146609468611*^9}, 3.5661466450210733`*^9, { 3.5661467117267904`*^9, 3.5661467122883916`*^9}, {3.566146817869377*^9, 3.5661468378998117`*^9}, {3.5661468705194693`*^9, 3.5661468715646715`*^9}, {3.5661469228107615`*^9, 3.5661469312971764`*^9}, {3.5661469641976337`*^9, 3.5661469774888573`*^9}, {3.5661477373577533`*^9, 3.5661477492774353`*^9}, {3.5663088744236484`*^9, 3.566308893580482*^9}}], Cell[BoxData[ Graphics3DBox[{SphereBox[{0, 0, 0}, 0.17819997281079103`], SphereBox[{0, 0, 0}, 0.1], SphereBox[{0, 0, 1.303373512862841}, 1.12517354005205], SphereBox[{-0.5125595746552611, 0.8877792252088051, -0.615115021979262}, 1.12517354005205], {Opacity[0.5], SphereBox[{-0.5125595746552611, -0.8877792252088051, -0.615115021979262}, 1.0173065142262625`], SphereBox[{1.0251191493105223`, 0, -0.615115021979262}, 1.0173065142262625`]}, {RGBColor[ NCache[ Rational[2, 3], 0.6666666666666666], 0, 0], SphereBox[{-0.7654064707701024, 1.3257228958158, -0.9185527719533515}, 0.06], SphereBox[{-0.7654064707701024, -1.3257228958158, \ -0.9185527719533515}, 0.06], SphereBox[{1.5308129415402048`, 0, -0.9185527719533515}, 0.06]}}, AspectRatio->Automatic, Boxed->False]], "Output", CellChangeTimes->{3.5663163013870935`*^9, 3.5663164024908714`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{"Projection", " ", "onto", " ", "LP"}], "-", "HH", "-", "plane"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.5661434217882133`*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"plot2", "=", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "P"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{"P", "+", "R"}]}], "}"}], ",", "R"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{"P", "+", "R"}]}], "}"}], ",", "0.08"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"wd", "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}]}], "}"}], ",", "Q"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"wd", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Dashing", "[", RowBox[{"{", "0.02", "}"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "/", "2"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}]}], "}"}], ",", "Q"}], "]"}]}], "}"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "/", "2"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Thickness", "[", "0.012", "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"wd", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}], "}"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Dashing", "[", RowBox[{"{", "0.02", "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "/", "2"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}], "}"}]}], "}"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Dashing", "[", RowBox[{"{", RowBox[{"0.02", ",", "0.04"}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", RowBox[{"P", "+", "R"}]}], "}"}]}], "}"}], "]"}]}], "}"}]}], " ", "}"}]}], "}"}], "]"}], "/.", "u"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plot2", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "3"}], ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "}"}]}], ",", RowBox[{"Axes", "\[Rule]", "True"}], ",", RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", RowBox[{"Frame", "\[Rule]", "True"}]}], "}"}]}], "]"}]}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.5661434217882133`*^9}}], Cell[BoxData[ GraphicsBox[{CircleBox[{0, 0}, 0.17819997281079103`], DiskBox[{0, 0}, 0.1], CircleBox[{0, 1.303373512862841}, 1.12517354005205], CircleBox[{0, 1.303373512862841}, 0.08], CircleBox[{1.0251191493105223`, -0.615115021979262}, 1.0173065142262625`], DiskBox[{1.5308129415402048`, -0.9185527719533515}, 0.1], {Dashing[{0.02}], CircleBox[{-0.5125595746552611, -0.615115021979262}, 1.0173065142262625`]}, DiskBox[{-0.7654064707701024, -0.9185527719533515}, 0.1], {Thickness[0.012], LineBox[{{0, 0}, {1.5308129415402048`, -0.9185527719533515}}], {Dashing[{0.02}], LineBox[{{0, 0}, {-0.7654064707701024, -0.9185527719533515}}]}, {Dashing[{0.02, 0.04}], LineBox[{{0, 0}, {0, 1.303373512862841}}]}}}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, PlotRange->{{-3, 3}, {-3, 3}}]], "Output", CellChangeTimes->{3.5663163014026937`*^9, 3.5663164025064716`*^9}] }, Open ]] }, WindowSize->{882, 817}, WindowMargins->{{40, Automatic}, {Automatic, 25}}, 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, 13418, 396, 173, "Input"], Cell[13978, 418, 10577, 320, 1052, "Input"], Cell[24558, 740, 1966, 40, 112, "Input"], Cell[26527, 782, 1637, 37, 92, "Input"], Cell[28167, 821, 987, 23, 50, "Text"], Cell[CellGroupData[{ Cell[29179, 848, 1294, 30, 92, "Input"], Cell[30476, 880, 578, 13, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[31091, 898, 1799, 45, 332, "Input"], Cell[32893, 945, 129, 2, 31, "Output"], Cell[33025, 949, 112, 1, 31, "Output"], Cell[33140, 952, 111, 1, 31, "Output"], Cell[33254, 955, 113, 2, 31, "Output"], Cell[33370, 959, 112, 1, 31, "Output"], Cell[33485, 962, 112, 1, 31, "Output"], Cell[33600, 965, 112, 1, 31, "Output"], Cell[33715, 968, 112, 1, 31, "Output"], Cell[33830, 971, 112, 1, 31, "Output"], Cell[33945, 974, 111, 1, 31, "Output"], Cell[34059, 977, 111, 1, 31, "Output"], Cell[34173, 980, 112, 1, 31, "Output"] }, Open ]], Cell[34300, 984, 751, 18, 52, "Input"], Cell[CellGroupData[{ Cell[35076, 1006, 5737, 165, 232, "Input"], Cell[40816, 1173, 897, 17, 379, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[41750, 1195, 5403, 144, 172, "Input"], Cell[47156, 1341, 922, 20, 430, "Output"] }, Open ]], Cell[48093, 1364, 475, 9, 31, "Input"], Cell[CellGroupData[{ Cell[48593, 1377, 5434, 152, 232, "Input"], Cell[54030, 1531, 955, 20, 411, "Output"] }, Open ]] } ] *) (* End of internal cache information *)