ABSTRACT
We present a method to self-consistently propagate stellar-mass $\hbox{$\hbox{${\rm M}$}_{\star }$}=\log (\hbox{${\rm M}$}/\hbox{${\rm M}_{\odot }$})$ and star-formation-rate $\hbox{${\Psi ...}$}=\log (\hbox{${\psi }$}/\hbox{${\rm M}_{\odot }$}\, {\rm yr}^{-1}$) uncertainties on to intercept (α), slope (β), and intrinsic-scatter (σ) estimates for a simple model of the main sequence of star-forming galaxies, where $\hbox{${\Psi }$}= \alpha + \beta \hbox{$\hbox{${\rm M}$}_{\star }$}+ \mathcal {N}(0,\sigma)$. To test this method and compare it with other published methods, we construct mock photometric samples of galaxies at z ∼ 5 based on idealized models combined with broad- and medium-band filters at wavelengths 0.8–5 μm. Adopting simple Ψ estimates based on dust-corrected ultraviolet luminosity can underestimate σ. We find that broad-band fluxes alone cannot constrain the contribution from emission lines, implying that strong priors on the emission-line contribution are required if no medium-band constraints are available. Therefore, at high redshifts, where emission lines contribute a higher fraction of the broad-band flux, photometric fitting is sensitive to Ψ variations on short (∼10 Myr) time-scales. Priors on age imposed with a constant (or rising) star formation history (SFH) do not allow one to investigate a possible dependence of σ on $\hbox{${\rm M}$}_{\star }$ at high redshifts. Delayed exponential SFHs have less constrained priors, but do not account for Ψ variations on short time-scales, a problem if σ increases due to stochasticity of star formation. A simple SFH with current star formation decoupled from the previous history is appropriate. We show that, for simple exposure-time calculations assuming point sources, with low levels of dust, we should be able to obtain unbiased estimates of the main sequence down to $\mathrm{ log}(\hbox{${\rm M}$}/\hbox{${\rm M}_{\odot }$})\sim 8$ at z ∼ 5 with the James Webb Space Telescope while allowing for stochasticity of star formation.
We present a new determination of the ultraviolet (UV) galaxy luminosity function (LF) at redshift z 7 and 8, and a first estimate at z 9. An accurate determination of the form and evolution of the ...galaxy LF during this era is of key importance for improving our knowledge of the earliest phases of galaxy evolution and the process of cosmic reionization. Our analysis exploits to the full the new, deepest Wide Field Camera 3/infrared imaging from our Hubble Space Telescope (HST) Ultra-Deep Field 2012 (UDF12) campaign, with dynamic range provided by including a new and consistent analysis of all appropriate, shallower/wider area HST survey data. Our new measurement of the evolving LF at z 7 to 8 is based on a final catalogue of 600 galaxies, and involves a step-wise maximum-likelihood determination based on the photometric redshift probability distribution for each object; this approach makes full use of the 11-band imaging now available in the Hubble Ultra-Deep Field (HUDF), including the new UDF12 F140W data, and the latest Spitzer IRAC imaging. The final result is a determination of the z 7 LF extending down to UV absolute magnitudes M
1500 = −16.75 (AB mag) and the z 8 LF down to M
1500 = −17.00. Fitting a Schechter function, we find M1500
* = −19.90+0.23
−0.28, log φ* = −2.96+0.18
−0.23 and a faint-end slope α = −1.90+0.14
−0.15 at z 7, and M1500* = −20.12+0.37
−0.48, log φ* = −3.35+0.28
−0.47 and α = −2.02+0.22
+0.23 at z 8. These results strengthen previous suggestions that the evolution at z > 7 appears more akin to 'density evolution' than the apparent 'luminosity evolution' seen at z 5 − 7. We also provide the first meaningful information on the LF at z 9, explore alternative extrapolations to higher redshifts, and consider the implications for the early evolution of UV luminosity density. Finally, we provide catalogues (including derived z
phot, M
1500 and photometry) for the most robust z ∼ 6.5-11.9 galaxies used in this analysis. We briefly discuss our results in the context of earlier work and the results derived from an independent analysis of the UDF12 data based on colour-colour selection.
ABSTRACT
The Hubble Frontier Fields represent the opportunity to probe the high-redshift evolution of the main sequence of star-forming galaxies to lower masses than possible in blank fields thanks ...to foreground lensing of massive galaxy clusters. We use the beagle SED-fitting code to derive stellar masses, $\rm{{M_{\star }}}=\log ({\it M}/{\rm{M_{\odot }}})$, SFRs, $\rm{{\Psi }}=\log (\rm{{\psi }}/{\rm{M_{\odot }}}\, {\rm{yr}}^{-1})$, and redshifts from galaxies within the astrodeep catalogue. We fit a fully Bayesian hierarchical model of the main sequence over 1.25 < z < 6 of the form $\rm{{\Psi }}= \rm{\alpha _\mathrm{9.7}}(z) + \rm{\beta }({\rm{M_{\star }}}-9.7) + \mathcal {N}(0,\rm{\sigma }^2)$ while explicitly modelling the outlier distribution. The redshift-dependent intercept at $\rm{{M_{\star }}}=9.7$ is parametrized as $\rm{\alpha _\mathrm{9.7}}(z) = \log {\it N}(1+{\it z})^{\rm{\gamma }} + 0.7$. Our results agree with an increase in normalization of the main sequence to high redshifts that follows the redshift-dependent rate of accretion of gas on to dark matter haloes with $\rm{\gamma }=2.40^{+0.18}_{-0.18}$. We measure a slope and intrinsic scatter of $\rm{\beta }=0.79^{+0.03}_{-0.04}$ and $\rm{\sigma }=0.26^{+0.02}_{-0.02}$. We find that the sampling of the SED provided by the combination of filters (Hubble + ground-based Ks-band + Spitzer 3.6 and 4.5 μm) is insufficient to constrain M⋆ and Ψ over the full dynamic range of the observed main sequence, even at the lowest redshifts studied. While this filter set represents the best current sampling of high-redshift galaxy SEDs out to z > 3, measurements of the main sequence to low masses and high redshifts still strongly depend on priors employed in SED fitting (as well as other fitting assumptions). Future data sets with JWST should improve this.
We present spectroscopic confirmation of 10 highly luminous (L≥ 2L
★) Lyα emitters (LAEs) in the redshift range 6.01 < z < 6.49 (nine galaxies and one active galactic nucleus), initially drawn from a ...sample of 14 z
phot≥ 6 Lyman-break galaxies (LBGs) selected from an area of 0.25 deg2 within the UK Infrared Deep Sky Survey (UKIDSS) Ultra-Deep Survey (UDS). Overall, our high rate of spectroscopic confirmation (≥71 per cent) and low rate of contamination provides a strong vindication of the photometric redshift analysis used to define the original sample. By considering star formation rate estimates based on the Lyα and ultraviolet continuum luminosity we conclude that our sample is consistent with a Lyα escape fraction of ≃25 per cent. Moreover, after careful consideration of the potential uncertainties and biases, we find that 40-50 per cent of our sample of L≥ 2L
★ galaxies at 6.0 < z < 6.5 display strong Lyα emission (rest-frame equivalent width ≥25 Å), a fraction which is a factor of ≃2 higher than previously reported for L≤L
★ galaxies at z≃ 6. Our results suggest that, as the epoch of reionization is approached, it is plausible that the LAE fraction amongst luminous (L≥ 2L
★) LBGs shows a similarly sharp increase to that observed in their lower luminosity (L≤L
★) counterparts.
We use the new ultra-deep, near-infrared imaging of the Hubble Ultra-Deep Field (HUDF) provided by our UDF12 Hubble Space Telescope (HST) Wide Field Camera 3/IR campaign to explore the rest-frame ...ultraviolet (UV) properties of galaxies at redshifts z > 6.5. We present the first unbiased measurement of the average UV power-law index, 〈β〉, (f
λ ∝ λβ) for faint galaxies at z 7, the first meaningful measurements of 〈β〉 at z 8, and tentative estimates for a new sample of galaxies at z 9. Utilizing galaxy selection in the new F140W (J
140) imaging to minimize colour bias, and applying both colour and power-law estimators of β, we find 〈β〉 = −2.1 ± 0.2 at z 7 for galaxies with M
UV −18. This means that the faintest galaxies uncovered at this epoch have, on average, UV colours no more extreme than those displayed by the bluest star-forming galaxies at low redshift. At z 8 we find a similar value, 〈β〉 = −1.9 ± 0.3. At z 9, we find 〈β〉 = −1.8 ± 0.6, essentially unchanged from z 6 to 7 (albeit highly uncertain). Finally, we show that there is as yet no evidence for a significant intrinsic scatter in β within our new, robust z 7 galaxy sample. Our results are most easily explained by a population of steadily star-forming galaxies with either solar metallicity and zero dust, or moderately sub-solar ( 10-20 per cent) metallicity with modest dust obscuration (A
V
0.1-0.2). This latter interpretation is consistent with the predictions of a state-of-the-art galaxy-formation simulation, which also suggests that a significant population of very-low metallicity, dust-free galaxies with β −2.5 may not emerge until M
UV > −16, a regime likely to remain inaccessible until the James Webb Space Telescope.
We present a new prospective analysis of deep multi-band imaging with the
James Webb
Space Telescope (JWST). In this work, we investigate the recovery of high-redshift 5 <
z
< 12 galaxies through ...extensive image simulations of accepted JWST programs, including the Early Release Science in the EGS field and the Guaranteed Time Observations in the HUDF. We introduced complete samples of ∼300 000 galaxies with stellar masses of log(
M
*
/
M
⊙
) > 6 and redshifts of 0 <
z
< 15, as well as galactic stars, into realistic mock NIRCam, MIRI, and HST images to properly describe the impact of source blending. We extracted the photometry of the detected sources, as in real images, and estimated the physical properties of galaxies through spectral energy distribution fitting. We find that the photometric redshifts are primarily limited by the availability of blue-band and near-infrared medium-band imaging. The stellar masses and star formation rates are recovered within 0.25 and 0.3 dex, respectively, for galaxies with accurate photometric redshifts. Brown dwarfs contaminating the
z
> 5 galaxy samples can be reduced to < 0.01 arcmin
−2
with a limited impact on galaxy completeness. We investigate multiple high-redshift galaxy selection techniques and find that the best compromise between completeness and purity at 5 <
z
< 10 using the full redshift posterior probability distributions. In the EGS field, the galaxy completeness remains higher than 50% at magnitudes
m
UV
< 27.5 and at all redshifts, and the purity is maintained above 80 and 60% at
z
≤ 7 and 10, respectively. The faint-end slope of the galaxy UV luminosity function is recovered with a precision of 0.1–0.25, and the cosmic star formation rate density within 0.1 dex. We argue in favor of additional observing programs covering larger areas to better constrain the bright end.
ABSTRACT
We study the ionizing photon production efficiency at the end of the Epoch of Reionization (z ∼ 5.4 − 6.6) for a sample of 30 Ly α emitters. This is a crucial quantity to infer the ionizing ...photon budget of the universe. These objects were selected to have reliable spectroscopic redshifts, assigned based on the profile of their Ly α emission line, detected in the MUSE deep fields. We exploit medium-band observations from the JWST Extragalactic Medium-band Survey (JEMS) to find the flux excess corresponding to the redshifted Hα emission line. We estimate the ultraviolet (UV) luminosity by fitting the full JEMS photometry, along with several HST photometric points, with Prospector. We find a median UV continuum slope of $\beta = -2.09^{+0.23}_{-0.21}$, indicating young stellar populations with little-to-no dust attenuation. Supported by this, we derive ξion,0 with no dust attenuation and find a median value of log$\frac{\xi _{ion,0}}{\text{Hz erg}^{-1}} = 25.44^{+0.21}_{-0.15}$. If we perform dust attenuation corrections and assume a Calzetti attenuation law, our values are lowered by ∼0.1 dex. Our results suggest Ly α emitters at the Epoch of Reionization have slightly enhanced ξion,0 compared to previous estimations from literature, in particular, when compared to the non-Ly α emitting population. This initial study provides a promising outlook on the characterization of ionizing photon production in the early universe. In the future, a more extensive study will be performed on the entire data set provided by the JWST Advanced Deep Extragalactic Survey (JADES). Thus, for the first time, allowing us to place constraints on the wider galaxy populations driving reionization.
We present the results of a study investigating the sizes and morphologies of redshift 4 < z < 8 galaxies in the CANDELS (Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey) GOODS-S ...(Great Observatories Origins Deep Survey southern field), HUDF (Hubble Ultra-Deep Field) and HUDF parallel fields. Based on non-parametric measurements and incorporating a careful treatment of measurement biases, we quantify the typical size of galaxies at each redshift as the peak of the lognormal size distribution, rather than the arithmetic mean size. Parametrizing the evolution of galaxy half-light radius as r
50 ∝ (1 + z)
n
, we find n = −0.20 ± 0.26 at bright UV-luminosities (0.3L
*(z = 3) < L < L
*) and n = −0.47 ± 0.62 at faint luminosities (0.12L
* < L < 0.3L
*). Furthermore, simulations based on artificially redshifting our z ∼ 4 galaxy sample show that we cannot reject the null hypothesis of no size evolution. We show that this result is caused by a combination of the size-dependent completeness of high-redshift galaxy samples and the underestimation of the sizes of the largest galaxies at a given epoch. To explore the evolution of galaxy morphology we first compare asymmetry measurements to those from a large sample of simulated single Sérsic profiles, in order to robustly categorize galaxies as either ‘smooth’ or ‘disturbed’. Comparing the disturbed fraction amongst bright (M
1500 ≤ −20) galaxies at each redshift to that obtained by artificially redshifting our z ∼ 4 galaxy sample, while carefully matching the size and UV-luminosity distributions, we find no clear evidence for evolution in galaxy morphology over the redshift interval 4 < z < 8. Therefore, based on our results, a bright (M
1500 ≤ −20) galaxy at z ∼ 6 is no more likely to be measured as ‘disturbed’ than a comparable galaxy at z ∼ 4, given the current observational constraints.
We report the results of a comprehensive study of the relationship between galaxy size, stellar mass and specific star formation rate (sSFR) at redshifts 1.3 < z < 1.5. Based on a mass-complete (M
≥ ...6 × 1010 M), spectroscopic sample from the UK Infrared Deep Sky Survey (UKIDSS) Ultradeep Survey, with accurate stellar-mass measurements derived from spectro-photometric fitting, we find that at z 1.4 the location of massive galaxies on the size-mass plane is determined primarily by their sSFR. At this epoch, we find that massive galaxies which are passive (sSFR ≤ 0.1 Gyr−1) follow a tight size-mass relation, with half-light radii a factor of f
g = 2.4 ± 0.2 smaller than their local counterparts. Moreover, amongst the passive sub-sample we find no evidence that the off-set from the local size-mass relation is a function of stellar population age. In contrast, we find that massive star-forming galaxies at this epoch lie closer to the local late-type size-mass relation and are only a factor of f
g = 1.6 ± 0.2 smaller than observed locally. Based on a sub-sample with dynamical-mass estimates, which consists of both passive and star-forming objects, we also derive an independent estimate of f
g = 2.3 ± 0.3 for the typical growth in half-light radius between z 1.4 and the present day. Focusing on the passive sub-sample, we conclude that to produce the necessary evolution predominantly via major mergers would require an unfeasible number of merger events and overpopulate the high-mass end of the local stellar-mass function. In contrast, we find that a scenario in which mass accretion is dominated by minor mergers can comfortably produce the necessary evolution, whereby an increase in stellar mass of only a factor of 2, accompanied by size growth of a factor of 3.5, is required to reconcile the size-mass relation at z 1.4 with that observed locally. Finally, we note that a significant fraction (44 ± 12 per cent) of the passive galaxies in our sample have a disc-like morphology, providing additional evidence that separate physical processes are responsible for the quenching of star formation and morphological transformation in massive galaxies.
ABSTRACT
We present the addition of nebular emission from the narrow-line regions (NLR) surrounding active galactic nuclei (AGNs) to beagle (BayEsian Analysis of GaLaxy sEds). Using a set of ...idealized spectra, we fit to a set of observables (emission-line ratios and fluxes) and test the retrieval of different physical parameters. We find that fitting to standard diagnostic-line ratios plus O ii λ3726, λ3729/O iii λ5007, H β/H α, O i λ6300/O ii λ3726, λ3729, and H α flux, degeneracies remain between dust-to-metal mass ratio ($\xi _\rm {d}^{\small NLR}$) and ionization parameter ($U_\rm {s}^{\small NLR}$) in the NLR gas, and between slope of the ionizing radiation ($\alpha _{\small PL}$, characterizing the emission from the accretion disc around the central black hole) and total accretion-disc luminosity ($L_\rm {acc}$). Since these degeneracies bias the retrieval of other parameters even at maximal signal-to-noise ratio (S/N), without additional observables, we suggest fixing $\alpha _{\small PL}$ and dust-to-metal mass ratios in both NLR and H ii regions. We explore the S/N in H β required for un-biased estimates of physical parameters, finding that S/N(H β) ∼ 10 is sufficient to identify an NLR contribution, but that higher S/N is required for un-biased parameter retrieval (∼20 for NLR-dominated systems, ∼30 for objects with approximately equal H β contributions from NLR and H ii regions). We also compare the predictions of our models for different line ratios to previously published models and data. By adding He ii λ4686-line measurements to a set of published line fluxes for a sample of 463 AGN NLR, we show that our models with $-4\lt \hbox{$\log U_{\small S}^{\small NLR}$}\lt -1.5$ can account for the full range of observed AGN properties in the local Universe.