Remaining systematics in white-dwarf core crystallization: quantum-ion and solid-structure effects on multicomponent plasma phase diagrams and their quantitative reconciliation with observed cooling signatures

White-dwarf core crystallization is governed by the liquid–solid phase diagram of dense carbon–oxygen plasma (with trace $^{22}$Ne and other impurities). The foundational molecular-dynamics phase diagram of Horowitz, Schneider and Berry (2010) came with four stated open issues: unquantified approximations in the interaction model (classical ions, neglected electron mass, neglected quantum effects), possible nonequilibrium bias in two-phase coexistence simulations, the unstudied role of a small neon abundance, and the need for white-dwarf luminosity functions computed with the new phase diagram. Later literature resolved the binary C/O phase diagram (three independent methods agree; azeotropic form with azeotrope near $x_O \approx 0.2$), superseded the nonequilibrium concern, and extensively mapped multicomponent and trace-species crystallization (including $^{22}$Ne distillation, which halts cooling in a delayed population). The surviving open core is: (1) a quantification of ion-quantum and solid-structure effects (short-range ordering, non-bcc and polycrystalline configurations, miscibility gaps and exsolution) on the phase diagrams and fractionation of crystallizing multicomponent plasmas, which to date has been probed only semi-analytically — first-principles quantum (path-integral) simulation exists only for the one-component plasma, not for mixtures; and (2) quantitative closure between the modern phase diagrams and observed crystallization signatures (Gaia pile-up, cluster luminosity functions), for which sub-Gyr to ~1 Gyr residual systematic cooling-age uncertainties persist (Pathak, Blouin & Herwig 2025) and the pile-up magnitude remains underestimated even with the latest phase diagram (Blouin et al. 2020).

Background

A white dwarf (WD) is the electron-degenerate remnant of a low- or intermediate-mass star; its core is a fully ionized plasma of carbon and oxygen (with traces of $^{22}$Ne, and O/Ne for the most massive remnants) at densities of about 10^6–10^7 g/cm^3. The ions in this plasma are strongly coupled: the Coulomb coupling parameter $\Gamma = Z^2 e^2/(a k_B T)$, the ratio of the electrostatic energy at the mean ion spacing $a$ to the thermal energy, exceeds ~175 at crystallization. As a WD cools, its core undergoes a first-order liquid–solid phase transition that releases latent heat and, because the solid and liquid have different C/O compositions, also gravitational energy from compositional redistribution; these effects delay WD cooling by up to ~1 Gyr and are central to the use of WDs as cosmic clocks (cosmochronology).

A phase diagram of the dense C/O mixture — the liquidus and solidus curves giving the coexisting liquid and solid compositions and the melting temperature as functions of the oxygen fraction $x_O$ — is the microphysical ingredient that determines when and how this happens. Horowitz, Schneider and Berry (2010) computed it with two-phase molecular dynamics (MD): a simulation box prepared half liquid and half solid, with ions interacting through a Yukawa potential $Z_i Z_j e^2 e^{-r/\lambda}/r$ (electron screening length $\lambda$), evolved until the two phases coexist. Their diagram is azeotropic (the liquidus and solidus touch at an azeotrope, a composition at which liquid and solid have identical composition, near $x_O \approx 0.2$) rather than spindle-shaped, the solid being oxygen-enriched. They used it, together with the observed crystallization peak in the luminosity function (number of WDs per luminosity interval) of the globular cluster NGC 6397 (Winget et al. 2009), to bound the core oxygen fraction to order 60% or less and the astrophysical $S$-factor (the low-energy normalization of the nuclear reaction cross section) of the $^{12}$C($\alpha,\gamma$)$^{16}$O reaction to $S_{300} \le 170$ keV b. The paper itself listed four unresolved issues: (1) the simulations treat the ions classically, neglect the electron mass (which shortens the screening length by ~20% at WD central densities) and neglect quantum effects, and the quantitative impact of these approximations on the phase diagram was unquantified; (2) in two-phase MD, equilibration of impurity diffusion across the liquid–solid interface is slow, so finite run times can leave residual nonequilibrium bias in the coexistence compositions; (3) the role of a small neon abundance (multicomponent effects beyond binary C/O) needed explicit three-component MD; (4) WD luminosity functions should be recomputed with this phase diagram, with careful attention to systematic errors, before the inferred limits on melting temperature, oxygen fraction and $S_{300}$ could be regarded as secure.

Subsequent work transformed this landscape. Three independent techniques now exist: semi-analytic phase diagrams built on Monte Carlo free-energy fits (Medin & Cumming 2010 lineage), two-phase MD (refined by Schneider et al. 2012 with 55296-ion boxes and an explicit nonequilibrium study), and Clapeyron (Gibbs–Duhem) integration (Blouin & Daligault 2021), which evolves the coexistence curve by numerically integrating the mixture Clausius–Clapeyron equation with Monte Carlo simulations of each phase performed separately in the isobaric semi-grand canonical ensemble (a statistical ensemble in which the total number of particles is fixed but species identity can exchange, at constant pressure). The last method avoids solid–liquid interfaces and finite-size/nonequilibrium issues and includes the electron background explicitly. All three methods agree on an azeotropic C/O diagram, and the 2022 review by Saumon, Blouin and Tremblay judges the binary C/O problem "mostly solved", identifying the remaining challenges as the classical-ion approximation (quantum effects), the assumed bcc (body-centered cubic) lattice structure of the solid, and the statistics of trace-species fractionation. Meanwhile Gaia (the European Space Agency's precision astrometry satellite) detected the crystallization pile-up — an excess accumulation of white dwarfs at the luminosities where crystallization delays cooling — in the field WD cooling sequence (Tremblay et al. 2019), luminosity-function peaks were measured in globular clusters (Bedin et al. 2019), and an ~8 Gyr extra cooling delay was found on the Q branch (an overdensity of high-mass white dwarfs on the Gaia colour–magnitude diagram) (Cheng et al. 2019), traced to $^{22}$Ne distillation — a crystallization-driven process in which $^{22}$Ne-poor, buoyant crystals float and transport $^{22}$Ne downward, releasing gravitational energy (Blouin et al. 2021; Bédard, Blouin and Cheng 2024). Yet Blouin et al. (2020) showed that even with an accurate new C/O phase diagram, models still underestimate the magnitude of the crystallization pile-up, and gigayear-level systematic uncertainties in cooling ages persist (Pathak, Blouin and Herwig 2025). Semi-analytic work by Baiko (2022, 2023, 2025) finds miscibility gaps and exsolution (spontaneous separation of a solid solution into two solids) in crystallized C/O and O/Ne mixtures, and shows that ion quantum effects can qualitatively restructure the C/$^{22}$Ne phase diagram — while first-principles quantum treatment by path-integral Monte Carlo (a simulation technique that represents quantum nuclei as discretized closed polymer rings) exists only for the one-component plasma, a single ionic species in a uniform neutralizing electron background (Baiko & Chugunov 2022), not for mixture phase diagrams or fractionation.

Problem Statement

The 2010 determination of the dense carbon–oxygen phase diagram underlying white-dwarf core crystallization posed four open issues. Two — nonequilibrium bias in two-phase coexistence simulations and the missing three-component (C–O–Ne) treatment — have been resolved by later literature, and the binary C/O phase diagram itself is now considered mostly solved by three concordant independent methods. The surviving problem is to complete the microphysical basis of white-dwarf core crystallization to the precision required for cosmochronology, specifically: (1) quantify how the approximations still shared by all current phase-diagram techniques — classical (non-quantum) ions and the assumed body-centered-cubic, randomly substituted solid lattice — affect the liquid–solid phase diagrams and, critically, the fractionation (composition difference between coexisting phases) of crystallizing multicomponent plasmas under white-dwarf conditions, including the miscibility-gap/exsolution behavior of the solid and the phase diagrams of trace species such as $^{22}$Ne whose quantum sensitivity may be qualitative; and (2) achieve quantitative reconciliation between the modern phase diagrams (with trace-species distillation included) and observed white-dwarf cooling signatures — the Gaia crystallization pile-up (the observed excess of white dwarfs at the luminosities where crystallization delays cooling), cluster luminosity functions, and the Q-branch cooling delay (the multi-gigayear extra cooling of massive white dwarfs) — with a careful assessment of observational and model systematic errors. The objective comprises both parts (1) and (2).

References

1. C. J. Horowitz, A. S. Schneider, D. K. Berry, Crystallization of Carbon-Oxygen Mixtures in White Dwarf Stars, Phys. Rev. Lett. 104, 231101 (2010). doi:10.1103/PhysRevLett.104.231101, https://arxiv.org/abs/1005.2441
2. D. E. Winget, S. O. Kepler, F. Campos, M. H. Montgomery, L. Girardi, P. Bergeron, K. Williams, The Physics of Crystallization from Globular Cluster White Dwarf Stars in NGC 6397, Astrophys. J. 693, L6 (2009). doi:10.1088/0004-637X/693/1/L6, https://doi.org/10.1088/0004-637X/693/1/L6
3. Z. Medin, A. Cumming, Crystallization of classical multicomponent plasmas, Phys. Rev. E 81, 036107 (2010). doi:10.1103/PhysRevE.81.036107, https://doi.org/10.1103/PhysRevE.81.036107
4. A. S. Schneider, C. J. Horowitz, J. Hughto, D. K. Berry, Phase diagram of carbon-oxygen plasma mixtures in white dwarf stars, J. Phys. Conf. Ser. 402, 012026 (2012). doi:10.1088/1742-6596/402/1/012026, https://doi.org/10.1088/1742-6596/402/1/012026
5. L. G. Althaus, E. García-Berro, J. Isern, A. H. Córsico, M. M. Miller Bertolami, New phase diagrams for dense carbon-oxygen mixtures and white dwarf evolution, Astron. Astrophys. 537, A33 (2011). doi:10.1051/0004-6361/201117902, https://doi.org/10.1051/0004-6361/201117902
6. M. E. Caplan, C. J. Horowitz, A. Cumming, Neon Cluster Formation and Phase Separation during White Dwarf Cooling, Astrophys. J. Lett. 902, L44 (2020). doi:10.3847/2041-8213/abbda0, https://arxiv.org/abs/2010.00036
7. P.-E. Tremblay, G. Fontaine, N. P. Gentile Fusillo, B. H. Dunlap, B. T. Gänsicke, M. A. Hollands, J. J. Hermes, T. R. Marsh, E. Cukanovaite, T. Cunningham, Core crystallization and pile-up in the cooling sequence of evolving white dwarfs, Nature 565, 202–205 (2019). doi:10.1038/s41586-018-0791-x, https://doi.org/10.1038/s41586-018-0791-x
8. S. Cheng, J. D. Cummings, B. Ménard, A Cooling Anomaly of High-mass White Dwarfs, Astrophys. J. 886, 100 (2019). doi:10.3847/1538-4357/ab4989, https://doi.org/10.3847/1538-4357/ab4989
9. L. R. Bedin, M. Salaris, J. Anderson, M. Libralato, D. Apai, D. Nardiello, R. M. Rich, A. Bellini, A. Dieball, P. Bergeron, A. J. Burgasser, A. P. Milone, A. F. Marino, The HST large programme on NGC 6752 – III. Detection of the peak of the white dwarf luminosity function, Mon. Not. R. Astron. Soc. 488, 3857–3865 (2019). doi:10.1093/mnras/stz1968, https://doi.org/10.1093/mnras/stz1968
10. S. Blouin, J. Daligault, D. Saumon, A. Bédard, P. Brassard, Toward precision cosmochronology: a new C/O phase diagram for white dwarfs, Astron. Astrophys. 640, L11 (2020). doi:10.1051/0004-6361/202038879, https://arxiv.org/abs/2007.13669
11. S. Blouin, J. Daligault, Direct evaluation of the phase diagrams of dense multicomponent plasmas by integration of the Clapeyron equations, Phys. Rev. E 103, 043204 (2021). doi:10.1103/PhysRevE.103.043204, https://arxiv.org/abs/2104.00599
12. S. Blouin, J. Daligault, D. Saumon, $^{22}$Ne Phase Separation as a Solution to the Ultramassive White Dwarf Cooling Anomaly, Astrophys. J. Lett. 911, L5 (2021). doi:10.3847/2041-8213/abf14b, https://arxiv.org/abs/2103.12892
13. S. Blouin, J. Daligault, Phase Separation in Ultramassive White Dwarfs, Astrophys. J. 919, 87 (2021). doi:10.3847/1538-4357/ac1513, https://arxiv.org/abs/2107.07094
14. D. Saumon, S. Blouin, P.-E. Tremblay, Current challenges in the physics of white dwarf stars, Phys. Rep. 988, 1–63 (2022). doi:10.1016/j.physrep.2022.09.001, https://arxiv.org/abs/2209.02846
15. A. S. Jermyn, J. Schwab, E. Bauer, F. X. Timmes, A. Y. Potekhin, Skye: A Differentiable Equation of State, Astrophys. J. 913, 72 (2021). doi:10.3847/1538-4357/abf48e, https://arxiv.org/abs/2104.00691
16. D. A. Baiko, Phase diagrams of binary ionic mixtures and white dwarf cooling, Mon. Not. R. Astron. Soc. 517, 3962–3973 (2022). doi:10.1093/mnras/stac2693, https://arxiv.org/abs/2209.11563
17. D. A. Baiko, A. I. Chugunov, Ab initio thermodynamics of one-component plasma for astrophysics of white dwarfs and neutron stars, Mon. Not. R. Astron. Soc. 510, 2628–2643 (2022). doi:10.1093/mnras/stab3613, https://arxiv.org/abs/2112.04822
18. D. A. Baiko, On excess entropy and latent heat in crystallizing white dwarfs, Mon. Not. R. Astron. Soc. Lett. 522, L26–L30 (2023). doi:10.1093/mnrasl/slad039, https://arxiv.org/abs/2303.11311
19. M. E. Caplan, S. Blouin, I. F. Freeman, Distillation of $^{56}$Fe in Ultramassive O–Ne White Dwarfs, Astrophys. J. 946, 78 (2023). doi:10.3847/1538-4357/acbfaa, https://arxiv.org/abs/2303.03409
20. M. E. Camisassa, D. A. Baiko, S. Torres, A. Rebassa-Mansergas, Exsolution process in white dwarf stars, Astron. Astrophys. 683, A101 (2024). doi:10.1051/0004-6361/202348344, https://arxiv.org/abs/2312.00604
21. A. Bédard, S. Blouin, S. Cheng, Buoyant crystals halt the cooling of white dwarf stars, Nature 627, 286–288 (2024). doi:10.1038/s41586-024-07102-y, https://arxiv.org/abs/2409.04419
22. M. Castro-Tapia, A. Cumming, Three-component Phase Separation for Ultramassive White Dwarf Models, Astrophys. J. 991, 64 (2025). doi:10.3847/1538-4357/adf745, https://arxiv.org/abs/2506.16529
23. D. A. Baiko, Short-range ordering and equilibrium structure of binary crystal mixtures of atomic nuclei in white dwarf cores, Phys. Rev. E 111, 045211 (2025). doi:10.1103/PhysRevE.111.045211, https://arxiv.org/abs/2504.11545
24. P. Pathak, S. Blouin, F. Herwig, Quantifying Systematic Uncertainties in White Dwarf Cooling Age Determinations, Astrophys. J. 983, 158 (2025). doi:10.3847/1538-4357/adc10f, https://arxiv.org/abs/2410.14014