Instantons and the quantum bound to chaos.

The rate at which information scrambles in a quantum system can be quantified using out-of-time-ordered correlators. A remarkable prediction is that the associated Lyapunov exponent [Formula: see text] that quantifies the scrambling rate in chaotic systems obeys a universal bound [Formula: see text]. Previous numerical and analytical studies have indicated that this bound has a quantum-statistical origin. Here, we use path-integral techniques to show that a minimal theory to reproduce this bound involves adding contributions from quantum thermal fluctuations (describing quantum tunneling and zero-point energy) to classical dynamics. By propagating a model quantum-Boltzmann-conserving classical dynamics for a system with a barrier, we show that the bound is controlled by the stability of thermal fluctuations around the barrier instanton (a delocalized structure which dominates the tunneling statistics). This stability requirement appears to be general, implying that there is a close relation between the formation of instantons, or related delocalized structures, and the imposition of the quantum-chaos bound.

Path Integral Simulations of Condensed-Phase Vibrational Spectroscopy.

Recent theoretical and algorithmic developments have improved the accuracy with which path integral dynamics methods can include nuclear quantum effects in simulations of condensed-phase vibrational spectra. Such methods are now understood to be approximations to the delocalized classical Matsubara dynamics of smooth Feynman paths, which dominate the dynamics of systems such as liquid water at room temperature. Focusing mainly on simulations of liquid water and hexagonal ice, we explain how the recently developed quasicentroid molecular dynamics (QCMD), fast-QCMD, and temperature-elevated path integral coarse-graining simulations (Te PIGS) methods generate classical dynamics on potentials of mean force obtained by averaging over quantum thermal fluctuations. These new methods give very close agreement with one another, and the Te PIGS method has recently yielded excellent agreement with experimentally measured vibrational spectra for liquid water, ice, and the liquid-air interface. We also discuss the limitations of such methods.

Quantum rates in dissipative systems with spatially varying friction.

We investigate whether making the friction spatially dependent on the reaction coordinate introduces quantum effects into the thermal reaction rates for dissipative reactions. Quantum rates are calculated using the numerically exact multi-configuration time-dependent Hartree method, as well as the approximate ring-polymer molecular dynamics (RPMD), ring-polymer instanton methods, and classical molecular dynamics. By conducting simulations across a wide range of temperatures and friction strengths, we can identify the various regimes that govern the reactive dynamics. At high temperatures, in addition to the spatial-diffusion and energy-diffusion regimes predicted by Kramer's rate theory, a (coherent) tunneling-dominated regime is identified at low friction. At low temperatures, incoherent tunneling dominates most of Kramer's curve, except at very low friction, when coherent tunneling becomes dominant. Unlike in classical mechanics, the bath's influence changes the equilibrium time-independent properties of the system, leading to a complex interplay between spatially dependent friction and nuclear quantum effects even at high temperatures. More specifically, a realistic friction profile can lead to an increase (or decrease) of the quantum (classical) rates with friction within the spatial-diffusion regime, showing that classical and quantum rates display qualitatively different behaviors. Except at very low frictions, we find that RPMD captures most of the quantum effects in the thermal reaction rates.

Comparison of Matsubara dynamics with exact quantum dynamics for an oscillator coupled to a dissipative bath.

We report the first numerical calculations in which converged Matsubara dynamics is compared directly with exact quantum dynamics with no artificial damping of the time-correlation functions (TCFs). The system treated is a Morse oscillator coupled to a harmonic bath. We show that, when the system-bath coupling is sufficiently strong, the Matsubara calculations can be converged by explicitly including up to M = 200 Matsubara modes, with the remaining modes included as a harmonic "tail" correction. The resulting Matsubara TCFs are in near-perfect agreement with the exact quantum TCFs, for non-linear as well as linear operators, at a temperature at which the TCFs are dominated by quantum thermal fluctuations. These results provide compelling evidence that incoherent classical dynamics can arise in the condensed phase at temperatures at which the statistics are dominated by quantum (Boltzmann) effects, as a result of smoothing of imaginary-time Feynman paths. The techniques developed here may also lead to efficient methods for benchmarking system-bath dynamics in the overdamped regime.

Improved torque estimator for condensed-phase quasicentroid molecular dynamics.

We describe improvements to the quasicentroid molecular dynamics (QCMD) path-integral method, which was developed recently for computing the infrared spectra of condensed-phase systems. The main development is an improved estimator for the intermolecular torque on the quasicentroid. When applied to qTIP4P/F liquid water and ice, the new estimator is found to remove an artificial 25 cm-1 red shift from the libration bands, to increase slightly the intensity of the OH stretch band in the liquid, and to reduce small errors noted previously in the QCMD radial distribution functions. We also modify the mass-scaling used in the adiabatic QCMD algorithm, which allows the molecular dynamics timestep to be quadrupled, thus reducing the expense of a QCMD calculation to twice that of Cartesian centroid molecular dynamics for qTIP4P/F liquid water at 300 K, and eight times for ice at 150 K.

On the "Matsubara heating" of overtone intensities and Fermi splittings.

Classical molecular dynamics (MD) and imaginary-time path-integral dynamics methods underestimate the infrared absorption intensities of overtone and combination bands by typically an order of magnitude. Plé et al. [J. Chem. Phys. 155, 2863 (2021)] have shown that this is because such methods fail to describe the coupling of the centroid to the Matsubara dynamics of the fluctuation modes; classical first-order perturbation theory (PT) applied to the Matsubara dynamics is sufficient to recover most of the lost intensity in simple models and gives identical results to quantum (Rayleigh-Schrödinger) PT. Here, we show numerically that the results of this analysis can be used as post-processing correction factors, which can be applied to realistic (classical MD or path-integral dynamics) simulations of infrared spectra. We find that the correction factors recover most of the lost intensity in the overtone and combination bands of gas-phase water and ammonia and much of it for liquid water. We then re-derive and confirm the earlier PT analysis by applying canonical PT to Matsubara dynamics, which has the advantage of avoiding secular terms and gives a simple picture of the perturbed Matsubara dynamics in terms of action-angle variables. Collectively, these variables "Matsubara heat" the amplitudes of the overtone and combination vibrations of the centroid to what they would be in a classical system with the oscillators (of frequency Ωi) held at their quantum effective temperatures [of ℏΩi coth(ßℏΩi/2)/2kB]. Numerical calculations show that a similar neglect of "Matsubara heating" causes path-integral methods to underestimate Fermi resonance splittings.

Testing the quasicentroid molecular dynamics method on gas-phase ammonia.

Quasicentroid molecular dynamics (QCMD) is a path-integral method for approximating nuclear quantum effects in dynamics simulations, which has given promising results for gas- and condensed-phase water. In this work, by simulating the infrared spectrum of gas-phase ammonia, we test the feasibility of extending QCMD beyond water. Overall, QCMD works as well for ammonia as for water, reducing or eliminating blue shifts from the classical spectrum without introducing the artificial red shifts or broadening associated with other imaginary-time path-integral methods. However, QCMD gives only a modest improvement over the classical spectrum for the position of the symmetric bend mode, which is highly anharmonic (since it correlates with the inversion pathway). We expect QCMD to have similar problems with large-amplitude degrees of freedom in other molecules but otherwise to work as well as for water.

Which quantum statistics-classical dynamics method is best for water?

There are a variety of methods for including nuclear quantum effects in dynamics simulations by combining quantum Boltzmann statistics with classical dynamics. Among them are thermostatted ring-polymer molecular dynamics (TRPMD), centroid molecular dynamics (CMD), quasi-centroid molecular dynamics (QCMD), and the linearised semi-classical initial value representation (LSC-IVR). Here we make a systematic comparison of these methods by calculating the infrared spectrum of water in the gas phase, and in the liquid and ice phases (using the q-TIP4P/F model potential). Some of these results are taken from previous work, and some of them are new (including the LSC-IVR calculations for ice, and extensions of all the spectra into the near-infrared region dominated by overtone and combination bands). Our results suggest that QCMD is the best method for reproducing fundamental transitions in the spectrum, and that LSC-IVR gives the best overall description of the spectrum (albeit with large errors in the bend fundamental band caused by zero-point-energy leakage). The TRPMD method gives damped spectra that line up with the QCMD spectra, and is by far the cheapest method.

Mean-field Matsubara dynamics: Analysis of path-integral curvature effects in rovibrational spectra.

It was shown recently that smooth and continuous "Matsubara" phase-space loops follow a quantum-Boltzmann-conserving classical dynamics when decoupled from non-smooth distributions, which was suggested as the reason that many dynamical observables appear to involve a mixture of classical dynamics and quantum Boltzmann statistics. Here we derive a mean-field version of this "Matsubara dynamics" which sufficiently mitigates its serious phase problem to permit numerical tests on a two-dimensional "champagne-bottle" model of a rotating OH bond. The Matsubara-dynamics rovibrational spectra are found to converge toward close agreement with the exact quantum results at all temperatures tested (200-800 K), the only significant discrepancies being a temperature-independent 22 cm-1 blue-shift in the position of the vibrational peak and a slight broadening in its line shape. These results are compared with centroid molecular dynamics (CMD) to assess the importance of non-centroid fluctuations. Above 250 K, only the lowest-frequency non-centroid modes are needed to correct small CMD red-shifts in the vibrational peak; below 250 K, more non-centroid modes are needed to correct large CMD red-shifts and broadening. The transition between these "shallow curvature" and "deep curvature" regimes happens when imaginary-time Feynman paths become able to lower their actions by cutting through the curved potential surface, giving rise to artificial instantons in CMD.

Approximating Matsubara dynamics using the planetary model: Tests on liquid water and ice.

Matsubara dynamics is the quantum-Boltzmann-conserving classical dynamics which remains when real-time coherences are taken out of the exact quantum Liouvillian [T. J. H. Hele et al., J. Chem. Phys. 142, 134103 (2015)]; because of a phase-term, it cannot be used as a practical method without further approximation. Recently, Smith et al. [J. Chem. Phys. 142, 244112 (2015)] developed a "planetary" model dynamics which conserves the Feynman-Kleinert (FK) approximation to the quantum-Boltzmann distribution. Here, we show that for moderately anharmonic potentials, the planetary dynamics gives a good approximation to Matsubara trajectories on the FK potential surface by decoupling the centroid trajectory from the locally harmonic Matsubara fluctuations, which reduce to a single phase-less fluctuation particle (the "planet"). We also show that the FK effective frequency can be approximated by a direct integral over these fluctuations, obviating the need to solve iterative equations. This modification, together with use of thermostatted ring-polymer molecular dynamics, allows us to test the planetary model on water (gas-phase, liquid, and ice) using the q-TIP4P/F potential surface. The "planetary" fluctuations give a poor approximation to the rotational/librational bands in the infrared spectrum, but a good approximation to the bend and stretch bands, where the fluctuation lineshape is found to be motionally narrowed by the vibrations of the centroid.

Is the simplest chemical reaction really so simple?

Modern computational methods have become so powerful for predicting the outcome for the H + H2 → H2 + H bimolecular exchange reaction that it might seem further experiments are not needed. Nevertheless, experiments have led the way to cause theorists to look more deeply into this simplest of all chemical reactions. The findings are less simple.

Rovibrational transitions of the methane-water dimer from intermolecular quantum dynamical computations.

Rovibrational quantum nuclear motion computations, with J = 0, 1, and 2, are reported for the intermolecular degrees of freedom of the methane-water dimer, where J is the quantum number describing the overall rotation of the complex. The computations provide the first explanation of the far-infrared spectrum of this complex published in J. Chem. Phys., 1994, 100, 863. All experimentally reported rovibrational transitions, up to J = 2, can be assigned to transitions between the theoretically computed levels. The deviation of the experimental and computed rovibrational transitions is 0.5 cm(-1) for the ortho and 2 cm(-1) for the para species with a variance of 0.005 cm(-1). In addition to a lower systematic error, the overall agreement of theory and experiment is also better for the ortho species (involving ortho-H2O). Most importantly, for this species all levels of the 24-fold tunneling splitting manifold corresponding to the zero-point vibration (ZPV) are involved in at least one experimentally reported transition. For the para species there are a few energy levels in the computed ZPV manifold that are not involved in the reported experimental transitions. Furthermore, computed energy levels are identified that correspond to the ZPV tunneling splitting manifold of the secondary minimum structure of the dimer, which presumably appear in rovibrational transitions in the same energy regime as the observed transitions, but have not been experimentally reported.

Calculating splittings between energy levels of different symmetry using path-integral methods.

It is well known that path-integral methods can be used to calculate the energy splitting between the ground and the first excited state. Here we show that this approach can be generalized to give the splitting patterns between all the lowest energy levels from different symmetry blocks that lie below the first-excited totally symmetric state. We demonstrate this property numerically for some two-dimensional models. The approach is likely to be useful for computing rovibrational energy levels and tunnelling splittings in floppy molecules and gas-phase clusters.

An alternative derivation of ring-polymer molecular dynamics transition-state theory.

In a previous article [T. J. H. Hele and S. C. Althorpe, J. Chem. Phys. 138, 084108 (2013)], we showed that the t → 0+ limit of ring-polymer molecular dynamics (RPMD) rate-theory is also the t → 0+ limit of a new type of quantum flux-side time-correlation function, in which the dividing surfaces are invariant to imaginary-time translation; in other words, that RPMD transition-state theory (RMPD-TST) is a t → 0+ quantum transition-state theory (QTST). Recently, Jang and Voth [J. Chem. Phys. 144, 084110 (2016)] rederived this quantum t → 0+ limit and claimed that it gives instead the centroid-density approximation. Here we show that the t → 0+ limit derived by Jang and Voth is in fact RPMD-TST.

Quantum tunneling splittings from path-integral molecular dynamics.

We illustrate how path-integral molecular dynamics can be used to calculate ground-state tunnelling splittings in molecules or clusters. The method obtains the splittings from ratios of density matrix elements between the degenerate wells connected by the tunnelling. We propose a simple thermodynamic integration scheme for evaluating these elements. Numerical tests on fully dimensional malonaldehyde yield tunnelling splittings in good overall agreement with the results of diffusion Monte Carlo calculations.

Non-equilibrium dynamics from RPMD and CMD.

We investigate the calculation of approximate non-equilibrium quantum time correlation functions (TCFs) using two popular path-integral-based molecular dynamics methods, ring-polymer molecular dynamics (RPMD) and centroid molecular dynamics (CMD). It is shown that for the cases of a sudden vertical excitation and an initial momentum impulse, both RPMD and CMD yield non-equilibrium TCFs for linear operators that are exact for high temperatures, in the t = 0 limit, and for harmonic potentials; the subset of these conditions that are preserved for non-equilibrium TCFs of non-linear operators is also discussed. Furthermore, it is shown that for these non-equilibrium initial conditions, both methods retain the connection to Matsubara dynamics that has previously been established for equilibrium initial conditions. Comparison of non-equilibrium TCFs from RPMD and CMD to Matsubara dynamics at short times reveals the orders in time to which the methods agree. Specifically, for the position-autocorrelation function associated with sudden vertical excitation, RPMD and CMD agree with Matsubara dynamics up to O(t4) and O(t1), respectively; for the position-autocorrelation function associated with an initial momentum impulse, RPMD and CMD agree with Matsubara dynamics up to O(t5) and O(t2), respectively. Numerical tests using model potentials for a wide range of non-equilibrium initial conditions show that RPMD and CMD yield non-equilibrium TCFs with an accuracy that is comparable to that for equilibrium TCFs. RPMD is also used to investigate excited-state proton transfer in a system-bath model, and it is compared to numerically exact calculations performed using a recently developed version of the Liouville space hierarchical equation of motion approach; again, similar accuracy is observed for non-equilibrium and equilibrium initial conditions.

Differential Cross Sections for the H + D2 → HD(v' = 3, j' = 4-10) + D Reaction above the Conical Intersection.

We report rovibrationally selected differential cross sections (DCSs) of the benchmark reaction H + D2 → HD(v' = 3, j' = 4-10) + D at a collision energy of 3.26 eV, which exceeds the conical intersection of the H3 potential energy surface at 2.74 eV. We use the PHOTOLOC technique in which a fluorine excimer laser at 157.64 nm photodissociates hydrogen bromide (HBr) molecules to generate fast H atoms and the HD product is detected in a state-specific manner by resonance-enhanced multiphoton ionization. Fully converged quantum wave packet calculations were performed for this reaction at this high collision energy without inclusion of the geometric phase (GP) effect, which takes into account coupling to the first excited state of the H3 potential energy surface. Multimodal structures can be observed in most of the DCSs up to j' = 10, which is predicted by theory and also well-reproduced by experiment. The theoretically calculated DCSs are in good overall agreement with the experimental measurements, which indicates that the GP effect is not large enough that its existence can be verified experimentally at this collision energy.

Boltzmann-conserving classical dynamics in quantum time-correlation functions: "Matsubara dynamics".

We show that a single change in the derivation of the linearized semiclassical-initial value representation (LSC-IVR or "classical Wigner approximation") results in a classical dynamics which conserves the quantum Boltzmann distribution. We rederive the (standard) LSC-IVR approach by writing the (exact) quantum time-correlation function in terms of the normal modes of a free ring-polymer (i.e., a discrete imaginary-time Feynman path), taking the limit that the number of polymer beads N → ∞, such that the lowest normal-mode frequencies take their "Matsubara" values. The change we propose is to truncate the quantum Liouvillian, not explicitly in powers of h(2) at h(0) (which gives back the standard LSC-IVR approximation), but in the normal-mode derivatives corresponding to the lowest Matsubara frequencies. The resulting "Matsubara" dynamics is inherently classical (since all terms O(h(2)) disappear from the Matsubara Liouvillian in the limit N → ∞) and conserves the quantum Boltzmann distribution because the Matsubara Hamiltonian is symmetric with respect to imaginary-time translation. Numerical tests show that the Matsubara approximation to the quantum time-correlation function converges with respect to the number of modes and gives better agreement than LSC-IVR with the exact quantum result. Matsubara dynamics is too computationally expensive to be applied to complex systems, but its further approximation may lead to practical methods.

Communication: Relation of centroid molecular dynamics and ring-polymer molecular dynamics to exact quantum dynamics.

We recently obtained a quantum-Boltzmann-conserving classical dynamics by making a single change to the derivation of the "Classical Wigner" approximation. Here, we show that the further approximation of this "Matsubara dynamics" gives rise to two popular heuristic methods for treating quantum Boltzmann time-correlation functions: centroid molecular dynamics (CMD) and ring-polymer molecular dynamics (RPMD). We show that CMD is a mean-field approximation to Matsubara dynamics, obtained by discarding (classical) fluctuations around the centroid, and that RPMD is the result of discarding a term in the Matsubara Liouvillian which shifts the frequencies of these fluctuations. These findings are consistent with previous numerical results and give explicit formulae for the terms that CMD and RPMD leave out.