Quantum dynamics in simple fluids

We use quantum-correction factors to calculate approximately the quantum velocity time-correlation function (TCF) of supercritical Lennard-Jones argon from the classical TCF. We find that for this quite classical system, several different quantum-correction schemes yield essentially identical results for the real and imaginary parts of the quantum TCF, and also agree well with the recent forward-backward semiclassical dynamics (FBSD) results of Wright and Makri [J. Chem. Phys. 119, 1634 (2003)]. We also consider a more quantum-mechanical fluid of lighter atoms (neon) at a lower temperature. In this case different quantum-correction schemes give different results. FBSD calculations show that the harmonic quantum correction factor works the best for this system     

Evaluation of quantum correlation functions from classical data: Anharmonic models

The previously introduced method of evaluating quantum mechanical time correlation functions using as input only classical simulation data is generalized and applied to two anharmonic model systems, as a further test. The quantum correction approach utilizes the relation between a general quantum correlation function and its classical analog. For the tested models, we obtain numerical results of nonlinear correlation functions with comparable accuracy to that of the centroid molecular dynamics method, although the present method is much simpler to implement and not limited to real valued quantum correlation functions.     

Few-body quantum and many-body classical hyperspherical approaches to reactions and to cluster dynamics

The hyperspherical method is a widely used and successful approach for the quantum treatment of elementary chemical processes. It has been mostly applied to three-atomic systems, and current progress is here outlined concerning the basic theoretical framework for the extension to four-body bound state and reactive scattering problems. Although most applications only exploit the advantages of the hyperspherical coordinate systems for the formulation of the few-body problem, the full power of the technique implies representations explicitly involving quantum hyperangular momentum operators as dynamical quantities and hyperspherical harmonics as basis functions. In terms of discrete analogues of these harmonics one has a universal representation for the kinetic energy and a diagonal representation for the potential (hyperquantization algorithm). Very recently, advances have been made on the use of the approach in classical dynamics, provided that a hyperspherical formulation is given based on “classical” definitions of the hyperangular momenta and related quantities. The aim of the present paper is to offer a retrospective and prospective view of the hyperspherical methods both in quantum and classical dynamics. Specifically, regarding the general quantum hyperspherical approaches for three- and four-body systems, we first focus on the basis set issue, and then we present developments on the classical formulation that has led to applications involving the implementations of hyperspherical techniques for classical molecular dynamics simulations of simple nanoaggregates.     

Forward--backward semiclassical dynamics with information-guided noise reduction for a molecule in solution

The forward--backward semiclassical dynamics (FBSD) methodology is used to obtain expressions for time correlation functions of a system (atom or molecule) in solution. We use information-guided noise reduction (IGNoR) [Makri, N. Chem. Phys. Lett. 2004, 400, 446] to minimize the statistical error associated with the Monte Carlo integration of oscillatory functions. This is possible by reformulating the correlation function in terms of an oscillatory solvent-dependent contribution whose integral can be obtained analytically and a slowly varying function obtained via a grid-based iterative evaluation of solute properties. Knowledge of the exact integral of the oscillatory function, combined with correlated statistics, leads to partial cancellation of the Monte Carlo error. Application on a one-dimensional solute-solvent model shows a substantial improvement of convergence in the IGNoR-enhanced FBSD correlation function for a fixed number of Monte Carlo samples. The reduction of statistical error achieved by using the IGNoR methodology becomes more significant as the number of solvent particles increases.     

Semiclassical initial value series representation in the continuum limit: application to vibrational relaxation

A recently formulated continuum limit semiclassical initial value series representation (SCIVR) of the quantum dynamics of dissipative systems is applied to the study of vibrational relaxation of model harmonic and anharmonic oscillator systems. As is well known, the classical dynamics of dissipative systems may be described in terms of a generalized Langevin equation. The continuum limit SCIVR uses the Langevin trajectories as input, albeit with a quantum noise rather than a classical noise. Combining this development with the forward-backward form of the prefactor-free propagator leads to a tractable scheme for computing quantum thermal correlation functions. Here we present the first implementation of this continuum limit SCIVR series method to study two model problems of vibrational relaxation. Simulations of the dissipative harmonic oscillator system over a wide range of parameters demonstrate that at most only the first two terms in the SCIVR series are needed for convergence of the correlation function. The methodology is then applied to the vibrational relaxation of a dissipative Morse oscillator. Here, too, the SCIVR series converges rapidly as the first two terms are sufficient to provide the quantum mechanical relaxation with an estimated accuracy on the order of a few percent. The results in this case are compared with computations obtained using the classical Wigner approximation for the relaxation dynamics.     

Simulating quantum dynamics in classical environments

Methods for simulating the dynamics of composite systems, where part of the system is treated quantum mechanically and its environment is treated classically, are discussed. Such quantum–classical systems arise in many physical contexts where certain degrees of freedom have an essential quantum character while the other degrees of freedom to which they are coupled may be treated classically to a good approximation. The dynamics of these composite systems are governed by a quantum–classical Liouville equation for either the density matrix or the dynamical variables which are operators in the Hilbert space of the quantum subsystem and functions of the classical phase space variables of the classical environment. Solutions of the evolution equations may be formulated in terms of surface-hopping dynamics involving ensembles of trajectory segments interspersed with quantum transitions. The surface-hopping schemes incorporate quantum coherence and account for energy exchanges between the quantum and classical degrees of freedom. Various simulation algorithms are discussed and illustrated with calculations on simple spin-boson models but the methods described here are applicable to realistic many-body environments.     

An accurate and simple quantum model for liquid water

The path-integral molecular dynamics and centroid molecular dynamics methods have been applied to investigate the behavior of liquid water at ambient conditions starting from a recently developed simple point charge/flexible (SPC/Fw) model. Several quantum structural, thermodynamic, and dynamical properties have been computed and compared to the corresponding classical values, as well as to the available experimental data. The path-integral molecular dynamics simulations show that the inclusion of quantum effects results in a less structured liquid with a reduced amount of hydrogen bonding in comparison to its classical analog. The nuclear quantization also leads to a smaller dielectric constant and a larger diffusion coefficient relative to the corresponding classical values. Collective and single molecule time correlation functions show a faster decay than their classical counterparts. Good agreement with the experimental measurements in the low-frequency region is obtained for the quantum infrared spectrum, which also shows a higher intensity and a redshift relative to its classical analog. A modification of the original parametrization of the SPC/Fw model is suggested and tested in order to construct an accurate quantum model, called q-SPC/Fw, for liquid water. The quantum results for several thermodynamic and dynamical properties computed with the new model are shown to be in a significantly better agreement with the experimental data. Finally, a force-matching approach was applied to the q-SPC/Fw model to derive an effective quantum force field for liquid water in which the effects due to the nuclear quantization are explicitly distinguished from those due to the underlying molecular interactions. Thermodynamic and dynamical properties computed using standard classical simulations with this effective quantum potential are found in excellent agreement with those obtained from significantly more computationally demanding full centroid molecular dynamics simulations. The present results suggest that the inclusion of nuclear quantum effects into an empirical model for water enhances the ability of such model to faithfully represent experimental data, presumably through an increased ability of the model itself to capture realistic physical effects.     

Including quantum effects in the dynamics of complex (i.e., large) molecular systems

The development in the 1950s and 1960s of crossed molecular beam methods for studying chemical reactions at the single-collision molecular level stimulated the need and desire for theoretical methods to describe these and other dynamical processes in molecular systems. Chemical dynamics theory has made great strides in the ensuing decades, so that methods are now available for treating the quantum dynamics of small molecular systems essentially completely. For the large molecular systems that are of so much interest nowadays (e.g., chemical reactions in solution, in clusters, in nanostructures, in biological systems, etc.), however, the only generally available theoretical approach is classical molecular dynamics (MD) simulations. Much effort is currently being devoted to the development of approaches for describing the quantum dynamics of these complex systems. This paper reviews some of these approaches, especially the use of semiclassical approximations for adding quantum effects to classical MD simulations, also showing some new versions that should make these semiclassical approaches even more practical and accurate.     

An approach for generating trajectory-based dynamics which conserves the canonical distribution in the phase space formulation of quantum mechanics. II. Thermal correlation functions

We show the exact expression of the quantum mechanical time correlation function in the phase space formulation of quantum mechanics. The trajectory-based dynamics that conserves the quantum canonical distribution-equilibrium Liouville dynamics (ELD) proposed in Paper I is then used to approximately evaluate the exact expression. It gives exact thermal correlation functions (of even nonlinear operators, i.e., nonlinear functions of position or momentum operators) in the classical, high temperature, and harmonic limits. Various methods have been presented for the implementation of ELD. Numerical tests of the ELD approach in the Wigner or Husimi phase space have been made for a harmonic oscillator and two strongly anharmonic model problems, for each potential autocorrelation functions of both linear and nonlinear operators have been calculated. It suggests ELD can be a potentially useful approach for describing quantum effects for complex systems in condense phase.     

A semiclassical initial-value representation for quantum propagator and boltzmann operator

Starting from the position-momentum integral representation, we apply the correction operator method to the derivation of a uniform semiclassical approximation for the quantum propagator and then extend it to approximate the Boltzmann operator. In this approach, the involved classical dynamics is determined by the method itself instead of given beforehand. For the approximate Boltzmann operator, the corresponding classical dynamics is governed by a complex Hamiltonian, which can be described as a pair of real Hamiltonian systems. It is demonstrated that the semiclassical Boltzmann operator is exact for linear systems. A quantum propagator in the complex time is thus proposed and preliminary numerical results show that it is a reasonable approximation for calculating thermal correlation functions of general systems. © 2018 Wiley Periodicals, Inc.     

Mixed quantum-classical equilibrium

We present an analysis of the equilibrium limits of the two most widely used approaches for simulating the dynamics of molecular systems that combine both quantum and classical degrees of freedom. For a two-level quantum system connected to an infinite number of classical particles, we derive a simple analytical expression for the equilibrium mean energy attained by the self-consistent-field (Ehrenfest) method and show that it deviates substantially from Boltzmann. By contrast, "fewest switches" surface hopping achieves Boltzmann quantum state populations. We verify these analytical results with simulations.     

Time scales and other problems in linking simulations of simple chemical systems to more complex ones

We shall start with very small systems like H₂ and H₃, computed with very accurate methods (Hylleraas–CI ) or atomic systems up to Zn with accurate methods (CI ), then move to more complex ones, like C₆₀, but now with somewhat less accurate methods, specifically Hartree–Fock with density functionals, the latter for the correlation energy but not for the exchange energy. For even more complex tasks like geometry optimization of C₆₀, we have resorted to even simpler and parametrized methods, like local density functionals. Then, we could use quantum mechanics either to provide interaction potentials for classical molecular dynamics or to directly solve dynamical systems, in a quantum molecular dynamics approximation. Having demonstrated that we can use the computational output from small systems as input to larger ones, we discuss in detail a new model for liquid water, which is borne out entirely from ab initio methods and nicely links spectroscopic, thermodynamics, and other physicochemical data. Concerning time scales, we use classical molecular dynamics to determine friction coefficients, and with these we perform stochastic dynamic simulations. The use of simulation results from smaller systems to provide inputs for larger system simulations is the “global simulation” approach, which, today, with the easily available computers, is becoming more and more feasible. Projections on simulations in the 1996–1998 period are discussed, new computational areas are outlined, and a N⁴ complexity algorithm is compared to density functional approaches. © 1993 John Wiley & Sons, Inc.     

Symmetry-adapted correlation function for semiclassical quantization

We study a very simple method to incorporate quantum-mechanical symmetries, including the permutational symmetry on an equal footing with spatial symmetries, into the semiclassical calculation of correlation functions. This method is applied to the calculation of energy spectra to verify its validity by reproducing quantum energy levels for systems of bosons (symmetrized) and fermions (antisymmetrized). The mechanism of how the phase-space structure of classical dynamics is linked with the relevant quantum symmetry is discussed.     

Direct assessment of quantum nuclear effects on hydrogen bond strength by constrained-centroid ab initio path integral molecular dynamics

The impact of quantum nuclear effects on hydrogen (H-) bond strength has been inferred in earlier work from bond lengths obtained from path integral molecular dynamics (PIMD) simulations. To obtain a direct quantitative assessment of such effects, we use constrained-centroid PIMD simulations to calculate the free energy changes upon breaking the H-bonds in dimers of HF and water. Comparing ab initio simulations performed using PIMD and classical nucleus molecular dynamics (MD), we find smaller dissociation free energies with the PIMD method. Specifically, at 50 K, the H-bond in (HF)(2) is about 30% weaker when quantum nuclear effects are included, while that in (H(2)O)(2) is about 15% weaker. In a complementary set of simulations, we compare unconstrained PIMD and classical nucleus MD simulations to assess the influence of quantum nuclei on the structures of these systems. We find increased heavy atom distances, indicating weakening of the H-bond consistent with that observed by direct calculation of the free energies of dissociation.     

Finite temperature application of the corrected propagator method to reactive dynamics in a condensed-phase environment

The recently proposed mixed quantum-classical method is extended to applications at finite temperatures. The method is designed to treat complex systems consisting of a low-dimensional quantum part (the primary system) coupled to a dissipative bath described classically. The method is based on a formalism showing how to systematically correct the approximate zeroth-order evolution rule. The corrections are defined in terms of the total quantum Hamiltonian and are taken to the classical limit by introducing the frozen Gaussian approximation for the bath degrees of freedom. The evolution of the primary system is governed by the corrected propagator yielding the exact quantum dynamics. The method has been tested on a standard model system describing proton transfer in a condensed-phase environment: a symmetric double-well potential bilinearly coupled to a bath of harmonic oscillators. Flux correlation functions and thermal rate constants have been calculated at two different temperatures for a range of coupling strengths. The results have been compared to the fully quantum simulations of Topaler and Makri [J. Chem. Phys. 101, 7500 (1994)] with the real path integral method.     

Two more approaches for generating trajectory-based dynamics which conserves the canonical distribution in the phase space formulation of quantum mechanics

We show two more approaches for generating trajectory-based dynamics in the phase space formulation of quantum mechanics: "equilibrium continuity dynamics" (ECD) in the spirit of the phase space continuity equation in classical mechanics, and "equilibrium Hamiltonian dynamics" (EHD) in the spirit of the Hamilton equations of motion in classical mechanics. Both ECD and EHD can recover exact thermal correlation functions (of even nonlinear operators, i.e., nonlinear functions of position or momentum operators) in the classical, high temperature, and harmonic limits. Both ECD and EHD conserve the quasi-probability within the infinitesimal volume dx(t)dp(t) around the phase point (x(t), p(t)) along the trajectory. Numerical tests of both approaches in the Wigner phase space have been made for two strongly anharmonic model problems and a double well system, for each potential auto-correlation functions of both linear and nonlinear operators have been calculated. The results suggest EHD and ECD are two additional potential useful approaches for describing quantum effects for complex systems in condense phase.     

On the correspondence between classical and quantum mechanics in macromolecular systems: Too much classical chaos

In a recent study we found the classical dynamics of a polyethylene (PE) chain to exhibit low dimensional chaos at temperatures as low as a few Kelvin. These results strongly suggest that classical molecular dynamic simulations in polymer systems can grossly overestimate vibrational motion, which consequently results in disordered structures. In contrast, quantum mechanical calculations using Internal Coordinate Quantum Monte Carlo (an improved method with an initial conjecture for the correct wave function) indicate that the quantum ground state for a three-dimensional model PE chain is far more rigid than determined from molecular dynamics (MD) simulations, even at energies as low as a small fraction of the ground state energy. This result casts uncertainty on the reliability of MD estimates of dynamical or structural quantities relevant to the study of some macromolecular systems.     

Quantum effects in ice Ih

Quantum and classical simulations are carried out on ice Ih over a range of temperatures utilizing the TIP4P water model. The rigid-body centroid molecular dynamics method employed allows for the investigation of equilibrium and dynamical properties of the quantum system. The impact of quantization on the local structure, as measured by the radial and spatial distribution functions, as well as the energy is presented. The effects of quantization on the lattice vibrations, associated with the molecular translations and librations, are also reported. Comparison of quantum and classical simulation results indicates that shifts in the average potential energy are equivalent to rising the temperature about 80 K and are therefore non-negligible. The energy shifts due to quantization and the quantum mechanical uncertainties observed in ice are smaller than the values previously reported for liquid water. Additionally, we carry out a comparative study of melting in our classical and quantum simulations and show that there are significant differences between classical and quantum ice.