Part II: Compact Objects | Chapter 4

Gamma-Ray Bursts

The most luminous explosions in the Universe: the fireball model, relativistic beaming, afterglow emission, and the short vs long GRB dichotomy

Overview

Gamma-ray bursts (GRBs) are brief, intense flashes of gamma-ray radiation from cosmological distances. With isotropic-equivalent luminosities reaching $10^{53}$ erg/s and total energies of $10^{51}\text{--}10^{54}$ erg (isotropic equivalent), they are the most energetic electromagnetic events in the Universe. GRBs were first detected by the Vela nuclear test detection satellites in 1967 and declassified in 1973. The cosmological origin was established by the BeppoSAX satellite in 1997 through the first afterglow detection and redshift measurement.

In this chapter we derive the compactness problem and its resolution through relativistic expansion, develop the fireball model, analyze synchrotron afterglow emission from the external shock, and explore the physical distinction between short and long GRBs.

1. The Compactness Problem and Relativistic Motion

The key observational constraints on GRBs are their enormous luminosities and rapid variability (timescales as short as milliseconds). Together, these create the "compactness problem."

1.1 The Compactness Argument

A source of luminosity $L$ varying on timescale $\delta t$ has a size $R \lesssim c\,\delta t$. For a typical GRB with $L_\gamma \sim 10^{52}$ erg/s and $\delta t \sim 10$ ms:

$$R \lesssim c\,\delta t \sim 3 \times 10^8\;\text{cm}$$

The optical depth to pair production ($\gamma\gamma \to e^+e^-$) for photons above the threshold energy $m_e c^2 = 511$ keV is:

$$\tau_{\gamma\gamma} \sim \frac{\sigma_T L_\gamma}{4\pi R m_e c^3} \sim \frac{\sigma_T L_\gamma}{4\pi c^2 \delta t\, m_e c^3} \sim 10^{15}$$

This implies the source should be completely opaque to gamma rays! Yet we observe non-thermal spectra extending to high energies, meaning the photons escape freely. The resolution is relativistic bulk motion.

1.2 Resolution via Bulk Lorentz Factor

If the emitting material moves toward us with Lorentz factor $\Gamma \gg 1$, two effects reduce $\tau_{\gamma\gamma}$:

(1) The source size in the comoving frame is $R' \sim \Gamma^2 c\,\delta t$ (from relativistic time compression).

(2) The photon energies in the comoving frame are lower by a factor $\Gamma$, reducing the fraction of photons above pair threshold.

The net effect is:

$$\boxed{\tau_{\gamma\gamma} \propto \Gamma^{-(4+2\beta)}}$$

where $\beta$ is the high-energy photon index. Requiring $\tau_{\gamma\gamma} < 1$gives minimum Lorentz factors of $\Gamma \gtrsim 100\text{--}1000$. This makes GRBs the most relativistic macroscopic objects known.

2. The Fireball Model

The standard fireball model, developed by Rees and Meszaros (1992, 1994), describes how a compact energy release produces the observed GRB emission.

2.1 Initial Fireball Phase

Energy $E \sim 10^{51}$ erg deposited in a volume $R_0 \sim 10^7$ cm produces a radiation-dominated fireball with temperature:

$$T_0 = \left(\frac{E}{aR_0^3}\right)^{1/4} \sim 10^{10}\;\text{K} \sim 1\;\text{MeV}$$

This is well above the pair-production threshold, so the fireball is initially opaque and loaded with $e^+e^-$ pairs. The radiation pressure drives rapid expansion with terminal Lorentz factor:

$$\boxed{\Gamma_\infty = \frac{E}{M_b c^2} = \eta}$$

where $M_b$ is the baryonic mass entrained in the outflow and $\eta$ is the baryon loading parameter. For $\eta \gtrsim 100$, the fireball achieves the required ultra-relativistic velocities.

2.2 Internal Shocks

The central engine (likely an accreting black hole or magnetar) produces an unsteady outflow with variable Lorentz factor. Faster shells catch up with slower shells at the internal shock radius:

$$R_{\text{IS}} \sim 2\Gamma^2 c\,\delta t \sim 3 \times 10^{13}\left(\frac{\Gamma}{300}\right)^2\left(\frac{\delta t}{10\,\text{ms}}\right)\;\text{cm}$$

The collisions dissipate kinetic energy and accelerate particles, producing the observed prompt gamma-ray emission through synchrotron radiation and inverse Compton scattering. The efficiency of internal shocks is typically $\sim 10\text{--}20\%$.

3. Relativistic Beaming

A key consequence of ultra-relativistic motion is that radiation is concentrated into a narrow cone of opening angle $\sim 1/\Gamma$.

3.1 Beaming and the Jet Break

If the GRB outflow is collimated into a jet of half-opening angle $\theta_j$, the true energy release is reduced from the isotropic equivalent by the beaming factor:

$$\boxed{E_{\text{true}} = f_b\,E_{\text{iso}}, \qquad f_b = 1 - \cos\theta_j \approx \frac{\theta_j^2}{2}}$$

For typical jet angles $\theta_j \sim 5°\text{--}10°$, the beaming correction is $f_b \sim 10^{-2}\text{--}10^{-3}$, reducing the $\sim 10^{53}$ erg isotropic equivalent to $\sim 10^{51}$ erg true energy — remarkably similar to the energy of a supernova.

3.2 Superluminal Motion

A source moving at angle $\theta$ to the line of sight with speed $\beta c$ appears to move across the sky with apparent transverse velocity:

$$\beta_{\text{app}} = \frac{\beta\sin\theta}{1 - \beta\cos\theta}$$

The maximum apparent velocity is $\beta_{\text{app,max}} = \Gamma\beta$, achieved at $\theta = 1/\Gamma$. For $\Gamma = 300$, the apparent transverse velocity can reach $300c$, an effect directly observed in GRB afterglow radio imaging.

4. Afterglow Physics

When the relativistic ejecta decelerate against the external medium, they produce a broadband afterglow visible from radio to X-rays, lasting days to months.

4.1 Deceleration Radius

The ejecta begin to decelerate significantly when they have swept up a mass comparable to $E/(\Gamma^2 c^2)$:

$$R_{\text{dec}} = \left(\frac{3E}{4\pi n m_p c^2 \Gamma^2}\right)^{1/3} \sim 10^{17}\;\text{cm}$$

4.2 Synchrotron Spectrum

Shock-accelerated electrons with a power-law energy distribution $N(\gamma_e) \propto \gamma_e^{-p}$(with $p \approx 2.2\text{--}2.5$) produce synchrotron radiation in the shock-amplified magnetic field. The characteristic synchrotron frequency of an electron with Lorentz factor $\gamma_e$ is:

$$\nu_{\text{syn}} = \frac{3}{4\pi}\frac{eB'}{m_e c}\gamma_e^2\,\Gamma$$

The afterglow spectrum has three characteristic break frequencies: the self-absorption frequency $\nu_a$, the minimum injection frequency $\nu_m$, and the cooling frequency $\nu_c$. Between these breaks, the flux follows power laws with slopes determined by $p$:

$$\boxed{F_\nu \propto \begin{cases} \nu^2 & \nu < \nu_a \\ \nu^{1/3} & \nu_a < \nu < \nu_m \\ \nu^{-(p-1)/2} & \nu_m < \nu < \nu_c \\ \nu^{-p/2} & \nu > \nu_c \end{cases}}$$

5. Short vs Long GRBs: Two Progenitor Channels

The duration distribution of GRBs is bimodal, with a division at approximately 2 seconds.

5.1 Long GRBs: Collapsar Model

Long GRBs ($T_{90} > 2$ s) are associated with the death of massive stars. The collapsar model (Woosley 1993; MacFadyen & Woosley 1999) posits that the iron core of a rapidly rotating massive star collapses to form a black hole. An accretion disk forms from the infalling stellar material, and the accretion energy $\dot{E} \sim \eta \dot{M} c^2$ powers a relativistic jet along the rotation axis. The jet must drill through the stellar envelope (taking $\sim 10$ s), which naturally explains the long duration. Observationally, long GRBs are found exclusively in star-forming galaxies and are sometimes associated with broad-lined Type Ic supernovae (e.g., GRB 030329/SN 2003dh).

5.2 Short GRBs: Compact Binary Mergers

Short GRBs ($T_{90} < 2$ s) originate from the merger of compact binary systems (neutron star–neutron star or neutron star–black hole). The merger timescale is set by the orbital decay due to gravitational wave emission. The merger produces a hypermassive neutron star or black hole with an accretion disk, which launches a short-lived relativistic jet. The definitive confirmation came from GW170817/GRB 170817A: the near-simultaneous detection of gravitational waves from a binary neutron star merger and a short GRB by LIGO/Virgo and Fermi-GBM.

5.3 The Amati Relation

An empirical correlation exists between the peak energy of the prompt emission spectrum $E_{\text{peak}}$ and the isotropic equivalent energy $E_{\text{iso}}$:

$$E_{\text{peak}} \propto E_{\text{iso}}^{0.5}$$

Long GRBs follow this relation closely, while short GRBs are outliers. The physical origin of this correlation remains debated, but it may reflect a common emission mechanism.

Applications

GRBs as Probes of the High-Redshift Universe

GRBs have been detected out to $z = 9.4$ (GRB 090429B), making them the most distant electromagnetic sources known. Their bright afterglows serve as backlights for absorption spectroscopy of the intervening intergalactic medium, probing the reionization epoch, the chemical enrichment history of the early Universe, and the properties of host galaxies at redshifts inaccessible to other techniques.

Kilonova and r-Process Nucleosynthesis

The kilonova AT2017gfo associated with GW170817/GRB 170817A showed that neutron star mergers produce substantial quantities of r-process elements. The red kilonova component (powered by lanthanide-rich ejecta with high opacities) and blue component (lanthanide-poor) together indicate synthesis of $\sim 0.05\,M_\odot$ of heavy elements including gold and platinum.

Very-High-Energy Afterglows

The MAGIC and H.E.S.S. Cherenkov telescope arrays detected TeV gamma-ray emission from GRB afterglows for the first time in 2019 (GRB 190114C and GRB 180720B). The TeV photons are produced by inverse Compton scattering of synchrotron photons by the same shock- accelerated electrons (synchrotron self-Compton, SSC). The detection of the SSC component constrains the magnetic field strength and electron energy distribution in the afterglow shock, resolving long-standing degeneracies in afterglow modeling. These observations also set the strongest constraints on Lorentz invariance violation at energies approaching the Planck scale.

The BOAT: GRB 221009A

GRB 221009A, the "Brightest Of All Time" (BOAT), was detected on 9 October 2022 at a redshift of $z = 0.151$. With an isotropic-equivalent energy of $\sim 10^{55}$ erg, it was the most energetic GRB ever recorded and saturated most gamma-ray detectors. It produced detectable TeV emission and ionized the Earth's ionosphere. Statistical analysis suggests such an event occurs only once every $\sim 10{,}000$ years, making it a truly exceptional opportunity to study the most extreme explosion physics.

Historical Notes

GRBs were accidentally discovered by the Vela satellites, designed to monitor nuclear test ban treaty compliance, in 1967. The discovery was published by Klebesadel, Strong, and Olson in 1973. For 25 years, even the distance scale was unknown — were GRBs Galactic or cosmological? The BATSE instrument on the Compton Gamma Ray Observatory (1991–2000) showed that GRBs are isotropically distributed on the sky with a deficit of faint events, ruling out Galactic disk models. The breakthrough came in 1997 when BeppoSAX detected the first X-ray afterglow (GRB 970228), enabling optical identification and redshift measurement. The collapsar model was confirmed by the GRB 030329/SN 2003dh association. The multi-messenger detection of GW170817/GRB 170817A in 2017 confirmed the compact merger origin of short GRBs and inaugurated the era of multi-messenger astronomy.

The theoretical framework for GRBs was developed through contributions from many physicists. Bohdan Paczynski and James Goodman (1986) proposed the cosmological fireball model. Martin Rees and Peter Meszaros (1992, 1994) developed the internal-external shock framework that remains the standard model. Stan Woosley (1993) proposed the collapsar model for long GRBs, and Neil Gehrels, Chryssa Kouveliotou, and collaborators established the observational framework for afterglow studies. The Swift satellite (launched 2004) with its rapid-slewing X-ray telescope transformed GRB science by enabling prompt afterglow observations for hundreds of events per year.

The field continues to evolve rapidly. The detection of very-high-energy ($> 100$ GeV) afterglow emission by MAGIC and H.E.S.S. has confirmed the synchrotron self-Compton mechanism. The Einstein Probe satellite (launched 2024) is detecting X-ray transients with unprecedented speed and sensitivity. Future missions such as THESEUS and the Gamow Explorer aim to detect GRBs at the highest redshifts ($z > 10$) to probe the earliest star formation in the Universe. The joint detection of gravitational waves and electromagnetic emission from GRBs remains a primary science driver for next-generation multi-messenger facilities. The observation rate of such joint events is expected to increase dramatically with improvements in gravitational wave detector sensitivity during the O5 observing run and beyond.

Central Engine Physics

The central engine of a GRB must produce $\sim 10^{51}$ erg of collimated kinetic energy in $\sim 10\text{--}100$ seconds (long GRBs) or$\sim 0.1\text{--}2$ seconds (short GRBs). Two leading models exist:

Black Hole Accretion

A stellar-mass black hole ($3\text{--}10\,M_\odot$) surrounded by a massive accretion disk ($\sim 0.01\text{--}1\,M_\odot$) can power a jet through neutrino-antineutrino annihilation ($\nu\bar{\nu} \to e^+e^-$) along the rotation axis, or through the Blandford-Znajek mechanism extracting the black hole's spin energy. The jet power from BZ extraction is:

$$P_{\text{jet}} \sim 10^{50}\left(\frac{a}{0.9}\right)^2\left(\frac{M}{3\,M_\odot}\right)^2\left(\frac{B}{10^{15}\,\text{G}}\right)^2\;\text{erg/s}$$

Millisecond Magnetar

Alternatively, a rapidly rotating ($P \sim 1$ ms) magnetar with surface field $B \sim 10^{15}\text{--}10^{16}$ G can power a GRB through its spin-down luminosity. The rotational energy is:

$$E_{\text{rot}} = \frac{1}{2}I\Omega^2 \approx 2 \times 10^{52}\left(\frac{P}{1\,\text{ms}}\right)^{-2}\;\text{erg}$$

This energy is released over the spin-down timescale, producing a plateau phase in the X-ray afterglow that has been observed in many GRBs by the Swift satellite. The magnetar model can also explain extended emission episodes and the energy injection signatures seen in afterglow light curves.

The Prompt Emission Spectrum

The prompt gamma-ray emission spectrum is empirically described by the Band function (Band et al., 1993), a smoothly-broken power law:

$$N(E) \propto \begin{cases} E^\alpha\,\exp(-E/E_0) & E < (\alpha - \beta)E_0 \\ E^\beta & E > (\alpha - \beta)E_0 \end{cases}$$

with typical parameters $\alpha \approx -1$ (low-energy index),$\beta \approx -2.3$ (high-energy index), and peak energy $E_{\text{peak}} \approx 200\text{--}300$ keV. The physical origin of the Band spectrum remains debated: synchrotron radiation in a magnetized outflow, photospheric emission from the expanding fireball, or inverse Compton scattering all produce spectra that can approximate the Band function under appropriate conditions.

The spectral-energy correlations (Amati relation, Ghirlanda relation) connect the intrinsic properties of the prompt emission:

$$E_{\text{peak,rest}} \propto E_{\text{iso}}^{0.5}, \qquad E_{\text{peak,rest}} \propto E_\gamma^{0.7}$$

where $E_\gamma$ is the collimation-corrected energy. These correlations have been proposed as distance indicators for GRBs, which would extend the Hubble diagram to $z > 8$, though their use remains controversial due to selection effects and the lack of a clear physical derivation.

Polarization as a Diagnostic

The polarization of GRB prompt emission constrains the magnetic field geometry in the emitting region. Synchrotron radiation from an ordered magnetic field produces $\Pi \sim 60\text{--}70\%$ linear polarization, while a random field yields $\Pi \sim 0$. Measurements by POLAR and GAP find typical polarization of $\sim 10\text{--}30\%$, suggesting partially ordered fields. Afterglow polarization measurements show evolution from $\sim 1\text{--}3\%$(consistent with a shock-compressed random field) with a smooth rotation of the polarization angle across the jet break, confirming the collimated jet geometry.

Computational Exploration

The following simulation models the GRB fireball dynamics, computes synchrotron afterglow light curves, illustrates the compactness problem resolution, and generates the duration-hardness distribution of GRBs.

GRB Fireball Dynamics, Afterglow Light Curves, and Duration Distribution

Python

script.py244 lines

import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt

fig, axes = plt.subplots(2, 2, figsize=(14, 11))
fig.patch.set_facecolor('#0a0a0a')
for ax in axes.flat:
    ax.set_facecolor('#0f0f0f')
    ax.tick_params(colors='white')
    ax.xaxis.label.set_color('white')
    ax.yaxis.label.set_color('white')
    ax.title.set_color('white')
    for spine in ax.spines.values():
        spine.set_color('#333')

rng = np.random.default_rng(42)

# --- Plot 1: Fireball dynamics - Lorentz factor vs radius ---
ax1 = axes[0, 0]

# Fireball parameters
E_iso = 1e52  # erg
eta_values = [100, 300, 1000]  # baryon loading parameter
colors_eta = ['#44bbff', '#ffa500', '#ff4444']
R_0 = 1e7  # initial radius (cm)

r_range = np.logspace(7, 18, 1000)

for eta, col in zip(eta_values, colors_eta):
    Gamma = np.zeros_like(r_range)

# Acceleration phase: Gamma ~ r/R_0 until Gamma = eta
    R_sat = R_0 * eta  # saturation radius

for i, r in enumerate(r_range):
        if r < R_sat:
            Gamma[i] = r / R_0
        else:
            # Coasting phase until deceleration
            # Deceleration radius
            n_ext = 1.0  # cm^-3
            m_p = 1.67e-24
            R_dec = (3 * E_iso / (4 * np.pi * n_ext * m_p * 3e10**2 * eta**2))**(1.0/3.0)

if r < R_dec:
                Gamma[i] = eta
            else:
                # Blandford-McKee deceleration: Gamma ~ (R/R_dec)^{-3/2} * eta
                Gamma[i] = eta * (r / R_dec)**(-1.5)

Gamma = np.maximum(Gamma, 1.0)
    ax1.loglog(r_range, Gamma, color=col, linewidth=2.5, label=f'eta = {eta}')

# Mark key radii
ax1.axvline(x=1e7, color='#666', linestyle=':', linewidth=1)
ax1.text(1.2e7, 5e2, 'R_0', color='#888', fontsize=9)
ax1.axvline(x=3e13, color='#666', linestyle=':', linewidth=1)
ax1.text(4e13, 5e2, 'R_IS', color='#888', fontsize=9)
ax1.axvline(x=1e17, color='#666', linestyle=':', linewidth=1)
ax1.text(1.2e17, 5e2, 'R_dec', color='#888', fontsize=9)

ax1.set_xlabel('Radius (cm)', fontsize=12)
ax1.set_ylabel('Lorentz Factor', fontsize=12)
ax1.set_title('Fireball Dynamics', fontsize=14, fontweight='bold')
ax1.legend(fontsize=10, facecolor='#1a1a1a', edgecolor='#444', labelcolor='white')
ax1.set_xlim(1e7, 1e19)
ax1.set_ylim(1, 3e3)
ax1.grid(True, alpha=0.15)

# --- Plot 2: Afterglow light curves (multi-band) ---
ax2 = axes[0, 1]

# Simplified afterglow model: F_nu(t) power-law segments
# For ISM (constant density), slow cooling (nu_m < nu < nu_c):
# F_nu ~ t^{3(1-p)/4} * nu^{-(p-1)/2}
# For t < t_jet: standard decay
# After jet break: steeper decline

p = 2.3
t_obs = np.logspace(-2, 3, 500)  # days

# Band frequencies (relative to nu_m at 1 day)
bands = [
    ('X-ray (1 keV)', 1e4, '#44bbff'),
    ('Optical (R-band)', 1.0, '#ffa500'),
    ('Radio (8.5 GHz)', 1e-4, '#ff4444'),
]

t_jet = 3.0  # jet break time in days

for label, nu_ratio, col in bands:
    # Pre-jet-break
    if nu_ratio > 1:  # above nu_c
        alpha = (2 - 3*p) / 4  # temporal index
        beta = -p / 2  # spectral index
    elif nu_ratio > 0.01:  # between nu_m and nu_c
        alpha = 3*(1-p)/4
        beta = -(p-1)/2
    else:  # below nu_m (radio)
        alpha = 0.5
        beta = 1.0/3.0

F_norm = 1e-3 * nu_ratio**beta  # mJy at 1 day

F_ag = np.zeros_like(t_obs)
    for i, t in enumerate(t_obs):
        if t < 0.01:
            F_ag[i] = F_norm * 0.01  # rising
        elif t < t_jet:
            F_ag[i] = F_norm * t**alpha
        else:
            # Post jet break: steeper by ~t^{-1}
            F_ag[i] = F_norm * t_jet**alpha * (t/t_jet)**(alpha - 0.75)

# Radio: add scintillation and self-absorption effects
    if nu_ratio < 1e-3:
        F_ag[:50] *= np.linspace(0.01, 1, 50)  # delayed rise
        F_ag *= (1 + 0.3 * rng.normal(0, 1, len(F_ag)))
        F_ag = np.maximum(F_ag, 0)

ax2.loglog(t_obs, np.maximum(F_ag, 1e-8), color=col, linewidth=2.5, label=label)

# Mark jet break
ax2.axvline(x=t_jet, color='#888', linestyle='--', linewidth=1.5, alpha=0.5)
ax2.text(t_jet*1.3, 1e-1, 'Jet break', color='#aaa', fontsize=10)

ax2.set_xlabel('Time since GRB (days)', fontsize=12)
ax2.set_ylabel('Flux density (mJy)', fontsize=12)
ax2.set_title('GRB Afterglow Light Curves', fontsize=14, fontweight='bold')
ax2.legend(fontsize=10, facecolor='#1a1a1a', edgecolor='#444', labelcolor='white')
ax2.set_xlim(0.01, 1000)
ax2.set_ylim(1e-7, 10)
ax2.grid(True, alpha=0.15)

# --- Plot 3: Duration-hardness diagram ---
ax3 = axes[1, 0]

# Generate synthetic GRB population
# Long GRBs
N_long = 400
T90_long = 10**(rng.normal(1.5, 0.4, N_long))  # centered at ~30 s
HR_long = 10**(rng.normal(-0.1, 0.3, N_long))  # hardness ratio

# Short GRBs
N_short = 100
T90_short = 10**(rng.normal(-0.5, 0.3, N_short))  # centered at ~0.3 s
HR_short = 10**(rng.normal(0.4, 0.25, N_short))  # harder spectrum

ax3.scatter(T90_long, HR_long, s=8, color='#ffa500', alpha=0.5, label='Long GRBs')
ax3.scatter(T90_short, HR_short, s=15, color='#44bbff', alpha=0.7, label='Short GRBs')

# Division line
ax3.axvline(x=2.0, color='#ff4444', linestyle='--', linewidth=2, alpha=0.7, label='T90 = 2 s')

ax3.set_xscale('log')
ax3.set_yscale('log')
ax3.set_xlabel('T90 Duration (s)', fontsize=12)
ax3.set_ylabel('Hardness Ratio (H/S)', fontsize=12)
ax3.set_title('GRB Duration-Hardness Distribution', fontsize=14, fontweight='bold')
ax3.legend(fontsize=10, facecolor='#1a1a1a', edgecolor='#444', labelcolor='white')
ax3.set_xlim(0.01, 1000)
ax3.set_ylim(0.1, 10)
ax3.grid(True, alpha=0.15)

# Annotate progenitors
ax3.text(0.15, 5, 'NS-NS / NS-BH\nmergers', color='#44bbff', fontsize=10, ha='center')
ax3.text(50, 0.2, 'Collapsars\n(massive star death)', color='#ffa500', fontsize=10, ha='center')

# --- Plot 4: Compactness resolution ---
ax4 = axes[1, 1]

Gamma_range = np.logspace(0, 3.5, 200)
p_index = 2.0  # photon index

# Optical depth to pair production as function of Gamma
L_gamma = 1e52  # erg/s
delta_t = 0.01  # seconds
sigma_T = 6.65e-25
m_e = 9.109e-28
c = 3e10

tau_static = sigma_T * L_gamma / (4 * np.pi * (c * delta_t)**2 * m_e * c**2 / delta_t)
tau_GG = tau_static * Gamma_range**(-(4 + 2*p_index))

ax4.loglog(Gamma_range, tau_GG, color='#ffa500', linewidth=3, label=r'$	au_{gammagamma}(Gamma)$')
ax4.axhline(y=1.0, color='#ff4444', linestyle='--', linewidth=2, alpha=0.7, label=r'$	au = 1$ (transparency)')

# Mark minimum Gamma
idx_trans = np.argmin(np.abs(tau_GG - 1.0))
Gamma_min = Gamma_range[idx_trans]
ax4.axvline(x=Gamma_min, color='#44bbff', linestyle=':', linewidth=1.5, alpha=0.7)
ax4.text(Gamma_min * 1.3, 1e5, f'Gamma_min ~ {Gamma_min:.0f}', color='#44bbff', fontsize=11)

# Shaded regions
ax4.fill_between(Gamma_range, 1, 1e20, where=tau_GG > 1, alpha=0.1, color='#ff4444')
ax4.fill_between(Gamma_range, 1e-15, 1, where=tau_GG < 1, alpha=0.1, color='#44bbff')

ax4.text(3, 1e8, 'OPAQUE\n(no gamma-rays\nescape)', color='#ff6644', fontsize=11, ha='center')
ax4.text(500, 1e-5, 'TRANSPARENT\n(gamma-rays\nescape)', color='#44bbff', fontsize=11, ha='center')

ax4.set_xlabel('Bulk Lorentz Factor (Gamma)', fontsize=12)
ax4.set_ylabel('Pair-production optical depth', fontsize=12)
ax4.set_title('Compactness Problem Resolution', fontsize=14, fontweight='bold')
ax4.legend(fontsize=10, facecolor='#1a1a1a', edgecolor='#444', labelcolor='white')
ax4.set_xlim(1, 3000)
ax4.set_ylim(1e-10, 1e15)
ax4.grid(True, alpha=0.15)

# Print summary
print("=" * 65)
print("GAMMA-RAY BURSTS: KEY RESULTS")
print("=" * 65)

print("\nCompactness problem:")
print(f"  L_gamma = {L_gamma:.0e} erg/s, delta_t = {delta_t*1e3:.0f} ms")
print(f"  Naive source size: R ~ c*dt = {c*delta_t:.2e} cm")
print(f"  Pair optical depth (static): tau ~ {tau_static:.2e}")
print(f"  Minimum Lorentz factor: Gamma_min ~ {Gamma_min:.0f}")

print("\nFireball parameters:")
for eta in eta_values:
    R_sat = R_0 * eta
    n_ext = 1.0
    R_dec = (3*E_iso/(4*np.pi*n_ext*m_p*c**2*eta**2))**(1.0/3.0)
    print(f"  eta = {eta:4d}: R_sat = {R_sat:.2e} cm, R_dec = {R_dec:.2e} cm")

print("\nBeaming corrections:")
for theta_j_deg in [2, 5, 10, 20]:
    theta_j = theta_j_deg * np.pi / 180
    f_b = theta_j**2 / 2
    E_true = f_b * 1e53
    print(f"  theta_j = {theta_j_deg:2d} deg: f_b = {f_b:.4f}, E_true = {E_true:.2e} erg")

print("\nAfterglow temporal indices (p = 2.3):")
print(f"  nu_m < nu < nu_c: alpha = {3*(1-p)/4:.2f}")
print(f"  nu > nu_c:        alpha = {(2-3*p)/4:.2f}")
print(f"  Post jet break:   alpha ~ -2")

m_p = 1.67e-24
plt.tight_layout()
plt.savefig('output.png', dpi=150, bbox_inches='tight', facecolor='#0a0a0a')
print("\nPlot saved to output.png")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Chapter Summary

Gamma-ray bursts require ultra-relativistic bulk motion ($\Gamma \gtrsim 100$) to resolve the compactness problem. The fireball model describes how energy deposited near a compact object drives a baryon-loaded outflow that converts kinetic energy to radiation through internal shocks (prompt emission) and external shocks (afterglow).

Relativistic beaming collimates the emission into jets with typical opening angles of $5°\text{--}10°$, reducing the true energy to $\sim 10^{51}$ erg. The achromatic jet break in afterglow light curves constrains $\theta_j$.

Long GRBs originate from collapsars (the death of massive stars), while short GRBs arise from compact binary mergers — confirmed definitively by the multi-messenger detection of GW170817/GRB 170817A. GRBs serve as unique probes of the high-redshift Universe and the physics of ultra-relativistic outflows.

← Supernovae Part III: Galaxies & Cosmology →