-
Revisiting stellar properties of star-forming galaxies with stellar and nebular spectral modelling
Authors:
Leandro S. M. Cardoso,
Jean Michel Gomes,
Polychronis Papaderos,
Ciro Pappalardo,
Henrique Miranda,
Ana Paulino-Afonso,
José Afonso,
Patricio Lagos
Abstract:
Spectral synthesis is a powerful tool for interpreting the physical properties of galaxies by decomposing their spectral energy distributions into the main luminosity contributors (e.g. stellar populations 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. The objective of this work is to es…
▽ More
Spectral synthesis is a powerful tool for interpreting the physical properties of galaxies by decomposing their spectral energy distributions into the main luminosity contributors (e.g. stellar populations 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. The objective of this work is to estimate the relations between stellar properties of SF galaxies from SDSS DR7 by simultaneously fitting the stellar and nebular continua with FADO and comparing them to the results derived using STARLIGHT, a representative of purely stellar population synthesis codes. Differences between codes regarding average mass, mean age and mean metallicity values can go as high as $\sim$0.06 dex for the overall population of galaxies and $\sim$0.12 dex for SF galaxies (galaxies with EW(H$α$)>3 Å), with the most prominent difference between both codes in the light-weighted mean stellar age. A closer look into the average light- and mass-weighted star formation histories of intensively SF galaxies (EW(H$α$)>75 Å) suggests that STARLIGHT is underestimating the average light-weighted age of intensively SF galaxies by up to $\sim$0.17 dex and overestimating the light-weighted metallicity by up to $\sim$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$α$) and decreasing total stellar mass. This work finds indirect evidence that a purely stellar population synthesis approach negatively impacts the inferred stellar properties of galaxies with relatively high star formation rates. In turn, this can bias interpretations of fundamental relations such as the mass-age or mass-metallicity.
△ Less
Submitted 30 August, 2022;
originally announced August 2022.
-
Self-consistent population spectral synthesis with FADO: II. Star formation history of galaxies in spectral synthesis methods
Authors:
Ciro Pappalardo,
Leandro S. M. Cardoso,
Jean Michel Gomes,
Polychronis Papaderos,
José Afonso,
Iris Breda,
Andrew Humphrey,
Tom Scott,
Stergios Amarantidis,
Israel Matute,
Rodrigo Carvajal,
Silvio Lorenzoni,
Patricio Lagos,
Ana Paulino-Afonso,
Henrique Miranda
Abstract:
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…
▽ More
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.
△ Less
Submitted 15 September, 2021; v1 submitted 17 May, 2021;
originally announced May 2021.
-
Stellar age gradients and inside-out star formation quenching in galaxy bulges
Authors:
Iris Breda,
Polychronis Papaderos,
Jean Michel Gomes,
José Manuel Vílchez,
Bodo L. Ziegler,
Michaela Hirschmann,
Leandro S. M. Cardoso,
Patricio Lagos,
Fernando Buitrago
Abstract:
Radial age gradients hold the cumulative record of the multitude of physical processes driving the build-up of stellar populations and the ensuing star formation (SF) quenching process in galaxy bulges, therefore potentially sensitive discriminators between competing theoretical concepts on bulge formation and evolution. Based on spectral modeling of integral field spectroscopy data from the CALIF…
▽ More
Radial age gradients hold the cumulative record of the multitude of physical processes driving the build-up of stellar populations and the ensuing star formation (SF) quenching process in galaxy bulges, therefore potentially sensitive discriminators between competing theoretical concepts on bulge formation and evolution. Based on spectral modeling of integral field spectroscopy data from the CALIFA survey, we derive mass- and light-weighted stellar age gradients ($\nabla$(t,B)L,M) within the photometrically determined bulge radius (RB) of a representative sample of local face-on late-type galaxies that span 2.6 dex in stellar mass. Our analysis documents a trend for decreasing $\nabla$(t,B)L,M with increasing M,T, with high-mass bulges predominantly showing negative age gradients and vice versa. The inversion from positive to negative $\nabla$(t,B)L,M occurs at logM,T ~ 10, which roughly coincides with the transition from lower-mass bulges whose gas excitation is powered by SF to bulges classified as Composite, LINER or Seyfert. We discuss two limiting cases for the origin of radial age gradients in massive LTG bulges. The first assumes that the stellar age in the bulge is initially spatially uniform, thus the observed age gradients arise from an inside-out SF quenching (ioSFQ) front that is radially expanding with a mean velocity vq. In this case, the age gradients translate into a slow ioSFQ that lasts until z~2, suggesting mild negative feedback by SF or an AGN. If negative age gradients in massive bulges are not due to ioSFQ but primarily due to their inside-out formation process, then the standard hypothesis of quasi-monolithic bulge formation has to be discarded in favor of a scenario that involves gradual buildup of stellar mass over 2-3 Gyr through, e.g., inside-out SF and inward migration of SF clumps from the disk. In this case, rapid AGN-driven ioSFQ cannot be ruled out.
△ Less
Submitted 16 January, 2020;
originally announced January 2020.
-
Self-consistent population spectral synthesis with FADO - I. The importance of nebular emission in modelling star-forming galaxies
Authors:
Leandro S. M. Cardoso,
Jean Michel Gomes,
Polychronis Papaderos
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…
▽ More
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.
△ Less
Submitted 30 January, 2019;
originally announced January 2019.
-
Impact of an AGN featureless continuum on estimation of stellar population properties
Authors:
Leandro S. M. Cardoso,
Jean Michel Gomes,
Polychronis Papaderos
Abstract:
The effect of the featureless power-law (PL) continuum of an active galactic nucleus (AGN) on the estimation of physical properties of galaxies with optical population spectral synthesis (PSS) remains largely unknown. With this in mind, we fit synthetic galaxy spectra representing a wide range of galaxy star formation histories (SFHs) and including distinct PL contributions of the form…
▽ More
The effect of the featureless power-law (PL) continuum of an active galactic nucleus (AGN) on the estimation of physical properties of galaxies with optical population spectral synthesis (PSS) remains largely unknown. With this in mind, we fit synthetic galaxy spectra representing a wide range of galaxy star formation histories (SFHs) and including distinct PL contributions of the form $F_ν \propto ν^{-α}$ with the PSS code STARLIGHT to study to which extent various inferred quantities (e.g. stellar mass, mean age, and mean metallicity) match the input. The synthetic spectral energy distributions (SEDs) computed with our evolutionary spectral synthesis code include an AGN PL component with $0.5 \leq α\leq 2$ and a fractional contribution $0.2 \leq x_{\mathrm{AGN}} \leq 0.8$ to the monochromatic flux at 4020 Å. At the empirical AGN detection threshold $x_{\mathrm{AGN}}\simeq 0.26$ that we previously inferred in a pilot study on this subject, our results show that the neglect of a PL component in spectral fitting can lead to an overestimation by $\sim$2 dex in stellar mass and by up to $\sim$1 and $\sim$4 dex in the light- and mass-weighted mean stellar age, respectively, whereas the light- and mass-weighted mean stellar metallicity are underestimated by up to $\sim$0.3 and $\sim$0.6 dex, respectively. Other fitting set-ups including either a single PL or multiple PLs in the base reveal, on average, much lower unsystematic uncertainties of the order of those typically found when fitting purely stellar SEDs with stellar templates, however, reaching locally up to $\sim$1, 3 and 0.4 dex in mass, age and metallicity, respectively. Our results underscore the importance of an accurate modelling of the AGN spectral contribution in PSS fits as a minimum requirement for the recovery of the physical and evolutionary properties of stellar populations in active galaxies.
△ Less
Submitted 11 May, 2017;
originally announced May 2017.
-
Semi-empirical AGN detection threshold in spectral synthesis studies of Lyman-continuum-leaking early-type galaxies
Authors:
Leandro S. M. Cardoso,
Jean Michel Gomes,
Polychronis Papaderos
Abstract:
Various lines of evidence suggest that the cores of a large portion of early-type galaxies (ETGs) are virtually evacuated of warm ionised gas. This implies that the Lyman-continuum (LyC) radiation produced by an assumed active galactic nucleus (AGN) can escape from the nuclei of these systems without being locally reprocessed into nebular emission, which would prevent their reliable spectroscopic…
▽ More
Various lines of evidence suggest that the cores of a large portion of early-type galaxies (ETGs) are virtually evacuated of warm ionised gas. This implies that the Lyman-continuum (LyC) radiation produced by an assumed active galactic nucleus (AGN) can escape from the nuclei of these systems without being locally reprocessed into nebular emission, which would prevent their reliable spectroscopic classification as Seyfert galaxies with standard diagnostic emission-line ratios. The spectral energy distribution (SED) of these ETGs would then lack nebular emission and be essentially composed of an old stellar component and the featureless power-law (PL) continuum from the AGN. A question that arises in this context is whether the AGN component can be detected with current spectral population synthesis in the optical, specifically, whether these techniques effectively place an AGN detection threshold in LyC-leaking galaxies. To quantitatively address this question, we took a combined approach that involves spectral fitting with STARLIGHT of synthetic SEDs composed of stellar emission that characterises a 10 Gyr old ETG and an AGN power-law component that contributes a fraction $0\leq x_{\mathrm{AGN}} < 1$ of the monochromatic luminosity at $λ_0=$ 4020 Å. In addition to a set of fits for PL distributions $F_ν \propto ν^{-α}$ with the canonical $α=1.5$, we used a base of multiple PLs with $0.5 \leq α\leq 2$ for a grid of synthetic SEDs with a signal-to-noise ratio of 5-$10^3$. Our analysis indicates an effective AGN detection threshold at $x_{\mathrm{AGN}}\simeq 0.26$, which suggests that a considerable fraction of ETGs hosting significant accretion-powered nuclear activity may be missing in the AGN demographics.
△ Less
Submitted 28 September, 2016;
originally announced September 2016.
-
MUSE 3D Spectroscopy and Kinematics of the gigahertz peaked spectrum Radio Galaxy PKS 1934-63: Interaction, Recently Triggered AGN and Star Formation
Authors:
Nathan Roche,
Andrew Humphrey,
Patricio Lagos,
Polychronis Papaderos,
Marckelson Silva,
Leandro S. M. Cardoso,
Jean Michel Gomes
Abstract:
We observe the radio galaxy PKS 1934-63 (at $z=0.1825$) using MUSE (Multi Unit Spectroscopic Explorer) on the Very Large Telescope (VLT). The radio source is GigaHertz Peaked Spectrum and compact (0.13 kpc), implying an early stage of evolution ($\leq 10^4$ yr). Our data show an interacting pair of galaxies, projected separation 9.1 kpc, velocity difference $Δ(v)=216$ km $\rm s^{-1}$. The larger g…
▽ More
We observe the radio galaxy PKS 1934-63 (at $z=0.1825$) using MUSE (Multi Unit Spectroscopic Explorer) on the Very Large Telescope (VLT). The radio source is GigaHertz Peaked Spectrum and compact (0.13 kpc), implying an early stage of evolution ($\leq 10^4$ yr). Our data show an interacting pair of galaxies, projected separation 9.1 kpc, velocity difference $Δ(v)=216$ km $\rm s^{-1}$. The larger galaxy is a $\rm M_{*}\simeq 10^{11}M_{\odot}$ spheroidal with the emission-line spectrum of a high-excitation young radio AGN, e.g. strong [OI]6300 and [OIII]5007. Emission-line ratios indicate a large contribution to the line luminosity from high-velocity shocks ($\simeq 550$ km $\rm s^{-1}$) . The companion is a non-AGN disk galaxy, with extended $\rm Hα$ emission from which its star-formation rate is estimated as $\rm 0.61~M_{\odot}yr^{-1}$.
Both galaxies show rotational velocity gradients in $\rm Hα$ and other lines, with the interaction being prograde-prograde. The SE-NW velocity gradient of the AGN host is misaligned from the E-W radio axis, but aligned with a previously discovered central ultraviolet source, and a factor 2 greater in amplitude in $\rm Hα$ than in other (forbidden) lines (e.g. [OIII]5007). This could be produced by a fast rotating (100-150 km $\rm s^{-1}$) disk with circumnuclear star-formation. We also identify a broad component of [OIII]5007 emission, blueshifted with a velocity gradient aligned with the radio jets, and associated with outflow. However, the broad component of [OI]6300 is redshifted. In spectral fits, both galaxies have old stellar populations plus $\sim 0.1\%$ of very young stars, consistent with the galaxies undergoing first perigalacticon, triggering infall and star-formation from $\sim 40$ Myr ago followed by the radio outburst.
△ Less
Submitted 11 May, 2016; v1 submitted 1 April, 2016;
originally announced April 2016.