Molecular Dynamics

Etotal = Evalence + Enon-bond

Eele(Rij)= S QiQj/e Rij

Edisp(Rij)= Do{[(6/ z -6) e z(1-Rij/Ro) ]-[(z / z -6)(Ro/Rij)6]} Ehbond(Rij)= Do{ 5(Ro/ R)12 - 6 (Ro/Rij)10]}