We present near-infrared observations of 42 gravitationally lensed galaxies obtained in the framework of the KMOS Lensed Emission Lines and VElocity Review (KLEVER) Survey, a program aimed at investigating the spatially resolved properties of the ionised gas in 1.2<z<2.5 galaxies by means of a full coverage of the YJ, H and K near-infrared bands. Detailed metallicity maps and gradients are derived for a sub-sample of 28 galaxies from reconstructed source plane emission line maps, exploiting the variety of different emission line diagnostics provided by the broad wavelength coverage of the survey. About 85% of these galaxies are characterised by metallicity gradients shallower than 0.05 dex/kpc and 89% are consistent with flat slope within 3$\sigma$ (67% within 1$\sigma$), suggesting a mild evolution with cosmic time. In the context of cosmological simulations and chemical evolution models, the presence of efficient feedback mechanisms and/or extended star formation profiles on top of the classical "inside-out" scenario of mass assembly is generally required to reproduce the observed flatness of the metallicity gradients beyond z$\sim$1 . Three galaxies with significantly (> 3$\sigma$) "inverted" gradients are also found, showing an anti-correlation between metallicity and star formation rate density on local scales, possibly suggesting recent episodes of pristine gas accretion or strong radial flows in place. Nevertheless, the individual metallicity maps are characterised by a variety of different morphologies, with flat radial gradients sometimes hiding non-axisymmetric variations on kpc scales which are washed out by azimuthal averages, especially in interacting systems or in those undergoing local episodes of recent star formation.

The Hubble Frontier Fields offer an exceptionally deep window into the high-redshift universe, covering a substantially larger area than the Hubble Ultra-Deep field at low magnification and probing 1--2 mags deeper in exceptional high-magnification regions. We aim to leverage Atacama Large Millimetre Array (ALMA) band 6 ($\approx$263\,GHz) mosaics in the central portions of five Frontier Fields to characterize the infrared (IR) properties of $1582$ ultraviolet (UV)-selected Lyman-Break Galaxies (LBGs) at redshifts of $z {\sim}$2--8. We investigated individual and stacked fluxes and IR excess (IRX) values of the LBG sample as functions of stellar mass ($\mathrm{M}_{\bigstar}$), redshift, UV luminosity and slope $\beta$, and lensing magnification. Two (2) LBG candidates were individually detected above a significance of $4.1{-}\sigma$, while stacked samples of the remaining LBG candidates yielded no significant detections. We investigated our detections and upper limits in the context of the IRX-$\mathrm{M}_{\bigstar}$ and IRX-$\beta$ relations, probing at least one dex lower in stellar mass than past studies have done. Our upper limits exclude substantial portions of parameter space and they are sufficiently deep in a handful of cases to create mild tension with the typically assumed attenuation and consensus relations. We observe a clear and smooth trend between $\mathrm{M}_{\bigstar}$ and $\beta$, which extends to low masses and blue (low) $\beta$ values, consistent with expectations from previous works.

We report the first detection obtained with ALMA of the [N II] 122$\mu$m line emission from a galaxy group BRI 1202-0725 at $z=4.69$ consisting of a QSO and a submilimeter-bright galaxy (SMG). Combining with a detection of [N II] 205$\mu$m line in both galaxies, we constrain the electron densities of the ionized gas based on the line ratio of [NII]122/205. The derived electron densities are $26^{+12}_{-11}$ and $134^{+50}_{-39}$ cm$^{-3}$ for the SMG and the QSO, respectively. The electron density of the SMG is similar to that of the Galactic Plane and to the average of the local spirals. Higher electron densities by up to a factor of three could, however, be possible for systematic uncertainties of the line flux estimates. The electron density of the QSO is comparable to high-$z$ star-forming galaxies at $z=1.5-2.3$, obtained using rest-frame optical lines and with the lower limits suggested from stacking analysis on lensed starbursts at $z=1-3.6$ using the same tracer of [NII]. Our results suggest a large scatter of electron densities in global scale at fixed star formation rates for extreme starbursts. The success of the [N II] 122$\mu$m and 205$\mu$m detections at $z=4.69$ demonstrates the power of future systematic surveys of extreme starbursts at $z>4$ for probing the ISM conditions and the effects on surrounding environments.

[abridged] Characterizing the number counts of faint, dusty star-forming galaxies is currently a challenge even for deep, high-resolution observations in the FIR-to-mm regime. They are predicted to account for approximately half of the total extragalactic background light at those wavelengths. Searching for dusty star-forming galaxies behind massive galaxy clusters benefits from strong lensing, enhancing their measured emission while increasing spatial resolution. Derived number counts depend, however, on mass reconstruction models that properly constrain these clusters. We estimate the 1.1 mm number counts along the line of sight of three galaxy clusters, i.e. Abell 2744, MACSJ0416.1-2403 and MACSJ1149.5+2223, which are part of the ALMA Frontier Fields Survey. We perform detailed simulations to correct these counts for lensing effects. We use several publicly available lensing models for the galaxy clusters to derive the intrinsic flux densities of our sources. We perform Monte Carlo simulations of the number counts for a detailed treatment of the uncertainties in the magnifications and adopted source redshifts. We find an overall agreement among the number counts derived for the different lens models, despite their systematic variations regarding source magnifications and effective areas. Our number counts span ~2.5 dex in demagnified flux density, from several mJy down to tens of uJy. Our number counts are consistent with recent estimates from deep ALMA observations at a 3$\sigma$ level. Below $\approx$ 0.1 mJy, however, our cumulative counts are lower by $\approx$ 1 dex, suggesting a flattening in the number counts. In our deepest ALMA mosaic, we estimate number counts for intrinsic flux densities $\approx$ 4 times fainter than the rms level. This highlights the potential of probing the sub-10 uJy population in larger samples of galaxy cluster fields with deeper ALMA observations.

Most sub-mm emission line studies of galaxies to date have targeted sources with known redshifts where the frequencies of the lines are well constrained. Recent blind line scans circumvent the spectroscopic redshift requirement, which could represent a selection bias. Our aim is to detect emission lines present in continuum oriented observations. The detection of such lines provides spectroscopic redshift and yields properties of the galaxies. We perform a search for emission lines in the ALMA observations of five Frontier Fields clusters and assess the reliability of our detection by associating line candidates with detected galaxies in deep near-infrared imaging. We find 26 significant emission lines candidates, with observed line fluxes between 0.2-4.6 Jy km s$^{-1}$ and velocity dispersions (FWHM) of 25-600 km s$^{-1}$. Nine of these candidates lie nearby to near-infrared sources, boosting their reliability; in six cases the observed line frequency and strength are consistent with expectations given the photometric redshift and properties of the galaxy counterparts. We present redshift identifications, magnifications and molecular gas estimates for the galaxies with identified lines. We show that two of these candidates likely originate from starburst galaxies, one of which is a jellyfish galaxy, while another two are consistent with being main sequence galaxies based in their depletion times. This work highlights the degree to which serendipitous emission lines can be discovered in large mosaic continuum observations when deep ancillary data are available. The low number of high-significance line detections, however, confirms that such surveys are not as optimal as blind line scans. We stress that Monte Carlo simulations should be used to assess the line detections significances, since using the negative noise suffers from stochasticity and incurs larger uncertainties.

Dusty star-forming galaxies are among the most prodigious systems at high redshift (z>1), characterized by high star formation rates and huge dust reservoirs. The bright end of this population has been well characterized in recent years, but considerable uncertainties remain for fainter dusty star-forming galaxies, which are responsible for the bulk of star formation at high redshift and thus play a key role in galaxy growth and evolution. In this first paper of our series, we describe our methods for finding high redshift faint dusty galaxies using millimeter observations with ALMA. We obtained ALMA 1.1 mm mosaic images for three strong-lensing galaxy clusters from the Frontier Fields survey. The 2'x2' mosaics overlap with the deep HST WFC3/IR footprints and encompass the high magnification regions of each cluster. The combination of extremely high ALMA sensitivity and the magnification power of these clusters allows us to systematically probe the sub-mJy population of dusty star-forming galaxies over a large surveyed area. We present a description of the reduction and analysis of the ALMA continuum observations for the galaxy clusters Abell 2744 (z=0.308), MACSJ0416.1-2403 (z=0.396) and MACSJ1149.5+2223 (z=0.543), for which we reach observed rms sensitivities of 55, 59 and 71 $\mu$Jy beam$^{-1}$ respectively. We detect 12 dusty star-forming galaxies at S/N>5.0 across the three clusters, all of them presenting coincidence with NIR detected counterparts in the HST images. None of the sources fall close to the lensing caustics, thus they are not strongly lensed. The observed 1.1 mm flux densities for the total sample of galaxies range from 0.41 to 2.82 mJy, with observed effective radii spanning <0.05" to 0.37"$\pm$0.21". The lensing-corrected sizes of the detected sources appear to be in the same range as those measured in brighter samples, albeit with possibly larger dispersion.

We searched for z > 7 Lyman-break galaxies (LBGs) in the optical-to-mid-infrared Hubble Frontier Field and associated parallel field observations of the strong-lensing cluster MACS J0416-2403. We discovered 22 candidates, of which six lie at z > 9 and one lies at z > 10. Based on the Hubble and Spitzer photometry, all have secure photometric redshifts and a negligible probability of being at lower redshifts, according to their peak probability ratios, R. This substantial increase in the number of known high-redshift galaxies allows a solid determination of the luminosity function at z > 8. The number of high-z candidates in the parallel field is considerably higher than that in the Abell 2744 parallel field. Our candidates have median stellar masses of log(M_*) ~ 8.40^{+0.55}_{-0.31}~Msun, SFRs of ~ 1.6^{+0.5}_{-0.4} Msun yr^-1, and SFR-weighted ages of < 310^{+70}_{-140} Myr. Finally, we are able to put strong constraints on the z = 7,8,9 and 10 luminosity functions. One of the objects in the cluster field is a z ~ 10 candidate, with a magnification of mu ~ 20 +- 13. This object is likely the faintest z ~ 10 object known to date, allowing a first look into the extreme faint-end (L ~ 0.04L*) of the z ~ 10 luminosity function.

We used near-infrared integral field spectroscopic observations from the AMAZE and LSD programs to constrain the metallicity in a sample of 40 star forming galaxies at 3<z<5 (most of which at z~3.4). We measure metallicities by exploiting strong emission line diagnostics. We found that a significant fraction of star-forming galaxies at z~3.4 deviate from the Fundamental Metallicity Relation (FMR), with a metallicity up to a factor of ten lower than expected according to the FMR. This deviation does not correlate with the dynamical properties of the galaxy or with the presence of interactions. To investigate the origin of the metallicity deviations in more detail, we also infer information on the gas content, by inverting the Schmidt-Kennicutt relation. In agreement with recent CO observational data, we found that, in contrast with the steeply rising trend at 0<z<2, the gas fraction in massive galaxies remains constant, with indication of a marginal decline, at 2<z<3.5. When combined with the metallicity information, we infer that to explain both the low metallicity and gas content in z~3.4 galaxies, both prominent outflows and massive pristine gas inflows are needed. In ten galaxies we can also spatially resolve the metallicity distribution. We found that the metallicity generally anticorrelates with the distribution of star formation and with the gas surface density. We discuss these findings in terms of pristine gas inflows towards the center, and outflows of metal rich gas from the center toward the external regions. (Abridged)

We analyse 14 LBGs at z~2.8-3.8 constituting the only sample where both a spectroscopic measurement of their metallicity and deep IR observations (CANDELS+HUGS survey) are available. Fixing the metallicity of population synthesis models to the observed values, we determine best-fit physical parameters under different assumptions about the star-formation history and also consider the effect of nebular emission. For comparison we determine the UV slope of the objects, and use it to estimate their SFR_UV99 by correcting the UV luminosity following Meurer et al. (1999). A comparison between SFR obtained through SED-fitting (SFR_fit) and the SFR_UV99 shows that the latter are underestimated by a factor 2-10, regardless of the assumed SFH. Other SFR indicators (radio, far-IR, X-ray, recombination lines) coherently indicate SFRs a factor of 2-4 larger than SFR_UV99 and in closer agreement with SFR_fit. This discrepancy is due to the solar metallicity implied by the usual beta-A1600 conversion factor. We propose a refined relation, appropriate for sub-solar metallicity LBGs: A1600 = 5.32+1.99beta. This relation reconciles the dust-corrected UV with the SED-fitting and the other SFR indicators. We show that the fact that z~3 galaxies have sub-solar metallicity implies an upward revision by a factor of ~1.5-2 of the global SFRD, depending on the assumptions about the age of the stellar populations. We find very young best-fit ages (10-500 Myrs) for all our objects. From a careful examination of the uncertainties in the fit and the amplitude of the Balmer break we conclude that there is little evidence of the presence of old stellar population in at least half of the LBGs in our sample, suggesting that these objects are probably caught during a huge star-formation burst, rather than being the result of a smooth evolution.

We present a SINFONI integral field kinematical study of 33 galaxies at z~3 from the AMAZE and LSD projects which are aimed at studying metallicity and dynamics of high-redshift galaxies. The number of galaxies analyzed in this paper constitutes a significant improvement compared to existing data in the literature and this is the first time that a dynamical analysis is obtained for a relatively large sample of galaxies at z~3. 11 galaxies show ordered rotational motions (~30% of the sample), in these cases we estimate dynamical masses by modeling the gas kinematics with rotating disks and exponential mass distributions. We find dynamical masses in the range 2 \times 10^9 M\odot - 2 \times 10^11 M\odot with a mean value of ~ 2 \times 10^10 M\odot. By comparing observed gas velocity dispersion with that expected from models, we find that most rotating objects are dynamically "hot", with intrinsic velocity dispersions of the order of ~90 km s-1. The median value of the ratio between the maximum disk rotational velocity and the intrinsic velocity dispersion for the rotating objects is 1.6, much lower than observed in local galaxies value (~10) and slightly lower than the z~2 value (2 - 4). Finally we use the maximum rotational velocity from our modeling to build a baryonic Tully-Fisher relation at z~3. Our measurements indicate that z~3 galaxies have lower stellar masses (by a factor of ten on average) compared to local galaxies with the same dynamical mass. However, the large observed scatter suggests that the Tully-Fisher relation is not yet "in place" at these early cosmic ages, possibly due to the young age of galaxies. A smaller dispersion of the Tuly-Fisher relation is obtained by taking into account the velocity dispersion with the use of the S_0.5 indicator, suggesting that turbulent motions might have an important dynamical role.