(* 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[ 31314, 886] NotebookOptionsPosition[ 29546, 825] NotebookOutlinePosition[ 30188, 848] CellTagsIndexPosition[ 30145, 845] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{"BH3", " ", "molecule", " ", "with", " ", RowBox[{"B", "[", RowBox[{"1", "s2", " ", "3", "sp"}], "]"}]}], ",", " ", "overlap", ",", " ", RowBox[{"parametrized", " ", "ES", " ", "31.12", ".2012"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"Variables", ":", " ", "P"}], ",", "Q", ",", "R", ",", "S"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"with", " ", "core", " ", StyleBox["exchange", FontColor->RGBColor[1, 0, 0]]}], " ", "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"Clear", "[", RowBox[{ "k1", ",", "k2", ",", "s1", ",", "s2", ",", " ", "Z", ",", " ", "P", ",", "Q", ",", "R", ",", "S", ",", "Vne", ",", "Vee", ",", "Vnn", ",", "Veec", ",", "Vnec", ",", "Veecb", ",", "Vnecb", ",", "Vnncb", ",", "Veebb", ",", "Vnebb", ",", "Vnebb", ",", "d1", ",", "d2", ",", "d3"}], "]"}], ";"}]}]], "Input", CellChangeTimes->{{3.5659373241689606`*^9, 3.5659373285525684`*^9}, { 3.56595464158808*^9, 3.565954646127688*^9}, {3.5659551127557077`*^9, 3.56595511394131*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Z", "=", "5.0"}], ";", " ", RowBox[{"w3", "=", RowBox[{"Sqrt", "[", "3", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"R", "=", "1.11"}], ";", " ", RowBox[{"S", "=", "1.14"}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"B", " ", "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[{"BH", " ", "bond"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Tb", "=", RowBox[{"9.0", "*", RowBox[{"k2", "/", RowBox[{"Q", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vneb", "=", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"2", "/", "Q"}], ")"}]}], "*", RowBox[{"(", RowBox[{"1.5", "-", RowBox[{"0.5", "*", RowBox[{ RowBox[{"(", RowBox[{"S", "/", "Q"}], ")"}], "^", "2"}]}]}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veeb", "=", RowBox[{"3.0", "*", RowBox[{"s2", "/", "Q"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"core", "-", "bond"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vnecb", "=", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"2.0", "*", RowBox[{"Z", "/", "Q"}]}], ")"}]}], "*", RowBox[{"(", RowBox[{"1.5", "-", RowBox[{"0.5", "*", RowBox[{ RowBox[{"(", RowBox[{"R", "/", "Q"}], ")"}], "^", "2"}]}]}], ")"}]}], "-", RowBox[{"2.0", "/", RowBox[{"(", RowBox[{"R", "+", "S"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Veecb", "=", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "Q"}], ")"}], "*", RowBox[{"(", RowBox[{"1.5", "-", RowBox[{"0.5", "*", RowBox[{ RowBox[{"(", RowBox[{"R", "/", "Q"}], ")"}], "^", "2"}]}], "-", RowBox[{"0.3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", "/", "Q"}], ")"}], "^", "2"}]}]}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vex", " ", "=", " ", RowBox[{ RowBox[{"-", StyleBox["1.2", FontColor->RGBColor[1, 0, 0]]}], StyleBox["*", FontColor->RGBColor[1, 0, 0]], StyleBox[ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"P", "/", "Q"}], ")"}], "^", "3"}], "/", "Q"}], FontColor->RGBColor[1, 0, 0]]}]}], StyleBox[";", FontColor->RGBColor[1, 0, 0]]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Veecb", " ", "=", " ", RowBox[{"Veecb", " ", "+", " ", "Vex"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Vnncb", "=", RowBox[{"Z", "/", RowBox[{"(", RowBox[{"R", "+", "S"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"bond", "-", "bond"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d1", "=", RowBox[{ RowBox[{"(", RowBox[{"R", "+", "S"}], ")"}], "*", "w3"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d2", "=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"R", "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"R", "+", "S"}], ")"}], "^", "2"}], "+", RowBox[{"R", "*", RowBox[{"(", RowBox[{"R", "+", "S"}], ")"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d3", "=", RowBox[{"R", "*", "w3"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vnebb", " ", "=", " ", RowBox[{ RowBox[{"-", "4"}], "/", "d2"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Veebb", " ", "=", " ", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "Q"}], ")"}], "*", RowBox[{"(", RowBox[{"1.2", "-", RowBox[{"0.5", "*", RowBox[{ RowBox[{"(", RowBox[{"d3", "/", "Q"}], ")"}], "^", "2"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"3", "/", "16"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"d3", "/", "Q"}], ")"}], "^", "3"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "160"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"d3", "/", "Q"}], ")"}], "^", "5"}]}]}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Vnnbb", " ", "=", " ", RowBox[{"1", "/", "d1"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"sum", " ", "terms"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"T", " ", "=", " ", RowBox[{"Tc", "+", RowBox[{"3", "*", "Tb"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vne", " ", "=", " ", RowBox[{"Vnec", "+", RowBox[{"3", "*", "Vneb"}], "+", RowBox[{"3", "*", "Vnecb"}], "+", RowBox[{"3", "*", "Vnebb"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vee", " ", "=", " ", RowBox[{"Veec", "+", RowBox[{"3", "*", "Veeb"}], "+", RowBox[{"3", "*", "Veecb"}], "+", RowBox[{"3", "*", "Veebb"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vnn", " ", "=", " ", RowBox[{ RowBox[{"3", "*", "Vnncb"}], "+", RowBox[{"3", "*", "Vnnbb"}]}]}], ";"}]}], "Input", CellChangeTimes->{{3.565935567309475*^9, 3.565935568339077*^9}, { 3.565935810966303*^9, 3.565935812432706*^9}, {3.5659543970012503`*^9, 3.565954518566264*^9}, 3.565954614678033*^9, 3.565954679620947*^9, { 3.5659547998971577`*^9, 3.565954822688798*^9}, {3.5659551235353265`*^9, 3.565955133066943*^9}, {3.565955524752431*^9, 3.56595554663927*^9}, { 3.5659792122987747`*^9, 3.5659792139211774`*^9}, {3.565980634288072*^9, 3.565980647797696*^9}, {3.5659810477199984`*^9, 3.565981063039225*^9}, { 3.565981138590158*^9, 3.5659811433481665`*^9}}], Cell[BoxData[ RowBox[{ StyleBox["(*", FontSize->12], StyleBox[" ", FontSize->12], RowBox[{"Parameters", " ", "in", " ", "list", " ", "c"}], StyleBox[" ", FontSize->12], StyleBox["*)", FontSize->12]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{"c", "=", RowBox[{"{", RowBox[{ RowBox[{"k1", " ", "->", " ", "1.0"}], ",", " ", RowBox[{"k2", "->", " ", "1.0"}], ",", " ", RowBox[{"s1", " ", "->", " ", "0.3"}], ",", " ", RowBox[{"s2", " ", "->", " ", "0.3"}]}], "}"}]}], ";"}], " ", "*)"}], " ", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"c", "=", RowBox[{"{", RowBox[{ RowBox[{"k1", " ", "\[Rule]", " ", "1.01613950"}], ",", " ", RowBox[{"k2", " ", "\[Rule]", " ", "0.701"}], ",", " ", RowBox[{"s1", " ", "->", " ", "0.25742249"}], ",", RowBox[{"s2", " ", "->", " ", "0.22958899"}]}], "}"}]}], ";"}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"c", "=", RowBox[{"{", RowBox[{ RowBox[{"k1", " ", "\[Rule]", " ", "1.04011439"}], ",", " ", RowBox[{"k2", " ", "\[Rule]", " ", "0.75422468"}], ",", " ", RowBox[{"s1", " ", "->", " ", "0.21404979"}], ",", RowBox[{"s2", " ", "->", " ", "0.07673395"}]}], "}"}]}], ";"}], "*)"}], "\[IndentingNewLine]", " ", RowBox[{"(*", RowBox[{ RowBox[{"c", "=", RowBox[{"{", RowBox[{ RowBox[{"k1", " ", "\[Rule]", " ", "1.0"}], ",", " ", RowBox[{"k2", " ", "\[Rule]", " ", "0.75"}], ",", " ", RowBox[{"s1", " ", "\[Rule]", " ", "0.3"}], ",", RowBox[{"s2", " ", "\[Rule]", " ", "0.3"}]}], "}"}]}], ";"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"c", "=", RowBox[{"{", RowBox[{ RowBox[{"k1", " ", "\[Rule]", " ", "1.02"}], ",", " ", RowBox[{"k2", " ", "\[Rule]", " ", "0.683"}], ",", " ", RowBox[{"s1", " ", "->", " ", "0.21"}], ",", RowBox[{"s2", " ", "->", " ", "0.21"}]}], "}"}]}], ";"}], "\[IndentingNewLine]"}]}]], "Input", CellChangeTimes->{{3.5659549426218085`*^9, 3.5659549956619015`*^9}, { 3.565955028047559*^9, 3.565955042056383*^9}, 3.5659550753936415`*^9, { 3.56595513716975*^9, 3.565955147216168*^9}, {3.5659551787594233`*^9, 3.5659552181494923`*^9}, {3.5659553061648474`*^9, 3.5659553459917173`*^9}, {3.565955636277027*^9, 3.5659556838415103`*^9}, { 3.5659789796867657`*^9, 3.56597899349279*^9}, {3.565979097591773*^9, 3.5659790989489756`*^9}, {3.565979307724142*^9, 3.5659795037696867`*^9}, { 3.5659796254031*^9, 3.565979640691127*^9}, {3.565980062656268*^9, 3.5659802109657288`*^9}, {3.5659802575474105`*^9, 3.565980388946441*^9}, { 3.5659804562137594`*^9, 3.5659806055216217`*^9}, {3.5659807011029897`*^9, 3.565980704878196*^9}, {3.56598073534505*^9, 3.5659807355478497`*^9}, { 3.5659807668103046`*^9, 3.565980841503236*^9}, {3.5659809097065563`*^9, 3.565981017939546*^9}, {3.565981082133659*^9, 3.5659810978428864`*^9}, { 3.5659811906786494`*^9, 3.5659813434653177`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"func", "=", RowBox[{ RowBox[{"T", "+", "Vne", "+", "Vee", "+", "Vnn"}], " ", "/.", " ", "c"}]}], ";"}]], "Input", CellChangeTimes->{{3.565955151786976*^9, 3.5659551527073774`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"t", "=", RowBox[{"FindMinimum", "[", RowBox[{"func", ",", RowBox[{"{", RowBox[{"P", ",", "0.32"}], "}"}], ",", RowBox[{"{", RowBox[{"Q", ",", "1.985"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Method", " ", "\[Rule]", " ", "\"\\""}], ",", " ", RowBox[{"MaxIterations", " ", "\[Rule]", " ", "1000"}]}], "}"}]}], "]"}], " "}]], "Input", CellChangeTimes->{{3.5659358688580046`*^9, 3.5659358824300284`*^9}, { 3.56595515980539*^9, 3.565955160834992*^9}, {3.565979036080865*^9, 3.565979071290127*^9}, 3.565979615013482*^9, {3.565980412486882*^9, 3.565980424062103*^9}, 3.5659806610733194`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "27.003190714933073`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.32017612112738686`"}], ",", RowBox[{"Q", "\[Rule]", "2.0446568008725934`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.5659797015156336`*^9, {3.565980070191081*^9, 3.565980214631735*^9}, { 3.5659802616970177`*^9, 3.5659804284301105`*^9}, {3.565980465932576*^9, 3.5659806096088285`*^9}, 3.5659806669857297`*^9, 3.5659807093086042`*^9, 3.565980740149858*^9, {3.5659807711471124`*^9, 3.565980845543643*^9}, { 3.5659809168357687`*^9, 3.56598095181103*^9}, {3.5659809836038857`*^9, 3.565981023446356*^9}, {3.5659810674228325`*^9, 3.565981101805293*^9}, 3.565981147482174*^9, {3.5659811969030604`*^9, 3.5659813489565277`*^9}}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"u", "=", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}], ";"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Vne", " ", "/.", "u"}]], "Input"], Cell[BoxData[ RowBox[{"-", "77.35498726769413`"}]], "Output", CellChangeTimes->{ 3.565979701531234*^9, {3.5659800702066813`*^9, 3.5659802146629353`*^9}, { 3.565980261712618*^9, 3.5659804284457107`*^9}, {3.5659804659481764`*^9, 3.5659806096244287`*^9}, 3.5659806670013294`*^9, 3.565980709324204*^9, 3.5659807401654577`*^9, {3.565980771162712*^9, 3.565980845559243*^9}, { 3.5659809168513684`*^9, 3.5659809518266296`*^9}, {3.565980983619486*^9, 3.5659810234619555`*^9}, {3.5659810674228325`*^9, 3.5659811018208933`*^9}, 3.5659811474977736`*^9, {3.56598119691866*^9, 3.565981348987727*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"Vee", " ", "/.", "u"}], " ", "/.", " ", "c"}]], "Input"], Cell[BoxData["16.11680338503265`"], "Output", CellChangeTimes->{ 3.565979701531234*^9, {3.5659800702222815`*^9, 3.5659802146629353`*^9}, { 3.565980261728218*^9, 3.5659804284613104`*^9}, {3.5659804659637766`*^9, 3.565980609640029*^9}, 3.5659806670169296`*^9, 3.5659807093398037`*^9, 3.565980740181058*^9, {3.5659807711783123`*^9, 3.5659808455748434`*^9}, { 3.565980916866968*^9, 3.56598095184223*^9}, {3.5659809836350856`*^9, 3.5659810234775553`*^9}, {3.5659810674384327`*^9, 3.5659811018364935`*^9}, 3.5659811474977736`*^9, {3.5659811969342604`*^9, 3.565981348987727*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"3", "*", "Vex"}], " ", "/.", " ", "u"}]], "Input"], Cell[BoxData[ RowBox[{"-", "0.006760637218700009`"}]], "Output", CellChangeTimes->{ 3.565979701546834*^9, {3.5659800702378817`*^9, 3.5659802146785355`*^9}, { 3.565980261743818*^9, 3.5659804284769106`*^9}, {3.5659804659793763`*^9, 3.565980609640029*^9}, 3.56598066703253*^9, 3.565980709355404*^9, 3.565980740181058*^9, {3.5659807711939125`*^9, 3.565980845590443*^9}, { 3.5659809168825684`*^9, 3.56598095184223*^9}, {3.5659809836506853`*^9, 3.5659810234931555`*^9}, {3.565981067454033*^9, 3.5659811018520937`*^9}, 3.5659811475133734`*^9, {3.5659811969498606`*^9, 3.5659813490033274`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Vnn", " ", "/.", " ", "u"}]], "Input"], Cell[BoxData["7.436467025586168`"], "Output", CellChangeTimes->{ 3.565979701546834*^9, {3.5659800702378817`*^9, 3.565980214694135*^9}, { 3.5659802617594175`*^9, 3.565980428492511*^9}, {3.5659804659793763`*^9, 3.5659806096556287`*^9}, 3.5659806670481296`*^9, 3.565980709371004*^9, 3.565980740196658*^9, {3.5659807712095127`*^9, 3.565980845590443*^9}, { 3.5659809168981686`*^9, 3.56598095185783*^9}, {3.5659809836662855`*^9, 3.5659810235087557`*^9}, {3.565981067454033*^9, 3.5659811018676934`*^9}, 3.5659811475289736`*^9, {3.5659811969654603`*^9, 3.5659813490189276`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"Vne", "+", "Vee", "+", "Vnn"}], ")"}]}], "/", "T"}], " ", "/.", "u"}], "/.", "c"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"R", "+", "S"}], ")"}], "*", "0.529177"}], " ", "/.", "u"}]}], "Input"], Cell[BoxData["2.007637157793875`"], "Output", CellChangeTimes->{ 3.565979701562434*^9, {3.5659800702534814`*^9, 3.565980214709735*^9}, { 3.5659802617750177`*^9, 3.565980428492511*^9}, {3.565980465994976*^9, 3.565980609671229*^9}, 3.5659806670637293`*^9, 3.5659807093866043`*^9, 3.5659807402122583`*^9, {3.565980771225113*^9, 3.565980845606043*^9}, { 3.5659809169137683`*^9, 3.5659809518734303`*^9}, {3.5659809836662855`*^9, 3.565981023524356*^9}, {3.565981067469633*^9, 3.565981101883293*^9}, 3.5659811475289736`*^9, {3.5659811969810605`*^9, 3.565981349034528*^9}}], Cell[BoxData["1.19064825`"], "Output", CellChangeTimes->{ 3.565979701562434*^9, {3.5659800702534814`*^9, 3.565980214709735*^9}, { 3.5659802617750177`*^9, 3.565980428492511*^9}, {3.565980465994976*^9, 3.565980609671229*^9}, 3.5659806670637293`*^9, 3.5659807093866043`*^9, 3.5659807402122583`*^9, {3.565980771225113*^9, 3.565980845606043*^9}, { 3.5659809169137683`*^9, 3.5659809518734303`*^9}, {3.5659809836662855`*^9, 3.565981023524356*^9}, {3.565981067469633*^9, 3.565981101883293*^9}, 3.5659811475289736`*^9, {3.5659811969810605`*^9, 3.565981349034528*^9}}] }, Open ]], Cell[BoxData[{ RowBox[{ 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"}], "]"}], ",", "\n", "\t\t\t", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"R", ",", "0"}], "}"}], ",", "Q"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "R"}], "/", "2"}], ",", RowBox[{"w3", "*", RowBox[{"R", "/", "2"}]}]}], "}"}], ",", "Q"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "R"}], "/", "2"}], ",", RowBox[{ RowBox[{"-", "w3"}], "*", RowBox[{"R", "/", "2"}]}]}], "}"}], ",", "Q"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"R", "+", "S"}], ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R", "+", "S"}], ")"}]}], "/", "2"}], ",", RowBox[{"w3", "*", RowBox[{ RowBox[{"(", RowBox[{"R", "+", "S"}], ")"}], "/", "2"}]}]}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R", "+", "S"}], ")"}]}], "/", "2"}], ",", RowBox[{ RowBox[{"-", "w3"}], "*", RowBox[{ RowBox[{"(", RowBox[{"R", "+", "S"}], ")"}], "/", "2"}]}]}], "}"}], ",", "0.1"}], "]"}], ",", "\n", "\t\t\t", RowBox[{"{", RowBox[{ RowBox[{"Thickness", "[", "0.012", "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"R", "+", "S"}], ",", "0"}], "}"}]}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R", "+", "S"}], ")"}]}], "/", "2"}], ",", RowBox[{"w3", "*", RowBox[{ RowBox[{"(", RowBox[{"R", "+", "S"}], ")"}], "/", "2"}]}]}], "}"}]}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R", "+", "S"}], ")"}]}], "/", "2"}], ",", RowBox[{ RowBox[{"-", "w3"}], "*", RowBox[{ RowBox[{"(", RowBox[{"R", "+", "S"}], ")"}], "/", "2"}]}]}], "}"}]}], "}"}], "]"}]}], "}"}]}], "}"}], "]"}], "/.", " ", "u"}]}], ";"}], "\n"}], "\n", RowBox[{" "}]}], "Input", CellChangeTimes->{3.5659356141563573`*^9}], Cell[CellGroupData[{ Cell[BoxData[ StyleBox[ RowBox[{"Show", "[", RowBox[{"plot1", ",", " ", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", "->", "Automatic"}], ",", " ", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "3"}], ",", "3.5"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "}"}]}], ",", RowBox[{"Axes", "\[Rule]", "True"}], ",", RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", RowBox[{"Frame", "\[Rule]", "True"}]}], "}"}]}], "]"}], FontSize->12]], "Input", CellChangeTimes->{{3.565955368798957*^9, 3.5659553945390024`*^9}}], Cell[BoxData[ GraphicsBox[{CircleBox[{0, 0}, 0.32017612112738686`], DiskBox[{0, 0}, 0.1], CircleBox[{1.11, 0}, 2.0446568008725934`], CircleBox[{-0.555, 0.961288198200727}, 2.0446568008725934`], CircleBox[{-0.555, -0.961288198200727}, 2.0446568008725934`], DiskBox[{2.25, 0}, 0.1], DiskBox[{-1.125, 1.9485571585149868`}, 0.1], DiskBox[{-1.125, -1.9485571585149868`}, 0.1], {Thickness[0.012], LineBox[{{0, 0}, {2.25, 0}}], LineBox[{{0, 0}, {-1.125, 1.9485571585149868`}}], LineBox[{{0, 0}, {-1.125, -1.9485571585149868`}}]}}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, PlotRange->{{-3, 3.5}, {-3, 3}}]], "Output", CellChangeTimes->{ 3.565979701578034*^9, {3.5659800702690816`*^9, 3.565980214725335*^9}, { 3.565980261790618*^9, 3.5659804285237107`*^9}, {3.5659804660105762`*^9, 3.5659806096868286`*^9}, 3.5659806670793295`*^9, 3.565980709402204*^9, 3.5659807402278585`*^9, {3.565980771240713*^9, 3.565980845621643*^9}, { 3.5659809169293685`*^9, 3.56598095188903*^9}, {3.565980983697486*^9, 3.565981023539956*^9}, {3.5659810674852333`*^9, 3.5659811018988934`*^9}, 3.565981147544574*^9, {3.56598119699666*^9, 3.5659813490501275`*^9}}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"plot2", "=", RowBox[{ RowBox[{"Graphics3D", "[", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.3", "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"R", ",", "0", ",", "0"}], "}"}], ",", "Q"}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "R"}], "/", "2"}], ",", RowBox[{"w3", "*", RowBox[{"R", "/", "2"}]}], ",", "0"}], "}"}], ",", "Q"}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "R"}], "/", "2"}], ",", RowBox[{ RowBox[{"-", "w3"}], "*", RowBox[{"R", "/", "2"}]}], ",", "0"}], "}"}], ",", "Q"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Darker", "[", "Red", "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"R", "+", "S"}], ",", "0", ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R", "+", "S"}], ")"}]}], "/", "2"}], ",", RowBox[{"w3", "*", RowBox[{ RowBox[{"(", RowBox[{"R", "+", "S"}], ")"}], "/", "2"}]}], ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R", "+", "S"}], ")"}]}], "/", "2"}], ",", RowBox[{ RowBox[{"-", "w3"}], "*", RowBox[{ RowBox[{"(", RowBox[{"R", "+", "S"}], ")"}], "/", "2"}]}], ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", "P"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Darker", "[", RowBox[{"Red", ",", "1"}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", "0.1"}], "]"}]}], "}"}]}], "}"}], "]"}], "/.", " ", "u"}]}], ";"}]], "Input", CellChangeTimes->{{3.5659360172610655`*^9, 3.565936046214716*^9}, { 3.5659360762759686`*^9, 3.5659362057093964`*^9}, {3.5659362514954767`*^9, 3.5659362794663258`*^9}, {3.5659363202759976`*^9, 3.5659363211807995`*^9}, {3.565936369478484*^9, 3.565936386045713*^9}, 3.565936446386619*^9, {3.5659365030303183`*^9, 3.5659365147459393`*^9}, { 3.5659365969580836`*^9, 3.565936608127703*^9}, {3.5659366510121784`*^9, 3.565936656347388*^9}, {3.565936732116721*^9, 3.5659367877308187`*^9}, 3.5659368895833974`*^9, {3.5659369318126717`*^9, 3.5659369888307714`*^9}, { 3.5659370409504633`*^9, 3.5659370439456687`*^9}, {3.5659370898721495`*^9, 3.5659371221798058`*^9}, {3.565937209118759*^9, 3.5659372209747796`*^9}, 3.565937260349249*^9, {3.5659495099658947`*^9, 3.5659495328667345`*^9}, { 3.5659497088818436`*^9, 3.5659497198798633`*^9}, 3.565949766898346*^9, 3.5659497991124024`*^9, {3.5659498450388827`*^9, 3.5659498459904847`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ StyleBox["Show", FontSize->12], StyleBox["[", FontSize->12], RowBox[{ StyleBox["plot2", FontSize->12], StyleBox[",", FontSize->12], StyleBox[" ", FontSize->12], RowBox[{ StyleBox["{", FontSize->12], RowBox[{ StyleBox[ RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], FontSize->12], StyleBox[",", FontSize->12], RowBox[{ StyleBox["Boxed", FontSize->12], StyleBox[" ", FontSize->12], StyleBox["\[Rule]", FontSize->12], " ", "False"}]}], StyleBox["}", FontSize->12]}]}], StyleBox["]", FontSize->12]}]], "Input", CellChangeTimes->{{3.5659368348585014`*^9, 3.565936881393383*^9}, { 3.5659369367110806`*^9, 3.565936937194681*^9}, 3.56595541024823*^9}], Cell[BoxData[ Graphics3DBox[ {Opacity[ 0.3], {SphereBox[{1.11, 0, 0}, 2.0446568008725934`], SphereBox[{-0.555, 0.961288198200727, 0}, 2.0446568008725934`], SphereBox[{-0.555, -0.961288198200727, 0}, 2.0446568008725934`]}, {RGBColor[ NCache[ Rational[2, 3], 0.6666666666666666], 0, 0], SphereBox[{2.25, 0, 0}, 0.1], SphereBox[{-1.125, 1.9485571585149868`, 0}, 0.1], SphereBox[{-1.125, -1.9485571585149868`, 0}, 0.1], SphereBox[{0, 0, 0}, 0.32017612112738686`]}, {RGBColor[0, 0, 0], SphereBox[{0, 0, 0}, 0.1]}}, AspectRatio->Automatic, AutomaticImageSize->True, Boxed->False, ImageSize->{365.9423160364236, 371.29580732661987`}, ViewPoint->{-0.09445522150097252, -0.3526151908065938, 3.3640363758948317`}, ViewVertical->{-0.3968070306057606, 0.6570135349795441, 0.9241299912283574}]], "Output", CellChangeTimes->{ 3.5659797015936337`*^9, {3.5659800702846813`*^9, 3.5659802147409353`*^9}, { 3.565980261806218*^9, 3.5659804285393105`*^9}, {3.5659804660417767`*^9, 3.565980609702429*^9}, 3.5659806670949297`*^9, 3.565980709417804*^9, 3.565980740259058*^9, {3.565980771256313*^9, 3.5659808456528435`*^9}, { 3.5659809169449687`*^9, 3.5659809519046297`*^9}, {3.565980983713086*^9, 3.565981023555556*^9}, {3.565981067500833*^9, 3.5659811019144936`*^9}, 3.5659811475601735`*^9, {3.5659811970278606`*^9, 3.5659813490657277`*^9}}] }, Open ]] }, WindowSize->{766, 740}, WindowMargins->{{Automatic, 492}, {Automatic, -23}}, PrintingCopies->1, PrintingPageRange->{32000, 32000}, PrintingOptions->{"Magnification"->1., "PaperOrientation"->"Portrait", "PaperSize"->{595.2, 841.8}}, PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}}, CellLabelAutoDelete->True, Magnification->0.8999999761581421, 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, 1216, 29, 130, "Input"], Cell[1776, 51, 6335, 200, 650, "Input"], Cell[8114, 253, 237, 10, 29, "Input"], Cell[8354, 265, 2913, 65, 130, "Input"], Cell[11270, 332, 229, 6, 29, "Input"], Cell[CellGroupData[{ Cell[11524, 342, 692, 16, 29, "Input"], Cell[12219, 360, 811, 16, 28, "Output"] }, Open ]], Cell[13045, 379, 121, 4, 29, "Input"], Cell[CellGroupData[{ Cell[13191, 387, 57, 1, 29, "Input"], Cell[13251, 390, 612, 9, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[13900, 404, 92, 2, 29, "Input"], Cell[13995, 408, 597, 8, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[14629, 421, 85, 2, 29, "Input"], Cell[14717, 425, 612, 9, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[15366, 439, 62, 1, 29, "Input"], Cell[15431, 442, 595, 8, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[16063, 455, 336, 12, 50, "Input"], Cell[16402, 469, 590, 8, 28, "Output"], Cell[16995, 479, 583, 8, 28, "Output"] }, Open ]], Cell[17593, 490, 3936, 118, 190, "Input"], Cell[CellGroupData[{ Cell[21554, 612, 723, 20, 46, "Input"], Cell[22280, 634, 1227, 22, 320, "Output"] }, Open ]], Cell[23522, 659, 3746, 95, 130, "Input"], Cell[CellGroupData[{ Cell[27293, 758, 820, 34, 29, "Input"], Cell[28116, 794, 1414, 28, 347, "Output"] }, Open ]] } ] *) (* End of internal cache information *)