The LSST survey will provide unprecedented statistical power for measurements of dark energy. Consequently, controlling systematic uncertainties is becoming more important than ever. The LSST observing strategy will affect the statistical uncertainty and systematics control for many science cases; here, we focus on weak lensing systematics. The fact that the LSST observing strategy involves hundreds of visits to the same sky area provides new opportunities for systematics mitigation. We explore these opportunities by testing how different dithering strategies (pointing offsets and rotational angle of the camera in different exposures) affect additive weak lensing shear systematics on a baseline operational simulation, using the $\rho-$statistics formalism. Some dithering strategies improve systematics control at the end of the survey by a factor of up to $\sim 3-4$ better than others. We find that a random translational dithering strategy, applied with random rotational dithering at every filter change, is the most effective of those strategies tested in this work at averaging down systematics. Adopting this dithering algorithm, we explore the effect of varying the area of the survey footprint, exposure time, number of exposures in a visit, and exposure to the Galactic plane. We find that any change that increases the average number of exposures (in filters relevant to weak lensing) reduces the additive shear systematics. Some ways to achieve this increase may not be favorable for the weak lensing statistical constraining power or for other probes, and we explore the relative trade-offs between these options given constraints on the overall survey parameters.

We report studies on mitigation of optical effects of bright Low-Earth-Orbit (LEO) satellites on Vera C. Rubin Observatory and its Legacy Survey of Space and Time (LSST). These include options for pointing the telescope to avoid satellites, laboratory investigations of bright trails on the Rubin Observatory LSST Camera sensors, algorithms for correcting image artifacts caused by bright trails, experiments on darkening SpaceX Starlink satellites, and ground-based follow-up observations. Starlink satellites with no brightness mitigation are presently g ~ 5.1 mag, and an initial experiment "DarkSat" is g ~ 6.1 mag. Future Starlink darkening plans may reach g ~ 7 mag, a brightness level that enables non-linear image artifact correction to the same level as background noise. However, the satellite trails will still exist at S/N ~ 100, generating systematics that may impact data analysis and limiting some science. LEO satellite trails are wider than a point-spread function because satellites are slightly out of focus due to their finite distance; for Rubin Observatory's 8.4-m mirror and a satellite at 550 km, this widening effect is about 3 arcsec, which helps avoid saturation by decreasing the trail's peak surface brightness.

We test the impact of some systematic errors in weak lensing magnification measurements with the COSMOS 30-band photo-$z$ Survey flux limited to $i_{AB}<25.0$ using correlations of both source galaxy counts and fluxes. Systematic obscuration effects are measured by comparing counts and flux correlations. We use the ACS-HST catalogs to identify potential blending objects (close pairs) and perform the magnification analyses with and without blended objects. We find that blending effects start to be important ($\sim$ 0.04~mag obscuration) at angular scales smaller than 0.1 arcmin. Extinction and other systematic obscuration effects can be as large as 0.10~mag (U-band) but are typically smaller than 0.02~mag depending on the band. After applying these corrections, we measure a $3.9\sigma$ magnification signal that is consistent for both counts and flux magnification. The corresponding projected mass profiles of galaxies at redshift $z \simeq 0.6$ ($M_I \simeq -21$) is $\Sigma= 25\pm 6 M_{sun}h^3/pc^2$ at 0.1 Mpc/h, consistent with NFW type profile with $M_{200} \simeq 2 \times 10^{12} M_{sun} h/pc^2$. Tangential shear and flux-size magnification over the same lenses show similar mass profiles. We conclude that magnification from counts and fluxes using photometric redshifts has the potential to provide complementary weak lensing information in future wide field surveys once we carefully take into account systematic effects, such as obscuration and blending.

A dedicated simulator, Poisson_CCD, has been constructed which models astronomical CCDs by solving Poisson's equation numerically and simulating charge transport within the CCD. The potentials and free carrier densities within the CCD are self-consistently solved for, giving realistic results for the charge distribution within the CCD storage wells. The simulator has been used to model the CCDs which are being used to construct the LSST digital camera. The simulator output has been validated by comparing its predictions with several different types of CCD measurements, including astrometric shifts, brighter-fatter induced pixel-pixel covariances, saturation effects, and diffusion spreading. The code is open source and freely available.

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.

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

Cosmology is one of the four science pillars of LSST, which promises to be transformative for our understanding of dark energy and dark matter. The LSST Dark Energy Science Collaboration (DESC) has been tasked with deriving constraints on cosmological parameters from LSST data. Each of the cosmological probes for LSST is heavily impacted by the choice of observing strategy. This white paper is written by the LSST DESC Observing Strategy Task Force (OSTF), which represents the entire collaboration, and aims to make recommendations on observing strategy that will benefit all cosmological analyses with LSST. It is accompanied by the DESC DDF (Deep Drilling Fields) white paper (Scolnic et al.). We use a variety of metrics to understand the effects of the observing strategy on measurements of weak lensing, large-scale structure, clusters, photometric redshifts, supernovae, strong lensing and kilonovae. In order to reduce systematic uncertainties, we conclude that the current baseline observing strategy needs to be significantly modified to result in the best possible cosmological constraints. We provide some key recommendations: moving the WFD (Wide-Fast-Deep) footprint to avoid regions of high extinction, taking visit pairs in different filters, changing the 2x15s snaps to a single exposure to improve efficiency, focusing on strategies that reduce long gaps (>15 days) between observations, and prioritizing spatial uniformity at several intervals during the 10-year survey.

We present cosmological parameter measurements from the Deep Lens Survey (DLS) using galaxy-mass and galaxy-galaxy power spectra in the multipole range $\ell=250\sim2000$. We measure galaxy-galaxy power spectra from two lens bins centered at $z\sim0.27$ and $0.54$ and galaxy-mass power spectra by cross-correlating the positions of galaxies in these two lens bins with galaxy shapes in two source bins centered at $z\sim0.64$ and $1.1$. We marginalize over a baryonic feedback process using a single-parameter representation and a sum of neutrino masses, as well as photometric redshift and shear calibration systematic uncertainties. For a flat $\Lambda$CDM cosmology, we determine $S_8\equiv\sigma_8\sqrt{\Omega_m/0.3}=0.810^{+0.039}_{-0.031}$, in good agreement with our previous DLS cosmic shear and the Planck Cosmic Microwave Background (CMB) measurements. Without the baryonic feedback marginalization, $S_8$ decreases by $\sim0.05$ because the dark matter-only power spectrum lacks the suppression at the highest $\ell$'s due to Active Galactic Nuclei (AGN) feedback. Together with the Planck CMB measurement, we constrain the baryonic feedback parameter to $A_{baryon}=1.07^{+0.31}_{-0.39}$, which suggests an interesting possibility that the actual AGN feedback might be stronger than the recipe used in the OWLS simulations. The interpretation is limited by the validity of the baryonic feedback simulation and the one-parameter representation of the effect.

We explore the weak lensing E- and B-mode shear signals of a field of galaxy clusters using both large scale structure N-body simulations and multi-color Suprime-cam & Hubble Space Telescope observations. Using the ray-traced and observed shears along with photometric redshift catalogs, we generate mass maps of the foreground overdensities by optimally filtering the tangential shear that they induce on background galaxies. We then develop and test a method to approximate the foreground structure as a superposition of NFW-like halos by locating these overdensities and determining their mass and redshift, thereby modeling the background correlated shear field as a sum of lensings induced by the foreground clusters. We demonstrate that the B-mode maps and shear correlation functions, which are generated by similarly filtering the cross shear in this method, are in agreement with observations and are related to the estimated cluster masses and locations as well as the distribution of background sources. Using the foreground mass model, we identify several sources of weak lensing B-modes including leakage and edge effects, source clustering, and multiple lensing which can be observed in deep cosmic shear surveys.

The 3.2 gigapixel LSST camera, an array of 189 thick fully-depleted CCDs, will repeatedly image the southern sky and accomplish a wide variety of science goals. However, its trove of tens of billions of object images implies stringent requirements on systematic biases imprinted during shift-and-stare CCD observation. In order to correct for these biases which, without correction, violate requirements on weak lensing precision, we investigate CCD systematics using both simulations of charge transport as well as with a unique bench-top optical system matched to the LSST's fast f/1.2 beam. By illuminating single CCDs with realistic scenes of stars and galaxies and then analyzing these images with the LSST data management pipelines, we can characterize the survey's imaging performance well before the camera's first light. We present measurements of several CCD systematics under varying conditions in the laboratory, including the brightness-dependent broadening of star and galaxy images, charge transport anomalies in the silicon bulk as well as the edges, and serial deferred charge. Alongside these measurements, we also present the development and testing of physics-based models which inform corrections or mitigation strategies for these systematics. Optimization of the CCD survey operation under a variety of realistic observational conditions, including systematic effects from the optics, clocking, sky brightness, and image analysis, will be critical to achieve the LSST's goals of precision astronomy and cosmology.

We develop an improved sky background estimator which employs optimal filters for both spatial and pixel intensity distributions. It incorporates growth of masks around detected objects and a statistical estimate of the flux from undetected faint galaxies in the remaining sky pixels. We test this algorithm for underlying sky estimation and compare its performance with commonly used sky estimation codes on realistic simulations which include detected galaxies, faint undetected galaxies, and sky noise. We then test galaxy surface brightness recovery using GALFIT 3, a galaxy surface brightness profile fitting optimizer, yielding fits to S\'{e}rsic profiles. This enables robust sky background estimates accurate at the 4 parts-per-million level. This background sky estimator is more accurate and is less affected by surface brightness profiles of galaxies and the local image environment compared with other methods.

(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 Large Synoptic Survey Telescope (LSST) is a high etendue imaging facility that is being constructed atop Cerro Pachon in Northern Chile. It is scheduled to begin science operations in 2022. With an 8.4m (6.5m effective) aperture, a three-mirror design achieving a seeing-limited 9.6deg^2 field of view, and a 3.2 Gigapixel camera, the LSST has the deep-wide-fast imaging capability to carry out an 18,000deg^2 survey in six passbands (ugrizy) to a coadded depth of r~27.5 over 10 years using 90% of its observational time. The remaining 10% will be devoted to deeper and faster time-domain observations and smaller surveys. In total, each patch of the sky in the main survey will receive 800 visits allocated across the six passbands with 30s exposure visits. The huge volume of high-quality LSST data will provide a wide range of science opportunities and open a new era of precision cosmology with unprecedented statistical power and tight control of systematic errors. In this review, we give a brief account of the LSST cosmology program with an emphasis on dark energy investigations. The LSST will address dark energy physics and cosmology in general by exploiting diverse precision probes including large-scale structure, weak lensing, type Ia supernovae, galaxy clusters, and strong lensing. Combined with the cosmic microwave background data, these probes form interlocking tests on the cosmological model and the nature of dark energy. The LSST data products will be made available to the U.S. and Chilean scientific communities and to international partners with no proprietary period. Close collaborations with coeval surveys observing at a variety of wavelengths, resolutions, depths, and timescales will be a vital part of the LSST science program, which will not only enhance specific studies but also allow a more complete understanding of the universe through different windows.

The Large Synoptic Survey Telescope (LSST) will enable revolutionary studies of galaxies, dark matter, and black holes over cosmic time. The LSST Galaxies Science Collaboration has identified a host of preparatory research tasks required to leverage fully the LSST dataset for extragalactic science beyond the study of dark energy. This Galaxies Science Roadmap provides a brief introduction to critical extragalactic science to be conducted ahead of LSST operations, and a detailed list of preparatory science tasks including the motivation, activities, and deliverables associated with each. The Galaxies Science Roadmap will serve as a guiding document for researchers interested in conducting extragalactic science in anticipation of the forthcoming LSST era.

Reduction of images and science analysis from ground-based telescopes such as the LSST requires detailed knowledge of the PSF of the image, which includes components attributable to the instrument as well as components attributable to the atmosphere. Because the atmospheric component is constantly changing, the PSF is typically extracted from each image by measuring the size and shape of star images across the CCD, then building a fitting function over the focal plane which is used to model the PSF for analysis of extended sources such as galaxies. Since the stars in each CCD field have a range of brightnesses, accurate knowledge of the PSF for point sources of varying brightness is essential. It has been found that in thick, fully-depleted CCDs, the electrostatic repulsion of charge stored in the collecting wells gives rise to a larger and slightly more elliptical PSF for brighter stars. This "brighter-fatter" effect has been reported in some detail in the literature. In this work, we report direct and indirect measurements of this effect in prototype LSST sensors, and describe a detailed physics-based model of the electrostatics and charge transport within the CCD.

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 pre-cursor Large Synoptic Survey Telescope (LSST) survey with an emphasis on depth. Using five tomographic redshift bins, we study their auto- and cross-correlations to constrain cosmological parameters. We use a luminosity-dependent 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 9-year 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 WMAP9-only 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.

Using overlapping fields with space-based Hubble Space Telescope (HST) and ground-based Subaru Telescope imaging we identify a population of blended galaxies that are blended to such a large degree that they are detected as single objects in the ground-based monochromatic imaging, which we label as 'ambiguous blends'. For deep imaging data, such as the depth targeted with the Large Synoptic Survey Telescope (LSST), the ambiguous blend population is both large ($\sim 14$%) and has a distribution of ellipticities that is different from that of unblended objects in a way that will likely be important for the weak lensing measurements. Most notably, for a limiting magnitude of $i \sim 27$ we find that ambiguous blending results in a ~14% increase in shear noise (or ~12% decrease in the effective projected number density of lensed galaxies; neff) due to 1) larger intrinsic ellipticity dispersion, 2) a scaling with the galaxy number density $N_{gal}$ that is shallower than 1/$\sqrt{N_{gal}}$. For the LSST Gold Sample ($i < 25.3$) there is a ~7% increase in shear noise (or ~7% decrease in $n_{eff}$). More importantly than these increases in the shear noise, we find that the ellipticity distribution of ambiguous blends has an RMS 13% larger than that of non-blended galaxies. Given the need of future weak lensing surveys to constrain the ellipticity distribution of galaxies to better than a percent in order to mitigate cosmic shear multiplicative biases, the different ellipticity distribution of ambiguous blends could be a dominant systematic if unaccounted for.

We characterize the astrometric distortion at the edges of thick, fully-depleted CCDs in the lab using a bench-top simulation of LSST observing. By illuminating an array of forty thousand pinholes (30mu m diameter) at the object plane of a f/1.2 optical reimager, thousands of PSFs can be imaged over a 4Kx4K pixel CCD. Each high purity silicon pixel, 10mu m square by 100mu m deep, can then be individually characterized through a series of sub-pixel dithers in the X/Y plane. The unique character [response, position, shape] of each pixel as a function of flux, wavelength, back side bias, etc. can be investigated. We measure the magnitude and onset of astrometric error at the edges of the detector as a test of the experimental setup, using a LSST prototype CCD. We show that this astrometric error at the edge is sourced from non-uniformities in the electric field lines that define pixel boundaries. This edge distortion must be corrected in order to optimize the science output of weak gravitational lensing and large scale structure measurements for the LSST.

Thick fully depleted CCDs, while enabling wide spectral response, also present challenges in understanding the systematic errors due to 3D charge transport. This 2014 Workshop on Precision Astronomy with Fully Depleted CCDs covered progress that has been made in the testing and modeling of these devices made since a workshop by the same name in 2013. Presentations covered the science drivers, CCD characterization, laboratory measurements of systematics, calibration, and different approaches to modeling the response and charge transport. The key issue is the impact of these CCD sensor features on dark energy science, including astrometry and photometry. Successful modeling of the spatial systematics can enable first order correction in the data processing pipeline.

The focus of this report is on the opportunities enabled by the combination of LSST, Euclid and WFIRST, the optical surveys that will be an essential part of the next decade's astronomy. The sum of these surveys has the potential to be significantly greater than the contributions of the individual parts. As is detailed in this report, the combination of these surveys should give us multi-wavelength high-resolution images of galaxies and broadband data covering much of the stellar energy spectrum. These stellar and galactic data have the potential of yielding new insights into topics ranging from the formation history of the Milky Way to the mass of the neutrino. However, enabling the astronomy community to fully exploit this multi-instrument data set is a challenging technical task: for much of the science, we will need to combine the photometry across multiple wavelengths with varying spectral and spatial resolution. We identify some of the key science enabled by the combined surveys and the key technical challenges in achieving the synergies.

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.

We present a cosmic shear study from the Deep Lens Survey (DLS), a deep BVRz multi-band imaging survey of five 4 sq. degree fields with two National Optical Astronomy Observatory (NOAO) 4-meter telescopes at Kitt Peak and Cerro Tololo. For both telescopes, the change of the point-spread-function (PSF) shape across the focal plane is complicated, and the exposure-to-exposure variation of this position-dependent 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 PSF-induced systematics from gravitational lensing effects. We calibrate our shears and estimate the errors, utilizing an image simulator, which generates sheared ground-based 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 cosmology-dependent 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 non-tomographic analysis alone constrains the Omega_M-sigma_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 weak-lensing tomographic study will further tighten these constraints because explicit treatment of the redshift dependence of cosmic shear more efficiently breaks the Omega_M-sigma_8 degeneracy. Combining the current results with the Wilkinson Microwave Anisotropy Probe 7-year (WMAP7) likelihood data, we obtain Omega_M=0.278+-0.018 and sigma_8=0.815+-0.020.

Using about 450,000 galaxies in the Deep Lens Survey, we present a detection of the gravitational magnification of z > 4 Lyman Break Galaxies by massive foreground galaxies with 0.4 < z < 1.0, grouped by redshift. The magnification signal is detected at S/N greater than 20, and rigorous checks confirm that it is not contaminated by any galaxy sample overlap in redshift. The inferred galaxy mass profiles are consistent with earlier lensing analyses at lower redshift. We then explore the tomographic lens magnification signal by splitting our foreground galaxy sample into 7 redshift bins. Combining galaxy-magnification cross-correlations and galaxy angular auto-correlations, we develop a bias-independent estimator of the tomographic signal. As a diagnostic of magnification tomography, the measurement of this estimator rejects a flat dark matter dominated Universe at > 7.5{\sigma} with a fixed \sigma_8 and is found to be consistent with the expected redshift-dependence of the WMAP7 {\Lambda}CDM cosmology.

We examine the projected correlation of galaxies with mass from small scales (<few hundred kpc) where individual dark matter halos dominate, out to 15 Mpc where correlated large-scale structure dominates. We investigate these profiles as a function of galaxy luminosity and redshift. Selecting 0.8 million galaxies in the Deep Lens Survey, we use photometric redshifts and stacked weak gravitational lensing shear tomography out to radial scales of 1 degree from the centers of foreground galaxies. We detect correlated mass density from multiple halos and large-scale structure at radii larger than the virial radius, and find the first observational evidence for growth in the galaxy-mass correlation on 10 Mpc scales with decreasing redshift and fixed range of luminosity. For a fixed range of redshift, we find a scaling of projected halo mass with rest-frame luminosity similar to previous studies at lower redshift. We control systematic errors in shape measurement and photometric redshift, enforce volume completeness through absolute magnitude cuts, and explore residual sample selection effects via simulations.

We present DLSCL J0916.2+2951 (z=0.53), a newly discovered major cluster merger in which the collisional cluster gas has become dissociated from the collisionless galaxies and dark matter. We identified the cluster using optical and weak lensing observations as part of the Deep Lens Survey. Our follow-up observations with Keck, Subaru, Hubble Space Telescope, and Chandra show that the cluster is a dissociative merger and constrain the dark matter self-interaction cross-section {\sigma}_{DM}m_{DM}^{-1}\leq7 cm^2g^{-1}. The system is observed at least 0.7\pm0.2 Gyr since first pass-through, thus providing a picture of cluster mergers 2-5 times further progressed than similar systems observed to date. This improved temporal leverage has implications for our understanding of merging clusters and their impact on galaxy evolution.

The uncertainty in the redshift distributions of galaxies has a significant potential impact on the cosmological parameter values inferred from multi-band imaging surveys. The accuracy of the photometric redshifts measured in these surveys depends not only on the quality of the flux data, but also on a number of modeling assumptions that enter into both the training set and SED fitting methods of photometric redshift estimation. In this work we focus on the latter, considering two types of modeling uncertainties: uncertainties in the SED template set and uncertainties in the magnitude and type priors used in a Bayesian photometric redshift estimation method. We find that SED template selection effects dominate over magnitude prior errors. We introduce a method for parameterizing the resulting ignorance of the redshift distributions, and for propagating these uncertainties to uncertainties in cosmological parameters.

The weak-lensing science of the LSST project drives the need to carefully model and separate the instrumental artifacts from the intrinsic lensing signal. The dominant source of the systematics for all ground based telescopes is the spatial correlation of the PSF modulated by both atmospheric turbulence and optical aberrations. In this paper, we present a full FOV simulation of the LSST images by modeling both the atmosphere and the telescope optics with the most current data for the telescope specifications and the environment. To simulate the effects of atmospheric turbulence, we generated six-layer phase screens with the parameters estimated from the on-site measurements. For the optics, we combined the ray-tracing tool ZEMAX and our simulated focal plane data to introduce realistic aberrations and focal plane height fluctuations. Although this expected flatness deviation for LSST is small compared with that of other existing cameras, the fast f-ratio of the LSST optics makes this focal plane flatness variation and the resulting PSF discontinuities across the CCD boundaries significant challenges in our removal of the systematics. We resolve this complication by performing PCA CCD-by-CCD, and interpolating the basis functions using conventional polynomials. We demonstrate that this PSF correction scheme reduces the residual PSF ellipticity correlation below 10^-7 over the cosmologically interesting scale. From a null test using HST/UDF galaxy images without input shear, we verify that the amplitude of the galaxy ellipticity correlation function, after the PSF correction, is consistent with the shot noise set by the finite number of objects. Therefore, we conclude that the current optical design and specification for the accuracy in the focal plane assembly are sufficient to enable the control of the PSF systematics required for weak-lensing science with the LSST.

Over the past decade, sky surveys such as the Sloan Digital Sky Survey have proven the power of large data sets for answering fundamental astrophysical questions. This observational progress, based on a synergy of advances in telescope construction, detectors, and information technology, has had a dramatic impact on nearly all fields of astronomy, and areas of fundamental physics. The next-generation instruments, and the surveys that will be made with them, will maintain this revolutionary progress. The hardware and computational technical challenges and the exciting science opportunities are attracting scientists and engineers from astronomy, optics, low-light-level detectors, high-energy physics, statistics, and computer science. The history of astronomy has taught us repeatedly that there are surprises whenever we view the sky in a new way. This will be particularly true of discoveries emerging from a new generation of sky surveys. Imaging data from large ground-based active optics telescopes with sufficient etendue can address many scientific missions simultaneously. These new investigations will rely on the statistical precision obtainable with billions of objects. For the first time, the full sky will be surveyed deep and fast, opening a new window on a universe of faint moving and distant exploding objects as well as unraveling the mystery of dark energy.

As our capacity to study ever-expanding domains of our science has increased (including the time domain, non-electromagnetic phenomena, magnetized plasmas, and numerous sky surveys in multiple wavebands with broad spatial coverage and unprecedented depths), so have the horizons of our understanding of the Universe been similarly expanding. This expansion is coupled to the exponential data deluge from multiple sky surveys, which have grown from gigabytes into terabytes during the past decade, and will grow from terabytes into Petabytes (even hundreds of Petabytes) in the next decade. With this increased vastness of information, there is a growing gap between our awareness of that information and our understanding of it. Training the next generation in the fine art of deriving intelligent understanding from data is needed for the success of sciences, communities, projects, agencies, businesses, and economies. This is true for both specialists (scientists) and non-specialists (everyone else: the public, educators and students, workforce). Specialists must learn and apply new data science research techniques in order to advance our understanding of the Universe. Non-specialists require information literacy skills as productive members of the 21st century workforce, integrating foundational skills for lifelong learning in a world increasingly dominated by data. We address the impact of the emerging discipline of data science on astronomy education within two contexts: formal education and lifelong learners.

We present weak lensing mass estimates of seven shear-selected galaxy cluster candidates from the Deep Lens Survey. The clusters were previously identified as mass peaks in convergence maps of 8.6 sq. deg of R band imaging, and followed up with X-ray and spectroscopic confirmation, spanning a redshift range 0.19 - 0.68. Most clusters contained multiple X-ray peaks, yielding 17 total mass concentrations. In this paper, we constrain the masses of these X-ray sources with weak lensing, using photometric redshifts from the full set of BVRz' imaging to properly weight background galaxies according to their lensing distance ratios. We fit both NFW and singular isothermal sphere profiles, and find that the results are insensitive to the assumed profile. We also show that the results do not depend significantly on the assumed prior on the position of the mass peak, but that this may become an issue in future larger samples. The inferred velocity dispersions for the extended X-ray sources range from 250-800 km/s, with the exception of one source for which no lensing signal was found. This work further establishes shear selection as a viable technique for finding clusters, but also highlights some unresolved issues such as determination of the mass profile center without biasing the mass estimate, and fully accounting for line-of-sight projections. A follow-up paper will examine the mass-X-ray scaling relations of these clusters.

We present a regularized maximum likelihood weak lensing reconstruction of the Deep Lens Survey F2 field (4 deg^2). High signal-to-noise ratio peaks in our lensing significance map appear to be associated with possible projected filamentary structures. The largest apparent structure extends for over a degree in the field and has contributions from known optical clusters at three redshifts (z ~ 0.3, 0.43, 0.5). Noise in weak lensing reconstructions is known to potentially cause "false positives"; we use Monte Carlo techniques to estimate the contamination in our sample, and find that 10-25% of the peaks are expected to be false detections. For significant lensing peaks we estimate the total signal-to-noise ratio of detection using a method that accounts for pixel-to-pixel correlations in our reconstruction. We also report the detection of a candidate relative underdensity in the F2 field with a total signal-to-noise ratio of ~ 5.5.

Wide-angle surveys have been an engine for new discoveries throughout the modern history of astronomy, and have been among the most highly cited and scientifically productive observing facilities in recent years. This trend is likely to continue over the next decade, as many of the most important questions in astrophysics are best tackled with massive surveys, often in synergy with each other and in tandem with the more traditional observatories. We argue that these surveys are most productive and have the greatest impact when the data from the surveys are made public in a timely manner. The rise of the "survey astronomer" is a substantial change in the demographics of our field; one of the most important challenges of the next decade is to find ways to recognize the intellectual contributions of those who work on the infrastructure of surveys (hardware, software, survey planning and operations, and databases/data distribution), and to make career paths to allow them to thrive.

The Cold Dark Matter theory of gravitationally-driven hierarchical structure formation has earned its status as a paradigm by explaining the distribution of matter over large spans of cosmic distance and time. However, its central tenet, that most of the matter in the universe is dark and exotic, is still unproven; the dark matter hypothesis is sufficiently audacious as to continue to warrant a diverse battery of tests. While local searches for dark matter particles or their annihilation signals could prove the existence of the substance itself, studies of cosmological dark matter in situ are vital to fully understand its role in structure formation and evolution. We argue that gravitational lensing provides the cleanest and farthest-reaching probe of dark matter in the universe, which can be combined with other observational techniques to answer the most challenging and exciting questions that will drive the subject in the next decade: What is the distribution of mass on sub-galactic scales? How do galaxy disks form and bulges grow in dark matter halos? How accurate are CDM predictions of halo structure? Can we distinguish between a need for a new substance (dark matter) and a need for new physics (departures from General Relativity)? What is the dark matter made of anyway? We propose that the central tool in this program should be a wide-field optical imaging survey, whose true value is realized with support in the form of high-resolution, cadenced optical/infra-red imaging, and massive-throughput optical spectroscopy.

We describe features of the LSST science database that are amenable to scientific data mining, object classification, outlier identification, anomaly detection, image quality assurance, and survey science validation. The data mining research agenda includes: scalability (at petabytes scales) of existing machine learning and data mining algorithms; development of grid-enabled parallel data mining algorithms; designing a robust system for brokering classifications from the LSST event pipeline (which may produce 10,000 or more event alerts per night); multi-resolution methods for exploration of petascale databases; indexing of multi-attribute multi-dimensional astronomical databases (beyond spatial indexing) for rapid querying of petabyte databases; and more.

We present a weak-lensing analysis of the galaxy cluster CL J1226+3332 at z=0.89 using Hubble Space Telescope Advanced Camera for Surveys images. The cluster is the hottest (>10 keV), most X-ray luminous system at z>0.6 known to date. The relaxed X-ray morphology, as well as its high temperature, is unusual at such a high redshift. Our mass reconstruction shows that on a large scale the dark matter distribution is consistent with a relaxed system with no significant substructures. However, on a small scale the cluster core is resolved into two mass clumps highly correlated with the cluster galaxy distribution. The dominant mass clump lies close to the brightest cluster galaxy whereas the other less massive clump is located ~40" (~310 kpc) to the southwest. Although this secondary mass clump does not show an excess in the X-ray surface brightness, the gas temperature of the region is much higher (12~18 keV) than those of the rest. We propose a scenario in which the less massive system has already passed through the main cluster and the X-ray gas has been stripped during this passage. The elongation of the X-ray peak toward the southwestern mass clump is also supportive of this possibility. We measure significant tangential shears out to the field boundary (~1.5 Mpc), which are well described by an Navarro-Frenk-White profile with a concentration parameter of c200=2.7+-0.3 and a scale length of rs=78"+-19" (~600 kpc) with chi^2/d.o.f=1.11. Within the spherical volume r200=1.6 Mpc, the total mass of the cluster becomes M(r<r200)=(1.4+-0.2) x 10^15 solar mass. Our weak-lensing analysis confirms that CL1226+3332 is indeed the most massive cluster known to date at z>0.6.

Baryon acoustic oscillations (BAOs) and weak lensing (WL) are complementary probes of cosmology. We explore the distance and growth factor measurements from photometric BAO and WL techniques and investigate the roles of the distance and growth factor in constraining dark energy. We find for WL that the growth factor has a great impact on dark energy constraints but is much less powerful than the distance. Dark energy constraints from WL are concentrated in considerably fewer distance eigenmodes than those from BAO, with the largest contributions from modes that are sensitive to the absolute distance. Both techniques have some well determined distance eigenmodes that are not very sensitive to the dark energy equation of state parameters w_0 and w_a, suggesting that they can accommodate additional parameters for dark energy and for the control of systematic uncertainties. A joint analysis of BAO and WL is far more powerful than either technique alone, and the resulting constraints on the distance and growth factor will be useful for distinguishing dark energy and modified gravity models. The Large Synoptic Survey Telescope (LSST) will yield both WL and angular BAO over a sample of several billion galaxies. Joint LSST BAO and WL can yield 0.5% level precision on ten comoving distances evenly spaced in log(1+z) between redshift 0.3 and 3 with cosmic microwave background priors from Planck. In addition, since the angular diameter distance, which directly affects the observables, is linked to the comoving distance solely by the curvature radius in the Friedmann-Robertson-Walker metric solution, LSST can achieve a pure metric constraint of 0.017 on the mean curvature parameter Omega_k of the universe simultaneously with the constraints on the comoving distances.

Next generation probes of dark matter and dark energy require high precision reconstruction of faint galaxy shapes from hundreds of dithered exposures. Current practice is to stack the images. While valuable for many applications, this stack is a highly compressed version of the data. Future weak lensing studies will require analysis of the full dataset using the stack and its associated catalog only as a starting point. We describe a "Multi-Fit" algorithm which simultaneously fits individual galaxy exposures to a common profile model convolved with each exposure's point spread function at that position in the image. This technique leads to an enhancement of the number of usable small galaxies at high redshift and, more significantly, a decrease in systematic shear error.

We compare X-ray and weak-lensing masses for four galaxy clusters that comprise the top-ranked shear-selected cluster system in the Deep Lens Survey. The weak-lensing observations of this system, which is associated with A781, are from the Kitt Peak Mayall 4-m telescope, and the X-ray observations are from both Chandra and XMM-Newton. For a faithful comparison of masses, we adopt the same matter density profile for each method, which we choose to be an NFW profile. Since neither the X-ray nor weak-lensing data are deep enough to well constrain both the NFW scale radius and central density, we estimate the scale radius using a fitting function for the concentration derived from cosmological hydrodynamic simulations and an X-ray estimate of the mass assuming isothermality. We keep this scale radius in common for both X-ray and weak-lensing profiles, and fit for the central density, which scales linearly with mass. We find that for three of these clusters, there is agreement between X-ray and weak-lensing NFW central densities, and thus masses. For the other cluster, the X-ray central density is higher than that from weak-lensing by 2 sigma. X-ray images suggest that this cluster may be undergoing a merger with a smaller cluster. This work serves as an additional step towards understanding the possible biases in X-ray and weak-lensing cluster mass estimation methods. Such understanding is vital to efforts to constrain cosmology using X-ray or weak-lensing cluster surveys to trace the growth of structure over cosmic time.

Since type Ia Supernovae (SNe) explode in galaxies, they can, in principle, be used as the same tracer of the large-scale structure as their hosts to measure baryon acoustic oscillations (BAOs). To realize this, one must obtain a dense integrated sampling of SNe over a large fraction of the sky, which may only be achievable photometrically with future projects such as the Large Synoptic Survey Telescope. The advantage of SN BAOs is that SNe have more uniform luminosities and more accurate photometric redshifts than galaxies, but the disadvantage is that they are transitory and hard to obtain in large number at high redshift. We find that a half-sky photometric SN survey to redshift z = 0.8 is able to measure the baryon signature in the SN spatial power spectrum. Although dark energy constraints from SN BAOs are weak, they can significantly improve the results from SN luminosity distances of the same data, and the combination of the two is no longer sensitive to cosmic microwave background priors.

(Abridged) The Large Synoptic Survey Telescope (LSST) is currently by far the most ambitious proposed ground-based optical survey. Solar System mapping is one of the four key scientific design drivers, with emphasis on efficient Near-Earth Object (NEO) and Potentially Hazardous Asteroid (PHA) detection, orbit determination, and characterization. In a continuous observing campaign of pairs of 15 second exposures of its 3,200 megapixel camera, LSST will cover the entire available sky every three nights in two photometric bands to a depth of V=25 per visit (two exposures), with exquisitely accurate astrometry and photometry. Over the proposed survey lifetime of 10 years, each sky location would be visited about 1000 times. The baseline design satisfies strong constraints on the cadence of observations mandated by PHAs such as closely spaced pairs of observations to link different detections and short exposures to avoid trailing losses. Equally important, due to frequent repeat visits LSST will effectively provide its own follow-up to derive orbits for detected moving objects. Detailed modeling of LSST operations, incorporating real historical weather and seeing data from LSST site at Cerro Pachon, shows that LSST using its baseline design cadence could find 90% of the PHAs with diameters larger than 250 m, and 75% of those greater than 140 m within ten years. However, by optimizing sky coverage, the ongoing simulations suggest that the LSST system, with its first light in 2013, can reach the Congressional mandate of cataloging 90% of PHAs larger than 140m by 2020.

Starting around 2013, data from the Large Synoptic Survey Telescope (LSST) will be analyzed for a wide range of phenomena. By separately tracing the development of mass structure and rate of expansion of the universe, these data will address the physics of dark matter and dark energy, the possible existence of modified gravity on large scales, large extra dimensions, the neutrino mass, and possible self interaction of dark matter particles.

We present the first sample of galaxy clusters selected on the basis of their weak gravitational lensing shear. The shear induced by a cluster is a function of its mass profile and its redshift relative to the background galaxies being sheared; in contrast to more traditional methods of selecting clusters, shear selection does not depend on the cluster's star formation history, baryon content, or dynamical state. Because mass is the property of clusters which provides constraints on cosmological parameters, the dependence on these other parameters could induce potentially important biases in traditionally-selected samples. Comparison of a shear-selected sample with optically and X-ray selected samples is therefore of great importance. Here we present the first step toward a new shear-selected sample: the selection of cluster candidates from the first 8.6 deg$^2$ of the 20 deg$^2$ Deep Lens Survey (DLS), and tabulation of their basic properties such as redshifts and optical and X-ray counterparts.

We consider cosmological applications of galaxy number density correlations to be inferred from future deep and wide multi-band optical surveys. We mostly focus on very large scales as a probe of possible features in the primordial power spectrum. We find the proposed survey of the Large Synoptic Survey Telescope may be competitive with future all-sky CMB experiments over a broad range of scales. On very large scales the inferred power spectrum is robust to photometric redshift errors, and, given a sufficient number density of galaxies, to angular variations in dust extinction and photometric calibration errors. We also consider other applications, such as constraining dark energy with the two CMB-calibrated standard rulers in the matter power spectrum, and controlling the effect of photometric redshift errors to facilitate the interpretation of cosmic shear data. We find that deep photometric surveys over wide area can provide constraints that are competitive with spectroscopic surveys in small volumes.

The Smithsonian Hectospec Lensing Survey (SHELS) combines a large deep complete redshift survey with a weak lensing map from the Deep Lens Survey (Wittman et al. 2002; 2005). We use maps of the velocity dispersion based on systems identified in the redshift survey to compare the three-dimensional matter distribution with the two-dimensional projection mapped by weak lensing. We demonstrate directly that the lensing map images the three-dimensional matter distribution obtained from the kinematic data.

Small distortions in the observed shapes of distant galaxies, a cosmic shear due to gravitational lensing, can be used to simultaneously determine the distance-redshift relation, r(z), and the density contrast growth factor, g(z). Both of these functions are sensitive probes of the acceleration. Their simultaneous determination allows for a consistency test and provides sensitivity to physics beyond the standard dark energy paradigm.

We present an analysis of weak lensing observations for RXJ1347-1145 over a 43' X 43' field taken in B and R filters on the Blanco 4m telescope at CTIO. RXJ1347-1145 is a massive cluster at redshift z=0.45. Using a population of galaxies with 20<R<26, we detect a weak lensing signal at the p<0.0005 level, finding best-fit parameters of \sigma_v=1400^{+130}_{-140} km s^{-1} for a singular isothermal sphere model and r_{200} = 3.5^{+0.8}_{-0.2} Mpc with c = 15^{+64}_{-10} for a NFW model in an \Omega_m = 0.3, \Omega_\Lambda = 0.7 cosmology. In addition, a mass to light ratio M/L_R =90 \pm 20 M_\odot / L_{R\odot} was determined. These values are consistent with the previous weak lensing study of RXJ1347--1145 by Fischer and Tyson, 1997, giving strong evidence that systemic bias was not introduced by the relatively small field of view in that study. Our best-fit parameter values are also consistent with recent X-ray studies by Allen et al, 2002 and Ettori et al, 2001, but are not consistent with recent optical velocity dispersion measurements by Cohen and Kneib, 2002.

We report on the methodology and first results from the Deep Lens Survey transient search. We utilize image subtraction on survey data to yield all sources of optical variability down to 24th magnitude. Images are analyzed immediately after acquisition, at the telescope and in near-real time, to allow for followup in the case of time-critical events. All classes of transients are posted to the web upon detection. Our observing strategy allows sensitivity to variability over several decades in timescale. The DLS is the first survey to classify and report all types of photometric and astrometric variability detected, including solar system objects, variable stars, supernovae, and short timescale phenomena. Three unusual optical transient events were detected, flaring on thousand-second timescales. All three events were seen in the B passband, suggesting blue color indices for the phenomena. One event (OT 20020115) is determined to be from a flaring Galactic dwarf star of spectral type dM4. From the remaining two events, we find an overall rate of \eta = 1.4 events deg-2 day-1 on thousand-second timescales, with a 95% confidence limit of \eta < 4.3. One of these events (OT 20010326) originated from a compact precursor in the field of galaxy cluster Abell 1836, and its nature is uncertain. For the second (OT 20030305) we find strong evidence for an extended extragalactic host. A dearth of such events in the R passband yields an upper 95% confidence limit on short timescale astronomical variability between 19.5 < R < 23.4 of \eta_R < 5.2. We report also on our ensemble of astrometrically variable objects, as well as an example of photometric variability with an undetected precursor.

We report on preliminary results of our X-ray survey of the most massive clusters currently identified from the Deep Lens Survey (DLS). The DLS cluster sample is selected based on weak lensing shear, which makes it possible for the first time to study clusters in a baryon-independent way. In this article we present X-ray properties of a subset of the shear-selected cluster sample.

We report the weak lensing discovery, spectroscopic confirmation, and weak lensing tomography of a massive cluster of galaxies at $z=0.68$, demonstrating that shear selection of clusters works at redshifts high enough to be cosmologically interesting. The mass estimate from weak lensing, $11.1 +- 2.8 x 10^{14} (r/Mpc)$ solar masses within projected radius r, agrees with that derived from the spectroscopy ($\sigma_v = 980 km s^{-1}$), and with the position of an arc which is likely to be a strongly lensed background galaxy. The redshift estimate from weak lensing tomography is consistent with the spectroscopy, demonstrating the feasibility of baryon-unbiased mass surveys. This tomographic technique will be able to roughly identify the redshifts of any dark clusters which may appear in shear-selected samples, up to z ~ 1.

A large wide-field telescope and camera with optical throughput over 200 m^2 deg^2 -- a factor of 50 beyond what we currently have -- would enable the detection of faint moving or bursting optical objects: from Earth threatening asteroids to energetic events at the edge of the optical universe. An optimized design for LSST is a 8.4 m telescope with a 3 degree field of view and an optical throughput of 260 m^2 deg^2. With its large throughput and dedicated all-sky monitoring mode, the LSST will reach 24th magnitude in a single 10 second exposure, opening unexplored regions of astronomical parameter space. The heart of the 2.3 Gpixel camera will be an array of imager modules with 10 micron pixels. Once each month LSST will survey up to 14,000 deg^2 of the sky with many ~10 second exposures. Over time LSST will survey 30,000 deg^2 deeply in multiple bandpasses, enabling innovative investigations ranging from galactic structure to cosmology. This is a shift in paradigm for optical astronomy: from "survey follow-up" to "survey direct science." The resulting real-time data products and fifteen petabyte time-tagged imaging database and photometric catalog will provide a unique resource. A collaboration of ~80 engineers and scientists is gearing up to confront this exciting challenge.

We measure seeing-corrected ellipticities for 2 x 10^6 galaxies with magnitude R<23 in 12 widely separated fields totalling 75 deg^2 of sky. At angular scales >30\arcmin, ellipticity correlations are detected at high significance and exhibit nearly the pure "E-mode" behavior expected of weak gravitational lensing. Even when smoothed to the full field size of 2.5 degrees, which is ~25h^-1 Mpc at the lens distances, an rms shear variance of <\gamma^2>^1/2 = 0.0012 +- 0.0003 is detected. At smaller angular scales there is significant "B-mode" power, an indication of residual uncorrected PSF distortions. The >30\arcmin data constrain the power spectrum of matter fluctuations on comoving scales of ~10h^-1 Mpc to have \sigma_8 (\Omega_m/0.3)^{0.57} = 0.71^{+0.12}_{-0.16} (95% CL, \LambdaCDM, \Gamma=0.21), where the systematic error includes statistical and calibration uncertainties, cosmic variance, and a conservative estimate of systematic contamination based upon the detected B-mode signal. This normalization of the power spectrum is lower than previous weak-lensing results but generally consistent them, is at the lower end of the \sigma_8 range from various analyses of galaxy cluster abundances, and agrees with recent determinations from CMB and galaxy clustering. The large and dispersed sky coverage of our survey reduces random errors and cosmic variance, while the relatively shallow depth allows us to use existing redshift-survey data to reduce systematic uncertainties in the N(z) distribution to insignificance. Reanalysis of the data with more sophisticated algorithms will hopefully reduce the systematic (B-mode) contamination, and allow more precise, multidimensional constraint of cosmological parameters.

The Deep Lens Survey (DLS) is a deep BVRz' imaging survey of seven 2x2 degree fields, with all data to be made public. The primary scientific driver is weak gravitational lensing, but the survey is also designed to enable a wide array of other astrophysical investigations. A unique feature of this survey is the search for transient phenomena. We subtract multiple exposures of a field, detect differences, classify, and release transients on the Web within about an hour of observation. Here we summarize the scientific goals of the DLS, field and filter selection, observing techniques and current status, data reduction, data products and release, and transient detections. Finally, we discuss some lessons which might apply to future large surveys such as LSST.

The number of mass clusters and their distribution in redshift are very sensitive to the density of matter Omega_m and the equation of state of dark energy w. Using weak lens gravitational tomography one can detect clusters of dark matter, weigh them, image their projected mass distribution, and determine their 3-D location. The degeneracy curve in the Omega_m - w plane is nearly orthogonal to that from CMB or SN measurements. Thus, a combination of CMB data with weak lens tomography of clusters can yield precision measurements of Omega_m and w, independently of the SN observations. The Large Synoptic Survey Telescope (LSST) will repeatedly survey 30,000 square degrees of the sky in multiple wavelengths. LSST will create a 3-D tomographic assay of mass overdensities back to half the age of the universe by measuring the shear and color-redshift of billions of high redshift galaxies. By simultaneously measuring several functions of cosmic shear and mass cluster abundance, LSST will provide a number of independent constraints on the dark energy density and the equation of state.

We report the discovery of a cluster of galaxies via its weak gravitational lensing effect on background galaxies, the first spectroscopically confirmed cluster to be discovered through its gravitational effects rather than by its electromagnetic radiation. This fundamentally different selection mechanism promises to yield mass-selected, rather than baryon or photon-selected, samples of these important cosmological probes. We have confirmed this cluster with spectroscopic redshifts of fifteen members at z=0.276, with a velocity dispersion of 615 km/s. We use the tangential shear as a function of source photometric redshift to estimate the lens redshift independently and find z_l = 0.30 +- 0.08. The good agreement with the spectroscopy indicates that the redshift evolution of the mass function may be measurable from the imaging data alone in shear-selected surveys.

Weak gravitational lensing enables direct reconstruction of dark matter maps over cosmologically significant volumes. This research is currently telescope-limited. The Dark Matter Telescope (DMT) is a proposed 8.4 m telescope with a 3 degree field of view, with an etendue of 260 $(m. degree)^2$, ten times greater than any other current or planned telescope. With its large etendue and dedicated observational mode, the DMT fills a nearly unexplored region of parameter space and enables projects that would take decades on current facilities. The DMT will be able to reach 10-sigma limiting magnitudes of 27-28 magnitude in the wavelength range .3 - 1 um over a 7 square degree field in 3 nights of dark time. Here we review its unique weak lensing cosmology capabilities and the design that enables those capabilities.

Most of the matter in the universe is not luminous and can be observed directly only through its gravitational effect. An emerging technique called weak gravitational lensing uses background galaxies to reveal the foreground dark matter distribution on large scales. Light from very distant galaxies travels to us through many intervening overdensities which gravitationally distort their apparent shapes. The observed ellipticity pattern of these distant galaxies thus encodes information about the large-scale structure of the universe, but attempts to measure this effect have been inconclusive due to systematic errors. We report the first detection of this ``cosmic shear'' using 145,000 background galaxies to reveal the dark matter distribution on angular scales up to half a degree in three separate lines of sight. The observed angular dependence of this effect is consistent with that predicted by two leading cosmological models, providing new and independent support for these models.

One hundred years ago we did not know how stars generate energy, the age of the Universe was thought to be only millions of years, and our Milky Way galaxy was the only galaxy known. Today, we know that we live in an evolving and expanding Universe comprising billions of galaxies, all held together by dark matter. With the hot big-bang model, we can trace the evolution of the Universe from the hot soup of quarks and leptons that existed a fraction of a second after the beginning to the formation of galaxies a few billion years later, and finally to the Universe we see today 13 billion years after the big bang, with its clusters of galaxies, superclusters, voids, and great walls. The attractive force of gravity acting on tiny primeval inhomogeneities in the distribution of matter gave rise to all the structure seen today. A paradigm based upon deep connections between cosmology and elementary particle physics -- inflation + cold dark matter -- holds the promise of extending our understanding to an even more fundamental level and much earlier times, as well as shedding light on the unification of the forces and particles of nature. As we enter the 21st century, a flood of observations is testing this paradigm.

We construct a high resolution mass map of the z=0.39 cluster 0024+1654, based on parametric inversion of the associated gravitational lens. The lens creates eight well-resolved sub-images of a background galaxy, seen in deep imaging with HST. Excluding mass concentrations centered on visible galaxies, more than 98% of the remaining mass is represented by a smooth concentration of dark matter centered near the brightest cluster galaxies, with a 35 h^{-1} kpc soft core. The asymmetry in the mass distribution is <3% inside 107 ~h^{-1} kpc radius. The dark matter distribution we observe in CL0024 is far more smooth, symmetric, and nonsingular than in typical simulated clusters in either Omega=1 or Omega=0.3 CDM cosmologies. Integrated to 107 h^{-1} kpc radius, the rest-frame mass to light ratio is M/L_V = 276\pm 40 h (M/L_V)_solar.

New HST WFPC2 observations of the lensed double QSO 0957+561 will allow improved constraints on the lens mass distribution and hence will improve the derived value of H$_0$. We first present improved optical positions and photometry for the known components of this lens. The optical separation between the A and B quasar images agrees with VLBI data at the 10 mas level, and the optical center of the primary lensing galaxy G1 coincides with the VLBI source G' to within 10 mas. The best previous model for this lens (Grogin and Narayan 1996) is excluded by these data and must be reevaluated. Several new resolved features are found within 10\arcsec of G1, including an apparent fold arc with two bright knots. Several other small galaxies are detected, including two which may be multiple images of each other. We present positions and crude photometry of these objects.

Galaxy cluster mass distribution are potentially useful probes of $\Omega_0$ and the nature of the dark matter. Large clusters will distort the observed shapes of background galaxies through gravitational lensing allowing the measurement of the cluster mass distributions. In this paper we describe weak statistical lensing measurements of the most luminous X-ray cluster known, RXJ1347.5-1145 at z=0.45. We detect a shear signal in the background galaxies at a signal-to-noise ratio of 7.5 in the radial range $120 \le r \le 1360 h^{-1}$ kpc. A mass map of the cluster reveals an 11$\sigma$ peak in the cluster mass distribution consistent with the position of the central dominant galaxy and 3 $\sigma$ evidence for the presence of a subcluster at a projected radius of 1.3 - 1.7 h^{-1} Mpc from the cluster center. In the range $120 \le r \le 1360 h^{-1}$ kpc mass traces light, and the azimuthally averaged cluster mass and light profiles are consistent with singular isothermal spheres with M(r<1 Mpc) = $1.7 \pm 0.4 \times 10^{15} M_\odot$. Assuming an isotropic velocity distribution function, the implied velocity dispersion is $\sigma = 1500 \pm 160 km s^{-1}$. The rest-frame mass-to-light ratio is $M/L_B = 200 \pm 50 h M_\odot/L_{B\odot}$. The lensing mass estimate is almost twice as high as a previously determined X-ray mass estimate.

Multiply gravitationally lensed objects with known time delays can lead to direct determinations of H$_0$ independent of the distance ladder if the mass distribution of the lens is known. Currently, the double QSO 0957+561 is the only lensed object with a precisely known time delay. The largest remaining source of systematic error in the H$_0$ determination results from uncertainty in the mass distribution of the lens which is comprised of a massive galaxy (G1) and the cluster in which it resides. We have obtained V-band CCD images from CFHT in order to measure the mass distribution in the cluster from its gravitional distorting effect on the appearance of background galaxes. We use this data to constuct a two-dimensional mass map of the field. A mass peak is detected at the $4.5\sigma$ level, offset from, but consistent with, the position of G1. Simple tests reveal no significant substructure and the mass distribution is consistent with a spherical cluster. The peak in the number density map of bright galaxies is offset from G1 similarly to the mass peak. We constructed an azimuthally averaged mass profile centered on G1 out to 2 \arcmin ($400 h^{-1}$ kpc). It is consistent with an isothermal mass distribution with a small core ($r_c \approx 5 \arcsec = 17 h^{-1}$ kpc). The inferred mass within 1 Mpc is consistent with the dynamical mass estimate but $2\sigma$ higher than the upper limits from a ROSAT X-ray study. We discuss implications for H$_0$ in a future paper.

We present the results of a search for low surface brightness galaxies (hereafter LSBs) in the Coma cluster. Bernstein et al. report on deep CCD observations in R of a 7.5 by 7.5 arcminute region in the core of the Coma cluster, and we extend this work by finding and measuring 36 LSBs within this field. We report both R and Bj results. The average magnitude based on the best fit exponential to the images is 22.5 (R) and the typical exponential scale is 1.3 arcseconds. The range of exponential scales is 0.4 to 1.2 kpc (distance modulus 34.89), and the range of central surface brightnesses is 24 to 27.4 R magnitudes per square arcsecond. Many of these objects are similar in terms of scale length and central surface brightness to those found by others in nearby clusters such as Fornax (Bothun, Impey and Malin), as well as in the low luminosity end of the dwarfs cataloged in the review by Ferguson and Binggeli. We find no evidence for a dependence of color on central surface brightness or on distance from the D galaxies or the X-ray center of Coma. We also find that these LSBs make a small contribution to the overall mass of the cluster. We discuss these results in terms of possible scenarios of LSB formation and evolution.

In this Letter we present calibrated observations of the average mass of galaxies with 22< I < 25 in the Hubble Deep Field. We measure the mean mass profile using the statistical gravitational lens distortion of faint blue background galaxies, based on 2221 foreground-background pairs. The observed weak lensing distortion is calibrated via full 3-D simulations and other HST data on a foreground cluster of known mass. Inside a projected radius of 5 arcsec we find a 3\sigma detection of shear. Fitting to an isothermal mass distribution, we find a shear of 0.066^{+0.018}_{-0.019} at a radius of 2\arcsec\ (\sim 8 ~ h^{-1} kpc). The (1\sigma) limits on the parameters for the truncated isothermal model are \sigma_v =185^{+30}_{-35} km s^{-1}, r_{outer} \ge 15 ~ h^{-1} kpc, and r_{core}= 0.6^{+1.3}_{-0.6} ~ h^{-1} kpc. This corresponds to an average galaxy mass interior to 20 h^{-1} kpc of 5.9^{+2.5}_{-2.7}\times 10^{11} h^{-1} M_\odot. Inside 10 h^{-1} kpc, we find an average rest-frame V band mass-to-light ratio of 11.4^{+6.0}_{-6.7} h (M/L_V)_\odot.

Optical variability in galaxies at high redshift is a tracer of evolution in AGN activity, and should provide a useful constraint on models of galaxy evolution, AGN structure, and cosmology. We studied optical variability in multiple deep CCD and photographic surveys of blank fields for galaxies with $B_j = 20 - 25$ mag. Weakly variable objects are far more common than strongly variable ones. For objects near $B_j = 22$, $0.74\% \pm 0.2 \%$ vary by 0.026~mag RMS or more, over a decade. This is small compared with previous claims based on photographic surveys, and also small compared with the fraction of bright quasars ($\approx 5\%$ at $B_j = 20$~mag) or Seyferts ($\approx 1-2\%$ for $B_j < 18$). The fraction of objects that vary increases slowly with magnitude. Detection probabilities and error rates were checked by simulations and statistical analysis of fluctuations of sample sky spots.

A unique reconstruction of the image of a high redshift source galaxy responsible for multiple long arcs in the z = 0.4 cluster 0024+1654 is obtained by inverse lensing. Deep B and I imaging with the Hubble Space Telescope (Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by AURA, under NASA contract NAS 5-26555) enables high resolution of the arcs due to strong gravitational lensing of the background source. The gravitational lens in the foreground cluster is thus used to obtain a magnified view of the distant source. Four strongly lensed images of the source lead to a unique reconstruction. Each of the long arcs, when unlensed, leads to the same reconstructed source image exhibiting a beaded ring-like morphology. The U luminosity of the ring alone is equivalent to a normal galaxy. This is likely a galaxy in formation.

We present a digital catalog of images of 113 galaxies in this paper. These galaxies are all nearby, bright, large and well resolved. All images were recorded with charge coupled devices (CCDs) at the Palomar Observatory with the 1.5 meter telescope and at the Lowell Observatory with the 1.1 meter telescope. At Palomar we used the Thuan--Gunn g, r and i photometric bands to take 3 images each of 31 spiral galaxies; at Lowell we used the B_J and R bands (2 images per galaxy) of the photometric system by Gullixson et al. (1995) to observe 82 spirals and ellipticals. The galaxies were selected to span the Hubble classification classes. All data are photometrically calibrated with foreground stars removed. Important data on these galaxies published in the "Third Reference Catalog of Bright Galaxies" (RC3) are recorded in the FITS file headers. All files are available through anonymous FTP from ftp://astro.princeton.edu/, through WWW at http://astro.princeton.edu/~frei/galaxy_catalog.html, and Princeton University Press will soon publish the data on CD-ROM.

In this letter we present calibrated mass and light profiles of the rich cluster of galaxies Abell 1689 out to 1 $h^{-1}$ Mpc from the center. The high surface density of faint blue galaxies at high redshift, selected by their low surface brightness, are unique tools for mapping the projected mass distribution of foreground mass concentrations. The systematic gravitational lens distortions of $10^4$ of these background galaxies in 15\arcmin\ fields reveal detailed mass profiles for intervening clusters of galaxies, and are a direct measure of the growth of mass inhomogeneity. The mass is measured directly, avoiding uncertainties encountered in velocity or X-ray derived mass estimates. Mass in the rich cluster Abell 1689 follows smoothed light, outside 100 h$^{-1}$ kpc, with a rest-frame V band mass-to-light ratio of $400 \pm 60$ $h^{-1} (M/L_V)_\odot$. Near the cluster center, mass appears to be more smoothly distributed than light. Out to a radius of 1 $h^{-1}$ Mpc the total mass follows a steeper than isothermal profile. Comparing with preliminary high resolution N-body clustering simulations for various cosmogonies on these scales, these data are incompatible with hot dark matter, a poor fit to most mixed dark matter models, and favor open or $\Lambda > 0$ cold dark matter. Substructure is seen in both the mass and the light, but detailed correspondence is erased on scales less than 100 $h^{-1}$ kpc.

We determine the luminosity function (LF) of galaxies in the core of the Coma cluster for M_R<=-11.4 (assuming H_0=75 km/s/Mpc), a magnitude regime previously explored only in the Local Group. Objects are counted in a deep CCD image of Coma having RMS noise of 27.7 R mag~arcsec$^{-2}$. A correction for objects in the foreground or background of the Coma cluster---and the uncertainty in this correction---are determined from images of five other high-latitude fields, carefully matched to the Coma image in both resolution and noise level. Accurate counts of Coma cluster members are obtained as faint as R=25.5, or M_R=-9.4. The LF for galaxies is well fit by a power law dN/dL\propto L^\alpha, with \alpha=-1.42\pm0.05, over the range -19.4<M_R<-11.4; faintward of this range, the galaxies are unresolved and indistinguishable from globular clusters, but the data are consistent with an extrapolation of the power law. Surface brightness biases are minimized since galaxies are not subjected to morphological selection, and the limiting detection isophote is 27.6 R mag~arcsec$^{-2}$. We find the typical $M_R\approx-12$ Coma cluster galaxy to have an exponential scale length $\approx200$~pc, similar to Local Group galaxies of comparable magnitude. These extreme dwarf galaxies show a surface density increasing towards the giant elliptical NGC~4874 as $r^{-1.3}$, similar to the diffuse light and globular cluster distributions. The luminosity in the detected dwarf galaxies is at most a few percent of the total diffuse light of the giant galaxies in the cluster, and the contribution of the dwarfs to the mass of the cluster is likely negligible as well.

Luminous objects associated with the 7.06\arcsec\ separation double quasar Q2345+007 have, until now, escaped detection. In this letter we present the results of the deepest known imaging of the region surrounding the quasar. The total exposure times were 47600 seconds in B$_j$ (101 frames) and 32160 seconds in R (82 frames). The frames came from CFHT, CTIO, and KPNO. We detect a B$_j$ = 25.0 mag galaxy (B$_j$ -- R = 0.5 mag) in close proximity to the fainter QSO image, possibly at $z = 1.49$ given by several absorption features in the QSO spectra. Furthermore, there is a 2.9$\sigma$ enhancement in the number density of faint galaxies ($24 \le$ B$_j$ $\le 28$, B$_j-$ R $ \le 1.5$) near the quasar and another 3.1$\sigma$ enhancement further away. These discoveries support the hypothesis that Q2345+007 is being lensed by one or more distant mass concentrations and may imply that compact ``clusters'' already exist at early epochs. We discuss several lens models for the system and the cosmological implications.