Search NASA⌕ Search

SEARCH · Search NASA

Results for “Liquid-solid phase transition”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

High-pressure phase diagram of beryllium from ab initio free-energy calculations

In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Strain release by 3D atomic misfit in fivefold twinned icosahedral nanoparticles with amorphization and dislocations

Multiple twinning to form fivefold twinned nanoparticles in crystal growth is common and has attracted broad attention ranging from crystallography research to physical chemistry and materials science. Lattice-misfit strain and defects in multiple twinned nanoparticles (MTP) are key to understand and tailor their electronic properties. However, the structural defects and related strain distributions in MTPs are poorly understood in three dimensions (3D). Here, we show the 3D atomic misfit and strain relief mechanism in fivefold twinned icosahedral nanoparticles with amorphization and dislocations by using atomic resolution electron tomography. We discover a two-sided heterogeneity in variety of structural characteristics. A nearly ideal crystallographic fivefold face is always found opposite to a less ordered face, forming Janus-like icosahedral nanoparticles with two distinct hemispheres. The disordered amorphous domains release a large amount of strain. Molecular dynamics simulations further reveal the Janus-like icosahedral nanoparticles are prevalent in the MTPs formed in liquid-solid phase transition. This work provides insights on the atomistic models for the modelling of formation mechanisms of fivefold twinned structures and computational simulations of lattice distortions and defects. We anticipate it will inspire future studies on fundamental problems such as twin boundary migration and kinetics of structures in 3D at atomic level.

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

Zero-point motion of liquid and solid hydrogen

We present an inelastic neutron scattering study of liquid and solid hydrogen carried out using the wide Angular Range Chopper Spectrometer at Oak Ridge National Laboratory. From the observed dynamic structure factor, we obtained empirical estimates of the molecular mean-squared displacement and average translational kinetic energy. Here we find that the former quantity increases with temperature, indicating that a combination of thermal and quantum effects is important near the liquid-solid phase transition, contrary to previous measurements. We also find that the kinetic energy drops dramatically on melting of the crystals, a consequence of the large increase in molar volume together with the Heisenberg indeterminacy principle. Our results are compared with quantum Monte Carlo simulations based on different model potentials. In general, there is good agreement between our findings and theoretical predictions based on the Silvera-Goldman and Buck potentials.

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 OSTI2021

Shock compression of silver to 300 GPa: Wave profile measurements and melting transition

Recent in-situ x-ray diffraction (XRD) measurements on laser-shock compressed Ag foils demonstrated a face-centered-cubic to body-centered-cubic transformation at ~150 GPa and melting between 172 and 197 GPa [Phys. Rev. Lett. 124, 235701 (2020)]. As a complement to the XRD work, we conducted plate impact experiments to obtain shock velocity and wave profile measurements on thicker Ag samples shock-compressed to peak stresses between 30 and 300 GPa. The shock velocity–particle velocity results were fitted very well by a linear relation over the entire stress range, providing an accurate determination of the Ag Hugoniot (locus of the stress-volume states achieved under shock compression). For peak stresses below 187 GPa and above 210 GPa—corresponding to the solid and liquid phases, respectively—the wave profiles show clean single waves. No wave profile features related to the fcc-bcc transformation at ~150 GPa were observed, implying minimal volume change for the transformation. For stresses between 187 and 210 GPa, an initial jump was followed by a time-dependent increase in the particle velocity (20–80 ns risetime) to the peak state—corresponding to the solid-liquid mixed phase response. Unlike the solid and liquid response, the mixed-phase response cannot be readily analyzed analytically. Instead, numerical simulations incorporating an accurate multiphase equation of state for Ag—not currently available—are required to analyze the wave profiles measured at 187–210 GPa stresses. Furthermore, the present work shows the potential for using wave profile measurements to examine the melting transition under shock compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Structural study of hcp and liquid iron under shock compression up to 275 GPa

Here we combine nanosecond laser shock compression with in situ picosecond x-ray diffraction to provide structural data on iron up to 275 GPa. We constrain the extent of hcp-liquid coexistence, the onset of total melt, and the structure within the liquid phase. Our results indicate that iron, under shock compression, melts completely by 258(8) GPa. A coordination number analysis indicates that iron is a simple liquid at these pressure-temperature conditions. We also perform texture analysis between the ambient body-centered-cubic (bcc) α, and the hexagonal-closed-packed (hcp) high-pressure ε-phase. We rule out the Rong-Dunlop orientation relationship (OR) between the α and ε - phase s. However, we cannot distinguish between three other closely related ORs: Burger's, Mao-Bassett-Takahashi, and Potter's OR. The solid-liquid coexistence region is constrained from a melt onset pressure of 225(3) GPa from previously published sound speed measurements and full melt [246.5(1.8)–258(8) GPa] from x-ray diffraction measurements, with an associated maximum latent heat of melting of 623 J/g. This value is lower than recently reported theoretical estimates and suggests that the contribution to the earth's geodynamo energy budget from heat release due to freezing of the inner core is smaller than previously thought. Melt pressures for these nanosecond shock experiments are consistent with gas gun shock experiments that last for microseconds, indicating that the melt transition occurs rapidly.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Mixed ortho- H 2 and para- H 2 clusters studied by vibrational coherent anti-Stokes Raman spectroscopy

The search for macroscopic quantum effects, including superfluidity, in molecular hydrogen is mostly focused on its parahydrogen (p-H 2 ) nuclear spin modification because of weaker intermolecular interaction compared to orthohydrogen (o-H 2 ), both modifications being bosonic. In this work, mixed clusters of o-H 2 and p-H 2 containing similar to 10(4) molecules are prepared by supersonic expansion with helium and studied by vibrational coherent anti-Stokes Raman scattering (CARS) spectroscopy. At similar experimental conditions the neat p-H 2 clusters avoid freezing and remain fluid at 1-2 K, which is predicted to be the realm of their superfluid behavior [Phys. Rev. Lett. 101, 205301 (2008)]. Dependence of the vibrational frequencies and intensities of the main CARS peaks due to o-H 2 and p-H 2 versus the ratio of the o-H 2 and p-H 2 concentrations in the expanding gas suggests that o-H 2 and p-H 2 molecules are uniformly mixed in the interior of the clusters. A weak spectral feature at 4157 cm -1 that appears independent of the concentration ratio is assigned to the outer shell of the clusters enriched with p-H 2 molecules. Although the phase of the mixed clusters could not be unambiguously identified, the shift of the vibrational frequencies with respect to the bulk solid is consistent with the liquid state of the clusters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Shock compression of vanadium at extremes: Theory and experiment

Here, the equation of state (EOS) and shock compression of bulk vanadium were investigated using canonical ab initio molecular dynamic simulations, with experimental validation to 865 GPa from shock data collected at Sandia's Z Pulsed Power Facility. In simulations the phase space was sampled along isotherms ranging from 3000 K to 50 000 K, for densities between ρ = 3 and 15 g / cm 3 , with a focus on the liquid regime and the body-centered-cubic phase in the vicinity of the melting limit. The principal Hugoniot predicted from first principles is overall consistent with shock data, while it showed that current multiphase SESAME-type EOS for vanadium needed revision in the liquid regime. A more accurate SESAME EOS was developed using constraints from experiments and simulations. This work emphasizes the need to use a combined theoretical and experimental approach to develop high-fidelity EOS models for extreme conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Sound speed measurements in silver shock compressed to 300 GPa: Solid-state transition, melting, and liquid-state response

To gain insight into the thermodynamic response of shock compressed Ag states corresponding to the face-centered-cubic (fcc) to body-centered-cubic (bcc) transformation and melting observed using in-situ x-ray diffraction (XRD) measurements, we present longitudinal sound speed results and their analysis in Ag at peak stresses ranging from 60 to 300 GPa. The measured sound speeds increased linearly with density compression to 171 GPa, showing that the sound speeds (and longitudinal moduli) in the fcc and bcc phases are very similar. Between 171 and 218 GPa, the sound speed dropped significantly, consistent with the melting reported using XRD measurements in shock compressed Ag. From 218 to 300 GPa, the increasing sound speeds and Hugoniot states provide a determination of the liquid phase Ag response. In particular, determination of the Grüneisen parameter (Γ) showed that the density – Grüneisen parameter product (ρΓ) for liquid Ag is constant, but differs significantly from that for solid Ag at ambient conditions. Furthermore, the Mie- Grüneisen equation of state can be used to describe the Hugoniot and off-Hugoniot response of liquid Ag.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Dynamical defects in a two-dimensional Wigner crystal: Self-doping and kinetic magnetism

We study the quantum dynamics of interstitials and vacancies in a two-dimensional Wigner crystal (WC) using a semi-classical instanton method that is asymptotically exact at low density, i.e., in the r s → ∞ limit. Here, the dynamics of these point defects mediates magnetism with much higher energy scales than the exchange energies of the pure WC. Via exact diagonalization of the derived effective Hamiltonians in the single-defect sectors, we find the dynamical corrections to the defect energies. The resulting expression for the interstitial (vacancy) energy extrapolates to 0 at r s = r mit ≈ 70 (r s ≈ 30), suggestive of a self-doping instability to a partially melted WC for some range of r s below r mit . We thus propose a “metallic electron crystal” phase of the two-dimensional electron gas at intermediate densities between a low density insulating WC and a high density Fermi fluid.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Melting anisotropy in crystalline solids

Despite the long and extensive history of melting studies, experiments to determine the dependence of melting on crystal orientation are lacking. Using longitudinal sound speed measurements in aluminum single crystals shock compressed along < 100 > and < 110 > to 168 GPa, we address this need and show that the melting transition (both onset and completion) is strongly anisotropic. Furthermore, the present findings demonstrate the need to consider the role of deformation induced lattice defects on the melting transition and provide a general approach to examine melting anisotropy in crystalline solids.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Melting curves of atomic hydrogen and deuterium calculated using path-integral Monte Carlo

We calculate the melting line of atomic hydrogen and deuterium up to 900 GPa with path-integral Monte Carlo using a machine-learned interatomic potential. We improve upon previous simulations of melting by treating the electrons with reptation quantum Monte Carlo, and by performing solid and liquid simulations using isothermal-isobaric path-integral Monte Carlo. Here, the resulting melting line for atomic hydrogen is higher than previous estimates. There is a small but resolvable decrease in the melting temperature as pressure is increased, which can be attributed to quantum effects.

08 HYDROGEN↗
DOE OSTI2023

Rocket drops: The self-propulsion of supercooled freezing drops

Isolated supercooled water drops have been observed to move spontaneously while freezing in vacuum. This motion is caused by an increase of the evaporation rate during the early stages of freezing, which transfers momentum to the drops, similar to rocket propulsion. As in other cases of self-propulsion, symmetry-breaking is necessary, and occurs when ice nucleation occurs away from the center of the drop. The self-propulsion velocity was modeled analytically for a simplified case, and numerically for two experiments. Here, the model predicts peak velocities on the order of 1 m/s in vacuum, and a drop kinematics similar to that observed experimentally. In air, the self-propulsion velocity is expected to be much smaller but may be detectable for micron-sized drops.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Anharmonicity and Phase Diagram of Magnesium Oxide in the Megabar Regime

With density functional molecular dynamics simulations, we computed the phase diagram of MgO from 50 to 2000 GPa up to 20000 K. Via thermodynamic integration (TDI), we derive the Gibbs free energies of the B1, B2, and liquid phases and determine their phase boundaries. With TDI and a pseudo-quasi-harmonic approach, we show that anharmonic effects are important and stabilize the B1 phase in particular. As a result, the B1-B2 transition boundary in the pressure-temperature plane exhibits a steep slope. We predict the B1-B2-liquid triple point to occur at approximately T= 10000 K and P = 370 GPa, which is higher in pressure than was inferred with quasi-harmonic methods alone. We predict the principal shock Hugoniot curve to enter the B2 phase stability domain but only over a very small range of parameters. Here, this may render it difficult to observe this phase with shock experiments because of kinetic effects.

3-dimensional systems↗
DOE OSTI2025

Microscopic changes governing melting anisotropy: Real-time nanosecond x-ray diffraction

To understand the microscopic origins governing melting anisotropy, in-situ x-ray diffraction (XRD) measurements were obtained in aluminum single crystals shocked along $\langle$100$\rangle$ and $\langle$110$\rangle$ to stresses below and above the melting threshold for each orientation. XRD results and analysis for the two orientations showed significant differences in the microstructure prior to the melting threshold stresses, demonstrating the key role of deformation induced microstructure – in addition to temperature and pressure – on the melting transition. As a result, our findings make a strong case for reconsidering theoretical approaches to melting, specifically for dislocation mediated melting, in shock compressed solids.

crystal melting↗
DOE OSTI2025

Dynamic compression effects of H 2 ⁡O in a dynamic diamond anvil cell: Origin of metastable ice VII and its crystal growth kinetics

We report on the structural verification of metastable ice VII solidifying in the phase space of ice VI at 1.80 GPa at room temperature. Using time-resolved (TR) x-ray diffraction and TR ruby luminescence paired with high-speed microphotography utilizing a dynamic diamond anvil cell, an initial compression rate range from 0.12 to 95.84 GPa/s was explored. The solidification pressure of metastable ice VII has a potential sigmoidal dependence upon compression rate with a turnover compression rate of ∼80 GPa/s. The preferred crystallization of ice VII in the stability field of ice VI is due to the increased nucleation rate of ice VII over ice VI at 1.77 GPa that is driven by the surface energy difference between the liquid and solid phases along with the change in Gibbs free energy of solidification. The dynamic pressure-volume–compression behaviors of ice phases (VI and VII) show a lattice stiffening in both phases, especially during the compression loading. It is also found that the compression rate greatly affects the solid-solid phase transition between ice VI and VII but does not affect the liquid-solid transition between water and ice VI as much. Lastly, a third phase transition was found to occur after metastable ice VII transforms into high-density amorphous (HDA) ice, which could be a disordered hydrogen-bonded network configuration of ice VII forming out of HDA ice facilitated by the decoupling of the oxygen movement and reorientation of the H 2⁡ O molecule. These results demonstrate the complexity of a seemingly simple molecule H 2⁡ O, how it can readily change its static properties with the modification of (de)compression rate, and highlight the need to use multiple TR structural and spectroscopic probes at higher time resolutions to realize the most comprehensive understanding.

Chemical bonding↗
DOE OSTI2023

Convection-Induced Compositional Patterning at Grain Boundaries in Irradiated Alloys

We consider the stability of precipitates formed at grain boundaries (GBs) by radiation-induced segregation in dilute alloys subjected to irradiation. Here, the effects of grain size and misorientation of symmetric-tilt GBs are quantified using phase field modeling. A novel regime is identified where, at long times, GBs are decorated by precipitate patterns that resist coarsening. Maps of the chemical Péclet number indicate that arrested coarsening takes place when solute advection dominates over thermal diffusion right up to the precipitate/matrix interface, preventing interfacial local equilibrium and overriding capillary effects. This contrasts with liquid-solid mixtures where convection always accelerates coarsening.

36 MATERIALS SCIENCE↗