No Arabic abstract
Spectral population synthesis (PS) is a fundamental tool in extragalactic research that aims to decipher the assembly history of galaxies from their SED. However, until recently all PS codes were restricted to purely stellar fits, neglecting the essential contribution of nebular emission (NE). With the advent of FADO, the now possible self-consistent modelling of stellar and NE opens new routes to the exploration of galaxy SFHs. The main goal of this study is to quantitatively explore the accuracy to which FADO can recover physical and evolutionary properties of galaxies and compare its output with that from purely stellar PS codes. With this in mind, FADO and STARLIGHT were applied to synthetic SEDs that track the spectral evolution of stars and gas in extinction-free mock galaxies that form their stellar mass ($M_star$) according to different parametric SFHs. Spectral fits were computed for two different set-ups that approximate the spectral range of SDSS and CALIFA data. Our analysis indicates that FADO can recover the key physical and evolutionary properties of galaxies, such as $M_star$ and mass- and light-weighted mean age and metallicity, with an accuracy better than 0.2 dex. This is the case even in phases of strongly elevated sSFR and thus with considerable NE contamination. As for STARLIGHT, our analysis documents a moderately good agreement with theoretical values only for evolutionary phases for which NE drops to low levels. Indeed, fits with STARLIGHT during phases of high sSFR severely overestimate both $M_star$ and the mass-weighted stellar age, whereas strongly underestimate the light-weighted age and metallicity. The insights from this study suggest that the neglect of nebular continuum emission in STARLIGHT and similar purely stellar PS codes could systematically impact $M_star$ and SFH estimates for star-forming galaxies.
The field of galaxy evolution will make a great leap forward in the next decade as a consequence of the huge effort by the scientific community in multi-object spectroscopic facilities. To maximise the impact of such incoming data, the analysis methods must also step up, extracting reliable information from the available spectra. In this paper, we aim to investigate the limits and the reliability of different spectral synthesis methods in the estimation of the mean stellar age and metallicity. The main question this work aims to address is which signal-to-noise ratios (S/N) are needed to reliably determine the mean stellar age and metallicity from a galaxy spectrum and how this depends on the tool used to model the spectra. To address this question we built a set of realistic simulated spectra containing stellar and nebular emission, reproducing the evolution of a galaxy in two limiting cases: a constant star formation rate and an exponentially declining star formation. We degraded the synthetic spectra built from these two star formation histories (SFHs) to different S/N and analysed them with three widely used spectral synthesis codes, namely FADO, STECKMAP, and STARLIGHT. For S/N < 5 all three tools show a large diversity in the results. The FADO and STARLIGHT tools find median differences in the light-weighted mean stellar age of ~0.1 dex, while STECKMAP shows a higher value of ~0.2 dex. Detailed investigations of the best-fit spectrum for galaxies with overestimated mass-weighted quantities point towards the inability of purely stellar models to fit the observed spectra around the Balmer jump. Our results imply that when a galaxy enters a phase of high specific star formation rate the neglect of the nebular continuum emission in the fitting process has a strong impact on the estimation of its SFH when purely stellar fitting codes are used, even in presence of high S/N spectra.
We highlight and discuss the importance of accounting for nebular emission in the SEDs of high redshift galaxies, as lines and continuum emission can contribute significantly or subtly to broad-band photometry. Physical parameters such as the galaxy age, mass, star-formation rate, dust attenuation and others inferred from SED fits can be affected to different extent by the treatment of nebular emission. We analyse a large sample of Lyman break galaxies from z~3-6, and show some main results illustrating e.g. the importance of nebular emission for determinations of the mass-SFR relation, attenuation and age. We suggest that a fairly large scatter in such relations could be intrinsic. We find that the majority of objects (~60-70%) is better fit with SEDs accounting for nebular emission; the remaining galaxies are found to show relatively weak or no emission lines. Our modeling, and supporting empirical evidence, suggests the existence of two categories of galaxies, starbursts and post-starbursts (lower SFR and older galaxies) among the LBG population, and relatively short star-formation timescales.
Accounting for nebular emission when modeling galaxy spectral energy distributions (SEDs) is important, as both line and continuum emission can contribute significantly to the total observed flux. In this work, we present a new nebular emission model integrated within the Flexible Stellar Population Synthesis code that computes the total line and continuum emission for complex stellar populations using the photoionization code Cloudy. The self-consistent coupling of the nebular emission to the matched ionizing spectrum produces emission line intensities that correctly scale with the stellar population as a function of age and metallicity. This more complete model of galaxy SEDs will improve estimates of global gas properties derived with diagnostic diagrams, star formation rates based on H$alpha$, and stellar masses derived from NIR broadband photometry. Our models agree well with results from other photoionization models and are able to reproduce observed emission from H II regions and star-forming galaxies. Our models show improved agreement with the observed H II regions in the Ne III/O II plane and show satisfactory agreement with He II emission from $z=2$ galaxies when including rotating stellar models. Models including post-asymptotic giant branch stars are able to reproduce line ratios consistent with low-ionization emission regions (LIERs).
Galaxy surveys targeting emission lines are characterising the evolution of star-forming galaxies, but there is still little theoretical progress in modelling their physical properties. We predict nebular emission from star-forming galaxies within a cosmological galaxy formation model. Emission lines are computed by combining the semi-analytical model sag with the photoionisation code mapp. We characterise the interstellar medium (ISM) of galaxies by relating the ionisation parameter of gas in galaxies to their cold gas metallicity, obtaining a reasonable agreement with the observed ha, oii, oiii luminosity functions, and the the BPT diagram for local star-forming galaxies. The average ionisation parameter is found to increase towards low star-formation rates and high redshifts, consistent with recent observational results. The predicted link between different emission lines and their associated star-formation rates is studied by presenting scaling relations to relate them. Our model predicts that emission line galaxies have modest clustering bias, and thus reside in dark matter haloes of masses below $M_{rm halo} lesssim 10^{12} {[rm h^{-1} M_{odot}]}$. Finally, we exploit our modelling technique to predict galaxy number counts up to $zsim 10$ by targeting far-infrared (FIR) emission lines detectable with submillimetre facilities
Star forming galaxies exhibit a variety of physical conditions, from quiescent normal spirals to the most powerful dusty starbursts. In order to study these complex systems, we need a suitable tool to analyze the information coming from observations at all wavelengths. We present a new spectro-photometric model which considers in a consistent way starlight as reprocessed by gas and dust. We discuss preliminary results to interpret some observed properties of VLIRGs.