The missing cooling delay of crystallizing white dwarfs: reconciling the observed Gaia crystallization pile-up with cooling theory

Gaia Data Release 2 (DR2) photometry and parallaxes of the white dwarfs within 100 pc revealed a pile-up in the cooling sequence at the luminosities where core crystallization sets in, providing the first direct observational evidence of crystallization in strongly coupled Coulomb plasmas. The release of latent heat and of gravitational energy from oxygen phase separation slows the cooling by roughly a billion years, but the observed pile-up is narrower and of higher magnitude than evolving-model predictions, indicating an additional cooling delay of order 1 Gyr that standard microphysics does not account for. The open problem is to close the energy budget of crystallizing white dwarfs: identify and validate the physical ingredient(s) — the leading candidate being phase separation and distillation of the neutron-rich isotope $^{22}$Ne, together with improved multicomponent phase diagrams including impurity effects — that allow evolutionary and population models to quantitatively reproduce the observed pile-up amplitude, width, and shape across the white-dwarf mass range, with direct implications for the systematic-uncertainty budget of white-dwarf cooling ages (cosmochronology).

Background

A white dwarf is the degenerate remnant of a low- or intermediate-mass star, supported against gravity by electron degeneracy pressure and devoid of nuclear energy sources, so that it simply cools over billions of years; its cooling rate therefore acts as a clock, the basis of white-dwarf cosmochronology (age-dating of stellar populations such as the Galactic disk, open and globular clusters, and the halo). The interior is a dense plasma of carbon and oxygen ions immersed in a degenerate electron gas, characterized by the Coulomb coupling parameter $\Gamma$, the ratio of the electrostatic potential energy between neighboring ions to the thermal energy. As the star cools, $\Gamma$ rises and the ion liquid is predicted to undergo a first-order phase transition to a solid (crystallization), conventionally taken to set in at $\Gamma\simeq175$. Crystallization releases latent heat; in addition, because the solid preferentially retains $^{16}$O over $^{12}$C, the growing crystal core becomes oxygen-enriched and the excess carbon is displaced upward, releasing further gravitational energy (oxygen phase separation). Both effects delay cooling by roughly a billion years, and the neutron-rich isotope $^{22}$Ne (produced during helium burning in the progenitor) can in addition sediment and release gravitational energy of its own; at sufficiently high $^{22}$Ne abundance the crystallizing core enters a distillation regime, in which the solid phase is so depleted in $^{22}$Ne that the crystals become buoyant, float up and melt, and the displaced neutron-rich liquid sinks, releasing a large additional amount of gravitational energy. Gaia's precise parallaxes and photometry yielded the first unbiased colour–absolute-magnitude diagram for the local white-dwarf population, and Tremblay et al. (2019) identified a 'transverse' sequence in the 100 pc sample — a mass-dependent pile-up (overdensity) exactly where crystallization is predicted — as direct evidence of the phase transition, requiring both latent heat and oxygen phase separation to reproduce the luminosity function. However, the authors noted that the modelled bump is imperfect: the observed feature is narrower and of higher magnitude than the models, with the discrepancy depending on uncertain inputs — the adopted coupling parameter ($\Gamma=175$), the assumed homogeneous 50/50-by-mass $^{12}$C/$^{16}$O core profile, the envelope stratification (hydrogen layer mass $M_{\rm H}/M_{\rm WD}=10^{-4}$, helium layer $M_{\rm He}/M_{\rm WD}=10^{-2}$), and possible $^{22}$Ne sedimentation. Subsequent work sharpened the discrepancy rather than removing it: even with a newly computed, accurate carbon–oxygen phase diagram, the pile-up magnitude remains underestimated by current evolution models, and at the high-mass end the so-called Q branch (an overdensity of massive white dwarfs in the Gaia Hertzsprung–Russell diagram, kinematically dated) requires a small fraction of objects to experience an extra cooling delay of several gigayears that standard sources cannot supply. Because the pile-up is the observational manifestation of the crystallization energy budget, resolving it means determining the missing energy source and the correct dense-matter microphysics (multicomponent phase diagrams, trace-impurity effects, sedimentation energies at stellar densities — conditions that cannot be reproduced in terrestrial laboratories), and propagating the answer into the systematic uncertainties of the cooling ages on which white-dwarf cosmochronology relies.

Problem Statement

Close the energy budget of crystallizing white dwarfs as tested by the Gaia observations: determine the physical ingredient(s) — beyond latent heat and $^{16}$O/$^{12}$C phase separation — that evolutionary and population-synthesis models must include so that the predicted crystallization-induced pile-up in the Gaia white-dwarf luminosity function and colour–magnitude diagram quantitatively reproduces the observed feature (its amplitude, width, and shape) across the white-dwarf mass range. The observational target is the pile-up measured for the 100-pc Gaia white-dwarf sample by Tremblay et al. (2019; ref. 1) and its published follow-up refinements; its amplitude, width, and shape in the luminosity function and colour–magnitude diagram, together with their quoted uncertainties, define the agreement sought. Concretely, the problem is to identify and validate the missing cooling-delay mechanism for standard-composition C/O white dwarfs (order 1 Gyr at crystallization; the leading candidate is phase separation and distillation of $^{22}$Ne, but alternatives — revised multicomponent phase diagrams including trace-impurity effects, envelope and convective-coupling physics — must be either incorporated or excluded) and to show, via evolutionary and population modeling compared against the observed 100-pc Gaia pile-up, that the observed amplitude, width, and shape are quantitatively reproduced, with the population-model inputs (star-formation history, initial-mass function, initial–final mass relation linking progenitor mass to remnant mass, merger fraction) either pinned independently or shown not to control the comparison. The target is a single validated physical account of the crystallization cooling delay for the standard-composition population in which the discrepancy was identified; a mechanism validated only for ultramassive merger remnants lies outside that scope.

References

1. Tremblay, P.-E., Fontaine, G., Gentile Fusillo, N. P., Dunlap, B. H., Gänsicke, B. T., Hollands, M. A., Hermes, J. J., Marsh, T. R., Cukanovaite, E., & Cunningham, T. (2019). Core crystallization and pile-up in the cooling sequence of evolving white dwarfs. Nature 565, 202–205. doi:10.1038/s41586-018-0791-x, arXiv:1908.00370. https://doi.org/10.1038/s41586-018-0791-x
2. Cheng, S., Cummings, J. D., & Ménard, B. (2019). A cooling anomaly of high-mass white dwarfs. The Astrophysical Journal 886, 100. doi:10.3847/1538-4357/ab4989, arXiv:1905.12710. https://doi.org/10.3847/1538-4357/ab4989
3. Blouin, S., Daligault, J., Saumon, D., Bédard, A., & Brassard, P. (2020). Toward precision cosmochronology: a new C/O phase diagram for white dwarfs. Astronomy & Astrophysics 640, L11. doi:10.1051/0004-6361/202038879, arXiv:2007.13669. https://doi.org/10.1051/0004-6361/202038879
4. Bauer, E. B., Schwab, J., Bildsten, L., & Cheng, S. (2020). Multi-gigayear white dwarf cooling delays from clustering-enhanced gravitational sedimentation. The Astrophysical Journal 902, 93. doi:10.3847/1538-4357/abb5a5, arXiv:2009.04025. https://doi.org/10.3847/1538-4357/abb5a5
5. Caplan, M. E., Horowitz, C. J., & Cumming, A. (2020). Neon cluster formation and phase separation during white dwarf cooling. The Astrophysical Journal Letters 902, L44. doi:10.3847/2041-8213/abbda0, arXiv:2010.00036. https://doi.org/10.3847/2041-8213/abbda0
6. Blouin, S., Daligault, J., & Saumon, D. (2021). $^{22}$Ne phase separation as a solution to the ultramassive white dwarf cooling anomaly. The Astrophysical Journal Letters 911, L5. doi:10.3847/2041-8213/abf14b, arXiv:2103.12892. https://doi.org/10.3847/2041-8213/abf14b
7. Blouin, S., & Daligault, J. (2021). Phase separation in ultramassive white dwarfs. The Astrophysical Journal 919, 87. doi:10.3847/1538-4357/ac1513, arXiv:2107.07094. https://doi.org/10.3847/1538-4357/ac1513
8. Blouin, S., & Daligault, J. (2021). Direct evaluation of the phase diagrams of dense multicomponent plasmas by integration of the Clapeyron equations. Physical Review E 103, 043204. doi:10.1103/PhysRevE.103.043204, arXiv:2104.00599. https://doi.org/10.1103/PhysRevE.103.043204
9. Saumon, D., Blouin, S., & Tremblay, P.-E. (2022). Current challenges in the physics of white dwarf stars. Physics Reports 988, 1–63. doi:10.1016/j.physrep.2022.09.001, arXiv:2209.02846. https://doi.org/10.1016/j.physrep.2022.09.001
10. Bauer, E. B. (2023). Carbon–oxygen phase separation in Modules for Experiments in Stellar Astrophysics (MESA) white dwarf models. The Astrophysical Journal 950, 115. doi:10.3847/1538-4357/acd057, arXiv:2303.10110. https://doi.org/10.3847/1538-4357/acd057
11. Shen, K. J., Blouin, S., & Breivik, K. (2023). The Q branch cooling anomaly can be explained by mergers of white dwarfs and subgiant stars. The Astrophysical Journal Letters 955, L33. doi:10.3847/2041-8213/acf57b, arXiv:2308.04559. https://doi.org/10.3847/2041-8213/acf57b
12. Barrientos, M., Kilic, M., Bergeron, P., Blouin, S., Brown, W. R., & Andrews, J. J. (2024). Fundamental tests of white dwarf cooling physics with wide binaries. The Astrophysical Journal 973, 88. doi:10.3847/1538-4357/ad6303, arXiv:2407.18763. https://doi.org/10.3847/1538-4357/ad6303
13. Bédard, A., Blouin, S., & Cheng, S. (2024). Buoyant crystals halt the cooling of white dwarf stars. Nature 627, 286–288. doi:10.1038/s41586-024-07102-y, arXiv:2409.04419. https://doi.org/10.1038/s41586-024-07102-y
14. Montgomery, M. H., & Dunlap, B. H. (2024). Fluid mixing during phase separation in crystallizing white dwarfs. The Astrophysical Journal 961, 197. doi:10.3847/1538-4357/ad16dc, arXiv:2312.11647. https://doi.org/10.3847/1538-4357/ad16dc
15. Tremblay, P.-E., Bédard, A., O'Brien, M. W., Munday, J., Elms, A. K., Gentile Fusillo, N. P., & Sahu, S. (2024). The Gaia white dwarf revolution. New Astronomy Reviews 99, 101705. doi:10.1016/j.newar.2024.101705, arXiv:2402.14960. https://doi.org/10.1016/j.newar.2024.101705
16. Pathak, P., Blouin, S., & Herwig, F. (2025). Quantifying systematic uncertainties in white dwarf cooling age determinations. The Astrophysical Journal 983, 158. doi:10.3847/1538-4357/adc10f, arXiv:2410.14014. https://doi.org/10.3847/1538-4357/adc10f
17. Castro-Tapia, M., & Cumming, A. (2025). Three-component phase separation for ultramassive white dwarf models. The Astrophysical Journal 991, 64. doi:10.3847/1538-4357/adf745, arXiv:2506.16529. https://doi.org/10.3847/1538-4357/adf745
18. Bauer, E. B., Dotter, A., Conroy, C., Cunningham, T., Park, M., & Tremblay, P.-E. (2026). MESA Isochrones and Stellar Tracks (MIST). III. The white dwarf cooling sequence. The Astrophysical Journal Supplement Series 283, 41. doi:10.3847/1538-4365/ae401e, arXiv:2509.21717. https://doi.org/10.3847/1538-4365/ae401e
19. Xu, M., Smith, A. A., Prokhorov, L., Guarrera, V., & Barontini, G. (2026). Measuring impurity-induced shifts in Coulomb crystallization. arXiv:2603.14980. https://arxiv.org/abs/2603.14980