DoNeutS:
A Double Neutron Star Catalogue
Gravitational wave (GW) emissions are ripples in the spacetime manifold caused by the movement of massive objects, and were predicted by Einstein's theory of General Relativity (GR). The most intense sources of GW emissions are the merger of compact objects like neutron stars and black holes, realeasing a great amount of energy. These components lose energy while orbiting each other, and as a consequence they spiral-in reducing the orbital period with time. The discovery of PSR B1913+16 (Hulse & Taylor, 1975) provided the first evidence for the existence of this effect, by the observation of orbital decay during time, \( \dot{P}_b \) (Taylor & Weisberg, 1982 ). This discovery triggered a great interest in the scientific community and placed DNSs as the main target for groundbased interferometric detectors.
Since the produced signals are very weak, detection of GW emssions was a challenging task for decades. The Laser Interferometer Gravitationl-Wave Observatory (LIGO) was founded in 1992 by Kip Thorne and Ronald Drever from Caltech and Rainer Weiss from the MIT (Abbott et al., 2009). LIGO consists of two identical interferometers separeted by 3002 kilometers, located in fairly isolated areas of Washington (LIGO Hanford) and Lousiana (LIGO Livingston). The reason behind this large distance is found in their sensitivities. Since these instruments are very sensitive, they can detect the tiniest vibrations on Earth, from sources that are nearby until those which are thousands of kilometers away, and this can mask the GW signal in each detector. Locating the interferometers far away, it is possible to distinguish between GW signals and local noise. The original instruments were operational from 2002 to 2010, but no detection was made at that time. During these years the team made huge progress in detector engineering and both interferometers were redesigned from 2010 to 2014.
LIGO's interferometers are, in a sense, similar to Michelson interferometers, invented in the 1880's to (un)probe the existence of ether through the detection of interference patterns, with an ideal configuration to precisely measure phase changes of a wave traveling along two perpendicular arms with an "L-shape", since GW's cause the manifold to stretch in one direction and compress in a perpendicular direction. The layout is shown in the figure below. In summary, a powerful laser beam is equally splitted and directed towards both arms, bouncing off mirrors at the end of each arm. An additional mirror is placed in each arm, near the beam splitter and 4 km away from the end mirror, to create a Fabry Perot cavity. This causes the laser to bounce between the two mirrors about 300 times before being recombined. When gravitational waves pass through, an interference pattern is observed in the recombined signal. The recombined laser beams are directed towards photodetectors that measure this interference pattern created by the supperpositon of the two beams. By analyzing the interference it is possible to determine the characteristics of the GW and gain insights into the astrophysical events that produced them.
The Virgo interferometer, with a 3 km arms and built outside of Pisa in Italy, joined the observations of LIGO in 2017, after a period improving its sensitivity (Accadia et al., 2012). Its operation, combined with both LIGO detectors is crucial to localize the GW sources in the sky, since the localization is determined from the time delay between different observatories. More recently, the KAGRA observatory also started operating (Aso et al., 2013). It is constructed in Japan, inside the Kamioka mine, where is also located the neutrino detector Super Kamiokande. The underground location and use of cryogenic systems took this detector to another level, minimizing seismic and molecular vibrations that can mask the GW signal.
The first detection of a gravitational wave occurred in 2015 from a binary black hole merger, GW150914 (Abbott, B. P. et al, 2016 ), by the Advanced LIGO detectors. However, the true era of multimessenger and multiwavelength astronomy started two years latter with the first detection of the compact binary merger GW170817 (Abbott, B. P. et al, 2017 ), quickly associated with a short gamma-ray burst (GRB) emitted from the same sky location, \(1.7~s\) after the merger (Goldstein et al., 2017 ). The precise location of the event allowed a follow-up, leading to detections in X-ray, ultraviolet, optical, infrared and radio bands. The optical counterpart, AT2017gfo (Coulter et al., 2017), placed NGC 4993 as the host galaxy with a Tully-Fisher distance of \( \sim40~Mpc \) (Freedman et al., 2001 ), consistent with the inferred GW luminosity distance. The association of the gravitational wave with the sGRB, as well as the observations across the electromagnetic spectrum by different groups, the total and component masses of the system and the offset related to the host galaxy center provided strong evidences for the DNS nature of the GW170817 progenitor.
In April 2019, the second gravitational wave signal consistent with the coalescence of a DNS was detected, GW190425 (Abbott, B. P. et al., 2020 ). This event is identified as single-detector since at the time it was detected the LIGO-LHO was offline and the signal-to-noise ratio of Virgo was 2.5, below the threshold of 4.0 for significance estimation. The localization of a GW signal relies on measuring the time delay between observatories, but since GW190425 is a single-detector event, was not possible to set straight constraints on the sky map region. As a consequence, a follow-up became challenging and until now there is now firm detection of an electromagnetic counterpart. Nonetheless, it is important to note that the event itself is highly significant. While it is not possible to rule out the possibility of one or both components being BH's, due to the absence of sensitivity to probe matter effects, this hypothesis is unlikely.
Parameters that influences the evolution of waveform phase can be precisely constrained given the large number of cycles in the observable waveform (Cutler & Flanagan, 1994 ). The same does not hold for parameters which affects only wave polarization or amplitude. The best constrained quantity from a GW signal is the called chirp mass (\(\mathcal{M} \)).
To lowest order, it is possible to describe the gravitational radiation from the quadrupole Newtonian formalism of two dimensionless particles in an almost circular orbit (assuming that gravitational damping circularizes the orbit). The angular orbital frequency is then given by: \begin{equation} \Omega_b^2 = \frac{GM_{tot}}{a^3}, \end{equation} where \(G\) is the gravitational constant, \(M_{tot}\) is the system total mass and \(a\) is the orbital separation. The loss of energy leads to a decrease in \(a\) and, consequently, a decrease in orbital period, since \(P_{b} \equiv 2\pi /\Omega_b\). Combining the inspiral rate of system and the equation above, the signal frequency (\(\nu=1/P\)) is then obtained: \begin{equation} \nu^{-11/3} \frac{d\nu}{dt} = \frac{96}{5}\pi^{8/3}\mathcal{M}^{5/3}. \end{equation}
Therefore, from the observed signal frequency and its derivative is possible to accurately measure the chirp mass, which is a combination of total and reduced (\(\mu\)) masses of the system, \(\mathcal{M} = \mu^{3/5} M_{tot}^{2/5}\), and to determine the rate at which the GW frequency changes over time. These calculations are valid for the source-frame. When detected on Earth the signal is Doppler shifted due to relative motion, so the detected value will be \(\mathcal{M}_{det} = (1 + z)\mathcal{M}\), where \(z\) is the redshift of the source. All binaries with the same \(\mathcal{M}\) will have the same chirping sweep. Still, post-newtonian corrections to Eq. () do introduce a dependence on the individual masses of the system. This quantity is inferred by a matched filtering technique, where the wave's phase evolution is matched with model templates predicted under GR.