We combine data from the Sloan Digital Sky Survey and the Galaxy Evolution Explorer to simultaneously analyze optical spectra and ultraviolet photometry of 231643 galaxies with the starlight spectral ...synthesis code using state-of-the-art stellar population models. We present a new method to estimate GALEX photometry in the SDSS spectroscopic aperture, which proves quite reliable if applied to large samples. In agreement with previous experiments with CALIFA, we find that adding UV constraints leads to a moderate increase on the fraction of ∼ 10 7-10 8 yr populations and a concomitant decrease of younger and older components, yielding slightly older lumi-nosity weighted mean stellar ages. These changes are most relevant in the low-mass end of the blue cloud. An increase in dust attenuation is observed for galaxies dominated by young stars. We investigate the contribution of different stellar populations to the fraction of light in GALEX and SDSS bands across the UV-optical color-magnitude diagram. As an example application, we use this λ dependence to highlight differences between retired galaxies with and without emission lines. In agreement with an independent study by Herpich et al., we find that the former show an excess of intermediate age populations when compared to the later. Finally, we test the suitability of two different prescription for dust, finding that our dataset is best fitted using the attenuation law of starburst galaxies. However, results for the Milky Way extinction curve improve with decreasing τ V , especially for edge-on galaxies.
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, such as O iii/Hβ versus N ii/Hα (the BPT diagram), require the measurement of at least four emission lines. This paper aims to remedy this situation by transposing the usual divisory lines between star-forming (SF) galaxies and active galactic nuclei (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 Sloan Digital Sky Survey Data Release 7 to constitute a suitable sample of 280 000 ELGs, one-third of which are WLGs. (2) Galaxies with strong emission lines are classified using the widely applied criteria of Kewley et al., Kauffmann et al. and Stasińska et al. to distinguish SF galaxies and AGN hosts and Kewley et al. to distinguish Seyferts from LINERs. (3) We transpose these classification schemes to alternative diagrams keeping N ii/Hα as a horizontal axis, but replacing Hβ by a stronger line (Hα or O ii), or substituting the ionization-level sensitive O iii/Hβ ratio with the equivalent width of Hα(WHα). 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 per cent 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. In the course of this analysis, we were led to make the following points. (i) The Kewley et al. BPT line for galaxy classification is generally ill-used. (ii) Replacing O iii/Hβ by WHα in the classification introduces a change in the philosophy of the distinction between LINERs and Seyferts, but not in its results. Because the WHα versus N ii/Hα diagram can be applied to the largest sample of ELGs without loss of discriminating power between Seyferts and LINERs, we recommend its use in further studies. (iii) The dichotomy between Seyferts and LINERs is washed out by WLGs in the BPT plane, but it subsists in other diagnostic diagrams. This suggests that the right wing in the BPT diagram is indeed populated by at least two classes, tentatively identified with bona fide AGN and ‘retired’ galaxies that have stopped forming stars and are ionized by their old stellar populations.
We investigate the cosmic evolution of the absolute and specific star formation rate (SFR, sSFR) of galaxies as derived from a spatially resolved study of the stellar populations in a set of 366 ...nearby galaxies from the Calar Alto Legacy Integral Field Area (CALIFA) survey. The sample spans stellar masses from M⋆ ~ 109 to 1012M⊙ and a wide range of Hubble types. The analysis combines images obtained with the Galaxy Evolution Explorer (GALEX; far-ultraviolet and near-ultraviolet) and Sloan Digital Sky Survey (SDSS; u, g, r, i, z) with the 4000 Å break, Hβ, and MgFe′ indices measured from the CALIFA data cubes to constrain parametric models for the star formation history (SFH), which are then used to study the cosmic evolution of the SFR density (ρSFR), the sSFR, the main sequence of star formation (MSSF), and the stellar mass density (ρ⋆). Several SFH laws are used to fit the observational constrains. A delayed-τ model, SFR ∝ (t0 − t)exp(−(t0 − t)∕τ), provides the best results, in good agreement with those obtained from cosmological surveys. Our main results from this model are that (a) the mass currently in the inner (≤0.5 half-light radius, HLR) regions formed at earlier epochs than the mass in the outer (1–2 HLR) regions of galaxies. The time since the onset of the star formation is longer in the inner regions (t0 ~ 13−10 Gyr) than in the outer ones (t0 ~ 11−9 Gyr) for all the morphologies, while the e-folding timescale τ in the inner region is similar to or shorter than in the outer regions. These results confirm that galaxies of any Hubble type grow inside-out. (b) The sSFR declines rapidly as the Universe evolves, and faster for early- than for late-type galaxies, and for the inner than for the outer regions of galaxies. (c) The evolution of ρSFR and ρ⋆ agrees well with results from cosmological surveys, particularly with the recent results from the Galaxy And Mass Assembly (GAMA), the G10-Cosmological Evolution Survey (COSMOS), and the 3D Hubble Space Telescope (HST) survey. At low redshift, z ≤ 0.5, most star formation takes place in the outer regions of late spiral galaxies, while at z > 2, the inner regions of the progenitors of the current E and S0 are the main contributors to ρSFR. (d) Similarly, the inner regions of galaxies are the main contributor to ρ⋆ at z > 0.5, growing their mass faster than the outer regions, with a lookback time at 50% ρ⋆ of t50 ~ 9 and 6 Gyr for the inner and outer regions. (e) The MSSF follows a power law at high redshift, with the slope evolving with time but always remaining sub-linear, in good agreement with the Illustris simulation. (f) In agreement with galaxy surveys at different redshifts, the average SFH of CALIFA galaxies indicates that galaxies grow their mass mainly in a mode that is well represented by a delayed-τ model, with the peak at z ~ 2 and an e-folding time of ~3.9 Gyr.
We use the W
Hα versus N ii/Hα (WHAN) diagram introduced by us in previous work to provide a comprehensive emission-line classification of Sloan Digital Sky Survey galaxies. This classification is ...able to cope with the large population of weak line galaxies that do not appear in traditional diagrams due 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 low-ionization nuclear emission-line region (LINER) region of traditional diagnostic diagrams. These are galaxies hosting a weakly active galactic nucleus (wAGN) and 'retired galaxies' (RGs), i.e. galaxies that have stopped forming stars and are ionized by their hot low-mass evolved stars.
A useful criterion to distinguish true from fake AGN (i.e. the RGs) is the value of ξ, which measures the ratio of the extinction-corrected Hα luminosity with respect to the Hα luminosity expected from photoionization by stellar populations older than 108 yr. We find that ξ follows a markedly bimodal distribution, with a ξ≫ 1 population composed by systems undergoing star formation and/or nuclear activity, and a peak at ξ∼ 1 corresponding to the prediction of the RG model. We base our classification scheme not on ξ but on a more readily available and model-independent quantity which provides an excellent observational proxy for ξ: the equivalent width of Hα. Based on the bimodal distribution of W
Hα, we set the practical division between wAGN and RGs at W
Hα= 3 Å.
Five classes of galaxies are identified within the WHAN diagram:
pure star-forming galaxies:
and W
Hα > 3 Å;
strong AGN (i.e. Seyferts):
and W
Hα > 6 Å;
weak AGN:
and W
Hα between 3 and 6 Å;
RGs (i.e. fake AGN): W
Hα < 3 Å;
passive galaxies (actually, lineless galaxies): W
Hα and W
N ii < 0.5 Å.
A comparative analysis of star formation histories and of other physical and observational properties in these different classes of galaxies corroborates our proposed differentiation between RGs and wAGN in the LINER-like family. This analysis also shows similarities between strong and weak AGN on the one hand, and retired and passive galaxies on the other.
ABSTRACT
Stellar population studies of early-type galaxies (ETGs) based on their optical stellar continuum suggest that these are quiescent systems. However, emission lines and ultraviolet photometry ...reveal a diverse population. We use a new version of the starlight spectral synthesis code and state-of-the-art stellar population models to simultaneously fit Sloan Digital Sky Survey spectra and Galaxy Evolution Explorer photometry for a sample of 3453 galaxies at z < 0.1 with near ultraviolet (NUV) − r > 5 that are classified as elliptical by Galaxy Zoo. We reproduce far ultraviolet (FUV) magnitudes for 80 per cent of UV upturn galaxies selected using criteria from the literature, suggesting that additional stellar population ingredients such as binaries and extreme horizontal branch stars may have a limited contribution to the UV upturn. The addition of ultraviolet data leads to a broadening of the distributions of mean stellar ages, metallicities, and attenuation. Stellar populations younger than $1\,$ Gyr are required to reproduce the ultraviolet emission in 17 per cent of our sample. These systems represent 43 per cent of the sample at 5 < NUV − r < 5.5 and span the same stellar mass range as other ETGs in our sample. ETGs with young stellar components have larger H α equivalent widths (WH α) and larger dust attenuation. Emission line ratios and WH α indicate that the ionizing source in these systems is a mixture of young and old stellar populations. Their young stellar populations are metal-poor, especially for high-mass galaxies, indicating recent star formation associated with rejuvenation events triggered by external processes, such as minor mergers.
We use the Sloan Digital Sky Survey (SDSS) and Wide-field Sky Survey Explorer surveys to investigate the real nature of galaxies defined as low-ionization nuclear emission-line regions (LINERs) in ...the BPT diagram. After establishing a mid-infrared colour W2–W3 = 2.5 as the optimal separator between galaxies with and without star formation, we investigate the loci of different galaxy classes in the W
H α versus W2–W3 space. We find that: (1) a large fraction of LINER-like galaxies are emission-line retired galaxies, i.e. galaxies which have stopped forming stars and are powered by hot low-mass evolved stars (HOLMES). Their W2–W3 colours show no sign of star formation and their Hα equivalent widths, W
H α, are consistent with ionization by their old stellar populations. (2) Another important fraction have W2–W3 indicative of star formation. This includes objects located in the supposedly ‘pure AGN’ zone of the BPT diagram. (3) A smaller fraction of LINER-like galaxies have no trace of star formation from W2–W3 and a high W
H α, pointing to the presence of an AGN. (4) Finally, a few LINERs tagged as retired by their W
H α but with W2–W3 values indicative of star formation are late-type galaxies whose SDSS spectra cover only the old ‘retired’ bulge. This reinforces the view that LINER-like galaxies are a mixed bag of objects involving different physical phenomena and observational effects thrusted into the same locus of the BPT diagram.
ABSTRACT
We investigate the impact of the diffuse ionized gas (DIG) on abundance determinations in star-forming (SF) galaxies. The DIG is characterized using the H α equivalent width (WH α). From a ...set of 1 409 SF galaxies from the Mapping Nearby Galaxies at APO (MaNGA) survey, we calculate the fractional contribution of the DIG to several emission lines using high-S/N data from SF spaxels (instead of using noisy emission-lines in DIG-dominated spaxels). Our method is applicable to spectra with observed WH α ≳ 10 Å (which are not dominated by DIG emission). Since the DIG contribution depends on galactocentric distance, we provide DIG-correction formulae for both entire galaxies and single aperture spectra. Applying those to a sample of $\, \gt 90\, 000$ SF galaxies from the Sloan Digital Sky Survey, we find the following. (1) The effect of the DIG on strong-line abundances depends on the index used. It is negligible for the (O iii/H β)/(N ii/H α) index, but reaches ∼0.1 dex at the high-metallicity end for N ii/H α. (2) This result is based on the ∼kpc MaNGA resolution, so the real effect of the DIG is likely greater. (3) We revisit the mass–metallicity–star formation rate (SFR) relation by correcting for the DIG contribution in both abundances and SFR. The effect of DIG removal is more prominent at higher stellar masses. Using the N ii/Hα index, O/H increases with SFR at high stellar mass, contrary to previous claims.
Abstract
We use spatially resolved spectroscopy from the Calar Alto Legacy Integral Field Area (CALIFA) survey to study the nature of the line emitting gas in galaxies of different Hubble types, ...focusing on the separation of star-forming (SF) regions from those better characterized as diffuse ionized gas (DIG). The diagnosis is carried out in terms of the equivalent width of H α (WH α). Three nebular regimes are identified. Regions where WH α < 3 Å define what we call the hDIG, the component of the DIG where photoionization is dominated by hot, low-mass, evolved stars. Regions where WH α > 14 Å trace SF complexes. WH α values in the intermediate 3–14 Å range reflect a mixed regime (mDIG) where more than one process contributes. This three-tier scheme is inspired both by theoretical and empirical considerations. Its application to CALIFA galaxies of different types and inclinations leads to the following results: (i) the hDIG component is prevalent throughout ellipticals and S0’s as well as in bulges, and explains the strongly bimodal distribution of WH α both among and within galaxies. (ii) Early-type spirals have some hDIG in their discs, but this component becomes progressively less relevant for later Hubble types. (iii) hDIG emission is also present above and below galactic discs, as seen in several edge-on spirals in our sample. (iv) The SF/mDIG proportion grows steadily from early- to late-type spirals, and from inner to outer radii. (v) Besides circumventing basic inconsistencies in conventional DIG/SF separation criteria based on the H α surface brightness, our WH α-based method produces results in agreement with a classical excitation diagram analysis.
ABSTRACT
The H α and H β emission-line luminosities measured in a single integrated spectrum are affected in non-trivial ways by point-to-point variations in dust attenuation in a galaxy. This work ...investigates the impact of this variation when estimating global H α luminosities corrected for the presence of dust by a global Balmer decrement. Analytical arguments show that the dust-corrected H α luminosity is always underestimated when using the global H α/H β flux ratio to correct for dust attenuation. We measure this effect on 156 face-on star-forming galaxies from the Mapping Nearby Galaxies at APO (MaNGA) survey. At 1–2 kpc spatial resolution, the effect is small but systematic, with the integrated dust-corrected H α luminosity underestimated by 2–4 per cent (and typically not more than by 10 per cent), and depends on the specific star formation rate of the galaxy. Given the spatial resolution of MaNGA, these are lower limits for the effect. From Multi Unit Spectroscopic Explorer (MUSE) observations of NGC 628 with a resolution of 36 pc, we find the discrepancy between the globally and the point-by-point dust-corrected H α luminosity to be 14 ± 1 per cent, which may still underestimate the true effect. We use toy models and simulations to show that the true difference depends strongly on the spatial variance of the H α/H β flux ratio, and on the slope of the relation between H αluminosity and dust attenuation within a galaxy. Larger samples of higher spatial resolution observations are required to quantify the dependence of this effect as a function of galaxy properties.