(* 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[    136902,       2773]
NotebookOptionsPosition[    134210,       2689]
NotebookOutlinePosition[    134808,       2711]
CellTagsIndexPosition[    134765,       2708]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell["\<\
(* HCOOH Formic acid, with \"orthogonal\" \[Sigma],\[Pi] LPs, 18.01.2013, \
corr 25.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.5675116107397537`*^9, 3.5675116255753803`*^9}, {
   3.5675120779761744`*^9, 3.567512078787376*^9}, {3.568101252202099*^9, 
   3.568101261546515*^9}}],

Cell["\<\
Clear[ZC,ZN,ZO,sig1,sig2,sig3,sig4,sig5,sig6,sig7,sig8,sig9,k1,k2,k3,k4k5,k6,\
k7,k8,k9,
k10,k11,nc,w1,w2,w3,pnh,pch,vee,vne,vnn,xn,yn,zn,xc,yc,zc,c1,c2,c3,cs1,ss1,\
cs2,ss2,cs3,
ss3,pi,i,j,d,t];
Clear[Cr,Nr,Oa,NH,CH,CO,CN,OL,NM,OC,OO,OM,w1,w2,w3,pch,pnh];
ZC=6.0; ZN=8.0; ZO=8.0; sig1=0.3; sig2=0.3; sig3=0.3; sig4=0.3; sig5=0.3; \
sig6=0.3;
sig7=0.3; sig8=0.3; sig9=0.3;
nc=17; nn=6; pi=0.0;
k1=1.0; k2=0.995; k3=0.995; k4=1.255; k5=1.46; k6=1.293; k7=1.24; k8=1.85; \
k9=1.1; k10=1.1; 
k11=0.78; k12=0.99;  
bohr=0.529177; rad=57.29578;

(* 17 Variables: Cr,Nr,Or,NH,CH,CN,CO,OL,NM,OC,OO,OM,pnh,pch,w1,w2,w3 *)\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{
  3.5668390630515456`*^9, {3.5668394752510695`*^9, 3.5668394888386936`*^9}, {
   3.566840005402401*^9, 3.5668400293172426`*^9}, {3.566840906678384*^9, 
   3.566840923027212*^9}, {3.5668417775723133`*^9, 3.5668417895999346`*^9}, {
   3.566841995894697*^9, 3.5668420196535387`*^9}, {3.566842057031204*^9, 
   3.5668420654552193`*^9}, {3.5668422519691467`*^9, 3.5668422563371544`*^9}, 
   3.5668422887072115`*^9, {3.566842337504097*^9, 3.566842360326937*^9}, {
   3.5668424117758274`*^9, 3.5668424218690453`*^9}, {3.566842627885008*^9, 
   3.566842684060706*^9}, {3.566842771030859*^9, 3.5668427769432697`*^9}, {
   3.5668428297025623`*^9, 3.5668428906206694`*^9}, {3.5668429241763277`*^9, 
   3.566842930135538*^9}, {3.5668871695283947`*^9, 3.5668871699963956`*^9}, {
   3.5668872245808916`*^9, 3.566887225407693*^9}, {3.5668875348030367`*^9, 
   3.566887535037037*^9}, {3.566887597327946*^9, 3.566887637981618*^9}, 
   3.5668876998201265`*^9, {3.5668877336253853`*^9, 3.5668877352165885`*^9}, {
   3.566887814495928*^9, 3.5668878152135286`*^9}, {3.566887857926404*^9, 
   3.566887961448186*^9}, 3.5668880086070685`*^9, {3.5668880951092205`*^9, 
   3.5668881570569296`*^9}, {3.5668882853983545`*^9, 
   3.5668884070317683`*^9}, {3.566888509258748*^9, 3.5668885097423487`*^9}, {
   3.566888582719277*^9, 3.566888583795679*^9}, {3.566888668176227*^9, 
   3.566888671935834*^9}, {3.5668887042746906`*^9, 3.566888705070292*^9}, {
   3.566888771822809*^9, 3.566888808857274*^9}, {3.5668889463559155`*^9, 
   3.566888956839134*^9}, {3.5669685950612297`*^9, 3.5669685991952367`*^9}, {
   3.566968722700654*^9, 3.5669687698127365`*^9}, {3.5669688330864477`*^9, 
   3.566968848015674*^9}, {3.5669689606010714`*^9, 3.566968973471094*^9}, {
   3.5669881655485907`*^9, 3.5669881659853916`*^9}, {3.566996986024498*^9, 
   3.5669969910321074`*^9}, {3.5669973015818524`*^9, 
   3.5669973626247597`*^9}, {3.566997996329073*^9, 3.566998005049488*^9}, {
   3.56705755136038*^9, 3.5670575543243856`*^9}, {3.5670590414593973`*^9, 
   3.5670590456714053`*^9}, {3.567061276787324*^9, 3.567061372633892*^9}, {
   3.567061756144966*^9, 3.5670617569249673`*^9}, {3.567061904329626*^9, 
   3.567061915546046*^9}, {3.567062309945938*^9, 3.5670623507556105`*^9}, {
   3.567062381019663*^9, 3.5670625124810944`*^9}, {3.5670638527106485`*^9, 
   3.5670638629598665`*^9}, 3.5670687979573345`*^9, {3.567069322009055*^9, 
   3.5670693806495576`*^9}, {3.567069419290826*^9, 3.567069451957283*^9}, {
   3.5670694881649466`*^9, 3.5670695381162343`*^9}, {3.56708248268017*^9, 
   3.567082487516179*^9}, {3.567085712191843*^9, 3.56708586995492*^9}, {
   3.5670863390009437`*^9, 3.5670863409197474`*^9}, {3.5670863724006023`*^9, 
   3.5670863736018047`*^9}, {3.567086415363078*^9, 3.56708641629908*^9}, {
   3.5670864537547455`*^9, 3.567086586245778*^9}, {3.56708706580622*^9, 
   3.567087066742222*^9}, {3.5670871956452484`*^9, 3.567087197376851*^9}, {
   3.567087313846656*^9, 3.5670873147670574`*^9}, {3.567087346825114*^9, 
   3.5670873478703156`*^9}, {3.5670873954503994`*^9, 3.567087396292801*^9}, {
   3.5670874292868586`*^9, 3.5670874532017007`*^9}, {3.567087493090971*^9, 
   3.5670874935745716`*^9}, {3.5670979031336823`*^9, 
   3.5670979250517206`*^9}, {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.567510445168107*^9, 3.567510445464507*^9}, {3.5675104830761733`*^9, 
   3.567510486024578*^9}, {3.567516190147399*^9, 3.5675161919570017`*^9}, {
   3.567516252282308*^9, 3.5675162528595095`*^9}, {3.56751628775677*^9, 
   3.5675162893947735`*^9}, {3.5675163515452824`*^9, 
   3.5675163533704853`*^9}, {3.567516397939764*^9, 3.5675164719774942`*^9}, {
   3.567516528870794*^9, 3.5675165297443953`*^9}, 3.5675165885720987`*^9, {
   3.567516622346158*^9, 3.56751662343816*^9}, {3.5675166547474146`*^9, 
   3.5675166556522164`*^9}, {3.567516696259088*^9, 3.567516698209091*^9}, {
   3.567516730906749*^9, 3.5675167601879997`*^9}, {3.567516795974463*^9, 
   3.567516822166909*^9}, {3.567516867172988*^9, 3.5675169111650653`*^9}, {
   3.5675324036190743`*^9, 3.5675324740843983`*^9}}],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
   RowBox[{"Cr", "=", "0.26"}], ";", " ", 
   RowBox[{"Nr", "=", "0.22"}], ";", " ", 
   RowBox[{"Oa", "=", "0.186"}], ";", " ", 
   RowBox[{"NH", "=", "1.2"}], ";", " ", 
   RowBox[{"CH", "=", "1.25"}], ";", " ", 
   RowBox[{"CO", "=", "0.928"}], ";", " ", 
   RowBox[{"CN", "=", "1.04"}], ";", " ", 
   RowBox[{"OL", "=", "1.0"}], ";", "\[IndentingNewLine]", 
   RowBox[{"NM", "=", "1.1"}], ";", " ", 
   RowBox[{"OC", "=", "1.4"}], ";", " ", 
   RowBox[{"OM", "=", "1.05"}], ";", " ", 
   RowBox[{"w1", "=", "1.9"}], ";", " ", 
   RowBox[{"w2", "=", "1.0"}], ";", " ", 
   RowBox[{"w3", "=", "1.0"}], ";", " ", 
   RowBox[{"pch", "=", "0.58"}], ";", 
   RowBox[{"pnh", "=", "0.5"}], ";"}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.5670689879656677`*^9, 3.5670691959452333`*^9}, 
   3.5670698299303465`*^9, {3.567077287748646*^9, 3.5670773058290777`*^9}, {
   3.567077396527637*^9, 3.567077402346447*^9}, {3.5670803249623804`*^9, 
   3.5670804001545124`*^9}, 3.5670804667198296`*^9, 3.5670805071239004`*^9, {
   3.5670810452624454`*^9, 3.567081083794513*^9}, {3.5670811174593725`*^9, 
   3.5670811463974233`*^9}, {3.5670811875190954`*^9, 
   3.5670812332895756`*^9}, {3.5670812684832373`*^9, 3.567081269356839*^9}, {
   3.5670852872782965`*^9, 3.567085293970708*^9}}],

Cell["\<\
(* C He shell *)
Ekin = 2.25*k1/Cr^2;
vee=3.0*sig1/Cr;
vne=-3.0*ZC/Cr;

(* O He shell *)
Ekin = Ekin+2.25*k2/Nr^2;
vee=vee+3.0*sig2/Nr;
vne=vne-3.0*ZN/Nr;

(* O He shell *)
Ekin = Ekin+2.25*k3/Oa^2;
vee=vee+3.0*sig3/Oa;
vne=vne-3.0*ZO/Oa;\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668390630515456`*^9, 3.5668390682151546`*^9}, {
  3.566971938334302*^9, 3.566971944714713*^9}, {3.566972152507078*^9, 
  3.5669721651899004`*^9}, {3.5669970287685733`*^9, 3.5669972828462195`*^9}, {
  3.566997377709986*^9, 3.5669974183012576`*^9}, {3.5669979768290386`*^9, 
  3.566997987171857*^9}, {3.566998792414071*^9, 3.5669988056272945`*^9}, {
  3.567063874659887*^9, 3.567063885189905*^9}, {3.5670698395087633`*^9, 
  3.567069848260379*^9}, {3.567574529015103*^9, 3.567574529763904*^9}}],

Cell["\<\
(* bonding and lone pairs *)
Ekin = Ekin + \
2.25*(2*k4/NH^2+k5/CN^2+k6/CO^2+k7/CH^2+k8/NM^2+k9/OL^2+k10/OM^2+k11/OC^2+k12/\
OO^2);

vee=vee+3.0*(2*sig4/NH+sig5/CN+sig6/CO+sig7/CH+sig9/OL);\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668390630515456`*^9, 3.566839078760773*^9}, {
   3.566839318626795*^9, 3.5668393238528037`*^9}, {3.566839506201524*^9, 
   3.5668395108503323`*^9}, {3.566839541207986*^9, 3.566839568664034*^9}, {
   3.566839684650237*^9, 3.56683970887708*^9}, {3.5668397799352045`*^9, 
   3.5668398201988754`*^9}, {3.566839877872177*^9, 3.5668398909606*^9}, {
   3.5668404106911125`*^9, 3.566840417539525*^9}, {3.5668404815464373`*^9, 
   3.566840691460406*^9}, {3.566840732519678*^9, 3.5668408608923035`*^9}, {
   3.566840930000425*^9, 3.5668409412168446`*^9}, {3.566841374686006*^9, 
   3.566841375388007*^9}, {3.5668421835006266`*^9, 3.5668421845770283`*^9}, {
   3.566887398037596*^9, 3.566887402078003*^9}, {3.5668874322328568`*^9, 
   3.5668874455864797`*^9}, {3.566887493166563*^9, 3.5668874939309645`*^9}, {
   3.5668882488786907`*^9, 3.5668882511094947`*^9}, {3.56688847075788*^9, 
   3.566888479322295*^9}, {3.5668885903476906`*^9, 3.5668885954956994`*^9}, {
   3.566971858103361*^9, 3.566971918927868*^9}, {3.566987354128765*^9, 
   3.5669873624903803`*^9}, {3.566987488164201*^9, 3.566987535619484*^9}, {
   3.566987795843541*^9, 3.5669878015843515`*^9}, {3.566987856168847*^9, 
   3.56698786325126*^9}, {3.566988264733565*^9, 3.566988269819174*^9}, {
   3.566988366133743*^9, 3.5669883719525533`*^9}, {3.5669885431160536`*^9, 
   3.566988549309265*^9}, {3.566991946182942*^9, 3.5669919483981457`*^9}, 
   3.5669919866182127`*^9, {3.5669920226542764`*^9, 3.5669920505783253`*^9}, {
   3.566992103743219*^9, 3.5669921079396257`*^9}, 3.5669921962201815`*^9, {
   3.5669922303998413`*^9, 3.5669922340658474`*^9}, {3.566992291177548*^9, 
   3.5669922926751504`*^9}, {3.566992333890423*^9, 3.5669923348264246`*^9}, 
   3.566992409675356*^9, {3.566997472043352*^9, 3.566997587171554*^9}, {
   3.566997674297707*^9, 3.5669977903463106`*^9}, {3.566997868580448*^9, 
   3.5669978780652647`*^9}, {3.566998052270771*^9, 3.566998071786405*^9}, {
   3.5669981665409718`*^9, 3.5669981730305834`*^9}, {3.56699824522751*^9, 
   3.566998326394452*^9}, {3.5669983623837156`*^9, 3.5669984017425847`*^9}, {
   3.566998470788306*^9, 3.5669984830811276`*^9}, {3.566998534405218*^9, 
   3.5669985995821323`*^9}, {3.5669988153461113`*^9, 3.566998936995125*^9}, 
   3.566998967134378*^9, 3.5669990224832754`*^9, {3.5669990973790064`*^9, 
   3.56699912179305*^9}, {3.5669993047345705`*^9, 3.5669993241410046`*^9}, {
   3.566999416274767*^9, 3.566999426742385*^9}, {3.566999497301309*^9, 
   3.5669995032137194`*^9}, {3.566999549733001*^9, 3.566999551433404*^9}, {
   3.5670567651502*^9, 3.567056829110312*^9}, {3.5670569077344503`*^9, 
   3.5670570287126627`*^9}, {3.56705707303234*^9, 3.5670571842137356`*^9}, {
   3.5670572582982655`*^9, 3.567057292103525*^9}, {3.56705732891959*^9, 
   3.5670574906294737`*^9}, {3.5670578301328697`*^9, 3.56705795512029*^9}, {
   3.5670580197668033`*^9, 3.567058030437222*^9}, {3.5670580810281105`*^9, 
   3.5670581611654515`*^9}, {3.5670581932547083`*^9, 3.567058240616391*^9}, {
   3.5670582801312604`*^9, 3.5670582990540934`*^9}, {3.567058343560972*^9, 
   3.5670584018894744`*^9}, {3.56705845049916*^9, 3.5670585083128614`*^9}, {
   3.567058578731385*^9, 3.567058593645011*^9}, {3.5670586348134837`*^9, 
   3.5670586577143235`*^9}, {3.567058689007979*^9, 3.5670587629209085`*^9}, {
   3.567058832777831*^9, 3.5670588394702425`*^9}, {3.567058955628047*^9, 
   3.5670589616496572`*^9}, {3.5670592284569263`*^9, 3.567059354442747*^9}, {
   3.5670593898548098`*^9, 3.567059560456709*^9}, {3.5670596222796183`*^9, 
   3.5670596362884426`*^9}, {3.5670596726677065`*^9, 3.567059775347087*^9}, {
   3.567059838449198*^9, 3.5670598453756094`*^9}, {3.5670598764976645`*^9, 
   3.5670598805380716`*^9}, 3.567059911254525*^9, {3.5670599744502363`*^9, 
   3.567060007818695*^9}, {3.567060059127185*^9, 3.5670600818564253`*^9}, {
   3.5670601143824825`*^9, 3.567060314593234*^9}, {3.5670607181347427`*^9, 
   3.5670607311607656`*^9}, {3.567060785573661*^9, 3.567061232998047*^9}, {
   3.5670626282800975`*^9, 3.56706267537658*^9}, 3.567062795371991*^9, {
   3.567069851801585*^9, 3.567069879460434*^9}, {3.567070767429593*^9, 
   3.5670707823588195`*^9}, {3.5670805354535503`*^9, 
   3.5670806312689185`*^9}, {3.567083393113369*^9, 3.5670834586022844`*^9}, {
   3.567083501299559*^9, 3.5670835804696984`*^9}, {3.5670836156477604`*^9, 
   3.5670836215757704`*^9}, {3.567083656504232*^9, 3.5670836570034328`*^9}, 
   3.5670837376399746`*^9, {3.567247718249197*^9, 3.567247775142497*^9}}],

Cell["\<\
(*  1  2  3  4  5  6  7  8  9  10 11   12   13   14   15 16 17  *)
(*  C  O  O  Hn Hn Hc B  B  L2 L1 L1   Bp   Bp   Bp   Bp L1 L1  *)
as={\" \",C,O,O,H,L,H};    (* atom string *)  
rr={Cr,Nr,Oa,NH,NH,CH,CN,CO,OL,NM,NM,  OC,  OC,  OO,  OO,OM,OM};
oc={-2,-2,-2,-2,-2,-2,-2,-2,-2,-1,-1,-1/2,-1/2,-1/2,-1/2,-1,-1};
ch={ZC,ZN,ZO, 1, 0, 1, 0, 0, 0, 0, 0,   0,   0,   0,   0, 0, 0};
csw1=Cos[w1/2];  ssw1=Sin[w1/2];
csw2=Cos[w2];    ssw2=Sin[w2];
csw3=Cos[w3];    ssw3=Sin[w3];\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668390630515456`*^9, 3.566839078760773*^9}, {
   3.566839318626795*^9, 3.5668393238528037`*^9}, {3.566839506201524*^9, 
   3.5668395108503323`*^9}, {3.566839541207986*^9, 3.566839568664034*^9}, {
   3.566839684650237*^9, 3.56683970887708*^9}, {3.5668397799352045`*^9, 
   3.5668398201988754`*^9}, {3.566839877872177*^9, 3.5668398909606*^9}, {
   3.5668404106911125`*^9, 3.566840417539525*^9}, {3.5668404815464373`*^9, 
   3.566840691460406*^9}, {3.566840732519678*^9, 3.5668408608923035`*^9}, {
   3.566840930000425*^9, 3.5668409412168446`*^9}, {3.566841374686006*^9, 
   3.566841375388007*^9}, {3.5668421835006266`*^9, 3.5668421845770283`*^9}, {
   3.566887398037596*^9, 3.566887402078003*^9}, {3.5668874322328568`*^9, 
   3.5668874455864797`*^9}, {3.566887493166563*^9, 3.5668874939309645`*^9}, {
   3.5668882488786907`*^9, 3.5668882511094947`*^9}, {3.56688847075788*^9, 
   3.566888479322295*^9}, {3.5668885903476906`*^9, 3.5668885954956994`*^9}, {
   3.566971858103361*^9, 3.566971918927868*^9}, {3.566987354128765*^9, 
   3.5669873624903803`*^9}, {3.566987488164201*^9, 3.566987535619484*^9}, {
   3.566987795843541*^9, 3.5669878015843515`*^9}, {3.566987856168847*^9, 
   3.56698786325126*^9}, {3.566988264733565*^9, 3.566988269819174*^9}, {
   3.566988366133743*^9, 3.5669883719525533`*^9}, {3.5669885431160536`*^9, 
   3.566988549309265*^9}, {3.566991946182942*^9, 3.5669919483981457`*^9}, 
   3.5669919866182127`*^9, {3.5669920226542764`*^9, 3.5669920505783253`*^9}, {
   3.566992103743219*^9, 3.5669921079396257`*^9}, 3.5669921962201815`*^9, {
   3.5669922303998413`*^9, 3.5669922340658474`*^9}, {3.566992291177548*^9, 
   3.5669922926751504`*^9}, {3.566992333890423*^9, 3.5669923348264246`*^9}, 
   3.566992409675356*^9, {3.566997472043352*^9, 3.566997587171554*^9}, {
   3.566997674297707*^9, 3.5669977903463106`*^9}, {3.566997868580448*^9, 
   3.5669978780652647`*^9}, {3.566998052270771*^9, 3.566998071786405*^9}, {
   3.5669981665409718`*^9, 3.5669981730305834`*^9}, {3.56699824522751*^9, 
   3.566998326394452*^9}, {3.5669983623837156`*^9, 3.5669984017425847`*^9}, {
   3.566998470788306*^9, 3.5669984830811276`*^9}, {3.566998534405218*^9, 
   3.5669985995821323`*^9}, {3.5669988153461113`*^9, 3.566998936995125*^9}, 
   3.566998967134378*^9, 3.5669990224832754`*^9, {3.5669990973790064`*^9, 
   3.56699912179305*^9}, {3.5669993047345705`*^9, 3.5669993241410046`*^9}, {
   3.566999416274767*^9, 3.566999426742385*^9}, {3.566999497301309*^9, 
   3.5669995032137194`*^9}, {3.566999549733001*^9, 3.566999551433404*^9}, {
   3.5670567651502*^9, 3.567056829110312*^9}, {3.5670569077344503`*^9, 
   3.5670570287126627`*^9}, {3.56705707303234*^9, 3.5670571842137356`*^9}, {
   3.5670572582982655`*^9, 3.567057292103525*^9}, {3.56705732891959*^9, 
   3.5670574906294737`*^9}, {3.5670578301328697`*^9, 3.56705795512029*^9}, {
   3.5670580197668033`*^9, 3.567058030437222*^9}, {3.5670580810281105`*^9, 
   3.5670581611654515`*^9}, {3.5670581932547083`*^9, 3.567058240616391*^9}, {
   3.5670582801312604`*^9, 3.5670582990540934`*^9}, {3.567058343560972*^9, 
   3.5670584018894744`*^9}, {3.56705845049916*^9, 3.5670585083128614`*^9}, {
   3.567058578731385*^9, 3.567058593645011*^9}, {3.5670586348134837`*^9, 
   3.5670586577143235`*^9}, {3.567058689007979*^9, 3.5670587629209085`*^9}, {
   3.567058832777831*^9, 3.5670588394702425`*^9}, {3.567058955628047*^9, 
   3.5670589616496572`*^9}, {3.5670592284569263`*^9, 3.567059354442747*^9}, {
   3.5670593898548098`*^9, 3.567059560456709*^9}, {3.5670596222796183`*^9, 
   3.5670596362884426`*^9}, {3.5670596726677065`*^9, 3.567059775347087*^9}, {
   3.567059838449198*^9, 3.5670598453756094`*^9}, {3.5670598764976645`*^9, 
   3.5670598805380716`*^9}, 3.567059911254525*^9, {3.5670599744502363`*^9, 
   3.567060007818695*^9}, {3.567060059127185*^9, 3.5670600818564253`*^9}, {
   3.5670601143824825`*^9, 3.567060314593234*^9}, {3.5670607181347427`*^9, 
   3.5670607311607656`*^9}, {3.567060785573661*^9, 3.567061232998047*^9}, {
   3.5670626282800975`*^9, 3.56706267537658*^9}, 3.567062795371991*^9, {
   3.567069851801585*^9, 3.567069879460434*^9}, {3.567070767429593*^9, 
   3.5670707823588195`*^9}, {3.5670805354535503`*^9, 
   3.5670806312689185`*^9}, {3.567083393113369*^9, 3.5670834586022844`*^9}, {
   3.567083501299559*^9, 3.5670835804696984`*^9}, {3.5670836156477604`*^9, 
   3.5670836215757704`*^9}, {3.567083656504232*^9, 3.5670836570034328`*^9}, 
   3.5670837376399746`*^9, {3.567228361875599*^9, 3.5672284052748756`*^9}, {
   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.567510844107608*^9, 3.5675108447316084`*^9}, {
   3.5675109368653708`*^9, 3.567510937832572*^9}, 3.5675112986612062`*^9, 
   3.5675114921015453`*^9, {3.5675120605509443`*^9, 3.5675120611437454`*^9}, {
   3.5675140171219807`*^9, 3.567514049944438*^9}, {3.5680997047417803`*^9, 
   3.5680997073157854`*^9}}],

Cell["\<\
dcnx=Cr+Nr+2*CN;             dcox=(Cr+Oa+2*CO)*csw2;  dnhx=dcnx+(Nr+NH)*csw1;
dchx=(Cr+CH)*csw3;           dccnx=Cr+CN; 
pnhx=dcnx+(Nr+NH+pnh)*csw1;  pchx=(Cr+CH+pch)*csw3;
dcoy=(Cr+Oa+2*CO)*ssw2;      dnhy=(Nr+NH)*ssw1;       dchy=(Cr+CH)*ssw3;  
pnhy=(Nr+NH+pnh)*ssw1;       pchy=(Cr+CH+pch)*ssw3;
dccox=(Cr+CO)*csw2;          dccoy=(Cr+CO)*ssw2;      
dolx=(Cr+2*CO+2*Oa+OL)*csw2; doly=(Cr+2*CO+2*Oa+OL)*ssw2;
dom1x=dcox-(Oa+OM)*ssw2;     dom2x=dcox+(Oa+OM)*ssw2;
dom1y=dcoy+(Oa+OM)*csw2;     dom2y=dcoy-(Oa+OM)*csw2;\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668390630515456`*^9, 3.566839078760773*^9}, {
   3.566839318626795*^9, 3.5668393238528037`*^9}, {3.566839506201524*^9, 
   3.5668395108503323`*^9}, {3.566839541207986*^9, 3.566839568664034*^9}, {
   3.566839684650237*^9, 3.56683970887708*^9}, {3.5668397799352045`*^9, 
   3.5668398201988754`*^9}, {3.566839877872177*^9, 3.5668398909606*^9}, {
   3.5668404106911125`*^9, 3.566840417539525*^9}, {3.5668404815464373`*^9, 
   3.566840691460406*^9}, {3.566840732519678*^9, 3.5668408608923035`*^9}, {
   3.566840930000425*^9, 3.5668409412168446`*^9}, {3.566841374686006*^9, 
   3.566841375388007*^9}, {3.5668421835006266`*^9, 3.5668421845770283`*^9}, {
   3.566887398037596*^9, 3.566887402078003*^9}, {3.5668874322328568`*^9, 
   3.5668874455864797`*^9}, {3.566887493166563*^9, 3.5668874939309645`*^9}, {
   3.5668882488786907`*^9, 3.5668882511094947`*^9}, {3.56688847075788*^9, 
   3.566888479322295*^9}, {3.5668885903476906`*^9, 3.5668885954956994`*^9}, {
   3.566971858103361*^9, 3.566971918927868*^9}, {3.566987354128765*^9, 
   3.5669873624903803`*^9}, {3.566987488164201*^9, 3.566987535619484*^9}, {
   3.566987795843541*^9, 3.5669878015843515`*^9}, {3.566987856168847*^9, 
   3.56698786325126*^9}, {3.566988264733565*^9, 3.566988269819174*^9}, {
   3.566988366133743*^9, 3.5669883719525533`*^9}, {3.5669885431160536`*^9, 
   3.566988549309265*^9}, {3.566991946182942*^9, 3.5669919483981457`*^9}, 
   3.5669919866182127`*^9, {3.5669920226542764`*^9, 3.5669920505783253`*^9}, {
   3.566992103743219*^9, 3.5669921079396257`*^9}, 3.5669921962201815`*^9, {
   3.5669922303998413`*^9, 3.5669922340658474`*^9}, {3.566992291177548*^9, 
   3.5669922926751504`*^9}, {3.566992333890423*^9, 3.5669923348264246`*^9}, 
   3.566992409675356*^9, {3.566997472043352*^9, 3.566997587171554*^9}, {
   3.566997674297707*^9, 3.5669977903463106`*^9}, {3.566997868580448*^9, 
   3.5669978780652647`*^9}, {3.566998052270771*^9, 3.566998071786405*^9}, {
   3.5669981665409718`*^9, 3.5669981730305834`*^9}, {3.56699824522751*^9, 
   3.566998326394452*^9}, {3.5669983623837156`*^9, 3.5669984017425847`*^9}, {
   3.566998470788306*^9, 3.5669984830811276`*^9}, {3.566998534405218*^9, 
   3.5669985995821323`*^9}, {3.5669988153461113`*^9, 3.566998936995125*^9}, 
   3.566998967134378*^9, 3.5669990224832754`*^9, {3.5669990973790064`*^9, 
   3.56699912179305*^9}, {3.5669993047345705`*^9, 3.5669993241410046`*^9}, {
   3.566999416274767*^9, 3.566999426742385*^9}, {3.566999497301309*^9, 
   3.5669995032137194`*^9}, {3.566999549733001*^9, 3.566999551433404*^9}, {
   3.5670567651502*^9, 3.567056829110312*^9}, {3.5670569077344503`*^9, 
   3.5670570287126627`*^9}, {3.56705707303234*^9, 3.5670571842137356`*^9}, {
   3.5670572582982655`*^9, 3.567057292103525*^9}, {3.56705732891959*^9, 
   3.5670574906294737`*^9}, {3.5670578301328697`*^9, 3.56705795512029*^9}, {
   3.5670580197668033`*^9, 3.567058030437222*^9}, {3.5670580810281105`*^9, 
   3.5670581611654515`*^9}, {3.5670581932547083`*^9, 3.567058240616391*^9}, {
   3.5670582801312604`*^9, 3.5670582990540934`*^9}, {3.567058343560972*^9, 
   3.5670584018894744`*^9}, {3.56705845049916*^9, 3.5670585083128614`*^9}, {
   3.567058578731385*^9, 3.567058593645011*^9}, {3.5670586348134837`*^9, 
   3.5670586577143235`*^9}, {3.567058689007979*^9, 3.5670587629209085`*^9}, {
   3.567058832777831*^9, 3.5670588394702425`*^9}, {3.567058955628047*^9, 
   3.5670589616496572`*^9}, {3.5670592284569263`*^9, 3.567059354442747*^9}, {
   3.5670593898548098`*^9, 3.567059560456709*^9}, {3.5670596222796183`*^9, 
   3.5670596362884426`*^9}, {3.5670596726677065`*^9, 3.567059775347087*^9}, {
   3.567059838449198*^9, 3.5670598453756094`*^9}, {3.5670598764976645`*^9, 
   3.5670598805380716`*^9}, 3.567059911254525*^9, {3.5670599744502363`*^9, 
   3.567060007818695*^9}, {3.567060059127185*^9, 3.5670600818564253`*^9}, {
   3.5670601143824825`*^9, 3.567060314593234*^9}, {3.5670607181347427`*^9, 
   3.5670607311607656`*^9}, {3.567060785573661*^9, 3.567061232998047*^9}, {
   3.5670626282800975`*^9, 3.56706267537658*^9}, 3.567062795371991*^9, {
   3.567069851801585*^9, 3.567069879460434*^9}, {3.567070767429593*^9, 
   3.5670707823588195`*^9}, {3.5670805354535503`*^9, 
   3.5670806312689185`*^9}, {3.567083393113369*^9, 3.5670834586022844`*^9}, {
   3.567083501299559*^9, 3.5670835804696984`*^9}, {3.5670836156477604`*^9, 
   3.5670836215757704`*^9}, {3.567083656504232*^9, 3.5670836570034328`*^9}, {
   3.5670837376399746`*^9, 3.56708374668799*^9}, {3.5673117165829782`*^9, 
   3.5673117180649805`*^9}}],

Cell["\<\
(* 6 nuclei: coord *)
xn={0,-dcnx,dcox,-pnhx,-pnhx,pchx};
yn={0,0,dcoy,pnhy,-pnhy,-pchy};
zn={0,0,0,0,0,0};\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668390630515456`*^9, 3.566839078760773*^9}, {
   3.566839318626795*^9, 3.5668393238528037`*^9}, {3.566839506201524*^9, 
   3.5668395108503323`*^9}, {3.566839541207986*^9, 3.566839568664034*^9}, {
   3.566839684650237*^9, 3.56683970887708*^9}, {3.5668397799352045`*^9, 
   3.5668398201988754`*^9}, {3.566839877872177*^9, 3.5668398909606*^9}, {
   3.5668404106911125`*^9, 3.566840417539525*^9}, {3.5668404815464373`*^9, 
   3.566840691460406*^9}, {3.566840732519678*^9, 3.5668408608923035`*^9}, {
   3.566840930000425*^9, 3.5668409412168446`*^9}, {3.566841374686006*^9, 
   3.566841375388007*^9}, {3.5668421835006266`*^9, 3.5668421845770283`*^9}, {
   3.566887398037596*^9, 3.566887402078003*^9}, {3.5668874322328568`*^9, 
   3.5668874455864797`*^9}, {3.566887493166563*^9, 3.5668874939309645`*^9}, {
   3.5668882488786907`*^9, 3.5668882511094947`*^9}, {3.56688847075788*^9, 
   3.566888479322295*^9}, {3.5668885903476906`*^9, 3.5668885954956994`*^9}, {
   3.566971858103361*^9, 3.566971918927868*^9}, {3.566987354128765*^9, 
   3.5669873624903803`*^9}, {3.566987488164201*^9, 3.566987535619484*^9}, {
   3.566987795843541*^9, 3.5669878015843515`*^9}, {3.566987856168847*^9, 
   3.56698786325126*^9}, {3.566988264733565*^9, 3.566988269819174*^9}, {
   3.566988366133743*^9, 3.5669883719525533`*^9}, {3.5669885431160536`*^9, 
   3.566988549309265*^9}, {3.566991946182942*^9, 3.5669919483981457`*^9}, 
   3.5669919866182127`*^9, {3.5669920226542764`*^9, 3.5669920505783253`*^9}, {
   3.566992103743219*^9, 3.5669921079396257`*^9}, 3.5669921962201815`*^9, {
   3.5669922303998413`*^9, 3.5669922340658474`*^9}, {3.566992291177548*^9, 
   3.5669922926751504`*^9}, {3.566992333890423*^9, 3.5669923348264246`*^9}, 
   3.566992409675356*^9, {3.566997472043352*^9, 3.566997587171554*^9}, {
   3.566997674297707*^9, 3.5669977903463106`*^9}, {3.566997868580448*^9, 
   3.5669978780652647`*^9}, {3.566998052270771*^9, 3.566998071786405*^9}, {
   3.5669981665409718`*^9, 3.5669981730305834`*^9}, {3.56699824522751*^9, 
   3.566998326394452*^9}, {3.5669983623837156`*^9, 3.5669984017425847`*^9}, {
   3.566998470788306*^9, 3.5669984830811276`*^9}, {3.566998534405218*^9, 
   3.5669985995821323`*^9}, {3.5669988153461113`*^9, 3.566998936995125*^9}, 
   3.566998967134378*^9, 3.5669990224832754`*^9, {3.5669990973790064`*^9, 
   3.56699912179305*^9}, {3.5669993047345705`*^9, 3.5669993241410046`*^9}, {
   3.566999416274767*^9, 3.566999426742385*^9}, {3.566999497301309*^9, 
   3.5669995032137194`*^9}, {3.566999549733001*^9, 3.566999551433404*^9}, {
   3.5670567651502*^9, 3.567056829110312*^9}, {3.5670569077344503`*^9, 
   3.5670570287126627`*^9}, {3.56705707303234*^9, 3.5670571842137356`*^9}, {
   3.5670572582982655`*^9, 3.567057292103525*^9}, {3.56705732891959*^9, 
   3.5670574906294737`*^9}, {3.5670578301328697`*^9, 3.56705795512029*^9}, {
   3.5670580197668033`*^9, 3.567058030437222*^9}, {3.5670580810281105`*^9, 
   3.5670581611654515`*^9}, {3.5670581932547083`*^9, 3.567058240616391*^9}, {
   3.5670582801312604`*^9, 3.5670582990540934`*^9}, {3.567058343560972*^9, 
   3.5670584018894744`*^9}, {3.56705845049916*^9, 3.5670585083128614`*^9}, {
   3.567058578731385*^9, 3.567058593645011*^9}, {3.5670586348134837`*^9, 
   3.5670586577143235`*^9}, {3.567058689007979*^9, 3.5670587629209085`*^9}, {
   3.567058832777831*^9, 3.5670588394702425`*^9}, {3.567058955628047*^9, 
   3.5670589616496572`*^9}, {3.5670592284569263`*^9, 3.567059354442747*^9}, {
   3.5670593898548098`*^9, 3.567059560456709*^9}, {3.5670596222796183`*^9, 
   3.5670596362884426`*^9}, {3.5670596726677065`*^9, 3.567059775347087*^9}, {
   3.567059838449198*^9, 3.5670598453756094`*^9}, {3.5670598764976645`*^9, 
   3.5670598805380716`*^9}, 3.567059911254525*^9, {3.5670599744502363`*^9, 
   3.567060007818695*^9}, {3.567060059127185*^9, 3.5670600818564253`*^9}, {
   3.5670601143824825`*^9, 3.567060314593234*^9}, {3.5670607181347427`*^9, 
   3.5670607311607656`*^9}, {3.567060785573661*^9, 3.567061232998047*^9}, {
   3.5670626282800975`*^9, 3.56706267537658*^9}, 3.567062795371991*^9, {
   3.567069851801585*^9, 3.567069879460434*^9}, {3.567070767429593*^9, 
   3.5670707823588195`*^9}, {3.5670805354535503`*^9, 
   3.5670806312689185`*^9}, {3.567083393113369*^9, 3.5670834586022844`*^9}, {
   3.567083501299559*^9, 3.5670835804696984`*^9}, {3.5670836156477604`*^9, 
   3.5670836215757704`*^9}, {3.567083656504232*^9, 3.5670836570034328`*^9}, {
   3.5670837376399746`*^9, 3.567083757514409*^9}}],

Cell["\<\
(* 17 clouds: coord *)        
xc={0,-dcnx,dcox,-dnhx,-dnhx,dchx,-dccnx,dccox,dolx,-dcnx,-dcnx,0,0,dcox,dcox,\
dom1x,dom2x};
yc={0,0,dcoy,dnhy,-dnhy,-dchy,0,dccoy,doly,0,0,0,0,dcoy,dcoy,dom1y,dom2y};
zc={0,0,0,0,0,0,0,0,0,NM,-NM,OC,-OC,OO,-OO,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.5672479122355375`*^9, 
   3.567247913155939*^9}, {3.5672481119782887`*^9, 3.567248113819092*^9}, {
   3.567248190337226*^9, 3.567248240850115*^9}}],

Cell["\<\
(* epot terms *)
vne=vne-(3-(pnh/NH)^2)/NH;
vne=vne-(3-(pch/CH)^2)/CH;
vnn = 0.0;\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668390630515456`*^9, 3.566839078760773*^9}, {
   3.566839318626795*^9, 3.5668393238528037`*^9}, {3.566839506201524*^9, 
   3.5668395108503323`*^9}, {3.566839541207986*^9, 3.566839568664034*^9}, {
   3.566839684650237*^9, 3.56683970887708*^9}, {3.5668397799352045`*^9, 
   3.5668398201988754`*^9}, {3.566839877872177*^9, 3.5668398909606*^9}, {
   3.5668404106911125`*^9, 3.566840417539525*^9}, {3.5668404815464373`*^9, 
   3.566840691460406*^9}, {3.566840732519678*^9, 3.5668408608923035`*^9}, {
   3.566840930000425*^9, 3.5668409412168446`*^9}, {3.566841374686006*^9, 
   3.566841375388007*^9}, {3.5668421835006266`*^9, 3.5668421845770283`*^9}, {
   3.566887398037596*^9, 3.566887402078003*^9}, {3.5668874322328568`*^9, 
   3.5668874455864797`*^9}, {3.566887493166563*^9, 3.5668874939309645`*^9}, {
   3.5668882488786907`*^9, 3.5668882511094947`*^9}, {3.56688847075788*^9, 
   3.566888479322295*^9}, {3.5668885903476906`*^9, 3.5668885954956994`*^9}, {
   3.566971858103361*^9, 3.566971918927868*^9}, {3.566987354128765*^9, 
   3.5669873624903803`*^9}, {3.566987488164201*^9, 3.566987535619484*^9}, {
   3.566987795843541*^9, 3.5669878015843515`*^9}, {3.566987856168847*^9, 
   3.56698786325126*^9}, {3.566988264733565*^9, 3.566988269819174*^9}, {
   3.566988366133743*^9, 3.5669883719525533`*^9}, {3.5669885431160536`*^9, 
   3.566988549309265*^9}, {3.566991946182942*^9, 3.5669919483981457`*^9}, 
   3.5669919866182127`*^9, {3.5669920226542764`*^9, 3.5669920505783253`*^9}, {
   3.566992103743219*^9, 3.5669921079396257`*^9}, 3.5669921962201815`*^9, {
   3.5669922303998413`*^9, 3.5669922340658474`*^9}, {3.566992291177548*^9, 
   3.5669922926751504`*^9}, {3.566992333890423*^9, 3.5669923348264246`*^9}, 
   3.566992409675356*^9, {3.566997472043352*^9, 3.566997587171554*^9}, {
   3.566997674297707*^9, 3.5669977903463106`*^9}, {3.566997868580448*^9, 
   3.5669978780652647`*^9}, {3.566998052270771*^9, 3.566998071786405*^9}, {
   3.5669981665409718`*^9, 3.5669981730305834`*^9}, {3.56699824522751*^9, 
   3.566998326394452*^9}, {3.5669983623837156`*^9, 3.5669984017425847`*^9}, {
   3.566998470788306*^9, 3.5669984830811276`*^9}, {3.566998534405218*^9, 
   3.5669985995821323`*^9}, {3.5669988153461113`*^9, 3.566998936995125*^9}, 
   3.566998967134378*^9, 3.5669990224832754`*^9, {3.5669990973790064`*^9, 
   3.56699912179305*^9}, {3.5669993047345705`*^9, 3.5669993241410046`*^9}, {
   3.566999416274767*^9, 3.566999426742385*^9}, {3.566999497301309*^9, 
   3.5669995032137194`*^9}, {3.566999549733001*^9, 3.566999551433404*^9}, {
   3.5670567651502*^9, 3.567056829110312*^9}, {3.5670569077344503`*^9, 
   3.5670570287126627`*^9}, {3.56705707303234*^9, 3.5670571842137356`*^9}, {
   3.5670572582982655`*^9, 3.567057292103525*^9}, {3.56705732891959*^9, 
   3.5670574906294737`*^9}, {3.5670578301328697`*^9, 3.56705795512029*^9}, {
   3.5670580197668033`*^9, 3.567058030437222*^9}, {3.5670580810281105`*^9, 
   3.5670581611654515`*^9}, {3.5670581932547083`*^9, 3.567058240616391*^9}, {
   3.5670582801312604`*^9, 3.5670582990540934`*^9}, {3.567058343560972*^9, 
   3.5670584018894744`*^9}, {3.56705845049916*^9, 3.5670585083128614`*^9}, {
   3.567058578731385*^9, 3.567058593645011*^9}, {3.5670586348134837`*^9, 
   3.5670586577143235`*^9}, {3.567058689007979*^9, 3.5670587629209085`*^9}, {
   3.567058832777831*^9, 3.5670588394702425`*^9}, {3.567058955628047*^9, 
   3.5670589616496572`*^9}, {3.5670592284569263`*^9, 3.567059354442747*^9}, {
   3.5670593898548098`*^9, 3.567059560456709*^9}, {3.5670596222796183`*^9, 
   3.5670596362884426`*^9}, {3.5670596726677065`*^9, 3.567059775347087*^9}, {
   3.567059838449198*^9, 3.5670598453756094`*^9}, {3.5670598764976645`*^9, 
   3.5670598805380716`*^9}, 3.567059911254525*^9, {3.5670599744502363`*^9, 
   3.567060007818695*^9}, {3.567060059127185*^9, 3.5670600818564253`*^9}, {
   3.5670601143824825`*^9, 3.567060314593234*^9}, {3.5670607181347427`*^9, 
   3.5670607311607656`*^9}, {3.567060785573661*^9, 3.567061232998047*^9}, {
   3.5670626282800975`*^9, 3.56706267537658*^9}, 3.567062795371991*^9, {
   3.567069851801585*^9, 3.567069879460434*^9}, {3.567070767429593*^9, 
   3.5670707823588195`*^9}, {3.5670805354535503`*^9, 
   3.5670806312689185`*^9}, {3.567083393113369*^9, 3.5670834586022844`*^9}, {
   3.567083501299559*^9, 3.5670835804696984`*^9}, {3.5670836156477604`*^9, 
   3.5670836215757704`*^9}, {3.567083656504232*^9, 3.5670836570034328`*^9}, {
   3.5670837376399746`*^9, 3.5670837651584225`*^9}, {3.5675138108896184`*^9, 
   3.567513811342019*^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.125/OC-0.125/OO;  (* selfinteraction corrected *)\
\>", "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.5681008975666757`*^9, 
  3.568100982477625*^9}}],

Cell["\<\
(*nn*)
For[i = 1, i < nn, i++,
  For[j = i+1, j < nn+1, j++,
    vnn = vnn + \
ch[[i]]*ch[[j]]/Sqrt[(xn[[i]]-xn[[j]])^2+(yn[[i]]-yn[[j]])^2+(zn[[i]]-zn[[j]])\
^2]]]\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668390630515456`*^9, 3.5668390890567913`*^9}, {
  3.566841524165468*^9, 3.566841530218279*^9}, {3.5670590745470557`*^9, 
  3.567059089304682*^9}}],

Cell["\<\
(*nw*)
For[i = 1, i < nc+1, i++,
  For[j = 1, j < nn+1, j++,
    If[i \[NotEqual] j,
      vne = vne + \
oc[[i]]*ch[[j]]/Sqrt[(xc[[i]]-xn[[j]])^2+(yc[[i]]-yn[[j]])^2+(zc[[i]]-zn[[j]])\
^2]]]]\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668390630515456`*^9, 3.5668390979332066`*^9}, 
   3.5668412394337683`*^9, {3.566841289743857*^9, 3.566841311927096*^9}, 
   3.5668415580955276`*^9, {3.5668416231788425`*^9, 3.5668416346760626`*^9}, {
   3.5669003862430363`*^9, 3.5669003862430363`*^9}, {3.5670591033447065`*^9, 
   3.5670591078063145`*^9}}],

Cell[CellGroupData[{

Cell["\<\
Epot=vne+vee+vnn;
func=Ekin+Epot;
t = FindMinimum[func,{Cr,0.26},{Nr,0.19},{Oa,0.19},{NH,0.88},{CH,1.25},{CO,0.\
92},{CN,0.97},{OL,0.905},{NM,0.93},{OC,1.45},{OO,1.22},{OM,0.85},{w1,2.28},{\
w2,1.197},{w3,0.85},{pch,0.58},{pnh,0.32},{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.5675105559751015`*^9, 
   3.5675106443492565`*^9}, {3.5675109760214396`*^9, 3.567511011776702*^9}, {
   3.567511104113264*^9, 3.567511174110587*^9}, {3.567511372371335*^9, 
   3.567511453023477*^9}, {3.5675122037435956`*^9, 3.5675122283916388`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "188.55380346547008`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"Cr", "\[Rule]", "0.25559559017708283`"}], ",", 
     RowBox[{"Nr", "\[Rule]", "0.18966729854432163`"}], ",", 
     RowBox[{"Oa", "\[Rule]", "0.18771946482544016`"}], ",", 
     RowBox[{"NH", "\[Rule]", "1.0383282480711433`"}], ",", 
     RowBox[{"CH", "\[Rule]", "1.237275196687085`"}], ",", 
     RowBox[{"CO", "\[Rule]", "0.8999252425636789`"}], ",", 
     RowBox[{"CN", "\[Rule]", "1.0264269313212875`"}], ",", 
     RowBox[{"OL", "\[Rule]", "0.9157712114866523`"}], ",", 
     RowBox[{"NM", "\[Rule]", "0.8768151106038778`"}], ",", 
     RowBox[{"OC", "\[Rule]", "1.3099009653202152`"}], ",", 
     RowBox[{"OO", "\[Rule]", "1.1059447491867593`"}], ",", 
     RowBox[{"OM", "\[Rule]", "0.8572259933635882`"}], ",", 
     RowBox[{"w1", "\[Rule]", "2.2645186072642343`"}], ",", 
     RowBox[{"w2", "\[Rule]", "1.2135867010034747`"}], ",", 
     RowBox[{"w3", "\[Rule]", "0.8484612239084653`"}], ",", 
     RowBox[{"pch", "\[Rule]", "0.562749228287675`"}], ",", 
     RowBox[{"pnh", "\[Rule]", "0.5659407637751243`"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.567346286077984*^9, 3.5675104559789257`*^9, 3.5675104963829966`*^9, 
   3.5675106613688865`*^9, 3.5675108565096292`*^9, 3.5675109484093904`*^9, 
   3.56751103923275*^9, 3.567511182456602*^9, 3.5675113144328337`*^9, 
   3.567511462258693*^9, 3.5675114994959583`*^9, 3.5675117143083363`*^9, 
   3.567511977028397*^9, {3.567512071954564*^9, 3.5675120852613873`*^9}, 
   3.5675122599192944`*^9, 3.5675138867993517`*^9, 3.567514066121667*^9, 
   3.5675161988990145`*^9, 3.56751625919312*^9, 3.567516295634784*^9, 
   3.567516359860097*^9, {3.5675164066133795`*^9, 3.567516477297103*^9}, 
   3.567516535578806*^9, 3.5675165969181137`*^9, 3.56751662914777*^9, 
   3.567516661517827*^9, 3.5675167055567045`*^9, {3.567516739595964*^9, 
   3.5675167661160107`*^9}, {3.567516802822875*^9, 3.567516829077721*^9}, {
   3.5675168733817987`*^9, 3.5675169171554756`*^9}, 3.56753241237069*^9, {
   3.5675324502475567`*^9, 3.5675324800436087`*^9}, 3.5680997197802067`*^9, 
   3.5681010076248693`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Ekin /. t[[2]]
vne /. t[[2]]
vee /. t[[2]]
vnn /. t[[2]]
-Epot/Ekin /. t[[2]]
w1*rad /. t[[2]]
w2*rad /. t[[2]]
w3*rad /. t[[2]]
(Cr+Oa+2*CO)*bohr /. t[[2]]
(Cr+Nr+2*CN)*bohr /. t[[2]]
(Cr+CH+pch)*bohr /. t[[2]]
(Nr+NH+pnh)*bohr /. t[[2]]\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668390630515456`*^9, 3.566839112846833*^9}, {
   3.566842810779729*^9, 3.566842811091729*^9}, {3.5669879418285975`*^9, 
   3.5669879480686083`*^9}, {3.5670621290636206`*^9, 3.567062152463662*^9}, {
   3.567063905688341*^9, 3.5670639108675504`*^9}, {3.5670689307603674`*^9, 
   3.5670689312127686`*^9}, {3.5670696234639845`*^9, 
   3.5670696694216647`*^9}, {3.5670795267235785`*^9, 
   3.5670795309823856`*^9}, {3.567079591604092*^9, 3.567079592524494*^9}, {
   3.567079628810158*^9, 3.567079717433913*^9}, {3.567082843118803*^9, 
   3.567082848547613*^9}, {3.567082895987296*^9, 3.567082897079298*^9}, {
   3.567085039347061*^9, 3.5670851279708166`*^9}, 3.567345843318006*^9}],

Cell[BoxData["188.5538034744355`"], "Output",
 CellChangeTimes->{
  3.5673462861091833`*^9, 3.567510456010126*^9, 3.567510496398597*^9, 
   3.567510661384486*^9, 3.5675108565252295`*^9, 3.5675109484249907`*^9, 
   3.5675110392639503`*^9, 3.567511182472202*^9, 3.567511314448434*^9, 
   3.567511462274293*^9, 3.5675114995115585`*^9, 3.567511714323936*^9, 
   3.5675119770439973`*^9, {3.5675120719701643`*^9, 3.5675120852769876`*^9}, 
   3.5675122599504943`*^9, 3.5675138868305516`*^9, 3.567514066137267*^9, 
   3.5675161989146147`*^9, 3.56751625920872*^9, 3.567516295650384*^9, 
   3.567516359875697*^9, {3.567516406628979*^9, 3.567516477312703*^9}, 
   3.5675165355944057`*^9, 3.5675165969337134`*^9, 3.56751662916337*^9, 
   3.567516661533427*^9, 3.5675167055723042`*^9, {3.5675167396115637`*^9, 
   3.567516766131611*^9}, {3.567516802838475*^9, 3.567516829093321*^9}, {
   3.567516873397399*^9, 3.567516917171076*^9}, 3.56753241238629*^9, {
   3.5675324502631564`*^9, 3.567532480059209*^9}, 3.568099719811407*^9, 
   3.568101007640469*^9}],

Cell[BoxData[
 RowBox[{"-", "585.8600165292419`"}]], "Output",
 CellChangeTimes->{
  3.5673462861091833`*^9, 3.567510456010126*^9, 3.567510496398597*^9, 
   3.567510661384486*^9, 3.5675108565252295`*^9, 3.5675109484249907`*^9, 
   3.5675110392639503`*^9, 3.567511182472202*^9, 3.567511314448434*^9, 
   3.567511462274293*^9, 3.5675114995115585`*^9, 3.567511714323936*^9, 
   3.5675119770439973`*^9, {3.5675120719701643`*^9, 3.5675120852769876`*^9}, 
   3.5675122599504943`*^9, 3.5675138868305516`*^9, 3.567514066137267*^9, 
   3.5675161989146147`*^9, 3.56751625920872*^9, 3.567516295650384*^9, 
   3.567516359875697*^9, {3.567516406628979*^9, 3.567516477312703*^9}, 
   3.5675165355944057`*^9, 3.5675165969337134`*^9, 3.56751662916337*^9, 
   3.567516661533427*^9, 3.5675167055723042`*^9, {3.5675167396115637`*^9, 
   3.567516766131611*^9}, {3.567516802838475*^9, 3.567516829093321*^9}, {
   3.567516873397399*^9, 3.567516917171076*^9}, 3.56753241238629*^9, {
   3.5675324502631564`*^9, 3.567532480059209*^9}, 3.568099719811407*^9, 
   3.568101007640469*^9}],

Cell[BoxData["136.4558809314361`"], "Output",
 CellChangeTimes->{
  3.5673462861091833`*^9, 3.567510456010126*^9, 3.567510496398597*^9, 
   3.567510661384486*^9, 3.5675108565252295`*^9, 3.5675109484249907`*^9, 
   3.5675110392639503`*^9, 3.567511182472202*^9, 3.567511314448434*^9, 
   3.567511462274293*^9, 3.5675114995115585`*^9, 3.567511714323936*^9, 
   3.5675119770439973`*^9, {3.5675120719701643`*^9, 3.5675120852769876`*^9}, 
   3.5675122599504943`*^9, 3.5675138868305516`*^9, 3.567514066137267*^9, 
   3.5675161989146147`*^9, 3.56751625920872*^9, 3.567516295650384*^9, 
   3.567516359875697*^9, {3.567516406628979*^9, 3.567516477312703*^9}, 
   3.5675165355944057`*^9, 3.5675165969337134`*^9, 3.56751662916337*^9, 
   3.567516661533427*^9, 3.5675167055723042`*^9, {3.5675167396115637`*^9, 
   3.567516766131611*^9}, {3.567516802838475*^9, 3.567516829093321*^9}, {
   3.567516873397399*^9, 3.567516917171076*^9}, 3.56753241238629*^9, {
   3.5675324502631564`*^9, 3.567532480059209*^9}, 3.568099719811407*^9, 
   3.568101007640469*^9}],

Cell[BoxData["72.29652865790015`"], "Output",
 CellChangeTimes->{
  3.5673462861091833`*^9, 3.567510456010126*^9, 3.567510496398597*^9, 
   3.567510661384486*^9, 3.5675108565252295`*^9, 3.5675109484249907`*^9, 
   3.5675110392639503`*^9, 3.567511182472202*^9, 3.567511314448434*^9, 
   3.567511462274293*^9, 3.5675114995115585`*^9, 3.567511714323936*^9, 
   3.5675119770439973`*^9, {3.5675120719701643`*^9, 3.5675120852769876`*^9}, 
   3.5675122599504943`*^9, 3.5675138868305516`*^9, 3.567514066137267*^9, 
   3.5675161989146147`*^9, 3.56751625920872*^9, 3.567516295650384*^9, 
   3.567516359875697*^9, {3.567516406628979*^9, 3.567516477312703*^9}, 
   3.5675165355944057`*^9, 3.5675165969337134`*^9, 3.56751662916337*^9, 
   3.567516661533427*^9, 3.5675167055723042`*^9, {3.5675167396115637`*^9, 
   3.567516766131611*^9}, {3.567516802838475*^9, 3.567516829093321*^9}, {
   3.567516873397399*^9, 3.567516917171076*^9}, 3.56753241238629*^9, {
   3.5675324502631564`*^9, 3.567532480059209*^9}, 3.568099719811407*^9, 
   3.568101007640469*^9}],

Cell[BoxData["1.999999999952452`"], "Output",
 CellChangeTimes->{
  3.5673462861091833`*^9, 3.567510456010126*^9, 3.567510496398597*^9, 
   3.567510661384486*^9, 3.5675108565252295`*^9, 3.5675109484249907`*^9, 
   3.5675110392639503`*^9, 3.567511182472202*^9, 3.567511314448434*^9, 
   3.567511462274293*^9, 3.5675114995115585`*^9, 3.567511714323936*^9, 
   3.5675119770439973`*^9, {3.5675120719701643`*^9, 3.5675120852769876`*^9}, 
   3.5675122599504943`*^9, 3.5675138868305516`*^9, 3.567514066137267*^9, 
   3.5675161989146147`*^9, 3.56751625920872*^9, 3.567516295650384*^9, 
   3.567516359875697*^9, {3.567516406628979*^9, 3.567516477312703*^9}, 
   3.5675165355944057`*^9, 3.5675165969337134`*^9, 3.56751662916337*^9, 
   3.567516661533427*^9, 3.5675167055723042`*^9, {3.5675167396115637`*^9, 
   3.567516766131611*^9}, {3.567516802838475*^9, 3.567516829093321*^9}, {
   3.567516873397399*^9, 3.567516917171076*^9}, 3.56753241238629*^9, {
   3.5675324502631564`*^9, 3.567532480059209*^9}, 3.568099719811407*^9, 
   3.568101007656069*^9}],

Cell[BoxData["129.74735992771798`"], "Output",
 CellChangeTimes->{
  3.5673462861091833`*^9, 3.567510456010126*^9, 3.567510496398597*^9, 
   3.567510661384486*^9, 3.5675108565252295`*^9, 3.5675109484249907`*^9, 
   3.5675110392639503`*^9, 3.567511182472202*^9, 3.567511314448434*^9, 
   3.567511462274293*^9, 3.5675114995115585`*^9, 3.567511714323936*^9, 
   3.5675119770439973`*^9, {3.5675120719701643`*^9, 3.5675120852769876`*^9}, 
   3.5675122599504943`*^9, 3.5675138868305516`*^9, 3.567514066137267*^9, 
   3.5675161989146147`*^9, 3.56751625920872*^9, 3.567516295650384*^9, 
   3.567516359875697*^9, {3.567516406628979*^9, 3.567516477312703*^9}, 
   3.5675165355944057`*^9, 3.5675165969337134`*^9, 3.56751662916337*^9, 
   3.567516661533427*^9, 3.5675167055723042`*^9, {3.5675167396115637`*^9, 
   3.567516766131611*^9}, {3.567516802838475*^9, 3.567516829093321*^9}, {
   3.567516873397399*^9, 3.567516917171076*^9}, 3.56753241238629*^9, {
   3.5675324502631564`*^9, 3.567532480059209*^9}, 3.568099719811407*^9, 
   3.568101007656069*^9}],

Cell[BoxData["69.53339663162086`"], "Output",
 CellChangeTimes->{
  3.5673462861091833`*^9, 3.567510456010126*^9, 3.567510496398597*^9, 
   3.567510661384486*^9, 3.5675108565252295`*^9, 3.5675109484249907`*^9, 
   3.5675110392639503`*^9, 3.567511182472202*^9, 3.567511314448434*^9, 
   3.567511462274293*^9, 3.5675114995115585`*^9, 3.567511714323936*^9, 
   3.5675119770439973`*^9, {3.5675120719701643`*^9, 3.5675120852769876`*^9}, 
   3.5675122599504943`*^9, 3.5675138868305516`*^9, 3.567514066137267*^9, 
   3.5675161989146147`*^9, 3.56751625920872*^9, 3.567516295650384*^9, 
   3.567516359875697*^9, {3.567516406628979*^9, 3.567516477312703*^9}, 
   3.5675165355944057`*^9, 3.5675165969337134`*^9, 3.56751662916337*^9, 
   3.567516661533427*^9, 3.5675167055723042`*^9, {3.5675167396115637`*^9, 
   3.567516766131611*^9}, {3.567516802838475*^9, 3.567516829093321*^9}, {
   3.567516873397399*^9, 3.567516917171076*^9}, 3.56753241238629*^9, {
   3.5675324502631564`*^9, 3.567532480059209*^9}, 3.568099719811407*^9, 
   3.568101007656069*^9}],

Cell[BoxData["48.61324762359017`"], "Output",
 CellChangeTimes->{
  3.5673462861091833`*^9, 3.567510456010126*^9, 3.567510496398597*^9, 
   3.567510661384486*^9, 3.5675108565252295`*^9, 3.5675109484249907`*^9, 
   3.5675110392639503`*^9, 3.567511182472202*^9, 3.567511314448434*^9, 
   3.567511462274293*^9, 3.5675114995115585`*^9, 3.567511714323936*^9, 
   3.5675119770439973`*^9, {3.5675120719701643`*^9, 3.5675120852769876`*^9}, 
   3.5675122599504943`*^9, 3.5675138868305516`*^9, 3.567514066137267*^9, 
   3.5675161989146147`*^9, 3.56751625920872*^9, 3.567516295650384*^9, 
   3.567516359875697*^9, {3.567516406628979*^9, 3.567516477312703*^9}, 
   3.5675165355944057`*^9, 3.5675165969337134`*^9, 3.56751662916337*^9, 
   3.567516661533427*^9, 3.5675167055723042`*^9, {3.5675167396115637`*^9, 
   3.567516766131611*^9}, {3.567516802838475*^9, 3.567516829093321*^9}, {
   3.567516873397399*^9, 3.567516917171076*^9}, 3.56753241238629*^9, {
   3.5675324502631564`*^9, 3.567532480059209*^9}, 3.568099719811407*^9, 
   3.568101007656069*^9}],

Cell[BoxData["1.1870316110293102`"], "Output",
 CellChangeTimes->{
  3.5673462861091833`*^9, 3.567510456010126*^9, 3.567510496398597*^9, 
   3.567510661384486*^9, 3.5675108565252295`*^9, 3.5675109484249907`*^9, 
   3.5675110392639503`*^9, 3.567511182472202*^9, 3.567511314448434*^9, 
   3.567511462274293*^9, 3.5675114995115585`*^9, 3.567511714323936*^9, 
   3.5675119770439973`*^9, {3.5675120719701643`*^9, 3.5675120852769876`*^9}, 
   3.5675122599504943`*^9, 3.5675138868305516`*^9, 3.567514066137267*^9, 
   3.5675161989146147`*^9, 3.56751625920872*^9, 3.567516295650384*^9, 
   3.567516359875697*^9, {3.567516406628979*^9, 3.567516477312703*^9}, 
   3.5675165355944057`*^9, 3.5675165969337134`*^9, 3.56751662916337*^9, 
   3.567516661533427*^9, 3.5675167055723042`*^9, {3.5675167396115637`*^9, 
   3.567516766131611*^9}, {3.567516802838475*^9, 3.567516829093321*^9}, {
   3.567516873397399*^9, 3.567516917171076*^9}, 3.56753241238629*^9, {
   3.5675324502631564`*^9, 3.567532480059209*^9}, 3.568099719811407*^9, 
   3.568101007656069*^9}],

Cell[BoxData["1.3219459281365367`"], "Output",
 CellChangeTimes->{
  3.5673462861091833`*^9, 3.567510456010126*^9, 3.567510496398597*^9, 
   3.567510661384486*^9, 3.5675108565252295`*^9, 3.5675109484249907`*^9, 
   3.5675110392639503`*^9, 3.567511182472202*^9, 3.567511314448434*^9, 
   3.567511462274293*^9, 3.5675114995115585`*^9, 3.567511714323936*^9, 
   3.5675119770439973`*^9, {3.5675120719701643`*^9, 3.5675120852769876`*^9}, 
   3.5675122599504943`*^9, 3.5675138868305516`*^9, 3.567514066137267*^9, 
   3.5675161989146147`*^9, 3.56751625920872*^9, 3.567516295650384*^9, 
   3.567516359875697*^9, {3.567516406628979*^9, 3.567516477312703*^9}, 
   3.5675165355944057`*^9, 3.5675165969337134`*^9, 3.56751662916337*^9, 
   3.567516661533427*^9, 3.5675167055723042`*^9, {3.5675167396115637`*^9, 
   3.567516766131611*^9}, {3.567516802838475*^9, 3.567516829093321*^9}, {
   3.567516873397399*^9, 3.567516917171076*^9}, 3.56753241238629*^9, {
   3.5675324502631564`*^9, 3.567532480059209*^9}, 3.568099719811407*^9, 
   3.568101007671669*^9}],

Cell[BoxData["1.0877868327580067`"], "Output",
 CellChangeTimes->{
  3.5673462861091833`*^9, 3.567510456010126*^9, 3.567510496398597*^9, 
   3.567510661384486*^9, 3.5675108565252295`*^9, 3.5675109484249907`*^9, 
   3.5675110392639503`*^9, 3.567511182472202*^9, 3.567511314448434*^9, 
   3.567511462274293*^9, 3.5675114995115585`*^9, 3.567511714323936*^9, 
   3.5675119770439973`*^9, {3.5675120719701643`*^9, 3.5675120852769876`*^9}, 
   3.5675122599504943`*^9, 3.5675138868305516`*^9, 3.567514066137267*^9, 
   3.5675161989146147`*^9, 3.56751625920872*^9, 3.567516295650384*^9, 
   3.567516359875697*^9, {3.567516406628979*^9, 3.567516477312703*^9}, 
   3.5675165355944057`*^9, 3.5675165969337134`*^9, 3.56751662916337*^9, 
   3.567516661533427*^9, 3.5675167055723042`*^9, {3.5675167396115637`*^9, 
   3.567516766131611*^9}, {3.567516802838475*^9, 3.567516829093321*^9}, {
   3.567516873397399*^9, 3.567516917171076*^9}, 3.56753241238629*^9, {
   3.5675324502631564`*^9, 3.567532480059209*^9}, 3.568099719811407*^9, 
   3.568101007671669*^9}],

Cell[BoxData["0.9493098349235608`"], "Output",
 CellChangeTimes->{
  3.5673462861091833`*^9, 3.567510456010126*^9, 3.567510496398597*^9, 
   3.567510661384486*^9, 3.5675108565252295`*^9, 3.5675109484249907`*^9, 
   3.5675110392639503`*^9, 3.567511182472202*^9, 3.567511314448434*^9, 
   3.567511462274293*^9, 3.5675114995115585`*^9, 3.567511714323936*^9, 
   3.5675119770439973`*^9, {3.5675120719701643`*^9, 3.5675120852769876`*^9}, 
   3.5675122599504943`*^9, 3.5675138868305516`*^9, 3.567514066137267*^9, 
   3.5675161989146147`*^9, 3.56751625920872*^9, 3.567516295650384*^9, 
   3.567516359875697*^9, {3.567516406628979*^9, 3.567516477312703*^9}, 
   3.5675165355944057`*^9, 3.5675165969337134`*^9, 3.56751662916337*^9, 
   3.567516661533427*^9, 3.5675167055723042`*^9, {3.5675167396115637`*^9, 
   3.567516766131611*^9}, {3.567516802838475*^9, 3.567516829093321*^9}, {
   3.567516873397399*^9, 3.567516917171076*^9}, 3.56753241238629*^9, {
   3.5675324502631564`*^9, 3.567532480059209*^9}, 3.568099719811407*^9, 
   3.568101007671669*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
plot1=Graphics[{Circle[{xc[[1]],yc[[1]]},Cr],      \
Circle[{xc[[2]],yc[[2]]},Nr],Circle[{xc[[3]],yc[[3]]},Oa],Circle[{xc[[4]],yc[[\
4]]},NH],
Circle[{xc[[5]],yc[[5]]},NH],Circle[{xc[[6]],yc[[6]]},CH],Circle[{xc[[7]],yc[[\
7]]},CN],Circle[{xc[[8]],yc[[8]]},CO],Circle[{xc[[9]],yc[[9]]},OL],Circle[{xc[\
[10]],yc[[10]]},NM],Circle[{xc[[11]],yc[[11]]},NM],Circle[{xc[[12]],yc[[12]]},\
OC],Circle[{xc[[13]],yc[[13]]},OC],Circle[{xc[[14]],yc[[14]]},OO],Circle[{xc[[\
15]],yc[[15]]},OO],Circle[{xc[[16]],yc[[16]]},OM],Circle[{xc[[17]],yc[[17]]},\
OM],Disk[{xn[[1]],yn[[1]]},0.08], Disk[{xn[[2]],yn[[2]]},0.08], \
Disk[{xn[[3]],yn[[3]]},0.08],Disk[{xn[[4]],yn[[4]]},0.08], \
(*Disk[{xn[[5]],yn[[5]]},0.08],*) Disk[{xn[[6]],yn[[6]]},0.08],
{Thickness[0.008],
Line[{{xn[[2]],yn[[2]]},{xn[[1]],yn[[1]]}}],
Line[{{xn[[3]]+0.02,yn[[3]]-0.07},{xn[[1]]+0.03,yn[[1]]-0.07}}],
Line[{{xn[[3]]-0.02,yn[[3]]+0.06},{xn[[1]]-0.02,yn[[1]]+0.06}}],
Line[{{xn[[2]],yn[[2]]},{xn[[4]],yn[[4]]}}],
(*Line[{{xn[[2]],yn[[2]]},{xn[[5]],yn[[5]]}}],*)
Line[{{xn[[6]],yn[[6]]},{xn[[1]],yn[[1]]}}]}} ] /. 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.5675116826870804`*^9, 
   3.5675117055879207`*^9}}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{0, 0}, 0.25559559017708283`], 
   CircleBox[{-2.4981167513639795`, 0}, 0.18966729854432163`], 
   CircleBox[{0.7843483199037025, 2.1015683080714793`}, 0.18771946482544016`],
    CircleBox[{-3.0195425457907588`, 1.1117950366024378`}, 
    1.0383282480711433`], 
   CircleBox[{-3.0195425457907588`, -1.1117950366024378`}, 
    1.0383282480711433`], 
   CircleBox[{0.9869942279020033, -1.1200471330977853`}, 1.237275196687085], 
   CircleBox[{-1.2820225214983703`, 0}, 1.0264269313212875`], 
   CircleBox[{0.40404098920022924`, 1.0825799157309375`}, 0.8999252425636789],
    CircleBox[{1.1701963736967702`, 3.135402410861584}, 0.9157712114866523], 
   CircleBox[{-2.4981167513639795`, 0}, 0.8768151106038778], 
   CircleBox[{-2.4981167513639795`, 0}, 0.8768151106038778], 
   CircleBox[{0, 0}, 1.3099009653202152`], 
   CircleBox[{0, 0}, 1.3099009653202152`], 
   CircleBox[{0.7843483199037025, 2.1015683080714793`}, 1.1059447491867593`], 
   CircleBox[{0.7843483199037025, 2.1015683080714793`}, 1.1059447491867593`], 
   CircleBox[{-0.19463616399049166`, 2.4669453611061445`}, 
    0.8572259933635882], 
   CircleBox[{1.7633328037978966`, 1.7361912550368142`}, 0.8572259933635882], 
   DiskBox[{0, 0}, 0.08], DiskBox[{-2.4981167513639795`, 0}, 0.08], 
   DiskBox[{0.7843483199037025, 2.1015683080714793`}, 0.08], 
   DiskBox[{-3.2598496976882316`, 1.6241829959157847`}, 0.08], 
   DiskBox[{1.3590493615167112`, -1.5422575918613213`}, 0.08], 
   {Thickness[0.008], LineBox[{{-2.4981167513639795`, 0}, {0, 0}}], 
    LineBox[{{0.8043483199037025, 2.0315683080714795`}, {0.03, -0.07}}], 
    LineBox[{{0.7643483199037024, 2.1615683080714794`}, {-0.02, 0.06}}], 
    LineBox[{{-2.4981167513639795`, 0}, {-3.2598496976882316`, 
      1.6241829959157847`}}], 
    LineBox[{{1.3590493615167112`, -1.5422575918613213`}, {0, 0}}]}},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->Automatic]], "Output",
 CellChangeTimes->{
  3.5673462861403837`*^9, 3.567510456041326*^9, 3.5675104964297967`*^9, 
   3.5675106614156866`*^9, 3.567510856556429*^9, 3.567510948456191*^9, 
   3.5675110392951508`*^9, 3.5675111825034018`*^9, 3.5675113144952335`*^9, 
   3.567511462321093*^9, 3.567511499542759*^9, 3.567511714355136*^9, 
   3.567511977075197*^9, {3.567512072001364*^9, 3.5675120853081875`*^9}, 
   3.567512259981694*^9, 3.567513886846152*^9, 3.5675140661684666`*^9, 
   3.5675161989614143`*^9, 3.5675162592399206`*^9, 3.5675162956815844`*^9, 
   3.5675163599068975`*^9, {3.567516406644579*^9, 3.5675164773439035`*^9}, 
   3.5675165356256056`*^9, 3.5675165969649134`*^9, 3.5675166291789703`*^9, 
   3.5675166615646267`*^9, 3.567516705603504*^9, {3.567516739642764*^9, 
   3.5675167661628103`*^9}, {3.567516802869675*^9, 3.5675168291245213`*^9}, {
   3.567516873444199*^9, 3.567516917202276*^9}, 3.56753241243309*^9, {
   3.567532450278756*^9, 3.5675324801060085`*^9}, 3.568099719858207*^9, 
   3.568101007687269*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
plot2=Graphics[{Circle[{xc[[1]],zc[[1]]},Cr],      \
Circle[{xc[[2]],zc[[2]]},Nr],Circle[{xc[[3]],zc[[3]]},Oa],Circle[{xc[[4]],zc[[\
4]]},NH],
Circle[{xc[[5]],zc[[5]]},NH],Circle[{xc[[6]],zc[[6]]},CH],Circle[{xc[[7]],zc[[\
7]]},CN],Circle[{xc[[8]],zc[[8]]},CO],Circle[{xc[[9]],zc[[9]]},OL],Circle[{xc[\
[10]],zc[[10]]},NM],Circle[{xc[[11]],zc[[11]]},NM],Circle[{xc[[12]],zc[[12]]},\
OC],Circle[{xc[[13]],zc[[13]]},OC],Circle[{xc[[14]],zc[[14]]},OO],Circle[{xc[[\
15]],zc[[15]]},OO],Circle[{xc[[16]],zc[[16]]},OM],Circle[{xc[[17]],zc[[17]]},\
OM],Disk[{xn[[1]],zn[[1]]},0.08], Disk[{xn[[2]],zn[[2]]},0.08], \
Disk[{xn[[3]],zn[[3]]},0.08],Disk[{xn[[4]],zn[[4]]},0.08], \
(*Disk[{xn[[5]],zn[[5]]},0.08],*) Disk[{xn[[6]],zn[[6]]},0.08]} ] /. t[[2]];

(* xz-plane *)
Show[plot2,{AspectRatio \[Rule] Automatic,Axes -> True,GridLines -> \
Automatic, PlotRange \[Rule] Automatic, Frame -> True}]\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668391703953342`*^9, 3.5668391881949654`*^9}, {
   3.5668412398393693`*^9, 3.5668412401201696`*^9}, {3.566841311459095*^9, 
   3.5668413117242956`*^9}, {3.5668419284557786`*^9, 
   3.5668419511070185`*^9}, {3.566842124953724*^9, 3.566842139087349*^9}, 
   3.5670689313531685`*^9, {3.5670705170023537`*^9, 3.567070548654809*^9}, {
   3.567076321179748*^9, 3.567076324596154*^9}, {3.5670763612250185`*^9, 
   3.567076444107964*^9}, {3.5670764854792366`*^9, 3.5670764972884574`*^9}, {
   3.567076691212398*^9, 3.5670767228336535`*^9}, {3.5670850083498063`*^9, 
   3.5670850144962173`*^9}, {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.56751178485166*^9, 3.5675117914972715`*^9}}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{0, 0}, 0.25559559017708283`], 
   CircleBox[{-2.4981167513639795`, 0}, 0.18966729854432163`], 
   CircleBox[{0.7843483199037025, 0}, 0.18771946482544016`], 
   CircleBox[{-3.0195425457907588`, 0}, 1.0383282480711433`], 
   CircleBox[{-3.0195425457907588`, 0}, 1.0383282480711433`], 
   CircleBox[{0.9869942279020033, 0}, 1.237275196687085], 
   CircleBox[{-1.2820225214983703`, 0}, 1.0264269313212875`], 
   CircleBox[{0.40404098920022924`, 0}, 0.8999252425636789], 
   CircleBox[{1.1701963736967702`, 0}, 0.9157712114866523], 
   CircleBox[{-2.4981167513639795`, 0.8768151106038778}, 0.8768151106038778], 
   CircleBox[{-2.4981167513639795`, -0.8768151106038778}, 0.8768151106038778],
    CircleBox[{0, 1.3099009653202152`}, 1.3099009653202152`], 
   CircleBox[{0, -1.3099009653202152`}, 1.3099009653202152`], 
   CircleBox[{0.7843483199037025, 1.1059447491867593`}, 1.1059447491867593`], 
   CircleBox[{0.7843483199037025, -1.1059447491867593`}, 1.1059447491867593`],
    CircleBox[{-0.19463616399049166`, 0}, 0.8572259933635882], 
   CircleBox[{1.7633328037978966`, 0}, 0.8572259933635882], 
   DiskBox[{0, 0}, 0.08], DiskBox[{-2.4981167513639795`, 0}, 0.08], 
   DiskBox[{0.7843483199037025, 0}, 0.08], 
   DiskBox[{-3.2598496976882316`, 0}, 0.08], 
   DiskBox[{1.3590493615167112`, 0}, 0.08]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->Automatic]], "Output",
 CellChangeTimes->{
  3.567346286155984*^9, 3.567510456056926*^9, 3.567510496445397*^9, 
   3.567510661431287*^9, 3.567510856556429*^9, 3.567510948456191*^9, 
   3.5675110393107505`*^9, 3.567511182534602*^9, 3.5675113145108337`*^9, 
   3.5675114623366933`*^9, 3.567511499558359*^9, 3.567511714370736*^9, 
   3.5675119770907974`*^9, {3.567512072016964*^9, 3.5675120853237877`*^9}, 
   3.5675122599972944`*^9, 3.5675138868773518`*^9, 3.5675140661684666`*^9, 
   3.5675161989614143`*^9, 3.5675162592555203`*^9, 3.5675162956971846`*^9, 
   3.5675163599224973`*^9, {3.567516406660179*^9, 3.5675164773595037`*^9}, 
   3.5675165356256056`*^9, 3.5675165969805136`*^9, 3.56751662919457*^9, 
   3.567516661580227*^9, 3.5675167056191044`*^9, {3.567516739642764*^9, 
   3.5675167661784105`*^9}, {3.567516802885275*^9, 3.5675168291401215`*^9}, {
   3.5675168734597993`*^9, 3.567516917217876*^9}, 3.56753241243309*^9, {
   3.5675324502943563`*^9, 3.5675324801060085`*^9}, 3.568099719873807*^9, 
   3.5681010077028694`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
plot3=Graphics3D[{
Sphere[{xc[[1]],yc[[1]],zc[[1]]},Cr],      
Sphere[{xc[[2]],yc[[2]],zc[[2]]},Nr],
Sphere[{xc[[3]],yc[[3]],zc[[3]]},Oa],
{Opacity[0.6],
Sphere[{xc[[4]],yc[[4]],zc[[4]]},NH],
Sphere[{xc[[5]],yc[[5]],zc[[5]]},NH],
Sphere[{xc[[6]],yc[[6]],zc[[6]]},CH],
Sphere[{xc[[7]],yc[[7]],zc[[7]]},CN],
Sphere[{xc[[8]],yc[[8]],zc[[8]]},CO],
{Lighter[Blue],
Sphere[{xc[[9]],yc[[9]],zc[[9]]},OL]},
{Lighter[Blue,0.6],
Sphere[{xc[[16]],yc[[16]],zc[[16]]},OM],
Sphere[{xc[[17]],yc[[17]],zc[[17]]},OM]}},
{Opacity[0.4],{Lighter[Green],
Sphere[{xc[[10]],yc[[10]],zc[[10]]},NM],
Sphere[{xc[[11]],yc[[11]],zc[[11]]},NM],
Sphere[{xc[[12]],yc[[12]],zc[[12]]},OC],
Sphere[{xc[[13]],yc[[13]],zc[[13]]},OC],
Sphere[{xc[[14]],yc[[14]],zc[[14]]},OO],
Sphere[{xc[[15]],yc[[15]],zc[[15]]},OO]}},
{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]],zn[[2]]},{xn[[1]],yn[[1]],zn[[1]]}},0.05],
Cylinder[{{xn[[3]],yn[[3]]-0.14,zn[[3]]},{xn[[1]],yn[[1]]-0.14,zn[[1]]}},0.05],
Cylinder[{{xn[[3]],yn[[3]]+0.14,zn[[3]]},{xn[[1]],yn[[1]]+0.14,zn[[1]]}},0.05],
Cylinder[{{xn[[2]],yn[[2]],zn[[2]]},{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[[1]],yn[[1]],zn[[1]]}},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.5675118112001057`*^9, 3.56751183603535*^9}}],

Cell[BoxData[
 Graphics3DBox[{SphereBox[{0, 0, 0}, 0.25559559017708283`], 
   SphereBox[{-2.4981167513639795`, 0, 0}, 0.18966729854432163`], 
   SphereBox[{0.7843483199037025, 2.1015683080714793`, 0}, 
    0.18771946482544016`], 
   {Opacity[0.6], 
    SphereBox[{-3.0195425457907588`, 1.1117950366024378`, 0}, 
     1.0383282480711433`], 
    SphereBox[{-3.0195425457907588`, -1.1117950366024378`, 0}, 
     1.0383282480711433`], 
    SphereBox[{0.9869942279020033, -1.1200471330977853`, 0}, 
     1.237275196687085], 
    SphereBox[{-1.2820225214983703`, 0, 0}, 1.0264269313212875`], 
    SphereBox[{0.40404098920022924`, 1.0825799157309375`, 0}, 
     0.8999252425636789], 
    {RGBColor[
      NCache[
       Rational[1, 3], 0.3333333333333333], 
      NCache[
       Rational[1, 3], 0.3333333333333333], 1], 
     SphereBox[{1.1701963736967702`, 3.135402410861584, 0}, 
      0.9157712114866523]}, 
    {RGBColor[0.6, 0.6, 1.], 
     SphereBox[{-0.19463616399049166`, 2.4669453611061445`, 0}, 
      0.8572259933635882], 
     SphereBox[{1.7633328037978966`, 1.7361912550368142`, 0}, 
      0.8572259933635882]}}, 
   {Opacity[0.4], 
    {RGBColor[
      NCache[
       Rational[1, 3], 0.3333333333333333], 1, 
      NCache[
       Rational[1, 3], 0.3333333333333333]], 
     SphereBox[{-2.4981167513639795`, 0, 0.8768151106038778}, 
      0.8768151106038778], 
     SphereBox[{-2.4981167513639795`, 0, -0.8768151106038778}, 
      0.8768151106038778], 
     SphereBox[{0, 0, 1.3099009653202152`}, 1.3099009653202152`], 
     SphereBox[{0, 0, -1.3099009653202152`}, 1.3099009653202152`], 
     SphereBox[{0.7843483199037025, 2.1015683080714793`, 1.1059447491867593`},
       1.1059447491867593`], 
     SphereBox[{0.7843483199037025, 
      2.1015683080714793`, -1.1059447491867593`}, 1.1059447491867593`]}}, 
   {RGBColor[0, 0, 0], SphereBox[{0, 0, 0}, 0.08], 
    SphereBox[{-2.4981167513639795`, 0, 0}, 0.08], 
    SphereBox[{0.7843483199037025, 2.1015683080714793`, 0}, 0.08], 
    SphereBox[{-3.2598496976882316`, 1.6241829959157847`, 0}, 0.08], 
    SphereBox[{1.3590493615167112`, -1.5422575918613213`, 0}, 0.08], 
    CylinderBox[{{-2.4981167513639795`, 0, 0}, {0, 0, 0}}, 0.05], 
    CylinderBox[{{0.7843483199037025, 1.9615683080714792`, 0}, {0, -0.14, 0}},
      0.05], CylinderBox[{{0.7843483199037025, 2.2415683080714794`, 0}, {
      0, 0.14, 0}}, 0.05], 
    CylinderBox[{{-2.4981167513639795`, 0, 0}, {-3.2598496976882316`, 
       1.6241829959157847`, 0}}, 0.05], 
    CylinderBox[{{1.3590493615167112`, -1.5422575918613213`, 0}, {0, 0, 0}}, 
     0.05]}},
  AspectRatio->Automatic,
  AutomaticImageSize->True,
  Boxed->False,
  ImageSize->{436.39571126404326`, 430.1736774502665},
  ViewPoint->{-0.023648508190708933`, -0.4672200995998511, 3.351290218198101},
  ViewVertical->{-0.4113999549274383, 0.649182900240934, 
   0.7577649280257959}]], "Output",
 CellChangeTimes->{
  3.5673462864055843`*^9, 3.5675104563221264`*^9, 3.5675104964609966`*^9, 
   3.5675106614624867`*^9, 3.567510856572029*^9, 3.567510948471791*^9, 
   3.56751103932635*^9, 3.567511182550202*^9, 3.567511314526434*^9, 
   3.5675114623522935`*^9, 3.567511499573959*^9, 3.567511714401936*^9, 
   3.5675119771063976`*^9, {3.567512072032564*^9, 3.5675120853237877`*^9}, 
   3.5675122600128946`*^9, 3.567513886892952*^9, 3.567514066184067*^9, 
   3.5675161992110147`*^9, 3.56751625927112*^9, 3.5675162957127843`*^9, 
   3.567516359938097*^9, {3.5675164066757793`*^9, 3.5675164773751035`*^9}, 
   3.567516535641206*^9, 3.567516596996114*^9, 3.56751662921017*^9, 
   3.567516661595827*^9, 3.5675167056347046`*^9, {3.5675167396583643`*^9, 
   3.5675167661784105`*^9}, {3.567516802900875*^9, 3.567516829155721*^9}, {
   3.5675168734597993`*^9, 3.5675169172334757`*^9}, 3.56753241244869*^9, {
   3.5675324503099566`*^9, 3.5675324801216087`*^9}, 3.5680997202014074`*^9, 
   3.568101007718469*^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.`", "0.`", "0.`"},
     {"O", 
      RowBox[{"-", "1.3219459281365367`"}], "0.`", "0.`"},
     {"O", "0.41505909088168164`", "1.1121016125603413`", "0.`"},
     {"H", 
      RowBox[{"-", "1.7250374834735653`"}], "0.8594802852297272`", "0.`"},
     {"L", 
      RowBox[{"-", "1.7250374834735653`"}], 
      RowBox[{"-", "0.8594802852297272`"}], "0.`"},
     {"H", "0.7191776639793287`", 
      RowBox[{"-", "0.8161272456883985`"}], "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.5673462864211845`*^9, 3.5675104563377266`*^9, 3.5675104964609966`*^9, 
   3.567510661524887*^9, 3.5675108565876293`*^9, 3.5675109484873905`*^9, 
   3.56751103938875*^9, 3.567511182597002*^9, 3.567511314588834*^9, 
   3.567511462414694*^9, 3.567511499573959*^9, 3.5675117144643364`*^9, 
   3.567511977121998*^9, {3.5675120720481644`*^9, 3.567512085339388*^9}, 
   3.5675122600752945`*^9, 3.567513886955352*^9, 3.567514066199667*^9, 
   3.567516199226615*^9, 3.5675162592867203`*^9, 3.5675162957127843`*^9, 
   3.567516359938097*^9, {3.5675164066913795`*^9, 3.5675164773751035`*^9}, 
   3.567516535656806*^9, 3.567516596996114*^9, 3.56751662922577*^9, 
   3.567516661595827*^9, 3.5675167056347046`*^9, {3.567516739673964*^9, 
   3.5675167661940107`*^9}, {3.567516802900875*^9, 3.567516829155721*^9}, {
   3.567516873475399*^9, 3.5675169172334757`*^9}, 3.5675324124642897`*^9, {
   3.5675324503255568`*^9, 3.567532480137209*^9}, 3.5680997202014074`*^9, 
   3.5681010077808695`*^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.3219459281365367`", "1.1870316110293102`", 
      "1.9272936154326317`", "1.9272936154326317`", "1.087786832758007`"},
     {"1.3219459281365367`", "0.`", "2.0625121654801926`", 
      "0.9493098349235607`", "0.9493098349235607`", "2.1982377486174935`"},
     {"1.1870316110293102`", "2.0625121654801926`", "0.`", 
      "2.154954960686962`", "2.909836512118946`", "1.9520641988126795`"},
     {"1.9272936154326317`", "0.9493098349235607`", "2.154954960686962`", 
      "0.`", "1.7189605704594544`", "2.9634183445318194`"},
     {"1.9272936154326317`", "0.9493098349235607`", "2.909836512118946`", 
      "1.7189605704594544`", "0.`", "2.444599593609482`"},
     {"1.087786832758007`", "2.1982377486174935`", "1.9520641988126795`", 
      "2.9634183445318194`", "2.444599593609482`", "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.5673462864367847`*^9, 3.5675104563533263`*^9, 3.567510496476597*^9, 
   3.567510661524887*^9, 3.5675108565876293`*^9, 3.5675109484873905`*^9, 
   3.5675110394043503`*^9, 3.567511182612602*^9, 3.567511314620034*^9, 
   3.5675114624302936`*^9, 3.5675114995895586`*^9, 3.5675117144643364`*^9, 
   3.5675119771375976`*^9, {3.5675120720481644`*^9, 3.5675120853549876`*^9}, 
   3.5675122600908947`*^9, 3.567513886970952*^9, 3.5675140662152667`*^9, 
   3.567516199226615*^9, 3.5675162592867203`*^9, 3.5675162957283845`*^9, 
   3.567516359953697*^9, {3.5675164066913795`*^9, 3.567516477390703*^9}, 
   3.5675165356724057`*^9, 3.5675165970117135`*^9, 3.56751662922577*^9, 
   3.5675166616114273`*^9, 3.5675167056503043`*^9, {3.5675167396895638`*^9, 
   3.567516766209611*^9}, {3.5675168029164753`*^9, 3.567516829171321*^9}, {
   3.567516873490999*^9, 3.567516917249076*^9}, 3.56753241247989*^9, {
   3.5675324503255568`*^9, 3.567532480152809*^9}, 3.5680997202170076`*^9, 
   3.568101007812069*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  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", "]"}], "]"}]}], "}"}], ",", "Nr"}], "]"}], ",",
         "\n", 
        RowBox[{"Sphere", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "3", "]"}], "]"}], ",", 
            RowBox[{"yc", "[", 
             RowBox[{"[", "3", "]"}], "]"}], ",", 
            RowBox[{"zc", "[", 
             RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "Oa"}], "]"}], ",",
         "\n", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Lighter", "[", "Blue", "]"}], ",", "\n", 
          RowBox[{"Sphere", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xc", "[", 
               RowBox[{"[", "9", "]"}], "]"}], ",", 
              RowBox[{"yc", "[", 
               RowBox[{"[", "9", "]"}], "]"}], ",", 
              RowBox[{"zc", "[", 
               RowBox[{"[", "9", "]"}], "]"}]}], "}"}], ",", "OL"}], "]"}]}], 
         "}"}], ",", "\n", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Lighter", "[", 
           RowBox[{"Blue", ",", "0.6"}], "]"}], ",", "\n", 
          RowBox[{"Sphere", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xc", "[", 
               RowBox[{"[", "16", "]"}], "]"}], ",", 
              RowBox[{"yc", "[", 
               RowBox[{"[", "16", "]"}], "]"}], ",", 
              RowBox[{"zc", "[", 
               RowBox[{"[", "16", "]"}], "]"}]}], "}"}], ",", "OM"}], "]"}], 
          ",", "\n", 
          RowBox[{"Sphere", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xc", "[", 
               RowBox[{"[", "17", "]"}], "]"}], ",", 
              RowBox[{"yc", "[", 
               RowBox[{"[", "17", "]"}], "]"}], ",", 
              RowBox[{"zc", "[", 
               RowBox[{"[", "17", "]"}], "]"}]}], "}"}], ",", "OM"}], "]"}]}],
          "}"}], ",", "\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[{"(*", 
           RowBox[{
            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[{"[", "2", "]"}], "]"}], ",", 
                RowBox[{"yn", "[", 
                 RowBox[{"[", "2", "]"}], "]"}], ",", 
                RowBox[{"zn", "[", 
                 RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", 
              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[{"[", "3", "]"}], "]"}], ",", 
                RowBox[{
                 RowBox[{"yn", "[", 
                  RowBox[{"[", "3", "]"}], "]"}], "-", "0.14"}], ",", 
                RowBox[{"zn", "[", 
                 RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"xn", "[", 
                 RowBox[{"[", "1", "]"}], "]"}], ",", 
                RowBox[{
                 RowBox[{"yn", "[", 
                  RowBox[{"[", "1", "]"}], "]"}], "-", "0.14"}], ",", 
                RowBox[{"zn", "[", 
                 RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", 
            "0.05"}], "]"}], ",", "\n", 
          RowBox[{"Cylinder", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
                RowBox[{"xn", "[", 
                 RowBox[{"[", "3", "]"}], "]"}], ",", 
                RowBox[{
                 RowBox[{"yn", "[", 
                  RowBox[{"[", "3", "]"}], "]"}], "+", "0.14"}], ",", 
                RowBox[{"zn", "[", 
                 RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"xn", "[", 
                 RowBox[{"[", "1", "]"}], "]"}], ",", 
                RowBox[{
                 RowBox[{"yn", "[", 
                  RowBox[{"[", "1", "]"}], "]"}], "+", "0.14"}], ",", 
                RowBox[{"zn", "[", 
                 RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", 
            "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[{"[", "4", "]"}], "]"}], ",", 
                RowBox[{"yn", "[", 
                 RowBox[{"[", "4", "]"}], "]"}], ",", 
                RowBox[{"zn", "[", 
                 RowBox[{"[", "4", "]"}], "]"}]}], "}"}]}], "}"}], ",", 
            "0.05"}], "]"}], ",", "\n", 
          RowBox[{"(*", 
           RowBox[{
            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[{"[", "6", "]"}], "]"}], ",", 
                RowBox[{"yn", "[", 
                 RowBox[{"[", "6", "]"}], "]"}], ",", 
                RowBox[{"zn", "[", 
                 RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"xn", "[", 
                 RowBox[{"[", "1", "]"}], "]"}], ",", 
                RowBox[{"yn", "[", 
                 RowBox[{"[", "1", "]"}], "]"}], ",", 
                RowBox[{"zn", "[", 
                 RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", 
            "0.05"}], "]"}]}], "}"}]}], "}"}], "]"}], " ", "/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\n"}]], "Input",
 CellChangeTimes->{{3.567343131222042*^9, 3.5673432051973724`*^9}, {
   3.5673432772226987`*^9, 3.567343288626319*^9}, 3.5673433408552103`*^9, {
   3.5673433993709135`*^9, 3.56734341448734*^9}, {3.5673449337584085`*^9, 
   3.5673449778908863`*^9}, {3.5673450098397417`*^9, 
   3.5673450314145803`*^9}, {3.5673457004529552`*^9, 
   3.5673457120281754`*^9}, {3.5675118638189983`*^9, 3.567511888498242*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"(*", " ", 
    RowBox[{
     RowBox[{
      RowBox[{
      "Two", " ", "lone", " ", "pairs", " ", "of", " ", "O", " ", "atom"}], 
      ";", " ", 
      RowBox[{"one", " ", "of", " ", 
       SubscriptBox["\[Pi]", "x"]}]}], ",", " ", 
     RowBox[{
     "the", " ", "other", " ", "of", " ", "\[Sigma]", " ", "symmetry"}]}], 
    " ", "*)"}], "Text"], "\[IndentingNewLine]", 
  RowBox[{"Show", "[", 
   RowBox[{"plot4", ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", 
      RowBox[{"Boxed", "->", "False"}]}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.567343131222042*^9, 3.5673432051973724`*^9}, {
   3.5673432772226987`*^9, 3.567343288626319*^9}, 3.5673433408552103`*^9, {
   3.5673433993709135`*^9, 3.5673435539359846`*^9}}],

Cell[BoxData[
 Graphics3DBox[{SphereBox[{0, 0, 0}, 0.25559559017708283`], 
   SphereBox[{-2.4981167513639795`, 0, 0}, 0.18966729854432163`], 
   SphereBox[{0.7843483199037025, 2.1015683080714793`, 0}, 
    0.18771946482544016`], 
   {RGBColor[
     NCache[
      Rational[1, 3], 0.3333333333333333], 
     NCache[
      Rational[1, 3], 0.3333333333333333], 1], 
    SphereBox[{1.1701963736967702`, 3.135402410861584, 0}, 
     0.9157712114866523]}, 
   {RGBColor[0.6, 0.6, 1.], 
    SphereBox[{-0.19463616399049166`, 2.4669453611061445`, 0}, 
     0.8572259933635882], 
    SphereBox[{1.7633328037978966`, 1.7361912550368142`, 0}, 
     0.8572259933635882]}, 
   {RGBColor[0, 0, 0], SphereBox[{0, 0, 0}, 0.08], 
    SphereBox[{-2.4981167513639795`, 0, 0}, 0.08], 
    SphereBox[{0.7843483199037025, 2.1015683080714793`, 0}, 0.08], 
    SphereBox[{-3.2598496976882316`, 1.6241829959157847`, 0}, 0.08], 
    SphereBox[{1.3590493615167112`, -1.5422575918613213`, 0}, 0.08], 
    CylinderBox[{{-2.4981167513639795`, 0, 0}, {0, 0, 0}}, 0.05], 
    CylinderBox[{{0.7843483199037025, 1.9615683080714792`, 0}, {0, -0.14, 0}},
      0.05], CylinderBox[{{0.7843483199037025, 2.2415683080714794`, 0}, {
      0, 0.14, 0}}, 0.05], 
    CylinderBox[{{-2.4981167513639795`, 0, 0}, {-3.2598496976882316`, 
       1.6241829959157847`, 0}}, 0.05], 
    CylinderBox[{{1.3590493615167112`, -1.5422575918613213`, 0}, {0, 0, 0}}, 
     0.05]}},
  AspectRatio->Automatic,
  Boxed->False]], "Output",
 CellChangeTimes->{
  3.567346286467984*^9, 3.5675104563845263`*^9, 3.567510496492197*^9, 
   3.5675106615560865`*^9, 3.5675108566188297`*^9, 3.567510948518591*^9, 
   3.5675110394199505`*^9, 3.567511182628202*^9, 3.567511314635634*^9, 
   3.5675114624458933`*^9, 3.567511499620759*^9, 3.567511714495536*^9, 
   3.5675119771687975`*^9, {3.5675120720793643`*^9, 3.5675120853705873`*^9}, 
   3.5675122601064944`*^9, 3.5675138869865522`*^9, 3.567514066230867*^9, 
   3.567516199257815*^9, 3.5675162593179207`*^9, 3.5675162957595844`*^9, 
   3.5675163599692974`*^9, {3.567516406722579*^9, 3.5675164774063034`*^9}, 
   3.5675165356880054`*^9, 3.5675165970273137`*^9, 3.5675166292569704`*^9, 
   3.567516661627027*^9, 3.5675167056815042`*^9, {3.567516739705164*^9, 
   3.567516766225211*^9}, {3.567516802947675*^9, 3.567516829186921*^9}, {
   3.5675168735065994`*^9, 3.5675169172802763`*^9}, 3.56753241249549*^9, {
   3.5675324503567567`*^9, 3.567532480168409*^9}, 3.568099720232608*^9, 
   3.568101007827669*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"plot5", "=", 
    RowBox[{
     RowBox[{"Graphics3D", "[", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Sphere", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "1", "]"}], "]"}], ",", 
            RowBox[{"yc", "[", 
             RowBox[{"[", "1", "]"}], "]"}], ",", 
            RowBox[{"zc", "[", 
             RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "Cr"}], "]"}], ",", 
        RowBox[{"Sphere", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "2", "]"}], "]"}], ",", 
            RowBox[{"yc", "[", 
             RowBox[{"[", "2", "]"}], "]"}], ",", 
            RowBox[{"zc", "[", 
             RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "Nr"}], "]"}], ",", 
        RowBox[{"Sphere", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "3", "]"}], "]"}], ",", 
            RowBox[{"yc", "[", 
             RowBox[{"[", "3", "]"}], "]"}], ",", 
            RowBox[{"zc", "[", 
             RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "Oa"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Opacity", "[", "0.4", "]"}], ",", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"Lighter", "[", 
             RowBox[{"Green", ",", "0.3"}], "]"}], ",", 
            RowBox[{"Sphere", "[", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
                RowBox[{"xc", "[", 
                 RowBox[{"[", "10", "]"}], "]"}], ",", 
                RowBox[{"yc", "[", 
                 RowBox[{"[", "10", "]"}], "]"}], ",", 
                RowBox[{"zc", "[", 
                 RowBox[{"[", "10", "]"}], "]"}]}], "}"}], ",", "NM"}], "]"}],
             ",", 
            RowBox[{"Sphere", "[", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
                RowBox[{"xc", "[", 
                 RowBox[{"[", "11", "]"}], "]"}], ",", 
                RowBox[{"yc", "[", 
                 RowBox[{"[", "11", "]"}], "]"}], ",", 
                RowBox[{"zc", "[", 
                 RowBox[{"[", "11", "]"}], "]"}]}], "}"}], ",", "NM"}], "]"}],
             ",", 
            RowBox[{"Sphere", "[", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
                RowBox[{"xc", "[", 
                 RowBox[{"[", "12", "]"}], "]"}], ",", 
                RowBox[{"yc", "[", 
                 RowBox[{"[", "12", "]"}], "]"}], ",", 
                RowBox[{"zc", "[", 
                 RowBox[{"[", "12", "]"}], "]"}]}], "}"}], ",", "OC"}], "]"}],
             ",", 
            RowBox[{"Sphere", "[", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
                RowBox[{"xc", "[", 
                 RowBox[{"[", "13", "]"}], "]"}], ",", 
                RowBox[{"yc", "[", 
                 RowBox[{"[", "13", "]"}], "]"}], ",", 
                RowBox[{"zc", "[", 
                 RowBox[{"[", "13", "]"}], "]"}]}], "}"}], ",", "OC"}], "]"}],
             ",", 
            RowBox[{"Sphere", "[", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
                RowBox[{"xc", "[", 
                 RowBox[{"[", "14", "]"}], "]"}], ",", 
                RowBox[{"yc", "[", 
                 RowBox[{"[", "14", "]"}], "]"}], ",", 
                RowBox[{"zc", "[", 
                 RowBox[{"[", "14", "]"}], "]"}]}], "}"}], ",", "OO"}], "]"}],
             ",", 
            RowBox[{"Sphere", "[", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
                RowBox[{"xc", "[", 
                 RowBox[{"[", "15", "]"}], "]"}], ",", 
                RowBox[{"yc", "[", 
                 RowBox[{"[", "15", "]"}], "]"}], ",", 
                RowBox[{"zc", "[", 
                 RowBox[{"[", "15", "]"}], "]"}]}], "}"}], ",", "OO"}], 
             "]"}]}], "}"}]}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Darker", "[", 
           RowBox[{"Red", ",", "1"}], "]"}], ",", 
          RowBox[{"Sphere", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xn", "[", 
               RowBox[{"[", "1", "]"}], "]"}], ",", 
              RowBox[{"yn", "[", 
               RowBox[{"[", "1", "]"}], "]"}], ",", 
              RowBox[{"zn", "[", 
               RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
          ",", 
          RowBox[{"Sphere", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xn", "[", 
               RowBox[{"[", "2", "]"}], "]"}], ",", 
              RowBox[{"yn", "[", 
               RowBox[{"[", "2", "]"}], "]"}], ",", 
              RowBox[{"zn", "[", 
               RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
          ",", 
          RowBox[{"Sphere", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xn", "[", 
               RowBox[{"[", "3", "]"}], "]"}], ",", 
              RowBox[{"yn", "[", 
               RowBox[{"[", "3", "]"}], "]"}], ",", 
              RowBox[{"zn", "[", 
               RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
          ",", 
          RowBox[{"Sphere", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xn", "[", 
               RowBox[{"[", "4", "]"}], "]"}], ",", 
              RowBox[{"yn", "[", 
               RowBox[{"[", "4", "]"}], "]"}], ",", 
              RowBox[{"zn", "[", 
               RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
          ",", 
          RowBox[{"(*", 
           RowBox[{
            RowBox[{"Sphere", "[", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
                RowBox[{"xn", "[", 
                 RowBox[{"[", "5", "]"}], "]"}], ",", 
                RowBox[{"yn", "[", 
                 RowBox[{"[", "5", "]"}], "]"}], ",", 
                RowBox[{"zn", "[", 
                 RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "0.08"}], 
             "]"}], ","}], "*)"}], 
          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[{"yn", "[", 
                 RowBox[{"[", "2", "]"}], "]"}], ",", 
                RowBox[{"zn", "[", 
                 RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"xn", "[", 
                 RowBox[{"[", "1", "]"}], "]"}], ",", 
                RowBox[{"yn", "[", 
                 RowBox[{"[", "1", "]"}], "]"}], ",", 
                RowBox[{"zn", "[", 
                 RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", 
            "0.05"}], "]"}], ",", 
          RowBox[{"Cylinder", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
                RowBox[{"xn", "[", 
                 RowBox[{"[", "3", "]"}], "]"}], ",", 
                RowBox[{
                 RowBox[{"yn", "[", 
                  RowBox[{"[", "3", "]"}], "]"}], "-", "0.14"}], ",", 
                RowBox[{"zn", "[", 
                 RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"xn", "[", 
                 RowBox[{"[", "1", "]"}], "]"}], ",", 
                RowBox[{
                 RowBox[{"yn", "[", 
                  RowBox[{"[", "1", "]"}], "]"}], "-", "0.14"}], ",", 
                RowBox[{"zn", "[", 
                 RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", 
            "0.05"}], "]"}], ",", 
          RowBox[{"Cylinder", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
                RowBox[{"xn", "[", 
                 RowBox[{"[", "3", "]"}], "]"}], ",", 
                RowBox[{
                 RowBox[{"yn", "[", 
                  RowBox[{"[", "3", "]"}], "]"}], "+", "0.14"}], ",", 
                RowBox[{"zn", "[", 
                 RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"xn", "[", 
                 RowBox[{"[", "1", "]"}], "]"}], ",", 
                RowBox[{
                 RowBox[{"yn", "[", 
                  RowBox[{"[", "1", "]"}], "]"}], "+", "0.14"}], ",", 
                RowBox[{"zn", "[", 
                 RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", 
            "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[{"[", "4", "]"}], "]"}], ",", 
                RowBox[{"yn", "[", 
                 RowBox[{"[", "4", "]"}], "]"}], ",", 
                RowBox[{"zn", "[", 
                 RowBox[{"[", "4", "]"}], "]"}]}], "}"}]}], "}"}], ",", 
            "0.05"}], "]"}], ",", 
          RowBox[{"(*", 
           RowBox[{
            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[{"[", "1", "]"}], "]"}], ",", 
                RowBox[{"yn", "[", 
                 RowBox[{"[", "1", "]"}], "]"}], ",", 
                RowBox[{"zn", "[", 
                 RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", 
            "0.05"}], "]"}]}], "}"}]}], "}"}], "]"}], "/.", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\n", "\[IndentingNewLine]", 
  StyleBox[
   RowBox[{"(*", " ", 
    RowBox[{
     RowBox[{"Overlapping", " ", "\[Pi]", " ", "system"}], ",", " ", 
     RowBox[{
     "forcing", " ", "the", " ", "N", " ", "atom", " ", "and", " ", "NH", " ",
       "bonds", " ", "into", " ", "plane", " ", "of", " ", "all", " ", "6", 
      " ", "atoms"}]}], " ", "*)"}], "Text"]}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{"plot5", ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", 
     RowBox[{"Boxed", "\[Rule]", "False"}]}], "}"}]}], 
  "]"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.567343604714074*^9, 3.5673436277865143`*^9}, {
   3.5673436837438126`*^9, 3.5673437262850876`*^9}, {3.5673438309300714`*^9, 
   3.5673438393072863`*^9}, {3.567343886060568*^9, 3.5673439230794334`*^9}, {
   3.567344408021885*^9, 3.56734451319727*^9}, 3.5673456889869347`*^9, 
   3.5673459664490223`*^9, {3.5673460936360455`*^9, 3.5673461706689806`*^9}, {
   3.567511905205871*^9, 3.5675119287307124`*^9}}],

Cell[BoxData[
 Graphics3DBox[{SphereBox[{0, 0, 0}, 0.25559559017708283`], 
   SphereBox[{-2.4981167513639795`, 0, 0}, 0.18966729854432163`], 
   SphereBox[{0.7843483199037025, 2.1015683080714793`, 0}, 
    0.18771946482544016`], 
   {Opacity[0.4], 
    {RGBColor[0.3, 1., 0.3], 
     SphereBox[{-2.4981167513639795`, 0, 0.8768151106038778}, 
      0.8768151106038778], 
     SphereBox[{-2.4981167513639795`, 0, -0.8768151106038778}, 
      0.8768151106038778], 
     SphereBox[{0, 0, 1.3099009653202152`}, 1.3099009653202152`], 
     SphereBox[{0, 0, -1.3099009653202152`}, 1.3099009653202152`], 
     SphereBox[{0.7843483199037025, 2.1015683080714793`, 1.1059447491867593`},
       1.1059447491867593`], 
     SphereBox[{0.7843483199037025, 
      2.1015683080714793`, -1.1059447491867593`}, 1.1059447491867593`]}}, 
   {RGBColor[0, 0, 0], SphereBox[{0, 0, 0}, 0.08], 
    SphereBox[{-2.4981167513639795`, 0, 0}, 0.08], 
    SphereBox[{0.7843483199037025, 2.1015683080714793`, 0}, 0.08], 
    SphereBox[{-3.2598496976882316`, 1.6241829959157847`, 0}, 0.08], 
    SphereBox[{1.3590493615167112`, -1.5422575918613213`, 0}, 0.08], 
    CylinderBox[{{-2.4981167513639795`, 0, 0}, {0, 0, 0}}, 0.05], 
    CylinderBox[{{0.7843483199037025, 1.9615683080714792`, 0}, {0, -0.14, 0}},
      0.05], CylinderBox[{{0.7843483199037025, 2.2415683080714794`, 0}, {
      0, 0.14, 0}}, 0.05], 
    CylinderBox[{{-2.4981167513639795`, 0, 0}, {-3.2598496976882316`, 
       1.6241829959157847`, 0}}, 0.05], 
    CylinderBox[{{1.3590493615167112`, -1.5422575918613213`, 0}, {0, 0, 0}}, 
     0.05]}},
  AspectRatio->Automatic,
  AutomaticImageSize->True,
  Boxed->False,
  ImageSize->{388.9407521890828, 382.5739513228262},
  ViewPoint->{0.06107555653723071, -0.30484743135150033`, 3.3694714451961256`},
  ViewVertical->{-0.5011232630071526, 0.7091910856939899, 
   0.6560171231862861}]], "Output",
 CellChangeTimes->{
  3.5673462864835844`*^9, 3.5675104564157267`*^9, 3.567510496523397*^9, 
   3.5675106615716867`*^9, 3.56751085663443*^9, 3.567510948534191*^9, 
   3.567511039451151*^9, 3.567511182659402*^9, 3.5675113146512337`*^9, 
   3.5675114624614935`*^9, 3.5675114996363587`*^9, 3.567511714511136*^9, 
   3.567511977199998*^9, {3.5675120720949645`*^9, 3.5675120853861876`*^9}, 
   3.5675122601220946`*^9, 3.5675138870177517`*^9, 3.5675140662464666`*^9, 
   3.5675161992734146`*^9, 3.5675162593335204`*^9, 3.5675162957751846`*^9, 
   3.567516360000497*^9, {3.567516406738179*^9, 3.5675164774219036`*^9}, 
   3.5675165357036057`*^9, 3.5675165970585136`*^9, 3.56751662927257*^9, 
   3.5675166616426272`*^9, 3.5675167056971045`*^9, {3.567516739720764*^9, 
   3.567516766240811*^9}, {3.567516802963275*^9, 3.5675168292181215`*^9}, {
   3.567516873522199*^9, 3.567516917295876*^9}, 3.5675324125110903`*^9, {
   3.5675324503723564`*^9, 3.567532480184009*^9}, 3.568099720263808*^9, 
   3.5681010078432693`*^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"}], "]"}], 
         ",", 
         RowBox[{"Sphere", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"xc", "[", 
              RowBox[{"[", "2", "]"}], "]"}], ",", 
             RowBox[{"yc", "[", 
              RowBox[{"[", "2", "]"}], "]"}], ",", 
             RowBox[{"zc", "[", 
              RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "Nr"}], "]"}], 
         ",", 
         RowBox[{"Sphere", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"xc", "[", 
              RowBox[{"[", "3", "]"}], "]"}], ",", 
             RowBox[{"yc", "[", 
              RowBox[{"[", "3", "]"}], "]"}], ",", 
             RowBox[{"zc", "[", 
              RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "Oa"}], "]"}], 
         ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"Opacity", "[", "0.6", "]"}], ",", 
           RowBox[{"Sphere", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"xc", "[", 
                RowBox[{"[", "4", "]"}], "]"}], ",", 
               RowBox[{"yc", "[", 
                RowBox[{"[", "4", "]"}], "]"}], ",", 
               RowBox[{"zc", "[", 
                RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "NH"}], "]"}], 
           ",", 
           RowBox[{"Sphere", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"xc", "[", 
                RowBox[{"[", "5", "]"}], "]"}], ",", 
               RowBox[{"yc", "[", 
                RowBox[{"[", "5", "]"}], "]"}], ",", 
               RowBox[{"zc", "[", 
                RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "NH"}], "]"}], 
           ",", 
           RowBox[{"Sphere", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"xc", "[", 
                RowBox[{"[", "6", "]"}], "]"}], ",", 
               RowBox[{"yc", "[", 
                RowBox[{"[", "6", "]"}], "]"}], ",", 
               RowBox[{"zc", "[", 
                RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], 
           ",", 
           RowBox[{"Sphere", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"xc", "[", 
                RowBox[{"[", "7", "]"}], "]"}], ",", 
               RowBox[{"yc", "[", 
                RowBox[{"[", "7", "]"}], "]"}], ",", 
               RowBox[{"zc", "[", 
                RowBox[{"[", "7", "]"}], "]"}]}], "}"}], ",", "CN"}], "]"}], 
           ",", 
           RowBox[{"Sphere", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"xc", "[", 
                RowBox[{"[", "8", "]"}], "]"}], ",", 
               RowBox[{"yc", "[", 
                RowBox[{"[", "8", "]"}], "]"}], ",", 
               RowBox[{"zc", "[", 
                RowBox[{"[", "8", "]"}], "]"}]}], "}"}], ",", "CO"}], "]"}]}],
           "}"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"Darker", "[", 
            RowBox[{"Red", ",", "1"}], "]"}], ",", 
           RowBox[{"Sphere", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"xn", "[", 
                RowBox[{"[", "1", "]"}], "]"}], ",", 
               RowBox[{"yn", "[", 
                RowBox[{"[", "1", "]"}], "]"}], ",", 
               RowBox[{"zn", "[", 
                RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}],
            ",", 
           RowBox[{"Sphere", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"xn", "[", 
                RowBox[{"[", "2", "]"}], "]"}], ",", 
               RowBox[{"yn", "[", 
                RowBox[{"[", "2", "]"}], "]"}], ",", 
               RowBox[{"zn", "[", 
                RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}],
            ",", 
           RowBox[{"Sphere", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"xn", "[", 
                RowBox[{"[", "3", "]"}], "]"}], ",", 
               RowBox[{"yn", "[", 
                RowBox[{"[", "3", "]"}], "]"}], ",", 
               RowBox[{"zn", "[", 
                RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}],
            ",", 
           RowBox[{"Sphere", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"xn", "[", 
                RowBox[{"[", "4", "]"}], "]"}], ",", 
               RowBox[{"yn", "[", 
                RowBox[{"[", "4", "]"}], "]"}], ",", 
               RowBox[{"zn", "[", 
                RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}],
            ",", 
           RowBox[{"(*", 
            RowBox[{
             RowBox[{"Sphere", "[", 
              RowBox[{
               RowBox[{"{", 
                RowBox[{
                 RowBox[{"xn", "[", 
                  RowBox[{"[", "5", "]"}], "]"}], ",", 
                 RowBox[{"yn", "[", 
                  RowBox[{"[", "5", "]"}], "]"}], ",", 
                 RowBox[{"zn", "[", 
                  RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "0.08"}], 
              "]"}], ","}], "*)"}], 
           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[{"yn", "[", 
                  RowBox[{"[", "2", "]"}], "]"}], ",", 
                 RowBox[{"zn", "[", 
                  RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", 
               RowBox[{"{", 
                RowBox[{
                 RowBox[{"xn", "[", 
                  RowBox[{"[", "1", "]"}], "]"}], ",", 
                 RowBox[{"yn", "[", 
                  RowBox[{"[", "1", "]"}], "]"}], ",", 
                 RowBox[{"zn", "[", 
                  RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", 
             "0.05"}], "]"}], ",", 
           RowBox[{"Cylinder", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"{", 
                RowBox[{
                 RowBox[{"xn", "[", 
                  RowBox[{"[", "3", "]"}], "]"}], ",", 
                 RowBox[{
                  RowBox[{"yn", "[", 
                   RowBox[{"[", "3", "]"}], "]"}], "-", "0.14"}], ",", 
                 RowBox[{"zn", "[", 
                  RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", 
               RowBox[{"{", 
                RowBox[{
                 RowBox[{"xn", "[", 
                  RowBox[{"[", "1", "]"}], "]"}], ",", 
                 RowBox[{
                  RowBox[{"yn", "[", 
                   RowBox[{"[", "1", "]"}], "]"}], "-", "0.14"}], ",", 
                 RowBox[{"zn", "[", 
                  RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", 
             "0.05"}], "]"}], ",", 
           RowBox[{"Cylinder", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"{", 
                RowBox[{
                 RowBox[{"xn", "[", 
                  RowBox[{"[", "3", "]"}], "]"}], ",", 
                 RowBox[{
                  RowBox[{"yn", "[", 
                   RowBox[{"[", "3", "]"}], "]"}], "+", "0.14"}], ",", 
                 RowBox[{"zn", "[", 
                  RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", 
               RowBox[{"{", 
                RowBox[{
                 RowBox[{"xn", "[", 
                  RowBox[{"[", "1", "]"}], "]"}], ",", 
                 RowBox[{
                  RowBox[{"yn", "[", 
                   RowBox[{"[", "1", "]"}], "]"}], "+", "0.14"}], ",", 
                 RowBox[{"zn", "[", 
                  RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", 
             "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[{"[", "4", "]"}], "]"}], ",", 
                 RowBox[{"yn", "[", 
                  RowBox[{"[", "4", "]"}], "]"}], ",", 
                 RowBox[{"zn", "[", 
                  RowBox[{"[", "4", "]"}], "]"}]}], "}"}]}], "}"}], ",", 
             "0.05"}], "]"}], ",", 
           RowBox[{"(*", 
            RowBox[{
             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[{"[", "1", "]"}], "]"}], ",", 
                 RowBox[{"yn", "[", 
                  RowBox[{"[", "1", "]"}], "]"}], ",", 
                 RowBox[{"zn", "[", 
                  RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", 
             "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.5675119422403364`*^9, 3.567511964267575*^9}}],

Cell[BoxData[
 Graphics3DBox[{SphereBox[{0, 0, 0}, 0.25559559017708283`], 
   SphereBox[{-2.4981167513639795`, 0, 0}, 0.18966729854432163`], 
   SphereBox[{0.7843483199037025, 2.1015683080714793`, 0}, 
    0.18771946482544016`], 
   {Opacity[0.6], 
    SphereBox[{-3.0195425457907588`, 1.1117950366024378`, 0}, 
     1.0383282480711433`], 
    SphereBox[{-3.0195425457907588`, -1.1117950366024378`, 0}, 
     1.0383282480711433`], 
    SphereBox[{0.9869942279020033, -1.1200471330977853`, 0}, 
     1.237275196687085], 
    SphereBox[{-1.2820225214983703`, 0, 0}, 1.0264269313212875`], 
    SphereBox[{0.40404098920022924`, 1.0825799157309375`, 0}, 
     0.8999252425636789]}, 
   {RGBColor[0, 0, 0], SphereBox[{0, 0, 0}, 0.08], 
    SphereBox[{-2.4981167513639795`, 0, 0}, 0.08], 
    SphereBox[{0.7843483199037025, 2.1015683080714793`, 0}, 0.08], 
    SphereBox[{-3.2598496976882316`, 1.6241829959157847`, 0}, 0.08], 
    SphereBox[{1.3590493615167112`, -1.5422575918613213`, 0}, 0.08], 
    CylinderBox[{{-2.4981167513639795`, 0, 0}, {0, 0, 0}}, 0.05], 
    CylinderBox[{{0.7843483199037025, 1.9615683080714792`, 0}, {0, -0.14, 0}},
      0.05], CylinderBox[{{0.7843483199037025, 2.2415683080714794`, 0}, {
      0, 0.14, 0}}, 0.05], 
    CylinderBox[{{-2.4981167513639795`, 0, 0}, {-3.2598496976882316`, 
       1.6241829959157847`, 0}}, 0.05], 
    CylinderBox[{{1.3590493615167112`, -1.5422575918613213`, 0}, {0, 0, 0}}, 
     0.05]}},
  AspectRatio->Automatic,
  AutomaticImageSize->True,
  Boxed->False,
  ImageSize->{295.60962735440614`, 256.66852778600247`},
  ViewPoint->{0.30364390661022517`, 0.5601352233907192, 3.3232587785930106`},
  ViewVertical->{0.017263820249459494`, 1.0679468664157148`, 
   1.0510539576658122`}]], "Output",
 CellChangeTimes->{
  3.567346286530384*^9, 3.567510456431327*^9, 3.5675104965389967`*^9, 
   3.567510661587287*^9, 3.5675108566500297`*^9, 3.567510948549791*^9, 
   3.5675110394667506`*^9, 3.567511182675002*^9, 3.567511314666834*^9, 
   3.5675114624770937`*^9, 3.5675114996519585`*^9, 3.5675117145267363`*^9, 
   3.567511977262398*^9, {3.567512072110564*^9, 3.5675120854017878`*^9}, 
   3.5675122601532946`*^9, 3.567513887033352*^9, 3.567514066262067*^9, 
   3.567516199289015*^9, 3.56751625934912*^9, 3.5675162957907844`*^9, 
   3.567516360016097*^9, {3.5675164067537794`*^9, 3.5675164774531035`*^9}, 
   3.567516535719206*^9, 3.567516597074114*^9, 3.56751662928817*^9, 
   3.567516661673827*^9, 3.5675167057127047`*^9, {3.5675167397363644`*^9, 
   3.5675167662564106`*^9}, {3.567516802978875*^9, 3.5675168292337217`*^9}, {
   3.567516873537799*^9, 3.5675169173114758`*^9}, 3.56753241252669*^9, {
   3.5675324503879566`*^9, 3.567532480199609*^9}, 3.568099720279408*^9, 
   3.5681010078588696`*^9}]
}, Open  ]]
},
WindowToolbars->"EditBar",
WindowSize->{836, 964},
WindowMargins->{{Automatic, 527}, {Automatic, -3}},
PrintingCopies->1,
PrintingPageRange->{5, 6},
PrintingOptions->{"Magnification"->1.,
"PaperOrientation"->"Portrait",
"PaperSize"->{595.2, 841.8}},
PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}},
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (November 20, 2012)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[557, 20, 2803, 41, 31, "Input"],
Cell[3363, 63, 5876, 89, 235, "Input"],
Cell[9242, 154, 1319, 26, 72, "Input"],
Cell[10564, 182, 804, 23, 252, "Input"],
Cell[11371, 207, 4743, 68, 99, "Input"],
Cell[16117, 277, 5668, 80, 167, "Input"],
Cell[21788, 359, 5105, 72, 167, "Input"],
Cell[26896, 433, 4639, 66, 82, "Input"],
Cell[31538, 501, 4932, 69, 82, "Input"],
Cell[36473, 572, 4667, 67, 82, "Input"],
Cell[41143, 641, 605, 15, 133, "Input"],
Cell[41751, 658, 384, 11, 99, "Input"],
Cell[42138, 671, 571, 14, 116, "Input"],
Cell[CellGroupData[{
Cell[42734, 689, 2072, 31, 116, "Input"],
Cell[44809, 722, 2189, 38, 72, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[47035, 765, 988, 24, 218, "Input"],
Cell[48026, 791, 1041, 15, 31, "Output"],
Cell[49070, 808, 1058, 16, 31, "Output"],
Cell[50131, 826, 1041, 15, 31, "Output"],
Cell[51175, 843, 1041, 15, 31, "Output"],
Cell[52219, 860, 1041, 15, 31, "Output"],
Cell[53263, 877, 1042, 15, 31, "Output"],
Cell[54308, 894, 1041, 15, 31, "Output"],
Cell[55352, 911, 1041, 15, 31, "Output"],
Cell[56396, 928, 1042, 15, 31, "Output"],
Cell[57441, 945, 1042, 15, 31, "Output"],
Cell[58486, 962, 1042, 15, 31, "Output"],
Cell[59531, 979, 1042, 15, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[60610, 999, 2907, 47, 354, "Input"],
Cell[63520, 1048, 2981, 50, 361, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[66538, 1103, 1897, 30, 235, "Input"],
Cell[68438, 1135, 2461, 41, 303, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[70936, 1181, 4126, 74, 643, "Input"],
Cell[75065, 1257, 3879, 77, 445, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[78981, 1339, 3008, 67, 152, "Input"],
Cell[81992, 1408, 2007, 41, 116, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[84036, 1454, 2050, 47, 72, "Input"],
Cell[86089, 1503, 2372, 42, 116, "Output"]
}, Open  ]],
Cell[88476, 1548, 11202, 292, 472, "Input"],
Cell[CellGroupData[{
Cell[99703, 1844, 843, 22, 52, "Input"],
Cell[100549, 1868, 2485, 46, 289, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[103071, 1919, 13025, 337, 532, "Input"],
Cell[116099, 2258, 2890, 51, 398, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[119026, 2314, 12420, 321, 472, "Input"],
Cell[131449, 2637, 2745, 49, 272, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)