Intrinsic Ellipticity Correlation Research Articles

A unified linear intrinsic alignment model for elliptical and disc galaxies and the resulting ellipticity spectra

Alignments of disc galaxies were thought to result from tidal torquing, where tidal field of the cosmic large-scale structure exert torquing moments onto dark matter haloes, determining their angular momentum and ultimately the orientation of galactic discs. In this model, resulting intrinsic ellipticity correlations are typically present on small scales, but neither observations nor simulations have found empirical evidence; instead, simulations point at the possibility that alignments of disc galaxies follow a similar alignment model as elliptical galaxies, but with a weaker alignment amplitude. In our article we make the case for the theory of linear alignments resulting from tidal distortions of the galactic disc, investigate the physical properties of this model and derive the resulting angular ellipticity spectra, as they would appear as a contribution to weak gravitational lensing in surveys such as Euclid’s. We discuss in detail on the statistical and physical properties of tidally induced alignments in disc galaxies, as they are relevant for mitigation of alignment contamination in weak lensing data, comment on the consistency between the alignment amplitude in spiral and elliptical galaxies and finally, estimate their observability with a Euclid-like survey.

Testing tidal alignment models for anisotropic correlations of halo ellipticities with N-body simulations

ABSTRACT There is a growing interest of using the intrinsic alignment (IA) of galaxy images as a tool to extract cosmological information complimentary to galaxy clustering analysis. Recently, Okumura & Taruya derived useful formulas for the intrinsic ellipticity–ellipticity correlation, the gravitational shear–intrinsic ellipticity correlation, and the velocity–intrinsic ellipticity correlation functions based on the linear alignment (LA) model. In this paper, using large-volume N-body simulations, we measure these alignment statistics for dark-matter haloes in real and redshift space and compare them to the LA and non-linear alignment model predictions. We find that anisotropic features of baryon acoustic oscillations in the IA statistics can be accurately predicted by our models. The anisotropy due to redshift-space distortions (RSDs) is also well described in the large-scale limit. Our results indicate that one can extract the cosmological information encoded in the IA through the Alcock–Paczynski and RSD effects.

Anisotropies of galaxy ellipticity correlations in real and redshift space: angular dependence in linear tidal alignment model

ABSTRACT Investigating intrinsic alignments (IAs) of galaxy shapes is important not only to constrain cosmological parameters unbiasedly from gravitational lensing but also to extract cosmological information complimentary to galaxy clustering analysis. We derive simple and useful formulas for the various IA statistics, including the intrinsic ellipticity–ellipticity correlation, the gravitational shear–intrinsic ellipticity correlation, and the velocity–intrinsic ellipticity correlation functions. The angular dependence of each statistic is explicitly given, namely the angle between the line-of-sight direction and the separation vector of two points. It thus allows us to analyse anisotropies of baryon acoustic oscillations encoded in the IA statistics, and we can extract the maximum cosmological information using the Alcock–Paczynski and redshift-space distortion effects. We also provide these formulas for the intrinsic ellipticities decomposed into E and B modes.

Environmental dependence of ellipticity correlation functions of intrinsic alignments

In this work we investigate the environmental dependence of the intrinsic ellipticity correlations in cosmic shear surveys. We use the quadratic and linear alignment model to describe the contributions by spiral and elliptical galaxies, respectively. The density field is in both cases described by a Gaussian random field and ellipticity correlation functions that are conditional on the environment of the galaxies are constructed by sampling random values for the tidal tensor and inertial tensor. The covariance of the Gaussian random process from which the tensor entries are drawn is decomposed by means of a spherical Fourier-Bessel transformation of the density field. The dependence on environment is modelled by the number of positive eigenvalues of the tidal tensor, which allows a differentiation between voids, sheets, filaments and superclusters. We find that elliptical galaxies align strongest in elongated structures such as sheets and filaments with an amplitude almost an order of magnitude higher compared to the alignment in clusters or voids. In contrast to this, spiral galaxies align equally strong in all environments. Cross-alignments between different environments are smaller than the respective auto-correlations subject to the Cauchy-Schwarz inequality which is an effective bound on their amplitude. Furthermore, we find misalignment between inertial and tidal tensor to be stronger in anisotropic regions compared to clusters or voids. While the imprint of weak lensing on galaxy ellipticities is agnostic about the environment, using environment information can help to distinguish between lensing and the intrinsic alignment signal.

Angular spectra of the intrinsic galaxy ellipticity field, their observability and their impact on lensing in tomographic surveys

This paper describes intrinsic ellipticity correlations between galaxies, their statistical properties, their observability with future surveys and their interference with weak gravitational lensing measurements. Using an angular-momentum-based, quadratic intrinsic alignment model we derive correlation functions of the ellipticity components and project them to yield the four non-zero angular ellipticity spectra |$C^\epsilon _E(\ell )$|⁠, |$C^\epsilon _B(\ell )$|⁠, |$C^\epsilon _C(\ell )$| and |$C^\epsilon _S(\ell )$| in their generalization to tomographic surveys. For a Euclid-like survey, these spectra would have amplitudes smaller than the weak lensing effect on non-linear structures, but would constitute an important systematics. Computing estimation biases for cosmological parameters derived from an alignment-contaminated survey suggests biases of +5σw for the dark energy equation of state parameter w, |$-20\sigma _{\Omega _{\rm m}}$| for the matter density Ωm and |$-12\sigma _{\sigma _8}$| for the spectrum normalization σ8. Intrinsic alignments yield a signal that is easily observable with a survey similar to Euclid: while not independent, significances for estimates of each of the four spectra reach values of tens of σ if weak lensing and shape noise are considered as noise sources, which suggests relative uncertainties on alignment parameters at the percent level, implying that galaxy alignment mechanisms can be investigated by future surveys.

Intrinsic alignments of galaxies in the Illustris simulation

We study intrinsic alignments (IA) of galaxy image shapes within the Illustris cosmic structure formation simulations. We investigate how IA correlations depend on observable galaxy properties such as stellar mass, apparent magnitude, redshift and photometric type, and on the employed shape measurement method. The correlations considered include the matter density–intrinsic ellipticity (mI), galaxy density–intrinsic ellipticity (dI), gravitational shear–intrinsic ellipticity (GI) and intrinsic ellipticity–intrinsic ellipticity (II) correlations. We find stronger correlations for more massive and more luminous galaxies, as well as for earlier photometric types, in agreement with observations. Moreover, the correlations significantly depend on the choice of shape estimator, even if calibrated to serve as unbiased shear estimators. In particular, shape estimators that down-weight the outer parts of galaxy images produce much weaker IA signals on intermediate and large scales than methods employing flat radial weights. The expected contribution of IA to the observed ellipticity correlation in tomographic cosmic shear surveys may be below one percent or several percent of the full signal depending on the details of the shape measurement method. A comparison of our results to a tidal alignment model indicates that such a model is able to reproduce the IA correlations well on intermediate and large scales, provided the effect of varying galaxy density is correctly taken into account. We also find that the GI contributions to the observed ellipticity correlations could be inferred directly from measurements of galaxy density–intrinsic ellipticity correlations, except on small scales, where systematic differences between mI and dI correlations are large.

3D weak lensing: Modified theories of gravity

Weak lensing (WL) promises to be a particularly sensitive probe of both the growth of large scale structure (LSS) as well as the fundamental relation between matter density perturbations and metric perturbations, thus providing a powerful tool with which we may constrain modified theories of gravity (MG) on cosmological scales. Future deep, wide-field WL surveys will provide an unprecedented opportunity to constrain deviations from General Relativity (GR). Employing a three-dimensional (3D) analysis based on the spherical Fourier-Bessel (sFB) expansion, we investigate the extent to which MG theories will be constrained by a typical 3D WL survey configuration including noise from the intrinsic ellipticity distribution $\sigma_{\epsilon}$ of source galaxies. Here we focus on two classes of screened theories of gravity: i) $f(R)$ chameleon models and ii) environmentally dependent dilaton models. We use one-loop perturbation theory combined with halo models in order to accurately model the evolution of matter power-spectrum with redshift in these theories. Using a Fisher information matrix based approach, we show that for an all-sky spectroscopic survey, the parameter $f_{R_0}$ can be constrained in the range $f_{R_0}< 5\times 10^{-6}(9\times 10^{-6})$ for $n=1(2)$ with a 3$\sigma$ confidence level. This can be achieved by using relatively low order angular harmonics $\ell<100$. Including higher order harmonics $\ell>100$ can further tighten the constraints, making them comparable to current solar-system constraints. We also employ a Principal Component Analysis (PCA) in order to study the parameter degeneracies in the MG parameters. Our results can trivially be extended to other MG theories, such as the K-mouflage models. The confusion from intrinsic ellipticity correlation and modification of the matter power-spectrum at small scale due to feedback mechanisms is briefly discussed.

Intrinsic galaxy shapes and alignments – II. Modelling the intrinsic alignment contamination of weak lensing surveys

Intrinsic galaxy alignments constitute the major astrophysical systematic of forthcoming weak gravitational lensing surveys but also yield unique insights into galaxy formation and evolution. We build analytic models for the distribution of galaxy shapes based on halo properties extracted from the Millennium Simulation, differentiating between early- and late-type galaxies as well as central galaxies and satellites. The resulting ellipticity correlations are investigated for their physical properties and compared to a suite of current observations. The best-faring model is then used to predict the intrinsic alignment contamination of planned weak lensing surveys. We find that late-type galaxy models generally have weak intrinsic ellipticity correlations, marginally increasing towards smaller galaxy separation and higher redshift. The signal for early-type models at fixed halo mass strongly increases by three orders of magnitude over two decades in galaxy separation, and by one order of magnitude from z=0 to z=2. The intrinsic alignment strength also depends strongly on halo mass, but not on galaxy luminosity at fixed mass, or galaxy number density in the environment. We identify models that are in good agreement with all observational data, except that all models over-predict alignments of faint early-type galaxies. The best model yields an intrinsic alignment contamination of a Euclid-like survey between 0.5-10% at z>0.6 and on angular scales larger than a few arcminutes. Cutting 20% of red foreground galaxies using observer-frame colours can suppress this contamination by up to a factor of two.

Intrinsic ellipticity correlations of galaxies: models, likelihoods and interplay with weak lensing

Subject of this paper are the statistical properties of ellipticity alignments between galaxies evoked by their coupled angular momenta. Starting from physical angular momentum models, we bridge the gap towards ellipticity correlations, ellipticity spectra and derived quantities such as aperture moments, comparing the intrinsic signals with those generated by gravitational lensing, with the projected galaxy sample of EUCLID in mind. We investigate the dependence of intrinsic ellipticity correlations on cosmological parameters and show that intrinsic ellipticity correlations give rise to non-Gaussian likelihoods as a result of nonlinear functional dependencies. Comparing intrinsic ellipticity spectra to weak lensing spectra we quantify the magnitude of their contaminating effect on the estimation of cosmological parameters and find that biases on dark energy parameters are very small in an angular-momentum based model in contrast to the linear alignment model commonly used. Finally, we quantify whether intrinsic ellipticities can be measured in the presence of the much stronger weak lensing induced ellipticity correlations, if prior knowledge on a cosmological model is assumed.

Evolution of intrinsic ellipticity correlations due to peculiar motion

Topic of this paper is the time-evolution of intrinsic correlations of galaxy ellipticities due to peculiar motion. In our model, the galaxy ellipticities are determined from the angular momentum of their host haloes, which can be computed from the fluctuations statistics of a Gaussian random field. Subsequent peculiar motion distorts the ellipticity field and causes changes in the ellipticity correlations. Using analogies between this problem of shifted ellipticity tensors and the displacements of polarisation tensors in gravitational lensing of the cosmic microwave background we compute E-mode and B-mode spectra of the time-evolved ellipticity field, where the displacements are modelled with first and second order Lagrangian perturbation theory. For EUCLID, ellipticity correlations are decreased on large multipoles l >1000, amounting to up to 10% in the E-mode spectrum C_E^epsilon(l) and up to 60% in the B-mode spectrum C_B^epsilon(l) at l~3000 due to the dispersing effect of peculiar motion. E/B-mode conversion in analogy to CMB-lensing is present but small. We conclude that distortions of the ellipticity field due to peculiar motion is not affecting the prediction of ellipticity models on the scales relevant for lensing in the case of EUCLID's galaxy distribution, but should affect larger scales for surveys at lower redshifts.

Self-calibration for three-point intrinsic alignment autocorrelations in weak lensing surveys

The weak lensing signal (cosmic shear) has been shown to be strongly contaminated by the various types of galaxy intrinsic alignment (IA) correlations, which poses a barrier to precision weak lensing measurements. The redshift dependence of the IA signal has been used at the 2-point level to reduce this contamination by only measuring cross-correlations between large redshift bins, which significantly reduces the galaxy intrinsic ellipticity - intrinsic ellipticity (II) correlation. A self-calibration technique based on the redshift dependences of the IA correlations has also been proposed as a means to remove the 2-point IA contamination from the lensing signal. We explore here the redshift dependences of the IA and lensing bispectra in order to propose a self-calibration of the IA auto-correlations at the 3-point level (i.e. GGI, GII, and III), which can be well understood without the assumption of any particular IA model. We find that future weak lensing surveys will be able to measure the distinctive IA redshift dependence over ranges of $|\Delta z^P|\le 0.2$. Using conservative estimates of photo-z accuracy, we describe the 3-point self-calibration technique for the total IA signal, which can be accomplished through lensing tomography of photo-z bin size $\sim 0.01$. We find that the 3-point self-calibration can function at the accuracy of the 2-point technique with modest constraints in redshift separation. This allows the 3-point IA auto-correlation self-calibration technique proposed here to significantly reduce the contamination of the IA contamination to the weak lensing bispectrum.

Self-calibration technique for three-point intrinsic alignment correlations in weak lensing surveys

The intrinsic alignment (IA) of galaxies has been shown to be a significant barrier to precision cosmic shear measurements. (Zhang, 2010, ApJ, 720, 1090) proposed a self-calibration technique for the power spectrum to calculate the induced gravitational shear-galaxy intrinsic ellipticity correlation (GI) in weak lensing surveys with photo-z measurements which is expected to reduce the IA contamination by at least a factor of 10 for currently proposed surveys. We confirm this using an independent analysis and propose an expansion to the self-calibration technique for the bispectrum in order to calculate the dominant IA gravitational shear-gravitational shear-intrinsic ellipticity correlation (GGI) contamination. We first establish an estimator to extract the galaxy density-density-intrinsic ellipticity (ggI) correlation from the galaxy ellipticity-density-density measurement for a photo-z galaxy sample. We then develop a relation between the GGI and ggI bispectra, which allows for the estimation and removal of the GGI correlation from the cosmic shear signal. We explore the performance of these two methods, compare to other possible sources of error, and show that the GGI self-calibration technique can potentially reduce the IA contamination by up to a factor of 5-10 for all but a few bin choices, thus reducing the contamination to the percent level. The self-calibration is less accurate for adjacent bins, but still allows for a factor of three reduction in the IA contamination. The self-calibration thus promises to be an efficient technique to isolate both the 2-point and 3-point intrinsic alignment signals from weak lensing measurements.

Erratum: Intrinsic alignment-lensing interference as a contaminant of cosmic shear [Phys. Rev. D70, 063526 (2004)

Cosmic shear surveys have great promise as tools for precision cosmology, but can be subject to systematic errors including intrinsic ellipticity correlations of the source galaxies. The intrinsic alignments are believed to be small for deep surveys, but this is based on intrinsic and lensing distortions being uncorrelated. Here we show that the gravitational lensing shear and intrinsic shear need not be independent: correlations between the tidal field and the intrinsic shear cause the intrinsic shear of nearby galaxies to be correlated with the gravitational shear acting on more distant galaxies. We estimate the magnitude of this effect for two simple intrinsic alignment models: one in which the galaxy ellipticity is linearly related to the tidal field, and one in which it is quadratic in the tidal field as suggested by tidal torque theory. The first model predicts a gravitational-intrinsic (GI) correlation that can be much greater than the intrinsic-intrinsic (II) correlation for broad redshift distributions, and that remains when galaxies pairs at similar redshifts are rejected. The second model, in its simplest form, predicts no gravitational-intrinsic correlation. In the first model and assuming a normalization consistent with recently claimed detections of intrinsic correlations we find that the GI correlation term can exceed the usual II term by >1 order of magnitude and the intrinsic correlation induced B-mode by 2 orders of magnitude. These interference effects can suppress the lensing power spectrum for a single broad redshift bin by of order ~10% at z_s=1 and ~30% at z_s=0.5.

SELF-CALIBRATION OF GRAVITATIONAL SHEAR-GALAXY INTRINSIC ELLIPTICITY CORRELATION IN WEAK LENSING SURVEYS

The galaxy intrinsic alignment is a severe challenge to precision cosmic shear measurement. We propose to self-calibrate the induced gravitational shear-galaxy intrinsic ellipticity correlation (the GI correlation, \citealt{Hirata04b}) in weak lensing surveys with photometric redshift measurement. (1) We propose a method to extract the intrinsic ellipticity-galaxy density cross correlation (I-g) from the galaxy ellipticity-density measurement in the same redshift bin. (2) We also find a generic scaling relation to convert the extracted I-g correlation to the demanded GI correlation. We perform concept study under simplified conditions and demonstrate its capability to significantly reduce the GI contamination. We discuss the impact of various complexities on the two key ingredients of the self-calibration technique, namely the method to extract the I-g correlation and the scaling relation between the I-g and the GI correlation. We expect none of them is likely able to completely invalidate the proposed self-calibration technique.

INTRINSIC ELLIPTICITY CORRELATION OF SDSS LUMINOUS RED GALAXIES AND MISALIGNMENT WITH THEIR HOST DARK MATTER HALOS

We investigate the orientation correlation of giant elliptical galaxies by measuring the intrinsic ellipticity correlation function of 83,773 luminous red galaxies (LRGs) at redshifts 0.16 -- 0.47 from the Sloan Digital Sky Survey. We have accurately determined the correlation up to 30 $h^{-1}$Mpc. Luminosity dependence of the ellipticity correlation is also detected although the error bars are large, while no evidence is found for its redshift evolution between $z=0.2$ and $z=0.4$. Then we use a cosmological $N$-body simulation to examine misalignment between the central LRGs and their parent dark matter halos. Central and satellite galaxies are assigned to simulated halos by employing a halo occupation distribution model for the LRGs. The ellipticity correlation is predicted to have the same shape as but an amplitude about 4 times higher than our observation if the central LRGs are perfectly aligned with their host halos. This indicates that the central LRG galaxies are preferentially but not perfectly aligned with their host halos. With the assumption that there is a misalignment angle between a central LRG and its host halo which follows a Gaussian distribution with a zero mean and a width $\sigma_\theta$, we obtain a tight constraint on the misalignment parameter, $\sigma_\theta={35.4}^{+4.0}_{-3.3}$ deg. This type of intrinsic ellipticity correlation, if not corrected, can lead to contamination at 5% level to the shear power spectrum in weak lensing surveys of limiting magnitude $R_{AB}=24.5$ if the source central galaxies follow the same misalignment distribution as the LRGs.

Sources of contamination to weak lensing three-point statistics: constraints fromN-body simulations

We investigate the impact of the observed correlation between a galaxies shape and its surrounding density field on the measurement of third order weak lensing shear statistics. Using numerical simulations, we estimate the systematic error contribution to a measurement of the third order moment of the aperture mass statistic (GGG) from three-point intrinsic ellipticity correlations (III), and the three-point coupling between the weak lensing shear experienced by distant galaxies and the shape of foreground galaxies (GGI and GII). We find that third-order weak lensing statistics are typically more strongly contaminated by these physical systematics compared to second-order shear measurements, contaminating the measured three-point signal for moderately deep surveys with a median redshift z_m ~ 0.7 by ~ 15%. It has been shown that accurate photometric redshifts will be crucial to correct for this effect, once a model and the redshift dependence of the effect can be accurately constrained. To this end we provide redshift-dependent fitting functions to our results and propose a new tool for the observational study of intrinsic galaxy alignments. For a shallow survey with z_m ~ 0.4 we find III to be an order of magnitude larger than the expected cosmological GGG shear signal. Compared to the two-point intrinsic ellipticity correlation which is similar in amplitude to the two-point shear signal at these survey depths, third order statistics therefore offer a promising new way to constrain models of intrinsic galaxy alignments. Early shallow data from the next generation of very wide weak lensing surveys will be optimal for this type of study.

Intrinsic alignment-lensing interference as a contaminant of cosmic shear

Cosmic shear surveys have great promise as tools for precision cosmology, but can be subject to systematic errors including intrinsic ellipticity correlations of the source galaxies. The intrinsic alignments are believed to be small for deep surveys, but this is based on intrinsic and lensing distortions being uncorrelated. Here we show that the gravitational lensing shear and intrinsic shear need not be independent: correlations between the tidal field and the intrinsic shear cause the intrinsic shear of nearby galaxies to be correlated with the gravitational shear acting on more distant galaxies. We estimate the magnitude of this effect for two simple intrinsic-alignment models: one in which the galaxy ellipticity is linearly related to the tidal field, and one in which it is quadratic in the tidal field as suggested by tidal torque theory. The first model predicts a gravitational-intrinsic (GI) correlation that can be much greater than the intrinsic-intrinsic (II) correlation for broad redshift distributions, and that remains when galaxies pairs at similar redshifts are rejected. The second model, in its simplest form, predicts no gravitational-intrinsic correlation. In the first model, and assuming a normalization consistent with recently claimed detections of intrinsic correlations, we find that the GI correlation term can exceed the usual II term by >1 order of magnitude and the intrinsic correlation induced B-mode by 2 orders of magnitude. These interference effects can suppress the lensing power spectrum for a single broad redshift bin by of order ∼10% at zs=1 and ∼30% at zs=0.5. We conclude that, depending on the intrinsic-alignment model, the GI correlation may be the dominant contaminant of the lensing signal and can even affect cross spectra between widely separated bins. We describe two ways to constrain this effect, one based on density-shear correlations and one based on scaling of the cross correlation tomography signal with redshift.

Weak lensing with COMBO-17: estimation and removal of intrinsic alignments

We estimate and remove the contamination of weak gravitational lensing measurements by the intrinsic alignment of close pairs of galaxies. We do this by investigating both the aperture mass B mode statistic, and the shear correlations of close and distant pairs of galaxies. These can be used to quantify non-lensing effects in weak lensing surveys. We reanalyse the {COMBO-17} survey, and study published results from the Red-sequence Cluster Survey and the {VIRMOS-DESCART} survey, concluding that the intrinsic alignment effect is at the lower end of the range of theoretical predictions. We also revisit this theoretical issue, and show that misalignment of baryon and halo angular momenta may be an important effect which can reduce the intrinsic ellipticity correlations estimated from numerical simulations to the level that we and the {SuperCOSMOS} survey observe. We re-examine the cosmological parameter estimation from the {COMBO-17} survey, using the shear correlation function, and now marginalizing over the Hubble constant. Assuming no evolution in galaxy clustering, and marginalizing over the intrinsic alignment signal, we find the mass clustering amplitude is reduced by 0.03 to ??8(??m/0.27)0.6= 0.71 +/- 0.11, where ??m is the matter density parameter. We consider the forthcoming {SuperNova/Acceleration} Probe {(SNAP)} wide weak lensing survey, and the {Canada-France-Hawaii} Telescope Legacy wide synoptic survey; we expect them to be contaminated on scales {\textgreater}1 arcmin by intrinsic alignments at the level of {\textasciitilde}1 and {\textasciitilde}2 per cent, respectively. Division of the {SNAP} survey for lensing tomography significantly increases the contamination in the lowest redshift bin to {\textasciitilde}7 per cent and possibly higher. Removal of the intrinsic alignment effect by the downweighting of nearby galaxy pairs will therefore be vital for {SNAP.}

Intrinsic correlation of halo ellipticity and its implications for large-scale weak lensing surveys

We use a large set of state-of-the-art cosmological N-body simulations [512(3) particles] to study the intrinsic ellipticity correlation functions of haloes. With the simulations of different resolutions, we find that the ellipticity correlations converge once the haloes have more than 160 members. For haloes with fewer members, the correlations are underestimated, and the underestimation amounts to a factor of 2 when the haloes have only 20 particles. After correcting for the resolution effects, we show that the ellipticity correlations of haloes in the bigger box L = 300 h(-1) Mpc) agree very well with those obtained in the smaller box (L = 100 h(-1) Mpc). Combining these results from the different simulation boxes, we present accurate fitting formulae for the ellipticity correlation function c(11)(r) and for the projected correlation functions Sigma(11)(r(p)) and Sigma(22) (r(p)) over three orders of magnitude in halo mass. The latter two functions are useful for predicting the contribution of the intrinsic correlations to deep lensing surveys. With reasonable assumptions for the redshift distribution of galaxies and for the mass of galaxies, we find that the intrinsic ellipticity correlation can contribute significantly not only to shallow surveys but also to deep surveys. Our results indicate that previous similar studies significantly underestimated this contribution for their limited simulation resolutions.