
We investigate 3D density and weak lensing profiles of dark matter haloes
predicted by a cosmologyrescaling algorithm for $N$body simulations. We
extend the rescaling method of Angulo & White (2010) and Angulo & Hilbert
(2015) to improve its performance on intrahalo scales by using models for the
concentrationmassredshift relation based on excursion set theory. The
accuracy of the method is tested with numerical simulations carried out with
different cosmological parameters. We find that predictions for median density
profiles are more accurate than $\sim 5\,\%$ for haloes with masses of
$10^{12.0}  10^{14.5} h^{1}\,M_{\odot}$ for radii $0.05 < r/r_{200\text{m}} <
0.5$, and for cosmologies with $\Omega_\text{m} \in [0.15,\,0.40]$ and
$\sigma_8 \in [0.6,\,1.0]$. For larger radii, $0.5 < r/r_{200\text{m}} < 5$,
the accuracy degrades to $\sim20\,\%$, due to inaccurate modelling of the
cosmological and redshift dependence of the splashback radius. For changes in
cosmology allowed by current data, the residuals decrease to $\lesssim2\,\%$ up
to scales twice the virial radius. We illustrate the usefulness of the method
by estimating the mean halo mass of a mock galaxy group sample. We find that
the algorithm's accuracy is sufficient for current data. Improvements in the
algorithm, particularly in the modelling of baryons, are likely required for
interpreting future (dark energy task force stage IV) experiments.

Strong gravitational lenses with measured time delay are a powerful tool to
measure cosmological parameters, especially the Hubble constant ($H_0$). Recent
studies show that by combining just three multiplyimaged AGN systems, one can
determine $H_0$ to 3.8% precision. Furthermore, the number of timedelay lens
systems is growing rapidly, enabling, in principle, the determination of $H_0$
to 1% precision in the near future. However, as the precision increases it is
important to ensure that systematic errors and biases remain subdominant. For
this purpose, challenges with simulated datasets are a key component in this
process. Following the experience of the past challenge on time delay, where it
was shown that time delays can indeed be measured precisely and accurately at
the subpercent level, we now present the "Time Delay Lens Modeling Challenge"
(TDLMC). The goal of this challenge is to assess the present capabilities of
lens modeling codes and assumptions and test the level of accuracy of inferred
cosmological parameters given realistic mock datasets. We invite scientists to
model a set of simulated Hubble Space Telescope (HST) observations of 50 mock
lens systems. The systems are organized in rungs, with the complexity and
realism increasing going up the ladder. The goal of the challenge is to infer
$H_0$ for each rung, given the HST images, the time delay, and a stellar
velocity dispersion of the deflector, for a fixed background cosmology. The
TDLMC challenge will start with the mock data release on 2018 January 8th, with
a deadline for blind submission of 2018 August 8th. This first paper gives an
overview of the challenge including the data design, and a set of metrics to
quantify the modeling performance and challenge details. After the deadline,
the results of the challenge will be presented in a companion paper with all
challenge participants as coauthors.

Tidal gravitational forces can modify the shape of galaxies and clusters of
galaxies, thus correlating their orientation with the surrounding matter
density field. We study the dependence of this phenomenon, known as intrinsic
alignment (IA), on the mass of the dark matter haloes that host these bright
structures, analysing the Millennium and MillenniumXXL $N$body simulations.
We closely follow the observational approach, measuring the halo positionhalo
shape alignment and subsequently dividing out the dependence on halo bias. We
derive a theoretical scaling of the IA amplitude with mass in a dark matter
universe, and predict a powerlaw with slope $\beta_{\mathrm{M}}$ in the range
$1/3$ to $1/2$, depending on mass scale. We find that the simulation data agree
with each other and with the theoretical prediction remarkably well over three
orders of magnitude in mass, with the joint analysis yielding an estimate of
$\beta_{\mathrm{M}} = 0.36^{+0.01}_{0.01}$. This result does not depend on
redshift or on the details of the halo shape measurement. The analysis is
repeated on observational data, obtaining a significantly higher value,
$\beta_{\mathrm{M}} = 0.56^{+0.05}_{0.05}$. There are also small but
significant deviations from our simple model in the simulation signals at both
the high and lowmass end. We discuss possible reasons for these
discrepancies, and argue that they can be attributed to physical processes not
captured in the model or in the dark matteronly simulations.

We study intrinsic alignments (IA) of galaxy image shapes within the
Illustris cosmic structure formation simulations. We investigate how IA
correlations depend on observable galaxy properties such as stellar mass,
apparent magnitude, redshift, and photometric type, and on the employed shape
measurement method. The correlations considered include the matter
densityintrinsic ellipticity (mI), galaxy densityintrinsic ellipticity (dI),
gravitational shearintrinsic ellipticity (GI), and intrinsic
ellipticityintrinsic ellipticity (II) correlations. We find stronger
correlations for more massive and more luminous galaxies, as well as for
earlier photometric types, in agreement with observations. Moreover, the
correlations significantly depend on the choice of shape estimator, even if
calibrated to serve as unbiased shear estimators. In particular, shape
estimators that downweight the outer parts of galaxy images produce much
weaker IA signals on intermediate and large scales than methods employing flat
radial weights. The expected contribution of intrinsic alignments to the
observed ellipticity correlation in tomographic cosmic shear surveys may be
below one percent or several percent of the full signal depending on the
details of the shape measurement method. A comparison of our results to a tidal
alignment model indicates that such a model is able to reproduce the IA
correlations well on intermediate and large scales, provided the effect of
varying galaxy density is correctly taken into account. We also find that the
GI contributions to the observed ellipticity correlations could be inferred
directly from measurements of galaxy densityintrinsic ellipticity
correlations, except on small scales, where systematic differences between mI
and dI correlations are large.

Testing predictions of semianalytic models of galaxy evolution against
observations help to understand the complex processes that shape galaxies. We
compare predictions from the Garching and Durham models implemented on the
Millennium Run with observations of galaxygalaxy lensing (GGL) and
galaxygalaxygalaxy lensing (G3L) for various galaxy samples with stellar
masses in the range 0.5 < (M_* / 10^10 M_Sun) < 32 and photometric redshift
range 0.2 < z < 0.6 in the CanadaFranceHawaii Telescope Lensing Survey
(CFHTLenS). We find that the predicted GGL and G3L signals are in qualitative
agreement with CFHTLenS data. Quantitatively, the models succeed in reproducing
the observed signals in the highest stellar mass bin (16 < ( M_* / 10^10 M_Sun)
< 32) but show different degrees of tension for the other stellar mass samples.
The Durham models are strongly excluded at the 95% confidence level by the
observations as they largely overpredict the amplitudes of the GGL and G3L
signals, probably because they predict too many satellite galaxies in massive
halos.

We present a Bayesian reconstruction algorithm that infers the
threedimensional largescale matter distribution from the weak gravitational
lensing effects measured in the image shapes of galaxies. The algorithm assumes
that the prior probability distribution of the matter density is lognormal, in
contrast to many existing methods that assume normal (Gaussian) distributed
density fields. We compare the reconstruction results for both priors in a
suite of increasingly realistic tests on mock data. We find that in cases of
high noise levels (i.e. for low source galaxy densities and/or high shape
measurement uncertainties), both normal and lognormal priors lead to
reconstructions of comparable quality. In the lownoise regime, however, the
lognormal model produces significantly better reconstructions than the normal
model: The lognormal model 1) enforces nonnegative densities, while negative
densities are present when a normal prior is employed, 2) better traces the
extremal values and the skewness of the true underlying distribution, and 3)
yields a higher correlation between the reconstruction and the true density
distribution. Hence, the lognormal model is to be preferred over the normal
model, in particular since these differences become relevant for data from
current and futures surveys.

Strong gravitational lenses with measured time delays between the multiple
images allow a direct measurement of the timedelay distance to the lens, and
thus a measure of cosmological parameters, particularly the Hubble constant,
$H_{0}$. We present a blind lens model analysis of the quadruplyimaged quasar
lens HE 04351223 using deep Hubble Space Telescope imaging, updated timedelay
measurements from the COSmological MOnitoring of GRAvItational Lenses
(COSMOGRAIL), a measurement of the velocity dispersion of the lens galaxy based
on Keck data, and a characterization of the mass distribution along the line of
sight. HE 04351223 is the third lens analyzed as a part of the $H_{0}$ Lenses
in COSMOGRAIL's Wellspring (H0LiCOW) project. We account for various sources of
systematic uncertainty, including the detailed treatment of nearby perturbers,
the parameterization of the galaxy light and mass profile, and the regions used
for lens modeling. We constrain the effective timedelay distance to be
$D_{\Delta t} = 2612_{191}^{+208}~\mathrm{Mpc}$, a precision of 7.6%. From HE
04351223 alone, we infer a Hubble constant of $H_{0} =
73.1_{6.0}^{+5.7}~\mathrm{km~s^{1}~Mpc^{1}}$ assuming a flat $\Lambda$CDM
cosmology. The cosmographic inference based on the three lenses analyzed by
H0LiCOW to date is presented in a companion paper (H0LiCOW Paper V).

We present high signaltonoise galaxygalaxy lensing measurements of the
BOSS CMASS sample using 250 square degrees of weak lensing data from CFHTLenS
and CS82. We compare this signal with predictions from mock catalogs trained to
match observables including the stellar mass function and the projected and two
dimensional clustering of CMASS. We show that the clustering of CMASS, together
with standard models of the galaxyhalo connection, robustly predicts a lensing
signal that is 2040% larger than observed. Detailed tests show that our
results are robust to a variety of systematic effects. Lowering the value of
$S_{\rm 8}=\sigma_{\rm 8} \sqrt{\Omega_{\rm m}/0.3}$ compared to Planck2015
reconciles the lensing with clustering. However, given the scale of our
measurement ($r<10$ $h^{1}$ Mpc), other effects may also be at play and need
to be taken into consideration. We explore the impact of baryon physics,
assembly bias, massive neutrinos, and modifications to general relativity on
$\Delta\Sigma$ and show that several of these effects may be nonnegligible
given the precision of our measurement. Disentangling cosmological effects from
the details of the galaxyhalo connection, the effects of baryons, and massive
neutrinos, is the next challenge facing joint lensing and clustering analyses.
This is especially true in the context of large galaxy samples from Baryon
Acoustic Oscillation surveys with precise measurements but complex selection
functions.

Based on spectroscopy and multiband widefield observations of the
gravitationally lensed quasar HE 04351223, we determine the probability
distribution function of the external convergence $\kappa_\mathrm{ext}$ for
this system. We measure the under/overdensity of the line of sight towards the
lens system and compare it to the average line of sight throughout the
universe, determined by using the CFHTLenS as a control field. Aiming to
constrain $\kappa_\mathrm{ext}$ as tightly as possible, we determine
under/overdensities using various combinations of relevant informative weighing
schemes for the galaxy counts, such as projected distance to the lens,
redshift, and stellar mass. We then convert the measured under/overdensities
into a $\kappa_\mathrm{ext}$ distribution, using raytracing through the
Millennium Simulation. We explore several limiting magnitudes and apertures,
and account for systematic and statistical uncertainties relevant to the
quality of the observational data, which we further test through simulations.
Our most robust estimate of $\kappa_\mathrm{ext}$ has a median value
$\kappa^\mathrm{med}_\mathrm{ext} = 0.004$ and a standard deviation of
$\sigma_\kappa = 0.025$. The measured $\sigma_\kappa$ corresponds to $2.5\%$
uncertainty on the time delay distance, and hence the Hubble constant $H_0$
inference from this system. The median $\kappa^\mathrm{med}_\mathrm{ext}$ value
is robust to $\sim0.005$ (i.e. $\sim0.5\%$ on $H_0$) regardless of the adopted
aperture radius, limiting magnitude and weighting scheme, as long as the latter
incorporates galaxy number counts, the projected distance to the main lens, and
a prior on the external shear obtained from mass modeling. The availability of
a wellconstrained $\kappa_\mathrm{ext}$ makes \hequad\ a valuable system for
measuring cosmological parameters using strong gravitational lens time delays.

We present a tomographic cosmic shear study from the Deep Lens Survey (DLS),
which, providing a limiting magnitude r_{lim}~27 (5 sigma), is designed as a
precursor Large Synoptic Survey Telescope (LSST) survey with an emphasis on
depth. Using five tomographic redshift bins, we study their auto and
crosscorrelations to constrain cosmological parameters. We use a
luminositydependent nonlinear model to account for the astrophysical
systematics originating from intrinsic alignments of galaxy shapes. We find
that the cosmological leverage of the DLS is among the highest among existing
>10 sq. deg cosmic shear surveys. Combining the DLS tomography with the 9year
results of the Wilkinson Microwave Anisotropy Probe (WMAP9) gives
Omega_m=0.293_{0.014}^{+0.012}, sigma_8=0.833_{0.018}^{+0.011},
H_0=68.6_{1.2}^{+1.4} km/s/Mpc, and Omega_b=0.0475+0.0012 for LCDM, reducing
the uncertainties of the WMAP9only constraints by ~50%. When we do not assume
flatness for LCDM, we obtain the curvature constraint
Omega_k=0.010_{0.015}^{+0.013} from the DLS+WMAP9 combination, which however
is not well constrained when WMAP9 is used alone. The dark energy equation of
state parameter w is tightly constrained when Baryonic Acoustic Oscillation
(BAO) data are added, yielding w=1.02_{0.09}^{+0.10} with the DLS+WMAP9+BAO
joint probe. The addition of supernova constraints further tightens the
parameter to w=1.03+0.03. Our joint constraints are fully consistent with the
final Planck results and also the predictions of a LCDM universe.

We present weak lensing constraints on the ellipticity of galaxyscale matter
haloes and the galaxyhalo misalignment. Using data from the
CanadaFranceHawaii Telescope Lensing Survey (CFHTLenS), we measure the
weightedaverage ratio of the aligned projected ellipticity components of
galaxy matter haloes and their embedded galaxies, $f_\mathrm{h}$, split by
galaxy type. We then compare our observations to measurements taken from the
Millennium Simulation, assuming different models of galaxyhalo misalignment.
Using the Millennium Simulation we verify that the statistical estimator used
removes contamination from cosmic shear. We also detect an additional signal in
the simulation, which we interpret as the impact of intrinsic shapeshear
alignments between the lenses and their largescale structure environment.
These alignments are likely to have caused some of the previous observational
constraints on $f_\mathrm{h}$ to be biased high. From CFHTLenS we find
$f_\mathrm{h}=0.04 \pm 0.25$ for earlytype galaxies, which is consistent with
current models for the galaxyhalo misalignment predicting $f_\mathrm{h}\simeq
0.20$. For latetype galaxies we measure $f_\mathrm{h}=0.69_{0.36}^{+0.37}$
from CFHTLenS. This can be compared to the simulated results which yield
$f_\mathrm{h}\simeq 0.02$ for misaligned latetype models.

Depending on the exact experimental conditions, the thermodynamic properties
of physical systems can be related to one or more thermostatistical ensembles.
Here, we survey the notion of thermodynamic temperature in different
statistical ensembles, focusing in particular on subtleties that arise when
ensembles become nonequivalent. The 'mother' of all ensembles, the
microcanonical ensemble, uses entropy and internal energy (the most
fundamental, dynamically conserved quantity) to derive temperature as a
secondary thermodynamic variable. Over the past century, some confusion has
been caused by the fact that several competing microcanonical entropy
definitions are used in the literature, most commonly the volume and surface
entropies introduced by Gibbs. It can be proved, however, that only the volume
entropy satisfies exactly the traditional form of the laws of thermodynamics
for a broad class of physical systems, including all standard classical
Hamiltonian systems, regardless of their size. This mathematically rigorous
fact implies that negative 'absolute' temperatures and Carnot efficiencies $>1$
are not achievable within a standard thermodynamical framework. As an important
offspring of microcanonical thermostatistics, we shall briefly consider the
canonical ensemble and comment on the validity of the Boltzmann weight factor.
We conclude by addressing open mathematical problems that arise for systems
with discrete energy spectrum.

The recent experimental realization of exotic matter states in isolated
quantum systems and the ensuing controversy about the existence of negative
absolute temperatures demand a careful analysis of the conceptual foundations
underlying microcanonical thermostatistics. Here, we provide a detailed
comparison of the most commonly considered microcanonical entropy definitions,
focussing specifically on whether they satisfy or violate the zeroth, first and
second law of thermodynamics. Our analysis shows that, for a broad class of
systems that includes all standard classical Hamiltonian systems, only the
Gibbs volume entropy fulfills all three laws simultaneously. To avoid
ambiguities, the discussion is restricted to exact results and analytically
tractable examples.

We respond to invalid criticism in the recent Comment by Schneider et al.
[arXiv:1407.4127v1]

In their paper [arXiv:1403.4299v1], Frenkel and Warren claim that the Gibbs
temperature does not characterize thermal equilibrium correctly. We point out
the main logical errors in their argument.

The microphysical properties of the DM particle can, in principle, be
constrained by the properties and abundance of substructures in DM halos, as
measured through strong gravitational lensing. Unfortunately, there is a lack
of accurate theoretical predictions for the lensing signal of substructures,
mainly because of the discreteness noise inherent to Nbody simulations. Here
we present RecursiveTCM, a method that is able to provide lensing predictions
with an arbitrarily low discreteness noise, without any free parameters or
smoothing scale. This solution is based on a novel way of interpreting the
results of Nbody simulations, where particles simply trace the evolution and
distortion of Lagrangian phasespace volume elements. We discuss the advantages
of this method over the widely used cloudincells and adaptivekernel
smoothing density estimators. Applying the new method to a clustersized DM
halo simulated in warm and cold DM scenarios, we show how the expected
differences in their substructure population translate into differences in the
convergence and magnification maps. We anticipate that our method will provide
the highprecision theoretical predictions required to interpret and fully
exploit strong gravitational lensing observations.

A considerable body of experimental and theoretical work claims the existence
of negative absolute temperatures in spin systems and ultracold quantum gases.
Here, we clarify that such findings can be attributed to the use of a popular
yet inconsistent entropy definition, which violates fundamental thermodynamic
relations and fails to produce sensible results for simple analytically
tractable classical and quantum systems. Within a mathematically consistent
thermodynamic formalism, based on an entropy concept originally derived by
Gibbs, absolute temperature remains positive even for systems with bounded
spectrum. We address spurious arguments against the Gibbs formalism and comment
briefly on heat engines with efficiencies greater than one.

High precision cosmological distance measurements towards individual objects
such as time delay gravitational lenses or type Ia supernovae are affected by
weak lensing perturbations by galaxies and groups along the line of sight. In
time delay gravitational lenses, "external convergence," kappa, can dominate
the uncertainty in the inferred distances and hence cosmological parameters. In
this paper we attempt to reconstruct kappa, due to line of sight structure,
using a simple halo model. We use mock catalogues from the Millennium
Simulation, and calibrate and compare our reconstructed P(kappa) to raytraced
kappa "truth" values; taking into account realistic observational
uncertainties. We find that the reconstruction of kappa provides an improvement
in precision of ~50% over galaxy number counts. We find that the lowestkappa
lines of sight have the best constrained P(kappa). In anticipation of large
future samples of lenses, we find that selecting the third of the systems with
the highest precision kappa estimates gives a sample of unbiased time delay
distance measurements with just ~1% uncertainty due to line of sight external
convergence effects. Photometric data are sufficient to preselect the
bestconstrained lines of sight, and can be done before investment in
lightcurve monitoring. Conversely, we show that selecting lines of sight with
high external shear could, with the reconstruction model presented, induce
biases of up to 1% in time delay distance. We find that a major potential
source of systematic error is uncertainty in the high mass end of the stellar
masshalo mass relation; this could introduce ~2% biases on the timedelay
distance if completely ignored. We suggest areas for the improvement of this
general analysis framework (including more sophisticated treatment of high mass
structures) that should allow yet more accurate cosmological inferences to be
made.

In order to use strong gravitational lens time delays to measure precise and
accurate cosmological parameters the effects of mass along the line of sight
must be taken into account. We present a method to achieve this by constraining
the probability distribution function of the effective line of sight
convergence k_ext. The method is based on matching the observed overdensity in
the weighted number of galaxies to that found in mock catalogs with k_ext
obtained by raytracing through structure formation simulations. We explore
weighting schemes based on projected distance, mass, luminosity, and redshift.
This additional information reduces the uncertainty of k_ext from sigma_k $0.06
to ~0.04 for very overdense lines of sight like that of the system B1608+656.
For more common lines of sight, sigma_k is reduced to ~<0.03, corresponding to
an uncertainty of ~<3% on distance. This uncertainty has comparable effects on
cosmological parameters to that arising from the mass model of the deflector
and its immediate environment. Photometric redshifts based on g, r, i and K
photometries are sufficient to constrain k_ext almost as well as with
spectroscopic redshifts. As an illustration, we apply our method to the system
B1608+656. Our most reliable k_ext estimator gives sigma_k=0.047 down from
0.065 using only galaxy counts. Although deeper multiband observations of the
field of B1608+656 are necessary to obtain a more precise estimate, we conclude
that griK photometry, in addition to spectroscopy to characterize the immediate
environment, is an effective way to increase the precision of timedelay
cosmography.

We use weak gravitational lensing to analyse the dark matter halos around
satellite galaxies in galaxy groups in the CFHTLenS dataset. This dataset is
derived from the CFHTLSWide survey, and encompasses 154 sq. deg of
highquality shape data. Using the photometric redshifts, we divide the sample
of lens galaxies with stellar masses in the range 10^9 Msun to 10^10.5 Msun
into those likely to lie in highdensity environments (HDE) and those likely to
lie in lowdensity environments (LDE). Through comparison with galaxy
catalogues extracted from the Millennium Simulation, we show that the sample of
HDE galaxies should primarily (~61%) consist of satellite galaxies in groups,
while the sample of LDE galaxies should consist of mostly (~87%) nonsatellite
(field and central) galaxies. Comparing the lensing signals around samples of
HDE and LDE galaxies matched in stellar mass, the lensing signal around HDE
galaxies clearly shows a positive contribution from their host groups on their
lensing signals at radii of ~5001000 kpc, the typical separation between
satellites and group centres. More importantly, the subhalos of HDE galaxies
are less massive than those around LDE galaxies by a factor 0.65 +/ 0.12,
significant at the 2.9 sigma level. A natural explanation is that the halos of
satellite galaxies are stripped through tidal effects in the group environment.
Our results are consistent with a typical tidal truncation radius of ~40 kpc.

We present a cosmic shear study from the Deep Lens Survey (DLS), a deep BVRz
multiband imaging survey of five 4 sq. degree fields with two National Optical
Astronomy Observatory (NOAO) 4meter telescopes at Kitt Peak and Cerro Tololo.
For both telescopes, the change of the pointspreadfunction (PSF) shape across
the focal plane is complicated, and the exposuretoexposure variation of this
positiondependent PSF change is significant. We overcome this challenge by
modeling the PSF separately for individual exposures and CCDs with principal
component analysis (PCA). We find that stacking these PSFs reproduces the final
PSF pattern on the mosaic image with high fidelity, and the method successfully
separates PSFinduced systematics from gravitational lensing effects. We
calibrate our shears and estimate the errors, utilizing an image simulator,
which generates sheared groundbased galaxy images from deep Hubble Space
Telescope archival data with a realistic atmospheric turbulence model. For
cosmological parameter constraints, we marginalize over shear calibration
error, photometric redshift uncertainty, and the Hubble constant. We use
cosmologydependent covariances for the Markov Chain Monte Carlo analysis and
find that the role of this varying covariance is critical in our parameter
estimation. Our current nontomographic analysis alone constrains the
Omega_Msigma_8 likelihood contour tightly, providing a joint constraint of
Omega_M=0.262+0.051 and sigma_8=0.868+0.071. We expect that a future DLS
weaklensing tomographic study will further tighten these constraints because
explicit treatment of the redshift dependence of cosmic shear more efficiently
breaks the Omega_Msigma_8 degeneracy. Combining the current results with the
Wilkinson Microwave Anisotropy Probe 7year (WMAP7) likelihood data, we obtain
Omega_M=0.278+0.018 and sigma_8=0.815+0.020.

The gravitational lensing effect provides various ways to study the mass
environment of galaxies. We investigate how galaxygalaxy(galaxy) lensing can
be used to test models of galaxy formation and evolution. We consider two
semianalytic galaxy formation models based on the Millennium Run Nbody
simulation: the Durham model by Bower et al. (2006) and the Garching model by
Guo et al. (2011). We generate mock lensing observations for the two models,
and then employ Fast Fourier Transform methods to compute second and
thirdorder aperture statistics in the simulated fields for various galaxy
samples. We find that both models predict qualitatively similar aperture
signals, but there are large quantitative differences. The Durham model
predicts larger amplitudes in general. In both models, red galaxies exhibit
stronger aperture signals than blue galaxies. Using these aperture measurements
and assuming a linear deterministic bias model, we measure relative bias ratios
of red and blue galaxy samples. We find that a linear deterministic bias is
insufficient to describe the relative clustering of model galaxies below ten
arcmin angular scales. Dividing galaxies into luminosity bins, the aperture
signals decrease with decreasing luminosity for brighter galaxies, but increase
again for fainter galaxies. This increase is likely an artifact due to too many
faint satellite galaxies in massive group and cluster halos predicted by the
models. Our study shows that galaxygalaxy(galaxy) lensing is a sensitive
probe of galaxy evolution.

We demonstrate the possibility of detecting tidal stripping of dark matter
subhalos within galaxy groups using weak gravitational lensing. We have run
raytracing simulations on galaxy catalogues from the Millennium Simulation to
generate mock shape catalogues. The raytracing catalogues assume a halo model
for galaxies and groups, using various models with different distributions of
mass between galaxy and group halos to simulate different stages of group
evolution. Using these mock catalogues, we forecast the lensing signals that
will be detected around galaxy groups and satellite galaxies, as well as test
two different methods for isolating the satellites' lensing signals. A key
challenge is to determine the accuracy to which group centres can be
identified. We show that with current and ongoing surveys, it will possible to
detect stripping in groups of mass 10^1210^15 Msun.

We study the ability of future weak lensing (WL) surveys to constrain
primordial nonGaussianity of the local type. We use a large ensemble of
simulated WL maps with survey specifications relevant to Euclid and LSST. The
simulations assume Cold Dark Matter cosmologies that vary certain parameters
around fiducial values: the nonGaussianity parameter f_NL, the matter density
parameter Omega_m, the amplitude of the matter power spectrum sigma_8, the
spectral index of the primordial power spectrum n_s, and the darkenergy
equationofstate parameter w_0. We assess the sensitivity of the cosmic shear
correlation functions, the thirdorder aperture mass statistics, and the
abundance of shear peaks to these parameters. We find that each of the
considered probes provides unmarginalized constraints of Delta f_NL ~ 20 on
f_NL. Marginalized constraints from any individual WL probe are much weaker due
to strong correlations between parameters. However, the parameter errors can be
substantially reduced by combining information from different WL probes.
Combining all WL probes yields the following marginal (68% confidence level)
uncertainties: Delta f_NL ~ 50, Delta Omega_m ~ 0.002, Delta sigma_8 ~ 0.004,
Delta n_s ~ 0.007, and Delta w_0 ~ 0.03. We examine the bias induced by
neglecting f_NL on the constraints on the other parameters. We find sigma_8 and
w_0 to be the most affected. Moreover, neglecting nonGaussianity leads to a
severe underestimation of the uncertainties in the other cosmological
parameters.

We present a new method to extract cosmological constraints from weak lensing
(WL) peak counts, which we denote as `the hierarchical algorithm'. The idea of
this method is to combine information from WL maps sequentially smoothed with a
series of filters of different size, from the largest down to the smallest,
thus increasing the cosmological sensitivity of the resulting peak function. We
compare the cosmological constraints resulting from the peak abundance measured
in this way and the abundance obtained by using a filter of fixed size, which
is the standard practice in WL peak studies. For this purpose, we employ a
large set of WL maps generated by raytracing through Nbody simulations, and
the Fisher matrix formalism. We find that if lowS/N peaks are included in the
analysis (S/N ~ 3), the hierarchical method yields constraints significantly
better than the singlesized filtering. For a large future survey such as
Euclid or LSST, combined with information from a CMB experiment like Planck,
the results for the hierarchical (singlesized) method are: \Delta n=0.0039
(0.004); \Delta \Omega m=0.002 (0.0045); \Delta \sigma 8=0.003 (0.006); \Delta
w=0.019 (0.0525). This forecast is conservative, as we assume no knowledge of
the redshifts of the lenses, and consider a single broad bin for the redshifts
of the sources. If only peaks with S/N >= 6 are considered, then there is
little difference between the results of the two methods. We also examine the
statistical properties of the hierarchical peak function: Its covariance matrix
has offdiagonal terms for bins with S/N <= 6 and aperture mass of M < 3 x
1e+14 Ms/h, the higher bins being largely uncorrelated and therefore well
described by a Poisson distribution.