Accurate weak lensing mass estimates of clusters are needed in order to calibrate mass proxies for the cosmological exploitation of galaxy cluster surveys. Such measurements require accurate knowledge of the redshift distribution of the weak lensing source galaxies. In this context, we investigate the accuracy of photometric redshifts (photo-$z$s) computed by the 3D-HST team for the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey fields, which provide a relevant photometric reference data set for deep weak lensing studies. Through the comparison to spectroscopic redshifts and photo-$z$s based on very deep data from the Hubble Ultra Deep Field, we identify catastrophic redshift outliers in the 3D-HST/CANDELS catalogue. These would significantly bias weak lensing results if not accounted for. We investigate the cause of these outliers and demonstrate that the interpolation of spectral energy distribution (SED) templates and a well-selected combination of photometric data can reduce the net impact for weak lensing studies.

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.

The thermodynamic properties of the hot plasma in galaxy clusters retains information on the processes leading to the formation and evolution of the gas in their deep, dark matter potential wells. These processes are dictated not only by gravity but also by gas physics, e.g. AGN feedback and turbulence. In this work, we study the thermodynamic properties, e.g. density, temperature, pressure, and entropy, of the most massive and the most distant ($z > 1.2$) SPT-selected clusters, and compare them with those of the nearby clusters ($z<0.1$) to constrain their evolution as a function of time and radius. We find that thermodynamic properties in the outskirts of high redshift clusters are remarkably similar to the low redshift clusters, and their evolution follows the prediction of the self-similar model. Their intrinsic scatter is larger, indicating that the physical properties that lead to the formation and virialization of cluster outskirts show evolving variance. On the other hand, thermodynamic properties in the cluster cores deviates significantly from self-similarity indicating that the processes that regulate the core are already in place in these very high redshift clusters. This result is supported by the unevolving physical scatter of all thermodynamic quantities in cluster cores.

The Lya emitter (LAE) fraction, X_LAE, is a potentially powerful probe of the evolution of the intergalactic neutral hydrogen gas fraction. However, uncertainties in the measurement of X_LAE are still debated. Thanks to deep data obtained with MUSE, we can measure the evolution of X_LAE homogeneously over a wide redshift range of z~3-6 for UV-faint galaxies (down to M_1500~-17.75). This is significantly fainter than in former studies, and allows us to probe the bulk of the population of high-z star-forming galaxies. We construct a UV-complete photo-redshift sample following UV luminosity functions and measure the Lya emission with MUSE using the second data release from the MUSE HUDF Survey. We derive the redshift evolution of X_LAE for M_1500 in [-21.75;-17.75] for the first time with a equivalent width range EW(Lya)>=65 A and find low values of X_ LAE<~30% at z<~6. For M_1500 in [-20.25;-18.75] and EW(Lya)<~25 A, our X_LAE values are consistent with those in the literature within 1sigma at z<~5, but our median values are systematically lower than reported values over the whole redshift range. In addition, we do not find a significant dependence of X_LAE on M_1500 for EW(Lya)>~50 A at z~3-4, in contrast with previous work. The differences in X_LAE mainly arise from selection biases for Lyman Break Galaxies (LBGs) in the literature: UV-faint LBGs are more easily selected if they have strong Lya emission, hence X_LAE is biased towards higher values. Our results suggest either a lower increase of X_LAE towards z~6 than previously suggested, or even a turnover of X_LAE at z~5.5, which may be the signature of a late or patchy reionization process. We compared our results with predictions from a cosmological galaxy evolution model. We find that a model with a bursty star formation (SF) can reproduce our observed X_LAE much better than models where SF is a smooth function of time.

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

While low-luminosity galaxies dominate number counts at all redshifts, their contribution to cosmic Reionization is poorly understood due to a lack of knowledge of their physical properties. We isolate a sample of 35 z~4-5 continuum-faint Lyman-alpha emitters from deep VLT/MUSE spectroscopy and directly measure their Halpha emission using stacked Spitzer/IRAC Ch. 1 photometry. Based on Hubble Space Telescope imaging, we determine that the average UV continuum magnitude is fainter than -16 (~0.01 L_star), implying a median Lyman-alpha equivalent width of 249 Angstroms. By combining the Halpha measurement with the UV magnitude we determine the ionizing photon production efficiency, xi_ion, a first for such faint galaxies. The measurement of log (xi_ion [Hz/erg]) = 26.28 (+0.28; -0.40) is in excess of literature measurements of both continuum- and emission line-selected samples, implying a more efficient production of ionizing photons in these lower-luminosity, Lyman-alpha-selected systems. We conclude that this elevated efficiency can be explained by stellar populations with metallicities between 4e-4 and 0.008, with light-weighted ages less than 3 Myr.

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

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 present the concept of BlueMUSE, a blue-optimised, medium spectral resolution, panoramic integral field spectrograph based on the MUSE concept and proposed for the Very Large Telescope. With an optimised transmission down to 350 nm, a larger FoV (1.4 x 1.4 arcmin$^2$) and a higher spectral resolution compared to MUSE, BlueMUSE will open up a new range of galactic and extragalactic science cases allowed by its specific capabilities, beyond those possible with MUSE. For example a survey of massive stars in our galaxy and the Local Group will increase the known population of massive stars by a factor $>$100, to answer key questions about their evolution. Deep field observations with BlueMUSE will also significantly increase samples of Lyman-alpha emitters, spanning the era of Cosmic Noon. This will revolutionise the study of the distant Universe: allowing the intergalactic medium to be detected unambiguously in emission, enabling the study of the exchange of baryons between galaxies and their surroundings. By 2030, at a time when the focus of most of the new large facilities (ELT, JWST) will be on the infra-red, BlueMUSE will be a unique facility, outperforming any ELT instrument in the Blue/UV. It will have a strong synergy with ELT, JWST as well as ALMA, SKA, Euclid and Athena.

We present the first measurements of the Lyman-continuum photon production efficiency $\xi_{\textrm{ion,0}}$ at $z\sim4$-5 for galaxies fainter than 0.2 $L^*$ ($-$19 mag). $\xi_{\textrm{ion,0}}$ quantifies the production rate of ionizing photons with respect to the UV luminosity density assuming a fiducial escape fraction of zero. Extending previous measurements of $\xi_{\textrm{ion,0}}$ to the faint population is important, as ultra-faint galaxies are expected to contribute the bulk of the ionizing emissivity. We probe $\xi_{\textrm{ion,0}}$ to such faint magnitudes by taking advantage of 200-hour depth Spitzer/IRAC observations from the GREATS program and $\approx$300 3<$z$<6 galaxies with spectroscopic redshifts from the MUSE GTO Deep + Wide programs. Stacked IRAC [3.6]$-$[4.5] colors are derived and used to infer the H$\alpha$ rest-frame equivalent widths, which range from 403\r{A} to 2818\r{A}. The derived $\xi_{\textrm{ion,0}}$ is $\log_{10}(\xi_{\textrm{ion,0}} / \textrm{Hz erg}^{-1}) = 25.36 \pm 0.08$ over $-$20.5 < M$_{\textrm{UV}}$ < $-$17.5, similar to those derived for brighter galaxy samples at the same redshift and therefore suggesting that $\xi_{\textrm{ion}}$ shows no strong dependence on $M_{UV}$. The $\xi_{\textrm{ion,0}}$ values found in our sample imply that the Lyman-continuum escape fraction for $M_{\textrm{UV}} \approx -19$ star-forming galaxies cannot exceed $\approx$8-20\% in the reionization era.

We present the results obtained with VLT/MUSE on the faint-end of the Lyman-alpha luminosity function (LF) based on deep observations of four lensing clusters. The precise aim of the present study is to further constrain the abundance of Lyman-alpha emitters (LAEs) by taking advantage of the magnification provided by lensing clusters. We blindly selected a sample of 156 LAEs, with redshifts between $2.9 \le z \le 6.7$ and magnification-corrected luminosities in the range $ 39 \lesssim \log L_{Ly_{\alpha}}$ [erg s$^{-1}$] $\lesssim 43$. The price to pay to benefit from magnification is a reduction of the effective volume of the survey, together with a more complex analysis procedure. To properly take into account the individual differences in detection conditions (including lensing configurations, spatial and spectral morphologies) when computing the LF, a new method based on the 1/Vmax approach was implemented. The LAE LF has been obtained in four different redshift bins with constraints down to $\log L_{Ly_{\alpha}} = 40.5$. From our data only, no significant evolution of LF mean slope can be found. When performing a Schechter analysis including data from the literature to complete the present sample a steep faint-end slope was measured varying from $\alpha = -1.69^{+0.08}_{-0.08}$ to $\alpha = -1.87^{+0.12}_{-0.12}$ between the lowest and the highest redshift bins. The contribution of the LAE population to the star formation rate density at $z \sim 6$ is $\lesssim 50$% depending on the luminosity limit considered, which is of the same order as the Lyman-break galaxy (LBG) contribution. The evolution of the LAE contribution with redshift depends on the assumed escape fraction of Lyman-alpha photons, and appears to slightly increase with increasing redshift when this fraction is conservatively set to one. (abridged)

The existence of a spatially resolved Star-Forming Main Sequence (rSFMS) and a spatially resolved Mass-Metallicity Relation (rMZR) is now well established for local galaxies. Moreover, gradients with metallicity decreasing with radius seem to be common in local disc galaxies. These observations suggest that galaxy formation is a self-regulating process, and provide constraints for galaxy evolution models. Studying the evolution of these relations at higher redshifts is still however very challenging. In this paper, we analyse three gravitationally lensed galaxies at z = 0.6, 0.7 and 1, observed with MUSE and SINFONI. These galaxies are highly magnified by galaxy clusters, which allow us to observe resolved scaling relations and metallicity gradients on physical scales of a couple of hundred parsecs, comparable to studies of local galaxies. We confirm that the rSFMS is already in place at these redshifts on sub-kpc scales, and establish, for the first time, the existence of the rMZR at higher redshifts. We develop a forward-modelling approach to fit 2D metallicity gradients of multiply imaged lensed galaxies in the image plane, and derive gradients of -0.027+/-0.003, -0.019+/-0.003 and -0.039+/-0.060 dex/kpc. Despite the fact that these are clumpy galaxies, typical of high redshift discs, the metallicity variations in the galaxies are well described by global linear gradients, and we do not see any difference in metallicity associated with the star-forming clumps.

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

We 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 an updated strong-lensing analysis of the massive cluster Abell 370 (A370), continuing the work first presented in Lagattuta et al. (2017). In this new analysis, we take advantage of the deeper imaging data from the Hubble Space Telescope (HST) Frontier Fields program, as well as a large spectroscopic mosaic obtained with the Multi-Unit Spectroscopic Explorer (MUSE). Thanks to the extended coverage of this mosaic, we probe the full 3D distribution of galaxies in the field, giving us a unique picture of the extended structure of the cluster and its surroundings. Our final catalog contains 584 redshifts, representing the largest spectroscopic catalog of A370 to date. Constructing the model, we measure a total mass distribution that is quantitatively similar to our previous work -- though to ensure a low rms error in the model fit, we invoke a significantly large external shear term. Using the redshift catalog, we search for other bound groups of galaxies, which may give rise to a more physical interpretation of this shear. We identify three structures in narrow redshift ranges along the line of sight, highlighting possible infalling substructures into the main cluster halo. We also discover additional substructure candidates in low-resolution imaging at larger projected radii. More spectroscopic coverage of these regions (pushing close to the A370 virial radius) and more extended, high-resolution imaging will be required to investigate this possibility, further advancing the analysis of these interesting developments.

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

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.

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 report that LAMOST-HVS1 is a massive hyper-runaway subgiant star with mass of 8.3 Msun and super-Solar metallicity, ejected from the inner stellar disk of the Milky Way $\sim$ 33 Myr ago with the intrinsic ejection velocity of $568^{+19}_{-17}$ km/s (corrected for the streaming motion of the disk), based on the proper motion data from Gaia Data Release 2 (DR2) and high-resolution spectroscopy. The extremely large ejection velocity indicates that this star was not ejected by the supernova explosion of the binary companion. Rather, it was probably ejected by a 3- or 4-body dynamical interaction with more massive objects in a high-density environment. Such a high-density environment may be attained at the core region of a young massive cluster with mass of $\gtrsim 10^4$ Msun. The ejection agent that took part in the ejection of LAMOST-HVS1 may well be an intermediate mass black hole ($\gtrsim$ 100 Msun), a very massive star ($\gtrsim$ 100 Msun), or multiple ordinary massive stars ($\gtrsim$ 30 Msun). Based on the flight time and the ejection location of LAMOST-HVS1, we argue that its ejection agent or its natal star cluster is currently located near the Norma spiral arm. The natal star cluster of LAMOST-HVS1 may be an undiscovered young massive cluster near the Norma spiral arm.

Measuring time delays from strongly lensed supernovae (SNe) is emerging as a novel and independent tool for estimating the Hubble constant $(H_0)$. This is very important given the recent discord in the value of $H_0$ from two methods that probe different distance ranges. The success of this technique will rely of our ability to discover strongly lensed SNe with measurable time delays. Here, we present the magnifications and the time delays for the multiply-imaged galaxies behind the Hubble Frontier Fields (HFF) galaxy clusters, by using recently published lensing models. Continuing on our previous work done for Abell 1689 (A1689) and Abell 370, we also show the prospects of observing strongly lensed SNe behind the HFF clusters with the upcoming James Webb Space Telescope (JWST). With four 1-hour visits in one year, the summed expectations of all six HFF clusters are $\sim0.5$ core-collapse (CC) SNe and $\sim0.06$ Type Ia SNe (SNe Ia) in F115W band, while with F150W the expectations are higher, $\sim0.9$ CC SNe and $\sim0.06$ SNe Ia. These estimates match those expected by only surveying A1689, proving that the performance of A1689 as gravitational telescope is superior. In the five HFF clusters presented here, we find that F150W will be able to detect SNe Ia (SNe IIP) exploding in 93 (80) pairs multiply-imaged galaxies with time delays of less than 5 years.

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

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

The physical origin of the near-ultraviolet MgII emission remains an under-explored domain, contrary to more typical emission lines detected in the spectra of star-forming galaxies. We explore the nebular and physical properties for a sample of 381 galaxies between 0.70 < z < 2.34 drawn from the MUSE Hubble Ultra Deep Survey. The spectra of these galaxies show a wide variety of profiles of the MgII 2796,2803 resonant doublet, from absorption to emission. We present a study on the main drivers for the detection of MgII emission in galaxy spectra. By exploiting photoionization models we verified that the emission-line ratios observed in galaxies with MgII in emission are consistent with nebular emission from HII regions. From a simultaneous analysis of MUSE spectra and ancillary HST information via spectral energy distribution (SED) fitting, we find that galaxies with MgII in emission have lower stellar masses, smaller sizes, bluer spectral slopes and lower optical depth than those with absorption. This leads us to suggest that MgII emission is a potential tracer of physical conditions not merely related to those of the ionized gas. We show that these differences in MgII emission/absorption can be explained in terms of a higher dust and neutral gas content in the interstellar medium (ISM) of galaxies showing MgII in absorption, confirming the extreme sensitivity of MgII to the presence of the neutral ISM. We conclude with an analogy between the MgII doublet and the Ly-alpha line, due to their resonant nature. Further investigations with current and future facilities, including JWST, are promising as the detection of MgII emission and its potential connection with Ly-alpha could provide new insights on the ISM content in the early Universe.

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 analyse a sample of 8 highly magnified galaxies at redshift 0.6<z<1.5 observed with MUSE, exploring the resolved properties of these galaxies at sub-kiloparsec scales. Combining multi-band HST photometry and MUSE spectra, we derive the stellar mass, global star formation rates, extinction and metallicity from multiple nebular lines, concluding that our sample is representative of z~1 star-forming galaxies. We derive the 2D kinematics of these galaxies from the [OII] emission and model it with a new method that accounts for lensing effects and fits multiple images simultaneously. We use these models to calculate the 2D beam-smearing correction and derive intrinsic velocity dispersion maps. We find them to be fairly homogeneous, with relatively constant velocity dispersions between 15 - 80 km/s and Gini coefficient of <0.3. We do not find any evidence for higher (or lower) velocity dispersions at the positions of bright star-forming clumps. We derive resolved maps of dust attenuation and attenuation-corrected star formation rates from emission lines for two objects in the sample. We use this information to study the relation between resolved star formation rate and velocity dispersion. We find that these quantities are not correlated, and the high-velocity dispersions found for relatively low star-forming densities seems to indicate that, at sub-kiloparsec scales, turbulence in high-z discs is mainly dominated by gravitational instability rather than stellar feedback.

We present the deepest study to date of the Lya luminosity function (LF) in a blank field using blind integral field spectroscopy from MUSE. We constructed a sample of 604 Lya emitters (LAEs) across the redshift range 2.91 < z < 6.64 using automatic detection software in the Hubble Ultra Deep Field. We calculate accurate total Lya fluxes capturing low surface brightness extended Lya emission now known to be a generic property of high-redshift star-forming galaxies. We simulated realistic extended LAEs to characterise the selection function of our samples, and performed flux-recovery experiments to test and correct for bias in our determination of total Lya fluxes. We find an accurate completeness correction accounting for extended emission reveals a very steep faint-end slope of the LF, alpha, down to luminosities of log10 L erg s^-1< 41.5, applying both the 1/Vmax and maximum likelihood estimators. Splitting the sample into three broad redshift bins, we see the faint-end slope increasing from -2.03+1.42-inf at z ~ 3.44 to -2.86+0.76-inf at z ~ 5.48, however no strong evolution is seen between the 68% confidence regions in L*-alpha parameter space. Using the Lya line flux as a proxy for star formation activity, and integrating the observed LFs, we find that LAEs' contribution to the cosmic SFRD rises with redshift until it is comparable to that from continuum-selected samples by z ~ 6. This implies that LAEs may contribute more to the star-formation activity of the early Universe than previously thought - any additional interglactic medium correction would act to further boost the Lya luminosities. Finally, assuming fiducial values for the escape of Lya and LyC radiation, and the clumpiness of the IGM, we integrated the maximum likelihood LF at 5.00 < z < 6.64 and find we require only a small extrapolation beyond the data (< 1 dex in L) for LAEs alone to maintain an ionised IGM at z ~ 6.

We present an analysis of MUSE observations obtained on the massive Frontier Fields cluster Abell 2744. This new dataset covers the entire multiply-imaged region around the cluster core. We measure spectroscopic redshifts for HST-selected continuum sources together with line emitters blindly detected in the datacube. The combined catalog consists of 514 spectroscopic redshifts (with 414 new identifications), including 156 cluster members and 326 magnified background sources. We use this redshift information to perform a strong-lensing analysis of all multiple images previously found in the deep Frontier Field images, and add three new MUSE-detected multiply-imaged systems with no obvious HST counterpart. The combined strong lensing constraints include a total of 60 systems producing 188 images altogether, out of which 29 systems and 83 images are spectroscopically confirmed, making Abell 2744 one of the most well-constrained clusters to date. A parametric mass model including two cluster-scale components in the core and several group-scale substructures at larger radii accurately reproduces all the spectroscopic multiple systems, reaching an rms of 0.67" in the image plane. Overall, the large number of spectroscopic redshifts gives us a robust model and we estimate the systematics on the mass density and magnification within the cluster core to be typically ~9%.

Strong gravitational lensing by galaxy clusters is a fundamental tool to study dark matter and constrain the geometry of the Universe. Recently, the Hubble Space Telescope Frontier Fields programme has allowed a significant improvement of mass and magnification measurements but lensing models still have a residual root mean square between 0.2 arcsec and few arcsec- onds, not yet completely understood. Systematic errors have to be better understood and treated in order to use strong lensing clusters as reliable cosmological probes. We have analysed two simulated Hubble-Frontier-Fields-like clusters from the Hubble Frontier Fields Comparison Challenge, Ares and Hera. We use several estimators (relative bias on magnification, den- sity profiles, ellipticity and orientation) to quantify the goodness of our reconstructions by comparing our multiple models, optimized with the parametric software LENSTOOL , with the input models. We have quantified the impact of systematic errors arising, first, from the choice of different density profiles and configurations and, secondly, from the availability of con- straints (spectroscopic or photometric redshifts, redshift ranges of the background sources) in the parametric modelling of strong lensing galaxy clusters and therefore on the retrieval of cosmological parameters. We find that substructures in the outskirts have a significant im- pact on the position of the multiple images, yielding tighter cosmological contours. The need for wide-field imaging around massive clusters is thus reinforced. We show that competitive cosmological constraints can be obtained also with complex multimodal clusters and that photometric redshifts improve the constraints on cosmological parameters when considering a narrow range of (spectroscopic) redshifts for the sources.

We present a strong lensing analysis on the massive cluster Abell 370 (A370; z = 0.375), using a combination of deep multi-band Hubble Space Telescope (HST) imaging and Multi-Unit Spectroscopic Explorer (MUSE) spectroscopy. From only two hours of MUSE data, we are able to measure 120 redshifts in the Southern BCG area, including several multiply-imaged lens systems. In total, we increase the number of multiply-imaged systems with a secure redshift from 4 to 15, nine of which are newly discovered. Of these, eight are located at z > 3, greatly extending the redshift range of spectroscopically-confirmed systems over previous work. Using these systems as constraints, we update a parametric lens model of A370, probing the mass distribution from cluster to galaxy scales. Overall, we find that a model with only two cluster- scale dark matter halos (one for each BCG) does a poor job of fitting these new image constraints. Instead, two additional mass clumps -- a central "bar" of mass located between the BCGs, and another clump located within a "crown" of galaxies in the Northern part of the cluster field -- provide significant improvements to the fit. Additional physical evidence suggests these clumps are indeed real features of the system, but with relatively few image constraints in the crown region, this claim is difficult to evaluate from a modeling perspective. Additional MUSE observations of A370 covering the entire strong-lensing region will greatly help these efforts, further improving our understanding of this intriguing cluster.

We present a high-precision mass model of the galaxy cluster MACSJ1149.6+2223, based on a strong-gravitational-lensing analysis of Hubble Space Telescope Frontier Fields (HFF) imaging data and spectroscopic follow-up with Gemini/GMOS and VLT/MUSE. Our model includes 12 new multiply imaged galaxies, bringing the total to 22, comprised of 65 individual lensed images. Unlike the first two HFF clusters, Abell 2744 and MACSJ0416.1-2403, MACSJ1149 does not reveal as many multiple images in the HFF data. Using the Lenstool software package and the new sets of multiple images, we model the cluster with several cluster-scale dark-matter halos and additional galaxy-scale halos for the cluster members. Consistent with previous analyses, we find the system to be complex, composed of five cluster-scale halos. Their spatial distribution and lower mass, however, makes MACSJ1149 a less powerful lens. Our best-fit model predicts image positions with an RMS of 0.91". We measure the total projected mass inside a 200~kpc aperture as ($1.840\pm 0.006$)$\times 10^{14}$M$_{\odot}$, thus reaching again 1% precision, following our previous HFF analyses of MACSJ0416.1-2403 and Abell 2744. In light of the discovery of the first resolved quadruply lensed supernova, SN Refsdal, in one of the multiply imaged galaxies identified in MACSJ1149, we use our revised mass model to investigate the time delays and predict the rise of the next image between November 2015 and January 2016.