Molecular dynamics with «-body potentials

In principle, we could find the minimum-energy crystal lattice from electronic structure calculations, determine the appropriate A-body interaction potential in the presence of lattice defects, and use molecular dynamics methods to calculate ab initio dynamic macroscale material properties. Some of the problems associated with this approach are considered by Wallace [1]. Because of these problems it is useful to establish a bridge between the micro-... 

A full-scale treatment of crystal growth, however, requires methods adapted for larger scales on top of these quantum-mechanical methods, such as effective potential methods like the embedded atom method (EAM) [11] or Stillinger-Weber potentials [10] with three-body forces necessary. The potentials are obtained from quantum mechanical calculations and then used in Monte Carlo or molecular dynamics methods, to be discussed below. 

Dang LX, Chang TM (1997) Molecular dynamics study of water clusters, liquid, and liquid-vapor interface of water with many-body potentials. J Chem Phys 106(19) 8149—8159... 

Medium-range interactions can be defined as those which dominate the dynamics when atoms interact with energies within a few eV of their molecular binding energies. These forces determine a majority of the physical and chemical properties of surface reactions which are of interest, and so their incorporation in computer simulations can be very important. Unfortunately, they are usually many-body in nature, and can require complicated functional forms to be adequately represented. This means that severe approximations are often required when one is interested in performing molecular dynamics simulations. Recently, several potentials have been semi-empirically developed which have proven to be sufficiently simple to be useful in computer simulations while still capturing the essentials of chemical bonding. 

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... 

Abstract. The physical nature of nonadditivity in many-particle systems and the methods of calculations of many-body forces are discussed. The special attention is devoted to the electron correlation contributions to many-body forces and their role in the Be r and Li r cluster formation. The procedure is described for founding a model potential for metal clusters with parameters fitted to ab initio energetic surfaces. The proposed potential comprises two-body, three-body, and four body interation energies each one consisting of exchange and dispersion terms. Such kind of ab initio model potentials can be used in the molecular dynamics simulation studies and in the cinalysis of binding in small metal clusters. 

After the formulation of defect thermodynamics, it is necessary to understand the nature of rate constants and transport coefficients in order to make practical use of irreversible thermodynamics in solid state kinetics. Even the individual jump of a vacancy is a complicated many-body problem involving, in principle, the lattice dynamics of the whole crystal and the coupling with the motion of all other atomic structure elements. Predictions can be made by simulations, but the relevant methods (e.g., molecular dynamics, MD, calculations) can still be applied only in very simple situations. What are the limits of linear transport theory and under what conditions do the (local) rate constants and transport coefficients cease to be functions of state When do they begin to depend not only on local thermodynamic parameters, but on driving forces (potential gradients) as well Various relaxation processes give the answer to these questions and are treated in depth later. 

Two primary aspects to the practical implementation of molecular dynamics are (i) the numerical integration of the equations of motion along with the boundary conditions and any constraints on the system and (ii) the choice of the interatomic potential. For a single-component system, the potential energy can be written as an expansion in terms of -body potentials ... 

The molecular dynamics simulation was performed using the MOTECC suite of programs [54] in the context of a microcanonical statistical ensemble. The system considered is a cube, with periodic boundary conditions, which contains 343 water molecules. The molecular dynamic simulation of water performed at ambient conditions revealed good agreement with experimental measurements. The main contribution to the total potential energy comes from the two-body term, while the many-body polarisation term contribution amounts to 23% of the total potential energy. Some of the properties calculated during the simulation are reported in Table 3. 

Computer simulation of molecular dynamics is concerned with solving numerically the simultaneous equations of motion for a few hundred atoms or molecules that interact via specified potentials. One thus obtains the coordinates and velocities of the ensemble as a function of time that describe the structure and correlations of the sample. If a model of the induced polarizabilities is adopted, the spectral lineshapes can be obtained, often with certain quantum corrections [425,426]. One primary concern is, of course, to account as accurately as possible for the pairwise interactions so that by carefully comparing the calculated with the measured band shapes, new information concerning the effects of irreducible contributions of inter-molecular potential and cluster polarizabilities can be identified eventually. Pioneering work has pointed out significant effects of irreducible long-range forces of the Axilrod-Teller triple-dipole type [10]. Very recently, on the basis of combined computer simulation and experimental CILS studies, claims have been made that irreducible three-body contributions are observable, for example, in dense krypton [221]. 

Because of their importance to nucleation kinetics, there have been a number of attempts to calculate free energies of formation of clusters theoretically. The most important approaches for the current discussion are harmonic models, " Monte Carlo studies, and molecular dynamics calcula-tions. In the harmonic model the cluster is assumed to be composed of constituent atoms with harmonic intermolecular forces. The most recent calculations, which use the harmonic model, have taken the geometries of the clusters to be those determined by the minimum in the two-body additive Lennard-Jones potential surface. The oscillator frequencies have been obtained by diagonalizing the Lennard-Jones force constant matrix. In the harmonic model the translational and rotational modes of the clusters are treated classically, and the vibrational modes are treated quantum mechanically. The harmonic models work best at low temjjeratures where anharmonic-ity effects are least important and the system is dominated by a single structure. 

---