No Arabic abstract
Developing analysis pipelines based on statistics beyond two-point functions is critical for extracting a maximal amount of cosmological information from current and upcoming weak lensing surveys. In this paper, we study the impact of the intrinsic alignment of galaxies (IA) on three promising probes measured from aperture mass maps -- the lensing peaks, minima and full PDF, in comparison and in combination with the shear two-point correlation functions ($gamma$-2PCFs). Our two-dimensional IA infusion method converts the light-cone-projected mass sheets into projected tidal tensors, which are then linearly coupled to an intrinsic ellipticity component with a strength controlled by the coupling parameter $A_{rm IA}$. We validate our method with the $gamma$-2PCFs statistics, recovering well the analytical calculations from the linear alignment model of citet{BridleKing} in a full tomographic setting, and for different $A_{rm IA}$ values. We next use our method to infuse at the galaxy catalogue level a non-linear IA model that includes the density-weighting term introduced in citet{Blazek2015}, and compute the impact on the three aperture mass map statistics. We find that large snr peaks are maximally affected, with deviations reaching 30% (10%) for a {it Euclid}-like (KiDS-like) survey. Modelling the signal in a $w$CDM cosmology universe with $N$-body simulations, we forecast the cosmological bias caused by unmodelled IA for 100 deg$^2$ of {it Euclid}-like data, finding very large offsets in $w_0$ (5-10$sigma_{rm stat}$), $Omega_{rm m}$ (4-6$sigma_{rm stat}$), and $S_8 equiv sigma_8sqrt{Omega_{rm m}/0.3}$ ($sim$3$sigma_{rm stat}$). The method presented in this paper offers a compelling avenue to account for IA in beyond-two-point weak lensing statistics, with a flexibility comparable to that of current $gamma$-2PCFs IA analytical models.
We consider the effect of galaxy intrinsic alignments (IAs) on dark energy constraints from weak gravitational lensing. We summarise the latest version of the linear alignment model of IAs, following the brief note of Hirata & Seljak (2010) and further interpretation in Laszlo et al. (2011). We show the cosmological bias on the dark energy equation of state parameters w0 and wa that would occur if IAs were ignored. We find that w0 and wa are both catastrophically biased, by an absolute value of just greater than unity under the Fisher matrix approximation. This contrasts with a bias several times larger for the earlier IA implementation. Therefore there is no doubt that IAs must be taken into account for future Stage III experiments and beyond. We use a flexible grid of IA and galaxy bias parameters as used in previous work, and investigate what would happen if the universe used the latest IA model, but we assumed the earlier version. We find that despite the large difference between the two IA models, the grid flexibility is sufficient to remove cosmological bias and recover the correct dark energy equation of state. In an appendix, we compare observed shear power spectra to those from a popular previous implementation and explain the differences.
We constrain cosmological parameters from a joint cosmic shear analysis of peak-counts and the two-point shear correlation functions, as measured from the Dark Energy Survey (DES-Y1). We find the structure growth parameter $S_8equiv sigma_8sqrt{Omega_{rm m}/0.3} = 0.766^{+0.033}_{-0.038}$, which at 4.8% precision, provides one of the tightest constraints on $S_8$ from the DES-Y1 weak lensing data. In our simulation-based method we determine the expected DES-Y1 peak-count signal for a range of cosmologies sampled in four $w$CDM parameters ($Omega_{rm m}$, $sigma_8$, $h$, $w_0$). We also determine the joint covariance matrix with over 1000 realisations at our fiducial cosmology. With mock DES-Y1 data we calibrate the impact of photometric redshift and shear calibration uncertainty on the peak-count, marginalising over these uncertainties in our cosmological analysis. Using dedicated training samples we show that our measurements are unaffected by mass resolution limits in the simulation, and that our constraints are robust against uncertainty in the effect of baryon feedback. Accurate modelling for the impact of intrinsic alignments on the tomographic peak-count remains a challenge, currently limiting our exploitation of cross-correlated peak counts between high and low redshift bins. We demonstrate that once calibrated, a fully tomographic joint peak-count and correlation functions analysis has the potential to reach a 3% precision on $S_8$ for DES-Y1. Our methodology can be adopted to model any statistic that is sensitive to the non-Gaussian information encoded in the shear field. In order to accelerate the development of these beyond-two-point cosmic shear studies, our simulations are made available to the community upon request.
One of the most pernicious theoretical systematics facing upcoming gravitational lensing surveys is the uncertainty introduced by the effects of baryons on the power spectrum of the convergence field. One method that has been proposed to account for these effects is to allow several additional parameters (that characterize dark matter halos) to vary and to fit lensing data to these halo parameters concurrently with the standard set of cosmological parameters. We test this method. In particular, we use this technique to model convergence power spectrum predictions from a set of cosmological simulations. We estimate biases in dark energy equation of state parameters that would be incurred if one were to fit the spectra predicted by the simulations either with no model for baryons, or with the proposed method. We show that neglecting baryonic effect leads to biases in dark energy parameters that are several times the statistical errors for a survey like the Dark Energy Survey. The proposed method to correct for baryonic effects renders the residual biases in dark energy equation of state parameters smaller than the statistical errors. These results suggest that this mitigation method may be applied to analyze convergence spectra from a survey like the Dark Energy Survey. For significantly larger surveys, such as will be carried out by the Large Synoptic Survey Telescope, the biases introduced by baryonic effects are much more significant. We show that this mitigation technique significantly reduces the biases for such larger surveys, but that a more effective mitigation strategy will need to be developed in order ensure that the residual biases in these surveys fall below the statistical errors.
Correlations of galaxy ellipticities with large-scale structure, due to galactic tidal interactions, provide a potentially significant contaminant to measurements of cosmic shear. However, these intrinsic alignments are still poorly understood for galaxies at the redshifts typically used in cosmic shear analyses. For spiral galaxies, it is thought that tidal torquing is significant in determining alignments resulting in zero correlation between the intrinsic ellipticity and the gravitational potential in linear theory. Here, we calculate the leading-order correction to this result in the tidal-torque model from non-linear evolution, using second-order perturbation theory, and relate this to the contamination from intrinsic alignments to the recently-measured cross-correlation between galaxy ellipticities and the CMB lensing potential. On the scales relevant for CMB lensing observations, the squeezed limit of the gravitational bispectrum dominates the correlation. Physically, the large-scale mode that sources CMB lensing modulates the small-scale power and hence the intrinsic ellipticity, due to non-linear evolution. We find that the angular cross-correlation from tidal torquing has a very similar scale dependence as in the linear alignment model, believed to be appropriate for elliptical galaxies. The amplitude of the cross-correlation is predicted to depend strongly on the formation redshift, being smaller for galaxies that formed at higher redshift when the bispectrum of the gravitational potential was smaller. Finally, we make simple forecasts for constraints on intrinsic alignments from the correlation of forthcoming cosmic shear measurements with current CMB lensing measurements. We note that cosmic variance can be significantly reduced in measurements of the difference in the intrinsic alignments for elliptical and spiral galaxies if these can be separated (e.g., using colour).
We present cosmological constraints from a cosmic shear analysis of the fourth data release of the Kilo-Degree Survey (KiDS-1000), doubling the survey area with nine-band optical and near-infrared photometry with respect to previous KiDS analyses. Adopting a spatially flat $Lambda$CDM model, we find $S_8 = sigma_8 (Omega_{rm m}/0.3)^{0.5} = 0.759^{+0.024}_{-0.021}$ for our fiducial analysis, which is in $3sigma$ tension with the prediction of the Planck Legacy analysis of the cosmic microwave background. We compare our fiducial COSEBIs (Complete Orthogonal Sets of E/B-Integrals) analysis with complementary analyses of the two-point shear correlation function and band power spectra, finding results to be in excellent agreement. We investigate the sensitivity of all three statistics to a number of measurement, astrophysical, and modelling systematics, finding our $S_8$ constraints to be robust and dominated by statistical errors. Our cosmological analysis of different divisions of the data pass the Bayesian internal consistency tests, with the exception of the second tomographic bin. As this bin encompasses low redshift galaxies, carrying insignificant levels of cosmological information, we find that our results are unchanged by the inclusion or exclusion of this sample.