Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase stability modeling”

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

Realization of two-dimensional discrete time crystals with anisotropic Heisenberg coupling

A discrete time crystal (DTC) is an out-of-equilibrium phase of matter that spontaneously breaks discrete time-translation symmetry. Previous studies have been limited to a set of models with Ising-like couplings - and mostly only in one dimension - thus precluding our understanding of the existence (or not) of DTCs in models with more realistic interactions. In this work, by combining the latest generation of IBM quantum processors with state-of-the-art tensor network methods, we demonstrate the existence of a DTC in a two-dimensional system governed by anisotropic Heisenberg interactions. We uncover a rich phase diagram encompassing spin-glass, ergodic, and time-crystalline phases, and identify the interplay of initialization, interaction anisotropy, and driving protocols in stabilizing the DTC phase. By extending the study of Floquet matter beyond simplified models, we lay the groundwork for exploring how driven systems bridge the gap between quantum coherence and emergent non-equilibrium thermodynamics.

Phase transitions and critical phenomena↗
DOE OSTI2024

Phase Pattern Formation in Grain Boundaries

Here, in this Letter, we derive conditions that predict the existence of two-phase periodic-pattern grain boundary structures that are stable against coarsening. While previous research has established that elastic effects can lead to phase pattern formation on crystal surfaces, the possibility of stable grain boundary structures composed of alternating grain boundary phases has not been previously analyzed. Our theory identifies the specific combination of grain boundary and materials properties that enable the emergence of patterned grain boundary states and shows that the dislocation content of grain boundary phase junctions, absent in surface phenomena, weakens the stability of the patterned structures. The predictions of the theory are tested using a model copper grain boundary that exhibits multiple phases and two-phase pattern formation. We discuss how, similarly to surfaces, elastic effects associated with grain boundary phase junctions have profound implications for how grain boundary phases transform.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Energetics of hydroxylbastnäsite solid solutions, La 1-$\chi$ Nd $\chi$ CO 3 OH

Bastnäsites (LnCO 3 (F,OH)) are a group of common rare earth elements (REE)-bearing minerals and are one of the primary global sources of REE. Due to the chemical similarities among REE, bastnӓsites tend to occur as solid solutions instead of end members in REE containing ores. To better understand the processes and the mechanisms of formation of such deposits, it is essential to determine the thermodynamic properties of bastnӓsites, including hydroxylbastnӓsite (LnCO 3 OH) solid solutions. In this work, we performed detailed structural and calorimetric investigations on synthetic hexagonal La–Nd hydroxylbastnӓsite (La 1-$\chi$ NdCO 3 OH, $\chi$ = 0, 0.25, 0.5, 0.75, 1) solid solutions. X-ray diffraction confirms the crystal structure of the solid solution series in the $\ P$6 space group, and pair distribution function (PDF) analysis reveals local bonding environments characterized by three different types of 9-coordinated metal-oxygen polyhedra. Unit cell parameters of La1–xNdxCO3OH exhibit a nearly linear relation with the Nd content x, suggesting a random distribution of La and Nd in the structure. Their standard enthalpies of formation (Δ$\ H$° f ) were determined by high temperature oxide melt drop solution calorimetry, from which the enthalpies of mixing (Δ$\ H$ mix ) were derived. The Δ$\ H$ mix can be fitted by a regular solution model with an interaction parameter of 12.58 ± 0.16 kJ/mol, suggesting enthalpic metastability of La1–xNdxCO3OH relative to the two endmembers. Combining entropy and enthalpy, we further estimated the Gibbs free energies of mixing (Δ$\ G$ mix ) at relevant temperatures, revealing favorable temperatures under which the intermediate La 1-$\chi$ NdCO 3 OH phases can be stabilized. Such entropy-driven stabilization, as is consistent with our geochemical modeling results, may explain the enhancement of thermal stability of the solid solutions in nature. Additionally, the temperature range constrained from this study may be used to estimate the thermal history of REE bastnӓsite deposit.

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

Compositional Tuning of Magnetic Properties in a Series of Transition Metal Site-Deficient UCo x Bi 2 and UNi x Bi 2 Phases

Varying the electronic structure of topological materials through aliovalent substitution is a primary approach to tuning their physical properties. Unlike substitution, metal site deficiency intrinsic to some structure types, including HfCuSi 2 -type, has rarely been employed for controlling the properties of topological phases. In this report, we describe the synthesis and characterization of two new series of compounds, UCoxBi 2 and UNixBi 2 , which demonstrate the variation of transition metal content through synthetic conditions. Magnetic measurements reveal the dependence between the extent of transition metal incorporation and the magnetism of the resulting phase. DFT calculations demonstrated the ability to model their formation and predict the stability ranges of transition metal-site deficient compounds.

Crystal structure↗
DOE OSTI2025

In Situ/Operando Probing of Dynamic Phase Structures of Alumina‐Supported Ultrasmall Copper‐Gold Alloy Nanoparticles Under Reaction Conditions

Abstract The ability to control phase structures and surface sites of ultrasmall alloy nanoparticles under reaction conditions is essential for preparing catalysts by design. This is, however, challenging due to limited understanding of the atomic‐scale phases and their correlation with the ensemble‐averaged structures and activities of catalysts during catalytic reactions. We reveal here a dynamic structural stability of alumina‐supported ultrasmall and equiatomic copper‐gold alloy nanoparticles under reaction conditions as a model system in the in situ/operando study. In situ atomic‐scale morphological tracking under oxygen reveals temperature‐dependent dynamic crystalline‐amorphous dual‐phase structures, showing dynamic stability over an elevated temperature range. This atomic‐scale dynamic phase stability coincides with a “conversion plateau” observed for carbon monoxide oxidation on the catalyst. It is substantiated by the stable lattice ordering/disordering structures and surface sites with oscillatory characteristics shown by operando ensemble‐average structural tracking of the catalyst during the oxidation reaction. The understanding of the atomic‐scale dynamic phase structures in correlation with the ensemble‐average dynamic ordering/disordering phase structures and surface sites provides fresh insights into the unique synergy of the supported alloy nanoparticles. This understanding has implications for the design and structural tuning of active and stable ultrasmall alloy catalysts under elevated temperatures.

Chemistry↗
DOE OSTI2023

Impedance Modeling and Stability Analysis for the PLL-less and Voltage Sensor-less Grid-tied Converters

The newly developed control methodology, which omits the phase-locked loop (PLL) and voltage sensor, has garnered considerable attention for its rapid dynamic response and reduced hardware expenses. Nonetheless, its stability concerns have not yet been thoroughly scrutinized. In this paper, the impedance modeling of the PLL-less and voltage sensorless control strategy has been established. The impedance model has been contrasted with the traditional PLL-based control scheme. Our findings demonstrate that the PLL-less and voltage sensorless control approach lacks the negative resistance behavior that is frequently observed in the traditional scheme, which results in unstable operations. Furthermore, the role of virtual resistance has been examined. To substantiate the aforementioned findings, simulation results are presented.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

A stress-sensitive precipitate nucleation model beyond classical nucleation theory

The dynamic evolution of precipitates and second phases dictates the strength and stability of most engineering alloys. By design, or as a consequence of thermo-mechanical aging, engineering metals and alloys often form precipitates of second phases when subjecting to diverse thermal and mechanical loads. Precipitation is governed by several factors, including the alloy’s composition, processing/operating temperature, and stresses — either as a result of external loads or from residual stresses. However, state-of-the-art models for precipitate nucleation (i.e., classical nucleation theory) typically lacks consistent method to capture the effects of externally applied and/or internal stresses on nucleation; thereby severely limiting the applicability of these models to complex materials systems and to representative loading scenarios. Here, in this work, we extend upon classical nucleation theory to account for the effect of stresses on precipitation kinetics and thermodynamics. This is achieved via the use of an Eshelbian micromechanics framework keeping track of (i) the stress build up resulting from second phase formation as a function of mechanical load and, (ii) the effects of dislocations on precipitate formation. This new model is applied to σ precipitate in Fe–Cr binary alloys and M 23 C 6 precipitate in 316H stainless steel (SS). Simulations demonstrate the important role of both the remotely applied loads and dislocation pile ups on precipitate nucleation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Liquid Phase Modeling in Porous Media: Adsorption of Methanol and Ethanol in H-MFI in Condensed Water

Zeolites are used in the chemical and separation industries for their exceptional selectivity, adsorption capacity, regenerability, and stability in gas and liquid phase processing. Here, we developed an explicit solvation method for predicting solvent/condensed phase effects on adsorption free energies in microporous media such as zeolites based on the hybrid quantum mechanical/molecular mechanical free energy perturbation (QM/MM-FEP) technique. Our explicit solvation method for zeolite systems, called eSZS, aims to capture site-specific interactions during the adsorption process at the Brønsted acid sites of H-MFI zeolite while still considering the diverse configuration space of the solvent molecules. This strategy is ideal for chemical reactions or adsorbates that interact with the microporous medium in few distinct adsorbate/transition state configurations, i.e., the harmonic or similar approximations are acceptable for the adsorbate/transition state while such approximations break down for the solvent molecules that require extensive configuration space sampling. In this way, our approach effectively overcomes the limitations of implicit solvation models and classical force field methods for describing solvation effects on chemical reactions within porous materials such as zeolites. Specifically, in this study, we investigated various aspects of our hybrid QM/MM approach, including QM cluster size dependencies in a periodic electrostatically embedded cluster model (PEECM), rules for link atoms at the QM/MM boundary, and functional and basis set considerations for converged and reasonably accurate gas and aqueous phase methanol and ethanol adsorption free energy predictions in H-MFI. For gas phase adsorption of methanol and ethanol in H-MFI at a Brønsted acid site in T12 position, we compute adsorption free energies at 298 K of −0.61 and −0.75 eV, respectively, using a PEECM containing 50 Si and 1 Al atom with ωB97x-D/def2-TZVP level of theory. For solvent effect calculations, we sample the aqueous phase using grand canonical Monte Carlo (GCMC) simulations to (1) obtain a mean field of electrostatic interactions in the reaction system and (2) perform a rigorous free energy perturbation calculation. Similar to the experimentally and computationally observed endergonic solvation effects observed for hydrocarbon adsorption on metal surfaces, we also observe that a condensed aqueous environment destabilizes methanol and ethanol at these acid sites in H-MFI at 298 K. Specifically, the computed solvation free energies of adsorption (ΔΔG solv ) for methanol and ethanol are +0.44 and +0.54 eV, respectively. From this study, it is evident that adsorbates (methanol and ethanol) are competing with water for adsorption space inside the H-MFI zeolite, leading to an endergonic solvation effect. Here, we expect that the endergonic, aqueous solvent effect during adsorption in microporous zeolites is highly tunable by changing the pore size and hydrophobicity of the microporous material as this will affect the water density inside the pore structure.

Adsorption↗
DOE OSTI2023

Intermittent Defect Fluctuations in Oxide Heterostructures

Abstract The heterogeneous nature, local presence, and dynamic evolution of defects typically govern the ionic and electronic properties of a wide variety of functional materials. While the last 50 years have seen considerable efforts into development of new methods to identify the nature of defects in complex materials, such as the perovskite oxides, very little is known about defect dynamics and their influence on the functionality of a material. Here, the discovery of the intermittent behavior of point defects (oxygen vacancies) in oxide heterostructures employing X‐ray photon correlation spectroscopy is reported. Local fluctuations between two ordered phases in strained SrCoO x with different degrees of stability of the oxygen vacancies are observed. Ab‐initio‐informed phase‐field modeling reveals that fluctuations between the competing ordered phases are modulated by the oxygen ion/vacancy interaction energy and epitaxial strain. The results demonstrate how defect dynamics, evidenced by measurement and modeling of their temporal fluctuations, give rise to stochastic properties that now can be fully characterized using coherent X‐rays, coupled for the first time to multiscale modeling in functional complex oxide heterostructures. The study and its findings open new avenues for engineering the dynamical response of functional materials used in neuromorphic and electrochemical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Causal horizons, geodesic completeness and stability in slow contraction cosmology

We show that cosmological models with a semi-infinite phase of slow contraction (ekpyrosis) possess a combination of properties that can address several fundamental problems in cosmology, otherwise faced in contracting de Sitter phases or standard big bang expansion. In particular, flat or open slow contraction admits a stable past attractor that asymptotes to Minkowski space and is past geodesically complete, as well as a stable, flat, homogeneous, and isotropic future attractor with negligible Weyl curvature (and, therefore, negligible gravitational entropy). In bouncing cosmologies, this contracting attractor is terminated by a smooth, non-singular bounce that transforms the attractor properties at the end of contraction into the initial conditions for the subsequent expanding phase. Cosmologies incorporating a slow contraction phase have no particle horizon and therefore avoid the causal horizon problem. The past Minkowski attractor also generates an initial spectrum of vacuum-like quantum fluctuations on all wavelengths. Moreover, because the averaged expansion rate along past-directed geodesics is non-positive, models incorporating a semi-infinite phase of slow contraction also evade the Borde–Guth–Vilenkin theorem. By contrast, contracting de Sitter space possesses a finite particle horizon and becomes unstable in the presence of scalar fields, matter, or radiation.

Cosmology and Nongalactic Astrophysics (astro-ph.C↗
DOE OSTI2020

Analytic fits to atom-in-jellium shear modulus predictions

Fits to AJ shear modulus calculations from ρ0 to the limit of the corresponding AJ EOS table (usually 1000 ρ0), adjusted to reproduce STP shear modulus where possible. Fitting minimized the fractional difference, so low-pressure points would be reproduced accurately despite the wide-ranging tabulation. AJ calculations were at T=0. Isochoric variation G(T) has been small for all cases spot-tested so far, and may have either sign. The fitted equation usually matches the AJ data to within a few percent. Between numerical noise in the AJ calculation and probably-physical structure not captured by the equation, the deviation could be up to 20% in some places in most models, and 30% in a few. AJ does not capture structural phase transitions or phases stabilized by directional bonds, where the shear modulus may vary by a greater amount. Where possible, STP ρ,G were used as parameters, and low pressure AJ points were de-weighted or removed if necessary. Otherwise, such as where AJ fails to capture solid phases with a significantly different shear modulus, the AJ data were fitted as far down in pressure as possible, G0 was also fitted if necessary, and ρ0 was also adjusted if needed to keep G0>0. The resulting models are not intended for use at low pressure, though some are probably adequate for practical purposes

36 MATERIALS SCIENCE↗
DOE OSTI2024

Correction to “Frank–Kasper Phases of Diblock Copolymer Melts Studied with the DPD Model: SCF Results”

In a recent study, (1) we performed the self-consistent field (SCF) calculations of a dissipative particle dynamics (DPD) model of conformationally asymmetric diblock copolymer A-B melts to study the stability of various Frank–Kasper (FK) phases formed by such A-B melts. While all the equations in ref (1) are correct, we did not calculate the dimensionless (mean-field) Helmholtz free energies per chain βf c = βu c,χ + βu c,κ – s c /k B and the dimensionless stresses dβf c /dθ of ordered phases accurately. For example, the dimensionless (nonbonded) internal energy per chain due to the system compressibility βu c,κ = (N/2κ)∑ q βû 0 (|q|)(ϕ̂ A (q) + ϕ̂ B (q) – δ q,0 )(ϕ̂ A (−q) + ϕ̂ B (−q) – δ– q,0 ) should be calculated in the reciprocal space as given (note that βû 0 (|q|) is obtained analytically, and that since ϕ A (r) is real 𝜙̂ A (−𝐪)= $\overline{𝜙̂ A}$(𝐪), i.e., the complex conjugate of ϕ̂ A (q)); this corresponds to numerically evaluating βu c,κ = (N/2κ)∫dr(ϕ A (r) + ϕ B (r) – 1)ψ(r) with ψ(r) ≡ (∫dr′βu 0 (|r–r′|)(ϕ A (r′) + ϕ B (r′) – 1)/V in the real space via the (composite) trapezoidal rule. While the trapezoidal rule is less accurate than the Romberg integration (2) in general, the opposite occurs when the integration domain is periodic; this was pointed out by Matsen in ref (3) and explained in detail in refs (4and5). With our ignorance of this exceptional convergence property of the trapezoidal rule for periodic functions, we calculated all convolutions having the form of ψ(r) in the reciprocal space via the fast Fourier transforms (6) but those having the form of βu c,κ in the real space via Romberg integration; the latter leads to inaccurate results and unfortunately incorrect phase diagrams reported in ref (1), which are corrected here.

He, Juntong↗
DOE OSTI2023

Gas-Phase Composition as a Predictive Metric for Calendar Life Behavior of Next-Generation Silicon Anodes

The expansion of renewable technologies and electrification of the transportation sector is driving increased demand for next-generation battery materials that provide higher power and energy density with superior cycling and calendar life stability. Silicon (Si) has a theoretical capacity nearly 10x that of graphite, and is therefore a promising anode material candidate to meet these rigorous performance demands. While leading Si anode battery demonstrations are approaching target metrics for cycle life, a series of complex and interrelated modes of reactivity lead to reduced calendar life and therefore challenge practical adoption of these materials. Deconvoluting the degradation processes that impact Si calendar life is critical to informing the rational and accelerated design of improved Si materials. In the present work, we employ novel sampling techniques and GC-MS-FID characterization to measure gas-phase composition during initial Si cycling, which we tie to selective mechanisms of Si passivation. We utilize a tiered analysis approach to identify and quantify the gas-phase reaction products associated with three advanced Si material candidates under practical operating conditions. Ex situ analysis of Si powders (pure chemical reactivity) is coupled with nondestructive in situ sampling of Si electrodes in a practical pouch-cell format (coupled chemical and electrochemical reactivity). We link the observed gas-phase species evolution to electrochemical behavior and measured calendar life of the three Si materials. Further, we evaluate the voltage-resolved evolution of gas-phase species for one such Si nanomaterial, where nonmonotonic gas generation implies competition between passivating reaction pathways. The measured gas-phase compositional data serves as a critical input for our advanced electrochemical SEI models to identify favorable vs unfavorable reaction pathways to stabilize Si. In addition to bolstering a fundamental understanding of Si reactivity, the present approach informs specific and quantifiable gas-phase metrics tied to calendar life improvements in Si, which can streamline and accelerate the process of next-generation material development.

DIRECT ENERGY CONVERSION,ENERGY STORAGE↗
DOE OSTI2025

Charting the chemical space of Zintl phases with graph neural networks and bonding insights

A large number of Zintl phases have been discovered by solid-state chemists driven by empirical knowledge, chemical intuition and in some cases, through serendipitous accidents. These discoveries have only scratched the surface, given the vast compositional and structural diversity that Zintl phases can accommodate. The large chemical space of Zintl phases, as well as intermetallic compounds in general, remain under-explored. Here, we use graph neural networks and the upper bound energy minimization approach to efficiently scan a large chemical space of >90 000 hypothetical Zintl phases and accurately discover 1810 new thermodynamically stable phases with 90% precision, as validated with first-principles calculations. We show that our approach is more than 2× more accurate in predicting DFT stability than M3GNet (40% precision) on the same dataset. Using a random forest model and SHAP analysis, we demonstrate the critical role of ionic bonding in the thermodynamic stability of Zintl phases. Our results not only expand the known chemical landscape of Zintl phases but also highlight the efficacy of machine learning frameworks combined with domain knowledge in uncovering chemically meaningful insights across complex intermetallics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Controlling the Phase Transition in CsPbI 3 Nanowires

Cesium lead iodide (CsPbI 3 ) is a promising semiconductor with a suitable band gap for optoelectronic devices. CsPbI 3 has a metastable perovskite phase that undergoes a phase transition into an unfavorable nonperovskite phase in an ambient environment. This phase transition changes the optoelectronic properties of CsPbI 3 and hinders its potential for device applications. Therefore, it is of central importance to understand the kinetics of such instability and develop strategies to control and stabilize the perovskite phase. Here, we use ultralong CsPbI 3 nanowires as a model platform to investigate the phase transition kinetics. Our results depict the role of environmental stressors (moisture and temperature) in controlling the phase transition dynamics of CsPbI 3 , which can serve as guiding principles for future phase transition studies and the design of related photovoltaics. Furthermore, we demonstrate the controllability of phase propagation on individual nanowires by varying the moisture level and temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Thermal equation of state of ice-VII revisited by single-crystal X-ray diffraction

Abstract Ice-VII is a high-pressure polymorph of H2O ice and an important mineral widely present in many planetary environments, such as in the interiors of large icy planetary bodies, within some cold subducted slabs, and in diamonds of deep origin as mineral inclusions. However, its stability at high pressures and high temperatures and thermoelastic properties are still under debate. In this study, we synthesized ice-VII single crystals in externally heated diamond-anvil cells and conducted single-crystal X-ray diffraction experiments up to 78 GPa and 1000 K to revisit the high-pressure and high-temperature phase stability and thermoelastic properties of ice-VII. No obvious unit-cell volume discontinuity or strain anomaly of the high-pressure ice was observed up to the highest achieved pressures and temperatures. The volume-pressure-temperature data were fitted to a high-temperature Birch-Murnaghan equation of state formalism, yielding bulk modulus KT0 = 21.0(4) GPa, its first pressure derivative KT0′ = 4.45(6), dK/dT = –0.009(4) GPa/K, and thermal expansion relation αT = 15(5) × 10–5 + 15(8) × 10–8 × (T – 300) K–1. The determined phase stability and thermoelastic properties of ice-VII can be used to model the inner structure of icy cosmic bodies. Combined with the thermoelastic properties of diamonds, we can reconstruct the isomeke P-T paths of ice-VII inclusions in diamond from depth, offering clues on the water-rich regions in Earth’s deep mantle and the formation environments of those diamonds.

Geochemistry & Geophysics↗
DOE OSTI2025

First‐Order Empirical Interpolation Method for Real‐Time Solution of Parametric Time‐Dependent Nonlinear PDEs

ABSTRACT We present a model reduction approach for the real‐time solution of time‐dependent nonlinear partial differential equations (PDEs) with parametric dependencies. A major challenge in constructing efficient and accurate reduced‐order models for nonlinear PDEs is the efficient treatment of nonlinear terms. We address this by unifying the implementation of hyperreduction methods to deal with nonlinear terms. Furthermore, we introduce a first‐order empirical interpolation method (EIM) to provide an efficient approximation of the nonlinear terms in time‐dependent PDEs. We demonstrate the effectiveness of our approach on the Allen–Cahn equation, which models phase separation, and the Buckley–Leverett equation, which describes two‐phase fluid flow in porous media. Numerical results highlight the accuracy, efficiency, and stability of the proposed method compared with both the Galerkin–Newton approach and hyper‐reduced models using the standard EIM.

Nguyen, Ngoc Cuong [Center for Computational Engin↗