
We measure the gravitational lensing signal around satellite galaxies in a
sample of galaxy clusters at $z<0.15$ by combining highquality imaging data
from the CanadaFranceHawaii Telescope with a large sample of
spectroscopicallyconfirmed 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
physicallymotivated 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 subhalotostellarmass 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 bestfitting
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
totaltostellar 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.

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 KiloDegree 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 signaltonoise 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 volumelimited 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 highredshift 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.

Accurate shape measurements are essential to infer cosmological parameters
from large area weak gravitational lensing studies. The compact
diffractionlimited pointspread function (PSF) in spacebased observations is
greatly beneficial, but its chromaticity for a broad band observation can lead
to new subtle effects that could hitherto be ignored: the PSF of a galaxy is no
longer uniquely defined and spatial variations in the colours of galaxies
result in biases in the inferred lensing signal. Taking Euclid as a reference,
we show that this colourgradient bias (CG bias) can be quantified with high
accuracy using available multicolour Hubble Space Telescope (HST) data. In
particular we study how noise in the HST observations might impact such
measurements and find this to be negligible. We determine the CG bias using HST
observations in the F606W and F814W filters and observe a correlation with the
colour, in line with expectations, whereas the dependence with redshift is
weak. The biases for individual galaxies are generally well below 1%, which may
be reduced further using morphological information from the Euclid data. Our
results demonstrate that CG bias should not be ignored, but it is possible to
determine its amplitude with sufficient precision, so that it will not
significantly bias the weak lensing measurements using Euclid data.

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 galaxymatter crosscorrelation
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 realspace 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
KiDS450 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 largescale structure
probe combination.

We demonstrate that deep goodseeing VLT/HAWKI $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 singleorbit 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
highredshift galaxy samples. Incorporating results from the 3DHST 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{HAWKI}
$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 simulationbased priors on the
concentration, making it one of the most massive galaxy clusters known in the
\mbox{$z\gtrsim 0.7$} Universe. We also crosschecked the HAWKI measurements
through an analysis of overlapping HST/ACS images, yielding fully consistent
estimates of the lensing signal.

We measure the projected galaxy clustering and galaxygalaxy lensing signals
using the Galaxy And Mass Assembly (GAMA) survey and KiloDegree Survey (KiDS)
to study galaxy bias. We use the concept of nonlinear 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 nonlinearity and
stochasticity in halo occupation models can be inferred. Our observations
unveil the physical reason for the nonlinearity and stochasticity, further
confirmed using hydrodynamical simulations, with the stochasticity mostly
originating from the nonPoissonian 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 nonlinearity 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.

Observations of relaxed, massive and distant clusters can provide important
tests of standard cosmological models e.g. using the gas mass fraction. We
study the very luminous, high redshift ($z=0.902$) galaxy cluster
ClJ120958.9+495352 using XMMNewton data and measure the temperature profile
and cooling time to investigate the dynamical status with respect to the
presence of a cool core as well as global cluster properties. We use HST weak
lensing data to estimate its total mass and determine the gas mass fraction. We
perform a spectral analysis using an XMMNewton observation of 15ks cleaned
exposure time. As the treatment of the background is crucial, we use two
different approaches to account for the background emission to verify our
results. We account for pointspreadfunction effects and deproject our results
to estimate the gas mass fraction of the cluster. We measure weak lensing
galaxy shapes from mosaic HST imaging and select background galaxies
photometrically in combination with imaging data from the William Herschel
Telescope. The Xray luminosity of ClJ120958.9+495352 in the 0.12.4keV band
estimated from our XMMNewton data is $L_X =
(13.4_{1.0}^{+1.2})\times10^{44}$erg/s and thus it is one of the most Xray
luminous clusters known at similarly high redshift. We find clear indications
for the presence of a cool core from the temperature profile and the central
cooling time, which is very rare at such high redshifts. Based on the weak
lensing analysis we estimate a cluster mass of
$M_\mathrm{500}/10^{14}M_\odot=4.4^{+2.2}_{2.0}(\mathrm{stat.})\pm0.6(\mathrm{sys.})$
and a gas mass fraction of $f_\mathrm{gas,2500} = 0.11_{0.03}^{+0.06}$ in good
agreement with previous findings for high redshift and local clusters.

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 subhalotostellar 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.

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 (KiDS450). We measure high
signaltonoise 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
massconcentration 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 highSNR 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 NFWlike
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 signaltonoise 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 shearposition 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 singlemassive halos. This
demonstrates the gain of information from lowS/N peaks which correspond to the
projection of several smallmass halos along the lineofsight. Our results are
in good agreement with the tomographic shear twopoint correlation function
measurement in KiDS450. Combining shear peaks with nontomographic
measurements of the shear twopoint correlation functions yields an $\sim20\%$
improvement in the uncertainty on $S_8$ compared to the shear twopoint
correlation functions alone, highlighting the great potential of peaks as a
cosmological probe.

In recent years, multiple studies have reported substantial populations of
large, lowsurfacebrightness galaxies in local galaxy clusters. Various
theories that aim to explain the presence of such ultradiffuse galaxies (UDGs)
have since been proposed. A key question that will help to differentiate
between models is whether UDGs have counterparts in lowermass 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 spectroscopicallyselected groups from the Galaxy And
Mass Assembly (GAMA) survey. We make use of the overlapping imaging from the
ESO KiloDegree 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 rband. 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 massrichness 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.

We perform a combined analysis of cosmic shear tomography, galaxygalaxy
lensing tomography, and redshiftspace 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 2degree Field Lensing Survey
(2dFLenS) and the Baryon Oscillation Spectroscopic Survey (BOSS). We restrict
the galaxygalaxy lensing and multipole power spectrum measurements to the
overlapping regions with KiDS, and selfconsistently 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 galaxygalaxy 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.

The convolution of galaxy images by the pointspread function (PSF) is the
dominant source of bias for weak gravitational lensing studies, and an accurate
estimate of the PSF is required to obtain unbiased shape measurements. The PSF
estimate for a galaxy depends on its spectral energy distribution (SED),
because the instrumental PSF is generally a function of the wavelength. In this
paper we explore various approaches to determine the resulting `effective' PSF
using broadband data. Considering the Euclid mission as a reference, we find
that standard SED template fitting methods result in biases that depend on
source redshift, although this may be remedied if the algorithms can be
optimised for this purpose. Using a machinelearning algorithm we show that, at
least in principle, the required accuracy can be achieved with the current
survey parameters. It is also possible to account for the correlations between
photometric redshift and PSF estimates that arise from the use of the same
photometry. We explore the impact of errors in photometric calibration, errors
in the assumed wavelength dependence of the PSF model and limitations of the
adopted template libraries. Our results indicate that the required accuracy for
Euclid can be achieved using the data that are planned to determine photometric
redshifts.

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 timedependent equation of state (in the form of
the $w_0w_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_0w_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
NavarroFrenkWhite (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 colddarkmatteronly 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 34 sigma detection of a nonzero 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 stellartohalo 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 stellartohalo 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 stellartohalo mass
relation in these denser environments. We find no significant differences
compared to the relation from the full sample, which suggests that the
stellartohalo mass relation does not vary strongly with local density.
Furthermore, we find that the stellartohalo 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
stellartohalo 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.

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 largescale 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 Euclidlike 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 subpercent 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,
highimagingquality 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 present the tomographic crosscorrelation between galaxy lensing measured
in the Kilo Degree Survey (KiDS450) with overlapping lensing measurements of
the cosmic microwave background (CMB), as detected by Planck 2015. We compare
our joint probe measurement to the theoretical expectation for a flat
$\Lambda$CDM cosmology, assuming the bestfitting cosmological parameters from
the KiDS450 cosmic shear and Planck CMB analyses. We find that our results are
consistent within $1\sigma$ with the KiDS450 cosmology, with an amplitude
rescaling parameter $A_{\rm KiDS} = 0.86 \pm 0.19$. Adopting a Planck
cosmology, we find our results are consistent within $2\sigma$, with $A_{\it
Planck} = 0.68 \pm 0.15$. We show that the agreement is improved in both cases
when the contamination to the signal by intrinsic galaxy alignments is
accounted for, increasing $A$ by $\sim 0.1$. This is the first tomographic
analysis of the galaxy lensing  CMB lensing crosscorrelation signal, and is
based on five photometric redshift bins. We use this measurement as an
independent validation of the multiplicative shear calibration and of the
calibrated source redshift distribution at high redshifts. We find that
constraints on these two quantities are strongly correlated when obtained from
this technique, which should therefore not be considered as a standalone
competitive calibration tool.

We measure the crosscorrelation between FermiLAT gammaray photons and over
1000 deg$^2$ of weak lensing data from the CanadaFranceHawaii 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 crosscorrelations and the first application of
spectral binning to crosscorrelations 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 nondetection of a crosscorrelation signal, we derive
constraints on weakly interacting massive particle (WIMP) dark matter. We
compute exclusion limits on the dark matter annihilation crosssection
$\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 crosscorrelation signal,
tomography does not significantly improve the constraining power of the
analysis. Assuming a strong contribution to the gammaray flux due to
smallscale clustering of dark matter and accounting for known astrophysical
sources of gamma rays, we exclude the thermal relic crosssection for particle
masses of $m_\rm{DM}\lesssim 20$ GeV.

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
`selfcalibration' that partially corrects for the effect of noise bias. We
also discuss the `weight bias' that may arise in optimallyweighted
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 `selfcalibration' 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.

Verlinde (2016) proposed that the observed excess gravity in galaxies and
clusters is the consequence of Emergent Gravity (EG). In this theory the
standard gravitational laws are modified on galactic and larger scales due to
the displacement of dark energy by baryonic matter. EG gives an estimate of the
excess gravity (described as an apparent dark matter density) in terms of the
baryonic mass distribution and the Hubble parameter. In this work we present
the first test of EG using weak gravitational lensing, within the regime of
validity of the current model. Although there is no direct description of
lensing and cosmology in EG yet, we can make a reasonable estimate of the
expected lensing signal of low redshift galaxies by assuming a background
LambdaCDM cosmology. We measure the (apparent) average surface mass density
profiles of 33,613 isolated central galaxies, and compare them to those
predicted by EG based on the galaxies' baryonic masses. To this end we employ
the ~180 square degrees overlap of the KiloDegree Survey (KiDS) with the
spectroscopic Galaxy And Mass Assembly (GAMA) survey. We find that the
prediction from EG, despite requiring no free parameters, is in good agreement
with the observed galaxygalaxy lensing profiles in four different stellar mass
bins. Although this performance is remarkable, this study is only a first step.
Further advancements on both the theoretical framework and observational tests
of EG are needed before it can be considered a fully developed and solidly
tested theory.

We present predictions for the galaxygalaxy lensing profile from the EAGLE
hydrodynamical cosmological simulation at redshift z=0.18, in the spatial range
0.02 < R/(Mpc/h) < 2, and for five logarithmically equispaced stellar mass
bins in the range 10.3 < $\log_{10}$(Mstar/ $M_{\odot}$) < 11.8. We compare
these excess surface density profiles to the observed signal from background
galaxies imaged by the Kilo Degree Survey around spectroscopically confirmed
foreground galaxies from the GAMA survey. Exploiting the GAMA galaxy group
catalogue, the profiles of central and satellite galaxies are computed
separately for groups with at least five members to minimise contamination.
EAGLE predictions are in broad agreement with the observed profiles for both
central and satellite galaxies, although the signal is underestimated at
R$\approx$0.52 Mpc/h for the highest stellar mass bins. When central and
satellite galaxies are considered simultaneously, agreement is found only when
the selection function of lens galaxies is taken into account in detail.
Specifically, in the case of GAMA galaxies, it is crucial to account for the
variation of the fraction of satellite galaxies in bins of stellar mass induced
by the fluxlimited nature of the survey. We report the inferred
stellartohalo mass relation and we find good agreement with recent published
results. We note how the precision of the galaxygalaxy lensing profiles in the
simulation holds the potential to constrain finegrained aspects of the
galaxydark matter connection.

Euclid is a European Space Agency medium class mission selected for launch in
2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to
understand the origin of the accelerated expansion of the universe. Euclid will
explore the expansion history of the universe and the evolution of cosmic
structures by measuring shapes and redshifts of galaxies as well as the
distribution of clusters of galaxies over a large fraction of the sky. Although
the main driver for Euclid is the nature of dark energy, Euclid science covers
a vast range of topics, from cosmology to galaxy evolution to planetary
research. In this review we focus on cosmology and fundamental physics, with a
strong emphasis on science beyond the current standard models. We discuss five
broad topics: dark energy and modified gravity, dark matter, initial
conditions, basic assumptions and questions of methodology in the data
analysis. This review has been planned and carried out within Euclid's Theory
Working Group and is meant to provide a guide to the scientific themes that
will underlie the activity of the group during the preparation of the Euclid
mission.

Using a sample of 98 galaxy clusters recently imaged in the near infrared
with the ESO NTT, WIYN and WHT telescopes, supplemented with 33 clusters from
the ESO archive, we measure how the stellar mass of the most massive galaxies
in the universe, namely Brightest Cluster Galaxies (BCG), increases with time.
Most of the BCGs in this new sample lie in the redshift range $0.2<z<0.6$,
which has been noted in recent works to mark an epoch over which the growth in
the stellar mass of BCGs stalls. From this sample of 132 clusters, we create a
subsample of 102 systems that includes only those clusters that have estimates
of the cluster mass. We combine the BCGs in this subsample with BCGs from the
literature, and find that the growth in stellar mass of BCGs from 10 billion
years ago to the present epoch is broadly consistent with recent semianalytic
and semiempirical models. As in other recent studies, tentative evidence
indicates that the stellar mass growth rate of BCGs may be slowing in the past
3.5 billion years. Further work in collecting larger samples, and in better
comparing observations with theory using mock images is required if a more
detailed comparison between the models and the data is to be made.