Through the lens of the LSST Science Collaborations' experience, this paper advocates for new and improved ways to fund large, complex collaborations at the interface of data science and astrophysics as they work in preparation for and on peta-scale, complex surveys, of which LSST is a prime example. We advocate for the establishment of programs to support both research and infrastructure development that enables innovative collaborative research on such scales.

Galaxy color gradients - i.e., spectral energy distributions that vary across the galaxy profile - will impact galaxy shape measurements when the modeled point spread function (PSF) corresponds to that for a galaxy with spatially uniform color. This paper describes the techniques and results of a study of the expected impact of galaxy color gradients on weak lensing measurements with the Large Synoptic Survey Telescope (LSST) when the PSF size depends on wavelength. The bias on cosmic shear measurements from color gradients is computed both for parametric bulge+disk galaxy simulations and for more realistic chromatic galaxy surface brightness profiles based on HST V- and I-band images in the AEGIS survey. For the parametric galaxies, and for the more realistic galaxies derived from AEGIS galaxies with sufficient SNR that color gradient bias can be isolated, the predicted multiplicative shear biases due to color gradients are found to be at least a factor of 2 below the LSST full-depth requirement on the total systematic uncertainty in the redshift-dependent shear calibration. The analysis code and data products are publicly available (https://github.com/sowmyakth/measure_cg_bias).

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.

The Large Synoptic Survey Telescope (LSST) Dark Energy Science Collaboration (DESC) will use five cosmological probes: galaxy clusters, large scale structure, supernovae, strong lensing, and weak lensing. This Science Requirements Document (SRD) quantifies the expected dark energy constraining power of these probes individually and together, with conservative assumptions about analysis methodology and follow-up observational resources based on our current understanding and the expected evolution within the field in the coming years. We then define requirements on analysis pipelines that will enable us to achieve our goal of carrying out a dark energy analysis consistent with the Dark Energy Task Force definition of a Stage IV dark energy experiment. This is achieved through a forecasting process that incorporates the flowdown to detailed requirements on multiple sources of systematic uncertainty. Future versions of this document will include evolution in our software capabilities and analysis plans along with updates to the LSST survey strategy.

One of the scientific goals of the Large Synoptic Survey Telescope (LSST) is to measure the evolution of dark energy by measuring subtle distortions of galaxy shapes due to weak gravitational lensing caused by the evolving dark matter distribution. Understanding the point spread function (PSF) for LSST is a crucial step to accurate measurements of weak gravitational lensing. Atmospheric contributions dominate the LSST PSF. Simulations of Kolmogorov turbulence models are commonly used to characterize and correct for these atmospheric effects. In order to validate these simulations, we compare the predicted atmospheric behavior to empirical data. The simulations are carried out in GalSim, an open-source software package for simulating images of astronomical objects and PSFs. Atmospheric simulations are run by generating large phase screens at varying altitude and evolving them over long time scales. We compare the turbulence strength and temporal behavior of atmospheres generated from simulations to those from reconstructed telemetry data from the Gemini Planet Imager (GPI). GPI captures a range of spatial frequencies by sampling the atmosphere with 18-cm subapertures. The LSST weak lensing analysis will measure correlations of galaxy ellipticity, requiring very accurate knowledge of the magnitude and correlations of PSF shape parameters. Following from the first analysis, we use simulations and sequential short exposure observations from the Differential Speckle Survey Instrument (DSSI) to study the behavior of PSF parameters - e.g., ellipticity and size - as a function of exposure time. These studies could help inform discussions of possible variable exposure times for LSST visits for example, to provide more uniform depth of visits.

(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.

Wavelength-dependent point spread functions (PSFs) violate an implicit assumption in current galaxy shape measurement algorithms that deconvolve the PSF measured from stars (which have stellar spectral energy distributions (SEDs)) from images of galaxies (which have galactic SEDs). Since the absorption length of silicon depends on wavelength, CCDs are a potential source of PSF chromaticity. Here we develop two toy models to estimate the sensitivity of the cosmic shear survey from the Large Synoptic Survey Telescope to chromatic effects in CCDs. We then compare these toy models to simulated estimates of PSF chromaticity derived from the LSST photon simulator PhoSim. We find that even though sensor contributions to PSF chromaticity are subdominant to atmospheric contributions, they can still significantly bias cosmic shear results if left uncorrected, particularly in the redder filter bands and for objects that are off-axis in the field of view.

Current and future imaging surveys will measure cosmic shear with statistical precision that demands a deeper understanding of potential systematic biases in galaxy shape measurements than has been achieved to date. We use analytic and computational techniques to study the impact on shape measurements of two atmospheric chromatic effects for ground-based surveys such as the Dark Energy Survey and the Large Synoptic Survey Telescope (LSST): (i) atmospheric differential chromatic refraction and (ii) wavelength dependence of seeing. We investigate the effects of using the point spread function (PSF) measured with stars to determine the shapes of galaxies that have different spectral energy distributions than the stars. We find that both chromatic effects lead to significant biases in galaxy shape measurements for current and future surveys, if not corrected. Using simulated galaxy images, we find a form of chromatic `model bias' that arises when fitting a galaxy image with a model that has been convolved with a stellar, instead of galactic, point spread function. We show that both forms of atmospheric chromatic biases can be predicted (and corrected) with minimal model bias by applying an ordered set of perturbative PSF-level corrections based on machine-learning techniques applied to six-band photometry. Catalog-level corrections do not address the model bias. We conclude that achieving the ultimate precision for weak lensing from current and future ground-based imaging surveys requires a detailed understanding of the wavelength dependence of the PSF from the atmosphere, and from other sources such as optics and sensors. The source code for this analysis is available at https://github.com/DarkEnergyScienceCollaboration/chroma .

We employ robust weak gravitational lensing measurements to improve cosmological constraints from measurements of the galaxy cluster mass function and its evolution, using X-ray selected clusters detected in the ROSAT All-Sky Survey. Our lensing analysis constrains the absolute mass scale of such clusters at the 8 per cent level, including both statistical and systematic uncertainties. Combining it with the survey data and X-ray follow-up observations, we find a tight constraint on a combination of the mean matter density and late-time normalization of the matter power spectrum, $\sigma_8(\Omega_m/0.3)^{0.17}=0.81\pm0.03$, with marginalized, one-dimensional constraints of $\Omega_m=0.26\pm0.03$ and $\sigma_8=0.83\pm0.04$. For these two parameters, this represents a factor of two improvement in precision with respect to previous work, primarily due to the reduced systematic uncertainty in the absolute mass calibration provided by the lensing analysis. Our new results are in good agreement with constraints from cosmic microwave background (CMB) data, both WMAP and Planck (plus WMAP polarization), under the assumption of a flat $\Lambda$CDM cosmology with minimal neutrino mass. Consequently, we find no evidence for non-minimal neutrino mass from the combination of cluster data with CMB, supernova and baryon acoustic oscillation measurements, regardless of which all-sky CMB data set is used (and independent of the recent claimed detection of B-modes on degree scales). We also present improved constraints on models of dark energy (both constant and evolving), modifications of gravity, and primordial non-Gaussianity. Assuming flatness, the constraints for a constant dark energy equation of state from the cluster data alone are at the 15 per cent level, improving to $\sim 6$ per cent when the cluster data are combined with other leading probes.

In light of the tension in cosmological constraints reported by the Planck team between their SZ-selected cluster counts and Cosmic Microwave Background (CMB) temperature anisotropies, we compare the Planck cluster mass estimates with robust, weak-lensing mass measurements from the Weighing the Giants (WtG) project. For the 22 clusters in common between the Planck cosmology sample and WtG, we find an overall mass ratio of $\left< M_{Planck}/M_{\rm WtG} \right> = 0.688 \pm 0.072$. Extending the sample to clusters not used in the Planck cosmology analysis yields a consistent value of $\left< M_{Planck}/M_{\rm WtG} \right> = 0.698 \pm 0.062$ from 38 clusters in common. Identifying the weak-lensing masses as proxies for the true cluster mass (on average), these ratios are $\sim 1.6\sigma$ lower than the default mass bias of 0.8 assumed in the Planck cluster analysis. Adopting the WtG weak-lensing-based mass calibration would substantially reduce the tension found between the Planck cluster count cosmology results and those from CMB temperature anisotropies, thereby dispensing of the need for "new physics" such as uncomfortably large neutrino masses (in the context of the measured Planck temperature anisotropies and other data). We also find modest evidence (at 95 per cent confidence) for a mass dependence of the calibration ratio and discuss its potential origin in light of systematic uncertainties in the temperature calibration of the X-ray measurements used to calibrate the Planck cluster masses. Our results exemplify the critical role that robust absolute mass calibration plays in cluster cosmology, and the invaluable role of accurate weak-lensing mass measurements in this regard.

We present improved methods for using stars found in astronomical exposures to calibrate both star and galaxy colors as well as to adjust the instrument flat field. By developing a spectroscopic model for the SDSS stellar locus in color-color space, synthesizing an expected stellar locus, and simultaneously solving for all unknown zeropoints when fitting to the instrumental locus, we increase the calibration accuracy of stellar locus matching. We also use a new combined technique to estimate improved flat-field models for the Subaru SuprimeCam camera, forming `star flats' based on the magnitudes of stars observed in multiple positions or through comparison with available SDSS magnitudes. These techniques yield galaxy magnitudes with reliable color calibration (< 0.01 - 0.02 mag accuracy) that enable us to estimate photometric redshift probability distributions without spectroscopic training samples. We test the accuracy of our photometric redshifts using spectroscopic redshifts z_s for ~5000 galaxies in 27 cluster fields with at least five bands of photometry, as well as galaxies in the COSMOS field, finding sigma((z_p - z_s)/(1 + z_s)) ~ 0.03 for the most probable redshift z_p. We show that the full posterior probability distributions for the redshifts of galaxies with five-band photometry exhibit good agreement with redshifts estimated from thirty-band photometry in the COSMOS field. The growth of shear with increasing distance behind each galaxy cluster shows the expected redshift-distance relation for a flat Lambda-CDM cosmology. Photometric redshifts and calibrated colors are used in subsequent papers to measure the masses of 51 galaxy clusters from their weak gravitational shear. We make our Python code for stellar locus matching available at http://big-macs-calibrate.googlecode.com; the code requires only a catalog and filter functions.

We report weak-lensing masses for 51 of the most X-ray luminous galaxy clusters known. This cluster sample, introduced earlier in this series of papers, spans redshifts 0.15 < z_cl < 0.7, and is well suited to calibrate mass proxies for current cluster cosmology experiments. Cluster masses are measured with a standard `color-cut' lensing method from three-filter photometry of each field. Additionally, for 27 cluster fields with at least five-filter photometry, we measure high-accuracy masses using a new method that exploits all information available in the photometric redshift posterior probability distributions of individual galaxies. Using simulations based on the COSMOS-30 catalog, we demonstrate control of systematic biases in the mean mass of the sample with this method, from photometric redshift biases and associated uncertainties, to better than 3%. In contrast, we show that the use of single-point estimators in place of the full photometric redshift posterior distributions can lead to significant redshift-dependent biases on cluster masses. The performance of our new photometric redshift-based method allows us to calibrate `color-cut` masses for all 51 clusters in the present sample to a total systematic uncertainty of ~7% on the mean mass, a level sufficient to significantly improve current cosmology constraints from galaxy clusters. Our results bode well for future cosmological studies of clusters, potentially reducing the need for exhaustive spectroscopic calibration surveys as compared to other techniques, when deep, multi-filter optical and near-IR imaging surveys are coupled with robust photometric redshift methods.

This is the first in a series of papers in which we measure accurate weak-lensing masses for 51 of the most X-ray luminous galaxy clusters known at redshifts 0.15<z<0.7, in order to calibrate X-ray and other mass proxies for cosmological cluster experiments. The primary aim is to improve the absolute mass calibration of cluster observables, currently the dominant systematic uncertainty for cluster count experiments. Key elements of this work are the rigorous quantification of systematic uncertainties, high-quality data reduction and photometric calibration, and the "blind" nature of the analysis to avoid confirmation bias. Our target clusters are drawn from RASS X-ray catalogs, and provide a versatile calibration sample for many aspects of cluster cosmology. We have acquired wide-field, high-quality imaging using the Subaru and CFHT telescopes for all 51 clusters, in at least three bands per cluster. For a subset of 27 clusters, we have data in at least five bands, allowing accurate photo-z estimates of lensed galaxies. In this paper, we describe the cluster sample and observations, and detail the processing of the SuprimeCam data to yield high-quality images suitable for robust weak-lensing shape measurements and precision photometry. For each cluster, we present wide-field color optical images and maps of the weak-lensing mass distribution, the optical light distribution, and the X-ray emission, providing insights into the large-scale structure in which the clusters are embedded. We measure the offsets between X-ray centroids and Brightest Cluster Galaxies in the clusters, finding these to be small in general, with a median of 20kpc. For offsets <100kpc, weak-lensing mass measurements centered on the BCGs agree well with values determined relative to the X-ray centroids; miscentering is therefore not a significant source of systematic uncertainty for our mass measurements. [abridged]

Current and future imaging surveys will measure cosmic shear with a statistical precision that demands a deeper understanding of potential systematic biases in galaxy shape measurements than has been achieved to date. We investigate the effects of using the point spread function (PSF) measured with stars to determine the shape of a galaxy that has a different spectral energy distribution (SED) than the star. We demonstrate that a wavelength dependent PSF size, for example as may originate from atmospheric seeing or the diffraction limit of the primary aperture, can introduce significant shape measurement biases. This analysis shows that even small wavelength dependencies in the PSF may introduce biases, and hence that achieving the ultimate precision for weak lensing from current and future imaging surveys will require a detailed understanding of the wavelength dependence of the PSF from all sources, including the CCD sensors.

We discuss some of the key science questions that are bringing particle physicists and astrophysicists together, and comment on some of the cultural and funding issues that have arisen as these two communities become increasingly intertwined.