Search NASA⌕ Search

SEARCH · Search NASA

Results for “High potential cathodes”

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.

857 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Humins-Derived Hard Carbon as a Low-Cost Material for Sodium-Ion Battery Anodes

The growing demand for sodium-ion batteries (SIBs) in grid storage underscores the need for electrode materials that balance performance and cost, including sustainable and robust carbon sources. The equitable and abundant distribution of materials for SIBs, along with their superior low-temperature performance, safety, and fast-charging capability, further distinguishes them from lithium-ion batteries (LIBs). Hard carbon (HC) is the state-of-the-art anode for SIBs, but current commercial HC production is localized mostly to one region of the world, raising concerns about supply chain vulnerability, critical material dependency, and environmental aspects. Here, the first demonstration of humins, an abundant biorefinery byproduct, as a precursor for HC anodes for sodium-ion storage is reported. Humins were carbonized at 1100 degrees C-1300 degrees C, and the resulting materials were subjected to comprehensive materials and electrochemical characterization. Among the temperatures studied, humins-derived hard carbon synthesized at 1200 degrees C delivers an initial reversible capacity 270mA h g-1 , with stable cycling performance up to 500 cycles and excellent rate capability, representing the optimal performance. This study establishes humins as a promising and low-cost carbon source that provides a route to mitigate supply chain risks and valorizes an underutilized biorefinery waste stream for high-performance SIB anodes.

25 ENERGY STORAGE↗
DOE OSTI2025

Symmetry‐Breaking Strategy Yields Dopant‐Free Small Molecule Hole Transport Materials for Inorganic Perovskite Solar Cells with 20.58% Efficiency and Outstanding Stability

Abstract Inorganic perovskites are known for their excellent photothermal stability; however, the photothermal stability of all‐inorganic n‐i‐p perovskite solar cells (PSCs) is compromised due to ion diffusion and free radical‐induced degradation caused by the use of doped spiro‐OMeTAD hole transport materials (HTMs). In this study, two isomeric donor–acceptor–donor (D–A–D) type small molecules, namely HBT and HiBT, were developed and used as dopant‐free HTMs, using 2,1,3‐benzothiadiazole or benzo[d][1,2,3]thiadiazole as acceptor moieties. The HiBT molecule, with its symmetry‐breaking features, exhibits a large dipole moment, enhanced coordination‐active sites, and a well‐aligned energy level structure, all of which contribute to passivating perovskite surface defects and improving free charge separation. As a result, inorganic CsPbI 3 PSCs with HiBT HTM achieved an impressive power conversion efficiency (PCE) of 20.58%, the highest reported for dopant‐free HTM‐based inorganic PSCs. Moreover, the enhanced hydrophobic properties of HiBT molecules, coupled with their ability to passivate perovskite surface defects, contribute to significantly improved device stability. The unencapsulated devices based on HiBT HTM retained over 83% and 80% of their initial efficiency after being stored at 85 °C for 50 days and undergoing maximum power point (MPP) tracking at 85 °C for 1100 h, respectively. These results highlight that the symmetry‐breaking strategy is an exceptionally effective approach for designing efficient, dopant‐free small molecule HTMs, significantly contributing to both the high efficiency and enhanced stability of all‐inorganic PSCs.

Chemistry↗
DOE OSTI2026

Precise Modeling of a Complex Solenoidal Magnetic Field Using a Combination of Analytic Functions and a PINN

We demonstrate an iterative approach to modeling a sparsely measured magnetic field in a large-bore solenoid. This approach uses a hybrid of traditional and machine learning techniques. The traditional technique is a linear least-squares fit using a series solution to Laplace's equation, while the machine learning technique involves the training of a physics-informed neural network (PINN) on the least-squares fit residuals. We use a newly defined activation function "DELTAsnake," a modification to the snake activation function proposed by Ziyin et al. that allows for stronger curvature and non-monotonicity. The combined model approximately obeys Maxwell's equations to a level sufficient for producing high quality physics simulations and analysis. Our approach is applied to a highly realistic calculation of the expected magnetic field in the Mu2e experiment's Detector Solenoid which includes a simple model for the expected statistical measurement uncertainties. Using ten toy measurement simulations, we demonstrate the capabilities of our model in comparison to the least-squares method alone; the least-squares method alone results in a reduced chi-squared statistic of ${2.15 \pm 0.01}$, while our approach improves the reduced chi-square to ${1.034 \pm 0.005}$. Furthermore, for an average toy simulation, we show that the range of the RMS of the three field component residuals reduces from ${0.07-0.37}$ Gauss to ${0.05-0.07}$ Gauss. We find that this novel method is robust against a realistic systematic uncertainty deriving from Hall probe calibration bias and can be used to significantly reduce the number of measurements required to achieve an accurate model.

Kampa, Cole [Caltech] (ORCID:0000000192972920)↗
DOE OSTI2026

Conditional diffusion machine-learning framework for mapping valence electron distribution from convergent beam electron diffraction

Quantitative convergent beam electron diffraction (CBED) enables determination of aspherical valence electron distributions through refinement of low-order structure factors, which are highly sensitive to chemical bonding and charge density variations. However, conventional quantitative CBED (QCBED) requires solving a highly nonlinear inverse problem with many coupled parameters, and computationally intensive dynamical diffraction calculations, making it time-consuming and difficult to apply to complex systems. More broadly, reconstructing charge density and orbital electron distribution from diffraction data has long been a central challenge in both x-ray and electron crystallography. Here, in this study, we introduce an artificial-intelligence (AI)-based framework that replaces traditional refinement with a data-driven inverse solver. Using a large synthetic CBED dataset generated by Bloch-wave simulations, we train a conditional diffusion model to directly infer crystal structural parameters and multipole density formalism parameters, and hence valence electron distributions, from CBED patterns alone. By learning from forward simulations across realistic parameter space, the model effectively solves the inverse problem. Compared with direct regression approaches, the diffusion-based framework provides posterior parameter distributions for rigorous uncertainty quantification while preserving quantitative fidelity and reducing analysis time by orders of magnitude. By eliminating the need for external single-crystal x-ray diffraction data and complex nonlinear refinement, this approach enables practical, high-throughput, and in situ quantitative CBED, enabling real-time mapping of valence electron distributions and their correlation with functional responses in quantum and energy materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Modeling a generic TRISO-fueled heat pipe microreactor using SCALE: Depletion, transportation criticality, and shielding

This paper demonstrates the applicability of the SCALE code system to tristructural-isotropic (TRISO)-fueled heat pipe microreactors through depletion, transportation criticality, and shielding analyses of a generic reference design. The study conducted supports US Nuclear Regulatory Commission code readiness efforts for advanced non–light-water reactor technologies and is intended as a code capability demonstration rather than as an optimization of a specific microreactor design. The modeled reactor employs high-assay low-enriched uranium (HALEU) uranium oxycarbide (UCO) TRISO fuel and beryllium oxide (BeO) reflectors and operates at 7.5 MWth with a nominal lifetime of about 3 effective full power years. Representative cases for fresh and irradiated cores were selected to exercise SCALE methods relevant to reactor operation and post-irradiation transport. The discharged-core decay heat is approximately 6% of operating power immediately after shutdown. Transportation criticality calculations show that internal water ingress is the dominant reactivity effect, with fully flooded fresh core and irradiated core configurations remain above the subcriticality criterion, even with the available control mechanisms. Shielding calculations for a simplified transportation package indicate that normal-condition dose rates are governed mainly by shielding thickness and cooling time, whereas the breached hypothetical accident case is governed primarily by cooling time. Overall, the study shows that SCALE supports depletion, transportation criticality, and shielding evaluations efficiently for TRISO-fueled heat pipe microreactors within a single code system.

Criticality↗
DOE OSTI2026

Atomic faulting drives exceptional toughness in low thermal expansion chromium alloys

Endowing functional properties with mechanical responses in traditional metals has been a frontier topic, akin to transforming base metal into gold. Chromium and its alloys, with their functional deficiencies and limited ductility, serve as typical examples. Herein, we report a Cr96Fe4Ge1.3B1 alloy that unifies low thermal expansion (LTE, αl = 1.79 × 10-6 K-1, 200 − 315 K) with exceptional toughness (240.2 J·cm-3). The enhancement in mechanical responses is primarily attributed to layered Cr2B intermetallic precipitates, which ameliorate interfacial cohesion and simultaneously refine the grain structure. The weakened interlayer interactions within the Cr-B layers facilitate the nucleation and movement of numerous tiny stacking faults in precipitates, efficiently alleviating strain energy and resulting in marked work-hardening ability. Additionally, antiferromagnetic fluctuations in the BCC matrix contribute to the unique LTE behavior. This paves the way for the design of high-performance alloys featuring layered-symmetry precipitates.

Yu, Chengyi [University of Science and Technology ↗
DOE OSTI2026

NOx mediated formation of ε – UO 3 during thermal decomposition of uranyl nitrate hexahydrate under static air conditions

The thermal decomposition of uranyl nitrate hexahydrate in air at temperatures of 650 – 800 °C is expected to yield α – U 3 O 8 through intermediate UO 3 phases. Here, n this study, Raman spectroscopy complimented by powder X-ray diffraction revealed the unexpected formation of ε – UO 3 as an accompanying phase during thermal decomposition of UNH at 700 °C under static air conditions. Experiments demonstrated that ε – UO 3 does not form during the initial ramp up stage nor during the high temperature heating period; but instead forms during the cooling stage between 250 – 400 °C. Additional experiments revealed that UNH initially decomposes through an amorphous UO 3 intermediate prior to α – U 3 O 8 formation. Attempts to bypass the U 3 O 8 precursor in the formation of ε – UO 3 were unsuccessful, supporting the necessity of U 3 O 8 (α – U 3 O 8 in this study) in the formation of ε – UO 3 . The observed phase evolution is proposed to result from the NOx species generated during UNH thermal decomposition, which create localized oxidizing conditions within the furnace under static air conditions. Furthermore, a predominantly phase pure ε – UO 3 was synthesized under static air conditions through the addition of an extra plateau at 250 °C during the cooling step. These findings demonstrate the importance of gas-phase chemistry and cooling conditions in uranium oxide phase evolution and provide additional insight into ε – UO 3 formation pathways.

Epsilon uranium trioxide↗
DOE OSTI2026

SAM Finite Volume Method Development Status Update: GCR Application, Restart, and MultiApp

The System Analysis Module (SAM) is being developed as a modern system analysis code for advanced non-light-water-reactor safety analysis under the U.S. DOE NEAMS program. Previous feasibility studies have demonstrated that a staggered-grid finite volume method (SG-FVM), implemented under the MOOSE framework, can deliver more than an order of magnitude speedup over the existing continuous Galerkin finite element method (CG-FEM) solver for liquid-cooled, incompressible but thermally expandable flow systems. This work extends the previous effort to compressible, gas-cooled reactor applications, where pressure couples directly into the mass equation adding additional nonlinearity into the equation system. New code capabilities are implemented for pebble bed high-temperature gas-cooled reactor (PB-HTGR) analysis, including a pebble bed CoreChannel component, built-in pebble bed effective thermal conductivity model and channel-to-channel crossflow model. The capabilities are tested, benchmarked, and demonstrated for problems with increased level of model and physical complexities, including the HTTU effective thermal conductivity test, the SANA passive cooling test, and a demonstration case using the GPBR200 reactor design covering steady-state operation, DLOFC and PLOFC transients. Across all cases, the SG-FVM solver demonstrated strong robustness and efficiency, and the solutions agree well with reference results and data. The finding of this work proves that SG-FVM is a viable and efficient solver pathway for compressible, gas-cooled reactor system analysis in SAM. In addition, work has been done to successfully support SAM-FVM recover/restart code feature that is essential to reactor safety analysis applications, and MultiApp code feature that is essential to multi-scale and multi-physics simulations. In summary, this work continued from previous feasibility studies, and further demonstrated that the SG-FVM will serve as a strong foundation for SAM’s advanced solver algorithm for future deployment.

Zou, Ling↗
DOE OSTI2026

Pinning ångström-size solid ionic channels for rare-earth element separation

High-purity rare-earth elements are essential for modern technologies, yet current solvent extraction processes are energy-intensive and environmentally harmful because of inadequate selectivity and ligand toxicity. Although combining size exclusion and binding affinity can improve lanthanide separation, the role of long-range confinement remains underexplored. Here we report lanthanide separation in aqueous systems using extremely confined manganese oxide solid ionic channels with optimized layer spacing. Different lanthanides induce distinct solid-state phase transformations in manganese oxide, creating a strong driving force for separation. Two lanthanide groups, differing by ~1.4 Å in spacing, were identified and confirmed to be stable by density functional theory. The narrower confinement of heavier Group II lanthanides improves cross-group separation by increasing the dehydration barrier for lighter Group I lanthanides without inducing strong binding. Here, we further developed a strategy to pin the confinement dimensions and enhance same-group separation, increasing enrichment factors for La–Nd and La–Pr pairs from 1.6 ± 0.1 and 1.5 ± 0.1 to 5.4 ± 0.1 and 4.2 ± 0.1, respectively.

Chemical engineering↗
DOE OSTI2025

Analysis of Liquids Using the Submerged Plasma for Isotopic Detection and Elemental Resolution (SPIDER)

We developed a sensor called the Submerged Plasma for Isotopic Detection and Elemental Resolution (SPIDER) probe, which uses an atmospheric pressure glow discharge below the surface of liquids to excite species in the liquid. Through emission spectroscopy of molten salts, liquid metals, and heavy water, we demonstrated the SPIDER probe’s high resolution, accuracy, and versatility. We successfully identified trace concentrations of transition and rare-earth metals in molten salts and detected the isotopic shift of the H β → D β emission line. Our analysis revealed unconventional spectral alkali line shapes, indicating two competing excitation modes: film explosion and droplet vaporization. The film explosion mode, characterized by dense plasma, exhibited self-reversal and broadband continuum emission, while the droplet vaporization mode, associated with diffusive plasma, produced narrow-line emissions. Furthermore, by analyzing circuit transients alongside individual plasma events, we observed that film explosions generate higher currents, likely due to a shorter plasma length as the current preferentially flows through the thin liquid layer. Altogether, our results highlight the SPIDER probe’s efficacy and flexibility, making it well-suited for online material quantification of liquids in extreme environments.

AES↗
DOE OSTI2026

Terahertz‐Nanoscale Visualization of the Microscopic Spin‐Charge Architecture of Colossal Magnetoresistive Switching

Resolving sub-10 nm spin switching and the associated terahertz (THz) electrodynamics during the colossal magnetoresistance (CMR) transition is a definitive frontier in reaching the fundamental spatial, temporal, and energy-dissipation limits of spin-electronics. Yet, simultaneous control of high magnetic field, cryogenic environment, and nanometer resolution has remained an elusive benchmark for THz nanoscopy, leaving the local THz dynamics of these transitions largely unexplored. Here, we overcome these limitations by utilizing a custom-built cryogenic magneto-THz scattering-type scanning near-field optical microscopy (cm-THz-sSNOM) to resolve the near-field THz spectroscopic evolution of the magnetic field-driven CMR transition in a manganite single crystal. Our measurements provide a nanoscale visualization of the THz conductivity, capturing the moment that magnetic-field-induced spin switching triggers the transition from an antiferromagnetic insulator to a ferromagnetic metal. An ellipsoidal near-field model reveals a multi-scale transition initiated by 1–2 nm isolated spin-flip sites at low magnetic fields, which coalesce into ∼15 nm conducting regions as the threshold field is approached. These results provide an nano-THz view of CMR switching, establishing an analysis framework for mapping spin–charge–lattice–orbit–coupled dynamics at spatial scales that transcend the nominal sSNOM resolution.

Haeuser, Samuel [Iowa State University, Ames, IA (↗
DOE OSTI2026

Tokamak Energy’s pre-concept design for a fusion power plant: an overview of ST-E1

Climate change and rapidly rising energy demand, driven in part by artificial intelligence and data-centre growth, create an urgent need for stable, low-carbon, and abundant power. Fusion is a promising long-term solution, yet its commercialisation faces a fundamental paradox in today’s investment environment: pilot plants are essential to de-risk physics, engineering, and operations, but their limited lifetime energy output and high upfront costs make them difficult to finance. This paper presents Tokamak Energy’s response: ST-E1, a pre-concept design for a low-aspect-ratio tokamak power plant engineered specifically to overcome this challenge. ST-E1 is designed from the outset for phased operation—pilot and commercial phases, with an upgrade phase in between—with emphasis on commercial viability, maintainability, nuclear engineering, modularity, and upgradability. A key design principle is the deliberate separation of long-lived assets, such as the magnet cage and vacuum vessel, from replaceable in-vessel systems. This provides an attractive and credible investment approach to generate operational data and de-risk key technologies while preserving most capital-intensive assets for later commercial phases. The architecture supports continuous optimisation toward high net electric power (targeting 800–1000 MW net electric), a normalised capital expenditure of $\$$ 12–14k/kW of net electric power, and high availability (targeting > 80%). A tokamak core with a 5 m major radius, aspect ratio of 2.3, and on-plasma axis toroidal field of 5.25 T was selected to meet these objectives. This paper summarises the ST-E1 design philosophy, principal features, and development methodology. It introduces a Focus Collection of 11 papers detailing the pre-concept design of the entire tokamak and corresponding plant.

ST-E1↗
DOE OSTI2026

Spin‐Polarized Interfaces Redirect CO 2 Reduction From CO to Formate

Achieving high product selectivity in electrocatalytic carbon dioxide reduction (CO 2 RR) remains a critical challenge due to competition between multiple proton-coupled electron-transfer pathways on catalyst surfaces. Meanwhile, chirality-induced spin selectivity (CISS), which enables spin-polarized electron transport through chiral interfaces, has recently emerged as a promising strategy to modulate interfacial electrochemical reactions. Although the CISS effect has been shown to enhance selectivity and efficiency in the spin-sensitive oxygen evolution reaction (OER), its role in regulating CO 2 RR pathways and in stabilizing intermediates remains largely unexplored. Here, chiral molecules (R- and S-1,1′-bi-2-naphthyl-2,2′-diyl hydrogen phosphate, BNP) were integrated with SnO 2 to construct chiral-modified catalysts (R-BNP/SnO 2 and S-BNP/SnO 2 ). Compared with bare SnO 2 and racemic BNP-modified SnO 2 (Rac-BNP/SnO 2 ), the chiral catalysts exhibited a pronounced shift in product selectivity from CO toward formate production. Importantly, in-situ attenuated total reflectance surface-enhanced infrared absorption spectroscopy (ATR-SEIRAS) reveals that the chiral interface selectively stabilizes the O-bound *OCHO intermediate associated with the formate pathway and modulates interfacial water structure and hydrogen-bonding dynamics. These findings demonstrate that spin-polarized interfacial electron transfer can regulate CO 2 RR pathway selectivity by modulating the stabilization of key intermediates. More broadly, this work establishes chiral spin-selective interfaces as a new strategy for regulating competitive electrocatalytic reaction pathways.

14 SOLAR ENERGY↗
DOE OSTI2026

Emissivity measurements of CuCrZr alloy

The Facility for Rare Isotope Beams (FRIB) heavy-ion accelerator, in user operation since 2022, produces rare isotope beams via interactions of high-intensity stable ion beams with a graphite production target. Approximately 20–40% of the primary beam power is deposited in the target, while the remaining 60–80% is absorbed by the beam dump. The minichannel beam dump (MCBD), currently operated at 20 kW and designed for operation up to 50 kW, uses CuCrZr alloy absorber plates. Thermal validation and thermal cycling tests of the MCBD were conducted at the Applied Research Laboratory (ARL) at Pennsylvania State University. Temperature measurements were obtained from an infrared (IR) camera. Since accurate temperature determination requires reliable emissivity values, the emissivity of CuCrZr was measured using the IR camera validated against thermocouple reference temperatures up to approximately 650 °C. The measurements were conducted under a vacuum level of approximately 10 -5 torr to minimize emissivity variations due to surface oxidation. The emissivity of CuCrZr was determined to be 0.057 ± 0.009 using a constant fit to the measured data over the surface temperature range from 100–650 °C.

Accelerator Subsystems and Technologies↗
DOE OSTI2025

AUTOIGNITION DELAY TIMES FOR REFORMATE GAS MIXTURES FROM METHANE GAS ENGINES

Methane slip is a prominent issue in natural gas reciprocating engines that are used in transportation and marine applications. The incomplete combustion that results in methane slip can be resolved with the introduction of hydrogen within the combustion mixture to improve methane oxidation and further enable combustion within the engine crevices where methane has previously remained unreacted. Steam methane reforming (SMR) is a common method used to produce hydrogen and can be used to design an onboard device to reduce methane slip from reciprocating engines. The development of this reformer device requires the validation of high-fidelity chemical kinetic models at the low temperatures of the crevice volumes of these engines. In this work, auto-ignition data is obtained using a shock tube at lean (φ—0.714 or λ—1.4) and stoichiometric (φ, λ = 1) equivalence ratios spanning a temperature range of 1042–1234 K at the 80-bar operating pressure of the test engine. Blends of methane, hydrogen, and reformate products from the SMR reaction are shock-heated in synthetic air, with the ignition delay time measured using an OH* chemiluminescence detector at 310 nm and a CH* detector at 430 nm. The experimental results are compared to several state-of-the-art chemical kinetic mechanisms from the literature. In general, most of the mechanisms show very good agreement with experiments at higher temperatures, with simulation results showing little deviation from experiments at lower temperatures. A sensitivity analysis was conducted, and the results reveal that the reaction H2 + CH3O2 = H + CH3O2H has a very significant role in determining low-temperature ignition delay times (IDTs) of SMR mixtures. These findings provide valuable insights into the chemical kinetics governing methane reformate combustion and contribute to the optimization of onboard reformer designs aimed at mitigating methane slip in natural gas-fueled engines.

Fraze, Matthew↗
DOE OSTI2026

NanoPSD: A software for automatic detection of Nano-Particle Shape Distribution in electron microscopy images

Accurate quantification of the size and morphology of nanoparticles from electron microscopy (EM) images is essential to understand growth mechanisms, surface reactivity, and functional behavior in nanoscale materials. Manual analysis remains slow, subjective, and difficult to reproduce in large datasets. We introduce NanoPSD (Nano-Particle Shape Distribution), an open-source and fully automated framework for quantitative particle detection and morphology analysis from EM images. NanoPSD integrates adaptive contrast enhancement, polarity-agnostic scale-bar detection, Optical Character Recognition (OCR)-based calibration, and classical segmentation via Otsu thresholding with morphological refinement. Particle contours are used to extract geometric descriptors, including equivalent circular diameter, aspect ratio, circularity, and solidity, enabling automated classification into spherical, rod-like, and aggregate morphologies. The framework supports both single-image and batch processing, generating publication-quality visualizations, LaTeX-ready tables, and structured comma-separated values (CSV) datasets. As a demonstration, we applied NanoPSD to plasma-synthesized nanoparticle samples diagnosed via transmission electron microscopy (TEM). The code produced statistically robust size and morphology distributions spanning a few to tens of nanometers with minimal user supervision. The pipeline demonstrates high reproducibility and scalability, processing large image collections with consistent calibration and output formatting. Its modular design enables seamless integration of future deep-learning-based segmentation models, providing a pathway toward intelligent, data-driven electron microscopy analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Decorated Honeycomb Lattice and Successive Magnetic Transitions of the Delafossite Derivative K 4 Ni 5 Te 3 O 16

Here, we report the synthesis and characterization of the delafossite derivative K 4 Ni 5 Te 3 O 16 , which hosts a decorated honeycomb Ni 2+ lattice. Single-crystal X-ray diffraction, electron microscopy, and neutron diffraction establish a noncentrosymmetric Pmm2 crystal structure that combines regular honeycomb units with elongated motifs formed by Ni 2+ trimers along extended edges. Magnetic susceptibility reveals strong antiferromagnetic interactions, while heat capacity measurements identify two successive magnetic transitions at 30 and 12 K. Neutron diffraction shows that these transitions correspond to a progression from partial to full long-range magnetic order on distinct Ni sites, reflecting competition among the various exchange interactions. The resulting magnetic point group, mm2.1′, inherits the lattice polarity and permits high-order magnetoelectric coupling. Consistent with this symmetry, a quadratic magnetic-field-induced polarization is experimentally observed below the magnetic ordering temperature, likely arising from exchange-driven magnetostriction coupled to the polar lattice. These results establish K 4 Ni 5 Te 3 O 16 as an interesting platform for engineering B-site ordering in delafossite derivatives to realize complex quantum magnetism and functional responses.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Multitask graph neural networks for elastoplastic response prediction in dual-phase polycrystals

Microstructure-sensitive prediction of elastoplastic response remains a recurring bottleneck in multiscale damage and fatigue modeling, where large ensembles of statistically distinct polycrystals are required to quantify variability and extreme-value behavior. In this work, we develop a multitask graph neural network (GNN) surrogate that maps dual-phase ferrite–martensite polycrystal microstructures to Statistical Volume Element (SVE)-level elastoplastic Quantities of Interest (QoIs). Each SVE is represented as a grain-adjacency graph, with node features encoding phase, geometry, and crystallographic orientation, and edge features encoding relative misorientation. A message-passing graph convolution generates node embeddings, which are pooled into a graph representation and passed to a multitask regression head that jointly predicts 10 scalar QoIs and vector-valued stress–strain responses in orthogonal loading directions across multiple martensite volume fractions and SVE sizes. Results show high accuracy for scalar QoIs and strong agreement for full stress–strain trajectories, with population envelopes reproducing both median behavior and finite-SVE variability across compositions and partition scales. A unified model trained on pooled volume-fraction data preserves most within-regime accuracy relative to regime-specific models while also capturing the broader cross-regime variation reflected in the pooled test set. Distributional comparisons further demonstrate that the surrogate preserves heterogeneity under SVE partitioning, enabling statistically consistent block-wise random-field construction for mesoscale analyses. Overall, the proposed grain-graph surrogate provides a practical pathway to accelerate ensemble-based studies of SVE-level constitutive variability in dual-phase polycrystals.

Crystal plasticity↗