Ab Initio Molecular Dynamics Simulation of the Aqueous Ru2+/Ru3+ Redox Reaction: The Marcus Perspective

Redox Properties of Flavin in BLUF and LOV Photoreceptor Proteins from Hybrid QM/MM Molecular Dynamics Simulation

Flavins play an important role in many oxidation and reduction processes in biological systems. For example, flavin adenine dinucleotide (FAD) and flavin mononucleotide (FMN) are common cofactors found in enzymatic proteins that use the special redox properties of these flavin molecules for their catalytic or photoactive functions. The redox potential of the flavin is strongly affected by its (protein) environment; however, the underlying molecular interactions of this effect are still unknown. Using hybrid quantum mechanics/molecular mechanics (QM/MM) simulation techniques, we have studied the redox properties of flavin in the gas phase, aqueous solution, and two different protein environments, in particular, a BLUF and a LOV photoreceptor domain. By mapping the changes in electrostatic potential and solvent structure, we gain insight into how specific polarization of the flavin by its environment tunes the reduction potential. We find also that accurate calculation of the reduction potentials of these systems by using the hybrid QM/MM approach is hampered by a too limited sampling of the counterion configurations and by artifacts at the QM/MM boundary. We make suggestions for how these issues can be overcome.

Electron transfer reaction of TEMPO-based organic radical batteries in different solvent environments: comparing quantum and classical approaches

In this study, we delve into the complex electron transfer reactions associated with the redox-active (2,2,6,6-tetramethylpiperidin-1-yl)oxyl (TEMPO), a common component in organic radical batteries (ORBs). Our approach estimates quantum electron-transfer (ET) energies using Density Functional Theory (DFT) calculations by sampling from structures simulated classically. This work presents a comparative study of reorganization energies in ET reactions across different solvents. Furthermore, we investigate how changes in the electrolyte environment can modify the reorganization energy and, consequently, impact ET dynamics. We also explore the relationship between classical and quantum vertical energies using linear regression models. Importantly, this comparison between quantum and classical vertical energies underscores the role of quantum effects, like charge delocalization, in offering added stabilization post-redox reactions. These effects are not adequately represented by the classical vertical energy distribution. Our study shows that, although we find a significant correlation between the vertical energies computed by DFT and the classical force field, the regression parameters depend on the solvent, highlighting that classical methods should be benchmarked by DFT before applying them to novel electrolyte materials.

Hybrid Functional and Plane Waves based Ab Initio Molecular Dynamics Study of the Aqueous Fe2+ /Fe3+ Redox Reaction

Kohn-Sham density functional theory and plane wave basis set based ab initio molecular dynamics (AIMD) simulation is a powerful tool for studying complex reactions in solutions, such as electron transfer (ET) reactions involving Fe²⁺ /Fe³⁺ ions in water. In most cases, such simulations are performed using density functionals at the level of Generalized Gradient Approximation (GGA). The challenge in modelling ET reactions is the poor quality of GGA functionals in predicting properties of such open-shell systems due to the inevitable self-interaction error (SIE). While hybrid functionals can minimize SIE, standard plane-wave based AIMD at that level of theory is typically 150 times slower than GGA for systems containing ∼100 atoms. Among several approaches reported to speed-up AIMD simulations with hybrid functionals, the noise-stabilized MD (NSMD) procedure, together with the use of localized orbitals to compute the required exchange integrals, is an attractive option. In this work, we demonstrate the application of the NSMD approach for studying the Fe²⁺ /Fe³⁺ redox reaction in water. It is shown here that long AIMD trajectories at the level of hybrid density functionals can be obtained using this approach. Redox properties of the aqueous Fe²⁺ /Fe³⁺ system computed from these simulations are compared with the available experimental data for validation.

Implementation and Validation of Constrained Density Functional Theory Forces in the CP2K Package

Constrained density functional theory (CDFT) is a powerful tool for the prediction of electron transfer parameters in condensed phase simulations at a reasonable computational cost. In this work we present an extension to CDFT in the popular mixed Gaussian/plane wave electronic structure package CP2K, implementing the additional force terms arising from a constraint based on Hirshfeld charge partitioning. This improves upon the existing Becke partitioning scheme, which is prone to give unphysical atomic charges. We verify this implementation for a variety of systems: electron transfer in (H₂O)₂⁺ in a vacuum, electron tunnelling between oxygen vacancy centers in solid MgO, and electron self-exchange in aqueous Ru²⁺–Ru³⁺. We find good agreement with previous plane-wave CDFT results for the same systems, but at a significantly lower computational cost, and we discuss the general reliability of condensed phase CDFT calculations.

Probing the thermodynamics and kinetics of ethylene carbonate reduction at the electrode-electrolyte interface with molecular simulations

Understanding the formation of the solid-electrolyte interphase (SEI) in lithium-ion batteries is an ongoing area of research due to its high degree of complexity and the difficulties encountered by experimental studies. Herein, we investigate the initial stage of SEI growth, the reduction reaction of ethylene carbonate (EC), from both a thermodynamic and a kinetic approach with theory and molecular simulations. We employed both the potential distribution theorem and the Solvation Method based on Density (SMD) to EC solvation for the estimation of reduction potentials of Li⁺, EC, and Li⁺-solvating EC (s-EC) as well as reduction rate constants of EC and s-EC. We find that solvation effects greatly influence these quantities of interest, particularly the Li⁺/Li reference electrode potential in EC solvent. Furthermore, we also compute the inner- and outer-sphere reorganization energies for both EC and s-EC at the interface of liquid EC and a hydroxyl-terminated graphite surface, where total reorganization energies are predicted to be 76.6 and 88.9 kcal/mol, respectively. With the computed reorganization energies, we estimate reduction rate constants across a range of overpotentials and show that EC has a larger electron transfer rate constant than s-EC at equilibrium, despite s-EC being more thermodynamically favorable. Overall, this manuscript demonstrates how ion solvation effects largely govern the prediction of reduction potentials and electron transfer rate constants at the electrode-electrolyte interface.

The Influence of Spectator Cations on Solvent Reorganization Energy Is a Short-Range Effect

In this manuscript, we use classical molecular dynamics simulation to explore the origin of specific cation effects on the rates of bulk-phase aqueous electron transfer (ET) reactions. We consider 0.6 M solutions of Cl^- and a series of different cations: Li⁺, Na⁺, K⁺, Rb⁺, and Cs⁺. We evaluate the collective electrostatic fluctuations that drive Marcus-like ET and find that they are essentially unaffected by changes in the cationic species. This finding implies that the structure making/breaking properties of various cations do not exert a significant influence on bulk-phase ET reactions. We evaluate the role of ion pairing in these specific cation effects and find, unsurprisingly, that model redox anions that are more highly charged tend to pair more effectively with spectator cations than their monovalent counterparts. We demonstrate that this ion pairing significantly affects local electrostatic fluctuations for the anionic redox species and thus conclude that ion pairing is one of the likely sources of rate-dependent cation effects in aqueous ET reactions.

Frontiers in molecular simulation of solvated ions, molecules and interfaces

This themed collection is a collection of articles on frontiers in molecular simulation of solvated ions, molecules and interfaces.

Microscopic Picture of the Solvent Reorganization During Electron Transfer to Flavin in Water

The redox potential of molecular species is largely modulated by its molecular environment so that a change of the environment will lead to a different redox potential. However, a detailed molecular picture of reorganization of the environment upon reduction is still unclear. To unravel the details of the solvent reorganization during electron transfer, we have performed density functional theory-based molecular dynamics (DFT-MD) and hybrid quantum mechanics/molecular mechanics (QM/MM) simulations of the reduction of lumiflavin. Previously, we have calculated the reduction free energy curves of the redox half reactions of lumiflavin in water as a function of the instantaneous gap energy (ΔE) ( J. Chem. Theory Comput. 2013 , 9 , 3889 - 3899 ). In this work, we focus on finding the changes in the solvent environment that correlate with this ΔE reaction coordinate. Comparing the QM/MM simulations, in which the solvent is modeled with an empirical force field, with the (full) DFT-MD simulations, we find that the response through electronic polarization plays a significant role in the latter case. Also a small charge transfer between flavin and solvent is observed in the full DFT treatment. As a result, we find only in the case of the QM/MM model a strong correlation between ΔE and the (pairwise computed) electrostatic potential (ESP) at the flavin due to the solvent. By analyzing the contribution of the ESP at the flavin per solvent molecule, we cannot only distinguish between the different modes of hydration by solvent molecules that coordinate at the hydrophilic and hydrophobic sides of the flavin molecule but also quantify their contribution to the reorganization free energy by measuring the ESP fluctuations per solvent molecule.

Electron transfer in extended systems: characterization by periodic density functional theory including the electronic coupling

We describe a new computer implementation of electron transfer (ET) theory in extended systems treated by periodic density functional theory (DFT), including the calculation of the electronic coupling transition element V_AB. In particular, the development opens up the full characterization of electron transfer in the solid state. The approach is valid for any single-determinant wavefunction with localized character representing the electronic structure of the system, from Hartree-Fock (HF) theory, to density functional theory (DFT), hybrid DFT theory, DFT+U theory, and constrained DFT (cDFT) theory. The implementation in CP2K reuses the high-performance functions of the code. The computational cost is equivalent to only one iteration of an HF calculation. We present test calculations for electron transfer in a number of systems, including a 1D-model of ferric oxide, hematite Fe₂O₃, rutile TiO₂, and finally bismuth vanadate BiVO₄.

Gibbs free-energy differences between polymorphs via a diabat approach

Polymorph free-energy differences are critical to several applications. A recently proposed diabat interpolation framework estimated free-energy differences between polymorphs by quadratic interpolation of diabats. This work extends the Zwanzig-Bennett relation to the NPT ensemble so that the diabats directly give Gibbs free-energy differences. We also demonstrate how the approach can be used in cases where the diabats are not parabolic. We illustrate the diabat method for Gibbs free-energy difference of zirconium (BCC and HCP phases) and compare it with the conventional lattice switch Monte Carlo approach.

Quantum modeling of ultrafast photoinduced charge separation

Phenomena involving electron transfer are ubiquitous in nature, photosynthesis and enzymes or protein activity being prominent examples. Their deep understanding thus represents a mandatory scientific goal. Moreover, controlling the separation of photogenerated charges is a crucial prerequisite in many applicative contexts, including quantum electronics, photo-electrochemical water splitting, photocatalytic dye degradation, and energy conversion. In particular, photoinduced charge separation is the pivotal step driving the storage of sun light into electrical or chemical energy. If properly mastered, these processes may also allow us to achieve a better command of information storage at the nanoscale, as required for the development of molecular electronics, optical switching, or quantum technologies, amongst others. In this Topical Review we survey recent progress in the understanding of ultrafast charge separation from photoexcited states. We report the state-of-the-art of the observation and theoretical description of charge separation phenomena in the ultrafast regime mainly focusing on molecular- and nano-sized solar energy conversion systems. In particular, we examine different proposed mechanisms driving ultrafast charge dynamics, with particular regard to the role of quantum coherence and electron-nuclear coupling, and link experimental observations to theoretical approaches based either on model Hamiltonians or on first principles simulations.

Reactive trajectories of the Ru2+/3+ self-exchange reaction and the connection to Marcus' theory

Outer sphere electron transfer between two ions in aqueous solution is a rare event on the time scale of first principles molecular dynamics simulations. We have used transition path sampling to generate an ensemble of reactive trajectories of the self-exchange reaction between a pair of Ru²⁺ and Ru³⁺ ions in water. To distinguish between the reactant and product states, we use as an order parameter the position of the maximally localised Wannier center associated with the transferring electron. This allows us to align the trajectories with respect to the moment of barrier crossing and compute statistical averages over the path ensemble. We compare our order parameter with two typical reaction coordinates used in applications of Marcus theory of electron transfer: the vertical gap energy and the solvent electrostatic potential at the ions.

A method to calculate redox potentials relative to the normal hydrogen electrode in nonaqueous solution by using density functional theory-based molecular dynamics

We demonstrate the redox potential calculation relative to the normal hydrogen electrode (NHE) in nonaqueous solution using a density functional theory-based molecular dynamics (DFT-MD) simulation. The calculation of the NHE in nonaqueous solution consists of two processes: the first step is the equilibrated simulation for a proton in nonaqueous solution to determine the space for inserting a proton in solution, and the second step is the thermodynamic integration method to calculate the solvation energy of the proton in the nonaqueous solution. In this work, we apply the method for a cation and an anion, i.e., copper(ii)/copper(i) and iodine/iodide in acetonitrile solution, and show that the errors in the calculated redox potential from experiments are within 0.21 V.

Computational electrochemistry: prediction of liquid-phase reduction potentials

This article reviews recent developments and applications in the area of computational electrochemistry. Our focus is on predicting the reduction potentials of electron transfer and other electrochemical reactions and half-reactions in both aqueous and nonaqueous solutions. Topics covered include various computational protocols that combine quantum mechanical electronic structure methods (such as density functional theory) with implicit-solvent models, explicit-solvent protocols that employ Monte Carlo or molecular dynamics simulations (for example, Car-Parrinello molecular dynamics using the grand canonical ensemble formalism), and the Marcus theory of electronic charge transfer. We also review computational approaches based on empirical relationships between molecular and electronic structure and electron transfer reactivity. The scope of the implicit-solvent protocols is emphasized, and the present status of the theory and future directions are outlined.

Self-assembled monolayers of C60-triphenylamine dyads as photo-switched interfacial layers for potential application in photovoltaic cells

C60-Triphenylamine dyads were synthesized for incorporation as photoswitched interfacial layers in organic photovoltaic (OPV) cells. Self-assembled monolayers (SAMs) of these dyads on gold (through S-Au and C60-Au interactions) were prepared through one or two adsorption processes, and their packing densities were fully characterized. Analysis using quartz crystal microbalance (QCM) and electrochemical impedance spectroscopy (EIS) measurements indicated that all SAMs exhibit dense coverage on the gold surfaces. Electrochemical desorption in KOH confirmed that the cis-1 dyad is anchored to the gold surface through its thiol group. Impedance measurements in the absence and presence of UV irradiation were performed to observe the photoswitched properties of these surface confined dyads. Upon UV light exposure of the SAMs, the charge-transfer resistance decreased when Fe(CN)6(3-/4-) was used as the probe redox couple and increased with Ru(NH3)6(3+/2+), confirming the generation of positive charges on the surface upon UV irradiation.

Fragment-based Quantum Mechanical/Molecular Mechanical Simulations of Thermodynamic and Kinetic Process of the Ru2+-Ru3+ Self-Exchange Electron Transfer

A fragment-based fractional number of electron (FNE) approach, is developed to study entire electron transfer (ET) processes from the electron donor region to the acceptor region in condensed phase. Both regions are described by the density-fragment interaction (DFI) method while FNE as an efficient ET order parameter is applied to simulate the electron transfer process. In association with the QM/MM energy expression, the DFI-FNE method is demonstrated to describe ET processes robustly with the Ru2+-Ru3+ self-exchange ET as a proof-of-concept example. This method allows for systematic calculations of redox free energies, reorganization energies, and electronic couplings, and the absolute ET rate constants within the Marcus regime.

Charge transfer in model peptides: obtaining Marcus parameters from molecular simulation

Charge transfer within and between biomolecules remains a highly active field of biophysics. Due to the complexities of real systems, model compounds are a useful alternative to study the mechanistic fundamentals of charge transfer. In recent years, such model experiments have been underpinned by molecular simulation methods as well. In this work, we study electron hole transfer in helical model peptides by means of molecular dynamics simulations. A theoretical framework to extract Marcus parameters of charge transfer from simulations is presented. We find that the peptides form stable helical structures with sequence dependent small deviations from ideal PPII helices. We identify direct exposure of charged side chains to solvent as a cause of high reorganization energies, significantly larger than typical for electron transfer in proteins. This, together with small direct couplings, makes long-range superexchange electron transport in this system very slow. In good agreement with experiment, direct transfer between the terminal amino acid side chains can be dicounted in favor of a two-step hopping process if appropriate bridging groups exist.

Real-world predictions from ab initio molecular dynamics simulations

In this review we present the techniques of ab initio molecular dynamics simulation improved to its current stage where the analysis of existing processes and the prediction of further chemical features and real-world processes are feasible. For this reason we describe the relevant developments in ab initio molecular dynamics leading to this stage. Among them, parallel implementations, different basis set functions, density functionals, and van der Waals corrections are reported. The chemical features accessible through AIMD are discussed. These are IR, NMR, as well as EXAFS spectra, sampling methods like metadynamics and others, Wannier functions, dipole moments of molecules in condensed phase, and many other properties. Electrochemical reactions investigated by ab initio molecular dynamics methods in solution, on surfaces as well as complex interfaces, are also presented.

Redox potentials and pKa for benzoquinone from density functional theory based molecular dynamics

The density functional theory based molecular dynamics (DFTMD) method for the computation of redox free energies presented in previous publications and the more recent modification for computation of acidity constants are reviewed. The method uses a half reaction scheme based on reversible insertion/removal of electrons and protons. The proton insertion is assisted by restraining potentials acting as chaperones. The procedure for relating the calculated deprotonation free energies to Brønsted acidities (pK(a)) and the oxidation free energies to electrode potentials with respect to the normal hydrogen electrode is discussed in some detail. The method is validated in an application to the reduction of aqueous 1,4-benzoquinone. The conversion of hydroquinone to quinone can take place via a number of alternative pathways consisting of combinations of acid dissociations, oxidations, or dehydrogenations. The free energy changes of all elementary steps (ten in total) are computed. The accuracy of the calculations is assessed by comparing the energies of different pathways for the same reaction (Hess's law) and by comparison to experiment. This two-sided test enables us to separate the errors related with the restrictions on length and time scales accessible to DFTMD from the errors introduced by the DFT approximation. It is found that the DFT approximation is the main source of error for oxidation free energies.

The diabatic picture of electron transfer, reaction barriers, and molecular dynamics

Diabatic states have a long history in chemistry, beginning with early valence bond pictures of molecular bonding and extending through the construction of model potential energy surfaces to the modern proliferation of methods for computing these elusive states. In this review, we summarize the basic principles that define the diabatic basis and demonstrate how they can be applied in the specific context of constrained density functional theory. Using illustrative examples from electron transfer and chemical reactions, we show how the diabatic picture can be used to extract qualitative insight and quantitative predictions about energy landscapes. The review closes with a brief summary of the challenges and prospects for the further application of diabatic states in chemistry.

Charge constrained density functional molecular dynamics for simulation of condensed phase electron transfer reactions

We present a plane-wave basis set implementation of charge constrained density functional molecular dynamics (CDFT-MD) for simulation of electron transfer reactions in condensed phase systems. Following the earlier work of Wu and Van Voorhis [Phys. Rev. A 72, 024502 (2005)], the density functional is minimized under the constraint that the charge difference between donor and acceptor is equal to a given value. The classical ion dynamics is propagated on the Born-Oppenheimer surface of the charge constrained state. We investigate the dependence of the constrained energy and of the energy gap on the definition of the charge and present expressions for the constraint forces. The method is applied to the Ru2+-Ru3+ electron self-exchange reaction in aqueous solution. Sampling the vertical energy gap along CDFT-MD trajectories and correcting for finite size effects, a reorganization free energy of 1.6 eV is obtained. This is 0.1-0.2 eV lower than a previous estimate based on a continuum model for solvation. The smaller value for the reorganization free energy can be explained by the fact that the Ru-O distances of the divalent and trivalent Ru hexahydrates are predicted to be more similar in the electron transfer complex than for the separated aqua ions.

Calculating solution redox free energies with ab initio quantum mechanical/molecular mechanical minimum free energy path method

A quantum mechanical/molecular mechanical minimum free energy path (QM/MM-MFEP) method was developed to calculate the redox free energies of large systems in solution with greatly enhanced efficiency for conformation sampling. The QM/MM-MFEP method describes the thermodynamics of a system on the potential of mean force surface of the solute degrees of freedom. The molecular dynamics (MD) sampling is only carried out with the QM subsystem fixed. It thus avoids "on-the-fly" QM calculations and thus overcomes the high computational cost in the direct QM/MM MD sampling. In the applications to two metal complexes in aqueous solution, the new QM/MM-MFEP method yielded redox free energies in good agreement with those calculated from the direct QM/MM MD method. Two larger biologically important redox molecules, lumichrome and riboflavin, were further investigated to demonstrate the efficiency of the method. The enhanced efficiency and uncompromised accuracy are especially significant for biochemical systems. The QM/MM-MFEP method thus provides an efficient approach to free energy simulation of complex electron transfer reactions.

Water at an electrochemical interface--a simulation study

The results of molecular dynamics simulations of the properties of water in an aqueous ionic solution close to an interface with a model metallic electrode are described. In the simulations the electrode behaves as an ideally polarizable hydrophilic metal, supporting image-charge interactions with charged species, and it is maintained at a constant electrical potential with respect to the solution so that the model is a textbook representation of an electrochemical interface through which no current is passing. We show how water is strongly attracted to and ordered at the electrode surface. This ordering is different to the structure that might be imagined from continuum models of electrode interfaces. Further, this ordering significantly affects the probability of ions reaching the surface. We describe the concomitant motion and configurations of the water and ions as functions of the electrode potential, and we analyze the length scales over which ionic atmospheres fluctuate. The statistics of these fluctuations depend upon surface structure and ionic strength. The fluctuations are large--sufficiently so that the mean ionic atmosphere is a poor descriptor of the aqueous environment near a metal surface. The importance of this finding for a description of electrochemical reactions is examined by calculating, directly from the simulation, Marcus free-energy profiles for transfer of charge between the electrode and a redox species in the solution and comparing the results with the predictions of continuum theories. Significant departures from the electrochemical textbook descriptions of the phenomenon are found and their physical origins are characterized from the atomistic perspective of the simulations.

Development and application of ab initio QM/MM methods for mechanistic simulation of reactions in solution and in enzymes

Determining the free energies and mechanisms of chemical reactions in solution and enzymes is a major challenge. For such complex reaction processes, combined quantum mechanics/molecular mechanics (QM/MM) method is the most effective simulation method to provide an accurate and efficient theoretical description of the molecular system. The computational costs of ab initio QM methods, however, have limited the application of ab initio QM/MM methods. Recent advances in ab initio QM/MM methods allowed the accurate simulation of the free energies for reactions in solution and in enzymes and thus paved the way for broader application of the ab initio QM/MM methods. We review here the theoretical developments and applications of the ab initio QM/MM methods, focusing on the determination of reaction path and the free energies of the reaction processes in solution and enzymes.

Electrostatic fluctuations in cavities within polar liquids and thermodynamics of polar solvation

We present the results of numerical simulations of the fluctuations of the electrostatic potential and electric field inside cavities created in the fluid of dipolar hard spheres. We found that the thermodynamics of polar solvation changes dramatically when the cavity size becomes about 4-5 times larger than the size of the liquid particle. The range of small cavities can be reasonably understood within the framework of current solvation models. On the contrary, the regime of large cavities is characterized by a significant softening of the cavity interface resulting in a decay of the fluctuation variances with the cavity size much faster than anticipated by both the continuum electrostatics and microscopic theories. For instance, the variance of the electrostatic potential at the cavity center decays with the cavity radius R0 approximately as 1R{0};{4-6} instead of the 1R{0} scaling expected from the standard electrostatics. Our results suggest that cores of nonpolar molecular assemblies in polar liquids lose solvation strength much faster than is traditionally anticipated.

Ab initio quantum mechanical/molecular mechanical simulation of electron transfer process: fractional electron approach

Electron transfer (ET) reactions are one of the most important processes in chemistry and biology. Because of the quantum nature of the processes and the complicated roles of the solvent, theoretical study of ET processes is challenging. To simulate ET processes at the electronic level, we have developed an efficient density functional theory (DFT) quantum mechanical (QM)/molecular mechanical (MM) approach that uses the fractional number of electrons as the order parameter to calculate the redox free energy of ET reactions in solution. We applied this method to study the ET reactions of the aqueous metal complexes Fe(H(2)O)(6)(2+/3+) and Ru(H(2)O)(6)(2+/3+). The calculated oxidation potentials, 5.82 eV for Fe(II/III) and 5.14 eV for Ru(II/III), agree well with the experimental data, 5.50 and 4.96 eV, for iron and ruthenium, respectively. Furthermore, we have constructed the diabatic free energy surfaces from histogram analysis based on the molecular dynamics trajectories. The resulting reorganization energy and the diabatic activation energy also show good agreement with experimental data. Our calculations show that using the fractional number of electrons (FNE) as the order parameter in the thermodynamic integration process leads to efficient sampling and validate the ab initio QM/MM approach in the calculation of redox free energies.