Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equation of State”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

BCS-BEC crossover of the strongly interacting 6 Li − 40 K mixture

Here, we present quantum Monte Carlo calculations of the properties of a two-component mass-imbalanced Fermi gas, corresponding to the 6 Li − 40 K mixture. We compute the equation of state of the unpolarized system as a function of the scattering length with particular attention paid to the unitary limit, where the effect of the effective range of the interaction is explored. In order to investigate differences from the equal-mass case we also compute the pair-distribution function and the momentum distribution over a range of interaction strengths, which can provide information about the structure of the system. Finally, we compute the heavy-light quasiparticle spectrum for the full crossover regime. Our theoretical predictions, based on quantum Monte Carlo calculations, should inform future theoretical studies and can be tested by experiments with ultracold fermionic gases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Two-component atomic Fermi superfluid with spin-orbital coupling in thin-spherical-shell geometry

We present a theory of two-component atomic Fermi superfluid with tunable pairing interaction in a thin spherical shell subject to spin-orbit coupling (SOC). By incorporating SOC into the Fermi superfluid in the BCS-Bose Einstein condensation (BEC) crossover, we obtain the energy spectrum and equations of state. While the order parameter and chemical potential are suppressed by SOC on the BCS side, the former remains positive but the latter may be pushed to negative values by SOC.Meanwhile, the compressibility exhibits kinks as the pairing interaction or SOC varies, indicating singular behavior of higher derivatives of the free energy despite the presence of the superfluid order parameter. The minimum of the energy dispersion indicates a decoupling of the energy gap from the order parameter, and the compressibility kinks occur when the energy gap approaches 0. We found the particle-hole mixing in the Fermi superfluid essential for the compressibility kinks since a Fermi gas with the same SOC but no pairing interaction only exhibits complicated dispersions but no singular behavior. Therefore, our results show that a combination of superfluid, SOC, and geometry can lead to interesting phenomena. Here, we also discuss possible experimental realizations and implications.

Fermi gases↗
DOE OSTI2020

Path integral Monte Carlo and density functional molecular dynamics simulations of warm dense MgSiO 3

In order to provide a comprehensive theoretical description of MgSiO 3 at extreme conditions, we combine results from path integral Monte Carlo and density functional molecular dynamics simulations and generate a consistent equation of state for this material. We consider a wide range of temperature and density conditions from 10 4 to 10 8 K and from 0.321 to 64.2 g/cm -3 (0.1- to 20-fold the ambient density). We study how the L and K shell electrons are ionized with increasing temperature and pressure. We derive the shock Hugoniot curve and compare with experimental results. Our Hugoniot curve is in good agreement with the experiments, and we predict a broad compression maximum that is dominated by the K shell ionization of all three nuclei while the peak compression ratio of 4.70 is obtained when the Si and Mg nuclei are ionized. To conclude, we analyze the heat capacity and structural properties of the liquid.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Measurement of the sound velocity and Grüneisen parameter of polystyrene at inertial confinement fusion conditions

The principal Hugoniot, sound velocity, and Grüneisen parameter of polystyrene were measured at conditions relevant to shocks in inertial confinement fusion implosions, from 100 to 1000 GPa. The sound velocity is in good agreement with quantum molecular dynamics calculations and all tabular equation of state models at pressures below 200 GPa. Above 200 GPa, the experimental results agree with two of the examined tables, but do not agree with the most recent table developed for design of ICF experiments. The Grüneisen parameter increases with density below ~3.1 g/cm 3 and approaches the asymptotic value for an ideal gas after complete dissociation. This behavior is in good agreement with quantum molecular dynamics results and previous work but is not represented by any of the tabular models. Finally, the discrepancy between tabular models and experimental measurement of the sound velocity and Grüneisen parameter is sufficient to impact simulations of ICF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Electronic energy gap closure and metal-insulator transition in dense liquid hydrogen

Here, using quantum Monte Carlo (QMC) calculations, we investigate the insulator-metal transition observed in liquid hydrogen at high pressure. Below the critical temperature of the transition from the molecular to the atomic liquid, the fundamental electronic gap closure occurs abruptly, with a small discontinuity reflecting the weak first-order transition in the thermodynamic equation of state. Above the critical temperature, molecular dissociation sets in while the gap is still open. When the gap closes, the decay of the off-diagonal reduced density matrix shows that the liquid enters a gapless, but localized, phase: there is a crossover between the insulating and the metallic liquids. Compared to different density functional theory (DFT) functionals, our QMC calculations provide larger values for the fundamental gap and the electronic density of states close to the band edges, indicating that optical properties from DFT potentially benefit from error cancellations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Deep learning of accurate force field of ferroelectric HfO 2

The discovery of ferroelectricity in HfO 2 -based thin films opens up new opportunities for using this silicon-compatible ferroelectric to realize low-power logic circuits and high-density nonvolatile memories. The functional performances of ferroelectrics are intimately related to their dynamic responses to external stimuli such as electric fields at finite temperatures. Molecular dynamics is an ideal technique for investigating dynamical processes on large length and time scales, though its applications to new materials are often hindered by the limited availability and accuracy of classical force fields. Here we present a deep neural network–based interatomic force field of HfO 2 learned from ab initio data using a concurrent learning procedure. The model potential is able to predict structural properties such as elastic constants, equation of states, phonon dispersion relationships, and phase transition barriers of various hafnia polymorphs with accuracy comparable with density functional theory calculations. Here, the validity of this model potential is further confirmed by the reproduction of experimental sequences of temperature-driven ferroelectric-paraelectric phase transitions of HfO 2 with isobaric-isothermal ensemble molecular dynamics simulations. We suggest a general approach to extend the model potential of HfO 2 to related material systems including dopants and defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

High-pressure synthesis and thermodynamic stability of PdH 1±ε up to 8 GPa

Palladium hydride alloys are superconductors and hydrogen storage materials. One synthesis route is compression of Pd to high pressure in a hydrogen-rich environment. Here we report the evolution of the unit cell volume of PdH x synthesized by compressing Pd in a pure H 2 medium to pressures from 0.2 to 8 GPa in a diamond anvil cell at room temperature. The volume of the face-centered cubic unit cell changes nonmonotonically with pressure, increasing upon compression from 0.2 to 1 GPa and decreasing upon compression from 1 to 8 GPa. Volume is reversible upon decompression and is independent of whether the sample was heated to 600 K at low pressure (P < 2 GPa). The x-ray diffraction data show no evidence for a phase transition between 0.2 and 8 GPa. The volume maximum at 1 GPa must be caused by progressive hydrogenation from 0 to 1 GPa. Assuming a pressure-volume-composition equation of state derived from previously published data, the [H]:[Pd] ratio in this study increases to a maximum value of x = 1 ± 0.02 at 2 ± 0.5 GPa and remains stable upon further compression to and from 8 GPa. These results add to a mounting body of evidence that PdH 1±ϵ is in thermodynamic equilibrium with pure H 2 at room temperature from 2 GPa to at least 8 GPa. Here, the simplest interpretation is that H atoms occupy all octahedral sites and no tetrahedral sites in face-centered cubic PdH 1.0 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Shock-compressed silicon: Hugoniot and sound speed up to 2100 GPa

In this work, the high-pressure equation of state and isentropic sound speed of fluid silicon were studied using laser-driven shock waves. Principal Hugoniot measurements to 2100 GPa were performed using impedance matching to an α-quartz reference. Sound speed was determined by time correlating imposed shock-velocity perturbations in both the sample (Si) and reference material (α-quartz). A change in shock velocity versus particle velocity (u s –u p ) slope on the fluid silicon principal Hugoniot was detected at 200 GPa. Density functional theory-based quantum molecular dynamics simulations suggest that an increase in ionic coordination and 50% increase in average ionization is coincident with the experimentally observed change in slope.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Ab initio anharmonic thermodynamic properties of cubic Ca Si O 3 perovskite

We present an ab initio study of the thermodynamic properties of cubic CaSiO 3 perovskite (CaPv) over the pressure and temperature range of the Earth's lower mantle. We compute the anharmonic phonon dispersions throughout the Brillouin zone by utilizing the phonon quasiparticle approach, which characterizes the intrinsic temperature dependence of phonon frequencies and, in principle, captures full anharmonicity. Such temperature-dependent phonon dispersions are used to calculate ab initio free energy in the thermodynamic limit (N → ∞) within the framework of the phonon gas model. Accurate free energy calculations enable us to investigate cubic CaPv's thermodynamic properties, e.g., thermal expansivity, Grüneisen parameter, bulk modulus, heat capacity and thermal equation of state, where anharmonic effects are demonstrated. Furthermore, the present methodology provides an important theoretical approach for exploring phase boundaries, thermodynamic, and thermoelastic properties of strongly anharmonic materials at high pressures and temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Thermodynamic properties of ε -Fe with thermal electronic excitation effects on vibrational spectra

The thermodynamic properties of hexagonal-close-packed iron (ε–Fe) are essential for investigating the internal structure and dynamic properties of planetary cores. Despite their importance to planetary sciences, experimental investigations of ε–Fe at relevant conditions are still challenging. Therefore, ab initio calculations are crucial to elucidating the thermodynamic properties of this system. Here, we use a free energy calculation scheme based on the phonon gas model compatible with temperature-dependent phonon frequencies. We investigate the effects of electronic thermal excitations, which introduces a temperature dependence on phonon frequencies, and the implication for the thermodynamic properties of ε–Fe at extreme pressure (P) and temperature (T) conditions. We disregard phonon-phonon interactions, i.e., anharmonicity and their effect on phonon frequencies. Nevertheless, the current scheme is also applicable to T -dependent anharmonic frequencies. We conclude that the impact of thermal electronic excitations on vibrational properties is not significant up to ~4000 K at 200 GPa but should not be ignored at higher temperatures or pressures. However, the static free energy F st must always include the effect of thermal excitation fully in a continuum of T. Furthermore, our results for isentropic equations of state show good agreement with data from recent ramp compression experiments up to 1400 GPa conducted at the National Ignition Facility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

In situ synchrotron diffraction of pressure-induced phase transition in DyPO 4 under variable hydrostaticity

In situ synchrotron x-ray diffraction was conducted on polycrystalline DyPO 4 to elucidate the details of the pressure-induced transition from the xenotime polymorph to the monazite polymorph. We used three different pressure-transmitting media (neon, a 16:3:1 methanol-ethanol-water mixture, and potassium chloride) to investigate the effect of hydrostaticity on the phase behavior. Specifically, our data clearly show a hydrostatic onset pressure of the xenotime-monazite transition of 9.1 GPa, considerably lower than the 15.3 GPa previously determined by Raman spectroscopy. Based on (quasi)hydrostatic data taken in a neon environment, third-order Birch-Murnaghan equation-of-state fits give a xenotime bulk modulus of 144 GPa and a monazite bulk modulus of 180 GPa (both with pressure derivatives of 4.0). Structural data and axial compressibilities show that DyPO 4 is sensitive to shear and has an anisotropic response to pressure. More highly deviatoric conditions cause the onset of the transition to shift to pressures at least as low as 7.0 GPa. We attribute early transition to shear-induced distortion of the PO4 tetrahedra. Finally, our characterization of the high-pressure behavior of DyPO 4 under variable hydrostaticity is critical for advancing rare earth orthophosphate fiber coating applications in ceramic matrix composites and may inform future tailoring of phase composition for controlled shear and pressure applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Shock melting and the hcp-bcc phase boundary of Mg under dynamic loading

The high-temperature and high-pressure response of magnesium has been investigated through shock-release experiments performed up to shock melting. The longitudinal and bulk sound speeds of Mg are reported along the Hugoniot from 25 to 56 GPa and used to determine the elastic properties and Grüneisen parameter. The hexagonal close-packed (hcp)–body-centered cubic (bcc) phase transition is marked by a reduction in the determined shear wave speed. Thermal softening is observed to begin between 40 and 44 GPa, with incipient melt at 55.5 GPa, in close agreement with previous diffraction measurements under laser shock loading. Examination of the release profiles showed that two different responses were observed, depending on the peak stress. When plotted in a stress-energy phase diagram, the two responses are observed to form separate lines that intersect the static hcp-bcc phase transition and incipient melt. The results indicate that the hcp-bcc transition occurs on the Hugoniot at 28.4 GPa and place the hcp-bcc-liquid triple point ~20 GPa. In conclusion, a multiphase equation of state is developed which places the melt boundary below the previously reported static measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atom-in-jellium predictions of the shear modulus at high pressure

Atom-in-jellium calculations of the Einstein frequency in condensed matter and of the equation of state were used to predict the variation of shear modulus from zero pressure to ~10 7 g/cm 3 , for several elements relevant to white dwarf stars and other self-gravitating systems. This is by far the widest range reported electronic structure calculation of shear modulus, spanning from ambient through the one-component plasma to extreme relativistic conditions. The predictions were based on a relationship between the Debye temperature and shear modulus, which we assess to be accurate at the o(10%) level, and is the first known use of atom-in-jellium theory to calculate a shear modulus. Finally, we assessed the overall accuracy of the method by comparing with experimental measurements and more detailed electronic structure calculations at lower pressures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Phase transitions in high-purity zirconium under dynamic compression

In this work, we present results from ramp compression experiments on high-purity Zr that show the α → ω, ω → β, as well as reverse β → ω phase transitions. Simulations with a multi-phase equation of state and phenomenological kinetic model match the experimental wave profiles well. While the dynamic α ! ! transition occurs ~ 9 GPa above the equilibrium phase boundary, the ω → β transition occurs within 0.9 GPa of equilibrium. We estimate that the dynamic compression path intersects the equilibrium ω – β line at P = 29:2 GPa, and T = 490 K. The thermodynamic path in the interior of the sample lies ~ 100 K above the isentrope at the point of the ω → β transition. Approximately half of this dissipative temperature rise is due to plastic work, and half is due to the non-equilibrium α → ω transition. The inferred rate of the α → ω transition is several orders of magnitude higher than that measured in dynamic diamond anvil cell (DDAC) experiments in an overlapping pressure range. We discuss a model for the influence of shear stress on the nucleation rate. We find that the shear stress s ji has the same effect on the nucleation rate as a pressure increase δP = cϵ ij s ji /(ΔV/V); where c is a geometric constant ~ 1 and, ϵ ij are the transformation shear strains. The small fractional volume change ΔV/V ≈ 0:1 at the α → ω transition amplifies the effect of shear stress, and we estimate that for this case δP is in the range of several GPa. Correcting our transition rate to a hydrostatic rate brings it approximately into line with the DDAC results, suggesting that shear stress plays a significant role in the transformation rate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Hugoniot, sound speed, and phase transitions of single-crystal sapphire for pressures 0.2–2.1 TPa

Sapphire (Al 2 O 3 ) is a major constituent of the Earth's mantle and has significant contributions to the field of high-pressure physics. Constraining its Hugoniot over a wide pressure range and identifying the location of shock-driven phase transitions allows for development of a multiphase equation of state and enables its use as an impedance-matching standard in shock physics experiments. In this paper we present measurements of the principal Hugoniot and sound velocity from direct impact experiments using magnetically launched flyers on the Z machine at Sandia National Laboratories. The Hugoniot was constrained for pressures from 0.2–2.1 TPa and a four-segment piecewise linear shock-velocity–particle-velocity fit was determined. First-principles molecular dynamics simulations were conducted and agree well with the experimental Hugoniot. Sound-speed measurements identified the onset of melt between 450 and 530 GPa, and the Hugoniot fit refined the onset to 525 ± 13 GPa. A phase diagram which incorporates literature diamond-anvil cell data and melting measurements is presented.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Elastic properties of the W 0.75 R e 0.25 alloy at high pressure up to 183 GPa

The high pressure equation of state for the W 0.75 R e 0.25 alloy is experimentally determined up to 183 GPa with synchrotron angle-dispersive powder x-ray diffraction in the diamond-anvil cell and to ∼925 GPa with density-functional theory. W-Re alloys are used in many industrial high-temperature applications and as a confining gasket material in high-pressure diamond-anvil cell research. The inclusion of 25 wt. % Re achieves the highest performance in terms of strength and ductility while also maintaining the body-centered-cubic (bcc) crystal structure, yet to date there has been no investigation into its elastic behavior at high pressure. We present the experimentally and theoretically determined volumetric and elastic pressure response and systematically compare these results to other W-Re alloys, finding that the bulk modulus of W-Re alloys varies nonlinearly with Re content and W 0.75 Re 0.25 becomes less incompressible than W at 85 GPa. Published by the American Physical Society 2025

Alloys↗
DOE OSTI2020

Neutrino signal from proto-neutron star evolution: Effects of opacities from charged-current–neutrino interactions and inverse neutron decay

In this work, we investigate the impact of charged-current–neutrino processes on the formation and evolution of neutrino spectra during the deleptonization of proto-neutron stars. To this end we develop the full kinematics of these reaction rates consistent with the nuclear equation of state, including weak magnetism contributions. This allows us to systematically study the impact of inelastic contributions and weak magnetism on the $ν_e$ and $\bar{ν}_e$ luminosities and average energies. Furthermore, we explore the role of the inverse neutron decay, also known as the direct Urca process, on the emitted spectra of $\bar{ν}_e$. This process is commonly considered in the cooling scenario of cold neutron stars but has so far been neglected in the evolution of hot proto-neutron stars. We find that the inverse neutron decay becomes the dominating opacity source for low-energy $\bar{ν}_e$. Accurate three-flavor Boltzmann neutrino transport enables us to relate the magnitude of neutrino fluxes and spectra to details of the treatment of weak processes. This allows us to quantify the corresponding impact on the conditions relevant for the nucleosynthesis in the neutrino-driven wind, which is ejected from the proto-neutron star surface during the deleptonization phase.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Pions in hot dense matter and their astrophysical implications

Here we study the role of pions in hot dense matter encountered in astrophysics. We find that strong interactions enhance the number density of negatively charged pions and that this enhancement can be calculated reliably for a relevant range of density and temperature by using the virial expansion. We assess the influence of pions and muons on the equation of state (EOS) and weak-interaction rates in hot dense matter. We find that thermal pions increase the proton fraction and soften the EOS. We also find that charged current weak reactions involving pions and muons $\bar{ν}_ μ + μ^– → π$ and $ν_μ + π^– → μ^–$ make an important contribution to the opacity of muon neutrinos. This could influence the dynamics of core-collapse supernovae and neutron-star mergers. Finally, we note that pion-nucleon reactions can alter the evolution of the proton fraction when weak interactions are not in equilibrium.

79 ASTRONOMY AND ASTROPHYSICS↗