Search for: All records

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. 

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. 

Abstract Stellar evolution models calculate convective boundaries using either the Schwarzschild or Ledoux criterion, but confusion remains regarding which criterion to use. Here we present a 3D hydrodynamical simulation of a convection zone and adjacent radiative zone, including both thermal and compositional buoyancy forces. As expected, regions that are unstable according to the Ledoux criterion are convective. Initially, the radiative zone adjacent to the convection zone is Schwarzschild unstable but Ledoux stable due to a composition gradient. Over many convective overturn timescales, the convection zone grows via entrainment. The convection zone saturates at the size originally predicted by the Schwarzschild criterion, although in this final state the Schwarzschild and Ledoux criteria agree. Therefore, the Schwarzschild criterion should be used to determine the size of stellar convection zones, except possibly during short-lived evolutionary stages in which entrainment persists.