This simulation renders a small cell of a white dwarf's carbon-oxygen core as it crystallizes. Because the Coulomb coupling parameter Γ scales as Z5/3, the higher-charge oxygen ions freeze into the lattice before carbon does at the same temperature. With phase separation switched on, that preference is modeled explicitly: the rising crystallization front captures oxygen ions first, leaving a carbon-enriched liquid buoyant above it — the denser oxygen-rich solid sinking under gravity releases extra energy on top of ordinary latent heat, which is the physical mechanism proposed to explain the pile-up of white dwarfs Gaia observes at intermediate luminosities. Drag the core-temperature slider down (or let Auto-cool do it) and watch the front rise, the readouts track live.
The Coulomb coupling of an ion species scales as Γ_i = Γ_e·Z_i^(5/3). Oxygen has Z=8 versus carbon's Z=6, so at any given temperature oxygen's Γ is roughly 1.4-1.6× higher — it crosses the Γ≈175 crystallization threshold first, while the carbon in the same volume is still a liquid.
Sinking oxygen-rich crystal releases gravitational potential energy in addition to the latent heat of freezing. Together they act as an extra heat source deep in the star, temporarily slowing the drop in luminosity — enough to measurably delay cooling by roughly a billion years for typical white dwarfs, which shows up as an excess pile-up of stars at a particular luminosity in large surveys like Gaia.