Toloza Castillo, Odette Fabiola

ABSTRACT White dwarfs with metal-polluted atmospheres have been studied widely in the context of the accretion of rocky debris from evolved planetary systems. One open question is the geometry of accretion and how material arrives and mixes in the white dwarf surface layers. Using the three-dimensional (3D) radiation hydrodynamics code co5bold, we present the first transport coefficients in degenerate star atmospheres that describe the advection–diffusion of a passive scalar across the surface plane. We couple newly derived horizontal diffusion coefficients with previously published vertical diffusion coefficients to provide theoretical constraints on surface spreading of metals in white dwarfs. Our grid of 3D simulations probes the vast majority of the parameter space of convective white dwarfs, with pure-hydrogen atmospheres in the effective temperature range of 6000–18 000 K and pure-helium atmospheres in the range of 12 000–34 000 K. Our results suggest that warm hydrogen-rich atmospheres (DA; ${\gtrsim} 13\, 000$ K) and helium-rich atmospheres (DB and DBA; ${\gtrsim} 30\, 000$ K) are unable to efficiently spread the accreted metals across their surface, regardless of the time dependence of accretion. This result may be at odds with the current non-detection of surface abundance variations in white dwarfs with debris discs. For cooler hydrogen- and helium-rich atmospheres, we predict a largely homogeneous distribution of metals across the surface within a vertical diffusion time-scale. This is typically less than 0.1 per cent of disc lifetime estimates, a quantity that is revisited in this paper using the overshoot results. These results have relevance for studies of the bulk composition of evolved planetary systems and models of accretion disc physics.

Abstract We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). The new auto_diff module implements automatic differentiation in MESA, an enabling capability that alleviates the need for hard-coded analytic expressions or finite-difference approximations. We significantly enhance the treatment of the growth and decay of convection in MESA with a new model for time-dependent convection, which is particularly important during late-stage nuclear burning in massive stars and electron-degenerate ignition events. We strengthen MESA’s implementation of the equation of state, and we quantify continued improvements to energy accounting and solver accuracy through a discussion of different energy equation features and enhancements. To improve the modeling of stars in MESA, we describe key updates to the treatment of stellar atmospheres, molecular opacities, Compton opacities, conductive opacities, element diffusion coefficients, and nuclear reaction rates. We introduce treatments of starspots, an important consideration for low-mass stars, and modifications for superadiabatic convection in radiation-dominated regions. We describe new approaches for increasing the efficiency of calculating monochromatic opacities and radiative levitation, and for increasing the efficiency of evolving the late stages of massive stars with a new operator-split nuclear burning mode. We close by discussing major updates to MESA’s software infrastructure that enhance source code development and community engagement.