1 INTRODUCTION
Weak gravitational lensing is a powerful probe of dark matter in the universe (e.g. Kaiser & Squires Reference Kaiser and Squires1993). Following initial investigations by Blain (Reference Blain2002) and Morales (Reference Morales2006), de Burgh-Day et al. (Reference de Burgh-Day, Taylor, Webster and Hopkins2015) have developed a new method to measure the weak lensing signal in individual galaxies called Direct Shear Mapping (DSM). The primary scientific application considered by de Burgh-Day et al. (Reference de Burgh-Day, Taylor, Webster and Hopkins2015) is the measurement of mass and mass distribution in dark matter halos around individual low-redshift galaxies. In particular, since the dark matter halo properties can be measured for individual galaxies, DSM will enable the measurement of the dispersion in the galaxy luminous matter to dark matter ratio, as a function of other galaxy observables.
The possibility of measuring individual galaxy dark matter halo masses through DSM is an exciting prospect, however the measurement itself is challenging, and potentially observationally expensive. We have consequently developed the approach presented here for identifying the most robust candidates for such a measurement. We have also used this approach to investigate the possibility of measuring the scatter in the M *–Mh relation using shear statistics, and we estimate the size of the statistical sample that would be required to make this measurement.
In this paper, the probability of weak lensing shear has been estimated as a function of the redshifts of the source and lensing galaxies, and a catalogue of candidate galaxy pairs is selected from the Galaxy and Mass Assembly Phase 1 Survey (GAMA I) Data Release 2 (DR2) catalogue (Driver et al. Reference Driver2011, Liske et al., 2015). We also find that the distribution of shears in a galaxy sample is influenced by the relationship between stellar mass and halo circular velocity, and the scatter in this relation. With enough shear measurements, it may be possible to constrain this relationship, and to measure the scatter.
DSM uses spatially resolved velocity field information for an object to obtain a shear measurement from the velocity map. DSM assumes intrinsic rotational symmetry in the the velocity map, and searches for departures from this symmetry. This requires either radio data cubes or spatially resolved optical spectroscopy. To identify prospective targets, it is desirable to first obtain an estimate of the shear signal present in a galaxy.
While galaxy–galaxy lensing has been used to measure halo masses in the past, those studies stack many galaxy–galaxy pairs statistically, to obtain average halo masses (Brainerd, Blandford, & Smail Reference Brainerd, Blandford and Smail1996; Hudson et al. Reference Hudson, Gwyn, Dahle and Kaiser1998; Wilson et al. Reference Wilson, Kaiser, Luppino and Cowie2001; Mandelbaum et al. Reference Mandelbaum, Seljak, Kauffmann, Hirata and Brinkmann2006). In addition, measurements of galaxy halo shapes have been made from stacked galaxy–galaxy weak lensing measurements (Hoekstra, Yee, & Gladders Reference Hoekstra, Yee and Gladders2004; van Uitert et al. Reference van Uitert, Hoekstra, Schrabback, Gilbank, Gladders and Yee2012).
To test our target selection algorithm, the weak-lensing statistics of a sample of galaxies in the Galaxy and Mass Assembly Data Release 2 (GAMA-DR2) catalogue were investigated, using the stellar mass estimates from Taylor et al. (Reference Taylor2011). The purpose of our lensing frequency algorithm is to estimate of the distribution of shear signals present in a dataset. This algorithm enables novel measurements to be made, and will improve the success rate of any survey intended to measure weak lensing via the DSM method.
The rest of this paper is organised as follows. In Section 2, relevant weak lensing theory is described. In Section 3, a theoretical estimate of the probability of weak lensing at low redshift is made, following Mortlock & Webster (Reference Mortlock and Webster2000). In Section 4, the inputs, structure, and outputs of the lensing frequency algorithm are outlined. In Section 5, results of the application of the algorithm to a dataset obtained from GAMA-DR2 are presented, along with an investigation of the possibility of using shear statistics to measure the scatter in the M *–Mh relation. Conclusions and a summary are presented in Section 6. Throughout the paper, we assume a flat Concordance cosmology, with Ωm = 0.3, ΩΛ = 0.7, H 0 = 70 km s−1 Mpc−1, and h 70 = H 0/100 = 0.7.
2 WEAK LENSING
In this section, we will outline the relevant weak lensing theory. Gravitational lensing is the deflection of light from some source on its path to the observer by an intervening mass. The deflection angle is given by
$$\begin{equation}
{\bf \alpha }({\bf \xi }) = \int {\rm d}^2{\bf \xi }^\prime \,\frac{4\,G\,\Sigma ({\bf \xi }^\prime )}{c^2}\,\frac{{\bf \xi }-{\bf \xi }^\prime }{|{\bf \xi }-{\bf \xi }^\prime |^2},
\end{equation}$$
where G is the gravitational constant, Σ is the projected mass distribution of the deflector, c is the speed of light, and ξ is the distance from the deflector (i.e. the impact parameter). Since the angle of incidence of the light has been altered, the source will appear to be in a different location to its true position. The true position of the source can be found by solving the lens equation
$$\begin{equation}
{\bf \beta } = {\bf \theta } - \frac{D_{\rm ds}}{D_{\rm s}}{\bf \alpha },
\end{equation}$$
where θ is the apparent angular separation of the deflector and source, β is the true angular separation of the deflector and source, D ds is the angular diameter distance between the deflector and source, and D s is the angular diameter distance between the observer and source. The angular coordinates β and θ can be related to the corresponding physical coordinates in the source and lens planes as
$$\begin{eqnarray}
{\bf \eta } &=& D_{\rm s}{\bf \beta }\nonumber \\
{\bf \xi } &=& D_{\rm d}{\bf \theta }.
\end{eqnarray}$$
In the case of a circularly symmetric lens, and perfect alignment between the observer, lens and source, the lens equation can be solved to obtain the Einstein radius, a characteristic length scale
$$\begin{equation}
\theta _{\rm E}^2 = \frac{4\pi GM(<\theta _{\rm E}D_{\rm d})}{c^{2}}\frac{D_{\rm ds}}{D_{\rm s}D_{\rm d}},
\end{equation}$$
where M(< θED d) is the mass enclosed within the Einstein radius, and D d is the observer–deflector angular diameter distance.
In the weak lensing regime, the light from the source passes well outside the Einstein radius, and the source is singly-imaged. In this case, one can assume that so long as the length scale of the source is much less than that of the deflector, then the lensing will be linear. Thus, it can be represented by a first-order Taylor expansion, allowing Equation (2) to be re-expressed as a linear coordinate mapping between the lensed and un-lensed coordinate systems
$$\begin{equation}
{\bf \beta } = \mathbf {A}\,{\bf \theta },
\end{equation}$$
where A is the Jacobian of transformation;
$$\begin{equation}
\bf {A} = \frac{\partial \beta}{\partial \theta} = \left(\begin{array}{cc}1-\kappa -\gamma _1 & -\gamma _2\\
-\gamma _2 & 1-\kappa +\gamma _1\\
\end{array}\right).
\end{equation}$$
Here, κ is the convergence, and γ1 and γ2 are the two components of the shear vector γ = (γ1, γ2). Equation (6) can be inverted to obtain
$$\begin{equation}
\kappa = 1-(A_{11}+A_{22})/2
\end{equation}$$
$$\begin{equation}
\gamma _1 = -(A_{11}-A_{22})/2
\end{equation}$$
$$\begin{equation}
\gamma _2 = -A_{21} = -A_{12} .
\end{equation}$$
The shear vector γ can be rewritten in polar co-ordinates as a function of a shear magnitude, γ, and an angle, ϕ
$$\begin{equation}
\gamma _1 = \gamma \cos \phi ; \gamma _2 = \gamma \sin \phi ,
\end{equation}$$
where ϕ is the angle of the shear vector and 
$\gamma = \sqrt{\gamma _1^2+\gamma _2^2}$. In this work, we are interested in the value of γ for a given lens–source system, which is a function of the lens projected mass density and the lens–source angular separation.
We assume a Singular Isothermal Sphere (SIS) lens profile throughout this paper, whose projected surface density has the form
$$\begin{equation}
\Sigma ({\bf \xi }) = \frac{\sigma _v^2}{2G}\frac{1}{|{\bf \xi }|},
\end{equation}$$
where σv is the halo velocity dispersion. Although the SIS profile is a primitive lens model, it is sufficient for this work and allows for a simple shear estimation. The Einstein radius for an SIS lens is given by
$$\begin{equation}
\theta _{\rm E}^2 = \frac{4\pi \sigma _v^2}{c^2}\frac{D_{\rm ds}}{D_{\rm s}},
\end{equation}$$
and from e.g. Lasky & Fluke (Reference Lasky and Fluke2009), the shear components for an SIS are
$$\begin{equation}
\gamma _1 = \frac{D_{\rm d}\theta _{\rm E}}{2|{\bf \xi }|^3}(\xi _2^2-\xi _1^2);\,\,\gamma _2 = -\frac{D_{\rm d}\theta _{\rm E}}{|{\bf \xi }|^3}\xi _1\xi _2,
\end{equation}$$
so that
$$\begin{equation}
\gamma = \frac{D_{\rm d}\theta _{\rm E}}{2|{\bf \xi }|}.
\end{equation}$$
Thus, we have an expression for the shear magnitude as a function of lens–source projected separation and lens velocity dispersion, the latter of which can be related to the halo mass.
3 AN ESTIMATION OF THE PROBABILITY OF LENSING
In this section, we describe the process by which a theoretical estimate may be made of the probability of a given source being weakly lensed, as a function of the source redshift. We begin by introducing a similar calculation for strong lensing in analytically solvable cosmologies from Mortlock & Webster (Reference Mortlock and Webster2000). We then adapt this work to the weak lensing case, and for a more realistic Concordance cosmology. In the following sections, we will assume two different lens populations: a population of halos housing elliptical galaxies, as in Mortlock & Webster (Reference Mortlock and Webster2000), and a Press–Schechter (Press & Schechter Reference Press and Schechter1974) population of halos. In both cases, we assume an SIS halo. We will discuss the general steps for obtaining the expression for the weak lensing optical depth, and will then discuss the elliptical galaxy and Press–Schechter halo cases individually.
Given an estimate for the spatial distribution and mass of lensing galaxies in a given volume, it is possible to make an estimate of the distribution of shears in the volume. This can in turn be used to make an estimate of the probability distribution of weak shears across the sky, assuming the distribution of lensing galaxies is isotropic. Mortlock & Webster (Reference Mortlock and Webster2000) have used these arguments to make an estimate of the probability of lensing of quasars by elliptical galaxies for three simplified cosmologies. They define the lensing optical depth, τ, as the fraction of the source plane within which the lens equation has multiple solutions. It can be used as an estimator for strong lensing probability. The contribution to the total optical depth by one lensing galaxy is
$$\begin{equation}
\tau _{{\rm g}} = \frac{\pi \beta _{\rm crit}^{2}}{4\pi },
\end{equation}$$
i.e. the fraction of the sky covered by its lensing cross-section. Here, βcrit is the angular distance from the deflector where a background source transitions between being singly or multiply imaged.
In the case of weak lensing, rather than being interested in the region where the source is multiply imaged, we are interested in the region where the source is singly imaged, yet sheared sufficiently such that it is still a measurable effect. This region takes the form of an annulus about the lens (assuming a spherically symmetric lens), the inner bound of which is the Einstein radius, θE, and the outer bound of which is a function of some limiting shear value, γlim. Thus, τ is re-written as a function of a new area, a(γlim), defined as the area covered by this annulus:
$$\begin{equation}
\tau _{{\rm g}}(\gamma _{\rm lim}) = \frac{a(\gamma _{\rm lim})}{4\pi },
\end{equation}$$
and
$$\begin{equation}
a(\gamma _{\rm lim}) = \pi \left[\left(\frac{1}{2\gamma _{\rm lim}} - 1\right)\theta _{\rm E}\right] ^{2},
\end{equation}$$
and substituting in Equation (12) this becomes
$$\begin{equation}
\tau _{{\rm g}}(\gamma _{\rm lim}) = \pi ^{2} \left(\frac{1}{2\gamma _{\rm lim}} - 1\right)^{2}\left(\frac{\sigma }{c}\right)^{4}\left(\frac{D_{\rm ds}}{D_{\rm s}}\right)^{2},
\end{equation}$$
where σ is the velocity dispersion of the deflector, c is the speed of light, D ds is the deflector-source distance, and D s is the observer-source distance.
If it is assumed that the population of lensing objects are non-evolving and have uniform volume density at all redshifts (which is reasonable at low redshift), then the differential number of objects at redshift z with velocity dispersion σ is
$$\begin{equation}
\frac{{\rm d}^{2}N_{\rm d}}{{\rm d}z_{\rm d}{\rm d}\sigma } = \frac{{\rm d}V_{0}}{{\rm d}z_{\rm d}}\frac{{\rm d}n_{\rm d}}{{\rm d}\sigma }.
\end{equation}$$
Here, dV 0/dz is the co-moving volume element at redshift z,
$$\begin{equation}
\frac{{\rm d}V_{0}}{{\rm d}z} = \frac{c}{H_{0}}\frac{(1+z)^{2}D_{A}(z)^{2}}{E(z)}{\rm d}\Omega ,
\end{equation}$$
where c is the speed of light, DA(z) is the angular diameter distance at redshift z, and
$$\begin{equation}
E(z) = \sqrt{(1+z)^{3}\Omega _{{\rm M}}+\Omega _{\Lambda }},
\end{equation}$$
assuming a flat universe.
The optical depth to a redshift z s is then obtained by integrating over the optical depths of the entire population of deflectors up to z s
$$\begin{equation}
\tau (z_{{\rm s}},\gamma _{\rm lim}) = \int ^{z_{\rm s}}_0\int ^{\infty }_0 \frac{{\rm d}^{2}N_{\rm d}}{{\rm d}z_{\rm d}{\rm d}\sigma } \tau _{\rm g}(\gamma _{\rm lim}) {\rm d}\sigma {\rm d}z_{\rm d}.
\end{equation}$$
We will now discuss using this relation to estimate the weak lensing optical depth as a function of source redshift for a population of elliptical galaxies, n g, and a population of dark matter halos, n h.
3.1. A population of halos housing elliptical galaxies
Mortlock & Webster (Reference Mortlock and Webster2000) gives the local co-moving number density of elliptical galaxies as
$$\begin{equation}
\frac{{\rm d} n_{\rm g}}{{\rm d} \sigma _{\parallel }} = \frac{\delta n_{*}}{\sigma _{*}} \left(\frac{\sigma _{\parallel }}{\sigma _{*}}\right)^{\delta (1+\alpha )-1}\exp \left[-\left(\frac{\sigma _{\parallel }}{\sigma _{*}}\right)^\delta \right],
\end{equation}$$
where σ∥ is the observed line-of-sight velocity dispersion, α = −1.07 ± 0.05, n * = (0.0019 ± 0.003)h 370 Mpc−3, σ* = 225 ± 20 km s−1, and δ = 3.7 ± 1. The values of α and n * are drawn from Efstathiou, Ellis, & Peterson (Reference Efstathiou, Ellis and Peterson1988), and the values of σ* and δ from de Vaucouleurs & Olson (Reference de Vaucouleurs and Olson1982). Substituting this into Equation (22), the inner integral can be solved analytically, to give
$$\begin{equation}
\tau (z_{\rm s},\gamma _{\rm lim}) = C_{\rm g} \int ^{z_{\rm s}}_{0} \frac{(1+z_{\rm g})^{2}}{E(z_{\rm g})} \left(\frac{D_{\rm d}D_{\rm ds}}{D_{\rm s}}\right)^{2}{\rm d}z_{\rm g},
\end{equation}$$
where
$$\begin{equation}
C_{\rm g} = \pi ^{2}n_{*}\left(\frac{\sigma _\parallel }{c}\right)^{4} \left(\frac{1}{2\gamma _{\rm lim}}-1\right)^2 \Gamma \left(1+\alpha +\frac{4}{\delta }\right).
\end{equation}$$
For a given limiting shear value, γlim, to obtain the probability of measuring a shear of at leastthe limiting value, one must solve Equation (24) numerically as a function of the source redshift, z s.
The dashed grey line in Figure 1 shows the lensing optical depth as a function of source redshift for a limiting shear of γlim = 0.02 based on a population of elliptical galaxies. The coloured lines in Figure 1 are discussed in the next section.
The lensing optical depth as a function of source redshift and a limiting shear of γlim = 0.02. The coloured solid lines show the optical depths obtained when using a population of lenses drawn from a Press–Schechter halo mass function. The three lines show the effect of different minimum halo masses. The minimum halo mass is usually taken to be 1010M⊙ (e.g. in the Millennium simulation). The dashed grey line shows the optical depth obtained when using a population of lenses drawn from an elliptical galaxy population. It is reassuring to see that the Press–Schechter curve which best matches the dashed curve is that which uses the commonly used minimum halo mass of 1010M⊙.
3.2. A Press–Schechter halo population
The Press–Schechter mass function is given by
$$\begin{equation}
\frac{{\rm d}n}{{\rm d}M} = \sqrt{\frac{2}{\pi }}\frac{\rho _m}{M}\times -\frac{{\rm d}\ln \sigma (M)}{{\rm d}M}\frac{\delta _{crit}}{\sigma (M)}\exp ^{-\frac{\delta ^2_{crit}}{2\sigma ^2(M)}}
\end{equation}$$
and gives the co-moving number density of dark matter halos as a function of halo mass (Press & Schechter Reference Press and Schechter1974).
Substituting this into Equation (22),
$$\begin{equation}
\tau (z_{{\rm s}},\gamma _{\rm lim}) = \int ^{z_{\rm s}}_0\int ^{\infty }_{M_{\rm min}} \frac{{\rm d}^{2}n_{\rm d}}{{\rm d}z_{\rm d}{\rm d}M} \tau _{\rm d}(\gamma _{\rm lim}) {\rm d}M{\rm d}z_{\rm d},
\end{equation}$$
where M min is the mass of the smallest halo capable of containing a galaxy, and
$$\begin{equation}
\tau _{{\rm g}}(\gamma _{\rm lim}) = \pi ^{2} \eta \left(\frac{1}{2\gamma _{\rm lim}} - 1\right)^{2}\left(\frac{(M\,h_{70}^{-1})^{0.316}}{c}\right)^{4}\left(\frac{D_{\rm ds}}{D_{\rm s}}\right)^{2},
\end{equation}$$
where η = 6.14656 × 10−7 and we have used the relation
$$\begin{equation}
V_{\rm c} = 2.8 \times 10^{-2} \left(M\,h_{70}^{-1}\right) ^{0.316},
\end{equation}$$
from Klypin, Trujillo-Gomez, & Primack (Reference Klypin, Trujillo-Gomez and Primack2011). V c is the halo asymptotic circular velocity, and we have assumed V c ≃ σ. The lower bound of the integral over M in Equation (27) has been truncated at M min because we are only interested in halos large enough to contain at least one galaxy. The value of M min at z ≃ 0 is usually taken to be log10(M min/M⊙) ≃ 10 (Barkana & Loeb Reference Barkana and Loeb2001).
In this case, the integral over mass is not analytically solvable, however the two integrals are still separable, giving the following expression:
$$\begin{eqnarray}
\tau (z_{{\rm s}},\gamma _{\rm lim}) & = C_h \displaystyle\int ^{\infty }_{M_{\rm min}} \frac{{\rm d}n_h}{{\rm d}M} \left(M\,h_{70}^{-1}\right) ^{1.264} {\rm d}M\nonumber\\
& \times \displaystyle\int ^{z_{\rm s}}_0 \frac{(1+z_h)^{2}}{E(z_h)} \left(\frac{D_{\rm d}D_{\rm ds}}{D_{\rm s}}\right)^{2}{\rm d}z_h,
\end{eqnarray}$$
where
$$\begin{equation}
C_h = \pi ^2\eta ^2 \left(\frac{1}{2\gamma _{\rm lim}} - 1\right)^{2}\left(\frac{\sigma ^*}{c}\right)^{4},
\end{equation}$$
and σ* ≃ 321 km s−1 is a halo characteristic velocity dispersion, corresponding to a characteristic halo mass M* defined as where σ(M*) = δc(z), where δc(z) is the critical density and δc(z = 0) ≃ 1.686. At z = 0, log10(M*/M⊙) ≃ 13 (Barkana & Loeb Reference Barkana and Loeb2001).
To obtain the probability of a shear of at least γlim from halos, as a function of source redshift, the mass and redshift integrals in Equation (30) must be solved numerically. In order to illustrate the effect of varying the minimum halo mass, the mass integral in Equation (30) has been solved for three values of the minimum halo mass: log10(M min/M⊙) = 9.5,10, and 10.5. The results of this are shown by the three coloured curves in Figure 1. As expected, the probability curve which best matches that of an elliptical population is for the case log10(M min/M⊙) = 10. Not surprisingly, the higher the minimum halo mass the lower the probability of lensing, since the majority of halos are of lower mass.
From Figure 1, roughly 1 in 1 000 sources at z ~ 0.2 will be sheared by at least γ = 0.02. In the calculation of halo lensing probability, four major assumptions have been made:
• The lensing cross-sections of each galaxy in the population do not overlap. This assumption breaks down if: (1) The redshift is high, since the number of lenses contributing to the optical depth increases with increasing source redshift. (2) M min is small, since this leads to more lenses and an increasing cross-section. (3) γlim is small, since the cross-section of each lens is larger for a smaller limiting shear.
• Both assumed halo populations do not evolve with redshift. This assumption is only correct over small redshift ranges. The elliptical and halo populations used in Figure 1 are for z = 0.
• The Press–Schechter halo population is dependent on the minimum halo mass, i.e. the smallest halo which will house a galaxy. This is discussed in detail in Section 5.1.1.
• For this calculation, we have assumed a uniform distribution of lens redshifts, which is not the case in real data. If the total lens mass distribution is specified however, any inhomogeneity in the galaxy distribution will average out over a large sample of source galaxies. Thus, when comparing the results of our theoretical calculation to lensing probabilities in real data, it is acceptable to use the redshifts of the source galaxies in the sample and treat them as homogenous. Such a comparison is discussed in greater detail in Section 5.1.1.
In spite of minor limitations, the results shown in Figure 1 are strong motivation for a more through calculation and estimation of the probability of detecting weak lensing, since we are interested primarily in relatively low-redshift weak lensing, where the DSM technique will be useful.
4 THE LENSING FREQUENCY ALGORITHM
We now describe an algorithm for estimating the lensing signal per galaxy in a catalogue. We can use this algorithm to estimate the probability of any given object in the catalogue being lensed by some value, and to identify potentially suitable targets for lensing studies.
There are three factors which will determine the strength of the shear signal a lens imposes on a source: the mass of the lens, the angular separation of the lens and source, and the redshifts of the lens and source. Different combinations of these can lead to systems in which the lens appears large, and almost obscuring the source, but has a shear of the same magnitude as well separated systems with a small but dense lens. The simplest approach to estimating the strength of the shear present in a potential lens–source system is to simply inspect their separations and redshifts. Clearly, this is prone to misidentifications, since systems may be falsely rejected using this method if the lens is particularly dense, and systems may be falsely selected if the lens appears large and close to a background source, but is of a very low density.
In the lensing frequency algorithm, the shear probability is estimated utilising information about the stellar mass of the lens, and the redshifts and angular separations of each lens–source pair. There are a number of steps involved in this process, beginning with the stellar masses of the galaxies in the catalogue, assuming an SIS density distribution for each galaxy, estimating the halo circular velocities of each galaxy, and ultimately returning an estimate of the shear each object imposes on its nearest neighbours. These measurements can then be binned in shear to obtain an estimate of the probability distribution function for shear, and to flag particularly promising targets for follow up and direct measurement of the shear. In this investigation, we focus on galaxy–galaxy lensing only, however this process is equally applicable to galaxy-group lensing, provided good group masses are available.
4.1. Estimation of the circular velocity, V c
The most uncertain step in the estimation of the estimated shears for each object is the calculation of the asymptotic circular velocity of the lens, Vc. This step is important as it takes us from the stellar mass to the total mass of each lens. Three approaches to this calculation were used:
1. Use a power-law relation between halo mass and halo circular velocity obtained from the Bolshoi simulations (Klypin et al. Reference Klypin, Trujillo-Gomez and Primack2011). In this case, it is necessary to compute the halo masses from the stellar masses, which is done using the relation derived by Moster et al. (Reference Moster, Somerville, Maulbetsch, van den Bosch, Macciò, Naab and Oser2010).
2. Assume the circular velocity at the outermost regions of the disk is approximately equal to the circular velocity in the halo. In this case, the baryonic Tully–Fisher relation (McGaugh et al. Reference McGaugh, Schombert, Bothun and de Blok2000) can be used to obtain Vc from the known stellar masses of the objects.
3. Use an empirically derived relation between V c and σ0 from Courteau et al. (Reference Courteau, McDonald, Widrow and Holtzman2007). This relation has a significant scatter.
There are significant uncertainties in all of the methods discussed here, and the discrepancies in the results obtained are at times significant (in particular for very high and very low mass objects). However, there is no single best method, and each method approaches the problem from a different starting point. These three methods may bracket the reality, and so by utilising all three methods simultaneously a good representation of possible values is obtained. Therefore, the shear resulting from all three methods are computed.
The steps taken by the algorithm are as follows:
1. The objects are sorted by redshift and the angular size distance for each object is computed by integrating over the Friedman equation:
(32)
$$\begin{equation}
D = \frac{c}{H_0} \int [\Omega _{\rm m} (1+z)^{3} + \Omega _\Lambda ]^{-1} {\rm d}z.
\end{equation}$$
2. Beginning with the lowest-redshift galaxy, assume it is a lens, and perform the following steps:
a. Compute the angular separation between the object and its neighbours out to a specified projected radius θmax, excluding those at a lower redshift:
(33)where Δ[Rcos D] = (R 1 − R 2)cos (D 1 + D 2), ΔD = D 1 − D 2, and Ri, Di are an objects right ascension and declinationFootnote 1;
$$\begin{equation}
\theta _{\rm sep} = \sqrt{\Delta [R\cos D]^2 + \Delta [D]^2},
\end{equation}$$
b. Sort the neighbours by angular distance from the object;
c. Compute the halo mass from the stellar mass (Moster et al. Reference Moster, Somerville, Maulbetsch, van den Bosch, Macciò, Naab and Oser2010)
(34)where Mh is the halo mass, m is the stellar mass, (m/Mh)0 is a normalisation, M 1 is a characteristic mass where m(Mh)/Mh is equal to (m/Mh)0, and ϕ and η are two slopes which indicate the behaviour of the relation at the low and high mass ends. The values for the free parameters used in this work are the best-fit values from Moster et al. (Reference Moster, Somerville, Maulbetsch, van den Bosch, Macciò, Naab and Oser2010):
$$\begin{equation}
\hspace*{-4.0pt} \frac{m(M_h)}{M_h} = 2 \left(\frac{m}{M_h}\right)_0\left[ \left(\frac{M_h}{M_1}\right)^{-\phi } + \left(\frac{M_h}{M_1}\right)^\eta \right]^{-1},
\end{equation}$$
(35)
$$\begin{eqnarray}
\log {(M_1)} & =& 11.884^{+0.030}_{-0.023} \nonumber\\
\left(\frac{m}{M_h}\right)_0 & =& 0.02820^{+0.00061}_{-0.00053} \nonumber\\
\phi & =& 1.057^{+0.054}_{-0.046} \nonumber\\
\eta & =& 0.556^{+0.010}_{-0.004} ;
\end{eqnarray}$$
d. Compute the circular velocity in three different ways:
i. Using a power law relation derived from the Bolshoi simulation (Klypin et al. Reference Klypin, Trujillo-Gomez and Primack2011), which uses a ΛCDM cosmology:
(36)where M vir is the virial mass and it is assumed that M vir = Mh.
$$\begin{equation}
V_{\rm c} = 2.8 \times 10^{-2} M_{\rm vir} ^{0.316},
\end{equation}$$
ii. From the stellar mass using the baryonic Tully–Fisher relation (McGaugh et al. Reference McGaugh, Schombert, Bothun and de Blok2000)
(37)where A = 26.25 h 70−2 M⊙ km−4 s4.
$$\begin{equation}
V_{\rm c} = \left(\frac{M_*}{{A}}\right)^{0.25},
\end{equation}$$
iii. Using an empirically derived relation between V c and σ0 (Courteau et al. Reference Courteau, McDonald, Widrow and Holtzman2007):
(38)
$$\begin{equation}
V_{\rm c} = \sqrt{2}\sigma _0.
\end{equation}$$
These will henceforth be referred to as the ΛCDM, bTF, and σ0 methods respectively.
e. The Einstein radius of the object is computed for each of the three V c values (where it is assumed that due to the virial theorem one can write Vc ~ σv), and the shear for each V c value is computed;
3. Move to the next highest redshift galaxy, assume it is a lens, and repeat the above steps.
4. Continue iterating, increasing in redshift until the second-highest redshift has been reached.
5 APPLICATION TO THE GAMA SURVEY DATA RELEASE 2 CATALOGUE
The target selection procedure outlined in Section 4 has been applied to the GAMA I DR2 catalogue (Driver et al. Reference Driver2011, Liske et al., Reference Liske2015). The GAMA Survey is part of a larger project aiming to exploit the latest generation of ground-based and space-borne survey facilities to study cosmology and galaxy formation and evolution (Driver et al. Reference Driver2009). Phase I of the GAMA Survey is a magnitude limited spectroscopic survey measuring galaxy spectra and redshifts in three equatorial regions centred at 9h, 12h, and 14.5h (called G09, G12, and G15, respectively), each with an area of 12 × 4 deg2 (Baldry et al. Reference Baldry2010). The fields were observed to a limiting r-band apparent magnitude of r app < 19.4, r app < 19.8, and r app < 19.4 mag respectively. The target galaxies are distributed over a redshift range 0 < z ≲ 0.5 with a median redshift of z ≃ 0.17. The GAMA DR2 catalogue contains all GAMA I main survey objects down to r < 19.0 mag (for G09 and G12) and r < 19.4 mag (for G15) including spectral redshifts (Baldry et al. Reference Baldry2014). The catalogue contains a total of 72 225 objects, of which 71 599 have derived stellar masses (Taylor et al. Reference Taylor2011). GAMA survey data is available on the GAMA websiteFootnote 2.
The lensing frequency algorithm has been applied to the GAMA I DR2 catalogue in two ways:
1. As the entire catalogue (with some minor cuts detailed below) to investigate the probability of any object being sheared by at least some value, and to identify conceivably suitable targets for potential follow up observation and shear measurement with DSM;
2. As a smaller sample to investigate the possibility of measuring the scatter in the M *–Mh relation, and to determine what sample size of galaxies is required to make this measurement.
These two analyses are described in detail in the following subsections.
5.1. Identifying weak lensing candidates in the GAMA survey Catalogue
In this section, we look at the entire GAMA 1 DR2 catalogue (with minor cuts described in the next paragraph), to investigate the probability of detecting shear of at least γ for any given object in the catalogue, and identify candidate targets for follow up observation and shear measurement with DSM.
After removing objects with undefined or uncertain redshift, and selecting for stellar mass in the range 8 < log10(M */M⊙) < 12, a dataset of 69 434 objects was obtained. The data was passed to the lensing frequency algorithm in two ways:
• As the original set of galaxies (i.e. 69 434 objects). The shears present in the dataset were calculated via the three halo circular velocity estimates (ΛCDM, bTF, and σ0). The three M * − V c relations are shown in Figure 2 as dotted (ΛCDM), dashed (bTF), and solid (σ0) lines. A histogram of shears present in the GAMA DR2 sample is shown in Figure 3 as solid lines. The probability of a given galaxy being sheared by a particular value is obtained by normalising by the total number of galaxies in the sample.
Figure 2.The M * − V c relation for three methods of calculation. ‘K’ denotes the ΛCDM method (Klypin et al. Reference Klypin, Trujillo-Gomez and Primack2011), ‘M’ the bTF method (McGaugh et al. Reference McGaugh, Schombert, Bothun and de Blok2000), and ‘C’ the σ0 method (Courteau et al. Reference Courteau, McDonald, Widrow and Holtzman2007). The dotted, dashed, and solid lines show the relation with no scatter introduced. The blue, maroon, and green points show a synthetic dataset obtained by making 100 realisations of the original dataset, and introducing a scatter of σΛCDM = 0.15 dex, σbTF = 0.14 h −270 dex, and σσ0 = 0.08 dex for the ΛCDM (Moster et al. Reference Moster, Somerville, Maulbetsch, van den Bosch, Macciò, Naab and Oser2010), bTF (McGaugh et al. Reference McGaugh, Schombert, Bothun and de Blok2000), and σ0 methods respectively. The three methods agree well in the range 8 < log10(M */M⊙) < 12 (to within ~ 15%), where most galaxies are situated.
Figure 3.Histogram of estimated shears for galaxies in the GAMA survey. The solid, coarsely-binned lines correspond to the shears present in the original GAMA DR2 sample. The dotted lines correspond to the shears present in a synthetic dataset, obtained by producing 100 relisations of the original dataset, and introducing a scatter in the M * − V c relation. For every galaxy in the GAMA survey, there are 100 galaxies in the synthetic dataset. Therefore, for ease of comparison of the histograms of the synthetic and real datasets, the number of sheared objects in the synthetic dataset have been divided by 100. One can see that the three methods agree well within each of the real and synthetic datasets, with the synthetic datasets sitting slightly above the real datasets.
• As a larger (synthetic) population of galaxies, to account for the scatter in the M * − V c relation. The synthetic population is produced by making 100 realisations of each input galaxy, with a scatter introduced to the realisations (giving a total dataset of 7 029 800 objects). The scatter in the population arises from the intrinsic scatter in each M * − V c relation and so is different for each of the three relations used. A scatter of σΛCDM = 0.15 dex was introduced into Equation (34) for the ΛCDM method (Moster et al. Reference Moster, Somerville, Maulbetsch, van den Bosch, Macciò, Naab and Oser2010); a scatter of σbTF = 0.14 h −270 dex was introduced into Equation (37) for the bTF method (McGaugh et al. Reference McGaugh, Schombert, Bothun and de Blok2000); and a scatter of σσ0 = 0.08 dex was introduced into Equation (38) for the σ0 method. The σ0 values and associated scatter used in the σ0 method were obtained by fitting for the stellar masses using velocity dispersions obtained from the Sloan Digital Sky Survey Data Release 10 (SDSSDR10; Ahn et al. Reference Ahn2014). The M * − V c relation of the synthetic population of galaxies is shown in Figure 2 as blue (ΛCDM), maroon (bTF), and green (σ0) points. One can see that in the range 8 < log10(M */M⊙) < 12, where most galaxies are situated, the three methods agree well. The shears present in the synthetic dataset were computed for every point, and the resulting number of sheared objects, normalised by 100, are presented in Figure 3, as dotted lines. The synthetic data is normalised by 100 to allow for better comparison to the real data, since for every galaxy in the real dataset there are 100 galaxies in the synthetic dataset. One can see that the three methods agree well within within each of the real and synthetic datasets, with the synthetic datasets sitting slightly above the real datasets. The implications of this upward shift are discussed in Section 5.2. The probability of a given galaxy being sheared by a particular value is obtained by normalising by the total number of galaxies in the sample.
5.1.1. Discussion
The number of galaxies with shears above the cut-off value of 0.02 varied between the three methods of calculating V c, with the ΛCDM method giving the largest estimate of measurable shears in the data, and the bTF method giving the smallest estimate. The probability of a shear of at least γ = 0.02 was P(⩾ γ) = 0.018, 0.005, and 0.007 for the ΛCDM, bTF, and σ0 methods respectively. The bTF method was chosen for selecting candidate targets because it gives the smallest shear probability, and so will be the least likely to overestimate the shears present in a sample. Using the bTF method, the number of objects in GAMA with an estimated shear of γ ⩾ 0.02 is 393. The bTF method always estimates a smaller shear than the ΛCDM and σ0 methods for any given object, so any of the 393 objects selected by the bTF method would also have been selected by the other two methods (although of course the converse is not true).
These objects were extracted from the sample, and matched to objects in the SDSS DR10 catalogue. The lens–source pairs in the sample were ranked by each of the authors by eye, based on their morphology, surface brightness, inclination angle, separation, and environment. The motivation for selecting on this criterion stems from the requirements of the DSM algorithm; that the source galaxy be undisturbed and stably-rotating, neither face-on nor edge-on and bright enough for observation with an Integral Field Unit (IFU). For a more detailed discussion of the requirements of the DSM algorithm, see de Burgh-Day et al. (Reference de Burgh-Day, Taylor, Webster and Hopkins2015). Two highly-ranked example lens–source pairs are shown in Figure 4.
Two example galaxy pairs from the GAMA DR2 Sample identified with the target selection algorithm. In both cases, the galaxy in the crosshairs is the source galaxy (i.e. the galaxy being lensed). The left-hand images show thumbnails of the galaxy pairs from the SDSS DR10 Finding Chart Tool, while the middle and right-hand images show the J-band images from the United Kingdom infrared telescope Infrared Deep Survey (UKIDSS), with the residuals from 2D Sérsic fits, taken from GAMA’s online Single Object Viewer tool. The top pair is at RA = 213.705 deg, DEC = 1.623 deg, and has an estimated shear of γ = 0.023. The lens and source redshifts are z = 0.128 and z = 0.186 respectively. The bottom pair is at RA = 213.705 deg, DEC = 1.623 deg, and has an estimated shear of γ = 0.053. The lens and source redshifts are z = 0.088 and z = 0.190 respectively.
Figure 5 shows the distribution of galaxies with significant estimated shear, with lens stellar mass plotted as a function of lens redshift. As might be expected, the selection favours more massive lensing galaxies, with significant numbers of candidates at all redshifts. The quality of the shear measurements however depends on the characteristics of the source galaxy.
Lens stellar mass as a function of lens redshift for the GAMA DR2 sample, showing the distribution of shears with these parameters. Grey contours show the density of all galaxies in the sample. Blue squares show galaxies with a shear in the range 0.01 < γ < 0.05, and red circles show those with a shear in the range γ > 0.5. The shears were estimated using the bTF method.
The normalised cumulative sum was taken of the data shown in Figure 3, the results of which are shown in Figure 6, giving the probability of measuring a shear of at least γlim. One can again see the upward shift in the distributions for synthetic data, relative to those for the real data. Also shown are two dashed lines corresponding to our theoretical estimate of the probability of lensing from Section 3. The probability of at least a given shear being observed in a galaxy increases as the minimum shear decreases. Since the lensing cross-section of each galaxy goes as γ−2lim, this is not surprising, as smaller minimum shears rapidly increase the fraction of the sky in which measurable lensing will occur. The method which gives the lowest overall probability is the bTF method. Since we expect to be able to measure shears as small as γ = 0.02 with DSM, we can expect to have at least a one in 200 chance of measuring shear in a randomly chosen galaxy (from the bTF method), and up to a roughly one in 50 chance (from the ΛCDM method).
Estimated probability of measuring a shear of at least γ for galaxies in the GAMA survey. The dotted lines are the probabilities with scatter introduced into the M * − V c relation (i.e. the synthetic dataset), and the solid lines are the probabilities with no scatter introduced. Again, one can see that the synthetic datasets sit higher than the real datasets. The dashed lines are included for comparison purposes, and correspond to the probabilities as a function of γ derived from the theoretical calculation of probability in Section 3, using the redshifts of the sheared objects in the GAMA survey.
It is obvious that the ΛCDM method results in much higher probabilities than the bTF and σ0 methods. It is interesting to note that the ΛCDM method is the one which uses only simulation and theoretically-derived relations to obtain values of V c from M *. In contrast, the bTF and σ0 methods utilise empirically derived relationships between V c and M *. One can easily identify the origin of the higher numbers of lensed galaxies from the ΛCDM method; the break in the M * − V c relation for this method. While giving good agreement in the intermediate mass range, this results in a much larger corresponding halo circular velocity. These higher V c galaxies will have a much larger lensing cross section, resulting in a larger number of lensed objects. Hence, while there are relatively few galaxies with log10(M */M⊙) > 12, they contribute strongly to the total lensing probability. This raises the question: how does this over-abundance of higher-V c galaxies affect the results of ΛCDM-based cosmological simulations? Further consideration of this topic is beyond the scope of this paper, and is left to future work.
As mentioned already, in Figure 6, we have included two dashed lines corresponding to our theoretical estimate of the probability of lensing from Section 3. The probability calculated in Section 3 is the probability that a source at redshift zs is lensed by at least some value γlim. The probability obtained from the lensing frequency algorithm however is the probability that a single galaxy from the GAMA DR2 sample, chosen at random, is lensed by at least γlim. In order to compare the two, we need to calculate the theoretical probability that one galaxy, chosen at random from a GAMA-like population of galaxies, is lensed by at least γlim. From Section 3, the probability that some object in the GAMA DR2 sample is lensed by at least γlim will be dependent on its redshift only (once a lens mass distribution has been assumed). Thus, the probability that any galaxy in the sample has been lensed by at least γlim will be the sum of the individual probabilities of each one having been lensed. Then, the probability that one galaxy chosen at random has been lensed by at least γlim will be the probability of any having been lensed, divided by the total number of galaxies in the sample. Figure 7 shows the histogram of the redshifts of the GAMA DR2 sample after a redshift quality cut has been made. While the distribution of redshifts in the GAMA sample is not uniform, it is sufficiently smoothly varying that the assumption made in Section 3 (that the lenses are uniformly distributed in redshift) will suffice. Thus, we have arrived at a method by which we can compare the results of Section 3 and Section 5.1.
Distribution of redshifts in the GAMA sample, after a redshift quality cut has been made. The mean and median of this distribution are 0.18 (solid white line) and 0.17 (dashed white line) respectively, and the range of the distribution is 0 ≲ z ≲ 0.65.
From Section 3, the theoretical probability of measuring a shear of at least γ = 0.02 in any one galaxy in the GAMA survey is P(⩾ γ) ≃ 0.002 and 0.003 for the halos housing elliptical galaxies and Press–Schechter halo populations respectively. Comparing to the probabilities obtained from the GAMA DR2 sample (P(⩾ γ) = 0.018, 0.005, and 0.007 for the ΛCDM, bTF, and σ0 methods respectively), we see the difference is a factor of ~ 6 for the ΛCDM method, and a factor of ~ 2 for the bTF and σ0 methods. This is an acceptable level of agreement, given the assumptions made in the calculations in Section 4. The level of agreement between the two approaches can be seen in Figure 6.
The theoretical estimate of P(> γlim) drawn from the Press–Schechter halo population is dependent on the value of the minimum halo mass in which a galaxy will form. As can be seen in Figure 1, a smaller (larger) value of M min will result in a higher (lower) probability of measuring a shear of at least γlim as a function of redshift. It can be seen in Figure 8 that adjusting M min will alter the level of agreement between the theoretical estimate and the probabilities derived from the GAMA DR2 sample. Each dashed grey line represents a minimum halo mass incremented by 0.5 log10(M⊙). The line with the highest corresponding probability is that with log10(M min/M⊙) = 8, and that with the lowest corresponding probability log10(M min/M⊙) = 12. The probabilities from the lensing frequency algorithm are as in Figure 6. It is interesting to note that for the bTF and σ0 methods, the best agreement with the theoretical approach is obtained for a minimum halo mass in the range 9 ≲ log10(M min/M⊙) ≲ 10, whereas the currently accepted value is log10(M min/M⊙) = 10.
Comparison of the results from the lensing frequency algorithm with the theoretically derived Press–Schechter probabilities, with varying M min. The highest theoretically-derived probability corresponds to a minimum halo mass of log10(M min/M⊙) = 8, while the lowest corresponds to log10(M min/M⊙) = 12. The results of the lensing frequency algorithm are plotted as in Figure 6. The bTF and σ0 methods best agree with a Press–Schechter halo population with a minimum halo mass in the range 9 ≲ log10(M min/M⊙) ≲ 10.
5.2. Measuring the scatter in the M *–Mh relation
Since halos of given mass can have different halo concentrations, spin parameters, and merger histories, we expect them to house galaxies with a range of masses. This manifests as a scatter in the stellar mass to halo mass relation. The value of the scatter in the M *–Mh relation is not well constrained, as directly measuring the masses of galaxies and their host halos is not trivial. Abundance matching techniques can be useful for describing the relationship between stellar mass and halo mass, however they cannot constrain the scatter in the M *–Mh relation (Conroy & Wechsler Reference Conroy and Wechsler2009). Leauthaud et al. (Reference Leauthaud2012) used traditional weak lensing techniques to constrain the M *–Mh relation, however as with abundance matching they are unable to place any meaningful constraints on the scatter in the relation. An alternative approach has been to utilise satellite galaxies to constrain the galaxy luminosity-halo mass relation. Early work involved stacking the central galaxies to obtain a statistical measure of the kinematics of the satellite galaxies (Erickson, Gottesman, & Hunter Reference Erickson, Gottesman and Hunter1987; Zaritsky et al. Reference Zaritsky, Smith, Frenk and White1993; Zaritsky & White Reference Zaritsky and White1994; Zaritsky et al. Reference Zaritsky, Smith, Frenk and White1997). However, recent work has avoided the need for stacking, largely by utilising the larger datasets offered by the 2dF Galaxy Redshift Survey (2dFGRS; Colless et al. Reference Colless2001) and SDSS (McKay et al. Reference McKay2002; Brainerd & Specian Reference Brainerd and Specian2003; Prada et al. Reference Prada2003; van den Bosch, Yang, & Mo Reference van den Bosch, Yang, Mo, Dettmar, Klein and Salucci2004a; van den Bosch et al. Reference van den Bosch, Norberg, Mo and Yang2004b; Conroy et al. Reference Conroy2005, Reference Conroy2007; More et al. Reference More, van den Bosch, Cacciato, Mo, Yang and Li2009). These techniques do not utilise weak lensing, and do not attempt to measure the scatter in the M *–Mh relation.
In this section, we test whether DSM measurements of a population of galaxies could be used to measure the scatter in the M *–Mh relation. This is a novel measurement which is difficult with traditional weak lensing techniques, but is made possible with DSM because it can measure individual shears with far greater accuracy around individual galaxies. We then briefly describe how to fit for the scatter in the M *–Mh relation, and measure the scatter in a set of simulated shear datasets with a known scatter incorporated.
As was noted in Section 5.1, and can be seen in Figures 3 and 6, scatter in the M *–Mh relation results in a shift in the distribution of shears present in a population of galaxies towards larger shears. If the measurement error in the shear is sufficiently small, or the population of galaxies with known stellar masses is large enough, it is possible to measure the scatter in the M *–Mh relation by comparing the distribution of shears to those obtained from a M *–Mh relation with zero scatter.
To perform this measurement with DSM, velocity maps from from an intermediate redshift survey would be required. The velocity maps can be obtained from several components of the galaxy, such as HI in radio wavelengths, or IFU observations of stellar velocities and gas emission in optical wavelengths. IFU maps obtained from bright, Hα emitting galaxies are the most practical and easily obtainable in the immediate term however, and so we will focus on these galaxies in this analysis.
We assume that the scatter in the M *–Mh relation is lognormal, which is the standard form assumed in the literature (e.g. Behroozi, Conroy, & Wechsler Reference Behroozi, Conroy and Wechsler2010, Moster et al. Reference Moster, Somerville, Maulbetsch, van den Bosch, Macciò, Naab and Oser2010, Behroozi, Wechsler, & Conroy Reference Behroozi, Wechsler and Conroy2013). The mean of a lognormal distribution is a function of the size of the scatter in the distribution, and is given by exp (μ + σ2/2), where σ is the scatter in the distribution, and μ is the natural logarithm of the mean of the underlying normal distribution. It can be seen that the mean increases with increasing scatter. Thus, if a lognormal scatter is introduced into the M *–Mh relation for a population of galaxies, the result is a larger number of galaxies with higher mass halos. It is this property of the lognormal scatter in the M *–Mh relation which leads to an increased probability of larger shears.
To compare the shear distributions obtained from populations with and without scatter in the M *–Mh relation, the lensing frequency algorithm was applied to a selection of mock catalogues generated from a subsample of the complete GAMA DR2 sample. This subsample was obtained by performing the same cuts on the GAMA DR2 as in the previous section, along with the following additional cuts: 0.1 < z < 0.15, r < 17.5, and keeping only galaxies with H-α emission. The resulting dataset contained 2 861 galaxies. When computing shears in the mock catalogues, a maximum lens–source separation of D max = 0.2 Mpc was used for step (2a) of the lensing frequency algorithm.
A successful measurement of the scatter in the M *–Mh relation will be limited by survey sample size and the size of the scatter. To investigate the range in which these parameters would enable a measurement of the scatter, a set of simulated datasets of varying size was created, with a range of scatters in the M *–Mh relation.
To produce the mock catalogues, larger populations of galaxies were generated from the subsample of GAMA DR2 galaxies by producing Monte-Carlo realisations of each galaxy. The resulting datasets contained N= 1 000, 15 000, 50 000, and 150 000 galaxies. The distribution of shears for a tight M *–Mh relation was computed for the simulated datasets. A scatter, 
$\sigma _{M_*\text{--}M_h}$, was then included in the M *–Mh relation and the shears were computed again. This process was repeated for a selection of values of 
$\sigma _{M_*\text{--}M_h}$ between 0.1 and 0.5. This range of scatter is chosen to bracket the current literature values (e.g. Yang, Mo, & van den Bosch Reference Yang, Mo and van den Bosch2009, More et al. Reference More, van den Bosch, Cacciato, Mo, Yang and Li2009, Behroozi et al. Reference Behroozi, Conroy and Wechsler2010, Moster et al. Reference Moster, Somerville, Maulbetsch, van den Bosch, Macciò, Naab and Oser2010). The resulting ‘true’ shears for each dataset were then ‘observed’ with a range of measurement errors, so that the ‘observed’ shear in the tight datasets contained a shear measurement error only, while the ‘observed’ shear in the scattered datasets contained a shear measurement error, and additional scatter from the M *–Mh relation. The key questions are: (1) Can we identify the scatter from the M *–Mh relation over the scatter from measurement error? (2) How many objects do we need to do so? (3) What is the uncertainty on the measurement? To answer these questions, one first needs to establish that for a given number of galaxies, the distribution of shears arising from a tight M *–Mh relation can be distinguished from the distribution with scatter in the M *–Mh relation. The distributions of shears with and without scatter in the M *–Mh relation are shown in Figures 9–11. Figure 9 shows the observed shear as a function of the predicted shear for N = 1 000, 15 000, 50 000, and 150 000, for a shear measurement error of σγobs = 0.02 and M *–Mh relation scatter of 
$\sigma _{M_*\text{--}M_h} =0.3$ dex. The solid blue and red lines show the mean in the tight and scattered distributions respectively, and the dashed blue and red lines show ± 1σ from the mean in the tight and scattered distributions respectively. The shaded background shows the relative excess (red) or shortfall (blue) of galaxies with scatter relative to galaxies without scatter, in 2D bins of size (Δγobs × Δγpred) = (0.02 × 0.002). Figure 10 shows the observed shear as a function of the predicted shear for N = 50 000, shear measurement error of σγobs = 0.02, and scatters of 
$\sigma _{M_*\text{--}M_h} =0.1$ and 0.5 dex. Lines and shading are as in Figure 9. Figure 11 shows the observed shear for σγobs = 0.005, 0.01, 0.03, and 0.1 for N = 50 000 and 
$\sigma _{M_*\text{--}M_h}= 0.3$ dex. Lines and shading are as in Figure 9.
Distributions of observed shears as a function of predicted shears with and without scatter in the M *–Mh relation, for N = 1 000, 15 000, 50 000, and 150 000, respectively (clockwise from top left), with σγobs = 0.02 and 
$\sigma _{M_*\text{--}M_h} =0.3$ dex. The solid blue and red lines show the mean in the tight and scattered distributions respectively. The dashed blue and red lines show ± 1σ from the mean in the tight and scattered distributions respectively. The shaded background (and colourbar at the top of each plot) shows the relative excess (red) or shortfall (blue) of galaxies with scatter relative to galaxies without scatter, n binscatt − n bintight, in 2D bins of size (Δγobs × Δγpred) = (0.02 × 0.002).
Distributions of observed shears as a function of predicted shears with and without scatter in the M *–Mh relation, for 
$\sigma _{M_*\text{--}M_h} = 0.1$ dex (left plot), and 0.5 dex (right plot), with σγobs = 0.02 and N = 15 000. All lines and shaded regions are as in Figure 9.
Distributions of observed shears as a function of predicted shears with and without scatter in the M *–Mh relation, for σγobs = 0.005, 0.01, 0.03, and 0.1, with 
$\sigma _{M_*\text{--}M_h} = 0.3$ dex and N = 15 000. All lines and shaded regions are as in Figure 9.
The distributions of observed shears with and without scatter in the M *–Mh relation were compared using the Two Sample Kolmogorov–Smirnov (2SKS) test. For it to be possible for the scatter in the M *–Mh relation to be measurable from comparing the distribution of observed to predicted shears with and without scatter, we require the distributions to fail the hypothesis that they are drawn from the same distribution, under the 2SKS test. That is, we require the p-value of the test to be small. The resulting 2SKS scores and p-values (denoted S 2SKS and P 2SKS) for the combinations of n and 
$\sigma _{M_*\text{--}M_h}$ considered are presented in Table 1.
The Two Sample Kolmogorov–Smirnov test score, S 2SKS, and p-value, P 2SKS, for the scattered and tight datasets, for each combination of N, σγobs and 
$\sigma _{M_*\text{--}M_h}$ investigated. n is the number of simulated datapoints for each real galaxy, N is the total number of simulated datapoints, and σM * − Mh is the scatter in the M *–Mh relation.
The next step is to fit for the scatter in M *–Mh relation in simulated datasets. To do this, we assume that the observed shears are distributed according to both the scatter in the M *–Mh relation, and an observation error, so that the likelihood of an observed shear given the true underlying shear is given by
$$\begin{eqnarray}
L(\gamma _{\rm obs}|\gamma _{\rm true}, \sigma _{M_*-M_h}, \sigma _{\gamma _{\rm obs}}) &=& \int N(\tau - \gamma _{\rm obs},0,\sigma _{\gamma _{\rm obs}}) \nonumber\\
&&\times\, M(\tau ,\gamma _{\rm true},\sigma _{M_*-M_h}){\rm d}\tau,\quad
\end{eqnarray}$$
where N is a Gaussian distribution representing the measurement error in the DSM method, and M is a lognormal distribution representing the scatter in the M *–Mh relation. N is given by
$$\begin{equation}
N(x,\mu ,\sigma ) = \frac{1}{\sqrt{2\pi }\sigma }\exp {\left[-\frac{1}{2}\left( \frac{x-\mu }{\sigma } \right)^2\right]},
\end{equation}$$
where μ is the mean of the distribution and σ is the standard deviation. The lognormal, M, is a distribution whose logarithm is normally distributed. It is given by
$$\begin{equation}
M(x,\mu ,\sigma ) = \frac{1}{\sqrt{2\pi }\sigma x}\exp {\left[-\frac{1}{2}\left( \frac{\log (x/\mu )}{\sigma } \right)^2\right]},
\end{equation}$$
where μ and σ are the mean and standard deviation of the underlying normal distribution. The total likelihood of a given value of σM * − Mh for fixed σγobs is then given by
$$\begin{equation}
\mathcal {L}(\sigma _{M_*-M_h}) = \sum _i L(\gamma _{{\rm obs},i}|\gamma _{{\rm true},i}, \sigma _{M_*-M_h} ,\sigma _{\gamma _{\rm obs}}).
\end{equation}$$
We have fitted for the scatter in the M *–Mh relation for an assumed true scatter of σM * − Mh = 0.3 dex by maximising 
$\log [\mathcal {L}(\sigma _{M_*\text{--}M_h})]$ for a range of values of N and σγobs to investigate the behaviour of the uncertainty in the fit with these parameters. A Fisher Matrix analysis was used to estimate the standard error in the maximum likelihood for each σγobs and N. Figure 12 shows the standard error, 
$\Delta \sigma _{M_*\text{--}M_h}$, as a function of these parameters. We find that for an assumed scatter of σM * − Mh = 0.3 dex, to obtain a robust fit with a measurement error of σγobs = 0.02 a dataset of N ~ 50 000 DSM measurements is required, a result which is consistent with the results of the 2SKS test.
The standard error in the maximum value of 
$\log [\mathcal {L}(\sigma _{M_*-M_h})]$ as a function of the shear measurement error, σγobs for values of N in the range 1 000 < N < 150 000.
5.2.1. Discussion
For fixed σγobs and σM * − Mh, the p-value for N = 15 000 is too large to rule out the hypothesis that the two samples of data come from the same distribution, while the p-value for N = 50 000 is sufficient to rule out this hypothesis, and the p-value for N = 150 000 can easily do so. Similarly, for fixed N and σM * − Mh, the p-value for σγobs = 0.03 is too large to confidently rule out the hypothesis that the two samples of data come from the same distribution, while the p-value for σγobs = 0.02 is sufficient to rule out this hypothesis, and the p-value for σγobs = 0.01 can easily do so. Not surprisingly, a larger scatter in the M *–Mh relation results in a smaller p-value, since a larger scatter directly increases the difference between the datasets with and without scatter in them. From the above considerations, we conclude that in order to measure a scatter of σM * − Mh = 0.3 dex, a sample of ≳ 50 000 galaxies with a shear measurement error of σγobs ≲ 0.02 would be required. For larger values of the scatter, larger values of σγobs and smaller values of N would be sufficient.
This result is confirmed in our fits for a scatter in the M *–Mh relation of 0.3 dex for a range of values of N and σγobs, presented in Figure 12. For a true M *–Mh scatter of 0.3 dex, a shear measurement error of σγobs = 0.02, and a sample of N = 50 000 shear measurements, we recover a fitted scatter of 0.308 ± 0.02 dex.
DSM is expected to achieve measurement errors of σγobs ~ 0.02, and so shear measurement error is not seen to be a limiting factor in measuring scatter in the M *–Mh relation. While a sample of ≳ 50 000 galaxies with spatial and spectral resolution does not yet exist, there are several surveys beginning in the near to intermediate future which will provide datasets of a sufficient size to perform this experiment, for example, surveys with the Hector instrument (Lawrence et al. Reference Lawrence2012) on the Anglo-Australian Telescope (AAT), and The Hobby–Eberly Telescope Dark Energy Experiment (HETDEX; Hill et al., Reference Hill, Kodama, Yamada and Aoki2008), or surveys on the Square Kilometre Array (SKA). We conclude that while it would not be possible to utilise DSM to measure the scatter in the M *–Mh relation with existing IFU survey data, it will be possible with data from upcoming surveys.
6 CONCLUSIONS
We have made an analytical estimate of the frequency of a source being weakly lensed given a uniformly distributed population of lenses, following Mortlock & Webster (Reference Mortlock and Webster2000). We have adapted their work for the weak lensing case, in which we consider the probability of the source being lensed by at least some limiting value γlim. The results of this analysis suggest the probability of detecting weak lensing greater than a limiting value of γlim = 0.02 in a realistically observable redshift range (z ≲ 1) is non-negligible. Given this, we have created a lensing frequency algorithm which searches an input dataset for all lens–source pairs with an estimated shear greater than a limiting value of γlim. Our algorithm has been applied to a dataset extracted from the GAMA survey catalogue, and the number of objects with an estimated shear of at least γ = 0.02 in the sample was found to be ~ 393. These targets can be matched to objects in the SDSS DR10 Catalogue, and a subsample of good targets can be chosen from this selection for follow up observations.
A scatter in the M *–Mh relation results in a shift towards higher measured shears for a given population of galaxies. Given this, we have investigated the feasibility of measuring the scatter in the M *–Mh relation using shear statistics. We find that for a given shear measurement error, our ability to differentiate between a distribution of shears from a ‘tight’ M *–Mh relation, and one with scatter, is dependent on the size of the scatter, the number of objects in the sample, and the shear measurement error. For a scatter of 0.3 dex in the M *–Mh relation, we find that a sample size of ~ 50 000 galaxies would be needed to measure the scatter, for a measurement error on the shear of 0.02 (a value consistent with the shear measurement accuracy achievable with DSM). We attempt to fit for the scatter in the M *–Mh relation for a set of simulated datasets. The result of this is shown in Figure 12. For a true scatter of 0.3 dex, a shear measurement error of 0.02 and ~ 50 000 shear measurements, we recover a scatter of 0.308 ± 0.02 dex. It should be noted that the technique we have demonstrated here is based on a relatively untested algorithm, however it is adequate as an illustrative example and our results are promising for future measurements. While there are no existing IFU survey catalogues of a sufficient size to apply this technique, there are several surveys beginning in the near to intermediate future which will provide datasets of a sufficient size to perform this experiment, for example HETDEX which aims to observe ~ 106 galaxies, or surveys on the SKA, such as the ‘billion galaxy survey’ which aims to observe ~ 109 galaxies.
ACKNOWLEDGEMENTS
We thank the referee for their constructive feedback. Their input into this publication has greatly improved it. This research was conducted as part of the Australian Research Council Centre for Excellence for All-Sky Astrophysics (CAASTRO), through project number CE110001020. We acknowledge financial support from The University of Melbourne, and the Australian Astronomical Observatory (AAO). Special thanks goes to the members of the GAMA survey team for their advice, input and helpful discussions. GAMA is a joint European–Australasian project based around a spectroscopic campaign using the Anglo-Australian Telescope. The GAMA input catalogue is based on data taken from the Sloan Digital Sky Survey and the UKIRT Infrared Deep Sky Survey. Complementary imaging of the GAMA regions is being obtained by a number of independent survey programs including GALEX MIS, VST KiDS, VISTA VIKING, WISE, Herschel-ATLAS, GMRT and ASKAP providing UV to radio coverage. GAMA is funded by the STFC (UK), the ARC (Australia), the AAO, and the participating institutions. The GAMA website is http://www.gama-survey.org/. Funding for SDSS-III has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, and the US Department of Energy Office of Science. The SDSS-III website is http://www.sdss3.org/.
