We propose an approximate method for calculating Kubo-transformed real-time correlation functions involving position-dependent operators, based on path integral (Parrinello-Rahman) molecular dynamics. The method gives the exact quantum mechanical correlation function at time zero, exactly satisfies the quantum mechanical detailed balance condition, and for correlation functions of the form C(Ax)(t) and C(xB)(t) it gives the exact result for a harmonic potential. It also works reasonably well at short times for more general potentials and correlation functions, as we illustrate with some example calculations. The method provides a consistent improvement over purely classical molecular dynamics that is most apparent in the low-temperature regime.

We further develop the ring polymer molecular dynamics (RPMD) method for calculating chemical reaction rates [I. R. Craig and D. E. Manolopoulos, J. Chem. Phys. 122, 084106 (2005)]. We begin by showing how the rate coefficient we obtained before can be calculated in a more efficient way by considering the side functions of the ring-polymer centroids, rather than averaging over the side functions of the individual ring-polymer beads. This has two distinct advantages. First, the statistics of the phase-space average over the ring-polymer coordinates and momenta are greatly improved. Second, the resulting flux-side correlation function converges to its long-time limit much more rapidly. Indeed the short-time limit of this flux-side correlation function already provides a "quantum transition state theory" approximation to the final rate coefficient. In cases where transition state recrossing effects are negligible, and the transition state dividing surface is put in the right place, the RPMD rate is therefore obtained almost instantly. We then go on to show that the long-time limit of the new flux-side correlation function, and hence the fully converged RPMD reaction rate, is rigorously independent of the choice of the transition state dividing surface. This is especially significant because the optimum dividing surface can often be very difficult to determine for reactions in complex chemical systems.

We show how the ring-polymer molecular dynamics method can be adapted to calculate approximate Kubo-transformed flux-side correlation functions, and hence rate coefficients for condensed phase reactions. An application of the method to the standard model for a chemical reaction in solution--a quartic double-well potential linearly coupled to a bath of harmonic oscillators--is found to give results of comparable accuracy to those of the classical Wigner model and the centroid molecular dynamics method. However, since the present method does not require that one evaluate the Wigner transform of a thermal flux operator or that one perform a separate path integral calculation for each molecular dynamics time step, we believe it will prove easier to apply to more general problems than either of these alternative techniques. We also present a (logarithmic) discretization scheme for the Ohmic bath in the system-bath model that gives converged results with just nine bath modes--a surprisingly small number for a model of a condensed phase reaction. Finally, we present some calculations of the transmission through an Eckart barrier which show that the present method provides a satisfactory (although not perfect) description of the deep quantum tunneling regime. Part of the reason for the success of the method is that it gives the exact quantum-mechanical rate constant for the transmission through a parabolic barrier, as we demonstrate analytically in the Appendix.

We have used the ring polymer molecular dynamics method to study the Azzouz-Borgis model for proton transfer between phenol (AH) and trimethylamine (B) in liquid methyl chloride. When the A-H distance is used as the reaction coordinate, the ring polymer trajectories are found to exhibit multiple recrossings of the transition state dividing surface and to give a rate coefficient that is smaller than the quantum transition state theory value by an order of magnitude. This is to be expected on kinematic grounds for a heavy-light-heavy reaction when the light atom transfer coordinate is used as the reaction coordinate, and it clearly precludes the use of transition state theory with this reaction coordinate. As has been shown previously for this problem, a solvent polarization coordinate defined in terms of the expectation value of the proton transfer distance in the ground adiabatic quantum state provides a better reaction coordinate with less recrossing. These results are discussed in light of the wide body of earlier theoretical work on the Azzouz-Borgis model and the considerable range of previously reported values for its proton and deuteron transfer rate coefficients.

The recently proposed multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) approach to evaluating reactive quantum dynamics is applied to two model condensed-phase proton transfer reactions. The models consist of a one-dimensional double-well "system" that is bilinearly coupled to a "bath" of harmonic oscillators parameterized to represent a condensed-phase environment. Numerically exact quantum-mechanical flux correlation functions and thermal rate constants are obtained for a broad range of temperatures and system-bath coupling strengths, thus demonstrating the efficacy of the ML-MCTDH approach. Particular attention is focused on the regime where low temperatures are combined with weak system-bath coupling. Under such conditions it is found that long propagation times are often required and that quantum coherence effects may prevent a rigorous determination of the rate constant.