Φ L J ( r ) = 4 ϵ [ ( σ r ) 12 − ( σ r ) 6 ] {\displaystyle \displaystyle \Phi _{LJ}(r)=4\epsilon \left[\left({\frac {\sigma }{r}}\right)^{12}-\left({\frac {\sigma }{r}}\right)^{6}\right]}
V ( r ) = { Φ L J ( r ) − Φ L J ( R c ) , si r ≤ R c 0 if r > R c {\displaystyle V(r)={\begin{cases}\Phi _{LJ}(r)-\Phi _{LJ}(R_{c}),&{\mbox{si }}r\leq R_{c}\\0&{\mbox{if }}r>R_{c}\end{cases}}}
a i = ∂ 2 r i ∂ t 2 = F i m i {\displaystyle a_{i}={\frac {\partial ^{2}r_{i}}{\partial t^{2}}}={\frac {F_{i}}{m_{i}}}}
F i = m i a i {\displaystyle F_{i}=m_{i}a_{i}\;}