Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reaction rate constants”

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.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Fundamental Studies of Bifunctional Catalysts for Tandem Reactions

During the past few years of our funding from DOE Basic Energy Sciences, we have focused on the catalytic conversion of biomass-derived 5-hydroxymethyl furfural (HMF) as a platform molecule for sustainable synthesis of value-added chemicals. The combination of the diene and aldehyde functionalities in HMF enables catalytic production of acetalized HMF derivatives with diol or epoxy reactants to allow reversible synthesis of norcantharimide derivatives upon Diels-Alder reaction with maleimides. Whereas the electron-withdrawing nature of aldehyde group inhibits the reactivity of diene group in HMF to undergo Diels-Alder reactions with dienophiles, by converting the aldehyde group in HMF to an electron-donating hydroxyl group, HMF can append maleimide-based chemicals through Diels-Alder coupling. Acetalization of HMF not only produces an electron-rich diene for Diels-Alder reaction, but it also alters the reversibility of acetal formation and thereby allows for the controlled release of chemicals by hydrolysis. Therefore, we synthesized various norcantharimide derivatives from biomass-derived HMF by acetalization over an acid catalyst, Amberlyst-15, followed by Diels-Alder reaction with maleimide. The norcantharimides release the starting materials by retro Diels-Alder reaction that is triggered by acetal hydrolysis under acidic (≤pH 3) conditions. In subsequent work, we explored the synthesis and properties of functional polyurethanes and polyesters with tunable properties from biomass-derived (HMF)-Acetone-HMF (HAH) monomers. HAH can be selectively hydrogenated over Cu and Ru catalysts to produce partially-hydrogenated HAH (PHAH) and fully-hydrogenated HAH (FHAH). HAH contains functional groups, including hydroxyl, furan, enone, and ketone functionalities, that can be exploited to further tune the polymer properties. The π-electron conjugation between the enone and furan groups in HAH monomer stabilizes the furan group to prevent the uncontrolled degradation of furans during reactions, such as etherification. Therefore, HAH is a renewable monomer that can be used to synthesize polymers, having high molecular weight, symmetric functionalities, an inexpensive production price, and being derived from renewable resources without competition with food resources. We demonstrated the synthesis of functional polyurethanes and a renewable polyester from HAH-derived monomers, and 4,4’-methylenebis(phenyl isocyanate) (MDI). We showed that it is essential to be able to selective hydrogenate C=C double bonds in HAH to be able to achieve effective Diels Alder reactions of the furan ring. Therefore, we studied the hydrogenation at temperatures from 313 – 393 K of HAH over Pd, Ru, and Cu based catalysts. HAH was selectively hydrogenated to produce partially-hydrogenated monomers (PHAH) over Cu and Ru catalysts and to fully-hydrogenated HAH monomers (FHAH) over the Ru catalyst. Pd based catalysts yielded a mixture of partially and fully hydrogenated monomers. Reaction kinetics models were employed to quantify the kinetic behavior for hydrogenation over Ru, Cu, and Pd catalysts. A 5-step pathway exhibited over Pd and Ru catalysts consists of both series and parallel reaction steps, where HAH is both converted to fully hydrogenated products sequentially via series reactions of partially hydrogenated intermediates, as well as converted directly in parallel reactions to form the fully hydrogenated products. In contrast, a 3-step pathway over the Cu catalyst consists only of the consecutive reaction steps, where the final product was formed via series reactions of intermediate products. Additionally, reaction over the Cu catalyst did not hydrogenate the furan rings of the HAH molecule and yielded a different final product than those hydrogenation over Pd and Ru catalysts. Using the results from our detailed reaction kinetics studies, we developed universal conditions for the maximum production of the various hydrogenated products for both batch and plug flow systems. In subsequent work, we employed reaction schemes to optimize the yields of products from selective hydrogenations of HAH in isopropanol solvent and in the presence of liquid water. Reaction schemes consisting of 7, 9, and 11 steps were examined to describe the rates of formation of the observed products and reaction intermediates for hydrogenation of HAH over Ru and Pd catalysts, and a 3-step scheme was studied over Cu catalysts. Rate constants and activation energies were calculated using these reaction schemes, and we then apply these reaction schemes to explore the effects of water addition on the hydrogenation pathways. The effects of water addition to isopropanol (IPA) solvents on the hydrogenation of HAH were markedly different over Pd, Ru, and Cu catalysts. Over the Pd catalyst, the addition of water to IPA increased hydrogenation rates and promoted hydrogenation of furan rings. The addition of water to IPA yielded significant carbon losses over the Ru catalyst, and slowed hydrogenation steps over Cu, while significantly inhibiting hydrogenation of the ketone group. This behavior opened routes toward increased production rates of PHAH=O (a partially hydrogenated form of HAH containing a C=O bond), a product in which the diene groups of the furan rings were not hydrogenated. Addition of water also allowed increased feed concentrations of HAH that were previously not possible in pure IPA solvents.

09 BIOMASS FUELS↗
DOE OSTI2025

Direct Evidence for Buffer-Enhanced Proton-Coupled Electron Transfer Generation of a High-Valent Metal-Oxo Complex

Here, the oxidation of metal-aquo and -hydroxo complexes to generate the high-valent metal-oxo species used in oxidative catalysis is often kinetically slow due to sluggish proton transfer between ligated −H 2 O/–OH in the proton-coupled electron transfer (PCET) chemistry. In this research, a ruthenium water oxidation catalyst anchored to a conductive tin-doped indium oxide (ITO) thin film, abbreviated ITO|Ru II –OH 2 , was characterized by spectroscopic and electrochemical methods in acetate or phosphate buffers. The deprotonated intermediate, Ru II –OH, was observed spectroscopically in the PCET half-reaction ITO(e – )|Ru III –OH + H + → ITO|Ru II –OH 2 indicating an underlying stepwise ET-PT mechanism. In contrast, at elevated buffer concentrations, this intermediate was absent, and a 2–4 order of magnitude increase in the proton transfer rate constant was observed. Kinetic data for this PCET reaction measured as a function of the driving force provided the reorganization energy λ = 1.05 eV and was assigned to a concerted electron–proton transfer (EPT) mechanism. In addition, the standard heterogeneous rate constants for two PCET equilibria, Ru III –OH + H + + e – ⇌ Ru II –OH 2 and Ru IV = O + H + + e – ⇌ Ru III –OH were enhanced by these same buffers. Collectively, the data show that the added buffers can enhance the kinetics and thermodynamics for PCET reactions relevant to oxidative catalysis.

Catalysts↗
DOE OSTI2021

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 (LIBs) 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 kinetic approach with theory and molecular simulation. We employed both the potential distribution theorem (PDT) and the SMD implicit solvent model 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. Further, 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. Furthermore, this manuscript demonstrates how ion solvation effects largely govern the prediction of reduction potentials and electron transfer rate constants at the electrode-electrolyte interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Synergistic catalysis by Mn promoted ceria for molecular oxygen assisted epoxidation

Mesoporous manganese doped cerium oxide catalysts were prepared by an inverse micelle method with different manganese loadings. These materials exhibited superior activity for epoxidation of alkenes with molecular oxygen. Highest active 10MnCe material showed an activity of 80 % conversion and >99 % selectivity for cyclooctene epoxidation. The fluorite crystal structure which facilitates the high oxygen mobility, low temperature reducibility, homogeneous distribution of manganese on ceria, nanoparticle nature, higher number of oxygen vacancies, low valent states of manganese and cerium, and the ability to generate reactive oxygen species of the manganese doped cerium oxide could be correlated to the high catalytic activity. The reaction followed first order kinetics with a rate constant of 0.31 h –1 at 100 °C, and the activation energy turned out to be 10.5 kJ/mol. In conclusion, the mechanistic study revealed that the reaction proceeds via a unique mechanism which involves in situ generated superoxide and singlet oxygen species..

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

An experimental, theoretical and kinetic modeling study of the N 2 O-H 2 system: Implications for N 2 O + H

The reaction of N 2 O with His the key step in consumption of nitrous oxide in thermal processes. The major product channel is N 2 + OH, while NH + NO constitute minor products. In addition, a pathway involving HNNO, initiated by N 2 O + H (+M)$\rightleftarrows$HNNO (+M) (R3, R4), has been inferred from experiment and theory by Burke and coworkers. At longer reaction times, the reaction may reach partial equilibration, and in addition to k 3 and k 4 the importance of this channel depends on the thermodynamic properties of HNNO and its consumption reactions, mainly HNNO + H. In the present work, we re-examined the thermochemistry of HNNO and calculated rate constants and branching fractions for the HNNO + H reaction. Experiments on the N 2 O-H 2 system were conducted in a high-pressure flow reactor at 100 atm as a function of temperature (600-925 K) and stoichiometry and explained in terms of an updated chemical kinetic model. The results support the importance of the HNNO pathway, which results in inhibition of N 2 O consumption and formation of NH 3 . In addition, selected literature results on the N 2 O-H 2 system are re-examined and the implications for the other product channels of N 2 O + H, in particular NH + NO, are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Association of Cl with C 2 H 2 by unified variable-reaction-coordinate and reaction-path variational transition-state theory

Significance A key issue in chemical kinetics is advancing the theoretical framework to handle reactions beyond the domain of textbook transition-state theory by including––for example––anharmonicity, barrierless transition states, transition states in series, and the effect of conformational flexibility on equilibrium constants. Other key issues are validating affordable electronic structure methods for direct dynamics and the direct calculation of high-pressure limiting rate constants, which are often obtainable experimentally only by extrapolation. The present article addresses all these issues for a prototype radical–molecule association reaction and thereby demonstrates how to combine improved theoretical methods to provide rate constants in cases where experimental data are uncertain or missing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Universal free-energy landscape produces efficient and reversible electron bifurcation

For decades, it was unknown how electron-bifurcating systems in nature prevented energy-wasting short-circuiting reactions that have large driving forces, so synthetic electron-bifurcating molecular machines could not be designed and built. The underpinning free-energy landscapes for electron bifurcation were also enigmatic. We predict that a simple and universal free-energy landscape enables electron bifurcation, and we show that it enables high-efficiency bifurcation with limited short-circuiting (the EB scheme). The landscape relies on steep free-energy slopes in the two redox branches to insulate against short-circuiting using an electron occupancy blockade effect, without relying on nuanced changes in the microscopic rate constants for the short-circuiting reactions. The EB scheme thus unifies a body of observations on biological catalysis and energy conversion, and the scheme provides a blueprint to guide future campaigns to establish synthetic electron bifurcation machines.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

The reaction NH 2 + CH 2 O: Kinetic measurements and computational studies

The reaction of amidogen with formaldehyde is relevant to astrochemistry and couples nitrogen and carbon chemistry in flames. The first measurements of the temperature dependence of the rate constant k were made, over 409–643 K, and are summarized as k = 7.1 × 10 −12 exp(−18.2 kJ mol −1 /RT) cm 3 molecule −1 s −1 with a confidence interval of ±20%. These results address a controversy over the presence of a significant barrier in the addition/elimination pathway leading to the production of formamide between theoretical models that include the zero-point energy of all modes at the transition state and a model that excludes transitional modes. The latter would lead to a negligible barrier, which is inconsistent with the experiments. Contrary to earlier claims, CCSDT(Q)-corrected energies combined with transition state theory yield quantitative accord with the measurements.

Astrochemistry↗
DOE OSTI2024

Quasi-Classical Trajectory Calculation of Rate Constants Using an Ab Initio Trained Machine Learning Model (aML-MD) with Multifidelity Data

Machine learning (ML) provides a great opportunity for the construction of models with improved accuracy in classical molecular dynamics (MD). However, the accuracy of a ML trained model is limited by the quality and quantity of the training data. Generating large sets of accurate ab initio training data can require significant computational resources. Furthermore, inconsistent or incompatible data with different accuracies obtained using different methods may lead to biased or unreliable ML models that do not accurately represent the underlying physics. Recently, transfer learning showed its potential for avoiding these problems as well as for improving the accuracy, efficiency, and generalization of ML models using multifidelity data. In this work, ab initio trained ML-based MD (aML-MD) models are developed through transfer learning using DFT and multireference data from multiple sources with varying accuracy within the Deep Potential MD framework. Further, the accuracy of the force field is demonstrated by calculating rate constants for the H + HO 2 → H 2 + 3 O 2 reaction using quasi-classical trajectories. We show that the aML-MD model with transfer learning can accurately predict the rate constants while reducing the computational cost by more than five times compared to the use of more expensive quantum chemistry training data sets. Hence, the aML-MD model with transfer learning shows great potential in using multifidelity data to reduce the computational cost involved in generating the training set for these potentials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Monolayer Molecular Functionalization Enabled by Acid–Base Interaction for High-Performance Photochemical CO 2 Reduction

We report the development of a hybrid catalyst consisting of carbon nitride (CN x ) and cobalt phthalocyanine tetracarboxylic acid (CoPc-COOH), which converts CO 2 to CO with high reaction rate (1067 μmol/g·h) and high selectivity (over 98%), under simulated solar irradiation. The carboxylic acid substituents on the phthalocyanine ligands play a critical role as they bind to the amine groups of CN x to enable nearly ideal monolayer coverage of the molecular co-catalyst on the semiconductor surface and promote catalytic activity from the molecular complex. Specifically, the CN x /CoPc-COOH hybrid material achieves a reaction rate 16 times higher than a CN x material containing unsubstituted CoPc molecules. We further show that activation and deactivation of the CN x /CoPc-COOH composite, which are associated with the reduction and decomposition of CoPc-COOH, respectively, both proceed at a nearly constant rate regardless of the CO 2 reduction reaction rate. Here, the decoupling of charge carrier injection and CO 2 reduction catalysis has important mechanistic implications for future performance optimization and materials design of photocatalysts for CO 2 reduction.

14 SOLAR ENERGY↗
DOE OSTI2025

Ab Initio Molecular Dynamics Study of the Reduction of Acetone by the Hydrated Electron

We have investigated the reaction dynamics of the reduction of acetone by the solvated electron in water using ab initio molecular dynamics simulations at 298 and 373 K. The rate constants derived from the simulations are consistent with experimental observations that the reaction has a low activation energy. Detailed analyses carried out to shed light on the mechanism of the reaction show that solvent reorganization plays a key role as a reaction coordinate, as expected from the Marcus electron transfer theory. Furthermore, constrained density functional theory calculations indicate that the electronic coupling is large, placing the reaction in the adiabatic limit. Indeed, the activation energies and rate constants of the simulations are in accordance with the predictions of adiabatic Marcus theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

From electronic structure to model application of key reactions for gasoline/alcohol combustion: Hydrogen-atom abstraction by $CH_3O\dot{O}$ radicals

Hydrogen atom abstraction by methyl peroxy ($CH_3O\dot{O}$) radicals can play an important role in gasoline/ethanol interacting chemistry for fuels that produce high concentrations of methyl radicals. Detailed kinetic reactions for hydrogen atom abstraction by $CH_3O\dot{O}$ radicals from the components of FGF-LLNL (a gasoline surrogate) including cyclopentane, toluene, 1-hexene, n -heptane, and isooctane have been systematically studied in this work. Here, the M06-2X/6-311 ++ G(d,p) level of theory was used to obtain the optimized structure and vibrational frequency for all stationary points and the low-frequency torsional modes. The 1-D hindered rotor treatment for low-frequency torsional modes was treated at M06-2X/6-31G level of theory. The UCCSD(T)-F12a/cc-pVDZ-F12 and QCISD(T)/CBS level of theory were used to calculate single point energies for all species. High pressure limiting rate constants for all hydrogen atom abstraction channels were performed using conventional transition state theory with unsymmetric tunneling corrections. Individual rate constants are reported in the temperature range from 298.15 to 2000 K. Our computed results show that the abstraction of allylic hydrogen atoms from 1-hexene is the fastest at low temperatures. When the temperature increases, the hydrogen atom abstraction reaction channel at the primary alkyl site gradually becomes dominant. Thermodynamics properties for all stable species and high-pressure limiting rate constants for each reaction pathway obtained in this work were incorporated into the latest gasoline surrogate/ethanol model to investigate the influence of the rate constants calculated here on model predicted ignition delay times.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Revisiting trends in the exchange current for hydrogen evolution

Nørskov and collaborators proposed a simple kinetic model to explain the volcano relation for the hydrogen evolution reaction on transition metal surfaces such that j 0 = k 0 f(ΔG H ) where j 0 is the exchange current density, f(ΔG H ) is a function of the hydrogen adsorption free energy ΔG H as computed from density functional theory, and k 0 is a universal rate constant. Herein, focusing on the hydrogen evolution reaction in acidic medium, we revisit the original experimental data and find that the fidelity of this kinetic model can be significantly improved by invoking metal-dependence on k 0 such that the logarithm of k 0 linearly depends on the absolute value of ΔG H . Here, we further confirm this relationship using additional experimental data points obtained from a critical review of the available literature. Our analyses show that the new model decreases the discrepancy between calculated and experimental exchange current density values by up to four orders of magnitude. Furthermore, we show the model can be further improved using machine learning and statistical inference methods that integrate additional material properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Mesoporous Fe-doped MgO nanoparticles as a heterogeneous photo-Fenton-like catalyst for degradation of salicylic acid in wastewater

Mesoporous Fe-doped MgO nanoparticles were synthesized using a facile sol-gel method and utilized for photo-Fenton-like degradation of salicylic acid (SA). The MgO surface dissolution facilitated an increase in the pH under the reaction conditions that allowed the Fe-MgO catalyst to be active without detectable iron leaching. Under simulated solar radiation, SA was completely degraded with an initial rate constant of 0.048 min -1 at the optimal reaction conditions of 500 ppm loading of 5% Fe-MgO, 20 mM H 2 O 2 concentration, and 50 ppm SA concentration. The catalyst was stable over 5 reaction cycles. The Fe-MgO catalyst was shown to have a surface area up to 171 m 2 /g and contain hematite (Fe 2 O 3 ) nanoparticles with octahedrally coordinated iron catalytic centers, as inferred from diffuse reflectance UV–vis measurements. Post-reaction catalyst characterization showed that some Fe 2+ was present in the catalyst due to the redox cycle during the chain initiation and propagation steps of the reaction.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Isopentane Disproportionation in Lewis Acidic Chloroaluminate Ionic Liquid

Lewis acidic chloroaluminate ionic liquid diluted in dichloromethane catalyzes the disproportionation of alkanes, such as isopentane, via carbenium ions in a reaction readily initiated by carbenium ion precursors such as tert-butyl chloride (TBC). The carbenium ion-AlCl 4 – ion-pairs stabilized by the ionic liquid-dichloromethane solution are the key intermediates in two distinctive kinetic regimes, i.e., a transient regime (0–5 min) and a steady-state regime (after 5 min). The transient regime constitutes the majority of isopentane conversion and is governed by the initial carbenium ion concentration. In the steady-state regime, disproportionation occurs at a considerably lower rate, affected by the carbenium ion concentration, the concentration of the ionic liquid, and the reaction temperature. The formation of alkenes observed in the 1 H NMR spectra of the reacting substrates, along with the DFT calculations, suggests that deprotonation of carbenium ion-pairs reduces their concentration, decreasing, in turn, the reaction rate. Kinetic modeling indicates that the transient regime is significantly controlled by the hydride transfer (k HT ) and the deprotonation rate constants (k DP ), while the steady-state regime is additionally influenced by the alkene protonation rate constant (k –DP ). The overall activation energy of the reaction at the steady state, expressed as E a,steady-state regime = E a,HT – E a,DP + E a,–DP , was 54 kJ/mol. The reaction mechanism and the kinetics highlight the potential of Lewis acid-catalyzed conversions of hydrocarbons under remarkably mild conditions.

Deprotonation↗
DOE OSTI2024

Kinetic Analysis of Proton-Coupled Electron Transfer at an Electrode-Immobilized Complex

An alkyne-terminated cobalt complex, [Co(Cp)(dppe ≡H )(Cl)] + (Cp = cyclopentadienyl; dppe ≡H = 1,2-bis-(di-(4-ethynyl-phenyl)phosphino)ethane), (Co ≡H ) was immobilized onto a glassy carbon electrode using two attachment strategies: Cu(I) catalyzed azide–alkyne click chemistry and reductive electropolymerization. The modified electrodes prepared through reductive electropolymerization exhibit current densities and peak resolutions for the electrochemical reduction of the cobalt species, which are amenable to electroanalytical quantification of coupled chemical reactions. Through peak shift analysis of cyclic voltammograms recorded in the presence of 4-chloroanilinium tetrafluoroborate, we quantified the proton transfer rate constant for the stepwise proton-coupled electron transfer reaction that reduces the electrode-immobilized [Co(Cp)(dppe ≡H )(Cl)] + to [H–Co(Cp)(dppe ≡H )] + (k PT app = (9.3 ± 1.8) × 10 5 M –1 s –1 ). In conclusion, the extraction of kinetic parameters for an elementary proton-coupled electron transfer reaction of an electrode-immobilized complex represents the first experimental measurement of its type and lays crucial groundwork for kinetic analyses of hybrid catalyst–electrode architectures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Experimental and computational studies of the kinetics of the reaction of hydrogen peroxide with the amidogen radical

In this work, the pulsed-laser photolysis/laser-induced fluorescence method is used to study the kinetics of the reaction of NH 2 with H 2 O 2 to yield a second-order rate constant of (2.42 ± 0.55) × 10 –14 cm 3 molecule –1 s –1 at 412 K in 10–22 mbar of Ar bath gas. There are no prior measurements for comparison. To check this value and enable reliable extrapolation to other temperatures, we also compute thermal rate constants for this process over the temperature range 298–3000 K via multi-structural canonical variational transition-state theory with small-curvature multidimensional tunneling (MS-CVT/SCT). The CVT/SCT rate constants are derived using a dual-level direct dynamics approach utilizing single-point CCSD(T)-F12b/cc-pVQZ-F12 energies—corrected for core-valence and scalar relativistic effects—and M06-2X/MG3S geometries, gradients, and Hessians—for all stationary and non-stationary points along the reaction path. The multistructural method with torsional anharmonicity, based on a coupled torsional potential, is then employed to calculate correction factors for the rate constants, accounting for the comprehensive effects of torsional anharmonicity on the kinetics of this reaction system. The final MS-CVT/SCT rate constants are found to be in good agreement with our measurements and can be expressed in modified Arrhenius form as 2.13 × 10 –15 (T/298 K) 4.02 exp(–513 K/T) cm 3 molecule –1 s –1 over the temperature range of 298–3000 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Electron transfer in a crystalline cytochrome with four hemes

Diffusion of electrons over distances on the order of 100 μm has been observed in crystals of a small tetraheme cytochrome (STC) from Shewanella oneidensis [J. Huang et al . J. Am. Chem. Soc. 142, 10459–10467 (2020)]. Electron transfer between hemes in adjacent subunits of the crystal is slower and more strongly dependent on temperature than had been expected based on semiclassical electron-transfer theory. Here, in this work, we explore explanations for these findings by molecular-dynamics simulations of crystalline and monomeric STC. New procedures are developed for including time-dependent quantum mechanical energy differences in the gap between the energies of the reactant and product states and for evaluating fluctuations of the electronic-interaction matrix element that couples the two hemes. Rate constants for electron transfer are calculated from the time- and temperature-dependent energy gaps, coupling factors, and Franck–Condon-weighted densities of states using an expression with no freely adjustable parameters. Back reactions are considered, as are the effects of various protonation states of the carboxyl groups on the heme side chains. Interactions with water are found to dominate the fluctuations of the energy gap between the reactant and product states. The calculated rate constant for electron transfer from heme IV to heme Ib in a neighboring subunit at 300 K agrees well with the measured value. However, the calculated activation energy of the reaction in the crystal is considerably smaller than observed. We suggest two possible explanations for this discrepancy. The calculated rate constant for transfer from heme I to II within the same subunit of the crystal is about one-third that for monomeric STC in solution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗