(* 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[ 131499, 2473] NotebookOptionsPosition[ 129090, 2397] NotebookOutlinePosition[ 129695, 2419] CellTagsIndexPosition[ 129652, 2416] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell["\<\ (*CH2CH2 Ethene, with \"orthogonal\" \[Sigma],\[Pi] LPs, 22.01.2013 *)\ \>", "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}, {3.5672691891313963`*^9, 3.567269197539811*^9}, 3.567334117644643*^9, {3.567762620897001*^9, 3.567762631193019*^9}, { 3.5677650115095997`*^9, 3.567765035986043*^9}, {3.5678408184988356`*^9, 3.5678408188888364`*^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,Cpi,NH,CH,CO,CN,CC,NM,OC,OO,OM,w1,w2,w3,pch,pnh]; ZC=6.0; sig1=0.3; sig4=0.5; sig5=0.3; nc=11; nn=6; pi=0.0; k1=1.02; k4=1.058; k5=1.12; k11=1.12; bohr=0.529177; rad=57.29578; (* 6 Variables: Cr,CH,CC,Cpi,pch,w1 *)\ \>", "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}, {3.567247634991851*^9, 3.567247676659524*^9}, { 3.56724778291131*^9, 3.5672477892761216`*^9}, {3.5672489000760727`*^9, 3.5672489007000737`*^9}, {3.5672492510454893`*^9, 3.567249259407104*^9}, { 3.567250434479168*^9, 3.5672504428095827`*^9}, {3.567250921215623*^9, 3.567250922058024*^9}, {3.567251110053954*^9, 3.5672511110211563`*^9}, 3.5672511638272486`*^9, {3.56725131811152*^9, 3.5672513271595354`*^9}, { 3.5672514083888783`*^9, 3.5672514174368944`*^9}, 3.5672514476697474`*^9, { 3.567251504469447*^9, 3.5672515050778484`*^9}, {3.5672534213240137`*^9, 3.567253422462816*^9}, {3.567253494862543*^9, 3.567253499308551*^9}, { 3.567268602539166*^9, 3.5672686237708035`*^9}, {3.5672686825673065`*^9, 3.5672686838621087`*^9}, {3.567268721052574*^9, 3.567268760146243*^9}, { 3.5672688031555185`*^9, 3.567268961854597*^9}, {3.5672690063302755`*^9, 3.5672690070946765`*^9}, {3.567269054035159*^9, 3.5672690547839603`*^9}, { 3.5673365018124304`*^9, 3.567336653163896*^9}, {3.567336686672755*^9, 3.56733668948076*^9}, {3.5673367723637056`*^9, 3.5673367928465414`*^9}, { 3.5673447921413593`*^9, 3.567344794621764*^9}, 3.5673448590030775`*^9, { 3.5677631198970776`*^9, 3.5677631210046797`*^9}, {3.5677635845750937`*^9, 3.56776369079568*^9}, {3.5677648244964714`*^9, 3.567764830518082*^9}, { 3.5677672065270557`*^9, 3.567767225808689*^9}, {3.567767298146016*^9, 3.567767316897249*^9}, {3.567767358440122*^9, 3.5677673642277327`*^9}, { 3.567767899159472*^9, 3.5677678998146734`*^9}, {3.5677731944066515`*^9, 3.567773206496673*^9}, {3.567773268359581*^9, 3.567773269295583*^9}, 3.567775337129215*^9, 3.567775412461747*^9, 3.567775449028211*^9, { 3.5677754839566727`*^9, 3.567775484206273*^9}, {3.567775531848757*^9, 3.567775553017994*^9}, {3.5677757367395163`*^9, 3.5677757538215466`*^9}, { 3.5677759866831555`*^9, 3.5677760378824453`*^9}, {3.5677762298719826`*^9, 3.567776244270808*^9}, {3.567776758837712*^9, 3.567776760491315*^9}, { 3.5677768013165865`*^9, 3.567776808367799*^9}, {3.567776846790666*^9, 3.5677768725151114`*^9}, 3.5677769237768016`*^9, {3.567776962215269*^9, 3.5677769890317163`*^9}, {3.5678310873237467`*^9, 3.5678310880257483`*^9}, {3.567831193731534*^9, 3.567831194527135*^9}, { 3.5678315027368765`*^9, 3.5678315306297255`*^9}, {3.5678321861116767`*^9, 3.5678321912752857`*^9}, 3.567832518485861*^9, {3.5678326470612864`*^9, 3.5678327286806297`*^9}, 3.567832782173124*^9, {3.567832936550995*^9, 3.567832938173398*^9}, {3.5678333274876814`*^9, 3.5678333279244823`*^9}, 3.5678333613085413`*^9, {3.567833468262329*^9, 3.567833544780463*^9}, { 3.567833664214273*^9, 3.567833693479924*^9}, 3.5678337663632526`*^9, { 3.56783380486412*^9, 3.567833837858178*^9}, {3.5678406463057337`*^9, 3.5678407057574377`*^9}, {3.567840841602476*^9, 3.5678408427880783`*^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.0*2.25*k1/Cr^2; vee=2.0*3.0*sig1/Cr; vne=-2.0*3.0*ZC/Cr;\ \>", "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}, {3.567763734928158*^9, 3.5677637671422143`*^9}, 3.5677638284971223`*^9, {3.567765306490518*^9, 3.5677653068649187`*^9}, {3.567775854238923*^9, 3.5677758995414023`*^9}}], Cell["\<\ (* bonding and lone pairs *) Ekin = Ekin + 2.25*(4*k4/CH^2+k5/CC^2+k11/Cpi^2); vee=vee+3.0*(4*sig4/CH+sig5/CC);\ \>", "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.567247718249197*^9, 3.567247775142497*^9}, { 3.5677637786862345`*^9, 3.5677637898246546`*^9}, {3.567763845891153*^9, 3.5677638701647954`*^9}, {3.567765140677827*^9, 3.5677651535010495`*^9}, { 3.5677651939987206`*^9, 3.5677651954027233`*^9}}], Cell["\<\ (* 1 2 3 4 5 6 7 8 9 10 11 *) (* C C Hc Hc Hc Hc CC L1 L1 L1 L1 *) as={\" \",C,C,H,H,H,H}; (* atom string *) rr={Cr,Cr,CH,CH,CH,CH,CC, Cpi ,Cpi, Cpi, Cpi}; oc={-2,-2,-2,-2,-2,-2,-2,-1/2,-1/2,-1/2,-1/2}; ch={ZC,ZC,1, 1, 1, 1}; csw1=Cos[w1/2]; ssw1=Sin[w1/2];\ \>", "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}, { 3.567247444921117*^9, 3.5672474899115963`*^9}, 3.5672475334980726`*^9, { 3.567247593043377*^9, 3.5672476082534037`*^9}, {3.5672478183077726`*^9, 3.567247880489482*^9}, {3.5672480027468967`*^9, 3.567248042745367*^9}, { 3.5672487779746585`*^9, 3.567248812669119*^9}, {3.5673358974985685`*^9, 3.5673358994953723`*^9}, {3.5677631586319456`*^9, 3.5677631590999465`*^9}, {3.5677632093788347`*^9, 3.5677633189534273`*^9}, {3.5677633719311204`*^9, 3.5677633964231634`*^9}, {3.567763882239217*^9, 3.567763944405326*^9}, { 3.567776064511692*^9, 3.5677761884383097`*^9}, {3.567832567547947*^9, 3.567832627514452*^9}}], Cell["\<\ dccx=Cr+CC; dchx=dccx+(Cr+CH)*csw1; dchy=(Cr+CH)*ssw1; pchx=dccx+(Cr+CH+pch)*csw1; pchy=(Cr+CH+pch)*ssw1; dcpi=0.5*Cr+Cpi; dpix=dccx-0.15; \ \>", "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}, {3.5673117165829782`*^9, 3.5673117180649805`*^9}, {3.5677634771845055`*^9, 3.5677635233137865`*^9}, {3.567763966167364*^9, 3.5677640107054424`*^9}, { 3.567764041234696*^9, 3.567764216859804*^9}, {3.5677642599470797`*^9, 3.567764312675173*^9}, {3.567831018153226*^9, 3.5678310302120466`*^9}, { 3.567831122595409*^9, 3.567831125450214*^9}, 3.5678315904402304`*^9, { 3.5678316279738965`*^9, 3.5678316535267415`*^9}, {3.567831973405303*^9, 3.5678319885061297`*^9}, {3.567832056210249*^9, 3.5678320591898537`*^9}, { 3.5678322168125305`*^9, 3.5678322234737425`*^9}, {3.5678323546231728`*^9, 3.567832357945979*^9}, {3.5678329016537333`*^9, 3.5678329050389395`*^9}, 3.567833006563918*^9, {3.5678330685116267`*^9, 3.567833091568467*^9}, { 3.567833132643339*^9, 3.56783313292414*^9}, {3.5678332038574643`*^9, 3.5678332045438657`*^9}, {3.5678333832733793`*^9, 3.567833383990981*^9}, { 3.5678335651384993`*^9, 3.5678335950437517`*^9}, {3.567833631173415*^9, 3.5678336318910165`*^9}}], Cell["\<\ (* 6 nuclei: coord *) xn={dccx,-dccx,pchx,pchx,-pchx,-pchx}; yn={0,0,pchy,-pchy,pchy,-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}, {3.567764337884817*^9, 3.567764477364662*^9}}], Cell["\<\ (* 11 clouds: coord *) xc={dccx,-dccx,dchx,dchx,-dchx,-dchx,0,dpix,dpix,-dpix,-dpix}; yc={0,0,dchy,-dchy,dchy,-dchy,0,0,0,0,0}; zc={0,0,0,0,0,0,0,dcpi,-dcpi,dcpi,-dcpi};\ \>", "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}, {3.5672479122355375`*^9, 3.567247913155939*^9}, {3.5672481119782887`*^9, 3.567248113819092*^9}, { 3.567248190337226*^9, 3.567248240850115*^9}, {3.5677645199995365`*^9, 3.5677647291647043`*^9}, {3.5677652477096148`*^9, 3.56776525051762*^9}, { 3.567767048701578*^9, 3.5677670726632204`*^9}, {3.5678310393536625`*^9, 3.5678310528476863`*^9}, {3.56783125443124*^9, 3.5678312778624816`*^9}, { 3.5678313501062083`*^9, 3.567831355238618*^9}, {3.5678314371543617`*^9, 3.5678314455003757`*^9}, {3.5678318942819643`*^9, 3.567831919585209*^9}}], Cell["\<\ (* epot terms *) vne=vne-4*(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}, {3.56776474933554*^9, 3.567764757104353*^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]]] vee=vee-0.25/Cpi;\ \>", "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}, {3.5677623897825947`*^9, 3.5677624289854636`*^9}, {3.567764779147192*^9, 3.567764795542821*^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},{CH,1.25},{CC,0.99},{Cpi,1.28},{w1,2.04},{pch,\ 0.55},{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}, { 3.5672486648744593`*^9, 3.56724869903852*^9}, {3.5677648551505256`*^9, 3.567764927175852*^9}, {3.567773235552724*^9, 3.5677732525567536`*^9}, { 3.567773333789296*^9, 3.5677733541317315`*^9}, {3.567775588461256*^9, 3.567775614716102*^9}, {3.567831732884081*^9, 3.567831752368515*^9}, { 3.5678320708742743`*^9, 3.5678320800470905`*^9}, {3.567832243051777*^9, 3.567832256826601*^9}, 3.567832371330802*^9, 3.5678328367888193`*^9, 3.567832892792918*^9, 3.567832982383876*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "78.13664027771551`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Cr", "\[Rule]", "0.26062448315986025`"}], ",", RowBox[{"CH", "\[Rule]", "1.2409871716713938`"}], ",", RowBox[{"CC", "\[Rule]", "0.9871003356287535`"}], ",", RowBox[{"Cpi", "\[Rule]", "1.2072995482981752`"}], ",", RowBox[{"w1", "\[Rule]", "2.0272652301267966`"}], ",", RowBox[{"pch", "\[Rule]", "0.5414621448498486`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.567773937931557*^9, 3.567774337011458*^9, 3.567774376370327*^9, 3.5677745394062138`*^9, 3.5677746263763666`*^9, 3.56777469668569*^9, 3.56777489628804*^9, 3.56777508398757*^9, 3.5677751633917093`*^9, 3.5677752738243036`*^9, 3.567775307161562*^9, 3.567775347628033*^9, 3.567775421681363*^9, 3.567775454503821*^9, 3.567775490103083*^9, { 3.5677755367939653`*^9, 3.567775561613609*^9}, 3.567775624294519*^9, 3.5677757593907566`*^9, 3.5677759371998687`*^9, 3.567776200387931*^9, 3.567776602931038*^9, {3.567776644458311*^9, 3.567776690712392*^9}, 3.5677767648437223`*^9, 3.567776812845007*^9, {3.5677768525158763`*^9, 3.5677768787863226`*^9}, 3.567776927973209*^9, {3.567776967472478*^9, 3.567776993618124*^9}, 3.5677848690954437`*^9, {3.567831064734907*^9, 3.5678310933921576`*^9}, 3.5678311344670296`*^9, 3.5678311996595445`*^9, 3.567831289032101*^9, 3.567831361837429*^9, 3.567831450710785*^9, { 3.567831508274886*^9, 3.5678315345453324`*^9}, 3.5678315954166393`*^9, { 3.5678316411715193`*^9, 3.5678316582379494`*^9}, 3.5678317687485437`*^9, 3.567832019893385*^9, 3.567832100061926*^9, 3.567832272910229*^9, 3.567832379162016*^9, 3.5678325254122725`*^9, {3.567832655329301*^9, 3.5678327336414385`*^9}, 3.5678327871183324`*^9, 3.5678328415468283`*^9, { 3.567832914601756*^9, 3.567832943446207*^9}, 3.5678330109631257`*^9, { 3.567833075094838*^9, 3.567833097278077*^9}, 3.5678331386025496`*^9, 3.567833209645075*^9, 3.56783333216769*^9, {3.5678333661601496`*^9, 3.56783338918579*^9}, {3.5678334726459365`*^9, 3.5678336000981607`*^9}, 3.5678336366022243`*^9, {3.5678336691282816`*^9, 3.567833697941532*^9}, 3.567833771511262*^9, {3.5678338174533424`*^9, 3.5678338420701857`*^9}, { 3.5678406904538107`*^9, 3.5678407102658453`*^9}, 3.567840850681692*^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]] (2*Cr+2*CC)*bohr /. t[[2]] (Cr+CH+pch)*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}, 3.567345843318006*^9, { 3.5677649519798956`*^9, 3.5677649783907423`*^9}, {3.567765362744217*^9, 3.567765377377043*^9}}], Cell[BoxData["78.07252628961751`"], "Output", CellChangeTimes->{ 3.567773937947157*^9, 3.5677743370270576`*^9, 3.567774376370327*^9, 3.5677745394062138`*^9, 3.5677746263919663`*^9, 3.56777469670129*^9, 3.5677748963036404`*^9, 3.56777508398757*^9, 3.5677751634073095`*^9, 3.567775273839904*^9, 3.567775307177162*^9, 3.567775347643633*^9, 3.567775421681363*^9, 3.5677754545194206`*^9, 3.5677754901186833`*^9, { 3.5677755368095655`*^9, 3.567775561629209*^9}, 3.567775624310119*^9, 3.5677757594063563`*^9, 3.5677759372154684`*^9, 3.567776200387931*^9, 3.567776602946638*^9, {3.567776644458311*^9, 3.567776690712392*^9}, 3.5677767648593225`*^9, 3.5677768128606067`*^9, {3.5677768525314765`*^9, 3.5677768788019223`*^9}, 3.5677769279888086`*^9, {3.567776967488078*^9, 3.567776993633724*^9}, 3.5677848690954437`*^9, {3.5678310647505074`*^9, 3.5678310934077578`*^9}, 3.5678311344670296`*^9, 3.567831199675144*^9, 3.567831289047701*^9, 3.567831361837429*^9, 3.567831450710785*^9, { 3.5678315082904863`*^9, 3.5678315345453324`*^9}, 3.5678315954322395`*^9, { 3.5678316411715193`*^9, 3.5678316582535496`*^9}, 3.567831768764144*^9, 3.567832019893385*^9, 3.567832100077526*^9, 3.5678322729258294`*^9, 3.567832379177616*^9, 3.5678325254278727`*^9, {3.567832655329301*^9, 3.5678327336414385`*^9}, 3.567832787133932*^9, 3.567832841562428*^9, { 3.567832914601756*^9, 3.5678329434618073`*^9}, 3.5678330109631257`*^9, { 3.5678330751104383`*^9, 3.5678330972936773`*^9}, 3.5678331386025496`*^9, 3.5678332096606746`*^9, 3.5678333321832895`*^9, {3.56783336617575*^9, 3.5678333892013903`*^9}, {3.5678334726615367`*^9, 3.5678336001137605`*^9}, 3.5678336366178246`*^9, {3.5678336691282816`*^9, 3.5678336979571323`*^9}, 3.567833771511262*^9, {3.567833817468942*^9, 3.5678338420857854`*^9}, { 3.5678406904538107`*^9, 3.5678407102814455`*^9}, 3.5678408506972923`*^9}], Cell[BoxData[ RowBox[{"-", "252.20657086708005`"}]], "Output", CellChangeTimes->{ 3.567773937947157*^9, 3.5677743370270576`*^9, 3.567774376370327*^9, 3.5677745394062138`*^9, 3.5677746263919663`*^9, 3.56777469670129*^9, 3.5677748963036404`*^9, 3.56777508398757*^9, 3.5677751634073095`*^9, 3.567775273839904*^9, 3.567775307177162*^9, 3.567775347643633*^9, 3.567775421681363*^9, 3.5677754545194206`*^9, 3.5677754901186833`*^9, { 3.5677755368095655`*^9, 3.567775561629209*^9}, 3.567775624310119*^9, 3.5677757594063563`*^9, 3.5677759372154684`*^9, 3.567776200387931*^9, 3.567776602946638*^9, {3.567776644458311*^9, 3.567776690712392*^9}, 3.5677767648593225`*^9, 3.5677768128606067`*^9, {3.5677768525314765`*^9, 3.5677768788019223`*^9}, 3.5677769279888086`*^9, {3.567776967488078*^9, 3.567776993633724*^9}, 3.5677848690954437`*^9, {3.5678310647505074`*^9, 3.5678310934077578`*^9}, 3.5678311344670296`*^9, 3.567831199675144*^9, 3.567831289047701*^9, 3.567831361837429*^9, 3.567831450710785*^9, { 3.5678315082904863`*^9, 3.5678315345453324`*^9}, 3.5678315954322395`*^9, { 3.5678316411715193`*^9, 3.5678316582535496`*^9}, 3.567831768764144*^9, 3.567832019893385*^9, 3.567832100077526*^9, 3.5678322729258294`*^9, 3.567832379177616*^9, 3.5678325254278727`*^9, {3.567832655329301*^9, 3.5678327336414385`*^9}, 3.567832787133932*^9, 3.567832841562428*^9, { 3.567832914601756*^9, 3.5678329434618073`*^9}, 3.5678330109631257`*^9, { 3.5678330751104383`*^9, 3.5678330972936773`*^9}, 3.5678331386025496`*^9, 3.5678332096606746`*^9, 3.5678333321832895`*^9, {3.56783336617575*^9, 3.5678333892013903`*^9}, {3.5678334726615367`*^9, 3.5678336001137605`*^9}, 3.5678336366178246`*^9, {3.5678336691282816`*^9, 3.5678336979571323`*^9}, 3.567833771511262*^9, {3.567833817468942*^9, 3.5678338420857854`*^9}, { 3.5678406904538107`*^9, 3.5678407102814455`*^9}, 3.5678408506972923`*^9}], Cell[BoxData["62.434301682186515`"], "Output", CellChangeTimes->{ 3.567773937947157*^9, 3.5677743370270576`*^9, 3.567774376370327*^9, 3.5677745394062138`*^9, 3.5677746263919663`*^9, 3.56777469670129*^9, 3.5677748963036404`*^9, 3.56777508398757*^9, 3.5677751634073095`*^9, 3.567775273839904*^9, 3.567775307177162*^9, 3.567775347643633*^9, 3.567775421681363*^9, 3.5677754545194206`*^9, 3.5677754901186833`*^9, { 3.5677755368095655`*^9, 3.567775561629209*^9}, 3.567775624310119*^9, 3.5677757594063563`*^9, 3.5677759372154684`*^9, 3.567776200387931*^9, 3.567776602946638*^9, {3.567776644458311*^9, 3.567776690712392*^9}, 3.5677767648593225`*^9, 3.5677768128606067`*^9, {3.5677768525314765`*^9, 3.5677768788019223`*^9}, 3.5677769279888086`*^9, {3.567776967488078*^9, 3.567776993633724*^9}, 3.5677848690954437`*^9, {3.5678310647505074`*^9, 3.5678310934077578`*^9}, 3.5678311344670296`*^9, 3.567831199675144*^9, 3.567831289047701*^9, 3.567831361837429*^9, 3.567831450710785*^9, { 3.5678315082904863`*^9, 3.5678315345453324`*^9}, 3.5678315954322395`*^9, { 3.5678316411715193`*^9, 3.5678316582535496`*^9}, 3.567831768764144*^9, 3.567832019893385*^9, 3.567832100077526*^9, 3.5678322729258294`*^9, 3.567832379177616*^9, 3.5678325254278727`*^9, {3.567832655329301*^9, 3.5678327336414385`*^9}, 3.567832787133932*^9, 3.567832841562428*^9, { 3.567832914601756*^9, 3.5678329434618073`*^9}, 3.5678330109631257`*^9, { 3.5678330751104383`*^9, 3.5678330972936773`*^9}, 3.5678331386025496`*^9, 3.5678332096606746`*^9, 3.5678333321832895`*^9, {3.56783336617575*^9, 3.5678333892013903`*^9}, {3.5678334726615367`*^9, 3.5678336001137605`*^9}, 3.5678336366178246`*^9, {3.5678336691282816`*^9, 3.5678336979571323`*^9}, 3.567833771511262*^9, {3.567833817468942*^9, 3.5678338420857854`*^9}, { 3.5678406904538107`*^9, 3.5678407102814455`*^9}, 3.5678408506972923`*^9}], Cell[BoxData["33.56310261756055`"], "Output", CellChangeTimes->{ 3.567773937947157*^9, 3.5677743370270576`*^9, 3.567774376370327*^9, 3.5677745394062138`*^9, 3.5677746263919663`*^9, 3.56777469670129*^9, 3.5677748963036404`*^9, 3.56777508398757*^9, 3.5677751634073095`*^9, 3.567775273839904*^9, 3.567775307177162*^9, 3.567775347643633*^9, 3.567775421681363*^9, 3.5677754545194206`*^9, 3.5677754901186833`*^9, { 3.5677755368095655`*^9, 3.567775561629209*^9}, 3.567775624310119*^9, 3.5677757594063563`*^9, 3.5677759372154684`*^9, 3.567776200387931*^9, 3.567776602946638*^9, {3.567776644458311*^9, 3.567776690712392*^9}, 3.5677767648593225`*^9, 3.5677768128606067`*^9, {3.5677768525314765`*^9, 3.5677768788019223`*^9}, 3.5677769279888086`*^9, {3.567776967488078*^9, 3.567776993633724*^9}, 3.5677848690954437`*^9, {3.5678310647505074`*^9, 3.5678310934077578`*^9}, 3.5678311344670296`*^9, 3.567831199675144*^9, 3.567831289047701*^9, 3.567831361837429*^9, 3.567831450710785*^9, { 3.5678315082904863`*^9, 3.5678315345453324`*^9}, 3.5678315954322395`*^9, { 3.5678316411715193`*^9, 3.5678316582535496`*^9}, 3.567831768764144*^9, 3.567832019893385*^9, 3.567832100077526*^9, 3.5678322729258294`*^9, 3.567832379177616*^9, 3.5678325254278727`*^9, {3.567832655329301*^9, 3.5678327336414385`*^9}, 3.567832787133932*^9, 3.567832841562428*^9, { 3.567832914601756*^9, 3.5678329434618073`*^9}, 3.5678330109631257`*^9, { 3.5678330751104383`*^9, 3.5678330972936773`*^9}, 3.5678331386025496`*^9, 3.5678332096606746`*^9, 3.5678333321832895`*^9, {3.56783336617575*^9, 3.5678333892013903`*^9}, {3.5678334726615367`*^9, 3.5678336001137605`*^9}, 3.5678336366178246`*^9, {3.5678336691282816`*^9, 3.5678336979571323`*^9}, 3.567833771511262*^9, {3.567833817468942*^9, 3.5678338420857854`*^9}, { 3.5678406904538107`*^9, 3.5678407102814455`*^9}, 3.5678408506972923`*^9}], Cell[BoxData["2.000821210624851`"], "Output", CellChangeTimes->{ 3.567773937947157*^9, 3.5677743370270576`*^9, 3.567774376370327*^9, 3.5677745394062138`*^9, 3.5677746263919663`*^9, 3.56777469670129*^9, 3.5677748963036404`*^9, 3.56777508398757*^9, 3.5677751634073095`*^9, 3.567775273839904*^9, 3.567775307177162*^9, 3.567775347643633*^9, 3.567775421681363*^9, 3.5677754545194206`*^9, 3.5677754901186833`*^9, { 3.5677755368095655`*^9, 3.567775561629209*^9}, 3.567775624310119*^9, 3.5677757594063563`*^9, 3.5677759372154684`*^9, 3.567776200387931*^9, 3.567776602946638*^9, {3.567776644458311*^9, 3.567776690712392*^9}, 3.5677767648593225`*^9, 3.5677768128606067`*^9, {3.5677768525314765`*^9, 3.5677768788019223`*^9}, 3.5677769279888086`*^9, {3.567776967488078*^9, 3.567776993633724*^9}, 3.5677848690954437`*^9, {3.5678310647505074`*^9, 3.5678310934077578`*^9}, 3.5678311344670296`*^9, 3.567831199675144*^9, 3.567831289047701*^9, 3.567831361837429*^9, 3.567831450710785*^9, { 3.5678315082904863`*^9, 3.5678315345453324`*^9}, 3.5678315954322395`*^9, { 3.5678316411715193`*^9, 3.5678316582535496`*^9}, 3.567831768764144*^9, 3.567832019893385*^9, 3.567832100077526*^9, 3.5678322729258294`*^9, 3.567832379177616*^9, 3.5678325254278727`*^9, {3.567832655329301*^9, 3.5678327336414385`*^9}, 3.567832787133932*^9, 3.567832841562428*^9, { 3.567832914601756*^9, 3.5678329434618073`*^9}, 3.5678330109631257`*^9, { 3.5678330751104383`*^9, 3.5678330972936773`*^9}, 3.5678331386025496`*^9, 3.5678332096606746`*^9, 3.5678333321832895`*^9, {3.56783336617575*^9, 3.5678333892013903`*^9}, {3.5678334726615367`*^9, 3.5678336001137605`*^9}, 3.5678336366178246`*^9, {3.5678336691282816`*^9, 3.5678336979571323`*^9}, 3.567833771511262*^9, {3.567833817468942*^9, 3.5678338420857854`*^9}, { 3.5678406904538107`*^9, 3.5678407102814455`*^9}, 3.5678408506972923`*^9}], Cell[BoxData["116.15374262699432`"], "Output", CellChangeTimes->{ 3.567773937947157*^9, 3.5677743370270576`*^9, 3.567774376370327*^9, 3.5677745394062138`*^9, 3.5677746263919663`*^9, 3.56777469670129*^9, 3.5677748963036404`*^9, 3.56777508398757*^9, 3.5677751634073095`*^9, 3.567775273839904*^9, 3.567775307177162*^9, 3.567775347643633*^9, 3.567775421681363*^9, 3.5677754545194206`*^9, 3.5677754901186833`*^9, { 3.5677755368095655`*^9, 3.567775561629209*^9}, 3.567775624310119*^9, 3.5677757594063563`*^9, 3.5677759372154684`*^9, 3.567776200387931*^9, 3.567776602946638*^9, {3.567776644458311*^9, 3.567776690712392*^9}, 3.5677767648593225`*^9, 3.5677768128606067`*^9, {3.5677768525314765`*^9, 3.5677768788019223`*^9}, 3.5677769279888086`*^9, {3.567776967488078*^9, 3.567776993633724*^9}, 3.5677848690954437`*^9, {3.5678310647505074`*^9, 3.5678310934077578`*^9}, 3.5678311344670296`*^9, 3.567831199675144*^9, 3.567831289047701*^9, 3.567831361837429*^9, 3.567831450710785*^9, { 3.5678315082904863`*^9, 3.5678315345453324`*^9}, 3.5678315954322395`*^9, { 3.5678316411715193`*^9, 3.5678316582535496`*^9}, 3.567831768764144*^9, 3.567832019893385*^9, 3.567832100077526*^9, 3.5678322729258294`*^9, 3.567832379177616*^9, 3.5678325254278727`*^9, {3.567832655329301*^9, 3.5678327336414385`*^9}, 3.567832787133932*^9, 3.567832841562428*^9, { 3.567832914601756*^9, 3.5678329434618073`*^9}, 3.5678330109631257`*^9, { 3.5678330751104383`*^9, 3.5678330972936773`*^9}, 3.5678331386025496`*^9, 3.5678332096606746`*^9, 3.5678333321832895`*^9, {3.56783336617575*^9, 3.5678333892013903`*^9}, {3.5678334726615367`*^9, 3.5678336001137605`*^9}, 3.5678336366178246`*^9, {3.5678336691282816`*^9, 3.5678336979571323`*^9}, 3.567833771511262*^9, {3.567833817468942*^9, 3.5678338420857854`*^9}, { 3.5678406904538107`*^9, 3.5678407102814455`*^9}, 3.5678408506972923`*^9}], Cell[BoxData["1.3205345528642047`"], "Output", CellChangeTimes->{ 3.567773937947157*^9, 3.5677743370270576`*^9, 3.567774376370327*^9, 3.5677745394062138`*^9, 3.5677746263919663`*^9, 3.56777469670129*^9, 3.5677748963036404`*^9, 3.56777508398757*^9, 3.5677751634073095`*^9, 3.567775273839904*^9, 3.567775307177162*^9, 3.567775347643633*^9, 3.567775421681363*^9, 3.5677754545194206`*^9, 3.5677754901186833`*^9, { 3.5677755368095655`*^9, 3.567775561629209*^9}, 3.567775624310119*^9, 3.5677757594063563`*^9, 3.5677759372154684`*^9, 3.567776200387931*^9, 3.567776602946638*^9, {3.567776644458311*^9, 3.567776690712392*^9}, 3.5677767648593225`*^9, 3.5677768128606067`*^9, {3.5677768525314765`*^9, 3.5677768788019223`*^9}, 3.5677769279888086`*^9, {3.567776967488078*^9, 3.567776993633724*^9}, 3.5677848690954437`*^9, {3.5678310647505074`*^9, 3.5678310934077578`*^9}, 3.5678311344670296`*^9, 3.567831199675144*^9, 3.567831289047701*^9, 3.567831361837429*^9, 3.567831450710785*^9, { 3.5678315082904863`*^9, 3.5678315345453324`*^9}, 3.5678315954322395`*^9, { 3.5678316411715193`*^9, 3.5678316582535496`*^9}, 3.567831768764144*^9, 3.567832019893385*^9, 3.567832100077526*^9, 3.5678322729258294`*^9, 3.567832379177616*^9, 3.5678325254278727`*^9, {3.567832655329301*^9, 3.5678327336414385`*^9}, 3.567832787133932*^9, 3.567832841562428*^9, { 3.567832914601756*^9, 3.5678329434618073`*^9}, 3.5678330109631257`*^9, { 3.5678330751104383`*^9, 3.5678330972936773`*^9}, 3.5678331386025496`*^9, 3.5678332096606746`*^9, 3.5678333321832895`*^9, {3.56783336617575*^9, 3.5678333892013903`*^9}, {3.5678334726615367`*^9, 3.5678336001137605`*^9}, 3.5678336366178246`*^9, {3.5678336691282816`*^9, 3.5678336979571323`*^9}, 3.567833771511262*^9, {3.567833817468942*^9, 3.5678338420857854`*^9}, { 3.5678406904538107`*^9, 3.5678407102814455`*^9}, 3.5678408506972923`*^9}], Cell[BoxData["1.0811476640938469`"], "Output", CellChangeTimes->{ 3.567773937947157*^9, 3.5677743370270576`*^9, 3.567774376370327*^9, 3.5677745394062138`*^9, 3.5677746263919663`*^9, 3.56777469670129*^9, 3.5677748963036404`*^9, 3.56777508398757*^9, 3.5677751634073095`*^9, 3.567775273839904*^9, 3.567775307177162*^9, 3.567775347643633*^9, 3.567775421681363*^9, 3.5677754545194206`*^9, 3.5677754901186833`*^9, { 3.5677755368095655`*^9, 3.567775561629209*^9}, 3.567775624310119*^9, 3.5677757594063563`*^9, 3.5677759372154684`*^9, 3.567776200387931*^9, 3.567776602946638*^9, {3.567776644458311*^9, 3.567776690712392*^9}, 3.5677767648593225`*^9, 3.5677768128606067`*^9, {3.5677768525314765`*^9, 3.5677768788019223`*^9}, 3.5677769279888086`*^9, {3.567776967488078*^9, 3.567776993633724*^9}, 3.5677848690954437`*^9, {3.5678310647505074`*^9, 3.5678310934077578`*^9}, 3.5678311344670296`*^9, 3.567831199675144*^9, 3.567831289047701*^9, 3.567831361837429*^9, 3.567831450710785*^9, { 3.5678315082904863`*^9, 3.5678315345453324`*^9}, 3.5678315954322395`*^9, { 3.5678316411715193`*^9, 3.5678316582535496`*^9}, 3.567831768764144*^9, 3.567832019893385*^9, 3.567832100077526*^9, 3.5678322729258294`*^9, 3.567832379177616*^9, 3.5678325254278727`*^9, {3.567832655329301*^9, 3.5678327336414385`*^9}, 3.567832787133932*^9, 3.567832841562428*^9, { 3.567832914601756*^9, 3.5678329434618073`*^9}, 3.5678330109631257`*^9, { 3.5678330751104383`*^9, 3.5678330972936773`*^9}, 3.5678331386025496`*^9, 3.5678332096606746`*^9, 3.5678333321832895`*^9, {3.56783336617575*^9, 3.5678333892013903`*^9}, {3.5678334726615367`*^9, 3.5678336001137605`*^9}, 3.5678336366178246`*^9, {3.5678336691282816`*^9, 3.5678336979571323`*^9}, 3.567833771511262*^9, {3.567833817468942*^9, 3.5678338420857854`*^9}, { 3.5678406904538107`*^9, 3.5678407102814455`*^9}, 3.5678408506972923`*^9}] }, Open ]], Cell[CellGroupData[{ Cell["\<\ plot1=Graphics[{Circle[{xc[[1]],yc[[1]]},Cr], Circle[{xc[[2]],yc[[2]]},Cr], Circle[{xc[[3]],yc[[3]]},CH], Circle[{xc[[4]],yc[[4]]},CH], Circle[{xc[[5]],yc[[5]]},CH], Circle[{xc[[6]],yc[[6]]},CH], Circle[{xc[[7]],yc[[7]]},CC], Circle[{xc[[8]],yc[[8]]},Cpi], Circle[{xc[[9]],yc[[9]]},Cpi], Circle[{xc[[10]],yc[[10]]},Cpi], Circle[{xc[[11]],yc[[11]]},Cpi], 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], {Thickness[0.008], Line[{{xn[[2]],yn[[2]]+0.06},{xn[[1]],yn[[1]]+0.06}}], Line[{{xn[[2]],yn[[2]]-0.06},{xn[[1]],yn[[1]]-0.06}}], Line[{{xn[[3]],yn[[3]]},{xn[[1]],yn[[1]]}}], Line[{{xn[[1]],yn[[1]]},{xn[[4]],yn[[4]]}}], Line[{{xn[[2]],yn[[2]]},{xn[[5]],yn[[5]]}}], Line[{{xn[[2]],yn[[2]]},{xn[[6]],yn[[6]]}}]}} ] /. t[[2]]; (* xy-plane *) 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}, {3.567248356181117*^9, 3.5672483876931725`*^9}, {3.567248623440787*^9, 3.567248633690005*^9}, 3.5672602418467417`*^9, {3.5672607929177094`*^9, 3.567260914441923*^9}, { 3.567260963628809*^9, 3.56726096965042*^9}, {3.5672610513009634`*^9, 3.5672610632349844`*^9}, {3.5672611071334615`*^9, 3.5672611271014967`*^9}, {3.5672643309709063`*^9, 3.5672643349645133`*^9}, {3.567264382794197*^9, 3.5672643891278086`*^9}, { 3.5672644516371183`*^9, 3.5672644630095387`*^9}, {3.5672644999192033`*^9, 3.567264505067212*^9}, {3.5672645512120934`*^9, 3.5672645514772935`*^9}, { 3.567264587497757*^9, 3.567264597388174*^9}, {3.56726465747948*^9, 3.567264679085518*^9}, {3.5672647158235826`*^9, 3.5672647395200243`*^9}, { 3.567264820562166*^9, 3.5672648250393744`*^9}, {3.5673459008509073`*^9, 3.567345907824119*^9}, 3.5673460056986914`*^9, {3.5677676200213814`*^9, 3.5677678327901554`*^9}, {3.567776562168166*^9, 3.5677766862975845`*^9}}], Cell[BoxData[ GraphicsBox[{CircleBox[{1.2477248187886139`, 0}, 0.26062448315986025`], CircleBox[{-1.2477248187886139`, 0}, 0.26062448315986025`], CircleBox[{2.0417485385104577`, 1.27450535285033}, 1.2409871716713938`], CircleBox[{2.0417485385104577`, -1.27450535285033}, 1.2409871716713938`], CircleBox[{-2.0417485385104577`, 1.27450535285033}, 1.2409871716713938`], CircleBox[{-2.0417485385104577`, -1.27450535285033}, 1.2409871716713938`], CircleBox[{0, 0}, 0.9871003356287535], CircleBox[{1.0977248187886137`, 0}, 1.2072995482981752`], CircleBox[{1.0977248187886137`, 0}, 1.2072995482981752`], CircleBox[{-1.0977248187886137`, 0}, 1.2072995482981752`], CircleBox[{-1.0977248187886137`, 0}, 1.2072995482981752`], DiskBox[{1.2477248187886139`, 0}, 0.08], DiskBox[{-1.2477248187886139`, 0}, 0.08], DiskBox[{2.328063435547241, 1.7340758415027397`}, 0.08], DiskBox[{2.328063435547241, -1.7340758415027397`}, 0.08], DiskBox[{-2.328063435547241, 1.7340758415027397`}, 0.08], DiskBox[{-2.328063435547241, -1.7340758415027397`}, 0.08], {Thickness[0.008], LineBox[{{-1.2477248187886139`, 0.06}, {1.2477248187886139`, 0.06}}], LineBox[{{-1.2477248187886139`, -0.06}, {1.2477248187886139`, -0.06}}], LineBox[{{2.328063435547241, 1.7340758415027397`}, { 1.2477248187886139`, 0}}], LineBox[{{1.2477248187886139`, 0}, { 2.328063435547241, -1.7340758415027397`}}], LineBox[{{-1.2477248187886139`, 0}, {-2.328063435547241, 1.7340758415027397`}}], LineBox[{{-1.2477248187886139`, 0}, {-2.328063435547241, -1.7340758415027397`}}]}}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, PlotRange->Automatic]], "Output", CellChangeTimes->{ 3.567773937962757*^9, 3.567774337042658*^9, 3.5677743764015274`*^9, 3.5677745394374137`*^9, 3.5677746264231663`*^9, 3.56777469671689*^9, 3.5677748963348403`*^9, 3.56777508401877*^9, 3.5677751634229097`*^9, 3.5677752738711033`*^9, 3.567775307192762*^9, 3.567775347659233*^9, 3.5677754217125635`*^9, 3.567775454535021*^9, 3.567775490134283*^9, { 3.5677755368251657`*^9, 3.567775561644809*^9}, 3.567775624341319*^9, 3.567775759421956*^9, 3.5677759372310686`*^9, 3.567776200419131*^9, 3.5677766029622383`*^9, {3.5677766444895105`*^9, 3.5677766907435923`*^9}, 3.5677767648905225`*^9, 3.567776812876207*^9, {3.567776852547076*^9, 3.567776878817522*^9}, 3.567776928020009*^9, {3.5677769675036783`*^9, 3.5677769936493244`*^9}, 3.567784869126644*^9, {3.567831064766107*^9, 3.5678310934233575`*^9}, 3.56783113449823*^9, 3.5678311996907444`*^9, 3.567831289063301*^9, 3.5678313618686295`*^9, 3.5678314507419853`*^9, { 3.5678315083060865`*^9, 3.567831534560932*^9}, 3.56783159546344*^9, { 3.5678316412027197`*^9, 3.56783165826915*^9}, 3.5678317687953444`*^9, 3.567832019908985*^9, 3.567832100093126*^9, 3.5678322729570293`*^9, 3.567832379208816*^9, 3.567832525443473*^9, {3.5678326553605013`*^9, 3.5678327336570387`*^9}, 3.5678327871495323`*^9, 3.567832841593628*^9, { 3.5678329146329565`*^9, 3.5678329434774075`*^9}, 3.5678330109943256`*^9, { 3.5678330751260386`*^9, 3.5678330973092775`*^9}, 3.56783313863375*^9, 3.5678332096918745`*^9, 3.5678333321988897`*^9, {3.5678333661913495`*^9, 3.56783338921699*^9}, {3.5678334726771364`*^9, 3.56783360012936*^9}, 3.5678336366490245`*^9, {3.5678336691594815`*^9, 3.5678336979727325`*^9}, 3.5678337715268617`*^9, {3.567833817484542*^9, 3.5678338421013856`*^9}, { 3.567840690485011*^9, 3.5678407102970457`*^9}, 3.5678408507128925`*^9}] }, Open ]], Cell[CellGroupData[{ Cell["\<\ plot2=Graphics[{Circle[{xc[[1]],yc[[1]]},Cr], Circle[{xc[[2]],zc[[2]]},Cr], Circle[{xc[[3]],zc[[3]]},CH], Circle[{xc[[4]],zc[[4]]},CH], Circle[{xc[[5]],zc[[5]]},CH], Circle[{xc[[6]],zc[[6]]},CH], Circle[{xc[[7]],zc[[7]]},CC], Circle[{xc[[8]],zc[[8]]},Cpi], Circle[{xc[[9]],zc[[9]]},Cpi], Circle[{xc[[10]],zc[[10]]},Cpi], Circle[{xc[[11]],zc[[11]]},Cpi], 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], {Thickness[0.008], Line[{{xn[[2]],zn[[2]]},{xn[[1]],zn[[1]]}}], Line[{{xn[[3]],zn[[3]]},{xn[[1]],zn[[1]]}}], Line[{{xn[[1]],zn[[1]]},{xn[[4]],zn[[4]]}}], Line[{{xn[[2]],zn[[2]]},{xn[[5]],zn[[5]]}}], Line[{{xn[[2]],zn[[2]]},{xn[[6]],zn[[6]]}}]}} ] /. t[[2]]; (* xz-plane *) Show[plot2,{AspectRatio \[Rule] Automatic,Axes -> True,GridLines -> \ Automatic, PlotRange \[Rule] {{-3.5,3.5},{-3.5,3.5}}, 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}, {3.5672484126064167`*^9, 3.567248437972061*^9}, { 3.567248597435541*^9, 3.5672486093851624`*^9}, {3.5672599486910267`*^9, 3.567260127888541*^9}, {3.5672601804762335`*^9, 3.5672601882294474`*^9}, 3.5672602282279177`*^9, {3.56734591933694*^9, 3.567345928369355*^9}, 3.5673459980702777`*^9, {3.5677679513759637`*^9, 3.567768010531267*^9}, 3.5677848584562254`*^9}], Cell[BoxData[ GraphicsBox[{CircleBox[{1.2477248187886139`, 0}, 0.26062448315986025`], CircleBox[{-1.2477248187886139`, 0}, 0.26062448315986025`], CircleBox[{2.0417485385104577`, 0}, 1.2409871716713938`], CircleBox[{2.0417485385104577`, 0}, 1.2409871716713938`], CircleBox[{-2.0417485385104577`, 0}, 1.2409871716713938`], CircleBox[{-2.0417485385104577`, 0}, 1.2409871716713938`], CircleBox[{0, 0}, 0.9871003356287535], CircleBox[{1.0977248187886137`, 1.3376117898781053`}, 1.2072995482981752`], CircleBox[{1.0977248187886137`, -1.3376117898781053`}, 1.2072995482981752`], CircleBox[{-1.0977248187886137`, 1.3376117898781053`}, 1.2072995482981752`], CircleBox[{-1.0977248187886137`, -1.3376117898781053`}, 1.2072995482981752`], DiskBox[{1.2477248187886139`, 0}, 0.08], DiskBox[{-1.2477248187886139`, 0}, 0.08], DiskBox[{2.328063435547241, 0}, 0.08], DiskBox[{2.328063435547241, 0}, 0.08], DiskBox[{-2.328063435547241, 0}, 0.08], DiskBox[{-2.328063435547241, 0}, 0.08], {Thickness[0.008], LineBox[{{-1.2477248187886139`, 0}, {1.2477248187886139`, 0}}], LineBox[{{2.328063435547241, 0}, {1.2477248187886139`, 0}}], LineBox[{{1.2477248187886139`, 0}, {2.328063435547241, 0}}], LineBox[{{-1.2477248187886139`, 0}, {-2.328063435547241, 0}}], LineBox[{{-1.2477248187886139`, 0}, {-2.328063435547241, 0}}]}}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, PlotRange->{{-3.5, 3.5}, {-3.5, 3.5}}]], "Output", CellChangeTimes->{ 3.5677739379939575`*^9, 3.567774337058258*^9, 3.567774376417127*^9, 3.5677745394530134`*^9, 3.5677746264387665`*^9, 3.56777469673249*^9, 3.5677748963504405`*^9, 3.56777508403437*^9, 3.5677751634541097`*^9, 3.5677752738867035`*^9, 3.567775307208362*^9, 3.5677753476748333`*^9, 3.5677754217281632`*^9, 3.567775454550621*^9, 3.5677754901498833`*^9, { 3.5677755368407655`*^9, 3.567775561660409*^9}, 3.5677756243569193`*^9, 3.5677757594375563`*^9, 3.567775937246669*^9, 3.567776200419131*^9, 3.567776602993438*^9, {3.5677766445051107`*^9, 3.5677766907591925`*^9}, 3.5677767648905225`*^9, 3.5677768128918066`*^9, {3.5677768525626764`*^9, 3.5677768788331223`*^9}, 3.567776928020009*^9, {3.5677769675192785`*^9, 3.567776993664924*^9}, 3.5677848691422443`*^9, {3.5678310647817073`*^9, 3.5678310934389577`*^9}, 3.5678311345138297`*^9, 3.5678311997063446`*^9, 3.5678312890945015`*^9, 3.567831361884229*^9, 3.5678314507575855`*^9, { 3.5678315083216867`*^9, 3.5678315345765324`*^9}, 3.5678315954790397`*^9, { 3.56783164121832*^9, 3.5678316583003497`*^9}, 3.567831768810944*^9, 3.5678320199245853`*^9, 3.5678321001243258`*^9, 3.5678322729726295`*^9, 3.567832379224416*^9, 3.567832525459073*^9, {3.5678326553761015`*^9, 3.567832733672639*^9}, 3.5678327871651325`*^9, 3.567832841609228*^9, { 3.5678329146485567`*^9, 3.567832943493007*^9}, 3.567833011009926*^9, { 3.5678330751572385`*^9, 3.5678330973248773`*^9}, 3.56783313864935*^9, 3.5678332097074747`*^9, 3.56783333221449*^9, {3.5678333662069497`*^9, 3.5678333892325897`*^9}, {3.5678334726927366`*^9, 3.5678336001449604`*^9}, 3.5678336366646247`*^9, {3.5678336691750813`*^9, 3.5678336979883327`*^9}, 3.5678337715424614`*^9, {3.567833817500142*^9, 3.5678338421169853`*^9}, { 3.5678406905006113`*^9, 3.567840710312646*^9}, 3.5678408507284927`*^9}] }, Open ]], Cell[CellGroupData[{ Cell["\<\ plot3=Graphics3D[{ Sphere[{xc[[1]],yc[[1]],zc[[1]]},Cr], Sphere[{xc[[2]],yc[[2]],zc[[2]]},Cr], {Opacity[0.6], Sphere[{xc[[3]],yc[[3]],zc[[3]]},CH], Sphere[{xc[[4]],yc[[4]],zc[[4]]},CH], Sphere[{xc[[5]],yc[[5]],zc[[5]]},CH], Sphere[{xc[[6]],yc[[6]],zc[[6]]},CH], Sphere[{xc[[7]],yc[[7]],zc[[7]]},CC], {Opacity[0.4],{Lighter[Green], Sphere[{xc[[8]],yc[[8]],zc[[8]]},Cpi], Sphere[{xc[[9]],yc[[9]],zc[[9]]},Cpi], Sphere[{xc[[10]],yc[[10]],zc[[10]]},Cpi], Sphere[{xc[[11]],yc[[11]],zc[[11]]},Cpi]}}, {Darker[Red,1], Sphere[{xn[[1]],yn[[1]],zn[[1]]},0.08], Sphere[{xn[[2]],yn[[2]],zn[[2]]},0.08], Sphere[{xn[[3]],yn[[3]],zn[[3]]},0.08], Sphere[{xn[[4]],yn[[4]],zn[[4]]},0.08], Sphere[{xn[[5]],yn[[5]],zn[[5]]},0.08], Sphere[{xn[[6]],yn[[6]],zn[[6]]},0.08], Cylinder[{{xn[[2]],yn[[2]]-0.1,zn[[2]]},{xn[[1]],yn[[1]]-0.1,zn[[1]]}},0.05], Cylinder[{{xn[[2]],yn[[2]]+0.1,zn[[2]]},{xn[[1]],yn[[1]]+0.1,zn[[1]]}},0.05], Cylinder[{{xn[[3]],yn[[3]],zn[[3]]},{xn[[1]],yn[[1]],zn[[1]]}},0.05], Cylinder[{{xn[[1]],yn[[1]],zn[[1]]},{xn[[4]],yn[[4]],zn[[4]]}},0.05], Cylinder[{{xn[[2]],yn[[2]],zn[[2]]},{xn[[5]],yn[[5]],zn[[5]]}},0.05], Cylinder[{{xn[[6]],yn[[6]],zn[[6]]},{xn[[2]],yn[[2]],zn[[2]]}},0.05]}}}] /. \ 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}, {3.5672484905285535`*^9, 3.5672485617270784`*^9}, { 3.567249135979687*^9, 3.567249154356519*^9}, {3.5672516676925335`*^9, 3.567251794832757*^9}, 3.5672519818002853`*^9, 3.56725201830435*^9, { 3.5672520987224913`*^9, 3.5672521037769003`*^9}, {3.567252239341138*^9, 3.5672522670467863`*^9}, {3.567252311678465*^9, 3.567252453279914*^9}, { 3.567252636627036*^9, 3.5672526483738565`*^9}, {3.56725271883918*^9, 3.567252725500392*^9}, 3.5672527761848807`*^9, {3.567252962511608*^9, 3.567252984897648*^9}, {3.567253050573763*^9, 3.5672531005094504`*^9}, { 3.5672531611155567`*^9, 3.567253313949025*^9}, {3.567253349579488*^9, 3.567253359189105*^9}, {3.567345074080655*^9, 3.567345086623077*^9}, { 3.5677674587014985`*^9, 3.5677675753897033`*^9}, {3.567768078883787*^9, 3.5677681632643356`*^9}, {3.5677731116825066`*^9, 3.567773141506959*^9}, 3.567773425798258*^9, {3.5677734894151692`*^9, 3.567773527884837*^9}, { 3.567773581907732*^9, 3.5677735889745445`*^9}, {3.5677749731493754`*^9, 3.567774985972598*^9}, {3.567775151722889*^9, 3.5677751566524973`*^9}, { 3.567775293979539*^9, 3.5677753012491517`*^9}}], Cell[BoxData[ Graphics3DBox[{SphereBox[{1.2477248187886139`, 0, 0}, 0.26062448315986025`], SphereBox[{-1.2477248187886139`, 0, 0}, 0.26062448315986025`], {Opacity[0.6], SphereBox[{2.0417485385104577`, 1.27450535285033, 0}, 1.2409871716713938`], SphereBox[{2.0417485385104577`, -1.27450535285033, 0}, 1.2409871716713938`], SphereBox[{-2.0417485385104577`, 1.27450535285033, 0}, 1.2409871716713938`], SphereBox[{-2.0417485385104577`, -1.27450535285033, 0}, 1.2409871716713938`], SphereBox[{0, 0, 0}, 0.9871003356287535], {Opacity[0.4], {RGBColor[ NCache[ Rational[1, 3], 0.3333333333333333], 1, NCache[ Rational[1, 3], 0.3333333333333333]], SphereBox[{1.0977248187886137`, 0, 1.3376117898781053`}, 1.2072995482981752`], SphereBox[{1.0977248187886137`, 0, -1.3376117898781053`}, 1.2072995482981752`], SphereBox[{-1.0977248187886137`, 0, 1.3376117898781053`}, 1.2072995482981752`], SphereBox[{-1.0977248187886137`, 0, -1.3376117898781053`}, 1.2072995482981752`]}}, {RGBColor[0, 0, 0], SphereBox[{1.2477248187886139`, 0, 0}, 0.08], SphereBox[{-1.2477248187886139`, 0, 0}, 0.08], SphereBox[{2.328063435547241, 1.7340758415027397`, 0}, 0.08], SphereBox[{2.328063435547241, -1.7340758415027397`, 0}, 0.08], SphereBox[{-2.328063435547241, 1.7340758415027397`, 0}, 0.08], SphereBox[{-2.328063435547241, -1.7340758415027397`, 0}, 0.08], CylinderBox[{{-1.2477248187886139`, -0.1, 0}, { 1.2477248187886139`, -0.1, 0}}, 0.05], CylinderBox[{{-1.2477248187886139`, 0.1, 0}, { 1.2477248187886139`, 0.1, 0}}, 0.05], CylinderBox[{{2.328063435547241, 1.7340758415027397`, 0}, { 1.2477248187886139`, 0, 0}}, 0.05], CylinderBox[{{1.2477248187886139`, 0, 0}, { 2.328063435547241, -1.7340758415027397`, 0}}, 0.05], CylinderBox[{{-1.2477248187886139`, 0, 0}, {-2.328063435547241, 1.7340758415027397`, 0}}, 0.05], CylinderBox[{{-2.328063435547241, -1.7340758415027397`, 0}, {-1.2477248187886139`, 0, 0}}, 0.05]}}}, AspectRatio->Automatic, Boxed->False, ImageSize->{347., 374.}, ViewPoint->{0.989438596541002, -2.8747493209240313`, 1.4855394998187677`}, ViewVertical->{0.08520446686896745, -0.32341989505724345`, 0.9724489502043494}]], "Output", CellChangeTimes->{ 3.5677739380095572`*^9, 3.5677743370894585`*^9, 3.567774376432727*^9, 3.5677745394686136`*^9, 3.5677746264543667`*^9, 3.56777469676369*^9, 3.56777489636604*^9, 3.5677750840499697`*^9, 3.5677751634697094`*^9, 3.5677752739023037`*^9, 3.567775307239562*^9, 3.5677753476904335`*^9, 3.5677754217281632`*^9, 3.567775454550621*^9, 3.5677754901498833`*^9, { 3.5677755368407655`*^9, 3.567775561676009*^9}, 3.5677756243725195`*^9, 3.5677757594531565`*^9, 3.567775937262269*^9, 3.567776200434731*^9, 3.567776603009038*^9, {3.567776644520711*^9, 3.567776690790392*^9}, 3.567776764906122*^9, 3.567776812907407*^9, {3.5677768525782766`*^9, 3.5677768788487225`*^9}, 3.567776928035609*^9, {3.5677769675348787`*^9, 3.567776993680524*^9}, 3.5677848694074445`*^9, {3.5678310647817073`*^9, 3.5678310934389577`*^9}, 3.56783113452943*^9, 3.5678311997063446`*^9, 3.5678312891101017`*^9, 3.567831361899829*^9, 3.5678314507887855`*^9, { 3.5678315083216867`*^9, 3.5678315345921326`*^9}, 3.5678315954946394`*^9, { 3.56783164123392*^9, 3.56783165831595*^9}, 3.567831768826544*^9, 3.567832019940185*^9, 3.567832100139926*^9, 3.5678322729882298`*^9, 3.567832379240016*^9, 3.567832525459073*^9, {3.567832655391701*^9, 3.5678327336882386`*^9}, 3.5678327871651325`*^9, 3.5678328416248283`*^9, { 3.5678329146641564`*^9, 3.567832943508607*^9}, 3.567833011041126*^9, { 3.567833075172838*^9, 3.5678330973404775`*^9}, 3.56783313866495*^9, 3.567833209723075*^9, 3.56783333223009*^9, {3.5678333662069497`*^9, 3.56783338924819*^9}, {3.567833472708337*^9, 3.5678336001449604`*^9}, 3.567833636680225*^9, {3.5678336691750813`*^9, 3.5678336980039325`*^9}, 3.5678337715580616`*^9, {3.5678338175157423`*^9, 3.5678338421325855`*^9}, { 3.5678406905006113`*^9, 3.567840710328246*^9}, 3.5678408507440925`*^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[{"[", RowBox[{"l", "+", "1"}], "]"}], "]"}], ",", RowBox[{"x", "[", RowBox[{"[", "l", "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", "l", "]"}], "]"}], ",", RowBox[{"z", "[", RowBox[{"[", "l", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"l", ",", "6"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", StyleBox[ RowBox[{"(*", " ", RowBox[{ "Cartesian", " ", "Coordinates", " ", "of", " ", "stationary", " ", "structure"}], " ", "*)"}], "Text"]}], "\[IndentingNewLine]", RowBox[{"u", "=", RowBox[{"MatrixForm", "[", "s", "]"}]}]}], "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}, 3.567249100786025*^9, {3.567336054232044*^9, 3.567336054762445*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"C", "0.6602672764321024`", "0.`", "0.`"}, {"C", RowBox[{"-", "0.6602672764321024`"}], "0.`", "0.`"}, {"H", "1.2319576246325825`", "0.9176330515788953`", "0.`"}, {"H", "1.2319576246325825`", RowBox[{"-", "0.9176330515788953`"}], "0.`"}, {"H", RowBox[{"-", "1.2319576246325825`"}], "0.9176330515788953`", "0.`"}, {"H", RowBox[{"-", "1.2319576246325825`"}], RowBox[{"-", "0.9176330515788953`"}], "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.5677739380095572`*^9, 3.567774337136258*^9, 3.5677743764951277`*^9, 3.5677745395310135`*^9, 3.5677746265167665`*^9, 3.5677746968104897`*^9, 3.567774896412841*^9, 3.5677750841123705`*^9, 3.5677751635165095`*^9, 3.5677752739647036`*^9, 3.567775307286362*^9, 3.567775347706033*^9, 3.567775421743763*^9, 3.5677754545662208`*^9, 3.5677754901654835`*^9, { 3.5677755368563657`*^9, 3.5677755616916094`*^9}, 3.5677756244349194`*^9, 3.5677757594531565`*^9, 3.567775937262269*^9, 3.5677762004503307`*^9, 3.5677766030714383`*^9, {3.5677766445831113`*^9, 3.5677766908371925`*^9}, 3.5677767649217224`*^9, 3.567776812907407*^9, {3.5677768525782766`*^9, 3.5677768788487225`*^9}, 3.567776928051209*^9, {3.5677769675348787`*^9, 3.567776993696124*^9}, 3.5677848694542446`*^9, {3.5678310647973075`*^9, 3.5678310934545574`*^9}, 3.5678311345918303`*^9, 3.5678311997219443`*^9, 3.5678312891569014`*^9, 3.5678313619466295`*^9, 3.5678314508355856`*^9, { 3.5678315083372865`*^9, 3.5678315345921326`*^9}, 3.5678315955570393`*^9, { 3.5678316412807198`*^9, 3.56783165836275*^9}, 3.5678317688889437`*^9, 3.567832019940185*^9, 3.567832100186726*^9, 3.5678322730350294`*^9, 3.567832379302416*^9, 3.567832525474673*^9, {3.567832655391701*^9, 3.567832733703839*^9}, 3.5678327871807327`*^9, 3.567832841687228*^9, { 3.567832914726557*^9, 3.567832943508607*^9}, 3.5678330110723257`*^9, { 3.5678330752196383`*^9, 3.5678330973404775`*^9}, 3.56783313866495*^9, 3.567833209754275*^9, 3.56783333223009*^9, {3.5678333662225494`*^9, 3.56783338924819*^9}, {3.567833472708337*^9, 3.5678336001605606`*^9}, 3.5678336367426248`*^9, {3.5678336691906815`*^9, 3.567833698019532*^9}, 3.5678337715580616`*^9, {3.5678338175157423`*^9, 3.5678338421481857`*^9}, { 3.5678406905162115`*^9, 3.567840710328246*^9}, 3.5678408507440925`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"(*", " ", RowBox[{"Distance", " ", "matrix", " ", "\[CapitalARing]"}], " ", "*)"}], "Text"], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Dij", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"x", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"y", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "6"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "6"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", "Dij", "]"}]}]}]], "Input", CellChangeTimes->{{3.5673121185644836`*^9, 3.5673123031908083`*^9}, { 3.567312519375988*^9, 3.5673125338996134`*^9}, 3.56731260009053*^9, { 3.5673126601194353`*^9, 3.5673126617262383`*^9}, {3.567312701537508*^9, 3.56731271390833*^9}, {3.5673127984604783`*^9, 3.567312798959679*^9}, { 3.56731301023085*^9, 3.5673130212756696`*^9}, {3.567313495251302*^9, 3.5673135176529417`*^9}, {3.5673140616882973`*^9, 3.567314081515932*^9}, { 3.5673141162571926`*^9, 3.5673141194083986`*^9}, {3.567314158034066*^9, 3.567314184585313*^9}, {3.567334145038291*^9, 3.567334156660311*^9}, { 3.5673342928329506`*^9, 3.567334322613403*^9}, {3.5673355964492397`*^9, 3.567335606152457*^9}, {3.5673360873353024`*^9, 3.56733610304453*^9}, { 3.5673361383629923`*^9, 3.567336138908993*^9}, {3.567336196395094*^9, 3.567336202213904*^9}, {3.567336253865595*^9, 3.56733627378683*^9}, { 3.567336312006897*^9, 3.5673363381057425`*^9}, {3.5673364031422567`*^9, 3.567336415684679*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0.`", "1.3205345528642047`", "1.0811476640938469`", "1.0811476640938469`", "2.1029896560751915`", "2.1029896560751915`"}, {"1.3205345528642047`", "0.`", "2.1029896560751915`", "2.1029896560751915`", "1.0811476640938469`", "1.0811476640938469`"}, {"1.0811476640938469`", "2.1029896560751915`", "0.`", "1.8352661031577906`", "2.463915249265165`", "3.072308582314186`"}, {"1.0811476640938469`", "2.1029896560751915`", "1.8352661031577906`", "0.`", "3.072308582314186`", "2.463915249265165`"}, {"2.1029896560751915`", "1.0811476640938469`", "2.463915249265165`", "3.072308582314186`", "0.`", "1.8352661031577906`"}, {"2.1029896560751915`", "1.0811476640938469`", "3.072308582314186`", "2.463915249265165`", "1.8352661031577906`", "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.567773938025157*^9, 3.567774337136258*^9, 3.5677743764951277`*^9, 3.5677745395466137`*^9, 3.5677746265323668`*^9, 3.56777469682609*^9, 3.5677748964284406`*^9, 3.5677750841279707`*^9, 3.5677751635321093`*^9, 3.5677752739647036`*^9, 3.5677753073019624`*^9, 3.567775347721633*^9, 3.567775421759363*^9, 3.5677754545818205`*^9, 3.5677754901810837`*^9, { 3.5677755368719654`*^9, 3.5677755616916094`*^9}, 3.5677756244349194`*^9, 3.5677757594687567`*^9, 3.5677759372778687`*^9, 3.5677762004503307`*^9, 3.5677766031026382`*^9, {3.567776644598711*^9, 3.5677766908527927`*^9}, 3.5677767649373226`*^9, 3.567776812923007*^9, {3.5677768525938764`*^9, 3.5677768788643227`*^9}, 3.5677769280668087`*^9, {3.5677769675504785`*^9, 3.567776993696124*^9}, 3.5677848694698443`*^9, {3.567831064812907*^9, 3.5678310934701576`*^9}, 3.56783113460743*^9, 3.567831199737544*^9, 3.5678312891725016`*^9, 3.5678313619622293`*^9, 3.567831450851186*^9, { 3.567831508352886*^9, 3.567831534607733*^9}, 3.5678315955726395`*^9, { 3.56783164129632*^9, 3.5678316583783503`*^9}, 3.567831768904544*^9, 3.567832019955785*^9, 3.567832100217926*^9, 3.56783227306623*^9, 3.567832379318016*^9, 3.567832525490273*^9, {3.567832655407301*^9, 3.567832733703839*^9}, 3.567832787196333*^9, 3.5678328417028284`*^9, { 3.567832914742157*^9, 3.567832943524207*^9}, 3.567833011103526*^9, { 3.567833075235238*^9, 3.567833097356077*^9}, 3.56783313868055*^9, 3.5678332097698746`*^9, 3.5678333322456903`*^9, {3.5678333662381496`*^9, 3.56783338926379*^9}, {3.567833472723937*^9, 3.567833600176161*^9}, 3.5678336367426248`*^9, {3.5678336692062817`*^9, 3.567833698019532*^9}, 3.567833771573662*^9, {3.5678338175313425`*^9, 3.5678338421481857`*^9}, { 3.567840690531811*^9, 3.5678407103438463`*^9}, 3.567840850759692*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"plot4", "=", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", "\n", RowBox[{ RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "Cr"}], "]"}], ",", " ", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "Cr"}], "]"}], ",", "\n", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.6", "]"}], ",", "\n", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.4", "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Lighter", "[", "Green", "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "8", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "8", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "8", "]"}], "]"}]}], "}"}], ",", "Cpi"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "9", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "9", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "9", "]"}], "]"}]}], "}"}], ",", "Cpi"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "10", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "10", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "10", "]"}], "]"}]}], "}"}], ",", "Cpi"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "11", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "11", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "11", "]"}], "]"}]}], "}"}], ",", "Cpi"}], "]"}]}], "}"}]}], "}"}], ",", "\n", RowBox[{"{", RowBox[{ RowBox[{"Darker", "[", RowBox[{"Red", ",", "1"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\n", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", "\n", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "4", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", "\n", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "5", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", "\n", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], "-", "0.1"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], "-", "0.1"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], "+", "0.1"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], "+", "0.1"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.06"}], "]"}]}], "}"}]}], "}"}]}], "}"}], "]"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}]], "Input", CellChangeTimes->{{3.567774789755453*^9, 3.567774875586804*^9}, 3.5677749531189404`*^9, {3.5677750023370266`*^9, 3.5677750700723457`*^9}, { 3.5677752442934513`*^9, 3.5677752481466584`*^9}, 3.5677763957626743`*^9, 3.5677764786924195`*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{"overlapped", " ", "\[Pi]", " ", "system"}], " ", "*)"}], "\n", RowBox[{"Show", "[", RowBox[{"plot4", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", RowBox[{"Boxed", "->", "False"}]}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.567774789755453*^9, 3.567774875586804*^9}, 3.5677749531189404`*^9, {3.5677750023370266`*^9, 3.5677750700723457`*^9}, { 3.5677752442934513`*^9, 3.5677752481466584`*^9}, 3.5677763957626743`*^9, { 3.5677764786924195`*^9, 3.567776510235675*^9}}], Cell[BoxData[ Graphics3DBox[{SphereBox[{1.2477248187886139`, 0, 0}, 0.26062448315986025`], SphereBox[{-1.2477248187886139`, 0, 0}, 0.26062448315986025`], {Opacity[0.6], {Opacity[0.4], {RGBColor[ NCache[ Rational[1, 3], 0.3333333333333333], 1, NCache[ Rational[1, 3], 0.3333333333333333]], SphereBox[{1.0977248187886137`, 0, 1.3376117898781053`}, 1.2072995482981752`], SphereBox[{1.0977248187886137`, 0, -1.3376117898781053`}, 1.2072995482981752`], SphereBox[{-1.0977248187886137`, 0, 1.3376117898781053`}, 1.2072995482981752`], SphereBox[{-1.0977248187886137`, 0, -1.3376117898781053`}, 1.2072995482981752`]}}, {RGBColor[0, 0, 0], SphereBox[{1.2477248187886139`, 0, 0}, 0.08], SphereBox[{-1.2477248187886139`, 0, 0}, 0.08], SphereBox[{2.328063435547241, 1.7340758415027397`, 0}, 0.08], SphereBox[{2.328063435547241, -1.7340758415027397`, 0}, 0.08], SphereBox[{-2.328063435547241, 1.7340758415027397`, 0}, 0.08], SphereBox[{-2.328063435547241, -1.7340758415027397`, 0}, 0.08], CylinderBox[{{2.328063435547241, 1.7340758415027397`, 0}, { 1.2477248187886139`, 0, 0}}, 0.05], CylinderBox[{{1.2477248187886139`, 0, 0}, { 2.328063435547241, -1.7340758415027397`, 0}}, 0.05], CylinderBox[{{-1.2477248187886139`, 0, 0}, {-2.328063435547241, 1.7340758415027397`, 0}}, 0.05], CylinderBox[{{-1.2477248187886139`, -0.1, 0}, { 1.2477248187886139`, -0.1, 0}}, 0.05], CylinderBox[{{-1.2477248187886139`, 0.1, 0}, { 1.2477248187886139`, 0.1, 0}}, 0.05], CylinderBox[{{-2.328063435547241, -1.7340758415027397`, 0}, {-1.2477248187886139`, 0, 0}}, 0.06]}}}, AspectRatio->Automatic, Boxed->False, ImageSize->{336.1915991035843, 364.}, ViewPoint->{0.6058715144117951, -2.6530321935182886`, 2.0110544219836273`}, ViewVertical->{-0.0028421488046460244`, 0.03527332099563409, 0.9998368924798626}]], "Output", CellChangeTimes->{ 3.5677748964440403`*^9, 3.56777508415917*^9, 3.5677751635477095`*^9, 3.567775273995904*^9, 3.567775307317562*^9, 3.567775347737233*^9, 3.5677754217749634`*^9, 3.5677754545974207`*^9, 3.567775490196684*^9, { 3.5677755368875656`*^9, 3.5677755617228093`*^9}, 3.5677756244661193`*^9, 3.567775759499956*^9, 3.5677759373090687`*^9, 3.567776200481531*^9, 3.5677766031182384`*^9, {3.567776644614311*^9, 3.567776690883992*^9}, 3.5677767649529223`*^9, 3.5677768129542065`*^9, {3.5677768526250763`*^9, 3.5677768788955226`*^9}, 3.567776928082409*^9, {3.5677769675816784`*^9, 3.5677769937273245`*^9}, 3.5677848695010443`*^9, {3.5678310648441076`*^9, 3.567831093485758*^9}, 3.56783113463863*^9, 3.5678311997531443`*^9, 3.5678312892037015`*^9, 3.5678313619778295`*^9, 3.5678314508667855`*^9, { 3.5678315083684864`*^9, 3.5678315346389327`*^9}, 3.5678315955882397`*^9, { 3.56783164132752*^9, 3.56783165839395*^9}, 3.567831768920144*^9, 3.567832019986985*^9, 3.567832100249126*^9, 3.5678322730818295`*^9, 3.5678323793336163`*^9, 3.5678325255058727`*^9, {3.567832655422901*^9, 3.5678327337350388`*^9}, 3.567832787211933*^9, 3.5678328417184286`*^9, { 3.5678329147577567`*^9, 3.5678329435554075`*^9}, 3.567833011150326*^9, { 3.5678330752508383`*^9, 3.5678330973872776`*^9}, 3.56783313871175*^9, 3.567833209801075*^9, 3.56783333226129*^9, {3.56783336625375*^9, 3.5678333892793903`*^9}, {3.5678334727551365`*^9, 3.5678336001917605`*^9}, 3.5678336367738247`*^9, {3.567833669221882*^9, 3.5678336980507326`*^9}, 3.5678337716048613`*^9, {3.5678338175625424`*^9, 3.5678338421793857`*^9}, { 3.567840690547411*^9, 3.567840710375046*^9}, 3.5678408507908926`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"plot6", "=", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "Cr"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "Cr"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.6", "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "7", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "7", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "7", "]"}], "]"}]}], "}"}], ",", "CC"}], "]"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"Darker", "[", RowBox[{"Red", ",", "1"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], "-", "0.1"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], "-", "0.1"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], "+", "0.1"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], "+", "0.1"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "4", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "5", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}]}], "}"}]}], "}"}], "]"}], "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", "\n", StyleBox[ RowBox[{"(*", " ", RowBox[{ "the", " ", "sigma", " ", "skeleton", " ", "of", " ", "the", " ", "molecule"}], " ", "*)"}], "Text"], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plot6", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", RowBox[{"Boxed", "\[Rule]", "False"}]}], "}"}]}], "]"}]}]}]], "Input", CellChangeTimes->{{3.567344022482808*^9, 3.567344060827675*^9}, { 3.5673441070661564`*^9, 3.5673441366438084`*^9}, {3.5673441903391027`*^9, 3.5673442015243225`*^9}, 3.56734427418925*^9, {3.567344526535293*^9, 3.567344545083726*^9}, 3.567345677474115*^9, 3.5673459753410378`*^9, { 3.5677739536719847`*^9, 3.567774072247793*^9}, {3.567774133930301*^9, 3.567774203350423*^9}, {3.5677742495577044`*^9, 3.567774257950519*^9}, { 3.5677742905857763`*^9, 3.567774324812237*^9}, 3.5677743668231106`*^9, { 3.567774432717626*^9, 3.5677745302177973`*^9}, {3.5677745831642904`*^9, 3.567774614442345*^9}, {3.567774646687602*^9, 3.5677746566404195`*^9}, { 3.5677746894628773`*^9, 3.5677746897904778`*^9}}], Cell[BoxData[ Graphics3DBox[{SphereBox[{1.2477248187886139`, 0, 0}, 0.26062448315986025`], SphereBox[{-1.2477248187886139`, 0, 0}, 0.26062448315986025`], {Opacity[0.6], SphereBox[{2.0417485385104577`, 1.27450535285033, 0}, 1.2409871716713938`], SphereBox[{2.0417485385104577`, -1.27450535285033, 0}, 1.2409871716713938`], SphereBox[{-2.0417485385104577`, 1.27450535285033, 0}, 1.2409871716713938`], SphereBox[{-2.0417485385104577`, -1.27450535285033, 0}, 1.2409871716713938`], SphereBox[{0, 0, 0}, 0.9871003356287535]}, {RGBColor[0, 0, 0], SphereBox[{1.2477248187886139`, 0, 0}, 0.08], SphereBox[{-1.2477248187886139`, 0, 0}, 0.08], SphereBox[{2.328063435547241, 1.7340758415027397`, 0}, 0.08], SphereBox[{2.328063435547241, -1.7340758415027397`, 0}, 0.08], SphereBox[{-2.328063435547241, 1.7340758415027397`, 0}, 0.08], SphereBox[{-2.328063435547241, -1.7340758415027397`, 0}, 0.08], CylinderBox[{{-1.2477248187886139`, -0.1, 0}, { 1.2477248187886139`, -0.1, 0}}, 0.05], CylinderBox[{{-1.2477248187886139`, 0.1, 0}, { 1.2477248187886139`, 0.1, 0}}, 0.05], CylinderBox[{{1.2477248187886139`, 0, 0}, { 2.328063435547241, -1.7340758415027397`, 0}}, 0.05], CylinderBox[{{1.2477248187886139`, 0, 0}, { 2.328063435547241, 1.7340758415027397`, 0}}, 0.05], CylinderBox[{{-1.2477248187886139`, 0, 0}, {-2.328063435547241, 1.7340758415027397`, 0}}, 0.05], CylinderBox[{{-2.328063435547241, -1.7340758415027397`, 0}, {-1.2477248187886139`, 0, 0}}, 0.05]}}, AspectRatio->Automatic, AutomaticImageSize->True, Boxed->False, ImageSize->{371.1275346379089, 317.9804574897849}, ViewPoint->{1.020892622353606, -1.594001508917872, 2.80480613290679}, ViewVertical->{-0.006637185673381674, 0.46928170722035956`, 2.3907718241138327`}]], "Output", CellChangeTimes->{ 3.5677743766823273`*^9, 3.567774539562214*^9, 3.5677746265635667`*^9, 3.56777469684169*^9, 3.567774896522041*^9, 3.56777508423717*^9, 3.56777516362571*^9, 3.567775274058304*^9, 3.567775307395562*^9, 3.5677753477528334`*^9, 3.5677754217905636`*^9, 3.567775454628621*^9, 3.5677754902122836`*^9, {3.567775536903166*^9, 3.567775561738409*^9}, 3.5677756245285196`*^9, 3.5677757595155563`*^9, 3.567775937324669*^9, 3.5677762004971313`*^9, 3.5677766031338387`*^9, {3.567776644645511*^9, 3.5677766908995924`*^9}, 3.567776764968522*^9, 3.5677768129698067`*^9, { 3.5677768526406765`*^9, 3.5677768789111223`*^9}, 3.567776928098009*^9, { 3.5677769675972786`*^9, 3.5677769937429247`*^9}, 3.5677848695322447`*^9, { 3.567831064859708*^9, 3.567831093501358*^9}, 3.56783113465423*^9, 3.5678311997687445`*^9, 3.5678312892193017`*^9, 3.5678313620090294`*^9, 3.5678314508823853`*^9, {3.5678315083840866`*^9, 3.567831534654533*^9}, 3.56783159560384*^9, {3.56783164134312*^9, 3.56783165842515*^9}, 3.5678317689357443`*^9, 3.5678320200025854`*^9, 3.5678321002647257`*^9, 3.5678322730974293`*^9, 3.5678323793492165`*^9, 3.567832525521473*^9, { 3.5678326554385014`*^9, 3.567832733750639*^9}, 3.567832787227533*^9, 3.5678328417340283`*^9, {3.567832914773357*^9, 3.5678329435710073`*^9}, 3.567833011165926*^9, {3.5678330752820387`*^9, 3.567833097402878*^9}, 3.56783313872735*^9, 3.5678332098166747`*^9, 3.56783333229249*^9, { 3.5678333662849493`*^9, 3.5678333892949905`*^9}, {3.5678334727707367`*^9, 3.5678336002229605`*^9}, 3.567833636789425*^9, {3.5678336692530823`*^9, 3.5678336980663323`*^9}, 3.5678337716204615`*^9, {3.567833817578142*^9, 3.5678338421949854`*^9}, {3.567840690563011*^9, 3.567840710390646*^9}, 3.567840850806493*^9}] }, Open ]] }, WindowToolbars->"EditBar", WindowSize->{900, 964}, WindowMargins->{{Automatic, 46}, {Automatic, 20}}, PrintingCopies->1, PrintingPageRange->{32000, 32000}, PrintingOptions->{"Magnification"->1., "PaperOrientation"->"Portrait", "PaperSize"->{595.2, 841.8}}, PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}}, FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (November 20, 2012)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[557, 20, 2779, 40, 31, "Input"], Cell[3339, 62, 6616, 97, 184, "Input"], Cell[9958, 161, 1319, 26, 72, "Input"], Cell[11280, 189, 784, 15, 82, "Input"], Cell[12067, 206, 4874, 69, 82, "Input"], Cell[16944, 277, 5476, 79, 133, "Input"], Cell[22423, 358, 5706, 84, 150, "Input"], Cell[28132, 444, 4692, 67, 82, "Input"], Cell[32827, 513, 5269, 73, 82, "Input"], Cell[38099, 588, 4639, 66, 65, "Input"], Cell[42741, 656, 611, 15, 116, "Input"], Cell[43355, 673, 384, 11, 82, "Input"], Cell[43742, 686, 571, 14, 116, "Input"], Cell[CellGroupData[{ Cell[44338, 704, 2150, 32, 82, "Input"], Cell[46491, 738, 2372, 38, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[48900, 781, 1002, 22, 150, "Input"], Cell[49905, 805, 1929, 26, 31, "Output"], Cell[51837, 833, 1947, 27, 31, "Output"], Cell[53787, 862, 1930, 26, 31, "Output"], Cell[55720, 890, 1929, 26, 31, "Output"], Cell[57652, 918, 1929, 26, 31, "Output"], Cell[59584, 946, 1930, 26, 31, "Output"], Cell[61517, 974, 1930, 26, 31, "Output"], Cell[63450, 1002, 1930, 26, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[65417, 1033, 2742, 53, 490, "Input"], Cell[68162, 1088, 3632, 59, 297, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[71831, 1152, 1996, 43, 473, "Input"], Cell[73830, 1197, 3431, 56, 375, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[77298, 1258, 4215, 71, 507, "Input"], Cell[81516, 1331, 4275, 75, 389, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[85828, 1411, 3008, 67, 152, "Input"], Cell[88839, 1480, 2912, 52, 102, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[91788, 1537, 2050, 47, 72, "Input"], Cell[93841, 1586, 3228, 53, 102, "Output"] }, Open ]], Cell[97084, 1642, 11278, 288, 452, "Input"], Cell[CellGroupData[{ Cell[108387, 1934, 618, 13, 52, "Input"], Cell[109008, 1949, 3762, 65, 379, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[112807, 2019, 12547, 311, 512, "Input"], Cell[125357, 2332, 3717, 62, 333, "Output"] }, Open ]] } ] *) (* End of internal cache information *)