The treatment of mixing in core helium burning models – I. Implications for asteroseismology

Abstract

The detection of mixed oscillation modes offers a unique insight into the internal structure of core helium burning (CHeB) stars. The stellar structure during CHeB is very uncertain because the growth of the convective core, and/or the development of a semiconvection zone, is critically dependent on the treatment of convective boundaries. In this study we calculate a suite of stellar structure models and their non-radial pulsations to investigate why the predicted asymptotic g-mode ℓ = 1 period spacing ΔΠ1 is systematically lower than is inferred from Kepler field stars. We find that only models with large convective cores, such as those calculated with our newly proposed ‘maximal-overshoot’ scheme, can match the average ΔΠ1 reported. However, we also find another possible solution that is related to the method used to determine ΔΠ1: mode trapping can raise the observationally inferred ΔΠ1 well above its true value. Even after accounting for these two proposed resolutions to the discrepancy in average ΔΠ1, models still predict more CHeB stars with low ΔΠ1 ( ≲ 270 s) than are observed. We establish two possible remedies for this: (i) there may be a difficulty in determining ΔΠ1 for early CHeB stars (when ΔΠ1 is lowest) because of the effect that the sharp composition profile at the hydrogen burning shell has on the pulsations, or (ii) the mass of the helium core at the flash is higher than predicted. Our conclusions highlight the need for the reporting of selection effects in asteroseismic population studies in order to safely use this information to constrain stellar evolution theory.