Search

Star clusters stand at the crossroads between galaxies and single stars. Resolving the formation of star clusters in cosmological simulations represents an ambitious and challenging goal, since modelling their internal properties requires very high resolution. This paper is the third of a series within the SImulating the Environment where Globular clusters Emerged (SIEGE) project, where we conduct zoom-in cosmological simulations with sub-parsec resolution that include the feedback of individual stars, aimed to model the formation of star clusters in high-redshift proto-galaxies. We investigate the role of three fundamental quantities in shaping the intrinsic properties of star clusters, i. e., i) pre-supernova stellar feedback (continuous or instantaneous ejection of mass and energy through stellar winds); ii) star formation efficiency, defined as the fraction of gas converted into stars per freefall time, for which we test 2 different values (epsi_ff=0.1 and 1), and iii) stellar initial mass function (IMF, standard vs top-heavy). All our simulations are run down to z=10.5, which is sufficient for investigating some structural properties of the emerging clumps and clusters. [Abridged] The prescription for a continuous, low-intensity feedback, along with the adoption of epsi_ff=1, produces star clusters with maximum stellar density values up to 10^4 M_sun pc^(-2), in good agreement with the surface density-size relation observed in local young star clusters (YSCs). Therefore, a realistic stellar wind description and a high star formation effiency are the key ingredients that allow us to achieve realistic star clusters characterised by properties comparable to those of local YSCs. In contrast, the other models produce too diffuse clusters, in particular the one with a top-heavy IMF., Comment: 20 pages, 14 pages

We aim at characterizing the X-ray-to-optical/near-infrared broad-band emission of luminous QSOs in the first Gyr of cosmic evolution to understand whether they exhibit differences compared to the lower-\textit{z} QSO population. Our goal is also to provide for these objects a reliable and uniform catalog of SED fitting derivable properties such as bolometric and monochromatic luminosities, Eddington ratios, dust extinction, strength of the hot dust emission. We characterize the X-ray/UV emission of each QSO using average SEDs from luminous Type 1 sources and calculate bolometric and monochromatic luminosities. Finally we construct a mean SED extending from the X-rays to the NIR bands. We find that the UV-optical emission of these QSOs can be modelled with templates of $z\sim$2 luminous QSOs. We observe that the bolometric luminosities derived adopting some bolometric corrections at 3000 \AA\ ($BC_{3000\text{\AA}}$) largely used in the literature are slightly overestimated by 0.13 dex as they also include reprocessed IR emission. We estimate a revised value, i.e. $BC_{3000\text{\AA}}=3.3 $ which can be used for deriving $L_\text{bol}$ in \textit{z} $\geq$ 6 QSOs. A sub-sample of 11 QSOs is provided with rest-frame NIR photometry, showing a broad range of hot dust emission strength, with two sources exhibiting low levels of emission. Despite potential observational biases arising from non-uniform photometric coverage and selection biases, we produce a X-ray-to-NIR mean SED for QSOs at \textit{z} $\gtrsim$ 6, revealing a good match with templates of lower-redshift, luminous QSOs up to the UV-optical range, with a slightly enhanced contribution from hot dust in the NIR.

We aim to gain new insights into the controversial dynamical status of MACS J0329-0211 (MACS0329), a massive cluster at z=0.4503, with a new analysis using a large sample of member galaxies as kinematic tracers. Our analysis is based on extensive spectroscopic data for more than 1700 galaxies obtained with the VIMOS and MUSE spectrographs at the Very Large Telescope (VLT), in combination with B and Rc Suprime-Cam photometry from the Subaru archive. According to our member selection procedure, we define a sample of 430 MACS0329 galaxies within 6 Mpc, corresponding to about 3 times the virial radius. We estimate the global velocity dispersion, sigmaV=841 km/s, and present the velocity dispersion profile. We estimate a mass M200=9.2E14 in units of solar masses, using 227 galaxies within R200=1.71 Mpc, for which sigmaV200=1018 km/s. The spatial distribution of the red galaxies traces a SE-NW elongated structure, without signs of a velocity gradient. This structure likely originates from the main phase of cluster assembly. The distribution of the blue galaxies is less concentrated, more rounded and shows signs of substructure, all characteristics indicating a recent infall of groups from the field. We detect two loose clumps of blue galaxies in the south and southwest at a distance about R200 from the cluster center. The strong spatial segregation among galaxy populations is not accompanied by kinematical difference. Thanks to our extensive catalog of spectroscopic redshift, we are able to study galaxy systems that are intervening along the line of sight. We can identify two foreground galaxy systems (GrG1 at z=0.31 and GrG2 at z=0.38) and one background system (GrG3 at z=0.47). We point out that the second brightest galaxy projected onto the MACS0329 core is in fact the dominant galaxy of the foreground group GrG2., Comment: Astronomy & Astrophysics accepted, 17 pages, 20 figures

One of the main open questions in the field of luminous ($L_{\rm bol}>10^{47}\,\rm erg\,s^{-1}$) quasars (QSOs) at $z \gtrsim 6$ is the rapid formation ($< 1\,$Gyr) of their supermassive black holes (SMBHs). For this work we analysed the relation between the X-ray properties and other properties describing the physics and growth of both the accretion disc and the SMBH in QSOs at the Epoch of Reionization (EoR). The sample consists of 21 $z>6$ QSOs, which includes 16 sources from the rapidly grown QSOs from the HYPERION sample and five other luminous QSOs with available high-quality archival X-ray data. We discovered a strong and statistically significant ($>3\sigma$) relation between the X-ray continuum photon index ($\Gamma$) and the $\rm C\,IV$ disc wind velocity ($v_{\rm C\,IV}$) in $z>6$ luminous QSOs, whereby the higher the $v_{\rm C\,IV}$, the steeper the $\Gamma$. This relation suggests a link between the disc-corona configuration and the kinematics of disc winds. Furthermore, we find evidence at $>2-3\sigma$ level that $\Gamma$ and $v_{\rm C\,IV}$ are correlated to the growth rate history of the SMBH. Although additional data are needed to confirm it, this result may suggest that, in luminous $z>6$ QSOs, the SMBH predominantly grows via fast accretion rather than via initial high seed BH mass., Comment: Accepted for publication in Astronomy & Astrophysics

We present JWST/NIRSpec PRISM follow-up of candidate galaxies at z=9-11 selected from deep JWST/NIRCam photometry in GLASS-JWST Early Release Science data. We spectroscopically confirm six sources with secure redshifts at z = 9.52-10.43, each showing multiple emission lines. An additional object is likely at z = 10.66, based on its Lya-break and a single emission feature, while one source is a lower redshift interloper. The sample includes the first JWST-detected candidate at z=10, GHZ1/GLASS-z10, which we confirm at z = 9.875, and the X-ray detected AGN GHZ9 confirmed at z = 10.145. Three objects in our sample, including GHZ9, have EW(CIII])>20A and occupy a region compatible with AGN emission in the EW(CIII]) vs CIV/CIII] diagram. The spectroscopic sample confirms a high abundance of galaxies at z > 9. We measure a number density of z=10 galaxies in the GLASS-JWST ERS field that is a factor of >3 higher than other JWST-based estimates at demagnified rest-frame magnitudes of -21 < Muv < -19. We find that the positions of these galaxies in redshift and angular space are not consistent with all of them being part of a unique progenitor of present-day galaxy clusters. The high density of objects in the GLASS region can be explained either by clustering on large scales or by a superposition of different forming structures of which we observe only the brightest members. By considering all the spectroscopic z=10 sources in the Abell-2744 field, we identify two potential galaxy proto-clusters centered around GHZ9 and JD1, with relative separations between their members of 1-2 pMpc. The potential AGN nature of three of the sources in our sample lends support to a scenario in which the high abundance of bright sources determined by JWST surveys at cosmic dawn may be affected by AGN contribution to their UV luminosity., Comment: Submitted to A&A

We present a detailed JWST/NIRSpec and NIRCam analysis of a gravitationally-lensed galaxy ($\rm \mu=17-21$) at redshift 6.14 magnified by the Hubble Frontier Field galaxy cluster MACS J0416. The target galaxy is overall a typical compact and UV-faint ($\rm M_{UV}=-17.8$) Lyman-$\alpha$ emitter; yet, the large magnification allows the detailed characterisation of structures on sub-galactic (down to few parsec) scales. Prominent optical $\rm H\alpha$, $\rm H\beta$ and [OIII]$\lambda\lambda4959,5007$ lines are spatially resolved with the high spectral resolution grating (G395H, R~2700), with large equivalent widths, EW($\rm H\beta$+[OIII])$\gtrsim1000$ \AA, and elevated ionising photon production efficiencies $\rm log(\xi_{ion}/erg^{-1}Hz)=25.2-25.7$. NIRCam deep imaging reveals the presence of compact rest-UV bright regions along with individual star clusters of sizes $\rm R_{eff}=3-8~pc$ and masses $\rm M\sim2\cdot10^5-5\cdot10^{6}~M_\odot$ These clusters are characterised by steep UV slopes, $\rm\beta_{UV}\lesssim-2.5$, in some cases associated with a dearth of line emission, indicating possible leaking of the ionising radiation, as also supported by a Lyman-$\rm \alpha$ emission peaking at $\rm \sim100~km~s^{-1}$ from the systemic redshift. While the entire system is characterised by low-metallicity, $\sim0.1~Z_\odot$, the NIRSpec-IFU map also reveals the presence of a low-luminosity, metal-poor region with $\rm Z\lesssim2\%~Z_\odot$, barely detected in NIRCam imaging; this region is displaced by $\rm >200~pc$ from one of the UV brightest structures of the system, and it would have been too faint to detect if not for the large magnification of the system., Comment: 19 pages (11 figures, 2 tables) + appendix (3 pages, 4 figures, 1 table). Submitted to A&A; comments are welcome

We report the serendipitous discovery of a faint (M_UV > -12.2), low-metallicity (Z ~ 0.02 Zsun), ionizing source (dubbed T2c) with a spectroscopic redshift of z=6.146. T2c is part of a larger structure amplified by the Hubble Frontier Field galaxy cluster MACSJ0416, and was observed with JWST/NIRSpec IFU. Stacking the short-wavelength NIRCam data reveals no stellar continuum detection down to a magnitude limit of m_UV ~ 31.0 (3 sigma). However, prominent Hb, [OIII]4959,5007, and Ha emissions are detected, with equivalent widths exceeding 200A, 800A, and 1300A (3 sigma), respectively. The corresponding intrinsic (magnification-corrected x23 +/- 3) ultraviolet and optical rest-frame magnitudes exceed 34.4 and 33.9 (corresponding to M_uv and M_opt fainter than -12.2 and -12.8, at lambda_rest ~ 2000A and ~5000A, respectively), suggesting a stellar mass lower than a few 10^4 Msun under an instantaneous burst scenario. The inferred ionizing photon production efficiency (xi_ion) is high, xi_ion >~ 26.08(25.86) 3(5)sigma, assuming no dust attenuation and no Lyman continuum leakage, indicating the presence of massive stars despite the low mass of the object. The very poor sampling of the initial mass function at such low mass star-forming complex suggests that the formation of very massive stars might be favored in very low metallicity environments. T2c is surrounded by Balmer and weak oxygen emission on a spatial scale of a few hundred parsecs after correcting for lensing effects. This system resembles an HII region potentially powered by currently undetected, extremely efficient, low-metallicity star complexes or clusters. We propose that massive O-type stars populate this low-mass and metallicity high-redshift satellites, likely caught in an early and short formation phase, contributing to the ionization of the surrounding medium., Comment: 9 pages, 5 figures, 1 table. Submitted to A&A. Comments are welcome

J1316+2614 at z=3.613 is the UV-brightest ($M_{UV}$ = -24.7) and strongest Lyman continuum (LyC, $f_{esc}^{LyC} \approx$ 90%) emitting star-forming galaxy known, showing also signatures of inflowing gas from its blue-dominated Ly$\alpha$ profile. Here, we present high-resolution imaging with the HST and VLT of the LyC, Ly$\alpha$, rest-UV, and optical emission of J1316+2614. Detailed analysis of the LyC and UV light distributions reveals compact yet resolved profiles, with LyC and UV morphologies showing identical half-light radii of $\simeq$ 220 pc. The continuum-subtracted Ly$\alpha$ emission reveals an extended filamentary structure of $\simeq$ 6.0 kpc oriented south-north with only weak/residual flux within the stellar core, suggesting a Ly$\alpha$ "hole". J1316+2614 presents remarkably high SFR and stellar mass surface densities of log($\Sigma_{SFR}$ [$M_{\odot}$/yr/kpc^2]) = 3.47$\pm$0.11 and log($\Sigma_{M}$ [$M_{\odot}$/pc^2]) = 4.20$\pm$0.06, respectively, which are among the highest observed in star-forming galaxies. Our findings indicate that J1316+2614 is a powerful, young, and compact starburst, leaking significant LyC photons due to the lack of gas and dust within the starburst. We explore the conditions for gas expulsion using a simple energetic balance and find that, given the strong binding force in J1316+2614, a high star formation efficiency ($\epsilon_{SF} \geq 0.7$) is necessary to remove the gas and explain its exposed nature. Our results thus suggest a close link between high $\epsilon_{SF}$ and high $f_{esc}^{LyC}$. This high efficiency can also naturally explain the remarkably high SFR, UV-luminosity, and efficient mass growth of J1316+2614, where at least 62% of its mass formed in the last 6 Myr. J1316+2614 may exemplify an intense, feedback-free starburst with a high $\epsilon_{SF}$, similar to those proposed for UV-bright galaxies at high redshifts. (ABRIDGED), Comment: 15 pages and 9 figures. Accepted for publication in A&A

The Extremely Large Telescopes (ELTs), thanks to their large apertures and cutting-edge Multi-Conjugate Adaptive Optics (MCAO) systems, promise to deliver sharper and deeper data even than the JWST. SHARP is a concept study for a near-IR (0.95-2.45 $\mu$m) spectrograph conceived to fully exploit the collecting area and the angular resolution of the upcoming generation of ELTs. In particular, SHARP is designed for the 2nd port of MORFEO@ELT. Composed of a Multi-Object Spectrograph, NEXUS, and a multi-Integral Field Unit, VESPER, MORFEO-SHARP will deliver high angular ($\sim$30 mas) and spectral (R$\simeq$300, 2000, 6000, 17000) resolution, outperforming NIRSpec@JWST (100 mas). SHARP will enable studies of the nearby Universe and the early Universe in unprecedented detail. NEXUS is fed by a configurable slit system deploying up to 30 slits with $\sim$2.4 arcsec length and adjustable width, over a field of about 1.2"$\times$1.2" (35 mas/pix). Each slit is fed by an inversion prism able to rotate by an arbitrary angle the field that can be seen by the slit. VESPER is composed of 12 probes of 1.7"$\times$1.5" each (spaxel 31 mas) probing a field 24"$\times$70". SHARP is conceived to exploit the ELTs apertures reaching the faintest flux and the sharpest angular resolution by joining the sensitivity of NEXUS and the high spatial sampling of VESPER to MORFEO capabilities. This article provides an overview of the scientific design drivers, their solutions, and the resulting optical design of the instrument achieving the required optical performance., Comment: Proceedings of the SPIE Astronomical Telescopes and Instrumentation 2024, Volume 13096, Paper No. 130965I, 11 pp, 11 figs

The study of gas-phase metallicity and its spatial distribution at high redshift is crucial to understand the processes that shaped the growth and evolution of galaxies in the early Universe. Here we study the spatially resolved metallicity in three systems at $z\sim6-8$, namely A2744-YD4, BDF-3299, and COSMOS24108, with JWST NIRSpec IFU low-resolution ($R\sim100$) spectroscopic observations. These are among the highest-$z$ sources in which metallicity gradients have been probed so far. Each of these systems hosts several spatial components in the process of merging within a few kpc, identified from the rest-frame UV and optical stellar continuum and ionised gas emission line maps. The sources have heterogeneous properties, with stellar masses log($M_*/M_\odot) \sim 7.6-9.3$, star formation rates (SFRs) $\sim1-15$ $M_\odot$ yr$^{-1}$, and gas-phase metallicities 12+log(O/H) $\sim 7.7-8.3$, which exhibit a large scatter within each system. Their properties are generally consistent with those of the highest-$z$ samples to date ($z\sim3-10$), though the sources in A2744-YD4 and COSMOS24108 are at the high end of the mass-metallicity relation (MZR) defined by the $z\sim3-10$ sources. Moreover, the targets in this work follow the predicted slope of the MZR at $z\sim 6-8$ from most cosmological simulations. The gas-phase metallicity gradients are consistent with being flat in the main sources of each system. Flat metallicity gradients are thought to arise from gas mixing processes on galaxy scales, such as mergers or galactic outflows and SN winds driven by intense stellar feedback, which wash out any gradient formed in the galaxy. The existence of flat gradients at $z\sim6-8$ sets also important constraints on cosmological simulations and chemical evolution models, whose predictions on the cosmic evolution of metallicity gradients differ significantly, but are mostly limited to $z<3$ so far., Comment: 15 pages (+5 appendix); 7 figures (+8 in appendix); in production on A&A

Observations at UV and optical wavelengths have revealed that galaxies at z~1-4 host star-forming regions, dubbed "clumps", which are believed to form due to the fragmentation of gravitationally unstable, gas-rich disks. However, the detection of the parent molecular clouds that give birth to such clumps is still possible only in a minority of galaxies, mostly at z~1. We investigated the [CII] and dust morphology of a z~3.4 lensed galaxy hosting four clumps detected in the UV continuum. We aimed to observe the [CII] emission of individual clumps that, unlike the UV, is not affected by dust extinction, to probe their nature and cold gas content. We conducted ALMA observations probing scales down to ~300 pc and detected three [CII] clumps. One (dubbed "NE") coincides with the brightest UV clump, while the other two ("SW" and "C") are not detected in the UV continuum. We do not detect the dust continuum. We converted the [CII] luminosity of individual clumps into molecular gas mass and found Mmol~10^8 Msun. By complementing it with the star formation rate (SFR) estimate from the UV continuum, we estimated the gas depletion time (tdep) of clumps and investigated their location in the Schmidt-Kennicutt plane. While the NE clump has a short tdep=0.16 Gyr, comparable with high-redshift starbursts, the SW and C clumps instead have longer tdep>0.65 Gyr and are likely probing the initial phases of star formation. The lack of dust continuum detection is consistent with the blue UV continuum slope estimated for this galaxy (beta~-2.5) and it indicates that dust inhomogeneities do not significantly affect the detection of UV clumps in this target. We pushed the observation of the cold gas content of individual clumps up to z~3.4 and showed that the [C II] line emission is a promising tracer of molecular clouds at high redshift, allowing the detection of clumps with a large range of depletion times., Comment: Accepted for publication in A&A

We present a novel approach to measuring the expansion rate and the geometry of the Universe, which combine time-delay cosmography in lens galaxy clusters with pure samples of 'cosmic chronometers' (CCs) by probing the member galaxies. The former makes use of the measured time delays between the multiple images of time-varying sources strongly lensed by galaxy clusters, while the latter exploits the most massive and passive cluster member galaxies to measure the differential time evolution of the Universe. We applied two different statistical techniques, adopting realistic errors on the measured quantities, to assess the accuracy and the gain in precision on the values of the cosmological parameters. We demonstrate that the proposed combined method allows for a robust and accurate measurement of the value of the Hubble constant. In addition, this provides valuable information on the other cosmological parameters thanks to the complementarity between the two different probes in breaking parameter degeneracies. Finally, we showcase the immediate observational feasibility of the proposed joint method by taking advantage of the existing high-quality spectro-photometric data for several lens galaxy clusters., Comment: 6 pages, 2 figures. Accepted for publication in A&A Letters

We present a simple and promising new method to measure the expansion rate and the geometry of the universe that combines observations related to the time delays between the multiple images of time-varying sources, strongly lensed by galaxy clusters, and Type Ia supernovae, exploding in galaxies belonging to the same lens clusters. By means of two different statistical techniques that adopt realistic errors on the relevant quantities, we quantify the accuracy of the inferred cosmological parameter values. We show that the estimate of the Hubble constant is robust and competitive, and depends only mildly on the chosen cosmological model. Remarkably, the two probes separately produce confidence regions on the cosmological parameter planes that are oriented in complementary ways, thus providing in combination valuable information on the values of the other cosmological parameters. We conclude by illustrating the immediate observational feasibility of the proposed joint method in a well-studied lens galaxy cluster, with a relatively small investment of telescope time for monitoring from a 2 to 3m class ground-based telescope., Comment: 6 pages, 2 figures, 1 table. Published in A&A Letters

We study the ultraviolet (UV) continuum $\beta$ slope of a sample of 166 clumps, individual star-forming regions observed in high redshift galaxies. They are hosted by 67 galaxies with redshift between 2 and 6.2, strongly lensed by the Hubble Frontier Fields cluster of galaxies MACS J0416.1-2403. The $\beta$ slope is sensitive to a variety of physical properties, such as the metallicity, the age of the stellar population, the dust attenuation throughout the galaxy, the stellar initial mass function (IMF), and the star-formation history (SFH). The aim of this study is to compare the $\beta$ values of individual clumps with those measured on the entire galaxy, to investigate possible physical differences between these regions and their hosts. We found a median value of $\beta \sim -2.4$, lower than that of integrated galaxies. This result confirms that clumps are sites of intense star formation, populated by young, massive stars, whose spectrum strongly emits in the UV. This is also consistent with the assumption that the dust extinction at the location of the clumps is lower than the average extinction of the galaxy, or that clumps have a different IMF or SFH. We made use of the correlations, discovered for high-redshift galaxies, of the $\beta$ value with those of redshift and UV magnitude, $M_{UV}$, finding that clumps follow the same relations, extended to much fainter magnitudes ($M_{UV}<-13$). We also find evidence of eight clumps with extremely blue ($\beta \lesssim -2.7$) slopes, which could be the signpost of low-metallicity stars and constrain the emissivity of ionizing photons at high redshift., Comment: 16 pages, 9 figures. Accepted for publication in MNRAS

VANDELS is a deep spectroscopic survey, performed with the VIMOS instrument at VLT, aimed at studying in detail the physical properties of high-redshift galaxies. VANDELS targeted about 2100 sources at 1

We use zoom-in, hydrodynamical, cosmological $N$-body simulations tracing the formation of the first stellar clumps from the SImulating the Environments where Globular clusters Emerged (SIEGE) project, to study key structural properties of dark matter haloes when the Universe was only $0.92$ Gyr old. The very high-resolution (maximum physical resolution 0.3 h$^{-1}$ pc at $z=6.14$, smallest dark-matter particle mass $164\,M_{\odot}$) allows us to reach the very low mass end of the stellar-to-halo mass relation ($M_{\rm vir}=10^{7.5-9.5}\,M_{\odot}$) to study the processes that mould dark matter haloes during the first stages of structure formation. We investigate the role of baryonic cooling and stellar feedback, modeled from individual stars, in shaping haloes, and of environmental effects as accretion of dark matter along cosmic filaments and mergers. We find that the onset of star formation (typically for $\log M_{\rm vir}/M_{\odot}\simeq7.6$) causes the inner cusp in the haloes density profile to flatten into a core with constant density and size proportionally to the halo virial mass. Even at these mass scales, we confirm that baryons make haloes that have formed stars rounder in the central regions than haloes that have not formed stars yet, with median minor-to-major $\langle q \rangle$ and intermediate-to-major $\langle s \rangle$ axes 0.66 and 0.84, respectively. Our morphological analysis shows that, at $z=6.14$, haloes are largely prolate in the outer parts, with the major axis aligned along filaments of the cosmic web or towards smaller sub-haloes, with the degree of elongation having no significant dependence on the halo mass., Comment: 20 pages, 11 figures. Accepted for publication in MNRAS

We present JWST/NIRSpec integral field spectroscopy (IFS) of a lensed Population III candidate stellar complex (dubbed Lensed And Pristine 1, LAP1), with a lensing-corrected stellar mass ~<10^4 Msun, absolute luminosity M_UV > -11.2 (m_UV > 35.6), confirmed at redshift 6.639 +/- 0.004. The system is strongly amplified (\mu >~ 100) by straddling a critical line of the Hubble Frontier Field galaxy cluster MACS J0416. Despite the stellar continuum is currently not detected in the Hubble and JWST/NIRCam and NIRISS imaging, arclet-like shapes of Lyman and Balmer lines, Lya, Hg, Hb and Ha are detected with NIRSpec IFS with signal-to-noise ratios SNR=5-13 and large equivalent widths (>300-2000A), along with a remarkably weak [OIII]4959-5007 at SNR ~ 4. LAP1 shows a large ionizing photon production efficiency, log(\xi_{ion}[erg~Hz^{-1}])>26. From the metallicity indexes R23 = ([OIII]4959-5007 + [OII]3727) / Hb ~< 0.74 and R3 = ([OIII]5007 / Hb) = 0.55 +/- 0.14, we derive an oxygen abundance 12+log(O/H) ~< 6.3. Intriguingly, the Ha emission is also measured in mirrored sub-components where no [OIII] is detected, providing even more stringent upper limits on the metallicity if in-situ star formation is ongoing in this region (12+log(O/H) < 6, or Z < 0.002 Zsun). The formal stellar mass limit of the sub-components would correspond to ~10^{3} Msun or M_UV fainter than -10. Alternatively, such a metal-free pure line emitting region could be the first case of a fluorescing HI gas region, induced by transverse escaping ionizing radiation from a nearby star-complex. The presence of large equivalent-width hydrogen lines and the deficiency of metal lines in such a small region, make LAP1 the most metal poor star-forming region currently known in the reionization era and a promising site that may host isolated, pristine stars., Comment: 15 pages, 10 figures, 1 table. Submitted to A&A. Comments welcome

We investigate the production efficiency of ionizing photons ($\xi_{ion}^*$) of 1174 galaxies with secure redshift at z=2-5 from the VANDELS survey to determine the relation between ionizing emission and physical properties of bright and massive sources. We constrain $\xi_{ion}^*$ and galaxy physical parameters by means of spectro-photometric fits performed with the BEAGLE code. The analysis exploits the multi-band photometry in the VANDELS fields, and the measurement of UV rest-frame emission lines (CIII]$\lambda 1909$, HeII$\lambda 1640$, OIII]$\lambda 1666$) from deep VIMOS spectra. We find no clear evolution of $\xi_{ion}^*$ with redshift within the probed range. The ionizing efficiency slightly increases at fainter $M_{UV}$, and bluer UV slopes, but these trends are less evident when restricting the analysis to a complete subsample at log(M$_{star}$/M$_{\odot}$)$>$9.5. We find a significant trend of increasing $\xi_{ion}^*$ with increasing EW(Ly$\alpha$), with an average log($\xi_{ion}^*$/Hz erg$^{-1}$)$>$25 at EW$>$50\AA, and a higher ionizing efficiency for high-EW CIII]$\lambda 1909$ and OIII]$\lambda 1666$ emitters. The most significant correlations are found with respect to stellar mass, specific star-formation rate (sSFR) and SFR surface density ($\Sigma_{SFR}$). The relation between $\xi_{ion}^*$ and sSFR shows a monotonic increase from log($\xi_{ion}^*$/Hz erg$^{-1}$) $\sim$24.5 at log(sSFR)$\sim$-9.5$yr^{-1}$ to $\sim$25.5 at log(sSFR)$\sim$-7.5$yr^{-1}$, a low scatter and little dependence on mass. The objects above the main-sequence of star-formation consistently have higher-than-average $\xi_{ion}^*$. A clear increase of $\xi_{ion}^*$ with $\Sigma_{SFR}$ is also found, with log($\xi_{ion}^*$/Hz erg$^{-1}$)$>$25 for objects at $\Sigma_{SFR}>$10 M$_{\odot}/yr/kpc^2$.(Abridged), Comment: 11 pages, 9 figures, accepted for publication in A&A

The existence of luminous quasars (QSO) at the Epoch of Reionization (EoR; i.e. z>6) powered by supermassive black holes (SMBH) with masses $\gtrsim10^9~M_\odot$ challenges models of early SMBH formation. To shed light on the nature of these sources we started a multiwavelength programme based on a sample of 18 HYPerluminous quasars at the Epoch of ReionizatION (HYPERION). These are the luminous QSOs whose SMBH must have had the fastest mass growth during the Universe first Gyr. In this paper we present the HYPERION sample and report on the first of the 3 years planned observations of the 2.4 Ms XMM-Newton Multi-Year Heritage program on which HYPERION is based. The goal of this program is to accurately characterize the X-ray nuclear properties of QSOs at the EoR. Through a joint X-ray spectral analysis of 10 sources, in the rest-frame $\sim2-50$ keV range, we report a steep average photon index ($\Gamma\sim2.4\pm0.1$). Absorption is not required. The average $\Gamma$ is inconsistent at $\geq4\sigma$ level with the canonical 1.8-2 value measured in QSO at z<6. This spectral slope is also much steeper than that reported in lower-z QSOs with similar luminosity or accretion rate, thus suggesting a genuine redshift evolution. Alternatively, we can interpret this result as the presence of an unusually low-energy cutoff $E_{cut}\sim20$ keV on a standard $\Gamma=1.9$ power-law. We also report on mild indications that HYPERION QSOs show higher soft X-ray emission at 2 keV compared to the UV one at 2500A than expected by lower-z luminous AGN. We speculate that a redshift-dependent coupling between the corona and accretion disc or intrinsically different coronal properties may account for the steep spectral slopes, especially in the presence of powerful winds. The reported slopes, if confirmed at lower luminosities, may have an important impact on future X-ray AGN studies in the early Universe., Comment: 21 pages (including appendix), 12 figures, 4 tables. Accepted for pubblication in A&A

Galaxy-scale strong lenses in galaxy clusters provide a unique tool to investigate their inner mass distribution and the sub-halo density profiles in the low-mass regime, which can be compared with the predictions from cosmological simulations. We search for galaxy-galaxy strong-lensing systems in HST multi-band imaging of galaxy cluster cores from the CLASH and HFF programs by exploring the classification capabilities of deep learning techniques. Convolutional neural networks are trained utilising highly-realistic simulations of galaxy-scale strong lenses injected into the HST cluster fields around cluster members. To this aim, we take advantage of extensive spectroscopic information on member galaxies in 16 clusters and the accurate knowledge of the deflection fields in half of these from high-precision strong lensing models. Using observationally-based distributions, we sample magnitudes, redshifts and sizes of the background galaxy population. By placing these sources within the secondary caustics associated with cluster galaxies, we build a sample of ~3000 galaxy-galaxy strong lenses which preserve the full complexity of real multi-colour data and produce a wide diversity of strong lensing configurations. We study two deep learning networks processing a large sample of image cutouts in three HST/ACS bands, and we quantify their classification performance using several standard metrics. We find that both networks achieve a very good trade-off between purity and completeness (85%-95%), as well as good stability with fluctuations within 2%-4%. We characterise the limited number of false negatives and false positives in terms of the physical properties of the background sources and cluster members. We also demonstrate the neural networks' high degree of generalisation by applying our method to HST observations of 12 clusters with previously known galaxy-scale lensing systems., Comment: 17 pages, 13 figures, to be published on A&A

The wavelength dependence of the Kormendy relation (KR) is well characterised at low redshift but poorly studied at intermediate redshifts. The KR provides information on the evolution of the population of early-type galaxies (ETGs). Therefore, by studying it, we may shed light on the assembly processes of these objects and their size evolution. As studies at different redshifts are generally conducted in different rest-frame wavebands, it is important to investigate whether the KR is dependent on wavelength. Knowledge of such a dependence is fundamental to correctly interpreting the conclusions we might draw from these studies. We analyse the KRs of the three Hubble Frontier Fields clusters, Abell S1063 (z = 0.348), MACSJ0416.1-2403 (z = 0.396), and MACS J1149.5+2223 (z = 0.542), as a function of wavelength. This is the first time the KR of ETGs has been explored consistently over such a large range of wavelengths at intermediate redshifts. We exploit very deep HST photometry, ranging from the observed B-band to the H-band, and MUSE integral field spectroscopy. We improve the structural parameter estimation we performed in a previous work by means of a newly developed Python package called morphofit. With its use on cluster ETGs, we find that the KR slopes increase smoothly with wavelength from the optical to the near-infrared (NIR) bands in all three clusters, with the intercepts becoming fainter at lower redshifts due to the passive ageing of the ETG stellar populations. The slope trend is consistent with previous findings at lower redshifts. The slope increase with wavelength implies that smaller ETGs are more centrally concentrated than larger ETGs in the NIR with respect to the optical regime. As different bands probe different stellar populations in galaxies, the slope increase also implies that smaller ETGs have stronger internal gradients with respect to larger ETGs., Comment: Accepted for publication in A&A as letter to the editor, 6 pages, 1 figure, 1 table

The redshift drift is a small, dynamic change in the redshift of objects following the Hubble flow. Its measurement provides a direct, real-time, model-independent mapping of the expansion rate of the Universe. It is fundamentally different from other cosmological probes: instead of mapping our (present-day) past light-cone, it directly compares different past light-cones. Being independent of any assumptions on gravity, geometry or clustering, it directly tests the pillars of the Lambda-CDM paradigm. Recent theoretical studies have uncovered unique synergies with other cosmological probes, including the characterization of the physical properties of dark energy. At the time of the original proposal by Sandage (1962) the expected change in the redshift of objects at cosmological distances appeared to be exceedingly small for reasonable observing times and beyond technological capabilities. In the last decades progress in the spectrographs (e.g. ESPRESSO), in the collecting area of telescopes and in the samples of cosmic beacons, enabled by new datasets and new machine-learning-based selections, have drastically changed the situation, bringing the Redshift Drift Grail within reach. As a consequence, this measurement is a flagship objective of the Extremely Large Telescope (ELT), specifically of its high-resolution spectrograph, ANDES., Comment: Proceedings of the Conference "Past, Present and Future of Astrophysical Spectroscopy", Trieste, 6-10 June 2022; accepted for publication; 8 pages, 4 figures

We report the identification of very massive stars (VMS; mass $> 100$\,\msun) that may be segregated in the center of the young massive star cluster at $z$=2.37 hosted in the lensed galaxy called {\tt Sunburst} galaxy. This result is based on two pieces of evidence: (1) VLT/MUSE spectra of several multiple images of the same star cluster show key spectral signatures of VMS, such as the \heii\ broad emission, \nivblue\ emission, and an \niv\ P-Cygni profile. In particular, \heii\ is broad ($\sim1610\pm300$ \kms), with an equivalent width of 3\AA,\ and asymmetric profile. These features require an extremely young ($\sim2.5$ Myr) stellar population component in which the masses of the stars exceed 100~\msun. When a Salpeter initial mass function and BPASS models for normal massive stars are assumed, the observed spectral features require $\sim$400 VMS. (2) The same star cluster is detected at a signal-to-noise ratio of~$\sim100$ in the Lyman continuum domain ($\lambda < 900$\AA). The Lyman continuum emission emerges from a region with a radius that is at least twice smaller than what is observed at 1700\AA~(independently of magnification) and is located in the center of the cluster. After delensing, the effective radii in absolute scales are R$_{\tt eff}[{\tt LyC}]\sim4.7 \pm 1.5$ pc and R$_{\tt eff}[1700]= 7.8 \pm 1.4$ pc. The Lyman continuum radiation is mainly produced by hot and massive stars, which implies that their spatial distribution (including that of VMS) is preferentially more confined in the central parts of the cluster. Approximately 400 VMS hosted by a cluster of $\sim 10^7$ \msun\ produce $\sim$15\% of the escaping Lyman continuum photons, and the remaining photons are produced by other massive early-type stars., Comment: 10 pages, 8 figures, Accepted to publication in A&A

In spite of the advent of extremely large telescopes in the UV/optical/NIR range, the current generation of 8-10m facilities is likely to remain competitive at ground-UV wavelengths for the foreseeable future. The Cassegrain U-Band Efficient Spectrograph (CUBES) has been designed to provide high-efficiency (>40%) observations in the near UV (305-400 nm requirement, 300-420 nm goal) at a spectral resolving power of R>20,000, although a lower-resolution, sky-limited mode of R ~ 7,000 is also planned. CUBES will offer new possibilities in many fields of astrophysics, providing access to key lines of stellar spectra: a tremendous diversity of iron-peak and heavy elements, lighter elements (in particular Beryllium) and light-element molecules (CO, CN, OH), as well as Balmer lines and the Balmer jump (particularly important for young stellar objects). The UV range is also critical in extragalactic studies: the circumgalactic medium of distant galaxies, the contribution of different types of sources to the cosmic UV background, the measurement of H2 and primordial Deuterium in a regime of relatively transparent intergalactic medium, and follow-up of explosive transients. The CUBES project completed a Phase A conceptual design in June 2021 and has now entered the Phase B dedicated to detailed design and construction. First science operations are planned for 2028. In this paper, we briefly describe the CUBES project development and goals, the main science cases, the instrument design and the project organization and management., Comment: Proceedings for the HACK100 conference, Trieste, June 2022. arXiv admin note: substantial text overlap with arXiv:2208.01672

GDS J033218.92-275302.7 (here GS-14) is a $z\sim5.5$ galaxy detected in [CII] as part of the ALPINE survey with unusual UV spectral features that have been interpreted as signatures of either a double stellar population or of an active galactic nucleus (AGN). We exploited the multi-wavelength coverage of GS-14 to investigate the properties and the origin of its emission. We performed UV-to-NIR SED-fitting, with single/double stellar population and/or AGN component. We analyzed the VIMOS spectrum, which shows highly-ionized emission lines (Ovi, Nv, and Niv). The line properties have been compared with those observed in galaxies and AGN, and with the predictions from radiation transfer models for star-forming galaxies, AGN, and shocks. The SED-fitting provides a total stellar mass of $M_*=(4 \pm 1) \times 10^{10} M_\odot$, an age of the main stellar population of $\sim670 Myr$ and a recent short (8 Myr) burst of star formation (SF) of $\sim 90 M_\odot yr^{-1}$. The Nv line has a characteristic P-Cygni profile, which suggests a $\sim 3 Myr$ old population of stars with a mass of $\sim 5 \times10^{7} M_\odot$. The Nv profile also shows evidence for an additional component of nebular emission. The comparison of the line ratios with theoretical models allows us to associate the emission with SF or AGN, but the strong radiation field required to ionize the Ovi is more commonly related to AGN activity. We found evidence for an old and already evolved stellar population at $z\sim 5.5$ and show that the galaxy is experiencing a second short burst of SF. In addition, GS-14 carries signatures of obscured AGN activity. The AGN could be responsible for the short depletion time of this galaxy, thus making GS-14 one of the two ALPINE sources with hints of an active nucleus and an interesting target for future follow-ups., Comment: Accepted in A&A, 10 pages, 4 figures

Star cluster formation in the early universe and their contribution to reionization remains to date largely unconstrained. Here we present JWST/NIRCam imaging of the most highly magnified galaxy known at z ~ 6, the Sunrise arc. We identify six young massive star clusters (YMCs) with measured radii spanning ~ 20 pc down to ~ 1 pc (corrected for lensing magnification), estimated stellar masses of ~ $10^{(6-7)}$ Msun, and with ages 1-30 Myr based on SED fitting to photometry measured in 8 filters extending to rest-frame 7000A. The resulting stellar mass surface densities are higher than 1000 Msun pc$^{-2}$ (up to a few $10^5$ Msun pc$^{-2}$) and their inferred dynamical ages qualify the majority of these systems as gravitationally-bound stellar clusters. The star cluster ages map the progression of star formation along the arc, with to evolved systems (>~ 10 Myr old) followed by very young clusters. The youngest stellar clusters (< 5 Myr) show evidence of prominent Hbeta + [OIII]4959,5007 emission, based on photometry, with equivalent widths larger than 1000 A rest-frame, and are hosted in a 200 pc sized star-forming complex. Such a region dominates the ionizing photon production, with a high efficiency log($\xi_{ion}$ [Hz erg$^{-1}$]) ~ 25.7. A significant fraction of the recently formed stellar mass of the galaxy (> 10-30 %) occurred in these YMCs. We speculate that such sources of ionizing radiation boost the ionizing photon production efficiency which eventually carve ionized channels that might favor the escape of Lyman continuum radiation. The survival of some of the clusters would make them the progenitors of massive and relatively metal-poor globular clusters in the local Universe., Comment: 18 pages, 8 figures, 3 tables. ApJ, Accepted

The origin of the Lyman-${\alpha}$ (Ly${\alpha}$) emission in galaxies is a long-standing issue: despite several processes known to originate this line (e.g. AGN, star formation, cold accretion, shock heating), it is difficult to discriminate among these phenomena based on observations. Recent studies have suggested that the comparison of the ultraviolet (UV) and optical properties of these sources could solve the riddle. For this reason, we investigate the rest-frame UV and optical properties of A2895b, a strongly lensed Ly${\alpha}$-emitter at redshift z ~ 3.7. From this study, we find that our target is a compact (r ~ 1.2 pkpc) star-forming (star formation rate ~ 11 M$_{\odot}$/yr) galaxy having a young stellar population. Interestingly, we measure a high ratio of the H${\beta}$ and the UV continuum monochromatic luminosities (L(H${\beta}$)/L(UV) ~ 100). Based on tracks of theoretical stellar models (Starburst99, BPASS), we can only partially explain this result by assuming a recent (< 10 Myr), bursty episode of star-formation and considering models characterised by binary stars, a top-heavy initial-mass function (IMF) and sub-solar metallicities (Z < 0.01 Z$_{\odot}$). These assumptions also explain the observed low (C/O) abundance of our target (~ 0.23(C/O)$_{\odot}$). By comparing the UV and optical datasets, we find that the Ly${\alpha}$ and UV continuum are more extended (x2) than the Balmer lines, and that the peak of the Ly${\alpha}$ is offset (~ 0.6 pkpc). The multi-wavelength results of our analysis suggest that the observed Ly${\alpha}$ emission originates from a recent star-formation burst, likely taking place in an off-centre clump., Comment: 18 pages (including Appendix), 5 figures. Paper accepted for publication on MNRAS. Accepted 2022 November 2. Received 2022 October 22; in original form 2022 August 3

We present the discovery and analysis of J1316+2614 at z=3.6130, a UV-bright star-forming galaxy ($M_{\rm UV} \simeq -24.7$) with large escape of Lyman continuum (LyC) radiation. J1316+2614 is a young ($\simeq 10$ Myr) star-forming galaxy with $SFR \simeq 500 M_{\odot}$ yr$^{-1}$ and a starburst mass of log($M_{\star}/M_{\odot}) \simeq 9.7$. It shows a very steep UV continuum, $\beta_{\rm UV} \simeq -2.59 \pm 0.05$, consistent with residual dust obscuration, $E(B-V)\simeq 0$. LyC emission is detected with high significance ($\simeq 17 \sigma$) down to $830$\r{A}, for which a very high relative (absolute) LyC escape fraction $f_{\rm esc} \rm (LyC) \simeq 0.92$ ($\simeq 0.87$) is inferred. The contribution of a foreground or AGN contamination to the LyC signal is discussed but is unlikely. J1316$+$2614 is the most powerful ionizing source known among the star-forming galaxy population, both in terms of production ($Q_{\rm H} \approx 10^{56}$ s$^{-1}$) and escape of ionizing photons ($f_{\rm esc} \rm (LyC) \approx 0.9$). Nebular emission in Ly$\alpha$, H$\beta$, and other rest-frame optical lines are detected, but these are weak ($EW_{0} \rm [H\beta] \simeq 35$\r{A}), with their strengths reduced roughly by $\simeq 90\%$. J1316+2614 is the first case known where the effect of large escape of ionizing photons on the strength of nebular lines and continuum emission is clearly observed. Gas inflows are detected in J1316+2614 from the blue-dominated peak Ly$\alpha$ emission (with a blue-to-red peak line ratio $I_{\rm blue}/I_{\rm red} \simeq 3.7$) and redshifted ISM absorption ($\simeq 100$ km s$^{-1}$). Our results suggest that J1316+2614 is undergoing a gas compaction event, possibly representing a short-lived phase in the evolution of massive and compact galaxies, where strong gas inflows have triggered an extreme star formation episode and nearly $100\%$ LyC photons are escaping., Comment: Accepted for publication in MNRAS. 18 pages, 10 figures, 5 tables

We present a detailed strong lensing analysis of the massive and distant ($z=0.870$) galaxy cluster ACT-CL J0102$-$4915 (ACT0102, also known as El Gordo), taking advantage of new spectroscopic data from the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope, and archival imaging from the Hubble Space Telescope. Thanks to the MUSE data, we measure secure redshifts for 374 single objects, including 23 multiply lensed galaxies, and 167 cluster members of ACT0102. The observed positions of 56 multiple images, along with their new spectroscopic redshift measurements, are used as constraints for our strong lensing model. Remarkably, some multiple images are detected out to a large projected distance of $\approx 1$ Mpc from the brightest cluster galaxy, allowing us to estimate a projected total mass value of $1.84_{-0.04}^{+0.03} \times 10^{15}\, \rm M_{\odot}$ within that radius. We find that we need two extended cluster mass components, the mass contributions from the cluster members and the additional lensing effect of a foreground ($z=0.633$) group of galaxies, to predict the positions of all multiple images with a root mean square offset of $0.75"$. The main cluster-scale mass component is centered very closely to the brightest cluster galaxy and the other extended mass component is located in the north-west region of the cluster. These two mass components have very similar values of mass projected within 300 kpc from their centers, namely $2.29_{-0.10}^{+0.09}\times10^{14}\,\rm M_{\odot}$ and $2.10_{-0.09}^{+0.08}\times10^{14}\,\rm M_{\odot}$, in agreement with the major merging scenario of ACT0102. We make publicly available the lens model, including the magnification maps and posterior distributions of the model parameter values, as well as the full spectroscopic catalogue containing all redshift measurements obtained with MUSE., Comment: 12 pages, 9 figures, 4 tables + appendix. Minor changes, A&A in press. MUSE redshift catalogue and lens model are publicly available

The combination of multi-band imaging from HST with MUSE integral field spectroscopy, obtained at the VLT, has recently driven remarkable progress in strong lensing (SL) modeling of galaxy clusters. From a few tens of multiple images with photometric redshifts per cluster, a new generation of high-precision SL models have recently been developed, by exploiting in some cases over a hundred of spectroscopically confirmed multiple images and cluster member galaxies. A further step forward is expected with JWST observations of SL clusters (from hundreds to possibly a thousand of multiple images). In this context, we present a new, state-of-the-art SL model of the galaxy cluster MACS J0416.1-2403, utilizing 237 spectroscopically confirmed multiple images, which is the largest sample of secure multiply lensed sources utilized to date. This model incorporates stellar kinematics information of 64 cluster galaxies and the hot-gas mass distribution of the cluster determined from Chandra X-ray observations. The observed positions of the many multiple images are reproduced with a remarkable accuracy of 0.43 arcsec. To further assess the reliability of this lens model and to highlight the improvement over previously published models, we show the extended surface brightness reconstruction of several lensed galaxies through a newly developed forward modeling software. The comparison with other SL models of the same cluster demonstrates that this new model is better suited to accurately reproduce the positions, shapes and fluxes of the observed multiple images. Besides a robust characterization of the total mass distribution of the cluster, our model can provide accurate and precise magnification maps that are key to studying the intrinsic physical properties of faint, high-redshift lensed sources. The model is made publicly available through our newly developed Strong Lensing Online Tool (or SLOT)., Comment: Accepted for publication in Astronomy & Astrophysics (16 pages, 13 figures)

SDSS J2222+2745, at z = 0.489, is one of the few currently known lens clusters with multiple images of a background (z = 2.801) quasar with measured time delays. We combine imaging from the Hubble Space Telescope (HST) with recent Multi Unit Spectroscopic Explorer (MUSE) spectroscopic data to securely identify 34 cluster members and 12 multiple images from 3 background sources. We measure the stellar velocity dispersions of 13 cluster galaxies, enabling an independent estimate of the contribution of the sub-halo mass component to the lens total mass. The projected total mass distribution of the lens cluster is best modelled with a single large-scale mass component, a galaxy-scale component, anchored by the MUSE kinematic information, and an external shear. The best-fit strong lensing model yields a root mean square separation between the model-predicted and observed positions of the multiple images of 0".29. When analysing the impact of systematic uncertainties, stemming from modelling assumptions and used observables, we find that the projected total mass profile, relative weight of the sub-halo mass component, and critical lines are consistent, within the statistical uncertainties. The predicted magnification and time delay values are, instead, more sensitive to the local details of the lens total mass distribution, and vary significantly among lens models that are similarly good at reproducing the observed multiple image positions. Due to its complex morphology, the low number of point-like multiple images, and current model degeneracies, it becomes clear that additional information (from the observed surface brightness distribution of lensed sources and the measured time delays) needs to be included in the modelling of SDSS J2222+2745 for accurate and precise cosmological measurements. The full MUSE secure spectroscopic catalogue presented in this work is made publicly available., Comment: 17 pages, 9 Figures, 5 Tables. Accepted for publication in Astronomy & Astrophysics. The full MUSE spectroscopic catalogue presented in this work is made publicly available at www.fe.infn.it/astro/lensing

In the era of Extremely Large Telescopes, the current generation of 8-10m facilities are likely to remain competitive at ground-UV wavelengths for the foreseeable future. The Cassegrain U-Band Efficient Spectrograph (CUBES) has been designed to provide high-efficiency (>40%) observations in the near UV (305-400 nm requirement, 300-420 nm goal) at a spectral resolving power of R>20,000 (with a lower-resolution, sky-limited mode of R ~ 7,000). With the design focusing on maximizing the instrument throughput (ensuring a Signal to Noise Ratio (SNR) ~20 per high-resolution element at 313 nm for U ~18.5 mag objects in 1h of observations), it will offer new possibilities in many fields of astrophysics, providing access to key lines of stellar spectra: a tremendous diversity of iron-peak and heavy elements, lighter elements (in particular Beryllium) and light-element molecules (CO, CN, OH), as well as Balmer lines and the Balmer jump (particularly important for young stellar objects). The UV range is also critical in extragalactic studies: the circumgalactic medium of distant galaxies, the contribution of different types of sources to the cosmic UV background, the measurement of H2 and primordial Deuterium in a regime of relatively transparent intergalactic medium, and follow-up of explosive transients. The CUBES project completed a Phase A conceptual design in June 2021 and has now entered the detailed design and construction phase. First science operations are planned for 2028., Comment: SPIE proceedings, SPIE Astronomical Telescopes + Instrumentation 2022, Montr\'eal, Canada; 20 pages, 13 figures, 2 tables

We investigate the blue and optical rest-frame sizes (lambda~2300A-4000A) of three compact star-forming regions in a galaxy at z=4 strongly lensed (x30, x45, x100) by the Hubble Frontier Field galaxy cluster A2744 using GLASS-ERS JWST/NIRISS imaging at 1.15um, 1.50mu and 2.0mu with PSF < 0.1". In particular, the Balmer break is probed in detail for all multiply-imaged sources of the system. With ages of a few tens of Myr, stellar masses in the range (0.7-4.0) x 10^6 Msun and optical/ultraviolet effective radii spanning the interval 3 < R_eff < 20 pc, such objects are currently the highest redshift (spectroscopically-confirmed) gravitationally-bound young massive star clusters (YMCs), with stellar mass surface densities resembling those of local globular clusters. Optical (4000A, JWST-based) and ultraviolet (1600A, HST-based) sizes are fully compatible. The contribution to the ultraviolet underlying continuum emission (1600A) is ~30%, which decreases by a factor of two in the optical for two of the YMCs (~4000A rest-frame), reflecting the young ages (<30 Myr) inferred from the SED fitting and supported by the presence of high-ionization lines secured with VLT/MUSE. Such bursty forming regions enhance the sSFR of the galaxy, which is ~10 Gyr^-1. This galaxy would be among the extreme analogs observed in the local Universe having high star formation rate surface density and high occurrence of massive stellar clusters in formation., Comment: 10 pages, 4 figures, 1 table. Submitted to ApJL

We present the serendipitous discovery of a late T-type brown dwarf candidate in JWST NIRCam observations of the Early Release Science Abell 2744 parallel field. The discovery was enabled by the sensitivity of JWST at 4~$\mu$m wavelengths and the panchromatic 0.9--4.5~$\mu$m coverage of the spectral energy distribution. The unresolved point source has magnitudes F115W = 27.95$\pm$0.15 and F444W = 25.84$\pm$0.01 (AB), and its F115W$-$F444W and F356W$-$F444W colors match those expected for other, known T dwarfs. We can exclude it as a reddened background star, high redshift quasar, or very high redshift galaxy. Comparison with stellar atmospheric models indicates a temperature of $T_{eff}$ $\approx$ 600~K and surface gravity $\log{g}$ $\approx$ 5, implying a mass of 0.03~M$_{\odot}$ and age of 5~Gyr. We estimate the distance of this candidate to be 570--720~pc in a direction perpendicular to the Galactic plane, making it a likely thick disk or halo brown dwarf. These observations underscore the power of JWST to probe the very low-mass end of the substellar mass function in the Galactic thick disk and halo., Comment: ApJL accepted

Star-forming galaxies can exhibit strong morphological differences between the rest-frame far-UV and optical, reflecting inhomogeneities in star-formation and dust attenuation. We exploit deep, high resolution NIRCAM 7-band observations to take a first look at the morphology of galaxies in the epoch of reionization ($z>7$), and its variation in the rest-frame wavelength range between Lyman $\alpha$ and 6000-4000\AA, at $z=7-12$. We find no dramatic variations in morphology with wavelength -- of the kind that would have overturned anything we have learned from the Hubble Space Telescope. No significant trends between morphology and wavelengths are detected using standard quantitative morphology statistics. We detect signatures of mergers/interactions in 4/19 galaxies. Our results are consistent with a scenario in which Lyman Break galaxies -- observed when the universe is only 400-800 Myrs old - are growing via a combination of rapid galaxy-scale star formation supplemented by accretion of star forming clumps and interactions., Comment: 15 pages, 4 figures, ApJL in press

Wide Field Slitless Spectroscopy (WFSS) provides a powerful tool for detecting strong line emission in star forming galaxies (SFGs) without the need for target pre-selection. As part of the GLASS-JWST-ERS program, we leverage the near-infrared wavelength capabilities of NIRISS ($1-2.2\mu$m) to observe rest-optical emission lines out to $\rm{z}\sim 3.4$, to a depth and with a spatial resolution higher than ever before (H$\alpha$ to z<2.4; [OIII]+H$\beta$ to z<3.4). In this letter we constrain the rest-frame [OIII]$\lambda5007$ equivalent width (EW) distribution for a sample of 76 $1<\rm{z}<3.4$ SFGs in the Abell 2744 Hubble Frontier Field and determine an abundance fraction of extreme emission line galaxies with EW$>750$A in our sample to be $12\%$. We determine a strong correlation between the measured H$\beta$ and [OIII]$\lambda5007$ EWs, supporting that the high [OIII]$\lambda5007$ EW objects require massive stars in young stellar populations to generate the high energy photons needed to doubly ionise oxygen. We extracted spectra for objects up to 2 mag fainter in the near-infrared than previous WFSS studies with the Hubble Space Telescope. Thus, this work clearly highlights the potential of JWST/NIRISS to provide high quality WFSS datasets in crowded cluster environments., Comment: 9 pages, 3 figures, submitted to APJL

We present the first rest-frame optical size-luminosity relation of galaxies at $z>7$, using the NIRCam imaging data obtained by the GLASS James Webb Space Telescope Early Release Science (GLASS-JWST-ERS) program, providing the deepest extragalactic data of the ERS campaign. Our sample consist of 19 photometrically selected bright galaxies with $m_\text{F444W}\leq27.8$ at $77$. Given the limited sample size and magnitude range, we first fix the slope to that observed for larger samples in rest-frame UV using HST samples. The median size $r_0$ at the reference luminosity $M=-21$ decreases slightly from rest-frame optical ($600\pm80$ pc) to UV ($450\pm130$ pc). We then re-fit the size-luminosity relation allowing the slope to vary. The slope is consistent with $\beta\sim0.2$ for all bands except F150W, where we find a marginally steeper slope of $\beta=0.53\pm0.15$. The steep UV slope is mainly driven by the smallest and faintest galaxies. If confirmed by larger samples, it implies that the UV size-luminosity relation breaks toward the faint end as suggested by lensing studies., Comment: Accepted by ApJL, 9 pages, 3 figures, 2 tables

We present a new strong lensing (SL) model of the Hubble Frontier Fields galaxy cluster Abell 2744, at z=0.3072, by exploiting archival Hubble Space Telescope (HST) multi-band imaging and Multi Unit Spectroscopic Explorer (MUSE) follow-up spectroscopy. The lens model considers 90 spectroscopically confirmed multiple images (from 30 background sources), which represents the largest secure sample for this cluster field prior to the recently acquired James Webb Space Telescope observations. The inclusion of the sub-structures within several extended sources as model constraints allows us to accurately characterize the inner total mass distribution of the cluster and the position of the cluster critical lines. We include the lensing contribution of 225 cluster members, 202 of which are spectroscopically confirmed. We also measure the internal velocity dispersion of 85 cluster galaxies to independently estimate the role of the subhalo mass component in the lens model. We investigate the effect of the cluster environment on the total mass reconstruction of the cluster core with two different mass parameterizations. We consider the mass contribution from three external clumps, either based on previous weak-lensing studies, or extended HST imaging of luminous members around the cluster core. In the latter case, the observed positions of the multiple images are better reproduced, with a remarkable accuracy of 0.37", a factor of $\sim2$ smaller than previous lens models. We develop and make publicly available a Strong Lensing Online Tool (SLOT) to exploit the predictive power and the full statistical information of this and future models, through a simple graphical interface. We plan to apply our high-precision SL model to the first analysis of the GLASS-JWST-ERS program, specifically to measure the intrinsic physical properties of high-$z$ galaxies from robust magnification maps., Comment: 14 pages, 9 figures, 4 tables

We present our lens mass model of SMACS J0723, the first strong gravitational lens observed by the James Webb Space Telescope (JWST). We use data from the Hubble Space Telescope and Multi Unit Spectroscopic Explorer (MUSE) to build our 'pre-JWST' lens model, and refine it with newly available JWST near-infrared imaging in our JWST model. To reproduce the positions of all multiple lensed images with good accuracy, the adopted mass parameterization consists of one cluster-scale component, accounting mainly for the dark matter distribution, the galaxy cluster members and an external shear component. The pre-JWST model has, as constraints, 19 multiple images from six background sources, of which four have secure spectroscopic redshift measurements from this work. The JWST model has more than twice the number of constraints, 30 additional multiple images from another eleven lensed sources. Both models can reproduce very well the multiple image positions with a $\delta_{rms}$ of $0.39''$ and $0.51''$, for the pre-JWST and JWST models, respectively. The total mass estimates within a radius of 128~kpc (~ the Einstein radius) are $7.9_{-0.2}^{+0.3}\times 10^{13}\rm M_{\odot}$ and $8.7_{-0.2}^{+0.2}\times 10^{13}\rm M_{\odot}$, for the pre-JWST and JWST models, respectively. We predict with our mass models the redshifts of the newly detected JWST sources, which are crucial information for systems without spectroscopic measurements for further studies and follow-up observations. Interestingly, one family detected with JWST is found to be at a very high redshift, $z>7.5$ (68% confidence level) and with one image having lensing magnification of $|\mu|=9.5_{-0.8}^{+0.9}$, making it an interesting case for future studies. The lens models, including magnification maps and redshifts estimated from the model are made publicly available, along with the full spectroscopic redshift catalogue from MUSE., Comment: 5 pages, 2 figures, 1 table, plus appendix. Published in A&A Letters. Lens model and redshift catalogue available at: http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/666/L9

We introduce a new set of zoom-in cosmological simulations with sub-pc resolution, intended to model extremely faint, highly magnified star-forming stellar clumps, detected at z=6.14 thanks to gravitational lensing. The simulations include feedback from individual massive stars (in both the pre-supernova and supernova phases), generated via stochastic, direct sampling of the stellar initial mass function. We adopt a modified 'delayed cooling' feedback scheme, specifically created to prevent artificial radiative loss of the energy injected by individual stars in very dense gas (n~10^3-10^5 cm^{-3}). The sites where star formation ignites are characterised by maximum densities of the order of 10^5 cm^{-3} and gravitational pressures P/k>10^7 K/cm^3, corresponding to the values of the local, turbulent regions where the densest stellar aggregates form. The total stellar mass at z=6.14 is 3.4x10^7 M_sun, in satisfactory agreement with the observed stellar mass of the observed systems. The most massive clumps have masses of ~10^6 M_sun and half-mass sizes of ~100 pc. These sizes are larger than the observed ones, including also other samples of lensed high-redshift clumps, and imply an average density one order of magnitude lower than the observed one. In the size-mass plane, our clumps populate a sequence that is intermediate between the ones of observed high-redshift clumps and local dSph galaxies., Comment: MNRAS, accepted for publication. 24 pages, 17 figures

The GLASS James Webb Space Telescope Early Release Science (hereafter GLASS-JWST-ERS) Program will obtain and make publicly available the deepest extragalactic data of the ERS campaign. It is primarily designed to address two key science questions, namely, "what sources ionized the universe and when?" and "how do baryons cycle through galaxies?", while also enabling a broad variety of first look scientific investigations. In primary mode, it will obtain NIRISS and NIRSpec spectroscopy of galaxies lensed by the foreground Hubble Frontier Field cluster, Abell 2744. In parallel, it will use NIRCam to observe two fields that are offset from the cluster center, where lensing magnification is negligible, and which can thus be effectively considered blank fields. In order to prepare the community for access to this unprecedented data, we describe the scientific rationale, the survey design (including target selection and observational setups), and present pre-commissioning estimates of the expected sensitivity. In addition, we describe the planned public releases of high-level data products, for use by the wider astronomical community., Comment: 23 pages, 10 figures; ApJ in press

We investigate the ionizing properties of the pair of bright Ly$\alpha$ emitting galaxies BDF521 and BDF2195 at z=7.012 in order to constrain their contribution to the formation of the BDF "reionized bubble". We obtain constraints on UV emission lines (CIV$\lambda 1548$ doublet, HeII$\lambda 1640$, OIII]$\lambda 1660$ doublet, and CIII]$\lambda 1909$ doublet) from deep VLT-XSHOOTER observations and compare them to those available for other high-redshift objects, and to models with mixed stellar and AGN emission. We use this spectroscopic information together with the photometry available in the field to constrain the physical properties of the two objects using the spectro-photometric fitting code BEAGLE. We do not detect any significant emission at the expected position of UV lines, with 3$\sigma$ upper limits of EW$\lesssim$2-7 AA rest-frame. We find that the two objects have lower CIII] emission than expected on the basis of the correlation between the Ly$\alpha$ and CIII] EWs. The EW limits on CIV and HeII emission exclude pure AGN templates at $\sim2-3\sigma$ significance, and only models with a $\lesssim$40% AGN contribution are compatible with the observations. The two objects are found to be relatively young ($\sim$20-30 Myrs) and metal-poor ($\lesssim 0.3 Z_{\odot}$) with stellar masses of a few $10^9M_{\odot}$. Their production rate of hydrogen ionizing photons per intrinsic UV luminosity is log($\xi_{ion}^*$/Hz erg$^{-1}$)=25.02-25.26, consistent with values typically found in high-redshift galaxies, but more than twice lower than values measured in $z>$7 galaxies with strong CIII] and/or optical line emission ($\simeq$25.6-25.7). The two BDF emitters have no evidence of higher than average ionizing capabilities and are not capable of reionizing their surroundings by their own means under realistic assumptions on the escape fraction of ionizing photons. (abridged), Comment: 9 pages, 4 figures, accepted for publication in A&A

We study the physical properties (size, stellar mass, luminosity, star formation rate) and scaling relations for a sample of 166 star-forming clumps with redshift $z \sim 2-6.2$. They are magnified by the Hubble Frontier Field galaxy cluster MACS~J0416 and have robust lensing magnification ($2\lesssim \mu \lesssim 82$) computed by using our high-precision lens model, based on 182 multiple images. Our sample extends by $\sim 3$ times the number of spectroscopically-confirmed lensed clumps at $z \gtrsim 2$. We identify clumps in ultraviolet continuum images and find that, whenever the effective spatial resolution (enhanced by gravitational lensing) increases, they fragment into smaller entities, likely reflecting the hierarchically-organized nature of star formation. Kpc-scale clumps, most commonly observed in field, are not found in our sample. The physical properties of our sample extend the parameter space typically probed by $z \gtrsim 1$ field observations and simulations, by populating the low mass (M$_\star \lesssim 10^7$ M$_\odot$), low star formation rate (SFR $\lesssim 0.5$ M$_\odot$ yr$^{-1}$), and small size (R$_\mathrm{eff} \lesssim 100$ pc) regime. The new domain probed by our study approaches the regime of compact stellar complexes and star clusters. In the mass-size plane, our sample spans the region between galaxies and globular clusters, with a few clumps in the region populated by young star clusters and globular-clusters. For the bulk of our sample, we measure star-formation rates which are higher than those observed locally in compact stellar systems, indicating different conditions for star formation at high redshift than in the local Universe., Comment: 25 pages, 20 figures, 3 table. Accepted to publication in MNRAS

We present a study designed to measure the average LyC escape fraction ($\langle f_{\rm esc}\rangle$) of star-forming galaxies at z=3.5. We assemble a sample of 148 galaxies from the VANDELS survey at $3.35\leq z_{\rm spec}\leq3.95$, selected to minimize line-of-sight contamination of their photometry. For this sample, we use ultra-deep, ground-based, $U-$band imaging and HST $V-$band imaging to robustly measure the distribution of $\mathcal{R_{\rm obs}}$ $=(L_{\rm LyC}/L_{\rm UV})_{\rm obs}$. We then model the distribution as a function of $\langle f_{\rm esc}\rangle$, carefully accounting for attenuation by dust, and the IGM (and CGM). A maximum likelihood fit to the $\mathcal{R_{\rm obs}}$ distribution returns a best-fitting value of $\langle f_{\rm esc}\rangle =0.07\pm0.02$, a result confirmed using an alternative Bayesian inference technique (both exclude $\langle f_{\rm esc}\rangle=0.0$ at $> 3\sigma$). By splitting our sample in two, we find evidence that $\langle f_{\rm esc}\rangle$ is positively correlated with Ly$\alpha$ equivalent width, with high and low sub-samples returning best fits of $\langle f_{\rm esc}\rangle=0.12^{+0.06}_{-0.04}$ and $\langle f_{\rm esc} \rangle=0.02^{+0.02}_{-0.01}$, respectively. In contrast, we find evidence that $\langle f_{\rm esc}\rangle$ is anti-correlated with intrinsic UV luminosity and UV dust attenuation; with low UV luminosity and dust attenuation sub-samples returning best fits in the range $0.10 \leq \langle f_{\rm esc}\rangle \leq 0.22$. We do not find evidence for a clear correlation between $f_{\rm esc}$ and galaxy stellar mass, suggesting it is not a primary indicator of leakage. Although larger samples are needed to further explore these trends, they suggest that it is entirely plausible that the low dust and metallicity galaxies found at z > 6 will display the $\langle f_{\rm esc}\rangle\geq0.1$ required to drive reionization., Comment: 15 pages (+ appendix), 12 figures, 1 table, Accepted for publication in Monthly Notices of the Royal Astronomical Society

Giant star-forming regions (clumps) are widespread features of galaxies at $z \approx 1-4$. Theory predicts that they can play a crucial role in galaxy evolution if they survive to stellar feedback for > 50 Myr. Numerical simulations show that clumps' survival depends on the stellar feedback recipes that are adopted. Up to date, observational constraints on both clumps' outflows strength and gas removal timescale are still uncertain. In this context, we study a line-emitting galaxy at redshift $z \simeq 3.4$ lensed by the foreground galaxy cluster Abell 2895. Four compact clumps with sizes $\lesssim$ 280 pc and representative of the low-mass end of clumps' mass distribution (stellar masses $\lesssim 2\times10^8\ {\rm M}_\odot$) dominate the galaxy morphology. The clumps are likely forming stars in a starbursting mode and have a young stellar population ($\sim$ 10 Myr). The properties of the Lyman-$\alpha$ (Ly$\alpha$) emission and nebular far-ultraviolet absorption lines indicate the presence of ejected material with global outflowing velocities of $\sim$ 200-300 km/s. Assuming that the detected outflows are the consequence of star formation feedback, we infer an average mass loading factor ($\eta$) for the clumps of $\sim$ 1.8 - 2.4 consistent with results obtained from hydro-dynamical simulations of clumpy galaxies that assume relatively strong stellar feedback. Assuming no gas inflows (semi-closed box model), the estimates of $\eta$ suggest that the timescale over which the outflows expel the molecular gas reservoir ($\simeq 7\times 10^8\ \text{M}_\odot$) of the four detected low-mass clumps is $\lesssim$ 50 Myr., Comment: 19 pages, 11 figures + Appendix. Published in MNRAS

We present Lyman continuum (LyC) radiation escape fraction $f_{\rm{esc}}$ measurements for 183 spectroscopically confirmed star-forming galaxies in the redshift range $3.11 < z < 3.53$ in the \textit{Chandra} Deep Field South. We use ground-based imaging to measure $f_{\rm{esc}}$, and use ground- and space-based photometry to derive galaxy physical properties using spectral energy distribution (SED) fitting. We additionally derive [O III]+H$\beta$ equivalent widths (that fall in the observed K band) by including nebular emission in the SED fitting. After removing foreground contaminants, we report the discovery of 11 new candidate LyC leakers, with absolute LyC escape fractions, $f_{\rm{esc}}$ in the range $0.14-0.85$. From non-detections, we place $1\sigma$ upper limits of $f_{\rm{esc}}<0.12$, where the Lyman-break selected galaxies have $f_{\rm{esc}} < 0.11$ and `blindly' discovered galaxies with no prior photometric selection have $f_{\rm{esc}}<0.13$. We find a slightly higher $1\sigma$ limit of $f_{\rm{esc}}<0.20$ for extreme emission line galaxies with rest-frame [O III]+H$\beta$ equivalent widths $>300$A. For candidate LyC leakers, we find a weak negative correlation between $f_{\rm{esc}}$ and galaxy stellar masses, no correlation between $f_{\rm{esc}}$ specific star-formation rates (sSFRs) and a positive correlation between $f_{\rm{esc}}$ and EW$_0$([O III]+H$\beta$). The weak/no correlations between stellar mass and sSFRs may be explained by misaligned viewing angles and/or non-coincident timescales of starburst activity and periods of high $f_{\rm{esc}}$. Alternatively, escaping radiation may predominantly occur in highly localised star-forming regions, or $f_{\rm{esc}}$ measurements may be impacted by stochasticity of the intervening neutral medium, obscuring any global trends with galaxy properties. These hypotheses have important consequences for models of reionisation., Comment: 19 pages, 11 figures, 1 appendix, published in Monthly Notices of the Royal Astronomical Society

Using the CLASH-VLT survey, we assembled an unprecedented sample of 1234 spectroscopically confirmed members in Abell~S1063, finding a dynamically complex structure at z_cl=0.3457 with a velocity dispersion \sigma_v=1380 -32 +26 km s^-1. We investigate cluster environmental and dynamical effects by analysing the projected phase-space diagram and the orbits as a function of galaxy spectral properties. We classify cluster galaxies according to the presence and strength of the [OII] emission line, the strength of the H$\delta$ absorption line, and colours. We investigate the relationship between the spectral classes of galaxies and their position in the projected phase-space diagram. We analyse separately red and blue galaxy orbits. By correlating the observed positions and velocities with the projected phase-space constructed from simulations, we constrain the accretion redshift of galaxies with different spectral types. Passive galaxies are mainly located in the virialised region, while emission-line galaxies are outside r_200, and are accreted later into the cluster. Emission-lines and post-starbursts show an asymmetric distribution in projected phase-space within r_200, with the first being prominent at Delta_v/sigma <~-1.5$, and the second at Delta_v/ sigma >~ 1.5, suggesting that backsplash galaxies lie at large positive velocities. We find that low-mass passive galaxies are accreted in the cluster before the high-mass ones. This suggests that we observe as passives only the low-mass galaxies accreted early in the cluster as blue galaxies, that had the time to quench their star formation. We also find that red galaxies move on more radial orbits than blue galaxies. This can be explained if infalling galaxies can remain blue moving on tangential orbits., Comment: 24 pages, 28 figures, 3 tables, Accepted for publication in Astronomy and Astrophysics

We present a strong lensing analysis of the galaxy cluster PSZ1 G311.65-18.48 (z=0.443) using multi-band observations with Hubble Space Telescope, complemented with VLT/MUSE spectroscopic data. The MUSE observations provide redshift estimates for the lensed sources and help reducing the mis-identification of the multiple images. Spectroscopic data are also used to measure the inner velocity dispersions of 15 cluster galaxies and calibrate the scaling relations to model the subhalo cluster component. The model is based on 62 multiple images grouped in 17 families belonging to 4 different sources. The majority of them are multiple images of compact stellar knots belonging to a single star-forming galaxy at z=2.3702. This source is strongly lensed by the cluster to form the Sunburst Arc system. To accurately reproduce all the multiple images, we build a parametric mass model, which includes both cluster-scale and galaxy-scale components. The resulting model has a r.m.s. separation between the model-predicted and the observed positions of the multiple images of only 0.14''. We conclude that PSZ1 G311.65-18.48 has a relatively round projected shape and a large Einstein radius (29'' for z_s = 2.3702), which could indicate that the cluster is elongated along the line of sight. The Sunburst Arc source is located at the intersection of a complex network of caustics, which explains why parts of the arc are imaged with unprecedented multiplicity (up to 12 times).

We investigate the strongly lensed (\mu x10-100) Lyman continuum (LyC) galaxy, dubbed Sunburst, at z=2.37, taking advantage of a new accurate model of the lens. A characterization of the intrinsic (delensed) properties of the galaxy yields a size of ~3 sq.kpc, a luminosity Muv=-20.3,and a stellar mass M~10^9 Msun;16% of the ultraviolet light is located in a 3 Myr old gravitationally-bound young massive star cluster (YMC) with an effective radius of Re~8 pc and a dynamical mass of ~10^7 Msun (similar to the stellar mass), from which LyC radiation is detected (\lambda < 912A). The inferred outflowing gas velocity (>300 km/s) exceeds the escape velocity of the star cluster. The resulting escape fraction of the ionizing radiation emerging from the Sunburst galaxy is >6-12%, whilst it is >46-93% if inferred from the YMC. 12 additional likely star clusters with 3~30%, which is consistent with the high \Gamma derived in local galaxies experiencing extreme gas physical conditions. The presence of the YMC influences the morphology (nucleation), photometry (photometric jumps) and spectroscopic output (nebular emission) of the entire galaxy. The de-lensed LyC and UV (1600A) magnitudes of the YMC are ~30.6 and ~26.9, whilst the galaxy has m1600~24.8. A relatively large rest-frame equivalent width of EWrest(Hb+[OIII]4959-5007)~450A emerges from the galaxy with the YMC contributing to ~30%. If O-type stars are mainly forged in star clusters, then such engines were the key ionizing agents during reionization and the increasing occurrence of high EW lines (Hb+[OIII]) observed at z>6.5 might be an indirect signature of a high \Gamma at reionization.Future facilities (like VLT/MAVIS or ELT), will probe bound clusters on moderately magnified (\mu<5-10) galaxies across cosmic epochs up to reionization[ABRIDGED], Comment: 17 pages (12 main body), 10 Figures and 1 Table. Accepted for publication in Astronomy & Astrophysics. Figure 6 shows the stellar cluster formation efficiency of the Sunburst galaxy; Figure C.1 shows the MUSE narrow field mode image of the counter-arc

VANDELS is an ESO Public Spectroscopic Survey designed to build a sample of high signal to noise, medium resolution spectra of galaxies at redshift between 1 and 6.5. Here we present the final Public Data Release of the VANDELS Survey, comprising 2087 redshift measurements. We give a detailed description of sample selection, observations and data reduction procedures. The final catalogue reaches a target selection completeness of 40% at iAB = 25. The high Signal to Noise ratio of the spectra (above 7 in 80% of the spectra) and the dispersion of 2.5{\AA} allowed us to measure redshifts with high precision, the redshift measurement success rate reaching almost 100%. Together with the redshift catalogue and the reduced spectra, we also provide optical mid-IR photometry and physical parameters derived through SED fitting. The observed galaxy sample comprises both passive and star forming galaxies covering a stellar mass range 8.3< Log(M*/Msolar)<11.7. All catalogues and spectra are accessible through the survey database (http://vandels.inaf.it) where all information can be queried interactively, and via the ESO Archive (https://www.eso.org/qi/)., Comment: 16 pages, 13 figures, accepted for publication in Astronomy & Astrophysics