We study the projected spatial offset between the ultraviolet continuum and Ly$\alpha$ emission for a sample of 65 lensed and unlensed galaxies in the Epoch of Reionization (5<z<7), the first such study at these redshifts, in order to understand the potential for these offsets to confuse estimates of the Ly$\alpha$ properties of galaxies observed in slit spectroscopy. While we find that ~40% of galaxies in our sample show significant projected spatial offsets ($|\Delta_{Ly\alpha-UV}|$), we find a modest average offset of 0.61$\pm$0.08 kpc. A small fraction of our sample, ~10%, exhibits offsets of 2-4 kpc, sizes that are larger than the effective radii of typical galaxies at these redshifts. An internal comparison and a comparison to studies at lower redshift yielded no significant evidence of evolution of $|\Delta_{Ly\alpha-UV}|$ with redshift. In our own sample, UV-bright galaxies showed offsets a factor of three greater than their fainter counterparts, 0.89$\pm$0.18 vs. 0.27$\pm$0.05 kpc, respectively. We argue that offsets are likely not the result of merging processes, but are rather due to anisotropic internal processes internal to the galaxies as a result of stellar feedback, which facilitates Ly$\alpha$ fluorescence and/or backscattering from nearby or outflowing gas. The reduction in the Ly$\alpha$ flux due to offset effects for various observational setups was quantified through mock observations of simple simulations. It was found that the loss of Ly$\alpha$ photons for galaxies with average offsets is not, if corrected for, a limiting factor for all but the narrowest slit widths (<0.4''). However, for the largest offsets, if such offsets are mostly perpendicular to the slit major axis, slit losses were found to be extremely severe in cases where slit widths of $\le$1'' were employed, such as those planned for James Webb Space Telescope/NIRSpec observations. (abridged)

When did galaxies start forming stars? What is the role of distant galaxies in galaxy formation models and the epoch of reionization? What are the conditions in typical star-forming galaxies at redshifts >~4? Why is galaxy evolution dependent on environment? The Spitzer Space Telescope has been a crucial tool for addressing these questions. Accurate knowledge of stellar masses, ages and star formation rates requires measuring rest-frame optical (and ultraviolet) light, which only Spitzer can probe at high redshifts for a sufficiently large sample of typical galaxies. Many of these science goals are the main science drivers for the James Webb Space Telescope, and Spitzer afforded us their first exploration.

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

We present 36 spectroscopically confirmed intrinsically UV-faint Ly$\alpha$ emitting galaxies from follow-up observations with Keck/DEIMOS of gravitationally lensed high-redshift candidates. Candidates were selected to be between $5\lesssim z \lesssim 7$ from photometric data using \textit{HST} and \textit{Spitzer} imaging surveys. We used photometric redshift information to perform an integrated photometric redshift probability cut $>1\%$ between $5<z<7$ to construct a sample of 198 high-redshift objects. Our high-redshift sample spans intrinsic UV luminosities from a few $L^*$ down to $0.001L^*$. We identified 19 high-confidence detections of Ly$\alpha$ and an additional 17 likely detections. We divided our sample into lower-redshift ($z\sim5.5$) and higher-redshift ($z\sim6.5$) bins and ran Monte Carlo trials, incorporating the strength of the Ly$\alpha$ emission and the photometric redshift of the non-detections. Considering only objects where Ly$\alpha$ could be detected at EW(Ly$\alpha$)$>$25{\AA} at $3\sigma$ at the fiducial depth of our survey, and only those galaxies with EW(Ly$\alpha$)$>$25{\AA} as true LAEs, and finally, only objects with $m_{AB}<26.8$, we found the LAE fraction to be flat, or modestly increase from 0.26$\pm0.04$ to 0.30$\pm0.04$. These values relative to those for lower-redshift samples are consistent with a rising LAE fraction with redshift out to $z\sim6$, but at $z\sim6.5$ there is some tension between our results and results from surveys at intrinsically brighter luminosities. We conclude intrinsically fainter galaxies have Ly$\alpha$ emission, and there is a steep drop in the LAE fraction from our high-redshift sample at $z\sim6.5$ and from similar galaxies at $z\sim7.5$. This likely indicates we are witnessing the tail end of the epoch of reionization, as such a drop is not expected due to changes of intrinsic galaxy properties between these redshifts.

We constrain the distribution of spatially offset Lyman-alpha emission (Ly$\alpha$) relative to rest-frame ultraviolet emission in $\sim300$ high redshift ($3<z<5.5$) Lyman-break galaxies (LBGs) exhibiting Ly$\alpha$ emission from VANDELS, a VLT/VIMOS slit-spectroscopic survey of the CANDELS Ultra Deep Survey and Chandra Deep Field South fields (${\simeq0.2}~\mathrm{deg}^2$ total). Because slit spectroscopy compresses two-dimensional spatial information into one spatial dimension, we use Bayesian inference to recover the underlying Ly$\alpha$ spatial offset distribution. We model the distribution using a 2D circular Gaussian, defined by a single parameter $\sigma_{r,\mathrm{Ly}\alpha}$, the standard deviation expressed in polar coordinates. Over the entire redshift range of our sample ($3<z<5.5$), we find $\sigma_{r,\mathrm{Ly}\alpha}=1.70^{+0.09}_{-0.08}$ kpc ($68\%$ conf.), corresponding to $\sim0.25$ arcsec at $\langle z\rangle=4.5$. We also find that $\sigma_{r,\mathrm{Ly}\alpha}$ decreases significantly with redshift. Because Ly$\alpha$ spatial offsets can cause slit-losses, the decrease in $\sigma_{r,\mathrm{Ly}\alpha}$ with redshift can partially explain the increase in the fraction of Ly$\alpha$ emitters observed in the literature over this same interval, although uncertainties are still too large to reach a strong conclusion. If $\sigma_{r,\mathrm{Ly}\alpha}$ continues to decrease into the reionization epoch, then the decrease in Ly$\alpha$ transmission from galaxies observed during this epoch might require an even higher neutral hydrogen fraction than what is currently inferred. Conversely, if spatial offsets increase with the increasing opacity of the IGM, slit losses may explain some of the drop in Ly$\alpha$ transmission observed at $z>6$. Spatially resolved observations of Ly$\alpha$ and UV continuum at $6<z<8$ are needed to settle the issue.

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

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 a large spectroscopic campaign with Keck/MOSFIRE targeting Lyman-alpha emission (Ly$\alpha$) from intrinsically faint Lyman-break Galaxies (LBGs) behind 12 efficient galaxy cluster lenses. Gravitational lensing allows us to probe the more abundant faint galaxy population to sensitive Ly$\alpha$ equivalent width limits. During the campaign we targeted 70 LBG candidates with MOSFIRE Y-band, selected photometrically to cover Ly$\alpha$ over the range $7<z<8.2$. We detect $S/N>5$ emission lines in 2 of these galaxies and find that they are likely Ly$\alpha$ at $z=7.148\pm0.001$ and $z=7.161\pm0.001$. We present new lens models for 4 of the galaxy clusters, using our previously published lens models for the remaining clusters to determine the magnification factors for the source galaxies. Using a Bayesian framework that employs large scale reionization simulations of the intergalactic medium (IGM) as well as realistic properties of the interstellar medium and circumgalactic medium, we infer the volume-averaged neutral hydrogen fraction, $\overline{x}_{\mathrm{HI}}$, in the IGM during reionization to be $\overline{x}_{\mathrm{HI}}=0.88^{+0.05}_{-0.10}$ at $z=7.6\pm0.6$. Our result is consistent with a late and rapid reionization scenario inferred by Planck.

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.

Golovich et al. 2017b presents an optical imaging and spectroscopic survey of 29 radio relic merging galaxy clusters. In this paper, we study this survey to identify substructure and quantify the dynamics of the mergers. Using a combined photometric and spectroscopic approach, we identify the minimum number of substructures in each system to describe the galaxy populations and estimate the line of sight velocity difference between likely merging subclusters. We find that the line-of-sight velocity components of the mergers are typically small compared with the maximum three dimensional relative velocity (usually $<1000$ km s$^{-1}$ and often consistent with zero). This suggests that the merger axes of these systems are generally in or near the plane of the sky matching findings in magneto-hydrodynamical simulations. In 28 of the 29 systems we identify substructures in the galaxy population aligned with the radio relic(s) and presumed associated merger induced shock. From this ensemble, we identify eight systems to include in a `gold' sample that is prime for further observation, modeling, and simulation study. Additional papers will present weak lensing mass maps and dynamical modeling for each merging system, ultimately leading to new insight into a wide range of astrophysical phenomena at some of the largest scales in the universe.

We present a gravitational lensing model of MACS J1149.5+2223 using ultra-deep Hubble Frontier Fields imaging data and spectroscopic redshifts from HST grism and VLT/MUSE spectroscopic data. We create total mass maps using 38 multiple images (13 sources) and 608 weak lensing galaxies, as well as 100 multiple images of 31 star-forming regions in the galaxy that hosts Supernova Refsdal. We find good agreement with a range of recent models within the HST field of view. We present a map of the ratio of projected stellar mass to total mass ($f_{\star}$), and find that the stellar mass fraction for this cluster peaks on the primary BCG. Averaging within a radius of 0.3 Mpc, we obtain a value of $\langle f_{\star} \rangle = 0.012^{+0.004}_{-0.003}$, consistent with other recent results for this ratio in cluster environments, though with a large global error (up to $\delta f_{\star} = 0.005$) primarily due to the choice of an IMF. We compare values of $f_{\star}$ and measures of star formation efficiency for this cluster to other Hubble Frontier Fields clusters studied in the literature, finding that MACS1149 has a higher stellar mass fraction than these other clusters, but a star formation efficiency typical of massive clusters.

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.

We present deep spectroscopic observations of a Lyman-break galaxy candidate (hereafter MACS1149-JD) at $z\sim9.5$ with the $\textit{Hubble}$ Space Telescope ($\textit{HST}$) WFC3/IR grisms. The grism observations were taken at 4 distinct position angles, totaling 34 orbits with the G141 grism, although only 19 of the orbits are relatively uncontaminated along the trace of MACS1149-JD. We fit a 3-parameter ($z$, F160W mag, and Ly$\alpha$ equivalent width) Lyman-break galaxy template to the three least contaminated grism position angles using an MCMC approach. The grism data alone are best fit with a redshift of $z_{\mathrm{grism}}=9.53^{+0.39}_{-0.60}$ ($68\%$ confidence), in good agreement with our photometric estimate of $z_{\mathrm{phot}}=9.51^{+0.06}_{-0.12}$ ($68\%$ confidence). Our analysis rules out Lyman-alpha emission from MACS1149-JD above a $3\sigma$ equivalent width of 21 \AA{}, consistent with a highly neutral IGM. We explore a scenario where the red $\textit{Spitzer}$/IRAC $[3.6] - [4.5]$ color of the galaxy previously pointed out in the literature is due to strong rest-frame optical emission lines from a very young stellar population rather than a 4000 \AA{} break. We find that while this can provide an explanation for the observed IRAC color, it requires a lower redshift ($z\lesssim9.1$), which is less preferred by the $\textit{HST}$ imaging data. The grism data are consistent with both scenarios, indicating that the red IRAC color can still be explained by a 4000 \AA{} break, characteristic of a relatively evolved stellar population. In this interpretation, the photometry indicate that a $340^{+29}_{-35}$ Myr stellar population is already present in this galaxy only $\sim500~\mathrm{Myr}$ after the Big Bang.

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.

Modern data empower observers to describe galaxies as the spatially and biographically complex objects they are. We illustrate this through case studies of four, $z\sim1.3$ systems based on deep, spatially resolved, 17-band + G102 + G141 Hubble Space Telescope grism spectrophotometry. Using full spectrum rest-UV/-optical continuum fitting, we characterize these galaxies' observed $\sim$kpc-scale structures and star formation rates (SFRs) and reconstruct their history over the age of the universe. The sample's diversity---passive to vigorously starforming; stellar masses $\log M_*/M_\odot=10.5$ to $11.2$---enables us to draw spatio-temporal inferences relevant to key areas of parameter space (Milky Way- to super-Andromeda-mass progenitors). Specifically, we find signs that bulge mass-fractions ($B/T$) and SF history shapes/spatial uniformity are linked, such that higher $B/T$s correlate with "inside-out growth" and central specific SFRs that peaked above the global average for all starforming galaxies at that epoch. Conversely, the system with the lowest $B/T$ had a flat, spatially uniform SFH with normal peak activity. Both findings are consistent with models positing a feedback-driven connection between bulge formation and the switch from rising to falling SFRs ("quenching"). While sample size forces this conclusion to remain tentative, this work provides a proof-of-concept for future efforts to refine or refute it: JWST, WFIRST, and the 30-m class telescopes will routinely produce data amenable to this and more sophisticated analyses. These samples---spanning representative mass, redshift, SFR, and environmental regimes---will be ripe for converting into thousands of sub-galactic-scale empirical windows on what individual systems actually looked like in the past, ushering in a new dialog between observation and theory.

Within one billion years of the Big Bang, intergalactic hydrogen was ionized by sources emitting ultraviolet and higher energy photons. This was the final phenomenon to globally affect all the baryons (visible matter) in the Universe. It is referred to as cosmic reionization and is an integral component of cosmology. It is broadly expected that intrinsically faint galaxies were the primary ionizing sources due to their abundance in this epoch. However, at the highest redshifts ($z>7.5$; lookback time 13.1 Gyr), all galaxies with spectroscopic confirmations to date are intrinsically bright and, therefore, not necessarily representative of the general population. Here, we report the unequivocal spectroscopic detection of a low luminosity galaxy at $z>7.5$. We detected the Lyman-$\alpha$ emission line at $\sim 10504$ {\AA} in two separate observations with MOSFIRE on the Keck I Telescope and independently with the Hubble Space Telescope's slit-less grism spectrograph, implying a source redshift of $z = 7.640 \pm 0.001$. The galaxy is gravitationally magnified by the massive galaxy cluster MACS J1423.8+2404 ($z = 0.545$), with an estimated intrinsic luminosity of $M_{AB} = -19.6 \pm 0.2$ mag and a stellar mass of $M_{\star} = 3.0^{+1.5}_{-0.8} \times 10^8$ solar masses. Both are an order of magnitude lower than the four other Lyman-$\alpha$ emitters currently known at $z > 7.5$, making it probably the most distant representative source of reionization found to date.

Exploiting the data of the Grism Lens-Amplified Survey from Space (GLASS), we characterize the spatial distribution of star formation in 76 high star forming galaxies in 10 clusters at 0.3< z <0.7. All these galaxies are likely restricted to first infall. In a companion paper we contrast the properties of field and cluster galaxies, whereas here we correlate the properties of H{\alpha} emitters to a number of tracers of the cluster environment to investigate its role in driving galaxy transformations. H{\alpha} emitters are found in the clusters out to 0.5 virial radii, the maximum radius covered by GLASS. The peak of the H{\alpha} emission is offset with respect to the peak of the UV-continuum. We decompose this offsets into a radial and tangential component. The radial compo- nent points away from the cluster center in 60% of the cases, with 95% confidence. The decompositions agree with cosmological simulations, i.e. the H{\alpha} emission offset correlates with galaxy velocity and ram-pressure stripping signatures. Trends between H{\alpha} emitter properties and surface mass density distributions and X-ray emissions emerge only for unrelaxed clusters. The lack of strong correlations with the global environment does not allow us to identify a unique environmental effect originating from the cluster center. In contrast, correla- tions between H{\alpha} morphology and local number density emerge. We conclude that local effects, uncorrelated to the cluster-centric radius, play a more important role in shaping galaxy properties.

We present the results of ALMA spectroscopic follow-up of a $z=6.765$ Lyman-$\alpha$ emitting galaxy behind the cluster RXJ1347-1145. We report the detection of [CII]158$\mu$m line fully consistent with the Lyman-$\alpha$ redshift and with the peak of the optical emission. Given the magnification of $\mu=5.0 \pm 0.3$ the intrinsic (corrected for lensing) luminosity of the [CII] line is $L_{[CII]} =1.4^{+0.2}_{-0.3} \times 10^7L_{\odot}$, which is ${\sim}5$ times fainter than other detections of $z\sim 7$ galaxies. The result indicates that low $L_{[CII]}$ in $z\sim 7$ galaxies compared to the local counterparts might be caused by their low metallicities and/or feedback. The small velocity off-set ($\Delta v = 20_{-40}^{+140} \rm km/s$) between the Lyman-$\alpha$ and [CII] line is unusual, and may be indicative of ionizing photons escaping.

Exploiting the slitless spectroscopy taken as part of the Grism Lens-Amplified Survey from Space (GLASS), we present an extended analysis of the spatial distribution of star formation in 76 galaxies in 10 clusters at 0.3< z <0.7. We use 85 foreground and background galaxies in the same redshift range as a field sample. The samples are well matched in stellar mass (10^8-10^11 M_sun) and star formation rate (0.5-50 M_sun/yr). We visually classify galaxies in terms of broad-band morphology, Halpha morphology and likely physical process acting on the galaxy. Most Halpha emitters have a spiral morphology (41+/-8% in clusters, 51+/-8% in the field), followed by mergers/interactions (28+/-8%, 31+/-7%, respectively) and early-type galaxies (remarkably as high as 29+/-8% in clusters and 15+/-6% in the field). A diversity of Halpha morphologies is detected, suggesting a diversity of physical processes. In clusters, 30+/-8% of the galaxies present a regular morphology, mostly consistent with star formation diffused uniformly across the stellar population (mostly in the disk component, when present). The second most common morphology (28+/-8%) is asymmetric/jellyfish, consistent with ram pressure stripping or other non-gravitational processes in 18+/-8% of the cases. Ram pressure stripping appears significantly less prominent in the field (2+/-2%), where the most common morphology/mechanism appears to be consistent with minor gas rich mergers or clump accretion. This work demonstrates that while environment specific mechanisms affect galaxy evolution at this redshift, they are diverse and their effects subtle. A full understanding of this complexity requires larger samples and detailed and spatially resolved physical models.

We report the detection of Ly$\alpha$ emission at $\sim9538$\AA{} in the Keck/DEIMOS and \HST WFC3 G102 grism data from a triply-imaged galaxy at $z=6.846\pm0.001$ behind galaxy cluster MACS J2129.4$-$0741. Combining the emission line wavelength with broadband photometry, line ratio upper limits, and lens modeling, we rule out the scenario that this emission line is \oii at $z=1.57$. After accounting for magnification, we calculate the weighted average of the intrinsic Ly$\alpha$ luminosity to be $\sim1.3\times10^{42}~\mathrm{erg}~\mathrm{s}^{-1}$ and Ly$\alpha$ equivalent width to be $74\pm15$\AA{}. Its intrinsic UV absolute magnitude at 1600\AA{} is $-18.6\pm0.2$ mag and stellar mass $(1.5\pm0.3)\times10^{7}~M_{\odot}$, making it one of the faintest (intrinsic $L_{UV}\sim0.14~L_{UV}^*$) galaxies with Ly$\alpha$ detection at $z\sim7$ to date. Its stellar mass is in the typical range for the galaxies thought to dominate the reionization photon budget at $z\gtrsim7$; the inferred Ly$\alpha$ escape fraction is high ($\gtrsim 10$\%), which could be common for sub-$L^*$ $z\gtrsim7$ galaxies with Ly$\alpha$ emission. This galaxy offers a glimpse of the galaxy population that is thought to drive reionization, and it shows that gravitational lensing is an important avenue to probe the sub-$L^*$ galaxy population.

We present a strong and weak gravitational lens model of the galaxy cluster MACSJ0416.1-2403, constrained using spectroscopy from the Grism Lens-Amplified Survey from Space (GLASS) and Hubble Frontier Fields (HFF) imaging data. We search for emission lines in known multiply imaged sources in the GLASS spectra, obtaining secure spectroscopic redshifts of 31 multiple images belonging to 16 distinct source galaxies. The GLASS spectra provide the first spectroscopic measurements for 6 of the source galaxies. The weak lensing signal is acquired from 884 galaxies in the F606W HFF image. By combining the weak lensing constraints with 15 multiple image systems with spectroscopic redshifts and 9 multiple image systems with photometric redshifts, we reconstruct the gravitational potential of the cluster on an adaptive grid. The resulting total mass density map is compared with a stellar mass density map obtained from the deep Spitzer Frontier Fields imaging data to study the relative distribution of stellar and total mass in the cluster. We find that the projected stellar mass to total mass ratio, $f_{\star}$, varies considerably with the stellar surface mass density. The mean projected stellar mass to total mass ratio is $\langle f_{\star} \rangle= 0.009 \pm 0.003 $ (stat.), but with a systematic error as large as $0.004-0.005$, dominated by the choice of the IMF. We find agreement with several recent measurements of $f_{\star}$ in massive cluster environments. The lensing maps of convergence, shear, and magnification are made available to the broader community in the standard HFF format.

We study the stellar population properties of the IRAC-detected $6 \lesssim z \lesssim 10$ galaxy candidates from the Spitzer UltRa Faint SUrvey Program (SURFS UP). Using the Lyman Break selection technique, we find a total of 16 new galaxy candidates at $6 \lesssim z \lesssim 10$ with $S/N \geq 3$ in at least one of the IRAC $3.6\mu$m and $4.5\mu$m bands. According to the best mass models available for the surveyed galaxy clusters, these IRAC-detected galaxy candidates are magnified by factors of $\sim 1.2$--$5.5$. We find that the IRAC-detected $6 \lesssim z \lesssim 10$ sample is likely not a homogeneous galaxy population: some are relatively massive (stellar mass as high as $4 \times 10^9\,M_{\odot}$) and evolved (age $\lesssim 500$ Myr) galaxies, while others are less massive ($M_{\text{stellar}}\sim 10^8\,M_{\odot}$) and very young ($\sim 10$ Myr) galaxies with strong nebular emission lines that boost their rest-frame optical fluxes. We identify two Ly$\alpha$ emitters in our sample from the Keck DEIMOS spectra, one at $z_{\text{Ly}\alpha}=6.76$ (in RXJ1347) and one at $z_{\text{Ly}\alpha}=6.32$ (in MACS0454). We show that IRAC $[3.6]-[4.5]$ color, when combined with photometric redshift, can be used to identify galaxies likely with strong nebular emission lines within certain redshift windows.

[abbreviated] We present a census of Ly\alpha\ emission at $z\gtrsim7$ utilizing deep near infrared HST grism spectroscopy from the first six completed clusters of the Grism Lens-Amplified Survey from Space (GLASS). In 24/159 photometrically selected galaxies we detect emission lines consistent with Ly\alpha\ in the GLASS spectra. Based on the distribution of signal-to-noise ratios and on simulations we expect the completeness and the purity of the sample to be 40-100% and 60-90%, respectively. For the objects without detected emission lines we show that the observed (not corrected for lensing magnification) 1$\sigma$ flux limits reaches $5\times10^{-18}$erg/s/cm$^{2}$ per position angle over the full wavelength range of GLASS (0.8-1.7$\mu$m). Based on the conditional probability of Ly\alpha\ emission measured from the ground at $z\sim7$ we would have expected 12-18 Ly\alpha\ emitters. This is consistent with the number of detections, within the uncertainties, confirming the drop in Ly\alpha\ emission with respect to $z\sim6$. These candidates include a promising source at $z=8.1$. The spatial extent of Ly\alpha\ in a deep stack of the most convincing Ly\alpha\ emitters with $\langle z\rangle=7.2$ is consistent with that of the rest-frame UV continuum. Extended Ly$\alpha$ emission, if present, has a surface brightness below our detection limit, consistent with the properties of lower redshift comparison samples. From the stack we estimate upper limits on rest-frame UV emission line ratios and find $f_\textrm{CIV} / f_\textrm{Ly${\alpha}$} \lesssim 0.32$ and $f_\textrm{CIII]} / f_\textrm{Ly$\alpha$} \lesssim 0.23$ in good agreement with other values published in the literature.

We give an overview of the Grism Lens Amplified Survey from Space (GLASS), a large Hubble Space Telescope program aimed at obtaining grism spectroscopy of the fields of ten massive clusters of galaxies at redshift z=0.308-0.686, including the Hubble Frontier Fields (HFF). The Wide Field Camera 3 yields near infrared spectra of the cluster cores, covering the wavelength range 0.81-1.69mum through grisms G102 and G141, while the Advanced Camera for Surveys in parallel mode provides G800L spectra of the infall regions of the clusters. The WFC3 spectra are taken at two almost orthogonal position angles in order to minimize the effects of confusion. After summarizing the scientific drivers of GLASS, we describe the sample selection as well as the observing strategy and data processing pipeline. We then utilize MACSJ0717.5+3745, a HFF cluster and the first one observed by GLASS, to illustrate the data quality and the high-level data products. Each spectrum brighter than H_AB=23 is visually inspected by at least two co-authors and a redshift is measured when sufficient information is present in the spectra. Furthermore, we conducted a thorough search for emission lines through all the GLASS WFC3 spectra with the aim of measuring redshifts for sources with continuum fainter than H_AB=23. We provide a catalog of 139 emission-line based spectroscopic redshifts for extragalactic sources, including three new redshifts of multiple image systems (one probable, two tentative). In addition to the data itself we also release software tools that are helpful to navigate the data.

We present a strong and weak lensing reconstruction of the massive cluster Abell 2744, the first cluster for which deep Hubble Frontier Field (HFF) images and spectroscopy from the Grism Lens-Amplified Survey from Space (GLASS) are available. By performing a targeted search for emission lines in multiply imaged sources using the GLASS spectra, we obtain 5 high-confidence spectroscopic redshifts and 2 tentative ones. We confirm 1 strongly lensed system by detecting the same emission lines in all 3 multiple images. We also search for additional line emitters blindly and use the full GLASS spectroscopic catalog to test reliability of photometric redshifts for faint line emitters. We see a reasonable agreement between our photometric and spectroscopic redshift measurements, when including nebular emission in photometric redshift estimations. We introduce a stringent procedure to identify only secure multiple image sets based on colors, morphology, and spectroscopy. By combining 7 multiple image systems with secure spectroscopic redshifts (at 5 distinct redshift planes) with 18 multiple image systems with secure photometric redshifts, we reconstruct the gravitational potential of the cluster pixellated on an adaptive grid, using a total of 72 images. The resulting mass map is compared with a stellar mass map obtained from the deep Spitzer Frontier Fields data to study the relative distribution of stars and dark matter in the cluster. We find that the stellar to total mass ratio varies substantially across the cluster field, suggesting that stars do not trace exactly the total mass in this interacting system. The maps of convergence, shear, and magnification are made available in the standard HFF format.

We present spectroscopy and laser guide star adaptive optics (LGSAO) images of the doubly imaged lensed quasar SDSS J1206+4332. We revise the deflector redshift proposed previously to $z_{d}=0.745,$ and measure for the first time its velocity dispersion $\sigma=(290\pm30)$ km/s. The LGSAO data show the lensed quasar host galaxy stretching over the astroid caustic thus forming an extra pair of merging images, which was previously thought to be an unrelated galaxy in seeing limited data. Owing to the peculiar geometry, the lens acts as a natural coronagraph on the broad-line region of the quasar so that only narrow [O III] emission is found in the fold arc. We use the data to reconstruct the source structure and deflector potential, including nearby perturbers. We reconstruct the point-spread function (PSF) from the quasar images themselves, since no additional point source is present in the field of view. From gravitational lensing and stellar dynamics, we find the slope of the total mass density profile to be $\gamma^{\prime}=-\log\rho/\log r =1.93\pm0.09.$ We discuss the potential of SDSS J1206+4332 for measuring time delay distance (and thus H$_0$ and other cosmological parameters), or as a standard ruler, in combination with the time delay published by the COSMOGRAIL collaboration. We conclude that this system is very promising for cosmography. However, in order to achieve competitive precision and accuracy, an independent characterization of the PSF is needed. Spatially resolved kinematics of the deflector would reduce the uncertainties further. Both are within the reach of current observational facilities.

We present new results from our search for z~7 galaxies from deep spectroscopic observations of candidate z-dropouts in the CANDELS fields. Despite the extremely low flux limits achieved by our sensitive observations, only 2 galaxies have robust redshift identifications, one from its Lyalpha emission line at z=6.65, the other from its Lyman-break, i.e. the continuum discontinuity at the Lyalpha wavelength consistent with a redshift 6.42, but with no emission line. In addition, for 23 galaxies we present deep limits in the Lyalpha EW derived from the non detections in ultra-deep observations. Using this new data as well as previous samples, we assemble a total of 68 candidate z~7 galaxies with deep spectroscopic observations, of which 12 have a line detection. With this much enlarged sample we can place solid constraints on the declining fraction of Ly$\alpha$ emission in z~7 Lyman break galaxies compared to z~6, both for bright and faint galaxies. Applying a simple analytical model, we show that the present data favor a patchy reionization process rather than a smooth one.

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$).

SURFSUP is a joint Spitzer and HST Exploration Science program using 10 galaxy clusters as cosmic telescopes to study z >~ 7 galaxies at intrinsically lower luminosities, enabled by gravitational lensing, than blank field surveys of the same exposure time. Our main goal is to measure stellar masses and ages of these galaxies, which are the most likely sources of the ionizing photons that drive reionization. Accurate knowledge of the star formation density and star formation history at this epoch is necessary to determine whether these galaxies indeed reionized the universe. Determination of the stellar masses and ages requires measuring rest frame optical light, which only Spitzer can probe for sources at z >~ 7, for a large enough sample of typical galaxies. Our program consists of 550 hours of Spitzer/IRAC imaging covering 10 galaxy clusters with very well-known mass distributions, making them extremely precise cosmic telescopes. We combine our data with archival observations to obtain mosaics with ~30 hours exposure time in both 3.6$\mu$m and 4.5$\mu$m in the central 4 arcmin x 4 arcmin field and ~15 hours in the flanking fields. This results in 3-$\sigma$ sensitivity limits of ~26.6 and ~26.2AB magnitudes for the central field in the IRAC 3.6 and 4.5$\mu$m bands, respectively. To illustrate the survey strategy and characteristics we introduce the sample, present the details of the data reduction and demonstrate that these data are sufficient for in-depth studies of z >~ 7 sources (using a z=9.5 galaxy behind MACSJ1149.5+2223 as an example). For the first cluster of the survey (the Bullet Cluster) we have released all high-level data mosaics and IRAC empirical PSF models. In the future we plan to release these data products for the entire survey.

High precision cosmological distance measurements towards individual objects such as time delay gravitational lenses or type Ia supernovae are affected by weak lensing perturbations by galaxies and groups along the line of sight. In time delay gravitational lenses, "external convergence," kappa, can dominate the uncertainty in the inferred distances and hence cosmological parameters. In this paper we attempt to reconstruct kappa, due to line of sight structure, using a simple halo model. We use mock catalogues from the Millennium Simulation, and calibrate and compare our reconstructed P(kappa) to ray-traced kappa "truth" values; taking into account realistic observational uncertainties. We find that the reconstruction of kappa provides an improvement in precision of ~50% over galaxy number counts. We find that the lowest-kappa lines of sight have the best constrained P(kappa). In anticipation of large future samples of lenses, we find that selecting the third of the systems with the highest precision kappa estimates gives a sample of unbiased time delay distance measurements with just ~1% uncertainty due to line of sight external convergence effects. Photometric data are sufficient to pre-select the best-constrained lines of sight, and can be done before investment in light-curve monitoring. Conversely, we show that selecting lines of sight with high external shear could, with the reconstruction model presented, induce biases of up to 1% in time delay distance. We find that a major potential source of systematic error is uncertainty in the high mass end of the stellar mass-halo mass relation; this could introduce ~2% biases on the time-delay distance if completely ignored. We suggest areas for the improvement of this general analysis framework (including more sophisticated treatment of high mass structures) that should allow yet more accurate cosmological inferences to be made.

We present the first results of our spectroscopic follow-up of 6.5 < z < 10 candidate galaxies behind clusters of galaxies. We report the spectroscopic confirmation of an intrinsically faint Lyman break galaxy (LBG) identified as a z 850LP-band dropout behind the Bullet Cluster. We detect an emission line at {\lambda} = 9412 {\AA} at >5{\sigma} significance using a 16 hr long exposure with FORS2 VLT. Based on the absence of flux in bluer broadband filters, the blue color of the source, and the absence of additional lines, we identify the line as Ly{\alpha} at z = 6.740 \pm 0.003. The integrated line flux is f = (0.7 \pm 0.1 \pm 0.3) \times 10^{-17} erg^{-1} s^{-1} cm^{-2} (the uncertainties are due to random and flux calibration errors, respectively) making it the faintest Ly{\alpha} flux detected at these redshifts. Given the magnification of {\mu} = 3.0 \pm 0.2 the intrinsic (corrected for lensing) flux is f^int = (0.23 \pm 0.03 \pm 0.10 \pm 0.02) \times 10^{-17} erg^{-1} s^{-1} cm^{-2} (additional uncertainty due to magnification), which is ~2-3 times fainter than other such measurements in z ~ 7 galaxies. The intrinsic H 160W-band magnitude of the object is m^int(H_160W)=27.57 \pm 0.17, corresponding to 0.5 L* for LBGs at these redshifts. The galaxy is one of the two sub-L* LBG galaxies spectroscopically confirmed at these high redshifts (the other is also a lensed z = 7.045 galaxy), making it a valuable probe for the neutral hydrogen fraction in the early universe.

We report spectroscopic confirmation and high-resolution infrared imaging of a z=2.79 triply-imaged galaxy behind the Bullet Cluster. This source, a Spitzer-selected luminous infrared galaxy (LIRG), is confirmed via polycyclic aromatic hydrocarbon (PAH) features using the Spitzer Infrared Spectrograph (IRS) and resolved with HST WFC3 imaging. In this galaxy, which with a stellar mass of M*=4e9 Msun is one of the two least massive ones studied with IRS at z>2, we also detect H_2 S(4) and H_2 S(5) pure rotational lines (at 3.1 sigma and 2.1 sigma) - the first detection of these molecular hydrogen lines in a high-redshift galaxy. From the molecular hydrogen lines we infer an excitation temperature T=377+68-84 K. The detection of these lines indicates that the warm molecular gas mass is 6(+36-4)% of the stellar mass and implies the likely existence of a substantial reservoir of cold molecular gas in the galaxy. Future spectral observations at longer wavelengths with facilities like the Herschel Space Observatory, the Large Millimeter Telescope, and the Atacama Pathfinder EXperiment (APEX) thus hold the promise of precisely determining the total molecular gas mass. Given the redshift, and using refined astrometric positions from the high resolution imaging, we also update the magnification estimate and derived fundamental physical properties of this system. The previously published values for total infrared luminosity, star formation rate, and dust temperature are confirmed modulo the revised magnification; however we find that PAH emission is roughly a factor of five stronger than would be predicted by the relations between the total infrared and PAH luminosity reported for SMGs and starbursts in Pope et al. (2008).

We measure the average mass properties of a sample of 41 strong gravitational lenses at moderate redshift (z ~ 0.4 - 0.9), and present the lens redshift for 6 of these galaxies for the first time. Using the techniques of strong and weak gravitational lensing on archival data obtained from the Hubble Space Telescope, we determine that the average mass overdensity profile of the lenses can be fit with a power-law profile (Delta_Sigma prop. to R^{-0.86 +/- 0.16}) that is within 1-sigma of an isothermal profile (Delta_Sigma prop. to R^{-1}) with velocity dispersion sigma_v = 260 +/- 20 km/s. Additionally, we use a two-component de Vaucouleurs+NFW model to disentangle the total mass profile into separate luminous and dark matter components, and determine the relative fraction of each component. We measure the average rest frame V-band stellar mass-to-light ratio (Upsilon_V = 4.0 +/- 0.6 h M_sol/L_sol) and virial mass-to-light ratio (tau_V = 300 +/- 90 h M_sol/L_sol) for our sample, resulting in a virial-to-stellar mass ratio of M_vir/M_* = 75 +/- 25. Finally, we compare our results to a previous study using low redshift lenses, to understand how galaxy mass profiles evolve over time. We investigate the evolution of M_vir/M_*(z) = alpha(1+z)^{beta}, and find best fit parameters of alpha = 51 +/- 36 and beta = 0.9 +/- 1.8, constraining the growth of virial to stellar mass ratio over the last ~7 Gigayears. We note that, by using a sample of strong lenses, we are able to constrain the growth of M_vir/M_*(z) without making any assumptions about the IMF of the stellar population.

We present a tomographic cosmological weak lensing analysis of the HST COSMOS Survey. Applying our lensing-optimized data reduction, principal component interpolation for the ACS PSF, and improved modelling of charge-transfer inefficiency, we measure a lensing signal which is consistent with pure gravitational modes and no significant shape systematics. We carefully estimate the statistical uncertainty from simulated COSMOS-like fields obtained from ray-tracing through the Millennium Simulation. We test our pipeline on simulated space-based data, recalibrate non-linear power spectrum corrections using the ray-tracing, employ photometric redshifts to reduce potential contamination by intrinsic galaxy alignments, and marginalize over systematic uncertainties. We find that the lensing signal scales with redshift as expected from General Relativity for a concordance LCDM cosmology, including the full cross-correlations between different redshift bins. For a flat LCDM cosmology, we measure sigma_8(Omega_m/0.3)^0.51=0.75+-0.08 from lensing, in perfect agreement with WMAP-5, yielding joint constraints Omega_m=0.266+0.025-0.023, sigma_8=0.802+0.028-0.029 (all 68% conf.). Dropping the assumption of flatness and using HST Key Project and BBN priors only, we find a negative deceleration parameter q_0 at 94.3% conf. from the tomographic lensing analysis, providing independent evidence for the accelerated expansion of the Universe. For a flat wCDM cosmology and prior w in [-2,0], we obtain w<-0.41 (90% conf.). Our dark energy constraints are still relatively weak solely due to the limited area of COSMOS. However, they provide an important demonstration for the usefulness of tomographic weak lensing measurements from space. (abridged)

The gravitational potential of clusters of galaxies acts as a cosmic telescope allowing us to find and study galaxies at fainter limits than otherwise possible and thus probe closer to the epoch of formation of the first galaxies. We use the Bullet Cluster 1E0657-56 (z = 0.296) as a case study, because its high mass and merging configuration makes it one of the most efficient cosmic telescopes we know. We develop a new algorithm to reconstruct the gravitational potential of the Bullet Cluster, based on a non-uniform adaptive grid, combining strong and weak gravitational lensing data derived from deep HST/ACS F606W-F775W-F850LP and ground-based imaging. We exploit this improved mass map to study z~5-6 Lyman Break Galaxies (LBGs), which we detect as dropouts. One of the LBGs is multiply imaged, providing a geometric confirmation of its high redshift, and is used to further improve our mass model. We quantify the uncertainties in the magnification map reconstruction in the intrinsic source luminosity, and in the volume surveyed, and show that they are negligible compared to sample variance when determining the luminosity function of high-redshift galaxies. With shallower and comparable magnitude limits to HUDF and GOODS, the Bullet cluster observations, after correcting for magnification, probe deeper into the luminosity function of the high redshift galaxies than GOODS and only slightly shallower than HUDF. We conclude that accurately focused cosmic telescopes are the most efficient way to sample the bright end of the luminosity function of high redshift galaxies and - in case they are multiply imaged - confirm their redshifts.

We constrain the physical nature of dark matter using the newly identified massive merging galaxy cluster MACSJ0025.4-1222. As was previously shown by the example of the Bullet Cluster (1E0657-56), such systems are ideal laboratories for detecting isolated dark matter, and distinguishing between cold dark matter (CDM) and other scenarios (e.g. self-interacting dark matter, alternative gravity theories). MACSJ0025.4-1222 consists of two merging subclusters of similar richness at z=0.586. We measure the distribution of X-ray emitting gas from Chandra X-ray data and find it to be clearly displaced from the distribution of galaxies. A strong (information from highly distorted arcs) and weak (using weakly distorted background galaxies) gravitational lensing analysis based on Hubble Space Telescope observations and Keck arc spectroscopy confirms that the subclusters have near-equal mass. The total mass distribution in each of the subclusters is clearly offset (at >4sigma significance) from the peak of the hot X-ray emitting gas (the main baryonic component), but aligned with the distribution of galaxies. We measure the fractions of mass in hot gas (0.09^{+0.07}_{-0.03}) and stars (0.010^{+0.007}_{-0.004}), consistent with those of typical clusters, finding that dark matter is the dominant contributor to the gravitational field. Under the assumption that the subclusters experienced a head-on collision in the plane of the sky, we obtain an order-of-magnitude estimate of the dark matter self-interaction cross-section of sigma/m < 4cm^2/g, re-affirming the results from the Bullet Cluster on the collisionless nature of dark matter.

The galaxy cluster RX J1347-1145 is one of the most X-ray luminous and most massive clusters known. Its extreme mass makes it a prime target for studying issues addressing cluster formation and cosmology. In this paper we present new high-resolution HST/ACS and Chandra X-ray data. The high resolution and sensitivity of ACS enabled us to detect and quantify several new multiply imaged sources, we now use a total of eight for the strong lensing analysis. Combining this information with shape measurements of weak lensing sources in the central regions of the cluster, we derive a high-resolution, absolutely-calibrated mass map. This map provides the best available quantification of the total mass of the central part of the cluster to date. We compare the reconstructed mass with that inferred from the new Chandra X-ray data, and conclude that both mass estimates agree extremely well in the observed region, namely within 400 / h_70 kpc of the cluster center. In addition we study the major baryonic components (gas and stars) and hence derive the dark matter distribution in the center of the cluster. We find that the dark matter and baryons are both centered on the BCG within the uncertainties (alignment is better than <10 kpc). We measure the corresponding 1-D profiles and find that dark matter distribution is consistent with both NFW and cored profiles, indicating that a more extended radial analysis is needed to pinpoint the concentration parameter, and hence the inner slope of the dark matter profile.