Sven Erik Matzen

Software Architect | Cloud & Security Expert | AI-enabled Solutions

Crystals in the Sky: White Dwarfs and the Clock That Dates Our Galaxy

🎧 Listen to this article

Astrophysics · 2026-11-05

EU label: fully AI-generated content Fully AI-generated article (no prior review).

The Hook: A clock whose error only Gaia could see

Imagine you have to determine the age of a building, and the only measuring instrument you have is a cup of coffee. You know how hot coffee is poured. You measure how warm it still is. From the cooling curve you work backward to how long the cup has been sitting there. This works — remarkably well, in fact — as long as you truly have one thing under control: the cooling law.

This is precisely the method astrophysics uses to determine the age of our galaxy. The cups are white dwarfs: burnt-out stellar cores that no longer run nuclear fusion and simply radiate away for billions of years. If you know how fast a white dwarf cools, you can read off from its present brightness how long it has been dead — and therefore when the first generation of stars in that region of the Milky Way died. The technique is called white dwarf cosmochronology.

And here is the point that carries this article. The cooling models had contained the decisive ingredient — a phase transition inside the star — as a prediction for decades. What was missing was the resolution to go and look. Only when a satellite named Gaia measured the distances to hundreds of thousands of these objects with previously unthinkable precision did two things become clear: the predicted effect is real — and for a subpopulation of high-mass white dwarfs, the models do not reproduce the observed delay. The discrepancy is not measured in millions of years but in billions.

White dwarfs crystallize: they transition from liquid to solid, like water turning to ice, and in doing so they release heat that slows the clock. This article tells three things: how the clock works, what Gaia actually showed in 2019, and what follows methodologically — even for someone who will never professionally lay eyes on a white dwarf, but who does work with models whose errors only become visible once measurement resolution improves.


Part 1: What a white dwarf actually is

The end of an ordinary star

A star like the Sun lives by fusing light nuclei into heavier ones in its interior. The radiation pressure from these reactions holds gravity in check. Eventually the fuel is spent — and then gravity wins. The outer envelopes are shed, and what remains is the naked, extremely compressed core.

NASA describes this remnant as a "roughly Earth-sized stellar cinder" that gradually cools over billions of years. And it states explicitly that this is also the fate of our Sun.

The numbers are absurd. According to NASA, an Earth-sized white dwarf has an average density of about 1 × 10⁹ kg/m³, while Earth itself reaches a mean of only 5.4 × 10³ kg/m³ — making the white dwarf roughly 200,000 times as dense. A typical white dwarf, per NASA, is about half as massive as the Sun yet only slightly bigger than Earth.

For the center, Tremblay and colleagues (2019) give values around 10⁷ grams per cubic centimeter, that is 10¹⁰ kg/m³ — ten times more than NASA's average, because the core is naturally denser than the object as a whole. For a sense of scale, a simple conversion: a teaspoon (about 5 cm³) of this core material would amount to 5 × 10⁷ grams, roughly 50 tonnes — more than a fully loaded 40-tonne semi-trailer truck.

Why the thing does not collapse further

What holds a white dwarf together at all — without fusion, without radiation pressure? The answer is pure quantum mechanics: electron degeneracy pressure. The Pauli principle forbids two electrons from occupying the same quantum state. Compress matter far enough and the electrons must retreat into ever higher momentum states. The resulting pressure does not depend on temperature — it exists even at absolute zero. That is why a white dwarf can become arbitrarily cold without collapsing.

This pressure has an upper limit, though. Above a certain mass it can no longer withstand gravity. NASA puts it this way: Subrahmanyan Chandrasekhar found this limit to be 1.4 times the mass of the Sun — known as the Chandrasekhar limit.

Chandrasekhar received the 1983 Nobel Prize in Physics, shared equally with William Alfred Fowler. The Nobel committee's citation for Chandrasekhar's half reads verbatim: "for his theoretical studies of the physical processes of importance to the structure and evolution of the stars." Notably, the official citation names no single work and does not name the mass limit by name — it honors the body of work in its breadth. This is a small but instructive detail for anyone who checks sources: the popular shorthand "Nobel Prize for the Chandrasekhar limit" is not what the citation says.

The composition: carbon and oxygen

What matters for everything that follows is what the core is made of. A star of intermediate mass fuses helium into carbon and on into oxygen, but does not get as far as iron. So what typically remains is a carbon-oxygen core (C/O), surrounded by thin helium and hydrogen layers. Added to this are traces of heavier elements — including, and this becomes the pivot of the story, the neon isotope ²²Ne.

The electrons in this core are degenerate. The ions — the carbon and oxygen nuclei — are not. They form a strongly coupled Coulomb plasma: a swarm of charged particles whose electrostatic interaction dominates thermal motion. And such a plasma can freeze.


Part 2: The clock — cooling as age measurement

The principle

Because a white dwarf no longer has an energy source, its evolution is essentially a thermodynamic cooling problem: stored heat is radiated away, luminosity drops, color reddens. If you know the cooling curve for a given mass and composition, then present-day luminosity is a direct measure of the cooling age.

Add the lifetime of the progenitor star and you get the total age. And if you look for the oldest, faintest white dwarfs in a population, you get a lower bound on that population's age. The review by Fontaine, Brassard & Bergeron (2001) in the Publications of the Astronomical Society of the Pacific, titled "The Potential of White Dwarf Cosmochronology" (PASP 113, 409), states exactly this logic: the oldest white dwarfs in the solar neighborhood — the remnants of the first generation of intermediate-mass disk stars — are still visible and can be used, in conjunction with cooling theory, to estimate the age of the disk. The article additionally points to a possible very old population of white dwarfs in the galactic halo that could yield an independent halo age.

The decisive phrase there is the fine print: "in conjunction with cooling theory". The method is only as good as the cooling model. It does not measure age, it measures brightness — and translates that with a model.

What the method delivers

How robust is this in practice? Kilic, Munn, Harris, von Hippel, Liebert, Williams, Jeffery and DeGennaro published an analysis of nearby white dwarfs in the Astrophysical Journal in 2017 (volume 837, page 162) titled "The Ages of the Thin Disk, Thick Disk, and the Halo from Nearby White Dwarfs". The results, as given in the abstract:

Population Local 40 pc sample Deep proper-motion catalog
Thin disk 6.8–7.0 Gyr 7.4–8.2 Gyr
Thick disk 8.7 ± 0.1 Gyr 9.5–9.9 Gyr
Inner galactic halo 12.5 (+1.4 / −3.4) Gyr –

Two further figures from the same work: the onset of star formation in the thin disk lagged that of the thick disk by 1.6 (+0.3 / −0.4) Gyr, and the ratio of thick-disk to thin-disk white dwarfs is "roughly 20 %" in the local 40 pc sample — a statement about that sample, not about the galaxy as a whole.

Two things stand out. First: the method works — it yields a plausible, layered chronology of the Milky Way. Second: the two samples within the same paper give different ages for the same thin disk — 6.8–7.0 versus 7.4–8.2 Gyr, two intervals that do not overlap. The paper does not directly explain this difference; it discusses as relevant factors the proper-motion limit of the 40 pc sample, the tangential-velocity cut of the deep survey, and the assumption about scale height (constant versus age-dependent), and explicitly attributes an earlier discrepancy to the assumption of a constant scale height.

This is not a scandal but science operating normally. It does show, however, where the method hurts — not in the photometry, but in the model and in the sample definition.


Part 3: The phase transition — when a star freezes

The prediction

For decades it had been theoretically clear: if a strongly coupled Coulomb plasma cools far enough, it must undergo a first-order phase transition. The ions arrange themselves into a crystal lattice. And like every first-order phase transition, this one releases latent heat — the same physics that keeps a glass of ice water at 0 °C while the ice melts, rather than immediately warming further.

For a white dwarf this means: when crystallization sets in, additional energy flows in that must be radiated away. The star lingers longer at a given luminosity. Cooling is delayed.

In the abstract by Tremblay and colleagues (2019), this prediction is summarized as follows: it has been predicted that a first-order phase transition occurs during cooling, leading to the crystallization of the non-degenerate carbon and oxygen ions in the core, which releases a considerable amount of latent heat and "delays the cooling process by about one billion years".

Important for interpretation: in that sentence, this roughly one billion years is the theoretical prediction, not a value measured by Gaia. In the same breath, the abstract notes that no direct observational evidence of the effect had been reported up to that point.

The effect was expected to become observable as an accumulation — a "pile-up" — of white dwarfs at particular brightnesses and colors: because they linger there, we find more of them there than uniform cooling would predict. Just as more cars stand at a traffic light than on open road.

The evidence (2019)

In 2019, Nature (volume 565, issue 7738, pages 202–205, DOI 10.1038/s41586-018-0791-x) published "Core crystallization and pile-up in the cooling sequence of evolving white dwarfs" by Pier-Emmanuel Tremblay, Gilles Fontaine, Nicola Pietro Gentile Fusillo, Bart H. Dunlap, Boris T. Gänsicke, Mark A. Hollands, J. J. Hermes, Thomas R. Marsh, Elena Cukanovaite and Tim Cunningham.

The result, as stated in the abstract: the authors report the presence of a pile-up in the cooling sequence of evolving white dwarfs within 100 parsecs of the Sun, determined using photometry and parallax data from the Gaia satellite. From modelling they infer that this pile-up arises from the release of latent heat as the cores crystallize — the abstract's wording is "we infer that this pile-up arises from the release of latent heat". So it is a model inference from an observed overdensity, not direct calorimetry.

And then comes the part that goes beyond mere confirmation. The authors state that, in addition to the release of latent heat, they find "strong evidence" that cooling is further slowed by the liberation of gravitational energy from element sedimentation in the crystallizing cores.

These are therefore two separate energy mechanisms. When a C/O mixture crystallizes, both species do not freeze out simultaneously. In the accompanying press coverage, Tremblay describes the process as the oxygen crystallizing first and then sinking to the core — a process similar to sedimentation on a river bed on Earth. This separation shifts mass inward, and in doing so gravitational energy is released, in addition to the latent heat.

In the same coverage, Tremblay places the finding historically: this is the first direct evidence that white dwarfs crystallize, that is, transition from liquid to solid. Fifty years ago it was predicted that, due to crystallization, a pile-up in the number of white dwarfs should be observed at certain luminosities and colors, and only now has that been observed.

Further figures are given there on the temperature scale and the sequence of events: when the core has cooled to about 10 million degrees, enough energy has been released that the fluid begins to solidify, forming a metallic core with a carbon-enriched mantle. All white dwarfs will crystallize at some point in their evolution, though more massive ones go through the process sooner; billions of white dwarfs in our galaxy have already completed the process and are essentially crystal spheres in the sky. The Sun itself will become a crystal white dwarf in about 10 billion years.

That crystalline structures of carbon and oxygen are expected inside many white dwarfs is also documented in NASA material: it describes a hypothesized crust about 50 km thick, at the bottom of which lies a crystalline lattice of carbon and oxygen atoms.

A note of caution, since the phrase "diamonds in the sky" is popular in accessible accounts: I am of the opinion that this image conceals more than it explains. A lattice of carbon and oxygen ions in a degenerate electron background is not a diamond in the mineralogical sense — it is a Coulomb lattice under conditions that have no counterpart on Earth. That is my own assessment and not part of the cited works; the NASA source speaks only of a "crystalline lattice of carbon and oxygen atoms".


Part 4: The Q branch — where the clock really goes wrong

Had the story ended here, it would be a clean success: prediction, measurement, confirmation, model improved. But a second finding came the same year — and it is uncomfortable.

Cheng, Cummings and Ménard published "A Cooling Anomaly of High-mass White Dwarfs" in the Astrophysical Journal in 2019 (volume 886, page 100, DOI 10.3847/1538-4357/ab4989). The subject is an overdensity in the Hertzsprung–Russell diagram that the authors call the "Q branch" and which, per the abstract, is located at the high-mass end of the recently identified crystallization branch.

The method is elegant because it sets two different age indicators, resting on quite different assumptions, against each other:

  1. Photometric isochrone age — the model-based cooling age from color and brightness.
  2. Dynamical age — estimated from transverse velocity. Stars are dynamically "heated" over time by encounters in the galactic disk; their velocity dispersion grows with age. Whoever is fast is statistically old.

Comparing these two quantities is the real point of the paper. (The image of "two independent clocks" is my phrasing, not the authors'; they speak of comparing the dynamical age inferred from transverse velocity with the photometric isochrone age.)

If the cooling model is right, both procedures must agree on average. The authors compared them for more than one thousand high-mass white dwarfs with masses from 1.08 to 1.23 solar masses drawn from Gaia DR2.

They do not agree. Per the abstract, the number density and velocity distribution on the Q branch cannot be explained by the cooling delay of crystallization alone, suggesting the existence of an extra cooling delay. Quantitatively: about 6 % of these high-mass white dwarfs must experience an extra cooling delay of 8 billion years on the Q branch — in addition to crystallization and merger delays. The authors explicitly call this "a challenge for WD cooling models" and note as a possible mechanism that ²²Ne settling in C/O-core white dwarfs could account for this extra cooling delay.

It is worth dwelling on that number for a moment. Eight billion years. That is more than half the lifetime of the universe, for a class of objects we use as a chronometer for that very universe. And note the structure of the claim: it is not that all white dwarfs are misdated by 8 Gyr, but that a minority of about 6 % within a narrow mass window of 1.08–1.23 M☉ is. The clock does not simply run slow — a small fraction of the clocks nearly stop for a while. For population statistics that is the more awkward variant, because a systematic offset could be corrected, whereas a rare, extreme subpopulation distorts the distribution.


Part 5: Distillation — a candidate for the solution

Where could billions of years of extra delay come from? The leading candidate is a trace element: ²²Ne.

²²Ne forms in the star's earlier history out of the original nitrogen. What matters is its neutron-rich nature: A > 2Z holds, so it carries more neutrons than its nuclear charge would "require". It is precisely from this that the large release of gravitational energy during phase separation in the crystallizing C/O core follows.

Blouin, Daligault and Saumon presented "²²Ne Phase Separation as a Solution to the Ultramassive White Dwarf Cooling Anomaly" in the Astrophysical Journal Letters in 2021 (volume 911, article L5, DOI 10.3847/2041-8213/abf14b). The proposed mechanism is not simple settling but distillation — and it is beautiful:

During crystallization, crystals form that are depleted in ²²Ne. Depleted here means: less dense than the surrounding liquid. So they float up, reach warmer layers and melt again. What is left behind and displaced downward is ²²Ne-rich liquid. A cycle of rising and melting thereby pumps ²²Ne systematically inward — a distillation process that concentrates heavy material at the center and in doing so releases a considerable amount of gravitational energy.

The figures from that paper (Table 1), exactly as given there:

Case Composition Extra delay Δτ
²²Ne-rich central core, 1.0 M☉ X(²²Ne)=0.035, X(O)=0.50 9.0 Gyr
²²Ne-rich central core, 1.1 M☉ X(²²Ne)=0.035, X(O)=0.50 7.1 Gyr
²²Ne-rich central core, 1.2 M☉ X(²²Ne)=0.035, X(O)=0.50 4.6 Gyr
²²Ne-rich central core, 1.2 M☉ X(²²Ne)=0.035, X(O)=0.60 3.6 Gyr
²²Ne-rich shell, 0.8 M☉ X(²²Ne)=0.014, X(O)=0.60 2.0 Gyr
²²Ne-rich shell, 1.0 M☉ X(²²Ne)=0.014, X(O)=0.60 1.6 Gyr

In summary, for the case of a ²²Ne-rich central core at X(²²Ne)=0.035, the paper gives an extra delay of 4 to 9 billion years, which in the body of the text (Section IV) is described as explaining "most of the ultramassive white dwarf cooling anomaly". In the abstract, the wording is that the mechanism can "largely resolve" the anomaly.

And now the part that tends to vanish from popular summaries: the authors weaken their own claim themselves, and do so repeatedly.

  • Their phrasing is conditional: the mechanism can largely resolve the anomaly if the delayed population consists of white dwarfs with moderately above-average ²²Ne abundances. That is not a statement of fact.
  • They point out that Bauer and colleagues (2020) argue the missing delay can be reduced to about 6 Gyr by changing assumptions in Cheng et al.'s analysis — thick disk age, age–velocity dispersion relation, star formation history.
  • They name two opposing corrections to the tracks used by Cheng et al.: for masses below 1.10 M☉, C/O phase separation is omitted, which underestimates time on the Q branch by up to about 1 Gyr; for masses at or above 1.10 M☉, O/Ne tracks were used, which overestimates by about 1.5 Gyr for objects with C/O cores.
  • They name an open problem with their own model: normal-mass white dwarfs with ²²Ne-rich central cores should show delays of several billion years that are not observed. As explanations they offer that few white dwarfs are that ²²Ne-rich, or that the central oxygen content is considerably higher (X(O) ≈ 0.85).
  • And they close with the statement that detailed cooling sequences and population synthesis simulations will be needed to confirm their conclusions.

It is important not to set the two numbers 8 and 6 against each other incorrectly. The 8 Gyr is the extra delay required by Cheng et al. for about 6 % of objects. The roughly 6 Gyr is not a second, competing measurement of the same quantity but the amount that remains to be explained after assumptions are changed. So the size of the anomaly hangs on the assumptions — and additionally on the track corrections mentioned above of up to ~1 Gyr (below 1.10 M☉) and ~1.5 Gyr (at or above 1.10 M☉).

The honest state of play is therefore: there is an observed anomaly whose amount to be explained is assumption-dependent, and a physically plausible, quantitatively fitting candidate mechanism that its own authors describe as not yet confirmed. That is not a deficiency — that is what good science looks like in its open state.


Part 6: What this means methodologically

Here the article leaves astrophysics and becomes general. Because the pattern of this story is not astronomical, it is epistemic — and it confronts anyone who works with models.

The pattern

Phase What happens Analogy in everyday technical work
1. Model established Cooling theory yields plausible ages; crystallization is included as a prediction but unobserved A capacity model with a term nobody ever measured
2. Resolution improves Gaia delivers parallaxes for hundreds of thousands of objects New monitoring, finer metrics, more instances
3. Structure becomes visible Pile-up, Q branch — patterns previously buried in noise The latency distribution suddenly shows two humps
4. Explanation A confirmed Crystallization: predicted, then evidenced as a pile-up The suspected cause really does apply
5. Residual error remains Extra delay for ~6 % of objects in 1.08–1.23 M☉ After the fix, a remainder is left that does not fit
6. New mechanism ²²Ne distillation as a candidate, described by its authors as unconfirmed A second, independent cause, not yet verified

The most valuable point is row 5. After a successful fix, the temptation to write off the residual error as noise is enormous. Cheng and colleagues did not do that — and precisely because of this, a second physical process became visible as a candidate at all.

Two measurement paths are better than one

The methodologically strongest element of the Q branch work is not a better model but a second measurement path with different assumptions. The photometric age comes from cooling theory, the dynamical age from kinematics. Neither of them is "true" — but their difference is diagnostic. A model error that creeps identically into both measurements is invisible; one that affects only one of the two becomes visible immediately.

Anyone who builds systems knows this principle: a measured value without an independent cross-check is an assertion. Two differently derived values whose difference you observe are a sensor for model error. The same idea lies behind clocks that carry their own uncertainty rather than delivering a number without error bars — see Clocks That Know Their Own Uncertainty: Google Spanner, TrueTime, and Mastering Time in the Cloud.

The distortion caused by the minority

Especially instructive is that not all objects are affected, but about 6 % in a narrow mass window. A uniform error is calibratable. A rare subpopulation with extreme behavior is not — it shifts means, inflates dispersions, and produces a tail in age distributions that is easily misread as a "particularly old population". Anyone working with percentiles in the cloud knows the same problem from the other side: the mean is blind to exactly the minority that causes the trouble — see The Long Tail of Latency: Tail Latency and Why Averages Lie in the Cloud.


Part 7: What remains open

Honesty about limits belongs to the topic, so explicitly:

  1. The amount to be explained is assumption-dependent. Cheng et al. (2019) require 8 Gyr for about 6 % of objects; per the argument of Bauer et al. (2020) cited by Blouin et al., about 6 Gyr remains if assumptions about the thick disk, the age–velocity dispersion relation and the star formation history are changed.
  2. The mechanism is a candidate, not a consensus. The authors themselves call for detailed cooling sequences and population synthesis to confirm it.
  3. The composition assumptions are uncertain. The delay depends sensitively on X(²²Ne) and X(O), as the table in Part 5 shows.
  4. A counter-observation is unresolved. Normal-mass white dwarfs with ²²Ne-rich central cores should show multi-billion-year delays that are not observed.
  5. The disk ages are model-dependent. The two samples in Kilic et al. (2017) give 6.8–7.0 and 7.4–8.2 Gyr respectively for the same thin disk.

I am of the opinion that this finding does not devalue white dwarf cosmochronology but rather makes it grow up: a method that knows and quantifies its own systematics is more trustworthy than one that appears to work perfectly because its resolution is too low to reveal problems. That assessment is my evaluation, not a statement of the cited works.


Fact Check

Claims checked: 41 (18 hard facts, 6 quotations, 5 source attributions, 12 synthesis claims) plus a check for internal consistency. Corrections or weakenings: 12. No claim was entirely FALSE; none remained UNVERIFIABLE.

The three most important corrections:

  1. The hook in the draft claimed the cooling law had been "systematically wrong for decades". No source supports this: crystallization had been part of the models as a prediction for decades, and Gaia showed for the first time that the models do not reproduce the size of the effect for a subpopulation of high-mass white dwarfs — weakened accordingly.
  2. The explanation that ²²Ne is heavier than its surroundings because buoyancy in a degenerate plasma is governed by the mass-to-charge ratio was wrongly attributed to Blouin et al. The paper attributes the large energy release to the neutron-rich nature of ²²Ne (A > 2Z) and describes the phase separation in terms of density; the text was corrected accordingly.
  3. The juxtaposition "8 versus 6 Gyr" was skewed: the 6 Gyr is not a competing measurement but the amount still remaining to be explained after assumptions are changed. The section was rewritten, and "The Central Takeaway" as well as Part 7 were adjusted to match.

Further weakenings concerned: the 20 % ratio (applies only to the local 40 pc sample), "contradict each other" → "give different ages" for Kilic et al. along with removal of the article's own selection-effect explanation, the diamond comparison (marked as the author's own assessment), the truck comparison (50 t exceeds a 40-tonne semi-trailer), the attribution of the two Blouin quotations to their locations (abstract vs. Section IV), the phrase "two independent clocks" (the article's image, not the source's), the referent of the fifty-year-old prediction (it concerns the pile-up in number, not the size of the delay), and the density figures (NASA's average versus Tremblay et al.'s central value, to avoid an apparent contradiction).

Automated AI verification against primary sources; no human expert review.


The Central Takeaway

Three levels:

  1. Physically: white dwarfs crystallize — Tremblay et al. (2019) report a pile-up within 100 pc using Gaia data and infer from modelling that latent heat is released, with strong evidence for an additional contribution from gravitational energy via element sedimentation.
  2. Methodologically: the cooling clock with which we date Milky Way populations has a known, not yet conclusively explained systematic: per Cheng et al. (2019), about 6 % of white dwarfs with 1.08–1.23 M☉ would have to experience an extra delay of 8 Gyr; the amount to be explained hangs on the assumptions and reduces to roughly 6 Gyr per Bauer et al. (2020). The leading candidate — ²²Ne distillation — is marked as unconfirmed by its own authors.
  3. Epistemically: progress came not from a better model but from a second measurement path with different assumptions (dynamical versus photometric age) and from not writing off the residual error as noise after a successful fix.

A concrete prompt to act on this week: Take a number you regularly rely on in your work — a capacity estimate, a cost forecast, a risk assessment, a runtime assumption — and ask it two questions. First: which model translates measurement into claim here, and when was that model last checked against reality? Second: is there a second, differently derived way to estimate the same quantity — and have I ever looked at the difference? If the answer to the second question is "no", you do not have a measurement; you have an assertion with decimal places.

Reflection question: Astronomy spent decades computing with a clock whose predicted error had never been verified — not out of negligence, but because the resolution to see it was missing. Which assumption in your own models stands unrefuted today only because you are not yet looking closely enough?


Cross-References in the Vault


Sources and Further Reading


Created as part of the daily learning workflow. Area of interest: Astrophysics. Estimated reading time: ~30 minutes.

← All articles