pith. sign in

arxiv: 2606.02932 · v1 · pith:FQAVHQONnew · submitted 2026-06-01 · 🌌 astro-ph.HE

On a re-examination of neutron star cooling in transient sources -- No shallow heating required?

Pith reviewed 2026-06-28 12:46 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords neutron star coolinglow-mass X-ray binariesshallow heatingthermonuclear burningcrust envelope boundaryaccretion outburstseffective temperaturetransient sources
0
0 comments X

The pith

Neutron star cooling after seven of eight outbursts matches predictions from envelope heat leakage alone, without shallow heating.

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

The paper tests whether cooling of neutron stars in low-mass X-ray binaries after accretion outbursts requires an unknown extra heat source in the crust. It replaces that source with a revised boundary condition at the envelope-crust interface that lets a fraction of thermonuclear heat from burning accreted hydrogen and helium flow downward, with the fraction set by the accretion rate. Simulations using this condition reproduce the observed effective temperatures for seven out of eight episodes from seven different sources. A sympathetic reader cares because the result removes the need for an ad hoc heating term whose physical origin has remained unclear and ties the crust temperature directly to observable accretion properties. For one source the match remains only qualitative and extra physics is still needed.

Core claim

By equipping the nscool code with a new boundary condition at the envelope/crust transition that incorporates thermonuclear heating leakage from the envelope into the crust depending on mass accretion rate, the cooling curves of neutron stars in seven out of eight outburst episodes from seven sources can be explained without invoking ad hoc shallow heating. The model evolves the fully relativistic crust and core. For the remaining source, qualitative features of the cooling curve are captured but a good quantitative fit still requires additional physics.

What carries the argument

A new boundary condition at the envelope/crust transition that incorporates thermonuclear heating leakage from the envelope into the crust as a function of the mass accretion rate.

If this is right

  • Observed post-outburst effective temperatures are reproduced using only thermonuclear burning and standard heat transport.
  • Shallow heating is not required to explain cooling in seven of the eight cases examined.
  • The amount of heat entering the crust is set by the accretion-rate-dependent leakage rather than an independent proportionality constant.
  • For the source EXO, the boundary condition alone is insufficient for a quantitative fit and additional physics must be included.

Where Pith is reading between the lines

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

  • The same boundary condition could be applied to other accretion-powered neutron-star systems to test whether their long-term temperatures also follow from envelope leakage.
  • If the accretion-rate dependence holds, it would link variations in outburst history directly to differences in subsequent cooling without separate heating mechanisms.
  • High-precision temperature measurements at late times in sources with extreme accretion rates could directly test the functional form of the leakage.

Load-bearing premise

The new boundary condition at the envelope/crust transition correctly captures the leakage of thermonuclear heating from the envelope into the crust as a function of mass accretion rate.

What would settle it

Cooling data from one of the modeled sources, or from a new outburst with a markedly different accretion rate, that deviates from the predicted temperature evolution while the accretion history is known.

Figures

Figures reproduced from arXiv: 2606.02932 by Dany Page, Martin Nava-Callejas, Yuri Cavecchi.

Figure 1
Figure 1. Figure 1: Best fitting cooling curves for Scenario A. [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
read the original abstract

Context: For the typical modeling of neutron stars cooling after an accretion episode in Low-Mass X-ray Binaries, an extra heating source of unknown physical origin, \textit{the shallow heating}, is invoked in order to account for the inferred high effective temperatures of the star up to hundreds of days after the end of accretion. The amount of the shallow heating generated in the crust is usually taken to be proportional to the accretion rate, although the proportionality constant may change from source to source. Aims: In this paper, we intend to model the effective temperature data of eight outburst episodes from seven different sources (\mxb, \xte, \exo, \igr, \swift, \rxs\ and \ks) without {\it ad hoc} shallow heating but accounting for the presence of thermonuclear heating due to the burning of the accreted H/He. Methods: We employ the fully relativistic code \texttt{nscool}, which simulates both the crust and core of a neutron star, equipped with a new boundary condition at the envelope/crust transition which considers thermonuclear heating leakage from the envelope into the crust and depends on the mass accretion rate. Results: We find that the neutron star cooling for seven out of eight of these outbursts can be well explained with this new boundary condition and without the requirement of {\it ad hoc} shallow heating. While the qualitative features of \exo\ cooling curve can be explained, a good fitting still requires additional physics.

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 / 1 minor

Summary. The manuscript re-examines post-outburst cooling of neutron stars in seven LMXBs using the nscool code. It introduces a new boundary condition at the envelope/crust transition that incorporates thermonuclear heating leakage as a function of mass accretion rate, claiming this eliminates the need for ad hoc shallow heating in seven of eight outbursts. For the eighth source (EXO), only qualitative features are reproduced and additional physics is required for a good fit.

Significance. If the boundary condition is shown to be derived from first-principles envelope calculations without source-specific tuning or adjustable constants, the result would be significant: it replaces an empirical shallow-heating term with a physically motivated, accretion-rate-dependent boundary condition, reducing the number of free parameters in cooling models. The use of the established, fully relativistic nscool code is a strength for reproducibility. The work directly addresses a long-standing modeling issue in transient neutron-star cooling.

major comments (2)
  1. [Methods] Methods: The new boundary condition is stated to depend on accretion rate and to capture thermonuclear heating leakage, but no explicit functional form, derivation from envelope structure equations, or demonstration that it contains no adjustable constants fitted to the cooling data is provided. This is load-bearing for the central claim that shallow heating is unnecessary, as the result would otherwise relocate rather than eliminate ad hoc elements.
  2. [Results] Results: The assertion that seven of eight outbursts 'can be well explained' lacks any quantitative fit statistics (e.g., χ^{2} values, residuals, or uncertainties on model parameters) or details on data selection and error treatment. Without these, the goodness-of-fit cannot be assessed, particularly since the eighth case (EXO) still requires extra physics.
minor comments (1)
  1. Abstract and throughout: Source names appear as LaTeX placeholders (\mxb, \xte, etc.); replace with actual designations (e.g., MXB 1659-29) for clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive comments. We address each major comment below and will revise the manuscript to incorporate additional details and quantitative information as requested.

read point-by-point responses
  1. Referee: [Methods] Methods: The new boundary condition is stated to depend on accretion rate and to capture thermonuclear heating leakage, but no explicit functional form, derivation from envelope structure equations, or demonstration that it contains no adjustable constants fitted to the cooling data is provided. This is load-bearing for the central claim that shallow heating is unnecessary, as the result would otherwise relocate rather than eliminate ad hoc elements.

    Authors: We agree that the explicit functional form, derivation, and confirmation of no fitted constants are necessary to support the claim. The boundary condition is constructed from envelope calculations that incorporate thermonuclear heating leakage depending on accretion rate. In the revised version we will add a dedicated methods subsection (or appendix) providing the explicit functional form, its derivation from the envelope structure equations, and a demonstration that it introduces no adjustable constants tuned to the cooling curves themselves. This will clarify that the approach is physically motivated rather than relocating ad hoc elements. revision: yes

  2. Referee: [Results] Results: The assertion that seven of eight outbursts 'can be well explained' lacks any quantitative fit statistics (e.g., χ^{2} values, residuals, or uncertainties on model parameters) or details on data selection and error treatment. Without these, the goodness-of-fit cannot be assessed, particularly since the eighth case (EXO) still requires extra physics.

    Authors: We acknowledge that quantitative fit statistics and data-treatment details were omitted from the submitted manuscript. The revised version will include χ² values (or equivalent goodness-of-fit metrics), residual plots, parameter uncertainties, and explicit descriptions of data selection and error treatment for the seven outbursts that are reproduced. For the EXO case we will expand the discussion of the additional physics required and why a quantitative fit is not achieved with the present boundary condition. revision: yes

Circularity Check

0 steps flagged

New boundary condition replaces shallow heating with accretion-rate dependence; no reduction to fitted quantity by construction

full rationale

The paper introduces a new boundary condition at the envelope/crust transition that depends on mass accretion rate to encode thermonuclear heating leakage, then uses it in nscool to model cooling curves for eight outbursts. Seven are explained without additional shallow heating. No quoted equations or text show the BC functional form or constants being fitted to the target cooling data itself, nor any self-citation chain that bears the central claim, nor renaming of a known result. The derivation remains self-contained against the external cooling observations once the BC is accepted as an input motivated by envelope physics. This matches the normal non-circular case (score 0-2).

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the accuracy of the nscool code and the new boundary condition for thermonuclear heating leakage; no new free parameters or invented entities are introduced beyond standard neutron star cooling assumptions.

axioms (2)
  • domain assumption The nscool code provides an accurate fully relativistic simulation of neutron star crust and core cooling with standard microphysical inputs.
    The paper employs this code as the simulation engine for all results.
  • domain assumption Thermonuclear heating from accreted H/He burning leaks into the crust according to a boundary condition that depends on the mass accretion rate.
    This is the key modeling choice replacing shallow heating, stated in the methods.

pith-pipeline@v0.9.1-grok · 5804 in / 1398 out tokens · 38191 ms · 2026-06-28T12:46:22.106687+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

122 extracted references · 56 canonical work pages · 20 internal anchors

  1. [1]

    P., Abbott, R., Abbott, T

    Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017, Phys. Rev. Lett., 119, 161101

  2. [2]

    R., & Ravenhall , D

    Akmal , A., Pandharipande , V. R., & Ravenhall , D. G. 1998, , 58, 1804

  3. [3]

    O., Sivakoff, G

    Bahramian, A., Heinke, C. O., Sivakoff, G. R., et al. 2013, The Astrophysical Journal, 780, 127

  4. [4]

    Brown , E. F. & Cumming , A. 2009, , 698, 1020

  5. [5]

    F., Cumming , A., Fattoyev , F

    Brown , E. F., Cumming , A., Fattoyev , F. J., et al. 2018, , 120, 182701

  6. [6]

    K., Heger , A., et al

    Cavecchi , Y., Galloway , D. K., Heger , A., et al. 2026, , 998, 25

  7. [7]

    1993, Nuclear Physics A, 555, 59

    Chen, J., Clark, J., Dav \'e , R., & Khodel, V. 1993, Nuclear Physics A, 555, 59

  8. [8]

    2016, Monthly Notices of the Royal Astronomical Society, 461, 4049

    Degenaar, N., Altamirano, D., Parker, M., et al. 2016, Monthly Notices of the Royal Astronomical Society, 461, 4049

  9. [9]

    2021, , 508, 882

    Degenaar , N., Page , D., van den Eijnden , J., et al. 2021, , 508, 882

  10. [10]

    & Wijnands, R

    Degenaar, N. & Wijnands, R. 2011, Monthly Notices of the Royal Astronomical Society: Letters, 414, L50

  11. [11]

    2015, , 451, 2071

    Degenaar , N., Wijnands , R., Bahramian , A., et al. 2015, , 451, 2071

  12. [12]

    F., et al

    Degenaar, N., Wijnands, R., Brown, E. F., et al. 2013, The Astrophysical Journal, 775, 48

  13. [13]

    Feigelson , E. D. & Babu , G. J. 2012, Modern Statistical Methods for Astronomy

  14. [14]

    K., Homan , J., Wijnands , R., et al

    Fridriksson , J. K., Homan , J., Wijnands , R., et al. 2011, , 736, 162

  15. [15]

    Y., Hanawa , T., Iben , I., J., & Richardson , M

    Fujimoto , M. Y., Hanawa , T., Iben , I., J., & Richardson , M. B. 1984, , 278, 813

  16. [16]

    Y., Hanawa , T., Iben , Icko, J., & Richardson , M

    Fujimoto , M. Y., Hanawa , T., Iben , Icko, J., & Richardson , M. B. 1987, , 315, 198

  17. [17]

    K., Muno , M

    Galloway , D. K., Muno , M. P., Hartman , J. M., Psaltis , D., & Chakrabarty , D. 2008, , 179, 360

  18. [18]

    H., Pethick , C

    Gudmundsson , E. H., Pethick , C. J., & Epstein , R. I. 1982, , 259, L19

  19. [19]

    F., Schatz , H., M \"o ller , P., & Kratz , K.-L

    Gupta , S., Brown , E. F., Schatz , H., M \"o ller , P., & Kratz , K.-L. 2007, , 662, 1188

  20. [20]

    Gusakov , M. E. & Chugunov , A. I. 2020, , 124, 191101

  21. [21]

    & Zdunik , J

    Haensel , P. & Zdunik , J. L. 2008, , 480, 459

  22. [22]

    & Fujimoto , M

    Hanawa , T. & Fujimoto , M. Y. 1984, , 36, 199

  23. [23]

    F., Schatz, H., et al

    Jain, R., Brown, E. F., Schatz, H., et al. 2025, Phys. Rev. Lett., 135, 162701

  24. [24]

    H., Ingram , A., Middleton , M., & Drake , J

    Knight , A. H., Ingram , A., Middleton , M., & Drake , J. 2022, , 510, 4736

  25. [25]

    R., et al

    Lin , D., Strader , J., Carrasco , E. R., et al. 2018, Nature Astronomy, 2, 656

  26. [26]

    2012, The Astrophysical Journal, 748, 82

    Linares, M., Altamirano, D., Chakrabarty, D., Cumming, A., & Keek, L. 2012, The Astrophysical Journal, 748, 82

  27. [27]

    L., Cackett , E

    Merritt , R. L., Cackett , E. M., Brown , E. F., et al. 2016, , 833, 186

  28. [28]

    1990, , 362, 572

    Miralda-Escude , J., Haensel , P., & Paczynski , B. 1990, , 362, 572

  29. [29]

    2024, RAS Techniques and Instruments, 3, 800

    Nava-Callejas , M., Cavecchi , Y., & Page , D. 2024, RAS Techniques and Instruments, 3, 800

  30. [30]

    2025, The Astrophysical Journal, 991, 62

    Nava-Callejas, M., Page, D., & Cavecchi, Y. 2025, The Astrophysical Journal, 991, 62

  31. [31]

    S., Page , D., Wijnands , R., & Degenaar , N

    Ootes , L. S., Page , D., Wijnands , R., & Degenaar , N. 2016, , 461, 4400

  32. [32]

    S., Vats , S., Page , D., et al

    Ootes , L. S., Vats , S., Page , D., et al. 2019 a , , 487, 1447

  33. [33]

    S., Wijnands , R., & Page , D

    Ootes , L. S., Wijnands , R., & Page , D. 2019 b , , 630, A95

  34. [34]

    2006, , 441, 1115

    \"O zel , F. 2006, , 441, 1115

  35. [35]

    O zel , F., Gould , A., & G \

    \"O zel , F., Gould , A., & G \"u ver , T. 2012, , 748, 5

  36. [36]

    2016, NSCool: Neutron star cooling code , Astrophysics Source Code Library, record ascl:1609.009

    Page , D. 2016, NSCool: Neutron star cooling code , Astrophysics Source Code Library, record ascl:1609.009

  37. [37]

    2022, The Astrophysical Journal, 933, 216

    Page, D., Homan, J., Nava-Callejas, M., et al. 2022, The Astrophysical Journal, 933, 216

  38. [38]

    & Reddy , S

    Page , D. & Reddy , S. 2013, , 111, 241102

  39. [39]

    S., Wijnands , R., Degenaar , N., Ootes , L., & Page , D

    Parikh , A. S., Wijnands , R., Degenaar , N., Ootes , L., & Page , D. 2018, , 476, 2230

  40. [40]

    S., Wijnands, R., Degenaar, N., et al

    Parikh, A. S., Wijnands, R., Degenaar, N., et al. 2017, Monthly Notices of the Royal Astronomical Society, 466, 4074

  41. [41]

    S., Wijnands , R., Homan , J., et al

    Parikh , A. S., Wijnands , R., Homan , J., et al. 2020, , 638, L2

  42. [42]

    S., Wijnands , R., Ootes , L

    Parikh , A. S., Wijnands , R., Ootes , L. S., et al. 2019, , 624, A84

  43. [43]

    Potekhin , A. Y. & Chabrier , G. 2021, , 645, A102

  44. [44]

    Y., Chabrier, G., & Yakovlev, D

    Potekhin, A. Y., Chabrier, G., & Yakovlev, D. G. 1997, Astron. Astrophys., 323, 415

  45. [45]

    Y., Chugunov , A

    Potekhin , A. Y., Chugunov , A. I., Shchechilin , N. N., & Gusakov , M. E. 2025, Journal of High Energy Astrophysics, 45, 116

  46. [46]

    & Reddy , S

    Roggero , A. & Reddy , S. 2016, , 94, 015803

  47. [47]

    Schwenk, A., Friman, B., & Brown, G. E. 2003, Nuclear Physics A, 713, 191

  48. [48]

    N., & Pons , J

    Turlione , A., Aguilera , D. N., & Pons , J. A. 2015, , 577, A5

  49. [49]

    2017, Journal of Astrophysics and Astronomy, 38, 49

    Wijnands , R., Degenaar , N., & Page , D. 2017, Journal of Astrophysics and Astronomy, 38, 49

  50. [50]

    P., Miller , J

    Wijnands , R., Muno , M. P., Miller , J. M., et al. 2002, , 566, 1060

  51. [51]

    L., Haensel , P., Paczynski , B., & Miralda-Escude , J

    Zdunik , J. L., Haensel , P., Paczynski , B., & Miralda-Escude , J. 1992, , 384, 129

  52. [52]

    , keywords =

    A Solution to the Tension of Burning on Neutron Stars and Nuclear Physics. , keywords =. doi:10.3847/1538-4357/ae2ff3 , archivePrefix =. 2506.11966 , primaryClass =

  53. [53]

    Cooling of Accretion-Heated Neutron Stars

    Cooling of Accretion-Heated Neutron Stars. Journal of Astrophysics and Astronomy , keywords =. doi:10.1007/s12036-017-9466-5 , archivePrefix =. 1709.07034 , primaryClass =

  54. [54]

    and Chamel, N

    Allard, V. and Chamel, N. , date-added =. Gapless Neutron Superfluidity Can Explain the Late Time Cooling of Transiently Accreting Neutron Stars , url =. Phys. Rev. Lett. , month =. 2024 , bdsk-url-1 =. doi:10.1103/PhysRevLett.132.181001 , issue =

  55. [55]

    and Fattoyev, Farrukh J

    Cumming, Andrew and Brown, Edward F. and Fattoyev, Farrukh J. and Horowitz, C. J. and Page, Dany and Reddy, Sanjay , doi =. Lower limit on the heat capacity of the neutron star core , url =. Phys. Rev. C , month =. 2017 , bdsk-url-1 =

  56. [56]

    doi:10.1088/0004-637X/698/2/1020 , eprint =

    , keywords =. doi:10.1088/0004-637X/698/2/1020 , eprint =

  57. [57]
  58. [58]

    Quiescent thermal emission from neutron stars in LMXBs

    doi:10.1051/0004-6361/201322690 , eid =. arXiv , author =:1309.3909 , journal =

  59. [59]

    doi:10.1016/j.jheap.2024.11.017 , eprint =

    Journal of High Energy Astrophysics , keywords =. doi:10.1016/j.jheap.2024.11.017 , eprint =

  60. [60]

    doi:10.1086/161851 , journal =

  61. [61]

    doi:10.1086/165124 , journal =

  62. [62]

    doi:10.1086/169295 , journal =

  63. [63]

    doi:10.1086/170858 , journal =

  64. [64]

    arXiv , author =:2001.08367 , journal =

    doi:10.1093/ptep/ptaa010 , eid =. arXiv , author =:2001.08367 , journal =

  65. [65]

    arXiv , author =:2105.13958 , journal =

    doi:10.3847/1538-4357/ac2821 , eid =. arXiv , author =:2105.13958 , journal =

  66. [66]

    arXiv , author =:2208.14622 , journal =

    doi:10.3847/1538-4357/ac8dfe , eid =. arXiv , author =:2208.14622 , journal =

  67. [67]

    doi:10.1086/183840 , journal =

  68. [68]

    arXiv , author =:astro-ph/9706148 , journal =

  69. [69]

    doi:10.1086/517869 , eprint =

    , keywords =. doi:10.1086/517869 , eprint =

  70. [70]

    doi:10.1051/0004-6361:20078578 , eprint =

    , keywords =. doi:10.1051/0004-6361:20078578 , eprint =

  71. [71]

    arXiv , author =:2004.04195 , journal =

    doi:10.1103/PhysRevLett.124.191101 , eid =. arXiv , author =:2004.04195 , journal =

  72. [72]

    , keywords =

    Long-term temperature evolution of neutron stars undergoing episodic accretion outbursts. , keywords =. doi:10.1051/0004-6361/201936035 , archivePrefix =. 1906.02554 , primaryClass =

  73. [73]

    Statistical Methods for Astronomy

    Modern Statistical Methods for Astronomy. doi:10.48550/arXiv.1205.2064 , adsurl =

  74. [74]

    ascl , author =:1609.009 , howpublished =

  75. [75]

    A luminous X-ray outburst from an intermediate-mass black hole in an off-centre star cluster

    A luminous X-ray outburst from an intermediate-mass black hole in an off-centre star cluster. Nature Astronomy , keywords =. doi:10.1038/s41550-018-0493-1 , archivePrefix =. 1806.05692 , primaryClass =

  76. [76]

    Continued cooling of the accretion-heated neutron star crust in the X-ray transient IGR J17480-2446 located in the globular cluster Terzan 5

    Continued cooling of the accretion-heated neutron star crust in the X-ray transient IGR J17480-2446 located in the globular cluster Terzan 5. , keywords =. doi:10.1093/mnras/stz1406 , archivePrefix =. 1805.00610 , primaryClass =

  77. [77]

    , keywords =

    Constraining the properties of dense neutron star cores: the case of the low-mass X-ray binary HETE J1900.1-2455. , keywords =. doi:10.1093/mnras/stab2202 , archivePrefix =. 2112.09711 , primaryClass =

  78. [78]

    Beznogov and Nathalie Degenaar and Rudy Wijnands and Aastha S

    Dany Page and Jeroen Homan and Martin Nava-Callejas and Yuri Cavecchi and Mikhail V. Beznogov and Nathalie Degenaar and Rudy Wijnands and Aastha S. Parikh , doi =. A ``Hyperburst'' in the MAXI J0556--332 Neutron Star: Evidence for a New Type of Thermonuclear Explosion , url =. The Astrophysical Journal , month =. 2022 , bdsk-url-1 =

  79. [79]

    doi:10.1093/rasti/rzae055 , eprint =

    RAS Techniques and Instruments , keywords =. doi:10.1093/rasti/rzae055 , eprint =

  80. [80]

    Thermonuclear Heating of Accreting Neutron Stars , url =

    Nava-Callejas, Mart. Thermonuclear Heating of Accreting Neutron Stars , url =. 2025 , bdsk-url-1 =. doi:10.3847/1538-4357/adf8e4 , journal =

Showing first 80 references.