Search

Support by grant 2020750 from the United States- Israel Binational Science Foundation (BSF) and grant 2109066 from the United States National Science Foundation (NSF), and by the Ministry of Science and Technology, Israel. J.C. acknowl- edges funding from the ‘FirstGalaxies’ Advanced Grant from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (Grant agreement No. 789056). E.C.L. acknowledges support of an STFC Webb Fellowship (ST/W001438/1). K.K. acknowledges the support by JSPS KAKENHI Grant Number JP17H06130 and the NAOJ ALMA Scientific Research Grant Number 2017–06B. D.E. acknowledges support from a Beatriz Galindo senior fellowship (BG20/00224) from the Spanish Ministry of Science and Innovation, projects PID2020-114414GB-100 and PID2020-113689GB-I00 financed by MCIN/AEI/10.13039/501100011033, project P20-00334 financed by the Junta de Andaluc´ıa, and project A-FQM-510-UGR20 of the FEDER/Junta de Andaluc´ıa-Consejer´ıa de Transformaci ´on Econ ´omica, Industria, Conocimiento y Universidades. G.E.M. ac- knowledges financial support from the Villum Young Investiga- tor grant 37440 and 13160 and the The Cosmic Dawn Center (DAWN), funded by the Danish National Research Foundation under grant No. 140. F.E.B. acknowledges support from ANID- Chile BASAL CATA FB210003, FONDECYT Regular 1200495 and 1190818, and Millennium Science Initiative Program–ICN12 009. K.K.K. acknowledges support from the Knut and Alice Wallenberg Foundation., We report the discovery of a triply imaged active galactic nucleus (AGN), lensed by the galaxy cluster MACS J0035.4−2015 (z d = 0.352). The object is detected in Hubble Space Telescope imaging taken for the RELICS program. It appears to have a quasi-stellar nucleus consistent with a point-source, with a de-magnified radius of re ≲ 100 pc. The object is spectroscopically confirmed to be an AGN at z spec = 2.063 ± 0.005 showing broad rest-frame UV emission lines, and detected in both X-ray observations with Chandra and in ALCS ALMA band 6 (1.2 mm) imaging. It has a relatively faint rest-frame UV luminosity for a quasar-like object, MUV, 1450 = −19.7 ± 0.2. The object adds to just a few quasars or other X-ray sources known to be multiply lensed by a galaxy cluster. Some diffuse emission from the host galaxy is faintly seen around the nucleus, and there is a faint object nearby sharing the same multiple-imaging symmetry and geometric redshift, possibly an interacting galaxy or a star-forming knot in the host. We present an accompanying lens model, calculate the magnifications and time delays, and infer the physical properties of the source. We find the rest-frame UV continuum and emission lines to be dominated by the AGN, and the optical emission to be dominated by the host galaxy of modest stellar mass M✶ ≃ 109.2 M⊙. We also observe some variation in the AGN emission with time, which may suggest that the AGN used to be more active. This object adds a low-redshift counterpart to several relatively faint AGN recently uncovered at high redshifts with HST and JWST., Horizon 2020 Framework Programme H2020, Consejería de Transformación Económica, Industria, Conocimiento y Universidades 13160, 37440, Science and Technology Facilities Council ST/W001438/1 STFC, European Research Council ERC, Fondo Nacional de Desarrollo Científico y Tecnológico 1190818, 1200495, ICN12_009 FONDECYT, Ministerio de Ciencia e Innovación MCIN/AEI/10.13039/501100011033, P20-00334, PID2020-113689GB-I00, PID2020-114414GB-100 MICINN, Horizon 2020 789056, European Regional Development Fund ERDF, Junta de Andalucía A-FQM-510-UGR20

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 detailed spectroscopic analysis of a galaxy at z ≃ 4.88 that is, by chance, magnified ∼30× by gravitational lensing. Only three sources at z ≳ 5 are known with such high magnification. This particular source has been shown to exhibit widespread, high equivalent width ${{\rm C\, \small {IV}}}\ \lambda 1549\, \mathring{\rm A}$ emission, implying it is a unique example of a metal-poor galaxy with a hard radiation field, likely representing the galaxy population responsible for cosmic reionization. Using ultraviolet (UV) nebular line ratio diagnostics, Very Large Telescope (VLT)/X-shooter observations rule out strong active galactic nuclei (AGN) activity, indicating a stellar origin of the hard radiation field instead. We present a new detection of ${[{\rm Ne\, \small {III}}]}\ \lambda 3870\, \mathring{\rm A}$ and use the [${\rm Ne\, \small {III}}$]/[${\rm O\, \small {II}}$] line ratio to constrain the ionization parameter and gas-phase metallicity. Closely related to the commonly used [${\rm O\, \small {III}}$]/[${\rm O\, \small {II}}$] ratio, our [${\rm Ne\, \small {III}}$]/[${\rm O\, \small {II}}$] measurement shows this source is similar to local ‘Green Pea’ galaxies and Lyman-continuum leakers. It furthermore suggests this galaxy is more metal poor than expected from the fundamental metallicity relation, possibly as a consequence of excess gas accretion diluting the metallicity. Finally, we present the highest redshift detection of ${{\rm Mg\, \small {II}}}\ \lambda 2796\, \mathring{\rm A}$, observed at high equivalent width in emission, in contrast to more evolved systems predominantly exhibiting ${\rm Mg\, \small {II}}$ absorption. Strong ${\rm Mg\, \small {II}}$ emission has been observed in most z ∼ 0 Lyman-continuum leakers known and has recently been proposed as an indirect tracer of escaping ionizing radiation. In conclusion, this strongly lensed galaxy, observed just $300\, \mathrm{Myr}$ after reionization ends, enables testing of observational diagnostics proposed to constrain the physical properties of distant galaxies in the James Webb Space Telescope (JWST)/Extremely Large Telescope (ELT) era.

ALMA observations of z > 6 galaxies hav e rev ealed abnormally high [O III ] 88 μm /[C II ] 158 μm ratios and [C II ] 158 μm deficits compared to local galaxies. The origin of this behaviour is unknown. Numerous solutions have been proposed including differences in C and O abundance ratios, observational bias, and differences in ISM properties, including ionization parameter, gas density, or photodissociation region (PDR) covering fraction. In order to elucidate the underlying physics that drives this high- redshift phenomenon, we employ SPHINX 20 , a state-of-the-art, cosmological radiation-hydrodynamics simulation, that resolves detailed ISM properties of thousands of galaxies in the epoch of reionization which has been post-processed with CLOUDY to predict emission lines. We find that the observed z > 6 [O III ] 88 μm -SFR and [C II ] 158 μm -SFR relations can only be reproduced when the C/O abundance ratio is ∼8 ×lower than Solar and the total metal production is ∼4 ×higher than that of a Kroupa IMF. This implies that high-redshift galaxies are potentially primarily enriched by low-metallicity core-collapse supernovae with a more top-heavy IMF. As AGB stars and type-Ia supernova begin to contribute to the galaxy metallicity, both the [C II ] 158 μm -SFR and [C II ] 158 μm luminosity functions are predicted to converge to observed values at z ∼4.5. While we demonstrate that ionization parameter, LyC escape fraction, ISM gas density, and CMB attenuation all drive galaxies towards higher [O III ] 88 μm /[C II ] 158 μm , observed values at z > 6 can only be reproduced with substantially lower C/O abundances compared to Solar. The combination of [C II ] 158 μm and [O III ] 88 μm can be used to predict the values of ionization parameter, ISM gas density, and LyC escape fraction and we provide estimates of these quantities for nine observed z > 6 galaxies. Finally, we demonstrate that [O I ] 63 μm can be used as a replacement for [C II ] 158 μm in high-redshift galaxies where [C II ] 158 μm is unobserved and argue that more observation time should be used to target [O I ] 63 μm at z > 6. Future simulations will be needed to self-consistently address the numerous uncertainties surrounding a varying IMF at high redshift and the associated metal returns.

The 100% detection rate of Lyα emission in a sample of four luminous z ∼ 8 galaxies with red Spitzer/IRAC colors suggests objects with unusual ionizing capabilities that created early ionized bubbles in a neutral era. Whether such bubbles reflect enhanced ionizing properties (nature) or an overdense environment (nurture), however, remains unclear. Here we aim to distinguish between these hypotheses via a search for Lyα emission in five fainter galaxies drawn from the CANDELS-GOODS fields using a similar IRAC excess and UV magnitudes that should reflect reduced clustering effects. Using Keck/MOSFIRE we tentatively detect >4σ line emission in only two targets at redshifts z Lyα = 7.1081 and 7.9622 with rest-frame EWs of 16–17 Å, ∼1.5× weaker compared to their brighter counterparts. Thus, we find a reduced rate for Lyα emission of 0.40 − 0.25 + 0.30 compared to 1.00 − 0.44 + 0.00 for more luminous examples. The lower rate agrees with predictions from simulations of a mostly neutral intergalactic medium and an intrinsic EW0,Lyα distribution for z ∼ 6 galaxies. However, even with an extreme EW0,Lyα model, it is challenging to match the detection rate for the luminous objects. Spectral energy distribution fitting of our fainter sample indicates young and star-forming systems, albeit with less extreme star formation rates and ionization parameters compared to their luminous counterparts. The enhanced Lyα rate in luminous galaxies is thus likely a byproduct of both extreme ionizing properties as well as environmental effects. Further studies with JWST may be required to resolve the physical nature of this puzzling population.

We present new observations with the Atacama Large Millimeter/submillimeter Array for a gravitationally-lensed galaxy at $z=9.1$, MACS1149-JD1. [O III] 88-$\mu$m emission is detected at 10$\sigma$ with a spatial resolution of $\sim0.3$ kpc in the source plane, enabling the most distant morpho-kinematic study of a galaxy. The [O III] emission is distributed smoothly without any resolved clumps and shows a clear velocity gradient with $\Delta V_{\rm obs}/2\sigma_{\rm tot}=0.84\pm0.23$, where $\Delta V_{\rm obs}$ is the observed maximum velocity difference and $\sigma_{\rm tot}$ is the velocity dispersion measured in the spatially-integrated line profile, suggesting a rotating system. Assuming a geometrically thin self-gravitating rotation disk model, we obtain $V_{\rm rot}/\sigma_V=0.67_{-0.26}^{+0.73}$, where $V_{\rm rot}$ and $\sigma_V$ are the rotation velocity and velocity dispersion, respectively, still consistent with rotation. The resulting disk mass of $0.65_{-0.40}^{+1.37}\times10^{9}$ M$_\odot$ is consistent with being associated with the stellar mass identified with a 300 Myr-old stellar population independently indicated by a Balmer break in the spectral energy distribution. We conclude that the most of the dynamical mass is associated with the previously-identified mature stellar population that formed at $z\sim15$., Comment: ApJL accepted

We present improved results of the measurement of the correlation between galaxies and the intergalactic medium transmission at the end of reionization. We have gathered a sample of 13 spectroscopically confirmed Lyman-break galaxies (LBGs) and 21 Lyman-α emitters (LAEs) at angular separations 20 arcsec ≲ θ ≲ 10 arcmin (∼0.1–4 pMpc at z ∼ 6) from the sightlines to eight background z ≳ 6 quasars. We report for the first time the detection of an excess of Lyman-α transmission spikes at ∼10–60 cMpc from LAEs (3.2σ) and LBGs (1.9σ). We interpret the data with an improved model of the galaxy–Lyman-α transmission and two-point cross-correlations, which includes the enhanced photoionization due to clustered faint sources, enhanced gas densities around the central bright objects and spatial variations of the mean free path. The observed LAE(LBG)–Lyman-α transmission spike two-point cross-correlation function (2PCCF) constrains the luminosity-averaged escape fraction of all galaxies contributing to reionization to $\langle f_{\rm esc} \rangle _{M_{\rm UV}\lt -12} = 0.14_{-0.05}^{+0.28}\, (0.23_{-0.12}^{+0.46})$. We investigate if the 2PCCF measurement can determine whether bright or faint galaxies are the dominant contributors to reionization. Our results show that a contribution from faint galaxies ($M_{\rm UV} \gt -20 \, (2\sigma)$) is necessary to reproduce the observed 2PCCF and that reionization might be driven by different sub-populations around LBGs and LAEs at z ∼ 6.

Context. The properties of galaxies at redshift z > 6 hold the key to our understanding of the early stages of galaxy evolution and can potentially identify the sources of the ultraviolet radiation that give rise to the epoch of reionisation. The far-infrared cooling line of [C II] at 158 μm is known to be bright and correlate with the star formation rate (SFR) of low-redshift galaxies, and hence is also suggested to be an important tracer of star formation and interstellar medium properties for very high-redshift galaxies. Aims. With the aim to study the interstellar medium properties of gravitationally lensed galaxies at z > 6, we search for [C II] and thermal dust emission in a sample of 52 z ∼ 6 galaxies observed by the ALMA Lensing Cluster Survey. Methods. We perform our analysis using LINESTACKER, stacking both [C II] and continuum emission. The target sample is selected from multiple catalogues, and the sample galaxies have spectroscopic redshift or low-uncertainty photometric redshifts (σz Results. Our analyses find no detection of either [C II] or continuum. An upper limit on L[CII] is derived, implying that [C II] remains marginally consistent for low-SFR z > 6 galaxies but likely is under-luminous compared to the local L[CII]-SFR relationship. We discuss potential biases and possible physical effects that may be the cause of the non-detection. Further, the upper limit on the dust continuum implies that less than half of the star formation is obscured.

A fundamental quest of modern astronomy is to locate the earliest galaxies and study how they influenced the intergalactic medium a few hundred million years after the Big Bang. The abundance of star-forming galaxies is known to decline from redshifts of about 6 to 10, but a key question is the extent of star formation at even earlier times, corresponding to the period when the first galaxies might have emerged. Here we present spectroscopic observations of MACS1149-JD1, a gravitationally lensed galaxy observed when the Universe was less than four per cent of its present age. We detect an emission line of doubly ionized oxygen at a redshift of $9.1096\pm0.0006$, with an uncertainty of one standard deviation. This precisely determined redshift indicates that the red rest-frame optical colour arises from a dominant stellar component that formed about 250 million years after the Big Bang, corresponding to a redshift of about 15. Our results indicate the it may be possible to detect such early episodes of star formation in similar galaxies with future telescopes., Comment: To appear in the 17 May issue of Nature

Infrared and nebular lines provide some of our best probes of the physics regulating the properties of the interstellar medium (ISM) at high-redshift. However, interpreting the physical conditions of high-redshift galaxies directly from emission lines remains complicated due to inhomogeneities in temperature, density, metallicity, ionisation parameter, and spectral hardness. We present a new suite of cosmological, radiation-hydrodynamics simulations, each centred on a massive Lyman-break galaxy that resolves such properties in an inhomogeneous ISM. Many of the simulated systems exhibit transient but well defined gaseous disks that appear as velocity gradients in [CII]~158.6$\mu$m emission. Spatial and spectral offsets between [CII]~158.6$\mu$m and [OIII]~88.33$\mu$m are common, but not ubiquitous, as each line probes a different phase of the ISM. These systems fall on the local [CII]-SFR relation, consistent with newer observations that question previously observed [CII]~158.6$\mu$m deficits. Our galaxies are consistent with the nebular line properties of observed $z\sim2-3$ galaxies and reproduce offsets on the BPT and mass-excitation diagrams compared to local galaxies due to higher star formation rate (SFR), excitation, and specific-SFR, as well as harder spectra from young, metal-poor binaries. We predict that local calibrations between H$\alpha$ and [OII]~3727$\AA$ luminosity and galaxy SFR apply up to $z>10$, as do the local relations between certain strong line diagnostics (R23 and [OIII]~5007$\AA$/H$\beta$) and galaxy metallicity. Our new simulations are well suited to interpret the observations of line emission from current (ALMA and HST) and upcoming facilities (JWST and ngVLA)., Comment: 21 pages, 16 figures, MNRAS accepted

We report the discovery of the first spectroscopically resolved C II /C II* 1334, 1335A doublet in the Lyman-break galaxy J0215-0555 at z = 5.754. The separation of the resonant and fluorescent emission channels was possible thanks to the large redshift of the source and long integration time, as well as the small velocity width of the feature, 0.6 +- 0.2A. We model this emission and find that at least two components are required to reproduce the combination of morphologies of C II* emission, C II absorption and emission, and Lyman-alpha emission from the object. We suggest that the close alignment between the fluorescence and Lyman-alpha emission could indicate an ionisation escape channel within the object. While the faintness of such a C II /C II* doublet makes it prohibitively difficult to pursue for similar systems with current facilities, we suggest it can become a valuable porosity diagnostic in the era of JWST and the upcoming generations of ELTs., 5 pages, 3 figures; accepted for publication in MNRAS Letters

The presence of spectroscopically confirmed Balmer breaks in galaxy spectral energy distributions (SEDs) at $z>9$ provides one of the best probes of the assembly history of the first generations of stars in our Universe. Recent observations of the gravitationally lensed source, MACS 1149_JD1 (JD1), indicate that significant amounts of star formation likely occurred at redshifts as high as $z\simeq15$. The inferred stellar mass, dust mass, and assembly history of JD1, or any other galaxy at these redshifts that exhibits a strong Balmer break, can provide a strong test of our best theoretical models from high-resolution cosmological simulations. In this work, we present the results from a cosmological radiation-hydrodynamics simulation of the region surrounding a massive Lyman-break galaxy. For two of our most massive systems, we show that dust preferentially resides in the vicinity of the young stars thereby increasing the strength of the measured Balmer break such that the simulated SEDs are consistent with the photometry of JD1 and two other $z>9$ systems (GN-z10-3 and GN-z9-1) that have proposed Balmer breaks at high redshift. We find strong variations in the shape and luminosity of the SEDs of galaxies with nearly identical stellar and halo masses, indicating the importance of morphology, assembly history, and dust distribution in making inferences on the properties of individual galaxies at high redshifts. Our results stress the importance that dust may play in modulating the observable properties of galaxies, even at the extreme redshifts of $z>9$., 16 pages, 13 Figures, Accepted 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 present a study of the intersection between the populations of star forming galaxies selected as either Lyman break galaxies (LBGs) or Lyman-alpha emitters (LAEs) in the redshift range 2.9 − 6.7 and 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 an LAE counterpart (408 galaxies), and LAE galaxies without an LBG counterpart (46 galaxies). All these galaxies are intrinsically faint because of the lensing nature of the sample (M1500 ≥ −20.5). The fraction of LAEs among all selected star-forming galaxies increases with redshift up to z ∼ 6 and decreases for higher redshifts, in agreement with previous findings. The evolution of LAE/LBG populations with UV magnitude and Lyα luminosity shows that the LAE selection is able to identify intrinsically UV faint galaxies with M1500 ≥ −15 that are typically missed in the deepest lensing photometric surveys. The LBG population seems to fairly represent the total population of star-forming galaxies down to M1500 ∼ −15. Galaxies with M1500 Lyα uv, whereas the opposite trend is observed within our sample for faint galaxies with M1500 > −17, including galaxies only detected by their Lyα emission, with a large scatter. These trends, previously observed in other samples of star-forming galaxies at high-z, are seen here for very faint M1500 ∼ −15 galaxies; that is, much fainter than in previous studies. The present results show no clear evidence for an intrinsic difference between the properties of the two populations selected as LBG and/or LAE. The observed trends could be explained by a combination of several phenomena, like the existence of different star-formation regimes, the dust content, the relative distribution and morphology of dust and stars, or the stellar populations.

We examine the absorption line spectra of a sample of 31 luminous (M_UV=-23) Lyman break galaxies at redshift z~6 using data taken with the FOCAS and OSIRIS spectrographs on the Subaru and GTC telescopes. For two of these sources we present longer exposure data taken at higher spectral resolution from ESO's X-shooter spectrograph. Using these data, we demonstrate the practicality of stacking our lower resolution data to measure the depth of various interstellar and stellar absorption lines to probe the covering fraction of low ionization gas and the gas-phase and stellar metallicities near the end of the era of cosmic reionization. From maximum absorption line depths of SiII1260 and CII1334, we infer a mean covering fraction of >0.85+/-0.16 for our sample. This is larger than that determined using similar methods for lower luminosity galaxies at slightly lower redshifts, suggesting that the most luminous galaxies appear to have a lower escape fraction than fainter galaxies, and therefore may not play a prominent role in concluding reionization. Using various interstellar absorption lines we deduce gas-phase metallicities close to solar indicative of substantial early enrichment. Using selected stellar absorption lines, we model our spectra with a range of metallicities using techniques successfully employed at lower redshift and deduce a stellar metallicity of 0.4 +0.3/-0.1 solar, consistent with the stellar mass - stellar metallicity relation recently found at z~3-5. We discuss the implications of these metallicity estimates for the typical ages of our luminous galaxies and conclude our results imply initial star formation at redshifts z~10, consistent with independent analyses of earlier objects., 13 pages, 10 figures, 3 tables, Accepted for Publication in ApJ

Identifying low-redshift galaxies that emit Lyman Continuum radiation (LyC leakers) is one of the primary, indirect methods of studying galaxy formation in the epoch of reionization. However, not only has it proved challenging to identify such systems, it also remains uncertain whether the low-redshift LyC leakers are truly "analogues" of the sources that reionized the Universe. Here, we use high-resolution cosmological radiation hydrodynamics simulations to examine whether simulated galaxies in the epoch of reionization share similar emission line properties to observed LyC leakers at $z\sim3$ and $z\sim0$. We find that the simulated galaxies with high LyC escape fractions ($f_{\rm esc}$) often exhibit high O32 and populate the same regions of the R23-O32 plane as $z\sim3$ LyC leakers. However, we show that viewing angle, metallicity, and ionisation parameter can all impact where a galaxy resides on the O32-$f_{\rm esc}$ plane. Based on emission line diagnostics and how they correlate with $f_{\rm esc}$, lower-metallicity LyC leakers at $z\sim3$ appear to be good analogues of reionization-era galaxies. In contrast, local [SII]-deficient galaxies do not overlap with the simulated high-redshift LyC leakers on the SII-BPT diagram; however, this diagnostic may still be useful for identifying leakers. We use our simulated galaxies to develop multiple new diagnostics to identify LyC leakers using IR and nebular emission lines. We show that our model using only [CII]$_{\rm 158\mu m}$ and [OIII]$_{\rm 88\mu m}$ can identify potential leakers from non-leakers from the local Dwarf Galaxy Survey. Finally, we apply this diagnostic to known high-redshift galaxies and find that MACS1149_JD1 at $z=9.1$ is the most likely galaxy to be actively contributing to the reionization of the Universe., Comment: 18 pages, 11 figures, Accepted to MNRAS

We introduce a new method for determining the influence of galaxies and active galactic nuclei (AGN) on the physical state of the intergalactic medium (IGM) at high redshift and illustrate its potential via a first application to the field of the $z=6.42$ QSO J1148+5251. By correlating the spatial positions of spectroscopically-confirmed Lyman break galaxies (LBGs) with fluctuations in the Lyman alpha forest seen in the high signal-to-noise spectrum of a background QSO, we provide a statistical measure of the typical escape fraction of Lyman continuum photons close to the end of cosmic reionisation. Here we use Keck DEIMOS spectroscopy to locate 7 colour-selected LBGs in the redshift range $5.3\lesssim z\lesssim 6.4$ and confirm a faint $z=5.701$ AGN. We then examine the spatial correlation between this sample and Ly$\alpha$/Ly$\beta$ transmission fluctuations in a Keck ESI spectrum of the QSO. Interpreting the statistical HI proximity effect as arising from faint galaxies clustered around the detected LBGs, we translate the observed mean Ly$\alpha$ transmitted flux around an average detected LBG into a constraint on the mean escape fraction $\langle f_{\rm esc}\rangle\geq0.08$ at $z\simeq6$. We also report evidence of the individual transverse HI proximity effect of a $z=6.177$ luminous LBG via a Ly$\beta$ transmission spike and two broad Ly$\alpha$ transmission spikes around the $z=5.701$ AGN. We discuss the possible origin of such associations which suggest that while faint galaxies are primarily driving reionisation, luminous galaxies and AGN may provide important contributions to the UV background or thermal fluctuations of the IGM at $z\simeq6$. Although a limited sample, our results demonstrate the potential of making progress using this method in resolving one of the most challenging aspects of the contribution of galaxies and AGN to cosmic reionisation., Comment: 21 pages, 16 figures, the version accepted in MNRAS

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., Comment: 21 pages, plus 13 pages of additional figures and tables. 13 figures, 4 tables. Accepted for publication in MNRAS

Physical properties of galaxies at z>7 are of interest for understanding both the early phases of star formation and the process of cosmic reionization. Chemical abundance measurements offer valuable information on the integrated star formation history, and hence ionizing photon production, as well as the rapid gas accretion expected at such high redshifts. We use reported measurements of [O III] 88$��$m emission and star formation rate to estimate gas-phase oxygen abundances in five galaxies at z=7.1-9.1 using the direct T_e method. We find typical abundances 12+log(O/H) = 7.9 ($\sim$0.2 times the solar value) and an evolution of 0.9$\pm$0.5 dex in oxygen abundance at fixed stellar mass from z$\simeq$8 to 0. These results are compatible with theoretical predictions, albeit with large (conservative) uncertainties in both mass and metallicity. We assess both statistical and systematic uncertainties to identify promising means of improvement with the Atacama Large Millimeter Array (ALMA) and the James Webb Space Telescope (JWST). In particular we highlight [O III] 52$��$m as a valuable feature for robust metallicity measurements. Precision of 0.1-0.2 dex in T_e-based O/H abundance can be reasonably achieved for galaxies at z$\approx$5-8 by combining [O III] 52$��$m with rest-frame optical strong lines. It will also be possible to probe gas mixing and mergers via resolved T_e-based abundances on kpc scales. With ALMA and JWST, direct metallicity measurements will thus be remarkably accessible in the reionization epoch., 12 pages, 3 figures, 4 tables. Accepted for publication in ApJ

Quasar proximity zones at z > 5.5 correspond to overdense and overionized environments. Galaxies found inside proximity zones can therefore display features that would otherwise be masked by absorption in the intergalactic medium. We demonstrate the utility of this quasar-galaxy synergy by reporting the discovery of the first three “proximate Lyα emitters” (LAEs) within the proximity zone of quasar J0836+0054 at z = 5.795 (Aerith A, B, and C). Aerith A, located behind the quasar with an impact parameter D^{\perp} = 278 \pm 8 pkpc, provides the first detection of an Lyα transverse proximity effect. We model the transmission and show that it constrains the onset of J0836ʼs quasar phase to 0.2Myr < 28Myr < t in the past. The second object, Aerith B at a distance D < 912 pkpc from the quasar, displays a bright and broad double-peaked Lyα emission line. The peak separation implies a low ionizing f_{esc} \leqslant 1%. We fit the Lyα line with an outflowing shell model, finding a typical central density log N_{HI}/cm^{-2} = 19.3_{-0.2}^{+0.8}, outflow velocity v_{out} = 16_{-11}^{+4}km s^{-1}, and gas temperature log T/K = 3.8_{-0.7}^{+0.8} compared to 2 < z < 3 analog LAEs. We detect object Aerith C via an Lyα emission line at z = 5.726. This corresponds with the edge of the quasar’s proximity zone (Dz < 0.02), suggesting that the proximity zone is truncated by a density fluctuation. Via the analyses conducted here, we illustrate how proximate LAEs offer unique insight into the ionizing properties of both quasars and galaxies during hydrogen reionization.

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), Accepted for publication in Astronomy & Astrophysics. 17 pages, 11 figures

Recent studies have revealed intense UV metal emission lines in a modest sample of z>7 Lyman-alpha emitters, indicating a hard ionizing spectrum is present. If such high ionization features are shown to be common, it may indicate that extreme radiation fields play a role in regulating the visibility of Lyman-alpha in the reionization era. Here we present deep near-infrared spectra of seven galaxies with Lyman-alpha emission at 5.47 photometric targets. In nine sources we do not detect UV metal lines. However in the z=8.683 galaxy EGSY8p7, we detect a 4.6 sigma emission line in the narrow spectral window expected for NV 1243. The feature is unresolved (FWHM7 have now been shown to have intense UV line emission, suggesting that extreme radiation fields are commonplace among the Lyman-alpha population. Future observations with JWST will eventually clarify the origin of these features and explain their role in the visibility of Lyman-alpha in the reionization era., Comment: 15 pages, 3 tables and 7 figures. Accepted for publication in MNRAS

The ubiquity of Lyman alpha (Ly$\alpha$) emission in a sample of four bright [O III]-strong star-forming galaxies with redshifts above 7 has led to the suggestion that such luminous sources represent a distinct population compared to their fainter, more numerous, counterparts. The presence of Ly$\alpha$ emission within the reionization era could indicate that these sources created early ionized bubbles due to their unusually strong radiation, possibly due to the presence of active galactic nuclei. To test this hypothesis we have secured long integration spectra with XSHOOTER on the VLT for three $z\simeq$7 sources selected to have similar luminosities and prominent excess fluxes in the IRAC 3.6 or 4.5$\mu$m band, usually attributed to strong [O III] emission. We secured additional spectroscopy for one of these galaxies at $z$=7.15 using MOSFIRE at the Keck telescope. For this, the most well-studied source in our sample with the strongest IRAC excess, we detect prominent nebular emission from He II and NV indicative of a non-thermal source. For the other two sources at $z$=6.81 and $z$=6.85, for which no previous optical/near infrared spectroscopy was initially available, Ly$\alpha$ is seen in one and CIII] emission in the other. Although a modest sample, our results further support the hypothesis that the phenomenon of intense [O III] emission is associated preferentially with sources lying in early ionized bubbles. However, even though one of our sources at $z$=7.15 clearly indicates the presence of non-thermal radiation, such ionized bubbles may not uniquely arise in this manner. We discuss the unique advantages of extending such challenging diagnostic studies with JWST., Comment: 11 pages, 7 figures, submitted to 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., new version 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 $$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$ \approx$1.6) and high dust contents (8.1 $, 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., Accepted for publication in ApJ

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., To appear in the Astrophysical Journal. 20 pages, 6 figures, 3 tables. arXiv admin note: text overlap with arXiv:1402.6743

[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

The different algorithms appropriate for point source photometry on data from the SPIRE instrument on-board the Herschel Space Observatory, within the Herschel Interactive Processing Environment (HIPE) are compared. Point source photometry of a large ensemble of standard calibration stars and dark sky observations is carried out using the 4 major methods within HIPE: SUSSEXtractor, DAOphot, the SPIRE Timeline Fitter and simple Aperture Photometry. Colour corrections and effective beam areas as a function of the assumed source spectral index are also included to produce a large number of photometric measurements per individual target, in each of the 3 SPIRE bands (250, 350, 500μm), to examine both the accuracy and repeatability of each of the 4 algorithms. It is concluded that for flux densities down to the level of 30mJy that the SPIRE Timeline Fitter is the method of choice. However, at least in the 250 and 350μm bands, all 4 methods provide photometric repeatability better than a few percent down to at approximately 100mJy. The DAOphot method appears in many cases to have a systematic offset of ∼8 % in all SPIRE bands which may be indicative of a sub-optimal aperture correction. In general, aperture photometry is the least reliable method, i.e. largest scatter between observations, especially in the longest wavelength band. At the faintest fluxes

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 $��$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, 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., 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}$$, Comment: 14 pages, 12 figures, accepted for publication in ApJ

International audience; We report on the detailed analysis of a gravitationally lensed Y-band dropout, A2744_YD4, selected from deep Hubble Space Telescope imaging in the Frontier Field cluster Abell 2744. Band 7 observations with the Atacama Large Millimeter/submillimeter Array (ALMA) indicate the proximate detection of a significant 1 mm continuum flux suggesting the presence of dust for a star-forming galaxy with a photometric redshift of $z\simeq 8$. Deep X-SHOOTER spectra confirms the high-redshift identity of A2744_YD4 via the detection of Lyα emission at a redshift z = 8.38. The association with the ALMA detection is confirmed by the presence of [O iii] 88 μm emission at the same redshift. Although both emission features are only significant at the 4σ level, we argue their joint detection and the positional coincidence with a high-redshift dropout in the Hubble Space Telescope images confirms the physical association. Analysis of the available photometric data and the modest gravitational magnification ($\mu \simeq 2$) indicates A2744_YD4 has a stellar mass of ∼2 × 10(9) ${M}_{\odot }$, a star formation rate of ∼20 ${M}_{\odot }$ yr(−)(1) and a dust mass of ∼6 × 10(6) ${M}_{\odot }$. We discuss the implications of the formation of such a dust mass only $\simeq 200\,\mathrm{Myr}$ after the onset of cosmic reionization.

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., 6 pages; 4 colour figures

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, Submitted to ApJ, 10 pages, 6 figures

We present a method for selecting $z>4$ dusty, star forming galaxies (DSFGs) using Herschel/SPIRE 250/350/500 $��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 $��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\,��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., Accepted for publication in ApJ

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)-19.7 with fixed alpha=-1.74 and Phi*=1.10x10^{-3}Mpc^{-3}. The corresponding star formation rate density should be rho_{SFR}, 4 pages, 2 figures, accepted for publication in A&A Letters

arXiv:1510.07084v1.-- Infante, Leopoldo et al., We searched for z ≲ 7 Lyman-break galaxies 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 6 lie at z ≲ 9 and 1 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, . This substantial increase in the number of known high-redshift galaxies allows a solid determination of the luminosity function (LF) 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 M, star formation rates (SFRs) of M yr, and SFR-weighted ages of . Finally, we are able to put strong constraints on the z = 7, 8, 9, and 10 LFs. One of the objects in the cluster field is a z ≃ 10 candidate, with a magnification of μ ∼ 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.04 L∗) of the z ∼ 10 LF (It is named >Tayna> in the Aymara language)., We acknowledge support from Basal-CATA PFB-06/2007 (LI, FEB, NL, SK, PT), FP7-SPACE-2012-ASTRODEEP-312725 (XS), NSFC grants 11103017 and 11233002 (XS), NASA through Hubble Fellowship grant HST-HF2-51334.01-A awarded by STScI (AZ), CONICYT-Chile grants FONDECYT 1141218 (FEB), gemini-conicyt 32120003 (FEB, NL), and “EMBIGGEN” Anillo ACT1101 (FEB), the Ministry of Economy, Development, and Tourism’s Millennium Science Initiative through grant IC120009, awarded to The Millennium Institute of Astrophysics, MAS (FEB), Spanish consolider project CAD2010-00064 and AYA2012-39475-C02-01 (JMD), the Brazilian funding agency-FAPESP pos-doc fellowship - process number 2014/11806-9 (AM) and CONICYT-Chile grant FONDECYT 3140542 (PT).

We investigate the nature and the physical properties of z, Y and J-dropout galaxies selected behind the lensing cluster A2667. This field is part of our project aimed at identifying z~7-10 candidates accessible to spectroscopic studies, based on deep photometry with ESO/VLT HAWK-I and FORS2 (zYJH and Ks-band images, AB(3 sigma)~26-27) on a sample of lensing clusters extracted from our multi-wavelength combined surveys with SPITZER, HST, and Herschel. In this paper we focus on the complete Y and J-dropout sample, as well as the bright z-dropouts fulfilling the selection criteria by Capak et al. (2011). 10 candidates are selected within the common field of ~33 arcmin2 (effective area once corrected for contamination and lensing dilution). All of them are detected in H and Ks bands in addition to J and/or IRAC 3.6/4.5, with H(AB)~23.4 to 25.2, and have modest magnification factors. Although best-fit photometric redshifts place all these candidates at high-z, the contamination by low-z interlopers is estimated at 50-75% level based on previous studies, and the comparison with the blank-field WIRCAM Ultra-Deep Survey (WUDS). The same result is obtained when photometric redshifts include a luminosity prior, allowing us to remove half of the original sample as likely z~1.7-3 interlopers with young stellar pulations and strong extinction. Two additional sources among the remaining sample could be identified at low-z based on a detection at 24 microns and on the HST z_850 band. These low-z interlopers are not well described by current templates given the large break, and cannot be easily identified based solely on optical and near-IR photometry. Given the estimated dust extinction and high SFRs, some of them could be also detected in the IR or sub-mm bands. After correction for likely contaminants, the observed counts at z>7.5 seem to be in agreement with an evolving LF. (abridged), 18 pages, 11 figures. Accepted for publication in A&A

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\,$��$m, with SPIRE/Herschel at 250, 350 and 500\,$��$m and with LABOCA/APEX at 870\,$��$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., Accepted for publication in A&A

Context. The Hubble Space Telescope (HST) Frontier Fields (HFFs) project started at the end of 2013 with the aim of providing extremely deep images of six massive galaxy clusters. One of the main goals of this program is to push several telescopes to their limits to provide the best current view of the earliest stages of the Universe. The analysis of the initial data has already demonstrated the huge capabilities of the program.Aims. We present a detailed analysis of z ~ 8 objects behind the HFFs lensing cluster, MACSJ0416.1-2403, combining 0.3−1.6 μ m imaging from HST, ground-based K s imaging from VLT HAWK-I, and 3.6 μ m and 4.5 μ m Spitzer Space Telescope. The images probe to 5σ depths of ≈29 AB for HST, 25.6 AB for HAWK-I, and ≈0.310 and 0.391 μ Jy at 3.6 and 4.5 μ m, respectively. With these datasets, we assess the photometric properties of z ~ 8 galaxies in this field, as well as their distribution in luminosity, to unprecedented sensitivity.Methods. We applied the classical Lyman break (LB) technique, which combines non detection criteria in bands blueward of the Lyman break at z ~ 8 and color-selection in bands redward of the break. To avoid contamination by mid-z interlopers, we required a strong break between optical and near-infrared data. We determined the photometric properties of z ~ 8 selected candidates using spectral energy distribution (SED)-fitting with standard library templates. The luminosity function at z ~ 8 is computed using a Monte-Carlo 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 ~ 8 galaxy in this cluster, which is now firmly detected as a mid-z interloper with a strong ≈1.5 mag Balmer break (between F606W and F125W ). Using the SED of this interloper, we estimated the contamination rate of our MACSJ0416.1−2403 sample, and that of previous samples in Abell 2744 that were based on HFF data, we highlight the dangers of pushing the LB technique too close to the photometry limits.Results. Our selection reliably recovers four objects with m F 160 W ranging from 26.0 to 27.9 AB that are located in modest-amplification regions (μ m and 4.5 μ m bands, which could be associated to the Balmer break or emission lines at z ~ 8. The SED-fitting analysis suggests that all of these objects favor high-z solutions with no reliable secondary solutions. The candidates generally have star formation rates around ~10 M ⊙ /yr and sizes ranging from 0.2 to 0.5 kpc, which agrees well with previous observations and expectations for objects in the early Universe. The sample size and luminosity distribution are consistent with previous findings.

We have used deep optical, near-IR, and IR observations from VLT, Spitzer, Herschel, and LABOCA in strong lensing clusters to study distant galaxies. In searches for optical-dropout galaxies (i.e. for z ≳ 7 candidates) we have found several galaxies with very unusual SEDs characterised by a strong spectral break, presumably indicative of high-z, although the objects are detected even in the Herschel bands between 160 and 500 μm and at 870 μm. The latter indicates, from simple estimates of the bolometric luminosity and from the IR SED, that these ob jects are most likely at z ~ 2–2.5.The resulting SEDs imply very high IR/UV ratios, indicative of very large attenuation. Despite this, the large spectral break observed between the optical and near-IR data is difficult to understand with currently know spectral templates from galaxies, EROs, SMGs, and others, both empirical and theoretical ones.

The Hubble Frontier Fields (HFF) program combines the capabilities of the Hubble Space Telescope (HST) with the gravitational lensing of massive galaxy clusters to probe the distant Universe to an unprecedented depth. Here, we present the results of the first combined HST and Spitzer observations of the cluster Abell 2744. We combine the full near-infrared data with ancillary optical images to search for gravitationally lensed high-redshift (z > 6) galaxies. We report the detection of 15 I814-dropout candidates at z ~ 6-7 and one Y105-dropout at z ~ 8 in a total survey area of 1.43 arcmin^2 in the source plane. The predictions of our lens model allow us to also identify five multiply-imaged systems lying at redshifts between z ~ 6 and z ~ 8. Thanks to constraints from the mass distribution in the cluster, we were able to estimate the effective survey volume corrected for completeness and magnification effects. This was in turn used to estimate the rest-frame ultraviolet luminosity function (LF) at z ~ 6-8. Our LF results are generally in agreement with the most recent blank field estimates, confirming the feasibility of surveys through lensing clusters. Although based on a shallower observations than what will be achieved in the final dataset including the full ACS observations, the LF presented here extends down to Muv ~ -18.5 at z ~ 7 with one identified object at Muv ~ -15 thanks to the highly-magnified survey areas. This early study forecasts the power of using massive galaxy clusters as cosmic telescopes and its complementarity to blank fields., Comment: 9 pages, 7 figures, 3 tables. accepted for publication in ApJ

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 report the discovery of a double-peaked Lyman-alpha profile in a galaxy at $z=6.803$, A370p\_z1, in the parallel Frontier Field of Abell 370. The velocity separation between the blue and red peaks of the Lyman-$\alpha$ profile ($\Delta v=101_{-19}^{+38} (\pm48)\,\text{km s}^{-1}$) suggests an extremely high escape fraction of ionising photons $> 59(51)\% (2\sigma)$. The spectral energy distribution indicates a young ($50$ Myr), star-forming ($12\pm 6 \, \rm{M}_\odot \rm{yr}^{-1}$) galaxy with an IRAC excess implying strong [OIII]+H$\beta$ emission. On the basis of the high escape fraction measured, we demonstrate that A370p\_z1 was solely capable of creating an ionised bubble sufficiently large to account for the blue component of its Lyman-alpha profile. We discuss whether A370p\_z1 may be representative of a larger population of luminous $z\simeq7$ double-peaked Lyman-alpha emitting sources with high escape fractions that self-ionised their surroundings without contributions from associated UV-fainter sources., Comment: Accepted for publication in MNRAS. 7 pages, 4 figures, 3 tables

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

Contact. This paper presents the results obtained with the Multi-Unit Spectroscopic Explorer (MUSE) at the ESO Very Large Telescope on the faint end of the Lyman-alpha luminosity function (LF) based on deep observations of four lensing clusters. The goal of our project is to set strong constraints on the relative contribution of the Lyman-alpha emitter (LAE) population to cosmic reionization. Aims. The precise aim of the present study is to further constrain the abundance of LAEs by taking advantage of the magnification provided by lensing clusters to build a blindly selected sample of galaxies which is less biased than current blank field samples in redshift and luminosity. By construction, this sample of LAEs is complementary to those built from deep blank fields, whether observed by MUSE or by other facilities, and makes it possible to determine the shape of the LF at fainter levels, as well as its evolution with redshift. Methods. We selected a sample of 156 LAEs with redshifts between 2.9 ≤ z ≤ 6.7 and magnification-corrected luminosities in the range 39 ≲ log LLyα [erg s−1] ≲43. To properly take into account the individual differences in detection conditions between the LAEs when computing the LF, including lensing configurations, and spatial and spectral morphologies, the non-parametric 1/Vmax method was adopted. 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 determine the effective volume Vmax for each galaxy. In this paper we present a complete procedure for the determination of the LF based on IFU detections in lensing clusters. This procedure, including some new methods for masking, effective volume integration and (individual) completeness determinations, has been fully automated when possible, and it can be easily generalized to the analysis of IFU observations in blank fields. Results. As a result of this analysis, the Lyman-alpha LF has been obtained in four different redshift bins: 2.9 < z < 6, 7, 2.9 < z < 4.0, 4.0 < z < 5.0, and 5.0 < z < 6.7 with constraints down to log LLyα = 40.5. From our data only, no significant evolution of LF mean slope can be found. When performing a Schechter analysis also including data from the literature to complete the present sample towards the brightest luminosities, a steep faint end slope was measured varying from α = −1.69−0.08+0.08 to α = −1.87−0.12+0.12 between the lowest and the highest redshift bins. Conclusions. The contribution of the LAE population to the star formation rate density at z ∼ 6 is ≲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. Depending on the intersection between the LAE/LBG populations, the contribution of the observed galaxies to the ionizing flux may suffice to keep the universe ionized at z ∼ 6., Astronomy & Astrophysics, 628, ISSN:0004-6361, ISSN:1432-0746