zhangxin@mail.neu.edu.cn
Prospects for weighing neutrinos in interacting dark energy models using joint observations of gravitational waves and -ray bursts
Abstract
Cosmological observations can be used to weigh neutrinos, but this method is model-dependent, with results relying on the cosmological model considered. If we consider interactions between dark energy and dark matter, the neutrino mass constraints differ from those derived under the standard model. On the contrary, gravitational wave (GW) standard siren observations can measure absolute cosmological distances, helping to break parameter degeneracies inherent in traditional cosmological observations, thereby improving constraints on neutrino mass. This paper examines the constraints on neutrino mass within interacting dark energy (IDE) models and explores how future GW standard siren observations could enhance these results. For multi-messenger GW observations, we consider the joint observations of binary neutron star mergers by third-generation ground-based GW detectors and short -ray burst observations by missions similar to the THESEUS satellite project. Using current cosmological observations (CMB+BAO+SN), we obtain an upper limit on the neutrino mass in the IDE models of 0.15 (or 0.16) eV. With the inclusion of GW data, the upper limit on the neutrino mass improves to 0.14 eV. This indicates that in the context of IDE models, the improvement in neutrino mass constraints from GW observations is relatively limited. However, GW observations significantly enhance the constraints on other cosmological parameters, such as matter density parameter, the Hubble constant, and coupling strength between dark energy and dark matter.
I Introduction
Experiments involving solar and atmospheric neutrinos have demonstrated that neutrinos possess mass and exhibit substantial mixing among different species [1]. Despite this, directly measuring the absolute mass scale of neutrinos remains a formidable challenge in particle physics experiments. Neutrino mass has implications for the cosmic microwave background (CMB) and the universeβs large-scale structure [2], which makes cosmological observations a crucial approach for determining their absolute mass.
Current principal cosmological observations fall into two categories. Firstly, there are observations related to the universeβs expansion history, such as baryon acoustic oscillations (BAO) and type Ia supernovae (SN). Secondly, observations concerning the cosmic structureβs growth, such as redshift space distortions, weak gravitational lensing, and galaxy clusters counts, are also pivotal. These observations help constrain the total neutrino mass, as extensively discussed in Refs. [3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 15, 16, 17, 18, 19, 20, 21, 23, 24, 25, 26, 27, 29, 30, 31, 33, 34, 37, 38, 39, 40, 43, 44, 22, 28, 35, 36, 47, 32, 45, 42, 46, 48, 49, 14, 41].
Recently, Ref. [44] explored the cosmological constraints on total neutrino mass within the framework of vacuum energy interacting with cold dark matter (ICDM), employing current observations. That study suggested that considering a direct interaction might alter the upper limit of the total neutrino mass, indicating that interactions between vacuum and cold dark matter could influence the cosmological measurement of neutrino mass. Nevertheless, current cosmological observations do not yet tightly constrain the total neutrino mass [50, 51]. However, gravitational-wave (GW) standard siren observations show promise as a powerful cosmological probe to aid in constraining the total neutrino mass.
The era of multi-messenger astronomy was inaugurated by the observation of the binary neutrino star (BNS) merger event GW170817 and its electromagnetic (EM) counterparts [52, 53]. Analyzing the GW waveform from such events directly yields the absolute luminosity distance, known as a standard siren. Additionally, identifying the EM counterpart of a GW source allows for the measurement of its redshift. Establishing the relationship between cosmic distances and redshifts paves the way for probing the universeβs expansion history. Thus, GW standard sirens are expected to significantly influence the estimation of cosmological parameters.
In the coming decades, third-generation (3G) ground-based GW detectors, such as Cosmic Explorer (CE) [54] and Einstein Telescope (ET) [55], are set to revolutionize the measurement of GW signals. These detectors are anticipated to form a network that will significantly enhance the capability to detect GW events. Ref. [56] presented a preliminary discussion on using future GW standard siren observations to weigh neutrinos within interacting dark energy (IDE) models. This analysis is based on an assumption that 1000 standard sirens could be detected over a 10-year period by either the ET or CE. While this estimate aligns with projections in other studies [29, 43, 44], variations in the redshift distribution could impact cosmological parameter estimations.
This work aims to explore the potential of 3G standard siren observations to measure neutrino mass in IDE models. We will analyze the detection strategy employing a network comprising the ET and two CEsβone in the US with a 40-km arm length and another in Australia with a 20-km arm length, collectively referred to as ET2CE. Additionally, we consider observations of -ray bursts (GRBs) from a detector similar to THESEUS, enhancing our capability to ascertain cosmological parameters.
In our analysis of IDE models, we focus exclusively on the ICDM models, characterized by a vacuum energy equation of state parameter . In these models, the energy conservation equations for vacuum energy and cold dark matter are described as follows:
(1) |
(2) |
where and represent the densities of dark energy and cold dark matter, respectively, is the Hubble parameter, the dot denotes the derivative with respect to cosmic time , and is the rate of energy transfer. In the field of IDE, one common assumption is that is proportional to the density of either dark matter or dark energy, i.e., or [57, 58, 59, 60, 61, 62, 63, 64, 65, 66], where is a dimensionless coupling parameter and facilitates computational ease. Another assumption posits as or [67, 68, 38, 40, 69, 70, 71], with the Hubble constant used for dimensional consistency. In this work, we consider two forms of the energy transfer rate: and . According to Eqs. (1) and (2), a positive implies that cold dark matter decays into vacuum energy, a negative suggests that vacuum energy decays into cold dark matter, and indicates no interaction between them.
II Methodology
The parameter space vector of the ICDM model is , where and are the present density of baryons and cold dark matter, respectively, is the Thomson scattering optical depth due to reionization, (multiplied by 100) is the radio between the comoving sound horizon and angular diameter distance at the decoupling epoch, is the coupling parameter in the ICDM model, is the amplitude of primordial scalar perturbation power spectrum, and is its power-law spectral index. When neutrino mass is considered in the ICDM model, one extra free parameter should be involved in the calculation.
Because we add neutrino mass into the ICDM model, the model considered in this paper is called the ICDM+ model. For convenience, in this paper, we use ICDM1+ and ICDM2+ to denote the corresponding and models that involve neutrino mass, respectively. Thus, there are eight independent parameters in total for these ICDM+ (ICDM1+ and ICDM2+) models.
It should be mentioned that, when we consider the situation of vacuum energy interacting with cold dark matter, then the vacuum energy is not a pure background any more, and in this case, we must consider the vacuum energy perturbations. Note that, in the calculation of the dark energy perturbation evolution, a problem of perturbation instability appears; namely, the cosmological perturbations of dark energy will be divergent in some parts of the parameter space, which ruins the IDE cosmology at the perturbation level. Under such a circumstance, to avoid the perturbation instability problem in the ICDM model, in this work, we treat the vacuum energy perturbations based on the extended parameterized post-Friedmann (PPF) approach [72, 73, 74] (for the original version of the PPF method, see Refs. [75, 76]). For more information about the calculation of cosmological perturbations in the IDE models and the PPF approach, see Refs. [77, 78, 65, 29].
In this paper, we simulate the GW standard siren data from the 3G GW detectors and future GRB detector and use them to constrain neutrino mass in the ICDM models. We will investigate whether the GW standard sirens can improve the constraint on neutrino mass.
To show the constraining capability of the simulated GW standard siren data, we consider two data combinations for comparison: CMB+BAO+SN (abbreviated as CBS) and CBS+GW. For the CMB data, we use the CMB likelihood, including the TT, TE, and EE spectra at , low- temperature commander likelihood, and low- SimAll EE likelihood, from the Planck 2018 data release [51]. For the BAO data, we use the measurements from the 6dFGS () [79], SDSS-MG () [80], and BOSS-DR12 ( 0.38, 0.51, and 0.61) [81]. For the SN data, we use the Pantheon sample, which is comprised of 1048 data points [82].
For the GW data, we use the simulation method described in Ref. [83]. Here, we provide only a brief overview. It is important to note that our approach includes a comprehensive calculation of the redshift distribution of GW-GRB events. This differs from the method used in previous studies, which assumed the detection of 1000 bright sirens over a 10-year observation period, as discussed in Refs. [91, 56, 84, 85, 86, 87, 88, 89, 90, 92, 93, 94, 95, 96].
Based on the star formation rate [97, 98, 99], the merger rate in the observer frame is
(3) |
where is the comoving volume element. is the BNS merger rate in the source frame that is given by
(4) |
where is the cosmic star formation rate, for which we adopt the Madau-Dickinson model of Ref. [100], is the universeβs age at the time of merger, is the time delay, and and are the minimum and maximum delay time for a massive binary to evolve merger [97], respectively. Here, Myr and is the Hubble time. The overall normalization is fixed by requiring that the value of agrees with the local rate estimated from the O1 LIGO and O2 LIGO/Virgo observation run [101], which is also in accordance with the latest O3 observation run [102]. is the distribution of the time delay , and we follow and adopt the exponential form [97]
(5) |
with an e-fold time of Myr for Myr. We simulate a catalog of BNS merger for 10-year observation. For each simulated GW source, the location , polarization angle , coalescence phase , and cosine of the inclination angle are drawn from uniform distributions. The mass of neutron star (NS) is assumed to be a Gaussian distribution, and the center value of the NS mass is , with a standard deviation of , where is the solar mass.
Under the stationary phase approximation (SPA) [103], the Fourier transform of the frequency-domain waveform for a detector network (with detectors) can be written as [104, 105]
(6) |
where is the diagonal matrix with , is the GW propagation direction, and is the location of the -th detector. () is
(7) |
where is the one-sided noise power spectral density of the -th detector, and is the frequency domain GW waveform.
We consider the waveform in the inspiralling stage for the non-spinning BNS system in this study, and we adopt the restricted Post-Newtonian (PN) approximation and calculate the waveform to the 3.5 PN order [106, 107]. The Fourier transform of the GW waveform of the -th detector can be expressed as
(8) |
where the detailed forms of and can be found in Refs. [106, 105]. is the Fourier amplitude, which can be written as
(9) |
where is the luminosity distance of the GW source, and are the antenna response functions of the -th GW detector, is the observed chirp mass, is the total mass of binary system with component masses and , and is the symmetric mass ratio. In this study, we adopt the GW waveform in the frequency domain, and then the time is replaced by [106, 105], where is the coalescence time.
To obtain the detection of GW events, we need to calculate the signal-to-noise ratio (SNR) for each GW event. The combined SNR for the detection network of detectors is given by
(10) |
The inner product is defined as
(11) |
where is the lower cutoff frequency ( Hz for ET and Hz for CE), is the frequency at the last stable orbit with , and are column matrices of the same dimension, and the star represents the conjugate transpose operator. Here, we adopt the SNR threshold to be 12 in the simulation, and we assume a running period of 10 years and duty cycle of 100% for the GW detector network.
For a GRB detected in coincidence with a GW signal, we require the peak flux to be higher than the flux limit of the satellite. Based on the fitting results of GRB170817A [108], we adopt the Gassian structured jet profile model
(12) |
where is the isotropically equivalent luminosity of short GRB observed at different viewing angle , and we assume the jetβs direction is aligned with the binary orbital angular momentum, namely, . is the on-axis isotropic luminosity defined by , and is the characteristic angle of the core. The detection probability of a short GRB is determined by . is the intrinsic luminosity function, given by
(13) |
where is the isotropic rest frame peak luminosity in the keV energy range, is a characteristic value separating the two regimes, and and are the characteristic slopes describing these regimes. In this paper, we use the values erg sec-1, , and , which is in accordance with Ref. [109]. In addition, we also assume a standard low-end cutoff in luminosity of erg sec-1. Here, we term the on-axis isotropic luminosity as the peak luminosity . According to the relation between luminosity and flux for GRB [110, 111], we can convert the flux limit to the luminosity [83], where is the 50-300 keV band for THESEUS-like telescope [112]. For the THESEUS-like telescope, we make the assumption of an 80% duty cycle. From the GW catalogue which has SNR larger than the threshold 12, then we can calculate the probability of the GRB detection for every GW event according to the probability distribution .
For the total uncertainty of the luminosity distance , we consider the instrument error , the weak-lensing error , and the peculiar velocity error , and thus the total error is given by
(14) |
We first estimate the instrument error by using the Fisher information matrix. The Fisher information matrix of a network including independent GW detectors can be written as
(15) |
where denotes nine parameters (, , , , , , , , ) for a GW event. The instrumental error of the parameter can be calculated by , where is the total Fisher information matrix for the network of interferometers.
The error caused by the peculiar velocity of the GW source is given by [115]
(17) |
where is the Hubble parameter, is the speed of light in a vacuum, and is the peculiar velocity of the GW source; we roughly use .
In this work, we wish to study the constraint on neutrino mass in the ICDM+ models from GW standard sirens. As the first attempt to explore the impacts of GW-GRB joint observation on the neutrino mass, we do not consider all the different cases of 3G GW observations, i.e., single ET, single CE, the CE-CE network (one in the US and another in Australia), and the ET-CE-CE network (one ET detector and two CE-like detectors). Instead, we only consider one case of ET-CE-CE network (abbreviated as GW hereafter) as a concrete example to analyze.
In the following, we use these observational data to place constraints on these ICDM+ models. We will use two data combinations, i.e., CBS and CBS+GW, to constrain the neutrino mass and other cosmological parameters. We employ the modified version of the publicly available Markov-Chain Monte Corlo (MCMC) package CosmoMC [116] to perform the calculations, and we use the PPF package [72, 73] for the ICDM models to handle the perturbations of vacuum energy.
III Results and discussion
In this section, we report the fitting results of these ICDM+ models and discuss the impacts of GW standard siren observation on constraining neutrino mass. We found that only approximately 640 standard sirens can be detected for the ET2CE+THESEUS network in 10 years, whose redshift distribution is given in Ref. [83]. The fitting results are presented in Tables 1 and 2 as well as Figs. 1 and 2. In Table 1, we quote errors for the parameters, but for the parameters that cannot be well-constrained, e.g., neutrino mass , we quote the 95.4% CL upper limits.
CDM+ β () | ICDM1+ β () | ICDM2+ β () | ||||||||
---|---|---|---|---|---|---|---|---|---|---|
Parameter | CBS | CBS+GW | CBS | CBS+GW | CBS | CBS+GW | ||||
β¦ | β¦ | |||||||||
CDM+ β () | ICDM1+ β () | ICDM2+ β () | ||||||||
---|---|---|---|---|---|---|---|---|---|---|
Parameter | CBS | CBS+GW | CBS | CBS+GW | CBS | CBS+GW | ||||
First, we discuss the effect of simulated GW data of ET2CE on constraining the total neutrino mass. In the case of CBS constraints, from Table 1, we have eV for CDM+, eV for ICDM1+, and eV for ICDM2+. After adding the GW data of ET2CE, namely, using the data combination CBS+GW, the constraint results become eV for CDM+, eV for ICDM1+, and eV for ICDM2+. We found that, compared with the CDM+ model, the constraint on neutrino mass becomes slighter looser in these ICDM+ models. We show the two-dimensional marginalized posterior contours (68.3% and 95.4% CL) in the β plane for the CDM+ and ICDM+ models in Fig. 1 (left panel for CBS and right panel for CBS+GW). From this figure, we can clearly see that once the interaction is considered in the model, the upper limit of becomes larger.
Moreover, we also found that the GW data help reduce the upper limits of neutrino mass by 16.26%, 4.67%, and 10.26% for the CDM+, ICDM1+, and ICDM2+ models, respectively. When the GW data are considered, the constraints on become slightly tighter in the CDM+ and ICDM+ models (also see Fig. 2). Thus, in the ICDM+ models, compared to the current observations, the GW standard siren observations from ET2CE can slightly improve the constraint on neutrino mass, which is in accordance with the conclusion on the IDE models in the previous study [56].
Next, we discuss how the GW data help improve the constraints on other parameters, i.e., , , and . In Table 2, we show the absolute and relative errors of , , and from the CBS and CBS+GW data combinations. Comparing the results from two data combinations, we found that the constraints on , , and are improved by 73.53%β83.95%, 5.56%β30%, and 90%β92.68%, respectively, when considering the ET2CE data. Obviously, the GW data can indeed effectively improve the constraints on the parameters and (also see Figs. 1 and 2).
Finally, we present the constraint results of the coupling constant . By using CBS data, we have for the ICDM1+ model and for ICDM2+ model. It is shown that a positive value of is slightly favored and is only at and levels, respectively. By using CBS+GW data, we have for the ICDM1+ model and for the ICDM2+ model. Thus, we found that the error is slightly shrunk and now is favored at and levels, respectively. For these two ICDM+ models, one can clearly see that, in the ICDM+ cases (with and ), is favored within the significance range, no matter whether the GW data are taken into account. Thus, in the two ICDM+ models, there is no evidence of a nonzero interaction. Moreover, comparing the values of from the two data combinations, we also found that the accuracy of is increased by 23.08% () and 42.11% () when adding the GW data of ET2CE in the cosmological fit (also see Fig. 2). This indicates that the GW data of the ET2CE can also improve the constraint accuracies of coupling parameter .
IV Conclusion
In this paper, we investigated how GW standard sirens from 3G ground-based GW detectors can constrain the total neutrino mass () in interacting vacuum energy models with energy transfer forms and . We used an extended version of the PPF approach to compute the vacuum energy perturbations within these interacting scenarios. We focused on the synergy between 3G GW detectors and a GRB detector for multi-messenger observations, specifically employing the ET2CE network as a representative for GW discussions. To evaluate the impact of GW data on the constraints of , we also considered existing CMB+BAO+SN data for comparison and combination.
Our analysis revealed that GW data can reduce the upper limits of by 16.26%, 4.67%, and 10.26% for the CDM+, ICDM1+, and ICDM2+ models, respectively. Thus, GW standard siren data from ET2CE offer a modest improvement in constraining compared to CBS alone. For the derived parameters and , incorporating GW data into the cosmological fit substantially enhances the precision: the accuracy of improved by 73.53%β83.95% and that of by 90%β92.68%. These significant enhancements underscore the value of including GW data from ET2CE. Moreover, based on the combined constraints from CBS and CBS+GW data, we found that GW data from ET2CE also significantly refine the accuracy of the coupling strength .
Acknowledgements.
We thank Shang-Jie Jin for helpful discussions. This work was supported by the National Natural Science Foundation of China (12305069, 11947022, 11975072, 11875102, 11835009), the National SKA Program of China (2022SKA0110200, 2022SKA0110203), the National 111 Project (B16009), and the Program of the Education Department of Liaoning Province (JYTMS20231695).References
- [1] J. Lesgourgues and S. Pastor, Phys. Rept. 429 (2006), 307-379 [arXiv:astro-ph/0603494 [astro-ph]].
- [2] K. N. Abazajian et al. [Topical Conveners: K.N. Abazajian, J.E. Carlstrom, A.T. Lee], Astropart. Phys. 63 (2015), 66-80 [arXiv:1309.5383 [astro-ph.CO]].
- [3] W. Hu, D. J. Eisenstein and M. Tegmark, Phys. Rev. Lett. 80 (1998), 5255-5258 [arXiv:astro-ph/9712057 [astro-ph]].
- [4] B. A. Reid, L. Verde, R. Jimenez and O. Mena, JCAP 01 (2010), 003 [arXiv:0910.0008 [astro-ph.CO]].
- [5] H. Li and X. Zhang, Phys. Lett. B 713 (2012), 160-164 [arXiv:1202.4071 [astro-ph.CO]].
- [6] Y. H. Li, S. Wang, X. D. Li and X. Zhang, JCAP 02 (2013), 033 [arXiv:1207.6679 [astro-ph.CO]].
- [7] X. Wang, X. L. Meng, T. J. Zhang, H. Shan, Y. Gong, C. Tao, X. Chen and Y. F. Huang, JCAP 11 (2012), 018 [arXiv:1210.2136 [astro-ph.CO]].
- [8] B. Audren, J. Lesgourgues, S. Bird, M. G. Haehnelt and M. Viel, JCAP 01 (2013), 026 [arXiv:1210.2194 [astro-ph.CO]].
- [9] S. Riemer-SΓΈrensen, D. Parkinson and T. M. Davis, Phys. Rev. D 89 (2014), 103505 [arXiv:1306.4153 [astro-ph.CO]].
- [10] J. F. Zhang, Y. H. Li and X. Zhang, Eur. Phys. J. C 74 (2014), 2954 [arXiv:1404.3598 [astro-ph.CO]].
- [11] X. Y. Zhou and J. H. He, Commun. Theor. Phys. 62 (2014), 102-108 [arXiv:1406.6822 [astro-ph.CO]].
- [12] J. F. Zhang, J. J. Geng and X. Zhang, JCAP 10 (2014), 044 [arXiv:1408.0481 [astro-ph.CO]].
- [13] N. Palanque-Delabrouille, C. Yèche, J. Lesgourgues, G. Rossi, A. Borde, M. Viel, E. Aubourg, D. Kirkby, J. M. LeGoff and J. Rich, et al. JCAP 02 (2015), 045 [arXiv:1410.7244 [astro-ph.CO]].
- [14] C. Q. Geng, C. C. Lee and J. L. Shen, Phys. Lett. B 740 (2015), 285-290 [arXiv:1411.3813 [astro-ph.CO]].
- [15] J. F. Zhang, M. M. Zhao, Y. H. Li and X. Zhang, JCAP 04 (2015), 038 [arXiv:1502.04028 [astro-ph.CO]].
- [16] C. Q. Geng, C. C. Lee, R. Myrzakulov, M. Sami and E. N. Saridakis, JCAP 01 (2016), 049 [arXiv:1504.08141 [astro-ph.CO]].
- [17] Y. Chen and L. Xu, Phys. Lett. B 752 (2016), 66-75 [arXiv:1507.02008 [astro-ph.CO]].
- [18] X. Zhang, Phys. Rev. D 93 (2016) no.8, 083011 [arXiv:1511.02651 [astro-ph.CO]].
- [19] Q. G. Huang, K. Wang and S. Wang, Eur. Phys. J. C 76 (2016) no.9, 489 [arXiv:1512.05899 [astro-ph.CO]].
- [20] Y. Chen, B. Ratra, M. Biesiada, S. Li and Z. H. Zhu, Astrophys. J. 829 (2016) no.2, 61 [arXiv:1603.07115 [astro-ph.CO]].
- [21] M. Moresco, R. Jimenez, L. Verde, A. Cimatti, L. Pozzetti, C. Maraston and D. Thomas, JCAP 12 (2016), 039 [arXiv:1604.00183 [astro-ph.CO]].
- [22] E. Giusarma, M. Gerbino, O. Mena, S. Vagnozzi, S. Ho and K. Freese, Phys. Rev. D 94 (2016) no.8, 083522 [arXiv:1605.04320 [astro-ph.CO]].
- [23] J. Lu, M. Liu, Y. Wu, Y. Wang and W. Yang, Eur. Phys. J. C 76 (2016) no.12, 679 [arXiv:1606.02987 [astro-ph.CO]].
- [24] S. Wang, Y. F. Wang, D. M. Xia and X. Zhang, Phys. Rev. D 94 (2016) no.8, 083519 [arXiv:1608.00672 [astro-ph.CO]].
- [25] M. M. Zhao, Y. H. Li, J. F. Zhang and X. Zhang, Mon. Not. Roy. Astron. Soc. 469 (2017) no.2, 1713-1724 [arXiv:1608.01219 [astro-ph.CO]].
- [26] S. Kumar and R. C. Nunes, Phys. Rev. D 94 (2016) no.12, 123511 [arXiv:1608.02454 [astro-ph.CO]].
- [27] L. Xu and Q. G. Huang, Sci. China Phys. Mech. Astron. 61 (2018) no.3, 039521 [arXiv:1611.05178 [astro-ph.CO]].
- [28] S. Vagnozzi, E. Giusarma, O. Mena, K. Freese, M. Gerbino, S. Ho and M. Lattanzi, Phys. Rev. D 96 (2017) no.12, 123503 [arXiv:1701.08172 [astro-ph.CO]].
- [29] R. Y. Guo, Y. H. Li, J. F. Zhang and X. Zhang, JCAP 05 (2017), 040 [arXiv:1702.04189 [astro-ph.CO]].
- [30] X. Zhang, Sci. China Phys. Mech. Astron. 60 (2017) no.6, 060431 [arXiv:1703.00651 [astro-ph.CO]].
- [31] E. K. Li, H. Zhang, M. Du, Z. H. Zhou and L. Xu, JCAP 08 (2018), 042 [arXiv:1703.01554 [astro-ph.CO]].
- [32] W. Yang, R. C. Nunes, S. Pan and D. F. Mota, Phys. Rev. D 95 (2017) no.10, 103522 [arXiv:1703.02556 [astro-ph.CO]].
- [33] S. Wang, Y. F. Wang and D. M. Xia, Chin. Phys. C 42 (2018) no.6, 065103 [arXiv:1707.00588 [astro-ph.CO]].
- [34] M. M. Zhao, J. F. Zhang and X. Zhang, Phys. Lett. B 779 (2018), 473-478 [arXiv:1710.02391 [astro-ph.CO]].
- [35] S. Vagnozzi, S. Dhawan, M. Gerbino, K. Freese, A. Goobar and O. Mena, Phys. Rev. D 98 (2018) no.8, 083501 [arXiv:1801.08553 [astro-ph.CO]].
- [36] E. Giusarma, S. Vagnozzi, S. Ho, S. Ferraro, K. Freese, R. Kamen-Rubio and K. B. Luk, Phys. Rev. D 98 (2018) no.12, 123526 [arXiv:1802.08694 [astro-ph.CO]].
- [37] R. Y. Guo, J. F. Zhang and X. Zhang, Chin. Phys. C 42 (2018) no.9, 095103 [arXiv:1803.06910 [astro-ph.CO]].
- [38] L. Feng, H. L. Li, J. F. Zhang and X. Zhang, Sci. China Phys. Mech. Astron. 63 (2020) no.2, 220401 [arXiv:1903.08848 [astro-ph.CO]].
- [39] J. F. Zhang, B. Wang and X. Zhang, Sci. China Phys. Mech. Astron. 63 (2020) no.8, 280411 [arXiv:1907.00179 [astro-ph.CO]].
- [40] L. Feng, D. Z. He, H. L. Li, J. F. Zhang and X. Zhang, Sci. China Phys. Mech. Astron. 63 (2020) no.9, 290404 [arXiv:1910.03872 [astro-ph.CO]].
- [41] Z. Liu and H. Miao, Int. J. Mod. Phys. D 29 (2020) no.13, 2050088 [arXiv:2002.05563 [astro-ph.CO]].
- [42] W. Yang, E. Di Valentino, O. Mena and S. Pan, Phys. Rev. D 102 (2020) no.2, 023535 [arXiv:2003.12552 [astro-ph.CO]].
- [43] M. Zhang, J. F. Zhang and X. Zhang, Commun. Theor. Phys. 72 (2020) no.12, 125402 [arXiv:2005.04647 [astro-ph.CO]].
- [44] H. L. Li, J. F. Zhang and X. Zhang, Commun. Theor. Phys. 72 (2020) no.12, 125401 [arXiv:2005.12041 [astro-ph.CO]].
- [45] W. Yang, E. Di Valentino, S. Pan and O. Mena, Phys. Dark Univ. 31 (2021), 100762 [arXiv:2007.02927 [astro-ph.CO]].
- [46] H. R. Amiri, A. Salehi and A. H. Noroozi, Eur. Phys. J. C 81 (2021) no.5, 479
- [47] I. Tanseri, S. Hagstotz, S. Vagnozzi, E. Giusarma and K. Freese, JHEAp 36 (2022), 1-26 [arXiv:2207.01913 [astro-ph.CO]].
- [48] P. Bazvand, A. Salehi and R. Sepahvand, Eur. Phys. J. Plus 138 (2023) no.5, 448
- [49] Y. H. Pang, X. Zhang and Q. G. Huang, Chin. Phys. C 48 (2024) no.6, 065102 [arXiv:2312.07188 [astro-ph.CO]].
- [50] N. Aghanim et al. [Planck], Astron. Astrophys. 594 (2016), A11 [arXiv:1507.02704 [astro-ph.CO]].
- [51] N. Aghanim et al. [Planck], Astron. Astrophys. 641 (2020), A6 [arXiv:1807.06209 [astro-ph.CO]].
- [52] B. P. Abbott et al. [LIGO Scientific and Virgo], Phys. Rev. Lett. 119 (2017) no.16, 161101 [arXiv:1710.05832 [gr-qc]].
- [53] B. P. Abbott et al. [LIGO Scientific, Virgo, Fermi GBM, INTEGRAL, IceCube, AstroSat Cadmium Zinc Telluride Imager Team, IPN, Insight-Hxmt, ANTARES, Swift, AGILE Team, 1M2H Team, Dark Energy Camera GW-EM, DES, DLT40, GRAWITA, Fermi-LAT, ATCA, ASKAP, Las Cumbres Observatory Group, OzGrav, DWF (Deeper Wider Faster Program), AST3, CAASTRO, VINROUGE, MASTER, J-GEM, GROWTH, JAGWAR, CaltechNRAO, TTU-NRAO, NuSTAR, Pan-STARRS, MAXI Team, TZAC Consortium, KU, Nordic Optical Telescope, ePESSTO, GROND, Texas Tech University, SALT Group, TOROS, BOOTES, MWA, CALET, IKI-GW Follow-up, H.E.S.S., LOFAR, LWA, HAWC, Pierre Auger, ALMA, Euro VLBI Team, Pi of Sky, Chandra Team at McGill University, DFN, ATLAS Telescopes, High Time Resolution Universe Survey, RIMAS, RATIR and SKA South Africa/MeerKAT], Astrophys. J. Lett. 848 (2017) no.2, L12 [arXiv:1710.05833 [astro-ph.HE]].
- [54] B. P. Abbott et al. [LIGO Scientific], Class. Quant. Grav. 34 (2017) no.4, 044001 [arXiv:1607.08697 [astro-ph.IM]].
- [55] M. Punturo, M. Abernathy, F. Acernese, B. Allen, N. Andersson, K. Arun, F. Barone, B. Barr, M. Barsuglia and M. Beker, et al. Class. Quant. Grav. 27 (2010), 194002 doi:10.1088/0264-9381/27/19/194002
- [56] S. J. Jin, R. Q. Zhu, L. F. Wang, H. L. Li, J. F. Zhang and X. Zhang, Commun. Theor. Phys. 74 (2022) no.10, 105404 [arXiv:2204.04689 [astro-ph.CO]].
- [57] M. Li, X. D. Li, S. Wang, Y. Wang and X. Zhang, JCAP 12 (2009), 014 [arXiv:0910.3855 [astro-ph.CO]].
- [58] Y. Li, J. Ma, J. Cui, Z. Wang and X. Zhang, Sci. China Phys. Mech. Astron. 54 (2011), 1367-1377 [arXiv:1011.6122 [astro-ph.CO]].
- [59] T. F. Fu, J. F. Zhang, J. Q. Chen and X. Zhang, Eur. Phys. J. C 72 (2012), 1932 [arXiv:1112.2350 [astro-ph.CO]].
- [60] Z. Zhang, S. Li, X. D. Li, X. Zhang and M. Li, JCAP 06 (2012), 009 [arXiv:1204.6135 [astro-ph.CO]].
- [61] S. Wang, Y. Z. Wang, J. J. Geng and X. Zhang, Eur. Phys. J. C 74 (2014) no.11, 3148 [arXiv:1406.0072 [astro-ph.CO]].
- [62] J. L. Cui, L. Yin, L. F. Wang, Y. H. Li and X. Zhang, JCAP 09 (2015), 024 [arXiv:1503.08948 [astro-ph.CO]].
- [63] L. Feng and X. Zhang, JCAP 08 (2016), 072 [arXiv:1607.05567 [astro-ph.CO]].
- [64] L. Feng, J. F. Zhang and X. Zhang, Phys. Dark Univ. 23 (2019), 100261 [arXiv:1712.03148 [astro-ph.CO]].
- [65] L. Feng, Y. H. Li, F. Yu, J. F. Zhang and X. Zhang, Eur. Phys. J. C 78 (2018) no.10, 865 [arXiv:1807.03022 [astro-ph.CO]].
- [66] H. L. Li, L. Feng, J. F. Zhang and X. Zhang, Sci. China Phys. Mech. Astron. 62 (2019) no.12, 120411 [arXiv:1812.00319 [astro-ph.CO]].
- [67] J. Zhang, L. Zhao and X. Zhang, Sci. China Phys. Mech. Astron. 57 (2014), 387-392 [arXiv:1306.1289 [astro-ph.CO]].
- [68] H. L. Li, J. F. Zhang, L. Feng and X. Zhang, Eur. Phys. J. C 77 (2017) no.12, 907 [arXiv:1711.06159 [astro-ph.CO]].
- [69] L. F. Wang, J. H. Zhang, D. Z. He, J. F. Zhang and X. Zhang, Mon. Not. Roy. Astron. Soc. 514 (2022) no.1, 1433-1440 [arXiv:2102.09331 [astro-ph.CO]].
- [70] Z. W. Zhao, L. F. Wang, J. G. Zhang, J. F. Zhang and X. Zhang, JCAP 04 (2023), 022 [arXiv:2210.07162 [astro-ph.CO]].
- [71] T. N. Li, S. J. Jin, H. L. Li, J. F. Zhang and X. Zhang, Astrophys. J. 963 (2024) no.1, 52 [arXiv:2310.15879 [astro-ph.CO]].
- [72] Y. H. Li, J. F. Zhang and X. Zhang, Phys. Rev. D 90 (2014) no.6, 063005 [arXiv:1404.5220 [astro-ph.CO]].
- [73] Y. H. Li, J. F. Zhang and X. Zhang, Phys. Rev. D 90 (2014) no.12, 123007 [arXiv:1409.7205 [astro-ph.CO]].
- [74] Y. H. Li and X. Zhang, JCAP 09 (2023), 046 [arXiv:2306.01593 [astro-ph.CO]].
- [75] W. Hu, Phys. Rev. D 77 (2008), 103524 [arXiv:0801.2433 [astro-ph]].
- [76] W. Fang, W. Hu and A. Lewis, Phys. Rev. D 78 (2008), 087303 [arXiv:0808.3125 [astro-ph]].
- [77] Y. H. Li, J. F. Zhang and X. Zhang, Phys. Rev. D 93 (2016) no.2, 023002 [arXiv:1506.06349 [astro-ph.CO]].
- [78] X. Zhang, Sci. China Phys. Mech. Astron. 60 (2017) no.5, 050431 [arXiv:1702.04564 [astro-ph.CO]].
- [79] F. Beutler, C. Blake, M. Colless, D. H. Jones, L. Staveley-Smith, L. Campbell, Q. Parker, W. Saunders and F. Watson, Mon. Not. Roy. Astron. Soc. 416 (2011), 3017-3032 [arXiv:1106.3366 [astro-ph.CO]].
- [80] A. J. Ross, L. Samushia, C. Howlett, W. J. Percival, A. Burden and M. Manera, Mon. Not. Roy. Astron. Soc. 449 (2015) no.1, 835-847 [arXiv:1409.3242 [astro-ph.CO]].
- [81] S. Alam et al. [BOSS], Mon. Not. Roy. Astron. Soc. 470 (2017) no.3, 2617-2652 [arXiv:1607.03155 [astro-ph.CO]].
- [82] D. M. Scolnic, D. O. Jones, A. Rest, Y. C. Pan, R. Chornock, R. J. Foley, M. E. Huber, R. Kessler, G. Narayan and A. G. Riess, et al. Astrophys. J. 859 (2018) no.2, 101 [arXiv:1710.00845 [astro-ph.CO]].
- [83] T. Han, S. J. Jin, J. F. Zhang and X. Zhang, Eur. Phys. J. C 84, no.7, 663 (2024) [arXiv:2309.14965 [astro-ph.CO]].
- [84] R. G. Cai and T. Yang, Phys. Rev. D 95 (2017) no.4, 044024 [arXiv:1608.08008 [astro-ph.CO]].
- [85] R. G. Cai, T. B. Liu, X. W. Liu, S. J. Wang and T. Yang, Phys. Rev. D 97 (2018) no.10, 103005 [arXiv:1712.00952 [astro-ph.CO]].
- [86] L. F. Wang, X. N. Zhang, J. F. Zhang and X. Zhang, Phys. Lett. B 782 (2018), 87-93 [arXiv:1802.04720 [astro-ph.CO]].
- [87] X. N. Zhang, L. F. Wang, J. F. Zhang and X. Zhang, Phys. Rev. D 99 (2019) no.6, 063510 [arXiv:1804.08379 [astro-ph.CO]].
- [88] X. Zhang, Sci. China Phys. Mech. Astron. 62 (2019) no.11, 110431 [arXiv:1905.11122 [astro-ph.CO]].
- [89] J. F. Zhang, H. Y. Dong, J. Z. Qi and X. Zhang, Eur. Phys. J. C 80 (2020) no.3, 217 [arXiv:1906.07504 [astro-ph.CO]].
- [90] J. F. Zhang, M. Zhang, S. J. Jin, J. Z. Qi and X. Zhang, JCAP 09 (2019), 068 [arXiv:1907.03238 [astro-ph.CO]].
- [91] H. L. Li, D. Z. He, J. F. Zhang and X. Zhang, JCAP 06 (2020), 038 [arXiv:1908.03098 [astro-ph.CO]].
- [92] S. J. Jin, D. Z. He, Y. Xu, J. F. Zhang and X. Zhang, JCAP 03 (2020), 051 [arXiv:2001.05393 [astro-ph.CO]].
- [93] S. J. Jin, L. F. Wang, P. J. Wu, J. F. Zhang and X. Zhang, Phys. Rev. D 104 (2021) no.10, 103507 [arXiv:2106.01859 [astro-ph.CO]].
- [94] P. J. Wu, Y. Shao, S. J. Jin and X. Zhang, JCAP 06 (2023), 052 [arXiv:2202.09726 [astro-ph.CO]].
- [95] S. J. Jin, S. S. Xing, Y. Shao, J. F. Zhang and X. Zhang, Chin. Phys. C 47 (2023) no.6, 065104 [arXiv:2301.06722 [astro-ph.CO]].
- [96] J. G. Zhang, Z. W. Zhao, Y. Li, J. F. Zhang, D. Li and X. Zhang, Sci. China Phys. Mech. Astron. 66 (2023) no.12, 120412 [arXiv:2307.01605 [astro-ph.CO]].
- [97] S. Vitale, W. M. Farr, K. Ng and C. L. Rodriguez, Astrophys. J. Lett. 886 (2019) no.1, L1 [arXiv:1808.00901 [astro-ph.HE]].
- [98] E. Belgacem, Y. Dirian, S. Foffa, E. J. Howell, M. Maggiore and T. Regimbau, JCAP 08 (2019), 015 [arXiv:1907.01487 [astro-ph.CO]].
- [99] T. Yang, JCAP 05 (2021), 044 [arXiv:2103.01923 [astro-ph.CO]].
- [100] P. Madau and M. Dickinson, Ann. Rev. Astron. Astrophys. 52 (2014), 415-486 [arXiv:1403.0007 [astro-ph.CO]].
- [101] A. Eichhorn, T. Koslowski and A. D. Pereira, Universe 5 (2019) no.2, 53 [arXiv:1811.12909 [gr-qc]].
- [102] R. Abbott et al. [KAGRA, VIRGO and LIGO Scientific], Phys. Rev. X 13 (2023) no.1, 011048 [arXiv:2111.03634 [astro-ph.HE]].
- [103] X. Zhang, T. Liu and W. Zhao, Phys. Rev. D 95 (2017) no.10, 104027 [arXiv:1702.08752 [gr-qc]].
- [104] L. Wen and Y. Chen, Phys. Rev. D 81 (2010), 082001 [arXiv:1003.2504 [astro-ph.CO]].
- [105] W. Zhao and L. Wen, Phys. Rev. D 97 (2018) no.6, 064031 [arXiv:1710.05325 [astro-ph.CO]].
- [106] C. Cutler, T. A. Apostolatos, L. Bildsten, L. S. Finn, E. E. Flanagan, D. Kennefick, D. M. Markovic, A. Ori, E. Poisson and G. J. Sussman, et al. Phys. Rev. Lett. 70 (1993), 2984-2987 [arXiv:astro-ph/9208005 [astro-ph]].
- [107] B. S. Sathyaprakash and B. F. Schutz, Living Rev. Rel. 12 (2009), 2 [arXiv:0903.0338 [gr-qc]].
- [108] E. J. Howell, K. Ackley, A. Rowlinson and D. Coward, [arXiv:1811.09168 [astro-ph.HE]].
- [109] D. Wanderman and T. Piran, Mon. Not. Roy. Astron. Soc. 448 (2015) no.4, 3026-3037 [arXiv:1405.5878 [astro-ph.HE]].
- [110] P. Meszaros and A. Meszaros, Astrophys. J. 449 (1995), 9-16 [arXiv:astro-ph/9503087 [astro-ph]].
- [111] A. Meszaros, J. Ripa and F. Ryde, Astron. Astrophys. 529 (2011), A55 [arXiv:1101.5040 [astro-ph.CO]].
- [112] G. Stratta, L. Amati, R. Ciolfi and S. Vinciguerra, Mem. Soc. Ast. It. 89 (2018) no.2, 205-212 [arXiv:1802.01677 [astro-ph.IM]].
- [113] C. M. Hirata, D. E. Holz and C. Cutler, Phys. Rev. D 81 (2010), 124046 [arXiv:1004.3988 [astro-ph.CO]].
- [114] L. Speri, N. Tamanini, R. R. Caldwell, J. R. Gair and B. Wang, Phys. Rev. D 103 (2021) no.8, 083526 [arXiv:2010.09049 [astro-ph.CO]].
- [115] B. Kocsis, Z. Frei, Z. Haiman and K. Menou, Astrophys. J. 637 (2006), 27-37 [arXiv:astro-ph/0505394 [astro-ph]].
- [116] A. Lewis and S. Bridle, Phys. Rev. D 66 (2002), 103511 [arXiv:astro-ph/0205436 [astro-ph]].