Search NASA⌕ Search

SEARCH · Search NASA

Results for “multicomponent system”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Model Data Archive Associated with Manuscript "Fire-altered Carbon Pools Create Disturbance Memory in Stream Dissolved Organic Carbon"

This data package supports the publication “Fire-altered Carbon Pools Create Disturbance Memory in Stream Dissolved Organic Carbon” by Li et al. (2026). The package contains processed model inputs, configuration files, restart files, simulation outputs, scripts, and visualization products used to evaluate post-fire dissolved organic carbon (DOC) dynamics in the Naches River Watershed, Washington, USA, following the 2021 Schneider Springs Fire. The modeling workflow couples ELM-BGC, the biogeochemistry-enabled Energy Exascale Earth System Model Land Model; ATS, the Advanced Terrestrial Simulator for integrated surface-subsurface hydrology; and PFLOTRAN, a reactive transport model for multicomponent aqueous geochemistry. Together, these models simulate how wildfire-induced changes in vegetation, litter, coarse woody debris, and soil organic matter influence DOC production, transport, and reaction from burned hillslopes to stream networks. The archive includes preprocessed meteorological, geospatial, hydrologic, and biogeochemical forcing data; ELM-BGC-derived DOC source terms; ATS mesh files; PFLOTRAN reactive-transport inputs; model configuration files; spin-up and transient restart files; watershed-scale diagnostic outputs; stream concentration time series; and figures or visualization files used to inspect and reproduce key results. File types include Hierarchical Data Format 5 (HDF5) files for gridded forcing and model-coupling data, model input and configuration files for ELM-BGC, ATS, and PFLOTRAN, restart and simulation-output files generated by the modeling workflow, tabular or time-series diagnostic outputs, scripts for post-processing and figure generation, and image or visualization products associated with the manuscript. Use of the package depends on the intended task. Re-running the simulations requires the relevant modeling software, including ELM-BGC, ATS, and PFLOTRAN as ATS's geochemical engine. Inspecting outputs and reproducing figures requires Python with scientific plotting libraries such as Matplotlib, and three-dimensional model outputs may be viewed with ParaView. Geographic information system files or maps may be inspected with ArcGIS Pro or comparable GIS software. The data package is intended to enable traceability, reuse, and partial reproduction of the coupled land-to-watershed hydro-biogeochemical modeling workflow used to test how wildfire disturbance affects terrestrial carbon pools and downstream DOC dynamics.

ATS↗
DOE OSTI2021

Structural and chemical disorder enhance point defect diffusion and atomic transport in Ni 3 Al-based γ' phase

The ordered L1 2 γ' Ni 3 Al intermetallic inclusions are routinely used as strengthening constituents in Ni-based superalloys and recent high-entropy alloys. As there are different metallic elements relative to the alloy matrix, the composition of γ' phase inclusions usually suffers some variations. In this work, we studied the elemental partition and its influence on atomic transport properties in a model (NiCo) 3 Al γ' phase with chemical and structural disorder. Combining the Monte Carlo and molecular static techniques, we first determined the elemental distribution in such compositionally-complexed γ' phase. The results suggest that Co tends to segregate, which profoundly influences defect energetics and diffusion properties. By modeling thermally activated diffusion by molecular dynamics, we found that chemical disorder induced by Co can enhance vacancy diffusion while significantly suppress diffusion via interstitial atom migration, as compared to the stoichiometric γ' Ni 3 Al. We further reveal that these observations are closely related to the defect energetics in these systems. Furthermore, our results explain the effects of disorder, long- and short-range orders on the diffusion properties in multicomponent intermetallic alloys, which helps to understand the stability and, thus, γ'-phase-assisted strengthening at different conditions including annealing and irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Glass formation during combinatorial sputtering in binary alloys

Glass formation is a complex phenomenon influenced by thermodynamic and kinetic aspects, which are often controlled by extrinsic contributions. While bulk metallic glasses are typically multicomponent alloys, binary alloys offer a simplified approach to studying glass formation. In this study, we fabricated 57 binary alloy systems through combinatorial sputtering, where each alloy system is represented in 66 different alloys. We developed an automated analysis to determine structure and composition using X-ray diffraction and energy-dispersive X-ray spectroscopy for over 3700 alloys. We found that ∼17 % of the alloys form glasses under the estimated cooling rate during sputtering of ∼10 8 K/s. Data analysis revealed that commonly used factors like atomic size ratio and heat of mixing are ineffective in predicting glass formation. However, the crystal structure mismatch of the alloys’ elements emerged as the strongest indicator of glass formation under sputtering conditions of binary alloys. Here, the differences in glass formation under slow cooling rates used for bulk glass formation and the here observed glass formation under rapid cooling rates are discussed.

Binary alloys↗
DOE OSTI2023

Predicting short-range order evolution in WTaCrVHf refractory high-entropy alloys

Short-range order (SRO) in multicomponent concentrated alloys affects their mechanical response. Hence, is paramount to understand how composition modifies the chemical ordering in the system to design materials with optimal properties. Here, in this work, we present a methodology to predict the SRO and thermodynamic properties in chemically complex systems and apply it to the WTaCrVHf quinary alloy. We observe that the addition of Hf significantly modifies the SRO, mainly at intermediate to low temperatures, matching experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Detection of non‐native species formed during fibrillization of the myocilin olfactomedin domain

Abstract Glaucoma is a group of neurodegenerative diseases that together are the leading cause of irreversible blindness worldwide. Myocilin‐associated glaucoma is an inherited form of this disease, caused by intracellular aggregation of misfolded mutant myocilin. In vitro, the myocilin C‐terminal olfactomedin domain (OLF), the relevant domain for glaucoma pathogenesis, can be driven to form amyloid‐like fibrils under mild conditions. Here we characterize a species present during in vitro fibrillization. Purified OLF was subjected to fibrillization at concentrations required for downstream electron microscopy imaging and NMR spectroscopy. Additional biophysical techniques, including analytical ultracentrifugation and X‐ray crystallography, were employed to further characterize the multicomponent mixture. Negative stain transmission electron microscopy (TEM) shows a non‐native species reminiscent of known prefibrillar oligomers from other amyloid systems, NMR indicates a minor population of partially misfolded species is present in solution, and cryo‐EM imaging shows two‐dimensional protein arrays. The predominant soluble species remaining in solution after the fibril reaction is natively folded, as evidenced by X‐ray crystallography. In summary, after incubating OLF under fibrillization‐promoting conditions, there is a heterogeneous mixture consisting of soluble folded protein, mature amyloid‐like fibrils, and partially misfolded intermediate species that at present belie additional molecular detail. The characterization of OLF fibrillar species illustrates the challenges associated with developing a comprehensive understanding of the fibrillization process for large, non‐model amyloidogenic proteins.

Scelsi, Hailee F. [School of Chemistry and Biochem↗
DOE OSTI2023

PFLOTRAN 5

PFLOTRAN leverages massively parallel, high performance computing to simulate large-scale non-isothermal multiphase flow, multicomponent reactive transport and electrical resistivity tomography (ERT) problems in the subsurface environment. Researchers have employed PFLOTRAN to simulate these Earth system processes on leadership class supercomputers for over two decades. The code is designed to predict the future estate of environmental systems and better inform stakeholders in the regulatory decision making process (e.g., fate of contaminants, long-term stewardship for nuclear waste, impact of climate change, etc.). A diverse team of scientists oversees PFLOTRAN development and maintenance under an open-source licensing agreement and manages contributions from an international community of researchers.

Hammond, Glenn↗
DOE OSTI2022

Closures for multicomponent reacting flows based on dispersion analysis

This work presents algebraic closure models associated with advective transport and nonlinear reactions in a Reynolds-averaged Navier-Stokes context for a system of species subject to binary reactions and transport by advection and diffusion. Expanding upon analysis originally developed for non-reactive transport in the context of Taylor dispersion of scalars, this work extends the modified gradient diffusion model explicated by Peters [N. Peters, Turbulent Combustion, Cambridge Monographs on Mechanics (Cambridge University Press, Cambridge, 2000)] and based on work by Corrsin [S. Corrsin, The reactant concentration spectrum in turbulent mixing with a first-order reaction, J. Fluid Mech. 11, 407 (1961)] beyond single-component transport phenomena and involving nonlinear reactions. The presented model forms, from this weakly nonlinear extension of the original dispersion theory, lead to an analytic expression for the eddy diffusivity matrix that explicitly captures the influence of the reaction kinetics on the closure operators. Furthermore, we demonstrate that the derived model form directly translates between flow topologies through a priori and a posteriori testing of a binary species system subject to homogeneous isotropic turbulence. Using two- and three-dimensional direct numerical simulations involving laminar and turbulent flows, it is shown that this framework improves prediction of mean quantities compared to previous results. Lastly, the presented model form, collapses to the earlier gradient diffusion and its modified version derived by Corrsin in the limits of nonreactive species and linear reactions, respectively.

42 ENGINEERING↗
DOE OSTI2021

Self-sorting in supramolecular assemblies

Supramolecular self-assembly enables living organisms to form highly functional hierarchical structures with individual components self-organized across multiple length scales. This has inspired work on multicomponent supramolecular materials to understand factors behind co-assembly versus self-sorting of molecules. We report here on a supramolecular system comprised of negatively charged peptide amphiphile (PA) molecules, in which only a tiny fraction of the molecules (0.7 mol%) were covalently conjugated to one of two different fluorophores, half to fluorescein isothiocyanate (FTIC) and the other half to tetramethylrhodamine (TAMRA). Confocal microscopy of the system revealed self-sorting of the two different fluorescent PA molecules, where TAMRA PA is concentrated in micron-scale domains while FITC PA remains dispersed throughout the sample. From Förster resonance energy transfer and fluorescence recovery experiments, we conclude that conjugation of the negatively charged FITC to PA significantly disrupts its co-assembly with the 99.3 mol% of unlabeled molecules, which are responsible for formation of micron-scale domains. Conversely, conjugation of the zwitterionic TAMRA causes no such disruption. Interestingly, this dissimilar behavior between FITC and TAMRA PA causes them to self-sort at large length scales in the supramolecular system, mediated not by specific interactions among the individual fluorophores but instead by their different propensities to co-assemble with the majority component. We also found that greater ionic strength in the aqueous environment of the system promotes mixing by lowering the electrostatic barriers involved in self-sorting. Furthermore, our results demonstrate great thermodynamic subtlety in the driving forces that mediate self-sorting versus co-assembly in supramolecular peptide assemblies.

36 MATERIALS SCIENCE↗
DOE OSTI2025

FY24 Task 4: Studies of Phosphate and Fluoride Solubility for Dissolution of High Phosphate Tank Waste

Knowledge gaps have been identified in phosphate solubility in almost every single- and multi-component system, and particularly for aluminate in phosphate/hydroxide, where uncertainties in predicted values of aluminum solubility are greater than 50%. Solubility data of multicomponent, aqueous electrolytes containing sodium hydroxide, sodium fluoride, sodium phosphate, sodium nitrite, sodium nitrate, and dissolved gibbsite are also sparse. For example, for solutions of (i) sodium nitrate and sodium phosphate, data are only available at 30 and 50 °C, and for (ii) dissolved gibbsite in sodium hydroxide and sodium phosphate, only two data points are available at 20 and 40 °C. There is no data available on mixtures of sodium hydroxide, sodium phosphate, sodium fluoride, and dissolved gibbsite. These knowledge gaps were identified in a technical review of waste solubility data and the impact of dilution on solution stabilities and will be addressed in this work to predict aluminum hydroxide and sodium phosphate solubility in multicomponent electrolytes upon dilution, and upon variation of temperature. Results will provide the technical basis to develop accurate models for (i) gibbsite solubility and mass transfer of aluminum between solid and liquid forms following the sluicing of sludge and saltcake with water; and (ii) further blending of these suspensions with bismuth from bismuth phosphate waste, and zirconium and uranium left from the fuel decladding. This work will be essential to developing a disposition path for retrieval solutions from bismuth phosphate wastes. This effort would also support sludge washing to further reduce phosphate concentration if that process were to be added back to the flowsheet in the future.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Thermodynamic modeling of adsorption at the liquid-solid interface

Adsorptive separation techniques are significantly energy efficient in comparison to conventional thermal separation techniques such as distillation. Despite extensive research and development activities undertaken for mixed-gas adsorption, the use of adsorption techniques for the separation of multicomponent liquid mixtures is still limited. A major barrier is the lack of accurate adsorption thermodynamic models, which form the scientific foundation of process simulation of such systems, making the translation to industrial scale challenging. In this work, we have rigorously computed the surface excess of adsorption for six binary liquid mixtures on silica gel at 303 K using the frameworks of the generalized Langmuir isotherm (gL) and the adsorbed solution theory (AST). The six binary liquid mixtures were formed by the pair-wise combinations of four components: benzene, 1,2-dichloroethane, cyclohexane, and n-heptane. We have based our calculations by considering simultaneous equilibria of three phases: saturated vapor phase, bulk liquid phase, and adsorbed phase. Further, the composition of the corresponding saturated vapor phase was estimated by the Nonrandom Two-Liquid activity coefficient model and experimental vapor-liquid equilibria data. The activity coefficients of the adsorbed phase, the central issue of multicomponent adsorption thermodynamics, were calculated using the adsorption Nonrandom Two-Liquid activity coefficient model. Devoid of simplifying assumptions, gL and AST provide rigorous thermodynamic frameworks for adsorption equilibria of multicomponent liquid mixtures.

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

Efficient Control of Atom Arrangement in Ternary Metal Chalcogenide Nanoparticles Using Precursor Oxidation State

Controlling both the concentration and the distribution of elements in a given material is often crucial to extracting and optimizing synergistic properties of the various constituents. An interesting class of such multielement materials is metal chalcogenide nanoparticles, which exhibit a wide range of composition-dependent optoelectronic properties including both bandgap-mediated processes and localized surface plasmon resonance properties, each of which is useful in applications ranging from energy conversion to sensing. Because metal chalcogenide nanoparticles can support several different metal elements in a variety of chalcogen lattices, this material class has particularly benefited from the ability to control both atom concentration and atom arrangement to tailor final particle properties. The primary method to access complex, multimetallic chalcogenide particles is via a postsynthetic cation exchange strategy. One-pot syntheses have been less explored to access these complex particles, although this route is desirable for economy and scalability. In this work, we compare the composition and morphology outcomes from cation exchange and one-pot preparation approaches using a Cu/Ag/Se system, which is already known to exhibit both binary and ternary metal chalcogenide phases. We show that at similar concentrations of the two metal cations, initial reaction conditions for the one-pot method yield multicomponent nanoparticles, whereas cation exchange yields homogeneous ternary metal chalcogenide structures. We then show that by tuning the precursor oxidation state for the one-pot method, this approach can be used to access homogeneous ternary metal chalcogenide particles that are similar in atom arrangement to the particles obtained using cation exchange. Taken together, our results demonstrate reliable synthetic methods that yield a variety of controlled compositions and composition morphologies in the Cu/Ag/Se system. Importantly, we demonstrate that this entire collection of architectures can all be accessed via a one-pot method simply by modifying metal precursor chemistry. The mechanistic insights gained and the resulting streamlined syntheses outlined indicate pathways to easily scaled, highly tailorable syntheses for rapid translation into downstream technologies.

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

Mesoporous Thin Film Architectures: Addressing Material Demands through Molecular Self-Assembly

Mesoporous thin films spark interest across a wide range of disciplines due to their tunable nanostructures, large internal surface areas, and strong compatibility with planar optical, electronic, and microfluidic devices. While attention in the porous materials community has shifted toward macroporous or disordered nanoporous systems, a resurgence in mesoporous thin film research is underway, driven by new molecular self-assembly methods, advanced materials chemistry, and improved characterization techniques. The integration of high-χN block copolymer design, kinetically persistent micelle templating, and postdeposition processing protocols now allows control over structural parameters such as pore size, wall thickness, porosity, and connectivity. These advances have overcome many of the thermodynamic and processing constraints that previously limited widespread adoption. Rather than serving only as high-surface-area supports, mesoporous thin films are engineered as active interfaces where responsive chemistries and nanoscale confinement act in tandem. Embedding switchable ligands, thermoresponsive polymers, redox mediators, or ion-selective groups directly within the pore walls enables real-time control over transport, optical, and electrochemical properties. These capabilities open up new directions in adaptive coatings, gated membranes, and fast-response biosensors. To further expand their functional scope, mesoporous films are integrated into hierarchical and multicomponent architectures. Techniques such as triblock terpolymer templating, crack-directed assembly, and nanoimprint lithography allow for control over spatial organization on the micron and submicron scale and pore system orientation. This enables programmable anisotropy, enhanced molecular diffusion, and wavelength-selective photonic behavior, essential for next-generation sensing, catalysis, and energy applications. Such structural and functional complexity requires equally sophisticated characterization. Multimodal and in situ techniques can track material dynamics under operational conditions. Recent progress includes extended-range ellipsometric porosimetry (EP) for hierarchical architectures, vacuum EP for interface energetics, time-resolved EP for diffusion kinetics, and correlative AFM-SAXS mapping. The introduction of advanced neutron-based spectroscopies, particularly quasielastic neutron scattering (QENS), promises to provide real-time access to ion transport dynamics and segmental motion under nanoscale confinement, offering a path toward deeper mechanistic understanding of structure-performance correlations in mesoporous systems. This Account reflects the technical advances made and the interdisciplinary collaborations that have shaped our collective vision. The particular dimensions of mesopores enable us to subtly tune interactions at the molecular, interfacial, and mesoscopic levels that permit us to harness nanoconfinement. What emerges is a versatile, modular platform capable of chemical gating, energy transduction, and sensing with a level of tunability unmatched by other porous materials. We highlight critical challenges including the need for more robust large-area processing, a deeper understanding of dynamic behavior under cycling, and better integration with device-level architectures. Our strategies support the transition of mesoporous thin films into active high-performance components in next-generation energy, environmental, and biomedical systems.

oxides↗
DOE OSTI2024

The Reservoir Temperature Estimator (RTEst): A multicomponent geothermometry tool

The Reservoir Temperature Estimator (RTEst) is a multicomponent geothermometry tool for estimating reservoir geochemical parameters including reservoir temperature, CO 2 fugacity, mass of water lost or gained, and a reaction factor. It estimates these parameters and their associated uncertainties by minimizing an objective function that is the weighted sum of squares of the saturation indices of a user-selected set of minerals believed to be equilibrated with the reservoir fluid. RTEst accomplishes these estimates by combining the geochemical modeling capabilities of The Geochemist’s Workbench® with the optimization/parameter estimation resources of PEST®. An included interface aids the user in selecting plausible mineral phases to comprise the objective function and calculates their weighting factors. The working principles of RTEst are described and its efficacy is illustrated by presenting results of its application to various geothermal fields with known conditions. These examples show RTEst can account for the alteration of ascending reservoir fluid by mineral (calcite) re-equilibration with changes in temperature, reconstruct waters with CO 2 loss, correct for the deficit of water and other volatiles (CO 2,gas , H 2 S gas ) from boiling, and determine the amount of mixing of thermal and non-thermal waters. RTEst can use data with basis species below detection limit, missing, or unreliable either by assuming equilibrium with a controlling mineral (fixed-analyte method) or by treating the analyte concentration as an optimization parameter. The inverse of variance weighting method included in RTEst provides more representative results than either the normalization or unit weighting methods. Finally, the ability of RTEst to calculate reservoir temperatures, gas fugacity, and mixing fractions demonstrates its usefulness as a tool for evaluating geothermal systems.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Atomistic‐Level Effects of Noncovalent Interactions and Crystalline Packing for Organic Material Structural Integrity upon Exposure to Gamma Radiation

Developing an atomistic understanding of ionizing radiation induced changes to organic materials is necessary for intentional design of greener and more sustainable materials for radiation shielding and detection. Cocrystals are promising for these purposes, but a detailed understanding of how the specific intermolecular interactions within the lattice upon exposure to radiation affect the structural stability of the organic crystalline material is unknown. This study evaluates atomistic-level effects of γ radiation on both single- and multicomponent organic crystalline materials and how specific noncovalent interactions and packing within the crystalline lattice enhance structural stability. Dose studies were performed on all crystalline systems and evaluated via experimental and computational methods. Changes in crystallinity were evaluated by p-XRD and free radical formation was analyzed via EPR spectroscopy. Type of intermolecular interactions and packing within the crystal lattice was delineated and related to the specific free radical species formed and the structural integrity of each material. Periodic DFT and HOMO-LUMO surface mapping calculations provided atomistic-level identifications of the most probable sites for the radicals formed upon exposure to γ radiation and relate intermolecular interactions and molecular packing within the crystalline lattice to experimental results.

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

Functionalized benzylamines from commercial kraft lignin

Benzylamines are key intermediates in pharmaceuticals, agrochemicals, and polymers, but their conventional production relies on benzyl chloride - a petroleum-derived compound with high toxicity and energy demands. Lignin, accounting for up to 30% of plant biomass, is the largest renewable source of aromatic carbon on Earth. However, its highly complex and recalcitrant structure poses a major barrier to efficient conversion into high-value chemicals. Here, in this study, we developed a catalytic approach to convert commercial kraft lignin into phenolic benzylamines through selective depolymerization and subsequent functionalization. We systematically evaluated the effects of three alcohol solvents, formic acid (FA), and a ruthenium-on‑carbon (Ru/C) catalyst on monophenol yield and selectivity. Up to 6.5 wt% monophenol yield was achieved using methanol (MeOH), FA, and Ru/C at 300 °C for 2 h. Quantum thermodynamic simulations based on the COSMO-RS model confirmed the superior solvation and reactivity of the MeOH + FA system, rationalizing observed product yield. The purified monophenolic products, primarily guaiacol and alkyl guaiacols, were then converted into functionalized benzylamines with >90% yield via a multicomponent Mannich reaction under mild conditions. Techno economic analysis (TEA) and life cycle assessment (LCA) underscore the importance of improving lignin depolymerization yields and expanding biorefinery scale. Solvent-only configurations outperform other options in both cost and emissions, with the methanol-only case performing the best ($\$$105 /kg and 26 kg CO 2e /kg) at a large-scale facility. This study establishes a scalable, bio-based pathway for producing benzylamines from commercial kraft lignin, advancing lignin valorization and offering a sustainable alternative to produce petrochemical-based benzylamines.

Catalytic reduction↗
DOE OSTI2023

Large-scale computational polymer solubility predictions and applications to dissolution-based plastic recycling

Dissolution-based plastic recycling is a promising approach to separate and recover high quality pure polymer resins from multicomponent plastic waste by exploiting differences in polymer solubility. The design of a dissolution-based polymer recycling process requires the selection of appropriate solvent systems and operating temperatures to dissolve only target polymers. Determining these parameters experimentally is challenging due to the wide range of solvents and temperatures possible for a given set of target polymers. In this work, we report a computational scheme that employs molecular dynamics simulations and the Conductor-like Screening Model for Realistic Solvents to predict polymer solubilities. Using this scheme, we established a computational solubility database for 8 common polymers and 1007 solvents at multiple temperatures and measured selected solubilities experimentally to validate computational predictions. Analysis of functional groups within this large database then provides chemical heuristics relating the molecular structures of good and non-solvents for selected polymers. We further developed a tool that automates the selection of solvents for all possible sequences in which target polymers can be selectively dissolved to guide the design of dissolution-based plastic recycling processes. Here, we demonstrate the application of these methods via multiple experimental case studies of representative dissolution-based polymer recycling processes in which pure polymer resins were successfully recovered from physical mixtures of polymers.

42 ENGINEERING↗
DOE OSTI2025

Thermodynamic Modeling of Complex Solid Solutions in the Lu-H-N System via Graph Neural Network Accelerated Monte Carlo Simulations

Metal hydrides are important across diverse applications, such as hydrogen storage, batteries, gas sensors, nuclear reactions, and high-temperature superconductivity. Previous computational studies of metal hydrides under extreme pressures, e.g., 𝑂⁡(10 2 ) ⁢GPa, usually treat them as stoichiometric compounds without considering interstitial lattice disorder. As pressures become more moderate in the 𝑂⁡(10 0 ) ⁢GPa and below range, hydrogen disorder at interstitial lattice sites becomes prominent, e.g., in the N-doped Lu hydride that was recently claimed superconducting near 1 GPa. Further adding compositional complexity from alloying and/or multielement interstitial occupation makes elucidating pressure- and temperature-dependent observables intractable by first-principles calculations alone. We therefore propose a lattice graph neural-network surrogate modeling approach to predict configuration- and pressure-dependent equation-of-state properties. Their efficiency permits Monte Carlo simulations to calculate Gibbs energies and pressure-dependent phase diagrams, thereby revealing insights into the synthesis conditions required for achieving desired phase equilibria. We demonstrate this concept for the compositionally complex cubic Lu(H,N,Va) 3 system where three constituents (hydrogen, nitrogen and vacancy) have disordered multielement interstitial occupancies and insights into pressure-dependent phase equilibria are critically needed, e.g., N-doping levels can significantly lower dehydrogenation temperatures and provide a new strategy to optimize hydrogen-storage alloys. This work can improve the thermodynamic understanding of the Lu-H-N system and help rational synthesis of N-doped Lu hydrides, but more generally demonstrates an efficient approach to model pressure-dependent thermodynamics of multicomponent solid solutions.

Monte Carlo methods↗