ABSTRACT
We use galaxies from the illustristng, massiveblack-ii, and illustris-1 hydrodynamic simulations to investigate the behaviour of large scale galaxy intrinsic alignments. Our analysis spans ...four redshift slices over the approximate range of contemporary lensing surveys z = 0−1. We construct comparable weighted samples from the three simulations, which we then analyse using an alignment model that includes both linear and quadratic alignment contributions. Our data vector includes galaxy–galaxy, galaxy–shape, and shape–shape projected correlations, with the joint covariance matrix estimated analytically. In all of the simulations, we report non-zero IAs at the level of several σ. For a fixed lower mass threshold, we find a relatively strong redshift dependence in all three simulations, with the linear IA amplitude increasing by a factor of ∼2 between redshifts z = 0 and z = 1. We report no significant evidence for non-zero values of the tidal torquing amplitude, A2, in TNG, above statistical uncertainties, although MBII favours a moderately negative A2 ∼ −2. Examining the properties of the TATT model as a function of colour, luminosity and galaxy type (satellite or central), our findings are consistent with the most recent measurements on real data. We also outline a novel method for constraining the TATT model parameters directly from the pixelized tidal field, alongside a proof-of-concept exercise using TNG. This technique is shown to be promising, although comparison with previous results obtained via other methods is non-trivial.
ABSTRACT
In cosmology, we routinely choose between models to describe our data, and can incur biases due to insufficient models or lose constraining power with overly complex models. In this paper, ...we propose an empirical approach to model selection that explicitly balances parameter bias against model complexity. Our method uses synthetic data to calibrate the relation between bias and the χ2 difference between models. This allows us to interpret χ2 values obtained from real data (even if catalogues are blinded) and choose a model accordingly. We apply our method to the problem of intrinsic alignments – one of the most significant weak lensing systematics, and a major contributor to the error budget in modern lensing surveys. Specifically, we consider the example of the Dark Energy Survey Year 3 (DES Y3), and compare the commonly used non-linear alignment (NLA) and tidal alignment and tidal torque (TATT) models. The models are calibrated against bias in the Ωm–S8 plane. Once noise is accounted for, we find that it is possible to set a threshold Δχ2 that guarantees an analysis using NLA is unbiased at some specified level Nσ and confidence level. By contrast, we find that theoretically defined thresholds (based on, e.g. p-values for χ2) tend to be overly optimistic, and do not reliably rule out cosmological biases up to ∼1–2σ. Considering the real DES Y3 cosmic shear results, based on the reported difference in χ2 from NLA and TATT analyses, we find a roughly $30{{\ \rm per\ cent}}$ chance that were NLA to be the fiducial model, the results would be biased (in the Ωm–S8 plane) by more than 0.3σ. More broadly, the method we propose here is simple and general, and requires a relatively low level of resources. We foresee applications to future analyses as a model selection tool in many contexts.
ABSTRACT
Galaxy intrinsic alignments (IAs) have long been recognized as a significant contaminant to weak lensing-based cosmological inference. In this paper we seek to quantify the impact of a ...common modelling assumption in analytic descriptions of IAs: that of spherically symmetric dark matter haloes. Understanding such effects is important as the current generation of IA models are known to be limited, particularly on small scales, and building an accurate theoretical description will be essential for fully exploiting the information in future lensing data. Our analysis is based on a catalogue of 113 560 galaxies between z = 0.06 and 1.00 from massiveblack-ii, a hydrodynamical simulation of box length $100 \, h^{-1}$ Mpc. We find satellite anisotropy contributes at the level of $\ge 30\!-\!40{{\ \rm per\ cent}}$ to the small-scale alignment correlation functions. At separations larger than $1 \, h^{-1}$ Mpc the impact is roughly scale independent, inducing a shift in the amplitude of the IA power spectra of $\sim 20{{\ \rm per\ cent}}$. These conclusions are consistent across the redshift range and between the massiveblack-ii and the illustris simulations. The cosmological implications of these results are tested using a simulated likelihood analysis. Synthetic cosmic shear data are constructed with the expected characteristics (depth, area, and number density) of a future LSST-like survey. Our results suggest that modelling alignments using a halo model based upon spherical symmetry could potentially induce cosmological parameter biases at the ∼1.5σ level for S8 and w.
Abstract
We investigate the expected cosmological constraints from a combination of cosmic shear and large-scale galaxy clustering using realistic photometric redshift distributions. Introducing a ...systematic bias in the lensing distributions (of 0.05 in redshift) produces a >2σ bias in the recovered matter power spectrum amplitude and dark energy equation of state for preliminary Stage III surveys. We demonstrate that cosmological error can be largely removed by marginalizing over biases in the assumed weak-lensing redshift distributions. Furthermore, the cosmological constraining power is retained despite removing much of the information on the lensing redshift biases. This finding relies upon high-quality redshift estimates for the clustering sample, but does not require spectroscopy. All galaxies in this analysis can thus be assumed to come from a single photometric survey. We show that this internal constraint on redshift biases arises from complementary degeneracy directions between cosmic shear and the combination of galaxy clustering and shear–density cross-correlations. Finally we examine a case where the assumed redshift distributions differ from the truth by more than a simple uniform bias. We find that the effectiveness of this self-calibration method will depend on the survey details and the nature of the uncertainties on the estimated redshift distributions.
In this work, we use a suite of simulated images based on Year 1 of the Dark Energy Survey to explore the impact of galaxy neighbours on shape measurement and shear cosmology. The hoopoe image ...simulations include realistic blending, galaxy positions, and spatial variations in depth and point spread function properties. Using the im3shape maximum-likelihood shape measurement code, we identify four mechanisms by which neighbours can have a non-negligible influence on shear estimation. These effects, if ignored, would contribute a net multiplicative bias of m ~0.03–0.09 in the Year One of the Dark Energy Survey (DES Y1) im3shape catalogue, though the precise impact will be dependent on both the measurement code and the selection cuts applied. This can be reduced to percentage level or less by removing objects with close neighbours, at a cost to the effective number density of galaxies neff of 30 percent. We use the cosmological inference pipeline of DES Y1 to explore the cosmological implications of neighbour bias and show that omitting blending from the calibration simulation for DES Y1 would bias the inferred clustering amplitude S8 ≡ σ8(Ωm/0.3)0.5 by 2σ towards low values. Lastly, we use the hoopoe simulations to test the effect of neighbour-induced spatial correlations in the multiplicative bias. We find the impact on the recovered S8 of ignoring such correlations to be subdominant to statistical error at the current level of precision.
We perform a joint analysis of intrinsic alignments and cosmology using tomographic weak lensing, galaxy clustering and galaxy-galaxy lensing measurements from Year 1 (Y1) of the Dark Energy Survey. ...We define early- and late-type subsamples, which are found to pass a series of systematics tests, including for spurious photometric redshift error and point spread function correlations. We analyse these split data alongside the fiducial mixed Y1 sample using a range of intrinsic alignment models. In a fiducial Nonlinear Alignment Model (NLA) analysis, assuming a flat \lcdm~cosmology, we find a significant difference in intrinsic alignment amplitude, with early-type galaxies favouring $A_\mathrm{IA} = 2.38^{+0.32}_{-0.31}$ and late-type galaxies consistent with no intrinsic alignments at $0.05^{+0.10}_{-0.09}$. We find weak evidence of a diminishing alignment amplitude at higher redshifts in the early-type sample. The analysis is repeated using a number of extended model spaces, including a physically motivated model that includes both tidal torquing and tidal alignment mechanisms. In multiprobe likelihood chains in which cosmology, intrinsic alignments in both galaxy samples and all other relevant systematics are varied simultaneously, we find the tidal alignment and tidal torquing parts of the intrinsic alignment signal have amplitudes $A_1 = 2.66 ^{+0.67}_{-0.66}$, $A_2=-2.94^{+1.94}_{-1.83}$, respectively, for early-type galaxies and $A_1 = 0.62 ^{+0.41}_{-0.41}$, $A_2 = -2.26^{+1.30}_{-1.16}$ for late-type galaxies. In the full (mixed) Y1 sample the best constraints are $A_1 = 0.70 ^{+0.41}_{-0.38}$, $A_2 = -1.36 ^{+1.08}_{-1.41}$. For all galaxy splits and IA models considered, we report cosmological parameter constraints that are consistent with the results of Troxel et al. (2017) and Dark Energy Survey Collaboration (2017).
We present constraints on extensions of the minimal cosmological models dominated by dark matter and dark energy, ΛCDM and wCDM, by using a combined analysis of galaxy clustering and weak ...gravitational lensing from the first-year data of the Dark Energy Survey (DES Y1) in combination with external data. We consider four extensions of the minimal dark energy-dominated scenarios: (1) nonzero curvature Ωk, (2) number of relativistic species Neff different from the standard value of 3.046, (3) time-varying equation-of-state of dark energy described by the parameters w0 and wa (alternatively quoted by the values at the pivot redshift, wp, and wa), and (4) modified gravity described by the parameters μ0 and Σ0 that modify the metric potentials. We also consider external information from Planck cosmic microwave background measurements; baryon acoustic oscillation measurements from SDSS, 6dF, and BOSS; redshift-space distortion measurements from BOSS; and type Ia supernova information from the Pantheon compilation of datasets. Constraints on curvature and the number of relativistic species are dominated by the external data; when these are combined with DES Y1, we find Ωk=0.0020−0.0032+0.0037 at the 68% confidence level, and the upper limit Neff<3.28(3.55) at 68% (95%) confidence, assuming a hard prior Neff>3.0. For the time-varying equation-of-state, we find the pivot value (wp,wa)=(−0.91−0.23+0.19,−0.57−1.11+0.93) at pivot redshift zp=0.27 from DES alone, and (wp,wa)=(−1.01−0.04+0.04,−0.28−0.48+0.37) at zp=0.20 from DES Y1 combined with external data; in either case we find no evidence for the temporal variation of the equation of state. For modified gravity, we find the present-day value of the relevant parameters to be Σ0=0.43−0.29+0.28 from DES Y1 alone, and (Σ0,μ0)=(0.06−0.07+0.08,−0.11−0.46+0.42) from DES Y1 combined with external data. These modified-gravity constraints are consistent with predictions from general relativity.
ABSTRACT
As the statistical power of galaxy weak lensing reaches per cent level precision, large, realistic, and robust simulations are required to calibrate observational systematics, especially ...given the increased importance of object blending as survey depths increase. To capture the coupled effects of blending in both shear and photometric redshift calibration, we define the effective redshift distribution for lensing, nγ(z), and describe how to estimate it using image simulations. We use an extensive suite of tailored image simulations to characterize the performance of the shear estimation pipeline applied to the Dark Energy Survey (DES) Year 3 data set. We describe the multiband, multi-epoch simulations, and demonstrate their high level of realism through comparisons to the real DES data. We isolate the effects that generate shear calibration biases by running variations on our fiducial simulation, and find that blending-related effects are the dominant contribution to the mean multiplicative bias of approximately $-2{{\ \rm per\ cent}}$. By generating simulations with input shear signals that vary with redshift, we calibrate biases in our estimation of the effective redshift distribution, and demonstrate the importance of this approach when blending is present. We provide corrected effective redshift distributions that incorporate statistical and systematic uncertainties, ready for use in DES Year 3 weak lensing analyses.
We present the first joint analysis of cluster abundances and auto or cross-correlations of three cosmic tracer fields: galaxy density, weak gravitational lensing shear, and cluster density split by ...optical richness. From a joint analysis (4×2pt+N) of cluster abundances, three cluster cross-correlations, and the auto correlations of the galaxy density measured from the first year data of the Dark Energy Survey, we obtain Ω_{m}=0.305_{-0.038}^{+0.055} and σ_{8}=0.783_{-0.054}^{+0.064}. This result is consistent with constraints from the DES-Y1 galaxy clustering and weak lensing two-point correlation functions for the flat νΛCDM model. Consequently, we combine cluster abundances and all two-point correlations from across all three cosmic tracer fields (6×2pt+N) and find improved constraints on cosmological parameters as well as on the cluster observable-mass scaling relation. This analysis is an important advance in both optical cluster cosmology and multiprobe analyses of upcoming wide imaging surveys.
We perform a joint analysis of the counts of redMaPPer clusters selected from the Dark Energy Survey (DES) year 1 data and multiwavelength follow-up data collected within the 2500 deg2 South Pole ...Telescope (SPT) Sunyaev-Zel'dovich (SZ) survey. The SPT follow-up data, calibrating the richness-mass relation of the optically selected redMaPPer catalog, enable the cosmological exploitation of the DES cluster abundance data. To explore possible systematics related to the modeling of projection effects, we consider two calibrations of the observational scatter on richness estimates: a simple Gaussian model which account only for the background contamination (BKG), and a model which further includes contamination and incompleteness due to projection effects (PRJ). Assuming either a Λ CDM + ∑ mν or w CDM + ∑ mν cosmology, and for both scatter models, we derive cosmological constraints consistent with multiple cosmological probes of the low and high redshift Universe, and in particular with the SPT cluster abundance data. This result demonstrates that the DES Y1 and SPT cluster counts provide consistent cosmological constraints, if the same mass calibration data set is adopted. It thus supports the conclusion of the DES Y1 cluster cosmology analysis which interprets the tension observed with other cosmological probes in terms of systematics affecting the stacked weak lensing analysis of optically selected low–richness clusters. Finally, we analyze the first combined optically SZ selected cluster catalog obtained by including the SPT sample above the maximum redshift probed by the DES Y1 redMaPPer sample (z = 0.65). Besides providing a mild improvement of the cosmological constraints, this data combination serves as a stricter test of our scatter models: the PRJ model, providing scaling relations consistent between the two abundance and multiwavelength follow-up data, is favored over the BKG model.