Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase diagram”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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

Phase diagram of magnetic shape memory alloy Ni 50 Mn $50–x$ In $x$ , 0 < $x$ , 25 from first principles, via spin cluster expansion and phonon vibrational entropies

The metamagnetic shape memory Heusler alloy Ni 50 Mn $50–x$ In $x$ exhibits a rich phase diagram featuring competing magnetic states, coupled magnetic–structural phase transitions, and strong compositional sensitivity. Existing first-principles approaches struggletocapturetheintertwinedchemical, magnetic, andvibrationaleffectsinthesealloys, necessitating a more integrated modeling framework. We develop a spin cluster expansion (spin-CE) framework augmented by a quasi-harmonic phonon model to capture both configurational (chemical and magnetic) and vibrational contributions to the free energy of Ni 50 Mn $50–x$ In $x$ over the full compositional range 0 ≤ x ≤25. The spin-CE includes both chemical clusters and composition-dependent Ising spin interactions, with parameters fit to a first-principles density functional theory (DFT) dataset. Using this approach, we predict the complete magnetostructural phase diagram and transformation temperatures of Ni 50 Mn $50–x$ In $x$ across the composition space. We find that vibrational entropy alone is insufficient to reproduce the martensitic transformation in the magnetic shape memory alloy regime, highlighting the essential role of magnetism. Incorporating both magnetic and vibrational contributions allows us to reproduce all experimentally known phases, including the disappearance of the stable martensite phase at a critical In concentration and the Curie temperature of the austenite phase. The method also captures the transition with increasing In in martensite from antiferromagnetic to ferromagnetic order and predicts re-entrant ferromagnetism, though the latter occurs at higher In content than reported experimentally. We discuss possible sources of this discrepancy and highlight the broader applicability of the method to other magnetostructurally complex materials, where it may offer mechanistic insight and predictive design capabilities.

Cluster expansion↗
DOE OSTI2023

First-principles molten salt phase diagrams through thermodynamic integration

Precise prediction of phase diagrams in molecular dynamics simulations is challenging due to the simultaneous need for long time and large length scales and accurate interatomic potentials. Here, we show that thermodynamic integration from low-cost force fields to neural network potentials trained using density-functional theory (DFT) enables rapid first-principles prediction of the solid–liquid phase boundary in the model salt NaCl. We use this technique to compare the accuracy of several DFT exchange–correlation functionals for predicting the NaCl phase boundary and find that the inclusion of dispersion interactions is critical to obtain good agreement with experiment. Importantly, our approach introduces a method to predict solid–liquid phase boundaries for any material at an ab initio level of accuracy, with the majority of the computational cost at the level of classical potentials.

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

Phase diagram of the TIP4P/Ice water model by enhanced sampling simulations

In this work, we studied the phase diagram for the TIP4P/Ice water model using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice–liquid free energy differences from biased coexistence simulations that reversibly sample the melting and growth of layers of ice. We computed a total of 19 melting points for five different ice polymorphs, which are in excellent agreement with the melting lines obtained from the integration of the Clausius–Clapeyron equation. For proton-ordered and fully proton-disordered ice phases, the results are in very good agreement with previous calculations based on thermodynamic integration. For the partially proton-disordered ice III, we find a large increase in stability that is in line with previous observations using direct coexistence simulations for the TIP4P/2005 model. This issue highlights the robustness of the approach employed here for ice polymorphs with diverse degrees of proton disorder. Our approach is general and can be applied to the calculation of other complex phase diagrams.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Multiferroicity and phase diagram of ferro-rotational magnet RbFe(SO 4 ) 2

Ferro-rotational magnet RbFe(SO 4 ) 2 has attracted attention for its stable ferro-rotational phase and electric-field-controlled magnetic chirality. This work presents the multiferroic properties and H–T phase diagram of RbFe(SO 4 ) 2 , which have been underexplored. Our measurements of magnetic susceptibility, ferroelectric polarization, and dielectric constant under various magnetic fields reveal four distinct phases: (I) a ferroelectric and helical magnetic phase below 4 K and 6 T, (II) a paraelectric and collinear magnetic phase below 4 K and above 6 T, (III) a paraelectric and non-collinear magnetic phase below 4 K and above 9 T, and (IV) a paraelectric and paramagnetic above 4 K. This study clarifies the multiferroic behavior and H–T phase diagram of RbFe(SO 4 ) 2 , providing valuable insights into ferro-rotational magnets.

ferro-rotation↗
DOE OSTI2022

Ab Initio Phase Diagram of Chromium to 2.5 TPa

Chromium possesses remarkable physical properties such as hardness and corrosion resistance. Chromium is also a very important geophysical material as it is assumed that lighter Cr isotopes were dissolved in the Earth’s molten core during the planet’s formation, which makes Cr one of the main constituents of the Earth’s core. Unfortunately, Cr has remained one of the least studied 3d transition metals. In a very recent combined experimental and theoretical study (Anzellini et al., Scientific Reports, 2022), the equation of state and melting curve of chromium were studied to 150 GPa, and it was determined that the ambient body-centered cubic (bcc) phase of crystalline Cr remains stable in the whole pressure range considered. However, the importance of the knowledge of the physical properties of Cr, specifically its phase diagram, necessitates further study of Cr to higher pressure. In this work, using a suite of ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries, we obtain the theoretical phase diagram of Cr to 2.5 TPa. We calculate the melting curves of the two solid phases that are present on its phase diagram, namely, the lower-pressure bcc and the higher-pressure hexagonal close-packed (hcp) ones, and obtain the equation for the bcc-hcp solid–solid phase transition boundary. We also obtain the thermal equations of state of both bcc-Cr and hcp-Cr, which are in excellent agreement with both experimental data and QMD simulations. We argue that 2180 K as the value of the ambient melting point of Cr which is offered by several public web resources (“Wikipedia,” “WebElements,” “It’s Elemental,” etc.) is most likely incorrect and should be replaced with 2135 K, found in most experimental studies as well as in the present theoretical work.

equation of state↗
DOE OSTI2023

Activity coefficient acquisition with thermodynamics–informed active learning for phase diagram construction

This work explores the use of thermodynamics-informed Gaussian processes (GPs) and active learning (AL) to model activity coefficients and construct phase diagrams. Relying on synthetic data generated from an excess Gibbs energy model, GPs were found to accurately describe the activity coefficients of several binary mixtures across large composition and temperature ranges. Moreover, GPs could estimate their own uncertainty and identify composition/temperature regions where activity coefficient data provide the most information to the models. This was leveraged to build AL algorithms targeted at modeling phase equilibria. In many cases, a single active-learning-acquired data point was sufficient to describe the phase diagrams studied. Lastly, the ability of AL to greatly reduce the amount of data needed to obtain accurate models was further verified on experimental case studies, namely individual ion activity coefficients, the solid–liquid and vapor–liquid equilibrium of deep eutectic solvents, and phase equilibria in ternary mixtures.

25 ENERGY STORAGE↗
DOE OSTI2022

Complexity Phase Diagram for Interacting and Long-Range Bosonic Hamiltonians

Here, we classify phases of a bosonic lattice model based on the computational complexity of classically simulating the system. We show that the system transitions from being classically simulable to classically hard to simulate as it evolves in time, extending previous results to include on-site number-conserving interactions and long-range hopping. Specifically, we construct a complexity phase diagram with easy and hard “phases” and derive analytic bounds on the location of the phase boundary with respect to the evolution time and the degree of locality. We find that the location of the phase transition is intimately related to upper bounds on the spread of quantum correlations and protocols to transfer quantum information. Remarkably, although the location of the transition point is unchanged by on-site interactions, the nature of the transition point does change. Specifically, we find that there are two kinds of transitions, sharp and coarse, broadly corresponding to interacting and noninteracting bosons, respectively. Our Letter motivates future studies of complexity in many-body systems and its interplay with the associated physical phenomena.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2022

Essential role of magnetic frustration in the phase diagrams of doped cobaltites

Doped perovskite cobaltites (e.g., La 1–x Sr x CoO 3 ) have been extensively studied for their spin-state physics, electronic inhomogeneity, and insulator-metal transitions. Ferromagnetically interacting spin-state polarons emerge at low x in the phase diagram of these compounds, eventually yielding long-range ferromagnetism. The onset of long-range ferromagnetism (x≈0.18) is substantially delayed relative to polaron percolation (x≈0.05), however, generating a troubling inconsistency. Here, in this study, Monte Carlo simulations of a disordered classical spin model are used to establish that previously ignored magnetic frustration is responsible for this effect, enabling faithful reproduction of the magnetic phase diagram.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exploring the phase diagram of 3D artificial spin-ice

Artificial spin-ices consist of lithographic arrays of single-domain magnetic nanowires organised into frustrated lattices. These geometries are usually two-dimensional, allowing a direct exploration of physics associated with frustration, topology and emergence. Recently, three-dimensional geometries have been realised, in which transport of emergent monopoles can be directly visualised upon the surface. Here we carry out an exploration of the three-dimensional artificial spin-ice phase diagram, whereby dipoles are placed within a diamond-bond lattice geometry. We find a rich phase diagram, consisting of a double-charged monopole crystal, a single-charged monopole crystal and conventional spin-ice with pinch points associated with a Coulomb phase. In experimental demagnetised systems, broken symmetry forces formation of ferromagnetic stripes upon the surface, forbidding the lower energy double-charged monopole crystal. Instead, we observe crystallites of single magnetic charge, superimposed upon an ice background. The crystallites are found to form due to the distribution of magnetic charge around the 3D vertex, which locally favours monopole formation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Hot QCD phase diagram from holographic Einstein–Maxwell–Dilaton models

In this review, we provide an up-to-date account of quantitative holographic descriptions of the strongly coupled quark-gluon plasma (QGP) produced in heavy-ion collisions, based on the class of gauge-gravity Einstein-Maxwell-Dilaton (EMD) models. Holography is employed to tentatively map the QCD phase diagram at finite temperature onto a dual theory of charged, asymptotically AdS black holes in 5D. With a quantitative focus on the hot QCD phase diagram, the EMD models reviewed are adjusted to describe lattice results for the finite-temperature QCD equation of state, with 2+1 flavors and physical quark masses, at zero chemical potential and vanishing electromagnetic fields. The predictive power of EMD models is tested by quantitatively comparing their predictions for the hot QCD equation of state at nonzero baryon density and the corresponding state-of-the-art lattice QCD results. The shear and bulk viscosities predicted by these EMD models are also compared to the corresponding profiles favored by the latest phenomenological multistage models describing different heavy-ion data. We report preliminary results from a Bayesian analysis which provide systematic evidence that lattice results at finite temperature and zero baryon density strongly constrains the free parameters of EMD models. Remarkably, the set of parameters constrained by lattice results at zero chemical potential produces EMD models in quantitative agreement with lattice QCD results also at finite baryon density. We also review results for equilibrium and transport properties from magnetic EMD models, describing the QGP at finite temperatures and magnetic fields. Lastly, we provide a critical assessment of the main limitations and drawbacks of the holographic models reviewed in the present work, and point out some perspectives we believe are of fundamental importance for future developments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Exotic Kondo Phases: the non-Kramers Doniach phase diagram

This is the final report for the early career award “Exotic Kondo Phases: the non-Kramers Doniach phase diagram”. This project searched for exotic novel phases in heavy fermion materials and explored the possible phases of the two-channel Kondo model. Significant results include an enhanced understanding of the channel-symmetry breaking heavy Fermi liquid and its spinorial nature; a novel platform for realizing Majorana defects in topological defects in two-channel Kondo insulators; and a new non-equilibrium way to realize multi-channel Kondo models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Thermodynamic insights into the intricate magnetic phase diagram of EuAl 4

The tetragonal intermetallic compound EuAl 4 hosts an exciting variety of low-temperature phases. In addition to a charge density wave below 140 K, four ordered magnetic phases are observed below 15.4 K. Recently, a skyrmion phase was proposed based on Hall effect measurements under a c-axis magnetic field. Here, we present a detailed investigation of the phase transitions in EuAl 4 under c-axis magnetic field. Our dilatometry, heat-capacity, DC magnetometry, AC magnetic susceptibility, and resonant ultrasound spectroscopy measurements reveal three magnetic phase transitions not previously reported. The first key result is a detailed H ∥ [001] magnetic phase diagram mapping the seven phases we observe. Second, we identify a high-field phase, phase VII, which directly corresponds to the region were skyrmions have been previously suggested. Our results provide guidance for future studies exploring the complex magnetic interactions and spin structures in EuAl 4 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Determination of Site Occupancy in the M–Pd–Zn (M = Cu, Ag, and Au) γ-Brass Phase by CALculation of PHAse Diagrams Modeling and Rietveld Refinement

The Pd–Zn γ-brass phase provides exciting opportunities for synthesizing site-isolated catalysts with precisely controlled Pd active site ensembles. Introducing a third metallic element into the γ-brass lattice further perturbs the catalytic active site ensembles. Here, in this work, we introduce coinage metallic elements M (M = Cu, Ag, and Au) into the Pd–Zn γ-brass phase and investigate the site occupation factors of each element in the γ-brass lattice. The CALculation of PHAse Diagrams (CALPHAD) modeling approach supported by energetics predicted by the density functional theory and X-ray and neutron diffraction with Rietveld refinement were used to identify the SOF on each Wyckoff site for various M amounts alloyed into the Pd–Zn γ-brass phase. The present analysis unveils the strong preference for Pd occupying the outer tetrahedral (OT) site in the γ-brass lattice, while the coinage metallic elements tend to substitute for Zn on the octahedral (OH) site. The determination of site occupancy in the bulk M–Pd–Zn γ-brass phase provides opportunities to investigate and tailor potential catalytically active site ensembles in the γ-brass phase materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Phase Diagrams and Piezoelectric Properties of Wurtzite Al1-x-yScxGdyN Heterostructural Alloys

Ternary nitride alloys based on wurtzite AlN are a promising platform to realize functional materials, particularly ferroelectrics and optical emitters, that can smoothly integrate with conventional microelectronics. Here, a strategic design is presented to enable multifunctional materials by substituting multiple elements into AlN to create quaternary nitride alloys. By combining computational predictions and combinatorial thin film synthesis, the phase diagram of these quaternary Al-Sc-Gd-N alloys (or pseudo-ternary heterostructural AlN-ScN-GdN alloys) is successfully predicted as a function of effective temperature, and we experimentally grow Al1-x-xScxGdyN thin films for the first time. It is revealed that Al1-x-xScxGdyN crystallizes in a wurtzite-derived structure for X + y <~ 0.35, consistent with the calculated phase diagram. The computational investigation explores whether co-substitution induces cooperative effects on these alloys' piezoelectric and ferroelectric properties, finding that it is beneficial for reducing the polarization switching barrier. We calculate that Al1-x-xScxGdyN thin films should display ferroelectric switching. This is supported by our experimental measurements of a high optical bandgap, enhanced piezoelectric coefficient, and a change in the calculated polarization switching mechanism, and we achieve preliminary ferroelectric switching that experimentally realizes the prediction. Overall, our work sets the foundation toward quaternary wurtzite-nitride-based multifunctional materials, including piezoelectrics, ferroelectrics, and possibly even multiferroics.

36 MATERIALS SCIENCE↗