Arxiv:2009.03476V3 [Astro-Ph.IM] 4 May 2021 Counterpart Is Found, So That We Can Use Gws As Stan- Than Each Individual LISA Or Taiji Detector [28]

Sky localization of space-based gravitational wave detectors

Chao Zhang (张超),1, ∗ Yungui Gong (龚云贵),1, † Hang Liu (刘航),2, ‡ Bin Wang (王斌),3, 2, § and Chunyu Zhang (张春雨)1, ¶ 1School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei 430074, China 2School of Aeronautics and Astronautics, Shanghai Jiao Tong University, Shanghai 200240, China 3Center for Gravitation and Cosmology, Yangzhou University, Yangzhou 225009, China. Localizing the sky position of the gravitational wave source is a key scientific goal for gravitational wave observations. Employing the Fisher information matrix approximation, we compute the angular resolutions of LISA and TianQin, two planned space-based gravitational wave detectors and examine how detectors’ configuration properties, such as the orientation change of the detector plane, heliocentric or geocentric motion and the arm length etc., affect the accuracy of source localization for monochromatic sources at selected values of frequencies. We find that the amplitude modulation due to the annual changing orientation of the detector plane helps LISA get better accuracy in the sky localization and better sky coverage at frequencies below several mHz, and its effect on TianQin is negligible although the orientation of TianQin’s detector plane is fixed. At frequencies above roughly 30mHz, TianQin’s ability in the sky localization is better than LISA. Further we explore potential space detector networks for fast and accurate localization of the gravitational wave sources. The LISA-TianQin network has better ability in sky localization for sources with frequencies in the range 1-100 mHz and the network has larger sky coverage for the angular resolution than the individual detector.

I. INTRODUCTION all directions within a few seconds to minutes. It requires three or more ground-based GW detectors at widely sep- Since the first gravitational wave (GW) event arated sites to locate GW sources with the method of GW150914 observed by the Laser Interferometer timing triangulation approximation [23–25]. However, Gravitational-Wave Observatory (LIGO) Scientific Col- space-based GW detector can measure GWs for months laboration and the Virgo Collaboration [1,2], there have to years, the periodic Doppler shift due to the detector been tens of confirmed GW detections by ground-based motion in space results in amplitude and phase modula- GW observatories [1–13]. The ground-based GW ob- tions of the detected GW signals which encode informa- servatories can measure GWs within frequency range tion about the detector position and the angular position 10 103 Hz only, and are difficult to explore lower fre- of the source. Therefore, a single space-based GW detec- quency− band where a wealth of astrophysical signals re- tor is able to locate the source position. side. The proposed space-based observatories including Like LISA, Taiji is composed of a triangle of three LISA [14, 15], TianQin [16] and Taiji [17] are expected spacecrafts with a larger separation distance in a helio- to detect GWs in the low-frequency regime. centric orbit ahead instead of behind the Earth by about Accurately localizing GW sources is important, since 20◦ [17]. The similarities in configurations construe that the source position is correlated with the physical prop- similar to LISA, Taiji’s angular resolution depends simi- erties of the binary star system which are necessary larly on the type of signal and on how much other infor- to understand the formation and evolution of the bi- mation must be extracted. Similarities and complements nary. Moreover, the accurate GW source localization between LISA and Taiji imply that LISA-Taiji network may provide important information about the environ- can effectively help to accurately localize GW sources, ments where such relativistic objects reside. In particu- since the angular resolution measurements for the net- lar, the accurate knowledge of the GW source position work depend on the configuration angle and separation is essential for the follow-up observations of counterparts of the two constellations [26, 27]. The network was es- and the statistical identification of the host galaxy if no timated to improve the angular resolution over 10 times arXiv:2009.03476v3 [astro-ph.IM] 4 May 2021 counterpart is found, so that we can use GWs as stan- than each individual LISA or Taiji detector [28]. dard sirens to explore the universe’s expansion history Unlike LISA and Taiji, TianQin has a geocentric or- and understand the problem of Hubble tension [18–22]. bit configuration with three spacecrafts orbiting the One ground-based GW observatory cannot localize Earth and further rotating around the Sun together GW source, since it is sensible to GW signals from nearly with the Earth. The normal vector of TianQin’s detector plane points to the source RX J0806.3+1527 at ◦ ◦ (θtq = 4.7 , φtq = 120.5 ). TianQin is slightly more − ∗ chao [email protected] sensitive to GWs with higher frequency than LISA. The † Corresponding author. [email protected] precision of the parameter estimation and the sky lo- ‡ [email protected] calization of equal mass supermassive black hole binary § 5 7 wang [email protected] systems with masses in the range 10 10 M for Tian- ¶ [email protected] Qin was discussed in [29]. It was further− argued that the 2

LISA-TianQin network can improve the sky localization The organization of the paper is as follows. In Sec. II, of Galactic double white dwarf binaries up to 3 orders of we review the FIM method of signal analysis. In Sec. III, magnitude [30], if compared with single TianQin obser- we devise several fiducial detectors to discuss the effects vation. of different factors on the sky localization. The orbits For TianQin, considering its special configuration and for these detectors are presented in the Appendix. Then design, it is important to fully study its angular resolu- in Sec. IV we use these results to analyze the angular tion and uncover its dependencies. It is needed to be resolutions of TianQin, compare the result with LISA clarified in TianQin project how information about the and finally explore the combined LISA-TianQin network. source position encodes in Doppler shift, the translational We present our conclusions and discussions in the last motion of the detector relative to the source, detector’s section. changing orientation, the rotation period of spacecrafts II. SIGNAL ANALYSIS and the arm length of the detector. Different from LISA, whose amplitude modulation improves the sky localiza- For a GW signal tion accuracy below 1 mHz [31, 32], the detector plane of TianQin points to a fixed direction so there is no am- X A hij(t) = eijhA(t), (1) plitude modulation in TianQin which is more sensitive A to signals above several mHz. It is interesting to explore how much the sky localization accuracy of TianQin the output in the detector α is is affected in the absence of the amplitude modulation. X A On the other hand, the Doppler modulation becomes Hα(t) = Fα hA(t) +n ˆα(t), (2) stronger as the frequency of GWs increases and it has the A equatorial pattern [31–33], it is natural to ask how the where A = +, stands for the plus and cross polariza- Doppler modulation affects the sky localization accuracy A × of TianQin and how to compare it with other space-based tions, eij is the polarization tensor,n ˆα(t) is the detector A detectors at frequencies above 10 mHz. Moreover, there noise, the angular response function Fα for the polariza- are other questions to be answered, for example: How tion A is the equatorial pattern affects the sky localization accu- A X ij A racy? Is it possible that TianQin has better sky localiza- Fα = Dα eij, (3) tion accuracy than LISA and Taiji even if it has shorter i,j arm length? Neglecting differences in orbit configura- Dij tions, if we explore the potential TianQin-LISA network, and α is the detector tensor. In GR, there are two po- can it perform better in source localization than Taiji- larizations. For alternative theory of gravity, there may LISA network? We will carefully discuss these issues and exist up to six polarizations [34–42]. In this paper, we present sky maps of angular resolutions for monochro- consider GR only. For equal arm space-based interfero- matic sources to identify the sensitive regions of differ- metric detector with a single round trip light travel, the ent detectors. These discussions are helpful to improve detector tensor for a monochromatic GW with the fre- f ω the detector design. Throughout the work, we will em- quency propagating along the direction ˆ is ploy the Fisher information matrix approximation (FIM) ij 1 i j i j to give robust estimation of the sky localization of the Dα = [ûαuˆαT (f, uˆα ωˆ) vˆαvˆαT (f, vˆα ωˆ)], (4) 2 · − · source with a high signal to noise ratio (SNR) where the inverse of FIM gives the covariance matrix of the param- whereu ˆα andv ˆα are the unit vectors along the arms of eters. the detector α, T (f, uˆα ωˆ) is [43, 44] ·

1 f f f f T (f, uˆα wˆ) = sinc[ (1 uˆα ωˆ)] exp[ i (3 +u ˆα ωˆ)] + sinc[ (1 +u ˆα ωˆ)] exp[ i (1 +u ˆα ωˆ)] , (5) · 2{ 2f ∗ − · − 2f ∗ · 2f ∗ · − 2f ∗ · }

∗ sinc(x) = sin x/x, f = c/(2πL) is the transfer frequency where the power spectral density Sn,α(f) satisfies ∗ 1 0 of the detector, c is the speed of light and L is the arm n˜α(f)ñ (f) = δ(f f )Sn,α. For space-based interfer- h α i 2 − length of the detector. With the signal Hα, we define ometers, the noise power spectral density Sn(f) is [16, 45] SNR as S S −4 2 X x 4 a 10 Hz ρ = (H H ) Sn(f) = + (1 + ). (7) α α L2 (2πf)4L2 f α | ∞ (6) X Z 1 = 4 df H (f)H∗(f), For TianQin, the acceleration noise is √Sa = S f α α −15 −2 1/2 α 0 n,α( ) 10 m s /Hz , the displacement noise is √Sx = 3

using the Parseval’s theorem [31, 46], the FIM becomes 24 TianQin 10− LISA X ∂Hα ∂Hα 27 Γij = 10− ∂θ ∂θ α i j ) 1 10 30 ∞ − − Z X 4 ∗ = ∂iHα(f)∂jHα(f)df (10) (Hz 33 S f 10− n,α( 0) 0 n α S Z ∞ 10 36 X 2 ∗ − = ∂iHα(t)∂jHα(t)dt , Sn,α(f0) −∞ 39 α 10− where θi is the ith parameter and ∂iHα = ∂Hα/∂θi. The 5 3 1 10− 10− 10− covariance matrix of the parameters is Frequency (Hz) i j −1 σij = ∆θ ∆θ (Γ )ij. (11) ≈ FIG. 1. The noise curves for LISA and TianQin. For a detected source with a signiﬁcant SNR (a threshold of ρ 7), the angular uncertainty of the sky localization is evaluated≥ as 1/2 5 1 pm/Hz and the arm length is Lt = 1.7 10 km q ∗ × 2 [16]. Its transfer frequency is f = 0.28 Hz. ∆Ωs 2π cos θs σθsθs σφsφs σ . (12) t ≡ | | − θsφs For LISA, the acceleration noise is √Sa = 3 1/2 × The signal in the detector coordinate system is H(t) = 10−15 m s−2/Hz , the displacement noise is √S = x F +h (t) + F ×h (t) +n ˆ (t). In practice, for space-based 1/2 6 + × α 15 pm/Hz and the arm length is Ls = 2.5 10 km detectors, we often work in the heliocentric coordinate ∗ × [15]. Its transfer frequency is fs = 0.02 Hz. The noise system. The translational motion of the center of the de- curves for LISA and TianQin are shown in Fig.1. The tector around the Sun leads to an extra phase modulation ratios of the noise Sn(f0) between LISA and TianQin factor are 23 at f0 = 1 mHz, 1.2 at f0 = 10 mHz and 0.9 at iφD (t) 2πif0R cos θs cos(2πt/T −φs−φα)/c f0 = 100 mHz. e = e , (13)

Compact binaries containing stellar or intermediate where φα is the ecliptic longitude of the detector α at mass black holes, white dwarfs or neutron stars could t = 0, the period T of the rotation is 1 year and the emit monochromatic GWs in the early inspiral phase radius R of the orbit is 1 AU. Therefore, the signal in the which are detectable by space-based detectors and the heliocentric coordinate system is corresponding frequency evolution is negligible dur- + × iφD (t) ing the mission of the detector. In the lowest or- H(t) = [F (t)h+(t) + F (t)h×(t)]e . (14) der quadrupole approximation, the two polarizations of Since the phase modulation is very large in most cases monochromatic GWs with the frequency f are 0 and contains the information about the source location, this effect is very important for source localization and it 2 is the same for LISA, Taiji and TianQin. From Eq. (13), h+ = 1 + L~ ωˆ exp(2πif0t + iφ0), A · (8) we get ~ h× =2i L ωˆ exp(2πif0t + iφ0), 1 A · δf0 = ∂tφD(t) 2π where = 2M1M2/(rdL) is the amplitude, M1 and M2 A 2πR (15) are masses of the binary components, r is the distance = f0 cos θs sin (2πt/T φs φα) | − − | cT between them, dL is the luminosity distance between the v ~ f0 , source and the observer, L is the unit vector for the bi- ∼ c nary’s orbital angular momentum, and φ is the initial 0 so the periodic phase modulation spreads the measured phase. We focus mainly on the following seven param- power of the monochromatic signal of frequency f over a eters of the monochromatic GW signal included in Eq. 0 range f (1 v/c) and the spread depends on the direction (8) 0 of the source.± In addition to the Doppler modulation imposed by the θ = θs, φs, θL, φL, , φ0, f0 . (9) { A } rotation of the center of the detector around the Sun, the motion of the spacecrafts with respect to the detec- These parameters describe a monochromatic source by tor center, the orientation change of the detector plane, the source direction (θs, φs), the direction of the binary’s the arm length and the noises of the detector also affect orbital angular momentum (θL, φL), the amplitude , the accuracy of sky localization. In the next section, we A the initial phase φ0 and frequency f0. For monochro- analyze how much influences these effects have on the matic sources there is almost no frequency evolution, by angular resolution of the detector in detail. 4

III. THE EFFECT OF DIFFERENT CONSTELLATION ON SKY LOCALIZATION

LISA mission is proposed as an equilateral triangle SC3 constellation with sides of 2.5 106 km to detect low frequency GWs. The constellation× has an inclination angle of 60◦ with respect to the ecliptic plane and trails the Earth by about 20◦. The inclination angle ensures the spacecrafts to keep the geometry of an equilateral tri- L L angle throughout the mission. As shown in Fig.2, the 2 3 normal vector of the detector plane rotates around the normal vector of the ecliptic plane and forms a cone with 60◦ half opening angle in one year. Taiji has similar configuration to LISA except that its arm length is 3 106 ° 60 L km. TianQin is an equilateral triangle constellation× with 1 SC2 sides of 1.73 105 km designed to orbit the Earth with the SC1 period of 3.×65 days. The normal vector of the detector plane points to the direction of RX J0806.3+1527 with ◦ ◦ the latitude θtq = 4.7 and the longitude φtq = 120.5 FIG. 3. The orientation of three arms in the detector plane. as shown in Fig.2.− The centers of LISA, Taiji, and Tian- Qin all follow the Earth-like orbit. Three spacecrafts of the detectors are connected with each other by lasers in a form of equilateral triangle as shown in Fig.3, we label three spacecrafts as SC1/2/3, three arms as 1/2/3 and their lengths as L1/L2/L3. We refer to the Michelson interferometer formed by the arms 1 and 2 as detector I and the Michelson interferometer formed by the arms words, we assume that Sn(f0) = 1 for all fiducial detec- 2 and 3 as detector II. We use both detectors I and II tors so that the effect of each factor of the constellation α to determine the sky localization, so = I and II in Eq. on the source localization can be easily understood. Then (10). we just multiply the ratio of Sn(f0) to the result when we compare the sky localization of LISA and TianQin. Before we start to compute the localization accuracy of detectors, we must give the parameters of sources and ensure that they could be detected. Thus we fix the SNR of all sources for the detector R (will be defined below) to be 7, and then we derive the amplitudes of sources from this SNR. This means that the amplitudesA of sources being considered are different for different detectors. For a complimentary comparison, we present the results on the effect of different constellation by taking the same source (i.e., the sources have the same distance and masses) for different constellation in AppendixB. For this setting, perhaps some sources can not be detected by other detectors due to small SNR, we can increase the amplitudes of all sources by the same factor so that all sources can be detected by all detectors. Note that FIG. 2. The normal vectors of the detector planes for TianQin the relative angular resolution between detectors remains and LISA. Illustration from [47]. unchanged although the values of angular resolutions become smaller for each detector. The other parameters As analyzed above, there are four different factors are fixed as θL = 1.0 and φL = φ0 = 0. The mission which influence the detectors’ localization precision. To time or the observation time To is set to be one year. We disclose the effect of each factor on the source localiza- simulate 3600 sources uniformly distributed in the sky tion and compare its role in different detectors, we as- with π/2 < θs < π/2 and π < φs < π. We ana- sume the other three factors the same in different detec- lyze the− source localization accuracy− at three frequencies tors. Because we consider monochromatic sources, the 10−3 Hz, 10−2 Hz and 10−1 Hz. The three frequencies noise power at the source frequency f0 appears linearly stand for low, medium and high frequencies relative to ∗ in the Fisher matrix (10), so Sn(f0) can be eliminated the transfer frequency fs = c/(2πLs) = 0.02 Hz in Eq. by normalizing the result to a fixed SNR [31]. In other (5). 5

1.0 f0=0.001 Hz R 0.8 C C1 R1 0.6 C2

0.4 Probability 0.2

0.0 2 1 0 10− 10− 10 ∆ΩS

1.00 f0=0.01 Hz R C 0.75 C1 R1 C2 0.50

Probability 0.25 FIG. 4. The sky map of angular resolutions ∆ΩS of sources from diﬀerent directions in the unit of steradian (1 steradian 0.00 4 3 2 is 3000 square degrees) for the ﬁducial detectors R and C. 10− 10− 10− The horizontal axis represents the longitude φs and the ver- ∆ΩS tical axis represents the latitude θs. The left panel is for the 1.00 f0=0.1 Hz detector R and the right panel is for the detector C. From R C top to bottom, the frequencies of monochromatic sources are 0.75 C1 −3 −2 −1 R1 10 Hz, 10 Hz and 10 Hz. C2 0.50

Probability 0.25 A. The rotation eﬀect

0.00 7 6 5 4 10− 10− 10− 10− One major difference between TianQin and LISA/Taiji ∆ΩS is the orientation of the detector plane. Naively we expect that the modulation of the signal caused by the rota- FIG. 5. Cumulative histograms of sky localization estimations tion of the detector plane increases the accuracy of source ∆ΩS for different detectors. localization for LISA and Taiji, and the precision of sky localization for TianQin becomes less without this modulation effect. To evaluate the effect of this time-changing f0(Hz) R C C1 C2 R1 orientation of the detector plane on sources’ angular reso- Mean Value 10−1 2.4 10−5 2.7 10−5 1.5 10−6 1.2 10−6 1.1 10−6 lutions, we construct two fiducial GW detectors with the × × × × × 10−2 1.8 10−3 2.3 10−3 2.0 10−3 2.1 10−3 1.6 10−3 same arm length, noise curve and rotation period of the × × × × × spacecrafts. The first detector is like LISA except that 10−3 5.9 10−2 1.4 10−1 1.5 10−1 1.7 10−1 5.9 10−2 × × × × × its arm length is 3.7 109 m, and we call it detector R. Median Value × −1 −6 −5 −7 −7 −7 The second detector is like TianQin except that its arm 10 7.0 10 1.0 10 4.7 10 4.8 10 4.2 10 9 × × × × × length is 3.7 10 m and its rotation period around the 10−2 7.6 10−4 1.0 10−3 7.9 10−4 7.9 10−4 5.0 10−4 × × × × × × Earth is one year, and we call it detector C. The detailed 10−3 4.8 10−2 7.1 10−2 7.1 10−2 7.6 10−2 4.8 10−2 orbit equations for the detectors R and C are presented × × × × × in AppendixA. For the arm length L = 3.7 109 m, the TABLE I. The mean and median values of angular resolutions ∗ × transfer frequency is f = c/(2πL) = 0.013 Hz. ∆ΩS for different detectors. The results of angular resolutions at the frequencies −1 −2 −3 f0 = 10 Hz, 10 Hz and 10 Hz for the detectors R and C are shown in Figs.4,5 and6. We show the sky olution at the frequency f0 = 1 mHz with the amplitude map of the angular resolution for the detectors R and modulation caused by the detector’s time-varying orien- C in Fig.4, the cumulative histograms of the angular tation. However, the median and mean values of angular resolution in Fig.5 and the ratios of the mean and me- resolutions of the detector as usually presented in the lit- dian angular resolutions between the detectors C and R erature are incomplete because we don’t know where the in Fig.6. The mean and median values of angular reso- poorly localized regions are. Thus we also present the lutions are summarized in TableI. From Figs.5 and6, sky map of angular resolutions in Fig.4. we see that the effect of the detector’s time-varying ori- From Fig.4, we see that there are some regions in entation decreases as the frequency of GWs increases, the which the uncertainties of the angular resolutions are improvement is only a factor of 2.4 for the mean angular large, but the locations of these regions cannot be seen resolution and a factor of 1.5 for the median angular res- in the histogram. As the frequencies of the monochro- 6 matic GW sources increase, the uncertainties of the an- Mean C/R gular resolutions become smaller. For the detector R, at 2.25 −3 Median C/R the low frequency f0 = 10 Hz, sky localization uncer- 2.00 tainties ∆ΩS are around 0.01 to 0.1 steradians across the 1.75 sky. Comparing with sources from other directions, ∆ΩS are more than two orders of magnitude larger for sources 1.50 1.25 along the equatorial plane with the latitude θs = 0 in the heliocentric coordinate system at the medium and high 3 2 1 10− 10− 10− frequencies. For the detector C, the angular resolutions Frequency (Hz) ◦ for sources along the detector plane (φs is around 150 and 30◦) and the equatorial plane in the heliocentric− co- FIG. 6. The ratios of the mean and median values of the ordinate system are much worse. angular resolutions between the detectors R and C. In Fourier space, the Doppler effect spreads the measured power of the monochromatic signal of frequency f0 over a range f0(1 v/c), where the speed of the motion B. The effect of the arm length of the detector’s center± around the Sun v 2πR/T = 4 ∼ 3 10 m/s, so the sky localization is better when f0 increases× and the frequency spread caused by the Doppler Another difference between TianQin and LISA/Taiji is −3 −7 the distance between spacecrafts. When the wavelength effect at f0 = 10 Hz is δf = f0v/c 10 . The Doppler effect on the angular resolution∼ for monochro- of GWs is larger than the arm length, i.e., f0 < f∗, the matic sources is [32] transfer function (5) is almost equal to 1, so the transfer function is independent of GWs’ frequency and the effect of the arm length on the transfer function is neg- f > f 2 ligible. If 0 ∗, then the frequency dependence of the 1 6Rs 2 2 δΩ = csc (θs) cos(θs) transfer function deteriorates the response of the detec- 2 2 2 2 2 f πR ρ √1 6x To tor. To evaluate the influence of the detector’s arm length s − 1/2 on angular resolutions, we devise two more fiducial GW Rs (16) 12 x cot(θ ) sin(φ ) csc(θ ) + csc2(θ ) detectors with the same noise curve, rotation period of T s s s s − 0 the spacecrafts and orientation of the detector plane but 2 2 −1 0.02(1mHz/f0) (10/ρ) sin(θs) , with different arm lengths. The third detector is like the ≈ | | detector C except that its arm length is 1.7 108 m, and we call it detector C1. For the detector C1,× we assume a where the distance between the Earth and the Sun Rs = fiducial Earth with the right mass to provide the orbit. 1 AU, the observational time To = 1 year and x = 1/π. Since the mass of the Earth does not affect the perfor- It is easy to see that the sky localization becomes better mance of the detector, so there is no problem with this as frequency increases and it depends on the direction assumption for the purpose of the discussion on the ef- of the source. For sources along the equatorial plane, fect of the arm length. The fourth detector is like the . 8 the latitude θs = 0, the sky localization becomes worse, detector R except that its arm length is 1 7 10 m, and × 8 so the Doppler modulation has the equatorial pattern we call it detector R1. For the arm length Lt = 1.7 10 ∗ × [32]. On the other hand, the time-varying orientation m, the transfer frequency is ft = c/2πLt = 0.28 Hz. The of the detector’s plane spreads the power over a range detailed orbit equations for the detectors C1 and R1 are −8 f0 2/T = f0 6.3 10 with T = 1 year (i.e., the described in AppendixA. The effect of the arm length frequency± spread± δf =× 2/T = 6.3 10−8), and it is in- can be analyzed by comparing the angular resolutions for × −3 dependent of the source location. At f0 = 10 Hz, the either the detectors R and R1 or C and C1. frequency spread caused by the Doppler modulation is Figure7 shows the sky map of angular resolutions for δf 10−7, so the frequency spread caused by the ampli- the detectors R1 and C1. Because of the Doppler effect, tude∼ modulation is the same order as that caused by the the angular resolution is better at higher frequency. Since Doppler modulation, this explains why the angular reso- the frequencies f0 = 1 mHz and f0 = 10 mHz are smaller lution for the detector R is almost the same for sources than the transfer frequencies f ∗ = 13 mHz (for the arm 9 ∗ from any direction, especially good angular resolutions length L = 3.7 10 m) and ft = 0.28 Hz (for the arm × 8 along the equatorial plane. As the frequency of GWs length Lt = 1.7 10 m), so the angular resolutions are increases, the Doppler effect dominates. Therefore the similar for either× the detectors R and R1 or the detectors angular resolutions for the detector R and C are almost C and C1 at f0 = 1 mHz and f0 = 10 mHz as shown the same at medium and high frequencies, and the sky in Figs.5 and8. At the frequency f0 = 0.1 Hz, the localization becomes worse for sources along the equa- wavelength of GWs is comparable with the arm length 8 torial plane. For the detector C, if sources are located Lt = 1.7 10 m and less than the arm length L = along the detector’s plane, then the SNR is smaller and 3.7 109 m,× so angular resolutions of the detectors R1 the sky localization becomes worse. and× C1 are better than those of the detectors R and C, as 7

FIG. 7. The sky map of angular resolutions ∆ΩS for the detectors R1, C1 and C2. The left panel is for the detector R1, the middle panel is for the detector C1 and the right panel is for the detector C2. From top to bottom, the frequencies are 10−3 Hz, 10−2 Hz and 10−1 Hz.

of angular resolutions for the detector C2, respectively. 20 Mean value C/C1 The mean and median values of angular resolutions for Median value C/C1 15 Mean value R/R1 the detector C2 are given in Table.I. These results show Median value R/R1 that the inﬂuence of the rotation period of the spacecrafts 10 is almost negligible. Because of the cylindrical symmetry

5 of the rotation, the rotations of the spacecrafts around the ﬁxed axis have no eﬀect on angular resolutions in one 0 3 2 1 period. 10− 10− 10− Frequency (Hz) In the real geocentric orbit, the rotation period depends on the arm length. If we increase the arm length, the detector’s sensitivity increases and the detector is FIG. 8. The ratios of the mean and median values of angular more sensitive at lower frequency. resolutions between the detectors R and R1, and the detectors C and C1.

IV. LISA-TIANQIN NETWORK shown in Figs.5 and8. We also summarize the mean and median values of angular resolutions for these detectors At the frequency f0 = 1 mHz, the contribution of the in TableI. These results tell us that the effect of the amplitude modulation due to the time changing orien- arm length is almost the same for the heliocentric and tation of the detector plane with a period of one year geocentric constellations and the factor of improvement to the accuracy of the angular resolution is comparable on the angular resolution is about L/Lt = 22 at f0 = 0.1 to that of the Doppler modulation due to the motion of Hz. As expected, TianQin is more sensitive than LISA the center of the detector around the Sun, so the am- at frequencies f0 & 10 mHz. plitude modulation helps LISA and Taiji not only get better angular resolution, but also enlarge the sky coverage because the amplitude modulation is independent C. Rotation period of the spacecrafts of the sources’ directions. At higher frequencies when the wavelength of GWs is larger than the detector’s arm The rotation period of the spacecrafts for TianQin is length, the frequency dependent transfer function dete- 3.65 days and it is 1 year for LISA/Taiji. To evaluate riorates the SNR registered in the detector, so we expect the effect of the rotation period of the spacecrafts on that the accuracy of the sky localization for TianQin is sources’ angular resolutions, we devise another fiducial better because of its shorter arm length. Now we discuss GW detector which is like the detector C1 except that the sky localization estimations of LISA, TianQin and its rotation period is 3.65 days, and we call it detector the LISA-TianQin network. Since Taiji and LISA have C2. The detailed orbit equations for the detector C2 are similar constellation, so we discuss LISA only. We fix the given in AppendixA. Without specifying the noises, the SNR of all sources for LISA to be 7, then we derive the detector C2 is the same as TianQin. amplitudes of sources from this fixed SNR, and cal- Figures5 and7 show the histograms and the sky map culate the angularA resolutions of LISA, TianQin and the 8

LISA-TianQin network. 1.00 f0=0.001 Hz We show the histograms of angular resolutions in Fig. LISA 0.75 TianQin 9 and summarize the mean and median values of angu- Net −3 lar resolutions in TableII. At the frequency f0 = 10 Hz, LISA’s angular resolution is roughly 50 times bet- 0.50

S f Probability ter than TianQin because its noise n( 0) is about 23 0.25 times smaller and the amplitude modulation contributes another factor of 2.4, the angular resolution of the com- 0.00 3 2 1 0 1 2 bined network is almost the same as LISA. At the fre- 10− 10− 10− 10 10 10 −2 ∆ΩS quency f0 = 10 Hz, LISA’s angular resolution is about 1.00 1.4 times better than TianQin because TianQin’s noise f0=0.01 Hz LISA Sn(f0) is about 1.2 times smaller and its rotation ef- 0.75 TianQin fect is about 1.2 times smaller, and the network’s an- Net gular resolution is a little better than LISA. At the fre- 0.50 quency f0 = 0.1 Hz, the frequency dependent transfer Probability function reduces LISA’s sensitivity by a factor about 10, 0.25 so TianQin’s angular resolution becomes about 10 times 0.00 4 3 2 better than LISA. Therefore, the angular resolution of 10− 10− 10− the combined LISA-TianQin network spans over the fre- ∆ΩS −6 quency ranges 1-100 mHz and it reaches 10 steradians 1.00 f0=0.1 Hz at f0 = 0.1 Hz. LISA 0.75 TianQin In Fig. 10, we plot the sky map of angular resolu- Net tions for LISA, TianQin and their combined network. 0.50 At f0 = 1 mHz, the amplitude and Doppler modulations contribute to LISA’s angular resolutions, so LISA’s Probability 0.25 angular resolutions are almost the same across the sky.

At higher frequencies, the contribution of the amplitude 0.00 7 6 5 4 modulation is negligible and only the Doppler modula- 10− 10− 10− 10− ∆ΩS tion matters. For sources along the equatorial plane, the angular resolution is the worst for both LISA and FIG. 9. Cumulative histograms of sky localization uncertain- TianQin at the frequencies f0 = 10 mHz and f0 = 100 mHz. For TianQin, the worst angular resolution also oc- ties ∆ΩS for LISA, TianQin and their combined network. ◦ curs for sources from the directions with φs = 30 and ◦ φs = 150 . Fig. 10 shows that the combined network − enlarges the sky coverage in addition to the slight im- comes better as the frequency of the monochromatic GW provement on angular resolutions at f0 = 10 mHz and source increases, and this effect depends on the direction f0 = 100 mHz. of the source. For sources along the plane with θs = π/2 in the heliocentric coordinate system, the improvement by the Doppler modulation on the angular resolution is V. CONCLUSION the largest. For frequencies above approximately 5 mHz, the contribution of the amplitude modulation is negligi- The detector’s noises, arm length, time-changing orien- ble and only the Doppler modulation matters. Therefore, tation, motion around the Sun are four main factors that TianQin’s angular resolution is not affected much due to influence the accuracy of source localization for spaced- the lack of amplitude modulation. By pointing to the based GW observatories. The detector’s noises and arm specific source, TianQin has some blind spots, its angular resolution is not good for sources from the directions length affect the angular resolution through SNR. The ◦ ◦ with φs around 30 or 150 . Doppler modulation on the amplitude and phase of the − GW signal imposed by the translational motion of the de- At the frequency f0 = 1 mHz, the contribution of tector’s center around the Sun carries the position infor- the amplitude modulation is comparable to that of the mation of the source. For LISA and Taiji, the amplitude Doppler modulation, so the amplitude modulation helps of GW signal is further modulated by the annual rotation LISA and Taiji not only get better angular resolution, of the detector’s orientation. Both the Doppler and am- but also enlarge the sky coverage. At higher frequen- plitude modulations spread the power of the monochro- cies when the wavelength of GWs is larger than the de- matic GWs over a frequency range in Fourier space. The tector’s arm length, the frequency dependent transfer amplitude modulation improves the accuracy of the sky function deteriorates the SNR registered in the detec- localization below several mHz and the improvement is tor. For the monochromatic sources with the frequency ∗ ∗ the same for sources from different directions. The effect fs = c/(2πLs) < f0 < ft = c/(2πLt), TianQin’s an- of the Doppler modulation on the angular resolution be- gular resolution is better than LISA’s by a factor of 9

LISA TianQin Network

f0(Hz) Mean Median Mean Median Mean Median 10−1 2.3 10−5 7.4 10−6 2.1 10−6 8.3 10−7 1.8 10−6 7.1 10−7 × × × × × × 10−2 2.0 10−3 8.2 10−4 2.5 10−3 1.1 10−3 8.9 10−4 3.8 10−4 × × × × × × 10−3 6.2 10−2 5.03 10−2 3.7 1.7 5.7 10−2 4.9 10−2 × × × × TABLE II. The mean and median values of angular resolutions for monochromatic sources.

−3 FIG. 10. The sky map of angular resolutions ∆ΩS for LISA, TianQin and their combined network at the frequencies f0 = 10 −2 −1 Hz, f0 = 10 Hz and f0 = 10 Hz.

roughly Sn,tLs/(Sn,sLt), here Sn,s and Ls are LISA’s the blind spots can be avoided. noise curve and arm length, and Sn,t and Lt are Tian- Qin’s noise curve and arm length. Since LISA and Taiji have the best sensitivity at around 1 mHz, and TianQin ACKNOWLEDGMENTS is more sensitive at 10 mHz, the LISA-TianQin or TaiJi- TianQin network have better ability of sky localization This research was supported in part by the National for sources in the frequency range 1-100 mHz and the Key Research and Development Program of China under network has larger sky coverage for the angular resolu- Grant No. 2020YFC2201504, the National Natural Sci- tion than the individual detector. By assuming SNR=7 ence Foundation of China under Grants No. 11875136 for LISA, the network’s angular resolution is about 200 and No. 12075202, and the Major Program of the Na- square degrees at 1 mHz, 3 square degrees at 10 mHz and tional Natural Science Foundation of China under Grant 0.005 square degrees at 100 mHz. The understanding of No. 11690021. the eﬀects of the amplitude and Doppler modulations on Appendix A: DETECTOR’S ORBITS LISA/Taiji/TianQin can be used to optimize their constellations. In particular, the result may help TianQin 1. The orbits for the detectors C, C1, and C2 to improve the design so that the equatorial pattern and T In the ecliptic coordinate system, the orbit ~rn (t) = T T T (Xn (t),Yn (t),Zn (t)), (n = 1, 2, 3) is

T 0 0 X (t) =RT [cos θtq cos φtq sin(αT n α ) + sin φtq cos(αT n α )] n − T 0 − T 0 1 0 0 + eT RT [sin φtq(cos 2(αT n α ) 3) + cos θtq cos φtq sin 2(αT n α )] 2 − T 0 − − T 0 (A1) 1 2 0 2 0 0 + e RT [ 6 sin φtq cos(αT n α ) sin (αT n α ) + cos θtq cos φtq(3 cos 2(αT n α ) 1)] 4 T − − T 0 − T 0 − T 0 − 1 3 2 2 3 3 + R cos(αT αT 0) + eR(cos 2(αT αT 0) 3) e R cos(αT αT 0) sin (αT αT 0) + O(e , e ), − 2 − − − 2 − − 1 10

T 0 0 Y (t) =RT [cos θtq sin φtq sin(αT n α ) cos φtq cos(αT n α )] n − T 0 − − T 0 1 0 0 + eT RT [ cos φtq(cos 2(αT n α ) 3) + cos θtq sin φtq sin 2(αT n α )] 2 − − T 0 − − T 0 1 2 0 2 0 + e RT [6 cos φtq cos(αT n α ) sin (αT n α )+ (A2) 4 T − T 0 − T 0 0 0 cos θtq sin φtq sin(αT n α )(3 cos 2(αT n α ) 1)] − T 0 − T 0 − 1 1 2 3 3 + R sin(αT αT 0) + eR sin 2(αT αT 0) + e R sin(αT αT 0)(e cos 2(αT αT 0) 1) + O(e , e ), − 2 − 4 − − − 1

T 0 1 0 Zn (t) = RT sin θtq sin(αT n αT 0) eT RT sin θtq sin 2(αT n αT 0) − − − 2 − (A3) 1 2 0 0 e RT sin θtq sin(αT n α ) [3 cos 2(αT n α ) 1], − 4 T − T 0 × − T 0 −

◦ ◦ where (θtq = 4.7 , φtq = 120.5 ) is the location of the with the same distance and masses for diﬀerent de- − source RX J0806.3+1527, eT is the eccentricity of the tectors. To consider monochromatic GWs with the orbit which we set to be 0 in this paper, e = 0.0167 is frequency (10−1, 10−2, 10−3) Hz, we take the sources the eccentricity of the Earth’s orbit, R = 1 AU, αT 0 and to be equal mass binary systems with the total mass 0 2 4 3 4 α are initial phases, αT = 2πfmt + κ0, fm = 1/(365 (6, 3 10 , 10 ) M and the distance (2.3, 1.3 10 , 10 ) T 0 × × days), αT n = 2πfsct + 2π(n 1)/3. For the detector Mpc, respectively. The noise power Sn(f0) is chosen as − −41 −41 −41 C, the semimajor axis of the spacecraft’s orbit RT = (3.6 10 , 3.6 10 , 4.1 10 ) at the frequency 6 −×1 −2 −3× × 3.7 10 km, fsc = 1/(365 days). For the detector C1, (10 , 10 , 10 ), respectively. The sky map of angular × 5 RT = 1.7 10 km, fsc = 1/(365 days). For the detector resolutions for the detectors R, C, R1, C1 and C2 are × 5 C2, RT = 1.7 10 km, fsc = 1/(3.65 days). shown in Figs. 11 and 12 which are consistent with our × conclusions (SNR=7).

2. The orbits for the detectors R and R1

L L L L The orbit ~rn (t) = (Xn (t),Yn (t),Zn (t)) in the ecliptic coordinate system is

L Xn (t) = R(cos αLn + eL) cos cos θLn q 2 R 1 e sin αLn sin θLn, (A4) − − L L Yn (t) = R(cos αLn + eL) cos sin θLn q 2 +R 1 e sin αLn cos θLn, (A5) − L L Zn (t) = R(cos αLn + eL) sin , (A6)

2 4 2 1 where eL = (1 + √ σ + σ ) 2 1, = arctan[σ/(1 + 3 3 − σ/√3)], σ = L/(2R), R 1 AU, θLn = 2π(n 1)/3, ≈ − αLn + eL sin αLn = 2πfmt 2π(n 1)/3 αL0 and αL0 is the initial phase. The arm− length− L = 3−.7 109 m for the detector R and L = 1.7 108 m for the detector× R1. × FIG. 11. The sky map of angular resolutions ∆ΩS for the ﬁdu- cial detectors R and C. For the same frequency, the sources have the same amplitude. The left panel is for the detector R Appendix B: THE EFFECT DIFFERENT and the right panel is for the detector C. From top to bottom, CONSTELLATION FOR SOURCES WITH FIXED the frequencies of monochromatic sources are 10−3 Hz, 10−2 AMPLITUDE Hz and 10−1 Hz.

In this section, instead of ﬁxing the SNR, we ﬁx the amplitude of the sources, i.e., we consider sources

[1] B. P. Abbott et al. (LIGO Scientiﬁc and Virgo Collabora- Black Hole Merger, Phys. Rev. Lett. 116, 061102 (2016). tions), Observation of Gravitational Waves from a Binary 11

FIG. 12. The sky map of angular resolutions ∆ΩS for the detectors R1, C1 and C2. For the same frequency, the sources have the same amplitude The left panel is for the detector R1, the middle panel is for the detector C1 and the right panel is for the detector C2. From top to bottom, the frequencies are 10−3 Hz, 10−2 Hz and 10−1 Hz.

[2] B. P. Abbott et al. (LIGO Scientific and Virgo Collabo- [10] R. Abbott et al. (LIGO Scientific and Virgo Collabora- rations), GW150914: The Advanced LIGO Detectors in tions), GW190412: Observation of a binary-black-hole the Era of First Discoveries, Phys. Rev. Lett. 116, 131103 coalescence with asymmetric masses, Phys. Rev. D 102, (2016). 043015 (2020). [3] B. P. Abbott et al. (LIGO Scientific and Virgo Collabora- [11] R. Abbott et al. (LIGO Scientific and Virgo Collabora- tions), GW151226: Observation of Gravitational Waves tions), GW190814: Gravitational waves from the coales- from a 22-Solar-Mass Binary Black Hole Coalescence, cence of a 23 solar mass black hole with a 2.6 solar mass Phys. Rev. Lett. 116, 241103 (2016). compact object, Astrophys. J. Lett. 896, L44 (2020). [4] B. P. Abbott et al. (LIGO Scientific and Virgo Collabo- [12] R. Abbott et al. (LIGO Scientific and Virgo Collabora- rations), GW170104: Observation of a 50-Solar-Mass Bi- tions), GW190521: A Binary Black Hole Merger with nary Black Hole Coalescence at Redshift 0.2, Phys. Rev. a Total Mass of 150M , Phys. Rev. Lett. 125, 101102 Lett. 118, 221101 (2017);, Erratum: Phys.Rev.Lett. 121, (2020). 129901 (2018). [13] R. Abbott et al. (LIGO Scientific and Virgo Collabo- [5] B. P. Abbott et al. (LIGO Scientific and Virgo Collab- rations), GWTC-2: Compact Binary Coalescences Ob- orations), GW170814: A Three-Detector Observation of served by LIGO and Virgo During the First Half of the Gravitational Waves from a Binary Black Hole Coales- Third Observing Run, arXiv:2010.14527 [Phys. Rev. X cence, Phys. Rev. Lett. 119, 141101 (2017). (to be published)]. [6] B. P. Abbott et al. (LIGO Scientific and Virgo Collabora- [14] K. Danzmann, LISA: An ESA cornerstone mission for a tions), GW170817: Observation of Gravitational Waves gravitational wave observatory, Classical Quantum Grav- from a Binary Neutron Star Inspiral, Phys. Rev. Lett. ity 14, 1399 (1997). 119, 161101 (2017). [15] P. Amaro-Seoane et al. (LISA Collaboration), Laser In- [7] B. . P. . Abbott et al. (LIGO Scientific and Virgo Col- terferometer Space Antenna, arXiv:1702.00786 [astro- laborations), GW170608: Observation of a 19-solar-mass ph.IM]. Binary Black Hole Coalescence, Astrophys. J. 851, L35 [16] J. Luo et al. (TianQin Collaboration), TianQin: a space- (2017). borne gravitational wave detector, Classical Quantum [8] B. P. Abbott et al. (LIGO Scientific and Virgo Collabora- Gravity 33, 035010 (2016). tions), GWTC-1: A Gravitational-Wave Transient Cata- [17] W.-R. Hu and Y.-L. Wu, The Taiji program in space log of Compact Binary Mergers Observed by LIGO and for gravitational wave physics and the nature of gravity, Virgo during the First and Second Observing Runs, Phys. Natl. Sci. Rev. 4, 685 (2017). Rev. X 9, 031040 (2019). [18] B. F. Schutz, Determining the Hubble Constant from [9] B. P. Abbott et al. (LIGO Scientific and Virgo Collab- Gravitational Wave Observations, Nature (London) 323, orations), GW190425: Observation of a compact binary 310 (1986). coalescence with total mass 3.4M , Astrophys. J. Lett. [19] B. P. Abbott et al. (LIGO Scientific, Virgo, 1M2H, Dark 892, L3 (2020). ∼ Energy Camera GW-E, DES, DLT40, Las Cumbres Ob- 12

servatory, VINROUGE, MASTER Collaborations), A [33] C. Zhang, Y. Gong, B. Wang, and C. Zhang, Accuracy of gravitational-wave standard siren measurement of the estimation of parameters with space-borne gravitational Hubble constant, Nature (London) 551, 85 (2017). wave observatory, arXiv:2012.01043 [gr-qc]. [20] H.-Y. Chen, M. Fishbach, and D. E. Holz, A two per [34] D. M. Eardley, D. L. Lee, and A. P. Lightman, cent Hubble constant measurement from standard sirens Gravitational-wave observations as a tool for testing rel- within ﬁve years, Nature 562, 545 (2018). ativistic gravity, Phys. Rev. D 8, 3308 (1973). [21] W. J. Wolf and M. Lagos, Standard Sirens as a Novel [35] D. Liang, Y. Gong, S. Hou, and Y. Liu, Polarizations of Probe of Dark Energy, Phys. Rev. Lett. 124, 061101 gravitational waves in f(R) gravity, Phys. Rev. D 95, (2020). 104034 (2017). [22] A. G. Riess, S. Casertano, W. Yuan, L. M. Macri, and [36] S. Hou, Y. Gong, and Y. Liu, Polarizations of Gravita- D. Scolnic, Large Magellanic Cloud Cepheid Standards tional Waves in Horndeski Theory, Eur. Phys. J. C 78, Provide a 1% Foundation for the Determination of the 378 (2018). Hubble Constant and Stronger Evidence for Physics be- [37] Y. Gong and S. Hou, Gravitational wave polarizations in yond ΛCDM, Astrophys. J. 876, 85 (2019). f(R) gravity and scalar-tensor theory, EPJ Web Conf. [23] S. Fairhurst, Triangulation of gravitational wave sources 168, 01003 (2018). with a network of detectors, New J. Phys. 11, 123006 [38] Y. Gong, E. Papantonopoulos, and Z. Yi, Constraints (2009); Erratum: New J.Phys. 13, 069602 (2011). on scalar–tensor theory of gravity by the recent observa- [24] S. Fairhurst, Source localization with an advanced gravi- tional results on gravitational waves, Eur. Phys. J. C 78, tational wave detector network, Classical Quantum Grav- 738 (2018). ity 28, 105021 (2011). [39] Y. Gong, S. Hou, D. Liang, and E. Papantonopoulos, [25] K. Grover, S. Fairhurst, B. F. Farr, I. Mandel, C. Ro- Gravitational waves in Einstein-æther and generalized driguez, T. Sidery, and A. Vecchio, Comparison of Grav- TeVeS theory after GW170817, Phys. Rev. D 97, 084040 itational Wave Detector Network Sky Localization Ap- (2018). proximations, Phys. Rev. D 89, 042004 (2014). [40] Y. Gong and S. Hou, The polarizations of gravitational [26] W.-H. Ruan, C. Liu, Z.-K. Guo, Y.-L. Wu, and R.-G. Cai, waves, Universe 4, 85 (2018). The LISA-Taiji network: precision localization of massive [41] Y. Gong, S. Hou, E. Papantonopoulos, and D. Tzortzis, black hole binaries, Research 2021, 6014164 (2021). Gravitational waves and the polarizations in Hoˇrava [27] W.-H. Ruan, C. Liu, Z.-K. Guo, Y.-L. Wu, and R.-G. gravity after GW170817, Phys. Rev. D 98, 104017 (2018). Cai, The LISA-Taiji network, Nat.Astron. 4, 108 (2020). [42] S. Hou and Y. Gong, Gravitational waves in Einstein- [28] G. Wang, W.-T. Ni, W.-B. Han, S.-C. Yang, and X.- Æther theory and generalized TeVeS theory after Y. Zhong, Numerical simulation of sky localization for GW170817, Universe 4, 84 (2018). LISA-TAIJI joint observation, Phys. Rev. D 102, 024089 [43] F. B. Estabrook and H. D. Wahlquist, Response of (2020). Doppler spacecraft tracking to gravitational radiation, [29] W.-F. Feng, H.-T. Wang, X.-C. Hu, Y.-M. Hu, and Gen. Relat. Gravit. 6, 439 (1975). Y. Wang, Preliminary study on parameter estimation [44] N. J. Cornish and S. L. Larson, Space missions to de- accuracy of supermassive black hole binary inspirals for tect the cosmic gravitational wave background, Classical TianQin, Phys. Rev. D 99, 123002 (2019). Quantum Gravity 18, 3473 (2001). [30] S.-J. Huang, Y.-M. Hu, V. Korol, P.-C. Li, Z.-C. Liang, [45] X.-C. Hu, X.-H. Li, Y. Wang, W.-F. Feng, M.-Y. Zhou, Y. Lu, H.-T. Wang, S. Yu, and J. Mei, Science with the Y.-M. Hu, S.-C. Hu, J.-W. Mei, and C.-G. Shao, Funda- TianQin Observatory: Preliminary results on Galactic mentals of the orbit and response for TianQin, Classical double white dwarf binaries, Phys. Rev. D 102, 063021 Quantum Gravity 35, 095008 (2018). (2020). [46] A. Vecchio and E. D. L. Wickham, The Eﬀect of the [31] C. Cutler, Angular resolution of the LISA gravitational LISA response function on observations of monochro- wave detector, Phys. Rev. D 57, 7089 (1998). matic sources, Phys. Rev. D 70, 082002 (2004). [32] A. Blaut, Accuracy of estimation of parameters with [47] D. Liang, Y. Gong, A. J. Weinstein, C. Zhang, and LISA, Phys. Rev. D 83, 083006 (2011). C. Zhang, Frequency response of space-based interfer- ometric gravitational-wave detectors, Phys. Rev. D 99, 104027 (2019).