# 6.4: Gravitational Wave Sources

- Page ID
- 45411

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\id}{\mathrm{id}}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\kernel}{\mathrm{null}\,}\)

\( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\)

\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\)

\( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)

\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)

\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vectorC}[1]{\textbf{#1}} \)

\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)In general, gravitational waves are produced by moving masses. This can include the collapse of stars or objects in orbit. One caveat, though, is that the motion can't be spherically symmetric. A perfectly spherical star, then, would not produce gravitational waves when it collapses or explodes, but an asymmetric star would. While perfectly spherical stars don't actually exist, the degree to which a typical star deviates from perfect spherical symmetry is small enough that wouldn't expect very strong gravitational waves to come from it. For this reason, the primary source that we can expect to detect gravitational waves from is objects in orbit. While many objects in space orbit around something, only the most massive objects have any hope of producing gravitational waves large enough for us to detect them. Furthermore, the orbit has to be fast enough to produce gravitational waves with a frequency on the order of 100 Hz. This vastly reduces the possible orbital scenarios that we can detect. So far, for example, we have only been able to detect binary systems (that is, two objects orbiting around each other) that consist of black holes and/or neutron stars.

Let's analyze a simple binary system consisting of objects with mass \(M_1\) and \(M_2\) separated by a distance \(r\). As long as \(r\) is much greater than the Schwarszchild radius and the orbital speed is non-relativistic, the total energy can be approximated using Newtonian mechanics. The result (see Box 6.4.1), in natural units, is

\[E=-\frac{M_1M_2}{2r}.\label{eq:orbitalEnergy}\]

The gravitational waves, meanwhile, carry energy away, which makes the total energy of the binary system *more negative*.

Use Equation \ref{eq:orbitalEnergy} to argue that, as the binary system loses energy, the distance between the orbiting objects decreases.

**Answer**-
Fig. 6.4.1 shows a plot of the total energy as a function of \(r\).

You can see from the graph that energy loss corresponds to smaller values of \(r\).

It can be shown (though we will not do so here), that the rate of energy loss is

\[\frac{dE}{dt}=-\frac{32}{5r^5}\left(M_1M_2\right)^2\left(M_1+M_2\right).\label{eq:GWEnergyLoss}\]

Equation \ref{eq:GWEnergyLoss} can then be used to determine the rate at which the distance between them decreases. The result (see Box 6.4.2) is

\[\frac{dr}{dt}=-\frac{64}{5r^3}M_1M_2\left(M_1+M_2\right).\label{eq:GWRadialRate}\]

This result shows two important things. First, \(\frac{dr}{dt}\) is negative, which confirms that the distance between the orbiting objects decreases with time. Second, the rate at which the distance decreases gets larger as the orbiting objects get closer together.

Finally, the orbital speed increases as the distance between them decreases (see Box 6.4.3), which increases the frequency of the emitted gravitational waves. The result is a **chirp**. That is, there is a very rapid increase in the frequency of the gravitational wave in the moments just before the orbiting objects merge. Fig. 6.4.2 shows what this looks like for an actual gravitational wave event at three different detectors (Virgo is another detector that came online after the original detection by Hanford and Livingston).

The point of this is that we have a theoretical expectation for what a gravitational wave signal should look like, and the properties of that signal depend on the masses of the orbiting objects, the distance between them, and the distance from us to the binary system. Gravitational wave detectors continuously look for signals with this characteristic shape. When they do find one, they perform a parameter fit, which gives us information about the binary system. Then by looking at the time delay between the signal at the different detectors, we can determine where on the sky on the signal came from.

The total energy in Equation \ref{eq:orbitalEnergy} is a combination of potential energy and kinetic energy. The potential energy, in natural units, is

\[PE=-\frac{M_1M_2}{r}\nonumber\]

while the kinetic energy is

\[KE=\frac{1}{2}\mu v^2,\nonumber\]

where \(\mu=\frac{M_1M_2}{M_1+M_2}\) is called the **reduced mass** (reduced mass is used when two objects revolve around a point between them rather than one object orbiting around another stationary object). Meanwhile, the gravitational force of attraction is given by

\[F=\frac{M_1M_2}{r^2}.\nonumber\]

Finally, the centripretal force is

\[F=\mu\frac{v^2}{r}.\nonumber\]

Use these to derive Equation \ref{eq:orbitalEnergy}.

Take the derivative of Equation \ref{eq:orbitalEnergy} with respect to time, then equate it to Equation \ref{eq:GWEnergyLoss} in order to prove Equation \ref{eq:GWRadialRate}. (Hint: the \(r\) in Equation \ref{eq:orbitalEnergy} is a function of time, so you must use the chain rule when taking the derivative.)

Use information from Box 6.4.1 to argue that orbital speed increases as the distance between the orbiting objects decreases.