Search

We study primordial non-Gaussian signatures in the redshift-space halo field on non-linear scales, using a quasi-maximum likelihood estimator based on optimally compressed power spectrum and modal bispectrum statistics. We train and validate the estimator on a suite of halo catalogues constructed from the Quijote-PNG N-body simulations, which we release to accompany this paper. We verify its unbiasedness and near optimality, for the three main types of primordial non-Gaussianity (PNG): local, equilateral, and orthogonal. We compare the modal bispectrum expansion with a $k$-binning approach, showing that the former allows for faster convergence of numerical derivatives in the computation of the score-function, thus leading to better final constraints. We find, in agreement with previous studies, that the local PNG signal in the halo-field is dominated by the scale-dependent bias signature on large scales and saturates at $k \sim 0.2~h\,\mathrm{Mpc}^{-1}$, whereas the small-scale bispectrum is the main source of information for equilateral and orthogonal PNG. Combining power spectrum and bispectrum on non-linear scales plays an important role in breaking degeneracies between cosmological and PNG parameters; such degeneracies remain however strong for equilateral PNG. We forecast that PNG parameters can be constrained with $\Delta f_\mathrm{NL}^\mathrm{local} = 45$, $\Delta f_\mathrm{NL}^\mathrm{equil} = 570$, $\Delta f_\mathrm{NL}^\mathrm{ortho} = 110$, on a cubic volume of $1 \left({ {\rm Gpc}/{ {\rm h}}} \right)^3$, at $z = 1$, considering scales up to $k_\mathrm{max} = 0.5~h\,\mathrm{Mpc}^{-1}$., Comment: 17 pages, 10 figures. v2: minor updates to match published version

We investigate how much can be learnt about four types of primordial non-Gaussianity (PNG) from small-scale measurements of the halo field. Using the QUIJOTE-PNG simulations, we quantify the information content accessible with measurements of the halo power spectrum monopole and quadrupole, the matter power spectrum, the halo-matter cross spectrum and the halo bispectrum monopole. This analysis is the first to include small, non-linear scales, up to $k_\mathrm{max}=0.5 \mathrm{h/Mpc}$, and to explore whether these scales can break degeneracies with cosmological and nuisance parameters making use of thousands of N-body simulations. We perform all the halo measurements in redshift space with a single sample comprised of all halos with mass $>3.2 \times 10^{13}~h^{-1}M_\odot$. For local PNG, measurements of the scale dependent bias effect from the power spectrum using sample variance cancellation provide significantly tighter constraints than measurements of the halo bispectrum. In this case measurements of the small scales add minimal additional constraining power. In contrast, the information on equilateral and orthogonal PNG is primarily accessible through the bispectrum. For these shapes, small scale measurements increase the constraining power of the halo bispectrum by up to $\times4$, though the addition of scales beyond $k\approx 0.3 \mathrm{h/Mpc}$ improves constraints largely through reducing degeneracies between PNG and the other parameters. These degeneracies are even more powerfully mitigated through combining power spectrum and bispectrum measurements. However even with combined measurements and small scale information, equilateral non-Gaussianity remains highly degenerate with $\sigma_8$ and our bias model., Comment: Updated to accepted version

Future large-scale structure surveys are expected to improve current bounds on primordial non-Gaussianity (PNG), with a significant impact on our understanding of early universe physics. The level of such improvements will however strongly depend on the extent to which late-time nonlinearities erase the PNG signal on small scales. In this work, we show how much primordial information remains in the bispectrum of the nonlinear dark matter density field by implementing a new, simulation-based methodology for joint estimation of PNG amplitudes (fNL) and standard Lambda cold dark matter parameters. The estimator is based on optimally compressed statistics, which, for a given input density field, combine power spectrum and modal bispectrum measurements, and numerically evaluate their covariance and their response to changes in cosmological parameters. In this first analysis, we focus on the matter density field, and we train and validate the estimator using a large suite of N-body simulations (QUIJOTE-PNG), including different types of PNG (local, equilateral, orthogonal). We explicitly test the estimator’s unbiasedness, optimality, and stability with respect to changes in the total number of input realizations. While the dark matter power spectrum itself contains negligible PNG information, as expected, including it as an ancillary statistic increases the PNG information content extracted from the bispectrum by a factor of order 2. As a result, we prove the capability of our approach to optimally extract PNG information on nonlinear scales beyond the perturbative regime, up to k max = h 0.5 Mpc 1- . At the same time, we discuss the significant information on cosmological parameters contained on these scales., G.J., M.L., and M.B. were supported by the project "Combining Cosmic Microwave Background and Large Scale Structure data: an Integrated Approach for Addressing Fundamental Questions in Cosmology," funded by the MIUR Progetti di Ricerca di Rilevante Interesse Nazionale (PRIN) Bando 2017—grant 2017YJYZAH. D.K. is supported by the South African Radio Astronomy Observatory (SARAO) and the National Research Foundation (grant No. 75415). G.J. and M.L. also acknowledge support from the INDARK INFN Initiative (https://web.infn.it/CSN4/IS/Linea5/InDark), which provided access to CINECA supercomputing facilities (https://www.cineca.it). M.B. acknowledges the use of computational resources from the parallel computing cluster of the Open Physics Hub (https://site.unibo.it/openphysicshub/en) at the Physics and Astronomy Department in Bologna. L.V. acknowledges ERC (BePreSySe, grant agreement 725327), PGC2018-098866-B-I00 MCIN/AEI/10.13039/501100011033 y FEDER "Una manera de hacer Europa," and the "Center of Excellence Maria de Maeztu 2020-2023" award to the ICCUB (CEX2019-000918-M funded by MCIN/AEI/10.13039/501100011033). B.D.W. acknowledges the support of the ANR BIG4 project, grant ANR-16-CE23-0002 of the French Agence Nationale de la Recherche; and the Labex ILP (reference ANR-10-LABX-63) part of the Idex SUPER, and received financial state aid managed by the Agence Nationale de la Recherche, as part of the program Investissements d'avenir under the reference ANR-11-IDEX-0004-02. The Flatiron Institute is supported by the Simons Foundation.

The 21cm emission of neutral hydrogen is a potential probe of the matter distribution in the Universe after reionisation. Cosmological surveys of this line intensity will be conducted in the coming years by the SKAO and HIRAX experiments, complementary to upcoming galaxy surveys. We present the first forecasts of the cosmological constraints from the combination of the 21cm power spectrum and bispectrum. Fisher forecasts are computed for the constraining power of these surveys on cosmological parameters, the BAO distance functions and the growth function. We also estimate the constraining power on dynamical dark energy and modified gravity. Finally we investigate the constraints on the 21cm clustering bias, up to second order. We consider the effects on the 21cm correlators of the telescope beam, instrumental noise, foreground avoidance, the Alcock-Paczynski effect and theoretical errors in the modelling of the correlators. Adding Planck priors, and marginalising over nuisance parameters, HIRAX achieves sub-percent precision on the $\Lambda$CDM parameters, with SKAO delivering slightly lower precision. The modified gravity parameter $\gamma$ is constrained at 1% (HIRAX) and 5% (SKAO). For the dark energy parameters $w_0,w_a$, HIRAX delivers percent-level precision while SKAO constraints are weaker. HIRAX achieves sub-percent precision on the BAO distance functions $D_A,H$, while SKAO reaches 1-2% for $0.6\lesssim z\lesssim 1$. The growth rate $f$ is constrained at a few-percent level for the whole redshift range of HIRAX and for $0.6\lesssim z\lesssim 1$ by SKAO. The different performances arise mainly since HIRAX is a packed inteferometer that is optimised for BAO measurements, while SKAO is not optimised for interferometer cosmology and operates better in single-dish mode, where the telescope beam limits access to the smaller scales that are covered by an interferometer., Comment: 37 pages, 6 figures, version accepted by JCAP

Primordial non-Gaussianity (PNG) is one of the most powerful probes of the early Universe and measurements of the large scale structure of the Universe have the potential to transform our understanding of this area. However relating measurements of the late time Universe to the primordial perturbations is challenging due to the non-linear processes that govern the evolution of the Universe. To help address this issue we release a large suite of N-body simulations containing four types of PNG: \textsc{quijote-png}. These simulations were designed to augment the \textsc{quijote} suite of simulations that explored the impact of various cosmological parameters on large scale structure observables. Using these simulations we investigate how much information on PNG can be extracted by extending power spectrum and bispectrum measurements beyond the perturbative regime at $z=0.0$. This is the first joint analysis of the PNG and cosmological information content accessible with power spectrum and bispectrum measurements of the non-linear scales. We find that the constraining power improves significantly up to $k_\mathrm{max}\approx 0.3 h/{\rm Mpc}$, with diminishing returns beyond as the statistical probes signal-to-noise ratios saturate. This saturation emphasizes the importance of accurately modelling all the contributions to the covariance matrix. Further we find that combining the two probes is a powerful method of breaking the degeneracies with the $\Lambda$CDM parameters., Comment: The simulation products are available at https://quijote-simulations.readthedocs.io/en/latest/png.html

Clustering measurements of Gravitational Wave (GW) mergers in Luminosity Distance Space can be used in the future as a powerful tool for Cosmology. We consider tomographic measurements of the Angular Power Spectrum of mergers both in an Einstein Telescope-like detector network and in some more advanced scenarios (more sources, better distance measurements, better sky localization). We produce Fisher forecasts for both cosmological (matter and dark energy) and merger bias parameters. Our fiducial model for the number distribution and bias of GW events is based on results from hydrodynamical simulations. The cosmological parameter forecasts with Einstein Telescope are less powerful than those achievable in the near future via galaxy clustering observations with, e.g., Euclid. However, in the more advanced scenarios we see significant improvements. Moreover, we show that bias can be detected at high statistical significance. Regardless of the specific constraining power of different experiments, many aspects make this type of analysis interesting anyway. For example, compact binary mergers detected by Einstein Telescope will extend up to very high redshifts. Furthermore, Luminosity Distance Space Distortions in the GW analysis have a different structure with respect to Redshift-Space Distortions in galaxy catalogues. Finally, measurements of the bias of GW mergers can provide useful insight into their physical nature and properties., Comment: 33 pages, 13 figures. Revised for publication

After reionisation, the 21cm emission line of neutral hydrogen within galaxies provides a tracer of dark matter. Next-generation intensity mapping surveys, with the SKA and other radio telescopes, will cover large sky areas and a wide range of redshifts, facilitating their use as probes of primordial non-Gaussianity. {Previous works have shown that the bispectrum can achieve tight constraints on primordial non-Gaussianity with future surveys that are purposely designed for intensity mapping in interferometer mode}. Here we investigate the constraints attainable from surveys operating in single-dish mode, \rev{using the combined power spectrum and bispectrum signal}. In the case of the power spectrum, single-dish surveys typically outperform interferometer surveys. We find that the reverse holds for the bispectrum: single-dish surveys are not competitive with surveys designed for interferometer mode., Comment: 31 pages, 5 figures, Power spectrum constraints added; version accepted by Phys. Dark Univ

The paper outlines research and development needed over the next decade to enable three-dimension mapping of the Universe using the 21 cm line of neutral hydrogen.

PUMA is a proposal for an ultra-wideband, low-resolution and transit interferometric radio telescope operating at $200-1100\,\mathrm{MHz}$. Its design is driven by six science goals which span three science themes: the physics of dark energy (measuring the expansion history and growth of the universe up to $z=6$), the physics of inflation (constraining primordial non-Gaussianity and primordial features) and the transient radio sky (detecting one million fast radio bursts and following up SKA-discovered pulsars). We propose two array configurations composed of hexagonally close-packed 6m dish arrangements with 50% fill factor. The initial 5,000 element 'petite array' is scientifically compelling, and can act as a demonstrator and a stepping stone to the full 32,000 element 'full array'. Viewed as a 21cm intensity mapping telescope, the program has the noise equivalent of a traditional spectroscopic galaxy survey comprised of 0.6 and 2.5 billion galaxies at a comoving wavenumber of $k=0.5\,h\mathrm{Mpc}^{-1}$ spanning the redshift range $z = 0.3 - 6$ for the petite and full configurations, respectively. At redshifts beyond $z=2$, the 21cm technique is a uniquely powerful way of mapping the universe, while the low-redshift range will allow for numerous cross-correlations with existing and upcoming surveys. This program is enabled by the development of ultra-wideband radio feeds, cost-effective dish construction methods, commodity radio-frequency electronics driven by the telecommunication industry and the emergence of sufficient computing power to facilitate real-time signal processing that exploits the full potential of massive radio arrays. The project has an estimated construction cost of 55 and 330 million FY19 USD for the petite and full array configurations. Including R&D, design, operations and science analysis, the cost rises to 125 and 600 million FY19 USD, respectively.

Author(s): Schlegel, David J; Kollmeier, Juna A; Aldering, Greg; Bailey, Stephen; Baltay, Charles; Bebek, Christopher; BenZvi, Segev; Besuner, Robert; Blanc, Guillermo; Bolton, Adam S; Bouri, Mohamed; Brooks, David; Buckley-Geer, Elizabeth; Cai, Zheng; Crane, Jeffrey; Dey, Arjun; Doel, Peter; Fan, Xiaohui; Ferraro, Simone; Font-Ribera, Andreu; Gutierrez, Gaston; Guy, Julien; Heetderks, Henry; Huterer, Dragan; Infante, Leopoldo; Jelinsky, Patrick; Johns, Matthew; Karagiannis, Dionysios; Kent, Stephen M; Kim, Alex G; Kneib, Jean-Paul; Kronig, Luzius; Konidaris, Nick; Lahav, Ofer; Lampton, Michael L; Lang, Dustin; Leauthaud, Alexie; Liguori, Michele; Linder, Eric V; Magneville, Christophe; Martini, Paul; Mateo, Mario; McDonald, Patrick; Miller, Christopher J; Moustakas, John; Myers, Adam D; Mulchaey, John; Newman, Jeffrey A; Nugent, Peter E; Palanque-Delabrouille, Nathalie; Padmanabhan, Nikhil; Piro, Anthony L; Poppett, Claire; Prochaska, Jason X; Pullen, Anthony R; Rabinowitz, David; Ramirez, Solange; Rix, Hans-Walter; Ross, Ashley J; Samushia, Lado; Schaan, Emmanuel; Schubnell, Michael; Seljak, Uros; Seo, Hee-Jong; Shectman, Stephen A; Silber, Joseph; Simon, Joshua D; Slepian, Zachary; Soares-Santos, Marcelle; Tarle, Greg; Thompson, Ian; Valluri, Monica; Wechsler, Risa H; White, Martin; Wilson, Michael J; Yeche, Christophe; Zaritsky, Dennis | Abstract: MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at 2

We forecast ability of dedicated 21 cm intensity mapping experiments to constraint primordial non-Gaussianity using power spectrum and bispectrum. We model the signal including the non-linear biasing expansion using a generalized halo model approach. We consider the importance of foreground filtering scale and of the foreground wedge. We find that the current generation intensity mapping experiments like CHIME do not posses sufficient sensitivity to be competitive with the existing limits. On the other hand, upcoming experiments like HIRAX can improve the current constraints and the proposed PUMA experiment can substantially improve them, reaching sensitivities below $\sigma (f_{\rm NL}), Comment: 24 pages, 4 figures, minor corrections and references updated. JCAP version

We analyse the clustering of 22,361 quasars between redshift 2.2 < z < 2.9 observed with the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS), which are included in the ninth data release (DR9). We fit the clustering results with a {\Lambda}CDM model to calculate the linear bias of the quasar sample, b=3.74\pm 0.12. The measured value of bias is consistent with the findings of White et al. (2012), where they analyse almost the same quasar sample, although only in the range, s, Comment: 19 pages, 15 figures, to be published in MNRAS

The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a radio interferometer array currently in development, with an initial 256-element array to be deployed at the South African Radio Astronomy Observatory (SARAO) Square Kilometer Array (SKA) site in South Africa. Each of the 6m, $f/0.23$ dishes will be instrumented with dual-polarisation feeds operating over a frequency range of 400-800 MHz. Through intensity mapping of the 21 cm emission line of neutral hydrogen, HIRAX will provide a cosmological survey of the distribution of large-scale structure over the redshift range of $0.775 < z < 2.55$ over $\sim$15,000 square degrees of the southern sky. The statistical power of such a survey is sufficient to produce $\sim$7 percent constraints on the dark energy equation of state parameter when combined with measurements from the Planck satellite. Additionally, HIRAX will provide a highly competitive platform for radio transient and HI absorber science while enabling a multitude of cross-correlation studies. In this paper, we describe the science goals of the experiment, overview of the design and status of the sub-components of the telescope system, and describe the expected performance of the initial 256-element array as well as the planned future expansion to the final, 1024-element array., 24 pages, 6 figures. Updated to match published version