pith. sign in

arxiv: 2512.19593 · v2 · submitted 2025-12-22 · ⚛️ nucl-th · astro-ph.HE· nucl-ex

Neutron star crust and outer core equation of state from chiral effective field theory with quantified uncertainties

Pith reviewed 2026-05-16 20:14 UTC · model grok-4.3

classification ⚛️ nucl-th astro-ph.HEnucl-ex
keywords neutron star crustchiral effective field theoryequation of stateGaussian processnuclear matterbeta equilibriummany-body perturbation theoryN3LO
0
0 comments X

The pith

Chiral effective field theory at N3LO combined with a two-dimensional Gaussian process yields consistent equation-of-state tables for neutron-star inner crusts with quantified uncertainties.

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

The work trains a two-dimensional Gaussian process on many-body perturbation theory calculations of asymmetric nuclear matter using chiral two- and three-nucleon forces through N3LO. This interpolation supplies the energy per particle, pressure, and chemical potentials of beta-equilibrated neutron-rich matter together with EFT truncation uncertainties up to twice saturation density. The resulting phase diagram tracks the progression from neutron drip through proton drip into the uniform phase while incorporating surface and Coulomb corrections. Equations of state for the inner crust are then assembled that remain consistent with the uniform-matter results at N3LO. These tables supply controlled nuclear-physics inputs for neutron-star models without arbitrary matching procedures at the crust-core boundary.

Core claim

We construct EOSs for the inner crust of neutron stars that are consistent with the chiral EFT results for uniform matter at N3LO by training a two-dimensional Gaussian process on many-body perturbation theory results for asymmetric nuclear matter, thereby providing efficient access to thermodynamic quantities and EFT truncation uncertainties across the relevant density and proton-fraction range.

What carries the argument

Two-dimensional Gaussian process trained on many-body perturbation theory results from chiral two- and three-nucleon interactions up to N3LO, which interpolates the energy per particle and supplies thermodynamic derivatives together with EFT truncation uncertainties.

If this is right

  • The inner-crust EOS tables include propagated EFT truncation uncertainties for all thermodynamic quantities.
  • The phase diagram maps the neutron-drip, proton-drip, and uniform-matter regions with surface and Coulomb corrections included.
  • Chemical potentials and pressure in beta equilibrium are available with quantified uncertainties up to twice saturation density.
  • The same Gaussian-process representation can be used to evaluate any thermodynamic derivative needed for neutron-star structure calculations.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The tabulated EOS could be inserted directly into neutron-star cooling or glitch models to propagate nuclear-theory uncertainties into observable predictions.
  • The interpolation technique may be applied at higher densities to generate consistent outer-core EOS segments without new many-body calculations.
  • Observational constraints on neutron-star radii or tidal deformability could be reinterpreted as direct tests of the N3LO uncertainty bands reported here.

Load-bearing premise

The two-dimensional Gaussian process accurately interpolates the many-body perturbation theory results and faithfully represents the EFT truncation uncertainties across the full range of densities and proton fractions relevant to the inner crust.

What would settle it

A calculation of the pressure or chemical potential at the crust-core transition density that lies outside the reported N3LO uncertainty band would falsify the consistency of the constructed inner-crust EOS with the uniform-matter results.

Figures

Figures reproduced from arXiv: 2512.19593 by A. Schwenk, H. G\"ottling, K. Hebeler, L. Hoff.

Figure 1
Figure 1. Figure 1: FIG. 1. Expansion coe [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Expansion coe [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Corrections to the energy per particle from LO to N [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Optimized marginal variance ¯c [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Model-checking diagnostics for the GP constructed with the training set [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Same as Fig [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Energy per particle from NLO to N [PITH_FULL_IMAGE:figures/full_fig_p007_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Energy per particle (left panel), pressure (middle panel), and proton fraction (right panel) for matter in [PITH_FULL_IMAGE:figures/full_fig_p008_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Chemical potentials of neutrons [PITH_FULL_IMAGE:figures/full_fig_p008_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Phase coexistence in neutron star matter, obtained using the [PITH_FULL_IMAGE:figures/full_fig_p010_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Pressure of neutron star matter as function of baryon density [PITH_FULL_IMAGE:figures/full_fig_p011_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13. Total proton fraction within the Wigner-Seitz cell as a func [PITH_FULL_IMAGE:figures/full_fig_p012_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14. Expansion coe [PITH_FULL_IMAGE:figures/full_fig_p014_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: FIG. 15. Expansion coe [PITH_FULL_IMAGE:figures/full_fig_p014_15.png] view at source ↗
read the original abstract

We study the order-by-order expansion of the energy per particle of asymmetric nuclear matter up to twice saturation density in chiral effective field theory (EFT) within a Bayesian framework. For this, we develop a two-dimensional Gaussian process (2D GP) that is trained using many-body perturbation theory results based on chiral two- and three-nucleon interactions from leading to next-to-next-to-next-to-leading order (N$^3$LO). This allows for an efficient evaluation of the equation of state (EOS) and thermodynamic derivatives with EFT truncation uncertainties. After benchmarking our 2D GP against Bayesian uncertainties for pure neutron matter and symmetric matter, we study the energy per particle, pressure, and chemical potentials of neutron star matter in $\beta$-equilibrium including EFT uncertainties. We investigate the phase diagram of neutron-rich matter from neutron- to proton-drip and to the uniform phase, including surface and Coulomb corrections. Based on this, we construct EOSs for the inner crust of neutron stars that are consistent with the chiral EFT results for uniform matter at N$^3$LO.

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

1 major / 2 minor

Summary. The manuscript develops a two-dimensional Gaussian process trained on many-body perturbation theory calculations from chiral EFT interactions up to N³LO to model the energy per particle, pressure, and chemical potentials of asymmetric nuclear matter up to twice saturation density. It benchmarks the GP against Bayesian uncertainties for pure neutron matter and symmetric matter, extends the approach to β-equilibrated neutron star matter, investigates the phase diagram from neutron drip through the inner crust including surface and Coulomb corrections, and constructs inner-crust EOS tables claimed to be consistent with the N³LO uniform-matter results.

Significance. If the central construction and uncertainty propagation hold, the work supplies a practical set of crust and outer-core EOS tables with quantified EFT truncation uncertainties that can be directly used in neutron-star structure calculations. The GP approach enables efficient evaluation of thermodynamic derivatives while propagating truncation errors, which is a clear methodological advance over pointwise EFT calculations for inhomogeneous phases.

major comments (1)
  1. [phase diagram and inner-crust EOS construction] The phase-diagram and inner-crust construction section: surface and Coulomb corrections are added via standard liquid-drop or Thomas-Fermi prescriptions whose parameters are not varied inside the Bayesian framework and whose uncertainties are not propagated into the final EOS bands. Because the location of the neutron-drip to uniform transition and the pressure jump at that point are sensitive to these terms, the advertised consistency of the complete inner-crust EOS with N³LO uniform-matter results is not fully demonstrated; the reported EFT truncation bands may be smaller than the unquantified model dependence introduced by the surface/Coulomb modeling.
minor comments (2)
  1. [benchmarking] The benchmarking subsection would benefit from an explicit table comparing GP-predicted central values and uncertainty bands against the original MBPT points at a few representative densities and proton fractions.
  2. [methods] Notation for the two-dimensional GP kernel and the precise definition of the truncation-error model should be collected in one place to avoid scattered references across the methods section.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. We address the single major comment below and have revised the manuscript to clarify the scope of the uncertainty quantification and to include additional sensitivity information.

read point-by-point responses
  1. Referee: The phase-diagram and inner-crust construction section: surface and Coulomb corrections are added via standard liquid-drop or Thomas-Fermi prescriptions whose parameters are not varied inside the Bayesian framework and whose uncertainties are not propagated into the final EOS bands. Because the location of the neutron-drip to uniform transition and the pressure jump at that point are sensitive to these terms, the advertised consistency of the complete inner-crust EOS with N³LO uniform-matter results is not fully demonstrated; the reported EFT truncation bands may be smaller than the unquantified model dependence introduced by the surface/Coulomb modeling.

    Authors: We agree that the surface and Coulomb corrections rely on fixed-parameter liquid-drop prescriptions whose uncertainties are not varied or propagated inside the Bayesian framework used for the chiral EFT truncation errors. This is a genuine limitation of the present work. The advertised consistency refers specifically to the fact that the bulk energy per particle, pressure, and chemical potentials in the uniform phase are taken directly from the N³LO chiral EFT results (via the 2D GP), while the inhomogeneous corrections are added on top using standard, widely employed prescriptions. In the revised manuscript we have (i) added an explicit statement in the phase-diagram section clarifying the scope of the reported uncertainty bands, (ii) performed a limited sensitivity analysis by varying the surface tension coefficient within its conventional range and reported the resulting shifts in the neutron-drip transition density and pressure, and (iii) included these variations as supplementary uncertainty estimates in the released EOS tables. These additions allow readers to judge the relative size of the unquantified model dependence versus the propagated EFT truncation uncertainties. revision: partial

Circularity Check

0 steps flagged

No significant circularity; GP interpolation of independent MBPT inputs remains non-reductive

full rationale

The central construction trains a 2D Gaussian process on many-body perturbation theory results computed from chiral two- and three-nucleon interactions at orders up to N3LO; these MBPT calculations constitute external, order-by-order inputs. The GP then supplies uniform-matter energies, pressures and chemical potentials (with truncation uncertainties) that are inserted into the phase-diagram construction together with standard liquid-drop surface and Coulomb terms. No equation or claim equates a derived EOS quantity to a fitted parameter by construction, nor does any load-bearing step rest on a self-citation whose content is itself unverified. The advertised consistency is explicitly restricted to the uniform-matter sector at N3LO, leaving the external corrections as a separate modeling choice whose uncertainties are not folded into the reported bands. This yields a minor self-citation exposure at most (score 2) while the derivation chain stays self-contained against the stated benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the validity of chiral EFT up to N3LO for densities up to twice saturation, the accuracy of many-body perturbation theory for the training data, and the ability of the 2D GP to represent truncation errors. No new particles or forces are introduced.

axioms (2)
  • domain assumption Chiral effective field theory truncated at N3LO is a reliable description of nuclear forces up to twice saturation density
    Invoked throughout the study of the order-by-order expansion
  • domain assumption Many-body perturbation theory provides sufficiently accurate training data for the Gaussian process
    Basis for the 2D GP training

pith-pipeline@v0.9.0 · 5507 in / 1356 out tokens · 24987 ms · 2026-05-16T20:14:49.282558+00:00 · methodology

discussion (0)

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

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Astrophysics equation of state inference with Bayesian chiral effective field theory uncertainties

    nucl-th 2026-05 unverdicted novelty 5.0

    Bayesian EOS inference with χEFT uncertainty priors and LIGO/NICER data yields posteriors consistent with prior work, a stiffening above 3n0, negligible pQCD impact, and an inferred symmetry-energy slope L of 42.6-56.7 MeV.

  2. Conformal prediction for uncertainties in the neutron star equation of state

    nucl-th 2026-04 unverdicted novelty 4.0

    Conformalized quantile regression applied post hoc to neutron star posterior samples yields reliable uncertainty bands validated by empirical coverage studies.

Reference graph

Works this paper leans on

45 extracted references · 45 canonical work pages · cited by 2 Pith papers

  1. [1]

    Yasin, S

    H. Yasin, S. Schäfer, A. Arcones, and A. Schwenk, Equation of state effects in core-collapse supernovae, Phys. Rev. Lett.124, 092701 (2020)

  2. [2]

    A. S. Schneider, L. F. Roberts, C. D. Ott, and E. O’Connor, Equation of state effects in the core collapse of a 20-M ⊙ star, Phys. Rev. C100, 055802 (2019)

  3. [3]

    J. M. Lattimer, Neutron Stars and the Nuclear Matter Equation of State, Annu. Rev. Nucl. Part. Sci.71, 433 (2021)

  4. [4]

    Jacobi, F

    M. Jacobi, F. M. Guercilena, S. Huth, G. Ricigliano, A. Ar- cones, and A. Schwenk, Effects of nuclear matter properties in neutron star mergers, Mon. Not. Roy. Astron. Soc.527, 8812 (2023)

  5. [5]

    Epelbaum, H.-W

    E. Epelbaum, H.-W. Hammer, and U.-G. Meißner, Modern the- ory of nuclear forces, Rev. Mod. Phys.81, 1773 (2009)

  6. [6]

    Machleidt and D

    R. Machleidt and D. R. Entem, Chiral effective field theory and nuclear forces, Phys. Rep.503, 1 (2011)

  7. [7]

    Hebeler and A

    K. Hebeler and A. Schwenk, Chiral three-nucleon forces and neutron matter, Phys. Rev. C82, 014314 (2010)

  8. [8]

    I. Tews, T. Krüger, K. Hebeler, and A. Schwenk, Neutron Mat- ter at Next-to-Next-to-Next-to-Leading Order in Chiral Effec- tive Field Theory, Phys. Rev. Lett.110, 032504 (2013)

  9. [9]

    J. W. Holt, N. Kaiser, and W. Weise, Nuclear chiral dynamics and thermodynamics, Prog. Part. Nucl. Phys.73, 35 (2013)

  10. [10]

    Carbone, A

    A. Carbone, A. Polls, and A. Rios, Symmetric nuclear matter with chiral three-nucleon forces in the self-consistent Green’s functions approach, Phys. Rev. C88, 044302 (2013)

  11. [11]

    Hagen, T

    G. Hagen, T. Papenbrock, A. Ekström, K. Wendt, G. Baardsen, S. Gandolfi, M. Hjorth-Jensen, and C. J. Horowitz, Coupled- cluster calculations of nucleonic matter, Phys. Rev. C89, 014319 (2014)

  12. [12]

    Coraggio, J

    L. Coraggio, J. W. Holt, N. Itaco, R. Machleidt, L. E. Mar- cucci, and F. Sammarruca, Nuclear-matter equation of state with consistent two- and three-body perturbative chiral inter- actions, Phys. Rev. C89, 044321 (2014)

  13. [13]

    Wellenhofer, J

    C. Wellenhofer, J. W. Holt, N. Kaiser, and W. Weise, Nu- clear thermodynamics from chiral low-momentum interactions, Phys. Rev. C89, 064009 (2014)

  14. [14]

    Wellenhofer, J

    C. Wellenhofer, J. W. Holt, and N. Kaiser, Thermodynamics of isospin-asymmetric nuclear matter from chiral effective field theory, Phys. Rev. C92, 015801 (2015)

  15. [15]

    J. E. Lynn, I. Tews, J. Carlson, S. Gandolfi, A. Gezerlis, K. E. Schmidt, and A. Schwenk, Chiral Three-Nucleon Interactions in Light Nuclei, Neutron-αScattering, and Neutron Matter, Phys. Rev. Lett.116, 062501 (2016)

  16. [16]

    Drischler, K

    C. Drischler, K. Hebeler, and A. Schwenk, Asymmetric nu- clear matter based on chiral two- and three-nucleon interac- tions, Phys. Rev. C93, 054314 (2016)

  17. [17]

    Ekström, G

    A. Ekström, G. Hagen, T. D. Morris, T. Papenbrock, and P. D. Schwartz,∆isobars and nuclear saturation, Phys. Rev. C97, 024332 (2018)

  18. [18]

    Drischler, K

    C. Drischler, K. Hebeler, and A. Schwenk, Chiral interactions up to next-to-next-to-next-to-leading order and nuclear satura- tion, Phys. Rev. Lett.122, 042501 (2019)

  19. [19]

    Carbone and A

    A. Carbone and A. Schwenk, Ab initio constraints on ther- mal effects of the nuclear equation of state, Phys. Rev. C100, 025805 (2019)

  20. [20]

    B.-N. Lu, N. Li, S. Elhatisari, D. Lee, J. E. Drut, T. A. Lähde, E. Epelbaum, and U.-G. Meißner, Ab Initio Nuclear Thermo- dynamics, Phys. Rev. Lett.125, 192502 (2020)

  21. [21]

    Keller, C

    J. Keller, C. Wellenhofer, K. Hebeler, and A. Schwenk, Neutron matter at finite temperature based on chiral effective field theory interactions, Phys. Rev. C103, 055806 (2021)

  22. [22]

    Keller, K

    J. Keller, K. Hebeler, and A. Schwenk, Nuclear Equation of State for Arbitrary Proton Fraction and Temperature Based on Chiral Effective Field Theory and a Gaussian Process Emulator, Phys. Rev. Lett.130, 072701 (2023)

  23. [23]

    Marino, W

    F. Marino, W. Jiang, and S. J. Novario, Diagrammatic ab initio methods for infinite nuclear matter with modern chiral interac- tions, Phys. Rev. C110, 054322 (2024)

  24. [24]

    F. Alp, Y . Dietz, K. Hebeler, and A. Schwenk, Equation of state and Fermi liquid properties of dense matter based on chi- ral effective field theory interactions, Phys. Rev. C112, 055802 (2025)

  25. [25]

    C. E. Rasmussen and C. K. I. Williams,Gaussian Processes for Machine Learning(MIT Press, Cambridge, 2005)

  26. [26]

    J. A. Melendez, R. J. Furnstahl, D. R. Phillips, M. T. Pratola, and S. Wesolowski, Quantifying Correlated Truncation Errors in Effective Field Theory, Phys. Rev. C100, 044001 (2019)

  27. [27]

    Drischler, R

    C. Drischler, R. J. Furnstahl, J. A. Melendez, and D. R. Phillips, How Well Do We Know the Neutron-Matter Equation of State at the Densities Inside Neutron Stars? A Bayesian Approach with Correlated Uncertainties, Phys. Rev. Lett.125, 202702 (2020)

  28. [28]

    Drischler, J

    C. Drischler, J. A. Melendez, R. J. Furnstahl, and D. R. Phillips, Quantifying uncertainties and correlations in the nuclear-matter equation of state, Phys. Rev. C102, 054315 (2020)

  29. [29]

    D. R. Entem, R. Machleidt, and Y . Nosyk, High-quality two- nucleon potentials up to fifth order of the chiral expansion, Phys. Rev. C96, 024004 (2017)

  30. [30]

    R. J. Furnstahl, N. Klco, D. R. Phillips, and S. Wesolowski, Quantifying truncation errors in effective field theory, Phys. Rev. C92, 024005 (2015)

  31. [31]

    Hebeler, Three-nucleon forces: Implementation and appli- cations to atomic nuclei and dense matter, Phys

    K. Hebeler, Three-nucleon forces: Implementation and appli- cations to atomic nuclei and dense matter, Phys. Rep.890, 1 (2021)

  32. [32]

    L. S. Bastos and A. O’Hagan, Diagnostics for Gaussian Process Emulators, Technometrics51, 425 (2009)

  33. [33]

    G. Baym, C. J. Pethick, and P. Sutherland, The Ground state of matter at high densities: Equation of state and stellar models, Astrophys. J.170, 299 (1970). 14 0.1 0.2 0.3 n [fm□3] □2 0 2 4 6 c x = 0 x = 0.5 c2 0.1 0.2 0.3 n [fm□3] c3 0.1 0.2 0.3 n [fm□3] c4 FIG. 14. Expansion coefficients at NLO (c 2, left), N 2LO (c 3, middle), and N 3LO (c 4, right pan...

  34. [34]

    Keller, K

    J. Keller, K. Hebeler, C. J. Pethick, and A. Schwenk, Neutron star matter as a dilute solution of protons in neutrons, Phys. Rev. Lett.132, 232701 (2024)

  35. [35]

    Chamel and P

    N. Chamel and P. Haensel, Physics of Neutron Star Crusts, Liv- ing Rev. Rel.11, 10 (2008)

  36. [36]

    Tews, Spectrum of shear modes in the neutron-star crust: Estimating the nuclear-physics uncertainties, Phys

    I. Tews, Spectrum of shear modes in the neutron-star crust: Estimating the nuclear-physics uncertainties, Phys. Rev. C95, 015803 (2017)

  37. [37]

    A. W. Steiner, Neutron star inner crust: Nuclear physics input, Phys. Rev. C77, 035805 (2008)

  38. [38]

    J. M. Lattimer, C. J. Pethick, D. G. Ravenhall, and D. Q. Lamb, Physical properties of hot, dense matter: The general case, Nucl. Phys. A432, 646 (1985)

  39. [39]

    M. Wang, W. J. Huang, F. G. Kondev, G. Audi, and S. Naimi, The AME 2020 atomic mass evaluation (II). Tables, graphs and references, Chin. Phys. C45, 030003 (2021)

  40. [40]

    C. J. Pethick and D. G. Ravenhall, Matter at large neutron ex- cess and the physics of neutron-star crusts, Annu. Rev. Nucl. Part. Sci.45, 429 (1995)

  41. [41]

    Haensel, A

    P. Haensel, A. Y . Potekhin, and D. G. Yakovlev,Neutron Stars 1: Equation of State and Structure(Springer, New York, USA, 2007)

  42. [42]

    Grams, R

    G. Grams, R. Somasundaram, J. Margueron, and S. Reddy, Properties of the neutron star crust: Quantifying and correlating uncertainties with improved nuclear physics, Phys. Rev. C105, 035806 (2022)

  43. [43]

    Hebeler, J

    K. Hebeler, J. M. Lattimer, C. J. Pethick, and A. Schwenk, Equation of state and neutron star properties constrained by nu- clear physics and observation, Astrophys. J.773, 11 (2013)

  44. [44]

    Svensson, A

    I. Svensson, A. Tichai, K. Hebeler, and A. Schwenk, A Bayesian approach for many-body uncertainties in nuclear structure: Many-body perturbation theory for finite nuclei, (2025), arXiv:2507.09079 [nucl-th]

  45. [45]

    Rutherford, M

    N. Rutherford, M. Mendes, I. Svensson, A. Schwenk, A. L. Watts, K. Hebeler, J. Keller, C. Prescod-Weinstein, D. Choud- hury, G. Raaijmakers,et al., Constraining the Dense Matter Equation of State with New NICER Mass–Radius Measure- ments and New Chiral Effective Field Theory Inputs, Astro- phys. J. Lett.971, L19 (2024)