Neutron Star Cooling Timescale Calculator
Estimate neutron star core and surface temperature from age, or age from an observed surface temperature, using simplified Urca cooling laws.
🌡️ What is a Neutron Star Cooling Timescale Calculator?
A neutron star cooling timescale calculator estimates how a neutron star's internal and surface temperature change as it ages, using simplified physical models of the dominant heat-loss processes. A newborn neutron star emerges from core-collapse supernova at a blistering internal temperature above 10^11 K, but within seconds to minutes it cools to around 10^10 K, and over the following years to hundreds of thousands of years, neutrino emission from deep inside the core drains its thermal energy far faster than photon radiation from the surface ever could. This calculator implements the two standard neutrino cooling channels, modified Urca and direct Urca, and the envelope relation that links core temperature to the surface temperature an X-ray telescope would actually observe.
The modified Urca process is the default, slower cooling channel present in essentially every neutron star. It requires a bystander nucleon to conserve momentum during the neutron-to-proton beta decay reaction that produces a neutrino, which suppresses the reaction rate and gives a comparatively gentle temperature decline, roughly T_core proportional to t^-1/6. The direct Urca process is a much more efficient neutrino-producing reaction that only becomes kinematically allowed when the local proton fraction exceeds a density-dependent threshold, typically reached only in the cores of the most massive neutron stars. When active, it cools the star far faster, roughly T_core proportional to t^-1/4, and can make a massive neutron star effectively invisible in thermal X-rays within a few thousand years.
A common point of confusion is treating core temperature and observed surface temperature as the same number. They are not: heat has to conduct outward through a solid crust and a thin, insulating envelope before it can radiate away as thermal photons, so the surface temperature an observer measures is typically one to two orders of magnitude lower than the core temperature. This calculator applies a simplified envelope scaling relation, in the spirit of the classic Gudmundsson, Pethick and Epstein (1983) and Potekhin, Chabrier and Yakovlev (1997) envelope models, to connect the two.
This tool is useful for astronomy students studying compact object thermal evolution, for building intuition about why some neutron stars (like the Cassiopeia A neutron star) match textbook cooling curves closely while others run unexpectedly hot or cold, and for quickly estimating the expected X-ray energy range of a neutron star's thermal emission given its age or a candidate age given its measured temperature.
📐 Formula
📖 How to Use This Calculator
Steps
💡 Example Calculations
Example 1 - Cassiopeia A Neutron Star
Age mode: t = 330 years, Modified Urca (Slow)
Example 2 - Vela Pulsar
Age mode: t = 11,000 years, Modified Urca (Slow)
Example 3 - Massive Neutron Star, Direct Urca Contrast
Age mode: t = 1,000 years, Direct Urca (Fast)
Example 4 - Age from an Observed Temperature
Temperature mode: Teff = 1,950,000 K, Modified Urca (Slow)
❓ Frequently Asked Questions
🔗 Related Calculators
How does a neutron star cool over time?
A young neutron star cools mainly by emitting neutrinos from its core (Urca processes), which dominate for roughly the first 10^5 years. Later, photon emission from the surface takes over as the primary cooling channel. This calculator models the early neutrino-dominated era using simplified modified and direct Urca power laws.
What is the difference between modified and direct Urca cooling?
Modified Urca is the standard, slower neutrino cooling process available in essentially all neutron stars; it involves an extra bystander nucleon and gives core temperature T_core proportional to t^-1/6. Direct Urca is a much faster process that only operates when the local proton fraction exceeds a density-dependent threshold, typically only in the cores of the most massive neutron stars; it gives T_core proportional to t^-1/4, cooling the star far more quickly.
Why is the surface temperature lower than the core temperature?
Heat generated or stored in the neutron star core must conduct outward through a solid crust and thin envelope, roughly a few hundred meters thick, before it can radiate as thermal photons at the surface. This envelope acts as a thermal insulator, so the observed surface temperature is typically 10 to 100 times lower than the core temperature.
What is the significance of the Cassiopeia A neutron star?
The Cassiopeia A neutron star formed in a supernova around 1680 CE, making it one of the youngest and best-dated known neutron stars. Chandra X-ray Observatory monitoring since 1999 detected a measurable decline in its surface temperature, providing some of the first direct observational evidence for the onset of core neutron superfluidity accelerating neutrino cooling.
How does this calculator estimate age from a measured temperature?
The reverse mode inverts the same power-law relations used in the forward mode: it first converts the observed surface temperature to an implied core temperature using the envelope relation, then inverts the cooling law (t proportional to core temperature to the power of -6 for modified Urca, or -4 for direct Urca) to solve for age.
How accurate is this calculator?
This calculator uses simplified, order-of-magnitude power-law approximations calibrated to reproduce published thermal evolution curves for standard neutron star cooling, in the spirit of the review by Yakovlev and Pethick (2004). It captures the correct scaling and gives realistic ballpark numbers, but real cooling curves depend on the equation of state, magnetic field, pairing gaps, and envelope composition, which detailed cooling codes model far more precisely.
What is the peak photon energy output showing?
It applies Wien's displacement law in energy form, E_peak = 2.82 k_B T_eff, to estimate the characteristic photon energy of the neutron star's thermal blackbody spectrum. For a surface temperature around 10^6 K this falls in the soft X-ray band (roughly 100 to 500 eV), consistent with the energies detected from real thermally emitting neutron stars.
Why do some old pulsars appear hotter than the standard cooling curve predicts?
Several old pulsars show surface temperatures above the standard modified Urca prediction. Proposed explanations include rotochemical heating (frictional heating from deviations between the actual and beta-equilibrium composition as the star spins down), magnetic field decay, or crustal impurity heating. These are active areas of neutron star astrophysics research.
At what age does photon cooling take over from neutrino cooling?
For a standard modified Urca cooling neutron star, neutrino emission dominates for roughly the first 10^5 years, after which photon emission from the surface becomes the dominant energy loss channel and the temperature decline accelerates. This calculator's simplified power laws are most reliable within the earlier neutrino-dominated era.
Can this formula be used for a white dwarf?
No. White dwarfs cool almost entirely by photon emission and their thermal evolution follows a very different physical model (Mestel cooling, roughly T proportional to t^-1/2 to t^-5/7 depending on the era). Neutron star cooling is dominated by neutrino emission for most of its early life because neutron star cores are vastly denser and hotter.