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 "Refsdal" with their estimated time delays. Starting from our reference model in a flat $\Lambda$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 $\Omega_{\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\sigma$) uncertainties for $H_{0}$ and $\Omega_{\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.

We discovered Bowen emission arising from a strongly lensed ($\mu>20$) transient stellar object hosted in the Sunburst arc at z=2.37. Unique ultraviolet lines emerge at the location of the transient. In particular, narrow ($\sigma_v$ ~ 40 km/s) ionisation lines of several chemical species (Fe, C, Si) fluoresce after being exposed to Lyman series (e.g., Ly$\alpha$, Ly$\beta$) of H, HeI and HeII that pump selectively their atomic levels. Similarities with the Resonance-Enhanced Two Photon-Ionisation (RETPI) spectral features observed in the circum-stellar Weigelt blobs of Eta~Carinae are observed. Data from VLT/MUSE, X-Shooter and ESPRESSO observations (the latter placed at the focus of the four UTs) at increasing spectral resolution of R=2500, 11400 and R=70000, respectively, confirm the lasing action since at least 3.3 years (~ 1 year rest-frame), and probe the circum-stellar dense gas condensations in radiation-rich conditions. We discuss the physical origin of the transient event, which remains unclear. We expect such transient events (including also supernova or impostors) will be easily recognised with ELTs thanks to high angular resolution provided by adaptive optics and large collecting area, especially in modest ($\mu < 3$) magnification regime.

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 $\gamma_{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 \gamma_{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.

We discovered a strongly lensed (\mu >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.

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.

We carry out a detailed strong lensing analysis of a sub-sample of eight galaxy clusters of the CLASH survey, in the redshift range of $ z_{\rm cluster} = [0.23-0.59]$, using extensive spectroscopic information, primarily MUSE archival data complemented with CLASH-VLT redshift measurements. Different models are tested in each cluster depending on the complexity of its mass distribution and on the number of detected multiple images. Four clusters show more than five spectroscopically confirmed multiple image families. In this sample, we do not make use of families that are only photometrically identified, in order to reduce model degeneracies and systematics due to the potential misidentifications of some multiple images. We present spectroscopic confirmation of 27 multiply lensed sources, with no previous spectroscopic measurements, spanning over the redshift range of $z_{\rm src}=[0.7-6.1]$. Moreover, we confirm an average of $48$ galaxy members in the core of each cluster, thanks to the high efficiency and large field of view of MUSE. Despite having different properties (i.e., number of mass components, total mass, redshift, etc), the projected total mass and mass density profiles of all clusters have very similar shapes, when rescaled by independent measurements of $M_{200c}$ and $R_{200c}$. Specifically, we measure the mean value of the projected total mass of our cluster sample within 10 (20)% of $R_{200c}$ to be 0.13 (0.32) of $M_{200c}$, with a remarkably small scatter of 5 (6)%. Furthermore, the large number of high-z sources and the precise magnification maps derived in this work for four clusters add up to the sample of high-quality gravitational telescopes to be used to study the faint and distant Universe. The strong lensing models and the full redshift catalogues from MUSE are publicly available.

We search for bright Ly$\rm \alpha$ emitters among Spitzer SMUVS galaxies at z > 2.9 with homogeneous MUSE data. Although it only covers a small region of COSMOS, MUSE has the unique advantage of providing spectral information over the entire field, without the need of target pre-selection. This gives an unbiased detection of all the brightest Ly$\rm \alpha$ emitters among SMUVS sources, which by design are stellar-mass selected galaxies. Within the studied area, ~14% of the SMUVS galaxies at z > 2.9 have Ly$\rm \alpha$ fluxes F$\rm _\lambda$ > 7 x 10$^{-18}$ erg s$^{-1}$ cm$^{-2}$. These Ly$\rm \alpha$ emitters are characterized by three types of emission, 47% show a single line profile, 19% present a double peak or a blue bump and 31% show a red tail. One object (3%) shows both a blue bump and a red tail. We also investigate the spectral energy distribution (SED) properties of the SMUVS MUSE-detected galaxies and MUSE non-detections. After stellar-mass matching both populations, we find that MUSE detected galaxies have generally lower extinction than SMUVS-only objects, while there is no clear intrinsic difference in the mass and age distributions. For the MUSE-detected SMUVS galaxies, we compare the instantaneous SFR lower limit given by Ly$\rm \alpha$ flux with its past average derived from SED fitting, and find evidence for rejuvenation in some of our oldest objects. We also study the spectra of those Ly$\rm \alpha$ emitters which are not detected in SMUVS in the same field. We find different distributions of the emission line profiles, which could be ascribed to the fainter Ly$\rm \alpha$ luminosities of the MUSE-only sources and an intrinsically different mass distribution. Finally, we search for the presence of galaxy associations. MUSE's integral coverage is 20 times more likely to find associations than all other existing spectral data in COSMOS, biased by target pre-selection.

We present an improved determination of the total mass distribution of three CLASH/HFF massive clusters, MACS J1206.2-0847 (z=0.44), MACS J0416.1-2403 (z=0.40), Abell S1063 (z=0.35). We specifically reconstruct the sub-halo mass component with robust stellar kinematics information of cluster galaxies, in combination with precise strong lensing models based on large samples of spectroscopically identified multiple images. We use VLT/MUSE integral-field spectroscopy in the cluster cores to measure the stellar velocity dispersion, $\sigma$, of 40-60 member galaxies per cluster, covering 4-5 magnitudes to $m_{F160W}\simeq 21.5$. We verify the robustness and quantify the accuracy of the velocity dispersion measurements with extensive spectral simulations. With these data, we determine the normalization and slope of the galaxy $L\mbox{-}\sigma$ Faber-Jackson relation in each cluster and use these parameters as a prior for the scaling relations of the sub-halo population in the mass distribution modeling. When compared to our previous lens models, the inclusion of member galaxies' kinematics provides a similar precision in reproducing the positions of the multiple images. However, the inherent degeneracy between the central effective velocity dispersion, $\sigma_0$, and truncation radius, $r_{cut}$, of sub-halos is strongly reduced, thus significantly alleviating possible systematics in the measurements of sub-halo masses. The three independent determinations of the $\sigma_0\mbox{-}r_{cut}$ scaling relation in each cluster are found to be fully consistent, enabling a statistical determination of sub-halo sizes as a function of $\sigma_0$, or halo masses. We derive galaxy central velocity dispersion functions of the three clusters and found them in agreement with each other. Sub-halo mass functions determined with this method can be compared with those obtained from cosmological simulations.

We report the discovery of a galaxy overdensity around a Compton-thick Fanaroff-Riley type II (FRII) radio galaxy at z=1.7 in the deep multiband survey around the z=6.3 QSO SDSS J1030+0524. Based on a 6hr VLT/MUSE and on a 4hr LBT/LUCI observation, we identify at least eight galaxy members in this structure with spectroscopic redshift z=1.687-1.699, including the FRII galaxy at z=1.699. Most of the identified overdensity members are blue, compact galaxies that are actively forming stars at rates of $\sim$8-60 $M_{\odot}$ yr$^{-1}$. Based on a 500ks Chandra ACIS-I observation we found that the FRII nucleus hosts a luminous QSO ($L_{2-10keV}=1.3\times10^{44}$ erg s$^{-1}$, intrinsic and rest-frame) that is obscured by Compton-thick absorption ($N_H=1.5\pm0.6\times 10^{24}$ cm$^{-2}$). Our Chandra observation, the deepest so far for a distant FRII within a galaxy overdensity, revealed significant diffuse X-ray emission within the region covered by the overdensity. In particular, X-ray emission extending for $\sim$240 kpc is found around the Eastern lobe of the FRII. Four out of the six MUSE star forming galaxies in the overdensity are distributed in an arc-like shape at the edge of this diffuse X-ray emission. The probability of observing by chance four out of the six $z=1.7$ sources at the edge of the diffuse emission is negligible. In addition, these four galaxies have the highest specific star formation rates of the MUSE galaxies in the overdensity and lie above the main sequence of field galaxies of equal stellar mass at z=1.7. We propose that the diffuse X-rays originate from an expanding bubble of gas that is shock-heated by the FRII jet, and that star formation is promoted by the compression of the cold interstellar medium of the galaxies around the bubble, which may be remarkable evidence of positive AGN feedback on cosmological scales. [shortened version]

We present a VLT/X-Shooter spectroscopy of the Lyman continuum (LyC) emitting galaxy 'Ion2' at z=3.2121 and compare it to that of the recently discovered strongly lensed LyC-emitter at z=2.37, known as the 'Sunburst' arc. Three main results emerge from the X-Shooter spectrum: (a) the Lya has three distinct peaks with the central one at the systemic redshift, indicating a ionised tunnel through which both Lya and LyC radiation escape; (b) the large O32 oxygen index ([OIII]4959-5007 / [OII]3727-3729) of 9.18(-1.32/+1.82) is compatible to those measured in local (z~0.4) LyC leakers; (c) there are narrow nebular high-ionisation metal lines with \sigma_v < 20 km/s, which confirms the presence of young hot, massive stars. The HeII1640 appears broad, consistent with a young stellar component including Wolf-Rayet stars. Similarly, the Sunburst LyC-emitter shows a triple-peaked Lya profile and from VLT/MUSE spectroscopy the presence of spectral features arising from young hot and massive stars. The strong lensing magnification, (\mu > 20), suggests that this exceptional object is a gravitationally-bound star cluster observed at a cosmological distance, with a stellar mass M <~ 10^7 Msun and an effective radius smaller than 20 pc. Intriguingly, sources like Sunburst but without lensing magnification might appear as Ion2-like galaxies, in which unresolved massive star clusters dominate the ultraviolet emission. This work supports the idea that dense young star clusters can contribute to the ionisation of the IGM through holes created by stellar feedback.

We report a multi-wavelength study of the massive ($M_{\star} \gtrsim 10^{11} \rm{M}_{\odot}$), $z\sim 2$ star-forming galaxy GMASS 0953, which hosts an obscured AGN. We combined near-infrared observations of the GNIRS, SINFONI and KMOS spectrographs to study the kinematics of the [O III]$\lambda 5007$ and H$\alpha$ emission lines. Our analysis shows that GMASS 0953 may host an ionized disc extending up to 13 kpc, which rotates at a velocity of $V_{\rm{ion}} = 203^{+17}_{-20} {\ \rm kms}^{-1}$ at the outermost radius. Evidence of rotation on a smaller scale ($R \sim 1$ kpc) arises from the CO(J=6-5) line. The central velocity $V_{\rm{CO}} = 320^{+ 92}_{-53} {\ \rm kms}^{-1}$ traced by the molecular gas is higher than $V_{\rm{ion}}$, suggesting that the galaxy harbors a multi-phase disc with a rotation curve that peaks in the very central regions. The galaxy appears well located on the $z = 0$ baryonic Tully-Fisher relation. We also discuss the possibility that the [O III]$\lambda 5007$ and H$\alpha$ velocity gradients are due to a galactic-scale wind. Besides, we found evidence of an AGN-driven outflow traced by a broad blueshifted wing affecting the [O III]$\lambda 5007$ line, which presents a velocity offset $\Delta v = -535 \pm 152 {\ \rm kms}^{-1}$ from the systemic velocity. Because of the short depletion timescale ($\tau_{\rm{dep}}\sim 10^8$ yr) due to gas ejection and gas consumption by star formation activity, GMASS 0953 may likely evolve into a passive galaxy. However, the role of the AGN in depleting the gas reservoir of the galaxy is quite unclear because of the uncertainties affecting the outflow rate.

Based on HST and MUSE data, we probe the stellar and gas properties (i.e. kinematics, stellar mass, star formation rate) of the radio-loud brightest cluster galaxy (BCG) located at the centre of the X-ray luminous cool core cluster Abell 2667 (z = 0.2343). The bi-dimensional modelling of the BCG surface brightness profile reveals the presence of a complex system of substructures extending all around the galaxy. Clumps of different size and shape plunged into a more diffuse component constitute these substructures, whose intense 'blue' optical colour hints to the presence of a young stellar population. Our results depict the BCG as a massive (M_star ~ 1.38 x 10^11 M_sun) dispersion-supported spheroid (v_star < 150 km/s, sigma_0 ~ 216 km/s) hosting an active supermassive black hole (M_SMBH ~ 3.8 x 10^9 M_sun) whose optical features are typical of low ionisation nuclear emission line regions. Although the velocity pattern of the stars in the BCG is irregular, the stellar kinematics in the regions of the clumps show a positive velocity of ~ 100 km/s, similarly to the gas component. An analysis of the mechanism giving rise to the observed lines in the clumps through empirical diagnostic diagrams points out that the emission is composite, suggesting the contribution from both star formation and AGN. We conclude our analysis describing how scenarios of both chaotic cold accretion and merging with a gas-rich disc galaxy can efficaciously explain the phenomena the BCG is undergoing.

We report on a superdense star-forming region with an effective radius (R_e) smaller than 13 pc identified at z=6.143 and showing a star-formation rate density \Sigma_SFR~1000 Msun/yr/kpc2 (or conservatively >300 Msun/yr/kpc2). Such a dense region is detected with S/N>40 hosted by a dwarf extending over 440 pc, dubbed D1 (Vanzella et al. 2017b). D1 is magnified by a factor 17.4+/-5.0 behind the Hubble Frontier Field galaxy cluster MACS~J0416 and elongated tangentially by a factor 13.2+/-4.0 (including the systematic errors). The lens model accurately reproduces the positions of the confirmed multiple images with a r.m.s. of 0.35", and the tangential stretch is well depicted by a giant multiply-imaged Lya arc. D1 is part of an interacting star-forming complex extending over 800 pc. The SED-fitting, the very blue ultraviolet slope (\beta ~ -2.5, F(\lambda) ~ \lambda^\beta) and the prominent Lya emission of the stellar complex imply that very young (< 10-100 Myr), moderately dust-attenuated (E(B-V)<0.15) stellar populations are present and organised in dense subcomponents. We argue that D1 (with a stellar mass of 2 x 10^7 Msun) might contain a young massive star cluster of M < 10^6 Msun and Muv~-15.6 (or m_uv=31.1), confined within a region of 13 pc, and not dissimilar from some local super star clusters (SSCs). The ultraviolet appearance of D1 is also consistent with a simulated local dwarf hosting a SSC placed at z=6 and lensed back to the observer. This compact system fits into some popular globular cluster formation scenarios. We show that future high spatial resolution imaging (e.g., E-ELT/MAORY-MICADO and VLT/MAVIS) will allow us to spatially resolve light profiles of 2-8 pc.

We present a multi-wavelength study of the massive ($M_{200\textrm{c}} \approx 1$-$2 \times 10^{15} M_\odot$) galaxy clusters RXC J2248.7$-$4431, MACS J0416.1$-$2403, and MACS J1206.2$-$0847 at $z \approx 0.4$. Using the X-ray surface brightness of the clusters from deep Chandra data to model their hot gas, we are able to disentangle this mass term from the diffuse dark matter in our new strong-lensing analysis, with approximately $50$-$100$ secure multiple images per cluster, effectively separating the collisional and collisionless mass components of the clusters. At a radial distance of $10\%$ of $R_{200\textrm{c}}$ (approximately $200$ kpc), we measure a projected total mass of $(0.129 \pm 0.001)$, $(0.131 \pm 0.001)$ and $(0.137 \pm 0.001)\times M_{200\textrm{c}}$, for RXC J2248, MACS J0416 and MACS J1206, respectively. These values are surprisingly similar, considering the large differences in the merging configurations, and, as a consequence, in the mass models of the clusters. Interestingly, at the same radii, the hot gas over total mass fractions differ substantially, ranging from $0.082 \pm 0.001$ to $0.133 \pm 0.001$, reflecting the various dynamical states of the clusters. Moreover, we do not find a statistically significant offset between the positions of the peak of the diffuse dark matter component and of the BCG in the more complex clusters of the sample. We extend to this sample of clusters previous findings of a number of massive sub-halos higher than in numerical simulations. These results highlight the importance of a proper separation of the different mass components to study in detail the properties of dark matter in galaxy clusters.

We present the first weak lensing calibration of $\mu_{\star}$, a new galaxy cluster mass proxy corresponding to the total stellar mass of red and blue members, in two cluster samples selected from the SDSS Stripe 82 data: 230 redMaPPer clusters at redshift $0.1\leq z<0.33$ and 136 Voronoi Tessellation (VT) clusters at $0.1 \leq z < 0.6$. We use the CS82 shear catalog and stack the clusters in $\mu_{\star}$ bins to measure a mass-observable power law relation. For redMaPPer clusters we obtain $M_0 = (1.77 \pm 0.36) \times 10^{14}h^{-1} M_{\odot}$, $\alpha = 1.74 \pm 0.62$. For VT clusters, we find $M_0 = (4.31 \pm 0.89) \times 10^{14}h^{-1} M_{\odot}$, $\alpha = 0.59 \pm 0.54$ and $M_0 = (3.67 \pm 0.56) \times 10^{14}h^{-1} M_{\odot}$, $\alpha = 0.68 \pm 0.49$ for a low and a high redshift bin, respectively. Our results are consistent, internally and with the literature, indicating that our method can be applied to any cluster finding algorithm. In particular, we recommend that $\mu_{\star}$ be used as the mass proxy for VT clusters. Catalogs including $\mu_{\star}$ measurements will enable its use in studies of galaxy evolution in clusters and cluster cosmology.

Realizing Refsdal's original idea from 1964, we present estimates of the Hubble constant that are complementary to and potentially competitive with those of other cosmological probes. We use the observed positions of 89 multiple images, with extensive spectroscopic information, from 28 background sources and the measured time delays between the images S1-S4 and SX of supernova "Refsdal" ($z = 1.489$), which were obtained thanks to Hubble Space Telescope (HST) deep imaging and Multi Unit Spectroscopic Explorer (MUSE) data. We extend the strong lensing modeling of the Hubble Frontier Fields (HFF) galaxy cluster MACS J1149.5$+$2223 ($z = 0.542$), published by Grillo et al. (2016), and explore different $\Lambda$CDM models. Taking advantage of the lensing information associated to the presence of very close pairs of multiple images at various redshifts and to the extended surface brightness distribution of the SN Refsdal host, we can reconstruct the total mass density profile of the cluster very precisely. The combined dependence of the multiple image positions and time delays on the cosmological parameters allows us to infer the values of $H_{0}$ and $\Omega_{\rm m}$ with relative (1$\sigma$) statistical errors of, respectively, 6% (7%) and 31% (26%) in flat (general) cosmological models, assuming a conservative 3% uncertainty on the final time delay of image SX and, remarkably, no priors from other cosmological experiments. Our best estimate of $H_{0}$, based on the model described in this work, will be presented when the final time-delay measurement becomes available. Our results show that it is possible to utilize time delays in lens galaxy clusters as an important alternative tool for measuring the expansion rate and the geometry of the Universe.

We report on the serendipitous discovery of a z=4.0, M1500=-22.20 star-forming galaxy (Ion3) showing copious Lyman continuum (LyC) leakage (~60% escaping), a remarkable multiple peaked Lya emission, and significant Lya radiation directly emerging at the resonance frequency. This is the highest redshift confirmed LyC emitter in which the ionising and Lya radiation possibly share a common ionised cavity (with N_HI<10^17.2 cm^-2). Ion3 is spatially resolved, it shows clear stellar winds signatures like the P-Cygni NV1240 profile, and has blue ultraviolet continuum (\beta = -2.5 +/- 0.25, F_\lambda~ \lambda^\beta) with weak low-ionisation interstellar metal lines. Deep VLT/HAWKI Ks and Spitzer/IRAC 3.6um and 4.5um imaging show a clear photometric signature of the Halpha line with equivalent width of 1000A rest-frame emerging over a flat continuum (Ks-4.5um ~ 0). From the SED fitting we derive a stellar mass of 1.5x10^9 Msun, SFR of 140 Msun/yr and age of ~10 Myr, with a low dust extinction, E(B-V)< 0.1, placing the source in the starburst region of the SFR-M^* plane. Ion3 shows similar properties of another LyC emitter previously discovered (z=3.21, Ion2, Vanzella et al. 2016). Ion3 (and Ion2) represents ideal high-redshift reference cases to guide the search for reionising sources at z>6.5 with JWST.

We analyse the Kormendy relations (KRs) of the two Frontier Fields clusters, Abell S1063, at z = 0.348, and MACS J1149.5+2223, at z = 0.542, exploiting very deep Hubble Space Telescope photometry and VLT/MUSE integral field spectroscopy. With this novel dataset, we are able to investigate how the KR parameters depend on the cluster galaxy sample selection and how this affects studies of galaxy evolution based on the KR. We define and compare four different galaxy samples according to: (a) S\'ersic indices: early-type ('ETG'), (b) visual inspection: 'ellipticals', (c) colours: 'red', (d) spectral properties: 'passive'. The classification is performed for a complete sample of galaxies with m$_{\textit{F814W}} \le$ 22.5 ABmag (M$_{*}$ $\gtrsim 10^{10.0}$ M$_{\odot}$). To derive robust galaxy structural parameters, we use two methods: (1) an iterative estimate of structural parameters using images of increasing size, in order to deal with closely separated galaxies and (2) different background estimations, to deal with the Intracluster light contamination. The comparison between the KRs obtained from the different samples suggests that the sample selection could affect the estimate of the best-fitting KR parameters. The KR built with ETGs is fully consistent with the one obtained for ellipticals and passive. On the other hand, the KR slope built on the red sample is only marginally consistent with those obtained with the other samples. We also release the photometric catalogue with structural parameters for the galaxies included in the present analysis.

We present the results from a $\sim500$ ks Chandra observation of the $z=6.31$ QSO SDSS J1030+0524. This is the deepest X-ray observation to date of a $z\sim6$ QSO. The QSO is detected with a total of 125 net counts in the full ($0.5-7$ keV) band and its spectrum can be modeled by a single power-law model with photon index of $\Gamma = 1.81 \pm 0.18$ and full band flux of $f=3.95\times 10^{-15}$ erg s$^{-1}$ cm$^{-2}$. When compared with the data obtained by XMM-Newton in 2003, our Chandra observation in 2017 shows a harder ($\Delta \Gamma \approx -0.6$) spectrum and a 2.5 times fainter flux. Such a variation, in a timespan of $\sim2$ yrs rest-frame, is unexpected for such a luminous QSO powered by a $> 10^9 \: M_{\odot}$ black hole. The observed source hardening and weakening could be related to an intrinsic variation in the accretion rate. However, the limited photon statistics does not allow us to discriminate between an intrinsic luminosity and spectral change, and an absorption event produced by an intervening gas cloud along the line of sight. We also report the discovery of diffuse X-ray emission that extends for 30"x20" southward the QSO with a signal-to-noise ratio of $\sim$6, hardness ratio of $HR=0.03_{-0.25}^{+0.20}$, and soft band flux of $f_{0.5-2 \: keV}= 1.1_{-0.3}^{+0.3} \times 10^{-15}$ erg s$^{-1}$ cm$^{-2}$, that is not associated to a group or cluster of galaxies. We discuss two possible explanations for the extended emission, which may be either associated with the radio lobe of a nearby, foreground radio galaxy (at $z \approx 1-2$), or ascribed to the feedback from the QSO itself acting on its surrounding environment, as proposed by simulations of early black hole formation.

We use Multi Unit Spectroscopic Explorer (MUSE) observations of the galaxy cluster MACS J1149.5+2223 to explore the kinematics of the grand-design spiral galaxy Sp1149 hosting the SN Refsdal. Sp1149 lies at $z\simeq1.49$, has a stellar mass $M_*\simeq5\times10^9 \, \mathrm{M_\odot}$, a star-formation rate $\mathrm{SFR} \simeq1-6 \, \mathrm{M_\odot/yr}$ and represents a likely progenitor of a Milky-Way-like galaxy. All the four multiple images of Sp1149 in our data show strong OII-line emissions pointing to a clear rotation pattern. We take advantage of the gravitational lensing magnification effect ($\simeq 4 \times$) on the OII emission of the least distorted image to fit 3D kinematic models to the MUSE data-cube and derive the rotation curve and the velocity dispersion profile of Sp1149. We find that the rotation curve steeply rises, peaks at $R\simeq1$ kpc and then (initially) declines and flattens to an average $V_\mathrm{flat} = 128^{+29}_{-19}$ km/s. The shape of the rotation curve is well determined but the actual value of $V_\mathrm{flat}$ is quite uncertain because of the nearly face-on configuration of the galaxy. The intrinsic velocity dispersion due to gas turbulence is almost constant across the entire disc with an average of $27\pm5$ km/s. This value is consistent with $z=0$ measurements in the ionized gas component and a factor 2-4 lower than other estimates in different galaxies at similar redshifts. The average stellar-to-total mass fraction is of the order of one fifth. Our kinematic analysis returns the picture of a regular star-forming, mildly turbulent, rotation-dominated ($V / \sigma\simeq5$) spiral galaxy in a 4 Gyr old Universe.

We present a high resolution dissection of the two-dimensional total mass distribution in the core of the Hubble Frontier Fields galaxy cluster MACS J0416.1-2403, at z ~ 0.396. We exploit HST/WFC3 near-IR (F160W) imaging, VLT/MUSE spectroscopy, and Chandra data to separate the stellar, hot gas, and dark-matter mass components in the inner 300 kpc of the cluster. We combine the recent results of our refined strong lensing analysis, which includes the contribution of the intracluster gas, with the modeling of the surface brightness and stellar mass distributions of 193 cluster members, of which 144 are spectroscopically confirmed. We find that moving from 10 to 300 kpc from the cluster center the stellar to total mass fraction decreases from 12% to 1% and the hot gas to total mass fraction increases from 3% to 9%, resulting in a baryon fraction of approximately 10% at the outermost radius. We measure that the stellar component represents ~ 30%, near the cluster center, and 15%, at larger clustercentric distances, of the total mass in the cluster substructures. We subtract the baryonic mass component from the total mass distribution and conclude that within 30 kpc (~ 3 times the effective radius of the BCG) from the cluster center the surface mass density profile of the total mass and global (cluster plus substructures) dark-matter are steeper and that of the diffuse (cluster) dark-matter is shallower than a NFW profile. Our current analysis does not point to a significant offset between the cluster stellar and dark-matter components. This detailed and robust reconstruction of the inner dark-matter distribution in a larger sample of galaxy clusters will set a new benchmark for different structure formation scenarios.

We report the spectroscopic confirmation of 22 new multiply lensed sources behind the Hubble Frontier Field (HFF) galaxy cluster MACS~J0416.1$-$2403 (MACS 0416), using archival data from the Multi Unit Spectroscopic Explorer (MUSE) on the VLT. Combining with previous spectroscopic measurements of 15 other multiply imaged sources, we obtain a sample of 102 secure multiple images with measured redshifts, the largest to date in a single strong lensing system. The newly confirmed sources are largely low-luminosity Lyman-$\alpha$ emitters with redshift in the range [3.08-6.15]. With such a large number of secure constraints, and a significantly improved sample of galaxy members in the cluster core, we have improved our previous strong lensing model and obtained a robust determination of the projected total mass distribution of MACS 0416. We find evidence of three cored dark-matter halos, adding to the known complexity of this merging system. The total mass density profile, as well as the sub-halo population, are found in good agreement with previous works. We update and make public the redshift catalog of MACS 0416 from our previous spectroscopic campaign with the new MUSE redshifts. We also release lensing maps (convergence, shear, magnification) in the standard HFF format.

We present a new strong lensing analysis of the galaxy cluster MACS J1206.2-0847 (MACS 1206), at z=0.44, using deep spectroscopy from CLASH-VLT and VLT/MUSE archival data in combination with imaging from the Cluster Lensing and Supernova survey with Hubble. MUSE observations enable the spectroscopic identification of 23 new multiply imaged sources, extending the previous compilations by a factor of approximately five. In total, we use the positional measurements of 82 spectroscopic multiple images belonging to 27 families at z=1.0-6.1 to reconstruct the projected total mass distribution of MACS 1206. Remarkably, 11 multiple images are found within 50 kpc of the brightest cluster galaxy, making this an unprecedented set of constraints for the innermost projected mass distribution of a galaxy cluster. We thus find that, although dynamically relaxed, the smooth matter component (dark matter plus hot gas) of MACS 1206 shows a significant asymmetry, which closely follows the asymmetric distribution of the stellar component (galaxy members and intracluster light). We determine the value of the innermost logarithmic slope of the projected total mass density profile and find it to be close to the canonical Navarro-Frenk-White value. We demonstrate that this quantity is very robust against different parametrizations of the diffuse mass component; however, this is not the case when only one central image is used in the mass reconstruction. We also show that the mass density profile from our new strong lensing model is in very good agreement with dynamical and X-ray measurements at larger radii, where they overlap.

We present a novel approach for a combined analysis of X-ray and gravitational lensing data and apply this technique to the merging galaxy cluster MACS J0416.1$-$2403. The method exploits the information on the intracluster gas distribution that comes from a fit of the X-ray surface brightness, and then includes the hot gas as a fixed mass component in the strong lensing analysis. With our new technique, we can separate the collisional from the collision-less diffuse mass components, thus obtaining a more accurate reconstruction of the dark matter distribution in the core of a cluster. We introduce an analytical description of the X-ray emission coming from a set of dual Pseudo-Isothermal Elliptical (dPIE) mass distributions, which can be directly used in most lensing softwares. By combining \emph{Chandra} observations with Hubble Frontier Fields imaging and MUSE spectroscopy in MACS J0416.1$-$2403, we measure a projected gas over total mass fraction of approximately $10\%$ at $350$ kpc from the cluster center. Compared to the results of a more traditional cluster mass model (diffuse halos plus member galaxies), we find a significant difference in the cumulative projected mass profile of the dark matter component and that the dark matter to total mass fraction is almost constant, out to more than $350$ kpc. In the coming era of large surveys, these results show the need of multi-probe analyses for detailed dark matter studies in galaxy clusters.

We perform a maximum likelihood kinematic analysis of the two dynamically relaxed galaxy clusters MACS J1206.2-0847 at $z=0.44$ and RXC J2248.7-4431 at $z=0.35$ to determine the total mass profile in modified gravity models, using a modified version of the MAMPOSSt code of Mamon, Biviano and Bou\'e. Our work is based on the kinematic and lensing mass profiles derived using the data from the Cluster Lensing And Supernova survey with Hubble (hereafter CLASH) and the spectroscopic follow-up with the Very Large Telescope (hereafter CLASH-VLT). We assume a spherical Navarro-Frenk-White (NFW hereafter) profile in order to obtain a constraint on the fifth force interaction range $\lambda$ for models in which the dependence of this parameter on the enviroment is negligible at the scale considered (i.e. $\lambda=const$) and fixing the fifth force strength to the value predicted in $f(R)$ gravity. We then use information from lensing analysis to put a prior on the other NFW free parameters. In the case of MACSJ 1206 the joint kinematic+lensing analysis leads to an upper limit on the effective interaction range $\lambda\le1.61\,\mbox{Mpc}$ at $\Delta\chi^{2}=2.71$ on the marginalized distribution. For RXJ 2248 instead a possible tension with the $\Lambda$CDM model appears when adding lensing information, with a lower limit $\lambda\ge0.14\,\mbox{Mpc}$ at $\Delta\chi^{2}=2.71$. This is consequence of the slight difference between the lensing and kinematic data, appearing in GR for this cluster, that could be explained in terms of modifications of gravity. We discuss the impact of systematics and the limits of our analysis as well as future improvements of the results obtained.[...]

A massive galaxy cluster can serve as a magnifying glass for distant stellar populations, with strong gravitational lensing exposing details in the lensed background galaxies that would otherwise be undetectable. The MACS J0416.1-2403 cluster (hereafter MACS0416) is one of the most efficient lenses in the sky, and in 2014 it was observed with high-cadence imaging from the Hubble Space Telescope (HST). Here we describe two unusual transient events that appeared behind MACS0416 in a strongly lensed galaxy at redshift $z = 1.0054 \pm 0.0002$. These transients---designated HFF14Spo-NW and HFF14Spo-SE and collectively nicknamed "Spock"---were faster and fainter than any supernova (SN), but significantly more luminous than a classical nova. They reached peak luminosities of $\sim10^{41}$ erg s$^{-1}$ ($M_{\rm AB} < -14$ mag) in 5 rest-frame days, then faded below detectability in roughly the same time span. Models of the cluster lens suggest that these events may be spatially coincident at the source plane, but are most likely not temporally coincident. We find that HFF14Spo can be explained as a luminous blue variable (LBV), a recurrent nova (RN), or a pair of stellar microlensing events. To distinguish between these hypotheses will require a clarification of the positions of nearby critical curves, along with high-cadence monitoring of the field that could detect new transient episodes in the host galaxy.

Exploiting the powerful tool of strong gravitational lensing by galaxy clusters to study the highest-redshift Universe and cluster mass distributions relies on precise lens mass modelling. In this work, we present the first attempt at modelling line-of-sight mass distribution in addition to that of the cluster, extending previous modelling techniques that assume mass distributions to be on a single lens plane. We focus on the Hubble Frontier Field cluster MACS J0416.1-2403, and our multi-plane model reproduces the observed image positions with a rms offset of ~0.53". Starting from this best-fitting model, we simulate a mock cluster that resembles MACS J0416.1-2403 in order to explore the effects of line-of-sight structures on cluster mass modelling. By systematically analysing the mock cluster under different model assumptions, we find that neglecting the lensing environment has a significant impact on the reconstruction of image positions (rms ~0.3"); accounting for line-of-sight galaxies as if they were at the cluster redshift can partially reduce this offset. Moreover, foreground galaxies are more important to include into the model than the background ones. While the magnification factors of the lensed multiple images are recovered within ~10% for ~95% of them, those ~5% that lie near critical curves can be significantly affected by the exclusion of the lensing environment in the models (up to a factor of ~200). In addition, line-of-sight galaxies cannot explain the apparent discrepancy in the properties of massive subhalos between MACS J0416.1-2403 and N-body simulated clusters. Since our model of MACS J0416.1-2403 with line-of-sight galaxies only reduced modestly the rms offset in the image positions, we conclude that additional complexities, such as more flexible halo shapes, would be needed in future models of MACS J0416.1-2403.

We study the spectrophotometric properties of a highly magnified (\mu~40-70) pair of stellar systems identified at z=3.2222 behind the Hubble Frontier Field galaxy cluster MACS~J0416. Five multiple images (out of six) have been spectroscopically confirmed by means of VLT/MUSE and VLT/X-Shooter observations. Each image includes two faint (m_uv~30.6), young (<100 Myr), low-mass (<10^7 Msun), low-metallicity (12+Log(O/H)~7.7, or 1/10 solar) and compact (30 pc effective radius) stellar systems separated by ~300pc, after correcting for lensing amplification. We measured several rest-frame ultraviolet and optical narrow (\sigma_v <~ 25 km/s) high-ionization lines. These features may be the signature of very hot (T>50000 K) stars within dense stellar clusters, whose dynamical mass is likely dominated by the stellar component. Remarkably, the ultraviolet metal lines are not accompanied by Lya emission (e.g., CIV / Lya > 15), despite the fact that the Lya line flux is expected to be 150 times brighter (inferred from the Hbeta flux). A spatially-offset, strongly-magnified (\mu>50) Lya emission with a spatial extent <~7.6 kpc^2 is instead identified 2 kpc away from the system. The origin of such a faint emission can be the result of fluorescent Lya induced by a transverse leakage of ionizing radiation emerging from the stellar systems and/or can be associated to an underlying and barely detected object (with m_uv > 34 de-lensed). This is the first confirmed metal-line emitter at such low-luminosity and redshift without Lya emission, suggesting that, at least in some cases, a non-uniform covering factor of the neutral gas might hamper the Lya detection.

We present the VISTA-CFHT Stripe 82 (VICS82) survey: a near-infrared (J+Ks) survey covering 150 square degrees of the Sloan Digital Sky Survey (SDSS) equatorial Stripe 82 to an average depth of J=21.9 AB mag and Ks=21.4 AB mag (80% completeness limits; 5-sigma point source depths are approximately 0.5 mag brighter). VICS82 contributes to the growing legacy of multi-wavelength data in the Stripe 82 footprint. The addition of near-infrared photometry to the existing SDSS Stripe 82 coadd ugriz photometry reduces the scatter in stellar mass estimates to delta log(M_stellar)~0.3 dex for galaxies with M_stellar>10^9M_sun at z~0.5, and offers improvement compared to optical-only estimates out to z~1, with stellar masses constrained within a factor of approximately 2.5. When combined with other multi-wavelength imaging of the Stripe, including moderate-to-deep ultraviolet (GALEX), optical and mid-infrared (Spitzer IRAC) coverage, as well as tens of thousands of spectroscopic redshifts, VICS82 gives access to approximately 0.5 Gpc^3 of comoving volume. Some of the main science drivers of VICS82 include (a) measuring the stellar mass function of L^star galaxies out to z~1; (b) detecting intermediate redshift quasars at 2<z<3.5; (c) measuring the stellar mass function and baryon census of clusters of galaxies, and (d) performing optical/near-infrared-cosmic microwave background lensing cross-correlation experiments linking stellar mass to large-scale dark matter structure. Here we define and describe the survey, highlight some early science results and present the first public data release, which includes an SDSS-matched catalogue as well as the calibrated pixel data itself.

We set out to quantify the number density of quiescent massive compact galaxies at intermediate redshifts. We determine structural parameters based on i-band imaging using the CFHT equatorial SDSS Stripe 82 (CS82) survey (~170 sq. degrees) taking advantage of an exquisite median seeing of ~0.6''. We select compact massive (M > 5x10^10 M_sun) galaxies within the redshift range of 0.2<z<0.6. The large volume sampled allows to decrease the effect of cosmic variance that has hampered the calculation of the number density for this enigmatic population in many previous studies. We undertake an exhaustive analysis in an effort to untangle the various findings inherent to the diverse definition of compactness present in the literature. We find that the absolute number of compact galaxies is very dependent on the adopted definition and can change up to a factor of >10. We systematically measure a factor of ~5 more compacts at the same redshift than what was previously reported on smaller fields with HST imaging, which are more affected by cosmic variance. This means that the decrease in number density from z ~ 1.5 to z ~ 0.2 might be only of a factor of ~2-5, significantly smaller than what previously reported. This supports progenitor bias as the main contributor to the size evolution. This milder decrease is roughly compatible with the predictions from recent numerical simulations. Only the most extreme compact galaxies, with Reff < 1.5x( M/10^11 M_sun)^0.75 and M > 10^10.7 M_sun, appear to drop in number by a factor of ~20 and hence likely experience a noticeable size evolution.

Star forming galaxies have long been considered the dominant sources of the cosmic ultraviolet background radiation at early epochs. However, observing and characterizing the galaxy population with significant ionizing emission has proven to be challenging. In particular, the fraction of ionizing radiation that escapes the local environment to the intergalactic medium is poorly known. We investigate the relation between the escape fraction and galaxy luminosity. We combine the deep ultraviolet observations of Hubble Ultra Deep Field (UVUDF) with the deep Multi Unit Spectroscopic Explorer (MUSE) observations of the same field, collecting a sample of 165 faint star forming galaxies in the $3 < z < 4$ redshift range with deep rest-frame observations of the Lyman continuum. In our sample, we do not find any galaxy with significant emission of LyC radiation. We bin the galaxies in various redshift and brightness intervals and stack their images. From stacked images we estimate the relative escape fraction upper limits as a function of the luminosity. Thanks to the depth of the sample we measure meaningful 1$\sigma$ upper limits of $f_{esc,rel} < 0.07, 0.2$ and 0.6 at $L \sim L_{\rm z=3}^{*}, 0.5L_{\rm z=3}^{*}$ and $0.1L_{\rm z=3}^{*}$, respectively. We use our estimates and theoretical predictions from the literature to study a possible dependence of the escape fraction on galaxy luminosity by modelling the ionizing background with different prescriptions of $f_{\rm esc} (M_{\rm UV})$. We show that the understanding of the luminosity dependence hinges on the ability to constrain the escape fraction down to $M_{\rm UV} \sim -18$ mag in the future.

We present a detailed strong lensing (SL) mass reconstruction of the core of the galaxy cluster MACSJ 2129.4-0741 ($\rm z_{cl}=0.589$) obtained by combining high-resolution HST photometry from the CLASH survey with new spectroscopic observations from the CLASH-VLT survey. A background bright red passive galaxy at $\rm z_{sp}=1.36$, sextuply lensed in the cluster core, has four radial lensed images located over the three central cluster members. Further 19 background lensed galaxies are spectroscopically confirmed by our VLT survey, including 3 additional multiple systems. A total of 31 multiple images are used in the lensing analysis. This allows us to trace with high precision the total mass profile of the cluster in its very inner region ($\rm R<100$ kpc). Our final lensing mass model reproduces the multiple images systems identified in the cluster core with high accuracy of $0.4"$. This translates to an high precision mass reconstruction of MACS 2129, which is constrained at level of 2%. The cluster has Einstein parameter $\Theta_E=(29\pm4)"$, and a projected total mass of $\rm M_{tot}(<\Theta_E)=(1.35\pm0.03)\times 10^{14}M_{\odot}$ within such radius. Together with the cluster mass profile, we provide here also the complete spectroscopic dataset for the cluster members and lensed images measured with VLT/VIMOS within the CLASH-VLT survey.

We report on five compact, extremely young (<10Myr) and blue (\beta_UV<-2.5, F_\lambda =\lambda^\beta) objects observed with VLT/MUSE at redshift 3.1169, 3.235, in addition to three objects at z=6.145. These sources are magnified by the Hubble Frontier Field galaxy clusters MACS~J0416 and AS1063. Their de-lensed half light radii (Re) are between 16 to 140pc, the stellar masses are ~1-20 X 10^6 Msun, the magnitudes are m_uv=28.8 - 31.4 (-17<Muv<-15) and specific star formation rates can be as large as ~800Gyr^-1. Multiple images of these systems are widely separated in the sky (up to 50'') and individually magnified by factors 3-40. Remarkably, the inferred physical properties of two objects are similar to those expected in some globular cluster formation scenarios, representing the best candidate proto-globular clusters (proto-GC) discovered so far. Rest-frame optical high dispersion spectroscopy of one of them at z=3.1169 yields a velocity dispersion \sigma_v~20km/s, implying a dynamical mass dominated by the stellar mass. Another object at z=6.145, with de-lensed Muv ~ -15.3 (m_uv ~ 31.4), shows a stellar mass and a star-formation rate surface density consistent with the values expected from popular GC formation scenarios. An additional star-forming region at z=6.145, with de-lensed m_uv ~ 32, a stellar mass of 0.5 X 10^6 Msun and a star formation rate of 0.06 Msun/yr is also identified. These objects currently represent the faintest spectroscopically confirmed star-forming systems at z>3, elusive even in the deepest blank fields. We discuss how proto-GCs might contribute to the ionization budget of the universe and augment Lya visibility during reionization. This work underlines the crucial role of JWST in characterizing the rest-frame optical and near-infrared properties of such low-luminosity high-z objects.

We report the discovery of a multiply lensed Lyman-$\alpha$ blob (LAB) behind the galaxy cluster AS1063 using the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT). The background source is at $z=$ 3.117 and is intrinsically faint compared to almost all previously reported LABs. We used our highly precise strong lensing model to reconstruct the source properties, and we find an intrinsic luminosity of $L_{\rm Ly\alpha}$=$1.9\times10^{42}$ erg s$^{-1}$, extending to 33 kpc. We find that the LAB is associated with a group of galaxies, and possibly a protocluster, in agreement with previous studies that find LABs in overdensities. In addition to Lyman-$\alpha$ (Ly$\alpha$) emission, we find \ion{C}{IV}, \ion{He}{II}, and \ion{O}{III}] ultraviolet (UV) emission lines arising from the centre of the nebula. We used the compactness of these lines in combination with the line ratios to conclude that the \Lya nebula is likely powered by embedded star formation. Resonant scattering of the \Lya photons then produces the extended shape of the emission. Thanks to the combined power of MUSE and strong gravitational lensing, we are now able to probe the circumgalatic medium of sub-$L_{*}$ galaxies at $z\approx 3$.

In spite of their conjectured importance for the Epoch of Reionization, the properties of low-mass galaxies are currently still under large debate. In this article, we study the stellar and gaseous properties of faint, low-mass galaxies at z>3. We observed the Frontier Fields cluster Abell S1063 with MUSE over a 2 arcmin^2 field, and combined integral-field spectroscopy with gravitational lensing to perform a blind search for intrinsically faint Lya emitters (LAEs). We determined in total the redshift of 172 galaxies of which 14 are lensed LAEs at z=3-6.1. We increased the number of spectroscopically-confirmed multiple-image families from 6 to 17 and updated our gravitational-lensing model accordingly. The lensing-corrected Lya luminosities are with L(Lya) <= 10^41.5 erg/s among the lowest for spectroscopically confirmed LAEs at any redshift. We used expanding gaseous shell models to fit the Lya line profile, and find low column densities and expansion velocities. This is to our knowledge the first time that gaseous properties of such faint galaxies at z>=3 are reported. We performed SED modelling to broadband photometry from the U-band through the infrared to determine the stellar properties of these LAEs. The stellar masses are very low (10^{6-8} Msun), and are accompanied by very young ages of 1-100 Myr. The very high specific star formation rates (~100 Gyr^-1) are characteristic of starburst galaxies. The UV-continuum slopes beta are low in our sample, with beta<-2 for all galaxies with M_star<10^8 Msun. We conclude that our low-mass galaxies at 3<z<6 are forming stars at higher rates when correcting for stellar mass effects than seen locally or in more massive galaxies. The young stellar populations with high star-formation rates and low HI column densities lead to continuum slopes and LyC-escape fractions expected for a scenario where low mass galaxies reionise the Universe.

This overview article describes the legacy prospect and discovery potential of the Dark Energy Survey (DES) beyond cosmological studies, illustrating it with examples from the DES early data. DES is using a wide-field camera (DECam) on the 4m Blanco Telescope in Chile to image 5000 sq deg of the sky in five filters (grizY). By its completion the survey is expected to have generated a catalogue of 300 million galaxies with photometric redshifts and 100 million stars. In addition, a time-domain survey search over 27 sq deg is expected to yield a sample of thousands of Type Ia supernovae and other transients. The main goals of DES are to characterise dark energy and dark matter, and to test alternative models of gravity; these goals will be pursued by studying large scale structure, cluster counts, weak gravitational lensing and Type Ia supernovae. However, DES also provides a rich data set which allows us to study many other aspects of astrophysics. In this paper we focus on additional science with DES, emphasizing areas where the survey makes a difference with respect to other current surveys. The paper illustrates, using early data (from `Science Verification', and from the first, second and third seasons of observations), what DES can tell us about the solar system, the Milky Way, galaxy evolution, quasars, and other topics. In addition, we show that if the cosmological model is assumed to be Lambda+ Cold Dark Matter (LCDM) then important astrophysics can be deduced from the primary DES probes. Highlights from DES early data include the discovery of 34 Trans Neptunian Objects, 17 dwarf satellites of the Milky Way, one published z > 6 quasar (and more confirmed) and two published superluminous supernovae (and more confirmed).

We report on the discovery of extended Ly-alpha nebulae at z~3.3 in the Hubble Ultra Deep Field (HUDF, ~ 40 kpc X 80 kpc) and behind the Hubble Frontier Fields galaxy cluster MACSJ0416 (~ 40kpc), spatially associated with groups of star-forming galaxies. VLT/MUSE integral field spectroscopy reveals a complex structure with a spatially-varying double peaked Ly-alpha emission. Overall, the spectral profiles of the two Ly-alpha nebulae are remarkably similar, both showing a prominent blue emission, more intense and slightly broader than the red peak. From the first nebula, located in the HUDF, no X-ray emission has been detected, disfavoring the possible presence of AGNs. Spectroscopic redshifts have been derived for 11 galaxies within two arcsec from the nebula and spanning the redshift range 1.037<z<5.97. The second nebula, behind MACSJ0416, shows three aligned star-forming galaxies plausibly associated to the emitting gas. In both systems, the associated galaxies reveal possible intense rest-frame-optical nebular emissions lines [OIII]4959-5007+Hbeta with equivalent widths as high as 1500A rest-frame and star formation rates ranging from a few to tens of solar masses per year. A possible scenario is that of a group of young, star-forming galaxies sources of escaping ionising radiation that induce Ly-alpha fluorescence, therefore revealing the kinematics of the surrounding gas. Also Ly-alpha powered by star-formation and/or cooling radiation may resemble the double peaked spectral properties and the morphology observed here. If the intense blue emission is associated with inflowing gas, then we may be witnessing an early phase of galaxy or a proto-cluster (or group) formation.

We present VLT/X-Shooter and MUSE spectroscopy of an faint F814W=28.60+/-0.33 (Muv=-17.0), low mass (~<10^7 Msun) and compact (Reff=62pc) freshly star-forming galaxy at z=3.1169 magnified (16x) by the Hubble Frontier Fields galaxy cluster Abell S1063. Gravitational lensing allows for a significant jump toward low-luminosity regimes, in moderately high resolution spectroscopy (R=lambda/dlambda ~ 3000-7400). We measured CIV1548,1550, HeII1640, OIII]1661,1666, CIII]1907,1909, Hbeta, [OIII]4959,5007, emission lines with FWHM< 50 km/s and (de-lensed) fluxes spanning the interval 1.0x10^-19 - 2.0x10^-18 erg/s/cm2 at S/N=4-30. The double peaked Lya emission with Delta_v(red-blue) = 280(+/-7)km/s and de-lensed fluxes 2.4_(blue)|8.5_(red)x10^-18 erg/s/cm2 (S/N=38_(blue)|110_(red)) indicate a low column density of neutral hydrogen gas consistent with a highly ionized interstellar medium as also inferred from the large [OIII]5007/[OII]3727>10 ratio. We detect CIV1548,1550 resonant doublet in emission, each component with FWHM ~< 45 km/s, and redshifted by +51(+/-10)km/s relative to the systemic redshift. We interpret this as nebular emission tracing an expanding optically-thin interstellar medium. Both CIV1548,1550 and HeII1640 suggest the presence of hot and massive stars (with a possible faint AGN). The ultraviolet slope is remarkably blue, beta =-2.95 +/- 0.20 (F_lambda=lambda^beta), consistent with a dust-free and young ~<20 Myr galaxy. Line ratios suggest an oxygen abundance 12+log(O/H)<7.8. We are witnessing an early episode of star-formation in which a relatively low NHI and negligible dust attenuation might favor a leakage of ionizing radiation. This galaxy currently represents a unique low-luminosity reference object for future studies of the reionization epoch with JWST.

We present VIMOS-VLT spectroscopy of the Frontier Fields cluster MACS~J0416.1-2403 (z=0.397). Taken as part of the CLASH-VLT survey, the large spectroscopic campaign provided more than 4000 reliable redshifts, including ~800 cluster member galaxies. The unprecedented sample of cluster members at this redshift allows us to perform a highly detailed dynamical and structural analysis of the cluster out to ~2.2$r_{200}$ (~4Mpc). Our analysis of substructures reveals a complex system composed of a main massive cluster ($M_{200}$~0.9$\times 10^{15} M_{\odot}$) presenting two major features: i) a bimodal velocity distribution, showing two central peaks separated by $\Delta V_{rf}$~1100 km s$^{-1}$ with comparable galaxy content and velocity dispersion, ii) a projected elongation of the main substructures along the NE-SW direction, with a prominent subclump ~600 kpc SW of the center and an isolated BCG approximately halfway between the center and the SW clump. We also detect a low mass structure at z~0.390, ~10' S of the cluster center, projected at ~3Mpc, with a relative line-of-sight velocity of $\Delta V_{rf}$~-1700 km s$^{-1}$. The cluster mass profile that we obtain through our dynamical analysis deviates significantly from the "universal" NFW, being best fit by a Softened Isothermal Sphere model instead. The mass profile measured from the galaxy dynamics is found to be in relatively good agreement with those obtained from strong and weak lensing, as well as with that from the X-rays, despite the clearly unrelaxed nature of the cluster. Our results reveal overall a complex dynamical state of this massive cluster and support the hypothesis that the two main subclusters are being observed in a pre-collisional phase, in line with recent findings from radio and deep X-ray data. With this article we also release the entire redshift catalog of 4386 sources in the field of this cluster.

We present MUSE observations in the core of the HFF galaxy cluster MACS J1149.5+2223, where the first magnified and spatially-resolved multiple images of SN 'Refsdal' at redshift 1.489 were detected. Thanks to a DDT program with the VLT and the extraordinary efficiency of MUSE, we measure 117 secure redshifts with just 4.8 hours of total integration time on a single target pointing. We spectroscopically confirm 68 galaxy cluster members, with redshift values ranging from 0.5272 to 0.5660, and 18 multiple images belonging to 7 background, lensed sources distributed in redshifts between 1.240 and 3.703. Starting from the combination of our catalog with those obtained from extensive spectroscopic and photometric campaigns using the HST, we select a sample of 300 (164 spectroscopic and 136 photometric) cluster members, within approximately 500 kpc from the BCG, and a set of 88 reliable multiple images associated to 10 different background source galaxies and 18 distinct knots in the spiral galaxy hosting SN 'Refsdal'. We exploit this valuable information to build 6 detailed strong lensing models, the best of which reproduces the observed positions of the multiple images with a rms offset of only 0.26". We use these models to quantify the statistical and systematic errors on the predicted values of magnification and time delay of the next emerging image of SN 'Refsdal'. We find that its peak luminosity should should occur between March and June 2016, and should be approximately 20% fainter than the dimmest (S4) of the previously detected images but above the detection limit of the planned HST/WFC3 follow-up. We present our two-dimensional reconstruction of the cluster mass density distribution and of the SN 'Refsdal' host galaxy surface brightness distribution. We outline the roadmap towards even better strong lensing models with a synergetic MUSE and HST effort.

We use high-precision kinematic and lensing measurements of the total mass profile of the dynamically relaxed galaxy cluster MACS J1206.2-0847 at $z=0.44$ to estimate the value of the ratio $\eta=\Psi/\Phi$ between the two scalar potentials in the linear perturbed Friedmann-Lemaitre-Robertson-Walker metric.[...] Complementary kinematic and lensing mass profiles were derived from exhaustive analyses using the data from the Cluster Lensing And Supernova survey with Hubble (CLASH) and the spectroscopic follow-up with the Very Large Telescope (CLASH-VLT). Whereas the kinematic mass profile tracks only the time-time part of the perturbed metric (i.e. only $\Phi$), the lensing mass profile reflects the contribution of both time-time and space-space components (i.e. the sum $\Phi+\Psi$). We thus express $\eta$ as a function of the mass profiles and perform our analysis over the radial range $0.5\,Mpc\le r\le r_{200}=1.96\,Mpc$. Using a spherical Navarro-Frenk-White mass profile, which well fits the data, we obtain $\eta(r_{200})=1.01\,_{-0.28}^{+0.31}$ at the 68\% C.L. We discuss the effect of assuming different functional forms for mass profiles and of the orbit anisotropy in the kinematic reconstruction. Interpreting this result within the well-studied $f(R)$ modified gravity model, the constraint on $\eta$ translates into an upper bound to the interaction length (inverse of the scalaron mass) smaller than 2 Mpc. This tight constraint on the $f(R)$ interaction range is however substantially relaxed when systematic uncertainties in the analysis are considered. Our analysis highlights the potential of this method to detect deviations from general relativity, while calling for the need of further high-quality data on the total mass distribution of clusters and improved control on systematic effects.

We present a strong lensing system, composed of 4 multiple images of a source at z = 2.387, created by two lens galaxies, G1 and G2, belonging to the galaxy cluster MACS J1115.9+0129 at z = 0.353. We use observations taken as part of the Cluster Lensing and Supernova survey with Hubble, CLASH, and its spectroscopic follow-up programme at the Very Large Telescope, CLASH-VLT, to estimate the total mass distributions of the two galaxies and the cluster through strong gravitational lensing models. We find that the total projected mass values within the half-light radii, R_{e}, of the two lens galaxies are M_{T,G1}(< R_{e,G1}) = (3.6 +/- 0.4) x 10^{10}M_{Sun} and M_{T,G2}(< R_{e,G2}) = (4.2 +/- 1.6) x 10^{10}M_{Sun}. The effective velocity dispersion values of G1 and G2 are (122 +/- 7) km/s and (137 +/- 27) km/s, respectively. We remark that these values are relatively low when compared to those of ~200-300 km/s, typical of lens galaxies found in the field by previous surveys. By fitting the spectral energy distributions of G1 and G2, we measure projected luminous over total mass fractions within R_{e} of 0.11 +/- 0.03, for G1, and 0.73 +/- 0.32, for G2. The fact that the less massive galaxy, G1, is dark-matter dominated in its inner regions raises the question of whether the dark matter fraction in the core of early-type galaxies depends on their mass. Further investigating strong lensing systems will help us understand the influence that dark matter has on the structure and evolution of the inner regions of galaxies.

We perform a comprehensive study of the total mass distribution of the galaxy cluster RXCJ2248 ($z=0.348$) with a set of high-precision strong lensing models, which take advantage of extensive spectroscopic information on many multiply lensed systems. In the effort to understand and quantify inherent systematics in parametric strong lensing modelling, we explore a collection of 22 models where we use different samples of multiple image families, parametrizations of the mass distribution and cosmological parameters. As input information for the strong lensing models, we use the CLASH HST imaging data and spectroscopic follow-up observations, carried out with the VIMOS and MUSE spectrographs, to identify bona-fide multiple images. A total of 16 background sources, over the redshift range $1.0-6.1$, are multiply lensed into 47 images, 24 of which are spectroscopically confirmed and belong to 10 individual sources. The cluster total mass distribution and underlying cosmology in the models are optimized by matching the observed positions of the multiple images on the lens plane. We show that with a careful selection of a sample of spectroscopically confirmed multiple images, the best-fit model reproduces their observed positions with a rms of $0.3$ in a fixed flat $\Lambda$CDM cosmology, whereas the lack of spectroscopic information lead to biases in the values of the model parameters. Allowing cosmological parameters to vary together with the cluster parameters, we find (at $68\%$ confidence level) $\Omega_m=0.25^{+0.13}_{-0.16}$ and $w=-1.07^{+0.16}_{-0.42}$ for a flat $\Lambda$CDM model, and $\Omega_m=0.31^{+0.12}_{-0.13}$ and $\Omega_\Lambda=0.38^{+0.38}_{-0.27}$ for a universe with $w=-1$ and free curvature. Using toy models mimicking the overall configuration of RXCJ2248, we estimate the impact of the line of sight mass structure on the positional rms to be $0.3\pm 0.1$.(ABRIDGED)

Supernova "Refsdal," multiply imaged by cluster MACSJ1149.5+2223, represents a rare opportunity to make a true blind test of model predictions in extragalactic astronomy, on a time scale that is short compared to a human lifetime. In order to take advantage of this event, we produced seven gravitational lens models with five independent methods, based on Hubble Space Telescope (HST) Hubble Frontier Field images, along with extensive spectroscopic follow-up observations by HST, the Very Large and the Keck Telescopes. We compare the model predictions and show that they agree reasonably well with the measured time delays and magnification ratios between the known images, even though these quantities were not used as input. This agreement is encouraging, considering that the models only provide statistical uncertainties, and do not include additional sources of uncertainties such as structure along the line of sight, cosmology, and the mass sheet degeneracy. We then present the model predictions for the other appearances of SN "Refsdal." A future image will reach its peak in the first half of 2016, while another image appeared between 1994 and 2004. The past image would have been too faint to be detected in existing archival images. The future image should be approximately one third as bright as the brightest known image (i.e., H_AB~25.7 mag at peak and H_AB~26.7 mag six months before peak), and thus detectable in single-orbit HST images. We will find out soon whether our predictions are correct.

We report the observation and confirmation of the first group- and cluster-scale strong gravitational lensing systems found in Dark Energy Survey (DES) data. Through visual inspection of data from the Science Verification (SV) season, we identified 53 candidate systems. We then obtained spectroscopic follow-up of 21 candidates using the Gemini Multi-Object Spectrograph (GMOS) at the Gemini South telescope and the Inamori-Magellan Areal Camera and Spectrograph (IMACS) at the Magellan/Baade telescope. With this follow-up, we confirmed six candidates as gravitational lenses: Three of the systems are newly discovered, and the remaining three were previously known. Of the 21 observed candidates, the remaining 15 were either not detected in spectroscopic observations, were observed and did not exhibit continuum emission (or spectral features), or were ruled out as lensing systems. The confirmed sample consists of one group-scale and five galaxy cluster-scale lenses. The lensed sources range in redshift z ~ 0.80-3.2, and in i-band surface brightness i_{SB} ~ 23-25 mag/sq.-arcsec. (2" aperture). For each of the six systems, we estimate the Einstein radius and the enclosed mass, which have ranges ~ 5.0 - 8.6" and ~ 7.5 x 10^{12} - 6.4 x 10^{13} solar masses, respectively.

The analysis of galaxy properties and the relations among them and the environment, can be used to investigate the physical processes driving galaxy evolution. We study the cluster A209 by using the CLASH-VLT spectroscopic data combined with Subaru photometry, yielding to 1916 cluster members down to a stellar mass of 10^{8.6} Msun. We determine: i) the stellar mass function of star-forming and passive galaxies; ii) the intra-cluster light and its properties; iii) the orbits of low- and high-mass passive galaxies; and iv) the mass-size relation of ETGs. The stellar mass function of the star-forming galaxies does not depend on the environment, while the slope found for passive galaxies becomes flatter in the densest region. The color distribution of the intra-cluster light is consistent with the color of passive members. The analysis of the dynamical orbits shows that low-mass passive galaxies have tangential orbits, avoiding small pericenters around the BCG. The mass-size relation of low-mass passive ETGs is flatter than that of high mass galaxies, and its slope is consistent with that of field star-forming galaxies. Low-mass galaxies are also more compact within the scale radius of 0.65 Mpc. The ratio between stellar and number density profiles shows a mass segregation in the center. The comparative analysis of the stellar and total density profiles indicates that this effect is due to dynamical friction. Our results are consistent with a scenario in which the "environmental quenching" of low-mass galaxies is due to mechanisms such as harassment out to R200, starvation and ram-pressure stripping at smaller radii, as supported by the analysis of the mass function, of the dynamical orbits and of the mass-size relation of passive early-types in different regions. Our analyses support the idea that the intra-cluster light is formed through the tidal disruption of subgiant galaxies.

In the effort to understand the link between the structure of galaxy clusters and their galaxy populations, we focus on MACSJ1206.2-0847 at z~0.44 and probe its substructure in the projected phase space through the spectrophotometric properties of a large number of galaxies from the CLASH-VLT survey. Our analysis is mainly based on an extensive spectroscopic dataset of 445 member galaxies, mostly acquired with VIMOS@VLT as part of our ESO Large Programme, sampling the cluster out to a radius ~2R200 (4 Mpc). We classify 412 galaxies as passive, with strong Hdelta absorption (red and blue galaxies, and with emission lines from weak to very strong. A number of tests for substructure detection are applied to analyze the galaxy distribution in the velocity space, in 2D space, and in 3D projected phase-space. Studied in its entirety, the cluster appears as a large-scale relaxed system with a few secondary, minor overdensities in 2D distribution. We detect no velocity gradients or evidence of deviations in local mean velocities. The main feature is the WNW-ESE elongation. The analysis of galaxy populations per spectral class highlights a more complex scenario. The passive galaxies and red strong Hdelta galaxies trace the cluster center and the WNW-ESE elongated structure. The red strong Hdelta galaxies also mark a secondary, dense peak ~2 Mpc at ESE. The emission line galaxies cluster in several loose structures, mostly outside R200. The observational scenario agrees with MACS J1206.2-0847 having WNW-ESE as the direction of the main cluster accretion, traced by passive galaxies and red strong Hdelta galaxies. The red strong Hdelta galaxies, interpreted as poststarburst galaxies, date a likely important event 1-2 Gyr before the epoch of observation. The emission line galaxies trace a secondary, ongoing infall where groups are accreted along several directions.

We present a detailed mass reconstruction and a novel study on the substructure properties in the core of the CLASH and Frontier Fields galaxy cluster MACS J0416.1-2403. We show and employ our extensive spectroscopic data set taken with the VIMOS instrument as part of our CLASH-VLT program, to confirm spectroscopically 10 strong lensing systems and to select a sample of 175 plausible cluster members to a limiting stellar mass of log(M_*/M_Sun) ~ 8.6. We reproduce the measured positions of 30 multiple images with a remarkable median offset of only 0.3" by means of a comprehensive strong lensing model comprised of 2 cluster dark-matter halos, represented by cored elliptical pseudo-isothermal mass distributions, and the cluster member components. The latter have total mass-to-light ratios increasing with the galaxy HST/WFC3 near-IR (F160W) luminosities. The measurement of the total enclosed mass within the Einstein radius is accurate to ~5%, including systematic uncertainties. We emphasize that the use of multiple-image systems with spectroscopic redshifts and knowledge of cluster membership based on extensive spectroscopic information is key to constructing robust high-resolution mass maps. We also produce magnification maps over the central area that is covered with HST observations. We investigate the galaxy contribution, both in terms of total and stellar mass, to the total mass budget of the cluster. When compared with the outcomes of cosmological $N$-body simulations, our results point to a lack of massive subhalos in the inner regions of simulated clusters with total masses similar to that of MACS J0416.1-2403. Our findings of the location and shape of the cluster dark-matter halo density profiles and on the cluster substructures provide intriguing tests of the assumed collisionless, cold nature of dark matter and of the role played by baryons in the process of structure formation.

The Navarro-Frenk-White (NFW) density profile is often used to model gravitational lenses. For low values of the characteristic convergence ($\kappa_s \ll 1$) of this model - corresponding to galaxy and galaxy group mass scales - a high numerical precision is required in order to accurately compute several quantities in the strong lensing regime. An alternative for fast and accurate computations is to derive analytic approximations in this limit. In this work we obtain analytic solutions for several lensing quantities for elliptical (ENFW) and pseudo-elliptical (PNFW) NFW lens models on the typical scales where gravitational arcs are expected to be formed, in the $\kappa_s \ll 1$ limit, establishing their domain of validity. We derive analytic solutions for the convergence and shear for these models, obtaining explicit expressions for the iso-convergence contours and constant distortion curves (including the tangential critical curve). We also compute the deformation cross section, which is given in closed form for the circular NFW model and in terms of a one-dimensional integral for the elliptical ones. In addition, we provide a simple expression for the ellipticity of the iso-convergence contours of the pseudo-elliptical models and the connection of characteristic convergences among the PNFW and ENFW models. We conclude that the set of solutions derived here is generally accurate for $\kappa_s \lesssim 0.1$. For low ellipticities, values up to $\kappa_s \simeq 0.18$ are allowed. On the other hand, the mapping between PNFW and the ENFW models is valid up to $\kappa_s \simeq 0.4$. The solutions derived in this work can be used to speed up numerical codes and ensure their accuracy in the low $\kappa_s$ regime, including applications to arc statistics and other strong lensing observables. (Abridged)

The statistics of gravitational arcs in galaxy clusters is a powerful probe of cluster structure and may provide complementary cosmological constraints. Despite recent progresses, discrepancies still remain among modelling and observations of arc abundance, specially regarding the redshift distribution of strong lensing clusters. Besides, fast "semi-analytic" methods still have to incorporate the success obtained with simulations. In this paper we discuss the contribution of the magnification in gravitational arc statistics. Although lensing conserves surface brightness, the magnification increases the signal-to-noise ratio of the arcs, enhancing their detectability. We present an approach to include this and other observational effects in semi-analytic calculations for arc statistics. The cross section for arc formation ({\sigma}) is computed through a semi-analytic method based on the ratio of the eigenvalues of the magnification tensor. Using this approach we obtained the scaling of {\sigma} with respect to the magnification, and other parameters, allowing for a fast computation of the cross section. We apply this method to evaluate the expected number of arcs per cluster using an elliptical Navarro--Frenk--White matter distribution. Our results show that the magnification has a strong effect on the arc abundance, enhancing the fraction of arcs, moving the peak of the arc fraction to higher redshifts, and softening its decrease at high redshifts. We argue that the effect of magnification should be included in arc statistics modelling and that it could help to reconcile arcs statistics predictions with the observational data.

The Perturbative Approach (PA) introduced by \citet{alard07} provides analytic solutions for gravitational arcs by solving the lens equation linearized around the Einstein ring solution. This is a powerful method for lens inversion and simulations in that it can be used, in principle, for generic lens models. In this paper we aim to quantify the domain of validity of this method for three quantities derived from the linearized mapping: caustics, critical curves, and the deformation cross section (i.e. the arc cross section in the infinitesimal circular source approximation). We consider lens models with elliptical potentials, in particular the Singular Isothermal Elliptic Potential and Pseudo-Elliptical Navarro--Frenk--White models. We show that the PA is exact for this first model. For the second, we obtain constraints on the model parameter space (given by the potential ellipticity parameter $\varepsilon$ and characteristic convergence $\kappa_s$) such that the PA is accurate for the aforementioned quantities. In this process we obtain analytic expressions for several lensing functions, which are valid for the PA in general. The determination of this domain of validity could have significant implications for the use of the PA, but it still needs to be probed with extended sources.

We present the first results of the SOAR (Southern Astrophysical Research) Gravitational Arc Survey (SOGRAS). The survey imaged 47 clusters in two redshift intervals centered at $z=0.27$ and $z=0.55$, targeting the richest clusters in each interval. Images were obtained in the $g'$, $r'$ and $i'$ bands using the SOAR Optical Imager (SOI), with a median seeing of 0.83, 0.76 and 0.71 arcsec, respectively, in these filters. Most of the survey clusters are located within the Sloan Digital Sky Survey (SDSS) Stripe 82 region and all of them are in the SDSS footprint. Photometric calibration was therefore performed using SDSS stars located in our SOI fields. We reached for galaxies in all fields the detection limits of $g \sim 23.5$, $r \sim 23$ and $i \sim 22.5$ for a signal-to-noise ratio (S/N) = 3. As a by-product of the image processing, we generated a source catalogue with 19760 entries, the vast majority of which are galaxies, where we list their positions, magnitudes and shape parameters. We compared our galaxy shape measurements to those of local galaxies and concluded that they were not strongly affected by seeing. From the catalogue data, we are able to identify a red sequence of galaxies in most clusters in the lower $z$ range. We found 16 gravitational arc candidates around 8 clusters in our sample. They tend to be bluer than the central galaxies in the lensing cluster. A preliminary analysis indicates that $\sim 10%$ of the clusters have arcs around them, with a possible indication of a larger efficiency associated to the high-$z$ systems when compared to the low-$z$ ones. Deeper follow-up images with Gemini strengthen the case for the strong lensing nature of the candidates found in this survey.

Models with elliptical potentials (pseudo-elliptical) are often used in gravitational lensing applications. Nevertheless, they generally lead to nonphysical mass distributions in some regions. In this paper we revisit the physical limitations of the pseudo-elliptical Navarro-Frenk-White (PNFW) model, for a broad range of the potential ellipticity parameter \epsilon\ and characteristic convergence \kappa_s focusing on the behavior of the mass distribution close to the tangential critical curve, where tangential arcs are expected to be formed. We investigate the shape of the mass distribution on this region and the presence of negative convergence. We obtain a mapping from the PNFW to the NFW model with elliptical mass distribution (ENFW) and provide fitting formulae for connecting the parameters of both models. We compare the arc cross section for these models using the "infinitesimal circular source approximation". We find that the PNFW model is well-suited to model an elliptical mass distribution on a larger \epsilon\ - \kappa_s parameter space than previously expected. In particular values as large as \epsilon\ ~ 0.65 are allowed for small \kappa_s. However, if we require the PNFW model to reproduce the arc cross section of the ENFW well, the ellipticity is more restricted. We also find that the negative convergence regions occur far from the arc formation region and should therefore not be a problem for studies with gravitational arcs. The determination of a domain of validity for the PNFW model and the mapping to ENFW models could have implications for the use of PNFW models for the inverse modeling of lenses and for fast arc simulations.