Theoretical Analysis of Proton Relays in Electrochemical Proton-Coupled Electron Transfer

General Kinetic Model for pH Dependence of Proton-Coupled Electron Transfer: Application to an Electrochemical Water Oxidation System

The pH dependence of proton-coupled electron transfer (PCET) reactions, which are critical to many chemical and biological processes, is a powerful probe for elucidating their fundamental mechanisms. Herein, a general, multichannel kinetic model is introduced to describe the pH dependence of both homogeneous and electrochemical PCET reactions. According to this model, a weak pH dependence can arise from the competition among multiple sequential and concerted PCET channels involving different forms of the redox species, such as protonated and deprotonated forms, as well as different proton donors and acceptors. The contribution of each channel is influenced by the relative populations of the reactant species, which often depend strongly on pH, leading to complex pH dependence of PCET apparent rate constants. This model is used to explain the origins of the experimentally observed weak pH dependence of the electrochemical PCET apparent rate constant for a ruthenium-based water oxidation catalyst attached to a tin-doped In2O3 (ITO) surface. The weak pH dependence is found to arise from the intrinsic differences in the rate constants of participating channels and the dependence of their relative contributions on pH. This model predicts that the apparent maximum rate constant will become pH-independent at higher pH, which is confirmed by experimental measurements. Our analysis also suggests that the dominant channels are electron transfer at lower pH and sequential PCET via electron transfer followed by fast proton transfer at higher pH. This work highlights the importance of considering multiple competing channels simultaneously for PCET processes.

Exploring Proton-Coupled Electron Transfer at Multiple Scales

The coupling of electron and proton transfer is critical for chemical and biological processes spanning a wide range of length and time scales and often occurring in complex environments. Thus, diverse modeling strategies, including analytical theories, quantum chemistry, molecular dynamics, and kinetic modeling, are essential for a comprehensive understanding of such proton-coupled electron transfer reactions. Each of these computational methods provides one piece of the puzzle, and all these pieces must be viewed together to produce the full picture.

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.

Structural and Thermodynamic Effects on the Kinetics of C-H Oxidation by Multisite Proton-Coupled Electron Transfer in Fluorenyl Benzoates

Our recent experimental and theoretical investigations have shown that fluorene C-H bonds can be activated through a mechanism in which the proton and electron are transferred from the C-H bond to a separate base and oxidant in a concerted, elementary step. This multisite proton-coupled electron transfer (MS-PCET) mechanism for C-H bond activation was shown to be facilitated by shorter proton donor-acceptor distances. With the goal of intentionally modulating this donor-acceptor distance, we have now studied C-H MS-PCET in the 3-methyl-substituted fluorenyl benzoate (2-Flr-3-Me-BzO^-). This derivative was readily oxidized by ferrocenium oxidants by initial C-H MS-PCET, with rate constants that were 6-21 times larger than those for 2-Flr-BzO^- with the same oxidants. Structural comparisons by X-ray crystallography and by computations showed that addition of the 3-methyl group caused the expected steric compression; however, the relevant C···O^- proton donor-acceptor distance was longer, due to a twist of the carboxylate group. The structural changes induced by the 3-Me group increased the basicity of the carboxylate, weakened the C-H bond, and reduced the reorganization energy for C-H bond cleavage. Thus, the rate enhancement for 2-Flr-3-Me-BzO^- was due to effects on the thermochemistry and kinetic barrier, rather than from compression of the C···O^- proton donor-acceptor distance. These results highlight both the challenges of controlling molecules on the 0.1 Å length scale and the variety of parameters that affect PCET rate constants.

Bimodal Evans-Polanyi Relationships in Hydrogen Atom Transfer from C(sp3)-H Bonds to the Cumyloxyl Radical. A Combined Time-Resolved Kinetic and Computational Study

The applicability of the Evans-Polanyi (EP) relationship to HAT reactions from C(sp3)-H bonds to the cumyloxyl radical (CumO•) has been investigated. A consistent set of rate constants, kH, for HAT from the C-H bonds of 56 substrates to CumO•, spanning a range of more than 4 orders of magnitude, has been measured under identical experimental conditions. A corresponding set of consistent gas-phase C-H bond dissociation enthalpies (BDEs) spanning 27 kcal mol-1 has been calculated using the (RO)CBS-QB3 method. The log kH' vs C-H BDE plot shows two distinct EP relationships, one for substrates bearing benzylic and allylic C-H bonds (unsaturated group) and the other one, with a steeper slope, for saturated hydrocarbons, alcohols, ethers, diols, amines, and carbamates (saturated group), in line with the bimodal behavior observed previously in theoretical studies of reactions promoted by other HAT reagents. The parallel use of BDFEs instead of BDEs allows the transformation of this correlation into a linear free energy relationship, analyzed within the framework of the Marcus theory. The ΔG⧧HAT vs ΔG°HAT plot shows again distinct behaviors for the two groups. A good fit to the Marcus equation is observed only for the saturated group, with λ = 58 kcal mol-1, indicating that with the unsaturated group λ must increase with increasing driving force. Taken together these results provide a qualitative connection between Bernasconi's principle of nonperfect synchronization and Marcus theory and suggest that the observed bimodal behavior is a general feature in the reactions of oxygen-based HAT reagents with C(sp3)-H donors.

Advances and challenges for experiment and theory for multi-electron multi-proton transfer at electrified solid–liquid interfaces

Understanding microscopic mechanism of multi-electron multi-proton transfer reactions at complexed systems is important for advancing electrochemistry-oriented science in the 21st century.

Strategies for Enhancing the Rate Constant of C-H Bond Cleavage by Concerted Proton-Coupled Electron Transfer

Recently selective C-H bond cleavage under mild conditions with weak oxidants was reported for fluorenyl-benzoates. This mechanism is based on multi-site concerted proton-coupled electron transfer (PCET) involving intermolecular electron transfer to an outer-sphere oxidant coupled to intramolecular proton transfer to a well-positioned proton acceptor. The electron transfer driving force depends predominantly on the oxidant, and the proton transfer driving force depends mainly on the basicity of the carboxylate, which is influenced by the substituent on the benzoate fragment. Experiments showed that the rate constants are much more sensitive to the carboxylate basicity than to the redox potential of the oxidant. Herein a vibronically nonadiabatic PCET theory is used to explain how changing the driving force for the electron and proton transfer components of the reaction through varying the oxidant and the substituent, respectively, impacts the PCET rate constant. In addition to increasing the driving force for proton transfer, enhancing the basicity of the carboxylate also decreases the equilibrium proton donor-acceptor distance, thereby facilitating the sampling of shorter proton donor-acceptor distances. This additional effect arising from the strong dependence of proton transfer on the proton donor-acceptor distance provides an explanation for the greater sensitivity of the rate constant to the carboxylate basicity than to the redox potential of the oxidant. These fundamental insights have broad implications for developing new strategies to activate C-H bonds, specifically by designing systems with shorter equilibrium proton donor-acceptor distances.

Theoretical analysis of the inverted region in photoinduced proton-coupled electron transfer

Photoinduced proton-coupled electron transfer (PCET) plays a key role in a wide range of energy conversion processes, and understanding how to design systems to control the PCET rate constant is a significant challenge. Herein a theoretical formulation of PCET is utilized to identify the conditions under which photoinduced PCET may exhibit inverted region behavior. In the inverted region, the rate constant decreases as the driving force increases even though the reaction becomes more thermodynamically favorable. Photoinduced PCET will exhibit inverted region behavior when the following criteria are satisfied: (1) the overlap integrals corresponding to the ground reactant and the excited product proton vibrational wavefunctions become negligible for a low enough product vibronic state and (2) the reaction free energies associated with the lower excited product proton vibrational wavefunctions contributing significantly to the rate constant are negative with magnitudes greater than the reorganization energy. These criteria are typically not satisfied by harmonic or Morse potentials but are satisfied by more realistic asymmetric double well potentials because the proton vibrational states above the barrier correspond to more delocalized proton vibrational wavefunctions with nodal structures leading to destructive interference effects. Thus, this theoretical analysis predicts that inverted region behavior could be observed for systems with asymmetric double well potentials characteristic of hydrogen-bonded systems and that the hydrogen/deuterium kinetic isotope effect will approach unity and could even become inverse in this region due to the oscillatory nature of the highly excited vibrational wavefunctions. These insights may help guide the design of more effective energy conversion devices.

Concerted proton-electron transfer reactions in the Marcus inverted region

Electron transfer reactions slow down when they become very thermodynamically favorable, a counterintuitive interplay of kinetics and thermodynamics termed the inverted region in Marcus theory. Here we report inverted region behavior for proton-coupled electron transfer (PCET). Photochemical studies of anthracene-phenol-pyridine triads give rate constants for PCET charge recombination that are slower for the more thermodynamically favorable reactions. Photoexcitation forms an anthracene excited state that undergoes PCET to create a charge-separated state. The rate constants for return charge recombination show an inverted dependence on the driving force upon changing pyridine substituents and the solvent. Calculations using vibronically nonadiabatic PCET theory yield rate constants for simultaneous tunneling of the electron and proton that account for the results.

Proton Discharge on a Gold Electrode from Triethylammonium in Acetonitrile: Theoretical Modeling of Potential-Dependent Kinetic Isotope Effects

The discharge of protons on electrode surfaces, known as the Volmer reaction, is a ubiquitous reaction in heterogeneous electrocatalysis and plays an important role in renewable energy technologies. Recent experiments with triethylammonium (TEAH⁺) donating the proton to a gold electrode in acetonitrile demonstrate significantly different Tafel slopes for TEAH⁺ and its deuterated counterpart, TEAD⁺. As a result, the kinetic isotope effect (KIE) for the hydrogen evolution reaction changes considerably as a function of applied potential. Herein a vibronically nonadiabatic approach for proton-coupled electron transfer (PCET) at an electrode interface is extended to heterogeneous electrochemical processes and is applied to this system. This approach accounts for the key effects of the electrical double layer and spans the electronically adiabatic and nonadiabatic regimes, as found to be necessary for this reaction. The experimental Tafel plots for TEAH⁺ and TEAD⁺ are reproduced using physically reasonable parameters within this model. The potential-dependent KIE or, equivalently, isotope-dependent Tafel slope is found to be a consequence of contributions from excited electron-proton vibronic states that depend on both isotope and applied potential. Specifically, the contributions from excited reactant vibronic states are greater for TEAD⁺ than for TEAH⁺. Thus, the two reactions proceed by the same fundamental mechanism yet exhibit significantly different Tafel slopes. This theoretical approach may be applicable to a wide range of other heterogeneous electrochemical PCET reactions.

Promoting proton coupled electron transfer in redox catalysts through molecular design

Mini-review on using the secondary coordination sphere to facilitate multi-electron, multi-proton catalysis.

Solvent-dependent transition from concerted electron-proton to proton transfer in photoinduced reactions between phenols and polypyridine Ru complexes with proton-accepting sites

The bimolecular rate coefficients (kobsq) for quenching the metal-to-ligand charge transfer excited states of two Ru polypyridine complexes containing H-bond accepting sites by six p-substituted phenols exhibit abrupt deviations from the expected linear correlations of log kobsq with phenol's Hammett σp constant. This pattern is attributed to a transition of the quenching mechanism from a concerted electron-proton transfer (EPT) to a proton transfer (PT); the latter becomes predominant for the most acidic phenols in acetonitrile, but not in dichloromethane. This assertion is supported by a detailed thermochemical analysis, which also excludes the quenching pathways involving electron transfer from phenols with or without deprotonation of phenols to the solvent, either concerted or sequential. The transition from EPT to PT upon the σp increase is consistent/supported by the magnitudes of the measured and computed PhOH/OD kinetic isotope effects and by the observed reduction of the EPT product yields upon replacing the low σp methoxyphenol by the high σp nitrophenol. In addition to modulating the relative contribution of the EPT and PT quenching pathways, the solvent strongly affects the bimolecular rate coefficients for the EPT quenching proper. Unlike with H-atom transfer reactions, this kinetic solvent effect could not be quantitatively accounted for by the phenol-solvent H-bonding alone, which suggests a solvent effect on the H-bonding constants in the phenol-Ru complex precursor exciplexes and/or on the unimolecular EPT rate coefficients within these exciplexes.

Theoretical Study of C-H Bond Cleavage via Concerted Proton-Coupled Electron Transfer in Fluorenyl-Benzoates

Developing new strategies to activate and cleave C-H bonds is important for a broad range of applications. Recently a new approach for C-H bond activation using multi-site concerted proton-coupled electron transfer (PCET) involving intermolecular electron transfer to an oxidant coupled to intramolecular proton transfer was reported. For a series of oxidants reacting with 2-(9 H-fluoren-9-yl)benzoate, experimental studies revealed an atypical Brønsted α, defined as the slope of the logarithm of the PCET rate constant versus the logarithm of the equilibrium constant or the scaled driving force. Herein this reaction is modeled with a vibronically nonadiabatic PCET theory. Hydrogen tunneling, thermal sampling of the proton donor-acceptor mode, solute and solvent reorganization, and contributions from excited vibronic states are found to play important roles. The calculations qualitatively reproduce the experimental observation of a Brønsted α significantly less than 0.5 and explain this shallow slope in terms of exoergic processes between pairs of electron-proton vibronic states. These fundamental mechanistic insights may guide the design of more effective strategies for C-H bond activation and cleavage.

Theoretical study of proton coupled electron transfer reaction in the light state of the AppA BLUF photoreceptor

The BLUF (blue light sensor using flavin adenine dinucleotide) domain is widely studied as a prototype for proton coupled electron transfer (PCET) reactions in biological systems. In this work, the photo-induced concerted PCET reaction from the light state of the AppA BLUF domain is investigated. To model the simultaneous transfer of two protons in the reaction, two-dimensional potential energy surfaces for the double proton transfer are first calculated for the locally excited and charge transfer states, which are then used to obtain the vibrational wave function overlaps and the vibrational energy levels. Contributions to the PCET rate constant from each pair of vibronic states are then analyzed using the theory based on the Fermi's golden rule. We show that, the recently proposed light state structure of the BLUF domain with a tautomerized Gln63 residue is consistent with the concerted transfer of one electron and two protons. It is also found that, thermal fluctuations of the protein structure, especially the proton donor-acceptor distances, play an important role in determining the PCET reaction rate. © 2018 Wiley Periodicals, Inc.

Diabatic model for electrochemical hydrogen evolution based on constrained DFT configuration interaction

The accuracy of density functional theory (DFT) based kinetic models for electrocatalysis is diminished by spurious electron delocalization effects, which manifest as uncertainties in the predicted values of reaction and activation energies. In this work, we present a constrained DFT (CDFT) approach to alleviate overdelocalization effects in the Volmer-Heyrovsky mechanism of the hydrogen evolution reaction (HER). This method is applied a posteriori to configurations sampled along a reaction path to correct their relative stabilities. Concretely, the first step of this approach involves describing the reaction in terms of a set of diabatic states that are constructed by imposing suitable density constraints on the system. Refined reaction energy profiles are then recovered by performing a configuration interaction (CDFT-CI) calculation within the basis spanned by the diabatic states. After a careful validation of the proposed method, we examined HER catalysis on open-ended carbon nanotubes and discovered that CDFT-CI increased activation energies and decreased reaction energies relative to DFT predictions. We believe that a similar approach could also be adopted to treat overdelocalization effects in other electrocatalytic proton-coupled electron transfer reactions, e.g., in the oxygen reduction reaction.

Proton-coupled electron transfer reactions: analytical rate constants and case study of kinetic isotope effects in lipoxygenase

A general theory has been developed for proton-coupled electron transfer (PCET), which is vital to a wide range of chemical and biological processes. This theory describes PCET reactions in terms of nonadiabatic transitions between reactant and product electron-proton vibronic states and includes the effects of thermal fluctuations of the solvent or protein environment, as well as the proton donor-acceptor motion. Within the framework of this general PCET theory, a series of analytical rate constant expressions has been derived for PCET reactions in well-defined regimes. Herein, the application of this theory to PCET in the enzyme soybean lipoxygenase illustrates the regimes of validity for the various rate constant expressions and elucidates the fundamental physical principles dictating PCET reactions. Such theoretical studies provide significant physical insights that guide the interpretation of experimental data and lead to experimentally testable predictions. A combination of theoretical treatments with atomic-level simulations is essential to understanding PCET.

Energy-Efficient Hydrogen Evolution by Fe-S Electrocatalysts: Mechanistic Investigations

The intrinsic catalytic property of a Fe-S complex toward H₂ evolution was investigated in a wide range of acids. The title complex exhibited catalytic events at -1.16 and -1.57 V (vs Fc⁺/Fc) in the presence of trifluoromethanesulfonic acid (HOTf) and trifluoroacetic acid (TFA), respectively. The processes corresponded to the single reduction of the Fe-hydride-S-proton and Fe-hydride species, respectively. When anilinium acid was used, the catalysis occurred at -1.16 V, identical with the working potential of the HOTf catalysis, although the employment of anilinium acid was only capable of achieving the Fe-hydride state on the basis of the spectral and calculated results. The thermodynamics and kinetics of individual steps of the catalysis were analyzed by density functional theory (DFT) calculations and electroanalytical simulations. The stepwise CCE or CE (C, chemical; E, electrochemical) mechanism was operative from the HOTf or TFA source, respectively. In contrast, the involvement of anilinium acid most likely initiated a proton-coupled electron transfer (PCET) pathway that avoided the disfavored intermediate after the initial protonation. Via the PCET pathway, the heterogeneous electron transfer rate was increased and the overpotential was decreased by 0.4 V in comparison with the stepwise pathways. The results showed that the PCET-involved catalysis exhibited substantial kinetic and thermodynamic advantages in comparison to the stepwise pathway; thus, an efficient catalytic system for proton reduction was established.

NRVS Studies of the Peroxide Shunt Intermediate in a Rieske Dioxygenase and Its Relation to the Native FeII O2 Reaction

The Rieske dioxygenases are a major subclass of mononuclear nonheme iron enzymes that play an important role in bioremediation. Recently, a high-spin FeIII-(hydro)peroxy intermediate (BZDOp) has been trapped in the peroxide shunt reaction of benzoate 1,2-dioxygenase. Defining the structure of this intermediate is essential to understanding the reactivity of these enzymes. Nuclear resonance vibrational spectroscopy (NRVS) is a recently developed synchrotron technique that is ideal for obtaining vibrational, and thus structural, information on Fe sites, as it gives complete information on all vibrational normal modes containing Fe displacement. In this study, we present NRVS data on BZDOp and assign its structure using these data coupled to experimentally calibrated density functional theory calculations. From this NRVS structure, we define the mechanism for the peroxide shunt reaction. The relevance of the peroxide shunt to the native FeII/O2 reaction is evaluated. For the native FeII/O2 reaction, an FeIII-superoxo intermediate is found to react directly with substrate. This process, while uphill thermodynamically, is found to be driven by the highly favorable thermodynamics of proton-coupled electron transfer with an electron provided by the Rieske [2Fe-2S] center at a later step in the reaction. These results offer important insight into the relative reactivities of FeIII-superoxo and FeIII-hydroperoxo species in nonheme Fe biochemistry.

Concerted One-Electron Two-Proton Transfer Processes in Models Inspired by the Tyr-His Couple of Photosystem II

Nature employs a Tyr_Z-His pair as a redox relay that couples proton transfer to the redox process between P680 and the water oxidizing catalyst in photosystem II. Artificial redox relays composed of different benzimidazole-phenol dyads (benzimidazole models His and phenol models Tyr) with substituents designed to simulate the hydrogen bond network surrounding the Tyr_Z-His pair have been prepared. When the benzimidazole substituents are strong proton acceptors such as primary or tertiary amines, theory predicts that a concerted two proton transfer process associated with the electrochemical oxidation of the phenol will take place. Also, theory predicts a decrease in the redox potential of the phenol by ∼300 mV and a small kinetic isotope effect (KIE). Indeed, electrochemical, spectroelectrochemical, and KIE experimental data are consistent with these predictions. Notably, these results were obtained by using theory to guide the rational design of artificial systems and have implications for managing proton activity to optimize efficiency at energy conversion sites involving water oxidation and reduction.

Nonadiabatic rate constants for proton transfer and proton-coupled electron transfer reactions in solution: Effects of quadratic term in the vibronic coupling expansion

Rate constant expressions for vibronically nonadiabatic proton transfer and proton-coupled electron transfer reactions are presented and analyzed. The regimes covered include electronically adiabatic and nonadiabatic reactions, as well as high-frequency and low-frequency proton donor-acceptor vibrational modes. These rate constants differ from previous rate constants derived with the cumulant expansion approach in that the logarithmic expansion of the vibronic coupling in terms of the proton donor-acceptor distance includes a quadratic as well as a linear term. The analysis illustrates that inclusion of this quadratic term in the framework of the cumulant expansion framework may significantly impact the rate constants at high temperatures for proton transfer interfaces with soft proton donor-acceptor modes that are associated with small force constants and weak hydrogen bonds. The effects of the quadratic term may also become significant in these regimes when using the vibronic coupling expansion in conjunction with a thermal averaging procedure for calculating the rate constant. In this case, however, the expansion of the coupling can be avoided entirely by calculating the couplings explicitly for the range of proton donor-acceptor distances sampled. The effects of the quadratic term for weak hydrogen-bonding systems are less significant for more physically realistic models that prevent the sampling of unphysical short proton donor-acceptor distances. Additionally, the rigorous relation between the cumulant expansion and thermal averaging approaches is clarified. In particular, the cumulant expansion rate constant includes effects from dynamical interference between the proton donor-acceptor and solvent motions and becomes equivalent to the thermally averaged rate constant when these dynamical effects are neglected. This analysis identifies the regimes in which each rate constant expression is valid and thus will be important for future applications to proton transfer and proton-coupled electron transfer in chemical and biological processes.

Donor-Dependent Kinetics of Interfacial Proton-Coupled Electron Transfer

The effect of the proton donor on the kinetics of interfacial concerted proton-electron transfer (CPET) to polycrystalline Au was probed indirectly by studying the rate of hydrogen evolution from trialkylammonium donors with different steric profiles, but the same pKa. Detailed kinetic studies point to a mechanism for HER catalysis that involves rate-limiting CPET from the proton donor to the electrode surface, allowing this catalytic reaction to serve as a proxy for the rate of interfacial CPET. In acetonitrile electrolyte, triethylammonium (TEAH(+)) displays up to 20-fold faster CPET kinetics than diisopropylethylammonium (DIPEAH(+)) at all measured potentials. In aqueous electrolyte, this steric constraint is largely lifted, suggesting a key role for water in mediating interfacial CPET. In acetonitrile, TEAH(+) also displays a much larger transfer coefficient (β = 0.7) than DIPEAH(+) (β = 0.4), and TEAH(+) displays a potential-dependent H/D kinetic isotope effect that is not observed for DIPEAH(+). These results demonstrate that proton donor structure strongly impacts the free energy landscape for CPET to extended solid surfaces and highlight the crucial role of the proton donor in the kinetics of electrocatalytic energy conversion reactions.

Mechanism of the reduction of the native intermediate in the multicopper oxidases: insights into rapid intramolecular electron transfer in turnover

The multicopper oxidases (MCOs) are the family of enzymes that catalyze the 4-electron reduction of O2 to H2O coupled to the four 1-electron oxidations of substrate. In the catalytic cycle electrons are transferred intramolecularly over ∼13 Å from a Type 1 (T1) Cu site that accepts electrons from substrate to a trinuclear Cu cluster (TNC) where O2 is reduced to H2O at rapid rates consistent with turnover (560 s(-1)). The oxygen reduction mechanism for the MCOs is well-characterized, whereas the rereduction is less understood. Our initial study of Rhus vernicifera Laccase (Heppner et al. J. Am. Chem. Soc. 2013, 135, 12212) experimentally established that the native intermediate (NI), the species formed upon O-O bond cleavage, is reduced with an IET rate >700 s(-1) and is the catalytically relevant fully oxidized form of the enzyme, rather than the resting state. In this report, we present kinetic and spectroscopic results coupled to DFT calculations that evaluate the mechanism of the 3 e(-)/3 H(+) reduction of NI, where all three catalytically relevant intramolecular electron transfer (IET) steps are rapid and involve three different structural changes. These three rapid IET processes reflect the sophisticated mechanistic control of the TNC to enable rapid turnover. All three IET processes are fast due to the associated protonation of the bridging oxo and hydroxo ligands, generated by O-O cleavage, to form water products that are extruded from the TNC upon full reduction, thereby defining a unifying mechanism for oxygen reduction and rapid IET by the TNC in the catalytic cycle of the MCOs.

Proton-coupled electron transfer in molecular electrocatalysis: theoretical methods and design principles

Molecular electrocatalysts play an essential role in a wide range of energy conversion processes. The objective of electrocatalyst design is to maximize the turnover frequency and minimize the overpotential for the overall catalytic cycle. Typically, the catalytic cycle is dominated by key proton-coupled electron transfer (PCET) processes comprised of sequential or concerted electron and proton transfer steps. Theoretical methods have been developed to investigate the mechanisms, thermodynamics, and kinetics of PCET processes in electrocatalytic cycles. Electronic structure methods can be used to calculate the reduction potentials and pKa's and to generate thermodynamic schemes, free energy reaction pathways, and Pourbaix diagrams, which indicate the most stable species under certain conditions. These types of calculations have assisted in identifying the thermodynamically favorable mechanisms under specified experimental conditions, such as acid strength and overpotential. Such calculations have also revealed linear correlations among the thermodynamic properties, which can be used to predict the impact of modifying the ligands, substituents, or metal centers. The thermodynamic properties can be tuned with electron-withdrawing or electron-donating substituents. Ligand modification can exploit the role of noninnocent ligands. For example, ligand protonation typically decreases the overpotential. Calculations of rate constants for electron and proton transfer, as well as concerted PCET, have assisted in identifying the kinetically favorable mechanisms under specified conditions. The concerted PCET mechanism is thought to lower the overpotential required for catalysis by avoiding high-energy intermediates. Rate constant calculations have revealed that the concerted mechanism involving intramolecular proton transfer will be favored by designing more flexible ligands that facilitate the proton donor-acceptor motion while also maintaining a sufficiently short equilibrium proton donor-acceptor distance. Overall, theoretical methods have assisted in the interpretation of experimental data and the design of more effective molecular electrocatalysts.

Intermolecular proton shuttling in excited state proton transfer reactions: insights from theory

The mechanism of intermolecular proton shuttling involved in a prototypical excited state proton transfer reaction is disclosed using DFT and TD-DFT.

Electrostatic effects on proton coupled electron transfer in oxomanganese complexes inspired by the oxygen-evolving complex of photosystem II

The influence of electrostatic interactions on the free energy of proton coupled electron transfer in biomimetic oxomanganese complexes inspired by the oxygen-evolving complex (OEC) of photosystem II (PSII) are investigated. The reported study introduces an enhanced multiconformer continuum electrostatics (MCCE) model, parametrized at the density functional theory (DFT) level with a classical valence model for the oxomanganese core. The calculated pKa's and oxidation midpoint potentials (E(m)'s) match experimental values for eight complexes, indicating that purely electrostatic contributions account for most of the observed couplings between deprotonation and oxidation state transitions. We focus on pKa's of terminal water ligands in [Mn(II/III)(H2O)6](2+/3+) (1), [Mn(III)(P)(H2O)2](3-) (2, P = 5,10,15,20-tetrakis(2,6-dichloro-3-sulfonatophenyl)porphyrinato), [Mn2(IV,IV)(μ-O)2(terpy)2(H2O)2](4+) (3, terpy = 2,2':6',2″-terpyridine), and [Mn3(IV,IV,IV)(μ-O)4(phen)4(H2O)2](4+) (4, phen = 1,10-phenanthroline) and the pKa's of μ-oxo bridges and Mn E(m)'s in [Mn2(μ-O)2(bpy)4] (5, bpy = 2,2'-bipyridyl), [Mn2(μ-O)2(salpn)2] (6, salpn = N,N'-bis(salicylidene)-1,3-propanediamine), [Mn2(μ-O)2(3,5-di(Cl)-salpn)2] (7), and [Mn2(μ-O)2(3,5-di(NO2)-salpn)2] (8). The analysis of complexes 6-8 highlights the strong coupling between electron and proton transfers, with any Mn oxidation lowering the pKa of an oxo bridge by 10.5 ± 0.9 pH units. The model also accounts for changes in the E(m)'s by ligand substituents, such as found in complexes 6-8, due to the electron withdrawing Cl (7) and NO2 (8). The reported study provides the foundation for analysis of electrostatic effects in other oxomanganese complexes and metalloenzymes, where proton coupled electron transfer plays a fundamental role in redox-leveling mechanisms.

pH-dependent reduction potentials and proton-coupled electron transfer mechanisms in hydrogen-producing nickel molecular electrocatalysts

The nickel-based P2(Ph)N2(Bn) electrocatalysts comprised of a nickel atom and two 1,5-dibenzyl-3,7-diphenyl-1,5-diaza-3,7-diphosphacyclooctane ligands catalyze H2 production in acetonitrile. Recent electrochemical experiments revealed a linear dependence of the Ni(II/I) reduction potential on pH with a slope of 57 mV/pH unit, implicating a proton-coupled electron transfer (PCET) process with the same number of electrons and protons transferred. The combined theoretical and experimental studies herein provide an explanation for this pH dependence in the context of the overall proposed catalytic mechanism. In the proposed mechanisms, the catalytic cycle begins with a series of intermolecular proton transfers from an acid to the pendant amine ligand and electrochemical electron transfers to the nickel center to produce the doubly protonated Ni(0) species, a precursor to H2 evolution. The calculated Ni(II/I) reduction potentials of the doubly protonated species are in excellent agreement with the experimentally observed reduction potential in the presence of strong acid, suggesting that the catalytically active species leading to the peak observed in these cyclic voltammetry (CV) experiments is doubly protonated. The Ni(I/0) reduction potential was found to be slightly more positive than the Ni(II/I) reduction potential, indicating that the Ni(I/0) reduction occurs spontaneously after the Ni(II/I) reduction, as implied by the experimental observation of a single CV peak. These results suggest that the PCET process observed in the CV experiments is a two-electron/two-proton process corresponding to an initial double protonation followed by two reductions. On the basis of the experimental and theoretical data, the complete thermodynamic scheme and the Pourbaix diagram were generated for this catalyst. The Pourbaix diagram, which identifies the most thermodynamically stable species at each reduction potential and pH value, illustrates that this catalyst undergoes different types of PCET processes for various pH ranges. These thermodynamic insights will aid in the design of more effective molecular catalysts for H2 production.

Theoretical Design of Molecular Electrocatalysts with Flexible Pendant Amines for Hydrogen Production and Oxidation

The design of hydrogen oxidation and production electrocatalysts is important for the development of alternative renewable energy sources. The overall objective is to maximize the turnover frequency and minimize the overpotential. We use computational methods to examine a variety of nickel-based molecular electrocatalysts with pendant amines. Our studies focus on the proton-coupled electron transfer (PCET) process involving electron transfer between the complex and the electrode and intramolecular proton transfer between the nickel center and the nitrogen of the pendant amine. The concerted PCET mechanism, which tends to require a lower overpotential, is favored by a smaller equilibrium Ni-N distance and a more flexible pendant amine ligand, thereby decreasing the energetic penalty for the nitrogen to approach the nickel center for proton transfer. Our calculations provide predictions about designing catalysts that incorporate these properties. These design principles will be useful for developing the next generation of hydrogen catalysts.

A redox active proton transfer and hydrogen-bonding system of tetrathiafulvalene-dicarboxylic acid and pyridine bases

Hydrogen bonds are important in supramolecular chemistry and responsible for proton transfer processes. This manuscript reports new hydrogen-bonding systems of redox-active acid-base assemblies consisting of dimethylthio-tetrathiafulvalene dicarboxylic acid (H(2)L) coupled with pyridine (py), 2,2'-bipyridine (2,2'-bpy), and 4,4'-bipyridine (4,4'-bpy). The (1)H NMR chemical shifts of the pyridines and the redox potential shifts of the tetrathiafulvalene moiety indicate that proton-transfer and strong hydrogen-bonding interactions exist between the acid and the bases in aprotic solvents. The results of cyclic voltammetry show a two-step square reaction with a redox and proton transfer coupled mechanism. The nature of the hydrogen bonds was characterized by X-ray single crystal analysis. In contrast to the mono-carboxyl-py system, in this ortho-substituted dicarboxyl-py system one proton of the H(2)L transfers completely from the carboxyl to the py group. The conjugated intramolecular and intermolecular hydrogen bonds are responsible for the proton transfer.

Insights into proton-coupled electron transfer mechanisms of electrocatalytic H2 oxidation and production

The design of molecular electrocatalysts for H(2) oxidation and production is important for the development of alternative renewable energy sources that are abundant, inexpensive, and environmentally benign. Recently, nickel-based molecular electrocatalysts with pendant amines that act as proton relays for the nickel center were shown to effectively catalyze H(2) oxidation and production. We developed a quantum mechanical approach for studying proton-coupled electron transfer processes in these types of molecular electrocatalysts. This theoretical approach is applied to a nickel-based catalyst in which phosphorous atoms are directly bonded to the nickel center, and nitrogen atoms of the ligand rings act as proton relays. The catalytic step of interest involves electron transfer between the nickel complex and the electrode as well as intramolecular proton transfer between the nickel and nitrogen atoms. This process can occur sequentially, with either the electron or proton transferring first, or concertedly, with the electron and proton transferring simultaneously without a stable intermediate. The electrochemical rate constants are calculated as functions of overpotential for the concerted electron-proton transfer reaction and the two electron transfer reactions in the sequential mechanisms. Our calculations illustrate that the concerted electron-proton transfer standard rate constant will increase as the equilibrium distance between the nickel and nitrogen atoms decreases and as the pendant amines become more flexible to facilitate the contraction of this distance with a lower energy penalty. This approach identifies the favored mechanisms under various experimental conditions and provides insight into the impact of substituents on the nitrogen and phosphorous atoms.

Hydrogen-bond relays in concerted proton-electron transfers

Reaction mechanisms in which electron and proton transfers are coupled are central to a huge number of processes, both natural and synthetic. Moreover, most of the new approaches to address modern energy challenges involve proton-coupled electron transfer (PCET). Recent research has focused on the possibility that the two steps are concerted, that is, concerted proton-electron transfer (CPET) reactions, rather than stepwise pathways in which proton transfer precedes (PET) or follows (EPT) electron transfer. CPET pathways have the advantage of bypassing the high-energy intermediates of stepwise pathways, although this thermodynamic benefit may have a kinetic cost. Concerted processes require short distances between the group being oxidized and the proton acceptor (and vice versa for a reduction process), which usually involves the formation of a hydrogen bond. Unlike the electron in outer-sphere electron-transfer reactions, the distance a proton may travel in a CPET is therefore rather limited. The idea has recently emerged, however, that this distance may be substantially increased via a H-bond relay located between the electron-transfer-triggered proton source and the proton acceptor. Generally speaking, the relay is a group bearing a H atom able to accept a H-bond from the moiety being oxidized and, at the same time, to form a H-bond with the proton-accepting group without going through a protonated intermediate. Although these molecules do not retain all the properties of chains of water molecules engaged in Grotthuss-type transport of a proton, the OH group in these molecules does possess a fundamental property of water molecules: namely, it is both a hydrogen-bond acceptor and a hydrogen-bond donor. Despite centuries of study, the mechanisms of proton movement in water remain active experimental and theoretical research areas, but so far with no connection to CPET reactions. In this Account, we bring together recent results concerning (i) the oxidative response of molecules containing a H-bond relay and (ii) the oxidation of phenol with water (in water) as the proton acceptor. In the first case, a nondestructive electrochemical method (cyclic voltammetry) was used to investigate the oxidation of phenol molecules containing one H-bond relay and an amine proton acceptor compared with a similar amino phenol deprived of relay. In the second, the kinetics of phenol oxidation with water (in water) as proton acceptor is contrasted with that of conventional proton acceptors (such as hydrogen phosphate and pyridine) to afford evidence of the concerted nature of Grotthuss-type proton displacement with electron transfer. First indications were provided by the same electrochemical method, whereas a more complete kinetic characterization was obtained from laser flash photolysis. Older electrochemical results concerning the reduction of superoxide ion in the presence of water are also examined. The result is a timely picture of current insight into concerted mechanisms involving electron transfer coupled with proton transport over simple H-bond relays and over H-bond networks.

Probing quantum and dynamic effects in concerted proton-electron transfer reactions of phenol-base compounds

The oxidation of three phenols, which contain an intramolecular hydrogen bond to a pendent pyridine or amine group, has been shown, in a previous experimental study, to undergo concerted proton-electron transfer (CPET). In this reaction, the electron is transferred to an outer-sphere oxidant, and the proton is transferred from the oxygen to nitrogen atom. In the present study, this reaction is studied computationally using a version of Hammes-Schiffer's multistate continuum theory where CPET is formulated as a transmission frequency between neutral and cation vibrational-electronic states. The neutral and cation proton vibrational wave functions are computed from one-dimensional potential energy surfaces (PESs) for the transferring proton in a fixed heavy atom framework. The overlap integrals for these neutral/cation wave functions, considering several initial (i.e., neutral) and final (i.e., cation) vibrational states, are used to evaluate the relative rates of oxidation. The analysis is extended to heavy atom configurations with various proton donor-acceptor (i.e., O-N) distances to assess the importance of heavy atom "gating". Such changes in d(ON) dramatically affect the nature of the proton PESs and wave functions. Surprisingly, the most reactive configurations have similar donor-acceptor distances despite the large (~0.2 Å) differences in the optimized structures. These theoretical results qualitatively reproduce the experimental faster reactivity of the reaction of the pyridyl derivative 1 versus the CH(2)-pyridyl 2, but the computed factor of 5 is smaller than the experimental 10(2). The amine derivative is calculated to react similarly to 1, which does not agree with the experiments, likely due to some of the simplifying assumptions made in applying the theory. The computed kinetic isotope effects (KIEs) and their temperature dependence are in agreement with experimental results.