Buckingham 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

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

Example Coulomb–Buckingham potential curve.

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

${\displaystyle \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

${\displaystyle \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 ${\displaystyle r_{0}}$ is the minimum energy distance, ${\displaystyle \alpha }$ is a free dimensionless parameter and ${\displaystyle \varepsilon }$ is the depth of the minimum energy.

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

${\displaystyle \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 ${\displaystyle \Phi _{12}(r)}$ is the interatomic potential between atom i and atom j, ${\displaystyle q_{1}}$and ${\displaystyle q_{2}}$ are the charges magnitudes, ${\displaystyle r_{12}}$ is the distance between atoms, and ${\displaystyle A_{ij}}$,${\displaystyle B_{ij}}$ and ${\displaystyle C_{ij}}$ are constant parameters based on the type of atoms.[5]

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

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

${\displaystyle \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 ${\displaystyle D_{ij}}$ were chosen to have the following values for Silica glass:

• Buckingham potential on SklogWiki