Properties of the neutron star crust informed by nuclear structure data
Pith reviewed 2026-05-22 01:30 UTC · model grok-4.3
The pith
Bayesian neutron star models from nuclear data predict thicker crusts and higher crustal moments of inertia.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Starting from the full posterior distribution of nuclear-matter parameters obtained from nuclear structure data and employing the extended Thomas-Fermi method for the inner crust produces unified npeμ equations of state that satisfy mass, tidal-deformability, and radius constraints while yielding an increase in both the neutron-star surface thickness and the crustal moment of inertia relative to previous calculations.
What carries the argument
The extended Thomas-Fermi treatment of the inner crust inside a Bayesian sampling of Skyrme functionals whose parameters are drawn from the nuclear-structure posterior, thereby preserving correlations among bulk, surface, spin-orbit, and effective-mass terms.
If this is right
- The crustal moment of inertia rises, altering predictions for the angular-momentum reservoir available to power pulsar glitches.
- The neutron-star surface region becomes thicker than in earlier models.
- Only those equations of state that simultaneously respect the nuclear posterior and the astrophysical constraints from GW170817 and NICER are retained.
- Correlations among bulk, surface, and spin-orbit parameters survive the extrapolation to neutron-star densities.
Where Pith is reading between the lines
- Future precision measurements of nuclear surface properties could tighten the posterior and further shift the predicted crustal moment of inertia.
- The same consistent Bayesian pipeline could be applied to other neutron-star observables such as cooling curves or r-mode damping times.
- Systematic differences between glitch-inferred moments of inertia and the present predictions would test the limits of the Thomas-Fermi description at supranuclear densities.
Load-bearing premise
The posterior distribution of nuclear-matter parameters obtained at nuclear densities and asymmetries remains valid when extrapolated to the higher densities and greater isospin asymmetries of the neutron-star inner crust under the extended Thomas-Fermi approximation.
What would settle it
A measurement or glitch-based inference of the crustal moment of inertia that lies significantly below the range obtained from these models would indicate that the nuclear posterior does not hold under neutron-star crust conditions.
Figures
read the original abstract
We perform a Bayesian analysis of the neutron star (NS) equation of state (EoS) based on a wide set of Skyrme functionals, derived from previous nuclear physics inferences. The novelty of this approach lies in starting from the full multidimensional posterior distribution of nuclear matter parameters, consistent with a comprehensive set of static and dynamic nuclear structure observables. We construct unified EoSs for $npe\mu$ matter, where the inner crust of the NS is treated using an extended Thomas-Fermi method, providing for the first time a fully consistent Bayesian treatment of the correlation of bulk with surface as well as with spin-orbit and effective mass parameters. We then employ a standard Bayesian framework to identify those EoSs that satisfy astrophysical constraints from NS mass measurements, the tidal deformability from GW170817, and NICER mass-radius observations. We also examine NS observables, such as the crustal moment of inertia, which is crucial in understanding pulsar glitches. Compared to previous works, we observe an increase in both the NS surface thickness and the crustal moment of inertia.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript performs a Bayesian analysis of neutron-star equations of state by propagating the full multidimensional posterior of Skyrme parameters (constrained by static and dynamic nuclear observables) through an extended Thomas-Fermi treatment of the inner crust. Unified npeμ EoSs are constructed, filtered by NS mass, GW170817 tidal deformability, and NICER mass-radius data, and used to compute crust observables; the central result is an increase in both neutron-star surface thickness and crustal moment of inertia relative to earlier studies.
Significance. If the extrapolation of the nuclear posterior to the inner-crust regime remains valid, the work supplies a statistically consistent treatment of bulk-surface-spin-orbit correlations that could tighten predictions for the crustal moment of inertia relevant to glitch modeling. The use of the complete posterior rather than point estimates is a methodological strength that improves uncertainty propagation.
major comments (2)
- [§3] §3 (extended Thomas-Fermi crust construction): the headline increase in surface thickness and crustal MOI rests on applying the nuclear posterior—derived primarily from observables near saturation density and moderate asymmetry—to the low-density (∼0.01–0.5 ρ₀), highly neutron-rich inner-crust regime. No explicit sensitivity test or additional low-density constraint is shown to confirm that the ETF surface and gradient terms remain reliable under this extrapolation; this is load-bearing for the central claim.
- [Results section] Results section (comparison to previous works): the reported enhancement in crustal moment of inertia is presented without a quantitative decomposition of contributions from the Bayesian propagation versus the specific ETF functional form or the astrophysical filters; without this breakdown it is difficult to judge whether the increase is robust or an artifact of the chosen approximation.
minor comments (2)
- [Abstract] Abstract: the phrase 'for the first time a fully consistent Bayesian treatment' should be qualified by a brief statement of which prior Bayesian crust studies are being surpassed in scope.
- [Figures] Figure captions: ensure that the baseline models from earlier literature are explicitly identified so that the claimed increases can be directly traced.
Simulated Author's Rebuttal
We thank the referee for their thorough review and constructive feedback on our manuscript. We appreciate the recognition of the methodological strengths in using the full posterior distribution. Below we address each major comment in detail, indicating the revisions we will make to the manuscript.
read point-by-point responses
-
Referee: [§3] §3 (extended Thomas-Fermi crust construction): the headline increase in surface thickness and crustal MOI rests on applying the nuclear posterior—derived primarily from observables near saturation density and moderate asymmetry—to the low-density (∼0.01–0.5 ρ₀), highly neutron-rich inner-crust regime. No explicit sensitivity test or additional low-density constraint is shown to confirm that the ETF surface and gradient terms remain reliable under this extrapolation; this is load-bearing for the central claim.
Authors: We agree that the reliability of the extrapolation is important for the robustness of our results. The Skyrme parameters in our posterior are constrained by a wide range of nuclear observables that include information on surface properties and isovector terms, which influence the low-density behavior. The extended Thomas-Fermi approach is widely used for crust modeling precisely because it captures the relevant physics in this regime. Nevertheless, to strengthen the manuscript, we will add an explicit sensitivity test in the revised version by sampling subsets of the posterior with varied surface and gradient coefficients and reporting the resulting variations in surface thickness and crustal MOI. This will help quantify the uncertainty due to the extrapolation. revision: yes
-
Referee: [Results section] Results section (comparison to previous works): the reported enhancement in crustal moment of inertia is presented without a quantitative decomposition of contributions from the Bayesian propagation versus the specific ETF functional form or the astrophysical filters; without this breakdown it is difficult to judge whether the increase is robust or an artifact of the chosen approximation.
Authors: We acknowledge that a quantitative decomposition would aid in interpreting the origin of the observed increase in crustal moment of inertia. In the current work, the enhancement arises from the consistent treatment of correlations across the full posterior, but we did not isolate the individual contributions. In the revised manuscript, we will include additional comparisons: (i) results using the mean parameters versus the full posterior to highlight the effect of Bayesian propagation, and (ii) a brief discussion of how the ETF functional compares to other crust models in the literature. A complete isolation of all factors would require extensive additional calculations, but these additions will provide a clearer picture of the robustness. revision: partial
Circularity Check
Derivation chain is self-contained with independent inputs and outputs
full rationale
The paper takes a pre-existing multidimensional posterior on Skyrme parameters (constrained by static and dynamic nuclear observables) as input, constructs unified EoSs via extended Thomas-Fermi for the crust, and applies separate astrophysical filters from NS mass, GW170817, and NICER data. The reported increases in surface thickness and crustal moment of inertia are computed outputs of this pipeline. No equation or step reduces a claimed result to a fitted parameter by construction, and no load-bearing uniqueness theorem or ansatz is imported via self-citation. The derivation remains externally benchmarked against nuclear data and astrophysical constraints.
Axiom & Free-Parameter Ledger
free parameters (1)
- Skyrme functional parameters
axioms (2)
- domain assumption Skyrme energy density functionals provide a sufficient description of nuclear matter at the densities and asymmetries relevant to the neutron-star crust
- domain assumption The extended Thomas-Fermi method accurately captures surface and spin-orbit effects in the inner crust
Reference graph
Works this paper leans on
-
[1]
Properties of the neutron star crust informed by nuclear structure data
provides a relatively simple analytic representation of the nuclear EoS at zero temperature and for a given proton fraction, under the core assumption that matter is purely nucleonic. For astrophysical applications, it is complemented with electrons and muons in weak equilib- rium, forming cold-catalyzednpeµmatter. In this study, arXiv:2505.16929v2 [nucl-...
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[2]
P. Klausner, G. Col` o, X. Roca-Maza, and E. Vigezzi, Phys. Rev. C111, 014311 (2025), arXiv:2410.18598 [nucl-th]
- [3]
- [4]
- [5]
-
[6]
Adhikariet al.(PREX Collaboration), Phys
D. Adhikariet al.(PREX Collaboration), Phys. Rev. Lett.126, 172502 (2021)
work page 2021
-
[7]
Adhikariet al.(CREX Collaboration), Phys
D. Adhikariet al.(CREX Collaboration), Phys. Rev. Lett.129, 042501 (2022)
work page 2022
-
[8]
P.-G. Reinhard, X. Roca-Maza, and W. Nazarewicz, Phys. Rev. Lett.129, 232501 (2022)
work page 2022
-
[9]
P. Koliogiannis, E. Y¨ uksel, T. Ghosh, and N. Paar, The Astrophysical Journal Letters996, L18 (2025)
work page 2025
-
[10]
B. T. Reed, F. J. Fattoyev, C. J. Horowitz, and J. Piekarewicz, Phys. Rev. Lett.126, 172503 (2021)
work page 2021
-
[11]
F. J. Fattoyev, J. Piekarewicz, and C. J. Horowitz, Phys. Rev. Lett.120, 172702 (2018)
work page 2018
-
[12]
X. Roca-Maza, X. Vi˜ nas, M. Centelles, B. K. Agrawal, G. Col` o, N. Paar, J. Piekarewicz, and D. Vretenar, Phys. Rev. C92, 064304 (2015)
work page 2015
-
[13]
In- sights into the physics of neutron star interiors from pulsar glitches,
M. Antonelli, A. Montoli, and P. Pizzochero, “In- sights into the physics of neutron star interiors from pulsar glitches,” (World Scientific, 2022) pp. 219–281, arXiv:2301.12769 [astro-ph.HE]
-
[14]
The equation of state for dense nucleonic matter from a metamodeling. I. Foundational aspects
J. Margueron, R. Hoffmann Casali, and F. Gulminelli, Phys. Rev. C97, 025805 (2018), arXiv:1708.06894 [nucl- th]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[15]
Higher-order effects on the incompressibility of isospin asymmetric nuclear matter
L.-W. Chen, B.-J. Cai, C. M. Ko, B.-A. Li, C. Shen, and J. Xu, Phys. Rev. C80, 014322 (2009), arXiv:0905.4323 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2009
-
[16]
E. Chabanat, P. Bonche, P. Haensel, J. Meyer, and R. Schaeffer, Nuclear Physics A627, 710 (1997)
work page 1997
- [17]
-
[18]
N. N. Shchechilin, N. Chamel, J. M. Pearson, A. I. Chugunov, and A. Y. Potekhin, Phys. Rev. C109, 055802 (2024)
work page 2024
- [19]
- [20]
-
[21]
Skyrme Interaction and Nuclear Matter Constraints
M. Dutra, O. Louren¸ co, J. S. S´ a Martins, A. Delfino, J. R. Stone, and P. D. Stevenson, Phys. Rev. C85, 035201 (2012), arXiv:1202.3902 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2012
- [22]
-
[23]
G. Montefusco, M. Antonelli, and F. Gulminelli, Astronomy and Astrophysical694, A150 (2025), arXiv:2410.08008 [nucl-th]
-
[24]
A Massive Pulsar in a Compact Relativistic Binary
J. Antoniadis, P. C. C. Freire, N. Wex, T. M. Tauris, R. S. Lynch, M. H. van Kerkwijk, M. Kramer, C. Bassa, V. S. Dhillon, T. Driebe, J. W. T. Hessels, V. M. Kaspi, V. I. Kondratiev, N. Langer, T. R. Marsh, M. A. McLaugh- lin, T. T. Pennucci, S. M. Ransom, I. H. Stairs, J. van Leeuwen, J. P. W. Verbiest, and D. G. Whelan, Science 340, 448 (2013), arXiv:13...
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[25]
B. P. Abbott, others (LIGO Scientific Collaboration, and V. Collaboration), Physical Review X9, 011001 (2019), arXiv:1805.11579 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[26]
S. Vinciguerra, T. Salmi, A. L. Watts, D. Choudhury, T. E. Riley, P. S. Ray, S. Bogdanov, Y. Kini, S. Guillot, D. Chakrabarty, W. C. G. Ho, D. Huppenkothen, S. M. Morsink, Z. Wadiasingh, and M. T. Wolff, Astrophysical Journal961, 62 (2024), arXiv:2308.09469 [astro-ph.HE]
-
[27]
M. C. Miller, F. K. Lamb, A. J. Dittmann, S. Bogdanov, Z. Arzoumanian, K. C. Gendreau, S. Guillot, W. C. G. Ho, J. M. Lattimer, M. Loewenstein, S. M. Morsink, P. S. Ray, M. T. Wolff, C. L. Baker, T. Cazeau, S. Man- thripragada, C. B. Markwardt, T. Okajima, S. Pollard, I. Cognard, H. T. Cromartie, E. Fonseca, L. Guillemot, M. Kerr, A. Parthasarathy, T. T. ...
work page internal anchor Pith review Pith/arXiv arXiv 2021
-
[28]
D. Choudhury, T. Salmi, S. Vinciguerra, T. E. Riley, Y. Kini, A. L. Watts, B. Dorsman, S. Bogdanov, S. Guil- lot, P. S. Ray, D. J. Reardon, R. A. Remillard, A. V. Bilous, D. Huppenkothen, J. M. Lattimer, N. Ruther- ford, Z. Arzoumanian, K. C. Gendreau, S. M. Morsink, and W. C. G. Ho, The Astrophysical Journal Letters971, L20 (2024), arXiv:2407.06789 [astro-ph.HE]
-
[29]
L. Mauviard, S. Guillot, T. Salmi, D. Choudhury, B. Dorsman, D. Gonz´ alez-Caniulef, M. Hoogkamer, D. Huppenkothen, C. Kazantsev, Y. Kini, J.-F. Olive, P. Stammler, A. L. Watts, M. Mendes, N. Rutherford, A. Schwenk, I. Svensson, S. Bogdanov, M. Kerr, P. S. Ray, L. Guillemot, I. Cognard, and G. Theureau, The Astrophysical Journal995, 60 (2025)
work page 2025
-
[30]
Physics of Neutron Star Crusts
N. Chamel and P. Haensel, Living Reviews in Relativity 11, 10 (2008), arXiv:0812.3955 [astro-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2008
-
[32]
J. W. Negele and D. Vautherin, Nuclear Physics A207, 298 (1973)
work page 1973
-
[33]
S. K. Greif, G. Raaijmakers, K. Hebeler, A. Schwenk, and A. L. Watts, Mon. Not. Roy. Astron. Soc.485, 5363 16 (2019), arXiv:1812.08188 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv 2019
- [34]
- [35]
- [36]
-
[37]
G. Baym, C. Pethick, and P. Sutherland, Astrophys. J. 170, 299 (1971)
work page 1971
-
[38]
A unified equation of state of dense matter and neutron star structure
F. Douchin and P. Haensel, Astron. Astrophys.380, 151 (2001), arXiv:astro-ph/0111092 [astro-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2001
-
[39]
Neutron star radii and crusts: uncertainties and unified equations of state
M. Fortin, C. Providˆ encia, A. R. Raduta, F. Gulminelli, J. L. Zdunik, P. Haensel, and M. Bejger, Phys. Rev. C 94, 035804 (2016), arXiv:1604.01944 [astro-ph.SR]
work page internal anchor Pith review Pith/arXiv arXiv 2016
- [40]
-
[41]
T. Carreau, F. Gulminelli, and J. Margueron, European Physical Journal A55, 188 (2019), arXiv:1902.07032 [nucl-th]
- [42]
-
[43]
H. Dinh Thi, C. Mondal, and F. Gulminelli, Universe7, 373 (2021), arXiv:2109.09675 [astro-ph.HE]
- [44]
- [45]
-
[46]
M. Pi, X. Vinas, M. Barranco, A. Polls, and A. Perez- Canyellas, Astronomy and Astrophysics supplement se- ries64, 439 (1986)
work page 1986
-
[47]
M. Onsi, H. Przysiezniak, and J. M. Pearson, Phys. Rev. C55, 3139 (1997)
work page 1997
- [48]
-
[49]
F. m. c. Aymard, F. Gulminelli, and J. Margueron, Phys. Rev. C89, 065807 (2014)
work page 2014
-
[50]
H. Dinh Thi, T. Carreau, A. F. Fantina, and F. Gul- minelli, Astronomy and Astrophysics654, A114 (2021), arXiv:2109.13638 [nucl-th]
-
[51]
X. Vi˜ nas, C. Gonzalez-Boquera, B. K. Sharma, and M. Centelles, Acta Phys. Pol. B Proc. Suppl.10, 259 (2017)
work page 2017
-
[52]
Core-crust transition in neutron stars: predictivity of density developments
C. Ducoin, J. Margueron, C. Providˆ encia, and I. Vida˜ na, Phys. Rev. C83, 045810 (2011), arXiv:1102.1283 [nucl- th]
work page internal anchor Pith review Pith/arXiv arXiv 2011
- [53]
- [54]
- [55]
-
[56]
Isospin-dependent clusterization of Neutron-Star Matter
C. Ducoin, P. Chomaz, and F. Gulminelli, Nuclear Physics A789, 403 (2007), arXiv:nucl-th/0612044 [nucl- th]
work page internal anchor Pith review Pith/arXiv arXiv 2007
-
[57]
P. Haensel, A. Y. Potekhin, and D. G. Yakovlev,Neu- tron stars 1: Equation of state and structure, Vol. 326 (Springer, New York, USA, 2007)
work page 2007
-
[58]
A. Pastore, M. Shelley, S. Baroni, and C. A. Diget, Jour- nal of Physics G: Nuclear and Particle Physics44, 094003 (2017)
work page 2017
-
[59]
D. G. Ravenhall, C. J. Pethick, and J. M. Lattimer, Nuc. Phys. A407, 571 (1983)
work page 1983
-
[60]
J. M. Lattimer and D. F. Swesty, Nuc. Phys. A535, 331 (1991)
work page 1991
-
[61]
S. Anti´ c, D. Chatterjee, T. Carreau, and F. Gulminelli, Journal of Physics G Nuclear Physics46, 065109 (2019), arXiv:1901.03959 [nucl-th]
- [62]
-
[63]
A. Y. Potekhin, A. F. Fantina, N. Chamel, J. M. Pearson, and S. Goriely, Astronomy and Astrophysics560, A48 (2013), arXiv:1310.0049 [astro-ph.SR]
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[64]
T. Carreau, F. Gulminelli, N. Chamel, A. F. Fantina, and J. M. Pearson, Astronomy and Astrophysics635, A84 (2020), arXiv:1912.01265 [astro-ph.HE]
-
[65]
Forecasting neutron star temperatures: predictability and variability
D. Page and S. Reddy, Phys. Rev. Lett.111, 241102 (2013), arXiv:1307.4455 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[66]
B. Link, R. I. Epstein, and J. M. Lattimer, Phys. Rev. Lett.83, 3362 (1999), arXiv:astro-ph/9909146 [astro-ph]
work page internal anchor Pith review Pith/arXiv arXiv 1999
-
[67]
N. Andersson, K. Glampedakis, W. C. G. Ho, and C. M. Espinoza, Phys. Rev. Lett.109, 241103 (2012)
work page 2012
- [68]
-
[69]
A. Montoli, M. Antonelli, B. Haskell, and P. Pizzochero, Universe7(2021), 10.3390/universe7010008
-
[70]
Superfluidity and Superconductivity in Neutron Stars
N. Chamel, Journal of Astrophysics and Astronomy38, 43 (2017), arXiv:1709.07288 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[71]
L. Gavassino, M. Antonelli, and B. Haskell, Universe7, 28 (2021), arXiv:2012.10288 [astro-ph.HE]
- [72]
- [73]
-
[74]
K. Sekizawa, S. Kobayashi, and M. Matsuo, Phys. Rev. C105, 045807 (2022), arXiv:2112.14350 [nucl-th]
-
[75]
Superfluid fraction in the crystalline crust of a neutron star: role of BCS pairing
N. Chamel, Phys. Rev. C111, 045803 (2025), arXiv:2412.05599 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv 2025
- [76]
-
[77]
P. Papakonstantinou, J. Margueron, F. Gulminelli, and A. R. Raduta, Phys. Rev. C88, 045805 (2013)
work page 2013
-
[78]
Giant Pulsar Glitches and the Inertia of Neutron-Star Crusts
T. Delsate, N. Chamel, N. G¨ urlebeck, A. F. Fantina, J. M. Pearson, and C. Ducoin, Phys. Rev. D94, 023008 (2016), arXiv:1606.00016 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[79]
B. P. Abbott, others (LIGO Scientific Collaboration, and V. Collaboration), Phys. Rev. Lett.121, 161101 (2018), arXiv:1805.11581 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[80]
A. Sorensen, ,et al., Progress in Particle and Nuclear Physics134, 104080 (2024), arXiv:2301.13253 [nucl-th]
-
[81]
E. Giliberti, G. Cambiotti, M. Antonelli, and P. M. Piz- zochero, Monthly Notices of the Royal Astronomical So- ciety491, 1064 (2020), arXiv:1809.08542 [astro-ph.HE]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.