First-principles modeling of dense hydrogen is crucial in materials and planetary sciences. Despite its apparent simplicity, predicting the ionic and electronic structure of hydrogen is a formidable challenge, and it is connected with the insulator-to-metal transition, a century-old problem in condensed matter. Accurate simulations of liquid hydrogen are also essential for modeling gas giant planets. Here, we perform an exhaustive study of the equation of state of hydrogen using density functional theory (DFT) and quantum Monte Carlo simulations. We find that the pressure predicted by DFT may vary qualitatively when using different functionals. The predictive power of first-principles simulations is restored by validating each functional against higher-level wavefunction theories, represented by computationally intensive variational and diffusion Monte Carlo calculations. Our simulations provide evidence that hydrogen is denser at planetary conditions, compared to currently used equations of state. For Jupiter, this implies a lower bulk metallicity (i.e., a smaller mass of heavy elements). Our results further amplify the inconsistency between Jupiter’s atmospheric metallicity measured by the Galileo probe and the envelope metallicity inferred from interior models.
