scholarly journals ShapeFit: extracting the power spectrum shape information in galaxy surveys beyond BAO and RSD

2021 ◽  
Vol 2021 (12) ◽  
pp. 054
Author(s):  
Samuel Brieden ◽  
Héctor Gil-Marín ◽  
Licia Verde

Abstract In the standard (classic) approach, galaxy clustering measurements from spectroscopic surveys are compressed into baryon acoustic oscillations and redshift space distortions measurements, which in turn can be compared to cosmological models. Recent works have shown that avoiding this intermediate step and fitting directly the full power spectrum signal (full modelling) leads to much tighter constraints on cosmological parameters. Here we show where this extra information is coming from and extend the classic approach with one additional effective parameter, such that it captures, effectively, the same amount of information as the full modelling approach, but in a model-independent way. We validate this new method (ShapeFit) on mock catalogs, and compare its performance to the full modelling approach finding both to deliver equivalent results. The ShapeFit extension of the classic approach promotes the standard analyses at the level of full modelling ones in terms of information content, with the advantages of i) being more model independent; ii) offering an understanding of the origin of the extra cosmological information; iii) allowing a robust control on the impact of observational systematics.

2019 ◽  
Vol 490 (3) ◽  
pp. 4237-4253 ◽  
Author(s):  
Florent Leclercq ◽  
Wolfgang Enzi ◽  
Jens Jasche ◽  
Alan Heavens

ABSTRACT We propose a new, likelihood-free approach to inferring the primordial matter power spectrum and cosmological parameters from arbitrarily complex forward models of galaxy surveys where all relevant statistics can be determined from numerical simulations, i.e. black boxes. Our approach, which we call simulator expansion for likelihood-free inference (selfi), builds upon approximate Bayesian computation using a novel effective likelihood, and upon the linearization of black-box models around an expansion point. Consequently, we obtain simple ‘filter equations’ for an effective posterior of the primordial power spectrum, and a straightforward scheme for cosmological parameter inference. We demonstrate that the workload is computationally tractable, fixed a priori, and perfectly parallel. As a proof of concept, we apply our framework to a realistic synthetic galaxy survey, with a data model accounting for physical structure formation and incomplete and noisy galaxy observations. In doing so, we show that the use of non-linear numerical models allows the galaxy power spectrum to be safely fitted up to at least kmax = 0.5 h Mpc−1, outperforming state-of-the-art backward-modelling techniques by a factor of ∼5 in the number of modes used. The result is an unbiased inference of the primordial matter power spectrum across the entire range of scales considered, including a high-fidelity reconstruction of baryon acoustic oscillations. It translates into an unbiased and robust inference of cosmological parameters. Our results pave the path towards easy applications of likelihood-free simulation-based inference in cosmology. We have made our code pyselfi and our data products publicly available at http://pyselfi.florent-leclercq.eu.


2020 ◽  
Vol 500 (1) ◽  
pp. 259-271
Author(s):  
Alex Smith ◽  
Arnaud de Mattia ◽  
Etienne Burtin ◽  
Chia-Hsun Chuang ◽  
Cheng Zhao

ABSTRACT Accurate mock catalogues are essential for assessing systematics in the cosmological analysis of large galaxy surveys. Anisotropic two-point clustering measurements from the same simulation show some scatter for different lines of sight (LOS), but are on average equal, due to cosmic variance. This results in scatter in the measured cosmological parameters. We use the OuterRim N-body simulation halo catalogue to investigate this, considering the three simulation axes as LOS. The quadrupole of the two-point statistics is particularly sensitive to changes in the LOS, with subper cent level differences in the velocity distributions resulting in ∼1.5 σ shifts on large scales. Averaging over multiple LOS can reduce the impact of cosmic variance. We derive an expression for the Gaussian cross-correlation between the power spectrum multipole measurements, for any two LOS, including shot noise, and the corresponding reduction in variance in the average measurement. Quadrupole measurements are anticorrelated, and for three orthogonal LOS, the variance on the average measurement is reduced by more than 1/3. We perform a Fisher analysis to predict the corresponding gain in precision on the cosmological parameter measurements, which we compare against a set of 300 extended Baryon Oscillation Spectroscopic Survey emission-line galaxy EZmocks. The gain in fσ8, which measures the growth of structure, is also better than 1/3. Averaging over multiple LOS in future mock challenges will allow the redshift space distortion models to be constrained with the same systematic error, with less than three times the CPU time.


2019 ◽  
Vol 485 (4) ◽  
pp. 5059-5072 ◽  
Author(s):  
Phoebe Upton Sanderbeck ◽  
Vid Iršič ◽  
Matthew McQuinn ◽  
Avery Meiksin

ABSTRACT Spatial fluctuations in ultraviolet backgrounds can subtly modulate the distribution of extragalactic sources, a potential signal and systematic for large-scale structure surveys. While this modulation has been shown to be significant for 3D Ly α forest surveys, its relevance for other large-scale structure probes has been hardly explored, despite being the only astrophysical process that likely can affect clustering measurements on the scales of ≳Mpc. We estimate that the background fluctuations, modulating the amount of H i, have a fractional effect of (0.03–0.3) × (k/[10−2 Mpc−1])−1 on the power spectrum of 21 cm intensity maps at z = 1–3. We find a smaller effect for H α and Ly α intensity mapping surveys of (0.001–0.1) × (k/[10−2 Mpc−1])−1 and even smaller effect for more traditional surveys that correlate the positions of individual H α or Ly α emitters. We also estimate the effect of backgrounds on low-redshift galaxy surveys in general based on a simple model in which background fluctuations modulate the rate halo gas cools, modulating star formation: We estimate a maximum fractional effect on the power of ∼0.01 (k/[10−2 Mpc−1])−1 at z = 1. We compare sizes of these imprints to cosmological parameter benchmarks for the next generation of redshift surveys: We find that ionizing backgrounds could result in a bias on the squeezed triangle non-Gaussianity parameter fNL that can be larger than unity for power spectrum measurements with a SPHEREx-like galaxy survey, and typical values of intensity bias. Marginalizing over a shape of the form k−1PL, where PL is the linear matter power spectrum, removes much of this bias at the cost of ${\approx } 40{{\ \rm per\ cent}}$ larger statistical errors.


2020 ◽  
Vol 496 (3) ◽  
pp. 3862-3869 ◽  
Author(s):  
Anatoly Klypin ◽  
Francisco Prada ◽  
Joyce Byun

ABSTRACT Making cosmological inferences from the observed galaxy clustering requires accurate predictions for the mean clustering statistics and their covariances. Those are affected by cosmic variance – the statistical noise due to the finite number of harmonics. The cosmic variance can be suppressed by fixing the amplitudes of the harmonics instead of drawing them from a Gaussian distribution predicted by the inflation models. Initial realizations also can be generated in pairs with 180○ flipped phases to further reduce the variance. Here, we compare the consequences of using paired-and-fixed versus Gaussian initial conditions on the average dark matter clustering and covariance matrices predicted from N-body simulations. As in previous studies, we find no measurable differences between paired-and-fixed and Gaussian simulations for the average density distribution function, power spectrum, and bispectrum. Yet, the covariances from paired-and-fixed simulations are suppressed in a complicated scale- and redshift-dependent way. The situation is particularly problematic on the scales of Baryon acoustic oscillations where the covariance matrix of the power spectrum is lower by only $\sim 20{{\ \rm per\ cent}}$ compared to the Gaussian realizations, implying that there is not much of a reduction of the cosmic variance. The non-trivial suppression, combined with the fact that paired-and-fixed covariances are noisier than from Gaussian simulations, suggests that there is no path towards obtaining accurate covariance matrices from paired-and-fixed simulations – result, that is theoretically expected and accepted in the field. Because the covariances are crucial for the observational estimates of galaxy clustering statistics and cosmological parameters, paired-and-fixed simulations, though useful for some applications, cannot be used for the production of mock galaxy catalogues.


2021 ◽  
Vol 2021 (12) ◽  
pp. 044
Author(s):  
G. Parimbelli ◽  
G. Scelfo ◽  
S.K. Giri ◽  
A. Schneider ◽  
M. Archidiacono ◽  
...  

Abstract We investigate and quantify the impact of mixed (cold and warm) dark matter models on large-scale structure observables. In this scenario, dark matter comes in two phases, a cold one (CDM) and a warm one (WDM): the presence of the latter causes a suppression in the matter power spectrum which is allowed by current constraints and may be detected in present-day and upcoming surveys. We run a large set of N-body simulations in order to build an efficient and accurate emulator to predict the aforementioned suppression with percent precision over a wide range of values for the WDM mass, Mwdm, and its fraction with respect to the totality of dark matter, fwdm. The suppression in the matter power spectrum is found to be independent of changes in the cosmological parameters at the 2% level for k≲ 10 h/Mpc and z≤ 3.5. In the same ranges, by applying a baryonification procedure on both ΛCDM and CWDM simulations to account for the effect of feedback, we find a similar level of agreement between the two scenarios. We examine the impact that such suppression has on weak lensing and angular galaxy clustering power spectra. Finally, we discuss the impact of mixed dark matter on the shape of the halo mass function and which analytical prescription yields the best agreement with simulations. We provide the reader with an application to galaxy cluster number counts.


2020 ◽  
Vol 501 (1) ◽  
pp. 833-852
Author(s):  
Toshiki Kurita ◽  
Masahiro Takada ◽  
Takahiro Nishimichi ◽  
Ryuichi Takahashi ◽  
Ken Osato ◽  
...  

ABSTRACT We use a suite of N-body simulations to study intrinsic alignments (IA) of halo shapes with the surrounding large-scale structure in the ΛCDM model. For this purpose, we develop a novel method to measure multipole moments of the three-dimensional power spectrum of the E-mode field of halo shapes with the matter/halo distribution, $P_{\delta E}^{(\ell)}(k)$ (or $P^{(\ell)}_{{\rm h}E}$), and those of the auto-power spectrum of the E-mode, $P^{(\ell)}_{EE}(k)$, based on the E/B-mode decomposition. The IA power spectra have non-vanishing amplitudes over the linear to non-linear scales, and the large-scale amplitudes at k ≲ 0.1 h−1 Mpc are related to the matter power spectrum via a constant coefficient (AIA), similar to the linear bias parameter of galaxy or halo density field. We find that the cross- and auto-power spectra PδE and PEE at non-linear scales, k ≳ 0.1 h−1 Mpc, show different k-dependences relative to the matter power spectrum, suggesting a violation of the non-linear alignment model commonly used to model contaminations of cosmic shear signals. The IA power spectra exhibit baryon acoustic oscillations, and vary with halo samples of different masses, redshifts, and cosmological parameters (Ωm, S8). The cumulative signal-to-noise ratio for the IA power spectra is about 60 per cent of that for the halo density power spectrum, where the super-sample covariance is found to give a significant contribution to the total covariance. Thus our results demonstrate that the IA power spectra of galaxy shapes, measured from imaging and spectroscopic surveys for an overlapping area of the sky, can be used to probe the underlying matter power spectrum, the primordial curvature perturbations, and cosmological parameters, in addition to the standard galaxy density power spectrum.


2020 ◽  
Vol 634 ◽  
pp. A74 ◽  
Author(s):  
Fabien Lacasa

As galaxy surveys improve their precision thanks to lower levels of noise and the push toward small, non-linear scales, the need for accurate covariances beyond the classical Gaussian formula becomes more acute. Here I investigate the analytical implementation and impact of non-Gaussian covariance terms that I had previously derived for the galaxy angular power spectrum. Braiding covariance is such an interesting class of such terms and it gets contributions both from in-survey and super-survey modes, the latter proving difficult to calibrate through simulations. I present an approximation for braiding covariance which speeds up the process of numerical computation. I show that including braiding covariance is a necessary condition for including other non-Gaussian terms, namely the in-survey 2-, 3-, and 4-halo covariance. Indeed these terms yield incorrect covariance matrices with negative eigenvalues if considered on their own. I then move to quantify the impact on parameter constraints, with forecasts for a survey with Euclid-like galaxy density and angular scales. Compared with the Gaussian case, braiding and in-survey covariances significantly increase the error bars on cosmological parameters, in particular by 50% for the dark energy equation of state w. The error bars on the halo occupation distribution (HOD) parameters are also affected between 12% and 39%. Accounting for super-sample covariance (SSC) also increases parameter errors, by 90% for w and between 7% and 64% for HOD. In total, non-Gaussianity increases the error bar on w by 120% (between 15% and 80% for other cosmological parameters) and the error bars on HOD parameters between 17% and 85%. Accounting for the 1-halo trispectrum term on top of SSC, as has been done in some current analyses, is not sufficient for capturing the full non-Gaussian impact: braiding and the rest of in-survey covariance have to be accounted for. Finally, I discuss why the inclusion of non-Gaussianity generally eases up parameter degeneracies, making cosmological constraints more robust for astrophysical uncertainties. I released publicly the data and a Python notebook reproducing the results and plots of the article.


Author(s):  
Juan García-Bellido

We test the Higgs dilaton inflation model (HDM) using the latest cosmological datasets, including the cosmic microwave background temperature, polarization and lensing data from the Planck satellite (2015), the BICEP and Keck Array experiments, the type Ia supernovae from the JLA catalogue, the baryon acoustic oscillations from CMASS, LOWZ and 6dF, the weak lensing data from the CFHTLenS survey and the matter power spectrum measurements from the latest SDSS data release. We find that the values of all cosmological parameters allowed by the HDM are well within the Planck satellite (2015) constraints. In particular, we determine , , , and (at 95.5% c.l.). We also place new stringent constraints on the couplings of the HDM, ξ χ <0.00328 and (at 95.5% c.l.). We find that the HDM is only slightly better than the w 0 w a CDM model, with . Given that the HDM has two fewer parameters, we find Bayesian evidence favouring the HDM over the w 0 w a CDM model. We also study the critical Higgs inflation model, taking into account the running of both the self-coupling λ( μ ) and the non-minimal coupling to gravity ξ ( μ ). We find peaks in the curvature power spectrum at scales corresponding to the critical value μ that re-enter during the radiation era and collapse to form a broad distribution of clustered primordial black holes, which could constitute today the main component of dark matter. This article is part of the Theo Murphy meeting issue ‘Higgs cosmology’.


2008 ◽  
Vol 23 (01) ◽  
pp. 25-32 ◽  
Author(s):  
SABINO MATARRESE ◽  
MASSIMO PIETRONI

Renormalization Group techniques, successfully employed in quantum field theory and statistical physics, are applied to study the dynamics of structure formation in the Universe. A semi-analytic approach to the computation of the nonlinear power-spectrum of dark matter density fluctuations is proposed. The method can be applied down to zero redshift and to length scales where perturbation theory fails. Our predictions accurately fit the results of numerical simulations in reproducing the "acoustic oscillations" features of the power spectrum, which will be accurately measured in future galaxy surveys and will provide a probe to distinguish among different dark energy models.


2017 ◽  
Vol 606 ◽  
pp. A104 ◽  
Author(s):  
F. Couchot ◽  
S. Henrot-Versillé ◽  
O. Perdereau ◽  
S. Plaszczynski ◽  
B. Rouillé d’Orfeuil ◽  
...  

When combining cosmological and oscillations results to constrain the neutrino sector, the question of the propagation of systematic uncertainties is often raised. We address this issue in the context of the derivation of an upper bound on the sum of the neutrino masses (Σmν) with recent cosmological data. This work is performed within the ΛCDM model extended to Σmν, for which we advocate the use of three mass-degenerate neutrinos. We focus on the study of systematic uncertainties linked to the foregrounds modelling in cosmological microwave background (CMB) data analysis, and on the impact of the present knowledge of the reionisation optical depth. This is done through the use of different likelihoods built from Planck data. Limits on Σmν are derived with various combinations of data, including the latest baryon acoustic oscillations (BAO) and Type Ia supernovae (SNIa) results. We also discuss the impact of the preference for current CMB data for amplitudes of the gravitational lensing distortions higher than expected within the ΛCDM model, and add the Planck CMB lensing. We then derive a robust upper limit: Σmν< 0.17 eV at 95% CL, including 0.01eV of foreground systematics. We also discuss the neutrino mass repartition and show that today’s data do not allow one to disentangle normal from inverted hierarchy. The impact on the other cosmological parameters is also reported, for different assumptions on the neutrino mass repartition, and different high and low multipole CMB likelihoods.


Sign in / Sign up

Export Citation Format

Share Document