scilib documentation

LeanChemicalTheories / molecular_mechanics.translation_invariance_LJ

Translation invariance of Lennard-Jones potential #

This section defines the translation invariance of Lennard-Jonnes potential function that describes the intermolecular pair potentials in molecular simulations. The commonly used expression is:
$$ E = 4ε [(\frac{σ}{r})^{12} - (\frac{σ}{r})^6] $$ where:

Here we show that if we translate all the coordinates equally - the molecule is not changing conformations and hence the energy does not change. Using pairwise disntace removes the choice of origin and makes the model translational invariant.

Assumption #

The pairwise distance is given by Euclidean distance or L^2 norm where x and y are the coordinates of the two vectors representing molecular positions.

To-Do #

theorem LJ_invariance {n : } (x y t : fin n ) (r rt a b E Et : ) (a1 : r = has_dist.dist x y) (a2 : rt = has_dist.dist (x + t) (y + t)) (a3 : E = a / r ^ 12 - b / r ^ 6) (a4 : Et = a / rt ^ 12 - b / rt ^ 6) :
E = Et