SHARDS: Constraints on the dust attenuation law of star-forming galaxies at z~2

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 $\Delta\beta$~0.4.Comment: 13 pages, 11+2 figures, 3 tables. MNRAS, in pres

Differences and similarities of stellar populations in LAEs and LBGs at z~3.4-6.8

Lyman alpha emitters (LAEs) and Lyman break galaxies (LBGs) represent the most common groups of star-forming galaxies at high z, and the differences between their inherent stellar populations (SPs) are a key factor in understanding early galaxy formation and evolution. We have run a set of SP burst-like models for a sample of 1558 sources at 3.4 < z < 6.8 from the Survey for High-z Absorption Red and Dead Sources (SHARDS) over the GOODS-N field. This work focuses on the differences between the three different observational subfamilies of our sample: LAE–LBGs, no-Ly α LBGs, and pure LAEs. Single and double SP synthetic spectra were used to model the spectral energy distributions, adopting a Bayesian information criterion to analyze under which situations a second SP is required. We find that the sources are well modelled using a single SP in ∼79 per cent of the cases. The best models suggest that pure LAEs are typically young low-mass galaxies (⁠t∼26+41−25 Myr; Mstar∼5.6+12.0−5.5×108 M⊙⁠), undergoing one of their first bursts of star formation. On the other hand, no-Ly α LBGs require older SPs (t ∼ 71 ± 12 Myr), and they are substantially more massive (Mstar ∼ 3.5 ± 1.1 × 109 M⊙). LAE–LBGs appear as the subgroup that more frequently needs the addition of a second SP, representing an old and massive galaxy caught in a strong recent star-forming episode. The relative number of sources found from each subfamily at each z supports an evolutionary scenario from pure LAEs and single SP LAE–LBGs to more massive LBGs. Stellar mass functions are also derived, finding an increase of M* with cosmic time and a possible steepening of the low-mass slope from z ∼ 6 to z ∼ 5 with no significant change to z ∼ 4. Additionally, we have derived the SFR–Mstar relation, finding an SFR∝Mβstar behaviour with negligible evolution from z ∼ 4 to z ∼ 6

Unveiling the hidden universe with JWST: The contribution of dust-obscured galaxies to the stellar mass function at z∼3−8\mathbf{z\sim3-8}

The emergence of massive, optically-faint galaxies in infrared observations has revealed that our view of the high-redshift Universe was previously incomplete. With the advent of JWST, we can for the first time probe the rest-frame optical emission of galaxies at $z>3$ with high sensitivity and spatial resolution, thus moving towards a more complete census of the galaxy population at high redshifts. To this end, we present a sample of 148 massive, dusty galaxies from the JWST/CEERS survey, colour-selected using solely JWST bands. With deep JWST/NIRCam data from 1.15$\mu$m to 4.44$\mu$m and ancillary HST/ACS and WFC3 data, we determine the physical properties of our sample using spectral energy distribution fitting with BAGPIPES. We demonstrate that our selection method efficiently identifies massive ($\mathrm{\langle \log M_\star/M_\odot \rangle \sim 10}$) and dusty ($\mathrm{\langle A_V\rangle \sim 2.7\ mag}$) sources, with a majority at $z>3$ and predominantly lying on the galaxy main-sequence. The main results of this work are the stellar mass functions (SMF) of red, optically-faint galaxies from redshifts between $3<z<8$: these galaxies make up a significant fraction of the pre-JWST total SMF at $3<z<4$, and dominate the high-mass end of the pre-JWST SMF at $4<z<6$ and $6<z<8$, suggesting that our census of the galaxy population needs amendment at these epochs. While larger areas need to be surveyed in the future, our results suggest already that the integrated stellar mass density at $\mathrm{\log M_\star/M_\odot>9.25}$ may have been underestimated by $\sim$20-25% at $z\sim3-6$, and $\sim$110% at $z\sim6-8$.Comment: 19 pages, 10 figures, submitted to MNRA

A simultaneous search for High-z LAEs and LBGs in the SHARDS survey

We have undertaken a comprehensive search for both Lyman Alpha Emitters (LAEs) and Lyman Break Galaxies (LBGs) in the SHARDS Survey of the GOODS-N field. SHARDS is a deep imaging survey, made with the 10.4 m Gran Telescopio Canarias (GTC), employing 25 medium band filters in the range from 500 to 941 nm. This is the first time that both LAEs and LBGs are surveyed simultaneously in a systematic way in a large field. We draw a sample of 1558 sources; 528 of them are LAEs. Most of the sources (1434) show rest-frame UV continua. A minority of them (124) are pure LAEs with virtually no continuum detected in SHARDS. We study these sources from z ∼ 3.35 up to z ∼ 6.8, well into the epoch of reionization. Note that surveys done with just one or two narrow band filters lack the possibility to spot the rest-frame UV continuum present in most of our LAEs. We derive redshifts, Star Formation Rates (SFRs), Lyα Equivalent Widths (EWs) and Luminosity Functions (LFs). Grouping within our sample is also studied, finding 92 pairs or groups of galaxies at the same redshift separated by less than 60 comoving kpc. In addition, we relate 87 and 55 UV-selected objects with two known overdensities at z = 4.05 and z = 5.198, respectively. Finally, we show that surveys made with broad band filters are prone to introduce many unwanted sources (∼20% interlopers), which means that previous studies may be overestimating the calculated LFs, specially at the faint end

From naked spheroids to disky galaxies: how do massive disk galaxies shape their morphology?

We investigate the assembly history of massive disk galaxies and describe how they shape their morphology through cosmic time. Using SHARDS and HST data, we modeled the surface brightness distribution of 91 massive galaxies at redshift $0.14<z\leq 1$ in the wavelength range $0.5-1.6$ $\mu$m, deriving the uncontaminated spectral energy distributions of their bulges and disks separately. This spectrophotometric decomposition allows us to compare the stellar populations properties of each component in individual galaxies. We find that the majority of massive galaxies ($\sim85\%$) builds inside-out, growing their extended stellar disk around the central spheroid. Some bulges and disks could start forming at similar epochs, but these bulges grow more rapidly than their disks, assembling $80\%$ of their mass in $\sim0.7$ Gyr and $\sim3.5$ Gyr, respectively. Moreover, we infer that both older bulges and older disks are more massive and compact than younger stellar structures. In particular, we find that bulges display a bimodal distribution of mass-weighted ages, i.e., they form in two waves. In contrast, our analysis of the disk components indicates that they form at $z\sim1$ for both first and second-wave bulges. This translates to first-wave bulges taking longer in acquiring a stellar disk ($5.2$ Gyr) compared to second-wave less-compact spheroids ($0.7$ Gyr). We do not find distinct properties (e.g., mass, star formation timescale, and mass surface density) for the disks in both types of galaxies. We conclude that the bulge mass and compactness mainly regulate the timing of the stellar disk growth, driving the morphological evolution of massive disk galaxies.Comment: Accepted for publication in ApJ (18 pages, 11 figures

Star-forming galaxies at low-redshift in the SHARDS survey

International audienceContext. The physical processes driving the evolution of star formation (SF) in galaxies over cosmic time still present many open questions. Recent galaxy surveys allow now to study these processes in great detail at intermediate redshift (0 ≤ z ≤ 0.5). Aims. We build a complete sample of star-forming galaxies and analyze their properties, reaching systems with low stellar masses and low star formation rates (SFRs) at intermediate-to-low redshift. Methods. We use data from the SHARDS multiband survey in the GOODS-North field. Its depth (up to magnitude m 3σ ∼ 26.5) and its spectro-photometric resolution (R ∼ 50) provides us with an ideal dataset to search for emission line galaxies (ELGs). We develop a new algorithm to identify low-redshift (z < 0.36) ELGs by detecting the [OIII]5007 and Hα emission lines simultaneously. We fit the spectral energy distribution (SED) of the selected sample, using a model with two single stellar populations. Results. We find 160 star-forming galaxies for which we derive equivalent widths (EWs) and absolute fluxes of both emission lines. We detect EWs as low as 12 Å, with median values for the sample of ∼35 Å in [OIII]5007 and ∼56 Å in Hα, respectively. Results from the SED fitting show a young stellar population with low median metallicity (36% of the solar value) and extinction (A V ∼ 0.37), with median galaxy stellar mass ∼10 8.5 M. Gas-phase metallicities measured from available spectra are also low. ELGs in our sample present bluer colours in the UVJ plane than the median colour-selected star-forming galaxy in SHARDS. We suggest a new V-J colour criterion to separate ELGs from non-ELGs in blue galaxy samples. In addition, several galaxies present high densities of O-type stars, possibly producing galactic superwinds, which makes them interesting targets for follow-up spectroscopy. Conclusions. We have demonstrated the efficiency of SHARDS in detecting low-mass ELGs (∼2 magnitudes deeper than previous spectroscopic surveys in the same field). The selected sample accounts for 20% of the global galaxy population at this redshift and luminosity, and is characterized by young SF bursts with sub-solar metallicities and low extinction. However, robust fits to the full SEDs can only be obtained including an old stellar population, suggesting the young component is built up by a recent burst of SF in an otherwise old galaxy

A simultaneous search for High-z LAEs and LBGs in the SHARDS survey

We have undertaken a comprehensive search for both Lyman Alpha Emitters (LAEs) and Lyman Break Galaxies (LBGs) in the SHARDS Survey of the GOODS-N field. SHARDS is a deep imaging survey, made with the 10.4 m Gran Telescopio Canarias (GTC), employing 25 medium band filters in the range from 500 to 941 nm. This is the first time that both LAEs and LBGs are surveyed simultaneously in a systematic way in a large field. We draw a sample of 1558 sources; 528 of them are LAEs. Most of the sources (1434) show rest-frame UV continua. A minority of them (124) are pure LAEs with virtually no continuum detected in SHARDS. We study these sources from z ∼ 3.35 up to z ∼ 6.8, well into the epoch of reionization. Note that surveys done with just one or two narrow band filters lack the possibility to spot the rest-frame UV continuum present in most of our LAEs. We derive redshifts, Star Formation Rates (SFRs), Lyα Equivalent Widths (EWs) and Luminosity Functions (LFs). Grouping within our sample is also studied, finding 92 pairs or groups of galaxies at the same redshift separated by less than 60 comoving kpc. In addition, we relate 87 and 55 UV-selected objects with two known overdensities at z = 4.05 and z = 5.198, respectively. Finally, we show that surveys made with broad band filters are prone to introduce many unwanted sources (∼20% interlopers), which means that previous studies may be overestimating the calculated LFs, specially at the faint end

Differential attenuation in star-forming galaxies at 0.3 ≲ z ≲ 1.5 in the SHARDS/CANDELS field

We use a sample of 706 galaxies, selected as [O II]λ3727 ([O II]) emitters in the Survey for High-z Absorption Red and Dead Sources (SHARDS) on the CANDELS/GOODS-N field, to study the differential attenuation of the nebular emission with respect to the stellar continuum. The sample includes only galaxies with a counterpart in the infrared and log10(M*/M⊙) > 9, over the redshift interval 0.3 ≲ z ≲ 1.5. Our methodology consists in the comparison of the star formation rates inferred from [O II] and Hα emission lines with a robust quantification of the total star-forming activity (SFR TOT) that is independently estimated based on both infrared and ultraviolet (UV) luminosities. We obtain $f\, =\, E(B-V)_{\mathrm{stellar}}$/E(B - V)nebular = 0.69$^{0.71}_{0.69}$ and 0.55$^{0.56}_{0.53}$ for [O II] and Hα, respectively. Our resulting f-factors display a significant positive correlation with the UV attenuation and shallower or not-significant trends with the stellar mass, the SFRTOT, the distance to the main sequence, and the redshift. Finally, our results favour an average nebular attenuation curve similar in shape to the typical dust curve of local starbursts

Differential attenuation in star-forming galaxies at 0.3 ≲ z ≲ 1.5 in the SHARDS/CANDELS field

We use a sample of 706 galaxies, selected as [O II]λ3727 ([O II]) emitters in the Survey for High-z Absorption Red and Dead Sources (SHARDS) on the CANDELS/GOODS-N field, to study the differential attenuation of the nebular emission with respect to the stellar continuum. The sample includes only galaxies with a counterpart in the infrared and log10(M∗/M) > 9, over the redshift interval 0.3 z 1.5. Our methodology consists in the comparison of the star formation rates inferred from [O II] and Hα emission lines with a robust quantification of the total star-forming activity (SFR TOT) that is independently estimated based on both infrared and ultraviolet (UV) luminosities. We obtain f = E(B − V )stellar/E(B − V)nebular = 0.690.71 0.69 and 0.550.56 0.53 for [O II] and Hα, respectively. Our resulting f-factors display a significant positive correlation with the UV attenuation and shallower or not-significant trends with the stellar mass, the SFRTOT, the distance to the main sequence, and the redshift. Finally, our results favour an average nebular attenuation curve similar in shape to the typical dust curve of local starbursts