We present a new-generation tool to model and interpret spectral energy distributions (SEDs) of galaxies, which incorporates in a consistent way the production of radiation and its transfer through ...the interstellar and intergalactic media. This flexible tool, named beagle (for BayEsian Analysis of GaLaxy sEds), allows one to build mock galaxy catalogues as well as to interpret any combination of photometric and spectroscopic galaxy observations in terms of physical parameters. The current version of the tool includes versatile modelling of the emission from stars and photoionized gas, attenuation by dust and accounting for different instrumental effects, such as spectroscopic flux calibration and line spread function. We show a first application of the beagle tool to the interpretation of broad-band SEDs of a published sample of
$\mathord {\sim }10^4$
galaxies at redshifts 0.1 ≲ z ≲ 8. We find that the constraints derived on photometric redshifts using this multipurpose tool are comparable to those obtained using public, dedicated photometric-redshift codes and quantify this result in a rigorous statistical way. We also show how the post-processing of beagle output data with the python extension pyp-beagle allows the characterization of systematic deviations between models and observations, in particular through posterior predictive checks. The modular design of the beagle tool allows easy extensions to incorporate, for example, the absorption by neutral galactic and circumgalactic gas, and the emission from an active galactic nucleus, dust and shock-ionized gas. Information about public releases of the beagle tool will be maintained on http://www.jacopochevallard.org/beagle.
We present a new model of the nebular emission from star-forming galaxies in a wide range of chemical compositions, appropriate to interpret observations of galaxies at all cosmic epochs. The model ...relies on the combination of state-of-the-art stellar population synthesis and photoionization codes to describe the ensemble of H ii regions and the diffuse gas ionized by young stars in a galaxy. A main feature of this model is the self-consistent yet versatile treatment of element abundances and depletion on to dust grains, which allows one to relate the observed nebular emission from a galaxy to both gas-phase and dust-phase metal enrichment. We show that this model can account for the rest-frame ultraviolet and optical emission-line properties of galaxies at different redshifts and find that ultraviolet emission lines are more sensitive than optical ones to parameters such as C/O abundance ratio, hydrogen gas density, dust-to-metal mass ratio and upper cut-off of the stellar initial mass function. We also find that, for gas-phase metallicities around solar to slightly subsolar, widely used formulae to constrain oxygen ionic fractions and the C/O ratio from ultraviolet and optical emission-line luminosities are reasonable faithful. However, the recipes break down at non-solar metallicities, making them inappropriate to study chemically young galaxies. In such cases, a fully self-consistent model of the kind presented in this paper is required to interpret the observed nebular emission.
We examine variations of the stellar initial mass function (IMF) in extreme environments within the formalism derived by Hennebelle & Chabrier. We focus on conditions encountered in progenitors of ...massive early-type galaxies and starburst regions. We show that, when applying the concept of turbulent Jeans mass as the characteristic mass for fragmentation in a turbulent medium, the peak of the IMF in such environments is shifted toward smaller masses, leading to a bottom-heavy IMF, as suggested by various observations. In very dense and turbulent environments, we predict that the high-mass tail of the IMF can become even steeper than the standard Salpeter IMF, with a limit for the power-law exponent alpha Asymptotically = to -2.7, in agreement with recent observational determinations. This steepening is a direct consequence of the high densities and Mach values in such regions but also of the time dependence of the fragmentation process, as incorporated in the Hennebelle-Chabrier theory. We provide analytical parameterizations of these IMFs in such environments tobe used in galaxy evolution calculations. We also calculate the star-formation rates and the mass-to-light ratios expected under such extreme conditions and show that they agree well with the values inferred in starburst environments and massive high-redshift galaxies. This reinforces the paradigm of star formation as being a universal process, i.e., the direct outcome of gravitationally unstable fluctuations in a density field initially generated by large-scale, shock-dominated turbulence. This globally enables us to infer the variations of the stellar IMF and related properties for atypical galactic conditions.
ABSTRACT
Galaxies in the reionization era have been shown to have prominent O iii + H β emission. Little is known about the gas conditions and radiation field of this population, making it ...challenging to interpret the spectra emerging at z ≳ 6. Motivated by this shortcoming, we have initiated a large MMT spectroscopic survey identifying rest-frame optical emission lines in 227 intense O iii emitting galaxies at 1.3 < z < 2.4. This sample complements the MOSDEF and KBSS surveys, extending to much lower stellar masses ($10^7\!-\!10^8 \, \mathrm{M}_\odot$) and larger specific star formation rates (5–300 Gyr−1), providing a window on galaxies directly following a burst or recent upturn in star formation. The hydrogen ionizing production efficiency (ξion) is found to increase with the O iii equivalent width (EW), in a manner similar to that found in local galaxies. We describe how this relationship helps explain the anomalous success rate in identifying Ly α emission in z ≳ 7 galaxies with strong O iii + H β emission. We probe the impact of the intense radiation field on the ISM using O32 and Ne3O2, two ionization-sensitive indices. Both are found to scale with the O iii EW, revealing extreme ionization conditions not commonly seen in older and more massive galaxies. In the most intense line emitters, the indices have very large average values (O32 = 9.1, Ne3O2 = 0.5) that have been shown to be linked to ionizing photon escape. We discuss implications for the nature of galaxies most likely to have O32 values associated with significant LyC escape. Finally we consider the optimal strategy for JWST spectroscopic investigations of galaxies at z ≳ 10 where the strongest rest-frame optical lines are no longer visible with NIRSpec.
We present a simple, largely empirical but physically motivated model to interpret the mid- and far-infrared spectral energy distributions of galaxies consistently with the emission at ultraviolet, ...optical and near-infrared wavelengths. Our model relies on an existing angle-averaged prescription to compute the absorption of starlight by dust in stellar birth clouds and in the ambient interstellar medium (ISM) in galaxies. We compute the spectral energy distribution of the power reradiated by dust in stellar birth clouds as the sum of three components: a component of polycyclic aromatic hydrocarbons (PAHs); a mid-infrared continuum characterizing the emission from hot grains at temperatures in the range 130–250 K; and a component of grains in thermal equilibrium with adjustable temperature in the range 30–60 K. In the ambient ISM, we fix for simplicity the relative proportions of these three components to reproduce the spectral shape of diffuse cirrus emission in the Milky Way, and we include a component of cold grains in thermal equilibrium with adjustable temperature in the range 15–25 K. Our model is both simple and versatile enough that it can be used to derive statistical constraints on the star formation histories and dust contents of large samples of galaxies using a wide range of ultraviolet, optical and infrared observations. We illustrate this by deriving median likelihood estimates of the star formation rates, stellar masses, effective dust optical depths, dust masses and relative strengths of different dust components of 66 well-studied nearby star-forming galaxies from the Spitzer Infrared Nearby Galaxy Survey (SINGS). We explore how the constraints derived in this way depend on the available spectral information. From our analysis of the SINGS sample, we conclude that the mid- and far-infrared colours of galaxies correlate strongly with the specific star formation rate, as well as with other galaxy-wide quantities connected to this parameter, such as the ratio of infrared luminosity between stellar birth clouds and the ambient ISM, the contributions by PAHs and grains in thermal equilibrium to the total infrared emission, and the ratio of dust mass to stellar mass. Our model can be straightforwardly applied to interpret ultraviolet, optical and infrared spectral energy distributions from any galaxy sample.
ABSTRACT In this paper, we present the GALEX-SDSS-WISE Legacy Catalog (GSWLC), a catalog of physical properties (stellar masses, dust attenuations, and star formation rates SFRs) for ∼700,000 ...galaxies with Sloan Digital Sky Survey (SDSS) redshifts below 0.3. GSWLC contains galaxies within the Galaxy Evolution Explorer footprint, regardless of a UV detection, covering 90% of SDSS. The physical properties were obtained from UV/optical spectral energy distribution (SED) fitting following Bayesian methodology of Salim et al., with improvements such as blending corrections for low-resolution UV photometry, flexible dust attenuation laws, and emission-line corrections. GSWLC also includes mid-IR SFRs derived from IR templates based on 22 Wide-field Infrared Survey Explorer observations. These estimates are independent of UV/optical SED fitting, in order to separate possible systematics. The paper argues that the comparison of specific SFRs (sSFRs) is more informative and physically motivated than the comparison of SFRs. The sSFRs resulting from the UV/optical SED fitting are compared to the mid-IR sSFRs and to sSFRs from three published catalogs. For "main-sequence" galaxies with no active galactic nucleus (AGN) all sSFRs are in very good agreement (within 0.1 dex on average). In particular, the widely used aperture-corrected SFRs from the MPA/JHU catalog show no systematic offsets, in contrast to some integral field spectroscopy results. For galaxies below the main sequence (log sSFR ), mid-IR (s)SFRs based on fixed luminosity-SFR conversion are severely biased (up to 2 dex) because the dust is primarily heated by old stars. Furthermore, mid-IR (s)SFRs are overestimated by up to 0.6 dex for galaxies with AGNs, presumably due to nonstellar dust heating. UV/optical (s)SFRs are thus preferred to IR-based (s)SFRs for quenched galaxies and those that host AGNs.
ABSTRACT
The variety of star formation histories (SFHs) of z ≳ 6 galaxies provides important insights into early star formation, but has been difficult to systematically quantify. Some observations ...suggest that many z ∼ 6–9 galaxies are dominated by ≳200 Myr stellar populations, implying significant star formation at z ≳ 9, while others find that most reionization era galaxies are ≲10 Myr, consistent with little z ≳ 9 star formation. Here, we quantify the distribution of ages of UV-bright ($-22.5\lesssim M_{\rm \small UV}\lesssim -21$) galaxies colour-selected to lie at z ≃ 6.6–6.9, an ideal redshift range to systematically study the SFHs of reionization era galaxies with ground-based observatories and Spitzer. We infer galaxy properties with two SED modelling codes and compare results, finding that stellar masses are largely insensitive to the model, but the inferred ages can vary by an order of magnitude. We infer a distribution of ages assuming a simple, parametric SFH model, finding a median age of ∼30–70 Myr depending on SED model. We quantify the fractions of ≤10 and ≥250 Myr galaxies, finding that these systems comprise ∼15–30 per cent and ∼20–25 per cent of the population, respectively. With a flexible SFH model, the shapes of the SFHs are consistent with those implied by the simple model (e.g. young galaxies have rapidly rising SFHs). However, stellar masses can differ significantly, with those of young systems sometimes being more than an order of magnitude larger with the flexible SFH. We quantify the implications of these results for z ≳ 9 stellar mass assembly and discuss improvements expected from JWST.
Abstract
Nearby dwarf galaxies provide a unique laboratory in which to test stellar population models below Z⊙/2. Such tests are particularly important for interpreting the surprising high-ionization ...ultraviolet (UV) line emission detected at z > 6 in recent years. We present HST/COS UV spectra of 10 nearby metal-poor star-forming galaxies selected to show He ii emission in SDSS optical spectra. The targets span nearly a dex in gas-phase oxygen abundance (7.8 < 12 + log O/H < 8.5) and present uniformly large specific star formation rates (sSFR ∼102 Gyr−1). The UV spectra confirm that metal-poor stellar populations can power extreme nebular emission in high-ionization UV lines, reaching C iii equivalent widths comparable to those seen in systems at z ∼ 6–7. Our data reveal a marked transition in UV spectral properties with decreasing metallicity, with systems below 12 + log O/H ≲ 8.0 (Z/Z⊙ ≲ 1/5) presenting minimal stellar wind features and prominent nebular emission in He ii and C iv. This is consistent with nearly an order of magnitude increase in ionizing photon production beyond the He+-ionizing edge relative to H-ionizing flux as metallicity decreases below a fifth solar, well in excess of standard stellar population synthesis predictions. Our results suggest that often-neglected sources of energetic radiation such as stripped binary products and very massive O-stars produce a sharper change in the ionizing spectrum with decreasing metallicity than expected. Consequently, nebular emission in C iv and He ii powered by these stars may provide useful metallicity constraints in the reionization era.
ABSTRACT
Spitzer/Infrared Array Camera (IRAC) imaging has revealed that the brightest z ∼ 7−8 galaxies often exhibit young ages and strong nebular line emission, hinting at high ionizing efficiency ...among early galaxies. However, IRAC’s limited sensitivity has long hindered efforts to study the fainter, more numerous population often thought largely responsible for reionization. Here, we use Cosmic Evolution Early Release Science (CEERS) JWST/NIRCam data to characterize 116 ultraviolet (UV)-faint (median MUV = −19.5) z ∼ 6.5−8 galaxies. The spectral energy distributions are typically dominated by young (∼10–50 Myr), low-mass (M* ∼ 108 M⊙) stellar populations, and we find no need for extremely high stellar masses (∼1011 M⊙). Considering previous studies of UV-bright (MUV ∼ −22) z ∼ 7−8 galaxies, we find evidence for a strong (5–10 times) increase in specific star formation rate (sSFR) toward lower luminosities (median sSFR = 103 Gyr−1 in CEERS). The larger sSFRs imply a more dominant contribution from OB stars in the relatively numerous UV-faint population, perhaps suggesting that these galaxies are very efficient ionizing agents (median ξion = 1025.7 erg−1 Hz). In spite of the much larger sSFRs, we find little increase in O iii + H β equivalent widths towards fainter MUV (median ≈780 $\mathrm{\mathring{A}}$). If confirmed, this may indicate that a substantial fraction of our CEERS galaxies possess extremely low metallicities (≲3 per cent Z⊙) where O iii emission is suppressed. Alternatively, high ionizing photon escape fractions or bursty star formation histories can also weaken the nebular lines in a subset of our sample. While the majority of galaxies in our sample are very blue (median β = −2.0), we identify a significant tail of very dusty galaxies (β ∼ −1) at ≈0.5$L_\mathrm{UV}^\ast$ which may contribute significantly to the z ∼ 7−8 star formation rate density.
In this paper, we investigate the implications of the integrated galaxy-wide stellar initial mass function (IGIMF) approach in the framework of the semi-analytical model GAEA (GAlaxy Evolution and ...Assembly), which features a detailed treatment of chemical enrichment and stellar feedback. The IGIMF provides an analytic description of the dependence of the stellar IMF shape on the rate of star formation in galaxies. We find that our model with a universal IMF predicts a rather flat a/Fe-stellar mass relation. The model assuming the IGIMF, instead, is able to reproduce the observed increase of a-enhancement with stellar mass, in agreement with previous studies. This is mainly due to the fact that massive galaxies are characterized by larger star formation rates at high redshift, leading to stronger a-enhancement with respect to low-mass galaxies. At the same time, the IGIMF hypothesis does not affect significantly the trend for shorter star formation time-scales for more massive galaxies. We argue that in the IGIMF scenario the a/Fe ratios are good tracers of the highest star formation events. The final stellar masses and mass-to-light ratio of our model massive galaxies are larger than those estimated from the synthetic photometry assuming a universal IMF, providing a self-consistent interpretation of similar recent results, based on dynamical analysis of local early-type galaxies.