pith. machine review for the scientific record. sign in

arxiv: 2605.01593 · v1 · submitted 2026-05-02 · ❄️ cond-mat.other

Recognition: unknown

First-principles simulation of shocked H-He mixture along the principal Hugoniot

Ammar A. Ellaboudy , Valentin V. Karasiev , S. X. Hu
Authors on Pith no claims yet

Pith reviewed 2026-05-10 14:55 UTC · model grok-4.3

classification ❄️ cond-mat.other
keywords H-He mixtureshock Hugoniotreflectivitydemixingab initio molecular dynamicswarm dense matterJupiter interior
0
0 comments X

The pith

Ab initio simulations of shocked H-He mixtures match experimental reflectivity but show no discontinuities indicative of demixing.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper performs first-principles molecular dynamics simulations of an H-He mixture with 11 percent helium under shock compression to replicate conditions inside Jupiter. It obtains ionic configurations and the equation of state with the Tr2SCANL meta-GGA functional, then computes optical properties using the RS-KDT0 range-separated thermal hybrid functional within the Kubo-Greenwood formalism. The resulting reflectivity agrees well with laser-shock experiments across the temperature range, yet exhibits no jumps at the elevated temperatures where demixing had been inferred from data. A reader would care because this finding directly questions whether reflectivity measurements alone can reliably detect helium-hydrogen phase separation in planetary interiors.

Core claim

The calculated reflectivity for the mixed H-He system along the principal Hugoniot shows overall good agreement with experimental measurements; however, no discontinuity is observed at elevated temperatures, and the mixed-system predictions remain consistent with the data in the temperature range where the mixture is inferred to be demixed.

What carries the argument

Ab initio molecular dynamics with the Tr2SCANL meta-GGA functional for equation of state and ionic structure, combined with the RS-KDT0 range-separated thermal hybrid functional for optical properties evaluated via the Kubo-Greenwood formalism.

Load-bearing premise

The chosen exchange-correlation functionals accurately describe both thermodynamic and optical properties of the H-He mixture without systematic bias in the warm-dense-matter regime.

What would settle it

An experimental measurement that clearly resolves a sharp discontinuity in reflectivity at the temperatures and pressures where demixing was previously inferred would contradict the simulation results showing continuous reflectivity for the mixed phase.

Figures

Figures reproduced from arXiv: 2605.01593 by Ammar A. Ellaboudy, S. X. Hu, Valentin V. Karasiev.

Figure 1
Figure 1. Figure 1: FIG. 1. Shock pressure as a function of density along the prin [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Temperature as a function of pressure along the prin [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Reflectivity for a shock-compressed H–He mixture [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
read the original abstract

Recent laser-shock experiments on an H--He mixture containing 11~$\%$ helium (atomic fraction) have suggested the presence of an immiscibility region inside Jupiter. Reflectivity measurements were used as the primary diagnostic of H--He demixing, with discontinuities in the optical reflectivity proposed as a signature of phase separation under conditions relevant to Jupiter's interior. Here, we investigate shock-compressed H--He using \textit{ab initio} molecular dynamics simulations with optical properties evaluated within the Kubo--Greenwood formalism. The equation of state and ionic configurations were obtained using the thermal Tr$^2$SCANL meta-GGA exchange--correlation (XC) functional, while optical properties were computed using the recently developed RS-KDT0 range-separated thermal hybrid XC, which provides state-of-the-art accuracy for band-gap predictions in the warm dense matter regime. The calculated reflectivity shows overall good agreement with experimental measurements; however, no discontinuity is observed at elevated temperatures. Moreover, the reflectivity predictions for the mixed system are consistent with the experimental measurements in the temperature range where the mixture is inferred to be demixed. These results suggest that reflectivity alone may not provide a unique or sensitive diagnostic of H-He demixing at low helium concentrations under these conditions.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper reports ab initio MD simulations of the principal Hugoniot for an 11% He (atomic) H-He mixture, using Tr²SCANL meta-GGA for the EOS and ionic configurations and the RS-KDT0 range-separated thermal hybrid for Kubo-Greenwood optical conductivity. The calculated reflectivity agrees overall with laser-shock experiments but exhibits no discontinuity along the Hugoniot; the authors conclude that reflectivity is not a unique or sensitive diagnostic of H-He demixing at low helium concentrations under Jupiter-interior conditions.

Significance. If the simulations are unbiased, the work would weaken the experimental claim that reflectivity discontinuities signal H-He phase separation, with implications for models of giant-planet interiors. The direct, parameter-free comparison to independent reflectivity data and the adoption of state-of-the-art thermal hybrid functionals for both EOS and optics constitute clear strengths.

major comments (2)
  1. [§2] §2 (Computational Methods, Kubo-Greenwood subsection): The choice of RS-KDT0 is motivated by its reported band-gap accuracy in WDM, yet the manuscript provides no direct benchmark of reflectivity or optical conductivity for H-He mixtures at the relevant densities and temperatures (1000–6000 K). A systematic error in the gap or dipole matrix elements could shift the entire reflectivity curve or mask a weak discontinuity, directly affecting the central claim that the mixed-phase reflectivity matches experiment in the inferred demixing regime.
  2. [§3] §3 (Results, reflectivity along the Hugoniot): No convergence tests, error bars, or sensitivity analysis are reported for k-point sampling, supercell size, or number of bands in the Kubo-Greenwood calculations. Without these, the absence of a discontinuity and the quantitative agreement with experiment cannot be assessed for robustness, which is load-bearing for the conclusion that reflectivity alone is insensitive to demixing.
minor comments (2)
  1. [Abstract] The abstract states that 'no discontinuity is observed at elevated temperatures' without specifying the temperature window or the expected size of a discontinuity that would be detectable given the reported precision.
  2. [Figures] Figure captions and axis labels for the reflectivity plots should explicitly state the wavelength or photon energy at which reflectivity is evaluated, as this is essential for comparing to the experimental diagnostic.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of our manuscript and for highlighting important points regarding the computational methodology and the robustness of our results. We address each major comment below and indicate the revisions we will make to strengthen the paper.

read point-by-point responses

  1. Referee: [§2] §2 (Computational Methods, Kubo-Greenwood subsection): The choice of RS-KDT0 is motivated by its reported band-gap accuracy in WDM, yet the manuscript provides no direct benchmark of reflectivity or optical conductivity for H-He mixtures at the relevant densities and temperatures (1000–6000 K). A systematic error in the gap or dipole matrix elements could shift the entire reflectivity curve or mask a weak discontinuity, directly affecting the central claim that the mixed-phase reflectivity matches experiment in the inferred demixing regime.

    Authors: We acknowledge that the manuscript does not include a direct benchmark of RS-KDT0 for reflectivity or optical conductivity specifically in H-He mixtures at the cited conditions. The functional was selected on the basis of its documented accuracy for band gaps across a range of warm-dense-matter systems in the literature. While a dedicated benchmark for this mixture would be ideal, the quantitative agreement between our computed reflectivity and the independent experimental data provides supporting evidence for the reliability of the approach. A uniform systematic error in the gap would primarily shift the absolute reflectivity level rather than eliminate a discontinuity arising from phase separation, which would involve distinct changes in local electronic structure and composition. We will revise the methods section to include a more explicit discussion of the functional's applicability, expected uncertainties, and the rationale for why such an error is unlikely to mask the signature we are testing for. revision: partial

  2. Referee: [§3] §3 (Results, reflectivity along the Hugoniot): No convergence tests, error bars, or sensitivity analysis are reported for k-point sampling, supercell size, or number of bands in the Kubo-Greenwood calculations. Without these, the absence of a discontinuity and the quantitative agreement with experiment cannot be assessed for robustness, which is load-bearing for the conclusion that reflectivity alone is insensitive to demixing.

    Authors: We agree that explicit convergence tests and error estimates are necessary to substantiate the robustness of the reported reflectivity curves and the absence of discontinuities. Although such tests were performed during the study (using supercells of 128–256 atoms, dense k-point sampling including the Gamma point for large cells, and inclusion of all valence plus a sufficient number of conduction bands), they were not presented in detail. We will add a dedicated appendix (or supplementary section) containing the convergence data, showing that reflectivity values are converged to within approximately 5% in the relevant photon-energy range and that the lack of discontinuity persists across the tested parameters. Error bars derived from these tests will also be incorporated into the main-text figures where appropriate. revision: yes

Circularity Check

0 steps flagged

No significant circularity; direct ab initio vs. independent experiment comparison

full rationale

The paper computes EOS and ionic trajectories with Tr²SCANL meta-GGA, then optical conductivity via Kubo-Greenwood using the RS-KDT0 hybrid, and compares the resulting reflectivity directly to external laser-shock measurements on the 11% He mixture. No target observable (reflectivity or discontinuity) is used to fit any parameter, define any functional, or close any equation; the absence of discontinuity is an output of the simulation, not an input. Prior citations to functional development are present but are not load-bearing for the central claim, which rests on the numerical comparison to independent data rather than on any self-referential definition or fitted prediction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on standard approximations of density-functional theory and the accuracy of two specific exchange-correlation functionals whose performance was established in earlier benchmark studies rather than derived here.

axioms (2)
  • standard math Born-Oppenheimer approximation separating electronic and nuclear motion
    Invoked implicitly in all ab initio molecular dynamics simulations described.
  • domain assumption Validity of the Kubo-Greenwood formalism for computing optical conductivity in the warm-dense-matter regime
    Used to obtain reflectivity from the simulated electronic states.

pith-pipeline@v0.9.0 · 5531 in / 1224 out tokens · 70604 ms · 2026-05-10T14:55:17.439343+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

48 extracted references · 2 canonical work pages

  1. [1]

    S. M. Wahl, W. B. Hubbard, B. Militzer, T. Guillot, Y. Miguel, N. Movshovitz, Y. Kaspi, R. Helled, D. Reese, E. Galanti, S. Levin, J. E. Connerney, and S. J. Bolton, Comparingjupiterinteriorstructuremodelstojunograv- ity measurements and the role of a dilute core, Geophys- ical Research Letters44, 4649 (2017)

    2017

  2. [2]

    Debras and G

    F. Debras and G. Chabrier, New models of jupiter in the context of juno and galileo, The Astrophysical Journal 872, 100 (2019)

    2019

  3. [3]

    C. R. Mankovich and J. J. Fortney, Evidence for a di- chotomy in the interior structures of jupiter and saturn from helium phase separation, The Astrophysical Journal 889, 51 (2020)

    2020

  4. [4]

    Howard, S

    S. Howard, S. Müller, and R. Helled, Evolution of jupiter and saturn with helium rain, Astronomy & Astrophysics 689, A15 (2024)

    2024

  5. [5]

    A. Sur, R. Tejada Arevalo, Y. Su, and A. Burrows, Si- multaneous evolutionary fits for jupiter and saturn incor- porating fuzzy cores, The Astrophysical Journal Letters 980, L5 (2025)

    2025

  6. [6]

    Tejada Arevalo, Y

    R. Tejada Arevalo, Y. Su, A. Sur, and A. Burrows, Equa- tions of state, thermodynamics, and miscibility curves for jovian planet and giant exoplanet evolutionary mod- els, The Astrophysical Journal Supplement Series274, 34 (2024)

    2024

  7. [7]

    D. J. Stevenson, Thermodynamics and phase separation of dense fully ionized hydrogen-helium fluid mixtures, Phys. Rev. B12, 3999 (1975)

    1975

  8. [8]

    Lorenzen, B

    W. Lorenzen, B. Holst, and R. Redmer, Demixing of hydrogen and helium at megabar pressures, Phys. Rev. Lett.102, 115701 (2009)

    2009

  9. [9]

    M. A. Morales, S. Hamel, K. Caspersen, and E. Schwe- gler, Hydrogen-helium demixing from first principles: From diamond anvil cells to planetary interiors, Phys. Rev. B87, 174105 (2013)

    2013

  10. [10]

    Schöttler and R

    M. Schöttler and R. Redmer, Ab initio calculation of the miscibility diagram for hydrogen-helium mixtures, Phys. Rev. Lett.120, 115703 (2018)

    2018

  11. [11]

    R. C. Clay, M. Holzmann, D. M. Ceperley, and M. A. Morales, Benchmarking density functionals for hydrogen- helium mixtures with quantum monte carlo: Energetics, pressures, and forces, Phys. Rev. B93, 035121 (2016)

    2016

  12. [12]

    V. V. Karasiev, S. X. Hu, J. P. Hinz, R. M. N. Goshadze, S. Zhang, A. Bergermann, and R. Redmer, Inhibiting conduction by he mixing in interiors of jupiter and saturn (2026), arXiv:2601.23152 [astro-ph.EP]

    work page arXiv 2026

  13. [13]

    Chang, B

    X. Chang, B. Chen, Q. Zeng, H. Wang, K. Chen, Q. Tong, X. Yu, D. Kang, S. Zhang, F. Guo, Y. Hou, Z. Zhao, Y. Yao, Y. Ma, and J. Dai, Theoretical evidence of h–he demixing under jupiter and saturn conditions, Nature Communications15, 8543 (2024)

    2024

  14. [14]

    Brygoo, P

    S. Brygoo, P. Loubeyre, M. Millot, J. R. Rygg, P. M. Celliers, J. H. Eggert, R. Jeanloz, and G. W. Collins, Evi- denceofhydrogen-heliumimmiscibilityatjupiter-interior conditions, Nature (London)593, 10.1038/s41586-021- 03516-0 (2021)

    work page doi:10.1038/s41586-021- 2021

  15. [15]

    Soubiran, S

    F. Soubiran, S. Mazevet, C. Winisdoerffer, and G. Chabrier, Optical signature of hydrogen-helium demixing at extreme density-temperature conditions, Phys. Rev. B87, 165114 (2013)

    2013

  16. [16]

    Loubeyre, R

    P. Loubeyre, R. Le Toullec, and J. P. Pinceaux, Binary phase diagrams ofh 2-he mixtures at high temperature and high pressure, Phys. Rev. B36, 3723 (1987)

    1987

  17. [17]

    J. A. Schouten, A. de Kuijper, and J. P. J. Michels, Crit- ical line of he-h2 up to 2500 k and the influence of at- traction on fluid-fluid separation, Phys. Rev. B44, 6630 (1991)

    1991

  18. [18]

    Kohn and L

    W. Kohn and L. J. Sham, Self-consistent equations in- cluding exchange and correlation effects, Phys. Rev.140, A1133 (1965)

    1965

  19. [19]

    N. D. Mermin, Thermal properties of the inhomogeneous electron gas, Phys. Rev.137, A1441 (1965)

    1965

  20. [20]

    J. P. Perdew and K. Schmidt, Jacob’s ladder of density functional approximations for the exchange-correlation energy, AIP Conference Proceedings577, 1 (2001)

    2001

  21. [21]

    J. P. Perdew, K. Burke, and Y. Wang, Generalized gra- dient approximation for the exchange-correlation hole of a many-electron system, Physical Review B54, 16533 (1996)

    1996

  22. [22]

    V. V. Karasiev, L. Calderín, and S. B. Trickey, Im- portance of finite-temperature exchange correlation for warm dense matter calculations, Phys. Rev. E93, 063207 (2016)

    2016

  23. [23]

    V. V. Karasiev, T. Sjostrom, J. Dufty, and S. B. Trickey, Accurate homogeneous electron gas exchange-correlation free energy for local spin-density calculations, Phys. Rev. Lett.112, 076403 (2014)

    2014

  24. [24]

    V. V. Karasiev, J. W. Dufty, and S. B. Trickey, Nonem- pirical semilocal free-energy density functional for matter under extreme conditions, Phys. Rev. Lett.120, 076401 (2018)

    2018

  25. [25]

    V. V. Karasiev, D. I. Mihaylov, and S. X. Hu, Meta- ggaexchange-correlationfreeenergydensityfunctionalto increase the accuracy of warm dense matter simulations, Phys. Rev. B105, L081109 (2022)

    2022

  26. [26]

    K. P. Hilleke, V. V. Karasiev, S. B. Trickey, R. M. N. Goshadze, and S. X. Hu, Fully thermal meta-gga ex- change correlation free-energy density functional, Phys. Rev. Mater.9, L050801 (2025)

    2025

  27. [27]

    D. I. Mihaylov, V. V. Karasiev, and S. X. Hu, Ther- mal hybrid exchange-correlation density functional for improving the description of warm dense matter, Phys. 10 Rev. B101, 245141 (2020)

    2020

  28. [28]

    A. A. Ellaboudy, V. V. Karasiev, D. I. Mihaylov, K. P. Hilleke, and S. X. Hu, Range-separated thermal hy- brid exchange-correlation density functional for accurate band-gap calculations of warm dense matter, Phys. Rev. B112, 155154 (2025)

    2025

  29. [29]

    J. Sun, A. Ruzsinszky, and J. P. Perdew, Strongly con- strained and appropriately normed semilocal density functional, Phys. Rev. Lett.115, 036402 (2015)

    2015

  30. [30]

    A. P. Bartok and J. R. Yates, Regularized scan func- tional, The Journal of Chemical Physics150, 161101 (2019)

    2019

  31. [31]

    J. P. Perdew, W. Yang, K. Burke, Z. Yang, E. K. U. Gross, M. Scheffler, G. E. Scuseria, T. M. Henderson, I. Y. Zhang, A. Ruzsinszky, H. Peng, J. Sun, E. Trushin, and A. Görling, Understanding band gaps of solids in generalized kohn–sham theory, Proceedings of the Na- tional Academy of Sciences114, 2801 (2017)

    2017

  32. [32]

    Mori-Sánchez, A

    P. Mori-Sánchez, A. J. Cohen, and W. Yang, Localization and delocalization errors in density functional theory and implications for band-gap prediction, Phys. Rev. Lett. 100, 146401 (2008)

    2008

  33. [33]

    A. J. Cohen, P. Mori-Sánchez, and W. Yang, Challenges for density functional theory, Chemical Reviews112, 289 (2012)

    2012

  34. [34]

    J. P. Perdew, M. Ernzerhof, and K. Burke, Rationale for mixing exact exchange with density functional approx- imations, The Journal of Chemical Physics105, 9982 (1996)

    1996

  35. [35]

    J. Heyd, G. E. Scuseria, and M. Ernzerhof, Hybrid func- tionals based on a screened coulomb potential, The Jour- nal of Chemical Physics118, 8207 (2003)

    2003

  36. [36]

    Kubo, Statistical-mechanical theory of irreversible processes

    R. Kubo, Statistical-mechanical theory of irreversible processes. i. general theory and simple applications to magnetic and conduction problems, Journal of the Phys- ical Society of Japan12, 570 (1957)

    1957

  37. [37]

    D. A. Greenwood, The boltzmann equation in the the- ory of electrical conduction in metals, Proceedings of the Physical Society71, 585 (1958)

    1958

  38. [38]

    S. X. Hu, T. R. Boehly, and L. A. Collins, Properties of warm dense polystyrene plasmas along the principal hugoniot, Phys. Rev. E89, 063104 (2014)

    2014

  39. [39]

    Blanchet, V

    A. Blanchet, V. Recoules, F. Soubiran, and M. Tacu, Computation of transport properties of warm dense mat- ter using abinit, Physics of Plasmas31, 062703 (2024)

    2024

  40. [40]

    Kresse and J

    G. Kresse and J. Furthmüller, Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set, Phys. Rev. B54, 11169 (1996)

    1996

  41. [41]

    Kresse and J

    G. Kresse and J. Furthmüller, Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set, Computational Materials Science 6, 15 (1996)

    1996

  42. [42]

    P. E. Blöchl, Projector augmented-wave method, Phys. Rev. B50, 17953 (1994)

    1994

  43. [43]

    Kresse and D

    G. Kresse and D. Joubert, From ultrasoft pseudopoten- tials to the projector augmented-wave method, Phys. Rev. B59, 1758 (1999)

    1999

  44. [44]

    Baldereschi, Mean-value point in the brillouin zone, Phys

    A. Baldereschi, Mean-value point in the brillouin zone, Phys. Rev. B7, 5212 (1973)

    1973

  45. [45]

    Zhang and S

    S. Zhang and S. X. Hu, Species separation and hydro- gen streaming upon shock release from polystyrene under inertial confinement fusion conditions, Phys. Rev. Lett. 125, 105001 (2020)

    2020

  46. [46]

    Lorenzen, B

    W. Lorenzen, B. Holst, and R. Redmer, Metallization in hydrogen-helium mixtures, Phys. Rev. B84, 235109 (2011)

    2011

  47. [47]

    Militzer, Equation of state calculations of hydrogen- helium mixtures in solar and extrasolar giant planets, Phys

    B. Militzer, Equation of state calculations of hydrogen- helium mixtures in solar and extrasolar giant planets, Phys. Rev. B87, 014202 (2013)

    2013

  48. [48]

    Hamel, M

    S. Hamel, M. A. Morales, and E. Schwegler, Signature of helium segregation in hydrogen-helium mixtures, Phys. Rev. B84, 165110 (2011)

    2011