This paper presents multiwavelength photometric catalogues of the last two Hubble Frontier Fields (HFF), the massive galaxy clusters Abell 370 and RXC J2248.7-4431. The photometry ranges from imaging performed on the Hubble Space Telescope (HST) to ground based Very Large Telescope (VLT) and Spitzer/IRAC, in collaboration with the ASTRODEEP team, and using the ASTRODEEP pipeline. While the main purpose of this paper is to release the catalogues, we also perform, as a proof of concept, a brief analysis of z > 6 objects selected using drop-out method, as well as spectroscopically confirmed sources and multiple images in both clusters. While dropout methods yield a sample of high-z galaxies, the addition of longer wavelength data reveals that as expected the samples have substantial contamination at the ~30-45% level by dusty galaxies at lower redshifts. Furthermore, we show that spectroscopic redshifts are still required to unambiguously determine redshifts of multiply imaged systems. Finally, the now publicly available ASTRODEEP catalogues were combined for all HFFs and used to explore stellar properties of a large sample of 20,000 galaxies across a large photometric redshift range. The powerful magnification provided by the HFF clusters allows us an exploration of the properties of galaxies with intrinsic stellar masses as low as $M_* \gtrsim 10^7M_{\odot}$ and intrinsic star formation rates $\mbox{SFRs}\sim 0.1\mbox{-}1M_\odot/\mbox yr$ at z > 6.

Measurements of stellar properties of galaxies when the universe was less than one billion years old yield some of the only observational constraints of the onset of star formation. We present here the inclusion of \textit{Spitzer}/IRAC imaging in the spectral energy distribution fitting of the seven highest-redshift galaxy candidates selected from the \emph{Hubble Space Telescope} imaging of the Reionization Lensing Cluster Survey (RELICS). We find that for 6/8 \textit{HST}-selected $z\gtrsim8$ sources, the $z\gtrsim8$ solutions are still strongly preferred over $z\sim$1-2 solutions after the inclusion of \textit{Spitzer} fluxes, and two prefer a $z\sim 7$ solution, which we defer to a later analysis. We find a wide range of intrinsic stellar masses ($5\times10^6 M_{\odot}$ -- $4\times10^9$ $M_{\odot}$), star formation rates (0.2-14 $M_{\odot}\rm yr^{-1}$), and ages (30-600 Myr) among our sample. Of particular interest is Abell1763-1434, which shows evidence of an evolved stellar population at $z\sim8$, implying its first generation of star formation occurred just $< 100$ Myr after the Big Bang. SPT0615-JD, a spatially resolved $z\sim10$ candidate, remains at its high redshift, supported by deep \textit{Spitzer}/IRAC data, and also shows some evidence for an evolved stellar population. Even with the lensed, bright apparent magnitudes of these $z \gtrsim 8$ candidates (H = 26.1-27.8 AB mag), only the \textit{James Webb Space Telescope} will be able further confirm the presence of evolved stellar populations early in the universe.

We introduce the OSIRIS Lens-Amplified Survey (OLAS), a kinematic survey of gravitationally lensed galaxies at cosmic noon taken with Keck adaptive optics. In this paper we present spatially resolved spectroscopy and nebular emission kinematic maps for 17 star forming galaxies with stellar masses 8 < log($M_*$/$M_{\odot}$) < 9.8 and redshifts 1.2 < z < 2.3. OLAS is designed to probe the stellar mass ($M_*$) and specific star formation rate (sSFR) range where simulations suggest that stellar feedback is most effective at driving gaseous outflows that create galaxy-wide potential fluctuations which can generate dark matter cores. We compare our kinematic data with the trend between sSFR, $M_*$ and H$\alpha$ velocity dispersion, $\sigma$, from the Feedback In Realistic Environments (FIRE) simulations. Our observations reveal a correlation between sSFR and sigma at fixed $M_*$ that is similar to the trend predicted by simulations: feedback from star formation drives star-forming gas and newly formed stars into more dispersion dominated orbits. The observed magnitude of this effect is in good agreement with the FIRE simulations, in which feedback alters the central density profiles of low mass galaxies, converting dark matter cusps into cores over time. Our data support the scenario that stellar feedback drives gaseous outflows and potential fluctuations, which in turn drive dark matter core formation in dwarf galaxies.

The rest-frame UV emission from massive stars contains a wealth of information about the physical nature and conditions of star formation in galaxies. Using studies of the rest-frame UV, the past decade has witnessed the beginning of knowledge about the existence and properties of galaxies during the first few billion years after the Big Bang. This period of history corresponds to the formation of the first stars, the rapid formation of galaxy stellar populations, the reionization of the IGM, the production and dissemination of heavy elements, and the formation of the first black holes. Massive stars in these galaxies drive all of these events, and their light dominates the spectral energy distributions of galaxies. As we look to the 2020s, fundamental questions remain about the nature of these stellar populations and their evolution, from just before the peak of the cosmic star formation density (z~3), up to the epoch of reionization (z > 6). This next decade will provide transformative gains both in our ability to identify star-forming galaxies and accreting supermassive black holes at these early epochs with imaging surveys in the rest-frame UV (e.g., LSST, WFIRST). Ground-based, rest-frame UV spectroscopy on >20 m-class telescopes (e.g., GMT/TMT) offers the ability to investigate the astrophysical conditions in galaxies at the earliest cosmic times. This includes studies of the evolution in galaxy stellar populations, gas ionization (temperature, pressure), metallicity, and interstellar (and circumgalactic) gas kinematics and covering fractions. In this white paper, we describe the scientific prospects and the requirements for research in this area.

The epoch of reionization (6 < z < 10) marks the period in our universe when the first large galaxies grew to fruition, and began to affect the universe around them. Massive stars, and potentially accreting supermassive black holes, filled the universe with ionizing radiation, burning off the haze of neutral gas that had filled the intergalactic medium (IGM) since recombination (z~1000). The evolution of this process constrains key properties of these earliest luminous sources, thus observationally constraining reionization is a key science goal for the next decade. The measurement of Lyman-alpha emission from photometrically-identified galaxies is a highly constraining probe of reionization, as a neutral IGM will resonantly scatter these photons, reducing detectability. While significant work has been done with 8-10m telescopes, these observations require extremely large telescopes (ELTs); the flux limits available from today's 10m class telescopes are sufficient for only the brightest known galaxies (m < 26). Ultra-deep surveys with the Giant Magellan Telescope (GMT) and Thirty Meter Telescope (TMT) will be capable of detecting Lyman-alpha emission from galaxies 2-3 magnitudes fainter than today's deepest surveys. Wide-field fiber spectroscopy on the GMT combined with narrow-field AO-assisted slit spectroscopy on the TMT will be able to probe the expected size of ionized bubbles throughout the epoch of reionization, following up degree scale deep imaging surveys with the Wide Field Infrared Space Telescope. These data will provide the first resolved Lyman-alpha-based maps of the ionized intergalactic medium throughout the epoch of reionization, constraining models of both the temporal and spatial evolution of this phase change.

Large surveys of galaxy clusters with the Hubble and Spitzer Space Telescopes, including CLASH and the Frontier Fields, have demonstrated the power of strong gravitational lensing to efficiently deliver large samples of high-redshift galaxies. We extend this strategy through a wider, shallower survey named RELICS, the Reionization Lensing Cluster Survey. This survey, described here, was designed primarily to deliver the best and brightest high-redshift candidates from the first billion years after the Big Bang. RELICS observed 41 massive galaxy clusters with Hubble and Spitzer at 0.4-1.7um and 3.0-5.0um, respectively. We selected 21 clusters based on Planck PSZ2 mass estimates and the other 20 based on observed or inferred lensing strength. Our 188-orbit Hubble Treasury Program obtained the first high-resolution near-infrared images of these clusters to efficiently search for lensed high-redshift galaxies. We observed 46 WFC3/IR pointings (~200 arcmin^2) with two orbits divided among four filters (F105W, F125W, F140W, and F160W) and ACS imaging as needed to achieve single-orbit depth in each of three filters (F435W, F606W, and F814W). As previously reported by Salmon et al., we discovered 322 z ~ 6 - 10 candidates, including the brightest known at z ~ 6, and the most distant spatially-resolved lensed arc known at z ~ 10. Spitzer IRAC imaging (945 hours awarded, plus 100 archival) has crucially enabled us to distinguish z ~ 10 candidates from z ~ 2 interlopers. For each cluster, two HST observing epochs were staggered by about a month, enabling us to discover 11 supernovae, including 3 lensed supernovae, which we followed up with 20 orbits from our program. We delivered reduced HST images and catalogs of all clusters to the public via MAST and reduced Spitzer images via IRSA. We have also begun delivering lens models of all clusters, to be completed before the JWST GO call for proposals.

Detections and non-detections of Lyman alpha (Ly$\alpha$) emission from $z>6$ galaxies ($<1$ Gyr after the Big Bang) can be used to measure the timeline of cosmic reionization. Of key interest to measuring reionization's mid-stages, but also increasing observational challenge, are observations at z > 7, where Ly$\alpha$ redshifts to near infra-red wavelengths. Here we present a search for z > 7.2 Ly$\alpha$ emission in 53 intrinsically faint Lyman Break Galaxy candidates, gravitationally lensed by massive galaxy clusters, in the KMOS Lens-Amplified Spectroscopic Survey (KLASS). With integration times of ~7-10 hours, we detect no Ly$\alpha$ emission with S/N>5 in our sample. We determine our observations to be 80% complete for 5$\sigma$ spatially and spectrally unresolved emission lines with integrated line flux $>5.7\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$. We define a photometrically selected sub-sample of 29 targets at $z=7.9\pm0.6$, with a median 5$\sigma$ Ly$\alpha$ EW limit of 58A. We perform a Bayesian inference of the average intergalactic medium (IGM) neutral hydrogen fraction using their spectra. Our inference accounts for the wavelength sensitivity and incomplete redshift coverage of our observations, and the photometric redshift probability distribution of each target. These observations, combined with samples from the literature, enable us to place a lower limit on the average IGM neutral hydrogen fraction of $> 0.76 \; (68\%), \; > 0.46 \; (95\%)$ at z ~ 8, providing further evidence of rapid reionization at z~6-8. We show that this is consistent with reionization history models extending the galaxy luminosity function to $M_\textrm{UV} \lesssim -12$, with low ionizing photon escape fractions, $f_\textrm{esc} \lesssim 15\%$.

We present a new gravitational lens model of the Hubble Frontier Fields cluster Abell 370 ($z = 0.375$) using imaging and spectroscopy from Hubble Space Telescope and ground-based spectroscopy. We combine constraints from a catalog of 1344 weakly lensed galaxies and 39 multiply-imaged sources comprised of 114 multiple images, including a system of multiply-imaged candidates at $z=7.93 \pm 0.02$, to obtain a best-fit mass distribution using the cluster lens modeling code Strong and Weak Lensing United. As the only analysis of A370 using strong and weak lensing constraints from Hubble Frontier Fields data, our method provides an independent check on assumptions in other methods on the mass distribution. Convergence, shear, and magnification maps are made publicly available through the HFF website. We find that the model we produce is similar to models produced by other groups, with some exceptions due to the differences in lensing code methodology. In an effort to study how our total projected mass distribution traces light, we measure the stellar mass density distribution using Spitzer/Infrared Array Camera imaging. Comparing our total mass density to our stellar mass density in a radius of 0.3 Mpc, we find a mean projected stellar to total mass ratio of $\langle f* \rangle = 0.011 \pm 0.003$ (stat.) using the diet Salpeter initial mass function. This value is in general agreement with independent measurements of $\langle f* \rangle$ in clusters of similar total mass and redshift.

We present the results of CANDELSz7, an ESO large program aimed at confirming spectroscopically a homogeneous sample of z~6 and z~7 star forming galaxies. The candidates were selected in the GOODS-South, UDS and COSMOS fields using the official CANDELS catalogs based on H160-band detections. Standard color criteria, which were tailored depending on the ancillary multi-wavelength data available for each field, were applied to select more than 160 candidate galaxies at z~6 and z~7. Deep medium resolution FORS2 spectroscopic observations were then conducted with integration times ranging from 12 to 20 hours, to reach a Lyalpha flux limit of approximately 1-3x 10-18 erg/s/cm^2 at 3sigma. For about 40% of the galaxies we could determine a spectroscopic redshift, mainly through the detection of a single emission line that we interpret as Lyalpha emission, or for some of the brightest objects (H160< 25.5) from the presence of faint continuum and sharp drop that we interpret as a Lyman break. In this paper we present the redshifts and main properties of 65 newly confirmed high redshift galaxies. Adding previous proprietary and archival data we assemble a sample of ~260 galaxies that we use to explore the evolution of the Lyalpha fraction in Lyman break galaxies and the change in the shape of the emission line between z~6 and z~7. We also discuss the accuracy of the CANDELS photometric redshifts in this redshift range.

We present a lens model for the cluster SPT-CLJ0615$-$5746, which is the highest redshift ($z=0.972$) system in the Reionization of Lensing Clusters Survey (RELICS), making it the highest redshift cluster for which a full strong lens model is published. We identify three systems of multiply-imaged lensed galaxies, two of which we spectroscopically confirm at $z=1.358$ and $z=4.013$, which we use as constraints for the model. We find a foreground structure at $z\sim0.4$, which we include as a second cluster-sized halo in one of our models; however two different statistical tests find the best-fit model consists of one cluster-sized halo combined with three individually optimized galaxy-sized halos, as well as contributions from the cluster galaxies themselves. We find the total projected mass density within $r=26.7"$ (the region where the strong lensing constraints exist) to be $M=2.51^{+0.15}_{-0.09}\times 10^{14}$~M$_{\odot}$. If we extrapolate out to $r_{500}$, our projected mass density is consistent with the mass inferred from weak lensing and from the Sunyaev-Zel'dovich effect ($M\sim10^{15}$~M$_{\odot}$). This cluster is lensing a previously reported $z\sim10$ galaxy, which, if spectroscopically confirmed, will be the highest-redshift strongly lensed galaxy known.

The optical observations of Ic-4 supernova (SN) 2016coi/ASASSN-16fp, from $\sim 2$ to $\sim450$ days after explosion, are presented along with analysis of its physical properties. The SN shows the broad lines associated with SNe Ic-3/4 but with a key difference. The early spectra display a strong absorption feature at $\sim 5400$ \AA\ which is not seen in other SNe~Ic-3/4 at this epoch. This feature has been attributed to He I in the literature. Spectral modelling of the SN in the early photospheric phase suggests the presence of residual He in a C/O dominated shell. However, the behaviour of the He I lines are unusual when compared with He-rich SNe, showing relatively low velocities and weakening rather than strengthening over time. The SN is found to rise to peak $\sim 16$ d after core-collapse reaching a bolometric luminosity of Lp $\sim 3\times10^{42}$ \ergs. Spectral models, including the nebular epoch, show that the SN ejected $2.5-4$ \msun\ of material, with $\sim 1.5$ \msun\ below 5000 \kms, and with a kinetic energy of $(4.5-7)\times10^{51}$ erg. The explosion synthesised $\sim 0.14$ \msun\ of 56Ni. There are significant uncertainties in E(B-V)host and the distance however, which will affect Lp and MNi. SN 2016coi exploded in a host similar to the Large Magellanic Cloud (LMC) and away from star-forming regions. The properties of the SN and the host-galaxy suggest that the progenitor had $M_\mathrm{ZAMS}$ of $23-28$ \msun\ and was stripped almost entirely down to its C/O core at explosion.

Massive foreground galaxy clusters magnify and distort the light of objects behind them, permitting a view into both the extremely distant and intrinsically faint galaxy populations. We present here the z ~ 6 - 8 candidate high-redshift galaxies from the Reionization Lensing Cluster Survey (RELICS), a Hubble and Spitzer Space Telescope survey of 41 massive galaxy clusters spanning an area of ~200 arcmin^2. These clusters were selected to be excellent lenses and we find similar high-redshift sample sizes and magnitude distributions as CLASH. We discover 321 candidate galaxies with photometric redshifts between z ~ 6 to z ~ 8, including extremely bright objects with H-band magnitudes of m_AB ~ 23 mag. As a sample, the observed (lensed) magnitudes of these galaxies are among the brightest known at z> 6, comparable to much wider, blank-field surveys. RELICS demonstrates the efficiency of using strong gravitational lenses to produce high-redshift samples in the epoch of reionization. These brightly observed galaxies are excellent targets for follow-up study with current and future observatories, including the James Webb Space Telescope.

Strong gravitational lensing by galaxy clusters magnifies background galaxies, enhancing our ability to discover statistically significant samples of galaxies at z>6, in order to constrain the high-redshift galaxy luminosity functions. Here, we present the first five lens models out of the Reionization Lensing Cluster Survey (RELICS) Hubble Treasury Program, based on new HST WFC3/IR and ACS imaging of the clusters RXC J0142.9+4438, Abell 2537, Abell 2163, RXC J2211.7-0349, and ACT-CLJ0102-49151. The derived lensing magnification is essential for estimating the intrinsic properties of high-redshift galaxy candidates, and properly accounting for the survey volume. We report on new spectroscopic redshifts of multiply imaged lensed galaxies behind these clusters, which are used as constraints, and detail our strategy to reduce systematic uncertainties due to lack of spectroscopic information. In addition, we quantify the uncertainty on the lensing magnification due to statistical and systematic errors related to the lens modeling process, and find that in all but one cluster, the magnification is constrained to better than 20% in at least 80% of the field of view, including statistical and systematic uncertainties. The five clusters presented in this paper span the range of masses and redshifts of the clusters in the RELICS program. We find that they exhibit similar strong lensing efficiencies to the clusters targeted by the Hubble Frontier Fields within the WFC3/IR field of view. Outputs of the lens models are made available to the community through the Mikulski Archive for Space Telescopes

We demonstrate that deep good-seeing VLT/HAWK-I $K_\mathrm{s}$ images complemented with $g$+$z$-band photometry can yield a sensitivity for weak lensing studies of massive galaxy clusters at redshifts \mbox{$0.7\lesssim z \lesssim 1.1$}, which is almost identical to the sensitivity of HST/ACS mosaics of single-orbit depth. Key reasons for this good performance are the excellent image quality frequently achievable for $K_\mathrm{s}$ imaging from the ground, a highly effective photometric selection of background galaxies, and a galaxy ellipticity dispersion that is noticeably lower than for optically observed high-redshift galaxy samples. Incorporating results from the 3D-HST and UltraVISTA surveys we also obtained a more accurate calibration of the source redshift distribution than previously achieved for similar optical weak lensing data sets. Here we studied the extremely massive galaxy cluster RCS2$J$232727.7$-$020437 (\mbox{$z=0.699$}), combining deep VLT/\mbox{HAWK-I} $K_\mathrm{s}$ images (point spread function with a 0\farcs35 full width at half maximum) with LBT/LBC photometry. The resulting weak lensing mass reconstruction suggests that the cluster consists of a single overdensity, which is detected with a peak significance of $10.1\sigma$. We constrained the cluster mass to \mbox{$M_\mathrm{200c}/(10^{15} \mathrm{M}_\odot) =2.06^{+0.28}_{-0.26}(\mathrm{stat.})\pm 0.12 (\mathrm{sys.})$} assuming a spherical Navarro, Frenk \& White model and simulation-based priors on the concentration, making it one of the most massive galaxy clusters known in the \mbox{$z\gtrsim 0.7$} Universe. We also cross-checked the HAWK-I measurements through an analysis of overlapping HST/ACS images, yielding fully consistent estimates of the lensing signal.

Galaxy-cluster gravitational lenses can magnify background galaxies by a total factor of up to ~50. Here we report an image of an individual star at redshift z=1.49 (dubbed "MACS J1149 Lensed Star 1 (LS1)") magnified by >2000. A separate image, detected briefly 0.26 arcseconds from LS1, is likely a counterimage of the first star demagnified for multiple years by a >~3 solar-mass object in the cluster. For reasonable assumptions about the lensing system, microlensing fluctuations in the stars' light curves can yield evidence about the mass function of intracluster stars and compact objects, including binary fractions and specific stellar evolution and supernova models. Dark-matter subhalos or massive compact objects may help to account for the two images' long-term brightness ratio.

Multi-band photometric and multi-object spectroscopic surveys of merging galaxy clusters allow for the characterization of the distributions of constituent dark matter and galaxy populations, constraints on the dynamics of the merging subclusters, and an understanding of galaxy evolution of member galaxies. We present deep photometric observations from Subaru/SuprimeCam and a catalog of $\sim$5400 spectroscopic cluster members from Keck/DEIMOS across 29 merging galaxy clusters ranging in redshift from $z=0.07$ to $0.55$. The ensemble is compiled based on the presence of radio relics, which highlight cluster scale collisionless shocks in the intra-cluster medium. Together with the spectroscopic and photometric information, the velocities, timescales, and geometries of the respective merging events may be tightly constrained. In this preliminary analysis, the velocity distributions of 28 of the 29 clusters are shown to be well fit by single Gaussians. This indicates that radio relic mergers largely occur transverse to the line of sight and/or near apocenter. In this paper, we present our optical and spectroscopic surveys, preliminary results, and a discussion of the value of radio relic mergers for developing accurate dynamical models of each system.

A massive galaxy cluster can serve as a magnifying glass for distant stellar populations, with strong gravitational lensing exposing details in the lensed background galaxies that would otherwise be undetectable. The MACS J0416.1-2403 cluster (hereafter MACS0416) is one of the most efficient lenses in the sky, and in 2014 it was observed with high-cadence imaging from the Hubble Space Telescope (HST). Here we describe two unusual transient events that appeared behind MACS0416 in a strongly lensed galaxy at redshift $z = 1.0054 \pm 0.0002$. These transients---designated HFF14Spo-NW and HFF14Spo-SE and collectively nicknamed "Spock"---were faster and fainter than any supernova (SN), but significantly more luminous than a classical nova. They reached peak luminosities of $\sim10^{41}$ erg s$^{-1}$ ($M_{\rm AB} < -14$ mag) in 5 rest-frame days, then faded below detectability in roughly the same time span. Models of the cluster lens suggest that these events may be spatially coincident at the source plane, but are most likely not temporally coincident. We find that HFF14Spo can be explained as a luminous blue variable (LBV), a recurrent nova (RN), or a pair of stellar microlensing events. To distinguish between these hypotheses will require a clarification of the positions of nearby critical curves, along with high-cadence monitoring of the field that could detect new transient episodes in the host galaxy.

The CIII] and CIV rest-frame UV emission lines are powerful probes of the ionizations states of galaxies. They have furthermore been suggested as alternatives for spectroscopic redshift confirmation of objects at the epoch of reionization ($z>6$), where the most frequently used redshift indicator, Ly$\alpha$, is attenuated by the high fraction of neutral hydrogen in the inter-galactic medium. However, currently only very few confirmations of carbon UV lines at these high redshifts exist, making it challenging to quantify these claims. Here, we present the detection of CIV$\lambda\lambda$1548,1551\AA\ in \HST\ slitless grism spectroscopy obtained by GLASS of a Ly$\alpha$ emitter at $z=6.11$ multiply imaged by the massive foreground galaxy cluster RXJ2248. The CIV emission is detected at the 3--5$\sigma$ level in two images of the source, with marginal detection in two other images. We do not detect significant CIII]$\lambda\lambda$1907,1909\AA\ emission implying an equivalent width EW$_\textrm{CIII]}<20$\AA\ (1$\sigma$) and $\textrm{CIV/CIII}>0.7$ (2$\sigma$). Combined with limits on the rest-frame UV flux from the HeII$\lambda$1640\AA\ emission line and the OIII]$\lambda\lambda$1661,1666\AA\ doublet, we put constraints on the metallicity and the ionization state of the galaxy. The estimated line ratios and equivalent widths do not support a scenario where an AGN is responsible for ionizing the carbon atoms. SED fits including nebular emission lines imply a source with a mass of log(M/M$_\odot)\sim9$, SFR of around 10M$_\odot$/yr, and a young stellar population $<50$Myr old. The source shows a stronger ionizing radiation field than objects with detected CIV emission at $z<2$ and adds to the growing sample of low-mass (log(M/M$_\odot)\lesssim9$) galaxies at the epoch of reionization with strong radiation fields from star formation.

We present the first results of the KMOS Lens-Amplified Spectroscopic Survey (KLASS), a new ESO Very Large Telescope (VLT) large program, doing multi-object integral field spectroscopy of galaxies gravitationally lensed behind seven galaxy clusters selected from the HST Grism Lens-Amplified Survey from Space (GLASS). Using the power of the cluster magnification we are able to reveal the kinematic structure of 25 galaxies at $0.7 \lesssim z \lesssim 2.3$, in four cluster fields, with stellar masses $8 \lesssim \log{(M_\star/M_\odot)} \lesssim 11$. This sample includes 5 sources at $z>1$ with lower stellar masses than in any previous kinematic IFU surveys. Our sample displays a diversity in kinematic structure over this mass and redshift range. The majority of our kinematically resolved sample is rotationally supported, but with a lower ratio of rotational velocity to velocity dispersion than in the local universe, indicating the fraction of dynamically hot disks changes with cosmic time. We find no galaxies with stellar mass $<3 \times 10^9 M_\odot$ in our sample display regular ordered rotation. Using the enhanced spatial resolution from lensing, we resolve a lower number of dispersion dominated systems compared to field surveys, competitive with findings from surveys using adaptive optics. We find that the KMOS IFUs recover emission line flux from HST grism-selected objects more faithfully than slit spectrographs. With artificial slits we estimate slit spectrographs miss on average 60% of the total flux of emission lines, which decreases rapidly if the emission line is spatially offset from the continuum.

Using deep Hubble Frontier Fields imaging and slitless spectroscopy from the Grism Lens-Amplified Survey from Space, we analyze 2200 cluster and 1748 field galaxies at $0.2\leq z\leq0.7$ to determine the impact of environment on galaxy size and structure at $\log M_*/M_\odot>7.8$, an unprecedented limit at these redshifts. Based on simple assumptions-$r_e=f(M_*)$-we find no significant differences in half-light radii ($r_e$) between equal-mass cluster or field systems. More complex analyses-$r_e=f(M_*,U-V,n,z,\Sigma$)-reveal local density $(\Sigma$) to induce only a $7\% \pm 3\%$ ($95\%$ confidence) reduction in $r_e$ beyond what can be accounted for by $U-V$ color, Sersic index ($n$), and redshift ($z$) effects.Almost any size difference between galaxies in high- and low-density regions is thus attributable to their different distributions in properties other than environment. Indeed, we find a clear color-$r_e$ correlation in low-mass passive cluster galaxies ($\log M_*/M_\odot<9.8$) such that bluer systems have larger radii, with the bluest having sizes consistent with equal-mass star-forming galaxies. We take this as evidence that large-$r_e$ low-mass passive cluster galaxies are recently acquired systems that have been environmentally quenched without significant structural transformation (e.g., by ram pressure stripping or starvation).Conversely, $\sim20\%$ of small-$r_e$ low-mass passive cluster galaxies appear to have been in place since $z\sim3$. Given the consistency of the small-$r_e$ galaxies' stellar surface densities (and even colors) with those of systems more than ten times as massive, our findings suggest that clusters mark places where galaxy evolution is accelerated for an ancient base population spanning most masses, with late-time additions quenched by environment-specific mechanisms are mainly restricted to the lowest masses.

Gravitational lensing by clusters of galaxies offers a powerful probe of their structure and mass distribution. Deriving a lens magnification map for a galaxy cluster is a classic inversion problem and many methods have been developed over the past two decades to solve it. Several research groups have developed techniques independently to map the predominantly dark matter distribution in cluster lenses. While these methods have all provided remarkably high precision mass maps, particularly with exquisite imaging data from the Hubble Space Telescope (HST), the reconstructions themselves have never been directly compared. In this paper, we report the results of comparing various independent lens modeling techniques employed by individual research groups in the community. Here we present for the first time a detailed and robust comparison of methodologies for fidelity, accuracy and precision. For this collaborative exercise, the lens modeling community was provided simulated cluster images -- of two clusters Ares and Hera -- that mimic the depth and resolution of the ongoing HST Frontier Fields. The results of the submitted reconstructions with the un-blinded true mass profile of these two clusters are presented here. Parametric, free-form and hybrid techniques have been deployed by the participating groups and we detail the strengths and trade-offs in accuracy and systematics that arise for each methodology. We note in conclusion that lensing reconstruction methods produce reliable mass distributions that enable the use of clusters as extremely valuable astrophysical laboratories and cosmological probes.

In Hubble Space Telescope (HST) imaging taken on 10 November 2014, four images of supernova (SN) "Refsdal" (redshift z=1.49) appeared in an Einstein-cross--like configuration (images S1-S4) around an early-type galaxy in the cluster MACS J1149.5+2223 (z=0.54). Almost all lens models of the cluster have predicted that the SN should reappear within a year in a second host-galaxy image created by the cluster's potential. In HST observations taken on 11 December 2015, we find a new source at the predicted position of the new image of SN Refsdal approximately 8" from the previous images S1-S4. This marks the first time the appearance of a SN at a particular time and location in the sky was successfully predicted in advance! We use these data and the light curve from the first four observed images of SN Refsdal to place constraints on the relative time delay and magnification of the new image (SX), compared to images S1-S4. This enables us, for the first time, to test "blind" lens model predictions of both magnifications and time delays for a lensed SN. We find that the timing and brightness of the new image are consistent with the blind predictions of a fraction of the models. The reappearance illustrates the discriminatory power of this blind test and its utility to uncover sources of systematic uncertainty. From planned HST photometry, we expect to reach a precision of 1-2% on the time delay between S1-S4 and SX.

The galaxy cluster RX J0603.3+4214 at z=0.225 is one of the rarest clusters boasting an extremely large (~2 Mpc) radio-relic. Because of the remarkable morphology of the relic, the cluster is nicknamed "Toothbrush Cluster". Although the cluster's underlying mass distribution is one of the critical pieces of information needed to reconstruct the merger scenario responsible for the puzzling radio-relic morphology, its proximity to the Galactic plane b~10 deg has imposed significant observational challenges. We present a high-resolution weak-lensing study of the cluster with Subaru/Suprime Cam and Hubble Space Telescope imaging data. Our mass reconstruction reveals that the cluster is comprised of complicated dark matter substructures closely tracing the galaxy distribution, however in contrast with the relatively simple binary X-ray morphology. Nevertheless, we find that the cluster mass is still dominated by the two most massive clumps aligned north-south with a ~3:1 mass ratio (M_{200}=6.29_{-1.62}^{+2.24} x 10^{14} Msun and 1.98_{-0.74}^{+1.24} x 10^{14} Msun for the northern and southern clumps, respectively). The southern mass peak is ~2' offset toward the south with respect to the corresponding X-ray peak, which has a "bullet"-like morphology pointing south. Comparison of the current weak-lensing result with the X-ray, galaxy, and radio-relic suggests that perhaps the dominant mechanism responsible for the observed relic may be a high-speed collision of the two most massive subclusters, although the peculiarity of the morphology necessitates involvement of additional sub-clusters. Careful numerical simulations should follow in order to obtain more complete understanding of the merger scenario utilizing all existing observations.

We present the first year of Hubble Space Telescope imaging of the unique supernova (SN) 'Refsdal', a gravitationally lensed SN at z=1.488$\pm$0.001 with multiple images behind the galaxy cluster MACS J1149.6+2223. The first four observed images of SN Refsdal (images S1-S4) exhibited a slow rise (over ~150 days) to reach a broad peak brightness around 20 April, 2015. Using a set of light curve templates constructed from SN 1987A-like peculiar Type II SNe, we measure time delays for the four images relative to S1 of 4$\pm$4 (for S2), 2$\pm$5 (S3), and 24$\pm$7 days (S4). The measured magnification ratios relative to S1 are 1.15$\pm$0.05 (S2), 1.01$\pm$0.04 (S3), and 0.34$\pm$0.02 (S4). None of the template light curves fully captures the photometric behavior of SN Refsdal, so we also derive complementary measurements for these parameters using polynomials to represent the intrinsic light curve shape. These more flexible fits deliver fully consistent time delays of 7$\pm$2 (S2), 0.6$\pm$3 (S3), and 27$\pm$8 days (S4). The lensing magnification ratios are similarly consistent, measured as 1.17$\pm$0.02 (S2), 1.00$\pm$0.01 (S3), and 0.38$\pm$0.02 (S4). We compare these measurements against published predictions from lens models, and find that the majority of model predictions are in very good agreement with our measurements. Finally, we discuss avenues for future improvement of time delay measurements -- both for SN Refsdal and for other strongly lensed SNe yet to come.

We have acquired Hubble Space Telescope (HST) and Very Large Telescope near-infrared spectra and images of supernova (SN) Refsdal after its discovery as an Einstein cross in Fall 2014. The HST light curve of SN Refsdal matches the distinctive, slowly rising light curves of SN 1987A-like supernovae (SNe), and we find strong evidence for a broad H-alpha P-Cygni profile in the HST grism spectrum at the redshift (z = 1.49) of the spiral host galaxy. SNe IIn, powered by circumstellar interaction, could provide a good match to the light curve of SN Refsdal, but the spectrum of a SN IIn would not show broad and strong H-alpha absorption. From the grism spectrum, we measure an H-alpha expansion velocity consistent with those of SN 1987A-like SNe at a similar phase. The luminosity, evolution, and Gaussian profile of the H-alpha emission of the WFC3 and X-shooter spectra, separated by ~2.5 months in the rest frame, provide additional evidence that supports the SN 1987A-like classification. In comparison with other examples of SN 1987A-like SNe, SN Refsdal has a blue B-V color and a high luminosity for the assumed range of potential magnifications. If SN Refsdal can be modeled as a scaled version of SN 1987A, we estimate it would have an ejecta mass of 20+-5 solar masses. The evolution of the light curve at late times will provide additional evidence about the potential existence of any substantial circumstellar material (CSM). Using MOSFIRE and X-shooter spectra, we estimate a subsolar host-galaxy metallicity (8.3+-0.1 dex and <8.4 dex, respectively) near the explosion site.

Supernova "Refsdal," multiply imaged by cluster MACSJ1149.5+2223, represents a rare opportunity to make a true blind test of model predictions in extragalactic astronomy, on a time scale that is short compared to a human lifetime. In order to take advantage of this event, we produced seven gravitational lens models with five independent methods, based on Hubble Space Telescope (HST) Hubble Frontier Field images, along with extensive spectroscopic follow-up observations by HST, the Very Large and the Keck Telescopes. We compare the model predictions and show that they agree reasonably well with the measured time delays and magnification ratios between the known images, even though these quantities were not used as input. This agreement is encouraging, considering that the models only provide statistical uncertainties, and do not include additional sources of uncertainties such as structure along the line of sight, cosmology, and the mass sheet degeneracy. We then present the model predictions for the other appearances of SN "Refsdal." A future image will reach its peak in the first half of 2016, while another image appeared between 1994 and 2004. The past image would have been too faint to be detected in existing archival images. The future image should be approximately one third as bright as the brightest known image (i.e., H_AB~25.7 mag at peak and H_AB~26.7 mag six months before peak), and thus detectable in single-orbit HST images. We will find out soon whether our predictions are correct.

We present the first study of the spatial distribution of star formation in z~0.5 cluster galaxies. The analysis is based on data taken with the Wide Field Camera 3 as part of the Grism Lens-Amplified Survey from Space (GLASS). We illustrate the methodology by focusing on two clusters (MACS0717.5+3745 and MACS1423.8+2404) with different morphologies (one relaxed and one merging) and use foreground and background galaxies as field control sample. The cluster+field sample consists of 42 galaxies with stellar masses in the range 10^8-10^11 M_sun, and star formation rates in the range 1-20 M_sun/yr. Both in clusters and in the field, H{\alpha} is more extended than the rest-frame UV continuum in 60% of the cases, consistent with diffuse star formation and inside out growth. In ~20% of the cases, the H{\alpha} emission appears more extended in cluster galaxies than in the field, pointing perhaps to ionized gas being stripped and/or star formation being enhanced at large radii. The peak of the H{\alpha} emission and that of the continuum are offset by less than 1 kpc. We investigate trends with the hot gas density as traced by the X-ray emission, and with the surface mass density as inferred from gravitational lens models and find no conclusive results. The diversity of morphologies and sizes observed in H_alpha illustrates the complexity of the environmental process that regulate star formation. Upcoming analysis of the full GLASS dataset will increase our sample size by almost an order of magnitude, verifying and strengthening the inference from this initial dataset.

SN HFF14Tom is a Type Ia Supernova (SN) discovered at z = 1.3457 +- 0.0001 behind the galaxy cluster Abell 2744 (z = 0.308). In a cosmology-independent analysis, we find that HFF14Tom is 0.77 +- 0.15 magnitudes brighter than unlensed Type Ia SNe at similar redshift, implying a lensing magnification of mu_obs = 2.03 +- 0.29. This observed magnification provides a rare opportunity for a direct empirical test of galaxy cluster lens models. Here we test 17 lens models, 13 of which were generated before the SN magnification was known, qualifying as pure "blind tests". The models are collectively fairly accurate: 8 of the models deliver median magnifications that are consistent with the measured mu to within 1-sigma. However, there is a subtle systematic bias: the significant disagreements all involve models overpredicting the magnification. We evaluate possible causes for this mild bias, and find no single physical or methodological explanation to account for it. We do find that model accuracy can be improved to some extent with stringent quality cuts on multiply-imaged systems, such as requiring that a large fraction have spectroscopic redshifts. In addition to testing model accuracies as we have done here, Type Ia SN magnifications could also be used as inputs for future lens models of Abell 2744 and other clusters, providing valuable constraints in regions where traditional strong- and weak-lensing information is unavailable.

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

We present a reconstructions of galaxy-cluster-scale mass distributions from simulated gravitational lensing data sets including strong lensing, weak lensing shear, and measurements of quadratic image distortions -- flexion. The lensing data is constructed to make a direct comparison between mass reconstructions with and without flexion. We show that in the absence of flexion measurements, significant galaxy-group scale substructure can remain undetected in the reconstructed mass profiles, and that the resulting profiles underestimate the aperture mass in the substructure regions by $\sim25-40\%$. When flexion is included, subhaloes down to a mass of $\sim3\times10^{12}$ M$_\odot$ can be detected at an angular resolution smaller than 10\arcsec. Aperture masses from profiles reconstructed with flexion match the input distribution values to within an error of $\sim13\%$, including both statistical error and scatter. This demonstrates the important constraint that flexion measurements place on substructure in galaxy clusters and its utility for producing high-fidelity mass reconstructions.

In 1964, Refsdal hypothesized that a supernova whose light traversed multiple paths around a strong gravitational lens could be used to measure the rate of cosmic expansion. We report the discovery of such a system. In Hubble Space Telescope imaging, we have found four images of a single supernova forming an Einstein cross configuration around a redshift z=0.54 elliptical galaxy in the MACS J1149.6+2223 cluster. The cluster's gravitational potential also creates multiple images of the z=1.49 spiral supernova host galaxy, and a future appearance of the supernova elsewhere in the cluster field is expected. The magnifications and staggered arrivals of the supernova images probe the cosmic expansion rate, as well as the distribution of matter in the galaxy and cluster lenses.

We present spatially resolved gas-phase metallicity for a system of three galaxies at z=1.85 detected in the Grism Lens-Amplified Survey from Space (GLASS). The combination of HST's diffraction limit and strong gravitational lensing by the cluster MACS J0717+3745 results in a spatial resolution of ~200-300 pc, enabling good spatial sampling despite the intrinsically small galaxy sizes. The galaxies in this system are separated by 50-200 kpc in projection and are likely in an early stage of interaction, evidenced by relatively high specific star formation rates. Their gas-phase metallicities are consistent with larger samples at similar redshift, star formation rate, and stellar mass. We obtain a precise measurement of the metallicity gradient for one galaxy and find a shallow slope compared to isolated galaxies at high redshift, consistent with a flattening of the gradient due to gravitational interaction. An alternative explanation for the shallow metallicity gradient and elevated star formation rate is rapid recycling of metal-enriched gas, but we find no evidence for enhanced gas-phase metallicities which should result from this effect. Notably, the measured stellar masses log(M/Msun) = 7.2-9.1 probe to an order of magnitude below previous mass-metallicity studies at this redshift. The lowest mass galaxy has properties similar to those expected for Fornax at this redshift, indicating that GLASS is able to directly study the progenitors of local group dwarf galaxies on spatially resolved scales. Larger samples from the full GLASS survey will be ideal for studying the effects of feedback, and the time evolution of metallicity gradients. These initial results demonstrate the utility of HST spectroscopy combined with gravitational lensing for characterizing resolved physical properties of galaxies at high redshift.

We present a pilot study on the origin and assembly history of the ICL for four galaxy clusters at 0.44<z<0.57 observed with the Hubble Space Telescope from the Cluster Lensing and Supernova Survey with Hubble (CLASH) sample. Using this sample of clusters we set an empirical limit on the amount of scatter in ICL surface brightness profiles of such clusters at z=0.5 and constrain the progenitor population and formation mechanism of the ICL by measuring the ICL surface brightness profile, the ICL color and color gradient, and the total ICL luminosity within 10<r<110 kpc. The observed scatter is physical, which we associate with differences in ICL assembly process, formation epoch, and/or ICL content. Using stellar population synthesis models we transform the observed colors to metallicity. For three of the four clusters we find clear negative gradients that, on average, decrease from super solar in the central regions of the BCG to sub-solar in the ICL. Such negative color/metallicity gradients can arise from tidal stripping of L* galaxies and/or the disruption of dwarf galaxies, but not major mergers with the BCG. We also find that the ICL at 110 kpc has a color comparable to m*+2 red sequence galaxies and a total luminosity between 10<r<110 kpc of 4-8 L*. This suggests that the ICL is dominated by stars liberated from galaxies with L>0.2 L* and that neither dwarf disruption nor major mergers with the BCG alone can explain the observed level of luminosity and remain consistent with either the observed evolution in the faint end slope of the luminosity function or predictions for the number of BCG major mergers since z=1. Taken together, the results of this pilot study are suggestive of a formation history for these clusters in which the ICL is built-up by the stripping of >0.2 L* galaxies, and disfavor significant contribution to the ICL by dwarf disruption or major mergers with the BCG.

We present Spitzer/IRAC observations of nine $z'$-band dropouts highly magnified (2<mu<12) by the Bullet Cluster. We combine archival imaging with our Exploratory program (SURFS UP), which results in a total integration time of ~30 hr per IRAC band. We detect (>3sigma) in both IRAC bands the brightest of these high-redshift galaxies, with [3.6]=23.80+-0.28 mag, [4.5]=23.78+-0.25 mag, and (H-[3.6])=1.17+-0.32 mag. The remaining eight galaxies are undetected to [3.6]~26.4 mag and [4.5]~26.0 mag with stellar masses of ~5x10^7 M_sol. The detected galaxy has an estimated magnification of mu=12+-4, which implies this galaxy has an ultraviolet luminosity of L_1500~0.3 L*_{z=7} --- the lowest luminosity individual source detected in IRAC at z>7. By modeling the broadband photometry, we estimate the galaxy has an intrinsic star-formation rate of SFR~1.3 M_sol/yr and stellar mass of M~2x10^9 M_sol, which gives a specific star-formation rate of sSFR~0.7 Gyr^-1. If this galaxy had sustained this star-formation rate since z~20, it could have formed the observed stellar mass (to within a factor of ~2), we also discuss alternate star-formation histories and argue the exponentially-increasing model is unlikely. Finally, based on the intrinsic star-formation rate, we estimate this galaxy has a likely [C II] flux of <f_[C II]> = 10^{-17} erg/s/cm2.

The Grism Lens-Amplified Survey from Space (GLASS) is a Hubble Space Telescope (HST) Large Program, which will obtain 140 orbits of grism spectroscopy of the core and infall regions of 10 galaxy clusters, selected to be among the very best cosmic telescopes. Extensive HST imaging is available from many sources including the CLASH and Frontier Fields programs. We introduce the survey by analyzing spectra of faint multiply-imaged galaxies and $z\gtrsim6$ galaxy candidates obtained from the first seven orbits out of fourteen targeting the core of the Frontier Fields cluster MACS0717.5+3745. Using the G102 and G141 grisms to cover the wavelength range 0.8-1.7$\mu$m, we confirm 4 strongly lensed systems by detecting emission lines in each of the images. For the 9 $z\gtrsim6$ galaxy candidates clear from contamination, we do not detect any emission lines down to a seven-orbit 1$\sigma$ noise level of $\sim$5$\times$10$^{-18}$erg s$^{-1}$cm$^{-2}$. Taking lensing magnification into account, our flux sensitivity reaches $\sim$0.2-5$\times$10$^{-18}$erg s$^{-1}$cm$^{-2}$. These limits over an uninterrupted wavelength range rule out the possibility that the high-$z$ galaxy candidates are instead strong line emitters at lower redshift. These results show that by means of careful modeling of the background - and with the assistance of lensing magnification - interesting flux limits can be reached for large numbers of objects, avoiding pre-selection and the wavelength restrictions inherent to ground-based multi-slit spectroscopy. These observations confirm the power of slitless HST spectroscopy even in fields as crowded as a cluster core.

Strong lensing gravitational time delays are a powerful and cost effective probe of dark energy. Recent studies have shown that a single lens can provide a distance measurement with 6-7 % accuracy (including random and systematic uncertainties), provided sufficient data are available to determine the time delay and reconstruct the gravitational potential of the deflector. Gravitational-time delays are a low redshift (z~0-2) probe and thus allow one to break degeneracies in the interpretation of data from higher-redshift probes like the cosmic microwave background in terms of the dark energy equation of state. Current studies are limited by the size of the sample of known lensed quasars, but this situation is about to change. Even in this decade, wide field imaging surveys are likely to discover thousands of lensed quasars, enabling the targeted study of ~100 of these systems and resulting in substantial gains in the dark energy figure of merit. In the next decade, a further order of magnitude improvement will be possible with the 10000 systems expected to be detected and measured with LSST and Euclid. To fully exploit these gains, we identify three priorities. First, support for the development of software required for the analysis of the data. Second, in this decade, small robotic telescopes (1-4m in diameter) dedicated to monitoring of lensed quasars will transform the field by delivering accurate time delays for ~100 systems. Third, in the 2020's, LSST will deliver 1000's of time delays; the bottleneck will instead be the aquisition and analysis of high resolution imaging follow-up. Thus, the top priority for the next decade is to support fast high resolution imaging capabilities, such as those enabled by the James Webb Space Telescope and next generation adaptive optics systems on large ground based telescopes.

In order to use strong gravitational lens time delays to measure precise and accurate cosmological parameters the effects of mass along the line of sight must be taken into account. We present a method to achieve this by constraining the probability distribution function of the effective line of sight convergence k_ext. The method is based on matching the observed overdensity in the weighted number of galaxies to that found in mock catalogs with k_ext obtained by ray-tracing through structure formation simulations. We explore weighting schemes based on projected distance, mass, luminosity, and redshift. This additional information reduces the uncertainty of k_ext from sigma_k $0.06 to ~0.04 for very overdense lines of sight like that of the system B1608+656. For more common lines of sight, sigma_k is reduced to ~<0.03, corresponding to an uncertainty of ~<3% on distance. This uncertainty has comparable effects on cosmological parameters to that arising from the mass model of the deflector and its immediate environment. Photometric redshifts based on g, r, i and K photometries are sufficient to constrain k_ext almost as well as with spectroscopic redshifts. As an illustration, we apply our method to the system B1608+656. Our most reliable k_ext estimator gives sigma_k=0.047 down from 0.065 using only galaxy counts. Although deeper multi-band observations of the field of B1608+656 are necessary to obtain a more precise estimate, we conclude that griK photometry, in addition to spectroscopy to characterize the immediate environment, is an effective way to increase the precision of time-delay cosmography.

Merging clusters of galaxies are unique in their power to directly probe and place limits on the self-interaction cross-section of dark matter. Detailed observations of several merging clusters have shown the intracluster gas to be displaced from the centroids of dark matter and galaxy density by ram pressure, while the latter components are spatially coincident, consistent with collisionless dark matter. This has been used to place upper limits on the dark matter particle self-inteaction cross-section of order 1 cm^2/g. The cluster Abell 520 has been seen as a possible exception. We revisit A520 presenting new HST ACS mosaic images and a Magellan image set. We perform a detailed weak lensing analysis and show that the weak lensing mass measurements and morphologies of the core galaxy-filled structures are mostly in good agreement with previous works. There is however one significant difference -- we do not detect the previously claimed "dark core" that contains excess mass with no significant galaxy overdensity at the location of the X-ray plasma. This peak has been suggested to be indicative of a large self-interaction cross-section for dark matter (at least ~5 sigma larger than the upper limit of 0.7 cm^2/g determined by observations of the Bullet Cluster). We find no such indication and instead find that the mass distribution of A520, after subtraction of the X-ray plasma mass, is in good agreement with the luminosity distribution of the cluster galaxies. We conclude that A520 shows no evidence to contradict the collisionless dark matter scenario.

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

We use imaging obtained with the Hubble Space Telescope Wide Field Camera 3 to search for z_850 dropouts at z~7 and J_110 dropouts at z~9 lensed by the Bullet Cluster. In total we find 10 z_850 dropouts in our 8.27 arcmin^2 field. Using magnification maps from a combined weak and strong lensing mass reconstruction of the Bullet Cluster and correcting for estimated completeness levels, we calculate the surface density and luminosity function of our z_850 dropouts as a function of intrinsic (accounting for magnification) magnitude. We find results consistent with published blank field surveys, despite using much shallower data, and demonstrate the effectiveness of cluster surveys in the search for z~7 galaxies.

Here we present the weak lensing results of A1758, which is known to consist of four subclusters undergoing two separate mergers, A1758N and A1758S. Weak lensing results of A1758N agree with previous weak lensing results of clusters 1E0657-558 (Bullet cluster) and MACS J0025.4-1222, whose X-ray gas components were found to be largely separated from their clusters' gravitational potentials. A1758N has a geometry that is different from previously published mergers in that one of its X-ray peaks overlays the corresponding gravitational potential and the other X-ray peak is well separated from its cluster's gravitational potential. The weak lensing mass peaks of the two northern clusters are separated at the 2.5 sigma level. We estimate the combined mass of the clusters in A1758N to be 2.2+/-0.5*10^15 M_sun and r_200=2300(+100/-130) kpc. We also detect seven strong lensing candidates, two of which may provide information that would improve the mass measurements of A1758N.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

In this paper we address the consumption of energy by astronomers while performing their professional duties. Although we find that astronomy uses a negligible fraction of the US energy budget, the rate at which energy is consumed by an average astronomer is similar to that of a typical high-flying businessperson. We review some of the ways in which astronomers are already acting to reduce their energy consumption. In the coming decades, all citizens will have to reduce their energy consumption to conserve fossil fuel reserves and to help avert a potentially catastrophic change in the Earth's climate. The challenges are the same for astronomers as they are for everyone: decreasing the distances we travel and investing in energy-efficient infrastructure. The high profile of astronomy in the media, and the great public interest in our field, can play a role in promoting energy-awareness to the wider population. Our specific recommendations are therefore to 1) reduce travel when possible, through efficient meeting organization, and by investing in high-bandwidth video conference facilities and virtual-world software, 2) create energy-efficient observatories, computing centers and workplaces, powered by sustainable energy resources, and 3) actively publicize these pursuits.

Whereas considerable effort has been afforded in understanding the properties of galaxies, a full physical picture, connecting their baryonic and dark-matter content, super-massive black holes, and (metric) theories of gravity, is still ill-defined. Strong gravitational lensing furnishes a powerful method to probe gravity in the central regions of galaxies. It can (1) provide a unique detection-channel of dark-matter substructure beyond the local galaxy group, (2) constrain dark-matter physics, complementary to direct-detection experiments, as well as metric theories of gravity, (3) probe central super-massive black holes, and (4) provide crucial insight into galaxy formation processes from the dark matter point of view, independently of the nature and state of dark matter. To seriously address the above questions, a considerable increase in the number of strong gravitational-lens systems is required. In the timeframe 2010-2020, a staged approach with radio (e.g. EVLA, e-MERLIN, LOFAR, SKA phase-I) and optical (e.g. LSST and JDEM) instruments can provide 10^(2-4) new lenses, and up to 10^(4-6) new lens systems from SKA/LSST/JDEM all-sky surveys around ~2020. Follow-up imaging of (radio) lenses is necessary with moderate ground/space-based optical-IR telescopes and with 30-50m telescopes for spectroscopy (e.g. TMT, GMT, ELT). To answer these fundamental questions through strong gravitational lensing, a strong investment in large radio and optical-IR facilities is therefore critical in the coming decade. In particular, only large-scale radio lens surveys (e.g. with SKA) provide the large numbers of high-resolution and high-fidelity images of lenses needed for SMBH and flux-ratio anomaly studies.

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

We present evidence for a Spitzer-selected luminous infrared galaxy (LIRG) behind the Bullet Cluster. The galaxy, originally identified in IRAC photometry as a multiply imaged source, has a spectral energy distribution consistent with a highly extincted (A_V~3.3), strongly star-forming galaxy at z=2.7. Using our strong gravitational lensing model presented in Bradac et al. (2006), we find that the magnifications are 10 to 50 for the three images of the galaxy. The implied infrared luminosity is consistent with the galaxy being a LIRG, with a stellar mass of M_*~2e11 M_Sun and a star formation rate of ~90 M_Sun/yr. With lensed fluxes at 24 microns of 0.58 mJy and 0.39 mJy in the two brightest images, this galaxy presents a unique opportunity for detailed study of an obscured starburst with star fomation rate comparable to that of L* galaxies at z>2.

Lens modeling is the key to successful and meaningful automated strong galaxy-scale gravitational lens detection. We have implemented a lens-modeling "robot" that treats every bright red galaxy (BRG) in a large imaging survey as a potential gravitational lens system. Using a simple model optimized for "typical" galaxy-scale lenses, we generate four assessments of model quality that are used in an automated classification. The robot infers the lens classification parameter H that a human would have assigned; the inference is performed using a probability distribution generated from a human-classified training set, including realistic simulated lenses and known false positives drawn from the HST/EGS survey. We compute the expected purity, completeness and rejection rate, and find that these can be optimized for a particular application by changing the prior probability distribution for H, equivalent to defining the robot's "character." Adopting a realistic prior based on the known abundance of lenses, we find that a lens sample may be generated that is ~100% pure, but only ~20% complete. This shortfall is due primarily to the over-simplicity of the lens model. With a more optimistic robot, ~90% completeness can be achieved while rejecting ~90% of the candidate objects. The remaining candidates must be classified by human inspectors. We are able to classify lens candidates by eye at a rate of a few seconds per system, suggesting that a future 1000 square degree imaging survey containing 10^7 BRGs, and some 10^4 lenses, could be successfully, and reproducibly, searched in a modest amount of time. [Abridged]

High resolution HST/ACS images of the galaxy cluster RX J1347-1145 have enabled us to identify several new multiple image candidates in the cluster, including a 5 image system with a central image. The multiple images allow us to construct an accurate 2-dimensional mass map of the central part of the cluster. The modelling of the cluster mass includes the most prominent cluster galaxies modelled as truncated isothermal spheres and a smooth halo component that is described with 2 parametric profiles. The mass reconstruction is done using a Markov chain Monte Carlo method that provides us with a total projected mass density as well as estimates for the parameters of interest and their respective errors. The mass profile is in reasonable agreement with previous mass estimates based on the X-ray emission from the hot intra-cluster gas, however the X-ray mass estimates are systematically lower than what we obtain with gravitational lensing.

(Abridged) We compare recent results from X-ray, strong lensing, weak lensing, and optical observations with numerical simulations of the merging galaxy cluster 1E0657-56. X-ray observations reveal a bullet-like subcluster with a prominent bow shock, while lensing results show that the positions of the total mass peaks are consistent with the centroids of the collisionless galaxies (and inconsistent with the X-ray brightness peaks). Previous studies, based on older observational datasets, have placed upper limits on the self-interaction cross-section of dark matter per unit mass, sigma/m, using simplified analytic techniques. In this work, we take advantage of new, higher-quality observational datasets by running N-body simulations of 1E0657-56 that include the effects of self-interacting dark matter, and comparing the results with observations. Furthermore, the recent data allow for a new independent method of constraining sigma/m, based on the non-observation of an offset between the bullet subcluster mass peak and galaxy centroid. This new method places an upper limit (68% confidence) of sigma/m < 1.25 cm^2/g. If we make the assumption that the subcluster and the main cluster had equal mass-to-light ratios prior to the merger, we derive our most stringent constraint of sigma/m < 0.7 cm^2/g, which comes from the consistency of the subcluster's observed mass-to-light ratio with the main cluster's, and with the universal cluster value, ruling out the possibility of a large fraction of dark matter particles being scattered away due to collisions. Our limit is a slight improvement over the previous result from analytic estimates, and rules out most of the 0.5 - 5cm^2/g range invoked to explain inconsistencies between the standard collisionless cold dark matter model and observations.

We present new weak lensing observations of 1E0657-558 (z=0.296), a unique cluster merger, that enable a direct detection of dark matter, independent of assumptions regarding the nature of the gravitational force law. Due to the collision of two clusters, the dissipationless stellar component and the fluid-like X-ray emitting plasma are spatially segregated. By using both wide-field ground based images and HST/ACS images of the cluster cores, we create gravitational lensing maps which show that the gravitational potential does not trace the plasma distribution, the dominant baryonic mass component, but rather approximately traces the distribution of galaxies. An 8-sigma significance spatial offset of the center of the total mass from the center of the baryonic mass peaks cannot be explained with an alteration of the gravitational force law, and thus proves that the majority of the matter in the system is unseen.

The galaxy cluster 1E0657-56 (z = 0.296) is remarkably well-suited for addressing outstanding issues in both galaxy evolution and fundamental physics. We present a reconstruction of the mass distribution from both strong and weak gravitational lensing data. Multi-color, high-resolution HST ACS images allow detection of many more arc candidates than were previously known, especially around the subcluster. Using the known redshift of one of the multiply imaged systems, we determine the remaining source redshifts using the predictive power of the strong lens model. Combining this information with shape measurements of "weakly" lensed sources, we derive a high-resolution, absolutely-calibrated mass map, using no assumptions regarding the physical properties of the underlying cluster potential. This map provides the best available quantification of the total mass of the central part of the cluster. We also confirm the result from Clowe et al. (2004,2006a).

We have shown that the cluster-mass reconstruction method which combines strong and weak gravitational lensing data, developed in the first paper in the series, successfully reconstructs the mass distribution of a simulated cluster. In this paper we apply the method to the ground-based high-quality multi-colour data of RX J1347.5-1145, the most X-ray luminous cluster to date. A new analysis of the cluster core on very deep, multi-colour data analysis of VLT/FORS data reveals many more arc candidates than previously known for this cluster. The combined strong and weak lensing reconstruction confirms that the cluster is indeed very massive. If the redshift and identification of the multiple-image system as well as the redshift estimates of the source galaxies used for weak lensing are correct, we determine the enclosed cluster mass in a cylinder to M(<360 h^-1 kpc)= (1.2 +/- 0.3) 10^15 Msun. In addition the reconstructed mass distribution follows the distribution found with independent methods (X-ray measurements, SZ). With higher resolution (e.g. HST imaging data) more reliable multiple imaging information can be obtained and the reconstruction can be improved to accuracies greater than what is currently possible with weak and strong lensing techniques.

Weak gravitational lensing is considered to be one of the most powerful tools to study the mass and the mass distribution of galaxy clusters. However, the mass-sheet degeneracy transformation has limited its success. We present a novel method for a cluster mass reconstruction which combines weak and strong lensing information on common scales and can, as a consequence, break the mass-sheet degeneracy. We extend the weak lensing formalism to the inner parts of the cluster and combine it with the constraints from multiple image systems. We demonstrate the feasibility of the method with simulations, finding an excellent agreement between the input and reconstructed mass also on scales within and beyond the Einstein radius. Using a single multiple image system and photometric redshift information of the background sources used for weak and strong lensing analysis, we find that we are effectively able to break the mass-sheet degeneracy, therefore removing one of the main limitations on cluster mass estimates. We conclude that with high resolution (e.g. HST) imaging data the method can more accurately reconstruct cluster masses and their profiles than currently existing lensing techniques.

Weak gravitational lensing is considered to be one of the most powerful tools to study the mass and the mass distribution of galaxy clusters. However, weak lensing mass reconstructions are plagued by the so-called mass-sheet degeneracy--the surface mass density \kappa of the cluster can be determined only up to a degeneracy transformation \kappa \to \kappa' = \lambda \kappa + (1 -\lambda), where \lambda is an arbitrary constant. This transformation fundamentally limits the accuracy of cluster mass determinations if no further assumptions are made. We discuss here a possibility to break the mass-sheet degeneracy in weak lensing mass maps using distortion and redshift information of background galaxies. Compared to other techniques proposed in the past, it does not rely on any assumptions on cluster potential and does not make use of weakly constrained information (such as the source number counts, used in the magnification effect).Our simulations show that we are effectively able to break the mass-sheet degeneracy for supercritical lenses and that for undercritical lenses the mass-sheet degeneracy is very difficult to be broken, even under idealised conditions.

Weak gravitational lensing is considered to be one of the most powerful tools to study the mass and the mass distribution of galaxy clusters. However, weak lensing mass reconstructions are plagued by the so-called mass-sheet degeneracy--the surface mass density \kappa of the cluster can be determined only up to a degeneracy transformation \kappa \to \kappa' = \lambda \kappa + (1 -\lambda), where \lambda is an arbitrary constant. This transformation fundamentally limits the accuracy of cluster mass determinations if no further assumptions are made. We describe here a method to break the mass-sheet degeneracy in weak lensing mass maps using distortion and redshift information of background galaxies and illustrate this by two simple toy models. Compared to other techniques proposed in the past, it does not rely on any assumptions on cluster potential; it can be easily applied to non-parametric mass-reconstructions and no assumptions on boundary conditions have to be made. In addition it does not make use of weakly constrained information (such as the source number counts, used in the magnification effect). Our simulations show that we are effectively able to break the mass-sheet degeneracy for supercritical lenses, but that for undercritical lenses the mass-sheet degeneracy is very difficult to be broken, even under idealised conditions.

We present a study of the lens properties of quadruply imaged systems, lensed by numerically simulated galaxies. We investigate a simulated elliptical and disc galaxy drawn from high resolution simulations of galaxy formation in a concordance $\Lambda$CDM universe. The simulations include the effects of gas dynamics, star formation and feedback processes. Flux-ratio anomalies observed in strong gravitational lensing potentially provide an indicator for the presence of mass substructure in lens galaxies as predicted from CDM simulations. We particularly concentrate on the prediction that, for an ideal cusp caustic, the sum of the signed magnifications of the three highly magnified images should vanish when the source approaches the cusp. Strong violation of this cusp relation indicates the presence of substructure, regardless of the global, smooth mass model of the lens galaxy. We draw the following conclusions: (1) the level of substructure present in simulations produces violations of the cusp relation comparable to those observed, (2) higher-order catastrophes (e.g. swallowtails) can also cause changes of the order of 0.6 in the cusp relation as predicted by a smooth model, (3) the flux anomaly distribution depends on the image parity and flux and both the brightest minimum and saddle-point images are more affected by substructure than the fainter images. In addition, the brightest saddle point is demagnified w.r.t. the brightest minimum. Our results are fully numerical and properly include all mass scales, without making semi-analytic assumptions. They are ultimately limited by the mass resolution of single particles in the simulation determined by current computational limits, however show that our results are not affected by shot-noise due to the finite number of particles.

Gravitationally lensed systems with multiply imaged quasars are an excellent tool for studying the properties of distant galaxies. In particular, they provide the most accurate mass measures for the lensing galaxy. The system B1422+231 is a well studied example of a quadruply imaged quasar, with high-quality VLBI data available. Very accurate data on image positions, fluxes and deconvolved image sizes provide good constraints for lensing models. We discuss here the failure of smooth models in fitting the data. Since it is intuitively clear that the mass of a lens galaxy is not a smooth entity, we have investigated how deviation from a smooth model can influence lensing phenomena, especially the image flux ratios. To explore expectations about the level of substructure in galaxies and its influence on strong lensing, N-body simulations of a model galaxy are employed. By using the mass distribution of this model galaxy as a lens, synthetic data sets of different four image system configurations are generated. Their analysis can possibly provide evidence for the presence and strong influence of substructure in the primary lens galaxy.

In this work we investigate the gravitationally lensed system B1422+231. High--quality VLBI image positions, fluxes and shapes as well as an optical HST lens galaxy position are used. First, two simple and smooth models for the lens galaxy are applied to fit observed image positions and fluxes; no even remotely acceptable model was found. Such models also do not accurately reproduce the image shapes. In order to fit the data successfully, mass substructure has to be added to the lens, and its level is estimated. To explore expectations about the level of substructure in galaxies and its influence on strong lensing, N-body simulation results of a model galaxy are employed. By using the mass distribution of this model galaxy as a lens, synthetic data sets of different four image system configurations are generated and simple lens models are again applied to fit them. The difficulties in fitting these lens systems turn out to be similar to the case of some real gravitationally lensed systems, thus possibly providing evidence for the presence and strong influence of substructure in the primary lens galaxy.