Search

[Abridged] Our view of dust in primordial galaxies is limited towards a few tens of z~7 galaxies, pre-selected from UV-optical observations, and are thus not necessarily representative of the bulk of the sources at these redshifts. In this work, we aim at constraining the dust properties of galaxies at 69) at 6-19mag), lower mass sources (log Mstar<9), and at higher redshift (z>7). Fitting these 68, compared to the fit without FIR. We extend the LIR vs MUV relation down to MUV=-19mag and show a tentative breakdown of the relation at fainter UV magnitudes. The positions of the JADES z~6.5 sample on the IRX versus beta and IRX versus Mstar diagrams are consistent with those of the ALPINE (z~5.5) and REBELS (z~6.5) samples, suggesting that the dust composition and content of our mass-selected sample are similar to these UV-selected galaxies. Extending our analysis of the infrared properties to z>7 galaxies, we find a non-evolution of beta with redshift in the MUV range probed by our sample (-17.24+/-0.62) and highlight the fact that samples from the literature are not representative of the bulk of galaxy populations at z>6. We confirm a linear relation between AV and 1/sSFR with a flatter slope than previously reported due to the use of ALMA constraints. Our results suggest that rapid and significant dust production has already happened by z~7., Comment: Submitted to A&A

Early galaxy formation, initiated by the dark matter and gas assembly, evolves through frequent mergers and feedback processes into dynamically hot, chaotic structures. In contrast, dynamically cold, smooth rotating disks have been observed in massive evolved galaxies merely 1.4 billion years after the Big Bang, suggesting rapid morphological and dynamical evolution in the early Universe. Probing this evolution mechanism necessitates studies of young galaxies, yet efforts have been hindered by observational limitations in both sensitivity and spatial resolution. Here we report high-resolution observations of a strongly lensed and quintuply imaged, low-luminosity, young galaxy at $z=6.072$ (dubbed the Cosmic Grapes), 930 million years after the Big Bang. Magnified by gravitational lensing, the galaxy is resolved into at least 15 individual star-forming clumps with effective radii of $r_{\rm e}\simeq$ 10--60 parsec (pc), which dominate $\simeq$ 70\% of the galaxy's total flux. The cool gas emission unveils a smooth, underlying rotating disk characterized by a high rotational-to-random motion ratio and a gravitationally unstable state (Toomre $Q \simeq$ 0.2--0.3), with high surface gas densities comparable to local dusty starbursts with $\simeq10^{3-5}$ $M_{\odot}$/pc$^{2}$. These gas properties suggest that the numerous star-forming clumps are formed through disk instabilities with weak feedback effects. The clumpiness of the Cosmic Grapes significantly exceeds that of galaxies at later epochs and the predictions from current simulations for early galaxies. Our findings shed new light on internal galaxy substructures and their relation to the underlying dynamics and feedback mechanisms at play during their early formation phases, potentially explaining the high abundance of bright galaxies observed in the early Universe and the dark matter core-cusp problem., Comment: Submitted. 44 pages, 14 figures, 5 tables. Comments are welcome! See also the companion papers on arXiv. Valentino+2024: arXiv:2402.17845 Gim\'enez-Arteaga+2024: arXiv:2402.17875

We present JWST/NIRCam observations of a strongly-lensed, multiply-imaged galaxy at $z=6.072$, with magnification factors >~20 across the galaxy. We perform a spatially-resolved analysis of the physical properties at scales of ~200 pc, inferred from SED modelling of 5 NIRCam imaging bands on a pixel-by-pixel basis. We find young stars surrounded by extended older stellar populations. By comparing H$\alpha$+[NII] and [OIII]+H$\beta$ maps inferred from the image analysis with our additional NIRSpec IFU data, we find that the spatial distribution and strength of the line maps are in agreement with the IFU measurements. We explore different parametric SFH forms with Bagpipes on the spatially-integrated photometry, finding that a double power-law star formation history retrieves the closest value to the spatially-resolved stellar mass estimate, and other SFH forms suffer from the dominant outshining emission from the youngest stars, thus underestimating the stellar mass - up to ~0.5 dex-. On the other hand, the DPL cannot match the IFU measured emission lines. Additionally, the ionizing photon production efficiency may be overestimated in a spatially-integrated approach by ~0.15 dex, when compared to a spatially-resolved analysis. The agreement with the IFU measurements points towards the pixel-by-pixel approach as a way to mitigate the general degeneracy between the flux excess from emission lines and underlying continuum, especially when lacking photometric medium-band coverage and/or IFU observations. This study stresses the importance of studying galaxies as the complex systems that they are, resolving their stellar populations when possible, or using more flexible SFH parameterisations. This can aid our understanding of the early stages of galaxy evolution by addressing the challenge of inferring robust stellar masses and ionizing photon production efficiencies of high redshift galaxies., Comment: Submitted to A&A; 13 pages. See also the companion papers on arXiv today: Fujimoto+2024 and Valentino+2024

Faint, star-forming galaxies likely play a dominant role in cosmic reionisation. Strides have been made in recent years to characterise these populations at high redshifts ($z>3$). Now for the first time, with JWST photometry beyond 1$\,\mu m$ in the rest frame, we can derive accurate stellar masses and position these galaxies on the galaxy main sequence. We seek to assess the place of 96 individual Lyman-alpha emitters (LAEs) selected behind the A2744 lensing cluster with MUSE spectroscopy on the galaxy main sequence. We also compare derived stellar masses to Lyman-alpha luminosities and equivalent widths to better quantify the relationship between the Lyman-alpha emission and the host galaxy. These 96 LAEs lie in the redshift range $2.9

As the primary fuel for star formation, molecular gas plays a key role in galaxy evolution. A number of techniques have been used for deriving the mass of molecular reservoirs in the early Universe (e.g., [CII]158$\mu$m, [CI], dust continuum), but the standard approach of CO-based estimates has been limited to a small number of galaxies due to the intrinsic faintness of the line. We present Jansky Very Large Array (JVLA) observations of the $z\sim8.31$ galaxy MACS0416_Y1, targeting CO(2-1) and rest-frame radio continuum emission, which result in upper limits on both quantities. Adding our continuum limit to the published far-infrared (FIR) spectral energy distribution (SED), we find a small non-thermal contribution to the FIR emission, a low dust mass ($\rm\log_{10}(M_D/M_{\odot})\sim5$), and an abnormally high dust temperature ($\rm T_D\gtrsim90K$) that may indicate a recent starburst. Assuming a low metallicity ($Z/Z_{\odot}\sim0.25$), we find evidence for $M_{\rm H_2,CO}\lesssim10^{10}$M$_{\odot}$, in agreement with previous [CII] investigations ($M_{\rm H_2,[CII]}\sim10^{9.6}$M$_{\odot}$). Upcoming JWST observations of this source will result in a precise determination of $Z$, enabling better constraints and an unprecedented view of the gaseous reservoir in this primordial starburst galaxy., Comment: 6 pages, 3 figures, accepted for publication in MNRASL

The ALMA lensing cluster survey (ALCS) is a 96-hr large program dedicated to uncovering and characterizing intrinsically faint continuum sources and line emitters with the assistance of gravitational lensing. All 33 cluster fields were selected from HST/Spitzer treasury programs including CLASH, Hubble Frontier Fields, and RELICS, which also have Herschel and Chandra coverages. The total sky area surveyed reaches $\sim$133 arcmin$^2$ down to a depth of $\sim$60 $\mu$Jy beam$^{-1}$ (1$\sigma$) at 1.2 mm, yielding 141 secure blind detections of continuum sources and additional 39 sources aided by priors. We present scientific motivation, survey design, the status of spectroscopy follow-up observations, and number counts down to $\sim$7 $\mu$Jy. Synergies with JWST are also discussed., Comment: 6 pages, 4 figures, Proceedings of the 7th Chile-Cologne-Bonn-Symposium: Physics and Chemistry of Star Formation, V. Ossenkopf-Okada, R. Schaaf, I. Breloy (eds.)

Early JWST spectroscopic campaigns have confirmed the presence of strong [O III] line-emitting galaxies in the redshift interval $7

We present EMIR, a powerful near-infrared (NIR) camera and multi-object spectrograph (MOS) installed at the Nasmyth focus of the 10.4 m GTC. EMIR was commissioned in mid-2016 and is offered as a common-user instrument. It provides spectral coverage of 0.9 to 2.5 $\mu m$ over a field of view (FOV) of 6.67x6.67 squared arcmin in imaging mode, and 6.67x4 squared arcmin in spectroscopy. EMIR delivers up to 53 spectra of different objects thanks to a robotic configurable cold slit mask system that is located inside the cryogenic chamber, allowing rapid reconfiguration of the observing mask. The imaging mode is attained by moving all bars outside the FOV and then leaving an empty space in the GTC focal surface. The dispersing suite holds three large pseudo-grisms, formed by the combination of high-efficiency FuSi ion-etched ruled transmission grating sandwiched between two identical ZnSe prisms, plus one standard replicated grism. These dispersing units offer the spectral recording of an atmospheric window $J,H,K$ in a single shot with resolving powers of 5000, 4250, 4000, respectively for a nominal slit width of 0.6\arcsec, plus the combined bands $YJ$ or $HK$, also in a single shot, with resolution of $\sim$ 1000. The original Hawaii2 FPA detector, which is prone to instabilities that add noise to the signal, is being replaced by a new Hawaii2RG detector array, and is currently being tested at the IAC. This paper presents the most salient features of the instrument, with emphasis on its observing capabilities and the functionality of the configurable slit unit. Sample early science data is also shown., Comment: 10 pages, 11 figures

According to the current paradigm of galaxy formation, the first galaxies have been likely formed within large dark matter haloes. The fragmentation of these massive haloes led to the formation of galaxy protoclusters, which are usually composed of one to a few bright objects, surrounded by numerous fainter (and less massive) galaxies. These early structures could have played a major role in reionising the neutral hydrogen within the first billion years of the Universe; especially, if their number density is significant.Taking advantage of the unprecedented sensitivity reached by the \textit{James Webb Space Telescope (JWST)}, galaxy protoclusters can now be identified and studied in increasing numbers beyond $z\geq\ $6. Characterising their contribution to the UV photon budget could supply new insights into the reionisation process. We analyse the first JWST dataset behind SMACS0723-7327 to search for protoclusters at $z\geq6$, combining the available spectroscopic and photometric data. We then compare our findings with semi-analytical models and simulations. In addition to two bright galaxies ($\leq$26.5 AB in F277W), separated by $\sim$11\arcsec and spectroscopically confirmed at $z_{spec}=7.66$, we identify 6 additional galaxies with similar colors in a $\theta\sim20$\arcsec radius around these (corresponding to R$\sim60-90$ kpc in the source plane). Using several methods, we estimate the mass of the dark matter halo of this protocluster, $\sim$3.3$\times$10$^{11}$M$_{\odot}$ accounting for magnification, consistent with various predictions. The physical properties of all protocluster members are also in excellent agreement with what has been previously found at lower redshifts: star-formation main sequence and protocluster size. This detection adds to just a few protoclusters currently known in the first billion years of the universe., Comment: 7 pages, 6 Figures. Accepted for publication in A&A Letter

Massive galaxy clusters are interesting astrophysical and cosmological study objects, but are relatively rare. In the redshift range z = 0.25 to 0.5 which is, for example, a favourable region for gravitational lensing studies, about 100 such systems are known. Most of them have been studied in X-rays. In this paper we study the six remaining massive clusters in this redshift interval in the highly complete CLASSIX survey which have so far not been observed with sufficiently deep exposures in X-rays. With data from our new XMM-Newton observations we characterise their structures, derive X-ray properties such as the X-ray luminosity and intra-cluster medium temperature and estimate their gas and total masses. We find that one cluster, RXCJ1230.7+3439, is dynamically young with three distinct substructures in the cluster outskirts and RXCJ1310.9+2157/RXCJ1310.4+2151 is a double cluster system. Mass determination is difficult in the systems with substructure. We therefore discuss several methods of mass estimation including scaling relations. In summary we find that five of the six study targets are indeed massive clusters as expected, while the last cluster RXCJ2116.2-0309 is a close projection of a distant and a nearby cluster which has led to a previous overestimation of its mass. In the XMM-Newton observation fields we also find three low redshift clusters close to the targets which are also analysed and described here. In the field of RXCJ2116.2-0309 we discover serendipitously a highly variable X-ray source which has decreased its flux within a year by more than a factor of eight. This source is most probably an AGN., Comment: 15 pages with 26 figures, Astronomy and Astrophysics 2022 (in press)

Before the end of the epoch of reionization, the Hydrogen in the Universe was predominantly neutral. This leads to a strong attenuation of Ly$\alpha$ lines of $z\gtrsim6$ galaxies in the intergalactic medium. Nevertheless, Ly$\alpha$ has been detected up to very high redshifts ($z\sim9$) for several especially UV luminous galaxies. Here, we test to what extent the galaxy's local environment might impact the Ly$\alpha$ transmission of such sources. We present an analysis of dedicated Hubble Space Telescope (HST) imaging in the CANDELS/EGS field to search for fainter neighbours around three of the most UV luminous and most distant spectroscopically confirmed Ly$\alpha$ emitters: EGS-zs8-1, EGS-zs8-2 and EGSY-z8p7 at $z_\mathrm{spec}=7.73$, 7.48, and 8.68, respectively. We combine the multi-wavelength HST imaging with Spitzer data to reliably select $z\sim7-9$ galaxies around the central, UV-luminous sources. In all cases, we find a clear enhancement of neighbouring galaxies compared to the expected number in a blank field (by a factor $\sim 3-9\times$). Our analysis thus reveals ubiquitous overdensities around luminous Ly$\alpha$ emitting sources in the heart of the cosmic reionization epoch. We show that our results are in excellent agreement with expectations from the Dragons simulation, confirming the theoretical prediction that the first ionized bubbles preferentially formed in overdense regions. JWST follow-up observations of the neighbouring galaxies identified here will be needed to confirm their physical association and to map out the ionized regions produced by these sources., Comment: 13 pages, 14 figures, Submitted to MNRAS

We discuss the spectral energy distributions and physical properties of six galaxies whose photometric redshifts suggest they lie beyond a redshift $z\simeq$9. Each was selected on account of a prominent excess seen in the Spitzer/IRAC 4.5$\mu$m band which, for a redshift above $z=9.0$, likely indicates the presence of a rest-frame Balmer break and a stellar component that formed earlier than a redshift $z\simeq10$. In addition to constraining the earlier star formation activity on the basis of fits using stellar population models with BAGPIPES, we have undertaken the necessary, but challenging, follow-up spectroscopy for each candidate using various combinations of Keck/MOSFIRE, VLT/X-shooter, Gemini/FLAMINGOS2 and ALMA. Based on either Lyman-$\alpha$ or [OIII] 88 $\mu$m emission, we determine a convincing redshift of $z$=8.78 for GN-z-10-3 and a likely redshift of $z$=9.28 for the lensed galaxy MACS0416-JD. For GN-z9-1, we conclude the case remains promising for a source beyond $z\simeq$9. Together with earlier spectroscopic data for MACS1149-JD1, our analysis of this enlarged sample provides further support for a cosmic star formation history extending beyond redshifts $z\simeq$10. We use our best-fit stellar population models to reconstruct the past rest-frame UV luminosities of our sources and discuss the implications for tracing earlier progenitors of such systems with the James Webb Space Telescope., Comment: updated version with correction of a numerical error on MACS1149-JD1 age estimate. The predicted earlier luminosity for this galaxy has been corrected on Figure 11

We report the discovery of an intrinsically faint, quintuply-imaged, dusty galaxy MACS0600-z6 at a redshift $z=$6.07 viewed through the cluster MACSJ0600.1-2008 ($z$=0.46). A $\simeq4\sigma$ dust detection is seen at 1.2mm as part of the ALMA Lensing Cluster Survey (ALCS), an on-going ALMA Large program, and the redshift is secured via [C II] 158 $\mu$m emission described in a companion paper. In addition, spectroscopic follow-up with GMOS/Gemini-North shows a break in the galaxy's spectrum, consistent with the Lyman break at that redshift. We use a detailed mass model of the cluster and infer a magnification $\mu\gtrsim$30 for the most magnified image of this galaxy, which provides an unprecedented opportunity to probe the physical properties of a sub-luminous galaxy at the end of cosmic reionisation. Based on the spectral energy distribution, we infer lensing-corrected stellar and dust masses of $\rm{2.9^{+11.5}_{-2.3}\times10^9}$ and $\rm{4.8^{+4.5}_{-3.4}\times10^6}$ $\rm{M_{\odot}}$ respectively, a star formation rate of $\rm{9.7^{+22.0}_{-6.6} M_{\odot} yr^{-1}}$, an intrinsic size of $\rm{0.54^{+0.26}_{-0.14}}$ kpc, and a luminosity-weighted age of 200$\pm$100 Myr. Strikingly, the dust production rate in this relatively young galaxy appears to be larger than that observed for equivalent, lower redshift sources. We discuss if this implies that early supernovae are more efficient dust producers and the consequences for using dust mass as a probe of earlier star formation., Comment: submitted to MNRAS, revised version after including referee's comments

We present a study of the intersection between the populations of star forming galaxies (SFGs) selected as either Lyman break galaxies (LBGs) or Lyman-alpha emitters (LAEs) in the redshift range 2.9 - 6.7, within the same volume of universe sampled by the Multi-Unit Spectroscopic Explorer (MUSE) behind the Hubble Frontier Fields lensing cluster A2744. We define three samples of star-forming galaxies: LBG galaxies with an LAE counterpart (92 galaxies), LBG galaxies without LAE counterpart (408 galaxies) and LAE galaxies without an LBG counterpart (46 galaxies). All these galaxies are intrinsically faint due to the lensing nature of the sample (Muv $\ge$ -20.5). The fraction of LAEs among all selected SFGs increases with redshift up to z $\sim$ 6 and decreases for higher redshifts. The evolution of LAE/LBG populations with UV magnitude and Lya luminosity shows that the LAE selection is able to identify intrinsically UV faint galaxies with Muv $\ge$ -15 that are typically missed in the deepest lensing photometric surveys. The LBG population seems to fairly represent the total population of SFGs down to Muv$\sim$-15. Galaxies with Muv$<-17$ tend to have SFRLya$<$SFRuv, whereas the opposite trend is observed within our sample for faint galaxies with Muv$>-17$, including galaxies only detected by their Lya emission, with a large scatter. These trends, previously observed in other samples of SFGs at high-$z$, are seen here for very faint Muv$\sim -15$ galaxies, much fainter than in previous studies. There is no clear evidence, based on the present results, for an intrinsic difference on the properties of the two populations selected as LBG and/or LAE. The observed trends could be explained by a combination of several facts, like the existence of different star-formation regimes, the dust content, the relative distribution and morphology of dust and stars, or the stellar populations

The scatter in the relationship between the strength of [CII] 158$\mu$m emission and the star formation rate at high-redshift has been the source of much recent interest. Although the relationship is well-established locally, several intensely star-forming galaxies have been found whose [CII] 158$\mu$m emission is either weak, absent or spatially offset from the young stars. Here we present new ALMA data for the two most distant, gravitationally-lensed and spectroscopically-confirmed galaxies, A2744\_YD4 at $z=$8.38 and MACS1149\_JD1 at $z=$9.11, both of which reveal intense [OIII] 88$\mu$m emission. In both cases we provide stringent upper limits on the presence of [CII] 158$\mu$m with respect to [OIII] 88$\mu$m. We review possible explanations for this apparent redshift-dependent [CII] deficit in the context of our recent hydrodynamical simulations. Our results highlight the importance of using several emission line diagnostics with ALMA to investigate the nature of the interstellar medium in early galaxies., Comment: 5 pages, 5 figures. Accepted for publication in MNRAS

We present the results obtained with VLT/MUSE on the faint-end of the Lyman-alpha luminosity function (LF) based on deep observations of four lensing clusters. The precise aim of the present study is to further constrain the abundance of Lyman-alpha emitters (LAEs) by taking advantage of the magnification provided by lensing clusters. We blindly selected a sample of 156 LAEs, with redshifts between $2.9 \le z \le 6.7$ and magnification-corrected luminosities in the range $ 39 \lesssim \log L_{Ly_{\alpha}}$ [erg s$^{-1}$] $\lesssim 43$. The price to pay to benefit from magnification is a reduction of the effective volume of the survey, together with a more complex analysis procedure. To properly take into account the individual differences in detection conditions (including lensing configurations, spatial and spectral morphologies) when computing the LF, a new method based on the 1/Vmax approach was implemented. The LAE LF has been obtained in four different redshift bins with constraints down to $\log L_{Ly_{\alpha}} = 40.5$. From our data only, no significant evolution of LF mean slope can be found. When performing a Schechter analysis including data from the literature to complete the present sample a steep faint-end slope was measured varying from $\alpha = -1.69^{+0.08}_{-0.08}$ to $\alpha = -1.87^{+0.12}_{-0.12}$ between the lowest and the highest redshift bins. The contribution of the LAE population to the star formation rate density at $z \sim 6$ is $\lesssim 50$% depending on the luminosity limit considered, which is of the same order as the Lyman-break galaxy (LBG) contribution. The evolution of the LAE contribution with redshift depends on the assumed escape fraction of Lyman-alpha photons, and appears to slightly increase with increasing redshift when this fraction is conservatively set to one. (abridged), Comment: 32 pages, 22 figures

The aim of this paper is to introduce the WIRCam Ultra Deep Survey (WUDS), a near-IR photometric public survey carried out at the CFH Telescope in the field of the CFHTLS-D3 field (Groth Strip). WUDS includes four near-IR bands (Y, J, H and K_s) over a field of view of ~400 arcmin^2. The typical depth of WUDS data reaches between ~26.8 in Y and J, and ~26 in H and K_s (AB, 3 sigma in 1.3 arcsec aperture). The area and depth of this survey were specifically tailored to set strong constraints on the cosmic star formation rate and the luminosity function brighter or around L* in the z~6-10 redshift domain, although these data are also useful for a variety of extragalactic projects.This first paper is intended to present the properties of WUDS: catalog building, completeness and depth, number counts, photometric redshifts, and global properties of the galaxy population. We have also concentrated on the study of galaxy samples at z~[4.5-7] in this field. UV luminosity functions were derived at z~5 and z~6 taking advantage from the fact that WUDS covers a particularly interesting regime at intermediate luminosities, which allows a combined determination of M* and Phi* with increased accuracy. Our results on the luminosity function are consistent with a small evolution of both M* and Phi* between z=5 and z=6, irrespective of the method used to derive them, either photometric redshifts applied to blindly-selected dropout samples or the classical Lyman Break Galaxy color-preselected samples. Our results lend support to higher Phi* determinations at z=6 than usually reported. The selection and combined analysis of different galaxy samples at z>7 will be presented in a forthcoming paper. WUDS is intended to provide a robust database in the near-IR for the selection of targets for detailed spectroscopic studies, in particular for the EMIR/GTC GOYA Survey (Abridged), Comment: Accepted for publication in Astronomy & Astrophysics. 17 pages, 11 figures

[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., Comment: 19 pages, 14 figures, 3 tables. Accepted for publication in A&A

[abridged] The Hubble and Spitzer Space Telescope surveys of the Frontier Fields (FF) provide extremely deep images around six massive, strong-lensing clusters of galaxies. The ALMA FF survey aims to cover the same fields at 1.1mm, with maps reaching (unlensed) sensitivities of $<$70$\mu$Jy, in order to explore the properties of background dusty star-forming galaxies. We report on the multi-wavelength photometric analysis of all 12 significantly detected ($>$5$\sigma$) sources in the first three FF clusters observed by ALMA, based on data from Hubble and Spitzer, the VLT and Herschel. We measure the total photometry in all available bands and determine the photometric redshifts and the physical properties of the counterparts via SED-fitting. In particular, we carefully estimate the FIR photometry using 1.1mm priors to limit the misidentification of blended FIR counterparts, which strongly affect some flux estimates in previous FIR catalogs. We identify robust near-infrared (NIR) counterparts for all 11 sources with K$_s$ detection, the majority of which are quite red, with eight having $F814W-K_s\gtrsim 4$ and five having $F160W-[4.5]\gtrsim3$. From the FIR point of view, all our objects have $z_{phot}$$\sim$1--3, whereas based on the optical SED one object prefers a high-$z$ solution ($z\geq\ $7). Five objects among our sample have spectroscopic redshifts from the GLASS survey for which we can reproduce their SEDs with existing templates. This verification confirms the validity of our photometric redshift methodology. The mean redshift of our sample is $z_{phot}$=1.99$\pm$0.27. All 1.1mm selected objects are massive (10.0$<\log[M_{\star}(M_{\odot})]<$ 11.5), with high star formation rates ($<\log[SFR(M_{\odot}/yr)]> \approx$1.6) and high dust contents (8.1 $<\log[M_{dust} (M_{\odot})]<$8.8), consistent with previous ALMA surveys., Comment: 18 pages, 11 figures, accepted for publication in A&A

We present \ci\,(2--1) and multi-transition $^{12}$CO observations of a dusty star-forming galaxy, ACT\,J2029+0120, which we spectroscopically confirm to lie at $z$\,=\,2.64. We detect CO(3--2), CO(5--4), CO(7--6), CO(8--7), and \ci\,(2--1) at high significance, tentatively detect HCO$^{+}$(4--3), and place strong upper limits on the integrated strength of dense gas tracers (HCN(4--3) and CS(7--6)). Multi-transition CO observations and dense gas tracers can provide valuable constraints on the molecular gas content and excitation conditions in high-redshift galaxies. We therefore use this unique data set to construct a CO spectral line energy distribution (SLED) of the source, which is most consistent with that of a ULIRG/Seyfert or QSO host object in the taxonomy of the \textit{Herschel} Comprehensive ULIRG Emission Survey. We employ RADEX models to fit the peak of the CO SLED, inferring a temperature of T$\sim$117 K and $n_{\text{H}_2}\sim10^5$ cm$^{-3}$, most consistent with a ULIRG/QSO object and the presence of high density tracers. We also find that the velocity width of the \ci\ line is potentially larger than seen in all CO transitions for this object, and that the $L'_{\rm C\,I(2-1)}/L'_{\rm CO(3-2)}$ ratio is also larger than seen in other lensed and unlensed submillimeter galaxies and QSO hosts; if confirmed, this anomaly could be an effect of differential lensing of a shocked molecular outflow., Comment: Accepted for publication in ApJ

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., Comment: new version accepted for publication in A&A

We search for high-redshift dropout galaxies behind the Hubble Frontier Fields (HFF) galaxy cluster MACS J1149.5+2223, a powerful cosmic lens that has revealed a number of unique objects in its field. Using the deep images from the Hubble and Spitzer space telescopes, we find 11 galaxies at z>7 in the MACS J1149.5+2223 cluster field, and 11 in its parallel field. The high-redshift nature of the bright z~9.6 galaxy MACS1149-JD, previously reported by Zheng et al., is further supported by non-detection in the extremely deep optical images from the HFF campaign. With the new photometry, the best photometric redshift solution for MACS1149-JD reduces slightly to z=9.44 +/- 0.12. The young galaxy has an estimated stellar mass of (7 +/- 2)X10E8 Msun, and was formed at z=13.2 +1.9-1.6 when the universe was ~300 Myr old. Data available for the first four HFF clusters have already enabled us to find faint galaxies to an intrinsic magnitude of M(UV) ~ -15.5, approximately a factor of ten deeper than the parallel fields., Comment: To appear in the Astrophysical Journal. 20 pages, 6 figures, 3 tables. arXiv admin note: text overlap with arXiv:1402.6743

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., Comment: Accepted for publication in A&A. New version addressing last referee comments, 17 pages and 12 figures

In this paper we present the results of our search for and study of $z \gtrsim 6$ galaxy candidates behind the third Frontier Fields (FF) cluster, MACSJ0717.5+3745, and its parallel field, combining data from Hubble and Spitzer. We select 39 candidates using the Lyman Break technique, for which the clear non-detection in optical make the extreme mid-$z$ interlopers hypothesis unlikely. We also take benefit from $z \gtrsim 6$ samples selected using previous Frontier Fields datasets of Abell 2744 and MACS0416 to improve the constraints on the properties of very high-redshift objects. We compute the redshift and the physical properties, such emission lines properties, star formation rate, reddening, and stellar mass for all Frontier Fields objects from their spectral energy distribution using templates including nebular emission lines. We study the relationship between several physical properties and confirm the trend already observed in previous surveys for evolution of star formation rate with galaxy mass, and between the size and the UV luminosity of our candidates. The analysis of the evolution of the UV Luminosity Function with redshift seems more compatible with an evolution of density. Moreover, no robust $z\ge$8.5 object is selected behind the cluster field, and few $z$$\sim$9 candidates have been selected in the two previous datasets from this legacy survey, suggesting a strong evolution in the number density of galaxies between $z$$\sim$8 and 9. Thanks to the use of the lensing cluster, we study the evolution of the star formation rate density produced by galaxies with L$>$0.03L$^{\star}$, and confirm the strong decrease observed between $z$$\sim$8 and 9., Comment: 21 pages - Accepted for publication in ApJ - v2: small corrections

We describe the search for Lyman-break galaxies (LBGs) near the sub-millimeter bright starburst galaxy HFLS3 at $z$$=$6.34 and a study on the environment of this massive galaxy during the end of reionization.We performed two independent selections of LBGs on images obtained with the \textit{Gran Telescopio Canarias} (GTC) and the \textit{Hubble Space Telescope} (HST) by combining non-detections in bands blueward of the Lyman-break and color selection. A total of 10 objects fulfilling the LBG selection criteria at $z$$>$5.5 were selected over the 4.54 and 55.5 arcmin$^2$ covered by our HST and GTC images, respectively. The photometric redshift, UV luminosity, and the star-formation rate of these sources were estimated with models of their spectral energy distribution. These $z$$\sim$6 candidates have physical properties and number densities in agreement with previous results. The UV luminosity function at $z$$\sim$6 and a Voronoi tessellation analysis of this field shows no strong evidence for an overdensity of relatively bright objects (m$_{F105W}$$<$25.9) associated with \textit{HFLS3}. However, the over-density parameter deduced from this field and the surface density of objects can not excluded definitively the LBG over-density hypothesis. Moreover we identified three faint objects at less than three arcseconds from \textit{HFLS3} with color consistent with those expected for $z$$\sim$6 galaxies. Deeper data are needed to confirm their redshifts and to study their association with \textit{HFLS3} and the galaxy merger that may be responsible for the massive starburst., Comment: 14 pages, 12 figures, accepted for publication in ApJ

(Abridged) The HST Frontier Fields project started at the end of 2013 with the aim of providing extremely deep images of 6 massive galaxy clusters. One of the main goals of this program is to push several telescopes to their limits in order to provide the best current view of the earliest stages of the Universe. We present a detailed analysis of $z$$\sim$8 objects behind the HFFs lensing cluster, MACS0416-2403, combining 0.3-1.6 $\mu$m imaging from HST, ground-based $K_s$ imaging from VLT HAWK-I, and 3.6 $\mu$m and 4.5 $\mu$m Spitzer Space Telescope. We apply the classical Lyman Break technique, combining non-detection criteria and color-selection. In order to avoid contamination by mid-$z$ interlopers, we require a strong break between optical and near infrared data. The luminosity function at $z$$\sim$8 is computed using a MC method taking advantage of the SED-fitting results. A piece of cautionary information is gleaned from new deep optical photometry of a previously identified $z$$\sim$8 galaxy in this cluster, which is now firmly detected as a mid-$z$ interloper with a strong $\approx1.5$ mag Balmer break (between F606W and F125W). Using the SED of this interloper, we estimate the contamination rate of our MACS0416$-$2403 sample, as well as for previous samples in Abell 2744 based on HFF data. Our selection recovers 4 robust objects with m$_{F160W}$ ranging from 26.0 to 27.9 AB and located in modest amplification regions ($\mu<$2.4). Two of the objects display a secondary break between the IRAC 3.6 $\mu$m and 4.5 $\mu$m bands which could be associated to the Balmer break or emission lines at $z$$\sim$8. The candidates generally have star formation rates around $\sim$10 M$_{\odot}$/yr and sizes ranging from 0.2 to 0.5 kpc, in good agreement with previous observations and expectations for objects in the early Universe., Comment: 11 pages, 8 figures, accepted for publication in A&A

We discuss the restframe UV emission from the starbursting galaxy HFLS3 at z=6.34, discovered in Herschel/SPIRE data due to its red color in the submm wavelengths from 250-500 um. The apparent inst. SFR of HFLS3 inferred from the total FIR luminosity measured with over 15 photometric data points between 100 to 1000 um is 2900 Msun/yr. Keck/NIRC2 Ks band adaptive optics imaging data showed two potential NIR counterparts near HFLS3. Previously, the northern galaxy was taken to be in the foreground at z=2.1 while the southern galaxy was assumed to HFLS3's NIR counterpart. New HST/WFC3 and ACS imaging data show both optically bright galaxies are in the foreground at z<6. A new lensing model based on HST data and mm-wave continuum emission yields a magnification of 2.2+/-0.3. The lack of multiple imaging constrains the lensing magnification to be lower than either 2.7 or 3.5 at the 95% confidence level for the two scenarios, which attribute one or two components to HFLS3 in the source plane. Correcting for gravitational lensing, the inst. SFR is 1320 Msun/yr with the 95% confidence lower limit around 830 Msun/yr. Using models for the restframe UV to FIR SED, the ave. SFR over the last 100 Myr is around 660 Msun/yr. The dust and stellar masses of HFLS3 from the same SED models are 3x10^8 Msun and ~5x10^10 Msun, respectively, with large systematic uncertainties on assumptions related to the SED model. With HST/WFC3 images we also find diffuse NIR emission about 0.5" (~3 kpc) SW of HFLS3 that remains undetected in the ACS data. The emission has a photometric redshift consistent with either z~6 or a dusty galaxy template at z~2. If at the same redshift as HFLS3 the detected diffuse emission could be part of the complex merger system that could be triggering the starburst. Alternatively, it could be part of the foreground structure at z~2.1 that is responsible for lensing of HFLS3., Comment: Submitted to ApJ, 10 pages, 6 figures

We describe a search for submillimeter emission in the vicinity of one of the most distant, luminous galaxies known, HerMES FLS3 at z=6.34, exploiting it as a signpost to a potentially biased region of the early Universe, as might be expected in hierarchical structure formation models. Imaging to the confusion limit with the innovative, wide-field submillimeter bolometer camera, SCUBA-2, we are sensitive to colder and/or less luminous galaxies in the surroundings of HFLS3. We use the Millennium Simulation to illustrate that HFLS3 may be expected to have companions if it is as massive as claimed, but find no significant evidence from the surface density of SCUBA-2 galaxies in its vicinity, or their colors, that HFLS3 marks an over-density of dusty, star-forming galaxies. We cannot rule out the presence of dusty neighbours with confidence, but deeper 450-um imaging has the potential to more tightly constrain the redshifts of nearby galaxies, at least one of which likely lies at z>~5. If associations with HFLS3 can be ruled out, this could be taken as evidence that HFLS3 is less biased than a simple extrapolation of the Millennium Simulation may imply. This could suggest either that it represents a rare short-lived, but highly luminous, phase in the evolution of an otherwise typical galaxy, or that this system has suffered amplification due to a foreground gravitational lens and so is not as intrinsically luminous as claimed., Comment: 6 pages; 4 colour figures

We present in this letter the first analysis of a z~8 galaxy candidate found in the Hubble and Spitzer imaging data of Abell 2744, as part of the Hubble Frontier Fields legacy program. We applied the most commonly-used methods to select exceptionally high-z galaxies by combining non-detection and color-criteria using seven HST bands. We used GALFIT on IRAC images for fitting and subtracting contamination of bright nearby sources. The physical properties have been inferred from SED-fitting using templates with and without nebular emission. This letter is focussed on the brightest candidate we found (m$_{F160W}$=26.2) over the 4.9 arcmin$^2$ field of view covered by the WFC3. It shows a non-detection in the ACS bands and at 3.6{\mu}m whereas it is clearly detected at 4.5{\mu}m with rather similar depths. This break in the IRAC data could be explained by strong [OIII]+H{\beta} lines at z~8 which contribute to the 4.5{\mu}m photometry. The best photo-z is found at z~8.0$^{+0.2}_{-0.5}$, although solutions at low-redshift (z~1.9) cannot be completely excluded, but they are strongly disfavoured by the SED-fitting work. The amplification factor is relatively small at {\mu}=1.49$\pm$0.02. The Star Formation Rate in this object is ranging from 8 to 60 Mo/yr, the stellar mass is in the order of M$_{\star}$=(2.5-10) x 10$^{9}$Mo and the size is r~0.35$\pm$0.15 kpc. This object is one of the first z~8 LBG candidates showing a clear break between 3.6{\mu}m and 4.5{\mu}m which is consistent with the IRAC properties of the first spectroscopically confirmed galaxy at a similar redshift. Due to its brightness, the redshift of this object could potentially be confirmed by near infrared spectroscopy with current 8-10m telescopes. The nature of this candidate will be revealed in the coming months with the arrival of new ACS and Spitzer data, increasing the depth at optical and near-IR wavelengths., Comment: 4 pages, 2 figures, Accepted for publication in Astronomy and Astrophysics Letters

We present a method for selecting $z>4$ dusty, star forming galaxies (DSFGs) using Herschel/SPIRE 250/350/500 $\mu m$ flux densities to search for red sources. We apply this method to 21 deg$^2$ of data from the HerMES survey to produce a catalog of 38 high-$z$ candidates. Follow-up of the first 5 of these sources confirms that this method is efficient at selecting high-$z$ DSFGs, with 4/5 at $z=4.3$ to $6.3$ (and the remaining source at $z=3.4$), and that they are some of the most luminous dusty sources known. Comparison with previous DSFG samples, mostly selected at longer wavelengths (e.g., 850 $\mu m$) and in single-band surveys, shows that our method is much more efficient at selecting high-$z$ DSFGs, in the sense that a much larger fraction are at $z>3$. Correcting for the selection completeness and purity, we find that the number of bright ($S_{500\,\mu m} \ge 30$ mJy), red Herschel sources is $3.3 \pm 0.8$ deg$^{-2}$. This is much higher than the number predicted by current models, suggesting that the DSFG population extends to higher redshifts than previously believed. If the shape of the luminosity function for high-$z$ DSFGs is similar to that at $z\sim2$, rest-frame UV based studies may be missing a significant component of the star formation density at $z=4$ to $6$, even after correction for extinction., Comment: Accepted for publication in ApJ

Strong gravitational lenses are now being routinely discovered in wide-field surveys at (sub)millimeter wavelengths. We present Submillimeter Array (SMA) high-spatial resolution imaging and Gemini-South and Multiple Mirror Telescope optical spectroscopy of strong lens candidates discovered in the two widest extragalactic surveys conducted by the Herschel Space Observatory: the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). From a sample of 30 Herschel sources with S_500>100 mJy, 21 are strongly lensed (multiply imaged), 4 are moderately lensed (singly imaged), and the remainder require additional data to determine their lensing status. We apply a visibility-plane lens modeling technique to the SMA data to recover information about the masses of the lenses as well as the intrinsic (i.e., unlensed) sizes (r_half) and far-infrared luminosities (L_FIR) of the lensed submillimeter galaxies (SMGs). The sample of lenses comprises primarily isolated massive galaxies, but includes some groups and clusters as well. Several of the lenses are located at z_lens>0.7, a redshift regime that is inaccessible to lens searches based on Sloan Digital Sky Survey spectroscopy. The lensed SMGs are amplified by factors that are significantly below statistical model predictions given the 500um flux densities of our sample. We speculate that this may reflect a deficiency in our understanding of the intrinsic sizes and luminosities of the brightest SMGs. The lensed SMGs span nearly one decade in L_FIR (median L_FIR=7.9x10^12 L_sun) and two decades in FIR luminosity surface density (median Sigma_FIR=6.0x10^11 L_sun kpc^-2). The strong lenses in this sample and others identified via (sub-)mm surveys will provide a wealth of information regarding the astrophysics of galaxy formation and evolution over a wide range in redshift., Comment: 28 pages, 7 tables, 9 figures. Accepted to ApJ; revised following minor comments from the referee and the community at large

We report new constraints on the galaxy luminosity function at z~9 based on observations carried out with ESO/VLT FORS2, HAWK-I and X-Shooter around the lensing cluster A2667, as part of our project aimed at selecting z~7-10 candidates accessible to spectroscopy. Only one J-dropout source was selected in this field fulfilling the color and magnitude criteria. This source was recently confirmed as a mid-z interloper based on X-Shooter spectroscopy. The depth and the area covered by our survey are well suited to set strong constraints on the bright-end of the galaxy luminosity function and hence on the star formation history at very high redshift. The non-detection of reliable J-dropout sources over the ~36arcmin2 field of view towards A2667 was used to carefully determine the lens-corrected effective volume and the corresponding upper-limit on the density of sources. The strongest limit is obtained for Phi(M_{1500}=-21.4+/-0.50)<6.70x10^{-6}Mpc^{-3}mag^{-1} at z~9. A maximum-likelihood fit of the luminosity function using all available data points including the present new result yields M*>-19.7 with fixed alpha=-1.74 and Phi*=1.10x10^{-3}Mpc^{-3}. The corresponding star formation rate density should be rho_{SFR}<5.97x10^{-3}M_{solar}/yr/Mpc^{3} at z~9. These results are in good agreement with the most recent estimates already published in this range of redshift and for this luminosity domain. This new result confirms the decrease in the density of luminous galaxies at very high-redshift, hence providing strong constraints for the design of future surveys aiming to explore the very high-redshift Universe., Comment: 4 pages, 2 figures, accepted for publication in A&A Letters

The spectral energy distributions (SED) of dusty galaxies at intermediate redshift may look similar to very high redshift galaxies in the optical/near infrared (NIR) domain. This can lead to the contamination of high redshift galaxy searches based on broad band optical/NIR photometry by lower redshift dusty galaxies as both kind of galaxies cannot be distinguished. The contamination rate could be as high as 50%. {This work shows how the far infrared (FIR) domain can help to recognize likely low-z interlopers in an optical/NIR search for high-z galaxies.} We analyse the FIR SEDs of two galaxies proposed as very high redshift ($z>7$) dropout candidates based on deep Hawk-I/VLT observations. The FIR SEDs are sampled with PACS/Herschel at 100 and 160\,$\mu$m, with SPIRE/Herschel at 250, 350 and 500\,$\mu$m and with LABOCA/APEX at 870\,$\mu$m. We find that redshifts $>7$ would imply extreme FIR SEDs (with dust temperatures $>100$\,K and FIR luminosities $>10^{13}$\,$L_{\odot}$). At z$\sim$2, instead, the SEDs of both sources would be compatible with that of typical ULIRGs/SMGs. Considering all the data available for these sources from visible to FIR we re-estimate the redshifts and we find $z\sim$1.6--2.5. Due to the strong spectral breaks observed in these galaxies, standard templates from the literature fail to reproduce the visible-near IR part of the SEDs even when additional extinction is included. These sources resemble strongly dust obscured galaxies selected in Spitzer observations with extreme visible-to-FIR colors, and the galaxy GN10 at $z=4$. Galaxies with similar SEDs could contaminate other high redshift surveys., Comment: Accepted for publication in A&A