Searching for gravitational waves from a nearby globular cluster
The detected signal which comprised the first direct detection of gravitational waves lasted around two-tenths of a second. However, these short, transient signals are far from the only type of gravitational waves that can be searched for. This publication is about a search for what are known as continuous gravitational waves. Rather than the "chirp" from the first detection, imagine a faint steady tone, of constant or very slowly changing pitch (i.e., frequency), ringing continuously. This type of gravitational wave signal, known as a continuous-wave (CW) signal, can be emitted by a rapidly rotating neutron star with some form of surface asymmetry – such as a mountain on its crust, or an oscillation (called an r-mode) produced on its surface.
The search described here is called a directed search: it targets a known sky location, but searches over a wide range of star rotation speeds (or rates). In other words, it knows exactly where it's looking, but makes as few assumptions as possible about what it might find there. In this way, the analysis has a good sensitivity to potential gravitational waves while its computational cost remains under control.
The target for this directed search is the globular cluster NGC 6544, pictured here. Globular clusters are collections of stars very tightly bound by gravity, giving them very high stellar densities – a large number of stars packed into a very small region. They are typically found in the halo of a galaxy, orbiting the galactic core. They are normally composed of very old stars, which would not be good candidates for gravitational waves because their asymmetries tend to disappear over time. However, the extremely high concentration of stars in a globular cluster could result in a high rate of interactions between stars — either with other stars or debris — which could spin a star back up and/or introduce new asymmetries. NGC 6544 is a globular cluster that is close — gravitational wave strain decreases with distance — and it has an especially high concentration of stars in its center.
One of the key features of this search is a new algorithm for making Doppler corrections, called barycentric resampling. A CW signal is shifted in frequency by the Doppler effect due to the motion of the Earth through space---both its daily rotation and its yearly journey around the Sun. Correcting for these Doppler shifts can be computationally intensive. The barycentric resampling method shifts each data sample in time by the light travel-time (remember that gravitational waves travel at the speed of light) between the Earth and the Solar System barycenter, which arranges the data as if it were being recorded by a detector located at the barycenter. The barycenter, unlike the Earth, does not move through space in the reference frame of our Solar System, so using this method to mimic a detector situated at the barycenter eliminates the Doppler shift and is more than ten times faster computationally. Although this search used data from the sixth initial LIGO science run (S6), this method has become standard for continuous wave searches and is already seeing extensive use in the Advanced LIGO/Virgo era.
The search found no gravitational wave signals, but was able to set what are called upper limits on the strength of possible signals. Given the sensitivity of our data, we can place a limit on the strain of gravitational waves and say that, if there had existed waves at or above the limit, our search would have detected them. These upper limits appear in Figure 2 from the paper, reproduced here. Based on energy conservation, as well as an assumed age of, and distance to, any possible sources in the globular cluster, we can set a limit on strain which represents the loudest possible signal we could see from the cluster in this search, represented by the black line. The computed upper limits lie below the loudest signal limit, which means that this analysis is scientifically relevant: it is sensitive to potentially real signals, weaker than the strongest possible one. These limits on gravitational wave strain translate into limits on the ellipticity — a measure of the spherical asymmetry — of possible neutron stars, and limits on the amplitude of the r-mode oscillations on their surfaces.
- Neutron star: When a large enough star has burnt all its nuclear fuel and reaches the end of its life cycle, the fusion processes at its core cease and the star collapses inward. With a mass comparable to that of our sun condensed into a star the size of a small city, neutron stars are held up against gravity by neutron degeneracy pressure--the effect of the Pauli exclusion principle not allowing the neutrons that make up the star to occupy the same point in space.
- r-modes: Waves on the fluid surface of a young, high-temperature neutron star. These r-modes have a frequency which is comparable to the spin frequency of the star, so for young neutron stars they could be in the LIGO and Virgo frequency band. See here for an informative animation.
- Asymmetry: Continuous gravitational waves are not produced by spherically symmetric objects. An asymmetry is some deviation from spherical symmetry and is required in order to produce continuous gravitational waves. Examples include mountains on the surface of an otherwise perfect sphere.
- Strain: The fractional stretching and squeezing of spacetime in response to a gravitational wave. E.g., the change in length of one of the LIGO detector's arms divided by the length of that arm.
- Solar System barycenter: The center of mass of our local solar system. The solar system barycenter is the point around which all bodies in the Solar System orbit, so signals received at that point would experience no Doppler shift in the frame of our Solar System.
- Initial LIGO: The first epoch of LIGO, spanning 2002 – 2010. The search described in this publication used data from late 2010.
- Advanced LIGO/Virgo: September 2015 – present. The second epoch of LIGO, using upgraded detectors which are (to date) about three times more sensitive than initial LIGO, with further improvements in sensitivity planned. LIGO has been joined by the European Virgo detector, which has been upgraded as well.
Globular cluster NGC 6544. Image credit: Wikimedia.
Figure from the publication
For more information on how this figure was generated and its meaning, see the preprint at arXiv.org.
Upper limits on gravitational wave strain are the red circles. The blue curve is an expected upper limit curve calculated from the noise spectrum of our data. (A running median has been applied to the noise spectrum, which causes some discrepancies when very narrow noise artefacts are suppressed). The black line is a limit on gravitational wave strain based on energy conservation as well as the age of, and distance to, any sources in the globular cluster; it represents the loudest possible signal we could hope for.