A Denser Hydrogen Inferred from First-Principles Simulations Challenges Jupiter's Interior Models
Pith reviewed 2026-05-23 05:33 UTC · model grok-4.3
The pith
Simulations show hydrogen is denser at Jupiter conditions than standard models assume, implying lower bulk metallicity.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
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.
What carries the argument
Density functional theory functionals for the hydrogen equation of state validated against variational and diffusion Monte Carlo calculations.
Load-bearing premise
The density functionals selected after quantum Monte Carlo validation accurately capture the true equation of state of hydrogen across the full range of planetary pressures and temperatures without residual systematic bias.
What would settle it
An independent high-accuracy calculation or experiment measuring the density of liquid hydrogen at pressures of hundreds of GPa and temperatures of several thousand K that matches the lower densities from existing equations of state.
Figures
read the original abstract
First-principle 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 and quantum Monte Carlo simulations. We find that the pressure predicted by Density Functional Theory may vary qualitatively when using different functionals. The predictive power of first-principle 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents an exhaustive first-principles study of the equation of state of dense hydrogen using Density Functional Theory (DFT) and quantum Monte Carlo (QMC) methods. By validating various DFT functionals against variational and diffusion Monte Carlo calculations, the authors find that hydrogen is denser at the pressure-temperature conditions relevant to gas giant planets than predicted by currently used equations of state. This leads to the conclusion that Jupiter's bulk metallicity must be lower, exacerbating the discrepancy with its atmospheric metallicity measured by the Galileo probe.
Significance. If the central result holds, the work would require revisions to standard interior models of Jupiter and other gas giants by implying a smaller total mass of heavy elements. The explicit anchoring in QMC validation (rather than direct fitting to planetary data) and the parameter-free character of the density comparison are strengths that would elevate the credibility of first-principles constraints in planetary science.
major comments (1)
- [Abstract and functional validation paragraph] Abstract, paragraph on functional validation: the manuscript states that predictive power is restored by QMC validation at selected points, yet provides no quantitative error bars on the DFT-QMC density differences, no convergence tests with respect to system size or k-point sampling, and no explicit assessment of residual bias when extrapolating across the full planetary P-T range. This directly affects the load-bearing claim that hydrogen is systematically denser than current EOS.
minor comments (2)
- Add explicit statements of the finite-size corrections applied in the QMC calculations and how they propagate into the reported densities.
- Clarify in the methods section which specific DFT functionals were retained after the QMC comparison and the quantitative selection criterion used.
Simulated Author's Rebuttal
We thank the referee for their careful and constructive review. The single major comment raises valid points about the presentation of our validation procedures. We address it below and will revise the manuscript to incorporate the requested quantitative details.
read point-by-point responses
-
Referee: [Abstract and functional validation paragraph] Abstract, paragraph on functional validation: the manuscript states that predictive power is restored by QMC validation at selected points, yet provides no quantitative error bars on the DFT-QMC density differences, no convergence tests with respect to system size or k-point sampling, and no explicit assessment of residual bias when extrapolating across the full planetary P-T range. This directly affects the load-bearing claim that hydrogen is systematically denser than current EOS.
Authors: We agree that the current manuscript would benefit from more explicit quantification. In the revised version we will add (i) quantitative error bars on all reported DFT-QMC density differences, obtained directly from the statistical uncertainties of the QMC runs; (ii) a new appendix or supplementary section documenting convergence tests with respect to system size (including results up to 128–256 atoms) and k-point sampling for the DFT calculations; and (iii) an explicit discussion of possible residual bias when extrapolating across the planetary P-T range, including a sensitivity analysis of how plausible bias levels would affect the inferred Jupiter metallicity. These additions will be placed in the main text or supplementary material as appropriate to support the central claim. revision: yes
Circularity Check
No significant circularity
full rationale
The paper's central result follows from direct first-principles DFT simulations whose functionals are validated at selected points by independent variational and diffusion Monte Carlo calculations. No parameter is fitted to Jupiter interior data, no prediction is defined in terms of the target density or metallicity, and no load-bearing step reduces to a self-citation or ansatz imported from the same authors. The derivation chain therefore remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Density-functional approximations, once benchmarked against variational and diffusion Monte Carlo, yield reliable pressures and densities for dense liquid hydrogen.
- domain assumption Quantum Monte Carlo calculations provide a higher-accuracy reference for the hydrogen equation of state.
Forward citations
Cited by 1 Pith paper
-
Neural Wave Functions for High-Pressure Atomic Hydrogen
Neural quantum states yield Born-Oppenheimer and non-Born-Oppenheimer energies for high-pressure atomic hydrogen that match or beat prior projector Monte Carlo results up to 128 atoms while avoiding symmetry assumptio...
Reference graph
Works this paper leans on
-
[1]
for further details). Note that the heavy-element mass fraction (i.e., the envelope metallicity) is assumed to be the same in the inner and outer envelopes. The left panel of Figure 8 shows the inferred equatorial radius and J2 for the Jupiter model with the MH13 EoS (blue star). Using CEPAM [137], we calculated a similar Jupiter model with SCAN+vv10 (ora...
work page 2000
-
[2]
J. M. McMahon, M. A. Morales, C. Pierleoni, and D. M. Ceperley, The properties of hydrogen and helium under extreme conditions, Reviews of Modern Physics84, 1607 (2012)
work page 2012
- [3]
-
[4]
E. Wigner and H. B. Huntington, On the possibility of a metallic modification of hydrogen, The Journal of Chemical Physics3, 764 (1935)
work page 1935
-
[5]
V. L. Ginzburg, Nobel lecture: On superconductivity and superfluidity (what i have and have not managed to do) as well as on the “physical minimum” at the beginning of the xxi century, Rev. Mod. Phys.76, 981 (2004)
work page 2004
-
[6]
P. Dalladay-Simpson, R. T. Howie, and E. Gregoryanz, Evidence for a new phase of dense hydrogen above 325 gigapascals, Nature529, 63 (2016)
work page 2016
-
[7]
R. P. Dias and I. F. Silvera, Observation of the wigner- huntington transition to metallic hydrogen, Science 10.1126/science.aal1579 (2017)
-
[8]
P. Loubeyre, F. Occelli, and P. Dumas, Synchrotron infrared spectroscopic evidence of the probable transition to metal hydrogen, Nature577, 631 (2020)
work page 2020
-
[9]
S. Weir, A. Mitchell, and W. Nellis, Metallization of fluid molecular hydrogen at 140 gpa (1.4 mbar), Physical Review Letters76, 1860 (1996)
work page 1996
-
[10]
P. M. Celliers, M. Millot, S. Brygoo, R. S. McWilliams, D. E. Fratanduono, J. R. Rygg, A. F. Goncharov, P. Loubeyre, J. H. Eggert, J. L. Peterson, N. B. Meezan, S. Le Pape, G. W. Collins, R. Jeanloz, and R. J. Hemley, Insulator-metal transition in dense fluid deuterium, Science361, 677 (2018)
work page 2018
-
[11]
M. Knudson, M. Desjarlais, A. Becker, R. Lemke, K. Cochrane, M. Savage, D. Bliss, T. Mattsson, and R. Redmer, Direct observation of an abrupt insulator- to-metal transition in dense liquid deuterium, Science 348, 1455 (2015)
work page 2015
- [12]
-
[13]
R. Helled and J. J. Fortney, The interiors of Uranus and Neptune: current understanding and open questions, Philosophical Transactions of the Royal Society of London Series A378, 20190474 (2020), arXiv:2007.10783 [astro-ph.EP]
- [14]
-
[15]
D. Durante, M. Parisi, D. Serra, M. Zannoni, V. Notaro, P. Racioppa, D. R. Buccino, G. Lari, L. Gomez Casajus, L. Iess, W. M. Folkner, G. Tommei, P. Tortora, and S. J. Bolton, Jupiter’s Gravity Field Halfway Through the Juno Mission, Geophysical Research Letter47, e86572 (2020)
work page 2020
-
[16]
L. Iess, B. Militzer, Y. Kaspi, P. Nicholson, D. Durante, P. Racioppa, A. Anabtawi, E. Galanti, W. Hubbard, M. J. Mariani, P. Tortora, S. Wahl, and M. Zannoni, Measurement and implications of Saturn’s gravity field and ring mass, Science364, aat2965 (2019)
work page 2019
-
[17]
R. Helled and S. Howard, Giant planet interiors and atmospheres, arXiv e-prints , arXiv:2407.05853 (2024), arXiv:2407.05853 [astro-ph.EP]
- [18]
-
[19]
S. Howard, T. Guillot, M. Bazot, Y. Miguel, D. J. Stevenson, E. Galanti, Y. Kaspi, W. B. Hubbard, B. Militzer, R. Helled, N. Nettelmann, B. Idini, and S. Bolton, Jupiter’s interior from Juno: Equation-of- state uncertainties and dilute core extent, Astronomy and Astrophysics672, A33 (2023), arXiv:2302.09082 [astro-ph.EP]
-
[20]
Y. Miguel, M. Bazot, T. Guillot, S. Howard, E. Galanti, Y. Kaspi, W. B. Hubbard, B. Militzer, R. Helled, S. K. Atreya, J. E. P. Connerney, D. Durante, L. Kulowski, J. I. Lunine, D. Stevenson, and S. Bolton, Jupiter’s inhomogeneous envelope, Astronomy & Astrophysics 662, A18 (2022), arXiv:2203.01866 [astro-ph.EP]
-
[21]
S. e. a. Wahl, Comparing Jupiter interior structure models to Juno gravity measurements and the role of a dilute core, Geophys. Res. Lett.44, 80 (2017)
work page 2017
-
[22]
N. Nettelmann, N. Movshovitz, D. Ni, J. J. Fortney, E. Galanti, Y. Kaspi, R. Helled, C. R. Mankovich, and S. Bolton, Theory of Figures to the Seventh Order and the Interiors of Jupiter and Saturn, The Planetary Science Journal2, 241 (2021), arXiv:2110.15452 [astro- ph.EP]
-
[23]
B. Militzer, W. B. Hubbard, S. Wahl, J. I. Lunine, E. Galanti, Y. Kaspi, Y. Miguel, T. Guillot, K. M. Moore, M. Parisi, J. E. P. Connerney, R. Helled, H. Cao, C. Mankovich, D. J. Stevenson, R. S. Park, M. Wong, S. K. Atreya, J. Anderson, and S. J. Bolton, Juno Spacecraft Measurements of Jupiter’s Gravity Imply a Dilute Core, The Planetary Science Journal3...
work page 2022
-
[24]
D. J. Stevenson, Jupiter’s Interior as Revealed by Juno, Annual Review of Earth and Planetary Sciences48, 465 (2020)
work page 2020
-
[25]
M. D. Knudson, M. P. Desjarlais, R. W. Lemke, T. R. Mattsson, M. French, N. Nettelmann, and R. Redmer, Probing the interiors of the ice giants: Shock compression of water to 700 gpa and 3.8g/cm 3, Phys. Rev. Lett.108, 091102 (2012)
work page 2012
-
[26]
M. D. Knudson, D. L. Hanson, J. E. Bailey, C. A. Hall, J. R. Asay, and C. Deeney, Principal Hugoniot, reverberating wave, and mechanical reshock measurements of liquid deuterium to 400 GPa using plate impact techniques, Phys. Rev. B69, 144209 (2004)
work page 2004
-
[27]
D. G. Hicks, T. R. Boehly, P. M. Celliers, J. H. Eggert, S. J. Moon, D. D. Meyerhofer, and G. W. Collins, Laser- driven single shock compression of fluid deuterium from 45 to 220 GPa, Phys. Rev. B79, 014112 (2009)
work page 2009
- [28]
-
[29]
N. Nettelmann, A. Becker, B. Holst, and R. Redmer, Jupiter models with improved ab initio hydrogen 28 equation of state (h-reos.2), The Astrophysical Journal 750, 52 (2012)
work page 2012
-
[30]
B. Militzer and W. B. Hubbard, Ab initio equation of state for hydrogen-helium mixtures with recalibration of the giant-planet mass-radius relation, The Astrophysical Journal774, 148 (2013)
work page 2013
-
[31]
Miguel, Y., Guillot, T., and Fayon, L., Jupiter internal structure: the effect of different equations of state (corrigendum), A&A618, C2 (2018)
work page 2018
- [32]
-
[33]
G. Chabrier, S. Mazevet, and F. Soubiran, A new equation of state for dense hydrogen–helium mixtures, The Astrophysical Journal872, 51 (2019)
work page 2019
-
[34]
A. Becker, W. Lorenzen, J. J. Fortney, N. Nettelmann, M. Sch¨ ottler, and R. Redmer, Ab initio equations of state for hydrogen (h-reos. 3) and helium (he-reos. 3) and their implications for the interior of brown dwarfs, The Astrophysical Journal Supplement Series215, 21 (2014)
work page 2014
-
[35]
Burke, Perspective on density functional theory, The Journal of Chemical Physics136, 150901 (2012)
K. Burke, Perspective on density functional theory, The Journal of Chemical Physics136, 150901 (2012)
work page 2012
-
[36]
S. Azadi and W. M. C. Foulkes, Fate of density functional theory in the study of high-pressure solid hydrogen, Physical Review B88, 014115 (2013)
work page 2013
-
[37]
G. Mazzola, R. Helled, and S. Sorella, Phase diagram of hydrogen and a hydrogen-helium mixture at planetary conditions by quantum monte carlo simulations, Phys. Rev. Lett.120, 025701 (2018)
work page 2018
-
[38]
C. Pierleoni, M. A. Morales, G. Rillo, M. Holzmann, and D. M. Ceperley, Liquid–liquid phase transition in hydrogen by coupled electron–ion monte carlo simulations, Proceedings of the National Academy of Sciences113, 4953 (2016)
work page 2016
-
[39]
R. C. Clay, J. Mcminis, J. M. McMahon, C. Pierleoni, D. M. Ceperley, and M. A. Morales, Benchmarking exchange-correlation functionals for hydrogen at high pressures using quantum monte carlo, Phys. Rev. B89, 184106 (2014)
work page 2014
-
[40]
M. D. Knudson, M. P. Desjarlais, M. Preising, and R. Redmer, Evaluation of exchange-correlation functionals with multiple-shock conductivity measurements in hydrogen and deuterium at the molecular-to-atomic transition, Phys. Rev. B98, 174110 (2018)
work page 2018
-
[41]
M. Bonitz, J. Vorberger, M. Bethkenhagen, M. P. B¨ ohme, D. M. Ceperley, A. Filinov, T. Gawne, F. Graziani, G. Gregori, P. Hamann, S. B. Hansen, M. Holzmann, S. X. Hu, H. K¨ ahlert, V. V. Karasiev, U. Kleinschmidt, L. Kordts, C. Makait, B. Militzer, Z. A. Moldabekov, C. Pierleoni, M. Preising, K. Ramakrishna, R. Redmer, S. Schwalbe, P. Svensson, and T. Do...
-
[42]
S. Van De Bund, H. Wiebe, and G. J. Ackland, Isotope quantum effects in the metallization transition in liquid hydrogen, Physical Review Letters126, 225701 (2021)
work page 2021
- [43]
-
[44]
N. Schuch and F. Verstraete, Computational complexity of interacting electrons and fundamental limitations of density functional theory, Nature physics5, 732 (2009)
work page 2009
-
[45]
R. O. Jones, Density functional theory: Its origins, rise to prominence, and future, Rev. Mod. Phys.87, 897 (2015)
work page 2015
-
[46]
K. Lejaeghere, G. Bihlmayer, T. Bj¨ orkman, P. Blaha, S. Bl¨ ugel, V. Blum, D. Caliste, I. E. Castelli, S. J. Clark, A. Dal Corso,et al., Reproducibility in density functional theory calculations of solids, Science351, aad3000 (2016)
work page 2016
-
[47]
E. Bosoni, L. Beal, M. Bercx, P. Blaha, S. Bl¨ ugel, J. Br¨ oder, M. Callsen, S. Cottenier, A. Degomme, V. Dikan, K. Eimre, E. Flage-Larsen, M. Fornari, A. Garcia, L. Genovese, M. Giantomassi, S. P. Huber, H. Janssen, G. Kastlunger, M. Krack, G. Kresse, T. D. K¨ uhne, K. Lejaeghere, G. K. H. Madsen, M. Marsman, N. Marzari, G. Michalicek, H. Mirhosseini, T...
work page 2023
-
[48]
P. Hohenberg and W. Kohn, Inhomogeneous electron gas, Phys. Rev.136, B864 (1964)
work page 1964
-
[49]
W. Kohn and L. J. Sham, Self-consistent equations including exchange and correlation effects, Phys. Rev. 140, A1133 (1965)
work page 1965
-
[50]
R. M. Martin,Electronic Structure: Basic Theory and Practical Methods(Cambridge University Press, 2004)
work page 2004
-
[51]
J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett. 77, 3865 (1996)
work page 1996
-
[52]
G. Chabrier and F. Debras, A new equation of state for dense hydrogen–helium mixtures. ii. taking into account hydrogen–helium interactions, The Astrophysical Journal917, 4 (2021)
work page 2021
-
[53]
A. D. Becke, Density-functional exchange-energy approximation with correct asymptotic behavior, Physical review A38, 3098 (1988)
work page 1988
-
[54]
C. Lee, W. Yang, and R. G. Parr, Development of the colle-salvetti correlation-energy formula into a functional of the electron density, Physical review B37, 785 (1988)
work page 1988
-
[55]
M. Dion, H. Rydberg, E. Schr¨ oder, D. C. Langreth, and B. I. Lundqvist, Van der waals density functional for general geometries, Physical review letters92, 246401 (2004)
work page 2004
-
[56]
T. Thonhauser, S. Zuluaga, C. Arter, K. Berland, E. Schr¨ oder, and P. Hyldgaard, Spin signature of nonlocal correlation binding in metal-organic frameworks, Physical review letters115, 136402 (2015)
work page 2015
-
[57]
K. Lee, ´E. D. Murray, L. Kong, B. I. Lundqvist, and D. C. Langreth, Higher-accuracy van der waals density functional, Physical Review B—Condensed Matter and Materials Physics82, 081101 (2010)
work page 2010
-
[58]
J. Sun, A. Ruzsinszky, and J. P. Perdew, Strongly constrained and appropriately normed semilocal density functional, Physical review letters115, 036402 (2015)
work page 2015
-
[59]
H.-C. Yang, K. Liu, Z.-Y. Lu, and H.-Q. Lin, First- principles study of solid hydrogen: Comparison among 29 four exchange-correlation functionals, Phys. Rev. B 102, 174109 (2020)
work page 2020
-
[60]
M. Sch¨ ottler and R. Redmer, Ab initio calculation of the miscibility diagram for hydrogen-helium mixtures, Phys. Rev. Lett.120, 115703 (2018)
work page 2018
-
[61]
S. Azadi and G. J. Ackland, The role of van der waals and exchange interactions in high-pressure solid hydrogen, Physical Chemistry Chemical Physics19, 21829 (2017)
work page 2017
-
[62]
B. Lu, D. Kang, D. Wang, T. Gao, and J. Dai, Towards the same line of liquid–liquid phase transition of dense hydrogen from various theoretical predictions, Chinese Physics Letters36, 103102 (2019)
work page 2019
-
[63]
J. Hinz, V. V. Karasiev, S. Hu, M. Zaghoo, D. Mej´ ıa- Rodr´ ıguez, S. Trickey, and L. Calder´ ın, Fully consistent density functional theory determination of the insulator- metal transition boundary in warm dense hydrogen, Physical Review Research2, 032065 (2020)
work page 2020
- [64]
-
[65]
M. Zaghoo and I. F. Silvera, Conductivity and dissociation in liquid metallic hydrogen and implications for planetary interiors, Proceedings of the National Academy of Sciences114, 11873 (2017)
work page 2017
-
[66]
C. Ji, B. Li, W. Liu, J. S. Smith, A. Majumdar, W. Luo, R. Ahuja, J. Shu, J. Wang, S. Sinogeikin,et al., Ultrahigh-pressure isostructural electronic transitions in hydrogen, Nature573, 558 (2019)
work page 2019
-
[67]
M. A. Morales, J. M. McMahon, C. Pierleoni, and D. M. Ceperley, Nuclear quantum effects and nonlocal exchange-correlation functionals applied to liquid hydrogen at high pressure, Physical Review Letters110, 065702 (2013)
work page 2013
-
[68]
H. Y. Geng, Q. Wu, M. Marqu´ es, and G. J. Ackland, Thermodynamic anomalies and three distinct liquid- liquid transitions in warm dense liquid hydrogen, Phys. Rev. B100, 134109 (2019)
work page 2019
-
[69]
R. C. Clay III, M. Holzmann, D. M. Ceperley, and M. A. Morales, Benchmarking density functionals for hydrogen-helium mixtures with quantum monte carlo: Energetics, pressures, and forces, Physical Review B93, 035121 (2016)
work page 2016
-
[70]
A. J. Cohen, P. Mori-S´ anchez, and W. Yang, Challenges for density functional theory, Chemical Reviews112, 289 (2011)
work page 2011
-
[71]
L. O. Wagner, E. Stoudenmire, K. Burke, and S. R. White, Reference electronic structure calculations in one dimension, Physical Chemistry Chemical Physics14, 8581 (2012)
work page 2012
-
[72]
M. Motta, D. M. Ceperley, G. K.-L. Chan, J. A. Gomez, E. Gull, S. Guo, C. A. Jim´ enez-Hoyos, T. N. Lan, J. Li, F. Ma, A. J. Millis, N. V. Prokof’ev, U. Ray, G. E. Scuseria, S. Sorella, E. M. Stoudenmire, Q. Sun, I. S. Tupitsyn, S. R. White, D. Zgid, and S. Zhang (Simons Collaboration on the Many-Electron Problem), Towards the solution of the many-electro...
work page 2017
-
[73]
M. Dubeckyy, L. Mitas, and P. Jurecka, Noncovalent interactions by quantum monte carlo, Chemical reviews 116, 5188 (2016)
work page 2016
-
[74]
M. J. Gillan, D. Alf` e, and A. Michaelides, Perspective: How good is dft for water?, The Journal of Chemical Physics144, 130901 (2016)
work page 2016
-
[75]
W. M. C. Foulkes, L. Mitas, R. J. Needs, and G. Rajagopal, Quantum monte carlo simulations of solids, Reviews of Modern Physics73, 33 (2001)
work page 2001
- [76]
-
[77]
A. Szabo and N. S. Ostlund,Modern quantum chemistry: introduction to advanced electronic structure theory(Courier Corporation, 2012)
work page 2012
-
[78]
K. E. Riley, M. Pitonak, P. Jurecka, and P. Hobza, Stabilization and structure calculations for noncovalent interactions in extended molecular systems based on wave function and density functional theories, Chemical Reviews110, 5023 (2010), pMID: 20486691, http://dx.doi.org/10.1021/cr1000173
- [79]
- [80]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.