Flux-correlation approach to characterizing reaction pathways in quantum systems: a study of condensed-phase proton-coupled electron transfer

Recent progress in approximate quantum dynamics methods for the study of proton-coupled electron transfer reactions

Proton-coupled electron transfer (PCET) reactions are ubiquitous natural processes at the heart of energy conversion reactions in photosynthesis and respiration, DNA repair, and diverse enzymatic reactions. Theoretical formulation and computational method developments have eyed modeling of thermal and photoinduced PCET for the last three decades. The accumulation of these studies, collected in dozens of reviews, accounts, and perspectives, has firmly established the influence of quantum effects, including non-adiabatic electronic transitions, vibrational relaxation, zero-point energy, and proton tunneling, on the rate and mechanism of PCET reactions. Here, we focus on some recently-developed methods, spanning the last eight years, that can quantitatively capture these effects in the PCET context and provide efficient means for their qualitative description in complex systems. The theoretical background of each method and their accuracy with respect to exact results are discussed and the results of relevant PCET simulations based on each method are presented.

Nonadiabatic Dynamics via the Symmetrical Quasi-Classical Method in the Presence of Anharmonicity

The symmetrical quasi-classical (SQC) method recently proposed by Miller and Cotton allows one to simulate nonadiabatic dynamics based on an algorithm with classical-like scaling with respect to system size. This is made possible by casting the electronic degrees of freedom in terms of mapping variables that can be propagated in a classical-like manner. While SQC was shown to be rather accurate when applied to benchmark models with harmonic electronic potential energy surfaces, it was also found to become inaccurate and to suffer numerical instabilities when applied to anharmonic systems. In this paper, we propose an extended SQC (E-SQC) methodology for overcoming those discrepancies by describing the anharmonic nuclear modes, which are coupled to the electronic degrees of freedom, in terms of classical-like mapping variables. The accuracy of E-SQC relative to standard SQC is demonstrated on benchmark models with quartic and Morse potential energy surfaces.

Photoexcited Ag ejection from a low-temperature He cluster: a simulation study by nonadiabatic Ehrenfest ring-polymer molecular dynamics

Nonadiabatic ring-polymer molecular dynamics simulations were performed to understand the photoexcitation dynamics of a low-temperature Ag·He₅₀₀ cluster.

Kinetically-constrained ring-polymer molecular dynamics for non-adiabatic chemistries involving solvent and donor–acceptor dynamical effects

We investigate the performance of the recently developed kinetically-constrained ring polymer molecular dynamics (KC-RPMD) method for the description of model condensed-phase electron transfer (ET) reactions in which solvent and donor–acceptor dynamics play an important role. Comparison of KC-RPMD with results from Golden-Rule rate theories and numerically exact quantum dynamics calculations demonstrates that KC-RPMD accurately captures the combination of electronic- and nuclear-dynamical effects throughout the Marcus (intermediate solvent friction) and Zusman (large solvent friction) regimes of ET. It is also demonstrated that KC-RPMD accurately describes systems in which the magnitude of the diabatic coupling depends on the position of a dynamical donor–acceptor mode. In addition to these successes, however, we present an unsurprising failure of KC-RPMD to capture the enhancement of the ET rate in the low solvent friction regime associated with nuclear coherence effects. In this analysis, we re-visit several aspects of the original KC-RPMD formulation, including the form of the kinetic constraint and the choice of the mass of the auxiliary electronic variable. In particular, we introduce a Langevin bath for the auxiliary electronic variable to correct for its unphysically low coupling with the nuclear degrees of freedom. This work demonstrates that the KC-RPMD method is well suited for the direct simulation of non-adiabatic donor–acceptor chemistries associated with many complex systems, including those for which solvent dynamics plays an important role in the reaction mechanism.

Microscopic molecular superfluid response: theory and simulations

Since its discovery in 1938, superfluidity has been the subject of much investigation because it provides a unique example of a macroscopic manifestation of quantum mechanics. About 60 years later, scientists successfully observed this phenomenon in the microscopic world though the spectroscopic Andronikashvili experiment in helium nano-droplets. This reduction of scale suggests that not only helium but also para-H2 (pH2) can be a candidate for superfluidity. This expectation is based on the fact that the smaller number of neighbours and surface effects of a finite-size cluster may hinder solidification and promote a liquid-like phase. The first prediction of superfluidity in pH2 clusters was reported in 1991 based on quantum Monte Carlo simulations. The possible superfluidity of pH2 was later indirectly observed in a spectroscopic Andronikashvili experiment in 2000. Since then, a growing number of studies have appeared, and theoretical simulations have been playing a special role because they help guide and interpret experiments. In this review, we go over the theoretical studies of pH2 superfluid clusters since the experiment of 2000. We provide a historical perspective and introduce the basic theoretical formalism along with key experimental advances. We then present illustrative results of the theoretical studies and comment on the possible future developments in the field. We include sufficient theoretical details such that the review can serve as a guide for newcomers to the field.