Estimating stellar population and emission line properties in S-PLUS galaxies

TL;DR

This paper introduces a new photometric method using the AlStar code to accurately estimate stellar populations and emission line properties of galaxies from S-PLUS survey data, achieving results comparable to spectral analysis.

Contribution

The study presents an improved empirical prior in AlStar that enhances emission line estimation from photometry, enabling detailed galaxy property recovery without spectroscopy.

Findings

Stellar properties are recovered with better than 0.2 dex accuracy.

Emission line equivalent widths are estimated with dispersions of 0.05–0.2 dex.

Method performs well on both integrated and spatially resolved galaxy data.

Abstract

We present tests of a new method to simultaneously estimate stellar population and emission line (EL) properties of galaxies out of S-PLUS photometry. The technique uses the AlStar code, updated with an empirical prior which greatly improves its ability to estimate ELs using only the survey's 12 bands. The tests compare the output of (noise-perturbed) synthetic photometry of SDSS galaxies to properties derived from previous full spectral fitting and detailed EL analysis. For realistic signal-to-noise ratios, stellar population properties are recovered to better than 0.2 dex in masses, mean ages, metallicities and $\pm 0.2$ mag for the extinction. More importantly, ELs are recovered remarkably well for a photometric survey. We obtain input $-$ output dispersions of 0.05--0.2 dex for the equivalent widths of $[\mathrm{O}\,\rm{II}]$, $[\mathrm{O}\,\rm{III}]$, H$\beta$, H$\alpha$, $[\mathrm{N}\,\rm{II}]$, and $[\mathrm{S}\,\rm{II}]$, and even better for lines stronger than $\sim 5$ $\mathring{A}$. These excellent results are achieved by combining two empirical facts into a prior which restricts the EL space available for the fits: (1) Because, for the redshifts explored here, H$\alpha$ and $[\mathrm{N}\,\rm{II}]$ fall in a single narrow band (J0660), their combined equivalent width is always well recovered, even when $[\mathrm{N}\,\rm{II}]$/H$\alpha$ is not. (2) We know from SDSS that $W_{H\alpha+[\mathrm{N}\,\rm{II}]}$ correlates with $[\mathrm{N}\,\rm{II}]$/H$\alpha$, which can be used to tell if a galaxy belongs to the left or right wings in the classical BPT diagnostic diagram. Example applications to integrated light and spatially resolved data are also presented, including a comparison with independent results obtained with MUSE-based integral field spectroscopy.