In this white paper, we discuss future uses of the LSST facility after the planned 10-year survey is complete. We expect the LSST survey to profoundly affect the scientific landscape over the next ten years, and it is likely that unexpected discoveries may drive its future scientific program. We discuss various operations and instrument options that could be considered for an extended LSST mission beyond ten years.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

The precise determination of the instrumental response function versus wavelength is a central ingredient in contemporary photometric calibration strategies. This typically entails propagating narrowband illumination through the system pupil, and comparing the detected photon rate across the focal plane to the amount of incident light as measured by a calibrated photodiode. However, stray light effects and reflections/ghosting (especially on the edges of filter passbands) in the optical train constitute a major source of systematic uncertainty when using a flat-field screen as the illumination source. A collimated beam projector that projects a mask onto the focal plane of the instrument can distinguish focusing light paths from stray and scattered light, allowing for a precise determination of instrumental throughput. This paper describes the conceptual design of such a system, outlines its merits, and presents results from a prototype system used with the Dark Energy Camera wide field imager on the 4-meter Blanco telescope. A calibration scheme that blends results from flat-field images with collimated beam projector data to obtain the equivalent of an illumination correction at high spectral and angular resolution is also presented. In addition to providing a precise system throughput calibration, by monitoring the evolution of the intensity and behaviour of the ghosts in the optical system, the collimated beam projector can be used to track the evolution of the filter transmission properties and various anti-reflective coatings in the optical system.

We use SDSS imaging data in $ugriz$ passbands to study the shape of the point spread function (PSF) profile and the variation of its width with wavelength and time. We find that the PSF profile is well described by theoretical predictions based on von K\'{a}rm\'{a}n's turbulence theory. The observed PSF radial profile can be parametrized by only two parameters, the profile's full width at half maximum (FWHM) and a normalization of the contribution of an empirically determined "instrumental" PSF. The profile shape is very similar to the "double gaussian plus power-law wing" decomposition used by SDSS image processing pipeline, but here it is successfully modeled with two free model parameters, rather than six as in SDSS pipeline. The FWHM variation with wavelength follows the $\lambda^{\alpha}$ power law, where $\alpha \approx-0.3$ and is correlated with the FWHM itself. The observed behavior is much better described by von K\'{a}rm\'{a}n's turbulence theory with the outer scale parameter in the range 5$-$100 m, than by the Kolmogorov's turbulence theory. We also measure the temporal and angular structure functions for FWHM and compare them to simulations and results from literature. The angular structure function saturates at scales beyond 0.5$-$1.0 degree. The power spectrum of the temporal behavior is found to be broadly consistent with a damped random walk model with characteristic timescale in the range $\sim5-30$ minutes, though data show a shallower high-frequency behavior. The latter is well fit by a single power law with index in the range $-1.5$ to $-1.0$. A hybrid model is likely needed to fully capture both the low-frequency and high-frequency behavior of the temporal variations of atmospheric seeing.

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

We report measurements from which we determine the spatial structure of the lunar contribution to night sky brightness, taken at the LSST site on Cerro Pachon in Chile. We use an array of six photodiodes with filters that approximate the Large Synoptic Survey Telescope's {\it u, g, r, i, z,} and {\it y} bands. We use the sun as a proxy for the moon, and measure sky brightness as a function of zenith angle of the point on sky, zenith angle of the sun, and angular distance between the sun and the point on sky. We make a correction for the difference between the illumination spectrum of the sun and the moon. Since scattered sunlight totally dominates the daytime sky brightness, this technique allows us to cleanly determine the contribution to the (cloudless) night sky from backscattered moonlight, without contamination from other sources of night sky brightness. We estimate our uncertainty in the relative lunar night sky brightness vs. zenith and lunar angle to be 10\,\%. This information is useful in planning the optimal execution of the LSST survey, and perhaps for other astronomical observations as well. Although our primary objective is to map out the angular structure and spectrum of the scattered light from the atmosphere and particulates, we also make an estimate of the expected number of scattered lunar photons per pixel per second in LSST, and find values that are in overall agreement with previous estimates.

The Large Synoptic Survey Telescope (LSST) will use an active optics system (AOS) to maintain alignment and surface figure on its three large mirrors. Corrective actions fed to the LSST AOS are determined from information derived from 4 curvature wavefront sensors located at the corners of the focal plane. Each wavefront sensor is a split detector such that the halves are 1mm on either side of focus. In this paper we describe the extensions to published curvature wavefront sensing algorithms needed to address challenges presented by the LSST, namely the large central obscuration, the fast f/1.23 beam, off-axis pupil distortions, and vignetting at the sensor locations. We also describe corrections needed for the split sensors and the effects from the angular separation of different stars providing the intra- and extra-focal images. Lastly, we present simulations that demonstrate convergence, linearity, and negligible noise when compared to atmospheric effects when the algorithm extensions are applied to the LSST optical system. The algorithm extensions reported here are generic and can easily be adapted to other wide-field optical systems including similar telescopes with large central obscuration and off-axis curvature sensing.

We present a comprehensive methodology for the simulation of astronomical images from optical survey telescopes. We use a photon Monte Carlo approach to construct images by sampling photons from models of astronomical source populations, and then simulating those photons through the system as they interact with the atmosphere, telescope, and camera. We demonstrate that all physical effects for optical light that determine the shapes, locations, and brightnesses of individual stars and galaxies can be accurately represented in this formalism. By using large scale grid computing, modern processors, and an efficient implementation that can produce 400,000 photons/second, we demonstrate that even very large optical surveys can be now be simulated. We demonstrate that we are able to: 1) construct kilometer scale phase screens necessary for wide-field telescopes, 2) reproduce atmospheric point-spread-function moments using a fast novel hybrid geometric/Fourier technique for non-diffraction limited telescopes, 3) accurately reproduce the expected spot diagrams for complex aspheric optical designs, and 4) recover system effective area predicted from analytic photometry integrals. This new code, the photon simulator (PhoSim), is publicly available. We have implemented the Large Synoptic Survey Telescope (LSST) design, and it can be extended to other telescopes. We expect that because of the comprehensive physics implemented in PhoSim, it will be used by the community to plan future observations, interpret detailed existing observations, and quantify systematics related to various astronomical measurements. Future development and validation by comparisons with real data will continue to improve the fidelity and usability of the code.

Early earthquake warning is a rapidly developing capability that has significant ramifications for many fields, including astronomical observatories. In this work, we describe the susceptibility of astronomical facilities to seismic events, including large telescopes as well as second-generation ground-based gravitational-wave interferometers. We describe the potential warning times for observatories from current seismic networks and propose locations for future seismometers to maximize warning times.

We describe a camera beam simulator for the LSST which is capable of illuminating a 60mm field at f/1.2 with realistic astronomical scenes, enabling studies of CCD astrometric and photometric performance. The goal is to fully simulate LSST observing, in order to characterize charge transport and other features in the thick fully depleted CCDs and to probe low level systematics under realistic conditions. The automated system simulates the centrally obscured LSST beam and sky scenes, including the spectral shape of the night sky. The doubly telecentric design uses a nearly unit magnification design consisting of a spherical mirror, three BK7 lenses, and one beam-splitter window. To achieve the relatively large field the beam-splitter window is used twice. The motivation for this LSST beam test facility was driven by the need to fully characterize a new generation of thick fully-depleted CCDs, and assess their suitability for the broad range of science which is planned for LSST. Due to the fast beam illumination and the thick silicon design [each pixel is 10 microns wide and over 100 microns deep] at long wavelengths there can be effects of photon transport and charge transport in the high purity silicon. The focal surface covers a field more than sufficient for a 40x40 mm LSST CCD. Delivered optical quality meets design goals, with 50% energy within a 5 micron circle. The tests of CCD performance are briefly described.

Near-future astronomical survey experiments, such as LSST, possess system requirements of unprecedented fidelity that span photometry, astrometry and shape transfer. Some of these requirements flow directly to the array of science imaging sensors at the focal plane. Availability of high quality characterization data acquired in the course of our sensor development program has given us an opportunity to develop and test a framework for simulation and modeling that is based on a limited set of physical and geometric effects. In this paper we describe those models, provide quantitative comparisons between data and modeled response, and extrapolate the response model to predict imaging array response to astronomical exposure. The emergent picture departs from the notion of a fixed, rectilinear grid that maps photo-conversions to the potential well of the channel. In place of that, we have a situation where structures from device fabrication, local silicon bulk resistivity variations and photo-converted carrier patterns still accumulating at the channel, together influence and distort positions within the photosensitive volume that map to pixel boundaries. Strategies for efficient extraction of modeling parameters from routinely acquired characterization data are described. Methods for high fidelity illumination/image distribution parameter retrieval, in the presence of such distortions, are also discussed.

The science goals for ground-based large-area surveys, such as the Dark Energy Survey, Pan-STARRS, and the Large Synoptic Survey Telescope, require calibration of broadband photometry that is stable in time and uniform over the sky to precisions of a per cent or better. This performance will need to be achieved with data taken over the course of many years, and often in less than ideal conditions. This paper describes a strategy to achieve precise internal calibration of imaging survey data taken in less than photometric conditions, and reports results of an observational study of the techniques needed to implement this strategy. We find that images of celestial fields used in this case study with stellar densities of order one per arcmin-squared and taken through cloudless skies can be calibrated with relative precision of 0.5 per cent (reproducibility). We report measurements of spatial structure functions of cloud absorption observed over a range of atmospheric conditions, and find it possible to achieve photometric measurements that are reproducible to 1 per cent in images that were taken through cloud layers that transmit as little as 25 per cent of the incident optical flux (1.5 magnitudes of extinction). We find, however, that photometric precision below 1 per cent is impeded by the thinnest detectable cloud layers. We comment on implications of these results for the observing strategies of future surveys.

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.

Cosmic shear holds great promise for a precision independent measurement of $\Omega\rm_m$, the mass density of the universe relative to the critical density. The signal is expected to be weak, so a thorough understanding of systematic effects is crucial. An important systematic effect is the atmosphere: shear power introduced by the atmosphere is larger than the expected signal. Algorithms exist to extract the cosmic shear from the atmospheric component, though a measure of their success applied to a range of seeing conditions is lacking. To gain insight into atmospheric shear, Gemini South imaging in conjunction with ground condition and satellite wind data were obtained. We find that under good seeing conditions Point-Spread-Function (PSF) correlations persist well beyond the separation typical of high-latitude stars. Under these conditions, ellipticity residuals based on a simple PSF interpolation can be reduced to within a factor of a few of the shot-noise induced ellipticity floor. We also find that the ellipticity residuals are highly correlated with wind direction. Finally, we correct stellar shapes using a more sophisticated procedure and generate shear statistics from stars. Under all seeing conditions in our data set the residual correlations lie everywhere below the target signal level. For good seeing we find that the systematic error attributable to atmospheric turbulence is comparable in magnitude to the statistical error (shape noise) over angular scales relevant to present lensing surveys.

We investigate the use of a narrow-band DDO51 filter for photometric identification of cool white dwarfs. We report photometric observations of 30 known cool white dwarfs with temperatures ranging from 10,000 K down to very cool temperatures (<3500 K). Follow-up spectroscopic observations of a sample of objects selected using this filter and our photometric observations show that DDO51 filter photometry can help select cool white dwarf candidates for follow-up multi--object spectroscopy by rejecting 65% of main sequence stars with the same broad--band colors as the cool white dwarfs. This technique is not selective enough to efficiently feed single--object spectrographs. We present the white dwarf cooling sequence using this filter. Our observations show that very cool white dwarfs form a sequence in the r-DDO vs. r-z color--color diagram and demonstrate that significant improvements are needed in white dwarf model atmospheres.

We report our search for the cool white dwarfs belonging to the Galactic disk by extending the NOAO Deep Wide-Field Survey. Narrow-band DDO51 photometry of the Deep Wide-Field Survey's northern field was obtained using the 4m-Mayall Telescope and the MOSAIC imager to separate cool white dwarfs from other stellar types of similar T_eff. Follow-up spectroscopy of four white dwarf candidates from our photometric search resulted in the discovery of two new cool white dwarfs as companions to M dwarfs.

SBS 1150+599A is a blue stellar object at high galactic latitude discovered in the Second Byurakan Survey. New high-resolution images of SBS 1150+599A are presented, demonstrating that it is very likely to be an old planetary nebula in the galactic halo, as suggested by Tovmassian et al (2001). An H-alpha image taken with the WIYN 3.5-m telescope and its "tip/tilt" module reveals the diameter of the nebula to be 9.2", comparable to that estimated from spectra by Tovmassian et al. Lower limits to the central star temperature were derived using the Zanstra hydrogen and helium methods to determine that the star's effective temperature must be > 68,000K and that the nebula is optically thin. New spectra from the MMT and FLWO telescopes are presented, revealing the presence of strong [Ne V] lambda 3425, indicating that the central star temperature must be > 100,000K. With the revised diameter, new central star temperature, and an improved central star luminosity, we can constrain photoionization models for the nebula significantly better than before. Because the emission-line data set is sparse, the models are still not conclusive. Nevertheless, we confirm that this nebula is an extremely metal-poor planetary nebula, having a value for O/H that is less than 1/100 solar, and possibly as low as 1/500 solar.

We present V, B-V, and DDO C(41-42) and C(42-45) photometry for a sample of 75 red giants down to M_V = +2 in the relatively metal-rich Galactic globular cluster M71. The C(41-42) colors reveal a bimodal distribution of CN band strengths generally anticorrelated with CH band strength as measured by the C(42-45) color. Both DDO colors agree well with those found in 47 Tucanae -- a nearby globular cluster of similar metallicity -- and suggest nearly identical C and N abundance patterns among the giants of both clusters. A comparison with synthetic DDO colors demonstrates that little change in surface C or N abundance is required to match the colors of the M71 giants over the entire luminosity range observed. Apparently like 47 Tuc (a cluster of much greater mass and central concentration), M71 exhibits an abundance pattern which cannot be solely the result of internal mixing.

The white dwarfs are promising laboratories for the study of cosmochronology and stellar evolution. Through observations of the pulsating white dwarfs, we can measure their internal structures and compositions, critical to understanding post main sequence evolution, along with their cooling rates, allowing us to calibrate their ages directly. The most important set of white dwarf variables to measure are the oldest of the pulsators, the cool DAVs, which have not previously been explored through asteroseismology due to their complexity and instability. Through a time-series photometry data set spanning ten years, we explore the pulsation spectrum of the cool DAV, G29-38 and find an underlying structure of 19 (not including multiplet components) normal-mode, probably l=1 pulsations amidst an abundance of time variability and linear combination modes. Modelling results are incomplete, but we suggest possible starting directions and discuss probable values for the stellar mass and hydrogen layer size. For the first time, we have made sense out of the complicated power spectra of a large-amplitude DA pulsator. We have shown its seemingly erratic set of observed frequencies can be understood in terms of a recurring set of normal-mode pulsations and their linear combinations. With this result, we have opened the interior secrets of the DAVs to future asteroseismological modelling, thereby joining the rest of the known white dwarf pulsators.

The F0 V star 9 Aur A exhibits an irregular variability of amplitude $\approx $0.1 magnitude at optical wavelengths. The variations are too slow for it to be a $\delta$ Scuti-type star. There is no evidence for a close, interacting companion or ring of dust, either from infrared, ultraviolet, or speckle data. The photometric variability of 9 Aur A is similar to two other early F dwarf stars: $\gamma$ Doradus and HD 96008. 9 Aur B appears to be an M dwarf, 9 Aur C is an early- to mid-K dwarf star, and 9 Aur E, if it is a member of the system, probably is a normal white dwarf. 9 Aur D is most likely an unrelated and distant K giant. [See note added in press regarding a fourth member of this class of ``variables without a cause", and short term variations of the radial velocity of 9 Aur.]