The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

Weak lensing, which is the deflection of light by matter along the line of sight, has proven to be an efficient method for constraining models of structure formation and reveal the nature of dark energy. So far, most weak-lensing studies have focused on the shear field that can be measured directly from the ellipticity of background galaxies. However, within the context of forthcoming full-sky weak-lensing surveys such as Euclid, convergence maps (mass maps) offer an important advantage over shear fields in terms of cosmological exploitation. While it carry the same information, the lensing signal is more compressed in the convergence maps than in the shear field. This simplifies otherwise computationally expensive analyses, for instance, non-Gaussianity studies. However, the inversion of the non-local shear field requires accurate control of systematic effects caused by holes in the data field, field borders, shape noise, and the fact that the shear is not a direct observable (reduced shear). We present the two mass-inversion methods that are included in the official Euclid data-processing pipeline: the standard Kaiser & Squires method (KS), and a new mass-inversion method (KS+) that aims to reduce the information loss during the mass inversion. This new method is based on the KS method and includes corrections for mass-mapping systematic effects. The results of the KS+ method are compared to the original implementation of the KS method in its simplest form, using the Euclid Flagship mock galaxy catalogue. In particular, we estimate the quality of the reconstruction by comparing the two-point correlation functions and third- and fourth- order moments obtained from shear and convergence maps, and we analyse each systematic effect independently and simultaneously...

Forthcoming space-based observations will require high-quality point-spread function (PSF) models for weak gravitational lensing measurements. One approach to generating these models is using a wavefront model based on the known telescope optics. We present an empirical framework for validating such models to confirm that they match the actual PSF to within requirements by comparing the models to the observed light distributions of isolated stars. We apply this framework to Tiny Tim, the standard tool for generating model PSFs for the Hubble Space Telescope (HST), testing its models against images taken by HST's Advanced Camera for Surveys in the Wide Field Channel. We show that Tiny Tim's models, in the default configuration, differ significantly from the observed PSFs, most notably in their sizes. We find that the quality of Tiny Tim PSFs can be improved through fitting the full set of Zernike polynomial coefficients which characterise the optics, to the point where the practical significance of the difference between model and observed PSFs is negligible for most use cases, resulting in additive and multiplicative biases both of order approximately 4e-4. We also show that most of this improvement can be retained through using an updated set of Zernike coefficients, which we provide.

The data from the Euclid mission will enable the measurement of the photometric redshifts, angular positions, and weak lensing shapes for over a billion galaxies. This large dataset will allow for cosmological analyses using the angular clustering of galaxies and cosmic shear. The cross-correlation (XC) between these probes can tighten constraints and it is therefore important to quantify their impact for Euclid. In this study we carefully quantify the impact of XC not only on the final parameter constraints for different cosmological models, but also on the nuisance parameters. In particular, we aim at understanding the amount of additional information that XC can provide for parameters encoding systematic effects, such as galaxy bias or intrinsic alignments (IA). We follow the formalism presented in Euclid Collaboration: Blanchard et al. (2019) and make use of the codes validated therein. We show that XC improves the dark energy Figure of Merit (FoM) by a factor $\sim 5$, whilst it also reduces the uncertainties on galaxy bias by $\sim 17\%$ and the uncertainties on IA by a factor $\sim 4$. We observe that the role of XC on the final parameter constraints is qualitatively the same irrespective of the galaxy bias model used. We also show that XC can help in distinguishing between different IA models, and that if IA terms are neglected then this can lead to significant biases on the cosmological parameters. We find that the XC terms are necessary to extract the full information content from the data in future analyses. They help in better constraining the cosmological model, and lead to a better understanding of the systematic effects that contaminate these probes. Furthermore, we find that XC helps in constraining the mean of the photometric-redshift distributions, but it requires a more precise knowledge of this mean in order not to degrade the final FoM. [Abridged]

Context. Future weak lensing surveys, such as the Euclid mission, will attempt to measure the shapes of billions of galaxies in order to derive cosmological information. These surveys will attain very low levels of statistical error, and systematic errors must be extremely well controlled. In particular, the point spread function (PSF) must be estimated using stars in the field, and recovered with high accuracy. Aims. The aims of this paper are twofold. Firstly, we took steps toward a nonparametric method to address the issue of recovering the PSF field, namely that of finding the correct PSF at the position of any galaxy in the field, applicable to Euclid. Our approach relies solely on the data, as opposed to parametric methods that make use of our knowledge of the instrument. Secondly, we studied the impact of imperfect PSF models on the shape measurement of galaxies themselves, and whether common assumptions about this impact hold true in an Euclid scenario. Methods. We extended the recently proposed resolved components analysis approach, which performs super-resolution on a field of under-sampled observations of a spatially varying, image-valued function. We added a spatial interpolation component to the method, making it a true 2-dimensional PSF model. We compared our approach to PSFEx, then quantified the impact of PSF recovery errors on galaxy shape measurements through image simulations. Results. Our approach yields an improvement over PSFEx in terms of the PSF model and on observed galaxy shape errors, though it is at present far from reaching the required Euclid accuracy. We also find that the usual formalism used for the propagation of PSF model errors to weak lensing quantities no longer holds in the case of an Euclid-like PSF. In particular, different shape measurement approaches can react differently to the same PSF modeling errors.

We present a publicly-available toolkit of flight-proven hardware and software to retrieve 5 TB of data or small physical samples from a stratospheric balloon platform. Before launch, a capsule is attached to the balloon, and rises with it. Upon remote command, the capsule is released and descends via parachute, continuously transmitting its location. Software to predict the trajectory can be used to select a safe but accessible landing site. We dropped two such capsules from the SuperBIT telescope, in September 2019. The capsules took ~37 minutes to descend from ~30 km altitude. They drifted 32 km and 19 km horizontally, but landed within 300 m and 600 m of their predicted landing sites. We found them easily, and successfully recovered the data. We welcome interest from other balloon teams for whom the technology would be useful.

Stage IV weak lensing experiments will offer more than an order of magnitude leap in precision. We must therefore ensure that our analyses remain accurate in this new era. Accordingly, previously ignored systematic effects must be addressed. In this work, we evaluate the impact of the reduced shear approximation and magnification bias, on the information obtained from the angular power spectrum. To first-order, the statistics of reduced shear, a combination of shear and convergence, are taken to be equal to those of shear. However, this approximation can induce a bias in the cosmological parameters that can no longer be neglected. A separate bias arises from the statistics of shear being altered by the preferential selection of galaxies and the dilution of their surface densities, in high-magnification regions. The corrections for these systematic effects take similar forms, allowing them to be treated together. We calculated the impact of neglecting these effects on the cosmological parameters that would be determined from Euclid, using cosmic shear tomography. To do so, we employed the Fisher matrix formalism, and included the impact of the super-sample covariance. We also demonstrate how the reduced shear correction can be calculated using a lognormal field forward modelling approach. These effects cause significant biases in Omega_m, sigma_8, n_s, Omega_DE, w_0, and w_a of -0.53 sigma, 0.43 sigma, -0.34 sigma, 1.36 sigma, -0.68 sigma, and 1.21 sigma, respectively. We then show that these lensing biases interact with another systematic: the intrinsic alignment of galaxies. Accordingly, we develop the formalism for an intrinsic alignment-enhanced lensing bias correction. Applying this to Euclid, we find that the additional terms introduced by this correction are sub-dominant.

The Euclid mission will observe well over a billion galaxies out to $z\sim6$ and beyond. This will offer an unrivalled opportunity to investigate several key questions for understanding galaxy formation and evolution. The first step for many of these studies will be the selection of a sample of quiescent and star-forming galaxies, as is often done in the literature by using well known colour techniques such as the `UVJ' diagram. However, given the limited number of filters available for the Euclid telescope, the recovery of such rest-frame colours will be challenging. We therefore investigate the use of observed Euclid colours, on their own and together with ground-based u-band observations, for selecting quiescent and star-forming galaxies. The most efficient colour combination, among the ones tested in this work, consists of the (u-VIS) and (VIS-J) colours. We find that this combination allows users to select a sample of quiescent galaxies complete to above $\sim70\%$ and with less than 15$\%$ contamination at redshifts in the range $0.75<z<1$. For galaxies at high-z or without the u-band complementary observations, the (VIS-Y) and (J-H) colours represent a valid alternative, with $>65\%$ completeness level and contamination below 20$\%$ at $1<z<2$ for finding quiescent galaxies. In comparison, the sample of quiescent galaxies selected with the traditional UVJ technique is only $\sim20\%$ complete at $z<3$, when recovering the rest-frame colours using mock Euclid observations. This shows that our new methodology is the most suitable one when only Euclid bands, along with u-band imaging, are available.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

At a fraction the total cost of an equivalent orbital mission, scientific balloon-borne platforms, operating above 99.7% of the Earth's atmosphere, offer attractive, competitive, and effective observational capabilities -- namely space-like resolution, transmission, and backgrounds -- that are well suited for modern astronomy and cosmology. SuperBIT is a diffraction-limited, wide-field, 0.5 m telescope capable of exploiting these observing conditions in order to provide exquisite imaging throughout the near-IR to near-UV. It utilizes a robust active stabilization system that has consistently demonstrated a 1 sigma sky-fixed pointing stability at 48 milliarcseconds over multiple 1 hour observations at float. This is achieved by actively tracking compound pendulations via a three-axis gimballed platform, which provides sky-fixed telescope stability at < 500 milliarcseconds and corrects for field rotation, while employing high-bandwidth tip/tilt optics to remove residual disturbances across the science imaging focal plane. SuperBIT's performance during the 2019 commissioning flight benefited from a customized high-fidelity science-capable telescope designed with exceptional thermo- and opto-mechanical stability as well as tightly constrained static and dynamic coupling between high-rate sensors and telescope optics. At the currently demonstrated level of flight performance, SuperBIT capabilities now surpass the science requirements for a wide variety of experiments in cosmology, astrophysics and stellar dynamics.

We study the shape of eight dynamically relaxed galaxy clusters observed with the Hubble Space Telescope and Chandra X-Ray Observatory. Using strong and weak gravitational lensing, the shape of the Brightest Cluster Galaxy and the X-ray isophote we study the ellipticity of the cluster halo at four different radii. We find that the proxies probing the inner regions of the cluster are strongly coupled with the BCG shape correlated with both the shape predicted by strong gravitational lensing and the X-ray isophote. Conversely we find no such correlation between the shape as predicted by the weak lensing and the other three probes suggesting any coupling between the inner and outer regions is weak. We also present in this paper the public release of the HST weak lensing shape measurement code pyRRG, directly available from PyPi (https://pypi.org/project/pyRRG/). This python3.7 code, based on Rhodes et al. (2000) adopts an automated star-galaxy classifier based on a Random Forest and outputs scientifically useful products such as weak lensing catalogues suitable for the mass mapping algorithm Lenstool.

We provide predictions of the yield of $7<z<9$ quasars from the Euclid wide survey, updating the calculation presented in the Euclid Red Book in several ways. We account for revisions to the Euclid near-infrared filter wavelengths; we adopt steeper rates of decline of the quasar luminosity function (QLF; $\Phi$) with redshift, $\Phi\propto10^{k(z-6)}$, $k=-0.72$, and a further steeper rate of decline, $k=-0.92$; we use better models of the contaminating populations (MLT dwarfs and compact early-type galaxies); and we use an improved Bayesian selection method, compared to the colour cuts used for the Red Book calculation, allowing the identification of fainter quasars, down to $J_{AB}\sim23$. Quasars at $z>8$ may be selected from Euclid $OYJH$ photometry alone, but selection over the redshift interval $7<z<8$ is greatly improved by the addition of $z$-band data from, e.g., Pan-STARRS and LSST. We calculate predicted quasar yields for the assumed values of the rate of decline of the QLF beyond $z=6$. For the case that the decline of the QLF accelerates beyond $z=6$, with $k=-0.92$, Euclid should nevertheless find over 100 quasars with $7.0<z<7.5$, and $\sim25$ quasars beyond the current record of $z=7.5$, including $\sim8$ beyond $z=8.0$. The first Euclid quasars at $z>7.5$ should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, $7<z<8$, $M_{1450}<-25$, using 8m class telescopes to confirm candidates, but follow-up with JWST or E-ELT will be required to measure the faint-end slope. Contamination of the candidate lists is predicted to be modest even at $J_{AB}\sim23$. The precision with which $k$ can be determined over $7<z<8$ depends on the value of $k$, but assuming $k=-0.72$ it can be measured to a 1 sigma uncertainty of 0.07.

The Wide-Field InfraRed Space Telescope (WFIRST) will be capable of delivering precise astrometry for faint sources over the enormous field of view of its main camera, the Wide-Field Imager (WFI). This unprecedented combination will be transformative for the many scientific questions that require precise positions, distances, and velocities of stars. We describe the expectations for the astrometric precision of the WFIRST WFI in different scenarios, illustrate how a broad range of science cases will see significant advances with such data, and identify aspects of WFIRST's design where small adjustments could greatly improve its power as an astrometric instrument.

Our aim is to quantify the impact of systematic effects on the inference of cosmological parameters from cosmic shear. We present an end-to-end approach that introduces sources of bias in a modelled weak lensing survey on a galaxy-by-galaxy level. Residual biases are propagated through a pipeline from galaxy properties (one end) through to cosmic shear power spectra and cosmological parameter estimates (the other end), to quantify how imperfect knowledge of the pipeline changes the maximum likelihood values of dark energy parameters. We quantify the impact of an imperfect correction for charge transfer inefficiency (CTI) and modelling uncertainties of the point spread function (PSF) for Euclid, and find that the biases introduced can be corrected to acceptable levels.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

Joint survey processing (JSP) is the pixel level combination of LSST, Euclid, and WFIRST datasets. By combining the high spatial resolution of the space-based datasets with deep, seeing-limited, ground-based images in the optical bands, systematics like source confusion and astrometric mismatch can be addressed to derive the highest precision optical/infrared photometric catalogs. This white paper highlights the scientific motivation, computational and algorithmic needs to build joint pixel level processing capabilities, which the individual projects by themselves will not be able to support. Through this white paper, we request that the Astro2020 decadal committee recognize the JSP effort as a multi-agency project with the natural outcome being a collaborative effort among groups which are normally supported by a single agency. JSP will allow the U.S. (and international) astronomical community to manipulate the flagship data sets and undertake innovative science investigations ranging from solar system object characterization, exoplanet detections, nearby galaxy rotation rates and dark matter properties, to epoch of reionization studies. It will also result in the ultimate constraints on cosmological parameters and the nature of dark energy, with far smaller uncertainties and a better handle on systematics than by any one survey alone.

We speculate on the development and availability of new innovative propulsion techniques in the 2040s, that will allow us to fly a spacecraft outside the Solar System (at 150 AU and more) in a reasonable amount of time, in order to directly probe our (gravitational) Solar System neighborhood and answer pressing questions regarding the dark sector (dark energy and dark matter). We identify two closely related main science goals, as well as secondary objectives that could be fulfilled by a mission dedicated to probing the local dark sector: (i) begin the exploration of gravitation's low-acceleration regime with a man-made spacecraft and (ii) improve our knowledge of the local dark matter and baryon densities. Those questions can be answered by directly measuring the gravitational potential with an atomic clock on-board a spacecraft on an outbound Solar System orbit, and by comparing the spacecraft's trajectory with that predicted by General Relativity through the combination of ranging data and the in-situ measurement (and correction) of non-gravitational accelerations with an on-board accelerometer. Despite a wealth of new experiments getting online in the near future, that will bring new knowledge about the dark sector, it is very unlikely that those science questions will be closed in the next two decades. More importantly, it is likely that it will be even more urgent than currently to answer them. Tracking a spacecraft carrying a clock and an accelerometer as it leaves the Solar System may well be the easiest and fastest way to directly probe our dark environment.

The Wide Field Infrared Survey Telescope (WFIRST) Coronagraph Instrument (CGI) will be the first high-performance stellar coronagraph using active wavefront control for deep starlight suppression in space, providing unprecedented levels of contrast, spatial resolution, and sensitivity for astronomical observations in the optical. One science case enabled by the CGI will be taking images and(R~50)spectra of faint interplanetary dust structures present in the habitable zone of nearby sunlike stars (~10 pc) and within the snow-line of more distant ones(~20pc), down to dust density levels commensurate with that of the solar system zodiacal cloud. Reaching contrast levels below~10-7 for the first time, CGI will cross an important threshold in debris disks physics, accessing disks with low enough optical depths that their structure is dominated by transport phenomena than collisions. Hence, CGI results will be crucial for determining how exozodiacal dust grains are produced and transported in low-density disks around mature stars. Additionally, CGI will be able to measure the brightness level and constrain the degree of asymmetry of exozodiacal clouds around individual nearby sunlike stars in the optical, at the ~10x solar zodiacal emission level. This information will be extremely valuable for optimizing the observational strategy of possible future exo-Earth direct imaging missions, especially those planning to operate at optical wavelengths, such as Habitable Exoplanet Observatory (HabEx) and the Large Ultraviolet/Optical/Infrared Surveyor (LUVOIR).

ATLAS (Astrophysics Telescope for Large Area Spectroscopy) is a concept for a NASA probe-class space mission. It is the spectroscopic follow-up mission to WFIRST, boosting its scientific return by obtaining deep NIR & MIR slit spectroscopy for most of the galaxies imaged by the WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for ~200M galaxies out to z=7 and beyond, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS and WFIRST together will produce a definitive 3D map of the Universe over 2000 sq deg. ATLAS Science Goals are: (1) Discover how galaxies have evolved in the cosmic web of dark matter from cosmic dawn through the peak era of galaxy assembly. (2) Discover the nature of cosmic acceleration. (3) Probe the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy. (4) Discover the bulk compositional building blocks of planetesimals formed in the outer Solar System. These flow down to the ATLAS Scientific Objectives: (1A) Trace the relation between galaxies and dark matter with less than 10% shot noise on relevant scales at 1<z<7. (1B) Probe the physics of galaxy evolution at 1<z<7. (2) Obtain definitive measurements of dark energy and tests of General Relativity. (3) Measure the 3D structure and stellar content of the inner Milky Way to a distance of 25 kpc. (4) Detect and quantify the composition of 3,000 planetesimals in the outer Solar System. ATLAS is a 1.5m telescope with a FoV of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 1000, and a wavelength range of 1-4 microns. ATLAS has an unprecedented spectroscopic capability based on DMDs, with a spectroscopic multiplex factor ~6,000. ATLAS is designed to fit within the NASA probe-class space mission cost envelope.

High-multiplex and deep spectroscopic follow-up of upcoming panoramic deep-imaging surveys like LSST, Euclid, and WFIRST is a widely recognized and increasingly urgent necessity. No current or planned facility at a U.S. observatory meets the sensitivity, multiplex, and rapid-response time needed to exploit these future datasets. FOBOS, the Fiber-Optic Broadband Optical Spectrograph, is a near-term fiber-based facility that addresses these spectroscopic needs by optimizing depth over area and exploiting the aperture advantage of the existing 10m Keck II Telescope. The result is an instrument with a uniquely blue-sensitive wavelength range (0.31-1.0 um) at R~3500, high-multiplex (1800 fibers), and a factor 1.7 greater survey speed and order-of-magnitude greater sampling density than Subaru's Prime Focus Spectrograph (PFS). In the era of panoramic deep imaging, FOBOS will excel at building the deep, spectroscopic reference data sets needed to interpret vast imaging data. At the same time, its flexible focal plane, including a mode with 25 deployable integral-field units (IFUs) across a 20 arcmin diameter field, enables an expansive range of scientific investigations. Its key programmatic areas include (1) nested stellar-parameter training sets that enable studies of the Milky Way and M31 halo sub-structure, as well as local group dwarf galaxies, (2) a comprehensive picture of galaxy formation thanks to detailed mapping of the baryonic environment at z~2 and statistical linking of evolving populations to the present day, and (3) dramatic enhancements in cosmological constraints via precise photometric redshifts and determined redshift distributions. In combination with Keck I instrumentation, FOBOS also provides instant access to medium-resolution spectroscopy for transient sources with full coverage from the UV to the K-band.

In order for Wide-Field Infrared Survey Telescope (WFIRST) and other Stage IV dark energy experiments (e.g., Large Synoptic Survey Telescope; LSST, and Euclid) to infer cosmological parameters not limited by systematic errors, accurate redshift measurements are needed. This accuracy can be met by using spectroscopic subsamples to calibrate the photometric redshifts for the full sample. In this work we employ the Self Organizing Map (SOM) spectroscopic sampling technique, to find the minimal number of spectra required for the WFIRST weak lensing calibration. We use galaxies from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) to build the LSST+WFIRST lensing analog sample of ~36 k objects and train the LSST+WFIRST SOM. We find that 26% of the WFIRST lensing sample consists of sources fainter than the Euclid depth in the optical, 91% of which live in color cells already occupied by brighter galaxies. We demonstrate the similarity between faint and bright galaxies as well as the feasibility of redshift measurements at different brightness levels. Our results suggest that the spectroscopic sample acquired for calibration to the Euclid depth is sufficient for calibrating the majority of the WFIRST color-space. For the spectroscopic sample to fully represent the synthetic color-space of WFIRST, we recommend obtaining additional spectroscopy of ~0.2-1.2 k new sources in cells occupied by mostly faint galaxies. We argue that either the small area of the CANDELS fields and the small overall sample size or the large photometric errors might be the reason for no/less bright galaxies mapped to these cells. Acquiring the spectra of these sources will confirm the above findings and will enable the comprehensive calibration of the WFIRST color-redshift relation.

This white paper is the result of the Tri-Agency Working Group (TAG) appointed to develop synergies between missions and is intended to clarify what LSST observations are needed in order to maximally enhance the combined science output of LSST and Euclid. To facilitate LSST planning we provide a range of possible LSST surveys with clear metrics based on the improvement in the Dark Energy figure of merit (FOM). To provide a quantifiable metric we present five survey options using only between 0.3 and 3.8% of the LSST 10 year survey. We also provide information so that the LSST DDF cadence can possibly be matched to those of \emph{Euclid} in common deep fields, SXDS, COSMOS, CDFS, and a proposed new LSST deep field (near the Akari Deep Field South). Co-coordination of observations from the Large Synoptic Survey Telescope (LSST) and Euclid will lead to a significant number of synergies. The combination of optical multi-band imaging from LSST with high resolution optical and near-infrared photometry and spectroscopy from \emph{Euclid} will not only improve constraints on Dark Energy, but provide a wealth of science on the Milky Way, local group, local large scale structure, and even on first galaxies during the epoch of reionization. A detailed paper has been published on the Dark Energy science case (Rhodes et al.) by a joint LSST/Euclid working group as well as a white paper describing LSST/Euclid/WFIRST synergies (Jain et al.), and we will briefly describe other science cases here. A companion white paper argues the general science case for an extension of the LSST footprint to the north at airmass < 1.8, and we support the white papers for southern extensions of the LSST survey.

We propose an extension of the LSST survey to cover the northern sky to DEC < +30 (accessible at airmass <1.8). This survey will increase the LSST sky coverage by ~9,600 square degrees from 18,900 to 28,500 square degrees (a 50% increase) but use only 0.6-2.5% of the time depending on the synergies with other surveys. This increased area addresses a wide range of science cases that enhance all of the primary LSST science goals by significant amounts. The science enabled includes: increasing the area of the sky accessible for follow-up of multi-messenger transients including gravitational waves, mapping the milky way halo and halo dwarfs including discovery of RR Lyrae stars in the outer galactic halo, discovery of z>7 quasars in combination Euclid, enabling a second generation DESI and other spectroscopic surveys, and enhancing all areas of science by improving synergies with Euclid, WFIRST, and unique northern survey facilities. This white paper is the result of the Tri-Agency Working Group (TAG) appointed to develop synergies between missions and presents a unified plan for northern coverage. The range of time estimates reflects synergies with other surveys. If the modified DESC WFD survey, the ecliptic plane mini survey, and the north galactic spur mini survey are executed this plan would only need 0.6% of the LSST time, however if none of these are included the overall request is 2.5% of the 10 year survey life. In other words, the majority of these observations are already suggested as part of these other surveys and the intent of this white paper is to propose a unified baseline plan to carry out a broad range of objectives to facilitate a combination of multiple science objectives. A companion white paper gives Euclid specific science goals, and we support the white papers for southern extensions of the LSST survey.

The Complete Calibration of the Color-Redshift Relation (C3R2) survey is a multi-institution, multi-instrument survey that aims to map the empirical relation of galaxy color to redshift to i~24.5 (AB), thereby providing a firm foundation for weak lensing cosmology with the Stage IV dark energy missions Euclid and WFIRST. Here we present 3171 new spectroscopic redshifts obtained in the 2016B and 2017A semesters with a combination of DEIMOS, LRIS, and MOSFIRE on the Keck telescopes. The observations come from all of the Keck partners: Caltech, NASA, the University of Hawaii, and the University of California. Combined with the 1283 redshifts published in DR1, the C3R2 survey has now obtained and published 4454 high quality galaxy redshifts. We discuss updates to the survey design and provide a catalog of photometric and spectroscopic data. Initial tests of the calibration method performance are given, indicating that the sample, once completed and combined with extensive data collected by other spectroscopic surveys, should allow us to meet the cosmology requirements for Euclid, and make significant headway toward solving the problem for WFIRST. We use the full spectroscopic sample to demonstrate that galaxy brightness is weakly correlated with redshift once a galaxy is localized in the Euclid or WFIRST color space, with potentially important implications for the spectroscopy needed to calibrate redshifts for faint WFIRST and LSST sources.

Two decades after its discovery, cosmic acceleration remains the most profound mystery in cosmology and arguably in all of physics. Either the Universe is dominated by a form of dark energy with exotic physical properties not predicted by standard model physics, or General Relativity is not an adequate description of gravity over cosmic distances. WFIRST emerged as a top priority of Astro2010 in part because of its ability to address the mystery of cosmic acceleration through both high precision measurements of the cosmic expansion history and the growth of cosmic structures with multiple and redundant probes. We illustrate in this white paper how mission design changes since Astro2010 have made WFIRST an even more powerful dark energy facility and have improved the ability of WFIRST to respond to changes in the experimental landscape. WFIRST is the space-based probe of DE the community needs in the mid-2020s.

The Cosmic Dawn and Reionization epochs remain a fundamental but challenging frontier of astrophysics and cosmology. We advocate a large-scale, multi-tracer approach to develop a comprehensive understanding of the physics that led to the formation and evolution of the first stars and galaxies. We highlight the line intensity mapping technique to trace the multi-phase reionization topology on large scales, and measure reionization history in detail. Besides 21cm, we advocate for Lya tomography mapping during the epoch of Wouthuysen-Field coupling as an additional probe of the cosmic dawn era.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.)

Free-floating planets are the remnants of violent dynamical rearrangements of planetary systems. It is possible that even our own solar system ejected a large planet early in its evolution. WFIRST will have the ability to detect free-floating planets over a wide range of masses, but it will not be able to directly measure their masses. Microlensing parallax observations can be used to measure the masses of isolated objects, including free-floating planets, by observing their microlensing events from two locations. The intra-L2 separation between WFIRST and Euclid is large enough to enable microlensing parallax measurements, especially given the exquisite photometric precision that both spacecraft are capable of over wide fields. In this white paper we describe how a modest investment of observing time could yield hundreds of parallax measurements for WFIRST's bound and free-floating planets. We also describe how a short observing campaign of precursor observations by Euclid can improve WFIRST's bound planet and host star mass measurements.

The Cosmic Dawn Intensity Mapper (CDIM) will transform our understanding of the era of reionization when the Universe formed the first stars and galaxies, and UV photons ionized the neutral medium. CDIM goes beyond the capabilities of upcoming facilities by carrying out wide area spectro-imaging surveys, providing redshifts of galaxies and quasars during reionization as well as spectral lines that carry crucial information on their physical properties. CDIM will make use of unprecedented sensitivity to surface brightness to measure the intensity fluctuations of reionization on large-scales to provide a valuable and complementary dataset to 21-cm experiments. The baseline mission concept is an 83-cm infrared telescope equipped with a focal plane of 24 \times 20482 detectors capable of R = 300 spectro-imaging observations over the wavelength range of 0.75 to 7.5 {\mu}m using Linear Variable Filters (LVFs). CDIM provides a large field of view of 7.8 deg2 allowing efficient wide area surveys, and instead of moving instrumental components, spectroscopic mapping is obtained through a shift-and-stare strategy through spacecraft operations. CDIM design and capabilities focus on the needs of detecting faint galaxies and quasars during reionization and intensity fluctuation measurements of key spectral lines, including Lyman-{\alpha} and H{\alpha} radiation from the first stars and galaxies. The design is low risk, carries significant science and engineering margins, and makes use of technologies with high technical readiness level for space observations.

Galaxies form and evolve in the context of their local and large-scale environments. Their baryonic content that we observe with imaging and spectroscopy is intimately connected to the properties of their dark matter halos, and to their location in the "cosmic web" of large-scale structure. Very large spectroscopic surveys of the local universe (e.g., SDSS and GAMA) measure galaxy positions (location within large-scale structure), statistical clustering (a direct constraint on dark matter halo masses), and spectral features (measuring physical conditions of the gas and stars within the galaxies, as well as internal velocities). Deep surveys with the James Webb Space Telescope (JWST) will revolutionize spectroscopic measurements of redshifts and spectral properties for galaxies out to the epoch of reionization, but with numerical statistics and over cosmic volumes that are too small to map large-scale structure and to constrain halo properties via clustering. Here, we consider advances in understanding galaxy evolution that would be enabled by very large spectroscopic surveys at high redshifts: very large numbers of galaxies (outstanding statistics) over large co-moving volumes (large-scale structure on all scales) over broad redshift ranges (evolution over most of cosmic history). The required observational facility can be established as part of the probe portfolio by NASA within the next decade.

The nature of dark energy remains a profound mystery 20 years after the discovery of cosmic acceleration. A very high number density galaxy redshift survey over a wide area (HD GRS Wide) spanning the redshift range of 0.5<z<4 using the same tracer, carried out using massively parallel wide field multi-object slit spectroscopy from space, will provide definitive dark energy measurements with minimal observational systematics by design. The HD GRS Wide will illuminate the nature of dark energy, and lead to revolutionary advances in particle physics and cosmology. It will also trace the cosmic web of dark matter and provide key insight into large-scale structure in the Universe. The required observational facility can be established as part of the probe portfolio by NASA within the next decade.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

The outer regions of globular clusters can enable us to answer many fundamental questions concerning issues ranging from the formation and evolution of clusters and their multiple stellar populations to the study of stars near and beyond the hydrogen-burning limit and to the dynamics of the Milky Way. The outskirts of globular clusters are still uncharted territories observationally. A very efficient way to explore them is through high-precision proper motions of low-mass stars over a large field of view. The Wide Field InfraRed Survey Telescope (WFIRST) combines all these characteristics in a single telescope, making it the best observational tool to uncover the wealth of information contained in the clusters' outermost regions.

Astrophysical transients have been observed for millennia and have shaped our most basic assumptions about the Universe. In the last century, systematic searches have grown from detecting handfuls of transients per year to over 7000 in 2018 alone. As these searches have matured, we have discovered both large samples of "normal" classes and new, rare classes. Recently, a transient was the first object observed in both gravitational waves and light. Ground-based observatories, including LSST, will discover thousands of transients in the optical, but these facilities will not provide the high-fidelity near-infrared (NIR) photometry and high-resolution imaging of a space-based observatory. WFIRST can fill this gap. With its survey designed to measure the expansion history of the Universe with Type Ia supernovae, WFIRST will also discover and monitor thousands of other transients in the NIR, revealing the physics for these high-energy events. Small-scale GO programs, either as a supplement to the planned survey or as specific target-of-opportunity observations, would significantly expand the scope of transient science that can be studied with WFIRST.

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

This white paper highlights the crucial and urgent synergies required between WFIRST, Subaru Hyper Suprime-Cam or other >25m-class telescopes galaxy observations and SKA 21cm measurements to constrain the nature of reionization (ionization history and topology) and its sources.

The Wide Field Infrared Survey Telescope (WFIRST) is a 2.4m space telescope with a 0.281 deg^2 field of view for near-IR imaging and slitless spectroscopy and a coronagraph designed for > 10^8 starlight suppresion. As background information for Astro2020 white papers, this article summarizes the current design and anticipated performance of WFIRST. While WFIRST does not have the UV imaging/spectroscopic capabilities of the Hubble Space Telescope, for wide field near-IR surveys WFIRST is hundreds of times more efficient. Some of the most ambitious multi-cycle HST Treasury programs could be executed as routine General Observer (GO) programs on WFIRST. The large area and time-domain surveys planned for the cosmology and exoplanet microlensing programs will produce extraordinarily rich data sets that enable an enormous range of Archival Research (AR) investigations. Requirements for the coronagraph are defined based on its status as a technology demonstration, but its expected performance will enable unprecedented observations of nearby giant exoplanets and circumstellar disks. WFIRST is currently in the Preliminary Design and Technology Completion phase (Phase B), on schedule for launch in 2025, with several of its critical components already in production.

ATLAS (Astrophysics Telescope for Large Area Spectroscopy) Probe is a concept for a NASA probe-class space mission. It is the follow-up space mission to WFIRST, boosting its scientific return by obtaining deep IR slit spectroscopy for 70% of all galaxies imaged by a 2000 sq deg WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for 200M galaxies out to z < 7, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS Probe science spans four broad categories: (1) Revolutionizing galaxy evolution studies by tracing the relation between galaxies and dark matter from galaxy groups to cosmic voids and filaments, from the epoch of reionization through the peak era of galaxy assembly; (2) Opening a new window into the dark Universe by weighing the dark matter filaments using 3D weak lensing with spectroscopic redshifts, and obtaining definitive measurements of dark energy and modification of General Relativity using galaxy clustering; (3) Probing the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy; and (4) Exploring the formation history of the outer Solar System by characterizing Kuiper Belt Objects. ATLAS Probe is a 1.5m telescope with a field of view of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 1000 over 1-4 microns, and a spectroscopic multiplex factor >5,000. ATLAS is designed to fit within the NASA probe-class space mission cost envelope; it has a single instrument, a telescope aperture that allows for a lighter launch vehicle, and mature technology. ATLAS Probe will lead to transformative science over the entire range of astrophysics: from galaxy evolution to the dark Universe, from Solar System objects to the dusty regions of the Milky Way.

The Coronagraph Instrument (CGI) for NASA's Wide Field Infrared Survey Telescope (WFIRST) will constitute a dramatic step forward for high-contrast imaging, integral field spectroscopy, and polarimetry of exoplanets and circumstellar disks, aiming to improve upon the sensitivity of current ground-based direct imaging facilities by 2-3 orders of magnitude. Furthermore, CGI will serve as a pathfinder for future exo-Earth imaging and characterization missions by demonstrating wavefront control, coronagraphy, and spectral retrieval in a new contrast regime, and by validating instrument and telescope models at unprecedented levels of precision. To achieve this jump in performance, it is critical to draw on the experience of ground-based high-contrast facilities. We discuss several areas of relevant commonalities, including: wavefront control, post-processing of integral field unit data, and calibration and observing strategies.

The Wide Field Infrared Survey Telescope (WFIRST) Coronagraph Instrument (CGI) is a high-contrast imager and integral field spectrograph that will enable the study of exoplanets and circumstellar disks at visible wavelengths. Ground-based high-contrast instrumentation has fundamentally limited performance at small working angles, even under optimistic assumptions for 30m-class telescopes. There is a strong scientific driver for better performance, particularly at visible wavelengths. Future flagship mission concepts aim to image Earth analogues with visible light flux ratios of more than 10^10. CGI is a critical intermediate step toward that goal, with a predicted 10^8-9 flux ratio capability in the visible. CGI achieves this through improvements over current ground and space systems in several areas: (i) Hardware: space-qualified (TRL9) deformable mirrors, detectors, and coronagraphs, (ii) Algorithms: wavefront sensing and control; post-processing of integral field spectrograph, polarimetric, and extended object data, and (iii) Validation of telescope and instrument models at high accuracy and precision. This white paper, submitted to the 2018 NAS Exoplanet Science Strategy call, describes the status of key CGI technologies and presents ways in which performance is likely to evolve as the CGI design matures.

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.

The Wide-Field Infrared Survey Telescope (WFIRST) is expected to launch in the mid-2020s. With its wide-field near-infrared (NIR) camera, it will survey the sky to unprecedented detail. As part of normal operations and as the result of multiple expected dedicated surveys, WFIRST will produce several relatively wide-field (tens of square degrees) deep (limiting magnitude of 28 or fainter) fields. In particular, a planned supernova survey is expected to image 3 deep fields in the LSST footprint roughly every 5 days over 2 years. Stacking all data, this survey will produce, over all WFIRST supernova fields in the LSST footprint, ~12-25 deg^2 and ~5-15 deg^2 regions to depths of ~28 mag and ~29 mag, respectively. We suggest LSST undertake mini-surveys that will match the WFIRST cadence and simultaneously observe the supernova survey fields during the 2-year WFIRST supernova survey, achieving a stacked depth similar to that of the WFIRST data. We also suggest additional observations of these same regions throughout the LSST survey to get deep images earlier, have long-term monitoring in the fields, and produce deeper images overall. These fields will provide a legacy for cosmology, extragalactic, and transient/variable science.

We review the measurements of dark energy enabled by observations of the Deep Drilling Fields and the optimization of survey design for cosmological measurements. This white paper is the result of efforts by the LSST DESC Observing Strategy Task Force (OSTF), which represents the entire collaboration, and aims to make recommendations on observing strategy for the DDFs that will benefit all cosmological analyses with LSST. It is accompanied by the DESC-WFD white paper (Lochner et al.). We argue for altering the nominal deep drilling plan to have $>6$ month seasons, interweaving $gri$ and $zy$ observations every 3 days with 2, 4, 8, 25, 4 visits in $grizy$, respectively. These recommendations are guided by metrics optimizing constraints on dark energy and mitigation of systematic uncertainties, including specific requirements on total number of visits after Y1 and Y10 for photometric redshifts (photo-$z$) and weak lensing systematics. We specify the precise locations for the previously-chosen LSST deep fields (ELAIS-S1, XMM-LSS, CDF-S, and COSMOS) and recommend Akari Deep Field South as the planned fifth deep field in order to synergize with Euclid and WFIRST. Our recommended DDF strategy uses $6.2\%$ of the LSST survey time. We briefly discuss synergy with white papers from other collaborations, as well as additional mini-surveys and Target-of-Opportunity programs that lead to better measurements of dark energy.

We present the latest version of the ray-tracing simulation code Skylens, which can be used to develop image simulations that reproduce strong lensing observations by any mass distribution with a high level of realism. Improvements of the code with respect to previous versions include the implementation of the multi-lens plane formalism, the use of denoised source galaxies from the Hubble eXtreme Deep Field, and the simulation of substructures in lensed arcs and images, based on a morphological analysis of bright nearby galaxies. Skylens can simulate observations with virtually any telescope. We present examples of space- and ground-based observations of a galaxy cluster through the Wide Field Channel on the Advanced Camera for Surveys of the Hubble Space Telescope, the Near Infrared Camera of the James Webb Space Telescope, the Wide Field Imager of the Wide Field Infrared Survey Telescope, the Hyper Suprime Camera of the Subaru telescope, and the Visible Imaging Channel of the Euclid space mission.

We present a community-led assessment of the solar system investigations achievable with NASA's next-generation space telescope, the Wide Field InfraRed Survey Telescope (WFIRST). WFIRST will provide imaging, spectroscopic, and coronagraphic capabilities from 0.43-2.0 $\mu$m and will be a potential contemporary and eventual successor to JWST. Surveys of irregular satellites and minor bodies are where WFIRST will excel with its 0.28 deg$^2$ field of view Wide Field Instrument (WFI). Potential ground-breaking discoveries from WFIRST could include detection of the first minor bodies orbiting in the Inner Oort Cloud, identification of additional Earth Trojan asteroids, and the discovery and characterization of asteroid binary systems similar to Ida/Dactyl. Additional investigations into asteroids, giant planet satellites, Trojan asteroids, Centaurs, Kuiper Belt Objects, and comets are presented. Previous use of astrophysics assets for solar system science and synergies between WFIRST, LSST, JWST, and the proposed NEOCam mission are discussed. We also present the case for implementation of moving target tracking, a feature that will benefit from the heritage of JWST and enable a broader range of solar system observations.

Balloon-borne telescopes present unique thermal design challenges which are a combination of those present for both space and ground telescopes. At altitudes of 35-40 km, convection effects are minimal and difficult to characterize. Radiation and conduction are the predominant heat transfer mechanisms reducing the thermal design options. For long duration flights payload mass is a function of power consumption making it an important optimization parameter. SuperBIT, or the Super-pressure Balloon-borne Imaging Telescope, aims to study weak lensing using a 0.5m modified Dall-Kirkham telescope capable of achieving 0.02" stability and capturing deep exposures from visible to near UV wavelengths. To achieve the theoretical stratospheric diffraction-limited resolution of 0.25", mirror deformation gradients must be kept to within 20nm. The thermal environment must thus be stable on time scales of an hour and the thermal gradients must be minimized on the telescope. SuperBIT plans to implement two types of parameter solvers; one to validate the thermal design and the other to tightly control the thermal environment.

Balloon-borne astronomy is a unique tool that allows for a level of image stability and significantly reduced atmospheric interference without the often prohibitive cost and long development time-scale that are characteristic of space-borne facility-class instruments. The Super-pressure Balloon-borne Imaging Telescope (SuperBIT) is a wide-field imager designed to provide 0.02" image stability over a 0.5 degree field-of-view for deep exposures within the visible-to-near-UV (300-900 um). As such, SuperBIT is a suitable platform for a wide range of balloon-borne observations, including solar and extrasolar planetary spectroscopy as well as resolved stellar populations and distant galaxies. We report on the overall payload design and instrumentation methodologies for SuperBIT as well as telescope and image stability results from two test flights. Prospects for the SuperBIT project are outlined with an emphasis on the development of a fully operational, three-month science flight from New Zealand in 2020.

Faint dwarf satellite galaxies are important as tracers of small-scale structure, but remain poorly characterized outside the Local Group, due to the difficulty of identifying them consistently at larger distances. We review a recently proposed method for estimating the average satellite population around a given sample of nearby bright galaxies, using a combination of size and magnitude cuts (to select low-redshift dwarf galaxies preferentially) and clustering measurements (to estimate the fraction of true satellites in the cut sample). We test this method using the high-precision photometric redshift catalog of the COSMOS survey, exploring the effect of specific cuts on the clustering signal. The most effective of the size-magnitude cuts considered recover the clustering signal around low-redshift primaries (z < 0.15) with about two-thirds of the signal and 80\%\ of the signal-to-noise ratio obtainable using the full COSMOS photometric redshifts. These cuts are also fairly efficient, with more than one third of the selected objects being clustered satellites. We conclude that structural selection represents a useful tool in characterizing dwarf populations to fainter magnitudes and/or over larger areas than are feasible with spectroscopic surveys.} In reviewing the low-redshift content of the COSMOS field, we also note the existence of several dozen objects that appear resolved or partially resolved in the HST imaging, and are confirmed to be local (at distances of $\sim$250 Mpc or less) by their photometric or spectroscopic redshifts. This underlines the potential for future space-based surveys to reveal local populations of intrinsically faint galaxies through imaging alone.

The "brighter-fatter" (BF) effect is a phenomenon (originally discovered in charge coupled devices) in which the size of the detector point spread function (PSF) increases with brightness. We present, for the first time, laboratory measurements demonstrating the existence of the effect in a Hawaii-2RG HgCdTe near infrared (NIR) detector. We use the Precision Projector Laboratory, a JPL facility for emulating astronomical observations with UV/VIS/NIR detectors, to project about 17,000 point sources onto the detector to stimulate the effect. After calibrating the detector for nonlinearity with flat-fields, we find evidence that charge is nonlinearly shifted from bright pixels to neighboring pixels during exposures of point sources, consistent with the existence of a BF-type effect. The Wide Field Infrared Survey Telescope (WFIRST) by NASA will use similar detectors to measure weak gravitational lensing from the shapes of hundreds of million of galaxies in the NIR. The WFIRST PSF size must be calibrated to approximately 0.1 percent to avoid biased inferences of dark matter and dark energy parameters; therefore further study and calibration of the BF effect in realistic images will be crucial.

Accurate shape measurements are essential to infer cosmological parameters from large area weak gravitational lensing studies. The compact diffraction-limited point-spread function (PSF) in space-based observations is greatly beneficial, but its chromaticity for a broad band observation can lead to new subtle effects that could hitherto be ignored: the PSF of a galaxy is no longer uniquely defined and spatial variations in the colours of galaxies result in biases in the inferred lensing signal. Taking Euclid as a reference, we show that this colourgradient bias (CG bias) can be quantified with high accuracy using available multi-colour Hubble Space Telescope (HST) data. In particular we study how noise in the HST observations might impact such measurements and find this to be negligible. We determine the CG bias using HST observations in the F606W and F814W filters and observe a correlation with the colour, in line with expectations, whereas the dependence with redshift is weak. The biases for individual galaxies are generally well below 1%, which may be reduced further using morphological information from the Euclid data. Our results demonstrate that CG bias should not be ignored, but it is possible to determine its amplitude with sufficient precision, so that it will not significantly bias the weak lensing measurements using Euclid data.

Cosmic acceleration is the most surprising cosmological discovery in many decades. Testing and distinguishing among possible explanations requires cosmological measurements of extremely high precision probing the full history of cosmic expansion and structure growth and, ideally, compare and contrast matter and relativistic tracers of the gravity potential. This program is one of the defining objectives of the Wide-Field Infrared Survey Telescope (WFIRST), as set forth in the New Worlds, New Horizons report (NWNH) in 2010. The WFIRST mission has the ability to improve these measurements by 1-2 orders of magnitude compared to the current state of the art, while simultaneously extending their redshift grasp, greatly improving control of systematic effects, and taking a unified approach to multiple probes that provide complementary physical information and cross-checks of cosmological results. We describe in this annual report the activities of the Science Investigation Team (SIT) "Cosmology with the High Latitude Survey (HLS)" during the year 2017. This team was selected by NASA in December 2015 in order to address the stringent challenges of the WFIRST dark energy (DE) program through the Project's formulation phase. This SIT has elected to jointly address Galaxy Redshift Survey, Weak Lensing and Cluster Growth and thus fully embrace the fact that the imaging and spectroscopic elements of the HLS will be realized as an integrated observing program, and they jointly impose requirements on performance and operations. WFIRST is designed to be able to deliver a definitive result on the origin of cosmic acceleration. It is not optimized for Figure of Merit sensitivity but for control of systematic uncertainties and for having multiple techniques each with multiple cross-checks. Our SIT work focuses on understanding the potential systematics in the WFIRST DE measurements.

The Wide Field Infrared Survey Telescope (WFIRST) was the top ranked large space mission in the 2010 New Worlds, New Horizons decadal survey, and it was formed by merging the science programs of 3 different mission concepts, including the Microlensing Planet Finder (MPF) concept (Bennett \etal\ 2010). The WFIRST science program (Spergel \etal\ 2015) consists of a general observer program, a wavefront controlled technology program, and two targeted science programs: a program to study dark energy, and a statistical census of exoplanets with a microlensing survey, which uses nearly one quarter of WFIRST's observing time in the current design reference mission. The New Worlds, New Horizons (decadal survey) midterm assessment summarizes the science case for the WFIRST exoplanet microlensing survey with this statement: "WFIRST's microlensing census of planets beyond 1 AU will perfectly complement Kepler's census of compact systems, and WFIRST will also be able to detect free-floating planets unbound from their parent stars\rlap."

The Wide Field Infrared Survey Telescope (WFIRST) was the top ranked large space mission in the 2010 New Worlds, New Horizons decadal survey, and it was formed by merging the science programs of 3 different mission concepts, including the Microlensing Planet Finder (MPF) concept (Bennett \etal\ 2010). The WFIRST science program (Spergel \etal\ 2015) consists of a general observer program, a wavefront controlled technology program, and two targeted science programs: a program to study dark energy, and a statistical census of exoplanets with a microlensing survey, which uses nearly one quarter of WFIRST's observing time in the current design reference mission. The New Worlds, New Horizons (decadal survey) midterm assessment summarizes the science case for the WFIRST exoplanet microlensing survey with this statement: "WFIRST's microlensing census of planets beyond 1 AU will perfectly complement Kepler's census of compact systems, and WFIRST will also be able to detect free-floating planets unbound from their parent stars\rlap."

As astronomical observations from space benefit from improved sensitivity, the effectiveness of scientific programs is becoming limited by systematics that often originate in poorly understood image sensor behavior. Traditional, bottom-up detector characterization methods provide one way to model underlying detector physics, and generate ever more faithful numerical simulations, but this approach is vulnerable to preconceptions and over-simplification. The alternative top-down approach is laboratory emulation, which enables observation, calibration, and analysis scenarios to be tested without relying on a complete understanding of the underlying physics. This complements detector characterization and simulation efforts by testing their validity. We describe a laboratory facility and experimental testbed that supports the emulation of a wide range of mission concepts such as gravitational weak lensing measurements by WFIRST and high precision spectrophotometry of transiting exoplanets by JWST. An Offner relay projects readily customizable "scenes" (e.g. stars, galaxies, spectra) with very low optical aberration over the full area of a typical optical or near infrared image sensor. f/8 and slower focal ratios may be selected, spanning those of most proposed space missions and approximating the point spread function (PSF) size of seeing limited ground based surveys. Diffraction limited PSFs are projected over a wide field of view and wavelength range to deliver highly predictable image patterns down to sub-pixel scales with stable intensity and fine motion control. The testbed enables realistic validation of detector performance on science-like images, which aids mission design and survey strategy, as well as targeted investigations of various detector effects.

Knowledge of the number density of H$\alpha$ emitting galaxies is vital for assessing the scientific impact of the Euclid and WFIRST missions. In this work we present predictions from a galaxy formation model, Galacticus, for the cumulative number counts of H$\alpha$-emitting galaxies. We couple Galacticus to three different dust attenuation methods and examine the counts using each method. A $\chi^2$ minimisation approach is used to compare the model predictions to observed galaxy counts and calibrate the dust parameters. We find that weak dust attenuation is required for the Galacticus counts to be broadly consistent with the observations, though the optimum dust parameters return large values for $\chi^2$, suggesting that further calibration of Galacticus is necessary. The model predictions are also consistent with observed estimates for the optical depth and the H$\alpha$ luminosity function. Finally we present forecasts for the redshift distributions and number counts for two Euclid-like and one WFIRST-like survey. For a Euclid-like survey with redshift range $0.9\leqslant z\leqslant 1.8$ and H$\alpha+{\rm [NII]}$ blended flux limit of $2\times 10^{-16}{\rm erg}\,{\rm s}^{-1}\,{\rm cm}^{-2}$ we predict a number density between 3900--4800 galaxies per square degree. For a WFIRST-like survey with redshift range $1\leqslant z\leqslant 2$ and blended flux limit of $1\times 10^{-16}{\rm erg}\,{\rm s}^{-1}\,{\rm cm}^{-2}$ we predict a number density between 10400--15200 galaxies per square degree.

Euclid and the Large Synoptic Survey Telescope (LSST) are poised to dramatically change the astronomy landscape early in the next decade. The combination of high cadence, deep, wide-field optical photometry from LSST with high resolution, wide-field optical photometry and near-infrared photometry and spectroscopy from Euclid will be powerful for addressing a wide range of astrophysical questions. We explore Euclid/LSST synergy, ignoring the political issues associated with data access to focus on the scientific, technical, and financial benefits of coordination. We focus primarily on dark energy cosmology, but also discuss galaxy evolution, transient objects, solar system science, and galaxy cluster studies. We concentrate on synergies that require coordination in cadence or survey overlap, or would benefit from pixel-level co-processing that is beyond the scope of what is currently planned, rather than scientific programs that could be accomplished only at the catalog level without coordination in data processing or survey strategies. We provide two quantitative examples of scientific synergies: the decrease in photo-z errors (benefitting many science cases) when high resolution Euclid data are used for LSST photo-z determination, and the resulting increase in weak lensing signal-to-noise ratio from smaller photo-z errors. We briefly discuss other areas of coordination, including high performance computing resources and calibration data. Finally, we address concerns about the loss of independence and potential cross-checks between the two missions and potential consequences of not collaborating.

In the last decade, astronomers have found a new type of supernova called `superluminous supernovae' (SLSNe) due to their high peak luminosity and long light-curves. These hydrogen-free explosions (SLSNe-I) can be seen to z~4 and therefore, offer the possibility of probing the distant Universe. We aim to investigate the possibility of detecting SLSNe-I using ESA's Euclid satellite, scheduled for launch in 2020. In particular, we study the Euclid Deep Survey (EDS) which will provide a unique combination of area, depth and cadence over the mission. We estimated the redshift distribution of Euclid SLSNe-I using the latest information on their rates and spectral energy distribution, as well as known Euclid instrument and survey parameters, including the cadence and depth of the EDS. We also applied a standardization method to the peak magnitudes to create a simulated Hubble diagram to explore possible cosmological constraints. We show that Euclid should detect approximately 140 high-quality SLSNe-I to z ~ 3.5 over the first five years of the mission (with an additional 70 if we lower our photometric classification criteria). This sample could revolutionize the study of SLSNe-I at z>1 and open up their use as probes of star-formation rates, galaxy populations, the interstellar and intergalactic medium. In addition, a sample of such SLSNe-I could improve constraints on a time-dependent dark energy equation-of-state, namely w(a), when combined with local SLSNe-I and the expected SN Ia sample from the Dark Energy Survey. We show that Euclid will observe hundreds of SLSNe-I for free. These luminous transients will be in the Euclid data-stream and we should prepare now to identify them as they offer a new probe of the high-redshift Universe for both astrophysics and cosmology.

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

A key goal of the Stage IV dark energy experiments Euclid, LSST and WFIRST is to measure the growth of structure with cosmic time from weak lensing analysis over large regions of the sky. Weak lensing cosmology will be challenging: in addition to highly accurate galaxy shape measurements, statistically robust and accurate photometric redshift (photo-z) estimates for billions of faint galaxies will be needed in order to reconstruct the three-dimensional matter distribution. Here we present an overview of and initial results from the Complete Calibration of the Color-Redshift Relation (C3R2) survey, designed specifically to calibrate the empirical galaxy color-redshift relation to the Euclid depth. These redshifts will also be important for the calibrations of LSST and WFIRST. The C3R2 survey is obtaining multiplexed observations with Keck (DEIMOS, LRIS, and MOSFIRE), the Gran Telescopio Canarias (GTC; OSIRIS), and the Very Large Telescope (VLT; FORS2 and KMOS) of a targeted sample of galaxies most important for the redshift calibration. We focus spectroscopic efforts on under-sampled regions of galaxy color space identified in previous work in order to minimize the number of spectroscopic redshifts needed to map the color-redshift relation to the required accuracy. Here we present the C3R2 survey strategy and initial results, including the 1283 high confidence redshifts obtained in the 2016A semester and released as Data Release 1.

We study the nonlinearity (NL) in the conversion from charge to voltage in infrared detectors (HXRG) for use in precision astronomy. We present laboratory measurements of the NL function of a H2RG detector and discuss the accuracy to which it would need to be calibrated in future space missions to perform cosmological measurements through the weak gravitational lensing technique. In addition, we present an analysis of archival data from the infrared H1RG detector of the Wide Field Camera 3 in the Hubble Space Telescope that provides evidence consistent with the existence of a sensor effect analogous to the brighter-fatter effect found in Charge-Coupled Devices. We propose a model in which this effect could be understood as shifts in the effective pixel boundaries, and discuss prospects of laboratory measurements to fully characterize this effect.

Balloon-borne astronomy is unique in that it allows for a level of image stability, resolution, and optical backgrounds that are comparable to space-borne systems due to greatly reduced atmospheric interference, but at a fraction of the cost and over a significantly reduced development time-scale. Instruments operating within visible-to-near-UV bands ($300$ - $900$ um) can achieve a theoretical diffraction limited resolution of $0.01"$ from the stratosphere ($35$ - $40$ km altitude) without the need for extensive adaptive optical systems required by ground-based systems. The {\it Superpressure Balloon-borne Imaging Telescope} ("SuperBIT") is a wide-field imager designed to achieve 0.02$"$ stability over a 0.5$^\circ$ field-of-view, for deep single exposures of up to 5 minutes. SuperBIT is thus well-suited for many astronomical observations, from solar or extrasolar planetary observations, to resolved stellar populations and distant galaxies (whether to study their morphology, evolution, or gravitational lensing by foreground mass). We report SuperBIT's design and implementation, emphasizing its two-stage real-time stabilization: telescope stability to $1$ - $2"$ at the telescope level (a goal surpassed during a test flight in September 2015) and image stability down to $0.02"$ via an actuated tip-tilt mirror in the optical path (to be tested during a flight in 2016). The project is progressing toward a fully operational, three month flight from New Zealand by 2018

The next generation weak lensing surveys (i.e., LSST, Euclid and WFIRST) will require exquisite control over systematic effects. In this paper, we address shear calibration and present the most realistic forecast to date for LSST/Euclid/WFIRST and CMB lensing from a stage 4 CMB experiment (CMB S4). We use the CosmoLike code to simulate a joint analysis of all the two-point functions of galaxy density, galaxy shear and CMB lensing convergence. We include the full Gaussian and non-Gaussian covariances and explore the resulting joint likelihood with Monte Carlo Markov Chains. We constrain shear calibration biases while simultaneously varying cosmological parameters, galaxy biases and photometric redshift uncertainties. We find that CMB lensing from CMB S4 enables the calibration of the shear biases down to 0.2% - 3% in 10 tomographic bins for LSST (below the ~0.5% requirements in most tomographic bins), down to 0.4% - 2.4% in 10 bins for Euclid and 0.6% - 3.2% in 10 bins for WFIRST. For a given lensing survey, the method works best at high redshift where shear calibration is otherwise most challenging. This self-calibration is robust to Gaussian photometric redshift uncertainties and to a reasonable level of intrinsic alignment. It is also robust to changes in the beam and the effectiveness of the component separation of the CMB experiment, and slowly dependent on its depth, making it possible with third generation CMB experiments such as AdvACT and SPT-3G, as well as the Simons Observatory.

Unlike optical CCDs, near-infrared detectors, which are based on CMOS hybrid readout technology, typically suffer from electrical crosstalk between the pixels. The interpixel capacitance (IPC) responsible for the crosstalk affects the point-spread function (PSF) of the telescope, increasing the size and modifying the shape of all objects in the images while correlating the Poisson noise. Upcoming weak lensing surveys that use these detectors, such as WFIRST, place stringent requirements on the PSF size and shape (and the level at which these are known), which in turn must be translated into requirements on IPC. To facilitate this process, we present a first study of the effect of IPC on WFIRST PSF sizes and shapes. Realistic PSFs are forward-simulated from physical principles for each WFIRST bandpass. We explore how the PSF size and shape depends on the range of IPC coupling with pixels that are connected along an edge or corner; for the expected level of IPC in WFIRST, IPC increases the PSF sizes by $\sim$5\%. We present a linear fitting formula that describes the uncertainty in the PSF size or shape due to uncertainty in the IPC, which could arise for example due to unknown time evolution of IPC as the detectors age or due to spatial variation of IPC across the detector. We also study of the effect of a small anisotropy in the IPC, which further modifies the PSF shapes. Our results are a first, critical step in determining the hardware and characterization requirements for the detectors used in the WFIRST survey.

SPHEREx is a proposed SMEX mission selected for Phase A. SPHEREx will carry out the first all-sky spectral survey and provide for every 6.2" pixel a spectra between 0.75 and 4.18 $\mu$m [with R$\sim$41.4] and 4.18 and 5.00 $\mu$m [with R$\sim$135]. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. It is readily apparent, however, that many other questions in astrophysics and planetary sciences could be addressed with the SPHEREx data. The SPHEREx team convened a community workshop in February 2016, with the intent of enlisting the aid of a larger group of scientists in defining these questions. This paper summarizes the rich and varied menu of investigations that was laid out. It includes studies of the composition of main belt and Trojan/Greek asteroids; mapping the zodiacal light with unprecedented spatial and spectral resolution; identifying and studying very low-metallicity stars; improving stellar parameters in order to better characterize transiting exoplanets; studying aliphatic and aromatic carbon-bearing molecules in the interstellar medium; mapping star formation rates in nearby galaxies; determining the redshift of clusters of galaxies; identifying high redshift quasars over the full sky; and providing a NIR spectrum for most eROSITA X-ray sources. All of these investigations, and others not listed here, can be carried out with the nominal all-sky spectra to be produced by SPHEREx. In addition, the workshop defined enhanced data products and user tools which would facilitate some of these scientific studies. Finally, the workshop noted the high degrees of synergy between SPHEREx and a number of other current or forthcoming programs, including JWST, WFIRST, Euclid, GAIA, K2/Kepler, TESS, eROSITA and LSST.

Weak gravitational lensing (WL) is one of the most powerful techniques to learn about the dark sector of the universe. To extract the WL signal from astronomical observations, galaxy shapes must be measured and corrected for the point spread function (PSF) of the imaging system with extreme accuracy. Future WL missions (such as the Wide-Field Infrared Survey Telescope, WFIRST) will use a family of hybrid nearinfrared CMOS detectors (HAWAII-4RG) that are untested for accurate WL measurements. Like all image sensors, these devices are subject to conversion gain nonlinearities (voltage response to collected photo-charge) that bias the shape and size of bright objects such as reference stars that are used in PSF determination. We study this type of detector nonlinearity (NL) and show how to derive requirements on it from WFIRST PSF size and ellipticity requirements. We simulate the PSF optical profiles expected for WFIRST and measure the fractional error in the PSF size and the absolute error in the PSF ellipticity as a function of star magnitude and the NL model. For our nominal NL model (a quadratic correction), we find that, uncalibrated, NL can induce an error of 0.01 (fractional size) and 0.00175 (absolute ellipticity error) in the H158 bandpass for the brightest unsaturated stars in WFIRST. In addition, our simulations show that to limit the bias of the size and ellipticity errors in the H158 band to approximately 10% of the estimated WFIRST error budget, the parameter of our quadratic NL model must be calibrated to about 1% and 2.4%, respectively. We present a fitting formula that can be used to estimate WFIRST detector NL requirements once a true PSF error budget is established.

We present a weak lensing study of the galaxy cluster IDCS J1426.5+3508 at $z=1.75$, which is the highest redshift strong lensing cluster known and the most distant cluster for which a weak lensing analysis has been undertaken. Using F160W, F814W, and F606W observations with the Hubble Space Telescope, we detect tangential shear at $2\sigma$ significance. Fitting a Navarro-Frenk-White mass profile to the shear with a theoretical median mass-concentration relation, we derive a mass $M_{200,\mathrm{crit}}=2.3^{+2.1}_{-1.4}\times10^{14}$ M$_{\odot}$. This mass is consistent with previous mass estimates from the Sunyaev-Zel'dovich (SZ) effect, X-ray, and strong lensing. The cluster lies on the local SZ-weak lensing mass scaling relation observed at low redshift, indicative of minimal evolution in this relation.

Calibrating the photometric redshifts of >10^9 galaxies for upcoming weak lensing cosmology experiments is a major challenge for the astrophysics community. The path to obtaining the required spectroscopic redshifts for training and calibration is daunting, given the anticipated depths of the surveys and the difficulty in obtaining secure redshifts for some faint galaxy populations. Here we present an analysis of the problem based on the self-organizing map, a method of mapping the distribution of data in a high-dimensional space and projecting it onto a lower-dimensional representation. We apply this method to existing photometric data from the COSMOS survey selected to approximate the anticipated Euclid weak lensing sample, enabling us to robustly map the empirical distribution of galaxies in the multidimensional color space defined by the expected Euclid filters. Mapping this multicolor distribution lets us determine where - in galaxy color space - redshifts from current spectroscopic surveys exist and where they are systematically missing. Crucially, the method lets us determine whether a spectroscopic training sample is representative of the full photometric space occupied by the galaxies in a survey. We explore optimal sampling techniques and estimate the additional spectroscopy needed to map out the color-redshift relation, finding that sampling the galaxy distribution in color space in a systematic way can efficiently meet the calibration requirements. While the analysis presented here focuses on the Euclid survey, similar analysis can be applied to other surveys facing the same calibration challenge, such as DES, LSST, and WFIRST.

Radiation damage to space-based Charge-Coupled Device (CCD) detectors creates defects which result in an increasing Charge Transfer Inefficiency (CTI) that causes spurious image trailing. Most of the trailing can be corrected during post-processing, by modelling the charge trapping and moving electrons back to where they belong. However, such correction is not perfect -- and damage is continuing to accumulate in orbit. To aid future development, we quantify the limitations of current approaches, and determine where imperfect knowledge of model parameters most degrade measurements of photometry and morphology. As a concrete application, we simulate $1.5\times10^{9}$ "worst case" galaxy and $1.5\times10^{8}$ star images to test the performance of the Euclid visual instrument detectors. There are two separable challenges: If the model used to correct CTI is perfectly the same as that used to add CTI, $99.68$ % of spurious ellipticity is corrected in our setup. This is because readout noise is not subject to CTI, but gets over-corrected during correction. Second, if we assume the first issue to be solved, knowledge of the charge trap density within $\Delta\rho/\rho\!=\!(0.0272\pm0.0005)$ %, and the characteristic release time of the dominant species to be known within $\Delta\tau/\tau\!=\!(0.0400\pm0.0004)$ % will be required. This work presents the next level of definition of in-orbit CTI calibration procedures for Euclid.

We present first results from the third GRavitational lEnsing Accuracy Testing (GREAT3) challenge, the third in a sequence of challenges for testing methods of inferring weak gravitational lensing shear distortions from simulated galaxy images. GREAT3 was divided into experiments to test three specific questions, and included simulated space- and ground-based data with constant or cosmologically-varying shear fields. The simplest (control) experiment included parametric galaxies with a realistic distribution of signal-to-noise, size, and ellipticity, and a complex point spread function (PSF). The other experiments tested the additional impact of realistic galaxy morphology, multiple exposure imaging, and the uncertainty about a spatially-varying PSF; the last two questions will be explored in Paper II. The 24 participating teams competed to estimate lensing shears to within systematic error tolerances for upcoming Stage-IV dark energy surveys, making 1525 submissions overall. GREAT3 saw considerable variety and innovation in the types of methods applied. Several teams now meet or exceed the targets in many of the tests conducted (to within the statistical errors). We conclude that the presence of realistic galaxy morphology in simulations changes shear calibration biases by $\sim 1$ per cent for a wide range of methods. Other effects such as truncation biases due to finite galaxy postage stamps, and the impact of galaxy type as measured by the S\'{e}rsic index, are quantified for the first time. Our results generalize previous studies regarding sensitivities to galaxy size and signal-to-noise, and to PSF properties such as seeing and defocus. Almost all methods' results support the simple model in which additive shear biases depend linearly on PSF ellipticity.

This report describes the 2014 study by the Science Definition Team (SDT) of the Wide-Field Infrared Survey Telescope (WFIRST) mission. It is a space observatory that will address the most compelling scientific problems in dark energy, exoplanets and general astrophysics using a 2.4-m telescope with a wide-field infrared instrument and an optical coronagraph. The Astro2010 Decadal Survey recommended a Wide Field Infrared Survey Telescope as its top priority for a new large space mission. As conceived by the decadal survey, WFIRST would carry out a dark energy science program, a microlensing program to determine the demographics of exoplanets, and a general observing program utilizing its ultra wide field. In October 2012, NASA chartered a Science Definition Team (SDT) to produce, in collaboration with the WFIRST Study Office at GSFC and the Program Office at JPL, a Design Reference Mission (DRM) for an implementation of WFIRST using one of the 2.4-m, Hubble-quality telescope assemblies recently made available to NASA. This DRM builds on the work of the earlier WFIRST SDT, reported by Green et al. (2012) and the previous WFIRST-2.4 DRM, reported by Spergel et. (2013). The 2.4-m primary mirror enables a mission with greater sensitivity and higher angular resolution than the 1.3-m and 1.1-m designs considered previously, increasing both the science return of the primary surveys and the capabilities of WFIRST as a Guest Observer facility. The addition of an on-axis coronagraphic instrument to the baseline design enables imaging and spectroscopic studies of planets around nearby stars.

Current constraints on dark matter density profiles from weak lensing are typically limited to radial scales greater than 50-100 kpc. In this paper, we explore the possibility of probing the very inner regions of galaxy/halo density profiles by measuring stacked weak lensing on scales of only a few tens of kpc. Our forecasts focus on scales smaller than the equality radius (Req) where the stellar component and the dark matter component contribute equally to the lensing signal. We compute the evolution of Req as a function of lens stellar mass and redshift and show that Req=7-34 kpc for galaxies with the stellar mass of 10^{9.5}-10^{11.5} solar masses. Unbiased shear measurements will be challenging on these scales. We introduce a simple metric to quantify how many source galaxies overlap with their neighbours and for which shear measurements will be challenging. Rejecting source galaxies with close-by companions results in about a 20 per cent decrease in the overall source density. Despite this decrease, we show that Euclid and WFIRST will be able to constrain galaxy/halo density profiles at Req with signal-to-noise ratio >20 for the stellar mass of >10^{10} solar masses. Weak lensing measurements at Req, in combination with stellar kinematics on smaller scales, will be a powerful means by which to constrain both the inner slope of the dark matter density profile as well as the mass and redshift dependence of the stellar initial mass function.

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.

Understanding the relationship between galaxies hosting active galactic nuclei (AGN) and the dark matter halos in which they reside is key to constraining how black-hole fueling is triggered and regulated. Previous efforts have relied on simple halo mass estimates inferred from clustering, weak gravitational lensing, or halo occupation distribution modeling. In practice, these approaches remain uncertain because AGN, no matter how they are identified, potentially live a wide range of halo masses with an occupation function whose general shape and normalization are poorly known. In this work, we show that better constraints can be achieved through a rigorous comparison of the clustering, lensing, and cross-correlation signals of AGN hosts to a fiducial stellar-to-halo mass relation (SHMR) derived for all galaxies. Our technique exploits the fact that the global SHMR can be measured with much higher accuracy than any statistic derived from AGN samples alone. Using 382 moderate luminosity X-ray AGN at z<1 from the COSMOS field, we report the first measurements of weak gravitational lensing from an X-ray selected sample. Comparing this signal to predictions from the global SHMR, we find that, contrary to previous results, most X-ray AGN do not live in medium size groups ---nearly half reside in relatively low mass halos with Mh~10^12.5 Msun. The AGN occupation function is well described by the same form derived for all galaxies but with a lower normalization---the fraction of halos with AGN in our sample is a few percent. By highlighting the relatively "normal" way in which moderate luminosity X-ray AGN hosts occupy halos, our results suggest that the environmental signature of distinct fueling modes for luminous QSOs compared to moderate luminosity X-ray AGN is less obvious than previously claimed.

We present the results of a numerical study based on the analysis of the MUSIC-2 simulations, aimed at estimating the expected concentration-mass relation for the CLASH cluster sample. We study nearly 1400 halos simulated at high spatial and mass resolution, which were projected along many lines-of-sight each. We study the shape of both their density and surface-density profiles and fit them with a variety of radial functions, including the Navarro-Frenk-White, the generalised Navarro-Frenk-White, and the Einasto density profiles. We derive concentrations and masses from these fits and investigate their distributions as a function of redshift and halo relaxation. We use the X-ray image simulator X-MAS to produce simulated Chandra observations of the halos and we use them to identify objects resembling the X-ray morphologies and masses of the clusters in the CLASH X-ray selected sample. We also derive a concentration-mass relation for strong-lensing clusters. We find that the sample of simulated halos which resemble the X-ray morphology of the CLASH clusters is composed mainly by relaxed halos, but it also contains a significant fraction of un-relaxed systems. For such a sample we measure an average 2D concentration which is ~11% higher than found for the full sample of simulated halos. After accounting for projection and selection effects, the average NFW concentrations of CLASH clusters are expected to be intermediate between those predicted in 3D for relaxed and super-relaxed halos. Matching the simulations to the individual CLASH clusters on the basis of the X-ray morphology, we expect that the NFW concentrations recovered from the lensing analysis of the CLASH clusters are in the range [3-6], with an average value of 3.87 and a standard deviation of 0.61. Simulated halos with X-ray morphologies similar to those of the CLASH clusters are affected by a modest orientation bias.

We present a new determination of the concentration-mass relation for galaxy clusters based on our comprehensive lensing analysis of 19 X-ray selected galaxy clusters from the Cluster Lensing and Supernova Survey with Hubble (CLASH). Our sample spans a redshift range between 0.19 and 0.89. We combine weak lensing constraints from the Hubble Space Telescope (HST) and from ground-based wide field data with strong lensing constraints from HST. The result are reconstructions of the surface-mass density for all CLASH clusters on multi-scale grids. Our derivation of NFW parameters yields virial masses between 0.53 x 10^15 and 1.76 x 10^15 M_sol/h and the halo concentrations are distributed around c_200c ~ 3.7 with a 1-sigma significant negative trend with cluster mass. We find an excellent 4% agreement between our measured concentrations and the expectation from numerical simulations after accounting for the CLASH selection function based on X-ray morphology. The simulations are analyzed in 2D to account for possible biases in the lensing reconstructions due to projection effects. The theoretical concentration-mass (c-M) relation from our X-ray selected set of simulated clusters and the c-M relation derived directly from the CLASH data agree at the 90% confidence level.

In this report, we present a wide variety of ways in which information from multiple probes of dark energy may be combined to obtain additional information not accessible when they are considered separately. Fundamentally, because all major probes are affected by the underlying distribution of matter in the regions studied, there exist covariances between them that can provide information on cosmology. Combining multiple probes allows for more accurate (less contaminated by systematics) and more precise (since there is cosmological information encoded in cross-correlation statistics) measurements of dark energy. The potential of cross-correlation methods is only beginning to be realized. By bringing in information from other wavelengths, the capabilities of the existing probes of dark energy can be enhanced and systematic effects can be mitigated further. We present a mixture of work in progress and suggestions for future scientific efforts. Given the scope of future dark energy experiments, the greatest gains may only be realized with more coordination and cooperation between multiple project teams; we recommend that this interchange should begin sooner, rather than later, to maximize scientific gains.

The GRavitational lEnsing Accuracy Testing 3 (GREAT3) challenge is the third in a series of image analysis challenges, with a goal of testing and facilitating the development of methods for analyzing astronomical images that will be used to measure weak gravitational lensing. This measurement requires extremely precise estimation of very small galaxy shape distortions, in the presence of far larger intrinsic galaxy shapes and distortions due to the blurring kernel caused by the atmosphere, telescope optics, and instrumental effects. The GREAT3 challenge is posed to the astronomy, machine learning, and statistics communities, and includes tests of three specific effects that are of immediate relevance to upcoming weak lensing surveys, two of which have never been tested in a community challenge before. These effects include realistically complex galaxy models based on high-resolution imaging from space; spatially varying, physically-motivated blurring kernel; and combination of multiple different exposures. To facilitate entry by people new to the field, and for use as a diagnostic tool, the simulation software for the challenge is publicly available, though the exact parameters used for the challenge are blinded. Sample scripts to analyze the challenge data using existing methods will also be provided. See http://great3challenge.info and http://great3.projects.phys.ucl.ac.uk/leaderboard/ for more information.

The "Airships: A New Horizon for Science" study at the Keck Institute for Space Studies investigated the potential of a variety of airships currently operable or under development to serve as observatories and science instrumentation platforms for a range of space, atmospheric, and Earth science. The participants represent a diverse cross-section of the aerospace sector, NASA, and academia. Over the last two decades, there has been wide interest in developing a high altitude, stratospheric lighter-than-air (LTA) airship that could maneuver and remain in a desired geographic position (i.e., "station-keeping") for weeks, months or even years. Our study found considerable scientific value in both low altitude (< 40 kft) and high altitude (> 60 kft) airships across a wide spectrum of space, atmospheric, and Earth science programs. Over the course of the study period, we identified stratospheric tethered aerostats as a viable alternative to airships where station-keeping was valued over maneuverability. By opening up the sky and Earth's stratospheric horizon in affordable ways with long-term flexibility, airships allow us to push technology and science forward in a project-rich environment that complements existing space observatories as well as aircraft and high-altitude balloon missions.

Extremely accurate shape measurements of galaxy images are needed to probe dark energy properties with weak gravitational lensing surveys. To increase survey area with a fixed observing time and pixel count, images from surveys such as the Wide Field Infrared Survey Telescope (WFIRST) or Euclid will necessarily be undersampled and therefore distorted by aliasing. Oversampled, unaliased images can be obtained by combining multiple, dithered exposures of the same source with a suitable reconstruction algorithm. Any such reconstruction must minimally distort the reconstructed images for weak lensing analyses to be unbiased. In this paper, we use the IMage COMbination (IMCOM) algorithm of Rowe, Hirata, and Rhodes to investigate the effect of image combination on shape measurements (size and ellipticity). We simulate dithered images of sources with varying amounts of ellipticity and undersampling, reconstruct oversampled output images from them using IMCOM, and measure shape distortions in the output. Our simulations show that IMCOM creates no significant distortions when the relative offsets between dithered images are precisely known. Distortions increase with the uncertainty in those offsets but become problematic only with relatively poor astrometric precision. E.g. for images similar to those from the Astrophysics Focused Telescope Asset (AFTA) implementation of WFIRST, combining eight undersampled images (sampling ratio Q=1) with highly pessimistic uncertainty in astrometric registration (\sigma_d~10^{-3} pixels) yields an RMS shear error of O(10^{-4}). Our analysis pipeline is adapted from that of the Precision Projector Laboratory -- a joint project between NASA Jet Propulsion Laboratory and Caltech which characterizes image sensors using laboratory emulations of astronomical data.

The American Physical Society's Division of Particles and Fields initiated a long-term planning exercise over 2012-13, with the goal of developing the community's long term aspirations. The sub-group "Dark Energy and CMB" prepared a series of papers explaining and highlighting the physics that will be studied with large galaxy surveys and cosmic microwave background experiments. This paper summarizes the findings of the other papers, all of which have been submitted jointly to the arXiv.

Large sets of objects with spectroscopic redshift measurements will be needed for imaging dark energy experiments to achieve their full potential, serving two goals:_training_, i.e., the use of objects with known redshift to develop and optimize photometric redshift algorithms; and_calibration_, i.e., the characterization of moments of redshift (or photo-z error) distributions. Better training makes cosmological constraints from a given experiment stronger, while highly-accurate calibration is needed for photo-z systematics not to dominate errors. In this white paper, we investigate the required scope of spectroscopic datasets which can serve both these purposes for ongoing and next-generation dark energy experiments, as well as the time required to obtain such data with instruments available in the next decade. Large time allocations on kilo-object spectrographs will be necessary, ideally augmented by infrared spectroscopy from space. Alternatively, precision calibrations could be obtained by measuring cross-correlation statistics using samples of bright objects from a large baryon acoustic oscillation experiment such as DESI. We also summarize the additional work on photometric redshift methods needed to prepare for ongoing and future dark energy experiments.

The discovery of cosmic acceleration has inspired ambitious experimental and observational efforts to understand its origin. Many of these take the form of large astronomical surveys, sometimes using new, special-purpose instrumentation, and in some cases entirely new facilities. This Report summarizes some of the major ongoing and planned dark energy experiments, focusing on those in which the U.S. community has a leading or significant supporting role. It provides background information for other Reports from the "Dark Energy and CMB" working group of the APS Division of Particles and Fields long-term planning exercise. We provide 1-2 page summaries of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) and its SDSS-IV successor eBOSS, the Dark Energy Survey (DES), the Subaru Hyper-Suprime Camera (HSC) and Prime Focus Spectrograph (PFS), the Hobby-Eberley Telescope Dark Energy Experiment (HEDTEX), the Dark Energy Spectroscopic Instrument (DESI), the Large Synpotic Survey Telescope (LSST), and the Euclid and AFTA-WFIRST (Wide Field Infrared Survey Telescope) space missions. Over the next one to two decades, these experiments seek to improve the precision of measurements of the cosmic expansion history and cosmic structure growth by factors of 10-100, while simultaneously increasing their redshift range and tightening control and cross-checks of sytematic effects. These measurements will dramatically improve the empirical constraints on theoretical explanations of cosmic acceleration and could lead to breakthroughs in our understanding of the nature of gravity and spacetime or the forms of energy present in the universe.

Weak gravitational lensing observations are a key science driver for the NASA Wide Field Infrared Survey Telescope (WFIRST). To validate the performance of the WFIRST infrared detectors, we have performed a laboratory emulation of weak gravitational lensing measurements. Our experiments used a custom precision projector system to image a target mask composed of a grid of pinholes, emulating stellar point sources, onto a 1.7 micron cut-off Teledyne HgCdTe/H2RG detector. We used a 880nm LED illumination source and f/22 pupil stop to produce undersampled point spread functions similar to those expected from WFIRST. We also emulated the WFIRST image reconstruction strategy, using the IMage COMbination (IMCOM) algorithm to derive oversampled images from dithered, undersampled input images. We created shear maps for this data and computed shear correlation functions to mimic a real weak lensing analysis. After removing only 2nd order polynomial fits to the shear maps, we found that the correlation functions could be reduced to O(10^-6). This places a conservative upper limit on the detector-induced bias to the correlation function (under our test conditions). This bias is two orders of magnitude lower than the expected weak lensing signal. Restricted to scales relevant to dark energy analyses (sky separations > 0.5 arcmin), the bias is O(10^-7): comparable to the requirement for future weak lensing missions to avoid biasing cosmological parameter estimates. Our experiment will need to be upgraded and repeated under different configurations to fully characterize the shape measurement performance of WFIRST IR detectors.

The scaling between X-ray observables and mass for galaxy clusters and groups is instrumental for cluster based cosmology and an important probe for the thermodynamics of the intracluster gas. We calibrate a scaling relation between the weak lensing mass and X-ray spectroscopic temperature for 10 galaxy groups in the COSMOS field, combined with 55 higher mass clusters from the literature. The COSMOS data includes HST imaging and redshift measurements of 46 source galaxies per square arcmin, enabling us to perform unique weak lensing measurements of low mass systems. Our sample extends the mass range of the lensing calibrated M-T relation an order of magnitude lower than any previous study, resulting in a power-law slope of 1.48$^{+0.13}_{-0.09}$. The slope is consistent with the self-similar model, predictions from simulations, and observations of clusters. However, X-ray observations relying on mass measurements derived under the assumption of hydrostatic equilibrium have indicated that masses at group scales are lower than expected. Both simulations and observations suggest that hydrostatic mass measurements can be biased low. Our external weak lensing masses provides the first observational support for hydrostatic mass bias at group level, showing an increasing bias with decreasing temperature and reaching a level of 30-50% at 1 keV.

We present a simulation analysis of weak gravitational lensing flexion and shear measurement using shapelet decomposition, and identify differences between flexion and shear measurement noise in deep survey data. Taking models of galaxies from the Hubble Space Telescope Ultra Deep Field (HUDF) and applying a correction for the HUDF point spread function we generate lensed simulations of deep, optical imaging data from Hubble's Advanced Camera for Surveys (ACS), with realistic galaxy morphologies. We find that flexion and shear estimates differ in our measurement pipeline: whereas intrinsic galaxy shape is typically the dominant contribution to noise in shear estimates, pixel noise due to finite photon counts and detector read noise is a major contributor to uncertainty in flexion estimates, across a broad range of galaxy signal-to-noise. This pixel noise also increases more rapidly as galaxy signal-to-noise decreases than is found for shear estimates. We provide simple power law fitting functions for this behaviour, for both flexion and shear, allowing the effect to be properly accounted for in future forecasts for flexion measurement. Using the simulations we also quantify the systematic biases of our shapelet flexion and shear measurement pipeline for deep Hubble data sets such as Galaxy Evolution from Morphology and SEDs, Space Telescope A901/902 Galaxy Evolution Survey or the Cosmic Evolution Survey. Flexion measurement biases are found to be significant but consistent with previous studies.

We use measurements of the stellar mass function, galaxy clustering, and galaxy-galaxy lensing within the COSMOS survey to constrain the stellar-to-halo mass relation (SHMR) of star forming and quiescent galaxies over the redshift range z=[0.2,1.0]. For massive galaxies, M*>~10^10.6 Msol, our results indicate that star-forming galaxies grow proportionately as fast as their dark matter halos while quiescent galaxies are outpaced by dark matter growth. At lower masses, there is minimal difference in the SHMRs, implying that the majority low-mass quiescent galaxies have only recently been quenched of their star formation. Our analysis also affords a breakdown of all COSMOS galaxies into the relative numbers of central and satellite galaxies for both populations. At z=1, satellite galaxies dominate the red sequence below the knee in the stellar mass function. But the number of quiescent satellites exhibits minimal redshift evolution; all evolution in the red sequence is due to low-mass central galaxies being quenched of their star formation. At M*~10^10 Msol, the fraction of central galaxies on the red sequence increases by a factor of ten over our redshift baseline, while the fraction of quenched satellite galaxies at that mass is constant with redshift. We define a "migration rate" to the red sequence as the time derivative of the passive galaxy abundances. We find that the migration rate of central galaxies to the red sequence increases by nearly an order of magnitude from z=1 to z=0. These results imply that the efficiency of quenching star formation for centrals is increasing with cosmic time, while the mechanisms that quench the star formation of satellite galaxies in groups and clusters is losing efficiency.

The Astro2010 Decadal Survey recommended a Wide Field Infrared Survey Telescope (WFIRST) as its top priority for a new large space mission. The report of the WFIRST-AFTA Science Definition Team (SDT) presents a Design Reference Mission for WFIRST that employs one of the 2.4-m, Hubble-quality mirror assemblies recently made available to NASA. The 2.4-m primary mirror enables a mission with greater sensitivity and higher angular resolution than the smaller aperture designs previously considered for WFIRST, increasing both the science return of the primary surveys and the capabilities of WFIRST as a Guest Observer facility. The option of adding an on-axis, coronagraphic instrument would enable imaging and spectroscopic studies of planets around nearby stars. This short article, produced as a companion to the SDT report, summarizes the key points of the WFIRST-2.4 DRM. It highlights the remarkable opportunity that the 2.4-m telescope affords for advances in many fields of astrophysics and cosmology, including dark energy, the demographics and characterization of exoplanets, the evolution of galaxies and quasars, and the stellar populations of the Milky Way and its neighbors.

The Astro2010 Decadal Survey recommended a Wide Field Infrared Survey Telescope (WFIRST) as its top priority for a new large space mission. As conceived by the decadal survey, WFIRST would carry out a dark energy science program, a microlensing program to determine the demographics of exoplanets, and a general observing program utilizing its ultra wide field. In October 2012, NASA chartered a Science Definition Team (SDT) to produce, in collaboration with the WFIRST Project Office at GSFC and the Program Office at JPL, a Design Reference Mission (DRM) for an implementation of WFIRST using one of the 2.4-m, Hubble-quality mirror assemblies recently made available to NASA. This DRM builds on the work of the earlier WFIRST SDT, reported by Green et al. (2012). The 2.4-m primary mirror enables a mission with greater sensitivity and higher angular resolution than the 1.3-m and 1.1-m designs considered previously, increasing both the science return of the primary surveys and the capabilities of WFIRST as a Guest Observer facility. The option of adding an on-axis, coronagraphic instrument would enable imaging and spectroscopic studies of planets around nearby stars. This document presents the final report of the SDT.

This paper describes the definition of a typical next-generation space-based weak gravitational lensing experiment. We first adopt a set of top-level science requirements from the literature, based on the scale and depth of the galaxy sample, and the avoidance of systematic effects in the measurements which would bias the derived shear values. We then identify and categorise the contributing factors to the systematic effects, combining them with the correct weighting, in such a way as to fit within the top-level requirements. We present techniques which permit the performance to be evaluated and explore the limits at which the contributing factors can be managed. Besides the modelling biases resulting from the use of weighted moments, the main contributing factors are the reconstruction of the instrument point spread function (PSF), which is derived from the stellar images on the image, and the correction of the charge transfer inefficiency (CTI) in the CCD detectors caused by radiation damage.

Large-scale structures (LSS) out to z $< 3.0$ are measured in the Cosmic Evolution Survey (COSMOS) using extremely accurate photometric redshifts (photoz). The Ks-band selected sample (from Ultra-Vista) is comprised of 155,954 galaxies. Two techniques -- adaptive smoothing and Voronoi tessellation -- are used to estimate the environmental densities within 127 redshift slices. Approximately 250 statistically significant overdense structures are identified out to z $= 3.0$ with shapes varying from elongated filamentary structures to more circularly symmetric concentrations. We also compare the densities derived for COSMOS with those based on semi-analytic predictions for a $\Lambda$CDM simulation and find excellent overall agreement between the mean densities as a function of redshift and the range of densities. The galaxy properties (stellar mass, spectral energy distributions (SEDs) and star formation rates (SFRs)) are strongly correlated with environmental density and redshift, particularly at z $< 1.0 - 1.2$. Classifying the spectral type of each galaxy using the rest-frame b-i color (from the photoz SED fitting), we find a strong correlation of early type galaxies (E-Sa) with high density environments, while the degree of environmental segregation varies systematically with redshift out to z $\sim 1.3$. In the highest density regions, 80% of the galaxies are early types at z=0.2 compared to only 20% at z = 1.5. The SFRs and the star formation timescales exhibit clear environmental correlations. At z $> 0.8$, the star formation rate density (SFRD) is uniformly distributed over all environmental density percentiles, while at lower redshifts the dominant contribution is shifted to galaxies in lower density environments.

In this paper we present results from the weak lensing shape measurement GRavitational lEnsing Accuracy Testing 2010 (GREAT10) Galaxy Challenge. This marks an order of magnitude step change in the level of scrutiny employed in weak lensing shape measurement analysis. We provide descriptions of each method tested and include 10 evaluation metrics over 24 simulation branches. GREAT10 was the first shape measurement challenge to include variable fields; both the shear field and the Point Spread Function (PSF) vary across the images in a realistic manner. The variable fields enable a variety of metrics that are inaccessible to constant shear simulations including a direct measure of the impact of shape measurement inaccuracies, and the impact of PSF size and ellipticity, on the shear power spectrum. To assess the impact of shape measurement bias for cosmic shear we present a general pseudo-Cl formalism, that propagates spatially varying systematics in cosmic shear through to power spectrum estimates. We also show how one-point estimators of bias can be extracted from variable shear simulations. The GREAT10 Galaxy Challenge received 95 submissions and saw a factor of 3 improvement in the accuracy achieved by shape measurement methods. The best methods achieve sub-percent average biases. We find a strong dependence in accuracy as a function of signal-to-noise, and indications of a weak dependence on galaxy type and size. Some requirements for the most ambitious cosmic shear experiments are met above a signal-to-noise ratio of 20. These results have the caveat that the simulated PSF was a ground-based PSF. Our results are a snapshot of the accuracy of current shape measurement methods and are a benchmark upon which improvement can continue. This provides a foundation for a better understanding of the strengths and limitations of shape measurement methods.

The spatial variation of the colour of a galaxy may introduce a bias in the measurement of its shape if the PSF profile depends on wavelength. We study how this bias depends on the properties of the PSF and the galaxies themselves. The bias depends on the scales used to estimate the shape, which may be used to optimise methods to reduce the bias. Here we develop a general approach to quantify the bias. Although applicable to any weak lensing survey, we focus on the implications for the ESA Euclid mission. Based on our study of synthetic galaxies we find that the bias is a few times 10^-3 for a typical galaxy observed by Euclid. Consequently, it cannot be neglected and needs to be accounted for. We demonstrate how one can do so using spatially resolved observations of galaxies in two filters. We show that HST observations in the F606W and F814W filters allow us to model and reduce the bias by an order of magnitude, sufficient to meet Euclid's scientific requirements. The precision of the correction is ultimately determined by the number of galaxies for which spatially-resolved observations in at least two filters are available. We use results from the Millennium Simulation to demonstrate that archival HST data will be sufficient for the tomographic cosmic shear analysis with the Euclid dataset.

The first half of this paper explores the origin of systematic biases in the measurement of weak gravitational lensing. Compared to previous work, we expand the investigation of PSF instability and fold in for the first time the effects of non-idealities in electronic imaging detectors and imperfect galaxy shape measurement algorithms. Together, these now explain the additive A(l) and multiplicative M(l) systematics typically reported in current lensing measurements. We find that overall performance is driven by a product of a telescope/camera's *absolute performance*, and our *knowledge about its performance*. The second half of this paper propagates any residual shear measurement biases through to their effect on cosmological parameter constraints. Fully exploiting the statistical power of Stage IV weak lensing surveys will require additive biases A<1.8e-12 and multiplicative biases M<4.0e-3. These can be allocated between individual budgets in hardware, calibration data and software, using results from the first half of the paper. If instrumentation is stable and well-calibrated, we find extant shear measurement software from GREAT10 already meet requirements on galaxies detected at S/N=40. Averaging over a population of galaxies with a realistic distribution of sizes, it also meets requirements for a 2D cosmic shear analysis from space. If used on fainter galaxies or for 3D cosmic shear tomography, existing algorithms would need calibration on simulations to avoid introducing bias at a level similar to the statistical error. Requirements on hardware and calibration data are discussed in more detail in a companion paper. Our analysis is intentionally general, but is specifically being used to drive the hardware and ground segment performance budget for the design of the European Space Agency's recently-selected Euclid mission.

In December 2010, NASA created a Science Definition Team (SDT) for WFIRST, the Wide Field Infra-Red Survey Telescope, recommended by the Astro 2010 Decadal Survey as the highest priority for a large space mission. The SDT was chartered to work with the WFIRST Project Office at GSFC and the Program Office at JPL to produce a Design Reference Mission (DRM) for WFIRST. Part of the original charge was to produce an interim design reference mission by mid-2011. That document was delivered to NASA and widely circulated within the astronomical community. In late 2011 the Astrophysics Division augmented its original charge, asking for two design reference missions. The first of these, DRM1, was to be a finalized version of the interim DRM, reducing overall mission costs where possible. The second of these, DRM2, was to identify and eliminate capabilities that overlapped with those of NASA's James Webb Space Telescope (henceforth JWST), ESA's Euclid mission, and the NSF's ground-based Large Synoptic Survey Telescope (henceforth LSST), and again to reduce overall mission cost, while staying faithful to NWNH. This report presents both DRM1 and DRM2.

We report on the detection of gravitational lensing magnification by a population of galaxy groups, at a significance level of 4.9 sigma. Using X-ray selected groups in the COSMOS 1.64 deg^2 field, and high-redshift Lyman break galaxies as sources, we measure a lensing-induced angular cross-correlation between the samples. After satisfying consistency checks that demonstrate we have indeed detected a magnification signal, and are not suffering from contamination by physical overlap of samples, we proceed to implement an optimally weighted cross-correlation function to further boost the signal to noise of the measurement. Interpreting this optimally weighted measurement allows us to study properties of the lensing groups. We model the full distribution of group masses using a composite-halo approach, considering both the singular isothermal sphere and Navarro-Frenk-White profiles, and find our best fit values to be consistent with those recovered using the weak-lensing shear technique. We argue that future weak-lensing studies will need to incorporate magnification along with shear, both to reduce residual systematics and to make full use of all available source information, in an effort to maximize scientific yield of the observations.

Using observations in the COSMOS field, we report an intriguing correlation between the star formation activity of massive (~10^{11.4}\msol) central galaxies, their stellar masses, and the large-scale (~10 Mpc) environments of their group-mass (~10^{13.6}\msol) dark matter halos. Probing the redshift range z=[0.2,1.0], our measurements come from two independent sources: an X-ray detected group catalog and constraints on the stellar-to-halo mass relation derived from a combination of clustering and weak lensing statistics. At z=1, we find that the stellar mass in star-forming centrals is a factor of two less than in passive centrals at the same halo mass. This implies that the presence or lack of star formation in group-scale centrals cannot be a stochastic process. By z=0, the offset reverses, probably as a result of the different growth rates of these objects. A similar but weaker trend is observed when dividing the sample by morphology rather than star formation. Remarkably, we find that star-forming centrals at z~1 live in groups that are significantly more clustered on 10 Mpc scales than similar mass groups hosting passive centrals. We discuss this signal in the context of halo assembly and recent simulations, suggesting that star-forming centrals prefer halos with higher angular momentum and/or formation histories with more recent growth; such halos are known to evolve in denser large-scale environments. If confirmed, this would be evidence of an early established link between the assembly history of halos on large scales and the future properties of the galaxies that form inside them.

We present a method for attaining sub-arcsecond pointing stability during sub- orbital balloon flights, as designed for in the High Altitude Lensing Observatory (HALO) concept. The pointing method presented here has the potential to perform near-space quality optical astronomical imaging at 1-2% of the cost of space-based missions. We also discuss an architecture that can achieve sufficient thermomechanical stability to match the pointing stability. This concept is motivated by advances in the development and testing of Ultra Long Duration Balloon (ULDB) flights which promise to allow observation campaigns lasting more than three months. The design incorporates a multi-stage pointing architecture comprising: a gondola coarse azimuth control system, a multi-axis nested gimbal frame structure with arcsecond stability, a telescope de-rotator to eliminate field rotation, and a fine guidance stage consisting of both a telescope mounted angular rate sensor and guide CCDs in the focal plane to drive a fast-steering mirror. We discuss the results of pointing tests together with a preliminary thermo-mechanical analysis required for sub-arcsecond pointing at high altitude. Possible future applications in the areas of wide-field surveys and exoplanet searches are also discussed.

The COSMOS field has been the subject of a wide range of observations, with a number of studies focusing on reconstructing the 3D dark matter density field. Typically, these studies have focused on one given method or tracer. In this paper, we reconstruct the distribution of mass in the COSMOS field out to a redshift z=1 by combining Hubble Space Telescope weak lensing measurements with zCOSMOS spectroscopic measurements of galaxy clustering. The distribution of galaxies traces the distribution of mass with high resolution (particularly in redshift, which is not possible with lensing), and the lensing data empirically calibrates the mass normalisation (bypassing the need for theoretical models). Two steps are needed to convert a galaxy survey into a density field. The first step is to create a smooth field from the galaxy positions, which is a point field. We investigate four possible methods for this: (i) Gaussian smoothing, (ii) convolution with truncated isothermal sphere, (iii) fifth nearest neighbour smoothing and (iv) a muliti-scale entropy method. The second step is to rescale this density field using a bias prescription. We calculate the optimal bias scaling for each method by comparing predictions from the smoothed density field with the measured weak lensing data, on a galaxy-by-galaxy basis. In general, we find scale-independent bias for all the smoothing schemes, to a precision of 10%. For the nearest neighbour smoothing case, we find the bias to be 2.51\pm 0.25. We also find evidence for a strongly evolving bias, increasing by a factor of ~3.5 between redshifts 0<z<0.8. We believe this strong evolution can be explained by the fact that we use a flux limited sample to build the density field.

This document describes the exposure time calculator for the Wide-Field Infrared Survey Telescope (WFIRST) high-latitude survey. The calculator works in both imaging and spectroscopic modes. In addition to the standard ETC functions (e.g. background and S/N determination), the calculator integrates over the galaxy population and forecasts the density and redshift distribution of galaxy shapes usable for weak lensing (in imaging mode) and the detected emission lines (in spectroscopic mode). The source code is made available for public use.

In this paper we present results from the Mapping Dark Matter competition that expressed the weak lensing shape measurement task in its simplest form and as a result attracted over 700 submissions in 2 months and a factor of 3 improvement in shape measurement accuracy on high signal to noise galaxies, over previously published results, and a factor 10 improvement over methods tested on constant shear blind simulations. We also review weak lensing shape measurement challenges, including the Shear TEsting Programmes (STEP1 and STEP2) and the GRavitational lEnsing Accuracy Testing competitions (GREAT08 and GREAT10).

In the theory of structure formation, galaxies are biased tracers of the underlying matter density field. The statistical relation between galaxy and matter density field is commonly referred as galaxy bias. In this paper, we test the linear bias model with weak-lensing and galaxy clustering measurements in the 2 square degrees COSMOS field (Scoville et al. 2007). We estimate the bias of galaxies between redshifts z=0.2 and z=1, and over correlation scales between R=0.2 h^-1 Mpc and R=15 h^-1 Mpc. We focus on three galaxy samples, selected in flux (simultaneous cuts I_814W < 26.5 and K_s < 24), and in stellar-mass (10^9 < M_* < 10^10 h^-2 Msun and 10^10 < M^*< 10^11 h^-2 Msun). At scales R > 2 h^-1 Mpc, our measurements support a model of bias increasing with redshift. The Tinker et al. (2010) fitting function provides a good fit to the data. We find the best fit mass of the galaxy halos to be log(M_200 h^-1 Msun) = 11.7^+0.6_-1.3 and log(M_200 h^-1 Msun) = 12.4^+0.2_-2.9 respectively for the low and high stellar-mass samples. In the halo model framework, bias is scale-dependent with a change of slope at the transition scale between the one and the two halo terms. We detect a scale-dependence of bias with a turn-down at scale R=2.3\pm1.5 h^-1 Mpc, in agreement with previous galaxy clustering studies. We find no significant amount of stochasticity, suggesting that a linear bias model is sufficient to describe our data. We use N-body simulations to quantify both the amount of cosmic variance and systematic errors in the measurement.

We present the first weak gravitational lensing analysis of the completed Canada-France-Hawaii Telescope Legacy Survey (CFHTLS). We study the 64 square degrees W1 field, the largest of the CFHTLS-Wide survey fields, and present the largest contiguous weak lensing convergence "mass map" yet made. 2.66 million galaxy shapes are measured, using a KSB pipeline verified against high-resolution Hubble Space Telescope imaging that covers part of the CFHTLS. Our i'-band measurements are also consistent with an analysis of independent r'-band imaging. The reconstructed lensing convergence map contains 301 peaks with signal-to-noise ratio {\nu}>3.5, consistent with predictions of a {\Lambda}CDM model. Of these peaks, 126 lie within 3.0' of a BCG identified from multicolor optical imaging in an independent, red sequence survey. We also identify 7 counterparts for massive clusters previously seen in X-ray emission within 6 square degrees XMM-LSS survey. With photometric redshift estimates for the source galaxies, we use a tomographic lensing method to fit the redshift and mass of each convergence peak. Matching these to the optical observations, we confirm 85 groups/clusters with \chi^2 reduced < 3.0, at a mean redshift <z_c> = 0.36 and velocity dispersion <\sigma_c> = 658.8 km/s. Future surveys, such as DES, LSST, KDUST and EUCLID, will be able to apply these same techniques to map clusters in much larger volumes and thus tightly constrain cosmological models.

GRavitational lEnsing Accuracy Testing 2010 (GREAT10) is a public image analysis challenge aimed at the development of algorithms to analyze astronomical images. Specifically, the challenge is to measure varying image distortions in the presence of a variable convolution kernel, pixelization and noise. This is the second in a series of challenges set to the astronomy, computer science and statistics communities, providing a structured environment in which methods can be improved and tested in preparation for planned astronomical surveys. GREAT10 extends upon previous work by introducing variable fields into the challenge. The "Galaxy Challenge" involves the precise measurement of galaxy shape distortions, quantified locally by two parameters called shear, in the presence of a known convolution kernel. Crucially, the convolution kernel and the simulated gravitational lensing shape distortion both now vary as a function of position within the images, as is the case for real data. In addition, we introduce the "Star Challenge" that concerns the reconstruction of a variable convolution kernel, similar to that in a typical astronomical observation. This document details the GREAT10 Challenge for potential participants. Continually updated information is also available from http://www.greatchallenges.info.

Weak lensing is commonly measured using shear through galaxy ellipticities, or using the effect of magnification bias on galaxy number densities. Here, we report on the first detection of weak lensing magnification with a new, independent technique using the distribution of galaxy sizes and magnitudes. These data come for free in galaxy surveys designed for measuring shear. We present the magnification estimator and apply it to an X-ray selected sample of galaxy groups in the COSMOS HST survey. The measurement of the projected surface density \Sigma(r) is consistent with the shear measurements within the uncertainties, and has roughly 40% of the signal-to-noise of the latter. We discuss systematic issues and challenges to realizing the potential of this new probe of weak lensing.

Gravitational lensing can provide pure geometric tests of the structure of space-time, for instance by determining empirically the angular diameter distance-redshift relation. This geometric test has been demonstrated several times using massive clusters which produce a large lensing signal. In this case, matter at a single redshift dominates the lensing signal, so the analysis is straightforward. It is less clear how weaker signals from multiple sources at different redshifts can be stacked to demonstrate the geometric dependence. We introduce a simple measure of relative shear which for flat cosmologies separates the effect of lens and source positions into multiplicative terms, allowing signals from many different source-lens pairs to be combined. Applying this technique to a sample of groups and low-mass clusters in the COSMOS survey, we detect a clear variation of shear with distance behind the lens. This represents the first detection of the geometric effect using weak lensing by multiple, low-mass systems. The variation of distance with redshift is measured with sufficient precision to constrain the equation of state of the universe under the assumption of flatness, equivalent to a detection of a dark energy component Omega_X at greater than 99% confidence for an equation-of-state parameter -2.5 < w < -0.1. For the case w = -1, we find a value for the cosmological constant density parameter Omega_Lambda = 0.85+0.044-0.19 (68% C.L.), and detect cosmic acceleration (q_0 < 0) at the 98% C.L.. We consider the systematic uncertainties associated with this technique and discuss the prospects for applying it in forthcoming weak-lensing surveys.

We explore the sensitivity of weak gravitational lensing to second-order corrections to the spacetime metric within a cosmological adaptation of the parameterized post-Newtonian framework. Whereas one might expect nonlinearities of the gravitational field to introduce non-Gaussianity into the statistics of the lensing convergence field, we show that such corrections are actually always small within a broad class of scalar-tensor theories of gravity. We show this by first computing the weak lensing convergence within our parameterized framework to second order in the gravitational potential, and then computing the relevant post-Newtonian parameters for scalar-tensor gravity theories. In doing so we show that this potential systematic factor is generically negligible, thus clearing the way for weak lensing to provide a direct tracer of mass on cosmological scales for a wide class of gravity theories despite uncertainties in the precise nature of the departures from general relativity.

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

We present a detailed strong lensing, weak lensing and X-ray analysis of Abell 2744 (z = 0.308), one of the most actively merging galaxy clusters known. It appears to have unleashed `dark', `ghost', `bullet' and `stripped' substructures, each ~10^14 solar masses. The phenomenology is complex and will present a challenge for numerical simulations to reproduce. With new, multiband HST imaging, we identify 34 strongly-lensed images of 11 galaxies around the massive Southern `core'. Combining this with weak lensing data from HST, VLT and Subaru, we produce the most detailed mass map of this cluster to date. We also perform an independent analysis of archival Chandra X-ray imaging. Our analyses support a recent claim that the Southern core and Northwestern substructure are post-merger and exhibit morphology similar to the Bullet Cluster viewed from an angle. From the separation between X-ray emitting gas and lensing mass in the Southern core, we derive a new and independent constraint on the self-interaction cross section of dark matter particles sigma/m <~ 3 \pm 1 cm^2 g^-1. In the Northwestern substructure, the gas, dark matter, and galaxy components have become separated by much larger distances. Most curiously, the `ghost' clump (primarily gas) leads the `dark' clump (primarily dark matter) by more than 150 kpc. We propose an enhanced `ram-pressure slingshot' scenario which may have yielded this reversal of components with such a large separation, but needs further confirmation by follow-up observations and numerical simulations. A secondary merger involves a second `bullet' clump in the North and an extremely `stripped' clump to the West. The latter appears to exhibit the largest separation between dark matter and X-ray emitting baryons detected to date in our sky.

Point Spread Function (PSF) modeling is a central part of any astronomy data analysis relying on measuring the shapes of objects. It is especially crucial for weak gravitational lensing, in order to beat down systematics and allow one to reach the full potential of weak lensing in measuring dark energy. A PSF modeling pipeline is made of two main steps: the first one is to assess its shape on stars, and the second is to interpolate it at any desired position (usually galaxies). We focus on the second part, and compare different interpolation schemes, including polynomial interpolation, radial basis functions, Delaunay triangulation and Kriging. For that purpose, we develop simulations of PSF fields, in which stars are built from a set of basis functions defined from a Principal Components Analysis of a real ground-based image. We find that Kriging gives the most reliable interpolation, significantly better than the traditionally used polynomial interpolation. We also note that although a Kriging interpolation on individual images is enough to control systematics at the level necessary for current weak lensing surveys, more elaborate techniques will have to be developed to reach future ambitious surveys' requirements.

Current and upcoming wide-field, ground-based, broad-band imaging surveys promise to address a wide range of outstanding problems in galaxy formation and cosmology. Several such uses of ground-based data, especially weak gravitational lensing, require highly precise measurements of galaxy image statistics with careful correction for the effects of the point-spread function (PSF). In this paper, we introduce the SHERA (SHEar Reconvolution Analysis) software to simulate ground-based imaging data with realistic galaxy morphologies and observing conditions, starting from space-based data (from COSMOS, the Cosmological Evolution Survey) and accounting for the effects of the space-based PSF. This code simulates ground-based data, optionally with a weak lensing shear applied, in a model-independent way using a general Fourier space formalism. The utility of this pipeline is that it allows for a precise, realistic assessment of systematic errors due to the method of data processing, for example in extracting weak lensing galaxy shape measurements or galaxy radial profiles, given user-supplied observational conditions and real galaxy morphologies. Moreover, the simulations allow for the empirical test of error estimates and determination of parameter degeneracies, via generation of many noise maps. The public release of this software, along with a large sample of cleaned COSMOS galaxy images (corrected for charge transfer inefficiency), should enable upcoming ground-based imaging surveys to achieve their potential in the areas of precision weak lensing analysis, galaxy profile measurement, and other applications involving detailed image analysis.

A simple, yet general, formalism for the optimized linear combination of astrophysical images is constructed and demonstrated. The formalism allows the user to combine multiple undersampled images to provide oversampled output at high precision. The proposed method is general and may be used for any configuration of input pixels and point spread function; it also provides the noise covariance in the output image along with a powerful metric for describing undesired distortion to the image convolution kernel. The method explicitly provides knowledge and control of the inevitable compromise between noise and fidelity in the output image. We present a first prototype implementation of the method, outlining steps taken to generate an efficient algorithm. This implementation is then put to practical use in reconstructing fully-sampled output images using simulated, undersampled input exposures that are designed to mimic the proposed \emph{Wide Field InfraRed Survey Telescope} (\emph{WFIRST}). We examine results using randomly rotated and dithered input images, while also assessing better-known "ideal" dither patterns: comparing results we illustrate the use of the method as a survey design tool. Finally, we use the method to test the robustness of linear image combination when subject to practical realities such as missing input pixels and focal plane plate scale variations.

In December 2010, NASA created a Science Definition Team (SDT) for WFIRST, the Wide Field Infra-Red Survey Telescope, recommended by the Astro 2010 Decadal Survey as the highest priority for a large space mission. The SDT was chartered to work with the WFIRST Project Office at GSFC and the Program Office at JPL to produce a Design Reference Mission (DRM) for WFIRST. This paper describes an Interim DRM. The DRM will be completed in 2012.

Future orbiting observatories will survey large areas of sky in order to constrain the physics of dark matter and dark energy using weak gravitational lensing and other methods. Lossy compression of the resultant data will improve the cost and feasibility of transmitting the images through the space communication network. We evaluate the consequences of the lossy compression algorithm of Bernstein et al. (2010) for the high-precision measurement of weak-lensing galaxy ellipticities. This square-root algorithm compresses each pixel independently, and the information discarded is by construction less than the Poisson error from photon shot noise. For simulated space-based images (without cosmic rays) digitized to the typical 16 bits per pixel, application of the lossy compression followed by image-wise lossless compression yields images with only 2.4 bits per pixel, a factor of 6.7 compression. We demonstrate that this compression introduces no bias in the sky background. The compression introduces a small amount of additional digitization noise to the images, and we demonstrate a corresponding small increase in ellipticity measurement noise. The ellipticity measurement method is biased by the addition of noise, so the additional digitization noise is expected to induce a multiplicative bias on the galaxies' measured ellipticities. After correcting for this known noise-induced bias, we find a residual multiplicative ellipticity bias of m ~ -4x10^{-4}. This bias is small when compared to the many other issues that precision weak lensing surveys must confront, and furthermore we expect it to be reduced further with better calibration of ellipticity measurement methods.

Cosmic shear has been identified as the method with the most potential to constrain dark energy. To capitalise on this potential it is necessary to measure galaxy shapes with great accuracy, which in turn requires a detailed model for the image blurring, the Point Spread Function (PSF). In general the PSF varies with wavelength and therefore the PSF integrated over an observing filter depends on the spectrum of the object. For a typical galaxy the spectrum varies across the galaxy image, thus the PSF depends on the position within the image. We estimate the bias on the shear due to such colour gradients by modelling galaxies using two co-centered, co-elliptical Sersic profiles, each with a different spectrum. We estimate the effect of ignoring colour gradients and find the shear bias from a single galaxy can be very large depending on the properties of the galaxy. We find that halving the filter width reduces the shear bias by a factor of about 5. We show that, to first order, tomographic cosmic shear two point statistics depend on the mean shear bias over the galaxy population at a given redshift. For a single broad filter, and averaging over a small galaxy catalogue, we find a mean shear bias which is subdominant to the predicted statistical errors for future cosmic shear surveys. However, the true mean shear bias may exceed the statistical errors, depending on how accurately the catalogue represents the observed distribution of galaxies in the cosmic shear survey. We then investigate the bias on the shear for two-filter imaging and find that the bias is reduced by at least an order of magnitude. Lastly, we find that it is possible to calibrate galaxies for which colour gradients were ignored using two-filter imaging of a fair sample of noisy galaxies, if the galaxy model is known. For a S/N of 25 the number of galaxies required in each tomographic redshift bin is of order 1e4.

We present ACS, NICMOS, and Keck AO-assisted photometry of 20 Type Ia supernovae SNe Ia from the HST Cluster Supernova Survey. The SNe Ia were discovered over the redshift interval 0.623 < z < 1.415. Fourteen of these SNe Ia pass our strict selection cuts and are used in combination with the world's sample of SNe Ia to derive the best current constraints on dark energy. Ten of our new SNe Ia are beyond redshift $z=1$, thereby nearly doubling the statistical weight of HST-discovered SNe Ia beyond this redshift. Our detailed analysis corrects for the recently identified correlation between SN Ia luminosity and host galaxy mass and corrects the NICMOS zeropoint at the count rates appropriate for very distant SNe Ia. Adding these supernovae improves the best combined constraint on the dark energy density \rho_{DE}(z) at redshifts 1.0 < z < 1.6 by 18% (including systematic errors). For a LambdaCDM universe, we find \Omega_\Lambda = 0.724 +0.015/-0.016 (68% CL including systematic errors). For a flat wCDM model, we measure a constant dark energy equation-of-state parameter w = -0.985 +0.071/-0.077 (68% CL). Curvature is constrained to ~0.7% in the owCDM model and to ~2% in a model in which dark energy is allowed to vary with parameters w_0 and w_a. Tightening further the constraints on the time evolution of dark energy will require several improvements, including high-quality multi-passband photometry of a sample of several dozen z>1 SNe Ia. We describe how such a sample could be efficiently obtained by targeting cluster fields with WFC3 on HST.

Using data from the COSMOS survey, we perform the first joint analysis of galaxy-galaxy weak lensing, galaxy spatial clustering, and galaxy number densities. Carefully accounting for sample variance and for scatter between stellar and halo mass, we model all three observables simultaneously using a novel and self-consistent theoretical framework. Our results provide strong constraints on the shape and redshift evolution of the stellar-to-halo mass relation (SHMR) from z=0.2 to z=1. At low stellar mass, we find that halo mass scales as Mh M*^0.46 and that this scaling does not evolve significantly with redshift to z=1. We show that the dark-to-stellar ratio, Mh/M*, varies from low to high masses, reaching a minimum of Mh/M*~27 at M*=4.5x10^10 Msun and Mh=1.2x10^12 Msun. This minimum is important for models of galaxy formation because it marks the mass at which the accumulated stellar growth of the central galaxy has been the most efficient. We describe the SHMR at this minimum in terms of the "pivot stellar mass", M*piv, the "pivot halo mass", Mhpiv, and the "pivot ratio", (Mh/M*)piv. Thanks to a homogeneous analysis of a single data set, we report the first detection of mass downsizing trends for both Mhpiv and M*piv. The pivot stellar mass decreases from M*piv=5.75+-0.13x10^10 Msun at z=0.88 to M*piv=3.55+-0.17x10^10 Msun at z=0.37. Intriguingly, however, the corresponding evolution of Mhpiv leaves the pivot ratio constant with redshift at (Mh/M*)piv~27. We use simple arguments to show how this result raises the possibility that star formation quenching may ultimately depend on Mh/M* and not simply Mh, as is commonly assumed. We show that simple models with such a dependence naturally lead to downsizing in the sites of star formation. Finally, we discuss the implications of our results in the context of popular quenching models, including disk instabilities and AGN feedback.

We present the 'simage' software suite for the simulation of artificial extragalactic images, based empirically around real observations of the Hubble Ultra Deep Field (UDF). The simulations reproduce galaxies with realistic and complex morphologies via the modeling of UDF galaxies as shapelets. Images can be created in the B, V, i and z bands for both space- and ground-based telescopes and instruments. The simulated images can be produced for any required field size, exposure time, Point Spread Function (PSF), telescope mirror size, pixel resolution, field star density, and a variety of detector noise sources. It has the capability to create images with both a pre-determined number of galaxies or one calibrated to the number counts of pre-existing data sets such as the HST COSMOS survey. In addition, simple options are included to add a known weak gravitational lensing (both shear and flexion) to the simulated images. The software is available in Interactive Data Language (IDL) and can be freely downloaded for scientific, developmental and teaching purposes.

[NIRSS is one of three concepts that contributed to the Wide-Field Infrared Survey Telescope (WFIRST) mission advocated by the Decadal Survey.] Operating beyond the reaches of the Earth's atmosphere, free of its limiting absorption and thermal background, the Near-Infrared Sky Surveyor (NIRSS) will deeply map the entire sky at near-infrared wavelengths, thereby enabling new and fundamental discoveries ranging from the identification of extrasolar planets to probing the reionization epoch by identifying thousands of quasars at z>10. NIRSS will directly address the NASA scientific objective of studying cosmic origins by using a 1.5-meter telescope to reach full-sky 0.2 uJy (25.6 mag AB) sensitivities in four passbands from 1 to 4 microns in a 4-yr mission. At the three shorter passbands (1 - 2.5 microns), the proposed depth is comparable to the deepest pencil-beam surveys done to date and is 3000 times more sensitive than the only previous all-sky near-infrared survey, 2MASS. At the longest passband (3.5 micron), which is not feasible from the ground, NIRSS will be 500 times more sensitive than WISE. NIRSS fills a pivotal gap in our knowledge of the celestial sphere, is a natural complement to WISE, and is well matched to the next generation of deep (0.1 uJy), wide-area (>2 pi ster), ground-based optical surveys (LSST and Pan-Starrs). With the high thermal backgrounds of ground-based infrared observations, a near-infrared full sky survey at sub-uJy sensitivity is only feasible from space.

We investigate the (large-scale) bar fraction in a mass-complete sample of M > 10^10.5 Msun disk galaxies at 0.2 < z < 0.6 in the COSMOS field. The fraction of barred disks strongly depends on mass, disk morphology, and specific star formation rate (SSFR). At intermediate stellar mass (10^10.5 < M < 10^11 Msun) the bar fraction in early-type disks is much higher, at all redshifts, by a factor ~2, than that in late-type disks. This trend is reversed at higher stellar mass (M > 10^11 Msun), where the fraction of bars in early-type disks becomes significantly lower, at all redshifts, than that in late-type disks. The bar fractions for galaxies with low and high SSFRs closely follow those of the morphologically-selected early-type and late-type populations, respectively. This indicates a close correspondence between morphology and SSFR in disk galaxies at these earlier epochs. Interestingly, the total bar fraction in 10^10.5 < M < 10^11 Msun disks is built up by a factor of ~2 over the redshift interval explored, while for M > 10^11 Msun disks it remains roughly constant. This indicates that, already by z ~ 0.6, spectral and morphological transformations in the most massive disk galaxies have largely converged to the familiar Hubble sequence that we observe in the local Universe, while for intermediate mass disks this convergence is ongoing until at least z ~ 0.2. Moreover, these results highlight the importance of employing mass-limited samples for quantifying the evolution of barred galaxies. Finally, the evolution of the barred galaxy populations investigated does not depend on the large-scale environmental density (at least, on the scales which can be probed with the available photometric redshifts).

We compare the surface brightness-inclination relation for a sample of COSMOS pure disk galaxies at z~0.7 with an artificially redshifted sample of SDSS disks well matched to the COSMOS sample in terms of rest-frame photometry and morphology, as well as their selection and analysis. The offset between the average surface brightness of face-on and edge-on disks in the redshifted SDSS sample matches that predicted by measurements of the optical depth of galactic disks in the nearby universe. In contrast, large disks at z~0.7 have a virtually flat surface brightness-inclination relation, suggesting that they are more opaque than their local counterparts. This could be explained by either an increased amount of optically thick material in disks at higher redshift, or a different spatial distribution of the dust.

(Abridged) We examine the effects of charge transfer inefficiency (CTI) during CCD readout on galaxy shape measurements required by studies of weak gravitational lensing. We simulate a CCD readout with CTI such as that caused by charged particle radiation damage. We verify our simulations on data from laboratory-irradiated CCDs. Only charge traps with time constants of the same order as the time between row transfers during readout affect galaxy shape measurements. We characterize the effects of CTI on various galaxy populations. We baseline our study around p-channel CCDs that have been shown to have charge transfer efficiency up to an order of magnitude better than several models of n-channel CCDs designed for space applications. We predict that for galaxies furthest from the readout registers, bias in the measurement of galaxy shapes, Delta(e), will increase at a rate of 2.65 +/- 0.02 x 10^(-4) per year at L2 for accumulated radiation exposure averaged over the solar cycle. If uncorrected, this will consume the entire shape measurement error budget of a dark energy mission within about 4 years. Software mitigation techniques demonstrated elsewhere can reduce this by a factor of ~10, bringing the effect well below mission requirements. CCDs with higher CTI than the ones we studeied may not meet the requirements of future dark energy missions. We discuss ways in which hardware could be designed to further minimize the impact of CTI.

Cosmic shear requires high precision measurement of galaxy shapes in the presence of the observational Point Spread Function (PSF) that smears out the image. The PSF must therefore be known for each galaxy to a high accuracy. However, for several reasons, the PSF is usually wavelength dependent, therefore the differences between the spectral energy distribution of the observed objects introduces further complexity. In this paper we investigate the effect of the wavelength-dependence of the PSF, focusing on instruments in which the PSF size is dominated by the diffraction-limit of the telescope and which use broad-band filters for shape measurement. We first calculate biases on cosmological parameter estimation from cosmic shear when the stellar PSF is used uncorrected. Using realistic galaxy and star spectral energy distributions and populations and a simple three-component circular PSF we find that the colour-dependence must be taken into account for the next generation of telescopes. We then consider two different methods for removing the effect (i) the use of stars of the same colour as the galaxies and (ii) estimation of the galaxy spectral energy distribution using multiple colours and using a telescope model for the PSF. We find that both of these methods correct the effect to levels below the tolerances required for per-cent level measurements of dark energy parameters. Comparison of the two methods favours the template-fitting method because its efficiency is less dependent on galaxy redshift than the broad-band colour method and takes full advantage of deeper photometry.

We investigate data compression schemes for proposed all-sky diffraction-limited visible/NIR sky surveys aimed at the dark energy problem. We show that lossy square-root compression to 1 bit of noise per pixel, followed by standard lossless compression algorithms, reduces the images to 2.5-4 bits per pixel, depending primarily upon the level of cosmic-ray contamination of the images. Compression to this level adds noise equivalent to <10% penalty in observing time. We derive an analytic correction to flux biases inherent to the square-root compression scheme. Numerical tests on simple galaxy models confirm that galaxy fluxes and shapes are measured with systematic biases <~10^{-4} induced by the compression scheme, well below the requirements of supernova and weak gravitational lensing dark-energy experiments. An accompanying paper (Vanderveld 2009) bounds the shape biases using realistic simulated images of the high-Galactic-latitude sky. The square-root preprocessing step has advantages over simple (linear) decimation when there are many bright objects or cosmic rays in the field, or when the background level will vary.

Measurements of X-ray scaling laws are critical for improving cosmological constraints derived with the halo mass function and for understanding the physical processes that govern the heating and cooling of the intracluster medium. In this paper, we use a sample of 206 X-ray selected galaxy groups to investigate the scaling relation between X-ray luminosity (Lx) and halo mass (M00) where M200 is derived via stacked weak gravitational lensing. This work draws upon a broad array of multi-wavelength COSMOS observations including 1.64 square degrees of contiguous imaging with the Advanced Camera for Surveys (ACS) and deep XMM-Newton/Chandra imaging. The combined depth of these two data-sets allows us to probe the lensing signals of X-ray detected structures at both higher redshifts and lower masses than previously explored. Weak lensing profiles and halo masses are derived for nine sub-samples, narrowly binned in luminosity and redshift. The COSMOS data alone are well fit by a power law, M200 ~ Lx^a, with a slope of a=0.66+-0.14. These results significantly extend the dynamic range for which the halo masses of X-ray selected structures have been measured with weak gravitational lensing. As a result, tight constraints are obtained for the slope of the M-Lx relation. The combination of our group data with previously published cluster data demonstrates that the M-Lx relation is well described by a single power law, a=0.64+-0.03, over two decades in mass, 10^13.5-10^15.5 h72^-1 Msun. These results are inconsistent at the 3.7 level with the self-similar prediction of a=0.75. We examine the redshift dependence of the M-Lx relation and find little evidence for evolution beyond the rate predicted by self-similarity from z ~ 0.25 to z ~ 0.8.

We use ~8,600 >5e10 Msol COSMOS galaxies to study how the morphological mix of massive ellipticals, bulge-dominated disks, intermediate-bulge disks, bulge-less disks and irregular galaxies evolves from z=0.2 to z=1. The morphological evolution depends strongly on mass. At M>3e11 Msol, no evolution is detected in the morphological mix: ellipticals dominate since z=1, and the Hubble sequence has quantitatively settled down by this epoch. At the 1e11 Msol mass scale, little evolution is detected, which can be entirely explained with major mergers. Most of the morphological evolution from z=1 to z=0.2 takes place at masses 5e10 - 1e11 Msol, where: (i) The fraction of spirals substantially drops and the contribution of early-types increases. This increase is mostly produced by the growth of bulge-dominated disks, which vary their contribution from ~10% at z=1 to >30% at z=0.2 (cf. the elliptical fraction grows from ~15% to ~20%). Thus, at these masses, transformations from late- to early-types result in disk-less elliptical morphologies with a statistical frequency of only 30% - 40%. Otherwise, the processes which are responsible for the transformations either retain or produce a non-negligible disk component. (ii) The bulge-less disk galaxies, which contribute ~15% to the intermediate-mass galaxy population at z=1, virtually disappear by z=0.2. The merger rate since z=1 is too low to account for the disappearance of these massive bulge-less disks, which most likely grow a bulge via secular evolution.

Charge Transfer Inefficiency (CTI) due to radiation damage above the Earth's atmosphere creates spurious trailing in Hubble Space Telescope (HST) images. Radiation damage also creates unrelated warm pixels - but these happen to be perfect for measuring CTI. We model CTI in the Advanced Camera for Surveys (ACS)/Wide Field Channel (WFC) and construct a physically motivated correction scheme. This operates on raw data, rather than secondary science products, by returning individual electrons to pixels from which they were unintentionally dragged during readout. We apply our correction to images from the HST COSMOS survey, successfully reducing the CTI trails by a factor of ~30 everywhere in the CCD and at all flux levels. We quantify changes in galaxy photometry, astrometry and shape. The remarkable 97% level of correction is more than sufficient to enable a (forthcoming) reanalysis of downstream science products, and the collection of larger surveys.

We present a new survey strategy to discover and study high redshift Type Ia supernovae (SNe Ia) using the Hubble Space Telescope (HST). By targeting massive galaxy clusters at 0.9<z<1.5, we obtain a twofold improvement in the efficiency of finding SNe compared to an HST field survey and a factor of three improvement in the total yield of SN detections in relatively dust-free red-sequence galaxies. In total, sixteen SNe were discovered at z>0.95, nine of which were in galaxy clusters. This strategy provides a SN sample that can be used to decouple the effects of host galaxy extinction and intrinsic color in high redshift SNe, thereby reducing one of the largest systematic uncertainties in SN cosmology.

We present the results of the GREAT08 Challenge, a blind analysis challenge to infer weak gravitational lensing shear distortions from images. The primary goal was to stimulate new ideas by presenting the problem to researchers outside the shear measurement community. Six GREAT08 Team methods were presented at the launch of the Challenge and five additional groups submitted results during the 6 month competition. Participants analyzed 30 million simulated galaxies with a range in signal to noise ratio, point-spread function ellipticity, galaxy size, and galaxy type. The large quantity of simulations allowed shear measurement methods to be assessed at a level of accuracy suitable for currently planned future cosmic shear observations for the first time. Different methods perform well in different parts of simulation parameter space and come close to the target level of accuracy in several of these. A number of fresh ideas have emerged as a result of the Challenge including a re-examination of the process of combining information from different galaxies, which reduces the dependence on realistic galaxy modelling. The image simulations will become increasingly sophisticated in future GREAT challenges, meanwhile the GREAT08 simulations remain as a benchmark for additional developments in shear measurement algorithms.

We investigate if the discrepancy between estimates of the total baryon mass fraction obtained from observations of the cosmic microwave background (CMB) and of galaxy groups/clusters persists when a large sample of groups is considered. To this purpose, 91 candidate X-ray groups/poor clusters at redshift 0.1 < z < 1 are selected from the COSMOS 2 deg^2 survey, based only on their X-ray luminosity and extent. This sample is complemented by 27 nearby clusters with a robust, analogous determination of the total and stellar mass inside R_500. The total sample of 118 groups and clusters with z < 1 spans a range in M_500 of ~10^13--10^15 M_sun. We find that the stellar mass fraction associated with galaxies at R_500 decreases with increasing total mass as (M_500)^-0.37 \pm 0.04, independent of redshift. Estimating the total gas mass fraction from a recently derived, high quality scaling relation, the total baryon mass fraction (f_500^stars+gas=f_500^stars+f_500^gas) is found to increase by ~ 25% when M_500 increases from <M>=5 X 10^13 M_sun to <M> = 7 X 10^14 M_sun. After consideration of a plausible contribution due to intra--cluster light (11--22% of the total stellar mass), and gas depletion through the hierarchical assembly process (10% of the gas mass), the estimated values of the total baryon mass fraction are still lower than the latest CMB measure of the same quantity (WMAP5), at a significance level of 3.3\sigma for groups of <M>=5 X 10^13~M_sun. The discrepancy decreases towards higher total masses, such that it is 1\sigma at <M>= 7 X 10^14~M_sun. We discuss this result in terms of non--gravitational processes such as feedback and filamentary heating.

The GRavitational lEnsing Accuracy Testing 2008 (GREAT08) Challenge focuses on a problem that is of crucial importance for future observations in cosmology. The shapes of distant galaxies can be used to determine the properties of dark energy and the nature of gravity, because light from those galaxies is bent by gravity from the intervening dark matter. The observed galaxy images appear distorted, although only slightly, and their shapes must be precisely disentangled from the effects of pixelisation, convolution and noise. The worldwide gravitational lensing community has made significant progress in techniques to measure these distortions via the Shear TEsting Program (STEP). Via STEP, we have run challenges within our own community, and come to recognise that this particular image analysis problem is ideally matched to experts in statistical inference, inverse problems and computational learning. Thus, in order to continue the progress seen in recent years, we are seeking an infusion of new ideas from these communities. This document details the GREAT08 Challenge for potential participants. Please visit http://www.great08challenge.info for the latest information.

We use HST/ACS images and a photometric catalog of the COSMOS field to analyze morphologies of the host galaxies of approximately 400 AGN candidates at redshifts 0.3 < z < 1.0. We compare the AGN hosts with a sample of non-active galaxies drawn from the COSMOS field to match the magnitude and redshift distribution of the AGN hosts. We perform 2-D surface brightness modeling with GALFIT to yield host galaxy and nuclear point source magnitudes. X-ray selected AGN host galaxy morphologies span a substantial range that peaks between those of early-type, bulge-dominated and late-type, disk-dominated systems. We also measure the asymmetry and concentration of the host galaxies. Unaccounted for, the nuclear point source can significantly bias results of these measured structural parameters, so we subtract the best-fit point source component to obtain images of the underlying host galaxies. Our concentration measurements reinforce the findings of our 2-D morphology fits, placing X-ray AGN hosts between early- and late-type inactive galaxies. AGN host asymmetry distributions are consistent with those of control galaxies. Combined with a lack of excess companion galaxies around AGN, the asymmetry distributions indicate that strong interactions are no more prevalent among AGN than normal galaxies. In light of recent work, these results suggest that the host galaxies of AGN at these X-ray luminosities may be in a transition from disk-dominated to bulge-dominated, but that this transition is not typically triggered by major mergers.

We present a method to simulate color, 3-dimensional images taken with a space-based observatory by building off of the established shapelets pipeline. The simulated galaxies exhibit complex morphologies, which are realistically correlated between, and include, known redshifts. The simulations are created using galaxies from the 4 optical and near-infrared bands (B, V, i and z) of the Hubble Ultra Deep Field (UDF) as a basis set to model morphologies and redshift. We include observational effects such as sky noise and pixelization and can add astronomical signals of interest such as weak gravitational lensing. The realism of the simulations is demonstrated by comparing their morphologies to the original UDF galaxies and by comparing their distribution of ellipticities as a function of redshift and magnitude to wider HST COSMOS data. These simulations have already been useful for calibrating multicolor image analysis techniques and for better optimizing the design of proposed space telescopes.

We present the first catalog of 67 strong galaxy-galaxy lens candidates discovered in the 1.64 square degree Hubble Space Telescope COSMOS survey. Twenty of these systems display multiple images or strongly curved large arcs. Our initial search is performed by visual inspection of the data and is restricted, for practical considerations, to massive early-type lens galaxies with arcs found at radii smaller than ~5''. Simple mass models are constructed for the best lens candidates and our results are compared to the strong lensing catalogs of the SLACS survey and the CASTLES database. These new strong galaxy-galaxy lensing systems constitute a valuable sample to study the mass distribution of early-type galaxies and their associated dark matter halos. We further expect this sample to play an important role in the testing of software algorithms designed to automatically search for strong gravitational lenses. From our analysis a robust lower limit is derived for the expected occurrence of strong galaxy-galaxy systems in current and future space-based wide-field imaging surveys. We expect that such surveys should uncover a large number of strong lensing systems (more than 10 systems per square degree), which will allow for a detailed statistical analysis of galaxy properties and will likely lead to constraints on models of gravitational structure formation and cosmology. The sample of strong lenses is available here: http://cosmosstronglensing.uni-hd.de/

We study in detail the photometric redshift requirements needed for tomographic weak gravitational lensing in order to measure accurately the Dark Energy equation of state. In particular, we examine how ground-based photometry (u,g,r,i,z,y) can be complemented by space-based near-infrared (IR) photometry (J,H), e.g. on board the planned DUNE satellite. Using realistic photometric redshift simulations and an artificial neural network photo-z method we evaluate the Figure of Merit for the Dark Energy parameters $(w_0, w_a)$. We consider a DUNE-like broad optical filter supplemented with ground-based multi-band optical data from surveys like the Dark Energy Survey, Pan-STARRS and LSST. We show that the Dark Energy Figure of Merit would improved by a factor of 1.3 to 1.7 if IR filters are added on board DUNE. Furthermore we show that with IR data catastrophic photo-z outliers can be removed effectively. There is an interplay between the choice of filters, the magnitude limits and the removal of outliers. We draw attention to the dependence of the results on the galaxy formation scenarios encoded into the mock galaxies, e.g the galaxy reddening. For example, deep u band data could be as effective as the IR. We also find that about $10^5-10^6$ spectroscopic redshifts are needed for calibration of the full survey.

The aim of this paper is to investigate ways to optimize the accuracy of photometric redshifts for a SNAP like mission. We focus on how the accuracy of the photometric redshifts depends on the magnitude limit and signal-to-noise ratio, wave-length coverage, number of filters and their shapes and observed galaxy type. We use simulated galaxy catalogs constructed to reproduce observed galaxy luminosity functions from GOODS, and derive photometric redshifts using a template fitting method. By using a catalog that resembles real data, we can estimate the expected number density of galaxies for which photometric redshifts can be derived. We find that the accuracy of photometric redshifts is strongly dependent on the signal-to-noise (S/N) (i.e., S/N>10 is needed for accurate photometric redshifts). The accuracy of the photometric redshifts is also dependent on galaxy type, with smaller scatter for earlier type galaxies. Comparing results using different filter sets, we find that including the U-band is important for decreasing the fraction of outliers, i.e., ``catastrophic failures''. Using broad overlapping filters with resolution ~4gives better photometric redshifts compared to narrower filters (resolution >~5) with the same integration time. We find that filters with square response curves result in a slightly higher scatter, mainly due to a higher fraction of outliers at faint magnitudes. We also compare a 9-filter set to a 17-filter set, where we assume that the available exposure time per filter in the latter set is half that of the first set. We find that the 9-filter set gives more accurate redshifts for a larger number of objects and reaches higher redshift, while the 17-filter set is gives better results at bright magnitudes.

We assess the relative merits of weak lensing surveys, using overlapping imaging data from the ground-based Subaru telescope and the Hubble Space Telescope (HST). Our tests complement similar studies undertaken with simulated data. From observations of 230,000 matched objects in the 2 square degree COSMOS field, we identify the limit at which faint galaxy shapes can be reliably measured from the ground. Our ground-based shear catalog achieves sub-percent calibration bias compared to high resolution space-based data, for galaxies brighter than i'~24.5 and with half-light radii larger than 1.8". This selection corresponds to a surface density of ~15 galaxies per sq arcmin compared to ~71 per sq arcmin from space. On the other hand the survey speed of current ground-based facilities is much faster than that of HST, although this gain is mitigated by the increased depth of space-based imaging desirable for tomographic (3D) analyses. As an independent experiment, we also reconstruct the projected mass distribution in the COSMOS field using both data sets, and compare the derived cluster catalogs with those from X-ray observations. The ground-based catalog achieves a reasonable degree of completeness, with minimal contamination and no detected bias, for massive clusters at redshifts 0.2<z<0.5. The space-based data provide improved precision and a greater sensitivity to clusters of lower mass or at higher redshift.

To derive a new H$\alpha$ luminosity function and to understand the clustering properties of star-forming galaxies at $z \approx 0.24$, we have made a narrow-band imaging survey for H$\alpha$ emitting galaxies in the HST COSMOS 2 square degree field. We used the narrow-band filter NB816 ($\lambda_c = 8150$ \AA, $\Delta \lambda = 120$ \AA) and sampled H$\alpha$ emitters with $EW_{\rm obs}(\rm H\alpha + [N\textsc{ii}]) > 12$ \AA in a redshift range between $z=0.233$ and $z=0.251$ corresponding to a depth of 70 Mpc. We obtained 980 H$\alpha$ emitting galaxies in a sky area of 5540 arcmin$^2$, corresponding to a survey volume of $3.1 \times 10^4 {\rm Mpc^3}$. We derive a H$\alpha$ luminosity function with a best-fit Schechter function parameter set of $\alpha = -1.35^{+0.11}_{-0.13}$, $\log\phi_* = -2.65^{+0.27}_{-0.38}$, and $\log L_* ({\rm erg s^{-1}}) = 41.94^{+0.38}_{-0.23}$. The H$\alpha$ luminosity density is $2.7^{+0.7}_{-0.6} \times 10^{39}$ ergs s$^{-1}$ Mpc$^{-3}$. After subtracting the AGN contribution (15 %) to the H$\alpha$ luminosity density, the star formation rate density is evaluated as $1.8^{+0.7}_{-0.4} \times 10^{-2}$ $M_{\sun}$ yr$^{-1}$ Mpc$^{-3}$. The angular two-point correlation function of H$\alpha$ emitting galaxies of $\log L({\rm H\alpha}) > 39.8$ is well fit by a power law form of $w(\theta) = 0.013^{+0.002}_{-0.001} \theta^{-0.88 \pm 0.03}$, corresponding to the correlation function of $\xi(r) = (r/1.9{\rm Mpc})^{-1.88}$. We also find that the H$\alpha$ emitters with higher H$\alpha$ luminosity are more strongly clustered than those with lower luminosity.

We present imaging data and photometry for the COSMOS survey in 15 photometric bands between 0.3um and 2.4um. These include data taken on the Subaru 8.3m telescope, the KPNO and CTIO 4m telescopes, and the CFHT 3.6m telescope. Special techniques are used to ensure that the relative photometric calibration is better than 1% across the field of view. The absolute photometric accuracy from standard star measurements is found to be 6%. The absolute calibration is corrected using galaxy spectra, providing colors accurate to 2% or better. Stellar and galaxy colors and counts agree well with the expected values. Finally, as the first step in the scientific analysis of these data we construct panchromatic number counts which confirm that both the geometry of the universe and the galaxy population are evolving.

We present measurements of the two-point galaxy angular correlation function w(\theta) in the COSMOS field. Independent determinations of w(\theta) as a function of magnitude limit are presented for both the HST ACS catalog and also for the ground-based data from Subaru and the CFHT. Despite having significantly different masks, these three determinations agree well. At bright magnitudes (IAB<22), our data generally match very well with existing measurements and with mock catalogs based on semi-analytic galaxy formation calculations of Kitzbichler and White from the Millennium Simulation. The exception is that our result is at the upper end of the expected cosmic variance scatter for \theta > 10 arcmin, which we attribute to a particularly rich structure known to exist at z~0.8. For fainter samples, however, the level of clustering is somewhat higher than reported by some previous studies: in all three catalogues we find w(\theta=1')~0.014 at a median IAB magnitude of 24. At these very faintest magnitudes, our measurements agree well with the latest determinations from the Canada-France Legacy Survey. This level of clustering is approximately double what is predicted by the semi-analytic catalogs (at all angles). The semi-analytic results allow an estimate of cosmic variance, which is too small to account for the discrepancy. We therefore conclude that the mean amplitude of clustering at this level is higher than previously estimated.

We present a weak gravitational lensing analysis of 22 early-type strong lens galaxies, based on deep HST images obtained as part of the Sloan Lens ACS Survey. Using the most advanced techniques to control systematic uncertainties related to the variable PSF and charge transfer efficiency of the ACS, we detect weak lensing signal out to 300 kpc/h. We analyze blank control fields from the COSMOS survey in the same manner, inferring that the residual systematic uncertainty in the tangential shear is <0.3%. A joint strong and weak lensing analysis shows that the average total mass density profile is consistent with isothermal over two decades in radius (3-300 kpc/h, approximately 1-100 Reff). This finding extends by over an order of magnitude in radius previous results, based on strong lensing and/or stellar dynamics, that luminous and dark component ``conspire'' to form an isothermal mass distribution. In order to disentangle the contributions of luminous and dark matter, we fit a two-component mass model (R^1/4 + NFW) to the weak and strong lensing constraints. It provides a good fit to the data with only two free parameters; i) the average stellar mass-to-light ratio M_*/L_V=4.48 +- 0.46 hMo/Lo, in agreement with that expected for an old stellar population; ii) the average virial mass-to-light ratio M_vir/L_V = 246+101-87 hMo/Lo. [abridged]

Weak gravitational lensing is a promising probe of dark matter and dark energy requiring accurate measurement of the shapes of faint, distant galaxies. Such measures are hindered by the finite resolution and pixel scale of typical cameras. On the other hand, as imaging telescopes are practically limited to a fixed number of pixels and operational life-span, so the survey area increases with pixel size. We investigate the optimum choice of pixel scale in this trade-off for a space-based mission, using the full engineering model and survey strategy of the proposed SuperNova/Acceleration Probe as an example. Our methodology is to simulate realistic astronomical images of known shear and to evaluate the surface density of sources where the shear is accurately recovered using the Rhodes, Refregier, & Groth algorithm in the context of the derived dark matter power spectrum. In addition to considering single exposures, we also examine the benefits of sub-pixel dithering. Although some of our results depend upon the adopted shape measurement method, the relative trends, particularly those involving the surface density of resolved galaxies, are robust. Our approach provides a practical counterpart to studies which consider the effects of pixelation from analytic principles, which necessarily assume an idealized shape measurement method. We find that the statistical error on the mass power spectrum is minimized with a pixel scale equal to 75-80% of the FWHM of the point-spread function, and that dithering is marginally beneficial at larger pixel scales.

We describe the details of the Hubble Space Telescope (HST) Advanced Camera for Surveys / Wide Field Channel (ACS/WFC) observations of the COSMOS field, including the data calibration and processing procedures. We obtained a total of 583 orbits of HST ACS/WFC imaging in the F814W filter, covering a field that is 1.64 square degrees in area, the largest contiguous field ever imaged with HST. The median exposure depth across the field is 2028 seconds (one HST orbit), achieving a limiting point-source depth AB(F814W) = 27.2 (5 sigma). We also present details about the astrometric image registration, distortion removal and image combination using MultiDrizzle, as well as motivating the choice of our final pixel scale (30 milliarcseconds per pixel), based on the requirements for weak lensing science. The final set of images are publicly available through the archive sites at IPAC and STScI, along with further documentation on how they were produced.

We present the first identification of large-scale structures (LSS) at z $< 1.1$ in the Cosmic Evolution Survey (COSMOS). The structures are identified from adaptive smoothing of galaxy counts in the pseudo-3d space ($\alpha,\delta$,z) using the COSMOS photometric redshift catalog. The technique is tested on a simulation including galaxies distributed in model clusters and a field galaxy population -- recovering structures on all scales from 1 to 20\arcmin without {\it a priori} assumptions for the structure size or density profile. Our procedure makes {\bf no} {\it a priori} selection on galaxy spectral energy distribution (SED, for example the Red Sequence), enabling an unbiased investigation of environmental effects on galaxy evolution. The COSMOS photometric redshift catalog yields a sample of $1.5\times10^5$ galaxies with redshift accuracy, $\Delta z_{FWHM}/(1+z) \leq 0.1$ at z $< 1.1$ down to I$_{AB} \leq 25$ mag. Using this sample of galaxies, we identify 42 large-scale structures and clusters. abstract truncated for astroph 25 line limit -- see preprint

The Cosmic Evolution Survey (COSMOS) was initiated with an extensive allocation (590 orbits in Cycles 12-13) using the Hubble Space Telescope (HST) for high resolution imaging. Here we review the characteristics of the HST imaging with the Advanced Camera for Surveys (ACS) and parallel observations with NICMOS and WFPC2. A square field (1.8$\sq$\deg) has been imaged with single-orbit ACS I-F814W exposures with 50% completeness for sources 0.5\arcsec in diameter at I$_{AB} $ = 26.0 mag. The ACS imaging is a key part of the COSMOS survey, providing very high sensitivity and high resolution (0.09\arcsec FWHM, 0.05\arcsec pixels) imaging and detecting 1.2 million objects to a limiting magnitude of 26.5 (AB). These images yield resolved morphologies for several hundred thousand galaxies. The small HST PSF also provides greatly enhanced sensitivity for weak lensing investigations of the dark matter distribution.

With a primary goal of conducting precision weak lensing measurements from space, the COSMOS survey has imaged the largest contiguous area observed by the Hubble Space Telescope (HST) to date using the Advanced Camera for Surveys (ACS). This is the first paper in a series where we describe our strategy for addressing the various technical challenges in the production of weak lensing measurements from the COSMOS data. The COSMOS ACS catalog is constructed from 575 ACS/WFC tiles (1.64 deg^2) and contains a total 1.2x10^6 objects to a limiting magnitude of F814W=26.5. This catalog is made publicly available. The shapes of galaxies have been measured and corrected for the distortion induced by the time varying ACS Point Spread Function and for Charge Transfer Efficiency effects. Next, simulated images are used to derive the shear susceptibility factors that are necessary in order to transform shape measurements into unbiased shear estimators. Finally, for each galaxy, we derive a shape measurement error and utilize this quantity to extract the intrinsic shape noise of the galaxy sample. Interestingly, our results indicate that the intrinsic shape noise varies little with either size, magnitude or redshift. Representing a number density of 66 galaxies per arcmin^2, the final COSMOS weak lensing catalog contains 3.9x10^5 galaxies with accurate shape measurements. The properties of the COSMOS weak lensing catalog described throughout this paper will provide key input numbers for the preparation and design of next-generation wide field space missions.

ABRIDGED: We study the evolution since z~1 of the rest-frame B luminosity function of the early-type galaxies (ETGs) in ~0.7 deg^2 in the COSMOS field. In order to identify ALL progenitors of local ETGs we construct the sample of high-z galaxies using two complementary criteria: (i) A morphological selection based on the Zurich Estimator of Structural Types, and (ii) A photometric selection based on the galaxy properties in the (U-V)-M_V color-magnitude diagram. We furthermore constrain both samples so as to ensure that the selected progenitors of ETGs are compatible with evolving into systems which obey the mu_B-r_{hl} Kormendy relation. Assuming the luminosity evolution derived from studies of the fundamental plane for high-z ETGs, our analysis shows no evidence for a decrease in the number density of the most massive ETGs out to z~ 0.7: Both the morphologically- and the photometrically-selected sub-samples show no evolution in the number density of bright (~L>2.5L*) ETGs. Allowing for different star formation histories, and cosmic variance, we estimate a maximum decrease in the number density of massive galaxies at that redshift of ~30%. We observe, however, in both the photometrical and morphological samples, a deficit of up to ~2-3 of fainter ETGs over the same cosmic period. Our results argue against a significant contribution of recent dissipationless ``dry'' mergers to the formation of the most massive ETGs. We suggest that the mass growth in low luminosity ETGs can be explained with a conversion from z~0.7 to z=0 of blue, irregular and disk galaxies into low- and intermediate-mass ``red'' ETGs, possibly also through gas rich mergers.

(abridged) We examine the spatial and temporal stability of the HST ACS Wide Field Camera (WFC) point spread function (PSF) using the two square degree COSMOS survey. We show that stochastic aliasing of the PSF necessarily occurs during `drizzling'. This aliasing is maximal if the output pixel scale is equal to the input pixel scale of 0.05''. We show that this source of PSF variation can be significantly reduced by choosing a Gaussian drizzle kernel and by setting the output pixel size to 0.03''. We show that the PSF is temporally unstable, most likely due to thermal fluctuations in the telescope's focus. We find that the primary manifestation of this thermal drift in COSMOS images is an overall slow periodic focus change. Using a modified version of TinyTim, we create undistorted stars in a 30x30 grid across the ACS WFC CCDs. These PSF models are created for telescope focus values in the range -10 microns to +5 microns, thus spanning the allowed range of telescope focus values. We then use the approximately ten well measured stars in each COSMOS field to pick the best-fit focus value for each field. The TinyTim model stars are then used to perform PSF corrections for weak lensing allowing systematics due to incorrectly modeled PSFs to be greatly reduced. We have made the software for PSF modeling using our modified version of TinyTim available to the astronomical community. We show the effects of Charge Transfer Efficiency (CTE) degradation, which distorts the object in the readout direction, mimicking a weak lensing signal. We derive a parametric correction for the effect of CTE on the shapes of objects in the COSMOS field as a function of observation date, position within the ACS WFC field, and object flux.

We present a three dimensional cosmic shear analysis of the Hubble Space Telescope COSMOS survey, the largest ever optical imaging program performed in space. We have measured the shapes of galaxies for the tell-tale distortions caused by weak gravitational lensing, and traced the growth of that signal as a function of redshift. Using both 2D and 3D analyses, we measure cosmological parameters Omega_m, the density of matter in the universe, and sigma_8, the normalization of the matter power spectrum. The introduction of redshift information tightens the constraints by a factor of three, and also reduces the relative sampling (or "cosmic") variance compared to recent surveys that may be larger but are only two dimensional. From the 3D analysis, we find sigma_8*(Omega_m/0.3)^0.44=0.866+^0.085_-0.068 at 68% confidence limits, including both statistical and potential systematic sources of error in the total budget. Indeed, the absolute calibration of shear measurement methods is now the dominant source of uncertainty. Assuming instead a baseline cosmology to fix the geometry of the universe, we have measured the growth of structure on both linear and non-linear physical scales. Our results thus demonstrate a proof of concept for tomographic analysis techniques that have been proposed for future weak lensing surveys by a dedicated wide-field telescope in space.

Ordinary baryonic particles (such as protons and neutrons) account for only one-sixth of the total matter in the Universe. The remainder is a mysterious "dark matter" component, which does not interact via electromagnetism and thus neither emits nor reflects light. As dark matter cannot be seen directly using traditional observations, very little is currently known about its properties. It does interact via gravity, and is most effectively probed through gravitational lensing: the deflection of light from distant galaxies by the gravitational attraction of foreground mass concentrations. This is a purely geometrical effect that is free of astrophysical assumptions and sensitive to all matter -- whether baryonic or dark. Here we show high fidelity maps of the large-scale distribution of dark matter, resolved in both angle and depth. We find a loose network of filaments, growing over time, which intersect in massive structures at the locations of clusters of galaxies. Our results are consistent with predictions of gravitationally induced structure formation, in which the initial, smooth distribution of dark matter collapses into filaments then into clusters, forming a gravitational scaffold into which gas can accumulate, and stars can be built.

We have identified a large-scale structure at z~0.73 in the COSMOS field, coherently described by the distribution of galaxy photometric redshifts, an ACS weak-lensing convergence map and the distribution of extended X-ray sources in a mosaic of XMM observations. The main peak seen in these maps corresponds to a rich cluster with Tx= 3.51+0.60/-0.46 keV and Lx=(1.56+/-0.04) x 10^{44} erg/s ([0.1-2.4] keV band). We estimate an X-ray mass within $r500$ corresponding to M500~1.6 x 10^{14} Msun and a total lensing mass (extrapolated by fitting a NFW profile) M(NFW)=(6+/-3) x 10^15 Msun. We use an automated morphological classification of all galaxies brighter than I_AB=24 over the structure area to measure the fraction of early-type objects as a function of local projected density Sigma_10, based on photometric redshifts derived from ground-based deep multi-band photometry. We recover a robust morphology-density relation at this redshift, indicating, for comparable local densities, a smaller fraction of early-type galaxies than today. Interestingly, this difference is less strong at the highest densities and becomes more severe in intermediate environments. We also find, however, local "inversions'' of the observed global relation, possibly driven by the large-scale environment. In particular, we find direct correspondence of a large concentration of disk galaxies to (the colder side of) a possible shock region detected in the X-ray temperature map and surface brightness distribution of the dominant cluster. We interpret this as potential evidence of shock-induced star formation in existing galaxy disks, during the ongoing merger between two sub-clusters.

[ABRIGED] We present the optical identification of a sample of 695 X-ray sources detected in the first 1.3 deg^2 of the XMM-COSMOS survey, down to a 0.5-2 keV (2-10 keV) limiting flux of ~10^-15 erg cm-2 s-1 (~5x10^-15 erg cm^-2 s-1). We were able to associate a candidate optical counterpart to ~90% (626) of the X-ray sources, while for the remaining ~10% of the sources we were not able to provide a unique optical association due to the faintness of the possible optical counterparts (I_AB>25) or to the presence of multiple optical sources within the XMM-Newton error circles. We also cross-correlated the candidate optical counterparts with the Subaru multicolor and ACS catalogs and with the Magellan/IMACS, zCOSMOS and literature spectroscopic data; the spectroscopic sample comprises 248 objects (~40% of the full sample). Our analysis reveals that for ~80% of the counterparts there is a very good agreement between the spectroscopic classification, the morphological parameters as derived from ACS data, and the optical to near infrared colors. About 20% of the sources show an apparent mismatch between the morphological and spectroscopic classifications. All the ``extended'' BL AGN lie at redshift <1.5, while the redshift distribution of the full BL AGN population peaks at z~1.5. Our analysis also suggests that the Type 2/Type 1 ratio decreases towards high luminosities, in qualitative agreement with the results from X-ray spectral analysis and the most recent modeling of the X-ray luminosity function evolution.

(ABRIDGED) Motivated by the desire to reliably and automatically classify structure of thousands of COSMOS galaxies, we present ZEST, the Zurich Estimator of Structural Types. To classify galaxy structure, ZEST uses: (i) Five non-parametric diagnostics: asymmetry, concentration, Gini coefficient, 2nd-order moment of the brightest 20% of galaxy pixels, and ellipticity; and (ii) The exponent n of single--Sersic fits to the 2D surface brightness distributions. To fully exploit the wealth of information while reducing the redundancy present in these diagnostics, ZEST performs a principal component (PC) Analysis. We use a sample of ~56,000 I<24 COSMOS galaxies to show that the first three PCs fully describe the key aspects of the galaxy structure, i.e., to calibrate a three-dimensional classification grid of axis PC_1, PC_2, and PC_3. We demonstrate the robustness of the ZEST grid on the z=0 sample of Frei et al. (1996). The ZEST classification breaks most of the degeneracy between different galaxy populations that affects morphological classifications based on only some of the diagnostics included in ZEST. As a first application, we present the evolution since z~1 of the Luminosity Functions of COSMOS galaxies of early, disk and irregular galaxies and, for disk galaxies, of different bulge-to-disk ratios. Overall, we find that the LF up to a redshift z=1 is consistent with a pure-luminosity evolution (of about 0.95 magnitudes at z \~0.7). We highlight however two trends, that are in general agreement with a down-sizing scenario for galaxy formation: (1.) A deficit of a factor of about two at z~0.7 of MB>-20.5 structurally--classified early--type galaxies; and (2.) An excess of a factor of about three, at a similar redshift, of irregular galaxies.

The Shear TEsting Programme (STEP) is a collaborative project to improve the accuracy and reliability of weak lensing measurement, in preparation for the next generation of wide-field surveys. We review sixteen current and emerging shear measurement methods in a common language, and assess their performance by running them (blindly) on simulated images that contain a known shear signal. We determine the common features of algorithms that most successfully recover the input parameters. We achieve previously unattained discriminatory precision in our analysis, via a combination of more extensive simulations, and pairs of galaxy images that have been rotated with respect to each other, thus removing noise from their intrinsic ellipticities. The robustness of our simulation approach is also confirmed by testing the relative calibration of methods on real data. Weak lensing measurement has improved since the first STEP paper. Several methods now consistently achieve better than 2% precision, and are still being developed. However, the simulations can now distinguish all methods from perfect performance. Our main concern continues to be the potential for a multiplicative shear calibration bias: not least because this can not be internally calibrated with real data. We determine which galaxy populations are responsible and, by adjusting the simulated observing conditions, we also investigate the effects of instrumental and atmospheric parameters. We have isolated several previously unrecognised aspects of galaxy shape measurement, in which focussed development could provide further progress towards the sub-percent level of precision desired for future surveys. [ABRIDGED]

We study a sample of approximately 16,500 galaxies with I_AB <= 22.5 in the COSMOS field. Structural information on the galaxies is derived by fitting single Sersic models to their two-dimensional surface brightness distributions. We investigate the evolution of the number density of disk galaxies larger than 5 kpc between redshift z~1 and the present epoch. To this end, we use the measurements of the half-light radii to construct, as a function of redshift, the size function of both the total disk galaxy population and of disk galaxies split in four bins of bulge-to-disk ratio. Furthermore, we use a selected sample of roughly 1800 SDSS galaxies to calibrate our results with respect to the local universe. We find that: (i) The number density of disk galaxies with intermediate sizes (r_{1/2}~5-7 kpc) remains nearly constant from z~1 to today. (ii) The number density of the largest disks (r_{1/2}>7 kpc) decreases by a factor of about two out to z~1. (iii) There is a constancy in the number density of large bulgeless disks out to z~1; the deficit of large disks at early epochs seems to arise from a smaller number of bulged disks. Our results indicate that the bulk of the large disk galaxy population has completed its growth by z~1, and support the hypothesis that secular evolution processes produce - or at least add stellar mass to - the bulge components of disk galaxies.

We propose a tunable laser-based satellite-mounted spectrophotometric and absolute flux calibration system, to be utilized by ground- and space-based telescopes. As spectrophotometric calibration may play a significant role in the accuracy of photometric redshift measurement, and photometric redshift accuracy is important for measuring dark energy using SNIa, weak gravitational lensing, and baryon oscillations, a method for reducing such uncertainties is needed. We propose to improve spectrophotometric calibration, currently obtained using standard stars, by placing a tunable laser and a wide-angle light source on a satellite by early next decade (perhaps included in the upgrade to the GPS satellite network) to improve absolute flux calibration and relative spectrophotometric calibration across the visible and near-infrared spectrum. As well as fundamental astrophysical applications, the system proposed here potentially has broad utility for defense and national security applications such as ground target illumination and space communication.

The ability to accurately measure the shapes of faint objects in images taken with the Advanced Camera for Surveys(ACS) on the Hubble Space Telescope (HST) depends upon detailed knowledge of the Point Spread Function (PSF). We show that thermal fluctuations cause the PSF of the ACS Wide Field Camera (WFC) to vary over time. We describe a modified version of the TinyTim PSF modeling software to create artificial grids of stars across the ACS field of view at a range of telescope focus values. These models closely resemble the stars in real ACS images. Using ~10 bright stars in a real image, we have been able to measure HST's apparent focus at the time of the exposure. TinyTim can then be used to model the PSF at any position on the ACS field of view. This obviates the need for images of dense stellar fields at different focus values, or interpolation between the few observed stars. We show that residual differences between our TinyTim models and real data are likely due to the effects of Charge Transfer Efficiency (CTE) degradation. Furthermore, we discuss stochastic noise that is added to the shape of point sources when distortion is removed, and we present Multidrizzle parameters that are optimal for weak lensing science. Specifically, we find that reducing the Multidrizzle output pixel scale and choosing a Gaussian kernel significantly stabilizes the resulting PSF after image combination, while still eliminating cosmic rays/bad pixels, and correcting the large geometric distortion in the ACS. We discuss future plans, which include more detailed study of the effects of CTE degradation on object shapes and releasing our TinyTim models to the astronomical community.

The Shear TEsting Programme, STEP, is a collaborative project to improve the accuracy and reliability of all weak lensing measurements in preparation for the next generation of wide-field surveys. In this first STEP paper we present the results of a blind analysis of simulated ground-based observations of relatively simple galaxy morphologies. The most successful methods are shown to achieve percent level accuracy. From the cosmic shear pipelines that have been used to constrain cosmology, we find weak lensing shear measured to an accuracy that is within the statistical errors of current weak lensing analyses, with shear measurements accurate to better than 7%. The dominant source of measurement error is shown to arise from calibration uncertainties where the measured shear is over or under-estimated by a constant multiplicative factor. This is of concern as calibration errors cannot be detected through standard diagnostic tests. The measured calibration errors appear to result from stellar contamination, false object detection, the shear measurement method itself, selection bias and/or the use of biased weights. Additive systematics (false detections of shear) resulting from residual point-spread function anisotropy are, in most cases, reduced to below an equivalent shear of 0.001, an order of magnitude below cosmic shear distortions on the scales probed by current surveys. Our results provide a snapshot view of the accuracy of current ground-based weak lensing methods and a benchmark upon which we can improve. To this end we provide descriptions of each method tested and include details of the eight different implementations of the commonly used Kaiser, Squires and Broadhurst (1995) method (KSB+) to aid the improvement of future KSB+ analyses.

The Deep Extragalactic Exploratory Probe (DEEP) is a multi-phase research program dedicated to the study of the formation and evolution of galaxies and of large scale structure in the distant Universe. This paper describes the first five-year phase, denoted DEEP1. A series of ten DEEP1 papers will discuss a range of scientific topics (e.g., the study of photometric and spectral properties of a general distant galaxy survey, the evolution observed in galaxy populations of varied morphologies). The observational basis for these studies is the Groth Survey Strip field, a 127 square arcminute region which has been observed with the Hubble Space Telescope in both broad I-band and V-band optical filters and with the Low Resolution Imaging Spectrograph on the Keck Telescopes. Catalogs of photometric and structural parameters have been constructed for 11,547 galaxies and stars at magnitudes brighter than 29, and spectroscopy has been conducted for a magnitude-color weighted subsample of 818 objects. We evaluate three independent techniques for constructing an imaging catalog for the field from the HST data, and discuss the depth and sampling of the resultant catalogs. The selection of the spectroscopic subsample is discussed, and we describe the multifaceted approach taken to prioritizing objects of interest for a variety of scientific subprograms. A series of Monte Carlo simulations then demonstrates that the spectroscopic subsample can be adequately modeled as a simple function of magnitude and color cuts in the imaging catalog.

The Deep Extragalactic Evolutionary Probe (DEEP) is a series of spectroscopic surveys of faint galaxies, targeted at the properties and clustering of galaxies at redshifts z ~ 1. We present the redshift catalog of the DEEP 1 GSS pilot phase of this project, a Keck/LRIS survey in the HST/WFPC2 Groth Survey Strip. The redshift catalog and data, including reduced spectra, are publicly available through a Web-accessible database. The catalog contains 658 secure galaxy redshifts with a median z=0.65, and shows large-scale structure walls to z = 1. We find a bimodal distribution in the galaxy color-magnitude diagram which persists to z = 1. A similar color division has been seen locally by the SDSS and to z ~ 1 by COMBO-17. For red galaxies, we find a reddening of only 0.11 mag from z ~ 0.8 to now, about half the color evolution measured by COMBO-17. We measure structural properties of the galaxies from the HST imaging, and find that the color division corresponds generally to a structural division. Most red galaxies, ~ 75%, are centrally concentrated, with a red bulge or spheroid, while blue galaxies usually have exponential profiles. However, there are two subclasses of red galaxies that are not bulge-dominated: edge-on disks and a second category which we term diffuse red galaxies (DIFRGs). The distant edge-on disks are similar in appearance and frequency to those at low redshift, but analogs of DIFRGs are rare among local red galaxies. DIFRGs have significant emission lines, indicating that they are reddened mainly by dust rather than age. The DIFRGs in our sample are all at z>0.64, suggesting that DIFRGs are more prevalent at high redshifts; they may be related to the dusty or irregular extremely red objects (EROs) beyond z>1.2 that have been found in deep K-selected surveys. (abridged)

We present the results of a 200 ks Chandra observation of part of the Groth Strip region, using the ACIS-I instrument. We present a relatively simple method for the detection of point-sources and calculation of limiting sensitivities, which we argue is at least as sensitive and more self-consistent than previous methods presented in the literature. 158 distinct X-ray sources are included in our point-source catalogue in the ACIS-I area. The number counts show a relative dearth of X-ray sources in this region. For example at a flux limit of 1E-15 (cgs) around 20 per cent more soft band sources are detected in the HDF-N and almost 50 per cent more in the ELAIS-N1 field, which we have analysed by the same method for comparison. We find, however, that these differences are consistent with Poisson variations at 2 sigma significance, and therefore there is no evidence for cosmic variance based on these number counts alone. We determine the average spectra of the objects and find a marked difference between the soft-band selected sources, which have Gamma=1.9 typical of unobscured AGN, and the hard-band selected sources, which have Gamma=1.0. Reassuringly, the sample as a whole has a mean spectrum of Gamma=1.4\pm 0.1, the same as the X-ray background. Nonetheless, our results imply that the fraction of sources with significant obscuration is only ~25 per cent, much less than predicted by standard AGN population synthesis models. This is confirmed by direct spectral fitting, with only a handful of objects showing evidence for absorption. After accounting for absorption, all objects are consistent with mean intrinsic spectrum of Gamma=1.76 \pm 0.08, very similar to local Seyferts (abridged).

The Supernova / Acceleration Probe (SNAP) is a proposed space-based experiment designed to study the dark energy and alternative explanations of the acceleration of the Universe's expansion by performing a series of complementary systematics-controlled measurements. We describe a self-consistent reference mission design for building a Type Ia supernova Hubble diagram and for performing a wide-area weak gravitational lensing study. A 2-m wide-field telescope feeds a focal plane consisting of a 0.7 square-degree imager tiled with equal areas of optical CCDs and near infrared sensors, and a high-efficiency low-resolution integral field spectrograph. The SNAP mission will obtain high-signal-to-noise calibrated light-curves and spectra for several thousand supernovae at redshifts between z=0.1 and 1.7. A wide-field survey covering one thousand square degrees resolves ~100 galaxies per square arcminute. If we assume we live in a cosmological-constant-dominated Universe, the matter density, dark energy density, and flatness of space can all be measured with SNAP supernova and weak-lensing measurements to a systematics-limited accuracy of 1%. For a flat universe, the density-to-pressure ratio of dark energy can be similarly measured to 5% for the present value w0 and ~0.1 for the time variation w'. The large survey area, depth, spatial resolution, time-sampling, and nine-band optical to NIR photometry will support additional independent and/or complementary dark-energy measurement approaches as well as a broad range of auxiliary science programs. (Abridged)

A wide field space-based imaging telescope is necessary to fully exploit the technique of observing dark matter via weak gravitational lensing. This first paper in a three part series outlines the survey strategies and relevant instrumental parameters for such a mission. As a concrete example of hardware design, we consider the proposed Supernova/Acceleration Probe (SNAP). Using SNAP engineering models, we quantify the major contributions to this telescope's Point Spread Function (PSF). These PSF contributions are relevant to any similar wide field space telescope. We further show that the PSF of SNAP or a similar telescope will be smaller than current ground-based PSFs, and more isotropic and stable over time than the PSF of the Hubble Space Telescope. We outline survey strategies for two different regimes - a ``wide'' 300 square degree survey and a ``deep'' 15 square degree survey that will accomplish various weak lensing goals including statistical studies and dark matter mapping.

Weak gravitational lensing provides a unique method to directly map the dark matter in the universe and measure cosmological parameters. Current weak lensing surveys are limited by the atmospheric seeing from the ground and by the small field of view of existing space telescopes. We study how a future wide-field space telescope can measure the lensing power spectrum and skewness, and set constraints on cosmological parameters. The lensing sensitivity was calculated using detailed image simulations and instrumental specifications studied in earlier papers in this series. For instance, the planned SuperNova/Acceleration Probe (SNAP) mission will be able to measure the matter density parameter Omega_m and the dark energy equation of state parameter w with precisions comparable and nearly orthogonal to those derived with SNAP from supernovae. The constraints degrade by a factor of about 2 if redshift tomography is not used, but are little affected if the skewness only is dropped. We also study how the constraints on these parameters depend upon the survey geometry and define an optimal observing strategy.

We study the accuracy with which weak lensing measurements could be made from a future space-based survey, predicting the subsequent precisions of 3-dimensional dark matter maps, projected 2-dimensional dark matter maps, and mass-selected cluster catalogues. As a baseline, we use the instrumental specifications of the Supernova/Acceleration Probe (SNAP) satellite. We first compute its sensitivity to weak lensing shear as a function of survey depth. Our predictions are based on detailed image simulations created using `shapelets', a complete and orthogonal parameterization of galaxy morphologies. We incorporate a realistic redshift distribution of source galaxies, and calculate the average precision of photometric redshift recovery using the SNAP filter set to be Delta z=0.034. The high density of background galaxies resolved in a wide space-based survey allows projected dark matter maps with a rms sensitivity of 3% shear in 1 square arcminute cells. This will be further improved using a proposed deep space-based survey, which will be able to detect isolated clusters using a 3D lensing inversion techniques with a 1 sigma mass sensitivity of approximately 10^13 solar masses at z~0.25. Weak lensing measurements from space will thus be able to capture non-Gaussian features arising from gravitational instability and map out dark matter in the universe with unprecedented resolution.

Weak lensing provides a direct measure of the distribution of mass in the universe, and is therefore a uniquely powerful probe of dark matter. Weak lensing can also be used to measure the twin phenomenon of dark energy, via its effect upon the cosmological growth rate of structures. Essential for this technique are well-resolved images of background galaxies out to large distances. As a concrete example of the surveys that will become available by the end of the decade, we consider the planned Supernova/Acceleration Probe (SNAP) space telescope. Detailed simulations of space-based images, manufactured using the shapelets formalism, enable us to quantitatively predict the future sensitivity to weak lensing shear. The high number density of galaxies resolved from space will enable maps of dark matter to be produced in two and three dimensions, with a resolution superior to that from the ground. Such observations will also afford reduced systematics for high-precision measurements of weak lensing statistics. These will be used to set tight constraints on cosmological parameters. In particular, the parameter for equation of state of dark energy, w, will be measured using weak lensing with a precision comparable to and somewhat orthogonal to constraints from other methods.

Weak lensing by large-scale structure allows a direct measure of the dark matter distribution. We have used parallel images taken with the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope to measure weak lensing, or cosmic shear. We measure the shapes of 26036 galaxies in 1292 STIS fields and measure the shear variance at a scale of 0.51 arcminutes. The charge transfer efficiency (CTE) of STIS has degraded over time and introduces a spurious ellipticity into galaxy shapes during the readout process. We correct for this effect as a function of signal to noise and CCD position. We further show that the detected cosmic shear signal is nearly constant in time over the approximately four years of observation. We detect cosmic shear at the 5.1 sigma level, and our measurement of the shear variance is consistent with theoretical predictions in a LambdaCDM universe. This provides a measure of the normalization of the mass power spectrum sigma_8=(1.02 +- 0.16) (0.3/Omega_m)^{0.46} (0.21/Gamma)^{0.18}$. The one-sigma error includes noise, cosmic variance, systematics and the redshift uncertainty of the source galaxies. This is consistent with previous cosmic shear measurements, but tends to favor those with a high value of sigma_8. It is also consistent with the recent determination of sigma_8 from the Wilkinson Microwave Anisotropy Probe (WMAP) experiment.

We present the luminosity function of [OII]-emitting galaxies at a median redshift of z=0.9, as measured in the deep spectroscopic data in the STIS Parallel Survey (SPS). The luminosity function shows strong evolution from the local value, as expected. By using random lines of sight, the SPS measurement complements previous deep single field studies. We calculate the density of inferred star formation at this redshift by converting from [OII] to H-alpha line flux as a function of absolute magnitude and find rho_dot=0.043 +/- 0.014 Msun/yr/Mpc^3 at a median redshift z~0.9 within the range 0.46<z<1.415 (H_0 = 70 km/s/Mpc, Omega_M=0.3, Omega_Lambda=0.7. This density is consistent with a (1+z)^4 evolution in global star formation since z~1. To reconcile the density with similar measurements made by surveys targeting H-alpha may require substantial extinction correction.

In the first three years of operation STIS obtained slitless spectra of approximately 2500 fields in parallel to prime HST observations as part of the STIS Parallel Survey (SPS). The archive contains almost 300 fields at high galactic latitude (|b|>30) with spectroscopic exposure times greater than 3000 seconds. This sample contains 220 fields (excluding special regions and requiring a consistent grating angle) observed between 6 June 1997 and 21 September 2000, with a total survey area of about 160 square arcminutes. At this depth, the SPS detects an average of one emission line galaxy per three fields. We present the analysis of these data, and the identification of 131 low to intermediate redshift galaxies detected by optical emission lines. The sample contains 78 objects with emission lines that we infer to be redshifted [OII]3727 emission at 0.43<z<1.7. The comoving number density of these objects is comparable to that of H-alpha emitting galaxies in the NICMOS parallel observations. One quasar and three probable Seyfert galaxies are detected. Many of the emission-line objects show morphologies suggestive of mergers or interactions. The reduced data are available upon request from the authors.

The Galactic Exoplanet Survey Telescope (GEST) will observe a 2 square degree field in the Galactic bulge to search for extra-solar planets using a gravitational lensing technique. This gravitational lensing technique is the only method employing currently available technology that can detect Earth-mass planets at high signal-to-noise, and can measure the frequency of terrestrial planets as a function of Galactic position. GEST's sensitivity extends down to the mass of Mars, and it can detect hundreds of terrestrial planets with semi-major axes ranging from 0.7 AU to infinity. GEST will be the first truly comprehensive survey of the Galaxy for planets like those in our own Solar System.

Weak lensing by large-scale structure provides a direct measurement of matter fluctuations in the universe. We report a measurement of this `cosmic shear' based on 271 WFPC2 archival images from the Hubble Space Telescope Medium Deep Survey (MDS). Our measurement method and treatment of systematic effects were discussed in an earlier paper. We measure the shear variance on scales ranging from 0.7' to 1.4', with a detection significance greater than 3.8. This allows us to measure the normalization of the matter power spectrum to be sigma8 = (0.94 +/- 0.10 +/- 0.14) (0.3/Omega_m)^0.44 (0.21/Gamma)^0.15, in a LCDM universe. The first 1sigma error includes statistical errors only, while the latter also includes (gaussian) cosmic variance and the uncertainty in the galaxy redshift distribution. Our results are consistent with earlier cosmic shear measurements from the ground and from space. We compare our cosmic shear results and those from other groups to the normalization from cluster abundance and galaxy surveys. We find that the combination of four recent cosmic shear measurements are somewhat inconsistent with the recent normalization using these methods, and discuss possible explanations for the discrepancy.

Weak lensing by large-scale structure provides a unique method to directly measure matter fluctuations in the universe, and has recently been detected from the ground. Here, we report the first detection of this `cosmic shear' based on space-based images. The detection was derived from the Hubble Space Telescope (HST) Survey Strip (or Groth Strip), a 4' by 42' set of 28 contiguous WFPC2 pointings with I<27. The small size of the HST Point-Spread Function (PSF) affords both a lower statistical noise, and a much weaker sensitivity to systematic effects, a crucial limiting factor of cosmic shear measurements. Our method and treatment of systematic effects were discussed in an earlier paper (Rhodes, Refregier & Groth 2000). We measure an rms shear of 1.8% on the WFPC2 chip scale (1.27'), in agreement with the predictions of cluster-normalized CDM models. Using a Maximum Likelihood (ML) analysis, we show that our detection is significant at the 99.5% confidence level (CL), and measure the normalization of the matter power spectrum to be sigma8*Omega_m^(0.48) = 0.51 (+0.14,-0.17), in a LambdaCDM universe. These 68% CL errors include (Gaussian) cosmic variance, systematic effects and the uncertainty in the redshift distribution of the background galaxies. Our result is consistent with earlier lensing measurements from the ground, and with the normalization derived from cluster abundance. We discuss how our measurement can be improved with the analysis of a large number of independent WFPC2 fields.

The weak distortions produced by gravitational lensing in the images of background galaxies provide a method to measure directly the distribution of mass in the universe. However this technique requires high precision measurements of the lensing shear and cautious corrections for systematic effects. Kaiser, Squires, & Broadhurst (1995) proposed a method to calibrate the ellipticity-shear relation in the presence of Point Spread Function (PSF) anisotropies and camera distortions. We revisit the KSB method and show that both the PSF and the camera distortions can be corrected for using source moments, as opposed to ellipticities. We clarify the applicability of some of the approximations made in this method. We derive expressions for the corrections which only involve the galaxy moments. We derive an explicit relation between the shear and the average ellipticity. We discuss the shortcomings of the method, and test its validity using numerical simulations. As an application of the method, we repeat the analysis of the HST WFPC2 camera performed by Hoekstra et al. (1998). We confirm the presence of sizable (10%) PSF ellipticities at the edge of the WFPC2 chips. We also show that the PSF ellipticity varies by as much as 2% over time. We use these measurements to correct the shape of galaxies in the HST Survey Strip (``Groth'' Strip). By considering the dependence of the ellipticities on object size, we show that, after corrections, the residual systematic uncertainty for galaxies with radii greater than 0.15 arcsec, is about 0.4%, when averaged over each chip. We discuss how these results provide good prospects for measuring weak lensing by large-scale structure with deep HST surveys.