Pair potential

Last updated

In physics, a pair potential is a function that describes the potential energy of two interacting objects solely as a function of the distance between them. [1]

Contents

Some interactions, like Coulomb's law in electrodynamics or Newton's law of universal gravitation in mechanics naturally have this form for simple spherical objects. For other types of more complex interactions or objects it is useful and common to approximate the interaction by a pair potential, for example interatomic potentials in physics and computational chemistry that use approximations like the Lennard-Jones and Morse potentials.


Functional form

The total energy of a system of objects at positions , that interact through pair potential is given by

Equivalently, this can be expressed as

This expression uses the fact that interaction is symmetric between particles and . It also avoids self-interaction by not including the case where .

Potential range

A fundamental property of a pair potential is its range. It is expected that pair potentials go to zero for infinite distance as particles that are too far apart do not interact. In some cases the potential goes quickly to zero and the interaction for particles that are beyond a certain distance can be assumed to be zero, these are said to be short-range potentials. Other potentials, like the Coulomb or gravitational potential, are long range: they go slowly to zero and the contribution of particles at long distances still contributes to the total energy.

Computational cost

The total energy expression for pair potentials is quite simple to use for analytical and computational work. It has some limitations however, as the computational cost is proportional to the square of number of particles. This might be prohibitively expensive when the interaction between large groups of objects needs to be calculated.

For short-range potentials the sum can be restricted only to include particles that are close, reducing the cost to linearly proportional to the number of particles.

Infinitely periodic systems

In some cases it is necessary to calculate the interaction between an infinite number of particles arranged in a periodic pattern.

Beyond pair potentials

Pair potentials are very common in physics and computational chemistry and biology; exceptions are very rare. An example of a potential energy function that is not a pair potential is the three-body Axilrod-Teller potential. Another example is the Stillinger-Weber potential for silicon, which includes the angle in a triangle of silicon atoms as an input parameter. [2] [3]

Common pair potentials

Some commonly used pair potentials are listed below.

Related Research Articles

<span class="mw-page-title-main">Force</span> Influence that can change motion of an object

In physics, a force is an influence that can cause an object to change its velocity, i.e., to accelerate, meaning a change in speed or direction, unless counterbalanced by other forces. The concept of force makes the everyday notion of pushing or pulling mathematically precise. Because the magnitude and direction of a force are both important, force is a vector quantity. The SI unit of force is the newton (N), and force is often represented by the symbol F.

In physics, screening is the damping of electric fields caused by the presence of mobile charge carriers. It is an important part of the behavior of charge-carrying fluids, such as ionized gases, electrolytes, and charge carriers in electronic conductors . In a fluid, with a given permittivity ε, composed of electrically charged constituent particles, each pair of particles interact through the Coulomb force as

<span class="mw-page-title-main">Molecular dynamics</span> Computer simulations to discover and understand chemical properties

Molecular dynamics (MD) is a computer simulation method for analyzing the physical movements of atoms and molecules. The atoms and molecules are allowed to interact for a fixed period of time, giving a view of the dynamic "evolution" of the system. In the most common version, the trajectories of atoms and molecules are determined by numerically solving Newton's equations of motion for a system of interacting particles, where forces between the particles and their potential energies are often calculated using interatomic potentials or molecular mechanical force fields. The method is applied mostly in chemical physics, materials science, and biophysics.

Density-functional theory (DFT) is a computational quantum mechanical modelling method used in physics, chemistry and materials science to investigate the electronic structure of many-body systems, in particular atoms, molecules, and the condensed phases. Using this theory, the properties of a many-electron system can be determined by using functionals, i.e. functions of another function. In the case of DFT, these are functionals of the spatially dependent electron density. DFT is among the most popular and versatile methods available in condensed-matter physics, computational physics, and computational chemistry.

<span class="mw-page-title-main">Lennard-Jones potential</span> Model of intermolecular interactions

In computational chemistry, molecular physics, and physical chemistry the Lennard-Jones potential is an intermolecular pair potential. Out of all the intermolecular potentials, the Lennard-Jones potential is probably the one that has been the most extensively studied. It is considered an archetype model for simple yet realistic intermolecular interactions.

<span class="mw-page-title-main">Electrostatics</span> Study of stationary or slow-moving electric charges

Electrostatics is a branch of physics that studies slow-moving or stationary electric charges.

In particle, atomic and condensed matter physics, a Yukawa potential is a potential named after the Japanese physicist Hideki Yukawa. The potential is of the form:

<span class="mw-page-title-main">Pseudopotential</span>

In physics, a pseudopotential or effective potential is used as an approximation for the simplified description of complex systems. Applications include atomic physics and neutron scattering. The pseudopotential approximation was first introduced by Hans Hellmann in 1934.

<span class="mw-page-title-main">Radial distribution function</span> Description of particle density in statistical mechanics

In statistical mechanics, the radial distribution function, in a system of particles, describes how density varies as a function of distance from a reference particle.

The Bose–Hubbard model gives a description of the physics of interacting spinless bosons on a lattice. It is closely related to the Hubbard model that originated in solid-state physics as an approximate description of superconducting systems and the motion of electrons between the atoms of a crystalline solid. The model was introduced by Gersch and Knollman in 1963 in the context of granular superconductors. The model rose to prominence in the 1980s after it was found to capture the essence of the superfluid-insulator transition in a way that was much more mathematically tractable than fermionic metal-insulator models.

In computational chemistry and computational physics, the embedded atom model, embedded-atom method or EAM, is an approximation describing the energy between atoms and is a type of interatomic potential. The energy is a function of a sum of functions of the separation between an atom and its neighbors. In the original model, by Murray Daw and Mike Baskes, the latter functions represent the electron density. The EAM is related to the second moment approximation to tight binding theory, also known as the Finnis-Sinclair model. These models are particularly appropriate for metallic systems. Embedded-atom methods are widely used in molecular dynamics simulations.

In atomic, molecular, and optical physics and quantum chemistry, the molecular Hamiltonian is the Hamiltonian operator representing the energy of the electrons and nuclei in a molecule. This operator and the associated Schrödinger equation play a central role in computational chemistry and physics for computing properties of molecules and aggregates of molecules, such as thermal conductivity, specific heat, electrical conductivity, optical, and magnetic properties, and reactivity.

<span class="mw-page-title-main">Water model</span> Aspect of computational chemistry

In computational chemistry, a water model is used to simulate and thermodynamically calculate water clusters, liquid water, and aqueous solutions with explicit solvent. The models are determined from quantum mechanics, molecular mechanics, experimental results, and these combinations. To imitate a specific nature of molecules, many types of models have been developed. In general, these can be classified by the following three points; (i) the number of interaction points called site, (ii) whether the model is rigid or flexible, (iii) whether the model includes polarization effects.

First introduced by M. Pollak, the Coulomb gap is a soft gap in the single-particle density of states (DOS) of a system of interacting localized electrons. Due to the long-range Coulomb interactions, the single-particle DOS vanishes at the chemical potential, at low enough temperatures, such that thermal excitations do not wash out the gap.

<span class="mw-page-title-main">Light front quantization</span> Technique in computational quantum field theory

The light-front quantization of quantum field theories provides a useful alternative to ordinary equal-time quantization. In particular, it can lead to a relativistic description of bound systems in terms of quantum-mechanical wave functions. The quantization is based on the choice of light-front coordinates, where plays the role of time and the corresponding spatial coordinate is . Here, is the ordinary time, is one Cartesian coordinate, and is the speed of light. The other two Cartesian coordinates, and , are untouched and often called transverse or perpendicular, denoted by symbols of the type . The choice of the frame of reference where the time and -axis are defined can be left unspecified in an exactly soluble relativistic theory, but in practical calculations some choices may be more suitable than others.

In statistical mechanics, the cluster expansion is a power series expansion of the partition function of a statistical field theory around a model that is a union of non-interacting 0-dimensional field theories. Cluster expansions originated in the work of Mayer & Montroll (1941). Unlike the usual perturbation expansion which usually leads to a divergent asymptotic series, the cluster expansion may converge within a non-trivial region, in particular when the interaction is small and short-ranged.

Static force fields are fields, such as a simple electric, magnetic or gravitational fields, that exist without excitations. The most common approximation method that physicists use for scattering calculations can be interpreted as static forces arising from the interactions between two bodies mediated by virtual particles, particles that exist for only a short time determined by the uncertainty principle. The virtual particles, also known as force carriers, are bosons, with different bosons associated with each force.

<span class="mw-page-title-main">Metadynamics</span> Scientific computer simulation method

Metadynamics is a computer simulation method in computational physics, chemistry and biology. It is used to estimate the free energy and other state functions of a system, where ergodicity is hindered by the form of the system's energy landscape. It was first suggested by Alessandro Laio and Michele Parrinello in 2002 and is usually applied within molecular dynamics simulations. MTD closely resembles a number of newer methods such as adaptively biased molecular dynamics, adaptive reaction coordinate forces and local elevation umbrella sampling. More recently, both the original and well-tempered metadynamics were derived in the context of importance sampling and shown to be a special case of the adaptive biasing potential setting. MTD is related to the Wang–Landau sampling.

<span class="mw-page-title-main">Interatomic potential</span> Functions for calculating potential energy

Interatomic potentials are mathematical functions to calculate the potential energy of a system of atoms with given positions in space. Interatomic potentials are widely used as the physical basis of molecular mechanics and molecular dynamics simulations in computational chemistry, computational physics and computational materials science to explain and predict materials properties. Examples of quantitative properties and qualitative phenomena that are explored with interatomic potentials include lattice parameters, surface energies, interfacial energies, adsorption, cohesion, thermal expansion, and elastic and plastic material behavior, as well as chemical reactions.

The KTHNY-theory describes the melting of crystals in two dimensions (2D). The name is derived from the initials of the surnames of John Michael Kosterlitz, David J. Thouless, Bertrand Halperin, David R. Nelson, and A. Peter Young, who developed the theory in the 1970s. It is, beside the Ising model in 2D and the XY model in 2D, one of the few theories, which can be solved analytically and which predicts a phase transition at a temperature .

References

  1. Pei, Jun; Song, Lin Frank; Merz Jr., Kenneth M. (June 19, 2020). "Pair Potentials as Machine Learning Features". J. Chem. Theory Comput. 16 (8): 5385–5400. doi:10.1021/acs.jctc.9b01246. PMID   32559380. S2CID   219947826 . Retrieved 26 July 2022.
  2. Stillinger, Frank H.; Weber, Thomas A. (15 April 1985). "Computer simulation of local order in condensed phases of silicon". Physical Review B. 31 (8): 5262–5271. Bibcode:1985PhRvB..31.5262S. doi:10.1103/PhysRevB.31.5262. PMID   9936488 . Retrieved 26 July 2022.
  3. Stillinger, Frank H.; Weber, Thomas A. (15 January 1986). "Erratum: Computer simulation of local order in condensed phases of silicon [Phys. Rev. B 31, 5262 (1985)]". Physical Review B. 33 (2): 1451. Bibcode:1986PhRvB..33.1451S. doi: 10.1103/PhysRevB.33.1451 . PMID   9938428.