Search

The large scale limit of the galaxy power spectrum provides a unique window into the early Universe through a possible detection of scale dependent bias produced by primordial non Gaussianities. On such large scales, relativistic effects could become important and be confused for a primordial signal. In this Letter we provide the first consistent estimate of such effects in the observed galaxy power spectrum, and discuss their possible degeneracy with local primordial non Gaussianities. We also clarify the physical differences between the two signatures, as revealed by their different sensitivity to the large scale gravitational potential. Our results indicate that, while relativistic effects could easily account for 10% of the observed power spectrum, the subset of those with a similar scale dependence to a primordial signal can be safely ignored for current galaxy surveys, but it will become relevant for future observational programs., Comment: 4 pages, 2 figures, code available on GitHub at https://github.com/foglienimatteo/GaPSE.jl. v2: minor modifications to match published version

Primordial non-Gaussianities (PNG) leave unique signatures in the bispectrum of the large-scale structure. With upcoming galaxy surveys set to improve PNG constraints by at least one order of magnitude, it is important to account for any potential contamination. In our work we show how to include wide-angle effects into the 3-dimensional observed galaxy bispectrum. We compute the leading wide-angle corrections to the monopole, finding that they could mimic local PNG with an amplitude of $f_{\rm NL} = \mathcal{O}\left( 0.1 \right)$. We also compute the dipole induced by wide-angle effects, whose amplitude is a few-percent of the flat-sky monopole. We estimate that wide-angle effects in the monopole can be safely neglected for survey volumes of the order of $8~\mathrm{Gpc}^3 h^{-3}$, while the dipole can start being detected from surveys probing volumes larger than $50~\mathrm{Gpc}^3 h^{-3}$. Our formalism can be readily adapted to realistic survey geometries and to include relativistic effects, which may become relevant at high redshifts., Comment: 11+15 pages, 6 figures. v2: minor modifications to match published version

We present constraints on the amplitude of local Primordial Non-Gaussianities (PNG), $f_{\rm NL}$, using the quasar sample in the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey Data Release 16. We analyze the power spectrum monopole, testing for the presence of scale dependent galaxy bias induced by local PNG. Our analysis makes use of optimal redshift weights that maximize the response of the quasar sample to the possible presence of non zero PNG. We find $-4

Cosmology offers opportunities to test Dark Matter independently of its interactions with the Standard Model. We study the imprints of long-range forces acting solely in the dark sector on the distribution of galaxies, the so-called Large Scale Structure (LSS). We derive the strongest constraint on such forces from a combination of Planck and BOSS data. Along the way we consistently develop, for the first time, the Effective Field Theory of LSS in the presence of new dynamics in the dark sector. We forecast that future surveys will improve the current bound by an order of magnitude., Comment: 5+18 pages, 4 figures

The large scale limit of the galaxy power spectrum provides a unique window into the early Universe through a possible detection of scale dependent bias produced by primordial non Gaussianities. On such large scales, relativistic effects could become important and be confused for a primordial signal. In this Letter we provide the first consistent estimate of such effects in the observed galaxy power spectrum, and discuss their possible degeneracy with local primordial non Gaussianities. We also clarify the physical differences between the two signatures, as revealed by their different sensitivity to the large scale gravitational potential. Our results indicate that, while relativistic effects could easily account for 10% of the observed power spectrum, the subset of those with a similar scale dependence to a primordial signal can be safely ignored for current galaxy surveys, but it will become relevant for future observational programs., Comment: 4 pages, 2 figures, code available on GitHub at https://github.com/foglienimatteo/GaPSE.jl. v2: minor modifications to match published version

Primordial non-Gaussianities (PNG) leave unique signatures in the bispectrum of the large-scale structure. With upcoming galaxy surveys set to improve PNG constraints by at least one order of magnitude, it is important to account for any potential contamination. In our work we show how to include wide-angle effects into the 3-dimensional observed galaxy bispectrum. We compute the leading wide-angle corrections to the monopole, finding that they could mimic local PNG with an amplitude of $f_{\rm NL} = \mathcal{O}\left( 0.1 \right)$. We also compute the dipole induced by wide-angle effects, whose amplitude is a few-percent of the flat-sky monopole. Our formalism can be readily adapted to realistic survey geometries and to include relativistic effects, which may become relevant at high redshifts., Comment: 11+15 pages, 6 figures

Cosmic inflation provides a window to the highest energy densities accessible in nature, far beyond those achievable in any realistic terrestrial experiment. Theoretical insights into the inflationary era and its observational probes may therefore shed unique light on the physical laws underlying our universe. This white paper describes our current theoretical understanding of the inflationary era, with a focus on the statistical properties of primordial fluctuations. In particular, we survey observational targets for three important signatures of inflation: primordial gravitational waves, primordial non-Gaussianity and primordial features. With the requisite advancements in analysis techniques, the tremendous increase in the raw sensitivities of upcoming and planned surveys will translate to leaps in our understanding of the inflationary paradigm and could open new frontiers for cosmology and particle physics. The combination of future theoretical and observational developments therefore offer the potential for a dramatic discovery about the nature of cosmic acceleration in the very early universe and physics on the smallest scales., Comment: Contribution to Snowmass 2021; 103 pages, 8 figures, 378 endorsers; v2: minor changes

We initiate the exploration of the cosmology of dark fifth forces: new forces acting solely on Dark Matter. We focus on long range interactions which lead to an effective violation of the Equivalence Principle on cosmological scales. At the microscopic level, the dark fifth force can be realized by a light scalar with mass smaller than the Hubble constant today ($\lesssim 10^{-33}\,\text{eV}$) coupled to Dark Matter. We study the behavior of the background cosmology and linear perturbations in such a Universe. At the background level, the new force modifies the evolution of the Dark Matter energy density and thus the Hubble flow. At linear order, it modifies the growth of matter perturbations and generates relative density and velocity perturbations between Dark Matter and baryons that grow over time. We derive constraints from current CMB and BAO data, bounding the strength of the dark fifth force to be less than a percent of gravity. These are the strongest constraints to date. We present potential implications of this scenario for the Hubble tension and discuss how our results are modified if the light scalar mediator accounts for the observed density of the Dark Energy. Finally, we comment on the interplay between our constraints and searches for violations of the Equivalence Principle in the visible sector., Comment: 50 pages, 15 figures

Cosmic inflation provides a window to the highest energy densities accessible in nature, far beyond those achievable in any realistic terrestrial experiment. Theoretical insights into the inflationary era and its observational probes may therefore shed unique light on the physical laws underlying our universe. This white paper describes our current theoretical understanding of the inflationary era, with a focus on the statistical properties of primordial fluctuations. In particular, we survey observational targets for three important signatures of inflation: primordial gravitational waves, primordial non-Gaussianity and primordial features. With the requisite advancements in analysis techniques, the tremendous increase in the raw sensitivities of upcoming and planned surveys will translate to leaps in our understanding of the inflationary paradigm and could open new frontiers for cosmology and particle physics. The combination of future theoretical and observational developments therefore offer the potential for a dramatic discovery about the nature of cosmic acceleration in the very early universe and physics on the smallest scales., Comment: Contribution to Snowmass 2021; 99 pages, 7 figures, 362 endorsers

Many well-motivated extensions of the Standard Model predict the existence of new light species that may have been produced in the early universe. Prominent examples include axions, sterile neutrinos, gravitinos, dark photons, and more. The gravitational influence of light relics leaves imprints in the cosmic microwave background fluctuations, the large-scale structure of the universe and the primordial element abundances. In this paper, we detail the physics of cosmological light relics, and describe how measurements of their relic density and mass serve as probes of physics beyond the Standard Model. A measurement of the light relic density at the precision of upcoming cosmological surveys will point the way toward new physics or severely constrain the range of viable extensions to the Standard Model., Comment: Contribution to Snowmass 2021. 59 pages, 12 figures

We present a method to perform the exact convolution of the model prediction for bispectrum multipoles in redshift space with the survey window function. We extend a widely applied method for the power spectrum convolution to the bispectrum, taking advantage of a 2D-FFTlog algorithm. As a preliminary test of its accuracy, we consider the toy model of a spherical window function in real space. This setup provides an analytical evaluation of the 3-point function of the window, and therefore it allows to isolate and quantify possible systematic errors of the method. We find that our implementation of the convolution in terms of a mixing matrix shows differences at the percent level in comparison to the measurements from a very large set of mock halo catalogs. It is also able to recover unbiased constraints on halo bias parameters in a likelihood analysis of a set of numerical simulations with a total volume of $100\, h^{-3} \, {\rm Gpc}^3$. For the level of accuracy required by these tests, the multiplication with the mixing matrix is performed in the time of one second or less., Comment: 18+10 pages, 10 figures, 2 tables

We make a number of comments regarding modeling degeneracies in strong lensing measurements of the Hubble parameter $H_0$. The first point concerns the impact of weak lensing associated with different segments of the line of sight. We show that external convergence terms associated with the lens-source and observer-lens segments need to be included in cosmographic modeling, in addition to the usual observer-source term, to avoid systematic bias in the inferred value of $H_0$. Specifically, we show how an incomplete account of some line of sight terms biases stellar kinematics as well as ray tracing simulation methods to alleviate the mass sheet degeneracy. The second point concerns the use of imaging data for multiple strongly-lensed sources in a given system. We show that the mass sheet degeneracy is not fully resolved by the availability of multiple sources: some degeneracy remains because of differential external convergence between the different sources. Similarly, differential external convergence also complicates the use of multiple sources in addressing the approximate mass sheet degeneracy associated with a local ("internal") core component in lens galaxies. This internal-external degeneracy is amplified by the non-monotonicity of the angular diameter distance as a function of redshift. For a rough assessment of the weak lensing effects, we provide estimates of external convergence using the nonlinear matter power spectrum, paying attention to non-equal time correlators., Comment: 16 pages, 6 figures. Comments welcome!

Many well-motivated extensions of the Standard Model predict the existence of new light species that may have been produced in the early universe. Prominent examples include axions, sterile neutrinos, gravitinos, dark photons, and more. The gravitational influence of light relics leaves imprints in the cosmic microwave background fluctuations, the large-scale structure of the universe and the primordial element abundances. In this paper, we detail the physics of cosmological light relics, and describe how measurements of their relic density and mass serve as probes of physics beyond the Standard Model. A measurement of the light relic density at the precision of upcoming cosmological surveys will point the way toward new physics or severely constrain the range of viable extensions to the Standard Model., Comment: Contribution to Snowmass 2021. 59 pages, 12 figures

Cosmic inflation provides a window to the highest energy densities accessible in nature, far beyond those achievable in any realistic terrestrial experiment. Theoretical insights into the inflationary era and its observational probes may therefore shed unique light on the physical laws underlying our universe. This white paper describes our current theoretical understanding of the inflationary era, with a focus on the statistical properties of primordial fluctuations. In particular, we survey observational targets for three important signatures of inflation: primordial gravitational waves, primordial non-Gaussianity and primordial features. With the requisite advancements in analysis techniques, the tremendous increase in the raw sensitivities of upcoming and planned surveys will translate to leaps in our understanding of the inflationary paradigm and could open new frontiers for cosmology and particle physics. The combination of future theoretical and observational developments therefore offer the potential for a dramatic discovery about the nature of cosmic acceleration in the very early universe and physics on the smallest scales., Comment: Contribution to Snowmass 2021; 99 pages, 7 figures, 362 endorsers

We make a number of comments regarding modeling degeneracies in strong lensing measurements of the Hubble parameter $H_0$. The first point concerns the impact of weak lensing associated with different segments of the line of sight. We show that external convergence terms associated with the lens-source and observer-lens segments need to be included in cosmographic modeling, in addition to the usual observer-source term, to avoid systematic bias in the inferred value of $H_0$. Specifically, we show how an incomplete account of some line of sight terms biases stellar kinematics as well as ray tracing simulation methods to alleviate the mass sheet degeneracy. The second point concerns the use of imaging data for multiple strongly-lensed sources in a given system. We show that the mass sheet degeneracy is not fully resolved by the availability of multiple sources: some degeneracy remains because of differential external convergence between the different sources. Similarly, differential external convergence also complicates the use of multiple sources in addressing the approximate mass sheet degeneracy associated with a local ("internal") core component in lens galaxies. This internal-external degeneracy is amplified by the non-monotonicity of the angular diameter distance as a function of redshift. For a rough assessment of the weak lensing effects, we provide estimates of external convergence using the nonlinear matter power spectrum, paying attention to non-equal time correlators., Comment: 16 pages, 6 figures. v2: match JCAP version. Comments welcome!

Cosmic inflation provides a window to the highest energy densities accessible in nature, far beyond those achievable in any realistic terrestrial experiment. Theoretical insights into the inflationary era and its observational probes may therefore shed unique light on the physical laws underlying our universe. This white paper describes our current theoretical understanding of the inflationary era, with a focus on the statistical properties of primordial fluctuations. In particular, we survey observational targets for three important signatures of inflation: primordial gravitational waves, primordial non-Gaussianity and primordial features. With the requisite advancements in analysis techniques, the tremendous increase in the raw sensitivities of upcoming and planned surveys will translate to leaps in our understanding of the inflationary paradigm and could open new frontiers for cosmology and particle physics. The combination of future theoretical and observational developments therefore offer the potential for a dramatic discovery about the nature of cosmic acceleration in the very early universe and physics on the smallest scales., Comment: Contribution to Snowmass 2021; 103 pages, 8 figures, 378 endorsers; v2: minor changes

We summarize progress made in theoretical astrophysics and cosmology over the past decade and areas of interest for the coming decade. This Report is prepared as the TF09 "Astrophysics and Cosmology" topical group summary for the Theory Frontier as part of the Snowmass 2021 process., Comment: 57 pages

We summarize progress made in theoretical astrophysics and cosmology over the past decade and areas of interest for the coming decade. This Report is prepared as the TF09 "Astrophysics and Cosmology" topical group summary for the Theory Frontier as part of the Snowmass 2021 process., Comment: 57 pages

We summarize progress made in theoretical astrophysics and cosmology over the past decade and areas of interest for the coming decade. This Report is prepared as the TF09 "Astrophysics and Cosmology" topical group summary for the Theory Frontier as part of the Snowmass 2021 process., Comment: 57 pages

We initiate the exploration of the cosmology of dark fifth forces: new forces acting solely on Dark Matter. We focus on long range interactions which lead to an effective violation of the Equivalence Principle on cosmological scales today. At the microscopic level, the dark fifth force can be realized by a light scalar with mass smaller than the Hubble constant today ($\lesssim 10^{-33}\,\text{eV}$) coupled to Dark Matter. We study the behavior of the background cosmology and linear perturbations in such a Universe. At the background level, the new force modifies the evolution of the Dark Matter energy density and thus the Hubble flow. At linear order, it modifies the growth of matter perturbations and generates relative density and velocity perturbations between Dark Matter and baryons that grow over time. We derive constraints from current CMB and BAO data, bounding the strength of the dark fifth force to be less than a percent of gravity. These are the strongest constraints to date. We present potential implications of this scenario for the Hubble tension and discuss how our results are modified if the light scalar mediator accounts for the observed density of the Dark Energy. Finally, we comment on the interplay between our constraints and searches for violations of the Equivalence Principle in the visible sector., Comment: 54 pages, 15 figures. v2: revision of Section 6, now including atomic clock constraints; added Table 1 summarizing the main results; extended discussion of BAO; added several references. Matches version to appear in JCAP. v3: fixed some incorrect equations, all results and figures unaffected. Updated with latest experimental constraints from MICROSCOPE and atomic clocks

Many well-motivated extensions of the Standard Model predict the existence of new light species that may have been produced in the early universe. Prominent examples include axions, sterile neutrinos, gravitinos, dark photons, and more. The gravitational influence of light relics leaves imprints in the cosmic microwave background fluctuations, the large-scale structure of the universe and the primordial element abundances. In this paper, we detail the physics of cosmological light relics, and describe how measurements of their relic density and mass serve as probes of physics beyond the Standard Model. A measurement of the light relic density at the precision of upcoming cosmological surveys will point the way toward new physics or severely constrain the range of viable extensions to the Standard Model., Comment: Contribution to Snowmass 2021. 59 pages, 12 figures

We present a method to perform the exact convolution of the model prediction for bispectrum multipoles in redshift space with the survey window function. We extend a widely applied method for the power spectrum convolution to the bispectrum, taking advantage of a 2D-FFTlog algorithm. As a preliminary test of its accuracy, we consider the toy model of a spherical window function in real space. This setup provides an analytical evaluation of the 3-point function of the window, and therefore it allows to isolate and quantify possible systematic errors of the method. We find that our implementation of the convolution in terms of a mixing matrix shows differences at the percent level in comparison to the measurements from a very large set of mock halo catalogs. It is also able to recover unbiased constraints on halo bias parameters in a likelihood analysis of a set of numerical simulations with a total volume of $100\, h^{-3} \, {\rm Gpc}^3$. For the level of accuracy required by these tests, the multiplication with the mixing matrix is performed in the time of one second or less., Comment: 18+10 pages, 10 figures, 2 tables

Cosmic inflation provides a window to the highest energy densities accessible in nature, far beyond those achievable in any realistic terrestrial experiment. Theoretical insights into the inflationary era and its observational probes may therefore shed unique light on the physical laws underlying our universe. This white paper describes our current theoretical understanding of the inflationary era, with a focus on the statistical properties of primordial fluctuations. In particular, we survey observational targets for three important signatures of inflation: primordial gravitational waves, primordial non-Gaussianity and primordial features. With the requisite advancements in analysis techniques, the tremendous increase in the raw sensitivities of upcoming and planned surveys will translate to leaps in our understanding of the inflationary paradigm and could open new frontiers for cosmology and particle physics. The combination of future theoretical and observational developments therefore offer the potential for a dramatic discovery about the nature of cosmic acceleration in the very early universe and physics on the smallest scales., Comment: Contribution to Snowmass 2021; 103 pages, 8 figures, 378 endorsers; v2: minor changes

We make a number of comments regarding modeling degeneracies in strong lensing measurements of the Hubble parameter $H_0$. The first point concerns the impact of weak lensing associated with different segments of the line of sight. We show that external convergence terms associated with the lens-source and observer-lens segments need to be included in cosmographic modeling, in addition to the usual observer-source term, to avoid systematic bias in the inferred value of $H_0$. Specifically, we show how an incomplete account of some line of sight terms biases stellar kinematics as well as ray tracing simulation methods to alleviate the mass sheet degeneracy. The second point concerns the use of imaging data for multiple strongly-lensed sources in a given system. We show that the mass sheet degeneracy is not fully resolved by the availability of multiple sources: some degeneracy remains because of differential external convergence between the different sources. Similarly, differential external convergence also complicates the use of multiple sources in addressing the approximate mass sheet degeneracy associated with a local ("internal") core component in lens galaxies. This internal-external degeneracy is amplified by the non-monotonicity of the angular diameter distance as a function of redshift. For a rough assessment of the weak lensing effects, we provide estimates of external convergence using the nonlinear matter power spectrum, paying attention to non-equal time correlators., Comment: 16 pages, 6 figures. v3: corrected mistake in Sec.III, changing details in weak lensing bias due to kinematics and imaging combination

The line-of-sight peculiar velocities of galaxies contribute to their observed redshifts, breaking the translational invariance of galaxy clustering down to a rotational invariance around the observer. This becomes important when the line-of-sight direction varies significantly across a survey, leading to what are known as `wide angle' effects in redshift space distortions. Wide-angle effects will also be present in measurements of the momentum field, i.e. the galaxy density-weighted velocity field, in upcoming peculiar velocity surveys. In this work we study how wide-angle effects modify the predicted correlation function and power spectrum for momentum statistics, both in auto-correlation and in cross-correlation with the density field. Using both linear theory and the Zeldovich approximation, we find that deviations from the plane-parallel limit are large and could become important in data analysis for low redshift surveys. We point out that even multipoles in the cross-correlation between density and momentum are non-zero regardless of the choice of line of sight, and therefore contain new cosmological information that could be exploited. We discuss configuration-space, Fourier-space and spherical analyses, providing exact expressions in each case rather than relying on an expansion in small angles. We hope these expressions will be of use in the analysis of upcoming surveys for redshift-space distortions and peculiar velocities.

The line-of-sight peculiar velocities of galaxies contribute to their observed redshifts, breaking the translational invariance of galaxy clustering down to a rotational invariance around the observer. This becomes important when the line-of-sight direction varies significantly across a survey, leading to what are known as `wide angle' effects in redshift space distortions. Wide-angle effects will also be present in measurements of the momentum field, i.e. the galaxy density-weighted velocity field, in upcoming peculiar velocity surveys. In this work we study how wide-angle effects modify the predicted correlation function and power spectrum for momentum statistics, both in auto-correlation and in cross-correlation with the density field. Using both linear theory and the Zeldovich approximation, we find that deviations from the plane-parallel limit are large and could become important in data analysis for low redshift surveys. We point out that even multipoles in the cross-correlation between density and momentum are non-zero regardless of the choice of line of sight, and therefore contain new cosmological information that could be exploited. We discuss configuration-space, Fourier-space and spherical analyses, providing exact expressions in each case rather than relying on an expansion in small angles. We hope these expressions will be of use in the analysis of upcoming surveys for redshift-space distortions and peculiar velocities.

An observational program focused on the high redshift ($2

The long wavelength modes lost to bright foregrounds in the interferometric 21-cm surveys can partially be recovered using a forward modeling approach that exploits the non-linear coupling between small and large scales induced by gravitational evolution. In this work, we build upon this approach by considering how adding external galaxy distribution data can help to fill in these modes. We consider supplementing the 21-cm data at two different redshifts with a spectroscopic sample (good radial resolution but low number density) loosely modeled on DESI-ELG at $z=1$ and a photometric sample (high number density but poor radial resolution) similar to LSST sample at $z=1$ and $z=4$ respectively. We find that both the galaxy samples are able to reconstruct the largest modes better than only using 21-cm data, with the spectroscopic sample performing significantly better than the photometric sample despite much lower number density. We demonstrate the synergies between surveys by showing that the primordial initial density field is reconstructed better with the combination of surveys than using either of them individually. Methodologically, we also explore the importance of smoothing the density field when using bias models to forward model these tracers for reconstruction.

Line intensity mapping (LIM) is an emerging technique with a unique potential to probe a wide range of scales and redshifts. Realizing the full potential of LIM, however, relies on accurate modeling of the signal. We introduce an extended halo model for the power spectrum of intensity fluctuations of CO rotational lines and [CII] fine transition line in real space, modeling nonlinearities in matter fluctuations and biasing relation between the line intensity fluctuations and the underlying dark matter distribution. We also compute the stochastic contributions beyond the Poisson approximation using the halo model framework. To establish the accuracy of the model, we create the first cosmological-scale simulations of CO and [CII] intensity maps, \textsf{MithraLIMSims}, at redshifts $0.5 \leq z\leq6$, using halo catalogs from Hidden-Valley simulations, and painting halos according to mass-redshift-luminosity relations for each line. We show that at $z=1$ on scales $k_{\rm max} \lesssim 0.8 \ {\rm Mpc}^{-1}h$, the model predictions of clustering power (with only two free parameters) are in agreement with the measured power spectrum at better than 5\%. At higher redshift of $z=4.5$, this remarkable agreement extends to smaller scale of $ k_{\rm max} \lesssim 2 \ {\rm Mpc}^{-1}h$. Furthermore, we show that on large scales, the stochastic contributions to CO and CII power spectra are non-Poissonian, with amplitudes reproduced reasonably well by the halo model prescription. Lastly, we assess the performance of the theoretical model of the baryon acoustic oscillations (BAO) and show that hypothetical LIM surveys probing CO lines at $z=1$, that can be deployed within this decade, will be able to make a high significance measurement of the BAO. On a longer time scale, a space-based mission probing [CII] line can uniquely measure the BAO on a wide range of redshifts at an unprecedented precision., Comment: v3 matches the published article

Measurements of the clustering of galaxies in Fourier space, and at low wavenumbers, offer a window into the early Universe via the possible presence of scale dependent bias generated by Primordial Non Gaussianites. On such large scales a Newtonian treatment of density perturbations might not be sufficient to describe the measurements, and a fully relativistic calculation should be employed. The interpretation of the data is thus further complicated by the fact that relativistic effects break statistical homogeneity and isotropy and are potentially divergent in the Infra-Red (IR). In this work we compute for the first time the ensemble average of the most used Fourier space estimator in spectroscopic surveys, including all general relativistic (GR) effects, and allowing for an arbitrary choice of angular and radial selection functions. We show that any observable is free of IR sensitivity once all the GR terms, individually divergent, are taken into account, and that this cancellation is a consequence of the presence of the Weinberg adiabatic mode as a solution to Einstein's equations. We then study the importance of GR effects, including lensing magnification, in the interpretation of the galaxy power spectrum multipoles, finding that they are in general a small, less than ten percent level, correction to the leading redshift space distortions term. This work represents the baseline for future investigations of the interplay between Primordial Non Gaussianities and GR effects on large scales and in Fourier space., Comment: 44 pages, 10 figures

An observational program focused on the high redshift ($2

The long wavelength modes lost to bright foregrounds in the interferometric 21-cm surveys can partially be recovered using a forward modeling approach that exploits the non-linear coupling between small and large scales induced by gravitational evolution. In this work, we build upon this approach by considering how adding external galaxy distribution data can help to fill in these modes. We consider supplementing the 21-cm data at two different redshifts with a spectroscopic sample (good radial resolution but low number density) loosely modeled on DESI-ELG at $z=1$ and a photometric sample (high number density but poor radial resolution) similar to LSST sample at $z=1$ and $z=4$ respectively. We find that both the galaxy samples are able to reconstruct the largest modes better than only using 21-cm data, with the spectroscopic sample performing significantly better than the photometric sample despite much lower number density. We demonstrate the synergies between surveys by showing that the primordial initial density field is reconstructed better with the combination of surveys than using either of them individually. Methodologically, we also explore the importance of smoothing the density field when using bias models to forward model these tracers for reconstruction.

Line intensity mapping (LIM) is an emerging technique with a unique potential to probe a wide range of scales and redshifts. Realizing the full potential of LIM, however, relies on accurate modeling of the signal. We introduce an extended halo model for the power spectrum of intensity fluctuations of CO rotational lines and [CII] fine transition line in real space, modeling nonlinearities in matter fluctuations and biasing relation between the line intensity fluctuations and the underlying dark matter distribution. We also compute the stochastic contributions beyond the Poisson approximation using the halo model framework. To establish the accuracy of the model, we create the first cosmological-scale simulations of CO and [CII] intensity maps, \textsf{MithraLIMSims}, at redshifts $0.5 \leq z\leq6$, using halo catalogs from Hidden-Valley simulations, and painting halos according to mass-redshift-luminosity relations for each line. We show that at $z=1$ on scales $k_{\rm max} \lesssim 0.8 \ {\rm Mpc}^{-1}h$, the model predictions of clustering power (with only two free parameters) are in agreement with the measured power spectrum at better than 5\%. At higher redshift of $z=4.5$, this remarkable agreement extends to smaller scale of $ k_{\rm max} \lesssim 2 \ {\rm Mpc}^{-1}h$. Furthermore, we show that on large scales, the stochastic contributions to CO and CII power spectra are non-Poissonian, with amplitudes reproduced reasonably well by the halo model prescription. Lastly, we assess the performance of the theoretical model of the baryon acoustic oscillations (BAO) and show that hypothetical LIM surveys probing CO lines at $z=1$, that can be deployed within this decade, will be able to make a high significance measurement of the BAO. On a longer time scale, a space-based mission probing [CII] line can uniquely measure the BAO on a wide range of redshifts at an unprecedented precision., Comment: v3 matches the published article

Measurements of the clustering of galaxies in Fourier space, and at low wavenumbers, offer a window into the early Universe via the possible presence of scale dependent bias generated by Primordial Non Gaussianites. On such large scales a Newtonian treatment of density perturbations might not be sufficient to describe the measurements, and a fully relativistic calculation should be employed. The interpretation of the data is thus further complicated by the fact that relativistic effects break statistical homogeneity and isotropy and are potentially divergent in the Infra-Red (IR). In this work we compute for the first time the ensemble average of the most used Fourier space estimator in spectroscopic surveys, including all general relativistic (GR) effects, and allowing for an arbitrary choice of angular and radial selection functions. We show that any observable is free of IR sensitivity once all the GR terms, individually divergent, are taken into account, and that this cancellation is a consequence of the presence of the Weinberg adiabatic mode as a solution to Einstein's equations. We then study the importance of GR effects, including lensing magnification, in the interpretation of the galaxy power spectrum multipoles, finding that they are in general a small, less than ten percent level, correction to the leading redshift space distortions term. This work represents the baseline for future investigations of the interplay between Primordial Non Gaussianities and GR effects on large scales and in Fourier space., Comment: 44 pages, 11 figures

An observational program focused on the high redshift ($2

The long wavelength modes lost to bright foregrounds in the interferometric 21-cm surveys can partially be recovered using a forward modeling approach that exploits the non-linear coupling between small and large scales induced by gravitational evolution. In this work, we build upon this approach by considering how adding external galaxy distribution data can help to fill in these modes. We consider supplementing the 21-cm data at two different redshifts with a spectroscopic sample (good radial resolution but low number density) loosely modeled on DESI-ELG at $z=1$ and a photometric sample (high number density but poor radial resolution) similar to LSST sample at $z=1$ and $z=4$ respectively. We find that both the galaxy samples are able to reconstruct the largest modes better than only using 21-cm data, with the spectroscopic sample performing significantly better than the photometric sample despite much lower number density. We demonstrate the synergies between surveys by showing that the primordial initial density field is reconstructed better with the combination of surveys than using either of them individually. Methodologically, we also explore the importance of smoothing the density field when using bias models to forward model these tracers for reconstruction., Comment: 16 pages, 7 Figures

The QUIJOTE simulations are a set of 44,100 full N-body simulations spanning more than 7000 cosmological models in the {Omega(m), Omega(b), h, n(s), sigma(8), M-nu, w} hyperplane. At a single redshift, the simulations contain more than 8.5 trillion particles over a combined volume of 44,100 (h(-1) Gpc)(3); each simulation follows the evolution of 256(3), 512(3), or 1024(3) particles in a box of 1 h(-1) Gpc length. Billions of dark matter halos and cosmic voids have been identified in the simulations, whose runs required more than 35 million core hours. The QUIJOTE simulations have been designed for two main purposes: (1) to quantify the information content on cosmological observables and (2) to provide enough data to train machine-learning algorithms. In this paper, we describe the simulations and show a few of their applications. We also release the petabyte of data generated, comprising hundreds of thousands of simulation snapshots at multiple redshifts; halo and void catalogs; and millions of summary statistics, such as power spectra, bispectra, correlation functions, marked power spectra, and estimated probability density functions.

Fluctuations with wavelengths larger than the volume of a galaxy survey affect the measurement of the galaxy power spectrum within the survey itself. In the presence of local Primordial Non- Gaussianities (PNG), in addition to the super-sample matter density and tidal fluctuations, the large-scale gravitational potential also induces a modulation of the observed power spectrum. In this work we investigate this modulation by computing for the first time the response of the redshift-space galaxy power spectrum to the presence of a long wavelength gravitational potential, fully accounting for the stochastic contributions. For biased tracers new response functions arise due to couplings between the small-scale fluctuations in the density, velocity and gravitational fields, the latter through scale dependent bias operators, and the large-scale gravitational potential. We study the impact of the super-sample modes on the measurement of the amplitude of the primordial bispectrum of the local-shape, $f_{\rm NL}^{\rm loc}$, accounting for modulations of both the signal and the covariance of the galaxy power spectrum by the long modes. Considering DESI-like survey specifications, we show that in most cases super-sample modes cause little or no degradation of the constraints, and could actually reduce the errorbars on $f_{\rm NL}^{\rm loc}$ by (10 - 30)\%, if external information on the bias parameters is available., Comment: 15 pages, 4 figures

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

The time delay measured between the images of gravitationally lensed quasars probes a combination of the angular diameter distance to the source-lens system and the mass density profile of the lens. Observational campaigns to measure such systems have reported a determination of the Hubble parameter H0 that shows significant tension with independent determination based on the cosmic microwave background (CMB) and large scale structure (LSS). We show that lens mass models that exhibit a cored component, coexisting with a cusp, probe a degenerate direction in the lens model parameter space, being an approximate mass sheet transformation. This family of lens models has not been considered by the cosmographic analyses. Once added to the model, the cosmographic error budget should become dependent on stellar kinematics uncertainties. We propose that a core component coexisting with a cusp could bring the lensing measurements of H0 to accord with the CMB/LSS value., Comment: 11 pages, 4 figures. v3: added comparison to stellar kinematics constraints, fixed typo in Fig.3, dialed down reference to DM with respect to the core structure; match journal version

The time delay measured between the images of gravitationally lensed quasars probes a combination of the angular diameter distance to the source-lens system and the mass density profile of the lens. Observational campaigns to measure such systems have reported a determination of the Hubble parameter H0 that shows significant tension with independent determination based on the cosmic microwave background (CMB) and large scale structure (LSS). We show that lens mass models that exhibit a cored component, coexisting with a cusp, probe a degenerate direction in the lens model parameter space, being an approximate mass sheet transformation. This family of lens models has not been considered by the cosmographic analyses. Once added to the model, the cosmographic error budget should become dependent on stellar kinematics uncertainties. We propose that a core component coexisting with a cusp could bring the lensing measurements of H0 to accord with the CMB/LSS value., Comment: 11 pages, 4 figures. v3: added comparison to stellar kinematics constraints, fixed typo in Fig.3, dialed down reference to DM with respect to the core structure; match journal version

We present the one-loop 2-point function of biased tracers in redshift space computed with Lagrangian perturbation theory, including a full resummation of both long-wavelength (infrared) displacements and associated velocities. The resulting model accurately predicts the power spectrum and correlation function of halos and mock galaxies from two different sets of N-body simulations at the percent level for quasi-linear scales, including the damping of the baryon acoustic oscillation signal due to the bulk motions of galaxies. We compare this full resummation with other, approximate, techniques including the moment expansion and Gaussian streaming model. We discuss infrared resummation in detail and compare our Lagrangian formulation with the Eulerian theory augmented by an infrared resummation based on splitting the input power spectrum into "wiggle" and "no-wiggle" components. We show that our model is able to recover unbiased cosmological parameters in mock data encompassing a volume much larger than what will be available to future galaxy surveys. We demonstrate how to efficiently compute the resulting expressions numerically, making available a fast Python code capable of rapidly computing these statistics in both configuration and Fourier space., Comment: 37 pages, 10 figures, updated to match version accepted by JCAP

Fluctuations with wavelengths larger than the volume of a galaxy survey affect the measurement of the galaxy power spectrum within the survey itself. In the presence of local Primordial Non- Gaussianities (PNG), in addition to the super-sample matter density and tidal fluctuations, the large-scale gravitational potential also induces a modulation of the observed power spectrum. In this work we investigate this modulation by computing for the first time the response of the redshift-space galaxy power spectrum to the presence of a long wavelength gravitational potential, fully accounting for the stochastic contributions. For biased tracers new response functions arise due to couplings between the small-scale fluctuations in the density, velocity and gravitational fields, the latter through scale dependent bias operators, and the large-scale gravitational potential. We study the impact of the super-sample modes on the measurement of the amplitude of the primordial bispectrum of the local-shape, $f_{\rm NL}^{\rm loc}$, accounting for modulations of both the signal and the covariance of the galaxy power spectrum by the long modes. Considering DESI-like survey specifications, we show that in most cases super-sample modes cause little or no degradation of the constraints, and could actually reduce the errorbars on $f_{\rm NL}^{\rm loc}$ by (10 - 30)\%, if external information on the bias parameters is available., Comment: 15 pages, 4 figures

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

The time delay measured between the images of gravitationally lensed quasars probes a combination of the angular diameter distance to the source-lens system and the mass density profile of the lens. Observational campaigns to measure such systems have reported a determination of the Hubble parameter H0 that shows significant tension with independent determination based on the cosmic microwave background (CMB) and large scale structure (LSS). We show that lens mass models that exhibit a cored component, coexisting with a cusp, probe a degenerate direction in the lens model parameter space, being an approximate mass sheet transformation. This family of lens models has not been considered by the cosmographic analyses. Once added to the model, the cosmographic error budget should become dependent on stellar kinematics uncertainties. We propose that a core component coexisting with a cusp could bring the lensing measurements of H0 to accord with the CMB/LSS value., Comment: 11 pages, 4 figures. v3: added comparison to stellar kinematics constraints, fixed typo in Fig.3, dialed down reference to DM with respect to the core structure; match journal version

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale., Comment: 5 pages + references; Submitted to the Astro2020 call for science white papers. This version: fixed author list

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

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.