We present a stellar population analysis of the absorption line strength maps for 48 early-type galaxies from the SAURON sample. Using the line strength index maps of Hβ, Fe5015 and Mg b, measured in ...the Lick/IDS system and spatially binned to a constant signal-to-noise ratio, together with predictions from up-to-date stellar population models, we estimate the simple stellar population-equivalent (SSP-equivalent) age, metallicity and abundance ratio α/Fe over a two-dimensional field extending up to approximately one effective radius. A discussion of calibrations and differences between model predictions is given. Maps of SSP-equivalent age, metallicity and abundance ratio α/Fe are presented for each galaxy. We find a large range of SSP-equivalent ages in our sample, of which ∼40 per cent of the galaxies show signs of a contribution from a young stellar population. The most extreme cases of post-starburst galaxies, with SSP-equivalent ages of ≤3 Gyr observed over the full field-of-view, and sometimes even showing signs of residual star formation, are restricted to low-mass systems (σe≤ 100 km s−1 or ∼2 × 1010 M⊙). Spatially restricted cases of young stellar populations in circumnuclear regions can almost exclusively be linked to the presence of star formation in a thin, dusty disc/ring, also seen in the near-UV or mid-IR on top of an older underlying stellar population. The flattened components with disc-like kinematics previously identified in all fast rotators are shown to be connected to regions of distinct stellar populations. These range from the young, still star-forming circumnuclear discs and rings with increased metallicity preferentially found in intermediate-mass fast rotators, to apparently old structures with extended disc-like kinematics, which are observed to have an increased metallicity and mildly depressed α/Fe ratio compared to the main body of the galaxy. The slow rotators, often harbouring kinematically decoupled components (KDC) in their central regions, generally show no stellar population signatures over and above the well-known metallicity gradients in early-type galaxies and are largely consistent with old (≥10 Gyr) stellar populations. Using radially averaged stellar population gradients we find in agreement with Spolaor et al. a mass–metallicity gradient relation where low-mass fast rotators form a sequence of increasing metallicity gradient with increasing mass. For more massive systems (above ∼3.5 × 1010 M⊙) there is an overall downturn such that metallicity gradients become shallower with increased scatter at a given mass leading to the most massive systems being slow rotators with relatively shallow metallicity gradients. The observed shallower metallicity gradients and increased scatter could be a consequence of the competition between different star formation and assembly scenarios following a general trend of diminishing gas fractions and more equal-mass mergers with increasing mass, leading to the most massive systems being devoid of ordered motion and signs of recent star formation.
Following our study on the incidence, morphology and kinematics of the ionized gas in early-type galaxies, we now address the question of what is powering the observed nebular emission. To constrain ...the likely sources of gas excitation, we resort to a variety of ancillary data we draw from complementary information on the gas kinematics, stellar populations and galactic potential from the sauron data, and use the sauron-specific diagnostic diagram juxtaposing the O iiiλ5007/Hβ and N iλλ5197, 5200/Hβ line ratios. We find a tight correlation between the stellar surface brightness and the flux of the Hβ recombination line across our sample, which points to a diffuse and old stellar source as the main contributor of ionizing photons in early-type galaxies, with post-asymptotic giant branch (pAGB) stars being still the best candidate based on ionizing balance arguments. The role of AGN photoionization is confined to the central 2–3 arcsec of an handful of objects with radio or X-ray cores. OB-stars are the dominant source of photoionization in 10 per cent of the sauron sample, whereas for another 10 per cent the intense and highly ionized emission is powered by the pAGB population associated to a recently formed stellar subcomponent. Fast shocks are not an important source of ionization for the diffuse nebular emission of early-type galaxies since the required shock velocities can hardly be attained in the potential of our sample galaxies. Finally, in the most massive and slowly or non-rotating galaxies in our sample, which can retain a massive X-ray halo, the finding of a spatial correlation between the hot and warm phases of the interstellar medium (ISM) suggests that the interaction with the hot ISM provides an additional source of ionization besides old ultraviolet-bright stars. This is also supported by a distinct pattern towards lower values of the O iii/Hβ ratio. These results lead us to investigate the relative role of stellar and AGN photoionization in explaining the ionized gas emission observed in early-type galaxies by the Sloan Digital Sky Survey (SDSS). By simulating how our sample galaxies would appear if placed at further distance and targeted by the SDSS, we conclude that only in very few, if any, of the SDSS galaxies which display modest values for the equivalent width of the O iii line (less than ∼2.4 Å) and low-ionization nuclear emission-line region like O iii/Hβ values the nebular emission is truly powered by an AGN.
The hypothesis of a universal initial mass function (IMF) – motivated by observations in nearby stellar systems – has been recently challenged by the discovery of a systematic variation of the IMF ...with the central velocity dispersion, σ, of early-type galaxies (ETGs), towards an excess of low-mass stars in high-σ galaxies. This trend has been derived so far from integrated spectra, and remains unexplained at present. To test whether such trend depends on the local properties within a galaxy, we have obtained new, extremely deep, spectroscopic data, for three nearby ETGs, two galaxies with high σ (∼300 km s− 1), and one lower mass system, with σ ∼ 100 km s− 1. From the analysis of IMF-sensitive spectral features, we find that the IMF depends significantly on galactocentric distance in the massive ETGs, with the enhanced fraction of low-mass stars mostly confined to their central regions. In contrast, the low-σ galaxy does not show any significant radial gradient in the IMF, well described by a shallower distribution, relative to the innermost regions of massive galaxies, at all radii. Such a result indicates that the IMF should be regarded as a local (rather than global) property, and suggests a significant difference between the formation process of the core and the outer regions of massive ETGs.
Abstract
To investigate star formation and assembly processes of massive galaxies, we present here a spatially resolved stellar population analysis of a sample of 45 elliptical galaxies (Es) selected ...from the Calar Alto Legacy Integral Field Area survey. We find rather flat age and Mg/Fe radial gradients, weakly dependent on the effective velocity dispersion of the galaxy within half-light radius. However, our analysis shows that metallicity gradients become steeper with increasing galaxy velocity dispersion. In addition, we have homogeneously compared the stellar population gradients of our sample of Es to a sample of nearby relic galaxies, i.e. local remnants of the high-z population of red nuggets. This comparison indicates that, first, the cores of present-day massive galaxies were likely formed in gas-rich, rapid star formation events at high redshift (z ≳ 2). This led to radial metallicity variations steeper than observed in the local Universe, and positive Mg/Fe gradients. Secondly, our analysis also suggests that a later sequence of minor dry mergers, populating the outskirts of early-type galaxies (ETGs), flattened the pristine Mg/Fe and metallicity gradients. Finally, we find a tight age--Mg/Fe relation, supporting that the duration of the star formation is the main driver of the Mg/Fe enhancement in massive ETGs. However, the star formation time-scale alone is not able to fully explain our Mg/Fe measurements. Interestingly, our results match the expected effect that a variable stellar initial mass function would have on the Mg/Fe ratio.
Two-dimensional stellar kinematics of 48 representative elliptical (E) and lenticular (S0) galaxies obtained with the SAURON integral-field spectrograph reveal that early-type galaxies appear in two ...broad flavours, depending on whether they exhibit clear large-scale rotation or not. We define a new parameter , which involves luminosity-weighted averages over the full two-dimensional kinematic field as a proxy to quantify the observed projected stellar angular momentum per unit mass. We use it as a basis for a new kinematic classification: early-type galaxies are separated into slow and fast rotators, depending on whether they have λR values within their effective radius Re below or above 0.1, respectively. Slow and fast rotators are shown to be physically distinct classes of galaxies, a result which cannot simply be the consequence of a biased viewing angle. Fast rotators tend to be relatively low-luminosity galaxies with MB≳−20.5. Slow rotators tend to be brighter and more massive galaxies, but are still spread over a wide range of absolute magnitude. Three slow rotators of our sample, among the most massive ones, are consistent with zero rotation. Remarkably, all other slow rotators (besides the atypical case of NGC 4550) contain a large kpc-scale kinematically decoupled core (KDC). All fast rotators (except one galaxy with well-known irregular shells) show well-aligned photometric and kinemetric axes, and small velocity twists, in contrast with most slow rotators which exhibit significant misalignments and velocity twists. These results are supported by a supplement of 18 additional early-type galaxies observed with SAURON. In a companion paper (Paper X), we also show that fast and slow rotators are distinct classes in terms of their orbital distribution. We suggest that gas is a key ingredient in the formation and evolution of fast rotators, and that the slowest rotators are the extreme evolutionary end point reached deep in gravitational potential wells where dissipationless mergers had a major role in the evolution, and for which most of the baryonic angular momentum was expelled outwards. Detailed numerical simulations in a cosmological context are required to understand how to form large-scale KDCs within slow rotators, and more generally to explain the distribution of λR values within early-type galaxies and the distinction between fast and slow rotators.
At present, the main challenge to the interpretation of variations in gravity-sensitive line strengths as driven by a non-universal initial mass function (IMF) lies in understanding the effect of the ...other population parameters. Most notably, α/Fe-enhanced populations or even departures in the individual element abundances with respect to the solar-scaled ratio may lead to similar observational results. We combine various TiO-based, IMF-sensitive indicators in the optical and NIR spectral windows, along with the FeH-based Wing–Ford band to break this degeneracy. We obtain a significant radial trend of the IMF slope in XSG1, a massive early-type galaxy
(ETG), with velocity dispersion σ ∼ 300 km s−1, observed with the Very Large Telescope/X-shooter instrument. In addition, we constrain – for the first time – both the shape and normalization of the IMF, using only a stellar population analysis. We robustly rule out a single power law to describe the IMF, whereas a power law tapered off to a constant value at low masses (defined as a bimodal IMF) is consistent with all the observational spectroscopic data and with the stellar M/L constraints based on the Jeans anisotropic modelling method. The IMF in XSG1 is bottom-heavy in the central regions (corresponding to a bimodal IMF slope Γb ∼ 3, or a mass normalization mismatch parameter α ∼ 2), changing towards a standard Milky Way-like IMF (Γb ∼ 1.3; α ∼ 1) at around one half of the effective radius. This result, combined with previous observations of local IMF variations in massive ETGs, reflects the varying processes underlying the formation of the central core and the outer regions in this type of galaxies.
We present the emission-line fluxes and kinematics of 48 representative elliptical and lenticular galaxies obtained with our custom-built integral-field spectrograph, SAURON, operating on the William ...Herschel Telescope. Hβ, O iiiλλ4959,5007 and N iλλ5198,5200 emission lines were measured using a new procedure that simultaneously fits both the stellar spectrum and the emission lines. Using this technique we can detect emission lines down to an equivalent width of 0.1 Å set by the current limitations in describing galaxy spectra with synthetic and real stellar templates, rather than by the quality of our spectra. Gas velocities and velocity dispersions are typically accurate to within 14 and 20 km s−1, respectively, and at worse to within 25 and 40 km s−1. The errors on the flux of the O iii and Hβ lines are on average 10 and 20 per cent, respectively, and never exceed 30 per cent. Emission is clearly detected in 75 per cent of our sample galaxies, and comes in a variety of resolved spatial distributions and kinematic behaviours. A mild dependence on the Hubble type and galactic environment is observed, with higher detection rates in lenticular galaxies and field objects. More significant is the fact that only 55 per cent of the galaxies in the Virgo cluster exhibit clearly detected emission. The ionized-gas kinematics is rarely consistent with simple coplanar circular motions. However, the gas almost never displays completely irregular kinematics, generally showing coherent motions with smooth variations in angular momentum. In the majority of the cases, the gas kinematics is decoupled from the stellar kinematics, and in half of the objects this decoupling implies a recent acquisition of gaseous material. Over the entire sample however, the distribution of the mean misalignment values between stellar and gaseous angular momenta is inconsistent with a purely external origin. The distribution of kinematic misalignment values is found to be strongly dependent on the apparent flattening and the level of rotational support of galaxies, with flatter, fast rotating objects hosting preferentially corotating gaseous and stellar systems. In a third of the cases, the distribution and kinematics of the gas underscore the presence of non-axisymmetric perturbations of the gravitational potential. Consistent with previous studies, the presence of dust features is always accompanied by gas emission while the converse is not always true. A considerable range of values for the O iii/Hβ ratio is found both across the sample and within single galaxies. Despite the limitations of this ratio as an emission-line diagnostic, this finding suggests either that a variety of mechanisms is responsible for the gas excitation in E and S0 galaxies or that the metallicity of the interstellar material is quite heterogeneous.
We investigate the well-known correlations between the dynamical mass-to-light ratio (M/L) and other global observables of elliptical (E) and lenticular (S0) galaxies. We construct two-integral Jeans ...and three-integral Schwarzschild dynamical models for a sample of 25 E/S0 galaxies with SAURON integral-field stellar kinematics to about one effective (half-light) radius Re. They have well-calibrated I-band Hubble Space Telescope WFPC2 and large-field ground-based photometry, accurate surface brightness fluctuation distances, and their observed kinematics is consistent with an axisymmetric intrinsic shape. All these factors result in an unprecedented accuracy in the M/L measurements. We find a tight correlation of the form (M/L) = (3.80 ± 0.14) × (σe/200 km s−1)0.84±0.07 between the M/L (in the I band) measured from the dynamical models and the luminosity-weighted second moment σe of the LOSVD within Re. The observed rms scatter in M/L for our sample is 18 per cent, while the inferred intrinsic scatter is ∼13 per cent. The (M/L)-σe relation can be included in the remarkable series of tight correlations between σe and other galaxy global observables. The comparison of the observed correlations with the predictions of the Fundamental Plane (FP), and with simple virial estimates, shows that the ‘tilt’ of the FP of early-type galaxies, describing the deviation of the FP from the virial relation, is almost exclusively due to a real M/L variation, while structural and orbital non-homology have a negligible effect. When the photometric parameters are determined in the ‘classic’ way, using growth curves, and the σe is measured in a large aperture, the virial mass appears to be a reliable estimator of the mass in the central regions of galaxies, and can be safely used where more ‘expensive’ models are not feasible (e.g. in high-redshift studies). In this case the best-fitting virial relation has the form (M/L)vir = (5.0 ± 0.1) ×Reσ2e/(LG), in reasonable agreement with simple theoretical predictions. We find no difference between the M/L of the galaxies in clusters and in the field. The comparison of the dynamical M/L with the (M/L)pop inferred from the analysis of the stellar population, indicates a median dark matter fraction in early-type galaxies of ∼30 per cent of the total mass inside one Re, in broad agreement with previous studies, and it also shows that the stellar initial mass function varies little among different galaxies. Our results suggest a variation in M/L at constant (M/L)pop, which seems to be linked to the galaxy dynamics. We speculate that fast-rotating galaxies have lower dark matter fractions than the slow-rotating and generally more-massive ones. If correct, this would suggest a connection between the galaxy assembly history and the dark matter halo structure. The tightness of our correlation provides some evidence against cuspy nuclear dark matter profiles in galaxies.
We use the integral-field spectrograph SAURON to measure the stellar line-of-sight velocity distribution and absorption line strengths out to four effective radii (Re) in the early-type galaxies NGC ...3379 and 821. With our newly developed observing technique, we can now probe these faint regions in galaxies that were previously not accessible with traditional long-slit spectroscopy. We make optimal use of the large field-of-view and high throughput of the spectrograph: by adding the signal of all ∼1400 lenslets into one spectrum, we obtain sufficient signal-to-noise in a few hours of observing time to reliably measure the absorption line kinematics and line strengths out to large radius. We find that the line strength gradients previously observed within 1 Re remain constant out to at least 4 Re, which puts constraints on the merger histories of these galaxies. The stellar halo populations are old and metal poor. By constructing orbit-based Schwarzschild dynamical models, we find that dark matter is necessary to explain the observed kinematics in NGC 3379 and 821, with 30–50 per cent of the total matter being dark within 4 Re. The radial anisotropy in our best-fitting halo models is less than in our models without halo, due to differences in orbital structure. The halo also has an effect on the Mg b–Vesc relation: its slope is steeper when a dark matter halo is added to the model.