Search

We present the results of a set of cosmological N-body simulations with standard $\Lambda$CDM cosmology but characterized by a scale-dependent primordial non-Gaussianity of the local type featuring a power-law dependence of the $f_{\rm NL}^{\rm loc}(k)$ at large scales followed by a saturation to a constant value at smaller scales where non-linear growth leads to the formation of collapsed cosmic structures. Such models are built to ensure consistency with current Cosmic Microwave Background bounds on primordial non-Gaussianity yet allowing for large effects of the non-Gaussian statistics on the properties of non-linear structure formation. We show the impact of such scale-dependent non-Gaussian scenarios on a wide range of properties of the resulting cosmic structures, such as the non-linear matter power spectrum, the halo and sub-halo mass functions, the concentration-mass relation, the halo and void density profiles, and we highlight for the first time that some of these models might mimic the effects of Warm Dark Matter for several of such observables, Comment: 21 pages, 9 figures, 2 tables; to be submitted to JCAP

A major goal of cosmology is to understand the nature of the field(s) which drove primordial Inflation. Through future observations, the statistics of large-scale structure will allow us to probe primordial non-Gaussianity of the curvature perturbation at the end of Inflation. We show how a new correlation statistic can significantly improve these constraints over conventional methods. Next-generation radio telescope arrays are under construction which will map the density field of neutral hydrogen to high redshifts. These telescopes can operate as an interferometer, able to probe small scales, or as a collection of single dishes, combining signals to map the large scales. We show how to fuse these these operating modes in order to measure the squeezed bispectrum with higher precision and greater economy. This leads to constraints on primordial non-Gaussianity that will improve on measurements by Planck, and out-perform other surveys such as Euclid. We forecast that $\sigma(f_{\mathrm{NL}}^{\mathrm{loc}})\sim 3$, achieved by using a small subset, $\mathcal{O}(10^2-10^3)$, of the total number of accessible triangles. The proposed method can be used for any signal that peaks in squeezed configurations., Comment: 13 pages, 6 figures

We explore the information content of 21cm intensity maps in redshift space using the 1-loop Effective Field Theory power spectrum model and the bispectrum at tree level. The 21cm signal contains signatures of dark matter, dark energy and the growth of large-scale structure in the Universe. These signatures are typically analyzed via the 2-point correlation function or power spectrum. However, adding the information from the 3-point correlation function or bispectrum will be crucial to exploiting next-generation intensity mapping experiments. The bispectrum could offer a unique opportunity to break key parameter degeneracies that hinder the measurement of cosmological parameters and improve on the precision. We use a Fisher forecast analysis to estimate the constraining power of the HIRAX survey on cosmological parameters, dark energy and modified gravity., Comment: 16 + 4 pages, "Minor corrections, main results unchanged. Version accepted by Phys. Dark Univ"

Primordial non-Gaussianity of the local type induces a strong scale-dependent bias on the clustering of halos in the late-time Universe. This signature is particularly promising to provide constraints on the non-Gaussianity parameter $f_{\rm NL}$ from galaxy surveys, as the bias amplitude grows with scale and becomes important on large, linear scales. However, there is a well-known degeneracy between the real prize, the $f_{\rm NL}$ parameter, and the (non-Gaussian) assembly bias i.e., the halo formation history-dependent contribution to the amplitude of the signal, which could seriously compromise the ability of large-scale structure surveys to constrain $f_{\rm NL}$. We show how the assembly bias can be modeled and constrained, thus almost completely recovering the power of galaxy surveys to competitively constrain primordial non-Gaussianity. In particular, studying hydrodynamical simulations, we find that a proxy for the halo properties that determine assembly bias can be constructed from photometric properties of galaxies. Using a prior on the assembly bias guided by this proxy degrades the statistical errors on $f_{\rm NL}$ only mildly compared to an ideal case where the assembly bias is perfectly known. The systematic error on $f_{\rm NL}$ that the proxy induces can be safely kept under control., Comment: 30 pages, 13 figures. v2: minor updates to match accepted version

We study signatures of primordial non-Gaussianity (PNG) in the redshift-space halo field on non-linear scales, using a combination of three summary statistics, namely the halo mass function (HMF), power spectrum, and bispectrum. The choice of adding the HMF to our previous joint analysis of power spectrum and bispectrum is driven by a preliminary field-level analysis, in which we train graph neural networks on halo catalogues to infer the PNG $f_\mathrm{NL}$ parameter. The covariance matrix and the responses of our summaries to changes in model parameters are extracted from a suite of halo catalogues constructed from the Quijote-PNG N-body simulations. We consider the three main types of PNG: local, equilateral and orthogonal. Adding the HMF to our previous joint analysis of power spectrum and bispectrum produces two main effects. First, it reduces the equilateral $f_\mathrm{NL}$ predicted errors by roughly a factor $2$, while also producing notable, although smaller, improvements for orthogonal PNG. Second, it helps break the degeneracy between the local PNG amplitude, $f_\mathrm{NL}^\mathrm{local}$, and assembly bias, $b_{\phi}$, without relying on any external prior assumption. Our final forecasts for PNG parameters are $\Delta f_\mathrm{NL}^\mathrm{local} = 40$, $\Delta f_\mathrm{NL}^\mathrm{equil} = 210$, $\Delta f_\mathrm{NL}^\mathrm{ortho} = 91$, on a cubic volume of $1 \left(h^{-1}{\rm Gpc}\right)^3$, with a halo number density of $\bar{n}\sim 5.1 \times 10^{-5}~h^3\mathrm{Mpc}^{-3}$, at $z = 1$, and considering scales up to $k_\mathrm{max} = 0.5~h\,\mathrm{Mpc}^{-1}$., Comment: 17 pages, 11 figures. v3 (minor caption fix)

The power spectrum and bispectrum of dark matter tracers are key and complementary probes of the Universe. Next-generation surveys will deliver good measurements of the bispectrum, opening the door to improved cosmological constraints and the breaking of parameter degeneracies, from the combination of the power spectrum and bispectrum. Multi-tracer power spectra have been used to suppress cosmic variance and mitigate the effects of nuisance parameters and systematics. We present a bispectrum multi-tracer formalism that can be applied to next-generation survey data. Then we perform a simple Fisher analysis to illustrate qualitatively the improved precision on primordial non-Gaussianity that is expected to come from the bispectrum multi-tracer. In addition, we investigate the parametric dependence of conditional errors from multi-tracer power spectra and multi-tracer bispectra, on the differences between the biases and the number densities of two tracers. Our results suggest that optimal constraints arise from maximising the ratio of number densities, the difference between the linear biases, the difference between the quadratic biases, and the difference between the products $b_1\,b_\Phi$ for each tracer, where $b_\Phi$ is the bias for the primordial potential., Comment: 30 pages, 7 figures, JCAP version

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

In this white paper, we present the MegaMapper concept. The MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at $2

In this white paper, we present an experimental road map for spectroscopic experiments beyond DESI. DESI will be a transformative cosmological survey in the 2020s, mapping 40 million galaxies and quasars and capturing a significant fraction of the available linear modes up to z=1.2. DESI-II will pilot observations of galaxies both at much higher densities and extending to higher redshifts. A Stage-5 experiment would build out those high-density and high-redshift observations, mapping hundreds of millions of stars and galaxies in three dimensions, to address the problems of inflation, dark energy, light relativistic species, and dark matter. These spectroscopic data will also complement the next generation of weak lensing, line intensity mapping and CMB experiments and allow them to reach their full potential., Comment: Contribution to Snowmass 2021

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

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

Future Large Scale Structure surveys are expected to improve over 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 non-linearities erase the PNG signal on small scales. In this work, we show how much primordial information remains in the bispectrum of the non-linear dark matter density field by implementing a new, simulation-based, methodology for joint estimation of PNG amplitudes ($f_{\rm NL}$) and standard $\Lambda$CDM 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. 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 non-linear scales beyond the perturbative regime, up to $k_{\rm max} = 0.5~h\,{\rm Mpc}^{-1}$, obtaining marginalized $1$-$\sigma$ bounds of $\Delta f_{\rm NL}^{\rm local} \sim 16$, $\Delta f_{\rm NL}^{\rm equil} \sim 77$ and $\Delta f_{\rm NL}^{\rm ortho} \sim 40$ on a cubic volume of $1~(\mathrm{Gpc}/h)^3$ at $z=1$. At the same time, we discuss the significant information on cosmological parameters contained on these scales., Comment: 22 pages, 12 figures

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

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., Comment: 24 pages, 6 figures. Updated to match published version

We study the clustering of Gravitational Wave (GW) merger events and Supernovae IA (SN), as cosmic tracers in Luminosity Distance Space. We modify the publicly available CAMB code to numerically evaluate auto- and cross- power spectra for the different sources, including Luminosity Distance Space distortion effects generated by peculiar velocities and lensing convergence. We perform a multitracer Fisher analysis to forecast expected constraints on cosmological and GW bias coefficients, using outputs from hydrodinamical N-body simulations to determine the bias fiducial model and considering future observations from the Vera Rubin Observatory and Einstein Telescope (ET), both single and in a 3 detector network configuration. We find that adding SN to the GW merger dataset considerably improves the forecast, mostly by breaking significant parameter degeneracies, with final constraints comparable to those obtainable from a Euclid-like survey. GW merger bias is forecasted to be detectable with good significance even in the single ET case., Comment: Version accepted for publication on JCAP, 25 pages, 10 figures

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 Packed Ultra-wideband Mapping Array (PUMA) is a proposed low-resolution transit interferometric radio telescope operating over the frequency range 200 - 1100MHz. Its rich science portfolio will include measuring structure in the universe from redshift z = 0.3 to 6 using 21cm intensity mapping, detecting one million fast radio bursts, and monitoring thousands of pulsars. It will allow PUMA to advance science in three different areas of physics (the physics of dark energy, the physics of cosmic inflation and time-domain astrophysics). This document is a response to a request for information (RFI) by the Panel on Radio, Millimeter, and Submillimeter Observations from the Ground (RMS) of the Decadal Survey on Astronomy and Astrophysics 2020. We present the science case of PUMA, the development path and major risks to the project., Comment: 46 pages, 16 figures, 7 tables; response to the request for information (RFI) by the Panel on Radio, Millimeter, and Submillimeter Observations from the Ground (RMS) of the Astro2020 Decadal Survey regarding PUMA APC submission (arXiv:1907.12559); v2: updated with correct bbl file

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})<5$ for equilateral and orthogonal configurations and $\sigma( f_{\rm NL}) < 1$ for the local shape if good foreground control is achieved., Comment: 24 pages, 4 figures, minor corrections and references updated. JCAP version

Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and analyze such data. On the positive side, it has become clear that the outstanding issues are purely technical in nature and can be solved with sufficient development activity. Such activity will enable science across redshifts, from early galaxy evolution in the pre-reionization era to dark energy evolution at low redshift., Comment: 10 pages + references, 2 figures, 1 table; APC white paper submitted to the Astro2020 decadal survey. arXiv admin note: text overlap with arXiv:1810.09572

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., Comment: 10 pages + references, 3 figures, 3 tables; project white paper submitted to the Astro2020 decadal survey; further details in updated arXiv:1810.09572

MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at 2

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space., Comment: Science white paper submitted to the Astro2020 Decadal Survey

This white paper envisions a revolutionary post-DESI, post-LSST dark energy program based on intensity mapping of the redshifted 21cm emission line from neutral hydrogen at radio frequencies. The proposed intensity mapping survey has the unique capability to quadruple the volume of the Universe surveyed by optical programs, provide a percent-level measurement of the expansion history to $z \sim 6$, open a window to explore physics beyond the concordance $\Lambda$CDM model, and to significantly improve the precision on standard cosmological parameters. In addition, characterization of dark energy and new physics will be powerfully enhanced by cross-correlations with optical surveys and cosmic microwave background measurements. The rich dataset obtained by the proposed intensity mapping instrument will be simultaneously useful in exploring the time-domain physics of fast radio transients and pulsars, potentially in live "multi-messenger" coincidence with other observatories. The core dark energy/inflation science advances enabled by this program are the following: (i) Measure the expansion history of the universe over $z=0.3-6$ with a single instrument, extending the range deep into the pre-acceleration era, providing an unexplored window for new physics; (ii) Measure the growth rate of structure in the universe over the same redshift range; (iii) Observe, or constrain, the presence of inflationary relics in the primordial power spectrum, improving existing constraints by an order of magnitude; (iv) Observe, or constrain, primordial non-Gaussianity with unprecedented precision, improving constraints on several key numbers by an order of magnitude. Detailed mapping of the enormous, and still largely unexplored, volume of cosmic space will thus provide unprecedented information on fundamental questions of the vacuum energy and early-universe physics., Comment: This white-paper was prepared under auspices of the Department of Energy Cosmic Visions Dark Energy program; v1: initial release in Oct 2018; v2: Jul 2019 update: strawman design matches PUMA Decadal Submission, updated forecasting, new figures and science section; v3: fixed reference, missing author in metadata

We forecast constraints on primordial non-Gaussianity (PNG) and bias parameters from measurements of galaxy power spectrum and bispectrum in future radio continuum and optical surveys. In the galaxy bispectrum, we consider a comprehensive list of effects, including the bias expansion for non-Gaussian initial conditions up to second order, redshift space distortions, redshift uncertainties and theoretical errors. These effects are all combined in a single PNG forecast for the first time. Moreover, we improve the bispectrum modelling over previous forecasts, by accounting for trispectrum contributions. All effects have an impact on final predicted bounds, which varies with the type of survey. We find that the bispectrum can lead to improvements up to a factor $\sim 5$ over bounds based on the power spectrum alone, leading to significantly better constraints for local-type PNG, with respect to current limits from \textit{Planck}. Future radio and photometric surveys could obtain a measurement error of $\sigma(f_{\mathrm{NL}}^{\mathrm{loc}}) \approx 0.2$. In the case of equilateral PNG, galaxy bispectrum can improve upon present bounds only if significant improvements in the redshift determinations of future, large volume, photometric or radio surveys could be achieved. For orthogonal non-Gaussianity, expected constraints are generally comparable to current ones., Comment: 36 pages, 19 figures, 12 tables, MNRAS 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<40 Mpc/h. At large scales we observe an excess or plateau in the clustering correlation function. By fitting a model that incorporates a scale dependent additional term in the bias introduced by primordial non-Gaussianity of the local type, we calculate the amplitude of the deviation from the Gaussian initial conditions as 70