Search NASA⌕ Search

SEARCH · Search NASA

Results for “Phase diagrams”

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.

291 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Rethinking 𝛼−RuCl 3 : Parameters, models, and phase diagram

RuCl 3 was likely the first ever deliberately synthesized ruthenium compound, following the discovery of the 44 Ru element in 1844. For a long time it was known as an oxidation catalyst, with its physical properties being discrepant and confusing, until a decade ago when its allotropic form 𝛼−RuCl 3 rose to exceptional prominence. This “rediscovery” of 𝛼−RuCl 3 has not only reshaped the hunt for a material manifestation of the Kitaev spin liquid, but it has opened the floodgates of theoretical and experimental research in the many unusual phases and excitations that the anisotropic-exchange magnets as a class of compounds have to offer. Given its importance for the field of Kitaev materials, it is astonishing that the low-energy spin model that describes this compound and its possible proximity to the much-desired spin-liquid state is still a subject of significant debate ten years later. In the present study, we argue that the existing key phenomenological observations put strong natural constraints on the effective microscopic spin model of 𝛼−RuCl 3 , and specifically on its spin-orbit-induced anisotropic-exchange parameters that are responsible for the nontrivial physical properties of this material. These constraints allow one to focus on the relevant region of the multidimensional phase diagram of the 𝛼−RuCl 3 model, suggest an intuitive description of it via a different parametrization of the exchange matrix, offer a unifying view on the earlier assessments of its parameters, and bring closer together several approaches to the derivation of anisotropic-exchange models. We explore extended phase diagrams relevant to the 𝛼−RuCl 3 parameter space using quasiclassical, Luttinger-Tisza, exact diagonalization, and density-matrix renormalization-group methods, demonstrating a remarkably close quantitative accord between them on the general structure and hierarchy of the phases, with the zigzag, ferromagnetic, and incommensurate phases that are proximate to each other. As a result, one of the highlights is the detailed agreement on the nature of the incommensurate phases that realize two distinct counterrotating helical states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

The phase diagram of quantum chromodynamics in one dimension on a quantum computer

The quantum chromodynamics (QCD) phase diagram, which reveals the state of strongly interacting matter at different temperatures and densities, is key to answering open questions in physics, ranging from the behaviour of particles in neutron stars to the conditions of the early universe. However, classical simulations of QCD face significant computational barriers, such as the sign problem at finite matter densities. Quantum computing offers a promising solution to overcome these challenges. Here, we take an important step toward exploring the QCD phase diagram with quantum devices by preparing thermal states in one-dimensional non-Abelian gauge theories. We experimentally simulate the thermal states of SU(2) and SU(3) gauge theories at finite densities on a trapped-ion quantum computer using a variational method. This is achieved by introducing two features: Firstly, we add motional ancillae to the existing qubit register to efficiently prepare thermal probability distributions. Secondly, we introduce charge-singlet measurements to enforce colour-neutrality constraints. This work pioneers the quantum simulation of QCD at finite density and temperature for two and three colours, laying the foundation to explore QCD phenomena on quantum platforms.

Quantum information↗
DOE OSTI2026

Superconducting phase diagram of multilayer square-planar nickelates

The discovery of superconductivity in square-planar nickelates has offered a rich materials platform to explore the origins of high-temperature superconductivity. However, experimental investigations have largely been limited to the infinite-layer R NiO 2 ( R , rare earth) nickelates. For this work we constructed a phase diagram of multilayer square-planar Nd n+1 Ni n O 2n+2 compounds and found signatures of superconductivity for dimensionality n = 4 to 8. Upon decreasing n , the superconducting anisotropy evolves owing to 4ƒ electron effects, and electronic structure characteristics approach cuprate-like behavior. Magnetic fluctuations persist from within the superconducting regime and into the overdoped, nonsuperconducting regime. The superconducting regime overlaps with that of chemically doped infinite-layer nickelates, demonstrating underlying commonalities as well as differences across varying structural realizations of square-planar nickelates. Our work establishes this layered template for creating new nickel-based superconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Strong Kitaev Interaction in BaCo 2⁢ (AsO 4 ) 2

The inelastic neutron scattering results and their analysis unequivocally point to a dominant Kitaev interaction in the honeycomb-lattice cobaltate BaCo 2 ⁢(AsO 4 ) 2 . Our anisotropic-exchange model closely describes all available neutron scattering data in the material’s field-polarized phase. Furthermore, the density-matrix renormalization group results for our model are in close accord with the unusual double-zigzag magnetic order and the low in-plane saturation field of BaCo 2 ⁢(AsO 4 ) 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Metrology for femtosecond pulsed x-ray heating in diamond anvil cell experiments at the European XFEL: Revisiting the iron phase diagram up to 150 GPa

The development of pulsed intense x-ray sources, such as free electron laser, offers new avenues for high pressure experiments. Here, we study the feasibility and metrology of x-ray heating in diamond anvil cells at the European x-ray free electron laser. This method enables one to volumetrically heat the sample while inhibiting chemical migration and probing the crystallographic structure of the sample throughout the heating with a high repetition rate. We focus our study on iron, whose phase diagram is well established up to 100 GPa, to explore the possibilities and limitations of this technique. We volumetrically heat iron samples at starting pressures ranging from 10 to 138 GPa, using the x-ray beam pulsed at 4.5 MHz in a serial pump-and-probe experimental design. Experimental challenges arise from temperature gradients within the sample, changes in temperature at the 100 ns timescale, the difficulty of direct temperature estimates, the effect of thermal pressure, and the presence of metastable crystallites due to rapid cycles of heating and cooling. Hence, we develop a multi-crystal-like data processing method that allows us to account for sample heterogeneity in probed conditions. We then calibrate our measurements using known physical properties of iron under pressure. Thermal pressure in our experiments increases from 4% of the isochoric prediction at 10 GPa to 23% at 138 GPa, and we show that our data are in agreement with most previous observations of iron in this pressure range. The method can now be implemented at higher pressures and temperatures and on materials with unknown phase diagrams.

Materials science↗
DOE OSTI2026

Topological magnons on the ferromagnetic zigzag lattice

Motivated by the experimental identification of magnetic compounds consisting of zigzag chains, we analyze the band structure topology of magnons in ferromagnets on a zigzag lattice. We account for the general lattice geometry by including spatially anisotropic Heisenberg exchange interactions and by Dzyaloshinskii-Moriya interaction on inversion asymmetric bonds. Within the linear spin-wave theory, we find two magnon branches, whose band structure topology (i.e., Chern numbers) we map out in a comprehensive phase diagram. Notably, besides topologically trivial and gapless phases, we identify topologically nontrivial phases that support chiral edge magnons. We show that these edge states are robust against elastic defect scattering.

Subramanian, Skandan [ORNL] (ORCID:000900087255487↗
DOE OSTI2023

Diagnostic diagrams for ram pressure stripped candidates

ABSTRACT This paper presents a method for finding ram pressure stripped (RPS) galaxy candidates by performing a morphological analysis of galaxy images obtained from the Legacy survey. We consider a sample of about 600 galaxies located in different environments such as groups and clusters, tidally interacting pairs and the field. The sample includes 160 RPS previously classified in the literature into classes from J1 to J5, based on the increasing level of disturbances. Our morphological analysis was done using the astromorphlib software followed by the inspection of diagnostic diagrams involving combinations of different parameters like the asymmetry (A), concentration (C), Sérsic index (n), and bulge strength parameters $F(G,\, M_{20})$. We found that some of those diagrams display a distinct region in which galaxies classified as J3, J4, and J5 decouples from isolated galaxies. We call this region as the morphological transition zone and we also found that tidally interacting galaxies in pairs are predominant within this zone. Nevertheless, after visually inspecting the objects in the morphological transition zone to discard obvious contaminants, we ended up with 33 bona fide new RPS candidates in the studied nearby groups and clusters (Hydra, Fornax, and CLoGS sample), of which one-third show clear evidence of unwinding arms. Future works may potentially further increase significantly the samples of known RPS using such method.

Astronomy & Astrophysics↗
DOE OSTI2025

Spin-strain interactions under hydrostatic pressure in α-RuCl 3

We investigate the effects of hydrostatic pressure on 𝛼−RuCl 3 , a prototypical material for the Kitaev spin model on a honeycomb lattice with a possible spin-liquid ground state. Using ultrasound measurements at pressures up to 1.16 GPa, we reveal significant modifications of the acoustic properties and the 𝐻−𝑇 phase diagram of this material. Hydrostatic pressure suppresses the three-dimensional magnetic order and induces a dimerization transition at higher pressures. At low pressures, the sound attenuation exhibits a linear temperature dependence, while above 0.28 GPa, it becomes nearly temperature independent, suggesting a shift in the phonon scattering regime dominated by Majorana fermions. These findings provide new insights into spin-strain interactions in Kitaev magnets and deliver a detailed characterization of the 𝐻−𝑇 phase diagram of 𝛼−RuCl 3 under hydrostatic pressure.

Focused ion beam↗
DOE OSTI2024

Complex orders and chirality in the classical Kitaev-Γ model

It is well recognized that the low-energy physics of many Kitaev materials is governed by two dominant energy scales, the Ising-type Kitaev coupling 𝐾 and the symmetric off-diagonal Γ coupling. An understanding of the interplay between these two scales is therefore the natural starting point toward a quantitative description that includes subdominant perturbations that are inevitably present in real materials. This study focuses on the classical 𝐾−Γ model on the honeycomb lattice, with a specific emphasis on the region 𝐾< 0 and Γ > 0 , which is the most relevant for the available materials and which remains enigmatic in both quantum and classical limits, despite much effort. We employ large-scale Monte Carlo simulations on specially designed finite-size clusters and unravel the presence of a complex multisublattice magnetic order in a wide region of the phase diagram, whose structure is characterized in detail. We show that this order can be quantified in terms of a coarse-grained scalar-chirality order, featuring a counterrotating modulation on the two spin sublattices. Here, we also provide a comparison to previous studies and discuss the impact of quantum fluctuations on the phase diagram.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Computational discovery of a tetragonal MBene phase: diamond-shaped boron lattices, viable MAB precursors, and selective CO 2 reduction

MBenes are a class of two-dimensional transition-metal borides derived from layered MAB phases. Here, we use first-principles calculations to predict a new tetragonal MB (1 : 1) MBene phase featuring a diamond-shaped B–B lattice and systematically assess its phase stability, synthetic accessibility, and catalytic performance towards the CO 2 reduction reaction (CO 2 RR). Four tetragonal configurations are benchmarked against known hexagonal and orthorhombic phases, and a multi-tier screening identifies six robust members (CrB, FeB, MoB, WB, IrB, and PtB). Exploration of the I4/mmm MAB space reveals 11 viable precursors, including the experimentally synthesized Ir 2 ZnB 2 , which validates our screening approach. CO 2 RR free-energy diagrams show that CrB favors CH 3 OH while FeB, MoB, WB, IrB, and PtB preferentially yield HCOOH. The corresponding limiting potentials are competitive with, and in some cases superior to, those reported for state-of-the-art orthorhombic and hexagonal MBenes. These findings expand the MBene landscape by establishing a tetragonal phase with structural, synthetic, and catalytic promise.

Lu, Linguo [University of Puerto Rico, San Juan, P↗
DOE OSTI2025

Quasicrystal stability and nucleation kinetics from density functional theory

The aperiodic order of quasicrystals bridges the amorphous and crystalline regime, so it has remained unclear whether quasicrystals are metastable or stable phases of matter. Density functional theory is often used to evaluate thermodynamic stability, but quasicrystals are long-range aperiodic and their energies cannot be calculated using conventional ab initio methods. Here, in this work, we perform first-principles calculations on quasicrystal nanoparticles of increasing size, from which we can directly extrapolate their bulk and surface energies. Using this technique, we determine with high confidence that the icosahedral quasicrystals ScZn 7.33 and YbCd 5.7 are ground-state phases, thus revealing that translational symmetry is not a necessary condition for the zero-temperature stability of inorganic solids. Although we found the ScZn 7.33 quasicrystal to be thermodynamically stable, we show on a mixed thermodynamic and kinetic phase diagram that its solidification from the melt is limited by nucleation, which illustrates why even stable materials may be kinetically challenging to grow. Our techniques broadly open the door to first-principles investigations into the structure–bonding–stability relationships of aperiodic materials.

density functional theory↗
NASA NTRS2025

Mechanics of Preloaded Bolt Tensile Loading With Focus on Load Introduction Factor

The bolt tensile and joint separation loads are directly influenced by the locations at which the external loads enter the clamped members of a preloaded bolted joint (PBJ) and the associated load-paths through the joint. This physical load introduction mechanism affecting the bolt tensile loading is typically represented in the bolt tensile load equation, in part, by a load introduction factor (LIF), which was shown by H.M. Lee of Marshall Spaceflight Center to be a natural product of the bolt tensile load equation using a linear spring stiffness model. This LIF, being a function of load-path stiffness, has subsequently been denoted as the stiffness-based LIF (SBLIF), providing a framework to calculate the LIF using whatever load-path stiffness approximations are appropriate. Expanding upon the work of Lee, it is shown that the SBLIF and the joint stiffness factor are functions of the stiffnesses of the same load-paths and regions within a PBJ, and thus they should not be treated as independent variables. Mathematical expressions for the SBLIF are presented. Comparisons are shown between the analytically calculated SBLIF, the analytically calculated geometric LIF (GLIF), which is a simple clamped-member thickness ratio, the experimentally derived LIF, and the LIF determined by finite element analysis (FEA). Using experiment and FEA as a benchmark, the SBLIF, using traditional load-path stiffness approximations, enables a more accurate prediction of bolt tensile loading than the GLIF, although it can be unconservative near joint separation. The GLIF generally attributes more of the externally applied tensile load to the bolt than does the SBLIF, potentially resulting in heavier and/or more costly bolted joints. Mathematical relationships between the SBLIF and the GLIF are developed. Supplemental material is provided in the appendixes where the historical practice of using the joint compressive stiffness in place of the joint tensile stiffness is evaluated. The appendixes include step-by-step examples demonstrating the calculation of the SBLIF using traditional stiffness approximations and conclude with the development of the joint diagram in terms of the SBLIF, culminating into formulas for the key features of a joint diagram, which is useful for programming.

Load Path↗
DOE OSTI2026

Laser-heated diamond anvil cell synthesis and recovery of metastable MnSb 2 and YbZn 2 for post-synthesis transport studies

The creation and exploration of new materials under extreme pressure–temperature conditions has become increasingly reliant on laser-heated diamond anvil cell (LHDAC) techniques, which provide direct access to previously unexplored regions of multinary phase diagrams. Whereas numerous high-pressure phases have been identified in situ, systematic recovery and post-synthesis physical property characterization of these materials remain significant challenges. In this work, we describe the setup and implementation of an LHDAC-based synthesis and recovery workflow and demonstrate its application to metastable MnSb 2 and YbZn 2 phases. Synchrotron x-ray diffraction and spatial mapping confirm dominant formation of the targeted phases, whereas laboratory-based refinement quantifies phase fractions despite intrinsic microstrain and minor secondary phases. High-pressure transport measurements on recovered samples reveal pressure-tunable electronic instabilities in both systems. In MnSb 2 , pressure suppresses two high-temperature magnetic ordering anomalies, observed in transport, by ∼5 GPa and, for higher pressures, induces a new low-temperature feature that increases with further pressure increase. In hexagonal high-pressure YbZn 2 , an electronic reconstruction emerges at ∼11 GPa, characterized by semiconducting-like behavior from ∼30 to 300 K and a broad low-temperature coherence crossover near 30 K. Our results establish LHDAC synthesis not only as a structural discovery tool but also as an experimental platform for investigating correlated quantum states stabilized far from equilibrium thermodynamic conditions.

Huyan, S. [Iowa State University, Ames, IA (United↗
NASA NTRS2020

Evaluation of Fatigue Damage Accumulation Functions for Delamination Initiation and Propagation

The present report follows on the cohesive fatigue damage model methodology proposed in NASA-TP-2018-219838. In that report, an empirical function describing the incremental damage due to cyclic loading was used to calculate fatigue damage within a cohesive formulation. The form of the function was developed such that, when integrated at a constant stress amplitude from no damage to failure, it produces a life versus load response that is consistent with an S-N diagram. Therefore, the parameters of the damage model could be obtained by fitting the model predictions to an S-N diagram. The finite element analyses performed demonstrate that the cohesive fatigue accumulation function provides a link between the S-N diagram that describes crack initiation, and the Paris law that characterizes the rate of crack propagation. However, when the model was proposed, it was not known whether the form of the damage accumulation function associated with a desired S-N diagram is unique and, if not, if the link between S-N and the Paris law is unique and independent of the fatigue function selected. In the effort described herein, several alternative forms of the damage function that reproduce the desired features of S-N diagrams were found and evaluated. The effects of each of these functions on the predicted parameters of the Paris law and the propagation threshold are discussed. The results indicate that the predicted exponent m of the Paris law is indeed independent of the damage accumulation function. However, different functions predict different values for the pre-factor C of the Paris law. Therefore, the proper damage accumulation function must be selected by comparison with experiments. One of the new damage accumulation functions proposed herein was found to be particularly useful for analysis because of the ease with which the model parameters can be determined with a minimal amount of experimental information. The effectiveness of the proposed methodology and damage function was demonstrated by conducting analyses of a double cantilever beam test, a mixed-mode bending test, and a three-point bending test of a skin/doubler specimen. The results indicate that the same set of model parameters can provide accurate predictions of the rate of fatigue crack propagation for a variety of material interfaces, mode mixities, load levels, and stress ratios.

DKIN/Stiffener Debonding↗
DOE OSTI2025

CASM Monte Carlo: Calculations of the thermodynamic and kinetic properties of complex multicomponent crystals

Monte Carlo techniques play a central role in statistical mechanics approaches that connect macroscopic thermodynamic and kinetic properties to the electronic structure of a material. This paper describes the implementation of Monte Carlo techniques for the study of multicomponent crystalline materials within the Clusters Approach to Statistical Mechanics (CASM) software suite, and demonstrates their use in model systems to calculate free energies and kinetic coefficients, study phase transitions, and construct phase diagrams from first principles. Many crystal structures are complex, with multiple sublattices occupied by differing sets of chemical species, along with the presence of vacancies or interstitial species. This imposes constraints on concentration variables, the form of thermodynamic potentials, and the values of kinetic transport coefficients. The framework used by CASM to formulate thermodynamic potentials and kinetic transport coefficients accounting for arbitrarily complex crystal structures is presented and demonstrated with examples of increasing complexity. Additionally, an overview of the capabilities of the CASM software specific to Monte Carlo methods is given, and a new CASM software package is introduced, casm-flow, which helps automate the setup, submission, management, and analysis of Monte Carlo simulations.

Cluster expansion↗
NASA NTRS1978

ITOS Power System Design

A brief description of the ITOS power systems is presented. The descriptions are from viewgraphs shown to the participants of the Workshop. A block diagram of ITOS power supply subsystem is presented along with graphs representing: ITOS voltage limit versus temperature; and battery charge-discharge current profiles over an orbit period.

A F Obenschain↗
NASA NTRS1993

Hubble Space Telescope nickel hydrogen battery system briefing

The topics covered are presented in viewgraph form and include the following: the Hubble Space Telescope (HST) Mission; system constraints; battery specification; battery module; simplified block diagram; cell design summary; present status; voltage decay; system depth of discharge; pressure since launch; system capacity; eclipse time vs. trickle charge; capacity test objectives; and capacity during tests.

Nawrocki, David↗