Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Modeling Electrochemical Reaction”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Marcus–Hush–Chidsey kinetics at electrode–electrolyte interfaces

Electrochemical kinetics at electrode–electrolyte interfaces limit the performance of devices including fuel cells and batteries. While the importance of moving beyond Butler–Volmer kinetics and incorporating the effect of electronic density of states of the electrode has been recognized, a unified framework that incorporates these aspects directly into electrochemical performance models is still lacking. In this work, we explicitly account for the density functional theory-calculated density of states numerically in calculating electrochemical reaction rates for a variety of electrode–electrolyte interfaces. We first show the utility of this for two cases related to Li metal electrodeposition and stripping on a Li surface and a Cu surface (anode-free configuration). The deviation in reaction rates is minor for cases with flat densities of states such as Li, but is significant for Cu due to nondispersive d-bands creating large variation. Finally, we consider a semiconducting case of a solid-electrolyte interphase consisting of LiF and Li2CO3 and note the importance of the Fermi level at the interface pinned by the redox reaction occurring there. We identify the asymmetry in reaction rates as a function of discharge/charge naturally within this approach.

Kurchin, Rachel (ORCID:0000000221474809)↗
DOE OSTI2023

Isothermal Microcalorimetric Analysis of Li/CF x Battery Discharge

Lithium/carbon monofluoride (Li/CF x ) batteries generate substantial amounts of heat during discharge in part due to the large deviation of the loaded voltage from the measured and thermodynamically predicted open circuit voltage. Here, in this study, we further analyze this system by estimating the equilibrium voltage (V eq ) and its temperature dependence (dV eq /dT) over the entire discharge range. Based on these results, the ohmic and entropic heat contributions to the overall heat flow are calculated from experimental data. The ohmic heat flow is consistent with an electrochemical impedance spectroscopy model and electrode density measurements, indicating decreasing porosity of the electrode during discharge. Entropic heating increases during discharge as the cell reaction becomes increasingly entropically unfavorable. The total energy dissipated by the cell (electrical and thermal) remains similar over the entire discharge, evidence for a two-phase reaction with a constant, and rate-independent, ΔG rxn of about -465 kJ mol -1 . However, near the end of discharge (>90% DOD), ΔG rxn changes significantly indicating possible secondary reactions.

25 ENERGY STORAGE↗
DOE OSTI2025

Off-Equilibrium Reactivity of Boron-Enriched Metal Diboride Surfaces in Electroreduction Conditions

Boron-based materials, featuring B-dependent reactivity and diverse phases, are emerging as promising catalyst systems. However, the catalytic mechanism on many borides remains poorly understood due to complex surface reconstructions under reaction conditions. Here, we investigate the MoB 2 surface in conditions of hydrogen evolution reaction in acidic media, using grand canonical global optimization, grand canonical density functional theory, ab initio molecular dynamics, free energy surface sampling, and an analytical model for electrochemical barrier evaluation. We propose a boron-enrichment strategy to tune the surface reactivity of the hexagonal face of MoB 2 . We reveal the dynamic nature of the B-enriched surface under H coverage and kinetic trapping of the system in the metastable regime with an extensive examination of the deactivation pathways. The metastable center B site on B-enriched surfaces, featuring buckled-up configuration and a usual relaxation effect, is found to be highly active toward HER via the Volmer–Heyrovsky mechanism. In conclusion, this work demonstrates how off-equilibrium behaviors can arise from the interplay between adsorbate coverage and surface reconstruction on a seemingly simple surface, and we present a theoretical framework and computational workflows to address these behaviors, along with other realistic complexities, in kinetics simulations.

Adsorption↗
DOE OSTI2021

Intermediate Sr 2 Co 1.5 Fe 0.5 O 6-δ Tetragonal Structure between Perovskite and Brownmillerite as a Model Catalyst with Layered Oxygen Deficiency for Enhanced Electrochemical Water Oxidation

The generation of hydrogen in an environmentally benign way is highly essential to meet future energy demands. However, in the process of splitting water electrochemically, sluggish kinetics of the oxygen evolution reaction (OER) curtails its applicability, as it drags energy input. Herein, we synthesized SrCo-Fe-O oxides to optimize their OER activity by varying the Co/Fe ratio. Among them, Sr 2 Co 1.5 Fe 0.5 O 6-δ exhibited the best OER catalytic activity in the series, with an overpotential of 318 mV at 10 mA cm -2 and Tafel slope of 44.8 mV dec -1 . High-resolution neutron powder diffraction analysis identified an intermediate structure between the perovskite and brownmillerite, with alternating layers of disorderly orientated oxygen-deficient tetrahedra and fully stoichiometric octahedra. The unique stacking of tetrahedral and octahedral units facilitates desired interactions between the electrode surface and electrolyte. Theoretical calculations revealed that increased covalency of Co 3d and O 2p in Sr 2 Co 1.5 Fe 0.5 O 6-δ oxide is another primary contributor to its augmented water oxidation ability. As a model for developing catalysts with such an intermediate structure, the synergetic effect of oxygen vacancy and hybridization between Co 3d and O 2p assured the Sr 2 Co 1.5 Fe 0.5 O 6-δ oxide as a better catalyst for its enhanced OER activity.

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

Characterizing the Geometry and Quantifying the Impact of Nanoscopic Electrocatalyst/Semiconductor Interfaces under Solar Water Splitting Conditions

The materials that are receiving the most attention in photoelectrochemical water splitting are metallic nanoparticle electrocatalysts (np-EC) attached to the surface of a semiconductor (SC) light absorber. In these multicomponent systems, the interface between the semiconductor and electrocatalysts critically affects performance. However, the np-EC/SC interface remains poorly understood as it is complex on atomic scales, dynamic under reaction conditions, and inaccessible to direct experimental probes. This contribution sheds light on how the electrocatalyst/semiconductor interface evolves under reaction conditions by investigating the behavior of nickel electrocatalysts (as nanoparticles and films) deposited on silicon semiconductors. Rigorous electrochemical experiments, interfacial atomistic characterization, and computational modeling are combined to demonstrate critical links between the atomistic features of the interface and the overall performance. It is shown that electrolyte-induced atomistic changes to the interface lead to (1) modulation of the charge carrier fluxes and a dramatic decrease in the electron/hole recombination rates and (2) a change in the barrier height of the interface. Additionally, the critical roles of nonidealities and electrocatalyst coverage due to interfacial geometry are explored. Each of these factors must be considered to optimize the design of metal/semiconductor interfaces which are broadly applicable to photoelectrocatalysis and photovoltaic research.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mathematical Modeling of Hydroxide-Exchange-Membrane Water Electrolyzer

Water electrolyzers can transform intermittent renewable energy like solar energy and wind energy into the chemical energy of hydrogen with zero greenhouse-gas emissions. The hydroxide-exchange membrane electrolyzer (HEME) combines the capability to produce pressurized hydrogen with the advantage of being able to use low or non-platinum group metal (PGM) electrocatalysts in the alkaline environment.1 Hydroxide salts, for example, KOH, are added to the HEME water feed on both anode and cathode to improve its performance. However, the specific mechanism of performance improvement still needs to be further understood. In addition, at high current densities, bubble evolution can result in mass-transport limitations, a less well studied phenomena. Mathematical modeling is ideal to explore these issues as it is cost and time efficient and can deconvolute the physics, processes, and observed phenomena and study the applied-voltage breakdown. In this work, we extend our previously developed 1D two-phase continuum model2 to study the varies processes in the HEME and provide insights on performance optimizations. First, the oxygen evolution reaction (OER) and hydrogen evolution reaction (HER) kinetics at different hydroxide concentrations have been studied by rotating disk electrodes (RDE) and implemented in the model. Then, the model is calibrated and validated against experimental HEME polarization curves for different KOH concentrations as a liquid electrolyte. The model clearly shows a performance increase with increasing KOH concentrations, which is consistent with the experimental results. The reduced ohmic resistance and increased electrochemical active surface area (ECSA) are the two main reasons for performance increase. The large amount of hydroxide in the liquid electrolyte not only helps to distribute the reactant hydroxide throughout the catalyst layer (CL), which reduces ohmic loss, but also enables reaction at the interface between the liquid electrolyte and electrocatalyst, which increases the ECSA. Applied-voltage breakdown demonstrates that the electrolyzer performance is dominated by anode kinetics and ohmic loss. A comparison with the DI water feed shows a more uniform current distribution in the anode CL when KOH is added, which indicates a higher utilization of the CL. Second, we present modeling on the effects of bubble coverage. As gas evolves, part of the ECSA is minimized due to bubble coverage. To account for this effect, an empirical relationship between the fractional bubble coverage and the current density is implemented in the model.3 The model shows this bubble coverage effect is more pronounced at large current densities with DI water feed. Acknowledgements This work was funded under the HydroGEN Consortium by the Energy Efficiency and Renewable Energy, Hydrogen and Fuel Cell Technologies Office, of the U. S. Department of Energy under contract number DE-AC02-05CH11231. References R. Abbasi, B. P. Setzler, S. Lin, J. Wang, Y. Zhao, H. Xu, B. Pivovar, B. Tian, X. Chen, G. Wu and Y. Yan, 31, 1805876 (2019). L. N. Stanislaw, M. R. Gerhardt and A. Z. Weber, ECS Transactions, 92, 767 (2019). H. Vogt and R. J. Balzer, Electrochimica Acta, 50, 2073 (2005).

Liu, Jiangjin↗
DOE OSTI2020

Reduction of uranyl and uranyl-organic complexes mediated by magnetite and ilmenite: A combined electrochemical AFM and DFT study

Reduction of mobile aqueous uranyl complexes can be hindered by the strength of the uranyl-ligand complex and the resistivity of the complexes to electron transfer. Semiconducting solids, such as magnetite (Fe 3 O 4 ) and ilmenite (FeTiO 3 ), can facilitate otherwise-slow electron transfer between the reductant and oxidant, especially if these species are adsorbed on their surfaces as inner-sphere complexes. Electrons can then be shuttled through the surface or the structure. This process can be facilitated by either finding a new electron transfer path or by weakening or destroying the ligand-uranyl complex. Using in situ electrochemical atomic force microscopy (EC-AFM), we imaged reduction products of uranyl as mediated by magnetite and ilmenite in the absence or presence of the organic ligands, ethylenediaminetetraacetate (EDTA) and oxalate. In situ images of uranyl reduction show gradual nucleation and growth of reduced uranium on the mineral surfaces with decreasing reduction potentials. X-ray photoelectron spectroscopy (XPS) measurements indicate pentavalent uranium as the dominant uranium species resulting from uranyl reduction mediated by magnetite and ilmenite. The reduction potentials of uranyl-organic species relevant to our reaction system, (U VI O 2 ∙oxalate) 0 and (U VI O 2 ∙HEDTA) - , are derived using quantum–mechanical modeling and compared with the midpoint potentials obtained from cyclic voltammetry measurements. The spectroscopic, electrochemical, and computational data suggest that the catalysis of Fe-bearing minerals facilitates the reduction of uranyl-organic complexes. U(VI)-oxalate complexes are redox-active and participate directly in the reduction mediated by magnetite. In the presence of EDTA, reduction of uranyl may proceed in concert with ligand exchange between metals (i.e., U and Fe) on the mineral surface. Finally, the proposed mechanism is the breaking of uranyl-organic bonds upon interaction with the mineral surface followed by the formation of Fe-organic bonds, thereby allowing the reduction of hexavalent uranyl to pentavalent one.

58 GEOSCIENCES↗
DOE OSTI2021

Electrochemically Induced Strain Evolution in Pt–Ni Alloy Nanoparticles Observed by Bragg Coherent Diffraction Imaging

Strain is known to enhance the activity of the oxygen reduction reaction in catalytic platinum alloy nanoparticles, whose inactivity is the primary impediment to efficient fuel cells and metal–air batteries. Here, Bragg coherent diffraction imaging (BCDI) was employed to reveal the strain evolution during the voltammetric cycling in Pt–Ni alloy nanoparticles composed of Pt 2 Ni 3 , Pt 1 Ni 1 , and Pt 3 Ni 2 . Analysis of the 3D strain images using a core–shell model shows that the strain as large as 5% is induced on Pt-rich shells due to Ni dissolution. The composition dependency of the strain on the shells is in excellent agreement with that of the catalytic activity. The present study demonstrates that BCDI enables quantitative determination of the strain on alloy nanoparticles during electrochemical reactions, which provides a means to exploit surface strain to design a wide range of electrocatalysts.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Time-resolved electrical potential pump – X-ray photoelectron spectroscopy probe developments for investigating dynamic processes occurring at electrochemical interfaces

Electrode–electrolyte interfaces are of critical importance in several fields, including renewable energy, corrosion, and environmental chemistry. However, investigating these interfaces under operational conditions poses considerable challenges due to the limitations of the instrumentation employed. While recent advancements in in situ and operando techniques have enhanced our comprehension of the steady-state properties of solid-liquid interfaces, the dynamic behaviors of these systems remain inadequately explored. This study introduces a time-resolved X-ray photoelectron spectroscopy (XPS) technique designed to capture transient reaction intermediates and charging dynamics at electrified interfaces. The presented proof-of-principle study demonstrates that electrochemical processes, represented by an equivalent electrical circuit (EEC) model, can be probed and understood using square wave voltage pulses of a potentiostat synchronized to the modified data acquisition of an XPS setup. This method offers a valuable alternative to traditional pump–probe techniques, facilitating the investigation of a broader range of electrochemical systems. A dedicated software package for analyzing time- and energy-resolved XPS with a focus on extracting parameters of the EEC is geared towards benchmarking different EECs in future real-world electrochemical experiments.

Electrochemistry↗
DOE OSTI2020

Operando Revealing Dynamic Reconstruction of NiCo Carbonate Hydroxide for High-Rate Energy Storage

Transition metal carbonate hydroxide (CH) has been widely explored as a promising battery-type electrode for high-rate energy storage. However, its genuine active sites under realistic operating conditions remains elusive. In this work, by virtue of the cutting-edge operando X-ray absorption spectroscopy, we unfold the dynamic evolution of the local electronic and geometric structures for the CH electrode upon cycling. Our results reveal that the irreversible redox of Co cations at the early stage of cyclic voltammetry (CV) cycling induces a phase transition in the NiCo CH nanowires forming oxygen vacancies-enriched NiCo layered double hydroxide (LDH) nanosheets with superior high-rate energy storage ability. Theoretical modeling demonstrates that the unsaturated 5-coordinated Co sites exhibit optimal redox reaction energy barrier and immensely promote the performance. These findings not only elucidate the electrochemical susceptibility of the Co-based carbonate hydroxide under operando conditions, but also highlight an electrochemical coordination-engineering strategy for the rational design of high-performance electrode.

25 ENERGY STORAGE↗
DOE OSTI2020

Engineering Atomically Dispersed FeN 4 Active Sites for CO 2 Electroreduction

Atomically dispersed FeN 4 active sites have exhibited exceptional catalytic activity and selectivity for the electrochemical CO 2 reduction reaction (CO2RR) to CO. However, the understanding behind the intrinsic and morphological factors contributing to the catalytic properties of FeN 4 sites is still lacking. By using a Fe-N-C model catalyst derived from the ZIF-8, we deconvoluted three key morphological and structural elements of FeN 4 sites, including particle sizes of catalysts, Fe content, and Fe-N bond structures. Furthermore, their respective impacts on the CO2RR were comprehensively elucidated. Engineering the particle size and Fe doping is critical to control extrinsic morphological factors of FeN 4 sites for optimal porosity, electrochemically active surface areas, and the graphitization of the carbon support. In contrast, the intrinsic activity of FeN 4 sites was only tunable by varying thermal activation temperatures during the formation of FeN 4 sites, which impacted the length of the Fe-N bonds and the local strains. The structural evolution of Fe-N bonds was examined at the atomic level. First-principles calculations further elucidated the origin of intrinsic activity improvement associated with the optimal local strain of the Fe-N bond.

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

Heterogenous electromediated depolymerization of highly crystalline polyoxymethylene

Abstract Post-consumer plastic waste in the environment has driven the scientific community to develop deconstruction methods that yield valued substances from these synthetic macromolecules. Electrocatalysis is a well-established method for achieving challenging transformations in small molecule synthesis. Here we present the first electro-chemical depolymerization of polyoxymethylene—a highly crystalline engineering thermoplastic (Delrin®)—into its repolymerizable monomer, formaldehyde/1,3,5-trioxane, under ambient conditions. We investigate this electrochemical deconstruction by employing solvent screening, cyclic voltammetry, divided cell studies, electrolysis with redox mediators, small molecule model studies, and control experiments. Our findings determine that the reaction proceeds via a heterogeneous electro-mediated acid depolymerization mechanism. The bifunctional role of the co-solvent 1,1,1,3,3,3-hexafluoro-2-propanol (HFIP) is also revealed. This study demonstrates the potential of electromediated depolymerization serving as an important role in sustainable chemistry by merging the concepts of renewable energy and circular plastic economy.

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

Mapping Interfacial Solution Forces that Drive Fluorescent Nanoparticle Incorporation into Crystals

Interfacial solution structure governs processes ranging from catalytic and electrochemical reactions to particle aggregation and crystallization. Here, we design bio-inspired nanocomposites by mapping interfacial forces and controlling nanoparticle incorporation during crystal growth. Our analytical model contains measured kinetic barriers for surface approaches and equilibrium binding constants with the crystal surface. We validate this model using fluorescent silica nanoparticles and calcite. Our results show that surface chemistry dictates incorporation pathways: methoxy groups are least favorable, hydroxyls reduce kinetic barriers by interacting favorably with hydration layers, carboxylates bind strongly but must overcome a kinetic barrier, and amines outperform all kinetically and thermodynamically. This framework resolves the interplay between interfacial solution structures, particle forces and dynamics, and growth kinetics, enabling predicative control of nanocomposite formation and multiplexed mixed-particle systems.

Zhang, Mingyi [Pacific Northwest National Laborato↗
DOE OSTI2022

Electrochemical CO2 Reduction over Metal-/Nitrogen-Doped Graphene Single-Atom Catalysts Modeled Using the Grand-Canonical Density Functional Theory

Renewably driven, electrochemical conversion of carbon dioxide into value-added products is expected to be a critical tool in global decarbonization. However, theoretical studies based on the computational hydrogen electrode largely ignore the nonlinear effects of the applied potential on the calculated results, leading to inaccurate predictions of catalytic behavior or mechanistic pathways. Here, we use grand canonical density functional theory (GC-DFT) to model electrochemical CO2 reduction (CO2R) over metal- and nitrogen-doped graphene catalysts (MNCs) and explicitly include the effects of the applied potential. We used GC-DFT to compute the CO2 to CO reaction intermediate energies at -0.3, -0.7, and -1.2 VSHE catalyzed by MNCs each doped with 1 of the 10 3d block metals coordinated by four pyridinic nitrogen atoms. Our results predict that Sc-, Ti-, Co-, Cu-, and Zn-N4Cs effectively catalyze CO2R at moderate to large reducing potentials (-0.7 to -1.2 VSHE). ZnN4C is a particularly promising electrocatalyst for CO2R to CO both at low and moderate applied potentials based on our thermodynamic analysis. Our findings also explain the observed pH independence of CO production over FeN4C and predict that the rate-determining step of CO2R over FeN4C is not *CO2- formation but rather *CO desorption. Additionally, the GC-DFT-computed density of states analysis illustrates how the electronic states of MNCs and adsorbates change non-uniformly with applied potential, resulting in a significantly increased *CO2- stability relative to other intermediates and demonstrating that the formation of the adsorbed *CO2- anion is critical to CO2R activation. This work demonstrates how GC-DFT paves the way for physically realistic and accurate theoretical simulations of reacting electrochemical systems.

CO2 reduction↗
DOE OSTI2025

Uncovering the True Active Sites in Ni–N–C Catalysts for CO 2 Electroreduction

Understanding and designing active sites in single-atom catalysts (SACs) requires going beyond static models to capture their dynamic evolution under realistic electrochemical conditions. Here, in this work, we develop an integrated theoretical framework that accounts for operational conditions, by combining grand canonical density functional theory (GC-DFT) with machine-learning-accelerated sampling, to uncover structure–activity–stability relationships in Ni–N–C SACs for the CO 2 reduction reaction (CO 2 RR). A library of NiN x C 4–x (x = 0–4) motifs─representing coordination defects likely formed during high-temperature synthesis─was systematically evaluated. Under working conditions, these sites were found to undergo hydrogenation, and NiN 3 C 1_ H 1 was identified as the most probable active site. At reducing potentials, hydrogen adsorbs spontaneously at C–Ni bridge sites rather than Ni top sites, while subsurface hydrogen facilitates bent CO 2 adsorption crucial for activation. High CO 2 RR selectivity toward CO arises from site separation: Ni centers drive CO2RR, while the hydrogen evolution reaction (HER) occurs at the C–Ni bridge or N sites and from thermodynamic suppression of HER at moderate hydrogen coverage. At more negative potentials, a shift in the CO 2 RR rate-determining process (RDP) and Ni out-of-surface displacement induced by coadsorption of H and H 2 O jointly reduce activity and selectivity. Thus, both the high CO2RR selectivity of Ni–N–C catalysts and its reversal with more negative potentials can be rationalized by accounting for hydrogenated surfaces. This highlights the necessity of modeling realistic; in situ conditions. This framework provides generalizable insights into the dynamic behavior of active sites in SACs, offering guidance for the rational design of active and robust catalysts for a wide range of electrochemical reactions.

25 ENERGY STORAGE↗
DOE OSTI2024

Effect of H + Exchange and Surface Impurities on Bulk and Interfacial Electrochemistry of Garnet Solid Electrolytes

Contact loss and current constriction pose significant challenges at the Li metal interface of solid-state batteries. For garnet-structured Li 7 La 3 Zr 2 O 12 (LLZO), these effects are amplified by Li + /H + exchange and surface contamination reactions, which lead to conductivity losses and poor Li wetting. In this study, we utilize a variety of surface treatment processes across 37 cells to selectively induce proton exchange and contamination reactions in LLZO. The resulting bulk and surface chemistry is systematically characterized and correlated to changes in electrochemical properties. Additionally, we combine impedance analysis and finite element method modeling to deconvolute sources of impedance contributions at the Li metal interface. Specifically, we show that constriction impedance at the Li metal interface arises not solely from voids, but also from ionically-resistive surface contaminants. Further, these findings emphasize the connection between ionic conductivity and constriction, demonstrating that micron-scale ionically-resistive components increase constriction even with identical contact geometries. Finally, we leverage our comprehensive dataset to highlight unstable overpotential growth as a failure mechanism, additionally showing that the phase of a cell’s impedance is a sensitive indicator for the onset of interfacial instability. Overall, this study clarifies the impacts of proton exchange and surface contamination on electrochemical properties at the Li|solid electrolyte interface and elucidates insights that are generalizable to other solid-state battery systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗