We present new estimates of the luminosity function (LF) and star formation rate density (SFRD) for an Hα-selected sample at z ∼ 0.62 from the Deep And Wide Narrow-band (DAWN) survey. Our results are based on a new Hα sample in the extended COSMOS region (compared to Coughlin et al.) with the inclusion of flanking fields, resulting in a total area coverage of ∼1.5 deg2. A total of 241 Hα emitters were selected based on robust selection criteria using spectrophotometric redshifts and broadband color–color classification. Given that dust extinction is a dominant uncertainty in the estimation of LF and SFRD, we explore the effect of different dust-correction prescriptions by calculating the LF and SFRD using a constant dust extinction correction, AHα = 1 mag, a luminosity-dependent correction, and a stellar-mass-dependent correction. The resulting Hα LFs are well fitted using Schechter functions with best-fit parameters: L* = 1042.24 erg s−1, ϕ* = 10−2.85 Mpc−3, α = −1.62 for constant dust correction, L erg s−1, ϕ* = 10−2.8 Mpc−3, α = −1.39 for luminosity-dependent dust correction, and L* = 1042.36 erg s−1, ϕ* = 10−2.91 Mpc−3, α = −1.48, for stellar-mass-dependent dust correction. The deep and wide nature of the DAWN survey effectively samples Hα emitters over a wide range of luminosities, thereby providing better constraints on both the faint and bright ends of the LF. Also, the SFRD estimates ρSFR = 10−1.39 M⊙ yr−1 Mpc−3 (constant dust correction), ρSFR = 10−1.47 M⊙ yr−1 Mpc−3 (luminosity-dependent dust correction), and ρSFR = 10−1.46 M⊙ yr−1 Mpc−3 (stellar-mass-dependent dust correction) are in good agreement with the evolution of SFRD across redshifts (0 < z < 2) seen from previous Hα surveys.
ABSTRACT We investigate the relation between active galactic nucleus (AGN) and star formation (SF) activity at 0.5 < z < 3 by analysing 898 galaxies with X-ray luminous AGNs (LX > 1044 erg s−1) and a large comparison sample of ∼320 000 galaxies without X-ray luminous AGNs. Our samples are selected from a large (11.8 deg2) area in Stripe 82 that has multiwavelength (X-ray to far-IR) data. The enormous comoving volume (∼0.3 Gpc3) at 0.5 < z < 3 minimizes the effects of cosmic variance and captures a large number of massive galaxies (∼30 000 galaxies with M* > 1011 M⊙) and X-ray luminous AGNs. While many galaxy studies discard AGN hosts, we fit the SED of galaxies with and without X-ray luminous AGNs with Code Investigating GALaxy Emission and include AGN emission templates. We find that without this inclusion, stellar masses and star formation rates (SFRs) in AGN host galaxies can be overestimated, on average, by factors of up to ∼5 and ∼10, respectively. The average SFR of galaxies with X-ray luminous AGNs is higher by a factor of ∼3–10 compared to galaxies without X-ray luminous AGNs at fixed stellar mass and redshift, suggesting that high SFRs and high AGN X-ray luminosities may be fuelled by common mechanisms. The vast majority ($\gt 95 {{\ \rm per\ cent}}$) of galaxies with X-ray luminous AGNs at z = 0.5−3 do not show quenched SF: this suggests that if AGN feedback quenches SF, the associated quenching process takes a significant time to act and the quenched phase sets in after the highly luminous phases of AGN activity.
We present a joint analysis of the rest-frame ultraviolet (UV) luminosity functions of continuum-selected star-forming galaxies and galaxies dominated by active galactic nuclei (AGNs) at $z \sim$ 4. These 3,740 $z \sim$ 4 galaxies are selected from broad-band imaging in nine photometric bands over 18 deg$^2$ in the \textit{Spitzer}/HETDEX Exploratory Large Area Survey (SHELA) field. The large area and moderate depth of our survey provide a unique view of the intersection between the bright end of the galaxy UV luminosity function (M$_{AB}<-$22) and the faint end of the AGN UV luminosity function. We do not separate AGN-dominated galaxies from star-formation-dominated galaxies, but rather fit both luminosity functions simultaneously. These functions are best fit with a double power-law (DPL) for both the galaxy and AGN components, where the galaxy bright-end slope has a power-law index of $-3.80\pm0.10$, and the corresponding AGN faint-end slope is $\alpha_{AGN} = -1.49^{+0.30}_{-0.21}$. We cannot rule out a Schechter-like exponential decline for the galaxy UV luminosity function, and in this scenario the AGN luminosity function has a steeper faint-end slope of $-2.08^{+0.18}_{-0.11}$. Comparison of our galaxy luminosity function results with a representative cosmological model of galaxy formation suggests that the molecular gas depletion time must be shorter, implying that star formation is more efficient in bright galaxies at $z=4$ than at the present day. If the galaxy luminosity function does indeed have a power-law shape at the bright end, the implied ionizing emissivity from AGNs is not inconsistent with previous observations. However, if the underlying galaxy distribution is Schechter, it implies a significantly higher ionizing emissivity from AGNs at this epoch.
We present a sample of 30 massive (log$(M_{\ast}/M_\odot) >11$) $z=3-5$ quiescent galaxies selected from the \textit{Spitzer-}HETDEX Exploratory Large Area (SHELA) Survey and observed at 1.1mm with Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 observations. These ALMA observations would detect even modest levels of dust-obscured star-formation, on order of $\sim 20 \ M_\odot \textrm{yr}^{-1}$ at $z\sim4$ at a $1σ$ level, allowing us to quantify the amount of contamination from dusty star-forming sources in our quiescent sample. Starting with a parent sample of candidate massive quiescent galaxies from the Stevans et al. 2021 v1 SHELA catalog, we use the Bayesian \textsc{Bagpipes} spectral energy distribution fitting code to derive robust stellar masses ($M_*$) and star-formation rates (SFRs) for these sources, and select a conservative sample of 36 candidate massive ($M_* > 10^{11}M_\odot$) quiescent galaxies, with specific SFRs at $>2σ$ below the star-forming main sequence at $z\sim4$. Based on ALMA imaging, six of these candidate quiescent galaxies have the presence of significant dust-obscured star-formation, thus were removed from our final sample. This implies a $\sim 17\%$ contamination rate from dusty star-forming galaxies with our selection criteria using the v1 SHELA catalog. This conservatively-selected quiescent galaxy sample at $z=3-5$ will provide excellent targets for future observations to better constrain how massive galaxies can both grow and shut-down their star-formation in a relatively short time period.
Abstract The study of galaxy evolution hinges on our ability to interpret multiwavelength galaxy observations in terms of their physical properties. To do this, we rely on spectral energy distribution (SED) models, which allow us to infer physical parameters from spectrophotometric data. In recent years, thanks to wide and deep multiwave band galaxy surveys, the volume of high-quality data have significantly increased. Alongside the increased data, algorithms performing SED fitting have improved, including better modeling prescriptions, newer templates, and more extensive sampling in wavelength space. We present a comprehensive analysis of different SED-fitting codes including their methods and output with the aim of measuring the uncertainties caused by the modeling assumptions. We apply 14 of the most commonly used SED-fitting codes on samples from the CANDELS photometric catalogs at z ∼ 1 and z ∼ 3. We find agreement on the stellar mass, while we observe some discrepancies in the star formation rate (SFR) and dust-attenuation results. To explore the differences and biases among the codes, we explore the impact of the various modeling assumptions as they are set in the codes (e.g., star formation histories, nebular, dust and active galactic nucleus models) on the derived stellar masses, SFRs, and A V values. We then assess the difference among the codes on the SFR–stellar mass relation and we measure the contribution to the uncertainties by the modeling choices (i.e., the modeling uncertainties) in stellar mass (∼0.1 dex), SFR (∼0.3 dex), and dust attenuation (∼0.3 mag). Finally, we present some resources summarizing best practices in SED fitting.
The study of galaxy evolution hinges on our ability to interpret multi-wavelength galaxy observations in terms of their physical properties. To do this, we rely on spectral energy distribution (SED) models which allow us to infer physical parameters from spectrophotometric data. In recent years, thanks to the wide and deep multi-waveband galaxy surveys, the volume of high quality data have significantly increased. Alongside the increased data, algorithms performing SED fitting have improved, including better modeling prescriptions, newer templates, and more extensive sampling in wavelength space. We present a comprehensive analysis of different SED fitting codes including their methods and output with the aim of measuring the uncertainties caused by the modeling assumptions. We apply fourteen of the most commonly used SED fitting codes on samples from the CANDELS photometric catalogs at z~1 and z~3. We find agreement on the stellar mass, while we observe some discrepancies in the star formation rate (SFR) and dust attenuation results. To explore the differences and biases among the codes, we explore the impact of the various modeling assumptions as they are set in the codes (e.g., star formation histories, nebular, dust, and AGN models) on the derived stellar masses, SFRs, and A_V values. We then assess the difference among the codes on the SFR-stellar mass relation and we measure the contribution to the uncertainties by the modeling choices (i.e., the modeling uncertainties) in stellar mass (~0.1dex), SFR (~0.3dex), and dust attenuation (~0.3mag). Finally, we present some resources summarizing best practices in SED fitting.
We present an ultraviolet spectroscopic survey of strong H I absorbers in the intergalactic medium, probing their evolution over the last 6-7 Gyr at redshifts $0.24 \leq z \leq 0.84$. We measure column densities $N_{\rm HI} \,( {\rm cm}^{-2})$ from the pattern of Lyman-series absorption lines and flux decrement at the Lyman limit (LL) when available. We analyzed 220 H I absorbers in ultraviolet spectra of 102 active galactic nuclei (AGN) taken by the Cosmic Origins Spectrograph aboard the Hubble Space Telescope with G130M/G160M gratings (1134-1795 \AA). For 158 absorbers with $\log N_{\rm HI} \geq 15$, the mean frequency is $d {\cal N}/dz = 4.95 \pm 0.39$ over pathlength $\Delta z = 31.94$ ($0.24 \leq z \leq 0.84)$. We identify 8 Lyman Limit Systems (LLS, $\log N_{\rm HI} \geq 17.2$) and 54 partial systems (pLLS) with $16.0 \leq \log N_{\rm HI} < 17.2$. Toward 159 AGN between $0.01 < z_{\rm abs} < 0.84$ with $\Delta z \approx 48$, we find four damped Ly$\alpha$ absorbers (DLA) with $(d {\cal N}/dz)_{\rm DLA} = 0.083^{+0.066}_{-0.040}$ at $\langle z \rangle = 0.18$. The mean LLS frequency between $z = 0.24-0.48$ is $(d {\cal N}/dz)_{\rm LLS} = 0.36^{+0.20}_{-0.13}$ fitted to $N(z) = (0.25^{+0.13}_{-0.09})(1+z)^{1.14}$. For 54 pLLS we find $(d {\cal N}/dz)_{\rm pLLS} = 1.69\pm0.23$ at $\langle z \rangle = 0.39$, a frequency consistent with gaseous halo sizes $R \approx 100 h^{-1}~{\rm kpc}$ for ($0.3-3L^*$) galaxies. A maximum-likelihood analysis yields a distribution $f(N,z) = C_0 N^{-\beta} (1+z)^{\gamma}$ with $\beta = 1.48 \pm 0.05$ and $\gamma = 1.14^{+0.88}_{-0.89}$ for $15 \leq \log N_{\rm HI} \leq 17.5$. The far-UV opacity gradient is $d \tau_{\rm eff} / dz \approx (0.444)(1+z)^{1.14}$ over the range $15 \leq \log N_{\rm HI} \leq 17$, implying mean LyC optical depth $\tau_{\rm eff} \approx 0.3-0.5$ toward sources at $z = 1-2$.
