Large area surveys have detected significant samples of galaxy clusters that can be used to constrain cosmological parameters, provided that the masses of the clusters are measured robustly. To improve the calibration of cluster masses using weak gravitational lensing we present new results for 48 clusters at $0.05<z<0.15$, observed as part of the Multi Epoch Nearby Cluster Survey (MENeaCS), and reevaluate the mass estimates for 52 clusters from the Canadian Cluster Comparison Project (CCCP). Updated high-fidelity photometric redshift catalogues of reference deep fields are used in combination with advances in shape measurements and state-of-the-art cluster simulations, yielding an average systematic uncertainty in the lensing signal below 5%, similar to the statistical uncertainty for our cluster sample. We derive a scaling relation with Planck measurements for the full sample and find a bias in the Planck masses of $1-b=0.84 \pm 0.04$. We find no statistically significant trend of the mass bias with redshift or cluster mass, but find that different selections could change the bias by up to 1.5$\sigma$. We find a gas fraction of $0.139 \pm 0.014$ for 8 relaxed clusters in our sample, which can also be used to infer cosmological parameters.

Weak-lensing measurements of the masses of galaxy clusters are commonly based on the assumption of spherically symmetric density profiles. Yet, the cold dark matter model predicts the shapes of dark matter halos to be triaxial. Halo triaxiality, and the orientation of the major axis with respect to the line of sight, are expected to be the leading cause of intrinsic scatter in weak-lensing mass measurements. The shape of central cluster galaxies (Brightest Cluster Galaxies; BCGs) is expected to follow the shape of the dark matter halo. Here we investigate the use of BCG ellipticity as predictor of the weak-lensing mass bias in individual clusters compared to the mean. Using weak lensing masses $M^{\rm WL}_{500}$ from the Weighing the Giants project, and $M_{500}$ derived from gas masses as low-scatter mass proxy, we find that, on average, the lensing masses of clusters with the roundest / most elliptical 25% of BCGs are biased $\sim 20$% high / low compared to the average, as qualitatively predicted by the cold dark matter model. For cluster cosmology projects utilizing weak-lensing mass estimates, the shape of the BCG can thus contribute useful information on the effect of orientation bias in weak lensing mass estimates as well as on cluster selection bias.

In modern weak-lensing surveys, the common approach to correct for residual systematic biases in the shear is to calibrate shape measurement algorithms using simulations. These simulations must fully capture the complexity of the observations to avoid introducing any additional bias. In this paper we study the importance of faint galaxies below the observational detection limit of a survey. We simulate simplified Euclid VIS images including and excluding this faint population, and measure the shift in the multiplicative shear bias between the two sets of simulations. We measure the shear with three different algorithms: a moment-based approach, model fitting, and machine learning. We find that for all methods, a spatially uniform random distribution of faint galaxies introduces a shear multiplicative bias of the order of a few times $10^{-3}$. This value increases to the order of $10^{-2}$ when including the clustering of the faint galaxies, as measured in the Hubble Space Telescope Ultra-Deep Field. The magnification of the faint background galaxies due to the brighter galaxies along the line of sight is found to have a negligible impact on the multiplicative bias. We conclude that the undetected galaxies must be included in the calibration simulations with proper clustering properties down to magnitude 28 in order to reach a residual uncertainty on the multiplicative shear bias calibration of a few times $10^{-4}$, in line with the $2\times10^{-3}$ total accuracy budget required by the scientific objectives of the Euclid survey. We propose two complementary methods for including faint galaxy clustering in the calibration simulations.

Exploiting the full statistical power of future cosmic shear surveys will necessitate improvements to the accuracy with which the gravitational lensing signal is measured. We present a framework for calibrating shear with image simulations that demonstrates the importance of including realistic correlations between galaxy morphology, size and more importantly, photometric redshifts. This realism is essential so that selection and shape measurement biases can be calibrated accurately for a tomographic cosmic shear analysis. We emulate Kilo-Degree Survey (KiDS) observations of the COSMOS field using morphological information from {\it Hubble} Space Telescope imaging, faithfully reproducing the measured galaxy properties from KiDS observations of the same field. We calibrate our shear measurements from lensfit, and find through a range of sensitivity tests that lensfit is robust and unbiased within the allowed 2 per cent tolerance of our study. Our results show that the calibration has to be performed by selecting the tomographic samples in the simulations, consistent with the actual cosmic shear analysis, because the joint distributions of galaxy properties are found to vary with redshift. Ignoring this redshift variation could result in misestimating the shear bias by an amount that exceeds the allowed tolerance. To improve the calibration for future cosmic shear analyses, it will be essential to also correctly account for the measurement of photometric redshifts, which requires simulating multi-band observations.

We measure the gravitational lensing signal around satellite galaxies in a sample of galaxy clusters at $z<0.15$ by combining high-quality imaging data from the Canada-France-Hawaii Telescope with a large sample of spectroscopically-confirmed cluster members. We use extensive image simulations to assess the accuracy of shape measurements of faint, background sources in the vicinity of bright satellite galaxies. We find a small but significant bias, as light from the lenses makes the shapes of background galaxies appear radially aligned with the lens. We account for this bias by applying a correction that depends on both lens size and magnitude. We also correct for contamination of the source sample by cluster members. We use a physically-motivated definition of subhalo mass, namely the mass bound to the subhalo, $m_\mathrm{bg}$, similar to definitions used by common subhalo finders in numerical simulations. Binning the satellites by stellar mass we provide a direct measurement of the subhalo-to-stellar-mass relation, $\log m_\mathrm{bg}/\mathrm{M}_\odot = (11.54\pm0.05) + (0.95\pm0.10)\log[m_\star/(2\times10^{10}\mathrm{M}_\odot)]$. This best-fitting relation implies that, at a stellar mass $m_\star\sim3\times10^{10}\,\mathrm{M}_\odot$, subhalo masses are roughly 50 per cent of those of central galaxies, and this fraction decreases at higher stellar masses. We find some evidence for a sharp change in the total-to-stellar mass ratio around the clusters' scale radius, which could be interpreted as galaxies within the scale radius having suffered more strongly from tidal stripping, but remain cautious regarding this interpretation.

One-opposition near-Earth asteroids (NEAs) are growing in number, and they must be recovered to prevent loss and mismatch risk, and to improve their orbits, as they are likely to be too faint for detection in shallow surveys at future apparitions. We aimed to recover more than half of the one-opposition NEAs recommended for observations by the Minor Planet Center (MPC) using the Isaac Newton Telescope (INT) in soft-override mode and some fractions of available D-nights. During about 130 hours in total between 2013 and 2016, we targeted 368 NEAs, among which 56 potentially hazardous asteroids (PHAs), observing 437 INT Wide Field Camera (WFC) fields and recovering 280 NEAs (76% of all targets). Engaging a core team of about ten students and amateurs, we used the THELI, Astrometrica, and the Find_Orb software to identify all moving objects using the blink and track-and-stack method for the faintest targets and plotting the positional uncertainty ellipse from NEODyS. Most targets and recovered objects had apparent magnitudes centered around V~22.8 mag, with some becoming as faint as V~24 mag. One hundred and three objects (representing 28% of all targets) were recovered by EURONEAR alone by Aug 2017. Orbital arcs were prolonged typically from a few weeks to a few years; our oldest recoveries reach 16 years. The O-C residuals for our 1,854 NEA astrometric positions show that most measurements cluster closely around the origin. In addition to the recovered NEAs, 22,000 positions of about 3,500 known minor planets and another 10,000 observations of about 1,500 unknown objects (mostly main-belt objects) were promptly reported to the MPC by our team. Four new NEAs were discovered serendipitously in the analyzed fields, increasing the counting to nine NEAs discovered by the EURONEAR in 2014 and 2015.

The recent discovery of thousands of ultra diffuse galaxies (UDGs) in nearby galaxy clusters has opened a new window into the process of galaxy formation and evolution. Several scenarios have been proposed to explain the formation history of UDGs, and their ability to survive in the harsh cluster environments. A key requirement to distinguish between these scenarios is a measurement of their halo masses which, due to their low surface brightnesses, has proven difficult if one relies on stellar tracers of the potential. We exploit weak gravitational lensing, a technique that does not depend on these baryonic tracers, to measure the average subhalo mass of 784 UDGs selected in 18 clusters at $z\leq0.09$. Our sample of UDGs has a median stellar mass $\langle m_\star\rangle=2\times10^8\,\mathrm{M}_\odot$ and a median effective radius $\langle r_\mathrm{eff}\rangle=2.8$ kpc. We constrain the average mass of subhaloes within 30 kpc to $\log m_\mathrm{UDG}(r<30\,\mathrm{kpc})/\mathrm{M}_\odot\leq10.99$ at 95 per cent credibility, implying an effective virial mass $\log m_{200}/\mathrm{M}_\odot\leq11.80$, and a lower limit on the stellar mass fraction within 10 kpc of 1.0 per cent. Such mass is consistent with a simple extrapolation of the subhalo-to-stellar mass relation of typical satellite galaxies in massive clusters. However, our analysis is not sensitive to scatter about this mean mass; the possibility remains that extreme UDGs reside in haloes as massive as the Milky Way.

The Kilo-Degree Survey (KiDS) is an ongoing optical wide-field imaging survey with the OmegaCAM camera at the VLT Survey Telescope. It aims to image 1500 square degrees in four filters (ugri). The core science driver is mapping the large-scale matter distribution in the Universe, using weak lensing shear and photometric redshift measurements. Further science cases include galaxy evolution, Milky Way structure, detection of high-redshift clusters, and finding rare sources such as strong lenses and quasars. Here we present the third public data release (DR3) and several associated data products, adding further area, homogenized photometric calibration, photometric redshifts and weak lensing shear measurements to the first two releases. A dedicated pipeline embedded in the Astro-WISE information system is used for the production of the main release. Modifications with respect to earlier releases are described in detail. Photometric redshifts have been derived using both Bayesian template fitting, and machine-learning techniques. For the weak lensing measurements, optimized procedures based on the THELI data reduction and lensfit shear measurement packages are used. In DR3 stacked ugri images, weight maps, masks, and source lists for 292 new survey tiles (~300 sq.deg) are made available. The multi-band catalogue, including homogenized photometry and photometric redshifts, covers the combined DR1, DR2 and DR3 footprint of 440 survey tiles (447 sq.deg). Limiting magnitudes are typically 24.3, 25.1, 24.9, 23.8 (5 sigma in a 2 arcsec aperture) in ugri, respectively, and the typical r-band PSF size is less than 0.7 arcsec. The photometric homogenization scheme ensures accurate colors and an absolute calibration stable to ~2% for gri and ~3% in u. Separately released are a weak lensing shear catalogue and photometric redshifts based on two different machine-learning techniques.

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

Improvements in the accuracy of shape measurements are essential to exploit the statistical power of planned imaging surveys that aim to constrain cosmological parameters using weak lensing by large-scale structure. Although a range of tests can be performed using the measurements, the performance of the algorithm can only be quantified using simulated images. This yields, however, only meaningful results if the simulated images resemble the real observations sufficiently well. In this paper we explore the sensitivity of the multiplicative bias to the input parameters of Euclid-like image simulations.We find that algorithms will need to account for the local density of sources. In particular the impact of galaxies below the detection limit warrants further study, because magnification changes their number density, resulting in correlations between the lensing signal and multiplicative bias. Although achieving sub-percent accuracy will require further study, we estimate that sufficient archival Hubble Space Telescope data are available to create realistic populations of galaxies.

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 describe and test the pipeline used to measure the weak lensing shear signal from the Kilo Degree Survey (KiDS). It includes a novel method of `self-calibration' that partially corrects for the effect of noise bias. We also discuss the `weight bias' that may arise in optimally-weighted measurements, and present a scheme to mitigate that bias. To study the residual biases arising from both galaxy selection and shear measurement, and to derive an empirical correction to reduce the shear biases to $\lesssim 1\%$, we create a suite of simulated images whose properties are close to those of the KiDS survey observations. We find that the use of `self-calibration' reduces the additive and multiplicative shear biases significantly, although further correction via a calibration scheme is required, which also corrects for a dependence of the bias on galaxy properties. We find that the calibration relation itself is biased by the use of noisy, measured galaxy properties, which may limit the final accuracy that can be achieved. We assess the accuracy of the calibration in the tomographic bins used for the KiDS cosmic shear analysis, testing in particular the effect of possible variations in the uncertain distributions of galaxy size, magnitude and ellipticity, and conclude that the calibration procedure is accurate at the level of multiplicative bias $\lesssim 1\%$ required for the KiDS cosmic shear analysis.

Current optical imaging surveys for cosmology are covering large areas of sky. To exploit the statistical power of these surveys for weak lensing measurements requires shape measurement methods with subpercent systematic errors. We introduce a new weak lensing shear measurement algorithm, Shear Nulling after PSF Gaussianisation (SNAPG), designed to avoid the noise biases that affect most other methods. SNAPG operates on images that have been convolved with a kernel that renders the Point Spread Function (PSF) a circular Gaussian, and uses weighted second moments of the sources. The response of such second moments to a shear of the pre-seeing galaxy image can be predicted analytically, allowing us to construct a shear nulling scheme that finds the shear parameters for which the observed galaxies are consistent with an unsheared, isotropically oriented population of sources. The inverse of this nulling shear is then an estimate of the gravitational lensing shear. We identify the uncertainty of the estimated centre of each galaxy as the source of noise bias, and incorporate an approximate estimate of the centroid covariance into the scheme. We test the method on extensive suites of simulated galaxies of increasing complexity, and find that it is capable of shear measurements with multiplicative bias below 0.5 percent.

We present cosmological parameter constraints from a tomographic weak gravitational lensing analysis of ~450deg$^2$ of imaging data from the Kilo Degree Survey (KiDS). For a flat $\Lambda$CDM cosmology with a prior on $H_0$ that encompasses the most recent direct measurements, we find $S_8\equiv\sigma_8\sqrt{\Omega_{\rm m}/0.3}=0.745\pm0.039$. This result is in good agreement with other low redshift probes of large scale structure, including recent cosmic shear results, along with pre-Planck cosmic microwave background constraints. A $2.3$-$\sigma$ tension in $S_8$ and `substantial discordance' in the full parameter space is found with respect to the Planck 2015 results. We use shear measurements for nearly 15 million galaxies, determined with a new improved `self-calibrating' version of $lens$fit validated using an extensive suite of image simulations. Four-band $ugri$ photometric redshifts are calibrated directly with deep spectroscopic surveys. The redshift calibration is confirmed using two independent techniques based on angular cross-correlations and the properties of the photometric redshift probability distributions. Our covariance matrix is determined using an analytical approach, verified numerically with large mock galaxy catalogues. We account for uncertainties in the modelling of intrinsic galaxy alignments and the impact of baryon feedback on the shape of the non-linear matter power spectrum, in addition to the small residual uncertainties in the shear and redshift calibration. The cosmology analysis was performed blind. Our high-level data products, including shear correlation functions, covariance matrices, redshift distributions, and Monte Carlo Markov Chains are available at http://kids.strw.leidenuniv.nl.

The Planck catalogue of SZ sources limits itself to a significance threshold of 4.5 to ensure a low contamination rate by false cluster candidates. This means that only the most massive clusters at redshift z>0.5, and in particular z>0.7, are expected to enter into the catalogue, with a large number of systems in that redshift regime being expected around and just below that threshold. In this paper, we follow-up a sample of SZ sources from the Planck SZ catalogues from 2013 and 2015. In the latter maps, we consider detections around and at lower significance than the threshold adopted by the Planck Collaboration. To keep the contamination rate low, our 28 candidates are chosen to have significant WISE detections, in combination with non-detections in SDSS/DSS, which effectively selects galaxy cluster candidates at redshifts $z\gtrsim0.5$. By taking r- and z-band imaging with MegaCam@CFHT, we bridge the 4000A rest-frame break over a significant redshift range, thus allowing accurate redshift estimates of red-sequence cluster galaxies up to z~0.8. After discussing the possibility that an overdensity of galaxies coincides -by chance- with a Planck SZ detection, we confirm that 16 of the candidates have likely optical counterparts to their SZ signals, 13 (6) of which have an estimated redshift z>0.5 (z>0.7). The richnesses of these systems are generally lower than expected given the halo masses estimated from the Planck maps. However, when we follow a simplistic model to correct for Eddington bias in the SZ halo mass proxy, the richnesses are consistent with a reference mass-richness relation established for clusters detected at higher significance. This illustrates the benefit of an optical follow-up, not only to obtain redshift estimates, but also to provide an independent mass proxy that is not based on the same data the clusters are detected with, and thus not subject to Eddington bias.

The Kilo-Degree Survey (KiDS) is a multi-band imaging survey designed for cosmological studies from weak lensing and photometric redshifts. It uses the ESO VLT Survey Telescope with its wide-field camera OmegaCAM. KiDS images are taken in four filters similar to the SDSS ugri bands. The best-seeing time is reserved for deep r-band observations that reach a median 5-sigma limiting AB magnitude of 24.9 with a median seeing that is better than 0.7arcsec. Initial KiDS observations have concentrated on the GAMA regions near the celestial equator, where extensive, highly complete redshift catalogues are available. A total of 109 survey tiles, one square degree each, form the basis of the first set of lensing analyses, which focus on measurements of halo properties of GAMA galaxies. 9 galaxies per square arcminute enter the lensing analysis, for an effective inverse shear variance of 69 per square arcminute. Accounting for the shape measurement weight, the median redshift of the sources is 0.53. KiDS data processing follows two parallel tracks, one optimized for galaxy shape measurement (for weak lensing), and one for accurate matched-aperture photometry in four bands (for photometric redshifts). This technical paper describes how the lensing and photometric redshift catalogues have been produced (including an extensive description of the Gaussian Aperture and Photometry pipeline), summarizes the data quality, and presents extensive tests for systematic errors that might affect the lensing analyses. We also provide first demonstrations of the suitability of the data for cosmological measurements, and explain how the shear catalogues were blinded to prevent confirmation bias in the scientific analyses. The KiDS shear and photometric redshift catalogues, presented in this paper, are released to the community through http://kids.strw.leidenuniv.nl .

We study the evolution of the luminosity-to-halo 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) SDSS-DR10 in the ~450 deg^2 that overlaps with imaging data from the second Red-sequence 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 mass-concentration 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 luminosity-to-halo 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 25_{-14}^{+16} % between z=0.36 and 0.22 to an average value of 6.43+/-0.52 x 10^13 h70^-1 Msun. 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 80^{+39}_{-28} % between z=0.6 and 0.2

Masses of clusters of galaxies from weak gravitational lensing analyses of ever larger samples are increasingly used as the reference to which baryonic scaling relations are compared. In this paper we revisit the analysis of a sample of 50 clusters studied as part of the Canadian Cluster Comparison Project. We examine the key sources of systematic error in cluster masses. We quantify the robustness of our shape measurements and calibrate our algorithm empirically using extensive image simulations. The source redshift distribution is revised using the latest state-of-the-art photometric redshift catalogs that include new deep near-infrared observations. Nonetheless we find that the uncertainty in the determination of photometric redshifts is the largest source of systematic error for our mass estimates. We use our updated masses to determine b, the bias in the hydrostatic mass, for the clusters detected by Planck. Our results suggest 1-b=0.76+-0.05(stat)}+-0.06(syst)}, which does not resolve the tension with the measurements from the primary cosmic microwave background.

This paper presents the first results of a far-infrared search for protocluster-associated galaxy overdensities using the SPIRE instrument on-board the {\it Herschel} Space Observatory. Large ($\sim$400 arcmin$^{2}$) fields surrounding 26 powerful high-redshift radio galaxies ($2.0 < z < 4.1$; $L_{\rm 500 MHz} > 10^{28.5}$ WHz$^{-1}$) are mapped at 250, 350 and 500\mic to give a unique wide-field sample. On average the fields have a higher than expected, compared to blank fields, surface density of 500\mic sources within 6 comoving Mpc of the radio galaxy. The analysis is then restricted to potential protocluster members only, which are identified using a far-infrared colour selection; this reveals significant overdensities of galaxies in 2 fields, neither of which are previously known protoclusters. The probability of finding 2 overdensities of this size by chance, given the number of fields observed is $5 \times 10^{-4}$. Overdensities here exist around radio galaxies with $L_{\rm 500 MHz} \gtrsim 10^{29}$ WHz$^{-1}$ and $z < 3$. The radial extent of the average far-infrared overdensity is found to be $\sim$6 comoving Mpc. Comparison with predictions from numerical simulations shows that the overdensities are consistent with having masses $> 10^{14}$Msolar. However, the large uncertainty in the redshift estimation means that it is possible that these far-infrared overdensities consist of several structures across the redshift range searched.