An official website of the United States government
ABSTRACT A large fraction of white dwarfs (WDs) have metal-polluted atmospheres, which are produced by accreting material from remnant planetary systems. The composition of the accreted debris broadly resembles that of rocky Solar system objects. Volatile-enriched debris with compositions similar to long-period comets (LPCs) is rarely observed. We attempt to reconcile this dearth of volatiles with the premise that exo-Oort clouds (XOCs) occur around a large fraction of planet-hosting stars. We estimate the comet accretion rate from an XOC analytically, adapting the ‘loss cone’ theory of LPC delivery in the Solar system. We investigate the dynamical evolution of an XOC during late stellar evolution. Using numerical simulations, we show that 1–30 per cent of XOC objects remain bound after anisotropic stellar mass-loss imparting a WD natal kick of $${\sim}1 \, {\rm km \, s^{-1}}$$. We also characterize the surviving comets’ distribution function. Surviving planets orbiting a WD can prevent the accretion of XOC comets by the star. A planet’s ‘dynamical barrier’ is effective at preventing comet accretion if the energy kick imparted by the planet exceeds the comet’s orbital binding energy. By modifying the loss cone theory, we calculate the amount by which a planet reduces the WD’s accretion rate. We suggest that the scarcity of volatile-enriched debris in polluted WDs is caused by an unseen population of 10–$$100 \, \mathrm{au}$$ scale giant planets acting as barriers to incoming LPCs. Finally, we constrain the amount of volatiles delivered to a planet in the habitable zone of an old, cool WD.
more »
« less
3D Simulations Demonstrate Propagating Thermohaline Convection for Polluted White Dwarfs
Abstract Polluted white dwarfs (WDs) with small surface convection zones deposit significant concentrations of heavy elements to the underlying radiative interior, presumably driving thermohaline convection. Current models of polluted WDs frequently fail to account for this effect, although its inclusion can increase the inferred accretion rate by orders of magnitude. It has been argued that this instability cannot be treated as a continuous mixing process and thus should not be considered in these models. In this work, we study 3D simulations of a thermohaline-unstable layer propagating into an underlying stable region, approximating the polluted WD scenario. We find that although thermohaline convection works to reduce driving gradients somewhat, the front continues to propagate and the system remains unstable. Importantly, the turbulent flux of metals broadly dominates over the diffusive flux in quantitative agreement with existing mixing prescriptions implemented in some stellar evolution models (except slightly below the boundary of the propagating front, where recent prescriptions neglect overshoot-like effects). Thus, our results broadly support polluted WD models that include thermohaline mixing in their estimates of the settling rate.
more »
« less
- Award ID(s):
- 2023499
- PAR ID:
- 10667372
- Publisher / Repository:
- Institute of Physics
- Date Published:
- Journal Name:
- The Astrophysical Journal Letters
- Volume:
- 986
- Issue:
- 1
- ISSN:
- 2041-8205
- Page Range / eLocation ID:
- L10
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Polluted white dwarfs (WDs) offer a unique way to study the bulk compositions of exoplanetary material, but it is not always clear if this material originates from comets, asteroids, moons, or planets. We combineN-body simulations with an analytical model to assess the prevalence of extrasolar moons as WD polluters. Using a sample of observed polluted WDs, we find that the extrapolated parent body masses of the polluters are often more consistent with those of many solar system moons, rather than solar-like asteroids. We provide a framework for estimating the fraction of WDs currently undergoing observable moon accretion based on results from simulated WD planetary and moon systems. Focusing on a three-planet WD system of super-Earth to Neptune-mass bodies, we find that we could expect about one percent of such systems to be currently undergoing moon accretions as opposed to asteroid accretion.more » « less
-
Abstract Observations show an almost ubiquitous presence of extra mixing in low-mass upper giant branch stars. The most commonly invoked explanation for this is thermohaline mixing. One-dimensional stellar evolution models include various prescriptions for thermohaline mixing, but the use of observational data directly to discriminate between thermohaline prescriptions has thus far been limited. Here, we propose a new framework to facilitate direct comparison: using carbon-to-nitrogen measurements from the Sloan Digital Sky Survey-IV APOGEE survey as a probe of mixing and a fluid parameter known as the reduced density ratio from one-dimensional stellar evolution programs, we compare the observed amount of extra mixing on the upper giant branch to predicted trends from three-dimensional fluid dynamics simulations. Using this method, we are able to empirically constrain how mixing efficiency should vary with the reduced density ratio. We find the observed amount of extra mixing is strongly correlated with the reduced density ratio and that trends between reduced density ratio and fundamental stellar parameters are robust across choices for modeling prescription. We show that stars with available mixing data tend to have relatively low density ratios, which should inform the regimes selected for future simulation efforts. Finally, we show that there is increased mixing at low reduced density ratios, which is consistent with current hydrodynamical models of thermohaline mixing. The introduction of this framework sets a new standard for theoretical modeling efforts, as validation for not only the amount of extra mixing, but trends between the degree of extra mixing and fundamental stellar parameters is now possible.more » « less
-
Abstract We seek signatures of the current experimental 12 C α , γ 16 O reaction rate probability distribution function in the pulsation periods of carbon–oxygen white dwarf (WD) models. We find that adiabatic g-modes trapped by the interior carbon-rich layer offer potentially useful signatures of this reaction rate probability distribution function. Probing the carbon-rich region is relevant because it forms during the evolution of low-mass stars under radiative helium-burning conditions, mitigating the impact of convective mixing processes. We make direct quantitative connections between the pulsation periods of the identified trapped g-modes in variable WD models and the current experimental 12 C α , γ 16 O reaction rate probability distribution function. We find an average spread in relative period shifts of Δ P / P ≃ ±2% for the identified trapped g-modes over the ±3 σ uncertainty in the 12 C α , γ 16 O reaction rate probability distribution function—across the effective temperature range of observed DAV and DBV WDs and for different WD masses, helium shell masses, and hydrogen shell masses. The g-mode pulsation periods of observed WDs are typically given to six to seven significant figures of precision. This suggests that an astrophysical constraint on the 12 C α , γ 16 O reaction rate could, in principle, be extractable from the period spectrum of observed variable WDs.more » « less
-
Abstract Fingering convection (also known as thermohaline convection) is a process that drives the vertical transport of chemical elements in regions of stellar radiative zones where the mean molecular weight increases with radius. Recently, Harrington & Garaud used three-dimensional direct numerical simulations (DNS) to show that a vertical magnetic field can dramatically enhance the rate of chemical mixing by fingering convection. Furthermore, they proposed a so-called “parasitic saturation” theory to model this process. Here, we test their model over a broad range of parameter space, using a suite of DNS of magnetized fingering convection, varying the magnetic Prandtl number, magnetic field strength, and composition gradient. We find that the rate of chemical mixing measured in the simulations is not always predicted accurately by their existing model, in particular when the magnetic diffusivity is large. We then present an extension of the Harrington & Garaud model which resolves this issue. When applied to stellar parameters, it recovers the results of Harrington & Garaud except in the limit where fingering convection becomes marginally stable, where the new model is preferred. We discuss the implications of our findings for stellar structure and evolution.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.3847/2041-8213/addbd5
