Lennard-Jones model

The Lennard-Jones intermolecular pair potential is a special case of the Mie potential and takes its name from Sir John Edward Lennard-Jones The Lennard-Jones model consists of two 'parts'; a steep repulsive term, and smoother attractive term, representing the London dispersion forces. Apart from being an important model in itself, the Lennard-Jones potential frequently forms one of 'building blocks' of many force fields. It is worth mentioning that the 12-6 Lennard-Jones model is not the most faithful representation of the potential energy surface, but rather its use is widespread due to its computational expediency. For example, the repulsive term is maybe better described with the exp-6 potential. One of the first computer simulations using the Lennard-Jones model was undertaken by Rahman in 1964 in a study of liquid argon.

Functional form

The Lennard-Jones potential is given by

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Phi_{12}(r) = 4 \epsilon \left[ \left(\frac{\sigma}{r} \right)^{12}- \left( \frac{\sigma}{r}\right)^6 \right] }

where

• Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r := |\mathbf{r}_1 - \mathbf{r}_2|}
• Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Phi_{12}(r) } is the intermolecular pair potential between two particles or sites
• Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sigma } is the value of Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r} at which Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Phi_{12}(r)=0}
• Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \epsilon } is the well depth (energy)

In reduced units:

• Density: Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho^* := \rho \sigma^3 } , where Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho := N/V } (number of particles Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle N } divided by the volume Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle V } )
• Temperature: Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle T^* := k_B T/\epsilon } , where Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle T } is the absolute temperature and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle k_B } is the Boltzmann constant

The following is a plot of the Lennard-Jones model for the Rowley, Nicholson and Parsonage parameter set (Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \epsilon/k_B = } 119.8 K and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sigma=} 0.3405 nm). See argon for other parameter sets.

Special points

• Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Phi_{12}(\sigma) = 0 }
• Minimum value of $\Phi _{12}(r)$ at Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r = r_{min} } ;

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{r_{min}}{\sigma} = 2^{1/6} \simeq 1.12246 ... }

Critical point

The location of the critical point is

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle T_c^* = 1.326 \pm 0.002}

at a reduced density of

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho_c^* = 0.316 \pm 0.002} .

Vliegenthart and Lekkerkerker have suggested that the critical point is related to the second virial coefficient via the expression

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle B_2 \vert_{T=T_c}= -\pi \sigma^3}

Triple point

The location of the triple point as found by Mastny and de Pablo is

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle T_{tp}^* = 0.694}

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho_{tp}^* = 0.84} (liquid); Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho_{tp}^* = 0.96} (solid)

Approximations in simulation: truncation and shifting

The Lennard-Jones model is often used with a cutoff radius of Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle 2.5 \sigma} , beyond which Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Phi_{12}(r)} is set to zero. Setting the well depth Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \epsilon } to be 1 in the potential on arrives at Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Phi_{12}(r)\simeq -0.0163} , i.e. at this distance the potential is at less than 2% of the well depth. See Mastny and de Pablo for an analysis of the effect of this cutoff on the melting line. See Panagiotopoulos for critical parameters .

n-m Lennard-Jones potential

It is relatively common to encounter potential functions given by:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Phi_{12}(r) = c_{n,m} \epsilon \left[ \left( \frac{ \sigma }{r } \right)^n - \left( \frac{\sigma}{r} \right)^m \right]. }

with Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle n } and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle m } being positive integers and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle n > m } . Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle c_{n,m} } is chosen such that the minimum value of Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Phi_{12}(r) } being Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Phi_{min} = - \epsilon } . Such forms are usually referred to as n-m Lennard-Jones Potential. For example, the 9-3 Lennard-Jones interaction potential is often used to model the interaction between the atoms/molecules of a fluid and a continuous solid wall. On the '9-3 Lennard-Jones potential' page a justification of this use is presented. Another example is the n-6 Lennard-Jones potential, where Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle m} is fixed at 6, and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle n} is free to adopt a range of integer values. The potentials form part of the larger class of potentials known as the Mie potential.

See also

• 8-6 Lennard-Jones potential
• 9-3 Lennard-Jones potential
• 9-6 Lennard-Jones potential
• 10-4-3 Lennard-Jones potential
• n-6 Lennard-Jones potential

Radial distribution function

The following plot is of a typical radial distribution function for the monatomic Lennard-Jones liquid (here with Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sigma=3.73 {\mathrm {\AA}}} and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \epsilon=0.294} kcal/mol at a temperature of 111.06K):

Typical radial distribution function for the monatomic Lennard-Jones liquid.

Equation of state

Main article: Lennard-Jones equation of state

Virial coefficients

Main article: Lennard-Jones model: virial coefficients

Phase diagram

Main article: Phase diagram of the Lennard-Jones model

Perturbation theory

The Lennard-Jones model is also used in perturbation theories, for example see: Weeks-Chandler-Anderson perturbation theory.

Mixtures

• Binary Lennard-Jones mixtures
• Multicomponent Lennard-Jones mixtures

Related models

• Kihara potential
• Lennard-Jones model in 1-dimension (rods)
• Lennard-Jones model in 2-dimensions (disks)
• Lennard-Jones model in 4-dimensions
• Lennard-Jones sticks
• Mie potential
• Soft-core Lennard-Jones model
• Soft sphere potential
• Stockmayer potential