The Early Universe, Big Bang Nucleosynthesis, CMB & Astrobiology

Extrapolating cosmic expansion backward in time reveals that the early universe was an ultra-dense, ultra-hot relativistic plasma of elementary particles in thermal equilibrium.

Thermal Physics of the Radiation Era

For relativistic particles ($k T \gg m c^2$), the total energy density $\rho_r c^2$ is governed by the Stefan-Boltzmann law generalized to multiple particle species:

where $g_*(T) = \sum_{\text{bosons}} g_b \left(\frac{T_b}{T}\right)^4 + \frac{7}{8} \sum_{\text{fermions}} g_f \left(\frac{T_f}{T}\right)^4$ counts the effective relativistic degrees of freedom (the factor $7/8$ arises from Fermi-Dirac vs Bose-Einstein integrals). In the radiation era where $H^2 = \frac{8\pi G}{3}\rho_r$, the time-temperature relationship is:

Chronological Epochs of the Early Universe

Big Bang Nucleosynthesis (BBN) represents the earliest testable empirical milestone in cosmology, predicting the primordial elemental abundances synthesized during the first twenty minutes of the Universe.

1. Neutron-to-Proton Freeze-Out ($t \approx 1\text{ s}$)

At $T > 1\text{ MeV}$ ($t < 1\text{ s}$), neutrons and protons are kept in thermal chemical equilibrium via rapid weak interactions mediated by electron neutrinos:

The neutron-to-proton ratio is governed by the Boltzmann factor with neutron-proton mass difference $\Delta m = m_n - m_p = 1.293\text{ MeV}$:

The weak reaction rate scales as $\Gamma_{\text{weak}} = n_e \langle \sigma v \rangle \propto G_F^2 T^5$ (where $G_F$ is the Fermi coupling constant). Meanwhile, the Hubble expansion rate in the radiation era scales as $H \propto \sqrt{g_*} T^2$.

As the universe cools, $\Gamma_{\text{weak}}$ plummets much faster than $H$. Freeze-out occurs when the reaction rate drops below the cosmic expansion rate:

At freeze-out, the neutron-to-proton ratio is locked at:

2. Free Neutron Decay & The Deuterium Bottleneck

Nuclear fusion cannot proceed directly to ${}^4\text{He}$ because four-body collisions ($2p + 2n$) are impossibly rare. Synthesis must proceed through the two-body stepping stone of deuterium:

Although the binding energy of deuterium is $B_d = 2.225\text{ MeV}$, the immense ratio of photons to baryons ($\eta^{-1} = n_\gamma / n_b \approx 1.6 \times 10^9$) means that high-energy photons in the Planck tail constantly photo-dissociate deuterium back into free protons and neutrons! This delay is the deuterium bottleneck.

During this delay from $t = 1\text{ s}$ to $t \approx 300\text{ s}$ ($T$ dropping to $\approx 0.08\text{ MeV}$), free neutrons undergo radioactive beta decay ($n \to p + e^- + \bar{\nu}_e$) with mean lifetime $\tau_n = 879.4\text{ seconds}$:

3. The Primordial Helium-4 Mass Fraction ($Y_p$)

Once deuterium stabilizes at $T \sim 80\text{ keV}$, nuclear fusion cascades rapidly through two-body reactions:

Because ${}^4\text{He}$ has a colossal binding energy ($28.3\text{ MeV} = 7.07\text{ MeV/nucleon}$), virtually all available neutrons are rapidly locked into ${}^4\text{He}$ nuclei. Since each ${}^4\text{He}$ nucleus requires 2 neutrons and 2 protons, the primordial helium mass fraction $Y_p$ is:

Substituting $(n/p) \approx 1/7$:

This matches observations of pristine, metal-poor extragalactic H II regions ($Y_p = 0.245 \pm 0.003$). Because no stable nuclei exist with mass numbers $A=5$ or $A=8$ (neither ${}^5\text{He}, {}^5\text{Li}$ nor ${}^8\text{Be}$ are stable), and the universe cools too rapidly to ignite the triple-alpha process, BBN ceases after 20 minutes, leaving $75\%$ hydrogen, $25\%$ helium, and trace amounts of deuterium ($D/H \sim 2.5 \times 10^{-5}$) and lithium (${}^7\text{Li}/H \sim 1.6 \times 10^{-10}$).

For the first 380,000 years, the universe was an opaque, foggy plasma where photons were tightly coupled to free electrons via Thomson scattering. Photons could not travel freely, possessing a mean free path $\ell_{\text{mfp}} = \frac{1}{n_e \sigma_T}$ of only a few light-years.

Cosmological Hydrogen Recombination

As the universe expanded and cooled below $T \sim 3{,}000\text{ K}$, free electrons combined with protons to form neutral hydrogen: $p + e^- \rightleftharpoons \text{H} + \gamma$.

The ionization fraction $X_e = n_e / n_b = n_e / (n_p + n_{\text{H}})$ is governed by the cosmological Saha equation:

Because the photon-to-baryon ratio is so high ($\eta \approx 6 \times 10^{-10}$), recombination does not occur at $k T = 13.6\text{ eV}$ ($T \sim 150{,}000\text{ K}$), but is delayed until $k T \approx 0.3\text{ eV}$ ($T_{\text{rec}} \approx 3{,}000\text{ K}$), corresponding to redshift $z_{\text{rec}} \approx 1100$.

Photon Decoupling & The Surface of Last Scattering

As $X_e \to 0$, the optical depth for Thomson scattering drops below unity:

Photons decoupled from matter and streamed freely across the transparent universe. When we look out into deep space, we look back in time to this spherical boundary—the Surface of Last Scattering (SLS).

Discovery & The Blackbody Nature of the CMB

In 1965, Arno Penzias and Robert Wilson discovered an isotropic, unpolarized microwave noise at Bell Labs with temperature $\approx 3.5\text{ K}$ (Nobel Prize 1978), matching predictions by Ralph Alpher, Robert Herman, and George Gamow (1948). In 1990, the FIRAS instrument on NASA's COBE satellite measured the CMB spectrum, establishing it as the most perfect blackbody spectrum ever observed in nature, with modern temperature:

The photon number density today is $n_\gamma = \frac{2.404}{\pi^2}\left(\frac{kT}{\hbar c}\right)^3 \approx 411\text{ photons cm}^{-3}$.

While the CMB is isotropic to one part in $10^5$, minute temperature fluctuations $\Delta T(\theta, \phi) / T \sim 10^{-5}$ contain a pristine snapshot of the primordial density perturbations that seeded all cosmic structures.

Multipole Expansion & Angular Power Spectrum

Temperature fluctuations across the celestial sphere are expanded into spherical harmonics:

The variance of the expansion coefficients defines the angular power spectrum ($C_\ell$):

Multipole $\ell$ corresponds to angular separation $\theta \approx 180^\circ / \ell$. The power per logarithmic interval is plotted as $\mathcal{D}_\ell = \frac{\ell(\ell+1)}{2\pi} C_\ell$ in $\mu\text{K}^2$.

Physics of the Acoustic Peaks

Prior to recombination, dark matter gravity pulled gas into gravitational potential wells, while photon radiation pressure pushed back. This competition set up longitudinal standing sound waves in the relativistic baryon-photon fluid—acoustic oscillations.

• First Acoustic Peak ($\ell \approx 220, \theta \approx 0.8^\circ$): Corresponds to perturbation modes that had time to undergo exactly one maximum gravitational compression before decoupling. Its physical scale is the sound horizon $r_s \approx 147\text{ Mpc}$. Its angular location on the sky measures the geometry of space: $$\theta = \frac{r_s}{d_A} \implies \ell_{\text{peak}} \approx \frac{220}{\sqrt{1 - \Omega_k}}$$ Measuring $\ell = 220$ proves that our Universe is spatially flat: $\Omega_k = 0.0007 \pm 0.0019$!
• Second Peak ($\ell \approx 540$): Rarefaction mode (maximum decompression). The ratio of the first-to-second peak height measures the total baryonic mass density $\Omega_b h^2 = 0.02237 \pm 0.00015$ (baryons add gravitational inertia, enhancing odd compression peaks relative to even rarefaction peaks).
• Third Peak ($\ell \approx 800$): Second compression mode, measuring the cold dark matter density $\Omega_c h^2 = 0.1200 \pm 0.0012$.

The Planck 2018 mission established the cosmological concordance parameters to sub-percent precision: $\Omega_b = 0.049$, $\Omega_c = 0.266$, $\Omega_\Lambda = 0.685$, and age $t_0 = 13.787 \pm 0.020\text{ Gyr}$.

Despite the triumphs of the standard Big Bang model, three fundamental cosmological paradoxes remained unexplained.

Three Puzzles of Classical Big Bang Theory

1. The Horizon Problem: The particle horizon at recombination subtended an angle of only $\theta_{\text{hor}} \approx 1^\circ$ on the CMB sky. The sky contains $\approx 40{,}000$ causally disconnected regions that could never have exchanged light signals since the Big Bang. Yet their temperatures are identical to within one part in $10^5$! How did causally disconnected patches establish thermal equilibrium?
2. The Flatness Problem: The Friedmann equation gives $|\Omega(t) - 1| = \frac{|k| c^2}{a^2 H^2}$. In a matter- or radiation-dominated universe, $a^2 H^2 \propto a^{-1}$ or $a^{-2}$ decreases with time, so $|\Omega - 1|$ grows rapidly. To observe $|\Omega - 1| < 0.002$ today requires the early universe to have been fine-tuned at the Planck epoch to $|\Omega(t_P) - 1| < 10^{-60}$!
3. The Magnetic Monopole Problem: Grand Unified Theories (GUTs) predict that spontaneous symmetry breaking at $T_{\text{GUT}} \sim 10^{16}\text{ GeV}$ would produce copious topological defects, including supermassive magnetic monopoles ($m \sim 10^{16}\text{ GeV}$), with density exceeding critical density by $10^{12}$! None have ever been detected.

Guth's Cosmic Inflation Hypothesis (1981)

Alan Guth proposed that during the GUT epoch ($t \sim 10^{-36}\text{ s}$), a scalar field $\phi$ (the inflaton) was displaced from its potential minimum into a 'false vacuum' state with potential energy $V(\phi) \approx \text{const}$. The universe underwent a transient de Sitter exponential expansion:

Over a duration of $\Delta t \sim 10^{-32}\text{ s}$, the universe expanded by a factor of $e^N$ with $N \ge 60$ e-folds:

Inflation effortlessly resolves all three paradoxes:

• Horizon Problem Resolved: The entire observable universe was expanded from a sub-microscopic patch ($< 10^{-28}\text{ cm}$) that had ample time to achieve causal thermal equilibrium before inflation stretched it beyond the horizon.
• Flatness Problem Resolved: During inflation, $a^2 H^2 \propto e^{2 H t}$ grows exponentially, driving $|\Omega - 1| \propto e^{-2N} \to 0$. Just as blowing up a balloon flattens its local surface curvature, inflation drives space to near-perfect flatness ($k \to 0$).
• Monopole Problem Resolved: The primordial monopole density is diluted to less than one monopole per observable universe.
• Origin of Structure: Subatomic quantum vacuum fluctuations in the inflaton field $\delta\phi$ were stretched to macroscopic astronomical scales, generating the scale-invariant Gaussian density perturbations ($n_s \approx 0.965$) that collapsed into all modern galaxies.

The geometry and equation of state $w$ dictate both the observational limits of our view and the ultimate fate of the cosmos.

Cosmological Horizons

• The Particle Horizon ($d_{\text{part}}$): The maximum proper distance from which an observer at time $t$ could have received a light signal emitted at the beginning of time ($t=0$):
  $$d_{\text{part}}(t) = a(t) \int_0^t \frac{c\,dt'}{a(t')} = c \int_0^z \frac{dz'}{H(z')}$$
  Today, $d_{\text{part}}(t_0) \approx 14.3\text{ Gpc} \approx 46.5\text{ billion light-years}$, bounding the observable universe.
• The Event Horizon ($d_{\text{event}}$): The maximum proper distance from which a light signal emitted now can ever reach an observer in the infinite future ($t \to \infty$):
  $$d_{\text{event}}(t) = a(t) \int_t^\infty \frac{c\,dt'}{a(t')}$$
  In an accelerating universe dominated by $\Lambda$, $d_{\text{event}}$ approaches a finite constant: $d_{\text{event}} \approx c / H_0 \sqrt{\Omega_\Lambda} \approx 5.2\text{ Gpc} \approx 17\text{ Gly}$. Any galaxy currently situated beyond $17\text{ Gly}$ is forever lost to our future; light it emits today will never reach the Milky Way!

The Ultimate Fate of the Universe

1. The Big Freeze / Heat Death (Standard $\Lambda\text{CDM}$, $w = -1$): Dark energy accelerates space forever. In $\sim 100\text{ Gyr}$, all galaxies outside the Local Group are swept beyond the event horizon, leaving our merged galaxy (Milkomeda) isolated in an empty universe. In $\sim 10^{14}\text{ yr}$, star formation ceases as gas is exhausted. In $\sim 10^{40}\text{ yr}$, protons decay. In $\sim 10^{100}\text{ yr}$, supermassive black holes evaporate via Hawking radiation, leaving a cold, dilute sea of photons, electrons, and positrons at thermodynamic maximum entropy.
2. The Big Rip (Phantom Dark Energy, $w < -1$): If dark energy density grows with time ($w < -1$), the scale factor diverges to infinity at a finite time $t_{\text{rip}} = t_0 + \frac{2}{3|1+w|H_0 \sqrt{\Omega_\Lambda}}$. As the phantom energy density surges, it progressively tears apart galaxy clusters, galaxies, planetary orbits, planets, atoms, and finally spacetime itself.
3. The Big Crunch (Closed Decelerating Universe, $\Omega > 1, \Lambda \le 0$): Self-gravity halts expansion and reverses it into cosmic collapse, ending in a catastrophic high-density singularity.

Astrobiology investigates the origin, evolution, distribution, and future of life in the cosmos, bridging stellar astrophysics, planetary science, and cosmology.

1. The Circumstellar Habitable Zone (CHZ)

The Habitable Zone (colloquially the 'Goldilocks Zone') is the circumstellar orbital shell where an Earth-like planet with an atmosphere can maintain liquid water on its surface. The inner boundary is set by the runaway greenhouse effect (solar flux evaporates oceans, water vapor saturates the stratosphere and is lost via UV photolysis), while the outer boundary is set by the maximum greenhouse / runaway glaciation (CO$_2$ condensation and ice-albedo freeze):

For the Sun ($1 L_\odot$), the conservative habitable zone spans $\approx 0.95 - 1.67\text{ AU}$. For an M-dwarf ($L \sim 10^{-3} L_\odot$), the habitable zone lies exceedingly close ($r \sim 0.03 - 0.1\text{ AU}$), where exoplanets become tidally locked and exposed to intense stellar magnetic flares.

2. The Drake Equation (1961)

Frank Drake formulated a probabilistic framework to estimate the number $N$ of active, communicative technological civilizations in our Milky Way galaxy:

where:

• $R_*$: Average rate of star formation in the Milky Way ($\sim 1-2 M_\odot\text{ yr}^{-1}$).
• $f_p$: Fraction of stars that have planetary systems ($f_p \approx 1.0$, confirmed by Kepler).
• $n_e$: Number of planets per system in the habitable zone ($n_e \approx 0.2-0.5$).
• $f_l$: Fraction of habitable planets where life actually arises ($0 < f_l \le 1$).
• $f_i$: Fraction of life-bearing planets that develop intelligent civilizations ($0 < f_i \le 1$).
• $f_c$: Fraction of civilizations that develop detectable radio communication technology ($0 < f_c \le 1$).
• $L$: The average longevity of such a communicative civilization in years.

Notice that the product of the astrophysical terms is well constrained: $R_* f_p n_e \sim 0.5 - 1.0\text{ yr}^{-1}$. Thus, the equation simplifies to $N \approx (f_l f_i f_c) \times L$. If civilizations endure for $L \sim 10{,}000\text{ years}$, $N$ could be in the hundreds or thousands; if $L \sim 100\text{ years}$ (due to nuclear, environmental, or technological self-destruction), $N \approx 1$—we may be entirely alone.

3. The Fermi Paradox & The Great Filter

During a 1950 lunch at Los Alamos, Enrico Fermi famously asked: "Where is everybody?"

Given that our Galaxy is $\sim 13\text{ billion years}$ old and crossing it at $0.01c$ requires only $\sim 10\text{ million years}$ (a cosmic blink of an eye), an intelligent civilization should have colonized the entire Galaxy long ago. Proposed resolutions include:

• The Great Filter (Hanson 1998): A near-insurmountable evolutionary barrier exists somewhere between prebiotic chemistry and interstellar colonization. If the filter lies behind us (abiogenesis or eukaryotic complexity is exceedingly rare), humanity is unique. If the filter lies ahead of us (technological civilizations inevitably annihilate themselves), humanity faces imminent peril.
• The Rare Earth Hypothesis: Complex multicellular life requires an exceptionally rare confluence of planetary factors (plate tectonics, large stabilizing Moon, Jupiter shield against comets, quiet magnetic star).
• The Zoo Hypothesis / Technological Singularity: Extraterrestrial intelligence deliberately avoids contacting undeveloped planets, or transitions into post-biological digital matrices unobservable by radio telescopes.