Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solid–liquid interfaces”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Electrolyte Organization Leads to Potential-Dependence in Thermochemical Catalysis of Nonpolar Reactions

Electrochemical polarization is now known to play a key role in thermochemical catalysis at solid–liquid interfaces. However, existing frameworks cannot account for why even nonpolar, nonfaradaic reactions are sensitive to interfacial polarization. In order to uncover the molecular basis of this phenomenon, we herein study the potential-dependent reaction kinetics of ethylene and trans-2-butene hydrogenation at Pt–liquid interfaces. Measurements were performed in aqueous and ortho-difluorobenzene (o-DFB) solutions, spontaneously polarizing the Pt–liquid interfaces by, respectively, varying the pH or dissolving distinct metallocene redox buffers into solution. Here, we find that at comparable mechanistic regimes, the rates of both ethylene and trans-2-butene hydrogenation are maximized near the same electrochemical potential, E. Moreover, the potential-dependence, defined as $\frac{∂ln 𝑟}{∂𝐸}$, of trans-2-butene hydrogenation is approximately 2.2× greater than that of ethylene hydrogenation across the full potential range studied. These observations are all consistent with a model in which polarization of the Pt surface away from the local potential of zero free charge (E PZFC ) induces electrostatic organization of the polar solvent and charged ions near the interface, which impedes olefin adsorption and surface reaction because these surface reactions induce electrolyte displacement. Accordingly, interfacial polarization alters the free energy landscape and thus the rate of nonpolar heterogeneous catalysis by controlling the degree of electrostatic organization of polar and charged spectators at the interface, which do not in general need to be specifically chemisorbed onto the surface but could simply be close enough to the surface to be perturbed by the olefin adsorption. These results point toward electrochemical design handles, namely, the electrolyte, catalyst potential, and local E PZFC of the catalyst, with which to tune interfacial catalysis of thermochemical organic transformations.

adsorption↗
DOE OSTI2026

Shining a Light on Some Fundamental Research Opportunities in Semiconductor Photoelectrochemistry

Decades of research in semiconductor photoelectrochemistry have yielded a deep understanding of charge transfer, energetics, and stability at solid−liquid interfaces. Theoretical frameworks developed by Gerischer and contemporaries, together with extensive experimental validation, have clarified the key principles affecting the interfacial kinetics and energetics of semiconductor photoelectrodes. Nevertheless, significant opportunities remain for advances in fundamental understanding of semiconductor photoelectrochemistry. Exciting opportunities include exploiting advances in theory, synthesis, and instrumentation to determine the chemical identity and reactivity of surface states; exerting control of band-edge energetics through molecular-level surface modification processes; and systematically improving emerging photoelectrode protection strategies to enable long-term photoelectrode operation under both oxidative and reductive conditions. Advanced morphologies, such as nanowire and microwire arrays, present new pathways to combine efficient light absorption with effective charge collection and catalyst integration. Unique light−matter interactions during photoelectrochemical deposition of p-type semiconductors readily allow preparation at scale of complex 3D morphologies that are difficult, if not impossible, to access by other methods. Continued exploration of these avenues will expand the fundamental understanding of semiconductor−liquid interfaces and could additionally advance the realization of efficient, stable, and scalable systems for solar fuel generation and other emerging photoelectrochemical applications.

Bean, Paul J. L. [California Institute of Technolo↗
DOE OSTI2025

Quantification and prediction of solidification textures under additive manufacturing conditions

Crystallographic textures are a major determinant of the macroscale anisotropic properties of polycrystalline metallic alloys produced in a wide range of additive manufacturing (AM) processes. Here, we introduce a statistical method that can accurately quantify the degree of orientational order of textures despite the large random fluctuations in the orientation of individual grains inherent in AM processes. The method, demonstrated for laser and resolidification of AlSi thin films, extends Z-scoring to a dynamical regime to assess the statistical significance of observed textures compared to randomly generated ones at different stages of solidification. We further show that, combined with phase-field modeling, this method can be used to infer fundamental anisotropic properties of the solid-liquid interface that are essential for texture prediction, and are compared here to the results of atomistic simulations. In addition, phase-field modeling reveals that, even at rapid AM solidification rates, the observed 〈110〉-dominated textures in the AlSi thin films are controlled predominantly by the anisotropy of the interface free-energy and sheds light on the physical mechanism of grain competition. These results significantly enhance both the existing tools for the quantification and prediction of AM crystallographic textures and our basic understanding of their formation.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Electronic origin of reorganization energy in interfacial electron transfer

Electron transfer (ET) reactions underpin energy conversion and chemical transformations in both biological and a biological systems. The efficiency of any ET process relies on achieving a desired ET rate within an optimal driving force range. Marcus theory provides a microscopic framework for understanding the activation free energy—and therefore the rate—of ET in terms of a key parameter: the reorganization energy. For electrified solid–liquid interfaces, it has long been conventionally understood that only factors in the electrolyte phase are responsible for determining the reorganization energy and that the electronic density of states (DOS) of the electrode only serves to dictate the number of thermally accessible channels for ET. Here we show instead that the electrode DOS plays a central role in governing the reorganization energy, far outweighing its conventionally assumed role. Using atomically layered heterostructures, we tune the DOS of graphene and measure outer-sphere ET kinetics. We find the ensuing variation in ET rate arises from strong modulation in a reorganization energy associated with image potential localization in the electrode. Here we redefine the traditional paradigm of heterogeneous ET kinetics, revealing a deeper role of the electrode electronic structure in interfacial reactivity.

Electrochemistry↗
DOE OSTI2023

Gating ion and fluid transport with chiral solvent

Nonlinear spectroscopy and electrokinetic measurements reveal that the electrochemical properties of a solid–liquid interface in salt solutions in propylene carbonate are dependent on the solvent chirality.

Chemistry↗
DOE OSTI2025

Charges on a suspended silicon nitride membrane under a high-energy electron beam

Thin silicon nitride (SiN x ) membranes are widely used in gas and liquid phase transmission electron microscopy (TEM) and as phase plates to enhance imaging contrast. SiN x contains trap sites for both positive and negative charges, which can be manipulated by high-energy electron irradiation, external potential biasing, or light exposure. Charge accumulation on the membrane can significantly affect in situ TEM processes, including chemical and electrochemical reactions, nanoparticle dynamics, and catalytic activity, or introduce unwanted phase shifts when used as a phase plate. Here, in this study, charge accumulation on suspended SiN x membranes was investigated using off-axis electron holography combined with model-free charge analysis, supported by custom finite element analysis (FEA) simulations. An average residual positive charge density of approximately 2.8 × 10 −4 C m −2 was measured. Localized and stable regions of both positive and negative charges were identified on the membrane. The global positive and localized positive/negative charges give rise to strong electric fields and electroosmotic slip velocities at the membrane surface, which are sufficient to induce non-Brownian particle behavior and directional fluid flow, offering a physical explanation for previously observed anomalies in particle dynamics, nucleation, and growth during gas and liquid phase TEM experiments. These results provide a benchmark for understanding charge behavior at SiNx interfaces in gas and liquid phase TEM. Furthermore, the FEA simulations establish a framework for future investigations into charge distribution, electrostatic potentials, and electrical double layers at solid–liquid interfaces, particularly in complex geometries and chemically dynamic environments.

Suspended silicon↗
DOE OSTI2025

Synthesis challenges, thermodynamic stability, and growth kinetics of La–Si–P ternary compounds

Although many new compounds have been recently predicted with the help of machine learning, the successful experimental synthesis of these compounds remains challenging. Computational insights about the thermodynamic stability and phase formation kinetics among the ground state and competing metastable phases are highly desirable to rationalize and attempt to overcome synthesis challenges experimentally. In this work, we explore synthetic challenges within ternary La–Si–P compounds through feedback between experimental and computational studies. We discuss the experimental challenges in forming three computationally predicted ternary phases (La 2 SiP, La 5 SiP 3 , and La 2 SiP 3 ). To understand the synthetic challenges, we performed molecular dynamics (MD) simulations using an accurate and efficient artificial neural network machine learning (ANN-ML) interatomic potential. We study the phase stability and formation kinetics of these ternary phases in relation to the reported and synthesized La 2 SiP 4 phase. While the growth of the La 2 SiP 4 phase can be reproduced by our MD simulation, our results indicate that the rapid formation of a Si-substituted LaP crystalline phase is a major barrier to the synthesis of the predicted La 2 SiP, La 5 SiP 3 , and La 2 SiP 3 ternary compounds, agreeing well with experimental observations. Our simulations also suggest that there is a narrow temperature window in which the La 2 SiP 3 phase can be grown from the solid–liquid interface.

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

Melting temperature of bismuth to 55 GPa using synchrotron X-ray phase contrast imaging

The melting temperature of elemental bismuth under high pressure has been measured to 55 GPa using synchrotron X-ray phase-contrast imaging in the laser-heated diamond anvil cell. Imaging of solid-liquid interface formation, combined with radiometric temperature and X-ray diffraction measurements, reveals a pronounced reduction in melting boundary slope in Bi-V with pressure. The unusually steep initial slope is attributed to low configurational entropy of melting, arising from structural ordering and coordination matching in the cool liquid, while slope reduction is driven by entropy increase correlated with significant liquid structure changes with rising pressure and temperature. Finally, the data rule out kinetic effects on melting in shock compression experiments and demonstrate the need for improved theoretical phase diagrams.

Materials science↗
DOE OSTI2026

HydraGNN_Predictive_GFM_2026 - Ensemble of predictive graph foundation models for atomistic materials modeling

This release contains data and parameters of HydraGNN-based graph foundation models trained as a result of the work published in the pre-print "Exascale Multi-Task Graph Foundation Models for Imbalanced, Multi-Fidelity Atomistic Data" by M. Lupo Pasini et al. (https://arxiv.org/abs/2604.15380). We jointly train on 16 open first-principles datasets (544+ million structures covering 85+ elements) using a multi-task architecture with per-dataset heads and a scalable ADIOS2/DDStore data pipeline. On Frontier, we execute six large-scale DeepHyper hyperparameter optimization campaigns in FP64 and promote the top-performing message-passing models to sustained 2,048-node training, yielding a PaiNN-based lead model. The version of HydraGNN used to generate the outputs provided in this release is HydraGNN v5.0 (https://github.com/ORNL/HydraGNN/releases/tag/v5.0) The list of datasets used for the training of the graph foundation model is the following: 1) Alexandria [1] 2) ANI1x [2] 3) MPTrj [3] 4) Open Catalyst 2020 (OC20) [4] 5) Open Catalyst 2022 (OC22) [5] 6) Open Catalyst 2025 (OC25) [6] 7) Open Direct ir Capture 2023 (ODAC23) [7] 8) Open Materials 2024 (OMat24) [8] 9) Open Molecules 2025 (OMol25) [9] 10) OMol25-neutral (subset of OMol25 that contains only molecules with zero total charge) 11) OMol25-non-neutral (subset of OMol25 that contains only molecules with non-zero total charge) 12) Open Polymers 2026 (OPoly2026) [10] 13) Nabla2DFT [11] 14) QCML [12] 15) QM7X [reference 13] 16) transition1x [14] Dataset references: [1] J. Schmidt et al., “A dataset of 175k stable and metastable materials calculated with the PBEsol and SCAN functionals,” Scientific Data, vol. 9, p. 64, 2022. [2] J. S. Smith et al., “The ANI-1ccx and ANI-1x data sets, coupled-cluster and density functional theory properties for molecules,” Scientific Data, vol. 7, p. 134, 2020. [Online]. Available: https: //www.nature.com/articles/s41597-020-0473-z [3] A. Jain et al., “Commentary: The Materials Project: A materials genome approach to accelerating materials innovation,” APL Materials, vol. 1, no. 1, p. 011002, 07 2013. [Online]. Available: https://doi.org/10.1063/1.4812323 [4] L. Chanussot et al., “Open catalyst 2020 (oc20) dataset and community challenges,” ACS Catalysis, vol. 11, no. 10, pp. 6059–6072, 2021. [Online]. Available: https://doi.org/10.1021/acscatal.0c04525 [5] K. Tran et al., “Open catalyst 2022 (oc22) dataset and challenges for oxidation electrocatalysts,” ACS Catalysis, vol. 13, no. 5, pp. 3066–3084, 2023. [Online]. Available: https://doi.org/10.1021/acscatal.2c05426 [6] S. J. Sahoo et al., “The open catalyst 2025 (oc25) dataset and models for solid-liquid interfaces,” arXiv preprint arXiv:2509.17862, 2025. [Online]. Available: https://arxiv.org/abs/2509.17862 [7] A. Sriram et al., “The open DAC 2023 dataset and challenges for sorbent discovery in direct air capture,” ACS Central Science, vol. 10, no. 5, pp. 923–941, 2024. [8] L. Barroso-Luque et al., “Open materials 2024 (omat24) inorganic materials dataset and models,” 2024. [Online]. Available: https://arxiv.org/abs/2410.12771 [9] D. S. Levine et al., “The open molecules 2025 (OMol25) dataset, evaluations, and models,” 2025. [Online]. Available: https://arxiv.org/abs/2505.08762 [10] D. S. Levine et al., The open polymers 2026 (OPoly26) dataset and evaluations,” arXiv preprint arXiv:2512.23117, 2025. [Online]. Available: https://arxiv.org/abs/2512.23117 [11] K. Khrabrov et al., “Nabla2dft: A universal quantum chemistry dataset of drug-like molecules and a benchmark for neural network potentials,” in NeurIPS 2024 Datasets and Benchmarks Track, 2024. [Online]. Available: https://openreview.net/forum?id=ElUrNM9U8c [12] S. Ganscha et al., “The QCML dataset, quantum chemistry reference data from 33.5M DFT and 14.7B semi-empirical calculations,” Scientific Data, vol. 12, p. 406, 2025. [13] J. Hoja et al., “QM7-X, a comprehensive dataset of quantum-mechanical properties spanning the chemical space of small organic molecules,” Scientific Data, vol. 8, p. 43, 2021. [Online]. Available: https://www.nature.com/articles/s41597-021-00812-2 [14] M. Schreiner et al., “Transition1x - a dataset for building generalizable reactive machine learning potentials,” Scientific Data, vol. 9, p. 779, 2022. The folder "datasets_ADIOS2_format" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "datasets_ADIOS2_format" directory contains 2 sub-directories, one for the version "v1" of the datasets and one for the version "v2" of the datasets. The version "v1" of the datasets provides values of the total energy as they are extracted from the original data as it was released by the respective institutions. The version "v2" of the datasets provides values of the energy that have been realigned. The realignment was performed by training a linear regression model that predicts the total energy as a function of the chemical composition of the atomistic structure, and then subtract such prediction from the original value of the total energy. Both folders "v1" and "v2" contain 16 sub-directories, each corresponding to an ADIOS2-formatted dataset The folder "DeepHyper-results" contains the configurational files and model's parameters for all the 186 HPO trials that were successfully completed by the scalable hyperparameter optimization (HPO) runs on Frontier. The content of the folder "DeepHyper-results" I structured as follows: 1) task-list.txt: list of mpnn name, jobid, and deephyper task id 2) gfm_${MPNN}_${JOBID}_0.${TASKID}: run directory with checkpoint files 3) gfm_${MPNN}: deephyper summary directory (*.csv) for each specific MPNN type 4) deephyper-experiment-${JOBID}: output and error logs for each job The file "deephyper-sorted.csv" contains the details of each HydraGNN model built and tested by HPO, obtained by merging the (*.csv) filed from each HPO run executed. Out of all the HPO trials, we selected 10 to continue the training of the respective HydraGNN models. Due to limited computational budget available in the LRN070 allocation we could not complete the training till convergence for all these 10 selected models. The folder "models" contains multiple sub-folders, one per each HydraGNN model trained. Each model sub-folder contains the parameters of each HydraGNN model, with multiple checkpoint-restarts. The list of sub-folders are as follows: 1) multidataset_hpo-BEST1-fp64 2) multidataset_hpo-BEST2-fp64 3) multidataset_hpo-BEST3-fp64 4) multidataset_hpo-BEST4-fp64 5) multidataset_hpo-BEST5-fp64 6) multidataset_hpo-BEST6-fp64 7) multidataset_hpo-BEST7-fp64 8) multidataset_hpo-BEST8-fp64 9) multidataset_hpo-BEST9-fp64 10) multidataset_hpo-BEST10-fp64 Within each one of these folders, additional auxiliary log files are provided with descriptions about how the training proceeded. The lead PaiNN-model is contained inside "multidataset_hpo-BEST6-fp64". The file "mlp_branch_weights" contains the parameters of the multi-layer perceptron (MLP) used to reconcile the predictions of the 16 output decoding heads of the HydragNN architectures. The MLP takes in input the chemical composition of the atomistic structure and predicts averaging weights to linearly mix the predictions of each output decoding head toward consolidating them into a single one. The folder "1.1billion-structure-inference" contains 1.1 billion atomistic structures randomly generated. Each structures is associated with energy and forces predicted with the lead-PaiNN model combined with the MLP model for reconciliation of the multi-branch predictions generated by the 16 output decoding heads. The folder "1.1billion-structure-inference" contains 9,300 (*.tar.gz) subdirectories, one per Frontier compute node used to execute the inference at exascale. Once uncompressed, each (*.tar.gz) subdirectory contains an ADIOS2 (*.bp) file container, where each atomistic structure is stored as a PyTorch-Geometric Data object. The file "export_dataset_environment_variables.sh" contains the environment variables that need to be set before running the HydraGNN code to reproduce the results provided in this dataset release. The code that can be used to load the ADIOS2 files, load HydraGNN models, and run inference is available at: https://github.com/ORNL/HydraGNN/releases/tag/v5.0

36 MATERIALS SCIENCE↗
DOE OSTI2025

Controlling Interfacial Energetics and Charge Transfer Rates in 2D Semiconductors: Fundamental Studies en Route to Photoelectrochemical Energy Conversion Beyond the Shockley-Queisser Limit (Final Scientific/Technical Report)

Current photovoltaic and solar-to-fuel technologies do not fully utilize the energy of sunlight because excess photon energy above the semiconductor band gap is rapidly lost as heat through hot-carrier thermalization. Overcoming this loss mechanism is critical, as hot-carrier-based energy conversion systems are predicted to exceed the conventional efficiency limit of ~33%. This project advanced fundamental understanding of hot-carrier energy conversion in two-dimensional (2D) semiconductors, with a focus on monolayer MoS 2 . Using a combination of electrochemical microscopy and in situ ultrafast spectroscopic measurements, this research directly demonstrated hot-carrier extraction from monolayer MoS 2 photoelectrodes in proof-of-concept liquid junction solar cells. These measurements established that hot-carrier transfer can compete with ultrafast carrier cooling at solid–liquid interfaces, providing unambiguous experimental evidence that hot-carrier extraction is feasible in atomically thin semiconductors under operating photoelectrochemical conditions. Beyond demonstration, the project identified design rules for tuning hot-carrier extraction rates relative to cooling rates in 2D semiconductor photoelectrodes. The outcomes of this research provide foundational thermodynamic and kinetic insights for the rational design of next-generation hot-carrier-enabled solar energy conversion systems. These findings have broad implications for photoelectrochemical solar fuels production, electrocatalysis, and emerging energy conversion architectures that seek to harness nonequilibrium charge carriers for enhanced efficiency.

14 SOLAR ENERGY↗
DOE OSTI2026

Final Technical Report for DE-SC0021049: Manipulating interfacial reactivity with atomically layered heterostructures

This final technical report summarizes the work accomplished in this DOE Early Career Research Program project that has established moiré superlattice materials and two-dimensional (2D) heterostructures as a highly tunable platform for controlling heterogeneous charge transfer (ET) kinetics at solid-liquid interfaces. By precisely engineering van der Waals heterostructures of atomically thin 2D materials, particularly bilayer and trilayer graphene with controlled twist angles, this project demonstrated systematic control of interfacial charge transfer rates spanning three orders of magnitude. This research addresses fundamental questions about how electronic structure, charge localization, and atomic layer-dependent properties govern charge transfer at electrochemical interfaces, with broad implications for energy conversion, electrocatalysis, and next-generation electrochemical devices.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Validated Reactive Force Field Quantifies MXene Interfacial Properties, Mechanics, and Thermal Transport

MXenes combine rich surface chemistry, mechanical strength, and high conductivity for a multitude of emerging applications. Predictive modeling supports accelerated materials designs and has been limited by the absence of validated and transferable force fields. Here, we introduce an interpretable, reactive INTERFACE force field (IFF and IFF-R) for Ti 3 C 2 T x MXenes that is trained based on chemical knowledge and achieves quantitative agreement with experiments across lattice parameters (<0.5%), density (<0.2%), liquid contact angles, Raman spectra, and the in-plane elastic modulus (∼320 GPa). The models cover surface terminations from hydroxyl (−OH) to fluorine (−F) groups and are extensible to other chemistries. We introduce pH-resolved surface chemistry and identify dopamine adsorption mechanisms at MXene–aqueous interfaces supported by QCM-D and UV–Vis experiments. The data reveal coplanar and perpendicular binding modes and concentration-dependent multilayer assembly. We predict previously inaccessible properties, including termination-dependent cleavage energies, interlayer shear moduli and dynamic shear failure, nanoindentation and brittle fracture, anisotropic in-plane and out-of-plane thermal conductivities, including the role of defects. Agreement with available experimental data is consistently close and exceeds DFT accuracy across the benchmark properties examined. The IFF/IFF-R model is compatible with CHARMM, AMBER, OPLS, and CVFF force fields for simulations of MXenes with diverse surface terminations, electrolyte interfaces, biointerfaces, and polymer composites without additional parameters. Parameter sets, 3D models, and analysis scripts are provided for community use. The validated, reactive, and transferable IFF framework facilitates predictive design of MXene-based films, membranes, sensing interfaces, and composites.

MXene↗
DOE OSTI2025

Probing the pH Effect on Boehmite Particles in Water Using Vacuum Ultraviolet Single-Photon Ionization Mass Spectrometry

Boehmite has been widely used in theoretical research and industry, especially for hazardous material processing. For the liquid-phase treating process, the interfacial properties of boehmite are believed to be affected by pH conditions, which change its physicochemical behavior. However, molecular-level detection on cluster ions is challenging when using bulk approaches. Herein we employ in situ vacuum ultraviolet single-photon ionization mass spectrometry (VUV SPI-MS) coupled with a vacuum-compatible microreactor system for analysis at the liquid–vacuum interface (SALVI) to investigate the solute molecular composition of boehmite under different pH conditions for the first time. The mass spectral results show that more complex clustering of solute molecules exists at the solid–liquid (s–l) interface than conventionally perceived in a “simple” aqueous solution. Besides solute ions, such as boehmite molecules and fragments, the composition and appearance energies of these newly discovered solvated cluster ions are determined by VUV SPI-MS in different pH solutions. We offer new results for the pH-dependent effect of boehmite and provide insights into a more detailed solvation mechanism at the s–l interface. By comparing the key products under different pH conditions, fundamental understanding of boehmite dissolution is revealed to assist the engineering design of waste processing and storage solutions.

SALVI↗
DOE OSTI2025

A foundation model for atomistic materials chemistry

Atomistic simulations of matter, especially those that leverage first-principles (ab initio) electronic structure theory, provide a microscopic view of the world, underpinning much of our understanding of chemistry and materials science. Over the last decade or so, machine-learned force fields have transformed atomistic modeling by enabling simulations of ab initio quality over unprecedented time and length scales. However, early machine-learning (ML) force fields have largely been limited by (i) the substantial computational and human effort required to develop and validate potentials for each particular system of interest and (ii) a general lack of transferability from one chemical system to the next. Here, we show that it is possible to create a general-purpose atomistic ML model, trained on a public dataset of moderate size, that is capable of running stable molecular dynamics for a wide range of molecules and materials. We demonstrate the power of the MACE-MP-0 model-and its qualitative and at times quantitative accuracy-on a diverse set of problems in the physical sciences, including properties of solids, liquids, gases, chemical reactions, interfaces, and even the dynamics of a small protein. The model can be applied out of the box as a starting or "foundation" model for any atomistic system of interest and, when desired, can be fine-tuned on just a handful of application-specific data points to reach ab initio accuracy. Establishing that a stable force-field model can cover almost all materials changes atomistic modeling in a fundamental way: experienced users obtain reliable results much faster, and beginners face a lower barrier to entry. Foundation models thus represent a step toward democratizing the revolution in atomic-scale modeling that has been brought about by ML force fields.

Batatia, Ilyes↗
DOE OSTI2026

Mechanistic insights into CO 2 capture and electrochemical conversion in nonaqueous Na–CO 2 batteries

Developing efficient energy storage systems that capture and convert CO 2 is critical for mitigating carbon emissions. Here, we report a Na–CO 2 battery with ruthenium dioxide (RuO 2 ) cathode catalysts and propane-1,3-diamine (PDA) as an electrolyte additive to enhance CO 2 capture and conversion efficiency. The integration of CO 2 adsorption and electrochemical reduction facilitates activation of the inert CO 2 molecule and circumvents gas–solid–liquid ternary-phase reactions at the interface. We employed density functional theory (DFT) calculations to systematically unravel the reaction mechanisms and energetics governing CO 2 reduction, both with and without PDA. Our results reveal an energetically favorable pathway toward the formation of Na 2 CO 3 and C as final discharge products, rather than sodium oxalate (Na 2 C 2 O 4 ). The CO 2 –amine adduct facilitates charge transfer from PDA to CO 2 , which results in activation of CO 2 . The kinetics of CO 2 conversion and regeneration of PDA were found to be significantly enhanced on the RuO 2 surface compared to the bulk electrolyte. More importantly, pre-activation of CO 2 via the amine–CO 2 adduct lowers the total overpotential to 2.44 V, compared to 3.13 V without PDA. This study provides fundamental insights into CO 2 electroreduction in Na–CO 2 batteries and underscores the promise of electrolyte engineering for sustainable CO 2 utilization and high-performance energy storage.

25 ENERGY STORAGE↗
DOE OSTI2026

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The Role of Cooperative Interactions Among Surfaces, Solvents, and Reactive Intermediates on Catalysis at Liquid–Solid Interfaces (Final Report DE-SC0020224)

This project established quantitative links between inner‑sphere chemistry (active metal identity, coordination, and zeolite topology) and outer‑sphere organization (solvent identity, hydrogen‑bond networks, and pore condensation) that govern rates, activation barriers, and selectivities for alkene epoxidation and epoxide ring‑opening at solid–liquid and quasi‑liquid–solid interfaces. We deconvoluted contributions from covalent interactions at active sites and noncovalent, solvent‑mediated interactions within pores by pairing well‑defined metal substituted zeolites with controlled solvent environments. We then mapped those contributions onto measurable kinetics (ΔH‡, ΔS‡), adsorption thermodynamics (ITC), and in situ spectroscopy. The transferrable outcomes include a set of design rules that include the following understandings. First, tune silanol ((SiOH)x) density and pore topology to organize solvent networks that selectively stabilize transition states. Second, exploit activity‑coefficient‑normalized rates and adsorption– barrier correlations to diagnose when solvent reorganization rather than surface chemistry limits performance. Third, use partial pore condensation (e.g., acetonitrile, water but also generalizable to other solvents) to elicit liquid‑like stabilization effects even in nominally vapor‑phase reactors. Collectively, these results provide strategies to increase epoxidation rates, improve oxidant utilization (i.e., selectivities), and steer regioselectivity in zeolite‑based catalytic processes relevant to sustainable oxidation chemistry. These outcomes should be transferable to other classes of reactions that proceed in microporous materials and under confinement provided by organized solvents (e.g., electrochemical double layers).

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

Molecular Interlocking Multidimensional Modulations of Cathode‐Electrolyte Interface for Constructing High Energy Density Quasi‐Solid‐State Batteries

Gel polymers are regarded as a promising candidate electrolyte for lithium-metal quasi-solid-state batteries, primarily due to their high ionic conductivity and solid-liquid synergistic properties. However, challenges such as interfacial side reactions, limitations in Li + transport caused by interfacial issues, and leaching of transition metals from the cathode have yet to be effectively solved. Herein, a novel gel electrolyte modulation strategy based on electrostatic filler assembly is proposed to address the issues of ineffective capacity utilization and inadequate cycling stability of high-energy-density cathode materials in solid-state lithium-ion batteries. It constructs a 3D interpenetrating charge-bridge network that effectively tackles the phase-separation challenge between fillers and electrolytes at the molecular level. Meanwhile, the molecular interlocking structure effectively inhibits the electrolyte erosion. More critically, it optimizes and stabilizes the cathode-electrolyte interface film, which facilitates the conduction of Li + -ions through a size-sieving mechanism. Consequently, this strategy enables effective adaptation across diverse high-energy-density cathode materials with satisfactory capacity performance (170.4 mAh g −1 at 4.5 V/1 C for LiNi 0.6 Co 0.2 Mn 0.2 O 2 and 194.0 mAh g −1 at 4.3 V/1 C for LiNi 0.9 Co 0.08 Mn 0.02 O 2 ). In conclusion, this investigation offers a straightforward and effective reference for addressing the critical challenges of ionic transport and interface stabilization in the design of gel electrolytes.

cathode-electrolyte interface↗