\jyear

2021

[1]\fnmSheean \surJolicoeur

[1]\orgdivPhysics & Astronomy, \orgnameUniversity of the Western Cape, \orgaddress\stateCape Town, \postcode7535, \countrySouth Africa

2]\orgdivInstitute of Cosmology & Gravitation, \orgnameUniversity of Portsmouth, \orgaddress\statePortsmouth, \postcodePO1 3FX, \countryUnited Kingdom

3]\orgdivNational Institute for Theoretical & Computational Sciences (NITheCS), \orgaddress\stateCape Town, \postcode7535, \countrySouth Africa

Constraining primordial non-Gaussianity by combining next-generation galaxy and 21 cm intensity mapping surveys

[email protected] \fnmRoy \surMaartens [email protected] \fnmSimthembile \surDlamini [email protected] * [ [

Abstract

Surveys of the matter distribution contain ‘fossil’ information on possible non-Gaussianity that is generated in the primordial Universe. This primordial signal survives only on the largest scales where cosmic variance is strongest. By combining different surveys in a multi-tracer approach, we can suppress the cosmic variance and signficantly improve the precision on the level of primordial non-Gaussianity. We consider a combination of an optical galaxy survey, like the recently initiated DESI survey, together with a new and very different type of survey, a 21 cm intensity mapping survey, like the upcoming SKAO survey. A Fisher forecast of the precision on the local primordial non-Gaussianity parameter $f_{\mathrm{NL}}$ , shows that this multi-tracer combination, together with non-overlap single-tracer information, can deliver precision comparable to that from the CMB. Taking account of the largest systematic, i.e. foreground contamination in intensity mapping, we find that $\sigma(f_{\mathrm{NL}})\sim 4$ .

1 Introduction

Constraining the local primordial non-Gaussianity (PNG) parameter $f_{\mathrm{NL}}$ presents an important challenge for next-generation large-scale structure surveys. Local-type PNG affects the power spectrum of dark matter tracers by inducing a scale-dependent correction to the tracer bias (Matarrese:2008nc, ; Dalal:2007cu, ). This correction is strongly suppressed on small to medium scales, but on very large scales, it is $\propto f_{\mathrm{NL}}(H_{0}^{2}/k^{2})$ . This means that ultra-large survey volumes are required for high-precision constraints.

The Planck survey of the cosmic microwave background (CMB) gives the current state-of-the-art constraint $\sigma(f_{\mathrm{NL}})=5.1$ from the bispectrum (Akrami:2019izv, ). Future CMB surveys will lead to an improvement on the Planck precision, but not by much and not sufficient to access $\sigma(f_{\mathrm{NL}})\lesssim 1$ . This level of precision would enable us to rule out many single-field and multi-field inflationary scenarios (Alvarez:2014vva, ; dePutter:2016trg, ; Meerburg:2019qqi, ). The hope is that this goal can be achieved by using the extra modes in 3-dimensional surveys of the large-scale structure (see e.g. (Giannantonio:2011ya, ; Camera:2013kpa, ; Font-Ribera:2013rwa, ; Camera:2014bwa, ; Alonso:2015uua, ; Raccanelli:2015vla, )). However, extracting $f_{\mathrm{NL}}$ from these surveys faces formidable difficulties:

•

Observational systematics on very large scales Leistedt:2014zqa ; Fonseca:2020lmi ; Cunnington:2020wdu ; Rezaie:2021voi ; Liu:2020izx ; Spinelli:2021emp ; Riquelme:2022amo . For our simplified Fisher analysis, we neglect these systematics, except for the foreground contamination of 21 cm intensity mapping (see below).
•

Astrophysical complications entailed in the modelling of tracer clustering, in particular, the uncertainties in modelling the scale-dependent bias Barreira:2022sey ; Lazeyras:2022koc . We use a simplified model since a full treatment requires simulations, beyond the scope of this paper.
•

A theoretical systematic arising from the neglect of lensing magnification and other relativistic lightcone effects on the power spectrum, which can bias the best-fit $f_{\mathrm{NL}}$ Namikawa:2011yr ; Bruni:2011ta ; Jeong:2011as ; Camera:2014sba ; Kehagias:2015tda ; Lorenz:2017iez ; Contreras:2019bxy ; Jelic-Cizmek:2020pkh ; Bernal:2020pwq ; Wang:2020ibf ; Maartens:2020jzf ; Castorina:2021xzs ; Martinelli:2021ahc ; Viljoen:2021ocx . Although the bias on best-fit can be large, the error $\sigma(f_{\mathrm{NL}})$ is not significantly affected and we therefore omit these lightcone effects.
•

Cosmic (or sample) variance that dominates ultra-large scale measurements. We deal with this problem via a multi-tracer approach – which also helps to mitigate uncorrelated systematics.

Recent constraints on $f_{\mathrm{NL}}$ have used the BOSS galaxy Cabass:2022ymb and eBOSS quasar (Castorina:2019wmr, ) samples, with the tightest constraint to date of $\sigma(f_{\mathrm{NL}})=21$ from eBOSS (Mueller:2021tqa, ). The much larger volumes of upcoming surveys should facilitate significant improvement on this precision. But if we rely on individual surveys using the power spectrum of single tracers, it turns out that $\sigma(f_{\mathrm{NL}})\lesssim 1$ is not achievable – even if we neglect systematics and complexities in scale-dependent bias Alonso:2015uua . The fundamental problem is cosmic variance.

A way to evade cosmic variance is the multi-tracer technique (Seljak:2008xr, ; McDonald:2008sh, ; Bernstein:2011ju, ; Hamaus:2011dq, ; Abramo:2013awa, ; Abramo:2015iga, ; Witzemann:2018cdx, ; Wang:2020dtd, ; Abramo:2022qir, ; Alonso:2015sfa, ; Fonseca:2015laa, ). Forecasts indicate that multi-tracing next-generation surveys can surpass the CMB precision, and in some cases can achieve $\sigma(f_{\mathrm{NL}})\lesssim 1$ (Ferramacho:2014pua, ; Yamauchi:2014ioa, ; dePutter:2014lna, ; Alonso:2015sfa, ; Fonseca:2015laa, ; Fonseca:2016xvi, ; Schmittfull:2017ffw, ; Munchmeyer:2018eey, ; Fonseca:2018hsu, ; SKA:2018ckk, ; Gomes:2019ejy, ; Ballardini:2019wxj, ; Bermejo-Climent:2021jxf, ; Viljoen:2021ypp, ). The performance of the multi-tracer improves when the difference in tracer properties, especially the Gaussian clustering biases, is significant. It also helps to suppress uncorrelated systematics.

In this paper, our aim is not to identify survey combinations that can achieve $\sigma(f_{\mathrm{NL}})\lesssim 1$ – for this purpose, one typically requires the very high number densities of photometric surveys dePutter:2014lna ; Alonso:2015sfa . Instead, our aim is to investigate the improvements over single-tracer constraints when using a new type of spectroscopic large-scale structure survey – neutral hydrogen (HI) intensity mapping of the 21 cm emission line – in combination with spectroscopic galaxy surveys. Such a pair of spectroscopic samples has very different clustering biases and systematics, which accentuates the advantages of a multi-tracer analysis.

We use a simple Fisher forecast on mock surveys, which are based on next-generation surveys that are starting to observe or close to starting. The two surveys we have in mind are: the Dark Energy Spectroscopic Instrument (DESI), with the Bright Galaxy Sample (BGS) and Emission Line Galaxy (ELG) samples (DESI:2016fyo, ; Yahia-Cherif:2020knp, ), and the Square Kilometre Array Observatory (SKAO), with intensity mapping surveys in Bands 1 and 2 (SKA:2018ckk, ; Berti:2022ilk, ). Galaxy number count surveys are well established, but a cosmological HI intensity mapping survey has not yet been implemented. However, pilot surveys on the SKAO precursor telescope MeerKAT have already been used to:
(a) measure the cross-power between MeerKAT and WiggleZ galaxy surveys Cunnington:2022uzo , and
(b) start developing a pipeline for the planned SKAO surveys Santos:2017qgq ; Wang:2020lkn .

We combine pairs of these surveys at low and high redshifts, using the Fisher single-tracer constraints in the non-overlap volume and the multi-tracer constraints in the overlap volume. The results show a reasonable improvement for the low-redshift pair, and a significant improvement for the high-redshift pair, compared to the standard single-tracer forecasts for $\sigma(f_{\mathrm{NL}})$ . The combination of all low- and high-redshift surveys improve on Planck, with $\sigma(f_{\mathrm{NL}})\sim 4$ . This constraint is based on avoiding foreground contamination of HI intensity mapping on very large scales. Foreground contamination of HI intensity mapping has a strong effect on the precision: if we neglect this foreground contamination, we find that $\sigma(f_{\mathrm{NL}})\sim 1$ .

2 Multi-tracer power spectra

At first order in perturbations, the density (or temperature) contrast of a tracer $A$ is related to the matter density contrast $\delta$ by

\Delta_{A}(z,\bm{k})=\Big{[}b_{A}(z)+f(z)\mu^{2}\Big{]}\delta(z,\bm{k})\,,

(1)

where $b_{A}$ is the (Gaussian) clustering bias, $\mu=\bm{\hat{k}}\cdot\bm{{n}}$ is the projection along the line-of-sight direction $\bm{n}$ , and $f$ is the linear growth rate. We assume $b_{A}$ has a known redshift evolution $\beta_{A}$ , following Agarwal:2020lov , so that

b_{A}(z)=b_{A0}\,\beta_{A}(z)\,,

(2)

where the amplitude $b_{A0}$ is a free parameter. We highlight the point that a multi-tracer approach reduces the impact of this assumption on $\sigma(f_{\mathrm{NL}})$ Viljoen:2021ocx .

The Fourier power spectra at tree-level are given by

\big{\langle}\Delta_{A}(z,\bm{k})\Delta_{B}(z,\bm{k}^{\prime})\big{\rangle}=(2\pi)^{3}P_{AB}(z,\bm{k})\delta^{\mathrm{D}}(\bm{k}+\bm{k}^{\prime})\,,

(3)

where the linear matter power spectrum (computed using CLASS (Blas:2011rf, )) is

P_{AB}=\big{(}b_{A}+f\mu^{2}\big{)}\big{(}b_{B}+f\mu^{2}\big{)}P.

(4)

In the presence of local PNG, the bias acquires a scale-dependent correction Matarrese:2008nc ; Dalal:2007cu ; Lazeyras:2022koc :

\hat{b}_{A}(z,k)=b_{A}(z)+b_{A\phi}(z)\frac{f_{\rm NL}}{\mathcal{M}(z,k)}\,,\quad\mathcal{M}(z,k)=\frac{2}{3\Omega_{m0}H_{0}^{2}}\frac{D(z)}{g_{\rm in}}T(k)\,k^{2}\,.

(5)

Here $T(k)$ is the matter transfer function (normalized to 1 on very large scales), $D$ is the growth factor (normalized to 1 today) and $g_{\rm in}$ is the initial growth suppression function, defined deep in the matter era. For $\Lambda$ CDM Villa:2015ppa

g_{\rm in}=\frac{3}{5}\big{(}1+z\big{)}D(z)\bigg{[}1+\frac{2f(z)}{3\Omega_{m}(z)}\bigg{]}\,.

(6)

The non-Gaussian bias factor $b_{A\phi}$ is halo-dependent and should be determined with the aid of simulations (Barreira:2020kvh, ; Barreira:2021dpt, ; Barreira:2022sey, ). A simplified model reduces it to

b_{A\phi}(z)=2\delta_{\rm c}\big{[}b_{A}(z)-p_{A}\big{]}\,,

(7)

where the critical collapse density is given by $\delta_{\rm c}=1.686\,$ and $p_{A}$ are halo-dependent constants. In the simplest (universal) halo model, $p_{A}=1$ for all tracers $A$ . But this model is not consistent with simulations (Barreira:2020kvh, ; Barreira:2021dpt, ; Barreira:2022sey, ). An improved (but still over-simplified) model allows the constant $p_{A}$ to vary with tracer. For galaxies chosen by stellar mass, simulations indicate that the rough approximation (Barreira:2020kvh, )

p_{g}\approx 0.55\,,

(8)

is an improvement, and we use this for the DESI-like samples. For HI intensity mapping, we use the approximation (Barreira:2021dpt, )

p_{H}\approx 1.25\,.

(9)

The $f_{\mathrm{NL}}$ terms in the power spectra are of order $H_{0}^{2}/k^{2}$ on ultra-large scales, $k\lesssim k_{\rm eq}$ , where $T\approx 1$ . Local PNG therefore dominates the power on ultra-large scales – but these are also the scales where cosmic variance is largest. We deal with the cosmic variance via a multi-tracer analysis, which includes the information from all auto- and cross-power spectra of the two (or more) tracers (see section 4).

3 Modelling the surveys

We consider two mock spectroscopic surveys:

$A=g$ : galaxy survey, similar to DESI surveys DESI:2016fyo ; Yahia-Cherif:2020knp ;

$A=H$ : 21 cm HI intensity mapping (IM) survey, similar to SKAO surveys SKA:2018ckk ; Viljoen:2020efi ; Berti:2022ilk .

In both cases, we have a low- and a high-redshift survey. The sky and redshift coverage of the individual and overlapping surveys are given in Table 1.

Table 1: Sky area and redshift range of surveys.

Survey	Sample	$\Omega_{\mathrm{sky}}$	$t_{\mathrm{tot}}$	redshift
		$\big{[}10^{3}\,\mathrm{deg}^{2}\big{]}$	$\big{[}10^{3}\,\mathrm{hr}\big{]}$	range
$g$ (DESI-like)	BGS	14	-	0.00–0.50
	ELG	14	-	0.60–1.70
$H$ (SKAO-like)	Band 2	20	10	0.10–0.58
	Band 1	20	10	0.35–3.05
$g\times H$ (low $z$ )	BGS $\times$ Band 2	10	5	0.10–0.50
$g\times H$ (high $z$ )	ELG $\times$ Band 1	10	5	0.60–1.70

Refer to caption — Figure 1: Background comoving galaxy number densities (red, left-hand $y$ -axis) and HI IM temperature (blue, right-hand $y$ -axis).

The one-parameter model (2) for the Gaussian clustering bias of the galaxy surveys follows (Jelic-Cizmek:2020pkh, )

b_{g}(z)=\frac{b_{g0}}{D(z)}\,.

(10)

Fiducial values are $b_{g0}=1.34$ for the Bright Galaxy Sample (BGS) and $b_{g0}=0.84$ for the Emission Line Galaxy (ELG) sample. For the HI IM surveys we use a fit based on (Villaescusa-Navarro:2018vsg, ; Cunnington:2022ryj, ):

b_{H}(z)={b_{H0}}{\big{(}1+0.823\,z-0.0546\,z^{2}\big{)}}\quad\mbox{with fiducial value}~{}b_{H0}=0.842\,.

(11)

The background comoving number density $\bar{n}_{g}$ for the BGS and ELG samples is modelled following (Jelic-Cizmek:2020pkh, ), which leads to the smoothed curves shown in Figure 1. For HI IM, the background brightness temperature is modelled via the fit given in (Santos:2017qgq, ; Fonseca:2019qek, ):

\bar{T}_{H}(z)=0.0559+0.2324\,z-0.0241\,z^{2}~{}~{}\mathrm{mK}\,,

(12)

which is also shown in Figure 1.

3.1 Noise

In galaxy surveys, the shot noise (assumed to be Poissonian) is

P^{\rm shot}_{gg}(z)=\frac{1}{\bar{n}_{g}(z)}\,.

(13)

Then the total signal that we measure for the galaxy auto-power spectrum is

\tilde{P}_{gg}(z,k,\mu)=P_{gg}(z,k,\mu)+P^{\rm shot}_{gg}(z)\,,

(14)

where $P_{gg}$ is given by (4).

In IM surveys, there is shot noise, but also thermal (or instrumental) noise. The Poissonian shot noise (see Umeh:2021xqm for non-Poissonian corrections) is derived in a halo-model framework and given by Castorina:2016bfm ; Villaescusa-Navarro:2018vsg ,

P_{H\!H}^{\rm shot}(z)=\frac{1}{\bar{n}_{H}(z)}=\frac{1}{\bar{\rho}_{H}^{\,2}(z)}\int\mathrm{d}M\,{\mathcal{N}}_{h}(M,z)\,M_{H}^{2}(M,z)\,,

(15)

where $\bar{n}_{H}$ is the effective average comoving number density of HI, $\bar{\rho}_{H}$ is the average comoving HI density, ${\cal N}_{h}$ is the halo mass function (average comoving halo number density per mass) and $M_{H}$ is the average HI mass in a halo of mass $M$ . However, on the linear scales that we consider, the shot noise is much smaller than the thermal noise and can be safely neglected Castorina:2016bfm ; Villaescusa-Navarro:2018vsg .

The thermal noise in HI IM depends on the sky temperature in the radio band, the survey specifications and the array configuration (single-dish or interferometer). For the single-dish mode of SKAO-like surveys, the thermal noise power spectrum is Bull:2014rha ; Alonso:2017dgh ; Jolicoeur:2020eup :

P_{H\!H}^{\mathrm{therm}}(z)=\frac{\Omega_{\mathrm{sky}}}{2\nu_{21}t_{\mathrm{tot}}}\,\frac{(1+z)r(z)^{2}}{\mathcal{H}(z)}\,\left[\frac{T_{\mathrm{sys}}(z)}{\bar{T}_{H}(z)}\right]^{2}\,\frac{1}{N_{\rm d}}\,,

(16)

where $\nu_{21}=1420\,\mathrm{MHz}$ is the rest-frame frequency of the 21 cm emission, $t_{\mathrm{tot}}$ is the total observing time, and the number of dishes is $N_{\mathrm{d}}=197$ (with dish diameter $D_{\mathrm{d}}=15\,$ m). The system temperature is modelled as (Ansari:2018ury, ):

T_{\mathrm{sys}}(z)=T_{\rm d}(z)+T_{\rm sky}(z)=T_{\rm d}(z)+2.7+25\bigg{[}\frac{400\,\mathrm{MHz}}{\nu_{21}}(1+z)\bigg{]}^{2.75}~{}\mathrm{K},

(17)

where $T_{\rm d}$ is the dish receiver temperature (see Viljoen:2020efi ). The total signal is then

\tilde{P}_{H\!H}(z,k,\mu)=P_{H\!H}(z,k,\mu)+P^{\rm shot}_{H\!H}(z)+P_{H\!H}^{\mathrm{therm}}(z)\approx P_{H\!H}(z,k,\mu)+P_{H\!H}^{\mathrm{therm}}(z)\,.

(18)

The noise for all surveys is shown in Figure 2.

In the case of the cross-power spectrum $P_{gH}$ , cross-shot noise arises if the halos that host galaxies and HI overlap. Following the arguments given in Viljoen:2020efi ; Casas:2022vik , we assume that the cross-shot noise may be neglected. The total cross-power signal is then

\tilde{P}_{gH}=P_{gH}~{}~{}\mbox{with}~{}~{}P_{gH}^{\rm shot}\approx 0\,.

(19)

3.2 Intensity mapping beam and foregrounds

HI IM surveys like the SKAO surveys have poor angular resolution, which results in power loss on small transverse scales, i.e. for large $k_{\perp}=(1-\mu^{2})^{1/2}k$ . This effect is typically modelled by a Gaussian beam factor (Bull:2014rha, ):

{\cal D}_{\rm b}(z,k,\mu)=\exp\left[-{(1-\mu^{2})k^{2}r(z)^{2}\theta_{\rm b}(z)^{2}\over 16\ln 2}\right]\quad\mbox{with}\quad\theta_{\rm b}(z)=1.22\,{\lambda_{21}(1+z)\over D_{\rm d}}\,.

(20)

HI IM surveys are also contaminated by foregrounds that are much larger than the HI signal. Since these foregrounds are spectrally smooth, they can be separated from the non-smooth signal on small to medium scales. However, on very large radial scales, i.e. for small $k_{\|}=\mu k$ , the signal becomes smoother and therefore the separation fails. A comprehensive treatment includes simulations of foreground cleaning of the HI signal (e.g. Cunnington:2020wdu ; Spinelli:2021emp ). For a simplified Fisher forecast we can instead use a foreground avoidance approach, by excising the regions of Fourier space where the foregrounds are significant. This means removing large radial scales, which can be modelled via an exponential suppression factor (Bernal:2019jdo, ; Cunnington:2022ryj, ):

{\cal D}_{\rm fg}(k,\mu)={1-\exp\bigg{[}-\Big{(}\frac{\mu k}{k_{\rm fg}}\Big{)}^{2}\bigg{]}\,.}

(21)

We choose a typically used value:

k_{\rm fg}=0.01\,h\,{\rm Mpc}^{-1}\,.

(22)

Figure 3 illustrates the consequences of foreground avoidance for the monopoles of the HI power spectra,

P_{AB,0}\,(z,k)=\frac{1}{2}\int_{-1}^{+1}\mathrm{d}\mu\,P_{AB}\,(z,k,\mu)\;.

(23)

Note the consequence of foreground avoidance for the monopole of the HI power spectra $P_{AH,0}$ , which are suppressed on large scales. It is clear that large-scale (small $k$ ) power in $P_{AH,0}$ is lost and that the effect of local PNG on these large scales is suppressed by foreground avoidance. Note that $f_{\rm NL}$ reduces the HI power on very large scales at $z=0.3$ , since $b_{H}<1.25$ , which leads to $b_{H\phi}<0$ in (7). Also note that the beam effect on $P_{AH,0}$ , which tends to suppress small scale (large $k$ ) modes, is negligible at the low redshift, but becomes apparent at higher redshift.

In summary, the effects of the radio telescope beam and radio foreground avoidance lead to the following modifications of the power spectra $P_{AH}$ :

	$\displaystyle P_{H\!H}(z,{k},\mu)$	$\displaystyle\rightarrow$	$\displaystyle{\cal D}_{\rm fg}(k,\mu)\,{{\cal D}_{\rm b}(z,k,\mu)}^{2}\,P_{H\!H}(z,{k},\mu)\,,$		(24)
	$\displaystyle P_{gH}(z,{k},\mu)$	$\displaystyle\rightarrow$	$\displaystyle{\cal D}_{\rm fg}(k,\mu)\,{\cal D}_{\rm b}(z,k,\mu)\,P_{gH}(z,{k},\mu)\,.$		(25)

We note that these equations are consistent with the recent paper Cunnington:2023jpq on foreground cleaning via a transfer function, which shows that the effect of foregrounds is the same in the HI auto-power as in the cross-correlation power spectra with galaxies. However, we should also point out that the modelling of foreground removal is an ongoing project which has not yet reached maturity.

4 Fisher forecast

For a multi-tracer combination of two dark matter tracers $g$ and $H$ , the data vector of power spectra is

\bm{P}=\big{(}P_{gg}\,,\,P_{gH}\,,\,P_{H\!H}\big{)}\,.

(26)

We exclude the noise from the data vector, since the noise is independent of the cosmological and nuisance parameters. (The noise appears in the covariance, as given below.) The standard cosmological parameters are measured on medium to small scales and are effectively independent of the ultra-large scale parameter $f_{\mathrm{NL}}$ . Fixing their values could bias the best-fit value of $f_{\mathrm{NL}}$ but is unlikely to have any significant impact on $\sigma(f_{\mathrm{NL}})$ . We include in the Fisher analysis the cosmological parameters that directly affect the large-scale amplitude ( $\sigma_{8,0}$ ) and shape ( $n_{s}$ ) of the power spectrum. We therefore consider the following cosmological and nuisance parameters:

\vartheta_{\alpha}=\big{(}\sigma_{8,0},\,n_{s},\,f_{\mathrm{NL}}\mathchar 24635\relax\;\,b_{g0},\,b_{H0}\big{)}\,.

(27)

Here we assume that the degeneracies between $b_{A}$ and $\sigma_{8,0}$ , and between $\bar{T}_{H}$ and $\sigma_{8,0}$ , have been broken by other surveys that focus on medium to small scales.

The covariance for the multi-tracer power spectra is given by Barreira:2020ekm ; Karagiannis:2023 :

\mathrm{Cov}(\bm{P},\bm{P})=\frac{k_{\rm f}^{3}}{2\pi k^{2}\Delta k}\,\frac{2}{\Delta\mu}\,\begin{pmatrix}\tilde{P}_{gg}^{2}&&\tilde{P}_{gg}\tilde{P}_{gH}&&\tilde{P}_{gH}^{2}\\ &&&&\\ \tilde{P}_{gg}\tilde{P}_{gH}&&\frac{1}{2}\big{[}\tilde{P}_{gg}\tilde{P}_{H\!H}+\tilde{P}_{gH}^{2}\big{]}&&\tilde{P}_{H\!H}\tilde{P}_{gH}\\ &&&&\\ \tilde{P}_{gH}^{2}&&\tilde{P}_{H\!H}\tilde{P}_{gH}&&\tilde{P}_{H\!H}^{2}\end{pmatrix}\,.

(28)

Here $\Delta k$ and $\Delta\mu$ are the bin-widths for $k$ and $\mu$ and $k_{\mathrm{f}}$ is the fundamental mode, corresponding to the longest wavelength, which is determined by the comoving survey volume of the redshift bin centred at $z$ :

V(z)=\frac{\Omega_{\mathrm{sky}}}{3}\bigg{[}r\Big{(}z+\frac{\Delta z}{2}\Big{)}^{3}-r\Big{(}z-\frac{\Delta z}{2}\Big{)}^{3}\bigg{]}=\bigg{[}\frac{2\pi}{k_{\mathrm{f}}(z)}\bigg{]}^{3}\,.

(29)

Then the multi-tracer Fisher matrix in a redshift bin is

F_{\alpha\beta}^{\bm{P}}=\sum_{\mu=-1}^{+1}\,\sum_{k=k_{\mathrm{min}}}^{k_{\mathrm{max}}}\,\partial_{\alpha}\,\bm{P}\cdot\mathrm{Cov}(\bm{P},\bm{P})^{-1}\cdot\partial_{\beta}\,\bm{P}^{\mathrm{T}}\,,

(30)

where $\partial_{\alpha}=\partial\,/\,\partial\vartheta_{\alpha}$ . We choose the bin-widths and $k_{\mathrm{min}}$ following (Karagiannis:2018jdt, ; Yankelevich:2018uaz, ; Maartens:2019yhx, ):

\Delta\mu=0.04\,,~{}~{}\Delta k=k_{\mathrm{f}}\,,~{}~{}k_{\mathrm{min}}=k_{\mathrm{f}}\,.

(31)

The smallest scale (largest $k$ ) is chosen to ensure that linear perturbations remain accurate, since we do not require information from nonlinear scales Smith:2002dz :

k_{\mathrm{max}}=0.08(1+z)^{2/(2+n_{s})}\,h\,\text{Mpc}^{-1}\,.

(32)

The galaxy and HI IM surveys do not have the same sky area and redshift ranges (see Table 1). The overlap redshift ranges are obvious; for the sky areas, we assume a nominal overlap area of 10,000 deg². The multi-tracer applies in the overlap redshift range and overlap sky area. However, we can add the independent Fisher matrices obtained in the non-overlapping regions of the two individual surveys to the multi-tracer Fisher matrix in the overlapping region (Viljoen:2020efi, ). In detail, the non-overlap region is

•

the non-overlap sky area of each survey, across the full redshift range of each survey;
•

the overlap sky area, across the non-overlap parts of the redshift ranges of each survey.

Then the full Fisher matrix (denoted by $g\otimes H$ ) is

F_{\alpha\beta}^{g\otimes H}=F_{\alpha\beta}^{\bm{P}}(\text{overlap})+F_{\alpha\beta}^{g}(\text{non-overlap})+F_{\alpha\beta}^{H}(\text{non-overlap}).

(33)

Figure 4 shows the 1 $\sigma$ error contours computed from this Fisher matrix after marginalising over the bias parameters in (27). We use the fiducial values $\sigma_{8,0}=0.8102$ , $n_{s}=0.9665$ and $f_{\mathrm{NL}}=0$ , keeping other cosmological parameters fixed to their Planck 2018 best-fit values (Planck:2018vyg, ). Table 2 lists the marginalised $\sigma(f_{\mathrm{NL}})$ constraints, including the improvements (in parentheses) that follow when we ignore the HI IM foregrounds, i.e. when $k_{\rm fg}=0$ .

The results for BGS, Band 1 and Band 2 are consistent with Viljoen:2021ocx , which uses the angular power spectrum rather than the Fourier power spectrum used here.

Table 2: Cumulative marginalised error

\sigma(f_{\mathrm{NL}})

for: the single-tracer surveys; the low-

z

and high-

z

multi-tracer pairs – including single-tracer information from the non-overlap volumes as in (33); and the sum of the multi-tracer pairs.
Values in parenthesis correspond to ignoring foregrounds, i.e., with

k_{\mathrm{fg}}=0

Survey	$\sigma(f_{\mathrm{NL}})$
BGS	29.4
ELG	4.9
Band 2	95.7 (66.9)
Band 1	6.1 (1.9)
BGS $\otimes$ Band 2	13.0 (3.2)
ELG $\otimes$ Band 1	4.4 (1.8)
BGS $\otimes$ Band 2 + ELG $\otimes$ Band 1	4.0 (1.5)

5 Conclusion

We applied a simplified Fisher forecast analysis in order to gain insights into the improvements on $\sigma(f_{\mathrm{NL}})$ from a new multi-tracer combination of large-scale structure surveys that will be made possible by the upcoming HI intensity mapping surveys with the SKAO. For the spectroscopic galaxy samples, we used mock surveys based on DESI. This allowed for multi-tracer pairs at low and high redshifts, in order to see the level of improved precision from the higher volume at high redshifts.

We included all available Fisher information from these mock surveys, using the multi-tracer Fisher matrix on the overlap volume (sky area and redshift range) and the single-tracer Fisher matrices on the non-overlap volumes, as in (33). Our results for $\sigma(f_{\mathrm{NL}})$ are summarised in Table 2, which gives the single-tracer results for the 4 samples, the full multi-tracer results from (33) for the low- and high-redshift pairs, and finally the results from the Fisher sum of these pairs. The impact of foreground avoidance can be seen from the results (in parenthesis) when foregrounds are ignored, corresponding to $k_{\rm fg}=0$ in (21).

Figure 4 shows the $1\sigma$ contours for the 3 cosmological parameters ( $\sigma_{8,\,0}$ , $n_{s}$ , $f_{\mathrm{NL}}$ ), marginalising over the 2 bias nuisance parameters ( $b_{g0}$ , $b_{H0}$ ).

As expected, the foreground avoidance severely degrades the constraining power of the HI intensity surveys. Without foregrounds, the high-redshift Band 1 survey would perform as well as the multi-tracer combination with the ELG survey: $\sigma(f_{\mathrm{NL}})=1.8$ . The reduced precision on $f_{\rm NL}$ for single-tracer HI intensity constraints is consistent with the findings of Cunnington:2020wdu , where foreground cleaning is simulated. However, using the multi-tracer with a spectroscopic galaxy survey helps to reduce the effect of foregrounds on $\sigma(f_{\mathrm{NL}})$ , as shown by Table 2.

The multi-tracer analysis is also expected to mitigate other systematics in the galaxy and HI intensity samples. For example, HI intensity mapping is also affected by radio frequency interference, receiver noise, calibration errors and polarisation leakage. Incorporating these and the galaxy systematics requires extensive simulations, beyond the scope of our paper.

Acknowledgements
We thank Dionysis Karagiannis for very helpful comments. We are supported by the South African Radio Astronomy Observatory (SARAO) and the National Research Foundation (Grant No. 75415).

References

\bibcommenthead
(1) Matarrese, S., Verde, L.: The effect of primordial non-Gaussianity on halo bias. Astrophys. J. Lett. 677, 77–80 (2008) arXiv:0801.4826 [astro-ph]. https://doi.org/10.1086/587840
(2) Dalal, N., Dore, O., Huterer, D., Shirokov, A.: The imprints of primordial non-gaussianities on large-scale structure: scale dependent bias and abundance of virialized objects. Phys. Rev. D77, 123514 (2008) arXiv:0710.4560 [astro-ph]. https://doi.org/10.1103/PhysRevD.77.123514
(3) Akrami, Y., et al.: Planck 2018 results. IX. Constraints on primordial non-Gaussianity. Astron. Astrophys. 641, 9 (2020) arXiv:1905.05697 [astro-ph.CO]. https://doi.org/10.1051/0004-6361/201935891
(4) Alvarez, M., et al.: Testing Inflation with Large Scale Structure: Connecting Hopes with Reality (2014) arXiv:1412.4671 [astro-ph.CO]
(5) de Putter, R., Gleyzes, J., Doré, O.: Next non-Gaussianity frontier: What can a measurement with $\sigma$ (fNL) $\lesssim$ 1 tell us about multifield inflation? Phys. Rev. D 95(12), 123507 (2017) arXiv:1612.05248 [astro-ph.CO]. https://doi.org/10.1103/PhysRevD.95.123507
(6) Meerburg, P.D., et al.: Primordial Non-Gaussianity (2019) arXiv:1903.04409 [astro-ph.CO]
(7) Giannantonio, T., Porciani, C., Carron, J., Amara, A., Pillepich, A.: Constraining primordial non-Gaussianity with future galaxy surveys. Mon. Not. Roy. Astron. Soc. 422, 2854–2877 (2012) arXiv:1109.0958 [astro-ph.CO]. https://doi.org/10.1111/j.1365-2966.2012.20604.x
(8) Camera, S., Santos, M.G., Ferreira, P.G., Ferramacho, L.: Cosmology on Ultra-Large Scales with HI Intensity Mapping: Limits on Primordial non-Gaussianity. Phys. Rev. Lett. 111, 171302 (2013) arXiv:1305.6928 [astro-ph.CO]. https://doi.org/10.1103/PhysRevLett.111.171302
(9) Font-Ribera, A., McDonald, P., Mostek, N., Reid, B.A., Seo, H.-J., Slosar, A.: DESI and other dark energy experiments in the era of neutrino mass measurements. JCAP 05, 023 (2014) arXiv:1308.4164 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2014/05/023
(10) Camera, S., Santos, M.G., Maartens, R.: Probing primordial non-Gaussianity with SKA galaxy redshift surveys: a fully relativistic analysis. Mon. Not. Roy. Astron. Soc. 448(2), 1035–1043 (2015) arXiv:1409.8286 [astro-ph.CO]. https://doi.org/10.1093/mnras/stv040
(11) Alonso, D., Bull, P., Ferreira, P.G., Maartens, R., Santos, M.: Ultra large-scale cosmology in next-generation experiments with single tracers. Astrophys. J. 814(2), 145 (2015) arXiv:1505.07596 [astro-ph.CO]. https://doi.org/10.1088/0004-637X/814/2/145
(12) Raccanelli, A., Montanari, F., Bertacca, D., Doré, O., Durrer, R.: Cosmological Measurements with General Relativistic Galaxy Correlations. JCAP 05, 009 (2016) arXiv:1505.06179 [astro-ph.CO]. https://doi.org/%****␣paper.bbl␣Line␣250␣****10.1088/1475-7516/2016/05/009
(13) Leistedt, B., Peiris, H.V., Roth, N.: Constraints on Primordial Non-Gaussianity from 800 000 Photometric Quasars. Phys. Rev. Lett. 113(22), 221301 (2014) arXiv:1405.4315 [astro-ph.CO]. https://doi.org/10.1103/PhysRevLett.113.221301
(14) Fonseca, J., Liguori, M.: Measuring ultralarge scale effects in the presence of 21 cm intensity mapping foregrounds. Mon. Not. Roy. Astron. Soc. 504(1), 267–279 (2021) arXiv:2011.11510 [astro-ph.CO]. https://doi.org/10.1093/mnras/stab903
(15) Cunnington, S., Camera, S., Pourtsidou, A.: The degeneracy between primordial non-Gaussianity and foregrounds in 21 cm intensity mapping experiments. Mon. Not. Roy. Astron. Soc. 499(3), 4054–4067 (2020) arXiv:2007.12126 [astro-ph.CO]. https://doi.org/10.1093/mnras/staa2986
(16) Rezaie, M., et al.: Primordial non-Gaussianity from the completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey – I: Catalogue preparation and systematic mitigation. Mon. Not. Roy. Astron. Soc. 506(3), 3439–3454 (2021) arXiv:2106.13724 [astro-ph.CO]. https://doi.org/10.1093/mnras/stab1730
(17) Liu, R.H., Breysse, P.C.: Coupling parsec and gigaparsec scales: Primordial non-Gaussianity with multitracer intensity mapping. Phys. Rev. D 103(6), 063520 (2021) arXiv:2002.10483 [astro-ph.CO]. https://doi.org/10.1103/PhysRevD.103.063520
(18) Spinelli, M., Carucci, I.P., Cunnington, S., Harper, S.E., Irfan, M.O., Fonseca, J., Pourtsidou, A., Wolz, L.: SKAO HI intensity mapping: blind foreground subtraction challenge. Mon. Not. Roy. Astron. Soc. 509(2), 2048–2074 (2021) arXiv:2107.10814 [astro-ph.CO]. https://doi.org/10.1093/mnras/stab3064
(19) Riquelme, W., et al.: Primordial non-Gaussianity with Angular correlation function: Integral constraint and validation for DES (2022) arXiv:2209.07187 [astro-ph.CO]
(20) Barreira, A.: Can we actually constrain f_NL using the scale-dependent bias effect? An illustration of the impact of galaxy bias uncertainties using the BOSS DR12 galaxy power spectrum. JCAP 11, 013 (2022) arXiv:2205.05673 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2022/11/013
(21) Lazeyras, T., Barreira, A., Schmidt, F., Desjacques, V.: Assembly bias in the local PNG halo bias and its implication for f _NL constraints. JCAP 01, 023 (2023) arXiv:2209.07251 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2023/01/023
(22) Namikawa, T., Okamura, T., Taruya, A.: Magnification effect on the detection of primordial non-Gaussianity from photometric surveys. Phys. Rev. D 83, 123514 (2011) arXiv:1103.1118 [astro-ph.CO]. https://doi.org/%****␣paper.bbl␣Line␣425␣****10.1103/PhysRevD.83.123514
(23) Bruni, M., Crittenden, R., Koyama, K., Maartens, R., Pitrou, C., Wands, D.: Disentangling non-Gaussianity, bias and GR effects in the galaxy distribution. Phys. Rev. D 85, 041301 (2012) arXiv:1106.3999 [astro-ph.CO]. https://doi.org/10.1103/PhysRevD.85.041301
(24) Jeong, D., Schmidt, F., Hirata, C.M.: Large-scale clustering of galaxies in general relativity. Phys. Rev. D85, 023504 (2012) arXiv:1107.5427 [astro-ph.CO]. https://doi.org/10.1103/PhysRevD.85.023504
(25) Camera, S., Maartens, R., Santos, M.G.: Einstein’s legacy in galaxy surveys. Mon. Not. Roy. Astron. Soc. 451(1), 80–84 (2015) arXiv:1412.4781 [astro-ph.CO]. https://doi.org/10.1093/mnrasl/slv069
(26) Kehagias, A., Dizgah, A.M., Norena, J., Perrier, H., Riotto, A.: A Consistency Relation for the Observed Galaxy Bispectrum and the Local non-Gaussianity from Relativistic Corrections. JCAP 1508(08), 018 (2015) arXiv:1503.04467 [astro-ph.CO]. https://doi.org/%****␣paper.bbl␣Line␣500␣****10.1088/1475-7516/2015/08/018
(27) Lorenz, C.S., Alonso, D., Ferreira, P.G.: Impact of relativistic effects on cosmological parameter estimation. Phys. Rev. D 97(2), 023537 (2018) arXiv:1710.02477 [astro-ph.CO]. https://doi.org/10.1103/PhysRevD.97.023537
(28) Contreras, D., Johnson, M.C., Mertens, J.B.: Towards detection of relativistic effects in galaxy number counts using kSZ Tomography. JCAP 10, 024 (2019) arXiv:1904.10033 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2019/10/024
(29) Jelic-Cizmek, G., Lepori, F., Bonvin, C., Durrer, R.: On the importance of lensing for galaxy clustering in photometric and spectroscopic surveys. JCAP 04, 055 (2021) arXiv:2004.12981 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2021/04/055
(30) Bernal, J.L., Bellomo, N., Raccanelli, A., Verde, L.: Beware of commonly used approximations. Part II. Estimating systematic biases in the best-fit parameters. JCAP 10, 017 (2020) arXiv:2005.09666 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2020/10/017
(31) Wang, M.S., Beutler, F., Bacon, D.: Impact of Relativistic Effects on the Primordial Non-Gaussianity Signature in the Large-Scale Clustering of Quasars. Mon. Not. Roy. Astron. Soc. 499(2), 2598–2607 (2020) arXiv:2007.01802 [astro-ph.CO]. https://doi.org/10.1093/mnras/staa2998
(32) Maartens, R., Jolicoeur, S., Umeh, O., De Weerd, E.M., Clarkson, C.: Local primordial non-Gaussianity in the relativistic galaxy bispectrum. JCAP 04, 013 (2021) arXiv:2011.13660 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2021/04/013
(33) Castorina, E., di Dio, E.: The observed galaxy power spectrum in General Relativity. JCAP 01(01), 061 (2022) arXiv:2106.08857 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2022/01/061
(34) Martinelli, M., Dalal, R., Majidi, F., Akrami, Y., Camera, S., Sellentin, E.: Ultralarge-scale approximations and galaxy clustering: Debiasing constraints on cosmological parameters. Mon. Not. Roy. Astron. Soc. 510(2), 1964–1977 (2022) arXiv:2106.15604 [astro-ph.CO]. https://doi.org/10.1093/mnras/stab3578
(35) Viljoen, J.-A., Fonseca, J., Maartens, R.: Multi-wavelength spectroscopic probes: biases from neglecting light-cone effects. JCAP 12(12), 004 (2021) arXiv:2108.05746 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2021/12/004
(36) Cabass, G., Ivanov, M.M., Philcox, O.H.E., Simonović, M., Zaldarriaga, M.: Constraints on multifield inflation from the BOSS galaxy survey. Phys. Rev. D 106(4), 043506 (2022) arXiv:2204.01781 [astro-ph.CO]. https://doi.org/10.1103/PhysRevD.106.043506
(37) Castorina, E., et al.: Redshift-weighted constraints on primordial non-Gaussianity from the clustering of the eBOSS DR14 quasars in Fourier space. JCAP 09, 010 (2019) arXiv:1904.08859 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2019/09/010
(38) Mueller, E.-M., et al.: The clustering of galaxies in the completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Primordial non-Gaussianity in Fourier Space (2021) arXiv:2106.13725 [astro-ph.CO]
(39) Seljak, U.: Extracting primordial non-gaussianity without cosmic variance. Phys. Rev. Lett. 102, 021302 (2009) arXiv:0807.1770 [astro-ph]. https://doi.org/10.1103/PhysRevLett.102.021302
(40) McDonald, P., Seljak, U.: How to measure redshift-space distortions without sample variance. JCAP 10, 007 (2009) arXiv:0810.0323 [astro-ph]. https://doi.org/10.1088/1475-7516/2009/10/007
(41) Bernstein, G.M., Cai, Y.-C.: Cosmology without cosmic variance. Mon. Not. Roy. Astron. Soc. 416, 3009 (2011) arXiv:1104.3862 [astro-ph.CO]. https://doi.org/10.1111/j.1365-2966.2011.19249.x
(42) Hamaus, N., Seljak, U., Desjacques, V.: Optimal Constraints on Local Primordial Non-Gaussianity from the Two-Point Statistics of Large-Scale Structure. Phys. Rev. D 84, 083509 (2011) arXiv:1104.2321 [astro-ph.CO]. https://doi.org/10.1103/PhysRevD.84.083509
(43) Abramo, L.R., Leonard, K.E.: Why multi-tracer surveys beat cosmic variance. Mon. Not. Roy. Astron. Soc. 432, 318 (2013) arXiv:1302.5444 [astro-ph.CO]. https://doi.org/10.1093/mnras/stt465
(44) Abramo, L.R., Secco, L.F., Loureiro, A.: Fourier analysis of multitracer cosmological surveys. Mon. Not. Roy. Astron. Soc. 455(4), 3871–3889 (2016) arXiv:1505.04106 [astro-ph.CO]. https://doi.org/10.1093/mnras/stv2588
(45) Witzemann, A., Alonso, D., Fonseca, J., Santos, M.G.: Simulated multitracer analyses with HI intensity mapping. Mon. Not. Roy. Astron. Soc. 485(4), 5519–5531 (2019) arXiv:1808.03093 [astro-ph.CO]. https://doi.org/10.1093/mnras/stz778
(46) Wang, Y., Zhao, G.-B.: A brief review on cosmological analysis of galaxy surveys with multiple tracers (2020) arXiv:2009.03862 [astro-ph.CO]. https://doi.org/10.1088/1674-4527/20/10/158
(47) Abramo, L.R., Dinarte Ferri, J.a.V., Tashiro, I.L., Loureiro, A.: Fisher matrix for the angular power spectrum of multi-tracer galaxy surveys. JCAP 08, 073 (2022) arXiv:2204.05057 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2022/08/073
(48) Alonso, D., Ferreira, P.G.: Constraining ultralarge-scale cosmology with multiple tracers in optical and radio surveys. Phys. Rev. D 92(6), 063525 (2015) arXiv:1507.03550 [astro-ph.CO]. https://doi.org/10.1103/PhysRevD.92.063525
(49) Fonseca, J., Camera, S., Santos, M., Maartens, R.: Hunting down horizon-scale effects with multi-wavelength surveys. Astrophys. J. 812(2), 22 (2015) arXiv:1507.04605 [astro-ph.CO]. https://doi.org/10.1088/2041-8205/812/2/L22
(50) Ferramacho, L.D., Santos, M.G., Jarvis, M.J., Camera, S.: Radio galaxy populations and the multitracer technique: pushing the limits on primordial non-Gaussianity. Mon. Not. Roy. Astron. Soc. 442(3), 2511–2518 (2014) arXiv:1402.2290 [astro-ph.CO]. https://doi.org/10.1093/mnras/stu1015
(51) Yamauchi, D., Takahashi, K., Oguri, M.: Constraining primordial non-Gaussianity via a multitracer technique with surveys by Euclid and the Square Kilometre Array. Phys. Rev. D 90(8), 083520 (2014) arXiv:1407.5453 [astro-ph.CO]. https://doi.org/10.1103/PhysRevD.90.083520
(52) de Putter, R., Doré, O.: Designing an Inflation Galaxy Survey: how to measure $\sigma(f_{\rm NL})\sim 1$ using scale-dependent galaxy bias. Phys. Rev. D 95(12), 123513 (2017) arXiv:1412.3854 [astro-ph.CO]. https://doi.org/10.1103/PhysRevD.95.123513
(53) Fonseca, J., Maartens, R., Santos, M.G.: Probing the primordial Universe with MeerKAT and DES. Mon. Not. Roy. Astron. Soc. 466(3), 2780–2786 (2017) arXiv:1611.01322 [astro-ph.CO]. https://doi.org/10.1093/mnras/stw3248
(54) Schmittfull, M., Seljak, U.: Parameter constraints from cross-correlation of CMB lensing with galaxy clustering. Phys. Rev. D 97(12), 123540 (2018) arXiv:1710.09465 [astro-ph.CO]. https://doi.org/10.1103/PhysRevD.97.123540
(55) Münchmeyer, M., Madhavacheril, M.S., Ferraro, S., Johnson, M.C., Smith, K.M.: Constraining local non-Gaussianities with kinetic Sunyaev-Zeldovich tomography. Phys. Rev. D 100(8), 083508 (2019) arXiv:1810.13424 [astro-ph.CO]. https://doi.org/10.1103/PhysRevD.100.083508
(56) Fonseca, J., Maartens, R., Santos, M.G.: Synergies between intensity maps of hydrogen lines. Mon. Not. Roy. Astron. Soc. 479(3), 3490–3497 (2018) arXiv:1803.07077 [astro-ph.CO]. https://doi.org/10.1093/mnras/sty1702
(57) Bacon, D.J., et al.: Cosmology with Phase 1 of the Square Kilometre Array: Red Book 2018: Technical specifications and performance forecasts. Publ. Astron. Soc. Austral. 37, 007 (2020) arXiv:1811.02743 [astro-ph.CO]. https://doi.org/10.1017/pasa.2019.51
(58) Gomes, Z., Camera, S., Jarvis, M.J., Hale, C., Fonseca, J.: Non-Gaussianity constraints using future radio continuum surveys and the multitracer technique. Mon. Not. Roy. Astron. Soc. 492(1), 1513–1522 (2020) arXiv:1912.08362 [astro-ph.CO]. https://doi.org/10.1093/mnras/stz3581
(59) Ballardini, M., Matthewson, W.L., Maartens, R.: Constraining primordial non-Gaussianity using two galaxy surveys and CMB lensing. Mon. Not. Roy. Astron. Soc. 489(2), 1950–1956 (2019) arXiv:1906.04730 [astro-ph.CO]. https://doi.org/10.1093/mnras/stz2258
(60) Bermejo-Climent, J.R., Ballardini, M., Finelli, F., Paoletti, D., Maartens, R., Rubiño-Martín, J.A., Valenziano, L.: Cosmological parameter forecasts by a joint 2D tomographic approach to CMB and galaxy clustering. Phys. Rev. D 103(10), 103502 (2021) arXiv:2106.05267 [astro-ph.CO]. https://doi.org/%****␣paper.bbl␣Line␣1100␣****10.1103/PhysRevD.103.103502
(61) Viljoen, J.-A., Fonseca, J., Maartens, R.: Multi-wavelength spectroscopic probes: prospects for primordial non-Gaussianity and relativistic effects. JCAP 11, 010 (2021) arXiv:2107.14057 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2021/11/010
(62) Aghamousa, A., et al.: The DESI Experiment Part I: Science,Targeting, and Survey Design (2016) arXiv:1611.00036 [astro-ph.IM]
(63) Yahia-Cherif, S., Blanchard, A., Camera, S., Casas, S., Ilić, S., Markovic, K., Pourtsidou, A., Sakr, Z., Sapone, D., Tutusaus, I.: Validating the Fisher approach for stage IV spectroscopic surveys. Astron. Astrophys. 649, 52 (2021) arXiv:2007.01812 [astro-ph.CO]. https://doi.org/10.1051/0004-6361/201937312
(64) Berti, M., Spinelli, M., Viel, M.: Multipole expansion for 21cm Intensity Mapping power spectrum: forecasted cosmological parameters estimation for the SKA Observatory (2022) arXiv:2209.07595 [astro-ph.CO]. https://doi.org/10.1093/mnras/stad685
(65) Cunnington, S., et al.: HI intensity mapping with MeerKAT: power spectrum detection in cross-correlation with WiggleZ galaxies (2022) arXiv:2206.01579 [astro-ph.CO]. https://doi.org/10.1093/mnras/stac3060
(66) Santos, M.G., et al.: MeerKLASS: MeerKAT Large Area Synoptic Survey. In: MeerKAT Science: On the Pathway to the SKA (2017)
(67) Wang, J., et al.: HI intensity mapping with MeerKAT: calibration pipeline for multidish autocorrelation observations. Mon. Not. Roy. Astron. Soc. 505(3), 3698–3721 (2021) arXiv:2011.13789 [astro-ph.CO]. https://doi.org/10.1093/mnras/stab1365
(68) Agarwal, N., Desjacques, V., Jeong, D., Schmidt, F.: Information content in the redshift-space galaxy power spectrum and bispectrum. JCAP 03, 021 (2021) arXiv:2007.04340 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2021/03/021
(69) Blas, D., Lesgourgues, J., Tram, T.: The Cosmic Linear Anisotropy Solving System (CLASS) II: Approximation schemes. JCAP 1107, 034 (2011) arXiv:1104.2933 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2011/07/034
(70) Villa, E., Rampf, C.: Relativistic perturbations in $\Lambda$ CDM: Eulerian & Lagrangian approaches. JCAP 1601(01), 030 (2016) arXiv:1505.04782 [gr-qc]. https://doi.org/10.1088/1475-7516/2016/01/030
(71) Barreira, A., Cabass, G., Schmidt, F., Pillepich, A., Nelson, D.: Galaxy bias and primordial non-Gaussianity: insights from galaxy formation simulations with IllustrisTNG. JCAP 12, 013 (2020) arXiv:2006.09368 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2020/12/013
(72) Barreira, A.: The local PNG bias of neutral Hydrogen, H ${I}$ . JCAP 04(04), 057 (2022) arXiv:2112.03253 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2022/04/057
(73) Viljoen, J.-A., Fonseca, J., Maartens, R.: Constraining the growth rate by combining multiple future surveys. JCAP 09, 054 (2020) arXiv:2007.04656 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2020/09/054
(74) Villaescusa-Navarro, F., et al.: Ingredients for 21 cm Intensity Mapping. Astrophys. J. 866(2), 135 (2018) arXiv:1804.09180 [astro-ph.CO]. https://doi.org/10.3847/1538-4357/aadba0
(75) Cunnington, S.: Detecting the power spectrum turnover with HI intensity mapping. Mon. Not. Roy. Astron. Soc. 512(2), 2408–2425 (2022) arXiv:2202.13828 [astro-ph.CO]. https://doi.org/10.1093/mnras/stac576
(76) Fonseca, J., Viljoen, J.-A., Maartens, R.: Constraints on the growth rate using the observed galaxy power spectrum. JCAP 1912(12), 028 (2019) arXiv:1907.02975 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2019/12/028
(77) Umeh, O., Maartens, R., Padmanabhan, H., Camera, S.: The effect of finite halo size on the clustering of neutral hydrogen. JCAP 06, 027 (2021) arXiv:2102.06116 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2021/06/027
(78) Castorina, E., Villaescusa-Navarro, F.: On the spatial distribution of neutral hydrogen in the Universe: bias and shot-noise of the HI power spectrum. Mon. Not. Roy. Astron. Soc. 471(2), 1788–1796 (2017) arXiv:1609.05157 [astro-ph.CO]. https://doi.org/10.1093/mnras/stx1599
(79) Bull, P., Ferreira, P.G., Patel, P., Santos, M.G.: Late-time cosmology with 21cm intensity mapping experiments. Astrophys. J. 803(1), 21 (2015) arXiv:1405.1452 [astro-ph.CO]. https://doi.org/10.1088/0004-637X/803/1/21
(80) Alonso, D., Ferreira, P.G., Jarvis, M.J., Moodley, K.: Calibrating photometric redshifts with intensity mapping observations. Phys. Rev. D 96(4), 043515 (2017) arXiv:1704.01941 [astro-ph.CO]. https://doi.org/10.1103/PhysRevD.96.043515
(81) Jolicoeur, S., Maartens, R., De Weerd, E.M., Umeh, O., Clarkson, C., Camera, S.: Detecting the relativistic bispectrum in 21cm intensity maps. JCAP 06, 039 (2021) arXiv:2009.06197 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2021/06/039
(82) Ansari, R., et al.: Inflation and Early Dark Energy with a Stage II Hydrogen Intensity Mapping experiment (2018) arXiv:1810.09572 [astro-ph.CO]
(83) Casas, S., Carucci, I.P., Pettorino, V., Camera, S., Martinelli, M.: Constraining gravity with synergies between radio and optical cosmological surveys. Phys. Dark Univ. 39, 101151 (2023) arXiv:2210.05705 [astro-ph.CO]. https://doi.org/10.1016/j.dark.2022.101151
(84) Bernal, J.L., Breysse, P.C., Gil-Marin, H., Kovetz, E.D.: User’s guide to extracting cosmological information from line-intensity maps. Phys. Rev. D 100(12), 123522 (2019) arXiv:1907.10067 [astro-ph.CO]. https://doi.org/10.1103/PhysRevD.100.123522
(85) Cunnington, S., et al.: The foreground transfer function for HI intensity mapping signal reconstruction: MeerKLASS and precision cosmology applications (2023) arXiv:2302.07034 [astro-ph.CO]
(86) Barreira, A.: On the impact of galaxy bias uncertainties on primordial non-Gaussianity constraints. JCAP 12, 031 (2020) arXiv:2009.06622 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2020/12/031
(87) Karagiannis, D., Maartens, R., Fonseca, J., Camera, S., Clarkson, C.: Multi-tracer power spectra and bispectra I: Formalism (2023) arXiv:2301.nnnnn [astro-ph.CO]
(88) Karagiannis, D., Lazanu, A., Liguori, M., Raccanelli, A., Bartolo, N., Verde, L.: Constraining primordial non-Gaussianity with bispectrum and power spectrum from upcoming optical and radio surveys. Mon. Not. Roy. Astron. Soc. 478(1), 1341–1376 (2018) arXiv:1801.09280 [astro-ph.CO]. https://doi.org/10.1093/mnras/sty1029
(89) Yankelevich, V., Porciani, C.: Cosmological information in the redshift-space bispectrum. Mon. Not. Roy. Astron. Soc. 483(2), 2078–2099 (2019) arXiv:1807.07076 [astro-ph.CO]. https://doi.org/10.1093/mnras/sty3143
(90) Maartens, R., Jolicoeur, S., Umeh, O., De Weerd, E.M., Clarkson, C., Camera, S.: Detecting the relativistic galaxy bispectrum. JCAP 03(03), 065 (2020) arXiv:1911.02398 [astro-ph.CO]. https://doi.org/10.1088/1475-7516/2020/03/065
(91) Smith, R.E., Peacock, J.A., Jenkins, A., White, S.D.M., Frenk, C.S., Pearce, F.R., Thomas, P.A., Efstathiou, G., Couchmann, H.M.P.: Stable clustering, the halo model and nonlinear cosmological power spectra. Mon. Not. Roy. Astron. Soc. 341, 1311 (2003) arXiv:astro-ph/0207664 [astro-ph]. https://doi.org/10.1046/j.1365-8711.2003.06503.x
(92) Aghanim, N., et al.: Planck 2018 results. VI. Cosmological parameters. Astron. Astrophys. 641, 6 (2020) arXiv:1807.06209 [astro-ph.CO]. https://doi.org/10.1051/0004-6361/201833910. [Erratum: Astron.Astrophys. 652, C4 (2021)]