# Lennard-Jones Potential

Proposed by Sir John Edward Lennard-Jones, the Lennard-Jones potential describes the potential energy of interaction between two non-bonding atoms or molecules based on their distance of separation. The potential equation accounts for the difference between attractive forces (dipole-dipole, dipole-induced dipole, and London interactions) and repulsive forces.

### Introduction

Imagine two rubber balls separated by a large distance. Both objects are far enough apart that they are not interacting. The two balls can be brought closer together with minimal energy, allowing interaction. The balls can continuously be brought closer together until they are touching. At this point, it becomes difficult to further decrease the distance between the two balls. In order to bring the balls any closer together, increasing amounts of energy must be added. This is because eventually, as the balls begin to invade each other’s space, they repel each other; the force of repulsion is far greater than the force of attraction.

This scenario is similar to that which takes place in neutral atoms and molecules, and is best described by the **Lennard-Jones potential**.

### The Lennard-Jones Potential

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.The Lennard-Jones Potential is given by the following equation:

\[ V(r)= 4 \epsilon \left [ {\left (\dfrac{\sigma}{r} \right )}^{12}-{\left (\dfrac{\sigma}{r} \right )}^{6} \right] \tag{1}\]

or is sometimes expressed as

\[ V(r) = \frac{A}{r^{12}}- \dfrac{B}{r^6} \tag{2}\]

where

**V**is the intermolecular potential between the two atoms or molecules.- \(\epsilon\) is the well depth and a measure of how strongly the two particles attract each other.
- \(\sigma\) is the distance at which the intermolecular potential between the two particles is zero (See Figure 1.2). \(\sigma\) gives a measurement of how close two nonbonding particles can get and is thus referred to as the
*van der Waals radius*. It is equal to one-half of the internuclear distance between nonbonding particles. - \(r\) is the distance of separation between both particles (measured from the center of one particle to the center of the other particle).
- \(A= 4\epsilon \sigma^{12}\), \(B= 4\epsilon \sigma^{6}\)
- Minimum value of \( \Phi_{12}(r) \) at \( r = r_{min} \);

**Equation 2.1** describes both the attraction and repulsion between nonionic particles. The first part of the equation, (\(\sigma\)/r)^{12} describes the repulsive forces between particles while the latter part of the equation, (\(\sigma\)/r)^{6} denotes attraction.

#### Bonding Potential

The Lennard-Jones potential is a function of the distance between the centers of two particles. When two non-bonding particles are an infinite distance apart, the possibility of them coming together and interacting is minimal. For simplicity's sake, their bonding potential energy is considered zero. However, as the distance of separation decreases, the probability of interaction increases. The particles come closer together until they reach a region of separation where the two particles become bound; their bonding potential energy decreases from zero to a negative quantity. While the particles are bound, the distance between their centers continue to decrease until the particles reach an equilibrium, specified by the separation distance at which the minimum potential energy is reached.

If the two bound particles are further pressed together, past their equilibrium distance, repulsion begins to occur: the particles are so close together that their electrons are forced to occupy each other’s orbitals. Repulsion occurs as each particle attempts to retain the space in their respective orbitals. Despite the repulsive force between both particles, their bonding potential energy increases rapidly as the distance of separation decreases.

#### Stability and Force of Interactions

Like the bonding potential energy, the stability of an arrangement of atoms is a function of the Lennard-Jones separation distance. As the separation distance decreases below equilibrium, the potential energy becomes increasingly positive (indicating a repulsive force). Such a large potential energy is energetically unfavorable, as it indicates an overlapping of atomic orbitals. However, at long separation distances, the potential energy is negative and approaches zero as the separation distance increases to infinity (indicating an attractive force). This indicates that at long-range distances, the pair of atoms or molecules experiences a small stabilizing force. Lastly, as the separation between the two particles reaches a distance slightly greater than σ, the potential energy reaches a minimum value (indicating zero force). At this point, the pair of particles is most stable and will remain in that orientation until an external force is exerted upon it.

### Problems

- The \(\epsilon\)
- Calculate the intermolecular potential between two Argon (Ar) atoms separated by a distance of 4.0 Angstroms (use \(\epsilon\)=0.997 kJ/mol and \(\sigma\)=3.40 Angstroms).
- Two molecules, separated by a distance of 3.0 angstroms, are found to have a \(\sigma\) value of 4.10 angstroms. By decreasing the separation distance between both molecules to 2.0 angstroms, the intermolecular potential between the molecules becomes more negative. Do these molecules follow the Lennard-Jones potential? Why or why not?
- The second part of the Lennard-Jones equation is (\(\sigma\)/r)
^{6}and denotes attraction. Name at least three types of intermolecular interactions that represent attraction. - At what separation distance in the Lennard-Jones potential does a species have a repulsive force acting on it? An attractive force? No force?

### Solutions

1. Recall from Section 2.1 that the van der Waals radius is equal to one-half of the internuclear distance between nonbonding particles. Because \(\sigma\) gives a measure of how close two non-bonding particles can be, the van der Waals radius for Xenon (Xe) is given by:

r = \(\sigma\)/2 = 4.10Angstroms/2 = 2.05 Angstroms

2. To solve for the intermolecular potential between the two Argon atoms, we use equation 2.1 where V is the intermolecular potential between two non-bonding particles.

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

The data given are \(\epsilon\)**=0.997 kJ/mol, \(\sigma\)****=3.40 Angstroms**, and the distance of separation, **r=4.0 Angstroms. **We plug these values into equation 2.1 and solve as follows:

\[V = 4(0.997kJ/mol) \left[\left(\dfrac{3.40\;\text{Angstroms}}{4.0\;\text{Angstroms}}\right)^{12}-\left(\dfrac{3.40 \;\text{Angstroms}}{4.0 \;\text{Angstroms}}\right)^6\right]\]

V = 3.988(0.14-0.38)

V = 3.988(-0.24)

V = -0.96 kJ/mol

3. Recall that \(\sigma\) is the distance at which the bonding potential between two particles is zero. On a graph of the Lennard-Jones potential, then, this value gives the x-intersection of the graph. According to the Lennard-Jones potential, any value of \(r\) greater than \(\sigma\) should yield a negative bonding potential and any value of r smaller than \(\sigma\) should yield a positive bonding potential. In this scenario, as the separation between the two molecules decreases from 3.0 angstroms to 2.0 angstroms, the bonding potential is becomes more negative. In essence however, because the starting separation (3.0 angstroms) is already less than \(\sigma\) (4.0 angstroms), decreasing the separation even further (2.0 angstroms) should result in a more positive bonding potential. Therefore, these molecules do not follow the Lennard-Jones potential.

4. From section 2.1, dipole-dipole, dipole-induced dipole, and London interactions are all attractive forces.

5. See Figure C. A species will have a repulsive force acting on it when r is less than the equilibrium distance between the particles. A species will have an attractive force acting on it when r is greater than the equilibrium distance between the particles. Lastly, when \(r\) is equal to the equilibrium distance between both particles, the species will have no force acting upon it.

### References

- Atkins, Peter and de Paula, Julio.
*Physical Chemistry for the Life Sciences*. New York, N.Y. W. H. Freeman Company, 2006. (469-472). - Royer, Donald J.
*Bonding Theory*. New York, NY.:McGraw-Hill, Inc., 1968. - Chang, Raymond.
*Physical Chemistry for the Biosciences.*Sausalito, CA. University Science Books, 2005. (498-500) - ↑ John Edward Lennard-Jones "On the Determination of Molecular Fields. I. From the Variation of the Viscosity of a Gas with Temperature", Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character
**106**pp. 441-462 (1924) § 8 (ii) - ↑ John Edward Lennard-Jones "On the Determination of Molecular Fields. II. From the Equation of State of a Gas", Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character
**106**pp. 463-477 (1924) Eq. 2.05 - ↑ F. London "Zur Theorie und Systematik der Molekularkräfte", Zeitschrift für Physik A Hadrons and Nuclei
**63**pp. 245-279 (1930) - ↑ W. W. Wood and F. R. Parker "Monte Carlo Equation of State of Molecules Interacting with the Lennard‐Jones Potential. I. A Supercritical Isotherm at about Twice the Critical Temperature", Journal of Chemical Physics
**27**pp. 720- (1957) - ↑ L. A. Rowley, D. Nicholson and N. G. Parsonage "Monte Carlo grand canonical ensemble calculation in a gas-liquid transition region for 12-6 Argon", Journal of Computational Physics
**17**pp. 401-414 (1975) - ↑ J. M. Caillol " Critical-point of the Lennard-Jones fluid: A finite-size scaling study", Journal of Chemical Physics
**109**pp. 4885-4893 (1998) - ↑ V. L. Kulinskii "The critical compressibility factor of fluids from the global isomorphism approach", Journal of Chemical Physics
**139**184119 (2013) - ↑ G. A. Vliegenthart and H. N. W. Lekkerkerker "Predicting the gas–liquid critical point from the second virial coefficient", Journal of Chemical Physics
**112**pp. 5364-5369 (2000) - ↑ L. A. Bulavin and V. L. Kulinskii "Generalized principle of corresponding states and the scale invariant mean-field approach", Journal of Chemical Physics '
**133**134101 (2010) - ↑
^{10.0}^{10.1}Ethan A. Mastny and Juan J. de Pablo "Melting line of the Lennard-Jones system, infinite size, and full potential", Journal of Chemical Physics**127**104504 (2007) - ↑ John G. Kirkwood, Victor A. Lewinson, and Berni J. Alder "Radial Distribution Functions and the Equation of State of Fluids Composed of Molecules Interacting According to the Lennard-Jones Potential", Journal of Chemical Physics
**20**pp. 929- (1952) - ↑ Martin A. van der Hoef "Free energy of the Lennard-Jones solid", Journal of Chemical Physics
**113**pp. 8142-8148 (2000) - ↑ J. Karl Johnson, John A. Zollweg and Keith E. Gubbins "The Lennard-Jones equation of state revisited", Molecular Physics
**78**pp. 591-618 (1993) - ↑ E. M. Apfelbaum, V. S. Vorob’ev and G. A. Martynov "Regarding the Theory of the Zeno Line", Journal of Physical Chemistry A
**112**pp. 6042-6044 (2008) - ↑ V. V. Brazhkin, Yu. D. Fomin, A. G. Lyapin, V. N. Ryzhov, and E. N. Tsiok "Widom Line for the Liquid–Gas Transition in Lennard-Jones System", Journal of Physical Chemistry B Article ASAP (2011)
- ↑ Alauddin Ahmed and Richard J. Sadus "Effect of potential truncations and shifts on the solid-liquid phase coexistence of Lennard-Jones fluids", Journal of Chemical Physics
**133**124515 (2010) - ↑ A. Z. Panagiotopoulos "Molecular simulation of phase coexistence: Finite-size effects and determination of critical parameters for two- and three-dimensional Lennard-Jones fluids", International Journal of Thermophysics
**15**pp. 1057-1072 (1994) - ↑ Søren Toxvaerd and Jeppe C. Dyre "Communication: Shifted forces in molecular dynamics", Journal of Chemical Physics
**134**081102 (2011) - ↑ Spotswood D. Stoddard and Joseph Ford "Numerical Experiments on the Stochastic Behavior of a Lennard-Jones Gas System", Physical Review A
**8**pp. 1504-1512 (1973)

### Contributor

- Rabia Naeem