Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Composite models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

KRAS4a and KRAS4b show distinct lipid-dependent regulation of RAS-RAF membrane dynamics

KRAS4a and KRAS4b are important regulators of signaling, and their interactions with the plasma membrane are dynamic and influenced by lipid composition. KRAS 4a and 4b have nearly identical globular domains but differ in their membrane-associated hyper variable region (HVR). The functional distinctions between these isoforms remain unclear, particularly with regards to their dependence on specific lipids and the membrane environment. Previous work showed that the membrane orientation of KRAS4b affects its ability to bind to RAF kinase RBDCRD and that the KRAS–RBDCRD complex adopts different poses on the membrane as well as influences the size and composition of the lipid environment. To model differences between KRAS 4a and 4b protein–lipid interactions, we extended the Multiscale Machine-Learned Modeling Infrastructure (MuMMI) to incorporate continuum simulations in the grand canonical ensemble, enabling sampling across macroscopic, coarse-grained, and all-atom resolutions. Using this framework, we systematically altered PIP2 concentrations, KRAS 4a versus 4b, and RAF RBDCRD complexation to assess impacts on membrane–protein interactions and dynamics. Our results reveal that reducing PIP2 shifts and broadens the membrane orientational preference of both KRAS 4b and 4a, with stronger effects on 4b HVR localization versus 4a. We demonstrate that with depletion of the strong negatively charged PIP2 lipid, the less charged phosphatidylserine replaces PIP2. Our findings highlight similarities and distinctions in the dynamics and lipid dependency of KRAS isoforms and suggest that ordering of the local lipid composition by HVRs is a shared property and key modulator of RAS-mediated signaling at the plasma membrane.

Biological and medical sciences↗
DOE OSTI2026

Aviation Fuel Characterization at Operationally Relevant Conditions

To accelerate approval and potentially expand the allowable property range for aviation fuels, we are using high performance computing simulations to reveal fuel property effects on aviation combustor performance. These simulations are supported by fuel property measurements over temperatures and pressures that the fuel experiences in an aircraft engine and by validated chemical kinetics models for SAF combustion. Here we report density, viscosity, and surface tension results for conventional jet fuel and multiple synthetic fuels from -30 degrees Celsius to 200 degrees Celsius (-40 degrees Celsius for viscosity) and 1 atm to 70 atm including an assessment of method repeatability. Properties of surrogate mixtures are also investigated. Distillation, ICN, LHV, flashpoint, and Cp are also reported, and data are being used to develop models to predict fuel properties from composition (GCxGC).

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Atomistic models of In and Ga diffusion in Cu(In,Ga)Se 2

To improve the performance of Cu(In,Ga)Se 2 thin-film photovoltaic devices, a robust understanding of the dominant diffusion pathways of the alloy species In and Ga is needed. Here, the most probable defect complexes and mechanisms for In and Ga diffusion are identified with the aid of density functional theory. The binding energies and migration barriers for these complexes are calculated in bulk CuInSe 2 and CuGaSe 2 . Analytic models and kinetic lattice Monte Carlo simulations are employed to predict the diffusivity of In and Ga under variations in composition and temperature. Here, we find that a model based on coulombic interactions between group III antisites and vacancies on the Cu-sublattice produces results that match well with experiment.

14 SOLAR ENERGY↗
DOE OSTI2022

Fatty acid capped, metal oxo clusters as the smallest conceivable nanocrystal prototypes

Metal oxo clusters of the type M 6 O 4 (OH) 4 (OOCR) 12 (M = Zr or Hf) are valuable building blocks for materials science. Here, we synthesize a series of zirconium and hafnium oxo clusters with ligands that are typically used to stabilize oxide nanocrystals (fatty acids with long and/or branched chains). The fatty acid capped oxo clusters have a high solubility but do not crystallize, precluding traditional purification and single-crystal XRD analysis. We thus develop alternative purification strategies and we use X-ray total scattering and Pair Distribution Function (PDF) analysis as our main method to elucidate the structure of the cluster core. We identify the correct structure from a series of possible clusters (Zr3, Zr4, Zr6, Zr12, Zr10, and Zr26). Excellent refinements are only obtained when the ligands are part of the structure model. Further evidence for the cluster composition is provided by nuclear magnetic resonance (NMR), infrared spectroscopy (FTIR), thermogravimetry analysis (TGA), and mass spectrometry (MS). We find that hydrogen bonded carboxylic acid is an intrinsic part of the oxo cluster. Using our analytical tools, we elucidate the conversion from a Zr6 monomer to a Zr12 dimer (and vice versa), induced by carboxylate ligand exchange. Finally, we compare the catalytic performance of Zr12-oleate clusters with oleate capped, 5.5 nm zirconium oxide nanocrystals in the esterification of oleic acid with ethanol. The oxo clusters present a five times higher reaction rate, due to their higher surface area. Since the oxo clusters are the lower limit of downscaling oxide nanocrystals, we present them as appealing catalytic materials, and as atomically precise model systems. In addition, the lessons learned regarding PDF analysis are applicable to other areas of cluster science as well, from semiconductor and metal clusters, to polyoxometalates.

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

Stiff and strong, lightweight bi-material sandwich plate-lattices with enhanced energy absorption

Plate-based lattices are predicted to reach theoretical Hashin–Shtrikman and Suquet upper bounds on stiffness and strength. However, simultaneously attaining high energy absorption in these plate-lattices still remains elusive, which is critical for many structural applications such as shock wave absorber and protective devices. In this article, we present bi-material isotropic cubic + octet sandwich plate-lattices composed of carbon fiber-reinforced polymer (stiff) skins and elastomeric (soft) core. This bi-material configuration enhances their energy absorption capability while retaining stretching-dominated behavior. We investigate their mechanical properties through an analytical model and finite element simulations. Our results show that they achieve enhanced energy absorption approximately 2–2.8 times higher than their homogeneous counterparts while marginally compromising their stiffness and strength. When compared to previously reported materials, these materials achieve superior strength-energy absorption characteristics, making them an excellent candidate for stiff and strong, lightweight energy absorbing applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Interface pinning causes the hysteresis of the hydride transformation in binary metal hydrides

Hydriding and dehydriding transitions in bulk and nanocrystalline binary metal hydrides were studied using the Pd-H model system by measuring pressure-composition isotherms with in situ x-ray diffractometry. Nanocrystalline Pd showed a smaller pressure hysteresis, solvus hysteresis, and hysteresis in lattice parameter, compared to bulk Pd. The time-dependence of pressure equilibration was measured after dosing with aliquots of hydrogen, giving equilibration times that were much faster in the single-phase regions than in the two-phase plateaus. In the broad two-phase plateaus, the pressure relaxations were exponential functions of time. An explanation of hysteresis is developed that is based on a dissipative potential barrier that impedes the motion of the interface due to interactions between lattice defects and the two-phase interface. The exponential pressure relaxations and hysteresis are consistent for this mechanism. For a simple model of the pinning potential, the potential barrier maximum is an order of magnitude less than typical grain boundary energies. These pinning effects are substantially different in the nanocrystalline Pd, suggesting differences in the hydriding mechanism.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Understanding solidification of near eutectic alloy using Cellular Automata (CA)

The solidification microstructure of an alloy strongly influences the mechanical properties. Dendritic and eutectic solidification are two important pathways of alloy solidification. The formation of phases, and their size and distribution depend on which pathway to take. In contrast to the Scheil model that has been widely used in describing the solidification pathways of alloys manufactured from conventional casting approach, more sophisticated solidification models are required to describe the state-of-the art manufacturing processes, such as additive manufacturing where distinctly different phase morphology, size and distribution than conventional casting can form, due to the vast difference in cooling conditions. Therefore, this work is to develop a computational Cellular Automata framework which includes the modeling of nucleation and growth of dendritic and eutectic solidification, as well as their competition as a function of alloy composition, undercooling and cooling rate. Further, after individual solidification models were validated against analytical solutions, the models were then combined to predict the competition between dendritic and eutectic solidification in alloys with off-eutectic compositions. This work predicts increased cooling rates suppress the dendritic solidification and promote eutectic solidification. Although quantitative validation is needed, this finding is consistent with qualitative observation in literature.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Smart culture medium optimization for recombinant protein production: Experimental, modeling, and AI/ML-driven strategies

Recombinant protein production (RPP) is central to biotechnology, where recombinant proteins are used as either end products or catalysts in the synthesis of chemicals, fuels, and materials. Among the major cost drivers, culture medium plays a pivotal role in determining protein yield and quality. This review presents a comprehensive perspective on the critical stages of “smart” culture medium optimization: planning, screening, modeling, optimization, and validation. In the planning stage, we examine the nutritional and energetic roles of medium components, including carbon, nitrogen, amino acids, salts, and trace metals, and their impacts on culture parameters such as pH, oxidative state, and osmolality. We highlight the variability in trace metal content due to water sources, culture vessels, and raw materials, which can substantially influence RPP. The screening stage covers Design of Experiments (DoE) approaches, assessing their theoretical basis, implementation, and limitations. For modeling, we describe methods that integrate experimental data to develop predictive models for smart medium formulation. Model-based optimization strategies can then be employed to select optimal media compositions for a given application. The validation stage aims to evaluate model predictions and provide feedback for model training and refinement. Finally, we survey mechanistic and artificial intelligence/machine learning (AI/ML)-driven models as integrated, transformational tools for predictive modeling of bioprocess conditions, nutrient availability, cellular metabolism, and protein quality, with the goal of optimizing culture media to enhance protein yields while reducing costs and environmental impact. We conclude by addressing the challenges of translating laboratory-scale medium optimization to industrial-scale settings and exploring future AI/ML-driven approaches that may overcome current bottlenecks and accelerate medium design for RPP. Overall, this review provides a unified framework for advancing smart medium design in RPP.

Artificial Intelligence/Machine Learning (AI/ML)↗
DOE OSTI2024

On the applicability of the Redlich-Kister framework for viscosity estimation of molten halide salt mixtures

For molten halide salt mixtures already being utilized or under consideration for carbon-free energy production systems, it is crucial that their viscosity is well understood so that system thermal hydraulics can be reliably assessed. Because of the difficulty in accurately measuring molten halide viscosity and the sheer size of the matrix of possible higher order salt mixtures that may be of interest to the energy industry, there are several gaps in the quantified understanding of molten halide viscosity across this matrix. As such, both first-principles and semi-empirical modeling techniques may be crucial for rapidly assessing this broad, complex compositional domain. Herein, the Redlich-Kister framework is applied to assess the feasibility of broadly interpolating and estimating the viscosity of several pseudobinary and pseudoternary molten halide salt systems that may be of key interest to the energy industry. The framework is based on the assumption that an ideal component and a nonideal component collectively describe the viscosity as a function of composition and temperature for a given molten halide system. Three different ideal models were considered for the ideal component, including Grunburg-Nissan, Katti-Chaudhri, and Gambill methods. Regarding the pseudobinary interpolations, the Redlich-Kister models with either the Grunburg-Nissan or Katti-Chaudhri models as the ideal component resulted in either highly (average error less than 5%) or reasonably (average error less than 15%) accurate interpolations of pseudobinary halide viscosity; BeF 2 - or UF 4 -bearing salts tended to result in reasonably accurate interpolations, whereas other pseudobinary mixtures tended to show high accuracy. Regarding the pseudoternary extrapolations, the Redlich-Kister framework shows reasonable success at estimating the extent to which a pseudoternary system may indicate deviations from ideal Grunburg-Nissan mixing, where discrepancies with comparative experimental data generally stay within 30%. Finally, the primary reasons identified for such discrepancies are (1) inaccuracy in the underlying experimental data, (2) different complexation behavior in the higher order systems compared to the pseudobinary subsystems, and (3) extrapolation into temperatures too far out of the domain, which is valid for the underlying experimental data feeding the Redlich-Kister model.

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

Atomic-Level Features for Kinetic Monte Carlo Models of Complex Chemistry from Molecular Dynamics Simulations

The high computational cost of evaluating atomic interactions recently motivated the development of computationally inexpensive kinetic models, which can be parametrized from MD simulations of complex chemistry of thousands of species or other processes and accelerate the prediction of the chemical evolution by up to four order of magnitude. Such models go beyond the commonly employed potential energy surface fitting methods in that they are aimed purely at describing kinetic effects. So far, such kinetic models utilize molecular descriptions of reactions and have been constrained to only reproduce molecules previously observed in MD simulations. Therefore, these descriptions fail to predict the reactivity of unobserved molecules, for example in the case of large molecules or solids. In this work, we propose a new approach for the extraction of reaction mechanisms and reaction rates from MD simulations, namely the use of atomic-level features. Using the complex chemical network of hydrocarbon pyrolysis as example, it is demonstrated that kinetic models built using atomic features are able to explore chemical reaction pathways never observed in the MD simulations used to parametrize them, a critical feature to describe rare events. Atomic-level features are shown to construct reaction mechanisms and estimate reaction rates of unknown molecular species from elementary atomic events. Through comparisons of the model ability to extrapolate to longer simulation timescales and different chemical compositions than the ones used for parameterization, it is demonstrated that kinetic models employing atomic features retain the same level of accuracy and transferability as the use of features based on molecular species, while being more compact and parametrized with less data. We also find that atomic features can better describe the formation of large molecules enabling the simultaneous description of small molecules and condensed phases.

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

Higher Dimensionality in the Mg–Co–B System: Synthesis and Structure of Incommensurate Composite Mg 1+ε Co 4 B 4

Guided by high-temperature in situ X-ray diffraction, the discovery and synthesis of Mg 1+ε Co 4 B 4 (ε ≈ 0.272) using a MgH 2 hydride precursor is reported, along with a detailed crystal structure description and measurement of magnetic properties. The mismatch in lattice periodicities between Mg and Co–B substructures places Mg 1+ε Co 4 B 4 in the family of incommensurate composite crystals and prompted structural refinement in a (3 + 1)-dimensional model. The structure of Mg 1+ε Co 4 B 4 (P4 2 /ncm(00γ)s00s, a = 6.75847(7) Å, c = 3.94007(8) Å, q = (0, 0, 1.2721(3))) was refined from neutron powder diffraction and high-resolution powder X-ray diffraction data and confirmed by scanning transmission electron microscopy and electron diffraction. Mg 1+ε Co 4 B 4 is isostructural to Nd 1+ε Fe 4 B 4 and several related ternary borides with 0.07 ≤ ε ≤ 0.17, with Mg occupying the rare-earth site. Satellite reflections in the electron diffraction patterns hinted at positional modulation of the transition metal–boron substructure by Mg atoms, but this could not be refined from the neutron or X-ray diffraction data. Low-temperature magnetic measurements show no indications of long-range magnetic ordering or superconductivity down to 5 K. DFT calculations confirmed the absence of a magnetically ordered ground state and the stability of a 5:4 supercell (ε = 0.25) relative to the fully commensurate structure. Neutron diffraction and synthesis from elemental Mg demonstrated that Mg 1+ε Co 4 B 4 is not a hydrogen-stabilized phase. Mg 1+ε Co 4 B 4 represents the second compound reported in the Mg–Co–B system and the first superspace symmetry model of a Nd 1+ε Fe 4 B 4 -type incommensurate composite compound refined from powder diffraction data.

chemical structure↗
DOE OSTI2022

Functional Traits Resolve Mechanisms Governing the Assembly and Distribution of Nitrogen-Cycling Microbial Communities in the Global Ocean

Microorganisms drive much of the marine nitrogen (N) cycle, which jointly controls the primary production in the global ocean. However, our understanding of the microbial communities driving the global ocean N cycle remains fragmented. Focusing on “who is doing what, where, and how?”, this study draws a clear picture describing the global biogeography of marine N-cycling microbial communities by utilizing the Tara Oceans shotgun metagenomes. The marine N-cycling communities are highly variable taxonomically but relatively even at the functional trait level, showing clear functional redundancy properties. The functional traits and taxonomic groups are shaped by the same set of geo-environmental factors, among which, depth is the major factor impacting marine N-cycling communities, differentiating mesopelagic from epipelagic communities. Latitudinal diversity gradients and distance-decay relationships are observed for taxonomic groups, but rarely or weakly for functional traits. The composition of functional traits is strongly deterministic as revealed by null model analysis, while a higher degree of stochasticity is observed for taxonomic composition. Integrating multiple lines of evidence, in addition to drawing a biogeographic picture of marine N-cycling communities, this study also demonstrated an essential microbial ecological theory—determinism governs the assembly of microbial communities performing essential biogeochemical processes; the environment selects functional traits rather than taxonomic groups; functional redundancy underlies stochastic taxonomic community assembly.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Microstructure Scale Lithium-Ion Battery Modeling: Part I. On Through-Plane Heterogeneity, Impact of Mesh Representation, and Differences between Macro- and Microscale Models

Li-ion battery performance and degradation are strongly correlated with the electrode microstructures and can be modeled at different scales, each with their own limitations. Herein, we compare predictions achieved with a macro- and a micro-scale model, that is, respectively, neglecting or considering the microstructural heterogeneity of the composite electrodes, on virtual numerically generated and real microstructures. While both models are in relative agreement at the low charge rates, differences arise for fast charging scenarios and especially for the real, highly heterogenous, microstructures. The microscale model predicts that electrolyte concentration saturation and depletion, respectively, at the back of the cathode and of the anode are exacerbated, and that lithium plating occurs earlier for real microstructures. The present work also indicates that the mesh representation significantly impacts the microscale model predictions, and consequently that microscale models should add surface area as a parameter to consider explicitly surface roughness. This article is the first of a series, with subsequent entries further investigating in-plane heterogeneities, lithium plating, and the impact of microstructure representativity on model predictions.

25 ENERGY STORAGE↗
DOE OSTI2023

TRACER-Vertical Profile of Aerosol in Convective System (TRACER-VPACS) Field Campaign Report

We investigated physicochemical properties of aerosol at different altitudes during the U.S. Department of Energy (DOE)’s Atmospheric Radiation Measurement (ARM) user facility Tracking Aerosol Convection Interactions Experiment (TRACER). We used ARM’s tethered balloon system (TBS) and multi-modal offline analysis techniques available at DOE’s Environmental Molecular Sciences Laboratory (EMSL) at Pacific Northwest National Laboratory (PNNL) and at the Advanced Light Source (ALS) at the Lawrence Berkley National Laboratory (LBNL). The main scientific aim of this project was to understand the influence of the vertical profile of aerosol composition on the properties of deep convective clouds. The effect of aerosol properties on convective activity is still poorly understood and debated. For example, previous studies suggested that higher concentration of hygroscopic aerosol may lead to lower collection efficiency and eventually delay in onset of precipitation and total amount of precipitation. However, some observations suggested that increased aerosol concentration may result in vigorous convection and a higher amount of precipitation. In addition, the radiative forcing of the ambient aerosol is strongly influenced by the vertical distribution of aerosols and their chemical composition, and not well constrained in climate models. The estimation of the indirect effect of aerosol also has a significant uncertainty due to the spatial variation of the aerosol, especially at the cloud base. Detailed size-resolved chemical composition studies of aerosol and their microphysical properties, along with in situ measurements from ARM’s mobile facilities, help to improve our understanding of the role of aerosol in deep convective cloud formation over a broad range of environmental thermodynamic and aerosol conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Impact of dry intrusion events on the composition and mixing state of particles during the winter Aerosol and Cloud Experiment in the Eastern North Atlantic (ACE-ENA)

Abstract. Long-range transport of continental emissions has a far-reaching influence over remote regions, resulting in substantial change in the size, morphology, and composition of the local aerosol population and cloud condensation nuclei (CCN) budget. Here, we investigate the physicochemical properties of atmospheric particles collected on board a research aircraft flown over the Azores during the winter 2018 Aerosol and Cloud Experiment in the Eastern North Atlantic (ACE-ENA) campaign. Particles were collected within the marine boundary layer (MBL) and free troposphere (FT) after long-range atmospheric transport episodes facilitated by dry intrusion (DI) events. Chemical and physical properties of individual particles were investigated using complementary capabilities of computer-controlled scanning electron microscopy and X-ray spectromicroscopy to probe particle external and internal mixing state characteristics. Furthermore, real-time measurements of aerosol size distribution, cloud condensation nuclei (CCN) concentration, and back-trajectory calculations were utilized to help bring into context the findings from offline spectromicroscopy analysis. While carbonaceous particles were found to be the dominant particle type in the region, changes in the percent contribution of organics across the particle population (i.e., external mixing) shifted from 68 % to 43 % in the MBL and from 92 % to 46 % in FT samples during DI events. This change in carbonaceous contribution is counterbalanced by the increase in inorganics from 32 % to 57 % in the MBL and 8 % to 55 % in FT. The quantification of the organic volume fraction (OVF) of individual particles derived from X-ray spectromicroscopy, which relates to the multi-component internal composition of individual particles, showed a factor of 2.06 ± 0.16 and 1.11 ± 0.04 increase in the MBL and FT, respectively, among DI samples. We show that supplying particle OVF into the κ-Köhler equation can be used as a good approximation of field-measured in situ CCN concentrations. We also report changes in the κ values in the MBL from κMBL, non-DI=0.48 to κMBL, DI=0.41, while changes in the FT result in κFT, non-DI=0.36 to κFT, DI=0.33, which is consistent with enhancements in OVF followed by the DI episodes. Our observations suggest that entrainment of particles from long-range continental sources alters the mixing state population and CCN properties of aerosol in the region. The work presented here provides field observation data that can inform atmospheric models that simulate sources and particle composition in the eastern North Atlantic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Improving ChemCam LIBS long-distance elemental compositions using empirical abundance trends

The ChemCam instrument on the Curiosity rover provides chemical compositions of Martian rocks and soils using remote laser-induced breakdown spectroscopy (LIBS). The elemental calibration is stable as a function of distance for Ti, Fe, Mg, and Ca. The calibration shows small, systematically increasing abundance trends as a function of distance for Al, Na, K, and to some extent, Si. The distance effect is known to be due to a dependence with distance on the relative strengths of atomic transition lines. Emission lines representing transitions from relatively low energy levels remain intense at longer distances while emission lines representing transitions from higher energy levels decrease in intensity more rapidly as a function of distance. The multivariate algorithms used to determine elemental compositions rely on a large number of emission lines in many cases, so rather than trying to correct all emission lines, a study was made of the predicted compositions as a function of distance, in order to determine an empirical correction. Abundance trends can be well approximated by a linear trend with distance within the ranges of abundances and distances observed up to ~6 m. Data from 11 distinct geological members and data groups of the Murray formation in Gale crater, Mars, were used to form the model, selecting the members and data groups yielding the best statistics. The model was tested using data from several targets observed from two different distances, and using data from the Kimberley formation, the composition of which is significantly different from the Murray formation, showing that the model works on other compositions beyond those used to build the model. For long-distance observations up to ~6 m, corrections can be made back to an equivalent composition at the median distance of ChemCam observations (2.6 m). Finally, the model has been validated up to 6.2 m, although ChemCam is able to observe bedrock targets to >7 m, and iron meteorites to distances of >9 m.

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

Neural network kinetics for exploring diffusion multiplicity and chemical ordering in compositionally complex materials

Diffusion involving atom transport from one location to another governs many important processes and behaviors such as precipitation and phase nucleation. The inherent chemical complexity in compositionally complex materials poses challenges for modeling atomic diffusion and the resulting formation of chemically ordered structures. Here, we introduce a neural network kinetics (NNK) scheme that predicts and simulates diffusion-induced chemical and structural evolution in complex concentrated chemical environments. The framework is grounded on efficient on-lattice structure and chemistry representation combined with artificial neural networks, enabling precise prediction of all path-dependent migration barriers and individual atom jumps. To demonstrate the method, we study the temperature-dependent local chemical ordering in a refractory NbMoTa alloy and reveal a critical temperature at which the B2 order reaches a maximum. The atomic jump randomness map exhibits the highest diffusion heterogeneity (multiplicity) in the vicinity of this characteristic temperature, which is closely related to chemical ordering and B2 structure formation. The scalable NNK framework provides a promising new avenue to exploring diffusion-related properties in the vast compositional space within which extraordinary properties are hidden.

36 MATERIALS SCIENCE↗
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↗