Search NASA⌕ Search

SEARCH · Search NASA

Results for “proxy”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

How Does HF-DFT Achieve Chemical Accuracy for Water Clusters?

Bolstered by recent calculations of exact functional-driven errors (FEs) and density-driven errors (DEs) of semilocal density functionals in the water dimer binding energy , we investigate approximate FEs and DEs in neutral water clusters containing up to 20 monomers, charged water clusters, and alkali- and halide-water clusters. Our proxy for the exact density is r 2 SCAN 50, a 50% global hybrid of exact exchange with r 2 SCAN, which may be less correct than r 2 SCAN for the compact water monomer but importantly more correct for long-range electron transfers in the noncompact water clusters. We show that SCAN makes substantially larger FEs for neutral water clusters than r 2 SCAN, while both make essentially the same DEs. Unlike the case for barrier heights, these FEs are small in a relative sense and become large in an absolute sense only due to an increase in cluster size. SCAN@HF, short for SCAN evaluated on the Hartree–Fock (HF) density, produces a cancellation of errors that makes it chemically accurate for predicting the absolute binding energies of water clusters. Likewise, adding a long-range dispersion correction to r 2 SCAN@HF, as in the composite method HF-r 2 SCAN-DC4, makes its FE more negative than in r 2 SCAN@HF, permitting a near-perfect cancellation of FE and DE. r 2 SCAN by itself (and even more so, r 2 SCAN evaluated on the r 2 SCAN 50 density), is almost perfect for the energy differences between water hexamers, and thus probably also for liquid water away from the boiling point. Thus, the accuracy of composite methods like SCAN@HF and HF-r 2 SCAN-DC4 is not due to the HF density being closer to the exact density, but to a compensation of errors from its greater degree of localization. We also give an argument for the approximate reliability of this unconventional error cancellation for diverse molecular properties. Lastly, we confirm this unconventional error cancellation for the SCAN description of the water trimer via Kohn–Sham inversion of the CCSD(T) density.

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

Accelerating Embedding Potential Optimization by Reconstructing the Pseudo-Valence Electron Density

Density functional embedding theory (DFET) enables use of electronic structure methods with higher accuracy than density functional theory in a local region, with applications thus far ranging from (photo/electro)catalysis to reactions in solution. DFET partitions a large collection of atoms into smaller groups that interact via a shared embedding (interaction) potential V emb , determined via functional optimization. The optimized effective potential (OEP) process used to optimize V emb is time-consuming and becomes a computational bottleneck due to sharp, oscillating features of V emb near nuclei. Here, similar to pseudopotential theory, by reconstructing electron densities used in the OEP process from smoother pseudo-valence-only (PVO) electron densities as proxies for total densities of the full system and subsystems, we can retain accuracy in the embedded electronic structure calculations while potentially reducing the overhead of V emb construction, within the projector augmented-wave (PAW) formalism. We explore three different chemical reactions as exemplars to test PVO–DFET, namely, H 2 dissociative adsorption on a Cu(111) surface, H 2 O adsorption on a Pt(111) surface, and aqueous [Ca 2+ –SO 4 2– ] ion-pair formation. The PVO approximation works well for all three systems with minimal loss of accuracy (∼10–70 meV error relative to the original exact-derivative (ED) approach) while accelerating V emb generation for the Cu and Pt systems respectively by 20× and 5×. Given proper numerical convergence parameters, the spatial distributions of differences between PVO- and ED-based V emb outside the core regions are small, explaining the exceptional agreement between the two approaches. Finally, we anticipate that this more efficient PVO–DFET approximation will be useful whenever computation of V emb is much more expensive than subsequent embedded high-level electron correlation calculations.

approximation↗
DOE OSTI2019

Unified Approach to Implicit and Explicit Solvent Simulations of Electrochemical Reaction Energetics

One of the major open challenges in ab initio simulations of the electrochemical interface is the determination of electrochemical barriers under a constant driving force. Existing methods to do so include extrapolation techniques based on fully explicit treatments of the electrolyte, as well as implicit solvent models which allow for a continuous variation in electrolyte charge. Emerging hybrid continuum models have the potential to revolutionize the field, since they account for the electrolyte with little computational cost while retaining some explicit electrolyte, representing a “best of both worlds” method. In this work, we present a unified approach to determine reaction energetics from fully explicit, implicit, and hybrid treatments of the electrolyte based on a new multicapacitor model of the electrochemical interface. A given electrode potential can be achieved by a variety of interfacial structures; a crucial insight from this work is that the effective surface charge gives a good proxy of the local potential, the true driving force of electrochemical processes. In contrast, we show that the traditionally considered work function gives rise to multivalued functions depending on the simulation cell size. Furthermore, we show that the reaction energetics are largely insensitive to the countercharge distribution chosen in hybrid implicit/explicit models, which means that any of the myriad implicit electrolyte models can be equivalently applied. Furthermore, this work thus paves the way for the accurate treatment of ab initio reaction energetics of general surface electrochemical processes using both implicit and explicit electrolytes.

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

Investigating Electron Conductivity Regimes in the Bacterial Cytochrome Wire OmcS

The anaerobic bacterium Geobacter sulfurreducens produces extracellular, electronically conductive cytochrome polymer wires that are conductive over micron length scales. Structure models from cryo-electron microscopy data show OmcS wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. However, the mechanism by which this heme arrangement supports long-range electronic conduction remains unknown. Structure models from cryo-electron microscopy data show these wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. Existing computational models using static heme redox potentials and coupling energies fail to explain experimental observations, predicting conductances 10,000 to 100,000 times lower than measured values. Here, we investigate how dynamic disorder affects site energies, interheme coupling, and long-range electronic conductivity within these cytochrome wires. We introduce an approach to extract charge carrier site information directly from Kohn–Sham density functional theory, without employing projector schemes, and show that site and coupling energies are highly sensitive to changes in interheme geometry and the surrounding electrostatic environment. Unlike models that incorporate dynamic disorder as a thermally averaged quantity, our quantum charge carrier model incorporates proxies for dynamic disorder through decoherence corrections, yielding predicted diffusion coefficient closer to what is expected from experiment and comparable with other organic-based electronic materials. Based on these simulations, we propose that the instantaneous fluctuations of the local electrostatic environment can transiently lift energy degeneracies and delocalize charge carriers. Furthermore, these studies reveal how incorporating dynamic fluctuations associated with the environment resolves the discrepancy between theory and experiment in microbial cytochrome wires and highlight design principles for bioinspired, heme-based conductive materials.

Bioinorganic chemistry↗
DOE OSTI2024

Initial Laboratory Measurements Probing Hydrogen Interactions with Eagle Ford Shale and Pyrite: Potential Implications for Subsurface Hydrogen Storage

Hydrogen (H 2 ) has the potential to be a transformative technology as an enabler to a low-carbon future and promoter of renewable energy. When H 2 is injected and stored in the subsurface, it has the potential to interact with the caprock (usually shale) which overlies and seals the storage reservoir. Here, this study examines geochemical reactions or changes in surface morphology to Eagle Ford Shale, a proxy for caprock, upon exposure to H 2 at 50°C and 10.3 MPa. Reactions were also performed with N 2 to provide an experimental control. Scanning Electron Microscopy with Energy Dispersive Spectroscopy (SEM-EDS), Atomic Force Microscopy (AFM), and Optical Photothermal-Infrared (O-PTIR) Spectroscopy were applied to quantify changes on the microscale and nanoscale level. Fluid chemistry changes were monitored with Ion Chromatography (IC) and Inductively Coupled Plasma Mass Spectrometry (ICP-MS). Exposure of Eagle Ford Shale to H 2 gas alone did not result in any alterations to the shale chemically or any changes in the surface morphology. Exposure of Eagle Ford Shale to both H 2 and water as well as N 2 and water resulted in changes to the surface morphology because of gypsum dissolution and reprecipitation, thus indicating that H 2 is not necessary to promote changes. Pure pyrite was the most reactive with H 2 possibly resulting in a reduction to pyrrhotite. These initial studies suggest that the extent of reactions activated by hydrogen with caprock are minor under the temperature and pressure conditions that would represent underground hydrogen storage.

08 HYDROGEN↗
DOE OSTI2025

Fundamental Insights into Cathode Stability: Linking Compositional Tuning and Local Coordination in Complex Metal Oxides under Aqueous Transformations

Compositional tuning of complex metal oxides in Li-ion battery materials influences their performance as well as their end-of-life behavior, in particular, the tendency to release toxic metal cations in aqueous solution. We modeled ternary variants of a parent LiCoO 2 delafossite structure by varying the metal identity and relative amounts. This yielded ten model formulations of Li(A 4/6 B 1/6 C 1/6 )O 2 , where the material is enriched with the A metal and doped with B and C, with Ni, Mn, Co, Fe, Al, V, and Ti as constituent metals. To assess their stability in aqueous conditions, metal release energetics were calculated using a combination of Density Functional Theory calculations and thermodynamics. Metal release in ternary oxides is dictated by subtle variations in the coordination environment of the leaving group. To identify governing chemical features across diverse compositions with varying local coordination environments, we leverage random forest regression and descriptor importance analysis. A key result is that metal–oxygen orbital hybridization, quantified using a projected density-of-states-derived descriptor, H d/p , provides a physically grounded measure of interaction strength that governs metal release energetics. This refined perspective goes beyond conventional oxidation state considerations and offers more robust insights for materials science. Finally, we model defect surface-bound O 2 dimer formation as a proxy for reactive oxygen species (ROS) generation. The results show that Ni-rich compositions more readily stabilize spin-polarized O 2 dimers, corroborating experimental reports of an increased ROS-driven biological response. In conclusion, our results establish a compositional and electronic basis for metal release and surface oxygen reactivity that form a rationale for complex metal oxide design principles.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Super-Resolved Single-Molecule Tracking Studies of Rhodamine B Accumulation on Fresh and Aged Polyethylene Terephthalate

It is well-known that toxic organic micropollutants (OMs) accumulate on the surfaces of microplastics. However, much remains to be learned about the exact molecular level mechanisms of OM accumulation and how these evolve as the plastics age. In this work, super-resolved single-molecule tracking (SMT) is used for the first time to investigate the accumulation of Rhodamine B (RhB) dye on fresh and artificially aged polyethylene terephthalate (PET) surfaces. PET thin films serve as models for microplastics, while RhB serves as a proxy for the OMs they accumulate. Artificial aging of the films is accomplished by exposing them in a UV-ozone chamber. Water contact angle, spectroscopic ellipsometry, and carbonyl index measurements reveal a gradual decrease in film hydrophobicity, thickness, and carbonyl content with age. Atomic force microscopy (AFM) data reveal an increase in surface roughness and confirm that the films remain largely intact and continuous across the aging times explored. In SMT experiments, wide-field fluorescence videos acquired from the water/PET interface under 7.5 pM RhB reveal both mobile and immobile dye molecules. Measurements of the frame-to-frame displacements of the dye show that diffusion occurs by a desorption-mediated mechanism and that the diffusion rate varies with PET film age. The surface density of mobile dye molecules decreases with increasing PET age, while the population of immobile molecules becomes relatively larger, suggesting an age-dependent transformation of the mechanism(s) by which the dye is accumulated. SMT data reveal that both mobile and immobile molecules repeatedly adsorb over the same surface sites, consistent with the emergence of nanoscale PET surface heterogeneity also revealed by AFM. Estimates of the adsorption coefficients are obtained using a nearest-neighbor analysis, giving values from 9.9 × 10 5 to 2.1 × 10 6 M –1 for immobile molecules and from 1.8 × 10 5 to 2.5 × 10 5 M –1 for mobile molecules on fresh and 5 min aged PET, respectively. Here, anomalous age-dependent variations in the velocity of molecular motion on the PET surface and in the population of immobile molecules are shown to correlate with changes in the strength of RhB adsorption.

Adsorption↗
DOE OSTI2021

Beyond Thermodynamics: Assessing the Dynamical Softness of Hydrated Ions from First Principles

Ion (de)hydration is a key rate-determining step in interfacial processes from corrosion to electrochemical energy storage. However, predicting the kinetics of ion (de)hydration remains challenging, prompting the use of static proxies such as hydration energy and valence. While useful for assessing thermodynamic preferences, such descriptors cannot fully capture the dynamical softness of the hydration shell that dictates kinetics. Accordingly, we use first-principles molecular dynamics to analyze hydration shell softness for a diverse set of metal cations. Additionally, three dynamic metrics are introduced to intuitively describe the bond rigidity, shape deformability, and exchange fluidity of the solvation shell. Together, these metrics capture the relevant physics in the static descriptors, while offering a far more complete and efficient representation for the overall propensity for (de)hydration. Application to the hydrated ion set demonstrates a weak connection between dynamical softness and hydration energy, confirming that dynamical descriptors of hydration are key for correctly describing ion transfer processes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Unconventional Error Cancellation Explains the Success of Hartree–Fock Density Functional Theory for Barrier Heights

Energy barriers, which control the rates of chemical reactions, are seriously underestimated by computationally efficient semilocal approximations for the exchange-correlation energy. The accuracy of a semilocal density functional approximation is strongly boosted for reaction barrier heights by evaluating that approximation non-self-consistently on Hartree–Fock electron densities, which has been known for ~30 years. Here, the conventional explanation is that the Hartree–Fock theory yields the more accurate density. This work presents a benchmark Kohn–Sham inversion of accurate coupled-cluster densities for the reaction H 2 + F → HHF → H + HF and finds a strong, understandable cancellation between positive (excessively overcorrected) density-driven and large negative functional-driven errors (expected from stretched radical bonds in the transition state) within this Hartree–Fock density functional theory. This confirms earlier conclusions (Kaplan, A. D., et al. J. Chem. Theory Comput. 2023, 19, 532–543) based on 76 barrier heights and three less reliable, but less expensive, fully nonlocal density functional proxies for the exact density.

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

Human Coronavirus-229E Hijacks Key Host-Cell RNA-Processing Complexes for Replication

The recent rise in zoonotic coronavirus outbreaks underscores the urgency to understand virus-host interactions and develop potent antiviral therapeutics. Systems biology approaches, particularly proteomics have been invaluable in providing a global overview of such interactions. However, these conventional approaches rely on measuring protein abundance changes which don’t reflect functional shifts. In this study, we employed a high-throughput structural proteomics approach called limited proteolysis-based mass spectrometry (LiP-MS) to capture conformational changes, which we demonstrate are better proxies for functional alterations. We applied this tool to both immortalized and primary human lung cells following human coronavirus 229E (HCoV-229E) infection. We identified significant infection-induced structural changes within RNA processing complexes such as the spliceosome-C and NOP56-associated complex. These observations emphasize that HCoV-229E infection propagates a multi-pronged effort to obstruct the house keeping RNA processing functions in the host. Finally, we show that HCoV-229E replication can be attenuated by the targeted disruption of these complexes, indicating that the identified cellular factories are viable targets to prevent coronavirus infection.

coronavirus↗
DOE OSTI2024

Molecular Mechanisms of Sorbed Ion Effects during Boehmite Particle Aggregation

Classical theories of particle aggregation, such as Derjaguin–Landau–Verwey–Overbeek (DLVO), do not explain recent observations of ion-specific effects or the complex concentration dependence for aggregation. Thus, here, we probe the molecular mechanisms by which selected alkali nitrate ions (Na + , K + , and NO 3 –) influence aggregation of the mineral boehmite (γ-AlOOH) nanoparticles. Nanoparticle aggregation was analyzed using classical molecular dynamics (CMD) simulations coupled with the metadynamics rare event approach for stoichiometric surface terminations of two boehmite crystal faces. Calculated free energy landscapes reveal how electrolyte ions alter aggregation on different crystal faces relative to pure water. Consistent with experimental observations, we find that adding an electrolyte significantly reduces the energy barrier for particle aggregation (~3–4×). However, in this work, we show this is due to the ions disrupting interstitial water networks, and that aggregation between stoichiometric (010) basal–basal surfaces is more favorable than between (001) edge–edge surfaces (~5–6×) due to the higher interfacial water densities on edge surfaces. The interfacial distances in the interlayer between aggregated particles with electrolytes (~5–10 Å) are larger than those in pure water (a few Ångströms). Together, aggregation/disaggregation in salt solutions is predicted to be more reversible due to these lower energy barriers, but there is uncertainty on the magnitudes of the energies that lead to aggregation at the molecular scale. By analyzing the peak water densities of the first monolayer of interstitial water as a proxy for solvent ordering, we find that the extent of solvent ordering likely determines the structures of aggregated states as well as the energy barriers to move between them. Further, the results suggest a path for developing a molecular-level basis to predict the synergies between ions and crystal faces that facilitate aggregation under given solution conditions. Such fundamental understanding could be applied extensively to the aggregation and precipitation utilization in the biological, pharmaceutical, materials design, environmental remediation, and geological regimes.

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

Theory of Cation Solvation in the Helmholtz Layer of Li-Ion Battery Electrolytes

The solvation environments of Li + in conventional nonaqueous battery electrolytes, such as LiPF 6 in mixtures of ethylene carbaronate (EC) and ethyl methyl carbonate (EMC), are often used to rationalize transport properties and solid electrolyte interphase (SEI) formation. Solvation environments in the compact electrical double layer (EDL) next to the electrode, also known as the Helmholtz layer, determine (partially) what species can react to form the SEI, with bulk solvation environments often being used as a proxy. Here, we develop and test a theory of cation solvation in the Helmholtz layer of nonaqueous Li-ion battery electrolytes. First, we validate the theory against bulk and diffuse EDL atomistic molecular dynamics (MD) simulations of LiPF 6 EC/EMC mixtures as a function of surface charge, where we find the theory can qualitatively capture the solvation environments. Next, we turn to the Helmholtz layer, where we find the main effect of the solvation structures next to the electrode is an apparent reduction in the number of binding sites between Li + and the solvents, again where we find reasonable agreement with our developed theory. Finally, by solving a simplified version of the theory, we find that the probability of Li + binding to each solvent remains equal to the bulk probability, suggesting that the bulk solvation environments are a reasonable place to start when understanding battery electrolytes. Our developed formalism can be parametrized from bulk MD simulations and used to predict the solvation environments in the Helmholtz layer through reducing the number of available coordination sites, which can be used to determine what could react and form the SEI.

Helmholtz↗
DOE OSTI2025

Effect of Networking Density on the Patterning Performance of Molecular Layer Deposited Alucone Electron Beam/EUV Resists

One of the main barriers to continued device scaling in the era of extreme ultraviolet (EUV) lithography is the need for improved photoresist chemistries to address challenges such as poor EUV sensitivity, inadequate etch resistance, and pattern collapse. Metal-organic photoresists are a promising class of materials that can address many of these challenges, and among them, resists deposited via hybrid molecular layer deposition (MLD) have attracted interest for their unique advantages in thickness control, chemical homogeneity, and compatibility with vacuum processing. However, despite many successful demonstrations of patterning, little is known about how the molecular design of hybrid MLD resists affects their lithographic performance. In this work, we study the effect of the network structure, a common feature among all hybrid MLD resists, via a series of aluminum alkoxide ("alucone") negative tone resists with varying networking density. Their patterning mechanism is investigated via electron beam lithography (EBL)─a common proxy for EUV─and compared to their EUV-induced reactions studied via flood exposure and in situ characterization. In conclusion, we show that the resist with the least networking density demonstrates the best sensitivity and resolution, with the ability to resolve dense line/space gratings as small as 14 nm half pitch via EBL.

electron beam lithography↗
DOE OSTI2025

Influence of Gel Formation on the Heterogeneous Oxidation of Organic Aerosol

Recent work has shown that divalent cations can perturb the viscosity and phase state of oxygenated organic particles in unexpected ways, sometimes leading to the formation of gel states. Despite the prevalence of divalent cations and oxygenated organic molecules in atmospheric aerosol, particularly in marine environments, the influence of gel states on the reaction rates is not well established. Here, in this work, we measure and compare the impact of viscosity and gel formation on the ozonolysis chemistry of aerosol particles containing the unsaturated, oxygenated organic species ascorbic acid (AA). AA serves as an excellent proxy compound for highly oxidized molecules in the atmosphere. We measured the hygroscopic growth, phase state, viscosity, and water diffusion rate in binary particles containing AA and water and in ternary particles with an additional co-solute of either ammonium sulfate or CaCl 2 . These measurements reveal gel formation in particles mixed with CaCl 2 , as particles become rigid but allow for relatively rapid diffusion rates of water, consistent with previous observations of gel behavior. We measured the rate of ozonolysis of these particles under a range of relative humidity (RH) conditions, showing a clear decrease in the rate with decreasing RH, consistent with the increase in viscosity and slowed diffusion. Notably, we measured the rate of ozonolysis of particles existing in a gel state, showing that following an initial period of decay, the rate of reaction arrests, indicating that the organic material in the solid state is protected from the reaction. This work shows that viscous particles and gel particles exhibit different reactivity toward heterogeneous oxidants, with significant implications for how we understand the chemical evolution of aerosol particles in the atmosphere.

gels↗
DOE OSTI2019

Ab Initio Calculation of Equilibrium Isotopic Fractionations of Potassium and Rubidium in Minerals and Water

Here, we used first-principle approaches to calculate the equilibrium isotopic fractionation factors of potassium (K) and rubidium (Rb) in a variety of minerals of geological relevance (orthoclase, albite, muscovite, illite, sylvite, and phlogopite). We also used molecular dynamics simulations to calculate the equilibrium isotopic fractionation factors of K in water. Our results indicate that K and Rb form bonds of similar strengths and that the ratio between the equilibrium fractionations of K and Rb is approximately 3–4. Under low-temperature conditions relevant to weathering of continents or alteration of seafloor basalts (~25 °C), the K isotopic fractionation between solvated K+ and illite (a proxy for K-bearing clays) is +0.24‰, exceeding the current analytical precision, so equilibrium isotopic fractionation can induce measurable isotopic fractionations for this system at low temperature. These findings, however, cannot easily explain why the δ 41 K value of seawater is shifted by +0.6‰ relative to igneous rocks. Our results indicate that part of the observed fractionation is most likely due to kinetic effects. The narrow range of mean force constants for K and Rb in silicate minerals suggests that phase equilibrium is unlikely to create large K and Rb isotopic fractionations at magmatic temperatures (at least in silicate systems). Kinetic effects associated with diffusion can, however, produce large K and Rb isotopic fractionations in igneous rocks.

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

Source-Resolved Inversion of Elemental Carbon Emissions in California Using Log-Space Bayesian Inference

Elemental carbon (EC), operationally quantified by thermal-optical analysis, is widely used as a proxy for black carbon (BC) relevant to short-term climate forcing and public health. Current EC emission inventories remain highly uncertain, with persistent discrepancies between bottom-up and top-down estimates. In this study, we develop a source-resolved, log-space Bayesian inversion framework applied to estimate California’s statewide EC emissions in 2019. By integrating surface EC measurements from the EPA’s Air Quality System network with high-resolution source contributions simulated by a chemical transport model, we identify a one-third underestimation in the existing statewide EC inventory, requiring an increase of the total from a prior of 8.58 [5.49–13.75] Gg year–1 to a posterior estimate of 12.78 [10.71–15.37] Gg year–1. This discrepancy is primarily driven by substantial underestimations in the power and industrial and off-road mobile sectors. Furthermore, population-weighted exposure analysis reveals a marked sectoral divergence between emission mass and health burden: off-road mobile sources dominate both emissions and exposure, accounting for 31% of statewide exposure, while residential wood combustion contributes 26% of total exposure despite comprising only 19% of total emissions, due to its source proximity to population. These findings underscore the need to update sector-specific EC speciation profiles and demonstrate that mitigation strategies targeting off-road mobile sources and residential wood combustion are critical for reducing EC-related health impacts in California.

Zhang, Jie↗
DOE OSTI2026

Direct Ab Initio Simulation of the Synthesis of BaZrO 3 and the Microstructure Impacts on Proton Transport

Controlling and predicting the processing-structure-performance relationship in functional materials is a grand challenge in materials science, with important implications for a wide range of emerging applications; a high fidelity understanding of the performance impact of microstructures formed under synthesis conditions is required to develop advanced materials, such as solid-state fuel cells and electrolyzers. Using the ceramic BaZrO 3 as a case study, we directly simulate the synthesis and investigate how proton transport is dictated by microstructures. We develop a framework that couples density functional theory (DFT), machine-learning interatomic potential (MLIP) driven molecular dynamics, and grand canonical Monte Carlo to perform large-scale, microstructure-resolved, atomistic simulations of proton transport in experimentally representative polycrystalline structures. Our fully ab initio approach, using a MLIP as a proxy for DFT, allows us to quantify the competition between two distinct diffusion mechanisms: one associated with grain-boundary regions and another within grains. When the impacts of grain boundaries are taken into account, proton transport exhibits substantial deviation from the bulk oxide limit. This addresses long-standing discrepancies between theory and experiments. Our integrated approach provides atomistic insight into microstructure-dependent proton pathways in BaZrO 3 and establishes a general protocol for predicting processing-structure-performance relationships.

organic↗
DOE OSTI2025

Challenges in Product Selectivity for Electrocatalytic Reduction of Amine-Captured CO 2 : Implications for Reactive Carbon Capture

CO 2 is a potential feedstock for carbon-based fuels or materials, but is only available in dilute streams. Integrated processes for CO 2 capture and conversion directly valorize the CO 2 captured by sorbent materials, skipping the energetically expensive sorbent regeneration step. Amines are the most heavily studied liquid-phase sorbent materials for CO 2 capture from dilute streams. Amines react with CO 2 in a 2:1 ratio to form the corresponding ammonium carbamate. Ammonium carbamate [NH 4 ][H 2 NCO 2 ] was tested as the substrate using the highly selective and robust CO 2 -to-formate reduction electrocatalyst [( tBu POCOP)Ir(H)(NCCH 3 ) 2 ], where ( tBu POCOP) is the tridentate pincer ligand 2,6-bis(di tert -butyl-phosphonito). When ammonium carbamate was used as the substrate instead of CO 2 , only hydrogen was produced. An equivalent electrolysis with ammonium hexafluorophosphate with CO 2 also resulted in primarily hydrogen. Methyl carbamate and urea were also tested as substrates as proxies for carbamate that do not contain an equivalent of ammonium, and there was also negligible reduction to carbon-based products. These results indicate that the loss of selectivity observed for aminecaptured CO 2 , or ammonium carbamate, is likely due to the generation of the acidic ammonium equivalent as well as the greater challenge of reducing carbamate compared to CO 2 . This study illustrates that catalysts with high selectivity for concentrated CO 2 can favor hydrogen evolution and loss of carbon-based products when amine-captured CO 2 is used instead.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗