Semi-empirical analysis of Sloan Digital Sky Survey galaxies - I. Spectral synthesis method

Abstract

The study of stellar populations in galaxies is entering a new era with the availability of large and high-quality data bases of both observed galactic spectra and state-of-the-art evolutionary synthesis models. In this paper we investigate the power of spectral synthesis as a means to estimate the physical properties of galaxies. Spectral synthesis is nothing more than the decomposition of an observed spectrum in terms of a superposition of a base of simple stellar populations of various ages and metallicities, producing as output the star formation and chemical histories of a galaxy, its extinction and velocity dispersion. Our implementation of this method uses the recent models of Bruzual & Charlot and observed spectra in the 3650–8000 A range. The reliability of this approach is studied by three different means: (1) simulations, (2) comparison with previous work based on a different technique, and (3) analysis of the consistency of results obtained for a sample of galaxies from the Sloan Digital Sky Survey (SDSS). We find that spectral synthesis provides reliable physical parameters as long as one does not attempt a very detailed description of the star formation and chemical histories. Robust and physically interesting parameters are obtained by combining the (individually uncertain) strengths of each simple stellar population in the base. In particular, we show that, besides providing excellent fits to observed galaxy spectra, this method is able to recover useful information on the distributions of stellar ages and, more importantly, stellar metallicities. Stellar masses, velocity dispersion and extinction are also found to be accurately retrieved for realistic signal-to-noise ratios. We apply this synthesis method to a volume-limited sample of 50 362 galaxies from the SDSS Data Release 2, producing a catalogue of stellar population properties. Emission lines are also studied, their measurement being performed after subtracting the computed starlight spectrum from the observed one. A comparison with recent estimates of both observed and physical properties of these galaxies obtained by other groups shows good qualitative and quantitative agreement, despite substantial differences in the methods of analysis. The confidence in the present method is further strengthened by several empirical and astrophysically reasonable correlations between synthesis results and independent quantities. For instance, we report the existence of strong correlations between stellar and nebular metallicities, stellar and nebular extinctions, mean stellar age and equivalent width of Hα and 4000-A break, and between stellar mass and velocity dispersion.