The VANDELS survey: Dust attenuation in star-forming galaxies at z=3−4\mathbf{z=3-4}

We present the results of a new study of dust attenuation at redshifts $3 < z < 4$ based on a sample of $236$ star-forming galaxies from the VANDELS spectroscopic survey. Motivated by results from the First Billion Years (FiBY) simulation project, we argue that the intrinsic spectral energy distributions (SEDs) of star-forming galaxies at these redshifts have a self-similar shape across the mass range $8.2 \leq$ log$(M_{\star}/M_{\odot}) \leq 10.6$ probed by our sample. Using FiBY data, we construct a set of intrinsic SED templates which incorporate both detailed star formation and chemical abundance histories, and a variety of stellar population synthesis (SPS) model assumptions. With this set of intrinsic SEDs, we present a novel approach for directly recovering the shape and normalization of the dust attenuation curve. We find, across all of the intrinsic templates considered, that the average attenuation curve for star-forming galaxies at $z\simeq3.5$ is similar in shape to the commonly-adopted Calzetti starburst law, with an average total-to-selective attenuation ratio of $R_{V}=4.18\pm0.29$. We show that the optical attenuation ($A_V$) versus stellar mass ($M_{\star}$) relation predicted using our method is consistent with recent ALMA observations of galaxies at $2<z<3$ in the \emph{Hubble} \emph{Ultra} \emph{Deep} \emph{Field} (HUDF), as well as empirical $A_V - M_{\star}$ relations predicted by a Calzetti-like law. Our results, combined with other literature data, suggest that the $A_V - M_{\star}$ relation does not evolve over the redshift range $0<z<5$, at least for galaxies with log$(M_{\star}/M_{\odot}) \gtrsim 9.5$. Finally, we present tentative evidence which suggests that the attenuation curve may become steeper at log$(M_{\star}/M_{\odot}) \lesssim 9.0$.Comment: 16 pages, 12 figures, accepted for publication in MNRA

Constraining Lyman-alpha spatial offsets at 3<z<5.53<z<5.5 from VANDELS slit spectroscopy

We constrain the distribution of spatially offset Lyman-alpha emission (Ly$\alpha$) relative to rest-frame ultraviolet emission in $\sim300$ high redshift ($3<z<5.5$) Lyman-break galaxies (LBGs) exhibiting Ly$\alpha$ emission from VANDELS, a VLT/VIMOS slit-spectroscopic survey of the CANDELS Ultra Deep Survey and Chandra Deep Field South fields (${\simeq0.2}~\mathrm{deg}^2$ total). Because slit spectroscopy compresses two-dimensional spatial information into one spatial dimension, we use Bayesian inference to recover the underlying Ly$\alpha$ spatial offset distribution. We model the distribution using a 2D circular Gaussian, defined by a single parameter $\sigma_{r,\mathrm{Ly}\alpha}$, the standard deviation expressed in polar coordinates. Over the entire redshift range of our sample ($3<z<5.5$), we find $\sigma_{r,\mathrm{Ly}\alpha}=1.70^{+0.09}_{-0.08}$ kpc ($68\%$ conf.), corresponding to $\sim0.25$ arcsec at $\langle z\rangle=4.5$. We also find that $\sigma_{r,\mathrm{Ly}\alpha}$ decreases significantly with redshift. Because Ly$\alpha$ spatial offsets can cause slit-losses, the decrease in $\sigma_{r,\mathrm{Ly}\alpha}$ with redshift can partially explain the increase in the fraction of Ly$\alpha$ emitters observed in the literature over this same interval, although uncertainties are still too large to reach a strong conclusion. If $\sigma_{r,\mathrm{Ly}\alpha}$ continues to decrease into the reionization epoch, then the decrease in Ly$\alpha$ transmission from galaxies observed during this epoch might require an even higher neutral hydrogen fraction than what is currently inferred. Conversely, if spatial offsets increase with the increasing opacity of the IGM, slit losses may explain some of the drop in Ly$\alpha$ transmission observed at $z>6$. Spatially resolved observations of Ly$\alpha$ and UV continuum at $6<z<8$ are needed to settle the issue.Comment: Submitted to MNRA

The VANDELS survey: Dust attenuation in star-forming galaxies at z=3−4\mathbf{z=3-4}

We present the results of a new study of dust attenuation at redshifts $3 < z < 4$ based on a sample of $236$ star-forming galaxies from the VANDELS spectroscopic survey. Motivated by results from the First Billion Years (FiBY) simulation project, we argue that the intrinsic spectral energy distributions (SEDs) of star-forming galaxies at these redshifts have a self-similar shape across the mass range $8.2 \leq$ log$(M_{\star}/M_{\odot}) \leq 10.6$ probed by our sample. Using FiBY data, we construct a set of intrinsic SED templates which incorporate both detailed star formation and chemical abundance histories, and a variety of stellar population synthesis (SPS) model assumptions. With this set of intrinsic SEDs, we present a novel approach for directly recovering the shape and normalization of the dust attenuation curve. We find, across all of the intrinsic templates considered, that the average attenuation curve for star-forming galaxies at $z\simeq3.5$ is similar in shape to the commonly-adopted Calzetti starburst law, with an average total-to-selective attenuation ratio of $R_{V}=4.18\pm0.29$. We show that the optical attenuation ($A_V$) versus stellar mass ($M_{\star}$) relation predicted using our method is consistent with recent ALMA observations of galaxies at $2<z<3$ in the \emph{Hubble} \emph{Ultra} \emph{Deep} \emph{Field} (HUDF), as well as empirical $A_V - M_{\star}$ relations predicted by a Calzetti-like law. Our results, combined with other literature data, suggest that the $A_V - M_{\star}$ relation does not evolve over the redshift range $0<z<5$, at least for galaxies with log$(M_{\star}/M_{\odot}) \gtrsim 9.5$. Finally, we present tentative evidence which suggests that the attenuation curve may become steeper at log$(M_{\star}/M_{\odot}) \lesssim 9.0$.Comment: 16 pages, 12 figures, accepted for publication in MNRA

The role of galaxy mass on AGN emission: a view from the VANDELS survey

We present a comparative analysis of the properties of AGN emitting at radio and X-ray wavelengths. The study is performed on 907 X-ray AGN and 100 radio AGN selected on the CDFS and UDS fields and makes use of new and ancillary data available to the VANDELS collaboration. Our results indicate that the mass of the host galaxy is a fundamental quantity which determines the level of AGN activity at the various wavelengths. Indeed large stellar masses are found to be connected with AGN radio emission, as virtually all radio-active AGN reside within galaxies of M*>10^{10} Msun. Large stellar masses also seem to favour AGN activity in the X-ray, even though X-ray AGN present a mass distribution which is more spread out and with a non-negligible tail at M*<10^{9} Msun. Stellar mass alone is also observed to play a fundamental role in simultaneous radio and X-ray emission: the percentage of AGN active at both wavelengths increases from around 1% of all X-ray AGN residing within hosts of M*<10^{11} Msun to about 13% in more massive galaxies. In the case of radio-selected AGN, such a percentage moves from about 15% to about 45% (but up to 80% in the deepest fields). Neither cosmic epoch, nor radio luminosity, X-ray luminosity, Eddington ratio or star-formation rate of the hosts are found to be connected to an enhanced probability for joint radio+X-ray emission of AGN origin. Furthermore, only a loose relation is observed between X-ray and radio luminosity in those AGN which are simultaneously active at both frequencies.Comment: 18 pages, 18 figures, to appear on MNRA

The VANDELS survey: A strong correlation between Lyα\alpha equivalent width and stellar metallicity at 3≤z≤5\mathbf{3\leq z \leq 5}

We present the results of a new study investigating the relationship between observed Ly$\alpha$ equivalent width ($W_{\lambda}$(Ly$\alpha$)) and the metallicity of the ionizing stellar population ($Z_{\star}$) for a sample of $768$ star-forming galaxies at $3 \leq z \leq 5$ drawn from the VANDELS survey. Dividing our sample into quartiles of rest-frame $W_{\lambda}$(Ly$\alpha$) across the range -58 \unicode{xC5} \lesssim $W_{\lambda}$(Ly$\alpha$) \lesssim 110 \unicode{xC5} we determine $Z_{\star}$ from full spectral fitting of composite far-ultraviolet (FUV) spectra and find a clear anti-correlation between $W_{\lambda}$(Ly$\alpha$) and $Z_{\star}$. Our results indicate that $Z_{\star}$ decreases by a factor $\gtrsim 3$ between the lowest $W_{\lambda}$(Ly$\alpha$) quartile ($\langle$$W_{\lambda}$(Ly$\alpha$)\rangle=-18\unicode{xC5}) and the highest $W_{\lambda}$(Ly$\alpha$) quartile ($\langle$$W_{\lambda}$(Ly$\alpha$)\rangle=24\unicode{xC5}). Similarly, galaxies typically defined as Lyman Alpha Emitters (LAEs; $W_{\lambda}$(Ly$\alpha$) >20\unicode{xC5}) are, on average, metal poor with respect to the non-LAE galaxy population ($W_{\lambda}$(Ly$\alpha$) \leq20\unicode{xC5}) with $Z_{\star}$$_{\rm{non-LAE}}\gtrsim 2 \times$ $Z_{\star}$$_{\rm{LAE}}$. Finally, based on the best-fitting stellar models, we estimate that the increasing strength of the stellar ionizing spectrum towards lower $Z_{\star}$ is responsible for $\simeq 15-25\%$ of the observed variation in $W_{\lambda}$(Ly$\alpha$) across our sample, with the remaining contribution ($\simeq 75-85\%$) being due to a decrease in the HI/dust covering fractions in low $Z_{\star}$ galaxies.Comment: 10 pages, 6 figures, MNRAS accepte

Insights into the reionization epoch from cosmic-noon-Civ emitters in the VANDELS survey

Recently, intense emission from nebular CaIII and CIV emission lines have been observed in galaxies in the epoch of reionization (z&gt;6) and have been proposed as the prime way of measuring their redshift and studying their stellar populations. These galaxies might represent the best examples of cosmic reionizers, as suggested by recent low-z observations of Lyman continuum emitting galaxies, but it is hard to directly study the production and escape of ionizing photons at such high redshifts. The ESO spectroscopic public survey VANDELS offers the unique opportunity to find rare examples of such galaxies at cosmic noon (z~3), thanks to the ultra deep observations available. We have selected a sample of 39 galaxies showing CIV emission, whose origin (after a careful comparison to photoionization models) can be ascribed to star formation and not to active galactic nuclei. By using a multiwavelength approach, we determined their physical properties including metallicity and the ionization parameter and compared them to the properties of the parent population to understand what the ingredients are that could characterize the analogs of the cosmic reionizers. We find that CIV emitters are galaxies with high photon production efficiency and there are strong indications that they might also have a large escape fraction: given the visibility of CIV in the epoch of reionization, this could become the best tool to pinpoint the cosmic reioinzers

Insights into the reionization epoch from cosmic-noon-CIV emitters in the VANDELS survey

Recently, intense emission from nebular C III] and C IV emission lines have been observed in galaxies in the epoch of reionization ($z>6$) and have been proposed as the prime way of measuring their redshift and studying their stellar populations. These galaxies might represent the best examples of cosmic reionizers, as suggested by recent low-z observations of Lyman Continuum emitting galaxies, but it is hard to directly study the production and escape of ionizing photons at such high redshifts. The ESO spectroscopic public survey VANDELS offers the unique opportunity to find rare examples of such galaxies at cosmic noon ($z\sim 3$), thanks to the ultra deep observations available. We have selected a sample of 39 galaxies showing C IV emission, whose origin (after a careful comparison to photoionization models) can be ascribed to star formation and not to AGN. By using a multi-wavelength approach, we determine their physical properties including metallicity and ionization parameter and compare them to the properties of the parent population to understand what are the ingredients that could characterize the analogs of the cosmic reionizers. We find that C IV emitters are galaxies with high photons production efficiency and there are strong indications that they might have also large escape fraction: given the visibility of C IV in the epoch of reionization this could become the best tool to pinpoint the cosmic reioinzers.Comment: 15 pages, 10 figures, resubmitted to A&A after addressing the referee comment

Timing the earliest quenching events with a robust sample of massive quiescent galaxies at 2 &lt;z &lt;5

We present a sample of 151 massive ($M_* > 10^{10}\mathrm{M_\odot}$) quiescent galaxies at $2 < z < 5$, based on a sophisticated Bayesian spectral energy distribution fitting analysis of the CANDELS UDS and GOODS-South fields. Our sample includes a robust sub-sample of 61 objects for which we confidently exclude low-redshift and star-forming solutions. We identify 10 robust objects at $z>3$, of which 2 are at $z>4$. We report formation redshifts, demonstrating that the oldest objects formed at $z > 6$, however individual ages from our photometric data have significant uncertainties, typically $\sim0.5$ Gyr. We demonstrate that the UVJ colours of the quiescent population evolve with redshift at $z>3$, becoming bluer and more similar to post-starburst galaxies at lower redshift. Based upon this we construct a model for the time-evolution of quiescent galaxy UVJ colours, concluding that the oldest objects are consistent with forming the bulk of their stellar mass at $z\sim6-7$ and quenching at $z\sim5$. We report spectroscopic redshifts for two of our objects at $z=3.440$ and $3.396$, which exhibit extremely weak Ly$\alpha$ emission in ultra-deep VANDELS spectra. We calculate star-formation rates based on these line fluxes, finding that these galaxies are consistent with our quiescent selection criteria, provided their Ly$\alpha$ escape fractions are $>3$ and $>10$ per cent respectively. We finally report that our highest-redshift robust object exhibits a continuum break at $\lambda\sim7000$A in a spectrum from VUDS, consistent with our photometric redshift of $z_\mathrm{phot}=4.72^{+0.06}_{-0.04}$. If confirmed as quiescent this object would be the highest-redshift known quiescent galaxy. To obtain stronger constraints on the times of the earliest quenching events, high-SNR spectroscopy must be extended to $z\gtrsim3$ quiescent objects.Comment: 14 pages, 8 figures, MNRAS accepte

The VANDELS survey: the role of ISM and galaxy physical properties in the escape of Lyα emission in z ∼ 3.5 star-forming galaxies⋆

International audienceAims. We wish to investigate the physical properties of a sample of Lyα 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α photons. Methods. From all the Lyα emitting galaxies in the VANDELS Data Release 2 at 3.5 z 4.5, we select a sample of 52 galaxies which also have a precise systemic redshift determination from at least one nebular emission line (HeII or CIII]). For these galaxies, we derive different physical properties (stellar mass, age, dust extinction and star formation rate) from spectral energy distribution (SED) fitting of the exquisite multi-wavelength photometry available in the VANDELS fields, using a state-of-the-art spectral modeling tool, BEAGLE and the UV β slope from the observed photometry. We characterize the Lyα emission in terms of kinematics, EW, FWHM and spatial extension and then estimate the velocity of the neutral outflowing gas. Thanks to the ultra-deep VANDELS spectra (up to 80 hours on-source integration) this can be achieved for individual galaxies, without relying on stacks. We then investigate the correlations between the Lyα properties and the other measured properties, to study how they affect the shape and intensity of Lyα emission. Results. We reproduce some of the well known correlations between Lyα EW and stellar mass, dust extinction and UV β 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α tend to be more compact both in UV and in Lyα. Our data reveal a new interesting correlation between the Lyα 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α 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α kinematics. In support to our interpretation, we observe that galaxies with high HI column densities have much more extended Lyα 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