skip to main content


Title: Multi-gigayear White Dwarf Cooling Delays from Clustering-enhanced Gravitational Sedimentation
Abstract Cooling white dwarfs (WDs) can yield accurate ages when theoretical cooling models fully account for the physics of the dense plasma of WD interiors. We use MESA to investigate cooling models for a set of massive and ultramassive WDs (0.9–1.3 ) for which previous models have failed to match kinematic age indicators based on Gaia DR2. We find that the WDs in this population can be explained as C/O cores experiencing unexpectedly rapid 22 Ne sedimentation in the strongly liquid interior just prior to crystallization. We propose that this rapid sedimentation is due to the formation of solid clusters of 22 Ne in the liquid C/O background plasma. We show that these heavier solid clusters sink faster than individual 22 Ne ions and enhance the sedimentation heating rate enough to dramatically slow WD cooling. MESA models including our prescription for cluster formation and sedimentation experience cooling delays of ≈4 Gyr on the WD Q branch, alleviating tension between cooling ages and kinematic ages. This same model then predicts cooling delays coinciding with crystallization of 6 Gyr or more in lower-mass WDs (0.6–0.8 ). Such delays are compatible with, and perhaps required by, observations of WD populations in the local 100 pc WD sample and the open cluster NGC 6791. These results motivate new investigations of the physics of strongly coupled C/O/Ne plasma mixtures in the strongly liquid state near crystallization and tests through comparisons with observed WD cooling.  more » « less
Award ID(s):
1663688
NSF-PAR ID:
10414681
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
902
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
93
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We enhance the treatment of crystallization for models of white dwarfs (WDs) in the stellar evolution software Modules for Experiments in Stellar Astrophysics (MESA) by implementing carbon–oxygen (C/O) phase separation. The phase separation process during crystallization leads to transport of oxygen toward the centers of WDs, resulting in a more compact structure that liberates gravitational energy as additional heating that modestly slows WD cooling timescales. We quantify this cooling delay in MESA C/O WD models over the mass range 0.5–1.0M, finding delays of 0.5–0.8 Gyr for typical C/O interior profiles. MESA WD cooling timescales including this effect are generally comparable to other WD evolution models that make similar assumptions about input physics. When considering phase separation alongside22Ne sedimentation, however, we find that both MESA and BaSTI WD cooling models predict a more modest sedimentation delay than the latest LPCODE models, and this may therefore require a reevaluation of previously proposed solutions to some WD cooling anomalies that were based on LPCODE models of22Ne sedimentation. Our implementation of C/O phase separation in the open-source stellar evolution software MESA provides an important tool for building realistic grids of WD cooling models, as well as a framework for expanding on our implementation to explore additional physical processes related to phase transitions and associated fluid motions in WD interiors.

     
    more » « less
  2. Abstract

    Gaia's exquisite parallax measurements allowed for the discovery and characterization of the Q branch in the Hertzsprung–Russell diagram, where massive C/O white dwarfs (WDs) pause their dimming due to energy released during crystallization. Interestingly, the fraction of old stars on the Q branch is significantly higher than in the population of WDs that will become Q branch stars or that were Q branch stars in the past. From this, Cheng et al. inferred that ∼6% of WDs passing through the Q branch experience a much longer cooling delay than that of standard crystallizing WDs. Previous attempts to explain this cooling anomaly have invoked mechanisms involving supersolar initial metallicities. In this paper, we describe a novel scenario in which a standard composition WD merges with a subgiant star. The evolution of the resulting merger remnant leads to the creation of a large amount of26Mg, which, along with the existing22Ne, undergoes a distillation process that can release enough energy to explain the Q branch cooling problem without the need for atypical initial abundances. The anomalously high number of old stars on the Q branch may thus be evidence that mass transfer from subgiants to WDs leads to unstable mergers.

     
    more » « less
  3. Abstract

    When white dwarfs freeze, the plasma mixtures inside them undergo separation processes that can produce radical changes in the composition profile of the star. The abundance of neutron-rich elements, such as22Ne or56Fe, determines whether or not the first crystals are more or less dense than the surrounding fluid and thus whether they sink or float. These processes have now been studied for C–O–Ne and C–O–Fe mixtures, finding that distillation and precipitation processes are possible in white dwarfs. In this work, we calculate the phase diagram of more complicated O–Ne–Fe mixtures and make predictions for the internal structure of the separated white dwarf. There are two possible outcomes determined by a complicated interplay between the Ne abundance, the22Ne fraction, and the56Fe abundance. Either Fe distills to form an inner core because the first O–Ne solids are buoyant, or an O–Ne inner core forms and Fe accumulates in the liquid until Fe distillation begins and forms an Fe shell. In the case of an Fe shell, a Rayleigh–Taylor instability may arise and overturn the core. In either case, Fe distillation may only produce a cooling delay of order 0.1 Gyr, as these processes occur early at high white dwarf luminosities. Fe inner cores and shells may be detectable through asteroseismology and could enhance the yield of neutron-rich elements such as55Mn and58Ni in supernovae.

     
    more » « less
  4. Abstract

    The diffusion coefficients of neutron rich nuclei in crystallizing white dwarf (WD) stars are essential microphysics input for modeling the evolution of the composition profile. Recently, molecular dynamics simulations have been used to compute diffusion coefficients for realistic mixtures of C–O and O–Ne WDs with many trace nuclides that could be important sedimentary heat sources such as22Ne,23Na,25Mg, and27Mg. In this brief note, I repeat these simulations but now include56Fe. I find that for the large charge ratios involved in these mixtures the empirical law developed in our earlier work tends to under-predict diffusion coefficients in the moderately coupled regime by 30%–40%. As this formalism is presently implemented in the stellar evolution code MESA, it is important for authors studying mixtures containing heavy nuclides like56Fe to be aware of these systematics. However, the impact on astrophysics is expected to be small.

     
    more » « less
  5. ABSTRACT

    We study environmental quenching using the spatial distribution of current star formation and stellar population ages with the full SAMI Galaxy Survey. By using a star formation concentration index [C-index, defined as log10(r50, H α/r50, cont)], we separate our sample into regular galaxies (C-index ≥−0.2) and galaxies with centrally concentrated star formation (SF-concentrated; C-index <−0.2). Concentrated star formation is a potential indicator of galaxies currently undergoing ‘outside-in’ quenching. Our environments cover ungrouped galaxies, low-mass groups (M200 ≤ 1012.5M⊙), high-mass groups (M200 in the range 1012.5–14 M⊙) and clusters (M200 > 1014M⊙). We find the fraction of SF-concentrated galaxies increases as halo mass increases by 9 ± 2 per cent, 8 ± 3 per cent, 19 ± 4 per cent, and 29 ± 4 per cent for ungrouped galaxies, low-mass groups, high-mass groups, and clusters, respectively. We interpret these results as evidence for ‘outside-in’ quenching in groups and clusters. To investigate the quenching time-scale in SF-concentrated galaxies, we calculate light-weighted age (AgeL) and mass-weighted age (AgeM) using full spectral fitting, as well as the Dn4000 and HδA indices. We assume that the average galaxy age radial profile before entering a group or cluster is similar to ungrouped regular galaxies. At large radius (1–2 Re), SF-concentrated galaxies in high-mass groups have older ages than ungrouped regular galaxies with an age difference of 1.83 ± 0.38 Gyr for AgeL and 1.34 ± 0.56 Gyr for AgeM. This suggests that while ‘outside-in’ quenching can be effective in groups, the process will not quickly quench the entire galaxy. In contrast, the ages at 1–2 Re of cluster SF-concentrated galaxies and ungrouped regular galaxies are consistent (difference of 0.19 ± 0.21 Gyr for AgeL, 0.40 ± 0.61 Gyr for AgeM), suggesting the quenching process must be rapid.

     
    more » « less