PhysSandbox
Classical MechanicsWaves & SoundElectricity & MagnetismOptics & LightGravity & OrbitsLabs
🌙Astronomy & The Sky🌡️Thermodynamics🌍Biophysics, Fluids & Geoscience📐Math Visualization🔧Engineering🧪Chemistry

Related simulators

Continue with similar topics in this category — or all 51 in Astronomy & The Sky.

View category →
NewSchool

Supernova Light Curves

Launch Simulator

Schematic Ia rise/decay vs II-P plateau; standard-candle note.

NewSchool

Nuclear Binding Curve

Launch Simulator

Qualitative B/A vs A with fusion/fission context.

NewUniversity / research

CMB Power Spectrum (Acoustic Peaks)

Launch Simulator

Cosmic Microwave Background temperature D_ℓ vs ℓ with Sakharov peaks: tune Ω_b h², Ω_c h², n_s, A_s, τ, h and watch the parity flip between odd / even peaks, the Silk damping tail, and the Sachs–Wolfe plateau move. Pedagogical parametric ΛCDM model.

NewUniversity / research

Pair Production Threshold & σ(E_γ)

Launch Simulator

Pair creation channels γ + nucleus → e⁺e⁻ (Bethe–Heitler, σ ∝ Z²), the higher-threshold triplet γ + e⁻ → e⁻e⁺e⁻, and Breit–Wheeler γγ → e⁺e⁻ that limits TeV photons against the cosmic background. Live threshold marker, log-σ curve, and material presets (H, C, Al, Cu, Pb) make the 1.022 MeV / 4 m_e c² thresholds intuitive.

NewUniversity / research

Hydrogen Saha Equilibrium (recombination)

Launch Simulator

n_p n_e / n_H vs T and ionized fraction X; steep Boltzmann knee ~3000 K scale (density dependent).

NewUniversity / research

Neutron star TOV: toy M–R (polytrope)

Launch Simulator

Schwarzschild TOV for polytropes n = 1 or 3/2; K normalized so the ρ_c scan peaks near ~2 M☉ (schematic only).

PhysSandbox

Interactive physics, chemistry, and engineering simulators for students, teachers, and curious minds.

Physics

  • Classical Mechanics
  • Waves & Sound
  • Electricity & Magnetism

Science

  • Optics & Light
  • Gravity & Orbits
  • Astronomy & The Sky

More

  • Thermodynamics
  • Biophysics, Fluids & Geoscience
  • Math Visualization
  • Engineering
  • Chemistry

© 2026 PhysSandbox. Free interactive science simulators.

PrivacyTermsContact
Home/Astronomy & The Sky/Big Bang Nucleosynthesis (BBN)

Big Bang Nucleosynthesis (BBN)

Light-element abundance curves H, ⁴He, D, ³He, ⁷Li vs cosmic time / temperature. Weak freeze-out, neutron decay gap, deuterium bottleneck → Y_p ≈ 0.245. Slide η₁₀ and N_eff over the classic BBN curves; observed values overlaid.

BBN parameters

6.1
3.046

Shortcuts

  • •Slide η₁₀ to see how the surviving abundance of D drops as you add more baryons

Measured values

Y_p (He-4 mass)0.2486
D / H2.532×10⁻⁵
³He / H1.040×10⁻⁵
⁷Li / H1.86×10⁻¹⁰
(n/p)_freeze ≈ exp(−Q/T_F)0.199
(n/p) at BBN onset0.150
Y_p from 2(n/p)/(1+n/p)0.261
Ω_b h² (from η₁₀)0.0224

About this model

Light-element Big Bang nucleosynthesis (BBN) curves track mass fractions of H, ⁴He, D, ³He, and ⁷Li versus cosmic time and temperature. Weak freeze-out sets the n/p ratio; a neutron-decay gap follows; then the deuterium bottleneck delays fusion until T is low enough for D to survive photodissociation, after which ⁴He locks most neutrons into Y_p ≈ 0.245. Sliders η₁₀ (baryon-to-photon ratio ×10¹⁰) and N_eff (effective neutrino species) shift the classic abundance curves, with observed bands overlaid. The network is pedagogical: incomplete reaction rates, no full nuclear Monte Carlo, and no stellar astration modeling. Vary η₁₀ and N_eff to see how helium and deuterium respond — the classic probe of early-universe expansion and baryon density.

Who it's for: Cosmology and nuclear astrophysics courses covering BBN abundances and early-universe parameters.

Key terms

  • Big Bang nucleosynthesis
  • helium mass fraction
  • deuterium bottleneck
  • baryon-to-photon ratio
  • N_eff
  • weak freeze-out

How it works

Big-Bang nucleosynthesis (BBN) — within the first ~3 minutes the universe forged the bulk of its primordial light elements. Watch the n/p ratio drop from equilibrium toward freeze-out at T_F ≈ 0.8 MeV, decay slowly through neutron β-decay, and then suddenly lock into He-4 when the deuterium bottleneck breaks at T ≈ 0.07 MeV. Trace amounts of D, ³He and ⁷Li survive. The asymptotic abundances depend almost only on the baryon-to-photon ratio η₁₀ (and a little on the effective neutrino number N_eff) — moving the slider walks you across the famous BBN curves and explains why measuring just D/H in distant quasars or Y_p in HII regions pins down Ω_b h² independently of the CMB. Right-edge ticks mark the observationally inferred values; note the long-standing ⁷Li discrepancy.

Key equations

T(t) = (1.32 s·MeV²)^{1/2} / √t (radiation era)
(n/p)_eq = exp(−Q/T), Q = m_n − m_p = 1.293 MeV
Y_p ≈ 2(n/p) / (1 + n/p) at BBN onset
D/H ∝ η₁₀^{−1.6}, ³He/H ∝ η₁₀^{−0.6}, ⁷Li/H : U-shape vs η

Frequently asked questions

Why is most helium already set by BBN?
Once deuterium forms past the bottleneck, nearly all available neutrons are quickly assembled into ⁴He. That yields Y_p ≈ 0.25 with only weak dependence on later stellar processing for the primordial floor. Stars make more metals and some helium, but the bulk cosmic ⁴He fraction is a Big Bang relic.
What does increasing η₁₀ do to deuterium?
Higher baryon density means more efficient burning of D into heavier nuclei, so residual deuterium falls steeply with η₁₀. Helium Y_p rises only mildly. That is why measured primordial D/H is a sharp baryometer. A misconception is that all light elements scale the same way with η — they do not.
How does N_eff change the curves?
Extra radiation speeds early expansion, freezing a higher n/p ratio and typically raising Y_p while shifting freeze-out timing. The simulator’s N_eff slider shows that expansion-rate lever next to η₁₀. Full neutrino decoupling physics is simplified; treat the curves as standard pedagogical BBN trends with observed values marked for comparison.