(* 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[ 59003, 1575] NotebookOptionsPosition[ 57103, 1514] NotebookOutlinePosition[ 57751, 1536] CellTagsIndexPosition[ 57708, 1533] 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["5", "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}, { 3.5663654716175947`*^9, 3.566365471929596*^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", "]"}]}], ";", " ", RowBox[{"rad", "=", RowBox[{"180", "/", "\[Pi]"}]}], ";"}], "\[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[{"L", 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}, {3.580635060257883*^9, 3.5806350659830933`*^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.566316505903453*^9, 3.566316526277089*^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]", " ", "\"\\""}]}], "]"}]}]}], "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}, 3.5663654313383245`*^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.580635382803105*^9, 3.5806358068118496`*^9, 3.580635888743193*^9, 3.580635920520449*^9, 3.5806359970385838`*^9, 3.5806360453518686`*^9, 3.5806362106809587`*^9, {3.5806364952790585`*^9, 3.580636513687091*^9}, 3.5806365643403797`*^9, 3.5806366338853025`*^9, 3.580636674835374*^9, 3.5806367632407293`*^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"}]}], "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}, 3.5663654431007447`*^9, {3.5806347422981243`*^9, 3.580634799175824*^9}, { 3.5806348635259376`*^9, 3.5806348661935425`*^9}, {3.5806349155520287`*^9, 3.5806349190464354`*^9}, 3.5806349624769115`*^9, {3.5806352527768764`*^9, 3.5806352824949284`*^9}}], Cell[BoxData[ RowBox[{"-", "198.35190730103852`"}]], "Output", CellChangeTimes->{ 3.580635382803105*^9, 3.5806358068274493`*^9, 3.580635888743193*^9, 3.580635920536049*^9, 3.5806359970541835`*^9, 3.5806360453674684`*^9, 3.5806362106809587`*^9, {3.5806364952946587`*^9, 3.580636513702691*^9}, 3.58063656435598*^9, 3.5806366339009027`*^9, 3.580636674835374*^9, 3.5806367632407293`*^9}], Cell[BoxData["36.269025578002974`"], "Output", CellChangeTimes->{ 3.580635382803105*^9, 3.5806358068274493`*^9, 3.580635888743193*^9, 3.580635920536049*^9, 3.5806359970541835`*^9, 3.5806360453674684`*^9, 3.5806362106809587`*^9, {3.5806364952946587`*^9, 3.580636513702691*^9}, 3.58063656435598*^9, 3.5806366339009027`*^9, 3.580636674835374*^9, 3.5806367632407293`*^9}], Cell[BoxData["9.319335984458716`"], "Output", CellChangeTimes->{ 3.580635382803105*^9, 3.5806358068274493`*^9, 3.580635888743193*^9, 3.580635920536049*^9, 3.5806359970541835`*^9, 3.5806360453674684`*^9, 3.5806362106809587`*^9, {3.5806364952946587`*^9, 3.580636513702691*^9}, 3.58063656435598*^9, 3.5806366339009027`*^9, 3.580636674835374*^9, 3.5806367632407293`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"politzerratio", "=", RowBox[{ RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"Vne", "+", RowBox[{"2", "*", "Vnn"}]}], ")"}]}], "/", RowBox[{"(", RowBox[{"Vee", "-", "Vnn"}], ")"}]}], " ", "/.", " ", "u"}], " ", "/.", " ", "par"}]}]], "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}, 3.5663654431007447`*^9, {3.5806347422981243`*^9, 3.580634799175824*^9}, { 3.5806348635259376`*^9, 3.5806348661935425`*^9}, {3.5806349155520287`*^9, 3.5806349190464354`*^9}, 3.5806349624769115`*^9, {3.5806352527768764`*^9, 3.5806352772377195`*^9}, {3.5806353401058297`*^9, 3.5806353677334785`*^9}}], Cell[BoxData["6.668471438541948`"], "Output", CellChangeTimes->{ 3.5806353828187046`*^9, 3.5806358068430495`*^9, 3.580635888758793*^9, 3.580635920536049*^9, 3.5806359970697837`*^9, 3.5806360453674684`*^9, 3.580636210696559*^9, {3.5806364952946587`*^9, 3.5806365137182913`*^9}, 3.58063656435598*^9, 3.5806366339165025`*^9, 3.580636674850974*^9, 3.5806367632563295`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ 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[{"rad", "*", "\[Beta]"}], "/.", "u"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"angle", " ", "LP"}], "-", "N", "-", "LP"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"rad", "*", "\[Alpha]"}], " ", "/.", " ", "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}, 3.5663654431007447`*^9, {3.5806347422981243`*^9, 3.580634758522153*^9}, { 3.5806350839075246`*^9, 3.580635093563942*^9}}], Cell[BoxData[ RowBox[{"-", "2.`"}]], "Output", CellChangeTimes->{ 3.580635382834305*^9, 3.5806358068430495`*^9, 3.580635888774393*^9, 3.580635920551649*^9, 3.5806359970697837`*^9, 3.580636045383068*^9, 3.580636210712159*^9, {3.580636495310259*^9, 3.5806365137182913`*^9}, 3.58063656437158*^9, 3.5806366339165025`*^9, 3.5806366748665743`*^9, 3.5806367632719297`*^9}], Cell[BoxData["0.9447147056716415`"], "Output", CellChangeTimes->{ 3.580635382834305*^9, 3.5806358068430495`*^9, 3.580635888774393*^9, 3.580635920551649*^9, 3.5806359970697837`*^9, 3.580636045383068*^9, 3.580636210712159*^9, {3.580636495310259*^9, 3.5806365137182913`*^9}, 3.58063656437158*^9, 3.5806366339165025`*^9, 3.5806366748665743`*^9, 3.5806367632719297`*^9}], Cell[BoxData["103.34396659114164`"], "Output", CellChangeTimes->{ 3.580635382834305*^9, 3.5806358068430495`*^9, 3.580635888774393*^9, 3.580635920551649*^9, 3.5806359970697837`*^9, 3.580636045383068*^9, 3.580636210712159*^9, {3.580636495310259*^9, 3.5806365137182913`*^9}, 3.58063656437158*^9, 3.5806366339165025`*^9, 3.5806366748665743`*^9, 3.5806367632719297`*^9}], Cell[BoxData["120.96554554624991`"], "Output", CellChangeTimes->{ 3.580635382834305*^9, 3.5806358068430495`*^9, 3.580635888774393*^9, 3.580635920551649*^9, 3.5806359970697837`*^9, 3.580636045383068*^9, 3.580636210712159*^9, {3.580636495310259*^9, 3.5806365137182913`*^9}, 3.58063656437158*^9, 3.5806366339165025`*^9, 3.5806366748665743`*^9, 3.5806367632719297`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"d1", " ", "/.", "u"}], "\[IndentingNewLine]", RowBox[{"d3", " ", "/.", "u"}], "\[IndentingNewLine]", RowBox[{"d5", " ", "/.", "u"}], "\[IndentingNewLine]", RowBox[{"d2", " ", "/.", "u"}], "\[IndentingNewLine]", RowBox[{"d4", " ", "/.", "u"}]}], "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}, 3.5663654431007447`*^9, {3.5806347422981243`*^9, 3.580634758522153*^9}}], Cell[BoxData["1.8756856345332227`"], "Output", CellChangeTimes->{ 3.580635382849905*^9, 3.5806358068586493`*^9, 3.5806358887899933`*^9, 3.5806359205672493`*^9, 3.5806359970853834`*^9, 3.5806360453986683`*^9, 3.580636210727759*^9, {3.580636495325859*^9, 3.5806365137338915`*^9}, 3.5806365643871803`*^9, 3.5806366339477024`*^9, 3.5806366748821745`*^9, 3.5806367632719297`*^9}], Cell[BoxData["2.8009659614062543`"], "Output", CellChangeTimes->{ 3.580635382849905*^9, 3.5806358068586493`*^9, 3.5806358887899933`*^9, 3.5806359205672493`*^9, 3.5806359970853834`*^9, 3.5806360453986683`*^9, 3.580636210727759*^9, {3.580636495325859*^9, 3.5806365137338915`*^9}, 3.5806365643871803`*^9, 3.5806366339477024`*^9, 3.5806366748821745`*^9, 3.58063676328753*^9}], Cell[BoxData["2.366755647766626`"], "Output", CellChangeTimes->{ 3.580635382849905*^9, 3.5806358068586493`*^9, 3.5806358887899933`*^9, 3.5806359205672493`*^9, 3.5806359970853834`*^9, 3.5806360453986683`*^9, 3.580636210727759*^9, {3.580636495325859*^9, 3.5806365137338915`*^9}, 3.5806365643871803`*^9, 3.5806366339477024`*^9, 3.5806366748821745`*^9, 3.58063676328753*^9}], Cell[BoxData["2.019953840865406`"], "Output", CellChangeTimes->{ 3.580635382849905*^9, 3.5806358068586493`*^9, 3.5806358887899933`*^9, 3.5806359205672493`*^9, 3.5806359970853834`*^9, 3.5806360453986683`*^9, 3.580636210727759*^9, {3.580636495325859*^9, 3.5806365137338915`*^9}, 3.5806365643871803`*^9, 3.5806366339477024`*^9, 3.5806366748821745`*^9, 3.58063676328753*^9}], Cell[BoxData["2.5115384660834215`"], "Output", CellChangeTimes->{ 3.580635382849905*^9, 3.5806358068586493`*^9, 3.5806358887899933`*^9, 3.5806359205672493`*^9, 3.5806359970853834`*^9, 3.5806360453986683`*^9, 3.580636210727759*^9, {3.580636495325859*^9, 3.5806365137338915`*^9}, 3.5806365643871803`*^9, 3.5806366339477024`*^9, 3.5806366748821745`*^9, 3.58063676328753*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"Projection", " ", "onto", " ", "HOH"}], "-", "plane"}], ",", " ", RowBox[{"LonePairs", " ", "dashed"}]}], " ", "*)"}]], "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.5663626853591013`*^9, 3.5663626922387133`*^9}, { 3.5663637103901014`*^9, 3.566363711466503*^9}, {3.566365546279326*^9, 3.5663655869329977`*^9}, 3.566369639264944*^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[{"{", RowBox[{ RowBox[{"Dashing", "[", RowBox[{"{", RowBox[{"0.02", ",", "0.02"}], "}"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}], "*", RowBox[{"Cos", "[", RowBox[{"\[Alpha]", "/", "2"}], "]"}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}], ",", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "*", RowBox[{"Cos", "[", RowBox[{"\[Alpha]", "/", "2"}], "]"}], "*", RowBox[{"Sin", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}]}], "}"}], ",", "R"}], "]"}]}], "}"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}], "*", RowBox[{"Sin", "[", RowBox[{"\[Beta]", "-", RowBox[{"\[Pi]", "/", "2"}]}], "]"}]}], ",", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "-", RowBox[{"\[Pi]", "/", "2"}]}], "]"}]}]}], "}"}], ",", "Q"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"P", "+", "Q"}], ",", "0"}], "}"}], ",", "Q"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}], "*", RowBox[{"Sin", "[", RowBox[{"\[Beta]", "-", RowBox[{"\[Pi]", "/", "2"}]}], "]"}]}], ",", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "-", RowBox[{"\[Pi]", "/", "2"}]}], "]"}]}]}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", 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[{"P", "+", "Q", "+", "S"}], ",", "0"}], "}"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Dashing", "[", RowBox[{"{", RowBox[{"0.01", ",", "0.01"}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}], "*", RowBox[{"Cos", "[", RowBox[{"\[Alpha]", "/", "2"}], "]"}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}], ",", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "*", RowBox[{"Cos", "[", RowBox[{"\[Alpha]", "/", "2"}], "]"}], "*", RowBox[{"Sin", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}]}], "}"}]}], "}"}], "]"}]}], "}"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}], "*", RowBox[{"Sin", "[", RowBox[{"\[Beta]", "-", RowBox[{"\[Pi]", "/", "2"}]}], "]"}]}], ",", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "-", RowBox[{"\[Pi]", "/", "2"}]}], "]"}]}]}], "}"}]}], "}"}], "]"}]}], "}"}]}], "}"}], "]"}], "/.", "u"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plot1", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "2"}], ",", "2.5"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "2.5"}], ",", "2"}], "}"}]}], "}"}]}], ",", 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, {3.5663176430442986`*^9, 3.566317682730768*^9}, {3.5663204234522157`*^9, 3.566320431267829*^9}, { 3.5663204642618876`*^9, 3.5663204702366977`*^9}, {3.566362812374524*^9, 3.5663628459301834`*^9}, {3.5663629092506943`*^9, 3.566362925334323*^9}, { 3.5663630995242286`*^9, 3.5663631074646425`*^9}, {3.566363142408704*^9, 3.5663631503023176`*^9}, {3.566363335505843*^9, 3.5663634250188007`*^9}, { 3.5663635259509773`*^9, 3.566363546730214*^9}, 3.5663635877426863`*^9, { 3.5663641657237015`*^9, 3.5663642640662737`*^9}, 3.5663644285061626`*^9, 3.5663644816710563`*^9, {3.5663645292199397`*^9, 3.566364585130438*^9}, { 3.5663650531156597`*^9, 3.566365054753663*^9}, {3.566365147714226*^9, 3.5663651764494762`*^9}, {3.566365323167734*^9, 3.566365365833809*^9}}], Cell[BoxData[ GraphicsBox[{CircleBox[{0, 0}, 0.17819997281079103`], DiskBox[{0, 0}, 0.1], {Dashing[{0.02, 0.02}], CircleBox[{-0.39823926301175055`, 0.5037517326518361}, 1.12517354005205]}, CircleBox[{-0.2759186527435787, -1.1632302685263332`}, 1.0173065142262625`], CircleBox[{1.1955064870370535`, 0}, 1.0173065142262625`], DiskBox[{-0.41203000130891465`, -1.7370546148211736`}, 0.1], DiskBox[{1.7852527711363901`, 0}, 0.1], {Thickness[0.012], LineBox[{{0, 0}, {1.7852527711363901`, 0}}], {Dashing[{0.01, 0.01}], LineBox[{{0, 0}, {-0.39823926301175055`, 0.5037517326518361}}]}, LineBox[{{0, 0}, {-0.41203000130891465`, -1.7370546148211736`}}]}}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, PlotRange->{{-2, 2.5}, {-2.5, 2}}]], "Output", CellChangeTimes->{ 3.5806353828811045`*^9, 3.5806358068898497`*^9, 3.5806358888055935`*^9, 3.5806359205828495`*^9, 3.580635997116584*^9, 3.5806360454142685`*^9, 3.580636210758959*^9, {3.5806364953570585`*^9, 3.580636513765091*^9}, 3.5806365644027805`*^9, 3.5806366339633026`*^9, 3.5806366748977747`*^9, 3.5806367633031297`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"plot3", " ", "=", " ", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.5", "]"}], ",", RowBox[{"Lighter", "[", "Red", "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", " ", "P"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Darker", "[", RowBox[{"Red", ",", "1.0"}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", " ", "0.08"}], "]"}]}], "}"}], ",", RowBox[{"Opacity", "[", "0.6", "]"}], ",", RowBox[{"Lighter", "[", RowBox[{"Red", ",", "0.6"}], "]"}], ",", " ", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "*", RowBox[{"Sin", "[", RowBox[{"\[Alpha]", "/", "2"}], "]"}]}], ",", "0", ",", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "*", RowBox[{"Cos", "[", RowBox[{"\[Alpha]", "/", "2"}], "]"}]}]}], "}"}], ",", " ", "R"}], "]"}], ",", " ", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}], "*", RowBox[{"Sin", "[", RowBox[{"\[Alpha]", "/", "2"}], "]"}]}], ",", "0", ",", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "*", RowBox[{"Cos", "[", RowBox[{"\[Alpha]", "/", "2"}], "]"}]}]}], "}"}], ",", " ", "R"}], "]"}], ",", RowBox[{"Opacity", "[", "0.4", "]"}], ",", RowBox[{"Lighter", "[", RowBox[{"Cyan", ",", "0.4"}], "]"}], ",", " ", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}], "*", RowBox[{"Sin", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}], ",", " ", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}]}], "}"}], ",", " ", "Q"}], "]"}], ",", " ", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], "*", RowBox[{"Sin", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}], ",", " ", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}]}], "}"}], ",", " ", "Q"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Darker", "[", "Red", "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q", "+", "S"}], ")"}], "*", RowBox[{"Sin", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}], ",", " ", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q", "+", "S"}], ")"}]}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}]}], "}"}], ",", " ", "0.06"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q", "+", "S"}], ")"}]}], "*", RowBox[{"Sin", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}], ",", " ", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q", "+", "S"}], ")"}]}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}]}], "}"}], ",", " ", "0.06"}], "]"}]}], "}"}]}], "}"}], "]"}], " ", "/.", " ", "u"}]}], ";"}], " "}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plot3", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "->", " ", "Automatic"}], ",", " ", RowBox[{"Boxed", " ", "\[Rule]", " ", "False"}], ",", RowBox[{"SphericalRegion", "\[Rule]", "True"}]}], "}"}]}], "]"}]}], "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}, { 3.566317471849598*^9, 3.5663175599429526`*^9}, {3.5663175959010158`*^9, 3.56631760969144*^9}, {3.566321253732474*^9, 3.5663213099081726`*^9}, { 3.5663213595474596`*^9, 3.5663214564548297`*^9}, {3.5663214993705053`*^9, 3.566321505719716*^9}, {3.5663216132507052`*^9, 3.5663216271815295`*^9}, { 3.56632166145479*^9, 3.566321670892806*^9}, 3.5663217012504597`*^9, { 3.566321797674229*^9, 3.566321848187118*^9}, 3.566370359268608*^9, 3.5806358020694413`*^9, {3.580635868338358*^9, 3.5806358833299837`*^9}, 3.5806359140620375`*^9, {3.580635957008913*^9, 3.580635975104945*^9}, 3.5806360369902534`*^9, {3.5806360910443487`*^9, 3.580636112728387*^9}, { 3.5806361975457354`*^9, 3.580636407591305*^9}, {3.580636441427764*^9, 3.5806366255080876`*^9}, {3.580636657113743*^9, 3.5806366663177595`*^9}, { 3.580636733195077*^9, 3.5806367571099186`*^9}}], Cell[BoxData[ Graphics3DBox[{{ {RGBColor[1, NCache[ Rational[1, 3], 0.3333333333333333], NCache[ Rational[1, 3], 0.3333333333333333]], Opacity[0.5], SphereBox[{0, 0, 0}, 0.17819997281079103`]}, {RGBColor[0., 0., 0.], SphereBox[{0, 0, 0}, 0.08]}, {RGBColor[1., 0.6, 0.6], Opacity[0.6], SphereBox[{1.1342055348474336`, 0, 0.6421528780235041}, 1.12517354005205], SphereBox[{-1.1342055348474336`, 0, 0.6421528780235041}, 1.12517354005205]}, {RGBColor[0.4, 1., 1.], Opacity[0.4], SphereBox[{0, 0.9378428172666112, -0.7414085315459367}, 1.0173065142262625`], SphereBox[{0, -0.9378428172666112, -0.7414085315459367}, 1.0173065142262625`]}}, {RGBColor[ NCache[ Rational[2, 3], 0.6666666666666666], 0, 0], SphereBox[{0, 1.400482980703127, -1.107147179742172}, 0.06], SphereBox[{0, -1.400482980703127, -1.107147179742172}, 0.06]}}, AspectRatio->Automatic, Boxed->False, ImageSize->{574.1074710560829, 574.}, SphericalRegion->True, ViewPoint->{ 1.8072134919722702`, -2.8515488343281934`, -0.22945291428727055`}, ViewVertical->{0.14517716108713158`, -0.7720450087731019, 0.9352443990166104}]], "Output", CellChangeTimes->{ 3.5806360454454684`*^9, 3.580636210774559*^9, {3.5806364953726587`*^9, 3.580636513780691*^9}, 3.5806365644339805`*^9, 3.5806366339789023`*^9, 3.5806366749133744`*^9, 3.5806367633187294`*^9}] }, Open ]] }, WindowSize->{829, 847}, WindowMargins->{{255, Automatic}, {5, Automatic}}, PrintingCopies->1, PrintingPageRange->{32000, 32000}, PrintingOptions->{"Magnification"->1., "PaperOrientation"->"Portrait", "PaperSize"->{595.2, 841.8}}, PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}}, FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)", StyleDefinitions->FrontEnd`FileName[{"Utility"}, "Notepad.nb", CharacterEncoding -> "WindowsANSI"] ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[557, 20, 13470, 397, 194, "Input"], Cell[14030, 419, 10608, 320, 1052, "Input"], Cell[24641, 741, 1966, 40, 112, "Input"], Cell[26610, 783, 1584, 36, 92, "Input"], Cell[CellGroupData[{ Cell[28219, 823, 1301, 30, 72, "Input"], Cell[29523, 855, 855, 18, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[30415, 878, 950, 17, 92, "Input"], Cell[31368, 897, 403, 7, 30, "Output"], Cell[31774, 906, 386, 6, 30, "Output"], Cell[32163, 914, 385, 6, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[32585, 925, 1019, 20, 31, "Input"], Cell[33607, 947, 387, 6, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[34031, 958, 1279, 31, 152, "Input"], Cell[35313, 991, 385, 7, 30, "Output"], Cell[35701, 1000, 384, 6, 30, "Output"], Cell[36088, 1008, 384, 6, 30, "Output"], Cell[36475, 1016, 384, 6, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[36896, 1027, 742, 11, 112, "Input"], Cell[37641, 1040, 393, 6, 30, "Output"], Cell[38037, 1048, 390, 6, 30, "Output"], Cell[38430, 1056, 389, 6, 30, "Output"], Cell[38822, 1064, 389, 6, 30, "Output"], Cell[39214, 1072, 390, 6, 30, "Output"] }, Open ]], Cell[39619, 1081, 707, 12, 31, "Input"], Cell[CellGroupData[{ Cell[40351, 1097, 7121, 184, 252, "Input"], Cell[47475, 1283, 1171, 23, 381, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[48683, 1311, 6939, 163, 212, "Input"], Cell[55625, 1476, 1462, 35, 621, "Output"] }, Open ]] } ] *) (* End of internal cache information *)