Informatics Educational Institutions & Programs

In theoretical chemistry, the Buckingham potential is a formula proposed by Richard Buckingham which describes the Pauli exclusion principle and van der Waals energy $\Phi _{12}(r)$ for the interaction of two atoms that are not directly bonded as a function of the interatomic distance $r$ . It is a variety of interatomic potentials.

\Phi _{12}(r)=A\exp \left(-Br\right)-{\frac {C}{r^{6}}}

Here, $A$ , $B$ and $C$ are constants. The two terms on the right-hand side constitute a repulsion and an attraction, because their first derivatives with respect to $r$ are negative and positive, respectively.

Buckingham proposed this as a simplification of the Lennard-Jones potential, in a theoretical study of the equation of state for gaseous helium, neon and argon.^[1]

As explained in Buckingham's original paper and, e.g., in section 2.2.5 of Jensen's text,^[2] the repulsion is due to the interpenetration of the closed electron shells. "There is therefore some justification for choosing the repulsive part (of the potential) as an exponential function". The Buckingham potential has been used extensively in simulations of molecular dynamics.

Because the exponential term converges to a constant as $r$ → $0$ , while the $r^{-6}$ term diverges, the Buckingham potential becomes attractive as $r$ becomes small. This may be problematic when dealing with a structure with very short interatomic distances, as any nuclei that cross a certain threshold will become strongly (and unphysically) bound to one another at a distance of zero.^[2]

Modified Buckingham (Exp-Six) potential

The modified Buckingham potential, also called the "exp-six" potential, is used to calculate the interatomic forces for gases based on Chapman and Cowling collision theory.^[3] The potential has the form

$\Phi _{12}(r)={\frac {\epsilon }{1-6/\alpha }}\left[{\frac {6}{\alpha }}\exp \left[\alpha \left(1-{\frac {r}{r_{min}}}\right)\right]-\left({\frac {r_{min}}{r}}\right)^{6}\right]$

where $\Phi _{12}(r)$ is the interatomic potential between atom i and atom j, $\epsilon$ is the minimum potential energy, $\alpha$ is the measurement of the repulsive energy steepness which is the ratio $\sigma /r_{min}$ , $\sigma$ is the value of $r$ where $\Phi _{12}(r)$ is zero, and $r_{min}$ is the value of $r$ which can achieve the minimum interatomic potential $\epsilon$ . This potential function is only valid when $r>r_{max}$ , as the potential will decay towards $-\infty$ as $r\rightarrow 0$ . This is corrected by identifying $r_{max}$ , which is the value of $r$ at which the potential is maximized; when $r\leq {r_{max}}$ , the potential is set to infinity.

Coulomb–Buckingham potential

The Coulomb–Buckingham potential is an extension of the Buckingham potential for application to ionic systems (e.g. ceramic materials). The formula for the interaction is

\Phi _{12}(r)=A\exp \left(-Br\right)-{\frac {C}{r^{6}}}+{\frac {q_{1}q_{2}}{4\pi \varepsilon _{0}r}}

where A, B, and C are suitable constants and the additional term is the electrostatic potential energy.

The above equation may be written in its alternate form as

\Phi (r)=\varepsilon \left\{{\frac {6}{\alpha -6}}\exp \left(\alpha \left[1-{\frac {r}{r_{0}}}\right]\right)-{\frac {\alpha }{\alpha -6}}\left({\frac {r_{0}}{r}}\right)^{6}\right\}+{\frac {q_{1}q_{2}}{4\pi \varepsilon _{0}r}}

where $r_{0}$ is the minimum energy distance, $\alpha$ is a free dimensionless parameter and $\varepsilon$ is the depth of the minimum energy.

Beest Kramer van Santen (BKS) potential

The BKS potential is a force field that may be used to simulate the interatomic potential between Silica glass atoms.^[4] Rather than relying only on experimental data, the BKS potential is derived by combining ab initio quantum chemistry methods on small silica clusters to describe accurate interaction between nearest-neighbors, which is the function of accurate force field. The experimental data is applied to fit larger scale force information beyond nearest neighbors. By combining the microscopic and macroscopic information, the applicability of the BKS potential has been extended to both the silica polymorphs and other tetrahedral network oxides systems systems that have same cluster structure, such as aluminophosphates, carbon and silicon.

The form of this interatomic potential is the usual Buckingham form, with the addition of a Coulomb force term. The formula for the BKS potential is expressed as

\Phi _{12}(r)=\left[A_{12}\exp \left(-B_{12}r_{12}\right)-{\frac {C_{12}}{r_{12}^{6}}}\right]+{\frac {q_{1}q_{2}}{r_{12}}}

where $\Phi _{12}(r)$ is the interatomic potential between atom i and atom j, $q_{1}$ and $q_{2}$ are the charges magnitudes, $r_{12}$ is the distance between atoms, and $A_{ij}$ , $B_{ij}$ and $C_{ij}$ are constant parameters based on the type of atoms.^[5]

The BKS potential parameters for common atoms are shown below:^[5]

BKS parameters
i-j	A_ij(eV)	B_ij(Å⁻¹)	C_ij(eV•Å⁶)
O - O	1388.7730	2.76000	175.0000
O - Si	18,003.757	4.87318	133.5381
Si - Si	0	0	0
Al - O	16,008.5345	4.79667	130.5659
Al - Al	0	0	0
P - O	9,034.2080	5.19098	19.8793
P - P	0	0	0

An updated version of the BKS potential introduced a new repulsive term to prevent atom overlapping.^[6] The modified potential is taken as

$\Phi _{12}(r)=\left[A_{12}\exp \left(-B_{12}r_{12}\right)-{\frac {C_{12}}{r_{12}^{6}}}\right]+{\frac {q_{1}q_{2}}{r_{12}}}+{\frac {D_{12}}{r_{12}^{24}}}$

where the constant parameters $D_{ij}$ were chosen to have the following values for Silica glass:

Parameter value for Silica glass
	Si - Si	Si - O	O - O
D_ij(eV•Å²⁴)	3423200	29	113

References

^ Buckingham, R. A. (1938). "The Classical Equation of State of Gaseous Helium, Neon and Argon". Proceedings of the Royal Society A. 168 (933): 264–283. Bibcode:1938RSPSA.168..264B. doi:10.1098/rspa.1938.0173. JSTOR 97239.
^ ^a ^b F. Jensen, Introduction to Computational Chemistry, 2nd ed., Wiley, 2007,
^ Mason, Edward A. (2004-12-29). "Transport Properties of Gases Obeying a Modified Buckingham (Exp-Six) Potential". The Journal of Chemical Physics. 22 (2): 169–186. doi:10.1063/1.1740026. ISSN 0021-9606.
^ van Beest, B. W. H.; Kramer, G. J.; van Santen, R. A. (1990-04-16). "Force fields for silicas and aluminophosphates based onab initiocalculations". Physical Review Letters. 64 (16): 1955–1958. Bibcode:1990PhRvL..64.1955V. doi:10.1103/physrevlett.64.1955. ISSN 0031-9007. PMID 10041537.
^ ^a ^b Kramer, G. J.; Farragher, N. P.; van Beest, B. W. H.; van Santen, R. A. (1991-02-15). "Interatomic force fields for silicas, aluminophosphates, and zeolites: Derivation based onab initiocalculations". Physical Review B. 43 (6): 5068–5080. Bibcode:1991PhRvB..43.5068K. doi:10.1103/physrevb.43.5068. ISSN 0163-1829. PMID 9997885.
^ Carré, Antoine; Ispas, Simona; Horbach, Jürgen; Kob, Walter (2016-11-01). "Developing empirical potentials from ab initio simulations: The case of amorphous silica". Computational Materials Science. 124: 323–334. doi:10.1016/j.commatsci.2016.07.041. ISSN 0927-0256.

External links

Buckingham potential on SklogWiki

[1] Buckingham, R. A. (1938). "The Classical Equation of State of Gaseous Helium, Neon and Argon". Proceedings of the Royal Society A. 168 (933): 264–283. Bibcode:1938RSPSA.168..264B. doi:10.1098/rspa.1938.0173. JSTOR 97239.

[jensen-2] F. Jensen, Introduction to Computational Chemistry, 2nd ed., Wiley, 2007,

[3] Mason, Edward A. (2004-12-29). "Transport Properties of Gases Obeying a Modified Buckingham (Exp-Six) Potential". The Journal of Chemical Physics. 22 (2): 169–186. doi:10.1063/1.1740026. ISSN 0021-9606.

[4] van Beest, B. W. H.; Kramer, G. J.; van Santen, R. A. (1990-04-16). "Force fields for silicas and aluminophosphates based onab initiocalculations". Physical Review Letters. 64 (16): 1955–1958. Bibcode:1990PhRvL..64.1955V. doi:10.1103/physrevlett.64.1955. ISSN 0031-9007. PMID 10041537.

[:0-5] Kramer, G. J.; Farragher, N. P.; van Beest, B. W. H.; van Santen, R. A. (1991-02-15). "Interatomic force fields for silicas, aluminophosphates, and zeolites: Derivation based onab initiocalculations". Physical Review B. 43 (6): 5068–5080. Bibcode:1991PhRvB..43.5068K. doi:10.1103/physrevb.43.5068. ISSN 0163-1829. PMID 9997885.

[6] Carré, Antoine; Ispas, Simona; Horbach, Jürgen; Kob, Walter (2016-11-01). "Developing empirical potentials from ab initio simulations: The case of amorphous silica". Computational Materials Science. 124: 323–334. doi:10.1016/j.commatsci.2016.07.041. ISSN 0927-0256.

[1]

[2]

[3]

[4]

[5]

[6]