On the possibility of chemically driven convection in red giants. Implications for the He-core flash and mixing above the Red Giant Branch Bump

Marcelo M. Miller Bertolami; M. Miguel Ocampo

arxiv: 2604.12117 · v1 · submitted 2026-04-13 · 🌌 astro-ph.SR · physics.flu-dyn

On the possibility of chemically driven convection in red giants. Implications for the He-core flash and mixing above the Red Giant Branch Bump

M. Miguel Ocampo , Marcelo M. Miller Bertolami This is my paper

Pith reviewed 2026-05-10 15:24 UTC · model grok-4.3

classification 🌌 astro-ph.SR physics.flu-dyn

keywords chemically driven convectionred giantshelium core flashmean molecular weightstellar mixingRGB bumpthermohaline instability

0 comments

The pith

Chemically driven convection in red giants can occur with much smaller mean molecular weight inversions than standard criteria require, potentially impacting the helium core flash.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper derives an alternative criterion showing that chemically driven convection can be sustained by smaller mean molecular weight inversions than typically assumed in stellar models. This criterion indicates that while the inversion above the red giant branch bump is too weak and short-lived for steady convection, rapid carbon production during the helium core flash can maintain a steady chemically driven convective region. A reader would care because this affects predictions for mixing processes and the behavior of the helium flash in red giant stars.

Core claim

We demonstrate that the standard criterion adopted in stellar evolution calculations does not accurately distinguish between thermohaline and Rayleigh-Taylor instabilities. We derive an alternative criterion and show that chemically driven convection in stellar interiors might be viable under much smaller mean molecular weight inversions than it is normally assumed. We find that the inversion at the base of the convective envelope above the RGBB is too weak and short-lived to sustain steady-state convection. In contrast, rapid carbon production at the base of the He-flash-driven convective zone can maintain a steady chemically driven convective region.

What carries the argument

Alternative criterion for chemically driven convection based on the size and persistence of mean molecular weight inversions.

If this is right

The inversion above the RGBB cannot sustain steady convection.
Rapid carbon production maintains steady convection during the He-core flash.
This could significantly alter our understanding of the He-core flash.
Further study of this process is needed.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Incorporating this criterion into stellar evolution codes could improve accuracy in modeling red giant phases and surface abundances.
It may help explain certain observed chemical anomalies in red giant stars if the additional mixing is confirmed.
Self-consistent modeling including the feedback of mixing on inversion profiles would be a natural next step.

Load-bearing premise

The mean molecular weight inversion profiles and timescales from standard stellar evolution models are accurate enough to determine if convection can be sustained without the mixing altering those profiles.

What would settle it

Numerical simulations of the helium core flash that test whether a steady chemically driven convective region forms due to carbon production.

Figures

Figures reproduced from arXiv: 2604.12117 by Marcelo M. Miller Bertolami, M. Miguel Ocampo.

**Figure 1.** Figure 1: Schematic representation of the different regions in the (∇rad − ∇ad) − B space where eqs. 4 and 5 have one or three solutions. In the figure only the quadrant where W = (∇rad − ∇ad) < 0 and B < 0 is displayed. ltherm = (κT ν/NT 2 ) 1/4 where κT is the thermal diffusivity, ν the kinematic viscosity, and NT 2 the squared thermal Brunt-V¨ais¨al¨a frequency (P. Garaud 2018, 2020). At each layer of the star,… view at source ↗

**Figure 2.** Figure 2: Values of the dimensionless velocity ν = (∇−∇e −B) 1/2 for different values of U. The dashed cyan lines indicate the location where the intermediate solution (or the only solution for large values of U) fulfills ∇ = ∇ad + B (neutrality according to the Ledoux criterion). The criterion usually adopted in stellar evolution computations, given by Eq. 3, is shown as a solid blue line. Dashed white lines in the… view at source ↗

**Figure 3.** Figure 3: Values for the sub-adiabacity ∇ − ∇ad for U = 10−2 and U = 10−9 when the slow (left panels), intermediate (middle panels) or fast solution (right panels) is chosen in the region between Eqs. 11 (solid) and 12 (dashed). In the bottom panels, Eq. 11 is not displayed since it overlaps with the vertical axis at that scale. Note that the color bar range of the lower-left panel is linear to highlight the radiati… view at source ↗

**Figure 4.** Figure 4: ). The initial metallicity of the models was taken as Z = 0.01. Models were computed under very standard assumptions for convection, i.e., Schwarzschild Criterion plus MLT, with no additional convective boundary mixing unless otherwise stated. Thermohaline mixing was included in the evolution as in F. C. Wachlin et al. (2014) adopting the state-of-the-art prescription of J. M. Brown et al. (2013). Rece… view at source ↗

**Figure 5.** Figure 5: Physical conditions in the stellar interior as a function of the radius r after the bump in the RGB. Upper panel: Abundances of H (above its surface value δH=H-Hsur), 3He, and 13C. Medium panel: Layer-by-layer value of the critical value of −B for the onset of (fast) chemically driven convection. Lower panel: Actual value of the instability discriminant −B − D(W, U). The inner regions that are either stab… view at source ↗

**Figure 6.** Figure 6: Temporal evolution of the chemical instabilities after the RGBB with the mixing artificially increased to fast convective-like velocities. In a span of a few months, the region is homogenized and, thus, B becomes lower than the Bcrit defined by D(W, U) in a short period of time. Blue and yellow are used to highlight the regions in the curve that would be stable or unstable according to Eq. 12 respectively.… view at source ↗

**Figure 7.** Figure 7: Same as [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗

**Figure 8.** Figure 8: Same as [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗

**Figure 9.** Figure 9 [PITH_FULL_IMAGE:figures/full_fig_p012_9.png] view at source ↗

read the original abstract

Turbulent mixing remains one of the primary uncertainties in the modeling of stellar interiors. In stellar evolution simulations, regions where mixing occurs are typically identified using instability criteria. A particularly interesting situation arises when nuclear reactions produce inversions in the mean molecular weight within stellar interiors. Under these conditions, the material can become unstable to either thermohaline or a Rayleigh-Taylor instabilities. We demonstrate that the standard criterion adopted in stellar evolution calculations does not accurately distinguish between these two regimes. We derive an alternative criterion and show that chemically driven convection in stellar interiors might be viable under much smaller mean molecular weight inversions than it is normally assumed. We investigate whether inversions in the mean molecular weight can trigger chemically driven convection above the red giant branch bump (RGBB) or during the helium core flash. We find that the inversion at the base of the convective envelope above the RGBB is too weak and short-lived to sustain steady-state convection. In contrast, rapid carbon production at the base of the He-flash-driven convective zone can maintain a steady chemically driven convective region. This process could significantly alter our understanding of the He-core flash and warrants further study.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper gives a sharper criterion for when mean-molecular-weight inversions drive convection in red giants, but the helium-flash case needs self-consistent modeling to confirm.

read the letter

The paper's main contribution is a new analytic criterion that separates thermohaline from Rayleigh-Taylor behavior for mean-molecular-weight inversions, and the claim that this makes chemically driven convection possible during the helium flash. They do a clean job of showing why the standard threshold misses the distinction and then test it on the RGB bump and the flash. The RGBB case comes out negative, which matches what people expect, while the flash gets a positive because of the carbon production rate. The derivation is the part that holds up. It is grounded in the physics of the two instabilities and gives a concrete way to decide which one applies. The weak point is the helium flash section. They take the inversion profiles straight from ordinary stellar models and argue that carbon burning keeps the inversion going even with the new mixing. But those profiles do not include the mixing, so the amplitude and duration could change once the convection is active. That makes the steady-state claim conditional rather than demonstrated. The work is for stellar modelers who care about mixing during the flash and its effect on later stages. Anyone updating convective criteria in codes would find the distinction useful to consider. The equations and logic are presented without obvious internal contradictions, and the citations to prior work on thermohaline mixing look standard. I recommend sending it for peer review. The idea is specific enough that referees can check the derivation and the self-consistency issue directly.

Referee Report

2 major / 2 minor

Summary. The paper derives an alternative instability criterion for regions with mean-molecular-weight inversions produced by nuclear burning, arguing that it better distinguishes thermohaline from Rayleigh-Taylor regimes and permits chemically driven convection at smaller inversions than the standard criterion. Applying the new criterion to standard stellar-evolution profiles, the authors conclude that the μ inversion above the RGB bump is too weak and transient to drive steady convection, whereas rapid carbon production during the helium core flash can sustain a steady chemically driven convective zone at the base of the flash-driven convection.

Significance. If the new criterion is shown to be robust and the helium-flash conclusion survives self-consistent modeling, the result would affect mixing prescriptions in stellar-evolution codes and could revise predictions for the helium-core flash and subsequent evolution. The work highlights a regime where nuclear-driven μ inversions may drive additional transport not captured by current Ledoux or Schwarzschild implementations.

major comments (2)

The central claim for the helium-core flash (that rapid carbon production maintains a steady chemically driven region) rests on μ-inversion profiles and timescales taken directly from unmodified stellar-evolution calculations. Because the proposed mixing is not included in those calculations, the amplitude, spatial extent, and lifetime of the inversion are not guaranteed to remain above the derived threshold once the new transport operates; a self-consistent test is required to substantiate the steady-state conclusion.
The manuscript states that the standard criterion fails to distinguish thermohaline from Rayleigh-Taylor regimes but does not provide a side-by-side comparison of the new criterion against the Ledoux and Schwarzschild criteria (including the explicit form of the new threshold and the numerical values of the critical μ gradients) in the regions of interest.

minor comments (2)

The abstract and introduction would benefit from a concise statement of the quantitative threshold (e.g., the minimum |∇μ| or inversion amplitude) required by the new criterion versus the standard one.
Figure captions and axis labels should explicitly indicate whether the plotted μ profiles are taken from standard models or from any test calculations that include the proposed mixing.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and insightful comments on our manuscript. We have addressed each major point below and revised the paper to incorporate clarifications and additional comparisons where appropriate.

read point-by-point responses

Referee: The central claim for the helium-core flash (that rapid carbon production maintains a steady chemically driven region) rests on μ-inversion profiles and timescales taken directly from unmodified stellar-evolution calculations. Because the proposed mixing is not included in those calculations, the amplitude, spatial extent, and lifetime of the inversion are not guaranteed to remain above the derived threshold once the new transport operates; a self-consistent test is required to substantiate the steady-state conclusion.

Authors: We agree that the analysis uses μ-inversion profiles from standard stellar-evolution calculations without the proposed mixing included, which represents a limitation for claiming a fully steady state. However, the helium-core flash is driven by extremely rapid nuclear timescales (on the order of hours to days for carbon production), which are orders of magnitude shorter than the convective turnover or diffusion timescales associated with the chemically driven instability we derive. This rapid replenishment of the μ inversion supports the possibility of sustained convection even if some smoothing occurs. In the revised manuscript we have added an explicit discussion paragraph in Section 4.2 acknowledging this caveat, quantifying the relevant timescales, and identifying self-consistent modeling as an important direction for future work. We view the current result as a first indication that the new criterion permits such a region rather than a definitive demonstration of its steady-state properties. revision: partial
Referee: The manuscript states that the standard criterion fails to distinguish thermohaline from Rayleigh-Taylor regimes but does not provide a side-by-side comparison of the new criterion against the Ledoux and Schwarzschild criteria (including the explicit form of the new threshold and the numerical values of the critical μ gradients) in the regions of interest.

Authors: We appreciate this observation. Although the manuscript contrasts the regimes conceptually, we agree that an explicit side-by-side presentation strengthens the argument. The revised manuscript now includes a new subsection (Section 3.2) that tabulates the explicit mathematical forms of the Ledoux, Schwarzschild, and our new instability criteria. It also provides numerical values of the critical μ gradients evaluated at the base of the convective envelope above the RGB bump and at the base of the flash-driven convection zone during the helium core flash, demonstrating that our threshold permits instability at smaller inversions than the standard Ledoux criterion. revision: yes

Circularity Check

0 steps flagged

No circularity: new instability criterion derived independently and applied to external model profiles

full rationale

The paper's central step is the derivation of an alternative criterion to distinguish thermohaline versus Rayleigh-Taylor regimes under mean-molecular-weight inversions. This is presented as a first-principles result rather than a re-expression of any fitted parameter, self-cited uniqueness theorem, or ansatz imported from prior work by the same authors. The subsequent application to RGBB and He-flash profiles extracts inversion amplitudes and timescales directly from unmodified stellar-evolution calculations; these profiles are treated as external input and are not redefined or refitted within the paper to produce the claimed steady convection. No self-definitional loop, fitted-input prediction, or renaming of known results occurs. The analysis therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The work relies on standard stellar-structure equations and nuclear-reaction networks already present in the literature; no new free parameters, ad-hoc axioms, or invented entities are introduced in the abstract.

axioms (1)

domain assumption Standard stellar evolution codes correctly compute mean-molecular-weight profiles and their time evolution in the absence of the new mixing.
The conclusions about inversion strength and duration are drawn from existing models; any feedback from the proposed convection would invalidate the input profiles.

pith-pipeline@v0.9.0 · 5518 in / 1372 out tokens · 22676 ms · 2026-05-10T15:24:43.228363+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

3 extracted references · 3 canonical work pages

[1]

, " * write output.state after.block = add.period write newline

ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...

work page
[2]

write newline

" write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...

work page
[3]

3=r] ]

thebibliography [1] 20pt to REFERENCES 6pt =0pt \@twocolumntrue 12pt -12pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key o...

work page doi:10.1093/mnras/stu2238 2017

[1] [1]

, " * write output.state after.block = add.period write newline

ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...

work page

[2] [2]

write newline

" write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...

work page

[3] [3]

3=r] ]

thebibliography [1] 20pt to REFERENCES 6pt =0pt \@twocolumntrue 12pt -12pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key o...

work page doi:10.1093/mnras/stu2238 2017