Evolution of the luminositytohalo mass relation of LRGs from a combined SDSSDR10+RCS2 analysis
Key Words.:
gravitational lensing  dark matter haloesWe study the evolution of the luminositytohalo mass relation of Luminous Red Galaxies (LRGs). We select a sample of 52 000 LOWZ and CMASS LRGs from the Baryon Oscillation Spectroscopic Survey (BOSS) SDSSDR10 in the 450 deg that overlaps with imaging data from the second Redsequence Cluster Survey (RCS2), group them into bins of absolute magnitude and redshift and measure their weak lensing signals. The source redshift distribution has a median of 0.7, which allows us to study the lensing signal as a function of lens redshift. We interpret the lensing signal using a halo model, from which we obtain the halo masses as well as the normalisations of the massconcentration relations. We find that the concentration of haloes that host LRGs is consistent with dark matter only simulations once we allow for miscentering or satellites in the modelling. The slope of the luminositytohalo mass relation has a typical value of 1.4 and does not change with redshift, but we do find evidence for a change in amplitude: the average halo mass of LOWZ galaxies increases by % between and 0.22 to an average value of . If we extend the redshift range using the CMASS galaxies and assume that they are the progenitors of the LOWZ sample, we find that the average mass of LRGs increases by between and 0.2.
1 Introduction
Hierarchical models of structure formation predict that galaxies form in small dark matter haloes, which subsequently clump together and merge into larger ones (White & Rees 1978). At large scales, the evolution of structure is mainly determined by the properties of dark matter and dark energy. However, at smaller, galactic scales, baryonic physics cannot be ignored. Processes such as supernova and AGN feedback impact the relation between the observable (baryonic) properties of galaxies and their dark matter haloes. Hence by measuring these relations, we gain insight into the processes that affected them. Studying this with numerical simulations is notoriously difficult, although in recent years this field has rapidly advanced through the use of semianalytic models (e.g. Baugh 2006) and hydrodynamical simulations (e.g. Vogelsberger et al. 2014; Schaye et al. 2015). To test these simulations and guide them with further input, we need observations of the relation between the properties of galaxies and their dark matter haloes. This is also crucial for understanding the effect of baryonic physics on the dark matter power spectrum (e.g van Daalen et al. 2011; Semboloni et al. 2011), which is the main observable in weak lensing studies that aim to extract cosmological parameters, such as Euclid (Laureijs et al. 2011).
The properties of dark matter haloes around galaxies can be studied with weak gravitational lensing. As the photons emitted by distant galaxies traverse the Universe, they are deflected due to the curvature of space around intervening mass inhomogeneities in the foreground. Consequently, the observed shapes of these background galaxies slightly deform, a distortion that can be reliably measured out to projected separations of tens of Mpcs around the lenses (e.g. Mandelbaum et al. 2013). Since this completely covers the regime where the dark matter halo of any lens dominates, weak gravitational lensing offers an excellent tool to measure halo masses. The weak lensing signal of individual galaxies is too noisy to be detected, but by averaging the signal of many lenses of similar observable properties, e.g. in a certain luminosity range, we can learn about the average halo properties of such lens samples.
The relation between the properties of galaxies and their dark matter haloes has been studied before with weak lensing (e.g. Hoekstra et al. 2005; Mandelbaum et al. 2006b; Li et al. 2009; van Uitert et al. 2011; Brimioulle et al. 2013; Velander et al. 2014), but most of these studies focused on lenses at a limited redshift range. However, to study how galaxies evolve, one would like to measure how the luminositytohalo mass relation depends on lookback time. Recent imaging surveys such as the CanadaFranceHawaii Telescope Survey (CFHTLS) and the second Redsequence Cluster Survey (RCS2) contain sufficient statistical power to enable such studies. Redshift dependent constraints that are derived in a homogeneous way, as is done in this study, are particularly useful for numerical simulations, as they can potentially disentangle degeneracies among the model parameters and limit the space for finetuning to match lowredshift observations (for an example, see Figure 23 in Guo et al. 2011).
In this work, we study a particular type of galaxies: Luminous Red Galaxies (LRGs). They form an interesting subsample of the total population of galaxies, as they trace the highest density peaks in the Universe. These galaxies are thought to have formed around 2 during a relatively short and intense period of star formation, after which the formation of stars practically halted. Their luminosity evolution can therefore be approximately described by ‘passive evolution’, the evolution of a stellar population without forming new stars (e.g. Glazebrook et al. 2004; Cimatti et al. 2006; Roseboom et al. 2006; Cool et al. 2008; Banerji et al. 2010). This enables us to model the luminosity evolution, for example with stellar population synthesis models (e.g. Bruzual & Charlot 2003; Conroy et al. 2009, 2010; Maraston et al. 2009), and separate that from the halo mass evolution part. Lowlevel star formation and mergers may also contribute to the luminosity evolution of LRGs, but this is thought to mainly affect less massive LRGs (Scarlata et al. 2007; Pozzetti et al. 2010; Tojeiro & Percival 2010; Tojeiro et al. 2011; Tojeiro & Percival 2011; Tojeiro et al. 2012). How large the average impact is on the luminosity evolution, compared to the pure passive evolution scenario, is not clear. However, for massive and luminous LRGs, the luminosity evolution is thought to be well understood.
Also from an observational perspective, LRGs are advantageous to study. They are easily selected in multiband optical datasets and their redshifts can be relatively easily determined using the 4000Å break (Eisenstein et al. 2001). More than a million LRGs have been observed spectroscopically as part of the Baryon Oscillation Spectroscopic Survey (BOSS; Dawson et al. 2013), forming the LOWZ sample, which targets galaxies, and the CMASS sample, which targets galaxies. From a weak lensing perspective, the advantage of LRGs is that they are massive and therefore produce a large lensing signal that can be measured up to relatively high redshift. The overlap between the BOSS survey and the RCS2 therefore offers a perfect combined dataset to study the evolution of the luminositytohalo mass relation of LRGs.
The outline of this work is as follows. In Section 2 we describe the data that we use in this work, how we compute the luminosities and how we perform the lensing analysis. We interpret the lensing measurements with the halo model, which we describe in Section 3. The evolution of the luminositytohalo mass relation is presented and discussed in Section 4. The massconcentration relation is discussed in Section 5. We conclude in Section 6. Unless stated otherwise, we assume a WMAP7 cosmology (Komatsu et al. 2011) with , , , , and km s Mpc with the Hubble constant; all distances are quoted in physical (rather than comoving) units; and all apparent magnitudes have been corrected using the dust maps from Schlegel et al. (1998).
2 Data analysis
In this work we make use of data from the tenth data release (DR10; Ahn et al. 2014) from the Sloan Digital Sky Survey (SDSS; York et al. 2000) and from the second Redsequence Cluster Survey (RCS2; Gilbank et al. 2011). As in van Uitert et al. (2011, 2013); Cacciato et al. (2014), we take advantage that the SDSS contains more ancillary data on galaxies than is available in the RCS2 due to its photometry in five optical band and its spectroscopy for over a million of galaxies. However, the RCS2 imaging is 2 magnitudes deeper and achieved a median seeing of approximately , compared to for SDSS, making the RCS2 better suited for a weak lensing analysis of lenses at higher redshifts. The total overlap between the RCS2 and the DR10 amounts to roughly 450 square degrees.
A first combined analysis of the overlap between the ninth data release of SDSS (DR9; Ahn et al. 2012) and the RCS2 was presented in Cacciato et al. (2014), where the lensing signal of the DR9 galaxies with spectroscopy was studied using RCS2 galaxies as sources. In that work we did not study the redshift evolution of the lensing signal, as, in contrary to the current work, we studied a mixed sample of early and latetype galaxies, whose combined luminosity evolution was not well understood. Also, in DR9, the number of (highredshift) BOSS spectra was roughly half of that in DR10.
2.1 Lens sample
We use a subset of the total sample of overlapping DR10 galaxies with spectroscopy as our lenses, i.e. only the LRGs. We select all galaxies that have been targeted as part of BOSS. These are selected from the SDSS catalogues by requiring

BOSS_TARGET1 && 2

SPECPRIMARY == 1

ZWARNING_NOQSO == 0

TILEID 10324
for the LOWZ sample, and

BOSS_TARGET1 && 2

SPECPRIMARY == 1

ZWARNING_NOQSO == 0

(CHUNK != “boss1”) && (CHUNK != “boss2”)

for the CMASS (highz) sample. Additionally, we select all objects with reliable spectroscopy from the SDSS catalogues that satisfy the BOSS LOWZ target selection cuts:




SCIENCEPRIMARY==1

ZWARNING_NOQSO == 0

where , and indicate model magnitudes and cmodel magnitudes. Note that we replaced the BOSS selection criterion with to ensure that we have no stars. Finally, we also selected all objects that satisfied the CMASS selection cuts, but we found that all objects were already targeted and labeled as being BOSS galaxies, and it therefore did not increase the lens sample.
Even though the LOWZ and CMASS samples mainly consist of LRGs, the populations differ due to the different colour and magnitude selection cuts. Tojeiro et al. (2012) study which fraction of the CMASS LRGs are progenitors of the LOWZ sample, and find that this strongly depends on absolute magnitude, with the highest fractions found for the most luminous objects. A second but weaker trend is found with restframe colour. Therefore, we choose to analyse the LOWZ and the CMASS samples separately. We do investigate what we can conclude about the evolution of the luminositytohalo mass relation of LRGs if we consider the CMASS sample as progenitors of the LOWZ LRGs.
Luminosities
In order to study how the average halo mass of LRGs evolves, we want to compare LRGs at low redshifts to their predecessors at higher redshifts. Hence we need to obtain the luminosities of our LRGs, corrected for the redshift of their spectra through the passbands (i.e. the correction). We compute the correction using the KCORRECT v4_2 code (Blanton et al. 2003; Blanton & Roweis 2007), where we use the , , , and model magnitudes and the spectroscopic redshift as input. Furthermore, we correct for the intrinsic evolution of the luminosities (the correction), accounting for the difference between the observerframe absolute magnitude of a galaxy with and without an evolving spectrum.
The luminosities of LRGs are thought to evolve passively, which can be modeled using a stellar population synthesis code. We make use of one of the publicly available codes, GALAXEV (Bruzual & Charlot 2003), in the default configuration, i.e. adopting a Chabrier (2003) IMF and using the Padova 1994 tracks for the stellar evolution. We compute a range of instantaneousburst models, where we vary the formation time and the metallicity. In Figure 1, we show the evolution of the and colours of these models, together with the observed colours of the LRGs. The set of models that describe the data best are those that assume a metallicity of (). However, at , the observed colours are a bit too red, and at the colours are somewhat too red. Maraston et al. (2009) improved the modelling by including a very low metallicity component to the model that consisted of 3% in mass, and by using the Pickles (1998) empirical spectral library instead of the theoretical one. However, below a redshift of 0.5 the evolution correction is fairly insensitive to the details of the modelling (see Figure 2), while at higher redshifts it is not clear whether the changes from Maraston et al. (2009) improve the match due to the low number of objects at this redshift range used in that work. As we will discuss below, our results do not critically depend on the choice of the model, hence we do not deem it necessary to include the improvements from Maraston et al. (2009).
In Figure 2 we show the correction that these GALAXEV models predict, together with the correction for the LRGs that have been computed using KCORRECT. We find that at , the tracks agree well, but at higher redshifts the correction values from KCORRECT are somewhat lower than the GALAXEV model. In fact, the agreement at with the metallicity models is remarkable, but the validity of these models for our LRGs at low redshift is excluded based on the colour evolution in Figure 1. However, at the LRGs show an increasing scatter in their colours and become more compatible with the models.
We show the luminosity evolution of some of the GALAXEV models in Figure 2. We only show the models with and , as the models with different metallicities are excluded based on their colour evolution and corrections. Also, we only show models with a formation time of 10, 11 and 12 Gyrs, as most previous works on the luminosity evolution of LRGs have adopted a formation time in this range (e.g Wake et al. 2006; Maraston et al. 2009; Banerji et al. 2010; Carson & Nichol 2010; Liu et al. 2012). For our nominal luminosity evolution correction, we adopt the model that formed 11 Gyrs ago (at redshift 0). As none of the models exactly captures the trends in Figure 1 and 2, the evolution correction we use may have a small bias. However, we have also tried different evolution corrections, with corresponding models that broadly cover the observed colour evolution and correction values. We detail on this test in Section 4.1; the main result is that our results do not change significantly. This suggests that the systematic bias in the luminosities caused by an incorrect evolution correction is likely insignificant for this work.
LRGs have formed over a certain range of time and with some range of metallicities. Hence their actual luminosity evolution corrections may have some scatter compared to our nominal correction, as we found that the luminosity evolution correction is increasingly sensitive with redshift to these parameters. If this scatter is random with respect to our nominal correction, this causes an Eddington bias, as lenses are preferentially scattered to where there are fewer of them. In Appendix A, we estimate the impact this may have on our masses estimates. We find that it is significantly smaller than our statistical errors and can be safely ignored.
In Figure 3, we show the distribution of absolute magnitudes after including the correction and the ()correction. In the range the distribution of corrected absolute magnitudes is fairly flat. At redshifts we have a tail of fainter objects in our catalogues, which are likely different types of galaxies. Therefore, we exclude them from this analysis. At higher redshifts, we start loosing fainter objects due to incompleteness. Since we determine both the mean luminosity and the mean halo mass for a given lens sample, this should not bias the overall masstoluminosity relation.
2.2 Lensing measurement
The shapes of the background galaxies are measured on images from the RCS2. Details on the data reduction and the shape measurement process can be found in van Uitert et al. (2011), and some important improvements of our lensing analysis were discussed in Cacciato et al. (2014). It suffices to say that we measure the shapes of the galaxies with the KSB method (Kaiser et al. 1995; Luppino & Kaiser 1997; Hoekstra et al. 1998), using the implementation described by Hoekstra et al. (1998, 2000). This method was tested on a range of simulations as part of the Shear Testing Programme (STEP) 1 and 2 (the ‘HH’ method in Heymans et al. 2006 and Massey et al. 2007 respectively) where it was found to have a multiplicative bias of a few per cent at most and a negligible additive bias. Recently, Hoekstra et al. (2015) found that these results were driven by the overly simplistic nature of the STEP simulations; for more realistic simulations, KSB suffers from noise bias (Kacprzak et al. 2012; Melchior & Viola 2012; Refregier et al. 2012) as any other shape measurement method that is currently in use. We calibrate our KSB implementation on realistic image simulations generated with GalSim
(1) 
with FWHM the size of stars in an image. We use this to correct the shear measured in each RCS2 image. We do not need to correct for residual additive bias as that generally averages out in galaxygalaxy lensing due to symmetry in lenssource pair orientations. The (multiplicative) noise bias correction increases the average lensing signal by 1015%. Note that the correction is not very sensitive to the adopted width of the intrinsic ellipticity distribution, but that it is critically important to simulate galaxies up to 1.5 magnitudes deeper than the nominal magnitude limit of the survey (see Hoekstra et al. 2015).
The lensing signal is extracted by azimuthally averaging the tangential projections of the ellipticities of the source galaxies in concentric radial bins, i.e. by measuring the tangential shear as a function of projected separation:
(2) 
where is the difference between the mean projected surface density inside radius and the projected surface density at , and is the critical surface density:
(3) 
with , , and the angular diameter distance to the lens, the source, and between the lens and the source, respectively. All galaxies with an apparent magnitude of and a welldefined shape measurement are selected as sources.
We measure the lensing signal from the BOSS lenses in each 11 deg RCS2 pointing, including the sources from the neighbouring pointings (if present). We bootstrap over these patches to obtain the covariance matrix, which accounts for intrinsic shape noise, measurement noise, as well as the contribution from largescale structures. The offdiagonal elements are consistent with zero on the radial range of interest (Mpc), hence we only use the inverse of the diagonal as the errors on the measurement when fitting the models to the data.
As these patches overlap, the contribution from largescale structures might be somewhat underestimated at large scales; therefore, as a test, we also perform the lensing measurements on 22 deg nonoverlapping patches and use that in the bootstrap resampling. We find that for both methods, the signal and the covariance matrix barely change in the radial range that we use in this work. Since the total area decreases if we limit our analysis to 22 deg nonoverlapping patches only, and since it makes no difference to the signal and its error, we decide to use all 11 deg RCS2 pointings plus neighbours as basis for the bootstrapping.
In order to compute we need the distances to the lenses and sources. We compute for each lens separately using its spectroscopic redshift from SDSS. The lensing efficiencies are determined by averaging over the source redshift distribution, which is obtained by applying the same band selection to the photometric redshift catalogues of the COSMOS field from Ilbert et al. (2013). The procedure is described in more detail in Appendix C of Cacciato et al. (2014). Note that we previously used the photometric redshift catalogues from Ilbert et al. (2009) as the former was not yet publicly available. Comparing the average lensing efficiencies from the two catalogues, we find them to agree for low lens redshifts, but to be increasingly different at higher redshifts (up to 15% at z=0.7). At increasingly high lens redshift the lensing efficiencies are more sensitive to the form of the adopted source redshift distribution, which is somewhat different for the two catalogues. We discuss the robustness of the derived lensing efficiencies in more detail in Appendix B.
A fraction of the sources is physically associated to the lens galaxies, representing an overdensity of sources galaxies that are not lensed. We cannot remove them since we lack redshifts for our sources. Such a contamination in the source catalogue dilutes the lensing signal. This can be corrected for by measuring the excess source number density relative to the background as a function of projected separation, and boosting the lensing signal with this factor (e.g. Mandelbaum et al. 2006b; van Uitert et al. 2011). We follow the same procedure.
This boost correction itself is biased low as the galaxies associated to the lens, and the lens itself, block light from the background sky, suppressing the source number density. The effect is described in Simet & Mandelbaum (2014). As discussed in that work, a correction for this bias is obtained by multiplying the boost correction with a factor , where is the fraction of the sky that is obscured by the foreground galaxies. We compute this by using the ISOAREA_IMAGE keyword in SExtractor, which stores how many pixels a galaxy spans on the sky. is taken to be the sum of these values of all galaxies whose centroids fall inside a radial bin, divided by the total number of pixels in that bin (accounting for the effect of the survey masks and geometry). Before doing the correction, we subtract the average skybackground of from the one observed around the lenses, as we are only interested in the additional obscuration. We find that the correction is at most 5%, in the radial bins closest to the most luminous and lowredshift lenses. The correction decreases at larger radii, as well as for fainter, higher redshift LRGs, as expected.
The robustness of the lensing signal has been addressed in Appendix B of Cacciato et al. (2014). There, we show that the cross shear signal of our lens sample is consistent with zero. Also, the random shear signal, which is used to correct for the effect of residual systematics in the shape measurement catalogues, is smaller than the real signal for all the projected separations we use in this work. However, an overall multiplicative bias could still be present, either through an incorrect determination of the noise bias correction, or through the use of incorrect lensing efficiencies. In Appendix B we perform an internal consistency check of our measurement pipeline, which provides strong evidence that such a bias is unlikely to be significant.
3 Halo model
In this section we describe the model that we employ to provide a physical interpretation of our measurements. The halo model provides a useful framework to describe the stacked weak lensing signal around galaxies (see e.g. Mandelbaum et al. 2006b; Cacciato et al. 2009, 2014; Miyatake et al. 2013). It is based on a statistical description of dark matter properties, such as their average density profile, their abundance, and their large scale bias, complemented with a statistical description of the way galaxies of a given luminosity populate dark matter haloes of different masses (also known as halo occupation statistics). In its fundamental assumptions, the model is similar to the one presented in Seljak (2000), Cooray & Sheth (2002) and Cacciato et al. (2009).
Galaxygalaxy lensing probes the average matter distribution projected along the lineofsight at a given projected physical separation, , for a set of lenses. The quantity of interest is the excess surface mass density profile, , which is determined from the projected surface mass density, . Since we measure the average signal of many lenses, the projected matter density can be expressed in terms of the galaxydark matter crosscorrelation, :
(4) 
where the integral is along the line of sight, is the comoving distance from the observer, the comoving distance to the source and is the mean matter density at the redshift of the lens. The threedimensional comoving distance is related to through , with the comoving distance to the lens and the angular separation between lens and source (see Figure 1 in Cacciato et al. (2009)). Note that the galaxydark matter crosscorrelation is evaluated at the average redshift of the lens galaxies.
Under the assumption that each galaxy resides in a dark matter halo, can be computed using a statistical description of how galaxies are distributed over dark matter haloes of different masses (see e.g. van den Bosch et al. 2013). Specifically, it is fairly straightforward to obtain the twopoint correlation function, , by Fourier transforming the galaxydark matter powerspectrum, , i.e.
(5) 
with the wavenumber. The quantity can be expressed as a sum of a term that describes the small scales (onehalo, 1h) and one that describes the large scales (twohalo, 2h), each of which can be further subdivided based upon the type of the galaxies (central or satellite) that contribute to the power spectrum. This reads
(6) 
The terms in Equation 6 can be written in compact form as
(7) 
(8)  
where ‘x’ and ‘y’ are either ‘c’ (for central), ‘s’ (for satellite), or ‘m’ (for matter), describes the power spectrum of haloes of mass and , and it contains the large scale bias of haloes from Tinker et al. (2010) (but see van den Bosch et al. (2013) for a more sophisticated modelling of this term). is the halo mass function of Tinker et al. (2010). Furthermore, we have defined
(9) 
and
(10) 
where
(11) 
and
(12) 
is the parameter that describes the probability that the ‘central’ galaxy does not reside at the centre of the dark matter halo, whereas quantifies the amount of offcentering in terms of the halo scale radius, (see e.g. Skibba et al. 2011; More et al. 2014). In our fiducial model, we set , but we explore the impact of this assumption in Section 4. The functions and represent the average number of central and satellite galaxies in a halo of mass , defined as:
(13) 
(14) 
where
(15) 
We use a flat, noninformative prior for and , set (see Equation 18) and . Since LRGs are thought to be predominantly central galaxies (see e.g. Wake et al. 2008; Zheng et al. 2009; Parejko et al. 2013), we set to zero and only fit for and in our nominal runs. We test the impact of this assumption on the derived quantities in the result sections by additionally fitting for . We have tested that the result is fairly insensitive to the details of the modelling of .
is the number density of galaxies at redshift :
(16)  
The last equality is exact in the case of LRGs being only central galaxies. is the Fourier transform of the normalized density distribution of dark matter within a halo of mass , for which we assume a NavarroFrenkWhite (NFW) profile (Navarro et al. 1996) and a massconcentration relation from Duffy et al. (2008):
(17) 
with , , , and . Note that is a free parameter that allows the normalisation of this relation to vary. Specifically, we apply a noninformative flat prior on this parameter.
The average halo mass in a given luminosity bin, which is what we shall refer to as ‘effective’ halo mass in what follows, can then be computing taking into account the weight of the halo mass function:
(18) 
where is the mean redshift of the lens galaxies in a given luminosity bin. The distinction between the mass associated with the mean of the lognormal distribution and the halo mass inferred accounting for the mass function is of relevance because LRGs populate fairly massive haloes for which the mass function is steep (see e.g. Figure 7 in Leauthaud et al. (2015)).
At small scales one expects the baryonic mass of LRGs to contribute to the lensing signal. The smallest scale used in this study is 50 kpc, which is much larger than the typical extent of the baryonic content of a galaxy. Therefore, it is adequate to model the lensing signal of the LRGs itself as a point source of mass . This reads
(19) 
We use the value of the average stellar masses, , for the galaxies in the luminosity bins under investigation here. The stellar masses are obtained by matching our lens catalogue to the MPAJHU stellar mass catalogue
To summarize, our fiducial model for the lensing signal is the sum of three terms: one describing the lensing due to the baryonic mass; the second is responsible for the small (subMpc) scale signal mostly due to the dark matter density profile of haloes hosting central LRGs; and the last describes the large (a few Mpc) scale signal due to the clustering of dark matter haloes around LRGs. This reads:
(20) 
We simultaneously fit the halo model to the three luminosity bins and do this for each redshift slice separately. We have five free parameters in each fit: the three mean masses of the luminosity bins, the scatter and the normalisation of the massconcentration relation. The fit is performed using a Markov Chain Monte Carlo Method (MCMC). Details of its implementation can be found in Appendix C.
We fit the model to the measurements on scales between 0.05 and 2 Mpc. At these scales, both the measured lensing signal and the halo model predictions are fairly robust. At larger scales, the lensing signal becomes smaller and more susceptible to residual systematics. At scales smaller than 0.05 Mpc, lens light may bias the shape measurements. For the halo model, the overlap between the 1halo and 2halo term is notoriously hard to model because of, amongst others, halo exclusion and nonlinear biasing. This mainly impacts the few Mpc regime. Most of the information about the halo masses and concentrations is contained in the lensing signal within the virial radius, so we do not loose much statistical precision by limiting ourselves to these scales.
4 Luminositytohalo mass relation
To study how the luminositytohalo mass relation of LRGs evolves with redshift, we divide our sample in bins of ()corrected absolute magnitude and redshift as detailed in Table 1 and Figure 3. For , we only select lenses from the LOWZ sample, at higher redshifts we exclusively select CMASS galaxies. The average log stellar masses for the consecutive luminosity bins are 11.2, 11.5 and 11.7 []. For each bin we measure the average lensing signal, which is shown in Figure 4, together with the bestfit halo models and the model uncertainties (computed as detailed in Appendix C). We find values of 1.8, 1.6, 1.2 and 1.0, going from the lowest to the highest redshift slice. Hence the fits of the CMASS samples are good, but for the LOWZ samples the values are somewhat large, suggesting that either our error bars are underestimated, or that the model that we fit to the data is overly simplistic.
The errors on the lensing measurements account for intrinsic shape noise, measurement noise and the impact of largescale structures. We have, however, ignored some small sources of error, as their amplitude is much smaller than the statistical errors on the lensing signal: the error on the boost correction, which is typically a few percent at small scales; the error on the obscuration correction, which is even smaller; the error on determining the lensing efficiency and the error on the multiplicative bias calibration, whose magnitudes are unknown but are likely of the order of a few percent. Combined, they might increase the errors by as much as 10%, although the exact number is difficult to estimate reliably. If we increase our error bars by this amount, we would get values of 1.5, 1.4, 1.0 and 0.8, respectively. The fact, however, that we find reasonable values for the CMASS sample, but not for the LOWZ sample, suggests that a systematic underestimate of our errors is unlikely to be the dominant cause.
Even though a visual inspection of the covariance matrix lead us to believe that it is diagonal on scales 2 Mpc, there could be lowlevel offdiagonal terms present that, if included, would lower the values. This potentially affects the LOWZ results more, as the measurements have a higher signaltonoise ratio and the covariance matrix is less noisy. To test this, we recompute the values using the full covariance matrix which we obtained from bootstrapping (see Section 2.2) for the bestfit models. Note that we only include the covariance between radial bins of a lens sample, but not the covariance between the radial bins of the different luminosity samples. If present, they would lower the values even more. We find that of the first redshift slice reduces to 1.6, while it does not change for the other three redshift slices. Hence the effect is small and does not fully explain the high values.
Figure 4 shows that the signaltonoise ratio of the lensing measurements of the LOWZ samples is very high and would allow for a more sophisticated modelling. When we include a satellite term or a miscentering term in the halo model, however, the values do not improve, as the lensing signal alone cannot constrain the miscentering distribution parameters very well, and the expected number of satellites is low. Allowing for even more freedom in the fit might lead to overfitting of the CMASS results. Using different halo models for the different samples, or splitting the LOWZ sample up in more luminosity bins, reduces the homogeneity of the analysis, which is one of the key advantages of our work. Hence we choose to stick to the settings described above. In Section 4.1 we perform a sensitivity analysis and find that our results do not critically depend on various choices in the analysis, suggesting that the quantities we derive from the fits are robust.
M  N  

(1)  (2)  (3)  (4)  (5)  (6)  (7)  
0.150.29 (LOWZ)  
L1z1  [21.8,21.2]  2969  0.219  0.65  
L2z1  [22.4,21.8]  2606  0.226  0.99  1.8  
L3z1  [22.8,22.4]  300  0.234  1.58  
0.290.43 (LOWZ)  
L1z2  [21.8,21.2]  3771  0.351  0.66  
L2z2  [22.4,21.8]  4502  0.364  1.00  1.6  
L3z2  [22.8,22.4]  721  0.368  1.58  
0.430.55 (CMASS)  
L1z3  [21.8,21.2]  8530  0.499  0.61  
L2z3  [22.4,21.8]  4213  0.503  0.99  1.2  
L3z3  [22.8,22.4]  587  0.500  1.59  
0.550.70 (CMASS)  
L1z4  [21.8,21.2]  5256  0.596  0.64  
L2z4  [22.4,21.8]  5161  0.611  1.00  1.0  
L3z4  [22.8,22.4]  969  0.616  1.60 
In the halo model, we fit the mean and the scatter of the lognormal distribution that describes . The log of the mean has typical values of 14.5, 15 and 15.5 for the three luminosity bins, while the scatter ranges between 0.7 and 0.8. Neither evolves with redshift. Note that this is the scatter in the log of the halo mass and not in luminosity. The latter was fit in Cacciato et al. (2014) where it was found to have a value of , obtained by fitting the halo model to the lensing signal of all galaxies in the DR9 that overlap with RCS2. The scatter in halo mass is much larger, because the luminosityhalo mass relation flattens at higher luminosities; a small scatter in luminosity corresponds to a large one in halo mass (see e.g. Figure 3 and the discussion in More et al. (2009)).
The quantity of interest that we can compare to other works is the ‘effective’ halo mass, which is given in Table 1. We plot it as a function of luminosity in Figure 5. We find that the masses increase with luminosity and decrease with redshift. To quantify this, we parametrize the luminositytohalo mass relation by , using a pivot luminosity of . The bestfit powerlaw fits are shown in the same figure; the confidence contours of the fitted amplitude and slope are shown in Figure 6. We have listed the fit parameters in Table 2.
The slope of the luminositytohalo mass relation does not change significantly for our different redshift samples and has a typical value of 1.4. The amplitude, however, is about 4 higher for our lowest redshift slice compared to the highest one. On average, the masses of LOWZ galaxies increase by % between redshift 0.36 and 0.22; the masses of CMASS galaxies increase by % from redshift 0.6 to 0.5. If we assume that CMASS galaxies evolve into LOWZ galaxies and combine the results, we find an average increase of % (stat. errors) in at from 0.6 to 0.2. Fixing the slope to its average value of 1.4 and only fitting the amplitude changes this number to %.
Tojeiro et al. (2012) found that at brighter absolute magnitudes, a larger fraction of CMASS galaxies are the progenitors of LOWZ galaxies. If we discard the lowest luminosity bin, we find that, for a pivot luminosity of , the average halo mass increases with %. Considering only the brightest luminosity bin, the average halo mass increases with %.
4.1 Sensitivity analysis
To study how sensitive our results are to the adopted luminosity evolution correction, we perform the following test. We multiply our nominal correction with a factor, such that the resulting luminosity evolution curves roughly cover the range of reasonable models that are shown in the lower panel of Figure 2. The factors we choose are and , respectively. Next, we recompute the luminosities, repeat the lensing measurements (using the same cuts) and the halo model fits, and compare the resulting bestfit halo masses. We find that the bestfit halo masses of the individual luminosity bins do not significantly shift compared to the nominal results. We fit the luminositytohalo mass relation and list the parameters in Table 2. As is shown there, we find that the bestfit slopes are very similar. The bestfit amplitude shifts with 1 at most compared to our nominal results. For the modification factor, the resulting increase in halo mass of LRGs from to is %; for the modification factor, it is %. The halo masses and corresponding growth change somewhat because the lens selection shifts systematically, such that we analyse lens samples that to some extent are intrinsically different. Importantly, however, our results do not critically depend on the choice of the luminosity evolution correction. For future work that is expected to have an improved statistical precision, more detailed knowledge of this correction will be required.
To test how sensitive our results are to the assumption that all LRGs are located at the centre of their dark matter haloes, we do two halo model runs where we allow for more flexibility. First, we assume that a fraction of the LRGs is miscentered, following Equation 11. We use and as additional free parameters with a flat uninformative prior in the range [0,1]. The allowed miscentering distribution ranges from the lenses being all correctly centered ( = =0) to all being miscentered and located at the halo scale radius ( = =1). We find , , and for the low and high redshift bins of LOWZ and CMASS, respectively, with corresponding miscentering radii of , , and . The resulting powerlaw parameters are listed in Table 2. We find that the powerlaw parameters do not change significantly. The total increase in halo mass of LRGs corresponds to %, consistent with our nominal result of %. These constraints on the miscentering distribution are in broad agreement with previous galaxygalaxy lensing and clustering results of CMASS galaxies. For instance, Miyatake et al. (2013) report and , whereas More et al. (2014) find and .
As a related test, we study how our results change when we assume that a fraction of LRGs are satellites. We use a simple model with only one free parameter, i.e. . We follow this approach because we are not interested to determine the satellite HOD (as lensing alone is not very sensitive to that), but because we want to get an estimate how much our results could be affected by ignoring the contribution of satellites. The typical satellite fraction for LRGs is 10% and decreases for more massive LRGs (White et al. 2011; Parejko et al. 2013; More et al. 2014). We therefore put a prior on such that the resulting satellite fractions from our model are between 5% and 15%. The resulting luminositytohalo mass relation parameters are listed in Table 2. The bestfit slopes are consistent with our nominal run, but the amplitude decreases by 12. The total increase of halo mass is % over the full redshift range of our LRG sample, consistent with our nominal result. The normalisation decreases because the satellites are associated to more massive haloes, with correspondingly larger lensing signals. This lowers the required contribution to the total signal from central LRGs, and hence their mass. This does not happen in the miscentering run, where the lensing signal is merely smeared out, but the integrated signal (and hence the mass) stays the same.
nominal run  
0.150.29  1.8  
0.290.43  1.6  
0.430.55  1.2  
0.550.70  1.0  
correction  
0.150.29  1.9  
0.290.43  1.3  
0.430.55  1.0  
0.550.70  1.0  
correction  
0.150.29  1.9  
0.290.43  1.4  
0.430.55  1.3  
0.550.70  1.0  
miscentering run  
0.150.29  1.9  
0.290.43  1.8  
0.430.55  1.3  
0.550.70  1.1  
satellite fraction run  
0.150.29  1.9  
0.290.43  1.7  
0.430.55  1.2  
0.550.70  1.0 
4.2 Comparison to previous work
Several previous works have provided mass estimates of LRGs, using gravitational lensing, clustering, abundance matching, or a combination of these. The selection of the samples, the models fit to the data and the definitions of mass, generally differ between these studies, limiting the level of detail with which we can perform a comparison.
Lensing results
Mandelbaum et al. (2006a) measured the masses for a sample of over LRGs with spectroscopic redshifts from SDSSI/II using weak lensing. Bluer and fainter LRGs are discarded using colourmagnitude cuts, as well as LRGs that are likely satellites of larger systems, hence the selection is not identical to ours. The resulting sample is split into a faint and bright part using a cut at , and the mean luminosity of the two samples is and , respectively. The luminosities are computed in a similar manner as in our work. All LRGs are selected in the redshift range and have a mean effective redshift of 0.24. Masses are estimated using NFW fits plus a baryonic component, where the fitting range is restricted to small scales were the 2halo term can be neglected. The masses are defined as the enclosed mass in a sphere where the density is 180 times the mean background density, , instead of 200 times the mean density which we use; the difference between these definitions is only a few percent. The faint LRGs are found to reside in haloes of masses and the bright ones in haloes of masses . The measurements are shown in Figure 5. We find that the masses of our lowredshift slices are substantially larger than the ones from Mandelbaum et al. (2006a). The discrepancy may be caused by the scatter between luminosity and halo mass, which is not included in Mandelbaum et al. (2006a). Allowing for a nonzero scatter results in larger halo masses.
In Miyatake et al. (2013), the lensing signal of 4,807 CMASS galaxies with is measured in the overlap with CFHTLS using the publicly available CFHTLenS catalogues (Heymans et al. 2012). The lensing signal is fit together with the projected clustering signal using a halo model that is similar to the one we have adopted here. Halo masses are defined with respect to 200 times the background density as we do. The average halo mass of CMASS galaxies is found to be . To compare it with our results, we compute the average luminosity of CMASS galaxies with in our catalogue. We find a value of and assume that this number is representative for the average luminosity of the lenses in that work. This estimate is actually a bit too low, as Miyatake et al. (2013) also apply a cut on stellar mass to remove the least massive (hence faintest) objects, which we cannot mimic. However, as an indication, we find that if we remove the faintest 10% of our CMASS sample, the average luminosity only increases to , hence it is unlikely that the average luminosity is very far off. We compare the results in Figure 5 and find that our CMASS masses are somewhat lower but consistent.
Several papers have studied how the average mass of galaxies changes as a function of redshift (e.g. van Uitert et al. 2011; Choi et al. 2012; Leauthaud et al. 2012; Tinker et al. 2013; Hudson et al. 2015), using weak lensing measurements. Since these works do not target LRGs specifically, and since the modelling of the signal differs from our approach, we cannot compare the results in detail. However, both Tinker et al. (2013) and Hudson et al. (2015) include red, massive galaxies in their work, hence we can at least compare the recovered trends.
Tinker et al. (2013) use measurements of weak lensing, clustering and the stellar mass function of galaxies in COSMOS to constrain the stellartohalo mass relation. Halo masses are defined like ours. The relations they report predict the average as a function of halo mass, instead of the average halo mass at a given stellar mass, which is what we measure; these relations are not the same due to intrinsic stellar mass scatter, which is illustrated in their Figure 7. From the righthand panel of that figure we observe that for passive galaxies, at the average stellar masses of our samples, the mean halo mass is roughly 0.5 dex lower than what we find. There are many differences between the analyses that could contribute to this difference, such as systematic offsets between stellar mass estimates, the selection of the samples, the definition of mass and the modelling of the signal. Nonetheless, at a stellar mass of 11.4 (typical for LRGs), the average halo mass increases from to between redshifts of 0.88 and 0.36, an increase of almost 100%, similar to the average growth in halo mass that we find for our LRG sample from redshift 0.62 to 0.21.
Hudson et al. (2015) use the shape measurements from CFHTLenS to measure the lensing signal for blue and red galaxies in three redshift slices. The lenses are binned in luminosity rather than stellar mass to avoid an Eddington bias due to the larger observational errors on stellar mass compared to luminosity. The stellar mass is then determined using the mean stellarmasstoluminosity ratio. For the highest luminosity bin of red lenses, which has a mean stellar mass of , the average halo mass increases from to from redshift 0.67 to 0.29. The masses are defined with respect to instead of the mean density, resulting in masses that are 3040% smaller compared to ours. Furthermore, the intrinsic scatter between luminosity/stellar mass and halo mass is not accounted for in their modelling, which also leads to smaller masses. Finally, the selection of the lens samples differs. However, the 58% increase in average halo mass is comparable to what we find.
Clustering results
The clustering of LRGs is wellstudied in the literature and has been used to derive halo masses (e.g. Blake et al. 2008; Wake et al. 2008; Zheng et al. 2009; Sawangwit et al. 2011; Nikoloudakis et al. 2013; Parejko et al. 2013; Guo et al. 2014). Parejko et al. (2013) measure the clustering of galaxies with from the LOWZ sample and fit it with a halo model. The probability distribution of halo masses, as shown in their Figure 9, has a mean of . We plot it in Figure 5 and find that it is larger than our LOWZ measurements. Although not specified, we assume that the mass is defined as , as is mentioned in a companion paper (White et al. 2011), which is comparable to our definition. Their halo mass distribution is fairly broad, however, and our constraints may well fall inside their 68% confidence region.
Guo et al. (2014) measure the clustering of CMASS galaxies divided in three band magnitude selected samples, which we cannot directly compare to our measurements. They find that going from their faintest to their brightest bins, the peak host halo mass increases from to , which is quite comparable to our results. Their masses are defined with respect to the mean density, like ours. However, we measure an ‘effective’ mass and not the peak halo mass, and it is unclear how much these definitions differ.
Zheng et al. (2009) fit the clustering signal of SDSSI/II LRGs using an HOD approach. Their LRGs are divided in a faint and bright sample using their magnitude, hence we cannot directly compare. The distribution of halo masses (defined like ours) of these samples peaks at and , respectively, similar to the values we find for faint and bright LOWZ samples, but note, again, that we measure the effective mass and not the peak halo mass. The scaling of luminosity with host halo mass, which is given by , is consistent with our results.
Wake et al. (2008) measure the evolution of the clustering signal of galaxies from SDSS and the 2dFSDSS LRG and QSO Survey (2SLAQ; Cannon et al. 2006). They match the selections using colour and magnitude cuts, which complicates the comparison with our results. However, they find that the effective halo masses increase with 50% from to , consistent with our findings. Sawangwit et al. (2011) study three separate LRGs samples in SDSS, with a mean redshift of 0.35, 0.55 and 0.75. After applying additional selection criteria such that the space density of LRGs is similar to the SDSSI/II LRG sample, the effective halo mass, defined as the virial mass, is found to be , and , respectively. We assume that the average luminosity of these three samples is similar to that of the LOWZ sample and show their measurements in Figure 5. We find that their results agree fairly well with ours. Blake et al. (2008) and Nikoloudakis et al. (2013) study the clustering signal of different samples of LRGs, which is difficult to compare with our results. The typical halo masses of LRGs are found to be a few times , hence broadly in agreement.
Abundance matching results
Masaki et al. (2013) apply an abundance matching technique to Nbody simulations in order to construct mock LRG samples whose number density matches that of LRGs in SDSS. From these samples a mock lensing signal is constructed, whose shape agrees well with, but whose amplitude is 20% larger than, the measured lensing signal of SDSS LRGs presented in Mandelbaum et al. (2013). They find that their LRGs reside in haloes with a mean virial mass of . This definition of mass is 10% smaller than for masses and concentrations that are typical for LRGs at 0.3. We assume the average luminosity equals that of SDSSI/II LRGs and show the measurement in Figure 5. We find that this mass estimate is in fair agreement with our results.
4.3 Interpretation
We find that of LRGs increases by approximately 80% from redshift 0.6 to 0.2. This is not only because of dark matter accretion. Part of the growth can be attributed to the socalled pseudoevolution (Diemer et al. 2013). is defined with respect to a reference density, which is redshift dependent. Even if a halo is static, i.e. it does not accrete anything, increases towards lower redshift. The halo mass function that we use to determine uses masses that are defined with respect to . Hence the evolution of the halo mass function is a mix of physical evolution and pseudoevolution. Therefore, the evolution of is also a mix of the two.
Since the halo masses of our LRGs span a relatively narrow range, we can estimate the contribution from pseudoevolution to our observed increase in mass using the results from Diemer et al. (2013). For halo masses that are typical for our LRGs, increases by approximately 33% from to due to pseudoevolution, as is illustrated in Figure 7. If we fit the amplitude of the pseudoevolution curve to our measurements, we find that , providing weak support that the slope is steeper and requires additional dark matter accretion.
It is interesting to compare our growth rates to those obtained from simulations. Ideally, we would like to compare to hydrodynamical simulations, but those are not available for the mass range we are interested in. However, a comparison with darkmatter only simulations is interesting as well, as a good/poor agreement points towards the relevance of baryonic physics. We first compare to the two Millennium simulations, for which growth and merger rates have been derived in Fakhouri et al. (2010). They find that the growth rate is well described by:
(21) 
Using this equation, we find that a halo of mass grows by 38% from redshift 0.6 to 0.22, whilst a halo of mass grows by 43%. More recently, Wetzel & Nagai (2014) use hydrodynamical simulations to measure the amount of physical accretion (i.e. after accounting for pseudoevolution) for haloes with masses . However, for masses in the range , which are typical for LRGs, they use darkmatter only simulations and derive a growth of 10% between =0.5 to =0 at scales smaller than a few hundred kpc. Together with pseudoevolution, the total growth amounts to 43%, comparable to the results of Fakhouri et al. (2010).
Wetzel & Nagai (2014) only study isolated haloes; LRGs cluster strongly, so the typical halo in that work may not be very representative for LRGs. However, no environment selection cuts are made in Fakhouri et al. (2010) and the derived growth rates are similar, hence this seems to be unimportant. In fact, Fakhouri & Ma (2010) measure the growth rate as a function of environment and find a weak trend of a decreasing growth rate towards denser environments, but the statistics for massive haloes is poor due to the low number of haloes.
We find that the growth rates predicted from darkmatter only simulations are marginally consistent with our results. Our measurements suggest a larger growth, particularly towards more massive haloes. This could point toward the impact of baryons; particularly, AGN feedback may have an impact on the distribution of matter in massive haloes (Duffy et al. 2010; Velliscig et al. 2014) and their accretion history. However, we need to improve the statistics of our measurements before we can make stronger claims.
Our conclusion about the evolution of the luminositytohalo mass relation depends, however, on the validity of the assumption of pure passive evolution. The absence of a clear tilt in Figure 3 supports this view. Nonetheless, some residual star formation and/or mergers may also affect the luminosities, which could have an impact on our conclusions.
The amount of star formation that is allowed in LRGs is limited, based on their colour evolution (Wake et al. 2006; Maraston et al. 2009). Also, spectral analyses from SDSSI/II LRGs and CMASS galaxies point toward a very low fraction of galaxies that either form stars or have AGN activity (Greisel et al. 2013; Thomas et al. 2013). There are indications that intermediate mass earlytype galaxies have a low level of ongoing star formation (e.g Kaviraj et al. 2007; Schawinski et al. 2007; Salim & Rich 2010), but it seems unlikely that it is sufficiently strong to affect our conclusions.
Mergers also complicate a purely passive luminosity evolution scenario. The merging history of LRGs or massive earlytype galaxies has received considerable attention in recent years (e.g Tal et al. 2012; LópezSanjuan et al. 2012; Gabor & Davé 2012; Bédorf & Portegies Zwart 2013; Ruiz et al. 2014). We focus here on the results from Tojeiro & Percival (2010), as they estimate the amount of luminosity growth in LRGs due to mergers. Their analysis is based on the measured luminosity function and clustering strength of LRGs in the SDSS in the redshift range , from which they deduce that the average luminosity of LRGs increases by 1.56% Gyr due to mergers, depending on luminosity, such that the growth mainly happens for the faintest LRGs. For LRGs with , the evolution is consistent with passive evolution. In Table 4 of that work, luminosity growth rates from recent works are compared, suggesting that their values are fairly representative.
Tojeiro & Percival (2010) propose two type of mergers that could contribute to the luminosity growth: a merger between an LRG and a small companion, or a merger between two small companions, whose combined luminosity is sufficient to classify it as LRG. In the first scenario, the luminosities of LRGs increase towards lower redshifts; hence they are overestimated compared to the luminosities of LRGs that have evolved through pure passive evolution.
To estimate how that may impact our derived luminositytohalo mass relations and its evolution, we assume that the luminosities of our three LRGs samples increase through mergers with 2, 4 and 6% Gyr, from the bright to the faint bin, respectively, without affecting the halo masses, i.e. we assume that the masstolight ratio of the smaller companions equals zero. This provides us with an upper limit on the bias in the derived halo mass growth. We use the highest redshift slice as our reference, and lower the luminosities of the lower redshift slices, to mimic how the evolution would have looked in the absence of mergers. For example, between and , approximately 3.3 Gyr passed, so we lower the luminosity of the L1z1 bin by a factor . We apply this factor to the average passively evolved luminosities, i.e. we do not apply it to the individual luminosities before the selection of the samples. After this adjustment, we fit the luminositytohalo mass relation again. We find that the retrieved slopes of the luminositytohalo mass relation are within the error bars of our nominal results. The amplitudes move up by 2 for our redshift slice, and by 1 for our slice. The average growth in halo mass becomes %, consistent with our nominal result. We ignored here that mergers also cause growth in mass, hence the LRGs move diagonally rather than horizontally in the luminositytohalo mass plane. If the masstolight ratio of the smaller companions is similar to that of the LRG, they move along the luminositytohalo mass relation and no bias is caused; if this ratio is larger than that of LRGs, the actual growth in halo mass would be smaller than our nominal value.
The second channel for luminosity growth through mergers described by Tojeiro is through the merging of two faint galaxies, which are individually not bright enough to be selected as an LRG. As the formation history differs from that of the typical LRG that formed at high redshift without much activity afterwards, the properties of their haloes may be different, complicating the interpretation of the measured trends.
Apart from that, there are several other physical processes that complicate the interpretation. For example, star formation and mergers may be linked in LRGs (Kaviraj et al. 2011). In addition, some highredshift LRGs may not become lowredshift LRGs, as they may be tidally disrupted or merge with other galaxies. Furthermore, AGN activity could trigger star formation, leading to too blue colours to match the LRG colour selection at lower redshift. The subset of LRGs whose luminositytohalo mass evolution is most easy to interpret is probably the brightest one: Tojeiro et al. (2012) find that the brightest CMASS galaxies are most likely to evolve in SDSSI/II LRGs, and Tojeiro & Percival (2010) derive that the brightest LRGs experience the least amount of luminosity growth from mergers.
To summarise the above: although passive evolution accounts for the bulk of the luminosity evolution of LRGs, there are several processes that complicate this picture. This affects our ability to select and compare LRGs and their progenitors at different redshifts. In order to extract the full information content on the evolution of LRGs that is contained in our measurements, one should compare with numerical simulations that contain all relevant physical processes. No luminosity evolution correction needs to be applied in that case, since there is no need to match the samples at different redshifts, as long as the selection can be matched in observations and simulations. The range in luminosity and redshift of the LRG sample studied in this work would form an ideal observational test case for such a study.
5 Massconcentration relation
As discussed in Section 3, we assume a functional form for the massconcentration relation, i.e. the one from Duffy et al. (2008), but we allow the overall normalisation of this relation to vary in the fit. The resulting constraints for the LOWZ bins are shown in Figure 8, together with the nominal Duffy et al. (2008) model. The constraints on as a function of redshift are shown in Figure 9 and listed in Table 1.
We find that the derived massconcentration relations are lower than the reference model. The largest discrepancy is for our lowest redshift slice, where we find . For the other redshift slices, our results are about 1 below the relation from Duffy et al. (2008). There is no evidence for a trend of with redshift, although the errors of the higher redshift bins are still fairly large. It does suggest that the redshift scaling of this relation is well captured by its functional form (Equation 17).
In Figure 8, we also plot the relation from Macciò et al. (2008), based on dark matter only simulations using the WMAP3 cosmology (Hinshaw et al. 2009), and the relation from Dutton & Macciò (2014), who derived the massconcentration relation with Nbody simulations using the Planck cosmology (Planck Collaboration et al. 2014). All the relations are shown for a redshift of 0.23. These relations derived from dark matter only simulations consistently predict higher concentrations in the mass range that we study.
It is possible that the low normalisation results from choices in our modelling. Therefore, we also derive the normalisation of the massconcentration relation for the halo model runs where we either allow a fraction of the LRGs to be miscentered, where we allow a nonzero satellite fraction, and for the alternative luminosity evolution correction runs. The results are listed in Table 2. For the alternative luminosity evolution run, the constraints on do not change significantly. For the other two runs, the normalisations increase and become consistent with unity, as is shown in Figure 9. By adding either satellites or a miscentered component, we basically shift power from small scales to intermediate scales. Consequently, to fit the same data, the lensing signal of the central LRGs needs to become steeper, i.e. the concentrations need to increase. Note that the errors on increase as well, due to the increased freedom in the model (particularly for the miscentering run).
In Figure 8 we also show the results from Mandelbaum et al. (2008), who derived the massconcentration relation by combining lensing measurements for L type galaxies, galaxy groups traced by LRGs and the maxBCG cluster sample in the SDSS. The model used for their main results did not include a miscentering component, but they used a minimum scale in the fit of 500 kpc for the maxBCG sample to reduce the impact of miscentering. We only show the measurements that overlap with our range of halo masses. We find that our results are consistent with theirs. Towards lower masses, the massconcentration relation of Mandelbaum et al. (2008) is in agreement with theoretical expectations, but at larger masses, the normalisation is lower, as can been seen in their Figure 5. As their lenses span a broad range of mass, they can fit a more flexible massconcentration relation to their measurements, resulting in a broader confidence region compared to ours.
The massconcentration relation of CMASS galaxies is also derived in Miyatake et al. (2013). Their results are very similar to ours: using their fiducial model, in which all galaxies are located at the centre of dark matter haloes, they find a normalisation of with respect to the relation of Macciò et al. (2007). However, when they include a miscentering component identical to ours, the normalisation factor becomes consistent with unity.
The results from Miyatake et al. (2013) and the ones presented here suggest that the low normalisation of the massconcentration relation can be explained by allowing for satellites and/or a miscentering component in the modelling, although the constraints on the miscentering parameters are still very weak. To understand whether the effect is real, one either needs to move towards larger datasets whose corresponding statistical errors are much smaller, or alternatively, one can make use of group catalogues with spectroscopically identified members, to study the brightest group members and the satellites separately. The overlap between the Galaxy And Mass Assembly (GAMA; Driver et al. 2009, 2011; Robotham et al. 2011) and the Kilo Degree Survey (KiDS; Kuijken et al., in prep) could offer an excellent dataset for such a work.
However, it is interesting to consider alternative explanations, such as baryonic processes that can lower the concentrations. AGN feedback, for example, may affect the distribution of dark matter at the relevant radii, as has been studied in Velliscig et al. (2014). In that work, several baryonic processes such as cooling, supernova feedback and AGN feedback are implemented in hydrodynamical simulations to study their effect on the density profiles of haloes. For halo masses typical for LRGs, AGN feedback is the dominant process and affects the dark matter distribution out to several times the virial radius, such that the density is lower at small scales. A similar conclusion was reached in Duffy et al. (2010); whilst cooling leads to more concentrated haloes, this effect is counteracted by AGN feedback, which could cause the concentrations to fall 15% below the ones from dark matter only simulations.
6 Conclusion
We study the evolution of the luminositytohalo mass relation of LRGs, combining SDSS photometry and spectroscopy with the excellent imaging data from the RCS2, which enables us to measure the weak lensing signal up to 0.6. We use stellar population synthesis modelling to compute the correction to account for passive luminosity, which is thought to dominate the luminosity evolution of LRGs, enabling us to compare lowredshift LRGs with their predecessors at higher redshift. We split the LOWZ and CMASS galaxies in two redshift slices and three luminosity bins each, resulting in twelve lens samples, and fit a halo model to the lensing signals of the three luminosity bins simultaneously, for each redshift slice separately. The halo mass estimates that we obtain are broadly consistent with various literature results that are based on a variety of measurement techniques, but span a considerably larger combined range of luminosity and redshift.
We find a typical value of 1.4 for the slope of the luminositytohalo mass relation and no evidence that it changes with redshift. The amplitude of this relation, however, does increase significantly with redshift. We find that the average halo mass of LOWZ galaxies increases by % from to . The halo masses of CMASS galaxies grows by % from to . If CMASS galaxies are the predecessors of the LOWZ ones, the total growth of LRGs from to is %.
This growth in halo mass is somewhat larger than what is expected for pure pseudoevolution, i.e. the evolution in the definition of caused by the change of the mean background density as the universe expands, which by itself causes an apparent halo mass growth of 33% at typical LRG masses between =0.6 and =0.2. Our measurements provide weak support for additional dark matter accretion.
We have tested the sensitivity of these results against changes in the luminosity evolution and changes in the halo model. We find that the inferred slopes of the luminositytohalo mass relation do not significantly change. The amplitude of this relation decreases by at most 2, in the case where we include a satellite component in the halo model. The inferred average growth in halo mass of LRGs does not change by more than 1. Hence for this work, systematic errors are likely subdominant to the statistical ones. For future work that uses measurements with a higher statistical precision, advances in the modelling of the luminosity evolution and of the set up of the halo model, are necessary to avoid biases in the results. Such measurements will be highly valuable to constrain the impact of baryonic processes on the distribution of dark matter, which is essential for a correct and optimized exploitation of future cosmic shear surveys such as Euclid (Laureijs et al. 2011).
We also constrain the overall normalisation of the massconcentration relation for each of the four redshift slices. We find that for our lowest redshift slice, the bestfit relation is lower than what is expected from dark matter only simulations. However, if we allow for miscentering or for the contribution from satellites in the halo model, the normalisation increases and becomes consistent with the results from dark matteronly simulations.
Acknowledgements
We would like to thank Alexie Leauthaud and Claudia Maraston for useful discussions and suggestions. This work was supported by a grant from the German Space Agency DLR. HH, MC and RH acknowledge support from NWO VIDI grant number 639.042.814 and ERC FP7 grant 279396.
This work is based on observations obtained with MegaPrime/MegaCam, a joint project of CFHT and CEA/DAPNIA, at the CanadaFranceHawaii Telescope (CFHT) which is operated by the National Research Council (NRC) of Canada, the Institute National des Sciences de l’Univers of the Centre National de la Recherche Scientifique of France, and the University of Hawaii. We used the facilities of the Canadian Astronomy Data Centre operated by the NRC with the support of the Canadian Space Agency.
Funding for SDSSIII has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, and the U.S. Department of Energy Office of Science. The SDSSIII web site is http://www.sdss3.org/.
SDSSIII is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSSIII Collaboration including the University of Arizona, the Brazilian Participation Group, Brookhaven National Laboratory, University of Cambridge, Carnegie Mellon University, University of Florida, the French Participation Group, the German Participation Group, Harvard University, the Instituto de Astrofisica de Canarias, the Michigan State/Notre Dame/JINA Participation Group, Johns Hopkins University, Lawrence Berkeley National Laboratory, Max Planck Institute for Astrophysics, Max Planck Institute for Extraterrestrial Physics, New Mexico State University, New York University, Ohio State University, Pennsylvania State University, University of Portsmouth, Princeton University, the Spanish Participation Group, University of Tokyo, University of Utah, Vanderbilt University, University of Virginia, University of Washington, and Yale University.
Appendix A Luminosity evolution scatter
The evolution corrected luminosities of galaxies become increasingly uncertain with redshift. This is mainly due to the intrinsic variation in properties of the LRGs, as illustrated in the bottom panel of Figure 2, where the difference between the correction curves increases with redshift for the different models.
To obtain a rough estimate of how intrinsic scatter in the ()corrected absolute magnitudes affects our results, we do a simple test. We start with randomly drawing redshifts and ()corrected absolute magnitudes from the original distribution, i.e. the one that is shown in the bottom panel of Figure 3. We assume that these magnitudes are the intrinsic ones. Next, we assign a mass to each object using our nominal bestfit luminositytohalo mass relation from Table 2. Note that our conclusions do not sensitively depend on the choice of slope and offset. We compute an NFW profile for each object at 100 logarithmically spaced radial bins between 0.05 h Mpc and 1 h Mpc. We stack the NFW profiles of objects that fall inside a lens bin as defined in Figure 3, and fit an NFW to the resulting profile using the mean redshift of those objects.
Next, we simulate the intrinsic scatter in absolute magnitude by assuming that it can be described by a Gaussian whose width increases as . This particular choice covers the range of correction curves of the SSP models that are in reasonable agreement with the colours of the LRGs, as shown in Figure 2. We draw a random number from this Gaussian and add that to the ‘intrinsic’ magnitudes: these are the ‘observed’ magnitudes. We stack the NFW profiles, but now using the ‘observed’ magnitudes to select the lenses. Again, we fit an NFW profile to the stacked profile using the mean redshift of the observed objects. The impact of intrinsic luminosity scatter is then estimated from the ratio of this ‘observed’ mass to the intrinsic one. We find that the mass typically changes by a few percent. The largest difference is found for the L3z4 bin, where the difference is 10%. However, in all cases, the error is considerably smaller than the statistical errors. We therefore conclude that intrinsic scatter of the luminosities can be safely ignored.
Appendix B Systematic tests
Since this work aims at measuring the redshift dependence of the lensing signal, we have to be particularly careful with systematic effects that mimic a redshift scaling. One issue that requires attention is the computation of the mean lensing efficiency. For example, if we miss the high redshift tail of the source redshift distribution, the mean lensing efficiency would be biased increasingly low for higher lens redshifts, which would mimic a redshift dependence in the lensing signal. Motivated by the differences in mean lensing efficiency computed from the photometric redshift catalogue from Ilbert et al. (2009) and Ilbert et al. (2013) at high lens redshifts, we perform the following test.
We divide our source sample in six different samples, with magnitude cuts of , , and , , , . We split the LRGs in a lowredshift sample with and a highredshift one with , and measure their lensing signals using each source sample. Note that for each source sample, we separately determine the shear signal around random points, the mean lensing efficiencies using the Ilbert et al. (2013) photometric redshift catalogues, the contamination of source galaxies around the lenses, and the noise bias correction. We measure the weighted mean of the lensing signal times the projected separation, as that is roughly a constant with radius, over the range 0.15 10 h Mpc. We show these values for the two LRG samples in Figure 10.
We find that the weighted mean of the lensing signal of each lens sample is consistent for the various source samples. Clearly, the measurements are correlated for the source samples that overlap in apparent magnitude, but the four bins of 0.5 magnitude width are more or less independent (not completely as they are similarly affected by cosmic variance). This result strongly suggests that the measurement process is robust.
We have repeated this test using the Ilbert et al. (2009) photometric redshift catalogue instead. For the LRG sample the results are consistent, but for the highz one, some bins differ by 23. This points at a problem with the source redshift distribution used to compute the lensing efficiencies. Since the Ilbert et al. (2013) photometric redshift catalogue gives consistent results for the different source samples, even at high redshifts, this suggests that the source redshift distribution is sufficiently accurately determined with the latter, but not with the former (see also the discussion in Hoekstra et al. 2015) .
Appendix C Fitting Methodology
We use Bayesian inference techniques to determine the posterior probability distribution of the model parameters , given the data . According to Bayes’ theorem,
(22) 
where is the likelihood of the data given the model parameters, is the prior probability of these parameters, and
(23) 
is called the evidence. Since we do not intend to perform model selection, the evidence just acts as a normalisation constant which needs not to be calculated. Therefore, the posterior distribution is given by
(24) 
where
(25) 
denotes the model prediction, is the corresponding error, and for the low and high redshift bins of the LOWZ and CMASS samples, respectively. For our fiducial model, the set of model parameters is , where:

is the mean halo mass for the th luminosity bin;

is the scatter of ;

is the normalisation of the relation.
When exploring model variations (see Section 4.1) we employ and . For all model parameters, we adopt a flat, sufficiently wide prior such that the results are not biased.
We sample the posterior distribution of our model parameters given the data using a MonteCarlo Markov chain (MCMC). In particular, we implement the MetropolisHastings algorithm to construct the MCMC. At any point in the chain, a trial model is generated using a method specified below. The statistic for the trial model, , is calculated using Equation (25). This trial model is accepted to be a member of the chain with a probability given by
(26) 
where denotes the for the current model in the chain. We initialize the chain from a random position in our multidimensional parameter space and obtain a chain of 500,000 models. We discard the first models (the burnin period) allowing the chain to sample from a more probable part of the distribution. We use this chain of models to estimate the confidence levels on the parameters and on the lensing signal, as shown in Figure 4.
A proper choice of the proposal distribution is very important in order to achieve fast convergence and a reasonable acceptance rate for the trial models. The posterior distribution in a multidimensional parameter space, such as the one we are dealing with, will have degeneracies and in general can be very difficult to sample from. We have adopted the following strategy to overcome these difficulties. During the first half of the burnin stage, we chose an independent Gaussian proposal distribution for every model parameter, as is common for the MetropolisHastings algorithm. Halfway through the burnin stage, we perform a Fisher information matrix analysis at the bestfit model found thus far. The Fisher information matrix, given by
(27) 
is a symmetric matrix, where denotes the number of parameters in our model, and is the likelihood. The inverse of the Fisher matrix gives the covariance matrix, , of the posterior constraints on the model parameters. More importantly, the eigenvectors of the covariance matrix are an excellent guide to the degeneracies in the posterior distribution, and the corresponding eigenvalues set a scale for how wide the posterior ought to be in a given direction. Therefore, for the second half of the burnin period, we utilize this information and use a proposal distribution which is a multivariate Gaussian centered at the current value of the parameters and with a covariance given by the aforementioned matrix. In practice, the trial model () can be generated from the current model () using
(28) 
where is a vector consisting of standard normal deviates, the matrix is such that , and is a parameter that we have chosen to achieve an average acceptance rate of 30%. We repeat the Fisher matrix analysis once again at the end of the burnin period (using the bestfit model found thus far) and use the covariance matrix to define our proposal distribution to be used for the MCMC. This strategy has proven to be extremely efficient in sampling posterior distributions for similar studies (see e.g. Cacciato et al. 2013, 2014).
Appendix D kcorrected only results
In this work, we have attempted to account for the evolution of the luminosities of LRGs by applying a luminosity evolution correction. Unfortunately, this correction is uncertain and will also add scatter to the corrected luminosities as LRGs cover a range of intrinsic properties such as formation age and metallicity, which complicates the interpretation of the results. Therefore, it is also interesting to apply only the correction to the absolute magnitudes. The resulting luminosities are closer to the real luminosities, i.e. they have smaller scatter. Interpreting any trend in the luminositytohalo mass relation will be harder as both the luminosity and the halo mass may evolve simultaneously, but comparison to for example simulations should be more straightforward.
Hence we divide the lens sample into bins of corrected absolute magnitude and redshift. Details of the lens samples can be found in Table 3. The selection is also illustrated in Figure 3. For each subsample, we stack the lensing signal of all the lenses inside that bin and show that in Figure 11. We fit the halo model using the same setup as before and show the bestfit models, together with the model uncertainties, in Figure 11. The corresponding effective masses and normalisations of the massconcentration relation can be found in Table 3.
M  N  
(1)  (2)  (3)  (4)  (5)  (6)  
0.150.29 (LOWZ)  
L1z1  [21.8,21.2]  1110  0.204  0.68  
L2z1  [22.4,21.8]  3784  0.225  1.04  
L3z1  [22.8,22.4]  954  0.234  1.59  
0.290.43 (LOWZ)  
L2z2  [22.4,21.8]  4725  0.351  1.10  
L3z2  [22.8,22.4]  2901  0.368  1.62  
L4z2  [23.2,22.8]  981  0.374  2.31  
0.430.55 (CMASS)  
L2z3  [22.4,21.8]  7798  0.496  1.09  
L3z3  [22.8,22.4]  4197  0.506  1.62  
L4z3  [23.2,22.8]  1646  0.506  2.31  
0.550.70 (CMASS)  
L2z4  [22.4,21.8]  2532  0.584  1.17  
L3z4  [22.8,22.4]  4357  0.599  1.67  
L4z4  [23.2,22.8]  3365  0.615  2.35  
L5z4  [23.6,23.2]  1314  0.625  3.37 
Footnotes
 institutetext: ArgelanderInstitut für Astronomie, Auf dem Hügel 71, 53121 Bonn, Germany
 institutetext: Leiden Observatory, Leiden University, Niels Bohrweg 2, NL2333 CA Leiden, The Netherlands
 http://www.sdss3.org/dr9/algorithms/boss_galaxy_ts.php
 https://github.com/GalSimdevelopers/GalSim
 http://www.mpagarching.mpg.de/SDSS/DR7/
 (1) absolute magnitude range (after ()correction); (2) number of lenses; (3) mean redshift; (4) mean luminosity [10] (after ()correction); (5) bestfit halo mass [10]; (6) bestfit normalisation of the massconcentration relation; (7) reduced chisquared of the fit.
 (1) absolute magnitude range (after correction); (2) number of lenses; (3) mean redshift; (4) mean luminosity [10