
A survey that can cover the sky in optical bands over wide fields to faint
magnitudes with a fast cadence will enable many of the exciting science
opportunities of the next decade. The Large Synoptic Survey Telescope (LSST)
will have an effective aperture of 6.7 meters and an imaging camera with field
of view of 9.6 deg^2, and will be devoted to a tenyear imaging survey over
20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with
fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a
total pointsource depth of r~27.5. The LSST Science Book describes the basic
parameters of the LSST hardware, software, and observing plans. The book
discusses educational and outreach opportunities, then goes on to describe a
broad range of science that LSST will revolutionize: mapping the inner and
outer Solar System, stellar populations in the Milky Way and nearby galaxies,
the structure of the Milky Way disk and halo and other objects in the Local
Volume, transient and variable objects both at low and high redshift, and the
properties of normal and active galaxies at low and high redshift. It then
turns to farfield cosmological topics, exploring properties of supernovae to
z~1, strong and weak lensing, the largescale distribution of galaxies and
baryon oscillations, and how these different probes may be combined to
constrain cosmological models and the physics of dark energy.

We present optical Rband light curves of five SDSS double QSOs (SDSS
J0903+5028, SDSS J1001+5027, SDSS J1206+4332, SDSS J1353+1138, SDSS J1335+0118)
obtained from monitoring at the Nordic Optical Telescope (NOT) between
September 2005 and September 2007. We also present analytical and pixelated
modeling of the observed systems. For SDSS J1206+4332, we measured the time
delay to be 116 days, which, for a Singular Isothermal Ellipsoid model,
corresponds to a Hubble constant of 73 km/s/Mpc. Simultaneous pixeleted
modeling of five other systems for which a time delay has now been previously
measured at the NOT leads to H_0 = 61.5 km/s/Mpc. Finally, by comparing
lightcurves of the two images of each system, suitably shifted by the predicted
or observed timedelays, we found no evidence for microlensing variability over
the course of the monitoring period.

In this paper, we describe a procedure for modelling strong lensing galaxy
clusters with parametric methods, and to rank models quantitatively using the
Bayesian evidence. We use a publicly available Markov chain MonteCarlo (MCMC)
sampler ('Bayesys'), allowing us to avoid local minima in the likelihood
functions. To illustrate the power of the MCMC technique, we simulate three
clusters of galaxies, each composed of a clusterscale halo and a set of
perturbing galaxyscale subhalos. We raytrace three light beams through each
model to produce a catalogue of multiple images, and then use the MCMC sampler
to recover the model parameters in the three different lensing configurations.
We find that, for typical Hubble Space Telescope (HST)quality imaging data,
the total mass in the Einstein radius is recovered with ~15% error according
to the considered lensing configuration. However, we find that the mass of the
galaxies is strongly degenerated with the cluster mass when no multiple images
appear in the cluster centre. The mass of the galaxies is generally recovered
with a 20% error, largely due to the poorly constrained cutoff radius.
Finally, we describe how to rank models quantitatively using the Bayesian
evidence. We confirm the ability of strong lensing to constrain the mass
profile in the central region of galaxy clusters in this way. Ultimately, such
a method applied to strong lensing clusters with a very large number of
multiple images may provide unique geometrical constraints on cosmology. The
implementation of the MCMC sampler used in this paper has been done within the
framework of the Lenstool software package, which is publicly available.

We present a parametric strong lensing model of the cluster Abell 2218 based
on HST ACS data. We constrain the lens model using 8 previously known multiply
imaged systems, of which 7 have spectroscopically confirmed redshifts. In
addition, we propose five candidate multiply imaged systems and estimate their
redshifts using our mass model. The model parameters are optimized in the
source plane by a bayesian Monte Carlo Markov Chain as implemented in the the
publicly available software Lenstool. We find rms=0."12 for the scatter of the
sources in the source plane, which translates into rms=1."49 between the
predicted and measured image positions in the image plane. We find that the
projected mass distribution of Abell 2218 is bimodal, which is supported by an
analysis of the light distribution. We also find evidence for two structures in
velocity space, separated by ~1000 km/s, corresponding to the two large scale
dark matter clumps. We find that the lensing constraints can not be well
reproduced using only dark matter halos associated with the cluster galaxies,
but that the dark matter is required to be smoothly distributed in large scale
halos. At 100" the enclosed projected mass is 3.8e14 solar masses. At that
radius, the large scale halos contribute ~85% of the mass. We find that the
model is not very sensitive to the fainter galaxy sized halos, unless they
locally perturb a given multiply imaged system. Therefore, dark galaxy sized
substructure can be reliably constrained only if it locally perturbs one of the
systems. In an appendix we give a selfcontained description of the parametric
profile we use, the dual pseudo isothermal elliptical mass distribution (dPIE).
This profile is a two component pseudo isothermal mass distribution (PIEMD)
with both a core radii and a scale radii. (Abridged)

We investigate the strong lensing properties of the Sersic profile as an
alternative to the NFW profile, focusing on applications to lens modelling of
clusters. Given an underlying Sersic dark matter profile, we study whether an
NFW profile can provide an acceptable fit to strong lensing constraints in the
form of single or multiple measured Einstein radii. We conclude that although
an NFW profile that fits the lensing constraints can be found in many cases,
the derived parameters may be biased. In particular, we find that for n~2,
which corresponds to massive clusters, the mass at r_200 of the best fit NFW is
overestimated (by a factor of ~2) and the concentration is very low (c~2). The
differences are important enough to warrant the inclusion of Sersic profile for
future analysis of strong lensing clusters.

We present a survey of the extinction properties of ten lensing galaxies, in
the redshift range z = 0.04  1.01, using multiply lensed quasars imaged with
the ESO VLT in the optical and near infrared. The multiple images act as
'standard light sources' shining through different parts of the lensing galaxy,
allowing for extinction studies by comparison of pairs of images. We explore
the effects of systematics in the extinction curve analysis, including
extinction along both lines of sight and microlensing, using theoretical
analysis and simulations. In the sample, we see variation in both the amount
and type of extinction. Of the ten systems, seven are consistent with
extinction along at least one line of sight. The mean differential extinction
for the most extinguished image pair for each lens is A(V) = 0.56 + 0.04,
using Galactic extinction law parametrization. The corresponding mean R_V = 2.8
+ 0.4 is consistent with that of the Milky Way at R_V = 3.1, where R_V =
A(V)/E(BV). We do not see any strong evidence for evolution of extinction
properties with redshift. Of the ten systems, B1152+199 shows the strongest
extinction signal of A(V) = 2.43 + 0.09 and is consistent with a Galactic
extinction law with R_V = 2.1 + 0.1. Given the similar redshift distribution
of SN Ia hosts and lensing galaxies, a large space based study of multiply
imaged quasars would be a useful complement to future dark energy SN Ia
surveys, providing independent constraints on the statistical extinction
properties of galaxies up to z~1.