Large-scale Structure Surveys Research Articles

Quijote-PNG: Quasi-maximum Likelihood Estimation of Primordial Non-Gaussianity in the Nonlinear Dark Matter Density Field

Future large-scale structure surveys are expected to improve current bounds on primordial non-Gaussianity (PNG), with a significant impact on our understanding of early universe physics. The level of such improvements will however strongly depend on the extent to which late-time nonlinearities erase the PNG signal on small scales. In this work, we show how much primordial information remains in the bispectrum of the nonlinear dark matter density field by implementing a new, simulation-based methodology for joint estimation of PNG amplitudes (f NL) and standard Lambda cold dark matter parameters. The estimator is based on optimally compressed statistics, which, for a given input density field, combine power spectrum and modal bispectrum measurements, and numerically evaluate their covariance and their response to changes in cosmological parameters. In this first analysis, we focus on the matter density field, and we train and validate the estimator using a large suite of N-body simulations (Quijote-png), including different types of PNG (local, equilateral, orthogonal). We explicitly test the estimator’s unbiasedness, optimality, and stability with respect to changes in the total number of input realizations. While the dark matter power spectrum itself contains negligible PNG information, as expected, including it as an ancillary statistic increases the PNG information content extracted from the bispectrum by a factor of order 2. As a result, we prove the capability of our approach to optimally extract PNG information on nonlinear scales beyond the perturbative regime, up to . At the same time, we discuss the significant information on cosmological parameters contained on these scales.

Rescuing constraints on modified gravity using gravitational redshift in large-scale structure

The distribution of galaxies provides an ideal laboratory to test for deviations from General Relativity. In particular, redshift-space distortions are commonly used to constrain modifications to the Poisson equation, which governs the strength of dark matter clustering. Here, we show that these constraints rely on the validity of the weak equivalence principle, which has never been tested for the dark matter component. Relaxing this restrictive assumption leads to modifications in the growth of structure that are fully degenerate with modifications induced by the Poisson equation. This in turns strongly degrades the constraining power of redshift-space distortions. Such degeneracies can however be broken and tight constraints on modified gravity can be recovered by measuring gravitational redshift from the galaxy distribution, an effect that will be detectable by the coming generation of large-scale structure surveys.

One spectrum to cure them all: signature from early Universe solves major anomalies and tensions in cosmology

Acoustic peaks in the Cosmic Microwave Background (CMB) temperature spectrum as observed by the Planck satellite appear to be smoother than our expectation from the standard model lensing effect. This anomalous effect can be also mimicked by a spatially closed Universe with a very low value of Hubble constant that consequently aggravates the already existing discordance between cosmological observations. We reconstruct a signature from the early Universe, a particular form of oscillation in the primordial spectrum of quantum fluctuations with a characteristic frequency, that solves all these anomalies. Interestingly, we find this form of the primordial spectrum resolves or substantially subsides, various tensions in the standard model of cosmology in fitting different observations, namely Planck CMB, clustering and weak lensing shear measurements from several large scale structure surveys, local measurements of Hubble constant, and recently estimated age of the Universe from globular clusters. We support our findings phenomenologically, by proposing an analytical form of the primordial spectrum with similar features and demonstrate that it agrees remarkably well with various combinations of cosmological observations. We support further our findings theoretically, by introducing a single scalar field potential for inflation that can generate such a form of the primordial spectrum.

Molecular Chemistry for Dark Matter. II. Recombination, Molecule Formation, and Halo Mass Function in Atomic Dark Matter

Dissipative dark matter predicts rich observable phenomena that can be tested with future large-scale structure surveys. As a specific example, we study atomic dark matter, consisting of a heavy particle and a light particle charged under a dark electromagnetism. In particular, we calculate the cosmological evolution of atomic dark matter focusing on dark recombination and dark molecule formation. We have obtained the relevant interaction rate coefficients by rescaling the rates for normal hydrogen, and evolved the abundances for ionized, atomic, and molecular states using a modified version of Recfast++ (which we have released publicly at a a https://github.com/jamesgurian/RecfastJulia ). We also provide an analytical approximation for the final abundances. We then calculate the effects of atomic dark matter on the linear power spectrum, which enter through a dark photon diffusion and dark acoustic oscillations. At formation time, the atomic dark matter model suppresses halo abundances on scales smaller than the diffusion scale, just as warm dark matter models suppress the abundance below the free-streaming scale. The subsequent evolution with radiative cooling, however, will alter the halo mass function further.

Magrathea-Pathfinder: a 3D adaptive-mesh code for geodesic ray tracing in N-body simulations

We introduce MAGRATHEA-PATHFINDER, a relativistic ray-tracing framework that can reconstruct the past light cone of observers in cosmological simulations. The code directly computes the 3D trajectory of light rays through the null geodesic equations, with the weak-field limit as its only approximation. This approach offers high levels of versatility while removing the need for many of the standard ray-tracing approximations such as plane-parallel, Born, or multiple-lens. Moreover, the use of adaptive integration steps and interpolation strategies based on adaptive-mesh refinement grids allows MAGRATHEA-PATHFINDER to accurately account for the nonlinear regime of structure formation and fully take advantage of the small-scale gravitational clustering. To handle very large N-body simulations, the framework has been designed as a high-performance computing post-processing tool relying on a hybrid paralleliza-tion that combines MPI tasks with C++11 std::threads. In this paper, we describe how realistic cosmological observables can be computed from numerical simulation using ray-tracing techniques. We discuss in particular the production of simulated catalogues and sky maps that account for all the observational effects considering first-order metric perturbations (such as peculiar velocities, gravitational potential, integrated Sachs-Wolfe, time-delay, and gravitational lensing). We perform convergence tests of our gravitational lensing algorithms and conduct performance benchmarks of the null geodesic integration procedures. MAGRATHEA-PATHFINDER introduces sophisticated ray-tracing tools to make the link between the space of N-body simulations and light-cone observables. This should provide new ways of exploring existing cosmological probes and building new ones beyond standard assumptions in order to prepare for the next generation of large-scale structure surveys.

Towards an optimal estimation of cosmological parameters with the wavelet scattering transform

Optimal extraction of the non-Gaussian information encoded in the Large-Scale Structure (LSS) of the universe lies at the forefront of modern precision cosmology. We propose achieving this task through the use of the Wavelet Scattering Transform (WST), which subjects an input field to a layer of non-linear transformations that are sensitive to non-Gaussianity in spatial density distributions through a generated set of WST coefficients. In order to assess its applicability in the context of LSS surveys, we apply the WST on the 3D overdensity field obtained by the Quijote simulations, out of which we extract the Fisher information in 6 cosmological parameters. It is subsequently found to deliver a large improvement in the marginalized errors on all parameters, ranging between $1.2-4\times$ tighter than the corresponding ones obtained from the regular 3D cold dark matter + baryon power spectrum, as well as a $50 \%$ improvement over the neutrino mass constraint given by the marked power spectrum. Through this first application on 3D cosmological fields, we demonstrate the great promise held by this novel statistic and set the stage for its future application to actual galaxy observations.

The integrated Sachs–Wolfe effect in 4D Einstein–Gauss–Bonnet gravity

A novel 4-dimensional Einstein-Gauss-Bonnet (4D EGB) gravity has been proposed that asserts to bypass the Lovelock's theorem and to result in a non-trivial contribution to the gravitational dynamics in four-dimensional spacetime. In this work, we study the integrated Sachs-Wolfe (ISW) effect in the 4D EGB model. For this purpose, we calculate the evolution of the gravitational potential and the linear growth factor as a function of redshift for the 4D EGB model and compare it with the corresponding result obtained from the $\Lambda$-cold dark matter ($\Lambda$CDM) model. We also calculate the ISW-auto power spectrum and the ISW-cross power spectrum as functions of cosmic microwave background (CMB) multipoles for the 4D EGB model and compare those with the one obtained from the $\Lambda$CDM model. To do this, we use the strongest constraint on the coupling parameter proposed for the 4D EGB model. Additionally, to calculate the ISW effect for the 4D EGB model, we employ three large-scale structure surveys from different wavelengths. The results exhibit that the ISW effect in the 4D EGB model is higher than the one obtained from the $\Lambda$CDM model. Hence, we show that the 4D EGB model can amplify the amplitude of the ISW power spectrum, which can be considered as a relative advantage of the 4D EGB model comparing the $\Lambda$CDM one. Also, we indicate that the deviation from the $\Lambda$CDM model is directly proportional to the value of the dimensionless coupling parameter $\beta$.

CosmoPower: emulating cosmological power spectra for accelerated Bayesian inference from next-generation surveys

ABSTRACT We present CosmoPower, a suite of neural cosmological power spectrum emulators providing orders-of-magnitude acceleration for parameter estimation from two-point statistics analyses of Large-Scale Structure (LSS) and Cosmic Microwave Background (CMB) surveys. The emulators replace the computation of matter and CMB power spectra from Boltzmann codes; thus, they do not need to be re-trained for different choices of astrophysical nuisance parameters or redshift distributions. The matter power spectrum emulation error is less than $0.4{{\ \rm per\ cent}}$ in the wavenumber range $k \in [10^{-5}, 10] \, \mathrm{Mpc}^{-1}$ for redshift z ∈ [0, 5]. CosmoPower emulates CMB temperature, polarization, and lensing potential power spectra in the 5-σ region of parameter space around the Planck best-fitting values with an error ${\lesssim}10{{\ \rm per\ cent}}$ of the expected shot noise for the forthcoming Simons Observatory. CosmoPower is showcased on a joint cosmic shear and galaxy clustering analysis from the Kilo-Degree Survey, as well as on a Stage IV Euclid-like simulated cosmic shear analysis. For the CMB case, CosmoPower is tested on a Planck 2018 CMB temperature and polarization analysis. The emulators always recover the fiducial cosmological constraints with differences in the posteriors smaller than sampling noise, while providing a speed-up factor up to O(104) to the complete inference pipeline. This acceleration allows posterior distributions to be recovered in just a few seconds, as we demonstrate in the Planck likelihood case. CosmoPower is written entirely in python, can be interfaced with all commonly used cosmological samplers, and is publicly available at: https://github.com/alessiospuriomancini/cosmopower.

The Three Hundred project: Contrasting clusters galaxy density in hydrodynamical and dark matter only simulations

Cluster number count is a major cosmological probe for the next generation of cosmological large scale-structure surveys like the one expected from the Euclid satellite mission. Cosmological constraints will be mainly limited by the understanding of the selection function (SF), which characterize the probability of detecting a cluster of a given mass and redshift. The SF can be estimated by injecting realistic simulated clusters into the survey and re-applying the detection procedure. For this purpose we intend to use The Three Hundreds project, a 324 cluster sample simulated with full-physics hydrodynamical re-simulations. In this paper we concentrate on the study of the distribution of member galaxies in the cluster sample. First, we study possible resolution effects by comparing low and high resolution simulations. Finally, accounting for the latter we derive the density profiles of the member galaxies and discuss their evolution with cluster mass and redshift.

Deep learning simulations of the microwave sky

We present 500 high-resolution, full-sky millimeter-wave Deep Learning (DL) simulations that include lensed CMB maps and correlated foreground components. We find that these MillimeterDL simulations can reproduce a wide range of non-Gaussian summary statistics matching the input training simulations, while only being optimized to match the power spectra. The procedure we develop in this work enables the capability to mass produce independent full-sky realizations from a single expensive full-sky simulation, when ordinarily the latter would not provide enough training data. We also circumvent a common limitation of high-resolution DL simulations that they be confined to small sky areas, often due to memory or GPU issues; we do this by developing a "stitching" procedure that can faithfully recover the high-order statistics of a full-sky map without discontinuities or repeated features. In addition, since our network takes as input a full-sky lensing convergence map, it can in principle take a full-sky lensing convergence map from any large-scale structure (LSS) simulation and generate the corresponding lensed CMB and correlated foreground components at millimeter wavelengths; this is especially useful in the current era of combining results from both CMB and LSS surveys, which require a common set of simulations.

Magnification and evolution biases in large-scale structure surveys

Measurements of galaxy clustering in upcoming surveys such as those planned for the Euclid and Roman satellites, and the SKA Observatory, will be sensitive to distortions from lensing magnification and Doppler effects, beyond the standard redshift-space distortions. The amplitude of these contributions depends sensitively on magnification bias and evolution bias in the galaxy number density. Magnification bias quantifies the change in the observed number of galaxies gained or lost by lensing magnification, while evolution bias quantifies the physical change in the galaxy number density relative to the conserved case. These biases are given by derivatives of the number density, and consequently are very sensitive to the form of the luminosity function. We give a careful derivation of the magnification and evolution biases, clarifying a number of results in the literature. We then examine the biases for a variety of surveys, encompassinggalaxy surveys and line intensity mapping at radio and optical/near-infrared wavelengths.

Multi-tasking the growth of cosmological structures

Next-generation large-scale structure surveys will deliver a significant increase in the precision of growth data, allowing us to use ‘agnostic’ methods to study the evolution of perturbations without the assumption of a cosmological model. We focus on a particular machine learning tool, Gaussian processes, to reconstruct the growth rate f , the root mean square of matter fluctuations σ 8 , and their product f σ 8 . We apply this method to simulated data, representing the precision of upcoming Stage IV galaxy surveys. We extend the standard single-task approach to a multi-task approach that reconstructs the three functions simultaneously, thereby taking into account their inter-dependence. We find that this multi-task approach outperforms the precision of the single-task approach for future surveys. By contrast, the limited sensitivity of current data severely hinders the use of agnostic methods, since the Gaussian processes parameters need to be fine tuned in order to obtain robust reconstructions.

AbacusHOD: a highly efficient extended multitracer HOD framework and its application to BOSS and eBOSS data

ABSTRACT We introduce the AbacusHOD model and present two applications of AbacusHOD and the AbacusSummit simulations to observations. AbacusHOD is a Halo Occupation Distribution (HOD) framework written in Python that is particle-based, multitracer, highly generalized, and highly efficient. It is designed specifically with multitracer/cosmology analyses for next-generation large-scale structure surveys in mind, and takes advantage of the volume and precision offered by the new state-of-the-art AbacusSummit cosmological simulations. The model is also highly customizable and should be broadly applicable to any upcoming surveys and a diverse range of cosmological analyses. In this paper, we demonstrate the capabilities of the AbacusHOD framework through two example applications. The first example demonstrates the high efficiency and the large HOD extension feature set through an analysis of full-shape redshift-space clustering of BOSS galaxies at intermediate to small scales ($\lt 30\, h^{-1}$ Mpc), assessing the necessity of introducing secondary galaxy biases (assembly bias). We find strong evidence for using halo environment instead of concentration to trace secondary galaxy bias, a result which also leads to a moderate reduction in the ‘lensing is low’ tension. The second example demonstrates the multitracer capabilities of the AbacusHOD package through an analysis of the extended Baryon Oscillation Spectroscopic Survey cross-correlation measurements between three different galaxy tracers: luminous red galaxies, emission-line galaxies, and quasi-stellar objects. We expect the AbacusHOD framework, in combination with the AbacusSummit simulation suite, to play an important role in a simulation-based analysis of the upcoming Dark Energy Spectroscopic Instrument data sets.

Non-Gaussian stochastic gravitational waves from phase transitions

Cosmological phase transitions in the primordial universe can produce anisotropic stochastic gravitational wave backgrounds (GWB), similar to the cosmic microwave background (CMB). For adiabatic perturbations, the fluctuations in GWB follow those in the CMB, but if primordial fluctuations carry an isocurvature component, this need no longer be true. It is shown that in non-minimal inflationary and reheating settings, primordial isocurvature can survive in GWB and exhibit significant non-Gaussianity (NG) in contrast to the CMB, while obeying current observational bounds. While probing such NG GWB is at best a marginal possibility at LISA, there is much greater scope at future proposed detectors such as DECIGO and BBO. It is even possible that the first observations of inflation-era NG could be made with gravitational wave detectors as opposed to the CMB or Large-Scale Structure surveys.

Dark sector to restore cosmological concordance

We develop a new phenomenological model that addresses current tensions between observations of the early and late Universe. Our scenario features: (i) a decaying dark energy fluid (DDE), which undergoes a transition at $z \sim 5,000$, to raise today's value of the Hubble parameter -- addressing the $H_0$ tension, and (ii) an ultra-light axion (ULA), which starts oscillating at $z\gtrsim 10^4$, to suppress the matter power spectrum -- addressing the $S_8$ tension. Our Markov Chain Monte Carlo analyses show that such a Dark Sector model fits a combination of Cosmic Microwave Background (CMB), Baryon Acoustic Oscillations, and Large Scale Structure (LSS) data slightly better than the $\Lambda$CDM model, while importantly reducing both the $H_0$ and $S_8$ tensions with late universe probes ($\lesssim 3\sigma$). Combined with measurements from cosmic shear surveys, we find that the discrepancy on $S_8$ is reduced to the $1.4\sigma$ level, and the value of $H_0$ is further raised. Adding local supernovae measurements, we find that the $H_0$ and $S_8$ tensions are reduced to the $1.4\sigma$ and $1.2\sigma$ level respectively, with a significant improvement $\Delta\chi^2\simeq -18$ compared to the $\Lambda$CDM model. With this complete dataset, the DDE and ULA are detected at $\simeq 4\sigma$ and $\simeq2\sigma$, respectively. We discuss a possible particle physics realization of this model, with a dark confining gauge sector and its associated axion, although embedding the full details within microphysics remains an urgent open question. Our scenario will be decisively probed with future CMB and LSS surveys.

Galaxy number counts at second order in perturbation theory: a leading-order term comparison

The Galaxy number density is a key quantity to compare theoretical predictions to the observational data from current and future large scale structure surveys. The precision demanded by these stage IV surveys requires the use of second order cosmological perturbation theory. Based on the independent calculation published previously, we present the result of the comparison with the results of three other groups at leading order. Overall we find that the differences between the different approaches lie mostly on the definition of certain quantities, where the ambiguity of signs results in the addition of extra terms at second order in perturbation theory.

Baryonic effects in the Effective Field Theory of Large-Scale Structure and an analytic recipe for lensing in CMB-S4

Upcoming Large-Scale Structure surveys will likely become the next leading sources of cosmological information, making it crucial to have a precise understanding of the influence of baryons on cosmological observables. The Effective Field Theory of Large-Scale Structure (EFTofLSS) provides a consistent way to predict the clustering of dark matter and baryons on large scales, where their leading corrections in perturbation theory are given by a simple and calculable functional form even after the onset of baryonic processes. In this paper, we extend the two-fluid-like system up to two-loop order in perturbation theory. Along the way, we show that a new linear counterterm proportional to the relative velocity of the fluids could generically be present, but we show that its effects are expected to be small in our universe. Regardless, we show how to consistently perform perturbation theory in the presence of this new term. We find that the EFTofLSS at two-loop order can accurately account for the details of baryonic processes on large scales. We compare our results to a hydrodynamical N-body simulation at many redshifts and find that the counterterms associated with the leading corrections to dark matter and baryons start to differ between redshifts z ≈ 3 and z ≈ 2, signaling the onset of star-formation physics. We then use these fits to compute the lensing power spectrum, show that the understanding of baryonic processes will be important for analyzing CMB-S4 data, and show that the two-loop EFTofLSS accurately captures these effects for ℓ ≲ 2000. Our results are also potentially of interest for current and future weak lensing surveys.

Measuring anisotropic stress with relativistic effects

One of the main goal of large-scale structure surveys is to test the consistency of General Relativity at cosmological scales. In the $\Lambda$CDM model of cosmology, the relations between the fields describing the geometry and the content of our Universe are uniquely determined. In particular, the two gravitational potentials -- that describe the spatial and temporal fluctuations in the geometry -- are equal. Whereas large classes of dark energy models preserve this equality, theories of modified gravity generally create a difference between the potentials, known as anisotropic stress. Even though measuring this anisotropic stress is one of the key goals of large-scale structure surveys, there are currently no methods able to measure it directly. Current methods all rely on measurements of galaxy peculiar velocities (through redshift-space distortions), from which the time component of the metric is inferred, assuming that dark matter follows geodesics. If this is not the case, all the proposed tests fail to measure the anisotropic stress. In this letter, we propose a novel test which directly measures anisotropic stress, without relying on any assumption about the unknown dark matter. Our method uses relativistic effects in the galaxy number counts to provide a direct measurement of the time component of the metric. By comparing this with lensing observations our test provides a direct measurement of the anisotropic stress.

Dark matter neutrino scattering in the galactic centre with IceCube

While there is evidence for the existence of dark matter, its properties have yet to be discovered. Simultaneously, the nature of high-energy astrophysical neutrinos detected by IceCube remains unresolved. If dark matter and neutrinos are coupled to each other, they may exhibit a non-zero elastic scattering cross section. Such an interaction between an isotropic extragalactic neutrino flux and dark matter would be concentrated in the Galactic Centre, where the dark matter column density is greatest. This scattering would attenuate the flux of high-energy neutrinos, which could be observed in IceCube. Using the seven-year Medium Energy Starting Events, we perform an unbinned likelihood analysis, searching for a signal based on a possible dark matter-neutrino interaction scenario. We search for a suppression of the high-energy astrophysical neutrino flux in the direction of the Galactic Centre, and compare these constraints to complementary low-energy information from large scale structure surveys and the cosmic microwave background.

Does jackknife scale really matter for accurate large-scale structure covariances?

ABSTRACT The jackknife method gives an internal covariance estimate for large-scale structure surveys and allows model-independent errors on cosmological parameters. Using the SDSS-III BOSS CMASS sample, we study how the jackknife size and number of resamplings impact the precision of the covariance estimate on the correlation function multipoles and the error on the inferred baryon acoustic scale. We compare the measurement with the MultiDark Patchy mock galaxy catalogues, and we also validate it against a set of lognormal mocks with the same survey geometry. We build several jackknife configurations that vary in size and number of resamplings. We introduce the Hartlap factor in the covariance estimate that depends on the number of jackknife resamplings. We also find that it is useful to apply the tapering scheme to estimate the precision matrix from a limited number of resamplings. The results from CMASS and mock catalogues show that the error estimate of the baryon acoustic scale does not depend on the jackknife scale. For the shift parameter α, we find an average error of 1.6 per cent, 2.2 per cent and 1.2 per cent, respectively, from CMASS, Patchy, and lognormal jackknife covariances. Despite these uncertainties fluctuate significantly due to some structural limitations of the jackknife method, our α estimates are in reasonable agreement with published pre-reconstruction analyses. Jackknife methods will provide valuable and complementary covariance estimates for future large-scale structure surveys.