HUNTING LOW-FREQUENCY CONTINUOUS GRAVITATIONAL WAVES WITH THE EINSTEIN@HOME DISTRIBUTED COMPUTING SYSTEM

Gravitational waves (GWs) are ripples in the gravitational field, predicted to exist by Einstein’s General Theory of Relativity. A neutron star is the extremely dense remnant of the core of a massive star. Born in a supernova explosion, neutron stars are some of the most interesting objects in the universe. They have the mass of about the sun but the size of a small city. All the forces of nature — gravity, the strong and weak nuclear interactions, and electromagnetism — play huge roles in their properties and behavior.

Rapidly rotating neutron stars are expected to emit GWs continuously as they spin, if they are not perfectly spherical. The amount by which a neutron star deviates from being perfectly spherical can be described by a parameter called the ellipticity. The larger the ellipticity, the stronger the GWs emitted by the neutron star. Neutron stars are expected to have very small ellipticities: the astrophysics of their structure predicts that they are among the most spherical objects in the universe, with ellipticities of 10-5 or lower. The best human technology produces spheres that are close to perfect at that level (for example, the quartz sphere gyroscopes that flew on Gravity Probe B). For normal neutron star material, the maximum deformation that the neutron star crust could sustain before collapsing is about 10 cm. Nevertheless, by analyzing a long stretch of data from the LIGO detectors with carefully prepared software, we could pick even such a very weak signal out of the detector noise. At present, a few thousand neutron stars in our Galaxy have been identified as pulsars by their radio, X-ray or gamma-ray emissions, but our Galaxy is estimated to contain roughly 100 million neutron stars that are not observed through their electromagnetic emissions. With a search for continuous GWs that scans the entire sky, we have the potential to make new detections of neutron stars and learn about their structure.

Figure 1: This plot shows the upper limits on the GW strain amplitude for this search (black curve) and for four other all-sky searches (red, blue, purple and green). The curves represent the source strain amplitude at which 90% of simulated signals could be confidently detected.

GWs from rapidly spinning neutron stars will be detected at some specific (but unknown) signal frequency, close to a perfect sine wave. We nominally expect the GW signal frequency to be twice the neutron star spin frequency, but it could also be at any multiple of the spin frequency. The spin frequency is the number of revolutions of the neutron star per second, measured in Hertz (Hz). These all-sky searches need to look for unknown sources located everywhere in the sky, for a wide range of possible signal frequencies and frequency derivatives, where the frequency derivative is the rate of change of the frequency. The ideal data analysis strategy used to extract the faint continuous wave signals from the detector noise is given by the matched filtering method, applied to many months of data. Matched filtering is a data analysis method in which the data are correlated against a simulated waveform in order to try to identify that signal buried in the detector noise. However, this data analysis technique requires a huge amount of computing power, more than is available to the LIGO Collaboration, when data stretches of the order of months or years are used and a wide fraction of the parameter space is searched over. Hence, to search for these signals using a realistic amount of computing power, we use a “hierarchical” approach in which the entire data set is broken into shorter segments and afterwards the information from the different segments is combined. The huge computational effort is distributed through the volunteer computing project Einstein@Home, which relies on members of the public to donate their idle computing power to GW searches. Hundreds of thousands of host machines, that contributed a total of approximately 25000 CPU (Central Processing Unit) years, were needed to carry out the computations. We can use segments of a few days, allowing for very sensitive searches. As a bonus, it provides an opportunity to involve the public in the search for GWs!

This all-sky search for GWs from isolated neutron stars uses data from the LIGO detectors (LIGO Hanford in Washington and LIGO Livingston in Louisiana) collected in 2015 – 2016 during the first Advanced LIGO run(O1). The Advanced LIGO detectors have undergone many upgrades to make them more sensitive than the Initial LIGO detectors. The sensitivity improvement is especially significant at low frequencies. In fact, the previous Einstein@Home all-sky search, using Initial LIGO data, did not look at frequencies below 50 Hz due to lack of detector sensitivity in that range. In this search, we focused the immense computing resources of the Einstein@Home computing project on an all-sky search for GWs from spinning neutron stars in the low frequency range of Advanced LIGO data, from 20 Hz to 100 Hz. Under the assumption that the GW frequency is twice the neutron star frequency, this means we are looking for neutron stars spinning between 10 and 50 revolutions per second. Our goal is to perform a deep all-sky search in a region much of which had not previously been explored with LIGO.

Figure 2: This plot shows the ellipticity of a source at different distances (indicated by the colors) that would have been detected by this search, with 90% confidence, if it was emitting continuous gravitational waves.

The search considers neutron stars with frequency derivative ranging from about − 26 × 10-10 Hz/s to 3 × 10-10Hz/s. We mainly consider negative frequency derivatives, as we would expect the spin frequency of the neutron star to decrease over time as it loses energy. For example, if a neutron star is spinning at 50 Hz and emitting GWs at 100 Hz, after 12 years the signal frequency will fall to as low as 99 Hz. The observed frequency derivative can differ from the actual frequency derivative due to radial acceleration of the neutron star, which is why we consider a small range of positive frequency derivatives.

No GW signals have been detected in this search, presumably because there are no neutron stars emitting strongly enough in our frequency band to be detected at our sensitivity level. We quantify that conclusion by calculating upper limits on the GW strain amplitude. The upper limit is the smallest GW signal amplitude that we would have been able to detect with a confidence of 90%. Since we have not detected any signal, we can exclude the existence of signals of the type that we have searched for with amplitudes greater than our upper limit value. For example, at the highest frequency of this search near 100 Hz, where the detector is most sensitive, we can exclude the presence of signals with a GW strain amplitude of 1.8 × 10-25 or higher with a 90% confidence level. We can only compute upper limits with a finite statistical confidence level because detector noise is random and thus has to be described using probability and statistics. As shown in Figure 1, the upper limits set by this analysis (black curve) are ~1.5 to 2 times more constraining than those set by other all-sky searches at the same frequency. This improvement is primarily thanks to the significant computing resources of the Einstein@Home project, but also partly due to the methods and post-processing techniques used in the Einstein@Home search. One should also keep in mind that the other all-sky searches cover a broader range of frequency and frequency derivatives. The upper limits set by this analysis show a factor of 5 to 10 improvement with respect to the best results from the initial LIGO detectors.

If you have been running Einstein@Home on your computer, we thank you very much for your donated computations! If you haven’t been participating, but want to, you can get started by going to https://einsteinathome.org.

SEE MORE: