We study the radial structure of the stellar mass surface density ($mu$) and stellar population age as a function of the total stellar mass and morphology for a sample of 107 galaxies from the CALIFA survey. We use the fossil record to recover the st
ar formation history (SFH) in spheroidal and disk dominated galaxies with masses from 10$^9$ to 10$^{12}$ M$_odot$. We derive the half mass radius, and we find that galaxies are on average 15% more compact in mass than in light. HMR/HLR decreases with increasing mass for disk galaxies, but is almost constant in spheroidal galaxies. We find that the galaxy-averaged stellar population age, stellar extinction, and $mu$ are well represented by their values at 1 HLR. Negative radial gradients of the stellar population ages support an inside-out formation. The larger inner age gradients occur in the most massive disk galaxies that have the most prominent bulges; shallower age gradients are obtained in spheroids of similar mass. Disk and spheroidal galaxies show negative $mu$ gradients that steepen with stellar mass. In spheroidal galaxies $mu$ saturates at a critical value that is independent of the galaxy mass. Thus, all the massive spheroidal galaxies have similar local $mu$ at the same radius (in HLR units). The SFH of the regions beyond 1 HLR are well correlated with their local $mu$, and follow the same relation as the galaxy-averaged age and $mu$; suggesting that local stellar mass surface density preserves the SFH of disks. The SFH of bulges are, however, more fundamentally related to the total stellar mass, since the radial structure of the stellar age changes with galaxy mass even though all the spheroid dominated galaxies have similar radial structure in $mu$. Thus, galaxy mass is a more fundamental property in spheroidal systems while the local stellar mass surface density is more important in disks.
In a companion paper we have presented many products derived from the application of the spectral synthesis code STARLIGHT to datacubes from the CALIFA survey, including 2D maps of stellar population properties and 1D averages in the temporal and spa
tial dimensions. Here we evaluate the uncertainties in these products. Uncertainties due to noise and spectral shape calibration errors and to the synthesis method are investigated by means of a suite of simulations based on 1638 CALIFA spectra for NGC 2916, with perturbations amplitudes gauged in terms of the expected errors. A separate study was conducted to assess uncertainties related to the choice of evolutionary synthesis models. We compare results obtained with the Bruzual & Charlot models, a preliminary update of them, and a combination of spectra derived from the Granada and MILES models. About 100k CALIFA spectra are used in this comparison. Noise and shape-related errors at the level expected for CALIFA propagate to 0.10-0.15 dex uncertainties in stellar masses, mean ages and metallicities. Uncertainties in A_V increase from 0.06 mag in the case of random noise to 0.16 mag for shape errors. Higher order products such as SFHs are more uncertain, but still relatively stable. Due to the large number statistics of datacubes, spatial averaging reduces uncertainties while preserving information on the history and structure of stellar populations. Radial profiles of global properties, as well as SFHs averaged over different regions are much more stable than for individual spaxels. Uncertainties related to the choice of base models are larger than those associated with data and method. Differences in mean age, mass and metallicity are ~ 0.15 to 0.25 dex, and 0.1 mag in A_V. Spectral residuals are ~ 1% on average, but with systematic features of up to 4%. The origin of these features is discussed. (Abridged)
Fossil record methods based on spectral synthesis techniques have matured over the past decade, and their application to integrated galaxy spectra fostered substantial advances on the understanding of galaxies and their evolution. Yet, because of the
lack of spatial resolution, these studies are limited to a global view, providing no information about the internal physics of galaxies. Motivated by the CALIFA survey, which is gathering Integral Field Spectroscopy over the full optical extent of 600 galaxies, we have developed an end-to-end pipeline which: (i) partitions the observed data cube into Voronoi zones in order to, when necessary and taking due account of correlated errors, increase the S/N, (ii) extracts spectra, including propagated errors and bad-pixel flags, (iii) feeds the spectra into the STARLIGHT spectral synthesis code, (iv) packs the results for all galaxy zones into a single file, (v) performs a series of post-processing operations, including zone-to-pixel image reconstruction and unpacking the spectral and stellar population properties into multi-dimensional time, metallicity, and spatial coordinates. This paper provides an illustrated description of this whole pipeline and its products. Using data for the nearby spiral NGC 2916 as a show case, we go through each of the steps involved, presenting ways of visualizing and analyzing this manifold. These include 2D maps of properties such as the v-field, stellar extinction, mean ages and metallicities, mass surface densities, star formation rates on different time scales and normalized in different ways, 1D averages in the temporal and spatial dimensions, projections of the stellar light and mass growth (x,y,t) cubes onto radius-age diagrams, etc. The results illustrate the richness of the combination of IFS data with spectral synthesis, providing a glimpse of what is to come from CALIFA and future surveys. (Abridged)
We use the W_Ha versus [NII]/Ha (WHAN) diagram to provide a comprehensive emission-line classification of SDSS galaxies. This classification is able to cope with the large population of weak line galaxies that do not appear in traditional diagrams du
e to a lack of some of the diagnostic lines. A further advantage of the WHAN diagram is to allow the differentiation between two very distinct classes that overlap in the LINER region of traditional diagnostic diagrams. These are galaxies hosting a weakly active nucleus (wAGN) and retired galaxies (RGs), i.e. galaxies that have stopped forming stars and are ionized by their hot evolved low-mass stars. A useful criterion to distinguish true from fake AGN (i.e. the RGs) is the ratio (xi) of the extinction-corrected L_Ha with respect to the Ha luminosity expected from photoionization by stellar populations older than 100 Myr. This ratio follows a markedly bimodal distribution, with a xi >> 1 population composed by systems undergoing star-formation and/or nuclear activity, and a peak at xi ~ 1 corresponding to the prediction of the RG model. We base our classification scheme on the equivalent width of Ha, an excellent observational proxy for xi. Based on the bimodal distribution of W_Ha, we set the division between wAGN and RGs at W_Ha = 3 A. Five classes of galaxies are identified within the WHAN diagram: (a) Pure star forming galaxies: log [NII]/Ha < -0.4 and W_Ha > 3 A. (b) Strong AGN (i.e., Seyferts): log [NII]/Ha > -0.4 and W_Ha > 6 A. (c) Weak AGN: log [NII]/Ha > -0.4 and W_Ha between 3 and 6 A. (d) RGs: W_Ha < 3 A. (e) Passive galaxies (actually, line-less galaxies): W_Ha and W_[NII] < 0.5 A. A comparative analysis of star formation histories and of other properties in these different classes of galaxies corroborates our proposed differentiation between RGs and weak AGN in the LINER-like family. (Abridged)
A numerous population of weak line galaxies (WLGs) is often left out of statistical studies on emission line galaxies (ELGs) due to the absence of an adequate classification scheme, since classical diagnostic diagrams, like [OIII]/Hb vs [NII]/Ha (the
BPT diagram), require the measurement of at least 4 emission lines. This paper aims to remedy this situation by transposing the usual divisory lines between Star Forming (SF) and AGN hosts, and between Seyferts and LINERs to diagrams that are more economical in terms of line quality requirements. By doing this, we rescue from the classification limbo a substantial number of sources and modify the global census of ELGs. More specifically: (1) We use the SDSS DR7 to constitute a suitable sample of 280k ELGs, 1/3 of which are WLGs. (2) Galaxies with strong emission lines are classified using the widely applied criteria of Kewley et al (2001), Kauffmann et al (2003), Stasinska et al (2006) and Kewley et al (2006). (3) We transpose these classification schemes to alternative diagrams keeping [NII]/Ha as a horizontal axis, but replacing Hb by a stronger line (Ha or [OII]), or substituting [OIII]/Hb ratio with the equivalent width of Ha. Optimized equations for the transposed divisory lines are provided. (4) We show that nothing significant is lost in the translation, but that the new diagrams allow one to classify up to 50% more ELGs. (5) Introducing WLGs in the census of galaxies in the local Universe increases the proportion of metal-rich SF galaxies and especially LINERs. (abridged)
Optical studies of starbursts, AGN and their connections usually leave out galaxies whose emission lines are too weak to warrant reliable measurement and classification. Yet, weak line galaxies abound, and deserve a closer look. We show that these ga
laxies are either massive, metal rich star-forming systems, or, more often, LINERs. From our detailed stellar population analysis, we find that these LINERs have stopped forming stars long ago. Moreover, their ionizing radiation field is amazingly consistent with that expected from their old stellar populations alone. The black-hole in the centers of these massive, early-type galaxies is not active enough to overwhelm stellar ionization, and thus, despite their looks, they should not be called AGN.
Retrieving the Star Formation History (SFH) of a galaxy out of its integrated spectrum is the central goal of stellar population synthesis. Recent advances in evolutionary synthesis models have given new breath to this old field of research. Modern s
pectral synthesis techniques incorporating these advances now allow the fitting of galaxy spectra on an angstrom-by-angstrom basis. These detailed fits are useful for a number of studies, like emission line, stellar kinematics, and specially galaxy evolution. Applications of this semi-empirical approach to mega data sets are teaching us a lot about the lives of galaxies. The STARLIGHT spectral synthesis code is one of the tools which allows one to harness this favorable combination of plentifulness of data and models. To illustrate this, we show how SFHs vary across classical emission line diagnostic diagrams. Systematic trends are present along both the star-forming and active-galaxy sequences. We also briefly describe experiments with ne
We study the evolution of 82302 star-forming (SF) galaxies from the SDSS. Our main goals are to explore new ways of handling star formation histories (SFH) obtained with our publicly available spectral synthesis code STARLIGHT, and apply them to inve
stigate how SFHs vary as a function of nebular metallicity (Zneb). Our main results are: (1) A conventional correlation analysis shows how global properties such as luminosity, mass, dust content, mean stellar metallicity and mean stellar age relate to Zneb. (2) We present a simple formalism which compresses the results of the synthesis into time-dependent star formation rates (SFR) and mass assembly histories. (3) The current SFR derived from the population synthesis and that from H-alpha are shown to agree within a factor of two. Thus we now have a way to estimate SFR in AGN hosts, where the H-alpha method cannot be applied. (4) Fully time-dependent SFHs are derived for all galaxies and averaged over six Zneb bins spanning the entire SF wing in the [OIII]/H-beta X [NII]/H-alpha diagram. (5) We find that SFHs vary systematically along the SF sequence, such that low-Zneb systems evolve slower and are currently forming stars at a higher relative rate. (6) At any given time, the distribution of specific SFRs for galaxies within a Zneb-bin is broad and roughly log-normal. (7) The same results are found grouping galaxies in stellar mass (M*) or surface mass density (S*) bins. (8) The overall pattern of SFHs as a function of Zneb, M* or S* is robust against changes in selection criteria, choice of evolutionary synthesis models for the spectral fits, and differential extinction effects. (Abridged)
