Search

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

Context: Understanding galaxy evolution in dense environments, particularly proto-clusters, is crucial for studying mechanisms driving star formation and quenching. Aims: This study examines how two proto-cluster over-densities at 3 < z < 4 impact star formation rate (SFR), stellar mass, and morphology, focusing on quenched galaxies. Methods: We identified proto-cluster over-densities in the Chandra Deep Field South (CDFS) and Ultra Deep Survey (UDS) regions of the VANDELS survey. Using spectral energy distribution analysis, Bayesian methods (BEAGLE and BAGPIPES) helped derive best-fit parameters and U-V and V-J rest-frame colours (UVJ), classifying galaxies as quenched or star-forming based on UVJ diagrams and specific star formation rates (sSFR). TNG300 simulations aided interpretation. Results: Two of 13 proto-cluster over-densities host quenched galaxies with red U-V colours, low sSFR, and properties like massive passive galaxies. These quenched members are redder, older, more massive, and more compact. The highest-density peaks at z=3.55 and z=3.43 have dark matter halo masses consistent with proto-clusters and host AGNs, with five and three AGNs, respectively. Compared to field galaxies, these quenched members are in denser environments. TNG300 simulations suggest proto-clusters with quenched galaxies at high redshift evolve to contain more passive galaxies by z=1. Conclusions: The over-densities host massive quenched galaxies and AGNs in their densest peaks. Simulations reveal that sSFR for passive galaxies in proto-clusters was high at z=6, with median mass growth rates of 96% from z=6 to z=3. Conditions for mass assembly likely involve galaxy interactions and high gas accretion in dense environments. Black hole growth and AGN feedback appear to drive quenching at z=3, aligning with the properties of quenched galaxies observed in our study., Comment: 16 pages, 13 figures, 1 table. Published in A&A

In this work, we devise a model for estimating UV and optical line emission (i.e., CIII] $1909$A, H$\beta$, [OIII] $5007$A, H$\alpha$, [NII] $6583$A) tracing HII regions in the interstellar medium (ISM) of galaxies at $z\sim4-6$ from the ALMA Large Programme ALPINE. The aim is to investigate the impact of binary stars in the stellar population along with an abrupt quenching in the Star Formation History (SFH) on line emission. This is crucial for understanding the ISM's properties in early galaxies and identifying new star formation tracers in high-$z$ galaxies. The model simulates HII+Photodissociation Region (PDR) complexes through radiative transfer in 1D slabs, characterized by gas density ($n$), ionisation parameter ($U$), and metallicity ($Z$). It considers: (a) heating from star formation (SF), simulated with Starburst99 and BPASS to quantify binary stars impact; (b) constant, exponentially declining, and quenched SFH scenarios. For each galaxy, we select theoretical ratios from CLOUDY models between [CII] line emission, tracing PDRs, and nebular lines from HII regions, using these to derive expected optical/UV lines from observed [CII]. We find binary stars strongly impact line emission post-quenching, keeping UV photon flux higher for longer, maintaining free electron temperature and ionised column density in HII regions up to 5 Myr after quenching. We constrain ISM properties of our subsample, finding a low ionisation parameter $\log U{\approx}-3.8\pm 0.2$ and moderate/high densities $\log(n/\rm cm^{-3}){\approx}2.9\pm 0.6$. Finally, we derive UV/optical line luminosities-SFR relations for different burstiness parameters ($k_s$). In the fiducial BPASS model, relations show negligible SFH dependence but depend on $k_s$, while in the SB99 case, dependence is on SFH. We propose their use for characterising the burstiness of high-$z$ galaxies., Comment: 14 pages, 8 figures. Accepted for publication in A&A

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: 14 pages, 11 figures. Submitted to A&A

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

Galaxy protoclusters are sites of rapid growth, with a high density of massive galaxies driving elevated rates of star formation and accretion onto supermassive black holes. Here, we present new JWST/NIRSpec IFU observations of the J1000+0234 group at $z=4.54$, a dense region of a protocluster hosting a massive, dusty star forming galaxy (DSFG). The new data reveal two extended, high-equivalent-width (EW$_0>1000\r{A}$) [O III] nebulae that appear at both sides of the DSFG along its minor axis (namely O3-N and O3-S). On one hand, the spectrum of O3-N shows a broad and blueshifted component with a full width at half maximum (FWHM) of 1300 km/s, suggesting an outflow origin. On the other hand, O3-S stretches over 8.6 kpc, and has a velocity gradient that spans 800 km/s, but shows no evidence of a broad component. However, both sources seem to be powered by an active galactic nucleus (AGN), so we classified them as extended emission-line regions (EELRs). The strongest evidence comes from the detection of the high-ionization [Ne V] $\lambda 3427$ line toward O3-N, which paired with the lack of hard X-rays implies an obscuring column density above the Compton-thick regime. The [Ne V] line is not detected in O3-S, but we measure a He II $\lambda 4687$/H$\beta$=0.25, which is well above the expectation for star formation. Despite the remarkable alignment of O3-N and O3-S with two radio sources, we do not find evidence of shocks from a radio jet that could be powering the EELRs. We interpret this as O3-S being externally irradiated by the AGN, akin to the famous Hanny's Voorwerp object in the local Universe. In addition, classical line ratio diagnostics (e.g., [O III]/H$\beta$ vs [N II]/H$\alpha$) put the DSFG itself in the AGN region of the diagrams, and therefore suggest it to be the most probable AGN host. These results showcase the ability of JWST to unveil obscured AGN at high redshifts., Comment: 13 pages, 10 figures, and one table. Accepted for publication in A&A on November 12th 2024

Dusty star-forming galaxies (DSFGs) contribute significantly to the stellar buildup at cosmic noon. Major mergers and gas accretion are often invoked to explain DSFGs' prodigious star-formation rates (SFRs) and large stellar masses. We conducted a spatially-resolved morphological analysis of the rest-frame UV/NIR emission in three DSFGs at z~2.5. Initially discovered as CO emitters by NOEMA observations of a bright Herschel source, we observed them with the JWST/NIRCam as part of the PEARLS program. The NIRCam data reveal the galaxies' stellar populations and dust distributions on scales of 250 pc. Spatial variations in stellar mass, SFR, and dust extinction are determined in resolved maps obtained through pixel-based SED fitting. The CO emitters are massive, dusty starburst galaxies with SFRs=340-2500 Msun/yr, positioning them among the most active SFGs at 2

The ultraviolet continuum traces young stars while the near-infrared unveils older stellar populations and dust-obscured regions. Balmer emission lines provide insights on gas properties and young stellar objects but are highly affected by dust attenuation. The near-infrared Paschen lines suffer less dust attenuation and can be used to measure star formation rates (SFRs) in star-forming regions obscured by dust clouds. We select 13 sources between redshifts 1 and 3 observed with HST, JWST/NIRCam and NIRSpec based on the availability of at least one Balmer and one Paschen line with S/N > 5. With a newly-developed version of CIGALE, we fit their hydrogen line equivalent widths (EWs) and photometric data. We assess the impacts of the removal of spectroscopic data by comparing the quality of the fits of the spectro-photometric data to those with photometric data only. We compare the single (BC03) vs binary (BPASS) stellar populations models in the fitting process of spectro-photometric data. We derive the differential attenuation and explore different attenuation recipes by fitting spectro-photometric data with BC03. For each stellar model and for each input dataset (with and without EWs), we quantify the deviation on the SFRs and stellar masses from the "standard" choice. On average, the SFRs are overestimated and the stellar masses are underestimated when EWs are not included as input data. We find a major contribution of the H${\alpha}$ emission line to the broadband photometric measurements of our sources, and a trend of increasing contribution with specific SFR. Using the BPASS models has a significant impact on the derived SFRs and stellar masses. We show that a flexible attenuation recipe provides more accurate estimates of the dust attenuation parameters, especially the differential attenuation which agrees with the original value of Charlot & Fall (2000)., Comment: Accepted in Astronomy & Astrophysics, 19 pages, 15 figures

We present reduced images and catalogues of photometric and emission line data ($\sim$230,000 and $\sim$8,000 sources, respectively) for the WFC3 Infrared Spectroscopic Parallel (WISP) Survey. These data are made publicly available on the Mikulski Archive for Space Telescopes (MAST) and include reduced images from various facilities: ground-based $ugri$, HST WFC3, and Spitzer IRAC (Infrared Array Camera). Coverage in at least one additional filter beyond the WFC3/IR data are available for roughly half of the fields (227 out of 483), with $\sim$20% (86) having coverage in six or more filters from $u$-band to IRAC 3.6$\mu$m (0.35-3.6$\mu$m). For the lower spatial resolution (and shallower) ground-based and IRAC data, we perform PSF-matched, prior-based, deconfusion photometry (i.e., forced-photometry) using the TPHOT software to optimally extract measurements or upper limits. We present the methodology and software used for the WISP emission line detection and visual inspection. The former adopts a continuous wavelet transformation that significantly reduces the number of spurious sources as candidates before the visual inspection stage. We combine both WISP catalogues and perform SED fitting on galaxies with reliable spectroscopic redshifts and multi-band photometry to measure their stellar masses. We stack WISP spectra as functions of stellar mass and redshift and measure average emission line fluxes and ratios. We find that WISP emission line sources are typically `normal' star-forming galaxies based on the Mass-Excitation diagram ([OIII]/H$\beta$ vs. $M_\star$; $0.74

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

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 Kennicutt-Schmidt (KS) relation between the gas and the star formation rate (SFR) surface density ($\Sigma_{\rm gas}$-$\Sigma_{\rm SFR}$) is essential to understand star formation processes in galaxies. So far, it has been measured up to z~2.5 in main-sequence galaxies. In this letter, we aim to put constraints at z~4.5 using a sample of four massive main-sequence galaxies observed by ALMA at high resolution. We obtained ~0.3"-resolution [CII] and continuum maps of our objects, which we then converted into gas and obscured SFR surface density maps. In addition, we produced unobscured SFR surface density maps by convolving Hubble ancillary data in the rest-frame UV. We then derived the average $\Sigma_{\rm SFR}$ in various $\Sigma_{\rm gas}$ bins, and estimated the uncertainties using a Monte Carlo sampling. Our galaxy sample follows the KS relation measured in main-sequence galaxies at lower redshift and is slightly lower than predictions from simulations. Our data points probe the high end both in terms of $\Sigma_{\rm gas}$ and $\Sigma_{\rm gas}$, and gas depletion timescales (285-843 Myr) remain similar to z~2 objects. However, three of our objects are clearly morphologically disturbed, and we could have expected shorter gas depletion timescales (~100 Myr) similar to merger-driven starbursts at lower redshifts. This suggests that the mechanisms triggering starbursts at high redshift may be different than in the low- and intermediate-z Universe., Comment: 9 pages, 7 figures, 2 tables, accepted by A&A (letter)

VANDELS is a deep spectroscopic survey, performed with the VIMOS instrument at VLT, aimed at studying in detail the physical properties of high-redshift galaxies. VANDELS targeted about 2100 sources at 1

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

Supernova (SN) H0pe was discovered as a new transient in James Webb Space Telescope (JWST) NIRCam images of the galaxy cluster PLCK G165.7+67.0 taken as part of the "Prime Extragalactic Areas for Reionization and Lensing Science" (PEARLS) JWST GTO program (# 1176) on 2023 March 30 (AstroNote 2023-96; Frye et al. 2023). The transient is a compact source associated with a background galaxy that is stretched and triply-imaged by the cluster's strong gravitational lensing. This paper reports spectra in the 950-1370 nm observer frame of two of the galaxy's images obtained with Large Binocular Telescope (LBT) Utility Camera in the Infrared (LUCI) in longslit mode two weeks after the \JWST\ observations. The individual average spectra show the [OII] doublet and the Balmer and 4000 Angstrom breaks at redshift z=1.783+/-0.002. The CIGALE best-fit model of the spectral energy distribution indicates that SN H0pe's host galaxy is massive (Mstar~6x10^10 Msun after correcting for a magnification factor ~7) with a predominant intermediate age (~2 Gyr) stellar population, moderate extinction, and a magnification-corrected star formation rate ~13 Msun/yr, consistent with being below the main sequence of star formation. These properties suggest that H0pe might be a type Ia SN. Additional observations of SN H0pe and its host recently carried out with JWST (JWST-DD-4446; PI: B. Frye) will be able to both determine the SN classification and confirm its association with the galaxy analyzed in this work., Comment: 6 pages, 4 figures, Letter accepted for publication in Astronomy & Astrophysics

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

To better understand the ionizing properties of galaxies in the EoR, we investigate deep, rest-frame ultraviolet (UV) spectra of $\simeq 500$ star-forming galaxies at $3 \leq z \leq 5$ selected from the public ESO-VANDELS spectroscopic survey. The absolute ionizing photon escape fraction ($f_{\rm esc}^{\rm abs}$) is derived by combining absorption line measurements with estimates of the UV attenuation. The ionizing production efficiency ($\xi_{ion}$) is calculated by fitting the far-UV (FUV) stellar continuum of the VANDELS galaxies. We find that the $f_{\rm esc}^{\rm abs}$ and $\xi_{ion}$ parameters increase towards low-mass, blue UV-continuum slopes and strong Ly$\alpha$ emitting galaxies, and both are just slightly higher-than-average for the UV-faintest galaxies in the sample. Potential Lyman Continuum Emitters (LCEs) and selected Lyman Alpha Emitters (LAEs) show systematically higher $\xi_{ion}$ ($\log \xi_{ion}$ (Hz\erg) $\approx 25.38, 25.41$) than non-LCEs and non-LAEs ($\log \xi_{ion}$ (Hz\erg) $\approx 25.18, 25.14$) at similar UV magnitudes. This indicates very young underlying stellar populations ($\approx 10~{\rm Myr}$) at relatively low metallicities ($\approx 0.2~{\rm Z_{\odot}}$). The FUV non-ionizing spectra of potential LCEs is characterized by very blue UV slopes ($\leq -2$), enhanced Ly$\alpha$ emission ($\leq -25$A), strong UV nebular lines (e.g., high CIV1550/CIII]1908 $\geq 0.75$ ratios), and weak absorption lines ($\leq 1$A). The latter suggests the existence of low gas-column-density channels in the interstellar medium which enables the escape of ionizing photons. By comparing our VANDELS results against other surveys in the literature, our findings imply that the ionizing budget in the EoR was likely dominated by UV-faint, low-mass and dustless galaxies., Comment: 25 pages, 18 figures, 2 tables; accepted for publication in MNRAS on April 24, 2023

Gaseous flows inside and outside galaxies are key to understanding galaxy evolution, as they regulate their star formation activity across cosmic time. We study the ISM kinematics of 330 CIII or HeII emitters, using far-UV ISM absorption lines detected in the VANDELS spectra. These galaxies span a broad range of stellar masses M$_\ast$ from $10^8$ to $10^{11}$ M$_\odot$, and SFRs from 1 to 500 M$_\odot$/yr, in the redshift range between 2 and 5. We find that the bulk ISM velocity v$_{ism}$ is globally in outflow, with v$_{ism}$ of -60 $\pm$ 10 km/s for low ionization gas traced by SiII 1260 Angstrom, CII 1334, SiII 1526, and AlII 1670, and v$_{ism}$ of -160 $\pm$ 30 and -170 $\pm$ 30 km/s for higher ionization gas traced respectively by AlIII 1854-1862 and SiIV 1393-1402. Interestingly, BPASS models are able to better reproduce the stellar continuum around the SiIV doublet than other stellar population templates. For individual galaxies, $34\%$ of the sample has a positive ISM velocity shift, almost double the fraction reported at lower redshifts. Comparing v$_{ism}$ to the host galaxies properties, we find no significant correlations with M$_\ast$ or SFR, and only a marginally significant dependence (at $\sim 2\sigma$) on morphology-related parameters, with slightly higher velocities in galaxies of smaller size (probed by the equivalent radius), higher concentration, and higher SFR surface density. The outflows are consistent with models of accelerating, momentum-driven winds, with densities decreasing towards the outskirts. Our moderately lower ISM velocities compared to those found in similar studies at lower redshifts suggest that inflows and internal turbulence might play an increased role at $z>2$. We estimate mass outflow rates comparable to the SFRs of the galaxies, and an average escape velocity of 625 km/s, suggesting that most of the ISM will remain bound to the galaxy halo., Comment: 24 pages, 2 tables and 20 figures, accepted for publication in A&A. This replaces the submitted version posted before. Fig. A.2 replaced, Fig. A.1 updated, typos corrected. Language editing by A&A has been performed

The circumgalactic medium (CGM) is the location where the interplay between large-scale outflows and accretion onto galaxies occurs. Metals in different ionization states flowing between the circumgalactic and intergalactic mediums are affected by large galactic outflows and low-ionization state inflowing gas. Observational studies on their spatial distribution and their relation with galaxy properties may provide important constraints on models of galaxy formation and evolution. To provide new insights into the spatial distribution of the circumgalactic of star-forming galaxies, we select a sample of 238 close pairs at $1.5 < z <4.5$ ($\langle z\rangle\sim$2.6) from the VIMOS Ultra Deep Survey. We then generate composite spectra by co-adding spectra of $background$ galaxies that provide different sight-lines across the CGM to examine the spatial distribution of the gas located around these galaxies and investigate possible correlations between the strength of the low- and high-ionization absorption features with different galaxy properties. We detect C II, Si II, Si IV and C IV) up to separations $\langle b \rangle=$ 172 kpc and 146 kpc. Our $W_{0}$ radial profiles suggest a potential redshift evolution for the CGM gas content producing these absorptions. We find a correlation between C II and C IV with star formation rate, stellar mass and trends with galaxy size estimated by the effective radius and azimuthal angle. Galaxies with high star formation rate show stronger C IV absorptions compared with star-forming galaxies with low SFR and low stellar mass. These results could be explained by stronger outflows, softer radiation fields unable to ionize high-ionization state lines or by the galactic fountain scenario where metal-rich gas ejected from previous star-formation episodes fall back to the galaxy., Comment: Accepted for publication in A&A

We present the average [CII] $158\,\rm{\mu m}$ emission line sizes of UV-bright star-forming galaxies at $z\sim7$. Our results are derived from a stacking analysis of [CII] $158\,\rm{\mu m}$ emission lines and dust continua observed by ALMA, taking advantage of the large program Reionization Era Bright Emission Line Survey (REBELS). We find that the average [CII] emission at $z\sim7$ has an effective radius $r_e$ of $2.2\pm0.2\,\rm{kpc}$. It is $\gtrsim2\times$ larger than the dust continuum and the rest-frame UV emission, in agreement with recently reported measurements for $z\lesssim6$ galaxies. Additionally, we compared the average [CII] size with $4

In the local universe, galaxies in clusters show different properties compared to more isolated systems. Understanding how this difference originates and whether it is already in place at high redshift is still a matter of debate. Thanks to uniquely deep optical spectra from the VANDELS survey, we investigate environmental effects on the stellar mass-metallicity relation (MZR) for a sample of ~1000 star-forming galaxies at 2

The luminosity functions at z < 4 - 5 suggest that most galaxies have a relatively low stellar mass (logM_star = 10) and a low dust attenuation (A_FUV = 1.0). The physical properties of these objects are quite homogeneous. We used an approach where we combined their rest-frame far-infrared and submillimeter emissions and utilized the universe and the redshift as a spectrograph to increase the amount of information in a collective way. From a subsample of 27 ALMA-detected galaxies at z > 4.5, we built an infrared spectral energy distribution composite template. It was used to fit, with CIGALE, the 105 galaxies (detections and upper limits) in the sample from the FUV to the FIR. The derived physical parameters provide information to decipher the nature of the dust cycle and of the stellar populations in these galaxies. The derived IR composite template is consistent with the galaxies in the studied sample. A delayed star formation history with tau_main = 500 Myrs is slightly favored by the statistical analysis as compared to a delayed with a final burst or a continuous star formation history. The position of the sample in the star formation rate (SFR)- M_star diagram is consistent with previous papers. The redshift evolution of the log M_star versus A_FUV relation is in agreement with evolution in the redshift of this relation. This evolution is necessary to explain the cosmic evolution of the average dust attenuation of galaxies. Evolution is also observed in the L_dust/ L_FUV (IRX) versus UV slope beta_FUV diagram: younger galaxies have bluer beta_FUV. We modeled the shift of galaxies in the IRX versus the beta_FUV diagram with the mass-weighted age as a free parameter, and we provide an equation to make predictions., Comment: 40 pages, 13 figures

We present a study designed to measure the average LyC escape fraction ($\langle f_{\rm esc}\rangle$) of star-forming galaxies at z=3.5. We assemble a sample of 148 galaxies from the VANDELS survey at $3.35\leq z_{\rm spec}\leq3.95$, selected to minimize line-of-sight contamination of their photometry. For this sample, we use ultra-deep, ground-based, $U-$band imaging and HST $V-$band imaging to robustly measure the distribution of $\mathcal{R_{\rm obs}}$ $=(L_{\rm LyC}/L_{\rm UV})_{\rm obs}$. We then model the distribution as a function of $\langle f_{\rm esc}\rangle$, carefully accounting for attenuation by dust, and the IGM (and CGM). A maximum likelihood fit to the $\mathcal{R_{\rm obs}}$ distribution returns a best-fitting value of $\langle f_{\rm esc}\rangle =0.07\pm0.02$, a result confirmed using an alternative Bayesian inference technique (both exclude $\langle f_{\rm esc}\rangle=0.0$ at $> 3\sigma$). By splitting our sample in two, we find evidence that $\langle f_{\rm esc}\rangle$ is positively correlated with Ly$\alpha$ equivalent width, with high and low sub-samples returning best fits of $\langle f_{\rm esc}\rangle=0.12^{+0.06}_{-0.04}$ and $\langle f_{\rm esc} \rangle=0.02^{+0.02}_{-0.01}$, respectively. In contrast, we find evidence that $\langle f_{\rm esc}\rangle$ is anti-correlated with intrinsic UV luminosity and UV dust attenuation; with low UV luminosity and dust attenuation sub-samples returning best fits in the range $0.10 \leq \langle f_{\rm esc}\rangle \leq 0.22$. We do not find evidence for a clear correlation between $f_{\rm esc}$ and galaxy stellar mass, suggesting it is not a primary indicator of leakage. Although larger samples are needed to further explore these trends, they suggest that it is entirely plausible that the low dust and metallicity galaxies found at z > 6 will display the $\langle f_{\rm esc}\rangle\geq0.1$ required to drive reionization., Comment: 15 pages (+ appendix), 12 figures, 1 table, Accepted for publication in Monthly Notices of the Royal Astronomical Society

The [CII] 158$~\mu$m emission line represents so far one of the most profitable tools for the investigation of the high-redshift galaxies in the early Universe. Being one of the brightest cooling lines in the rest-frame far-infrared regime of star-forming galaxies, it has been successfully exploited as a tracer of star-formation rate (SFR) in local sources. The picture is more complex at higher redshifts, where its usability in this context is still under investigation. Recent results from the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey suggest that there is no (or weak) evolution of the L$\mathrm{_{[CII]}}$-SFR relation up to $z\sim6$ but their reliability is hampered by the presence of a large population of [CII] non-detected galaxies. In this work, we characterize the population of [CII] non-detections in ALPINE. By stacking their ALMA spectra, we obtain a signal detected at $\sim5.1\sigma$, resulting in a [CII] luminosity of $\mathrm{log(L_\mathrm{[CII]}}/\mathrm{L_{\odot}})$ $\sim7.8$. When combining this value with those from the [CII] detections, we find a $\mathrm{L_{[CII]}}$-SFR relation with a slope $b=1.14\pm0.11$, in agreement within the uncertainties both with the linear relation found in the local Universe, and with the previous findings from ALPINE at $z\sim5$. This suggests that the [CII] line can be considered a good tracer of star formation up to the distant Universe. Finally, we show that the galaxies of our sample that most deviate from the observed L$_\mathrm{[CII]}$-SFR relation could suffer from a less precise redshift estimation, perhaps artificially reducing their [CII] luminosity. In this respect, we claim that there is no evidence in favour of a deficit of [CII] content in high-z galaxies, in contrast with earlier studies., Comment: Accepted for publication in A&A. 10 pages, 5 figures

We present $10$ main-sequence ALPINE galaxies (log($M/M_{\odot}$) = 9.2-11.1 and ${\rm SFR}=23-190\,{\rm M_{\odot}\,yr^{-1}}$) at $z\sim4.5$ with optical [OII] measurements from Keck/MOSFIRE spectroscopy and Subaru/MOIRCS narrow-band imaging. This is the largest such multi-wavelength sample at these redshifts, combining various measurements in the ultra-violet, optical, and far-infrared including [CII]$_{158{\rm \mu m}}$ line emission and dust continuum from ALMA and H$\alpha$ emission from Spitzer photometry. For the first time, this unique sample allows us to analyze the relation between [OII] and total star-formation rate (SFR) and the interstellar medium (ISM) properties via [OII]/[CII] and [OII]/\halpha luminosity ratios at $z\sim4.5$. The [OII]$-$SFR relation at $z\sim4.5$ cannot be described using standard local descriptions, but is consistent with a metal-dependent relation assuming metallicities around $50\%$ solar. To explain the measured dust-corrected luminosity ratios of $L[OII]/L[CII] \sim 0.98^{+0.21}_{-0.22}$ and $L[OII]/LHa \sim -0.22^{+0.13}_{-0.15}$ for our sample, ionization parameters $\log(U)< -2$ and electron densities $\log(\rm n_e / {\rm [cm^{-3}]}) \sim 2.5-3$ are required. The former is consistent with galaxies at $z\sim2-3$, however lower than at $z>6$. The latter may be slightly higher than expected given the galaxies' specific SFR. The analysis of this pilot sample suggests that typical log($ M/M_{\odot})$ > 9 galaxies at $z\sim4.5$ to have broadly similar ISM properties as their descendants at $z\sim2$ and suggest a strong evolution of ISM properties since the Epoch of Reionization at $z>6$., Comment: 14 pages, 8 figures, Accepted for publication to MNRAS

We present Lyman continuum (LyC) radiation escape fraction $f_{\rm{esc}}$ measurements for 183 spectroscopically confirmed star-forming galaxies in the redshift range $3.11 < z < 3.53$ in the \textit{Chandra} Deep Field South. We use ground-based imaging to measure $f_{\rm{esc}}$, and use ground- and space-based photometry to derive galaxy physical properties using spectral energy distribution (SED) fitting. We additionally derive [O III]+H$\beta$ equivalent widths (that fall in the observed K band) by including nebular emission in the SED fitting. After removing foreground contaminants, we report the discovery of 11 new candidate LyC leakers, with absolute LyC escape fractions, $f_{\rm{esc}}$ in the range $0.14-0.85$. From non-detections, we place $1\sigma$ upper limits of $f_{\rm{esc}}<0.12$, where the Lyman-break selected galaxies have $f_{\rm{esc}} < 0.11$ and `blindly' discovered galaxies with no prior photometric selection have $f_{\rm{esc}}<0.13$. We find a slightly higher $1\sigma$ limit of $f_{\rm{esc}}<0.20$ for extreme emission line galaxies with rest-frame [O III]+H$\beta$ equivalent widths $>300$A. For candidate LyC leakers, we find a weak negative correlation between $f_{\rm{esc}}$ and galaxy stellar masses, no correlation between $f_{\rm{esc}}$ specific star-formation rates (sSFRs) and a positive correlation between $f_{\rm{esc}}$ and EW$_0$([O III]+H$\beta$). The weak/no correlations between stellar mass and sSFRs may be explained by misaligned viewing angles and/or non-coincident timescales of starburst activity and periods of high $f_{\rm{esc}}$. Alternatively, escaping radiation may predominantly occur in highly localised star-forming regions, or $f_{\rm{esc}}$ measurements may be impacted by stochasticity of the intervening neutral medium, obscuring any global trends with galaxy properties. These hypotheses have important consequences for models of reionisation., Comment: 19 pages, 11 figures, 1 appendix, published in Monthly Notices of the Royal Astronomical Society

We present a rest-frame UV-optical stacked spectrum representative of massive quiescent galaxies at $1.010.8$. The stack is constructed using VANDELS survey data, combined with new KMOS observations. We apply two independent full-spectral-fitting approaches, measuring a total metallicity, [Z/H]=$-0.13\pm0.08$ with Bagpipes, and [Z/H]=$0.04\pm0.14$ with Alf, a fall of $\sim0.2-0.3$ dex compared with the local Universe. We also measure an iron abundance, [Fe/H] =$-0.18\pm0.08$, a fall of $\sim0.15$ dex compared with the the local Universe. We measure the alpha enhancement via the magnesium abundance, obtaining [Mg/Fe]=$0.23\pm$0.12, consistent with similar-mass galaxies in the local Universe, indicating no evolution in the average alpha enhancement of log$(M_*/\rm{M_\odot})=11$ quiescent galaxies over the last $\sim8$ Gyr. This suggests the very high alpha enhancements recently reported for several bright $z\sim1-2$ quiescent galaxies are due to their extreme masses, log$(M_*/\rm{M_\odot})\gtrsim11.5$, rather than being typical of the $z\gtrsim1$ population. The metallicity evolution we observe with redshift (falling [Z/H], [Fe/H], constant [Mg/Fe]) is consistent with recent studies. We recover a mean stellar age of $2.5^{+0.6}_{-0.4}$ Gyr, corresponding to a formation redshift, $z_\rm{form}=2.4^{+0.6}_{-0.3}$. Recent studies have obtained varying average formation redshifts for $z\gtrsim1$ massive quiescent galaxies, and, as these studies report consistent metallicities, we identify different star-formation-history models as the most likely cause. Larger spectroscopic samples from upcoming ground-based instruments will provide precise constraints on ages and metallicities at $z\gtrsim1$. Combining these with precise JWST $z>2$ quiescent-galaxy stellar-mass functions will provide an independent test of formation redshifts derived from spectral fitting., Comment: 16 pages, 3 figures, accepted for publication in ApJ

Galaxy mergers are thought to be one of the main mechanisms of the mass assembly of galaxies. Recently, many works have suggested a possible increase in the fraction of major mergers in the early Universe, reviving the debate on which processes (e.g., cold accretion, star formation, mergers) most contribute to the mass build-up of galaxies through cosmic time. To estimate the importance of major mergers in this context, we make use of the new data collected by the ALMA Large Program to INvestigate [CII] at Early times (ALPINE), which observed the [CII] 158 $\mu$m emission line from a sample of 75 main-sequence star-forming galaxies at 4.4 < z < 5.9. We used, for the first time, the morpho-kinematic information provided by the [CII] emission to obtain the fraction of major mergers ($f_{MM}$) at z~5. By adopting different prescriptions for the merger timescales ($T_{MM}$), we converted this fraction into the merger rate per galaxy ($R_{MM}$) and per volume ($\Gamma_{MM}$). We then combined our results with those at lower redshifts from the literature, computing the cosmic evolution of the merger fraction. This is described by a rapid increase from z~0 to higher redshifts, a peak at z~3, and a slow decrease towards earlier epochs. Depending on the timescale prescription used, this fraction translates into a merger rate ranging between ~0.1 and ~4.0 Gyr$^{-1}$ at z~5. Finally, we compare the specific star formation and star-formation rate density with the analogous quantities from major mergers. Our new ALPINE data reveal the presence of a significant merging activity in the early Universe. However, whether this population of mergers can provide a relevant contribution to the galaxy mass assembly at these redshifts and through the cosmic epochs is strongly dependent on the assumption of the merger timescale., Comment: Accepted for publication in section 4. Extragalactic astronomy of Astronomy and Astrophysics

We study the mean properties of a large representative sample of 217 galaxies showing CIII] emission at $2

The dust content of normal galaxies and the dust mass density (DMD) at high-z (z>4) are unconstrained given the source confusion and the sensitivity limitations of previous observations. The ALMA Large Program to INvestigate [CII] at Early Times (ALPINE), which targeted 118 UV-selected star-forming galaxies at 4.4

We use VANDELS spectroscopic data overlapping with the $\simeq$7 Ms Chandra Deep Field South survey to extend studies of high-mass X-ray binary systems (XRBs) in 301 normal star-forming galaxies in the redshift range $3 < z < 5.5$. Our analysis evaluates correlations between X-ray luminosities ($L_X$), star formation rates (SFR) and stellar metallicities ($Z_\star$) to higher redshifts and over a wider range in galaxy properties than hitherto. Using a stacking analysis performed in bins of both redshift and SFR for sources with robust spectroscopic redshifts without AGN signatures, we find convincing evolutionary trends in the ratio $L_X$/SFR to the highest redshifts probed, with a stronger trend for galaxies with lower SFRs. Combining our data with published samples at lower redshift, the evolution of $L_X$/SFR to $z\simeq5$ proceeds as $(1 + z)^{1.03 \pm 0.02}$. Using stellar metallicities derived from photospheric absorption features in our spectroscopic data, we confirm indications at lower redshifts that $L_X$/SFR is stronger for metal-poor galaxies. We use semi-analytic models to show that metallicity dependence of $L_X$/SFR alone may not be sufficient to fully explain the observed redshift evolution of X-ray emission from high-mass XRBs, particularly for galaxies with SFR $<30$ $M_\odot$ yr$^{-1}$. We speculate that the discrepancy may arise due to reduced overall stellar ages in the early Universe leading to higher $L_X$/SFR for the same metallicity. We use our data to define the redshift-dependent contribution of XRBs to the integrated X-ray luminosity density and, in comparison with models, find that the contribution of high-mass XRBs to the cosmic X-ray background at $z>6$ may be $\gtrsim 0.25$ dex higher than previously estimated., Comment: 16 pages, 10 figures, Published in MNRAS

We present ASTRODEEP-GS43, a new multiwavelength photometric catalogue of the GOODS-South field, which builds and improves upon the previously released CANDELS catalogue. We provide photometric fluxes and corresponding uncertainties in 43 optical and infrared bands (25 wide and 18 medium filters), as well as photometric redshifts and physical properties of the 34930 CANDELS $H$-detected objects, plus an additional sample of 178 $H$-dropout sources, of which 173 are $Ks$-detected and 5 IRAC-detected. We keep the CANDELS photometry in 7 bands (CTIO $U$, Hubble Space Telescope WFC3 and ISAAC-$K$), and measure from scratch the fluxes in the other 36 (VIMOS, HST ACS, HAWK-I $Ks$, Spitzer IRAC, and 23 from Subaru SuprimeCAM and Magellan-Baade Fourstar) with state-of-the-art techniques of template-fitting. We then compute new photometric redshifts with three different software tools, and take the median value as best estimate. We finally evaluate new physical parameters from SED fitting, comparing them to previously published ones. Comparing to a sample of 3931 high quality spectroscopic redshifts, for the new photo-$z$'s we obtain a normalized median absolute deviation (NMAD) of 0.015 with 3.01$\%$ of outliers (0.011, 0.22$\%$ on the bright end at $I814$<22.5), similarly to the best available published samples of photometric redshifts, such as the COSMOS UltraVISTA catalogue. The ASTRODEEP-GS43 results are in qualitative agreement with previously published catalogues of the GOODS-South field, improving on them particularly in terms of SED sampling and photometric redshift estimates. The catalogue is available for download from the Astrodeep website., Comment: 15 pages, 14 figures, accepted for publication on A&A

We present results from the NIRVANDELS survey investigating the gas-phase metallicity ($\mathrm{Z}_{\mathrm{gas}}$, tracing O/H) and stellar metallicity ($Z_{\star}$, tracing Fe/H) of 33 star-forming galaxies at redshifts $2.95 < z < 3.80$. Based on a combined analysis of deep optical and near-IR spectra, tracing the rest-frame far ultraviolet and rest-frame optical respectively, we present the first simultaneous determination of the stellar and gas-phase mass-metallicity relationships (MZRs) at $z\simeq3.4$. In both cases, we find that metallicity increases with increasing stellar mass ($M_{\star}$), and that the power-law slope at $M_{\star} \lesssim 10^{10} \mathrm{M}_{\odot}$ of both MZRs scales as $Z \propto M_{\star}^{0.3}$. Comparing the stellar and gas-phase MZRs, we present direct evidence for super-solar O/Fe ratios (i.e., $\alpha$-enhancement) at $z>3$, finding $\mathrm{(O/Fe)}\simeq (2.54 \pm 0.38) \times \mathrm{(O/Fe)}_{\odot}$, with no clear dependence on $M_{\star}$., Comment: 18 pages, 11 figures, MNRAS accepted

Aim. Our aim is to analyse the variance of the Inter-Galactic Medium transmission (IGM) by studying this parameter in the rest-frame UV spectra of a large sample of high redshift galaxies. Method. We make use of the VIMOS Ultra Deep Survey and the VANDELS public survey to have an insight into the far UV spectrum of $2.7

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

The estimate of stellar metallicities (Z*) of high-z galaxies are of paramount importance in order to understand the complexity of dust effects and the reciprocal interrelations among stellar mass, dust attenuation, stellar age, and metallicity. Benefiting from uniquely deep FUV spectra of >500 star-forming galaxies at redshifts 2

Utilizing spectroscopic observations taken for the VIMOS Ultra-Deep Survey (VUDS), new observations from Keck/DEIMOS, and publicly available observations of large samples of star-forming galaxies, we report here on the relationship between the star formation rate (SFR) and the local environment ($\delta_{gal}$) of galaxies in the early universe ($2

We want to investigate whether we can use Lyalpha emission to obtain information on the environment properties and whether Lyalpha emitters show different characteristics as a function of their environment. We estimated local densities in the VANDELS Chandra Deep Field-South (CDFS) and UKIDSS Ultra Deep Survey (UDS) fields, by using a three-dimensional algorithm which works in the RA-dec-redshift space. We selected a sample of 131 Lyalpha-emitting galaxies (EW(Lyalpha)>0 A), unbiased with respect to environmental density, to study their location with respect to the over- or under-dense environment. We identify 13 (proto)cluster candidates in the CDFS and nine in the UDS at 2

Star formation rate (SFR) measurements at z>4 have relied mostly on rest-frame far-ultraviolet (FUV) observations. The corrections for dust attenuation based on IRX-$\beta$ relation are highly uncertain and are still debated in the literature. Hence, rest-frame far-infrared (FIR) observations are necessary to constrain the dust-obscured component of the SFR. In this paper, we exploit the rest-frame FIR continuum observations collected by the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) to directly constrain the obscured SFR in galaxies at 4.4

The Lyman-$\alpha$ (Ly$\alpha$) emission line has been ubiquitously used to confirm and study high redshift galaxies. We report on the line morphology as seen in the 2D spectra from the VIMOS Ultra Deep Survey in a sample of 914 Ly$\alpha$ emitters from a parent sample of 4192 star-forming galaxies at $2

We study the projected spatial offset between the ultraviolet continuum and Ly$\alpha$ emission for 65 lensed and unlensed galaxies in the Epoch of Reionization ($5\leq z\leq7$), the first such study at these redshifts, in order to understand the potential for these offsets to confuse estimates of the Ly$\alpha$ properties of galaxies observed in slit spectroscopy. While we find that ~40% of galaxies in our sample show significant projected spatial offsets ($|\Delta_{Ly\alpha-UV}|$), we find a modest average offset of 0.61$\pm$0.08 kpc. A small fraction of our sample, ~10%, exhibits offsets of 2-4 kpc, sizes that are larger than the effective radii of typical galaxies at these redshifts. An internal comparison and a comparison to studies at lower redshift yielded no significant evidence of evolution of $|\Delta_{Ly\alpha-UV}|$ with redshift. In our own sample, UV-bright galaxies showed offsets a factor of three greater than their fainter counterparts, 0.89$\pm$0.18 vs. 0.27$\pm$0.05 kpc, respectively. We argue that offsets are likely not the result of merging processes, but are rather due to internal anisotropic processes resulting from stellar feedback facilitates Ly$\alpha$ fluorescence and/or backscattering from nearby or outflowing gas. The reduction in the Ly$\alpha$ flux due to offset effects for various observational setups was quantified through mock observations of simple simulations. It was found that the loss of Ly$\alpha$ photons for galaxies with average offsets is not, if corrected for, a limiting factor for all but the narrowest slit widths (<0.4''). However, for the largest offsets, if such offsets are mostly perpendicular to the slit major axis, slit losses were found to be extremely severe in cases where slit widths of $\leq$1'' were employed, such as those planned for James Webb Space Telescope/NIRSpec observations. (abridged), Comment: 21 pages, 7 figures, 2 tables. Updated with the accepted MNRAS version that includes minor changes to the text and two tables

We present the detailed characterisation of a sample of 56 sources serendipitously detected in ALMA band 7, as part of the ALMA Large Program to INvestigate CII at Early Times (ALPINE) in COSMOS and ECDFS. These sources have been used to derive the total infrared luminosity function (LF) and to estimate the cosmic star formation rate density (SFRD) up to z=6. We have looked for counterparts in all the available multi-wavelength and photometric redshift catalogues, and in deeper near- and mid-IR source lists and maps, to identify optically dark sources with no matches in the public catalogues. Our ALMA blind survey allows us to push further the study of the nature and evolution of dusty galaxies at high-z, identifying luminous and massive sources to redshifts and faint luminosities never probed before by any far-infrared surveys. The ALPINE data are the first ones to sample the faint-end of the infrared LF, showing little evolution from z=2.5 to z=6, and a flat slope up to the highest redshifts. The SFRD obtained by integrating the luminosity function remains almost constant between z=2 and 6, and significantly higher than the optical/UV derivations, showing an important contribution of dusty galaxies and obscured star formation up to high-z. About 14 per cent of the ALPINE serendipitous continuum sources are optically+near-IR dark (six show a counterpart only in the mid-IR and no HST or near-IR identification, while two are detected as [CII] emitters at z=5). The six HST and near-IR dark galaxies with mid-IR counterpart contribute for about 17 per cent of the total SFRD at z=5 and dominate the high-mass end of the stellar mass function at z>3., Comment: 26 pages, 18 figures. Accepted for publication in A&A

We present the [CII]158$\mu$m line luminosity functions (LFs) at $z\sim4-6$ using the ALMA observations of 118 sources, which are selected to have UV luminosity $M_{1500A}<-20.2$ and optical spectroscopic redshifts in COSMOS and ECDF-S. Of the 118 targets, 75 have significant [CII] detections and 43 are upper limits. This is by far the largest sample of [CII] detections which allows us to set constraints to the volume density of [CII] emitters at $z\sim4-6$. But because this is a UV-selected sample, we are missing [CII]-bright but UV-faint sources making our constraints strict lower limits. Our derived LFs are statistically consistent with the $z\sim0$ [CII] LF at $10^{8.25} - 10^{9.75}L_\odot$. We compare our results with the upper limits of the [CII] LF derived from serendipitous sources in the ALPINE maps (Loiacono et al. 2020). We also infer the [CII] LFs based on published far-IR and CO LFs at $z\sim4-6$. Combining our robust lower limits with these additional estimates, we set further constraints to the true number density of [CII] emitters at $z\sim 4 - 6$. These additional LF estimates are largely above our LF at $L_{[CII]}>10^9L_{\odot}$, suggesting that UV-faint but [CII]-bright sources likely make a significant contributions to the [CII] emitter volume density. When we include all the LF estimates, we find that available model predictions underestimate the number densities of [CII] emitters at $z\sim4-6$. Finally, we set a constraint on the molecular gas mass density at $z\sim4-6$, with $\rho_{mol} \sim (2-7)\times10^7M_\odot$\,Mpc$^{-3}$. This is broadly consistent with previous studies., Comment: 17 pages, 9 Figures, Update to match with the published version. Accepted for the publication in ApJ

We present ALMA observations of a merging system at z ~ 4.57, observed as a part of the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey. Combining ALMA [CII] 158 micron and far-infrared continuum data with multi-wavelength ancillary data we find that the system is composed of two massive (Mstar >~ 10^10 Msun) star-forming galaxies experiencing a major merger (stellar mass ratio r_mass ~ 0.9) at close spatial (~13 kpc; projected) and velocity (delta_v < 300 km/s) separations, and two additional faint narrow [CII]-emitting satellites. The overall system belongs to a larger-scale protocluster environment and is coincident to one of its overdensity peaks. ALMA reveals also the presence of [CII] emission arising from a circumgalactic gas structure, extending up to a diameter-scale of ~30 kpc. Our morpho-spectral decomposition analysis shows that about 50% of the total flux resides between the individual galaxy components, in a metal-enriched gaseous envelope characterized by a disturbed morphology and complex kinematics. Similarly to observations of shock-excited [CII] emitted from tidal tails in local groups, our results can be interpreted as a possible signature of interstellar gas stripped by strong gravitational interactions, with a possible contribution from material ejected by galactic outflows and emission triggered by star formation in small faint satellites. Our findings suggest that mergers could be an efficient mechanism of gas mixing in the circumgalactic medium around high-z galaxies, and thus play a key role in the galaxy baryon cycle at early epochs., Comment: 12 pages, 7 figures. Accepted for publication in A&A

The molecular gas content of normal galaxies at z>4 is poorly constrained, because the commonly used molecular gas tracers become hard to detect. We use the [CII]158um luminosity, recently proposed as a molecular gas tracer, to estimate the molecular gas content in a large sample of main-sequence star-forming galaxies at z=4.4-5.9, with a median stellar mass of 10^9.7 Msun, drawn from the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey. The agreement between molecular gas masses derived from [CII] luminosity, dynamical mass, and rest-frame 850um luminosity, extrapolated from the rest-frame 158um continuum, supports [CII] as a reliable tracer of molecular gas in our sample. We find a continuous decline of the molecular gas depletion timescale from z=0 to z=5.9, which reaches a mean value of (4.6+/-0.8)x10^8 yr at z~5.5, only a factor of 2-3 shorter than in present-day galaxies. This suggests a mild enhancement of star formation efficiency toward high redshifts, unless the molecular gas fraction significantly increases. Our estimates show that the rise in molecular gas fraction as reported previously, flattens off above z~3.7 to achieve a mean value of 63%+/-3 over z=4.4-5.9. This redshift evolution of the gas fraction is in line with the one of the specific star formation rate. We use multi-epoch abundance matching to follow the gas fraction evolution over cosmic time of progenitors of z=0 Milky Way-like galaxies in 10^13 Msun halos and of more massive z=0 galaxies in 10^14 Msun halos. Interestingly, the former progenitors show a monotonic decrease of the gas fraction with cosmic time, while the latter show a constant gas fraction from z=5.9 to z~2 and a decrease at z<2. We discuss three possible effects, namely outflows, halt of gas supplying, and over-efficient star formation, which may jointly contribute to the gas fraction plateau of the latter massive galaxies., Comment: 16 pages, 9 figures, 2 tables. Submitted to Astronomy & Astrophysics

We present dust attenuation properties of spectroscopically confirmed star forming galaxies on the main sequence at a redshift of ~4.4-5.8. Our analyses are based on the far infrared continuum observations of 118 galaxies at rest-frame 158{\mu}m obtained with the Atacama Large Millimeter Array (ALMA) Large Program to INvestigate [CII] at Early times (ALPINE). We study the connection between the ultraviolet (UV) spectral slope ($\beta$), stellar mass (M_*), and infrared excess (IRX=L_IR/L_UV). Twenty-three galaxies are individually detected in the continuum at >3.5 sigma significance. We perform a stacking analysis using both detections and nondetections to study the average dust attenuation properties at z~4.4-5.8. The individual detections and stacks show that the IRX-$\beta$ relation at z~5 is consistent with a steeper dust attenuation curve than typically found at lower redshifts (z<4). The attenuation curve is similar to or even steeper than that of the extinction curve of the Small Magellanic Cloud (SMC). This systematic change of the IRX-$\beta$ relation as a function of redshift suggests an evolution of dust attenuation properties at z>4. Similarly, we find that our galaxies have lower IRX values, up to 1dex on average, at a fixed mass compared to previously studied IRX-M_* relations at z<4, albeit with significant scatter. This implies a lower obscured fraction of star formation than at lower redshifts. Our results suggest that dust properties of UV-selected star forming galaxies at z>4 are characterised by (i) a steeper attenuation curve than at z<4, and (ii) a rapidly decreasing dust obscured fraction of star formation as a function of redshift. Nevertheless, even among this UV-selected sample, massive galaxies (log M_*/$M_\odot$>10) at z~5-6 already exhibit an obscured fraction of star formation of ~45%, indicating a rapid build-up of dust during the epoch of reionization., Comment: 12 pages, 6 figures. accepted for publication in Astronomy & Astrophysics. Minor revisions

We explore X-ray emission from a sample of 18 He II 1640 emitting star-forming galaxies at z ~ 2.3-3.6 from the VANDELS survey in the Chandra Deep Field South, to set constraints on the role of X-ray sources in powering the He II emission. We find that 4 He II emitters have tentative detections with S/N ~ 2 and have X-ray luminosities, L_X = 1.5-4.9 x 10^41 erg/s. The stacked luminosity of all 18 He II emitters is 2.6 x 10^41 erg/s, and that of a subset of 13 narrow He II emitters (FHWM(He II) < 1000 km/s) is 3.1 x 10^41 erg/s. We also measure stacked L_X for non-He II emitters through bootstrapping of matched samples, and find L_X = 2.5 x 10^41 erg/s, which is not significantly different from L_X measured for He II emitters. The L_X per star-formation rate for He II emitters (log (L_X/SFR) ~ 40.0) and non-emitters (log (L_X/SFR) ~ 39.9) are also comparable and in line with the redshift evolution and metallicity dependence predicted by models. Due to the non-significant difference between the X-ray emission from galaxies with and without He II, we conclude that X-ray binaries or weak or obscured AGNs are unlikely to be the dominant producers of He II ionising photons in VANDELS star-forming galaxies at z ~ 3. Given the comparable physical properties of both He II emitters and non-emitters reported previously, alternative He II ionising mechanisms such as localised low-metallicity stellar populations, Pop-III stars, etc. may need to be explored., Comment: 12 pages, 5 figures, accepted for publication in MNRAS