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 present a generalization of surface photometry to the higher-order moments of the line-of-sight velocity distribution of galaxies observed with integral-field spectrographs. The generalization ...follows the approach of surface photometry by determining the best-fitting ellipses along which the profiles of the moments can be extracted and analysed by means of harmonic expansion. The assumption for the odd moments (e.g. mean velocity) is that the profile along an ellipse satisfies a simple cosine law. The assumption for the even moments (e.g. velocity dispersion) is that the profile is constant, as it is used in surface photometry. We test the method on a number of model maps and discuss the meaning of the resulting harmonic terms. We apply the method to the kinematic moments of an axisymmetric model elliptical galaxy and probe the influence of noise on the harmonic terms. We also apply the method to SAURON observations of NGC 2549, NGC 2974, NGC 4459 and NGC 4473 where we detect multiple co- and counter-rotating (NGC 2549 and NGC 4473, respectively) components. We find that velocity profiles extracted along ellipses of early-type galaxies are well represented by the simple cosine law (with 2 per cent accuracy), while possible deviations are carried in the fifth harmonic term which is sensitive to the existence of multiple kinematic components, and has some analogy to the shape parameter of photometry. We compare the properties of the kinematic and photometric ellipses and find that they are often very similar, but a study on a larger sample is necessary. Finally, we offer a characterization of the main velocity structures based only on the kinemetric parameters which can be used to quantify the features in velocity maps.
We make use of interferometric CO and H i observations, and optical integral-field spectroscopy from the ATLAS3D survey, to probe the origin of the molecular and ionized interstellar medium (ISM) in ...local early-type galaxies. We find that 36 ± 5 per cent of our sample of fast-rotating early-type galaxies have their ionized gas kinematically misaligned with respect to the stars, setting a strong lower limit on the importance of externally acquired gas (e.g. from mergers and cold accretion). Slow rotators have a flat distribution of misalignments, indicating that the dominant source of gas is external. The molecular, ionized and atomic gas in all the detected galaxies are always kinematically aligned, even when they are misaligned from the stars, suggesting that all these three phases of the ISM share a common origin. In addition, we find that the origin of the cold and warm gas in fast-rotating early-type galaxies is strongly affected by environment, despite the molecular gas detection rate and mass fractions being fairly independent of group/cluster membership. Galaxies in dense groups and the Virgo cluster nearly always have their molecular gas kinematically aligned with the stellar kinematics, consistent with a purely internal origin (presumably stellar mass loss). In the field, however, kinematic misalignments between the stellar and gaseous components indicate that at least 42 ± 5 per cent of local fast-rotating early-type galaxies have their gas supplied from external sources. When one also considers evidence of accretion present in the galaxies' atomic gas distributions, ≳46 per cent of fast-rotating field ETGs are likely to have acquired a detectable amount of ISM from accretion and mergers. We discuss several scenarios which could explain the environmental dichotomy, including preprocessing in galaxy groups/cluster outskirts and the morphological transformation of spiral galaxies, but we find it difficult to simultaneously explain the kinematic misalignment difference and the constant detection rate. Furthermore, our results suggest that galaxy mass may be an important independent factor associated with the origin of the gas, with the most massive fast-rotating galaxies in our sample (M
K
≲−24 mag; stellar mass of ≈8 × 1010 M⊙) always having kinematically aligned gas. This mass dependence appears to be independent of environment, suggesting it is caused by a separate physical mechanism.
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.
We study the evidence for a diversity of formation processes in early-type galaxies by presenting the first complete volume-limited sample of slow rotators with both integral-field kinematics from ...the ATLAS
3D
Project and high spatial resolution photometry from the
Hubble
Space Telescope. Analysing the nuclear surface brightness profiles of 12 newly imaged slow rotators, we classify their light profiles as core-less, and place an upper limit to the core size of about 10 pc. Considering the full magnitude and volume-limited ATLAS
3D
sample, we correlate the presence or lack of cores with stellar kinematics, including the proxy for the stellar angular momentum (
λ
Re
) and the velocity dispersion within one half-light radius (
σ
e
), stellar mass, stellar age,
α
-element abundance, and age and metallicity gradients. More than half of the slow rotators have core-less light profiles, and they are all less massive than 10
11
M
⊙
. Core-less slow rotators show evidence for counter-rotating flattened structures, have steeper metallicity gradients, and a larger dispersion of gradient values (ΔZ/H¯ = −0.42 ± 0.18) than core slow rotators (ΔZ/H¯ = −0.23 ± 0.07). Our results suggest that core and core-less slow rotators have different assembly processes, where the former, as previously discussed, are the relics of massive dissipation-less merging in the presence of central supermassive black holes. Formation processes of core-less slow rotators are consistent with accretion of counter-rotating gas or gas-rich mergers of special orbital configurations, which lower the final net angular momentum of stars, but support star formation. We also highlight core fast rotators as galaxies that share properties of core slow rotators (i.e. cores, ages,
σ
e
, and population gradients) and core-less slow rotators (i.e. kinematics,
λ
Re
, mass, and larger spread in population gradients). Formation processes similar to those for core-less slow rotators can be invoked to explain the assembly of core fast rotators, with the distinction that these processes form or preserve cores.
ABSTRACT
We present the stellar kinematics of 48 representative elliptical and lenticular galaxies obtained with our custom‐built integral‐field spectrograph SAURON operating on the William Herschel ...Telescope. The data were homogeneously processed through a dedicated reduction and analysis pipeline. All resulting SAURON data cubes were spatially binned to a constant minimum signal‐to‐noise ratio. We have measured the stellar kinematics with an optimized (penalized pixel‐fitting) routine which fits the spectra in pixel space, via the use of optimal templates, and prevents the presence of emission lines to affect the measurements. We have thus generated maps of the mean stellar velocity V, the velocity dispersion σ, and the Gauss–Hermite moments h3 and h4 of the line‐of‐sight velocity distributions. The maps extend to approximately one effective radius. Many objects display kinematic twists, kinematically decoupled components, central stellar discs, and other peculiarities, the nature of which will be discussed in future papers of this series.
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.
Most of the supermassive black hole mass (M•) estimates based on stellar kinematics use the assumption that galaxies are axisymmetric oblate spheroids or spherical. Here, we use fully general ...triaxial orbit-based models to explore the effect of relaxing the axisymmetric assumption on the previously studied galaxies M32 and NGC 3379. We find that M32 can only be modelled accurately using an axisymmetric shape viewed nearly edge-on, and our black hole mass estimate is identical to previous studies. When the observed 5° kinematical twist is included in our model of NGC 3379, the best shape is mildly triaxial, and we find that our best-fitting black hole mass estimate doubles with respect to the axisymmetric model. This particular black hole mass estimate is still within the errors of that of the axisymmetric model and consistent with the M•–σ relationship. However, this effect may have a pronounced impact on black hole demography, since roughly a third of the most massive galaxies are strongly triaxial.
In the past 15 yr, the triaxial Schwarzschild orbit-superposition code developed by van den Bosch and van de Ven in Leiden has been widely applied to study the dynamics of galaxies. Recently, a bug ...was reported in the orbit calculation of this code, specifically in the mirroring procedure that is used to speed up the computation. We have fixed the incorrect mirroring in the DYNAMITE code, which is the publicly-released successor of the Leiden triaxial Schwarzschild code. In this study, we provide a thorough quantification of how this bug has affected the results of dynamical analyses performed with this code. We compare results obtained with the original and corrected versions of DYNAMITE, and discuss the differences in the phase-space distribution of a single orbit and in the global stellar orbit distribution, in the mass estimate of the central black hole in the highly triaxial galaxy PGC 46832, and in the measurement of intrinsic shape and enclosed mass for more than 50 galaxies. Focusing on the typical scientific applications of the Schwarzschild method, in all our tests we find that differences are negligible with respect to the statistical and systematic uncertainties. We conclude that previous results with the Leiden triaxial Schwarzschild code are not significantly affected by the incorrect mirroring.
We combine SAURON integral field data of a representative sample of local early-type, red sequence galaxies with Spitzer/Infrared Array Camera imaging in order to investigate the presence of trace ...star formation in these systems. With the Spitzer data, we identify galaxies hosting low-level star formation, as traced by polycyclic aromatic hydrocarbon emission, with measured star formation rates that compare well to those estimated from other tracers. This star formation proceeds according to established scaling relations with molecular gas content, in surface density regimes characteristic of disc galaxies and circumnuclear starbursts. We find that star formation in early-type galaxies happens exclusively in fast-rotating systems and occurs in two distinct modes. In the first, star formation is a diffuse process, corresponding to widespread young stellar populations and high molecular gas content. The equal presence of co- and counter-rotating components in these systems strongly implies an external origin for the star-forming gas, and we argue that these star formation events may be the final stages of (mostly minor) mergers that build up the bulges of red sequence lenticulars. In the second mode of star formation, the process is concentrated into well-defined disc or ring morphologies, outside of which the host galaxies exhibit uniformly evolved stellar populations. This implies that these star formation events represent rejuvenations within previously quiescent stellar systems. Evidence for earlier star formation events similar to these in all fast-rotating early-type galaxies suggests that this mode of star formation may be common to all such galaxies, with a duty cycle of roughly 1/10, and likely contributes to the embedded, corotating inner stellar discs ubiquitous in this population.