Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Interface correction”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

The Baghdad Atlas: A relational database of inelastic neutron-scattering (n,n ' γ) data

A relational database has been developed based on the original (n,n'γ) work carried out by A. M. Demidov et al., at the Nuclear Research Institute in Baghdad, Iraq (Demidov et al., 1978) for 105 independent measurements comprising 76 elemental samples of natural composition and 29 isotopically-enriched samples. The information from this Atlas includes: γ-ray energies and relative intensities; nuclide and level data corresponding to the residual nucleus and meta data associated with the target sample that allows for the extraction of the flux-weighted (n,n'γ) cross sections for a given transition relative to a defined value. The optimized angular-distribution-corrected fast-neutron flux-weighted partial γ-ray cross section for the production of the 846.8-keV 21+→0gs+γ-ray transition in 56Fe, determined to be $\langle$σγ$\rangle$=143(29) mb, is used for this purpose. However, different values for the adopted cross section can be readily implemented to accommodate user preference based on revised determinations of this quantity. The Atlas (n,n'γ) data has been compiled into a series of CSV-style ASCII data sets and a suite of Python scripts have been developed to build and install the database locally. The database can then be accessed directly through the SQLite engine, or using alternative methods such as the Jupyter Notebook Python-browser interface. Several examples exploiting different interaction methodologies are distributed with the complete software package.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

The sol–gel autocombustion as a route towards highly CO 2 -selective, active and long-term stable Cu/ZrO 2 methanol steam reforming catalysts

The adaption of the sol-gel autocombustion method to the Cu/ZrO 2 system opens new pathways for the specific optimisation of the activity, long-term stability and CO 2 selectivity of methanol steam reforming (MSR) catalysts. Calcination of the same post-combustion precursor at 400 °C, 600 °C or 800 °C allows accessing Cu/ZrO 2 interfaces of metallic Cu with either amorphous, tetragonal or monoclinic ZrO 2 , influencing the CO 2 selectivity and the MSR activity distinctly different. While the CO 2 selectivity is less affected, the impact of the post-combustion calcination temperature on the Cu and ZrO 2 catalyst morphology is more pronounced. A porous and largely amorphous ZrO 2 structure in the sample, characteristic for sol-gel autocombustion processes, is obtained at 400 °C. This directly translates into superior activity and long-term stability in MSR compared to Cu/tetragonal ZrO 2 and Cu/monoclinic ZrO 2 obtained by calcination at 600 °C and 800 °C. The morphology of the latter Cu/ZrO 2 catalysts consists of much larger, agglomerated and non-porous crystalline particles. Based on aberration-corrected electron microscopy, we attribute the beneficial catalytic properties of the Cu/amorphous ZrO 2 material partially to the enhanced sintering resistance of copper particles provided by the porous support morphology.

36 MATERIALS SCIENCE↗
DOE OSTI2021

River and hyporheic zone water geochemical data from OsmoSamplers, East River Pumphouse, Colorado, Nov 2017-Sept 2018

This data package contains geochemical and river water temperature data from OsmoSamplers deployed at East River, Colorado pump house. There are three data files. The first one contains sulfate (SO4) and chloride (Cl) concentrations from riverwater and hyporheic zone (sediment water interface (SWI), 4 cm below seafloor (cmbsf), and 12 (cmbsf)). The second file contains methane concentrations, methane stable carbon isotope ratios and water isotopes from hyporheic zone porewaters. And the third is from a temperature logger during the deployment. Files are provided as .xls or .csv. This research was performed as part of the SBR funded project “Quantifying Subsurface Biogeochemical Variability in a High Altitude Watershed During Winter Isolation” DE-SC0018328. The novel part of this dataset is that we quantify these geochemical parameters throughout the winter isolation period, which is difficult to measure, as well as within the hyporheic zone. The data files were updated 3/1/2021 because the column label for the water isotopes was incorrect. Also, the csv files were also updated since they were not formatted correctly.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Importance of van der Waals interactions for ab initio studies of topological insulators

We investigate the lattice and electronic structures of the bulk and surface of the prototypical layered topological insulators Bi 2 Se 3 and Bi 2 Te 3 using ab initio density functional methods, and systematically compare the results of different methods of including van der Waals (vdW) interactions. We show that the methods utilizing semi-empirical energy corrections yield accurate descriptions of these materials, with the most precise results obtained by properly accounting for the long-range tail of the vdW interactions. The bulk lattice constants, distances between quintuple layers and the Dirac velocity of the topological surface states (TSS) are all in excellent agreement with experiment. In Bi 2 Te 3 , hexagonal warping of the energy dispersion leads to complex spin textures of the TSS at moderate energies, while in Bi 2 Se 3 these states remain almost perfectly helical away from the Dirac point, showing appreciable signs of hexagonal warping at much higher energies, above the minimum of the bulk conduction band. Furthermore, our results establish a framework for unified and systematic self-consistent first principles calculations of topological insulators in bulk, slab and interface geometries, and provides the necessary first step toward ab initio modeling of topological heterostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Isomers and band assignments in the cryogenic vibrational spectrum of the binary complex between water and protonated formic acid using two color, IR–IR photobleaching and H/D isotopic substitution

Protonated formic acid (PFA) is purported to be the active species in the catalytic activation of condensation reactions at the acidic interface of microdroplets. Here, we investigate the fundamental interaction between PFA and water with cryogenic ion vibrational spectroscopy of the binary PFA–H 2 O complexes generated via electrospray ionization followed by buffer gas cooling to about 20 K. The patterns displayed by the isomer-specific IR spectra of D 2 -tagged PFA–H 2 O indicate that two distinct, non-interconverting rotamers are present at low temperatures based on the cis and trans structures of the HCO 2 $H$$^{+}_{2}$ core ion. Both of these occur with the water molecule attached to the OH that is in a cis-configuration relative to the CH group (denoted E), but differ in the E vs Z (cis vs trans relative to the CH group) orientation of the spectator OH. This assignment scheme corrects a previous theoretical analysis that invoked a scenario in which structures with E- and Z-bound water molecules interconvert at low (20 K) temperatures. Isomer-specific bands arising from the OH stretches and water bending modes are deconvoluted using isotopomer-specific spectroscopy of the complexes with partial H/D exchange. The dependence of the nominal shared proton OH stretch frequency on the deuteration of the tethered water confirms strong coupling between this mode and the water bending fundamental.

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

Chemical Compensation Challenges in Processing Antiferroelectric PbZrO 3 Thin Films

The electric field-induced antipolar-to-polar phase transition in antiferroelectric materials is accompanied by large changes in the structural and functional response, making them attractive for many applications ranging from pulsed energy capacitors to high-strain-high (blocking-)force actuators, solid-state heating and cooling systems, and optoelectronic devices. PbZrO 3 , as an end member of the PZT, (x)PbZr 1–x Ti x O 3 , solid solution, has been the subject of many studies. However, processing of perovskite PbZrO 3 is extremely challenging due to phase instabilities induced by Pb loss at the temperatures necessary for perovskite crystallization. Here, we discuss the challenges associated with Pb loss, as well as its compensation through bulk Pb overstoichiometry and interfacial PbO additions, in chemical solution processed, highly oriented PbZrO 3 thin films. For both 042 o - and 001 o -oriented (o-orthorhombic distortion) PbZrO 3 thin films, the crystallization interfaces are the most important contributors to Pb loss and off-stoichiometric outcomes, and no single approach was sufficient to address these challenges. Pb-rich and Pb-deficient nonperovskite phases including ZrO x were only observable through microscopic characterization, and X-ray diffraction alone could not rule out the presence of such secondary phases. “Ideal” macroscopic antiferroelectric polarization-electric field double hysteresis curves were observed despite (at times) the large presence of secondary phases. However, reduced saturation polarization and increased phase transition electric fields were observed concomitantly with the presence of secondary phase(s) and tentatively assigned to the voltage drop across the ZrO x nanocrystals. Based on these results, it is imperative that the presence of secondary phases always be addressed (beyond X-ray diffraction spectra) in order to correctly evaluate the properties of antiferroelectric films.

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

Quantifying transport and electrocatalytic reaction processes in a gastight rotating cylinder electrode reactor via integration of Computational Fluid Dynamics modeling and experiments

Understanding the complexity of the multiple processes of mass, momentum, charge, and heat transport, and how these affect reaction kinetics at the electrode/electrolyte interface is one of the major challenges in the field of energy and catalysis. The rapid and rational scale-up of electrocatalytic systems to industrial scales require a detailed understanding of nonlinear transport-reaction processes, accessible only through the building of multi-physics models that capture with high fidelity the complexity of real-world devices. The gastight rotating cylinder electrode (RCE) reactor is a promising lab-scale tool that can decouple transport from intrinsic kinetics to generate data for first-principle models useful in the design of industrial, electrochemical reactors. Computational Fluid Dynamics (CFD) studies have previously been used to investigate the bulk flow in RCE reactors for simple corrosion and electroplating processes. However, the quantification of changes in local concentration within the viscous layer where catalysis takes place requires capturing the correct flow conditions inside the hydrodynamic boundary layer near the surface of the electrode. Further, this requires simulations with spatial resolution in the nm and μm scale and temporal resolutions between ms and s scales that are similar to the timescales for reactions on the electrode surface. In this study, experimental electrocatalysis is combined with CFD modeling to elucidate and parameterize the hydrodynamics in a gastight RCE reactor. CFD simulations of the electrochemical ferricyanide reduction reaction under mass transport limited conditions are used to evaluate the validity of the CFD model parameters by comparing calculated dimensionless mass transport descriptors to dimensionless correlations obtained experimentally. Justifications for assumptions and details of the simulation methods used in this study are presented to provide a detailed understanding of the effect that each model parameter has on the ability to accurately simulate electrocatalysis in RCE systems. The simulation methodology reported here is a first step towards the development of multi-scale models for the study of transport dependent electrocatalytic processes, such as the electrochemical transformation of CO 2 to fuels and chemicals.

42 ENGINEERING↗
DOE OSTI2026

Alchemy: A Model-Based Approach for 2D to 3D Autonomous Nuclear System Design

Engineering design of nuclear power plant (NPP) piping and equipment systems frequently bypasses crucial 2D system planning, instead moving straight to 3D modeling. This often leads to designs that exceed building envelope constraints, forcing expensive and time-consuming redesigns. When 2D modeling is employed, it typically involves labor-intensive manual workflows that convert 2D drawings into 3D models, resulting in inefficiencies and errors across design iterations. These workflows further suffer from poor software interoperability and dependence on proprietary software ecosystems, thereby contributing to schedule delays and cost overruns. This paper presents Alchemy, an autonomous framework that transforms 2D system definitions into Industry Foundation Classes (IFC)-compliant 3D building information models (BIMs) for expediting nuclear facility design at the conceptual preliminary phase. Using a model-based approach, the framework treats the 2D system diagram as the central reference model employed to automatically generate all subsequent outputs, ensuring consistency between the system definition and the resulting physical design. A web-based interface enables engineers to define hierarchical system topologies including associated equipment, geometric properties, and connectivity requirements. A two-phase equipment layout optimization algorithm automatically computes collision-free spatial configurations within predefined building envelopes. An artificial intelligence (AI)-assisted pipe routing module then generates orthogonal, collision-free routing paths, allowing the user to select either an A* search-based method or an Ant Colony Optimization (ACO)-based method. All outputs are authored natively in IFC format, relying on open-source technologies and standardized formats in order to ensure extensibility and eliminate proprietary software dependencies. The proposed framework is validated on two representative pressurized-water reactor (PWR)-based case studies, for which it autonomously generates IFC-compliant 3D models in minutes, drastically reducing workflows that typically require hours of manual effort. The generated model demonstrates topologically correct equipment placement, physically plausible spatial relationships, and collision-free pipe routing consistent with known PWR loop configurations. This work represents a foundational step toward digital engineering for nuclear facility preliminary design, with future ongoing development targeting design code compliance and expanded system complexity.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Structure of domain walls in chiral spin liquids

The chiral spin liquid is a canonical state of quantum spins combining topological and symmetry-breaking order, and possible experimental realizations have attracted growing interest. We examine the physics at interfaces between chiral spin liquid domains of opposite chirality. We show that a self-consistent mean-field description of spinons remains possible in the vicinity of a domain wall and use this to formulate a Ginzburg–Landau theory of the domain wall. The bulk of a chiral spin liquid contains gapped spinon excitations and gauge fluctuations, set by a finite spinon mass and a nonzero spinon Chern number. A third class of excitations consists of amplitude fluctuations of the spinon hoppings, which admit a geometric interpretation in terms of effective vielbein fields. These fluctuations are usually neglected because they are irrelevant for a homogeneous chiral spin liquid and are suppressed in standard large- N treatments. Going beyond the purely topological Chern–Simons limit, we incorporate these fluctuations into an effective field theoretic framework and show that they generate momentum-dependent corrections, including Chern–Simons-like linear terms and higher-order contributions, beyond the universal topological limit. We then analyze nontopological properties, including domain wall tension and edge velocity, and explain how they modify observables relative to the uniform case. These results connect measurable, nonuniversal quantities such as domain wall width, domain wall tension, and edge velocity to microscopic parameters and provide concrete targets for experiments.

domain walls↗
DOE OSTI2021

Prediction of protein assemblies, the next frontier: The CASP14-CAPRI experiment

Here we present the results for CAPRI Round 50, the fourth joint CASP-CAPRI protein assembly prediction challenge. The Round comprised a total of twelve targets, including six dimers, three trimers, and three higher-order oligomers. Four of these were easy targets, for which good structural templates were available either for the full assembly, or for the main interfaces (of the higher-order oligomers). Eight were difficult targets for which only distantly related templates were found for the individual subunits. Twenty-five CAPRI groups including eight automatic servers submitted ~1250 models per target. Twenty groups including six servers participated in the CAPRI scoring challenge submitted ~190 models per target. The accuracy of the predicted models was evaluated using the classical CAPRI criteria. The prediction performance was measured by a weighted scoring scheme that takes into account the number of models of acceptable quality or higher submitted by each group as part of their five top-ranking models. Compared to the previous CASP-CAPRI challenge, top performing groups submitted such models for a larger fraction (70–75%) of the targets in this Round, but fewer of these models were of high accuracy. Scorer groups achieved stronger performance with more groups submitting correct models for 70–80% of the targets or achieving high accuracy predictions. Servers performed less well in general, except for the MDOCKPP and LZERD servers, who performed on par with human groups. In addition to these results, major advances in methodology are discussed, providing an informative overview of where the prediction of protein assemblies currently stands.

protein complexes↗
DOE OSTI2026

Toward Trustworthy Autonomous Science: A Two-Year Community Roadmap

One year ago, the AISLE roadmap argued that autonomous laboratories operated as isolated islands and proposed a grassroots network organized around five critical dimensions. The field has since moved faster than that roadmap anticipated: multi-agent systems have produced experimentally validated hypotheses, self-driving laboratories have grown more interoperable and orchestrated, reasoning-trained and domain foundation models have raised the capability ceiling, and the Genesis Mission has placed autonomous experimentation at the center of U.S. federal science strategy, with industry emerging as a primary actor. Progress has met a sobering counter-current, including a corrected flagship discovery result, benchmarks showing that agents which rival experts on closed-ended questions still complete only a fraction of open-ended research, and fabricated citations surfacing at leading venues. We read this as the defining tension of the field: producing a candidate discovery is no longer the hard part, but verifying it is, and this asymmetry now limits autonomous science more than raw model capability. Accordingly, we update the roadmap around seven dimensions, revisiting the original five and elevating two former cross-cutting concerns, trust, verification, and reproducibility, and safety, security, and governance, to first-class status. We assess the original milestones (M1 through M14) as achieved, partially achieved, reframed, or open, add four new milestones (M15 through M18) for the elevated dimensions, and scope the path forward to a two-year horizon, with the first year concentrating on interfaces, protocol adoption, and the scaffolding of verification, and the second targeting federation, zero-trust coordination, and governance. Throughout, we position the grassroots network as the interoperability fabric that lets national programs, international initiatives, and commercial platforms connect rather than re-silo.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2026

ChatMPI: LLM-Driven MPI Code Generation for HPC Workloads

The Message Passing Interface (MPI) standard plays a crucial role in enabling scientific applications for parallel computing and is an essential component in high-performance computing (HPC). However, implementing MPI code manually—especially applying a proper domain decomposition and communication pattern—is a challenging and error-prone task. We present ChatMPI, an AI assistant for MPI parallelization of sequential C codes. In our analysis, we focus on testing six essential HPC workloads, which are based on Basic Linear Algebra Subprograms levels 1, 2, and 3 as well as sparse, stencil, and iterative operations. We analyze the process of creating ChatMPI by using the ChatHPC library. This lightweight large language model (LLM)–based infrastructure enables HPC experts to efficiently create and supervise trustworthy AI capabilities for critical HPC software tasks. We study the data required for training (fine-tuning) ChatMPI to generate parallel codes that not only use MPI syntax correctly but also apply HPC techniques to reduce memory communication and maximize performance by using proper work decomposition. With a relatively small training dataset composed of a few dozen prompts and fewer than 15 minutes of fine-tuning on one node equipped with two NVIDIA H100 GPUs, ChatMPI elevates trustworthiness for MPI code generation of current LLMs (e.g., Code Llama, ChatGPT-4o and ChatGPT 5). Additionally, we evaluate the performance of the MPI codes generated by ChatMPI in comparison with the ones generated by ChatGPT-4o and ChatGPT-5. The codes generated by ChatMPI provide up to a 4 × boost in performance by using better problem decomposition, communication patterns, and HPC techniques (e.g., communication avoiding).

Valero Lara, Pedro [ORNL] (ORCID:0000000214794310)↗
DOE OSTI2023

The effect of terrain-influenced winds on fire spread in QUIC-Fire

Here in this manuscript, we describe the implementation of the terrain-following version of QUIC-URB into QUIC-Fire and a demonstration of the impacts of terrain-influenced winds on QUIC-Fire-simulated fire spread. No changes to the underlying QUIC-Fire fire spread algorithm were made other than what was required to correctly account for the inclusion of terrain. This paper summarizes simulations used to understand how the QUIC-URB terrain-influenced winds affect upslope fire behavior without additional changes to the QUIC-Fire fire spread algorithm. Previously published FIRETEC results are compared to simulation results from the modified QUIC-Fire (incorporating terrain-influenced winds) that use the same topographies and fuels. QUIC-Fire results showed overall similar behaviors in terms of how the topographies affected fire shapes and trends in spread rates. Due to the terrain-following version of QUIC-URB being unable to generate flow separations at the crest of hills, fire spread rates in these regions across all non-flat topographies were over-predicted when compared to FIRETEC. Lateral fire growth showed similar trends with FIRETEC between topographies but did not capture the increase in spread due to a diagonal interface between grassland and forested fuel regions in the test domain. These simulations suggest possible refinements that are necessary to improve QUIC-Fire and thus guide ongoing efforts related to: how flame tilt angle is accounted for, the incorporation of non-local drag effects, and the inclusion of the wake-eddy parameterizations that are used in QUIC-URB.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Executive summary of error sources in dynamic surface temperature measurements

Obtaining bulk T requires measuring apparent sample interface T surface , knowledge of any window conditions T window , and knowledge of thermal transport from the sample to the window. To obtain bulk T with uncertainties below 5% requires relatively small uncertainties in each of these areas, and large uncertainties in one area require smaller uncertainties in others to maintain the error budget. For example, if T window is known to 10%, a 5% total uncertainty can be obtained if T surface is known to 1% and combined transport uncertainties are known to 20%. If T surface can only be measured to 2%, combined transport uncertainties must be 17% to reach the same overall uncertainty. As the form of transport is unknown at high pressure, and window temperatures are difficult to measure by their very nature, reducing surface T uncertainties is the practical first step. Here we will discuss various error sources in the measurement of apparent T surface , and which ones must be correctly obtained prior to the experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Neutron and photon in-vivo materials characterization at the evolving plasma-material interface in plasma-burning fusion environments

This report includes the first set of MOSS data collected in the context of fusion-relevant plasma-material interactions. The data collected show that the surface curvature and stress irreversibly change during low-energy deuterium ion irradiation. This holds for both tungsten thin films, bulk tungsten, and single-crystal silicon. Meanwhile, argon irradiation did not produce this same effect. The surface morphology of the exposed samples was virtually identical regardless of gas species. There were, however, some surface chemistry differences between argon-exposed and deuterium-exposed samples, mainly concerning the oxide states of the W4f XPS signal. The changes in surface chemistry could be due to hydrogen dynamics on the tungsten surfaces, though further investigation will be required with in-situ irradiation and analysis. The changes in surface stress could also not be explained by temperature alone, as the temperature reached during the ion irradiations was too low to induce non-elastic strain. With this in mind, deuterium retention in the materials could change the surface curvature and stress. To confirm this, further work will be needed to quantify the change in curvature and relate it to ion flux, fluence, and deuterium retention. Unfortunately, when this work was conducted, the experimental setup for thermal desorption spectroscopy was not operating correctly, and that work could not be included in this dissertation. Despite experimental difficulties, this work shows the first results of a promising technique to study deuterium dynamics in tungsten-based reconstituted thin films.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Unifying principles for catalytic hydrotreating processes (Final Technical Report)

This project builds on the hypothesis that the hydrotreating processes for the removal of oxygen and sulfur are fundamentally similar at the atomic-scale and existing knowledge from the treatment of petroleum derived feedstock can be leveraged for the design of novel catalysts for the upgrade of bio-oil. We tested this hypothesis by comparing computed potential energy diagrams for hydrodesulfurization (HDS) of thiophene over MoS 2 with hydrodeoxygenation (HDO) of furan over MoO 3 and concluded that certain aspects, such as catalyst promotion with transition metals, are valid strategies for both reactions. On the other hand, we also noticed significant differences in the mechanism for hydrogen (H 2 ) activation, which requires sites with metallic character. While MoS 2 is known to have metallic edge states that can catalyze H 2 dissociation, this elementary step is prohibitively slow on defect-free oxides. Only in the presence of vacancies or by creating metal/oxide interfaces can efficient H 2 activation sites during HDO be formed. The need for bifunctional catalyst when it comes to efficient and selective HDO or dehydrogenation reactions was further corroborated in joint experimental and theoretical studies of the Guerbet reaction for the coupling of biomass derived oxygenates over PdCu alloys, nitrate reduction over In-promoted Pd nanoparticles, and ethylene dehydroaromatization over Ga-exchanged ZSM-5 zeolites. All of these catalytic systems have in common that catalytic sites with distinct functional requirements are needed to create a working catalyst. Detailed computational studies were carried out for HDO of m -cresol and phenol on Ru-modified TiO 2 surfaces, which allowed us to attribute catalytic activity to the metal/oxide interface. A surprising finding was that heterolytic cleavage of the H-H bond across the Ru/TiO 2 interface was critically important, despite lower barriers for homolytic H 2 activation on Ru metal. The explanation lies in the high barriers for hydrogen spillover from Ru to TiO 2 , which becomes unnecessary in the heterolytic activation pathway. Moreover, we also reported that proton transfer steps between metal and oxide sites are mediated by weakly adsorbed surface water. During attempts to develop and validate a kinetic Monte Carlo (kMC) model for HDO reactions at the Ru/TiO 2 interface, it became clear that lateral interactions are paramount to describe realistic surface chemistry and without these interactions, the reduction and hydroxylation behavior from our simulations was inconsistent with reported experiments. To assess the importance of lateral interactions in popular computational catalyst design strategies relying on the identification of reactivity descriptors, which can be used along with Brønsted–Evans–Polanyi (BEP) and scaling relations as input to a microkinetic model (MKM) to make predictions for activity or selectivity trends, we compared predicted trends with those obtained from descriptor-based kMC models. We critically evaluated the benefits of kMC over MKM in terms of trend predictions and computational cost when using only a small set of input parameters. After confirming that in the absence of lateral interactions the kMC and MKM approaches yield identical trends and mechanistic information, we observed substantial differences between the two kinetic models when lateral interactions were introduced. The mean-field implementation applies coverage corrections directly to the descriptors, causing an artificial overprediction of the activity of strongly binding metals. In contrast, the cluster expansion in kMC implementation can differentiate among the highly active metals but it is very sensitive to the set of included interaction parameters. Considering that computational screening relies on a minimal set of descriptors, for which MKM makes reasonable trend predictions at a ca. three orders of magnitude lower computational cost than kMC, we concluded that the MKM approach does provide an overall better entry point for computational catalyst design. Overall, this project has led to 11 peer-reviewed publications, and their scientific is impact is well illustrated by their combined 702 citations.

08 HYDROGEN↗
DOE OSTI2024

IDAES-PSE 2.6.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications. IDAES-PSE 2.6.0 Release Highlights Upcoming Changes IDAES will be switching to the new Pyomo solver interface in the next release. Whilst this will hopefully be a smooth transition for most users, there are a few important changes to be aware of. The new solver interface uses a different version of the IPOPT writer (“ipopt_v2”) and thus any custom configuration options you might have set for IPOPT will not carry over and will need to be reset. By default, the new Pyomo linear presolver will be activated with ipopt_v2. Whilst are working to identify any bugs in the presolver, it is possible that some edge cases will remain. IDAES will begin deploying a new set of scaling tools and APIs over the next few releases that make use of the new solver writers. The old scaling tools and APIs will remain for backward compatibility but will begin to be deprecated. New Models, Tools and Features New Intersphinx extension automatically linking Jupyter notebook examples to project documentation New end-to-end diagnostics example demonstrated on a real problem New complementarity formulation for VLE with cubic equations of state, backward compatibility for old formulation New solver interface with presolve (ipopt_v2) in support of upcoming changes to the initialization and APIs methods, with default set to ipopt to maintain backwards compatibility; this will deprecate once all examples have been updated New forecaster and parameterized bidder methods within grid integration library Updated surrogates API and examples to support Keras 3, with backwards compatibility for older formats such as TensorFlow SavedModel (TFSM) Updated costing base dictionary to include the 2023 cost year index value Updated ProcessBlock to include information on the constructing block class Updated Flowsheet Visualizer to allow visualize() method to return value and functions Bug Fixes Fixed bug in the Modular Property Framework that would cause errors when trying to use phase-based material balances with phase equilibria. Fixed bug in Modular Properties Framework that caused errors when initializing models with non-vapor-liquid phase equilibria. Fixed typos flagged by June update to crate-ci/typos and removed DMF-related exceptions Minor corrections of units of measurement handling in power plant waste/transport costing expressions, control volume material holdup expressions, and BTX property package parameters Fixed throwing >7500 numpy deprecation warnings by replacing scalar value assignment with element extraction and item iteration calls Testing and Robustness Migrated slow tests (>10s) to integration, impacting test coverage but also yielding a nearly 30% decrease in local test runtime Pinned pint to avoid issues with older supported Python versions Pinned codecov versions to avoid tokenless upload behavior with latest version Bumped extensions to version 3.4.2 to allow pointing to non-standard install location Deprecations and Removals Python 3.8 is no longer supported. The supported Python versions are 3.9 through 3.12 The Data Management Framework (DMF) is no longer supported. Importing idaes.core.dmf will cause a deprecation warning to be displayed until the next release The SOFC Keras surrogates have been removed. The current version of the SOFC surrogate model in the examples repository is a PySMO Kriging model.

AS↗
DOE OSTI2026

From bulk to surface: Structure and dynamics of amorphous alumina from deep potential molecular dynamics

Understanding the atomic-scale structure and dynamics of amorphous oxide surfaces is essential for interpreting their chemical reactivity, mechanical stability, and interfacial behavior, yet direct experimental characterization remains challenging. We employ Deep Potential (DP) molecular dynamics to generate large-scale, ab initio -quality models of amorphous Al 2 O 3 bulk glasses and melt-quenched free surfaces, enabling a quantitative analysis of both structure and relaxation dynamics with statistical confidence inaccessible to direct ab initio simulation. The trained DP model reproduces experimental liquid and glass structure, captures the cooling-rate dependence of the bulk glass transition, and corrects systematic biases in the polyhedral populations predicted by widely used classical force fields. At the free surface, mass density recovers to bulk values over ~10 Å, while local coordination requires a slightly wider subsurface region to fully converge. The outermost layer is oxygen-enriched, exhibits altered polyhedral connectivity with contracted Al–O bonds, and hosts a broad population of under-coordinated motifs (notably AlO 3 and OAl 2 ) whose abundances are governed by glass stability. These under-coordinated surface motifs exhibit distinct vibrational signatures and occur as locally paired Lewis acid and Brønsted base sites consistent with bond-valence compensation, yet remain spatially dispersed rather than aggregating into extended clusters. Despite this pronounced structural heterogeneity, surface relaxation and the glass-transition temperature remain comparable to their bulk counterparts, suggesting that the disordered surface is kinetically stable once formed. Together, these results establish a molecular-level picture of amorphous alumina surfaces and demonstrate the capability of machine-learned potentials to resolve structure–property relationships in disordered oxide interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗