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

Emerging superconductivity rules in rare-earth and alkaline-earth metal hydrides

Hydrides of alkaline-earth and rare-earth metals have garnered significant interest in high-temperature superconductor research due to their excellent electron-phonon coupling and high T c upon pressurization. This study explores the electronic structures and electron-phonon coupling of metal hydrides XHn (n = 4,6), where X includes Ca, Mg, Sc, and Y. The involvement of d-orbital electrons alters the Fermi surface, leading to saddle-point nesting and a charge density wave (CDW) phase transition, which opens the superconducting gap. For instance, in YH6, the exchange coupling between Y-4d and H-1s holes in the phonon softening region results in T c values up to 230 K. The study suggests that factors, such as the origin of the CDW order, hydrogen concentration, and d-orbital contributions are crucial to superconductivity. This work proposes a new rule for high T c superconductors, emphasizing the importance of double gaps and electron-phonon interactions at exchange coupling sites, and predicts potential high-quality superconductors among rare-earth hydrides.

High-pressure structures of solid hydrogen: Insights from ab initio molecular dynamics simulations

Understanding the structural behavior of solid hydrogen under high pressures is crucial for uncovering its unique properties and potential applications. In this study, starting from the phase I of solid hydrogen-free-rotator hcp structure, we conduct extensive ab initio molecular dynamics calculations to simulate the cooling, heating, and equilibrium processes within a pressure range of 80-260 GPa. Without relying on any structure previously predicted, we identify the high-pressure phase structures of solid hydrogen as P21/c for phase II, P6522 for phase III, and BG1BG2BG3 six-layer structure for phase IV, which are different from those proposed previously using the structure-search method. The reasonability of these structures are validated by Raman spectra and x-ray diffraction patterns by comparison with the experimental results. Our results actually show pronounced changes in the c/a ratio between phases I, III, and IV, which hold no brief for the experimental interpretation of an isostructural hcp transformations for phases I-III-IV.

On: X-ray diffraction from the electron gas in monatomic metallic hydrogen

Solid hydrogen is expected to become a monatomic metal under sufficiently high compression. With hydrogen having only a single valence electron and no ion core, the nature of x-ray diffraction patterns from the electron gas of monatomic metallic hydrogen is uncertain, and it is unclear whether they may yield enough information for a crystal structure determination. With emphasis on the Cs-IV-type (I41/amd) structure predicted for hydrogen at ∼500 GPa, the electron density distributions, zero-point and thermal atomic motion, and x-ray diffraction intensities are determined from first-principles calculations for several candidate phases of metallic hydrogen. It is shown that the electron distribution is much more structured than might be expected from the commonly employed free-electron-gas picture, and in fact more modulated than what is obtained from the superposition of free-atom charge densities. We demonstrate that an identification of the crystal structure of monatomic metallic hydrogen from x-ray diffraction is fundamentally possible and discuss the possibility of single-crystal diffraction from metallic hydrogen. An atomic scattering factor for the hydrogen atom in monatomic metallic hydrogen is constructed to aid the quantitative analysis of diffraction intensities from future x-ray diffraction experiments.

Stability and superconductivity of freestanding two-dimensional transition metal boridene: M4/3B2

The small atomic mass of boron indicates strong electron-phonon coupling (EPC), so it may have a brilliant performance in superconductivity. Recently, a new 2D boride sheet with ordered metal vacancies and surface terminals (Mo4/3B2-x) was realized in experiments (Zhouet al2021Science373801). Here, the 2D monolayer freestanding Mo4/3B2is evidenced to be thermodynamically stable. Through electronic structure, phonon spectrum and EPC, monolayer Mo4/3B2is found to be an intrinsic phonon-mediated superconductor. The superconducting transition temperature (Tc) is determined to be 4.06 K by the McMillian-Allen-Dynes formula. Remarkably, theTcof monolayer Mo4/3B2can be increased to 6.78 K with an appropriate biaxial tensile strain (+5%). Moreover, we predict that other transition metal replacing Mo atoms is also stable and retaining the superconductivity. Such as monolayer W4/3B2is also a superconductor with theTcof 2.37 K. Our research results enrich the database of 2D monolayer superconductors and boron-related formed materials science.

Prediction of a Supersolid Phase in High-Pressure Deuterium

Supersolid is a mysterious and puzzling state of matter whose possible existence has stirred a vigorous debate among physicists for over 60 years. Its elusive nature stems from the coexistence of two seemingly contradicting properties, long-range order and superfluidity. We report computational evidence of a supersolid phase of deuterium under high pressure (p>800  GPa) and low temperature (T<1.0  K). In our simulations, that are based on bosonic path integral molecular dynamics, we observe a highly concerted exchange of atoms while the system preserves its crystalline order. The exchange processes are favored by the soft core interactions between deuterium atoms that form a densely packed metallic solid. At the zero temperature limit, Bose-Einstein condensation is observed as the permutation probability of N deuterium atoms approaches 1/N with a finite superfluid fraction. Our study provides concrete evidence for the existence of a supersolid phase in high-pressure deuterium and could provide insights on the future investigation of supersolid phases in real materials.

High-Temperature Superconducting Phases in Cerium Superhydride with a T_{c} up to 115 K below a Pressure of 1 Megabar

The discoveries of high-temperature superconductivity in H_{3}S and LaH_{10} have excited the search for superconductivity in compressed hydrides, finally leading to the first discovery of a room-temperature superconductor in a carbonaceous sulfur hydride. In contrast to rapidly expanding theoretical studies, high-pressure experiments on hydride superconductors are expensive and technically challenging. Here, we experimentally discovered superconductivity in two new phases, Fm3[over ¯]m-CeH_{10} (SC-I phase) and P6_{3}/mmc-CeH_{9} (SC-II phase) at pressures that are much lower (<100  GPa) than those needed to stabilize other polyhydride superconductors. Superconductivity was evidenced by a sharp drop of the electrical resistance to zero and decreased critical temperature in deuterated samples and in external magnetic field. SC-I has T_{c}=115  K at 95 GPa, showing an expected decrease in further compression due to the decrease of the electron-phonon coupling (EPC) coefficient λ (from 2.0 at 100 GPa to 0.8 at 200 GPa). SC-II has T_{c}=57  K at 88 GPa, rapidly increasing to a maximum T_{c}∼100  K at 130 GPa, and then decreasing in further compression. According to the theoretical calculation, this is due to a maximum of λ at the phase transition from P6_{3}/mmc-CeH_{9} into a symmetry-broken modification C2/c-CeH_{9}. The pressure-temperature conditions of synthesis affect the actual hydrogen content and the actual value of T_{c}. Anomalously low pressures of stability of cerium superhydrides make them appealing for studies of superhydrides and for designing new superhydrides with stability at even lower pressures.

The stochastic self-consistent harmonic approximation: calculating vibrational properties of materials with full quantum and anharmonic effects

The efficient and accurate calculation of how ionic quantum and thermal fluctuations impact the free energy of a crystal, its atomic structure, and phonon spectrum is one of the main challenges of solid state physics, especially when strong anharmonicy invalidates any perturbative approach. To tackle this problem, we present the implementation on a modular Python code of the stochastic self-consistent harmonic approximation (SSCHA) method. This technique rigorously describes the full thermodynamics of crystals accounting for nuclear quantum and thermal anharmonic fluctuations. The approach requires the evaluation of the Born-Oppenheimer energy, as well as its derivatives with respect to ionic positions (forces) and cell parameters (stress tensor) in supercells, which can be provided, for instance, by first principles density-functional-theory codes. The method performs crystal geometry relaxation on the quantum free energy landscape, optimizing the free energy with respect to all degrees of freedom of the crystal structure. It can be used to determine the phase diagram of any crystal at finite temperature. It enables the calculation of phase boundaries for both first-order and second-order phase transitions from the Hessian of the free energy. Finally, the code can also compute the anharmonic phonon spectra, including the phonon linewidths, as well as phonon spectral functions. We review the theoretical framework of the SSCHA and its dynamical extension, making particular emphasis on the physical inter pretation of the variables present in the theory that can enlighten the comparison with any other anharmonic theory. A modular and flexible Python environment is used for the implementation, which allows for a clean interaction with other packages. We briefly present a toy-model calculation to illustrate the potential of the code. Several applications of the method in superconducting hydrides, charge-density-wave materials, and thermoelectric compounds are also reviewed.

Synthesis of molecular metallic barium superhydride: pseudocubic BaH12

Following the discovery of high-temperature superconductivity in the La-H system, we studied the formation of new chemical compounds in the barium-hydrogen system at pressures from 75 to 173 GPa. Using in situ generation of hydrogen from NH3BH3, we synthesized previously unknown superhydride BaH12 with a pseudocubic (fcc) Ba sublattice in four independent experiments. Density functional theory calculations indicate close agreement between the theoretical and experimental equations of state. In addition, we identified previously known P6/mmm-BaH2 and possibly BaH10 and BaH6 as impurities in the samples. Ab initio calculations show that newly discovered semimetallic BaH12 contains H2 and H3- molecular units and detached H12 chains which are formed as a result of a Peierls-type distortion of the cubic cage structure. Barium dodecahydride is a unique molecular hydride with metallic conductivity that demonstrates the superconducting transition around 20 K at 140 GPa.

Identification of synthesisable crystalline phases of water – a prototype for the challenges of computational materials design

We discuss the identification of experimentally realisable crystalline phases of water to outline and contextualise some of the diverse building blocks of a computational materials design process.

Speed of sound from fundamental physical constants

Two dimensionless fundamental physical constants, the fine structure constant α and the proton-to-electron mass ratio [Formula: see text], are attributed a particular importance from the point of view of nuclear synthesis, formation of heavy elements, planets, and life-supporting structures. Here, we show that a combination of these two constants results in a new dimensionless constant that provides the upper bound for the speed of sound in condensed phases, v_u We find that [Formula: see text], where c is the speed of light in vacuum. We support this result by a large set of experimental data and first-principles computations for atomic hydrogen. Our result expands the current understanding of how fundamental constants can impose new bounds on important physical properties.

Energy Gap Closure of Crystalline Molecular Hydrogen with Pressure

We study the gap closure with pressure of crystalline molecular hydrogen. The gaps are obtained from grand-canonical quantum Monte Carlo methods properly extended to quantum and thermal crystals, simulated by coupled electron ion Monte Carlo methods. Nuclear zero point effects cause a large reduction in the gap (∼2  eV). Depending on the structure, the fundamental indirect gap closes between 380 and 530 GPa for ideal crystals and 330-380 GPa for quantum crystals. Beyond this pressure the system enters into a bad metal phase where the density of states at the Fermi level increases with pressure up to ∼450-500  GPa when the direct gap closes. Our work partially supports the interpretation of recent experiments in high pressure hydrogen.

Nonlocal Electronic Correlations in the Cohesive Properties of High-Pressure Hydrogen Solids

High-pressure hydrogen exhibits remarkable phenomena including the insulator-to-metal (IM) transition; however, a complete resolution of its phase diagram is still an elusive goal despite many efforts and much controversy. Theoretical modeling is typically based on density functional theory (DFT) with a mean-field description of electronic correlations, which is known to be rather limited in describing IM transitions. Herein, we show that nonlocal electron correlations play a central role in the relative stability of solid hydrogen phases, and that DFT-correcting for these correlations by the many-body dispersion (MBD) model reaches the accuracy of quantum Monte Carlo (QMC) simulations and predicts the same C2/c-24 → Cmca-12 → Cs(IV) IM transition. In contrast with the conventional assumption that many-body electronic correlations become localized in metallic systems because of exponential screening with interelectronic distance, we find that the anisotropy of the electronic response of hydrogen solids under pressure leads to longer-ranged many-body effects in metallic phases relative to insulating ones. This refines our understanding of phase diagram of hydrogen solids as well as anisotropic many-body correlations.

Synthesis of clathrate cerium superhydride CeH9 at 80-100 GPa with atomic hydrogen sublattice

Hydrogen-rich superhydrides are believed to be very promising high-T_c superconductors. Recent experiments discovered superhydrides at very high pressures, e.g. FeH₅ at 130 GPa and LaH₁₀ at 170 GPa. With the motivation of discovering new hydrogen-rich high-T_c superconductors at lowest possible pressure, here we report the prediction and experimental synthesis of cerium superhydride CeH₉ at 80–100 GPa in the laser-heated diamond anvil cell coupled with synchrotron X-ray diffraction. Ab initio calculations were carried out to evaluate the detailed chemistry of the Ce-H system and to understand the structure, stability and superconductivity of CeH₉. CeH₉ crystallizes in a P6₃/mmc clathrate structure with a very dense 3-dimensional atomic hydrogen sublattice at 100 GPa. These findings shed a significant light on the search for superhydrides in close similarity with atomic hydrogen within a feasible pressure range. Discovery of superhydride CeH₉ provides a practical platform to further investigate and understand conventional superconductivity in hydrogen rich superhydrides.

Hydrogen-rich superhydrides are promising high-temperature superconductors which have been observed only at pressures above 170 GPa. Here the authors show that CeH₉ can be synthesized at 80-100 GPa with laser heating, and is characterized by a clathrate structure with a dense 3-dimensional atomic hydrogen sublattice.

Hydrogen-Induced High-Temperature Superconductivity in Two-Dimensional Materials: The Example of Hydrogenated Monolayer MgB_{2}

Hydrogen-based compounds under ultrahigh pressure, such as the polyhydrides H_{3}S and LaH_{10}, superconduct through the conventional electron-phonon coupling mechanism to attain the record critical temperatures known to date. Here we exploit the intrinsic advantages of hydrogen to strongly enhance phonon-mediated superconductivity in a completely different system, namely, a two-dimensional material with hydrogen adatoms. We find that van Hove singularities in the electronic structure, originating from atomiclike hydrogen states, lead to a strong increase of the electronic density of states at the Fermi level, and thus of the electron-phonon coupling. Additionally, the emergence of high-frequency hydrogen-related phonon modes in this system boosts the electron-phonon coupling further. As a concrete example, we demonstrate the effect of hydrogen adatoms on the superconducting properties of monolayer MgB_{2}, by solving the fully anisotropic Eliashberg equations, in conjunction with a first-principles description of the electronic and vibrational states, and their coupling. We show that hydrogenation leads to a high critical temperature of 67 K, which can be boosted to over 100 K by biaxial tensile strain.

Nuclear Magnetic Resonance Spectroscopy as a Dynamical Structural Probe of Hydrogen under High Pressure

An unambiguous crystallographic structure solution for the observed phases II-VI of high pressure hydrogen does not exist due to the failure of standard structural probes at extreme pressure. In this work we propose that nuclear magnetic resonance spectroscopy provides a complementary structural probe for high pressure hydrogen. We show that the best structural models available for phases II, III, and IV of high pressure hydrogen exhibit markedly distinct nuclear magnetic resonance spectra which could therefore be used to discriminate amongst them. As an example, we demonstrate how nuclear magnetic resonance spectroscopy could be used to establish whether phase III exhibits polymorphism. Our calculations also reveal a strong renormalization of the nuclear magnetic resonance response in hydrogen arising from quantum fluctuations, as well as a strong isotope effect. As the experimental techniques develop, nuclear magnetic resonance spectroscopy can be expected to become a useful complementary structural probe in high pressure experiments.

Metallic hydrogen

Hydrogen is the simplest and most abundant element in the Universe. There are two pathways for creating metallic hydrogen under high pressures. Over 80 years ago Wigner and Huntington predicted that if solid molecular hydrogen was sufficiently compressed in the T  =  0 K limit, molecules would dissociate to form atomic metallic hydrogen (MH). We have observed this transition at a pressure of 4.95 megabars. MH in this form has probably never existed on Earth or in the Universe; it may be a room temperature superconductor and is predicted to be metastable. If metastable it will have an important technological impact. Liquid metallic hydrogen can also be produced at intermediate pressures and high temperatures and is believed to make up ~90% of the planet Jupiter. We have observed this liquid-liquid transition, also known as the plasma phase transition, at pressures of ~1-2 megabar and temperatures ~1000-2000 K. However, in this paper we shall focus on the Wigner-Huntington transition. We shall discuss the methods used to observe metallic hydrogen at extreme conditions of static pressure in the laboratory, extending our understanding of the phase diagram of the simplest atom in the periodic table.

Structure and Metallicity of Phase V of Hydrogen

A new phase V of hydrogen was recently claimed in experiments above 325 GPa and 300 K. Because of the extremely small sample size at such record pressures the measurements were limited to Raman spectroscopy. The experimental data on increase of pressure show decreasing Raman activity and darkening of the sample, which suggests band gap closure and impending molecular dissociation, but no definite conclusions could be reached. Furthermore, the available data are insufficient to determine the structure of phase V, which remains unknown. Introducing saddle-point ab initio random structure searching, we find several new structural candidates of hydrogen which could describe the observed properties of phase V. We investigate hydrogen metallization in the proposed candidate structures, and demonstrate that smaller band gaps are associated with longer bond lengths. We conclude that phase V is a stepping stone towards metallization.

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.

The role of van der Waals and exchange interactions in high-pressure solid hydrogen

Our study of the van der Waals interactions in solid molecular hydrogen structures indicates two candidates for phase III.

Anharmonic enhancement of superconductivity in metallic molecular Cmca - 4 hydrogen at high pressure: a first-principles study

First-principles calculations based on density-functional theory including anharmonicity within the variational stochastic self-consistent harmonic approximation are applied to understand how the quantum character of the proton affects the candidate metallic molecular Cmca  -  4 structure of hydrogen in the 400-450 GPa pressure range, where metallization of hydrogen is expected to occur. Anharmonic effects, which become crucial due to the zero-point motion, have a large impact on the hydrogen molecules by increasing the intramolecular distance by approximately a 6%. This induces two new electron pockets at the Fermi surface opening new scattering channels for the electron-phonon interaction. Consequently, the electron-phonon coupling constant and the superconducting critical temperature are approximately doubled by anharmonicity and Cmca  -  4 hydrogen becomes a superconductor above 200 K in all the studied pressure range. Contrary to many superconducting hydrides, where anharmoncity tends to lower the superconducting critical temperature, our results show that it can enhance superconductivity in molecular hydrogen.

Topological Surface States in Dense Solid Hydrogen

Metallization of dense hydrogen and associated possible high-temperature superconductivity represents one of the key problems of physics. Recent theoretical studies indicate that before becoming a good metal, compressed solid hydrogen passes through a semimetallic stage. We show that such semimetallic phases predicted to be the most stable at multimegabar (∼300  GPa) pressures are not conventional semimetals: they exhibit topological metallic surface states inside the bulk "direct" gap in the two-dimensional surface Brillouin zone; that is, metallic surfaces may appear even when the bulk of the material remains insulating. Examples include hydrogen in the Cmca-12 and Cmca-4 structures; Pbcn hydrogen also has metallic surface states but they are of a nontopological nature. The results provide predictions for future measurements, including probes of possible surface superconductivity in dense hydrogen.

Chemical accuracy from quantum Monte Carlo for the benzene dimer

We report an accurate study of interactions between benzene molecules using variational quantum Monte Carlo (VMC) and diffusion quantum Monte Carlo (DMC) methods. We compare these results with density functional theory using different van der Waals functionals. In our quantum Monte Carlo (QMC) calculations, we use accurate correlated trial wave functions including three-body Jastrow factors and backflow transformations. We consider two benzene molecules in the parallel displaced geometry, and find that by highly optimizing the wave function and introducing more dynamical correlation into the wave function, we compute the weak chemical binding energy between aromatic rings accurately. We find optimal VMC and DMC binding energies of -2.3(4) and -2.7(3) kcal/mol, respectively. The best estimate of the coupled-cluster theory through perturbative triplets/complete basis set limit is -2.65(2) kcal/mol [Miliordos et al., J. Phys. Chem. A 118, 7568 (2014)]. Our results indicate that QMC methods give chemical accuracy for weakly bound van der Waals molecular interactions, comparable to results from the best quantum chemistry methods.

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.

Quantum Monte Carlo study of the phase diagram of solid molecular hydrogen at extreme pressures

Establishing the phase diagram of hydrogen is a major challenge for experimental and theoretical physics. Experiment alone cannot establish the atomic structure of solid hydrogen at high pressure, because hydrogen scatters X-rays only weakly. Instead, our understanding of the atomic structure is largely based on density functional theory (DFT). By comparing Raman spectra for low-energy structures found in DFT searches with experimental spectra, candidate atomic structures have been identified for each experimentally observed phase. Unfortunately, DFT predicts a metallic structure to be energetically favoured at a broad range of pressures up to 400 GPa, where it is known experimentally that hydrogen is non-metallic. Here we show that more advanced theoretical methods (diffusion quantum Monte Carlo calculations) find the metallic structure to be uncompetitive, and predict a phase diagram in reasonable agreement with experiment. This greatly strengthens the claim that the candidate atomic structures accurately model the experimentally observed phases.

Molecular-Atomic Transition along the Deuterium Hugoniot Curve with Coupled Electron-Ion Monte Carlo Simulations

We have performed simulations of the principal deuterium Hugoniot curve using coupled electron-ion Monte Carlo calculations. Using highly accurate quantum Monte Carlo methods for the electrons, we study the region of maximum compression along the Hugoniot, where the system undergoes a continuous transition from a molecular fluid to a monatomic fluid. We include all relevant physical corrections so that a direct comparison to experiment can be made. Around 50 GPa we find a maximum compression of 4.85. This compression is approximately 5.5% higher than previous theoretical predictions and 15% higher than the most accurate experimental data. Thus first-principles simulations encompassing the most advanced techniques are in disagreement with the results of the best experiments.

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.