Researcher profile

C. Grillo

C. Grillo contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2026arXiv

A Pristine Star-Forming Complex at z=4.19

We report the discovery of a faint (M_1700 ~ -12.2), oxygen-deficient strongly-lensed ionizing source -- dubbed LAP2 (Lensed And Pristine 2) -- at a spectroscopic redshift of z=4.19. LAP2 appears to be isolated and lies very close to the caustic produced by the lensing galaxy cluster Abell 2744. It was observed with the James Webb Space Telescope (JWST) NIRSpec MSA in prism mode as part of the UNCOVER program. The NIRSpec spectrum reveals prominent Lya emission (7.1 sigma), clear Ha emission (6.2 sigma), tentative Hb emission (2.8 sigma), and no detectable [OIII]4959,5007 (~ 7 times fainter than Ha). The inferred [OIII] 2 sigma upper limit corresponds to an R3 index <0.85 (assuming the Ha/Hb ~ 2.86 case~B recombination ratio), which, under high-ionization conditions, implies a metallicity of Z < 6 X 10^(-3) Z_sun. The combination of faint ultraviolet luminosity, large rest-frame Ha equivalent width (~ 650 A), and extremely compact size (< 10 pc) suggests that LAP2 is being caught in an early, pristine formation phase consistent with an instantaneous-burst scenario, with an estimated stellar mass of at most a few 10^4 Msun. Deep VLT/MUSE observations further reveal copious Lya emission forming an arclet that straddles the critical line. LAP2 thus joins the rare class of extremely metal-poor star-forming complexes that JWST has started to unveil at redshifts 3 - 7, and provides a rare glimpse into a still very poorly explored faint-luminosity regime.

preprint2025arXiv

Parametric strong lensing model of the galaxy cluster Abell 2390 from Euclid and MUSE observations

We present a new high precision parametric strong lensing total mass reconstruction of the Euclid Early Release Observations (ERO) galaxy cluster Abell 2390, at redshift z = 0.231. We include in this analysis 35 multiple images from 13 background sources, of which 25 are spectroscopically confirmed thanks to observations from the MUSE, spanning a redshift range from z = 0.535 to z = 4.877. After fully reanalysing the MUSE spectroscopy, we combine it with archival spectroscopic catalogues allowing us to select 65 secure cluster members. This sample is further complemented with 114 photometric member galaxies, identified within the Euclid VIS and NISP imaging down to magnitude H = 23. We also measure the stellar velocity dispersions for 22 cluster members, in order to calibrate the Faber-Jackson relation and hence the scaling relations for the subhalo mass components. We test and compare 11 total mass parametrisations of the galaxy cluster with increasing complexity. To do so, we employ the new parametric strong lensing modelling code Gravity.jl. Our best fit total mass parametrisation is characterised by a single large scale halo, 179 subhalo components, and an external shear term. The reference model yields a mean scatter between the model predicted and observed positions of the multiple images of 0.32 arcseconds. We are able to quantify the systematics arising from our modelling choices by taking advantage of all the different explored total mass parametrisations. When our results are compared with those from other lensing studies, we notice an overall agreement in the reconstructed cluster total mass profile in the outermost strong lensing regime: the discrepancy in the innermost region of the cluster (few kiloparsec from the brightest cluster galaxy, where few or no strong lensing features are observed) could possibly be ascribed to the different data and modelling choices.

preprint2025arXiv

Pushing JWST to the extremes: search and scrutiny of bright galaxy candidates at z$\simeq$15-30

We designed customized Lyman-break color selection techniques to identify galaxy candidates in the redshift ranges $15 \leq z \leq 20$ and $20 \leq z \leq 28$. The selection was performed on the ASTRODEEP-JWST multi-band catalogs of the CEERS, Abell-2744, JADES, NGDEEP, and PRIMER survey fields, covering a total area of $\sim0.2$ sq. deg. We identify five candidates at $15 \leq z \leq 20$, while no objects are found based on the $z\gtrsim20$ color selection criteria. Despite exhibiting a $>$1.5 mag break, all the objects display multimodal redshift probability distributions across different SED-fitting codes and methodologies. The alternative solutions correspond to poorly understood populations of low-mass quiescent or dusty galaxies at z$\sim$3-7. This conclusion is supported by the analysis of five F200W-dropout objects that we find to be interlopers on the basis of NIRSpec PRISM spectra: four dusty star-forming galaxies at z$\sim$2.2-6.6, and a passive galaxy at z=4.91 with log$(M_{\rm star}/{\rm M}_{\odot}) \lesssim$ 9. We measured the UV luminosity function under different assumptions on the contamination level within our sample. We find that if even a fraction of the candidates is indeed at $z\gtrsim15$, the resulting UV LF points to a very mild evolution compared to estimates at $z<15$, implying a significant tension with existing theoretical models. In particular, confirming our bright ($M_{\text{UV}}<-21$) candidates would require substantial revisions to the theoretical framework. In turn, if all these candidates will be confirmed to be interlopers, we conclude that future surveys may need ten times wider areas to select $M_{\text{UV}}\lesssim-20$ galaxies at $z>15$. Observations in the F150W and F200W filters at depths comparable to those in the NIRCam LW bands are also required to mitigate contamination from rare red objects at z$\lesssim$8.

preprint2022arXiv

Does the Fe L-shell blend bias abundance measurements in intermediate-temperature clusters?

In intermediate-temperature (T = 2.5 - 4.5 keV) galaxy clusters, abundance measurements are almost-equally driven by Fe K and L transitions, at $\sim$ 6.7 keV and 0.9 - 1.3 keV, respectively. While K-shell-derived measurements are considered reliable, the resolution of the currently available instrumentation, as well as our current knowledge of the atomic processes, makes the modelling of the L-line complex challenging, resulting in potential biases for abundance measurements. In this work, we study systematics related to the modelling of the Fe L-line complex that may influence iron-abundance measurements in the intermediate-temperature range. To this aim, we select a sample of three bright galaxy clusters, with long XMM-Newton observations available and temperature in the 2.5 - 4.5 keV range. We fit spectra extracted from concentric rings with APEC and APEC+APEC models, by alternatively excluding the L and K bands, and derive the fractional difference of the metal abundances, $ΔZ/Z$, as indication of the consistency between K- and L-shell-derived measurements. The $ΔZ/Z$ distribution is then studied as a function of the cluster radius, ring temperature and X-ray flux. The L-induced systematics, measured through an individual fit of each MOS and pn spectrum, remain constant at a 5 - 6% value in the whole 2.5 - 4.5 keV temperature range. Conversely, a joint fit of MOS and pn spectra leads to a slight excess of 1 - 2% in the above estimate. No significant dependence on the ring X-ray flux is highlighted. The measured 5 - 8% value indicates a modest contribution of the systematics to the derived iron abundances, giving confidence for future measurements. To date, these findings represent the best-achievable estimate of the systematics in analysis, while future microcalorimeters will significantly improve our understanding of the atomic processes underlying the Fe L emissions.

preprint2022arXiv

Early results from GLASS-JWST. XIII. A faint, distant, and cold brown dwarf

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~$μ$m wavelengths and the panchromatic 0.9--4.5~$μ$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.

preprint2022arXiv

Exploring the physical properties of lensed star-forming clumps at $2\lesssim z \lesssim6$

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

preprint2022arXiv

The GLASS James Webb Space Telescope Early Release Science Program. I. Survey Design and Release Plans

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, &#34;what sources ionized the universe and when?&#34; and &#34;how do baryons cycle through galaxies?&#34;, 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.

preprint2021arXiv

ALMA Lensing Cluster Survey: an ALMA galaxy signposting a MUSE galaxy group at z=4.3 behind &#39;El Gordo&#39;

We report the discovery of a Multi Unit Spectroscopic Explorer (MUSE) galaxy group at z=4.32 lensed by the massive galaxy cluster ACT-CL J0102-4915 (aka El Gordo) at z=0.87, associated with a 1.2 mm source which is at a 2.07+/-0.88 kpc projected distance from one of the group galaxies. Three images of the whole system appear in the image plane. The 1.2 mm source has been detected within the Atacama Large Millimetre/submillimetre Array (ALMA) Lensing Cluster Survey (ALCS). As this ALMA source is undetected at wavelengths lambda < 2 microns, its redshift cannot be independently determined, however, the three lensing components indicate that it belongs to the same galaxy group at z=4.32. The four members of the MUSE galaxy group have low to intermediate stellar masses (~ 10^7-10^{10} Msun) and star formation rates (SFRs) of 0.4-24 Msun/yr, resulting in high specific SFRs (sSFRs) for two of them, which suggest that these galaxies are growing fast (with stellar-mass doubling times of only ~ 2x10^7 years). This high incidence of starburst galaxies is likely a consequence of interactions within the galaxy group, which is compact and has high velocity dispersion. Based on the magnification-corrected sub-/millimetre continuum flux density and estimated stellar mass, we infer that the ALMA source is classified as an ordinary ultra-luminous infrared galaxy (with associated dust-obscured SFR~200-300 Msun/yr) and lies on the star-formation main sequence. This reported case of an ALMA/MUSE group association suggests that some presumably isolated ALMA sources are in fact signposts of richer star-forming environments at high redshifts.

preprint2021arXiv

Galaxy cluster strong lensing cosmography: cosmological constraints from a sample of regular galaxy clusters

Cluster strong lensing cosmography is a promising probe of the background geometry of the Universe and several studies have emerged, thanks to the increased quality of observations using space and ground-based telescopes. For the first time, we use a sample of five cluster strong lenses to measure the values of cosmological parameters and combine them with those from classical probes. In order to assess the degeneracies and the effectiveness of strong-lensing cosmography in constraining the background geometry of the Universe, we adopt four cosmological scenarios. We find good constraining power on the total matter density of the Universe ($Ω_{\rm m}$) and the equation of state of the dark energy parameter $w$. For a flat $w$CDM cosmology, we find $Ω_{\rm m} = 0.30_{-0.11}^{+0.09}$ and $w=-1.12_{-0.32}^{+0.17}$ from strong lensing only. Interestingly, we show that the constraints from the Cosmic Microwave Background (CMB) are improved by factors of 2.5 and 4.0 on $Ω_{\rm m}$ and $w$, respectively, when combined with our posterior distributions in this cosmological model. In a scenario where the equation of state of dark energy evolves with redshift, the strong lensing constraints are compatible with a cosmological constant (i.e. $w=-1$). In a curved cosmology, our strong lensing analyses can accommodate a large range of values for the curvature of the Universe of $Ω_{\rm k}=0.28_{-0.21}^{+0.16}$. In all cosmological scenarios, we show that our strong lensing constraints are complementary and in good agreement with measurements from the CMB, baryon acoustic oscillations and Type Ia supernovae. Our results show that cluster strong lensing cosmography is a potentially powerful probe to be included in the cosmological analyses of future surveys.

preprint2021arXiv

High star cluster formation efficiency in the strongly lensed Sunburst Lyman-continuum galaxy at z=2.37

We investigate the strongly lensed (μ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 (λ< 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<Re<20 pc are identified in the galaxy from which we derive a cluster formation efficiency Γ>~30%, which is consistent with the high Γ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 Γat reionization.Future facilities (like VLT/MAVIS or ELT), will probe bound clusters on moderately magnified (μ<5-10) galaxies across cosmic epochs up to reionization[ABRIDGED]

preprint2021arXiv

The Young Supernova Experiment: Survey Goals, Overview, and Operations

Time domain science has undergone a revolution over the past decade, with tens of thousands of new supernovae (SNe) discovered each year. However, several observational domains, including SNe within days or hours of explosion and faint, red transients, are just beginning to be explored. Here, we present the Young Supernova Experiment (YSE), a novel optical time-domain survey on the Pan-STARRS telescopes. Our survey is designed to obtain well-sampled $griz$ light curves for thousands of transient events up to $z \approx 0.2$. This large sample of transients with 4-band light curves will lay the foundation for the Vera C. Rubin Observatory and the Nancy Grace Roman Space Telescope, providing a critical training set in similar filters and a well-calibrated low-redshift anchor of cosmologically useful SNe Ia to benefit dark energy science. As the name suggests, YSE complements and extends other ongoing time-domain surveys by discovering fast-rising SNe within a few hours to days of explosion. YSE is the only current four-band time-domain survey and is able to discover transients as faint $\sim$21.5 mag in $gri$ and $\sim$20.5 mag in $z$, depths that allow us to probe the earliest epochs of stellar explosions. YSE is currently observing approximately 750 square degrees of sky every three days and we plan to increase the area to 1500 square degrees in the near future. When operating at full capacity, survey simulations show that YSE will find $\sim$5000 new SNe per year and at least two SNe within three days of explosion per month. To date, YSE has discovered or observed 8.3% of the transient candidates reported to the International Astronomical Union in 2020. We present an overview of YSE, including science goals, survey characteristics and a summary of our transient discoveries to date.

preprint2020arXiv

An accurate strong lensing model of the Abell 2163 core

Abell 2163 at $z \simeq 0.201$ is one of the most massive galaxy clusters known, very likely in a post-merging phase. Data from several observational windows suggest a complex mass structure with interacting subsystems, which makes the reconstruction of a realistic merging scenario very difficult. A missing key element in this sense is unveiling the cluster mass distribution at high resolution. We perform such a reconstruction of the cluster inner total mass through a strong lensing model based on new spectroscopic redshift measurements. We use data from the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT) to confirm 12 multiple images of 4 sources with redshift values from 1.16 to 2.72. We also discover four new multiple images and identify 29 cluster members and 35 foreground and background sources. The resulting galaxy member and image catalogs are used to build five cluster total mass models. The fiducial model consists of 111 small-scale subhalos plus a diffuse component, which is centered $\sim2$ arcseconds away from the BCG belonging to the east Abell 2163 subcluster. We confirm that the latter is well represented by a single, large-scale mass component. Its strong elongation towards a second (west) subcluster confirms the existence of a preferential axis, corresponding to the merging direction. From the fiducial model, we extrapolate the cumulative projected total mass profile and measure a value of $M(<300\,$kpc$) = 1.43^{+0.07}_{-0.06}\times 10^{14}\,$M$_{\odot}$, which has a significantly reduced statistical error compared with previous estimates, thanks to the inclusion of the spectroscopic redshifts. Our strong lensing results are very accurate: the model-predicted positions of the multiple images are, on average, only $0.15$ arcseconds away from the observed ones.

preprint2020arXiv

Candidate Population III stellar complex at z=6.629 in the MUSE Deep Lensed Field

We discovered a strongly lensed (μ>40) Lya emission at z=6.629 (S/N~18) in the MUSE Deep Lensed Field (MDLF) targeting the Hubble Frontier Field galaxy cluster MACS~J0416. Dedicated lensing simulations imply that the Lya emitting region necessarily crosses the caustic. The arc-like shape of the Lya extends 3 arcsec on the observed plane and is the result of two merged multiple images, each one with a de-lensed Lya luminosity L<~2.8 x 10^(40) erg/s arising from a confined region (< 150 pc effective radius). A spatially unresolved HST counterpart is barely detected at S/N~2 after stacking the near-infrared bands, corresponding to an observed(intrinsic) magnitude m_(1500)>~30.8(>~35.0). The inferred rest-frame Lya equivalent width is EWo > 1120 A if the IGM transmission is T(IGM)<0.5. The low luminosities and the extremely large Lya EWo match the case of a Population~III star complex made of several dozens stars (~ 10^4 Msun) which irradiate a HII region crossing the caustic. While the Lya and stellar continuum are among the faintest ever observed at this redshift, the continuum and the Lya emissions could be affected by differential magnification, possibly biasing the EWo estimate. The aforementioned tentative HST detection tend to favor a large EWo, making such a faint Pop~III candidate a key target for the James Webb Space Telescope and Extremely Large Telescopes.

preprint2020arXiv

CLASH-VLT: a full dynamical reconstruction of the mass profile of Abell S1063 from 1 kpc out to the virial radius

Context: The shape of the mass density profiles of cosmological halos informs us of the nature of DM and DM-baryons interactions. Previous estimates of the inner slope of the mass density profiles of clusters of galaxies are in opposition to predictions derived from numerical simulations of CDM. Aims: We determine the inner slope of the DM density profile of a massive cluster of galaxies, Abell S1063 (RXC J2248.7$-$4431) at $z=0.35$, with a dynamical analysis based on an extensive spectroscopic campaign carried out with the VIMOS and MUSE spectrographs at the ESO VLT. This new data set provides an unprecedented sample of 1234 spectroscopic members, 104 of which are located in the cluster core (R < 200 kpc), extracted from the MUSE integral field spectroscopy. The latter also allows the stellar velocity dispersion profile of the brightest cluster galaxy (BCG) to be measured out to 40 kpc. Methods: We used an upgraded version of the MAMPOSSt technique to perform a joint maximum likelihood fit to the velocity dispersion profile of the BCG and to the velocity distribution of cluster member galaxies over a radial range from 1 kpc to the virial radius (~ 2.7 Mpc). Results: We find a value of $γ_{DM} =0.99 \pm 0.04$ for the inner logarithmic slope of the DM density profile after marginalizing over all the other parameters. The newly determined dynamical mass profile is found to be in excellent agreement with the mass density profiles obtained from the independent X-ray hydrostatic analysis based on deep Chandra data, as well as the strong and weak lensing analyses. Our value of γ_{DM} is in very good agreement with predictions from cosmological CDM simulations. We will extend our analysis to more clusters in future works. If confirmed on a larger cluster sample, our result makes this DM model more appealing than alternative models.

preprint2020arXiv

On the accuracy of time-delay cosmography in the Frontier Fields Cluster MACS J1149.5+2223 with supernova Refsdal

We study possible systematic effects on the values of the cosmological parameters measured through strong lensing analyses of the HFF galaxy cluster MACS J1149.5+2223. We use the observed positions of a large set of spectroscopically selected multiple images, including those of supernova &#34;Refsdal&#34; with their estimated time delays. Starting from our reference model in a flat $Λ$CDM cosmology, published in Grillo et al. (2018), we confirm the relevance of the longest measurable time delay, between SX and S1, and an approximately linear relation between its value and that of $H_{0}$. We perform true blind tests by considering a range of time delays around its original estimate of $345 \pm 10$ days, as an accurate measurement of this time delay was not known at the time of analysis and writing. We investigate separately the impact of a constant sheet of mass at the cluster redshift, of a power-law profile for the mass density of the cluster main halo and of some scatter in the cluster member scaling relations. Remarkably, we find that these systematic effects do not introduce a significant bias on the inferred values of $H_{0}$ and $Ω_{\rm m}$ and that the statistical uncertainties dominate the total error budget: a 3% uncertainty on the time delay of image SX translates into approximately 6% and 40% (including both statistical and systematic $1σ$) uncertainties for $H_{0}$ and $Ω_{\rm m}$, respectively. Furthermore, our model accurately reproduces the extended surface brightness distribution of the supernova host, covering more than $3 \times 10^{4}$ $HST$ pixels. We also present the interesting possibility of measuring the value of the equation-of-state parameter $w$ of the dark energy density, currently with a 30% uncertainty. We conclude that time-delay cluster lenses have the potential to become soon an alternative and competitive cosmological probe.

preprint2020arXiv

The MUSE Deep Lensed Field on the Hubble Frontier Field MACS~J0416

Context: A census of faint and tiny star forming complexes at high redshift is key to improving our understanding of reionizing sources, galaxy growth and the formation of globular clusters. Aims: We present the MUSE Deep Lensed Field (MDLF) program. Methods: We describe Deep MUSE observations of 17.1 hours integration on a single pointing over the Hubble Frontier Field galaxy cluster MACS~J0416. Results: We confirm spectroscopic redshifts for all 136 multiple images of 48 source galaxies at 0.9<z<6.2. Within those galaxies, we securely identify 182 multiple images of 66 galaxy components that we use to constrain our lens model. We identify 116 clumps belonging to background high-z galaxies; the majority of them are multiple images and span magnitude, size and redshift intervals of [-18,-10], [~400-3] parsec and 1<z<6.6, respectively, with the most magnified ones probing possible single gravitationally bound star clusters. The depth of the MDLF combined with lensing magnification lead us to reach a detection limit for unresolved emission lines of a few 10$^{-20}$ erg/s/cm2, after correction for lensing magnification. Ultraviolet high-ionization metal lines (and HeII1640) are detected with S/N>10 for individual objects down to de-lensed magnitude 28-30 suggesting that they are common in such faint sources. Conclusions:Deep MUSE observations, in combination with existing HST imaging, allowed us to:(1) confirm redshifts for extremely faint high-z sources;(2) peer into their internal clumps (down to 100-200 pc scale);(3) in some cases break down such clumps into bound star clusters (<20 pc scale);(4) double the number of constraints for the lens model,reaching an unprecedented set of 182 bona-fide multiple images and confirming up to 213 galaxy cluster members. These results demonstrate the power that JWST and future ELTs will have when combined to study gravitational telescopes.[abridged]