Showing 1–50 of 57 results for author: Coulton, W

The Impact of Non-Gaussian Primordial Tails on Cosmological Observables

Authors: William R. Coulton, Oliver H. E. Philcox, Francisco Villaescusa-Navarro

Abstract: Whilst current observational evidence favors a close-to-Gaussian spectrum of primordial perturbations, there exist many models of the early Universe that predict this distribution to have exponentially enhanced or suppressed tails. In this work, we generate realizations of the primordial potential with non-Gaussian tails via a phenomenological model; these are then evolved numerically to obtain ma… ▽ More Whilst current observational evidence favors a close-to-Gaussian spectrum of primordial perturbations, there exist many models of the early Universe that predict this distribution to have exponentially enhanced or suppressed tails. In this work, we generate realizations of the primordial potential with non-Gaussian tails via a phenomenological model; these are then evolved numerically to obtain maps of the cosmic microwave background (CMB) and large-scale structure (LSS). In the CMB maps, our added non-Gaussianity manifests as a localized enhancement of hot and cold spots, which would be expected to contribute to $N$-point functions up to large $N$. Such models are indirectly constrained by \textit{Planck} trispectrum bounds, which restrict the changes in the temperature fluctuations to $O(10μ\mathrm{K})$. In the late-time Universe, we find that tailed cosmologies lead to a halo mass function enhanced at high masses, as expected. Furthermore, significant scale-dependent bias in the halo-halo and halo-matter power spectrum is also sourced, which arises from the squeezed limit of large $N$-point functions that are implicitly generated through the enhancement of the tails. These results underscore that a detection of scale-dependent bias alone cannot be used to rule out single field inflation, but can be used together with other statistics to probe a wide range of primordial processes. △ Less

Submitted 21 June, 2024; originally announced June 2024.

Comments: Comments welcome!

GANSky -- fast curved sky weak lensing simulations using Generative Adversarial Networks

Authors: Supranta S. Boruah, Pier Fiedorowicz, Rafael Garcia, William R. Coulton, Eduardo Rozo, Giulio Fabbian

Abstract: Extracting non-Gaussian information from the next generation weak lensing surveys will require fast and accurate full-sky simulations. This is difficult to achieve in practice with existing simulation methods: ray-traced $N$-body simulations are computationally expensive, and approximate simulation methods (such as lognormal mocks) are not accurate enough. Here, we present GANSky, an interpretable… ▽ More Extracting non-Gaussian information from the next generation weak lensing surveys will require fast and accurate full-sky simulations. This is difficult to achieve in practice with existing simulation methods: ray-traced $N$-body simulations are computationally expensive, and approximate simulation methods (such as lognormal mocks) are not accurate enough. Here, we present GANSky, an interpretable machine learning method that uses Generative Adversarial Networks (GANs) to produce fast and accurate full-sky tomographic weak lensing maps. The input to our GAN are lognormal maps that approximately describe the late-time convergence field of the Universe. Starting from these lognormal maps, we use GANs to learn how to locally redistribute mass to achieve simulation-quality maps. This can be achieved using remarkably small networks ($\approx 10^3$ parameters). We validate the GAN maps by computing a number of summary statistics in both simulated and GANSky maps. We show that GANSky maps correctly reproduce both the mean and $χ^2$ distribution for several statistics, specifically: the 2-pt function, 1-pt PDF, peak and void counts, and the equilateral, folded and squeezed bispectra. These successes makes GANSky an attractive tool to compute the covariances of these statistics. In addition to being useful for rapidly generating large ensembles of artificial data sets, our method can be used to extract non-Gaussian information from weak lensing data with field-level or simulation-based inference. △ Less

Submitted 9 June, 2024; originally announced June 2024.

Comments: 14 pages, 10 Figures, comments welcome

Denoising Diffusion Delensing Delight: Reconstructing the Non-Gaussian CMB Lensing Potential with Diffusion Models

Authors: Thomas Flöss, William R. Coulton, Adriaan J. Duivenvoorden, Francisco Villaescusa-Navarro, Benjamin D. Wandelt

Abstract: Optimal extraction of cosmological information from observations of the Cosmic Microwave Background critically relies on our ability to accurately undo the distortions caused by weak gravitational lensing. In this work, we demonstrate the use of denoising diffusion models in performing Bayesian lensing reconstruction. We show that score-based generative models can produce accurate, uncorrelated sa… ▽ More Optimal extraction of cosmological information from observations of the Cosmic Microwave Background critically relies on our ability to accurately undo the distortions caused by weak gravitational lensing. In this work, we demonstrate the use of denoising diffusion models in performing Bayesian lensing reconstruction. We show that score-based generative models can produce accurate, uncorrelated samples from the CMB lensing convergence map posterior, given noisy CMB observations. To validate our approach, we compare the samples of our model to those obtained using established Hamiltonian Monte Carlo methods, which assume a Gaussian lensing potential. We then go beyond this assumption of Gaussianity, and train and validate our model on non-Gaussian lensing data, obtained by ray-tracing N-body simulations. We demonstrate that in this case, samples from our model have accurate non-Gaussian statistics beyond the power spectrum. The method provides an avenue towards more efficient and accurate lensing reconstruction, that does not rely on an approximate analytic description of the posterior probability. The reconstructed lensing maps can be used as an unbiased tracer of the matter distribution, and to improve delensing of the CMB, resulting in more precise cosmological parameter inference. △ Less

Submitted 6 June, 2024; v1 submitted 9 May, 2024; originally announced May 2024.

Comments: 12 pages, 10 figures. v2: typo in one of the equations fixed, references added

Signal-preserving CMB component separation with machine learning

Authors: Fiona McCarthy, J. Colin Hill, William R. Coulton, David W. Hogg

Abstract: Analysis of microwave sky signals, such as the cosmic microwave background, often requires component separation with multi-frequency methods, where different signals are isolated by their frequency behaviors. Many so-called "blind" methods, such as the internal linear combination (ILC), make minimal assumptions about the spatial distribution of the signal or contaminants, and only assume knowledge… ▽ More Analysis of microwave sky signals, such as the cosmic microwave background, often requires component separation with multi-frequency methods, where different signals are isolated by their frequency behaviors. Many so-called "blind" methods, such as the internal linear combination (ILC), make minimal assumptions about the spatial distribution of the signal or contaminants, and only assume knowledge of the frequency dependence of the signal. The ILC is a minimum-variance linear combination of the measured frequency maps. In the case of Gaussian, statistically isotropic fields, this is the optimal linear combination, as the variance is the only statistic of interest. However, in many cases the signal we wish to isolate, or the foregrounds we wish to remove, are non-Gaussian and/or statistically anisotropic (in particular for Galactic foregrounds). In such cases, it is possible that machine learning (ML) techniques can be used to exploit the non-Gaussian features of the foregrounds and thereby improve component separation. However, many ML techniques require the use of complex, difficult-to-interpret operations on the data. We propose a hybrid method whereby we train an ML model using only combinations of the data that $\textit{do not contain the signal}$, and combine the resulting ML-predicted foreground estimate with the ILC solution to reduce the error from the ILC. We demonstrate our methods on simulations of extragalactic temperature and Galactic polarization foregrounds, and show that our ML model can exploit non-Gaussian features, such as point sources and spatially-varying spectral indices, to produce lower-variance maps than ILC - eg, reducing the variance of the B-mode residual by factors of up to 5 - while preserving the signal of interest in an unbiased manner. Moreover, we often find improved performance when applying our model to foreground models on which it was not trained. △ Less

Submitted 4 April, 2024; originally announced April 2024.

Comments: 22 pages, 13 figures

Quijote-PNG: Optimizing the summary statistics to measure Primordial non-Gaussianity

Authors: Gabriel Jung, Andrea Ravenni, Michele Liguori, Marco Baldi, William R. Coulton, Francisco Villaescusa-Navarro, Benjamin D. Wandelt

Abstract: We apply a suite of different estimators to the Quijote-PNG halo catalogues to find the best approach to constrain Primordial non-Gaussianity (PNG) at non-linear cosmological scales, up to $k_{\rm max} = 0.5 \, h\,{\rm Mpc}^{-1}$. The set of summary statistics considered in our analysis includes the power spectrum, bispectrum, halo mass function, marked power spectrum, and marked modal bispectrum.… ▽ More We apply a suite of different estimators to the Quijote-PNG halo catalogues to find the best approach to constrain Primordial non-Gaussianity (PNG) at non-linear cosmological scales, up to $k_{\rm max} = 0.5 \, h\,{\rm Mpc}^{-1}$. The set of summary statistics considered in our analysis includes the power spectrum, bispectrum, halo mass function, marked power spectrum, and marked modal bispectrum. Marked statistics are used here for the first time in the context of PNG study. We perform a Fisher analysis to estimate their cosmological information content, showing substantial improvements when marked observables are added to the analysis. Starting from these summaries, we train deep neural networks (NN) to perform likelihood-free inference of cosmological and PNG parameters. We assess the performance of different subsets of summary statistics; in the case of $f_\mathrm{NL}^\mathrm{equil}$, we find that a combination of the power spectrum and a suitable marked power spectrum outperforms the combination of power spectrum and bispectrum, the baseline statistics usually employed in PNG analysis. A minimal pipeline to analyse the statistics we identified can be implemented either with our ML algorithm or via more traditional estimators, if these are deemed more reliable. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: 13 pages, 10 figures

A new "temperature inversion" estimator to detect CMB patchy screening by large-scale structure

Authors: Theo Schutt, Abhishek S. Maniyar, Emmanuel Schaan, William R. Coulton, Nishant Mishra

Abstract: Thomson scattering of cosmic microwave background (CMB) photons imprints various properties of the baryons around galaxies on the CMB. One such imprint, called patchy screening, is a direct probe of the gas density profile around galaxies. It usefully complements the information from the kinematic and thermal Sunyaev-Zel'dovich effects and does not require individual redshifts. In this paper, we d… ▽ More Thomson scattering of cosmic microwave background (CMB) photons imprints various properties of the baryons around galaxies on the CMB. One such imprint, called patchy screening, is a direct probe of the gas density profile around galaxies. It usefully complements the information from the kinematic and thermal Sunyaev-Zel'dovich effects and does not require individual redshifts. In this paper, we derive new estimators of patchy screening called the "temperature inversion" (TI) and "signed" estimators, analogous to the gradient inversion estimator of CMB lensing. Pedagogically, we clarify the relation between these estimators and the standard patchy screening quadratic estimator (QE). The new estimators trade optimality for robustness to biases caused by the dominant CMB lensing and foreground contaminants, allowing the use of smaller angular scales. We perform a simulated analysis to realistically forecast the expected precision of patchy screening measurements from four CMB experiments, ACT, SPT, Simons Observatory (SO) and CMB-S4, cross-correlated with three galaxy samples from BOSS, unWISE and the simulated Rubin LSST Data Challenge 2 catalog. Our results give further confidence in the first detection of this effect from the ACT$\times$unWISE data in the companion paper and show patchy screening will be a powerful observable for future surveys like SO, CMB-S4 and LSST. Implementations of the patchy screening QE and the TI and signed estimators are publicly available in our LensQuEst and ThumbStack software packages, available at https://github.com/EmmanuelSchaan/LensQuEst and https://github.com/EmmanuelSchaan/ThumbStack , respectively. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: 20 pages, 7 figures, submitted to PRD

The Atacama Cosmology Telescope: Detection of Patchy Screening of the Cosmic Microwave Background

Authors: William R. Coulton, Theo Schutt, Abhishek S. Maniyar, Emmanuel Schaan, Rui An, Zachary Atkins, Nicholas Battaglia, J Richard Bond, Erminia Calabrese, Steve K. Choi, Mark J. Devlin, Adriaan J. Duivenvoorden, Jo Dunkley, Simone Ferraro, Vera Gluscevic, J. Colin Hill, Matt Hilton, Adam D. Hincks, Arthur Kosowsky, Darby Kramer, Aleksandra Kusiak, Adrien La Posta, Thibaut Louis, Mathew S. Madhavacheril, Gabriela A. Marques , et al. (15 additional authors not shown)

Abstract: Spatial variations in the cosmic electron density after reionization generate cosmic microwave background anisotropies via Thomson scattering, a process known as the ``patchy screening" effect. In this paper, we propose a new estimator for the patchy screening effect that is designed to mitigate biases from the dominant foreground signals. We use it to measure the cross-correlation between \textit… ▽ More Spatial variations in the cosmic electron density after reionization generate cosmic microwave background anisotropies via Thomson scattering, a process known as the ``patchy screening" effect. In this paper, we propose a new estimator for the patchy screening effect that is designed to mitigate biases from the dominant foreground signals. We use it to measure the cross-correlation between \textit{unWISE} galaxies and patchy screening, the latter measured by the Atacama Cosmology Telescope and \textit{Planck} satellite. We report the first detection of the patchy screening effect, with the statistical significance of the cross-correlation exceeding $7σ$. This measurement directly probes the distribution of electrons around these galaxies and provides strong evidence that gas is more extended than the underlying dark matter. By comparing our measurements to electron profiles extracted from simulations, we demonstrate the power of these observations to constrain galaxy evolution models. Requiring only the 2D positions of objects and no individual redshifts or velocity estimates, this approach is complementary to existing gas probes, such as those based on the kinetic Sunyaev-Zeldovich effect. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: See Schutt et al for a detailed comparison of patchy screening estimators. 17 pages with 8 figures

Taming assembly bias for primordial non-Gaussianity

Authors: Emanuele Fondi, Licia Verde, Francisco Villaescusa-Navarro, Marco Baldi, William R. Coulton, Gabriel Jung, Dionysios Karagiannis, Michele Liguori, Andrea Ravenni, Benjamin D. Wandelt

Abstract: 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 betw… ▽ More 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. △ Less

Submitted 2 February, 2024; v1 submitted 16 November, 2023; originally announced November 2023.

Comments: 30 pages, 13 figures. v2: minor updates to match accepted version

Journal ref: JCAP02(2024)048

Small-scale signatures of primordial non-Gaussianity in k-Nearest Neighbour cumulative distribution functions

Authors: William R. Coulton, Tom Abel, Arka Banerjee

Abstract: Searches for primordial non-Gaussianity in cosmological perturbations are a key means of revealing novel primordial physics. However, robustly extracting signatures of primordial non-Gaussianity from non-linear scales of the late-time Universe is an open problem. In this paper, we apply k-Nearest Neighbor cumulative distribution functions, kNN-CDFs, to the \textsc{quijote-png} simulations to explo… ▽ More Searches for primordial non-Gaussianity in cosmological perturbations are a key means of revealing novel primordial physics. However, robustly extracting signatures of primordial non-Gaussianity from non-linear scales of the late-time Universe is an open problem. In this paper, we apply k-Nearest Neighbor cumulative distribution functions, kNN-CDFs, to the \textsc{quijote-png} simulations to explore the sensitivity of kNN-CDFs to primordial non-Gaussianity. An interesting result is that for halo samples with $M_h<10^{14}$ M$_\odot$/h, the kNN-CDFs respond to \textit{equilateral} PNG in a manner distinct from the other parameters. This persists in the galaxy catalogs in redshift space and can be differentiated from the impact of galaxy modelling, at least within the halo occupation distribution (HOD) framework considered here. kNN-CDFs are related to counts-in-cells and, through mapping a subset of the kNN-CDF measurements into the count-in-cells picture, we show that our results can be modeled analytically. A caveat of the analysis is that we only consider the HOD framework, including assembly bias. It will be interesting to validate these results with other techniques for modeling the galaxy--halo connection, e.g., (hybrid) effective field theory or semi-analytical methods. △ Less

Submitted 26 September, 2023; originally announced September 2023.

Cosmology from weak lensing peaks and minima with Subaru Hyper Suprime-Cam survey first-year data

Authors: Gabriela A. Marques, Jia Liu, Masato Shirasaki, Leander Thiele, Daniela Grandón, Kevin M. Huffenberger, Sihao Cheng, Joachim Harnois-Déraps, Ken Osato, William R. Coulton

Abstract: We present cosmological constraints derived from peak counts, minimum counts, and the angular power spectrum of the Subaru Hyper Suprime-Cam first-year (HSC Y1) weak lensing shear catalog. Weak lensing peak and minimum counts contain non-Gaussian information and hence are complementary to the conventional two-point statistics in constraining cosmology. In this work, we forward-model the three summ… ▽ More We present cosmological constraints derived from peak counts, minimum counts, and the angular power spectrum of the Subaru Hyper Suprime-Cam first-year (HSC Y1) weak lensing shear catalog. Weak lensing peak and minimum counts contain non-Gaussian information and hence are complementary to the conventional two-point statistics in constraining cosmology. In this work, we forward-model the three summary statistics and their dependence on cosmology, using a suite of $N$-body simulations tailored to the HSC Y1 data. We investigate systematic and astrophysical effects including intrinsic alignments, baryon feedback, multiplicative bias, and photometric redshift uncertainties. We mitigate the impact of these systematics by applying cuts on angular scales, smoothing scales, statistic bins, and tomographic redshift bins. By combining peaks, minima, and the power spectrum, assuming a flat-$Λ$CDM model, we obtain $S_{8} \equiv σ_8\sqrt{Ω_m/0.3}= 0.810^{+0.022}_{-0.026}$, a 35\% tighter constraint than that obtained from the angular power spectrum alone. Our results are in agreement with other studies using HSC weak lensing shear data, as well as with Planck 2018 cosmology and recent CMB lensing constraints from the Atacama Cosmology Telescope and the South Pole Telescope. △ Less

Submitted 21 August, 2023; originally announced August 2023.

The Atacama Cosmology Telescope: High-resolution component-separated maps across one-third of the sky

Authors: William R. Coulton, Mathew S. Madhavacheril, Adriaan J. Duivenvoorden, J. Colin Hill, Irene Abril-Cabezas, Peter A. R. Ade, Simone Aiola, Tommy Alford, Mandana Amiri, Stefania Amodeo, Rui An, Zachary Atkins, Jason E. Austermann, Nicholas Battaglia, Elia Stefano Battistelli, James A. Beall, Rachel Bean, Benjamin Beringue, Tanay Bhandarkar, Emily Biermann, Boris Bolliet, J Richard Bond, Hongbo Cai, Erminia Calabrese, Victoria Calafut , et al. (129 additional authors not shown)

Abstract: Observations of the millimeter sky contain valuable information on a number of signals, including the blackbody cosmic microwave background (CMB), Galactic emissions, and the Compton-$y$ distortion due to the thermal Sunyaev-Zel'dovich (tSZ) effect. Extracting new insight into cosmological and astrophysical questions often requires combining multi-wavelength observations to spectrally isolate one… ▽ More Observations of the millimeter sky contain valuable information on a number of signals, including the blackbody cosmic microwave background (CMB), Galactic emissions, and the Compton-$y$ distortion due to the thermal Sunyaev-Zel'dovich (tSZ) effect. Extracting new insight into cosmological and astrophysical questions often requires combining multi-wavelength observations to spectrally isolate one component. In this work, we present a new arcminute-resolution Compton-$y$ map, which traces out the line-of-sight-integrated electron pressure, as well as maps of the CMB in intensity and E-mode polarization, across a third of the sky (around 13,000 sq.~deg.). We produce these through a joint analysis of data from the Atacama Cosmology Telescope (ACT) Data Release 4 and 6 at frequencies of roughly 93, 148, and 225 GHz, together with data from the \textit{Planck} satellite at frequencies between 30 GHz and 545 GHz. We present detailed verification of an internal linear combination pipeline implemented in a needlet frame that allows us to efficiently suppress Galactic contamination and account for spatial variations in the ACT instrument noise. These maps provide a significant advance, in noise levels and resolution, over the existing \textit{Planck} component-separated maps and will enable a host of science goals including studies of cluster and galaxy astrophysics, inferences of the cosmic velocity field, primordial non-Gaussianity searches, and gravitational lensing reconstruction of the CMB. △ Less

Submitted 3 July, 2023; originally announced July 2023.

Comments: The Compton-y map and associated products will be made publicly available upon publication of the paper. The CMB T and E mode maps will be made available when the DR6 maps are made public

Signatures of a Parity-Violating Universe

Authors: William R. Coulton, Oliver H. E. Philcox, Francisco Villaescusa-Navarro

Abstract: What would a parity-violating universe look like? We present a numerical and theoretical study of mirror asymmetries in the late universe, using a new suite of $N$-body simulations: QUIJOTE-Odd. These feature parity-violating initial conditions, injected via a simple ansatz for the imaginary primordial trispectrum and evolved into the non-linear regime. We find that the realization-averaged power… ▽ More What would a parity-violating universe look like? We present a numerical and theoretical study of mirror asymmetries in the late universe, using a new suite of $N$-body simulations: QUIJOTE-Odd. These feature parity-violating initial conditions, injected via a simple ansatz for the imaginary primordial trispectrum and evolved into the non-linear regime. We find that the realization-averaged power spectrum, bispectrum, halo mass function, and matter PDF are not affected by our modifications to the initial conditions, deep into the non-linear regime, which we argue arises from rotational and translational invariance. In contrast, the parity-odd trispectrum of matter (measured using a new estimator), shows distinct signatures proportional to the parity-violating parameter, $p_{\rm NL}$, which sets the amplitude of the primordial trispectrum. We additionally find intriguing signatures in the angular momentum of halos, with the primordial trispectrum inducing a non-zero correlation between angular momentum and smoothed velocity field, proportional to $p_{\rm NL}$. Our simulation suite has been made public to facilitate future analyses. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Comments: 19 pages, 9 figures. Simulations available at https://quijote-simulations.readthedocs.io/en/latest/odd.html

Quijote-PNG: The Information Content of the Halo Mass Function

Authors: Gabriel Jung, Andrea Ravenni, Marco Baldi, William R. Coulton, Drew Jamieson, Dionysios Karagiannis, Michele Liguori, Helen Shao, Licia Verde, Francisco Villaescusa-Navarro, Benjamin D. Wandelt

Abstract: 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 neura… ▽ More 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_φ$, without relying on any external prior assumption. Our final forecasts for PNG parameters are $Δf_\mathrm{NL}^\mathrm{local} = 40$, $Δf_\mathrm{NL}^\mathrm{equil} = 210$, $Δ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}$. △ Less

Submitted 4 February, 2024; v1 submitted 17 May, 2023; originally announced May 2023.

Comments: 17 pages, 11 figures. v3 (minor caption fix)

Journal ref: Astrophys.J. 957 (2023) 1, 50

How to estimate Fisher information matrices from simulations

Authors: William R. Coulton, Benjamin D. Wandelt

Abstract: The Fisher information matrix is a quantity of fundamental importance for information geometry and asymptotic statistics. In practice, it is widely used to quickly estimate the expected information available in a data set and guide experimental design choices. In many modern applications, it is intractable to analytically compute the Fisher information and Monte Carlo methods are used instead. The… ▽ More The Fisher information matrix is a quantity of fundamental importance for information geometry and asymptotic statistics. In practice, it is widely used to quickly estimate the expected information available in a data set and guide experimental design choices. In many modern applications, it is intractable to analytically compute the Fisher information and Monte Carlo methods are used instead. The standard Monte Carlo method produces estimates of the Fisher information that can be biased when the Monte-Carlo noise is non-negligible. Most problematic is noise in the derivatives as this leads to an overestimation of the available constraining power, given by the inverse Fisher information. In this work we find another simple estimate that is oppositely biased and produces an underestimate of the constraining power. This estimator can either be used to give approximate bounds on the parameter constraints or can be combined with the standard estimator to give improved, approximately unbiased estimates. Both the alternative and the combined estimators are asymptotically unbiased so can be also used as a convergence check of the standard approach. We discuss potential limitations of these estimators and provide methods to assess their reliability. These methods accelerate the convergence of Fisher forecasts, as unbiased estimates can be achieved with fewer Monte Carlo samples, and so can be used to reduce the simulated data set size by several orders of magnitude. △ Less

Submitted 3 June, 2023; v1 submitted 15 May, 2023; originally announced May 2023.

Comments: Supporting code available at https://github.com/wcoulton/CompressedFisher

The Atacama Cosmology Telescope: DR6 Gravitational Lensing Map and Cosmological Parameters

Authors: Mathew S. Madhavacheril, Frank J. Qu, Blake D. Sherwin, Niall MacCrann, Yaqiong Li, Irene Abril-Cabezas, Peter A. R. Ade, Simone Aiola, Tommy Alford, Mandana Amiri, Stefania Amodeo, Rui An, Zachary Atkins, Jason E. Austermann, Nicholas Battaglia, Elia Stefano Battistelli, James A. Beall, Rachel Bean, Benjamin Beringue, Tanay Bhandarkar, Emily Biermann, Boris Bolliet, J Richard Bond, Hongbo Cai, Erminia Calabrese , et al. (134 additional authors not shown)

Abstract: We present cosmological constraints from a gravitational lensing mass map covering 9400 sq. deg. reconstructed from CMB measurements made by the Atacama Cosmology Telescope (ACT) from 2017 to 2021. In combination with BAO measurements (from SDSS and 6dF), we obtain the amplitude of matter fluctuations $σ_8 = 0.819 \pm 0.015$ at 1.8% precision, $S_8\equivσ_8({Ω_{\rm m}}/0.3)^{0.5}=0.840\pm0.028$ an… ▽ More We present cosmological constraints from a gravitational lensing mass map covering 9400 sq. deg. reconstructed from CMB measurements made by the Atacama Cosmology Telescope (ACT) from 2017 to 2021. In combination with BAO measurements (from SDSS and 6dF), we obtain the amplitude of matter fluctuations $σ_8 = 0.819 \pm 0.015$ at 1.8% precision, $S_8\equivσ_8({Ω_{\rm m}}/0.3)^{0.5}=0.840\pm0.028$ and the Hubble constant $H_0= (68.3 \pm 1.1)\, \text{km}\,\text{s}^{-1}\,\text{Mpc}^{-1}$ at 1.6% precision. A joint constraint with CMB lensing measured by the Planck satellite yields even more precise values: $σ_8 = 0.812 \pm 0.013$, $S_8\equivσ_8({Ω_{\rm m}}/0.3)^{0.5}=0.831\pm0.023$ and $H_0= (68.1 \pm 1.0)\, \text{km}\,\text{s}^{-1}\,\text{Mpc}^{-1}$. These measurements agree well with $Λ$CDM-model extrapolations from the CMB anisotropies measured by Planck. To compare these constraints to those from the KiDS, DES, and HSC galaxy surveys, we revisit those data sets with a uniform set of assumptions, and find $S_8$ from all three surveys are lower than that from ACT+Planck lensing by varying levels ranging from 1.7-2.1$σ$. These results motivate further measurements and comparison, not just between the CMB anisotropies and galaxy lensing, but also between CMB lensing probing $z\sim 0.5-5$ on mostly-linear scales and galaxy lensing at $z\sim 0.5$ on smaller scales. We combine our CMB lensing measurements with CMB anisotropies to constrain extensions of $Λ$CDM, limiting the sum of the neutrino masses to $\sum m_ν < 0.12$ eV (95% c.l.), for example. Our results provide independent confirmation that the universe is spatially flat, conforms with general relativity, and is described remarkably well by the $Λ$CDM model, while paving a promising path for neutrino physics with gravitational lensing from upcoming ground-based CMB surveys. △ Less

Submitted 11 April, 2023; originally announced April 2023.

Comments: 30 pages, 16 figures, prepared for submission to ApJ. Cosmological likelihood data is here: https://lambda.gsfc.nasa.gov/product/act/actadv_prod_table.html ; likelihood software is here: https://github.com/ACTCollaboration/act_dr6_lenslike . Also see companion papers Qu et al and MacCrann et al. Mass maps will be released when papers are published

The Atacama Cosmology Telescope: A Measurement of the DR6 CMB Lensing Power Spectrum and its Implications for Structure Growth

Authors: Frank J. Qu, Blake D. Sherwin, Mathew S. Madhavacheril, Dongwon Han, Kevin T. Crowley, Irene Abril-Cabezas, Peter A. R. Ade, Simone Aiola, Tommy Alford, Mandana Amiri, Stefania Amodeo, Rui An, Zachary Atkins, Jason E. Austermann, Nicholas Battaglia, Elia Stefano Battistelli, James A. Beall, Rachel Bean, Benjamin Beringue, Tanay Bhandarkar, Emily Biermann, Boris Bolliet, J Richard Bond, Hongbo Cai, Erminia Calabrese , et al. (133 additional authors not shown)

Abstract: We present new measurements of cosmic microwave background (CMB) lensing over $9400$ sq. deg. of the sky. These lensing measurements are derived from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) CMB dataset, which consists of five seasons of ACT CMB temperature and polarization observations. We determine the amplitude of the CMB lensing power spectrum at $2.3\%$ precision ($43σ$ sign… ▽ More We present new measurements of cosmic microwave background (CMB) lensing over $9400$ sq. deg. of the sky. These lensing measurements are derived from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) CMB dataset, which consists of five seasons of ACT CMB temperature and polarization observations. We determine the amplitude of the CMB lensing power spectrum at $2.3\%$ precision ($43σ$ significance) using a novel pipeline that minimizes sensitivity to foregrounds and to noise properties. To ensure our results are robust, we analyze an extensive set of null tests, consistency tests, and systematic error estimates and employ a blinded analysis framework. The baseline spectrum is well fit by a lensing amplitude of $A_{\mathrm{lens}}=1.013\pm0.023$ relative to the Planck 2018 CMB power spectra best-fit $Λ$CDM model and $A_{\mathrm{lens}}=1.005\pm0.023$ relative to the $\text{ACT DR4} + \text{WMAP}$ best-fit model. From our lensing power spectrum measurement, we derive constraints on the parameter combination $S^{\mathrm{CMBL}}_8 \equiv σ_8 \left({Ω_m}/{0.3}\right)^{0.25}$ of $S^{\mathrm{CMBL}}_8= 0.818\pm0.022$ from ACT DR6 CMB lensing alone and $S^{\mathrm{CMBL}}_8= 0.813\pm0.018$ when combining ACT DR6 and Planck NPIPE CMB lensing power spectra. These results are in excellent agreement with $Λ$CDM model constraints from Planck or $\text{ACT DR4} + \text{WMAP}$ CMB power spectrum measurements. Our lensing measurements from redshifts $z\sim0.5$--$5$ are thus fully consistent with $Λ$CDM structure growth predictions based on CMB anisotropies probing primarily $z\sim1100$. We find no evidence for a suppression of the amplitude of cosmic structure at low redshifts △ Less

Submitted 28 May, 2024; v1 submitted 11 April, 2023; originally announced April 2023.

Comments: 45+22 pages, 50 figures. v2 matches with published version in ApJ. Cosmological likelihood data and lensing maps are here: https://lambda.gsfc.nasa.gov/product/act/actadv_prod_table.html ; likelihood software is here: https://github.com/ACTCollaboration/act_dr6_lenslike . Also see companion papers Madhavacheril et al and MacCrann et al

Report number: FERMILAB-PUB-23-237-PPD

The Atacama Cosmology Telescope: Mitigating the impact of extragalactic foregrounds for the DR6 CMB lensing analysis

Authors: Niall MacCrann, Blake D. Sherwin, Frank J. Qu, Toshiya Namikawa, Mathew S. Madhavacheril, Irene Abril-Cabezas, Rui An, Jason E. Austermann, Nicholas Battaglia, Elia S. Battistelli, James A. Beall, Boris Bolliet, J. Richard Bond, Hongbo Cai, Erminia Calabrese, William R. Coulton, Omar Darwish, Shannon M. Duff, Adriaan J. Duivenvoorden, Jo Dunkley, Gerrit S. Farren, Simone Ferraro, Joseph E. Golec, Yilun Guan, Dongwon Han , et al. (25 additional authors not shown)

Abstract: We investigate the impact and mitigation of extragalactic foregrounds for the CMB lensing power spectrum analysis of Atacama Cosmology Telescope (ACT) data release 6 (DR6) data. Two independent microwave sky simulations are used to test a range of mitigation strategies. We demonstrate that finding and then subtracting point sources, finding and then subtracting models of clusters, and using a prof… ▽ More We investigate the impact and mitigation of extragalactic foregrounds for the CMB lensing power spectrum analysis of Atacama Cosmology Telescope (ACT) data release 6 (DR6) data. Two independent microwave sky simulations are used to test a range of mitigation strategies. We demonstrate that finding and then subtracting point sources, finding and then subtracting models of clusters, and using a profile bias-hardened lensing estimator, together reduce the fractional biases to well below statistical uncertainties, with the inferred lensing amplitude, $A_{\mathrm{lens}}$, biased by less than $0.2σ$. We also show that another method where a model for the cosmic infrared background (CIB) contribution is deprojected and high frequency data from Planck is included has similar performance. Other frequency-cleaned options do not perform as well, incurring either a large noise cost, or resulting in biased recovery of the lensing spectrum. In addition to these simulation-based tests, we also present null tests performed on the ACT DR6 data which test for sensitivity of our lensing spectrum estimation to differences in foreground levels between the two ACT frequencies used, while nulling the CMB lensing signal. These tests pass whether the nulling is performed at the map or bandpower level. The CIB-deprojected measurement performed on the DR6 data is consistent with our baseline measurement, implying contamination from the CIB is unlikely to significantly bias the DR6 lensing spectrum. This collection of tests gives confidence that the ACT DR6 lensing measurements and cosmological constraints presented in companion papers to this work are robust to extragalactic foregrounds. △ Less

Submitted 11 April, 2023; originally announced April 2023.

Comments: Companion paper to Qu et al and Madhavacheril et al

The Atacama Cosmology Telescope: Map-Based Noise Simulations for DR6

Authors: Zachary Atkins, Adriaan J. Duivenvoorden, William R. Coulton, Frank J. Qu, Simone Aiola, Erminia Calabrese, Grace E. Chesmore, Steve K. Choi, Mark J. Devlin, Jo Dunkley, Carlos Hervías-Caimapo, Yilun Guan, Adrien La Posta, Zack Li, Thibaut Louis, Mathew S. Madhavacheril, Kavilan Moodley, Sigurd Naess, Federico Nati, Michael D. Niemack, Lyman Page, Roberto Puddu, Maria Salatino, Cristóbal Sifón, Suzanne T. Staggs , et al. (3 additional authors not shown)

Abstract: The increasing statistical power of cosmic microwave background (CMB) datasets requires a commensurate effort in understanding their noise properties. The noise in maps from ground-based instruments is dominated by large-scale correlations, which poses a modeling challenge. This paper develops novel models of the complex noise covariance structure in the Atacama Cosmology Telescope Data Release 6… ▽ More The increasing statistical power of cosmic microwave background (CMB) datasets requires a commensurate effort in understanding their noise properties. The noise in maps from ground-based instruments is dominated by large-scale correlations, which poses a modeling challenge. This paper develops novel models of the complex noise covariance structure in the Atacama Cosmology Telescope Data Release 6 (ACT DR6) maps. We first enumerate the noise properties that arise from the combination of the atmosphere and the ACT scan strategy. We then prescribe a class of Gaussian, map-based noise models, including a new wavelet-based approach that uses directional wavelet kernels for modeling correlated instrumental noise. The models are empirical, whose only inputs are a small number of independent realizations of the same region of sky. We evaluate the performance of these models against the ACT DR6 data by drawing ensembles of noise realizations. Applying these simulations to the ACT DR6 power spectrum pipeline reveals a $\sim 20\%$ excess in the covariance matrix diagonal when compared to an analytic expression that assumes noise properties are uniquely described by their power spectrum. Along with our public code, $\mathtt{mnms}$, this work establishes a necessary element in the science pipelines of both ACT DR6 and future ground-based CMB experiments such as the Simons Observatory (SO). △ Less

Submitted 7 March, 2023; originally announced March 2023.

Comments: 41 pages (+10 appendix), 22 figures (+5 appendix), submitted to JCAP

Constraining cosmological vorticity modes with CMB secondary anisotropies

Authors: William R. Coulton, Kazuyuki Akitsu, Masahiro Takada

Abstract: Observational searches for large-scale vorticity modes in the late time Universe are underexplored. Within the standard $Λ$CDM model, this is well motivated given the observed properties of the cosmic microwave background (CMB). However, this means that searches for cosmic vorticity modes can serve as a powerful consistency test of our cosmological model. We show that through combining CMB measure… ▽ More Observational searches for large-scale vorticity modes in the late time Universe are underexplored. Within the standard $Λ$CDM model, this is well motivated given the observed properties of the cosmic microwave background (CMB). However, this means that searches for cosmic vorticity modes can serve as a powerful consistency test of our cosmological model. We show that through combining CMB measurements of the kinetic Sunyaev-Zel'dovich and the moving lens effects with galaxy survey data we can constrain vorticity fields independently from the large scale cosmic velocity field. This approach can provide stringent constraints on the largest scale modes and can be achieved by a simple change in the standard estimators. Alternatively if one assumes there are no cosmic vorticity modes, this estimator can be used to test for systematic biases in existing analyses of kinetic Sunyaev-Zel'dovich effect in a manner analogous to curl-lensing. △ Less

Submitted 26 January, 2023; originally announced January 2023.

Comments: 6 pages, 1 figure

Extragalactic CO emission lines in the CMB experiments: a forgotten signal and a foreground

Authors: Abhishek S. Maniyar, Athanasia Gkogkou, William R. Coulton, Zack Li, Guilaine Lagache, Anthony R. Pullen

Abstract: High resolution cosmic microwave background (CMB) experiments have allowed us to precisely measure the CMB temperature power spectrum down to very small scales (multipole $\ell \sim 3000$). Such measurements at multiple frequencies enable separating the primary CMB anisotropies with other signals like CMB lensing, thermal and kinematic Sunyaev-Zel'dovich effects (tSZ and kSZ), and cosmic infrared… ▽ More High resolution cosmic microwave background (CMB) experiments have allowed us to precisely measure the CMB temperature power spectrum down to very small scales (multipole $\ell \sim 3000$). Such measurements at multiple frequencies enable separating the primary CMB anisotropies with other signals like CMB lensing, thermal and kinematic Sunyaev-Zel'dovich effects (tSZ and kSZ), and cosmic infrared background (CIB). In this paper, we explore another signal of interest at these frequencies that should be present in the CMB maps: extragalactic CO molecular rotational line emissions, which are the most widely used tracers of molecular gas in the line intensity mapping experiments. Using the SIDES simulations adopted for top hat bandpasses at 150 and 220 GHz, we show that the cross-correlation of the CIB with CO lines has a contribution similar to the CIB-tSZ correlation and the kSZ power, thereby contributing a non-negligible amount to the total power at these scales. This signal, therefore, may significantly impact the recently reported $\geq 3σ$ detection of the kSZ power spectrum from the South Pole Telescope (SPT) collaboration, as the contribution of the CO lines is not considered in such analyses. Our results also provide a new way of measuring the CO power spectrum in cross-correlation with the CIB. Finally, these results show that the CO emissions present in the CMB maps will have to be accounted for in all the CMB auto-power spectrum and cross-correlation studies involving a LSS tracer. △ Less

Submitted 25 January, 2023; originally announced January 2023.

Comments: 8 pages main text and bibliography + 2 pages appendices, 6 figures. Comments very welcome

Quijote-PNG: Quasi-maximum likelihood estimation of Primordial Non-Gaussianity in the non-linear halo density field

Authors: Gabriel Jung, Dionysios Karagiannis, Michele Liguori, Marco Baldi, William R Coulton, Drew Jamieson, Licia Verde, Francisco Villaescusa-Navarro, Benjamin D. Wandelt

Abstract: 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 unbia… ▽ More 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 $Δf_\mathrm{NL}^\mathrm{local} = 45$, $Δf_\mathrm{NL}^\mathrm{equil} = 570$, $Δ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}$. △ Less

Submitted 18 May, 2023; v1 submitted 14 November, 2022; originally announced November 2022.

Comments: 17 pages, 10 figures. v2: minor updates to match published version

Journal ref: Astrophys.J. 948 (2023) 2, 135

Snowmass Theory Frontier: Astrophysics and Cosmology

Authors: Daniel Green, Joshua T. Ruderman, Benjamin R. Safdi, Jessie Shelton, Ana Achúcarro, Peter Adshead, Yashar Akrami, Masha Baryakhtar, Daniel Baumann, Asher Berlin, Nikita Blinov, Kimberly K. Boddy, Malte Buschmann, Giovanni Cabass, Robert Caldwell, Emanuele Castorina, Thomas Y. Chen, Xingang Chen, William Coulton, Djuna Croon, Yanou Cui, David Curtin, Francis-Yan Cyr-Racine, Christopher Dessert, Keith R. Dienes , et al. (62 additional authors not shown)

Abstract: 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. 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. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Comments: 57 pages

Biases to primordial non-Gaussianity measurements from CMB secondary anisotropies

Authors: William Coulton, Alexander Miranthis, Anthony Challinor

Abstract: Our view of the last-scattering surface in the cosmic microwave background (CMB) is obscured by secondary anisotropies, sourced by scattering, extragalactic emission and gravitational processes between recombination and observation. Whilst it is established that non-Gaussianity from the correlation between the integrated-Sachs-Wolfe (ISW) effect and gravitational lensing can significantly bias pri… ▽ More Our view of the last-scattering surface in the cosmic microwave background (CMB) is obscured by secondary anisotropies, sourced by scattering, extragalactic emission and gravitational processes between recombination and observation. Whilst it is established that non-Gaussianity from the correlation between the integrated-Sachs-Wolfe (ISW) effect and gravitational lensing can significantly bias primordial non-Gaussianity (PNG) searches, recent work by Hill (2018) has suggested that other combinations of secondary anisotropies can also produce significant biases. Building on that work, we use the WebSky and Sehgal et al.(2010) simulations to perform an extensive examination of possible biases to PNG measurements for the local, equilateral and orthogonal shapes. For a Planck-like CMB experiment, without foreground cleaning, we find significant biases from cosmic infrared background (CIB)-lensing and thermal Sunyaev-Zel'dovich (tSZ)-lensing bispectra for the local and orthogonal templates, and from CIB and tSZ bispectra for the equilateral template. For future experiments, such as the Simons Observatory, biases from correlations between the ISW effect and the tSZ and CIB will also become important. Finally we investigate whether foreground-cleaning techniques are able to suppress these biases sufficiently. We find that the majority of these biases are effectively suppressed by the internal-linear-combination method and the total bias to Planck-like and SO-like experiments is less than the $1\,σ$ statistical error. However, the small total bias arises from the cancellation of several $1\,σ$ biases for Planck-like experiments and $2\,σ$ biases for SO-like. As this cancellation is likely sensitive to the precise modelling, to ensure robustness against these biases explicit removal methods should be used, likely at the cost of decreased constraining power. △ Less

Submitted 25 August, 2022; originally announced August 2022.

Quijote PNG: The information content of the halo power spectrum and bispectrum

Authors: William R Coulton, Francisco Villaescusa-Navarro, Drew Jamieson, Marco Baldi, Gabriel Jung, Dionysios Karagiannis, Michele Liguori, Licia Verde, Benjamin D. Wandelt

Abstract: 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… ▽ More 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 $σ_8$ and our bias model. △ Less

Submitted 20 December, 2022; v1 submitted 30 June, 2022; originally announced June 2022.

Comments: Updated to accepted version

Quijote-PNG: Quasi-maximum likelihood estimation of Primordial Non-Gaussianity in the non-linear dark matter density field

Authors: Gabriel Jung, Dionysios Karagiannis, Michele Liguori, Marco Baldi, William R Coulton, Drew Jamieson, Licia Verde, Francisco Villaescusa-Navarro, Benjamin D. Wandelt

Abstract: 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… ▽ More 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 $Λ$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$-$σ$ bounds of $Δf_{\rm NL}^{\rm local} \sim 16$, $Δf_{\rm NL}^{\rm equil} \sim 77$ and $Δ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. △ Less

Submitted 3 June, 2022; originally announced June 2022.

Comments: 22 pages, 12 figures

Journal ref: Astrophys.J. 940 (2022) 1, 71

Quijote-PNG: Simulations of primordial non-Gaussianity and the information content of the matter field power spectrum and bispectrum

Authors: William R Coulton, Francisco Villaescusa-Navarro, Drew Jamieson, Marco Baldi, Gabriel Jung, Dionysios Karagiannis, Michele Liguori, Licia Verde, Benjamin D. Wandelt

Abstract: 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 hel… ▽ More 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 $Λ$CDM parameters. △ Less

Submitted 26 May, 2023; v1 submitted 3 June, 2022; originally announced June 2022.

Comments: The simulation products are available at https://quijote-simulations.readthedocs.io/en/latest/png.html

An exploration of the properties of cluster profiles for the thermal and kinetic Sunyaev-Zel'dovich effects

Authors: Billy K. K. Lee, William R. Coulton, Leander Thiele, Shirley Ho

Abstract: With the advent of high-resolution, low-noise CMB measurements, the ability to extract cosmological information from thermal Sunyaev-Zel'dovich effect and kinetic Sunyaev-Zel'dovich effect will be limited not by statistical uncertainties but rather by systematic and theoretical uncertainties. The theoretical uncertainty is driven by the lack of knowledge about the electron pressure and density. Th… ▽ More With the advent of high-resolution, low-noise CMB measurements, the ability to extract cosmological information from thermal Sunyaev-Zel'dovich effect and kinetic Sunyaev-Zel'dovich effect will be limited not by statistical uncertainties but rather by systematic and theoretical uncertainties. The theoretical uncertainty is driven by the lack of knowledge about the electron pressure and density. Thus we explore the electron pressure and density distributions in the IllustrisTNG hydrodynamical simulations, and we demonstrate that the cluster properties exhibit a strong dependence on the halo concentration -- providing some of the first evidence of cluster assembly bias in the electron pressure and density. Further, our work shows evidence for a broken power-law mass dependence, with lower pressure in lower mass halos than previous work and a strong evolution with mass of the radial correlations in the electron density and pressure. Both of these effects highlight the differing impact of active galactic nuclei and supernova feedback on the gas in galaxy groups compared to massive clusters. We verified that we see qualitatively similar features in the SIMBA hydro-dynamical simulations, suggesting these effects could be generic features. Finally, we provide a parametric formula for the electron pressure and density profile as a function of dark matter halo mass, halo concentration, and redshift. These fitting formulae can reproduce the distribution of density and pressure of clusters and will be useful in extracting cosmological information from upcoming CMB surveys. △ Less

Submitted 3 May, 2022; originally announced May 2022.

Comments: 18 pages, 7 figures

Inflation: Theory and Observations

Authors: Ana Achúcarro, Matteo Biagetti, Matteo Braglia, Giovanni Cabass, Robert Caldwell, Emanuele Castorina, Xingang Chen, William Coulton, Raphael Flauger, Jacopo Fumagalli, Mikhail M. Ivanov, Hayden Lee, Azadeh Maleknejad, P. Daniel Meerburg, Azadeh Moradinezhad Dizgah, Gonzalo A. Palma, Guilherme L. Pimentel, Sébastien Renaux-Petel, Benjamin Wallisch, Benjamin D. Wandelt, Lukas T. Witkowski, W. L. Kimmy Wu

Abstract: 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… ▽ More 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. △ Less

Submitted 29 September, 2022; v1 submitted 15 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021; 103 pages, 8 figures, 378 endorsers; v2: minor changes

Snowmass 2021 CMB-S4 White Paper

Authors: Kevork Abazajian, Arwa Abdulghafour, Graeme E. Addison, Peter Adshead, Zeeshan Ahmed, Marco Ajello, Daniel Akerib, Steven W. Allen, David Alonso, Marcelo Alvarez, Mustafa A. Amin, Mandana Amiri, Adam Anderson, Behzad Ansarinejad, Melanie Archipley, Kam S. Arnold, Matt Ashby, Han Aung, Carlo Baccigalupi, Carina Baker, Abhishek Bakshi, Debbie Bard, Denis Barkats, Darcy Barron, Peter S. Barry , et al. (331 additional authors not shown)

Abstract: This Snowmass 2021 White Paper describes the Cosmic Microwave Background Stage 4 project CMB-S4, which is designed to cross critical thresholds in our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. We provide an overview of the science case, the technical design, and project plan. This Snowmass 2021 White Paper describes the Cosmic Microwave Background Stage 4 project CMB-S4, which is designed to cross critical thresholds in our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. We provide an overview of the science case, the technical design, and project plan. △ Less

Submitted 15 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021. arXiv admin note: substantial text overlap with arXiv:1908.01062, arXiv:1907.04473

Snowmass2021 Cosmic Frontier: Cosmic Microwave Background Measurements White Paper

Authors: Clarence L. Chang, Kevin M. Huffenberger, Bradford A. Benson, Federico Bianchini, Jens Chluba, Jacques Delabrouille, Raphael Flauger, Shaul Hanany, William C. Jones, Alan J. Kogut, Jeffrey J. McMahon, Joel Meyers, Neelima Sehgal, Sara M. Simon, Caterina Umilta, Kevork N. Abazajian, Zeeshan Ahmed, Yashar Akrami, Adam J. Anderson, Behzad Ansarinejad, Jason Austermann, Carlo Baccigalupi, Denis Barkats, Darcy Barron, Peter S. Barry , et al. (107 additional authors not shown)

Abstract: This is a solicited whitepaper for the Snowmass 2021 community planning exercise. The paper focuses on measurements and science with the Cosmic Microwave Background (CMB). The CMB is foundational to our understanding of modern physics and continues to be a powerful tool driving our understanding of cosmology and particle physics. In this paper, we outline the broad and unique impact of CMB science… ▽ More This is a solicited whitepaper for the Snowmass 2021 community planning exercise. The paper focuses on measurements and science with the Cosmic Microwave Background (CMB). The CMB is foundational to our understanding of modern physics and continues to be a powerful tool driving our understanding of cosmology and particle physics. In this paper, we outline the broad and unique impact of CMB science for the High Energy Cosmic Frontier in the upcoming decade. We also describe the progression of ground-based CMB experiments, which shows that the community is prepared to develop the key capabilities and facilities needed to achieve these transformative CMB measurements. △ Less

Submitted 15 March, 2022; originally announced March 2022.

Comments: contribution to Snowmass 2021

Predicting the Thermal Sunyaev-Zel'dovich Field using Modular and Equivariant Set-Based Neural Networks

Authors: Leander Thiele, Miles Cranmer, William Coulton, Shirley Ho, David N. Spergel

Abstract: Theoretical uncertainty limits our ability to extract cosmological information from baryonic fields such as the thermal Sunyaev-Zel'dovich (tSZ) effect. Being sourced by the electron pressure field, the tSZ effect depends on baryonic physics that is usually modeled by expensive hydrodynamic simulations. We train neural networks on the IllustrisTNG-300 cosmological simulation to predict the continu… ▽ More Theoretical uncertainty limits our ability to extract cosmological information from baryonic fields such as the thermal Sunyaev-Zel'dovich (tSZ) effect. Being sourced by the electron pressure field, the tSZ effect depends on baryonic physics that is usually modeled by expensive hydrodynamic simulations. We train neural networks on the IllustrisTNG-300 cosmological simulation to predict the continuous electron pressure field in galaxy clusters from gravity-only simulations. Modeling clusters is challenging for neural networks as most of the gas pressure is concentrated in a handful of voxels and even the largest hydrodynamical simulations contain only a few hundred clusters that can be used for training. Instead of conventional convolutional neural net (CNN) architectures, we choose to employ a rotationally equivariant DeepSets architecture to operate directly on the set of dark matter particles. We argue that set-based architectures provide distinct advantages over CNNs. For example, we can enforce exact rotational and permutation equivariance, incorporate existing knowledge on the tSZ field, and work with sparse fields as are standard in cosmology. We compose our architecture with separate, physically meaningful modules, making it amenable to interpretation. For example, we can separately study the influence of local and cluster-scale environment, determine that cluster triaxiality has negligible impact, and train a module that corrects for mis-centering. Our model improves by 70 % on analytic profiles fit to the same simulation data. We argue that the electron pressure field, viewed as a function of a gravity-only simulation, has inherent stochasticity, and model this property through a conditional-VAE extension to the network. This modification yields further improvement by 7 %, it is limited by our small training set however. (abridged) △ Less

Submitted 10 June, 2022; v1 submitted 28 February, 2022; originally announced March 2022.

Comments: 11 pages, 5 figures; condensed version accepted at the Fourth Workshop on Machine Learning and the Physical Sciences (NeurIPS 2021) as "Equivariant and Modular DeepSets with Applications in Cluster Cosmology"; revised version accepted by MLST includes numerous clarifications and a new appendix comparing to CNN

The CAMELS project: public data release

Authors: Francisco Villaescusa-Navarro, Shy Genel, Daniel Anglés-Alcázar, Lucia A. Perez, Pablo Villanueva-Domingo, Digvijay Wadekar, Helen Shao, Faizan G. Mohammad, Sultan Hassan, Emily Moser, Erwin T. Lau, Luis Fernando Machado Poletti Valle, Andrina Nicola, Leander Thiele, Yongseok Jo, Oliver H. E. Philcox, Benjamin D. Oppenheimer, Megan Tillman, ChangHoon Hahn, Neerav Kaushal, Alice Pisani, Matthew Gebhardt, Ana Maria Delgado, Joyce Caliendo, Christina Kreisch , et al. (22 additional authors not shown)

Abstract: The Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project was developed to combine cosmology with astrophysics through thousands of cosmological hydrodynamic simulations and machine learning. CAMELS contains 4,233 cosmological simulations, 2,049 N-body and 2,184 state-of-the-art hydrodynamic simulations that sample a vast volume in parameter space. In this paper we present… ▽ More The Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project was developed to combine cosmology with astrophysics through thousands of cosmological hydrodynamic simulations and machine learning. CAMELS contains 4,233 cosmological simulations, 2,049 N-body and 2,184 state-of-the-art hydrodynamic simulations that sample a vast volume in parameter space. In this paper we present the CAMELS public data release, describing the characteristics of the CAMELS simulations and a variety of data products generated from them, including halo, subhalo, galaxy, and void catalogues, power spectra, bispectra, Lyman-$α$ spectra, probability distribution functions, halo radial profiles, and X-rays photon lists. We also release over one thousand catalogues that contain billions of galaxies from CAMELS-SAM: a large collection of N-body simulations that have been combined with the Santa Cruz Semi-Analytic Model. We release all the data, comprising more than 350 terabytes and containing 143,922 snapshots, millions of halos, galaxies and summary statistics. We provide further technical details on how to access, download, read, and process the data at \url{https://camels.readthedocs.io}. △ Less

Submitted 4 January, 2022; originally announced January 2022.

Comments: 18 pages, 3 figures. More than 350 Tb of data from thousands of simulations publicly available at https://www.camel-simulations.org

Effects of boosting on extragalactic components: methods and statistical studies

Authors: William R. Coulton, Sydney Feldman, Karime Maamari, Elena Pierpaoli, Siavash Yasini, Klaus Dolag

Abstract: In this work we examine the impact of our motion with respect to the CMB rest frame on statistics of CMB maps by examining the one-, two-, three- and four- point statistics of simulated maps of the CMB and Sunyaev-Zeldovich (SZ) effects. We validate boosting codes by comparing their outcomes for temperature and polarization power spectra up to $\ell \simeq 6000$. We derive and validate a new analy… ▽ More In this work we examine the impact of our motion with respect to the CMB rest frame on statistics of CMB maps by examining the one-, two-, three- and four- point statistics of simulated maps of the CMB and Sunyaev-Zeldovich (SZ) effects. We validate boosting codes by comparing their outcomes for temperature and polarization power spectra up to $\ell \simeq 6000$. We derive and validate a new analytical formula for the computation of the boosted power spectrum of a signal with a generic frequency dependence. As an example we show how this increases the boosting correction to the power spectrum of CMB intensity measurements by $\sim 30\%$ at 150 GHz. We examine the effect of boosting on thermal and kinetic SZ power spectra from semianalytical and hydrodynamical simulations; the boosting correction is generally small for both simulations, except when considering frequencies near the tSZ null. For the non-Gaussian statistics, in general we find that boosting has no impact with two exceptions. We find that, whilst the statistics of the CMB convergence field are unaffected, quadratic estimators that are used to measure this field can become biased at the $O(1)\%$ level by boosting effects. We present a simple modification to the standard estimators that removes this bias. Second, bispectrum estimators can receive a systematic bias from the Doppler induced quadrupole when there is anisotropy in the sky -- in practice this anisotropy comes from masking and inhomogenous noise. This effect is unobservable and already removed by existing analysis methods. △ Less

Submitted 19 December, 2021; originally announced December 2021.

Comments: 21 pages and 14 figures

Covariant transverse-traceless projection for secondary gravitational waves

Authors: Atsuhisa Ota, Hayley J. Macpherson, William R. Coulton

Abstract: Second-order tensor modes induced by nonlinear gravity are a key component of the cosmological background of gravitational waves. A detection of this background would allow us to probe the primordial power spectrum at otherwise inaccessible scales. Usually, the energy density of these gravitational waves is studied within perturbation theory in a particular gauge -- a connection between our physic… ▽ More Second-order tensor modes induced by nonlinear gravity are a key component of the cosmological background of gravitational waves. A detection of this background would allow us to probe the primordial power spectrum at otherwise inaccessible scales. Usually, the energy density of these gravitational waves is studied within perturbation theory in a particular gauge -- a connection between our physical spacetime and a fictitious background. It is a widely recognized issue that the second-order, scalar-induced gravitational waves are gauge dependent. This issue arises because they are not well-defined as tensors in the physical spacetime at second-order and are thus unphysical. In this paper, we propose the covariant transverse-traceless projection of the extrinsic curvature to study cosmological gravitational waves on a spatial hypersurface. We define a new energy density which is based purely on spacetime tensors, independent of perturbation theory, and thus is gauge invariant by definition. We show that, in the context of second-order perturbation theory, this new energy density contains only propagating modes in the constant-time hypersurface in the Newtonian gauge. We further show that we can recover the same gravitational waves after a transformation to the synchronous gauge, so long as we correctly identify the Newtonian hypersurface. △ Less

Submitted 1 July, 2022; v1 submitted 17 November, 2021; originally announced November 2021.

Comments: 14 pages, 2 figure, major revision

Inference of the optical depth to reionization from low multipole temperature and polarisation Planck data

Authors: Roger de Belsunce, Steven Gratton, William Coulton, George Efstathiou

Abstract: This paper explores methods for constructing low multipole temperature and polarisation likelihoods from maps of the cosmic microwave background anisotropies that have complex noise properties and partial sky coverage. We use Planck 2018 High Frequency Instrument (HFI) and updated SRoll2 temperature and polarisation maps to test our methods. We present three likelihood approximations based on quad… ▽ More This paper explores methods for constructing low multipole temperature and polarisation likelihoods from maps of the cosmic microwave background anisotropies that have complex noise properties and partial sky coverage. We use Planck 2018 High Frequency Instrument (HFI) and updated SRoll2 temperature and polarisation maps to test our methods. We present three likelihood approximations based on quadratic cross spectrum estimators: (i) a variant of the simulation-based likelihood (SimBaL) techniques used in the Planck legacy papers to produce a low multipole EE likelihood; (ii) a semi-analytical likelihood approximation (momento) based on the principle of maximum entropy; (iii) a density-estimation `likelihood-free' scheme (DELFI). Approaches (ii) and (iii) can be generalised to produce low multipole joint temperature-polarisation (TTTEEE) likelihoods. We present extensive tests of these methods on simulations with realistic correlated noise. We then analyse the Planck data and confirm the robustness of our method and likelihoods on multiple inter- and intra-frequency detector set combinations of SRoll2 maps. The three likelihood techniques give consistent results and support a low value of the optical depth to reoinization, tau, from the HFI. Our best estimate of tau comes from combining the low multipole SRoll2 momento (TTTEEE) likelihood with the CamSpec high multipole likelihood and is tau = 0.0627+0.0050-0.0058. This is consistent with the SRoll2 team's determination of tau, though slightly higher by 0.5 sigma, mainly because of our joint treatment of temperature and polarisation. △ Less

Submitted 24 August, 2021; v1 submitted 26 March, 2021; originally announced March 2021.

Comments: 20 pages, 15 figures; published in MNRAS

Journal ref: MNRAS 507, 1072-1091 (2021)

The parity-odd Intrinsic Bispectrum

Authors: William R. Coulton

Abstract: At linear order the only expected source of a curl-like, B mode, polarization pattern in the cosmic microwave background (CMB) is primordial gravitational waves. At second-order B modes are also produced from purely scalar, density, initial conditions. Unlike B modes from primordial gravitational waves, these B modes are expected to be non-Gaussian and not independent from the temperature and grad… ▽ More At linear order the only expected source of a curl-like, B mode, polarization pattern in the cosmic microwave background (CMB) is primordial gravitational waves. At second-order B modes are also produced from purely scalar, density, initial conditions. Unlike B modes from primordial gravitational waves, these B modes are expected to be non-Gaussian and not independent from the temperature and gradient-like polarization, E mode, CMB anisotropies. We find that the three point function between a second-order B mode and two first-order T/E modes is a powerful probe of second-order B modes and should be detectable by upcoming CMB experiments. We focus on the contribution to the three point function arising from non-linear evolution and scattering processes before the end of recombination as this provides new information on the universe at $z> 1000$. We find that this contribution can be separated from the other contributions and is measurable at $\sim 2.5 σ$ by CMB experiments with noise levels of $\sim 1 μ$Karcmin and delensing efficiencies $\ge 90\%$, such as the proposed PICO satellite. We show that approximately half of the total signal arises from non-linearly induced vector and tensor metric perturbations, as evaluated in the Newtonian gauge. This bispectrum is a unique probe of these perturbations in the CMB, as their contribution to the power spectrum is suppressed. An important feature of this bispectrum is that the detectability will increase with decreasing experimental noise, in the absence of primordial B modes, provided that delensing efficiencies improve in parallel. △ Less

Submitted 15 March, 2021; originally announced March 2021.

Comments: 8 pages, 4 figures

Fundamental limits on constraining primordial non-Gaussianity

Authors: Alba Kalaja, P. Daniel Meerburg, Guilherme L. Pimentel, William R. Coulton

Abstract: We study the cosmic variance limit on constraining primordial non-Gaussianity for a variety of theory-motivated shapes. We consider general arguments for 2D and 3D surveys, with a particular emphasis on the CMB. A scale-invariant $N$-point correlator can be measured with a signal-to-noise that naively scales with the square root of the number of observed modes. This intuition generally fails for t… ▽ More We study the cosmic variance limit on constraining primordial non-Gaussianity for a variety of theory-motivated shapes. We consider general arguments for 2D and 3D surveys, with a particular emphasis on the CMB. A scale-invariant $N$-point correlator can be measured with a signal-to-noise that naively scales with the square root of the number of observed modes. This intuition generally fails for two reasons. First, the signal-to-noise scaling is reduced due to the blurring of the last scattering surface at short distances. This blurring is caused by the combination of projection and damping, but the loss of signal is not due to exponential decay, as both signal and noise are equally damped. Second, the behavior of the $N$-point correlator in the squeezed and collapsed (for $N>3$) limits can enhance the scaling of the signal-to-noise with the resolution, even with a reduced range of momenta probing these limits. We provide analytic estimates for all $N$-point correlators. We show that blurring affects equilateral-like shapes much more than squeezed ones. We discuss under what conditions the optimistic scalings in the collapsed limit can be exploited. Lastly, we confirm our analytical estimates with numerical calculations of the signal-to-noise for local, orthogonal and equilateral bispectra, and local trispectra. We also show that adding polarization to intensity data enhances the scaling for equilateral-like spectra. △ Less

Submitted 10 May, 2021; v1 submitted 18 November, 2020; originally announced November 2020.

Comments: 39 pages, 6 figures. Added a paragraph on lensing effect. Matches journal version

The primordial information content of Rayleigh Anisotropies

Authors: William R. Coulton, Benjamin Beringue, P. Daniel Meerburg

Abstract: Anisotropies in the cosmic microwave background (CMB) are primarily generated by Thomson scattering of photons by free electrons. Around recombination, the Thomson scattering probability quickly diminishes as the free electrons combine with protons to form neutral hydrogen off which CMB photons can scatter through Rayleigh scattering. Unlike Thomson scattering, Rayleigh scattering is frequency dep… ▽ More Anisotropies in the cosmic microwave background (CMB) are primarily generated by Thomson scattering of photons by free electrons. Around recombination, the Thomson scattering probability quickly diminishes as the free electrons combine with protons to form neutral hydrogen off which CMB photons can scatter through Rayleigh scattering. Unlike Thomson scattering, Rayleigh scattering is frequency dependent resulting in the generation of anisotropies with a different spectral dependence. Unfortunately the Rayleigh scattering efficiency rapidly decreases with the expansion of the neutral universe, with the result that only a small percentage of photons are scattered by neutral hydrogen. Although the effect is very small, future CMB missions with higher sensitivity and improved frequency coverage are poised to measure Rayleigh scattering signal. The uncorrelated component of the Rayleigh anisotropies contains unique information on the primordial perturbations that could potentially be leveraged to expand our knowledge of the early universe. In this paper we explore whether measurements of Rayleigh scattering anisotropies can be used to constrain primordial non-Gaussianity (NG) and examine the hints of anomalies found by WMAP and \textit{Planck} satellites. We show that the additional Rayleigh information has the potential to improve primordial NG constraints by $30\%$, or more. Primordial bispectra that are not of the local type benefit the most from these additional scatterings, which we attribute to the different scale dependence of the Rayleigh anisotropies. Unfortunately this different scaling means that Rayleigh measurements can not be used to constrain anomalies or features on large scales. On the other hand, anomalies that may persist to smaller scales, such as the potential power asymmetry seen in WMAP and \textit{Planck}, could be improved by the addition of Rayleigh measurements. △ Less

Submitted 20 October, 2020; originally announced October 2020.

Comments: 12 pages, 5 figures

Journal ref: Phys. Rev. D 103, 043501 (2021)

CMB-S4: Forecasting Constraints on Primordial Gravitational Waves

Authors: CMB-S4 Collaboration, :, Kevork Abazajian, Graeme E. Addison, Peter Adshead, Zeeshan Ahmed, Daniel Akerib, Aamir Ali, Steven W. Allen, David Alonso, Marcelo Alvarez, Mustafa A. Amin, Adam Anderson, Kam S. Arnold, Peter Ashton, Carlo Baccigalupi, Debbie Bard, Denis Barkats, Darcy Barron, Peter S. Barry, James G. Bartlett, Ritoban Basu Thakur, Nicholas Battaglia, Rachel Bean, Chris Bebek , et al. (212 additional authors not shown)

Abstract: CMB-S4---the next-generation ground-based cosmic microwave background (CMB) experiment---is set to significantly advance the sensitivity of CMB measurements and enhance our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. Among the science cases pursued with CMB-S4, the quest for detecting p… ▽ More CMB-S4---the next-generation ground-based cosmic microwave background (CMB) experiment---is set to significantly advance the sensitivity of CMB measurements and enhance our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. Among the science cases pursued with CMB-S4, the quest for detecting primordial gravitational waves is a central driver of the experimental design. This work details the development of a forecasting framework that includes a power-spectrum-based semi-analytic projection tool, targeted explicitly towards optimizing constraints on the tensor-to-scalar ratio, $r$, in the presence of Galactic foregrounds and gravitational lensing of the CMB. This framework is unique in its direct use of information from the achieved performance of current Stage 2--3 CMB experiments to robustly forecast the science reach of upcoming CMB-polarization endeavors. The methodology allows for rapid iteration over experimental configurations and offers a flexible way to optimize the design of future experiments given a desired scientific goal. To form a closed-loop process, we couple this semi-analytic tool with map-based validation studies, which allow for the injection of additional complexity and verification of our forecasts with several independent analysis methods. We document multiple rounds of forecasts for CMB-S4 using this process and the resulting establishment of the current reference design of the primordial gravitational-wave component of the Stage-4 experiment, optimized to achieve our science goals of detecting primordial gravitational waves for $r > 0.003$ at greater than $5σ$, or, in the absence of a detection, of reaching an upper limit of $r < 0.001$ at $95\%$ CL. △ Less

Submitted 27 August, 2020; originally announced August 2020.

Comments: 24 pages, 8 figures, 9 tables, submitted to ApJ. arXiv admin note: text overlap with arXiv:1907.04473

Minimizing gravitational lensing contributions to the primordial bispectrum covariance

Authors: William R. Coulton, P. Daniel Meerburg, David G. Baker, Selim Hotinli, Adriaan J. Duivenvoorden, Alexander van Engelen

Abstract: The next generation of ground-based CMB experiments aim to measure temperature and polarization fluctuations up to $\ell_{\rm max} \approx 5000$ over half of the sky. Combined with Planck data on large scales, this will provide improved constraints on primordial non-Gaussianity. However, the impressive resolution of these experiments will come at a price. Besides signal confusion from galactic for… ▽ More The next generation of ground-based CMB experiments aim to measure temperature and polarization fluctuations up to $\ell_{\rm max} \approx 5000$ over half of the sky. Combined with Planck data on large scales, this will provide improved constraints on primordial non-Gaussianity. However, the impressive resolution of these experiments will come at a price. Besides signal confusion from galactic foregrounds, extra-galactic foregrounds and late-time gravitational effects, gravitational lensing will introduce large non-Gaussianity that can become the leading contribution to the bispectrum covariance through the connected 4-point function. Here, we compute this effect analytically for the first time on the full sky for both temperature and polarization. We compare our analytical results with those obtained directly from map-based simulations of the CMB sky for several levels of instrumental noise. Of the standard shapes considered in the literature, the local shape is most affected, resulting in a 35\% increase of the estimator standard deviation for an experiment like the Simons Observatory (SO) and a 110\% increase for a cosmic-variance limited experiment, including both temperature and polarization modes up to $\ell_{\rm max} = 3800$. Because of the nature of the lensing 4-point function, the impact on other shapes is reduced while still non negligible for the orthogonal shape. Two possible avenues to reduce the non-Gaussian contribution to the covariance are proposed. First by marginalizing over lensing contributions, such as the ISW-lensing 3pt function in temperature, and second by delensing the CMB. We show the latter method can remove almost all extra covariance, reducing the effect to below $<$5\% for local bispectra. At the same time, delensing would remove signal biases from secondaries induced by lensing, such as ISW-lensing. △ Less

Submitted 16 December, 2019; originally announced December 2019.

Comments: 22 pages, 9 figures. Abstract abridged

Journal ref: Phys. Rev. D 101, 123504 (2020)

Cosmology with the thermal-kinetic Sunyaev-Zel'dovich effect

Authors: William R. Coulton, Atsuhisa Ota, Alexander van Engelen

Abstract: Compton scattering of the cosmic microwave background (CMB) from hot ionized gas produces a range of effects, and the leading order effects are the kinetic and thermal Sunyaev Zel'dovich (kSZ and tSZ) effects. In the near future, CMB surveys will provide the precision to probe beyond the leading order effects. In this work we study the cosmological information content of the next order term which… ▽ More Compton scattering of the cosmic microwave background (CMB) from hot ionized gas produces a range of effects, and the leading order effects are the kinetic and thermal Sunyaev Zel'dovich (kSZ and tSZ) effects. In the near future, CMB surveys will provide the precision to probe beyond the leading order effects. In this work we study the cosmological information content of the next order term which combines the tSZ and kSZ effects, hereafter called the thermal-kinetic Sunyaev Zel'dovich (tkSZ) effect. As the tkSZ effect has the same velocity dependence as the kSZ effect, it will also have many of the useful properties of the kSZ effect. However, it also has its own, unique spectral dependence, which allows it to be isolated from all other CMB signals. We show that with currently-envisioned CMB missions the tkSZ effect can be detected and can be used to reconstruct large scale velocity fields, with no appreciable bias from either the kSZ effect or other extragalactic foregrounds. Furthermore, since the relativistic corrections arise from the well-studied pressure of ionized gas, rather than the gas number density as in the kSZ effect, the degeneracy due to uncertain gas physics will be significantly reduced. Finally, for a very low-noise experiment the tkSZ effect will be measurable at higher precision than the kSZ. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: 8 pages, 3 figures

Journal ref: Phys. Rev. Lett. 125, 111301 (2020)

Weak Lensing Minima and Peaks: Cosmological Constraints and the Impact of Baryons

Authors: William R. Coulton, Jia Liu, Ian G. McCarthy, Ken Osato

Abstract: We present a novel statistic to extract cosmological information in weak lensing data: the lensing minima. We also investigate the effect of baryons on the cosmological constraints from peak and minimum counts. Using the \texttt{MassiveNuS} simulations, we find that lensing minima are sensitive to non-Gaussian cosmological information and are complementary to the lensing power spectrum and peak co… ▽ More We present a novel statistic to extract cosmological information in weak lensing data: the lensing minima. We also investigate the effect of baryons on the cosmological constraints from peak and minimum counts. Using the \texttt{MassiveNuS} simulations, we find that lensing minima are sensitive to non-Gaussian cosmological information and are complementary to the lensing power spectrum and peak counts. For an LSST-like survey, we obtain $95\%$ credible intervals from a combination of lensing minima and peaks that are significantly stronger than from the power spectrum alone, by $44\%$, $11\%$, and $63\%$ for the neutrino mass sum $\sum m_ν$, matter density $Ω_m$, and amplitude of fluctuation $A_s$, respectively. We explore the effect of baryonic processes on lensing minima and peaks using the hydrodynamical simulations \texttt{BAHAMAS} and \texttt{Osato15}. We find that ignoring baryonic effects would lead to strong ($\approx 4 σ$) biases in inferences from peak counts, but negligible ($\approx 0.5 σ$) for minimum counts, suggesting lensing minima are a potentially more robust tool against baryonic effects. Finally, we demonstrate that the biases can in principle be mitigated without significantly degrading cosmological constraints when we model and marginalize the baryonic effects. △ Less

Submitted 9 October, 2019; originally announced October 2019.

Comments: 12 pages and 11 figures

Baryonic effects on the matter bispectrum

Authors: Simon Foreman, William Coulton, Francisco Villaescusa-Navarro, Alexandre Barreira

Abstract: The large-scale clustering of matter is impacted by baryonic physics, particularly AGN feedback. Modelling or mitigating this impact will be essential for making full use of upcoming measurements of cosmic shear and other large-scale structure probes. We study baryonic effects on the matter bispectrum, using measurements from a selection of state-of-the-art hydrodynamical simulations: IllustrisTNG… ▽ More The large-scale clustering of matter is impacted by baryonic physics, particularly AGN feedback. Modelling or mitigating this impact will be essential for making full use of upcoming measurements of cosmic shear and other large-scale structure probes. We study baryonic effects on the matter bispectrum, using measurements from a selection of state-of-the-art hydrodynamical simulations: IllustrisTNG, Illustris, EAGLE, and BAHAMAS. We identify a low-redshift enhancement of the bispectrum, peaking at $k\sim 3h\,{\rm Mpc}^{-1}$, that is present in several simulations, and discuss how it can be associated to the evolving nature of AGN feedback at late times. This enhancement does not appear in the matter power spectrum, and therefore represents a new source of degeneracy breaking between two- and three-point statistics. In addition, we provide physical interpretations for other aspects of these measurements, and make initial comparisons to predictions from perturbation theory, empirical fitting formulas, and the response function formalism. We publicly release our measurements (including estimates of their uncertainty due to sample variance) and bispectrum measurement code as resources for the community. △ Less

Submitted 25 September, 2020; v1 submitted 8 October, 2019; originally announced October 2019.

Comments: 29 pages, 24 figures. v2: MNRAS published version, with references and clarifications added (conclusions unchanged). See Fig. 13 for a summary of the main results. Simulation measurements are available at https://github.com/sjforeman/hydro_bispectrum , and bispectrum measurement code is available at https://github.com/sjforeman/bskit

The Simons Observatory: Astro2020 Decadal Project Whitepaper

Authors: The Simons Observatory Collaboration, Maximilian H. Abitbol, Shunsuke Adachi, Peter Ade, James Aguirre, Zeeshan Ahmed, Simone Aiola, Aamir Ali, David Alonso, Marcelo A. Alvarez, Kam Arnold, Peter Ashton, Zachary Atkins, Jason Austermann, Humna Awan, Carlo Baccigalupi, Taylor Baildon, Anton Baleato Lizancos, Darcy Barron, Nick Battaglia, Richard Battye, Eric Baxter, Andrew Bazarko, James A. Beall, Rachel Bean , et al. (258 additional authors not shown)

Abstract: The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021… ▽ More The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4. △ Less

Submitted 16 July, 2019; originally announced July 2019.

Comments: Astro2020 Decadal Project Whitepaper. arXiv admin note: text overlap with arXiv:1808.07445

Journal ref: Bull. Am. Astron. Soc. 51 (2019) 147

Inflation and Dark Energy from spectroscopy at $z > 2$

Authors: Simone Ferraro, Michael J. Wilson, Muntazir Abidi, David Alonso, Behzad Ansarinejad, Robert Armstrong, Jacobo Asorey, Arturo Avelino, Carlo Baccigalupi, Kevin Bandura, Nicholas Battaglia, Chetan Bavdhankar, José Luis Bernal, Florian Beutler, Matteo Biagetti, Guillermo A. Blanc, Jonathan Blazek, Adam S. Bolton, Julian Borrill, Brenda Frye, Elizabeth Buckley-Geer, Philip Bull, Cliff Burgess, Christian T. Byrnes, Zheng Cai , et al. (118 additional authors not shown)

Abstract: 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… ▽ More 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 $σ(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. △ Less

Submitted 21 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey

Messengers from the Early Universe: Cosmic Neutrinos and Other Light Relics

Authors: Daniel Green, Mustafa A. Amin, Joel Meyers, Benjamin Wallisch, Kevork N. Abazajian, Muntazir Abidi, Peter Adshead, Zeeshan Ahmed, Behzad Ansarinejad, Robert Armstrong, Carlo Baccigalupi, Kevin Bandura, Darcy Barron, Nicholas Battaglia, Daniel Baumann, Keith Bechtol, Charles Bennett, Bradford Benson, Florian Beutler, Colin Bischoff, Lindsey Bleem, J. Richard Bond, Julian Borrill, Elizabeth Buckley-Geer, Cliff Burgess , et al. (114 additional authors not shown)

Abstract: The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic… ▽ More The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: 5 pages + references; 1 figure; science white paper submitted to the Astro2020 decadal survey

Primordial Non-Gaussianity

Authors: P. Daniel Meerburg, Daniel Green, Muntazir Abidi, Mustafa A. Amin, Peter Adshead, Zeeshan Ahmed, David Alonso, Behzad Ansarinejad, Robert Armstrong, Santiago Avila, Carlo Baccigalupi, Tobias Baldauf, Mario Ballardini, Kevin Bandura, Nicola Bartolo, Nicholas Battaglia, Daniel Baumann, Chetan Bavdhankar, José Luis Bernal, Florian Beutler, Matteo Biagetti, Colin Bischoff, Jonathan Blazek, J. Richard Bond, Julian Borrill , et al. (153 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 14 March, 2019; v1 submitted 11 March, 2019; originally announced March 2019.

Comments: 5 pages + references; Submitted to the Astro2020 call for science white papers. This version: fixed author list

The bispectrum of polarized galactic foregrounds

Authors: William R. Coulton, David N. Spergel

Abstract: Understanding the properties of the galactic emission at millimetre wavelengths is important for studies of the cosmic microwave background (CMB). In this work we explore the bispectrum, the harmonic equivalent of the three point function, from galactic dust and synchrotron emission. We investigate these effects across a broad range of frequencies using the synchrotron dominated S-band Polarizatio… ▽ More Understanding the properties of the galactic emission at millimetre wavelengths is important for studies of the cosmic microwave background (CMB). In this work we explore the bispectrum, the harmonic equivalent of the three point function, from galactic dust and synchrotron emission. We investigate these effects across a broad range of frequencies using the synchrotron dominated S-band Polarization All Sky Survey (SPASS) maps at 2.3 GHz, the Planck satellite maps (30-857 GHz) and dust dominated Infrared Astronomical Satellite (IRAS) maps at 3 THz. We measure bispectra of total intensity fields, T, as well as the gradient, E, and curl modes, B, of the polarization field. We find that the synchrotron and galactic dust have strong temperature bispectra with significant contributions in the squeezed limit, which probes the correlations between two small scale modes with a large scale mode. Additionally, we find that the dust also has strong polarised bispectra that also peak in the squeezed configuration. We explore parity odd bispectra, such as BTT bispectra, and find strong parity odd bispectra for the galactic dust notably in BTT, BTE and BEE configurations. After masking bright sources, we find no evidence for polarised synchrotron bispectra and no evidence for cross bispectra between the dust and synchrotron emission. The strong foreground bispectra discussed here need to be carefully controlled to avoid biasing measurements of primordial non-Gaussianity. Finally we use these bispectra tools to test for residual foregrounds in the component separated Planck maps and find no evidence of residual foregrounds. These tools will be useful for characterizing residual foregrounds in component separated maps, particularly for experiments with less frequency coverage than the Planck satellite. △ Less

Submitted 6 October, 2019; v1 submitted 14 January, 2019; originally announced January 2019.

Comments: Revised version

Constraining Neutrino Mass with the Tomographic Weak Lensing Bispectrum

Authors: William R. Coulton, Jia Liu, Mathew S. Madhavacheril, Vanessa Böhm, David N. Spergel

Abstract: We explore the effect of massive neutrinos on the weak lensing shear bispectrum using the Cosmological Massive Neutrino Simulations. We find that the primary effect of massive neutrinos is to suppress the amplitude of the bispectrum with limited effect on the bispectrum shape. The suppression of the bispectrum amplitude is a factor of two greater than the suppression of the small scale power-spect… ▽ More We explore the effect of massive neutrinos on the weak lensing shear bispectrum using the Cosmological Massive Neutrino Simulations. We find that the primary effect of massive neutrinos is to suppress the amplitude of the bispectrum with limited effect on the bispectrum shape. The suppression of the bispectrum amplitude is a factor of two greater than the suppression of the small scale power-spectrum. For an LSST-like weak lensing survey that observes half of the sky with five tomographic redshift bins, we explore the constraining power of the bispectrum on three cosmological parameters: the sum of the neutrino mass $\sum m_ν$, the matter density $Ω_m$ and the amplitude of primordial fluctuations $A_s$. Bispectrum measurements alone provide similar constraints to the power spectrum measurements and combining the two probes leads to significant improvements than using the latter alone. We find that the joint constraints tighten the power spectrum $95\%$ constraints by $\sim 32\%$ for $\sum m_ν$, $13\%$ for $Ω_m$ and $57\%$ for $A_s$ . △ Less

Submitted 4 October, 2018; originally announced October 2018.

Comments: To be submitted to JCAP

Constraining neutrino mass with tomographic weak lensing peak counts

Authors: Zack Li, Jia Liu, José Manuel Zorrilla Matilla, William R. Coulton

Abstract: Massive cosmic neutrinos change the structure formation history by suppressing perturbations on small scales. Weak lensing data from galaxy surveys probe the structure evolution and thereby can be used to constrain the total mass of the three active neutrinos. However, much of the information is at small scales where the dynamics are nonlinear. Traditional approaches with second order statistics t… ▽ More Massive cosmic neutrinos change the structure formation history by suppressing perturbations on small scales. Weak lensing data from galaxy surveys probe the structure evolution and thereby can be used to constrain the total mass of the three active neutrinos. However, much of the information is at small scales where the dynamics are nonlinear. Traditional approaches with second order statistics thus fail to fully extract the information in the lensing field. In this paper, we study constraints on the neutrino mass sum using lensing peak counts, a statistic which captures non-Gaussian information beyond the second order. We use the ray-traced weak lensing mocks from the Cosmological Massive Neutrino Simulations (MassiveNuS), and apply LSST-like noise. We discuss the effects of redshift tomography, multipole cutoff $\ell_{\rm max}$ for the power spectrum, smoothing scale for the peak counts, and constraints from peaks of different heights. We find that combining peak counts with the traditional lensing power spectrum can improve the constraint on neutrino mass sum, $Ω_m$, and $A_s$ by 39%, 32%, and 60%, respectively, over that from the power spectrum alone. △ Less

Submitted 3 October, 2018; originally announced October 2018.

Comments: 12 pages, 13 figures, comments welcome

Journal ref: Phys. Rev. D 99, 063527 (2019)