Off-centre convective zones in mass accreting stellar models

We report the physical origin of transient off-centre convective zones (oCZs) that arise in mass accreting stellar models. Using detailed MESA simulations of binary evolution, we find that these oCZs are not numerical artefacts, but emerge due to a local increase in density near the retreating edge of the convective core. The density enhancement raises the local opacity, which amplifies the radiative temperature gradient ∇_̊m rad. If this gradient surpasses the Ledoux threshold ∇_̊m L, defined by both thermal and compositional stratification, the region becomes convectively unstable. The resulting oCZs are detached from the convective core and transient: mixing within the oCZ erases the local gradient in mean molecular weight and leaves a sharp ∇_μ discontinuity at the boundary, thus stabilising the adjacent layers. This mechanism naturally explains the presence and evolution of oCZs, as previously reported in massive interacting stars.