Results for “Yield Prediction” — Page 4

DOE OSTI2024

MARVEL Reactor Fuel Performance Report (Rev.2)

The Microreactor Applications Research Validation and EvaLuation (MARVEL) project is producing a high temperature liquid metal-cooled nuclear test bed at Idaho National Laboratory (INL) to ultimately improve the integration of microreactors to end-user applications. This ambitious effort seeks to design, authorize, construct, test, and operate the reactor within five years. In order to construct and operate the MARVEL reactor in a timely manner, the system will utilize materials and component designs which have already been used, qualified, or licensed from previous reactors. The MARVEL reactor will be located at the INL Transient Reactor Test (TREAT) facility in the north high-bay equipment pit and will use the existing 304 stainless steel-clad U-ZrH1.6 pin-type fuel system developed by General Atomics and purchased from TRIGA International. This fuel has been previously qualified under the United States Department of Energy’s (US DOE) Reduced Enrichment for Research and Test Reactors (RERTR) Program. Even though the regulator of the MARVEL reactor is the US DOE, the standards and overall approach recommended by the Nuclear Regulatory Commission is well-defined and utilized here. Following NUREG-1537 regulatory guidance, this report documents the authorization case for the MARVEL fuel system’s application to MARVEL and establishes stable and predictable fuel performance during the most thermophysically unfavorable conditions achievable in the MARVEL reactor. To that end, this report provides a comprehensive survey of the known thermophysical properties, performance, and quantitative relationships associated with the MARVEL reactor fuel element and uses this information to determine its mechanical integrity and risk of reaching unacceptable conditions during the most extreme accident scenarios predicted for the reactor using the most conservative assumptions available. The information contained herein is compiled from a combination of historical reports and peer reviewed scientific publication manuscripts. Known mechanisms under which the fuel is susceptible to failure are highlighted and compared to conditions that could exist in the MARVEL reactor during an unanticipated transient or accident scenario. The two scenarios considered for analysis in this report are (1) an unprotected loss of flow accident and (2) a hypothetical unprotected loss of coolant accident during the loss of flow accident. Preliminary 2D steady-state analyses herein indicate that both fuel-cladding chemical interactions and fuel-cladding mechanical interactions are negligible throughout the fuel’s operational cycle under both normal and high temperature accident scenario conditions. Although higher fidelity 3D time-dependent modeling and simulations are planned, the following may be concluded presently. The MARVEL fuel element maintains its geometric stability and structural integrity during the most extreme accident scenarios predicted for the MARVEL reactor. The hoop stress during the unprotected loss of flow accident reaches about -1.27 MPa; this negative stress indicates that it is compressive rather than tensile. The compressive stress is a result of the NaK pressure on the outside of the fuel element, caused by the restricted thermal expansion of the NaK coolant, exceeding the internal pressures generated inside of the fuel element. The hoop stress generated in the cladding during the unprotected loss of coolant accident reaches a maximum of about approximately 10 MPa, which is nearly an order of magnitude less than the predicted yield strength of the cladding under high-temperature accident scenario conditions. Calculations were compared with results from high performance computational simulations using BISON and are in very close agreement. A conservative MARVEL fuel meat peak temperature limit of 900 °C is recommended presently, which is about 180 °C higher than the peak fuel temperature predicted to occur during the most extreme accident. Based on the known properties and behavior of the MARVEL fuel element, the fuel successfully meets its design and safety requirements under normal and most extreme accident conditions with a large safety margin.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

MARVEL Reactor Fuel Performance Report

The Microreactor Applications Research Validation and EvaLuation (MARVEL) project is producing a high temperature liquid metal-cooled nuclear test bed at Idaho National Laboratory (INL) to ultimately improve the integration of microreactors to end-user applications. This ambitious effort seeks to design, authorize, construct, test, and operate the reactor within five years. In order to construct and operate the MARVEL reactor in a timely manner, the system will utilize materials and component designs which have already been used, qualified, or licensed from previous reactors. The MARVEL reactor will be located at the INL Transient Reactor Test (TREAT) facility in the north high-bay equipment pit and will use the existing 304 stainless steel-clad U-ZrH1.6 pin-type fuel system developed by General Atomics and purchased from TRIGA International. This fuel has been previously qualified under the United States Department of Energy’s (US DOE) Reduced Enrichment for Research and Test Reactors (RERTR) Program. Even though the regulator of the MARVEL reactor is the US DOE, the standards and overall approach recommended by the Nuclear Regulatory Commission is well-defined and utilized here. Following NUREG-1537 regulatory guidance, this report documents the authorization case for the MARVEL fuel system’s application to MARVEL and establishes stable and predictable fuel performance during the most thermophysically unfavorable conditions achievable in the MARVEL reactor. To that end, this report provides a comprehensive survey of the known thermophysical properties, performance, and quantitative relationships associated with the MARVEL reactor fuel element and uses this information to determine its mechanical integrity and risk of reaching unacceptable conditions during the most extreme accident scenarios predicted for the reactor using the most conservative assumptions available. The information contained herein is compiled from a combination of historical reports and peer reviewed scientific publication manuscripts. Known mechanisms under which the fuel is susceptible to failure are highlighted and compared to conditions that could exist in the MARVEL reactor during an unanticipated transient or accident scenario. The two scenarios considered for analysis in this report are (1) an unprotected loss of flow accident and (2) a hypothetical unprotected loss of coolant accident during the loss of flow accident. Preliminary 2D steady-state analyses herein indicate that both fuel-cladding chemical interactions and fuel-cladding mechanical interactions are negligible throughout the fuel’s operational cycle under both normal and high temperature accident scenario conditions. Although higher fidelity 3D time-dependent modeling and simulations are planned, the following may be concluded presently. The MARVEL fuel element maintains its geometric stability and structural integrity during the most extreme accident scenarios predicted for the MARVEL reactor. The hoop stress during the unprotected loss of flow accident reaches about -1.27 MPa; this negative stress indicates that it is compressive rather than tensile. The compressive stress is a result of the NaK pressure on the outside of the fuel element, caused by the restricted thermal expansion of the NaK coolant, exceeding the internal pressures generated inside of the fuel element. The hoop stress generated in the cladding during the unprotected loss of coolant accident reaches a maximum of about approximately 10 MPa, which is nearly an order of magnitude less than the predicted yield strength of the cladding under high-temperature accident scenario conditions. Calculations were compared with results from high performance computational simulations using BISON and are in very close agreement. A conservative MARVEL fuel meat peak temperature limit of 900 °C is recommended presently, which is about 180 °C higher than the peak fuel temperature predicted to occur during the most extreme accident. Based on the known properties and behavior of the MARVEL fuel element, the fuel successfully meets its design and safety requirements under normal and most extreme accident conditions with a large safety margin.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Comparison of interlaminar damage modeling strategies for hybrid composite/aluminum laminates subjected to low-velocity impact

Low-velocity impact of hybrid metal-composite structures was investigated experimentally and computationally. Composite laminates consisting of 2D woven glass fiber reinforced polymer (GFRP) and carbon fiber reinforced polymer (CFRP) were joined with a 6061-T6 aluminum plate using an epoxy adhesive. Two variations of the structure were studied; one consisting of all plies oriented at 0° and one consisting of all plies oriented at 45°. A drop tower was used to impact structures at a range of energies, including energies above and below the threshold at which the aluminum layer was perforated. Numerical simulations were implemented using Sierra/SM, an in-house transient dynamics finite element code developed at Sandia National Laboratories. A Hosford plasticity model was used to describe the response of the aluminum layer. A newly implemented orthotropic continuum damage mechanics (CDM) constitutive model was used to represent the composite laminate. This 3D-CDM model was compared to a cohesive zone model (2D-CDM/CZM) to investigate efficacy of aluminum perforation energy prediction, delamination prediction, and computational cost. Accuracy of each model was evaluated using the experimental results. Each showed good agreement with the tests for both the force and velocity histories, as well as the observed damage mechanisms. The 2D-CDM/CZM model was marginally more accurate in capturing both the composite and aluminum behavior — this model averaged error percentages of -11.2% and 10.8% for residual velocity and peak force, respectively. Meanwhile, the 3D-CDM model predictions yielded average error percentages of -35.5% (velocity) and 22.6% (force). However, the 3D-CDM model generally resulted in a decreased computational cost; the average run time was 14% shorter than the 2D-CDM/CZM model and 3x as many timesteps per hour were computed using the same computational resources. In conclusion, new experimental data on the impact and perforation resistance of metal-composite laminates is presented in addition to numerical predictions of the impact behavior.

Carbon fiber↗

DOE OSTI2022

High Yield, Economical and Environmentally Benign Production of Rare Earth Elements from Coal Ash [Abstract]

Fly ash stored in landfills and ponds across the United States is an attractive, abundant domestic resource for the cost-effective recovery of rare earth elements (REE) and other critical minerals (CM). Physical Sciences Inc. (PSI) and its team members, Winner Water Services (WWS) and University of Kentucky/Center for Applied Energy Research (UK/CAER) successfully executed a multiphase program that developed and demonstrated pilot plant operations for cost-effective and environmentally-friendly production of rare earth element oxide (REO) concentrates, and the critical minerals scandium and aluminum (in the forms of salts or oxide products). We constructed and operated a sub-scale (0.5 kg/day) micropilot facility to demonstrate the key physical and chemical processing operations, predict yields, and troubleshoot process bottlenecks. The project team then designed, constructed and operated two decoupled pilot plant operations: (1) an operational pilot plant for physical separation processes with a capacity of 0.4 metric tons per day (tpd) where we optimized processes to produce selected ash fractions as the feedstock for chemical processing as well as valuable byproducts such as cement substitute, cenospheres, magnetic ash and secondary fuel carbon, and (2) an operational pilot plant for chemical processing with a capacity of 0.5 tpd that developed optimized processes for the production of: (a) REO concentrates, (b) critical minerals recovery (Sc, Al), and (c) beneficiated ash as valuable byproduct suitable for cement applications.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2025

Investigating Electron Conductivity Regimes in the Bacterial Cytochrome Wire OmcS

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

Bioinorganic chemistry↗

DOE OSTI2022

High Yield, Economical and Environmentally Benign Production of Rare Earth Elements from Coal Ash (Phase II Final Summary Report)

Fly ash stored in landfills and ponds across the United States is an attractive, abundant domestic resource for the cost-effective recovery of rare earth elements (REE) and other critical minerals (CM). Physical Sciences Inc. (PSI) and its team members, Winner Water Services (WWS) and University of Kentucky/Center for Applied Energy Research (UK/CAER) successfully executed a multiphase program that developed technologies and their implementation in a pilot plant. We demonstrated plant operations for cost-effective and environmentally-friendly production of rare earth element oxide (REO) concentrates, and the critical minerals scandium and aluminum (in the forms of salts or oxide products), from coal ash. We also constructed and demonstrated a research-scale (0.5 kg/day) micropilot facility to validate the key physical and chemical processing operations, predict yields, and troubleshoot process bottlenecks. The project team then designed, constructed and operated two decoupled pilot plants: (1) an operational pilot plant for physical separation processes with capacity of 0.4 metric tons per day (tpd), where we optimized processes to produce selected ash fractions as the feedstock for chemical processing and as valuable byproducts such as cenospheres, magnetic ash, and secondary fuel carbon, and (2) an operational pilot plant for chemical ash processing with a capacity of 0.5 tpd that developed optimized processes for the production of: (a) REO concentrates, (b) critical minerals (Sc, Al), and (c) beneficiated ash as a valuable byproduct suitable for cement applications. In Phase I, the project team (with Equinox Chemicals in place of WWS) developed and demonstrated the feasibility of the physical and chemical separation processes, developed the design of a pilot plant, and began the development of a preliminary techno-economic model. In the baseline (initial) Phase II program, the project team developed and demonstrated the above pilot scale plant, producing salable REE concentrates, including Y and Sc (REYSc), plus commercially viable byproducts, using environmentally safe and high-yield physical and chemical enrichment processes. The team successfully demonstrated chemical pilot design, construction, shakedown, and operations of the plant. We produced the Phase II deliverable REYSc concentrate ((50 g of >60 wt.% purity REYSc salts on elemental basis), generated the feed for the Phase II follow-on program, identified processing challenges for future optimizations, and refined the techno-economic model. In the Phase II follow-on program, the project team: (1) developed and demonstrated processes to increase the REE amount by 3X (content basis) and convert the Phase II REE salt mixture to an oxide mixture, (2) produced/delivered >38 g of REO mixture with >85 wt.% purity (elemental basis); (3) developed processes to recover critical minerals scandium and aluminum from intermediate streams; (4) produced/delivered > 1 g of scandium salt mixture with >85 wt. % purity (elemental basis); (5) produced/delivered > 100 g of aluminum oxide type material with >70% wt. purity (elemental basis); and (6) updated the techno-economic model from the baseline Phase II program to assess CAPEX and OPEX of a commercial operation. This program has developed extensive databases on process chemistry, unit operations, plant engineering, and techno-economics that will enable further scale-up toward commercial plant design. Specific future developments will be focused on achieving dramatic savings in energy, reagent usage, and operating costs. The combined results will contribute significantly for maturing the technologies of REE recovery from coal byproducts and promote the establishment of domestic REE and CM supply chains.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2021

Quantifying the propagation of parametric uncertainty on flux balance analysis

Flux balance analysis (FBA) and associated techniques operating on stoichiometric genome-scale metabolic models play a central role in quantifying metabolic flows and constraining feasible phenotypes. At the heart of these methods lie two important assumptions: (i) the biomass precursors and energy requirements neither change in response to growth conditions nor environmental/genetic perturbations, and (ii) metabolite production and consumption rates are equal at all times (i.e., steady-state). Despite the stringency of these two assumptions, FBA has been shown to be surprisingly robust at predicting cellular phenotypes. In this paper, we formally assess the impact of these two assumptions on FBA results by quantifying how uncertainty in biomass reaction coefficients, and departures from steady-state due to temporal fluctuations could propagate to FBA results. In the first case, conditional sampling of parameter space is required to re-weigh the biomass reaction so as the molecular weight remains equal to 1 g mmol –1 , and in the second case, metabolite (and elemental) pool conservation must be imposed under temporally varying conditions. Results confirm the importance of enforcing the aforementioned constraints and explain the robustness of FBA biomass yield predictions.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Bayesian discovery of optimal reduced order models from mechanistic and experimental data: A case study of Pd penetration in TRISO fuels using BISON

TRistructural ISOtropic (TRISO) particles rely on a silicon carbide (SiC) layer as the primary structural material and barrier to metallic fission products (FPs) release. Accurate prediction of palladium (Pd) transport and penetration is therefore critical for qualifying TRISO fuels for advanced reactors. The empirical correlation for Pd penetration in BISON is derived from historical particle-fuel data, but cannot explain the large scatter in the experimental data that arises from varying experimental conditions. To aid fuel qualification, we previously developed a mechanistic reduced order model (ROM) using BISON that resolves these dependencies. Here, in this work we build on that mechanistic ROM and perform validation and quantify its uncertainty using Bayesian uncertainty quantification (UQ). calibration against a suite of in-pile and out-of-pile experiments spanning particle compositions, geometries, and operating conditions, and we benchmark it against the empirical correlation. Bayesian UQ identifies influential parameters, calibrates them to data, and yields predictive intervals. Results show that while the empirical correlation can be tuned to fit a single experiment type, it transfers poorly; the mechanistic ROM sustains accuracy with credible uncertainty across disparate conditions. This demonstrates a practical path—via Bayesian UQ applied to mechanistic ROMs—to leverage single-effect experiments for inferring in-reactor behavior and supporting TRISO fuel qualification.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Impact of duration and missing data on the long-term photovoltaic degradation rate estimation

Accurate quantification of photovoltaic (PV) system degradation rate (R D ) is essential for lifetime yield predictions. Although R D is a critical parameter, its estimation lacks a standardized methodology that can be applied on outdoor field data. The purpose of this paper is to investigate the impact of time period duration and missing data on R D by analyzing the performance of different techniques applied to synthetic PV system data at different linear R D patterns and known noise conditions. The analysis includes the application of different techniques to a 10-year synthetic dataset of a crystalline Silicon PV system, with emulated degradation levels and imputed missing data. Here, the analysis demonstrated that the accuracy of ordinary least squares (OLS), year-on-year (YOY), autoregressive integrated moving average (ARIMA) and robust principal component analysis (RPCA) techniques is affected by the evaluation duration with all techniques converging to lower R D deviations over the 10-year evaluation, apart from RPCA at high degradation levels. Moreover, the estimated R D is strongly affected by the amount of missing data. Filtering out the corrupted data yielded more accurate R D results for all techniques. It is proven that the application of a change-point detection stage is necessary and guidelines for accurate R D estimation are provided.

14 SOLAR ENERGY↗

DOE OSTI2020

Determination of the strong coupling constant $\alpha_{S}(m_\mathrm{Z})$ from measurements of inclusive W$^\pm$ and Z boson production cross sections in proton-proton collisions at $ \sqrt{\mathrm{s}} $ = 7 and 8 TeV

Twelve measurements of inclusive cross sections of W$^{±}$ and Z boson production, performed in proton-proton collisions at centre-of-mass energies of 7 and 8 TeV, are compared with perturbative quantum chromodynamics calculations at next-to-next- to-leading order (NNLO) accuracy obtained with the CT14, HERAPDF2.0, MMHT14, and NNPDF3.0 parton distribution functions (PDFs). Data and theory agree well for all PDF sets, taking into account the experimental and theoretical uncertainties. A novel procedure is employed to extract the strong coupling constant at the Z pole mass from a detailed comparison of all the experimental fiducial cross sections to the corresponding NNLO theoretical predictions, yielding $ {\alpha}_{\mathrm{S}}\left({m}_{\mathrm{Z}}\right)={0.1163}_{-0.0031}^{+0.0024} $ (CT14), $ {0.1072}_{-0.0040}^{+0.0043} $ (HERAPDF2.0), 0.1186 ± 0.0025 (MMHT14), and 0.1147 ± 0.0023 (NNPDF3.0). Using the results obtained with the CT14 and MMHT14 PDFs, which yield the most robust and stable α$_{S}$(m$_{Z}$) extractions, a value $ {\alpha}_{\mathrm{S}}\left({m}_{\mathrm{Z}}\right)={0.1175}_{-0.0028}^{+0.0025} $ is determined.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Using in-situ strain measurements to evaluate the accuracy of stress estimation procedures from fracture injection/shut-in tests

Fracture injection/shut-in tests are commonly used to measure the state of stress. Injection creates a hydraulic fracture (or in some cases, opens a preexisting fracture), and then the pressure after shut-in is monitored to identify fracture closure. Different interpretation procedures have been proposed for estimating closure, and the procedures sometimes yield significantly different results. In this study, direct, in-situ strain measurements are used to observe fracture reopening and closure. The tests were performed as part of the EGS Collab project, a mesoscale project performed at 1.25 and 1.5 km depth at the Sanford Underground Research Facility. The tests were instrumented with the SIMFIP tool, a double-packer probe with a high-resolution three-dimensional borehole displacement sensor. The measurements provide a direct observation of the fracture closure signature, enabling a high-fidelity estimate of the fracture closure stress (ie, the normal stress on the fracture). In two of the four tests, injection created an opening mode fracture, and so the closure stress can be interpreted as the minimum principal stress. In the other two tests, injection probably opened preexisting natural fractures, and so the closure stress can be interpreted as the normal stress on the fractures. The strain measurements are compared against different proposed methods for estimating closure stress from pressure transients. The shut-in transients are analyzed with two techniques that are widely used in the field of petroleum engineering – the ‘tangent’ method and the ‘compliance’ method. In three of the four tests, the tangent method significantly underestimates the closure stress. The compliance method is reasonably accurate in all four tests. Closure stress is also interpreted using two other commonly-used methods – ‘first deviation from linearity’ and the method of Hayashi and Haimson. In comparison with the SIMFIP data, these methods tend to overestimate the closure stress, evidently because they identify closure from early-time transient effects, such as near-wellbore tortuosity. In two of the tests, microseismic imaging provides an independent estimate of the size of the fracture created by injection. When combined with a simple mass balance calculation, the SIMFIP stress measurements yield predictions of fracture size that are reasonably consistent with the estimates from microseismic. In conclusion, the calculations imply an apparent fracture toughness 2-3x higher than typical laboratory-derived values.

58 GEOSCIENCES↗

DOE OSTI2025

Mechanistic and kinetic relevance of hydrogen and water in CO 2 hydrogenation on Cu-based catalysts

Here, we ally steady-state kinetics, kinetic isotope effects, and density functional theory (DFT) calculations to illustrate that Cu-based catalysts remain saturated by H-adatoms (H*) and molecular formic acid (HCOOH**) during CO 2 hydrogenation. High H* coverage under methanol synthesis conditions is evidenced by reverse water-gas shift (RWGS) rates that exhibit positive H 2 reaction orders only at P H2 ≲ 0.5 bar, above which methanol synthesis and RWGS rates exhibit first and zeroth order dependence on P H2 , respectively. HCOOH** also accumulates on the surface with increasing P CO2 as informed by the Langmuir-type dependence on P CO2 (0.25-23 bar) for both methanol synthesis and RWGS. As both HCOOH** and H* have one H-atom per site occupied, the two species share the same P H2 dependence and give rise to CO 2 reaction orders that are independent of P H2 . Surface coverages determined based on kinetic analyses are further corroborated with DFT-derived adsorption energies that show favorable HCOOH** adsorbate-adsorbate interactions as well as repulsive interactions for bidentate formate (HCOO**) on H*-saturated surfaces. Methanol selectivity remains invariant with P CO2 and P CO despite CO inhibiting reaction rates, thereby demonstrating methanol synthesis and RWGS occur on the same active site. In contrast, water preferentially inhibits methanol synthesis rates, increases methanol synthesis H 2 reaction order from 1.0 to 1.5, and alters the methanol synthesis H 2 /D 2 kinetic isotope effect; the inhibitory effect of H 2 O thus cannot be attributed to competitive adsorption alone and instead reflects a change in the rate-determining step for methanol synthesis. The disparate kinetics of methanol synthesis and RWGS evince a branching pathway where methanol is formed from formates and CO is formed from carboxylates. The presented work thus identifies the relevant surface species, underscores the distinct catalytic role of water in branching methanol synthesis and RWGS pathways, and, in doing so, details a mechanistic picture that yields predictable rates and reaction orders for both methanol synthesis and RWGS on Cu-based CO 2 hydrogenation catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Critical statistical assessment of data in metal additive manufacturing

Obtaining high quality data reflecting the relationships between the additive manufacturing (AM) process parameters, material microstructure and mechanical properties is crucial for the use of machine learning in AM. A database of over 4,000 data entries of metal AM was created thanks to a large number of literature studies on key process parameters and indicators of build quality. Meta-analysis reveals critical biases in the literature. Firstly, majority of studies report only high quality builds, these imbalances in reporting result in weak correlation between process parameters, properties and consolidation, limiting the ability of machine learning models to generalize beyond optimized conditions. Nevertheless, the trained models accurately predict yield strength ($R^2 = 0.85$), suggesting that certain process–property relationships are effectively captured within these models. Secondly, quantitative microstructural data are largely absent, limiting the learning of the microstructure-mechanical properties relationships. Finally, current process window identification is based largely on the consolidation, despite significant uncertainty in its measurement. It is important to identify the process map on the basis of not only the consolidation, but also mechanical behaviour under loading. Such a identification shows that 316 L and Inconel have much larger process map (i.e. highly printable) in comparison to the AlSi10Mg and Ti6Al4V.

Additive manufacturing↗

DOE OSTI2024

Surface Studies of β-1,3,5,7-Tetranitro-1,3,5,7-Tetrazoctane and Pentaerythritol Tetranitrate from Density Functional Tight-Binding Calculations and Implications on Crystal Shape

Here, we use density functional tight-binding (DFTB) theory to calculate the surface energies of two energetic crystals: monoclinic β-1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (β-HMX) and tetragonal pentaerythritol tetranitrate (PETN). The results are then employed to determine crystal shapes using the Bravais–Friedel–Donnay–Harker, attachment energy, and surface energy models. We find that energy-based models yield predictions in good agreement with experimental observations. Additionally, we propose a simple model that reframes surface energy as a measure of the lost intermolecular interactions during the formation of a surface from the bulk. The model accurately captures the results from the DFTB calculations and enables us to explain and predict surface energies as a function of the local molecular environment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Machine Learning Interatomic Potentials for Modeling Framework Flexibility and Water Uptake in NbOFFIVE-1-Ni Metal–Organic Framework

Metal–organic frameworks (MOFs), with their distinctive porous structures and tunable chemical properties, have shown immense promise in the separation and storage of gases. Currently, the accurate simulation of their adsorptive properties remains challenging, especially for systems where the molecules fit very tightly into the pores. Traditional simulation methods often approximate the frameworks as rigid and do not account for the framework flexibility seen in materials such as NbOFFIVE-1-Ni. First-principles molecular dynamics (FPMD) simulations offer the desired accuracy in modeling this flexibility but are limited by their extensive computational demands, rendering them impractical for long simulations. Conversely, classical force field-based simulations offer computational efficiency but lack the necessary accuracy. Here, to break this accuracy-efficiency trade-off, we have developed machine learning interatomic potentials trained on energies and forces from FPMD to model the framework flexibility of NbOFFIVE-1-Ni in the presence of water over nanosecond time scales. Furthermore, by integrating MLIP-driven molecular dynamics (MLIP-MD) with grand canonical Monte Carlo (GCMC) simulations, we further incorporated framework flexibility into adsorption predictions, yielding water adsorption isotherms that better align with experimental data compared to those of conventional GCMC simulations. These advances offer new opportunities for the design and optimization of MOFs in gas storage and separation applications.

adsorption↗

DOE OSTI2022

Electrochemical Hydrogenation of CO on Cu(100): Insights from Accurate Multiconfigurational Wavefunction Methods

Copper (Cu) remains the most efficacious electrocatalyst for electrochemical CO 2 reduction (CO 2 R). Its activity and selectivity are highly facet-dependent. We recently examined the commonly proposed rate-limiting CO hydrogenation step on Cu(111) via embedded correlated wavefunction (ECW) theory and demonstrated that only this higher-level theory yields predictions consistent with potential-dependent experimental kinetics. Here, to understand the differing activities of Cu(111) and Cu(100) in catalyzing CO 2 R, we explore CO hydrogenation on Cu(100) using ECW theory. We predict that the preferred pathway involves the reduction of adsorbed CO (*CO) to *COH via proton-coupled electron transfer (PCET) at working potentials, although *CHO also may form with a kinetically accessible but higher barrier. In contrast, our earlier work on Cu(111) concluded that *COH and *CHO formation via PCET are equally feasible. Furthermore, this work illustrates one possible origin of the facet dependence of CO 2 R mechanisms and products on Cu electrodes and sheds light on how the selectivity of CO 2 R electrocatalysts can be controlled by the surface morphology.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2025

Low Temperature CO 2 Hydrogenation on Unsupported Mo 2 C Catalysts

CO 2 hydrogenation to methanol, a key reaction for decarbonizing the fuel and chemical industries, requires catalyst formulations that hydrogenate CO 2 selectively to methanol at temperatures where methanol conversion is not significantly equilibrium limited (<423 K). Herein we report continuous CO 2 hydrogenation at low temperatures (348-408 K, H 2 /CO 2 = 0.1-50, 5-35 bar) with high selectivity to methanol (up to ca. 80%) over unsupported β-Mo 2 C catalysts. Active site density quantification via titration with trifluoroacetic acid at reaction temperatures enables an assessment of site-specific rates. Methanation and reverse water gas shift (RWGS) occur concurrently with methanol synthesis during CO 2 hydrogenation over Mo 2 C. Reaction pathway analysis, product cofeeds, and reversibility formalisms show that all products form through primary reaction pathways from CO 2 , but secondary reactions of CO contribute significantly to rates of methanation. Dependences of forward rates on reactant and product concentration determined by independently varying the CO 2 , H 2 , CO, H 2 O, CH 3 OH, and CH 4 pressure in conjunction with reversibility formalisms reveal that all products form through H-assisted CO 2 activation and involve partially hydrogenated CO 2 -derived intermediates. Here, these inferences were verified by quantitative agreement between measured site-time yields and site-time yields predicted by closed form kinetic rate expressions in an integral reactor model over widely varying conditions (85-2000 kPa H 2 , 80-1500 kPa CO 2 , 0-45 kPa H 2 O, 0-21 kPa CO, 0-25 kPa CH 3 OH, 0-75 kPa CH 4 , 5-87 mol Mo s s mol CO 2 -1 ). Coverages calculated based on the kinetic model reveal that the Mo 2 C surface is covered with bidentate CO- and CO 2 -derived intermediates of the stoichiometry H 2 CO 2 and H 2 CO, indicating that H 2 and CO x do not compete for surface occupancy but instead adsorb cooperatively to form partially hydrogenated intermediates. Hydrogenation of the CO-derived H 2 CO** intermediate favors methanation, while hydrogenation of CO 2 -derived H 2 CO 2 ** favors methanol synthesis. Together, these findings demonstrate the ability of unsupported Mo 2 C to catalyze the hydrogenation of CO 2 to methanol at low temperatures and provide insight into the reaction network and mechanisms involved in its formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Dynamic Bottlebrush Polymer Networks: Self-Healing in Super-Soft Materials

Here, we introduce a design strategy to expand the range of accessible mechanical properties in covalent adaptable networks (CANs) using bottlebrush polymer building blocks. Well-defined bottlebrush polymers with rubbery poly(4-methylcaprolactone) side-chains were crosslinked in formulations that include a bislactone and strong Lewis acid (tin ethylhexanoate). The resulting materials exhibit tunable stress relaxation rates at elevated temperatures (160–180 °C) due to dynamic ester crosslinks that undergo transesterification with residual hydroxy groups. Varying the crosslinker loading or bottlebrush backbone degree of polymerization yields predictable low-frequency shear moduli ca. 10–100 kPa, well below values typical of linear polymer CANs (1 MPa). These extensible networks can be stretched to strains as large as 350% before failure and undergo efficient self-healing to recover >85% of their original toughness upon repeated fracture and melt processing. In summary, molecular architecture creates new opportunities to tailor the mechanical properties of CANs in ways that are otherwise difficult to achieve.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗