The quadrupole radiation formula calculates the gravitational wave emission from accelerating mass distributions.The quadrupole radiation formula follows from Einstein's equations in the weak-field approximation. Gravity, like electromagnetism, can create waves, but with an important difference: an electric dipole (a pair of + and - charges) radiates dipole radiation, but in gravity, dipole radiation is forbidden because of the conservation of momentum. The first allowed term is quadrupole. It's like sound: a uniformly pulsating sphere is silent, but one that flattens and stretches makes noise. Einstein realized this in 1916, but only in the 1970s was the formula refined for astrophysical calculations.
The idea of gravitational waves came to Albert Einstein in 1916, but initially he derived the formula for the simplest case of a rotating dumbbell. In 1918, he corrected errors and showed that gravitational radiation is indeed of quadrupole nature. In the 1960s–70s, with the discovery of pulsars, the formula was refined and rigorously derived in the work of scientists such as Richard Feynman (explained through a 'sticky bead') and Kip Thorne (developed detection theory). Today it forms the foundation of gravitational-wave astronomy.
The formula is critically important for calculating the evolution of close binary systems consisting of neutron stars or black holes. It predicts the shape of the wave signal (chirp), by which LIGO and Virgo identify merger events. Its range of applicability is weak fields and slow motions (v << c); near black holes, a full numerical relativity calculation is needed. The formula also fails for isolated axisymmetric systems where the quadrupole does not change.
Thanks to it, astrophysicists calculate when black holes will merge, and LIGO detectors are tuned to search for exactly this pattern—without it, gravitational waves would be predicted but not recognized.
The famous binary pulsar PSR B1913+16, discovered by Hulse and Taylor, decays exactly according to the quadrupole formula, which earned them the Nobel Prize.
| h_{ij}^{\text{TT}} | variable | Transverse-traceless gravitational wave strain tensor — gravitational wave strain |
| D | variable | Distance from source to observer — distance metre |
| I_{ij} | variable | Mass quadrupole moment tensor — mass quadrupole moment kilogram metre squared |
| t | variable | Time — time second |
| G | constant | Constant of gravitation in Newton's law — gravitational constant = 6.67430e-11 cubic metre per kilogram per second squared |
| c | constant | Universal speed limit and speed of light in vacuum — speed of light = 299792458 metre per second |
| d^2/dt^2 | operator | Second derivative with respect to time — second time derivative |
| A | variable | Sunrise azimuth angle — azimuth angle radian |
| φ | variable | Latitude of observer — latitude radian |
| δ | variable | Solar declination — solar declination radian |
| arccos | operator | Inverse cosine function — arccosine |
| tan | operator | Tangent function — tangent |
| v_c | variable | Circular orbital speed — velocity metre per second |
| M | variable | Mass of the central body — mass kilogram |
| r | variable | Orbital radius (distance from center of mass) — radius metre |
| G | constant | Universal constant of gravitation, determining the strength of gravitational force. — gravitational constant = 6.67430e-11 cubic metre per kilogram per second squared |
| v_esc | variable | escape velocity — velocity metre per second |
| M | variable | mass of the gravitating body — mass kilogram |
| r | variable | distance from the center of mass — radius metre |
| G | constant | Universal constant of gravitation, determining the strength of gravitational force. — gravitational constant = 6.67430e-11 cubic metre per kilogram per second squared |















