Proton-coupled electron transfer reactions in solution: Molecular dynamics with quantum transitions for model systems

Different flavors of exact-factorization-based mixed quantum-classical methods for multistate dynamics

The exact factorization approach has led to the development of new mixed quantum-classical methods for simulating coupled electron-ion dynamics. We compare their performance for dynamics when more than two electronic states are occupied at a given time, and analyze: (1) the use of coupled versus auxiliary trajectories in evaluating the electron-nuclear correlation terms, (2) the approximation of using these terms within surface-hopping and Ehrenfest frameworks, and (3) the relevance of the exact conditions of zero population transfer away from nonadiabatic coupling regions and total energy conservation. Dynamics through the three-state conical intersection in the uracil radical cation as well as polaritonic models in one dimension are studied.

DECIDE: A Deterministic Mixed Quantum-Classical Dynamics Approach

Mixed quantum-classical dynamics provides an efficient way of simulating the dynamics of quantum subsystems coupled to many-body environments. Many processes, including proton-transfer reactions, electron-transfer reactions, and vibrational energy transport, for example, take place in such open systems. The most accurate algorithms for performing mixed quantum-classical simulations require very large ensembles of trajectories to obtain converged expectation values, which is computationally prohibitive for quantum subsystems containing even a few degrees of freedom. The recently developed “Deterministic evolution of coordinates with initial decoupled equations” (DECIDE) method has demonstrated high accuracy and low computational cost for a host of model systems; however, these applications relied on representing the equations of motion in subsystem and adiabatic energy bases. While these representations are convenient for certain systems, the position representation is convenient for many other systems, including systems undergoing proton- and electron-transfer reactions. Thus, in this review, we provide a step-by-step derivation of the DECIDE approach and demonstrate how to cast the DECIDE equations in a quantum harmonic oscillator position basis for a simple one-dimensional (1D) hydrogen bond model. After integrating the DECIDE equations of motion on this basis, we show that the total energy of the system is conserved for this model and calculate various quantities of interest. Limitations of casting the equations in an incomplete basis are also discussed.

Benchmarking the Surface Hopping Method to Include Nuclear Quantum Effects

We have benchmarked the surface hopping method to capture nuclear quantum effects in the spin-Boson model in the deep tunneling regime. The thermal populations and the rate constants calculated using the surface hopping method are compared with those calculated using Boltzmann theory and Fermi's golden rule, respectively. Additionally, we have proposed a simple kinetic model that partially includes nuclear quantum effects within Marcus theory, and the results of the surface hopping method are analyzed under the framework of this simple kinetic model. A broad range of parameters are investigated to identify the regimes for the successes and failures of the surface hopping method. This work shows that with the accurate treatment of decoherence and velocity reversal, surface hopping can generally capture the nuclear quantum effects in the deep tunneling and weak diabatic coupling regime.

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.

Case studies of the time-dependent potential energy surface for dynamics in cavities

The exact time-dependent potential energy surface driving the nuclear dynamics was recently shown to be a useful tool to understand and interpret the coupling of nuclei, electrons, and photons in cavity settings. Here, we provide a detailed analysis of its structure for exactly solvable systems that model two phenomena: cavity-induced suppression of proton-coupled electron-transfer and its dependence on the initial state, and cavity-induced electronic excitation. We demonstrate the inadequacy of simply using a weighted average of polaritonic surfaces to determine the dynamics. Such a weighted average misses a crucial term that redistributes energy between the nuclear and the polaritonic systems, and this term can in fact become a predominant term in determining the nuclear dynamics when several polaritonic surfaces are involved. Evolving an ensemble of classical trajectories on the exact potential energy surface reproduces the nuclear wavepacket quite accurately, while evolving on the weighted polaritonic surface fails after a short period of time. The implications and prospects for application of mixed quantum-classical methods based on this surface are discussed.

Effect of many modes on self-polarization and photochemical suppression in cavities

The standard description of cavity-modified molecular reactions typically involves a single (resonant) mode, while in reality, the quantum cavity supports a range of photon modes. Here, we demonstrate that as more photon modes are accounted for, physicochemical phenomena can dramatically change, as illustrated by the cavity-induced suppression of the important and ubiquitous process of proton-coupled electron-transfer. Using a multi-trajectory Ehrenfest treatment for the photon-modes, we find that self-polarization effects become essential, and we introduce the concept of self-polarization-modified Born-Oppenheimer surfaces as a new construct to analyze dynamics. As the number of cavity photon modes increases, the increasing deviation of these surfaces from the cavity-free Born-Oppenheimer surfaces, together with the interplay between photon emission and absorption inside the widening bands of these surfaces, leads to enhanced suppression. The present findings are general and will have implications for the description and control of cavity-driven physical processes of molecules, nanostructures, and solids embedded in cavities.

Exact Potential Energy Surface for Molecules in Cavities

We find and analyze the exact time-dependent potential energy surface driving the proton motion for a model of cavity-induced suppression of proton-coupled electron transfer. We show how, in contrast to the polaritonic surfaces, its features directly correlate to the proton dynamics and we discuss cavity modifications of its structure responsible for the suppression. The results highlight the interplay between nonadiabatic effects from coupling to photons and coupling to electrons and suggest caution is needed when applying traditional dynamics methods based on polaritonic surfaces.

Perspective: How to understand electronic friction

Electronic friction is a correction to the Born-Oppenheimer approximation, whereby nuclei in motion experience a drag in the presence of a manifold of electronic states. The notion of electronic friction has a long history and has been (re-)discovered in the context of a wide variety of different chemical and physical systems including, but not limited to, surface scattering events, surface reactions or chemisorption, electrochemistry, and conduction through molecular-(or nano-) junctions. Over the years, quite a few different forms of electronic friction have been offered in the literature. In this perspective, we briefly review these developments of electronic friction, highlighting the fact that we can now isolate a single, unifying form for (Markovian) electronic friction. We also focus on the role of electron-electron interactions for understanding frictional effects and offer our thoughts on the strengths and weaknesses of using electronic friction to model dynamics in general.

Perspective: Reaches of chemical physics in biology

Chemical physics as a discipline contributes many experimental tools, algorithms, and fundamental theoretical models that can be applied to biological problems. This is especially true now as the molecular level and the systems level descriptions begin to connect, and multi-scale approaches are being developed to solve cutting edge problems in biology. In some cases, the concepts and tools got their start in non-biological fields, and migrated over, such as the idea of glassy landscapes, fluorescence spectroscopy, or master equation approaches. In other cases, the tools were specifically developed with biological physics applications in mind, such as modeling of single molecule trajectories or super-resolution laser techniques. In this introduction to the special topic section on chemical physics of biological systems, we consider a wide range of contributions, all the way from the molecular level, to molecular assemblies, chemical physics of the cell, and finally systems-level approaches, based on the contributions to this special issue. Chemical physicists can look forward to an exciting future where computational tools, analytical models, and new instrumentation will push the boundaries of biological inquiry.

Proton-coupled electron transfer in biology: results from synergistic studies in natural and model systems

Proton-coupled electron transfer (PCET) underpins energy conversion in biology. PCET may occur with the unidirectional or bidirectional transfer of a proton and electron and may proceed synchronously or asynchronously. To illustrate the role of PCET in biology, this review presents complementary biological and model systems that explore PCET in electron transfer (ET) through hydrogen bonds [azurin as compared to donor-acceptor (D-A) hydrogen-bonded networks], the activation of C-H bonds [alcohol dehydrogenase and soybean lipoxygenase (SLO) as compared to Fe(III) metal complexes], and the generation and transport of amino acid radicals [photosystem II (PSII) and ribonucleotide reductase (RNR) as compared to tyrosine-modified photoactive Re(I) and Ru(II) complexes]. In providing these comparisons, the fundamental principles of PCET in biology are illustrated in a tangible way.

Charge-dependent model for many-body polarization, exchange, and dispersion interactions in hybrid quantum mechanical/molecular mechanical calculations

This work explores a new charge-dependent energy model consisting of van der Waals and polarization interactions between the quantum mechanical (QM) and molecular mechanical (MM) regions in a combined QMMM calculation. van der Waals interactions are commonly treated using empirical Lennard-Jones potentials, whose parameters are often chosen based on the QM atom type (e.g., based on hybridization or specific covalent bonding environment). This strategy for determination of QMMM nonbonding interactions becomes tedious to parametrize and lacks robust transferability. Problems occur in the study of chemical reactions where the "atom type" is a complex function of the reaction coordinate. This is particularly problematic for reactions, where atoms or localized functional groups undergo changes in charge state and hybridization. In the present work we propose a new model for nonelectrostatic nonbonded interactions in QMMM calculations that overcomes many of these problems. The model is based on a scaled overlap model for repulsive exchange and attractive dispersion interactions that is a function of atomic charge. The model is chemically significant since it properly correlates atomic size, softness, polarizability, and dispersion terms with minimal one-body parameters that are functions of the atomic charge. Tests of the model are examined for rare-gas interactions with neutral and charged atoms in order to demonstrate improved transferability. The present work provides a new framework for modeling QMMM interactions with improved accuracy and transferability.

Nonadiabatic dynamics of charge transfer in diatomic anion clusters

We have studied the photodissociation and recombination dynamics of the diatomic anions X(2)(-) and XY(-) designed to mimic I(2)(-) and ICl(-), respectively, by using a one-electron model in size-selected N(2)O clusters. The one-electron model is composed of two nuclei and an extra electron moving in a two-dimensional plane including the two nuclei. The main purpose of this study is to explain the salient features of various dynamical processes of molecular ions in clusters using a simple theoretical model. For heteronuclear diatomic anions, a mass disparity and asymmetric electron affinity between the X and Y atoms lead to different phenomena from the homonuclear case. The XY(-) anion shows efficient recombination for a smaller cluster size due to the effect of collision-mediated energy transfer and an inherent potential wall on excited state at asymptotic region, while the recombination for the X(2)(-) anion is due to rearrangement of solvent configuration and faster nonadiabatic transitions. The results of the present study illustrate the microscopic details of the electronically nonadiabatic processes which control the photodissociation dynamics of molecular ions in clusters.

The Possible Role of Proton-Coupled Electron Transfer (PCET) in Water Oxidation by Photosystem II

All higher life forms use oxygen and respiration as their primary energy source. The oxygen comes from water by solar-energy conversion in photosynthetic membranes. In green plants, light absorption in photosystem II (PSII) drives electron-transfer activation of the oxygen-evolving complex (OEC). The mechanism of water oxidation by the OEC has long been a subject of great interest to biologists and chemists. With the availability of new molecular-level protein structures from X-ray crystallography and EXAFS, as well as the accumulated results from numerous experiments and theoretical studies, it is possible to suggest how water may be oxidized at the OEC. An integrated sequence of light-driven reactions that exploit coupled electron-proton transfer (EPT) could be the key to water oxidation. When these reactions are combined with long-range proton transfer (by sequential local proton transfers), it may be possible to view the OEC as an intricate structure that is "wired for protons".

Quantum-classical dynamics of wave fields

An approach to the quantum-classical mechanics of phase space dependent operators, which has been proposed recently, is remodeled as a formalism for wave fields. Such wave fields obey a system of coupled nonlinear equations that can be written by means of a suitable non-Hamiltonian bracket. As an example, the theory is applied to the relaxation dynamics of the spin-boson model. In the adiabatic limit, a good agreement with calculations performed by the operator approach is obtained. Moreover, the theory proposed in this paper can take nonadiabatic effects into account without resorting to surface-hopping approximations. Hence, the results obtained follow qualitatively those of previous surface-hopping calculations and increase by a factor of (at least) 2, the time length over which nonadiabatic dynamics can be propagated with small statistical errors. Moreover, it is worth to note that the dynamics of quantum-classical wave fields proposed here is a straightforward non-Hamiltonian generalization of the formalism for nonlinear quantum mechanics that Weinberg introduced recently.