Molecular dynamics quantum rate

Progress in the theoretical description of reaction rates in solution of course correlates strongly with that in other theoretical disciplines, in particular those which have profited most from the enonnous advances in computing power such as quantum chemistry and equilibrium as well as non-equilibrium statistical mechanics of liquid solutions where Monte Carlo and molecular dynamics simulations in many cases have taken on the traditional role of experunents, as they allow the detailed investigation of the influence of intra- and intemiolecular potential parameters on the microscopic dynamics not accessible to measurements in the laboratory. No attempt, however, will be made here to address these areas in more than a cursory way, and the interested reader is referred to the corresponding chapters of the encyclopedia. 

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. 

These results suggest a computational strategy for the study of reactions in condensed phases. One starts from some realistic intermolecular potentials and performs a molecular-dynamics-Kramers-Grote-Hynes scheme that consists of the following steps.First, we fix the proton at the transition state and run a MD simulation. The friction kernel y(t) is calculated and along with Eqs. (7,8) enables the calculation of the Grote-Hynes rate. This scheme has also been used as a means of obtaining input for quantum calculations as well. ... 

Table II Comparison of the ratio k/k E °f quantum rate k over k, which is the TST result corrected for zero-point energy in the reactant well. Also shown are the Landail-Zener and centroid calculations67 and the molecular dynamics with quantum transition result.68...

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. 

The research of Paul Brumer and his colleagues addresses several fundamental problems in theoretical chemical physics. These include studies of the control of molecular dynamics with lasers.98 In particular, the group has demonstrated that quantum interference effects can be used to control the motion of molecules, opening up a vast new area of research. For example, one can alter the rate and yield of production of desirable molecules in chemical reactions, alter the direction of motion of electrons in semiconductors, and change the refractive indices of materials etc. by creating and manipulating quantum interferences. In essence, this approach, called coherent control, provides a method for manipulating chemistry at its most fundamental level.99... 

Abstract We use Nuclear Magnetic Resonance relaxometry (i.e. the frequency variation of the NMR relaxation rates) of quadrupolar nucleus ( Na) and H Pulsed Gradient Spin Echo NMR to determine the mobility of the counterions and the water molecules within aqueous dispersions of clays. The local ordering of isotropic dilute clay dispersions is investigated by NMR relaxometry. In contrast, the NMR spectra of the quadrupolar nucleus and the anisotropy of the water self-diffusion tensor clearly exhibit the occurrence of nematic ordering in dense aqueous dispersions. Multi-scale numerical models exploiting molecular orbital quantum calculations, Grand Canonical Monte Carlo simulations, Molecular and Brownian Dynamics are used to interpret the measured water mobility and the ionic quadrupolar relaxation measurements. 

In the oxygen VER experiments (3) the n = 1 vibrational state of a given oxygen molecule is prepared with a laser, and the population of that state, probed at some later time, decays exponentially. Since in this case tiojo kT, we are in the limit where the state space can be truncated to two levels, and 1/Ti k, 0. Thus the rate constant ki o is measured directly in these experiments. Our starting point for the theoretical discussion is then Equation (14). For reasons discussed in some detail elsewhere (6), for this problem we use the Egelstaff scheme in Equation (19) to relate the Fourier transform of the quantum force-force time-correlation function to the classical time-correlation function, which we then calculate from a classical molecular dynamics computer simulation. The details of the simulation are reported elsewhere (4) here we simply list the site-site potential parameters used therein e/k = 38.003 K, and a = 3.210 A, and the distance between sites is re = 0.7063 A. 

Solvent effects can significantly influence the function and reactivity of organic molecules.1 Because of the complexity and size of the molecular system, it presents a great challenge in theoretical chemistry to accurately calculate the rates for complex reactions in solution. Although continuum solvation models that treat the solvent as a structureless medium with a characteristic dielectric constant have been successfully used for studying solvent effects,2,3 these methods do not provide detailed information on specific intermolecular interactions. An alternative approach is to use statistical mechanical Monte Carlo and molecular dynamics simulation to model solute-solvent interactions explicitly.4 8 In this article, we review a combined quantum mechanical and molecular mechanical (QM/MM) method that couples molecular orbital and valence bond theories, called the MOVB method, to determine the free energy reaction profiles, or potentials of mean force (PMF), for chemical reactions in solution. We apply the combined QM-MOVB/MM method to... 

The phase space structure of classical molecular dynamics is extensively used in developing classical reaction rate theory. If the quanmm reaction dynamics can also be viewed from a phase-space perspective, then a quantum reaction rate theory can use a significant amount of the classical language and the quantum-classical correspondence in reaction rate theory can be closely examined. This is indeed possible by use of, for example, the Wigner function approach. For simplicity let us consider a Hamiltonian system with only one DOF. Generalization to many-dimensional systems is straightforward. The Wigner function associated with a density operator /)( / is defined by... 

Unfortunately, quantitatively reliable quantum chemical calculations of nucleation rates for atmospherically relevant systems would require the application of both high-level electronic structure methods and complicated anharmonic thermochemical analysis to large cluster structures. Such computations are therefore computationally too expensive for currently available computer systems, and will likely remain so for the foreseeable future. Instead, a synthesis of different approaches will probably be necessary. In the future, successful nucleation studies are likely to contain combinations of the best features of both classical (Monte Carlo and molecular dynamics) and quantum chemical methods, with the ultimate objective being a chemically accurate, complete configurational sampling. 

---