Lennard-Jones potential


The Lennard-Jones potential is a mathematically simple model that approximates the intermolecular potential energy between a pair of neutral atoms or molecules as metals or cyclic alkanes. This interatomic potential emerged from Max Born's treatise on lattice enthalpies and was elaborated in 1924 by John Lennard-Jones. A common expression of the L-J potential is
where ε is the depth of the potential well, σ is the finite distance at which the inter-particle potential is zero, r is the distance between the particles, and rm is the distance at which the potential reaches its minimum. At rm, the potential function has the value −ε. The distances are related as rm = 21/6σ ≈ 1.122σ. These parameters can be fitted to reproduce experimental data on lattice parameters, surface energies, or accurate quantum chemistry calculations. Due to its computational simplicity and interpretability, the Lennard-Jones potential is used extensively in computer simulations. It is also part of more detailed interatomic potentials such as CHARMM, AMBER, OPLS-AA, CVFF, DREIDING, and IFF.
The Lennard-Jones potential is also used with different exponents, for example, in 9-6 form:
in force fields such as CFF, PCFF, and IFF. The interpretation of rm as atomic diameter and ε as polarizability and a measure for internal cohesion is the same as for the 12-6 Lennard-Jones potential.

Explanation

The r12 term, which is the repulsive term, describes Pauli repulsion at short ranges due to overlapping electron orbitals, and the r6 term, which is the attractive long-range term, describes attraction at long ranges.
Differentiating the L-J potential with respect to r gives an expression for the net inter-molecular force between 2 molecules. This inter-molecular force may be attractive or repulsive, depending on the value of r. When r is very small, the molecules repel each other.
Whereas the functional form of the attractive term has a clear physical justification, the repulsive term has no theoretical justification. It is used because it approximates the Pauli repulsion well and is more convenient due to the relative computing efficiency of calculating r12 as the square of r6.
The L-J potential is a relatively good approximation. Due to its simplicity, it is often used to describe the properties of gases and to model dispersion and overlap interactions in molecular models. It is especially accurate for noble gas atoms and is a good approximation at long and short distances for neutral atoms and molecules.
The lowest-energy arrangement of an infinite number of atoms described by a Lennard-Jones potential is a hexagonal close-packing. On raising temperature, the lowest-free-energy arrangement becomes cubic close packing, and then liquid. Under pressure, the lowest-energy structure switches between cubic and hexagonal close packing. Real materials include body-centered cubic structures also.
The Lennard-Jones potential was improved by the Buckingham potential later proposed by Richard Buckingham, incorporating an extra parameter and the repulsive part is replaced by an exponential function:
Other more recent methods, such as the Stockmayer potential, describe the interaction of molecules more accurately. Quantum chemistry methods, Møller–Plesset perturbation theory, coupled cluster method, or full configuration interaction can give extremely accurate results, but require large computing cost.

Alternative expressions

There are many different ways to formulate the Lennard-Jones potential. Some common forms follow.

AB form

This form is a simplified formulation that is used by some simulation software:
where, and. Conversely, and. This is the form in which Lennard-Jones wrote the 12-6 potential.
A more mathematically general form, which contains an extra variable n is
where is the bonding energy of the molecule. The exponent n could be related to the spring constant k as
from where n can be calculated if k is known. Normally the harmonic states are known,, where. n can also be related to the group velocity in a crystal,
where a is the lattice distance, and m is the mass of an atom.

Truncated and shifted form

To save computing time and satisfy the when using periodic boundary conditions, the Lennard-Jones potential is often truncated at a cut-off distance of rc = 2.5σ, where
i.e., at rc = 2.5σ, the Lennard-Jones potential VLJ is about 1/60 of its minimum value, ε. Beyond, the truncated potential is set to zero.
To avoid a jump discontinuity at, the LJ potential must be shifted upward a little, so that the truncated potential would be zero exactly at the cut-off distance,.
For clarity, let denote the LJ potential as defined above, i.e.,
Then the truncated Lennard-Jones potential
is defined as follows
It can be easily verified that VLJtrunc = 0,
thus eliminating the jump discontinuity at r = rc. Although the value of the Lennard Jones potential at r = rc = 2.5σ is rather small, the effect of the truncation can be significant; for instance, on the gas–liquid critical point. The potential energy can often be corrected for this effect in a mean-field manner by adding so-called tail corrections.
The stability of crystal structure is extremely sensitive to the truncation: a wide variety of close packed stackings may have the lowest energy.

Extensions for polarizability

Lennard-Jones potentials have been combined with bonded potentials to describe virtual electrons and image charges in metals. The combined potential reproduces lattice parameters, surface energies, the complete image potentials, hydration energies, and mechanical properties without additional fit parameters in higher accuracy than density functional theory calculations.

Dimensionless (reduced) units

Dimensionless units can be defined based on the Lennard-Jones potential, which are convenient for molecular dynamics simulations. From a numerical point, the advantages of dimensionless units include computing values which are closer to unity, using simplified equations and being able to easily scale the results.
When the Lennard-Jones potential is used for molecular dynamics simulations, the most convenient dimensionless units are obtained by choosing length σ, mass m and energy ε as the scaling factors for the various physical properties.

Limitations

For simulations of charged particles, e.g. ions, an additional term has been proposed to account for dipole interaction, giving rise to a so-called 12-6-4 potential.
where C is a system-dependent constant.