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$\mu$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$\mu$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$\mu$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.

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 luminous galaxies do 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.

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, but it also remains uncertain whether the low-redshift LyC leakers are true "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, we find that identifying low-redshift galaxies based on [SII] deficiencies does not seem to produce true analogues. 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.

We present improved results of the measurement of the correlation between galaxies and the intergalactic medium (IGM) transmission at the end of reionisation. We have gathered a sample of $13$ spectroscopically confirmed Lyman-break galaxies (LBGs) and $21$ Lyman-$\alpha$ emitters (LAEs) at angular separations $20'' \lesssim \theta \lesssim 10'$ ($\sim 0.1-4$ pMpc at $z\sim 6$) from the sightlines to $8$ background $z\gtrsim 6$ quasars. We report for the first time the detection of an excess of Lyman-$\alpha$ transmission spikes at $\sim 10-60$ cMpc from LAEs ($3.2\sigma$) and LBGs ($1.9\sigma$). We interpret the data with an improved model of the galaxy-Lyman-$\alpha$ transmission and two-point cross-correlations which includes the enhanced photoionisation 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-$\alpha$ transmission spike two-point cross-correlation function (2PCCF) constrains the luminosity-averaged escape fraction of all galaxies contributing to reionisation to $\langle f_{\rm esc} \rangle_{M_{\rm UV}<-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 reionisation. Our results show that a contribution from faint galaxies ($M_{\rm UV} > -20 \, (2\sigma)$) is necessary to reproduce the observed 2PCCF and that reionisation might be driven by different sub-populations around LBGs and LAEs at $z\sim 6$.

Prior to the launch of JWST, Spitzer/IRAC photometry offers the only means of studying the rest-frame optical properties of z>7 galaxies. Many such high redshift galaxies display a red [3.6] - [4.5] micron colour, often referred to as the "IRAC excess", which has conventionally been interpreted as arising from intense [OIII]+Hbeta emission within the [4.5] micron bandpass. An appealing aspect of this interpretation is similarly intense line emission seen in star-forming galaxies at lower redshift as well as the redshift-dependent behaviour of the IRAC colours beyond z~7 modelled as the various nebular lines move through the two bandpasses. In this paper we demonstrate that, given the photometric uncertainties, established stellar populations with Balmer (4000 A, rest-frame) breaks, such as those inferred at z>9 where line emission does not contaminate the IRAC bands, can equally well explain the redshift-dependent behaviour of the IRAC colours in 7<z<9 galaxies. We discuss possible ways of distinguishing between the two hypotheses using ALMA measures of [OIII] 88 micron and dust continuum fluxes. Prior to further studies with JWST, we show that the distinction is important in determining the assembly history of galaxies in the first 500 Myr.

Quasar proximity zones at $z>5.5$ correspond to over-dense and over-ionized environments. Galaxies found inside proximity zones can therefore display features which would otherwise be masked by absorption in the IGM. We demonstrate the utility of this quasar-galaxy synergy by reporting the discovery of the first three `proximate Lyman-$\alpha$ emitters' (LAEs) within the proximity zone of quasar J0836 at $z=5.802$ (\textit{Aerith A, B} and \textit{C}). \textit{Aerith A}, located behind the quasar with an impact parameter $D_\perp = 278$ pkpc, provides the first detection of a Lyman-$\alpha$ transverse proximity effect. We model the transmission and show it constrains the onset of J0836's quasar phase to $0.2 \text{Myr}<t<20\text{Myr}$ in the past. The second object, \textit{Aerith B} at a distance $D=750$ pkpc from the quasar, displays a bright, broad double-peaked \lal emission line. Based on relations calibrated at $z\leq3$, the peak separation implies a low ionizing $f_{\text{esc}} \lesssim 1\%$, the most direct such constraint on a reionization-era galaxy. We fit the Ly-$\alpha$ line with an outflowing shell model, finding a completely typical central density $\text{log N}_{\text{HI}}/\text{cm}^{-2} = 19.3_{-0.2}^{+0.8}$, outflow velocity $v=16_{-11}^{+4}$ km s$^{-1}$, and gas temperature $\text{log} T/\text{K} = 3.8_{-0.7}^{+0.8}$ compared to $2<z<3$ analogue LAEs. Finally, we detect an emission line at $\lambda=8177$ \AA\ in object \textit{Aerith C} which, if it is \lal at $z=5.726$, would correspond closely with the end of the quasar's proximity zone ($\Delta z<0.02$ from the boundary) and suggests the quasar influences the IGM up to $\sim85$ cMpc away, making it the largest quasar proximity zone. Via the analyses conducted here, we illustrate how proximate LAEs offer unique insight into the ionizing properties of both quasars and galaxies during the epoch of reionization.

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)

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).

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.

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.

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.

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$.

[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.

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)

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.4<z<8.7 (including a newly-confirmed lensed galaxy at z=6.031) and three bright z>7 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 (FWHM<90 km/s) and is likely nebular in origin. A deep H-band spectrum of EGSY8p7 reveals non-detections of CIV, He II, and OIII]. The presence of NV requires a substantial flux of photons above 77 eV, pointing to a hard ionizing spectrum powered by an AGN or fast radiative shocks. Regardless of its origin, the intense radiation field of EGSY8p7 may aid the transmission of Lyman-alpha through what is likely a partially neutral IGM. With this new detection, five of thirteen known Lyman-alpha emitters at z>7 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.

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.

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.

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.

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.

[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.

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.

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 Array (ALMA) indicate the proximate detection of a significant 1mm continuum flux suggesting the presence of dust for a star-forming galaxy with a photometric redshift of $z\simeq8$. Deep X-SHOOTER spectra confirms the high redshift identity of A2744_YD4 via the detection of Lyman $\alpha$ emission at a redshift $z$=8.38. The association with the ALMA detection is confirmed by the presence of [OIII] 88$\mu$m emission at the same redshift. Although both emission features are only significant at the 4 $\sigma$ level, we argue their joint detection and the positional coincidence with a high redshift dropout in the HST images confirms the physical association. Analysis of the available photometric data and the modest gravitational magnification ($\mu\simeq2$) indicates A2744_YD4 has a stellar mass of $\sim$ 2$\times$10$^9$ M$_{\odot}$, a star formation rate of $\sim20$ M$_{\odot}$/yr and a dust mass of $\sim$6$\times$10$^{6}$ M$_{\odot}$. We discuss the implications of the formation of such a dust mass only $\simeq$200 Myr after the onset of cosmic reionisation.

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.

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.

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.

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

We 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.

(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.

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.

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.

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.

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.

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, 500um), 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 350um 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, <30mJy, SUSSEXtractor or DAOphot provide a better alternative to the Timeline Fitter.

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.

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.

Laporte et al. (2011) reported a very high redshift galaxy candidate: a lensed J-band dropout (A2667-J1). J1 has a photometric redshift of z=9.6-12, the probability density function for which permits no low or intermediate z solution. We here report new spectroscopic observations of this galaxy with VLT/XShooter, which show clear [OIII]5007AA, Ly-alpha, H-alpha, and H-beta emission and place the galaxy firmly at z=2.082. The oxygen lines contribute only ~25% to the H-band flux, and do not significantly affect the dropout selection of J1. After correcting the broadband fluxes for line emission, we identify two roughly equally plausible natures for A2667-J1: either it is young heavily reddened starburst, or a maximally old system with a very pronounced 4000AA break, upon which a minor secondary burst of star formation is superimposed. Fits show that to make a 3 sigma detection of this object in the B-band (V-band), imaging of depth AB=30.2 (29.5) would be required - despite the relatively bright NIR magnitude, we would need optical data of equivalent depth to the Hubble Ultra Deep Field to rule out the mid-z solution on purely photometric grounds. Assuming that this stellar population can be scaled to the NIR magnitudes of recent HST/WFC3 IR-selected galaxies, we conclude that infeasibly deep optical data AB~32 would be required for the same level of security. There is a population of galaxies at z~2 with continuum colours alone that mimic those of our z=7-12 candidates.

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.

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.

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)