Arxiv:2105.00746V2 [Gr-Qc]

arXiv:2105.00746v2 [gr-qc] 11 May 2021 ntuetnieadatohsclfrgon,jitob- joint foreground, the astrophysical from early and imprint the relic noise cosmological on instrument stochastic distinguish understanding the To our of impact Universe. deeply detection will The GW background GW [8–13]. stochastic the (SGWB) bi- for Virgo searching compact Advanced actively Although the and are from LIGO 7]. Advanced all coalescences, are [6, nary detections parame- sources GW GW and current the localization the on sky determination the ter The of the improve resolution 6]. will angular Earth [5, the future surrounding ground- near network the detector in join network to interferometer expected is based detector [2– KAGRA tested The were polarizations 4]. GW local- alternative well the were and LIGO directions ized source between the orientation detectors, misaligned Virgo and the of Virgo. benefits triple Advanced As and by LIGO observed Advanced GW170817 coincidently of and interferometers detections GW170814 first The the were km. a with 3000 arms of interferometric de- separation aligned are (closely) interferometers be Two to signed [1]. LA Hanford, Livingston, sites and two WA at detectors LIGO Advanced by served ¶ ∗ § ‡ † ie u:[email protected] Luo: Ziren [email protected] Xu: Peng [email protected] Han: Wen-Biao [email protected] Ni: Wei-Tou [email protected] [email protected], Wang: Gang h rvttoa aedtcin W594 a ob- was GW150914, detection, wave gravitational The 4 5 aghuIsiuefrAvne td,Uiest fChine of University Study, Advanced for Institute Hangzhou colo srnm n pc cec,Uiest fChines of University Science, Space and Astronomy of School agWang, Gang o h ieetlcto n retto opstoso L of th compositions investigate orientation we and location work, sensitivit different this comparable the In with for detectors advantages. dual significant The with 2030s. the in ol civ h etcoscreainadosretest the observe and complement sensitivity. cross-correlation better best their the to achieve localiz due could sky system binary the massive on the performance for best the shows network TAIJIm nlnto n edn h at y20 by Earth the leading and inclination n edn h at y20 by Earth the leading and h at y20 by Earth the 1 h pc-on rvttoa ae(W eetr,LS an LISA detectors, (GW) wave gravitational space-borne The hnhiAtooia bevtr,CieeAaeyo Sci of Academy Chinese Observatory, Astronomical Shanghai .INTRODUCTION I. 3 eateto hsc,Ntoa sn u nvriy Hsin University, Hua Tsing National Physics, of Department 7 azo etrfrTertclPyis azo Universit Lanzhou Physics, Theoretical for Center Lanzhou 2 noainAaeyfrPeiinMaueetSineadT and Science Measurement Precision for Academy Innovation tt e aoaoyo antcRsnneadAoi n Mo and Atomic and Resonance Magnetic of Laboratory Key State 6 nttt fMcais hns cdm fSine,Beiji Sciences, of Academy Chinese Mechanics, of Institute 1, 8 aj aoaoyfrGaiainlWv nvre(Beijing Universe Wave Gravitational for Laboratory Taiji ∗ nvriyo hns cdm fSine,Biig100049, Beijing Sciences, of Academy Chinese of University e-o Ni, Wei-Tou ◦ n lentv AJ ri ofiuain nldn TAIJ including configurations orbit TAIJI alternative and ) hns cdm fSine,Whn407,China 430071, Wuhan Sciences, of Academy Chinese lentv IATIInetworks LISA-TAIJI Alternative ◦ ,TII (+60 TAIJIc ), ,3, 2, † e-ioHan, Wen-Biao Dtd a 2 2021) 12, May (Dated: ◦ .I h he IATIIcngrtos h LISA- the configurations, LISA-TAIJI three the In ). ◦ nlnto n olctdwt IA,TII ( TAIJIm LISA), with co-located and inclination tram,teSCfrainpaeo h LISA/TAIJI the be of to plane formation designed S/C is interferom- the stable achieve designed arms, To are eter orbit. missions heliocentric TAIJI will the and be missions to LISA space- the The three be by of (S/C). formed to craft Each constellation scheduled triangular are a 2030s. LISA include [21] the including TianQin around and observation launched [20], GW TAIJI band [19], frequency are low band frequency middle [18]. activ- this in present in reviewed The missions briefly 17]. ob- GW [16, be interferome- of DECIGO to or coplanar ities BBO expected and the also co-located of two is ters SGWB the localization by The sky served the sources. of the resolution for will angular constellations the their between improve detection, separations the large of the (SNR) and ratio signal-to-noise the increase ain elyn nteErhlk eicnrcobtwith orbit heliocentric 120 constel- Earth-like three the on proposed deploying missions Both lations DECIGO GW band. and frequency for BBO middle the missions deci-Hz space the in next-generation observation the in planned are detectors independent more demanded. or highly two from servations AJ slaigteErhb 20 by Earth the leading the is assuming By TAIJI [22]. framework Clohessy–Wiltshire the n h at y20 by Earth the ing 2]eautdntokcpblte o h GBobser- SGWB Orlando the and for [26], capabilities Seto network mis- [25], evaluated Seto [27] LISA and single Omiya the sion. to the compared of improvement network localization LISA-TAIJI sky the investigated [24] ,4 5, 4, 1, h pc-on isostreigfrtemilli-Hz the for targeting missions space-borne The h utpeitreoee oprto salso is cooperation interferometer multiple The eAaeyo cecs aghu302,China 310124, Hangzhou Sciences, of Academy se ◦ catcG akrudwt noptimal an with background GW ochastic eaain 1,1] h ubro eetr will detectors of number The 15]. [14, separations cdm fSine,Biig104,China 100049, Beijing Sciences, of Academy e ‡ S ri (+60 orbit ISA egXu, Peng e ilfr ewr bevn GW observing network a form will ies he osbeLS-AJ networks LISA-TAIJI possible three e r nen atr,adLISA-TAIJIc and pattern, antenna ary to n oaiaindetermination polarization and ation AJ,aepandt elaunched be to planned are TAIJI, d ne,Saga 000 China 200030, Shanghai ences, ,Lnhu700,China 730000, Lanzhou y, h,Tia,303 ROC 30013, Taiwan, chu, ,7 4, 7, 6, g109,China 100190, ng ± cnlg (APM), echnology ◦ /Hangzhou), 60 ◦ Ruan , § nlnto n trailing and inclination eua Physics, lecular ◦ n ie Luo Ziren and China p(+60 Ip ihteelpi yemploying by ecliptic the with tal. et ◦ inclination ◦ 2]adWang and [23] n IAi trail- is LISA and ,6 8, 6, 4, − 60 ◦ ¶ tal. et tal. et 2 vation. Wang et al. [28] and Wang et al. [29] estimated as TAIJIp in this work. The preset 20◦ trailing/leading the impact of the joint LISA-TAIJI observation on cos- angle is a practical compromise from the launch vehicle, mological parameter determination. And Wang and Han telemetry capabilities, and the gravitational perturbation [30] demonstrated the observation constraints on the GW reduction [31]. polarizations from the joint observation. The TAIJI orbital configuration could also have other Considering the orbital configuration of the TAIJI mis- choices without (significantly) increasing the launch bud- sion is not fully determined, the merits of the alterna- get. The first alternative is that the constellation forma- tive LISA-TAIJI network are worth evaluating. In this tion is tuned to be 60◦ inclination compared to the ◦ − work, by presetting the LISA orbit is determined, we in- TAIJIp’s +60 and we label it as TAIJIm. Another case vestigate the performances of three possible LISA-TAIJI would be that the TAIJI is co-located and coplanar with networks for different TAIJI orbital selections as shown the LISA which is named as TAIJIc in this work. These in Fig. 1, a) TAIJIp which leading the Earth by 20◦ two orbital configurations are shown in the right panel of and the formation of the constellation is +60◦ inclined∼ as Fig. 1. For the LISA-TAIJIc network, their orientations LISA, b) TAIJIm which also leading the Earth by 20◦ of S/C formation would be co-aligned. For the LISA- and the plane of the S/C has a 60◦ inclination∼ com- TAIJIp, the angle between the two formation planes is pared to the LISA, and c) TAIJIc which− is co-located and 34.5◦, their separation angle will be 40◦ and distance ∼ 8 ∼ coplanar with LISA and trailing the Earth by 20◦. The is 1 10 km. And the angle of the orientations be- ∼ × ◦ deployment and observation for the TAIJI mission∼ from tween the LISA and TAIJIm is around 71 . these orbit choices are expected to be not too much dif- In this work, by employing the numerical orbits, ferent. However, the joint observations with LISA from we investigate the performances of three pairs of alternative TAIJI mission orbits could yield different per- LISA-TAIJI networks (LISA-TAIJIp, LISA-TAIJIm, and formances for the supermassive black hole (SMBH) bi- LISA-TAIJIc) on the detectability for SMBH binaries nary and SGWB observations. We evaluate networks’ an- and the SGWB. These three combinations should cover gular resolutions of sky localization for the SMBH binary, all the possible dual detector scenarios except the LISA observations for the alternative polarizations beyond gen- with a co-located TAIJI having a 60◦ inclination which − eral relativity (GR), and the overlap reduction function would be insipid because neither the large separation for the SGWB observations. In the three pairs combina- benefiting for the compact binary observation nor strong tion, the LISA-TAIJIm network demonstrates the best cross-correlation for the SGWB observation is expected. parameter determinations for the SMBH binary system, The numerical orbits for the TAIJIp and TAIJIc are from and LISA-TAIJIc shows an optimal capability for the our work in [24, 32], and the orbit for TAIJIm is newly SGWB observation. obtained from our optimization method in [32–37]. This paper is organized as follows. In Sec. II, we intro- duce three LISA-TAIJI network configurations and their B. Response formulation of TDI channel joint sensitivities. In Sec. III, we report and compare the results of parameter determinations on the SMBH binary from three LISA-TAIJI networks including the For the space-borne GW missions, time-delay interfer- angular resolution and the alternative polarization con- ometry (TDI) is essential to suppress the laser frequency straints. In Sec. IV, we investigate the overlap reduction noise and achieve targeting sensitivity. The sensitivities functions of the three networks for the SGWB observa- for the different TDI channels have been evaluated nu- tion. We recapitulate our conclusions in Sec. V. (We set merically in our recent works [38, 39]. With implement- G = c = 1 in this work except otherwise stated). ing the TDI, the GW response is combined from the response of each evolved single link. The response functions to the GW tensor polarizations from GR in Doppler measurement were formulated in [40–43]. And the response II. ALTERNATIVE LISA-TAIJI NETWORKS functions for the polarizations beyond the GR were developed in Tinto and da Silva Alves [44]. To keep the A. The LISA and TAIJI orbital configurations integrity of the work, we reiterate the formulas of the response of TDI to the six polarizations as utilized in [30]. The LISA mission is scheduled to be launched in the The GW propagation vector from a source locating at 2030s which includes three S/C forming a 2.5 106 km ecliptic longitude λ and latitude θ (in the solar-system triangle trailing the Earth by 20◦ [19]. The constellation× barycentric coordinates) will be plane has a 60◦ inclination with respect to the ecliptic plane as shown in Fig. 1. The TAIJI mission is proposed kˆ = (cos λ cos θ, sin λ cos θ, sin θ). (1) as a LISA-like orbital configuration with a 3 106 km − arm length [20]. An assumed orbit for the TAIJI× is that The polarization tensors of the GW signal for the +, , the constellation is in front of the Earth by 20◦ and has scalar breathing (b), scalar longitudinal (L), vector x and× the same 60◦ inclination as LISA as shown in the left plot y, combining with the factors of the source’s inclination of Fig. 1, and this TAIJI orbital configuration is labeled angle ι are 3

FIG. 1. The diagram of LISA and TAIJI mission orbital conﬁgurations. The left panel shows the LISA (trailing the Earth by 20◦ and +60◦ inclined with respect to the ecliptic plane) and the TAIJIp (leading the Earth by 20◦ and +60◦ inclined as the∼ LISA). The right panel shows the LISA and anther two optional TAIJI orbital choices which are∼ the TAIJIm (leading the Earth by 20◦ and with a 60◦ inclination) and TAIJIc (coplanar and co-located with LISA). The angle between the LISA and TAIJIp∼ formation planes− is 34.5◦, and the angle between the LISA and TAIJIm formation planes is 71◦. ∼ ∼

1 0 0 2 0 1 0 T 1+cos ι T e+ 1 0 1 0 1 , e× 1 1 0 0 1 i( cos ι), ≡ O · 0− 0 0 · O × 2 ≡ O · 0 0 0 · O × − 1 0 0  0 0 0 T 2 T 2 eb 1 0 1 0 1 sin ι, eL 1 0 0 0 1 sin ι, (2) ≡ O · 0 0 0 · O × ≡ O · 0 0 1 · O × 0 0 1 0 0 0  T T ex 1 0 0 0 1 sin ι cos ι, ey 1 0 0 1 1 i sin ι, ≡ O · 1 0 0 · O × ≡ O · 0 1 0 · O ×     with sin λ cos ψ cos λ sin θ sin ψ sin λ sin ψ cos λ sin θ cos ψ cos λ cos θ − − − − 1 = cos λ cos ψ sin λ sin θ sin ψ cos λ sin ψ sin λ sin θ cos ψ sin λ cos θ , (3) O − cos−θ sin ψ cos−θ cos ψ − sin θ  −   where ψ is the polarization angle. The response to the GW polarization p in the link from S/Ci to j will be

h nîj ep nîj ˆ ˆ yp,ij (f)= · · exp(2πif(Lij + k pi)) exp(2πifk pj) , (4) 2(1 nîj kˆ) × · − · − · h i

wheren ˆij is the unit vector from S/Ci to j, Lij is the the time shift single links, arm length from S/Ci to j, p is the position of the S/Ci i h FX p(f) =( ∆21 + ∆21∆13∆31)y in the solar-system barycentric ecliptic coordinates. , − p,12 h + ( 1 + ∆13∆31)yp,21 − (5) h + (∆31 ∆31∆12∆21)y − p,13 h + (1 ∆12∆21)y , − p,31

where ∆ij = exp(2πifLij). The GW responses in the The first-generation Michelson TDI configuration and Michelson optimal A, E, and T channels are obtained by its corresponding optimal channels are employed to rep- applying [42, 45] resent the performance of each mission. The response of the Michelson-X channel for a specific polarization p in Z X X 2Y+Z X+Y+Z A= − , E= − , T= , (6) the frequency domain will be the sum of the responses in √2 √6 √3 4 where Y and Z channels are obtained from cyclical per- The instantaneous sensitivities to the tensor polariza- mutation of the S/C indexes in the X channel. tions for the LISA and joint LISA-TAIJI networks for ψ = 0, ι = 0, and f = 10 mHz are shown in Fig. 2. As we can see in the upper left plot, the LISA has the op- C. The sensitivity of the networks timum sensitivity around the normal directions (ecliptic latitude 30◦ ) of the triangular formation plane con- 1. The noises in TDI sidering its± 60◦ inclination. As expected from Fig. 1, the antenna pattern of the TAIJIp is shifted by 40◦ ∼ Multiple noise sources will be involved in the process of along the ecliptic latitude with respect to the LISA’s. TDI combinations from the link measurements. For the And their joint sensitivity is shown by the upper right ◦ Michelson-X channels, the expression of measurements panel in Fig. 2. For the TAIJIm, due to its 40 sep- ◦ could be described as [43], aration and 60 inclination with respect to the LISA, its antenna pattern− is not only shifted by 40◦ along the X =[ 31 13 21η12 + 31 13η21 + 31η13 + η31] D D D D D D latitude, also inversed with respect to the ecliptic plane. [η21 + 21η12 + 21 12η31 + 21 12 31η13], And their joint sensitivity of the LISA-TAIJIm is shown − D D D D D D (7) in the lower left panel. As for the TAIJIc case, since the where is a time-delay operator, η(t)= η(t L ), TAIJIc is co-located and coplanar with LISA, the joint Dij Dij − ij ηji are the combined observables from S/Cj to S/Ci [46– LISA-TAIJIc enhanced the LISA’s sensitivity as shown in 48], and the specific expressions for this work are defined the lower right plot. One caveat is that the sensitivity of in [39]. By assuming the dominant laser frequency noises the TAIJI is slightly better than the LISA mission, fully are sufficiently suppressed, the acceleration noise and op- symmetry should not be expected for the joint sensitivity tical path noise in ηij become the primary noise sources plots in Fig. 2. after TDI process. The noise budgets for the acceleration noise Sacc are assumed to be the same for the LISA and TAIJI [19, 49], III. PARAMETER DETERMINATIONS FOR SMBH BINARY COALESCENCE m/s2 0.4mHz 2 f 4 S1/2 =3 10−15 1+ 1+ . acc × √ f 8mHz Hz s s As the most promising GW source for the LISA and (8) TAIJI missions, the SMBH binary is selected to demon- And the optical path noises Sop requirement for two mis- strate the performances of parameter determination from sions are treated slightly different as the LISA and three LISA-TAIJI networks. m 2mHz 4 S1/2 = 10 10−12 1+ , op,LISA × √Hzs f (9) A. Fisher information method 4 1/2 −12 m 2mHz Sop,TAIJI =8 10 1+ . The Fisher information matrix (FIM) is employed in × √Hz s f this investigation to determine the uncertainty of pa- And the power spectrum density (PSD) of a TDI channel rameters from GW observation [50–53, and references Sn,TDI is obtained by implementing the numerical algo- therein]. For a single mission with full six links, the FIM rithm developed in [38, 39]. is combined from three optimal channels (A, E, and T). And the FIM of the joint LISA-TAIJI network is obtained by summing up the FIM from two missions, 2. The joint sensitivities TAIJI ∂h˜TDI ∂h˜TDI The antenna pattern of an interferometer will change Γij = , (12) ∂ξi ∂ξj with the geometric angles Ω(λ,θ,ψ,ι) and the frequency. LISA A,E,T ! X X For a given Ω and frequency, the sensitivities of the

LISA’s A+E+T channel and joint LISA-TAIJI network with at a given mission time could be evaluated respectively ∞ g∗(f)h(f) by, (g h)TDI = 4Re df, (13) | 0 STDI(f) −1/2 Z 2 1/2 FTDI(f, Ω) ˜ SLISA(f, Ω) = | | , (10) where hTDI is the frequency domain GW waveform re-  Sn,TDI(f)  A,E,T sponded in a TDI channel, ξ is the i-th parameter to be X i   −1/2 determined, and STDI(f) is the noise PSD of the corre- TAIJI 2 1/2 FTDI(f, Ω) sponding TDI channel. Sjoint(f, Ω) = | | . (11)  Sn TDI(f)  Considering the source location estimation will be sig- LISA A E T , X X, , niﬁcantly aﬀected by the polarization content of the   5

FIG. 2. The instantaneous sensitivities on the sky map for the LISA mission and LISA-TAIJI networks at ψ = 0, ι = 0, and f = 10 mHz. The single LISA sensitivity is shown by the upper left panel which obtained by using Eq. (10), and the joint LISA-TAIJI network sensitivities are obtained by using Eq. (11). The sensitivity of the LISA-TAIJIp is shown in the upper right panel, the sensitivity of LISA-TAIJIm is shown in the lower left panel, and the sensitivity of the LISA-TAIJIc is shown by the lower right panel. The plots reflect the antenna pattern of the detectors considering the orientation of the S/C formations. The LISA-TAIJIm network achieves a better averaged sensitivity to the different sky directions than the other two networks. (One caveat is that the sensitivity of the TAIJI is slightly better than the LISA mission, fully symmetry on the plots should not be expected). source [2–4], only tensor polarizations from GR are in- parameters for the high SNR ρ 1 will be cluded to investigate the angular resolution of the sky lo- ≫ calization. And 9 parameters are utilized to describe the −1 σi (Γ ) , GW signal and TDI responses from the LISA or LISA- ≃ ii (16) q −1 TAIJI network, which are ecliptic longitude and latitude σij = cov(ξi, ξj ) Γ . ≃ ij (λ, θ), polarization angle ψ, source inclination ι, luminos- ity distance D, the coalescence time and phase (tc, φc), The uncertainty of sky localization for one source is eval- the total mass of binary M and mass ratio q. The GW uated by signal responded by TDI incorporating two polarizations (+ and ) could be described as ∆Ω 2π cos θ σ σ σ2 . (17) × ≃ | | λ θ − λθ q h˜GR,TDI(f) =(F+ + F×)h˜GR, (14) The Monte Carlo simulation is performed for parameter determination by 1000 sources. The (λ, θ) are randomly sampled in the sky sphere, ψ is sampled in [0, 2π] ˜ where hGR is the frequency domain waveform repre- uniformly, cos ι is sampled randomly in [ 1, 1], the merge sented by IMRPhenomPv2 [54]. When the constraints − 5 time tc is randomly in one year. The m1 = 10 M⊙, q = on the alternative GW polarizations are investigated, ad- 1/3 at redshift z = 2 is fixed as we used in [24, 30]. ditional 6 ppE (parametrized post-Einsteinian) param- Considering the SNR is mainly contributed from the bi- eters, (β,b,αb, αL, αx, αx), are employed to qualify the nary coalescing stage, the 30 days observation before the deviations from the GR as developed in [55], and the merge is simulated to perform the investigation. waveform will be explained in Eq. (18). The variance-covariance matrix of the parameters is calculated by B. Sky localization of the networks

ρ≫1 The cumulative histograms of SNR from the LISA and ∆ξ ∆ξ = Γ−1 + (ρ−1) Γ−1 . (15) h i j i ij O ≃ ij LISA-TAIJI networks are shown in the left panel of Fig. 3. Compared to the single LISA mission, three LISA- The standard deviations σi and correlation σij of the TAIJI networks achieve more than √2 SNR by imple- 6

1.0 1 LISA LISA-TAIJIp -TAIJIp 0.8 LISA-TAIJIm LISA-TAIJIm LISA-TAIJIc LISA-TAIJIc

0.6

fraction of events 0.4

ΔΔΔ

Cumulative 0.2

0.0 1 1 102 103 1 1 1 1 1 1 SNR ΔΩ

FIG. 3. The cumulative histograms of SNR (left panel) and angular resolutions of the sky localization (right panel) from the LISA and three LISA-TAIJI networks. Three LISA-TAIJI networks could achieve more than √2 times SNR of the LISA mission. In the three networks, the SNR from LISA-TAIJIm has the most concentrated distribution because of the most averaged antenna pattern as shown in Fig. 2. The distribution of SNR from LISA-TAIJIc has a long tail since the two missions share the same sensitive/insensitive areas. On the right plot, the resolution of sky localization by LISA-TAIJIc is approximately twice better than the single LISA as the attribute to the SNR increase. The angular resolution of the LISA-TAIJIp is better than LISA-TAIJIc beneﬁting from the long separation. The performance of LISA-TAIJIm is better than LISA-TAIJIp because the angle of the LISA and TAIJIm formation planes is 71◦ and their orientation could yield more collaborative antenna pattern.

2 2 2 menting the quadratic sum (ρjoint = ρLISA + ρTAIJI) con- C. Observation for GW polarizations sidering the TAIJI is more sensitive to LISA in the selected GW frequency band. In the three networks, the The detector responded GW signal is modiﬁed as fol- LISA-TAIJIc has a larger range of the SNR distribu- lows to incorporate alternative polarization beyond GR tion with a longer tail, because the co-aligned detectors [30, 55], are sensitive/insensitive to the same directions and leav- ing the common optimal/blind areas. The LISA-TAIJIm ˜ b+5 hppE,TDI(f)= (FTDI,+ + FTDI,×)(1 + Cβu2 ) shows the most concentrated SNRs values compared to +α F + α F the two other networks since their joint antenna pattern b TDI,b L TDI,L (18) b 2iβu2 is more averaged on the sky map as shown in Fig. 2. +αxFTDI,x + αyFTDI,y] h˜GR e ,

where C is the function of b and is deﬁned by Eq. (11) in the Erratum [56] of [55], h˜GR is the GW waveform repre- 1/3 sented by IMRPhenomPv2 [54], and u2 (π f) . In The angular resolutions of the sky localization from this investigation, we choose b = 3 which≡ correspondsM the LISA and LISA-TAIJI networks are shown in the to the massive graviton theory [57–63],− and set β =0.01 right panel of Fig. 3. For the LISA-TAIJIc network, the which is from the rough boundary constrained in [64]. uncertainties of sky localization are improved by more The other four ppE parameters tuning the amplitudes than 2 times compared to the single LISA mission which of the alternative polarizations are set to be zeros, (αb, should be due to the more than √2 times SNR from the αL, αx, αy)=(0, 0, 0, 0). Although this speciﬁc selection network. Compared to the LISA-TAIJIc, the joint obser- could not represent all other gravity theories, we have vation from LISA-TAIJIp demonstrate the more than 2 demonstrated that the measurement on the ppE param- orders improvement on the localization resolution which eters could also be similarly improved by the LISA-TAIJI should be mainly attributed to the long baseline sepa- network for other value choices [30]. rations between the LISA and TAIJIp. On the other The constraints on the ppE parameters α for the scalar side, the LISA-TAIJIm yield a better capability on lo- (upper panel) and vector (lower panel) polarizations are cating the source than LISA-TAIJIp because the TAI- shown in Fig. 4. Similar to the results achieved for JIm’s formation plane is 71◦ with respect to LISA’s, and the sky localization, the constraint on α from the LISA- its antenna pattern could better compensate the LISA’s TAIJIc is more than √2 times better than single LISA insensitive directions. which is an attribute from the increase of SNR. The 7

1 1

fraction of events fraction of events

LISA LISA LISA-TAIJIp LISA-TAIJIp LISA-TAIJIm LISA-TAIJIm LISA-TAIJIc LISA-TAIJIc 1 11 1 1 11 1 α α

1 1

fraction of events fraction of events

LISA LISA LISA-TAIJIp LISA-TAIJIp LISA-TAIJIm LISA-TAIJIm LISA-TAIJIc LISA-TAIJIc 1 1 11 1 1 1 11 1 α α

FIG. 4. The cumulative histograms of the constraints on the amplitudes of scalar and vector polarizations. The result from the LISA-TAIJIc is more than √2 times better than the single LISA which is an attribute from the increase of SNR. The LISA-TAIJIp, with long-baseline separation, could significantly improve the constraints on polarizations. And LISA-TAIJIm achieves the best constraints on the polarizations in three LISA-TAIJI configurations as the benefit of better antenna pattern cooperation.

LISA-TAIJIp, with large separation, could improve the function is introduced to indicate the cross-correlation constraints on the polarizations significantly. And LISA- between a pair of detectors [66]. Whelan et al. [67] cal- TAIJIm could achieve the best constraints on the polar- culated the overlap reduction functions for the two LIGO izations in the three LISA-TAIJI configurations as the detectors and GEO. Omiya and Seto [25], Seto [26], and benefit of better antenna pattern cooperation. Orlando et al. [27] specified the overlap function of the LISA-TAIJIp network for optimal TDI channels and alternative GW polarizations. Schmitz [17] reviewed the IV. OVERLAP REDUCTION FUNCTION OF detectability of the ground- and space-based detectors THE LISA-TAIJI NETWORKS for the stochastic GW background. For a single LISA-like mission with full six measure- The response of the detector network to the stochastic ment links, the optimal TDI channels could be treated as background GW signal will depend on the locations and three equivalent interferometers. The observation from orientations of the interferometers. Flanagan [65] evalu- these TDI channels could be used to discriminate the ated the sensitivities of the ground-based GW interferom- stochastic GW background from the instrument noise eters to the stochastic background. An overlap reduction [68]. And the motion of the detectors may also help 8

-TAIJIp -TAIJIp LISA-TAIJIm LISA-TAIJIm LISA-TAIJIc LISA-TAIJIc ( ( ( ( 1 1 1 11 1 1 1 11

scalar (b) LISA-TAIJIp -TAIJIp LISA-TAIJIm LISA-TAIJIm LISA-TAIJIc LISA-TAIJIc

Reduction Function ) ( ) ( ) ( ) ( 1 1 1 11 1 1 1 11 (

FIG. 5. The overlap reduction functions of three LISA-TAIJI networks for different polarization modes. The overlap function of the LISA-TAIJIc network is unity for the frequency lower than 10 mHz which is optimal for the SGWB observation, and it changes the sign during the two detectors’ characteristic frequencies gap [50 mHz, 60 mHz]. The overlap function of the LISA- 8 8 TAIJIp/TAIJIm are close to zero around 1.5 mHz due to their 1 10 km separations (fcrit = c/(2 1 10 km) 1.5 mHz). The LISA-TAIJIm network has a worse cross-correlation than the× LISA-TAIJIp because of the more× misaligned× orientation≃ . to resolve the background, especially for the anisotropic ization angle ψ is set to be zero, and the inclination ι is signal [16]. The LISA and TAIJI could form an ideal set to be optimal for each polarization mode. network to separate the cosmological SGWB signal from The overlap reduction functions for the different LISA- other stochastic processes such as the instrument noise TAIJI networks for different polarizations are shown in and astrophysical foreground. Fig. 5. Since the orientations of the LISA and TAIJIc To characterize the cross-correlation between LISA are aligned and locations of them are coplanar and at the and different TAIJI orbital configurations, their overlap same location, the overlap function of the LISA-TAIJIc reduction functions are calculated for different polariza- network is unity for the frequency lower than 10 mHz tions, which indicates the strong correlation between the LISA and TAIJIc detectors, and the network is optimal for the κ γ (f)= dn F a (f, n) F b (f, n) SGWB observation for all polarization modes. We also ab,p 4π TDI,p TDI,p A E T A E T notice that their overlap reduction functions change the Z X, , X, , (19) sign during the two detectors’ characteristic frequencies a gap [ c = 50 mHz, c = 60 mHz] (c is the speed where FTDI,p is the response function to the polarization 2LTAIJI 2LLISA mode p (tensor, vector, scalar breathing, and scalar lon- of the light in this section). gitudinal) in the TDI channel from the mission a, and Considering the 1 108 km separation between LISA × κ is the normalization factor to make γab = 1 when the and TAIJIp/TAIJIm, we can see the overlap reduction two detectors are co-aligned and co-localized. The polar- functions quickly approach to zero around a critical fre- 9

8 quency fcrit c/(2 1 10 km) 1.5 mHz [16], and ≃ × × ≃ 90 oscillate and decay with the frequency increase. The γab from the LISA-TAIJIp pair is higher than the value of LISA-TAIJIm the LISA-TAIJIm because the orientation of the LISA is 75 more aligned with the TAIIJIp (34.5◦) than the TAIJIm ◦ 60 (71 ) to make TAIJIp have a stronger correlation with BBO/DECIGO the LISA. Therefore, in the three LISA-TAIJI networks, the LISA-TAIJIm should be the relatively worst config- 45 LISA-TAIJIp uration for the SGWB detections. There will be a trade-off for a LISA-like GW detec- 30 tor network to observe the compact binary system and 15 SGWB. A long baseline for detectors deployment will plus-plus promote the accuracy of parameter estimation for the LISA-TAIJIc plus-minus SMBH binary. However, the frequency for the detectable Angle between constellation planes (deg) 0 SGWB band would be lowered referring to the critical 0 20 40 60 80 100 120 140 160 180 frequency, Separation angle of two missions (deg) c c 0.5 mHz f = = (20) crit ǫ ǫ FIG. 6. The angle of the formation planes varying with the 2d 2 2AU sin 2 ≃ sin 2 × separation angle of two constellations. The curve of the plus- where d is the distance between two detectors, ǫ is the plus indicates the two constellations having +60◦ inclinations separation angle formed by the two lines connecting the with respect to the ecliptic plane, and the curve of the plus- ◦ Sun and detectors. minus indicates that one constellation having a +60 inclination and another having a 60◦ inclination. On the other hand, the angle between the constellation − plane also changes with the separation angle as shown in Fig. 6. Considering the orientation of the plane with an angle closing to 90◦ should be helpful to resolve the source parameters, the composition of a +60◦ with a 60◦ in- orientation with the LISA in the three TAIJI orbital con- clination could more cooperative than the two− missions figurations. with the same inclination for a separation angle small One lesson from this evaluation of three LISA-TAIJI than 90◦. For the BBO or DECIGO mission, the con- networks is that the parameter resolution of the compact stellations are planned to be separated by 120◦, and two binary coalescences will be impacted by the SNR, the dis- options could be considered for their orientation deploy- tance of detector separation, and the cooperative orienta- ment, 83◦ or 51◦ with respect to another formation plane tions of the detectors. The next generation space-based as tagged in Fig. 6. GW detectors, both DECIGO and BBO, are proposed to be the LISA-like orbit with multiple constellations to de- V. CONCLUSIONS tect the relic GW left by the Big Bang, intermediate-mass black holes, etc in the deci-Hz frequency band [14, 15]. The parameter resolution improvements have been per- In this work, we investigate the performances of three formed for the compact binaries in [14] as the results alternative LISA-TAIJI networks on the sky localizations of the multiple interferometers and long-baseline. The and polarization observations from the SMBH binary and orientations combinations of the constellations are also the overlap reduction function for the stochastic gravi- worth to be evaluated for the targeting sources. tational wave background observation. For the SMBH binary system, compared to the single LISA mission, the Beyond the LISA-like orbital formation, various space co-located and coplanar LISA-TAIJIc network ordinarily missions are proposed to arrange S/C equally on a planet improve the SNR and parameter resolution by a factor orbit in order to observe GW in the micro-Hz band, for of √2 times. With a 1 108 km separation, the joint ob- instance, ASTROD-GW [69], Folkner mission [70], and servations from LISA× and TAIJIp significantly improve µAres [71], etc. The ASTROD-GW is initially proposed parameter determinations for the SMBH binary than the to deploy 3 S/C around the Lagrange points L3, L4, and LISA-TAIJIc. The LISA-TAIJIm network demonstrates L5 of the Sun-Earth system, and an extended deployment a better capability to determine the sky location and po- could be 6 S/C to form two triangular interferometers to larizations than the LISA-TAIJIp network as the bene- enhance the sensitivity to the SGWB [69]. The µAres will fits of more misaligned orientation and complementary place the two orthogonal triangle interferometers with antenna pattern. For the detectability for the stochastic respect to the Mars (or Earth/Venus) orbit. The trade- gravitational wave background, the LISA-TAIJIc would off of detectability from various deployments could also have optimal performance as the benefit of the coplanar be explored to balance the GW observations from the formation and co-location, the TAIJIm present the worst compact binary systems and the cosmological stochastic cross-correlation with the LISA due to the least aligned background. 10

ACKNOWLEDGMENTS of the Chinese Academy of Sciences under Grants No. XDA15021102. This work made use of the High Perfor- mance Computing Resource in the Core Facility for Ad- This work was supported by NSFC Nos. 12003059 vanced Research Computing at Shanghai Astronomical and 11773059, the Strategic Priority Research Program Observatory.

[1] B. P. Abbott et al. (LIGO Scientific, Virgo), Obser- grounds using data from Advanced LIGO’s and Ad- vation of Gravitational Waves from a Binary Black vanced Virgo’s first three observing runs, (2021), Hole Merger, Phys. Rev. Lett. 116, 061102 (2016), arXiv:2103.08520 [gr-qc]. arXiv:1602.03837 [gr-qc]. [14] J. Crowder and N. J. Cornish, Beyond LISA: [2] B. P. Abbott et al. (LIGO Scientific, Virgo), Exploring future gravitational wave missions, GW170814: A Three-Detector Observation of Phys. Rev. D 72, 083005 (2005), arXiv:gr-qc/0506015. Gravitational Waves from a Binary Black Hole [15] S. Kawamura et al., The Japanese space Coalescence, Phys. Rev. Lett. 119, 141101 (2017), gravitational wave antenna DECIGO, arXiv:1709.09660 [gr-qc]. Class. Quant. Grav. 23, S125 (2006). [3] B. P. Abbott et al. (LIGO Scientific, Virgo), GW170817: [16] J. D. Romano and N. J. Cornish, Detection meth- Observation of Gravitational Waves from a Binary Neu- ods for stochastic gravitational-wave backgrounds: tron Star Inspiral, Phys. Rev. Lett. 119, 161101 (2017), a unified treatment, Living Rev. Rel. 20, 2 (2017), arXiv:1710.05832 [gr-qc]. arXiv:1608.06889 [gr-qc]. [4] B. P. Abbott et al. (LIGO Scientific, [17] K. Schmitz, New Sensitivity Curves for Gravitational- Virgo), Tests of General Relativity with Wave Signals from Cosmological Phase Transitions, GW170817, Phys. Rev. Lett. 123, 011102 (2019), JHEP 01, 097, arXiv:2002.04615 [hep-ph]. arXiv:1811.00364 [gr-qc]. [18] W.-T. Ni, Mid-frequency gravitational wave detection [5] T. Akutsu et al. (KAGRA), Overview of KA- and sources, Int. J. Mod. Phys. D 29, 1902005 (2020), GRA : KAGRA science 10.1093/ptep/ptaa120 (2020), arXiv:2004.05590 [gr-qc]. arXiv:2008.02921 [gr-qc]. [19] P. Amaro-Seoane, H. Audley, S. Babak, and et al (LISA [6] B. P. Abbott et al. (KAGRA, LIGO Scientific, Virgo), Team), Laser Interferometer Space Antenna, arXiv e- Prospects for observing and localizing gravitational-wave prints , arXiv:1702.00786 (2017). transients with Advanced LIGO, Advanced Virgo and [20] W.-R. Hu and Y.-L. Wu, The Taiji Program in Space KAGRA, Living Rev. Rel. 23, 3 (2020). for gravitational wave physics and the nature of gravity, [7] B. F. Schutz, Networks of gravita- Natl. Sci. Rev. 4, 685 (2017). tional wave detectors and three figures of [21] J. Luo et al. (TianQin Team), TianQin: merit, Class. Quant. Grav. 28, 125023 (2011), a space-borne gravitational wave detec- arXiv:1102.5421 [astro-ph.IM]. tor, Class. Quant. Grav. 33, 035010 (2016), [8] B. P. Abbott et al. (LIGO Scientific, Virgo), Up- arXiv:1512.02076 [astro-ph.IM]. per Limits on the Stochastic Gravitational-Wave [22] S. V. Dhurandhar, K. Rajesh Nayak, S. Koshti, Background from Advanced LIGO’s First Ob- and J. Y. Vinet, Fundamentals of the LISA sta- serving Run, Phys. Rev. Lett. 118, 121101 (2017), ble flight formation, Class. Quant. Grav. 22, 481 (2005), [Erratum: Phys.Rev.Lett. 119, 029901 (2017)], arXiv:gr-qc/0410093. arXiv:1612.02029 [gr-qc]. [23] W.-H. Ruan, C. Liu, Z.-K. Guo, Y.-L. Wu, and R.-G. Cai, [9] B. P. Abbott et al. (LIGO Scientific, Virgo), The LISA-Taiji network, Nature Astron. 4, 108 (2020), Search for Tensor, Vector, and Scalar Polar- arXiv:2002.03603 [gr-qc]. izations in the Stochastic Gravitational-Wave [24] G. Wang, W.-T. Ni, W.-B. Han, S.-C. Yang, Background, Phys. Rev. Lett. 120, 201102 (2018), and X.-Y. Zhong, Numerical simulation of sky arXiv:1802.10194 [gr-qc]. localization for LISA-TAIJI joint observation, [10] B. P. Abbott et al. (LIGO Scientific, Virgo), Phys. Rev. D 102, 024089 (2020), arXiv:2002.12628. Search for the isotropic stochastic background [25] H. Omiya and N. Seto, Searching for anoma- using data from Advanced LIGO’s second ob- lous polarization modes of the stochastic serving run, Phys. Rev. D 100, 061101 (2019), gravitational wave background with LISA arXiv:1903.02886 [gr-qc]. and Taiji, Phys. Rev. D 102, 084053 (2020), [11] R. Abbott et al. (LIGO Scientific, Virgo, KAGRA), Up- arXiv:2010.00771 [gr-qc]. per Limits on the Isotropic Gravitational-Wave Back- [26] N. Seto, Gravitational Wave Background Search ground from Advanced LIGO’s and Advanced Virgo’s by Correlating Multiple Triangular Detectors in Third Observing Run, (2021), arXiv:2101.12130 [gr-qc]. the mHz Band, Phys.Rev.D 102, 123547 (2020), [12] R. Abbott et al. (LIGO Scientific, Virgo, KAGRA), arXiv:2010.06877 [gr-qc]. Constraints on cosmic strings using data from the [27] G. Orlando, M. Pieroni, and A. Ricciardone, Measur- third Advanced LIGO-Virgo observing run, (2021), ing Parity Violation in the Stochastic Gravitational arXiv:2101.12248 [gr-qc]. Wave Background with the LISA-Taiji network, (2020), [13] R. Abbott et al. (LIGO Scientific, Virgo, KA- arXiv:2011.07059 [astro-ph.CO]. GRA), Search for anisotropic gravitational-wave back- [28] R. Wang, W.-H. Ruan, Q. Yang, Z.-K. Guo, R.- 11

G. Cai, and B. Hu, Hubble parameter estimation via arXiv:1010.1302 [gr-qc]. dark sirens with the LISA-Taiji network, (2020), [45] T. A. Prince, M. Tinto, S. L. Larson, arXiv:2010.14732 [astro-ph.CO]. and J. W. Armstrong, The LISA optimal [29] L.-F. Wang, S.-J. Jin, J.-F. Zhang, and X. Zhang, sensitivity, Phys. Rev. D 66, 122002 (2002), Forecast for cosmological parameter estimation with arXiv:gr-qc/0209039 [gr-qc]. gravitational-wave standard sirens from the LISA-Taiji [46] M. Otto, G. Heinzel, and K. Danzmann, TDI and clock network, (2021), arXiv:2101.11882 [gr-qc]. noise removal for the split interferometry configuration [30] G. Wang and W.-B. Han, Observing grav- of LISA, Class. Quant. Grav. 29, 205003 (2012). itational wave polarizations with the LISA- [47] M. Otto, Time-Delay Interferometry Simulations for the TAIJI network, Phys. Rev. D 103, 064021 (2021), Laser Interferometer Space Antenna, Ph.D. thesis, Leib- arXiv:2101.01991 [gr-qc]. niz U., Hannover (2015). [31] LISA Study Team, LISA (Laser Interferometer Space [48] M. Tinto and O. Hartwig, Time-Delay Antenna): A Cornerstone Mission for the Observation Interferometry and Clock-Noise Cali- of Gravitational Waves, Tech. Rep. 11 (ESA-SCI, 2000) bration, Phys. Rev. D 98, 042003 (2018), system and Technology Study Report. arXiv:1807.02594 [gr-qc]. [32] G. Wang and W.-T. Ni, Numerical simulation [49] Z. Luo, Z. Guo, G. Jin, Y. Wu, and W. Hu, of time delay interferometry for TAIJI and new A brief analysis to Taiji: Science and technology, LISA, Res. Astron. Astrophys. 19, 058 (2019), Results in Physics 16, 102918 (2020). arXiv:1707.09127 [astro-ph.IM]. [50] C. Cutler and E.´ E. Flanagan, Gravitational waves [33] G. Wang and W.-T. Ni, Time- from merging compact binaries: How accurately delay Interferometry for ASTROD-GW, can one extract the binary’s parameters from the Chin. Astron. Astrophys. 36, 211 (2012), and refer- inspiral waveform?, Phys. Rev. D 49, 2658 (1994), ences therein. arXiv:gr-qc/9402014 [gr-qc]. [34] G. Wang and W.-T. Ni, Numermcal sim- [51] C. Cutler, Angular resolution of the LISA gravi- ulation of time delay interferometry for tational wave detector, Phys. Rev. D57, 7089 (1998), NGO/eLISA, Class. Quant. Grav. 30, 065011 (2013), arXiv:gr-qc/9703068 [gr-qc]. arXiv:1204.2125 [gr-qc]. [52] M. Vallisneri, Use and abuse of the Fisher [35] G. Wang and W.-T. Ni, Orbit optimization for ASTROD- information matrix in the assessment of GW and its time delay interferometry with two arms gravitational-wave parameter-estimation using CGC ephemeris, Chin. Phys. B22, 049501 (2013), prospects, Phys. Rev. D77, 042001 (2008), arXiv:1205.5175 [gr-qc]. arXiv:gr-qc/0703086 [GR-QC]. [36] S. V. Dhurandhar, W. T. Ni, and G. Wang, Nu- [53] K. A. Kuns, H. Yu, Y. Chen, and R. X. Ad- merical simulation of time delay interferometry for a hikari, Astrophysics and cosmology with a deci- LISA-like mission with the simplification of having hertz gravitational-wave detector: TianGO, (2019), only one interferometer, Adv. Space Res. 51, 198 (2013), arXiv:1908.06004 [gr-qc]. arXiv:1102.4965 [gr-qc]. [54] S. Khan, S. Husa, M. Hannam, F. Ohme, M. Pürrer, [37] G. Wang and W.-T. Ni, Orbit optimiza- X. J. Forteza, and A. Bohé, Frequency-domain grav- tion and time delay interferometry for in- itational waves from nonprecessing black-hole bina- clined ASTROD-GW formation with half-year ries. II. A phenomenological model for the ad- precession-period, Chin. Phys. B24, 059501 (2015), vanced detector era, Phys. Rev. D 93, 044007 (2016), arXiv:1409.4162 [gr-qc]. arXiv:1508.07253 [gr-qc]. [38] G. Wang, W.-T. Ni, and W.-B. Han, Revisiting time [55] K. Chatziioannou, N. Yunes, and N. Cornish, Model- delay interferometry for unequal-arm LISA and TAIJI, independent test of general relativity: An ex- (2020), arXiv:2008.05812 [gr-qc]. tended post-Einsteinian framework with complete [39] G. Wang, W.-T. Ni, W.-B. Han, and C.-F. Qiao, Algo- polarization content, Phys. Rev. D 86, 022004 (2012), rithm for TDI numerical simulation and sensitivity inves- arXiv:1204.2585 [gr-qc]. tigation, (2020), arXiv:2010.15544 [gr-qc]. [56] K. Chatziioannou, N. Yunes, and N. Cornish, Erratum: [40] F. B. Estabrook and H. D. Wahlquist, Response of Model-independent test of general relativity: An ex- Doppler spacecraft tracking to gravitational radiation., tended post-einsteinian framework with complete po- General Relativity and Gravitation 6, 439 (1975). larization content [phys. rev. d 86, 022004 (2012)], [41] H. Wahlquist, The Doppler response to grav- Phys. Rev. D 95, 129901 (2017). itational waves from a binary star source., [57] C. M. Will, Bounding the mass of the gravi- General Relativity and Gravitation 19, 1101 (1987). ton using gravitational wave observations of inspi- [42] M. Vallisneri, J. Crowder, and M. Tinto, Sensitivity ralling compact binaries, Phys. Rev. D 57, 2061 (1998), and parameter-estimation precision for alternate LISA arXiv:gr-qc/9709011. configurations, Class. Quant. Grav. 25, 065005 (2008), [58] C. M. Will and N. Yunes, Testing al- arXiv:0710.4369 [gr-qc]. ternative theories of gravity using [43] M. Vallisneri and C. R. Galley, Non-sky-averaged LISA, Class. Quant. Grav. 21, 4367 (2004), sensitivity curves for space-based gravitational-wave arXiv:gr-qc/0403100. observatories, Class. Quant. Grav. 29, 124015 (2012), [59] E. Berti, A. Buonanno, and C. M. Will, Testing general arXiv:1201.3684 [gr-qc]. relativity and probing the merger history of massive black [44] M. Tinto and M. E. da Silva Alves, LISA Sensitiv- holes with LISA, Class. Quant. Grav. 22, S943 (2005), ities to Gravitational Waves from Relativistic Met- arXiv:gr-qc/0504017. ric Theories of Gravity, Phys. Rev. D 82, 122003 (2010), 12

[60] A. Stavridis and C. M. Will, Bounding the terferometer gravitational wave observatory (LIGO) mass of the graviton with gravitational waves: to a stochastic background, and its dependence on Eﬀect of spin precessions in massive black the detector orientations, Phys. Rev. D 48, 2389 (1993), hole binaries, Phys. Rev. D 80, 044002 (2009), arXiv:astro-ph/9305029. arXiv:0906.3602 [gr-qc]. [66] N. Christensen, Measuring the stochastic gravitational [61] K. G. Arun and C. M. Will, Bounding the mass radiation background with laser interferometric anten- of the graviton with gravitational waves: Eﬀect nas, Phys. Rev. D 46, 5250 (1992). of higher harmonics in gravitational waveform [67] J. T. Whelan, W. G. Anderson, M. Casquette, templates, Class. Quant. Grav. 26, 155002 (2009), M. C. Diaz, I. S. Heng, M. McHugh, J. D. Ro- arXiv:0904.1190 [gr-qc]. mano, J. Torres Charlie W., R. M. Trejo, and [62] D. Keppel and P. Ajith, Constraining the A. Vecchio, Progress on stochastic background search mass of the graviton using coalescing black- codes for LIGO, Class. Quant. Grav. 19, 1521 (2002), hole binaries, Phys. Rev. D 82, 122001 (2010), arXiv:gr-qc/0110019. arXiv:1004.0284 [gr-qc]. [68] M. R. Adams and N. J. Cornish, Discriminating be- [63] K. Yagi and T. Tanaka, Constraining alternative tween a Stochastic Gravitational Wave Background theories of gravity by gravitational waves from and Instrument Noise, Phys. Rev. D 82, 022002 (2010), precessing eccentric compact binaries with LISA, arXiv:1002.1291 [gr-qc]. Phys. Rev. D 81, 064008 (2010), [Erratum: Phys.Rev.D [69] W.-T. Ni, ASTROD-GW: Overview and 81, 109902 (2010)], arXiv:0906.4269 [gr-qc]. Progress, Int. J. Mod. Phys. D 22, 1341004 (2013), [64] N. Cornish, L. Sampson, N. Yunes, and F. Pre- arXiv:1212.2816 [astro-ph.IM]. torius, Gravitational Wave Tests of General Rel- [70] J. Baker et al., Space Based Gravitational ativity with the Parameterized Post-Einsteinian Wave Astronomy Beyond LISA, (2019), Framework, Phys. Rev. D 84, 062003 (2011), arXiv:1907.11305 [astro-ph.IM]. arXiv:1105.2088 [gr-qc]. [71] A. Sesana et al., Unveiling the Gravitational [65] E. E. Flanagan, The Sensitivity of the laser in- Universe at mu-Hz Frequencies, (2019), arXiv:1908.11391 [astro-ph.IM].\