Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interface capturing”

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.

558 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Operando methods: A new era of electrochemistry

One of the grand challenges facing electrochemistry is to directly resolve the complex nature of (electro)catalyst active sites and capture real-time “movies” of reaction dynamics, i.e. “watching chemistry in action”. The need for such fundamental understanding has stimulated the development of operando/in situ methods, which have greatly enhanced our ability to identify activity descriptors of electrocatalysts and establish structure–property relationships of energy materials. Here, this review summarizes the frontiers of operando electrochemical liquid-cell scanning transmission electron microscopy and correlative synchrotron X-ray methods, which are complementary tools to comprehensively investigate reaction dynamics across multiple spatiotemporal scales. In an effort to encourage greater adoption of advanced operando methods by the general electrochemistry community, this review points out the need to benchmark electrochemistry in confined and heterogenous liquid environment with minimal beam-induced damage. We anticipate that multimodal operando methods will become indispensable for understanding interfacial reaction mechanisms for the broad chemistry and energy materials communities.

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

Continuum-Level Modeling of Li-Ion Battery SEI by Upscaling Atomistically Informed Reaction Mechanisms

Understanding and controlling solid-electrolyte interphase (SEI) formation to stabilize cell performance is a significant challenge for next-generation Li-ion battery technologies. In recent years, computational modeling has become an essential tool in providing fundamental insights into SEI properties and dynamics. However, neither atomistic nor continuum-level approaches alone can capture the complexities of SEI chemistry across all relevant length and time scales. In this work, a continuum-level model is developed that is informed by reaction mechanisms obtained from first-principle calculations. The atomistically informed continuum-level model is used to understand electrolyte degradation, including the decomposition of ethylene carbonate (EC), ethyl methyl carbonate (EMC), and fluoroethylene carbonate (FEC). The model presented here is the most chemically complex continuum-level SEI model in the literature to date. The SEI model is calibrated against experimental irreversible leakage currents and shows qualitative agreement with expected SEI growth trends. The model framework is expected to accelerate fundamental understanding of SEI formation, facilitate mechanism development feedback, and dynamically interact with experimental insights.

continuum-level model↗
DOE OSTI2026

Thermodynamically consistent incorporation of the Langmuir adsorption model into compressible fluctuating hydrodynamics

For a gas–solid interfacial system where chemical species undergo reversible adsorption, we develop a mesoscopic stochastic modeling method that simulates both gas-phase hydrodynamics and surface coverage dynamics by coupling the Langmuir adsorption model with compressible fluctuating hydrodynamics. To this end, we derive a thermodynamically consistent mass–energy update scheme that accounts for how the mass and energy variables in the gas and surface subsystems should be updated according to the changes in the number of molecules of each species in each subsystem due to adsorption and desorption events. By performing a stochastic analysis for the ideal Langmuir model and the full hydrodynamic system, we analytically confirm that our mass–energy update scheme captures thermodynamic equilibrium predicted by equilibrium statistical mechanics. We find that an internal energy correction term is needed, which is attributed to the difference in the mean kinetic energy of gas molecules colliding with the surface from that computed from the Maxwell–Boltzmann distribution. By performing an equilibrium simulation study for an ideal gas mixture of CO and Ar, with CO undergoing reversible adsorption, we validate our overall simulation method and implementation.

Adsorption↗
DOE OSTI2026

Quasi-In-Situ Analysis of Electrode Top Atomic Layers via High-Sensitivity Low-Energy Ion Scattering and Potential-Controlled Sample Transfer

Electrocatalytic reactions involve interfacial interactions between the surfaces of electrodes and reactive species at an electrolyte interface. There are presently no universal or unambiguous methods to directly assay the active top atomic layer composition that influences the reactivity of these electrodes under relevant operating conditions. Low-energy ion scattering (LEIS) spectroscopy is a surface characterization technique that yields compositional analysis of the outermost atomic layer of a material, but it must be performed in ultrahigh vacuum (UHV). Application of LEIS measurements to electrochemical materials that are removed from ambient liquid-phase environments thus leaves an open question as to whether the surface that is transferred to UHV is truly the surface that manifested during the electrochemical reaction. Toward the goal of preserving the active surface state, we developed a sample transfer workflow for LEIS enabling air-free removal and drying of an electrode from an electrochemical cell while maintaining control of the potential using an auxiliary electrode. The potential-controlled emersion method was demonstrated to give distinct potential-dependent surface compositions for a Cu−Pd alloy relative to removal after uncontrolled return to open-circuit potential. A Cu-enriched surface was found at anodic potential and a Pd-enriched surface at cathodic potential, suggesting that the approach can be used to retain representative atomic configurations during transfer. Since adsorbates will often persist from the reaction environment, conventional sample pretreatment methods for removal, including atomic O and atomic H exposure, were also contrasted. Both methods were found to differ with results from incidental low-dose depth profiling by the LEIS primary ion source, which removes adventitious species and surface atoms during the course of repeated measurements. These depth profiles were found to be sensitive to sample history and thus qualitatively informative, despite the possible changes induced by ion damage. The results exhibit (i) the need for complete control over the polarization state of the sample at all times (no excursions to open circuit during transfer) and (ii) the utility of low-dose depth profiling to capture changes in the near-surface composition.

Alloys↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Adatom-Driven Oxygen Intermixing during the Deposition of Oxide Thin Films by Molecular Beam Epitaxy

Thin film deposition from the vapor phase is a complex process involving adatom adsorption, movement, and incorporation into the growing film. While simulations have been able to capture key aspects of these processes to reproduce film morphology, these models typically ignore any subsurface diffusion that occurs. The results of these diffusion and intermixing processes are often observed experimentally, but cannot be explained by invoking bulk lattice diffusion at the deposition temperature. In this work, we present quantitative experimental data that reveals anion and cation intermixing over long length scales during the deposition of epitaxial Fe 2 O 3 and Cr 2 O 3 films and heterostructures by oxygen-plasma-assisted molecular beam epitaxy. We track this diffusion by incorporating well-defined tracer layers containing 18 O and/or 57 Fe, and measure their redistribution on the nanometer scale with atom probe tomography. Molecular dynamics simulations suggest potential intermixing events, which are then examined via nudged elastic band calculations. These reveal that adatoms on the film surface act to “pull up” subsurface O and Fe. Subsequent ring-like rotation mechanisms involving both adatom and subsurface anions then facilitate their mixing. A simple model incorporating these events qualitatively agrees with the experimental intermixing results. In addition to film deposition, these intermixing mechanisms may to be operant during other surface-mediated processes such as heterogeneous catalysis and corrosion.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

High-speed impact of micron-sized diesel drop trains—Splashing dynamics, secondary droplet formation, and effects of pre-existing film thickness

This work investigates surface impingement of monodispersed diesel drop trains using computational simulations. The three-dimensional simulations include high-speed impact of micrometer-sized fuel drops onto initially dry and wetted stainless steel substrates. The drop size and impact velocity are representative of fuel injection conditions in internal combustion engines (ICEs). The drop trains serve as a simplified representation of fuel spray. To accurately capture the film that naturally develops on initially dry substrates, a dynamic diesel-stainless steel contact angle model was developed using experiments of single drop impact under ICE representative temperature. Simulations are presented for a highly splashing case and the effects of a pre-existing film on splashing dynamics were investigated, including the temporal evolution of splashed mass and film thickness. It was concluded that for thin films, the effect of pre-existing film thickness is less significant than in single drop impingement. Secondary droplet characterization was performed on simulation results, leading to instantaneous and time-averaged distributions of secondary droplet size, velocity magnitude, and trajectory angle. It was found that for each drop impingement, approximately 58% of the splashed mass is from the impinging drop itself, while the remainder of the splashed mass is composed of film liquid. From a high resolution case, details of secondary droplet formation are observed and three distinct phases of secondary droplet formation are identified. The detailed analysis of drop train impingement under engine-relevant conditions serves as a first step toward a robust understanding of fuel impingement and the development of cleaner and more efficient ICEs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Dynamic Behavior of Natural Seep Vents: Analysis of Field and Laboratory Observations and Modeling (Final Scientific/Technical Report)

In this project, we have analyzed data collected by the U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) in a high pressure water tunnel (HPWT) and data from two research cruises to natural seeps in the Gulf of Mexico to adapt and validate a numerical model to predict the dynamics of natural seeps in the deep oceans. The HPWT data include video observations of the shrinkage rate of individual methane and natural gas bubbles under simulated deep-water conditions. Field data were collected during two cruises by the Gulf Integrated Spill Research (GISR) Consortium led by Texas A&M University and funded by the Gulf of Mexico Research Initiative (GoMRI). These data included in situ observations from a remotely operated vehicle (ROV) of gas bubbles at two natural seep sites in the Gulf and acoustic observations of the natural seep bubble flares in the ocean water column. The acoustic data were from multibeam echosounders, one mounted in a forward-looking orientation on the ROV and another mounted down-looking in the haul of the ship. All of these laboratory and field data were focused on the dynamics of natural gas bubbles at temperatures and pressures favorable for clathrate hydrate formation between the gas and water. Our analyses of this data focused on understanding the mechanisms responsible for gas bubble dissolution within the hydrate stability zone (HSZ) of the oceans. Ice-like hydrate shells may form on the bubble-water interface under these conditions, and it was unknown how this might affect the mass transfer of gas into the ocean. We were able to extract bubble shrinkage rates from the HPWT datasets. Using this data we determined that mass transfer coefficients with and without a hydrate shell match empirical values for bubbles in contaminated systems (so-called dirty bubbles contaminated by naturally occurring surfactants). We also showed that free gas, and not gas hydrate, is the dominant dissolving phase when the hydrate sub-cooling is below 11 degree Celsius (temperature difference between hydrate the hydrate formation temperature and ambient temperator) or the pressure is reducing as bubbles rise through the ocean water column. Using this mass transfer model, our numerical model of bubble dissolution matched the over 200 HPWT experiments with an average error of 10% for predicting the bubble size at the end of an experiment. From field data in the literature, we also observed that gas bubbles dissolve faster when they are initially released, following mass transfer coefficients for so-called clean-bubbles (those not yet contaminated by surfactants). Shortly after release within the HSZ, a hydrate shell forms on the bubble-water interface, and the mass transfer reduces to rates matching those of dirty bubbles. We correlated this transition time from clean to dirty bubble behavior with the initial bubble surface area and the hydrate sub-cooling. With this model for hydrate formation time and using the mass transfer coefficients deduced from the HPWT data, we validated our numerical model for predicting the rise heights of natural seep flares in the oceans. Flare heights are commonly observed in haul-mounted acoustic multibeam data. The numerical model predicts bubbles to rise high in the ocean water column owing to the slower mass transfer rates for dirty bubbles that accompany the majority of their rise time. We found that the numerical model predictions matched the observed flare heights within 5% to 10% accuracy when we compared the rise heights of the largest bubbles released from the seafloor with the bubbles acoustically visible in the multibeam data. Bubbles become acoustically transparent as they shrink to sizes of order 1 mm in diameter for the multibeam frequencies used in the field. The forward-looking multibeam on the ROV also provided data on the lateral spreading of bubbles in natural seep flares. Our analysis of this data showed that spreading follows a diffusion process, with the effective diffusivity correlating with the wobbling length scale of these ellipsoidal bubbles. When we apply this diffusivity in a random displacement model of bubble spreading, our numerical simulations match closely the lateral spread observed by the M3 in the ocean water column. Finally, we compared the seep model predictions for the acoustic properties of these natural seep plumes with that observed by the acoustic instruments in the field. The M3 and EM 302 observations were converted to relative values of target strength using a calibration we obtained in the laboratory for the M3 and using an algorithm from the manufacturer for the EM 302. Comparing the numerical seep model to these data, we obtain good agreement over the whole height of rise of these bubble flares. This further validates the numerical model. Overall, our validated seep model captures the key dynamics of gas bubbles released from natural seeps in the oceans and helps to predict the fate of methane in the water column.

03 NATURAL GAS↗
DOE OSTI2024

Capturing Surface Coverage Effects in Heterogeneous Catalysis

Adsorbate–adsorbate lateral interactions at relevant surface coverages have a significant effect on chemical kinetics, thereby influencing the activity of a heterogeneous catalyst. Coverage-dependent kinetic and thermodynamic parameters therefore must be included in studies of such complex systems to properly predict the turnover frequencies and kinetic trends. Thus, it becomes extremely important to accurately capture the strength of lateral interactions between neighboring species under realistic reaction conditions. In this Perspective, we discuss the various existing computational and experimental methods for determining adspecies coverage and configurational effects. The choice of the tools and methods employed in such studies depends on factors such as time, length scales, computational cost, the presence of solvents, and reaction conditions. The applications of each method and the respective challenges are also discussed here. As a result, we discuss the recent developments and future of the state-of-the-art for inclusion of surface coverage and configuration into a holistic picture for accurate predictions of catalytic behavior.

09 BIOMASS FUELS↗
DOE OSTI2022

Deciphering the influence of high-flux helium plasma exposure on tungsten nanostructure growth using real-time, in situ spectroscopic ellipsometry

In situ analysis of surfaces during high-flux plasma exposure represents a long-standing challenge in the study of plasma–material interactions. While post-mortem microscopy can provide a detailed picture of structural and compositional changes, in situ techniques can capture the dynamic evolution of the surface. In this study, we demonstrate how spectroscopic ellipsometry can be applied to the real-time characterization of W nanostructure (also known as “fuzz”) growth during exposure to low temperature, high-flux He plasmas. Strikingly, over a wide range of sample temperatures and helium fluences, the measured ellipsometric parameters (Ψ,Δ) collapse onto a single curve that can be directly correlated with surface morphologies characterized by ex situ helium ion microscopy. The initial variation in the (Ψ,Δ) parameters appears to be governed by small changes in surface roughness (<50 nm) produced by helium bubble nucleation and growth, followed by the emergence of 50 nm diameter W tendrils. This basic behavior appears to be reproducible over a wide parameter space, indicating that the spectroscopic ellipsometry may be of general practical use as a diagnostic to study surface morphologies produced by high-flux He implantation in refractory metals. An advantage of the methods outlined here is that they are applicable at low incident ion energies, even below the sputtering threshold. As an example of this application, we apply in situ ellipsometry to examine how W fuzz growth is affected both by varying ion energy and the temperature of the surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Symmetry in Seeded Metal Nanocrystal Growth

Symmetry underpins the structure and function of the world around us and is also captured in modern nanomaterials, where changing the symmetry of a nanocrystal or the interparticle spacing and orientation of nanocrystal building blocks in a superlattice can give new function. However, the synthesis and assembly of nanocrystals have been limited largely to simple compositions and structures. It remains a grand challenge to achieve nanocrystals with compositional and structural complexity while maintaining the monodispersity required for their use. This Account will illustrate through recent examples that seeded methods enable the synthesis of compositionally and structurally complex multimetallic crystals with defined and predictable symmetries for applications in plasmonics and catalysis. This outcome arises because the barrier for heterogeneous nucleation (i.e., seeded) is lower than that of homogeneous nucleation, where seeds can serve as preferential sites for the growth of complex structures and crystal phases. Our analysis begins by considering metal overgrowth from single-crystalline seeds of different shapes and symmetries, where the kinetics of adatom addition to seeds relative to their diffusion across seeds accounts for the expressed nanocrystal shapes. These results are then compared to overgrowth from seeds with different internal structures (i.e., planar defects), where the relationships between nanocrystal size and volumetric strain energy and surface energy are discussed. A major finding from this analysis is that often the underlying symmetry of seeds can be predictably transferred to the final crystals during overgrowth processes. Consequences of this finding are the predictable syntheses of crystals with different hierarchies akin to snowcrystals as well as nanocrystals with complex compositions (e.g., quaternary nanoparticles). Yet, there are subtle aspects to seeded growth that pave a path toward examples where nanocrystal symmetry has been reduced compared to the original seeds in a controlled manner. As we found, both the concentrations of metal precursors and capping agents can impact whether symmetry is transferred or reduced during overgrowth. Examples from our laboratory will be placed in context to other reported strategies for symmetry breaking. As will be argued, understanding what conditions favor symmetry preservation versus symmetry reduction during seeded crystal growth is central to accessing next-generation crystal forms. The Account concludes by outlining synthetic challenges associated with forming nanoscale heterostructures with precise 3-D placement of different materials within a given nanocrystal as well as facet control within different material domain and interface engineering. Furthermore, we envision meeting these challenges through regioselective and chemoselective seeded syntheses for which a foundation is outlined herein.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Fox Trails

1. This software utilizes python pandas to pull data from P6 databases or XER files. The software transforms the datasets into multiple main tables by joining, filtering, iteratively flattening hierarchical structured data, and pivoting datasets to give simple flat output tables. The activity table includes all of the information related to an activity including activity codes, global, EPS, and project codes, UDFs, and WBS information as separate columns. This includes the code id, code value and sequence number for all levels in hierarchical codes. The resource table is similar to the activity table and includes all of the information related to resources on activities including UPFs and resource codes. The resource time phased table takes the resource information and time phases it for the budget, forecast, late, and actual dates/units/costs that closely matches P6's user interface's values as it implements the resource curve and calendars. The wbs table contains the WBS structure broken out by levels and includes UDFs, codes, and notebook topics. The final P6 data table is the relationships table which simply contains the relationships. 2. When a user updates the tool with data (via giving it P6 project names with database username/password information or XER files) the system creates the data in #1, then creates a networkx graph with the activity data imbedded in the node data and the relationships added as edges. Each edge also has it's float calculated (working time distance between the predecessor and successor) and attached to the edge. Activities are also tagged as a potential start of a path based on their constraints, constraint dates, remaining start date, and activity status. When a user enters an activity ID into the UI, it runs a shortest path calculation on the network graph between each node tagged as potential start to the entered activity id based on the float tagged on the edge. Each path returned by the algorithm contains all of the nodes on the path in order, as well as the total float of the edges that make the path. This data is then collected and returned to the user in the form of a gantt chart with groupings for each path that includes the total float for each group. 3. Similar to 2, if the user passes through a reference dataset each activity set in the path is checked to see if it had a path in the reference dataset, if that path was the primary path between the start and end activities, and what has changed regarding logic and durations. These changes are color coded and summarized before sent to the user to be displayed by the UI for simple discovery. 4. Utilizing the data from #1, the user can submit desired grouping code(s) and filters to the system. The system will then pull the activities, resources, and relationships and create a gantt chart based on the groupings sent and filtered based on the filters sent. 5. The system will produce a gantt chart in a similar method to #4, but allows interactivity with the data. As the user interacts with the gantt chart, the software captures the changes and stores it with the user making the change so that project controls and implement those changes in P6.

Fox, Ben↗
DOE OSTI2026

Local Structural Coherence and Interfacial Charge Transfer in Cu 2 ⁢S/Mo⁢S 2 Heterostructure

Precise control over electronic coupling at nanoscale interfaces is critical for designing materials with tunable charge-transfer behavior and catalytic function. Heterostructures with locally coherent interfaces provide a platform for interrogating interfacial charge redistribution in coupled material systems. Here, we report Cu 2 ⁢S/Mo⁢S 2 heterostructures exhibiting nanoscale crystallographic alignment, which are synthesized through a rapid thermal transformation pathway. We employed electrochemical reduction reactions to probe interfacial charge transfer, revealing shifts in product distribution attributable to modified interfacial energetics, even in the absence of optimized catalytic performance. The observed formate Faradaic efficiency suggests that interfacial electronic modulation in the heterostructure shifts product selectivity towards formate, highlighting how interface-driven electronic modulation can direct reaction pathways and influence product selectivity. Optimizing catalyst loading, architecture, and reactor configuration will be critical for future improvement. Structural and compositional integrity were confirmed through powder x-ray diffraction, x-ray photoelectron spectroscopy, and high-resolution transmission electron microscopy. Electron transfer between Cu 2 ⁢S and Mo⁢S 2 domains was further evaluated by electrochemical impedance spectroscopy, while selected-area electron diffraction revealed local crystallographic alignment consistent with a local epitaxial relationship at the Cu 2 S/Mo⁢S 2 heterointerface. To illustrate the broader applicability of this approach, a Zn⁢S/Mo⁢S 2 heterostructure was also synthesized using the same microwave strategy, confirming the generalizability of interfacial engineering principles across metal sulfide-Mo⁢S 2 systems. Collectively, these findings demonstrate that controlled local epitaxial alignment serves as an effective design principle for tuning interfacial energetics and catalytic reactivity in complex heterostructure materials.

carbon capture & utilization↗
DOE OSTI2026

Free Energy and Flexibility Analysis of Autoinhibited Human BRAF

The RAF serine/threonine protein kinases function as direct effectors of RAS in the intracellular transmission of extracellular growth signals, and they are key targets for drug discovery, given the high incidence of oncogenic mutations in RAF and other components of this signaling pathway. In its inactive state, RAF is held in an autoinhibited conformation in the cytosol through a combination of intramolecular interactions and binding to a regulatory 14−3−3 protein dimer. Activation of RAF is initiated by its interaction with membrane-localized GTP-bound RAS, which induces conformational changes that release RAF from its autoinhibited state. However, the molecular mechanisms governing RAF activation remain incomplete, largely due to the challenges in experimentally capturing the intermediate conformational states in this process. To address this gap, we developed a comprehensive all-atom model of BRAF based on existing cryo-EM structures. Using this model, we performed extensive molecular dynamics simulations to evaluate the stability and free energy landscape of autoinhibited BRAF in solution. Our analysis reveals conformational flexibility within the autoinhibited complex, suggesting that this dynamic behavior may play a role in facilitating BRAF activation upon engagement with the membrane-bound RAS.

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

Quantifying the Impact of Metal Population Distribution in MFI-Supported Mo Catalysts for Methane Dehydroaromatization

Precise evaluation of intrinsic kinetic behavior in Mo/MFI catalysts for methane dehydroaromatization (MDA) is confounded by variations in Mo dispersion and speciation, which are influenced by metal loading and zeolite acidity. This work advances the utility of H 2 -temperature programmed reduction (H 2 -TPR) for characterizing Mo/MFI catalysts by enabling quantitative comparison of MoO x populations distinguished by reduction behavior and linked to initial catalytic performance. A comprehensive H 2 reduction pathway is established through systematic H 2 - TPR studies varying catalyst composition (1–10 wt % Mo, Si/Al = 15, 40, ∞), supplemented by UV-Raman spectroscopy, X-ray powder diffraction, N 2 physisorption, NH 3 -TPD, and advanced spectroscopic analysis (in situ XAS with principal component analysis/multivariate curve resolution—alternating least squares). Two low-temperature H 2 -TPR regions capture distinct Mo populations undergoing initial Mo(VI)→Mo(IV) reduction: a lower-temperature population (Mo-RI) associated primarily with highly dispersed, anchored MoO x species expected to predominantly reside within MFI channels, and a higher-temperature population (Mo-RII) corresponding to a broader set of MoO x species that becomes increasingly bulk-like/extrazeolitic at higher Mo loading. Quantification of these populations provides practical, kinetically relevant descriptors for comparing initial MDA rates across catalysts with varying Mo loading, Si/Al ratio, and MoO x heterogeneity. Application of lower-temperature Mo-RI estimates to kinetic measurements reveals a minimum threshold of ∼0.12 × 10 –3 mol Mo-RI/g cat , above which initial forward benzene rates normalized to this population converge despite differences in metal loading and zeolite Brønsted acidity. This threshold coincides with a transition toward a common C 2 -mediated benzene-forming regime, as indicated by approach-to-equilibrium analysis of methane-to-ethane, ethane-to-ethylene, and ethylene-to-benzene reaction steps. Above this threshold, normalized initial benzene rates are nearly invariant with increasing Mo-RII/Mo-RI population ratio, indicating that excess Mo-RII populations, including bulk-like/extrazeolitic MoO x domains present at higher loading, do not measurably suppress benzene formation associated with anchored and mostly channel-confined Mo population under the initial-rate conditions examined.

Mo/MFI↗
DOE OSTI2022

Ba +2 ion trapping using organic submonolayer for ultra-low background neutrinoless double beta detector

If neutrinos are their own antiparticles the otherwise-forbidden nuclear reaction known as neutrinoless double beta decay can occur. The very long lifetime expected for these exceptional events makes its detection a daunting task. In order to conduct an almost background-free experiment, the NEXT collaboration is investigating novel synthetic molecular sensors that may capture the Ba dication produced in the decay of certain Xe isotopes in a high-pressure gas experiment. The use of such molecular detectors immobilized on surfaces must be explored in the ultra-dry environment of a xenon gas chamber. Here, using a combination of highly sensitive surface science techniques in ultra-high vacuum, we demonstrate the possibility of employing the so-called Fluorescent Bicolor Indicator as the molecular component of the sensor. We unravel the ion capture process for these molecular indicators immobilized on a surface and explain the origin of the emission fluorescence shift associated to the ion trapping.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Hydrogen Storage for Flexible Fossil Fuel Power Generation: Integration of Underground Hydrogen Storage with Gas Turbine (Final Report)

As the nation continues to encourage, through market structures and financial incentives, the proliferation of intermittent renewable electricity, how to optimize the ever-changing electric grid and identify means to retain and improve resilience, while ensuring continued reductions in GHG emissions, will be critical. According to Bloomberg, wind & solar generated 10.5% of US electricity in 2020 and that percentage continues to grow. In support of expanding renewable energy use, and to address its intermittent nature, this project will develop the Hydrogen Storage for Flexible Fossil Fuel Power Generation platform that is dispatchable, reliable, repeatable and have the ability to produce zero or negative carbon power while interfacing with geology capable of CO2 and hydrogen storage. GTI Energy (GTIE) and team members Illinois State Geological Survey (ISGS), Mitsubishi Heavy Industries America (MHIA), Ameren Illinois, Hexagon Purus, and the Low Carbon Resources Initiative (LCRI) completed a Phase I Conceptual Study under contract DE-FE0032012 for Hydrogen Storage for Flexible Fossil Fuel Power Generation: Integration of Underground Hydrogen Storage with Gas Turbine. The Hydrogen Storage for Flexible Fossil Fuel Power Generation platform addresses the intermittent nature of the expanding use of Variable Renewable Energy (VRE) generation. The low cost of the electricity (COE) generated results in greater dispatch and more operation at higher power levels (higher efficiency), fewer short intervals, and fewer start/stop cycles. The reliable, resilient system can produce zero carbon power and store hydrogen. It will demonstrate hydrogen storage in geologic formations like those used in natural gas underground storage thus enabling large scale storage of hydrogen in sedimentary strata across the United States rather than in geographically restricted salt caverns. The Phase I study confirmed the system is feasible and generates power at lower cost than other low carbon approaches. The demonstration defines the pathway for broad commercial application and will accelerate the development of larger systems suitable for centralized utility scale electricity production. The study advanced the maturity of the H 2 storage-based system with flexible power generation by completing a Pre-FEED study (Phase II). The Pre-FEED focused on the selected Energy Farm on the University of Illinois Urbana-Champaign (UIUC) site that includes above ground and underground hydrogen storage, low-carbon hydrogen production (GTI’s Compact Hydrogen Generator, CHG) with underground CO2 sequestration, and a 40-MW class gas turbine. The Pre-FEED addressed the entire system and its interconnection to the natural gas and electric grid and mitigation of key risks, such as storage behavior, load-following, and system operation. During Phase 1 of the project, the team completed key tasks, which moved the entire demonstration project, specific components and approaches closer to commercialization. These Phase I Accomplishments include: Completing System Requirements Review; Completing System Layout and Modeling - Heat & Mass Balance and Process Flow Diagram; Completing modelling of 9 turbine performance cases; Evaluating rock strata for underground storage of hydrogen and sequestration of carbon dioxide; Completing initial modelling of underground storage of hydrogen and withdrawal with evaluation of loss and water production; Identifying roadable storage for above ground hydrogen storage; Identifying existing electrical infrastructure for receiving/delivering electricity; Identifying existing gas supply infrastructure for receiving natural gas; Document concept design/development plans in required reports. Conclusions: The 12-month Feasibility study in Phase I study was completed and confirmed the system is feasible and generates power at lower cost than other low carbon approaches and even lower cost than the reference NGCC plant without carbon capture when taking advantage of 45Q carbon credits. The study enabled the fidelity of the concept to be improved and allowed identification of the requirements for the system. Defining the individual system and component requirements was performed via the system requirements review with the whole team. These requirements were then incorporated into and iterated with our Heat & Mass Balance process model and process flow diagrams were generated to reflect the overall system. This information was then used to complete the TEA and show economic feasibility. Large scale non-salt geologic storage of hydrogen is an enabling technology for a hydrogen-fired turbine that can be retrofitted into large-scale electric generating units (EGU). Our demonstration will include 428 MWh or ~4 hours full load of hydrogen storage (above and underground). Carbon capture inherent to the CHG process can capture 90% CO 2 (with upgrades to >98%). This system provides a COE of 23% savings relative to an NGCC with a post combustion amine system. Our proposed storage system decouples carbon capture and hydrogen production from power production; therefore, we expect our proposed system’s efficiency and variable COE to be superior resulting in overall higher dispatch and reduced deep cycling. Our demonstration will be full to multi-day hydrogen storage and has the potential for longer (seasonal) duration commercially. The demonstration defines the pathway for broad commercial application and will accelerate the development of larger systems suitable for centralized utility scale electricity production.

03 NATURAL GAS↗
DOE OSTI2025

Gerischer Electrochemistry Today

Semiconductor photoelectrochemistry is a dynamic and interdisciplinary field at the forefront of research in solar fuels, energy conversion, and catalysis. Here, this Perspective captures the collective insights from the second Gerischer Electrochemistry Today Symposium, held at Colorado State University in Fort Collins, CO, in August 2024, which convened leading researchers, early-career scientists, and industry partners to define the critical next steps for the field. Through interactive sessions, technical talks, panel discussions, and training initiatives─including a Semiconductor Electrochemistry Bootcamp─the symposium emphasized three pillars of advancement: (i) facilitating the exchange of new ideas in semiconductor electrochemistry and charge separation; (ii) fostering the development of future researchers, research topics, and participation in the semiconductor workforce; and (iii) building community. This Energy Focus distills key themes from the meeting and identifies major knowledge gaps in the following areas: mechanisms of charge separation and recombination, role of defects and disorder, dynamic and operando characterization methods, interfacial chemistry and surface passivation, theoretical and modeling limitations, and standardization and benchmarking. The inclusive and collaborative structure of the symposium enabled the generation of this comprehensive report that will serve as a roadmap for fundamental and applied research in the rapidly evolving field of semiconductor electrochemistry over the next decade.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗