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

This paper proposes a refined criterion for chemically driven convection, demonstrating that while mean molecular weight inversions above the red giant branch bump are insufficient to sustain mixing, rapid carbon production during the helium core flash could trigger a steady convective region that significantly alters our understanding of this stellar event.

The Gist

The Big Picture: Stars are Like Layered Cakes (Sometimes)

Imagine a star as a giant, glowing cake. Usually, we think of stars as having layers that stay put: heavy stuff sinks to the bottom, and light stuff floats on top. This is how nature likes to work; it's stable.

However, inside a star, nuclear reactions are constantly cooking up new ingredients. Sometimes, these reactions create "heavy" elements right on top of "lighter" ones. It's like if you baked a cake and suddenly dropped a heavy layer of lead into the middle of the fluffy frosting.

In physics, this is called a mean molecular weight inversion. It's unstable. The heavy stuff wants to sink, and the light stuff wants to rise. This creates a tug-of-war that can cause the star's interior to mix, or "stir," violently.

The Old Rulebook vs. The New Discovery

For decades, astronomers used a standard rulebook (called the Ledoux criterion) to predict when this stirring would happen. Think of this rulebook like a traffic light. It says: "If the difference in weight between the layers is small, the traffic is green (stable). If the difference is huge, the traffic is red (unstable, and mixing happens)."

The authors of this paper, M. Miguel Ocampo and Marcelo M. Miller Bertolami, looked at the math behind that traffic light and realized it's broken.

They found that the old rulebook is too strict. It assumes that for mixing to happen, the "weight difference" needs to be massive. But their new calculations show that mixing can actually start with a tiny weight difference, provided the "stirring" happens fast enough.

The Analogy:
Imagine trying to push a heavy boulder up a hill.

• The Old View: You need a giant, muscular superhero to push it. If you are just a regular person, the boulder won't move.
• The New View: If you give that regular person a really fast, powerful jetpack (representing the physics of how heat moves), they can actually push that boulder up the hill even if they aren't super strong.

The paper proves that inside stars, the "jetpack" (fast convection) is often stronger than we thought, meaning mixing can happen much more easily than the old rules predicted.

The Two Test Cases: The "Bump" and The "Flash"

To test their new theory, the authors looked at two specific moments in a star's life:

1. The Red Giant Branch "Bump" (The Slow Stir)

As a star ages, it swells up into a Red Giant. At a certain point, it hits a "bump" in its evolution. Here, a nuclear reaction creates a tiny bit of extra hydrogen, making a small weight inversion.

• The Finding: The authors checked if this tiny bump could trigger the "fast mixing" they discovered.
• The Result: No. The weight difference here is too weak and disappears too quickly. It's like trying to start a forest fire with a single spark in the rain. The mixing happens, but it's slow and weak (like a gentle simmer), not the violent storm the old simulations suggested.
• Conclusion: The extra mixing seen in stars at this stage is likely caused by something else, perhaps weak magnetic fields acting like a whisk, rather than this new type of fast convection.

2. The Helium Core Flash (The Explosive Mix)

Later in the star's life, the core runs out of hydrogen and starts fusing helium. This happens in a sudden, violent burst called the "Helium Flash."

• The Finding: Here, the star is churning out Carbon at a furious rate. This creates a massive weight difference.
• The Result: Yes! The new rules show that this environment is perfect for the "fast mixing." The nuclear reactions are so strong they can keep the "heavy layer" on top of the "light layer" constantly, feeding the instability.
• The Implication: This suggests that during the Helium Flash, the star doesn't just have a small convective zone; it might develop a deep, fast, adiabatic (super-hot) mixing zone that reaches all the way toward the center.

Why This Changes Everything

If the authors are right, our understanding of how stars die and explode is wrong.

1. The "Flash" might be different: Standard models say the helium flash happens in a shell around the core, creating a series of smaller "sub-flashes." But if this new fast mixing happens, it could carry the heat and reaction products straight to the center. This might mean the flash happens all at once in the very center, and those little "sub-flashes" never happen.
2. Hydrodynamics vs. Reality: Computer simulations that try to model stars in 3D (like video game physics) have shown fast, violent mixing in these areas. But the 1D models (the standard "rulebook" models) couldn't explain why that was happening. This paper provides the missing link: Chemical gradients are strong enough to drive fast convection, even when the old rules said they weren't.

The Bottom Line

The universe is a bit more chaotic than we thought. Inside stars, the "heavy stuff on top of light stuff" doesn't need to be a huge imbalance to cause a stir. If the star is cooking up new elements fast enough, even a small imbalance can trigger a massive, fast-moving mixing event.

This discovery suggests that the violent "flashes" in dying stars might be much more dramatic and central than our current textbooks describe, and it explains why high-tech computer simulations have been seeing things that the simple math couldn't predict.

Technical Summary

1. Problem Statement

Turbulent mixing in stellar interiors remains a primary uncertainty in stellar evolution modeling. Standard 1D codes rely on the Ledoux instability criterion to identify regions where mixing occurs due to mean molecular weight ($\mu$) inversions (where heavier material sits atop lighter material).

• Current Assumption: The standard criterion used in stellar evolution codes is $\nabla_{rad} > \nabla_{ad} + B$, where $\nabla_{rad}$ is the radiative gradient, $\nabla_{ad}$ is the adiabatic gradient, and $B$ is the Ledoux term (representing the chemical gradient). This assumes that if the radiative gradient required to transport energy is unstable, mixing will occur.
• The Issue: This assumption fails to distinguish between two distinct physical regimes:
  1. Thermohaline Mixing: Slow, diffusive mixing driven by the difference in diffusion rates between heat and chemical species.
  2. Chemically Driven Convection: Fast, nearly adiabatic convection similar to standard thermal convection but driven by composition gradients.
• The Discrepancy: Hydrodynamical simulations (2D/3D) of events like the Red Giant Branch Bump (RGBB) and the Helium (He) Core Flash show fast, large-scale mixing plumes that standard thermohaline prescriptions cannot explain. Conversely, standard 1D codes often fail to predict these fast motions because they rely on the simplified Ledoux criterion (Eq. 3), which may be too restrictive for fast convection.

2. Methodology

The authors developed a theoretical framework and performed numerical experiments to re-evaluate mixing criteria:

• Extended Mixing Length Theory (MLT): They extended the standard MLT equations to include chemical gradients ($B \neq 0$). They derived a system of equations relating the actual temperature gradient ($\nabla$), the element gradient ($\nabla_e$), the adiabatic gradient, and the radiative gradient.
• Dimensionless Velocity Analysis: By combining the MLT equations, they derived a quartic equation for a dimensionless velocity $\nu$. This equation reveals that for specific ranges of parameters, three distinct solutions can coexist:
  1. Slow solution: Corresponds to thermohaline mixing (radiative-like gradients).
  2. Intermediate solution: A transition state.
  3. Fast solution: Corresponds to chemically driven convection (adiabatic-like gradients).
• Derivation of a New Criterion: They identified critical lines in the parameter space ($\nabla_{rad} - \nabla_{ad}$ vs. $B$) that separate the regime of single solutions (slow only) from the regime of three coexisting solutions. They derived an approximate analytical expression (Eq. 12) for the critical chemical gradient required to trigger the fast convective solution.
• Numerical Simulations: Using the LPCODE stellar evolution code, they modeled a $1 M_\odot$ star from the Zero Age Main Sequence through the RGBB and the He-core flash.
  • They tested the stability of layers at the RGBB and during the He-flash using the new criterion.
  • They performed "experiments" where they artificially enhanced mixing coefficients to test if chemical gradients could be sustained against homogenization.

3. Key Contributions

1. Refutation of the Standard Criterion: The paper demonstrates that the standard criterion ($\nabla_{rad} > \nabla_{ad} + B$) is insufficient. It fails to identify regions where fast, adiabatic chemically driven convection is possible, even when the chemical gradient is much smaller than previously thought necessary.
2. New Instability Criterion (Eq. 12): The authors derived a new, more accurate criterion for the onset of fast chemically driven convection. This criterion depends on the ratio of free-fall to thermal adjustment timescales ($U$) and the mixing length ($l_m$). It shows that fast convection can occur at much lower chemical gradients than the standard Ledoux limit if the mixing elements are large enough (comparable to the pressure scale height).
3. Resolution of the "RGBB Mixing Problem": The study addresses the long-standing discrepancy where observed abundances in RGBB stars require mixing efficiencies far higher than standard thermohaline theory predicts.
4. Reinterpretation of the He-Core Flash: The paper provides a physical mechanism for the fast mixing observed in hydrodynamical simulations of the He-core flash, linking it to chemically driven convection rather than just thermohaline instability.

4. Key Results

A. Red Giant Branch Bump (RGBB)

• Instability: At the RGBB, the $^3$He($^3$He, 2p)$^4$He reaction creates a small $\mu$-inversion. The new criterion suggests this region could be unstable to fast convection if the mixing length is large ($l_m \sim H_P$).
• Self-Quenching: However, numerical experiments show that the chemical gradient generated by the reaction is too weak and short-lived. Once fast mixing begins, it homogenizes the layer faster than the reaction can replenish the gradient.
• Conclusion: Chemically driven convection is unlikely to be the primary driver of mixing at the RGBB. The observed mixing is more likely explained by enhanced thermohaline mixing due to weak magnetic fields, rather than a transition to fast convection.

B. Helium (He) Core Flash

• Instability: During the He-core flash, rapid carbon production via the triple-$\alpha$ reaction creates a strong $\mu$-inversion below the thermal convective zone.
• Sustainability: Unlike the RGBB, the nuclear burning rate is high enough to sustain the chemical gradient against the homogenizing effect of fast mixing.
• Adiabatic Convection: The new criterion predicts that a steady-state, adiabatic, chemically driven convective zone can develop below the thermal flash zone.
• Implications: This zone would significantly alter the thermal structure of the core. It could shift the temperature maximum inward, potentially causing the He-flash to ignite at the center of the star. This would eliminate the "He sub-flashes" (secondary flashes) typically seen in standard 1D models.

5. Significance and Implications

• Stellar Evolution Models: The findings suggest that standard 1D stellar evolution codes are missing a critical mixing mechanism (fast chemically driven convection) in specific high-gradient environments.
• He-Core Flash Physics: The possibility of a central ignition and the disappearance of sub-flashes fundamentally changes our understanding of the late-stage evolution of low-mass stars and the resulting white dwarf properties.
• Hydrodynamical Simulations: The work offers a theoretical justification for the fast, large-scale plumes observed in 2D/3D hydrodynamical simulations of the He-flash and RGBB, which were previously difficult to reconcile with 1D thermohaline models.
• Future Directions: The authors call for a reanalysis of inverse chemical gradients across the entire mass range of stellar evolution, including Asymptotic Giant Branch (AGB) stars and white dwarf crystallization fronts, using the new instability criterion.

In summary, the paper argues that chemically driven convection is a viable, fast mixing mechanism in stellar interiors that is currently overlooked by standard stability criteria. While it may not explain RGBB mixing, it likely plays a decisive role in the dynamics of the He-core flash, potentially rewriting the standard model of this evolutionary phase.