Unified expression for the rate constant of the bridged electron transfer derived by renormalization

Theoretical Modeling of Electrochemical Proton-Coupled Electron Transfer

Proton-coupled electron transfer (PCET) plays an essential role in a wide range of electrocatalytic processes. A vast array of theoretical and computational methods have been developed to study electrochemical PCET. These methods can be used to calculate redox potentials and pK_a values for molecular electrocatalysts, proton-coupled redox potentials and bond dissociation free energies for PCET at metal and semiconductor interfaces, and reorganization energies associated with electrochemical PCET. Periodic density functional theory can also be used to compute PCET activation energies and perform molecular dynamics simulations of electrochemical interfaces. Various approaches for maintaining a constant electrode potential in electronic structure calculations and modeling complex interactions in the electric double layer (EDL) have been developed. Theoretical formulations for both homogeneous and heterogeneous electrochemical PCET spanning the adiabatic, nonadiabatic, and solvent-controlled regimes have been developed and provide analytical expressions for the rate constants and current densities as functions of applied potential. The quantum mechanical treatment of the proton and inclusion of excited vibronic states have been shown to be critical for describing experimental data, such as Tafel slopes and potential-dependent kinetic isotope effects. The calculated rate constants can be used as input to microkinetic models and voltammogram simulations to elucidate complex electrocatalytic processes.

Photosynthetic reaction center variants made via genetic code expansion show Tyr at M210 tunes the initial electron transfer mechanism

Photosynthetic reaction centers (RCs) from Rhodobacter sphaeroides were engineered to vary the electronic properties of a key tyrosine (M210) close to an essential electron transfer component via its replacement with site-specific, genetically encoded noncanonical amino acid tyrosine analogs. High fidelity of noncanonical amino acid incorporation was verified with mass spectrometry and X-ray crystallography and demonstrated that RC variants exhibit no significant structural alterations relative to wild type (WT). Ultrafast transient absorption spectroscopy indicates the excited primary electron donor, P*, decays via a ∼4-ps and a ∼20-ps population to produce the charge-separated state P⁺H_A ^- in all variants. Global analysis indicates that in the ∼4-ps population, P⁺H_A ^- forms through a two-step process, P*→ P⁺B_A ^-→ P⁺H_A ^-, while in the ∼20-ps population, it forms via a one-step P* → P⁺H_A ^- superexchange mechanism. The percentage of the P* population that decays via the superexchange route varies from ∼25 to ∼45% among variants, while in WT, this percentage is ∼15%. Increases in the P* population that decays via superexchange correlate with increases in the free energy of the P⁺B_A ^- intermediate caused by a given M210 tyrosine analog. This was experimentally estimated through resonance Stark spectroscopy, redox titrations, and near-infrared absorption measurements. As the most energetically perturbative variant, 3-nitrotyrosine at M210 creates an ∼110-meV increase in the free energy of P⁺B_A ^- along with a dramatic diminution of the 1,030-nm transient absorption band indicative of P⁺B_A ^- formation. Collectively, this work indicates the tyrosine at M210 tunes the mechanism of primary electron transfer in the RC.

Charge Transfer through Redox Molecular Junctions in Nonequilibrated Solvents

Molecular conduction operating in dielectric solvent environments is often described using kinetic rates based on the Marcus theory of electron transfer at a molecule-metal electrode interface. However, the successive nature of charge transfer in such a system implies that the solvent does not necessarily reach equilibrium in such processes. Here we generalize the theory to account for solvent nonequilibrium and consider a molecular junction consisting of an electronic donor-acceptor system coupled to two metallic electrodes and placed in a polarizable solvent. We determine the nonequilbrium distribution of the solvent by solving diffusion equations in the strong- and weak-friction limits and calculate the charge current and its fluctuating behavior. In extreme limits, the absence of the solvent or fast solvent relaxation, the charge-transfer statistics is Poissonian, while it becomes correlated by the dynamic solvent between these limits. A Kramers-like turnover of the nonequilibrium current as a function of the solvent damping is found. Finally, we propose a way to tune the solvent-induced damping using geometrical control of the solvent dielectric response in nanostructured solvent channels.

Electron Transfer Methods in Open Systems

Utilization of electron transfer methods for description of quantum transport is popular due to simplicity of the formulation and its ability to account for basic physics of electron exchange between the system and baths. At the same time, the necessity to go beyond simple golden rule-type expressions for rates was indicated in the literature and ad hoc formulations were proposed. Similarly, kinetic schemes for quantum transport beyond the usual second-order Lindblad/Redfield considerations were discussed. Here we utilize recently introduced the nonequilibrium Hubbard Green's function diagrammatic technique to analyze the construction of rates in open systems. We show that previous considerations for rates of second and fourth order can be obtained as a particular case of zero- and second-order Green's function diagrammatic series with bare diagrams. We discuss limitations of previous considerations, stress advantages of the Hubbard Green's function approach in constructing the rates, and indicate that standard dressing of the diagrams is a natural way to account for additional baths/degrees of freedom in the formulation of generalized expressions for the rates.

Multistate treatments of the electronic coupling in donor-bridge-acceptor systems: insights and caveats from a simple model

We use a simple one-dimensional delta function electronic structure model (dfm) to investigate the results of a pair of multistate diabatization techniques (i.e., based on n states, with n ≥ 2) for linear DBA and DBBA (donor-bridge-acceptor) electron-transfer systems. In particular, we focus on the physical meaning of the couplings obtained from multistate methods and their relationship to two-state (n = 2) coupling elements. On the basis of the simple dfm approach, which allows exact as well as finite basis set treatment and has no many-electron effects, we conclude that for orthogonal diabatic states, it is difficult to assign clear physical significance to multistate matrix elements for coupling beyond nearest-neighbor contacts. The implications of these results for more complex multistate many-electron treatments are discussed. It is emphasized that physically meaningful coupling elements must involve states that are orthogonal, either explicitly or implicitly.