Rate theories molecular dynamics

The dynamic motion due to rapid energy exchange for the desorption of Xe atoms from a Pd(lOO) surface will be illustrated. Figure 5.9a shows the rate of Xe desorption as predicted according to transition-state theory. Figure 5.9b compares computed molecular-dynamics rates and the transition-state rates. The open data points are the computed desorption rates for Xe atoms that are allowed to readsorb once they have passed the transition-state barrier. The filled data points ignore the possibility of readsorption. The open data points, computed from the more exact theory, always remain lower than the transition-state result. Transition-state theory and molecular dynamics predict very similar rate constants for the desorption of xenon from palladium. 

The Langmuir-Hinshelwood picture is essentially that of Fig. XVIII-14. If the process is unimolecular, the species meanders around on the surface until it receives the activation energy to go over to product(s), which then desorb. If the process is bimolecular, two species diffuse around until a reactive encounter occurs. The reaction will be diffusion controlled if it occurs on every encounter (see Ref. 211) the theory of surface diffusional encounters has been treated (see Ref. 212) the subject may also be approached by means of Monte Carlo/molecular dynamics techniques [213]. In the case of activated bimolecular reactions, however, there will in general be many encounters before the reactive one, and the rate law for the surface reaction is generally written by analogy to the mass action law for solutions. That is, for a bimolecular process, the rate is taken to be proportional to the product of the two surface concentrations. It is interesting, however, that essentially the same rate law is obtained if the adsorption is strictly localized and species react only if they happen to adsorb on adjacent sites (note Ref. 214). (The apparent rate law, that is, the rate law in terms of gas pressures, depends on the form of the adsorption isotherm, as discussed in the next section.)... 

The assumptions of transition state theory allow for the derivation of a kinetic rate constant from equilibrium properties of the system. That seems almost too good to be true. In fact, it sometimes is [8,18-21]. Violations of the assumptions of TST do occur. In those cases, a more detailed description of the system dynamics is necessary for the accurate estimate of the kinetic rate constant. Keck [22] first demonstrated how molecular dynamics could be combined with transition state theory to evaluate the reaction rate constant (see also Ref. 17). In this section, an attempt is made to explain the essence of these dynamic corrections to TST. 

Evans and Baranyai [51, 52] have explored what they describe as a nonlinear generalization of Prigogine s principle of minimum entropy production. In their theory the rate of (first) entropy production is equated to the rate of phase space compression. Since phase space is incompressible under Hamilton s equations of motion, which all real systems obey, the compression of phase space that occurs in nonequilibrium molecular dynamics (NEMD) simulations is purely an artifact of the non-Hamiltonian equations of motion that arise in implementing the Evans-Hoover thermostat [53, 54]. (See Section VIIIC for a critical discussion of the NEMD method.) While the NEMD method is a valid simulation approach in the linear regime, the phase space compression induced by the thermostat awaits physical interpretation even if it does turn out to be related to the rate of first entropy production, then the hurdle posed by Question (3) remains to be surmounted. 

In this chapter, we give a brief account of two related aspects of chemical reactions in solution the so-called "stochastic" theoretical approach to the rates of reactions and related features, and Molecular Dynamics (MD) computer simulations designed to test such theories and to otherwise provide insight on the reaction dynamics. 

We have reviewed above the GH approach to reaction rate constants in solution, together with simple models that give a deeper perspective on the reaction dynamics and various aspects of the generalized frictional influence on the rates. The fact that the theory has always been found to agree with Molecular Dynamics computer simulation results for realistic models of many and varied reaction types gives confidence that it may be used to analyze real experimental results. 

These experimental and numerical developments have posed a challenge to the theorist. Given the complexity of the phenomena involved, is it still possible to present a theory which provides the necessary concepts and insight needed for understanding rate processes in condensed phases Although classical molecular dynamics computations are almost routine, real time quantum molecular dynamics are still largely computationally inaccessible. Are there alternatives Do we understand quantum effects in rate theory These are the topics of this review article. 

Following Fey nman s original work, several authors pmsued extensions of the effective potential idea to construct variational approximations for the quantum partition function (see, e g., Refs. 7,8). The importance of the path centroid variable in quantum activated rate processes was also explored and revealed, which gave rise to path integral quantum transition state theory and even more general approaches. The Centroid Molecular Dynamics (CMD) method for quantum dynamics simulation was also formulated. In the CMD method, the position centroid evolves classically on the efiective centroid potential. Various analysis and numerical tests for realistic systems have shown that CMD captures the main quantum effects for several processes in condensed matter such as transport phenomena. 

Interestingly, in the experiments devoted solely to computational chemistry, molecular dynamics calculations had the highest representation (96-98). The method was used in simulations of simple liquids, (96), in simulations of chemical reactions (97), and in studies of molecular clusters (98). One experiment was devoted to the use of Monte Carlo methods to distinguish between first and second-order kinetic rate laws (99). One experiment used DFT theory to study two isomerization reactions (100). 

Fortunately, the reaction rates of many important processes can be obtained without a full molecular dynamics simulation. Most reaction rate theories for elementary processes build upon the ideas introduced in the so-called transition state theory [88-90]. We shall focus on this theory here, particularly because it (and its harmonic approximation, HTST) has been shown to yield reliable results for elementary processes at surfaces. 

By contrast, few such calculations have as yet been made for diffusional problems. Much more significantly, the experimental observables of rate coefficient or survival (recombination) probability can be measured very much less accurately than can energy levels. A detailed comparison of experimental observations and theoretical predictions must be restricted by the experimental accuracy attainable. This very limitation probably explains why no unambiguous experimental assignment of a many-body effect has yet been made in the field of reaction kinetics in solution, even over picosecond timescale. Necessarily, there are good reasons to anticipate their occurrence. At this stage, all that can be done is to estimate the importance of such effects and include them in an analysis of experimental results. Perhaps a comparison of theoretical calculations and Monte Carlo or molecular dynamics simulations would be the best that could be hoped for at this moment (rather like, though less satisfactory than, the current position in the development of statistical mechanical theories of liquids). Nevertheless, there remains a clear need for careful experiments, which may reveal such effects as discussed in the remainder of much of this volume. 

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. 

Our calculations of the activation free energy barrier for the cuprous-cupric electron transfer were not precise enough to permit a very accurate estimate of the absolute value of the exchange current for comparison with experiment. In principle, a determination of the absolute rate from the activation energy requires a calculation of the relevant correlation function [82] when the ion is in the transition region within the molecular dynamics model. We have not carried out such a calculation, but can obtain some information about the amplitude by comparing experiments with the transition state theory expression [84]... 

Molecular dynamic simulations recently made by Soddemann et al. [52, 53] offer a very precise insight into the behavior of the layered systems under shear. As we will briefly discuss now, a direct comparison of these simulations to the analytic theory presented above shows very good agreement between both approaches [53], In Fig. 13 we have plotted the strain rate, y, as a function of the tilt angle, 0o-... 

Transition state theory, especially with its recent developments, has proved a very powerful tool, vastly superior to collision theory. It has only recently been challenged by modem advances in molecular beams and molecular dynamics which look at the microscopic details of a collision, and which can be regarded as a modified collision theory. These developments along with computer techniques, and modem experimental advances in spectroscopy and lasers along with fast reaction techniques, are now revolutionizing the science of reaction rates. 

Ray Kapral came to Toronto from the United States in 1969. His research interests center on theories of rate processes both in systems close to equilibrium, where the goal is the development of a microscopic theory of condensed phase reaction rates,89 and in systems far from chemical equilibrium, where descriptions of the complex spatial and temporal reactive dynamics that these systems exhibit have been developed.90 He and his collaborators have carried out research on the dynamics of phase transitions and critical phenomena, the dynamics of colloidal suspensions, the kinetic theory of chemical reactions in liquids, nonequilibrium statistical mechanics of liquids and mode coupling theory, mechanisms for the onset of chaos in nonlinear dynamical systems, the stochastic theory of chemical rate processes, studies of pattern formation in chemically reacting systems, and the development of molecular dynamics simulation methods for activated chemical rate processes. His recent research activities center on the theory of quantum and classical rate processes in the condensed phase91 and in clusters, and studies of chemical waves and patterns in reacting systems at both the macroscopic and mesoscopic levels. 

---