Atom dynamics

The simulation of molecular (or atomic) dynamics on a computer was invented by the physicist George Vineyard, working at Brookhaven National Laboratory in New York State. This laboratory, whose biography has recently been published (Crease 1999), was set up soon after World War II by a group of American universities. 

Atom dynamics Group contribution and rigid bonds/angels Specific adsorption Dipolar hard sphere SPC, ST2, TIPS Polarizable H Bonds... 

Section II discusses the real wave packet propagation method we have found useful for the description of several three- and four-atom problems. As with many other wave packet or time-dependent quantum mechanical methods, as well as iterative diagonalization procedures for time-independent problems, repeated actions of a Hamiltonian matrix on a vector represent the major computational bottleneck of the method. Section III discusses relevant issues concerning the efficient numerical representation of the wave packet and the action of the Hamiltonian matrix on a vector in four-atom dynamics problems. Similar considerations apply to problems with fewer or more atoms. Problems involving four or more atoms can be computationally very taxing. Modern (parallel) computer architectures can be exploited to reduce the physical time to solution and Section IV discusses some parallel algorithms we have developed. Section V presents our concluding remarks. 

In tfiis chapter we address first the electrochemical application of the more familiar method of molecular (or atom) dynamics, and later turn to consider Monte Carlo methods, in each case giving a short introduction that should motivate the reader to pursue reading more specific works. Although the present research field is relatively new, the investigations are already too extensive to review in detail in a single chapter. For this reason, we discuss here the more extended research branches in the field and present a few representative examples. The application of simulations applied to nanostructuring problems is discussed in Chapter 36 liquid-liquid interfaces have been addressed by I. Benjamin (1997). 

The main idea behind classical molecular (or atom) dynamics (MD) is fairly simple. To illustrate this for the relatively simple case of an ensemble of atoms, let us consider a system of N particles, each having mass m, with Cartesian coordinates r,. The motion of this system of particles can be described by solving a set of equations of the type... 

The EAM approach appears to provide a formalism within which realistic potentials which describe atomic dynamics can be developed. It should also provide a method for realistically incorporating adsorbates into dynamics simulations. Both of these applications can be considered significant advances, and will help molecular dynamics simulations to continue to contribute to the understanding of technologically important processes. 

As has become clear in previous sections, atomic thermal parameters refined from X-ray or neutron diffraction data contain information on the thermodynamics of a crystal, because they depend on the atom dynamics. However, as diffracted intensities (in kinematic approximation) provide magnitudes of structure factors, but not their phases, so atomic displacement parameters provide the mean amplitudes of atomic motion but not the phase of atomic displacement (i.e., the relative motion of atoms). This means that vibrational frequencies are not directly available from a model where Uij parameters are refined. However, Biirgi demonstrated [111] that such information is in fact available from sets of (7,yS refined on the same molecular crystals at different temperatures. 

J. W. Halley, Studies of the Interdependence of Electronic and Atomic Dynamics and Structure at the Electrode-Electrolyte Interface, Electrochim. Acta. 41 2229 (1996). 

Nuclei provide a large number of spectroscopic probes for the investigation of solid state reaction kinetics. At the same time these probes allow us to look into the atomic dynamics under in-situ conditions. However, the experimental and theoretical methods needed to obtain relevant results in chemical kinetics, and particularly in atomic dynamics, are rather laborious. Due to characteristic hyperfine interactions, nuclear spectroscopies can, in principle, identify atomic particles and furthermore distinguish between different SE s of the same chemical component on different lattice sites. In addition to the analytical aspect of these techniques, nuclear spectroscopy informs about the microscopic motion of the nuclear probes. In Table 16-2 the time windows for the different methods are outlined. 

March. N.H.. and M.P. Tosi Atomic Dynamics in Liquids. Dover Public,id,ms Inc. Mineola. NY. 1991... 

From the perspective of this symposium, analysis of the atomic dynamics and electronic structure of surfaces constitutes an even more exotic topic than surface atomic geometry. In both cases attention has been focused on a small number of model systems, e.g., single crystal transition metal and semiconductor surfaces, using rather specialized experimental facilities. General reviews have appeared for both atomic surface dynamics (21) and spectroscopic measurements of the electronic structure of single-crystal surfaces (, 22). An important emerging trend in the latter area is the use of synchrotron radiation for studying surface electronic structure via photoemission spectroscopy ( 23) Moreover, the use of the very intense synchrotron radiation sources also will enable major improvements in the application of core-level photoemission for surface chemical analysis (13). 

The rearrangement of nuclei in an elementary chemical reaction takes place over a distance of a few angstrom (1 angstrom = 10 10 m) and within a time of about 10-100 femtoseconds (1 femtosecond = 10 15 s a femtosecond is to a second what one second is to 32 million years ), equivalent to atomic speeds of the order of 1 km/s. The challenges in molecular reaction dynamics are (i) to understand and follow in real time the detailed atomic dynamics involved in the elementary processes, (ii) to use this knowledge in the control of these reactions at the microscopic level, e.g., by means of external laser fields, and (iii) to establish the relation between such microscopic processes and macroscopic quantities like the rate constants of the elementary processes. 

N. H. March and M. P. Tosi, Atomic Dynamics in Liquids, Macmillan, London, 1976. 

For other methods of computation, see for instance M. Born, Problems of Atomic Dynamics, Series II, Lecture 4, Massachusetts Institute of Technology, 1926. 

This chapter begins a series of chapters devoted to electronic structure and transport properties. In the present chapter, the foundation for understanding band structures of crystalline solids is laid. The presumption is, of course, that said electronic structures are more appropriately described from the standpoint of an MO (or Bloch)-type approach, rather than the Heitler-London valence-bond approach. This chapter will start with the many-body Schrodinger equation and the independent-electron (Hartree-Fock) approximation. This is followed with Bloch s theorem for wave functions in a periodic potential and an introduction to reciprocal space. Two general approaches are then described for solving the extended electronic structure problem, the free-electron model and the LCAO method, both of which rely on the independent-electron approximation. Finally, the consequences of the independent-electron approximation are examined. Chapter 5 studies the tight-binding method in detail. Chapter 6 focuses on electron and atomic dynamics (i.e. transport properties), and the metal-nonmetal transition is discussed in Chapter 7. 

For the experienced practitioner of atomic physics there appears to be an enigma right at this point. What does nonlinear chaos theory have to do with linear quantum mechanics, so successful in the classification of atomic states and the description of atomic dynamics The answer, interestingly, is the enormous advances in atomic physics itself. Modern day experiments are able to control essentially isolated atoms and molecules to unprecedented precision at very high quantum numbers. Key elements here are the development of atomic beam techniques and the revolutionary effect of lasers. Given the high quantum numbers, Bohr s correspondence principle tells us that atoms are best understood on the basis of classical mechanics. The classical counterpart of most atoms and molecules, however, is chaotic. Hence the importance of understanding chaos in atomic physics. 

A refined picture of the interaction dynamics can be incorporated into models whose aim is to predict the extent of polymer erosion in LEO. For example, the empirical models mentioned in Sec. 3.1 have the possibility to make quantitative predictions of erosion yields, but they have suffered from a lack of information about the initiation steps and the reaction products. Both models have assumed that CO and CO2 carry carbon away from the surface, but quantitative predictions depend on the branching between these two products. The models can be extended to include contributions from all oxidation products that may carry mass away from the surface, as long as their identity and relative yields are known. In addition, the accuracy of the model predictions would be improved by knowledge of the initial trapping probability of incident oxygen atoms. Dynamical data, such as those described in Secs. 3.3 and 3.4, can bolster erosion models and thus expand their acceptance by the community. 

Equation [73] has the same form as the equations of motion for molecules with constrained internal coordinates, and we already know that such equations can be solved effectively using the SHAKE algorithm4 ° Equations [72] and [73] play a key role in the Car-Parrinello method and enable one to run the dynamics for both ionic and electronic degrees of freedom in parallel. With carefully chosen effective mass p and a small time step, the electronic state adjusts itself instanteously to the nuclear configuration (Born-Oppenheimer principle), and, therefore, the atomic dynamics is computed along the system s Born-Oppenheimer surface. Note that there is no need to carry out the costly matrix-diagonalization procedure for performing electronic structure calculations. 

---