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

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

Auteurs: Charlie Sharpe, Geraint F. Lewis en Luke A. Barnes
Datum: 20 februari 2026 (conceptversie)

1. Het Probleem

Big Bang Nucleosynthese (BBN) is een hoeksteen van de moderne kosmologie die de abundantie van oorspronkelijke elementen in het heelal verklaart. De berekening van de output van BBN vereist de integratie van gekoppelde differentiaalvergelijkingen, waarbij de initiële voorwaarden (startcondities) cruciaal zijn.

De auteurs wijzen erop dat een klassiek artikel, Wagoner et al. (1967), een fundamentele fout bevat in de afleiding van de relatie tussen de temperatuur ($T$) van het vroege heelal en de tijd ($t$) in het hoge-temperatuurlimiet ($T \gtrsim 10^{10}$ K). Hoewel de numerieke factor in hun bekende vergelijking correct lijkt te zijn, is de analytische uitdrukking dimensionaal incorrect. Deze fout is decennialang onopgemerkt gebleven en is doorgegeven in andere belangrijke publicaties en numerieke codes voor BBN-simulaties.

2. Methodologie

De auteurs analyseren de afleiding van Wagoner et al. (1967) stap voor stap:

• Analyse van de oorspronkelijke afleiding: Wagoner et al. leiden een differentiaalvergelijking af voor de expansie van het heelal ($R$) als functie van temperatuur, gebaseerd op de Friedmann-Lemaître-Robertson-Walker-metric. Ze combineren dit met de Friedmann-vergelijking om een relatie voor $T(t)$ te vinden.
• Identificatie van de fout: In het hoge-temperatuurlimiet zijn elektronen en neutrino's ultrarelativistisch. Wagoner et al. gebruiken de energiedichtheden correct ($\rho = \frac{9}{2}\rho_\gamma$ en $\rho_1 = \frac{11}{4}\rho_\gamma$), maar hun uiteindelijke formule (hun vergelijking A15) luidt:
  $$T_9 = (12\pi G a_r c^{-2})^{1/4} t^{-1/2} = 10.4 t^{-1/2}$$
  De auteurs tonen aan dat de eerste gelijkheid dimensionaal incorrect is. In cgs-eenheden heeft de centrale uitdrukking de eenheid $s^{-1}K^{-1}$, terwijl $T_9$ (temperatuur in eenheden van $10^9$ K) de eenheid Kelvin moet hebben.
• Correcte afleiding: De auteurs lossen de onderliggende differentiaalvergelijking (vergelijking 2 in het artikel) direct op, rekening houdend met de juiste energiedichtheden van straling, elektronen en neutrino's. Ze vergelijken dit ook met eerdere werken zoals Alpher et al. (1953), die een andere numerieke factor vonden omdat zij geen rekening hielden met elektronen en neutrino's.

3. Belangrijkste Bijdragen en Resultaten

• Correcte Formule: De auteurs stellen de correcte dimensionale uitdrukking voor de temperatuur-tijd relatie vast:
  $$T_9 = (48\pi a_r c^{-2} G)^{-1/4} t^{-1/2} = 10.4 t^{-1/2}$$
  Hoewel de numerieke factor ($10.4$) toevallig overeenkomt met die van Wagoner et al., is de analytische vorm in de haakjes nu dimensionaal correct (eenheden van Kelvin).
• Identificatie van Verspreide Fouten: Het artikel toont aan dat de foutieve analytische uitdrukking van Wagoner et al. is overgenomen in moderne BBN-software:
  • NUC1231 (Kawano, 1992): Een Fortran-implementatie die de formule van Wagoner direct overneemt als startconditie, hoewel de broncode de numerieke factor correct toepast.
  • AlterBBN2 (Arbey et al., 2020): Een C-implementatie die de foutieve formule in de broncode implementeert voor het bepalen van de starttijd. De auteurs merken op dat dit in de natuurlijke eenheden van AlterBBN waarschijnlijk geen invloed heeft op de uiteindelijke integratie, maar dat de code toch gecorrigeerd zal worden.

4. Betekenis en Conclusie

De significance van dit artikel ligt in het herstellen van een fundamentele theoretische fout in een klassiek werk dat de basis vormt voor BBN-berekeningen.

• Theoretische Zuiverheid: Het corrigeert een dimensionale inconsistentie die in een standaardreferentie (Wagoner et al. 1967) aanwezig was.
• Impact op Simulaties: Hoewel de numerieke factor in de meeste gevallen toevallig correct bleek te zijn (waardoor de resultaten van bestaande simulaties niet per se onjuist waren), is het essentieel dat de onderliggende formules correct zijn voor toekomstige precisieberekeningen en voor het voorkomen van fouten in nieuwe codes.
• Actie: De auteurs hebben de auteurs van AlterBBN gecontacteerd, en een correctie van de initiële condities in een toekomstige update van die software is bevestigd.

Kortom, dit onderzoek noteert een "kleine typografische fout" die zich heeft verspreid tot in de kern van de kosmologische modellering en biedt de wiskundig correcte basis voor het berekenen van de vorming van oorspronkelijke elementen.