Melting and High P-T Transitions of Hydrogen up to 300 GPa

A platform for planar dynamic compression of crystalline hydrogen toward the terapascal regime

We describe a method for laser-driven planar compression of crystalline hydrogen that starts with a sample of solid para-hydrogen (even-valued rotational quantum number j) having an entropy of 0.06 kB/molecule at 10 K and 2 atm, with Boltzmann constant kB. Starting with this low-entropy state, the sample is compressed using a small initial shock (<0.2 GPa), followed by a pressure ramp that approaches isentropic loading as the sample is taken to hundreds of GPa. Planar loading allows for quantitative compression measurements; the objective of our low-entropy compression is to keep the sample cold enough to characterize crystalline hydrogen toward the terapascal range.

Thermodynamics of diamond formation from hydrocarbon mixtures in planets

Hydrocarbon mixtures are extremely abundant in the Universe, and diamond formation from them can play a crucial role in shaping the interior structure and evolution of planets. With first-principles accuracy, we first estimate the melting line of diamond, and then reveal the nature of chemical bonding in hydrocarbons at extreme conditions. We finally establish the pressure-temperature phase boundary where it is thermodynamically possible for diamond to form from hydrocarbon mixtures with different atomic fractions of carbon. Notably, here we show a depletion zone at pressures above 200 GPa and temperatures below 3000 K-3500 K where diamond formation is thermodynamically favorable regardless of the carbon atomic fraction, due to a phase separation mechanism. The cooler condition of the interior of Neptune compared to Uranus means that the former is much more likely to contain the depletion zone. Our findings can help explain the dichotomy of the two ice giants manifested by the low luminosity of Uranus, and lead to a better understanding of (exo-)planetary formation and evolution.

Stable Solid Molecular Hydrogen above 900 K from a Machine-Learned Potential Trained with Diffusion Quantum Monte Carlo

We survey the phase diagram of high-pressure molecular hydrogen with path integral molecular dynamics using a machine-learned interatomic potential trained with quantum Monte Carlo forces and energies. Besides the HCP and C2/c-24 phases, we find two new stable phases both with molecular centers in the Fmmm-4 structure, separated by a molecular orientation transition with temperature. The high temperature isotropic Fmmm-4 phase has a reentrant melting line with a maximum at higher temperature (1450 K at 150 GPa) than previously estimated and crosses the liquid-liquid transition line around 1200 K and 200 GPa.

Thermodynamic modeling of fluid polyamorphism in hydrogen at extreme conditions

Fluid polyamorphism, the existence of multiple amorphous fluid states in a single-component system, has been observed or predicted in a variety of substances. A remarkable example of this phenomenon is the fluid–fluid phase transition (FFPT) in high-pressure hydrogen between insulating and conducting high-density fluids. This transition is induced by the reversible dimerization/dissociation of the molecular and atomistic states of hydrogen. In this work, we present the first attempt to thermodynamically model the FFPT in hydrogen at extreme conditions. Our predictions for the phase coexistence and the reaction equilibrium of the two alternative forms of fluid hydrogen are based on experimental data and supported by the results of simulations. Remarkably, we find that the law of corresponding states can be utilized to construct a unified equation of state combining the available computational results for different models of hydrogen and the experimental data.

The 2021 room-temperature superconductivity roadmap

Designing materials with advanced functionalities is the main focus of contemporary solid-state physics and chemistry. Research efforts worldwide are funneled into a few high-end goals, one of the oldest, and most fascinating of which is the search for an ambient temperature superconductor (A-SC). The reason is clear: superconductivity at ambient conditions implies being able to handle, measure and access a single, coherent, macroscopic quantum mechanical state without the limitations associated with cryogenics and pressurization. This would not only open exciting avenues for fundamental research, but also pave the road for a wide range of technological applications, affecting strategic areas such as energy conservation and climate change. In this roadmap we have collected contributions from many of the main actors working on superconductivity, and asked them to share their personal viewpoint on the field. The hope is that this article will serve not only as an instantaneous picture of the status of research, but also as a true roadmap defining the main long-term theoretical and experimental challenges that lie ahead. Interestingly, although the current research in superconductor design is dominated by conventional (phonon-mediated) superconductors, there seems to be a widespread consensus that achieving A-SC may require different pairing mechanisms.In memoriam, to Neil Ashcroft, who inspired us all.

Fluid-like elastic response of superionic NH3 in Uranus and Neptune

Nondipolar magnetic fields exhibited at Uranus and Neptune may be derived from a unique geometry of their icy mantle with a thin convective layer on top of a stratified nonconvective layer. The presence of superionic H₂O and NH₃ has been thought as an explanation to stabilize such nonconvective regions. However, a lack of experimental data on the physical properties of those superionic phases has prevented the clarification of this matter. Here, our Brillouin measurements for NH₃ show a two-stage reduction in longitudinal wave velocity (V _p) by ∼9% and ∼20% relative to the molecular solid in the temperature range of 1,500 K and 2,000 K above 47 GPa. While the first V _p reduction observed at the boundary to the superionic α phase was most likely due to the onset of the hydrogen diffusion, the further one was likely attributed to the transition to another superionic phase, denoted γ phase, exhibiting the higher diffusivity. The reduction rate of V _p in the superionic γ phase, comparable to that of the liquid, implies that this phase elastically behaves almost like a liquid. Our measurements show that superionic NH₃ becomes convective and cannot contribute to the internal stratification.

Evidence for supercritical behaviour of high-pressure liquid hydrogen

Hydrogen, the simplest and most abundant element in the Universe, develops a remarkably complex behaviour upon compression¹. Since Wigner predicted the dissociation and metallization of solid hydrogen at megabar pressures almost a century ago², several efforts have been made to explain the many unusual properties of dense hydrogen, including a rich and poorly understood solid polymorphism^1,3-5, an anomalous melting line⁶ and the possible transition to a superconducting state⁷. Experiments at such extreme conditions are challenging and often lead to hard-to-interpret and controversial observations, whereas theoretical investigations are constrained by the huge computational cost of sufficiently accurate quantum mechanical calculations. Here we present a theoretical study of the phase diagram of dense hydrogen that uses machine learning to 'learn' potential-energy surfaces and interatomic forces from reference calculations and then predict them at low computational cost, overcoming length- and timescale limitations. We reproduce both the re-entrant melting behaviour and the polymorphism of the solid phase. Simulations using our machine-learning-based potentials provide evidence for a continuous molecular-to-atomic transition in the liquid, with no first-order transition observed above the melting line. This suggests a smooth transition between insulating and metallic layers in giant gas planets, and reconciles existing discrepancies between experiments as a manifestation of supercritical behaviour.

Distinguishing the Structure of High-Pressure Hydrogen with Dielectric Constants

Identifying the structures of high-pressure hydrogen has been one of the central goals in high-pressure physics; however, it still presents a fundamental challenge because of the lack of an effective measure for distinguishing the structures. Herein, we address this issue by focusing on the potential candidates of phases II and III of high-pressure hydrogen. We find that the anisotropic dielectric constants of the different hydrogen solids and their responses to pressure behave differently depending on the atomic structures, corresponding to the different polarization responses of the structures to the external electric field. These findings are robust regardless of the quantum and thermal motion of hydrogen solids. Therefore, the anisotropic dielectric property can serve as a potential measure for probing the structures of high-pressure hydrogen as well as other high-pressure materials.

A Spectroscopic Study of the Insulator-Metal Transition in Liquid Hydrogen and Deuterium

The insulator-to-metal transition in dense fluid hydrogen is an essential phenomenon in the study of gas giant planetary interiors and the physical and chemical behavior of highly compressed condensed matter. Using direct fast laser spectroscopy techniques to probe hydrogen and deuterium precompressed in a diamond anvil cell and laser heated on microsecond timescales, an onset of metal-like reflectance is observed in the visible spectral range at P >150 GPa and T ≥ 3000 K. The reflectance increases rapidly with decreasing photon energy indicating free-electron metallic behavior with a plasma edge in the visible spectral range at high temperatures. The reflectance spectra also suggest much longer electronic collision time (≥1 fs) than previously inferred, implying that metallic hydrogen at the conditions studied is not in the regime of saturated conductivity (Mott-Ioffe-Regel limit). The results confirm the existence of a semiconducting intermediate fluid hydrogen state en route to metallization.

Polarized angular-dependent reflectivity and density-dependent profiles of shock-compressed xenon plasmas

New data for the reflectivity of shock-compressed xenon plasmas at pressures of 10-12 GPa at large incident angles are presented. In addition, measurements have been performed at different densities. These data allow to analyze the free-electron density profile across the shock wave front. Assuming a Fermi-like density profile, the width of the front layer is inferred. The reflectivity coefficients for the s- and p-polarized waves are calculated. The influence of atoms, which was taken into account on the level of the collision frequency, proves to be essential for the understanding of the reflection process. Subsequently, a unique density profile is sufficient to obtain good agreement with the experimental data at different incident angles and at all investigated optical laser frequencies. Reflectivity measurements for different densities allow to determine the dependence of shock-front density profiles on the plasma parameters. As a result, it was found that the width of the front layer increases with decreasing density.

Optical properties of high-pressure fluid hydrogen across molecular dissociation

Optical properties of compressed fluid hydrogen in the region where dissociation and metallization is observed are computed by ab initio methods and compared with recent experimental results. We confirm that at T > 3,000 K, both processes are continuous, while at T < 1,500 K, the first-order phase transition is accompanied by a discontinuity of the dc conductivity and the thermal conductivity, while both the reflectivity and absorption coefficient vary rapidly but continuously. Our results support the recent analysis of National Ignition Facility (NIF) experiments [Celliers PM, et al. (2018) Science 361:677-682], which assigned the inception of metallization to pressures where the reflectivity is ∼0.3. Our results also support the conclusion that the temperature plateau seen in laser-heated diamond-anvil cell (DAC) experiments at temperatures higher than 1,500 K corresponds to the onset of optical absorption, not to the phase transition.

Anharmonic and Anomalous Trends in the High-Pressure Phase Diagram of Silicon

A multifaceted first-principles approach utilizing density functional theory, evolutionary algorithms, and lattice dynamics was used to construct the phase diagram of silicon up to 4 TPa and 26 000  K. These calculations predicted that (i) an anomalous sequence of face-centered cubic to body-centered cubic to simple cubic crystalline phase transitions occur at pressures of 2.87 and 3.89 TPa, respectively, along the cold curve, (ii) the orthorhombic phases of Imma and Cmce-16 appear on the phase diagram only when the anharmonic contribution to the Gibbs free energy is taken into account, and (iii) a substantial change in the slope of the principal Hugoniot is observed if the anharmonic free energy of the cubic diamond phase is considered.

Insulator-metal transition in dense fluid deuterium

Dense fluid metallic hydrogen occupies the interiors of Jupiter, Saturn, and many extrasolar planets, where pressures reach millions of atmospheres. Planetary structure models must describe accurately the transition from the outer molecular envelopes to the interior metallic regions. We report optical measurements of dynamically compressed fluid deuterium to 600 gigapascals (GPa) that reveal an increasing refractive index, the onset of absorption of visible light near 150 GPa, and a transition to metal-like reflectivity (exceeding 30%) near 200 GPa, all at temperatures below 2000 kelvin. Our measurements and analysis address existing discrepancies between static and dynamic experiments for the insulator-metal transition in dense fluid hydrogen isotopes. They also provide new benchmarks for the theoretical calculations used to construct planetary models.

Metallization and molecular dissociation of dense fluid nitrogen

Diatomic nitrogen is an archetypal molecular system known for its exceptional stability and complex behavior at high pressures and temperatures, including rich solid polymorphism, formation of energetic states, and an insulator-to-metal transformation coupled to a change in chemical bonding. However, the thermobaric conditions of the fluid molecular–polymer phase boundary and associated metallization have not been experimentally established. Here, by applying dynamic laser heating of compressed nitrogen and using fast optical spectroscopy to study electronic properties, we observe a transformation from insulating (molecular) to conducting dense fluid nitrogen at temperatures that decrease with pressure and establish that metallization, and presumably fluid polymerization, occurs above 125 GPa at 2500 K. Our observations create a better understanding of the interplay between molecular dissociation, melting, and metallization revealing features that are common in simple molecular systems.

Nitrogen is a model system still presenting unknown behaviors at the pressures and temperatures typical of deep planets’ interiors. Here the authors explore, by pulsed laser heating in a diamond anvil cell and optical measurements, the metallization and non-molecular states of nitrogen in a previously unexplored domain above 1 Mbar and at 2000-7000K.

Strong Electron-Phonon and Band Structure Effects in the Optical Properties of High Pressure Metallic Hydrogen

The recent claim of having produced metallic hydrogen in the laboratory relies on measurements of optical spectra. Here, we present first-principles calculations of the reflectivity of hydrogen between 400 and 600 GPa in the I4_{1}/amd crystal structure, the one predicted at these pressures, based on both time-dependent density functional and Eliashberg theories, thus, covering the optical properties from the infrared to the ultraviolet regimes. Our results show that atomic hydrogen displays an interband plasmon at around 6 eV that abruptly suppresses the reflectivity, while the large superconducting gap energy yields a sharp decrease of the reflectivity in the infrared region approximately at 120 meV. The experimentally estimated electronic scattering rates in the 0.7-3 eV range are in agreement with our theoretical estimations, which show that the huge electron-phonon interaction of the system dominates the electronic scattering in this energy range. The remarkable features in the optical spectra predicted here encourage extending the optical measurements to the infrared and ultraviolet regions as our results suggest optical measurements can potentially identify high-pressure phases of hydrogen.