We examine the massive colliding cluster El Gordo, one of the most massive clusters at high redshift. We use a free-form lensing reconstruction method that avoids making assumptions about the mass distribution. We use data from the RELICS program and identify new multiply lensed system candidates. The new set of constraints and free-form method provides a new independent mass estimate of this intriguing colliding cluster. Our results are found to be consistent with earlier parametric models, indirectly confirming the assumptions made in earlier work. By fitting a double gNFW profile to the lens model, and extrapolating to the virial radius, we infer a total mass for the cluster of $M_{200c}=(1.08^{+0.65}_{-0.12})\times10^{15}$M$_{\odot}$. We estimate the uncertainty in the mass due to errors in the photometric redshifts, and discuss the uncertainty in the inferred virial mass due to the extrapolation from the lens model. We also find in our lens map a mass overdensity corresponding to the large cometary tail of hot gas, reinforcing its interpretation as a large tidal feature predicted by hydrodynamical simulations that mimic El Gordo. Finally, we discuss the observed relation between the plasma and the mass map, finding that the peak in the projected mass map may be associated with a large concentration of colder gas, exhibiting possible star formation. El Gordo is one of the first clusters that will be observed with JWST, which is expected to unveil new high redshift lensed galaxies around this interesting cluster, and provide a more accurate estimation of its mass.

For studies of galaxy formation and evolution, one of the major benefits of the James Webb Space Telescope is that space-based IFUs like those on its NIRSpec and MIRI instruments will enable spatially resolved spectroscopy of distant galaxies, including spectroscopy at the scale of individual star-forming regions in galaxies that have been gravitationally lensed. In the meantime, there is only a very small subset of lensed sources where work like this is possible even with the Hubble Space Telescope's Wide Field Camera 3 infrared channel grisms. We examine two of these sources, SDSS J1723+3411 and SDSS J2340+2947, using HST WFC3/IR grism data and supporting spatially-unresolved spectroscopy from several ground-based instruments to explore the size of spatial variations in observed strong emission line ratios like O32, R23, which are sensitive to ionization parameter and metallicity, and the Balmer decrement as an indicator of reddening. We find significant spatial variation in the reddening and the reddening-corrected O32 and R23 values which correspond to spreads of a few tenths of a dex in ionization parameter and metallicity. We also find clear evidence of a negative radial gradient in star formation in SDSS J2340+2947 and tentative evidence of one in SDSS J1723+3411, though its star formation is quite asymmetric. Finally, we find that reddening can vary enough spatially to make spatially-resolved reddening corrections necessary in order to characterize gradients in line ratios and the physical conditions inferred from them, necessitating the use of space-based IFUs for future work on larger, more statistically robust samples.

We study the highly magnified arc SGAS J122651.3+215220 caused by a star-forming galaxy at $z_s=2.93$ crossing the lensing caustic cast by the galaxy cluster SDSS J1226+2152 ($z_l=0.43$), using Hubble Space Telescope observations. We report in the arc several asymmetric surface brightness features whose angular separations are a fraction of an arcsecond from the lensing critical curve and appear to be highly but unequally magnified image pairs of underlying compact sources, with one brightest pair having clear asymmetry consistently across four filters. One explanation of unequal magnification is microlensing by intracluster stars, which induces independent flux variations in the images of individual or groups of source stars in the lensed galaxy. For a second possibility, intracluster dark matter subhalos invisible to telescopes effectively perturb lensing magnifications near the critical curve and give rise to persistently unequal image pairs. Our modeling suggests, at least for the most prominent identified image pair, that the microlensing hypothesis is in tension with the absence of notable asymmetry variation over a six-year baseline, while subhalos of $\sim 10^6$--$10^8\,M_\odot$ anticipated from structure formation with Cold Dark Matter typically produce stationary and sizable asymmetries. We judge that observations at additional times and more precise lens models are necessary to stringently constrain temporal variability and robustly distinguish between the two explanations. The arc under this study is a scheduled target of a Director's Discretionary Early Release Science program of the James Webb Space Telescope, which will provide deep images and a high-resolution view with integral field spectroscopy.

Extensive surveys with the \textit{Hubble Space Telescope} (HST) over the past decade, targeting some of the most massive clusters in the sky, have uncovered dozens of galaxy-cluster strong lenses. The massive cluster strong-lens scale is typically $\theta_{E}\sim10\arcsec$ to $\sim30-35\arcsec$, with only a handful of clusters known with Einstein radii $\theta_{E}\sim40\arcsec$ or above (for $z_{source}=2$, nominally). Here we report another very large cluster lens, RXC J0032.1+1808 ($z=0.3956$), the second richest cluster in the redMapper cluster catalog and the 85th most massive cluster in the Planck Sunyaev-Zel'dovich catalog. With our Light-Traces-Mass and fully parametric (dPIEeNFW) approaches, we construct strong lensing models based on 18 multiple images of 5 background galaxies newly identified in the \textit{Hubble} data mainly from the \textit{Reionization Lensing Cluster Survey} (RELICS), in addition to a known sextuply imaged system in this cluster. Furthermore, we compare these models to Lenstool and GLAFIC models that were produced independently as part of the RELICS program. All models reveal a large effective Einstein radius of $\theta_{E}\simeq40\arcsec$ ($z_{source}=2$), owing to the obvious concentration of substructures near the cluster center. Although RXC J0032.1+1808 has a very large critical area and high lensing strength, only three magnified high-redshift candidates are found within the field targeted by RELICS. Nevertheless, we expect many more high-redshift candidates will be seen in wider and deeper observations with \textit{Hubble} or \emph{JWST}. Finally, the comparison between several algorithms demonstrates that the total error budget is largely dominated by systematic uncertainties.

Recent observations have revealed the presence of strong CIII] emission (EW$_{\rm{CIII]}}>20$ \r{A}) in $z>6$ galaxies, the origin of which remains unclear. In an effort to understand the nature of these line emitters, we have initiated a survey targeting CIII] emission in gravitationally-lensed reionization era analogs identified in HST imaging of clusters from the RELICS survey. Here we report initial results on four galaxies selected to have low stellar masses (2-8$\times$10$^7$ M$_\odot$) and J$_{125}$-band flux excesses indicative of intense [OIII]+H$\beta$ emission (EW$_{\rm{[OIII]+H\beta}}$=500-2000 \r{A}), similar to what has been observed at $z>6$. We detect CIII] emission in three of the four sources, with the CIII] EW reaching values seen in the reionization era (EW$_{\rm{CIII]}}\simeq 17-22$ \r{A}) in the two sources with the strongest optical line emission (EW$_{\rm{[OIII]+H\beta}}\simeq 2000$ \r{A}). We have obtained a Magellan/FIRE near-infrared spectrum of the strongest CIII] emitter in our sample, revealing gas that is both metal poor and highly ionized. Using photoionization models, we are able to simultaneously reproduce the intense CIII] and optical line emission for extremely young (2-3 Myr) and metal poor (0.06-0.08 Z$_\odot$) stellar populations, as would be expected after a substantial upturn in the SFR of a low mass galaxy. The sources in this survey are among the first for which CIII] has been used as the primary means of redshift confirmation. We suggest that it should be possible to extend this approach to $z>6$ with current facilities, using CIII] to measure redshifts of objects with IRAC excesses indicating EW$_{\rm{[OIII]+H\beta}}\simeq 2000$ \r{A}, providing a method of spectroscopic confirmation independent of Ly$\alpha$.

Our ability to study the properties of the interstellar medium (ISM) in the earliest galaxies will rely on emission line diagnostics at rest-frame ultraviolet (UV) wavelengths. In this work, we identify metallicity-sensitive diagnostics using UV emission lines. We compare UV-derived metallicities with standard, well-established optical metallicities using a sample of galaxies with rest-frame UV and optical spectroscopy. We find that the He2-O3C3 diagnostic (He II 1640 / C III 1906,1909 vs. O III 1666 / C III 1906,1909) is a reliable metallicity tracer, particularly at low metallicity (12+log(O/H) < 8), where stellar contributions are minimal. We find that the Si3-O3C3 diagnostic (Si III 1883 / C III 1906,1909 vs. O III 1666 / C III 1906,1909) is a reliable metallicity tracer, though with large scatter (0.2-0.3 dex), which we suggest is driven by variations in gas-phase abundances. We find that the C4-O3C3 diagnostic (C IV 1548,1550 / O III 1666 vs. O III 1666 / C III 1906,1909) correlates poorly with optically-derived metallicities. We discuss possible explanations for these discrepant metallicity determinations, including the hardness of the ionizing spectrum, contribution from stellar wind emission, and non-solar-scaled gas-phase abundances. Finally, we provide two new UV oxygen abundance diagnostics, calculated from polynomial fits to the model grid surface in the He2-O3C3 and Si3-O3C3 diagrams.

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.

We present a new galaxy cluster lens modeling approach, hybrid-Lenstool, that is implemented in the publicly available modeling software Lenstool. hybrid-Lenstool combines a parametric approach to model the core of the cluster, and a non-parametric (free-form) approach to model the outskirts. hybrid-Lenstool optimizes both strong- and weak-lensing constraints simultaneously (Joint-Fit), providing a self-consistent reconstruction of the cluster mass distribution on all scales. In order to demonstrate the capabilities of the new algorithm, we tested it on a simulated cluster. hybrid-Lenstool yields more accurate reconstructed mass distributions than the former Sequential-Fit approach where the parametric and the non-parametric models are optimized successively. Indeed, we show with the simulated cluster that the mass density profile reconstructed with a Sequential-Fit deviates form the input by $2-3\sigma$ at all scales while the Joint-Fit gives a profile that is within $1-1.5\sigma$ of the true value. This gain in accuracy is consequential for recovering mass distributions exploiting cluster lensing and therefore for all applications of clusters as cosmological probes. Finally we found that the Joint-Fit approach yields shallower slope of the inner density profile than the Sequential-Fit approach, thus revealing possible biases in previous lensing studies.

We present a catalog of galaxy cluster candidates detected in 100 square degrees surveyed with the SPTpol receiver on the South Pole Telescope. The catalog contains 89 candidates detected with a signal-to-noise ratio greater than 4.6. The candidates are selected using the Sunyaev-Zel'dovich effect at 95 and 150 GHz. Using both space- and ground-based optical and infrared telescopes, we have confirmed 81 candidates as galaxy clusters. We use these follow-up images and archival images to estimate photometric redshifts for 66 galaxy clusters and spectroscopic observations to obtain redshifts for 13 systems. An additional 2 galaxy clusters are confirmed using the overdensity of near-infrared galaxies only, and are presented without redshifts. We find that 15 candidates (18% of the total sample) are at redshift of $z \geq 1.0$, with a maximum confirmed redshift of $z_{\rm{max}} = 1.38 \pm 0.10$. We expect this catalog to contain every galaxy cluster with $M_{500c} > 2.6 \times 10^{14} M_\odot h^{-1}_{70}$ and $z > 0.25$ in the survey area. The mass threshold is approximately constant above $z = 0.25$, and the complete catalog has a median mass of approximately $ M_{500c} = 2.7 \times 10^{14} M_\odot h^{-1}_{70}$. Compared to previous SPT works, the increased depth of the millimeter-wave data (11.2 and 6.5 $\mu$K-arcmin at 95 and 150 GHz, respectively) makes it possible to find more galaxy clusters at high redshift and lower mass.

We describe the observations and resultant galaxy cluster catalog from the 2770 deg$^2$ SPTpol Extended Cluster Survey (SPT-ECS). Clusters are identified via the Sunyaev-Zel'dovich (SZ) effect, and confirmed with a combination of archival and targeted follow-up data, making particular use of data from the Dark Energy Survey (DES). With incomplete followup we have confirmed as clusters 244 of 266 candidates at a detection significance $\xi \ge 5$ and an additional 204 systems at $4<\xi<5$. The confirmed sample has a median mass of $M_{500c} \sim {4.4 \times 10^{14} M_\odot h_{70}^{-1}}$, a median redshift of $z=0.49$, and we have identified 44 strong gravitational lenses in the sample thus far. Radio data are used to characterize contamination to the SZ signal; the median contamination for confirmed clusters is predicted to be $\sim$1% of the SZ signal at the $\xi>4$ threshold, and $<4\%$ of clusters have a predicted contamination $>10\% $ of their measured SZ flux. We associate SZ-selected clusters, from both SPT-ECS and the SPT-SZ survey, with clusters from the DES redMaPPer sample, and find an offset distribution between the SZ center and central galaxy in general agreement with previous work, though with a larger fraction of clusters with significant offsets. Adopting a fixed Planck-like cosmology, we measure the optical richness-to-SZ-mass ($\lambda-M$) relation and find it to be 28% shallower than that from a weak-lensing analysis of the DES data---a difference significant at the 4 $\sigma$ level---with the relations intersecting at $\lambda=60$ . The SPT-ECS cluster sample will be particularly useful for studying the evolution of massive clusters and, in combination with DES lensing observations and the SPT-SZ cluster sample, will be an important component of future cosmological analyses.

We present spatially-resolved echelle spectroscopy of an intervening MgII-FeII-MgI absorption-line system detected at $z_{\rm abs}=0.73379$ toward the giant gravitational arc PSZ1 G311.65-18.48. The absorbing gas is associated to an inclined disk-like star-forming galaxy, whose major axis is aligned with the two arc-segments reported here. We probe in absorption the galaxy's extended disk continuously, at $\approx 3$ kpc sampling, from its inner region out to $15\times$ the optical radius. We detect strong ($W_0^{2796}>0.3$ \r{A}) coherent absorption along $13$ independent positions at impact parameters $D=0$--$29$ kpc on one side of the galaxy, and no absorption at $D=28$--$57$ kpc on the opposite side (all de-lensed distances at $z_{\rm abs}$). We show that: (1) the gas distribution is anisotropic; (2) $W_0^{2796}$, $W_0^{2600}$, $W_0^{2852}$, and the ratio $W_0^{2600}\!/W_0^{2796}$, all anti-correlate with $D$; (3) the $W_0^{2796}$-$D$ relation is not cuspy and exhibits significantly less scatter than the quasar-absorber statistics; (4) the absorbing gas is co-rotating with the galaxy out to $D \lesssim 20$ kpc, resembling a `flat' rotation curve, but at $D\gtrsim 20$ kpc velocities decline below the expectations from a 3D disk-model extrapolated from the nebular [OII] emission. These signatures constitute unambiguous evidence for rotating extra-planar diffuse gas, possibly also undergoing enriched accretion at its edge. Arguably, we are witnessing some of the long-sought processes of the baryon cycle in a single distant galaxy expected to be representative of such phenomena.

We present hitherto the largest sample of gas-phase metallicity radial gradients measured at sub-kiloparsec resolution in star-forming galaxies in the redshift range of $z\in[1.2, 2.3]$. These measurements are enabled by the synergy of slitless spectroscopy from the Hubble Space Telescope near-infrared channels and the lensing magnification from foreground galaxy clusters. Our sample consists of 79 galaxies with stellar mass ranging from 10$^7$ to 10$^{10}$ $M_\odot$, instantaneous star-formation rate in the range of [1, 100] $M_\odot$/yr, and global metallicity [$\frac{1}{12}$, 2] solar. At 2-$\sigma$ confidence level, 14/79 galaxies in our sample show negative radial gradients, whereas 7/79 show inverted gradients. Combining ours and all other metallicity gradients obtained at similar resolution currently available in the literature, we measure a negative mass dependence of $\Delta\log({\rm O/H})/\Delta r~ [\mathrm{dex~kpc^{-1}}] = (-0.020\pm0.007) + (-0.014\pm0.008) \log(M_\ast/10^{9.4}M_\odot)$ with the intrinsic scatter being $\sigma=0.060\pm0.006$ over four orders of magnitude in stellar mass. Our result is consistent with strong feedback, not secular processes, being the primary governor of the chemo-structural evolution of star-forming galaxies during the disk mass assembly at cosmic noon. We also find that the intrinsic scatter of metallicity gradients increases with decreasing stellar mass and increasing specific star-formation rate. This increase in the intrinsic scatter is likely caused by the combined effect of cold-mode gas accretion and merger-induced starbursts, with the latter more predominant in the dwarf mass regime of $M_\ast\lesssim10^9 M_\odot$.

During the epoch of reionisation, neutral gas in the early Universe was ionized by hard ultraviolet radiation emitted by young stars in the first galaxies. To do so, ionizing ultraviolet photons must escape from the host galaxy. We present Hubble Space Telescope observations of the gravitationally lensed galaxy PSZ1-ARC G311.6602-18.4624, revealing bright, multiply-imaged ionizing photon escape from a compact star-forming region through a narrow channel in an optically thick gas. The gravitational lensing magnification shows how ionizing photons escape this galaxy, contributing to the re-ionization of the Universe. The multiple sight lines to the source probe absorption by intergalactic neutral hydrogen on scales of no more than a few hundred, perhaps even less than ten, parsec.

We present an analysis of the mass distribution inferred from strong lensing by SPT-CL J0356-5337, a cluster of galaxies at redshift z = 1.0359 revealed in the follow-up of the SPT-SZ clusters. The cluster has an Einstein radius of Erad=14 for a source at z = 3 and a mass within 500 kpc of M_500kpc = 4.0+-0.8x10^14Msol. Our spectroscopic identification of three multiply-imaged systems (z = 2.363, z = 2.364, and z = 3.048), combined with HST F606W-band imaging allows us to build a strong lensing model for this cluster with an rms of <0.3'' between the predicted and measured positions of the multiple images. Our modeling reveals a two-component mass distribution in the cluster. One mass component is dominated by the brightest cluster galaxy and the other component, separated by ~170 kpc, contains a group of eight red elliptical galaxies confined in a ~9'' (~70 kpc) diameter circle. We estimate the mass ratio between the two components to be between 1:1.25 and 1:1.58. In addition, spectroscopic data reveal that these two near-equal mass cores have only a small velocity difference of 300 km/s between the two components. This small radial velocity difference suggests that most of the relative velocity takes place in the plane of the sky, and implies that SPT-CL J0356-5337 is a major merger with a small impact parameter seen face-on. We also assess the relative contributions of galaxy-scale halos to the overall mass of the core of the cluster and find that within 800 kpc from the brightest cluster galaxy about 27% of the total mass can be attributed to visible and dark matter associated with galaxies, whereas only 73% of the total mass in the core comes from cluster-scale dark matter halos.

In the past decade, our understanding of how stars and galaxies formed during the first 5 billion years after the Big Bang has been revolutionized by observations that leverage gravitational lensing by intervening masses, which act as natural cosmic telescopes to magnify background sources. Previous studies have harnessed this effect to probe the distant universe at ultraviolet, optical, infrared and millimeter wavelengths. However, strong lensing studies of young, star-forming galaxies have never extended into X-ray wavelengths, which uniquely trace high-energy phenomena. Here we report an X-ray detection of star formation in a highly magnified, strongly lensed galaxy. This lensed galaxy, seen during the first third of the history of the Universe, is a low--mass, low--metallicity starburst with elevated X-ray emission, and is a likely analog to the first generation of galaxies. Our measurements yield insight into the role that X-ray emission from stellar populations in the first generation of galaxies may play in re-ionizing the Universe. This observation paves the way for future strong lensing-assisted X-ray studies of distant galaxies reaching orders of magnitude below the detection limits of current deep fields, and previews the depths that will be attainable with future X-ray observatories.

We report the first sub-kiloparsec spatial resolution measurements of strongly inverted gas-phase metallicity gradients in two dwarf galaxies at $z$$\sim$2. The galaxies have stellar masses $\sim$$10^9M_\odot$, specific star-formation rate $\sim$20 Gyr$^{-1}$, and global metallicity $12+\log({\rm O/H})\sim8.1$ (1/4 solar), assuming the Maiolino et al. (2008) strong line calibrations of OIII/Hb and OII/Hb. Their metallicity radial gradients are measured to be highly inverted, i.e., 0.122$\pm$0.008 and 0.111$\pm$0.017 dex/kpc, which is hitherto unseen at such small masses in similar redshift ranges. From the Hubble Space Telescope observations of the source nebular emission and stellar continuum, we present the 2-dimensional spatial maps of star-formation rate surface density, stellar population age, and gas fraction, which show that our galaxies are currently undergoing rapid mass assembly via disk inside-out growth. More importantly, using a simple chemical evolution model, we find that the gas fractions for different metallicity regions cannot be explained by pure gas accretion. Our spatially resolved analysis based on a more advanced gas regulator model results in a spatial map of net gaseous outflows, triggered by active central starbursts, that potentially play a significant role in shaping the spatial distribution of metallicity by effectively transporting stellar nucleosynthesis yields outwards. The relation between wind mass loading factors and stellar surface densities measured in different regions of our galaxies shows that a single type of wind mechanism, driven by either energy or momentum conservation, cannot explain the entire galaxy. These sources present a unique constraint on the effects of gas flows on the early phase of disk growth from the perspective of spatially resolved chemical evolution within individual systems.

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.

The ionizing photons produced by massive stars are key actors in galaxy evolution. Ionizing photon production and escape is poorly understood. Improved space-based, spatially-resolved, multiplexed spectroscopic capabilities covering observed wavelengths of 1000 to 3000 Angstrom, in concert with spectroscopy from the ELTs and JWST, would lead to definitive answers as to how ionizing photons are produced and leaked, what populations of galaxies are responsible for ionizing photon leakage, what determines whether escape is possible, and how ionizing galaxy populations evolve over cosmic time.

We infer the properties of massive star populations using the far-ultraviolet stellar continua of 61 star-forming galaxies: 42 at low-z observed with HST and 19 at z~2 from the Megasaura sample. We fit each stellar continuum with a linear combination of up to 50 single age and single metallicity Starburst99 models. From these fits, we derive light-weighted ages and metallicities, which agree with stellar wind and photospheric spectral features, and infer the spectral shapes and strengths of the ionizing continua. Inferred light-weighted stellar metallicities span 0.05-1.5 Z$_\odot$ and are similar to the measured nebular metallicities. We quantify the ionizing continua using the ratio of the ionizing flux at 900\AA\ to the non-ionizing flux at 1500\AA\ and demonstrate the evolution of this ratio with stellar age and metallicity using theoretical single burst models. These single burst models only match the inferred ionizing continua of half of the sample, while the other half are described by a mixture of stellar ages. Mixed age populations produce stronger and harder ionizing spectra than continuous star formation histories, but, contrary to previous studies that assume constant star formation, have similar stellar and nebular metallicities. Stellar population age and metallicity affect the far-UV continua in different and distinguishable ways; assuming a constant star formation history diminishes the diagnostic power. Finally, we provide simple prescriptions to determine the ionizing photon production efficiency ($\xi_{ion}$) from the stellar population properties. $\xi_{ion}$ has a range of log($\xi_{ion})=24.4-25.7$ Hz erg$^{-1}$ that depends on stellar age, metallicity, star formation history, and contributions from binary star evolution. These stellar population properties must be observationally determined to determine the number of ionizing photons generated by massive stars.

We present new, deep observations of the Phoenix cluster from the Chandra X-ray Observatory, the Hubble Space Telescope, and the Karl Jansky Very Large Array. These data provide an order of magnitude improvement in depth and/or angular resolution at X-ray, optical, and radio wavelengths, yielding an unprecedented view of the core of the Phoenix cluster. We find that the one-dimensional temperature and entropy profiles are consistent with expectations for pure-cooling hydrodynamic simulations and analytic descriptions of homogeneous, steady-state cooling flow models. In the inner ~10 kpc, the cooling time is shorter by an order of magnitude than any other known cluster, while the ratio of the cooling time to freefall time approaches unity, signaling that the ICM is unable to resist multiphase condensation on kpc scales. When we consider the thermodynamic profiles in two dimensions, we find that the cooling is highly asymmetric. The bulk of the cooling in the inner ~20 kpc is confined to a low-entropy filament extending northward from the central galaxy. We detect a substantial reservoir of cool (10^4 K) gas (as traced by the [OII] doublet), which is coincident with the low-entropy filament. The bulk of this cool gas is draped around and behind a pair of X-ray cavities, presumably bubbles that have been inflated by radio jets, which are detected for the first time on kpc scales. These data support a picture in which AGN feedback is promoting the formation of a multiphase medium via a combination of ordered buoyant uplift and locally enhanced turbulence. These processes ought to counteract the tendency for buoyancy to suppress condensation, leading to rapid cooling along the jet axis. The recent mechanical outburst has sufficient energy to offset cooling, and appears to be coupling to the ICM via a cocoon shock, raising the entropy in the direction orthogonal to the radio jets.

We present strong gravitational lensing models for 37 galaxy clusters from the SDSS Giant Arcs Survey. We combine data from multi-band Hubble Space Telescope WFC3imaging, with ground-based imaging and spectroscopy from Magellan, Gemini, APO, and MMT, in order to detect and spectroscopically confirm new multiply-lensed background sources behind the clusters. We report spectroscopic or photometric redshifts of sources in these fields, including cluster galaxies and background sources. Based on all available lensing evidence, we construct and present strong lensing mass models for these galaxy clusters.

We present the first spatially resolved analysis of rest-frame optical and UV imaging and spectroscopy for a lensed galaxy at z = 2.39 hosting a Seyfert active galactic nucleus (AGN). Proximity to a natural guide star has enabled high signal-to-noise VLT SINFONI + adaptive optics observations of rest-frame optical diagnostic emission lines, which exhibit an underlying broad component with FWHM ~ 700 km/s in both the Balmer and forbidden lines. Measured line ratios place the outflow robustly in the region of the ionization diagnostic diagrams associated with AGN. This unique opportunity - combining gravitational lensing, AO guiding, redshift, and AGN activity - allows for a magnified view of two main tracers of the physical conditions and structure of the interstellar medium in a star-forming galaxy hosting a weak AGN at cosmic noon. By analyzing the spatial extent and morphology of the Ly-alpha and dust-corrected H-alpha emission, disentangling the effects of star formation and AGN ionization on each tracer, and comparing the AGN induced mass outflow rate to the host star formation rate, we find that the AGN does not significantly impact the star formation within its host galaxy.

Strong gravitational lensing by clusters of galaxies probes the mass distribution at the core of each cluster and magnifies the universe behind it. MACS J0417.5-1154 at z=0.443 is one of the most massive clusters known based on weak lensing, X-ray, and Sunyaev-Zel'dovich analyses. Here we compute a strong lens model of MACS J0417 based on Hubble Space Telescope imaging observations collected, in part, by the Reionization Lensing Cluster Survey (RELICS), and recently reported spectroscopic redshifts from the MUSE instrument on the Very Large Telescope (VLT). We measure an Einstein radius of ThetaE=36'' at z = 9 and a mass projected within 200 kpc of M(200 kpc) = 1.78+0.01-0.03x10**14Msol. Using this model, we measure a ratio between the mass attributed to cluster-member galaxy halos and the main cluster halo of order 1:100. We assess the probability to detect magnified high-redshift galaxies in the field of this cluster, both for comparison with RELICS HST results and as a prediction for the James Webb Space Telescope (JWST) Guaranteed Time Observations upcoming for this cluster. Our lensing analysis indicates that this cluster has similar lensing strength to other clusters in the RELICS program. Our lensing analysis predicts a detection of at least a few z~6-8 galaxies behind this cluster, at odds with a recent analysis that yielded no such candidates in this field. Reliable strong lensing models are crucial for accurately predicting the intrinsic properties of lensed galaxies. As part of the RELICS program, our strong lensing model produced with the Lenstool parametric method is publicly available through the Mikulski Archive for Space Telescopes (MAST).

We present observations of two strong-lensing galaxy clusters located within the $90$ per cent credible sky localization maps released following LIGO-Virgo's discovery of the binary black hole (BH-BH) gravitational wave (GW) source GW170814. Our objectives were (1) to search for candidate electromagnetic (EM) counterparts to GW170814 under the hypothesis that it was strongly-lensed, and thus more distant and less massive than inferred by LIGO-Virgo, and (2) to demonstrate the feasibility of rapid target of opportunity observations to search for faint lensed transient point sources in crowded cluster cores located within GW sky localizations. Commencing $20$ hours after discovery, and continuing over $12$ nights, we observed Abell 3084 ($z=0.22$) and SMACSJ0304.3$-$4401 ($z=0.46$) with GMOS on the Gemini-South telescope, and Abell 3084 with MUSE on ESO's Very Large Telescope. We detect no candidate EM counterparts in these data. Calibration of our photometric analysis methods using simulations yield $5\sigma$ detection limits for transients in difference images of the cores of these clusters of $i=25$. This is the most sensitive photometric search to date for counterparts to GW sources, and rules out the possibility that GW170814 was lensed by these clusters with a kilonova-like EM counterpart. Based on the detector frame masses of the compact objects, and assuming that at least one Neutron Star (NS) is required in the merging system to produce a kilonova-like counterpart, implies that GW170814 was neither a NS-NS nor NS-BH merger at $z>8$ lensed by either of these clusters. Also, in the first ever emission line search for counterparts to GW sources, we detected no lines down to a $5\sigma$ detection limit of $5\times10^{-17}{\rm erg\,s^{-1}\,cm^{-2}}$.

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.

We present $Hubble\ Space\ Telescope$ ($HST$) imaging and grism spectroscopy of a strongly lensed LIRG at $z=0.816$, SGAS 143845.1$+$145407, and use the magnification boost of gravitational lensing to study the distribution of star formation throughout this galaxy. Based on the $HST$ imaging data, we create a lens model for this system; we compute the mass distribution and magnification map of the $z=0.237$ foreground lens. We find that the magnification of the lensed galaxy ranges between $2$ and $10$, with a total magnification (measured over all the images of the source) of $\mu=11.8^{+4.6}_{-2.4}$. We find that the total projected mass density within $\sim34$ kpc of the brightest cluster galaxy is $6.0^{+0.3}_{-0.7}\times10^{12}\,M_{\odot}$. Using the lens model we create a source reconstruction for SGAS 143845.1$+$145407, which paired with a faint detection of H$\alpha$ in the grism spectroscopy, allows us to finally comment directly on the distribution of star formation in a $z\sim1$ LIRG. We find widespread star formation across this galaxy, in agreement with the current understanding of these objects. However, we note a deficit of H$\alpha$ emission in the nucleus of SGAS 143845.1$+$145407, likely due to dust extinction.

Strong gravitational lensing (SL) is a powerful means to map the distribution of dark matter. In this work, we perform a SL analysis of the prominent X-ray cluster RXJ0152.7-1357 (z=0.83, also known as CL 0152.7-1357) in \textit{Hubble Space Telescope} images, taken in the framework of the Reionization Lensing Cluster Survey (RELICS). On top of a previously known $z=3.93$ galaxy multiply imaged by RXJ0152.7-1357, for which we identify an additional multiple image, guided by a light-traces-mass approach we identify seven new sets of multiply imaged background sources lensed by this cluster, spanning the redshift range [1.79-3.93]. A total of 25 multiple images are seen over a small area of ~0.4 $arcmin^2$, allowing us to put relatively high-resolution constraints on the inner matter distribution. Although modestly massive, the high degree of substructure together with its very elongated shape make RXJ0152.7-1357 a very efficient lens for its size. This cluster also comprises the third-largest sample of z~6-7 candidates in the RELICS survey. Finally, we present a comparison of our resulting mass distribution and magnification estimates with those from a Lenstool model. These models are made publicly available through the MAST archive.

We present first results from a galaxy population study in the highest redshift galaxy clusters identified in the 2500 deg$^2$ South Pole Telescope Sunyaev Zel'dovich effect (SPT-SZ) survey. The cluster selection is to first order independent of galaxy properties, making the SPT-SZ sample particularly well suited for cluster galaxy population studies. We carry out a 4-band imaging campaign with the {\it Hubble} and {\it Spitzer} Space Telescopes of the five $z\gtrsim 1.4$, S/N$_{SZE}>$5 clusters, that are among the rarest most massive clusters known at this redshift. All five show clear overdensities of red galaxies whose colors agree with the initial cluster redshift estimates. The highest redshift cluster in this sample, SPT-CLJ0459-4947 at $z\sim1.72$, is the most distant $M_{500}>10^{14}~M_{\odot}$ ICM-selected cluster discovered thus far, and is one of only three known clusters in this mass range at $z\gtrsim 1.7$, regardless of selection. Based on UVJ-like photometric classification of quiescent and star-forming galaxies, the passive fraction in the cluster central regions ($r/r_{500}<0.7$) is higher than in the field at the same redshift, with corresponding environmental quenching efficiencies typically in the range $\sim0.5-0.8$ for stellar masses $\log(M/M_{\odot})>10.85$. We have explored the impact of emission from star formation on the selection of this sample, concluding that all five clusters studied here would still have been detected with S/N$_{SZE}>$5, even if they had the same passive fraction as measured in the field. Our results thus point towards an efficient suppression of star formation in massive galaxies in the central regions of the most massive clusters, occurring already earlier than $z\sim1.5$. [Abridged]

We present a multi-wavelength analysis of the core of the massive galaxy cluster MACS\,J0417.5-1154 ($z = 0.441$; MACS\;J0417). Our analysis takes advantage of VLT/MUSE observations which allow the spectroscopic confirmation of three strongly-lensed systems. One of these, nick-named \emph{The Doughnut}, consists of three complete images of a complex ring galaxy at $z = 0.8718$ and a fourth, partial and radial image close to the Brightest Cluster Galaxy (BCG) only discernible thanks to its strong [OII] line emission. The best-fit mass model (rms of 0.38\arcsec) yields a two-dimensional enclosed mass of $M({\rm R < 200\,kpc}) = (1.77\pm0.03)\times10^{14}\,\msun$ and almost perfect alignment between the peaks of the BCG light and the dark matter of ($0.5\pm0.5$)\arcsec . Our finding that a significant misalignment results when the radial image of \emph{The Doughnut} is omitted serves as an important caveat for studies of BCG-dark matter offsets in galaxy clusters. Using \emph{Chandra} data to map the intra-cluster gas, we observe an offset between the gas and dark-matter peaks of ($1.7\pm0.5$)\arcsec, and excellent alignment of the X-ray peak with the location of optical emission line associated with the BCG. We interpret all observational evidence in the framework of on-going merger activity, noting specifically that the coincidence between the gas peak and the peak of blue light from the BCG may be evidence of dense, cold gas leading to direct star formation. We use the surface area $\sigma_{\mu}$ above a given magnification factor $\mu$ as a metric to estimate the lensing power of MACS\,J0417. We obtain $\sigma(\mu > 3) = 0.22$\,arcmin$^2$ confirming MACS\,J0417 as an efficient gravitational lens. Finally, we discuss the differences between our mass model and Mahler et al. (2018).

We present measurements of the X-ray observables of the intra-cluster medium (ICM), including luminosity $L_X$, ICM mass $M_{ICM}$, emission-weighted mean temperature $T_X$, and integrated pressure $Y_X$, that are derived from XMM-Newton X-ray observations of a Sunyaev-Zel'dovich Effect (SZE) selected sample of 59 galaxy clusters from the South Pole Telescope SPT-SZ survey that span the redshift range of $0.20 < z < 1.5$. We constrain the best-fit power law scaling relations between X-ray observables, redshift, and halo mass. The halo masses are estimated based on previously published SZE observable to mass scaling relations, calibrated using information that includes the halo mass function. Employing SZE-based masses in this sample enables us to constrain these scaling relations for massive galaxy clusters ($M_{500}\geq 3 \times10^{14}$ $M_\odot$) to the highest redshifts where these clusters exist without concern for X-ray selection biases. We find that the mass trends are steeper than self-similarity in all cases, and with $\geq 2.5{\sigma}$ significance in the case of $L_X$ and $M_{ICM}$. The redshift trends are consistent with the self-similar expectation, but the uncertainties remain large. Core-included scaling relations tend to have steeper mass trends for $L_X$. There is no convincing evidence for a redshift-dependent mass trend in any observable. The constraints on the amplitudes of the fitted scaling relations are currently limited by the systematic uncertainties on the SZE-based halo masses, however the redshift and mass trends are limited by the X-ray sample size and the measurement uncertainties of the X-ray observables.

The galaxy phase-space distribution in galaxy clusters provides insights into the formation and evolution of cluster galaxies, and it can also be used to measure cluster mass profiles. We present a dynamical study based on $\sim$3000 passive, non-emission line cluster galaxies drawn from 110 galaxy clusters. The galaxy clusters were selected using the Sunyaev-Zel'dovich effect (SZE) in the 2500~deg$^2$ SPT-SZ survey and cover the redshift range $0.2 < z < 1.3$. We model the clusters using the Jeans equation, while adopting NFW mass profiles and a broad range of velocity dispersion anisotropy profiles. The data prefer velocity dispersion anisotropy profiles that are approximately isotropic near the center and increasingly radial toward the cluster virial radius, and this is true for all redshifts and masses we study. The pseudo-phase-space density profile of the passive galaxies is consistent with expectations for dark matter particles and subhalos from cosmological $N$-body simulations. The dynamical mass constraints are in good agreement with external mass estimates of the SPT cluster sample from either weak lensing, velocity dispersions, or X-ray $Y_X$ measurements. However, the dynamical masses are lower (at the 2.2$\sigma$ level) when compared to the mass calibration favored when fitting the SPT cluster data to a $\Lambda$CDM model with external cosmological priors, including CMB anisotropy data from Planck. The discrepancy grows with redshift, where in the highest redshift bin the ratio of dynamical to SPT+Planck masses is $\eta=0.63^{+0.13}_{-0.08}\pm0.06$ (statistical and systematic), corresponding to a $2.6\sigma$ discrepancy.

We present a strong-lensing analysis of four massive galaxy clusters imaged with the Hubble Space Telescope in the Reionization Lensing Cluster Survey. We use a Light-Traces-Mass technique to uncover sets of multiply images and constrain the mass distribution of the clusters. These mass models are the first published for Abell S295 and MACS J0159.8-0849, and are improvements over previous models for Abell 697 and MACS J0025.4-1222. Our analysis for MACS J0025.4-1222 and Abell S295 shows a bimodal mass distribution supporting the merger scenarios proposed for these clusters. The updated model for MACS J0025.4-1222 suggests a substantially smaller critical area than previously estimated. For MACS J0159.8-0849 and Abell 697 we find a single peak and relatively regular morphology, revealing fairly relaxed clusters. Despite being less prominent lenses, three of these clusters seem to have lensing strengths, i.e. cumulative area above certain magnification, similar to the Hubble Frontier Fields clusters (e.g., A($\mu>5$) $\sim 1-3$ arcmin$^2$, A($\mu>10$) $\sim 0.5-1.5$ arcmin$^2$), which in part can be attributed to their merging configurations. We make our lens models publicly available through the Mikulski Archive for Space Telescopes. Finally, using Gemini-N/GMOS spectroscopic observations we detect a single emission line from a high-redshift $J_{125}\simeq25.7$ galaxy candidate lensed by Abell 697. While we cannot rule out a lower-redshift solution, we interpret the line as Ly$\alpha$ at $z=5.800\pm 0.001$, in agreement with its photometric redshift and dropout nature. Within this scenario we measure a Ly$\alpha$ rest-frame equivalent width of $52\pm22$ \AA, and an observed Gaussian width of $117\pm 15$ km/s.

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.

Galaxy clusters can act as gravitational lenses to magnify the universe behind them, allowing us to see deep into the early universe. The Hubble Space Telescope Frontier Fields program (Lotz et al. 2017) uses six galaxy clusters imaged by Hubble to discover and study galaxies at z ~ 5-10. Seven independent teams developed lens models and derived magnifications for each galaxy cluster, based on positional and redshift constraints from the best available data at the time. In this work we evaluate ten models for MACSJ0416.1-2403, which were made public in 2015, by contrasting them with new spectroscopic redshifts that were measured in 2016 (Caminha et al. 2017). We developed an independent comparison method that uses the source plane root-mean-square as a metric of lensing model performance. Our analysis quantifies the ability of models to predict unknown multiple images. We examine the source plane scatter of multiply-imaged systems and explore dependence of the scatter on the location and the redshift of the background sources. The analysis we present evaluates the performance of the different algorithms in the specific case of the MACSJ0416.1-2403 models.

We present spectroscopic confirmation of five galaxy clusters at $1.25 < \textit{z} < 1.5$, discovered in the $2500$ deg$^{2}$ South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) survey. These clusters, taken from a mass-limited sample with a nearly redshift independent selection function, have multi-wavelength follow-up imaging data from the X-ray to near-infrared, and currently form the most homogeneous massive high-redshift cluster sample known. We identify $44$ member galaxies, along with $25$ field galaxies, among the five clusters, and describe the full set of observations and data products from Magellan/LDSS3 multi-object spectroscopy of these cluster fields. We briefly describe the analysis pipeline, and present ensemble analyses of cluster member galaxies that demonstrate the reliability of the measured redshifts. We report $\textit{z} = 1.259, 1.288, 1.316, 1.401$ and $1.474$ for the five clusters from a combination of absorption-line (Ca II H$\&$K doublet - $3968,3934$ {\AA}) and emission-line ([OII] $3727,3729$ {\AA}) spectral features. Moreover, the calculated velocity dispersions yield dynamical cluster masses in good agreement with SZ masses for these clusters. We discuss the velocity and spatial distributions of passive and [OII]-emitting galaxies in these clusters, showing that they are consistent with velocity segregation and biases observed in lower redshift SPT clusters. We identify modest [OII] emission and pronounced CN and H$\delta$ absorption in a stacked spectrum of $28$ passive galaxies with Ca II H$\&$K-derived redshifts. This work increases the number of spectroscopically-confirmed SZ-selected galaxy clusters at $\textit{z} > 1.25$ from three to eight, further demonstrating the efficacy of SZ selection for the highest redshift massive clusters, and enabling detailed study of these systems.

We estimate total mass ($M_{500}$), intracluster medium (ICM) mass ($M_{\mathrm{ICM}}$) and stellar mass ($M_{\star}$) in a Sunyaev-Zel'dovich effect (SZE) selected sample of 91 galaxy clusters with masses $M_{500}\gtrsim2.5\times10^{14}M_{\odot}$ and redshift $0.2 < z < 1.25$ from the 2500 deg$^2$ South Pole Telescope SPT-SZ survey. The total masses $M_{500}$ are estimated from the SZE observable, the ICM masses $M_{\mathrm{ICM}}$ are obtained from the analysis of $Chandra$ X-ray observations, and the stellar masses $M_{\star}$ are derived by fitting spectral energy distribution templates to Dark Energy Survey (DES) $griz$ optical photometry and $WISE$ or $Spitzer$ near-infrared photometry. We study trends in the stellar mass, the ICM mass, the total baryonic mass and the cold baryonic fraction with cluster mass and redshift. We find significant departures from self-similarity in the mass scaling for all quantities, while the redshift trends are all statistically consistent with zero, indicating that the baryon content of clusters at fixed mass has changed remarkably little over the past $\approx9$ Gyr. We compare our results to the mean baryon fraction (and the stellar mass fraction) in the field, finding that these values lie above (below) those in cluster virial regions in all but the most massive clusters at low redshift. Using a simple model of the matter assembly of clusters from infalling groups with lower masses and from infalling material from the low density environment or field surrounding the parent halos, we show that the measured mass trends without strong redshift trends in the stellar mass scaling relation could be explained by a mass and redshift dependent fractional contribution from field material. Similar analyses of the ICM and baryon mass scaling relations provide evidence for the so-called "missing baryons" outside cluster virial regions.

Strong gravitational lensing by galaxy clusters has become a powerful tool for probing the high-redshift Universe, magnifying distant and faint background galaxies. Reliable strong lensing (SL) models are crucial for determining the intrinsic properties of distant, magnified sources and for constructing their luminosity function. We present here the first SL analysis of MACS J0308.9+2645 and PLCK G171.9-40.7, two massive galaxy clusters imaged with the Hubble Space Telescope in the framework of the Reionization Lensing Cluster Survey (RELICS). We use the Light-Traces-Mass modeling technique to uncover sets of multiply imaged galaxies and constrain the mass distribution of the clusters. Our SL analysis reveals that both clusters have particularly large Einstein radii ($\theta_E>30"$ for a source redshift of $z_s=2$), providing fairly large areas with high magnifications, useful for high-redshift galaxy searches ($\sim2$ arcmin$^{2}$ with $\mu>5$ to $\sim1$ arcmin$^{2}$ with $\mu>10$, similar to a typical \textit{Hubble Frontier Fields} cluster). We also find that MACS J0308.9+2645 hosts a promising, apparently bright (J$\sim23.2-24.6$ AB), multiply imaged high-redshift candidate at $z\sim6.4$. These images are amongst the brightest high-redshift candidates found in RELICS. Our mass models, including magnification maps, are made publicly available for the community through the Mikulski Archive for Space Telescopes.

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.

The fraction of ionizing photons that escape high-redshift galaxies sensitively determines whether galaxies reionized the early universe. However, this escape fraction cannot be measured from high-redshift galaxies because the opacity of the intergalactic medium is large at high redshifts. Without methods to indirectly measure the escape fraction of high-redshift galaxies, it is unlikely that we will know what reionized the universe. Here, we analyze the far-ultraviolet (UV) H I (Lyman series) and low-ionization metal absorption lines of nine low-redshift, confirmed Lyman continuum emitting galaxies. We use the H I covering fractions, column densities, and dust attenuations measured in a companion paper to predict the escape fraction of ionizing photons. We find good agreement between the predicted and observed Lyman continuum escape fractions (within $1.4\sigma$) using both the H I and ISM absorption lines. The ionizing photons escape through holes in the H I, but we show that dust attenuation reduces the fraction of photons that escape galaxies. This means that the average high-redshift galaxy likely emits more ionizing photons than low-redshift galaxies. Two other indirect methods accurately predict the escape fractions: the Ly$\alpha$ escape fraction and the optical [O III]/[O II] flux ratio. We use these indirect methods to predict the escape fraction of a sample of 21 galaxies with rest-frame UV spectra but without Lyman continuum observations. Many of these galaxies have low escape fractions ($f_{\rm esc} \le 1$\%), but 11 have escape fractions $>1$\%. The methods presented here will measure the escape fractions of high-redshift galaxies, enabling future telescopes to determine whether star-forming galaxies reionized the early universe.

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.

Discovery of strongly-lensed gravitational wave (GW) sources will unveil binary compact objects at higher redshifts and lower intrinsic luminosities than is possible without lensing. Such systems will yield unprecedented constraints on the mass distribution in galaxy clusters, measurements of the polarization of GWs, tests of General Relativity, and constraints on the Hubble parameter. Excited by these prospects, and intrigued by the presence of so-called "heavy black holes" in the early detections by LIGO-Virgo, we commenced a search for strongly-lensed GWs and possible electromagnetic counterparts in the latter stages of the second LIGO observing run (O2). Here, we summarise our calculation of the detection rate of strongly-lensed GWs, describe our review of BBH detections from O1, outline our observing strategy in O2, summarize our follow-up observations of GW170814, and discuss the future prospects of detection.

Every star-forming galaxy has a halo of metal-enriched gas extending out to at least 100 kpc, as revealed by the absorption lines this gas imprints on the spectra of background quasars. However, quasars are sparse and typically probe only one narrow pencil beam through the intervening galaxy. Close quasar pairs and gravitationally lensed quasars have been used to circumvent this inherently one-dimensional technique, but these objects are rare and the structure of the circum-galactic medium remains poorly constrained. As a result, our understanding of the physical processes that drive the re-cycling of baryons across the lifetime of a galaxy is limited. Here we report integral-field (tomographic) spectroscopy of an extended background source -a bright giant gravitational arc. We can thus coherently map the spatial and kinematic distribution of Mg II absorption -a standard tracer of enriched gas- in an intervening galaxy system at redshift 0.98 (i.e., ~8 Gyr ago). Our gravitational-arc tomography unveils a clumpy medium in which the absorption-strength decreases with increasing impact parameter, in good agreement with the statistics towards quasars; furthermore, we find strong evidence that the gas is not distributed isotropically. Interestingly, we detect little kinematic variation over a projected area of ~600 kpc squared, with all line-of-sight velocities confined to within a few tens of km/s of each other. These results suggest that the detected absorption originates from entrained recycled material, rather than in a galactic outflow.

The most distant galaxies known are at z~10-11, observed 400-500 Myr after the Big Bang. The few z~10-11 candidates discovered to date have been exceptionally small- barely resolved, if at all, by the Hubble Space Telescope. Here we present the discovery of SPT0615-JD, a fortuitous z~10 (z_phot=9.9+/-0.6) galaxy candidate stretched into an arc over ~2.5" by the effects of strong gravitational lensing. Discovered in the Reionization Lensing Cluster Survey (RELICS) Hubble Treasury program and companion S-RELICS Spitzer program, this candidate has a lensed H-band magnitude of 25.7+/-0.1 AB mag. With a magnification of \mu~4-7 estimated from our lens models, the de-lensed intrinsic magnitude is 27.6+/-0.3 AB mag, and the half-light radius is r_e<0.8 kpc, both consistent with other z>9 candidates. The inferred stellar mass (log [M* /M_Sun]=9.7^{+0.7}_{-0.5}) and star formation rate (\log [SFR/M_Sun yr^{-1}]=1.3^{+0.2}_{-0.3}) indicate that this candidate is a typical star-forming galaxy on the z>6 SFR-M* relation. We note that three independent lens models predict two counterimages, at least one of which should be of a similar magnitude to the arc, but these counterimages are not yet detected. Counterimages would not be expected if the arc were at lower redshift. However, the only spectral energy distributions capable of fitting the Hubble and Spitzer photometry well at lower redshifts require unphysical combinations of z~2 galaxy properties. The unprecedented lensed size of this z~10 candidate offers the potential for the James Webb Space Telescope to study the geometric and kinematic properties of a galaxy observed 500 Myr after the Big Bang.

We introduce Project MEGaSaURA: The Magellan Evolution of Galaxies Spectroscopic and Ultraviolet Reference Atlas. MEGaSaURA comprises medium-resolution, rest-frame ultraviolet spectroscopy of N=15 bright gravitationally lensed galaxies at redshifts of 1.68$<$z$<$3.6, obtained with the MagE spectrograph on the Magellan telescopes. The spectra cover the observed-frame wavelength range $3200 < \lambda_o < 8280$ \AA ; the average spectral resolving power is R=3300. The median spectrum has a signal-to-noise ratio of $SNR=21$ per resolution element at 5000 \AA . As such, the MEGaSaURA spectra have superior signal-to-noise-ratio and wavelength coverage compared to what COS/HST provides for starburst galaxies in the local universe. This paper describes the sample, the observations, and the data reduction. We compare the measured redshifts for the stars, the ionized gas as traced by nebular lines, and the neutral gas as traced by absorption lines; we find the expected bulk outflow of the neutral gas, and no systemic offset between the redshifts measured from nebular lines and the redshifts measured from the stellar continuum. We provide the MEGaSaURA spectra to the astronomical community through a data release.

We stack the rest-frame ultraviolet spectra of N=14 highly magnified gravitationally lensed galaxies at redshifts 1.6<z<3.6. The resulting new composite spans $900< \lambda_{rest} < 3000$ \AA, with a peak signal-to-noise ratio of 103 per spectral resolution element ($\sim$100 km/s). It is the highest signal-to-noise ratio, highest spectral resolution composite spectrum of $z\sim2$--3 galaxies yet published. The composite reveals numerous weak nebular emission lines and stellar photospheric absorption lines that can serve as new physical diagnostics, particularly at high redshift with the James Webb Space Telescope (JWST). We report equivalent widths to aid in proposing for and interpreting JWST spectra. We examine the velocity profiles of strong absorption features in the composite, and in a matched composite of $z\sim0$ COS/HST galaxy spectra. We find remarkable similarity in the velocity profiles at $z\sim 0$ and $z\sim2$, suggesting that similar physical processes control the outflows across cosmic time. While the maximum outflow velocity depends strongly on ionization potential, the absorption-weighted mean velocity does not. As such, the bulk of the high-ionization absorption traces the low-ionization gas, with an additional blueshifted absorption tail extending to at least $-2000$ km/s . We interpret this tail as arising from the stellar wind and photospheres of massive stars. Starburst99 models are able to replicate this high-velocity absorption tail. However, these theoretical models poorly reproduce several of the photospheric absorption features, indicating that improvements are needed to match observational constraints on the massive stellar content of star-forming galaxies at $z \sim 2$. We publicly release our composite spectra.

We present Atacama Large Millimeter/submillimeter Array observations of CO lines and dust continuum emission of the source RCSGA 032727--132609, a young $z=1.7$ low-metallicity starburst galaxy. The CO(3-2) and CO(6-5) lines, and continuum at rest-frame $450\,\mu m$ are detected and show a resolved structure in the image plane. We use the corresponding lensing model to obtain a source plane reconstruction of the detected emissions revealing intrinsic flux density of $S_{450\,\mu m}=23.5_{-8.1}^{+26.8}$ $\mu$Jy and intrinsic CO luminosities $L'_{\rm CO(3-2)}=2.90_{-0.23}^{+0.21}\times10^{8}$ ${\rm K\,km\,s^{-1}\,pc^{2}}$ and $L'_{\rm CO(6-5)}=8.0_{-1.3}^{+1.4}\times10^{7}$ ${\rm K\,km\,s^{-1}\,pc^{2}}$. We used the resolved properties in the source plane to obtain molecular gas and star-formation rate surface densities of $\Sigma_{\rm H2}=16.2_{-3.5}^{+5.8}\,{\rm M}_{\odot}\,{\rm pc}^{-2}$ and $\Sigma_{\rm SFR}=0.54_{-0.27}^{+0.89}\,{\rm M}_{\odot}\,{\rm yr}^{-1}\,{\rm kpc}^{-2}$ respectively. The intrinsic properties of RCSGA 032727--132609 show an enhanced star-formation activity compared to local spiral galaxies with similar molecular gas densities, supporting the ongoing merger-starburst phase scenario. RCSGA 032727--132609 also appears to be a low--density starburst galaxy similar to local blue compact dwarf galaxies, which have been suggested as local analogs to high-redshift low-metallicity starburst systems. Finally, the CO excitation level in the galaxy is consistent with having the peak at ${\rm J}\sim5$, with a higher excitation concentrated in the star-forming clumps.

We report the detection of extended Lyman-$\alpha$ emission from the host galaxy of SDSS~J2222+2745, a strongly lensed quasar at $z = 2.8$. Spectroscopic follow-up clearly reveals extended Lyman-$\alpha$ in emission between two images of the central active galactic nucleus (AGN). We reconstruct the lensed quasar host galaxy in the source plane by applying a strong lens model to HST imaging, and resolve spatial scales as small as $\sim$200 parsecs. In the source plane we recover the host galaxy morphology to within a few hundred parsecs of the central AGN, and map the extended Lyman-$\alpha$ emission to its physical origin on one side of the host galaxy at radii $\sim$0.5-2 kpc from the central AGN. There are clear morphological differences between the Lyman-$\alpha$ and rest-frame ultraviolet stellar continuum emission from the quasar host galaxy. Furthermore, the relative velocity profiles of quasar Lyman-$\alpha$, host galaxy Lyman-$\alpha$, and metal lines in outflowing gas reveal differences in the absorbing material affecting the AGN and host galaxy. These data indicate the presence of patchy local intervening gas in front of the central quasar and its host galaxy. This interpretation is consistent with the central luminous quasar being obscured across a substantial fraction of its surrounding solid angle, resulting in strong anisotropy in the exposure of the host galaxy to ionizing radiation from the AGN. This work demonstrates the power of strong lensing-assisted studies to probe spatial scales that are currently inaccessible by other means.

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.

Using the combined resolving power of the Hubble Space Telescope and gravitational lensing, we resolve star-forming structures in a z~2.5 galaxy on scales much smaller than the usual kiloparsec diffraction limit of HST. SGAS J111020.0+645950.8 is a clumpy, star forming galaxy lensed by the galaxy cluster SDSS J1110+6459 at z = 0.659, with a total magnification ~30x across the entire arc. We use a hybrid parametric/non-parametric strong lensing mass model to compute the deflection and magnification of this giant arc, reconstruct the light distribution of the lensed galaxy in the source plane, and resolve the star formation into two dozen clumps. We develop a forward-modeling technique to model each clump in the source plane. We ray trace the model to the image plane, convolve with the instrumental point spread function (PSF), and compare with the GALFIT model of the clumps in the image plane, which decomposes clump structure from more extended emission. This technique has the advantage, over ray tracing, by accounting for the asymmetric lensing shear of the galaxy in the image plane and the instrument PSF. At this resolution, we can begin to study star formation on a clump-by-clump basis, toward the goal of understanding feedback mechanisms and the buildup of exponential disks at high redshift.

For lensed galaxy SGAS J111020.0+645950.8 at redshift z=2.481, which is magnified by a factor of 28 +- 8, we analyze the morphology of star formation as traced by rest-frame ultraviolet emission, in both the highly-magnified source plane, and in simulations of how this galaxy would appear without lensing magnification. Were this galaxy not lensed but drawn from an HST deep field, we would conclude that almost all its star formation arises from an exponential disk (S\'ersic index of 1.0 +- 0.4) with an effective radius of r_e = 2.7 +- 0.3 kpc measured from two-dimensional fitting to F606W using Galfit, and r_e=1.9 +- 0.1 kpc measured by fitting a radial profile to F606W elliptical isophotes. At the normal spatial resolution of the deep fields, there is no sign of clumpy star formation within SGAS J111020.0+645950.8 . However, the enhanced spatial resolution enabled by gravitational lensing tells a very different story: much of the star formation arises in two dozen clumps with sizes of r=30--50 pc spread across the 7 kpc length of the galaxy. The color and spatial distribution of the diffuse component suggests that still smaller clumps are unresolved. Despite this clumpy, messy morphology, the radial profile is still well-characterized by an exponential profile. In this lensed galaxy, stars are forming in complexes with sizes well below 100 pc; such sizes are wholly unexplored by surveys of galaxy evolution at 1<z<3.

We present measurements of the surface density of star formation, the star-forming clump luminosity function, and the clump size distribution function, for the lensed galaxy SGAS J111020.0+645950.8 at a redshift of z=2.481. The physical size scales that we probe, radii r=30--50 pc, are considerably smaller scales than have yet been studied at these redshifts. The star formation surface density we find within these small clumps is consistent with surface densities measured previously for other lensed galaxies at similar redshift. Twenty-two percent of the rest-frame ultraviolet light in this lensed galaxy arises from small clumps, with r<100 pc. Within the range of overlap, the clump luminosity function measured for this lensed galaxy is remarkably similar to those of z~0 galaxies. In this galaxy, star-forming regions smaller than 100 pc---physical scales not usually resolved at these redshifts by current telescopes---are important locations of star formation in the distant universe. If this galaxy is representative, this may contradict the theoretical picture in which the critical size scale for star formation in the distant universe is of order 1 kiloparsec. Instead, our results suggest that current telescopes have not yet resolved the critical size scales of star-forming activity in galaxies over most of cosmic time.

Until now, systematic errors in strong gravitational lens modeling have been acknowledged but never been fully quantified. Here, we launch an investigation into the systematics induced by constraint selection. We model the simulated cluster Ares 362 times using random selections of image systems with and without spectroscopic redshifts and quantify the systematics using several diagnostics: image predictability, accuracy of model-predicted redshifts, enclosed mass, and magnification. We find that for models with $>15$ image systems, the image plane rms does not decrease significantly when more systems are added; however the rms values quoted in the literature may be misleading as to the ability of a model to predict new multiple images. The mass is well constrained near the Einstein radius in all cases, and systematic error drops to $<2\%$ for models using $>10$ image systems. Magnification errors are smallest along the straight portions of the critical curve, and the value of the magnification is systematically lower near curved portions. For $>15$ systems, the systematic error on magnification is $\sim2\%$. We report no trend in magnification error with fraction of spectroscopic image systems when selecting constraints at random; however, when using the same selection of constraints, increasing this fraction up to $\sim0.5$ will increase model accuracy. The results suggest that the selection of constraints, rather than quantity alone, determines the accuracy of the magnification. We note that spectroscopic follow-up of at least a few image systems is crucial, as models without any spectroscopic redshifts are inaccurate across all of our diagnostics.

(Abridged) We combine deep HST grism spectroscopy with a new Bayesian method to derive maps of gas-phase metallicity, nebular dust extinction, and star-formation rate for 10 star-forming galaxies at high redshift ($1.2<z<2.3$). Exploiting lensing magnification by the foreground cluster MACS1149.6+2223, we reach sub-kpc spatial resolution and push the stellar mass limit associated with such high-z spatially resolved measurements below $10^8M_\odot$ for the first time. Our maps exhibit diverse morphologies, indicative of various effects such as efficient radial mixing from tidal torques, rapid accretion of low-metallicity gas, etc., which can affect the gas and metallicity distributions in individual galaxies. Based upon an exhaustive sample of all existing sub-kpc metallicity gradients at high-z, we find that predictions given by analytical chemical evolution models assuming a relatively extended star-formation profile in the early disk formation phase can explain the majority of observed gradients, without involving galactic feedback or radial outflows. We observe a tentative correlation between stellar mass and metallicity gradient, consistent with the downsizing galaxy formation picture that more massive galaxies are more evolved into a later phase of disk growth, where they experience more coherent mass assembly at all radii and thus show shallower metallicity gradients. In addition, we compile a sample of homogeneously cross-calibrated integrated metallicity measurements spanning three orders of magnitude in stellar mass at $z\sim1.8$. We use this sample to study the mass-metallicity relation (MZR) and test the fundamental metallicity relation (FMR). The slope of the observed MZR can rule out the momentum-driven wind model at 3-$\sigma$ confidence level. We find no significant offset with respect to the FMR, taking into account the intrinsic scatter and measurement uncertainties.

SDSS J2222+2745 is a galaxy cluster at z=0.49, strongly lensing a quasar at z=2.805 into six widely separated images. In recent HST imaging of the field, we identify additional multiply lensed galaxies, and confirm the sixth quasar image that was identified by Dahle et al. (2013). We used the Gemini North telescope to measure a spectroscopic redshift of z=4.56 of one of the secondary lensed galaxies. These data are used to refine the lens model of SDSS J2222+2745, compute the time delay and magnifications of the lensed quasar images, and reconstruct the source image of the quasar host and a second lensed galaxy at z=2.3. This second galaxy also appears in absorption in our Gemini spectra of the lensed quasar, at a projected distance of 34 kpc. Our model is in agreement with the recent time delay measurements of Dahle et al. (2015), who found tAB=47.7+/-6.0 days and tAC=-722+/-24 days. We use the observed time delays to further constrain the model, and find that the model-predicted time delays of the three faint images of the quasar are tAD=502+/-68 days, tAE=611+/-75 days, and tAF=415+/-72 days. We have initiated a follow-up campaign to measure these time delays with Gemini North. Finally, we present initial results from an X-ray monitoring program with Swift, indicating the presence of hard X-ray emission from the lensed quasar, as well as extended X-ray emission from the cluster itself, which is consistent with the lensing mass measurement and the cluster velocity dispersion.

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.

We probe the spatial distribution of outflowing gas along four lines of sight separated by up to 6 kpc in a gravitationally-lensed star-forming galaxy at z=1.70. Using MgII and FeII emission and absorption as tracers, we find that the clumps of star formation are driving galactic outflows with velocities of -170 to -250 km/sec. The velocities of MgII emission are redshifted with respect to the systemic velocities of the galaxy, consistent with being back-scattered. By contrast, the FeII fluorescent emission lines are either slightly blueshifted or at the systemic velocity of the galaxy. Taken together, the velocity structure of the MgII and FeII emission is consistent with arising through scattering in galactic winds. Assuming a thin shell geometry for the out owing gas, the estimated masses carried out by these outfows are large (> 30 - 50 $\rm{M_{\odot} yr^{-1}}$), with mass loading factors several times the star-formation rate. Almost 20% to 50% of the blueshifted absorption probably escapes the gravitational potential of the galaxy. In this galaxy, the outflow is "locally sourced", that is, the properties of the outflow in each line of sight are dominated by the properties of the nearest clump of star formation; the wind is not global to the galaxy. The mass outflow rates and the momentum flux carried out by outflows in individual star forming knots of this object are comparable to that of starburst galaxies in the local Universe.

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.

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.

The sample of cosmological strong lensing systems has been steadily growing in recent years and with the advent of the next generation of space-based survey telescopes, the sample will reach into the thousands. The accuracy of strong lens models relies on robust identification of multiple image families of lensed galaxies. For the most massive lenses, often more than one background galaxy is magnified and multiply-imaged, and even in the cases of only a single lensed source, identification of counter images is not always robust. Recently, we have shown that the Gini coefficient in space-telescope-quality imaging is a measurement of galaxy morphology that is relatively well-preserved by strong gravitational lensing. Here, we investigate its usefulness as a diagnostic for the purposes of image family identification and show that it can remove some of the degeneracies encountered when using color as the sole diagnostic, and can do so without the need for additional observations since whenever a color is available, two Gini coefficients are as well.

We present an initial study of the mass and evolutionary state of a massive and distant cluster, RCS2 J232727.6-020437. This cluster, at z=0.6986, is the richest cluster discovered in the RCS2 project. The mass measurements presented in this paper are derived from all possible mass proxies: X-ray measurements, weak-lensing shear, strong lensing, Sunyaev Zel'dovich effect decrement, the velocity distribution of cluster member galaxies, and galaxy richness. While each of these observables probe the mass of the cluster at a different radius, they all indicate that RCS2 J232727.6-020437 is among the most massive clusters at this redshift, with an estimated mass of M_200 ~3 x10^15 h^-1 Msun. In this paper, we demonstrate that the various observables are all reasonably consistent with each other to within their uncertainties. RCS2 J232727.6-020437 appears to be well relaxed -- with circular and concentric X-ray isophotes, with a cool core, and no indication of significant substructure in extensive galaxy velocity data.

We measure C III] 1907,1909 A emission lines in eleven gravitationally--lensed star-forming galaxies at z~1.6--3, finding much lower equivalent widths than previously reported for fainter lensed galaxies (Stark et al. 2014). While it is not yet clear what causes some galaxies to be strong C III] emitters, CIII] emission is not a universal property of distant star-forming galaxies. We also examine C III] emission in 46 star-forming galaxies in the local universe, using archival spectra from GHRS, FOS, and STIS on HST, and IUE. Twenty percent of these local galaxies show strong C III] emission, with equivalent widths <-5 A. Three nearby galaxies show C III] emission equivalent widths as large as the most extreme emitters yet observed in the distant universe; all three are Wolf-Rayet galaxies. At all redshifts, strong C III] emission may pick out low-metallicity galaxies experiencing intense bursts of star formation. Such local C III] emitters may shed light on the conditions of star formation in certain extreme high-redshift galaxies.

We report first results from an ongoing monitoring campaign to measure time delays between the six images of the quasar SDSS\,J2222$+$2745, gravitationally lensed by a galaxy cluster. The time delay between A and B, the two most highly magnified images, is measured to be $\tau_{\rm AB} = 47.7 \pm 6.0$ days (95\% confidence interval), consistent with previous model predictions for this lens system. The strong intrinsic variability of the quasar also allows us to derive a time delay value of $\tau_{\rm CA} = 722 \pm 24$ days between image C and A, in spite of modest overlap between their light curves in the current data set. Image C, which is predicted to lead all the other lensed quasar images, has undergone a sharp, monotonic flux increase of 60-75\% during 2014. A corresponding brightening is firmly predicted to occur in images A and B during 2016. The amplitude of this rise indicates that time delays involving all six known images in this system, including those of the demagnified central images D-F, will be obtainable from further ground-based monitoring of this system during the next few years.

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.

We test the effects of varying the cosmological parameter values used in the strong lens modeling process for the six Hubble Frontier Fields (HFF) galaxy clusters. The standard procedure for generating high fidelity strong lens models includes careful consideration of uncertainties in the output models that result from varying model parameters within the bounds of available data constraints. It is not, however, common practice to account for the effects of cosmological parameter value uncertainties. The convention is to instead use a single fiducial "concordance cosmology" and generate lens models assuming zero uncertainty in cosmological parameter values. We find that the magnification maps of the individual HFF clusters vary significantly when lens models are computed using different cosmological parameter values taken from recent literature constraints from space- and ground-based experiments. Specifically, the magnification maps have average variances across the best fit models computed using different cosmologies that are comparable in magnitude to - and as much as 2.5 times larger than - the model fitting uncertainties in each best fit model. We also find that estimates of the mass profiles of the cluster cores themselves vary only slightly when different input cosmological parameters are used. We conclude that cosmological parameter uncertainty is a non-negligible source of uncertainty in lens model products for the HFF clusters, and that it is important that current and future work which relies on precision strong lensing models take care to account for this additional source of uncertainty.

We present a detailed gravitational lens model of the galaxy cluster RCS2 J232727.6-020437. Due to cosmological dimming of cluster members and ICL, its high redshift ($z=0.6986$) makes it ideal for studying background galaxies. Using new ACS and WFC3/IR HST data, we identify 16 multiple images. From MOSFIRE follow up, we identify a strong emission line in the spectrum of one multiple image, likely confirming the redshift of that system to $z=2.083$. With a highly magnified ($\mu\gtrsim2$) source plane area of $\sim0.7$ arcmin$^2$ at $z=7$, RCS2 J232727.6-020437 has a lensing efficiency comparable to the Hubble Frontier Fields clusters. We discover four highly magnified $z\sim7$ candidate Lyman-break galaxies behind the cluster, one of which may be multiply-imaged. Correcting for magnification, we find that all four galaxies are fainter than $0.5 L_{\star}$. One candidate is detected at ${>10\sigma}$ in both Spitzer/IRAC [3.6] and [4.5] channels. A spectroscopic follow-up with MOSFIRE does not result in the detection of the Lyman-alpha emission line from any of the four candidates. From the MOSFIRE spectra we place median upper limits on the Lyman-alpha flux of $5-14 \times 10^{-19}\, \mathrm{erg \,\, s^{-1} cm^{-2}}$ ($5\sigma$).

We present a revised lens model of MACS J1149+2223, in which the first resolved multiply-imaged lensed supernova was discovered. The lens model is based on the model of Johnson et al. (2014) with some modifications. We include more lensing constraints from the host galaxy of the newly discovered supernova, and increase the flexibility of the model in order to better reproduce the lensing signal in the vicinity of this galaxy. The revised model accurately reconstructs the positions of the lensed supernova, provides magnifications, and predicts the time delay between the instances of the supernova. Finally, we reconstruct the source image of the host galaxy, and position the supernova on one of its spiral arms. Products of this lens model are available to the community through MAST.

We present strong-lensing models, as well as mass and magnification maps, for the cores of the six HST Frontier Fields galaxy clusters. Our parametric lens models are constrained by the locations and redshifts of multiple image systems of lensed background galaxies. We use a combination of photometric redshifts and spectroscopic redshifts of the lensed background sources obtained by us (for Abell 2744 and Abell S1063), collected from the literature, or kindly provided by the lensing community. Using our results, we (1) compare the derived mass distribution of each cluster to its light distribution, (2) quantify the cumulative magnification power of the HFF clusters, (3) describe how our models can be used to estimate the magnification and image multiplicity of lensed background sources at all redshifts and at any position within the cluster cores, and (4) discuss systematic effects and caveats resulting from our modeling methods. We specifically investigate the effect of the use of spectroscopic and photometric redshift constraints on the uncertainties of the resulting models. We find that the photometric redshift estimates of lensed galaxies are generally in excellent agreement with spectroscopic redshifts, where available. However, the flexibility associated with relaxed redshift priors may cause the complexity of large-scale structure that is needed to account for the lensing signal to be underestimated. Our findings thus underline the importance of spectroscopic arc redshifts, or tight photometric redshift constraints, for high precision lens models. All products from our best-fit lens models (magnification, convergence, shear, deflection field) and model simulations for estimating errors are made available via the Mikulski Archive for Space Telescopes.

Star formation occurs on physical scales corresponding to individual star forming regions, typically of order ~100 parsecs in size, but current observational facilities cannot resolve these scales within field galaxies beyond the local universe. However, the magnification from strong gravitational lensing allows us to measure the properties of these discrete star forming regions within galaxies in the distant universe. New results from multi-wavelength spectroscopic studies of a sample of extremely bright, highly magnified lensed galaxies are revealing the complexity of star formation on sub-galaxy scales during the era of peak star formation in the universe. We find a wide range of properties in the rest-frame UV spectra of individual galaxies, as well as in spectra that originate from different star forming regions within the same galaxy. Large variations in the strengths and velocity structure of Lyman-alpha and strong P Cygni lines such as C IV, and MgII provide new insights into the astrophysical relationships between extremely massive stars, the elemental abundances and physical properties of the nebular gas those stars ionize, and the galactic-scale outflows they power.

We present the mass distribution at the core of SDSS J1531+3414, a strong-lensing cluster at z=0.335. We find that the mass distribution is well described by two cluster-scale halos with a contribution from cluster-member galaxies. New HST observations of SDSS J1531+3414 reveal a signature of ongoing star formation associated with the two central galaxies at the core of the cluster, in the form of a chain of star forming regions at the center of the cluster. Using the lens model presented here, we place upper limits on the contribution of a possible lensed image to the flux at the center region, and rule out that this emission is coming from a background source.

New Hubble Space Telescope ultraviolet and optical imaging of the strong-lensing galaxy cluster SDSS J1531+3414 (z=0.335) reveals two centrally dominant elliptical galaxies participating in an ongoing major merger. The interaction is at least somewhat rich in cool gas, as the merger is associated with a complex network of nineteen massive superclusters of young stars (or small tidal dwarf galaxies) separated by ~1 kpc in projection from one another, combining to an estimated total star formation rate of ~5 solar masses per year. The resolved young stellar superclusters are threaded by narrow H-alpha, [O II], and blue excess filaments arranged in a network spanning ~27 kpc across the two merging galaxies. This morphology is strongly reminiscent of the well-known "beads on a string" mode of star formation observed on kpc-scales in the arms of spiral galaxies, resonance rings, and in tidal tails between interacting galaxies. Nevertheless, the arrangement of this star formation relative to the nuclei of the two galaxies is difficult to interpret in a dynamical sense, as no known "beads on a string" systems associated with kpc-scale tidal interactions exhibit such lopsided morphology relative to the merger participants. In this Letter we present the images and follow-up spectroscopy, and discuss possible physical interpretations for the unique arrangement of the young stellar clusters. While we suggest that this morphology is likely to be dynamically short-lived, a more quantitative understanding awaits necessary multiwavelength follow-up, including optical integral field spectroscopy, ALMA sub-mm interferometry, and Chandra X-ray imaging.

We present optical and near-IR imaging and spectroscopy of SGAS J105039.6$+$001730, a strongly lensed galaxy at z $=$ 3.6252 magnified by $>$30$\times$, and derive its physical properties. We measure a stellar mass of log(M$_{*}$/M$_{\odot}$) $=$ 9.5 $\pm$ 0.35, star formation rates from [O II]$\lambda$$\lambda$3727 and H-$\beta$ of 55 $\pm$ 20 and 84 $\pm$ 17 M$_{\odot}$ yr$^{-1}$, respectively, an electron density of n$_{e} \leq$ 10$^{3}$ cm$^{-2}$, an electron temperature of T$_{e} \leq$ 14000 K, and a metallicity of 12+log(O/H) $=$ 8.3 $\pm$ 0.1. The strong C III]$\lambda$$\lambda$1907,1909 emission and abundance ratios of C, N, O and Si are consistent with well-studied starbursts at z $\sim$ 0 with similar metallicities. Strong P Cygni lines and He II$\lambda$1640 emission indicate a significant population of Wolf-Rayet stars, but synthetic spectra of individual populations of young, hot stars do not reproduce the observed integrated P Cygni absorption features. The rest-frame UV spectral features are indicative of a young starburst with high ionization, implying either 1) an ionization parameter significantly higher than suggest by rest-frame optical nebular lines, or 2) differences in one or both of the initial mass function and the properties of ionizing spectra of massive stars. We argue that the observed features are likely the result of a superposition of star forming regions with different physical properties. These results demonstrate the complexity of star formation on scales smaller than individual galaxies, and highlight the importance of systematic effects that result from smearing together the signatures of individual star forming regions within galaxies.

We present a detailed analysis of Hubble Space Telescope (HST), Wide Field Camera 3 (WFC3) G141 grism spectroscopy for seven star-forming regions of the highly magnified lensed starburst galaxy RCSGA 032727-132609 at z=1.704. We measure the spatial variations of the extinction in RCS0327 through the observed H$\gamma$/H$\beta$ emission line ratios, finding a constant average extinction of $\mathrm{E(B-V)_{gas}}=0.40\pm0.07$. We infer that the star formation is enhanced as a result of an ongoing interaction, with measured star formation rates derived from demagnified, extinction-corrected H$\beta$ line fluxes for the individual star-forming clumps falling >1-2 dex above the star formation sequence. When combining the HST/WFC3 [OIII]$\lambda$5007/H$\beta$ emission line ratio measurements with [NII]/H$\alpha$ line ratios from Wuyts et a. (2014), we find that the majority of the individual star-forming regions fall along the local "normal" abundance sequence. With the first detections of the He I $\lambda$5876 and He II $\lambda$4686 recombination lines in a distant galaxy, we probe the massive-star content of the star-forming regions in RCS0327. The majority of the star-forming regions have a He I $\lambda$5876 to H$\beta$ ratio consistent with the saturated maximum value, which is only possible if they still contain hot O-stars. Two regions have lower ratios, implying that their last burst of new star formation ended ~5 Myr ago. Together, the He I $\lambda$5876 and He II $\lambda$4686 to H$\beta$ line ratios provide indirect evidence for the order in which star formation is stopping in individual star-forming knots of this high redshift merger. We place the spatial variations of the extinction, star formation rate and ionization conditions in the context of the star formation history of RCS0327.

We examine the Mg II 2796, 2803 Angstrom, Lyman alpha, and nebular line emission in five bright star-forming galaxies at 1.66<z<1.91 that have been gravitationally lensed by foreground galaxy clusters. All five galaxies show prominent Mg II emission and absorption in a P Cygni profile. We find no correlation between the equivalent widths of Mg II and Lyman alpha emission. The Mg II emission has a broader range of velocities than do the nebular emission line profiles; the Mg II emission is redshifted with respect to systemic by 100 to 200 km/s. When present, Lyman alpha is even more redshifted. The reddest components of Mg II and Lyman alpha emission have tails to 500-600 km/s, implying a strong outflow. The lack of correlation in the Mg II and Lyman alpha equivalent widths, the differing velocity profiles, and the high ratios of Mg II to nebular line fluxes together suggest that the bulk of Mg II emission does not ultimately arise as nebular line emission, but may instead be reprocessed stellar continuum emission.

We present an analysis of the line-of-sight structure toward a sample of ten strong lensing cluster cores. Structure is traced by groups that are identified spectroscopically in the redshift range, 0.1 $\leq$ z $\leq$ 0.9, and we measure the projected angular and comoving separations between each group and the primary strong lensing clusters in each corresponding line of sight. From these data we measure the distribution of projected angular separations between the primary strong lensing clusters and uncorrelated large scale structure as traced by groups. We then compare the observed distribution of angular separations for our strong lensing selected lines of sight against the distribution of groups that is predicted for clusters lying along random lines of sight. There is clear evidence for an excess of structure along the line of sight at small angular separations ($\theta \leq 6'$) along the strong lensing selected lines of sight, indicating that uncorrelated structure is a significant systematic that contributes to producing galaxy clusters with large cross sections for strong lensing. The prevalence of line-of-sight structure is one of several biases in strong lensing clusters that can potentially be folded into cosmological measurements using galaxy cluster samples. These results also have implications for current and future studies -- such as the Hubble Space Telescope Frontier Fields -- that make use of massive galaxy cluster lenses as precision cosmological telescopes; it is essential that the contribution of line-of-sight structure be carefully accounted for in the strong lens modeling of the cluster lenses.

We present a detailed analysis of multi-wavelength HST/WFC3 imaging and Keck/OSIRIS near-IR AO-assisted integral field spectroscopy for a highly magnified lensed galaxy at z=1.70. This young starburst is representative of UV-selected star-forming galaxies (SFG) at z~2 and contains multiple individual star-forming regions. Due to the lensing magnification, we can resolve spatial scales down to 100pc in the source plane of the galaxy. The velocity field shows disturbed kinematics suggestive of an ongoing interaction, and there is a clear signature of a tidal tail. We constrain the age, reddening, SFR and stellar mass of the star-forming clumps from SED modelling of the WFC3 photometry and measure their H-alpha luminosity, metallicity and outflow properties from the OSIRIS data. With strong star formation driven outflows in four clumps, RCSGA0327 is the first high redshift SFG at stellar mass <10^10 M_sun with spatially resolved stellar winds. We compare the H-alpha luminosities, sizes and dispersions of the star-forming regions to other high-z clumps as well as local giant HII regions and find no evidence for increased clump star formation surface densities in interacting systems, unlike in the local Universe. Spatially resolved SED modelling unveils an established stellar population at the location of the largest clump and a second mass concentration near the edge of the system which is not detected in H-alpha emission. This suggests a picture of an equal-mass mixed major merger, which has not triggered a new burst of star formation or caused a tidal tail in the gas-poor component.

Absorption-line spectroscopy of multiply-lensed QSOs near a known foreground galaxy provides a unique opportunity to go beyond the traditional one-dimensional application of QSO probes and establish a crude three-dimensional (3D) map of halo gas around the galaxy that records the line-of-sight velocity field at different locations in the gaseous halo. Two intermediate-redshift galaxies are targeted in the field around the quadruply-lensed QSO HE0435-1223 at redshift z=1.689, and absorption spectroscopy along each of the lensed QSOs is carried out in the vicinities of these galaxies. One galaxy is a typical, star-forming L* galaxy at z=0.4188 and projected distance of rho=50 kpc from the lensing galaxy. The other is a super-L* barred spiral at z=0.7818 and rho=33 kpc. Combining known orientations of the quadruply-lensed QSO to the two foreground galaxies with the observed MgII absorption profiles along individual QSO sightlines has for the first time led to spatially resolved kinematics of tenuous halo gas on scales of 5-10 kpc at z>0.2. A MgII absorber is detected in every sightline observed through the halos of the two galaxies, and the recorded absorber strength is typical of what is seen in previous close QSO--galaxy pair studies. While the multi-sightline study confirms the unity covering fraction of MgII absorbing gas at rho < 50 kpc from star-forming disks, the galaxies also present two contrasting examples of complex halo gas kinematics. Different models, including a rotating disk, collimated outflows, and gaseous streams from either accretion or tidal/ram-pressure stripping, are considered for comparisons with the absorption-line observations, and infalling streams/stripped gas of width >~ 10 kpc are found to best describe the observed gas kinematics across multiple sightlines.

The process by which the mass density profile of certain galaxy clusters becomes centrally concentrated enough to produce high strong lensing (SL) cross-sections is not well understood. It has been suggested that the baryonic condensation of the intra-cluster medium (ICM) due to cooling may drag dark matter to the cores and thus steepen the profile. In this work, we search for evidence of ongoing ICM cooling in the first large, well-defined sample of strong lensing selected galaxy clusters in the range 0.1 < z < 0.6. Based on known correlations between the ICM cooling rate and both optical emission line luminosity and star formation, we measure, for a sample of 89 strong lensing clusters, the fraction of clusters that have [OII]3727 emission in their brightest cluster galaxy (BCG). We find that the fraction of line-emitting BCGs is constant as a function of redshift for z > 0.2 and shows no statistically significant deviation from the total cluster population. Specific star formation rates, as traced by the strength of the 4000 angstrom break, D_4000, are also consistent with the general cluster population. Finally, we use optical imaging of the SL clusters to measure the angular separation, R_arc, between the arc and the center of mass of each lensing cluster in our sample and test for evidence of changing [OII] emission and D_4000 as a function of R_arc, a proxy observable for SL cross-sections. D_4000 is constant with all values of R_arc, and the [OII] emission fractions show no dependence on R_arc for R_arc > 10" and only very marginal evidence of increased weak [OII] emission for systems with R_arc < 10". These results argue against the ability of baryonic cooling associated with cool core activity in the cores of galaxy clusters to strongly modify the underlying dark matter potential, leading to an increase in strong lensing cross-sections.

We report the detection of CO(1-0) line emission in the bright, lensed star-forming galaxies SPT-S 233227-5358.5 (z=2.73) and SPT-S 053816-5030.8 (z=2.78), using the Australia Telescope Compact Array (ATCA). Both galaxies were discovered in a large-area millimeter survey with the South Pole Telescope (SPT) and found to be gravitationally lensed by intervening structures. The measured CO intensities imply galaxies with molecular gas masses of (3.2 \pm 0.5)x10^10 (mu/15)^{-1}(X_CO/0.8) M_sun and (1.7 \pm 0.3)x10^10 (mu/20)^{-1}(X_CO/0.8) M_sun, and gas depletion timescales of 4.9x10^7 (X_CO/0.8) yr and 2.6x10^7 (X_CO/0.8) yr, respectively, where mu corresponds to the lens magnification and X_CO is the CO luminosity to gas mass conversion factor. In the case of SPT-S 053816-5030.8, we also obtained significant detections of the rest-frame 115.7 and 132.4 GHz radio continuum. Based on the radio to infrared spectral energy distribution and an assumed synchrotron spectral index, we find that 42 \pm 10 % and 55 \pm 13 % of the flux at rest-frame 115.7 and 132.4 GHz arise from free-free emission. We find a radio-derived intrinsic star formation rate (SFR) of 470 \pm 170 M_sun/yr, consistent within the uncertainties with the infrared estimate. Based on the morphology of this object in the source plane, the derived gas mass and the possible flattening of the radio spectral index towards low frequencies, we argue that SPT-S 053816-5030.8 exhibits properties compatible with a scaled-up local ultra-luminous infrared galaxy.

We present Atacama Large Millimeter/submillimeter Array (ALMA) 860 micrometer imaging of four high-redshift (z=2.8-5.7) dusty sources that were detected using the South Pole Telescope (SPT) at 1.4 mm and are not seen in existing radio to far-infrared catalogs. At 1.5 arcsec resolution, the ALMA data reveal multiple images of each submillimeter source, separated by 1-3 arcsec, consistent with strong lensing by intervening galaxies visible in near-IR imaging of these sources. We describe a gravitational lens modeling procedure that operates on the measured visibilities and incorporates self-calibration-like antenna phase corrections as part of the model optimization, which we use to interpret the source structure. Lens models indicate that SPT0346-52, located at z=5.7, is one of the most luminous and intensely star-forming sources in the universe with a lensing corrected FIR luminosity of 3.7 X 10^13 L_sun and star formation surface density of 4200 M_sun yr^-1 kpc^-2. We find magnification factors of 5 to 22, with lens Einstein radii of 1.1-2.0 arcsec and Einstein enclosed masses of 1.6-7.2x10^11 M_sun. These observations confirm the lensing origin of these objects, allow us to measure the their intrinsic sizes and luminosities, and demonstrate the important role that ALMA will play in the interpretation of lensed submillimeter sources.

In the past decade, our understanding of galaxy evolution has been revolutionized by the discovery that luminous, dusty, starburst galaxies were 1,000 times more abundant in the early Universe than at present. It has, however, been difficult to measure the complete redshift 2 distribution of these objects, especially at the highest redshifts (z > 4). Here we report a redshift survey at a wavelength of three millimeters, targeting carbon monoxide line emission from the star-forming molecular gas in the direction of extraordinarily bright millimetrewave-selected sources. High-resolution imaging demonstrates that these sources are strongly gravitationally lensed by foreground galaxies. We detect spectral lines in 23 out of 26 sources and multiple lines in 12 of those 23 sources, from which we obtain robust, unambiguous redshifts. At least 10 of the sources are found to lie at z > 4, indicating that the fraction of dusty starburst galaxies at high redshifts is greater than previously thought. Models of lens geometries in the sample indicate that the background objects are ultra-luminous infrared galaxies, powered by extreme bursts of star formation.

We present the discovery and a detailed multi-wavelength study of a strongly-lensed luminous infrared galaxy at z=0.816. Unlike most known lensed galaxies discovered at optical or near-infrared wavelengths this lensed source is red, r-Ks = 3.9 [AB], which the data presented here demonstrate is due to ongoing dusty star formation. The overall lensing magnification (a factor of 17) facilitates observations from the blue optical through to 500micron, fully capturing both the stellar photospheric emission as well as the re-processed thermal dust emission. We also present optical and near-IR spectroscopy. These extensive data show that this lensed galaxy is in many ways typical of IR-detected sources at z~1, with both a total luminosity and size in accordance with other (albeit much less detailed) measurements in samples of galaxies observed in deep fields with the Spitzer telescope. Its far-infrared spectral energy distribution is well-fit by local templates that are an order of magnitude less luminous than the lensed galaxy; local templates of comparable luminosity are too hot to fit. Its size (D~7kpc) is much larger than local luminous infrared galaxies, but in line with sizes observed for such galaxies at z~1. The star formation appears uniform across this spatial scale. In this source, the luminosity of which is typical of sources that dominate the cosmic infrared background, we find that star formation is spatially extended and well organised, quite unlike the compact merger-driven starbursts which are typical for sources of this luminosity at z~0.

We report the discovery of a unique gravitational lens system, SDSSJ2222+2745, producing five spectroscopically confirmed images of a z_s=2.82 quasar lensed by a foreground galaxy cluster at z_l=0.49. We also present photometric and spectroscopic evidence for a sixth lensed image of the same quasar. The maximum separation between the quasar images is 15.1". Both the large image separations and the high image multiplicity of the lensed quasar are in themselves exceptionally rare, and observing the combination of these two factors is an exceptionally unlikely occurrence in present datasets. This is only the third known case of a quasar lensed by a cluster, and the only one with six images. The lens system was discovered in the course of the Sloan Giant Arcs Survey, in which we identify candidate lenses in the Sloan Digital Sky Survey and target these for follow up and verification with the 2.56m Nordic Optical Telescope. Multi-band photometry obtained over multiple epochs from September 2011 to September 2012 reveal significant variability at the ~10-30% level in some of the quasar images, indicating that measurements of the relative time delay between quasar images will be feasible. In this lens system we also identify a bright (g = 21.5) giant arc corresponding to a strongly lensed background galaxy at z_s=2.30. We fit parametric models of the lens system, constrained by the redshift and positions of the quasar images and the redshift and position of the giant arc. The predicted time delays between different pairs of quasar images range from ~100 days to ~6 years.

We present APEX SABOCA 350micron and LABOCA 870micron observations of 11 representative examples of the rare, extremely bright (S_1.4mm > 15mJy), dust-dominated millimeter-selected galaxies recently discovered by the South Pole Telescope (SPT). All 11 sources are robustly detected with LABOCA with 40 < S_870micron < 130mJy, approximately an order of magnitude higher than the canonical submillimeter galaxy (SMG) population. Six of the sources are also detected by SABOCA at >3sigma, with the detections or upper limits providing a key constraint on the shape of the spectral energy distribution (SED) near its peak. We model the SEDs of these galaxies using a simple modified blackbody and perform the same analysis on samples of SMGs of known redshift from the literature. These calibration samples inform the distribution of dust temperature for similar SMG populations, and this dust temperature prior allows us to derive photometric redshift estimates and far infrared luminosities for the sources. We find a median redshift of <z> = 3.0, higher than the <z> = 2.2 inferred for the normal SMG population. We also derive the apparent size of the sources from the temperature and apparent luminosity, finding them to appear larger than our unlensed calibration sample, which supports the idea that these sources are gravitationally magnified by massive structures along the line of sight.

We take advantage of gravitational lensing amplification by Abell 1689 (z=0.187) to undertake the first space-based census of emission line galaxies (ELGs) in the field of a massive lensing cluster. Forty-three ELGs are identified to a flux of i_775=27.3 via slitless grism spectroscopy. One ELG (at z=0.7895) is very bright owing to lensing magnification by a factor of ~4.5. Several Balmer emission lines detected from ground-based follow-up spectroscopy signal the onset of a major starburst for this low-mass galaxy (M_* = 2 x 10^9 solar masses) with a high specific star formation rate (~20 /Gyr). From the blue emission lines we measure a gas-phase oxygen abundance consistent with solar (12+log(O/H)=8.8 +/- 0.2). We break the continuous line-emitting region of this giant arc into seven ~1kpc bins (intrinsic size) and measure a variety of metallicity dependent line ratios. A weak trend of increasing metal fraction is seen toward the dynamical center of the galaxy. Interestingly, the metal line ratios in a region offset from the center by ~1kpc have a placement on the blue HII region excitation diagram with f([OIII])/f(Hbeta) and f([NeIII])/f(Hbeta) that can be fit by an AGN. This asymmetrical AGN-like behavior is interpreted as a product of shocks in the direction of the galaxy's extended tail, possibly instigated by a recent galaxy interaction.

We present multi-wavelength imaging and near-IR spectroscopy for ten gravitationally lensed galaxies at 0.9<z<2.5 selected from a new, large sample of strong lens systems in the Sloan Digital Sky Survey (SDSS) DR7. We derive stellar masses from the rest-frame UV to near-IR spectral energy distributions, star formation rates (SFR) from the dust-corrected Ha flux, and metallicities from the [N II]/Ha flux ratio. We combine the lensed galaxies with a sample of sixty star-forming galaxies from the literature in the same redshift range for which measurements of [N II]/Ha have been published. Due to the lensing magnification, the lensed galaxies probe intrinsic stellar masses that are on average a factor of 11 lower than have been studied so far at these redshifts. They have specific star formation rates that are an order of magnitude higher than seen for main-sequence star-forming galaxies at z~2. We measure an evolution of 0.16+/-0.06 dex in the mass-metallicity relation between z~1.4 and z~2.2. In contrast to previous claims, the redshift evolution is smaller at low stellar masses. We do not see a correlation between metallicity and SFR at fixed stellar mass. The combined sample is in general agreement with the local fundamental relation between metallicity, stellar mass and SFR from Mannucci et al. (2010, 2011). Using the Kennicutt-Schmidt law to infer gas fractions, we investigate the importance of gas inflows and outflows on the shape of the mass-metallicity relation using simple analytical models. This suggests that the Maiolino et al.(2008) calibration of the [N II]/Ha flux ratio is biased high. We conclude that both an absolute metallicity calibration and direct measurements of the gas mass are needed to use the observed mass-metallicity relation to gain insight into the impact of gas flows on the chemical evolution of galaxies.

We present a comprehensive analysis of the rest-frame UV to near-IR spectral energy distributions and rest-frame optical spectra of four of the brightest gravitationally lensed galaxies in the literature: RCSGA 032727-132609 at z=1.70, MS1512-cB58 at z=2.73, SGAS J152745.1+065219 at z=2.76 and SGAS J122651.3+215220 at z=2.92. This includes new Spitzer imaging for RCSGA0327 as well as new spectra, near-IR imaging and Spitzer imaging for SGAS1527 and SGAS1226. Lensing magnifications of 3-4 magnitudes allow a detailed study of the stellar populations and physical conditions. We compare star formation rates as measured from the SED fit, the H-alpha and [OII] emission lines, and the UV+IR bolometric luminosity where 24 micron photometry is available. The SFR estimate from the SED fit is consistently higher than the other indicators, which suggests that the Calzetti dust extinction law used in the SED fitting is too flat for young star-forming galaxies at z~2. Our analysis finds similar stellar population parameters for all four lensed galaxies: stellar masses 3-7*10^9 M_sun, young ages ~ 100 Myr, little dust content E(B-V)=0.10-0.25, and star formation rates around 20-100 M_sun/yr. Compared to typical values for the galaxy population at z~2, this suggests we are looking at newly formed, starbursting systems that have only recently started the build-up of stellar mass. These results constitute the first detailed, uniform analysis of a sample of the growing number of strongly lensed galaxies known at z~2.

We present new HST/WFC3 imaging data of RCSGA 032727-132609, a bright lensed galaxy at z=1.7 that is magnified and stretched by the lensing cluster RCS2 032727-132623. Using this new high-resolution imaging, we modify our previous lens model (which was based on ground-based data) to fully understand the lensing geometry, and use it to reconstruct the lensed galaxy in the source plane. This giant arc represents a unique opportunity to peer into 100-pc scale structures in a high-redshift galaxy. This new source reconstruction will be crucial for a future analysis of the spatially-resolved rest-UV and rest-optical spectra of the brightest parts of the arc.

We present the discovery of an extremely bright and extended lensed source from the second Red Sequence Cluster Survey (RCS2). RCSGA 032727-132609 is spectroscopically confirmed as a giant arc and counter-image of a background galaxy at $z=1.701$, strongly-lensed by the foreground galaxy cluster RCS2 032727-132623 at $z=0.564$. The giant arc extends over $\sim 38$\,\arcsec and has an integrated $g$-band magnitude of 19.15, making it $\sim 20$ times larger and $\sim 4$ times brighter than the prototypical lensed galaxy MS1512-cB58. This is the brightest distant lensed galaxy in the Universe known to date. Its location in the `redshift desert' provides unique opportunities to connect between the large samples of galaxies known at $z\sim3$ and $z\sim1$. We have collected photometry in 9 bands, ranging from $u$ to $K_s$, which densely sample the rest-frame UV and optical light, including the age-sensitive 4000\AA\ break. A lens model is constructed for the system, and results in a robust total magnification of $2.04 \pm 0.16$ for the counter-image; we estimate an average magnification of $17.2 \pm 1.4$ for the giant arc based on the relative physical scales of the arc and counter-image. Fits of single-component spectral energy distribution (SED) models to the photometry result in a moderately young age, $t = 115 \pm 65$\,Myr, small amounts of dust, $E(B-V) \le 0.035$, and an exponentially declining star formation history with \textit{e}-folding time $\tau = 10-100$\,Myr. After correcting for the lensing magnification, we find a stellar mass of $\log(\mathrm{M}/\mathrm{M}_\odot)=10.0 \pm 0.1$. Allowing for episodic star formation, an underlying old burst could contain up to twice the mass inferred from single-component modeling. This stellar mass estimate is consistent with the average stellar mass of a sample of `BM' galaxies ($1.4 < z < 2.0$) studied by Reddy et al. (2006).

We study the mass distribution of a sample of 28 galaxy clusters using strong and weak lensing observations. The clusters are selected via their strong lensing properties as part of the Sloan Giant Arcs Survey (SGAS) from the Sloan Digital Sky Survey (SDSS). Mass modelling of the strong lensing information from the giant arcs is combined with weak lensing measurements from deep Subaru/Suprime-cam images to primarily obtain robust constraints on the concentration parameter and the shape of the mass distribution. We find that the concentration c_vir is a steep function of the mass, c_vir \propto M_vir^-0.59\pm0.12, with the value roughly consistent with the lensing-bias-corrected theoretical expectation for high mass (10^15 h^-1 M_sun) clusters. However, the observationally inferred concentration parameters appear to be much higher at lower masses (10^14 h^-1 M_sun), possibly a consequence of the modification to the inner density profiles provided by baryon cooling. The steep mass-concentration relation is also supported from direct stacking analysis of the tangential shear profiles. In addition, we explore the two-dimensional shape of the projected mass distribution by stacking weak lensing shear maps of individual clusters with prior information on the position angle from strong lens modelling, and find significant evidence for a large mean ellipticity with the best-fit value of e = 0.47 \pm 0.06 for the mass distribution of the stacked sample. We find that the luminous cluster member galaxy distribution traces the overall mass distribution very well, although the distribution of fainter cluster galaxies appears to be more extended than the total mass.

We present results from a spectroscopic program targeting 26 strong lensing cluster cores that were visually identified in the Sloan Digital Sky Survey (SDSS) and the Second Red-Sequence Cluster Survey (RCS-2). The 26 galaxy cluster lenses span a redshift range of 0.2 < z < 0.65, and our spectroscopy reveals 69 unique background sources with redshifts as high as z=5.200. We also identify redshifts for 262 cluster member galaxies and measure the velocity dispersions and dynamical masses for 18 clusters where we have redshifts for N \geq 10 cluster member galaxies. We include an accounting for the expected biases in dynamical masses of strong lensing selected clusters as predicted by results from numerical simulations and discuss possible sources of bias in our observations. The median dynamical mass of the 18 clusters with N \geq 10 spectroscopic cluster members is M_{Vir} = 7.84 x 10^14 M_sun h_{0.7}^{-1}, which is somewhat higher than predictions for strong lensing selected clusters in simulations. The disagreement is not significant considering the large uncertainty in our dynamical data, systematic uncertainties in the velocity dispersion calibration, and limitations of the theoretical modeling. Nevertheless our study represents an important first step toward characterizing large samples of clusters that are identified in a systematic way as systems exhibiting dramatic strong lensing features.

We report rest-frame optical Keck/NIRSPEC spectroscopy of the bright lensed galaxy RCSGA 032727-132609 at z=1.7037. From precise measurements of the nebular lines, we infer a number of physical properties: redshift, extinction, star formation rate, ionization parameter, electron density, electron temperature, oxygen abundance, and N/O, Ne/O, and Ar/O abundance ratios. The limit on [O III]~4363 A tightly constrains the oxygen abundance via the "direct" or electron temperature method, for the first time in an average-metallicity galaxy at z~2. We compare this result to several standard "bright-line" O abundance diagnostics, thereby testing these empirically-calibrated diagnostics in situ. Finally, we explore the positions of lensed and unlensed galaxies in standard diagnostic diagrams, and explore the diversity of ionization conditions and mass--metallicity ratios at z=2.

We measure the redshift distribution of a sample of 28 giant arcs discovered as a part of the Sloan Giant Arcs Survey (SGAS). Gemini/GMOS-North spectroscopy provides precise redshifts for 24 arcs, and "redshift desert" constraints for the remaining four. This is a direct measurement of the redshift distribution of a uniformly selected sample of bright giant arcs, which is an observable that can be used to inform efforts to predict giant arc statistics. Our primary giant arc sample has a median redshift z=1.821 and nearly two thirds of the arcs - 64% - are sources at z \gtrsim 1.4, indicating that the population of background sources that are strongly lensed into bright giant arcs resides primarily at high redshift. We also analyze the distribution of redshifts for 19 secondary strongly lensed background sources that are not visually apparent in SDSS imaging, but were identified in deeper follow-up imaging of the lensing cluster fields. Our redshift sample for the secondary sources is not spectroscopically complete, but combining it with our primary giant arc sample suggests that a large fraction of all background galaxies which are strongly lensed by foreground clusters reside at z \gtrsim 1.4. Kolmogorov-Smirnov (KS) tests indicate that our well-selected, spectroscopically complete primary giant arc redshift sample can be reproduced with a model distribution that is constructed from a combination of results from studies of strong lensing clusters in numerical simulations, and observational constraints on the galaxy luminosity function.

We have measured the Sunyaev Zel'dovich (SZ) effect for a sample of ten strong lensing selected galaxy clusters using the Sunyaev Zel'dovich Array (SZA). The SZA is sensitive to structures on spatial scales of a few arcminutes, while the strong lensing mass modeling constrains the mass at small scales (typically < 30"). Combining the two provides information about the projected concentrations of the strong lensing clusters. The Einstein radii we measure are twice as large as expected given the masses inferred from SZ scaling relations. A Monte Carlo simulation indicates that a sample randomly drawn from the expected distribution would have a larger median Einstein radius than the observed clusters about 3% of the time. The implied overconcentration has been noted in previous studies with smaller samples of lensing clusters. It persists for this sample, with the caveat that this could result from a systematic effect such as if the gas fractions of the strong lensing clusters are substantially below what is expected.

Knowledge of the supernova (SN) delay time distribution (DTD) - the SN rate versus time that would follow a hypothetical brief burst of star formation - can shed light on SN progenitors and physics. We compile recent measurements of the Type-Ia SN (SN Ia) rate in galaxy clusters at redshifts z=0-1.45. Together with the observed iron-to-stellar mass ratio in clusters, which constrains the time-integrated number of SN Ia events in clusters, we recover the DTD of SNe Ia in cluster environments. The DTD peaks at the shortest time-delay interval we probe, 0<t<2.2 Gyr, with a low tail out to delays of ~10 Gyr, and is remarkably consistent with several recent DTD reconstructions based on different methods, in different environments. We test DTD models from the literature, requiring that they simultaneously reproduce the observed cluster SN rates and the observed iron-to-stellar mass ratios. A power-law DTD of the form t^{-1.2+/-0.3}, extending to a Hubble time, can satisfy both constraints. Shallower power laws, such as t^{-1/2} cannot, assuming a single DTD, and a single star-formation burst (either brief or extended) at high z. This implies 50-85% of SNe Ia explode within 1 Gyr of star formation. DTDs from double-degenerate (DD) models, which generically have ~t^{-1} shapes over a wide range of timescales, match the data, but only if their predictions are scaled up by factors of 5-10. Single degenerate (SD) DTDs always give poor fits to the data, due to a lack of delayed SNe and overall low numbers of SNe. The observations also permit a combination of two SN Ia populations - prompt (e.g. SD) SNe Ia that explode within a few Gyr of star formation, and produce about 60% of the iron mass in clusters, and a DD population that contributes the events seen at z<1.4. Our results support the existence of a DD progenitor channel for SNe Ia, if the overall predicted numbers can be suitably increased.

We present observations of two strongly lensed $z\sim5$ Lyman-$\alpha$ Emitting (LAE) galaxies that were discovered in the Sloan Giant Arcs Survey (SGAS). We identify the two sources as SGAS J091541+382655, at $z=5.200$, and SGAS J134331+415455 at $z=4.994$. We measure their AB magnitudes at $(i,z)=(23.34\pm0.09,23.29\pm0.13$) mags and $(i,z)=(23.78\pm0.18,24.24^{+0.18}_{-0.16}$) mags, and the rest-frame equivalent widths of the Lyman-$\alpha$ emission at $25.3\pm4.1$\AA~and $135.6\pm20.3$\AA~for SGAS J091541+382655 and SGAS J134331+415455, respectively. Each source is strongly lensed by a massive galaxy cluster in the foreground, and the magnifications due to gravitational lensing are recovered from strong lens modeling of the foreground lensing potentials. We use the magnification to calculate the intrinsic, unlensed Lyman-$\alpha$ and UV continuum luminosities for both sources, as well as the implied star formation rates (SFR). We find SGAS J091541+382655 and SGAS J134341+415455 to be galaxies with (L$_{Ly-\alpha}$, L$_{UV})\leq(0.6$L$_{Ly-\alpha}^{*}, 2$L$_{UV}^{*}$) and (L$_{Ly-\alpha}$, L$_{UV})=(0.5$L$_{Ly-\alpha}^{*}, 0.9$L$_{UV}^{*}$), respectively. Comparison of the spectral energy distributions (SEDs) of both sources against stellar population models produces estimates of the mass in young stars in each galaxy: we report an upper limit of M$_{stars} \leq 7.9^{+3.7}_{-2.5} \times 10^{7}$ M$_{\sun} h_{0.7}^{-1}$ for SGAS J091531+382655, and a range of viable masses for SGAS J134331+415455 of $2\times10^{8}$ M$_{\sun} h_{0.7}^{-1} <$ M$_{stars} < 6\times10^{9}$ M$_{\sun} h_{0.7}^{-1}$.

Supernova (SN) rates are potentially powerful diagnostics of metal enrichment and SN physics, particularly in galaxy clusters with their deep, metal-retaining potentials and relatively simple star-formation histories. We have carried out a survey for supernovae (SNe) in galaxy clusters, at a redshift range 0.5<z<0.9, using the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope. We reimaged a sample of 15 clusters that were previously imaged by ACS, thus obtaining two to three epochs per cluster, in which we discovered five likely cluster SNe, six possible cluster SNe Ia, two hostless SN candidates, and several background and foreground events. Keck spectra of the host galaxies were obtained to establish cluster membership. We conducted detailed efficiency simulations, and measured the stellar luminosities of the clusters using Subaru images. We derive a cluster SN rate of 0.35 SNuB +0.17/-0.12 (statistical) \pm0.13 (classification) \pm0.01 (systematic) [where SNuB = SNe (100 yr 10^10 L_B_sun)^-1] and 0.112 SNuM +0.055/-0.039 (statistical) \pm0.042 (classification) \pm0.005 (systematic) [where SNuM = SNe (100 yr 10^10 M_sun)^-1]. As in previous measurements of cluster SN rates, the uncertainties are dominated by small-number statistics. The SN rate in this redshift bin is consistent with the SN rate in clusters at lower redshifts (to within the uncertainties), and shows that there is, at most, only a slight increase of cluster SN rate with increasing redshift. The low and fairly constant SN Ia rate out to z~1 implies that the bulk of the iron mass in clusters was already in place by z~1. The recently observed doubling of iron abundances in the intracluster medium between z=1 and 0, if real, is likely the result of redistribution of existing iron, rather than new production of iron.

We report the discovery of two strongly-lensed z ~ 3 Lyman Break Galaxies (LBGs) discovered as u-band dropouts as part of the SDSS Giant Arcs Survey (SGAS). The first, SGAS J122651.3+215220 at z = 2.9233 is lensed by one of several sub-clusters, SDSS J1226+2152, in a complex massive cluster at z = 0.43. Its (g, r, i) magnitudes are (21.14, 20.60, 20.51) which translate to surface brightnesses, mu_{g,r,i}, of (23.78, 23.11, 22.81). The second, SGAS J152745.1+065219, is an LBG at z = 2.7593 lensed by the foreground SDSS J1527+0652 at z = 0.39, with (g, r, z)=(20.90, 20.52, 20.58) and mu_{g,r,z}=(25.15, 24.52, 24.12). Moderate resolution spectroscopy confirms the redshifts suggested by photometric breaks, and shows both absorption and emission features typical of LBGs. Lens mass models derived from combined imaging and spectroscopy reveal that SGAS J122651.3+215220 is a highly magnified source (M ~40), while SGAS J152745.1+065219 is magnified by no more than M ~ 15. Compared to LBG survey results (Steidel et al. 2003), the luminosities and lensing-corrected magnitudes suggest that SGAS J122651.3+215220 is among the faintest 20% of LBGs in that sample. SGAS J152745.1+065219, on the other hand, appears to be more representative of the average LBG, similar to the "Cosmic Eye".

We derive radial mass profiles of four strong lensing selected clusters which show prominent giant arcs (Abell 1703, SDSS J1446+3032, SDSS J1531+3414, and SDSS J2111-0115), by combining detailed strong lens modeling with weak lensing shear measured from deep Subaru Suprime-cam images. Weak lensing signals are detected at high significance for all four clusters, whose redshifts range from z=0.28 to 0.64. We demonstrate that adding strong lensing information with known arc redshifts significantly improves constraints on the mass density profile, compared to those obtained from weak lensing alone. While the mass profiles are well fitted by the universal form predicted in N-body simulations of the LCDM model, all four clusters appear to be slightly more centrally concentrated (the concentration parameters c_vir \sim 8) than theoretical predictions, even after accounting for the bias toward higher concentrations inherent in lensing selected samples. Our results are consistent with previous studies which similarly detected a concentration excess, and increases the total number of clusters studied with the combined strong and weak lensing technique to ten. Combining our sample with previous work, we find that clusters with larger Einstein radii are more anomalously concentrated. We also present a detailed model of the lensing cluster A1703 with constraints from multiple image families, and find the dark matter inner density profile to be cuspy with the slope consistent with -1, in agreement with expectations.

We present the results of deep spectroscopy for the central region of the cluster lens SDSS J1004+4112 with the Subaru telescope. A secure detection of an emission line of the faint blue stellar object (component E) near the center of the brightest cluster galaxy (G1) confirms that it is the central fifth image of the lensed quasar system. In addition, we measure the stellar velocity dispersion of G1 to be sigma_* = 352+-13 km/s. We combine these results to obtain constraints on the mass M_BH of the putative black hole (BH) at the center of the inactive galaxy G1, and hence on the M_BH-sigma_* relation at the lens redshift z_l=0.68. From detailed mass modeling, we place an upper limit on the black hole mass, M_BH < 2.1x10^{10}M_sun at 1-sigma level (<3.1x10^{10}M_sun at 3-sigma), which is consistent with black hole masses expected from the local and redshift-evolved M_BH-sigma_* relations, M_BH~10^{9}-10^{10}M_sun.

We report a measurement of the supernova (SN) rates (Ia and core-collapse) in galaxy clusters based on the 136 SNe of the sample described in Cappellaro et al. (1999) and Mannucci et al. (2005). Early-type cluster galaxies show a type Ia SN rate (0.066 SNuM) similar to that obtained by Sharon et al. (2007) and more than 3 times larger than that in field early-type galaxies (0.019 SNuM). This difference has a 98% statistical confidence level. We examine many possible observational biases which could affect the rate determination, and conclude that none of them is likely to significantly alter the results. We investigate how the rate is related to several properties of the parent galaxies, and find that cluster membership, morphology and radio power all affect the SN rate, while galaxy mass has no measurable effect. The increased rate may be due to galaxy interactions in clusters, inducing either the formation of young stars or a different evolution of the progenitor binary systems. We present the first measurement of the core-collapse SN rate in cluster late-type galaxies, which turns out to be comparable to the rate in field galaxies. This suggests that no large systematic difference in the initial mass function exists between the two environments.

Large samples of high-redshift supernovae (SNe) are potentially powerful probes of cosmic star formation, metal enrichment, and SN physics. We present initial results from a new deep SN survey, based on re-imaging in the R, i', z' bands, of the 0.25 deg2 Subaru Deep Field (SDF), with the 8.2-m Subaru telescope and Suprime-Cam. In a single new epoch consisting of two nights of observations, we have discovered 33 candidate SNe, down to a z'-band magnitude of 26.3 (AB). We have measured the photometric redshifts of the SN host galaxies, obtained Keck spectroscopic redshifts for 17 of the host galaxies, and classified the SNe using the Bayesian photometric algorithm of Poznanski et al. (2007) that relies on template matching. After correcting for biases in the classification, 55% of our sample consists of Type Ia supernovae and 45% of core-collapse SNe. The redshift distribution of the SNe Ia reaches z ~ 1.6, with a median of z ~ 1.2. The core-collapse SNe reach z ~ 1.0, with a median of z ~ 0.5. Our SN sample is comparable to the Hubble Space Telescope/GOODS sample both in size and redshift range. The redshift distributions of the SNe in the SDF and in GOODS are consistent, but there is a trend (which requires confirmation using a larger sample) for more high-z SNe Ia in the SDF. This trend is also apparent when comparing the SN Ia rates we derive to those based on GOODS data. Our results suggest a fairly constant rate at high redshift that could be tracking the star-formation rate. Additional epochs on this field, already being obtained, will enlarge our SN sample to the hundreds, and determine whether or not there is a decline in the SN Ia rate at z >~ 1.

We present results of the HST Advanced Camera for Surveys spectroscopic ground-based redshift survey in the field of A1689. We measure 98 redshifts, increasing the number of spectroscopically confirmed objects by sixfold. We present two spectra from this catalog of the Sextet Arcs, images which arise from a strongly-lensed Lyman Break Galaxy (LBG) at a redshift of z=3.038. Gravitational lensing by the cluster magnifies its flux by a factor of ~16 and produces six separate images with a total r-band magnitude of r_625=21.7. The two spectra, each of which represents emission from different regions of the LBG, show H I and interstellar metal absorption lines at the systemic redshift. Significant variations are seen in Ly-alpha profile across a single galaxy, ranging from strong absorption to a combination of emission plus absorption. A spectrum of a third image close to the brightest arc shows Ly-alpha emission at the same redshift as the LBG, arising from either another spatially distinct region of the galaxy, or from a companion galaxy close to the LBG. Taken as a group, the Ly-alpha equivalent width in these three spectra decreases with increasing equivalent width of the strongest interstellar absorption lines. We discuss how these variations can be used to understand the physical conditions in the LBG. Intrinsically, this LBG is faint, ~0.1L*, and forming stars at a modest rate, ~4 solar masses per year. We also detect absorption line systems toward the Sextet Arcs at z=2.873 and z=2.534. The latter system is seen across two of our spectra.

On 2006 May 5, a four second duration, low energy, ~10^49 erg, Gamma-Ray Burst (GRB) was observed, spatially associated with a z=0.0894 galaxy. Here, we report the discovery of the GRB optical afterglow and observations of its environment using Gemini-south, Hubble Space Telescope (HST), Chandra, Swift and the Very Large Array. The optical afterglow of this GRB is spatially associated with a prominent star forming region in the Sc-type galaxy 2dFGRS S173Z112. Its proximity to a star forming region suggests that the progenitor delay time, from birth to explosion, is smaller than about 10 Myr. Our HST deep imaging rules out the presence of a supernova brighter than an absolute magnitude of about -11 (or -12.6 in case of ``maximal'' extinction) at about two weeks after the burst, and limits the ejected mass of radioactive Nickel 56 to be less than about 2x10^-4 solar mass (assuming no extinction). Although it was suggested that GRB 060505 may belong to a new class of long-duration GRBs with no supernova, we argue that the simplest interpretation is that the physical mechanism responsible for this burst is the same as for short-duration GRBs.

Supernova (SN) rates are a potentially powerful diagnostic of star formation history (SFH), metal enrichment, and SN physics, particularly in galaxy clusters with their deep, metal-retaining potentials, and simple SFH. However, a low-redshift cluster SN rate has never been published. We derive the SN rate in galaxy clusters at 0.06<z<0.19, based on type Ia supernovae (SNe Ia) that were discovered by the Wise Observatory Optical Transient Survey. As described in a separate paper, a sample of 140 rich Abell clusters was monitored, in which six cluster SNe Ia were found and confirmed spectroscopically. Here, we determine the SN detection efficiencies of the individual survey images, and combine the efficiencies with the known spectral properties of SNe Ia to calculate the effective visibility time of the survey. The cluster stellar luminosities are measured from the Sloan Digital Sky Survey (SDSS) database in the griz SDSS bands. Uncertainties are estimated using Monte-Carlo simulations in which all input parameters are allowed to vary over their known distributions. We derive SN rates normalized by stellar luminosity, in SNU units (SNe per century per 10^10 L_sun) in five photometric bandpasses, of 0.36+/-(0.22,0.14)+/-0.02 (B), 0.351+/-(0.210,0.139)+/-0.020 (g), 0.288+/-(0.172,0.114)+/-0.018 (r), 0.229+/-(0.137,0.091)+/-0.014 (i), 0.186+/-(0.111,0.074)+/-0.010 (z), where the quoted errors are statistical and systematic, respectively. The SN rate per stellar mass unit, derived using a color-luminosity-mass relation, is 0.098+/-(0.059,0.039)+/-0.009 SNe (century 10^10 M_sun)^-1. The low cluster SN rates we find are similar to, and consistent with, the SN Ia rate in local elliptical galaxies.

We describe our ongoing program designed to measure the SN-Ia rate in a sample of massive z=0.5-0.9 galaxy clusters. The SN-Ia rate is a poorly known observable, especially at high z, and in cluster environments. The SN rate and its redshift dependence can serve as powerful discrimiminants for a number of key issues in astrophysics and cosmology. Our observations will put clear constraints on the characteristic SN-Ia ``delay time'', the typical time between the formation of a stellar population and the explosion of some of its members as SNe-Ia. Such constraints can exclude entire categories of SN-Ia progenitor models, since different models predict different delays. These data will also help to resolve the question of the dominant source of the high metallicity in the intracluster medium (ICM) - SNe-Ia, or core-collapse SNe from an early stellar population with a top-heavy IMF, perhaps those population III stars responsible for the early re-ionization of the Universe. Since clusters are excellent laboratories for studying enrichment (they generally have a simple star-formation history, and matter cannot leave their deep potentials), the results will be relevant for understanding metal enrichment in general, and the possible role of first generation stars in early Universal enrichment. Observations obtained so far during cycles 14 and 15 yield many SNe in our cluster fields, but our follow-up campaign reveals most are not in cluster galaxies.

Over the past decade our physical understanding of gamma-ray bursts (GRBs) has progressed rapidly thanks to the discovery and observation of their long-lived afterglow emission. Long-duration (T < 2 s) GRBs are associated with the explosive deaths of massive stars (``collapsars''), which produce accompanying supernovae, while the short-duration (T > 2 s) GRBs arise from a different origin, which has been argued to be the merger of two compact objects, either neutron stars or black holes. Here we present observations of GRB060614, a 100-s long burst discovered by the Swift satellite, which require the invocation of a new explosive process: either a massive ``collapsar'' that powers a GRB without any associated supernova, or a new type of engine, as long-lived as the collapsar but without any such massive stellar host. We also discuss the properties of this burst's redshift z=0.125 host galaxy, which distinguish it from other long-duration GRBs and suggest that an entirely new type of GRB progenitor may be required.

We have identified three multiply imaged galaxies in Hubble Space Telescope images of the redshift z=0.68 cluster responsible for the large-separation quadruply lensed quasar, SDSS J1004+4112. Spectroscopic redshifts have been secured for two of these systems using the Keck I 10m telescope. The most distant lensed galaxy, at z=3.332, forms at least four images, and an Einstein ring encompassing 3.1 times more area than the Einstein ring of the lensed QSO images at z=1.74, due to the greater source distance. For a second multiply imaged galaxy, we identify Ly_alpha emission at a redshift of z=2.74. The cluster mass profile can be constrained from near the center of the brightest cluster galaxy, where we observe both a radial arc and the fifth image of the lensed quasar, to the Einstein radius of the highest redshift galaxy, ~110 kpc. Our preliminary modeling indicates that the mass approximates an elliptical body, with an average projected logarithmic gradient of ~-0.5. The system is potentially useful for a direct measurement of world models in a previously untested redshift range.

ACS observations of massive lensing clusters permit an order of magnitude increase in the numbers of multiply-lensed background galaxies identified behind a given cluster. We have developed a code to take the pixels belonging to any given image and generate counter-images with full resolution, so that multiple systems are convincingly and exhaustively identified. Over 130 images of 35 multiply lensed galaxies are found behind A1689, including many radial arcs and also tiny counter-images projected on the center of mass. The derived mass profile is found to flatten steadily towards the center, like an NFW profile, with a mean slope d(log(Sigma))/d(log(r)) ~ -0.55 +/- 0.1, over the range r<250 kpc/h, which is somewhat steeper than predicted for such a massive halo. We also clearly see the expected geometric increase of bend angles with redshift, however, given the low redshift of A1689, z=0.18, the dependence on cosmological parameters is weak, but using higher redshift clusters from our GTO program we may derive a more competitive constraint.

We analyse deep multi-colour Advanced Camera images of the largest known gravitational lens, A1689. Radial and tangential arcs delineate the critical curves in unprecedented detail and many small counter-images are found near the center of mass. We construct a flexible light deflection field to predict the appearance and positions of counter-images. The model is refined as new counter-images are identified and incorporated to improve the model, yielding a total of 106 images of 30 multiply lensed background galaxies, spanning a wide redshift range, 1.0$<$z$<$5.5. The resulting mass map is more circular in projection than the clumpy distribution of cluster galaxies and the light is more concentrated than the mass within $r<50kpc/h$. The projected mass profile flattens steadily towards the center with a shallow mean slope of $d\log\Sigma/d\log r \simeq -0.55\pm0.1$, over the observed range, r$<250kpc/h$, matching well an NFW profile, but with a relatively high concentration, $C_{vir}=8.2^{+2.1}_{-1.8}$. A softened isothermal profile ($r_{core}=20\pm2$\arcs) is not conclusively excluded, illustrating that lensing constrains only projected quantities. Regarding cosmology, we clearly detect the purely geometric increase of bend-angles with redshift. The dependence on the cosmological parameters is weak due to the proximity of A1689, $z=0.18$, constraining the locus, $\Omega_M+\Omega_{\Lambda} \leq 1.2$. This consistency with standard cosmology provides independent support for our model, because the redshift information is not required to derive an accurate mass map. Similarly, the relative fluxes of the multiple images are reproduced well by our best fitting lens model.

We present the results of several surveys for supernovae (SNe) in galaxy clusters. SNe discovered in deep, archival HST images were used to measure the cluster SN Ia rate to z=1. A search for SNe in nearby (0.06 < z < 0.2) Abell galaxy clusters yielded 15 SNe, 12 of which were spectroscopically confirmed. Of these, 7 are cluster SNe Ia, which we will use to measure the SN Ia rate in nearby clusters. This search has also discovered the first convincing examples of intergalactic SNe. We conclude with a brief description of ongoing and future cluster SN surveys.

We have searched for high-redshift supernova (SN) candidates in multiple deep Hubble Space Telescope (HST) archival images of nine galaxy-cluster fields. We detect six apparent SNe, with I814 between 21.6 and 28.4 mag. There is roughly 1 SN per deep (flux limit I814 > 26 mag), doubly-imaged, WFPC2 cluster field. Two SNe are associated with cluster galaxies (at redshifts z=0.18 and z=0.83), three are probably in galaxies not in the clusters (at z=0.49, z=0.60, and z=0.98), and one is at unknown z. After accounting for observational efficiencies and uncertainties (statistical and systematic) we derive the rate of type-Ia SNe within the projected central 500 kpc of rich clusters: R=0.20(+0.84)(-0.19) SNu in clusters at z=0.18 to 0.37, and R=0.41(+1.23)(-0.39) SNu in clusters at z=0.83 to 1.27 (95 per cent confidence interval; H_0=50; 1 SNu = 1 SN per century per 10^10 L_B_sun). Combining the two redshift bins, the SN rate at a mean redshift of z=0.41 is R(z=0.41) = 0.30(+0.58)(-0.28) SNu. The upper bounds argue against SNe-Ia being the dominant source of the large iron mass observed in the intra-cluster medium. We also compare our observed counts of field SNe (i.e., non-cluster SNe of all types) to recent model predictions. The observed field count is zero or one SN with I814 < 26 mag, and 1 to 3 SNe with I814 < 27 mag. These counts are about two times lower than some of the predictions. Since the counts at these magnitudes are likely dominated by type-II SNe, our observations may suggest obscuration of distant SNe-II, or a SN-II luminosity distribution devoid of a large high-luminosity tail.