Local mixing length theory with compositional effects:\ First application to asymptotic giant branch evolution

TL;DR

This paper introduces an extended mixing length theory, MLT#, that incorporates compositional effects to improve modeling of AGB star evolution, showing significant differences in chemical profiles and pulsation properties compared to standard models.

Contribution

The paper develops and applies MLT# with the Ledoux criterion, providing a simpler alternative to GNA theory for chemically driven convection in stellar evolution.

Findings

MLT# closely reproduces GNA results

Chemical profiles differ significantly between MLT and MLT#

Chemical stratification impacts pulsation spectra of GW Vir stars

Abstract

During the evolution of stars on the asymptotic giant branch (AGB), thermal pulses lead to the formation of strongly stratified layers in the outer regions of the CO core, which might lead to inversions in the chemical gradient. Such inversions would produce instabilities beyond the ones predicted by the Schwarzschild criterion and the standard use of mixing length theory (MLT). We used a set of MLT equations that consider the impact of the background chemical gradients. This extension of MLT is referred to in this work as MLT$\sharp$, to make a distinction between both prescriptions. We applied MLT$\sharp$ in tandem with the more general Ledoux instability criterion. We computed the evolution in the AGB phase and compared the chemical profiles resulting from MLT, MLT$\sharp$ and the double diffusive GNA theory. We continued the evolution through a post-AGB thermal pulse and performed a pulsational analysis of the resultant GW Vir models to asses $g$-mode pulsation periods. Finally, we tested our results with pulsation properties of known GW Vir stars derived from recent observations. We find that the much simpler MLT$\sharp$ set of equations closely reproduces the results from the GNA theory. As such, MLT$\sharp$ offers a simple way to include chemically driven convection in stellar evolution computations. Stellar evolution simulations show that Rayleigh-Taylor and thermohaline instabilities can play an important role during the TP-AGB. We obtained significantly different chemical profiles using a standard MLT approach compared to those resulting from our MLT$\sharp$ and GNA computations. Our adiabatic pulsational analysis shows that these differences in the chemical stratification leave clear mode-trapping signatures in the pulsation spectrum of the GW Vir models.