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 ten-year 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 point-source 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 far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale 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 R-band 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 time-delays, 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 Monte-Carlo (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 cluster-scale halo and a set of perturbing galaxy-scale subhalos. We ray-trace 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 ~1-5% 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 cut-off 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 self-contained 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(B-V). 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.