Search NASASearch

SEARCH · Search NASA

Results for “Polymer melt mixing”

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.

226 records · Page 11

Fingerprinting Uranium Oxides with Electron Energy Loss Spectroscopy Supported by Theoretical Computations

Uranium oxides occur in a variety of phases that differ in crystal structure and uranium oxidation states. Electron energy loss spectroscopy (EELS) is one of the few techniques that has sufficient spatial resolution and sensitivity to electronic structure to distinguish amongst phases at the nanoscale. However, beam-sensitive materials such as uranium oxides are subject to spectral modification due to interactions with the electron beam. Therefore, theory support is essential to reliably exclude the impact of beam damage and generate true reference datasets. Here we use a comparison of theoretical and experimental spectra to probe the impact of beam damage on O K-edge and U N-edge (N6,7 and N4,5) EELS spectra of various single-valent and mixed-valence uranium oxide bulk phases. Using a low-dose experimental set-up, we show that the O K-edge theoretical spectra are in excellent agreement with experiment for both peak positions and relative intensities of respective peaks. In contrast, U N-edge features are less distinguishing due to the partially localized nature of the U 5f orbitals and overlapping multiplet and spin–orbit coupling effects. This work demonstrates that O K-edge EELS is sufficiently diagnostic to distinguish a wide range of uranium oxides and that the experimental approach used here minimizes beam damage and allows valence state discrimination across the U(IV), U(V) and U(VI) series. When combined with imaging modes available in electron mi-croscopy, the work enables detailed investigation and characterization of uranium redox transformations at the nanoscale.

Carbone, Jacopo

Tailoring Ion Transport in Li 3‐3y Ho 1+y Cl 6‐x Br x via Transition‐Metal Free Structural Planes and Charge Carrier Distribution

Abstract Localized atomistic disorder in halide‐based solid electrolytes (SEs) can be leveraged to boost Li + mobility. In this study, Li + transport in structurally modified Li 3 HoCl 6 , via Br − introduction and Li + deficiency, is explored. The optimized Li 3‐3 y Ho 1+ y Cl 6‐ x Br x achieves an ionic conductivity of 3.8 mS cm −1 at 25 °C, the highest reported for holmium halide materials. 6,7 Li nuclear magnetic resonance and relaxometry investigations unveil enhanced ion dynamics with bromination, attaining a Li + motional rate neighboring 116 MHz. X‐ray diffraction analyses reveal mixed‐anion‐induced phase transitions with disproportionate octahedral expansions and distortions, creating Ho‐free planes with favorable energetics for Li + migration. Bond valence site energy analysis highlights preferred Li + transport pathways, particularly in structural planes devoid of Ho 3+ blocking effects. Molecular dynamics simulations corroborate enhanced Li + diffusion with Br − introduction into Li 3 HoCl 6 . Li‐Ho electrostatic repulsions in the (001) plane presumably drive Li + diffusion into the Ho‐free (002) layer, enabling rapid intraplanar Li + motion and exchange between the 2d and 4h sites. Li 3‐3 y Ho 1+ y Cl 6‐ x Br x also demonstrates good battery cycling stability. These findings offer valuable insights into the intricate correlations between structure and ion transport and will help guide the design of high‐performance fast ion conductors for all‐solid‐state batteries.

Chemistry

Simulation to a Newborn Supernova Remnant from a Low-mass Iron Core Star

Supernova remnant observations show a high degree of asymmetry, mixing, and inhomogeneity. These asymmetries are seeded during the early seconds of the explosion and are further enhanced and modified as the shock and ejecta move through the stellar progenitor and into the circumstellar medium. We present simulations of a 9.6 M⊙ zero-metallicity progenitor initialized after shock revival and evolved for several years when the ejecta is in the circumstellar medium. A suite of 1D and 2D simulations examines the effects of neutron-star wind and radioactive decay heating. In 1D, decay heating forms a low-density bubble that suppresses the reverse shock. While in 2D, the heating is localized to metal-rich pockets, inflating them and compressing the surrounding material into dense shells. In 3D, the neutron-star wind and decay heating modify the plume morphology, producing more large-scale structures. The extended plume morphology leads to an asymmetrical shock breakout. After breakout, the leading plumes cannot keep up with the shock front, resulting in deceleration and fragmentation by the reverse shock while retaining the large-scale asymmetry. The projected ejecta morphology and velocities are strongly viewing angle dependent. The relatively uniform metal-rich distribution does not resemble the strongly inhomogeneous ejecta structure of Cas A. The 160-isotope decay network shows that 24.4% of the radioactive heating comes from decay chains other than the canonical 56Ni chain. The low explosion energy, low 56Ni yield, and Ni/Fe ratio greater than unity suggest an observational signature similar to an electron capture supernova.

Neopane, Sudarshan [University of Tennessee (UT)]

Molecular origin of anisotropic shear elastoplasticity in chitin

Chitin nano- and mesoscale structures present in the exoskeleton of crustaceans exhibit exceptional longitudinal stiffness and toughness, rivaling or even exceeding that of many synthetic polymer architectures. Here, we reveal the origin of the asymmetric shear response in chitin multiscale architectures, marked by pronounced anisotropy in deformation. Under shear aligned with the molecular axis, chitin accommodates strain through coherent atomic rearrangements that enable elastic recovery. In contrast, shear applied perpendicular to this axis induces liquid-like plasticity via localized sliding. These results demonstrate the intrinsic mechanical anisotropy of chitin, underpinning a dual function in resisting repetitive loading while dissipating internal stress during high-strain events. Our findings establish the molecular basis of shear elastoplasticity in multiscale chitin structures, wherein axial elasticity supports energy storage in load-bearing regions, whereas transverse plasticity enables controlled energy dissipation. These atomic-scale insights lay a foundation for the predictive design of chitin-based materials with tunable strength-toughness profiles.

Wan, Zhangmin [University of British Columbia, Van

Hybrid Data‐Driven Discovery of High‐Performance Silver Selenide‐Based Thermoelectric Composites

Optimizing material compositions often enhances thermoelectric performances. However, the large selection of possible base elements and dopants results in a vast composition design space that is too large to systematically search using solely domain knowledge. To address this challenge, a hybrid data-driven strategy that integrates Bayesian optimization (BO) and Gaussian process regression (GPR) is proposed to optimize the composition of five elements (Ag, Se, S, Cu, and Te) in AgSe-based thermoelectric materials. Data is collected from the literature to provide prior knowledge for the initial GPR model, which is updated by actively collected experimental data during the iteration between BO and experiments. Within seven iterations, the optimized AgSe-based materials prepared using a simple high-throughput ink mixing and blade coating method deliver a high power factor of 2100 µW m −1 K −2 , which is a 75% improvement from the baseline composite (nominal composition of Ag 2 Se 1 ). In conclusion, the success of this study provides opportunities to generalize the demonstrated active machine learning technique to accelerate the development and optimization of a wide range of material systems with reduced experimental trials.

36 MATERIALS SCIENCE

Process control-enabled mitigation of microstructural and plastic heterogeneities in additively manufactured Grade 91 steel

Synergizing wire arc-directed energy deposition (WA-DED) additive manufacturing (AM) with particle-strengthened creep strength-enhanced ferritic (CSEF) steels enables fabrication and repair of critical power-plant components. Investigations focused on fusion-welded particle-strengthened CSEF steels, such as Grade 91 steel, have linked microstructurally heterogeneous regions—forming due to heat affected zones (HAZ)—with premature failure during elevated temperature service. Fusion-based AM, including WA-DED, likewise generates microstructurally and plastically heterogeneous regions due to spatiotemporally varying thermokinetics during deposition. However, works investigating such microstructural heterogeneities, their implications for mechanical behavior, and strategies to mitigate their formation remain scarce. This work identifies microstructurally and plastically heterogeneous regions within the WA-DED-processed Grade 91 steel. Spatial microhardness variations in the as-fabricated specimen correlate with the variation in the attributes of grain, martensitic microstructure, and precipitates across the fusion zone and HAZ. Digital image correlation-enabled tensile tests performed at 500 °C revealed pronounced deformation localization and a wave-like strain distribution, with wavelength close to the melt pool depth, indicating susceptibility of the as-fabricated components to premature creep failure. Such heterogeneity in microstructural and mechanical behavior was attributed to recurring solid-state phase transformations. Subsequently, an interlayer temperature control strategy was implemented, wherein maintaining interlayer temperature above the martensitic start temperature mitigated the heterogeneous microstructural and plastic response in the as-fabricated condition. Findings open pathways to achieving deformation-localization- and creep-resistant microstructures in WA-DED fabricated particle-strengthened CSEF steel components, reducing reliance on post-welding heat treatments—conventionally required to enhance creep resistance—and enabling on-demand, short lead-time fabrication of next-generation power-plant components.

Heat affected zones

Chemo‐Mechanical Coupling in Hydrogels: Dynamics in the Diffusion‐Limited Regime

Hydrogels are characterized by substantial volume changes in response to external stimuli, making them promising candidates for developing smart materials with enhanced adaptability and responsiveness. By integrating chemical reactions, hydrogels acquire dynamic and tunable responsiveness to external stimuli through chemo‐mechanical coupling, expanding their potential in emerging applications. However, capturing their transient behavior remains challenging due to the complex interplay of chemical reactions, solvent transport, and polymer network deformation. Classical theories capture equilibrium swelling but fail to describe time‐dependent phenomena. To address this, a time‐dependent continuum model is developed that explicitly couples these processes. Volume phase transition in hydrogels with homogeneous chemical reactions is investigated, then the effects of reaction kinetics on these transitions are analyzed. The coupling of these mechanisms is further explored through a study of transient mechanical instabilities. To illustrate the impact of distinct reaction and diffusion timescales, a photo‐active gripper is studied for robotic applications. Finally, a photo‐active microswimmer that exhibits non‐reciprocal motion is proposed to highlight the solvent diffusion for locomotion at the micro‐ and nano‐ scales. The work establishes that transient dynamics of chemo‐mechanical hydrogels generate functions not accessible by steady state models and provides a predictive platform for designing adaptive materials in emerging applications.

chemo-mechanical coupling

Extreme confinement of hydrogen gas within fullerenelike nanoporous carbon

Nanoporous carbons and carbon nanostructures can store hydrogen at cryogenic temperatures but lack the volumetric and gravimetric capacity to be industrially significant. Recent inelastic neutron scattering experiments suggest a highly dense phase of hydrogen at temperatures well above the melting point of solid hydrogen. However, it remains unclear how pore geometry and intermolecular interactions enable these dense phases to exist, with dispersion (van der Waals) or electrostatic/induction suggested to be the key effects in slit and curved pores but their relative contributions have yet to be quantified. In this paper, we perform benchmark electronic structure calculations allowing the interactions between planar and curved aromatic molecules with hydrogen to be accurately determined. Dispersion was found to dominate over electrostatic and inductive effects with some many-body charge transfer (Dobson type-A) effects needed to capture the most highly curved structures. Density functional methods that include type-A many-body effects were found to accurately describe the intermolecular interactions at a fraction of the cost of coupled-cluster simulations and these approaches were used to calculate the energies inside large carbon bowl and slit pores. The interaction energies inside the bowl pores were found to depend on the orientation of the hydrogen molecule. This rotational barrier, modeled as a quantum hindered rotor, could reproduce the peak splitting observed in inelastic neutron scattering experiments, with weak splitting arising from bowl-like fullerene pores and strong splitting from highly confining nanotubelike pores. Increasing the fraction of such curved pores in nanoporous carbons may therefore offer a pathway to enhance their hydrogen-storage capacity. Moreover, the preferential adsorption of ortho hydrogen on nanotubelike pores could enable the storage of high-density hydrogen without the need to remove heat produced during the ortho-para hydrogen conversion.

36 MATERIALS SCIENCE

Joint scheduling of energy, fast and primary frequency response reserves in integrated transmission–distribution networks

Inverter-based distributed energy resources (DERs) connected to distribution networks (DNs) can provide fast frequency support, but their reserve deliverability depends on feeder constraints and differs from synchronous primary frequency response (PFR). Existing transmission–distribution coordination studies usually treat reserve generically or neglect feeder-level feasibility, while frequency-security scheduling studies rarely represent distribution feeders explicitly. This paper develops a bi-level day-ahead scheduling framework for integrated transmission–distribution networks that jointly clears energy, transmission-side PFR, and distribution-side fast frequency response (FFR) under exogenous hourly inertia and largest-loss inputs from an external unit commitment (UC) schedule. The transmission problem is modeled with DC-optimal power flow (OPF) and closed-form second-order cone (SOC) frequency-security constraints, whereas each DN is represented by a reserve-aware branch-flow AC-OPF so that scheduled fast reserves remain deliverable during activation. The bi-level problem is reformulated through Karush–Kuhn–Tucker (KKT) conditions into a mixed-integer SOC program, and a penalty term is used to tighten the distribution-network relaxation. In the reduced test system, lower exogenous inertia increased the required primary response from 179.64 MW to 191.08 MW, distribution-side fast response reduced total frequency-response procurement by up to 4.9%, and neglecting distribution constraints overstated the combined distribution-side energy and reserve award by up to 18%. In the expanded study, the largest case was solved in 2.02 s with a 0.00% optimality gap. Time-domain simulations kept the frequency nadir above 59.0 Hz in all tested hours. These results demonstrate the value of fast-response modeling and distribution-feasible reserve delivery in coordinated market clearing.

Noh, Seung-Gil

Fracture Analysis of Cohesive Zone Models for Modeling Residual Stress Induced Delamination in Composite Structures

A fracture study of coupon-scale composite cylinders with embedded defects was conducted with an objective to assess and validate a modeling approach using two available cohesive material models. The study included experimental and simulation evaluations of initiation of crack growth and progression. Interrupted thermal experiments used acoustic emissions monitoring to identify the onset of crack progression during each cooling interval and ultrasonic scanning provided images of defect growth. Verification, validation, and uncertainty quantification (VVUQ) processes were performed in the assessment of the simulation predicted temperature at which crack propagation begins (quantity of interest). The Sobol sensitivity analysis identified the hoop direction elastic modulus in the carbon fiber reinforced polymer (CFRP) plies as the most influential parameter for simulations using both cohesive models, accounting for at least 70% of the variation in the temperature at crack propagation. The UQ temperature range for the Tvergaard-Hutchinson model was higher (more conservative) than the experimental acoustic measurement indicators of crack progression, while the temperature range for the Thouless-Parmigiani model enveloped the experimental data points for the primary defect size of 0.75 x 1 in. The simulations could not capture the stable crack growth indicated in the experiments. This is likely due to the models’ inability to represent anisotropic fracture toughness attributed to the structure of the orthotropic fiber weave in a woven composite laminate.

42 ENGINEERING

Flame stabilization in DME spray flames under engine-relevant conditions characterized by OH* chemiluminescence and formaldehyde laser-induced fluorescence

The transient and quasi-steady flame structures of Dimethyl Ether (DME) fuel sprays, produced by a single-hole injector (Spray D), were investigated using Planar Laser-Induced Fluorescence (PLIF) and chemiluminescence imaging in a constant-volume chamber under Engine Combustion Network (ECN) Spray A conditions (900 K ambient temperature, 60 bar ambient pressure, 1500 bar injection pressure, and 22.8 kg/m 3 ambient density). Low-temperature chemical reaction zones were visualized using formaldehyde (CH 2 O) PLIF with 355 nm excitation, while high-temperature flame regions were captured via chemiluminescence imaging of excited-state hydroxyl radicals (OH*). Both transient and quasi-steady flame structures clearly show the transition from CH 2 O to OH*, highlighting the progression from low- to high-temperature combustion, while the position of the flame is displaced for DME compared to reference hydrocarbon n-dodecane. Homogeneous reactor calculations with detailed chemistry and using adiabatic mixing for initial temperature show that CH 2 O peaks are significantly higher for DME at the same equivalence ratio, with a higher heat-release during the cool-flame regime with respect to the fuel heating value. Thus, the cool-flame dynamic as a precursor to high-temperature combustion and flame stabilization exhibit distinct behavior for DME relative to conventional hydrocarbons, and these phenomena are effectively resolved through the soot-free nature of DME and the high-speed, time-resolved diagnostics.

CH2O laser-induced fluorescence

Conserved Macromolecular Architecture of Poplar Secondary Cell Walls Revealed by ssNMR and Atomistic Modeling

The macromolecular architecture of plant secondary cell walls governs wood's mechanical and biochemical properties, yet its natural intra-species variability remains poorly characterized. Here, we combined 13C solid-state NMR (ssNMR), multivariate statistical analysis, and molecular modeling to profile nanoscale structure across 13 genetically diverse Populus trichocarpa genotypes grown in 13C-enriched atmospheres. SsNMR-derived phenotypes spanning composition, structure, mobility, and inter-polymer proximities reveal a conserved architecture, with a subtle yet coordinated variation organizing into dominant structural and secondary mobility axes. A representative atomistic model captures these features and reproduces experimental metrics. Molecular dynamics simulations support a weak but consistent positive correlation between cellulose abundance and crystalline-like order, with interior cellulose chains enriched in tg (trans-gauche) conformations without expanding crystalline cores. Together, experiment and simulation reveal a genetically buffered, broadly conserved nanoscale architecture across genotypes, where subtle fine-tuning of cellulose bundling and matrix packing balances mechanical performance with biological function.

09 BIOMASS FUELS

Dimensional Evolution Guides Property Control in the A n Cu 4– n TiS 4 Semiconductor Series

Through progressive reduction of the three-dimensional (3D) covalent network of Cu 4 TiS 4 , we isolate seven new members of the A n Cu 4–n TiS 4 family (A = alkali metal; n = 0–4), spanning 3D, 2D, 1D, and 0D structural fragments. The dimensional reduction is rational, as it preserves the edge-sharing connectivity between [CuS 4 ] 7– and [TiS 4 ] 4– tetrahedra across the series. This structural evolution is driven by the stepwise substitution of Cu with alkali metals, guiding the formation of fragments with reduced dimensionality. The effects of “n” and “A” on the crystal structures, stabilities, electronic structures, and optoelectronic properties are profound, demonstrating that the manipulation of alkali metal size and A n Cu 4–n TiS 4 stoichiometry enables predictable variations in structure and properties. For example, the n = 0 and n = 4 end members of the A n Cu 4–n TiS 4 family set the range of achievable band gaps with 2.00 eV for Cu 4 TiS 4 , 2.60 eV for Na 4 TiS 4 , and intermediate values for the n = 1–3 members. Notably, CsCu 3 TiS 4 exhibits exceptional air stability and congruent melting, with density functional theory (DFT) calculating moderate hole and electron effective masses in specific crystallographic directions (mh = 1.24m 0 , me = 0.87m 0 ). Additionally, A 3 CuTiS 4 (A = Na, K, Rb) displays direct band gap behavior and long photoluminescence lifetimes of 2.3–8.6 μs, and K 3 CuTiS 4 has a PLQY of 5.19%. These findings underscore the potential of the A n Cu 4–n TiS 4 family for applications in optoelectronics and demonstrate widely applicable design concepts that unveil rational stoichiometries within a given composition space to generate a series of crystal structures related through an evolving covalent dimensionality that corresponds to a predictable electronic structure and property progression.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Vanishing dynamic strength measured in a transient high-pressure phase of tin

This work reports two independent experimental estimations of the dynamic strength of the body-centered tetragonal (bct) 𝛾 phase of tin, which occurs at pressures above about 9 GPa and transforms promptly back to the ambient 𝛽 phase upon pressure release. Measuring strength in such a transient high-pressure phase is challenging. One measurement used free-surface Richtmyer-Meshkov instabilities generated with gas gun impact. Another set of measurements used ramp-release loading in a pulsed-power facility. Strength estimations came from comparison to simulations using a comprehensive multiphase modeling framework that includes equation of state, shear moduli, and strength separately for each phase, and treats mixed-phase regions. Both experimental methods found the 𝛾-phase strength to be very low, within experimental uncertainty of zero. The existing literature on other materials, by contrast, almost universally reports higher strength in transient high-pressure phases compared with ambient phases. Recent Molecular-Dynamics simulations on tin in the literature showed almost zero deviatoric strength during a deformation-induced bct → bct transformation in which one of the 𝛾-phase 𝑎 axes “flips” to the 𝑐 axis. This reorientation currently provides the most plausible explanation for the low observed strength in 𝛾-phase tin.

36 MATERIALS SCIENCE

Revealing the Roles of CuF 2 /NiF 2 Incorporation in the Electrochemical Performance of FeF 3 Cathodes in Solid‐State Batteries

Mixed metal fluorides have been considered as a promising candidate to lower the voltage hysteresis of conversion-type iron fluoride cathodes, but their cycling stability is limited due to transition metal dissolution and interphase growth in liquid electrolyte batteries. Here, we study the role of incorporating CuF 2 and NiF 2 in the electrochemical performance of FeF 3 cathode in halide-based solid-state batteries to test whether we can transfer the kinetic benefit of low voltage hysteresis to solid-state batteries while using solid electrolyte to eliminate transition metal dissolution and stabilize the interphase. Synchrotron X-ray absorption spectroscopy results indicated the redox reactions are attributed to Cu 0 /Cu + and Fe 0 /Fe 2+ in 25CuF 2 -75FeF 3 and Ni 0 /Ni 2+ and Fe 0 /Fe 3+ in 10NiF 2 -90FeF 3 . While no apparent improvement in electrode kinetics can be observed, the incorporation of CuF 2 and NiF 2 can largely improve the cycling stability of FeF 3 cathodes. In conclusion, the results demonstrate the advantages of using solid-state concept to improve the cycling stability of conversion-type cathodes.

25 ENERGY STORAGE

Development of an Amine Oxide Polyzwitterion Brush Martini Model with Polarizable Water and Ions

Abstract In order to accurately simulate the fouling process of proteins onto polyzwitterion brushes, models that accurately capture the hydration properties and chain conformations of such brushes must first be established. We developed a Martini coarse-grained (CG) model for amine oxide polyzwitterion (PNOMA) brushes, a promising class of antifouling materials, in polarizable water and ions by fitting to all-atom bond and angle distributions, monomer hydration free energy, monomer–monomer distance potential of mean force (PMF), and monomer–salt radial distribution functions (RDFs). Martini 2.2P was selected for compatibility with the established polarizable water and ion models. For comparison with PNOMA, we also constructed models for conventional sulfobetaine (PSBMA) and phosphorylcholine (PMPC) polyzwitterions and the polycation PMETAC using established nonbonded bead types from the literature and refitting bond and angle potentials. We simulated each polymer brush chemistry for varying grafting density and chain length, validating brush height scaling relations against experimental data. The CG models captured the relative hydration strengths among different polyzwitterion chemistries, and brush heights extrapolated to higher molecular weights are in agreement with experimental ellipsometry data. We find that chain swelling of the superhydrophilic PNOMA brushes lies between that of the traditional polyzwitterions PSBMA/PMPC and the polycation PMETAC. For PNOMA brushes in NaCl solution, simulated brush height decreases with salt concentration due to the selectively strong interactions between amine oxide and sodium ions.

Walker, Christopher C. [Oak Ridge National Laborat

Design Trade-Offs in Composite Fuel Cell Membranes: Effects of Reinforcement and Chemical Additives

Perfluorosulfonic acid (PFSA) membranes are critical components in proton exchange membrane fuel cells, where performance depends on balancing ionic conductivity, mechanical durability, and chemical stability. This study characterizes a composite membrane (NC700) featuring PFSA-impregnated expanded polytetrafluoroethylene (ePTFE) reinforcement and cerium-based radical scavengers, benchmarked against unreinforced NR211. Complementary techniques, including electron microscopy, X-ray scattering, infrared spectroscopy, thermogravimetric analysis, and dynamic mechanical analysis, identify the structural and compositional strategies employed in NC700. Water sorption isotherms reveal lower water uptake for NC700 across all conditions, attributed to reinforcement and cerium incorporation. Reinforcement reduces in-plane swelling from 11% to 2.1% at 90% RH, confirming strong swelling anisotropy, while maintaining mechanical properties at elevated temperatures. While the ionic conductivity of NC700 is approximately 10% lower than that of NR211, the reduced thickness yields a 40% decrease in calculated area-specific resistance, suggesting the composite architecture can favorably shift the conductivity-stability trade-off. The composite structure also reduces gas permeability, indicating potential for improved separator function alongside favorable transport properties. Systematic deconvolution of reinforcement and additive contributions shows that conductivity losses from cerium incorporation are largely offset by gains from the lower equivalent-weight polymer, providing quantitative relationships that may guide composite membrane design for fuel cells and other electrochemical applications.

25 ENERGY STORAGE

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

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

36 MATERIALS SCIENCE