The Wolynes theory of electronically nonadiabatic reaction rates [P. G. Wolynes, J. Chem. Phys. 87, 6559 (1987)] is based on a saddle point approximation to the time integral of a reactive flux autocorrelation function in the nonadiabatic (golden rule) limit. The dominant saddle point is on the imaginary time axis at t=iλℏ, and provided λ lies in the range -β/2≤λ≤β/2, it is straightforward to evaluate the rate constant using information obtained from an imaginary time path integral calculation. However, if λ lies outside this range, as it does in the Marcus inverted regime, the path integral diverges. This has led to claims in the literature that Wolynes theory cannot describe the correct behaviour in the inverted regime. Here we show how the imaginary time correlation function obtained from a path integral calculation can be analytically continued to λ<-β/2, and the continuation used to evaluate the rate in the inverted regime. Comparison with exact golden rule results for a spin-boson model and a more demanding (asymmetric and anharmonic) model of electronic predissociation shows that the theory is just as accurate in the inverted regime as it is in the normal regime.

We present a simple interpolation formula for the rate of an electron transfer reaction as a function of the electronic coupling strength. The formula only requires the calculation of Fermi Golden Rule and Born-Oppenheimer rates and so can be combined with any methods that are able to calculate these rates. We first demonstrate the accuracy of the formula by applying it to a one dimensional scattering problem for which the exact quantum mechanical, Fermi Golden Rule, and Born-Oppenheimer rates are readily calculated. We then describe how the formula can be combined with the Wolynes theory approximation to the Golden Rule rate, and the ring polymer molecular dynamics (RPMD) approximation to the Born-Oppenheimer rate, and used to capture the effects of nuclear tunnelling, zero point energy, and solvent friction on condensed phase electron transfer reactions. Comparison with exact hierarchical equations of motion (HEOM) results for a demanding set of spin-boson models shows that the interpolation formula has an error comparable to that of RPMD rate theory in the adiabatic limit, and that of Wolynes theory in non-adiabatic limit, and is therefore as accurate as any method could possibly be that attempts to generalise these methods to arbitrary electronic coupling strengths. arXiv:1909.09882v1 [physics.chem-ph]

We present a new method for calculating the product yield of a radical pair recombination reaction in the presence of a weak time-dependent magnetic field. This method successfully circumvents the computational difficulties presented by a direct solution of the Liouville-von Neumann equation for a long-lived radical pair containing many hyperfine-coupled nuclear spins. Using a modified formulation of Floquet theory, treating the time-dependent magnetic field as a perturbation, and exploiting the slow radical pair recombination, we show that one can obtain a good approximation to the product yield by considering only nearly degenerate sub-spaces of the Floquet space. Within a significant parameter range, the resulting method is found to give product yields in good agreement with exact quantum mechanical results for a variety of simple model radical pairs. Moreover it is considerably more efficient than the exact calculation, and it can be applied to radical pairs containing significantly more nuclear spins. This promises to open the door to realistic theoretical investigations of the effect of radiofrequency electromagnetic radiation on the photochemically induced radical pair recombination reactions in the avian retina which are believed to be responsible for the magnetic compass sense of migratory birds. Published by AIP Publishing. [http://dx.

We shall use this introduction to the Faraday Discussion on quantum effects in complex systems to review the recent progress that has been made in using imaginary time path integral methods to calculate chemical reaction rates.

We describe a fast implementation of the quasi-centroid molecular dynamics (QCMD) method in which the quasi-centroid potential of mean force is approximated as a separable correction to the classical interaction potential. This correction is obtained by first calculating quasi-centroid radial and angular distribution functions in a short path integral molecular dynamics simulation and then using iterative Boltzmann inversion to obtain an effective classical potential that reproduces these distribution functions in a classical NVT simulation. We illustrate this approach with example applications to the vibrational spectra of gas phase molecules, obtaining excellent agreement with QCMD reference calculations for water and ammonia and good agreement with the quantum mechanical vibrational spectrum of methane.