Repository logo
 

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


Change log

Authors

Drummond, ND 
Monserrat, Bartomeu 
Lloyd-Williams, Jonathan H 
López Ríos, P 
Pickard, Chris J 

Abstract

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.

Description

Keywords

cond-mat.mtrl-sci, cond-mat.mtrl-sci

Journal Title

Nat Commun

Conference Name

Journal ISSN

2041-1723
2041-1723

Volume Title

6

Publisher

Springer Science and Business Media LLC
Sponsorship
Engineering and Physical Sciences Research Council (EP/J017639/1)
Engineering and Physical Sciences Research Council (EP/K037870/1)
Engineering and Physical Sciences Research Council (EP/F032773/1)
We thank Dominik Jochym for help with the implementation of the BLYP density functional. Financial support was provided by the Engineering and Physical Sciences Research Council (EPSRC), U.K. This research used resources of the Oak Ridge Leadership Computing Facility at the Oak Ridge National Laboratory, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725. Additional calculations were performed on the Cambridge High Performance Computing Service facility Darwin and the N8 high-performance computing facility provided and funded by the N8 consortium and EPSRC (Grant No. EP/K000225/1). We thank Dominik Jochym for help with the mplementation of the BLYP density functional. Financial support was provided by the Engineering and Physical Sciences Research Council (EPSRC), U.K. This research used resources of the Oak Ridge Leadership Computing Facility at the Oak Ridge National Laboratory, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725. Additional calculations were performed on the Cambridge High Performance Computing Service facility Darwin and the N8 high-performance computing facility provided and funded by the N8 consortium and EPSRC (Grant No. EP/K000225/1).