Part IV: Frontiers | Chapter 3

Multi-Messenger Astronomy

Combining gravitational waves, electromagnetic radiation, neutrinos, and cosmic rays to reveal the most extreme phenomena in the Universe

Overview

Multi-messenger astronomy observes the Universe through four distinct channels: electromagnetic radiation (photons from radio to gamma-ray), gravitational waves (ripples in spacetime), neutrinos (nearly massless particles that escape the densest environments), and cosmic rays (ultra-high-energy charged particles). The landmark event GW170817 — the first joint detection of gravitational waves and electromagnetic emission from a binary neutron star merger — inaugurated this new era on 17 August 2017, yielding discoveries that no single messenger could provide alone.

In this chapter we derive the gravitational wave signal from compact binary mergers, analyze the electromagnetic counterparts of GW170817, discuss neutrino astronomy from supernovae and active galaxies, and examine ultra-high-energy cosmic ray physics.

1. Gravitational Waves from Compact Binaries

The inspiral of two compact objects (neutron stars or black holes) produces gravitational waves whose frequency and amplitude increase as the orbit shrinks — the characteristic "chirp" signal.

1.1 The Chirp Signal

For a circular binary with component masses $m_1$ and $m_2$at distance $d$, the gravitational wave strain amplitude is:

$$\boxed{h = \frac{4}{d}\left(\frac{G\mathcal{M}}{c^2}\right)^{5/3}\left(\frac{\pi f_{\text{GW}}}{c}\right)^{2/3}}$$

where the chirp mass is:

$$\mathcal{M} = \frac{(m_1 m_2)^{3/5}}{(m_1 + m_2)^{1/5}}$$

The gravitational wave frequency is twice the orbital frequency: $f_{\text{GW}} = 2f_{\text{orb}}$. As the binary loses energy to gravitational radiation, the frequency evolves as:

$$\boxed{\dot{f}_{\text{GW}} = \frac{96}{5}\pi^{8/3}\left(\frac{G\mathcal{M}}{c^3}\right)^{5/3}f_{\text{GW}}^{11/3}}$$

The chirp mass is the best-measured parameter from the GW signal, determined to sub-percent precision for strong detections. For GW170817: $\mathcal{M} = 1.188 \pm 0.004\,M_\odot$.

1.2 Energy and Luminosity

The gravitational wave luminosity (quadrupole formula) is:

$$L_{\text{GW}} = \frac{32}{5}\frac{G^4}{c^5}\frac{(m_1 m_2)^2(m_1+m_2)}{a^5}$$

At merger, the peak luminosity reaches $\sim 3.6 \times 10^{56}$ erg/s for a binary neutron star and $\sim 10^{57}$ erg/s for stellar-mass BBH — briefly exceeding the combined electromagnetic luminosity of all stars in the observable Universe.

2. GW170817: The Rosetta Stone Event

On 17 August 2017, LIGO and Virgo detected a 100-second gravitational wave signal from the inspiral and merger of two neutron stars at a distance of $\sim 40$ Mpc. The electromagnetic follow-up campaign across the spectrum revealed an extraordinary sequence of phenomena.

2.1 The Short Gamma-Ray Burst

Fermi-GBM detected a short GRB (GRB 170817A) arriving 1.7 seconds after the merger, confirming the long-hypothesized connection between short GRBs and neutron star mergers. The delay constrains the speed of gravity:

$$\boxed{-3 \times 10^{-15} \leq \frac{v_{\text{GW}} - c}{c} \leq 7 \times 10^{-16}}$$

This extraordinary constraint rules out many modified gravity theories and confirms that gravitational waves travel at the speed of light to parts per quadrillion.

2.2 The Kilonova

The optical/infrared transient AT2017gfo was identified as a kilonova — thermal emission from radioactive r-process ejecta. The light curve showed two components:

Blue component (first 2 days): lanthanide-poor ejecta with $\kappa \sim 1$ cm$^2$/g, producing emission peaking at $\sim 10^{42}$ erg/s.

Red component (days 3–14): lanthanide-rich ejecta with $\kappa \sim 10$ cm$^2$/g (due to the complex atomic structure of lanthanides), producing redder, longer-lasting emission.

The total ejecta mass was $\sim 0.05\,M_\odot$, confirming that neutron star mergers produce substantial quantities of r-process elements.

2.3 The Hubble Constant Measurement

Gravitational waves provide a direct measurement of the luminosity distance (they are "standard sirens"), while the host galaxy NGC 4993 provides the redshift. Combined:

$$H_0 = 70.0^{+12}_{-8}\;\text{km s}^{-1}\;\text{Mpc}^{-1}$$

With more events, standard sirens will provide an independent measurement of $H_0$, potentially resolving the Hubble tension between CMB and local distance ladder measurements.

3. Neutrino Astronomy

Neutrinos interact so weakly that they escape from the densest astrophysical environments, carrying information about the interior physics of their sources.

3.1 Supernova Neutrinos

Core-collapse supernovae release $\sim 3 \times 10^{53}$ erg in neutrinos (99% of the gravitational binding energy). The neutrino luminosity is $L_\nu \sim 3 \times 10^{53}/10$ erg/s $\sim 10^{52}$ erg/s over $\sim 10$ seconds. The expected event rate in a detector with $N_t$target protons at distance $d$ is:

$$N_{\text{events}} \approx \frac{E_\nu^{\text{tot}}}{4\pi d^2 \langle E_\nu\rangle}\,\sigma_{\bar{\nu}_e p}\,N_t$$

For SN 1987A at $d = 50$ kpc, Kamiokande-II detected 11 events, IMB detected 8, and Baksan detected 5 — confirming the theoretical picture. A Galactic supernova would produce $\sim 10{,}000$ events in Super-Kamiokande and $\sim 100{,}000$in Hyper-Kamiokande.

3.2 High-Energy Neutrinos from AGN

IceCube, a cubic-kilometer neutrino detector at the South Pole, has detected astrophysical neutrinos with energies up to $\sim 10$ PeV. In September 2017 (IceCube-170922A), a $\sim 290$ TeV neutrino was traced to the blazar TXS 0506+056, providing the first identification of an extragalactic neutrino source. The neutrino production mechanism involves proton-photon interactions in the relativistic jet:

$$p + \gamma \to \Delta^+ \to \begin{cases} n + \pi^+ \to n + \mu^+ + \nu_\mu \\ p + \pi^0 \to p + 2\gamma \end{cases}$$

4. Ultra-High-Energy Cosmic Rays

Cosmic rays are charged particles (mostly protons and heavier nuclei) accelerated to extraordinary energies, up to $\sim 10^{20}$ eV — the $10^8$ times the energy of the LHC.

4.1 The Hillas Criterion

A particle of charge $Ze$ can be confined and accelerated in a region of size $R$ with magnetic field $B$ up to a maximum energy:

$$\boxed{E_{\text{max}} = Ze\beta_s cBR \approx 10^{18}\,Z\,\beta_s\left(\frac{B}{1\,\mu\text{G}}\right)\left(\frac{R}{1\,\text{kpc}}\right)\;\text{eV}}$$

where $\beta_s$ is the shock velocity in units of $c$. Only a few source classes can accelerate particles to $> 10^{19}$ eV: AGN jets, GRB shocks, and galaxy cluster accretion shocks.

4.2 The GZK Cutoff

Protons with energies above $\sim 5 \times 10^{19}$ eV interact with CMB photons via pion production:

$$p + \gamma_{\text{CMB}} \to \Delta^+ \to p + \pi^0 \;\;\text{or}\;\; n + \pi^+$$

The threshold energy in the proton rest frame requires the CMB photon to appear as a$\sim 300$ MeV photon, achieved when $E_p \gtrsim 5 \times 10^{19}$ eV. This limits the propagation distance of the most energetic cosmic rays to $\sim 100$ Mpc (the GZK horizon), confirmed by the Pierre Auger Observatory.

5. Gravitational Wave Detector Physics

LIGO and Virgo are laser interferometers that measure the differential arm length change induced by a passing gravitational wave.

5.1 Strain Sensitivity

A gravitational wave of strain $h$ changes the differential arm length by:

$$\Delta L = \frac{1}{2}h L$$

For LIGO ($L = 4$ km) and a typical binary merger signal ($h \sim 10^{-21}$):$\Delta L \sim 2 \times 10^{-18}$ m — about one-thousandth the diameter of a proton. Achieving this extraordinary sensitivity requires overcoming quantum shot noise, thermal noise, seismic noise, and Newtonian gravity gradient noise.

5.2 The Detection Horizon

The detection range for a given source type scales as:

$$d_{\text{horizon}} \propto \mathcal{M}^{5/6} \cdot h_{\text{rms}}^{-1}$$

LIGO at design sensitivity (O4) can detect binary neutron star mergers to $\sim 190$ Mpc and binary black hole mergers to $\sim\text{several Gpc}$. Next-generation detectors (Einstein Telescope, Cosmic Explorer) will extend the BNS range to $z \sim 2$.

Applications

Testing General Relativity

Gravitational wave observations test GR in the strong-field, highly dynamical regime inaccessible to other experiments. The consistency of the observed waveforms with GR predictions (no-hair theorem, quasi-normal mode spectrum of the merger remnant) has been confirmed to remarkable precision. The speed of gravity measurement from GW170817 eliminated a wide class of modified gravity theories.

The Nanohertz Gravitational Wave Background

Pulsar timing arrays (NANOGrav, EPTA, PPTA) reported evidence in 2023 for a stochastic gravitational wave background at nanohertz frequencies, likely produced by the cosmic population of supermassive black hole binaries. This represents the lowest-frequency gravitational wave detection and opens a new window on SMBH assembly and galaxy evolution.

Historical Notes

Multi-messenger astronomy has roots in the coincident detection of solar neutrinos and sunlight (by definition) and the detection of neutrinos from SN 1987A. The LIGO project was conceived by Rainer Weiss, Kip Thorne, and Barry Barish in the 1980s–1990s. The first gravitational wave detection (GW150914, two merging black holes) on 14 September 2015 earned Weiss, Thorne, and Barish the 2017 Nobel Prize in Physics. GW170817, detected just two years later, produced observations across the entire electromagnetic spectrum by over 70 observatories worldwide, demonstrating the power of multi-messenger astronomy. IceCube's identification of TXS 0506+056 as a neutrino source in 2017 added the neutrino channel to the multi-messenger toolkit.

The theoretical prediction of gravitational waves dates to Einstein's 1916 paper, though Einstein himself doubted their physical reality. The question was settled by the observation of orbital decay in the Hulse-Taylor binary pulsar (PSR B1913+16, discovered 1974), which matches the GR prediction for gravitational wave emission to 0.2% precision. Joseph Weber's resonant bar detector claims in the 1960s, though never confirmed, stimulated the development of laser interferometric detectors. The 40-year journey from concept (Weiss 1972) to detection (2015) required solving extraordinary technical challenges in seismic isolation, mirror coating, laser stability, and quantum noise reduction. The current network of Advanced LIGO (Hanford and Livingston), Advanced Virgo, and KAGRA detects binary mergers at a rate of approximately one per week, building a catalog that constrains the mass spectrum of black holes, the neutron star equation of state, and the rate of cosmic compact binary mergers.

The observation of GW190814 revealed a $23\,M_\odot$ black hole merging with a $2.6\,M_\odot$ compact object — either the lightest black hole or the most massive neutron star ever observed. The nature of the secondary component depends on the maximum neutron star mass, which in turn depends on the nuclear equation of state. The lack of an electromagnetic counterpart (expected if the secondary were a neutron star disrupted by tidal forces) favors a black hole interpretation, but the question remains open. The accumulating gravitational wave catalog is mapping the mass distribution of compact objects with increasing precision, revealing features like the lower mass gap and the pair-instability mass gap that encode the physics of massive stellar evolution and death.

Future Multi-Messenger Science

The next decade promises dramatic advances in multi-messenger astronomy. Third-generation gravitational wave detectors (Einstein Telescope, Cosmic Explorer) will detect binary neutron star mergers to $z \sim 2$, enabling hundreds of standard siren measurements of $H_0$. LISA (launch ~2035) will detect supermassive black hole mergers, extreme mass-ratio inspirals, and thousands of Galactic binaries. IceCube-Gen2 will increase the neutrino detection volume tenfold. The Cherenkov Telescope Array (CTA) will survey the TeV gamma-ray sky with unprecedented sensitivity. Combined with wide-field optical surveys (LSST), these facilities will enable routine multi- messenger follow-up of the most violent events in the Universe.

Diffuse Neutrino and Gravitational Wave Backgrounds

Beyond individual sources, the cumulative emission from all sources in the Universe produces diffuse backgrounds. IceCube has measured the diffuse astrophysical neutrino flux with a spectrum approximately $E^{-2.5}$, though the dominant sources remain uncertain (AGN, starburst galaxies, and galaxy clusters all contribute). The stochastic gravitational wave background at nanohertz frequencies, detected by pulsar timing arrays in 2023, likely originates from the cosmic population of inspiraling supermassive black hole binaries. At higher frequencies, the unresolved superposition of compact binary mergers creates a background detectable by LIGO at design sensitivity. These diffuse backgrounds carry information about the entire population of sources across cosmic history, complementing the study of individual events.

The Galactic Neutrino Sky

In addition to extragalactic point sources, IceCube has detected diffuse neutrino emission from the Galactic plane, produced by cosmic ray interactions with interstellar gas ($p + p \to \pi^\pm + X$, followed by $\pi^\pm \to \mu^\pm + \nu_\mu$). This detection, announced in 2023, confirms that the Milky Way is a source of high-energy neutrinos and provides the first neutrino "image" of our Galaxy. The neutrino emission traces cosmic ray propagation and the distribution of target gas, complementing gamma-ray observations that suffer from the ambiguity between hadronic ($\pi^0 \to 2\gamma$) and leptonic (inverse Compton) emission mechanisms.

The Gravitational Wave Transient Catalog

Through the first four observing runs (O1–O4), LIGO/Virgo/KAGRA have detected over 200 compact binary mergers, revealing the mass spectrum of stellar-mass black holes and neutron stars.

Key Population Features

The black hole mass spectrum shows several notable features: a peak at $\sim 8\text{--}10\,M_\odot$, a possible gap in the range $3\text{--}5\,M_\odot$ between neutron stars and black holes (the "lower mass gap"), and events in the pair-instability mass gap ($50\text{--}130\,M_\odot$) such as GW190521 (a $85 + 66\,M_\odot$merger). The existence of BHs in the pair-instability gap suggests hierarchical mergers in dense stellar environments (globular clusters, nuclear star clusters).

The spin distribution of merging black holes provides clues to their formation channels. Dynamically assembled binaries (in dense clusters) have isotropic spin orientations, while isolated binary evolution preferentially produces aligned spins. The observed distribution shows a mixture of both channels, with effective spin parameters concentrated near zero but with a tail to positive values.

The merger rate density is measured as $\mathcal{R}_{\text{BBH}} \approx 17\text{--}45$ Gpc$^{-3}$ yr$^{-1}$for binary black holes and $\mathcal{R}_{\text{BNS}} \approx 10\text{--}1700$Gpc$^{-3}$ yr$^{-1}$ for binary neutron stars. These rates are consistent with population synthesis models of massive binary stellar evolution.

Extreme Mass-Ratio Inspirals (EMRIs)

A stellar-mass compact object ($\sim 10\,M_\odot$) spiraling into a supermassive black hole ($\sim 10^6\,M_\odot$) produces a long-lived gravitational wave signal in the millihertz band, detectable by LISA. EMRIs trace$\sim 10^5$ orbits before plunge, mapping the spacetime geometry of the SMBH with exquisite precision and testing the no-hair theorem of general relativity (verifying that the SMBH is described by the Kerr metric). LISA is expected to detect tens to hundreds of EMRIs per year, providing precision measurements of SMBH masses and spins across cosmic time and probing the stellar dynamics in galactic nuclei.

Computational Exploration

The following simulation generates a binary inspiral gravitational wave chirp signal, models the kilonova light curve from GW170817, illustrates the cosmic ray energy spectrum, and computes the LIGO noise curve with detection horizons.

GW Chirp Signal, Kilonova Light Curve, and Cosmic Ray Spectrum

Python

script.py259 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')

# Constants
G = 6.674e-11
c = 3e8
M_sun = 1.989e30
Mpc = 3.086e22

# --- Plot 1: GW chirp signal ---
ax1 = axes[0, 0]

m1 = 1.4 * M_sun
m2 = 1.3 * M_sun
M_total = m1 + m2
M_chirp = (m1 * m2)**0.6 / M_total**0.2
d = 40 * Mpc  # GW170817 distance

# Time to merger from frequency f
# t(f) = 5/(256 pi^{8/3}) * (G*Mc/c^3)^{-5/3} * f^{-8/3}
f_isco = c**3 / (6**1.5 * np.pi * G * M_total)  # ISCO frequency

# Generate chirp
t_merge = 0  # merger at t=0
t_start = -30  # seconds before merger
dt = 1e-4
t = np.arange(t_start, -0.01, dt)

# Time to coalescence from current time
tau = -t  # time remaining to merger
tau = np.maximum(tau, 1e-10)

# Frequency evolution
f_gw = (5/(256 * tau))**(3.0/8.0) / np.pi * (G * M_chirp / c**3)**(-5.0/8.0)
f_gw = np.minimum(f_gw, f_isco)

# Strain amplitude
h_plus = 4 / d * (G * M_chirp / c**2)**(5.0/3.0) * (np.pi * f_gw / c)**(2.0/3.0)

# Phase
phi = np.cumsum(2 * np.pi * f_gw * dt)

# Waveform
h_t = h_plus * np.cos(phi)

# Plot last 2 seconds
mask = t > -2
ax1.plot(t[mask], h_t[mask] * 1e21, color='#ffa500', linewidth=1)

# Envelope
ax1.plot(t[mask], h_plus[mask] * 1e21, '--', color='#ff4444', linewidth=1, alpha=0.5)
ax1.plot(t[mask], -h_plus[mask] * 1e21, '--', color='#ff4444', linewidth=1, alpha=0.5)

ax1.set_xlabel('Time to merger (s)', fontsize=12)
ax1.set_ylabel('Strain h (x 10^-21)', fontsize=12)
ax1.set_title('GW170817-like Chirp Signal', fontsize=14, fontweight='bold')
ax1.set_xlim(-2, 0)
ax1.grid(True, alpha=0.15)

# Annotate
ax1.text(-1.8, max(h_plus[mask]*1e21)*0.8, f'Mc = {M_chirp/M_sun:.3f} M_sun\nd = 40 Mpc',
         color='#ffdd00', fontsize=10)

# --- Plot 2: Kilonova light curve ---
ax2 = axes[0, 1]

t_days = np.logspace(-1, 2, 200)

# Blue kilonova component
M_ej_blue = 0.015 * M_sun * 1e3  # grams
v_ej_blue = 0.3 * 3e10  # cm/s
kappa_blue = 1.0  # cm^2/g
t_m_blue = np.sqrt(2 * kappa_blue * M_ej_blue / (v_ej_blue * 3e10))  # seconds
t_m_blue_d = t_m_blue / 86400  # days

# Simple model: L ~ L_peak * exp(-(t/t_m)^2 - t_m/t)
L_peak_blue = 2e41  # erg/s
L_blue = L_peak_blue * np.exp(-(t_days / t_m_blue_d)**2) * np.exp(-0.5 * t_m_blue_d / t_days)
L_blue *= (t_days / t_m_blue_d)**(-1.3)
L_blue = np.minimum(L_blue, L_peak_blue)

# Red kilonova component
M_ej_red = 0.04 * M_sun * 1e3
kappa_red = 10.0  # lanthanide-rich
v_ej_red = 0.15 * 3e10
t_m_red = np.sqrt(2 * kappa_red * M_ej_red / (v_ej_red * 3e10))
t_m_red_d = t_m_red / 86400

L_peak_red = 5e40
L_red = L_peak_red * np.exp(-(t_days / t_m_red_d)) * (t_days / 1.0)**(-1.0)
L_red = np.minimum(L_red, L_peak_red)

L_total = L_blue + L_red
L_sun = 3.828e33

ax2.loglog(t_days, L_blue / L_sun, '--', color='#44bbff', linewidth=2, label='Blue (lanthanide-poor)')
ax2.loglog(t_days, L_red / L_sun, '--', color='#ff4444', linewidth=2, label='Red (lanthanide-rich)')
ax2.loglog(t_days, L_total / L_sun, color='#ffa500', linewidth=3, label='Total')

# Simulated AT2017gfo data
rng = np.random.default_rng(42)
t_obs = [0.5, 1.0, 1.5, 2.0, 3.0, 4.0, 5.0, 7.0, 10.0, 14.0]
L_obs = np.interp(t_obs, t_days, L_total / L_sun)
L_obs *= 10**(rng.normal(0, 0.1, len(t_obs)))
ax2.scatter(t_obs, L_obs, color='#ffdd00', s=50, zorder=5, edgecolor='white',
            linewidth=0.5, label='AT2017gfo data')

ax2.set_xlabel('Time since merger (days)', fontsize=12)
ax2.set_ylabel('Luminosity (L_sun)', fontsize=12)
ax2.set_title('Kilonova Light Curve (GW170817)', fontsize=14, fontweight='bold')
ax2.legend(fontsize=9, facecolor='#1a1a1a', edgecolor='#444', labelcolor='white')
ax2.set_xlim(0.1, 100)
ax2.set_ylim(1e4, 1e9)
ax2.grid(True, alpha=0.15)

# --- Plot 3: Cosmic ray energy spectrum ---
ax3 = axes[1, 0]

E_range = np.logspace(9, 21, 300)  # eV

# Approximate all-particle spectrum
# Below knee (3e15 eV): E^-2.7
# Between knee and ankle (3e18): E^-3.1
# Above ankle: E^-2.6
# GZK suppression above 5e19

flux = np.zeros_like(E_range)
for i, E in enumerate(E_range):
    if E < 3e15:
        flux[i] = 1e24 * E**(-2.7)
    elif E < 3e18:
        flux[i] = 1e24 * (3e15)**(-2.7) * (E/3e15)**(-3.1)
    elif E < 5e19:
        flux[i] = 1e24 * (3e15)**(-2.7) * (3e18/3e15)**(-3.1) * (E/3e18)**(-2.6)
    else:
        flux[i] = 1e24 * (3e15)**(-2.7) * (3e18/3e15)**(-3.1) * (5e19/3e18)**(-2.6) *                   (E/5e19)**(-4.0)

# Plot E^2.6 * flux for better visualization
scaled_flux = E_range**2.6 * flux

ax3.loglog(E_range, scaled_flux, color='#ffa500', linewidth=3)

# Mark features
ax3.axvline(x=3e15, color='#44bbff', linestyle='--', linewidth=1.5, alpha=0.7)
ax3.text(4e15, max(scaled_flux)*0.3, 'Knee', color='#44bbff', fontsize=11, rotation=90)
ax3.axvline(x=3e18, color='#ffdd00', linestyle='--', linewidth=1.5, alpha=0.7)
ax3.text(4e18, max(scaled_flux)*0.3, 'Ankle', color='#ffdd00', fontsize=11, rotation=90)
ax3.axvline(x=5e19, color='#ff4444', linestyle='--', linewidth=1.5, alpha=0.7)
ax3.text(6e19, max(scaled_flux)*0.3, 'GZK', color='#ff4444', fontsize=11, rotation=90)

# Source labels
ax3.text(1e11, max(scaled_flux)*2, 'SNR acceleration', color='#aaa', fontsize=9)
ax3.text(1e18, max(scaled_flux)*2, 'Extragalactic', color='#aaa', fontsize=9)

ax3.set_xlabel('Energy (eV)', fontsize=12)
ax3.set_ylabel('E^2.6 x Flux (arb. units)', fontsize=12)
ax3.set_title('Cosmic Ray Energy Spectrum', fontsize=14, fontweight='bold')
ax3.set_xlim(1e9, 1e21)
ax3.grid(True, alpha=0.15)

# --- Plot 4: LIGO noise curve and detection horizons ---
ax4 = axes[1, 1]

# Simplified LIGO noise curve (amplitude spectral density)
f_det = np.logspace(0.5, 4, 500)  # Hz

# aLIGO design sensitivity (approximate)
# Seismic wall below ~10 Hz, thermal noise ~50-200 Hz, shot noise above
h_noise = np.zeros_like(f_det)
for i, f in enumerate(f_det):
    seismic = 1e-15 * (10/f)**10 if f > 5 else 1e-5
    thermal = 5e-24 * (100/f)**0.5
    shot = 3e-24 * (f/1000)**0.5
    h_noise[i] = np.sqrt(seismic**2 + thermal**2 + shot**2)

# Characteristic strain of sources
# BNS at various distances
f_bns = np.logspace(1, np.log10(f_isco/2), 100)
for d_Mpc, col, lab in [(40, '#ff4444', 'BNS @ 40 Mpc (GW170817)'),
                         (200, '#ffa500', 'BNS @ 200 Mpc'),
                         (1000, '#44bbff', 'BNS @ 1 Gpc')]:
    d_m = d_Mpc * Mpc
    h_char = 4/d_m * (G*M_chirp/c**2)**(5./3) * (np.pi*f_bns/c)**(2./3) * np.sqrt(f_bns)
    # Characteristic strain: h_c = h * sqrt(f * N_cycles)
    N_cyc = f_bns**2 / (96/5 * np.pi**(8./3) * (G*M_chirp/c**3)**(5./3) * f_bns**(11./3))
    h_c = h_char * np.sqrt(np.maximum(N_cyc, 1))
    ax4.loglog(f_bns, h_c, color=col, linewidth=2, label=lab)

# BBH (30+30 Msun)
M_chirp_bbh = (30*30)**0.6 / 60**0.2 * M_sun
f_isco_bbh = c**3 / (6**1.5 * np.pi * G * 60 * M_sun)
f_bbh = np.logspace(0.5, np.log10(f_isco_bbh), 80)
d_bbh = 1000 * Mpc
h_bbh = 4/d_bbh * (G*M_chirp_bbh/c**2)**(5./3) * (np.pi*f_bbh/c)**(2./3)
N_cyc_bbh = f_bbh**2 / (96/5 * np.pi**(8./3) * (G*M_chirp_bbh/c**3)**(5./3) * f_bbh**(11./3))
h_c_bbh = h_bbh * np.sqrt(np.maximum(N_cyc_bbh, 1)) * np.sqrt(f_bbh)
ax4.loglog(f_bbh, h_c_bbh, '--', color='#ff69b4', linewidth=2, label='BBH (30+30) @ 1 Gpc')

# Noise curve
ax4.loglog(f_det, h_noise * np.sqrt(f_det), color='white', linewidth=2.5, alpha=0.8, label='aLIGO design')
ax4.fill_between(f_det, h_noise * np.sqrt(f_det), 1e-15, alpha=0.1, color='white')

ax4.set_xlabel('Frequency (Hz)', fontsize=12)
ax4.set_ylabel('Characteristic Strain', fontsize=12)
ax4.set_title('LIGO Sensitivity and Source Signals', fontsize=14, fontweight='bold')
ax4.legend(fontsize=7, facecolor='#1a1a1a', edgecolor='#444', labelcolor='white', loc='upper right')
ax4.set_xlim(3, 5000)
ax4.set_ylim(1e-25, 1e-18)
ax4.grid(True, alpha=0.15)

# Print summary
print("=" * 65)
print("MULTI-MESSENGER ASTRONOMY: KEY RESULTS")
print("=" * 65)

print(f"\nGW170817 parameters:")
print(f"  Chirp mass: {M_chirp/M_sun:.3f} M_sun")
print(f"  Distance: 40 Mpc")
print(f"  GW frequency at ISCO: {f_isco:.0f} Hz")
print(f"  Speed of gravity: |v_GW - c|/c < 7e-16")

print(f"\nKilonova AT2017gfo:")
print(f"  Blue ejecta: ~0.015 M_sun, v ~ 0.3c, kappa ~ 1 cm^2/g")
print(f"  Red ejecta: ~0.04 M_sun, v ~ 0.15c, kappa ~ 10 cm^2/g")
print(f"  Total r-process: ~0.05 M_sun")

print(f"\nH_0 from GW170817: 70 +12/-8 km/s/Mpc")

print(f"\nSN 1987A neutrinos:")
print(f"  Total energy: ~3e53 erg")
print(f"  Events detected: ~24 (across 3 detectors)")
print(f"  Duration: ~10 seconds")

print(f"\nCosmic ray features:")
print(f"  Knee: ~3e15 eV (Galactic SNR limit)")
print(f"  Ankle: ~3e18 eV (Galactic->extragalactic transition)")
print(f"  GZK cutoff: ~5e19 eV (pion production on CMB)")

print(f"\nLIGO strain sensitivity: ~3e-24 /sqrt(Hz) at 100 Hz")
print(f"  BNS range (O4): ~190 Mpc")
print(f"  BBH range: several Gpc")

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

Gravitational waves from compact binary mergers carry the characteristic chirp signal with strain $h \propto \mathcal{M}^{5/3}f^{2/3}/d$. The chirp mass is the best-measured parameter, determined to sub-percent precision.

GW170817 confirmed the connection between neutron star mergers, short GRBs, kilonovae, and r-process nucleosynthesis. It constrained the speed of gravity to parts per quadrillion and provided the first standard siren measurement of $H_0$.

Neutrino astronomy probes the cores of supernovae and the jets of blazars, while cosmic ray observations reveal particle acceleration to $10^{20}$ eV with the GZK cutoff confirming interactions with the CMB. Together, these four messengers provide complementary views of the most extreme phenomena in the Universe.

← Astrobiology Observational Techniques →