Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Thermochemistry”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Regularized machine learning on molecular graph model explains systematic error in DFT enthalpies

Abstract A major goal of materials research is the discovery of novel and efficient heterogeneous catalysts for various chemical processes. In such studies, the candidate catalyst material is modeled using tens to thousands of chemical species and elementary reactions. Density Functional Theory (DFT) is widely used to calculate the thermochemistry of these species which might be surface species or gas-phase molecules. The use of an approximate exchange correlation functional in the DFT framework introduces an important source of error in such models. This is especially true in the calculation of gas phase molecules whose thermochemistry is calculated using the same planewave basis set as the rest of the surface mechanism. Unfortunately, the nature and magnitude of these errors is unknown for most practical molecules. Here, we investigate the error in the enthalpy of formation for 1676 gaseous species using two different DFT levels of theory and the ‘ground truth values’ obtained from the NIST database. We featurize molecules using graph theory. We use a regularized algorithm to discover a sparse model of the error and identify important molecular fragments that drive this error. The model is robust to rigorous statistical tests and is used to correct DFT thermochemistry, achieving more than an order of magnitude improvement.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Application of NEAMS Multiphysics Framework for Species Tracking in Molten Salt Reactors

This report from Idaho National Laboratory (INL) summarizes the key modeling and simulation activities conducted under the Department of Energy (DOE) Molten Salt Reactor (MSR) Campaign during the Fiscal Year 2023 (FY23). The focus of the work was to leverage state-of-the-art modeling capabilities from the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) codes to enable novel multiphysics and multiscale modeling and simulation of MSRs. Through collaboration with NEAMS code developers, advanced multiphysics analysis capabilities for MSR systems were demonstrated by coupling depletion, thermal-hydraulics, and thermochemistry into an innovative framework for chemical species transport in MSRs. As a result, the framework can track nuclides throughout their lifetimes in the core, from production (depletion) to advection throughout the salt volume (thermal-hydraulics) and off-gassing or precipitation outside of the salt (thermochemistry). This work supports the near-term deployment of MSRs by integrating the synergistic efforts between the DOE’s MSR Campaign and NEAMS program. The resulting framework will help better connect system design modelers with experimentalists to better understand and predict complex physical behaviors in MSRs. Researchers and MSR developers alike can now leverage these new modeling and simulation capabilities to perform novel analyses with applications including: • MSR dynamics during normal operational transients and accident scenarios • Off-gas system design and performance for fuel cycle and depletion analysis • Corrosion and active chemistry control for reactor component health and lifetime determination • Source term, decay heat and activity determination in accident scenarios • Special nuclear material accountancy and chemical forensic analysis for safeguards • Digital twin development of experiments and experimental reactor demonstrations • Measurement requirements for instrumentation and control design • Uncertainty and sensitivity analysis of missing data to inform future experimental data collection.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Free Energies of Proton-Coupled Electron Transfer Reagents and Their Applications

Here, we present an update and revision to our 2010 review on the topic of proton- coupled electron transfer (PCET) reagent thermochemistry. Over the past decade, the data and thermochemical formalisms presented in that review have been of value to multiple fields. Concurrently, there have been advances in the thermochemical cycles and experimental methods used to measure these values. This Review (i) summarizes those advancements, (ii) corrects systematic errors in our prior review that shifted many of the absolute values in the tabulated data, (iii) provides updated tables of thermochemical values, and (iv) discusses new conclusions and opportunities from the assembled data and associated techniques. We advocate for updated thermochemical cycles that provide greater clarity and reduce experimental barriers to the calculation and measurement of Gibbs free energies for the conversion of X to XHn in PCET reactions. In particular, we demonstrate the utility and generality of reporting potentials of hydrogenation, E°(V vs H 2 ), in almost any solvent and how these values are connected to more widely reported bond dissociation free energies (BDFEs). The tabulated data demonstrate that E°(V vs H 2 ) and BDFEs are generally insensitive to the nature of the solvent and, in some cases, even to the phase (gas versus solution). This Review also presents introductions to several emerging fields in PCET thermochemistry to give readers windows into the diversity of research being performed. Some of the next frontiers in this rapidly growing field are coordination-induced bond weakening, PCET in novel solvent environments, and reactions at material interfaces.

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

Active Thermochemical Tables: Enthalpies of Formation of Bromo- and Iodo-Methanes, Ethenes and Ethynes

Here, the thermochemistry of halocarbon species containing iodine and bromine is examined through an extensive interplay between new Feller–Peterson–Dixon (FPD) style composite methods and a detailed analysis of all available experimental and theoretical determinations using the thermochemical network that underlies the Active Thermochemical Tables (ATcT). From the computational viewpoint, a slower convergence of the components of composite thermochemistry methods is observed relative to species that solely contain first row elements, leading to a higher computational expense for achieving comparable levels of accuracy. Potential systematic sources of computational uncertainty are investigated, and, not surprisingly, spin-orbit coupling is found to be a critical component, particularly for iodine containing molecular species. The ATcT analysis of available experimental and theoretical determinations indicates that prior theoretical determinations have significantly larger uncertainties than originally reported, particularly in cases where molecular spin-orbit effects were ignored. Accurate and reliable heats of formation are reported for 38 halogen containing systems, based on combining the current computations with previous experimental and theoretical work via the ATcT approach.

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

“Best” Iterative Coupled-Cluster Triples Model? More Evidence for 3CC

To follow up on the unexpectedly good performance of several coupled-cluster models with approximate inclusion of 3-body clusters we performed a more complete assessment of the 3CC method for accurate computational thermochemistry in the standard HEAT framework. New spin-integrated implementation of the 3CC method applicable to closed- and open-shell systems utilizes a new automated toolchain for derivation, optimization, and evaluation of operator algebra in many-body electronic structure. We found that with a double-ζ basis set the 3CC correlation energies and their atomization energy contributions are almost always more accurate (with respect to the CCSDTQ reference) than the CCSDT model as well as the standard CCSD(T) model. The mean absolute errors in cc-pVDZ {3CC, CCSDT, and CCSD(T)} electronic (per valence electron) and atomization energies relative to the CCSDTQ reference for the HEAT data set, were {24, 70, 122} μE h /e and {0.46, 2.00, 2.58} kJ/mol, respectively. The mean absolute errors in the complete-basis-set limit {3CC, CCSDT, and CCSD(T)} atomization energies relative to the HEAT model reference, were {0.52, 2.00, and 1.07} kJ/mol, The significant and systematic reduction of the error by the 3CC method and its lower cost than CCSDT suggests it as a viable candidate for post- CCSD(T) thermochemistry applications, as well as the preferred alternative to CCSDT in general.

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

Elaborated thermochemical treatment of HF, CO, N 2 , and H 2 O: Insight into HEAT and its extensions

Here, empirical, highly accurate non-relativistic electronic total atomization energies (eTAEs) are established by combining experimental or computationally converged treatments of the nuclear motion and relativistic contributions with the total atomization energies of HF, CO, N 2 , and H 2 O obtained from the Active Thermochemical Tables. These eTAEs, which have estimated (2σ) uncertainties of less than 10 cm -1 (0.12 kJ mol -1 ), form the basis for an analysis of high-level ab initio quantum chemical calculations that aim at reproducing these eTAEs for the title molecules. The results are then employed to analyze the performance of the high-accuracy extrapolated ab initio thermochemistry, or High-Accuracy Extrapolated Ab Initio Thermochemistry (HEAT), family of theoretical methods. The method known as HEAT-345(Q), in particular, is found to benefit from fortuitous error cancellation between its treatment of the zero-point energy, extrapolation errors in the Hartree-Fock and coupled cluster contributions, neglect of post-(T) core-correlation, and the basis-set error involved in higher-level correlation corrections. In addition to shedding light on a longstanding curiosity of the HEAT protocol—where the cheapest HEAT-345(Q) performs comparably to the theoretically more complete HEAT-456QP procedure—this study lays the foundation for extended HEAT variants that offer substantial improvements in accuracy relative to the established approaches.

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

Thermochemical Modeling in Molten Fluoride Salts for Radionuclide Speciation

An important aspect of the licensing process for nuclear reactors is providing a reasonable assurance of safety to the general public. This includes modeling potential radionuclide releases from the reactor during normal operations and accident scenarios, which is known as the reactor’s source term. A new class of advanced (non-LWR) reactors are being developed which use molten salts as the coolant fluid. Because the molten salt coolant represents a credited barrier for radionuclide transport between the fuel and the environment, a necessary aspect of mechanistic source term (MST) analysis for the KP-FHR is modeling the thermochemistry of molten salts. Provided here is a review of the theory of the thermodynamic principles governing multicomponent phase equilibria, the background of molten salt thermochemistry research, and a summary of the thermochemical data relevant to the KP-FHR coolant salt, Li 2 BeF 4 , commonly referred to as “FLiBe”. A review of literature is followed by a brief introduction to methods that can be used to model the thermochemical behavior of molten salt mixtures. The methodology outlined is based on the use of a commercial thermodynamic modeling software called FactSage, which is one of only a few available softwares based on the modified quasichemical model(MQM), which is the recommended solution model for molten salts.

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

Systematically derived thermodynamic properties for alkane oxidation

Key combustion properties, such as ignition delay time, show strong sensitivity to the thermochemistry of the main species in the standard radical oxidation pathway at low temperatures (600-1000K). Significant uncertainties persist in current estimates of thermodynamic properties, particularly for the larger species of relevance to common practical fuels. In this article, we use advanced computational schemes to evaluate thermodynamic properties for the fuel, fuel radical, peroxy, hydroperoxy-alkyl radical, and hydroperoxide species for a set of 17 fuels containing up to 9 carbon atoms, and with various degrees of branching in the alkane backbone. The procedure, termed STAR-1D, combines conformer sampling to find the minimum geometry, B2PLYP-D3/cc-pVTZ harmonic frequency evaluations, and omega B97X-D/cc-pVTZ one-dimensional torsional mapping. It includes two physically-based scaling routines: a frequency-dependent scaling to B2PLYP-D3/cc-pVTZ anharmonic frequencies and a scaling of one-dimensional omega B97X-D/ccpVTZ torsional profiles to reproduce the product of the B2PLYP-D3/cc-pVTZ frequencies in the harmonic limit. Substantive comparisons with existing experimental databases, together with careful examinations of key theoretical assumptions, are used to explore the accuracy of the predictions. These computationally intensive explorations of 195 species were facilitated by automated thermochemistry software. Further, high accuracy 0 K heats of formation from a separate study are used with the STAR-1D computations to generate NASA polynomial representations. In addition to their intrinsic value, the present results also provide a reliable database for the optimization of group additivity or machine learning schemes for scaling to larger combustion systems. Towards this end a complementary extensive conformational analysis is carried out for the medium sized species and the thermodynamic properties of the lowest energy hydrogen-bonded and non-hydrogen bonded conformers are contrasted for larger hydroperoxy-alkyl radical species.

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

A workflow for automatic generation and efficient refinement of individual pressure-dependent networks

Manual chemical kinetic model construction often requires modelists to at least implicitly guess all possible decay pathways and their relative fluxes for all important species. We present a workflow and associated tool for automatic generation and re finement of pressure-dependent networks that should enable modelists to efficiently and comprehensively identify decay pathways and estimate respective parameters for chemical species. This tool combines the capabilities of the Reaction Mechanism Generator (RMG) software to generate possible reaction paths, determine thermochemistry, approximate rate coefficients and estimate frequencies with the capabilities of the Arkane software to use quantum chemical parameters to compute thermochemistry and pressure dependent rates. A flux-based algorithm is used to decide which isomers to add to the network. Isomers added to the network are reacted to form other channels. When enough of the flux is accounted for by the isomers and bimolecular product channels, the generation process terminates to yield a comprehensive network. Network sensitivity analysis is applied to the network to identify the most important wells and barriers that can be refined using quantum chemistry calculations. Iterative refinement can then be used to achieve accuracy. The comprehensive network can be easily reduced using energy and flux-based algorithms to the most important channels and isomers in the network.

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

The C 2 H 4 O isomers in the oxidation of ethylene

A combined rotational spectroscopy, thermochemistry, and kinetic modeling study explores the mechanism of acetaldehyde (ethanal), vinyl alcohol (ethenol), and ethylene oxide (oxirane) formation in the oxidation of ethylene (ethene). Multiplexed quantitative detection of oxidation intermediates and products exiting a SiO 2 / SiC microreactor at 1700 K is demonstrated with the BrightSpec W-band chirped-pulse Fourier transform millimeter-wave spectrometer. The broadband rotational spectrum contains transitions of formaldehyde (CH 2 O), methoxy (CH 3 O), methanol (CH 3 OH), ketene (CH 2 CO), acetaldehyde (CH 3 CHO), syn-vinyl alcohol (syn- CH 2 CHOH), anti-vinyl alcohol (anti-CH 2 CHOH), oxirane (c-CH 2 OCH 2 ), propyne (CH 3 CCH), and syn-propanal (syn-CH 3 CH 2 CHO). We focus on the three C 2 H 4 O species and deduce their concentration ratio [CH 3 CHO]:[CH 2 CHOH]:[c-CH 2 OCH 2 ] = 1:0.7(2):0.06(2) by comparing the observed line intensities to those simulated with the PGOPHER software. Detailed thermochemistry of the C 2 H 4 O isomers and conformers is provided by Active Thermochemical Tables. The observed excess concentrations of vinyl alcohol and oxirane relative to the more stable acetaldehyde compared to the equilibrium ratio at 1700 K (1:0.087:0.000024), point to direct chemical pathways to these higher energy isomers. The mechanism for the formation of the three C 2 H 4 O isomers is analyzed using a 0-D homogenous reactor kinetics simulation for ethylene oxidation. The ratios of the C 2 H 4 O isomers concentrations predicted by the kinetic model are compared to the experimental values.

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

Modeling and simulation of oxygen transport in high burnup LWR fuel

We have developed a formulation for oxygen transport in uranium dioxide nuclear fuel that accounts for the effects of irradiation. The overall simulation combines the evolving isotopic composition, thermochemistry, and oxygen transport in irradiated fuel. The driving forces for oxygen transport are computed from local thermodynamic equilibrium calculations and account for the effects of temperature gradients and composition, including fission products. The proposed method provides a mechanism for including complex thermodynamic models of nuclear fuel in modeling of mass redistribution, and alleviates difficulties associated with the common thermodiffusion formulation. Here, the transport model has been implemented within the nuclear fuel performance code BISON utilizing the thermochemistry code Thermochimica, with burnup calculations provided by the ORIGEN isotopic transmutation code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Developing Scaling Relationships for Molecular Electrocatalysis through Studies of Fe-Porphyrin-Catalyzed O 2 Reduction

The oxygen reduction reaction (ORR) is a multiproton/multielectron transformation in which dioxygen (O 2 ) is reduced to water or hydrogen peroxide and serves as the cathode reaction in most fuel cells. The ORR (O 2 + 4e – + 4H + → 2H 2 O) involves up to nine substrates and thus requires navigating a complicated reaction landscape, typically with several high-energy intermediates. Many catalysts can perform this reaction, though few operate with fast rates and at low overpotentials (close to the thermodynamic potential). Attempts to optimize these parameters, both in homogeneous and heterogeneous electrocatalytic systems, have focused on modifying catalyst design and understanding kinetic/thermodynamic relationships between catalytic intermediates. One such method for analyzing and predicting catalyst reactivity and efficiency has been the development of “molecular scaling relationships”. In this work, we share our experience deriving and utilizing molecular scaling relationships for soluble, iron-porphyrin-catalyzed O 2 reduction in organic solvents. These relationships correlate turnover frequencies (TOF max ) and effective overpotentials (η eff ), properties uniquely defined for homogeneous catalysts. Following a general introduction of scaling relationships for both homogeneous and heterogeneous electrocatalysis, we describe the components of such scaling relationships: (i) the overall thermochemistry of the reaction and (ii) the rate and rate law of the catalyzed reaction. We then show how connecting these thermodynamic and kinetic parameters reveals multiple molecular scaling relationships for iron-porphyrin-catalyzed O 2 reduction. For example, the log(TOF max ) responds steeply to changes in ηeff that result from different catalyst reduction potentials (18.5 decades in TOF max /V in η eff ) but much less dramatically to changes in η eff that arise from varying the pKa of the acid buffer (5.1 decades in TOF max /V in η eff ). Thus, a single scaling relationship is not always sufficient for describing molecular electrocatalysis. This is particularly evident when the catalyst identity and reaction conditions are coupled. Using these multiple scaling relationships, we demonstrate that the metrics of turnover frequency and effective overpotential can be predictably tuned to achieve faster rates at lowered overpotentials. This Account uses a collection of related stories describing our research on soluble iron-porphyrin-catalyzed ORR to show how molecular scaling relationships can be derived and used for any electrocatalytic reaction. Such scaling relationships are powerful tools that connect the thermochemistry, mechanism, and rate law for a catalytic system. We hope that this collection shows the utility and simplicity of the molecular scaling approach for understanding catalysis, for enabling direct comparisons between catalyst systems, and for optimizing catalytic processes.

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

Uncertainty Quantification and Error Propagation in the Enthalpy and Entropy of Surface Reactions Arising from a Single DFT Functional

Accounting for parametric uncertainty in models is essential for quantifying the models’ predictive ability. Recently, approaches have been introduced to estimate parametric uncertainty in kinetic models while accounting for correlations among energy parameters. However, correlations have been estimated indirectly and correlations in entropies have not been accounted for. For surface-catalyzed microkinetic models of >C2 (more than two carbon-containing) molecules, which consist of thousands of reaction steps and intermediate surface species, first-principles density functional theory (DFT) is costly, and thus, estimation of thermochemistry and reaction barriers requires surrogate methods of DFT, such as group additivity and Brønsted–Evans–Polanyi relationships, respectively. For such parametrization, model uncertainty is unclear. This work develops a framework to overcome these gaps using group additivity and a single DFT functional. We estimate correlations in parameters of kinetic models and quantify uncertainty for thermochemistry, reaction barriers, reaction paths, and ultimately reaction rates, accounting also for the contribution of entropic uncertainty. Furthermore, the approach is illustrated on propane combustion and ethane oxidative dehydrogenation reactions.

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

Proton-Coupled Defects Impact O–H Bond Dissociation Free Energies on Metal Oxide Surfaces

Proton-coupled electron transfer (PCET) reactions on metal oxides require coupling between proton transfer at the solid-liquid interface and electron transfer involving defects at or near the band edge. Herein, hybrid functional periodic density functional theory is used to elucidate the impact of proton-coupled defects on the bond dissociation free energies (BDFEs) of O-H bonds on anatase TiO 2 surfaces. These O-H BDFEs are directly related to interfacial PCET thermochemistry. Comparison between geometrically similar O-H bonds associated with different defect types, namely conduction d-band electrons or valence p-band holes, reveals that the BDFEs differ by ~81 kcal/mol (3.50 eV), comparable to the wide TiO 2 band gap. These differences are shown to be determined primarily by differences in electron transfer driving forces, which are analyzed by using band energies and inner-sphere reorganization energies within a Marcus theory framework. Furthermore, these fundamental insights about the impact of proton-coupled defects on PCET thermochemistry at semiconductor surfaces have broad implications for electrocatalysis.

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

A graph-based network for predicting chemical reaction pathways in solid-state materials synthesis

Abstract Accelerated inorganic synthesis remains a significant challenge in the search for novel, functional materials. Many of the principles which enable “synthesis by design” in synthetic organic chemistry do not exist in solid-state chemistry, despite the availability of extensive computed/experimental thermochemistry data. In this work, we present a chemical reaction network model for solid-state synthesis constructed from available thermochemistry data and devise a computationally tractable approach for suggesting likely reaction pathways via the application of pathfinding algorithms and linear combination of lowest-cost paths in the network. We demonstrate initial success of the network in predicting complex reaction pathways comparable to those reported in the literature for YMnO 3 , Y 2 Mn 2 O 7 , Fe 2 SiS 4 , and YBa 2 Cu 3 O 6.5 . The reaction network presents opportunities for enabling reaction pathway prediction, rapid iteration between experimental/theoretical results, and ultimately, control of the synthesis of solid-state materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Optimizing the regularization in size-consistent second-order Brillouin-Wigner perturbation theory

Despite its simplicity and relatively low computational cost, second-order Møller-Plesset perturbation theory (MP2) is well-known to overbind noncovalent interactions between polarizable monomers and some organometallic bonds. In such situations, the pairwise-additive correlation energy expression in MP2 is inadequate. Although energy-gap dependent amplitude regularization can substantially improve the accuracy of conventional MP2 in these regimes, the same regularization parameter worsens the accuracy for small molecule thermochemistry and density-dependent properties. Recently, we proposed a repartitioning of Brillouin-Wigner perturbation theory that is size-consistent to second order (BW-s2), and a free parameter ($α$) was set to recover the exact dissociation limit of H 2 in a minimal basis set. Alternatively $α$ can be viewed as a regularization parameter, where each value of $α$ represents a valid variant of BW-s2, which we denote as BW-s2($α$). In this work, we semi-empirically optimize $α$ for noncovalent interactions, thermochemistry, alkane conformational energies, electronic response properties, and transition metal datasets, leading to improvements in accuracy relative to the ab initio parameterization of BW-s2 and MP2. We demonstrate that the optimal $α$ parameter ($α$ = 4) is more transferable across chemical problems than energy-gap-dependent regularization parameters. This is attributable to the fact that the BW-s2($α$) regularization strength depends on all of the information encoded in the t amplitudes rather than just orbital energy differences. While the computational scaling of BW-s2($α$) is iterative $\mathcal{O}$($N^5$), this effective and transferable approach to amplitude regularization is a promising route to incorporate higher-order correlation effects at second-order cost.

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

Accelerating Thermochemical Equilibrium Calculations for Nuclear Reactor Applications

Thermochemical properties play a key role in modeling and simulation of several key phenomena in nuclear reactors. There has been an increasing interest in incorporating CALPHAD-based formulations in multiphysics simulations including for Molten Salt Reactors where knowledge of phase evolution of the salt and the chemical potentials of various elements are of utmost importance in source term analyses and redox control. However, the size of such simulations is often limited by the high computational cost of full thermodynamic equilibrium calculations. This work discusses the current efforts aimed at accelerating thermochemical equilibrium calculations for multiphysics simulations performed using the open-source finite element / finite volume code Multiphysics Object Oriented Simulation Environment (MOOSE) [1]. While several methods have been proposed for accelerating phase equilibrium calculations [2], most focus on relatively small systems and often rely on a- priori knowledge of the state-space of the system. Nuclear materials, however, are often multi-component systems owing to the evolution of composition under irradiation and an approach based on a-priori mapping of phase diagram is often not enough. This work is aimed at demonstrating an on-the-fly surrogate modeling framework that uses active learning to reduce the number of full equilibrium calculations that must be performed. By combining with efficient coupling approaches, the surrogate framework helps in reducing the computational cost of thermodynamic equilibrium informed multiphysics simulations of nuclear materials. The performance is benchmarked against full coupling with the thermochemistry library Thermochimica [3]. This work uses a machine learning based approach for constructing surrogate models to predict the stable phases in a multicomponent system. The surrogates were constructed using neural networks and Gaussian process classification. In this work, we compare the relative performance of the two methods. We also demonstrate the use of caching previous calculations by interpolating the values from nearest neighbors. References [1] Lindsay, A.D., et al. "2.0 – MOOSE: Enabling massively parallel multiphysics simulation", SoftwareX, 20 (2022): 101202. [2] Roos, W.A. and Zietsman J.H. "Accelerating complex chemical equilibrium calculations – A Review", Calphad, 77 (2022): 102380. [3] Piro, M.H.A., et al. "The thermochemistry library Thermochimica", Computational Materials Science, 67 (2013): 266-272.

36 MATERIALS SCIENCE↗