🔭 astrophysics

A data-driven estimate of the protosolar helium mass fraction

This paper presents a data-driven update to the protosolar helium mass fraction, incorporating macroscopic mixing effects and recent helioseismic constraints to derive a revised primordial value of approximately 0.276, while identifying uncertainties in the solar equation of state as the dominant source of error.

Original authors: G. Buldgen, M. Kunitomo, A. Noels, T. Guillot, R. Scuflaire, N. Grevesse

Published 2026-03-03

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: G. Buldgen, M. Kunitomo, A. Noels, T. Guillot, R. Scuflaire, N. Grevesse

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the Sun as a giant, cosmic soup pot that has been simmering for 4.6 billion years. Inside this pot, the ingredients (mostly hydrogen and helium) have been slowly changing. The scientists in this paper are trying to figure out exactly how much helium was in the pot when the Sun was first born (the "protosolar" amount).

Why does this matter? Because knowing the original recipe helps us understand how giant planets like Jupiter and Saturn formed, and it helps us predict how other stars in the universe behave.

Here is the story of their discovery, broken down into simple concepts:

1. The Old Way of Guessing (The "Perfect Stillness" Mistake)

For a long time, scientists tried to guess the Sun's original helium by looking at the Sun today and running backward in time using computer models. They assumed the Sun was a perfectly still pot where heavy stuff (helium) slowly sinks to the bottom due to gravity, like sand settling in a jar of water.

They calculated: "If we see this much helium at the top today, and we know how much sank, we can figure out how much was there at the start."

The Problem: This method assumed the Sun was perfectly calm. But the Sun isn't a still jar; it's a churning, turbulent storm.

2. The New Discovery: The "Stirring Spoon"

The authors of this paper realized that at the very bottom of the Sun's outer layer (the convective zone), there is a lot of turbulence. Think of this turbulence as a giant, invisible spoon constantly stirring the soup.

Without stirring: Helium sinks easily to the bottom.
With stirring: The turbulence keeps the helium mixed up, preventing it from sinking as much as we thought.

Because of this "stirring," the helium at the surface today is actually higher than it would be in a calm model. This changes the math completely. If the helium didn't sink as much as we thought, it means the Sun started with less helium than previous models suggested.

3. The "Chemical Clocks": Lithium and Beryllium

How do they know how strong the "stirring spoon" is? They use two special ingredients as clocks: Lithium and Beryllium.

These elements are very fragile. If the Sun gets too hot or mixes too vigorously deep inside, they get destroyed (burned up).
The Sun has much less Lithium and Beryllium on its surface today than it did when it was born.
By measuring exactly how much is missing, the scientists can calibrate the "stirring spoon." They can say, "Okay, the Sun must have been stirred this hard to destroy exactly this amount of Lithium."

Once they know the stirring speed, they can accurately calculate how much helium was prevented from sinking.

4. The Equation of State: The "Recipe Book" Problem

There is another hurdle. To calculate how the Sun behaves, scientists need a "recipe book" called the Equation of State. This book tells them how the gas inside the Sun reacts to heat and pressure.

Different recipe books give slightly different answers.
The paper points out that the biggest uncertainty in their final number isn't the stirring or the measurements; it's that we don't have the perfect recipe book for solar gas yet. It's like trying to bake a cake without knowing exactly how much the flour weighs.

5. The Final Result: A Lower Number

When they combined all these new ideas (the stirring spoon, the chemical clocks, and the latest measurements of the Sun's surface), they got a new answer for the Sun's original helium content.

Old Estimate: About 27.8% helium.
New Estimate: About 27.6% helium (give or take a tiny bit).

It sounds like a small difference, but in the world of astronomy, this is a huge shift. It means the Sun started with slightly less helium than we thought.

Why Should You Care?

Think of the Sun as the "control group" for the entire universe. If we get the Sun's recipe wrong, our models for how planets form, how stars age, and how the universe evolves are all slightly off.

By realizing that the Sun is "stirred" and using Lithium and Beryllium as clues, these scientists have sharpened our picture of the Sun's history. It's a reminder that even something as familiar as our Sun is full of hidden turbulence and secrets that we are only just beginning to understand.

In a nutshell: The Sun isn't a calm jar of soup; it's a churning pot. Because of this churning, less helium sank to the bottom than we thought, meaning the Sun started with slightly less helium than our old maps said. We used the "missing" Lithium and Beryllium to figure out exactly how much churning was happening.

1. Problem Statement

The protosolar helium mass fraction ( $Y_P$ ) is a fundamental parameter for modeling the evolution of the Sun, giant planets, and stellar populations. However, current estimates of $Y_P$ suffer from significant uncertainties due to two main factors:

Reliance on Simplified Physics: Traditional estimates often rely on Standard Solar Models (SSMs) that assume only microscopic diffusion (gravitational settling) and ignore macroscopic mixing processes at the base of the convective zone (BCZ).
Equation of State (EOS) Uncertainties: The inference of the current surface helium mass fraction ( $Y_{CZ}$ ) from helioseismology is highly sensitive to the Equation of State of solar plasma. Different EOS choices (e.g., OPAL vs. MHD vs. SAHA-S) yield significantly different $Y_{CZ}$ values.
Calibration Gaps: Previous studies often failed to simultaneously reproduce the observed depletion of light elements (Lithium and Beryllium) alongside the current $Y_{CZ}$ , leading to biased estimates of the initial helium abundance.

2. Methodology

The authors employ a data-driven approach combining helioseismic inversions, spectroscopic abundances, and advanced solar modeling to update the $Y_P$ estimate.

Solar Modeling:
- Models were computed using the Liège Stellar Evolution Code (CLES).
- Physics Ingredients: The study utilizes AAG21 solar abundances, OPAL and LANL/OPLIB opacities, the SAHA-S (v7) and FreeEOS equations of state, and NACRE II nuclear reaction rates.
- Transport Mechanisms: The models explicitly include microscopic diffusion (Thoul et al. 1994) and macroscopic turbulent diffusion at the BCZ (following Proffitt & Michaud 1991). The turbulence is parameterized by a diffusion coefficient $D \propto (\rho_{BCZ}/\rho)^X$ .
- Calibration: Models are calibrated to reproduce the solar radius, luminosity, and surface metallicity at the current solar age. Crucially, the models are tested against the observed depletion of Lithium ( $A(\text{Li}) = 0.96 \pm 0.05$ ) and Beryllium ( $A(\text{Be}) = 1.21 \pm 0.05$ ).
Helioseismic Constraints:
- The study utilizes $\Gamma_1$ inversions from Buldgen et al. (2024a) to determine the current helium mass fraction in the convective envelope ( $Y_{CZ}$ ).
- A meta-analysis of literature values for $Y_{CZ}$ (spanning 30 years of data from Basu & Antia, Di Mauro, Vorontsov, etc.) is performed to establish a robust range.
Derivation of $Y_P$ :
- The protosolar abundance is derived via the relation: $Y_P = Y_{CZ} - (Y_{CZ} - Y_P)$ .
- The term $(Y_{CZ} - Y_P)$ represents the helium depletion over the Sun's lifetime, which is determined by the efficiency of turbulent mixing at the BCZ.

3. Key Contributions & Findings

A. The Critical Role of Macroscopic Mixing

The study demonstrates that macroscopic mixing at the BCZ cannot be neglected.

In Standard Solar Models (SSMs) without turbulence, helium settling is overestimated.
Including turbulent diffusion significantly reduces the efficiency of helium settling.
Quantitative Impact: Using an SSM leads to an overestimation of the helium depletion ( $Y_{CZ} - Y_P$ ) by approximately 40% compared to models calibrated with light element depletion. This results in a significantly overestimated $Y_P$ if turbulence is ignored.

B. Lithium and Beryllium as Calibrators

The authors show that the observed depletion of Lithium and Beryllium serves as an efficient empirical calibrator for the mixing efficiency at the BCZ.

While overshooting alone cannot reproduce both Li and Be depletion simultaneously (and often over-depletes Li in young twins), a calibrated turbulent diffusion model can reproduce the observed values.
The efficiency of this mixing directly dictates the amount of helium retained in the convective envelope, thereby fixing the $Y_{CZ} - Y_P$ difference.

C. Revised Protosolar Helium Abundance ( $Y_P$ )

By combining the latest helioseismic $Y_{CZ}$ estimates with the turbulence-calibrated depletion values, the authors derive updated intervals for $Y_P$ :

Primary Result (Data-Driven):
Using the latest $Y_{CZ}$ estimates from Buldgen et al. (2024a) and strict constraints on Li/Be depletion:
$Y_P = 0.27575 \pm 0.00315$
This value is slightly lower than previous estimates (e.g., Serenelli & Basu 2010) but consistent within uncertainties. The reduction is due to a higher inferred $Y_{CZ}$ in recent studies combined with a lower depletion rate caused by turbulence.
Alternative Scenario (Traditional $Y_{CZ}$ ):
If the traditional $Y_{CZ}$ value from Basu & Antia (2004) is used, but macroscopic mixing is still included:
$Y_P = 0.2669 \pm 0.00415$
This highlights that even with older surface helium data, accounting for turbulence lowers the primordial estimate significantly compared to SSMs.
Conservative Range:
A very conservative estimate spanning all literature $Y_{CZ}$ values and depletion scenarios is $0.2552 \leq Y_P \leq 0.2792$ .

D. Dominant Sources of Uncertainty

The study identifies the Equation of State (EOS) as the dominant source of uncertainty.

The uncertainty in $Y_{CZ}$ derived from helioseismology is primarily driven by the EOS (e.g., differences between OPAL, MHD, and SAHA-S).
Secondary uncertainties arise from the inference techniques (inversion methods) and surface effect corrections.
The uncertainty in the transport physics (turbulence) is secondary to the EOS uncertainty but remains crucial for the depletion term.

4. Significance and Implications

Planetary Modeling: The revised $Y_P$ values provide more accurate initial conditions for modeling the formation and evolution of giant planets in the Solar System.
Stellar Physics: The work underscores that SSMs are insufficient for precise chemical evolution studies. Future models must include macroscopic mixing calibrated against light element depletion to accurately predict stellar ages and internal structures.
Future Directions: To further reduce uncertainties in $Y_P$ $Y_{P}$ , the paper calls for:
1. Improved Equations of State for solar plasma.
2. Better modeling of surface effects (using 3D hydrodynamical simulations).
3. More precise measurements of the current solar photospheric Beryllium abundance, which is a critical anchor for mixing efficiency.

In conclusion, this paper establishes that ignoring macroscopic mixing leads to a systematic overestimation of the protosolar helium abundance. By integrating helioseismic data with light element depletion constraints, the authors provide a more robust and precise determination of $Y_P$ , essential for high-precision astrophysics.