Debye cooling of ultramassive white dwarfs

Abstract

Recent observations of cold and massive white dwarfs are beginning to probe Debye cooling: an accelerated cooling phase driven by the vanishing heat capacity $C_V(T)$ of the crystallized interiors at low temperatures $T$, due to the quantum nature of phonon oscillations. We analyse this regime by calculating, both analytically and numerically, the instantaneous cooling time $\tau (L)\equiv -\text{d}t/\text{d}\ln L$, which is a proxy for the number density of white dwarfs at a given luminosity $L$. We find that cooling accelerates ($\tau$ shortens) at $L\lesssim 10^{-4}{\text{L}}_{\odot }$ due to a combination of Debye cooling and the simultaneous convective coupling of the degenerate core to the surface, with neither process causing significant acceleration on its own. Specifically, strong magnetic fields can inhibit convective coupling, leading to an almost flat $\tau \propto L^{1/7}$ instead. We parametrize the uncertainty in the solid microphysics by generalizing the Helmholtz free energy such that $C_V\propto T^{\delta }$ at low temperatures. The currently observed decreasing luminosity function in the Debye regime is tentatively consistent with $\tau (L)$ of non-magnetic white dwarfs with $\delta =3$, corresponding to the standard Debye model. Larger samples could probe, in addition to magnetism, deviations $\delta <3$, which diverge appreciably before the Galactic disc's age $t\approx 10\text{ Gyr}$ for masses $M\gtrsim 1.0{\text{M}}_{\odot }$. For $\delta >3$, in contrast, the cooling curves converge because the heat capacity becomes dominated by the degenerate electrons $C_V\propto T$.

arXiv comment9 pages, 6 figures. Submitted to MNRAS, comments are welcome