Context. Spectroscopic surveys of massive galaxy clusters reveal the properties of faint background galaxies thanks to the magnification provided by strong gravitational lensing. Aims. We present a systematic analysis of integral-field-spectroscopy observations of 12 massive clusters, conducted with the Multi Unit Spectroscopic Explorer (MUSE). All data were taken under very good seeing conditions (∼0″​​.6) in effective exposure times between two and 15 h per pointing, for a total of 125 h. Our observations cover a total solid angle of ∼23 arcmin2 in the direction of clusters, many of which were previously studied by the MAssive Clusters Survey, Frontier Fields (FFs), Grism Lens-Amplified Survey from Space and Cluster Lensing And Supernova survey with Hubble programmes. The achieved emission line detection limit at 5σ for a point source varies between (0.77–1.5) × 10−18 erg s−1 cm−2 at 7000 Å. Methods. We present our developed strategy to reduce these observational data, detect continuum sources and line emitters in the datacubes, and determine their redshifts. We constructed robust mass models for each cluster to further confirm our redshift measurements using strong-lensing constraints, and identified a total of 312 strongly lensed sources producing 939 multiple images. Results. The final redshift catalogues contain more than 3300 robust redshifts, of which 40% are for cluster members and ∼30% are for lensed Lyman-α emitters. Fourteen percent of all sources are line emitters that are not seen in the available HST images, even at the depth of the FFs (∼29 AB). We find that the magnification distribution of the lensed sources in the high-magnification regime (μ = 2–25) follows the theoretical expectation of N(z) ∝ μ−2. The quality of this dataset, number of lensed sources, and number of strong-lensing constraints enables detailed studies of the physical properties of both the lensing cluster and the background galaxies. The full data products from this work, including the datacubes, catalogues, extracted spectra, ancillary images, and mass models, are made available to the community.

We present spatially resolved maps of six individually-detected Lyman α haloes (LAHs) as well as a first statistical analysis of the Lyman α (Lyα) spectral signature in the circum-galactic medium of high-redshift star-forming galaxies (−17.5 >  MUV >  −21.5) using the Multi-Unit Spectroscopic Explorer. Our resolved spectroscopic analysis of the LAHs reveals significant intrahalo variations of the Lyα line profile. Using a three-dimensional two-component model for the Lyα emission, we measured the full width at half maximum (FWHM), the peak velocity shift, and the asymmetry of the Lyα line in the core and in the halo of 19 galaxies. We find that the Lyα line shape is statistically different in the halo compared to the core (in terms of width, peak wavelength, and asymmetry) for ≈40% of our galaxies. Similarly to object-by-object based studies and a recent resolved study using lensing, we find a correlation between the peak velocity shift and the width of the Lyα line both at the interstellar and circum-galactic scales. This trend has been predicted by radiative transfer simulations of galactic winds as a result of resonant scattering in outflows. While there is a lack of correlation between the spectral properties and the spatial scale lengths of our LAHs, we find a correlation between the width of the line in the LAH and the halo flux fraction. Interestingly, UV bright galaxies (MUV <  −20) show broader, more redshifted, and less asymmetric Lyα lines in their haloes. The most significant correlation found is for the FWHM of the line and the UV continuum slope of the galaxy, suggesting that the redder galaxies have broader Lyα lines. The generally broad and red line shapes found in the halo component suggest that the Lyα haloes are powered either by scattering processes through an outflowing medium, fluorescent emission from outflowing cold clumps of gas, or a mix of both. Considering the large diversity of the Lyα line profiles observed in our sample and the lack of strong correlation, the interpretation of our results is still broadly open and underlines the need for realistic spatially resolved models of the LAHs.

Exploiting the fundamentally achromatic nature of gravitational lensing, we present a lens model for the massive galaxy cluster SMACS J0723.3−7323 (SMACS J0723, 𝑧 = 0.388) that significantly improves upon earlier work. Building on strong-lensing constraints identified in prior Hubble Space Telescope (HST) observations, the mass model utilizes 21 multiple-image systems, 16 of which were newly discovered in Early Release Observation (ERO) data from the James Webb Space Telescope (JWST). The resulting lens model maps the cluster mass distribution to an RMS spatial precision of 1. 08 and is publicly available a) . Consistent with previous analyses, our study shows SMACS J0723.3−7323 to be well described by a single large-scale component centered on the location of the brightest cluster galaxy, however JWST data point to the need for two additional diffuse components west of the cluster, which in turn accounts for all the currently identified multiply imaged systems. A comparison of the galaxy distribution, the mass distribution, and gas distribution in the core of SMACS J0723 based on HST, JWST, and Chandra data reveals a fairly concentrated regular elliptical profile along with with tell-tale signs of recent merger activity, possibly proceeding aligned closely to our line of sight. The exquisite sensitivity of JWST's NIRCAM reveals in spectacular fashion both the extended intra-cluster-light distribution and numerous star-forming clumps in magnified background galaxies. The high-precision lens model derived here for SMACS J0723−7323 demonstrates impressively the power of combining HST and JWST data for unprecedented studies of structure formation and evolution in the distant Universe.

We present an analysis of H I Lyman-α emission in deep VLT/MUSE observations of two highly magnified and extended galaxies at z = 3.5 and 4.03, including a newly discovered, almost complete Einstein ring. While these Lyman-α haloes are intrinsically similar to the ones typically seen in other MUSE deep fields, the benefits of gravitational lensing allows us to construct exceptionally detailed maps of Lyman-α line properties at sub-kpc scales. By combining all multiple images, we are able to observe complex structures in the Lyman-α emission and uncover small (∼ 120 km s −1 in Lyman-α peak shift), but significant at > 4 σ, systematic variations in the shape of the Lyman-α line profile within each halo. Indeed, we observe a global trend for the line peak shift to become redder at large radii, together with a strong correlation between the peak wavelength and line width. This systematic intrahalo variation is markedly similar to the object-to-object variations obtained from the integrated properties of recent large samples. Regions of high surface brightness correspond to relatively small line shifts, which could indicate that Lyman-α emission escapes preferentially from regions where the line profile has been less severely affected by scattering of Lyman-α photons.

ABSTRACT Ly α emission lines are typically found to be redshifted with respect to the systemic redshifts of galaxies, likely due to resonant scattering of Ly α photons. Here, we measure the average velocity offset for a sample of 96 z ≈ 3.3 Ly α emitters (LAEs) with a median Ly α flux (luminosity) of ${\approx}10^{-17}~\rm erg~cm^{-2}~s^{-1}$ (${\approx}10^{42}~\rm erg~s^{-1}$) and a median star formation rate (SFR) of ≈1.3 $\rm M_{\odot }~yr^{-1}$ (not corrected for possible dust extinction), detected by the Multi-Unit Spectroscopic Explorer as part of our MUSEQuBES circumgalactic medium (CGM) survey. By postulating that the stacked CGM absorption profiles of these LAEs, probed by eight background quasars, must be centred on the systemic redshift, we measure an average velocity offset, Voffset = 171\pm 8 km s−1, between the Ly α emission peak and the systemic redshift. The observed Voffset is lower by factors of ≈1.4 and ≈2.6 compared to the velocity offsets measured for narrow-band-selected LAEs and Lyman break galaxies, respectively, which probe galaxies with higher masses and SFRs. Consistent with earlier studies based on direct measurements for individual objects, we find that the Voffset is correlated with the full width at half-maximum of the red peak of the Ly α line, and anticorrelated with the rest-frame equivalent width. Moreover, we find that Voffset is correlated with SFR with a sub-linear scaling relation, $V_{\rm offset}\propto \rm SFR^{0.16\pm 0.03}$. Adopting the mass scaling for main-sequence galaxies, such a relation suggests that Voffset scales with the circular velocity of the dark matter haloes hosting the LAEs.

We present the clump populations detected in 18 lensed galaxies at redshifts 1 to 8.5 within the lensing cluster field SMACS0723. The recent JWST Early Release Observations of this poorly known region of the sky have revealed numerous point-like sources within and surrounding their host galaxies, undetected in the shallower HST images. We use JWST multiband photometry and the lensing model of this galaxy cluster to estimate the intrinsic sizes and magnitudes of the stellar clumps. We derive optical restframe effective radii from &lt;10 to hundreds pc and masses ranging from ∼105 to 109 M⊙, overlapping with massive star clusters in the local universe. Clump ages range from 1 Myr to 1 Gyr. We compare the crossing time to the age of the clumps and determine that between 45 and 60  per cent of the detected clumps are consistent with being gravitationally bound. On average, the dearth of Gyr old clumps suggests that the dissolution time scales are shorter than 1 Gyr. We see a significant increase in the luminosity (mass) surface density of the clumps with redshift. Clumps in reionisation era galaxies have stellar densities higher than star clusters in the local universe. We zoom in into single galaxies at redshift &lt;6 and find for two galaxies, the Sparkler and the Firework, that their star clusters/clumps show distinctive colour distributions and location surrounding their host galaxy that are compatible with being accredited or formed during merger events. The ages of some of the compact clusters are between 1 and 4 Gyr, e.g. globular cluster precursors formed around 9-12 Gyr ago. Our study, conducted on a small sample of galaxies, shows the potential of JWST observations for understanding the conditions under which star clusters form in rapidly evolving galaxies.

Exploiting the fundamentally achromatic nature of gravitational lensing, we present a lens model for the massive galaxy cluster SMACS J0723.3−7323 (SMACS J0723; z = 0.388) that significantly improves upon earlier work. Building on strong-lensing constraints identified in prior Hubble Space Telescope (HST) observations, the mass model utilizes 21 multiple-image systems, 17 of which were newly discovered in Early Release Observation data from the JWST. The resulting lens model maps the cluster mass distribution to an rms spatial precision of 0.″32, and is publicly available. Consistent with previous analyses, our study shows SMACS J0723.3 to be well described by a single large-scale component centered on the location of the brightest cluster galaxy. However, satisfying all lensing constraints provided by the JWST data, the model points to the need for the inclusion of an additional, diffuse component west of the cluster. A comparison of the galaxy, mass, and gas distributions in the core of SMACS J0723 based on HST, JWST, and Chandra data reveals a concentrated regular elliptical profile along with tell-tale signs of a recent merger, possibly proceeding almost along our line of sight. The exquisite sensitivity of JWST’s NIRCam reveals in spectacular fashion both the extended intracluster light distribution and numerous star-forming clumps in magnified background galaxies. The high-precision lens model derived here for SMACS J0723 demonstrates the unprecedented power of combining HST and JWST data for studies of structure formation and evolution in the distant universe.

Aims. We present the Lensed Lyman-Alpha MUSE Arcs Sample (LLAMAS) selected from MUSE and HST observations of 17 lensing clusters. The sample consists of 603 continuum-faint (-23