ABSTRACT
We explore the production and escape of ionizing photons in young galaxies by investigating the ultraviolet and optical emission-line properties of models of ionization-bounded and ...density-bounded H ii regions, active-galactic-nucleus (AGN) narrow-line regions, and radiative shocks computed all using the same physically consistent description of element abundances and depletion on to dust grains down to very low metallicities. We compare these models with a reference sample of metal-poor star-forming galaxies and Lyman-continuum (LyC) leakers at various redshifts, which allows the simultaneous exploration of more spectral diagnostics than typically available at once for individual subsamples. We confirm that current single- and binary-star population synthesis models do not produce hard-enough radiation to account for the high-ionization emission of the most metal-poor galaxies. Introducing either an AGN or radiative-shock component brings models into agreement with observations. A published model including X-ray binaries is an attractive alternative to reproduce the observed rise in He iiλ4686/H β ratio with decreasing oxygen abundance in metal-poor star-forming galaxies, but not the high observed He iiλ4686/Hβ ratios of galaxies with large EW(Hβ). A source of harder ionizing radiation appears to be required in these extreme objects, such as an AGN or radiative-shock component, perhaps linked to an initial-mass-function bias towards massive stars at low metallicity. This would also account for the surprisingly high O i/O iii ratios of confirmed LyC leakers relative to ionization-bounded models. We find no simple by-eye diagnostic of the nature of ionizing sources and the escape of LyC photon, which require proper simultaneous fits of several lines to be discriminated against.
We present a new approach to investigate the content and spatial distribution of dust in structurally unresolved star-forming galaxies from the observed dependence of integrated spectral properties ...on galaxy inclination. Motivated by the observation that different stellar populations reside in different spatial components of nearby star-forming galaxies, we develop an innovative combination of generic models of radiative transfer in dusty media with a prescription for the spectral evolution of galaxies, via the association of different geometric components of galaxies with stars in different age ranges. We start by showing that a wide range of radiative transfer models all predict a quasi-universal relation between slope of the attenuation curve at any wavelength, from the ultraviolet to the near-infrared, and V-band attenuation optical depth in the diffuse interstellar medium (ISM), at all galaxy inclinations. This relation predicts steeper (shallower) dust attenuation curves than both the Calzetti and Milky Way curves at small (large) attenuation optical depths, which implies that geometry and orientation effects have a stronger influence on the shape of the attenuation curve than changes in the optical properties of dust grains. We use our new, combined radiative transfer and spectral evolution model to interpret the observed dependence of the Hα/Hβ ratio and ugrizYJH attenuation curve on inclination in a sample of about 23 000 nearby star-forming galaxies, which we correct for systematic biases by developing a general method based on importance sampling. From the exploration of the model parameter space by means of a Bayesian Markov chain Monte Carlo technique, we measure the central face-on B-band optical depth of this sample to be τB⊥ 1.8 ± 0.2 (corresponding to an angle-averaged
). We also quantify the enhanced optical depth towards newly formed stars in their birth clouds, finding this to be significantly larger in galaxies with bulges than in disc-dominated galaxies, while τ
B⊥ is roughly similar in both cases. This can arise if, for example, galaxies with significant bulges have higher central star formation efficiencies than their disc-dominated counterparts at a fixed specific star formation rate, and dustier stellar birth clouds because of the higher metallicity. We find that over 80 per cent of the attenuation in galaxies in our sample is characteristic of that affecting thin-disc stars in radiative transfer models. The median unattenuated V-band luminosity ratio of thick-disc to thin-disc stars is 0.1-0.2, in good agreement with the results from spatially resolved studies of nearby edge-on disc galaxies. Finally, we show that neglecting the effect of geometry and orientation on attenuation can severely bias the interpretation of galaxy spectral energy distributions, as the impact on broad-band colours can reach up to 0.3-0.4 mag at optical wavelengths and 0.1 mag at near-infrared ones. This paper also contains an original application of Gaussian random processes to extend the wavelength range of dust attenuation curves.
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.
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 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.
We present an empirical method of assessing the star formation rate (SFR) of star-forming galaxies based on their locations in the rest-frame color–color diagram (NUV − r) vs. (r − K). By using the ...Spitzer 24 μm sample in the COSMOS field (~16 400 galaxies with 0.2 ≤ z ≤ 1.3) and a local GALEX-SDSS-SWIRE sample (~700 galaxies with z ≤ 0.2), we show that the mean infrared excess ⟨IRX⟩ = ⟨ LIR/LUV ⟩ can be described by a single vector, NRK , that combines the two colors. The calibration between ⟨IRX⟩ and NRK allows us to recover the IR luminosity, LIR, with an accuracy of σ ~ 0.21 for the COSMOS sample and 0.27 dex for the local one. The SFRs derived with this method agree with the ones based on the observed (UV+IR) luminosities and on the spectral energy distribution (SED) fitting for the vast majority (~85%) of the star-forming population. Thanks to a library of model galaxy SEDs with realistic prescriptions for the star formation history, we show that we need to include a two-component dust model (i.e., birth clouds and diffuse ISM) and a full distribution of galaxy inclinations in order to reproduce the behavior of the ⟨IRX⟩ stripes in the NUVrK diagram. In conclusion, the NRK method, based only on the rest-frame UV/optical colors available in most of the extragalactic fields, offers a simple alternative of assessing the SFR of star-forming galaxies in the absence of far-IR or spectral diagnostic observations.
Aims. 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 selected a sample of 52 galaxies that also have a precise systemic redshift determination from at least one nebular emission line (HeII or CIII). For these galaxies, we derived different physical properties (stellar mass, age, dust extinction, and star formation rate) from spectral energy distribution (SED) fitting of the exquisite multiwavelength photometry available in the VANDELS fields, using the dedicated spectral modeling tool BEAGLE and the UV β slope from the observed photometry. We characterized the Lyα emission in terms of kinematics, equivalent width (EW), full width at half-maximum, and spatial extension and then estimated the velocity of the neutral outflowing gas. The ultra-deep VANDELS spectra (up to 80 h on-source integration) enable this for individual galaxies without the need to rely on stacks. We then investigated 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 galaxies with lower mass that are less dusty and bluer. 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 in both UV and in Lyα. Our data reveal an interesting correlation between the Lyα velocity offset and the shift of the interstellar absorption lines with respect to the systemic redshift, observed for the first time at high redshifts: galaxies with higher interstellar medium (ISM) outflow 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; this is 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.
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.
We tested the performance of photometric redshifts for galaxies in the Hubble Ultra Deep field down to 30th magnitude. We compared photometric redshift estimates from three spectral fitting codes ...from the literature (EAZY, BPZ and BEAGLE) to high quality redshifts for 1227 galaxies from the MUSE integral field spectrograph. All these codes can return photometric redshifts with bias |(zMUSE−pz) / (1 + zMUSE)| < 0.05 down to F775W = 30 and spectroscopic incompleteness is unlikely to strongly modify this statement. We have, however, identified clear systematic biases in the determination of photometric redshifts: in the 0.4 < z < 1.5 range, photometric redshifts are systematically biased low by as much as (zMUSE−pz) / (1 + zMUSE) = −0.04 in the median, and at z> 3 they are systematically biased high by up to (zMUSE−pz) / (1 + zMUSE) = 0.05, an offset that can in part be explained by adjusting the amount of intergalactic absorption applied. In agreement with previous studies we find little difference in the performance of the different codes, but in contrast to those we find that adding extensive ground-based and IRAC photometry actually can worsen photo-z performance for faint galaxies. We find an outlier fraction, defined through |(zMUSE−pz) / (1 + zMUSE)| > 0.15, of 8% for BPZ and 10% for EAZY and BEAGLE, and show explicitly that this is a strong function of magnitude. While this outlier fraction is high relative to numbers presented in the literature for brighter galaxies, they are very comparable to literature results when the depth of the data is taken into account. Finally, we demonstrate that while a redshift might be of high confidence, the association of a spectrum to the photometric object can be very uncertain and lead to a contamination of a few percent in spectroscopic training samples that do not show up as catastrophic outliers, a problem that must be tackled in order to have sufficiently accurate photometric redshifts for future cosmological surveys.