No Arabic abstract
We present the results of a study investigating the dust attenuation law at $zsimeq 5$, based on synthetic spectral energy distributions (SEDs) calculated for a sample of N=498 galaxies drawn from the First Billion Years (FiBY) simulation project. The simulated galaxies at $zsimeq 5$, which have M$_{1500} leq -18.0$ and $7.5 leq rm{log(M/M}_{odot}rm{)} leq 10.2$, display a mass-dependent $alpha$-enhancement, with a median value of $[alpha/rm{Fe}]_{z=5}~simeq~4~times~[alpha/rm{Fe}]_{Z_{odot}}$. The median Fe/H ratio of the simulated galaxies is $0.14pm0.05$ which, even including the effects of nebular continuum, produces steep intrinsic UV continuum slopes; $langle beta_{i} rangle = -2.4 pm 0.05$. Using a set of simple dust attenuation models, in which the wavelength-dependent attenuation is assumed to be of the form $A(lambda) propto lambda^{n}$, we explore the parameter values which best reproduce the observed $z=5$ luminosity function (LF) and colour-magnitude relation (CMR). We find that a simple model in which the absolute UV attenuation is a linearly increasing function of log stellar mass, and the dust attenuation slope ($n$) is within the range $-0.7 leq n leq-0.3$, can successfully reproduce the LF and CMR over a wide range of stellar population synthesis model (SPS) assumptions. This range of attenuation curves is consistent with a power-law fit to the Calzetti attenuation law in the UV ($n=-0.55$), and other similarly `grey star-forming galaxy attenuation curves recently derived at $zsimeq2$. In contrast, attenuation curves as steep as the Small Magellanic Cloud (SMC) extinction curve ($n=-1.24$) are formally ruled out. Finally, we show that our models are consistent with recent 1.3mm ALMA observations of the Hubble Ultra Deep Field (HUDF), and predict the form of the $zsimeq5$ IRX$-beta$ relation.
We make use of SHARDS, an ultra-deep (<26.5AB) galaxy survey that provides optical photo-spectra at resolution R~50, via medium band filters (FWHM~150A). This dataset is combined with ancillary optical and NIR fluxes to constrain the dust attenuation law in the rest-frame NUV region of star-forming galaxies within the redshift window 1.5<z<3. We focus on the NUV bump strength (B) and the total-to-selective extinction ratio (Rv), targeting a sample of 1,753 galaxies. By comparing the data with a set of population synthesis models coupled to a parametric dust attenuation law, we constrain Rv and B, as well as the colour excess, E(B-V). We find a correlation between Rv and B, that can be interpreted either as a result of the grain size distribution, or a variation of the dust geometry among galaxies. According to the former, small dust grains are associated with a stronger NUV bump. The latter would lead to a range of clumpiness in the distribution of dust within the interstellar medium of star-forming galaxies. The observed wide range of NUV bump strengths can lead to a systematic in the interpretation of the UV slope ($beta$) typically used to characterize the dust content. In this study we quantify these variations, concluding that the effects are $Deltabeta$~0.4.
A diverse range of dust attenuation laws is found in star-forming galaxies. In particular, Tress et al. (2018) studied the SHARDS survey to constrain the NUV bump strength (B) and the total-to selective ratio (Rv) of 1,753 star-forming galaxies in the GOODS-N field at 1.5<z<3. We revisit here this sample to assess the implications and possible causes of the correlation found between Rv and B. The UVJ bicolour plot and main sequence of star formation are scrutinised to look for clues into the observed trend. The standard boundary between quiescent and star-forming galaxies is preserved when taking into account the wide range of attenuation parameters. However, an additional degeneracy, regarding the effective attenuation law, is added to the standard loci of star-forming galaxies in the UVJ diagram. A simple phenomenological model with an age-dependent extinction (at fixed dust composition) is compatible with the observed trend between Rv and B, whereby the opacity decreases with the age of the populations, resulting in a weaker NUV bump when the overall attenuation is shallower (greyer). In addition, we compare the constraints obtained by the SHARDS sample with dust models from the literature, supporting a scenario where geometry could potentially drive the correlation between Rv and B
Dust attenuation in galaxies has been extensively studied nearby, however, there are still many unknowns regarding attenuation in distant galaxies. We contribute to this effort using observations of star-forming galaxies in the redshift range z = 0.05-0.15 from the DYNAMO survey. Highly star-forming DYNAMO galaxies share many similar attributes to clumpy, star-forming galaxies at high redshift. Considering integrated Sloan Digital Sky Survey observations, trends between attenuation and other galaxy properties for DYNAMO galaxies are well matched to star-forming galaxies at high redshift. Integrated gas attenuations of DYNAMO galaxies are 0.2-2.0 mags in the V-band, and the ratio of stellar E(B-V) and gas E(B-V) is 0.78-0.08 (compared to 0.44 at low redshift). Four highly star-forming DYNAMO galaxies were observed at H-alpha using the Hubble Space Telescope and at Pa-alpha using integral field spectroscopy at Keck. The latter achieve similar resolution (~0.8-1 kpc) to our HST imaging using adaptive optics, providing resolved observations of gas attenuations of these galaxies on sub-kpc scales. We find < 1.0 mag of variation in attenuation (at H-alpha) from clump to clump, with no evidence of highly attenuated star formation. Attenuations are in the range 0.3-2.2 mags in the V band, consistent with attenuations of low redshift star-forming galaxies. The small spatial variation on attenuation suggests that a majority of the star-formation activity in these four galaxies occurs in relatively unobscured regions and, thus, star-formation is well characterised by our H-alpha observations.
We present millimetre dust emission measurements of two Lyman Break Galaxies at z~3 and construct for the first time fully sampled infrared spectral energy distributions (SEDs), from mid-IR to the Rayleigh-Jeans tail, of individually detected, unlensed, UV-selected, main sequence (MS) galaxies at $z=3$. The SED modelling of the two sources confirms previous findings, based on stacked ensembles, of an increasing mean radiation field <U> with redshift, consistent with a rapidly decreasing gas metallicity in z > 2 galaxies. Complementing our study with CO[3-2] emission line observations, we measure the molecular gas mass (M_H2) reservoir of the systems using three independent approaches: 1) CO line observations, 2) the dust to gas mass ratio vs metallicity relation and 3) a single band, dust emission flux on the Rayleigh-Jeans side of the SED. All techniques return consistent M_H2 estimates within a factor of ~2 or less, yielding gas depletion time-scales (tau_dep ~ 0.35 Gyrs) and gas-to-stellar mass ratios (M_H2/M* ~ 0.5-1) for our z~3 massive MS galaxies. The overall properties of our galaxies are consistent with trends and relations established at lower redshifts, extending the apparent uniformity of star-forming galaxies over the last 11.5 billion years.
Lyman Limit systems (LLSs) trace the low-density circumgalactic medium and the most dense regions of the intergalactic medium, so their number density and evolution at high redshift, just after reionisation, are important to constrain. We present a survey for LLSs at high redshifts, $z_{rm LLS} =3.5$--5.4, in the homogeneous dataset of 153 optical quasar spectra at $z sim 5$ from the Giant Gemini GMOS survey. Our analysis includes detailed investigation of survey biases using mock spectra which provide important corrections to the raw measurements. We estimate the incidence of LLSs per unit redshift at $z approx 4.4$ to be $ell(z) = 2.6 pm 0.4$. Combining our results with previous surveys at $z_{rm LLS} <4$, the best-fit power-law evolution is $ell(z) = ell_* [(1+z)/4]^alpha$ with $ell_* = 1.46 pm 0.11$ and $alpha = 1.70 pm 0.22$ (68% confidence intervals). Despite hints in previous $z_{rm LLS} <4$ results, there is no indication for a deviation from this single power-law soon after reionization. Finally, we integrate our new results with previous surveys of the intergalactic and circumgalactic media to constrain the hydrogen column density distribution function, $f(N_{rm HI},X)$, over 10 orders of magnitude. The data at $z sim 5$ are not well described by the $f(N_{rm HI},X)$ model previously reported for $z sim 2$--3 (after re-scaling) and a 7-pivot model fitting the full $z sim 2$--5 dataset is statistically unacceptable. We conclude that there is significant evolution in the shape of $f(N_{rm HI},X)$ over this $sim$2 billion year period.