Hydrogen liquid-liquid transition from first principles and machine learning

Abstract

The molecular-to-atomic liquid-liquid transition (LLT) in high-pressure hydrogen is a fundamental topic touching domains from planetary science to materials modeling. Yet, the nature of the LLT is still under debate. To resolve it, numerical simulations must cover length and time scales spanning several orders of magnitude. We overcome these size and time limitations by constructing a fast and accurate machine-learning interatomic potential (MLIP) built on the MACE neural network architecture. The MLIP is trained on Perdew-Burke-Ernzerhof (PBE) density functional calculations and uses a modified loss function correcting for an energy bias in the molecular phase. Classical and path-integral molecular dynamics driven by this MLIP show that the LLT is always supercritical above the melting temperature. The position of the corresponding Widom line agrees with previous ab initio PBE calculations, which in contrast predicted a first-order LLT. According to our calculations, the crossover line becomes a first-order transition only inside the molecular crystal region. These results call for a reconsideration of the LLT picture previously drawn.