Connect internal cooling to an observed thermal spectrum

After early thermal relaxation, the interior of a static star can be approximately isothermal in the redshifted temperature $\widetilde T=e^{\Phi(r)}T(r)$, not in the local temperature itself. An integrated relativistic cooling model uses that temperature, a compatible total heat capacity, and luminosities referred to a distant observer. Before relaxation, radial heat transport must be solved rather than replacing the entire star by one temperature. Rotation and strong fields can further complicate that simplification.

An energy balance includes neutrino loss, photon loss, and any heating. Accretion, crustal reactions, magnetic dissipation, or rotational processes can add heat and change the interpretation of a cooling curve. If a toy model uses heat capacity proportional to temperature and neutrino luminosity proportional to a power q of temperature, dividing the balance gives a temperature derivative with power q minus one. The coefficient contains microphysics and structure; the exponent alone does not predict an absolute age.

For degenerate normal nucleonic matter, direct Urca has a leading emissivity scaling roughly as $T^6$ when its Fermi-momentum triangle is open. Modified Urca includes an additional spectator nucleon to conserve momentum and typically scales as $T^8$. The coefficients and thresholds depend on composition, interactions, and effective masses. These power laws refer to particular low-temperature regimes and should not be compared as if they had identical normalizations or occupied the same core volume.

Pairing reduces the available quasiparticle phase space and changes heat capacity and neutrino rates. Pair formation and breaking can create an additional neutrino channel near a superfluid transition, so saying that superfluidity only suppresses cooling is incomplete. The size and density dependence of gaps are uncertain. An insulating envelope with light elements or a magnetic field can change the mapping from interior to surface temperature even with identical core neutrino physics.

For a static spherical exterior with the surface outside the photon sphere, local surface temperature and luminosity are redshifted on their way to infinity. The apparent radiation radius is larger than the areal surface radius. The redshift relations apply to a consistently defined bolometric thermal emission model; an observed color temperature still needs atmospheric and instrumental modeling. A fitted blackbody normalization is therefore not automatically a geometric measurement of the neutron-star surface.

As emphasized in the primary cooling analysis by Yakovlev and Pethick, comparison with observed thermal emission requires the envelope composition, age information, mass-dependent reaction thresholds, and pairing assumptions. The physical framework remains useful, while its 2004 examples should be recognized as dated observations rather than a catalog of the latest measured stars. Modern data improve the likelihoods but do not remove the structural degeneracy between heat transport, neutrino processes, heating, and atmospheric inference.