Water molecule bond length

Classical Mechanics

Newton Equations

$F=ma$ where $F=-\nabla V$. The system is 1D, thus gradient will be differential, and $F=-{\frac {dV}{dx}}$.

Potential Function

Lennard--Jones potential with parameters for TIPS model:

$V(x)=4\epsilon \left[\left({\frac {\sigma }{r}}\right)^{12}-\left({\frac {\sigma }{r}}\right)^{6}\right]$

${\begin{aligned}A&=580.0\times 10^{3}kcalA^{12}/mol\\B&=525.0kcalA^{6}/mol\end{aligned}}$ where $A$ is Ångströms.

Integration

The finite differences (Euler method) are

${\begin{aligned}v(t)&={\frac {x(t+\Delta t)-x(t)}{\Delta t}}\\a(t)&={\frac {v(t+\Delta t)-v(t)}{\Delta t}}\end{aligned}}$ and

$a(t)=-{\frac {{\frac {d}{dx}}(V(x(t)))}{m}}$

Velocity Verlet Algorithm

A very good and easy to implement integration method is velocity Verlet:

${\begin{aligned}x(t+\Delta t)&=x(t)+v(t)\Delta t+{\frac {1}{2}}a\Delta t^{2}\\v(t+\Delta t)&=v(t)+{\frac {1}{2}}\left(a(t)+a(t+\Delta t)\right)\Delta t\end{aligned}}$

Temperature/ Initial distribution

The initial velocity of the hydrogen atom is chosen randomly from the Maxwell-Boltzmann distribution at given temperature $T$

$p(v)=\left({\frac {m}{2\pi k_{B}T}}\right)^{1/2}\exp \left[-{\frac {1}{2}}{\frac {mv^{2}}{k_{B}T}}\right]$

Issues

1D statement