Welcome to visit Communications in Theoretical Physics,
Gravitation Theory, Astrophysics and Cosmology

Prospects for searching for sterile neutrinos with gravitational wave and γ-ray burst joint observations

  • Lu Feng 1, 2, 5 ,
  • Tao Han 2, 5 ,
  • Jing-Fei Zhang 2 ,
  • Xin Zhang , 2, 3, 4,
Expand
  • 1College of Physical Science and Technology, Shenyang Normal University, Shenyang 110034, China
  • 2Liaoning Key Laboratory of Cosmology and Astrophysics, College of Sciences, Northeastern University, Shenyang 110819, China
  • 3MOE Key Laboratory of Data Analytics and Optimization for Smart Industry, Northeastern University, Shenyang 110819, China
  • 4National Frontiers Science Center for Industrial Intelligence and Systems Optimization, Northeastern University, Shenyang 110819, China

These authors contributed equally to this paper.

Author to whom any correspondence should be addressed.

Received date: 2024-09-05

  Revised date: 2024-11-20

  Accepted date: 2024-11-26

  Online published: 2025-02-27

Copyright

© 2025 Institute of Theoretical Physics CAS, Chinese Physical Society and IOP Publishing. All rights, including for text and data mining, AI training, and similar technologies, are reserved.
This article is available under the terms of the IOP-Standard License.

Cite this article

Lu Feng , Tao Han , Jing-Fei Zhang , Xin Zhang . Prospects for searching for sterile neutrinos with gravitational wave and γ-ray burst joint observations[J]. Communications in Theoretical Physics, 2025 , 77(6) : 065403 . DOI: 10.1088/1572-9494/ad9c3e

1. Introduction

On August 17th 2017, the first observation of gravitational waves (GW) from a binary neutron star (BNS) merger [1], together with the first joint observation of GW from a BNS merger and its electromagnetic (EM) counterpart [2, 3], marked the beginning of a new era in multi-messenger astronomy and cosmology. The measurement of the GW signal directly provides information on the absolute luminosity distance to the source, while its redshift can be determined by identifying the EM counterpart of the GW source. This establishes an absolute distance-redshift relation, known as the standard siren method, which is crucial for cosmological studies. To date, only one bright siren, GW170817, has been identified. This is insufficient to probe cosmological parameters using current standard sirens, necessitating the use of next-generation GW detectors.
In the future, third-generation (3G) ground-based GW detectors, such as the Einstein Telescope (ET) [4, 5] in Europe and the Cosmic Explorer (CE) [6, 7] in the United States will become operational, with sensitivities improved one order of magnitude over the current detectors, and much more BNS merger events will be observed at much deeper redshifts. Recently, GW standard sirens have been widely discussed in the literature [847]. It has been found that future observations of GW standard sirens from the ET and CE will play a crucial role in the estimation of cosmological parameters [14, 16, 22, 25, 26, 28, 37, 38, 43, 46]. In particular, the GW standard sirens can break the parameter degeneracies generated by the current EM cosmological observations, thereby improving constraints on neutrino mass, see [14, 37].
A recent forecast [46] demonstrated that joint observations of BNS by 3G GW detectors and short γ-ray burst (GRB) observations by missions similar to the THESEUS satellite project can improve the constraints on the total active neutrino mass. Therefore, it is crucial to investigate how the combined GW-GRB observations would affect constraints on the sterile neutrino parameters.
The existence of light sterile neutrinos has been suggested by anomalies in short-baseline (SBL) neutrino experiments [4857]. To explain the SBL neutrino oscillation data, sterile neutrinos with eV-scale masses are required [5860]. Cosmological observations play a crucial role in constraining the mass of active neutrinos (see [61108]). Since sterile neutrinos have implications for the evolution of the Universe, cosmology can provide an independent test for their existence. For related works on sterile neutrinos, see [109134].
Currently, one of the most significant challenges in cosmology is the 'Hubble tension' [135], which refers to the discrepancy between early and late Universe observations. In the past few years, people often considered models that include light sterile neutrinos to alleviate the Hubble constant crisis. This is because when using the cosmic microwave background data to constrain cosmological parameters, the effective number of neutrino species (Neff) is positively correlated with the Hubble constant; if Neff is larger, the derived Hubble constant will also be larger. See [118, 125128, 130, 131, 134], for related studies. However, in recent years, the results of cosmological observation fits have shown that the effect of using this method to alleviate the Hubble crisis is no longer significant. Nevertheless, due to the significant correlation between Neff and H0, precise measurements of the Hubble constant using gravitational wave standard sirens are very helpful for determining the parameters of sterile neutrinos.
Additionally, the main advantage of the standard siren method for measuring the Hubble constant is that it avoids relying on the cosmic distance ladder. Therefore, in the future, the GW standard sirens could become a promising cosmological probe, playing a crucial role in measuring cosmological parameters, including those related to sterile neutrinos.
In this paper, we present a forecast for the search for sterile neutrinos using joint GW-GRB observations. The primary aim of this work is to investigate the impact of future GW standard siren observations on the constraints of sterile neutrino parameters.
This work is organized as follows. In section 2, we introduce the methodology used in this work. In section 3, we give the constraint results and make some relevant discussions. The conclusion is given in section 4.

2. Methodology

2.1. Gravitational wave simulation

In this subsection, we introduce the method of simulating the joint GW standard sirens and GRB events. We consider the THESEUS-like GRB detector in synergy with the 3G GW observation. We use the simulation method as prescribed in [43, 46]. Here, we provide only a brief overview.
The BNS merger rate with redshift in the observer frame is [136138]
$\begin{eqnarray}{R}_{{\rm{m}}}(z)=\frac{{{ \mathcal R }}_{{\rm{m}}}(z)}{1+z}\frac{{\rm{d}}V(z)}{{\rm{d}}z},\end{eqnarray}$
where dV(z)/dz is the comoving volume element, the factor (1 + z)−1 converts the merger rate in the source frame to the observer frame, and ${{ \mathcal R }}_{{\rm{m}}}(z)$ is the BNS merger rate in the source frame, expressed as
$\begin{eqnarray}{{ \mathcal R }}_{{\rm{m}}}(z)={\int }_{{t}_{{\rm{\min }}}}^{{t}_{{\rm{\max }}}}{{ \mathcal R }}_{{\rm{f}}}[t(z)-{t}_{{\rm{d}}}]P({t}_{{\rm{d}}}){\rm{d}}{t}_{{\rm{d}}},\end{eqnarray}$
which is commonly used in the literature [137147]6. Here, td is the delay time between the formation of BNS system and merger, ${t}_{{\rm{\min }}}=20$ Myr is the minimum delay time, ${t}_{{\rm{\max }}}={t}_{{\rm{H}}}$ is the maximum delay time, t(z) is the age of the Universe at the time of merger, ${{ \mathcal R }}_{{\rm{f}}}$ is the cosmic star formation rate in the source frame for which we adopt the Madau–Dickinson model [149], P(td) is the time delay distribution of the td, and we adopt the exponential time delay model [136], which is given by
$\begin{eqnarray}P({t}_{{\rm{d}}})=\frac{1}{\tau }{\rm{\exp }}(-{t}_{{\rm{d}}}/\tau ),\end{eqnarray}$
with an e-fold time of τ = 0.1 Gyr for ${t}_{{\rm{d}}}\gt {t}_{{\rm{\min }}}$.
In our calculations, for BNS mergers, we consider the local comoving merger rate to be ${{ \mathcal R }}_{{\rm{m}}}(z=0)=920\,{\rm{G}}p{c}^{-3}\,y{r}^{-1}$, which is the estimated median from the O1 LIGO and the O2 LIGO/Virgo observation run [150] and is also consistent with the O3 observation run [151]. We simulate a catalog of BNS mergers for 10 years observation. For each source, the location (θ, φ), the polarization angle ψ, the cosine of the inclination angle ι, and the coalescence phase ψc are drawn from uniform distributions. Currently, there are multiple candidate models for the neutron star (NS) mass distribution. However, different mass distributions of NSs have less impact on the cosmological analysis [43]. For simplicity, we employ a Gaussian mass distribution. This distribution has a mean of 1.33M for the NS mass and a standard deviation of 0.09M, where M represents the solar mass [152, 153].
Under the stationary phase approximation [154], the Fourier transform of the frequency-domain GW waveform for a detector network (with N detectors) is given by [155, 156]
$\begin{eqnarray}\tilde{{\boldsymbol{h}}}(f)={{\rm{e}}}^{{\rm{-i}}{\boldsymbol{\Phi }}}{\boldsymbol{h}}(f),\end{eqnarray}$
with the h(f) is given by
$\begin{eqnarray}{\boldsymbol{h}}(f)={\left[\frac{{h}_{1}(f)}{\sqrt{{S}_{{\rm{n}},1}(f)}},\frac{{h}_{2}(f)}{\sqrt{{S}_{{\rm{n}},2}(f)}},\ldots ,\frac{{h}_{N}(f)}{\sqrt{{S}_{{\rm{n}},N}(f)}}\right]}^{{\rm{T}}},\end{eqnarray}$
where Φ is the N × N diagonal matrix with Φij = 2πfδij(n · rk), n is the propagation direction of GW, and rk is the location of the k-th detector. Here Sn,k(f) is the one-side noise power spectral density of the k-th detector, The Fourier transform of the GW waveform of k-th detector is given by
$\begin{eqnarray}\begin{array}{rcl}{h}_{k}(f) & = & {{ \mathcal A }}_{k}{f}^{-7/6}{\rm{\exp }}\{{\rm{i}}[2\pi f{t}_{{\rm{c}}}-\pi /4-2{\psi }_{c}+2{\rm{\Psi }}(f/2)]\\ & & -\,{\varphi }_{k,(2,0)})\},\end{array}\end{eqnarray}$
where the Fourier amplitude can be written as
$\begin{eqnarray}\begin{array}{rcl}{{ \mathcal A }}_{k} & = & \displaystyle \frac{1}{{d}_{{\rm{L}}}}\sqrt{{({F}_{+,k}(1+{\cos }^{2}\iota ))}^{2}+{(2{F}_{\times ,k}\cos \iota )}^{2}}\\ & & \times \sqrt{5\pi /96}{\pi }^{-7/6}{{ \mathcal M }}_{{\rm{chirp}}}^{5/6}.\end{array}\end{eqnarray}$
Here, the detailed forms of $\Psi$(f/2) and φk,(2,0) can be found in [155, 157], dL is the luminosity distance of the GW source, ${{ \mathcal M }}_{{\rm{chirp}}}=(1+z){\eta }^{3/5}M$ is the observed chirp mass, η = m1m2/M2 is the symmetric mass ratio, and M = m1 + m2 is the total mass of the binary system with component masses m1 and m2, F+,k and F×,k are the antenna response functions of the k-th GW detector, we adopt the GW waveform in the frequency-domain, in which the time t is replaced by ${t}_{{\rm{f}}}={t}_{{\rm{c}}}-(5/256){{ \mathcal M }}_{{\rm{chirp}}}^{-5/3}{(\pi f)}^{-8/3}$ [155, 157], where tc is the coalescence time.
In this work, we consider the waveform in the inspiralling stage for the non-spinning BNS system. Here we adopt the restricted Post-Newtonian approximation and calculate the waveform to the 3.5 PN order [157, 158].
After simulating the GW catalog, we need to calculate the signal-to-noise ratio (SNR) for each GW event. The SNR for the detection network of N independent interferometers can be calculated by
$\begin{eqnarray}\rho ={(\tilde{{\boldsymbol{h}}}| \tilde{{\boldsymbol{h}}})}^{1/2}.\end{eqnarray}$
The inner product is defined as
$\begin{eqnarray}({\boldsymbol{a}}| {\boldsymbol{b}})=2{\int }_{{f}_{{\rm{lower}}}}^{{f}_{{\rm{upper}}}}\{{\boldsymbol{a}}(f){{\boldsymbol{b}}}^{* }(f)+{{\boldsymbol{a}}}^{* }(f){\boldsymbol{b}}(f)\}{\rm{d}}f,\end{eqnarray}$
where * represents conjugate transpose, a and b are column matrices of the same dimension, the lower cutoff frequency is set to flower = 1 Hz for ET and flower = 5 Hz for CE, and fupper = 2/(63/22πMobs) is the frequency at the last stable orbit with Mobs = (m1 + m2)(1 + z). In this work, we adopt the SNR threshold to be 12 in our simulation.
For the short GRB model, we adopt the model of Gaussian structured jet profile based on the GW170817/GRB170817A [159] observation,
$\begin{eqnarray}{L}_{{\rm{iso}}}({\theta }_{{\rm{v}}})={L}_{{\rm{on}}}\exp \left(-\displaystyle \frac{{\theta }_{{\rm{v}}}^{2}}{2{\theta }_{{\rm{c}}}^{2}}\right),\end{eqnarray}$
where θv is the viewing angle, Liso(θv) is the isotropically equivalent luminosity of short GRB observed at different θv, Lon = Liso(0) is the on-axis isotropic luminosity, θc = 4. 7 is the characteristic angle of the core, and the direction of the jet is assumed to align with the binary orbital angular momentum, namely ι = θv.
For the distribution of the short GRB, we assume the empirical broken-power-law luminosity function7
$\begin{eqnarray}{\rm{\Phi }}(L)\propto \left\{\begin{array}{ll}{(L/{L}_{* })}^{\alpha }, & L\lt {L}_{* },\\ {(L/{L}_{* })}^{\beta }, & L\geqslant {L}_{* },\end{array}\right.\end{eqnarray}$
which is commonly used in the literature [38, 42, 43, 46, 137, 138, 146, 148]. Here, L is the isotropic rest frame luminosity in the 1 − 10000 keV energy range, L* is the characteristic luminosity that separates the low and high end of the luminosity function, and the slopes describing these regimes are given by α and β, respectively. Following [146], we adopt L* = 2 × 1052 erg sec−1, α = − 1.95, and β = − 3. We assume a standard low end cutoff in luminosity of ${L}_{{\rm{\min }}}=1{0}^{49}$ erg sec−1, and we also term the on-axis isotropic luminosity Lon as the peak luminosity L [43, 137, 138, 160]. For the THESEUS mission [161], a GRB detection is recorded if the value of observed flux is greater than the flux threshold PT = 0.2phs−1cm−2 in the 50-300 keV band. For the GRB detection, we assume a duty cycle of 80% and a sky coverage fraction of 0.5. From the GW catalogue which has passed the threshold 12, we can select the GW-GRB events according to the probability distribution Φ(L)dL.
For a network with N independent interferometers, the Fisher information matrix is given by
$\begin{eqnarray}{F}_{ij}=\left(\left.\frac{\partial \tilde{{\boldsymbol{h}}}}{\partial {\theta }_{i}}\right|\frac{\partial \tilde{{\boldsymbol{h}}}}{\partial {\theta }_{j}}\right),\end{eqnarray}$
where θi denotes nine GW parameters (dL, ${{ \mathcal M }}_{{\rm{chirp}}}$, η, θ, φ, ι, tc, ψc, ψ) for a GW event.
For the total uncertainty of the luminosity distance dL, we first consider the instrumental error ${\sigma }_{{d}_{{\rm{L}}}}^{{\rm{inst}}}$. The covariance matrix is equal to the inverse of the Fisher information matrix, thus the instrumental error of GW parameter θi is
$\begin{eqnarray}{\rm{\Delta }}{\theta }_{i}=\sqrt{{({F}^{-1})}_{ii}},\end{eqnarray}$
where Fij is the total Fisher information matrix for the network of N interferometers.
In addition, the weak-lensing error ${\sigma }_{{d}_{{\rm{L}}}}^{{\rm{lens}}}$ and the peculiar velocity error ${\sigma }_{{d}_{{\rm{L}}}}^{{\rm{pv}}}$ are also considered. The error caused by weak lensing is adopted from [162, 163],
$\begin{eqnarray}\begin{array}{rcl}{\sigma }_{{d}_{{\rm{L}}}}^{{\rm{lens}}}(z) & = & \left[1-\frac{0.3}{\pi /2}\arctan (z/0.073)\right]\times {d}_{{\rm{L}}}(z)\\ & & \times 0.066{\left[\frac{1-{(1+z)}^{-0.25}}{0.25}\right]}^{1.8}.\end{array}\end{eqnarray}$
For the error caused by the peculiar velocity of the GW source is given by [164]
$\begin{eqnarray}{\sigma }_{{d}_{{\rm{L}}}}^{{\rm{pv}}}(z)={d}_{{\rm{L}}}(z)\times \left[1+\frac{c{(1+z)}^{2}}{H(z){d}_{{\rm{L}}}(z)}\right]\frac{\sqrt{\langle {v}^{2}\rangle }}{c},\end{eqnarray}$
where c is the speed of light in vacuum, H(z) is the Hubble parameter, and $\sqrt{\langle {v}^{2}\rangle }$ is the peculiar velocity of the GW source with respect to the Hubble flow is roughly set to $\sqrt{\langle {v}^{2}\rangle }=500\,{\rm{km}}{{\rm{s}}}^{-1}$.
Hence, the total error of dL can be written as
$\begin{eqnarray}{\sigma }_{{d}_{{\rm{L}}}}=\sqrt{{({\sigma }_{{d}_{{\rm{L}}}}^{{\rm{inst}}})}^{2}+{({\sigma }_{{d}_{{\rm{L}}}}^{{\rm{lens}}})}^{2}+{({\sigma }_{{d}_{{\rm{L}}}}^{{\rm{pv}}})}^{2}}.\end{eqnarray}$

2.2. Other cosmological observations

For comparison, we also employ three current EM cosmological observations, i.e., the cosmic microwave background (CMB) data, the baryon acoustic oscillation (BAO) data, and the type Ia supernova (SN) data. The details of these data are listed as follows.
The CMB data: the CMB likelihood including the TT, TE, EE spectra at l ≥ 30, the low-l temperature commander likelihood, and the low-l SimAll EE likelihood, from the Planck 2018 data release [165].
The BAO data: the measurements from 6dFGs at zeff = 0.106 [166], the SDSS-MGS at zeff = 0.15 [167], and BOSS-DR12 at zeff = 0.38, zeff = 0.51, and zeff = 0.61 [168].
The SN data: the Pantheon sample comprised of 1048 data points from the Pantheon compilation [169].

2.3. Methods of constraining cosmological parameters

In this paper, we will consider both cases of massless and massive sterile neutrinos in the framework of standard model (ΛCDM) of cosmology. For the ΛCDM model, the base parameter set (including six free parameters) is
$\begin{eqnarray*}{\bf{P}}=\{{\omega }_{b},\,{\omega }_{c},\,100{\theta }_{{\rm{MC}}},\,\tau ,\,\mathrm{ln}(1{0}^{10}{A}_{{\rm{s}}}),\,{n}_{{\rm{s}}}\},\end{eqnarray*}$
where ωb ≡ Ωbh2 and ωc ≡ Ωch2 are the physical densities of baryon and cold dark matter, respectively, θMC is the ratio (multiplied by 100) between the sound horizon and the angular diameter distance at the time of last-scattering, τ is the optical depth to the reionization, As is the amplitude of the primordial curvature perturbation, and ns is the scalar spectral index.
When we consider massless sterile neutrinos (as the dark radiation) in the ΛCDM model, an additional parameter Neff (the effective number of relativistic species) need to be added in the model, and this case is called ΛCDM+Neff model in this paper. When the massive sterile neutrinos are considered in the ΛCDM model, two extra free parameters, the Neff and ${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}$ (the effective sterile neutrino mass) need to be added in the model, and this case is called ΛCDM+Neff+${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}$ model in this paper. Thus, the ΛCDM+Neff model has seven independent parameters, and the ΛCDM+Neff+${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}$ has eight independent parameters. Note that in both the massless and massive sterile neutrino cases the total mass of active neutrinos is fixed at ∑mν = 0.06 eV.
For the GW standard siren observation with N data point, the χ2 function is defined as
$\begin{eqnarray}{\chi }_{{\rm{GW}}}^{2}=\displaystyle \sum _{i=1}^{N}{\left[\frac{{d}_{{\rm{L}}}^{i}-{d}_{{\rm{L}}}({z}_{i};\overrightarrow{{\rm{\Omega }}})}{{\sigma }_{{d}_{{\rm{L}}}}^{i}}\right]}^{2},\end{eqnarray}$
where zi, ${d}_{{\rm{L}}}^{i}$, and ${\sigma }_{{d}_{{\rm{L}}}}^{i}$ are the i-th GW redshift, luminosity distance, and the measurement error of the luminosity distance, respectively, $\overrightarrow{{\rm{\Omega }}}$ denotes the set of cosmological parameters.
In this work, we present the first forecast for the search for sterile neutrinos using joint GW-GRB observation. We use the public Markov-chain Monte Carlo (MCMC) package CosmoMC [170] to constrain sterile neutrino and other cosmological parameters. To demonstrate the impact of simulated GW data on constraining sterile neutrino parameters, we will consider all the different cases of 3G GW observations, the single ET, the single CE, the CE-CE network (one CE in the United States with 40-km arm length and another one in Australia with 20-km arm length, abbreviated as 2CE hereafter), and the ET-CE-CE network (one ET detector and two CE-like detectors, abbreviated as ET2CE hereafter) to analysis. We utilize the sensitivity curves of ET from [171] and for CE from [172], as shown in figure 1. For the GW detector, in view of the high uncertainty of the duty cycle, we only calculate the ideal scenario assuming a 100% duty cycle for all detectors, as discussed in [173]. The specific parameters characterizing the geometry of GW detector (latitude φ, longitude λ, opening angle ζ, and arm bisector angle γ) are detailed in table 1. The number of GW standard sirens in the subsequent cosmological analysis are shown in table 2 and their redshift distributions are shown in figure 2.
Figure 1. Sensitivity curves of the 3G GW detectors considered in this work.
Table 1. The specific coordinate parameters considered in this work.
GW detector φ(deg) λ(deg) γ(deg) ζ(deg)
Einstein Telescope, Europe 40.443 9.457 0.000 60
Cosmic Explorer, USA 43.827 −112.825 45.000 90
Cosmic Explorer, Australia −34.000 145.000 90.000 90
Table 2. Numbers of GW standard sirens in cosmological analysis, triggered by THESEUS in synergy with ET, CE, 2CE, and ET2CE, respectively.
Detection strategy ET CE 2CE ET2CE
Number of GW standard sirens 400 538 600 640
Figure 2. Redshift distributions of BNS detected by THESEUS in synergy with ET, CE, 2CE, and ET2CE for a 10-year observation.
The primary goal of this work is to assess the influence of joint observations between 3G GW detectors and future GRB detectors on the cosmological measurement of sterile neutrino parameters. Such observations are crucial for alleviating the degeneracies in cosmological parameters that are commonly observed in traditional EM data. To elucidate this, we have conducted simulations to generate mock GW data, which we have subsequently integrated with the mainstream EM observations, i.e., CMB+BAO+SN data. Our analysis specifically focuses on the estimation errors and precision of sterile neutrino parameters derived from this combined dataset.
To avoid any inconsistencies in the cosmological parameters constrained by combining CMB+BAO+SN with GW mock data, we have adopted a strategic approach. This approach is designed to thoroughly investigate the capacity of GW mock data to break the parameter degeneracies present in conventional EM observations. For this purpose, we have utilized the best-fit values of the cosmological parameters derived from the CMB+BAO+SN dataset as the reference values for simulating the GW mock data corresponding to each cosmological model. This methodology enables a more accurate assessment of the role of GW data in enhancing the precision of parameter estimation and resolving the degeneracies encountered in traditional EM observations.
For simplicity, we use 'CBS' to denote the joint CMB+BAO+SN data combination. Thus, in our analysis, we use five data combinations: (1) CBS, (2) CBS+ET, (3) CBS+CE, (4) CBS+2CE, and (5) CBS+ET2CE. We will report the constraint results in the next section.

3. Results and discussion

In this section, we report the constraint results for the ΛCDM+Neff and ΛCDM+Neff+${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}$ models using the CBS, CBS+ET, CBS+CE, CBS+2CE, and CBS+ET2CE data combinations and analyze how the GW standard sirens affect the cosmological constraints on the sterile neutrino parameters. The fitting results are shown in figures 3 and 4 and in tables 3 and 4. In the tables, we quote ± 1σ errors for the parameters, but for the parameters that cannot be well constrained, e.g., the sterile neutrino parameters Neff and ${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}$, we quote the 2σ upper limits. For a parameter ξ, we use σ(ξ) and ϵ(ξ) = σ(ξ)/ξ to represent its absolute error and relative error, respectively.
Figure 3. Constraint results for the ΛCDM+Neff model from the CBS, CBS+ET, CBS+CE, CBS+2CE, and CBS+ET2CE data combinations. One-dimensional marginalized posterior distribution for Neff (left panel), and two-dimensional marginalized posterior contours (1σ and 2σ) in the H0Neff and ΩmNeff planes (right panel).
Figure 4. Two-dimensional marginalized posterior contours (1σ and 2σ) in the Neff–${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}$, H0–${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}$, Ωm–${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}$ planes for the ΛCDM+Neff+${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}$ model from the constraints of the CBS, CBS+ET, CBS+CE, CBS+2CE, and CBS+ET2CE data combinations.
Table 3. Fitting results of the ΛCDM+Neff model by using the CBS, CBS+ET, CBS+CE, CBS+2CE, and CBS+ET2CE data combinations. We quote  ± 1σ errors for the parameters, but for the parameters that cannot be well constrained, we quote the 2σ upper limits. Here, H0 is in units of kms−1 Mpc−1.
Model CBS CBS+ET CBS+CE CBS+2CE CBS+ET2CE
Ωbh2 0.02247 ± 0.00015 0.02248 ± 0.00011 0.02248 ± 0.00011 0.02248 ± 0.00011 0.02248 ± 0.00011
Ωch2 $0.122{0}_{-0.0028}^{+0.0017}$ 0.1217 ± 0.0017 0.1217 ± 0.0016 0.1217 ± 0.0016 0.1217 ± 0.0016
100θMC $1.0404{6}_{-0.00037}^{+0.00043}$ 1.04050 ± 0.00036 $1.0405{0}_{-0.00035}^{+0.00034}$ 1.04049 ± 0.00035 1.04050 ± 0.00034
τ $0.055{8}_{-0.0083}^{+0.0074}$ $0.055{8}_{-0.0082}^{+0.0075}$ $0.055{9}_{-0.0082}^{+0.0074}$ $0.055{9}_{-0.0082}^{+0.0075}$ $0.055{8}_{-0.0082}^{+0.0072}$
ns $0.969{7}_{-0.0057}^{+0.0046}$ 0.9696 ± 0.0030 $0.969{6}_{-0.0031}^{+0.0030}$ 0.9697 ± 0.0030 0.9697 ± 0.0029
${\rm{ln}}(1{0}^{10}{A}_{s})$ $3.05{1}_{-0.018}^{+0.016}$ 3.050 ± 0.016 3.050 ± 0.016 3.050 ± 0.016 3.050 ± 0.016
σ8 $0.817{7}_{-0.0101}^{+0.0083}$ $0.817{0}_{-0.0081}^{+0.0082}$ 0.8169 ± 0.0078 0.8171 ± 0.0079 0.8169 ± 0.0077
Ωm $0.307{7}_{-0.0061}^{+0.0060}$ 0.3071 ± 0.0038 $0.307{1}_{-0.0036}^{+0.0035}$ 0.3071 ± 0.0035 $0.307{0}_{-0.0035}^{+0.0034}$
H0 $68.67{0}_{-1.020}^{+0.640}$ $68.67{1}_{-0.055}^{+0.057}$ $68.67{1}_{-0.057}^{+0.056}$ $68.67{1}_{-0.054}^{+0.052}$ $68.67{2}_{-0.050}^{+0.052}$
Neff < 3.464 3.209 ± 0.060 3.209 ± 0.056 $3.21{0}_{-0.056}^{+0.055}$ $3.20{8}_{-0.055}^{+0.054}$

ΔNeff > 0 . . .  2.717σ 2.911σ 2.929σ 2.945σ
σm) 0.00605 0.00380 0.00355 0.00350 0.00345
σ(H0) 0.8300 0.0560 0.0565 0.0530 0.0510
ϵm) 1.966% 1.237% 1.156% 1.140% 1.124%
ϵ(H0) 1.209% 0.082% 0.082% 0.077% 0.074%
Table 4. Fitting results of the ΛCDM+${N}_{{\rm{eff}}}+{m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}$ model by using the CBS, CBS+ET, CBS+CE, CBS+2CE, and CBS+ET2CE data combinations. We quote  ± 1σ errors for the parameters, but for the parameters that cannot be well constrained, we quote the 2σ upper limits. Here, H0 is in units of kms−1Mpc−1 and ${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}$ is in units of eV.
Model CBS CBS+ET CBS+CE CBS+2CE CBS+ET2CE
Ωbh2 $0.0224{7}_{-0.00016}^{+0.00015}$ 0.02249  ±  0.00012 0.02249  ±  0.00012 0.02249  ±  0.00012 0.02249  ±  0.00012
Ωch2 $0.119{8}_{-0.0031}^{+0.0036}$ $0.119{3}_{-0.0018}^{+0.0033}$ $0.119{3}_{-0.0017}^{+0.0032}$ $0.119{5}_{-0.0017}^{+0.0031}$ $0.119{4}_{-0.0017}^{+0.0030}$
100θMC $1.0405{9}_{-0.00033}^{+0.00047}$ $1.0407{1}_{-0.00035}^{+0.00041}$ $1.0407{1}_{-0.00035}^{+0.00040}$ $1.0407{0}_{-0.00036}^{+0.00040}$ $1.0407{1}_{-0.00036}^{+0.00039}$
τ $0.056{1}_{-0.0083}^{+0.0073}$ $0.056{6}_{-0.0082}^{+0.0074}$ $0.056{6}_{-0.0082}^{+0.0074}$ $0.056{8}_{-0.0083}^{+0.0074}$ $0.056{9}_{-0.0083}^{+0.0074}$
ns $0.968{1}_{-0.0065}^{+0.0047}$ 0.9684  ±  0.0033 0.9684  ±  0.0033 0.9684  ±  0.0033 $0.968{4}_{-0.0032}^{+0.0033}$
${\rm{ln}}(1{0}^{10}{A}_{s})$ $3.04{8}_{-0.018}^{+0.016}$ $3.04{8}_{-0.017}^{+0.016}$ $3.04{8}_{-0.017}^{+0.016}$ 3.048  ±  0.016 3.048  ±  0.016
σ8 $0.79{6}_{-0.014}^{+0.024}$ $0.79{5}_{-0.012}^{+0.022}$ $0.79{5}_{-0.012}^{+0.022}$ $0.79{6}_{-0.012}^{+0.021}$ $0.79{6}_{-0.012}^{+0.022}$
Ωm 0.3114  ±  0.0064 $0.309{6}_{-0.0052}^{+0.0032}$ $0.309{6}_{-0.0050}^{+0.0032}$ $0.309{7}_{-0.0049}^{+0.0033}$ $0.309{6}_{-0.0048}^{+0.0033}$
H0 $68.15{0}_{-1.000}^{+0.500}$ $68.15{6}_{-0.053}^{+0.055}$ $68.15{6}_{-0.053}^{+0.054}$ $68.15{6}_{-0.053}^{+0.054}$ $68.15{6}_{-0.050}^{+0.054}$
Neff <  3.446 $3.14{8}_{-0.094}^{+0.039}$ $3.14{8}_{-0.090}^{+0.040}$ $3.15{0}_{-0.089}^{+0.043}$ $3.15{0}_{-0.089}^{+0.042}$
${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}$ <  0.5789 <  0.4842 <  0.4772 <  0.4321 <  0.4226

ΔNeff > 0 ... 1.085σ 1.133σ 1.169σ 1.169σ
σm) 0.00640 0.00420 0.00410 0.00410 0.00405
σ(H0) 0.7500 0.0540 0.0535 0.0535 0.0520
ϵm) 2.055% 1.357% 1.324% 1.324% 1.308%
ϵ(H0) 1.101% 0.079% 0.078% 0.078% 0.076%
In accordance with the cosmological parameter constraints obtained from the CBS data, the central values and uncertainties of the cosmological parameters from the combined CBS and GW observations are detailed in tables 3 and 4. Given that the CBS data is real, whereas the GW data is simulated, the constraints from the joint CBS and GW analysis represent a mixture of real and simulated data. Consequently, the central values should not be interpreted as from actual observational data. As such, this work emphasizes the significance of the uncertainties and precision in the derived cosmological constraints.

3.1. The case of massless sterile neutrino

The massless sterile neutrinos serve as the dark radiation, and thus in this case Neff is treated as a free parameter, the total relativistic energy density of radiation is given by
$\begin{eqnarray}{\rho }_{{\rm{r}}}=\left[1+{N}_{{\rm{eff}}}\frac{7}{8}{\left(\frac{4}{11}\right)}^{\frac{4}{3}}\right]{\rho }_{\gamma },\end{eqnarray}$
where ργ is the photon energy density. In the ΛCDM model, Neff = 3.046, and so ΔNeff = Neff − 3.046 > 0 indicates the presence of extra relativistic particle species in the early universe, and in this paper we take the fit results of ΔNeff > 0 as evidence of the existence of massless sterile neutrinos.
In table 3, we find that the CBS data provides only an upper limit, Neff < 3.464, indicating that the existence of massless sterile neutrinos is not favored by the CBS data. However, when GW standard sirens are included in the data combination, the results change significantly. After considering the GW data, we obtain Neff = 3.209 ± 0.060 for CBS+ET, Neff = 3.209 ± 0.056 for CBS+CE, ${N}_{{\rm{eff}}}=3.21{0}_{-0.056}^{+0.055}$ for CBS+2CE, and ${N}_{{\rm{eff}}}=3.20{8}_{-0.055}^{+0.054}$ for CBS+ET2CE, which indicates a detection of ΔNeff > 0 at the 2.717σ, 2.911σ, 2.929σ, and 2.945σ significance levels, respectively. Obviously, the GW data can indeed effectively improve the constraints on the massless sterile neutrino parameter Neff.
Note that here CBS refers to current real observational data, while GW denotes simulated data. When we use the actual CBS data to conduct constraints, we can only obtain an upper limit for Neff and not a definitive detection result. However, when simulated GW data is included, the errors in the parameter constraints are significantly reduced. With the central value of Neff remaining essentially unchanged and ΔNeff substantially decreased, we can achieve a result of ΔNeff > 0. Our results indicate that future GW observations can greatly assist in improving the cosmological detection of massless sterile neutrinos, with a significance level reaching up to 3σ. Of course, we must remember that our findings are based on simulated outcomes, and the central value of Neff is set to be consistent with the result from CBS.
In the right panel of figure 3, we show the two-dimensional posterior distribution contours (1σ and 2σ) in the NeffH0 and Neff–Ωm planes for the ΛCDM+Neff model using CBS, CBS+ET, CBS+CE, CBS+2CE, and CBS+ET2CE data combinations. We can see that Neff is in positive correlation with H0 by using the CBS data combination. However, after adding the GW data, this correlation becomes negligible, indicating that the degeneracy between Neff and H0 is effectively broken by the GW observations. In addition, we can also clearly see that when considering the GW data, the parameter space is greatly shrunk in each plane and the constraints on cosmological parameters of H0 and Ωm become much tighter.
In table 3, we also show absolute and relative errors of H0 and Ωm from the CBS, CBS+ET, CBS+CE, CBS+2CE, and CBS+ET2CE data combinations. Compared to the CBS data, we find that the accuracy of the H0 constraint improves by 93.253% for CBS+ET, 93.193% for CBS+CE, 93.614% for CBS+2CE, and 93.855% for CBS+ET2CE, respectively. Similarly, the accuracy of the Ωm constraint improves by 37.190% for CBS+ET, 41.322% for CBS+CE, 42.149% for CBS+2CE, and 42.975% for CBS+ET2CE, respectively. Obviously, the GW data can indeed effectively improve the constraints on the parameters H0 and Ωm.

3.2. The case of massive sterile neutrino

In this subsection, we investigated how GW standard sirens constrain the massive sterile neutrino parameters. Hence, the requirement of Neff > 3.046 still holds.
From table 4, we obtain Neff < 3.446 by using CBS data. After adding the GW data, the constraint results become ${N}_{{\rm{eff}}}=3.14{8}_{-0.094}^{+0.039}$ for CBS+ET, ${N}_{{\rm{eff}}}=3.14{8}_{-0.090}^{+0.039}$ for CBS+CE, ${N}_{{\rm{eff}}}=3.15{0}_{-0.089}^{+0.043}$ for CBS+2CE, and ${N}_{{\rm{eff}}}=3.15{0}_{-0.089}^{+0.042}$ for CBS+ET2CE, respectively. We find that Neff cannot be well constrained using the CBS data, but the addition of GW data can significantly improve the constraint on Neff, favoring ΔNeff > 0 at 1.085σ (CBS+ET), 1.133σ (CBS+CE), 1.169σ (CBS+2CE), and 1.169σ (CBS+ET2CE) statistical significance, respectively. For the mass of sterile neutrino, the CBS data give ${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}\lt 0.5789$ eV and further including the GW data leads to results of ${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}\lt 0.4842$ eV (CBS+ET), ${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}\lt 0.4772$ eV (CBS+CE), ${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}\lt 0.4321$ eV (CBS+2CE), and ${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}\lt 0.4226$ eV (CBS+ET2CE), respectively. Evidently, adding GW data tightens the constraint on ${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}$ significantly, which is in accordance with the conclusions in previous studies' active neutrinos mass by using the GW data [14, 37, 46]. Therefore, the GW data also plays an important role in constraining the mass of sterile neutrinos.
In figure 4, we show two-dimensional marginalized posterior contours (1σ and 2σ) in the Neff–${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}$, H0–${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}$, Ωm–${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}$ planes for the ΛCDM+Neff+${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}$ model. We can clearly see that when further considering the GW data, the parameter space is also greatly shrunk and the constraints on H0 and Ωm also become much tighter. From table 4, we find that the constraints on H0 and Ωm could be improved by 92.800% and 34.375%, respectively, when adding the ET data to the CBS data, 92.867% and 35.938% for the case of CE, 92.867% and 35.938% for the case of 2CE, and 93.067% and 36.719% for the case of ET2CE, respectively. These results are in accordance with the conclusions for both the case of massless sterile neutrinos in this study and active neutrinos mass in previous studies [14, 37, 46]. Therefore, the inclusion of GW data can significantly improve constraints on most cosmological parameters, particularly H0 and Ωm.

4. Conclusion

This work aims to forecast the search for sterile neutrinos using joint GW-GRB observations. We consider two cases of massless and massive sterile neutrinos, corresponding to the ΛCDM+Neff and ΛCDM+Neff+${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}$ models, respectively. We consider four GW detection observation strategies, i.e., ET, CE, the 2CE network, and the ET2CE network, to perform cosmological analysis. To evaluate the impact of GW data on the constraints of sterile neutrino parameters, we also considered existing CMB+BAO+SN data for comparison and combination.
For the ΛCDM+Neff model, in the case of using CMB+BAO+SN, only upper limits on Neff can be obtained. Further adding the GW data tightens the Neff significantly, and in this case the preference of ΔNeff > 0 at about 3σ level. Therefore, GW standard siren observations can greatly assist in the detection of massless sterile neutrinos.
For the ΛCDM+Neff+${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}$ model, only upper limits on Neff and ${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}$ can be derived by using the CMB+BAO+SN data. Further including GW data significantly improves the constraints, and we find that the GW data gives a rather tight upper limit on ${m}_{\nu ,{\rm{sterile}}}^{{\rm{eff}}}$ and favor ΔNeff > 0 at about 1.1σ level. These results also seem to favor massless sterile neutrinos.
Furthermore, we find that the GW data can significantly enhance the accuracy of constraints on the derived parameters H0 and Ωm. The accuracy of H0 improves by approximately 93% and the accuracy of Ωm increases by about 37% to 42%, when the GW data is included in the cosmological fit.

This work was supported by the National Natural Science Foundation of China under Grant Nos. 12305069, 11947022, 12473001, 11975072, 11875102, and 11835009, the National SKA Program of China under Grants Nos. 2022SKA0110200 and 2022SKA0110203, the Program of the Education Department of Liaoning Province under Grant No. JYTMS20231695, and the National 111 Project under Grant No. B16009.

1
Abbott B P et al 2017 GW170817: Observation of gravitational waves from a binary neutron star inspiral Phys. Rev. Lett. 119 161101

DOI

2
Abbott B P et al 2017 Multi-messenger observations of a binary neutron star merger Astrophys. J. Lett. 848 l12

DOI

3
Dìaz M C et al 2017 Observations of the first electromagnetic counterpart to a gravitational wave source by the TOROS collaboration Astrophys. J. Lett. 848 l29

DOI

4

5
Punturo M et al 2010 The Einstein telescope: A third-generation gravitational wave observatory Class. Quant. Grav. 27 194002

DOI

6

7
Benjamin P et al 2017 Exploring the sensitivity of next generation gravitational wave detectors Class. Quant. Grav 34 044001

DOI

8
Cai R-G, Yang T 2017 Estimating cosmological parameters by the simulated data of gravitational waves from the Einstein Telescope Phys. Rev. D 95 044024

DOI

9
Cai R-G, Yang T 2018 Standard sirens and dark sector with Gaussian process EPJ Web Conf. 168 01008

DOI

10
Liu T, Zhang X, Zhao W 2018 Constraining f(R) gravity in solar system, cosmology and binary pulsar systems Phys. Lett. B 777 286293

DOI

11
Cai R-G, Liu T-B, Liu X-W, Wang S-J, Yang T 2018 Probing cosmic anisotropy with gravitational waves as standard sirens Phys. Rev. D 97 103005

DOI

12
Berti E, Yagi K, Yunes N 2018 Extreme gravity tests with gravitational waves from compact binary coalescences: (I) inspiral-merger Gen. Rel. Grav. 50 46

DOI

13
Cai Y-F, Li C, Saridakis E N, Xue L 2018 f(T) gravity after GW170817 and GRB170817A Phys. Rev. D 97 103513

DOI

14
Wang L-F, Zhang X-N, Zhang J-F, Zhang X 2018 Impacts of gravitational-wave standard siren observation of the Einstein Telescope on weighing neutrinos in cosmology Phys. Lett. B 782 87 93

DOI

15
Zhao W, Wright B S, Li B 2018 Constraining the time variation of Newton's constant G with gravitational-wave standard sirens and supernovae J. Cosmol. Astropart. Phys. 2018 052

DOI

16
Zhang X-N, Wang L-F, Zhang J-F, Zhang X 2019 Improving cosmological parameter estimation with the future gravitational-wave standard siren observation from the Einstein Telescope Phys. Rev. D 99 063510

DOI

17
Du M, Yang W, Xu L, Pan S, Mota D F 2019 Future constraints on dynamical dark-energy using gravitational-wave standard sirens Phys. Rev. D 100 043535

DOI

18
He J-H 2019 Accurate method to determine the systematics due to the peculiar velocities of galaxies in measuring the Hubble constant from gravitational-wave standard sirens Phys. Rev. D 100 023527

DOI

19
Yang W, Pan S, Valentino E D, Wang B, Wang A 2020 Forecasting interacting vacuum-energy models using gravitational waves J. Cosmol. Astropart. Phys. 2020 050

DOI

20
Yang W, Vagnozzi S, Valentino E D, Nunes R C, Pan S, Mota D F 2019 Listening to the sound of dark sector interactions with gravitational wave standard sirens J. Cosmol. Astropart. Phys. 2019 037

DOI

21
Zhang X 2019 Gravitational wave standard sirens and cosmological parameter measurement Sci. China Phys. Mech. Astron. 62 110431

DOI

22
Zhang J-F, Dong H-Y, Qi J-Z, Zhang X 2020 Prospect for constraining holographic dark energy with gravitational wave standard sirens from the Einstein Telescope Eur. Phys. J. C 80 217

DOI

23
Bachega R R A, Costa A A, Abdalla E, Fornazier K S F 2020 Forecasting the interaction in dark matter-dark energy models with standard sirens from the Einstein Telescope J. Cosmol. Astropart. Phys. 2020 021

DOI

24
Wang L-F, Zhao Z-W, Zhang J-F, Zhang X 2020 A preliminary forecast for cosmological parameter estimation with gravitational-wave standard sirens from TianQin J. Cosmol. Astropart. Phys. 2020 012

DOI

25
Zhang J-F, Zhang M, Jin S-J, Qi J-Z, Zhang X 2019 Cosmological parameter estimation with future gravitational wave standard siren observation from the Einstein Telescope J. Cosmol. Astropart. Phys. 2019 068

DOI

26
Li H-L, He D-Z, Zhang J-F, Zhang X 2020 Quantifying the impacts of future gravitational-wave data on constraining interacting dark energy J. Cosmol. Astropart. Phys. 2020 038

DOI

27
Zhao Z-W, Wang L-F, Zhang J-F, Zhang X 2020 Prospects for improving cosmological parameter estimation with gravitational-wave standard sirens from Taiji Sci. Bull. 65 1340 1348

DOI

28
Jin S-J, He D-Z, Xu Y, Zhang J-F, Zhang X 2020 Forecast for cosmological parameter estimation with gravitational-wave standard siren observation from the Cosmic Explorer J. Cosmol. Astropart. Phys. 2020 051

DOI

29
Wang L-F, Jin S-J, Zhang J-F, Zhang X 2022 Forecast for cosmological parameter estimation with gravitational-wave standard sirens from the LISA-Taiji network Sci. China Phys. Mech. Astron. 65 210411

DOI

30
Qi J-Z, Jin S-J, Fan X-L, Zhang J-F, Zhang X 2021 Using a multi-messenger and multi-wavelength observational strategy to probe the nature of dark energy through direct measurements of cosmic expansion history J. Cosmol. Astropart. Phys. 2021 042

DOI

31
Jin S-J, Wang L-F, Wu P-J, Zhang J-F, Zhang X 2021 How can gravitational-wave standard sirens and 21-cm intensity mapping jointly provide a precise late-universe cosmological probe? Phys. Rev. D 104 103507

DOI

32
Zhu L-G, Xie L-H, Hu Y-M, Liu S, Li E-K, Napolitano N R, Tang B-T, Zhang J, Mei J 2022 Constraining the Hubble constant to a precision of about 1% using multi-band dark standard siren detections Sci. China Phys. Mech. Astron. 65 259811

DOI

33
Mendonça J, de Souza S, Sturani R, Alcaniz J 2022 Cosmography with standard sirens from the Einstein Telescope J. Cosmol. Astropart. Phys. 2022 025

DOI

34
Wang L-F, Shao Y, Zhang J-F, Zhang X 2022 Ultra-low-frequency gravitational waves from individual supermassive black hole binaries as standard sirens arXiv:2201.00607

35
Wu P-J, Shao Y, Jin S-J, Zhang X 2023 A path to precision cosmology: synergy between four promising late-universe cosmological probes J. Cosmol. Astropart. Phys. 06 052

DOI

36
Jin S-J, Li T-N, Zhang J-F, Zhang X 2023 Prospects for measuring the Hubble constant and dark energy using gravitational-wave dark sirens with neutron star tidal deformation J. Cosmol. Astropart. Phys. 2023 070

DOI

37
Jin S-J, Zhu R-Q, Wang L-F, Li H-L, Zhang J-F, Zhang X 2022 Impacts of gravitational-wave standard siren observations from Einstein Telescope and Cosmic Explorer on weighing neutrinos in interacting dark energy models Commun. Theor. Phys. 74 105404

DOI

38
Hou W-T, Qi J-Z, Han T, Zhang J-F, Cao S, Zhang X 2023 Prospects for constraining interacting dark energy models from gravitational wave and gamma ray burst joint observation J. Cosmol. Astropart. Phys. 2023 017

DOI

39
Song J-Y, Wang L-F, Li Y, Zhao Z-W, Zhang J-F, Zhao W, Zhang X 2024 Synergy between CSST galaxy survey and gravitational-wave observation: inferring the Hubble constant from dark standard sirens Sci. China Phys. Mech. Astron. 67 230411

DOI

40
Jin S-J, Xing S-S, Shao Y, Zhang J-F, Zhang X 2023 Joint constraints on cosmological parameters using future multi-band gravitational wave standard siren observations Chin. Phys. C 47 065104

DOI

41
Jin S-J, Zhang Y-Z, Song J-Y, Zhang J-F, Zhang X 2024 Taiji-TianQin-LISA network: precisely measuring the Hubble constant using both bright and dark sirens Sci. China Phys. Mech. Astron. 67 220412

DOI

42
Jin S-J, Zhu R-Q, Song J-Y, Han T, Zhang J-F, Zhang X 2023 Standard siren cosmology in the era of the 2.5-generation ground-based gravitational wave detectors: bright and dark sirens of LIGO Voyager and NEMO J. Cosmol. Astropart. Phys. 2024 050

DOI

43
Han T, Jin S-J, Zhang J-F, Zhang X 2024 A comprehensive forecast for cosmological parameter estimation using joint observations of gravitational waves and short -ray bursts Eur. Phys. J. C 84 663

DOI

44
Li T-N, Jin S-J, Li H-L, Zhang J-F, Zhang X 2024 Prospects for probing the interaction between dark energy and dark matter using gravitational-wave dark sirens with neutron star tidal deformation Astrophys. J. 963 52

DOI

45
Dong Y-Y, Song J-Y, Jin S-J, Zhang J-F, Zhang X 2024 Enhancing dark siren cosmology through multi-band gravitational wave synergetic observations arXiv:2404.18188

46
Feng L, Han T, Zhang J-F, Zhang X 2024 Prospects for weighing neutrinos in interacting dark energy models using joint observations of gravitational waves and γ-ray bursts Chin. Phys. C 48 095104

DOI

47
Bian L et al 2021 The gravitational-wave physics II: progress Sci. China Phys. Mech. Astron. 64 120401

DOI

48
Aguilar A et al 2001 Evidence for neutrino oscillations from the observation of ${\bar{\nu }}_{e}$ appearance in a ${\bar{\nu }}_{e}$ beam Phys. Rev. D 64 112007

DOI

49
Giunti C, Laveder M 2011 Statistical significance of the gallium anomaly Phys. Rev. C 83 065504

DOI

50
Akbar M, Quevedo H, Saifullah K, Sanchez A, Taj S 2011 Thermodynamic geometry of charged rotating BTZ black holes Phys. Rev. D 83 084031

DOI

51
Conrad J M, Ignarra C M, Karagiorgi G, Shaevitz M H, Spitz J 2013 Sterile neutrino fits to short baseline neutrino oscillation measurements Adv. High Energy Phys. 163897 2013

DOI

52
Aguilar-Arevalo A A et al 2012 A combined νμνe and ${\bar{\nu }}_{\mu }\to {\bar{\nu }}_{e}$ oscillation analysis of the MiniBooNE excesses arXiv:1207.4809

53
Giunti C, Laveder M, Li Y F, Liu Q Y, Long H W 2012 Update of short-baseline electron neutrino and antineutrino disappearance Phys. Rev. D 86 113014

DOI

54
Giunti C, Laveder M, Li Y F, Long H W 2013 Short-baseline electron neutrino oscillation length after troitsk Phys. Rev. D 87 013004

DOI

55
Kopp J, Machado P A N, Maltoni M, Schwetz T 2013 Sterile neutrino oscillations: the global picture J. High Energy Phys. 05 050

DOI

56
Giunti C, Laveder M, Li Y F, Long H W 2013 Pragmatic view of short-baseline neutrino oscillations Phys. Rev. D 88 073008

DOI

57
Gariazzo S, Giunti C, Laveder M 2013 Light sterile neutrinos in cosmology and short-baseline oscillation experiments J. High Energy Phys. 2013 JHEP11(2013)211

DOI

58
Abazajian K N et al 2012 Light sterile neutrinos: a white paper arXiv:1204.5379

59
Hannestad S, Tamborra I, Tram T 2012 Thermalisation of light sterile neutrinos in the early Universe J. Cosmol. Astropart. Phys. 2012 025

DOI

60
Conrad J M, Louis W C, Shaevitz M H 2012 The LSND and MiniBooNE oscillation searches at high Δm2 Ann. Rev. Nucl. Part. Sci. 63 4567

DOI

61
Hu W, Eisenstein D J, Tegmark M 1998 Weighing neutrinos with galaxy surveys Phys. Rev. Lett. 80 5255 5258

DOI

62
Reid B A, Verde L, Jimenez R, Mena O 2010 Robust neutrino constraints by combining low redshift observations with the CMB J. Cosmol. Astropart. Phys. 2010 003

DOI

63
Li H, Zhang X 2012 Constraining dynamical dark energy with a divergence-free parametrization in the presence of spatial curvature and massive neutrinos Phys. Lett. B 713 160 164

DOI

64
Wang X, Meng X-L, Zhang T-J, Shan H Y, Gong Y, Tao C, Chen X, Huang Y F 2012 Observational constraints on cosmic neutrinos and dark energy revisited J. Cosmol. Astropart. Phys. 2012 018

DOI

65
Hamann J, Hannestad S, Wong Y Y Y 2012 Measuring neutrino masses with a future galaxy survey J. Cosmol. Astropart. Phys. 2012 052

DOI

66
Li Y-H, Wang S, Li X-D, Zhang X 2013 Holographic dark energy in a Universe with spatial curvature and massive neutrinos: a full Markov Chain Monte Carlo exploration J. Cosmol. Astropart. Phys. 02 033

DOI

67
Riemer-Sørensen S, Parkinson D, Davis T M 2014 Combining Planck data with large scale structure information gives a strong neutrino mass constraint Phys. Rev. D 89 103505

DOI

68
Giusarma E, de Putter R, Ho S, Mena O 2013 Constraints on neutrino masses from Planck and Galaxy clustering data Phys. Rev. D 88 063515

DOI

69
Cahn R N, Dwyer D A, Freedman S J, Haxton W C, Kadel R W, Kolomensky Y G, Luk K B, McDonald P, Orebi Gann G D, Poon A W P 2013 White paper: measuring the neutrino mass hierarchy arXiv:1307.5487

70
Lesgourgues J, Pastor S 2014 Neutrino cosmology and Planck New J. Phys. 16 065002

DOI

71
Zhang J-F, Li Y-H, Zhang X 2014 Cosmological constraints on neutrinos after BICEP2 Eur. Phys. J. C 74 2954

DOI

72
Zhou X-Y, He J-H 2014 Weighing neutrinos in f(R) gravity in light of BICEP2 Commun. Theor. Phys. 62 102108

DOI

73
Costanzi M, Sartoris B, Viel M, Borgani S 2014 Neutrino constraints: what large-scale structure and CMB data are telling us? J. Cosmol. Astropart. Phys. 2014 081

DOI

74
Nathalie P-D et al 2014 Constraint on neutrino masses from SDSS-III/BOSS Lyα forest and other cosmological probes J. Cosmol. Astropart. Phys. 02 045

DOI

75
Zhang J-F, Zhao M-M, Li Y-H, Zhang X 2015 Neutrinos in the holographic dark energy model: constraints from latest measurements of expansion history and growth of structure J. Cosmol. Astropart. Phys. 04 038

DOI

76
Qian X, Vogel P 2015 Neutrino mass hierarchy Prog. Part. Nucl. Phys. 83 1 30

DOI

77
Patterson R B 2015 Prospects for measurement of the neutrino mass hierarchy Ann. Rev. Nucl. Part. Sci. 65 177 192

DOI

78
Allison R, Caucal P, Calabrese E, Dunkley J, Louis T 2015 Towards a cosmological neutrino mass detection Phys. Rev. D 92 123535

DOI

79
Geng C-Q, Lee C-C, Myrzakulov R, Sami M, Saridakis E N 2016 Observational constraints on varying neutrino-mass cosmology J. Cosmol. Astropart. Phys. 2016 049

DOI

80
Chen Y, Xu L 2016 Galaxy clustering, CMB and supernova data constraints on φ CDM model with massive neutrinos Phys. Lett. B 752 6675

DOI

81
Zhang X 2015 Impacts of dark energy on weighing neutrinos after Planck Phys. Rev. D 93 083011

DOI

82
Huang Q-G, Wang K, Wang S 2016 Constraints on the neutrino mass and mass hierarchy from cosmological observations Eur. Phys. J. C 76 489

DOI

83
Chen Y, Ratra B, Biesiada M, Li S, Zhu Z-H 2016 Constraints on non-flat cosmologies with massive neutrinos after Planck 2015 Astrophys. J. 829 61

DOI

84
Moresco M, Jimenez R, Verde L, Cimatti A, Pozzetti L, Maraston C, Thomas D 2016 Constraining the time evolution of dark energy, curvature and neutrino properties with cosmic chronometers J. Comsol. Astropart. Phys. 2016 039

DOI

85
Lu J, Liu M, Wu Y, Wang Y, Yang W 2016 Cosmic constraint on massive neutrinos in viable f(R) gravity with producing ΛCDM background expansion Eur. Phys. J. C 76 679

DOI

86
Hada R, Futamase T 2016 Constraints on neutrino masses from the lensing dispersion of Type Ia supernovae Astrophys. J. 828 112

DOI

87
Wang S, Wang Y-F, Xia D-M, Zhang X 2016 Impacts of dark energy on weighing neutrinos: mass hierarchies considered Phys. Rev. D 94 083519

DOI

88
Kumar S, Nunes R C 2016 Probing the interaction between dark matter and dark energy in the presence of massive neutrinos Phys. Rev. D 94 123511

DOI

89
Zhao M-M, Li Y-H, Zhang J-F, Zhang X 2016 Constraining neutrino mass and extra relativistic degrees of freedom in dynamical dark energy models using Planck 2015 data in combination with low-redshift cosmological probes: basic extensions to ΛCDM cosmology Mon. Not. Roy. Astron. Soc. 469 17131724

DOI

90
Böhringer H, Chon G 2016 Constraints on neutrino masses from the study of the nearby large-scale structure and galaxy cluster counts Mod. Phys. Lett. A 31 1640008

DOI

91
Xu L, Huang Q-G 2018 Detecting the neutrinos mass hierarchy from cosmological data Sci. China Phys. Mech. Astron. 61 039521

DOI

92
Vagnozzi S, Giusarma E, Mena O, Freese K, Gerbino M, Ho S, Lattanzi M 2017 Unveiling ν secrets with cosmological data: neutrino masses and mass hierarchy Phys. Rev. D 96 123503

DOI

93
Guo R-Y, Li Y-H, Zhang J-F, Zhang X 2017 Weighing neutrinos in the scenario of vacuum energy interacting with cold dark matter: application of the parameterized post-Friedmann approach J. Cosmol. Astropart. Phys. 05 040

DOI

94
Zhang X 2017 Weighing neutrinos in dynamical dark energy models Sci. China Phys. Mech. Astron. 60 060431

DOI

95
Chen L, Huang Q-G, Wang K 2017 New cosmological constraints with extended-Baryon Oscillation Spectroscopic Survey DR14 quasar sample Eur. Phys. J. C 77 762

DOI

96
Yang W, Nunes R C, Pan S, Mota D F 2017 Effects of neutrino mass hierarchies on dynamical dark energy models Phys. Rev. D 95 103522

DOI

97
Koksbang S M, Hannestad S 2017 Constraining dynamical neutrino mass generation with cosmological data J. Cosmol. Astropart. Phys. 2017 014

DOI

98
Li E-K, Zhang H, Du M, Zhou Z-H, Xu L 2018 Probing the neutrino mass hierarchy beyond ΛCDM model J. Cosmol. Astropart. Phys. 2018 042

DOI

99
Wang S, Wang Y-F, Xia D-M 2018 Constraints on the sum of neutrino masses using cosmological data including the latest extended Baryon Oscillation Spectroscopic Survey DR14 quasar sample Chin. Phys. C 42 065103

DOI

100
Zhao M-M, Zhang J-F, Zhang X 2018 Measuring growth index in a universe with massive neutrinos: a revisit of the general relativity test with the latest observations Phys. Lett. B 779 473 478

DOI

101
Boyle A, Komatsu E 2018 Deconstructing the neutrino mass constraint from galaxy redshift surveys J. Comsol. Astropart. Phys. 2018 035

DOI

102
Vagnozzi S, Dhawan S, Gerbino M, Freese K, Goobar A, Mena O 2018 Constraints on the sum of the neutrino masses in dynamical dark energy models with w(z)≥ − 1 are tighter than those obtained in ΛCDM Phys. Rev. D 98 083501

DOI

103
Guo R-Y, Zhang J-F, Zhang X 2018 Exploring neutrino mass and mass hierarchy in the scenario of vacuum energy interacting with cold dark matte. Chin. Phys. C 42 095103

DOI

104
Choudhury S R, Choubey S 2018 Updated bounds on sum of neutrino masses in various cosmological scenarios J. Comsol. Astropart. Phys. 2018 017

DOI

105
Feng L, Li H-L, Zhang J-F, Zhang X 2020 Exploring neutrino mass and mass hierarchy in interacting dark energy models Sci. China Phys. Mech. Astron. 63 220401

DOI

106
Zhang J-F, Wang B, Zhang X 2020 Forecast for weighing neutrinos in cosmology with SKA Sci. China Phys. Mech. Astron. 63 280411

DOI

107
Li H-L, Zhang J-F, Zhang X 2018 Constraints on neutrino mass in the scenario of vacuum energy interacting with cold dark matter after Planck Commun. Theor. Phys. 72 125401

DOI

108
Zhang M, Zhang J-F, Zhang X 2020 Impacts of dark energy on constraining neutrino mass after Planck 2018 Commun. Theor. Phys. 72 125402

DOI

109
de Holanda P C, Smirnov A Y 2011 Solar neutrino spectrum, sterile neutrinos and additional radiation in the Universe Phys. Rev. D 83 113011

DOI

110
Palazzo A 2013 Phenomenology of light sterile neutrinos: a brief review Mod. Phys. Lett. A 28 1330004

DOI

111
Hamann J, Hasenkamp J 2013 A new life for sterile neutrinos: resolving inconsistencies using hot dark matter J. Cosmol. Astropart. Phys. 10 044

DOI

112
Wyman M, Rudd D H, Vanderveld R A, Hu W 2014 Neutrinos help reconcile Planck measurements with the local Universe Phys. Rev. Lett. 112 051302

DOI

113
Battye R A, Moss A 2014 Evidence for massive neutrinos from cosmic microwave background and lensing observations Phys. Rev. Lett. 112 051303

DOI

114
Dvorkin C, Wyman M, Rudd D H, Hu W 2014 Neutrinos help reconcile Planck measurements with both the early and local Universe Phys. Rev. D 90 083503

DOI

115
Archidiacono M, Fornengo N, Gariazzo S, Giunti C, Hannestad S, Laveder M 2014 Light sterile neutrinos after BICEP-2 J. Cosmol. Astropart. Phys. 2014 031

DOI

116
Ko P, Tang Y 2014 νΛMDM: A model for sterile neutrino and dark matter reconciles cosmological and neutrino oscillation data after BICEP2 Phys. Lett. B 739 6267

DOI

117
Li Y-H, Zhang J-F, Zhang X 2014 Tilt of primordial gravitational wave spectrum in a universe with sterile neutrinos Sci. China Phys. Mech. Astron. 57 1455 1459

DOI

118
Zhang J-F, Li Y-H, Zhang X 2015 Sterile neutrinos help reconcile the observational results of primordial gravitational waves from Planck and BICEP2 Phys. Lett. B 740 359 363

DOI

119
Archidiacono M, Hannestad S, Hansen R S, Tram T 2015 Cosmology with self-interacting sterile neutrinos and dark matter - A pseudoscalar model Phys. Rev. D 91 065021

DOI

120
Bergström J, Gonzalez-Garcia M C, Niro V, Salvado J 2014 Statistical tests of sterile neutrinos using cosmology and short-baseline data J. High Energy Phys. 2014 JHEP10(2014)104

DOI

121
An F P et al 2014 Search for a light sterile neutrino at daya bay Phys. Rev. Lett. 113 141802

DOI

122
Zhang J-F, Geng J-J, Zhang X 2014 Neutrinos and dark energy after Planck and BICEP2: data consistency tests and cosmological parameter constraints J. Cosmol. Astropart. Phys. 10 044

DOI

123
Zhang J-F, Li Y-H, Zhang X 2014 Measuring growth index in a universe with sterile neutrinos Phys. Lett. B 739 102 105

DOI

124
Li Y-H, Zhang J-F, Zhang X 2015 Probing f(R) cosmology with sterile neutrinos via measurements of scale-dependent growth rate of structure Phys. Lett. B 744 213 217

DOI

125
Feng L, Zhang J-F, Zhang X 2017 A search for sterile neutrinos with the latest cosmological observations Eur. Phys. J. C 77 418

DOI

126
Zhao M-M, He D-Z, Zhang J-F, Zhang X 2017 Search for sterile neutrinos in holographic dark energy cosmology: Reconciling Planck observation with the local measurement of the Hubble constant Phys. Rev. D 96 043520

DOI

127
Feng L, Zhang J-F, Zhang X 2018 Searching for sterile neutrinos in dynamical dark energy cosmologies Sci. China Phys. Mech. Astron. 61 050411

DOI

128
Feng L, Zhang J-F, Zhang X 2019 Search for sterile neutrinos in a universe of vacuum energy interacting with cold dark matter Phys. Dark Univ. 23 100261

DOI

129
Knee A M, Contreras D, Scott D 2019 Cosmological constraints on sterile neutrino oscillations from Planck J. Cosmol. Astropart. Phys. 2019 039

DOI

130
Feng L, He D-Z, Li H-L, Zhang J-F, Zhang X 2020 Constraints on active and sterile neutrinos in an interacting dark energy cosmology Sci. China Phys. Mech. Astron. 63 290404

DOI

131
Feng L, Guo R-Y, Zhang J-F, Zhang X 2022 Cosmological search for sterile neutrinos after Planck 2018 Phys. Lett. B 827 136940

DOI

132
Di Valentino E, Gariazzo S, Giunti C, Mena O, Pan S, Yang W 2022 Minimal dark energy: key to sterile neutrino and Hubble constant tensions? Phys. Rev. D 105 103511

DOI

133
Chernikov P A, Ivanchik A V 2022 The influence of the effective number of active and sterile neutrinos on the determination of the values of cosmological parameters Astron. Lett. 48 689 701

DOI

134
Pan S, Seto O, Takahashi T, Toda Y 2023 Constraints on sterile neutrinos and the cosmological tensions Phys. Rev. D 110 083524

DOI

135
Verde L, Treu T, Riess A G 2019 Tensions between the early and the late Universe Nature Astron. 3 891 7

DOI

136
Vitale S, Farr W M, Ng K, Rodriguez C L 2019 Measuring the star formation rate with gravitational waves from binary black holes Astrophys. J. Lett. 886 L1

DOI

137
Yang T 2021 Gravitational-wave detector networks: standard sirens on cosmology and modified gravity theory J. Cosmol. Astropart. Phys. 05 044

DOI

138
Belgacem E, Dirian Y, Foffa S, Howell E J, Maggiore M, Regimbau T 2019 Cosmology and dark energy from joint gravitational wave-GRB observations J. Cosmol. Astropart. Phys. 2019 015

DOI

139
Chen Z-C, Huang F, Huang Q-G 2019 Stochastic gravitational-wave background from binary black holes and binary neutron stars and implications for LISA Astrophys. J. 871 97

DOI

140
Du M, Xu L 2022 How will our knowledge of short gamma-ray bursts affect the distance measurement of binary neutron stars? Sci. China Phys. Mech. Astron. 65 219811

DOI

141
de Souza J M S, Sturani R 2021 Cosmological model selection from standard siren detections by third-generation gravitational wave observatories Phys. Dark Univ. 32 100830

DOI

142
Regimbau T, Evans M, Christensen N, Katsavounidis E, Sathyaprakash B, Vitale S 2017 Digging deeper: Observing primordial gravitational waves below the binary black hole produced stochastic background Phys. Rev. Lett. 118 151105

DOI

143
Belgacem E, Dirian Y, Foffa S, Maggiore M 2018 Modified gravitational-wave propagation and standard sirens Phys. Rev. D 98 023510

DOI

144
Safarzadeh M, Berger E, Ng K K Y, Chen H-Y, Vitale S, Whittle C, Scannapieco E 2019 Measuring the delay time distribution of binary neutron stars. II. Using the redshift distribution from third-generation gravitational wave detectors network Astrophys. J. Lett. 878 L13

DOI

145
Song H-R, Ai S-K, Wang M-H, Xing N, Gao H, Zhang B 2019 Viewing angle constraints on S190425z and S190426c and the joint gravitational-wave/gamma-ray detection fractions for binary neutron star mergers Astrophys. J. Lett. 881 L40

DOI

146
Wanderman D, Piran T 2015 The rate, luminosity function and time delay of non-collapsar short GRBs Mon. Not. Roy. Astron. Soc. 448 3026 3037

DOI

147
Yu J, Song H, Ai S, Gao H, Wang F, Wang Y, Lu Y, Fang W, Zhao W 2021 Multimessenger detection rates and distributions of binary neutron star mergers and their cosmological implications Astrophys. J. 916 54

DOI

148
Regimbau T, Siellez K, Meacher D, Gendre B, Boër M 2015 Revisiting coincidence rate between gravitational wave detection and short gamma-ray burst for the advanced and third generation Astrophys. J. 799 69

DOI

149
Madau P, Dickinson M 2014 Cosmic star formation history Ann. Rev. Astron. Astrophys. 52 415 486

DOI

150
Eichhorn A, Koslowski T, Pereira A D 2019 Status of background-independent coarse-graining in tensor models for quantum gravity Universe 5 53

DOI

151
Abbott R et al 2023 Population of merging compact binaries inferred using gravitational waves through GWTC-3 Phys. Rev. X 13 011048

DOI

152
Abbott B P et al 2019 GWTC-1: A gravitational-wave transient catalog of compact binary mergers observed by LIGO and virgo during the first and second observing runs Phys. Rev. X 9 031040

DOI

153
Özel F, Freire P 2016 Masses, radii, and the equation of state of neutron stars Ann. Rev. Astron. Astrophys. 54 401 440

DOI

154
Zhang X, Liu T, Zhao W 2017 Gravitational radiation from compact binary systems in screened modified gravity Phys. Rev. D 95 104027

DOI

155
Zhao W, Wen L 2018 Localization accuracy of compact binary coalescences detected by the third-generation gravitational-wave detectors and implication for cosmology Phys. Rev. D 97 064031

DOI

156
Wen L, Chen Y 2010 Geometrical expression for the angular resolution of a network of gravitational-wave detectors Phys. Rev. D 81 082001

DOI

157
Curt Cutler et al 1993 The last three minutes: issues in gravitational wave measurements of coalescing compact binaries Phys. Rev. Lett. 70 2984 2987

DOI

158
Sathyaprakash B S, Schutz B F 2009 Physics, astrophysics and cosmology with gravitational waves Living Rev. Rel. 12 2

DOI

159
Howell E J, Ackley K, Rowlinson A, Coward D 2018 Joint gravitational wave—gamma-ray burst detection rates in the aftermath of GW170817 Mon. Not. Royal Astron. Soc. 485 1435 1447

DOI

160
Tan W-W, Yu Y-W 2020 The jet structure and the intrinsic luminosity function of short gamma-ray bursts Astrophys. J. 902 83

DOI

161
Stratta G, Amati L, Ciolfi R, Vinciguerra S 2018 THESEUS in the era of multi-messenger astronomy arXiv:1802.01677

162
Speri L, Tamanini N, Caldwell R R, Gair J R, Wang B 2021 Testing the quasar Hubble diagram with LISA standard sirens Phys. Rev. D 103 083526

DOI

163
Hirata C M, Holz D E, Cutler C 2010 Reducing the weak lensing noise for the gravitational wave Hubble diagram using the non-Gaussianity of the magnification distribution Phys. Rev. D 81 124046

DOI

164
Kocsis B, Frei Z, Haiman Z, Menou K 2006 Finding the electromagnetic counterparts of cosmological standard sirens Astrophys. J. 637 27 37

DOI

165
Aghanim N et al 2020 Planck 2018 results. VI. Cosmological parameters Astron. Astrophys. 641 A6

[Erratum: Astron. Astrophys. 652, C4 (2021)]

DOI

166
Beutler F, Blake C, Colless M, Jones D H, Lister S-S, Campbell L, Parker Q, Saunders W, Watson F 2011 The 6dF galaxy survey: Baryon acoustic oscillations and the local hubble constant Mon. Not. Roy. Astron. Soc. 416 3017 3032

DOI

167
Ross A J, Samushia L, Howlett C, Percival W J, Burden A, Manera M 2015 The clustering of the SDSS DR7 main Galaxy sample - I. A 4 per cent distance measure at z = 0: 15 Mon. Not. Roy. Astron. Soc. 449 835 847

DOI

168
Alam S et al 2017 The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological analysis of the DR12 galaxy sample Mon. Not. Roy. Astron. Soc. 470 2617 2652

DOI

169
Scolnic D M et al 2018 The complete light-curve sample of spectroscopically confirmed SNe Ia from Pan-STARRS1 and cosmological constraints from the combined pantheon sample Astrophys. J. 859 101

DOI

170
Lewis A, Bridle S 2002 Cosmological parameters from CMB and other data: A Monte Carlo approach Phys. Rev. D 66 103511

DOI

171

172

173
Zhu J-P et al 2023 Kilonovae and optical sfterglows from binary neutron star mergers. II. optimal search strategy for serendipitous observations and target-of-opportunity observations of gravitational wave triggers Astrophys. J. 942 88

DOI

Outlines

/