Featured image of post The Dawn of Gravitational-Wave Astronomy: Giant Laser Interferometers Capturing Ripples in Spacetime and the Mysteries of Cosmic Genesis

The Dawn of Gravitational-Wave Astronomy: Giant Laser Interferometers Capturing Ripples in Spacetime and the Mysteries of Cosmic Genesis

A miracle 100 years after Einstein's prediction. The incredible measurement precision of LIGO/Virgo/KAGRA and the future of multi-messenger astronomy.

The Dawn of Gravitational-Wave Astronomy: Giant Laser Interferometers Capturing Ripples in Spacetime and the Mysteries of Cosmic Genesis

Humanity’s “eyes” for observing the universe have, throughout recorded history, relied on electromagnetic waves (visible light, radio waves, X-rays, etc.). However, in 2015, we acquired “ears” to listen to the entirely new heartbeat of the cosmos. Those are gravitational waves. In this article, we will thoroughly explore and explain the monumental achievement in the history of physics that is the direct detection of gravitational waves, realized a century after Einstein’s prediction. We will also delve into the pinnacle of human extreme engineering that made it possible, and the future of cosmology being pioneered by multi-messenger astronomy.


Chapter 1: Einstein’s Doubts and the Theory of Gravitational Waves

The concept of gravitational waves naturally derives from the theory of general relativity, completed by Albert Einstein in 1915. In general relativity, gravity is described as the “distortion of spacetime.” When an object with mass accelerates, the distortion of the spacetime around it propagates through space like ripples at the speed of light—this phenomenon is known as “Gravitational Waves.”

Weak-Field Approximation of the Einstein Field Equations and Derivation of the Wave Equation

$$ R_{\mu\nu} - \frac{1}{2}g_{\mu\nu}R = \frac{8\pi G}{c^4} T_{\mu\nu} $$$$ g_{\mu\nu} = \eta_{\mu\nu} + h_{\mu\nu} \quad (|h_{\mu\nu}| \ll 1) $$$$ \bar{h}_{\mu\nu} \equiv h_{\mu\nu} - \frac{1}{2}\eta_{\mu\nu}h $$$$ \Box \bar{h}_{\mu\nu} = -\frac{16\pi G}{c^4} T_{\mu\nu} $$

Here, $\Box = \eta^{\alpha\beta}\partial_\alpha\partial_\beta = -\frac{1}{c^2}\frac{\partial^2}{\partial t^2} + \nabla^2$ is the d’Alembertian operator. In a vacuum ($T_{\mu\nu}=0$), this becomes the wave equation $\Box \bar{h}_{\mu\nu} = 0$, rigorously demonstrating that spacetime distortions are waves propagating at the speed of light $c$. Furthermore, by adopting the Transverse-Traceless (TT) gauge, the physical degrees of freedom are restricted to only two independent polarization modes, $h_+$ and $h_\times$.

Rigorous Derivation of the Quadrupole Formula

$$ \bar{h}_{\mu\nu}(t, \vec{x}) = \frac{4G}{c^4} \int \frac{T_{\mu\nu}(t - |\vec{x} - \vec{x}'|/c, \vec{x}')}{|\vec{x} - \vec{x}'|} d^3x' $$

We perform a multipole expansion assuming the distance to the observation point $r = |\vec{x}|$ is sufficiently larger than the size of the wave source ($r \gg |\vec{x}'|$). By repeatedly using the law of conservation of energy and momentum $\partial^\nu T_{\mu\nu} = 0$, the spatial integral of the spatial components $T_{ij}$ can be converted into the second time derivative of the moments of the energy density $T_{00}$ (i.e., mass density $\rho c^2$).

$$ \int T_{ij} d^3x = \frac{1}{2} \frac{d^2}{dt^2} \int T_{00} x_i x_j d^3x $$$$ h_{ij}^{TT}(t, r) = \frac{2G}{c^4 r} \left[ \ddot{I}_{ij}(t - r/c) \right]^{TT} $$

The generation of gravitational waves necessitates a time variation in the deviation from spherical symmetry of a mass distribution (the quadrupole moment). Radiation from a monopole (due to mass conservation) or a dipole (due to momentum conservation, or because the time derivative of the dipole moment is the total momentum, which is conserved) is forbidden. The coefficient $\frac{2G}{c^4}$ is an extremely minuscule value of about $1.65 \times 10^{-44} \text{ s}^2/\text{kg m}$, which is the fundamental reason why the detection of gravitational waves became humanity’s ultimate challenge, taking over 100 years.

Are Gravitational Waves a Physical Reality or a Coordinate Artifact? Historical Debate and Feynman’s Sticky Bead Argument

Einstein himself harbored doubts about the existence of gravitational waves throughout his life. Although he theoretically predicted them in 1916, in 1936 he attempted to write a paper with Nathan Rosen arguing that “gravitational waves do not exist due to the non-linearity of general relativity” (he later recognized his error and corrected it after it was pointed out by reviewer Howard Robertson). Among physicists of the time, there was a fierce debate: “Are gravitational waves merely a mathematical artifact arising from the choice of coordinate system, carrying no physical energy?”

The decisive thought experiment that put an end to this debate was the “Sticky bead argument” presented by Richard Feynman at the Chapel Hill Conference in 1957. Imagine a stick with beads threaded onto it with some friction. When a gravitational wave passes, the expansion and contraction of spacetime in orthogonal directions (tidal forces in the TT gauge) will cause relative acceleration between the beads and the stick. Because there is friction, this movement will generate heat energy. Feynman elegantly argued that since physical energy in the form of heat is produced, gravitational waves must be a “physical reality” that reliably carries energy. Subsequently, Hermann Bondi and others mathematically rigorously proved that gravitational waves do carry energy.


Chapter 2: A Century from Indirect Evidence to Direct Detection

Even after the existence of gravitational waves became theoretically certain, their direct detection remained a pipe dream. However, astronomical observations provided “indirect evidence” of their existence first.

The Hulse-Taylor Binary Pulsar and Orbital Decay

In 1974, Russell Hulse and Joseph Taylor discovered the neutron star binary system “PSR B1913+16” using the Arecibo radio telescope. This binary system orbits its common center of mass with a period of about 7.75 hours. As a result of precisely observing the arrival timing of radio pulses from the pulsar over many years, they found that the orbital period was shortening (the orbit was decaying) by about 76 microseconds per year.

$$ P = \frac{G}{45c^5} \langle \dddot{I}_{ij} \dddot{I}^{ij} \rangle $$

Assuming a Keplerian orbit with orbital eccentricity $e$, the rate of change of the period $T$, denoted as $\dot{T}$, can be theoretically derived. The observed amount of orbital decay perfectly matched (within a 0.2% error margin) the “energy loss due to gravitational wave radiation” predicted by general relativity. This became the first indirect evidence for the existence of gravitational waves, and Hulse and Taylor were awarded the Nobel Prize in Physics in 1993 for this achievement.

The Illusion of Joseph Weber’s Resonant Bar Detectors

The first serious challenge toward direct detection was initiated in the 1960s by Joseph Weber at the University of Maryland. He used huge aluminum cylinders (Weber bars) 2 meters long, 1 meter in diameter, and weighing about 1.5 tons. The principle was that when a gravitational wave passes near the resonant frequency of the cylinder, minute elastic vibrations would be excited in the cylinder due to tidal forces.

In 1969, Weber announced that he had “detected gravitational waves,” shocking the global physics community. However, although other research institutions built similar resonant bar detectors for follow-up testing, no one could reproduce Weber’s signals. The thermal noise (Brownian motion) of the aluminum completely drowned out the tiny gravitational wave signals, meaning the sensitivity was absolutely insufficient with the technology of the time. Weber’s claims were ultimately rejected, but his passion and challenge became an important cornerstone that paved the way for subsequent laser interferometer detectors.

Chapter 3: The Extreme Engineering of the Michelson Interferometer and the Noise Budget Curve

Scientists who felt the limitations of resonant bar types placed Michelson interferometers using lasers at the center stage of gravitational wave detection. When a gravitational wave passes, spacetime has the property (as a tensor wave) of stretching in a specific direction and shrinking in the orthogonal direction. The interferometer captures this minute differential phase change of $L_x - L_y$. However, to achieve the target strain sensitivity $h \sim 10^{-21} - 10^{-22}$, it was necessary to push the interferometer’s “noise budget” down to the absolute limit.

LIGO’s 4km Arms and Fabry-Perot Cavities

The American LIGO (Laser Interferometer Gravitational-Wave Observatory) consists of L-shaped interferometers with 4km arms built in Hanford, Washington, and Livingston, Louisiana. However, even with an optical path length of 4km, the expected expansion and contraction of space due to gravitational waves $\Delta L = h \times L$ is an agonizingly small $10^{-18}$ meters (less than 1/1000th of a proton).

To capture this minuscule change, “Fabry-Perot cavities” are incorporated into LIGO’s arms. A partially transmissive mirror (ITM) and a fully reflective mirror (ETM) are placed at both ends of the arm, allowing the laser light to bounce back and forth on average hundreds of times (finesse $\mathcal{F} \approx 450$) within the arm. This significantly amplifies the phase shift by stretching the effective optical path length close to the scale of the gravitational wave’s wavelength. Furthermore, an extremely complex optical system known as a “Dual-Recycled Fabry-Perot-Michelson Interferometer” is constructed by adding a “Power Recycling Mirror (PRM)” that pushes light returning toward the source from the beam splitter back into the interferometer, and a “Signal Recycling Mirror (SRM)” that optimizes the bandwidth of the signal component.

Quantum Noise: The Dilemma of Shot Noise and Radiation Pressure Noise

What limits the sensitivity of the interferometer in the high-frequency band (> 200 Hz) is “Shot noise,” caused by the discreteness of photons. The phase uncertainty due to the Poisson fluctuation of the number of photons reaching the photodetector decreases inversely proportional to the square root of the laser power $P$ ($\Delta \phi \propto 1/\sqrt{P}$). Therefore, in LIGO, an initially stabilized Nd:YAG laser of several tens of watts is boosted to hundreds of kilowatts inside the interferometer by power recycling.

However, increasing the laser power in turn makes “Radiation pressure noise” prominent in the low-frequency band (< 50 Hz). The fluctuation in the reaction force when a massive amount of photons collides with the mirror causes the mirror to shake randomly. This increases proportionally to the square root of the laser power ($\Delta x \propto \sqrt{P}$).

These two noises are direct consequences of Heisenberg’s uncertainty principle regarding the mirror’s position and momentum, $\Delta x \Delta p \ge \hbar/2$, and the theoretical lower limit of sensitivity determined by their intersection is called the “Standard Quantum Limit (SQL)”. In the noise budget curve of a gravitational wave detector, the SQL forms an insurmountable V-shaped valley.

Surpassing the Standard Quantum Limit with Squeezed Light

Introduced to shatter this SQL is the pinnacle of quantum optics, “Squeezed vacuum states” (squeezed light). This is a technology that compresses (squeezes) the fluctuation that affects observation—either the “phase fluctuation” or the “amplitude (radiation pressure) fluctuation” of light—while sacrificing the other to satisfy the uncertainty principle.

A squeezed vacuum state generated by an optical parametric oscillator using a nonlinear optical crystal (OPO) is injected from the output port (dark port) of the interferometer. Furthermore, “Frequency-dependent squeezing” is implemented in Advanced LIGO’s latest upgrade (A+) and KAGRA. This technology uses a long filter cavity to rotate the angle of the squeezed light’s ellipse for each frequency, optimally compressing phase fluctuation at high frequencies and amplitude fluctuation at low frequencies. This has successfully reduced quantum noise below the SQL limit simultaneously across the entire frequency band.


Chapter 4: Vibration Isolation Engineering and the Fluctuation-Dissipation Theorem of Thermal Noise

Dominating the low-to-mid frequency band (10 Hz to 100 Hz) of the interferometer are physical disturbances on Earth, namely seismic noise and thermal noise. To achieve a precision of 1/10,000th of an atomic nucleus, these must be eliminated to the utmost limit.

Seismic Noise and the Transfer Function of Multi-Stage Pendulums

Minute ground vibrations (microseisms) have a spectral density of around $10^{-7}/f^2 \text{ m}/\sqrt{\text{Hz}}$ depending on the frequency $f$, which is more than ten orders of magnitude larger than a gravitational wave signal.

To isolate this seismic noise, LIGO employs passive vibration isolation using “Multiple-stage pendulums.” A single-stage pendulum acts as a low-pass filter, attenuating external disturbances proportionally to $(f_0/f)^2$ in the band above its resonant frequency $f_0$. LIGO’s end mirrors (test masses) are suspended by 4-stage pendulums (quad suspensions). As a result, the transfer function attenuates with a fierce steepness of $(f_0/f)^8$ at high frequencies.

Moreover, by combining multi-degree-of-freedom active damping systems using hydraulics and piezoelectric elements (which measure ground shaking with seismometers and apply anti-phase forces via feedforward and feedback control to cancel it out), vibrations from the ground are practically completely isolated in the band above 10Hz.

Thermal Noise and the Fluctuation-Dissipation Theorem

Even if vibration isolation is perfect, as long as matter is not at absolute zero, the atoms making up the mirrors themselves randomly vibrate due to thermal energy $k_B T$. This is called “Thermal noise.”

$$ S_x(f) = \frac{k_B T}{\pi^2 f^2} \text{Re} [Z(f)] \approx \frac{k_B T}{\pi f} \frac{V_0}{E} \phi(f) $$

Here, $Z(f)$ is the mechanical impedance of the system, $V_0$ is the effective volume, $E$ is Young’s modulus, and $\phi(f)$ is the mechanical loss angle of the material.

Particularly severe around 100Hz are “coating thermal noise” from the dielectric multilayer coating deposited on the mirror’s reflective surface, and “suspension thermal noise” from the fibers suspending the mirrors. LIGO dramatically reduces suspension thermal noise by monolithically welding and suspending high-purity fused silica mirrors with extremely low mechanical loss using silica fibers of the same material.

KAGRA: The Underground Environment of Kamioka Mine and Cryogenic Cooling of Sapphire Mirrors

The ultimate approach to further lower the thermal noise $S_x(f)$ is to lower the temperature $T$ itself. The Japanese Large-scale Cryogenic Gravitational Wave Telescope “KAGRA” chose this path.

KAGRA is the only facility in the world combining the following two innovative technologies:

  1. Low Seismic Noise in an Underground Environment: Built over 200m underground in the Kamioka Mine, Gifu Prefecture. Compared to the surface, the seismic background noise is extremely quiet, about 1/100th, which directly leads to improved sensitivity in the low-frequency band.
  2. Cryogenic Sapphire Mirrors: For the test masses, “single-crystal sapphire” was adopted, which has dramatically high thermal conductivity and extremely small mechanical loss $\phi(f)$ at low temperatures. This is cooled down to 20K (minus 253 degrees Celsius) using cryogenic refrigerators and ultra-fine pure copper heat links.

Cryogenic operation is considered an essential technology for next-generation (third-generation) gravitational wave telescopes (Einstein Telescope, Cosmic Explorer). While confronting immense technical difficulties unique to cryogenic temperatures, such as optical asymmetry due to sapphire’s birefringence, minute vibrations transmitted from the cooling system (noise injection through heat links), and residual gas adsorption on the mirror surface (frosting phenomenon), KAGRA plays a crucial role as a demonstrator pioneering the forefront of humanity’s capabilities.

Chapter 5: September 14, 2015 - The Full Picture of the Historic Detection GW150914 and the Mathematics of Waveform Analysis

The moment when a century of theoretical exploration and decades of extreme engineering challenges bore fruit arrived suddenly. On September 14, 2015, at 09:50:45 UTC, both Advanced LIGO detectors in Hanford and Livingston recorded the exact same waveform showing a magnificent match. This was “GW150914,” the first gravitational wave directly detected in human history.

Binary Black Hole Merger and Mass Defect

As a result of data analysis, it was revealed that this signal was emitted when two black holes with masses 36 and 29 times that of the Sun spiraled towards each other in space about 1.3 billion light-years away from Earth (redshift $z \approx 0.09$), finally merging into a single massive black hole of 62 solar masses.

What is noteworthy is the mass defect. It should have been 36 + 29 = 65, but the mass after the merger was 62 solar masses. Where did the energy corresponding to the lost “3 solar masses” go? Following Einstein’s $E=mc^2$, all of it was converted into pure gravitational wave energy and released into space. For a brief moment just before the merger, the peak luminosity of the gravitational waves emitted by this binary system reached approximately $3.6 \times 10^{49}$ watts ($\sim 200 \text{ M}_\odot c^2 / \text{s}$), a staggering figure that was over 50 times the total power output of all the light from all the stars in the observable universe.

Post-Newtonian Expansion of the Chirp Signal and Matched Filtering

The waveform of GW150914 was a typical “Chirp Signal”—a waveform whose frequency and amplitude rapidly increase with time.

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

Here, $\mathcal{M}$ is a parameter called the “Chirp mass,” and using the masses $m_1, m_2$ of the two black holes, it is defined as $\mathcal{M} = \frac{(m_1 m_2)^{3/5}}{(m_1 + m_2)^{1/5}}$. From the observed frequency change $\dot{f}$ of the gravitational wave, this chirp mass is directly read with extremely high precision (for GW150914, $\mathcal{M} \approx 30 M_\odot$).

This waveform is broadly modeled in three phases:

  1. Inspiral Phase: The stage where two black holes approach while orbiting each other. A waveform model that calculates the above post-Newtonian approximation to a very high order (such as 3.5PN) is applied.
  2. Merger Phase: The moment when event horizons touch and violently merge. Because the gravitational field becomes extremely strong and nonlinearities dominate, the waveform can only be predicted by Numerical Relativity using supercomputers.
  3. Ringdown Phase: The stage where the highly distorted Kerr black hole after the merger settles into a spherical (more accurately, oblate) shape while radiating excess energy as gravitational waves. It is described as Quasinormal modes based on black hole perturbation theory, resulting in an exponentially decaying sine wave.
$$ \rho^2 = 4 \int_0^\infty \frac{|\tilde{s}(f) \tilde{h}^*(f)|}{S_n(f)} df $$

Through massive parallel computing using millions of templates, GW150914 was detected with a decisive significance of an SNR of 24.

Testing General Relativity in Strong Gravitational Fields

GW150914 not only proved the “reality of binary black holes” for the first time, but also enabled the first “verification of general relativity under extreme strong gravitational field and high-dynamics environments.” The observed waveform from inspiral to ringdown matched the predictions of the Einstein field equations perfectly. An upper limit on the mass of the graviton ($m_g < 1.2 \times 10^{-22} \text{ eV}/c^2$) was established, the propagation speed of gravity was proven to match the speed of light, and extremely strict constraints were placed on alternative theories of gravity.


Chapter 6: The Dawn of Multi-Messenger Astronomy and the Future of Cosmology

The detection of gravitational waves is a monumental pillar of physics in its own right, but its true value lies in collaboration with other observation methods. Light, radio waves, X-rays, neutrinos, and gravitational waves. The era of “Multi-Messenger Astronomy,” which observes the same celestial phenomena from multiple angles using various “messengers,” has begun.

GW170817: Neutron Star Merger and Simultaneous Observation of Electromagnetic Counterparts

Its biggest highlight is “GW170817,” observed on August 17, 2017. This was not a black hole, but a gravitational wave from the merger of two neutron stars. Unlike a black hole merger, when neutron stars collide, a massive amount of matter (neutron-rich matter) is scattered into space, accompanied by intense electromagnetic radiation.

$$ -3 \times 10^{-15} < \frac{v_{GW}-c}{c} < +7 \times 10^{-16} $$

This result killed off in one fell swoop many modified gravity theories (such as some tensor-scalar theories) that predicted the speed of gravitational waves to be different from the speed of light, which had been proposed to explain dark energy.

Kilonovae and Deciphering the Origin of Heavy Elements (Gold and Platinum)

A few hours later, ground-based optical telescopes captured the glow of a “Kilonova,” the debris of the merger phenomenon. This is a phenomenon where the fragments of neutron stars emit light through radioactive decay while expanding. Detailed spectral observations confirmed that elements heavier than iron (r-process elements) were synthesized in massive quantities during the merger process.

Until then, the main origins of heavy elements in the universe, such as gold, platinum, and uranium, had long been shrouded in mystery (it was thought that supernova explosions alone could not provide enough neutron density to explain the quantities). The observation of GW170817 presented incontrovertible evidence that the gold and platinum that make our rings shine were created by the cosmic catastrophe of a “neutron star collision” in the distant past.

Early Universe Inflation and Primordial Gravitational Waves

One of the ultimate targets that gravitational wave astronomy sets its sights on is “Primordial Gravitational Waves.” The “Inflation theory” posits that the universe underwent exponential rapid expansion immediately after its birth, prior to the Big Bang. During this dramatic expansion, quantum fluctuations of space were stretched to macroscopic scales and are believed to have been fixed as tensor-type fluctuations shaking the entire universe—namely, primordial gravitational waves.

Primordial gravitational waves are expected to leave traces in the polarization pattern (B-mode polarization) of the Cosmic Microwave Background (CMB), and also drift through space as a direct Stochastic Gravitational-Wave Background. If these can be detected, it will be direct evidence for the inflation theory and the greatest key to unlocking the laws of quantum gravity in the extreme energy regime of particle physics (Grand Unified Theory, Planck scale).

Prospects for the LISA Space Telescope and Next-Generation Ground Detectors

Current ground detectors (LIGO, Virgo, KAGRA) target the frequency band from 10Hz to several kHz (mergers of stellar-mass black holes and neutron stars). However, the universe is teeming with even lower-frequency (slower period) gravitational waves. Examples include the mergers of supermassive black holes of millions to billions of solar masses at the centers of galaxies, and Extreme Mass Ratio Inspirals (EMRI) of compact stars.

To capture these, plans are underway to move beyond the limits of seismic noise on the ground and build huge interferometers in space. This is the “LISA (Laser Interferometer Space Antenna)” project, driven primarily by the European Space Agency (ESA). LISA is a space interferometer of an astounding scale, linking three spacecraft flying in a triangular formation 2.5 million kilometers apart in a heliocentric orbit with laser links (scheduled for launch in the mid-2030s). Its frequency band will be $10^{-4}$ Hz to $10^{-1}$ Hz, allowing it to comprehensively cover the merger history of supermassive black holes across the universe and close in on the mysteries of galaxy formation and evolution.

Simultaneously on the ground, concepts for third-generation detectors (Europe’s Einstein Telescope, America’s Cosmic Explorer) with arm lengths extending from 10km to 40km are progressing. If these are realized, they will be able to capture every black hole merger occurring out to the edge of the observable universe (redshift $z>10$).


Conclusion: From Einstein’s Legacy, and Beyond

The direct detection of gravitational waves was an extraordinary feat exactly 100 years after the theory’s prediction. It marks a historic turning point where humanity became able to not only “see” the universe but also “listen” to it.

Laser interferometers battling the fluctuations of microscopic quantum noise and thermal noise, quieting the tremors of the Earth, and capturing the extreme distortions of spacetime. Behind them lie the tenacity and wisdom of thousands of scientists and engineers over several generations. The profound emotion of the moment when the quadrupole formula and post-Newtonian expansion, which were merely strings of mathematical equations, perfectly matched the heartbeat of the real universe proves the profundity of physics as a discipline and the triumph of human intellect.

We are now merely standing at the entrance to gravitational-wave astronomy. The advancement of the international observation network by LIGO, Virgo, and KAGRA, the construction of next-generation ground detectors, and the launch of space interferometers including LISA. The symphony of multi-messengers played by gravitational waves, electromagnetic waves, and neutrinos will surely continue to tell us the deepest, most violent, and most beautiful secrets of the universe. Humanity, having finished solving the final homework left by Einstein, is now stepping powerfully towards the frontiers of an unknown cosmology that even Einstein could not have imagined.

comments powered by Disqus