Gravitational Wave Strain Amplitude Calculator
Calculate gravitational wave strain amplitude for compact binaries or rotating neutron stars. Uses the quadrupole approximation.
🌊 What is a Gravitational Wave Strain Amplitude Calculator?
Gravitational wave strain amplitude h is the dimensionless quantity that characterizes how strongly a gravitational wave distorts spacetime as it passes through a detector. Concretely, h = Delta-L / L, where Delta-L is the change in arm length and L is the original arm length. For GW150914, LIGO's 4 km arms moved by about 4 x 10^-18 m, less than one-thousandth the diameter of a proton. This calculator implements the two leading GW source models in the quadrupole approximation: compact binary inspiral (BBH, BNS, or BH-NS) and continuous emission from a rotating neutron star with non-zero equatorial ellipticity.
For compact binaries, the strain is h = (4/r)(GMc/c^2)(pi*G*Mc*f_GW/c^3)^(2/3), where Mc is the chirp mass -- the particular combination (m1*m2)^(3/5)/(m1+m2)^(1/5) that dominates the inspiral waveform. Chirp mass can be measured to sub-percent accuracy from the rate of frequency evolution alone, making it the most precisely determined parameter in GW astronomy. GW150914 (Mc ~ 28.3 M_sun, 410 Mpc) and GW170817 (Mc ~ 1.19 M_sun, 40 Mpc) are the landmark events that opened the era of multi-messenger GW astrophysics.
For rotating neutron stars, a non-axisymmetric mass distribution (ellipticity epsilon) generates continuous monochromatic gravitational waves at exactly twice the rotation frequency. The strain is h = 4*G*I*epsilon*(2*pi*f_rot)^2 / (c^4*r), where I is the moment of inertia. Even epsilon = 10^-6 (a mountain only ~0.1 mm high on a 10 km star) gives strain near LIGO's detection floor for pulsars within a few kpc. Current searches have placed astrophysically interesting upper limits on epsilon for dozens of known pulsars.
The GW luminosity -- the power radiated as gravitational waves -- reaches staggering levels during compact binary mergers. GW150914 briefly radiated ~3.6 x 10^49 W, more than 10 times the total electromagnetic luminosity of all stars in the observable universe. Even the quieter neutron star binary GW170817 emitted ~10^46 W during its final inspiral seconds. These numbers illustrate why LIGO's precision interferometers, capable of measuring motions smaller than an atomic nucleus, represent one of the most demanding engineering achievements in physics.
📐 Formula
📖 How to Use This Calculator
Steps
💡 Example Calculations
Example 1 - GW150914 Binary Black Hole Merger
GW150914: m₁ = 36 M☉, m₂ = 29 M☉, f = 35 Hz, r = 410 Mpc
Example 2 - GW170817 Binary Neutron Star Merger
GW170817: m₁ = 1.36 M☉, m₂ = 1.17 M☉, f = 100 Hz, r = 40 Mpc
Example 3 - Crab Pulsar Continuous GW (upper limit)
Crab Pulsar: M = 1.4 M☉, R = 10 km, ε = 10⁻⁴, f_rot = 30 Hz, r = 2 kpc
❓ Frequently Asked Questions
🔗 Related Calculators
What is gravitational wave strain amplitude h?
The strain amplitude h is the fractional change in length Delta-L/L induced by a gravitational wave as it passes through a detector. It is dimensionless: h = 10^-21 means two mirrors 4 km apart oscillate by about 4 x 10^-18 m. LIGO measures strain by laser interferometry with sensitivity near 10^-23 per root-Hz in its most sensitive frequency band.
What is the formula for GW strain from a compact binary?
For a compact binary (BBH, BNS, or BH-NS) in circular orbit the leading-order quadrupole formula gives h = (4/r)(GMc/c^2)(pi*G*Mc*f_GW/c^3)^(2/3), where Mc is the chirp mass, f_GW is the gravitational wave frequency (twice the orbital frequency), and r is the luminosity distance to the source.
What is chirp mass and why is it important?
Chirp mass Mc = (m1 m2)^(3/5)/(m1+m2)^(1/5) is the combination of component masses that controls both the strain amplitude and the rate of frequency evolution (chirp rate df/dt = (96/5) pi^(8/3) (G Mc/c^3)^(5/3) f^(11/3)). It can be measured extremely precisely from the observed chirp even when individual masses are uncertain. GW150914 had Mc = 28.3 M_sun.
How does distance affect gravitational wave strain?
Strain h falls as 1/r with luminosity distance, just like electromagnetic flux falls as 1/r^2 in power (because h is an amplitude, not a flux). Doubling the distance halves the strain. This means LIGO's detectable volume grows as r^3 with sensitivity, and a factor-of-10 strain improvement in the detector opens a factor-of-1,000 larger survey volume.
What is the formula for GW strain from a rotating neutron star?
A neutron star with non-zero equatorial ellipticity epsilon = (Ixx - Iyy)/Izz emits continuous gravitational waves at twice its rotation frequency. The strain is h = 4 G I epsilon (2 pi f_rot)^2 / (c^4 r), where I is the moment of inertia (canonical value ~10^38 kg m^2 for a 1.4 M_sun, 10 km NS). Even epsilon = 10^-6 (a tiny mountain) gives detectable strain from nearby pulsars.
What is GW luminosity and how large is it?
For compact binaries, GW luminosity L = (32/5) c^5/G * (G Mc pi f / c^3)^(10/3). At peak, GW150914 radiated about 3.6 x 10^49 W, briefly outshining all stars in the observable universe combined. Even the quieter GW170817 (BNS) emitted roughly 10^46 W at 100 Hz. For rotating neutron stars the luminosity is much lower, typically 10^20 to 10^30 W.
Why is GW frequency twice the orbital or rotation frequency?
Gravitational waves are quadrupole radiation. A binary system has two arms of the mass quadrupole rotating at the orbital frequency f_orb; because the quadrupole completes a full cycle of symmetry twice per orbit, GW emission occurs at 2 f_orb. Similarly, a non-axisymmetric rotating star presents the same quadrupole orientation twice per rotation, giving f_GW = 2 f_rot.
What is the LIGO sensitivity needed to detect these signals?
LIGO Advanced (O3 run) achieved a design sensitivity noise floor of about 3 x 10^-24 per root-Hz near 100 Hz. A signal with h ~ 10^-21 stands well above the noise with SNR above 10 in a 1-second integration. At 410 Mpc (GW150914 distance), only extremely massive mergers are detectable; third-generation detectors like Einstein Telescope aim for sensitivity of 10^-25 per root-Hz.
What values of neutron star ellipticity are physically realistic?
Theoretical models suggest the NS crust can support mountains with epsilon up to roughly 10^-7 to 10^-5 depending on the breaking strain of nuclear pasta. Strong internal magnetic fields (magnetars) can distort the star to epsilon ~ 10^-4 to 10^-3. Current LIGO continuous-wave searches have set upper limits of epsilon below 10^-7 for many known millisecond pulsars.
What are the Einstein Telescope and Cosmic Explorer?
Einstein Telescope (ET) and Cosmic Explorer (CE) are proposed third-generation ground-based GW detectors planned for the 2030s. ET aims for 10 km triangular underground arms and would improve strain sensitivity by a factor of 10 over Advanced LIGO, detecting BBH mergers to redshift z ~ 20. CE proposes 40 km L-shaped surface arms with even greater reach.
How is strain measured from the example event GW150914?
GW150914 was a BBH merger of 36 and 29 solar masses at 410 Mpc. The observed peak strain was h ~ 1.0 x 10^-21 at about 150 Hz. This calculator returns h ~ 1.12 x 10^-21 at 35 Hz (an earlier inspiral point); the peak near merger is higher because the frequency and chirp mass are higher. Use the GW150914 preset to explore the parameter space.