(* 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[     29428,        704]
NotebookOptionsPosition[     28272,        664]
NotebookOutlinePosition[     28785,        683]
CellTagsIndexPosition[     28742,        680]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"O2", " ", 
    RowBox[{"with", ":", "  ", 
     RowBox[{"12.05", ".2013"}]}]}], " ", "*)"}], "\n", 
  RowBox[{
   RowBox[{
    RowBox[{"Clear", "[", 
     RowBox[{
     "a", ",", "z", ",", "sig1", ",", "sig2", ",", "sig3", ",", "sig4", ",", 
      "k1", ",", "k2", ",", "k3", ",", "k4", ",", "k5", ",", "nc", ",", "vee",
       ",", "vne", ",", "vnn", ",", "R1", ",", "R2", ",", "R3", ",", "R4", 
      ",", "R5", ",", "\n", "xc", ",", "yc", ",", "zc", ",", "xn", ",", "yn", 
      ",", "zn", ",", "ow", ",", "cn", ",", "rr", ",", "d1", ",", "d2", ",", 
      "d3", ",", "d4", ",", "d5", ",", "i", ",", "j", ",", "d", ",", "t", ",",
       "sq", ",", "Ekin", ",", "Epot"}], "]"}], ";"}], "\n", 
   RowBox[{
    RowBox[{"z", "=", "7.0"}], ";"}], "\n", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"pi", "=", "0.0"}], ";"}], "*)"}], " ", 
   RowBox[{"(*", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"distance", " ", "of", " ", "bananas", " ", "from", " ", "C"}], 
      "-", 
      RowBox[{"C", " ", "axis"}]}], ",", " ", 
     RowBox[{"pi", " ", "=", 
      RowBox[{"0", " ", "is", " ", "regular", " ", "tetraedric"}]}]}], " ", 
    "*)"}], "\n", 
   RowBox[{
    RowBox[{"nc", "=", "9"}], ";", "  ", 
    RowBox[{"(*", " ", 
     RowBox[{"number", " ", "of", " ", "clouds"}], " ", "*)"}], "\n", 
    RowBox[{"sig1", "=", "0.3"}], ";", " ", 
    RowBox[{"sig2", "=", "0.3"}], ";", " ", 
    RowBox[{"sig3", "=", "sig2"}], ";"}], " ", 
   RowBox[{"(*", " ", 
    RowBox[{
     RowBox[{"screening", " ", 
      RowBox[{"const", ".", " ", "from"}], " ", "e"}], "-", 
     RowBox[{"e", " ", "interaction", " ", "in", " ", "doubly", " ", 
      RowBox[{"occ", ".", " ", "clouds"}]}]}], " ", "*)"}], "\n", 
   RowBox[{
    RowBox[{"k1", "=", "0.8"}], ";", " ", 
    RowBox[{"k2", "=", "0.9"}], ";", " ", 
    RowBox[{"k3", "=", "0.9"}], ";", " ", 
    RowBox[{"k4", "=", "1.0"}], ";", " ", 
    RowBox[{"k5", "=", "1.0"}], ";", 
    RowBox[{"(*", " ", 
     RowBox[{
      RowBox[{
      "parameters", " ", "for", " ", "kinetic", " ", "energy", " ", "of", " ",
        "clouds"}], ";", " ", 
      RowBox[{"k", "=", 
       RowBox[{"1.0", " ", 
        RowBox[{"Kimball", "'"}], "s", " ", "lowest", " ", "value"}]}]}], " ",
      "*)"}], "\n", 
    RowBox[{"bohr", "=", "0.529177"}], ";"}], " ", "\n", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"R1", "=", "0.18"}], ";", " ", 
     RowBox[{"R2", "=", "1.0"}], ";", " ", 
     RowBox[{"R3", "=", "0.96"}], ";", " ", 
     RowBox[{"R4", "=", "1.0"}], ";", " ", 
     RowBox[{"R5", "=", "1.0"}], ";"}], "*)"}]}]}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5773415047046976`*^9, 3.5773415137059135`*^9}, {
   3.5773416513449545`*^9, 3.577341669752987*^9}, {3.5773417239318824`*^9, 
   3.577341893956581*^9}, 3.5773419407566633`*^9, {3.5773420175555983`*^9, 
   3.577342056165666*^9}, {3.5773420929973307`*^9, 3.5773421056489525`*^9}, {
   3.5773421396414127`*^9, 3.5773423470593767`*^9}, {3.577342382081438*^9, 
   3.577342441158742*^9}, {3.5773425028880506`*^9, 3.577342523682887*^9}, {
   3.5773425980794177`*^9, 3.5773426626791315`*^9}, {3.5773426933487854`*^9, 
   3.5773426941131864`*^9}, {3.577342738354864*^9, 3.577342889035529*^9}, {
   3.5773429412644205`*^9, 3.57734295234044*^9}, {3.577342991824109*^9, 
   3.5773430231333647`*^9}, {3.5773430561274223`*^9, 3.577343143737176*^9}, {
   3.5773431993668737`*^9, 3.577343258334977*^9}, 3.5773434426249013`*^9, {
   3.577343480042567*^9, 3.5773435125686235`*^9}, {3.5773436063903885`*^9, 
   3.5773436272008247`*^9}, {3.5773436617392855`*^9, 
   3.5773437808610945`*^9}, {3.577343868193248*^9, 3.577343879659268*^9}, {
   3.5773439214829416`*^9, 3.577344118838888*^9}, {3.5773441503509436`*^9, 
   3.577344170740179*^9}, {3.5773446920694942`*^9, 3.577344693286297*^9}, 
   3.5773526894956493`*^9, {3.5773528388503113`*^9, 3.577352838897112*^9}, {
   3.577353457328998*^9, 3.577353503364679*^9}, {3.5773556679320807`*^9, 
   3.5773557165573664`*^9}, 3.5773607112650003`*^9, {3.577363145384075*^9, 
   3.5773631641821084`*^9}, 3.5773713130222435`*^9, {3.5773735851350346`*^9, 
   3.577373695474028*^9}, {3.577373964730501*^9, 3.5773739677725067`*^9}, 
   3.5773740596254673`*^9, {3.577374092245125*^9, 3.5773741045067463`*^9}, {
   3.5773741853616886`*^9, 3.577374230414568*^9}, 3.5773742637986264`*^9, {
   3.5773746139412413`*^9, 3.5773746339092765`*^9}, {3.5773750495252066`*^9, 
   3.5773750633000307`*^9}, 3.5773754561555204`*^9, {3.5773757271123962`*^9, 
   3.577375728063998*^9}, 3.577375963078411*^9, {3.5773761205918875`*^9, 
   3.5773762048008356`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"C", " ", "He"}], "-", "shells"}], " ", "*)"}], "\n", 
  RowBox[{Cell["\<\
(* N,N He shell*)
T = 2*2.25*k1/R1^2;
vee=2*3.0*sig1/R1;
vne=-2*3.0*z/R1;

a=Table[{0,0,0},{nc}];
b=Table[{0,0,0},{nc}];
c=Table[{0,0,0},{nc}];
(* bonding and lone pairs *)
T = T + 2.25*(2*k2/R2^2+k3/R3^2)+4*1.125*k4/R4^2);
vee=vee+3.0*(2*sig2/R2+sig3/R3);
c[[1,1]]=-2;      c[[1,2]]=z;  c[[1,3]]=R1;
c[[2,1]]=-2;      c[[2,2]]=z;  c[[2,3]]=R1;
c[[3,1]]=-2;      c[[3,2]]=0;  c[[3,3]]=R2;
c[[4,1]]=-2;      c[[4,2]]=0;  c[[4,3]]=R2;
c[[5,1]]=-2;      c[[5,2]]=0;  c[[5,3]]=R3;
c[[6,1]]=-1;      c[[6,2]]=0;  c[[6,3]]=R4;
c[[7,1]]=-1;      c[[7,2]]=0;  c[[7,3]]=R4;
c[[8,1]]=-1;      c[[8,2]]=0;  c[[8,3]]=R4;
c[[9,1]]=-1;      c[[9,2]]=0;  c[[9,3]]=R4;
d1=R1+R3;              b[[1,1]]=d1; 
d2=R1+R3+R1+R2;        b[[2,1]]=-d1;
                       a[[1,1]]=d1;
                       a[[2,1]]=-d1;
                       a[[3,1]]=d2;
\t\t       a[[4,1]]=-d2;
                       a[[5,1]]=0.0;             
                       a[[6,1]]=d1;     a[[6,2]]=R4;
                       a[[7,1]]=d1;                 a[[7,3]]=-R4;
                       a[[8,1]]=-d1;    a[[8,2]]=R4;
\t\t       a[[9,1]]=-d1;                 a[[9,3]]=-R4
For[ i=1,i<nc,i++,
  For[ j=i+1,j<=nc,j++,
    d=Sqrt[(a[[i,1]]-a[[j,1]])^2+(a[[i,2]]-a[[j,2]])^2+(a[[i,3]]-a[[j,3]])^2];
    If[ d<>0, vee=vee+c[[i,1]]*c[[j,1]]/d;
vnn=z*z/(2*d1);
For[ i=1,i<=nc,i++,
  For[ j=i+1,j<=nc,j++,
    If[ i<>j,
      d=Sqrt[(a[[i,1]]-b[[j,1]])^2+(a[[i,2]]-b[[j,2]]^2)+(a[[i,3]]-b[[j,3]])^\
2];
      If[ d<>0, vne=vne+c[[i,1]]*c[[j,2]]/d;
Ekin=T;
Epot=vne+vee+vnn;
func=Ekin+Epot;\
\>", "Input"], "\n", 
   RowBox[{"t", " ", "=", " ", 
    RowBox[{"FindMinimum", "[", 
     RowBox[{"func", ",", 
      RowBox[{"{", 
       RowBox[{"R1", ",", "0.21"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"R2", ",", "0.8"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"R3", ",", "0.9"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"R4", ",", "0.9"}], "}"}], ",", 
      RowBox[{"Method", " ", "\[Rule]", "\"\<Newton\>\""}], ",", " ", 
      RowBox[{"MaxIterations", " ", "->", " ", "500"}]}], 
     "]"}]}]}]}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5773415047046976`*^9, 3.5773415137059135`*^9}, {
   3.5773416513449545`*^9, 3.577341669752987*^9}, {3.5773417239318824`*^9, 
   3.577341893956581*^9}, 3.5773419407566633`*^9, {3.5773420175555983`*^9, 
   3.577342056165666*^9}, {3.5773420929973307`*^9, 3.5773421056489525`*^9}, {
   3.5773421396414127`*^9, 3.5773423470593767`*^9}, {3.577342382081438*^9, 
   3.577342441158742*^9}, {3.5773425028880506`*^9, 3.577342523682887*^9}, {
   3.5773425980794177`*^9, 3.5773426626791315`*^9}, {3.5773426933487854`*^9, 
   3.5773426941131864`*^9}, {3.577342738354864*^9, 3.577342889035529*^9}, {
   3.5773429412644205`*^9, 3.57734295234044*^9}, {3.577342991824109*^9, 
   3.5773430231333647`*^9}, {3.5773430561274223`*^9, 3.577343143737176*^9}, {
   3.5773431993668737`*^9, 3.577343258334977*^9}, 3.5773434426249013`*^9, {
   3.577343480042567*^9, 3.5773435125686235`*^9}, {3.5773436063903885`*^9, 
   3.5773436272008247`*^9}, {3.5773436617392855`*^9, 
   3.5773437808610945`*^9}, {3.577343868193248*^9, 3.577343879659268*^9}, {
   3.5773439214829416`*^9, 3.577344118838888*^9}, {3.5773441503509436`*^9, 
   3.577344170740179*^9}, {3.5773446920694942`*^9, 3.577344693286297*^9}, 
   3.5773526894956493`*^9, {3.577352838897112*^9, 3.577352838928312*^9}, 
   3.577362343573867*^9, {3.577362462570876*^9, 3.577362524440585*^9}, 
   3.5773625762794757`*^9, {3.5773626257315626`*^9, 3.577362806021079*^9}, {
   3.5773631418428693`*^9, 3.577363185554146*^9}, 3.5773632603562775`*^9, {
   3.577371455731294*^9, 3.5773714994893713`*^9}, {3.5773724898259106`*^9, 
   3.577372494162718*^9}, {3.5773730043772144`*^9, 3.5773730238460484`*^9}, {
   3.5773731837619295`*^9, 3.577373238128025*^9}, {3.577373271434083*^9, 
   3.577373393987899*^9}, {3.577373430913163*^9, 3.5773734337835684`*^9}, {
   3.577373480833251*^9, 3.5773734946236753`*^9}, {3.5773739993625617`*^9, 
   3.5773740373018284`*^9}, {3.5773743739348197`*^9, 
   3.5773743842152376`*^9}, {3.5773744685645857`*^9, 3.577374501527444*^9}, {
   3.5773746732993455`*^9, 3.577374707073405*^9}, {3.577375258362773*^9, 
   3.577375352243738*^9}, 3.577375468307942*^9, {3.577375746784031*^9, 
   3.577375845345004*^9}, {3.577375895483492*^9, 3.5773759185091324`*^9}, 
   3.577376045072155*^9, 3.5773761022150555`*^9, {3.577376336184266*^9, 
   3.577376368523123*^9}, 3.5773799042893333`*^9, {3.577379940106996*^9, 
   3.577380154014572*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "112.224999145135`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"R1", "\[Rule]", "0.1791044776119398`"}], ",", 
     RowBox[{"R2", "\[Rule]", "7.725107497224924`*^7"}], ",", 
     RowBox[{"R3", "\[Rule]", "1.1206040075582859`*^8"}], ",", 
     RowBox[{"R4", "\[Rule]", "7.36327469994904`*^7"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.5773760631369867`*^9, {3.5773761074254646`*^9, 3.5773761538199463`*^9}, {
   3.5773761949728184`*^9, 3.577376209590044*^9}, 3.577376415260805*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", "func", " ", "*)"}], "\n", 
  RowBox[{
   RowBox[{"vne", "  ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}], " ", "\n", 
   RowBox[{"vee", "  ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}], " ", "\n", 
   RowBox[{"vnn", "  ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}], "\n", 
   RowBox[{"Vir", " ", "=", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "Epot"}], "/", "Ekin"}], "  ", "/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}]}], "\n", 
   RowBox[{
    RowBox[{"2", "*", "d1", "*", "bohr"}], "  ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}]}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5773415047046976`*^9, 3.5773415137059135`*^9}, {
   3.5773416513449545`*^9, 3.577341669752987*^9}, {3.5773417239318824`*^9, 
   3.577341893956581*^9}, 3.5773419407566633`*^9, {3.5773420175555983`*^9, 
   3.577342056165666*^9}, {3.5773420929973307`*^9, 3.5773421056489525`*^9}, {
   3.5773421396414127`*^9, 3.5773423470593767`*^9}, {3.577342382081438*^9, 
   3.577342441158742*^9}, {3.5773425028880506`*^9, 3.577342523682887*^9}, {
   3.5773425980794177`*^9, 3.5773426626791315`*^9}, {3.5773426933487854`*^9, 
   3.5773426941131864`*^9}, {3.577342738354864*^9, 3.577342889035529*^9}, {
   3.5773429412644205`*^9, 3.57734295234044*^9}, {3.577342991824109*^9, 
   3.5773430231333647`*^9}, {3.5773430561274223`*^9, 3.577343143737176*^9}, {
   3.5773431993668737`*^9, 3.577343258334977*^9}, 3.5773434426249013`*^9, {
   3.577343480042567*^9, 3.5773435125686235`*^9}, {3.5773436063903885`*^9, 
   3.5773436272008247`*^9}, {3.5773436617392855`*^9, 
   3.5773437808610945`*^9}, {3.577343868193248*^9, 3.577343879659268*^9}, {
   3.5773439214829416`*^9, 3.577344086609232*^9}, {3.5773444013153844`*^9, 
   3.5773444098953996`*^9}, {3.5773527214601054`*^9, 3.577352754235763*^9}, {
   3.577352839084312*^9, 3.577352839115512*^9}, {3.5773562861299667`*^9, 
   3.577356291496376*^9}, {3.5773717712106485`*^9, 3.5773718737028284`*^9}}],

Cell[BoxData[
 RowBox[{"-", "234.50000000000065`"}]], "Output",
 CellChangeTimes->{
  3.577376063152587*^9, {3.5773761074410644`*^9, 3.577376153835546*^9}, {
   3.577376194988418*^9, 3.577376209605644*^9}, 3.5773764152764053`*^9}],

Cell[BoxData["10.050000636232953`"], "Output",
 CellChangeTimes->{
  3.577376063152587*^9, {3.5773761074410644`*^9, 3.577376153835546*^9}, {
   3.577376194988418*^9, 3.577376209605644*^9}, 3.5773764152764053`*^9}],

Cell[BoxData["2.1863209301049814`*^-7"], "Output",
 CellChangeTimes->{
  3.577376063152587*^9, {3.5773761074410644`*^9, 3.577376153835546*^9}, {
   3.577376194988418*^9, 3.577376209605644*^9}, 3.5773764152764053`*^9}],

Cell[BoxData["1.999999992382574`"], "Output",
 CellChangeTimes->{
  3.577376063152587*^9, {3.5773761074410644`*^9, 3.577376153835546*^9}, {
   3.577376194988418*^9, 3.577376209605644*^9}, 3.5773764152764053`*^9}],

Cell[BoxData["1.1859957357109015`*^8"], "Output",
 CellChangeTimes->{
  3.577376063152587*^9, {3.5773761074410644`*^9, 3.577376153835546*^9}, {
   3.577376194988418*^9, 3.577376209605644*^9}, 3.5773764152764053`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"rr", "=", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
     "R1", ",", "R1", ",", "R2", ",", "R2", ",", "R3", ",", "R4", ",", "R4", 
      ",", "R4", ",", "R4"}], "}"}], " ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5773415047046976`*^9, 3.5773415137059135`*^9}, {
   3.5773416513449545`*^9, 3.577341669752987*^9}, {3.5773417239318824`*^9, 
   3.577341893956581*^9}, 3.5773419407566633`*^9, {3.5773420175555983`*^9, 
   3.577342056165666*^9}, {3.5773420929973307`*^9, 3.5773421056489525`*^9}, {
   3.5773421396414127`*^9, 3.5773423470593767`*^9}, {3.577342382081438*^9, 
   3.577342441158742*^9}, {3.5773425028880506`*^9, 3.577342523682887*^9}, {
   3.5773425980794177`*^9, 3.5773426626791315`*^9}, {3.5773426933487854`*^9, 
   3.5773426941131864`*^9}, {3.577342738354864*^9, 3.577342889035529*^9}, {
   3.5773429412644205`*^9, 3.57734295234044*^9}, {3.577342991824109*^9, 
   3.5773430231333647`*^9}, {3.5773430561274223`*^9, 3.577343143737176*^9}, {
   3.5773431993668737`*^9, 3.577343258334977*^9}, 3.5773434426249013`*^9, {
   3.577343480042567*^9, 3.5773435125686235`*^9}, {3.5773436063903885`*^9, 
   3.5773436272008247`*^9}, {3.5773436617392855`*^9, 
   3.5773437808610945`*^9}, {3.577343868193248*^9, 3.577343879659268*^9}, {
   3.5773439214829416`*^9, 3.577344086609232*^9}, {3.5773444013153844`*^9, 
   3.5773444098953996`*^9}, {3.5773527214601054`*^9, 3.577352754235763*^9}, {
   3.577352839084312*^9, 3.577352839115512*^9}, {3.5773562861299667`*^9, 
   3.577356291496376*^9}, {3.5773717712106485`*^9, 3.577371904793683*^9}, {
   3.577372601444106*^9, 3.577372628525754*^9}, 3.577375878432662*^9}],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
   RowBox[{"projection", " ", "on", " ", "xy"}], "-", 
   RowBox[{"plane", " ", "of", " ", "molecule"}]}], " ", "*)"}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5773413317003937`*^9, 3.5773413362088013`*^9}, {
   3.577352839115512*^9, 3.577352839131112*^9}, {3.5773715990175457`*^9, 
   3.5773716034635534`*^9}, {3.5773716711988726`*^9, 
   3.5773716785932856`*^9}, {3.5773719148557005`*^9, 
   3.5773720320119066`*^9}, {3.5773720625879602`*^9, 3.577372089326407*^9}, {
   3.5773721384508934`*^9, 3.577372247229884*^9}, {3.5773723618276854`*^9, 
   3.577372407098965*^9}, {3.577372640428575*^9, 3.57737267169103*^9}, {
   3.5773727054338894`*^9, 3.5773727803452206`*^9}, {3.577372855490553*^9, 
   3.5773728739765854`*^9}, {3.5773729151762576`*^9, 3.577372935331493*^9}, 
   3.577373068493327*^9}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"plot1", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics", "[", 
       RowBox[{"{", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "i", "]"}], "]"}], ",", 
            RowBox[{"yc", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"rr", "[", 
           RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "}"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "9"}], "}"}]}], "]"}], " ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\n", 
 RowBox[{"Show", "[", 
  RowBox[{"plot1", ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", 
     RowBox[{"Axes", " ", "->", " ", "True"}], ",", 
     RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", 
     RowBox[{"PlotRange", " ", "\[Rule]", " ", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "}"}]}], ",", " ", 
     RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], 
  "]"}], "\n"}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5773413317003937`*^9, 3.5773413362088013`*^9}, {
   3.577352839115512*^9, 3.577352839131112*^9}, {3.5773715990175457`*^9, 
   3.5773716034635534`*^9}, {3.5773716711988726`*^9, 
   3.5773716785932856`*^9}, {3.5773719148557005`*^9, 
   3.5773720320119066`*^9}, {3.5773720625879602`*^9, 3.577372089326407*^9}, {
   3.5773721384508934`*^9, 3.577372247229884*^9}, {3.5773723618276854`*^9, 
   3.577372407098965*^9}, {3.577372640428575*^9, 3.57737267169103*^9}, {
   3.5773727054338894`*^9, 3.5773727803452206`*^9}, {3.577372855490553*^9, 
   3.5773728739765854`*^9}, {3.5773729151762576`*^9, 3.577372935331493*^9}, 
   3.577373068493327*^9, 3.577375863940237*^9}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{1.1206040093493307`*^8, 0}, 0.1791044776119398], 
   CircleBox[{-1.1206040093493307`*^8, 0}, 0.1791044776119398], 
   CircleBox[{1.8931147608628678`*^8, 0}, 7.725107497224924*^7], 
   CircleBox[{-1.8931147608628678`*^8, 0}, 7.725107497224924*^7], 
   CircleBox[{0, 0}, 1.1206040075582859`*^8], 
   CircleBox[{1.1206040093493307`*^8, 7.36327469994904*^7}, 
    7.36327469994904*^7], 
   CircleBox[{1.1206040093493307`*^8, -7.36327469994904*^7}, 
    7.36327469994904*^7], 
   CircleBox[{-1.1206040093493307`*^8, 0}, 7.36327469994904*^7], 
   CircleBox[{-1.1206040093493307`*^8, 0}, 7.36327469994904*^7]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-4, 4}, {-3, 3}}]], "Output",
 CellChangeTimes->{
  3.577376063183787*^9, {3.5773761074722643`*^9, 3.5773761538511457`*^9}, {
   3.5773761950196185`*^9, 3.577376209621244*^9}, 3.5773764153076053`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
   RowBox[{
    RowBox[{"projection", " ", "on", " ", "xz"}], "-", "plane"}], ",", " ", 
   RowBox[{"perpendicular", " ", "to", " ", "molecular", " ", "plane"}]}], 
  " ", "*)"}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5773413317003937`*^9, 3.5773413362088013`*^9}, {
  3.577352839115512*^9, 3.577352839131112*^9}, {3.5773715990175457`*^9, 
  3.5773716034635534`*^9}, {3.5773716711988726`*^9, 3.5773716785932856`*^9}, {
  3.5773719148557005`*^9, 3.5773720320119066`*^9}, {3.5773720625879602`*^9, 
  3.577372089326407*^9}, {3.5773721384508934`*^9, 3.577372247229884*^9}, {
  3.5773723618276854`*^9, 3.577372407098965*^9}, {3.577372640428575*^9, 
  3.57737267169103*^9}, {3.5773727054338894`*^9, 3.5773727803452206`*^9}, {
  3.577372855490553*^9, 3.5773728739765854`*^9}, {3.5773729151762576`*^9, 
  3.577372935331493*^9}, {3.577373068493327*^9, 3.577373077494543*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"plot2", "=", 
   RowBox[{"plot1", "=", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"Graphics", "[", 
        RowBox[{"{", 
         RowBox[{"Circle", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"xc", "[", 
              RowBox[{"[", "i", "]"}], "]"}], ",", 
             RowBox[{"zc", "[", 
              RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
           RowBox[{"rr", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "}"}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "9"}], "}"}]}], "]"}], " ", "/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}]}]}], ";", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{"Graphics", "[", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "1", "]"}], "]"}], ",", 
            RowBox[{"zc", "[", 
             RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"rr", "[", 
           RowBox[{"[", "1", "]"}], "]"}]}], "]"}], ",", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "2", "]"}], "]"}], ",", 
            RowBox[{"zc", "[", 
             RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"rr", "[", 
           RowBox[{"[", "1", "]"}], "]"}]}], "]"}], ",", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "3", "]"}], "]"}], ",", 
            RowBox[{"zc", "[", 
             RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"rr", "[", 
           RowBox[{"[", "2", "]"}], "]"}]}], "]"}], ",", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "4", "]"}], "]"}], ",", 
            RowBox[{"zc", "[", 
             RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"rr", "[", 
           RowBox[{"[", "2", "]"}], "]"}]}], "]"}], ",", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "5", "]"}], "]"}], ",", 
            RowBox[{"zc", "[", 
             RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"rr", "[", 
           RowBox[{"[", "3", "]"}], "]"}]}], "]"}], ",", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "6", "]"}], "]"}], ",", 
            RowBox[{"zc", "[", 
             RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"rr", "[", 
           RowBox[{"[", "4", "]"}], "]"}]}], "]"}], ",", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "7", "]"}], "]"}], ",", 
            RowBox[{"zc", "[", 
             RowBox[{"[", "7", "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"rr", "[", 
           RowBox[{"[", "4", "]"}], "]"}]}], "]"}], ",", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "8", "]"}], "]"}], ",", 
            RowBox[{"zc", "[", 
             RowBox[{"[", "8", "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"rr", "[", 
           RowBox[{"[", "4", "]"}], "]"}]}], "]"}], ",", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "9", "]"}], "]"}], ",", 
            RowBox[{"zc", "[", 
             RowBox[{"[", "9", "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"rr", "[", 
           RowBox[{"[", "4", "]"}], "]"}]}], "]"}], ",", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "10", "]"}], "]"}], ",", 
            RowBox[{"zc", "[", 
             RowBox[{"[", "10", "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"rr", "[", 
           RowBox[{"[", "5", "]"}], "]"}]}], "]"}], ",", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "11", "]"}], "]"}], ",", 
            RowBox[{"zc", "[", 
             RowBox[{"[", "11", "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"rr", "[", 
           RowBox[{"[", "5", "]"}], "]"}]}], "]"}], ",", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "12", "]"}], "]"}], ",", 
            RowBox[{"zc", "[", 
             RowBox[{"[", "12", "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"rr", "[", 
           RowBox[{"[", "5", "]"}], "]"}]}], "]"}], ",", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "13", "]"}], "]"}], ",", 
            RowBox[{"zc", "[", 
             RowBox[{"[", "13", "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"rr", "[", 
           RowBox[{"[", "5", "]"}], "]"}]}], "]"}]}], "}"}], " ", "]"}], "  ",
      "/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}], ";"}], "*)"}], "\n", "\n", 
  RowBox[{"Show", "[", 
   RowBox[{"plot2", ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", 
      RowBox[{"Axes", " ", "->", " ", "True"}], ",", 
      RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", 
      RowBox[{"PlotRange", " ", "\[Rule]", " ", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "}"}]}], ",", " ", 
      RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}]}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5773413317003937`*^9, 3.5773413362088013`*^9}, {
   3.577352839115512*^9, 3.577352839131112*^9}, {3.5773715990175457`*^9, 
   3.5773716034635534`*^9}, {3.5773716711988726`*^9, 
   3.5773716785932856`*^9}, {3.5773719148557005`*^9, 
   3.5773720320119066`*^9}, {3.5773720625879602`*^9, 3.577372089326407*^9}, {
   3.5773721384508934`*^9, 3.577372247229884*^9}, {3.5773723618276854`*^9, 
   3.577372407098965*^9}, {3.577372640428575*^9, 3.57737267169103*^9}, {
   3.5773727054338894`*^9, 3.5773727803452206`*^9}, {3.577372855490553*^9, 
   3.5773728739765854`*^9}, {3.5773729151762576`*^9, 3.577372935331493*^9}, {
   3.577373068493327*^9, 3.577373077494543*^9}, 3.577375870741849*^9}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{1.1206040093493307`*^8, 0}, 0.1791044776119398], 
   CircleBox[{-1.1206040093493307`*^8, 0}, 0.1791044776119398], 
   CircleBox[{1.8931147608628678`*^8, 0}, 7.725107497224924*^7], 
   CircleBox[{-1.8931147608628678`*^8, 0}, 7.725107497224924*^7], 
   CircleBox[{0, 0}, 1.1206040075582859`*^8], 
   CircleBox[{1.1206040093493307`*^8, 0}, 7.36327469994904*^7], 
   CircleBox[{1.1206040093493307`*^8, 0}, 7.36327469994904*^7], 
   CircleBox[{-1.1206040093493307`*^8, 7.36327469994904*^7}, 
    7.36327469994904*^7], 
   CircleBox[{-1.1206040093493307`*^8, -7.36327469994904*^7}, 
    7.36327469994904*^7]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-4, 4}, {-3, 3}}]], "Output",
 CellChangeTimes->{
  3.577376063183787*^9, {3.5773761074878645`*^9, 3.577376153866746*^9}, {
   3.5773761950352182`*^9, 3.5773762096368437`*^9}, 3.577376415323205*^9}]
}, Open  ]]
},
WindowToolbars->"EditBar",
WindowSize->{1009, 964},
WindowMargins->{{Automatic, 422}, {-9, Automatic}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}},
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[557, 20, 4636, 94, 252, "Input"],
Cell[CellGroupData[{
Cell[5218, 118, 4658, 99, 960, "Input"],
Cell[9879, 219, 561, 13, 33, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10477, 237, 2123, 41, 132, "Input"],
Cell[12603, 280, 230, 4, 31, "Output"],
Cell[12836, 286, 213, 3, 31, "Output"],
Cell[13052, 291, 217, 3, 31, "Output"],
Cell[13272, 296, 212, 3, 31, "Output"],
Cell[13487, 301, 216, 3, 31, "Output"]
}, Open  ]],
Cell[13718, 307, 1742, 29, 31, "Input"],
Cell[15463, 338, 859, 15, 31, "Input"],
Cell[CellGroupData[{
Cell[16347, 357, 2050, 50, 92, "Input"],
Cell[18400, 409, 937, 19, 293, "Output"]
}, Open  ]],
Cell[19352, 431, 931, 16, 31, "Input"],
Cell[CellGroupData[{
Cell[20308, 451, 7010, 189, 192, "Input"],
Cell[27321, 642, 935, 19, 293, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)