Chemical Sciences: A Manual for CSIR-UGC National Eligibility Test for Lectureship and JRF/Lennard-Jones potential: Difference between revisions
[unreviewed revision] | [unreviewed revision] |
citation cleanup (see WP:WCC) |
|||
Line 1: | Line 1: | ||
{{Cleanup-reason|date=May 2009|math formatting needs fixed}} |
{{Cleanup-reason|date=May 2009|math formatting needs fixed}} |
||
[[Image: |
[[Image:Argon dimer potential and Lennard-Jones.png|thumb|320px|right|Lennard-Jones potential for argon dimer]] |
||
The '''Lennard-Jones potential''' (also referred to as the '''L-J potential''', '''6-12 potential''' or, less commonly, '''12-6 potential''') is a mathematically simple model that describes the interaction between a pair of neutral atoms or molecules. It was proposed in 1924 by [[John Lennard-Jones]].<ref>Lennard-Jones |
The '''Lennard-Jones potential''' (also referred to as the '''L-J potential''', '''6-12 potential''' or, less commonly, '''12-6 potential''') is a mathematically simple model that describes the interaction between a pair of neutral atoms or molecules. It was proposed in 1924 by [[John Lennard-Jones]].<ref>{{Citation |last=Lennard-Jones |first=J. E. |year=1924 |title=On the Determination of Molecular Fields |journal=Proc. R. Soc. Lond. A |volume=106 |issue=738 |pages=463–477 |doi=10.1098/rspa.1924.0082 }}.</ref> |
||
The expression of the L-J potential is |
The expression of the L-J potential is |
||
Line 13: | Line 13: | ||
The lowest energy arrangement of an infinite number of atoms described by a Lennard-Jones potential is a hexagonal [[close-packing]]. |
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 [[Thermodynamic free energy|free energy]] arrangement becomes cubic [[close packing]] and then liquid. Under pressure the lowest energy structure switches between cubic and hexagonal close packing.<ref>Barron |
On raising temperature, the lowest [[Thermodynamic free energy|free energy]] arrangement becomes cubic [[close packing]] and then liquid. Under pressure the lowest energy structure switches between cubic and hexagonal close packing.<ref>{{Citation |last=Barron |first=T. H. K. |last2=Domb |first2=C. |title=On the Cubic and Hexagonal Close-Packed Lattices |journal=Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences |year=1955 |volume=227 |issue=1171 |pages=447–465 |doi=10.1098/rspa.1955.0023 }}.</ref> |
||
Other more recent methods, such as the [[Stockmayer equation]] and the so-called [[multi equation]], 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 [[Algorithmic efficiency|computational cost]]. |
Other more recent methods, such as the [[Stockmayer equation]] and the so-called [[multi equation]], 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 [[Algorithmic efficiency|computational cost]]. |
||
Line 22: | Line 22: | ||
:<math>V(r) = \epsilon \left[ \left(\frac{r_\mathrm{min}}{r}\right)^{12} - 2\left(\frac{r_\mathrm{min}}{r}\right)^{6} \right],</math> |
:<math>V(r) = \epsilon \left[ \left(\frac{r_\mathrm{min}}{r}\right)^{12} - 2\left(\frac{r_\mathrm{min}}{r}\right)^{6} \right],</math> |
||
where ''r''<sub>min</sub> = {{radic|2|6}}''σ'' is the distance at the minimum of the potential. At ''r''<sub>min</sub>, the potential function has the value -ε. |
where ''r''<sub>min</sub> = {{radic|2|6}}''σ'' is the distance at the minimum of the potential. At ''r''<sub>min</sub>, the potential function has the value -ε. |
||
The simplest formulation, often used internally by simulation software, is: |
The simplest formulation, often used internally by simulation software, is: |
||
Line 129: | Line 128: | ||
|} |
|} |
||
It can be easily verified that ''V''<sub>comp</sub>(''r''<sub>c</sub>) = 0, |
It can be easily verified that ''V''<sub>comp</sub>(''r''<sub>c</sub>) = 0, |
||
thus eliminating the jump discontinuity at ''r'' = ''r''<sub>c</sub>. It should be noted that, although the value of the (unshifted) Lennard Jones potential at ''r'' = ''r''<sub>c</sub> = 2.5''σ'' is rather small, the effect of the truncation can be significant, for instance on the gas–liquid critical point<ref>Smit |
thus eliminating the jump discontinuity at ''r'' = ''r''<sub>c</sub>. It should be noted that, although the value of the (unshifted) Lennard Jones potential at ''r'' = ''r''<sub>c</sub> = 2.5''σ'' is rather small, the effect of the truncation can be significant, for instance on the gas–liquid critical point.<ref>{{Citation |last=Smit |first=B. |title=Phase diagrams of Lennard-Jones fluids |journal=Journal of Chemical Physics |volume=96 |issue=11 |pages=8639 |year=1992 |doi=10.1063/1.462271 }}.</ref> Fortunately, the potential energy can be corrected for this effect in a mean field manner by adding |
||
so-called tail corrections<ref>Frenkel |
so-called tail corrections.<ref>{{Citation |last=Frenkel |first=D. |lastauthoramp=yes |last2=Smit |first2=B. |title=Understanding Molecular Simulation |edition=Second |publisher=Academic Press |location=San Diego |year=2002 |isbn=0122673514 }}.</ref> |
||
==See also== |
==See also== |
||
Line 141: | Line 140: | ||
==References== |
==References== |
||
{{Reflist}} |
|||
<references/> |
|||
[[Category:Thermodynamics]] |
[[Category:Thermodynamics]] |
Revision as of 23:14, 8 May 2010
The Lennard-Jones potential (also referred to as the L-J potential, 6-12 potential or, less commonly, 12-6 potential) is a mathematically simple model that describes the interaction between a pair of neutral atoms or molecules. It was proposed in 1924 by John Lennard-Jones.[1]
The 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, and r is the distance between the particles. These parameters can be fitted to reproduce experimental data or accurate quantum chemistry calculations. The r−12 term describes Pauli repulsion at short ranges due to overlapping electron orbitals and the r−6 term describes attraction at long ranges (van der Waals force, or dispersion force).
The Lennard-Jones potential is an approximation. The form of the repulsion term has no theoretical justification; the repulsion force should depend exponentially on the distance, but the repulsion term of the L-J formula is more convenient due to the ease and efficiency of computing r12 as the square of r6. Its physical origin is related to the Pauli principle: when the electronic clouds surrounding the atoms start to overlap, the energy of the system increases abruptly. The exponent 12 was chosen exclusively because of ease of computation.
The attractive long-range potential, however, is derived from dispersion interactions. The L-J potential is a relatively good approximation and due to its simplicity is often used to describe the properties of gases, and to model dispersion and overlap interactions in molecular models. It is particularly accurate for noble gas atoms and is a good approximation at long and short distances for neutral atoms and molecules. On the graph, Lennard-Jones potential for argon dimer is shown. Small deviation from the accurate empirical potential due to incorrect short range part of the repulsion term can be seen.
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.[2]
Other more recent methods, such as the Stockmayer equation and the so-called multi equation, 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 computational cost.
Alternative expressions
The Lennard-Jones potential function is also often written as
where rmin = 6√2σ is the distance at the minimum of the potential. At rmin, the potential function has the value -ε.
The simplest formulation, often used internally by simulation software, is:
where A = 4εσ12 and B = 4εσ6; conversely, σ = 6√A/B and ε = B2/(4A).
Molecular dynamics simulation: Truncated potential
To save computational time, the Lennard-Jones (LJ) potential is often truncated at the cut-off distance of rc = 2.5σ, where
-
(1)
i.e., at rc = 2.5σ, the LJ potential V is about 1/60th of its minimum value ε (depth of potential well).
Beyond , the computational potential is set to zero. On the other hand, to avoid a jump discontinuity at , as shown in Eq.(1), the LJ potential is shifted upward a little so that the computational potential would be zero exactly at the cut-off distance .
For clarity, let denote the LJ potential as defined above, i.e.,
(2)
The computational potential is defined as follows [3]
(3)
It can be easily verified that Vcomp(rc) = 0, thus eliminating the jump discontinuity at r = rc. It should be noted that, although the value of the (unshifted) 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.[4] Fortunately, the potential energy can be corrected for this effect in a mean field manner by adding so-called tail corrections.[5]
See also
- Molecular mechanics
- Embedded atom model
- Morse potential
- Force field (chemistry)
- Force field implementation
- Lennard-Jones model on SklogWiki.
References
- ↑ Lennard-Jones, J. E. (1924), "On the Determination of Molecular Fields", Proc. R. Soc. Lond. A, 106 (738): 463–477, doi:10.1098/rspa.1924.0082.
- ↑ Barron, T. H. K.; Domb, C. (1955), "On the Cubic and Hexagonal Close-Packed Lattices", Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences, 227 (1171): 447–465, doi:10.1098/rspa.1955.0023.
- ↑ softmatter:Lennard-Jones Potential, Soft matter, Materials Digital Library Pathway
- ↑ Smit, B. (1992), "Phase diagrams of Lennard-Jones fluids", Journal of Chemical Physics, 96 (11): 8639, doi:10.1063/1.462271.
- ↑ Frenkel, D.; Smit, B. (2002), Understanding Molecular Simulation (Second ed.), San Diego: Academic Press, ISBN 0122673514
{{citation}}
: Unknown parameter|lastauthoramp=
ignored (|name-list-style=
suggested) (help).