Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermodynamic assessment”

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

Assessing cathode property prediction via exchange-correlation functionals with and without long-range dispersion corrections

In this work, we benchmark calculated interlayer spacings, average topotactic voltages, thermodynamic stabilities, and band gaps in layered lithium transition-metal oxides (TMOs) and their de-lithiated counterparts, which are used in lithium-ion batteries as positive electrode materials, against available experimental data. Specifically, we examine the accuracy of properties calculated within density functional theory (DFT) using eight different treatments of electron exchange-correlation: the strongly constrained and appropriately normed (SCAN) and Perdew–Burke–Ernzerhof (PBE) density functionals, Hubbard-U-corrected SCAN and PBE (i.e., SCAN+U and PBE+U), and SCAN(+U) and PBE(+U) with added long-range dispersion (D) interactions (i.e., DFT(+U)+D). van der Waals interactions are included respectively via the revised Vydrov-Van Voorhis (rVV10) for SCAN(+U) and the DFT-D3 for PBE(+U). We find that SCAN-based functionals predict larger voltages due to an underestimation of stability of the MO2 systems, while also predicting smaller interlayer spacings compared to their PBE-based counterparts. Furthermore, adding dispersion corrections to PBE has a greater effect on voltage predictions and interlayer spacings than with SCAN, indicating that DFT-SCAN – despite being a ground-state theory – fortuitously captures some short and medium-range dispersion interactions better than PBE. While SCAN-based and PBE-based functionals yield qualitatively similar band gap predictions, there is no significant quantitative improvement of SCAN-based functionals over the corresponding PBE-based versions. Finally, we expect SCAN-based functionals to yield more accurate property predictions than the respective PBE-based functionals for most TMOs, given SCAN's stronger theoretical underpinning and better predictions of systematic trends in interlayer spacings, intercalation voltages, and band gaps obtained in this work.

25 ENERGY STORAGE↗
DOE OSTI2025

Fundamental Insights into Cathode Stability: Linking Compositional Tuning and Local Coordination in Complex Metal Oxides under Aqueous Transformations

Compositional tuning of complex metal oxides in Li-ion battery materials influences their performance as well as their end-of-life behavior, in particular, the tendency to release toxic metal cations in aqueous solution. We modeled ternary variants of a parent LiCoO 2 delafossite structure by varying the metal identity and relative amounts. This yielded ten model formulations of Li(A 4/6 B 1/6 C 1/6 )O 2 , where the material is enriched with the A metal and doped with B and C, with Ni, Mn, Co, Fe, Al, V, and Ti as constituent metals. To assess their stability in aqueous conditions, metal release energetics were calculated using a combination of Density Functional Theory calculations and thermodynamics. Metal release in ternary oxides is dictated by subtle variations in the coordination environment of the leaving group. To identify governing chemical features across diverse compositions with varying local coordination environments, we leverage random forest regression and descriptor importance analysis. A key result is that metal–oxygen orbital hybridization, quantified using a projected density-of-states-derived descriptor, H d/p , provides a physically grounded measure of interaction strength that governs metal release energetics. This refined perspective goes beyond conventional oxidation state considerations and offers more robust insights for materials science. Finally, we model defect surface-bound O 2 dimer formation as a proxy for reactive oxygen species (ROS) generation. The results show that Ni-rich compositions more readily stabilize spin-polarized O 2 dimers, corroborating experimental reports of an increased ROS-driven biological response. In conclusion, our results establish a compositional and electronic basis for metal release and surface oxygen reactivity that form a rationale for complex metal oxide design principles.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Computational discovery of a novel double transition metal nitride MXene and its applications as an anchoring and catalytic material in Li–Se batteries

Double transition metal (DTM) nitride MXenes offer enhanced electrical conductivity and tunable catalytic properties compared to conventional carbide-based MXenes. In this study, we employed first-principles density functional theory (DFT) calculations to discover and computationally validate a novel DTM nitride MXene, Nb 2 TiN 2 , derived from its MAX phase precursor and investigated its potential as an anchoring material (AM) for Li–Se batteries. This newly proposed MXene expands the compositional landscape of DTM nitrides and opens new avenues for functional material design. We performed a comprehensive analysis of the thermodynamic and electronic properties of Nb 2 TiAlN 2 , and the MAX phase precursor to Nb 2 TiN 2 to assess its structural stability and exfoliation potential. Exfoliation energy calculations confirmed the feasibility of synthesizing Nb 2 TiN 2 from Nb 2 TiAlN 2 . We then explored the functionalized form, Nb 2 TiN 2 S 2 , evaluating its capability to serve as an effective anchoring material (AM) in Li–Se batteries by analyzing the reaction mechanisms and kinetics of the selenium reduction reaction (SeRR). Our results indicate that Nb 2 TiN 2 S 2 exhibits a strong binding affinity for lithium polyselenides (Li 2 Se n ), effectively suppressing the shuttle effect. Gibbs free energy calculations for the rate-limiting step of the SeRR reveal favorable kinetics and reduced reaction barriers. Overall, this study provides a detailed evaluation of the structural and electronic properties of a newly proposed DTM nitride MXene and its S-functionalized derivative and the catalyzing effect of Nb 2 TiN 2 S 2 in accelerating the reaction kinetics in Li–Se batteries. These findings underscore the potential importance of the further exploration of MXenes to address current challenges in high-performance Li–Se batteries.

25 ENERGY STORAGE↗
DOE OSTI2023

Can we observe the ion-neutral drift velocity in prestellar cores?

Given the low-ionization fraction of molecular clouds, ambipolar diffusion is thought to be an integral process in star formation. However, chemical and radiative-transfer effects, observational challenges, and the fact that the ion-neutral drift velocity is inherently very small render a definite detection of ambipolar diffusion extremely non-trivial. Here, we study the ion-neutral drift velocity in a suite of chemodynamical, non-ideal magnetohydrodynamic (MHD), two-dimensional axisymmetric simulations of prestellar cores where we alter the temperature, cosmic-ray ionization rate, visual extinction, mass-to-flux ratio, and chemical evolution. Subsequently, we perform a number of non-local thermodynamic equilibrium (non-LTE) radiative-transfer calculations considering various idealized and non-idealized scenarios in order to assess which factor (chemistry, radiative transfer, and/or observational difficulties) is the most challenging to overcome in our efforts to detect the ion-neutral drift velocity. We find that temperature has a significant effect in the amplitude of the drift velocity with the coldest modelled cores (T = 6 K) exhibiting drift velocities comparable to the sound speed. Against expectations, we find that in idealized scenarios (where two species are perfectly chemically co-evolving) the drift velocity ‘survives’ radiative-transfer effects and can in principle be observed. However, we find that observational challenges and chemical effects can significantly hinder our view of the ion-neutral drift velocity. Finally, we propose that HCN and HCNH + ⁠, being chemically co-evolving, could be used in future observational studies aiming to measure the ion-neutral drift velocity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

CALPHAD Uncertainty Quantification and TDBX

CALPHAD uncertainty quantification (UQ) is the foundation of materials design with quantified confidence. We report a framework and software packages to enable CALPHAD UQ assessment and calculation using commercial CALPHAD software (Thermo-Calc). This Bayesian inference framework is coupled with a Markov chain Monte Carlo algorithm to establish uncertainty traces with a given thermodynamic database file (TDB) and corresponding experimental data points. This general framework is demonstrated with the Ni–Cr binary system. The algorithm is firstly validated on synthetic data with known ground truth. Then it is applied to real experimental data to generate posterior traces. We develop a file format named TDBX, which provides a single source of truth by combining the original TDB content and the traces for each assessed Gibbs energy parameter. CALPHAD UQ calculations are performed based on the TDBX file, from which uncertainties for phase boundaries, enthalpy curves, and solidification range are collected as examples of basic design parameters. This TDBX file with corresponding scripts are made open-source. Finally, the combination of CALPHAD UQ assessments and calculations connected by TDBX supports uncertainty-assisted modeling, enabling the integrated application of modern design with uncertainty methodologies to computational materials design.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A unified finite strain gradient-enhanced micropolar continuum approach for modeling quasi-brittle failure of cohesive-frictional materials

Here, in this work, a novel framework for modeling quasi-brittle crack propagation and shear band dominated failure of cohesive-frictional materials like concrete, mortar, rock, tough ceramics, energetic materials, but also granular materials like sands or powders in terms of a unified continuum approach is proposed. It is based on a combination of the gradient-enhanced continuum with gradients of internal variables for representing quasi-brittle cracking, and the micropolar continuum, accounting for the deformation of the microstructure. For developing the gradient-enhanced micropolar framework, the set of balance equations and the kinematic relations are derived, and the constitutive relations are established in a general manner. The framework is formulated in a geometrically exact setting, based on the thermodynamically sound theory of hyperelasto-plasticity, and the numerical implementation by means of the finite element method is discussed. For assessing the approach, realizations of this new approach in terms of constitutive models for particular materials are developed. They are applied to numerical benchmark examples, investigating various loading conditions, and the obtained results are validated by means of a comparison with experiments from the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dynamic Barriers to Crystallization of Calcium Barium Carbonates

Metastable carbonates play important roles in geochemistry, biomineralization and serve as model systems for nonclassical theories of nucleation and growth. Balcite (Ca 0.5 Ba 0.5 CO 3 ) is a remarkable carbonate phase that is isostructural with a high-temperature modification of calcite (CaCO 3 ), yet can be synthesized at ambient conditions. Here, we investigate crystallization pathways in the Ba-Ca-CO 3 -H 2 O system, with a focus on the transformation of amorphous calcium barium carbonate (ACBC) to balcite and subsequent decomposition into the equilibrium calcite (CaCO3) and witherite (BaCO 3 ) phases. Density functional theory calculations show that balcite is an unstable solid solution (Ca 1-x BaxCO 3 , R 3¯ m) in the range 0.17 < x < 0.5, but is accessible through the amorphous ACBC precursor for x ≲ 0.5, and predict its decomposition into calcite and witherite. We confirm this pathway experimentally but found demixing to proceed slowly and remain incomplete even after 9 months. Nucleation kinetics of balcite from ACBC was assessed using a microfluidic assay, where increasing barium content led to a surprising increase in the balcite nucleation rate, despite decreasing thermodynamic driving force. We attribute crystallization rates that dramatically accelerate with time to changes in interfacial structure and composition during coarsening of the amorphous precipitate. By carefully quantifying the thermodynamic and kinetic contributions in the multistep crystallization of a metastable carbonate, we produce insights that allow us to better interpret the formation and persistence of metastable minerals in natural and synthetic environments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Timescales and Processes of Methane Hydrate Formation and Breakdown, With Application to Geologic Systems

Abstract Gas hydrate is an ice‐like form of water and low molecular weight gas stable at temperatures of roughly −10°C to 25°C and pressures of ~3 to 30 MPa in geologic systems. Natural gas hydrates sequester an estimated one sixth of Earth's methane and are found primarily in deepwater marine sediments on continental margins, but also in permafrost areas and under continental ice sheets. When gas hydrate is removed from its stability field, its breakdown has implications for the global carbon cycle, ocean chemistry, marine geohazards, and interactions between the geosphere and the ocean‐atmosphere system. Gas hydrate breakdown can also be artificially driven as a component of studies assessing the resource potential of these deposits. Furthermore, geologic processes and perturbations to the ocean‐atmosphere system (e.g., warming temperatures) can cause not only dissociation, but also more widespread dissolution of hydrate or even formation of new hydrate in reservoirs. Linkages between gas hydrate and disparate aspects of Earth's near‐surface physical, chemical, and biological systems render an assessment of the rates and processes affecting the persistence of gas hydrate an appropriate Centennial Grand Challenge. This paper reviews the thermodynamic controls on methane hydrate stability and then describes the relative importance of kinetic, mass transfer, and heat transfer processes in the formation and breakdown (dissociation and dissolution) of gas hydrate. Results from numerical modeling, laboratory, and some field studies are used to summarize the rates of hydrate formation and breakdown, followed by an extensive treatment of hydrate dynamics in marine and cryospheric gas hydrate systems.

58 GEOSCIENCES↗
DOE OSTI2024

Porous Organic Cage Membranes for Molecular Gas Separations (Final Technical Report for DE-SC0021357)

This proposal aims at demonstrating the development of a novel family of membranes, composed of porous organic cages (POC) which offer the possibility of displaying high separation performance for challenging molecular gas separations relevant to natural gas purification, and olefin/paraffin separation. The proposed POCs synthesized in membrane form will display the most desirable properties of polymers (facile processability and flexibility) and inorganic materials (hierarchically ordered pores with molecular sieving properties) leading to highly selective and permeable membranes. POCs should display distinctive structural, compositional, adsorption and transport properties than those of conventional porous materials, opening the doors for a new research direction in membrane science, and gas separations. Our preliminary results demonstrate the feasibility of preparing POC crystals with controlled size, and continuous POC membranes with remarkable high permeances, and separation ability for CO 2 /CH 4 , N 2 /CH 4 and C 3 H 6 /C 3 H 8 separations serving as a solid foundation for our proposed work. Fundamentally, this proposal aims at elucidating separation mechanisms of different gas mixtures related to natural gas composition, and olefin/paraffin separation over porous organic cage membranes. The proposed research will result in fundamental understanding of adsorption and transport properties of industrially relevant gas molecules through novel microporous membranes, and may lead to the development of a cost effective membrane technology for natural gas purification, and olefin/paraffin separation surpassing the conventional benchmark technology distillation. Furthermore, we aim at demonstrating selective water transport through POC membranes, which can be positively impactful in numerous industrial applications in which water is present. The ability to fabricate thin, chemically and mechanically stable POC membranes for societal relevant gas separations constitute a new and distinctive direction in membrane science. Our proposed work aims at addressing some of the challenges recognized in the Research Agenda for Transforming Separation Science . Specifically: (a) advancing understanding of complex mixtures on separation performance; (b) exploring thermodynamic and kinetic mechanisms through the elucidation of separation mechanisms, and (c) study potential stability issues of the membranes to be assessed by evaluating the long term membrane stability and performance at various temperatures and pressures. The team is uniquely qualified to execute the proposed work. The PI has solid expertise in the rational molecular engineering design of porous crystalline membranes for molecular gas separations. The PNNL collaborator has extensive experience in the synthesis, characterization, and functional applications of microporous crystals, with particular emphasis on gas adsorption.

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

Design principles for intrinsically circular polymers with tunable properties

This perspective discusses a set of design principles for next-generation kinetically trapped, intrinsically circular polymers (iCPs) that are inherently, selectively, and expediently depolymerizable to their monomer state once their kinetic barriers of deconstruction are overcome, thereby enabling not only the ideal shortest chemical circularity but also tunable performance properties. After describing four elements of the design principles - thermodynamics and kinetics, strategies to overcome trade-offs and unify conflicting properties, predictive modeling, and supply-chain life-cycle assessment and techno-economic analysis, which are illustrated with state-of-the-art examples - it concludes with presenting key challenges and opportunities for sustainable development of iCPs.

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

Interplay between London Dispersion, Hubbard U , and Metastable States for Uranium Compounds

High-throughput computational studies of lanthanide and actinide chemistry with density-functional theory are complicated by the need for Hubbard U corrections, which ensure localization of the f-electrons, but can lead to metastable states. This work presents a systematic investigation of the effects of both Hubbard U value and metastable states on the predicted structural and thermodynamic properties of four uranium compounds central to the field of nuclear fuels: UC, UN, UO 2 , and UCl 3 . We also assess the impact of the exchange-hole dipole moment (XDM) dispersion correction on the computed properties. Overall, the choice of Hubbard U value and inclusion of a dispersion correction cause larger variations in the computed geometric properties than result from metastable states. Further, the weak dependence of structure optimization on metastable states should simplify future high-throughput calculations on actinides. Conversely, addition of the dispersion correction is found to offset the repulsion introduced by the Hubbard U term and provides greatly improved agreement with experiment for both cell volumes and heats of formation. The XDM dispersion correction is largely invariant to the chosen U value, making it a robust dispersion correction for actinide systems.

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

Sequential Design of Experiments for Pilot Testing of Novel Solvent System

The CCSI2 program is supporting a six-month test campaign at the National Carbon Capture Center (NCCC) for evaluation of a novel water-lean solvent. This presentation describes CCSI2’s efforts in process modeling of the solvent system for both coal and natural gas-based flue gas sources and initial uncertainty quantification (UQ) work to estimate parametric uncertainty in key sub-models of interest (e.g., thermodynamics, mass transfer, reaction kinetics). Moreover, perspective is provided on how UQ and sequential design of experiments (SDoE) tools are used to assess the impact of model uncertainty on projected process performance, use this information to optimize data collection during the campaign, and refine process models through data collection. This framework is expected to reduce the overall model uncertainty, and thus risk associated with scale-up as the process moves towards commercialization.

Morgan, Joshua↗
DOE OSTI2026

Advancing Mass Timber Buildings: Novel Methods Improve Thermal Assessment and Material Use

For nearly a century, thermal demand calculations for buildings have relied on simplified models developed to match the technical constraints of their era. The first standards, introduced in Germany and Austria in 1929, established climate zones and material conductivity coefficients that, with only incremental updates, still underpin many current assessments. Yet, methods such as Hot box testing, originally designed for lightweight insulation, continue to be applied for mass timber buildings, overlooking thermodynamic characteristics confer real-world advantages. Recent research at Oak Ridge National Laboratory incorporates updated methodologies, aligned with ASHRAE Standard 55 (ASHRAE, 2023) accounting for factors such as thermal inertia, inner surface temperatures, emissivity, solar gains, and dynamic outdoor conditions. These factors better reflect observed heating and cooling loads and highlight opportunities for efficient use of materials in mass timber construction. This work provides a framework for designing comfortable, resilient, and resource-efficient buildings while aligning with performance expectations in energy codes.

Pickett, Robert [International Mass Timber Allianc↗
DOE OSTI2020

Development and Validation of Physics-Based Sub-Models of High-Pressure Supercritical Fuel Injection at Diesel Conditions Insert

The project sought to develop and validate accurate, physics-based, numerical submodels and implement them in computational fluid dynamics (CFD) software codes to enable accurate predictions of diesel sprays at pressures higher than the current systems. The combined experimental-computational effort has four main objectives: (1) acquire spatially and temporally resolved scalar measurements of fuel-oxidizer mixing in the near field of the fuel jet, (2) develop software modules to calculate thermo-physical properties of diesel surrogates and their mixtures with oxidizer of interest in real-gas operating regime, (3) integrate the real-fluid model into a commercial and an open-source CFD code to simulate fuel-air mixing in high-pressure diesel sprays, and (4) assess the robustness, accuracy, and uncertainty of the integrated CFD solvers in computing diesel sprays. Project accomplishments are consistent with the original goals of the project. Although we acquired spatially and temporally resolved scalar measurements, they were limited to refractive index difference and not to thermodynamic properties such as temperature and equivalence ratio. However, we successfully developed the theoretical background needed to acquire such measurements at diesel conditions and demonstrated the capability of the RSD for scalar measurements in a canonical turbulent jet flow replicating features of the diesel jet. Additional work will be needed to implement these procedures at high-pressure, high-temperature engine-like conditions of the CPFR. Mainly, we did not anticipate the need to find the liquid boundary and its statistical features or the state relationships between refractive index and thermodynamic properties, assuming the widely known Dale-Gladstone equation would suffice. Still, efforts devoted to these tasks were essential to fulfil the goals of the project and serves as the foundation on which many of the conclusions and key contributions of this project are based. Further assessment and application of the developed CFD codes is recommended to fully realize the benefits of this research.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

High-Temperature Corrosion of Chromia-Forming Ni-Based Alloys in CO 2 Containing Impurities

Future technologies require structural materials resistant to environmental degradation in high-temperature CO 2 -rich environments. In this work, we exposed several commercially available Ni-based alloys (230, 263, 282, 617, 625, and 740H) to atmospheric pressure gases intended to simulate the compositions expected in future direct-fired supercritical CO 2 power cycles. The alloys were exposed to 95% CO 2 + 4% H 2 O + 1% O 2 and the same gas containing 0.1% SO 2 at temperatures of 600, 650, 700, 750, and 800 °C for 2500 h. With minor exceptions, chromia scales formed on all alloys at all temperatures in the SO 2 -free gas, yielding parabolic growth rates that followed an Arrhenius temperature dependence. Behavior in the SO 2 -containing gas was more complex. Generally, the alloys performed well at temperatures of 650, 750, and 800 °C. While some alloys further performed relatively well across the whole temperature range, several of the alloys experienced chromia failure resulting in non-protective duplex oxide scales and high oxidation rates, at temperatures of 600 and 700 °C. Deviation from protective behavior was associated with internal sulfide formation and, additionally for the case of 600 °C, external sulfate formation. Extensive carburization accompanied growth of the non-protective duplex oxide scales, which made it more difficult for the alloy to recover after initiation of the sulfur-induced accelerated corrosion process. The thermodynamic and kinetic factors influencing the accelerated corrosion in the presence of sulfur are discussed. The results suggest that caution is required when assessing compatibility of Ni-based alloys for CO 2 -based systems when sulfur-based impurities are expected.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dependence of predicted bulk properties of hexagonal hydroxyapatite on exchange-correlation functional

Hydroxyapatite (HA) is the main component of human bones and teeth. HA has also been widely applied in various technological fields due to its unique properties. Reliable computational simulations for this critical material often require input parameters obtained from the first-principles density functional theory (DFT) calculations with appropriate exchange–correlation functionals. Previous DFT calculations are insufficient to verify the reliabilities of various functionals, particularly as they fail to assess predictions for multiple properties of the material. In this paper, we first select 18 different functionals to calculate geometric, elastic, electronic, and thermodynamic properties of hexagonal HA bulk crystal. We find that the results from optB86b-vdW and optB88-vdW functionals with dispersion corrections have the overall best agreement with available experimental data. Then, we choose optB88-vdW functional, as well as PBE functional without dispersion corrections as a comparison, to perform extensive first-principles DFT phonon calculations under the quasiharmonic approximation. Various thermodynamic properties (including phonon contributions to internal energy, entropy, and Helmholtz free energy) and thermal parameters (including thermal expansions of volume, thermal expansion coefficients, heat capacities, isothermal bulk moduli, etc.) versus temperature are consequently obtained. By comparing these quantities, we report that the results from optB88-vdW functional have significantly better agreement with available experimental data than those from PBE functional although the latter has been widely used in previous DFT calculations for HA-based material systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modeling and Analysis of a Thermophotovoltaic Integrated Self-Powered Furnace

This work investigates the energy efficiency and carbon reduction potential of self-powered residential building heating equipment using a thermodynamic modeling approach. An integrated thermophotovoltaic power module and residential scale furnace system (40,000 Btu/h) were modeled and studied in detail to assess the influence of different design configurations on primary energy efficiency. Operational characteristics such as total power generation, electrical efficiency, and heat recovery were examined in a self-powered system configuration. A sensitivity analysis was conducted to determine the influence of the electric grid’s carbon dioxide footprint (carbon intensity) and the cost of electricity on the environmental, as well as the economic, benefit associated with the self-powered configuration. Compared with a traditional furnace powered by an electric grid at a carbon intensity of 0.5 kg CO 2eq /kWh EL , the self-powered furnace was shown to decrease the annual carbon dioxide emissions by approximately 550 kg (~75% reduction), while also saving more than USD 200 in utility expenses, annually. Additionally, the carbon emission reduction potential of blending different concentrations of hydrogen in natural gas fuel was also studied.

14 SOLAR ENERGY↗
DOE OSTI2021

Combining artificial intelligence and physics-based modeling to directly assess atomic site stabilities: from sub-nanometer clusters to extended surfaces

The performance of functional materials is dictated by chemical and structural properties of individual atomic sites. In catalysts, for instance, the thermodynamic stability of constituting atomic sites is a key descriptor from which more complex properties, such as molecular adsorption energies and reaction rates, can be derived. In this study, we present a widely applicable machine learning (ML) approach to instantaneously compute the stability of individual atomic sites in structurally and electronically complex nano-materials. Conventionally, we determine such site stabilities using computationally intensive first-principles calculations. With our approach, we predict the stability of atomic sites in sub-nanometer metal clusters of 3–55 atoms with mean absolute errors in the range of 0.11–0.14 eV. To extract physical insights from the ML model, we introduce a genetic algorithm (GA) for feature selection. This algorithm distills the key structural and chemical properties governing the stability of atomic sites in size-selected nanoparticles, allowing for physical interpretability of the models and revealing structure–property relationships. The results of the GA are generally model and materials specific. In the limit of large nanoparticles, the GA identifies features consistent with physics-based models for metal–metal interactions. By combining the ML model with the physics-based model, we predict atomic site stabilities in real time for structures ranging from sub-nanometer metal clusters (3–55 atom) to larger nanoparticles (147 to 309 atoms) to extended surfaces using a physically interpretable framework. Finally, we present a proof of principle showcasing how our approach can determine stable and active nanocatalysts across a generic materials space of structure and composition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗