Neutron star heating vs. HST observations

This paper demonstrates that no single heating mechanism can explain the unexpectedly high surface temperatures of several old neutron stars, but a combined model of rotochemical heating and vortex creep successfully reproduces the observations for PSR J0437−4715 and PSR B0950+08 while remaining consistent with upper limits for other sources.

The Gist

Imagine a neutron star as a cosmic campfire. When it's first born in a supernova explosion, it's a roaring inferno, blazing at temperatures hotter than the center of the sun. But like any fire, it's supposed to run out of fuel and die down. According to standard physics, once a neutron star gets old (billions of years), it should have cooled down so much that it's practically invisible to our telescopes—like a cold, dead ember that no longer glows.

However, astronomers using the Hubble Space Telescope looked at five very old neutron stars and found something strange: four of them were still glowing with a warm, ultraviolet light. They were too hot to be just "dead embers." This paper asks: What is keeping these cosmic campfires warm?

The authors tested three different "heaters" that might be working inside these stars, and then combined them to see if they could explain the observations. Here is the breakdown using simple analogies:

The Three Potential Heaters

1. Rotochemical Heating (The "Squeezed Spring"):
   As a neutron star spins, it bulges at the equator. As it slows down over millions of years, it gets slightly more spherical. This change squeezes the star's core, changing the pressure. Imagine a spring that is slowly being compressed; eventually, the pressure builds up until it snaps back, releasing energy. In the star's core, this "snap" triggers nuclear reactions that release heat.

   • The Catch: For this to work efficiently, the star needs to be spinning very fast initially, and the particles inside need to be in a special "superfluid" state (like a frictionless liquid). If the particles are in this state, they act like a dam, holding back the reactions until the pressure gets huge, then releasing a massive burst of heat.

2. Vortex Creep (The "Rubbing Hands"):
   Inside the star's crust, there is a superfluid that spins faster than the solid crust on the outside. As the star slows down, the superfluid tries to keep spinning, creating tiny whirlpools (vortices). These whirlpools get stuck on the crust's atomic lattice, like a gear getting stuck in a machine. Eventually, they slip and slide, creating friction.

   • The Analogy: Think of rubbing your hands together to generate heat. The friction between the spinning superfluid and the solid crust generates warmth. This depends heavily on how fast the star is slowing down right now.

3. Crustal Heating (The "Compressed Squeezie"):
   Some neutron stars (called millisecond pulsars) were "rejuvenated" by stealing matter from a companion star. This extra weight squashed the star's crust. As the star continues to spin down, the crust compresses even more, triggering nuclear reactions deep inside the rock-like layers.

   • The Catch: The authors found this heater is too weak to explain the warmth of the hottest stars they observed.

The Great Detective Work

The team ran computer simulations to see which heater (or combination) could explain the temperatures of the five specific stars they observed:

• PSR J0437−4715: A very old, fast-spinning star that is surprisingly hot.
• PSR B0950+08: An old, slower-spinning star that is also warm.
• Three others: Stars that were not detected, meaning they are very cold (or at least, colder than a certain limit).

The Results:

• No single heater worked for everyone.
  • If you used only the "Rubbing Hands" (Vortex Creep) heater, you could explain the warmth of the slow star (B0950), but it wasn't strong enough to heat up the fast star (J0437).
  • If you used only the "Squeezed Spring" (Rotochemical) heater with the special "superfluid" conditions, you could explain the fast star (J0437), but it required the slow star to have started spinning impossibly fast in the past, which doesn't fit the data.

The Winning Combination:
The authors found that you need both heaters working together to explain the whole picture:

1. For the fast star (J0437): The "Squeezed Spring" (Rotochemical heating) is the main driver. The star must have started spinning incredibly fast (faster than a millisecond) and has a special internal structure (large energy gaps in the superfluid) that allows it to store up heat and release it now.
2. For the slow star (B0950): The "Rubbing Hands" (Vortex Creep) is the main driver. The friction from the slowing spin keeps it warm.
3. For the others: This combined model predicts that the three non-detected stars should be just barely cold enough to be invisible, but very close to the limit of detection.

The Bottom Line

The paper concludes that neutron stars aren't just cooling down passively. They are complex machines where different internal "engines" kick in depending on how fast they spin and what their internal ingredients are. To explain why some old stars are still glowing, we need a mix of friction from spinning and pressure-induced nuclear reactions, provided the star started its life spinning at a breakneck speed.

The authors suggest that if we look at these stars again with more sensitive telescopes, we should find that the "invisible" ones are actually just barely glowing, confirming this dual-heater theory.

Technical Summary

1. Problem Statement

Standard cooling models predict that isolated neutron stars (NSs) should cool rapidly via neutrino emission in their first $10^5$ years and subsequently via photon emission, reaching surface temperatures ($T_s$) below $10^4$ K within $10^7$ years, rendering them undetectable by current telescopes. However, Hubble Space Telescope (HST) observations have detected thermal ultraviolet (UV) emission from several old pulsars with characteristic ages ($\tau$) significantly exceeding this limit:

• Millisecond Pulsars (MSPs): PSR J0437−4715 ($\tau > 10^9$ yr) and PSR J2124−3358 ($\tau > 10^9$ yr).
• Classical Pulsars (CPs): PSR B0950+08 ($\tau \sim 10^7$ yr) and PSR J0108−1431 ($\tau \sim 10^8$ yr).

These detections imply surface temperatures of $T_s \sim 10^5$ K, suggesting the presence of internal heating mechanisms that counteract passive cooling. Conversely, the non-detection of PSR J2144−3933 provides a stringent upper limit on heating efficiency. The central problem is to identify which heating mechanisms (or combinations thereof) can simultaneously explain the high temperatures of specific pulsars while respecting the upper limits of others.

2. Methodology

The authors performed numerical simulations of the thermal evolution of five specific pulsars, comparing theoretical cooling/heating curves against HST observational data.

Observational Data:

• Updated surface temperature measurements and upper limits were derived from HST (and ROSAT for J0437) data, utilizing blackbody and hydrogen atmosphere models.
• Key parameters included rotation period ($P$), magnetic field ($B$), and characteristic age ($\tau$).

Physical Models:
The thermal evolution was computed by integrating the energy balance equation:
$$\dot{T}^\infty_c = \frac{1}{C} (L^\infty_H - L^\infty_\nu - L^\infty_\gamma)$$
where $L^\infty_H$ is heating power, $L^\infty_\nu$ is neutrino luminosity, and $L^\infty_\gamma$ is photon luminosity. The study focused on three specific heating mechanisms:

1. Rotochemical Heating:

   • Mechanism: Continuous spin-down creates chemical imbalances ($\eta$) in the core, driving non-equilibrium beta reactions (Urca processes) that release heat.
   • Variations: Modeled for Normal Matter (Modified Urca and Direct Urca) and Cooper-paired Matter (Superfluid neutrons/superconducting protons).
   • Key Feature: In Cooper-paired matter, reactions are suppressed until the chemical imbalance exceeds a threshold energy gap ($\Delta_{thr}$). Once exceeded, heating is significantly enhanced. The study assumed uniform, isotropic gaps and introduced reduction factors to mimic non-uniform gaps and suppress thermal oscillations.

2. Vortex Creep:

   • Mechanism: Friction between quantized vortices in the superfluid inner crust and the nuclear lattice as the star spins down.
   • Parameter: Treated the excess angular momentum ($J$) as a free parameter, constrained by the non-detection of PSR J2144.

3. Crustal Heating:

   • Mechanism: Pycnonuclear reactions in the deep crust of recycled MSPs triggered by compression due to spin-down.
   • Constraint: Modeled as an upper limit based on original predictions, acknowledging recent Equation of State (EoS) updates might reduce this effect.

Simulation Setup:

• Initial Conditions: MSPs started with $P_0 = 1$ ms and $T_0 = 10^9$ K; CPs started with $P_0 = 5$ ms and $T_0 = 10^{11}$ K.
• Spin-down: Assumed constant magnetic dipole braking.
• Scenarios: The authors tested 11 distinct cases (Table 3) combining different reaction types (Murca/Durca), matter states (Normal/Superfluid), and heating mechanisms.

3. Key Contributions

• Comprehensive Comparison: This is the first study to simultaneously model and contrast all three major heating mechanisms (rotochemical, vortex creep, crustal) against a unified set of HST observations for both MSPs and CPs.
• Refined Rotochemical Modeling: The authors rigorously applied the effects of Cooper pairing with large energy gaps ($\sim 1.5$ MeV) to rotochemical heating, demonstrating how this mechanism can be "switched on" only for specific initial rotation conditions.
• Oscillation Suppression: The study addresses the issue of thermal/chemical oscillations in models with large uniform gaps by introducing reaction rate reduction factors, allowing for stable quasi-steady-state predictions.
• Joint Mechanism Proposal: The paper proposes that a single heating mechanism is insufficient for the entire sample, necessitating a hybrid model.

4. Results

The simulations yielded the following findings:

• Passive Cooling: Fails completely; temperatures drop below $10^4$ K, inconsistent with observations.
• Rotochemical Heating (Normal Matter):
  • Can explain the temperatures of CPs (B0950, J0108) via Modified Urca reactions.
  • Fails to explain the high temperature of MSP J0437.
• Rotochemical Heating (Cooper-paired Matter):
  • With a large gap ($\Delta_{thr} \approx 1.5$ MeV), this mechanism successfully explains the high temperature of PSR J0437.
  • Requirement: Requires an initial rotation period $P_0 \lesssim 1.8$ ms. This is plausible for MSPs (like J0437) but implausible for CPs (which typically have $P_0 > 10$ ms), meaning this mechanism cannot explain the heating of B0950.
• Vortex Creep:
  • With an excess angular momentum parameter $J \sim 3 \times 10^{43}$ erg s, this mechanism explains the temperature of PSR B0950.
  • Fails to explain the high temperature of J0437 (the heating power is insufficient).
• Crustal Heating:
  • Provides a heating contribution for MSPs but is too weak to account for the observed temperature of J0437 on its own.
• The Hybrid Solution (Case XI):
  • A model combining Rotochemical Heating (with Cooper pairing, $\Delta \approx 1.5$ MeV) and Vortex Creep ($J \approx 3 \times 10^{43}$ erg s) successfully reproduces the observations for all five pulsars.
  • Mechanism: Rotochemical heating dominates in MSPs (activating only if $P_0$ is fast enough), while Vortex Creep dominates in CPs.
  • Prediction: This model predicts that the temperatures of the non-detected or marginally detected pulsars (J2124, J0108, J2144) should be close to their current observational upper limits.

5. Significance

• Equation of State Constraints: The success of the rotochemical heating model for J0437 implies the existence of a large Cooper pairing gap ($\sim 1.5$ MeV) in the NS core, providing a rare observational constraint on the superfluid properties of dense matter.
• Initial Conditions: The results suggest that the initial rotation period of recycled pulsars (MSPs) was likely very fast ($P_0 \lesssim 2$ ms), consistent with binary evolution models but distinct from classical pulsars.
• Future Observations: The paper provides a clear, testable prediction: deeper UV observations or observations with wider wavelength coverage for PSR J2124, PSR J0108, and PSR J2144 should detect thermal emission near their current upper limits. Confirming this would strongly validate the combined heating model and the specific physics of the NS interior.
• Exclusion of Alternatives: The study effectively rules out single-mechanism explanations (e.g., magnetic field decay or dark matter accretion) for the entire sample, narrowing the focus to internal thermodynamic processes driven by spin-down.