🔭 astrophysics

Big Bang Nucleosynthesis Initial Conditions: Revisiting Wagoner et al. (1967)

This paper corrects a previously overlooked error in the temperature-time relationship for the early universe found in the classic 1967 Wagoner et al. study, thereby providing the accurate initial conditions necessary for Big Bang Nucleosynthesis calculations.

Original authors: Charlie Sharpe, Geraint F. Lewis, Luke A. Barnes

Published 2026-02-20

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Charlie Sharpe, Geraint F. Lewis, Luke A. Barnes

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 universe as a giant, expanding balloon. In the very first moments after it popped into existence (the Big Bang), it was incredibly hot and dense. Scientists have spent decades trying to figure out exactly what happened in those first few minutes, specifically how the universe cooked up the first ingredients: hydrogen, helium, and a little bit of lithium. This process is called Big Bang Nucleosynthesis (BBN).

To understand this "cosmic cooking," scientists need a recipe. But a recipe isn't just about the ingredients; it's also about the timing. You need to know exactly how hot the oven is at every second of the cooking process. If you get the temperature wrong, your cake burns or stays raw.

The Classic Recipe Book

In 1967, a famous physicist named Wagoner and his team wrote a "classic recipe book" for this cosmic cooking. Their paper is still used today as a foundation for understanding how the early universe worked. They provided a formula to tell us: "If the universe is this hot, then this many seconds have passed since the Big Bang."

The Typo That Got Lost in Translation

Here is the problem: The authors of this new paper (Sharpe, Lewis, and Barnes) discovered that Wagoner's 1967 recipe had a mathematical typo.

Think of it like a cooking instruction that says: "Set the oven to 400 degrees, but make sure the timer is set to 'square root of 400' seconds."

The Good News: The number they wrote down (10.4) was actually correct. If you just used that number, your cake would turn out fine.
The Bad News: The formula they wrote to get that number was broken. It was like writing the instruction as "Temperature = (Gravity × Light Speed) divided by Time." If you tried to plug in real numbers into that broken formula, the units wouldn't match up (it would be like trying to measure your height in "seconds"). It was a dimensional disaster.

Why Does This Matter?

You might think, "Well, if the final number was right, who cares about the broken math?"

The issue is that this broken formula didn't just sit in a dusty library. It got copied into modern computer programs that scientists use to simulate the universe.

Program A (NUC1231): This program saw the broken formula, but the programmers were smart. They ignored the broken math part and just hard-coded the correct number (10.4) into the computer. So, this program is actually safe.
Program B (AlterBBN): This program saw the broken formula and tried to follow the instructions literally. It tried to calculate the starting time using the broken math.

The "Cosmic Stopwatch" Analogy

Imagine you are trying to start a race.

The Correct Way: You look at the clock, see it's 10 seconds, and say, "Go!"
The Broken Way: You look at a sign that says "Start time = 10 seconds," but the sign also has a confusing, wrong equation written underneath it.
- If you are a careful runner (like Program A), you ignore the confusing equation and just start at 10 seconds.
- If you are a literal runner (like Program B), you try to solve the confusing equation. Because the equation is broken, you might start the race at the wrong time, or in the wrong place.

In the case of the universe, starting the "race" (the simulation) at the wrong time means you are calculating the temperature of the universe at the wrong moment. While the authors of this paper found that Program B was lucky enough to start so early that it didn't ruin the final result, it's still a dangerous mistake to leave in the code.

The Fix

The authors of this paper have done a "correction notice." They have:

Identified the error: They showed exactly where the 1967 math went wrong (the units didn't match).
Provided the right formula: They gave the correct mathematical expression that actually works with the laws of physics.
Warned the community: They told the people writing the computer codes (like the creators of AlterBBN) to fix their software so they aren't using the broken math anymore.

The Bottom Line

This paper is a bit of "cosmic housekeeping." It's like finding a typo in a famous textbook that everyone has been using for 50 years. Even though the typo didn't change the final answer in this specific case, it's important to fix the math so that future scientists don't get confused or make mistakes when they try to solve even harder problems. It ensures that our understanding of the universe's "first meal" is built on solid, correct foundations.

1. Problem Statement

The paper addresses a long-standing, subtle error in the foundational literature of Big Bang Nucleosynthesis (BBN). Specifically, the authors identify an incorrect analytical expression for the relationship between the temperature of the early universe ( $T$ ) and time ( $t$ ) in the high-temperature limit ( $T \gtrsim 10^{10}$ K).

This error originates from the classic 1967 paper by Wagoner, Fowler, and Hoyle. While the numerical coefficient provided in their work happens to be correct, the analytical formula they derived is dimensionally inconsistent. Despite this, the incorrect formula has been propagated into subsequent key publications and, more critically, into modern numerical codes used to simulate primordial element formation. This creates a risk of theoretical inconsistency in BBN calculations, even if some codes inadvertently bypass the error through hardcoded numerical values.

2. Methodology

The authors employed a rigorous theoretical re-derivation and code audit approach:

Re-derivation of Dynamics: The authors revisited the derivation steps taken by Wagoner et al. (1967). They started with the Friedmann–Lemaître–Robertson–Walker (FLRW) metric and the work-energy equation.
- They accounted for the scaling of baryon mass density ( $\rho_b \propto R^{-3}$ ) and neutrino energy density ( $\rho_\nu \propto R^{-4}$ ).
- They formulated the first-order differential equation for the scale factor $R(T)$ and combined it with the Friedmann equation to derive the differential equation for temperature evolution $dT/dt$.
High-Temperature Limit Analysis: The analysis focused on the regime where $T \gtrsim 10^{10}$ $T ≳ 1 0^{10}$ K. In this limit, electrons are highly relativistic. The authors carefully calculated the total energy density ( $\rho$ $ρ$ ) and the effective energy density of electrons and radiation ( $\rho_1$ $ρ_{1}$ ), accounting for the statistical weights (spin states and species count) of photons, electrons, and neutrinos.
- They established that $\rho = \frac{9}{2}\rho_\gamma$ and $\rho_1 = \frac{11}{4}\rho_\gamma$ .
Dimensional Analysis: The authors performed a dimensional check on Wagoner's final equation (Equation A15 in the original paper) using CGS units to verify physical consistency.
Code Audit: The authors examined the source code of two prominent BBN simulation packages:
1. NUC1231 (Kawano 1992): A Fortran implementation.
2. AlterBBN2 (Arbey et al. 2020): A C implementation.
  They traced how these codes implemented the initial conditions derived from Wagoner et al. (1967).

3. Key Contributions

Identification of the Error: The primary contribution is pinpointing that while Wagoner et al. (1967) correctly calculated the numerical factor ($10.4$), their analytical expression was flawed.
- Wagoner's Incorrect Formula: $T_9 = (12\pi G a_r c^{-2})^{1/4} t^{-1/2}$
- The Flaw: In CGS units, the term inside the parenthesis has units of $s^{-1}K^{-1}$ . Consequently, the entire expression yields units of $s^{-1/2}K^{-1/2}$ , which is dimensionally inconsistent with temperature ( $K$ ).
Correction of the Formula: The authors derived the correct analytical expression by directly solving the differential equation for the high-temperature limit.
- Corrected Formula: $T_9 = (48\pi a_r c^{-2} G)^{-1/4} t^{-1/2}$
- This expression is dimensionally consistent (yielding units of Kelvin) and retains the correct numerical factor of $10.4$ when $t$ is in seconds.
Resolution of Historical Discrepancy: The paper clarifies why Wagoner's result differed from Alpher et al. (1953) (who found a factor of 15.2). Wagoner included electrons and neutrinos in the energy density, whereas Alpher did not. The authors confirm that Wagoner's numerical inclusion of these particles was correct, but the algebraic derivation of the constant was botched.

4. Results

Analytical Correction: The paper provides the mathematically rigorous, dimensionally correct expression for $T(t)$ in the early universe:
$T_9 = 10.4 t^{-1/2}$
where $T_9$ is temperature in $10^9$ K and $t$ is time in seconds. The authors emphasize that the numerical coefficient $10.4$ is correct, but the physical constant grouping in Wagoner's original text is wrong.
Impact on Simulation Codes:
- NUC1231: The code reproduces the incorrect analytical form in its documentation (Equation D.6) but, upon inspection of the source code, uses a hardcoded numerical constant (const1=0.09615). Thus, the simulation results remain numerically accurate despite the theoretical error in the documentation.
- AlterBBN2: The code cites the incorrect formula and implements it directly in the source code (sqrt(12.*pi*G*sigma_SB)/pow(Ti,2.)). However, the authors note that because AlterBBN uses natural (energy) units, the resulting starting time is extremely early (deep within the high-temperature limit). Consequently, the error does not significantly impact the final integration results for this specific code, though the theoretical formulation remains flawed.

5. Significance

Restoring Theoretical Rigor: Big Bang Nucleosynthesis is a pillar of modern cosmology. Ensuring that the initial conditions for these calculations are theoretically sound is essential for the integrity of the field. Correcting a 50-year-old error prevents the propagation of dimensional inconsistencies in future theoretical work.
Code Maintenance and Updates: The paper serves as a critical alert to developers of cosmological simulation software. It highlights that while some codes may be "lucky" enough to produce correct numerical results due to hardcoded constants, the underlying theoretical implementation must be corrected to ensure robustness and transparency.
Immediate Action: The authors have communicated their findings to the authors of AlterBBN, who have agreed to correct the initial condition implementation in a forthcoming update. This ensures that future iterations of widely used BBN tools will rely on the correct physical derivation.

In summary, Sharpe, Lewis, and Barnes have performed a necessary "sanity check" on the foundations of BBN, correcting a dimensional error in a seminal paper that, while not altering the numerical outcomes of specific simulations due to coincidental numerical fixes, was theoretically unsound and required correction for the sake of scientific accuracy.