# The Impact of Magnification and Size Bias on Weak Lensing Power Spectrum and Peak Statistics

## Abstract

The weak lensing power spectrum is a powerful tool to probe cosmological parameters. Additionally, lensing peak counts contain cosmological information beyond the power spectrum. Both of these statistics can be affected by the preferential selection of source galaxies in patches of the sky with high magnification, as well as by the dilution in the source galaxy surface density in such regions. If not accounted for, these biases introduce systematic errors for cosmological measurements. Here we quantify these systematic errors, using convergence maps from a suite of ray-tracing N-body simulations. At the cut-off magnitude of on-going and planned major weak lensing surveys, the logarithmic slope of the cumulative number counts is in the range . At , expected in the band for LSST, the inferred values of , and are biased by many (where denotes the marginalized error) and therefore the biases will need to be carefully modeled. We also find that the parameters are biased differently in the (, , ) parameter space when the power spectrum and when the peak counts are used. In particular, derived from the power spectrum is less affected than derived from peak counts, while the opposite is true for the best-constrained combination of (with from the power spectrum and from peak counts). This suggests that the combination of the power spectrum and peak counts can help mitigate the impact of magnification and size biases.

###### pacs:

PACS codes: 98.80.-k, 95.36.+x, 95.30.Sf, 98.62.Sb## I Introduction

By measuring the distortions of background galaxy shapes by foreground masses (galaxies, galaxy clusters, and large-scale structures), weak gravitational lensing (WL) surveys probe the mass density fluctuations throughout the cosmic span (see recent reviews by (1); (2); (3); (4)). WL observations, in conjunction with cosmological simulations, can be used to place precise constraints on cosmological parameters. Recent WL surveys, such as COSMOS (5) and CFHTLenS (6), have measured the shear power spectrum and have already placed useful constraints on (the matter density of the universe), (the amplitude of the primordial power spectrum on a scale of 8 comoving Mpc), and (the dark energy equation of state).

Because of the statistical nature of WL surveys, it is important to have an unbiased sample of source galaxies, fairly sampling the foreground density fluctuations across the sky. In this paper, we investigate possible sources of bias in flux-limited surveys, arising from a preferential selection of source galaxies in patches of the sky with high magnification, as well as by the dilution in the source galaxy surface density in such regions (known as magnification bias; hereafter MB). MB has been studied extensively in the past for its impact on galaxy–quasar and galaxy–galaxy correlation functions in 2D (7); (8); (9); (10); (11); (12); (13); (14); (15), and in 3D (16); (17); (18), and on the statistics of the Lyman- forest (19). An additional size bias (hereafter SB) can be present in surveys in which the selection of the source galaxies depends on their angular sizes. If not accounted for, these biases represent a systematic error for cosmological measurements. In the context of WL, the impact of MB and SB have been studied for the power spectrum (20); (21) and for high peaks caused by individual NFW halos (22).

Ref. (20) has shown that ignoring MB and SB in the
shear power spectrum can cause deviations in cosmological
parameter estimation for a DETF (23) Stage III experiment, such
as the Dark Energy Survey.^{1}^{2}^{3}

In this paper, we first show that MB is indeed significant for the power spectrum, extending earlier results (20) to explicitly compute the biases on cosmological parameters. We then focus on the impact of MB on peak counts. Lensing peaks were first considered as a cosmological probe in early ray-tracing simulations a decade ago (24). Peak counts have received increasing attention in recent years (25); (26); (27); (28); (29); (30); (31); (32); (33) as a way to access cosmological information from the strong non-Gaussianities in the lensing fields. In particular, these studies have shown that the number and height-distribution of peaks have high cosmological sensitivity, and can improve cosmological constraints by a factor of two, compared to using the power spectrum alone.

Peak counts are a simple and robust statistic, defined by recording local maxima in a 2D shear or convergence () map, smoothed by suitable filters. Ref. (28) investigated the physical origin of the individual peaks, by tracing their contributing light rays back in time across their N-body simulation boxes. They found that high peaks (with amplitudes 3.5, where is the r.m.s. of the convergence ) are typically created by individual massive halos. It has been shown that MB increases the signal-to-noise and therefore the total number of such peaks (22). By comparison, low peaks () are typically caused by a combination of (cosmology-independent) shape noise and a (cosmology-dependent) constellation of 4–8 lower-mass halos. These halos have masses of a few , and are offset by arcmin from the line of sight to the center of the peak. The low peaks are especially promising, as they carry the majority of cosmological information, and are relatively insensitive to baryonic cooling that affect the halo cores (32). We therefore extend the earlier results of (22) for high peaks, where increases in both the peak heights and number of high peaks were seen, to the low peaks, and to explicitly compute the biases on cosmological parameters.

To study the impact of MB, we build a simple numerical model to derive cosmological parameters (and their error bars) using either the power spectrum or peak counts measured in our simulations. We then apply magnification bias to a set of “true” convergence maps (which faithfully represent the projected dark matter distribution in a fiducial flat CDM cosmology) to create mock “biased” maps, mimicking an observed dataset. For each of these “biased” datasets, we find the best-fit set of the three cosmological parameters ( and ), using the “true” maps for the model fitting. Finally, we quantify the difference between the inferred cosmology and the true fiducial cosmology, as a function of the strength and sign of the magnification and size bias (determined by the slope of the galaxy luminosity function and the galaxy size distribution).

The rest of this paper is organized as follows: in § II, we introduce the formalism of magnification bias, and discuss its principal ingredient, the galaxy luminosity function. We then describe our computation methods in § III, including the convergence maps created with our ray-tracing N-body simulations, computing the power spectra and the peak distributions from these maps, determining the cosmology-dependence of these quantities, and finally applying biases to the maps to create mock observations. We present our main results in § IV, where we fit the mock data, and show that MB and SB will indeed alter the derived cosmological parameters by many . Finally, in § VI, we summarize our conclusions and the implications of this work.

## Ii Magnification Bias

Gravitational lensing causes a bias by modulating the apparent surface density of galaxies on the sky, through two competing effects (34). First, lensing can magnify (or demagnify) individual source galaxies in the background, increasing (or decreasing) their total flux. In a flux-limited WL surveys, some otherwise excluded faint galaxies can therefore make it into (or drop out of) the sample because of this (de)magnification. Second, a similar (de)magnification applies to the patch of the sky around the galaxy, geometrically diluting (or enhancing) the apparent surface density of galaxies in this region. These two effects counteract each other, and the net bias depends on the slope of the intrinsic (unlensed) galaxy luminosity function at the survey flux limit. In addition to these effects, lensing can increase (or decrease) the apparent angular size of spatially resolved individual galaxies. If either the survey selection, or a derived statistic such as WL shear, depends on the apparent size, then this can introduce an additional size bias.

### ii.1 Formalism

To quantify the effect of MB, we follow the discussion in Appendix A of (17). Including the effect of lensing on both the flux and on the geometrical surface density, we have the relation

(1) |

where is the observed (lensed) galaxy number density at position , as viewed by the observer, is the intrinsic (unlensed) galaxy number density, is the slope of the cumulative number counts evaluated at , and is the convergence. This equation assumes the weak lensing limit (), neglects the correspondingly small difference between lensed and unlensed directions on the sky, and also assumes that galaxy number density fluctuations on the angular scales of interest are small, as well. Under these assumptions, the above equation is valid to first order in , and . Finally, if we assume a survey with a sharp magnitude cut-off at :

(2) |

### ii.2 Galaxy Luminosity Functions and WL Surveys

The magnitude of MB depends on the galaxy luminosity function through the logarithmic slope . Observed luminosity functions are well described by a Schechter function (35),

(3) | |||||

where is the number density of galaxies with magnitude between and , is a characteristic number density (in Mpc), and is a characteristic magnitude. It consists of a power-law with index at the faint end and a exponential cut-off at the bright end. The cumulative galaxy number density can be written as

(4) |

Note that this equation holds at a given redshift. In
Fig. 1, we show calculated using
equation (2), as a function of cut-off magnitudes in the G,
R, I and Z bands at redshifts and . In this calculation,
we used the measurements of and by
(36); (37), which are all redshift-dependent.
^{4}

Table 1 lists the magnitude limits and the corresponding values of for several current and future WL surveys. We note that in order to measure the shape of the galaxies, it is necessary to adopt a brighter magnitude than for the point sources. For surveys where the magnitude limit was available only for point sources, we adopted a one magnitude brighter value for . For simplicity, for broad multi-band filters (R+I+Z), we have calculated using the central I band. Table 1 shows that surveys with have , assuming a mean redshift . While the effect of MB almost disappears () for galaxies at . Ref. (20) has shown that the –dependence is much weaker than the –dependence. For LSST, we expect the effective galaxy number density after applying lensing cuts to peak at a lower redshift () than the raw sample (). To illustrate the effect of MB, we adopt as our fiducial value, corresponding to for the conservative cut (see Figure 7 in ref. (38)). We will show in § IV that MB will significantly affect the power spectrum and the peak counts at this value.

Magnitude limit | ref | |||
---|---|---|---|---|

LSST | I 24.8 | 0.19 | 0.38 | (39) |

Euclid | R+I+Z 24.5 | 0.20 | 0.43 | (40) |

COSMOS | I 25 | 0.19 | 0.35 | (41) |

CFHTLS | I 24.7 | 0.19 | 0.39 | (42) |

DES | I 24.3 | 0.21 | 0.46 | (43) |

DUNE | R+I+Z 24.5 | 0.20 | 0.43 | (44) |

KiDS | R 25.2 (*) | 0.24 | 0.69 | (45) |

HSC | I 26.2 (*) | 0.18 | 0.32 | (46) |

### ii.3 Size Bias

If a survey has a cut in galaxy size , in addition to a flux cut, then equation (1) is modified to

(5) |

where in the case of a sharp cut, the new term is the logarithmic slope of the galaxy size distribution,

(6) |

This equation assumes that the size and flux cuts are independent, and also that the slopes and only weakly depend on and . Under these simple assumptions, the effects of size and magnification bias are equivalent, and only the combination matters. A more sophisticated treatment will eventually be necessary (and will depend on the details of the survey, including how galaxy sizes affect measurement errors). Here we simply note that at the limiting magnitudes of , the observed angular size distribution has a slope of (21). Therefore the additional effect of size bias is equivalent to increasing the value of by , i.e. the relevant fiducial value for LSST with a flux + size cut is changed from to . This means that the sign of the effect changes when we add size bias to magnification bias, as the effect of galaxy density dilution dominates over individual galaxy magnification.

## Iii Methodology

#### N-body Simulations

The N-body simulations and lensing maps were created with the Inspector Gadget lensing simulation pipeline on the New York Blue IBM BlueGene supercomputer. The N-body simulations are the same as the ones used in our earlier work (27); (28); (30); (32); (33). We refer readers to these papers for more detailed information. Here we briefly describe the basis of the simulations and the parameters used.

This work uses in total 35 different N-body simulations, covering 7 different cosmological models (1 fiducial cosmology plus 6 variations), each with 5 independent realizations of the same input primordial power spectrum. We chose our fiducial cosmological model to be , , Hubble constant , with a primordial matter power spectrum with and a spectral index of , using the best fit values from the seven-year results by the WMAP satellite (47). We vary each of the 3 parameters ( and ) one at a time (a higher value and a lower value than in the fiducial model), while keeping the other 2 parameters at the fiducial values. The 6 non-fiducial models have values of (while to keep a spatially flat universe), and . The combinations are listed in Table 2.

Fiducial | 0.798 | -1.0 | 0.26 |
---|---|---|---|

High- | 0.850 | -1.0 | 0.26 |

Low- | 0.750 | -1.0 | 0.26 |

High- | 0.798 | -0.8 | 0.26 |

Low- | 0.798 | -1.2 | 0.26 |

High- | 0.798 | -1.0 | 0.29 |

Low- | 0.798 | -1.0 | 0.23 |

The N-body simulations were generated using a modified version of the
Gadget-2 code^{5}^{6}

#### Ray-Tracing and Lensing Maps

To construct convergence maps, we perform ray-tracing. First, we output 3D boxes at redshifts corresponding to every 80 Mpc (comoving). We then divide the 3D box into many parallel pieces and project each slice onto a 2D plane perpendicular to the observer’s line of sight, using the the triangular shaped cloud (TSC) scheme (49). In the next step, we convert the surface density to gravitational potential at each plane using Poisson’s equation. Each 2D plane has a resolution of pixels. We then follow light rays from , traveling backward through the projection planes. The deflection angle and WL convergence and shear are calculated at each plane for each light ray. These depend on the first and second derivatives of the gravitational potential, respectively. Between the planes, the light rays travel in straight lines. Finally, for each of the 7 cosmological models, we create 1,000 convergence maps of 12 deg each in size. This is done by mixing simulations of different realizations, and randomly rotating and shifting the simulation data cubes.

We add galaxy ellipticity noise to our maps, due to variations in the intrinsic shapes of galaxies, and their random orientations on the sky. This shape noise is added to the raw convergence maps using a redshift-dependent expression for the noise in one component of the shear (50):

(7) |

For each pixel, we add drawn from a random Gaussian distribution centered at zero with variance (51)

(8) |

where is the number of galaxies per , and is the solid angle of a pixel in units of . In the case of LSST, we expect arcmin(38) for galaxies that are usable for shape measurements, and it follows that . This is much larger than the WL signal, whose r.m.s. value (at ) for noise-free maps is . To average out the random galaxy noise, we perform smoothing on individual maps with a Gaussian kernel:

(9) | |||||

(10) |

where is the smoothed value at pixel
and is the Gaussian kernel with a smoothing scale
arcmin.^{7}

For simplicity, we use only convergence maps for source galaxies at the single redshift , as the –dependence of MB has shown to be weak (20). Future work should employ tomography with multiple redshifts, and fold into the analysis the actual –distribution of the source galaxies. In total, we have 7,000 convergence maps; we call these the “true” maps, since they do not include any magnification bias. We use this set of maps to predict the cosmology-dependent observables (power spectra or peak counts), which will be described in detail in § III.0.5.

#### Power Spectra and Peak Counts

The power spectrum is the most widely used statistic in current WL surveys, and has already been shown to be affected significantly by MB (20). We revisit the impact of MB on the power spectrum in order to cross-check our simulation results, and to explicitly compute the resulting biases on the cosmological parameters.

We first compute the power spectra for spherical harmonic index in the range , with 1000 equally spaced (linear) bins. This covers the range of angles from our pixel size ( arcsec) to the linear size of our maps ( deg). In our previous work (27); (30) we compared our numerical power spectrum with the semi-analytical power spectrum obtained using the Limber approximation (52) and integrating the nonlinear 3D matter power spectrum along the line of sight (53). Our power spectrum loses power on large scales below due to our finite box size, and on small scales above due to spatial resolution; there is excellent agreement with the semi-analytic predictions between these scales.

Peak counting is done by simply scanning through the pixels on a convergence map, and identifying local maxima (i.e. pixels with a higher value of than its surrounding 8 pixels). We then record the number of peaks as a function of their central value.

#### Applying Bias to the Convergence Maps

On each of the 1000 maps in our fiducial cosmology, we apply different levels of MB, ranging from to 1.0, with a step size . To do this, on each fiducial map, we take into account the factor in eq. (1) and add when smoothing the map. Eq. (9) becomes (with dependence suppressed for and ):

(11) |

This is the smoothed at each pixel, weighted by the galaxy number densities modified by MB. Note that we assume the intrinsic (unlensed) galaxy number density to be a constant - this ignores the effects of shot noise arising from a discrete sampling of the field by a finite number of galaxies, as well as the clustering of galaxies. Other than applying MB, the same procedures are then followed to add noise, smooth the maps, count peaks, or compute power spectra, on the “bias” maps, as for the “true” maps.

#### Predictions in Other Cosmologies

In this subsection, we describe how we interpolate (and extrapolate) the peak counts and power spectra for other cosmologies, using our set of simulations in the 7 different cosmologies listed in Table 2.

First, for individual convergence maps, we histogram the peaks into 200 equally spaced bins ranging from to 0.19 (this choice for the number of bins will be justified in § V(vii) below). We then calculate the mean peak distribution (average of the 1,000 maps) in each of the 7 cosmology models. To predict the peak distribution for an arbitrary combination of cosmological parameters, we treat each bin individually, and use a Taylor expansion:

(12) | |||||

Here denotes the total number of peaks in the bin (=1, 2…200), averaged over 1000 maps. and are the differences of the desired parameters () from the fiducial parameters ().

The same method was followed for the power spectrum, by simply replacing the peak counts with , the total power in the bin.

In the body of our paper below, we chose to use the fiducial and the “high” models as defined in Table 2, to compute the cosmology derivatives in equation (12) by a simple finite difference. We call these “forward derivatives”. Given that we also have “low” models for each parameter, ideally we could use all three models to refine these predictions, either by including second-order terms in the Taylor expansion, or using two-sided linear derivatives. In practice, we chose to avoid a second–order expansion, in order to be able to perform an analytical minimization (see next subsection). We have attempted to use a two-sided derivative, but have found that this caused numerical problems (the discontinuity in the derivative can cause the fitting procedure, described below, to become stuck). We therefore use the forward derivatives in the bulk of this paper. We will discuss the differences in our results if we use “backward derivatives” instead in § V.

#### Finding the Best-Fit Cosmology

To fit a cosmology to one of our biased maps (or more generally to an arbitrary peak count distribution), we minimize a , defined as

(13) |

Here is the difference between the peak distribution in a given single map () and the model () in the bin, and is the unbiased estimator of the inverse covariance matrix (54); (55). Summation is implied over repeated indices , . We make the simple assumption that the peak counts depend linearly on the three parameters. It then becomes possible to write down an analytical solution to the best fitted parameters. By defining

(14) | |||||

(15) |

where ) is a three-component vector and denotes one of the three parameters, we can rewrite

(16) | |||||

(17) |

Setting , we obtain

(18) |

which is symmetric in and , and hence the two terms can be written combined as

(19) |

and the difference between the best fit and the fiducial model is simply

(20) |

To check these analytical calculations and to eliminate potential
numerical errors from matrix inversion, we also directly minimized
eq. (13) using the numerical scipy routine
“optimize.minimize”.^{8}

The same fitting procedure was performed using the power spectrum, by simply replacing the peak count with the power spectrum in the above equations. In the case of the power spectrum model (as for the peaks), we used the covariance matrix derived using noisy maps, to include the higher power at small induced by the galaxy shape noise. However, to measure the power spectrum derivatives with respect to cosmological parameters, we computed using the noiseless maps directly (since noise adds linearly). We choose to use the noisy maps directly, but only with , as cutting off at (corresponding to arcmin) is equivalent to smoothing but has the advantage of faster computation.

The above procedure, applied to each of the 1,000 individual “bias” maps, returns a set of 1,000 best-fit parameters for each specific value of . We then use the distribution of these best-fits to find the average bias in the cosmology parameters (corresponding to the mean best fit), confidence levels, and the goodness-of-fit values.

## Iv Results

### iv.1 Power Spectrum

The impact of MB on the power spectrum is illustrated in Fig. 2. The levels of bias we chose to show are , 0.4 and 0.8. The value is close to that expected in LSST; is the special case when MB effect disappears completely (); and corresponds to in (20), close to the value expected in the presence of an additional size bias. For comparison, we also show the impact on of varying each cosmological parameter.

For , the observations suffer a negative bias magnitude of . In this case, the effect of diluting a patch of sky wins over the number density increase due to magnification. At all bins, the power is reduced, as the result of the decreasing fluctuations. For , we have and expect the MB effect to be absent. This is verified by the lack of any difference between the power spectrum in the and the fiducial (unbiased) models, and merely serves as a test of our numerical code. For , the power is increased on all scales; this behavior has the opposite sign of the case, and is consistent with the expectations from . For cross-check, we calculate for using shear maps. Our results (Fig. 3) are very close to the ones obtained by ref. (20) (their Fig. 1) in the range (note that our case is equivalent to their case, as they also included the reduced shear correction). However, we noticed that the amplitude of is 10 times smaller than if we use convergence maps (as in this work).

### iv.2 Peak Counts

Fig. 4 shows the impact of MB on peak counts. For the pure MB case of , the height of any positive peak is reduced due to the negative overall bias. The case continues to show no effect from MB. Finally, for the MB + SB case of (or ), all peaks are boosted to a higher value, and consequently the whole distribution is shifted toward the right. The peak counts change in a direction opposite to the case, and with a larger amplitude, as expected. We note that for this large positive bias, the abundance of the (or ) peaks increases (as discussed in ref. (22)), but the number of the low peaks is reduced.

A positive MB effect () also reduces the total number of peaks (the number in brackets in Fig. 4). By directly comparing an example of the “bias” maps against its original “true” version, we found that out of the peaks in total, peaks disappeared after MB, while only new peaks were created. By visual examination of the maps, we found that peak disappearance and creation tends to happen in complex regions, where many peaks are interconnected through filament-like structures. As an illustration of this, in Fig. 5 we show a typical “high” peak. The shape of this peak is fairly round, likely due to one single massive halo. High peaks like this normally remain a peak after MB. In contrast, Fig. 6 shows a typical low peak that disappears after MB is applied. The original low peak merges into the neighboring, somewhat higher-amplitude peak at the lower left corner – this can be attributed to the lensing bias creating a “ridge” between the two original peaks. The opposite phenomenon happens when the overall MB is negative (), where we see an increase in total number of peaks, due to the bias “destroying” ridges and causing a net increase in the number of low peaks.

We have found that MB results in a monotonic increase or decrease for
all peaks before smoothing, depending on the sign of
. Fig. 7 shows the change in values
for all individual pixels, as well as for the peaks, for the
case. The peaks that survived the MB (the pixels that are peaks in
both “true” and “bias” maps) tend to have a smaller increase in
their value than other random pixels. We speculate that these
are the local dominating peaks that could not gain a higher value due
to the lack of higher peaks around them. Interestingly,
Fig. 7 also show a clear cutoff at below which no peaks are seen.^{9}

Fig. 4 shows that the changes caused by variations in cosmological parameters tend to be more symmetric in the two wings of the peak distribution. For example, at , we see fewer high– peaks, as well as fewer low– peaks. This shows that no single cosmological parameter can mimic the changes caused by MB – however, a linear combination of the three parameters may still resemble such change and can be degenerate with the effects of MB (as we will see below).

Examining the changes due to and , we see a clear degeneracy between these two parameters. This previously known issue (e.g. (57); (58); (27); (25)) is similar to that from cluster counts – both and can change the number of massive halos; therefore, we can obtain the same number of massive halos (hence the same peak distribution) for a higher value of , as long as we decrease . A product of the two parameters in the form of is much more tightly constrained by a fixed number of halos. The value of depends on the relevant mass scale being measured, and varies from 0.3 to 0.6 (57); (59); (60); (61); (6). From our error ellipse, we found for the power spectrum and for peak counts, by minimizing for the 1,000 fitted fiducial maps ( and are the differences between the fitted values for an individual map and the fiducial parameters).

### iv.3 Cosmological Parameters

We are now ready to show that without taking into account the effect of magnification bias, WL surveys can deliver cosmological parameters that are biased from the true values by many times their statistical error – for both the power spectrum and peak counts.

Fig. 8 shows the average deviation of fitted parameters using the power spectrum, in units of their standard deviation (). We have computed this cosmology bias for the range of . The standard deviation is calculated over the 1,000 fiducial maps. Each fitted parameter is marginalized over the other parameters, and scaled from our simulation (12 deg) to LSST’s planned sky coverage of 20,000 deg. The shaded region indicates the values of where the deviation of the derived parameter is within 1, 2 and 3 (dark to light), for galaxy density arcmin. For , (although interestingly, as shown in the figure, the bias is not monotonic in ) and at 1. We choose to plot instead of and individually, because the former is much more tightly constrained, as discussed in § IV.2.

In Fig. 9, we show the deviations of cosmological parameters inferred from the peak counts (). For , we find , times larger in magnitude than from the power spectrum; at 1, which, on the other hand, is much lower than from the power spectrum. For , we see deviations at similar magnitude but in opposite directions to the case, for both the power spectrum and peak counts.

The biases are again shown in two dimensions in Fig. 10, where the Monte Carlo error ellipses, enclosing 68% of the best-fits, are explicitly shown for the fiducial unbiased maps and biased maps (). In conclusion, WL observations in a survey as large as LSST will need to take MB into account, by including it in the modeling when fitting the observations. Combining information from both the power spectrum and the peak counts will be useful, as these two observables are impacted by MB in different ways, and their combination can help mitigate the biases. The value of (or other parameters describing higher-order lensing corrections) could be potentially additional parameters in a fitting procedure, simultaneously with the cosmological parameters. We expect that MB has a smaller impact on the current surveys, mainly due to their smaller sky coverage (e.g. COSMOS: 2 deg, CFHTLenS: 150 deg). After scaling by their sky coverage, we found the deviations to be of order for COSMOS, and for CFHTLenS.

The observed galaxy number density will also affect the level of MB. In Figures 8 and 9, we also show the parameter biases for and 45 arcmin. For the power spectrum, the slope near tends to be steeper for larger . This means deeper surveys with higher galaxy number densities (hence smaller galaxy noise) are more sensitive to MB when galaxy noise (eq. 8) is smaller. For peak counts, MB impacts the derived for shallower surveys ( arcmin) in opposite direction to surveys with higher number density.

## V Discussions

In this work, we made several assumptions and simplifications, which we must highlight here:

(i) We assumed the number of peaks and the power spectrum depends linearly on cosmology. For example, in our analysis, we used “forward” derivatives for , built with the fiducial and the three “high” cosmologies from Table 2 for a finite difference. We can also use the three “low” cosmologies to obtain “backward” finite-difference derivatives. When we do so, we find the resulting deviations to have similar magnitude (Table 3) to the ones from “forward” derivatives, except for a significantly lower value for for peak counts.

power spectrum | forward | backward | forward | backward |
---|---|---|---|---|

-21.0 | -23.2 | 61.6 | 67.0 | |

0.9 | 0.4 | -6.1 | -4.1 | |

7.6 | 8.6 | -28.3 | -30.9 | |

-25.0 | -27.3 | 36.0 | 37.4 | |

peak counts | forward | backward | forward | backward |

1.4 | 1.5 | -2.9 | -3.2 | |

-3.1 | -4.4 | 8.7 | 4.5 | |

-2.6 | -2.7 | 4.8 | 3.1 | |

-3.0 | -3.1 | 4.7 | 0.2 |

We also attempted to use a spline interpolation, using all three data-points for each parameter to describe the cosmology-dependence. This enables us to utilize all 7 cosmologies simultaneously, but we lose the advantage of the analytical method to obtain the best-fits (eq. 20). We used the numerical method to find the best-fits with spline interpolation, and found mean biases similar to those from linear interpolations. However, the error ellipses from spline interpolation were considerably smaller, and suspiciously coincident with our simulated range of model parameters. This is likely due to the spline tails that curve dramatically outside our parameter region, and hence artificially force the fit to stay within our simulated range of model parameters for each map. To solve this issue, we will need to have a larger grid of simulation parameters, which will also help us understand the dependence of peak counts on cosmology more accurately.

(ii) We used convergence maps only at a single redshift. This is motivated by the fact that the effect of MB depends weakly on . At low redshift, is mainly dependent on the slope of this power-law tail, and MB will have similar level of impact for all galaxies. For galaxies at higher redshift, , when redshifted to the rest frame of the galaxy, moves closer to the exponential part of the luminosity function, so we expect to increase to a larger value. A redshift-dependent correction to MB that folds in the correct -distribution of high- galaxies will eventually be necessary.

(iii) We ignored all instrumental and measurement errors. In reality, the point spread function (PSF) deconvolution and the measurement of galaxy shapes accurately is a difficult task, and has received thorough discussions (e.g.(62); (63)). Ref. (33) used simulated shear maps with realistic galaxy properties and has taken into account distortions from both the atmosphere and optical errors expected for LSST. They have shown that, though peak significance is reduced, the addition of these errors does not significantly degrade the cosmological constraints, compared to considering shape noise only. While our basic conclusion, that MB is significant, likely remains valid in the presence of such errors, the detailed modeling of MB will need to incorporate these additional sources of error.

(iv) In this paper, we choose to work with convergence maps, as they are computationally simpler. Using galaxies with sizes larger than the PSF, the convergence field can potentially be inferred by combining galaxy size and flux measurements, as lensing modifies these two quantities by different factors of and , respectively, in the weak lensing limit (64); (65); (66); (67). In current practice, reduced shear maps are obtained by measuring the shapes of individual galaxies. The observer can deduce the aperture mass (), a smoothed form of convergence, by applying a convolution over tangential components of shear (68). Ref. (31) has shown that both shear and convergence statistics give similar constraints when compared at the same scale, but once again, the impact of lensing bias should be modeled directly on the shear field.

(v) Although WL surveys may implement a sharp flux cut, the size bias is likely to be more complicated, with an effective weighting on galaxies that depends monotonically on their size, but in a gradual fashion, rather than a step function. In the idealized case of a sharp size cut, our analysis remains applicable, with a suitable re-interpretation of as a stand-in for , where is the logarithmic slope of the size distribution (at the size cut). In this simplified case, the bias induced by the size cut is likely larger than the one induced by the flux cut. For example, (21) showed that, for a survey with magnitude cut and size cut , the impact of MB becomes positive and , which is equivalent to our cases with . From Fig. 10, we see the derived parameters remain many away from the true parameters, but in the opposite direction. This demonstrates that size cut is likely to be important, and also that it is necessary in future work to investigate the effects of size bias in more detail.

(vi) We have tested the impact of MB on three parameters, , and . When additional cosmological parameters are considered (e.g. , , , ), the impact of MB may be more severe, since a combination involving the new parameters could mimic the MB better. In order to test this, we need to run more N-body simulations with other parameters varying to build a more complete cosmological model.

(vii) Optimizing the number of bins has not been the focus of this work. However, the choice of the number of bins has an effect on the error sizes. As shown in Fig. 11, for peak counts, the values of derived parameters and marginalized errors only converge at bins. Once the number of bins exceeds this value, we see a roughly constant plateau extending to bins (beyond which the results become unreliable, due to having too few realizations of maps and the sample covariance matrix becoming singular near 1000 bins). Therefore, we chose to use 200 bins in this work.

## Vi Summary

In this paper, we have studied the effect of magnification bias on peak statistics, using convergence maps from ray-tracing N-body simulations. Using maps in a suite of simulations, we can predict the convergence power spectrum or peak count distribution as a function of , and . Using this tool, we found the biases in cosmological parameters, when convergence maps in the fiducial cosmology, modified by magnification bias, were used to find the best-fit cosmology, without taking MB into account in the fits.

Near the flux limit of future WL survey, such as LSST, the galaxy number counts have a logarithmic slope of . This causes a bias in the inferred value of by and of by when using the power spectrum, and by for and for when using peak counts. These results are scaled to WL observations expected from LSST. However, for recent surveys, such as COSMOS and CFHTLenS, the deviations are generally negligible ( and , respectively), due to their smaller sky coverage.

We conclude that it is necessary that cosmological simulations consider MB effects, when they are used to match observations. We have found that inferred from the power spectrum is less impacted by MB, but peak count is a less biased method to infer . Future work on magnification biases should incorporate the many improvements we have emphasized that are necessary, including (i) the redshift-dependence of the bias, (ii) the impact on shear maps with realistic measurement errors and the peak statistics derived from these maps, (iii) more complex biases induced by the size-dependent measurement errors cut on galaxies, and (iv) additionally, the potential of using magnification bias and size bias as a signal to tighten the constraints on convergence field (64); (65); (66); (67). Our results suggest that lensing biases can be mitigated by combining the power spectrum and the peak counts, which produce biases in very different directions in cosmological parameter space.

###### Acknowledgements.

We thank Kevin Huffenberger for useful discussions. This research utilized resources at the New York Center for Computational Sciences, a cooperative effort between Brookhaven National Laboratory and Stony Brook University, supported in part by the State of New York. This work is supported in part by the U.S. Department of Energy under Contract No. DE-AC02-98CH10886 and by the NSF under grant AST-1210877. The simulations were created on the IBM Blue Gene/L and /P New York Blue computer and the maps were created analyzed on the LSST/Astro Linux cluster at BNL.### Footnotes

- http://www.darkenergysurvey.org
- http://www.lsst.org
- http://sci.esa.int/euclid
- Eq.1 in ref. (36) describes the redshift evolution of and . The parameters can be found in table 3 & 4 of ref. (36) for the G band, and in table 9 (“case 3”, with a constant ) of ref. (37) for the R, I, Z bands.
- http://www.mpa-garching.mpg.de/gadget/
- http://camb.info/
- We note that in previous papers of this series (27); (28); (30); (32), a different definition of was adopted. When using the more commonly used definition (eq. 10), our smoothing scale of arcmin is equivalent to their arcmin.
- http://scipy.org/
- This could potentially be a cosmological probe, in analogy with the cosmology-dependent minimum in the probability distribution of in random directions on the sky (56).

### References

- A. Refregier, ARA&A 41, 645 (2003), [arXiv:arXiv:astro-ph/0307212].
- P. Schneider, ArXiv Astrophysics e-prints (2005), [arXiv:arXiv:astro-ph/0509252].
- H. Hoekstra and B. Jain, Annual Review of Nuclear and Particle Science 58, 99 (2008), [arXiv:0805.0139].
- M. Bartelmann, Classical and Quantum Gravity 27, 233001 (2010).
- T. Schrabback et al., A&A 516, A63 (2010), [arXiv:0911.0053].
- M. Kilbinger et al., MNRAS 430, 2200 (2013), [arXiv:1212.3338].
- R. L. Webster, P. C. Hewett, M. E. Harding and G. A. Wegner, Nature 336, 358 (1988).
- W. Fugmann, A&A 204, 73 (1988).
- R. Narayan, ApJL 339, L53 (1989).
- P. Schneider, A&A 221, 221 (1989).
- J. Verner Villumsen, ArXiv Astrophysics e-prints (1995), [arXiv:arXiv:astro-ph/9512001].
- J. V. Villumsen, W. Freudling and L. N. da Costa, ApJ 481, 578 (1997), [arXiv:arXiv:astro-ph/9606084].
- R. Moessner, B. Jain and J. V. Villumsen, MNRAS 294, 291 (1998), [arXiv:arXiv:astro-ph/9708271].
- N. Kaiser, ApJ 498, 26 (1998), [arXiv:arXiv:astro-ph/9610120].
- M. Loverde, L. Hui and E. Gaztañaga, Phys. Rev. D75, 043519 (2007), [arXiv:arXiv:astro-ph/0611539].
- T. Matsubara, ApJL 537, L77 (2000), [arXiv:arXiv:astro-ph/0004392].
- L. Hui, E. Gaztañaga and M. Loverde, Phys. Rev. D76, 103502 (2007), [arXiv:0706.1071].
- L. Hui, E. Gaztañaga and M. Loverde, Phys. Rev. D77, 063526 (2008), [arXiv:0710.4191].
- M. Loverde, S. Marnerides, L. Hui, B. Ménard and A. Lidz, Phys. Rev. D82, 103507 (2010), [arXiv:1004.1165].
- F. Schmidt, E. Rozo, S. Dodelson, L. Hui and E. Sheldon, ApJ 702, 593 (2009), [arXiv:0904.4703].
- F. Schmidt, E. Rozo, S. Dodelson, L. Hui and E. Sheldon, Physical Review Letters 103, 051301 (2009), [arXiv:0904.4702].
- F. Schmidt and E. Rozo, ApJ 735, 119 (2011), [arXiv:1009.0757].
- A. Albrecht et al., arXiv:astro-ph/0609591.
- B. Jain and L. Van Waerbeke, ApJL 530, L1 (2000), [arXiv:astro-ph/9910459].
- J. P. Dietrich and J. Hartlap, MNRAS 402, 1049 (2010), [arXiv:0906.3512].
- M. Maturi, C. Angrick, F. Pace and M. Bartelmann, A&A 519, A23 (2010), [arXiv:0907.1849].
- J. M. Kratochvil, Z. Haiman and M. May, Phys. Rev. D81, 043519 (2010), [arXiv:0907.0486].
- X. Yang et al., Phys. Rev. D84, 043529 (2011), [arXiv:1109.6333].
- L. Marian, R. E. Smith, S. Hilbert and P. Schneider, MNRAS 423, 1711 (2012), [arXiv:1110.4635].
- J. M. Kratochvil et al., Phys. Rev. D85, 103513 (2012), [arXiv:1109.6334].
- S. Pires, A. Leonard and J.-L. Starck, MNRAS 423, 983 (2012), [arXiv:1203.2877].
- X. Yang, J. M. Kratochvil, K. Huffenberger, Z. Haiman and M. May, Phys. Rev. D87, 023511 (2013), [arXiv:1210.0608].
- D. Bard et al., ApJ 774, 49 (2013), [arXiv:1301.0830].
- E. L. Turner, J. P. Ostriker and J. R. Gott, III, ApJ 284, 1 (1984).
- P. Schechter, ApJ 203, 297 (1976).
- A. Gabasch et al., A&A 421, 41 (2004), [arXiv:arXiv:astro-ph/0403535].
- A. Gabasch et al., A&A 448, 101 (2006), [arXiv:arXiv:astro-ph/0510339].
- C. Chang et al., ArXiv e-prints (2013), [arXiv:1305.0793].
- LSST Science Collaboration et al., ArXiv e-prints (2009), [arXiv:0912.0201].
- R. Laureijs, ArXiv e-prints (2009), [arXiv:0912.0914].
- N. Scoville et al., ApJS 172, 38 (2007), [arXiv:arXiv:astro-ph/0612306].
- C. Heymans et al., MNRAS 427, 146 (2012), [arXiv:1210.0032].
- The Dark Energy Survey Collaboration, ArXiv Astrophysics e-prints (2005), [arXiv:arXiv:astro-ph/0510346].
- A. Refregier, Experimental Astronomy 23, 17 (2009), [arXiv:0802.2522].
- J. T. A. de Jong, G. A. Verdoes Kleijn, K. H. Kuijken and E. A. Valentijn, Experimental Astronomy 35, 25 (2013), [arXiv:1206.1254].
- M. Takada, Subaru Hyper Suprime-Cam Project, in American Institute of Physics Conference Series, edited by N. Kawai and S. Nagataki, , American Institute of Physics Conference Series Vol. 1279, pp. 120–127, 2010.
- E. Komatsu et al., ApJS 192, 18 (2011), [arXiv:1001.4538].
- A. Lewis, A. Challinor and A. Lasenby, ApJ 538, 473 (2000), [arXiv:arXiv:astro-ph/9911177].
- R. Hockney and J. Eastwood, Chap 4, 113 (1988).
- Y.-S. Song and L. Knox, Phys. Rev. D70, 063510 (2004), [arXiv:arXiv:astro-ph/0312175].
- L. van Waerbeke, MNRAS 313, 524 (2000), [arXiv:arXiv:astro-ph/9909160].
- D. N. Limber, ApJ 117, 134 (1953).
- R. E. Smith et al., MNRAS 341, 1311 (2003), [arXiv:arXiv:astro-ph/0207664].
- J. Hartlap, P. Simon and P. Schneider, A&A 464, 399 (2007), [arXiv:arXiv:astro-ph/0608064].
- T. W. Anderson, An Introduction to Multivariate Statistical Analysis, third ed. (Wiley, New York, NY, 2003).
- E. V. Linder, JCAP 3, 19 (2008), [arXiv:0711.0743].
- E. Rozo et al., ApJ 708, 645 (2010), [arXiv:0902.3702].
- S. Basilakos and M. Plionis, ApJL 714, L185 (2010), [arXiv:1003.2559].
- H. Mantz, K. Jacobs and K. Mecke, Journal of Statistical Mechanics: Theory and Experiment 12, 15 (2008).
- J. P. Henry, A. E. Evrard, H. Hoekstra, A. Babul and A. Mahdavi, ApJ 691, 1307 (2009), [arXiv:0809.3832].
- A. Vikhlinin et al., ApJ 692, 1060 (2009), [arXiv:0812.2720].
- N. Kaiser, G. Squires and T. Broadhurst, ApJ 449, 460 (1995), [arXiv:arXiv:astro-ph/9411005].
- H. Bonnet and Y. Mellier, A&A 303, 331 (1995).
- A. Vallinotto, S. Dodelson and P. Zhang, Phys. Rev. D84, 103004 (2011), [arXiv:1009.5590].
- F. Schmidt et al., ApJL 744, L22 (2012), [arXiv:1111.3679].
- B. Casaponsa et al., MNRAS 430, 2844 (2013), [arXiv:1209.1646].
- A. Heavens, J. Alsing and A. H. Jaffe, MNRAS 433, L6 (2013), [arXiv:1302.1584].
- P. Schneider, MNRAS 283, 837 (1996), [arXiv:arXiv:astro-ph/9601039].