• We present the first weak lensing calibration of $\mu_{\star}$, a new galaxy cluster mass proxy corresponding to the total stellar mass of red and blue members, in two cluster samples selected from the SDSS Stripe 82 data: 230 redMaPPer clusters at redshift $0.1\leq z<0.33$ and 136 Voronoi Tessellation (VT) clusters at $0.1 \leq z < 0.6$. We use the CS82 shear catalog and stack the clusters in $\mu_{\star}$ bins to measure a mass-observable power law relation. For redMaPPer clusters we obtain $M_0 = (1.77 \pm 0.36) \times 10^{14}h^{-1} M_{\odot}$, $\alpha = 1.74 \pm 0.62$. For VT clusters, we find $M_0 = (4.31 \pm 0.89) \times 10^{14}h^{-1} M_{\odot}$, $\alpha = 0.59 \pm 0.54$ and $M_0 = (3.67 \pm 0.56) \times 10^{14}h^{-1} M_{\odot}$, $\alpha = 0.68 \pm 0.49$ for a low and a high redshift bin, respectively. Our results are consistent, internally and with the literature, indicating that our method can be applied to any cluster finding algorithm. In particular, we recommend that $\mu_{\star}$ be used as the mass proxy for VT clusters. Catalogs including $\mu_{\star}$ measurements will enable its use in studies of galaxy evolution in clusters and cluster cosmology.
  • We study projected underdensities in the cosmic galaxy density field known as 'troughs', and their overdense counterparts, which we call 'ridges'. We identify these regions using a bright sample of foreground galaxies from the photometric Kilo-Degree Survey (KiDS), specifically selected to mimic the spectroscopic Galaxy And Mass Assembly survey (GAMA). From an independent sample of KiDS background galaxies, we measure the weak gravitational lensing profiles of the troughs/ridges. We quantify their lensing strength $A$ as a function of galaxy density percentile rank $P$ and overdensity $\delta$, and find that the skewness in the galaxy density distribution is reflected in the total mass distribution measured by weak lensing. We interpret our results using the mock galaxy catalogue from the Marenostrum Institut de Ci\`encies de l'Espai (MICE) Grand Challenge lightcone simulation, and find a good agreement with our observations. Using signal-to-noise weights derived from the Scinet LIghtCone Simulations (SLICS) mock catalogue we optimally stack the lensing signal of KiDS troughs with an angular radius $\theta_A$ = {5,10,15,20} arcmin, resulting in {16.8,14.9,10.13,7.55} $\sigma$ detections. Finally, we select troughs using a volume-limited sample of galaxies, split into two redshift bins between 0.1 < z < 0.3. For troughs/ridges with transverse comoving radius $R_A$ = 1.9 Mpc/h, we find no significant difference between the comoving $A'(P)$ and $A'(\delta)$ relation of the low- and high-redshift sample. Using the MICE and SLICS mocks we predict that trough and ridge evolution could be detected with gravitational lensing using deeper and wider lensing surveys, such as those from the Large Synoptic Survey Telescope and Euclid.
  • We present cosmological parameter constraints from a joint analysis of three cosmological probes: the tomographic cosmic shear signal in $\sim$450 deg$^2$ of data from the Kilo Degree Survey (KiDS), the galaxy-matter cross-correlation signal of galaxies from the Galaxies And Mass Assembly (GAMA) survey determined with KiDS weak lensing, and the angular correlation function of the same GAMA galaxies. We use fast power spectrum estimators that are based on simple integrals over the real-space correlation functions, and show that they are practically unbiased over relevant angular frequency ranges. We test our full pipeline on numerical simulations that are tailored to KiDS and retrieve the input cosmology. By fitting different combinations of power spectra, we demonstrate that the three probes are internally consistent. For all probes combined, we obtain $S_8\equiv \sigma_8 \sqrt{\Omega_{\rm m}/0.3}=0.800_{-0.027}^{+0.029}$, consistent with Planck and the fiducial KiDS-450 cosmic shear correlation function results. Marginalising over wide priors on the mean of the tomographic redshift distributions yields consistent results for $S_8$ with an increase of $28\%$ in the error. The combination of probes results in a $26\%$ reduction in uncertainties of $S_8$ over using the cosmic shear power spectra alone. The main gain from these additional probes comes through their constraining power on nuisance parameters, such as the galaxy intrinsic alignment amplitude or potential shifts in the redshift distributions, which are up to a factor of two better constrained compared to using cosmic shear alone, demonstrating the value of large-scale structure probe combination.
  • We demonstrate that deep good-seeing VLT/HAWK-I $K_\mathrm{s}$ images complemented with $g$+$z$-band photometry can yield a sensitivity for weak lensing studies of massive galaxy clusters at redshifts \mbox{$0.7\lesssim z \lesssim 1.1$}, which is almost identical to the sensitivity of HST/ACS mosaics of single-orbit depth. Key reasons for this good performance are the excellent image quality frequently achievable for $K_\mathrm{s}$ imaging from the ground, a highly effective photometric selection of background galaxies, and a galaxy ellipticity dispersion that is noticeably lower than for optically observed high-redshift galaxy samples. Incorporating results from the 3D-HST and UltraVISTA surveys we also obtained a more accurate calibration of the source redshift distribution than previously achieved for similar optical weak lensing data sets. Here we studied the extremely massive galaxy cluster RCS2$J$232727.7$-$020437 (\mbox{$z=0.699$}), combining deep VLT/\mbox{HAWK-I} $K_\mathrm{s}$ images (point spread function with a 0\farcs35 full width at half maximum) with LBT/LBC photometry. The resulting weak lensing mass reconstruction suggests that the cluster consists of a single overdensity, which is detected with a peak significance of $10.1\sigma$. We constrained the cluster mass to \mbox{$M_\mathrm{200c}/(10^{15} \mathrm{M}_\odot) =2.06^{+0.28}_{-0.26}(\mathrm{stat.})\pm 0.12 (\mathrm{sys.})$} assuming a spherical Navarro, Frenk \& White model and simulation-based priors on the concentration, making it one of the most massive galaxy clusters known in the \mbox{$z\gtrsim 0.7$} Universe. We also cross-checked the HAWK-I measurements through an analysis of overlapping HST/ACS images, yielding fully consistent estimates of the lensing signal.
  • We conduct a comprehensive study of the effects of incorporating galaxy morphology information in photometric redshift estimation. Using machine learning methods, we assess the changes in the scatter and catastrophic outlier fraction of photometric redshifts when galaxy size, ellipticity, S\'{e}rsic index and surface brightness are included in training on galaxy samples from the SDSS and the CFHT Stripe-82 Survey (CS82). We show that by adding galaxy morphological parameters to full $ugriz$ photometry, only mild improvements are obtained, while the gains are substantial in cases where fewer passbands are available. For instance, the combination of $grz$ photometry and morphological parameters almost fully recovers the metrics of $5$-band photometric redshifts. We demonstrate that with morphology it is possible to determine useful redshift distribution $N(z)$ of galaxy samples without any colour information. We also find that the inclusion of quasar redshifts and associated object sizes in training improves the quality of photometric redshift catalogues, compensating for the lack of a good star-galaxy separator. We further show that morphological information can mitigate biases and scatter due to bad photometry. As an application, we derive both point estimates and posterior distributions of redshifts for the official CS82 catalogue, training on morphology and SDSS Stripe-82 $ugriz$ bands when available. Our redshifts yield a 68th percentile error of $0.058(1+z)$, and a catastrophic outlier fraction of $5.2$ per cent. We further include a deep extension trained on morphology and single $i$-band CS82 photometry.
  • We measure the projected galaxy clustering and galaxy-galaxy lensing signals using the Galaxy And Mass Assembly (GAMA) survey and Kilo-Degree Survey (KiDS) to study galaxy bias. We use the concept of non-linear and stochastic galaxy biasing in the framework of halo occupation statistics to constrain the parameters of the halo occupation statistics and to unveil the origin of galaxy biasing. The bias function $\Gamma_{\text{gm}}(r_{\text{p}})$, where $r_{\text{p}}$ is the projected comoving separation, is evaluated using the analytical halo model from which the scale dependence of $\Gamma_{\text{gm}}(r_{\text{p}})$, and the origin of the non-linearity and stochasticity in halo occupation models can be inferred. Our observations unveil the physical reason for the non-linearity and stochasticity, further confirmed using hydrodynamical simulations, with the stochasticity mostly originating from the non-Poissonian behaviour of satellite galaxies in the dark matter haloes and their spatial distribution, which does not follow the spatial distribution of dark matter in the halo. The observed non-linearity is mostly due to the presence of the central galaxies, as was noted from previous theoretical work on the same topic. We also see that overall, more massive galaxies reveal a stronger scale dependence, and out to a larger radius. Our results show that a wealth of information about galaxy bias is hidden in halo occupation models. These models should therefore be used to determine the influence of galaxy bias in cosmological studies.
  • The spatial distribution of satellite galaxies encodes rich information of the structure and assembly history of galaxy clusters. In this paper, we select a redMaPPer cluster sample in SDSS Stripe 82 region with $0.1 \le z \le 0.33$, $20<\lambda<100$ and $P_{\rm cen}>0.7$. Using the high-quality weak lensing data from CS82 Survey, we constrain the mass profile of this sample. Then we compare directly the mass density profile with the satellite number density profile. We find that the total mass and number density profiles have the same shape, both well fitted by an NFW profile. The scale radii agree with each other within 1$\sigma$ error ($r_{\rm s,gal}=0.34_{-0.03}^{+0.04}$Mpc vs $r_{\rm s}=0.37_{-0.10}^{+0.15}$Mpc ).
  • We present a new measurement of the mass-concentration relation and the stellar-to-halo mass ratio over the halo mass range $5\times 10^{12}$ to $2\times 10^{14}M_{\odot}$. To achieve this, we use weak lensing measurements from the CFHT Stripe 82 Survey (CS82), combined with the central galaxies from the redMaPPer cluster catalogue and the LOWZ/CMASS galaxy sample of the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Tenth Data Release. The stacked lensing signals around these samples are modelled as a sum of contributions from the central galaxy, its dark matter halo, and the neighboring halos, as well as a term for possible centering errors. We measure the mass-concentration relation: $c_{200c}(M)=A(\frac{M_{200c}}{M_0})^{B}$ with $A=5.24\pm1.24, B=-0.13\pm0.10$ for $0.2<z<0.4$ and $A=6.61\pm0.75, B=-0.15\pm0.05$ for $0.4<z<0.6$. These amplitudes and slopes are completely consistent with predictions from recent simulations. We also measure the stellar-to-halo mass ratio for our samples, and find results consistent with previous measurements from lensing and other techniques.
  • This paper is the first of a series of papers constraining cosmological parameters with weak lensing peak statistics using $\sim 450~\rm deg^2$ of imaging data from the Kilo Degree Survey (KiDS-450). We measure high signal-to-noise ratio (SNR: $\nu$) weak lensing convergence peaks in the range of $3<\nu<5$, and employ theoretical models to derive expected values. These models are validated using a suite of simulations. We take into account two major systematic effects, the boost factor and the effect of baryons on the mass-concentration relation of dark matter haloes. In addition, we investigate the impacts of other potential astrophysical systematics including the projection effects of large scale structures, intrinsic galaxy alignments, as well as residual measurement uncertainties in the shear and redshift calibration. Assuming a flat $\Lambda$CDM model, we find constraints for $S_{\rm 8}=\sigma_{\rm 8}(\Omega_{\rm m}/0.3)^{0.5}=0.746^{+0.046}_{-0.107}$ according to the degeneracy direction of the cosmic shear analysis and $\Sigma_{\rm 8}=\sigma_{\rm 8}(\Omega_{\rm m}/0.3)^{0.38}=0.696^{+0.048}_{-0.050}$ based on the derived degeneracy direction of our high-SNR peak statistics. The difference between the power index of $S_{\rm 8}$ and in $\Sigma_{\rm 8}$ indicates that combining the two probes has the potential to break the degeneracy in $\sigma_{\rm 8}$ and $\Omega_{\rm m}$. Our results are consistent with the cosmic shear tomographic correlation analysis of the same dataset and $\sim 2\sigma$ lower than the Planck 2016 results.
  • We study the statistics of peaks in a weak lensing reconstructed mass map of the first 450 square degrees of the Kilo Degree Survey. The map is computed with aperture masses directly applied to the shear field with an NFW-like compensated filter. We compare the peak statistics in the observations with that of simulations for various cosmologies to constrain the cosmological parameter $S_8 = \sigma_8 \sqrt{\Omega_{\rm m}/0.3}$, which probes the ($\Omega_{\rm m}, \sigma_8$) plane perpendicularly to its main degeneracy. We estimate $S_8=0.750\pm0.059$, using peaks in the signal-to-noise range $0 \leq {\rm S/N} \leq 4$, and accounting for various systematics, such as multiplicative shear bias, mean redshift bias, baryon feedback, intrinsic alignment, and shear-position coupling. These constraints are $\sim25\%$ tighter than the constraints from the high significance peaks alone ($3 \leq {\rm S/N} \leq 4$) which typically trace single-massive halos. This demonstrates the gain of information from low-S/N peaks which correspond to the projection of several small-mass halos along the line-of-sight. Our results are in good agreement with the tomographic shear two-point correlation function measurement in KiDS-450. Combining shear peaks with non-tomographic measurements of the shear two-point correlation functions yields an $\sim20\%$ improvement in the uncertainty on $S_8$ compared to the shear two-point correlation functions alone, highlighting the great potential of peaks as a cosmological probe.
  • In the hierarchical formation model, galaxy clusters grow by accretion of smaller groups or isolated galaxies. During the infall into the centre of a cluster, the properties of accreted galaxies change. In particular, both observations and numerical simulations suggest that its dark matter halo is stripped by the tidal forces of the host. We use galaxy-galaxy weak lensing to measure the average mass of dark matter haloes of satellite galaxies as a function of projected distance to the centre of the host, for different stellar mass bins. Assuming that the stellar component of the galaxy is less disrupted by tidal stripping, stellar mass can be used as a proxy of the infall mass. We study the stellar to halo mass relation of satellites as a function of the cluster-centric distance to measure tidal stripping. We use the shear catalogues of the DES science verification archive, the CFHTLenS and the CFHT Stripe 82 (CS82) surveys, and we select satellites from the redMaPPer catalogue of clusters. For galaxies located in the outskirts of clusters, we find a stellar to halo mass relation in good agreement with the theoretical expectations from \citet{moster2013} for central galaxies. In the centre of the cluster, we find that this relation is shifted to smaller halo mass for a given stellar mass. We interpret this finding as further evidence for tidal stripping of dark matter haloes in high density environments.
  • In recent years, multiple studies have reported substantial populations of large, low-surface-brightness galaxies in local galaxy clusters. Various theories that aim to explain the presence of such ultra-diffuse galaxies (UDGs) have since been proposed. A key question that will help to differentiate between models is whether UDGs have counterparts in lower-mass host haloes, and what their abundance as a function of halo mass is. In this study we extend our previous study of UDGs in galaxy clusters to galaxy groups. We measure the abundance of UDGs in 325 spectroscopically-selected groups from the Galaxy And Mass Assembly (GAMA) survey. We make use of the overlapping imaging from the ESO Kilo-Degree Survey (KiDS), from which we can identify galaxies with mean surface brightnesses within their effective radii down to ~25.5 mag arcsec$^{-2}$ in the r-band. We are able to measure a significant overdensity of UDGs (with sizes r_eff > 1.5 kpc) in galaxy groups down to M200=10^12 Msun, a regime where approximately only 1 in 10 groups contains a UDG that we can detect. We combine measurements of the abundance of UDGs in haloes that cover three orders of magnitude in halo mass, finding that their numbers scale quite steeply with halo mass; N_UDG (R<R200) $\propto$ M200^(1.11+/-0.07). To better interpret this, we also measure the mass-richness relation for brighter galaxies down to $M^*_r$+2.5 in the same GAMA groups, and find a much shallower relation of N_Bright (R<R200) $\propto$ M200^(0.78+/-0.05). This shows that UDGs are relatively more abundant, compared to bright galaxies, in massive clusters than in groups. We discuss implications, but whether this difference is related to a higher destruction rate of UDGs in groups, or whether massive haloes have a positive effect on their formation, is not yet clear.
  • Galaxy systems with large magnitude gaps - defined as the difference in magnitude between the central galaxy and the brightest satellite in the central region - have been claimed to have earlier formation histories. In this study we measure the mass concentration, as an indicator of early formation, of ensembles of galaxy systems divided by redshift and magnitude gaps in the $r$ band. We use cross-correlation weak lensing measurements with NFW parametric mass profiles to measure masses and concentrations of these ensembles from a catalogue of systems built from the SDSS Coadd by the redMaPPer algorithm. The lensing shear data come from the CFHT Stripe 82 (CS82) survey, and consists of $i$-band images of the SDSS Stripe 82 region. We find that the stack made up of systems with larger magnitude gaps has a high probability of being more concentrated on average, in the lowest redshift slice ($0.2<z<0.4$), both when dividing in quartiles ($P=0.98$) and terciles ($P=0.85$). These results lend credibility to the claim that the magnitude gap is an indicator of earlier formed systems.
  • We perform a combined analysis of cosmic shear tomography, galaxy-galaxy lensing tomography, and redshift-space multipole power spectra (monopole and quadrupole) using 450 deg$^2$ of imaging data by the Kilo Degree Survey (KiDS) overlapping with two spectroscopic surveys: the 2-degree Field Lensing Survey (2dFLenS) and the Baryon Oscillation Spectroscopic Survey (BOSS). We restrict the galaxy-galaxy lensing and multipole power spectrum measurements to the overlapping regions with KiDS, and self-consistently compute the full covariance between the different observables using a large suite of $N$-body simulations. We methodically analyze different combinations of the observables, finding that galaxy-galaxy lensing measurements are particularly useful in improving the constraint on the intrinsic alignment amplitude (by 30%, positive at $3.5\sigma$ in the fiducial data analysis), while the multipole power spectra are useful in tightening the constraints along the lensing degeneracy direction (e.g. factor of two stronger matter density constraint in the fiducial analysis). The fully combined constraint on $S_8 \equiv \sigma_8 \sqrt{\Omega_{\rm m}/0.3} = 0.742 \pm 0.035$, which is an improvement by 20% compared to KiDS alone, corresponds to a $2.6\sigma$ discordance with Planck, and is not significantly affected by fitting to a more conservative set of scales. Given the tightening of the parameter space, we are unable to resolve the discordance with an extended cosmology that is simultaneously favored in a model selection sense, including the sum of neutrino masses, curvature, evolving dark energy, and modified gravity. The complementarity of our observables allow for constraints on modified gravity degrees of freedom that are not simultaneously bounded with either probe alone, and up to a factor of three improvement in the $S_8$ constraint in the extended cosmology compared to KiDS alone.
  • We build a background cluster candidate catalog from the Next Generation Virgo Cluster Survey, using our detection algorithm RedGOLD. The NGVS covers 104$deg^2$ of the Virgo cluster in the $u^*,g,r,i,z$-bandpasses to a depth of $ g \sim 25.7$~mag (5$\sigma$). Part of the survey was not covered or has shallow observations in the $r$--band. We build two cluster catalogs: one using all bandpasses, for the fields with deep $r$--band observations ($\sim 20 \ deg^2$), and the other using four bandpasses ($u^*,g,i,z$) for the entire NGVS area. Based on our previous CFHT-LS W1 studies, we estimate that both of our catalogs are $\sim100\%$($\sim70\%$) complete and $\sim80\%$ pure, at $z\le 0.6$($z\lesssim1$), for galaxy clusters with masses of $M\gtrsim10^{14}\ M_{\odot}$. We show that when using four bandpasses, though the photometric redshift accuracy is lower, RedGOLD detects massive galaxy clusters up to $z\sim 1$ with completeness and purity similar to the five-band case. This is achieved when taking into account the bias in the richness estimation, which is $\sim40\%$ lower at $0.5\le z<0.6$ and $\sim20\%$ higher at $0.6<z< 0.8$, with respect to the five-band case. RedGOLD recovers all the X-ray clusters in the area with mass $M_{500} > 1.4 \times 10^{14} \rm M_{\odot}$ and $0.08<z<0.5$. Because of our different cluster richness limits and the NGVS depth, our catalogs reach to lower masses than the published redMaPPer cluster catalog over the area, and we recover $\sim 90-100\%$ of its detections.
  • We set out to quantify the number density of quiescent massive compact galaxies at intermediate redshifts. We determine structural parameters based on i-band imaging using the CFHT equatorial SDSS Stripe 82 (CS82) survey (~170 sq. degrees) taking advantage of an exquisite median seeing of ~0.6''. We select compact massive (M > 5x10^10 M_sun) galaxies within the redshift range of 0.2<z<0.6. The large volume sampled allows to decrease the effect of cosmic variance that has hampered the calculation of the number density for this enigmatic population in many previous studies. We undertake an exhaustive analysis in an effort to untangle the various findings inherent to the diverse definition of compactness present in the literature. We find that the absolute number of compact galaxies is very dependent on the adopted definition and can change up to a factor of >10. We systematically measure a factor of ~5 more compacts at the same redshift than what was previously reported on smaller fields with HST imaging, which are more affected by cosmic variance. This means that the decrease in number density from z ~ 1.5 to z ~ 0.2 might be only of a factor of ~2-5, significantly smaller than what previously reported. This supports progenitor bias as the main contributor to the size evolution. This milder decrease is roughly compatible with the predictions from recent numerical simulations. Only the most extreme compact galaxies, with Reff < 1.5x( M/10^11 M_sun)^0.75 and M > 10^10.7 M_sun, appear to drop in number by a factor of ~20 and hence likely experience a noticeable size evolution.
  • We measured stacked weak lensing cluster masses for a sample of 1325 galaxy clusters detected by the RedGOLD algorithm in the Canada-France-Hawaii Telescope Legacy Survey W1 and the Next Generation Virgo Cluster Survey at $0.2<z<0.5$, in the optical richness range $10<\lambda<70$. After a selection of our best richness subsample ($20<\lambda<50$), this is the most comprehensive lensing study of a $\sim 100\%$ complete and $\sim 90\%$ pure optical cluster catalogue in this redshift range, with a total of 346 clusters in $\sim164~deg^2$. We test three different mass models, and our best model includes a basic halo model, with a Navarro Frenk and White profile, and correction terms that take into account cluster miscentering, non-weak shear, the two-halo term, the contribution of the Brightest Cluster Galaxy, and an a posteriori correction for the intrinsic scatter in the mass-richness relation. With this model, we obtain a mass-richness relation of $\log{M_{\rm 200}/M_{\odot}}=(14.48\pm0.04)+(1.14\pm0.23)\log{(\lambda/40)}$ (statistical uncertainties). This result is consistent with other published lensing mass-richness relations. When compared to X-ray masses and mass proxies, we find that on average weak lensing masses are $\sim 10\%$ higher than those derived in the X-ray in the range $2\times10^{13}M_{\rm \odot}<E(z) M^{X}_{\rm 200}<2\times10^{14}M_{\rm \odot}$, in agreement with most previous results and simulations. We also give the coefficients of the scaling relations between the lensing mass and X-ray mass proxies, $L_X$ and $T_X$, and compare them with previous results.
  • RedGOLD searches for red-sequence galaxy overdensities while minimizing contamination from dusty star-forming galaxies. It imposes an NFW profile and calculates cluster detection significance and richness. We optimize these latter two parameters using both simulations and X-ray detected cluster catalogs, and obtain a catalog $\sim 80\%$ pure up to $z \sim 1$, and $\sim 100\%$ ($\sim 70\%$) complete at $z\le 0.6$ ( $z\lesssim1$) for galaxy clusters with $M \gtrsim 10^{14}\ {\rm M_{\odot}}$ at the CFHT-LS Wide depth. In the CFHT-LS W1, we detect 11 cluster candidates per $\rm deg^2$ out to $z\sim1.1$. When we optimize both completeness and purity, RedGOLD obtains a cluster catalog with higher completeness and purity than other public catalogs, obtained using CFHT-LS W1 observations, for $M \gtrsim 10^{14}\ {\rm M_{\odot}}$. We use X-ray detected cluster samples to extend the study of the X-ray temperature-optical richness relation to a lower mass threshold, and find a mass scatter at fixed richness of $\sigma_{lnM|\lambda}=0.39\pm0.07$ and $\sigma_{lnM|\lambda}=0.30\pm0.13$ for the Gozaliasl et al. (2014) and Mehrtens et al. (2012) samples. When considering similar mass ranges as previous work, we recover a smaller scatter in mass at fixed richness. We recover $93\%$ of the redMaPPer detections, and find that its richness estimates is on average $\sim 40-50\%$ larger than ours at $z>0.3$. RedGOLD recovers X-ray cluster spectroscopic redshifts at better than $5\%$ up to $z\sim1$, and the centers within a few tens of arcseconds.
  • We test extensions to the standard cosmological model with weak gravitational lensing tomography using 450 deg$^2$ of imaging data from the Kilo Degree Survey (KiDS). In these extended cosmologies, which include massive neutrinos, nonzero curvature, evolving dark energy, modified gravity, and running of the scalar spectral index, we also examine the discordance between KiDS and cosmic microwave background measurements from Planck. The discordance between the two datasets is largely unaffected by a more conservative treatment of the lensing systematics and the removal of angular scales most sensitive to nonlinear physics. The only extended cosmology that simultaneously alleviates the discordance with Planck and is at least moderately favored by the data includes evolving dark energy with a time-dependent equation of state (in the form of the $w_0-w_a$ parameterization). In this model, the respective $S_8 = \sigma_8 \sqrt{\Omega_{\rm m}/0.3}$ constraints agree at the $1\sigma$ level, and there is `substantial concordance' between the KiDS and Planck datasets when accounting for the full parameter space. Moreover, the Planck constraint on the Hubble constant is wider than in LCDM and in agreement with the Riess et al. (2016) direct measurement of $H_0$. The dark energy model is moderately favored as compared to LCDM when combining the KiDS and Planck measurements, and remains moderately favored after including an informative prior on the Hubble constant. In both of these scenarios, marginalized constraints in the $w_0-w_a$ plane are discrepant with a cosmological constant at the $3\sigma$ level. Moreover, KiDS constrains the sum of neutrino masses to 4.0 eV (95% CL), finds no preference for time or scale dependent modifications to the metric potentials, and is consistent with flatness and no running of the spectral index. The analysis code is public at https://github.com/sjoudaki/kids450
  • We constrain the average halo ellipticity of ~2 600 galaxy groups from the Galaxy And Mass Assembly (GAMA) survey, using the weak gravitational lensing signal measured from the overlapping Kilo Degree Survey (KiDS). To do so, we quantify the azimuthal dependence of the stacked lensing signal around seven different proxies for the orientation of the dark matter distribution, as it is a priori unknown which one traces the orientation best. On small scales, the major axis of the brightest group/cluster member (BCG) provides the best proxy, leading to a clear detection of an anisotropic signal. In order to relate that to a halo ellipticity, we have to adopt a model density profile. We derive new expressions for the quadrupole moments of the shear field given an elliptical model surface mass density profile. Modeling the signal with an elliptical Navarro-Frenk-White (NFW) profile on scales < 250 kpc, which roughly corresponds to half the virial radius, and assuming that the BCG is perfectly aligned with the dark matter, we find an average halo ellipticity of e_h=0.38 +/- 0.12. This agrees well with results from cold-dark-matter-only simulations, which typically report values of e_h ~ 0.3. On larger scales, the lensing signal around the BCGs does not trace the dark matter distribution well, and the distribution of group satellites provides a better proxy for the halo's orientation instead, leading to a 3--4 sigma detection of a non-zero halo ellipticity at scales between 250 kpc and 750 kpc. Our results suggest that the distribution of stars enclosed within a certain radius forms a good proxy for the orientation of the dark matter within that radius, which has also been observed in hydrodynamical simulations.
  • We study the stellar-to-halo mass relation of central galaxies in the range 9.7<log_10(M_*/h^-2 M_sun)<11.7 and z<0.4, obtained from a combined analysis of the Kilo Degree Survey (KiDS) and the Galaxy And Mass Assembly (GAMA) survey. We use ~100 deg^2 of KiDS data to study the lensing signal around galaxies for which spectroscopic redshifts and stellar masses were determined by GAMA. We show that lensing alone results in poor constraints on the stellar-to-halo mass relation due to a degeneracy between the satellite fraction and the halo mass, which is lifted when we simultaneously fit the stellar mass function. At M_sun>5x10^10 h^-2 M_sun, the stellar mass increases with halo mass as ~M_h^0.25. The ratio of dark matter to stellar mass has a minimum at a halo mass of 8x10^11 h^-1 M_sun with a value of M_h/M_*=56_-10^+16 [h]. We also use the GAMA group catalogue to select centrals and satellites in groups with five or more members, which trace regions in space where the local matter density is higher than average, and determine for the first time the stellar-to-halo mass relation in these denser environments. We find no significant differences compared to the relation from the full sample, which suggests that the stellar-to-halo mass relation does not vary strongly with local density. Furthermore, we find that the stellar-to-halo mass relation of central galaxies can also be obtained by modelling the lensing signal and stellar mass function of satellite galaxies only, which shows that the assumptions to model the satellite contribution in the halo model do not significantly bias the stellar-to-halo mass relation. Finally, we show that the combination of weak lensing with the stellar mass function can be used to test the purity of group catalogues.
  • We investigate possible signatures of halo assembly bias for spectroscopically selected galaxy groups from the GAMA survey using weak lensing measurements from the spatially overlapping regions of the deeper, high-imaging-quality photometric KiDS survey. We use GAMA groups with an apparent richness larger than 4 to identify samples with comparable mean host halo masses but with a different radial distribution of satellite galaxies, which is a proxy for the formation time of the haloes. We measure the weak lensing signal for groups with a steeper than average and with a shallower than average satellite distribution and find no sign of halo assembly bias, with the bias ratio of $0.85^{+0.37}_{-0.25}$, which is consistent with the $\Lambda$CDM prediction. Our galaxy groups have typical masses of $10^{13} M_{\odot}/h$, naturally complementing previous studies of halo assembly bias on galaxy cluster scales.
  • We measure the cross-correlation between Fermi-LAT gamma-ray photons and over 1000 deg$^2$ of weak lensing data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), the Red Cluster Sequence Lensing Survey (RCSLenS), and the Kilo Degree Survey (KiDS). We present the first measurement of tomographic weak lensing cross-correlations and the first application of spectral binning to cross-correlations between gamma rays and weak lensing. The measurements are performed using an angular power spectrum estimator while the covariance is estimated using an analytical prescription. We verify the accuracy of our covariance estimate by comparing it to two internal covariance estimators. Based on the non-detection of a cross-correlation signal, we derive constraints on weakly interacting massive particle (WIMP) dark matter. We compute exclusion limits on the dark matter annihilation cross-section $\langle\sigma_\rm{ann} v \rangle$, decay rate $\Gamma_\rm{dec}$, and particle mass $m_\rm{DM}$. We find that in the absence of a cross-correlation signal, tomography does not significantly improve the constraining power of the analysis. Assuming a strong contribution to the gamma-ray flux due to small-scale clustering of dark matter and accounting for known astrophysical sources of gamma rays, we exclude the thermal relic cross-section for particle masses of $m_\rm{DM}\lesssim 20$ GeV.
  • We present The-wiZZ, an open source and user-friendly software for estimating the redshift distributions of photometric galaxies with unknown redshifts by spatially cross-correlating them against a reference sample with known redshifts. The main benefit of The-wiZZ is in separating the angular pair finding and correlation estimation from the computation of the output clustering redshifts allowing anyone to create a clustering redshift for their sample without the intervention of an "expert". It allows the end user of a given survey to select any sub-sample of photometric galaxies with unknown redshifts, match this sample's catalog indices into a value-added data file, and produce a clustering redshift estimation for this sample in a fraction of the time it would take to run all the angular correlations needed to produce a clustering redshift. We show results with this software using photometric data from the Kilo-Degree Survey (KiDS) and spectroscopic redshifts from the Galaxy and Mass Assembly (GAMA) survey and the Sloan Digital Sky Survey (SDSS). The results we present for KiDS are consistent with the redshift distributions used in a recent cosmic shear analysis from the survey. We also present results using a hybrid machine learning-clustering redshift analysis that enables the estimation of clustering redshifts for individual galaxies. The-wiZZ can be downloaded at http://github.com/morriscb/The-wiZZ/.
  • We determine the accuracy of galaxy redshift distributions as estimated from photometric redshift probability distributions $p(z)$. Our method utilises measurements of the angular cross-correlation between photometric galaxies and an overlapping sample of galaxies with spectroscopic redshifts. We describe the redshift leakage from a galaxy photometric redshift bin $j$ into a spectroscopic redshift bin $i$ using the sum of the $p(z)$ for the galaxies residing in bin $j$. We can then predict the angular cross-correlation between photometric and spectroscopic galaxies due to intrinsic galaxy clustering when $i \neq j$ as a function of the measured angular cross-correlation when $i=j$. We also account for enhanced clustering arising from lensing magnification using a halo model. The comparison of this prediction with the measured signal provides a consistency check on the validity of using the summed $p(z)$ to determine galaxy redshift distributions in cosmological analyses, as advocated by the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). We present an analysis of the photometric redshifts measured by CFHTLenS, which overlaps the Baryon Oscillation Spectroscopic Survey (BOSS). We also analyse the Red-sequence Cluster Lensing Survey (RCSLenS), which overlaps both BOSS and the WiggleZ Dark Energy Survey. We find that the summed $p(z)$ from both surveys are generally biased with respect to the true underlying distributions. If unaccounted for, this bias would lead to errors in cosmological parameter estimation from CFHTLenS by less than $\sim 4\%$. For photometric redshift bins which spatially overlap in 3-D with our spectroscopic sample, we determine redshift bias corrections which can be used in future cosmological analyses that rely on accurate galaxy redshift distributions.