Primordial Helium Mass Fraction

Primordial helium mass fraction: what this calculator estimates

Big Bang nucleosynthesis set the first light-element abundances, and this calculator focuses on the two outputs most often compared with observations: primordial helium-4 mass fraction (Y_p) and primordial deuterium abundance, written as D/H. It turns the baryon-to-photon ratio in convenient units (η₁₀), the neutron mean lifetime (τ_n), and any additional relativistic energy density beyond the Standard Model expressed as an effective neutrino contribution (ΔN_ν) into a quick estimate of how much helium and deuterium the early universe should have produced.

The result is a fitting-formula estimate, not a full reaction-network run. That makes it useful for classroom work, quick parameter scans, and checking how sensitive Y_p or D/H is to a small shift in η₁₀, τ_n, or ΔN_ν. For precision comparisons with data, however, you still want a dedicated BBN code and a careful uncertainty budget, because the details of nuclear rates and cosmological assumptions can matter once you start asking for decimal-place accuracy.

Why primordial helium mass fraction tracks the neutron-to-proton balance

In primordial helium calculations, the key story is that most neutrons surviving until nuclei can assemble are quickly captured into helium-4. A convenient way to picture the result is to relate Y_p to the neutron-to-proton ratio at the moment the deuterium bottleneck breaks. Let

x = n/p

Assuming most neutrons go into ⁴He, the helium mass fraction is approximately

The ratio x is governed by weak freeze-out and by the time available for free neutrons to decay before helium forms. That is why τ_n and the expansion rate matter so much: a longer neutron lifetime leaves more neutrons available, while extra radiation speeds up the expansion and can preserve neutrons by shortening the wait. The baryon density still matters, but mostly through the timing of nuclear assembly rather than by directly changing the neutron-to-proton balance.

Helium-4 fitting formula used by this calculator

This calculator uses a common linearized fit around standard cosmological values, with a baseline helium fraction and small corrections from η₁₀, τ_n, and S, where S summarizes any extra relativistic energy density:

Y_p = 0.2485 + 0.0016(η₁₀ − 6) + 0.0002(τ_n − 880) + 0.013(S − 1)

with

S = √(1 + 7ΔN_ν/43)

Here, η₁₀ is the baryon-to-photon ratio scaled as η₁₀ = 10¹⁰η, τ_n is the neutron mean lifetime in seconds, and ΔN_ν is the extra radiation term expressed in neutrino units. Because S = √(1 + 7ΔN_ν/43), increasing ΔN_ν raises the expansion rate, leaves less time for neutrons to decay, and usually pushes Y_p upward. Around the standard cosmological point, each term acts as a small correction rather than a complete reshaping of the prediction.

How primordial deuterium (D/H) behaves alongside helium

Primordial deuterium reacts even more strongly to the baryon density than helium does. When η₁₀ is higher, nuclear reactions run longer and more completely, so deuterium is burned into helium more efficiently and D/H falls. Extra radiation can leave a little more deuterium behind by speeding up the expansion, but η₁₀ is still the dominant lever in most realistic parameter sweeps.

This calculator therefore treats deuterium as the sharper diagnostic and helium as the smoother one. If you are trying to understand which input is driving a result, D/H will usually show the strongest response to changes in the baryon density, while Y_p is the better place to look for shifts in neutron survival and expansion rate.

Interpreting primordial helium mass fraction results

• Y_p (helium-4 mass fraction) is a mass fraction, so values should lie between 0 and 1. Standard-BBN benchmarks typically fall near 0.24–0.25, and a change of 0.001 in Y_p is already meaningful if you are comparing to observational or published values.
• D/H is a number ratio and is often of order a few × 10⁻⁵ in standard cosmology. If your output is far from that range, check the units on η₁₀, τ_n, and ΔN_ν, and make sure you are still inside the fit's intended neighborhood.
• Sensitivity guide: for this calculator, η₁₀ drives D/H most strongly, ΔN_ν has the clearest effect on Y_p, and τ_n shifts Y_p more gently. A useful habit is to change one input at a time so you can see which abundance moves first and how the two outputs respond differently.

Worked example: a near-standard primordial helium prediction

Suppose you enter η₁₀ = 6.1, τ_n = 880 s, and ΔN_ν = 0. Those inputs sit close to the standard reference point, so the calculator should only nudge the answer away from its baseline rather than sending it into an extreme regime.

1. Compute S: S = √(1 + 7×0/43) = 1.
2. Plug into the helium fit:
   • Baseline: 0.2485
   • Baryon term: 0.0016(6.1−6) = 0.00016
   • Neutron lifetime term: 0.0002(880−880) = 0
   • Radiation term: 0.013(1−1) = 0

So Y_p ≈ 0.2485 + 0.00016 = 0.24866. In percent terms, that means about 24.866% of the baryonic mass is in helium-4 in this simplified estimate, while the deuterium estimate stays close to the standard few × 10⁻⁵ level because the baryon density has only shifted slightly.

Primordial helium and deuterium effects at a glance

Assumptions, validity range, and limitations for primordial helium estimates

• Approximate fits: The calculator uses compact fitting relations that mimic the behavior of full BBN calculations near standard parameters. Once you move far away from that region, the linearized terms can drift and nonlinear effects can become more important, especially if several inputs shift at the same time.
• Parameter range: The outputs are most reliable when η₁₀ is in the few-to-10 range, τ_n stays near present experimental values, and ΔN_ν is small to moderate. Very large ΔN_ν can make the simple S mapping and the helium fit less trustworthy.
• Model dependence: Real BBN predictions also depend on nuclear reaction rates, neutrino physics, and the thermal history of the early universe. This calculator does not propagate those uncertainties, so it should be used as a clean pedagogical estimate rather than a final word on primordial abundances.
• Notation: ΔN_ν here is a convenient extra-radiation parameter. Different papers use N_eff, ΔN_eff, or closely related definitions, so compare like with like when you read a paper, a data table, or a published cosmology fit.
• Not a measurement tool: Astronomical determinations of primordial Y_p and D/H include astrophysical systematics such as chemical evolution, ionization corrections, and source selection. The calculator gives the theoretical side of the comparison only, which is helpful for intuition but not a substitute for an observational analysis.

References for primordial helium and deuterium context

For background on primordial helium and deuterium, standard Big Bang nucleosynthesis reviews, particle-data summaries, and modern numerical BBN papers are the right starting points. If you need publication-grade abundances, pair a current BBN code with up-to-date nuclear rates and the exact parameter definitions used in your source, then compare the theoretical output to the observational quantity you actually intend to test.

How to use this primordial helium mass fraction calculator

1. Enter Baryon-to-photon ratio η10 (10⁻¹⁰) as a positive value in the same 10⁻¹⁰ scaling used by the field.
2. Enter Neutron lifetime τₙ (s) in seconds, using the value you want to test against the fit.
3. Enter Extra relativistic species ΔNν as a nonnegative value that represents the extra radiation input.
4. Run the calculation, then try a nearby η10, τn, or ΔNν value to see whether Yp or D/H is the quantity that shifts most for your scenario.