Bonding changes in hot fluid hydrogen at megabar pressures

High temperature melting of dense molecular hydrogen from machine-learning interatomic potentials trained on quantum Monte Carlo

We present results and discuss methods for computing the melting temperature of dense molecular hydrogen using a machine learned model trained on quantum Monte Carlo data. In this newly trained model, we emphasize the importance of accurate total energies in the training. We integrate a two phase method for estimating the melting temperature with estimates from the Clausius-Clapeyron relation to provide a more accurate melting curve from the model. We make detailed predictions of the melting temperature, solid and liquid volumes, latent heat, and internal energy from 50 to 180 GPa for both classical hydrogen and quantum hydrogen. At pressures of roughly 173 GPa and 1635 K, we observe molecular dissociation in the liquid phase. We compare with previous simulations and experimental measurements.

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.

In-situ abiogenic methane synthesis from diamond and graphite under geologically relevant conditions

Diamond and graphite are fundamental sources of carbon in the upper mantle, and their reactivity with H₂-rich fluids present at these depths may represent the key to unravelling deep abiotic hydrocarbon formation. We demonstrate an unexpected high reactivity between carbons' most common allotropes, diamond and graphite, with hydrogen at conditions comparable with those in the Earth's upper mantle along subduction zone thermal gradients. Between 0.5-3 GPa and at temperatures as low as 300 °C, carbon reacts readily with H₂ yielding methane (CH₄), whilst at higher temperatures (500 °C and above), additional light hydrocarbons such as ethane (C₂H₆) emerge. These results suggest that the interaction between deep H₂-rich fluids and reduced carbon minerals may be an efficient mechanism for producing abiotic hydrocarbons at the upper mantle.

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.

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

High P-T Raman spectra of hydrogen in the vibron and lattice mode regions were measured up to 300 GPa and 900 K using externally heated diamond anvil cell techniques. A new melting line determined from the disappearance of lattice mode excitations was measured directly for the first time above 140 GPa. The results differ from theoretical predictions and extrapolations from lower pressure melting relations. In addition, discontinuities in Raman frequencies are observed as a function of pressure and temperature indicative of phase transition at these conditions. The appearance of a new Raman feature near 2700  cm^{-1} at ∼300  GPa and 370 K indicates the transformation to a new crystalline phase. Theoretical calculations of the spectrum suggest the new phase is the proposed Cmca-4 metallic phase. The transition pressure is close to that of a recently reported transition observed on dynamic compression.

Predicted reentrant melting of dense hydrogen at ultra-high pressures

The phase diagram of hydrogen is one of the most important challenges in high-pressure physics and astrophysics. Especially, the melting of dense hydrogen is complicated by dimer dissociation, metallization and nuclear quantum effect of protons, which together lead to a cold melting of dense hydrogen when above 500 GPa. Nonetheless, the variation of the melting curve at higher pressures is virtually uncharted. Here we report that using ab initio molecular dynamics and path integral simulations based on density functional theory, a new atomic phase is discovered, which gives an uplifting melting curve of dense hydrogen when beyond 2 TPa, and results in a reentrant solid-liquid transition before entering the Wigner crystalline phase of protons. The findings greatly extend the phase diagram of dense hydrogen, and put metallic hydrogen into the group of alkali metals, with its melting curve closely resembling those of lithium and sodium.

Structure and dynamics of liquid CS2: Going from ambient to elevated pressure conditions

Molecular dynamics simulation studies were performed to investigate the structural and dynamic properties of liquid carbon disulfide (CS₂) from ambient to elevated pressure conditions. The results obtained have revealed structural changes at high pressures, which are related to the more dense packing of the molecules inside the first solvation shell. The calculated neutron and X-ray structure factors have been compared with available experimental diffraction data, also revealing the pressure effects on the short-range structure of the liquid. The pressure effects on the translational, reorientational, and residence dynamics are very strong, revealing a significant slowing down when going from ambient pressure to 1.2 GPa. The translational dynamics of the linear CS₂ molecules have been found to be more anisotropic at elevated pressures, where cage effects and librational motions are reflected on the shape of the calculated time correlation functions and their corresponding spectral densities.

Optical Properties of Fluid Hydrogen at the Transition to a Conducting State

We use fast transient transmission and emission spectroscopies in the pulse laser heated diamond anvil cell to probe the energy-dependent optical properties of hydrogen at pressures of 10-150 GPa and temperatures up to 6000 K. Hydrogen is absorptive at visible to near-infrared wavelengths above a threshold temperature that decreases from 3000 K at 18 GPa to 1700 K at 110 GPa. Transmission spectra at 2400 K and 141 GPa indicate that the absorptive hydrogen is semiconducting or semimetallic in character, definitively ruling out a first-order insulator-metal transition in the studied pressure range.

Evidence for a new phase of dense hydrogen above 325 gigapascals

Almost 80 years ago it was predicted that, under sufficient compression, the H-H bond in molecular hydrogen (H2) would break, forming a new, atomic, metallic, solid state of hydrogen. Reaching this predicted state experimentally has been one of the principal goals in high-pressure research for the past 30 years. Here, using in situ high-pressure Raman spectroscopy, we present evidence that at pressures greater than 325 gigapascals at 300 kelvin, H2 and hydrogen deuteride (HD) transform to a new phase--phase V. This new phase of hydrogen is characterized by substantial weakening of the vibrational Raman activity, a change in pressure dependence of the fundamental vibrational frequency and partial loss of the low-frequency excitations. We map out the domain in pressure-temperature space of the suggested phase V in H2 and HD up to 388 gigapascals at 300 kelvin, and up to 465 kelvin at 350 gigapascals; we do not observe phase V in deuterium (D2). However, we show that the transformation to phase IV' in D2 occurs above 310 gigapascals and 300 kelvin. These values represent the largest known isotropic shift in pressure, and hence the largest possible pressure difference between the H2 and D2 phases, which implies that the appearance of phase V of D2 must occur at a pressure of above 380 gigapascals. These experimental data provide a glimpse of the physical properties of dense hydrogen above 325 gigapascals and constrain the pressure and temperature conditions at which the new phase exists. We speculate that phase V may be the precursor to the non-molecular (atomic and metallic) state of hydrogen that was predicted 80 years ago.

Phase boundary of hot dense fluid hydrogen

We investigated the phase transformation of hot dense fluid hydrogen using static high-pressure laser-heating experiments in a laser-heated diamond anvil cell. The results show anomalies in the heating efficiency that are likely to be attributed to the phase transition from a diatomic to monoatomic fluid hydrogen (plasma phase transition) in the pressure range between 82 and 106 GPa. This study imposes tighter constraints on the location of the hydrogen plasma phase transition boundary and suggests higher critical point than that predicted by the theoretical calculations.

Raman spectroscopy of hot hydrogen above 200 GPa

It has been theorized that at high pressure the increased energy of the zero-point oscillations in hydrogen would destabilize the lattice and form a ground fluid state at 0 K (ref. 1). Theory has also suggested that this fluid state, representing a new state of matter, might have unusual properties governed by quantum effects, such as superfluidity or superconductivity. Here, by combining Raman spectroscopy and in situ high-temperature, high-pressure techniques, we demonstrate that above 200 GPa a new phase transition occurs as temperature is increased, for example 480 K at 255 GPa. If the transformation is interpreted as melting, it would be the lowest melting temperature of any material at these high pressures. We also find a new triple point between phases I and IV and the new phase, and demonstrate that hydrogen retains its molecular character around this point. These data may require a significant revision of the phase diagram of hydrogen above 200 GPa.

Molecular to atomic phase transition in hydrogen under high pressure

The metallization of high-pressure hydrogen, together with the associated molecular to atomic transition, is one of the most important problems in the field of high-pressure physics. It is also currently a matter of intense debate due to the existence of conflicting experimental reports on the observation of metallic hydrogen on a diamond-anvil cell. Theoretical calculations have typically relied on a mean-field description of electronic correlation through density functional theory, a theory with well-known limitations in the description of metal-insulator transitions. In fact, the predictions of the pressure-driven dissociation of molecules in high-pressure hydrogen by density functional theory is strongly affected by the chosen exchange-correlation functional. In this Letter, we use quantum Monte Carlo calculations to study the molecular to atomic transition in hydrogen. We obtain a transition pressure of 447(3) GPa, in excellent agreement with the best experimental estimate of the transition 450 GPa based on an extrapolation to zero band gap from experimental measurements. Additionally, we find that C2/c is stable almost up to the molecular to atomic transition, in contrast to previous density functional theory (DFT) and DFT+quantum Monte Carlo studies which predict large stability regimes for intermediary molecular phases.

New iron hydrides under high pressure

The Fe-H system has been investigated by combined x-ray diffraction studies and total energy calculations at pressures up to 136 GPa. The experiments involve laser annealing of hydrogen-embedded iron in a diamond anvil cell. Two new FeHx compounds, with x∼2 and x=3, are discovered at 67 and 86 GPa, respectively. Their crystal structures are identified (unit cell and Fe positional parameters from x-ray diffraction, H positional parameters from ab initio calculations) as tetragonal with space group I4/mmm for FeH(∼2) and as simple cubic with space group Pm3m for FeH3. Large metastability regimes are observed that allowed to measure the P(V) equation of state at room temperature of FeH, FeH(∼2), and FeH3.

Dissociation of high-pressure solid molecular hydrogen: a quantum Monte Carlo and anharmonic vibrational study

A theoretical study is reported of the molecular-to-atomic transition in solid hydrogen at high pressure. We use the diffusion quantum Monte Carlo method to calculate the static lattice energies of the competing phases and a density-functional-theory-based vibrational self-consistent field method to calculate anharmonic vibrational properties. We find a small but significant contribution to the vibrational energy from anharmonicity. A transition from the molecular Cmca-12 direct to the atomic I41/amd phase is found at 374 GPa. The vibrational contribution lowers the transition pressure by 91 GPa. The dissociation pressure is not very sensitive to the isotopic composition. Our results suggest that quantum melting occurs at finite temperature.

High-pressure measurements of hydrogen phase IV using synchrotron infrared spectroscopy

Phase IV of dense solid hydrogen has been identified by its infrared spectrum using high-pressure synchrotron radiation techniques. The spectrum exhibits a sharp vibron band at higher frequency and lower intensity than that for phase III, indicating the stability of molecular H(2) with decreased intermolecular interactions and charge transfer between molecules. A low-frequency vibron having a strong negative pressure shift indicative of strongly interacting molecules is also observed. The character of the spectrum is consistent with an anisotropic, mixed layer structure related to those recently predicted theoretically. Phase IV was found to be stable from 220 GPa (300 K) to at least 340 GPa (near 200 K), with the I-III-IV triple point located. Infrared transmission observed to the lowest photon energies measured places constraints on the electronic properties of the phase.

Evidence of a liquid-liquid phase transition in hot dense hydrogen

We use pulsed-laser heating of hydrogen at static pressures in the megabar pressure region to search for the plasma phase transition to liquid atomic metallic hydrogen. We heat our samples substantially above the melting line and observe a plateau in a temperature vs. laser power curve that otherwise increases with power. This anomaly in the heating curve appears correlated with theoretical predictions for the plasma phase transition.

Nuclear quantum effects and nonlocal exchange-correlation functionals applied to liquid hydrogen at high pressure

Using first-principles molecular dynamics, we study the influence of nuclear quantum effects (NQEs) and nonlocal exchange-correlation density functionals (DFs) near molecular dissociation in liquid hydrogen. NQEs strongly influence intramolecular properties, such as bond stability, and are thus an essential part of the dissociation process. Moreover, by including DFs that account for either the self-interaction error or dispersion interactions, we find a much better description of molecular dissociation and metallization than previous studies based on classical protons and/or local or semilocal DFs. We obtain excellent agreement with experimentally measured optical properties along Hugoniot curves for precompressed states, and while we still find a first-order liquid-liquid transition at low temperatures, transition pressures are increased by more than 100 GPa.

A fresh look at dense hydrogen under pressure. I. an introduction to the problem, and an index probing equalization of H-H distances

In the first of a series of four papers on hydrogen under pressure, and its transitions from an initiating molecular state, we begin by defining carefully the problem, and setting the distance scale of interactions of protons and electrons in molecular aggregates of the first of the elements. Following a review of the experimental situation, in particular the phase diagram of hydrogen, in as much as it is known, and the behavior of its vibrons and rotons, we move onto the setting up of a numerical laboratory for probing the underlying physics and chemistry of interactions in hydrogen as the pressure increases. The laboratory consists of the preferred static structures emerging from calculations on the system in the range of 1 atm to 500 GPa, those of Pickard and Needs. The intermolecular (inter-pair) H···H separations naturally decrease with increasing pressure, first rapidly so, then more slowly. The intramolecular (intra-pair) H-H distances vary over a much smaller scale (0.05 Å) as the pressure increases, first decreasing, then increasing, and finally decreasing. We define an equalization function to gauge the approach to equality of the first neighbor and shortest next neighbor H (proton) separations in this numerical laboratory. And we find that metallization is likely to occur before bond equalization.

Synchrotron infrared measurements of dense hydrogen to 360 GPa

Diamond-anvil-cell techniques have been developed to confine and measure hydrogen samples under static conditions to pressures above 300 GPa from 12 to 300 K using synchrotron infrared and optical absorption techniques. A decreasing absorption threshold in the visible spectrum is observed, but the material remains transparent at photon energies down to 0.1 eV at pressures to 360 GPa over a broad temperature range. The persistence of the strong infrared absorption of the vibron characteristic of phase III indicates the stability of the paired state of hydrogen. There is no evidence for the predicted metallic state over these conditions, in contrast to recent reports, but electronic properties consistent with semimetallic behavior are observed.

Mixed molecular and atomic phase of dense hydrogen

We used Raman and visible transmission spectroscopy to investigate dense hydrogen (deuterium) up to 315 (275) GPa at 300 K. At around 200 GPa, we observe the phase transformation, which we attribute to phase III, previously observed only at low temperatures. This is succeeded at 220 GPa by a reversible transformation to a new phase, IV, characterized by the simultaneous appearance of the second vibrational fundamental and new low-frequency phonon excitations and a dramatic softening and broadening of the first vibrational fundamental mode. The optical transmission spectra of phase IV show an overall increase of absorption and a closing band gap which reaches 1.8 eV at 315 GPa. Analysis of the Raman spectra suggests that phase IV is a mixture of graphenelike layers, consisting of elongated H2 dimers experiencing large pairing fluctuations, and unbound H2 molecules.

Raman Spectroscopy at High Pressures

Raman spectroscopy is one of the most informative probes for studies of material properties under extreme conditions of high pressure. The Raman techniques have become more versatile over the last decades as a new generation of optical filters and multichannel detectors become available. Here, recent progress in the Raman techniques for high-pressure research and its applications in numerous scientific disciplines including physics and chemistry of materials under extremes, earth and planetary science, new materials synthesis, and high-pressure metrology will be discussed.