Empirical interatomic interactions

Our knowledge of the properties of orbitals indicates that some of the 3d orbitals might be combined with the 45 and 4p orbitals to form bond orbitals in metals, the other 3d orbitals being unsuited to bond formation, but does not suffice to give a theoretical derivation of the number of d orbitals in each of these classes. Empirical evidence, outlined below, indicates that about 2.44 d orbitals (on the average) show only weak interatomic interactions, and that the remaining 2.56 d orbitals combine with the 5 orbital and the p orbitals to form hybrid bond orbitals. 

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

Derissen JL, Smit PH (1978) Intermolecular interactions in crystals of carboxylic acids. IV. Empirical interatomic potential functions. Acta Cryst A 34 842-853... 

Computationally efficient ab initio quantum mechanical calculations within the framework of DFT play a significant role in the study of plasma-surface interactions. First, they are used to parametrize classical force fields for MD simulations. Second, they provide the quantitative accuracy needed in the development of a chemical reaction database for KMC simulations over long time scales upon identification of a surface chemical reaction through MD simulation, DFT can be used to address in quantitative detail the reaction energetics and kinetics. Third, DFT-based chemical reaction analysis and comparison with the corresponding predictions of the empirical interatomic potential used in the MD simulations provides further... 

Several empirical and semiempirical interatomic potentials have been developed for the Si H systembased on extensions and modifications of well-known potentials for Si including up to three-body interactions (StilUnger and Weber, 1985 Biswas and Hamann, 1985 Biswas et al., 1987 Mousseau and Lewis, 1991 Baskes, 1992). Recent atomic-scale simulation work of plasma-surface interactions in the PECVD of Si thin films has been based on an empirical description of interatomic interactions in the Si H system according to Tersoff s (1986, 1988, 1989) potential for Si, as extended by Ohira and co-workers (1994, 1995, 1996) to incorporate Si-H, H-H, and the corresponding three-body interactions. The extension of the potential to include the presence of hydrogen adopted the Tersoff parametrization to fit results of ab initio calculations for the structure and energetics of Sil 1., x <4, species in the gas phase (Ohira et al., 1994,1995,1996). A similar form of... 

In spite of the parametrization effort involved in the construction of an empirical interatomic potential, capturing the complexity of plasma-surface interactions poses severe demands on the accuracy and transfer-ability of such a model of interatomic interactions. Thus, careful testing is required prior to using any empirical model for the study of plasma-surface... 

The classical equations of motion used in molecular mechanics (MM) are only slightly more difficult to solve than simple additive bond energy equations hence MM calculations are fast and not very demanding of computer resources. In molecular mechanics, one determines the structure of a molecule from a knowledge of the force field, a collection of empirical force constants governing, in principle, all classical mechanical interatomic interactions within the molecule. In practice, it is not feasible for a parameter set to include all possible interactions within a complicated molecule. One hopes that all significant interactions have been included in the force field. 

In many applications involving the use of interatomic potentials in physics and materials science, it is not necessary to know the precise form of the force field between the interacting particles. Even if the exact functional form of the potential energy were known, its mathematical complexity would restrict it from being used in simple analytical work, finding utility only in detailed computer calculations. Empirical atomic interactions are based on a simple analytical model, which provides a mathematically tractable, analytical expression for the pairwise interaction between two atoms or ions. 

P, H. Smit, J. 1., Derissen, and F, B. van Duijneveldt, Mol. Phys., 37, 521 (1979). Intermolecular Interactions in Crystals of Carboxylic Adds III. Non-empirical Interatomic Potential Functions. 

Armed with the empirical knowledge that each element in the periodic table has a characteristic spectmm, and that heating materials to a sufficiently high temperature dismpts all interatomic interactions, Bunsen and Kirchoff invented the spectroscope, an instrument that atomizes substances in a flame and then records their emission spectmm. Using this instmment, the elemental composition of several compounds and minerals were deduced by measuring the wavelength of radiation that they emit. In addition, this new science led to the discovery of elements, notably caesium and mbidium. 

A molecular dynamics calculation was performed for thorium mononitride ThN(cr) in the temperature range from 300 to 2800 K to evaluate the thermophysical properties, viz. the lattice parameter, linear thermal expansion coefficient, compressibility, heat capacity (C° ), and thermal conductivity. A Morse-type function added to the Busing-Ida type potential was employed as the potential function for interatomic interactions. The interatomic potential parameters were semi-empirically determined by fitting to the experimental variation of the lattice parameter with temperature. 

Although MD calculations resorting of interatomic potentials have been successful in many instances, the major shortcoming associated with this type of study, is the reliance of the quantitative precision of the predicted property upon the accuracy of the empirical potential used to model interatomic interactions when interatomic distances are substantially different from those used to fit, the model potential. Ab initio MD circumvents entirely this problem and will play a decisive role in the study erf" mantle phases under pressure. In the next section we outline the theory behind a new ab initio constant pressure MD with variable cell shape (VCS). The following section illustrates its use as an efficient structural optimizer for two mantle phases MgSi03-perovskite and C2/c enstatite. Although these were 0 K calculations, finite temperature studies are similarly possible, the current limitation being simply computational power. 

Application of empirical potentials in study of the interatomic interaction. 

The exponent n is an important quantity in the empirical description of the interatomic interactions. This quantity uniquely determines the dependence of energy on distance, E r) Series of the type of (15.18) have been calculated and tabulated (see [7]). For n < 3 these series diverge. As n ooA approaches the number of nearest neighbors, which is 12 for the face-centered crystal lattice. 

It is the method with empirical potential model. As the model, the approximation that any molecule can behave as if obeying Lennard-Jones potential seems to be satisfactory. This (one-center) LJ model is valid only for rare gases and simple spherical molecules. But this model may also be valid for other simple molecules as a zeroth approximation. We may also use two-center LJ model where interatomic interactions are concentrated to flie major two atoms in the molecule. We expect that these one-center and two-center LJ models will play a role of Occam s razor. 

A theoretical calculation was made of the diffusivity of O in crystalline material. This was based upon constrained-path energy minimization and jump-rate theory by using an empirical interatomic potential which had been newly developed for modeling interactions between Si and O atoms. The calculations predicted that an O atom jumped, on (110) planes, from one bond-center site to another. The saddle-point configuration was farther away from the starting configuration than was the midpoint of the path. The O diffusivity was predicted to be given by ... 

Because the accuracy of the predictions of molecular stmc-ture, conformation or dynamic properties of polymers is intimately related to the interatomic interaction potentials, in the past decade, a variety of FFs have been developed for a wide range of synthetic and natural polymers (see, e.g.. References 135-137). In particular, many organic, metaloorganic, and inorganic compounds were explicitly parameterized in liquid or solid phases to empirical and ab initio information within the so-called Class II FFs, thus providing a starting point for approximation of the potential energy landscape of macromolecules. Accurate correlated electronic stmcture methods were used to benchmark particular intermolecular interactions. 

---