-
Studying large-scale structure probes of modified gravity with COLA
Authors:
Bartolomeo Fiorini,
Kazuya Koyama,
Albert Izard
Abstract:
We study the effect of two Modified Gravity (MG) theories, $f(R)$ and nDGP, on three probes of large-scale structure, the real space power spectrum estimator $Q_0$, bispectrum and voids, and validate fast approximate COLA simulations against full $N$-body simulations for the prediction of these probes. We find that using the first three even multipoles of the redshift space power spectrum to estim…
▽ More
We study the effect of two Modified Gravity (MG) theories, $f(R)$ and nDGP, on three probes of large-scale structure, the real space power spectrum estimator $Q_0$, bispectrum and voids, and validate fast approximate COLA simulations against full $N$-body simulations for the prediction of these probes. We find that using the first three even multipoles of the redshift space power spectrum to estimate $Q_0$ is enough to reproduce the MG boost factors of the real space power spectrum for both halo and galaxy catalogues. By analysing the bispectrum and reduced bispectrum of Dark Matter (DM), we show that the strong MG signal present in the DM bispectrum is mainly due to the enhanced power spectrum. We warn about adopting screening approximations in simulations as this neglects non-linear contributions that can source a significant component of the MG bispectrum signal at the DM level, but we argue that this is not a problem for the bispectrum of galaxies in redshift space where the signal is dominated by the non-linear galaxy bias. Finally, we perform void-finding on our galaxy mock catalogues by the ZOBOV watershed algorithm. To apply a linear model for Redshift-Space Distortion (RSD) in the void-galaxy cross-correlation function, we first examine the effects of MG on the void profiles entering into the RSD model. We find relevant MG signals in the integrated-density, velocity dispersion and radial velocity profiles in the nDGP theory. Fitting the RSD model for the linear growth rate, we recover the linear theory prediction in an nDGP model, which is larger than the $Λ$CDM prediction at the $3 σ$ level. In $f(R)$ theory we cannot naively compare the results of the fit with the linear theory prediction as this is scale-dependent, but we obtain results that are consistent with the $Λ$CDM prediction.
△ Less
Submitted 21 November, 2022; v1 submitted 2 August, 2022;
originally announced August 2022.
-
Dark Energy Survey Year 3 Results: Galaxy mock catalogs for BAO analysis
Authors:
I. Ferrero,
M. Crocce,
I. Tutusaus,
A. Porredon,
L. Blot,
P. Fosalba,
A. Carnero Rosell,
S. Avila,
A. Izard,
J. Elvin-Poole,
K. C. Chan,
H. Camacho,
R. Rosenfeld,
E. Sanchez,
P. Tallada-Crespí,
J. Carretero,
I. Sevilla-Noarbe,
E. Gaztanaga,
F. Andrade-Oliveira,
J. De Vicente,
J. Mena-Fernández,
A. J. Ross,
D. Sanchez Cid,
A. Ferté,
A. Brandao-Souza
, et al. (61 additional authors not shown)
Abstract:
The calibration and validation of scientific analysis in simulations is a fundamental tool to ensure unbiased and robust results in observational cosmology. In particular, mock galaxy catalogs are a crucial resource to achieve these goals in the measurement of baryon acoustic oscillation (BAO) in the clustering of galaxies. Here we present a set of 1952 galaxy mock catalogs designed to mimic the D…
▽ More
The calibration and validation of scientific analysis in simulations is a fundamental tool to ensure unbiased and robust results in observational cosmology. In particular, mock galaxy catalogs are a crucial resource to achieve these goals in the measurement of baryon acoustic oscillation (BAO) in the clustering of galaxies. Here we present a set of 1952 galaxy mock catalogs designed to mimic the Dark Energy Survey (DES) Year 3 BAO sample over its full photometric redshift range $0.6<z_{\rm photo}<1.1$. The mocks are based upon 488 ICE-COLA fast $N$-body simulations of full-sky light cones and were created by populating halos with galaxies, using a hybrid halo occupation distribution - halo abundance matching model. This model has ten free parameters, which were determined, for the first time, using an automatic likelihood minimization procedure. We also introduced a novel technique to assign photometric redshift for simulated galaxies, following a two-dimensional probability distribution with VIMOS Public Extragalactic Redshift Survey (VIPERS) data. The calibration was designed to match the observed abundance of galaxies as a function of photometric redshift, the distribution of photometric redshift errors, and the clustering amplitude on scales smaller than those used for BAO measurements. An exhaustive analysis was done to ensure that the mocks reproduce the input properties. Finally, mocks were tested by comparing the angular correlation function $w(θ)$, angular power spectrum $C_\ell$, and projected clustering $ξ_p(r_\perp)$ to theoretical predictions and data. The impact of volume replication in the estimate of the covariance is also investigated. The success in accurately reproducing the photometric redshift uncertainties and the galaxy clustering as a function of redshift render this mock creation pipeline as a benchmark for future analyses of photometric galaxy surveys.
△ Less
Submitted 13 December, 2021; v1 submitted 9 July, 2021;
originally announced July 2021.
-
Fast generation of mock galaxy catalogues in modified gravity models with COLA
Authors:
Bartolomeo Fiorini,
Kazuya Koyama,
Albert Izard,
Hans A. Winther,
Bill S. Wright,
Baojiu Li
Abstract:
We investigate the viability of producing galaxy mock catalogues with COmoving Lagrangian Acceleration (COLA) simulations in Modified Gravity (MG) models employing the Halo Occupation Distribution (HOD) formalism. In this work, we focus on two theories of MG: $f(R)$ gravity with the chameleon mechanism, and a braneworld model (nDGP) that incorporates the Vainshtein mechanism. We use a suite of ful…
▽ More
We investigate the viability of producing galaxy mock catalogues with COmoving Lagrangian Acceleration (COLA) simulations in Modified Gravity (MG) models employing the Halo Occupation Distribution (HOD) formalism. In this work, we focus on two theories of MG: $f(R)$ gravity with the chameleon mechanism, and a braneworld model (nDGP) that incorporates the Vainshtein mechanism. We use a suite of full $N$-body simulations in MG as a benchmark to test the accuracy of COLA simulations. At the level of Dark Matter (DM), we show that COLA accurately reproduces the matter power spectrum up to $k \sim 1 h {\rm Mpc}^{-1}$, while it is less accurate in reproducing the velocity field. To produce halo catalogues, we find that the ROCKSTAR halo-finder does not perform well with COLA simulations. On the other hand, using a simple Friends-of-Friends (FoF) finder and an empirical mass conversion from FoF to spherical over-density masses, we are able to produce halo catalogues in COLA that are in good agreement with those in $N$-body simulations. To consider the effects of the MG fifth force on the halo profile, we derive simple fitting formulae for the concentration-mass and the velocity dispersion-mass relations that we calibrate using ROCKSTAR halo catalogues in $N$-body simulations. We then use these results to extend the HOD formalism to modified gravity simulations in COLA. We use an HOD model with five parameters that we tune to obtain galaxy catalogues in redshift space. We find that despite the great freedom of the HOD model, MG leaves characteristic imprints in the redshift space power spectrum multipoles and these features are well captured by the COLA galaxy catalogues.
△ Less
Submitted 23 August, 2021; v1 submitted 9 June, 2021;
originally announced June 2021.
-
WFIRST: The Essential Cosmology Space Observatory for the Coming Decade
Authors:
O. Doré,
C. Hirata,
Y. Wang,
D. Weinberg,
T. Eifler,
R. J. Foley,
C. He Heinrich,
E. Krause,
S. Perlmutter,
A. Pisani,
D. Scolnic,
D. N. Spergel,
N. Suntzeff,
G. Aldering,
C. Baltay,
P. Capak,
A. Choi,
S. Deustua,
C. Dvorkin,
S. M. Fall,
X. Fang,
A. Fruchter,
L. Galbany,
S. Ho,
R. Hounsell
, et al. (24 additional authors not shown)
Abstract:
Two decades after its discovery, cosmic acceleration remains the most profound mystery in cosmology and arguably in all of physics. Either the Universe is dominated by a form of dark energy with exotic physical properties not predicted by standard model physics, or General Relativity is not an adequate description of gravity over cosmic distances. WFIRST emerged as a top priority of Astro2010 in p…
▽ More
Two decades after its discovery, cosmic acceleration remains the most profound mystery in cosmology and arguably in all of physics. Either the Universe is dominated by a form of dark energy with exotic physical properties not predicted by standard model physics, or General Relativity is not an adequate description of gravity over cosmic distances. WFIRST emerged as a top priority of Astro2010 in part because of its ability to address the mystery of cosmic acceleration through both high precision measurements of the cosmic expansion history and the growth of cosmic structures with multiple and redundant probes. We illustrate in this white paper how mission design changes since Astro2010 have made WFIRST an even more powerful dark energy facility and have improved the ability of WFIRST to respond to changes in the experimental landscape. WFIRST is the space-based probe of DE the community needs in the mid-2020s.
△ Less
Submitted 1 April, 2019;
originally announced April 2019.
-
Comparing approximate methods for mock catalogues and covariance matrices III: Bispectrum
Authors:
Manuel Colavincenzo,
Emiliano Sefusatti,
Pierluigi Monaco,
Linda Blot,
Martin Crocce,
Martha Lippich,
Ariel G. Sánchez,
Marcelo A. Alvarez,
Aniket Agrawal,
Santiago Avila,
Andrés Balaguera-Antolínez,
Richard Bond,
Sandrine Codis,
Claudio Dalla Vecchia,
Antonio Dorta,
Pablo Fosalba,
Albert Izard,
Francisco-Shu Kitaura,
Marcos Pellejero-Ibanez,
George Stein,
Mohammadjavad Vakili,
Gustavo Yepes
Abstract:
We compare the measurements of the bispectrum and the estimate of its covariance obtained from a set of different methods for the efficient generation of approximate dark matter halo catalogs to the same quantities obtained from full N-body simulations. To this purpose we employ a large set of three-hundred realisations of the same cosmology for each method, run with matching initial conditions in…
▽ More
We compare the measurements of the bispectrum and the estimate of its covariance obtained from a set of different methods for the efficient generation of approximate dark matter halo catalogs to the same quantities obtained from full N-body simulations. To this purpose we employ a large set of three-hundred realisations of the same cosmology for each method, run with matching initial conditions in order to reduce the contribution of cosmic variance to the comparison. In addition, we compare how the error on cosmological parameters such as linear and nonlinear bias parameters depends on the approximate method used for the determination of the bispectrum variance. As general result, most methods provide errors within 10% of the errors estimated from N-body simulations. Exceptions are those methods requiring calibration of the clustering amplitude but restrict this to two-point statistics. Finally we test how our results are affected by being limited to a few hundreds measurements from N-body simulation, and therefore to the bispectrum variance, by comparing with a larger set of several thousands realisations performed with one approximate method.
△ Less
Submitted 8 October, 2018; v1 submitted 25 June, 2018;
originally announced June 2018.
-
Comparing approximate methods for mock catalogues and covariance matrices II: Power spectrum multipoles
Authors:
Linda Blot,
Martin Crocce,
Emiliano Sefusatti,
Martha Lippich,
Ariel G. Sánchez,
Manuel Colavincenzo,
Pierluigi Monaco,
Marcelo A. Alvarez,
Aniket Agrawal,
Santiago Avila,
Andrés Balaguera-Antolínez,
Richard Bond,
Sandrine Codis,
Claudio Dalla Vecchia,
Antonio Dorta,
Pablo Fosalba,
Albert Izard,
Francisco-Shu Kitaura,
Marcos Pellejero-Ibanez,
George Stein,
Mohammadjavad Vakili,
Gustavo Yepes
Abstract:
We study the accuracy of several approximate methods for gravitational dynamics in terms of halo power spectrum multipoles and their estimated covariance matrix. We propagate the differences in covariances into parameter constrains related to growth rate of structure, Alcock-Paczynski distortions and biasing. We consider seven methods in three broad categories: algorithms that solve for halo densi…
▽ More
We study the accuracy of several approximate methods for gravitational dynamics in terms of halo power spectrum multipoles and their estimated covariance matrix. We propagate the differences in covariances into parameter constrains related to growth rate of structure, Alcock-Paczynski distortions and biasing. We consider seven methods in three broad categories: algorithms that solve for halo density evolution deterministically using Lagrangian trajectories (ICE-COLA, Pinocchio and PeakPatch), methods that rely on halo assignment schemes onto dark-matter overdensities calibrated with a target N-body run (Halogen, Patchy) and two standard assumptions about the full density PDF (Gaussian and Lognormal). We benchmark their performance against a set of three hundred N-body simulations, running similar sets of approximate simulations with matched initial conditions, for each method. We find that most methods reproduce the monopole to within $5\%$, while residuals for the quadrupole are sometimes larger and scale dependent. The variance of the multipoles is typically reproduced within $10\%$. Overall, we find that covariances built from approximate simulations yield errors on model parameters within $10\%$ of those from the N-body based covariance.
△ Less
Submitted 18 February, 2019; v1 submitted 25 June, 2018;
originally announced June 2018.
-
Comparing approximate methods for mock catalogues and covariance matrices I: correlation function
Authors:
Martha Lippich,
Ariel G. Sánchez,
Manuel Colavincenzo,
Emiliano Sefusatti,
Pierluigi Monaco,
Linda Blot,
Martin Crocce,
Marcelo A. Alvarez,
Aniket Agrawal,
Santiago Avila,
Andrés Balaguera-Antolínez,
Richard Bond,
Sandrine Codis,
Claudio Dalla Vecchia,
Antonio Dorta,
Pablo Fosalba,
Albert Izard,
Francisco-Shu Kitaura,
Marcos Pellejero-Ibanez,
George Stein,
Mohammadjavad Vakili,
Gustavo Yepes
Abstract:
This paper is the first in a set that analyses the covariance matrices of clustering statistics obtained from several approximate methods for gravitational structure formation. We focus here on the covariance matrices of anisotropic two-point correlation function measurements. Our comparison includes seven approximate methods, which can be divided into three categories: predictive methods that fol…
▽ More
This paper is the first in a set that analyses the covariance matrices of clustering statistics obtained from several approximate methods for gravitational structure formation. We focus here on the covariance matrices of anisotropic two-point correlation function measurements. Our comparison includes seven approximate methods, which can be divided into three categories: predictive methods that follow the evolution of the linear density field deterministically (ICE-COLA, Peak Patch, and Pinocchio), methods that require a calibration with N-body simulations (Patchy and Halogen), and simpler recipes based on assumptions regarding the shape of the probability distribution function (PDF) of density fluctuations (log-normal and Gaussian density fields). We analyse the impact of using covariance estimates obtained from these approximate methods on cosmological analyses of galaxy clustering measurements, using as a reference the covariances inferred from a set of full N-body simulations. We find that all approximate methods can accurately recover the mean parameter values inferred using the N-body covariances. The obtained parameter uncertainties typically agree with the corresponding N-body results within 5% for our lower mass threshold, and 10% for our higher mass threshold. Furthermore, we find that the constraints for some methods can differ by up to 20% depending on whether the halo samples used to define the covariance matrices are defined by matching the mass, number density, or clustering amplitude of the parent N-body samples. The results of our configuration-space analysis indicate that most approximate methods provide similar results, with no single method clearly outperforming the others.
△ Less
Submitted 13 May, 2019; v1 submitted 25 June, 2018;
originally announced June 2018.
-
WFIRST Science Investigation Team "Cosmology with the High Latitude Survey" Annual Report 2017
Authors:
Olivier Doré,
Christopher Hirata,
Yun Wang,
David Weinberg,
Ivano Baronchelli,
Andrew Benson,
Peter Capak,
Ami Choi,
Tim Eifler,
Shoubaneh Hemmati,
Shirley Ho,
Albert Izard,
Bhuvnesh Jain,
Mike Jarvis,
Alina Kiessling,
Elisabeth Krause,
Elena Massara,
Dan Masters,
Alex Merson,
Hironao Miyatake,
Andres Plazas Malagon,
Rachel Mandelbaum,
Lado Samushia,
Chaz Shapiro,
Melanie Simet
, et al. (24 additional authors not shown)
Abstract:
Cosmic acceleration is the most surprising cosmological discovery in many decades. Testing and distinguishing among possible explanations requires cosmological measurements of extremely high precision probing the full history of cosmic expansion and structure growth and, ideally, compare and contrast matter and relativistic tracers of the gravity potential. This program is one of the defining obje…
▽ More
Cosmic acceleration is the most surprising cosmological discovery in many decades. Testing and distinguishing among possible explanations requires cosmological measurements of extremely high precision probing the full history of cosmic expansion and structure growth and, ideally, compare and contrast matter and relativistic tracers of the gravity potential. This program is one of the defining objectives of the Wide-Field Infrared Survey Telescope (WFIRST), as set forth in the New Worlds, New Horizons report (NWNH) in 2010. The WFIRST mission has the ability to improve these measurements by 1-2 orders of magnitude compared to the current state of the art, while simultaneously extending their redshift grasp, greatly improving control of systematic effects, and taking a unified approach to multiple probes that provide complementary physical information and cross-checks of cosmological results. We describe in this annual report the activities of the Science Investigation Team (SIT) "Cosmology with the High Latitude Survey (HLS)" during the year 2017. This team was selected by NASA in December 2015 in order to address the stringent challenges of the WFIRST dark energy (DE) program through the Project's formulation phase. This SIT has elected to jointly address Galaxy Redshift Survey, Weak Lensing and Cluster Growth and thus fully embrace the fact that the imaging and spectroscopic elements of the HLS will be realized as an integrated observing program, and they jointly impose requirements on performance and operations. WFIRST is designed to be able to deliver a definitive result on the origin of cosmic acceleration. It is not optimized for Figure of Merit sensitivity but for control of systematic uncertainties and for having multiple techniques each with multiple cross-checks. Our SIT work focuses on understanding the potential systematics in the WFIRST DE measurements.
△ Less
Submitted 10 April, 2018;
originally announced April 2018.
-
ICE-COLA: fast simulations for weak lensing observables
Authors:
Albert Izard,
Pablo Fosalba,
Martin Crocce
Abstract:
Approximate methods to full N-body simulations provide a fast and accurate solution to the development of mock catalogues for the modeling of galaxy clustering observables. In this paper we extend ICE-COLA (Izard et al. 2016), based on an optimized implementation of the approximate COLA method, to produce weak lensing maps and halo catalogues in the light cone using an integrated and self consiste…
▽ More
Approximate methods to full N-body simulations provide a fast and accurate solution to the development of mock catalogues for the modeling of galaxy clustering observables. In this paper we extend ICE-COLA (Izard et al. 2016), based on an optimized implementation of the approximate COLA method, to produce weak lensing maps and halo catalogues in the light cone using an integrated and self consistent approach. We show that despite the approximate dynamics, the catalogues thus produced enable an accurate modeling of weak lensing observables one decade beyond the characteristic scale where the growth becomes non-linear. In particular, we compare ICE-COLA to the MICE-GC N-body simulation for some fiducial cases representative of upcoming surveys and find that, for sources at redshift $z=1$, their convergence power spectra agree to within one percent up to high multipoles (i.e., of order $1000$). The corresponding shear two point functions, $ξ_{+}$ and $ξ_{-}$, yield similar accuracy down to $2$ and $20$ arcmin respectively, while tangential shear around a $z=0.5$ lens sample is accurate down to $4$ arcmin. We show that such accuracy is stable against an increased angular resolution of the weak lensing maps. Hence, this opens the possibility of using approximate methods for the joint modeling of galaxy clustering and weak lensing observables and their covariance in ongoing and future galaxy surveys.
△ Less
Submitted 25 September, 2017; v1 submitted 19 July, 2017;
originally announced July 2017.
-
ICE-COLA: Towards fast and accurate synthetic galaxy catalogues optimizing a quasi $N$-body method
Authors:
Albert Izard,
Martin Crocce,
Pablo Fosalba
Abstract:
Next generation galaxy surveys demand the development of massive ensembles of galaxy mocks to model the observables and their covariances, what is computationally prohibitive using $N$-body simulations. COLA is a novel method designed to make this feasible by following an approximate dynamics but with up to 3 orders of magnitude speed-ups when compared to an exact $N$-body. In this paper we invest…
▽ More
Next generation galaxy surveys demand the development of massive ensembles of galaxy mocks to model the observables and their covariances, what is computationally prohibitive using $N$-body simulations. COLA is a novel method designed to make this feasible by following an approximate dynamics but with up to 3 orders of magnitude speed-ups when compared to an exact $N$-body. In this paper we investigate the optimization of the code parameters in the compromise between computational cost and recovered accuracy in observables such as two-point clustering and halo abundance. We benchmark those observables with a state-of-the-art $N$-body run, the MICE Grand Challenge simulation (MICE-GC). We find that using 40 time steps linearly spaced since $z_i \sim 20$, and a force mesh resolution three times finer than that of the number of particles, yields a matter power spectrum within $1\%$ for $k \lesssim 1\,h {\rm Mpc}^{-1}$ and a halo mass function within $5\%$ of those in the $N$-body. In turn the halo bias is accurate within $2\%$ for $k \lesssim 0.7\,h {\rm Mpc}^{-1}$ whereas, in redshift space, the halo monopole and quadrupole are within $4\%$ for $k \lesssim 0.4\,h {\rm Mpc}^{-1}$. These results hold for a broad range in redshift ($0 < z < 1$) and for all halo mass bins investigated ($M > 10^{12.5} \, h^{-1} \, {\rm M_{\odot}}$). To bring accuracy in clustering to one percent level we study various methods that re-calibrate halo masses and/or velocities. We thus propose an optimized choice of COLA code parameters as a powerful tool to optimally exploit future galaxy surveys.
△ Less
Submitted 6 April, 2016; v1 submitted 15 September, 2015;
originally announced September 2015.
-
nIFTy Cosmology: Galaxy/halo mock catalogue comparison project on clustering statistics
Authors:
Chia-Hsun Chuang,
Cheng Zhao,
Francisco Prada,
Emiliano Munari,
Santiago Avila,
Albert Izard,
Francisco-Shu Kitaura,
Marc Manera,
Pierluigi Monaco,
Steven Murray,
Alexander Knebe,
Claudia G. Scoccola,
Gustavo Yepes,
Juan Garcia-Bellido,
Felipe A. Marin,
Volker Muller,
Ramin Skibba,
Martin Crocce,
Pablo Fosalba,
Stefan Gottlober,
Anatoly A. Klypin,
Chris Power,
Charling Tao,
Victor Turchaninov
Abstract:
We present a comparison of major methodologies of fast generating mock halo or galaxy catalogues. The comparison is done for two-point and the three-point clustering statistics. The reference catalogues are drawn from the BigMultiDark N-body simulation. Both friend-of-friends (including distinct halos only) and spherical overdensity (including distinct halos and subhalos) catalogs have been used w…
▽ More
We present a comparison of major methodologies of fast generating mock halo or galaxy catalogues. The comparison is done for two-point and the three-point clustering statistics. The reference catalogues are drawn from the BigMultiDark N-body simulation. Both friend-of-friends (including distinct halos only) and spherical overdensity (including distinct halos and subhalos) catalogs have been used with the typical number density of a large-volume galaxy surveys. We demonstrate that a proper biasing model is essential for reproducing the power spectrum at quasilinear and even smaller scales. With respect to various clustering statistics a methodology based on perturbation theory and a realistic biasing model leads to very good agreement with N-body simulations. However, for the quadrupole of the correlation function or the power spectrum, only the method based on semi-N-body simulation could reach high accuracy (1% level) at small scales, i.e., r<25 Mpc/h or k>0.15 h/Mpc. Full N-body solutions will remain indispensable to produce reference catalogues. Nevertheless, we have demonstrated that the far more efficient approximate solvers can reach a few percent accuracy in terms of clustering statistics at the scales interesting for the large-scale structure analysis after calibration with a few reference N-body calculations. This makes them useful for massive production aimed at covariance studies, to scan large parameter spaces, and to estimate uncertainties in data analysis techniques, such as baryon acoustic oscillation reconstruction, redshift distortion measurements, etc.
△ Less
Submitted 16 June, 2015; v1 submitted 24 December, 2014;
originally announced December 2014.