(* 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[ 88995, 1603] NotebookOptionsPosition[ 86709, 1531] NotebookOutlinePosition[ 87306, 1553] CellTagsIndexPosition[ 87263, 1550] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell["\<\ (*NH2COH Formamide, with \"orthogonal\" \[Sigma],\[Pi] LPs, 13.01.2013, not \ final optimization: slight overlaps *)\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{ 3.5668390630515456`*^9, {3.5668394752510695`*^9, 3.5668394888386936`*^9}, { 3.566840005402401*^9, 3.5668400293172426`*^9}, {3.566840906678384*^9, 3.566840923027212*^9}, {3.5668417775723133`*^9, 3.5668417895999346`*^9}, { 3.566841995894697*^9, 3.5668420196535387`*^9}, {3.566842057031204*^9, 3.5668420654552193`*^9}, {3.5668422519691467`*^9, 3.5668422563371544`*^9}, 3.5668422887072115`*^9, {3.566842337504097*^9, 3.566842360326937*^9}, { 3.5668424117758274`*^9, 3.5668424218690453`*^9}, {3.566842627885008*^9, 3.566842684060706*^9}, {3.566842771030859*^9, 3.5668427769432697`*^9}, { 3.5668428297025623`*^9, 3.5668428906206694`*^9}, {3.5668429241763277`*^9, 3.566842930135538*^9}, {3.5668871695283947`*^9, 3.5668871699963956`*^9}, { 3.5668872245808916`*^9, 3.566887225407693*^9}, {3.5668875348030367`*^9, 3.566887535037037*^9}, {3.566887597327946*^9, 3.566887637981618*^9}, 3.5668876998201265`*^9, {3.5668877336253853`*^9, 3.5668877352165885`*^9}, { 3.566887814495928*^9, 3.5668878152135286`*^9}, {3.566887857926404*^9, 3.566887961448186*^9}, 3.5668880086070685`*^9, {3.5668880951092205`*^9, 3.5668881570569296`*^9}, {3.5668882853983545`*^9, 3.5668884070317683`*^9}, {3.566888509258748*^9, 3.5668885097423487`*^9}, { 3.566888582719277*^9, 3.566888583795679*^9}, {3.566888668176227*^9, 3.566888671935834*^9}, {3.5668887042746906`*^9, 3.566888705070292*^9}, { 3.566888771822809*^9, 3.566888808857274*^9}, {3.5668889463559155`*^9, 3.566888956839134*^9}, {3.5669685950612297`*^9, 3.5669685991952367`*^9}, { 3.566968722700654*^9, 3.5669687698127365`*^9}, {3.5669688330864477`*^9, 3.566968848015674*^9}, {3.5669689606010714`*^9, 3.566968973471094*^9}, { 3.5669881655485907`*^9, 3.5669881659853916`*^9}, {3.566996986024498*^9, 3.5669969910321074`*^9}, {3.5669973015818524`*^9, 3.5669973626247597`*^9}, {3.566997996329073*^9, 3.566998005049488*^9}, { 3.56705755136038*^9, 3.5670575543243856`*^9}, {3.5670590414593973`*^9, 3.5670590456714053`*^9}, {3.567061276787324*^9, 3.567061372633892*^9}, { 3.567061756144966*^9, 3.5670617569249673`*^9}, {3.567061904329626*^9, 3.567061915546046*^9}, {3.567062309945938*^9, 3.5670623507556105`*^9}, { 3.567062381019663*^9, 3.5670625124810944`*^9}, {3.5670638527106485`*^9, 3.5670638629598665`*^9}, 3.5670687979573345`*^9, {3.567078692613513*^9, 3.56707870788594*^9}}], Cell["\<\ Clear[ZC,ZN,ZO,sig1,sig2,sig3,sig4,sig5,sig6,sig7,sig8,sig9,k1,k2,k3,k4k5,k6,\ k7,k8,k9,k10,k11,nc,w1,w2,w3,pnh,pch,vee,vne,vnn, xn,yn,zn,xc,yc,zc,c1,c2,c3,cs1,ss1,cs2,ss2,cs3,ss3,pi,i,j,d,t]; Clear[Cr,Nr,Oa,NH,CH,CO,CN,OL,NM,OC,OM,w1,w2,w3,pch,pnh]; ZC=6.0; ZN=7.0; ZO=8.0; sig1=0.3; sig2=0.3; sig3=0.3; sig4=0.3; sig5=0.3; \ sig6=0.3; sig7=0.3; sig8=0.3; sig9=0.3; nc=15; nn=6; pi=0.0; k1=1.02; k2=1.02; k3=1.02; k4=1.2; k5=1.25; k6=1.56; k7=1.26; k8=1.4; k9=1.1; \ k10=1.1; k11=1.0; bohr=0.529177; rad=57.29578; (* 16 Variables: Cr,Nr,Or,NH,CH,CN,CO,OL,NM,OC,OM,pnh,pch,w1,w2,w3 *)\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{ 3.5668390630515456`*^9, {3.5668394752510695`*^9, 3.5668394888386936`*^9}, { 3.566840005402401*^9, 3.5668400293172426`*^9}, {3.566840906678384*^9, 3.566840923027212*^9}, {3.5668417775723133`*^9, 3.5668417895999346`*^9}, { 3.566841995894697*^9, 3.5668420196535387`*^9}, {3.566842057031204*^9, 3.5668420654552193`*^9}, {3.5668422519691467`*^9, 3.5668422563371544`*^9}, 3.5668422887072115`*^9, {3.566842337504097*^9, 3.566842360326937*^9}, { 3.5668424117758274`*^9, 3.5668424218690453`*^9}, {3.566842627885008*^9, 3.566842684060706*^9}, {3.566842771030859*^9, 3.5668427769432697`*^9}, { 3.5668428297025623`*^9, 3.5668428906206694`*^9}, {3.5668429241763277`*^9, 3.566842930135538*^9}, {3.5668871695283947`*^9, 3.5668871699963956`*^9}, { 3.5668872245808916`*^9, 3.566887225407693*^9}, {3.5668875348030367`*^9, 3.566887535037037*^9}, {3.566887597327946*^9, 3.566887637981618*^9}, 3.5668876998201265`*^9, {3.5668877336253853`*^9, 3.5668877352165885`*^9}, { 3.566887814495928*^9, 3.5668878152135286`*^9}, {3.566887857926404*^9, 3.566887961448186*^9}, 3.5668880086070685`*^9, {3.5668880951092205`*^9, 3.5668881570569296`*^9}, {3.5668882853983545`*^9, 3.5668884070317683`*^9}, {3.566888509258748*^9, 3.5668885097423487`*^9}, { 3.566888582719277*^9, 3.566888583795679*^9}, {3.566888668176227*^9, 3.566888671935834*^9}, {3.5668887042746906`*^9, 3.566888705070292*^9}, { 3.566888771822809*^9, 3.566888808857274*^9}, {3.5668889463559155`*^9, 3.566888956839134*^9}, {3.5669685950612297`*^9, 3.5669685991952367`*^9}, { 3.566968722700654*^9, 3.5669687698127365`*^9}, {3.5669688330864477`*^9, 3.566968848015674*^9}, {3.5669689606010714`*^9, 3.566968973471094*^9}, { 3.5669881655485907`*^9, 3.5669881659853916`*^9}, {3.566996986024498*^9, 3.5669969910321074`*^9}, {3.5669973015818524`*^9, 3.5669973626247597`*^9}, {3.566997996329073*^9, 3.566998005049488*^9}, { 3.56705755136038*^9, 3.5670575543243856`*^9}, {3.5670590414593973`*^9, 3.5670590456714053`*^9}, {3.567061276787324*^9, 3.567061372633892*^9}, { 3.567061756144966*^9, 3.5670617569249673`*^9}, {3.567061904329626*^9, 3.567061915546046*^9}, {3.567062309945938*^9, 3.5670623507556105`*^9}, { 3.567062381019663*^9, 3.5670625124810944`*^9}, {3.5670638527106485`*^9, 3.5670638629598665`*^9}, 3.5670687979573345`*^9, {3.567069322009055*^9, 3.5670693806495576`*^9}, {3.567069419290826*^9, 3.567069451957283*^9}, { 3.5670694881649466`*^9, 3.5670695381162343`*^9}, {3.56708248268017*^9, 3.567082487516179*^9}, {3.567085712191843*^9, 3.56708586995492*^9}, { 3.5670863390009437`*^9, 3.5670863409197474`*^9}, {3.5670863724006023`*^9, 3.5670863736018047`*^9}, {3.567086415363078*^9, 3.56708641629908*^9}, { 3.5670864537547455`*^9, 3.567086586245778*^9}, {3.56708706580622*^9, 3.567087066742222*^9}, {3.5670871956452484`*^9, 3.567087197376851*^9}, { 3.567087313846656*^9, 3.5670873147670574`*^9}, {3.567087346825114*^9, 3.5670873478703156`*^9}, {3.5670873954503994`*^9, 3.567087396292801*^9}, { 3.5670874292868586`*^9, 3.5670874532017007`*^9}, {3.567087493090971*^9, 3.5670874935745716`*^9}, {3.5670979031336823`*^9, 3.5670979250517206`*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{"Cr", "=", "0.26"}], ";", " ", RowBox[{"Nr", "=", "0.22"}], ";", " ", RowBox[{"Oa", "=", "0.186"}], ";", " ", RowBox[{"NH", "=", "1.2"}], ";", " ", RowBox[{"CH", "=", "1.25"}], ";", " ", RowBox[{"CO", "=", "0.928"}], ";", " ", RowBox[{"CN", "=", "1.04"}], ";", " ", RowBox[{"OL", "=", "1.0"}], ";", "\[IndentingNewLine]", RowBox[{"NM", "=", "1.1"}], ";", " ", RowBox[{"OC", "=", "1.4"}], ";", " ", RowBox[{"OM", "=", "1.05"}], ";", " ", RowBox[{"w1", "=", "1.9"}], ";", " ", RowBox[{"w2", "=", "1.0"}], ";", " ", RowBox[{"w3", "=", "1.0"}], ";", " ", RowBox[{"pch", "=", "0.58"}], ";", RowBox[{"pnh", "=", "0.5"}], ";"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.5670689879656677`*^9, 3.5670691959452333`*^9}, 3.5670698299303465`*^9, {3.567077287748646*^9, 3.5670773058290777`*^9}, { 3.567077396527637*^9, 3.567077402346447*^9}, {3.5670803249623804`*^9, 3.5670804001545124`*^9}, 3.5670804667198296`*^9, 3.5670805071239004`*^9, { 3.5670810452624454`*^9, 3.567081083794513*^9}, {3.5670811174593725`*^9, 3.5670811463974233`*^9}, {3.5670811875190954`*^9, 3.5670812332895756`*^9}, {3.5670812684832373`*^9, 3.567081269356839*^9}, { 3.5670852872782965`*^9, 3.567085293970708*^9}}], Cell["\<\ (* C He shell *) Ekin = 2.25*k1/Cr^2; vee=3.0*sig1/Cr; vne=-3.0*ZC/Cr; (* N He shell *) Ekin = Ekin+2.25*k2/Nr^2; vee=vee+3.0*sig2/Nr; vne=vne-3.0*ZN/Nr; (* O He shell *) Ekin = Ekin+2.25*k3/Oa^2; vee=vee+3.0*sig3/Oa; vne=vne-3.0*ZO/Oa;\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5668390630515456`*^9, 3.5668390682151546`*^9}, { 3.566971938334302*^9, 3.566971944714713*^9}, {3.566972152507078*^9, 3.5669721651899004`*^9}, {3.5669970287685733`*^9, 3.5669972828462195`*^9}, { 3.566997377709986*^9, 3.5669974183012576`*^9}, {3.5669979768290386`*^9, 3.566997987171857*^9}, {3.566998792414071*^9, 3.5669988056272945`*^9}, { 3.567063874659887*^9, 3.567063885189905*^9}, {3.5670698395087633`*^9, 3.567069848260379*^9}}], Cell["\<\ (* bonding and lone pairs *) Ekin = Ekin + \ 2.25*(2*k4/NH^2+k5/CN^2+k6/CO^2+k7/CH^2+k8/NM^2+k9/OL^2+k10/OM^2+k11/OC^2); vee=vee+3.0*(2*sig4/NH+sig5/CN+sig6/CO+sig7/CH+sig9/OL);\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5668390630515456`*^9, 3.566839078760773*^9}, { 3.566839318626795*^9, 3.5668393238528037`*^9}, {3.566839506201524*^9, 3.5668395108503323`*^9}, {3.566839541207986*^9, 3.566839568664034*^9}, { 3.566839684650237*^9, 3.56683970887708*^9}, {3.5668397799352045`*^9, 3.5668398201988754`*^9}, {3.566839877872177*^9, 3.5668398909606*^9}, { 3.5668404106911125`*^9, 3.566840417539525*^9}, {3.5668404815464373`*^9, 3.566840691460406*^9}, {3.566840732519678*^9, 3.5668408608923035`*^9}, { 3.566840930000425*^9, 3.5668409412168446`*^9}, {3.566841374686006*^9, 3.566841375388007*^9}, {3.5668421835006266`*^9, 3.5668421845770283`*^9}, { 3.566887398037596*^9, 3.566887402078003*^9}, {3.5668874322328568`*^9, 3.5668874455864797`*^9}, {3.566887493166563*^9, 3.5668874939309645`*^9}, { 3.5668882488786907`*^9, 3.5668882511094947`*^9}, {3.56688847075788*^9, 3.566888479322295*^9}, {3.5668885903476906`*^9, 3.5668885954956994`*^9}, { 3.566971858103361*^9, 3.566971918927868*^9}, {3.566987354128765*^9, 3.5669873624903803`*^9}, {3.566987488164201*^9, 3.566987535619484*^9}, { 3.566987795843541*^9, 3.5669878015843515`*^9}, {3.566987856168847*^9, 3.56698786325126*^9}, {3.566988264733565*^9, 3.566988269819174*^9}, { 3.566988366133743*^9, 3.5669883719525533`*^9}, {3.5669885431160536`*^9, 3.566988549309265*^9}, {3.566991946182942*^9, 3.5669919483981457`*^9}, 3.5669919866182127`*^9, {3.5669920226542764`*^9, 3.5669920505783253`*^9}, { 3.566992103743219*^9, 3.5669921079396257`*^9}, 3.5669921962201815`*^9, { 3.5669922303998413`*^9, 3.5669922340658474`*^9}, {3.566992291177548*^9, 3.5669922926751504`*^9}, {3.566992333890423*^9, 3.5669923348264246`*^9}, 3.566992409675356*^9, {3.566997472043352*^9, 3.566997587171554*^9}, { 3.566997674297707*^9, 3.5669977903463106`*^9}, {3.566997868580448*^9, 3.5669978780652647`*^9}, {3.566998052270771*^9, 3.566998071786405*^9}, { 3.5669981665409718`*^9, 3.5669981730305834`*^9}, {3.56699824522751*^9, 3.566998326394452*^9}, {3.5669983623837156`*^9, 3.5669984017425847`*^9}, { 3.566998470788306*^9, 3.5669984830811276`*^9}, {3.566998534405218*^9, 3.5669985995821323`*^9}, {3.5669988153461113`*^9, 3.566998936995125*^9}, 3.566998967134378*^9, 3.5669990224832754`*^9, {3.5669990973790064`*^9, 3.56699912179305*^9}, {3.5669993047345705`*^9, 3.5669993241410046`*^9}, { 3.566999416274767*^9, 3.566999426742385*^9}, {3.566999497301309*^9, 3.5669995032137194`*^9}, {3.566999549733001*^9, 3.566999551433404*^9}, { 3.5670567651502*^9, 3.567056829110312*^9}, {3.5670569077344503`*^9, 3.5670570287126627`*^9}, {3.56705707303234*^9, 3.5670571842137356`*^9}, { 3.5670572582982655`*^9, 3.567057292103525*^9}, {3.56705732891959*^9, 3.5670574906294737`*^9}, {3.5670578301328697`*^9, 3.56705795512029*^9}, { 3.5670580197668033`*^9, 3.567058030437222*^9}, {3.5670580810281105`*^9, 3.5670581611654515`*^9}, {3.5670581932547083`*^9, 3.567058240616391*^9}, { 3.5670582801312604`*^9, 3.5670582990540934`*^9}, {3.567058343560972*^9, 3.5670584018894744`*^9}, {3.56705845049916*^9, 3.5670585083128614`*^9}, { 3.567058578731385*^9, 3.567058593645011*^9}, {3.5670586348134837`*^9, 3.5670586577143235`*^9}, {3.567058689007979*^9, 3.5670587629209085`*^9}, { 3.567058832777831*^9, 3.5670588394702425`*^9}, {3.567058955628047*^9, 3.5670589616496572`*^9}, {3.5670592284569263`*^9, 3.567059354442747*^9}, { 3.5670593898548098`*^9, 3.567059560456709*^9}, {3.5670596222796183`*^9, 3.5670596362884426`*^9}, {3.5670596726677065`*^9, 3.567059775347087*^9}, { 3.567059838449198*^9, 3.5670598453756094`*^9}, {3.5670598764976645`*^9, 3.5670598805380716`*^9}, 3.567059911254525*^9, {3.5670599744502363`*^9, 3.567060007818695*^9}, {3.567060059127185*^9, 3.5670600818564253`*^9}, { 3.5670601143824825`*^9, 3.567060314593234*^9}, {3.5670607181347427`*^9, 3.5670607311607656`*^9}, {3.567060785573661*^9, 3.567061232998047*^9}, { 3.5670626282800975`*^9, 3.56706267537658*^9}, 3.567062795371991*^9, { 3.567069851801585*^9, 3.567069879460434*^9}, {3.567070767429593*^9, 3.5670707823588195`*^9}, {3.5670805354535503`*^9, 3.5670806312689185`*^9}, {3.567083393113369*^9, 3.5670834586022844`*^9}, { 3.567083501299559*^9, 3.5670835804696984`*^9}, {3.5670836156477604`*^9, 3.5670836215757704`*^9}, {3.567083656504232*^9, 3.5670836570034328`*^9}, 3.5670837376399746`*^9}], Cell["\<\ (* 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 *) (* C N O Hn Hn Hc B B L2 L1 L1 Bp Bp L1 L1 *) as={C,N,O,H,H,H}; (* atom string *) rr={Cr,Nr,Oa,NH,NH,CH,CN,CO,OL,NM,NM,OC,OC,OM,OM}; oc={-2,-2,-2,-2,-2,-2,-2,-2,-2,-1,-1,-1,-1,-1,-1}; ch={ZC,ZN,ZO, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0}; csw1=Cos[w1/2]; ssw1=Sin[w1/2]; csw2=Cos[w2]; ssw2=Sin[w2]; csw3=Cos[w3]; ssw3=Sin[w3];\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5668390630515456`*^9, 3.566839078760773*^9}, { 3.566839318626795*^9, 3.5668393238528037`*^9}, {3.566839506201524*^9, 3.5668395108503323`*^9}, {3.566839541207986*^9, 3.566839568664034*^9}, { 3.566839684650237*^9, 3.56683970887708*^9}, {3.5668397799352045`*^9, 3.5668398201988754`*^9}, {3.566839877872177*^9, 3.5668398909606*^9}, { 3.5668404106911125`*^9, 3.566840417539525*^9}, {3.5668404815464373`*^9, 3.566840691460406*^9}, {3.566840732519678*^9, 3.5668408608923035`*^9}, { 3.566840930000425*^9, 3.5668409412168446`*^9}, {3.566841374686006*^9, 3.566841375388007*^9}, {3.5668421835006266`*^9, 3.5668421845770283`*^9}, { 3.566887398037596*^9, 3.566887402078003*^9}, {3.5668874322328568`*^9, 3.5668874455864797`*^9}, {3.566887493166563*^9, 3.5668874939309645`*^9}, { 3.5668882488786907`*^9, 3.5668882511094947`*^9}, {3.56688847075788*^9, 3.566888479322295*^9}, {3.5668885903476906`*^9, 3.5668885954956994`*^9}, { 3.566971858103361*^9, 3.566971918927868*^9}, {3.566987354128765*^9, 3.5669873624903803`*^9}, {3.566987488164201*^9, 3.566987535619484*^9}, { 3.566987795843541*^9, 3.5669878015843515`*^9}, {3.566987856168847*^9, 3.56698786325126*^9}, {3.566988264733565*^9, 3.566988269819174*^9}, { 3.566988366133743*^9, 3.5669883719525533`*^9}, {3.5669885431160536`*^9, 3.566988549309265*^9}, {3.566991946182942*^9, 3.5669919483981457`*^9}, 3.5669919866182127`*^9, {3.5669920226542764`*^9, 3.5669920505783253`*^9}, { 3.566992103743219*^9, 3.5669921079396257`*^9}, 3.5669921962201815`*^9, { 3.5669922303998413`*^9, 3.5669922340658474`*^9}, {3.566992291177548*^9, 3.5669922926751504`*^9}, {3.566992333890423*^9, 3.5669923348264246`*^9}, 3.566992409675356*^9, {3.566997472043352*^9, 3.566997587171554*^9}, { 3.566997674297707*^9, 3.5669977903463106`*^9}, {3.566997868580448*^9, 3.5669978780652647`*^9}, {3.566998052270771*^9, 3.566998071786405*^9}, { 3.5669981665409718`*^9, 3.5669981730305834`*^9}, {3.56699824522751*^9, 3.566998326394452*^9}, {3.5669983623837156`*^9, 3.5669984017425847`*^9}, { 3.566998470788306*^9, 3.5669984830811276`*^9}, {3.566998534405218*^9, 3.5669985995821323`*^9}, {3.5669988153461113`*^9, 3.566998936995125*^9}, 3.566998967134378*^9, 3.5669990224832754`*^9, {3.5669990973790064`*^9, 3.56699912179305*^9}, {3.5669993047345705`*^9, 3.5669993241410046`*^9}, { 3.566999416274767*^9, 3.566999426742385*^9}, {3.566999497301309*^9, 3.5669995032137194`*^9}, {3.566999549733001*^9, 3.566999551433404*^9}, { 3.5670567651502*^9, 3.567056829110312*^9}, {3.5670569077344503`*^9, 3.5670570287126627`*^9}, {3.56705707303234*^9, 3.5670571842137356`*^9}, { 3.5670572582982655`*^9, 3.567057292103525*^9}, {3.56705732891959*^9, 3.5670574906294737`*^9}, {3.5670578301328697`*^9, 3.56705795512029*^9}, { 3.5670580197668033`*^9, 3.567058030437222*^9}, {3.5670580810281105`*^9, 3.5670581611654515`*^9}, {3.5670581932547083`*^9, 3.567058240616391*^9}, { 3.5670582801312604`*^9, 3.5670582990540934`*^9}, {3.567058343560972*^9, 3.5670584018894744`*^9}, {3.56705845049916*^9, 3.5670585083128614`*^9}, { 3.567058578731385*^9, 3.567058593645011*^9}, {3.5670586348134837`*^9, 3.5670586577143235`*^9}, {3.567058689007979*^9, 3.5670587629209085`*^9}, { 3.567058832777831*^9, 3.5670588394702425`*^9}, {3.567058955628047*^9, 3.5670589616496572`*^9}, {3.5670592284569263`*^9, 3.567059354442747*^9}, { 3.5670593898548098`*^9, 3.567059560456709*^9}, {3.5670596222796183`*^9, 3.5670596362884426`*^9}, {3.5670596726677065`*^9, 3.567059775347087*^9}, { 3.567059838449198*^9, 3.5670598453756094`*^9}, {3.5670598764976645`*^9, 3.5670598805380716`*^9}, 3.567059911254525*^9, {3.5670599744502363`*^9, 3.567060007818695*^9}, {3.567060059127185*^9, 3.5670600818564253`*^9}, { 3.5670601143824825`*^9, 3.567060314593234*^9}, {3.5670607181347427`*^9, 3.5670607311607656`*^9}, {3.567060785573661*^9, 3.567061232998047*^9}, { 3.5670626282800975`*^9, 3.56706267537658*^9}, 3.567062795371991*^9, { 3.567069851801585*^9, 3.567069879460434*^9}, {3.567070767429593*^9, 3.5670707823588195`*^9}, {3.5670805354535503`*^9, 3.5670806312689185`*^9}, {3.567083393113369*^9, 3.5670834586022844`*^9}, { 3.567083501299559*^9, 3.5670835804696984`*^9}, {3.5670836156477604`*^9, 3.5670836215757704`*^9}, {3.567083656504232*^9, 3.5670836570034328`*^9}, 3.5670837376399746`*^9, {3.567228361875599*^9, 3.5672284052748756`*^9}}], Cell["\<\ dcnx=Cr+Nr+2*CN; dcox=(Cr+Oa+2*CO)*csw2; dnhx=dcnx+(Nr+NH)*csw1; \ dchx=(Cr+CH)*csw3; dccnx=Cr+CN; pnhx=dcnx+(Nr+NH+pnh)*csw1; pchx=(Cr+CH+pch)*csw3; dcoy=(Cr+Oa+2*CO)*ssw2; dnhy=(Nr+NH)*ssw1; dchy=(Cr+CH)*ssw3; pnhy=(Nr+NH+pnh)*ssw1; pchy=(Cr+CH+pch)*ssw3; dccox=(Cr+CO)*csw2; dccoy=(Cr+CO)*ssw2; dolx=(Cr+2*CO+2*Oa+OL)*csw2; doly=(Cr+2*CO+2*Oa+OL)*ssw2; dom1x=dcox-(Oa+OM)*ssw2; dom2x=dcox+(Oa+OM)*ssw2; dom1y=dcoy+(Oa+OM)*csw2; dom2y=dcoy-(Oa+OM)*csw2;\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5668390630515456`*^9, 3.566839078760773*^9}, { 3.566839318626795*^9, 3.5668393238528037`*^9}, {3.566839506201524*^9, 3.5668395108503323`*^9}, {3.566839541207986*^9, 3.566839568664034*^9}, { 3.566839684650237*^9, 3.56683970887708*^9}, {3.5668397799352045`*^9, 3.5668398201988754`*^9}, {3.566839877872177*^9, 3.5668398909606*^9}, { 3.5668404106911125`*^9, 3.566840417539525*^9}, {3.5668404815464373`*^9, 3.566840691460406*^9}, {3.566840732519678*^9, 3.5668408608923035`*^9}, { 3.566840930000425*^9, 3.5668409412168446`*^9}, {3.566841374686006*^9, 3.566841375388007*^9}, {3.5668421835006266`*^9, 3.5668421845770283`*^9}, { 3.566887398037596*^9, 3.566887402078003*^9}, {3.5668874322328568`*^9, 3.5668874455864797`*^9}, {3.566887493166563*^9, 3.5668874939309645`*^9}, { 3.5668882488786907`*^9, 3.5668882511094947`*^9}, {3.56688847075788*^9, 3.566888479322295*^9}, {3.5668885903476906`*^9, 3.5668885954956994`*^9}, { 3.566971858103361*^9, 3.566971918927868*^9}, {3.566987354128765*^9, 3.5669873624903803`*^9}, {3.566987488164201*^9, 3.566987535619484*^9}, { 3.566987795843541*^9, 3.5669878015843515`*^9}, {3.566987856168847*^9, 3.56698786325126*^9}, {3.566988264733565*^9, 3.566988269819174*^9}, { 3.566988366133743*^9, 3.5669883719525533`*^9}, {3.5669885431160536`*^9, 3.566988549309265*^9}, {3.566991946182942*^9, 3.5669919483981457`*^9}, 3.5669919866182127`*^9, {3.5669920226542764`*^9, 3.5669920505783253`*^9}, { 3.566992103743219*^9, 3.5669921079396257`*^9}, 3.5669921962201815`*^9, { 3.5669922303998413`*^9, 3.5669922340658474`*^9}, {3.566992291177548*^9, 3.5669922926751504`*^9}, {3.566992333890423*^9, 3.5669923348264246`*^9}, 3.566992409675356*^9, {3.566997472043352*^9, 3.566997587171554*^9}, { 3.566997674297707*^9, 3.5669977903463106`*^9}, {3.566997868580448*^9, 3.5669978780652647`*^9}, {3.566998052270771*^9, 3.566998071786405*^9}, { 3.5669981665409718`*^9, 3.5669981730305834`*^9}, {3.56699824522751*^9, 3.566998326394452*^9}, {3.5669983623837156`*^9, 3.5669984017425847`*^9}, { 3.566998470788306*^9, 3.5669984830811276`*^9}, {3.566998534405218*^9, 3.5669985995821323`*^9}, {3.5669988153461113`*^9, 3.566998936995125*^9}, 3.566998967134378*^9, 3.5669990224832754`*^9, {3.5669990973790064`*^9, 3.56699912179305*^9}, {3.5669993047345705`*^9, 3.5669993241410046`*^9}, { 3.566999416274767*^9, 3.566999426742385*^9}, {3.566999497301309*^9, 3.5669995032137194`*^9}, {3.566999549733001*^9, 3.566999551433404*^9}, { 3.5670567651502*^9, 3.567056829110312*^9}, {3.5670569077344503`*^9, 3.5670570287126627`*^9}, {3.56705707303234*^9, 3.5670571842137356`*^9}, { 3.5670572582982655`*^9, 3.567057292103525*^9}, {3.56705732891959*^9, 3.5670574906294737`*^9}, {3.5670578301328697`*^9, 3.56705795512029*^9}, { 3.5670580197668033`*^9, 3.567058030437222*^9}, {3.5670580810281105`*^9, 3.5670581611654515`*^9}, {3.5670581932547083`*^9, 3.567058240616391*^9}, { 3.5670582801312604`*^9, 3.5670582990540934`*^9}, {3.567058343560972*^9, 3.5670584018894744`*^9}, {3.56705845049916*^9, 3.5670585083128614`*^9}, { 3.567058578731385*^9, 3.567058593645011*^9}, {3.5670586348134837`*^9, 3.5670586577143235`*^9}, {3.567058689007979*^9, 3.5670587629209085`*^9}, { 3.567058832777831*^9, 3.5670588394702425`*^9}, {3.567058955628047*^9, 3.5670589616496572`*^9}, {3.5670592284569263`*^9, 3.567059354442747*^9}, { 3.5670593898548098`*^9, 3.567059560456709*^9}, {3.5670596222796183`*^9, 3.5670596362884426`*^9}, {3.5670596726677065`*^9, 3.567059775347087*^9}, { 3.567059838449198*^9, 3.5670598453756094`*^9}, {3.5670598764976645`*^9, 3.5670598805380716`*^9}, 3.567059911254525*^9, {3.5670599744502363`*^9, 3.567060007818695*^9}, {3.567060059127185*^9, 3.5670600818564253`*^9}, { 3.5670601143824825`*^9, 3.567060314593234*^9}, {3.5670607181347427`*^9, 3.5670607311607656`*^9}, {3.567060785573661*^9, 3.567061232998047*^9}, { 3.5670626282800975`*^9, 3.56706267537658*^9}, 3.567062795371991*^9, { 3.567069851801585*^9, 3.567069879460434*^9}, {3.567070767429593*^9, 3.5670707823588195`*^9}, {3.5670805354535503`*^9, 3.5670806312689185`*^9}, {3.567083393113369*^9, 3.5670834586022844`*^9}, { 3.567083501299559*^9, 3.5670835804696984`*^9}, {3.5670836156477604`*^9, 3.5670836215757704`*^9}, {3.567083656504232*^9, 3.5670836570034328`*^9}, { 3.5670837376399746`*^9, 3.56708374668799*^9}}], Cell["\<\ (* 6 nuclei: coord *) xn={0,-dcnx,dcox,-pnhx,-pnhx,pchx}; yn={0,0,dcoy,pnhy,-pnhy,-pchy}; zn={0,0,0,0,0,0};\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5668390630515456`*^9, 3.566839078760773*^9}, { 3.566839318626795*^9, 3.5668393238528037`*^9}, {3.566839506201524*^9, 3.5668395108503323`*^9}, {3.566839541207986*^9, 3.566839568664034*^9}, { 3.566839684650237*^9, 3.56683970887708*^9}, {3.5668397799352045`*^9, 3.5668398201988754`*^9}, {3.566839877872177*^9, 3.5668398909606*^9}, { 3.5668404106911125`*^9, 3.566840417539525*^9}, {3.5668404815464373`*^9, 3.566840691460406*^9}, {3.566840732519678*^9, 3.5668408608923035`*^9}, { 3.566840930000425*^9, 3.5668409412168446`*^9}, {3.566841374686006*^9, 3.566841375388007*^9}, {3.5668421835006266`*^9, 3.5668421845770283`*^9}, { 3.566887398037596*^9, 3.566887402078003*^9}, {3.5668874322328568`*^9, 3.5668874455864797`*^9}, {3.566887493166563*^9, 3.5668874939309645`*^9}, { 3.5668882488786907`*^9, 3.5668882511094947`*^9}, {3.56688847075788*^9, 3.566888479322295*^9}, {3.5668885903476906`*^9, 3.5668885954956994`*^9}, { 3.566971858103361*^9, 3.566971918927868*^9}, {3.566987354128765*^9, 3.5669873624903803`*^9}, {3.566987488164201*^9, 3.566987535619484*^9}, { 3.566987795843541*^9, 3.5669878015843515`*^9}, {3.566987856168847*^9, 3.56698786325126*^9}, {3.566988264733565*^9, 3.566988269819174*^9}, { 3.566988366133743*^9, 3.5669883719525533`*^9}, {3.5669885431160536`*^9, 3.566988549309265*^9}, {3.566991946182942*^9, 3.5669919483981457`*^9}, 3.5669919866182127`*^9, {3.5669920226542764`*^9, 3.5669920505783253`*^9}, { 3.566992103743219*^9, 3.5669921079396257`*^9}, 3.5669921962201815`*^9, { 3.5669922303998413`*^9, 3.5669922340658474`*^9}, {3.566992291177548*^9, 3.5669922926751504`*^9}, {3.566992333890423*^9, 3.5669923348264246`*^9}, 3.566992409675356*^9, {3.566997472043352*^9, 3.566997587171554*^9}, { 3.566997674297707*^9, 3.5669977903463106`*^9}, {3.566997868580448*^9, 3.5669978780652647`*^9}, {3.566998052270771*^9, 3.566998071786405*^9}, { 3.5669981665409718`*^9, 3.5669981730305834`*^9}, {3.56699824522751*^9, 3.566998326394452*^9}, {3.5669983623837156`*^9, 3.5669984017425847`*^9}, { 3.566998470788306*^9, 3.5669984830811276`*^9}, {3.566998534405218*^9, 3.5669985995821323`*^9}, {3.5669988153461113`*^9, 3.566998936995125*^9}, 3.566998967134378*^9, 3.5669990224832754`*^9, {3.5669990973790064`*^9, 3.56699912179305*^9}, {3.5669993047345705`*^9, 3.5669993241410046`*^9}, { 3.566999416274767*^9, 3.566999426742385*^9}, {3.566999497301309*^9, 3.5669995032137194`*^9}, {3.566999549733001*^9, 3.566999551433404*^9}, { 3.5670567651502*^9, 3.567056829110312*^9}, {3.5670569077344503`*^9, 3.5670570287126627`*^9}, {3.56705707303234*^9, 3.5670571842137356`*^9}, { 3.5670572582982655`*^9, 3.567057292103525*^9}, {3.56705732891959*^9, 3.5670574906294737`*^9}, {3.5670578301328697`*^9, 3.56705795512029*^9}, { 3.5670580197668033`*^9, 3.567058030437222*^9}, {3.5670580810281105`*^9, 3.5670581611654515`*^9}, {3.5670581932547083`*^9, 3.567058240616391*^9}, { 3.5670582801312604`*^9, 3.5670582990540934`*^9}, {3.567058343560972*^9, 3.5670584018894744`*^9}, {3.56705845049916*^9, 3.5670585083128614`*^9}, { 3.567058578731385*^9, 3.567058593645011*^9}, {3.5670586348134837`*^9, 3.5670586577143235`*^9}, {3.567058689007979*^9, 3.5670587629209085`*^9}, { 3.567058832777831*^9, 3.5670588394702425`*^9}, {3.567058955628047*^9, 3.5670589616496572`*^9}, {3.5670592284569263`*^9, 3.567059354442747*^9}, { 3.5670593898548098`*^9, 3.567059560456709*^9}, {3.5670596222796183`*^9, 3.5670596362884426`*^9}, {3.5670596726677065`*^9, 3.567059775347087*^9}, { 3.567059838449198*^9, 3.5670598453756094`*^9}, {3.5670598764976645`*^9, 3.5670598805380716`*^9}, 3.567059911254525*^9, {3.5670599744502363`*^9, 3.567060007818695*^9}, {3.567060059127185*^9, 3.5670600818564253`*^9}, { 3.5670601143824825`*^9, 3.567060314593234*^9}, {3.5670607181347427`*^9, 3.5670607311607656`*^9}, {3.567060785573661*^9, 3.567061232998047*^9}, { 3.5670626282800975`*^9, 3.56706267537658*^9}, 3.567062795371991*^9, { 3.567069851801585*^9, 3.567069879460434*^9}, {3.567070767429593*^9, 3.5670707823588195`*^9}, {3.5670805354535503`*^9, 3.5670806312689185`*^9}, {3.567083393113369*^9, 3.5670834586022844`*^9}, { 3.567083501299559*^9, 3.5670835804696984`*^9}, {3.5670836156477604`*^9, 3.5670836215757704`*^9}, {3.567083656504232*^9, 3.5670836570034328`*^9}, { 3.5670837376399746`*^9, 3.567083757514409*^9}}], Cell["\<\ (* 15 clouds: coord *) xc={0,-dcnx,dcox,-dnhx,-dnhx,dchx,-dccnx,dccox,dolx,-dcnx,-dcnx,dccox,dccox,\ dom1x,dom2x}; yc={0,0,dcoy,dnhy,-dnhy,-dchy,0,dccoy,doly,0,0,dccoy,dccoy,dom1y,dom2y}; zc={0,0,0,0,0,0,0,0,0,NM,-NM,OC,-OC,0,0};\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5668390630515456`*^9, 3.566839078760773*^9}, { 3.566839318626795*^9, 3.5668393238528037`*^9}, {3.566839506201524*^9, 3.5668395108503323`*^9}, {3.566839541207986*^9, 3.566839568664034*^9}, { 3.566839684650237*^9, 3.56683970887708*^9}, {3.5668397799352045`*^9, 3.5668398201988754`*^9}, {3.566839877872177*^9, 3.5668398909606*^9}, { 3.5668404106911125`*^9, 3.566840417539525*^9}, {3.5668404815464373`*^9, 3.566840691460406*^9}, {3.566840732519678*^9, 3.5668408608923035`*^9}, { 3.566840930000425*^9, 3.5668409412168446`*^9}, {3.566841374686006*^9, 3.566841375388007*^9}, {3.5668421835006266`*^9, 3.5668421845770283`*^9}, { 3.566887398037596*^9, 3.566887402078003*^9}, {3.5668874322328568`*^9, 3.5668874455864797`*^9}, {3.566887493166563*^9, 3.5668874939309645`*^9}, { 3.5668882488786907`*^9, 3.5668882511094947`*^9}, {3.56688847075788*^9, 3.566888479322295*^9}, {3.5668885903476906`*^9, 3.5668885954956994`*^9}, { 3.566971858103361*^9, 3.566971918927868*^9}, {3.566987354128765*^9, 3.5669873624903803`*^9}, {3.566987488164201*^9, 3.566987535619484*^9}, { 3.566987795843541*^9, 3.5669878015843515`*^9}, {3.566987856168847*^9, 3.56698786325126*^9}, {3.566988264733565*^9, 3.566988269819174*^9}, { 3.566988366133743*^9, 3.5669883719525533`*^9}, {3.5669885431160536`*^9, 3.566988549309265*^9}, {3.566991946182942*^9, 3.5669919483981457`*^9}, 3.5669919866182127`*^9, {3.5669920226542764`*^9, 3.5669920505783253`*^9}, { 3.566992103743219*^9, 3.5669921079396257`*^9}, 3.5669921962201815`*^9, { 3.5669922303998413`*^9, 3.5669922340658474`*^9}, {3.566992291177548*^9, 3.5669922926751504`*^9}, {3.566992333890423*^9, 3.5669923348264246`*^9}, 3.566992409675356*^9, {3.566997472043352*^9, 3.566997587171554*^9}, { 3.566997674297707*^9, 3.5669977903463106`*^9}, {3.566997868580448*^9, 3.5669978780652647`*^9}, {3.566998052270771*^9, 3.566998071786405*^9}, { 3.5669981665409718`*^9, 3.5669981730305834`*^9}, {3.56699824522751*^9, 3.566998326394452*^9}, {3.5669983623837156`*^9, 3.5669984017425847`*^9}, { 3.566998470788306*^9, 3.5669984830811276`*^9}, {3.566998534405218*^9, 3.5669985995821323`*^9}, {3.5669988153461113`*^9, 3.566998936995125*^9}, 3.566998967134378*^9, 3.5669990224832754`*^9, {3.5669990973790064`*^9, 3.56699912179305*^9}, {3.5669993047345705`*^9, 3.5669993241410046`*^9}, { 3.566999416274767*^9, 3.566999426742385*^9}, {3.566999497301309*^9, 3.5669995032137194`*^9}, {3.566999549733001*^9, 3.566999551433404*^9}, { 3.5670567651502*^9, 3.567056829110312*^9}, {3.5670569077344503`*^9, 3.5670570287126627`*^9}, {3.56705707303234*^9, 3.5670571842137356`*^9}, { 3.5670572582982655`*^9, 3.567057292103525*^9}, {3.56705732891959*^9, 3.5670574906294737`*^9}, {3.5670578301328697`*^9, 3.56705795512029*^9}, { 3.5670580197668033`*^9, 3.567058030437222*^9}, {3.5670580810281105`*^9, 3.5670581611654515`*^9}, {3.5670581932547083`*^9, 3.567058240616391*^9}, { 3.5670582801312604`*^9, 3.5670582990540934`*^9}, {3.567058343560972*^9, 3.5670584018894744`*^9}, {3.56705845049916*^9, 3.5670585083128614`*^9}, { 3.567058578731385*^9, 3.567058593645011*^9}, {3.5670586348134837`*^9, 3.5670586577143235`*^9}, {3.567058689007979*^9, 3.5670587629209085`*^9}, { 3.567058832777831*^9, 3.5670588394702425`*^9}, {3.567058955628047*^9, 3.5670589616496572`*^9}, {3.5670592284569263`*^9, 3.567059354442747*^9}, { 3.5670593898548098`*^9, 3.567059560456709*^9}, {3.5670596222796183`*^9, 3.5670596362884426`*^9}, {3.5670596726677065`*^9, 3.567059775347087*^9}, { 3.567059838449198*^9, 3.5670598453756094`*^9}, {3.5670598764976645`*^9, 3.5670598805380716`*^9}, 3.567059911254525*^9, {3.5670599744502363`*^9, 3.567060007818695*^9}, {3.567060059127185*^9, 3.5670600818564253`*^9}, { 3.5670601143824825`*^9, 3.567060314593234*^9}, {3.5670607181347427`*^9, 3.5670607311607656`*^9}, {3.567060785573661*^9, 3.567061232998047*^9}, { 3.5670626282800975`*^9, 3.56706267537658*^9}, 3.567062795371991*^9, { 3.567069851801585*^9, 3.567069879460434*^9}, {3.567070767429593*^9, 3.5670707823588195`*^9}, {3.5670805354535503`*^9, 3.5670806312689185`*^9}, {3.567083393113369*^9, 3.5670834586022844`*^9}, { 3.567083501299559*^9, 3.5670835804696984`*^9}, {3.5670836156477604`*^9, 3.5670836215757704`*^9}, {3.567083656504232*^9, 3.5670836570034328`*^9}, { 3.5670837376399746`*^9, 3.5670837651584225`*^9}}], Cell["\<\ (* epot terms *) vne=vne-2*(3-(pnh/NH)^2)/NH; vne=vne-(3-(pch/CH)^2)/CH; vnn = 0.0;\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5668390630515456`*^9, 3.566839078760773*^9}, { 3.566839318626795*^9, 3.5668393238528037`*^9}, {3.566839506201524*^9, 3.5668395108503323`*^9}, {3.566839541207986*^9, 3.566839568664034*^9}, { 3.566839684650237*^9, 3.56683970887708*^9}, {3.5668397799352045`*^9, 3.5668398201988754`*^9}, {3.566839877872177*^9, 3.5668398909606*^9}, { 3.5668404106911125`*^9, 3.566840417539525*^9}, {3.5668404815464373`*^9, 3.566840691460406*^9}, {3.566840732519678*^9, 3.5668408608923035`*^9}, { 3.566840930000425*^9, 3.5668409412168446`*^9}, {3.566841374686006*^9, 3.566841375388007*^9}, {3.5668421835006266`*^9, 3.5668421845770283`*^9}, { 3.566887398037596*^9, 3.566887402078003*^9}, {3.5668874322328568`*^9, 3.5668874455864797`*^9}, {3.566887493166563*^9, 3.5668874939309645`*^9}, { 3.5668882488786907`*^9, 3.5668882511094947`*^9}, {3.56688847075788*^9, 3.566888479322295*^9}, {3.5668885903476906`*^9, 3.5668885954956994`*^9}, { 3.566971858103361*^9, 3.566971918927868*^9}, {3.566987354128765*^9, 3.5669873624903803`*^9}, {3.566987488164201*^9, 3.566987535619484*^9}, { 3.566987795843541*^9, 3.5669878015843515`*^9}, {3.566987856168847*^9, 3.56698786325126*^9}, {3.566988264733565*^9, 3.566988269819174*^9}, { 3.566988366133743*^9, 3.5669883719525533`*^9}, {3.5669885431160536`*^9, 3.566988549309265*^9}, {3.566991946182942*^9, 3.5669919483981457`*^9}, 3.5669919866182127`*^9, {3.5669920226542764`*^9, 3.5669920505783253`*^9}, { 3.566992103743219*^9, 3.5669921079396257`*^9}, 3.5669921962201815`*^9, { 3.5669922303998413`*^9, 3.5669922340658474`*^9}, {3.566992291177548*^9, 3.5669922926751504`*^9}, {3.566992333890423*^9, 3.5669923348264246`*^9}, 3.566992409675356*^9, {3.566997472043352*^9, 3.566997587171554*^9}, { 3.566997674297707*^9, 3.5669977903463106`*^9}, {3.566997868580448*^9, 3.5669978780652647`*^9}, {3.566998052270771*^9, 3.566998071786405*^9}, { 3.5669981665409718`*^9, 3.5669981730305834`*^9}, {3.56699824522751*^9, 3.566998326394452*^9}, {3.5669983623837156`*^9, 3.5669984017425847`*^9}, { 3.566998470788306*^9, 3.5669984830811276`*^9}, {3.566998534405218*^9, 3.5669985995821323`*^9}, {3.5669988153461113`*^9, 3.566998936995125*^9}, 3.566998967134378*^9, 3.5669990224832754`*^9, {3.5669990973790064`*^9, 3.56699912179305*^9}, {3.5669993047345705`*^9, 3.5669993241410046`*^9}, { 3.566999416274767*^9, 3.566999426742385*^9}, {3.566999497301309*^9, 3.5669995032137194`*^9}, {3.566999549733001*^9, 3.566999551433404*^9}, { 3.5670567651502*^9, 3.567056829110312*^9}, {3.5670569077344503`*^9, 3.5670570287126627`*^9}, {3.56705707303234*^9, 3.5670571842137356`*^9}, { 3.5670572582982655`*^9, 3.567057292103525*^9}, {3.56705732891959*^9, 3.5670574906294737`*^9}, {3.5670578301328697`*^9, 3.56705795512029*^9}, { 3.5670580197668033`*^9, 3.567058030437222*^9}, {3.5670580810281105`*^9, 3.5670581611654515`*^9}, {3.5670581932547083`*^9, 3.567058240616391*^9}, { 3.5670582801312604`*^9, 3.5670582990540934`*^9}, {3.567058343560972*^9, 3.5670584018894744`*^9}, {3.56705845049916*^9, 3.5670585083128614`*^9}, { 3.567058578731385*^9, 3.567058593645011*^9}, {3.5670586348134837`*^9, 3.5670586577143235`*^9}, {3.567058689007979*^9, 3.5670587629209085`*^9}, { 3.567058832777831*^9, 3.5670588394702425`*^9}, {3.567058955628047*^9, 3.5670589616496572`*^9}, {3.5670592284569263`*^9, 3.567059354442747*^9}, { 3.5670593898548098`*^9, 3.567059560456709*^9}, {3.5670596222796183`*^9, 3.5670596362884426`*^9}, {3.5670596726677065`*^9, 3.567059775347087*^9}, { 3.567059838449198*^9, 3.5670598453756094`*^9}, {3.5670598764976645`*^9, 3.5670598805380716`*^9}, 3.567059911254525*^9, {3.5670599744502363`*^9, 3.567060007818695*^9}, {3.567060059127185*^9, 3.5670600818564253`*^9}, { 3.5670601143824825`*^9, 3.567060314593234*^9}, {3.5670607181347427`*^9, 3.5670607311607656`*^9}, {3.567060785573661*^9, 3.567061232998047*^9}, { 3.5670626282800975`*^9, 3.56706267537658*^9}, 3.567062795371991*^9, { 3.567069851801585*^9, 3.567069879460434*^9}, {3.567070767429593*^9, 3.5670707823588195`*^9}, {3.5670805354535503`*^9, 3.5670806312689185`*^9}, {3.567083393113369*^9, 3.5670834586022844`*^9}, { 3.567083501299559*^9, 3.5670835804696984`*^9}, {3.5670836156477604`*^9, 3.5670836215757704`*^9}, {3.567083656504232*^9, 3.5670836570034328`*^9}, { 3.5670837376399746`*^9, 3.5670837651584225`*^9}}], Cell["\<\ (*ww*) For[i = 1, i < nc, i++, For[j = i+1, j < nc+1, j++, vee = vee + \ oc[[i]]*oc[[j]]/Sqrt[(xc[[i]]-xc[[j]])^2+(yc[[i]]-yc[[j]])^2+(zc[[i]]-zc[[j]])\ ^2]]]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5668390630515456`*^9, 3.566839083924382*^9}, { 3.5668410500738354`*^9, 3.5668411059531336`*^9}, {3.5668412393245683`*^9, 3.5668412393713684`*^9}, {3.5668412896190567`*^9, 3.5668413118646955`*^9}, { 3.566841703815384*^9, 3.5668417040025845`*^9}}], Cell["\<\ (*nn*) For[i = 1, i < nn, i++, For[j = i+1, j < nn+1, j++, vnn = vnn + \ ch[[i]]*ch[[j]]/Sqrt[(xn[[i]]-xn[[j]])^2+(yn[[i]]-yn[[j]])^2+(zn[[i]]-zn[[j]])\ ^2]]]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5668390630515456`*^9, 3.5668390890567913`*^9}, { 3.566841524165468*^9, 3.566841530218279*^9}, {3.5670590745470557`*^9, 3.567059089304682*^9}}], Cell["\<\ (*nw*) For[i = 1, i < nc+1, i++, For[j = 1, j < nn+1, j++, If[i \[NotEqual] j, vne = vne + \ oc[[i]]*ch[[j]]/Sqrt[(xc[[i]]-xn[[j]])^2+(yc[[i]]-yn[[j]])^2+(zc[[i]]-zn[[j]])\ ^2]]]]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5668390630515456`*^9, 3.5668390979332066`*^9}, 3.5668412394337683`*^9, {3.566841289743857*^9, 3.566841311927096*^9}, 3.5668415580955276`*^9, {3.5668416231788425`*^9, 3.5668416346760626`*^9}, { 3.5669003862430363`*^9, 3.5669003862430363`*^9}, {3.5670591033447065`*^9, 3.5670591078063145`*^9}}], Cell[CellGroupData[{ Cell["\<\ Epot=vne+vee+vnn; func=Ekin+Epot; t = FindMinimum[func,{Cr,0.26},{Nr,0.22},{Oa,0.186},{NH,1.11},{CH,1.12},{CO,0.\ 92},{CN,1.03},{OL,0.905},{NM,0.765},{OC,1.36},{OM,0.86},{w1,2.15},{w2,1.0},{\ w3,1.12},{pch,0.51},{pnh,0.5},{Method -> \"Newton\"}, {MaxIterations -> 800}]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5668390630515456`*^9, 3.566839104968819*^9}, { 3.5668402974973135`*^9, 3.5668403782742558`*^9}, {3.5668414049812593`*^9, 3.5668414375229163`*^9}, 3.566841756605877*^9, {3.5668429499163733`*^9, 3.566842978277223*^9}, {3.5668880399631233`*^9, 3.5668880687763743`*^9}, { 3.56697197120356*^9, 3.5669719847755833`*^9}, {3.5669875623735313`*^9, 3.566987566460738*^9}, 3.566987810538767*^9, {3.566987882142893*^9, 3.5669878871349015`*^9}, {3.566987967677843*^9, 3.5669880000323*^9}, { 3.5669882132846746`*^9, 3.5669882156558785`*^9}, 3.5669882806455927`*^9, { 3.5669883912029867`*^9, 3.5669884039482093`*^9}, 3.5669885600420837`*^9, { 3.5670614496356277`*^9, 3.5670616095515084`*^9}, {3.56706164472957*^9, 3.567061671780018*^9}, {3.5670617748181987`*^9, 3.567061891771604*^9}, { 3.5670620082414083`*^9, 3.5670620167278233`*^9}, {3.5670640067297187`*^9, 3.567064010270925*^9}, {3.5670640525313993`*^9, 3.567064055854205*^9}, 3.567064094339473*^9, {3.5670689484819984`*^9, 3.5670689607748203`*^9}, 3.5670698952320614`*^9, {3.56708497059774*^9, 3.567084977087352*^9}, { 3.567085213599367*^9, 3.5670852679342623`*^9}, {3.567085494899061*^9, 3.567085499485469*^9}, {3.5670859700758953`*^9, 3.5670859708870974`*^9}, { 3.567086006080759*^9, 3.5670861172933545`*^9}, {3.5670866938391666`*^9, 3.567086753244071*^9}, {3.5670871411699524`*^9, 3.5670871501399684`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "168.06916898729622`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Cr", "\[Rule]", "0.2597231856403463`"}], ",", RowBox[{"Nr", "\[Rule]", "0.22405135815353058`"}], ",", RowBox[{"Oa", "\[Rule]", "0.1916417225081855`"}], ",", RowBox[{"NH", "\[Rule]", "1.1619297579240278`"}], ",", RowBox[{"CH", "\[Rule]", "1.2532410401053755`"}], ",", RowBox[{"CO", "\[Rule]", "0.9166962134031293`"}], ",", RowBox[{"CN", "\[Rule]", "1.0305980699531003`"}], ",", RowBox[{"OL", "\[Rule]", "0.8932960411034288`"}], ",", RowBox[{"NM", "\[Rule]", "0.7472908218886891`"}], ",", RowBox[{"OC", "\[Rule]", "1.5679221435272637`"}], ",", RowBox[{"OM", "\[Rule]", "0.8346241451827114`"}], ",", RowBox[{"w1", "\[Rule]", "2.1564408390806467`"}], ",", RowBox[{"w2", "\[Rule]", "0.9700759023513378`"}], ",", RowBox[{"w3", "\[Rule]", "1.1571282820932742`"}], ",", RowBox[{"pch", "\[Rule]", "0.572080058406644`"}], ",", RowBox[{"pnh", "\[Rule]", "0.5183973604881592`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.567230369521126*^9, 3.567244942692322*^9, 3.567245196036767*^9, 3.5672453996015244`*^9, 3.5672454374471908`*^9, 3.567245473732855*^9, 3.5672456740372066`*^9, 3.567245720447288*^9, 3.5672457602585583`*^9, 3.567245810989847*^9, 3.567245974041333*^9, 3.567246007487792*^9, 3.5672462443898087`*^9, 3.567246316617935*^9, 3.5672463763972397`*^9}] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Ekin /. t[[2]] vne /. t[[2]] vee /. t[[2]] vnn /. t[[2]] -Epot/Ekin /. t[[2]] w1*rad /. t[[2]] w2*rad /. t[[2]] w3*rad /. t[[2]] (Cr+Oa+2*CO)*bohr /. t[[2]] (Cr+Nr+2*CN)*bohr /. t[[2]] (Cr+CH+pch)*bohr /. t[[2]] (Nr+NH+pnh)*bohr /. t[[2]]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5668390630515456`*^9, 3.566839112846833*^9}, { 3.566842810779729*^9, 3.566842811091729*^9}, {3.5669879418285975`*^9, 3.5669879480686083`*^9}, {3.5670621290636206`*^9, 3.567062152463662*^9}, { 3.567063905688341*^9, 3.5670639108675504`*^9}, {3.5670689307603674`*^9, 3.5670689312127686`*^9}, {3.5670696234639845`*^9, 3.5670696694216647`*^9}, { 3.5670795267235785`*^9, 3.5670795309823856`*^9}, {3.567079591604092*^9, 3.567079592524494*^9}, {3.567079628810158*^9, 3.567079717433913*^9}, { 3.567082843118803*^9, 3.567082848547613*^9}, {3.567082895987296*^9, 3.567082897079298*^9}, {3.567085039347061*^9, 3.5670851279708166`*^9}}], Cell[BoxData["168.06916898729625`"], "Output", CellChangeTimes->{3.5672303695367255`*^9, 3.567244942707922*^9, 3.567245196052367*^9, 3.5672453996171246`*^9, 3.567245437462791*^9, 3.567245473748455*^9, 3.5672456740528064`*^9, 3.567245720447288*^9, 3.567245760274158*^9, 3.5672458110054474`*^9, 3.5672459740569334`*^9, 3.567246007503392*^9, 3.5672462444054084`*^9, 3.567246316617935*^9, 3.56724637641284*^9}], Cell[BoxData[ RowBox[{"-", "538.3611060012901`"}]], "Output", CellChangeTimes->{3.5672303695367255`*^9, 3.567244942707922*^9, 3.567245196052367*^9, 3.5672453996171246`*^9, 3.567245437462791*^9, 3.567245473748455*^9, 3.5672456740528064`*^9, 3.567245720447288*^9, 3.567245760274158*^9, 3.5672458110054474`*^9, 3.5672459740569334`*^9, 3.567246007503392*^9, 3.5672462444054084`*^9, 3.567246316617935*^9, 3.56724637641284*^9}], Cell[BoxData["130.60649695673055`"], "Output", CellChangeTimes->{3.5672303695367255`*^9, 3.567244942707922*^9, 3.567245196052367*^9, 3.5672453996171246`*^9, 3.567245437462791*^9, 3.567245473748455*^9, 3.5672456740528064`*^9, 3.567245720447288*^9, 3.567245760274158*^9, 3.5672458110054474`*^9, 3.5672459740569334`*^9, 3.567246007503392*^9, 3.5672462444054084`*^9, 3.567246316617935*^9, 3.56724637641284*^9}], Cell[BoxData["71.6162710699671`"], "Output", CellChangeTimes->{3.5672303695367255`*^9, 3.567244942707922*^9, 3.567245196052367*^9, 3.5672453996171246`*^9, 3.567245437462791*^9, 3.567245473748455*^9, 3.5672456740528064`*^9, 3.567245720447288*^9, 3.567245760274158*^9, 3.5672458110054474`*^9, 3.5672459740569334`*^9, 3.567246007503392*^9, 3.5672462444054084`*^9, 3.567246316617935*^9, 3.56724637641284*^9}], Cell[BoxData["1.9999999999999996`"], "Output", CellChangeTimes->{3.5672303695367255`*^9, 3.567244942707922*^9, 3.567245196052367*^9, 3.5672453996171246`*^9, 3.567245437462791*^9, 3.567245473748455*^9, 3.5672456740528064`*^9, 3.567245720447288*^9, 3.567245760274158*^9, 3.5672458110054474`*^9, 3.5672459740569334`*^9, 3.567246007503392*^9, 3.5672462444054084`*^9, 3.567246316617935*^9, 3.56724637641284*^9}], Cell[BoxData["123.55495989898013`"], "Output", CellChangeTimes->{3.5672303695367255`*^9, 3.567244942707922*^9, 3.567245196052367*^9, 3.5672453996171246`*^9, 3.567245437462791*^9, 3.567245473748455*^9, 3.5672456740528064`*^9, 3.567245720447288*^9, 3.567245760274158*^9, 3.5672458110054474`*^9, 3.5672459740569334`*^9, 3.567246007503392*^9, 3.5672462444054084`*^9, 3.567246316617935*^9, 3.56724637642844*^9}], Cell[BoxData["55.581255484423735`"], "Output", CellChangeTimes->{3.5672303695367255`*^9, 3.567244942707922*^9, 3.567245196052367*^9, 3.5672453996171246`*^9, 3.567245437462791*^9, 3.567245473748455*^9, 3.5672456740528064`*^9, 3.567245720447288*^9, 3.567245760274158*^9, 3.5672458110054474`*^9, 3.5672459740569334`*^9, 3.567246007503392*^9, 3.5672462444054084`*^9, 3.567246316617935*^9, 3.56724637642844*^9}], Cell[BoxData["66.29856748259418`"], "Output", CellChangeTimes->{3.5672303695367255`*^9, 3.567244942707922*^9, 3.567245196052367*^9, 3.5672453996171246`*^9, 3.567245437462791*^9, 3.567245473748455*^9, 3.5672456740528064`*^9, 3.567245720447288*^9, 3.567245760274158*^9, 3.5672458110054474`*^9, 3.5672459740569334`*^9, 3.567246007503392*^9, 3.5672462444054084`*^9, 3.567246316617935*^9, 3.56724637642844*^9}], Cell[BoxData["1.2090410322393712`"], "Output", CellChangeTimes->{3.5672303695367255`*^9, 3.567244942707922*^9, 3.567245196052367*^9, 3.5672453996171246`*^9, 3.567245437462791*^9, 3.567245473748455*^9, 3.5672456740528064`*^9, 3.567245720447288*^9, 3.567245760274158*^9, 3.5672458110054474`*^9, 3.5672459740569334`*^9, 3.567246007503392*^9, 3.5672462444054084`*^9, 3.567246316617935*^9, 3.56724637642844*^9}], Cell[BoxData["1.346739951488356`"], "Output", CellChangeTimes->{3.5672303695367255`*^9, 3.567244942707922*^9, 3.567245196052367*^9, 3.5672453996171246`*^9, 3.567245437462791*^9, 3.567245473748455*^9, 3.5672456740528064`*^9, 3.567245720447288*^9, 3.567245760274158*^9, 3.5672458110054474`*^9, 3.5672459740569334`*^9, 3.567246007503392*^9, 3.5672462444054084`*^9, 3.567246316617935*^9, 3.56724637642844*^9}], Cell[BoxData["1.1033574791548966`"], "Output", CellChangeTimes->{3.5672303695367255`*^9, 3.567244942707922*^9, 3.567245196052367*^9, 3.5672453996171246`*^9, 3.567245437462791*^9, 3.567245473748455*^9, 3.5672456740528064`*^9, 3.567245720447288*^9, 3.567245760274158*^9, 3.5672458110054474`*^9, 3.5672459740569334`*^9, 3.567246007503392*^9, 3.5672462444054084`*^9, 3.567246316617935*^9, 3.56724637642844*^9}], Cell[BoxData["1.0077532890936167`"], "Output", CellChangeTimes->{3.5672303695367255`*^9, 3.567244942707922*^9, 3.567245196052367*^9, 3.5672453996171246`*^9, 3.567245437462791*^9, 3.567245473748455*^9, 3.5672456740528064`*^9, 3.567245720447288*^9, 3.567245760274158*^9, 3.5672458110054474`*^9, 3.5672459740569334`*^9, 3.567246007503392*^9, 3.5672462444054084`*^9, 3.567246316617935*^9, 3.56724637642844*^9}] }, Open ]], Cell[CellGroupData[{ Cell["\<\ plot1=Graphics[{Circle[{xc[[1]],yc[[1]]},Cr], \ Circle[{xc[[2]],yc[[2]]},Nr],Circle[{xc[[3]],yc[[3]]},Oa],Circle[{xc[[4]],yc[[\ 4]]},NH], Circle[{xc[[5]],yc[[5]]},NH],Circle[{xc[[6]],yc[[6]]},CH],Circle[{xc[[7]],yc[[\ 7]]},CN],Circle[{xc[[8]],yc[[8]]},CO],Circle[{xc[[9]],yc[[9]]},OL],Circle[{xc[\ [10]],yc[[10]]},NM],Circle[{xc[[11]],yc[[11]]},NM],Circle[{xc[[12]],yc[[12]]},\ OC],Circle[{xc[[13]],yc[[13]]},OC],Circle[{xc[[14]],yc[[14]]},OM],Circle[{xc[[\ 15]],yc[[15]]},OM],Disk[{xn[[1]],yn[[1]]},0.08], \ Disk[{xn[[2]],yn[[2]]},0.08], \ Disk[{xn[[3]],yn[[3]]},0.08],Disk[{xn[[4]],yn[[4]]},0.08], \ Disk[{xn[[5]],yn[[5]]},0.08], Disk[{xn[[6]],yn[[6]]},0.08]} ] /. t[[2]]; Show[plot1,{AspectRatio \[Rule] Automatic,Axes -> True,GridLines -> \ Automatic, PlotRange \[Rule] Automatic, Frame -> True}]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5668391703953342`*^9, 3.5668391835773573`*^9}, { 3.5668412395117683`*^9, 3.566841239776969*^9}, {3.5668412898218565`*^9, 3.566841290087057*^9}, {3.56684186652367*^9, 3.566841909423745*^9}, 3.5670689312751684`*^9, {3.5670699737625995`*^9, 3.567070022403485*^9}, { 3.5670703462132535`*^9, 3.5670703473676558`*^9}, {3.56707040660096*^9, 3.567070410220166*^9}, 3.567071022770842*^9, {3.567071058073704*^9, 3.567071077339738*^9}, {3.567071109569394*^9, 3.56707112992743*^9}, { 3.567071166259894*^9, 3.5670711898627357`*^9}, {3.56707606019129*^9, 3.567076240699207*^9}, {3.567076605131447*^9, 3.5670766430707135`*^9}, { 3.5670849912053766`*^9, 3.567085001423394*^9}}], Cell[BoxData[ GraphicsBox[{CircleBox[{0, 0}, 0.2597231856403463], CircleBox[{-2.5449706837000776`, 0}, 0.22405135815353058`], CircleBox[{1.291429196139679, 1.884761711991131}, 0.1916417225081855], CircleBox[{-3.200397178358574, 1.2212124156850317`}, 1.1619297579240278`], CircleBox[{-3.200397178358574, -1.2212124156850317`}, 1.1619297579240278`], CircleBox[{0.6081672572221739, -1.3853495355430014`}, 1.2532410401053755`], CircleBox[{-1.2903212555934467`, 0}, 1.0305980699531003`], CircleBox[{0.6649556762928184, 0.9704620296599826}, 0.9166962134031293], CircleBox[{1.9046760746972449`, 2.77975792251382}, 0.8932960411034288], CircleBox[{-2.5449706837000776`, 0}, 0.7472908218886891], CircleBox[{-2.5449706837000776`, 0}, 0.7472908218886891], CircleBox[{0.6649556762928184, 0.9704620296599826}, 1.5679221435272637`], CircleBox[{0.6649556762928184, 0.9704620296599826}, 1.5679221435272637`], CircleBox[{0.44483311166516193`, 2.4648450688802646`}, 0.8346241451827114], CircleBox[{2.138025280614196, 1.3046783551019976`}, 0.8346241451827114], DiskBox[{0, 0}, 0.08], DiskBox[{-2.5449706837000776`, 0}, 0.08], DiskBox[{1.291429196139679, 1.884761711991131}, 0.08], DiskBox[{-3.4455458033171533`, 1.6779814766359766`}, 0.08], DiskBox[{-3.4455458033171533`, -1.6779814766359766`}, 0.08], DiskBox[{0.8381266667787244, -1.9091760938454152`}, 0.08]}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, PlotRange->Automatic]], "Output", CellChangeTimes->{3.5672303695679255`*^9, 3.567244942739122*^9, 3.567245196099167*^9, 3.5672453996639247`*^9, 3.567245437493991*^9, 3.5672454737796545`*^9, 3.5672456740840063`*^9, 3.567245720478488*^9, 3.567245760305358*^9, 3.567245811036647*^9, 3.567245974088134*^9, 3.567246007534592*^9, 3.5672462444366083`*^9, 3.5672463166647353`*^9, 3.5672463764440403`*^9}] }, Open ]], Cell[CellGroupData[{ Cell["\<\ plot2=Graphics[{Circle[{xc[[1]],zc[[1]]},Cr], \ Circle[{xc[[2]],zc[[2]]},Nr],Circle[{xc[[3]],zc[[3]]},Oa],Circle[{xc[[4]],zc[[\ 4]]},NH], Circle[{xc[[5]],zc[[5]]},NH],Circle[{xc[[6]],zc[[6]]},CH],Circle[{xc[[7]],zc[[\ 7]]},CN],Circle[{xc[[8]],zc[[8]]},CO],Circle[{xc[[9]],zc[[9]]},OL],Circle[{xc[\ [10]],zc[[10]]},NM],Circle[{xc[[11]],zc[[11]]},NM],Circle[{xc[[12]],zc[[12]]},\ OC],Circle[{xc[[13]],zc[[13]]},OC],Circle[{xc[[14]],zc[[14]]},OM],Circle[{xc[[\ 15]],zc[[15]]},OM],Disk[{xn[[1]],zn[[1]]},0.08], \ Disk[{xn[[2]],zn[[2]]},0.08], \ Disk[{xn[[3]],zn[[3]]},0.08],Disk[{xn[[4]],zn[[4]]},0.08], \ Disk[{xn[[5]],zn[[5]]},0.08], Disk[{xn[[6]],zn[[6]]},0.08]} ] /. t[[2]]; Show[plot2,{AspectRatio \[Rule] Automatic,Axes -> True,GridLines -> \ Automatic, PlotRange \[Rule] Automatic, Frame -> True}]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5668391703953342`*^9, 3.5668391881949654`*^9}, { 3.5668412398393693`*^9, 3.5668412401201696`*^9}, {3.566841311459095*^9, 3.5668413117242956`*^9}, {3.5668419284557786`*^9, 3.5668419511070185`*^9}, {3.566842124953724*^9, 3.566842139087349*^9}, 3.5670689313531685`*^9, {3.5670705170023537`*^9, 3.567070548654809*^9}, { 3.567076321179748*^9, 3.567076324596154*^9}, {3.5670763612250185`*^9, 3.567076444107964*^9}, {3.5670764854792366`*^9, 3.5670764972884574`*^9}, { 3.567076691212398*^9, 3.5670767228336535`*^9}, {3.5670850083498063`*^9, 3.5670850144962173`*^9}}], Cell[BoxData[ GraphicsBox[{CircleBox[{0, 0}, 0.2597231856403463], CircleBox[{-2.5449706837000776`, 0}, 0.22405135815353058`], CircleBox[{1.291429196139679, 0}, 0.1916417225081855], CircleBox[{-3.200397178358574, 0}, 1.1619297579240278`], CircleBox[{-3.200397178358574, 0}, 1.1619297579240278`], CircleBox[{0.6081672572221739, 0}, 1.2532410401053755`], CircleBox[{-1.2903212555934467`, 0}, 1.0305980699531003`], CircleBox[{0.6649556762928184, 0}, 0.9166962134031293], CircleBox[{1.9046760746972449`, 0}, 0.8932960411034288], CircleBox[{-2.5449706837000776`, 0.7472908218886891}, 0.7472908218886891], CircleBox[{-2.5449706837000776`, -0.7472908218886891}, 0.7472908218886891], CircleBox[{0.6649556762928184, 1.5679221435272637`}, 1.5679221435272637`], CircleBox[{0.6649556762928184, -1.5679221435272637`}, 1.5679221435272637`], CircleBox[{0.44483311166516193`, 0}, 0.8346241451827114], CircleBox[{2.138025280614196, 0}, 0.8346241451827114], DiskBox[{0, 0}, 0.08], DiskBox[{-2.5449706837000776`, 0}, 0.08], DiskBox[{1.291429196139679, 0}, 0.08], DiskBox[{-3.4455458033171533`, 0}, 0.08], DiskBox[{-3.4455458033171533`, 0}, 0.08], DiskBox[{0.8381266667787244, 0}, 0.08]}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, PlotRange->Automatic]], "Output", CellChangeTimes->{3.5672303695835257`*^9, 3.567244942754722*^9, 3.567245196114767*^9, 3.5672453996795244`*^9, 3.567245437509591*^9, 3.5672454737796545`*^9, 3.5672456740996065`*^9, 3.567245720494088*^9, 3.567245760320958*^9, 3.567245811052247*^9, 3.5672459741037335`*^9, 3.5672460075501924`*^9, 3.5672462444522085`*^9, 3.567246316680335*^9, 3.5672463764596405`*^9},ImageCache->GraphicsData["CompressedBitmap", "\<\ eJztfQmQpVf110u/tbtneiaz792TyezZGDGEqWELqUBJGBZDASEEwkyGRBIS QkqgnEgIFEuh0KIYC4KJFgKFYEQkSIEiUv8/gqUBt0CNCq7gGve4f36/c8+5 95zvnW953a9fd6iheJN+933fvffs55577rmvvv2BO++45/YH7jpz+8Ir7r/9 vjvvOvPehRvvvT9val/Sal1yZ/65bqGFv7NWS/7J8v9nd+EftNT/fZT+025l FxZbJ889Sd+O0L+d1pPnWvy/k63FC1l4tKuaw//OPZlxN+qVk4utC1l2kHuX 5pOLF1rH/EePVPReMqFeOQjHpC270Fo8RwPEOV5YPEkTybInW4uLF47qR0+e az2ZfzvIk8kwInf55DkatOU9f8zv2ptFN0x6QX5S8E7Rv23GCHqnb9QWeiXK teZL3+0LNuOzmMm5Vpjl/tL3ekPvXVhcxDuXNX4lIkNwB1yAePRz/jog8HB3 nNt4xAL50uwZcYcZ1eiSug9UPxZ7OSl8E3umXpiKVdQOxG4+Y69fPT7z55Eo BPnrzKOxi/z7Ig95PD4XGC2Joj/dpgi+wun3qDN+10434D5AIqKfXbjAJM9/ Y/IfidMR+uG9clkZVbaOVvR/pGJe3XoFR21RC10QBREhLzYHkWwFhjzk6LaS R+s6ttPoNdfe+u8N+A/+aN1CpiA0tNpsFwKNp5y23sW2i22/5W1QKdklF+Xi YtvFttWUi/1O27FVahuCnX7esoZmyMsK0ybz9Okls3/ooYeGel5OmzfTunev qeGti3O1Y2wKbRufeeaZl37/+9+/85FHHml1HvrgBz/Yan/7la985fdf+tKX 4tOawr/8w32f/vSn8fDCL3/5S/Z4nxvArqN/p3f8+te/fu0TTzwBaL78pje9 CaA9OxgMntm4cSP+fuTOO+8UJLzy29/+NmDFpzWFf8MP05++77778PAvFxYW MMiPrrvuuide+9rX4tc3ffnLX87y32UVtOYQMaB/Z0B1IAKwPH3kyJFf79ix Q8MAcAfPPvtsmNL0dT/60Y803rIcXwD/sbe//e14HHiVldyaAznwaf/I008/ DXCfuuaaa0BvgAt+RntrX0AO+gJvv/2xxx4DnEAL0AO+AA4O1oDYp5ahyFFn lSAP5rQH6kC8ATmgAdTXPPXUU5c53fArQApeAW+D0O//2Mc+BsR4r2gETIff WzoE1Z4Y9AM9BvgSNAYNQT9A7NEvsEcP0H3s/e9/P6CFdIMFICXzNRAvITBG COni9d+lbunP7KabboJYZVdccUWWs17r9/BvuYhlOeHot3xeWT7H1k34rRsm MGCU61DQJBiOgyHAEbgDYwHNZBq852fS81AjIAvEkVTH7qZspeJHk2CrQ3bO 4A3wCI3pafhgZrrQnBAyUZGEk52NdIcfXJoErEE2euB7KAnAC7jRmpWAO5fA ldfACvje2t4I3KFo1CQ490AYGjoOkg4KgQs9JHuvbwyvA0owA7QqzOW2GoAn Yf7226mBIJhatyFkmyxiYAGBGG+Wk4ZsIbTBCIs9YhaqB4vnAoMAjIBVW5vX AExsWgAL7CxpwkFDmLYGUkFCwb/4gOqXrgGw9oY2TA0LCTAiK/96sFhrgESQ K/LNNjaBySymxwXJnsR0YBysCFqzDSHZEdoAASAhptvgPHdVnP1kqLPbYbp1 DWFiEwZGA8PRamluDXDcLvw7BUabawgJYwEuLniULNa6NQAJYxjYJRWHb01h Ym4Fl0J7k/GeXQMwbU9cA1yTJGxsCNPeJIGAibh1Zg3AxMYeOMa8SDNc2hCm sNZsYz0CCSTncHoNwLTVwgScS5CoFiY2ZYAF7xJMgzUA0xbLexTL2dwQpoUE U9SU/bUDE9a2iMnQt6Yw7U+8F3VErwlMK2JfedbQdqQVNLbqIOFIBCgaNbgX kZy0fSWJmcLqaVtDSMJyhPQ2IKFvHee5SfMZ6zMsZLFQp29NYbo8tMG+kpdA 39YATGx3xJ6YlWodTAcdbp1aAzDxGJBleJqke3eMBhPwEP3tNQST0XFNYeKg jXjeROO8hZ45MWE45hwZ2tkQDk49A6+B5wSGSYOw3hGZXQ1BOGJFZlVIwPMH K8Ds0bcR5w/Tj/ezBx54YHj+K2Ia1zmMs7vhrDkHjTCez5q+XR1nOhmszybM IV5DXLOn4fylbc+ePeSE4tukFRDPHyEwxI3o26jzzw0fHEeCfdLz50WScSr2 jjj/vA1ralqITnr+vCAyhmnfSPPvoA3vg//Ah1evDghGhJcAAv6SPlYJBKM7 RwOhK22aEsEpnyAcvBSFKwGXgr7NLw0OTY7VgsPQYzQ4eroNNMFWOP57ZRGU FTFrfYcQC0sHAB8QAgQ5Hic9GVr0HVqMBkq/CIqxlldMmLc4DABNi2A3favK n6vrj5cRxoIeXx2YohUbE0zGqq4STNg/J126XJjaTn/HJgwTB24QzKGcM3y7 bPkwiWpYTZhMuGAMMEHXPCVevYZJ8ge8E22TyB/gkJW2KRLaWhK83J/Rr0cV vAN+qnAIboKg6pX5OECFiqLd0iKoox2io94oOehxGqVN6UO5taVvB7gN08z/ om+3cVu+YMlyLUDfDuJfTjFKWV3pBNoE8QyTSFGt5eKZRRLsKaktYYnPeOZF kzm65qWkjxtOxiVgBKz07fLlw2miZxpOURXeGbdJ0LXtyM8Y4DUrrWF4hw7E TSLXiEE1VmAMoBoroEGdhGVrO2K5HJh6jlgenjBMU44IHlw+TEYEVwkm462O ASbjTa8STEbUxwCT2agagmlFlsIMifFRxwAJIqaIOqb+JrgmZpjAbZQSMSaY sC7B+mQ1Oc5ou+XAxHEDo+0OrThMQ2ET+Oc4+nJ4+aCA3WKgfhVAEaU9BlBM tGIVQBHpHQMoxuVYBVBkMT8GUIyOXAVQxIIuGxTKC1l5BTYEgHjaKwTAmOzj yk3bW2usqFkcAkXW9uMGZeWloVMEJfHAcwqOIZKID/lbAMoYBcV4wuMCJeym 9+SsMT7APhgJFGhT0LeNYdFEiws6cRtcv9iNl2u+OhpviXjgmCUcFkk9hx84 FeAHKjAambp2gB+PoY3iguG5DlxhLP6BJ2zXrG+OlrEZ5bGihY9gAhZ0C9jg z4EfYHpD1z3EPORUE34GuuAvzawC+IFZ6cz8OOUkRFa64A94onIOj+LfW9Nv YBB4iQCFfvNSrUtQIPt7K+5h1jgFg+LEoBzkODpc8JDt1AHJwBDAB2D3spVX fl1g4jhLdx/iji9WXxBfTA/UpN2SkFnRAfDC31B/Xi5wxUIVizrJeRzHQtVs wCyd0SmBQnLpIca0+OUNMaAgjuPlCVfAa+Y3BngN/pYj2GfPngWsNDeeF8iJ Nlpienm2FXAaFhxDUMVs9C9Fk3OaIZW8kMSeIGp0IhgCSwUz8hZ6zztBzLCZ kOIYYDMhz6XAdiTNi+pdnD9/XsaLNTDQVgGT8QnHAJMJeS4FJs4n1mfL8Ptx RcacHYj1j5eDNRxfLAxNhzrYe8E4wBfxAZVP4R+igCnBa4oIY3mXgghOEDfb rUUmvfnmm0HkOqGEEoMHV4qN4Jv1ADXAB/JhtfEOnmPPF9Mo4Ak/TEEXtbpy Fhvt4GxCXlMtyUhDH4T6pSKNgTOZWzUaLEp1uatjtbe3iRT2nqfxDCQAUg2s AB7gA1IGuMiaBFhmgB+0gaXxjLwDtsV3sDDm4gkv76+ZLIah/bVS14YBkKN9 NYSKK2tR7UfL8WRCYB6eAp2pyo2UMQFu8V7kmZD23gP8wAOmmOV4keo25a65 zcw5FHHQnH04d33+V7/6FcAgUKbo3wGg/+sveclL/sorXvGKBx9++OFfzc/P P/zggw/i16t+9rOflXOZSS6NttarQ9R2uNfDYnhuqGKOlL8JeVY91GDR1XGA TPx8uAyBtFrSGR1N0RbOW3WACqDoHxw7dowWjP/vkksu6ZDi6f6d5z3vefgN k/nfnU7nv6xbt641Rb/1/+/U1BTGIHWzadMmPHf+wx/+MEnLUQZGjafPOnrA tB2tF9A4bUDhrsB8spAj0oRnB8UhpbYKSFNdpsfqsaVsfIcTrdNv/MpXvoLh /u2WLVsEQ8DOl9785jcD22BUN1s1oGAaKMBz2YkTJ/Ce7geMfPbzn/88G5ra wJP+2ct6CikvfQENGIMdMWuHQ6XU6mD66A8gAHUHyoc1am8p/Bqe6/++b33r W8KP/316eho8GvluNo0VE8KD9prxMIV+wNPgZfA/KCcYKTxuuvQw0nGUKVeD AWjgU7wPo+Lp4m5CqJTHg77dXz6S0ZtLQWg4B9QDO4I1gcz/2ev1wGKnfvjD H3I3xjf0VshhOj2o02I3aFtwXqlyYx1IjcutlnWAkLw1z6fkxBxgCWoWhgs0 nC8fxqwUloLQEC3qAqGQW3CocJWceAWXYjqlJ745f0pzOTgcmIzHI5r6Oo6O M6m4ARcdcbtBZcDo+S6cpidxBVosetPpOMN4mZWFKxOaZoaF17vABrAC7ADT pFE3OwLoYTiAQoW3oLbQx4XLL78cLB+PYdXZd6YmQIx1HPu2tiMdK/Lkpetw W0BlF216x5vw4ZEjcFoH42NO8OABT8X0jQwdG6KILVTfNHuNZQyiDpHHVEAK VRkhanwvjDibKPHomTNn0AWgJ02+y3neM2Ac/hcfHlMg3zOfOZiQwtqBMjNS ehTyh9/RD9pAhYjxfDxmGhPE2JMIju6JcE2LevF5SywZInG9qnTcn8kT0JQa 7fKKlKq7i0bpZKdPn8663S7OPqKO6DZqbqMgIBbn9G0HP4pCdLlMZO12O7v+ +usX0DyFEoHBvxrQA0j7xYu5m0+pvUA7Ojt79uwGfgxZv7feemt26tQpGjaX Neo3xzoOvPZ4Cvgpf4u+7eI2vAbHJ2/n1YxNLD7JbXfffTdCROiyNc9tmF1u 1imMchht3eXzOns84G+oC7gNUOx8ctzEqz1eDzakAy0Dfw4faDGvlB3PyORh 8jYGpFiKh2JIKCWIS+imD0kD4qAD8Qi4G6ykF4tmE25nUut4liSvaUW3sGNG BeEwHiQPf+9wHvXsvT7bshy6cN+gBWgCRSLKw+RBeURhNQ73EW4LmQ2v1l5A yYzJ08qCT4h3iqgOwr1Od6EJp8sf6mNZipVi7MirPO4Rg4tPgB8iIT3eYmyZ pc7SKTHkJcMBBBXgvggVjJ/qFSPiQkAQB1CMIoBbS6nQFwOvK2HGap+eTxJg 7utKm0KC6Cww1U0hKa9eXniOdsrE8OMTtpVVEJKsv7zkeSGXWuTQ5EPnQ2s6 UQT6GGVKOzfXs5STKvYKsP/pvn374Jnib7ZExh3xqMSyB1sPCmuXy6FSF4Bh CEMdby0QpmzqoP6xd7/73aT5CmWjKO9z2gExIDNWFtZOAPrzIqGhref5AsC3 tzG2KenlWFyLYrADonkv40KJ8oKE3oaOv7oEIhsHrP6rnTt3fuaee+7hI/WY YeQmjzCMJEQFQBjyYzeVIpoqcP7gxS9+MaKWh37xi194ZZeD1qEKkEDaf16/ fv0f/MQnPrHZkoMcmRkHoGAZ+iCXBBOk0jXwQe+F4DXFXvGcVDuFAhAHiKwZ psUQmgpgYYo9KTaO9/Eb9CvmQkHcXD/lQwNQNP+t5z//+f9m69atsah401Ow OfdiWrASol/MRp5HFNYn0IN4Vy8BHa+iC5aUiqqR4B5leP0EkP/r7OwsUIp5 bEzqC1xXWk1yNgmK9hGg7oDVrQmpwJIYrCzHHIim6BCrlgXlOpBQG9qlAjpk CQDhuYCjGIRBEwaGvfCKodetvpkTwetZvmog9O5L3BkDiF7+xvb0biRNeNoN +kjpetd0euGZAGoH7wEc0e5e2SGPROF1qhn+1Te84Q1Q1cDSH3nve9+LadCO TVgp4ivUXdyj2JhUKUp+gxZAMq2V2GDjb0j/H/jsZz+rQ8xKJnskd+zPYlhg Ce+EFKH+j6+99tpfHDp06F/u2rULbiwgev3Xv/51msLBGsJtssjP3vKWt0iV AbPe8gjH/hF0ZCTcBmvc8BH2BJa9I1LlNOuKpZWICGGOU6DM1qBX1TSY5N4N 3/3udxGbwGr2f+XLLSkDTd4BmcXWQPYuhT+e4TsPJCtpzkKFJkxL6n1zLWU4 Lr/Zvn07WOTFP/jBD8JEB5i8JLtAT2VcKRwAwXp30/7pn3jXu96FubbCbmAb PckcgvoZ2poKSqyjiUCGKyDE+HoeDXcmGkY/kUt9y3MG+V5oziMgh95EujQP iJSZPdn1DgE5kIl3/uxb3/pWWL3/0e/30efVP/3pTzn+XEyw09aKpzGUfyQV zEEWiMlcGkn7RxJiCr3M4isIB+UIYOKNEvmvGMFUIQuGYRp/g7P+0YEDB6A+ 3FnuidT0hVTpEbxL+w08pgirWcp5hN5lBZ28koB0cnQwF9GM3hHLchr39JJJ K+hphwUDomeLG0joAkOB6//mC17wAqgxuFykwrYMG0E8J8LnrSODaejoFQm4 LTAehaDlIhf04G27dhNsMjE8DjQX8rTiAajQ+7S8gkHf8YUvfMGbc6B4hVpW df7gfpKr4+2N7nVEyaM+Lyog3nGxIDrz1KlTmB9xIi9vZc6lpA+w9otpFZIv o/uSgmrknQS44gaipMbJTSV//O677waa6dnAsR2RVckjK9/yiP6xhD54USZR WfGR1b5/1Gvs9xUnRIDIclHnibC7LNmbOgRanHAApJ7YUPyV6RZ7Hd3rEXtP kt+i4wUGzG644YZAPkdAy4lt8TWVGEWw7R5FZyUNeIEXuVUDDCiX8jBcEfdm r9Nb9zEHYIy47uOxQV7QgsYPmmFQLA6p5wH6DJKsIcBJvBeQNi1LSwn/ulPc UUPeuUQ2yHL4Vk5eo3Y98tJzU9hmUZtAdFUJxXgde+ARdZ1DfEXUyOhtZ1ZM VLwLSOVSDVkPEvIXkkiK86McQGfXIq7jyfVtp6kUioF3ARqUrnhgulA4bZnJ rsqliW9jguKUP63tNTRk9wRWE4FJ0p8ea+51MOjRkNdGEM9ocniRrzeG8M7A UfEeQdcnCS+UrO0bjZEFj/xpyf9jYmIawLSALZJE0s5pBHgPUkuM4GmoMNc+ KIK+RGzwfDQLs4mhYObf86lPfUraxN2nPnm1D6phTNKkU848ttWQbl0S7+g8 ecWj9zoCUUE6REfBCrQRzDE2aJq4YcbkE4E03kUF/UxpX7ZKWqNGdyOolYFc 0CDCKZpWSn6H1IiOkJe0XD4su0J4++uvf/3rRaFxPEJ3qpllOr0llRrR2Var kMknVnojVh/f0oxUWBTFtNMKUpkCDh6p2DeGJ4mVlNZnBs0eK7ORMZt/HtlY u5uSucr9BBxgcZCQHxXpIClVJWeiHT2SOJGkIXeKWELoNBHaHn/b294GyZgp 63IqeeAdrRIwd46FCbfqRHwZk0i4uYZcs0kpYnFdWmB+b9LDtKwtI1cgZB9a gRxebjcF/T374aW9eqTyEmgY/4d//vOfY21J0rw+aaNoj6ZSG4UzQadjiQ2x Y0oSytEi4CKKQt4Sl3tO1+2Enoh6zjoxVlwUyvZSRdE2y+xNRfIR9Zho8ESG NrVKLFTcsvMCq4FoPawsoA7lWYnOVPhvZvfOo1dwwgbGn8w/mPbn7rjjDrkK gCGCUozu4FSSC9l4iHkbvEBQW0r94hDoNkbUe07/3aTxosveTWIV3U1Ps2xz nlNaKc4t0cwXPxYh41d6O1v7HPHwKEkkIL+ysCnt0bCbJhzXGF58nGmomVXH C0QWZ6xKUsa+o2mRsqGG8qm9gjVtywVEuY4zknf6yZNKT/lsdziBxxheCFUo 0+mEUChTis94oQ2mpnEHS6kZHF1kVP2T/fv342/PnHKbUTYeLTnUKBeKfOCj H/0oaOedzeqX6qiOtqiJnjHaaDRnp1zMTN79lDOcVwfYy+jxiEpLOjKhndRz nFVDUmJUCuaSBJST0rgxHimZFnBhkC+vseDR0UDn0ZFj/qe/8Y1v/LeZmRla zCj7FnNcmIimwk94bvg8TVKwQ7F/9Af8qSEi4ThGb4iUJSeXI4yGtbJkXDwk 7HRkT1lHmUzw2+qJCJtGGxtlRGQf0+wKe8+xPEKLYeMTMglnwwNBJWgj1EzT 9RK5mY7QoX/p1a9+NSGwPJZp27x6YV75w+DjDCVr6/e01LHHKUPJf2cdWSTx 8jwPdfNc1Lxe4vH6GgoyFYx99MRrPj0XLUcFBYUZRbwrKGimW07BjmTnitVV 3kssEc1derEZIx7lRHRD5EIloZRaUGA/EWtD8Ycxy3sXFxfJRniOxi5Hh3QS j0afe10zwhnRqyCcOaXoEW5hWMmIWYQIeu4Nt5k9Ro+AbBINoTM2HV7RN08H evnMTVNZ+bm/8aIXvegfX3bZZWrPi8aAuUBMXUgIciIJiL55JOR9AqM9GRfG DMzWkLCfyF6597+QWCump3vzYhIWDYGEzD0SKgcunhPxEvE5SmisJt+MLkU6 8E0Lo3dKe6pcazWl4l972cteho18RUXKWnjDV7/6VQFhNqmB77385S//+eHD h12e310ui4b5GhISGbaPIpOzhpBGPXmEnHeYnVEnCVNGABw0GV3iUVQl38R0 9wD7QM7l4psWTMGTp66aUpE1qnQl3RdSGfAtds8RbagYiY9GFvewt8dRKd45 p5kiVWn2bGyh+3DAivw+79JQXosbYfMMNdPSLCyZRljCUC6Boq+DLLs56J04 YlqatQnnnaJNTjiIZBYLSxjpb5pIzT6uwGXM8MBXRPKdf5ZcEvlZCF5BU8/Y D49cEQ5gCplgZwV1jQRWBIDM6oTnZQK0FdR165061DWmJcylK/P7y6961atk TzeMbxzHYVoOpetKTzwps3TpOciYKkcQWFAMgmbhWLC0IvzpCe3wKr5CFXPf JjzqEZiXfMY6eU7qvohM/hUxBzhFCrrqa0u8Y5cbkuAVzhqRl/T3jx8/jkQe 2Me28+gIdBWZU5a3KOqeNfaiPCbu1nWeq9iFMON2nLY6uqpQGo5603LPu5+a 6eocQ3EF14RXeYPKuNT5h97zcvIqKOsdH/jUe97zHpgU5Lu4Fd/LyWouW4p9 V9jh2jZVvCSmAvQcGazAnSerZoxBM5oa/9ej6YHy8TxeM0sdVpQmbURo6i3v vSPoG6IO4Ew59P7Nm266CTEKSWbGL2IM+w6wDiEjM3mHzbnNlBIObX2jXrJg OGOKed8BqoKInu87rJbqiEgBca/qHnuexnaVpnUQisP8p01AQgiWhWptMSTt FbbY4CgmznERKy1GSypWQKXDOHlHZ5y6GhE1XQc4r2w+8yCe+/1f+9rXZFNe nzD1+vXcSAOWt+bsN5M5LCEkQ9KVOSacCelW5OMY9zbAO4DLCX2nN/p1GMwr KDCXho0qIajPgcR99G+yfP+r119/fcDGtISI5NwsoMd/8R2CQkISGCduj4gW ZgVolDIT1Oz18HN/9P777/+PGzZs0Pm6BhGekmJCGqfGuz6njpD8DpAMg0jf PAk8mBgyhnU9b/TSpA8wf9I5HnOy866TcnBiVfaey89cWpUUiLpd867mbTi0 4kCjTZJZIKaS6y2ZVfgvvks+aSwCechSWdStV0Pfq9PP+ZKSGAkV98/27t3L 0lJrjebLFawZuiGhzcLfIzQX8TL3v4viZEGSREQ5OwZMyc4v/qYonlfgiCsY aCuK4zWlZSzmHKEJfawHA5Ky5flpPMo+Opjvn+/ZswfC9ZYvfvGLbnRhKs1T IgnM7JHi4t3xBoq2pXoanEIvKa2yHpC3Zx1gPIJzdMckA3hWvFckOMHFT+LU TAw+emTmonU6ewm/zyVmBfURKAEcXKAFSHqEb99Vi4phby3E0n569dVXU0W/ fACvkMOcIwKDxGIxtpV/5LtNoSQ8tg3CvWCuciojS6sAczSygb50yAFnIH54 6tQpjPUXX/Oa1/zOyZMn1SXuHXF4Ja/VW4TPOPqqguBeHHhYvVcEJfgd46N6 pGd3R0BAnQJ5DmgGKvVhU3kulBYLbSZWzxT/+Pve9z6c+aGChw/5VeDmHFA5 LcTEhLN0EiYMS+w1rfMqRZ/GZQeUCdDhea5qZzZyVjcpEjmoiSgvlvzgMjAv qA504NQX+iqezPEOgngxHI/cnuU2p8CG5Nsnt0ku8sh9NM0rk+o1s1Y5q/JR MR1R2Fg8EdbdAP9nV111FTyyaBcqSO3s1U6bNlaiUHfgqkfPnDmjyNOR3T+g lJLBHJWj8pnixqKqE4XhVP7fQM6vov2Oz33uc4AQIP2L3bt3gztUentc/zHC zNFXz49l4npJFObdbjPimtBueTnymGdFoKsrvgFiOqod0v4lRhP3gVRegGTu wBWES1gjxwYkb9nD2hJDQYYQK8WQnHFHLNkKx9Rjhj+7Ve71v/kUopooXNtJ kHANOXx/38c//nH0HZMZmc/wnURAJZvEJcKsw7MVVDbiWhrJqqeyl/3uUDmu HqVYO8HhlU2fS8wPI0bw8WpKEopiAmzeijgknSYuofb6xFAxp3TgoIBRBSbC c6g5KTmaAFEqF0r9FkzPSzpS1I7Y9Hb8eOtAznNjSNg+Ur7KQ47c4ZXnYqXt bat4gXzD8NyfMVzjovi2pOaiExiWViZKDNj1aSu3fI46ehTnzlKGVBaUBVBw d/QxK7zD6JClfqnCZzCM/avIC4MgAg0pbBuOqskUAbZ4p3KCSKebGjLucEBU tjeSjfWRWAm8DxRPEQtTfga2dL5z4403AjMs0l4aiPHER2QiEyz2zvjWMZEK pMd27woPPsyJBTCekw/2l/HuRucVRlDx1KZXTpxPbeI7nbTm3h/jQiDh58Hb Hn/8cUTCZQQ5sSm3EairY6mq4RK4y3hhjBuQUFwUMXBlM+WH0LWcNwPvI4r4 jdOnT5NW3V7dRRh1PWasKQUOgXJFAQHwKwAXw4KuaAOICsmGLYGia2F43CvB oE4axu2ECs4zhs1bYmjOk2pIUnUNd0F76ovHMxtaHjdy2QqsTAV93FM7rlTL 3lWn06WgihRjBjzMRKBUsWu9M8xT7d/4ne98B8whVVCB5qL/0ZQhnfiGy5BY TkQ56I02X2hzKZIiMZZYt29zaV+zpq/8QxXpWuHwC7j471555ZWooaAPe9TO qSIfZzms6QU8hxMRCvcsU0VAvjdb7pb2AgAjsijoFHHQCTiIVljR3ryr6oSD NlgJ/+tt27ZJ3R9VL4Z8xAKOZ7RnIdc5QV3AZYY3/oV3vOMdoDc795IjNDp7 6kxruEtFpnQm5yxdWutkfY0RpHSnBATRrovI49115bCDJqjXIHH2LPdI/tPc 3BzwhleFDzoNZ5gOZXnbreooR/HwwtgY88LiyZwvWfByJoXupDEyPHeLfTgw sfPw8rhYrkLQaIomsexdjuUgavqHPvKRjxQu0kJ6oxSZYsGYBv+AUaVws1yi hVNDEAJ9YhjMDAdrKUrV2Kl2Of9GJlJgxyB54t+5Yhq81CCRsgroB1wF8PE3 BBkFgdo4DtjqoiQYduLxDZVkQpiKfmsLitQUqBqOFOugCFmOHqrRGmbteci1 kFQkvZpFyYzz3MQ4fnlMTCuOAg4eldzIsneda89w3vAvvO51r4PWYpDkjJPu OiqJhEqqrgt0Qgvr69EGScyaOApe5rTnKPztEydOxFV2r+FUuUYnDr/C2xcH G9yimU508VxptzMfPn/+vBaDjLUIpokufu9PfvITYGE0JFqeWuvuQtNFelhf dcGjiBpDucZaeeX17vreLXRRu/EmP/AK5QsKxuqP7GI88MlPfhKrcffGuhHd 16brKchO3B/u2zmKS0pz5HxRCcpA9hDERIpLLOG0tQJGzh2WlRTmIjXyoC/+ zG233Yb32qkLuEIQbzyiurBHepRKLbJc7QG10mV7So4ccck0Wtyv9kY/dfQy 5vPw4h4WMua/qMW9vvmPZyP5akthGbPZUjBfRHKvcCmH4yXTQu6EAZtc++Mf /1hSBVQ4MO6JhldncCryd1/4whdKMT+5Tga/AkF/6p3vfCemle6K8fISxxC+ Gc51axYDjKuKcupXX22oVr5xlaZ2zeGEhbuuwsjFKxBZE+sAcNVmjgGjn7gu eiXqkBQ0IrSG5N3odyWBQNIwVGap2dhRpC+c4J8WqQUlJBQACsJ8or4jTck7 hNFJCjJivykHLCTNEk20V4ipIQeYDNGauP/QDY8DZ85hjdXVVURiyIeNh7kR MkNT4MI60nvlxIcXXa0O2hB8wwmplKRPJejBkCCTJBXoLZh4LYjay4u+QKGM iNrLozQVpHdg7w5cvP03v/mNTsiM4HspM7OOQFeQ3UsGNhG1hmSXVVAN2e1F l2p3OW7Q846tbNALGwgrcthTTuRjfx6mmKheva/nZD3OFE9hCCYhrB/9wAc+ INmUnP2MP6GtgACIuKpbbGVFETw6QQFRcWteDJXmVwz97zZv3gwNB2RCzets dGORZhy5rdiaNwF/7yBLQ0JL0dFSQns3fQYZjhk4gB/wcnk9gCDRFvCdtwRk 9wyP/sOjR49KtledhBs/KOBo1uwqZukostgbU3J/KlFeiplTpbTkglMkMQZd VFJGXP4GDpkVO4W5Q/OgPyTNISqEEw3/4dJLL431HT0bqGpcxO3JCoqbYJBK HYhBiMpkK3MMqCLjxtx6GvhmVgoYSjaibCpJAqX4eTqpwcmHo2toELWR/isI 7Ry+3mSAzT8yhEagJAhgajJtfZ0DMNJNagrdEU+qY4Nxt535C9BBgtENIipI uMEvEvUTyPXZXS+04VUE8ejsLRuHaV+fYwWrHQNPFel1xTthOYFSJxfruKmg TLO1sTpsTbFfgn2U2rtlmeJmRRf6mDPLo0ydIshSyiL+ljKmKGry8u9973uB Csadw0cWfFLIvZBCKdaDbZIeQE9ElaGRovjojLWc57O5p0ed/FkvM94Y/obp VubkUEX+rHsJ7rwDB2tr7QvG4qcVJDXhJM6CLiMh2pHoJHkjcocr0KtDRIBX 4j6yzVa4CDESVO6S8C69ZnI5sRxK1pBrd8FNxNnrknGIxqaClJ6wm/EbHl4w h/u8wwtNr+blwwvGrWF328Q5mQ0qr++dc9iVC32IWyILcJ4KInYQUGUh+jqY LVSV1RJ7HZGSovqZkoZxvSKlKisCqSNIsW+pzZ3isRKPgh6lzbgNKWiq/lQc P6m9N9grB8sJU4+eOXMGISYjyFn13cI8ZXMMgHMfxXeFPIpKBuhIj4N3Y4og a5Cdg19RxhkdJilFlWmKefJc9RT7vshskJQDDs6YtVHFWS/v8KcZuiH1jBGt OLDnlRTzqOcclpo2tRQU9aRcZOnFxl79Qo5kgKyIZ+IQFv4WrhV9V069Qq6N CkgZ3HsHsjgDQus5jCXJcF7otoKEngo1U6g7rsfKQvJE6Zt3PvqAwx4eCb26 ebySlOCUIaE6HxdvYyk9c1lwolidAldYDaKsB9DsHeNseJQ2cp4XDvRSHFlN F2sFS5lCsprMRsaKVpDUOz4wfLKinqRGYiqOvHspzt68DEt6h5B4m60Y5K66 zJmD215MCudpsRvw5i996Utuxmd5YYqhKQhKvRJRPUdGvcPx6tw9sEDLNq8m Z8Wh97KjX3EuDYsZmAPSFdUqzPwrqlV454yduvqj3ya90eF8jn9AYpCyBNuh asjURmxKSqsVy9UXVLjBvBegVqGHeLKIKe4cH6/3GjsO5EPUJch4hkY5erWC 9juzrqgv41Rmauv6Xa735N1wWVGBZDiw0urhvXjzXZZq+Hr7NSPWCBIlKMTk bBENqU4TV8eH4tEYL8rmIXFPucAOL+Mrogw8BxM/rCgEVXuv9nw5O2MBhIVQ +Oa869W9u9x5jje8nWD6tBSVliJCwu5GiXmV2UaktPSH3G30w9Wg9DCRhMwF UkpczIOYioriQR4vGiEeKgjlUxeBFYnHV1G39k5vdniLdZVlm8FztzxkVxDV +CkqzCu79pGxvIqYbQdFI5bc+8ob3/jGwu1oJNY4YSKU5WZZ6iHvnFSy0+vu clYzbXX12rhvcwKvonaic/7dJaQWU+hc6F5g2yOkpx/LayfOmDBJFjIJ4tmR LHk8XjHnwoGrkeio4npy0g59qwKY2Mu55zOf+YwsBDEk8meprYKItXWkGxbA NIvRCiIaBqmoXCrmRRyi0sqltAr0ZC/UKx3EamfPqsr8Waqi55Vo9iLlI9Yq FU6wRcNJvLBgkqH4qgRQELsPiGKUqU+uWOoFbI1uqCMYE8dE7yqKBddeTr4/ gYAqrCJtFaVmjQhUlJoVXimsI5xFrhtjxIwgrKbWRoNli0TSvEx93i2SKciw gYxd7SpF56eibLBxn704Xl3ZYKalqdJTQcvaO815bQr9gksANYadDFHrrHjV u0MWau/fb9q06Q9/6EMfknZZbqSEeg5vD1XNALCFouv4RHKX30fRMUY0M7eM 9M0SKgsx3DiuBzCn55ii7Z7Jq6vUPZMIUemT8BKx9vpzFj4xcwirgfkqaGbQ Ul5xvQ9swW8Q0IbLxLmHWDp6czUbpp6XS2SK4ni70Z4LpKr2Fqup1VaM2OEo EmWe4/w2NKMmbB/CMfp+Bo+atdekB2p2xfM0XOnB4dWIKb8Moa9xJWYxqneV 41G8aEdYRkkj7X7Isg/06zlTUff0FMEww/RSW7E4occGBrTtDlcr7ysK68Yi KfUFMybqUhGxqb1OnT1QsyjMUjZABQ3NdnL55SR9c3cTa0LsRxBdle8XbwIs HBemb7Ivevbs2ej/DYZ5BB8jmz2nf57/8DHIlne9m4WJT2l7CdPmSrdEsMrL nUygpYKMtRewz6f+4NNSBQlvaebB1HNYpeKCJ+9uL+gvmEzSTap4QOFiwHiA mcY5lgQm5kJxyWYpR2KqyA2cfr1EaCWchcoBFp6tae6RP9WZyBgkGLrpqUDL dUlXRpemIrhbe2M70ZIKI/J+iAk55x962ssGLafadDFWIOUQcbojOnRhK7pv LgXEeK2QrEKYPJK4CgqCMpvUdXd/+vbbb8e6APPlbMVidzIVTm0pRmRp9tss HXRGv4xMz9ddosa8aLY4KjYUay9lZ3k0gXWV/6gvJdNUqr2RPWB+pqjLjKxl oQ4I9FYAtQ8sSv1OgC1qjpjlYJpWTMySabH3r4upF6oCdaRvCavwclBrPnS1 JQkbiBLqcyVGjyNvbUYoo8AqCFV7jyQ/Z+6RZDdbGwi8o24yqbx0PUxxyJ6Z Zbq3IcTEitdvtkLJC6lNQ+g5kPhK8nWqrgDVhEeLHPAg28MExIkK2I3Fe++9 N4XLIitSt5uT+oxl/NrOVOpuAeUcFnNNXUXWRe116vycWXQMEuXJNYDiqakN VSBd6U1HXv5nAHMgFgvUlAQWrYHmk7ToYxYe8Mr/j+qQnzM3C7INRFc47iKZ tTS6StwgXnrggQfkoiLhL3Wh4tC06i7i5TaT3lqR/1R7MToF0MmUzSbKw9jr hUftdejBwPY0lbUTQPae1sp8HR5AArSgETACdSgbdhg4hPU7IshS8qRiCSqb uTQ859SgDa9ilBQBG3injPVEVN5GJhElttGSRx3zG7JQUERPMgRiKkjINslc IiBqf97RdBX32T+qE6EkXUd74N4SqvyK85720iTQrK7tNK5M0CjTGjQMKcV5 H37wwQcRlSSh2JmYHWhCt7qYthPPk8RmssmcHSKHVeS8jbfUD21dPRMCQLmd uKs0Jr4EvTMQXQri4tniTHfX0JTvQdK7g4Bq3tF+Hjk5Dos1HhxUlchKi2J9 sbh3WWE5RbsSbxbdQ1xR0FHEyd5xMk4DAQLQD3I4kHKMNSXG50icVIXCIzBn 4j16qyBVqlNkBjMKk512j9NlhetI+3ZSMiKYQ+t1LolKXYTuB3KC4U/edddd uMirOOG9RSITNlh5m+oe845e9OjKMVnk1sQ4aUC0ue8REyAr6sVrPeKyAZKg FOYRVZKqIxkjJOsd4vKpYAAuiECyI/gU7ypPkxa2aBMVDVjk8HnacTBJJVLV WQSIpaZXPCcrzg5zAGYDXQwWk0qjVAwa+lCUIuf8oWQQqIIAJQxYcVZ7Iwl9 mWXXCgKXnThxgr55MVKPsDuSYGNYyrfbMIwFfABtLCTlZS14BOZtclkh4zxG vEB+fdKQ0Yytc+Dj1RtOXgJGBI9AYDpt1Q7FQ+SCL2K20NYFPHhGCmzQUmVD Yloxl1Jun/FIJxz/3hVXXIGoZSyqtT7JuyxhsrzvKZrEehwtQZOeBr6jHYcV /0+7TeVR5DIGmsqBGlXM651s06ZNMVGj6YXpHPnDPIAuImt5Sbae1MITj1WU rRdrY20o2+PEyN7lTrMOSDNJOPAudvRQ5+Smb37zm+BUcBLQR322Q+UoUxtp Q+ImOZ4CASPTHGIssQS/lN+DfEoxZfQYCLmOVmN5Nyiv8udvvvlmGR6aEdjG Kwgz/blbbrkFgYvDDYn1UDhNTInjTa9L5wgEho5Gs7zEQ1eif2A2VOAkPHmU Usf6gFHkg0CGyA6QE+uRiD1gDAI0QohEq0MQ2FSCAXX5BDoZln+2OSLNiTpS iUJfeAAQ0Ca3c/GUp7HYQBuqUeAELuQIKJRDygTysUYkifYJ88FKg2SB021q L0XnwB3Gje7pJuc5dU0MivGgvJFLEnbogR+QENc4IqFoS9Km8VYqr6oBrxTx utgeOfSM99AO1IRXZ4FxaBq5K08eF9uiL70QSCWNSuWRdEU2RKtRInefJtbL ciUplER6LGQJ717ZiDhmexlLQehc7BjxsqD2qnM2fchSxKukJ7eU0ociVDAd ovbJKV9wnlc1AUEnYESCh8AqI0tfQuiWKN6ULAoQLuFJsk9hwdEtljQja+Ot D9UlwdEKbx3Gok4EakaEIUMvp/JwClA8+tpryjcli4O4DF3csrWUEgNtKoAe qZKOdlBpr/Mqp/zK6TLxudEaFc52h42bnplQCWlxkbNtGMP4iN+jMSxHv6VA S+vkYutC7riVIt5dUQB4HMGE2QZbs4NWe7d4MMG0vQNNBUMCTbetXFk5GTxd IYXodPKLdls20aYBz+jdUBP2YlrIOW6T8l1HC2ZXCb6QR7DDZ1ezI3PFEEEw swutxZPnWk9mmZv16VEjPNfBQhZ2E4oFlGGQai8J5+lDgvA6TLu42l6QwGxw 8npGnBypZSILDArNBMmio3jkCOV+JVAAOwUHSEUQ7EYDr+Okcr2qqlRPkrnU H8hMrLOzlJEtex1vQhYq5LGL+upgnYT4RpZb0+z666/fRs3tLPd3shtuuIG+ 7eBHT58+neVIRmXHLBfLBX701ltvzXJJCJsj3Jb7FRkul0E4bArlx/LXc1cf l4jgYGyWo3iWHp1G9Im6wE+If+byCGNHj+eGDiMxL9pZ7eI2AIDh89/4PCSh GH/Rt5Pchr5Q0w+3O887Uz+Mf7vL52j2BqHcoVsgvdB0/Gjt3dsc6ofSg26R it27y5naHMRU1WPEEkohRJ57V0qNyYemo34bspIVO8zmVA47e9AhMeG0acbi bNLKmK/E4j2wvbrK2kddDvGmLfZhYUG8zc6QXixSvQ66C+09I+vdt6tqqkFz yIeiWN6Cxku92Zf60PFJlTbjRYbIkGMepEkq5msCEePCOoce4ZAC60ifAAPR jiEveWvvyeaaezDlMCKwBnGvab4cntq7sb00Et4hAl51SWcPveysSEaXVILZ Vz4jk5iqMdyPQg8Mo4Ww3BTJvGGmEQT2hulEm1c1yEtfCbPswf+BtQWONb28 2FTTq6uVkY4p9mFx15HtcPFFPRXCa74dXJxeNs7Klx0W2KMK1TWRh9jG+VBg XPiSQATwir+xyqIlsFfdhA8265Au9ISUwkUIC85lsfJKEaO110V7RTU5oRaw 6WVX6SHqTqKALB4wprfq9qLRS0FrULcz4Cewp5TPFplG7Bbtkk7ExYCL+TD4 SUJNeCWbmZmRquBQM3gelAuHU0e/uNk753eA+1IQCpLjGtvrS12DoldoNLaX nO4lsS8F08HP7gLTQBXQBOalRJjZ4BLgbzAjoTD3EKnASUjeG1quCS688748 ZZOufcCZpioKJ9eGQMmSJyLo5X6LVSeJDB561QUVMQx1ZAnoOuDInZfS4hwv j/nAnoypnIG4+LzcmROHS/CcGB7ADhygLdgl2pSF9y11WaU6ozeyV/lwCDM0 CQ8fPMfaK30L+Ki9uLft6C0PH3zZH4CTUBM++DsG3gKL0zJPbgaVmjkZb+C7 NXpUqDE6CIciNprzDM+79vbcgCNawdTdjOtVxPfQE0zPDLoBi8B/gghoJx/t EhUObtg6ubkAXUskTa6DkuymWHnyIb8GoVf3aikC1/TGWvEH1O2zVWcRjRLy MBe8IYrOS2QXTCVR2XaKKZaspHpSch5EkpgxuvNCjuXejFVZI2yOxLaqq2CF cnJ9a9lcvIusDzrjszWRsndm/eLVeK+A2yikpcBdcjcqDynxLkQdXL/dO8Xh geyBUAGWydNZCljevZ9cC0pS/qru7vTqH4wBLHNCYilgefdbHk7cpEvPGM5t 6hQvAzYTsVkKbN71jqoq2VCopG5e3v3UY4DTlA5aOpz+BYd8RKx4cWG5HXRB Nmp7OSBzAMEo10NFkEu9nkJNb3yG7vALIZShe/jKtf+QLy3b2YeXD6ep3HR5 hG10Fnb2UYq3zHHmnqzCJW3HC9mUAC5mZgyAm9j4EIGXCDhnDBbvMSOXiLft wQ7exWRewLcEC6J2xoAF4yeOCwuFq7LEpTZXZakjUVIGtxWuaWzLBSz0hBfb LEGLLAvHgBZjKMaFFi5JIitT2VRCXxBmvoMJWIJTGDBCbT256kkcR8kEKXec hpAj3sqykUNZK8uxdbqNfWnJCZS8QAxB93Y6+UdeVk4JzLJGXFswl7etHACF 9fK4OHoEUCTw99wDpVMEJXHDSsJR42IsmRDFunvjAmBcvsPqiIdZpayCeDzn VdXYjO9zHoDCrvsqMJOkm6xV/3AEUMbo8Juw0cqD4q5Szcn+MaxSzS7lKsFk 5jAGmAyOVsosupAYDhlD2MSkfR2Ms58odUxUcQwwmajnyqtpFya3EOkyYDIh z1WCyYQSxwCTsUGrBJOpfTEGmEwq+qRhqqqTuhSY1FnjmNGylG2wMcBkoh1j gMkknWmY+C6UlJ7N/8PFwk1ztMcAr1sWdcnwksOqEywkCR1XJuOG0Cx7srW4 eGGCABrBGwNBjeBpgnIWfA7qImX1xTlQDtHj1FebUmtRPAbfDnAbJvMQp+He xm3nz5/HgT6q9XAQbV07ChB5ju/3nQQ2vSr7XlrKiP0Z8dCMw7HOeHcxvnlJ sisEpzHtY4DTuB4aTlEDiWtyurbK72seN7zM1SZl97Llw2v6G4a3w3CG/+WK oWl25hhANaRYDqhth1WOKVAnYbUYJuNdjwEmsyE0aZjYwzG7cfuXAZN3O9Dx 1YHJhNbGAJPx/lcJJpMAOgaYjIEdgmlF1rzeTtfCSJAMBWskpet4nPRkiOLd 5rFMUCQU0Oz45vhA4V1gE4mYHwkUk1kaUXLV6sBhSLIMOGIYctJweIWE9y0N DqMzVgkOQ4/R4IjbYKIxrp4wCDMOKfaODoIuQLlKIBhOGA0EaovpXddMeP6c tGEOcewZcf5y9ym+Dc1/RSwdz9qUdG46a04+pBTRLVu2kHt/dZzpZLDO9XjM snH3aPOHzGbirz9vwlzDNYbMVtKuhvPnnNZYTRO/nVid+ZvqrE3nz6nGujox fp80COpS4VjMdOdoIGjyrQYIPIaJ1e9oCAInupuVA4MxcW7iQ6vGjG1vCAfH 0E0NZS81eNJ2YaMjIU1h4vCkqd/t5f5OGibOMzNXpWwbDSYTUWgE04rYPy4L A48h1mPZ2hASDgPWXt19VZz9ZKjDp+5NJH9LQ5goVjSlLjdfVT7j/kyUZHND SDjyYPTaWoBpa8RwU0gWHHp6eeGThoQl3lQFu7QhTPOhTQrIlF4LN2mYWDOb ylkbG8K0L3FrrCHplUKaNEw7kiREXHvjeu/uTfiIEjizBmDaaWVCFXyth4nX N3LCkL55VwJPGib2m+UWDlXXsx4mfleK3JIzur4JTCtiVXkFJkfQyM57lVMr KCsn9cgP9HL5J21VuVR6sdLOTEOwVCHbeF3BBue5STOdUlrxmLVX2rHcAHSE 6WAI6DqHNQDWvOVA8s+8CoilFjrIYTyI5xXVnDRMCwndj3DVbcDXtGRdcDfo hNkTXOUAf29eA5BxtRqR+ULd+VrImDgIvUWCebUdJw3W5c68mt4au9FHiVcv UYMl2+kpi4hSByaxnR4W4D3hr6i8MzT7VZ6CXqezkFItARztVSLUQKaclwut xXNUxXIiuREc7cEUMV2EqcsgC9auI7ZcLhXY1RgyZEOEamApZ+h36L0uSvvh K+oDZnL/y7X0W4eqBCKnCOX5cgfnVWguVOe7EKpgIbWscamx5aAttPUwBswe zIQUjATlvaI3nAUDaySPx1KeXg0wF4cXFluLT2YRo5MB09zTITc2SEEeSDFZ fq84VSfJjxQGxwd/o22hkeR3c7AvkEAQC52cnHAcc8bgomcAWa6bgFqA1w3K hyX6jOkmC0d4pRI3KC9F7r2CTp7uA+FPivabTCqRBzyf9AXPAmBdVxzQqePu fWQAymVIQ3XEvcJMk7ZkHnzB853BHDFXzBm/FwtG4/fkFc7Izcu6wiYpMyY0 XYeXP+7V71kLUAeTNSjWLgcYoByeBc92nAPLoW1Gyh4KM4DoeAarJMlzwTQ8 s7IWMMCeFXSSnEtu8WFkfT45mPUQkC/IRP24awHOpu9WL6m9HlenrWPaskaL 57Uze8ewVPCHlxnsrZMutl1s+21quygXF9sutjWQiwx/3xIaWpf8fyRdWQU= \ \>"]] }, Open ]], Cell[CellGroupData[{ Cell["\<\ plot3=Graphics3D[{ Sphere[{xc[[1]],yc[[1]],zc[[1]]},Cr], Sphere[{xc[[2]],yc[[2]],zc[[2]]},Nr], Sphere[{xc[[3]],yc[[3]],zc[[3]]},Oa], {Opacity[0.6], Sphere[{xc[[4]],yc[[4]],zc[[4]]},NH], Sphere[{xc[[5]],yc[[5]],zc[[5]]},NH], Sphere[{xc[[6]],yc[[6]],zc[[6]]},CH], Sphere[{xc[[7]],yc[[7]],zc[[7]]},CN], Sphere[{xc[[8]],yc[[8]],zc[[8]]},CO], {Lighter[Blue], Sphere[{xc[[9]],yc[[9]],zc[[9]]},OL], Sphere[{xc[[10]],yc[[10]],zc[[10]]},NM], Sphere[{xc[[11]],yc[[11]],zc[[11]]},NM], Sphere[{xc[[14]],yc[[14]],zc[[14]]},OM], Sphere[{xc[[15]],yc[[15]],zc[[15]]},OM]}}, {Opacity[0.4],{Lighter[Yellow], Sphere[{xc[[12]],yc[[12]],zc[[12]]},OC], Sphere[{xc[[13]],yc[[13]],zc[[13]]},OC]}}, {Darker[Red,1], Sphere[{xn[[1]],yc[[1]],zn[[1]]},0.08], Sphere[{xn[[2]],yc[[2]],zn[[2]]},0.08], Sphere[{xn[[3]],yc[[3]],zn[[3]]},0.08], Sphere[{xn[[4]],yc[[4]],zn[[4]]},0.08], Sphere[{xn[[5]],yc[[5]],zn[[5]]},0.08], Sphere[{xn[[6]],yc[[6]],zn[[6]]},0.08]}}] /. t[[2]]; Show[plot3,{AspectRatio \[Rule] Automatic,Boxed->False}]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5668391703953342`*^9, 3.566839177774147*^9}, { 3.5669295568515253`*^9, 3.5669295979107976`*^9}, {3.5669296458028817`*^9, 3.5669298487904377`*^9}, {3.5669298934845166`*^9, 3.5669299137177525`*^9}, {3.566929958552231*^9, 3.566929984853877*^9}, { 3.5669300180975356`*^9, 3.5669300244935465`*^9}, {3.5669681367168245`*^9, 3.566968174578091*^9}, {3.5669682084145503`*^9, 3.5669682238429775`*^9}, { 3.56696839045127*^9, 3.566968403352493*^9}, {3.5669886931259174`*^9, 3.5669887406436005`*^9}, {3.56698877437086*^9, 3.5669887786296673`*^9}, { 3.5669888178481364`*^9, 3.5669888516221952`*^9}, {3.5669888927126675`*^9, 3.5669888937422695`*^9}, {3.56698893940355*^9, 3.566988940932352*^9}, { 3.566989011850077*^9, 3.5669890235968976`*^9}, {3.5669890790549955`*^9, 3.566989179831172*^9}, {3.5669892163196363`*^9, 3.5669892198764424`*^9}, { 3.566989533171793*^9, 3.5669895449498134`*^9}, 3.567068931415569*^9, { 3.5670704486430335`*^9, 3.567070500840725*^9}, {3.5670768428602643`*^9, 3.567076863427901*^9}, {3.567076934439225*^9, 3.5670771373487816`*^9}, { 3.567077191402877*^9, 3.567077229716544*^9}, {3.5670776494196815`*^9, 3.5670776603397007`*^9}, 3.567077706489382*^9, 3.5670779441365986`*^9, { 3.5670781500413604`*^9, 3.567078154253368*^9}, {3.567078187263026*^9, 3.5670782620127573`*^9}, {3.5670783597557287`*^9, 3.567078441453072*^9}, { 3.5670784741819296`*^9, 3.5670785488436604`*^9}, {3.567079302263584*^9, 3.567079352339672*^9}, 3.5670794136233797`*^9, {3.567084448847624*^9, 3.567084513977738*^9}, {3.567084586159065*^9, 3.5670845877970676`*^9}, { 3.567084647170772*^9, 3.5670846660780053`*^9}, {3.5670850221714306`*^9, 3.567085027023039*^9}}], Cell[BoxData[ Graphics3DBox[{SphereBox[{0, 0, 0}, 0.2597231856403463], SphereBox[{-2.5449706837000776`, 0, 0}, 0.22405135815353058`], SphereBox[{1.291429196139679, 1.884761711991131, 0}, 0.1916417225081855], {Opacity[0.6], SphereBox[{-3.200397178358574, 1.2212124156850317`, 0}, 1.1619297579240278`], SphereBox[{-3.200397178358574, -1.2212124156850317`, 0}, 1.1619297579240278`], SphereBox[{0.6081672572221739, -1.3853495355430014`, 0}, 1.2532410401053755`], SphereBox[{-1.2903212555934467`, 0, 0}, 1.0305980699531003`], SphereBox[{0.6649556762928184, 0.9704620296599826, 0}, 0.9166962134031293], {RGBColor[ NCache[ Rational[1, 3], 0.3333333333333333], NCache[ Rational[1, 3], 0.3333333333333333], 1], SphereBox[{1.9046760746972449`, 2.77975792251382, 0}, 0.8932960411034288], SphereBox[{-2.5449706837000776`, 0, 0.7472908218886891}, 0.7472908218886891], SphereBox[{-2.5449706837000776`, 0, -0.7472908218886891}, 0.7472908218886891], SphereBox[{0.44483311166516193`, 2.4648450688802646`, 0}, 0.8346241451827114], SphereBox[{2.138025280614196, 1.3046783551019976`, 0}, 0.8346241451827114]}}, {Opacity[0.4], {RGBColor[1, 1, NCache[ Rational[1, 3], 0.3333333333333333]], SphereBox[{0.6649556762928184, 0.9704620296599826, 1.5679221435272637`}, 1.5679221435272637`], SphereBox[{0.6649556762928184, 0.9704620296599826, -1.5679221435272637`}, 1.5679221435272637`]}}, {RGBColor[0, 0, 0], SphereBox[{0, 0, 0}, 0.08], SphereBox[{-2.5449706837000776`, 0, 0}, 0.08], SphereBox[{1.291429196139679, 1.884761711991131, 0}, 0.08], SphereBox[{-3.4455458033171533`, 1.2212124156850317`, 0}, 0.08], SphereBox[{-3.4455458033171533`, -1.2212124156850317`, 0}, 0.08], SphereBox[{0.8381266667787244, -1.3853495355430014`, 0}, 0.08]}}, AspectRatio->Automatic, AutomaticImageSize->True, Boxed->False, ImageSize->{361.1608400917048, 323.7429743739751}, ViewPoint->{-0.22029197196915629`, 0.1662356257434774, 3.3725119960972134`}, ViewVertical->{-0.16726004335463357`, 0.9574806600091574, 0.6332968355980947}]], "Output", CellChangeTimes->{3.567230369599126*^9, 3.567244942785922*^9, 3.5672451961303673`*^9, 3.5672453996951246`*^9, 3.5672454375407915`*^9, 3.567245473810855*^9, 3.5672456741152067`*^9, 3.567245720525288*^9, 3.567245760352158*^9, 3.5672458110678473`*^9, 3.5672459741193333`*^9, 3.567246007565792*^9, 3.5672462444678087`*^9, 3.5672463166959352`*^9, 3.5672463764908404`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"as", ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"x", "=", RowBox[{ RowBox[{"bohr", "*", "xn"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "=", RowBox[{ RowBox[{"bohr", "*", "yn"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"z", "=", RowBox[{ RowBox[{"bohr", "*", "zn"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"s", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"as", "[", RowBox[{"[", "l", "]"}], "]"}], ",", RowBox[{"x", "[", RowBox[{"[", "l", "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", "l", "]"}], "]"}], ",", RowBox[{"z", "[", RowBox[{"[", "l", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"l", ",", "6"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ "Cartesian", " ", "Coordinates", " ", "of", " ", "stationary", " ", "structure"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{"u", "=", RowBox[{"MatrixForm", "[", "s", "]"}]}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.5672284907474256`*^9, 3.56722854442712*^9}, { 3.5672285770935774`*^9, 3.567228591570403*^9}, 3.567228630320871*^9, { 3.56722880051717*^9, 3.567228812825591*^9}, {3.5672288915745296`*^9, 3.5672289256605897`*^9}, {3.5672289905567036`*^9, 3.5672290107431393`*^9}, {3.567229046592002*^9, 3.5672290563576193`*^9}, { 3.5672291676482143`*^9, 3.567229187179449*^9}, {3.5672292331995296`*^9, 3.567229305692857*^9}, {3.5672293946910133`*^9, 3.567229403739029*^9}, { 3.567229509632015*^9, 3.567229553296492*^9}, {3.567229926074747*^9, 3.5672299410507727`*^9}, {3.567229996883271*^9, 3.5672300622005854`*^9}, { 3.56723009870465*^9, 3.567230152836745*^9}, {3.567230250492916*^9, 3.567230251787719*^9}, {3.56723029233219*^9, 3.5672303016454067`*^9}, { 3.5672449305087004`*^9, 3.567244935375909*^9}, {3.5672450518145137`*^9, 3.567245059255727*^9}, {3.5672451122178197`*^9, 3.5672451847423472`*^9}, { 3.5672453557654476`*^9, 3.5672453904911084`*^9}, 3.567245430770379*^9, { 3.5672454639672375`*^9, 3.567245468038845*^9}, {3.5672455930730643`*^9, 3.567245597597072*^9}, {3.5672457069220643`*^9, 3.567245710744071*^9}, { 3.5672457475133357`*^9, 3.5672457532229457`*^9}, {3.5672457998514276`*^9, 3.5672458038294344`*^9}, {3.5672459644785166`*^9, 3.5672459999217787`*^9}, {3.567246112413576*^9, 3.567246114847181*^9}, { 3.5672461687920756`*^9, 3.567246233594589*^9}, {3.567246288413086*^9, 3.5672462895674877`*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"C", "0.`", "0.`", "0.`"}, {"N", RowBox[{"-", "1.346739951488356`"}], "0.`", "0.`"}, {"O", "0.683394627725607`", "0.9973725484663308`", "0.`"}, {"H", RowBox[{"-", "1.8233035915619613`"}], "0.8879492038617962`", "0.`"}, {"H", RowBox[{"-", "1.8233035915619613`"}], RowBox[{"-", "0.8879492038617962`"}], "0.`"}, {"H", "0.443517355145965`", RowBox[{"-", "1.0102920778128353`"}], "0.`"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{3.5672457205720882`*^9, 3.567245760383358*^9, 3.567245811114647*^9, 3.567245974166134*^9, 3.5672460076125927`*^9, 3.5672462445146084`*^9, 3.567246316742735*^9, 3.5672463765220404`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"s", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]], "Input", CellChangeTimes->{{3.567246294918297*^9, 3.567246308911522*^9}, { 3.567246367739225*^9, 3.5672463679576254`*^9}}], Cell[BoxData["0.683394627725607`"], "Output", CellChangeTimes->{3.5672463167895355`*^9, 3.5672463765376406`*^9}] }, Open ]] }, WindowToolbars->"EditBar", WindowSize->{830, 964}, WindowMargins->{{Automatic, 525}, {Automatic, 0}}, PrintingCopies->1, PrintingPageRange->{5, 6}, PrintingOptions->{"Magnification"->1., "PaperOrientation"->"Portrait", "PaperSize"->{595.2, 841.8}}, PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}}, 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, 2599, 38, 48, "Input"], Cell[3159, 60, 3877, 59, 218, "Input"], Cell[7039, 121, 1319, 26, 72, "Input"], Cell[8361, 149, 758, 23, 252, "Input"], Cell[9122, 174, 4686, 67, 82, "Input"], Cell[13811, 243, 4950, 71, 167, "Input"], Cell[18764, 316, 5053, 71, 167, "Input"], Cell[23820, 389, 4639, 66, 82, "Input"], Cell[28462, 457, 4771, 67, 82, "Input"], Cell[33236, 526, 4617, 66, 82, "Input"], Cell[37856, 594, 489, 12, 99, "Input"], Cell[38348, 608, 384, 11, 99, "Input"], Cell[38735, 621, 571, 14, 116, "Input"], Cell[CellGroupData[{ Cell[39331, 639, 1756, 26, 116, "Input"], Cell[41090, 667, 1497, 28, 72, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[42624, 700, 955, 24, 235, "Input"], Cell[43582, 726, 422, 6, 31, "Output"], Cell[44007, 734, 438, 7, 31, "Output"], Cell[44448, 743, 422, 6, 31, "Output"], Cell[44873, 751, 420, 6, 31, "Output"], Cell[45296, 759, 422, 6, 31, "Output"], Cell[45721, 767, 422, 6, 31, "Output"], Cell[46146, 775, 422, 6, 31, "Output"], Cell[46571, 783, 421, 6, 31, "Output"], Cell[46995, 791, 422, 6, 31, "Output"], Cell[47420, 799, 421, 6, 31, "Output"], Cell[47844, 807, 422, 6, 31, "Output"], Cell[48269, 815, 422, 6, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[48728, 826, 1572, 26, 218, "Input"], Cell[50303, 854, 1917, 32, 328, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[52257, 891, 1466, 25, 218, "Input"], Cell[53726, 918, 22969, 381, 327, 1364, 26, "CachedBoxData", "BoxData", \ "Output"] }, Open ]], Cell[CellGroupData[{ Cell[76732, 1304, 2803, 54, 490, "Input"], Cell[79538, 1360, 2631, 55, 339, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[82206, 1420, 2903, 64, 172, "Input"], Cell[85112, 1486, 1210, 29, 116, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[86359, 1520, 218, 5, 31, "Input"], Cell[86580, 1527, 113, 1, 31, "Output"] }, Open ]] } ] *) (* End of internal cache information *)