ABSTRACT
We present a fast methodology to produce mock observations of the thermal and kinetic Sunyaev–Zel’dovich (SZ) effects based on the dark matter only N-body simulations coupled with the ...analytic intracluster medium model. The methods employ two different approaches: halo-based pasting (HP) and particle-based pasting (PP). The former pastes gas density and pressure on to haloes and requires only a halo catalogue, and the latter considers the contribution from field particles as well, i.e. particles that do not belong to any haloes and thus utilize the full particle information. Therefore, the PP algorithm incorporates secondary effects beyond the HP algorithm: asphericity of haloes and contribution from diffuse gas. In particular, such a diffuse component is the dominant source of the kinetic SZ effect. As validation of our methods, we have produced 108 all-sky maps with HP and 108 flat-sky maps, which cover $5 \times 5 \, \mathrm{deg}^2$ with both HP and PP, and measured power spectra of the maps. Our method can produce a mock map within a few hours, even for all-sky coverage with a parallel computational environment. The power spectra of HP maps are consistent with the halo model prediction of the thermal SZ effect. On the other hand, the power spectra of PP maps are suppressed due to the halo asphericity but can reproduce better the theoretical prediction for the kinetic SZ effect. We discuss the utility of baryon-pasted mock SZ maps for estimating the covariance matrix of SZ statistics and modelling the selection and projection effects for cluster cosmology.
An accurate theoretical template of the redshift-space galaxy power spectrum, if applicable out to nonlinear scales, enables us to extract more stringent and robust constraints on cosmological ...parameters from the measured galaxy clustering. In this work, we develop a simulation-based template, so-called emulator, for the redshift-space power spectrum of dark matter halos. Using the redshift-space halo power spectra measured from the dark quest N -body simulation suite that covers 101 flat-geometry w -cold dark matter ( w CDM ) cosmologies around the Planck Λ CDM model, we feed these data into a feed-forward neural network to build the fast and accurate emulation of the power spectrum from the linear to nonlinear scales up to k ≃ 0.6 h Mpc−1. Our emulator achieves about 1% and 5% fractional accuracies in predicting the monopole and quadrupole moments of the power spectrum, respectively, for halos of ∼ 1013 h−1 M⊙ that correspond to host halos of the Sloan Digital Sky Survey (SDSS) LOWZ- and CMASS (constant mass)-like galaxies, where the achieved accuracies are sufficient compared to the statistical errors of SDSS volume. The validation and performance of the emulator are given by the comparison of the emulator predictions with the power spectra directly measured from the simulations for validation sets that are not used in the training. We demonstrate that the emulator outputs can be used to make model predictions for the redshift-space power spectrum of galaxies by employing user-fed models for the halo-galaxy connection, such as the halo occupation distribution. The emulator allows us to easily incorporate the Finger-of-God effect due to the virial motions of galaxies and the Alcock-Paczyński distortions. Our code can compute the redshift-space galaxy power spectrum in a CPU subseconds and is ready to perform the emulator-based cosmological analysis for the exiting and upcoming galaxy redshift surveys.
We present 108 full-sky gravitational lensing simulation data sets generated by performing multiple-lens plane ray-tracing through high-resolution cosmological N-body simulations. The data sets ...include full-sky convergence and shear maps from redshifts z = 0.05 to 5.3 at intervals of 150 h − 1 Mpc comoving radial distance (corresponding to a redshift interval of Δ z 0.05 at the nearby universe), enabling the construction of a mock shear catalog for an arbitrary source distribution up to z = 5.3. The dark matter halos are identified from the same N-body simulations with enough mass resolution to resolve the host halos of the Sloan Digital Sky Survey (SDSS) CMASS and luminous red galaxies (LRGs). Angular positions and redshifts of the halos are provided by a ray-tracing calculation, enabling the creation of a mock halo catalog to be used for galaxy-galaxy and cluster-galaxy lensing. The simulation also yields maps of gravitational lensing deflections for a source redshift at the last scattering surface, and we provide 108 realizations of lensed cosmic microwave background (CMB) maps in which the post-Born corrections caused by multiple light scattering are included. We present basic statistics of the simulation data, including the angular power spectra of cosmic shear, CMB temperature and polarization anisotropies, galaxy-galaxy lensing signals for halos, and their covariances. The angular power spectra of the cosmic shear and CMB anisotropies agree with theoretical predictions within 5% up to = 3000 (or at an angular scale θ > 0.5 arcmin). The simulation data sets are generated primarily for the ongoing Subaru Hyper Suprime-Cam survey, but are freely available for download at http://cosmo.phys.hirosaki-u.ac.jp/takahasi/allsky_raytracing/.
Abstract
Emission line galaxies (ELGs), more generally star-forming galaxies, are valuable tracers of large-scale structure and therefore main targets of upcoming wide-area spectroscopic galaxy ...surveys. We propose a fixed-aperture shape estimator of each ELG for extracting the intrinsic alignment (IA) signal, and assess the performance of the method using image simulations of ELGs generated from the IllustrisTNG simulation including observational effects such as the sky background noise. We show that our method enables a significant detection of the IA power spectrum with the linear-scale coefficient
A
IA
≃ (13–15) ± 3.0 up to
z
= 2, even from the small simulation volume
∼
0.009
(
h
−
1
Gpc
)
3
, in contrast to the null detection with the standard method. Thus the ELG IA signal, measured with our method, opens up opportunities to exploit cosmology and galaxy physics in high-redshift universe.
We study the accuracy with which cosmological parameters can be determined from a real space power spectrum of matter density contrast at weakly nonlinear scales using analytical approaches. From ...power spectra measured in N -body simulations and using the Markov chain Monte Carlo technique, the best-fitting cosmological input parameters are determined with several analytical methods as a theoretical template, such as the standard perturbation theory, the regularized perturbation theory, and the effective field theory. We show that at redshift 1, all two-loop level calculations can fit the measured power spectrum down to scales k ∼ 0.2 h Mpc − 1 , and cosmological parameters are successfully estimated in an unbiased way. Introducing the figure of bias (FoB) and figure of merit (FoM) parameter, we determine the validity range of those models and then evaluate their relative performances. With one free parameter, namely the damping scale, the regularized perturbation theory is found to be able to provide the largest FoM parameter while keeping the FoB in the acceptance range.
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.
We investigate the splashback features of dark-matter halos based on cosmic density and velocity fields. Besides the density correlation function binned by the halo orientation angle, which was used ...in the literature, we introduce, for the first time, the corresponding velocity statistic, alignment momentum correlation function, to take into account the asphericity of halos. Using large-volume, high-resolution N-body simulations, we measure the alignment statistics of density and velocity. On halo scales, x∼R200m∼1 h−1 Mpc, we detect a sharp steepening in the momentum correlation associated with the physical halo boundary, or the splashback feature, which is found more prominent than in the density correlation. We also find that the splashback radius determined from the density correlation becomes ∼3.5% smaller than that from the momentum correlation, with their correlation coefficient being 0.605. Moreover, the orientation-dependent splashback feature due to halo asphericity is measured when the density profile is determined by dark-matter particles, which can be used as a test of collisional cold dark matter since the halo shape is predicted to be rounder in such a model.
ABSTRACT
Cross-correlation analysis of the thermal Sunyaev–Zel’dovich (tSZ) effect and weak gravitational lensing (WL) provides a powerful probe of cosmology and astrophysics of the intracluster ...medium. We present the measurement of the cross-correlation of tSZ and WL from Planck and Subaru Hyper-Suprime Cam. The combination enables us to study cluster astrophysics at high redshift. We use the tSZ-WL cross-correlation and the tSZ autopower spectrum measurements to place a tight constraint on the hydrostatic mass bias, which is a measure of the degree of non-thermal pressure support in galaxy clusters. With the prior on cosmological parameters derived from the analysis of the cosmic microwave background anisotropies by Planck and taking into account foreground contributions both in the tSZ autopower spectrum and the tSZ-WL cross-correlation, the hydrostatic mass bias is estimated to be $26.9^{+8.9}_{-4.4} {{\ \rm per\ cent}}$ ($68{{\ \rm per\ cent}}$ CL), which is consistent with recent measurements by mass calibration techniques.
We constrain the scaling relation between optical richness (λ) and halo mass (M) for a sample of Sloan Digital Sky Survey (SDSS) red-sequence Matched-filter Probabilistic Percolation (redMaPPer) ...galaxy clusters within the context of the Planck cosmological model. We use a forward modeling approach where we model the probability distribution of optical richness for a given mass, . To model the abundance and the stacked lensing profiles, we use an emulator specifically built to interpolate the halo mass function and the stacked lensing profile for an arbitrary set of halo mass and redshift, which is calibrated based on a suite of high-resolution N-body simulations. We apply our method to 8312 SDSS redMaPPer clusters with 20 ≤ λ ≤ 100 and 0.10 ≤ zλ ≤ 0.33 and show that the lognormal distribution model for , with four free parameters, well reproduces the measured abundances and lensing profiles simultaneously. The constraints are characterized by the mean relation, , with and (68% CL), where the pivot mass scale Mpivot = 3 × 1014 h−1 M , and the scatter with and . We find that a large scatter in halo masses is required at the lowest-richness bins (20 ≤ λ 30) in order to reproduce the measurements. Without such a large scatter, the model prediction for the lensing profiles tends to overestimate the measured amplitudes. This might imply a possible contamination of intrinsically low-richness clusters due to the projection effects. Such a low-mass halo contribution is significantly reduced when applying our method to the sample of 30 ≤ λ ≤ 100.
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 cosmological constraints from peak and ...minimum counts. Using the 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{{\ \rm per\ cent}}$ credible intervals from a combination of lensing minima and peaks that are significantly stronger than from the power spectrum alone, by $44{{\ \rm per\ cent}}$, $11{{\ \rm per\ cent}}$, and $63{{\ \rm per\ cent}}$ for the neutrino mass sum ∑mν, matter density Ωm, and amplitude of fluctuation As, respectively. We explore the effect of baryonic processes on lensing minima and peaks using the hydrodynamical simulations BAHAMAS and Osato15. We find that ignoring baryonic effects would lead to strong (≈4σ) biases in inferences from peak counts, but negligible (≈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.