Search

Stellar feedback is an essential step in the baryon cycle of galaxies, but it remains unconstrained beyond Cosmic Noon. We study the ionized gas kinematics, dynamical mass and gas-flow properties of a sample of 16 sub-$L^{\star}$ star-forming galaxies at $4\leq z\leq7.6$, using high-resolution JWST/NIRSpec observations. The emission lines are resolved, with velocity dispersions ($\sigma_{\rm gas}{\rm~(km~s^{-1})}\simeq38-96$) comparable to more massive galaxies at Cosmic Noon. From $\sigma_{\rm gas}$ and the galaxy size ($r_e=400-960~$pc), we estimate the dynamical mass to be $\log M_{\rm dyn}/M_{\odot}=9.25-10.25$. Stellar-to-dynamical mass ratios are low ($\log M_{\star}/M_{\rm dyn}\in[-0.5,-2]$) and decrease with increasing SFR surface density ($\Sigma_{\rm SFR}$). We estimate the gas surface densities assuming a star-formation law, but the gas masses do not balance the baryon-to-dynamical mass ratios, which would require a decrease in the star-formation efficiency. We find evidence of ionized outflows in five out of the sixteen galaxies, based on the need of broad components to reproduce the emission-line wings. We only observe outflows from galaxies undergoing recent bursts of star formation ${\rm SFR_{10}/SFR_{100}\geq1}$, with elevated $\Sigma_{\rm SFR}$ and low $M_{\star}/M_{\rm dyn}$. This links high gas surface densities to increased outflow incidence and lower $M_{\star}/M_{\rm dyn}$. With moderate outflow velocities ($v_{\rm flow}{\rm~(km~s^{-1})}=150-250$) and mass outflow rates ($\dot{M}_{\rm flow}/{\rm M_{\odot} yr^{-1}}=0.2-5$), these high-redshift galaxies appear more efficient at removing baryons than low-redshift galaxies with similar $M_{\star}$, showing mass loading-factors of $\dot{M}_{\rm flow}/{\rm SFR}=0.04-0.4$. For their given dynamical mass, the outflow velocities exceed the escape velocities, meaning that they may eventually enrich the Circumgalactic Medium., Comment: 15 pages, 11 figures, 3 tables; submitted to MNRAS

We present an analysis of the UV continuum slope, beta, using a sample of 733 galaxies selected from a mixture of JWST ERS/GTO/GO observational programs and with z > 4. We consider spectroscopic data obtained with the low resolution PRISM/CLEAR NIRSpec configuration. Studying the correlation of beta with M_UV we find a decreasing trend of beta = (-0.056 +- 0.017) M_UV - (3.01 +- 0.34), consistent with brighter galaxies having redder beta as found in previous works. However, analysing the trend in separate redshift bins, we find that at high redshift the relation becomes much flatter, consistent with a flat slope. Furthermore, we find that beta decreases with redshift with an evolution as beta = (-0.075 +- 0.010) z - (1.496 +- 0.056), consistent with most previous results that show a steepening of the spectra going at higher z. We then select a sample of galaxies with extremely blue slopes (beta < -2.6): such slopes are steeper than what is predicted by stellar evolution models, even for dust free, young, metal poor populations, when the contribution of nebular emission is included. We select 51 extremely blue galaxies (XBGs) and we investigate the possible physical origin of their steep slopes, comparing them to a sub-sample of redder galaxies (matched in redshift and M_UV). We find that XBGs have younger stellar populations, stronger ionization fields, lower dust attenuation, and lower but not pristine metallicity (~ 10% solar) compared to red galaxies. However, these properties alone cannot explain the extreme beta values. By using indirect inference of Lyman continuum escape, using the most recent models, we estimate escape fractions f_esc > 10% in at least 25% of XBGs, while all the red sources have smaller f_esc. A reduced nebular continuum contribution as due to either a high escape fraction or to a bursty star-formation history is likely the origin of the extremely blue slopes., Comment: Submitted to A&A (12 pages, 6 figures, and 2 tables)

Investigating the ionizing emission of star-forming galaxies is critical to understanding their contribution to reionization and their impact on the surrounding environment. The number of ionizing photons available to reionize the intergalactic medium (IGM) depends not only on the abundance of galaxies but also on their efficiency in producing ionizing photons ($\xi_{ion}$). We aim to estimate the $\xi_{ion}$ using Balmer lines in a sample of 731 galaxies at $4\leq z \leq 10$ selected from different JWST surveys. We used the available HST and JWST photometry to perform a SED fitting in the sample to determine their physical properties. We used the BAGPIPES code and assumed a delayed exponential model for the star formation history. We used the NIRSpec spectra from prism or grating configurations to estimate Balmer luminosities and then constrained $\xi_{ion}$ values after dust correction. We find a mean value of 10$^{25.23}$Hz erg$^{-1}$ for $\xi_{ion}$ in the sample with an observed scatter of 0.42dex. We find an increase in the median values of $\xi_{ion}$ which confirms the redshift evolution of $\xi_{ion}$ found in other works. Regarding the relation with physical properties, we find a decrease of $\xi_{ion}$ with increasing stellar mass, indicating that low-mass galaxies are efficient producers of ionizing photons. We also find an increase of $\xi_{ion}$ with increasing specific star formation rate (sSFR) and increasing UV absolute magnitude, which indicates that faint galaxies and with high sSFR are also efficient producers. We also investigated the relation of $\xi_{ion}$ with the EW([OIII]$\lambda$5007) and find that galaxies with the higher EW([OIII]) are the more efficient producers of ionizing photons. Similarly, we find that galaxies with higher O32 ratios and lower gas-phase metallicities (based on the R23 calibration) show higher $\xi_{ion}$ values., Comment: Submitted to A&A. Comments are welcome

We present the first spectroscopic confirmation of a dust-obscured dwarf galaxy, CEERS-14821. The analysis is performed combining JWST NIRCam broad-band photometry and NIRSpec/PRISM spectroscopic data. From the detection of multiple rest-frame optical lines, we derive that CEERS-14821 is located at $z=4.883\pm0.003$. Moreover, from a secure detection of the $H_{\alpha}$ and $H_{\beta}$ we derived that the galaxy has a dust extinction ranging from Av=2.2 to Av=3.3, depending on the assumed reddening law. This value is extremely large given that we estimated a low stellar mass around log(M/Mo)=8.0-8.2. Moreover, using different metallicity tracers, we verify that the galaxy is also metal-rich, with 12+log(O/H)>8.3. This is well above the expectation from both the mass-metallicity relation and the fundamental mass-metalliticy relation. CEERS-14821 is going through a burst of star formation, there are no indications of a strong contribution from an active galactic nuclei (f(AGN)<0.5 with respect to the total dust luminosity). Based on the rest-frame optical images, this source has a size compatible with galaxies of similar stellar mass and redshift. Finally, with the current data, it seems that there are galaxies closely interacting with CEERS-14821., Comment: 15 pages, 12 figures. Accepted by A&A on the 13th of December 2024

We present indirect constraints on the absolute escape fraction of ionizing photons ($f_{\rm esc}^{\rm LyC}$) of the system GN 42912 which comprises two luminous galaxies ($M_{\rm UV}$ magnitudes of -20.89 and -20.37) at $z\sim7.5$, GN 42912-NE and GN 42912-SW, to determine their contribution to the ionizing photon budget of the Epoch of Reionization (EoR). The high-resolution James Webb Space Telescope NIRSpec and NIRCam observations reveal the two galaxies are separated by only ~0.1$"$ (0.5 kpc) on the sky and have a 358 km s$^{-1}$ velocity separation. GN 42912-NE and GN 42912-SW are relatively massive for this redshift (log($M_\ast/M_\odot$) $\sim$ 8.4 and 8.9, respectively), with gas-phase metallicities of 18 per cent and 23 per cent solar, O$_{32}$ ratios of 5.3 and $>5.8$, and $\beta$ slopes of $-1.92$ and $-1.51$, respectively. We use the Mg II$\lambda\lambda$2796,2803 doublet to constrain $f_{\rm esc}^{\rm LyC}$. Mg II has an ionization potential close to that of neutral hydrogen and, in the optically thin regime, can be used as an indirect tracer of the LyC leakage. We establish realistic conservative upper limits on $f_{\rm esc}^{\rm LyC}$ of 8.5 per cent for GN 42912-NE and 14 per cent for GN 42912-SW. These estimates align with $f_{\rm esc}^{\rm LyC}$ trends observed with $\beta$, O$_{32}$, and the H$\beta$ equivalent width at $z<4$. The small inferred ionized region sizes ($<0.3$ pMpc) around both galaxies indicate they have not ionized a significant fraction of the surrounding neutral gas. While these $z>7$ $f_{\rm esc}^{\rm LyC}$ constraints do not decisively determine a specific reionization model, they support a minor contribution from these two relatively luminous galaxies to the EoR., Comment: 18 pages, submitted to MNRAS

We present a set of photometric catalogs primarily aimed at providing the community with a comprehensive database for the study of galaxy populations in the high redshift Universe. The set gathers data from eight JWST NIRCam observational programs, targeting the Abell 2744 (GLASS-JWST, UNCOVER, DDT2756 and GO3990), EGS (CEERS), COSMOS and UDS (PRIMER), and GOODS North and South (JADES and NGDEEP) deep fields, for a total area of $\sim$0.2 sq. degrees. Photometric estimates are obtained by means of well-established techniques, including tailored improvements designed to enhance the performance on the specific dataset. We also include new measurements from HST archival data, thus collecting 16 bands spanning from 0.44 to 4.44 $\mu$m. A grand total of $\sim$530 thousand sources is detected on stacks of NIRCam 3.56 and 4.44 $\mu$m mosaics. We assess the photometric accuracy by comparing fluxes and colors against archival catalogs. We also provide photometric redshift estimates, statistically validated against a large set of robust spectroscopic data. The catalogs are publicly available on the Astrodeep website., Comment: 28 pages, 16 figures; accepted for publication on Astronomy & Astrophysics

The Cosmic Dawn Survey (DAWN survey) provides multiwavelength (UV/optical to mid-IR) data across the combined 59 deg$^{2}$ of the Euclid Deep and Auxiliary fields (EDFs and EAFs). Here, the first public data release (DR1) from the DAWN survey is presented. DR1 catalogues are made available for a subset of the full DAWN survey that consists of two Euclid Deep fields: Euclid Deep Field North (EDF-N) and Euclid Deep Field Fornax (EDF-F). The DAWN survey DR1 catalogues do not include $Euclid$ data as they are not yet public for these fields. Nonetheless, each field has been covered by the ongoing Hawaii Twenty Square Degree Survey (H20), which includes imaging from CFHT MegaCam in the new $u$ filter and from Subaru Hyper Suprime-Cam (HSC) in the $griz$ filters. Each field is further covered by $Spitzer$/IRAC 3.6-4.5$\mu$m imaging spanning 10 deg$^{2}$ and reaching $\sim$25 mag AB (5$\sigma$). All present H20 imaging and all publicly available imaging from the aforementioned facilities are combined with the deep $Spitzer$/IRAC data to create source catalogues spanning a total area of 16.87 deg$^{2}$ in EDF-N and 2.85 deg$^{2}$ in EDF-F for this first release. Photometry is measured using The Farmer, a well-validated model-based photometry code. Photometric redshifts and stellar masses are computed using two independent codes for modeling spectral energy distributions: EAZY and LePhare. Photometric redshifts show good agreement with spectroscopic redshifts ($\sigma_{\rm NMAD} \sim 0.5, \eta < 8\%$ at $i < 25$). Number counts, photometric redshifts, and stellar masses are further validated in comparison to the COSMOS2020 catalogue. The DAWN survey DR1 catalogues are designed to be of immediate use in these two EDFs and will be continuously updated. Future data releases will provide catalogues of all EDFs and EAFs and include $Euclid$ data.

Euclid will provide deep NIR imaging to $\sim$26.5 AB magnitude over $\sim$59 deg$^2$ in its deep and auxiliary fields. The Cosmic DAWN survey complements the deep Euclid data with matched depth multiwavelength imaging and spectroscopy in the UV--IR to provide consistently processed Euclid selected photometric catalogs, accurate photometric redshifts, and measurements of galaxy properties to a redshift of $z\sim 10$. In this paper, we present an overview of the survey, including the footprints of the survey fields, the existing and planned observations, and the primary science goals for the combined data set., Comment: 16 pages, 10 figures, submitted to A&A; Updated references; Updated author list

The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance., Comment: Accepted for publication in the A&A special issue`Euclid on Sky'

Extreme emission line galaxies (EELGs) are typically characterized by high equivalent widths (EWs) which are driven by elevated specific star formation rates (sSFR) in low-mass galaxies with subsolar metallicities and little dust. Such extreme systems are rare in the local universe, but the number density of EELGs increases with redshift. Such starburst galaxies are currently presumed to be the main drivers of hydrogen reionization over 5.50.4-1 which indicates that star formation may be the dominant source of ionization. Based on the photometric fluxes, we find an increase of EW([OIII]+H$\beta$) with sSFR and $\Sigma_{SFR}$, and a decrease with age and stellar mass. The sample of EELGs can reach $\Sigma_{SFR}>$10Msun yr$^{-1}$kpc$^{-2}$ which indicate they are strong candidates of LyC leakers. Another indirect indicator is the high values of O32>5 that can be reached for some galaxies in the sample., Comment: Accepted for publication in A&A. We updated the manuscript following referee's suggestions

Here we present high spectral resolution $\textit{JWST}$ NIRSpec observations of GN42437, a low-mass (log(M$_\ast/M_\odot)=7.9$), compact ($r_e < 500$pc), extreme starburst galaxy at $z=5.59$ with 13 emission line detections. GN42437 has a low-metallicity (5-10% Z$_\odot$) and its rest-frame H$\alpha$ equivalent width suggests nearly all of the observed stellar mass formed within the last 3 Myr. GN42437 has an extraordinary 7$\sigma$ significant [Ne V] 3427 $\mathring{\rm A}$ detection. The [Ne V] line has a rest-frame equivalent width of $11\pm2\mathring{\rm A}$, [Ne V]/H$\alpha =0.04\pm0.007$, [Ne V]/[Ne III] 3870$\mathring{\rm A} = 0.26\pm0.04$, and [Ne V]/He II 4687 $\mathring{\rm A} = 1.2\pm0.5$. Ionization from massive stars, shocks, or high-mass X-ray binaries cannot simultaneously produce these [Ne V] and low-ionization line ratios. Reproducing the complete nebular structure requires both massive stars and accretion onto a black hole. We do not detect broad lines nor do the traditional diagnostics indicate that GN42437 has an accreting black hole. Thus, the very-high-ionization emission lines powerfully diagnose faint narrow-line black holes at high-redshift. We approximate the black hole mass in a variety of ways as log(M$_{\rm BH}/M_\odot) \sim 5-7$. This black hole mass is consistent with local relations between the black hole mass and the observed velocity dispersion, but significantly more massive than the stellar mass would predict. Very-high-ionization emission lines may reveal samples to probe the formation and growth of the first black holes in the universe., Comment: 18 pages, submitted to MNRAS. Comments encouraged

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

The galaxy integrated star-formation rate (SFR) surface density ($\Sigma_{\rm SFR}$) has been proposed as a valuable diagnostic of the mass accumulation in galaxies as being more tightly related to the physics of star-formation (SF) and stellar feedback than other SF indicators. In this paper, we assemble a statistical sample of 230 galaxies observed with JWST in the GLASS and CEERS spectroscopic surveys to estimate Balmer line based dust attenuations and SFRs, and UV rest-frame effective radii. We study the evolution of galaxy SFR and $\Sigma_{\rm SFR}$ in the first 1.5 Billion years of our Universe, finding that $\Sigma_{\rm SFR}$ is mildly increasing with redshift with a linear slope of $0.16 \pm 0.06$. We also explore the dependence of SFR and $\Sigma_{\rm SFR}$ on stellar mass, showing that a SF 'Main-Sequence' and a $\Sigma_{\rm SFR}$ `Main-Sequence' are in place out to z=10, with a similar slope compared to the same relations at lower redshifts. We find that the specific SFR (sSFR) and $\Sigma_{\rm SFR}$ are correlated with the [OIII]5007/[OII]3727 ratio and with indirect estimates of the escape fraction of Lyman continuum photons, hence they likely play an important role in the evolution of ionization conditions and in the escape of ionizing radiation. We also search for spectral outflow signatures in a subset of galaxies observed at high resolution, finding an outflow incidence of $2/11$ ($=20\%^{32\%}_{9\%}$) at $z<6$, but no evidence at $z>6$ ($<26\%$). Finally, we find a positive correlation between A$_V$ and $\Sigma_{\rm SFR}$, and a flat trend as a function of sSFR, indicating that there is no evidence of a drop of A$_V$ in extremely star-forming galaxies between z=4 and 10. This might be at odds with a dust-clearing outflow scenario, which might instead take place at redshifts $z\geq 10$, as suggested by some theoretical models., Comment: Accepted for publication in A&A ; 19 pages, 9 figures, 1 table. Compared to the previous version, we have corrected a typo in Fig. 5 (the x and y axis labels were inverted)

The resonant scattering interaction between Ly$\alpha$ photons and neutral hydrogen implies that a partially neutral IGM can significantly impact the detectability of Ly$\alpha$ emission in galaxies. The redshift evolution of the Ly$\alpha$ equivalent width distribution of galaxies thus offers a key probe of the degree of ionization during the Epoch of Reionization (EoR). Previous in-depth investigations at $z$ $\geq$ 7 were limited by ground-based instrument capabilities. We present an extensive study of Ly$\alpha$ emission from galaxies at 4 < $z$ < 8.5, observed from the CEERS and JADES surveys in the JWST NIRSpec/PRISM configuration. The sample consists of 235 galaxies, among which we identify 65 as Ly$\alpha$ emitters. We first measure Ly$\alpha$ escape fractions from Balmer lines, and explore the correlations with the inferred galaxies' physical properties, which are similar to those found at lower redshift. We also investigate the possible connection between the escape of Ly$\alpha$ photons and the inferred escape fractions of LyC photons obtained from indirect indicators. We then analyze the redshift evolution of the Ly$\alpha$ emitter fraction, finding lower average values at $z$ = 5 and 6 compared to ground-based observations. At $z$ = 7 we find a very large difference in Ly$\alpha$ visibility between the EGS and GOODS-South fields, possibly due to the presence of early reionized regions in the EGS. Such large variance is also expected in the Cosmic Dawn II radiation-hydrodynamical simulation. Our findings suggest a scenario in which the ending phase of the EoR is characterized by $\sim$ 1 pMpc ionized bubbles around a high fraction of moderately bright galaxies. Finally, we characterize such two ionized regions found in the EGS at $z$ = 7.18 and $z$ = 7.49 by estimating the radius of the ionized bubble that each of the spectroscopically-confirmed members could have created., Comment: Submitted to A&A

The JWST North Ecliptic Pole (NEP) Time Domain Field (TDF) is a $>$14 arcmin diameter field optimized for multi-wavelength time-domain science with JWST. It has been observed across the electromagnetic spectrum both from the ground and from space, including with the Hubble Space Telescope (HST). As part of HST observations over 3 cycles (the "TREASUREHUNT" program), deep images were obtained with ACS/WFC in F435W and F606W that cover almost the entire JWST NEP TDF. Many of the individual pointings of these programs partially overlap, allowing an initial assessment of the potential of this field for time-domain science with HST and JWST. The cumulative area of overlapping pointings is ~88 arcmin$^2$, with time intervals between individual epochs that range between 1 day and 4$+$ years. To a depth of $m_{AB}$ $\simeq$ 29.5 mag (F606W), we present the discovery of 12 transients and 190 variable candidates. For the variable candidates, we demonstrate that Gaussian statistics are applicable, and estimate that ~80 are false positives. The majority of the transients will be supernovae, although at least two are likely quasars. Most variable candidates are AGN, where we find 0.42% of the general $z$ $<$ 6 field galaxy population to vary at the $~3\sigma$ level. Based on a 5-year timeframe, this translates into a random supernova areal density of up to ~0.07 transients per arcmin$^2$ (~245 deg$^{-2}$) per epoch, and a variable AGN areal density of ~1.25 variables per arcmin$^2$ (~4500 deg$^{-2}$) to these depths., Comment: 32 pages, 11 figures, 5 tables, 1 Appendix

The Epoch of Reionization (EoR) began when galaxies grew in abundance and luminosity, so their escaping Lyman continuum (LyC) radiation started ionizing the surrounding neutral intergalactic medium (IGM). Despite significant recent progress, the nature and role of cosmic reionizers are still unclear: in order to define them, it would be necessary to directly measure their LyC escape fraction ($f_{esc}$). However, this is impossible during the EoR due to the opacity of the IGM. Consequently, many efforts at low and intermediate redshift have been made to determine measurable indirect indicators in high-redshift galaxies so that their $f_{esc}$ can be predicted. This work presents the analysis of the indirect indicators of 62 spectroscopically confirmed star-forming galaxies at $6 \leq z \leq 9$ from the Cosmic Evolution Early Release Science (CEERS) survey, combined with 12 sources with public data from other JWST-ERS campaigns. From the NIRCam and NIRSpec observations, we measured their physical and spectroscopic properties. We discovered that on average $6

We analyze the Near Infrared ($\sim0.8-1\mu$m) rest-frame morphologies of galaxies with $\log M_*/M_\odot>9$ in the redshift range $010.5$) at $z\sim5$, and bulge-dominated galaxies also exist at these early epochs, confirming a rich and evolved morphological diversity of galaxies $\sim1$ Gyr after the Big Bang. Finally, we find that the morphology-quenching relation is already in place for massive galaxies at $z>3$, with massive quiescent galaxies ($\log M_*/M_\odot>10.5$) being predominantly bulge-dominated., Comment: Submitted to A&A, comments welcome

Mid-infrared observations are powerful in identifying heavily obscured Active Galactic Nuclei (AGN) which have weak emission in other wavelengths. Data from the Mid-Infrared Instrument (MIRI) onboard JWST provides an excellent opportunity to perform such studies. We take advantage of the MIRI imaging data from the Cosmic Evolution Early Release Science Survey (CEERS) to investigate the AGN population in the distant universe. We estimate the source properties of MIRI-selected objects by utilizing spectral energy distribution (SED) modelling, and classify them into star-forming galaxies (SF), SF-AGN mixed objects, and AGN. The source numbers of these types are 418, 111, and 31, respectively, from 4 MIRI pointings covering $\sim 9$ arcmin$^2$. The sample spans a redshift range of $\approx 0$--5. We derive the median SEDs for all three source types, respectively, and publicly release them. The median MIRI SED of AGN is similar to the typical SEDs of hot dust-obscured galaxies and Seyfert 2s, for which the mid-IR SEDs are dominantly from AGN-heated hot dust. Based on our SED-fit results, we estimate the black-hole accretion density (BHAD; i.e., total BH growth rate per comoving volume) as a function of redshift. At $z<3$, the resulting BHAD agrees with the X-ray measurements in general. At $z>3$, we identify a total of 27 AGN and SF-AGN mixed objects, leading to that our high-$z$ BHAD is substantially higher than the X-ray results ($\sim 0.5$ dex at $z \approx 3$--5). This difference indicates MIRI can identify a large population of heavily obscured AGN missed by X-ray surveys at high redshifts., Comment: ApJL in press, as part of the CEERS Focus Issue. 13 pages and 7 figures (median SEDs in Fig. 5 are available upon request)

We take advantage of the NIRCam photometric observations available as part of the Cosmic Evolution Early Release Science survey (CEERS) to identify and analyse very red sources in an effort to discover very dusty star forming galaxies. We select red galaxies as objects with a S/N>3 at 4.4 $\mu$m and a S/N<2 in all JWST and HST filters at $\lambda\leq2\mu$m, which corresponds to [F200W]-[F444W]>1.2 considering CEERS depths. This selection is ideal to identify very dusty (Av>1 mag) galaxies with stellar masses between $10^6$ to $10^{10}\, \rm M_{\odot}$ at z<5, more massive dusty galaxies at z=5-18 and galaxies at z>18 due to the Lyman absorption, independently of their dust extinction. Our sample of F200W-dropouts contains no strong candidates at z>6.5, instead it consists almost completely (~81%) of z<2 low-mass galaxies, with a median stellar mass of $10^{7.3} \rm M_{\odot}$. These galaxies show an exceptional dust extinction with median value of Av=4.9 mag, completely unexpected given their low stellar mass. The remaining galaxies, which are at z<6.5, show similar large dust extinction (Av>1), but they are generally more massive $>10^{7.5}\rm M_{\odot}$., Comment: 30 pages, 1 table, 19 figures, accepted for publication in A&A

Visual inspections of the first optical rest-frame images from JWST have indicated a surprisingly high fraction of disk galaxies at high redshifts. Here, we alternatively apply self-supervised machine learning to explore the morphological diversity at $z \geq 3$. Our proposed data-driven representation scheme of galaxy morphologies, calibrated on mock images from the TNG50 simulation, is shown to be robust to noise and to correlate well with the physical properties of the simulated galaxies, including their 3D structure. We apply the method simultaneously to F200W and F356W galaxy images of a mass-complete sample ($M_*/M_\odot>10^9$) at $ 3 \leq z \leq 6$ from the first JWST/NIRCam CEERS data release. We find that the simulated and observed galaxies do not exactly populate the same manifold in the representation space from contrastive learning. We also find that half the galaxies classified as disks -- either CNN-based or visually -- populate a similar region of the representation space as TNG50 galaxies with low stellar specific angular momentum and non-oblate structure. Although our data-driven study does not allow us to firmly conclude on the true nature of these galaxies, it suggests that the disk fraction at $z \geq 3$ remains uncertain and possibly overestimated by traditional supervised classifications. Deeper imaging and spectroscopic follow-ups as well as comparisons with other simulations will help to unambiguously determine the true nature of these galaxies, and establish more robust constraints on the emergence of disks at very high redshift., Comment: 28 pages, 22 figures. Accepted for publication in ApJ

Most of the hydrogen in the intergalactic medium (IGM) was rapidly ionized at high-redshifts. While observations have established that reionization occurred, observational constraints on the emissivity of ionizing photons at high-redshift remains elusive. Here, we present a new analysis of the Low-redshift Lyman Continuum Survey (LzLCS) and archival observations, a combined sample of 89 star-forming galaxies at z~0.3 with Hubble Space Telescope observations of their ionizing continua (or Lyman Continuum, LyC). We find a strong (6$\sigma$ significant) inverse correlation between the continuum slope at 1550\r{A} (defined as F$_\lambda\propto\lambda^{\beta}$) and both the LyC escape fraction (f$_{esc}$) and f$_{esc}$ times the ionizing photon production efficiency ($\xi_{ ion}$). On average, galaxies with redder continuum slopes have smaller f$_{esc}$ than galaxies with bluer slopes due to higher dust attenuation. More than 5% (20%) of the LyC emission escapes galaxies with $\beta$<-2.1 (-2.6). We find strong correlations between $\beta$ and the gas-phase ionization ([OIII]/[OII] flux ratio; at 7.5$\sigma$ significance), galaxy stellar mass (at 5.9$\sigma$), the gas-phase metallicity (at 4.6$\sigma$), and the observed FUV absolute magnitude (M$_{UV}$ at 3.4$\sigma$). Using previous observations of $\beta$ at high-redshift, we estimate the evolution of f$_{esc}$ with both $z$ and M$_{UV}$. The LzLCS suggest that fainter and lower mass galaxies dominate the ionizing photon budget at higher redshift, possibly due to their rapidly evolving metal and dust content. Finally, we use our correlation between $\beta$ and f$_{ esc}\times\xi_{ion}$ to predict the ionizing emissivity of galaxies during the epoch of reionization. Our estimated emissivities match IGM observations, and suggest that star-forming galaxies emit sufficient LyC photons into the IGM to exceed recombinations near redshifts of 7-8., Comment: 14 pages plus appendix. Submitted to MNRAS. Comments encouraged

(Abridged) The Euclid mission is expected to discover thousands of z>6 galaxies in three Deep Fields, which together will cover a ~40 deg2 area. However, the limited number of Euclid bands and availability of ancillary data could make the identification of z>6 galaxies challenging. In this work, we assess the degree of contamination by intermediate-redshift galaxies (z=1-5.8) expected for z>6 galaxies within the Euclid Deep Survey. This study is based on ~176,000 real galaxies at z=1-8 in a ~0.7 deg2 area selected from the UltraVISTA ultra-deep survey, and ~96,000 mock galaxies with 25.3$\leq$H<27.0, which altogether cover the range of magnitudes to be probed in the Euclid Deep Survey. We simulate Euclid and ancillary photometry from the fiducial, 28-band photometry, and fit spectral energy distributions (SEDs) to various combinations of these simulated data. Our study demonstrates that identifying z>6 with Euclid data alone will be very effective, with a z>6 recovery of 91(88)% for bright (faint) galaxies. For the UltraVISTA-like bright sample, the percentage of z=1-5.8 contaminants amongst apparent z>6 galaxies as observed with Euclid alone is 18%, which is reduced to 4(13)% by including ultra-deep Rubin (Spitzer) photometry. Conversely, for the faint mock sample, the contamination fraction with Euclid alone is considerably higher at 39%, and minimized to 7% when including ultra-deep Rubin data. For UltraVISTA-like bright galaxies, we find that Euclid (I-Y)>2.8 and (Y-J)<1.4 colour criteria can separate contaminants from true z>6 galaxies, although these are applicable to only 54% of the contaminants, as many have unconstrained (I-Y) colours. In the most optimistic scenario, these cuts reduce the contamination fraction to 1% whilst preserving 81% of the fiducial z>6 sample. For the faint mock sample, colour cuts are infeasible., Comment: 27 pages, 14 figures, 9 tables; version accepted for publication by A&A

Compact star formation appears to be generally common in dusty star-forming galaxies (SFGs). However, its role in the framework set by the scaling relations in galaxy evolution remains to be understood. In this work we follow up on the galaxy sample from the GOODS-ALMA 2.0 survey, an ALMA blind survey at 1.1mm covering a continuous area of 72.42arcmin$^2$ using two array configurations. We derived physical properties, such as star formation rates, gas fractions, depletion timescales, and dust temperatures for the galaxy sample built from the survey. There exists a subset of galaxies that exhibit starburst-like short depletion timescales, but they are located within the scatter of the so-called main sequence of SFGs. These are dubbed starbursts in the main sequence and display the most compact star formation and they are characterized by the shortest depletion timescales, lowest gas fractions, and highest dust temperatures of the galaxy sample, compared to typical SFGs at the same stellar mass and redshift. They are also very massive, accounting for $\sim 60\%$ of the most massive galaxies in the sample ($\log (M_{\rm{*}}/M_{\odot}) > 11.0$). We find trends between the areas of the ongoing star formation regions and the derived physical properties for the sample, unveiling the role of compact star formation as a physical driver of these properties. Starbursts in the main sequence appear to be the extreme cases of these trends. We discuss possible scenarios of galaxy evolution to explain the results drawn from our galaxy sample. Our findings suggest that the star formation rate is sustained in SFGs by gas and star formation compression, keeping them within the main sequence even when their gas fractions are low and they are presumably on the way to quiescence., Comment: Accepted by A&A. 26 pages, 13 figures

Submillimeter/millimeter observations of dusty star-forming galaxies with the Atacama Large Millimeter/submillimeter Array (ALMA) have shown that dust continuum emission generally occurs in compact regions smaller than the stellar distribution. However, it remains to be understood how systematic these findings are. Studies often lack homogeneity in the sample selection, target discontinuous areas with inhomogeneous sensitivities, and suffer from modest $uv$ coverage coming from single array configurations. GOODS-ALMA is a 1.1mm galaxy survey over a continuous area of 72.42arcmin$^2$ at a homogeneous sensitivity. In this version 2.0, we present a new low resolution dataset and its combination with the previous high resolution dataset from the survey, improving the $uv$ coverage and sensitivity reaching an average of $\sigma = 68.4\mu$Jy beam$^{-1}$. A total of 88 galaxies are detected in a blind search (compared to 35 in the high resolution dataset alone), 50% at $S/N_{peak} \geq 5$ and 50% at $3.5 \leq S/N_{peak} \leq 5$ aided by priors. Among them, 13 out of the 88 are optically dark or faint sources ($H$- or $K$-band dropouts). The sample dust continuum sizes at 1.1mm are generally compact, with a median effective radius of $R_{e} = 0"10 \pm 0"05$ (a physical size of $R_{e} = 0.73 \pm 0.29$kpc at the redshift of each source). Dust continuum sizes evolve with redshift and stellar mass resembling the trends of the stellar sizes measured at optical wavelengths, albeit a lower normalization compared to those of late-type galaxies. We conclude that for sources with flux densities $S_{1.1mm} > 1$mJy, compact dust continuum emission at 1.1mm prevails, and sizes as extended as typical star-forming stellar disks are rare. The $S_{1.1mm} < 1$mJy sources appear slightly more extended at 1.1mm, although they are still generally compact below the sizes of typical star-forming stellar disks., Comment: Accepted by A&A. 30 pages, 18 figures

NASA's Great Observatories have opened up the electromagnetic spectrum from space, providing sustained access to wavelengths not accessible from the ground. Together, Hubble, Compton, Chandra, and Spitzer have provided the scientific community with an agile and powerful suite of telescopes with which to attack broad scientific questions, and react to a rapidly changing scientific landscape. As the existing Great Observatories age, or are decommissioned, community access to these wavelengths will diminish, with an accompanying loss of scientific capability. This report, commissioned by the NASA Cosmic Origins, Physics of the Cosmos and Exoplanet Exploration Program Analysis Groups (PAGs), analyzes the importance of multi-wavelength observations from space during the epoch of the Great Observatories, providing examples that span a broad range of astrophysical investigations., Comment: A report by the Great Observatories Science Analysis Group (SAG-10), commissioned by NASA's Cosmic Origins, Physics of the Cosmos, and Exoplanet Exploration Program Analysis Groups. 87 pages, 23 figures

VANDELS is an ESO Public Spectroscopic Survey designed to build a sample of high signal to noise, medium resolution spectra of galaxies at redshift between 1 and 6.5. Here we present the final Public Data Release of the VANDELS Survey, comprising 2087 redshift measurements. We give a detailed description of sample selection, observations and data reduction procedures. The final catalogue reaches a target selection completeness of 40% at iAB = 25. The high Signal to Noise ratio of the spectra (above 7 in 80% of the spectra) and the dispersion of 2.5{\AA} allowed us to measure redshifts with high precision, the redshift measurement success rate reaching almost 100%. Together with the redshift catalogue and the reduced spectra, we also provide optical mid-IR photometry and physical parameters derived through SED fitting. The observed galaxy sample comprises both passive and star forming galaxies covering a stellar mass range 8.3< Log(M*/Msolar)<11.7. All catalogues and spectra are accessible through the survey database (http://vandels.inaf.it) where all information can be queried interactively, and via the ESO Archive (https://www.eso.org/qi/)., Comment: 16 pages, 13 figures, accepted for publication in Astronomy & Astrophysics

We present the first results from the 2mm Mapping Obscuration to Reionization (MORA) survey, the largest ALMA contiguous blank-field survey to-date with a total area of 184 sq. arcmin and the only at 2mm to search for dusty star-forming galaxies (DSFGs). We use the 13 sources detected above 5sigma to estimate the first ALMA galaxy number counts at this wavelength. These number counts are then combined with the state-of-the-art galaxy number counts at 1.2mm and 3mm and with a backward evolution model to place constraints on the evolution of the IR luminosity function and dust-obscured star formation in the last 13 billion years. Our results suggest a steep redshift evolution on the space density of DSFGs and confirm the flattening of the IR luminosity function at faint luminosities, with a slope of $\alpha_{LF} = -0.42^{+0.02}_{-0.04}$. We conclude that the dust-obscured component, which peaks at z=2-2.5, has dominated the cosmic history of star formation for the past ~12 billion years, back to z~4. At z=5, the dust-obscured star formation is estimated to be ~35% of the total star formation rate density and decreases to 25%-20% at z=6-7, implying a minor contribution of dust-enshrouded star formation in the first billion years of the Universe. With the dust-obscured star formation history constrained up to the end of the epoch of reionization, our results provide a benchmark to test galaxy formation models, to study the galaxy mass assembly history, and to understand the dust and metal enrichment of the Universe at early times., Comment: Accepted for publication in The Astrophysical Journal

In this paper we study the properties of the six optically dark galaxies detected in the 69 arcmin^2 GOODS-ALMA 1.1mm continuum survey. While none of them are listed in the deepest H-band based CANDELS catalog in the GOODS-South field down to H=28.16AB, we were able to de-blend two of them from their bright neighbor and measure an $H$-band flux for them. We note that AGS4 and AGS15 have H=25.23, 27.11AB respectively. Their extreme proximity (0.50", 0.27") to a bright optical source and their extreme faintness prevented them from being included in the H-band catalog. We present the spectroscopic scan follow-up of five of the six sources with ALMA band 4. All are detected in the 2mm continuum with signal-to-noise ratios higher than eight. One emission line is detected in AGS4 (\nu_{obs} =151.44GHz with a S/N=8.58) and AGS17 (\nu_{obs} =154.78GHz with a S/N=10.23), which we interpret in both cases as being due to the CO(6-5) line at z^{AGS4}_{spec}=3.556 and z^{AGS4}_{spec}=3.467, respectively. These redshifts match both the probability distribution of the photometric redshifts derived from the UV to near-infrared spectral energy distributions (SEDs) and the far-infrared SEDs for typical dust temperatures of galaxies at these redshifts. We present evidence that nearly 70% (4/6 of galaxies) of the optically dark galaxies belong to the same overdensity of galaxies at z~3.5. overdensity The most massive one, AGS24 (M_{\star} = 10^{11.32^{+0.02}_{-0.19}} M_{\odot}), is the most massive galaxy without an active galactic nucleus (AGN) at $z$\,>\,3 in the GOODS-ALMA field. It falls in the very center of the peak of the galaxy surface density, which suggests that the surrounding overdensity is a proto-cluster in the process of virialization and that AGS24 is the candidate progenitor of the future brightest cluster galaxy (BCG)., Comment: 17 pages, 13 figures, 3 tables, accepted by A&A

We investigate the properties of a sample of 35 galaxies, detected with ALMA at 1.1 mm in the GOODS-ALMA field (area of 69 arcmin$^2$, resolution = 0.60", RMS $\simeq$ 0.18 mJy beam$^{-1}$). Using the UV-to-radio deep multiwavelength coverage of the GOODS-South field, we fit the spectral energy distributions of these galaxies to derive their key physical properties. The galaxies detected by ALMA are among the most massive at $z$ = 2-4 (M$_{\star,med}$ = 8.5$ \times$ 10$^{10}$ M$_\odot$) and are either starburst or located in the upper part of the galaxy star-forming main sequence. A significant portion of our galaxy population ($\sim$ 40%), located at $z\sim$ 2.5-3, exhibits abnormally low gas fractions. The sizes of these galaxies, measured with ALMA, are compatible with the trend between $H$-band size and stellar mass observed for $z\sim2$ elliptical galaxies suggesting that they are building compact bulges. We show that there is a strong link between star formation surface density (at 1.1 mm) and gas depletion time: the more compact a galaxy's star-forming region is, the shorter its lifetime will be (without gas replenishment). The identified compact sources associated with relatively short depletion timescales ($\sim$100 Myr), are the ideal candidates to be the progenitors of compact elliptical galaxies at $z$ $\sim$ 2., Comment: 18 pages, 8 figures, Accepted for publication in A&A

In this paper, we extend the source detection in the GOODS-ALMA field (69 arcmin$^2$, rms sensitivity $\sigma$ $\simeq$ 0.18 mJy.beam$^{-1}$), to deeper levels than presented in Franco et al. (2018). Using positional information at 3.6 and 4.5 $\mu$m (from Spitzer-IRAC), we explore the presence of galaxies detected at 1.1 mm with ALMA below our original blind detection limit of 4.8-$\sigma$ at which the number of spurious sources starts to dominate over that of real sources. In this Supplementary Catalog, we find a total of 16 galaxies, including 2 galaxies with no counterpart in HST images (also known as optically-dark galaxies) down to a 5$\sigma$ limiting depth of H = 28.2 AB (HST/WFC3 F160W). This brings the total sample of GOODS-ALMA 1.1 mm sources to 35 galaxies. Galaxies in the new sample cover a wider dynamic range in redshift ($z$ = 0.65 - 4.73), are on average twice as large (1.3 vs 0.65 kpc) and and have lower stellar mass (M$_{\star}^{\rm SC}$ = 7.6$\times$10$^{10}$M$_\odot$ vs M$_{\star}^{\rm MC}$ = 1.2$\times$10$^{11}$M$_\odot$). Although exhibiting larger physical sizes, these galaxies have still far-infrared sizes significantly more compact than inferred from their optical emission. We show that the astrometry of the HST image does not only suffer from a global astrometric shift, as already discussed in previous papers, but also from local shifts. These distortions were artificially introduced in the process of building the mosaic of the GOODS-South HST image. By comparing the positions of almost 400 galaxies detected by HST, Pan-STARRS and ALMA, we create a distortion map which can be used to correct for these astrometric issues., Comment: 18 pages, 13 figures. Accepted for publication in A&A

While most of the inter-galactic medium (IGM) today is permeated by ionized hydrogen, it was largely filled with neutral hydrogen for the first 700 million years after the Big Bang. The process that ionized the IGM (cosmic reionization) is expected to be spatially inhomogeneous, with fainter galaxies playing a significant role. However, we still have only a few direct constraints on the reionization process. Here we report the first spectroscopic confirmation of two galaxies and very likely a third galaxy in a group (hereafter EGS77) at redshift z = 7.7, merely 680 Myrs after the Big Bang. The physical separation among the three members is < 0.7 Mpc. We estimate the radius of ionized bubble of the brightest galaxy to be about 1.02 Mpc, and show that the individual ionized bubbles formed by all three galaxies likely overlap significantly, forming a large yet localized ionized region, which leads to the spatial inhomogeneity in the reionization process. It is striking that two of three galaxies in EGS77 are quite faint in the continuum, thanks to our selection of reionizing sources using their Lyman-alpha line emission. Indeed, one is the faintest spectroscopically confirmed galaxy yet discovered at such high redshifts. Our observations provide direct constraints in the process of cosmic reionization, and allow us to investigate the properties of sources responsible for reionizing the universe., Comment: Submitted after addressing referee's comments

Aims. Our aim is to estimate the intergalactic medium transmission towards UV-selected star-forming galaxies at redshift 4 and above and study the effect of the dust attenuation on these measurements. Methods. The ultra-violet spectrum of high redshift galaxies is a combination of their intrinsic emission and the effect of the Inter-Galactic medium (IGM) absorption along their line of sight. Using data coming from the unprecedented deep spectroscopy from the VANDELS ESO public survey carried out with the VIMOS instrument we compute both the dust extinction and the mean transmission of the IGM as well as its scatter from a set of 281 galaxies at z>3.87. Because of a degeneracy between the dust content of the galaxy and the IGM, we first estimate the stellar dust extinction parameter E(B-V) and study the result as a function of the dust prescription. Using these measurements as constraint for the spectral fit we estimate the IGM transmission Tr(Lyalpha). Both photometric and spectroscopic SED fitting are done using the SPectroscopy And photometRy fiTting tool for Astronomical aNalysis (SPARTAN) that is able to fit the spectral continuum of the galaxies as well as photometric data. Results. Using the classical Calzetti's attenuation law we find that E(B-V) goes from 0.11 at z=3.99 to 0.08 at z=5.15. These results are in very good agreement with previous measurements from the literature. We estimate the IGM transmission and find that the transmission is decreasing with increasing redshift from Tr(Lyalpha)=0.53 at z=3.99 to 0.28 at z=5.15. We also find a large standard deviation around the average transmission that is more than 0.1 at every redshift. Our results are in very good agreement with both previous measurements from AGN studies and with theoretical models., Comment: Accepted for publication in A&A

We investigate the physical properties of a sample of 52 Ly{\alpha} emitting galaxies in the VANDELS survey, with particular focus on the role of kinematics and neutral hydrogen column density in the escape and spatial distribution of Ly{\alpha} photons. For these galaxies, we derive different physical properties and characterize the Ly{\alpha} emission in terms of kinematics, EW, FWHM and spatial extension and then estimate the velocity of the neutral outflowing gas. We reproduce some of the well known correlations between Ly{\alpha} EW and stellar mass, dust extinction and UV \b{eta} slope, in the sense that the emission line appears brighter in lower mass, less dusty and bluer galaxies. We do not find any correlation with the SED-derived star formation rate, while we find that galaxies with brighter Ly{\alpha} tend to be more compact both in UV and in Ly{\alpha}. Our data reveal a new interesting correlation between the Ly{\alpha} velocity and the offset of the inter-stellar absorption lines with respect to the systemic redshift, in the sense that galaxies with larger inter-stellar medium (ISM) out-flow velocities show smaller Ly{\alpha} velocity shifts. We interpret this relation in the context of the shell-model scenario, where the velocity of the ISM and the HI column density contribute together in determining the Ly{\alpha} kinematics. In support to our interpretation, we observe that galaxies with high HI column densities have much more extended Ly{\alpha} spatial profiles, a sign of increased scattering. However, we do not find any evidence that the HI column density is related to any other physical properties of the galaxies, although this might be due in part to the limited range of parameters that our sample spans., Comment: submitted to A&A

Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.)

Recent near-IR power-spectra and panchromatic Extragalactic Background Light measurements provide upper limits on the near-IR surface brightness (SB>31 mag/arcsec^2) that may come from Pop III stars and accretion disks around resulting stellar-mass black holes (BHs) in the epoch of First Light (z=7-17). Physical parameters for zero metallicity Pop III stars at z>7 can be estimated from MESA stellar evolution models through helium-depletion, and for BH accretion disks from quasar microlensing results and multicolor accretion models. Second-generation stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions in their AGB stage. The near-IR SB constraints can be used to calculate the number of caustic transits behind lensing clusters that JWST and the 25~39 m ground-based telescopes may detect for both Pop III stars and stellar mass BH accretion disks. Because Pop III stars and stellar mass BH accretion disks have sizes of a few x 10^-11 arcsec at z>7, typical caustic magnifications can be mu=10^4~10^5, with rise times of hours and decline times of < 1 year for cluster transverse velocities of v<1000 km/s. Microlensing by intracluster medium objects can modify transit magnifications, and lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Pop III stars. To observe Pop III caustic transits directly may require monitoring 3~30 lensing clusters to AB< 29 mag over a decade or more. Such a program must be started with JWST in Cycle 1, and -- depending on the role of microlensing in the Intra Cluster Light -- should be continued for decades with the GMT and TMT, where JWST and the ground-based telescopes each will play a unique and strongly complementary role., Comment: 6 pages + references, 3 figures, science white paper submitted to the Astro2020 decadal survey. arXiv admin note: substantial text overlap with arXiv:1801.03584

We present the results of CANDELSz7, an ESO large program aimed at confirming spectroscopically a homogeneous sample of z~6 and z~7 star forming galaxies. The candidates were selected in the GOODS-South, UDS and COSMOS fields using the official CANDELS catalogs based on H160-band detections. Standard color criteria, which were tailored depending on the ancillary multi-wavelength data available for each field, were applied to select more than 160 candidate galaxies at z~6 and z~7. Deep medium resolution FORS2 spectroscopic observations were then conducted with integration times ranging from 12 to 20 hours, to reach a Lyalpha flux limit of approximately 1-3x 10-18 erg/s/cm^2 at 3sigma. For about 40% of the galaxies we could determine a spectroscopic redshift, mainly through the detection of a single emission line that we interpret as Lyalpha emission, or for some of the brightest objects (H160< 25.5) from the presence of faint continuum and sharp drop that we interpret as a Lyman break. In this paper we present the redshifts and main properties of 65 newly confirmed high redshift galaxies. Adding previous proprietary and archival data we assemble a sample of ~260 galaxies that we use to explore the evolution of the Lyalpha fraction in Lyman break galaxies and the change in the shape of the emission line between z~6 and z~7. We also discuss the accuracy of the CANDELS photometric redshifts in this redshift range., Comment: Accepted for publication in Astronomy & Astrophysics

This paper describes the observations and the first data release (DR1) of the ESO public spectroscopic survey "VANDELS, a deep VIMOS survey of the CANDELS CDFS and UDS fields". VANDELS' main targets are star-forming galaxies at 2.4

VANDELS is a uniquely-deep spectroscopic survey of high-redshift galaxies with the VIMOS spectrograph on ESO's Very Large Telescope (VLT). The survey has obtained ultra-deep optical (0.48 < lambda < 1.0 micron) spectroscopy of ~2100 galaxies within the redshift interval 1.0 < z < 7.0, over a total area of ~0.2 sq. degrees centred on the CANDELS UDS and CDFS fields. Based on accurate photometric redshift pre-selection, 85% of the galaxies targeted by VANDELS were selected to be at z>=3. Exploiting the red sensitivity of the refurbished VIMOS spectrograph, the fundamental aim of the survey is to provide the high signal-to-noise ratio spectra necessary to measure key physical properties such as stellar population ages, masses, metallicities and outflow velocities from detailed absorption-line studies. Using integration times calculated to produce an approximately constant signal-to-noise ratio (20 < t_int < 80 hours), the VANDELS survey targeted: a) bright star-forming galaxies at 2.4 < z < 5.5, b) massive quiescent galaxies at 1.0 < z < 2.5, c) fainter star-forming galaxies at 3.0 < z < 7.0 and d) X-ray/Spitzer-selected active galactic nuclei and Herschel-detected galaxies. By targeting two extragalactic survey fields with superb multi-wavelength imaging data, VANDELS will produce a unique legacy data set for exploring the physics underpinning high-redshift galaxy evolution. In this paper we provide an overview of the VANDELS survey designed to support the science exploitation of the first ESO public data release, focusing on the scientific motivation, survey design and target selection., Comment: 19 pages, 9 figures, accepted for publication in MNRAS

We present a 69 arcmin$^2$ ALMA survey at 1.1mm, GOODS-ALMA, matching the deepest HST-WFC3 H-band part of the GOODS-South field. We taper the 0"24 original image with a homogeneous and circular synthesized beam of 0"60 to reduce the number of independent beams - thus reducing the number of purely statistical spurious detections - and optimize the sensitivity to point sources. We extract a catalogue of galaxies purely selected by ALMA and identify sources with and without HST counterparts down to a 5$\sigma$ limiting depth of H=28.2 AB (HST/WFC3 F160W). ALMA detects 20 sources brighter than 0.7 mJy in the 0"60 tapered mosaic (rms sensitivity =0.18 mJy/beam) with a purity greater than 80%. Among these detections, we identify three sources with no HST nor Spitzer-IRAC counterpart, consistent with the expected number of spurious galaxies from the analysis of the inverted image; their definitive status will require additional investigation. An additional three sources with HST counterparts are detected either at high significance in the higher resolution map, or with different detection-algorithm parameters ensuring a purity greater than 80%. Hence we identify in total 20 robust detections. Our wide contiguous survey allows us to push further in redshift the blind detection of massive galaxies with ALMA with a median redshift of $z$=2.92 and a median stellar mass of M$_{\star}$ = 1.1 $\times 10^{11}$M$_\odot$. Our sample includes 20% HST-dark galaxies (4 out of 20), all detected in the mid-infrared with IRAC. The near-infrared based photometric redshifts of two of them $z\sim$4.3 and 4.8) suggest that these sources have redshifts $z$>4. At least 40% of the ALMA sources host an X-ray AGN, compared to 14% for other galaxies of similar mass and redshift. The wide area of our ALMA survey provides lower values at the bright end of number counts than single-dish telescopes, Comment: 26 pages, 23 figures, accepted to A&A

JWST was conceived and built to answer one of the most fundamental questions that humans can address empirically: "How did the Universe make its first stars?". Our First Lights At REionization (FLARE) project transforms the quest for the epoch of reionization from the static to the time domain. It targets the complementary question: "What happened to those first stars?". It will be answered by observations of the most luminous events: supernovae and accretion on to black holes formed by direct collapse from the primordial gas clouds. These transients provide direct constraints on star-formation rates (SFRs) and the truly initial Initial Mass Function (IMF), and they may identify possible stellar seeds of supermassive black holes (SMBHs). Furthermore, our knowledge of the physics of these events at ultra-low metallicity will be much expanded. JWST's unique capabilities will detect these most luminous and earliest cosmic messengers easily in fairly shallow observations. However, these events are very rare at the dawn of cosmic structure formation and so require large area coverage. Time domain astronomy can be advanced to an unprecedented depth by means of a shallow field of JWST reaching 27 mag AB in 2 and 4.4 microns over a field as large as 0.1 square degree visited multiple times each year. Such a survey may set strong constraints or detect massive Pop III SNe at redshifts beyond 10, pinpointing the redshift of the first stars, or at least their death. Based on our current knowledge of superluminous supernovae (SLSNe), such a survey will find one or more SLSNe at redshifts above 6 in five years and possibly several direct collapse black holes. Although JWST is not designed as a wide field survey telescope, we show that such a wide field survey is possible with JWST and is critical in addressing several of its key scientific goals.

The fraction of Lyman-$\alpha$ emitters among the galaxy population has been found to increase from $z\sim0$ to $z\sim6$ and drop dramatically at $z>6$. This drop has been interpreted as an effect of an increasingly neutral intergalactic medium with increasing redshift, while a LyC escape fraction evolving with redshift. We report the result of a large VLT/FORS2 program aiming to confirm spectroscopically a large galaxy sample at $z\geq6$ that has been selected in several independent fields through the Lyman Break technique. Combining those data with archival data, we create a large and homogeneous sample of $z\sim6$ galaxies ($N=127$), complete in terms of Ly$\alpha$ detection at $>95\%$ for EW(Ly$\alpha)\geq25\AA$. We use this sample to derive a new measurement of the LAE fraction at $z\sim6$ and derive the physical properties of these galaxies through spectral energy distribution fitting. We find a median LAE fraction at $z\sim6$ lower than in previous studies. The main difference between LAEs and non-LAEs is that the latter are significantly dustier. Using predictions of our SED fitting code accounting for nebular emission, we find an effective Ly$\alpha$ escape fraction $f^{eff}_{esc}(Ly\alpha)=0.23^{+0.36}_{-0.17}$ remarkably consistent with the value derived by comparing UV luminosity function with Ly$\alpha$ luminosity function. We conclude that the drop in the LAE fraction from $z\sim6$ to $z>6$ is less dramatic than previously found and the effect of an increasing IGM neutral fraction is possibly observed at $5

We analyze a sample of $z$-dropout galaxies in the CANDELS GOODS South and UDS fields that have been targeted by a dedicated spectroscopic campaign aimed at detecting their Ly$\alpha$ line. Deep IRAC observations at 3.6 and 4.5 $\mu$m are used to determine the strength of optical emission lines affecting these bands at z$\sim$6.5-6.9 in order to i) investigate possible physical differences between Ly$\alpha$ emitting and non-emitting sources; ii) constrain the escape fraction of ionizing photons; iii) provide an estimate of the specific star-formation rate at high redshifts. We find evidence of strong [OIII]+H$\beta$ emission in the average (stacked) SEDs of galaxies both with and without Ly$\alpha$ emission. The blue IRAC [3.6]-[4.5] color of the stack with detected Ly$\alpha$ line can be converted into a rest-frame equivalent width EW([OIII]+H$\beta$)=1500$^{+530}_{-440}\AA$ assuming a flat intrinsic stellar continuum. This strong optical line emission enables a first estimate of f$_{esc}\lesssim$20% on the escape fraction of ionizing photons from Ly$\alpha$ detected objects. The objects with no Ly$\alpha$ line show less extreme EW([OIII]+H$\beta$)=520$^{+170}_{-150}\AA$ suggesting different physical conditions of the HII regions with respect to Ly$\alpha$-emitting ones, or a larger f$_{esc}$. The latter case is consistent with a combined evolution of f$_{esc}$ and the neutral hydrogen fraction as an explanation of the lack of bright Ly$\alpha$ emission at z$>$6. A lower limit on the specific star formation rate, SSFR$>$9.1$Gyr^{-1}$ for $M_{star}=2 \times 10^9 M_{\odot}$ galaxies at these redshifts can be derived from the spectroscopically confirmed sample., Comment: 10 pages, 7 figures; The Astrophysical Journal in press; matched to the published version

We present a multi-wavelength photometric catalog in the COSMOS field as part of the observations by the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). The catalog is based on Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) and Advanced Camera for Surveys (ACS) observations of the COSMOS field (centered at RA: $10^h00^m28^s$, Dec:$+02^{\circ}12^{\prime}21^{\prime\prime}$). The final catalog has 38671 sources with photometric data in forty two bands from UV to the infrared ($\rm \sim 0.3-8\,\mu m$). This includes broad-band photometry from the HST, CFHT, Subaru, VISTA and Spitzer Space Telescope in the visible, near infrared and infrared bands along with intermediate and narrow-band photometry from Subaru and medium band data from Mayall NEWFIRM. Source detection was conducted in the WFC3 F160W band (at $\rm 1.6\,\mu m$) and photometry is generated using the Template FITting algorithm. We further present a catalog of the physical properties of sources as identified in the HST F160W band and measured from the multi-band photometry by fitting the observed spectral energy distributions of sources against templates., Comment: 28 Pages, 22 Figures and 10 Tables. Accepted for Publication in the ApJS

Recent analysis of strongly-lensed sources in the Hubble Frontier Fields indicates that the rest-frame UV luminosity function of galaxies at $z=$6--8 rises as a power law down to $M_\mathrm{UV}=-15$, and possibly as faint as -12.5. We use predictions from a cosmological radiation hydrodynamic simulation to map these luminosities onto physical space, constraining the minimum dark matter halo mass and stellar mass that the Frontier Fields probe. While previously-published theoretical studies have suggested or assumed that early star formation was suppressed in halos less massive than $10^9$--$10^{11} M_\odot$, we find that recent observations demand vigorous star formation in halos at least as massive as (3.1, 5.6, 10.5)$\times10^9 M_\odot$ at $z=(6,7,8)$. Likewise, we find that Frontier Fields observations probe down to stellar masses of (8.1, 18, 32)$\times10^6 M_\odot$; that is, they are observing the likely progenitors of analogues to Local Group dwarfs such as Pegasus and M32. Our simulations yield somewhat different constraints than two complementary models that have been invoked in similar analyses, emphasizing the need for further observational constraints on the galaxy-halo connection., Comment: 8 pages, 5 figures, re-submitted to MNRAS after incorporating referee's comments

We present the results of the first, deep ALMA imaging covering the full 4.5 sq arcmin of the Hubble Ultra Deep Field (HUDF) as previously imaged with WFC3/IR on HST. Using a mosaic of 45 pointings, we have obtained a homogeneous 1.3mm image of the HUDF, achieving an rms sensitivity of 35 microJy, at a resolution of 0.7 arcsec. From an initial list of ~50 >3.5sigma peaks, a rigorous analysis confirms 16 sources with flux densities S(1.3) > 120 microJy. All of these have secure galaxy counterparts with robust redshifts ( = 2.15), and 12 are also detected at 6GHz in new deep JVLA imaging. Due to the wealth of supporting data in this unique field, the physical properties of the ALMA sources are well constrained, including their stellar masses (M*) and UV+FIR star-formation rates (SFR). Our results show that stellar mass is the best predictor of SFR in the high-z Universe; indeed at z > 2 our ALMA sample contains 7 of the 9 galaxies in the HUDF with M* > 2 x 10^10 Msun and we detect only one galaxy at z > 3.5, reflecting the rapid drop-off of high-mass galaxies with increasing redshift. The detections, coupled with stacking, allow us to probe the redshift/mass distribution of the 1.3-mm background down to S(1.3) ~ 10 micro-Jy. We find strong evidence for a steep `main sequence' for star-forming galaxies at z ~ 2, with SFR \propto M* and a mean specific SFR = 2.2 /Gyr. Moreover, we find that ~85% of total star formation at z ~ 2 is enshrouded in dust, with ~65% of all star formation at this epoch occurring in high-mass galaxies (M* > 2 x 10^10 Msun), for which the average obscured:unobscured SF ratio is ~200. Finally, we combine our new ALMA results with the existing HST data to revisit the cosmic evolution of star-formation rate density; we find that this peaks at z ~ 2.5, and that the star-forming Universe transits from primarily unobscured to primarily obscured thereafter at z ~ 4., Comment: 26 pages, 15 figures, updated to match version accepted for publication in MNRAS

The Frontier Fields are a director's discretionary time campaign with HST and the Spitzer Space Telescope to see deeper into the universe than ever before. The Frontier Fields combine the power of HST and Spitzer with the natural gravitational telescopes of massive high-magnification clusters of galaxies to produce the deepest observations of clusters and their lensed galaxies ever obtained. Six clusters - Abell 2744, MACSJ0416.1-2403, MACSJ0717.5+3745, MACSJ1149.5+2223, Abell S1063, and Abell 370 - were selected based on their lensing strength, sky darkness, Galactic extinction, parallel field suitability, accessibility to ground-based facilities, HST, Spitzer and JWST observability, and pre-existing ancillary data. These clusters have been targeted by the HST ACS/WFC and WFC3/IR with coordinated parallels of adjacent blank fields for over 840 HST orbits. The Spitzer Space Telescope has dedicated > 1000 hours of director's discretionary time to obtain IRAC 3.6 and 4.5 micron imaging to ~26.5, 26.0 ABmag 5-sigma point-source depths in the six cluster and six parallel Frontier Fields. The Frontier Field parallel fields are the second-deepest observations thus far by HST with ~29th ABmag 5-sigma point source depths in seven optical - near-infrared bandpasses. Galaxies behind the Frontier Field cluster lenses experience typical magnification factors of a few, with small regions near the critical curves magnified by factors 10-100. Therefore, the Frontier Field cluster HST images achieve intrinsic depths of ~30-33 magnitudes over very small volumes. Early studies of the Frontier Fields have probed galaxies fainter than any seen before during the epoch of reionization 6 < z < 10, mapped out the cluster dark matter to unprecedented resolution, and followed lensed transient events., Comment: submitted to ApJ; 18 pages; see http://www.stsci.edu/hst/campaigns/frontier-fields/ and http://ssc.spitzer.caltech.edu/warmmission/scheduling/approvedprograms/ddt/frontier/ for data and more information

Galaxies at high redshifts provide a valuable tool to study cosmic dawn, and therefore it is crucial to reliably identify these galaxies. Here, we present an unambiguous and first simultaneous detection of both the Lyman-alpha emission and the Lyman break from a z = 7.512+/- 0.004 galaxy, observed in the Faint Infrared Grism Survey (FIGS). These spectra, taken with G102 grism on Hubble Space Telescope (HST), show a significant emission line detection (6 sigma) in multiple observational position angles (PA), with total integrated Ly{\alpha} line flux of 1.06+/- 0.12 e10-17erg s-1cm-2. The line flux is nearly a factor of four higher than the previous MOSFIRE spectroscopic observations of faint Ly{\alpha} emission at {\lambda} = 1.0347{\mu}m, yielding z = 7.5078+/- 0.0004. This is consistent with other recent observations implying that ground-based near-infrared spectroscopy underestimates total emission line fluxes, and if confirmed, can have strong implications for reionization studies that are based on ground-based Lyman-{\alpha} measurements. A 4-{\sigma} detection of the NV line in one PA also suggests a weak Active Galactic Nucleus (AGN), potentially making this source the highest-redshift AGN yet found. Thus, this observation from the Hubble Space Telescope clearly demonstrates the sensitivity of the FIGS survey, and the capability of grism spectroscopy to study the epoch of reionization., Comment: Published in ApJL; matches published version