Searching for and characterising the high redshift star-forming galaxy population in the early universe

This thesis focuses on the study of LAEs and LBGs at z=3.4-6.8, a period of time between ~0.8-2 Gyr after the Big Bang, hosting the end of the reionisation epoch. To reach such deep galaxies, we have made use of the SHARDS survey of the GOODS-N field, with 200 h of exposure time with GTC. The 25 consecutive medium-band filters in the optical/NIR have provided us with very good SEDs of these galaxies, which we have further completed with ancillary data from other previous deep observations of the field, allowing a comprehensive analysis of the SEDs. In the first part of this work, we develop our own robust selection criteria, paying special attention to the interlopers rejection. The final selected sample consists of 1558 sources, divided in 528 LAEs and 1030 non-LAEs/LBGs. The comparison of this sample with previous broad-band studies has revealed that the latter are prone to introduce a non-negligible amount of low redshift interlopers due to the lack of good spectral resolution in their SEDs. From the SEDs in our sample, we have derived photometric redshifts, Ly alpha EWs, SFRs and LFs. In order to calculate ages and stellar masses of our galaxies, we have carried out a SP synthesis modeling of the SEDs, using burst-like SFHs. We focus the SP analysis on differences between subfamilies of LAEs and LBGs, as well as on the need of one or two SPs to properly describe each SED. The results derived from the model fitting, joined to the relative fraction of each subfamily found at each z in this work, support an evolutionary relation from low mass LAEs to more massive LBGs. With the model-derived masses, we have built SMFs and estimated the SMD at different redshifts. We have also performed a search for galaxies in close groups, reporting 92 of them, some with up to six members. Many of these groups show tail-like structures that point to a gravitational bound between their galaxies. Additionally, we have studied a previously discovered z=5.198 overdensity in GOODS-N, finding 44 new candidates within our sample. Some of them have been further confirmed through MOS observations. If most of the remaining new candidates are confirmed, this proto-cluster would be the most populated one up to date beyond z=5 in such a tight redshift range. The final part of this thesis is dedicated to an unsupervised ML SEDs clustering of both our sample and the SC4K-COSMOS high-z LAEs sample. The results point to the Ly alpha line and the UV slope as the most relevant parameters driving the classes differentiation, separating sources with strong Ly alpha emission and blue UV slope from those with steeper slope and fainter line emission. The study of the median ages and stellar masses of the groups supports the evolutionary hypothesis between LAEs and LBGs subgroups. Moreover, the ML classification has shown to be very effective for the detection of some few remaining interlopers in both samples

The evolution of rest-frame UV properties, Lya EWs and the SFR-Stellar mass relation at z~2-6 for SC4K LAEs

We explore deep rest-frame UV to FIR data in the COSMOS field to measure the individual spectral energy distributions (SED) of the ~4000 SC4K (Sobral et al. 2018) Lyman-alpha (Lya) emitters (LAEs) at z~2-6. We find typical stellar masses of 10$^{9.3\pm0.6}$ M$_{\odot}$ and star formation rates (SFR) of SFR$_{SED}=4.5^{+8.8}_{-2.5}$ M$_{\odot}$/yr and SFR$_{Lya}=5.9^{+6.3}_{-2.6}$ M$_{\odot}$/yr, combined with very blue UV slopes of beta=-2.0$^{+0.3}_{-0.5}$, but with significant variations within the population. M$_{UV}$ and beta are correlated in a similar way to UV-selected sources, but LAEs are consistently bluer. This suggests that LAEs are the youngest and/or most dust-poor subset of the UV-selected population. We also study the Lya rest-frame equivalent width (EW$_0$) and find 45 "extreme" LAEs with EW$_0>240$ A (3 $\sigma$), implying a low number density of $(7\pm1)\times10^{-7}$ Mpc$^{-3}$. Overall, we measure little to no evolution of the Lya EW$_0$ and scale length parameter ($w_0$) which are consistently high (EW$_0=140^{+280}_{-70}$ A, $w_0=129^{+11}_{-11}$ A) from z~6 to z~2 and below. However, $w_0$ is anti-correlated with M$_{UV}$ and stellar mass. Our results imply that sources selected as LAEs have a high Lya escape fraction (f$_{esc, Lya}$) irrespective of cosmic time, but f$_{esc, Lya}$ is still higher for UV-fainter and lower mass LAEs. The least massive LAEs (

A First Look at the Abundance Pattern -- O/H, C/O, Ne/O, and Fe/O -- in z>7z>7 Galaxies with JWST/NIRSpec

We analyze the rest-frame near-UV and optical nebular spectra of three z > 7 galaxies from the Early Release Observations taken with the Near-Infrared Spectrograph (NIRSpec) on the James Webb Space Telescope (JWST). These three high-z galaxies show the detection of several strong-emission nebular lines, including the temperature-sensitive [O III] $\lambda$4363 line, allowing us to directly determine the nebular conditions and gas-phase abundances for O/H, C/O, Ne/O, and Fe/O. We derive O/H abundances and ionization parameters that are generally consistent with other recent analyses. The lowest-mass galaxy has a large O/H uncertainty, which as a significant effect on anchoring the mass-metallicity relationship (i.e., slope) and tests of its redshift evolution. We also detect the C III] $\lambda$$\lambda$1907,1909 emission in a z > 8 galaxy from which we determine the most distant C/O abundance to date. This valuable detection provides the first test of C/O redshift evolution out to high-redshift. For neon, we use the high-ionization [Ne III] $\lambda$3869 line to measure the first Ne/O abundances at z>7, finding no evolution in this $\alpha$-element ratio. To investigate the Fe abundance, we explore the tentative detection of weak [Fe II] and [Fe III] lines in a z>8 galaxy, which would indicate a rapid build up of metals. Importantly, we demonstrate that properly flux-calibrated and higher S/N spectra are crucial to robustly determine the abundance pattern in z>7 galaxies with NIRSpec/JWST.Comment: Submitted to The Astrophysical Journal Letters. Comments welcom

The evolution of the UV luminosity and stellar mass functions of Lyman-alpha emitters from z~2 to z~6

We measure the evolution of the rest-frame UV luminosity function (LF) and the stellar mass function (SMF) of Lyman-alpha (Lya) emitters (LAEs) from z~2 to z~6 by exploring ~4000 LAEs from the SC4K sample. We find a correlation between Lya luminosity (LLya) and rest-frame UV (M_UV), with best-fit M_UV=-1.6+-0.2 log10(LLya/erg/s)+47+-12 and a shallower relation between LLya and stellar mass (Mstar), with best-fit log10( Mstar/Msun)=0.9+-0.1 log10(LLya/erg/s)-28+-4.0. An increasing LLya cut predominantly lowers the number density of faint M_UV and low Mstar LAEs. We estimate a proxy for the full UV LFs and SMFs of LAEs with simple assumptions of the faint end slope. For the UV LF, we find a brightening of the characteristic UV luminosity (M_UV*) with increasing redshift and a decrease of the characteristic number density (Phi*). For the SMF, we measure a characteristic stellar mass (Mstar*/Msun) increase with increasing redshift, and a Phi* decline. However, if we apply a uniform luminosity cut of log10 (LLya/erg/s) >= 43.0, we find much milder to no evolution in the UV and SMF of LAEs. The UV luminosity density (rho_UV) of the full sample of LAEs shows moderate evolution and the stellar mass density (rho_M) decreases, with both being always lower than the total rho_UV and rho_M of more typical galaxies but slowly approaching them with increasing redshift. Overall, our results indicate that both rho_UV and rho_M of LAEs slowly approach the measurements of continuum-selected galaxies at z>6, which suggests a key role of LAEs in the epoch of reionisation

Extreme damped Lyman-α\alpha absorption in young star-forming galaxies at z=9−11z=9-11

The onset of galaxy formation is thought to be initiated by the infall of neutral, pristine gas onto the first protogalactic halos. However, direct constraints on the abundance of neutral atomic hydrogen (HI) in galaxies have been difficult to obtain at early cosmic times. Here we present spectroscopic observations with JWST of three galaxies at redshifts $z=8.8 - 11.4$, about $400-600$ Myr after the Big Bang, that show strong damped Lyman-$\alpha$ absorption ($N_{\rm HI} > 10^{22}$ cm$^{-2}$) from HI in their local surroundings, an order of magnitude in excess of the Lyman-$\alpha$ absorption caused by the neutral intergalactic medium at these redshifts. Consequently, these early galaxies cannot be contributing significantly to reionization, at least at their current evolutionary stages. Simulations of galaxy formation show that such massive gas reservoirs surrounding young galaxies so early in the history of the universe is a signature of galaxy formation in progress.Comment: Submitte

Probing the earliest phases in the formation of massive galaxies with simulated HST+JWST imaging data from Illustris

We use the Illustris-1 simulation to explore the capabilities of the $\textit{Hubble}$ and $\textit{James Webb Space Telescope}$ data to analyze the stellar populations in high-redshift galaxies, taking advantage of the combined depth, spatial resolution, and wavelength coverage. For that purpose, we use simulated broad-band ACS, WFC3 and NIRCam data and 2-dimensional stellar population synthesis (2D-SPS) to derive the integrated star formation history (SFH) of massive (M$_{\ast}>10^{10}\,$M$_{\odot}$) simulated galaxies at $110^{11}\,$M$_{\odot}$ galaxy. In particular, we explore the potential of HST and JWST datasets reaching a depth similar to those of the CANDELS and ongoing CEERS observations, respectively, and concentrate on determining the capabilities of this dataset for characterizing the first episodes in the SFH of local M$_{\ast}>10^{11}\,$M$_{\odot}$ galaxies by studying their progenitors at $z>1$. The 2D-SPS method presented in this paper has been calibrated to robustly recover the cosmic times when the first star formation episodes occurred in massive galaxies, i.e., the first stages in their integrated SFHs. In particular, we discuss the times when the first 1% to 50% of their total stellar mass formed in the simulation. We demonstrate that we can recover these ages with typical median systematic offset of less than 5% and scatter around 20%-30%. According to our measurements on Illustris data, we are able to recover that local M$_{\ast}>10^{11}\,$M$_{\odot}$ galaxies would have started their formation by $z=16$, forming the first 5% of their stellar mass present at $z \sim 1$ by $z=4.5$, 10% by $z=3.7$, and 25% by $z=2.7$.Comment: 28 pages, 13 figures, 4 tables. ApJ in press. Summary of changes from original submission: the major change is that we now include in Sec. 6 the comparison of the results obtained for our sample of massive 1 < z < 4 progenitors with those obtained by considering all massive galaxies at 1 < z < 4 in the simulated images. Several figures and sections have been update

Mapping dusty galaxy growth at z>5z>5 with FRESCO: Detection of Hα\alpha in submm galaxy HDF850.1 and the surrounding overdense structures

We report the detection of a 13$\sigma$ H$\alpha$ emission line from HDF850.1 at $z=5.188\pm0.001$ using the FRESCO NIRCam F444W grism observations. Detection of H$\alpha$ in HDF850.1 is noteworthy, given its high far-IR luminosity, substantial dust obscuration, and the historical challenges in deriving its redshift. HDF850.1 shows a clear detection in the F444W imaging data, distributed between a northern and southern component, mirroring that seen in [CII] from the Plateau de Bure Interferometer. Modeling the SED of each component separately, we find that the northern component has a higher mass, star formation rate (SFR), and dust extinction than the southern component. The observed H$\alpha$ emission appears to arise entirely from the less-obscured southern component and shows a similar $\Delta$v$\sim$+130 km/s velocity offset to that seen for [CII] relative to the source systemic redshift. Leveraging H$\alpha$-derived redshifts from FRESCO observations, we find that HDF850.1 is forming in one of the richest environments identified to date at $z>5$, with 100 $z=5.17-5.20$ galaxies distributed across 10 structures and a $\sim$(15 cMpc)$^3$ volume. Based on the evolution of analogous structures in cosmological simulations, the $z=5.17-5.20$ structures seem likely to collapse into a single $>$10$^{14}$ $M_{\odot}$ cluster by $z\sim0$. Comparing galaxy properties forming within this overdensity with those outside, we find the masses, SFRs, and $UV$ luminosities inside the overdensity to be clearly higher. The prominence of H$\alpha$ line emission from HDF850.1 and other known highly-obscured $z>5$ galaxies illustrates the potential of NIRCam-grism programs to map both the early build-up of IR-luminous galaxies and overdense structures.Comment: Submitted to MNRAS. 20 pages, 10 figures and 8 tables (including appendices

A z=1.85 galaxy group in CEERS: evolved, dustless, massive intra-halo light and a brightest group galaxy in the making

We present CEERS JWST/NIRCam imaging of a massive galaxy group at z=1.85, to explore the early JWST view on massive group formation in the distant Universe. The group contains >16 members (including 6 spectros. confirmations) down to log10(Mstar/Msun)=8.5, including the brightest group galaxy (BGG) in the process of actively assembling at this redshift. The BGG is comprised of multiple merging components extending ~3.6" (30kpc) across the sky. The BGG contributes 69% of the group's total galactic stellar mass, with one of the merging components containing 76% of the total mass of the BGG and a SFR>1810Msun/yr. Most importantly, we detect intra-halo light (IHL) in several HST and JWST/NIRCam bands, allowing us to construct a state-of-the-art rest-frame UV-NIR Spectral Energy Distribution of the IHL for the first time at this high redshift. This allows stellar population characterisation of both the IHL and member galaxies, as well as the morphology distribution of group galaxies vs. their star-formation activity when coupled with Herschel data. We create a stacked image of the IHL, giving us a sensitivity to extended emission of 28.5 mag/arcsec2 at rest-frame 1um. We find that the IHL is extremely dust poor (Av~0), containing an evolved stellar population of log10(t50/yr)=8.8, corresponding to a formation epoch for 50% of the stellar material 0.63Gyr before z=1.85. There is no evidence of ongoing star-formation in the IHL. The IHL in this group at z=1.85 contributes ~10% of the total stellar mass, comparable with what is observed in local clusters. This suggests that the evolution of the IHL fraction is more self-similar with redshift than predicted by some models, challenging our understanding of IHL formation during the assembly of high-redshift clusters. JWST is unveiling a new side of group formation at this redshift, which will evolve into Virgo-like structures in the local Universe.Comment: 14 pages + appendix, 11 figures, 4 tables. Accepted to A&A on 15th May 202

CEERS: 7.7 μ{\mu}m PAH Star Formation Rate Calibration with JWST MIRI

We test the relationship between UV-derived star formation rates (SFRs) and the 7.7 ${\mu}$m polycyclic aromatic hydrocarbon (PAH) luminosities from the integrated emission of galaxies at z ~ 0 - 2. We utilize multi-band photometry covering 0.2 - 160 ${\mu}$m from HST, CFHT, JWST, Spitzer, and Herschel for galaxies in the Cosmic Evolution Early Release Science (CEERS) Survey. We perform spectral energy distribution (SED) modeling of these data to measure dust-corrected far-UV (FUV) luminosities, $L_{FUV}$, and UV-derived SFRs. We then fit SED models to the JWST/MIRI 7.7 - 21 ${\mu}$m CEERS data to derive rest-frame 7.7 ${\mu}$m luminosities, $L_{770}$, using the average flux density in the rest-frame MIRI F770W bandpass. We observe a correlation between $L_{770}$ and $L_{FUV}$, where log $L_{770}$ is proportional to (1.27+/-0.04) log $L_{FUV}$. $L_{770}$ diverges from this relation for galaxies at lower metallicities, lower dust obscuration, and for galaxies dominated by evolved stellar populations. We derive a "single-wavelength" SFR calibration for $L_{770}$ which has a scatter from model estimated SFRs (${{\sigma}_{{\Delta}SFR}}$) of 0.24 dex. We derive a "multi-wavelength" calibration for the linear combination of the observed FUV luminosity (uncorrected for dust) and the rest-frame 7.7 ${\mu}$m luminosity, which has a scatter of ${{\sigma}_{{\Delta}SFR}}$ = 0.21 dex. The relatively small decrease in ${\sigma}$ suggests this is near the systematic accuracy of the total SFRs using either calibration. These results demonstrate that the rest-frame 7.7 ${\mu}$m emission constrained by JWST/MIRI is a tracer of the SFR for distant galaxies to this accuracy, provided the galaxies are dominated by star-formation with moderate-to-high levels of attenuation and metallicity.Comment: 20 pages, 11 figures, 2 tables, submitted to Ap

A Milky Way-like barred spiral galaxy at a redshift of 3

International audienceThe majority of massive disk galaxies in the local Universe show a stellar barred structure in their central regions, including our Milky Way. Bars are supposed to develop in dynamically cold stellar disks at low redshift, as the strong gas turbulence typical of disk galaxies at high redshift suppresses or delays bar formation. Moreover, simulations predict bars to be almost absent beyond $z = 1.5$ in the progenitors of Milky Way-like galaxies. Here we report observations of ceers-2112, a barred spiral galaxy at redshift $z_{\rm phot} \sim 3$, which was already mature when the Universe was only 2 Gyr old. The stellar mass ($M_{\star} = 3.9 \times 10^9 M_{\odot}$) and barred morphology mean that ceers-2112 can be considered a progenitor of the Milky Way, in terms of both structure and mass-assembly history in the first 2 Gyr of the Universe, and was the closest in mass in the first 4 Gyr. We infer that baryons in galaxies could have already dominated over dark matter at $z \sim 3$, that high-redshift bars could form in approximately 400 Myr and that dynamically cold stellar disks could have been in place by redshift $z = 4-5$ (more than 12 Gyrs ago)