Concerted Electron-Nuclear Motion in Proton-Coupled Electron Transfer-Driven Grotthuss-Type Proton Translocation
Not Available
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Not Available
The development of next generation energy conversion and catalytic systems requires fundamental understanding of the interplay between photo-excitation and the resulting proton and electron transfer processes that occur in solution. To address these challenges requires accurate and efficient methods to compute ground and excited states, as well as the ability to treat the dynamics of energy transfer in the presence of solvent fluctuations. Our team is developing accurate and efficient theoretical models for solution phase reactions. These developments provide an improved understanding of photo-initiated excitations, electron transfer, and proton coupled electron transfer processes. Our research team brings together expertise in electron and nuclear dynamics (classical and quantum), electronic structure, and analysis of solvation networks. We are developing and validating techniques for modeling condensed phase charge and energy transfer and to simulate the methods to monitor them using optical and vibrational spectroscopy.
Stereocontrol is of critical importance in organic synthesis. In this study, we demonstrate how heterogeneous electrochemical hydrogenation enables diastereocontrol simply by tuning electrochemical hydrogenation mechanisms without altering the adsorption conformation of a reactant on the electrode. We use 4-hydroxy-1-tetralone (4-OH-tetralone) as a model reactant, where diastereomers can be produced during the hydrogenation of the carbonyl group. In traditional thermocatalytic hydrogenation, H 2 first dissociates on the catalyst surface to form surface-adsorbed hydrogen (H*), and therefore, H* is always added to the organic reactant from the catalyst side via hydrogen atom transfer (HAT). Thus, in order to flip the diastereoselectivity, the adsorbed reactant itself must be physically flipped. In contrast, electrochemical hydrogenation can occur either via HAT, where H is added from the electrode surface, or via proton-coupled electron transfer (PCET), where H is added from the solution side of the adsorbed reactant. Thus, without changing the adsorption conformation of the reactant, opposite diastereomers can be obtained by switching the hydrogenation mechanism (HAT vs PCET). In this work, using a combination of experimental and computational studies, we demonstrate two examples of flipping diastereoselectivity by different electrochemical hydrogenation mechanisms. In the first case, we achieve opposite diastereoselectivities using metals that adopt different hydrogenation mechanisms (HAT vs PCET). In the second case, we flip the diastereoselectivity by varying the applied potential, which switches one hydrogenation mechanism to the other on the same metal electrode. In each case, our results offer an atomic-level understanding of the preferred hydrogenation mechanism that enables the corresponding diastereoselectivity.
Biology uses relatively few electron-transfer cofactors, tuning their potentials, electronic couplings, and reorganization energies to carry out the required chemistry. It is remarkable that the potential ordering of two-electron transfer active flavins can be normal (first oxidation at low potential and second oxidation at high potential) or inverted, and the gap between the potentials can be as large as one volt. Analysis based on structural bioinformatics and electrostatics indicates that the ordering of the flavin redox potential is influenced by protein electrostatics. In all 36 flavoproteins examined, the introduction of a negative charge near the flavin in silico increases the extent of potential inversion (by lowering the electrochemical potential of the second electron-transfer step); the introduction of a positive charge near the flavin favors normally ordered potentials. We also find that the addition of positive charges increases the electrochemical potential for the naturally occurring one-electron transition in flavodoxins (between deprotonated hydroquinone and neutral semiquinone) and also increases the second one-electron transition in bifurcating flavins (between anionic semiquinone and fully oxidized flavin). Finally, we find that proximity of a proton acceptor, notably conserved arginine, supports proton-coupled electron transfer because it may act as a proton acceptor, promoting potential inversion. This key arginine residue may enable two-electron transfer chemistry by promoting the proton-coupled electron transfer process over the pure electron transfer process, suggesting how a protein's flavin environment may influence one- or two-electron chemistry in flavoproteins.
We report the results of modeling CO 2 reduction (CO 2 R) to CO over Ag(110) and Cu(211) surfaces at different applied potentials using grand-canonical density functional theory (GC-DFT), a method specifically designed to accurately model electrochemical systems. In addition to demonstrating GC-DFT’s ability to accurately model electrochemical processes, we also compare it with the computational hydrogen electrode (CHE) approach. GC-DFT predicts that the geometries of these reacting systems strongly depend on the applied potential and the Helmholtz free energies vary nonlinearly with the applied potential, which contradicts a central assumption of the CHE approach. The CHE approach neglects the change in the number of electrons on the electrode surface at different applied potentials, which reduces its accuracy as the potential changes from the potential of zero charge. Our results further demonstrate that the grand free energies of the reaction intermediates not only depend on the value of the applied potential but also on the metal surface type, adsorption site, and adsorbate. GC-DFT’s ability to predict the effect of the applied potential on adsorbate geometry enables it to evaluate different possible reaction mechanisms at different applied potentials. For instance, GC-DFT predicts that the first step of CO2R likely switches from proton-coupled electron transfer to sequential electron transfer and then proton transfer at more reducing potentials, a result that cannot be determined by the CHE because it assumes that all electron transfers are coupled to proton transfers and neglects the effect of the applied potential on the adsorbate geometry.
Computational tools provide a unique opportunity to study and design optimal materials by enhancing our ability to comprehend the connections between their atomistic structure and functional properties. However, designing materials with tailored functionalities is complicated due to the necessity to integrate various computational-chemistry software (not necessarily compatible with one another), the heterogeneous nature of the generated data, and the need to explore vast chemical and parameter spaces. The latter is especially important to avoid bias in scattered data points-based models and derive statistical trends only accessible by systematic datasets. Here, we introduce a robust high-throughput multi-scale computational infrastructure coined MISPR (Materials Informatics for Structure–Property Relationships) that seamlessly integrates classical molecular dynamics (MD) simulations with density functional theory (DFT). By enabling high-performance data analytics and coupling between different methods and scales, MISPR addresses critical challenges arising from the needs of automated workflow management and data provenance recording. The major features of MISPR include automated DFT and MD simulations, error handling, derivation of molecular and ensemble properties, and creation of output databases that organize results from individual calculations to enable reproducibility and transparency. In this work, we describe fully automated DFT workflows implemented in MISPR to compute various properties such as nuclear magnetic resonance chemical shift, binding energy, bond dissociation energy, and redox potential with support for multiple methods such as electron transfer and proton-coupled electron transfer reactions. The infrastructure also enables the characterization of large-scale ensemble properties by providing MD workflows that calculate a wide range of structural and dynamical properties in liquid solutions. MISPR employs the methodologies of materials informatics to facilitate understanding and prediction of phenomenological structure–property relationships, which are crucial to designing novel optimal materials for numerous scientific applications and engineering technologies.
The [FeFe]-hydrogenases are enzymes that catalyze the reversible activation of H 2 coupled to the reduction–oxidation of electron carriers. Members of the different taxonomic groups of [FeFe]-hydrogenases display a wide range of preference, or bias, for H 2 oxidation or H 2 production reactions, despite sharing a common catalytic cofactor, or H-cluster. Identifying the properties that control reactivity remains an active area of investigation, and models have emerged that include diversity in the catalytic site coordination environments and compositions of electron transfer chains. The kinetics of proton-coupled electron transfer at the H-cluster might be expected to be a point of control of reactivity. To test this hypothesis, systematic changes were made to the conserved cysteine residue that functions in proton exchange with the H-cluster in the three model enzymes: CaI, CpII, and CrHydA1. CaI and CpII both employ electron transfer accessory clusters but differ in bias, whereas CrHydA1 lacks accessory clusters having only the H-cluster. Changing from cysteine to either serine (more basic) or aspartate (more acidic) modifies the sidechain pKa and thus the barrier for the proton exchange step. The reaction rates for H 2 oxidation or H 2 evolution were surveyed and measured for model [FeFe]-hydrogenases, and the results show that the initial proton-transfer step in [FeFe]-hydrogenase is tightly coupled to the control of reactivity; a change from cysteine to more basic serine favored H 2 oxidation in all enzymes, whereas a change to more acidic aspartate caused a shift in preference toward H 2 evolution. Overall, the changes in reactivity profiles were profound, spanning 10 5 in ratio of the H 2 oxidation-to-H 2 evolution rates. The fact that the change in reactivity follows a common trend implies that the effect of changing the proton-transfer residue pKa may also be framed as an effect on the scaling relationship between the H-cluster di(thiolmethyl)amine (DTMA) ligand pKa and E m values of the H-cluster. Experimental observations that support this relationship, and how it relates to catalytic function in [FeFe]-hydrogenases, are discussed.
Proton-coupled electron transfer (PCET) reactions on metal oxides require coupling between proton transfer at the solid-liquid interface and electron transfer involving defects at or near the band edge. Herein, hybrid functional periodic density functional theory is used to elucidate the impact of proton-coupled defects on the bond dissociation free energies (BDFEs) of O-H bonds on anatase TiO 2 surfaces. These O-H BDFEs are directly related to interfacial PCET thermochemistry. Comparison between geometrically similar O-H bonds associated with different defect types, namely conduction d-band electrons or valence p-band holes, reveals that the BDFEs differ by ~81 kcal/mol (3.50 eV), comparable to the wide TiO 2 band gap. These differences are shown to be determined primarily by differences in electron transfer driving forces, which are analyzed by using band energies and inner-sphere reorganization energies within a Marcus theory framework. Furthermore, these fundamental insights about the impact of proton-coupled defects on PCET thermochemistry at semiconductor surfaces have broad implications for electrocatalysis.
In photoredox reactions, charge recombination (CR) limits quantum yields, hindering the efficient conversion of light energy into catalytic activity. To address this, we drew inspiration from redox relays in photosystem II (PSII) and developed a new series of iridium(III) complexes featuring covalently attached benzimidazole-phenol-pyridine (BIP-Py) groups to facilitate intramolecular multiproton-coupled electron transfer (MPCET). Herein, we evaluate the effects of MPCET through an extended and well-defined hydrogen-bond network to improve photocatalytic activity and mitigate rapid charge recombination. Infrared spectroelectrochemistry reveals pyridine protonation upon phenol oxidation, while visible spectroelectrochemistry and transient absorption spectroscopy confirm the electro- and photochemical formation of chargeseparated states (CSS) involving oxidized BIP, resulting from intramolecular proton-coupled electron transfer (PCET). The application of the BIP-Py platform in a photocatalytic Nhydroxyphthalimide ester reduction reaction resulted in a ∼106-fold reduction in CR rate and a quantum yield enhancement of up to 157%. Our findings suggest that incorporating MPCET-based redox relays into photocatalyst frameworks is an effective strategy to enhance the efficiency of photocatalytic systems.
Simulating the nuclear–electronic quantum dynamics of large-scale molecular systems in the condensed phase is key for studying biologically and chemically important processes such as proton transfer and proton-coupled electron transfer reactions. Herein, the real-time nuclear–electronic orbital time-dependent density functional theory (RT-NEO-TDDFT) approach is combined with a hybrid quantum mechanical/molecular mechanical (QM/MM) strategy to enable the accurate description of coupled nuclear–electronic quantum dynamics in the presence of heterogeneous environments such as solvent or proteins. The densities of the electrons and quantum protons are propagated in real time, while the other nuclei are propagated classically on the instantaneous electron–proton vibronic surface. This approach is applied to phenol bound to lysozyme, intramolecular proton transfer in malonaldehyde, and nonequilibrium excited-state intramolecular proton transfer in o-hydroxybenzaldehyde. Furthermore, these examples illustrate that the RT-NEO-TDDFT framework, coupled with an atomistic representation of the environment, allows the simulation of condensed-phase systems that exhibit significant nuclear quantum effects.
Heterogenization of molecular electrocatalysts offers an attractive way to improve the catalytic selectivity and efficiency of CO 2 conversion to liquid fuels. Herein, we employ density functional theory to compare the mechanism of CO 2 RR by a cobalt(II) tetra(amino)phthalocyanine (Co(II)Pc(NH 2 ) 4 ) electrocatalyst with and without the presence of fullerene support. Our DFT calculations suggest that the CO 2 reduction mechanism is initiated by a metal-based electron reduction followed by subsequent CO 2 nucleophilic addition, electron transfer, proton transfer, water dissociation, and proton-coupled electron transfer steps that lead to CO and methanol formation. We show that graphitic interactions between the Co(II)Pc(NH 2 ) 4 electrocatalyst and C 60 support selectively improve the CO 2 RR to methanol at mild potentials. The undesirable hydrogen evolution reaction (HER) was also investigated for both electrocatalysts and proceeds via the protonation of the cobalt metal center over the nitrogen atom in the inner ring. The competition between the HER and the CO 2 RR was improved in favor of CO and methanol formation using the Co(II)Pc(NH 2 ) 4 @C 60 electrocatalyst. Overall, our results suggest C 60 as a promising graphitic support for molecular electrocatalysts integration for CO 2 catalysis.
A cobalt hangman porphyrin system with a phenyl spacer between the porphyrin ring and an internal carboxylic acid group as well as its non-hangman analogue were synthesized and utilized for the study of the proton-coupled electron transfer (PCET) kinetics attendant to electrocatalytic hydrogen evolution. Cyclic voltammetry (CV) together with simulations show that a short distance between the proton relay and the redox active cobalt center as well as the increased proton donating strength results in superior catalytic activity. Here, the mechanism of hydrogen generation is at the nexus of proton transfer–electron transfer (PTET) and concerted proton–electron transfer (CPET), as opposed to an ETPT mechanism that is characteristic of hangman systems with longer proton relay networks.
Here, the proton-coupled electron transfer (PCET) mechanism for the reaction M ox -OH + e − + H + → M red -OH 2 was de-termined through kinetic resolution of the independent electron transfer (ET) and proton transfer (PT) steps. The reaction of interest was triggered by visible light excitation of [Ru II (tpy)(bpy′)H 2 O] 2+ , Ru II -OH 2 , where tpy is 2,2′:6′,2″-terpyridine and bpy′ is 4,4′-diaminopropylsilatrane-2,2′-bipyridine, anchored to In 2 O 3 :Sn (ITO) thin films in aqueous solutions. Interfacial kinetics for the PCET reduction reaction were quantified by nanosecond transient absorption spectroscopy as a function of solution pH and applied potential. Data acquired from pH = 5-10 revealed a stepwise electron-transfer proton-transfer (ET-PT) mechanism, while kinetic measurements made below the pK a (Ru III -OH/OH 2 ) = 1.3 were used to study the analogous interfacial reaction where electron transfer was the only mechanistic step. Analysis of this data with a recently reported multi-channel kinetic model was used to construct a PCET zone diagram and supported the assignment of an ET-PT mechanism from pH = 5-10. Ultimately, this study represents a unique example amongst M ox -OH/M red -OH 2 reactivity where the protona-tion and oxidation state of the intermediate was kinetically and spectrally resolved to firmly establish the PCET mechanism.
Upon photoinitiated electron transfer, charge recombination limits the quantum yield of photoredox reactions for which the rates for the forward reaction and back electron transfer are competitive. Taking inspiration from a proton-coupled electron transfer (PCET) process in Photosystem II, here a benzimidazole-phenol (BIP) has been covalently attached to the 2,2'-bipyridyl ligand of [Ir(dF(CF 3 )ppy) 2 (bpy)][PF 6 ] (dF(CF 3 )ppy = 2-(2,4-difluorophenyl)-5-(trifluoromethyl)pyridine; bpy = 2,2'-bipyridyl). Excitation of the [Ir(dF(CF 3 )ppy) 2 (BIP-bpy)][PF 6 ] photocatalyst results in intramolecular PCET to form a charge-separated state with oxidized BIP. Subsequent reduction of methyl viologen dication (MV 2+ ), a substrate surrogate, by the reducing moiety of the charge separated species demonstrates that the inclusion of BIP significantly slows the charge recombination rate. The effect of ~24-fold slower charge recombination in a photocatalytic phthalimide ester reduction resulted in a greater than 2-fold increase in reaction quantum efficiency.
Nitrogenase catalyzes biological nitrogen fixation, the conversion of atmospheric N 2 into bioavailable ammonia. The three nitrogenase isozymes—Mo-nitrogenase, V-nitrogenase, and Fe-nitrogenase—utilize catalytic cofactors distinguished by their metal composition (Fe 7 M, M = Mo, V, or Fe; denoted FeM-co). Their catalytic cycles involve stepwise addition of 8[e−/H+] to FeM-co, generating intermediates designated E n , where n is the number of [e − /H + ] delivered. The electron-transfer has been extensively characterized, but the proton delivery has not. Here, we investigate [e − /H + ] delivery during early-stage conversions, primarily E 0 → E 1 (H), for each of the three nitrogenases, using as reductants γ-ray-generated thermolyzed, mobile electrons at 77 K, and radiation-generated solvent radicals during subsequent annealing to higher temperatures. Our results show E 0 → E 1 (H) conversion differs among the three MFe-proteins. The FeMo-co of MoFe-protein accepts an electron (ET) during 77 K γ-irradiation, but proton transfer (PT) to generate E 1 (H) is only enabled by conformational or thermodynamic activation upon cryoannealing to ∼200 K(ET/PT). For VFe-protein, E 1 (H) forms during annealing at-and-above 210 K by electron-transfer to FeV-co from radicals through proton-coupled electron transfer (PCET), which too is enabled by activated proton transfer. FeFe-protein differs in directly exhibiting delivery of protons at 77 K, which together with the mobile electrons react to form E 1 (H). This could well occur by PCET at 77 K, but does not preclude the possibility of sequential 77 K electron/proton transfer (ET/PT). In addition, 450 nm photolysis reveals the E 1 (H) state of FeV-co, like that of FeFe-co, contains a hydride bound to a formally oxidized cofactor. The mechanistic differences observed here provide a contribution towards understanding the sources of catalytic differences among the three nitrogenase isozymes.
An azadithiolate bridged CN — bound pentacarbonyl bis-iron complex, mimicking the active site of [Fe–Fe] H 2 ase is synthesized. The geometric and electronic structure of this complex is elucidated using a combination of EXAFS analysis, infrared and Mössbauer spectroscopy and DFT calculations. The electrochemical investigations show that complex 1 effectively reduces H + to H 2 between pH 0–3 at diffusion-controlled rates (10 11 M —1 s —1 ) i.e. 10 8 s —1 at pH 3 with an overpotential of 140 mV. Electrochemical analysis and DFT calculations suggests that a CN— ligand increases the p K a of the cluster enabling hydrogen production from its Fe(I)–Fe(0) state at pHs much higher and overpotential much lower than its precursor bis-iron hexacarbonyl model which is active in its Fe(0)–Fe(0) state. The formation of a terminal Fe–H species, evidenced by spectroelectrochemistry in organic solvent, via a rate determining proton coupled electron transfer step and protonation of the adjacent azadithiolate, lowers the kinetic barrier leading to diffusion controlled rates of H 2 evolution. The stereo-electronic factors enhance its catalytic rate by 3 order of magnitude relative to a bis-iron hexacarbonyl precursor at the same pH and potential.
Titanium dioxide (TiO 2 ) has long been employed as (photo)electrodes for reactions relevant to energy storage and renewable energy synthesis. Proton-coupled electron transfer (PCET) reactions with equimolar amounts of protons and electrons at the TiO 2 surface or within the bulk structure lie at the center of these reactions. Because a proton and an electron are thermochemically equivalent to an H-atom, these reactions are essentially H-atom transfer reactions. Thermodynamics of H-atom transfer has a complex dependence on the synthetic protocol and chemical history of the electrode, the reaction medium, and many others; together, these complications preclude the understanding of the H-atom transfer thermochemistry with atomic-level structural knowledge. Herein, we report our success in employing open-circuit potential (E OCP ) measurements to quantitatively determine the H-atom transfer thermochemistry at structurally well-defined Ti-oxo clusters within a colloidally stabilized metal–organic framework (MOFs), Ti-MIL-125. The free energy to transfer H-atom, Ti 3+ O–H bond dissociation free energy (BDFE), was measured to be 68(2) kcal mol -1 . To the best of our understanding, this is the first report on using E OCP measurements to quantify thermochemistry on any MOFs. The proton topology, the structural change upon the redox reaction, and BDFE values were further quantitatively corroborated using computational simulations. Furthermore, comparisons of the E OCP -derived BDFEs of Ti-MIL-125 to similar parameters in the literature suggest that E OCP should be the preferred method for quantitatively accurate BDFE calculations. Here, the reported success in employing E OCP for nanosized Ti- MIL-125 should lay the ground for thermochemical measurements of other colloidal systems, which are otherwise challenging. Implications of these measurements on Ti-MIL-125 as an H-atom acceptor in chemical reactions and comparisons with other MOFs/metal oxides are discussed.
Abstract Several renewable energy schemes aim to use the chemical bonds in abundant molecules like water and ammonia as energy reservoirs. Because the O-H and N-H bonds are quite strong (>100 kcal/mol), it is necessary to identify substances that dramatically weaken these bonds to facilitate proton-coupled electron transfer processes required for energy conversion. Usually this is accomplished through coordination-induced bond weakening by redox-active metals. However, coordination-induced bond weakening is difficult with earth’s most abundant metal, aluminum, because of its redox inertness under mild conditions. Here, we report a system that uses aluminum with a redox non-innocent ligand to achieve significant levels of coordination-induced bond weakening of O-H and N-H bonds. The multisite proton-coupled electron transfer manifold described here points to redox non-innocent ligands as a design element to open coordination-induced bond weakening chemistry to more elements in the periodic table.