Exciton dynamics from the mapping approach to surface hopping: comparison with Förster and Redfield theories.

We compare the recently introduced multi-state mapping approach to surface hopping (MASH) with the Förster and Redfield theories of excitation energy transfer. Whereas Förster theory relies on weak coupling between chromophores, and Redfield theory assumes the electronic excitations to be weakly coupled to fast chromophore vibrations, MASH is free from any perturbative or Markovian approximations. We illustrate this with an example application to the rate of energy transfer in a Frenkel-exciton dimer, showing that MASH interpolates correctly between the opposing regimes in which the Förster and Redfield results are reliable. We then compare the three methods for a realistic model of the Fenna-Matthews-Olson complex with a structured vibrational spectral density and static disorder in the excitation energies. In this case there are no exact results for comparison so we use MASH to assess the validity of Förster and Redfield theories. We find that Förster theory is the more accurate of the two on the picosecond timescale, as has been shown previously for a simpler model of this particular light-harvesting complex. We also explore various ways to sample the initial electronic state in MASH and find that they all give very similar results for exciton dynamics.

A multi-state mapping approach to surface hopping.

We describe a multiple electronic state adaptation of the mapping approach to surface hopping introduced recently by Mannouch and Richardson [J. Chem. Phys. 158, 104111 (2023)]. Our modification treats populations and coherences on an equal footing and is guaranteed to give populations in any electronic basis that tend to the correct quantum-classical equilibrium values in the long-time limit (assuming ergodicity). We demonstrate its accuracy by comparison with exact benchmark results for three- and seven-state models of the Fenna-Matthews-Olson complex, obtaining electronic populations and coherences that are significantly more accurate than those of fewest switches surface hopping and at least as good as those of any other semiclassical method we are aware of. Since these results were obtained by adapting the scheme of Mannouch and Richardson, we go on to compare our results with theirs for a variety of problems with two electronic states. We find that their method is sometimes more accurate, especially in the Marcus inverted regime. However, in other situations, the accuracies are comparable, and since our scheme can be used with multiple electronic states it can be applied to a wider variety of electronically nonadiabatic systems.

Vibrational strong coupling in liquid water from cavity molecular dynamics.

We assess the cavity molecular dynamics method for the calculation of vibrational polariton spectra using liquid water as a specific example. We begin by disputing a recent suggestion that nuclear quantum effects may lead to a broadening of polariton bands, finding instead that they merely result in anharmonic red shifts in the polariton frequencies. We go on to show that our simulated cavity spectra can be reproduced to graphical accuracy with a harmonic model that uses just the cavity-free spectrum and the geometry of the cavity as input. We end by showing that this harmonic model can be combined with the experimental cavity-free spectrum to give results in good agreement with optical cavity measurements. Since the input to our harmonic model is equivalent to the input to the transfer matrix method of applied optics, we conclude that cavity molecular dynamics cannot provide any more insight into the effect of vibrational strong coupling on the absorption spectrum than this transfer matrix method, which is already widely used by experimentalists to corroborate their cavity results.

A smooth basis for atomistic machine learning.

Machine learning frameworks based on correlations of interatomic positions begin with a discretized description of the density of other atoms in the neighborhood of each atom in the system. Symmetry considerations support the use of spherical harmonics to expand the angular dependence of this density, but there is, as of yet, no clear rationale to choose one radial basis over another. Here, we investigate the basis that results from the solution of the Laplacian eigenvalue problem within a sphere around the atom of interest. We show that this generates a basis of controllable smoothness within the sphere (in the same sense as plane waves provide a basis with controllable smoothness for a problem with periodic boundaries) and that a tensor product of Laplacian eigenstates also provides a smooth basis for expanding any higher-order correlation of the atomic density within the appropriate hypersphere. We consider several unsupervised metrics of the quality of a basis for a given dataset and show that the Laplacian eigenstate basis has a performance that is much better than some widely used basis sets and competitive with data-driven bases that numerically optimize each metric. Finally, we investigate the role of the basis in building models of the potential energy. In these tests, we find that a combination of the Laplacian eigenstate basis and target-oriented heuristics leads to equal or improved regression performance when compared to both heuristic and data-driven bases in the literature. We conclude that the smoothness of the basis functions is a key aspect of successful atomic density representations.

Nuclear quantum effects in thermal conductivity from centroid molecular dynamics.

We show that the centroid molecular dynamics (CMD) method provides a realistic way to calculate the thermal diffusivity a = λ/ρcV of a quantum mechanical liquid such as para-hydrogen. Once a has been calculated, the thermal conductivity can be obtained from λ = ρcVa, where ρ is the density of the liquid and cV is the constant-volume heat capacity. The use of this formula requires an accurate quantum mechanical heat capacity cV, which can be obtained from a path integral molecular dynamics simulation. The thermal diffusivity can be calculated either from the decay of the equilibrium density fluctuations in the liquid or by using the Green-Kubo relation to calculate the CMD approximation to λ and then dividing this by the corresponding approximation to ρcV. We show that both approaches give the same results for liquid para-hydrogen and that these results are in good agreement with the experimental measurements of the thermal conductivity over a wide temperature range. In particular, they correctly predict a decrease in the thermal conductivity at low temperatures-an effect that stems from the decrease in the quantum mechanical heat capacity and has eluded previous para-hydrogen simulations. We also show that the method gives equally good agreement with the experimental measurements for the thermal conductivity of normal liquid helium.

Fast quasi-centroid molecular dynamics.

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.

Spin relaxation in radical pairs from the stochastic Schrödinger equation.

We show that the stochastic Schrödinger equation (SSE) provides an ideal way to simulate the quantum mechanical spin dynamics of radical pairs. Electron spin relaxation effects arising from fluctuations in the spin Hamiltonian are straightforward to include in this approach, and their treatment can be combined with a highly efficient stochastic evaluation of the trace over nuclear spin states that is required to compute experimental observables. These features are illustrated in example applications to a flavin-tryptophan radical pair of interest in avian magnetoreception and to a problem involving spin-selective radical pair recombination along a molecular wire. In the first of these examples, the SSE is shown to be both more efficient and more widely applicable than a recent stochastic implementation of the Lindblad equation, which only provides a valid treatment of relaxation in the extreme-narrowing limit. In the second, the exact SSE results are used to assess the accuracy of a recently proposed combination of Nakajima-Zwanzig theory for the spin relaxation and Schulten-Wolynes theory for the spin dynamics, which is applicable to radical pairs with many more nuclear spins. We also analyze the efficiency of trace sampling in some detail, highlighting the particular advantages of sampling with SU(N) coherent states.

A general non-adiabatic quantum instanton approximation.

We present a general quantum instanton approach to calculating reaction rates for systems with two electronic states and arbitrary values of the electronic coupling. This new approach, which we call the non-adiabatic quantum instanton (NAQI) approximation, reduces to Wolynes theory in the golden rule limit and to a recently proposed projected quantum instanton method in the adiabatic limit. As in both of these earlier theories, the NAQI approach is based on making a saddle point approximation to the time integral of a reactive flux autocorrelation function, although with a generalized definition of the projection operator onto the product states. We illustrate the accuracy of the approach by comparison with exact rates for one dimensional scattering problems and discuss its applicability to more complex reactions.

Confirming the role of nuclear tunneling in aqueous ferrous-ferric electron transfer.

We revisit the well-known aqueous ferrous-ferric electron transfer reaction in order to address recent suggestions that nuclear tunneling can lead to significant deviation from the linear response assumption inherent in the Marcus picture of electron transfer. A recent study of this reaction by Richardson and co-workers [Phys. Chem. Chem. Phys. 22, 10687 (2020)] has found a large difference between their new path-integral method, golden-rule quantum transition state theory (GR-QTST), and the saddle point approximation of Wolynes (Wolynes theory). They suggested that this difference could be attributed to the existence of multiple tunneling pathways, leading Wolynes theory to significantly overestimate the rate. This was used to argue that the linear response assumptions of Marcus theory may break down for liquid systems when tunneling is important. If true, this would imply that the commonly used method for studying such systems, where the problem is mapped onto a spin-boson model, is invalid. However, we have recently shown that size inconsistency in GR-QTST can lead to poor predictions of the rate in systems with many degrees of freedom. We have also suggested an improved method, the path-integral linear golden-rule (LGR) approximation, which fixes this problem. Here, we demonstrate that the GR-QTST results for ferrous-ferric electron transfer are indeed dominated by its size consistency error. Furthermore, by comparing the LGR and Wolynes theory results, we confirm the established picture of nuclear tunneling in this system. Finally, by comparing our path-integral results to those obtained by mapping onto the spin-boson model, we reassess the importance of anharmonic effects and the accuracy of this commonly used mapping approach.

An improved path-integral method for golden-rule rates.

We present a simple method for the calculation of reaction rates in the Fermi golden-rule limit, which accurately captures the effects of tunneling and zero-point energy. The method is based on a modification of the recently proposed golden-rule quantum transition state theory (GR-QTST) of Thapa, Fang, and Richardson [J. Chem. Phys. 150, 104107 (2019)]. While GR-QTST is not size consistent, leading to the possibility of unbounded errors in the rate, our modified method has no such issue and so can be reliably applied to condensed phase systems. Both methods involve path-integral sampling in a constrained ensemble; the two methods differ, however, in the choice of constraint functional. We demonstrate numerically that our modified method is as accurate as GR-QTST for the one-dimensional model considered by Thapa and co-workers. We then study a multidimensional spin-boson model, for which our method accurately predicts the true quantum rate, while GR-QTST breaks down with an increasing number of boson modes in the discretization of the spectral density. Our method is able to accurately predict reaction rates in the Marcus inverted regime without the need for the analytic continuation required by Wolynes theory.

Quantum mechanical spin dynamics of a molecular magnetoreceptor.

Radical pair recombination reactions are known to be sensitive to extremely weak magnetic fields and can therefore be said to function as molecular magnetoreceptors. The classic example is a carotenoid-porphyrin-fullerene (C•+PF•-) radical pair that has been shown to provide a "proof-of-principle" for the operation of a chemical compass [K. Maeda et al., Nature 453, 387 (2008)]. Previous simulations of this radical pair have employed semiclassical approximations, which are routinely applicable to its 47 coupled electronic and nuclear spins. However, calculating the exact quantum mechanical spin dynamics presents a significant challenge and has not been possible until now. Here, we use a recently developed method to perform numerically converged simulations of the C•+PF•- quantum mechanical spin dynamics, including all coupled spins. A comparison of these quantum mechanical simulations with various semiclassical approximations reveals that, while it is not perfect, the best semiclassical approximation does capture essentially all of the relevant physics in this problem.

Path integral methods for reaction rates in complex systems.

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. As a result of this progress, it is now routinely possible to calculate accurate rate constants including quantum mechanical zero point energy and tunnelling effects for arbitrarily complex (anharmonic and multi-dimensional) systems. This is true in the adiabatic (Born-Oppenheimer) limit, in the non-adiabatic (Fermi Golden Rule) limit, and everywhere between these two limits in the normal Marcus regime. Quantum mechanical effects on reaction rates can be enormous, even at room temperature, and the problem of including these effects in simulations of a wide variety of chemical reactions in complex systems has now effectively been solved.

How quantum is radical pair magnetoreception?

Currently the most likely mechanism of the magnetic compass sense in migratory songbirds relies on the coherent spin dynamics of pairs of photochemically formed radicals in the retina. Spin-conserving electron transfer reactions are thought to result in radical pairs whose near-degenerate electronic singlet and triplet states interconvert coherently as a result of hyperfine, exchange, and dipolar couplings and, crucially for a compass sensor, Zeeman interactions with the geomagnetic field. In this way, the yields of the reaction products can be influenced by magnetic interactions a million times smaller than kBT. The question we ask here is whether one can only account for the coherent spin dynamics using quantum mechanics. We find that semiclassical approximations to the spin dynamics of radical pairs only provide a satisfactory description of the anisotropic product yields when there is no electron spin-spin coupling, a situation unlikely to be consistent with a magnetic sensing function. Although these methods perform reasonably well for shorter-lived radical pairs with stronger electron-spin coupling, the accurate simulation of anisotropic magnetic field effects relevant to magnetoreception seems to require full quantum mechanical calculations.

Radical pair intersystem crossing: Quantum dynamics or incoherent kinetics?

Magnetic field effects on radical pair reactions arise due to the interplay of coherent electron spin dynamics and spin relaxation effects, a rigorous treatment of which requires the solution of the Liouville-von Neumann equation. However, it is often found that simple incoherent kinetic models of the radical pair singlet-triplet intersystem crossing provide an acceptable description of experimental measurements. In this paper, we outline the theoretical basis for this incoherent kinetic description, elucidating its connection to exact quantum mechanics. We show, in particular, how the finite lifetime of the radical pair spin states, as well as any additional spin-state dephasing, leads to incoherent intersystem crossing. We arrive at simple expressions for the radical pair spin state interconversion rates to which the functional form proposed recently by Steiner et al. [J. Phys. Chem. C 122, 11701 (2018)] can be regarded as an approximation. We also test the kinetic master equation against exact quantum dynamical simulations for a model radical pair and for a series of PTZ•+-Phn-PDI•- molecular wires.

An analysis of isomorphic RPMD in the golden rule limit.

We analyze the golden rule limit of the recently proposed isomorphic ring polymer (iso-RP) method. This method aims to combine an exact expression for the quantum mechanical partition function of a system with multiple electronic states with a pre-existing mixed quantum-classical (MQC) dynamics approximation, such as fewest switches surface hopping. Since the choice of the MQC method adds a degree of flexibility, we simplify the analysis by assuming that the dynamics used correctly reproduces the exact golden rule rate for a nonadiabatic (e.g., electron transfer) reaction in the high temperature limit. Having made this assumption, we obtain an expression for the iso-RP rate in the golden rule limit that is valid at any temperature. We then compare this rate with the exact rate for a series of simple spin-boson models. We find that the iso-RP method does not correctly predict how nuclear quantum effects affect the reaction rate in the golden rule limit. Most notably, it does not capture the quantum asymmetry in a conventional (Marcus) plot of the logarithm of the reaction rate against the thermodynamic driving force, and it also significantly overestimates the correct quantum mechanical golden rule rate for activationless electron transfer reactions. These results are analyzed and their implications discussed for the applicability of the iso-RP method to more general nonadiabatic reactions.

Electron spin relaxation in radical pairs: Beyond the Redfield approximation.

Relaxation processes can have a large effect on the spin selective electron transfer reactions of radical pairs. These processes are often treated using phenomenological relaxation superoperators or with some model for the microscopic relaxation mechanism treated within Bloch-Redfield-Wangsness theory. Here, we demonstrate that an alternative perturbative relaxation theory, based on the Nakajima-Zwanzig equation, has certain advantages over Redfield theory. In particular, the Nakajima-Zwanzig equation does not suffer from the severe positivity problem of Redfield theory in the static disorder limit. Combining the Nakajima-Zwanzig approach consistently with the Schulten-Wolynes semiclassical method, we obtain an efficient method for modeling the spin dynamics of radical pairs containing many hyperfine-coupled nuclear spins. This is then used to investigate the spin-dependent electron transfer reactions and intersystem crossing of dimethyljulolidine-naphthalenediimide (DMJ-NDI) radical ion pairs. By comparing our simulations with experimental data, we find evidence for a field-independent contribution to the triplet quantum yields of these reactions which cannot be explained by electron spin relaxation alone.

On the calculation of quantum mechanical electron transfer rates.

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 tunneling, zero point energy, and solvent friction on condensed phase electron transfer reactions. Comparison with exact hierarchical equations of motion 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 the nonadiabatic limit, and is therefore as accurate as any method could possibly be that attempts to generalize these methods to arbitrary electronic coupling strengths.

The quantum needle of the avian magnetic compass.

Migratory birds have a light-dependent magnetic compass, the mechanism of which is thought to involve radical pairs formed photochemically in cryptochrome proteins in the retina. Theoretical descriptions of this compass have thus far been unable to account for the high precision with which birds are able to detect the direction of the Earth's magnetic field. Here we use coherent spin dynamics simulations to explore the behavior of realistic models of cryptochrome-based radical pairs. We show that when the spin coherence persists for longer than a few microseconds, the output of the sensor contains a sharp feature, referred to as a spike. The spike arises from avoided crossings of the quantum mechanical spin energy-levels of radicals formed in cryptochromes. Such a feature could deliver a heading precision sufficient to explain the navigational behavior of migratory birds in the wild. Our results (i) afford new insights into radical pair magnetoreception, (ii) suggest ways in which the performance of the compass could have been optimized by evolution, (iii) may provide the beginnings of an explanation for the magnetic disorientation of migratory birds exposed to anthropogenic electromagnetic noise, and (iv) suggest that radical pair magnetoreception may be more of a quantum biology phenomenon than previously realized.

Simple and Accurate Method for Central Spin Problems.

We describe a simple quantum mechanical method that can be used to obtain accurate numerical results over long timescales for the spin correlation tensor of an electron spin that is hyperfine coupled to a large number of nuclear spins. This method does not suffer from the statistical errors that accompany a Monte Carlo sampling of the exact eigenstates of the central spin Hamiltonian obtained from the algebraic Bethe ansatz, or from the growth of the truncation error with time in the time-dependent density matrix renormalization group (TDMRG) approach. As a result, it can be applied to larger central spin problems than the algebraic Bethe ansatz, and for longer times than the TDMRG algorithm. It is therefore an ideal method to use to solve central spin problems, and we expect that it will also prove useful for a variety of related problems that arise in a number of different research fields.

Analytic continuation of Wolynes theory into the Marcus inverted regime.

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 tsp=iλspℏ, and provided λsp lies in the range -ß/2≤λsp≤ß/2, it is straightforward to evaluate the rate constant using information obtained from an imaginary time path integral calculation. However, if λsp 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 λsp<-ß/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.