Interatomic potentials many-body

Meath W J and Koulis M 1991 On the construction and use of reliable two- and many-body interatomic and intermolecular potentials J. Moi. Struct. (Theochem) 226 1... 

Classical molecular dynamics (MD) implementing predetermined potentials, either empirical or derived from independent electronic structure calculations, has been used extensively to investigate condensed-matter systems. An important aspect in any MD simulation is how to describe or approximate the interatomic interactions. Usually, the potentials that describe these interactions are determined a priori and the full interaction is partitioned into two-, three-, and many-body contributions, long- and short-range terms, etc., for which suitable analytical functional forms are devised. Despite the many successes with classical MD, the requirement to devise fixed potentials results in several serious problems... 

The form of the potential for the system under study was discussed in many publications [28,202,207,208]. Effective pair potentials are widely used in theoretical estimates and numerical calculations. When a many-particle interatomic potential is taken into account, the quantitative description of experimental data improves. For example, the consideration of three-body interactions along with two-particle interactions made it possible to quantitatively describe the stratification curve for interstitial hydrogen in palladium [209]. Let us describe the pair interaction of all the components (hydrogen and metal atoms in the a. and (j phases) by the Lennard Jones potential cpy(ry) = 4 zi [(ff )12- / )6], where Sy and ai are the parameters of the corresponding potentials. All the distances ry, are considered within c.s. of radius r (1 < r < R), where R is the largest radius of the radii of interaction Ry between atoms / and /). 

Periodicity is an important attribute of crystals with significant implications for their properties. Another important property of these systems is the fact that the amplitudes of atomic motions about their equilibrium positions are small enough to allow a harmonic approximation of the interatomic potential. The resulting theory of atomic motion in harmonic crystals constitutes the simplest example for many-body dynamics, which is discussed in this section. 

For larger Au NPs many theoretical calculations have been made using empirical interatomic potentials. A number of different models have been developed to represent the many-body character of bonding in metals, for example, Finnis-Sinclair, Gupta, and glue models. Here, we discuss the embedded atom method (EAM), which has many similarities with the models mentioned above but can be considered as more... 

Pair-additive interactions continued to be used in most materials-related simulations for over 20 years after Vineyard s work despite well-known deficiencies in their ability to model surface and bulk properties of most materials. Quantitative simulation of materials properties was therefore very limited. A breakthrough in materials-related atomistic simulation occurred in the 1980s, however, with the development of several many-body analytic potential energy functions that allow accurate quantitative predictions of structures and dynamics of materials.These methods demonstrated that even relatively simple analytic interatomic potential functions can capture many of the details of chemical bonding, provided the functional form is carefully derived from sound physical principles. 

Classically-based simulations do not work for metals because of the very delocalized bonding characteristic of nearly-free electrons in the substrate. Pair potentials also fail to woik for semiconductors that have strongly directional covalent bonding. For example, a pair potential would predict a close-packed structure for Si instead of the experimentally-observed diamond structure. Only three-body or n-body interatomic PES s can yield the diamond lattice structure, although that structure may contain zero-energy lattice defects [88Terl], Since n-body interatomic potentials contain too many free parameters, in practice potential ejq)ansions usually incorporate terms with no more than three bodies. The two most popular three-body interatomic PES s are the Stilhnger-Weber (SW) and Tersoff potentials. 

The fundamental difficulty with sub-micron modeling is the proliferation of length-scale niches, the ladder becomes a bit of a slide. At the level of 5-500 nanometers reside classical dynamic methods for atoms, which include many body interatomic potentials as well as the embedded atom method. At this length scale, Nafion looks like a crosslinked comb polymer in a polar solvent. 

The first technique used for simulating the behaviors of CNTs was MD method. This method uses realistic force fields (many-body interatomic potential functions) to determination the total energy of a system of particles. Whit the calculation of the total potential energy and force fields of a system, the realistic calculations of the behavior and the properties of a system of atoms and molecules can be acquired. Although the main aspect of both MD and MC simulations methods is based on second Newton s law, MD methods are deterministic approaches, in comparison to the MC methods that are stochastic ones. 

Embedded atom interatomic potentials can be formally derived from DFT [188-194]. This class of interatomic potentials generally works well for bulk metallic systems (i.e., accurately reproduces mechanical and stractural properties), due to its relationship to the band model of electronic structure, which is related to the overlap of local electron densities to form the/ree electron gas of the metallic system. By parameterizing the form of the electron density about an atomic center, one can build a model for the electron overlap in these systems and, hence, the many-body features of alloy systems. The resultant expression for... 

The basis on which interatomic potential methods are built is that the energy of a system can be expressed as a sum over many-body interaction terms, where the number of bodies runs from 1 through to infinity ... 

The problem with empirical many-body potentials lies in the choice of the parameters used to define the interatomic force fields. [89] Since empirical potentials are usually derived by fitting bulk properties, their application to finite systems can lead to incorrect answers. However, successful examples of empirical potential functions which incorporate duster size dependent effects, such as the sp hybridization in Be dusters, [90] have been reported. 

The potential energy of the system is constructed as a sum of individual bond energies. The interactions are truncated using a cutolf function of the interatomic distance ry. The expressions for the repulsive pair potential Vji(rij) and the attractive pair potential F (ry) have been taken from the original Tersoff potential, but the bond order term by modulating the strength of the attractive potential contribution is expressed by a neural network. This many-body term depends on the local environment of the bonds. There is one separate NN for each bond in the system. For each of these bond, each atom bonded either to atom i or j provides an input vector for the NN of the bond ij. As discussed in the previous section, a major... 

The functional form of the Tersoff potential model is motivated by the structural chanistry of covalent systems. The primary basis for this model is that the strength of a bond (i.e. bond order) depends on the local environment (i.e. coordination number). Most predominantly covalent systans have open structures, as an atom with fewer neighbours will form stronger bonds as compared to a close-packed structure where an atom has many neighbours. Tersoff states that the energy is modelled as a sum of pair-like interactions, where, however the coefficient of the attractive term in the pair-like potential (which plays the role of bond order) depends on the local environment, giving a many-body potential. The empirical interatomic potential function for multi-component systems as proposed by Tersoff is as shown below. 

X 10 ppm amagaH. The coefficients a, G and g arise from two-, three- and four-body interactions, respectively. At low densities the shielding constant depends linearly on density, whereas at high pressures many-body collisions become important as well and cause deviation from linearity. The second virial coefficient arises from the Xe-Xe pair interactions (with the potential V(r), r is the interatomic separation) and can be presented... 

Grand ganonical Monte Carlo simulations using realistic interatomic potentials were performed for a significant number of metallic systems, allowing us to draw a number of interesting conclusions. One of the novel features of the work is the exploration of the electrochemical phenomena of UPD and OPD in terms of lattice models that consider the many-body interactions typical of metallic systems. Thus, without the need to assume a particular type of interaction potential between the particles, phase transition phenomena in metallic monolayers could be studied. These studies comprised the formation... 

---