Einstein@Home

Einstein@Home is a volunteer computing project that searches for signals from spinning neutron stars in data from gravitational-wave detectors, from large radio telescopes, and from a gamma-ray telescope. Neutron stars are detected by their pulsed radio and gamma-ray emission as radio and/or gamma-ray pulsars. They also might be observable as continuous gravitational wave sources if they are rapidly spinning and non-axisymmetrically deformed. The project was officially launched on 19 February 2005 as part of the American Physical Society's contribution to the World Year of Physics 2005 event.
Einstein@Home searches data from the LIGO gravitational-wave detectors. The project conducts the most sensitive all-sky searches for continuous gravitational waves. While no such signal has yet been detected, the upper limits set by Einstein@Home analyses provide astrophysical constraints on the galactic population of spinning neutron stars in our Milky Way galaxy.
Einstein@Home also searches radio telescope data from the Arecibo Observatory, and has in the past analyzed data from Parkes Observatory. On 12 August 2010, the first discovery by Einstein@Home of a previously undetected radio pulsar J2007+2722, found in data from the Arecibo Observatory, was published in Science. This was the first data-based discovery by a volunteer computing project. As of December 2023, Einstein@Home had discovered 55 radio pulsars.
The project also analyses data from the Fermi Gamma-ray Space Telescope to discover gamma-ray pulsars. On 26 November 2013, the first Einstein@Home results of the Fermi data analysis was published: the discovery of four young gamma-ray pulsars in data from Fermi's Large Area Telescope. As of December 2023, Einstein@Home has discovered 39 previously unknown gamma-ray pulsars in data from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. The Einstein@Home search makes use of novel and more efficient data-analysis methods and discovered pulsars missed in other analyses of the same data.
The project runs on the Berkeley Open Infrastructure for Network Computing software platform and uses free software released under the GNU General Public License, version 2. Einstein@Home is hosted by the Max Planck Institute for Gravitational Physics and the University of Wisconsin–Milwaukee. The project is supported by the Max Planck Society, the American Physical Society, and the US National Science Foundation. The Einstein@Home project director is Bruce Allen.
Einstein@Home uses the power of volunteer computing in solving the computationally intensive problem of analyzing a large volume of data. Such an approach was pioneered by the SETI@home project, which is designed to look for signs of extraterrestrial life by analyzing radio wave data. Einstein@Home runs through the same software platform as SETI@home, the Berkeley Open Infrastructure for Network Computing. As of December 2023, more than 492,000 volunteers in 226 countries had participated in the project, making it the third-most-popular active BOINC application. Users regularly contribute about 7.7 petaFLOPS of computational power, which would rank Einstein@Home among the top 105 on the TOP500 list of supercomputers.

Scientific objectives

The Einstein@Home project was originally created to perform all-sky searches for previously unknown continuous gravitational-wave sources using data from the Laser Interferometer Gravitational-Wave Observatory detector instruments in Washington and Louisiana, USA. The best understood potential CW sources are rapidly spinning neutron stars which are expected to emit gravitational waves due to a deviation from Rotational symmetry. Besides validating Einstein's theory of General Relativity, direct detection of gravitational waves would also constitute an important new astronomical tool. As most neutron stars are electromagnetically invisible, gravitational-wave observations might also reveal completely new populations of neutron stars. A CW detection could potentially be extremely helpful in neutron-star astrophysics and would eventually provide unique insights into the nature of matter at high densities, because it provides a way of examining the bulk motion of the matter.
Since March 2009, part of the Einstein@Home computing power has also been used to analyze data taken by the PALFA Consortium at the Arecibo Observatory in Puerto Rico. This search effort is designed to find radio pulsars in tight binary systems. It is expected that there is one radio pulsar detectable from Earth in an orbital system with a period of less than one hour. A similar search has also been performed on two archival data sets from the Parkes Multi-beam Pulsar Survey. The Einstein@Home radio pulsar search employs mathematical methods developed for the search for gravitational waves.
Since July 2011, Einstein@Home is also analyzing data from the Large Area Telescope, the primary instrument on Fermi Gamma-ray Space Telescope to search for pulsed gamma-ray emission from spinning neutron stars. Some neutron stars are only detectable by their pulsed gamma-ray emission, which originates in a different area of the neutron star magnetosphere than the radio emission. Identifying the neutron star's rotation rate is computationally difficult, because for a typical gamma-ray pulsar only thousands of gamma-ray photons will be detected by the LAT over the course of billions of rotations. The Einstein@Home analysis of the LAT data makes use of methods initially developed for the detection of continuous gravitational waves.

Gravitational-wave data analysis and results

Einstein@Home has carried out many analysis runs using data from the LIGO instruments. Since its first search run in 2005, the sensitivity of the LIGO detectors has been improved in a series of steps and upgrades. This is continuing with the current Advanced LIGO detectors. At the same time, Einstein@Home search algorithms have also improved. Together these have increased the search sensitivity by several orders of magnitude.
Einstein@Home's first analysis used data from the "third science run" of LIGO. Processing of the S3 data set was conducted between 22 February 2005 and 2 August 2005. This analysis employed 60 segments from the LIGO Hanford 4-km detector, totaling ten hours of data each. Each 10-hour segment was analyzed for CW signals by the volunteers' computers using a matched-filtering technique. When all matched-filtering results were returned, the results from different segments were then combined in a "post-processing step" on Einstein@Home servers via a coincidence scheme to further enhance search sensitivity. Results were published on the Einstein@Home webpages.
Work on the S4 data set was started via interlacing with the S3 calculations and finished in July 2006. This analysis used 10 segments of 30 hours each from the LIGO Hanford 4-km detector and 7 segments of 30 hours each from the LIGO Livingston 4-km detector. Besides the S4 data being more sensitive, a more sensitive coincidence combination scheme was also applied in the post-processing. The results of this search have led to the first scientific publication of Einstein@Home in Physical Review D.
Einstein@Home gained considerable attention in the international volunteer computing community when an optimized application for the S4 data set analysis was developed and released in March 2006 by project volunteer Akos Fekete, a Hungarian programmer. Fekete improved the official S4 application and introduced SSE, 3DNow! and SSE3 optimizations into the code improving performance by up to 800%. Fekete was recognized for his efforts and was afterward officially involved with the Einstein@Home team in the development of the new S5 application. As of late July 2006, this new official application had become widely distributed among Einstein@Home users. The app created a large surge in the project's total performance and productivity, as measured by floating point speed, which over time has increased by approximately 50% compared to non-optimized S4 applications.
The first Einstein@Home analysis of the early LIGO S5 data set, where the instruments initially reached their design sensitivity, began on 15 June 2006. This search used 22 segments of 30 hours each from the LIGO Hanford 4-km detector and six segments of 30 hours from the LIGO Livingston 4-km detector. This analysis run, employing the search methodology as Einstein@Home, was very similar to the previous S4 analysis. However, the search results were more sensitive due to the use of more data of better quality compared to S4. Over large parts of the search parameter space, these results, which also appeared in Physical Review D, are the most exhaustive published to date.
The second Einstein@Home search of LIGO S5 data constituted a further major improvement regarding search sensitivity. As opposed to previous searches, the ensuing results were already combined on the volunteers' computers via a Hough transform technique. This method matched-filtered results from 84 data segments of 25 hours each, parameters from which came from both 4-km LIGO Hanford and Livingston instruments.
On 7 May 2010, a new Einstein@Home search, which uses a significantly improved search method, launched. This program analyzed 205 data segments of 25 hours each, using data from both 4-km LIGO Hanford and Livingston instruments. It employed a technique which exploited global parameter-space correlations to efficiently combine the matched-filtering results from the different segments.
Results from an Einstein@Home all-sky search for continuous gravitational waves in LIGO S5 data were published on 13 February 2013. In the most sensitive frequency band of the search, the presence of periodic gravitational waves with strain amplitude larger than 7.6×10⁻²⁵ could be excluded at 90% confidence. Overall, the search was 3 times as sensitive as previous Einstein@Home searches in LIGO S5 data. Details of the two-stage follow-up procedure for signal candidates used in this study were published on 25 June 2014.
A search for high-frequency continuous gravitational waves in LIGO S5 data by Einstein@Home, published on 26 September 2016, was the only such search in LIGO data. No signal candidates were identified. The search excluded neutron stars with spin frequencies between 625 Hertz and 770 Hertz and with ellipticities greater than 2.8×10⁻⁷ closer than 100 parsec to Earth.
Data from LIGO 6th science run were analyzed by Einstein@Home and the results were published on 18 November 2016. No signal was found and the search set the most stringent upper limits for an all-sky search for continuous gravitational waves at the time of publication. In the most sensitive frequency band between 170.5 Hertz and 171 Hertz there were no continuous gravitational waves with a strain amplitude of more than 5.5×10⁻²⁵ detected. At frequencies of 230 Hertz, the search results exclude neutron stars with ellipticities greater than 10⁻⁶ within 100 parsecs of Earth.
Einstein@Home conducted a directed search for continuous gravitational waves from the central object in the supernova remnant Cassiopeia A. It used data from the LIGO S6 run and searched over a range of frequencies from 50 Hertz to 1000 Hertz, because the spin frequency of the central object is unknown. No signal was found. The upper limits on gravitational-wave emission from Cassiopeia A were the most stringent at the time of publication, about a factor two lower than previous upper limits.
On 28 December 2016 results from a follow-up of the all-sky search for continuous gravitational waves in LIGO S6 data were published. Out of a total of 3.8 × 10¹⁰ signal candidates from the earlier search, the 16 million most promising were analyzed using a four-stage hierarchical process. No candidate was found to be consistent with an astrophysical source of continuous gravitational waves. In the frequency band between 170.5 Hertz and 171 Hertz the upper limit on the strain amplitude was 4.3×10⁻²⁵, a factor 1.3 lower than in the previous search.
Searches for continuous gravitational waves are limited by the available computing power. Within the project, research on improving the sensitivity of the searches with new methods is conducted. In late 2017 two publications were published, describing improved methods of candidate clustering in the hierarchical searches and new "veto" methods that distinguish between astrophysical continuous gravitational waves and detector artifacts mimicking them.
Both these new methods were employed in the first Einstein@Home all-sky search for continuous gravitational waves in Advanced LIGO data from the first observing run, the results of which were published on 8 December 2017. The first part of the search investigated the lower end of the LIGO frequency band between 20 Hertz and 100 Hertz. No signals were found. The most stringent upper limit on the gravitational-wave strain amplitude set by the search was 1.8×10⁻²⁵ at a frequency of 100 Hertz.
An Einstein@Home study on how to optimally use the limited computing power for directed searches was published on 31 January 2018. It describes the design of searches for continuous gravitational waves over a wide frequency range from three supernova remnants.
The results from the directed Einstein@Home search for continuous gravitational waves from the central objects of the supernova remnants Vela Jr., Cassiopeia A, and G347.3 was published on 29 July 2019. It covered a frequency range from 20 Hertz to 1500 Hertz and used data from LIGO's first observing run O1. No signal was found and the most stringent upper limit at the time of publication were set, improving earlier results by a factor of two for all three targets.
A follow-up of the Einstein@Home search for continuous gravitational waves from the central objects of the supernova remnants Vela Jr., Cassiopeia A, and G347.3 was published on 29 June 2020. It investigated the most promising 10,000 candidates from the previous search and followed them up in two stretches of data from LIGO's second observing run. A single candidate associated with G347.3 remained as a possible signal after the follow-up, but was not conclusively confirmed based on gravitational-wave data. Archival X-ray data were searched for pulsations at the putative rotation frequency of the neutron star and its integer multiples. No signal was found. It is expected that data from LIGO's third observing run will suffice to shed light on the nature of this potential candidate.
On 8 March 2021 results from an Einstein@Home all-sky search for continuous gravitational waves in LIGO O2 data were published. It used an eight-stage follow-up process and covered a frequency range from 20 Hertz to 585 Hertz and reached the highest sensitivity for any all-sky survey below 500 Hertz. Six candidates were found after all follow-up stages. They are consistent with and caused by validation hardware injections in the LIGO instruments. No other signal was found. The most stringent upper limit was set in a 0.5 Hertz band at 163 Hertz at a gravitational-wave strain amplitude of 1.3×10⁻²⁵. The results begin to probe neutron star astrophysics and population properties. They exclude neutron stars with rotation frequencies above 200 Hertz with ellipticities larger than 10⁻⁷ closer than 100 parsec.
Results from a dedicated Einstein@Home search for continuous gravitational waves from the central object of the supernova remnant G347.3 was published on 5 August 2021. In the analysed frequency range between 20 Hertz and 400 Hertz no signal was found. The derived upper limits correspond to ellipticities of less than 10⁻⁶ for most of the frequency band. In the most sensitive frequency band at 166 Hertz the upper limit on gravitational-wave strain is 7.0×10⁻²⁶.
In July 2023, the results of an all-sky search for continuous gravitational waves in the public LIGO O3 data were published. The search was the most sensitive at that time for gravitational waves with frequencies between 2o Hertz and 800 Hertz and with spin-downs of up to −2.6×10⁻⁹ Hz s⁻¹. No astrophysical gravitational-wave signal was identified, and all candidate signals could be attributed to artificial signals injected into the LIGO data for validation purposes. The results exclude the existence of isolated neutron stars spinning at rotational frequencies of more than 200 Hertz with ellipticities larger than 5×10⁻⁸ closer than 100 parsec.