Science with the TianQin Observatory: Preliminary Results on Stochastic Gravitational-Wave Background Zheng-Cheng Liang,1, 2, ∗ Yi-Ming Hu,2, y Yun Jiang,2, 3, z Jun Cheng,2, x Jian-dong Zhang,2, { and Jianwei Mei2, ∗∗ 1MOE Key Laboratory of Fundamental Physical Quantities Measurements, Hubei Key Laboratory of Gravitation and Quantum Physics, School of Physics, Huazhong University of Science and Technology, Wuhan 430074, China 2MOE Key Laboratory of TianQin Mission, TianQin Research Center for Gravitational Physics & School of Physics and Astronomy, Frontiers Science Center for TianQin, CNSA Research Center for Gravitational Waves, Sun Yat-sen University (Zhuhai Campus), Zhuhai 519082, China 3Max-Planck-Institut f¨urKernphysik, Saupfercheckweg 1, 69117 Heidelberg, Germany (Dated: July 20, 2021) In this work, we study the prospect of detecting the stochastic gravitational-wave background with the TianQin observatory. We consider both astrophysical-origin and cosmological-origin sources, in- cluding stellar-mass binary black holes, binary neutron stars, Galactic white dwarves, inflation, first order phase transition, and cosmic defects. For the detector configurations, we considered TianQin, TianQin I+II and TianQin + LISA. We studied the detectability of stochastic gravitational-wave backgrounds with the assumed methods of both cross-correlation and null channel, and present the corresponding power-law integrated sensitivity curves. We introduce the definition of the \joint foreground" with a network of detectors. With the joint foreground, the number of resolved double white dwarves in the Galaxy will be increased by 5% ∼ 22% compared with simple combination of individual detectors. The astrophysical background from the binary black holes and the binary neutron stars under the theoretical models are predicted to be detectable with signal-to-noise ratio of around 10 after five years operation. As for the cosmological sources, their models are highly uncertain, and we only roughly estimate the detection capability under certain cases. I. INTRODUCTION the event rate, mass distribution and formation mecha- nisms [16{18]. As for the DWDs in the Galaxy, one can The gravitational-wave (GW) emission from a large study the spatial structure by its anisotropy [19]. In ad- number of independent and unresolvable sources can not dition, some physical pictures of the early Universe may be detected individually by the GW detectors. Nev- be hidden in the cosmological SGWBs [20{22]. ertheless, their incoherent superposition would form a Specifically, the collection of Galaxy DWDs could ex- stochastic GW background (SGWB), which could be de- ceed the noise level of space-borne GW detectors [23]. tected as a bulk [1{3]. Many different mechanisms can Instead of forming a background, such signal is generally contribute to the SGWB, by which the SGWB can be classified as a foreground [24{27]. The capability to de- roughly divided into two categories: astrophysical-origin tect the other isotropic SGWB would not be significantly and cosmological-origin [4,5]. weakened when the foreground is suitably modeled [28]. The astrophysical SGWB mainly contains the GW Since the strength of foreground is comparable with the emission from plenty of compact binaries [6]. For detector noise, in the process of data analysis for other space-borne GW detectors, such source could be signals, the all-sky integrated foreground could be treated double white dwarfs (DWDs) [7,8], massive black as a part of the noise [29{34]. hole binarys (MBHBs) [9], stellar-mass binary black There exist a few programs for the space-borne GW holes (SBBHs) [10], and extreme-mass-ratio inspirals detectors [6, 35{38]. In this work, we constrain our focus (EMRIs) [11, 12]. On the other hand, the cosmological on TianQin [35, 39, 40] (and LISA [6] when a network is SGWB for space-borne GW detectors contains the GWs considered). TianQin is expected to explore GW astron- from physical processes linked to the early Universe, like omy as well as fundemantal physics during its operation the inflation [13], first order phase transition (PT) [14] [8{12, 41, 42]. For such detectors, the laser interferom- and cosmic defects [15]. A detection of the SGWB would etry is used to combine three satellites into a six-link have important implications in either astrophysics or cos- arXiv:2107.08643v1 [astro-ph.CO] 19 Jul 2021 triangle GW detector, with which one can build up three mology on the fundamental physics. When exploring the independent data channels. A link is formed when the SGWB of extragalactic compact binaries, one can limit laser emitted from one satellite is received by the other, which is further used to construct an interferometer. When extracting the GW from the resolved sources such as SBBHs, a common practice is to use the matched ∗ [email protected] filter method [43, 44], where a template bank with various y [email protected] z [email protected] waveforms are used to compare with the data. However, x [email protected] due to the stochastic nature, the SGWB has no definite { [email protected] waveforms. Thus, in the realm of SGWB detection, the ∗∗ [email protected] optimal filter method is used instead [45]. 2 For the space-borne GW missions, the laser phase noise where h:::i denotes ensemble average [61]. In the plane is usually orders of magnitude higher than other noises, wave expansion, the metric perturbation can be ex- while the time delay interferometer (TDI) can beat down pressed as [62, 63]: such laser phase noise [46{50]. The outcome, however, is Z 1 Z ^ i2πf(t−k^·~x=c) that instead of the Michelson channel, analysis should hab(t; ~x) = df dΩk^ehab(f; k)e ; (4) be performed on the TDI channels, among which the A −∞ S2 and E channels are conventional channels, while the T ^ P ^ P ^ P ^ where ehab(f; k) = P =+;× ehP (f; k)eab(k), with eab(k) channel is signal-insensitive, and is also referred to as being the polarization tensor[58, 64, 65]. the null channel or noise monitoring channel. The detector response can be expressed as the con- In the search of SGWB, one could cross-correlate the volution of metric perturbations h (t; ~x) and the im- outputs from multiple channels [51{53]. At the first ab pulse response of the detector Fab(t; ~x)[66]: h(t) = glance, such method could in principle be applied to the ab F (t; ~x) ∗ hab(t; ~x), through which the output in the A/E channels. However, the GW signals projected to frequency-domain is derived: the A and E channels belong to orthogonal polarizations. Z Therefore, cross correlation method can not be applied X P ^ ^ −i2πfk^·~x=c eh(f) = dΩk^ F (f; k)ehP (f; k)e ; 2 on the TDI channels [54]. Fortunately, the T channel S P =+;× can provide information about the noise spectrum, and (5) the SGWB can be identified in the A and E channels where F P (f; k^) = eP (k^)F ab(f; k^) and with the comparison of the T channel through the null ab Z 1 Z channel method [55, 56]. We calculate the responses ab −i2πfk^·~x=c 3 ab F (f; k^) = e d τ d y F (τ; ~y) of SGWB under different configurations with both the −∞ cross-correlation and null channel methods, as well as ^ × e−i2πf(τ−k·~y=c): (6) the spectra from various types of sources. The outline of the paper is as follows. In Sec.II, we re- Taking the Michelson channel as an example, we have view the fundamental for the SGWB to provide a theoret- 1 F ab(f; k^) = (uaub T (f; u^ ; k^) − uaub T (f; u^ ; k^)); (7) ical framework for the following calculation. In Sec. III, M 2 1 1 1 2 2 2 we introduce the methodology to realize the detection for whereu ^i are the unit vectors for the two arms of the the SGWB and show the power-law integrated (PI) sen- a Michelson channel and ui denotes the a-component of the sitivity curves for the different configurations. In Sec.IV, ^ we analyze the SGWBs from different origins. We sum- unit vectoru ^i. T (f; u;^ k) is the timing transfer function marize the conclusions in Sec.V. for each arm [67, 68]. The expectation of the cross-correlation between the channels I and J is defined as [69]: II. BASIC DEFINITIONS ∗ 0 1 0 hehI (f)eh (f )i = δ(f − f )ΓIJ (f)Sh(f); (8) J 2 A. Fundamentals of SGWB where Sh(f) is the one-sided power spectrum density (PSD) for SGWB, and the overlap reduction function The strength of SGWB is described by the ratio (ORF) [70] for I and J is: of SGWB energy density per logarithmic frequency Z 1 X P ^ P ∗ ^ −i2πfk^·∆~x=c bin to the critical density of the Universe ρc ≡ ΓIJ (f) = dΩk^FI (f; k)FJ (f; k)e ; 2 2 8π S2 3H0 c =(8πG)[57, 58]: P =+;× (9) 1 dρgw where ∆~x = ~xI − ~xJ , and we define ~xI;J as the location Ωgw(f) = : (1) 1 ρc d(ln f) where the interference happens. The ORF reflects the correlation of response between channels I and J, consid- Here dρgw is the SGWB energy density in the frequency ering a plane wave of frequency f and unit wave vector band [f, f + df], H0 is the Hubble constant, where c is k^. When applying to the same channel, i.e., I = J, the the light speed and G is the gravitational constant. In ORF reduces back to the transfer function. −1 −1 this work, we adopt H0 = 67:4 km s Mpc [59].