Model-independent test of the parity symmetry of gravity with gravitational waves

Gravitational wave (GW) data can be used to test the parity symmetry of gravity by investigating the difference between left-hand and right-hand circular polarization modes. In this article, we develop a waveform-independent method to decompose the circular polarizations of GWs produced during the inspiralling stage of compact binaries, with the help of stationary phase approximation. Applying this method to the mock data, we measure the arrival time difference between left-hand and right-hand polarizations, to test the parity symmetry of gravity by constraining the velocity birefringence of GWs. If a nearly edge-on binary neutron-stars with observed electromagnetic counterparts at 40 Mpc is detected by the 2nd-generation detector network, the arrival time difference can be constrained at the accuracy of $O(0.1s)$ for the GWs at $f\sim 100$Hz, and one could derive the model-independent test on the parity symmetry in gravity: lower limit velocity difference between two circular polarizations can be constrained to be within $\mathcal{O}(10^{-17})$, and the lower limit of the energy scale of parity violation can be constrained within $\mathcal{O}(10^4{\rm eV})$.

Introduction Although Einstein's General Relativity (GR) has been considered to be the most successful theory of gravity since it was proposed, it faces the difficulties in both theoretically (e.g. singularity, quantization, etc), and observationally (e.g. dark matter, dark energy, etc). Therefore, testing GR in various circumstance is an important topic since its birth [1][2][3][4]. The discovery of GW compact binary coalescence source GW170817 [5], and its electromagnetic (EM) counterparts in different frequency bands [6], opens the new window of multi-messenger GW astronomy, which also provides an excellent opportunity to test GR in the strong gravitational fields [7][8][9][10][11]. Numerous works have been carried out, which can be separated into two classes: One is the model-dependent method. For a specific theory of gravity, one calculates the GW waveforms, and constrains its deviation from that of GR (see, for instance, [12]). Due to the large number of modified gravity models, the modelindependent test becomes more attractive. One example is the parameterized post-Einsteinian framework developed by Yunes et al. [13].
Parity symmetry implies that a directional flipping to the left and right does not change the laws of physics. It is well known that nature is parity violating. Since the first discovery of parity violation in weak interactions [14], the experimental tests become more necessary in the other interactions, including gravity. In most previous works, test of parity symmetry in gravity has focused on the Chern-Simon modified gravity (see for instance the review article [15]). Many parity-violating gravities with different action forms have been proposed for different motivations [16] (see Appendix A for the details). If the parity symmetry is violated, it is expected that a GW behaves asymmetrically in its two cir-cular polarization modes [17]. The observable effects include the "amplitude birefringence" [18][19][20][21][22] and the "velocity birefringence" [23][24][25][26]. Therefore, reconstructing circular polarizations from observed GW signals is crucial. In this article, we develop a waveform-independent method to reconstruct the circular polarizations of GWs with the help of stationary phase approximation (SPA). By mocking the potential observation of the upcoming 2nd-generation GW detector network, we test the reliability of the method, and apply it to constrain modelindependently the velocity birefringence of GWs. The foremost advantage is that this method is simple (without tedious calculation), clean (with least assumption), independent of GW waveform, and is applicable to the existing detector network. Circular polarizations of GWs In general, it is convenient to describe GWs with complex oscillating functions h s (t) [27]. The detector response is a linear combination of the real part of two wave polarizations [28], i.e. d I (t + τ I ) = s=+,× F s I h R s (t) + n I (t + τ I ), where I = 1, 2, 3, · · · labels the I-th detector. τ I is the relative time delay with respect to a reference time t taken at the center of the Earth. F s I are the detector's beampattern functions. For a network with I ≥ 2, in principle, an unbiased estimator for polarization mode h R s (t) can be solved directly from the data d I by introducing the Moore-Penrose psudo-inverse matrix A, which is composed of the detectors' antenna beam-pattern functions F s I [29,30]. For a given source direction, the time-delay corrected data from a network of N d GW detectors can be written in the frequency domains as  [29], where σ 2 if is the noise variance of i-th detector at the fth frequency bin. Throughout this paper, tilde denotes the quantity in the frequency domain. The estimator for the vector {h R f } is given by the Moore-Penrose inverse [29], i.e.h By the inverse Fourier transformation, the estimators for the GW signal in time domain h R s (t) can be derived. The circular polarizations of GWs are defined by the complex h s (t) reconstructed from the observables h R s (t). In Appendix B, we demonstrate that this reconstruction is achievable, if considering the SPA, which is applicable for the GWs produced during the inspiralling stage of coalescing compact binaries. Consider h + as an example, it can be shown that the Fourier components of h + (t) and h R + (t) satisfy the following relations, that is, one could deriveh + (f ) fromh R + (f ). This plays a crucial role in the method proposed in this article, that reconstructs the GW circular polarizations from observation. However, we should emphasize that, this conclusion is applicable only for the inspiraling stage of GW events, where SPA is appropriate. The basic reason for this reconstruction to work is that, in SPA, in each small time span, the amplitude of a GW can be considered as a constant, it is therefore possible to derive both amplitude and phase information of GWs, instead of the combination of them in a general case. Equivalently, GWs can be decomposed as left-hand (L) and right-hand (R) circular polarizations, which are defined ash R/L = (h + ± ih × )/ √ 2 [27]. Using the inverse Fourier transformation, we can also obtain the time-domain function h R (t) and h L (t), which are both complex functions.
In this decomposition, the antenna beam-pattern functions of the I-th detectors F s I should be known in advance, which depend on the sky direction of the GW signal (RA,Dec), and the polarization angle ψ s [31]. The former is assumed to be obtained from its EM counterparts, while the latter remains unknown without a template fitting. On the other hand, a change of ψ s corresponds to the rotation of the celestial sphere frame coordinate system along z-axis [31]. Considering a rotation of coordinate system along z-axis by an angle ψ, the left-hand and right-hand modes in the new coordinate system are given by h L/R = h L/R e ±2ψ [27], which means that this rotation of coordinate system is completely equivalent to a change of the phase of h R and h L with the same value but opposite signs. For this reason, if we care only about the phase evolution, instead of their exact values of phases, the polarization angle ψ s can be arbitrarily chosen in the reconstruction, as confirmed in the simulation analysis.
Note that for GWs emitted by the coalescence of compact binaries, the amplitude ratio of left-hand and righthand polarization modes is determined by the inclination angle ι of the GW event. For the face-on sources with ι = 0 • , or 180 • , GW is circularly polarized, i.e. only left-hand or right-hand polarization exists. While for the edge-on sources with ι = 90 • , the amplitude of two circular polarizations are comparable.
We test the reliability of this decomposition method by simulations. In our analysis, GR is considered as the fiducial theory of gravity. We consider the a GW signal produced by the coalescence of BNSs, which are accompanied by observable EM emissions in various frequency bands. Without lose of generality, we assume a GW event produced by the coalescence of BNS located at (RA=0, Dec=0). Similar to GW170817, the luminosity distance is adopted as d L = 40 Mpc. Both masses of NSs are chosen as 1.4M and their tidal deformability parameter are assumed to be 425, referring to Fig.5 of Ref. [5]. For the inclination angle ι, we consider both the face-on case with ι = 0 • and edge-on case with ι = 90 • . Since we are only interested in GWs in the inspiraling stage, Tay-lorF2 model is used to calculate GW waveforms in time domian.
To mimic the realistic case, we consider the detector network including AdvLIGO [32], AdvVirgo [33], and KAGRA [34], as their designed sensitivities. The gaussian noise for each detector is randomly generated using the PyCBC package [35]. Mock data set for each detector is a linear combination of signal and noise.
FIG. 1: Time-frequency representations of left-hand (left panels) and right-hand (right panels) polarization signals reconstructed from the simulated GW data of signal + noise for the 2nd-generation detector network. Upper panels show the results of a face-on GW event, and the lower panels show those of an edge-on event.
Following the procedure described above, we reconstruct the circular polarizations of GWs for two extreme cases of the inclination angles. Similar to LIGO/Virgo Collaboration (LVC) in [5,36], we present the results in time-frequency representations of the reconstructed GW signals [5,36], as shown in Fig.1. We have applied the Q-transformation to the reconstructed complex strain serials, and optimized the Q value in the range of (100,110)Hz to maximize signal-to-noise ratio (SNR) in the diagram [37]. For the case with face-on burst, only the left-hand polarization mode, manifests itself as the significant signal in the upper left panel in Fig.1. While for the case with edge-on source, in the lower panels of Fig.1 we observe both polarizations with similar power as expected.
Testing the chirality of gravity In the general parityviolating gravity, GWs propagating in the Friedmann-Robertson-Walker universe satisfy (see Appendix A for the details), where prime denotes the derivative with respect the conformal time η, and H ≡ a /a. a is the scale factor, and the present value is adopted as a 0 = 1. k is the wave number, which relates to the GW frequency by f = k/2πa. Nonzero ν A and µ A represent the deviation from GR [38][39][40]. Term ν A determines the amplitude evolution of GWs, and the term µ A represents velocity of GWs. The case with ν R = ν L means different dampings of left-hand and right-hand polarizations, which is the effect of "amplitude birefringence" [18]. For an individual frequency, the effect of amplitude birefringence completely degenerates with the value of inclination angle [15], which is unknown in advance for the model-independent measurement. Therefore, in principle, this effect can be tested by comparing the amplitude ratio of left-hand and righthand polarizations among different frequencies, which is left as a future work. The case with µ R = µ L = 0 represents the violation of Lorentz symmetry in gravity, which has been tightly constrained by comparing the arrival times of GW170817 and GRB170817a [41]. In this paper, we consider only the parity-violating case. µ R = µ L means the velocities of GW polarizations are different, that is, there exists the "velocity birefringence". In general, µ A can be parameterized as µ A = αρ A (k/aΛ) β , with ρ R = 1 and ρ L = −1, Λ is the energy-scale of the theory, α and β are the coefficients, which depend on the theory of modified gravity. As shown in Appendix A, in general α is a function of the cosmic time. In the local universe, we can ignore the time-dependence, and treat α as a constant, which is absorbed by the definition of Λ as discussed below. Since the measurements of GWs using laser interferometers are not sensitive to the GW amplitude, in this article we restrict attention to the effect of velocity birefringence. The parametrization of µ A can be equivalently written as the modified dispersion relation which follows the group velocity of GWs, i.e.
Note that since the sign of v A /c − 1 is determined by ρ A , ρ R and ρ L have opposite signs. If one polarization mode is superluminal, then the other is subluminal. For a given GW signal at redshift z emitting both left and right circular polarizations, different propagation velocities will result in the difference in signal arrival times, which is given by [42] . ∆t e is the emitting time difference of the two modes. For the comparison of GW with their EM counterparts, the uncertainty of ∆t e is the main problem for the GW velocity measurement, which strongly depends on the theoretical models of GW and EM bursts. For instance, LVC has assumed that the emitting time difference between GW signal of GW170817 and EM signal of GRB170817a is smaller than 10s predicted in some models [41], and use this to constrain the Lorentz symmetry of gravity [43]. Therefore, the constraint of GW velocity by comparing the arrival times of GW and EM signals is model-dependent. Fortunately, we do not need to make such assumptions in our method. Since circular polarizations of GWs are the spin 2 modes, independent of theory of gravity, both modes are produced by the instantaneous acceleration of mass quadrupole of the systems. For this reason, for a fixed wavenumber, i.e. k R = k L , we have ∆t e = 0, and arrival time deference becomes This formula gives a direct relation between the arrival time difference and the energy scale of parity violation Λ. For the specific case with β = 1, i.e. the gravities with lowest order of parity-violating terms (see Appendix A), the constraint of t R−L directly corresponds to the constraint on the energy scale Λ of parity violation. In realistic observations, the differences of arrival times of the two GW polarizations are observable, as long as the left-hand and right-hand polarization can be reconstructed. For the cases with low noise level, e.g. for the 3rd-generation detector network, the SNR for each frequency channel can be large enough, and we can read out the arrival times of both modes from time-frequency representation, and calculate directly the time difference t R−L . However, in the case with the 2nd-generation detectors, the signal for individual frequency channel is too noisy, we have to combine the data within a frequency bin (f lower , f upper ) to amplify the SNR. To realize it, for each polarization mode, we calibrate all frequencies in the bands to compensate the emitting-time difference for different frequencies. To the lowest Newtonian order, the emitting time of GW at frequency f is given by [44] t c − t = 2.18s (1.21M /M c ) where t c is the time at which f becomes infinity and M c is the chirp mass of the binary. Utilizing this formula, for the simulated data, we calibrate all the frequency bands in the span (90, 110)Hz to f = 100Hz, and superimpose them to obtain the arrival times. The results are presented in Fig.2. From this figure, we obtain the arrival times of left-hand mode t L = −2.186 +0.085 −0.088 s and right-hand mode t R = −2.182 +0.088 −0.092 s, by measuring the FWHM of the signal in Fig.2. The difference between them is derived directly t R−L = 0.004 +0.122 −0.127 s. We observe that, the uncertainty of arrival time difference can be achieved at the level of O(0.1s) for this particular GW source configuration. Although these results are derived from an individual simulation, the stability of the conclusions is confirmed by repeating the analysis above but adopting different realizations of detector noises. It is important to emphasize, in our analysis, the uncertainties of t R and t L depend not only on the errors of arrival times in GW measurement, but also on the time resolution of the time-frequency representation in Fig.1. Therefore, the uncertainty of time difference t R−L derived above might not be the optimal result, which is expected to be significantly reduced if a better way can be used to read out the arriving times of GW signals. The value of |t R−L | can be directly translated to the constraint on parity-violating gravity. The analysis above indicates that if a nearly face-one BNSs at 40 Mpc observed by 2nd-generation network, |t R−L | 0.12s at frequency f ∼ 100Hz is expected to be obtained. The first constraint is for the velocity difference between polarization modes. Using the relation v R /v L −1 = t R−L /d L , and d L = 40 Mpc, we found that |v R /v L − 1| < 3.1 × 10 −17 . The constraint on the parity-violating energy scale Λ is also expected to be obtained. From the relation in Eq. (6) we have k/Λ = (|t R−L |/T β ) 1/β . As in general, we consider the case with β = 1, i.e. the gravity includes the lowest-order parity-violating terms in GWs. Considering the LCDM model with cosmological parameters H 0 = 70.0km/s/Mpc, Ω Λ = 0.7, Ω m = 0.3 [40], we have T β = 4.1 × 10 15 s, it follows that Λ > 1.4 × 10 4 eV.
GW170817 At this writing, GW170817 is the unique GW event with observed EM counterparts. Therefore, the reconstruction of two circular polarizations for this event is carried out and the results are presented in Appendix C. For a given GW event, to test the parity symmetry, the foremost requirement is that the amplitudes of left-hand and right-hand modes are comparable, which is mainly determined by the inclination angle ι. In GR, the amplitude ratio is given by |h L |/|h R | = [(1 + cos ι)/(1 − cos ι)] 2 . The condition |h L |/|h R | ∈ (1/3, 3) requires that ι ∈ (74 • , 106 • ). Therefore, only the nearly edge-on events can be used here. However, for the event GW170817, observation gives that constraint ι ≥ 152 • [5], which follows that |h L |/|h R | < 0.004. So, for this event, the right-hand mode is completely dominant, which is also proved in our analysis in Appendix C. For this reason, the method introduced above cannot be applied in GW170817.
Conclusions The advent of multi-messenger GW astronomy opens a new window to test the characteristics of gravity in strong gravitational-field. In this article, we develop a method to reconstruct the circular polarization modes of GWs, emitted during the inspiralling stage of compact binaries, with the help of the source information obtained from the observations of EM counterparts. By simulating the mock GW data, we test the reliability of the decomposition method for various cases.
By measuring the arrival time difference of left-hand and right-hand polarizations of a GW, one can modelindependently test the parity symmetry of gravity. For the general parity-violating gravities, since the GW observation is more sensitive to the arrival time rather than the amplitude, as an example of application, we test the velocity birefringence of GWs by means of mock data, and find that if a nearly edge-on BNS event at 40 Mpc is observed by the 2nd-generation GW detector network, the arrival-time difference can be constrained at the accuracy of O(0.1s) for the GWs at f ∼ 100Hz. It follows that the fractional velocity difference of two modes can be constrained at the level of O(10 −17 ). For the general theories of gravity with lowest-order parity-violating terms, this result implies the expected constraint on the energy scale of parity violating Λ O(10 4 ) eV. Although not exactly fair, it is instructive to compare this constraint with the existing constraint of Chern-Simons modified gravity [45]. We find that this constraint is 17 orders better than the existing constraint in solar system, which is Λ 10 −13 eV [46], and 6 orders better than that derived from binary pulsars, which is Λ 10 −2 eV [47]. We consider parity-violating gravity, whose action is of the form where R is the Ricci scalar, L PV is a parity-violating Lagrangian, and L other is the Lagrangian for scalar field, which may be coupled non-minimally to gravity, and all the other matters and modification terms of gravity.
In the flat Friedmann-Robertson-Walker universe, GW is the tensor perturbation of the metric, i.e.
where a(η) is the conformal scale factor, η is the conformal time and χ i is the comving coordinates. The quantity h ij stands for the GW perturbation, which we take to be transverse and traceless gauge, δ ij h ij = 0 and The equation of motion of GW is determined by the tensor quadratic action, which reads is the standard Lagrangian obtained from the Einstein-Hilbert term R. In the viewpoint of effective fields [17], the operatorsḣ 2 ij and (∂ k h ij ) 2 are the only quadratic operators with two derivatives. Note that, in principle this Lagrangian can be generalized by adding a time-dependent coefficient in each term, which might be caused by the Horndeski or Galileon scalar field, through the perturbation of extrinsic curvature of the spatial slices. Since they are not relevant to parity violation, we neglect them here. The first possible corrections to the tensor mode come from terms with three derivatives. There are only two possible terms [17], where ijk is the antisymmetric symbol. At least in the unitary gauge, the standard quadratic action is modified by the addition of [17] L where a dot denotes the derivative with respect to the cosmic time t, c 1 and c 2 are dimensionless coefficients, which could depend on time, and Λ is the scale that suppresses these higher dimension operators. The equation of motion for GWs can be derived by varying the quadratic action with respect to h ij . Decomposing the GW in the circular polarization basis, in the frequency domain the equation of motion can be written as [26]h where A = R or L, standing for the right-hand or lefthand polarization mode respectively, and Note that, Eq.(A6) is the unifying description for lowenergy effective description of generic parity-violating GWs. To our knowledge, all the known parity-violating theories of gravity in the literature, even if including more than three-derivative terms, can be casted into this form [16].
Chern-Simons (CS) modified gravity with Pontryagin term coupled with a scalar field corresponds to the case with c 1 = c 2 , i.e. ν A ∝ ρ A (k/aΛ) and µ A = 0 [18]. Due to the disappearing of µ A term, only the amplitude birefringence effect exists in CS modified gravity. This theory has higher-derivative field equations, which induces the dangerous Ostrogradsky ghosts. For this reason, CS modified gravity can only be treated as a lowenergy truncation of a fundamental theory.
Ghost-free parity-violating theories of gravity have recently been explored. One of the theories has the Lagrangian L PV (see Eqs.(3.1), (3.2) and (3.4) in [48]), which includes the scalar field and its first derivatives. GW in this theory corresponds to the case with nonzero functions of c 1 = c 2 [26,49].
Another ghost-free parity-violating theory contains second derivatives of the scalar field, and the Lagrangian L PV is given by Eqs.(3.12), (3.13), (3.14) and (3.18) in [48], which corresponds to the case with c 1 = 0, c 2 = 0 [26,49]. For both ghost-free theories, we have the nonzero ν A ∝ ρ A (k/aΛ) and µ A ∝ ρ A (k/aΛ), which indicates that both amplitude birefringence and velocity birefringence exist for GWs. The GWs in the general ghost-free parity-violating theory including higher derivatives are also investigated in the recent work [50].
Horava-Lifshitz (HL) gravity is power-counting renormalizable theory because of the presence of high-order spatial derivative operators [51,52]. The exclusion of high-order time-derivative operators, on the other hand, renders the theory unitary, whereby it is expected to be ultraviolet (UV) complete. In HL theory, the parity symmetry can be violated by introducing the third-and fifthorder spatial derivative operators. The theory with terms of third-order spatial derivatives, i.e. three-dimensional gravitational CS term, corresponds to c 1 = 0 and c 2 = 0, which is equivalent to ν A = 0, µ A ∝ ρ A (k/aΛ) [24,25].
Parity-violating extension of the symmetric teleparallel equivalent of GR theory is a non-Riemannian formulation of gravity, which allows for a wider variety of consistent extensions than the metric formulation of gravity theory, and has been studied explicitly in [53]. If considering only the lowest two-derivative parity-violating terms, the equation of motion of GW is given by Eq.(A6) with ν A = 0 and µ A ∝ ρ A (k/aΛ) −1 . However, if considering the three-derivative parity-violating terms, the field equation of GW corresponds to ν A ∝ ρ A (k/aΛ) and µ A ∝ ρ A (k/aΛ). In principle, there are non-vanishing terms of all orders in derivatives. However, as proved in [53], these higher-order terms are straightforward and somewhat trivial generalizations of second and thirdderivative gravity.
Another parity-violating theory is to consider the fifthorder derivative operator L [26], which corresponds to the equation of motion of Eq.(A6) with ν A ∝ ρ A (k/aΛ) 3 and µ A ∝ ρ A (k/aΛ) 3 . It is expected that, the parity-violating theories with other fifthorder or higher-order derivative operators can also be written in the form of (A6) with power-law form of ν A and/or µ A .
We observe that in all these theories, except for CS modified gravity, the velocities of left-hand and righthand polarizations are different, i.e., the effect of velocity birefringence exists. In the real observations of GW detectors, the observables are the real parts of h + (t) and h × (t). We consider h + as an example, which can be written as the following (general) form and the real part of this wave is where t r is the retarded time. The terms A and Φ represent the amplitude and phase of GWs respectively, which are both functions of t r . For GWs produced during the inspiralling stage of compact binaries, the conditions d ln A/dt dΦ/dt and |d 2 Φ/dt 2 | (dΦ/dt) 2 are satisfied. Therefore, we can use SPA to obtain the waveform in Fourier domainh R + (f ) as follows [44], where Ψ + (t) ≡ 2πf t r − Φ(t r ) − π/4. In this equation t * is defined as the time at which dΦ/dt = 2πf , and Ψ + (t * ) is the value of Ψ + (t) at t = t * . Similarly, assuming SPA is applicable, we can also obtain the h + (t) in Fourier domainh + (f ) as follows, andh By comparing Eqs.(B3), (B4), (B6), (B7), we obtain the following relation: which indicates that one can obtain the Fourier componenth + (f ) fromh R + (f ). However, we should emphasize that, this relation is applicable only to the GW produced during the inspiraling stage of compact binaries, where SPA is correct.
Once we obtainh + (f ) andh × (f ), the Fourier components of h R and h L are derived straightforwardly, Using the inverse Fourier transformation, we can also obtain the time-domain functions h R (t) and h L (t).
the data with total duration of 2048 seconds and sampling frequency 4096Hz. Note that, we do not use the AdvVirgo data, in which the SNR of GW signal is too weak.  [5]. The polarization angle can be randomly chosen as mentioned above, so we adopt it as ψ s = 0 • . Before analysis, we filter all time serials with a 20-500Hz bandpass filter to suppress large fluctuations outside the detector's most sensitive frequency band.
We repeat the reconstruction method as before, but here the mock data is replaced with the real data of LIGO detectors. Following the decomposition proceeding, we first construct the unbiased estimators forh R + andh R × , and translate them toh + andh × , respectively. Based on these, the left-hand and right-hand polarizations in frequency domain and time domain are derived straightforwardly. We show the results in the time-frequency representation, which are given in the upper panels of Fig.3. From both diagrams, we find the GW signal is weak (As expected, they are even weaker than those in simulations), and it seems difficult to identify them directly by eyes. Therefore, we should use the frequency superimposition technique to amplify signal and get the arrival times of both GW polarizations, where we adopted M c = 1.188M derived in GR framework to approximately displace the true value of chirp mass [5]. In the analysis, we scan the arrival times of the frequency bins with 5Hz in length from (20,25)Hz, (25,30)Hz to (195,200)Hz. For both polarizations, we observe the maximum SNR at the bin of (140, 145)Hz. Similarly, we change the frequency bins with 10Hz in length, and find the best one is at (135, 145)Hz. The arrival times of circular polarizations are presented in Fig.3 (low panels).
For the arrival time of right-hand polarization, we observe the significant peaks ( 4σ) for both frequency bins (140,145)Hz and (135,145)Hz, which follows the consistent results, i.e., t R = −0.911 +0.060 −0.062 s for the former bin, and t R = −0.901 +0.063 −0.063 s for the latter bin. Note that, for each bin, we have calibrated the arrival times of all the frequency bands in the bin to a fixed frequency at 140Hz. Since here we have used the blind method for the signal search, and consider the signal only in a small frequency bin, the SNR is much lower than that derived from template fitting (∼ 32σ). On the other hand, the SNR of arrival time for left-hand polarization is too low for both frequency bins. These results are consistent with what we anticipate: Since the inclination angle of GW170817 is ι ≥ 152 • [5], the amplitude ratio of circular polarization modes |h L |/|h R | < 0.004. Therefore, the observed GW signal of this event is completely dominated by the right-hand mode.