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

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

High-z galaxy gas-phase metallicities are usually determined through observations of strong optical emission lines with calibrations tied to the local universe. Recent debate has questioned if these calibrations are valid in the high-z universe. We investigate this by analysing a sample of 16 galaxies at z~2 available in the literature, and for which the metallicity can be robustly determined using oxygen auroral lines. The sample spans a redshift range of 1.4 < z < 3.6, has metallicities of 7.4-8.4 in 12+log(O/H) and stellar masses 10^7.5-10^11 Msun. We test commonly used strong line diagnostics (R23, O3, O2, O32, N2, O3N2 and Ne3O2 ) as prescribed by four different sets of empirical calibrations, as well as one fully theoretical calibration. We find that none of the strong line diagnostics (or calibration set) tested perform consistently better than the others. Amongst the line ratios tested, R23 and O3 deliver the best results, with accuracies as good as 0.01-0.04 dex and dispersions of ~0.2 dex in two of the calibrations tested. Generally, line ratios involving nitrogen predict higher values of metallicity, while results with O32 and Ne3O2 show large dispersions. The theoretical calibration yields an accuracy of 0.06 dex, comparable to the best strong line methods. We conclude that, within the metallicity range tested in this work, the locally calibrated diagnostics can still be reliably applied at z~2.

We present spatially resolved analysis of a lensed galaxy, SDSS1958+5950 at $z = 2.225$, from the Cambridge Sloan Survey of Wide Arcs in the Sky (CASSOWARY). We use our new high resolution imaging data to construct a robust lens model for the galaxy group at $z = 0.214$. We employ the updated lens model to combine the Integral Field Spectrographic observations on two highly distorted images of the lensed target. We adopt a forward-modeling approach to deconvolve the effects of point spread function from the combined source-plane reconstruction. The approach is adapted to the lens model magnification and enables a resolution of $\sim$170 pc in the galaxy-source plane. We propose an ongoing merger as the origin of the lensed system on the basis of its source-plane morphology, kinematics and rest-frame emission line ratios. Using our novel technique of adaptive coadded source plane reconstruction, we are able to detect different components in the velocity gradient that were not seen in previous studies of this object, plausibly belonging to different components in the merging system.

Galaxies at low-redshift typically possess negative gas-phase metallicity gradients (centres more metal-rich than their outskirts). Whereas, it is not uncommon to observe positive metallicity gradients in higher-redshift galaxies ($z \gtrsim 0.6$). Bridging these epochs, we present gas-phase metallicity gradients of 84 star-forming galaxies between $0.08 < z < 0.84$. Using the galaxies with reliably determined metallicity gradients, we measure the median metallicity gradient to be negative ($-0.039^{+0.007}_{-0.009}$ dex/kpc). Underlying this, however, is significant scatter: $(8\pm3)\%\ [7]$ of galaxies have significantly positive metallicity gradients, $(38 \pm 5)\%\ [32]$ have significantly negative gradients, $(31\pm5)\%\ [26]$ have gradients consistent with being flat. (The remaining $(23\pm5)\%\ [19]$ have unreliable gradient estimates.) We notice a slight trend for a more negative metallicity gradient with both increasing stellar mass and increasing star formation rate (SFR). However, given the potential redshift and size selection effects, we do not consider these trends to be significant. Indeed, once we normalize the SFR relative to that of the main sequence, we do not observe any trend between the metallicity gradient and the normalized SFR. This is contrary to recent studies of galaxies at similar and higher redshifts. We do, however, identify a novel trend between the metallicity gradient of a galaxy and its size. Small galaxies ($r_d < 3$ kpc) present a large spread in observed metallicity gradients (both negative and positive gradients). In contrast, we find no large galaxies ($r_d > 3$ kpc) with positive metallicity gradients, and overall there is less scatter in the metallicity gradient amongst the large galaxies. These large (well-evolved) galaxies may be analogues of present-day galaxies, which also show a common negative metallicity gradient.

The Lyman alpha (lya) line of Hydrogen is a prominent feature in the spectra of star-forming galaxies, usually redshifted by a few hundreds of km/s compared to the systemic redshift. This large offset hampers follow-up surveys, galaxy pair statistics and correlations with quasar absorption lines when only lya is available. We propose diagnostics that can be used to recover the systemic redshift directly from the properties of the lya line profile. We use spectroscopic observations of Lyman-Alpha Emitters (LAEs) for which a precise measurement of the systemic redshift is available. Our sample contains 13 sources detected between z~3 and z~6 as part of various Multi Unit Spectroscopic Explorer (MUSE) Guaranteed Time Observations (GTO). We also include a compilation of spectroscopic lya data from the literature spanning a wide redshift range (z~0-8). First, restricting our analysis to double-peaked lya spectra, we find a tight correlation between the velocity offset of the red peak with respect to the systemic redshift, Vpeak, and the separation of the peaks. Secondly, we find a correlation between Vpeak and the full width at half maximum of the lya line. Fitting formulas, to estimate systemic redshifts of galaxies with an accuracy of +-100 km/s when only the lya emission line is available, are given for the two methods.

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

We present a method to recover the gas-phase metallicity gradients from integral field spectroscopic (IFS) observations of barely resolved galaxies. We take a forward modelling approach and compare our models to the observed spatial distribution of emission line fluxes, accounting for the degrading effects of seeing and spatial binning. The method is flexible and is not limited to particular emission lines or instruments. We test the model through comparison to synthetic observations and use downgraded observations of nearby galaxies to validate this work. As a proof of concept we also apply the model to real IFS observations of high-redshift galaxies. From our testing we show that the inferred metallicity gradients and central metallicities are fairly insensitive to the assumptions made in the model and that they are reliably recovered for galaxies with sizes approximately equal to the half width at half maximum of the point-spread function. However, we also find that the presence of star forming clumps can significantly complicate the interpretation of metallicity gradients in moderately resolved high-redshift galaxies. Therefore we emphasize that care should be taken when comparing nearby well-resolved observations to high-redshift observations of partially resolved galaxies.

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

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

We present the results obtained with MUSE on the core of the lensing cluster A1689. Integral-field observations with MUSE provide a unique view of the central region, allowing us to conduct a complete census on both cluster galaxies and lensed background sources, identified based on their spectral features without preselection. We investigate the multiple-image configuration for all known sources in the field. Previous to our survey, 28 different lensed galaxies displaying 46 multiple images were known in the MUSE field of view, most of them based on photometric redshifts and lensing considerations. Among them, we spectroscopically confirm 12 images based on their emission-lines, corresponding to 7 different lensed galaxies between z = 0.95 and 5.0. In addition, 14 new galaxies have been spectroscopically identified in this area, with redshifts ranging between 0.8 and 6.2. All background sources within the MUSE field of view correspond to multiple-imaged systems lensed by A1689. 17 sources in total are found at z > 3 based on their Lyman-alpha emission, with Lyman-alpha luminosities ranging between 40.5 < log(Ly{\alpha}) < 42.5 after correction for magnification. This sample is particularly sensitive to the slope of the LF toward the faintest-end. The density of sources obtained in this survey is consistent with a steep value of {\alpha} < -1.5, although this result still needs further investigation. These results illustrate the efficiency of MUSE in the characterization of lensing clusters on one hand, and the study of faint and distant populations of galaxies on the other hand. In particular, our current survey of lensing clusters should provide a unique census of sources responsible for the reionization in a representative volume at z ~ 4-7.

We report the detection of extended Ly alpha emission around individual star-forming galaxies at redshifts z = 3-6 in an ultradeep exposure of the Hubble Deep Field South obtained with MUSE on the ESO-VLT. The data reach a limiting surface brightness (1sigma) of ~1 x 10^-19 erg s^-1 cm^-2 arcsec^-2 in azimuthally averaged radial profiles, an order of magnitude improvement over previous narrowband imaging. Our sample consists of 26 spectroscopically confirmed Ly alpha-emitting, but mostly continuum-faint (m_AB >~ 27) galaxies. In most objects the Ly alpha emission is considerably more extended than the UV continuum light. While 5 of the faintest galaxies in the sample show no significantly detected Ly alpha haloes, the derived upper limits suggest that this is just due to insufficient S/N. Ly alpha haloes therefore appear to be (nearly) ubiquitous even for low-mass (~10^8-10^9 M_sun) star-forming galaxies at z>3. We decompose the Ly alpha emission of each object into a compact `continuum-like' and an extended halo component, and infer sizes and luminosities of the haloes. The extended Ly alpha emission approximately follows an exponential surface brightness distribution with a scale length of a few kpc. While these haloes are thus quite modest in terms of their absolute sizes, they are larger by a factor of 5-15 than the corresponding rest-frame UV continuum sources as seen by HST. They are also much more extended, by a factor ~5, than Ly alpha haloes around low-redshift star-forming galaxies. Between ~40% and >90% of the observed Ly alpha flux comes from the extended halo component, with no obvious correlation of this fraction with either the absolute or the relative size of the Ly alpha halo. Our observations provide direct insights into the spatial distribution of at least partly neutral gas residing in the circumgalactic medium of low to intermediate mass galaxies at z > 3.

Spatially resolved studies of high redshift galaxies, an essential insight into galaxy formation processes, have been mostly limited to stacking or unusually bright objects. We present here the study of a typical (L$^{*}$, M$_\star$ = 6 $\times 10^9$ $M_\odot$) young lensed galaxy at $z=3.5$, observed with MUSE, for which we obtain 2D resolved spatial information of Ly$\alpha$ and, for the first time, of CIII] emission. The exceptional signal-to-noise of the data reveals UV emission and absorption lines rarely seen at these redshifts, allowing us to derive important physical properties (T$_e\sim$15600 K, n$_e\sim$300 cm$^{-3}$, covering fraction f$_c\sim0.4$) using multiple diagnostics. Inferred stellar and gas-phase metallicities point towards a low metallicity object (Z$_{\mathrm{stellar}}$ = $\sim$ 0.07 Z$_\odot$ and Z$_{\mathrm{ISM}}$ $<$ 0.16 Z$_\odot$). The Ly$\alpha$ emission extends over $\sim$10 kpc across the galaxy and presents a very uniform spectral profile, showing only a small velocity shift which is unrelated to the intrinsic kinematics of the nebular emission. The Ly$\alpha$ extension is $\sim$4 times larger than the continuum emission, and makes this object comparable to low-mass LAEs at low redshift, and more compact than the Lyman-break galaxies and Ly$\alpha$ emitters usually studied at high redshift. We model the Ly$\alpha$ line and surface brightness profile using a radiative transfer code in an expanding gas shell, finding that this model provides a good description of both observables.