Gravitational Wave Searches for Aligned-Spin Binary Neutron Stars Using Nonspinning Templates
We study gravitational wave searches for merging binary neutron stars (NSs). We use nonspinning template waveforms towards the signals emitted from aligned-spin NS-NS binaries, in which the spins of the NSs are aligned with the orbital angular momentum. We use the TaylorF2 waveform model, which can generate inspiral waveforms emitted from aligned-spin compact binaries. We employ the single effective spin parameter to represent the effect of two component spins () on the wave function. For a target system, we choose a binary consisting of the same component masses of and consider the spins up to , We investigate fitting factors of the nonspinning templates to evaluate their efficiency in gravitational wave searches for the aligned-spin NS-NS binaries. We find that the templates can achieve the fitting factors exceeding only for the signals in the range of . Therefore, we demonstrate the necessity of using aligned-spin templates not to lose the signals outside that range. We also show how much the recovered total mass can be biased from the true value depending on the spin of the signal.
pacs:04.30.–w, 04.80.Nn, 95.55.Ym
Recently, three gravitational wave (GW) signals have been detected by the two Advanced LIGO detectors GW1 ; GW2 ; GW3 . According to the observing dates, those signals were named as GW150914, GW151226, and GW170104. Detailed analyses showed that they were emitted from merging binary black holes GW1PE1 ; GW1PE2 ; GW12PE . On the other hand, although GW signals from black hole-neutron star (NS) binaries or NS-NS binaries have not yet been detected GWBNS , those signals are also expected to be captured in the near future by the advanced detector network ALIGO ; AVirgo ; KAGRA (for an overview of NS-NS merger rates, see Aba10 ; GWBNS ). They will tell us more information about the nature of NSs, such as equation of states and mass distribution Lat07 ; Pra13 ; Lee11 ; Lee14 .
In this work, we investigate fitting factors of nonspinning waveform templates for GW signals from aligned-spin NS-NS binaries based on the Advanced LIGO detectors. The fitting factor represents the best-match of template waveforms to a signal waveform Apo95 . It is generally used to evaluate the efficiency of templates in GW searches. For example, if , that means a loss of detection rates of , because the detection rate is proportional to the cube root of the signal-to-noise ratio (), , and . The studies on the efficiency of nonspinning templates for aligned-spin signals have been carried out in the past works Pri14 ; Dal14 ; Cap16 ; Cho16b , but those works have only considered binary black holes or black hole-NS binaries. Several authors have also performed such a study towards the precessing-spin signals Bro12 ; Aji14 ; Har14 ; Dal15 . Here, we focus on aligned-spin NS-NS binaries with the spins of and calculate fitting factors and biases using nonspinning templates. In particular, we give a detailed interpretation for the physical property of the results.
For merging NS-NS binaries, the GW frequencies are most likely to get out of the sensitivity frequency band of Advanced LIGO before the systems reach the ultra-relativistic regime. Therefore, the post-Newtonian approximation is valid for modeling the gravitational waveforms emitted from those binaries Bla14 . Here, we adopt a simple post-Newtonian waveform model, TaylorF2, in which the wave function is defined in the Fourier domain.
Ii GW search analysis
In signal processing, the matched filtering can be the most efficient method for signals of known shape buried in stationary Gaussian noise. Since the inspiral waveforms emitted from merging NS-NS binaries can be modeled almost accurately by the post-Newtonian approximation, the matched filter can be used in the GW search analysis. The match between a detector data stream , consisting of a GW signal and stationary Gaussian noise, and a waveform is defined as
where the tilde denotes the Fourier transform of the time-domain waveform, is the power spectral density (PSD) of the detector noise, and the low frequency cutoff () can be determined by the PSD curve. In this work, we use the zero-detuned, high-power noise PSD of Advanced LIGO apsd and assume Hz. If the normalized waveform is defined as , the signal-to-noise ratio can be obtained by .
In order to describe aligned-spin binary systems in circular orbits, the wave function should incorporate five extrinsic parameters (luminosity distance of the binary, two angles defining the sky position of the binary with respect to the detector, orbital inclination, and wave polarization), two arbitrary constants (coalescence time and coalescence phase ), and four intrinsic parameters (two masses and two spins). On the other hand, the extrinsic parameters only scale the wave amplitude, hence do not affect our analysis Cho14 . Therefore, we do not consider those extrinsic parameters in this work.
We define the overlap () by the match between a normalized signal and a normalized template Cho13 :
Here, the maximization over and can be easily performed by using certain analytic techniques All12 . Since we consider only nonspinning waveforms as templates, the overlaps can be distributed in the two-dimensional mass parameter space as
The fitting factor (FF) is defined as the best-match between and a set of Apo95 :
where represents the physical parameters considered in the template space. The connection of the fitting factor to the overlap can be given as
If we use a complete template waveform model that can produce exactly the same shape as the signal waveform, we obtain the fitting factor equal to 1, and the parameter values recovered by the templates are the same as the true values. However, typically, a waveform model cannot be complete; thus, the recovered values are likely to be biased from the true values systematically, and the fitting factor should be lower than 1. Therefore, in our overlap distributions, the recovered masses can be biased because we use nonspinning template waveforms towards the spinning signals. The bias can be easily determined by the distance from the true value () to the recovered value ():
The bias corresponds to a systematic error in the GW parameter estimation. As the efficiency of a template waveform model for the search is evaluated by the fitting factor, its validity for the parameter estimation can be examined by the bias Cut07 ; Fav14 ; Cho15d .
We determine the fitting factor and the bias in the following way Cho15c ; Cho15d ; Cho16b . First, we repeat a grid search around until we find the crude location of in the overlap surface. Note that, we use the chirp mass and the symmetric mass ratio , where , because those parameters are much more efficient than the component masses in our analysis Cho15a . Next, we estimate the size of the contour , where is the maximum value in the overlap surface. Finally, we find (almost) the exact location of by performing a grid search in the region of and choose the overlap value at that location as the fitting factor.
The waveform function of TaylorF2 is given as
where is the wave amplitude that consists of the masses and the extrinsic parameters. Since we use the normalized waveforms, the amplitude does not affect our analysis. The wave phase is expressed as
where is given by the post-Newtonian expansion with the expansion parameter . The coefficients of the phase equation are expressed as functions of and the two component spins with , being the spin angular momentum of the compact object. We consider the post-Newtonian expansion up to 3.5 pN order where the spin terms are included up to 2.5 pN order Aru09 .
For a target system, we choose a NS-NS binary consisting of the same component masses of . We consider the spins up to , which corresponds to the spin of the fastest-spinning millisecond pulsar observed so far Hes06 . Note that, the spins of known pulsars in double neutron star systems are below 0.04 Kra09 . In order to produce our template waveforms, we vary only the mass parameters fixing the spins to be . We then calculate two-dimensional overlap distributions by using Eq. (3) and obtain fitting factors and biases by using Eqs. (5) and (6), respectively.
In Fig. 1, we show the fitting factors of nonspinning templates (cf. Fig 6 of Aji11a ). One can easily find that the fitting factor contours are symmetric about the line of , and they are coincide with the lines of constant total spin. Generally, in aligned-spin systems, the effect of the two component spins can be represented by a single effective parameter Poi95 . Several definitions for the single parameter have been introduced by Refs. Dam01 ; Aji11a ; Aji11 ; San10 ; Kha16 We adopt the simplest one defined by , which has been introduced to model the phenomenological template families San10 ; Aji11 and used in the recent studies similar to this work Cho16b ; Cho17 . For equal mass systems, the lines of constant total spin are coincide with the lines of constant effective spin. Therefore, the pattern of the contours shows that the fitting factors mainly depend on the effective spin rather than the component spins.
On the other hand, the fitting factors show a significant discrepancy between the positive and the negative effective spins. Typically, a positively (negatively) aligned-spin signal tends to be recovered by a nonspinning template having a larger (smaller) value of than the true value of , i.e., , and the size of the bias () depends on the spin of the signal. On the other hand, the parameter value of is physically restricted to the range of (Note that the unphysical value of implies complex-valued masses), and this restriction results in a rapid fall-off of the fitting factor when reaches the physical boundaries (for more details, see Dal14 ; Cho16b ). In the region of , as the magnitude of the effective spin increases, decreases but does not reach the boundary in our spin range because this boundary is sufficiently far from the true value (). Therefore, the fitting factor gradually decreases with increasing . However, in the region of , is always equal to , i.e., because , and consequently the fitting factors have much lower values. In conclusion, the nonspinning templates are efficient for GW searches, i.e., , for equal mass NS-NS binaries only with the spins of . Note that, if we choose a binary having sufficiently asymmetric masses so that does not reach the boundaries and in the regions of and , respectively, then the fitting factors gradually decrease in both the regions (e.g., see Fig. 2 of Cho16b ).
One can infer how the recovered total mass can be biased by understanding the orbital motion of spinning binary systems. When the spin is negatively aligned with the orbital angular momentum, the spin-orbit coupling makes the binary’s phase evolution slightly faster, hence hastens the onset of the plunge phase, as compared to its nonspinning counterpart Cam06 . For that reason, the negative effective spin decreases the length of the waveform, as compared to the nonspinning case, and such a waveform best matches the one produced by a higher mass nonspinning binary. In a positively aligned-spin system, the spin-orbit coupling makes exactly the opposite effect. Therefore, a positive (negative) aligned-spin binary can be recovered by a lower (higher) mass nonspinning template. In Fig. 2, we show the biases for the total mass of the system, . As described above, one can find the positive biases in the region of . However, there are almost no biases in the region of . This is due to the constraint on the space. Typically, as increases from 0, the position of tends to move to the upper left hand side from the position of in the plane. However, for equal mass binaries, since is always unbiased due to the physical boundary in this region, its companion is also unbiased although the space is unrestricted here.
Iv summary and discussion
We studied the efficiency of nonspinning templates in GW searches for aligned-spin NS-NS binaries. We assumed that the target binary system had the same component masses of and calculated fitting factors and biases of nonspinning templates for aligned-spin signals with the spins in the range of . We found that both the fitting factor and the bias strongly depend on the effective spin () rather than the component spins, and we confirmed that unequal mass binaries also show a similar behavior. In particular, we found a significant discrepancy between the positive and the negative effective spins. In the region of , as the magnitude of the effective spin increases, the fitting factor decreases gradually, and the bias () increases steadily up to at . On the contrary, in the region of , the fitting factors have much lower values due to the effect of the physical boundary of , and almost no biases are shown for all the signals in that region.
We demonstrated that nonspinning templates can achieve the fitting factors exceeding for aligned-spin NS-NS binaries in the range of ; thus, they are efficient for GW searches only for the signals in that range. In order not to lose GW signals outside that range, one should take into account the spin parameters in the template waveforms. The detection pipeline of Advanced LIGO for NS-NS binaries has been using aligned-spin templates, where the magnitude of the component object’s spin is limited to GWBNS .
HSC was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIP) ( No. 2016R1C1B2010064 and No. 2015R1A2A2A01004238). CHL was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIP) (No. 2015R1A2A2A01004238 and No. 2016R1A5A1013277). This work used the computing resources at the KISTI Global Science Experimental Data Hub Center (GSDC).
- (1) B. P. Abbott et al., (LIGO Scientific Collaboration and Virgo Collaboration) Phys. Rev. Lett. 116, 061102 (2016).
- (2) B. P. Abbott et al., (LIGO Scientific Collaboration and Virgo Collaboration)Phys. Rev. Lett. 116, 241103 (2016).
- (3) B. P. Abbott et al., (LIGO Scientific Collaboration and Virgo Collaboration)Phys. Rev. Lett. 118, 221101 (2017).
- (4) B. P. Abbott et al., (LIGO Scientific Collaboration and Virgo Collaboration)Phys. Rev. Lett. 116, 241102 (2016).
- (5) B. P. Abbott et al., (LIGO Scientific Collaboration and Virgo Collaboration) Phys. Rev. X 6, 041014 (2016).
- (6) B. P. Abbott et al., (LIGO Scientific Collaboration and Virgo Collaboration) Phys. Rev. X 6, 041015 (2016).
- (7) J. Abadie et al., (LIGO Scientific Collaboration and Virgo Collaboration) Class. Quantum Grav. 27, 173001 (2010).
- (8) B. P. Abbott et al.(LIGO Scientific Collaboration and Virgo Collaboration), Astrophys. J. 832, L21 (2016).
- (9) J. Abadie et al.(LIGO Scientific Collaboration), Class. Quantum Grav. 32, 074001 (2015).
- (10) F. Acernese et al., Class. Quantum Grav. 32, 024001 (2015).
- (11) Y. Aso et al.(The KAGRA Collaboration), Phys. Rev. D 88, 043007 (2013).
- (12) J. M. Lattimer and M. Prakash, Phy. Repts. 442, 109 (2007).
- (13) M. Prakash, arXiv:1307.0397 (2013).
- (14) C.-H. Lee, H.-S. Cho, Y. M. Kim and H.-Y. Park, J. Korean Phys. Soc. 59, 2118 (2011).
- (15) C.-H. Lee and H.-S. Cho, Nuclear Physics A 928, 296 (2014).
- (16) T. A. Apostolatos, Phys. Rev. D 52, 605 (1995).
- (17) T. Dal Canton et al., Phys. Rev. D 90, 082004 (2014).
- (18) S. Privitera et al., Phys. Rev. D 89, 024003 (2014).
- (19) C. Capano, I. Harry, S. Privitera, and A. Buonanno, Phys. Rev. D 93, 124007 (2016).
- (20) H.-S. Cho, Phys. Rev. D 94, 124045 (2016).
- (21) D. A. Brown, I. Harry, A. Lundgren, and A. H. Nitz, Phys. Rev. D 86, 084017 (2012).
- (22) P. Ajith, N. Fotopoulos, S. Privitera, A. Neunzert, N. Mazumder, and A. J. Weinstein, Phys. Rev. D 89, 084041 (2014).
- (23) I. W. Harry, A. H. Nitz, D. A. Brown, A. P. Lundgren, E. Ochsner, and D. Keppel, Phys. Rev. D 89, 024010 (2014).
- (24) T. Dal Canton, A. P. Lundgren, and A. B. Nielsen, Phys. Rev. D 91, 062010 (2015).
- (25) L. Blanchet, Living Rev. Relativity 17, 2 (2014).
- (26) Advanced LIGO anticipated sensitivity curves, https://dcc.ligo.org/LIGO-T0900288/public.
- (27) H. -S. Cho and C. -H. Lee, Class. Quantum Grav. 31, 235009 (2014).
- (28) H. -S. Cho, E. Ochsner, R. O’Shaughnessy, C. Kim, and C. -H. Lee, Phys. Rev. D 87, 024004 (2013).
- (29) B. Allen, W. G. Anderson, P. R. Brady, D. A. Brown and J. D. E. Creighton, Phys. Rev. D 85, 122006 (2012).
- (30) C. Cutler and M. Vallisneri, Phys. Rev. D 76, 104018 (2007).
- (31) M. Favata, Phys. Rev. Lett. 112, 101101 (2014).
- (32) H.-S. Cho, Class. Quantum Grav. 32, 235007 (2015).
- (33) H.-S. Cho, Class. Quantum Grav. 32, 215023 (2015).
- (34) H.-S. Cho, J. Korean Phys. Soc. 66, 1637 (2015).
- (35) K. G. Arun, A. Buonanno, G. Faye, and E. Ochsner, Phys. Rev. D 79, 104023 (2009).
- (36) J. W. T. Hessels, S. M. Ransom, I. H. Stairs, P. C. C. Freire, V. M. Kaspi, and F. Camilo, Science 311, 1901 (2006).
- (37) M. Kramer and N. Wex, Class. Quantum Grav. 26, 073001 (2009).
- (38) P. Ajith, Phys. Rev. D 84, 084037 (2011).
- (39) E. Poisson and C. M. Will, Phys. Rev. D 52, 848 (1995).
- (40) T. Damour, Phys. Rev. D 64, 124013 (2001).
- (41) L. Santamaria et al., Phys. Rev. D 82, 064016 (2010).
- (42) P. Ajith et al., Phys. Rev. Lett. 106, 241101 (2011).
- (43) S. Khan et al., Phys. Rev. D 93, 044007 (2016).
- (44) H.-S. Cho, J. Korean Phys. Soc. 70, 735 (2017).
- (45) M. Campanelli, C. Lousto, and Y. Zlochower, Phys. Rev. D74, 041501 (2006).