Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermodynamic stability”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Effect of nitrogen doping and pressure on the stability of $\mathrm{LuH_3}$

A report on the near-ambient superconductivity in a nitrogen-doped lutetium hydride has stimulated great interest in this material [Dasenbrock-Gammon et al. , Nature (London) 615, 244 (2023)]. While its superconductivity is still a subject of debate, the structure of the claimed cubic phase remains uncertain. In this work, we study the effect of nitrogen doping and pressure on the energetic and dynamic stability of cubic $\mathrm{LuH_3}$. Our findings indicate that both pressure and nitrogen doping can enhance the stability of the cubic $\mathrm{LuH_3}$ phase. Here, we propose a $\mathrm{Lu_8}$$\mathrm{H_{21}}$$\mathrm{Lu}$ structure that exhibits a stable phonon, reasonable thermodynamic stability at 1 GPa, and an x-ray diffraction pattern similar to the experimental data. However, we do not observe electron-phonon coupling in the zone-center phonon modes of these crystal structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Gallium‐Doping Effects on Structure, Lithium‐Conduction, and Thermochemical Stability of Li 7‐3 x Ga x La 3 Zr 2 O 12 Garnet‐Type Electrolytes

Abstract One of the most promising electrolytes for all‐solid‐state lithium batteries is Li 7 La 3 Zr 2 O 12 . Previously, their thermodynamic stability, Li‐ion conductivity, and structural features induced by Ga‐doping have not been empirically determined or correlated. Here, their interplay was examined for Li 7−3 x Ga x La 3 Zr 2 O 12 with target x Ga=0, 0.25, 0.50, 0.75, and 1.00 atoms per formula unit (apfu). Formation enthalpies, obtained with calorimetry and found to be exothermic at all compositions, linearly decreased in stability with increased x Ga. At dilute x Ga substitution, the formation enthalpy curve shifted stepwise endothermically, and the conductivity increased to a maximum, coinciding with 0.529 Ga apfu. This correlated with percolation threshold analysis (0.558 Ga apfu). Further substitution (0.787 Ga apfu) produced a large decrease in the stability and conductivity due to a large increase in point defects and blocked Li‐migration pathways. At x Ga=1.140 apfu, a small exothermic shift was related to defect cluster organization extending the Li hopping distance and decreased Li‐ion conductivity.

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

High-throughput and data-driven search for stable optoelectronic AMSe 3 materials

The rapid advancement in emerging optoelectronic technologies demands highly efficient, affordable, and ecofriendly materials. In this context, ternary chalcogenides, especially ternary selenides, show early promise as a material class due to their stability and remarkable electronic, optical, and transport properties. In this work, we integrate first-principles-based high-throughput computations with machine learning (ML) techniques to predict the thermodynamic stability and optoelectronic properties of 920 valency-satisfied selenide compounds. Through investigating polymorphism, our study reveals the edge-sharing orthorhombic Pnma phase (NH 4 CdCl 3 -type) as the most stable structure for most ternary selenides. High-fidelity supervised ML models are trained and tested to accelerate stability and band gap predictions. These data-driven models pin down the most influential features that dominantly control key material characteristics. The multistep high-throughput computations identify the ternary selenides with optimal direct band gaps, light carrier masses, and strong optical absorption edges. The extensive materials screening considering phase stability, toxicity, and defect tolerance, finally identifies the seven most suitable candidates for photovoltaic applications. Two of these final compounds, SrZrSe 3 and SrHfSe 3 , have already been synthesized in a single-phase form, with the latter showing an optically suitable band gap, aligning well with our findings. The non-adiabatic molecular dynamics reveal sufficiently long photoexcited charge carrier lifetimes (on the order of nanoseconds) in some of these selected selenide materials, indicating their exciting characteristics. Overall, our study suggests a robust in silico framework that can be extended to screen large datasets of various material classes for identifying promising photoactive candidates.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Designing better electrolytes

Electrolytes and the associated interphases constitute the critical components to support the emerging battery chemistries that promise tantalizing energy but involve drastic phase and structure complications. Designing better electrolytes and interphases holds the key to the success of these batteries. As the only component that interfaces with every other component in the device, an electrolyte must satisfy multiple criteria simultaneously. These include transporting ions while insulating electrons between the electrodes and maintaining stability against electrodes of extreme chemical natures: the strongly oxidative cathode and the strongly reductive anode. In most advanced batteries, the two electrodes operate at potentials far beyond the thermodynamic stability limits of electrolytes, so the stability therein has to be realized kinetically through an interphase formed from the sacrificial reactions between electrolyte and electrodes.

Meng, Y. Shirley↗
DOE OSTI2023

Structure–Thermodynamic Relationship of a Polysaccharide Gel (Alginate) as a Function of Water Content and Counterion Type (Na vs Ca)

Biofilms are the predominant mode of microbial life on Earth, and so a deep understanding of microbial communities-and their impacts on environmental processes-requires a firm understanding of biofilm properties. Because of the importance of biofilms to their microbial inhabitants, microbes have evolved different ways of engineering and reconfiguring the matrix of extracellular polymeric substances (EPS) that constitute the main non-living component of biofilms. This ability makes it difficult to distinguish between the biotic and abiotic origins of biofilm properties. An important route toward establishing this distinction has been the study of simplified models of the EPS matrix. This study builds on such efforts by using atomistic simulations to predict the nanoscale (≤10 nm scale) structure of a model EPS matrix and the sensitivity of this structure to interpolymer interactions and water content. To accomplish this, we use replica exchange molecular dynamics (REMD) simulations to generate all-atom configurations of ten 3.4 kDa alginate polymers at a range of water contents and Ca–Na ratios. Simulated systems are solvated with explicitly modeled water molecules, which allows us to capture the discrete structure of the hydrating water and to examine the thermodynamic stability of water in the gels as they are progressively dehydrated. Our primary findings are that (i) the structure of the hydrogels is highly sensitive to the identity of the charge-compensating cations, (ii) the thermodynamics of water within the gels (specific enthalpy and free energy) are, surprisingly, only weakly sensitive to cation identity, and (iii) predictions of the differential enthalpy and free energy of hydration include a short-ranged enthalpic term that promotes hydration and a longer-ranged (presumably entropic) term that promotes dehydration, where short and long ranges refer to distances shorter or longer than ~0.6 nm between alginate strands.

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

Epitaxial Stabilization and Oxygen Evolution Reaction Activity of Metastable Columbite Iridium Oxide

Non-rutile polymorphs of binary iridium oxide such as columbite IrO 2 (α-IrO 2 ) are promising candidates for highly active acid-stable oxygen evolution reaction (OER) catalysts, yet their synthesis has been challenging due to the dominant thermodynamic stability of rutile IrO 2 (R-IrO 2 ). Here, we report the growth of α-IrO 2 via epitaxial thin films using pulsed laser deposition. We observe that, in competition with R-IrO 2 (100), the films can be optimized to be predominantly (100)-oriented α-IrO 2 . Surprisingly, the activity of α-IrO 2 shows a large discrepancy of ~0.2 V in the overpotential compared to its predicted activity, which is resolved via theoretical calculations to be a crystal orientation effect. This demonstrates that the electrocatalytic activity can be significantly varied upon crystal orientation, a parameter that is difficult to control in conventional polycrystalline systems but accessible in epitaxial thin films. In total, this study demonstrates epitaxial thin film growth as a powerful technique, which can overcome large energetic instabilities on the order of ~300 meV to stabilize metastable material structures inaccessible by bulk synthesis. As a result, this provides unique opportunities to effectively identify the atomic structure of active catalysts by combining investigations of metastable materials with theoretical predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The critical role of intrinsic defects and many-body interactions on the stability of MnBi2Te4

Intrinsic antisite defects pose a major challenge to understanding and predicting the exotic properties of the layered topological magnetic insulator MnBi2Te4 (MBT). In this work, we study the origin of the abundance of intrinsic defects in MBT, including many-body defect–defect interactions and many-body electronic correlations. Until now, ab initio methods have struggled to explain thermodynamic stability and properties influenced by defect behavior in MBT. We model native Mn–Bi antisite defects in MBT at finite temperatures using a cluster expansion that includes defect–defect interactions. To overcome the limitations of conventional density functional theory (DFT), we introduce a hybrid approach that incorporates high-accuracy quantum Monte Carlo (QMC) calculations, introducing missing correlations. This strategy allows for accurate estimation of defect energetics and finite-temperature properties. We compute the configurational free energy, defect concentration, and configurational heat capacity, revealing a second-order order–disorder phase transition near the experimental synthesis temperature. Our study provides the first theoretical insight into the thermodynamics of intrinsic defects in MBT. The negative free energy relative to pristine MBT at synthesis temperatures indicates that Mn–Bi antisite formation is thermodynamically spontaneous. We also present a broadly applicable general framework for correcting low-level theoretical theories using highly accurate many-body corrections from QMC.

Ghaffar, Abdul [ORNL] (ORCID:0000000241190168)↗
DOE OSTI2026

Bridging Transition Metal and Anion Redox Processes in Li-Rich Sulfide Cathodes

Li-ion batteries are essential for decarbonizing global transport and energy, but their scalability is constrained by limited supplies of critical cathode elements, such as Ni, Mn, Co, and P. To address this, we previously introduced high-energydensity Li-ion cathodes composed of Al, Fe, and S, which are elements already produced globally at industrial scale and batterygrade purity. These cathodes leverage sulfide anion redox, involving nonbonding S 3p states and localized distortions that form and break S−S bonds, enabling high capacity. Here, we expand this chemical space by incorporating Cu into cathodes Li 2.2d−z Cu z Al 0.2 Fe 0.6 S 2 (0 ≤ z ≤ 0.4), where highly covalent Cu−S interactions stabilize holes on Cu as Cu >1+ . This Cu redox extends charge compensation that was previously restricted to localized, electronically isolated S−S bonds. Cu also limits capacity, which we attribute to structural destabilization of the delithiated phase, despite the thermodynamic stability of Cu >1+ . By describing the effects of Cu on charge compensation and phase stability, we present a sulfide anion redox mechanism for next-generation multielectron redox Li-ion cathodes, where highly covalent transition metal states participate in otherwise electronically isolated redox processes involving anion nonbonding states.

Anions↗
DOE OSTI2022

Structure, Stability, and Electronic Properties of Boron Suboxide: A Density Functional Theory Study

Boron suboxide (B 6 O) is a boron-rich compound derived from the α-rhombohedral boron lattice with extreme hardness and unusual semiconducting properties. In this work, density functional theory (DFT) was used to show that unit cell volume, mechanical strength, band gaps, and thermodynamic stabilities of B 6 O were influenced by the interstitial elements and point defects at the icosahedral sites. While the hexagonal unit cell volume (HUCV) varies with interstitial occupancy, it is the icosahedral defect that weakens the intrinsic bulk modulus of B 6 O. Using the hybrid HSE functional, we confirmed that the perfect B 6 O bulk is a p-type semiconductor with a direct band gap of 2.8 eV. Furthermore, by screening α-boron compounds systematically, we found that a simple octet rule may offer a consistent explanation for the variations in the computed electronic structures. The formation free energies calculated over a wide range of temperatures (0–2500 K) and pressures (0–80 GPa) predict that formations of interstitial defects become favorable only at higher temperatures (ca. 1800 K) in bulk B 6 O lattices. The nudged elastic band (NEB) method was employed to identify the minimum energy pathways for the diffusions of dislocated B and O atoms. The diffusion of icosahedral B atoms has an energy barrier of 0.16 eV. More complex B diffusion paths involving the reorganization of icosahedral boron atoms incur higher barriers (>1 eV). In contrast, the diffusion of interstitial O atoms is facile with a barrier of 0.4 eV. Lastly, successive O insertions into the α-B lattice were performed using DFT to generate a basic understanding of the oxidation process. Furthermore, these calculations provide fundamental atomistic insights into the growth of B 6 O crystals and control of their point defects.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Insights on the corrosion thermodynamics of chromium in molten LiF-NaF-KF eutectic salts

The thermodynamic stability of chromium in the Cr 0 , Cr 2+ , and Cr 3+ oxidation states considering various F - ion coordination compounds was investigated in molten LiF–NaF-KF (FLiNaK) eutectic salt. Potential fluoride ion activity (F - and CrF 3 - ) diagrams were constructed to predict the most stable Cr oxidation states as a function of anion activity, solvation state of chromium ions and potential at 600 °C. The Gibbs free energies of these compounds were estimated by utilizing a combined cyclic voltammetry – Nernst theory analysis approach in FLiNaK salt. To verify the constructed diagrams, X-ray diffraction was utilized after exposure of Cr at various applied potentials to determine whether compounds detected in solidified FLiNaK salts agreed with thermodynamic calculations. This work aims to identify key thermodynamic factors that are significant to chromium corrosion for molten salt nuclear reactor applications. The F - stability region spans the entire region where spontaneous corrosion of Cr occurs. Cr may be oxidized to Cr 2+ and Cr 3+ in the presence HF (due to moisture as an impurity) spontaneously except some conditions of p 1/2 H 2 /a HF etc. This situation is not changed qualitatively for various states of solvation between oxidized Cr solute in F - solvent and this aspect is substantially similar for the two cases (pair 1: Cr 0 /CrF 3 - / CrF 6 3- ; pair 2: Cr 0 /CrF 4 2- /CrF 5 2- ) considered in this work.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Anticorrelation between electron-phonon coupling strength and stability of ternary metal diborides

This paper endeavors to expand the exploration of prospective superconductivity within the 𝑃⁢6/𝑚⁢𝑚⁢𝑚−MgB 2 structural framework by substituting Mg atoms with a random selection of two distinct cations from a pool of 24 candidates. Employing a high-throughput screening methodology to scrutinize zone-center electron-phonon couplings, we found an inverse relationship between thermodynamic stability and electron-phonon coupling strength in the resulting ternary metal borides. Notably, within the MgNiB 4 system, we demonstrated that electron-phonon superconductivity is stronger in less dynamically stable configurations. Further, the qualitative reasons for the emerging instabilities when electron-phonon coupling becomes strong are discussed. This insight indicates that the quest for superior superconducting materials should concentrate on the intricate balance between electronic and phononic characteristics in structures that are on the brink of instability.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2026

Thermodynamic Properties of Uranium-Zirconium Carbide Solid Solutions with Oxygen

Uranium-zirconium carbide solid solutions are potential candidates for fuel in space nuclear reactors. To support development of advanced robust fuels for space applications, (U y Zr 1-y )C x O z solid solutions were synthesized via carbothermic reduction with uranium, carbon, and oxygen contents raging from y=0.05 – 0.3, x=0.95 – 0.98, and z = 0.02 – 0.07. These solid solutions were subsequently characterized by X-ray diffraction and elemental combustion analysis. High-temperature oxidative drop solution calorimetry in molten sodium molybdate solvent at 1073 K was employed to determine the formation, oxidation, and mixing enthalpies of (U y Zr 1-y )C x O z solid solutions, and thereby for evaluating their thermodynamic stability. The solid solutions are enthalpically less stable than the UC and ZrC end members. They also become less stable and more reactive towards oxidation with increasing uranium content. Incorporation of oxygen into solid solutions increases their solid solution-stability and improves their oxidation resistance.

Calorimetry↗
DOE OSTI2025

Computational Analysis of the Energetic Stability of High-Entropy Structures of a Prototypical Lanthanide-Based Metal–Organic Framework

High-entropy materials are characterized by their complex compositions, typically comprising five or more elements in near-equiatomic proportions. Applying this concept to metal ions in metal−organic frameworks (MOFs) has paved the way for exploring a new class of high-entropy MOFs. While the compositional strategy of high-entropy materials leverages configurational entropy to aid thermodynamic stability, it also poses significant analytical challenges due to the vast compositional landscape and diverse phases that these materials can adopt. We present a computational study of several complexities associated with selecting potential high-entropy versions of a prototype lanthanidebased MOF. We compute the energetics of metal mixing of these heterometallic MOFs using density functional theory (DFT) and machine learning interatomic potential (MLIP) methods. The use of MLIP methods allows a systematic exploration of the convex hull of thermodynamically stable MOF structures containing up to 5 distinct metals.

Chemical structure↗
DOE OSTI2025

A framework to evaluate machine learning crystal stability predictions

The rapid adoption of machine learning in various scientific domains calls for the development of best practices and community agreed-upon benchmarking tasks and metrics. We present Matbench Discovery as an example evaluation framework for machine learning energy models, here applied as pre-filters to first-principles computed data in a high-throughput search for stable inorganic crystals. We address the disconnect between (1) thermodynamic stability and formation energy and (2) retrospective and prospective benchmarking for materials discovery. Alongside this paper, we publish a Python package to aid with future model submissions and a growing online leaderboard with adaptive user-defined weighting of various performance metrics allowing researchers to prioritize the metrics they value most. To answer the question of which machine learning methodology performs best at materials discovery, our initial release includes random forests, graph neural networks, one-shot predictors, iterative Bayesian optimizers and universal interatomic potentials. We highlight a misalignment between commonly used regression metrics and more task-relevant classification metrics for materials discovery. Accurate regressors are susceptible to unexpectedly high false-positive rates if those accurate predictions lie close to the decision boundary at 0 eV per atom above the convex hull. The benchmark results demonstrate that universal interatomic potentials have advanced sufficiently to effectively and cheaply pre-screen thermodynamic stable hypothetical materials in future expansions of high-throughput materials databases.

Riebesell, Janosh↗
DOE OSTI2022

Tuning Enzyme Thermostability via Computationally Guided Covalent Stapling and Structural Basis of Enhanced Stabilization

Enhancing the thermostability of enzymes without impacting their catalytic function represents an important yet challenging goal in protein engineering and biocatalysis. We recently introduced a novel method for enzyme thermostabilization that relies on the computationally guided installation of genetically encoded thioether “staples” into a protein via cysteine alkylation with the noncanonical amino acid O-2-bromoethyl tyrosine (O2beY). Here, we demonstrate the functionality of an expanded set of electrophilic amino acids featuring chloroacetamido, acrylamido, and vinylsulfonamido side-chain groups for protein stapling using this strategy. Using a myoglobin-based cyclopropanase as a model enzyme, our studies show that covalent stapling with pchloroacetamido-phenylalanine (pCaaF) provides higher stapling efficiency and enhanced stability (thermodynamic and kinetic) compared to the other stapled variants and the parent protein. Interestingly, molecular simulations of conformational flexibility of the cross-links show that the pCaaF staple allows fewer energetically feasible conformers than the other staples, and this property may be a broader indicator of stability enhancement. Using this strategy, pCaaF-stapled variants with significantly enhanced stability against thermal denaturation (ΔT m ' = +27 °C) and temperature-induced heme loss (ΔT 50 = +30 °C) were obtained while maintaining high levels of catalytic activity and stereoselectivity. Crystallographic analyses of singly and doubly stapled variants provide key insights into the structural basis for stabilization, which includes both direct interactions of the staples with protein residues and indirect interactions through adjacent residues involved in heme binding. This work expands the toolbox of protein stapling strategies available for protein stabilization.

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

High throughput exploration of the oxidation landscape in high entropy alloys

We report high entropy alloys (HEAs) have gained interest for structural applications in extreme environments. With a potentially vast chemical and phase space, there are significant opportunities to discover superior performing alloys. Crucial for most high-temperature applications is understanding and mitigating the oxidation behavior of these chemically complex alloys. Most experimental and computational HEA studies have focused on a limited set of compositions and only a fraction these compositions have been characterized for oxidation. We present a high-throughput framework that utilizes density-functional theory (DFT) in concert with a combined machine-learning model and grand-canonical linear programming for assessing phase stability, phase-fraction, chemical activity and high-temperature survivability of arbitrary HEAs. This framework considers temperature dependent contributions to the Gibbs energy of the competing phases arising from short-range order and vibrational entropy. We demonstrate the effectiveness of the framework by assessing the thermodynamic stability, oxidation behavior, chemical activity, and phase decomposition of body-centered cubic Mo-W-Ta-Ti-Zr refractory HEAs. A total of 51 compositions were analyzed and ranked in order of their survivability based on the Pareto-front analysis. Oxidation was performed at 1373 K on fours samples in air showing the difference in oxidation behavior determined experimentally through scale thickness and their mass changes. The insights on oxidation behavior presented in this work will enable the fast assessment of technologically useful HEAs needed for future structural application in extreme conditions

36 MATERIALS SCIENCE↗
DOE OSTI2020

The significant impact of aerosol vertical structure on lower atmosphere stability and its critical role in aerosol–planetary boundary layer (PBL) interactions

The aerosol–planetary boundary layer (PBL) interaction was proposed as an important mechanism to stabilize the atmosphere and exacerbate surface air pollution. Despite the tremendous progress made in understanding this process, its magnitude and significance still have large uncertainties and vary largely with aerosol distribution and meteorological conditions. In this study, we focus on the role of aerosol vertical distribution in thermodynamic stability and PBL development by jointly using micropulse lidar, sun photometer, and radiosonde measurements taken in Beijing. Despite the complexity of aerosol vertical distributions, cloud-free aerosol structures can be largely classified into three types: well-mixed, decreasing with height, and inverse structures. The aerosol–PBL relationship and diurnal cycles of the PBL height and PM 2.5 associated with these different aerosol vertical structures show distinct characteristics. The vertical distribution of aerosol radiative forcing differs drastically among the three types, with strong heating in the lower, middle, and upper PBL, respectively. Such a discrepancy in the heating rate affects the atmospheric buoyancy and stability differently in the three distinct aerosol structures. Absorbing aerosols have a weaker effect of stabilizing the lower atmosphere under the decreasing structure than under the inverse structure. As a result, the aerosol–PBL interaction can be strengthened by the inverse aerosol structure and can be potentially neutralized by the decreasing structure. Moreover, aerosols can both enhance and suppress PBL stability, leading to both positive and negative feedback loops. This study attempts to improve our understanding of the aerosol–PBL interaction, showing the importance of the observational constraint of aerosol vertical distribution for simulating this interaction and consequent feedbacks.

54 ENVIRONMENTAL SCIENCES↗