Heavy-atom tunnelling in singlet oxygen deactivation predicted by instanton theory with branch-point singularities.

The reactive singlet state of oxygen (O2) can decay to the triplet ground state nonradiatively in the presence of a solvent. There is a controversy about whether tunnelling is involved in this nonadiabatic spin-crossover process. Semiclassical instanton theory provides a reliable and practical computational method for elucidating the reaction mechanism and can account for nuclear quantum effects such as zero-point energy and multidimensional tunnelling. However, the previously developed instanton theory is not directly applicable to this system because of a branch-point singularity which appears in the flux correlation function. Here we derive a new instanton theory for cases dominated by the singularity, leading to a new picture of tunnelling in nonadiabatic processes. Together with multireference electronic-structure theory, this provides a rigorous framework based on first principles that we apply to calculate the decay rate of singlet oxygen in water. The results indicate a new reaction mechanism that is 27 orders of magnitude faster at room temperature than the classical process through the minimum-energy crossing point. We find significant heavy-atom tunnelling contributions as well as a large temperature-dependent H2O/D2O kinetic isotope effect of approximately 20, in excellent agreement with experiment.

Exact tunneling splittings from symmetrized path integrals.

We develop a new simulation technique based on path-integral molecular dynamics for calculating ground-state tunneling splitting patterns from ratios of symmetrized partition functions. In particular, molecular systems are rigorously projected onto their J = 0 rotational state by an "Eckart spring" that connects two adjacent beads in a ring polymer. Using this procedure, the tunneling splitting can be obtained from thermodynamic integration at just one (sufficiently low) temperature. Converged results are formally identical to the values that would have been obtained by solving the full rovibrational Schrödinger equation on a given Born-Oppenheimer potential energy surface. The new approach is showcased with simulations of hydronium and methanol, which are in good agreement with wavefunction-based calculations and experimental measurements. The method will be of particular use for the study of low-barrier methyl rotations and other floppy modes, where instanton theory is not valid.

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.

Nonadiabatic instanton rate theory beyond the golden-rule limit.

Fermi's golden rule (GR) describes the leading-order behavior of the reaction rate as a function of the diabatic coupling. Its asymptotic (ℏ → 0) limit is the semiclassical golden-rule instanton rate theory, which rigorously approximates nuclear quantum effects, lends itself to efficient numerical computation, and gives physical insight into reaction mechanisms. However, the golden rule by itself becomes insufficient as the strength of the diabatic coupling increases, so higher-order terms must be additionally considered. In this work, we give a first-principles derivation of the next-order term beyond the golden rule, represented as a sum of three components. Two of them lead to new instanton pathways that extend the GR case and, among other factors, account for effects of recrossing on the full rate. The remaining component derives from the equilibrium partition function and accounts for changes in potential energy around the reactant and product wells due to diabatic coupling. The new semiclassical theory demands little computational effort beyond a GR instanton calculation. It makes it possible to rigorously assess the accuracy of the GR approximation and sets the stage for future work on general semiclassical nonadiabatic rate theories.

Instanton theory for Fermi's golden rule and beyond.

Instanton theory provides a semiclassical approximation for computing quantum tunnelling effects in complex molecular systems. It is typically applied to proton-transfer reactions for which the Born-Oppenheimer approximation is valid. However, many processes in physics, chemistry and biology, such as electron transfers, are non-adiabatic and are correctly described instead using Fermi's golden rule. In this work, we discuss how instanton theory can be generalized to treat these reactions in the golden-rule limit. We then extend the theory to treat fourth-order processes such as bridge-mediated electron transfer and apply the method to simulate an electron moving through a model system of three coupled quantum dots. By comparison with benchmark quantum calculations, we demonstrate that the instanton results are much more reliable than alternative approximations based on superexchange-mediated effective coupling or a classical sequential mechanism. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.

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.