Pith. Sign in

REVIEW

Minimal Cooling of Neutron Stars: A New Paradigm

PaperPDF from arXiv

data files

ReviewThe objections and verdict, as JSON Paper breakdownThe authors' claims, conditions and measurements, as JSON Signed receiptProof this review has not been altered
X Bluesky LinkedIn Reddit Hacker News
Dany Page, James M. Lattimer, Madappa Prakash, Andrew W. Steiner

Astrophysics · Submitted 29 Mar 2004 · 139 references

Pith's verdict: The paper that named the null hypothesis — minimal cooling with PBF — and did the four-parameter survey carefully enough that the field still uses it as the baseline.

T0 · AI review

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on October 29, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 4 major objections · 6 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDastro-ph/0403657 (version 2)
Pith IDKZ52LIQH, Pith's permanent record of this paper and its review
CitationAndrew W. Steiner, Dany Page, James M. Lattimer, Madappa Prakash. "Minimal Cooling of Neutron Stars: A New Paradigm." arXiv:astro-ph/0403657 (2004). https://pith.science/paper/KZ52LIQH
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

The paper argues that a cooling model with Cooper-pair neutrino emission and no direct Urca process can account for every clearly detected thermal neutron star, leaving only a few luminosity upper limits as candidates for enhanced cooling.

The reading

The paper proposes a new baseline for how isolated neutron stars cool. It replaces the old standard-cooling picture with a minimal-cooling scenario: neutrino losses come from modified Urca and nucleon bremsstrahlung processes plus neutrino emission from the breaking and formation of Cooper pairs in superfluid neutrons and protons, but never from any direct Urca process involving nucleons or exotic constituents. The authors treat the prediction as a four-parameter family set by the dense-matter equation of state, the pairing gaps, the envelope composition, and the stellar mass, and they compute cooling curves across wide ranges of all four. They then compare those curves with the inferred temperatures, luminosities, and ages of thermally emitting neutron stars. The claim is that no clearly detected star requires enhanced cooling, so the minimal scenario can serve as the null hypothesis against which any future claim of exotic matter or direct Urca emission must be judged.

What carries the argument

The mechanism that carries the argument is neutrino emission from Cooper-pair breaking and formation (PBF). When a neutron star cools through the critical temperature of a superfluid pairing gap, Cooper pairs continually break and re-form, and each event can emit a neutrino-antineutrino pair; the resulting luminosity has a $T^7$ dependence and can exceed the modified Urca and bremsstrahlung rates near $T_c$. The paper folds this into a four-parameter cooling calculation: the equation of state, the $^1S_0$ and $^3P_2$ pairing gaps, the light-element content of the envelope, and the neutron-star mass. The PBF channel is what lets the minimal scenario produce temperatures low enough to match the young pulsars without invoking a direct Urca process.

Extended reading notes

Core claim

The central result is that a cooling model with no direct Urca channel, but with Cooper-pair breaking and formation added, is at least marginally compatible with every neutron star whose thermal emission has been clearly detected. The authors argue that the pair-breaking process dominates modified Urca near the pairing critical temperature and can make the star cool faster than the old standard scenario, so the observed low temperatures of young pulsars such as Vela do not by themselves demand exotic cooling. Within this baseline, the equation of state and stellar mass matter little, while the neutron $^3P_2$ gap and the envelope composition carry most of the uncertainty. The paper identifies two upper limits, for PSR J0205+6449 in 3C58 and RX J0007.0+7302 in CTA 1, that fall below all of its minimal models and therefore point toward some enhanced neutrino emission.

Load-bearing premise

The load-bearing premise is the size and density dependence of the neutron $^3P_2$ pairing gap, which the authors themselves call the most uncertain input and handle by testing a vanishing gap and three hand-chosen cases.

Editorial extensions

If this is right

  • If the minimal scenario is right, no detected neutron star requires a direct Urca process, so the nuclear symmetry energy must keep the direct Urca threshold above the central density of the observed stars.
  • The equation of state and stellar mass become secondary parameters for cooling, while the neutron $^3P_2$ gap and envelope composition control most of the spread in predicted temperatures.
  • Vela and PSR 1706-44 can be explained without exotic matter only if the $^3P_2$ gap is large enough and, for Vela, if the star is low mass with a light-element envelope and pairing through most of the core.
  • PSR J0205+6449 in 3C58 and RX J0007.0+7302 in CTA 1 become the leading candidates for physics beyond minimal cooling, because their luminosity upper limits lie below every minimal model.
  • The older stars PSR 0656+14, PSR 1055-52, Geminga, and RX J0720.4-3125 would require an internal heating source only if their ages or luminosities sit near the upper end of their estimated ranges.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A population study of young supernova remnants could test the PBF prediction directly: the minimal scenario predicts a characteristic change in slope of the temperature-age relation at the pairing critical temperature, not a single power law.
  • If 3C58 and CTA 1 are confirmed to contain neutron stars, their low limits would favor a mass threshold for enhanced cooling, which could be checked by measuring neutron-star masses in those remnants.
  • Time-dependent envelope composition offers an alternative to internal heating for the old stars: a star whose light-element envelope is gradually replaced by heavy elements can mimic a heating episode in the luminosity-age plane.
  • The same four-parameter family makes a falsifiable split between cooling dominated by neutron pairing and cooling dominated by proton pairing; future measurements of surface temperatures in the $10^8$-$10^9$ K range could identify which gap is controlling the neutrino losses.
Original abstract

A new classification of neutron star cooling scenarios, involving either ``minimal'' cooling or ``enhanced'' cooling is proposed. The minimal cooling scenario replaces and extends the so-called standard cooling scenario to include neutrino emission from the Cooper pair breaking and formation process. This emission dominates that due to the modified Urca process for temperatures close to the critical temperature for superfluid pairing. Minimal cooling is distinguished from enhanced cooling by the absence of neutrino emission from any direct Urca process, due either to nucleons or to exotica. Within the minimal cooling scenario, theoretical cooling models can be considered to be a four parameter family involving the equation of state of dense matter, superfluid properties of dense matter, the composition of the neutron star envelope, and the mass of the neutron star. Consequences of minimal cooling are explored through extensive variations of these parameters. Results are compared with the inferred properties of thermally-emitting neutron stars in order to ascertain if enhanced cooling occurs in any of them. All stars for which thermal emissions have been clearly detected are at least marginally consistent with the lack of enhanced cooling. The two pulsars PSR 0833-45 (Vela) and PSR 1706-44 would require enhanced cooling in case their ages and/or temperatures are on the lower side of their estimated values whereas the four stars PSR 0656+14, PSR 1055-52, Geminga, and RX J0720.4-3125 may require some source of internal heating in case their age and/or luminosity are on the upper side of their estimated values. The new upper limits on the thermal luminosity of PSR J0205+6449 and RX J0007.0+7302 are indicative of the occurrence of some enhanced neutrino emission beyond the minimal scenario.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 7,846 input tokens.

What the abstract itself states, as probabilities: claims first 0.12 · states limitation 0.57 · human participants 0.01 · controlled comparison 0.38 · randomized assignment 0.01 · states quantitative result 0.37 · compares with prior results 0.78 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v2. Results and the main claim are underlined; .

A new classification of neutron star cooling scenarios, involving either ``minimal'' cooling or ``enhanced'' cooling is proposed. The minimal cooling scenario replaces and extends the so-called standard cooling scenario to include neutrino emission from the Cooper pair breaking and formation process. This emission dominates that due to the modified Urca process for temperatures close to the critical temperature for superfluid pairing. Minimal cooling is distinguished from enhanced cooling by the absence of neutrino emission from any direct Urca process, due either to nucleons or to exotica. Within the minimal cooling scenario, theoretical cooling models can be considered to be a four parameter family involving the equation of state of dense matter, superfluid properties of dense matter, the composition of the neutron star envelope, and the mass of the neutron star. Consequences of minimal cooling are explored through extensive variations of these parameters. Results are compared with the inferred properties of thermally-emitting neutron stars in order to ascertain if enhanced cooling occurs in any of them. All stars for which thermal emissions have been clearly detected are at least marginally consistent with the lack of enhanced cooling. The two pulsars PSR 0833-45 (Vela) and PSR 1706-44 would require enhanced cooling in case their ages and/or temperatures are on the lower side of their estimated values whereas the four stars PSR 0656+14, PSR 1055-52, Geminga, and RX J0720.4-3125 may require some source of internal heating in case their age and/or luminosity are on the upper side of their estimated values. The new upper limits on the thermal luminosity of PSR J0205+6449 and RX J0007.0+7302 are indicative of the occurrence of some enhanced neutrino emission beyond the minimal scenario.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 A new classification of neutron star cooling scenarios, involving either ``minimal'' cooling or ``enhanced'' cooling is proposed. contribution 1.00 0.85 0.13 asserted 0.98
s1 The minimal cooling scenario replaces and extends the so-called standard cooling scenario to include neutrino emission from the Cooper pair breaking and formation process. contribution 0.96 0.64 0.07 asserted 1.00
s2 This emission dominates that due to the modified Urca process for temperatures close to the critical temperature for superfluid pairing. background 0.39 0.25 0.03 asserted 1.00
s3 Minimal cooling is distinguished from enhanced cooling by the absence of neutrino emission from any direct Urca process, due either to nucleons or to exotica. contribution 0.77 0.24 0.07 asserted 1.00
s4 Within the minimal cooling scenario, theoretical cooling models can be considered to be a four parameter family involving the equation of state of dense matter, superfluid properties of dense matter, the composition of the neutron star envelope, and the mass of the neutron star. method 0.80 0.62 0.05 asserted 0.90
s5 Consequences of minimal cooling are explored through extensive variations of these parameters. method 0.84 0.84 0.03 asserted 0.84
s6 Results are compared with the inferred properties of thermally-emitting neutron stars in order to ascertain if enhanced cooling occurs in any of them. method 0.90 0.83 0.05 asserted 0.65
s7 All stars for which thermal emissions have been clearly detected are at least marginally consistent with the lack of enhanced cooling. result 0.94 0.83 0.51 hedged possibility 0.90
s8 The two pulsars PSR 0833-45 (Vela) and PSR 1706-44 would require enhanced cooling in case their ages and/or temperatures are on the lower side of their estimated values whereas the four stars PSR 0656+14, PSR 1055-52, Geminga, and RX J0720.4-3125 may require some source of internal heating in case their age and/or luminosity are on the upper side of their estimated values. result 0.66 0.77 0.08 hedged possibility 0.91
s9 The new upper limits on the thermal luminosity of PSR J0205+6449 and RX J0007.0+7302 are indicative of the occurrence of some enhanced neutrino emission beyond the minimal scenario. result 0.81 0.74 0.13 hedged possibility 0.96

What kind of paper the classifier reports this is: theoretical result 0.63 · evidence from analysis of existing data 0.62.

What was examined, and what was not

sources read
  • arXiv listing abstract v2
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. — Inferred temperature T∞ (top) and lu￾minosity L∞ (bottom) versus age for neutron stars with thermal emission. Data from
Figure 2. One feature notable in Figures 1 and 2 is the sizes of the error boxes, particularly in the age di￾mension. These uncertainties represent an inher￾ent difficulty in using these observations to firmly constrain the details of neutron star cooling. For this reason, instead of attempting to detail prop￾erties of the equation of state, superconductiv￾ity, and/or neutrino emissivities from the observa￾tions, our approach will be to model a reasonably broad range of acceptable physical inputs in order to determine ranges of parameters that might be excluded by the present data. 3. INPUT…
Figure 3. — Symmetry energy (top panel) and proton
Figure 5. — Mass versus radius (top panel) and ver
Figure 6. — Fermi momenta of neutrons (upper curves) and protons (lower curves) versus baryon density for the APR and BPAL EOS’s (BPAL21 and BPAL31 give nearly the same values). BPAL21 BPAL31 APR WFF3 n p
Figure 7. — Ratio of the Landau effective mass m∗ to the vacuum mass m for neutrons and protons vs. density for the four EOS’s. whose momenta lie close to the Fermi momentum: this is the Cooper theorem (Cooper 1956). As a result of this instability, the ground state of the system is reorganized with a gap ∆ in the energy spectrum around the value of the Fermi energy; no particle can have an energy between EF − ∆ and EF + ∆. This instability usually appears as a second or￾der phase transition, the gap acting as the or￾der parameter of the transition which has a corre￾sponding critical temperatu
Figure 8. — Neutron 1S0 pairing critical temper￾ature Tc vs neutron fermi momentum kF from the calculations of Ainsworth, Wambach & Pines (1989) (labeled as “AWP II” & “AWP III”: two slightly different results), Wambach, Ainsworth & Pines (1993) (“WAP”), Chen, et al. (1993) (“CCDK”), Shulze et al. (1996) (“SCLBL”), and Schwenk, Friman & Brown (2003) (“SFB”). The dotted curve shows the results of Shulze et al. (1996) in the case where medium polarization is not included. Medium polarization effects reduce the 1S0 gap by about a factor three. The vertical dotted line shows the location of the c
Figure 9. — Proton 1S0 pairing critical temperature Tc vs proton fermi momentum kF from the cal￾culations of Takatsuka (1973) (labeled as “T”), Amundsen & Østgaard (1985a) (“AO”), Chao, Clark & Yang (1972) (“CCY ms” & “CCY ps”: two slightly different results), Chen, et al. (1993) (“CCDK”), Baldo et al. (1992) (“BCLL”), El￾garøy et al. (1996) (“EEHO”), Niskanen & Sauls (1981) (“NS”). Only the calculation of Niskanen & Sauls (1981) include medium polarization. The shaded region shows the estimates of Ainsworth, Wambach & Pines (1991) of the range of values in which Tc should lie due to medium po
Figure 10. — Neutron 3P2 −3F2 pairing critical tem￾perature Tc vs neutron fermi momentum kF . The three thin continuous lines show the results of: (1) Amundsen & Østgaard (1985b), (2) Takat￾suka (1972b), and (3) Elgarøy et al. (1996). The three dotted lines show results obtained assum￾ing m∗ n = mn to illustrate the strong reduction due to small m∗ n : (0) Hoffberg et al. (1970), (1) Amundsen & Østgaard (1985b), and (2) Takat￾suka (1972b). The three thick continuous lines, (a), (b), and (c), bracket the results of Baldo et al. (1998). The vertical dotted line shows the lo￾cation of the…
Figure 12. — Specific heat of neutrons in the core of a 1.4 M⊙ star built using the APR EOS, at six different temperatures, from 0.99 down to 0.1 times 109 K illustrating the effects of pairing. No gap is present for results shown by the dotted lines. The n 3P2 gap ”a” is assumed for continu￾ous curves. This gap has a maximum Tc of 1×109 K at ρ = 3.61 × 1014 gm cm−3 . the continuous curves show the jump of cv by a factor 2.188 (see
Figure 13. — Control functions for pairing effects in
Figure 14. One must finally mention that our calculations are based on the assumption of spherical symme￾try in the interior and that the only asymmetries considered, due to the presence of a magnetic field, are within the envelope and hence included into this outer boundary condition. However, this as￾sumption is questionable in some magnetic field configurations where the field is confined to the stellar crust. As shown by Geppert, Kuecker & Page (2004), the crust is highly non-isothermal in such cases and this can affect the thermal evo￾lution because the resulting photon luminosity is lowe
Figure 15. — Neutrino and photon cooling eras for two models of non magnetized envelopes formed by heavy iron-like elements (labeled “H”) and a maximum amount of light elements (labeled “L”). The effective temperature (panel A), the cen￾tral temperature (panel B) and neutrino (contin￾uous lines) and photon (dotted lines) luminosities (panel C), all redshifted to infinity, are shown as a function of time. Pairing effects are not included in these calculations. H H −15 −13 −11 −9 Log M = −7 ∆ L L L
Figure 16. — Effect on the cooling of various amounts
Figure 18. — Comparison of the cooling of a 1.4 M⊙ star, built using the EOS of APR, without and with nucleon pairing. In the model with pair￾ing neutrino emission by the PBF process has been either artificially turned off (continuous line) or al￾lowed (dotted line). Neutron 1S0 pairing is from AWP, 3P2 pairing from our case “c” and p 1S0 pairing from AO, as labeled in Figures 8, 10, and 9, respectively. The envelope is assumed to be composed of heavy elements. Without pairing 0 t 0−ν tν−γ time Temperature With pairing T
Figure 19. — Schematic representation of the power
Figure 20. — Upper panel: neutrino luminosities vs temperature from the PBF process for four dif￾ferent n 3P2 gaps labeled as in
Figure 21. — Upper panel: neutrino luminosities vs temperature from the PBF process for 3 different p 1S0 gaps labeled as in
Figure 22. — Comparison of luminosities from various processes during th
Figure 23. — Comparison of cooling trajectories with
Figure 24. — Cooling with adjusted modified Urca, for F = 1/10, 1, 10 and 100 as indicated, with and without nucleon pairing. Assumed pairing gaps are from our model “a” for n 3P2 and from AO for p 1S0 (and n 1S0 pairing from AWP for which the effect is very small). The envelope is assumed to be composed of heavy elements. 30
Figure 25. confirms that in the absence of pair￾ing, there is almost no mass effect, both during the neutrino and the photon cooling era. Simi￾larly, when n, but not p, pairing is included, the mass dependence is also small, though larger than with no pairing at all. When the p gap is included, the main variation with mass occurs in the pho￾ton cooling era in which more massive stars cool more slowly. This is a direct consequence of the lesser suppression of the proton specific heat with increasing mass, since the p 1S0 gap vanishes at high density and there is an increasingly larger…
Figure 26. — Effects of the EOS: cooling of 1.4 M⊙ stars built using the four chosen EOS’s. EOS’s are labeled in the cases with pairing, whereas for similar cases without pairing the tra￾jectories are too similar to be separately labeled. Pairing gaps: n 1S0 from AWP, n 3P2 from our model “a” and p 1S0 from AO. The envelope is assumed to be composed of heavy elements. 5 4 2 3 1 6
Figure 27. — Cooling of a 1.1 M⊙ star (continuous lines) and with various p 1S0 gaps (as labeled in the inset): 1 - NS, 2 - T, 3 - AO, 4 - BCLL, 5 - CCY ms, and 6 - CCDK (see
Figure 28. — Comparison of predictions of the minimal cooling scenario with
Figure 29. — Same as the central panels of Figure 28, but for a 1.1
Figure 30. — Global comparison of the upper lu￾minosity limits for sources lacking apparent ther￾mal emissions with cooling trajectories, satisfy￾ing the minimal cooling scenario. The area with the darkest gray shading contains models with heavy-element dominated envelopes, whereas that with the intermediate shading contains models with ligh-element envelopes. The region with the lightest shading contains models with intermedi￾ate compositions. 40

Showing 27 of the paper's 30 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:astro-ph/0403657 by the authors.

References
How later papers used it

140 later papers cite it

Evidence: substantial. 129 of the citing papers share no author with it, and 15 of those carry a confident label from Pith's citation classifier.

Pith read 19 of the 233 citing sentences, choosing only those whose citation markers point at this paper alone, and wrote the paraphrases from them. The summary, the counts, the evidence level, the links and the quotes come from code; the summary counts the papers listed below. A later paper is listed under a label when Pith's citation classifier gives one of its sentences that label at probability 0.8 or more (0.95 for support), Pith's reading finds that the citing paper itself did it, and a second check confirms the paraphrase against the sentence. This section lists uses only. It does not list criticism or disagreement, so their absence here says nothing about how the paper has been received. Papers that share an author with this one are listed apart, under "By the same authors". A quote is the citing authors' own sentence, shown when their paper's licence is CC BY, CC BY-SA or CC0. JSON

Builds on it 1 paper · 1 by the same authors

By the same authors

  • Went beyond the minimal cooling model first envisaged in this paper by varying the equation of state and including the direct Urca process, allowing each neutron star's mass to vary.

    Simultaneous Fitting of Neutron Star Structure and Cooling Data (2018) (same authors) · 0.97 · paraphrased: the licence does not allow quoting

Uses its method or tool 2 of 9 papers with this label · 1 by the same authors

By the same authors

Uses its data or values 3 of 7 papers with this label
How the literature uses it

73 arXiv papers cite this in 117 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers whose main result relates to this one's

Found among the papers closest to this one by title and main claim. Pith's classifier compared each pair's main results and a second model family confirmed the relation; checked against two further families, shown conflicts agreed 118 times in 120 and shown same results 218 times in 239. The quoted sentence is the other paper's main claim.

States the same result 1 paper

Papers citing it · 90 sentences without a confident label

Uses its method or tool 3 papers
Uses its data or values 3 papers
  • The triple of neutron $^1S_0$ and $^3P_2$-$^3F_2$ and proton $^1S_0$ gaps were fixed to the value “WAP-b-T73” of the NSCool code, where the acronyms refer to Ref. [ref] (WAP), model b of Ref. [ref] (b), and Ref. [ref] (T73).

    Axion cooling of neutron stars. II. Beyond hadronic axions · Results of simulations · 0.86

  • Meanwhile, the coldest NS on record, PSR J2144$-$3933, has its effective temperature bounded at $<$ 33000 Kelvin by Hubble Space Telescope (HST) observations [ref] ; its spin-down age is $3×10^8$ years, implying that if passive cooling were in place its temperature would be a few hundred Kelvin in “minimal cooling" models [ref] , those that incorporate only modified Urca processes of neutrino cooling.

    Reheated Sub-40000 Kelvin Neutron Stars at the JWST, ELT, and TMT · Introduction · 0.90

  • The energy deposited in the neutron star is converted into a kinetic heating of the apparent blackbody temperature [equation] A $𝒪(10^9 year)$-old neutron star is expected to cool to $𝒪(100 K)$ [ref] .

    Dark Kinetic Heating of Neutron Stars from Contact Interactions with Relativistic Targets · Kinetic Heating · 0.88

Cites it as background 18 papers

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.