Revisiting stellar properties of star-forming galaxies with stellar and nebular spectral modelling

Context. Spectral synthesis is a powerful tool for interpreting the physical properties of galaxies by decomposing their spectral energy distributions (SEDs) into the main luminosity contributors (e.g. stellar populations of distinct age and metallicity or ionised gas). However, the impact nebular emission has on the inferred properties of star-forming (SF) galaxies has been largely overlooked over the years, with unknown ramifications to the current understanding of galaxy evolution. Aims. The objective of this work is to estimate the relations between stellar properties (e.g. total mass, mean age, and mean metallicity) of SF galaxies by simultaneously fitting the stellar and nebular continua and comparing them to the results derived through the more common purely stellar spectral synthesis approach. Methods. The main galaxy sample from SDSS DR7 was analysed with two distinct population synthesis codes: FADO, which estimates self-consistently both the stellar and nebular contributions to the SED, and the original version of STARLIGHT, as representative of purely stellar population synthesis codes. Results. Differences between codes regarding average mass, mean age and mean metallicity values can go as high as similar to 0.06 dex for the overall population of galaxies and similar to 0.12 dex for SF galaxies (galaxies with EW(H alpha) > 3 angstrom), with the most prominent difference between both codes in the two populations being in the light-weighted mean stellar age. FADO presents a broader range of mean stellar ages and metallicities for SF galaxies than STARLIGHT, with the latter code preferring metallicity solutions around the solar value (Z(circle dot) = 0.02). A closer look into the average light- and mass-weighted star formation histories of intensively SF galaxies (EW(H alpha) > 75 angstrom) reveals that the light contributions of simple stellar populations (SSPs) younger than <= 10(7) (10(9)) years in STARLIGHT are higher by similar to 5.41% (9.11%) compared to FADO. Moreover, FADO presents higher light contributions from SSPs with metallicity <= Z(circle dot)/200 (Z(circle dot)/50) of around 8.05% (13.51%) when compared with STARLIGHT. This suggests that STARLIGHT is underestimating the average light-weighted age of intensively SF galaxies by up to similar to 0.17 dex and overestimating the light-weighted metallicity by up to similar to 0.13 dex compared to FADO (or vice versa). The comparison between the average stellar properties of passive, SF and intensively SF galaxy samples also reveals that differences between codes increase with increasing EW(H alpha) and decreasing total stellar mass. Moreover, comparing SF results from FADO in a purely stellar mode with the previous results qualitatively suggests that differences between codes are primarily due to mathematical and statistical differences and secondarily due to the impact of the nebular continuum modelling approach (or lack thereof). However, it is challenging to adequately quantify the relative role of each factor since they are likely interconnected. Conclusions. This work finds indirect evidence that a purely stellar population synthesis approach negatively impacts the inferred stellar properties (e.g. mean age and mean metallicity) of galaxies with relatively high star formation rates (e.g. dwarf spirals, 'green peas', and starburst galaxies). In turn, this can bias interpretations of fundamental relations such as the mass-age or mass-metallicity, which are factors worth bearing in mind in light of future high-resolution spectroscopic surveys at higher redshifts (e.g. MOONS and 4MOST-4HS).

Automated summary

This study aims to estimate the relations between stellar properties of star-forming galaxies by simultaneously fitting the stellar and nebular continua and comparing them to the results derived through purely stellar spectral synthesis. The results show differences in average mass, mean age, and mean metallicity values between different population synthesis codes, highlighting the impact of the nebular continuum modeling approach on the inferred properties of galaxies.