Transient Off-centre Convective Zones (oCZs)

• Transient off-centre convective zones (oCZs) are thin, detached convective shells forming outside the main convective region due to local density, opacity, and radiative gradient enhancements.
• They are short-lived phenomena observed in mass-accreting stars and massive stellar envelopes, impacting mixing processes, wave transmission, and asteroseismic signatures.
• Key formation mechanisms include mass-accretion induced opacity bumps and convective penetration, underscoring challenges in modeling nonlocal boundary physics in stellar interiors.

Searching arXiv for the cited papers and closely related work on off-centre convective zones, convective penetration, and wave transmission in layered radiative/convective structures. Transient off-centre convective zones (oCZs) are thin, detached convective shells that arise away from the stellar centre, embedded within otherwise radiative regions and separated from the main convective core or envelope by Ledoux- or Schwarzschild-stable layers. In the most direct recent formulation, oCZs in mass-accreting stellar models are “thin, short-lived convective shells that form just outside (at larger mass coordinate than) the main hydrogen-burning convective core” and are produced by a local increase in density, opacity, and radiative temperature gradient near the retreating core boundary (Miszuda, 18 Aug 2025). More broadly, related work shows that off-centre convective structures also arise from opacity bumps in hot massive stars, from convective penetration beyond formal stability boundaries, and from alternating radiative–convective stratifications that affect wave transmission, mixing, and mode coherence (Cantiello et al., 2010, Anders et al., 2021, Cai et al., 2021). The term therefore denotes both a structural phenomenon in stellar evolution calculations and a dynamical interface problem involving overshooting, penetrative convection, gravity-wave coupling, and composition-gradient feedback.

1. Definition and structural identity

In mass-accreting stellar models, oCZs appear as “thin convective shells located above the convective core in mass coordinate,” detached from it by “a radiative layer that is locally stable according to the Ledoux criterion” (Miszuda, 18 Aug 2025). They are chemically almost homogeneous within each shell, while very narrow intershell regions retain steep mean-molecular-weight gradients. This detached geometry is central: the core remains deeper and continuously present during core hydrogen burning, whereas the oCZs are shells at larger radius or mass coordinate that appear only during specific structural readjustments (Miszuda, 18 Aug 2025).

The defining instability condition in the mass-accretion context is the Ledoux criterion,

$$\nabla_{\rm rad} > \nabla_{\rm L}, \qquad \nabla_{\rm L} = \nabla_{\rm ad} + \frac{\varphi}{\delta}\,\nabla_{\mu},$$

with $\nabla_{\mu} = d\ln\mu/d\ln P$ supplying the stabilizing compositional term (Miszuda, 18 Aug 2025). In this framework, oCZs form where the radiative temperature gradient temporarily exceeds the Ledoux threshold inside a composition-stratified region, but not immediately at the convective-core boundary because the innermost $\mu$-gradient is too strong there (Miszuda, 18 Aug 2025).

A broader stellar-structure usage is supported by related studies. The iron convection zone (FeCZ) in hot massive OB stars is explicitly described as an “off-centre, subsurface convection zone” produced by the iron opacity bump and occupying a thin shell in the outer envelope (Cantiello et al., 2010). Convective penetration studies likewise show that the dynamically mixed region can extend beyond the Schwarzschild boundary into a formally stable layer, creating an effectively convective region outside the local instability boundary (Anders et al., 2021). This suggests that “oCZ” denotes not a single microphysical mechanism, but a class of detached convective regions arising from local opacity physics, boundary mixing, or structural readjustment.

2. Formation mechanisms in stellar interiors

The most explicit physical origin presently identified is the mass-accretion mechanism. During and after Roche-lobe overflow, the convective core of the accreting star expands, rejuvenates, and later partly retreats. Structural readjustment near the retreating core edge creates a narrow region where density and hydrogen abundance increase locally relative to adjacent layers. Because the radiative gradient scales as

$\nabla_{\rm rad} \propto \frac{\kappa L P}{m T^4},$

the density and hydrogen enhancement raises the Rosseland mean opacity $\kappa$, which in turn produces a local peak in $\nabla_{\rm rad}$ (Miszuda, 18 Aug 2025). If that peak crosses the Ledoux threshold, a detached convective shell forms (Miszuda, 18 Aug 2025).

The sequence reported for mass gainers is specific: mass accretion deposits fresh H-rich material; the core expands and partially mixes with surrounding layers; a strong $\mu$-discontinuity is established at the core edge; a narrow high-density layer forms just outside the core; this generates an opacity peak and radiative-gradient spike; convection starts in a shell where the $\mu$-gradient is present but still shallow enough that $\nabla_{\rm rad}$ can exceed $\nabla_{\rm L}$ (Miszuda, 18 Aug 2025). Convective mixing then homogenizes the shell interior while sharpening $\nabla_{\mu} = d\ln\mu/d\ln P$0 at its edges, which stabilizes adjacent layers and can segment broader unstable regions into multiple thinner oCZs (Miszuda, 18 Aug 2025).

Other astrophysical realizations differ in trigger but preserve the off-centre geometry. In hot massive stars, the FeCZ is caused by “a peak in the opacity due to iron recombination” near the surface, so that a thin shell in the outer envelope becomes convective while deeper and shallower layers remain radiative (Cantiello et al., 2010). In convective-penetration models, the local radiative conductivity profile near the Schwarzschild boundary controls whether a penetration zone develops above or below the formal boundary, with the penetration parameter $\nabla_{\mu} = d\ln\mu/d\ln P$1 governing how far the adiabatic stratification extends into the stable side (Anders et al., 2021). This suggests that oCZs emerge whenever a local flux/opacity/composition configuration produces a confined region in which convective energy transport or convective mixing remains dynamically favored despite surrounding stability.

3. Detachment, transience, and boundary physics

The detached character of oCZs is controlled by composition gradients. In the mass-accretion models, the layer immediately above the convective core is stabilized by a sufficiently strong positive $\nabla_{\mu} = d\ln\mu/d\ln P$2, so even an elevated $\nabla_{\mu} = d\ln\mu/d\ln P$3 does not trigger convection there (Miszuda, 18 Aug 2025). Slightly farther out, where the $\nabla_{\mu} = d\ln\mu/d\ln P$4-gradient is shallower but the opacity peak is strongest, the Ledoux condition is violated and a separate shell turns convective (Miszuda, 18 Aug 2025). This produces the characteristic sequence: convective core, radiative $\nabla_{\mu} = d\ln\mu/d\ln P$5-gradient layer, detached convective shell (Miszuda, 18 Aug 2025).

The transient nature of these shells is also a direct consequence of mixing. Once formed, an oCZ “quickly mixes its interior chemically, flattening both $\nabla_{\mu} = d\ln\mu/d\ln P$6 and $\nabla_{\mu} = d\ln\mu/d\ln P$7 inside,” while leaving neighboring layers largely unmixed because they remain Ledoux-stable (Miszuda, 18 Aug 2025). The shell therefore erases the very density and composition anomalies that created the local opacity peak, while simultaneously steepening $\nabla_{\mu} = d\ln\mu/d\ln P$8 at its boundaries. As the opacity peak weakens, $\nabla_{\mu} = d\ln\mu/d\ln P$9 falls below $\mu$0 and the shell disappears (Miszuda, 18 Aug 2025). The oCZs are thus “short-lived episodes” associated with late stages of Roche-lobe overflow and core rejuvenation rather than permanent structures (Miszuda, 18 Aug 2025).

Boundary dynamics in related work reinforce this picture. Reviews of overshooting emphasize that convective boundaries are not sharp because finite-velocity plumes penetrate into adjacent stable layers; the relevant distinction is between a penetration regime at large Péclet number and an overshoot regime at small Péclet number (Dintrans, 2010). Three-dimensional direct numerical simulations recover primarily the overshoot regime, with a thin thermal boundary layer and intermittent plume penetration rather than a broad, globally adiabatic extension (Dintrans, 2010). By contrast, Anders et al. identify penetration zones in which the mean stratification becomes approximately adiabatic beyond the Schwarzschild boundary, and parameterize the extent with

$\mu$1

showing that stellar-like cases with $\mu$2 can yield extensions of up to $\mu$3–$\mu$4 of a mixing length beyond the formal boundary (Anders et al., 2021). A plausible implication is that some transient oCZs in evolutionary calculations represent the time-dependent interplay of composition stabilization, plume penetration, and local radiative-gradient structure rather than purely local instability in a static background.

4. Dynamical coupling: overshoot, waves, and layered transport

oCZs are not dynamically isolated shells. Overshooting studies show that convective elements entering an adjacent stable layer excite internal gravity waves and that the nature of the interface depends on thermal diffusion, plume filling factor, and dimensionality (Dintrans, 2010). In 3D, plume filling factors are small, extensive adiabatic penetration is disfavored, and the interface behaves more like a thin wave-generating overshoot layer (Dintrans, 2010).

Subsurface off-centre convection in hot stars provides a concrete example of this coupling. In FeCZ simulations, convection in a thin shell beneath the photosphere excites gravity waves that propagate upward into the radiative layer above, and these wave-induced motions are proposed as a candidate for observed photospheric microturbulence (Cantiello et al., 2010). The same convective shell can also act as a dynamo when rotation and shear are present, with magnetic fields reaching equipartition and amplitudes “up to 3 kG,” potentially producing localized magnetic spots at the surface (Cantiello et al., 2010).

Wave transmission through embedded convective layers has also been analyzed in a general rotating Boussinesq framework with alternating radiative and convective layers (Cai et al., 2021). A transient oCZ embedded in a radiative region corresponds to a stable–convective–stable three-layer system, and such an arrangement can display “enhanced wave transmission” by resonant propagation or resonant tunneling when layer thicknesses and clamping-layer properties satisfy specific conditions (Cai et al., 2021). Efficient transmission occurs only when the total number of alternating layers is odd, and the theory predicts that embedded oCZs can function as wave filters or waveguides rather than simple barriers (Cai et al., 2021).

This suggests that transient oCZs can modulate the transport of angular momentum and wave energy through radiative interiors in a time-dependent way. When a shell appears, splits, or vanishes, the star’s layered transmission properties may change sharply, particularly for gravity and gravito-inertial waves whose frequencies lie near resonant windows (Cai et al., 2021).

5. Consequences for mixing, seismology, and observables

The immediate structural consequence of oCZs in mass-accreting stars is a reshaping of the hydrogen and $\mu$5 profiles around the core into a stepped configuration, with each step marking a shell that has homogenized its own material (Miszuda, 18 Aug 2025). The paper explicitly identifies effects on “buoyancy frequency profiles and thus the propagation cavities for g-modes and p-modes,” and argues that sharp $\mu$6 and density discontinuities associated with oCZs contribute to the “asteroseismic fingerprints” of post-interaction stars (Miszuda, 18 Aug 2025).

Convective-penetration models imply a related seismic effect. A MESA solar model with a penetration-zone prescription places a nearly adiabatic, chemically mixed layer beneath the Schwarzschild base of the convection zone and produces an acoustic glitch with a $\mu$7 increase in sound speed below the convection-zone base (Anders et al., 2021). Although presented as a proof of concept, this demonstrates that off-centre or penetrative convective extensions can leave measurable signatures in stratification-sensitive oscillation data (Anders et al., 2021).

In hot massive stars, the FeCZ has been linked to a wide range of surface and wind phenomena: non-thermal line broadening, line-profile variability, discrete absorption components, wind clumping, non-thermal X-rays, and stochastically excited pulsations (Cantiello et al., 2010). The proposed causal chain is subsurface off-centre convection $\mu$8 gravity waves and magnetic fields $\mu$9 photospheric turbulence and wind structure (Cantiello et al., 2010). These phenomena are not claims about mass-accreting oCZs specifically, but they demonstrate that thin off-centre convective shells can have consequences far outside their local thermal domain.

Non-static convection zones can also limit mode coherence. In DAV white dwarfs, the mass of the surface convection zone scales as $\nabla_{\rm rad} \propto \frac{\kappa L P}{m T^4},$0, so small pulsational temperature changes move the convective boundary and alter g-mode phase coherence through cavity-size changes and Doppler shifts of reflected waves (Montgomery et al., 2019). The paper generalizes this mechanism qualitatively to deeper time-dependent convective boundaries. This suggests that transient oCZs, when located within a g-mode cavity, may also broaden mode widths or modulate amplitudes through time-dependent reflection conditions (Montgomery et al., 2019).

6. Numerical modeling, robustness, and unresolved issues

The mass-accretion oCZs reported in recent work were computed with MESA r23.05.1 in a binary setup with an initially $\nabla_{\rm rad} \propto \frac{\kappa L P}{m T^4},$1 primary and $\nabla_{\rm rad} \propto \frac{\kappa L P}{m T^4},$2 secondary, conservative mass transfer, solar composition $\nabla_{\rm rad} \propto \frac{\kappa L P}{m T^4},$3, and Ledoux convection with time-dependent convection, semiconvection, thermohaline mixing, and exponentially decaying overshooting above the H-burning core (Miszuda, 18 Aug 2025). The calculations also imposed a minimum diffusion coefficient of $\nabla_{\rm rad} \propto \frac{\kappa L P}{m T^4},$4 and excluded rotational mixing (Miszuda, 18 Aug 2025).

The crucial methodological result is robustness. The oCZs persist under spatial- and temporal-resolution changes and under several alternative mixing treatments, including shell overshoot, shell over- and undershoot, convective premixing, higher $\nabla_{\rm rad} \propto \frac{\kappa L P}{m T^4},$5, and stronger semiconvection (Miszuda, 18 Aug 2025). Even where internal structure differs, “off-centre convective shells are still present” (Miszuda, 18 Aug 2025). The authors therefore argue that these shells are “not numerical artefacts” but a genuine physical response to mass accretion and boundary mixing (Miszuda, 18 Aug 2025).

At the same time, related literature highlights limitations in current physical modeling. Overshooting reviews emphasize that mixing-length theory treats convective boundaries as local and sharp, whereas realistic plume dynamics are nonlocal and diffusion-dependent (Dintrans, 2010). Convective-penetration simulations show that penetration zones can take “thousands of overturn times to develop,” so long simulations or accelerated evolutionary techniques are required (Anders et al., 2021). Wave-transmission models assume plane-parallel Boussinesq geometry, sharp interfaces, linear waves, and no magnetic fields (Cai et al., 2021). FeCZ simulations use local Cartesian boxes with deliberately unrealistic flux ratios and reduced density contrasts, so their field strengths and velocities are qualitative rather than directly stellar (Cantiello et al., 2010).

A recurrent unresolved issue is how to map between thin, composition-regulated shells in 1D stellar evolution and the multi-dimensional interface physics of penetration, overshoot, and wave excitation. The data support a coherent picture in which transient oCZs are real stellar structures, but the precise extent of their mixing, their lifetime relative to secular evolution, and their coupling to rotation, magnetism, and low-Prandtl-number turbulence remain open problems within the cited work (Miszuda, 18 Aug 2025, Dintrans, 2010, Anders et al., 2021).

7. Conceptual synthesis

Taken together, the cited studies define transient oCZs as detached convective regions whose existence depends on local departures from simple one-zone stability logic. In mass-accreting stars, they arise from a mass-transfer-induced density and opacity enhancement near a retreating core boundary, cross the Ledoux threshold only in a narrow shell, and vanish when their own mixing removes the anomaly that created them (Miszuda, 18 Aug 2025). In massive-star envelopes, an opacity bump can create a shallow subsurface off-centre shell with dynamical consequences disproportionate to its thickness (Cantiello et al., 2010). At convective boundaries more generally, penetration and overshoot blur the distinction between formally radiative and dynamically convective regions (Dintrans, 2010, Anders et al., 2021).

The unifying interpretation is that an oCZ is best understood as a local convective shell embedded in a larger radiative environment, with its onset governed by the competition between radiative-gradient enhancement and stabilizing stratification, and its evolution governed by mixing, diffusion, and nonlocal plume dynamics. This suggests that transient oCZs are not merely peculiar features of specific stellar codes, but a generic manifestation of how stars respond when opacity structure, composition gradients, or boundary energetics briefly favor convection away from the centre and away from an existing main convective region.