Efficient estimators of Fourier-space statistics for large number of objects rely on fast Fourier transforms (FFTs), which are affected by aliasing from unresolved small-scale modes due to the finite ...FFT grid. Aliasing takes the form of a sum over images, each of them corresponding to the Fourier content displaced by increasing multiples of the sampling frequency of the grid. These spurious contributions limit the accuracy in the estimation of Fourier-space statistics, and are typically ameliorated by simultaneously increasing grid size and discarding high-frequency modes. This results in inefficient estimates for e.g. the power spectrum when desired systematic biases are well under per cent level. We show that using interlaced grids removes odd images, which include the dominant contribution to aliasing. In addition, we discuss the choice of interpolation kernel used to define density perturbations on the FFT grid and demonstrate that using higher order interpolation kernels than the standard Cloud-In-Cell algorithm results in significant reduction of the remaining images. We show that combining fourth-order interpolation with interlacing gives very accurate Fourier amplitudes and phases of density perturbations. This results in power spectrum and bispectrum estimates that have systematic biases below 0.01 per cent all the way to the Nyquist frequency of the grid, thus maximizing the use of unbiased Fourier coefficients for a given grid size and greatly reducing systematics for applications to large cosmological data sets.
We present a new N-body simulation from the Marenostrum Institut de Ciències de l'Espai (MICE) collaboration, the MICE Grand Challenge (MICE-GC), containing about 70 billion dark matter particles in ...a (3 Gpc h
−1)3 comoving volume. Given its large volume and fine spatial resolution, spanning over five orders of magnitude in dynamic range, it allows an accurate modelling of the growth of structure in the universe from the linear through the highly non-linear regime of gravitational clustering. We validate the dark matter simulation outputs using 3D and 2D clustering statistics, and discuss mass-resolution effects in the non-linear regime by comparing to previous simulations and the latest numerical fits. We show that the MICE-GC run allows for a measurement of the BAO feature with per cent level accuracy and compare it to state-of-the-art theoretical models. We also use sub-arcmin resolution pixelized 2D maps of the dark matter counts in the lightcone to make tomographic analyses in real and redshift space. Our analysis shows the simulation reproduces the Kaiser effect on large scales, whereas we find a significant suppression of power on non-linear scales relative to the real space clustering. We complete our validation by presenting an analysis of the three-point correlation function in this and previous MICE simulations, finding further evidence for mass-resolution effects. This is the first of a series of three papers in which we present the MICE-GC simulation, along with a wide and deep mock galaxy catalogue built from it. This mock is made publicly available through a dedicated web portal, http://cosmohub.pic.es.
This is the second in a series of three papers in which we present an end-to-end simulation from the MICE collaboration, the MICE Grand Challenge (MICE-GC) run. The N-body contains about 70 billion ...dark-matter particles in a (3 h
−1 Gpc)3 comoving volume spanning five orders of magnitude in dynamical range. Here, we introduce the halo and galaxy catalogues built upon it, both in a wide (5000 deg2) and deep (z < 1.4) lightcone and in several comoving snapshots. Haloes were resolved down to few 1011 h
−1 M⊙. This allowed us to model galaxies down to absolute magnitude M
r
< −18.9. We used a new hybrid halo occupation distribution and abundance matching technique for galaxy assignment. The catalogue includes the spectral energy distributions of all galaxies. We describe a variety of halo and galaxy clustering applications. We discuss how mass resolution effects can bias the large-scale two-pt clustering amplitude of poorly resolved haloes at the ≲5 per cent level, and their three-pt correlation function. We find a characteristic scale-dependent bias of ≲6 per cent across the BAO feature for haloes well above M
⋆ ∼ 1012 h
−1 M⊙ and for luminous red galaxy like galaxies. For haloes well below M
⋆ the scale dependence at 100 h
−1 Mpc is ≲2 per cent. Lastly, we discuss the validity of the large-scale Kaiser limit across redshift and departures from it towards non-linear scales. We make the current version of the lightcone halo and galaxy catalogue (MICECATv1.0) publicly available through a dedicated web portal to help develop and exploit the new generation of astronomical surveys.
In Paper I of this series, we presented a new N-body light-cone simulation from the MICE Collaboration, the MICE Grand Challenge (MICE-GC), containing about 70 billion dark-matter particles in a (3 h
...−1 Gpc)3 comoving volume, from which we built halo and galaxy catalogues using a Halo Occupation Distribution and Halo Abundance Matching technique, as presented in the companion Paper II. Given its large volume and fine mass resolution, the MICE-GC simulation also allows an accurate modelling of the lensing observables from upcoming wide and deep galaxy surveys. In the last paper of this series (Paper III), we describe the construction of all-sky lensing maps, following the ‘Onion Universe’ approach, and discuss their properties in the light-cone up to z = 1.4 with sub-arcminute spatial resolution. By comparing the convergence power spectrum in the MICE-GC to lower mass-resolution (i.e. particle mass ∼1011 h
−1 M⊙) simulations, we find that resolution effects are at the 5 per cent level for multipoles ℓ ∼ 103 and 20 per cent for ℓ ∼ 104. Resolution effects have a much lower impact on our simulation, as shown by comparing the MICE-GC to recent numerical fits by Takahashi. We use the all-sky lensing maps to model galaxy lensing properties, such as the convergence, shear, and lensed magnitudes and positions, and validate them thoroughly using galaxy shear auto and cross-correlations in harmonic and configuration space. Our results show that the galaxy lensing mocks here presented can be used to accurately model lensing observables down to arcminute scales. Accompanying this series of papers, we make a first public data release of the MICE-GC galaxy mock, the MICECAT v1.0, through a dedicated web-portal for the MICE simulations, http://cosmohub.pic.es, to help developing and exploiting the new generation of astronomical surveys.
We study the impact of setting initial conditions in numerical simulations using the standard procedure based on the Zel'dovich approximation (ZA). As it is well known from the perturbation theory, ...ZA initial conditions have incorrect second- and higher-order growth and therefore excite long-lived transients in the evolution of the statistical properties of density and velocity fields. We also study the improvement brought by using more accurate initial conditions based on second-order Lagrangian perturbation theory (2LPT). We show that 2LPT initial conditions reduce transients significantly and thus are much more appropriate for numerical simulations devoted to precision cosmology. Using controlled numerical experiments with ZA and 2LPT initial conditions, we show that simulations started at redshift zi= 49 using the ZA underestimate the power spectrum in the non-linear regime by about 2, 4 and 8 per cent at z= 0, 1, and 3, respectively, whereas the mass function of dark matter haloes is underestimated by 5 per cent at m= 1015 M⊙ h−1 (z= 0) and 10 per cent at m= 2 × 1014 M⊙ h−1 (z= 1). The clustering of haloes is also affected to the few per cent level at z= 0. These systematics effects are typically larger than statistical uncertainties in recent mass function and power spectrum fitting formulae extracted from numerical simulations. At large scales, the measured transients in higher-order correlations can be understood from first principle calculations based on perturbation theory.
We present a method to build mock galaxy catalogues starting from a halo catalogue that uses halo occupation distribution (HOD) recipes as well as the subhalo abundance matching (SHAM) technique. ...Combining both prescriptions we are able to push the absolute magnitude of the resulting catalogue to fainter luminosities than using just the SHAM technique and can interpret our results in terms of the HOD modelling. We optimize the method by populating with galaxies friends-of-friends dark matter haloes extracted from the Marenostrum Institut de Ciències de l'Espai dark matter simulations and comparing them to observational constraints. Our resulting mock galaxy catalogues manage to reproduce the observed local galaxy luminosity function and the colour–magnitude distribution as observed by the Sloan Digital Sky Survey. They also reproduce the observed galaxy clustering properties as a function of luminosity and colour. In order to achieve that, the algorithm also includes scatter in the halo mass–galaxy luminosity relation derived from direct SHAM and a modified Navarro–Frenk–White mass density profile to place satellite galaxies in their host dark matter haloes. Improving on general usage of the HOD that fits the clustering for given magnitude limited samples, our catalogues are constructed to fit observations at all luminosities considered and therefore for any luminosity subsample. Overall, our algorithm is an economic procedure of obtaining galaxy mock catalogues down to faint magnitudes that are necessary to understand and interpret galaxy surveys.
The Dark Energy Survey Data Release 2 Abbott, T. M. C.; Amon, A.; Avila, S. ...
The Astrophysical journal. Supplement series,
08/2021, Volume:
255, Issue:
2
Journal Article
Peer reviewed
Open access
Abstract
We present the second public data release of the Dark Energy Survey, DES DR2, based on optical/near-infrared imaging by the Dark Energy Camera mounted on the 4 m Blanco telescope at Cerro ...Tololo Inter-American Observatory in Chile. DES DR2 consists of reduced single-epoch and coadded images, a source catalog derived from coadded images, and associated data products assembled from 6 yr of DES science operations. This release includes data from the DES wide-area survey covering ∼5000 deg
2
of the southern Galactic cap in five broad photometric bands,
grizY
. DES DR2 has a median delivered point-spread function FWHM of
g
= 1.11″,
r
= 0.95″,
i
= 0.88″,
z
= 0.83″, and
Y
= 0.″90, photometric uniformity with a standard deviation of < 3 mmag with respect to Gaia DR2
G
band, a photometric accuracy of ∼11 mmag, and a median internal astrometric precision of ∼27 mas. The median coadded catalog depth for a 1.″95 diameter aperture at signal-to-noise ratio = 10 is
g
= 24.7,
r
= 24.4,
i
= 23.8,
z
= 23.1, and
Y
= 21.7 mag. DES DR2 includes ∼691 million distinct astronomical objects detected in 10,169 coadded image tiles of size 0.534 deg
2
produced from 76,217 single-epoch images. After a basic quality selection, benchmark galaxy and stellar samples contain 543 million and 145 million objects, respectively. These data are accessible through several interfaces, including interactive image visualization tools, web-based query clients, image cutout servers, and Jupyter notebooks. DES DR2 constitutes the largest photometric data set to date at the achieved depth and photometric precision.
We present the methodology for a joint cosmological analysis of weak gravitational lensing from the fourth data release of the ESO Kilo-Degree Survey (KiDS-1000) and galaxy clustering from the ...partially overlapping Baryon Oscillation Spectroscopic Survey (BOSS) and the 2-degree Field Lensing Survey (2dFLenS). Cross-correlations between BOSS and 2dFLenS galaxy positions and source galaxy ellipticities have been incorporated into the analysis, necessitating the development of a hybrid model of non-linear scales that blends perturbative and non-perturbative approaches, and an assessment of signal contributions by astrophysical effects. All weak lensing signals were measured consistently via Fourier-space statistics that are insensitive to the survey mask and display low levels of mode mixing. The calibration of photometric redshift distributions and multiplicative gravitational shear bias has been updated, and a more complete tally of residual calibration uncertainties was propagated into the likelihood. A dedicated suite of more than 20 000 mocks was used to assess the performance of covariance models and to quantify the impact of survey geometry and spatial variations of survey depth on signals and their errors. The sampling distributions for the likelihood and the
χ
2
goodness-of-fit statistic have been validated, with proposed changes for calculating the effective number of degrees of freedom. The prior volume was explicitly mapped, and a more conservative, wide top-hat prior on the key structure growth parameter
S
8
=
σ
8
(Ω
m
/0.3)
1/2
was introduced. The prevalent custom of reporting
S
8
weak lensing constraints via point estimates derived from its marginal posterior is highlighted to be easily misinterpreted as yielding systematically low values of
S
8
, and an alternative estimator and associated credible interval are proposed. Known systematic effects pertaining to weak lensing modelling and inference are shown to bias
S
8
by no more than 0.1 standard deviations, with the caveat that no conclusive validation data exist for models of intrinsic galaxy alignments. Compared to the previous KiDS analyses,
S
8
constraints are expected to improve by 20% for weak lensing alone and by 29% for the joint analysis.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
ABSTRACT We search for excess γ-ray emission coincident with the positions of confirmed and candidate Milky Way satellite galaxies using six years of data from the Fermi Large Area Telescope (LAT). ...Our sample of 45 stellar systems includes 28 kinematically confirmed dark-matter-dominated dwarf spheroidal galaxies (dSphs) and 17 recently discovered systems that have photometric characteristics consistent with the population of known dSphs. For each of these targets, the relative predicted γ-ray flux due to dark matter annihilation is taken from kinematic analysis if available, and estimated from a distance-based scaling relation otherwise, assuming that the stellar systems are DM-dominated dSphs. LAT data coincident with four of the newly discovered targets show a slight preference (each 2 local) for γ-ray emission in excess of the background. However, the ensemble of derived γ-ray flux upper limits for individual targets is consistent with the expectation from analyzing random blank-sky regions, and a combined analysis of the population of stellar systems yields no globally significant excess (global significance ). Our analysis has increased sensitivity compared to the analysis of 15 confirmed dSphs by Ackermann et al. The observed constraints on the DM annihilation cross section are statistically consistent with the background expectation, improving by a factor of ∼2 for large DM masses ( and ) and weakening by a factor of ∼1.5 at lower masses relative to previously observed limits.
This overview paper describes the legacy prospect and discovery potential of the Dark Energy Survey (DES) beyond cosmological studies, illustrating it with examples from the DES early data. DES is ...using a wide-field camera (DECam) on the 4 m Blanco Telescope in Chile to image 5000 sq deg of the sky in five filters (grizY). By its completion, the survey is expected to have generated a catalogue of 300 million galaxies with photometric redshifts and 100 million stars. In addition, a time-domain survey search over 27 sq deg is expected to yield a sample of thousands of Type Ia supernovae and other transients. The main goals of DES are to characterize dark energy and dark matter, and to test alternative models of gravity; these goals will be pursued by studying large-scale structure, cluster counts, weak gravitational lensing and Type Ia supernovae. However, DES also provides a rich data set which allows us to study many other aspects of astrophysics. In this paper, we focus on additional science with DES, emphasizing areas where the survey makes a difference with respect to other current surveys. The paper illustrates, using early data (from ‘Science Verification’, and from the first, second and third seasons of observations), what DES can tell us about the Solar system, the Milky Way, galaxy evolution, quasars and other topics. In addition, we show that if the cosmological model is assumed to be Λ+cold dark matter, then important astrophysics can be deduced from the primary DES probes. Highlights from DES early data include the discovery of 34 trans-Neptunian objects, 17 dwarf satellites of the Milky Way, one published z > 6 quasar (and more confirmed) and two published superluminous supernovae (and more confirmed).
PDF
