Galactic archaeology based on star counts is instrumental to reconstruct the past mass assembly of Local Group galaxies. The development of new observing techniques and data reduction, coupled with ...the use of sensitive large field of view cameras, now allows us to pursue this technique in more distant galaxies exploiting their diffuse low surface brightness (LSB) light. As part of the ATLAS3D project, we have obtained with the MegaCam camera at the Canada–France–Hawaii Telescope extremely deep, multiband images of nearby early-type galaxies (ETGs). We present here a catalogue of 92 galaxies from the ATLAS3D sample, which are located in low- to medium-density environments. The observing strategy and data reduction pipeline, which achieve a gain of several magnitudes in the limiting surface brightness with respect to classical imaging surveys, are presented. The size and depth of the survey are compared to other recent deep imaging projects. The paper highlights the capability of LSB-optimized surveys at detecting new prominent structures that change the apparent morphology of galaxies. The intrinsic limitations of deep imaging observations are also discussed, among those, the contamination of the stellar haloes of galaxies by extended ghost reflections, and the cirrus emission from Galactic dust. The detection and systematic census of fine structures that trace the present and past mass assembly of ETGs are one of the prime goals of the project. We provide specific examples of each type of observed structures – tidal tails, stellar streams and shells – and explain how they were identified and classified. We give an overview of the initial results. The detailed statistical analysis will be presented in future papers.
We study the volume-limited and nearly mass-selected (stellar mass M
stars 6 × 109 M) ATLAS3D sample of 260 early-type galaxies (ETGs, ellipticals Es and lenticulars S0s). We construct detailed ...axisymmetric dynamical models (Jeans Anisotropic MGE), which allow for orbital anisotropy, include a dark matter halo and reproduce in detail both the galaxy images and the high-quality integral-field stellar kinematics out to about 1R
e, the projected half-light radius. We derive accurate total mass-to-light ratios (M/L)
e
and dark matter fractions f
DM, within a sphere of radius
centred on the galaxies. We also measure the stellar (M/L)stars and derive a median dark matter fraction f
DM = 13 per cent in our sample. We infer masses M
JAM ≡ L × (M/L)
e
2 × M
1/2, where M
1/2 is the total mass within a sphere enclosing half of the galaxy light. We find that the thin two-dimensional subset spanned by galaxies in the
coordinates system, which we call the Mass Plane (MP) has an observed rms scatter of 19 per cent, which implies an intrinsic one of 11 per cent. Here,
is the major axis of an isophote enclosing half of the observed galaxy light, while σ
e
is measured within that isophote. The MP satisfies the scalar virial relation
within our tight errors. This show that the larger scatter in the Fundamental Plane (FP) (L, σ
e
, R
e) is due to stellar population effects including trends in the stellar initial mass function (IMF). It confirms that the FP deviation from the virial exponents is due to a genuine (M/L)
e
variation. However, the details of how both R
e and σ
e
are determined are critical in defining the precise deviation from the virial exponents. The main uncertainty in masses or M/L estimates using the scalar virial relation is in the measurement of R
e. This problem is already relevant for nearby galaxies and may cause significant biases in virial mass and size determinations at high redshift. Dynamical models can eliminate these problems. We revisit the (M/L)
e
-σ
e
relation, which describes most of the deviations between the MP and the FP. The best-fitting relation is
(r band). It provides an upper limit to any systematic increase of the IMF mass normalization with σ
e
. The correlation is more shallow and has smaller scatter for slow rotating systems or for galaxies in Virgo. For the latter, when using the best distance estimates, we observe a scatter in (M/L)
e
of 11 per cent, and infer an intrinsic one of 8 per cent. We perform an accurate empirical study of the link between σ
e
and the galaxies circular velocity V
circ within 1R
e (where stars dominate) and find the relation max (V
circ) 1.76 × σ
e
, which has an observed scatter of 7 per cent. The accurate parameters described in this paper are used in the companion Paper XX (Cappellari et al.) of this series to explore the variation of global galaxy properties, including the IMF, on the projections of the MP.
Much of our knowledge of galaxies comes from analysing the radiation emitted by their stars, which depends on the present number of each type of star in the galaxy. The present number depends on the ...stellar initial mass function (IMF), which describes the distribution of stellar masses when the population formed, and knowledge of it is critical to almost every aspect of galaxy evolution. More than 50 years after the first IMF determination, no consensus has emerged on whether it is universal among different types of galaxies. Previous studies indicated that the IMF and the dark matter fraction in galaxy centres cannot both be universal, but they could not convincingly discriminate between the two possibilities. Only recently were indications found that massive elliptical galaxies may not have the same IMF as the Milky Way. Here we report a study of the two-dimensional stellar kinematics for the large representative ATLAS(3D) sample of nearby early-type galaxies spanning two orders of magnitude in stellar mass, using detailed dynamical models. We find a strong systematic variation in IMF in early-type galaxies as a function of their stellar mass-to-light ratios, producing differences of a factor of up to three in galactic stellar mass. This implies that a galaxy's IMF depends intimately on the galaxy's formation history.
The ATLAS3D project is a multiwavelength survey combined with a theoretical modelling effort. The observations span from the radio to the millimetre and optical, and provide multicolour imaging, ...two-dimensional kinematics of the atomic (H i), molecular (CO) and ionized gas (Hβ, O iii and N i), together with the kinematics and population of the stars (Hβ, Fe5015 and Mg b), for a carefully selected, volume-limited (1.16 × 105 Mpc3) sample of 260 early-type (elliptical E and lenticular S0) galaxies (ETGs). The models include semi-analytic, N-body binary mergers and cosmological simulations of galaxy formation. Here we present the science goals for the project and introduce the galaxy sample and the selection criteria. The sample consists of nearby (D < 42 Mpc, |δ− 29°| < 35°, |b| > 15°) morphologically selected ETGs extracted from a parent sample of 871 galaxies (8 per cent E, 22 per cent S0 and 70 per cent spirals) brighter than MK
< −21.5 mag (stellar mass M
★≳ 6 ×109 M⊙). We analyse possible selection biases and we conclude that the parent sample is essentially complete and statistically representative of the nearby galaxy population. We present the size-luminosity relation for the spirals and ETGs and show that the ETGs in the ATLAS3D sample define a tight red sequence in a colour-magnitude diagram, with few objects in the transition from the blue cloud. We describe the strategy of the SAURON integral field observations and the extraction of the stellar kinematics with the ppxf method. We find typical 1σ errors of ΔV≈ 6 km s−1, Δσ≈ 7 km s−1, Δh
3≈Δh
4≈ 0.03 in the mean velocity, the velocity dispersion and Gauss-Hermite (GH) moments for galaxies with effective dispersion σe≳ 120 km s−1. For galaxies with lower σe (≈40 per cent of the sample) the GH moments are gradually penalized by ppxf towards zero to suppress the noise produced by the spectral undersampling and only V and σ can be measured. We give an overview of the characteristics of the other main data sets already available for our sample and of the ongoing modelling projects.
We use interferometric 12CO(1-0) observations to compare and contrast the extent, surface brightness profiles and kinematics of the molecular gas in CO-rich ATLAS3D early-type galaxies (ETGs) and ...spiral galaxies. We find that the molecular gas extent is smaller in absolute terms in ETGs than in late-type galaxies, but that the size distributions are similar once scaled by the galaxies optical/stellar characteristic scalelengths. Amongst ETGs, we find that the extent of the gas is independent of its kinematic misalignment (with respect to the stars), but does depend on the environment, with Virgo cluster ETGs having less extended molecular gas reservoirs, further emphasizing that cluster ETGs follow different evolutionary pathways from those in the field. Approximately half of ETGs have molecular gas surface brightness profiles that follow the stellar light profile. These systems often have relaxed gas out to large radii, suggesting they are unlikely to have had recent merger/accretion events. A third of the sample galaxies show molecular gas surface brightness profiles that fall off slower than the light, and sometimes show a truncation. These galaxies often have a low mass, and either have disturbed molecular gas or are in the Virgo cluster, suggesting that recent mergers, ram pressure stripping and/or the presence of hot gas can compress/truncate the gas. The remaining galaxies have rings, or composite profiles, that we argue can be caused by the effects of bars. We investigated the kinematics of the molecular gas using position-velocity diagrams, and compared the observed kinematics with dynamical model predictions, and the observed stellar and ionized gas velocities. We confirm that the molecular gas reaches beyond the turnover of the circular velocity curve in 70 per cent of our CO-rich ATLAS3D ETGs, validating previous work on the CO Tully-Fisher relation. In general we find that in most galaxies the molecular gas is dynamically cold, and the observed CO rotation matches well model predictions of the circular velocity. In the galaxies with the largest molecular masses, dust obscuration and/or population gradients can cause model predictions of the circular velocity to disagree with observations of the molecular gas rotation; however, these effects are confined to the most star forming systems. Bars and non-equilibrium conditions can also make the gas deviate from circular orbits. In both these cases, one expects the model circular velocity to be higher than the observed CO velocity, in agreement with our observations. Molecular gas is a better direct tracer of the circular velocity than the ionized gas, justifying its use as a kinematic tracer for Tully-Fisher and similar analyses.
We use the ATLAS3D sample of 260 early-type galaxies to study the apparent kinematic misalignment angle, Ψ, defined as the angle between the photometric and kinematic major axes. We find that 71 per ...cent of nearby early-type galaxies are strictly aligned systems (Ψ≤ 5°), an additional 14 per cent have 5° < Ψ≤ 10° and 90 per cent of galaxies have Ψ≤ 15°. Taking into account measurement uncertainties, 90 per cent of galaxies can be considered aligned to better than 5°, suggesting that only a small fraction of early-type galaxies (∼10 per cent) are not consistent with the axisymmetry within the projected half-light radius. We identify morphological features such as bars and rings (30 per cent), dust structures (16 per cent), blue nuclear colours (6 per cent) and evidence of interactions (8 per cent) visible on ATLAS3D galaxies. We use kinemetry to analyse the mean velocity maps and separate galaxies into two broad types of regular and non-regular rotators. We find 82 per cent of regular rotators and 17 per cent of non-regular rotators, with two galaxies that we were not able to classify due to the poor data quality. The non-regular rotators are typically found in dense regions and are massive. We characterize the specific features in the mean velocity and velocity dispersion maps. The majority of galaxies do not have any specific features, but we highlight here the frequency of the kinematically distinct cores (7 per cent of galaxies) and the aligned double peaks in the velocity dispersion maps (4 per cent of galaxies). We separate galaxies into five kinematic groups based on the kinemetric features, which are then used to interpret the (Ψ-ε) diagram. Most of the galaxies that are misaligned have complex kinematics and are non-regular rotators. In addition, some show evidence of the interaction and might not be in equilibrium, while some are barred. While the trends are weak, there is a tendency that large values of Ψ are found in galaxies at intermediate environmental densities and among the most massive galaxies in the sample. Taking into account the kinematic alignment and the kinemetric analysis, the majority of early-type galaxies have velocity maps more similar to that of the spiral discs than to that of the remnants of equal-mass mergers. We suggest that the most common formation mechanism for early-type galaxies preserves the axisymmetry of the disc progenitors and their general kinematic properties. Less commonly, the formation process results in a triaxial galaxy with much lower net angular momentum.
In Paper I of this series we introduced a volume-limited parent sample of 871 galaxies from which we extracted the ATLAS3D sample of 260 early-type galaxies (ETGs). In Papers II and III we classified ...the ETGs using their stellar kinematics, in a way that is nearly insensitive to the projection effects, and we separated them into fast and slow rotators. Here we look at galaxy morphology and note that the edge-on fast rotators generally are lenticular galaxies. They appear like spiral galaxies with the gas and dust removed, and in some cases are flat ellipticals (E5 or flatter) with discy isophotes. Fast rotators are often barred and span the same full range of bulge fractions as spiral galaxies. The slow rotators are rounder (E4 or rounder, except for counter-rotating discs) and are generally consistent with being genuine, namely spheroidal-like, elliptical galaxies. We propose a revision to the tuning-fork diagram by Hubble as it gives a misleading description of ETGs by ignoring the large variation in the bulge sizes of fast rotators. Motivated by the fact that only one third (34 per cent) of the ellipticals in our sample are slow rotators, we study for the first time the kinematic morphology-density T-Σ relation using fast and slow rotators to replace lenticulars and ellipticals. We find that our relation is cleaner than using classic morphology. Slow rotators are nearly absent at the lowest density environments
per cent and generally constitute a small fraction f(SR) ≈ 4 per cent of the total galaxy population in the relatively low-density environments explored by our survey, with the exception of the densest core of the Virgo cluster f(SR) ≈ 20 per cent. This contrasts with the classic studies that invariably find significant fractions of (misclassified) ellipticals down to the lowest environmental densities. We find a clean log-linear relation between the fraction f(Sp) of spiral galaxies and the local galaxy surface density Σ3, within a cylinder enclosing the three nearest galaxies. This holds for nearly four orders of magnitude in the surface density down to Σ3≈ 0.01 Mpc−2, with f(Sp) decreasing by 10 per cent per dex in Σ3, while f(FR) correspondingly increases. The existence of a smooth kinematic T-Σ relation in the field excludes processes related to the cluster environment, like e.g. ram-pressure stripping, as main contributors to the apparent conversion of spirals into fast rotators in low-density environments. It shows that the segregation is driven by local effects at the small-group scale. This is supported by the relation becoming shallower when using a surface density estimator Σ10 with a cluster scale. Only at the largest densities in the Virgo core does the f(Sp) relation break down and steepen sharply, while the fraction of slow rotators starts to significantly increase. This suggests that a different mechanism is at work there, possibly related to the stripping of the gas from spirals by the hot intergalactic medium in the cluster core and the corresponding lack of cold accretion.
Abstract
We observed 12 nearby H i-detected early-type galaxies (ETGs) of stellar mass ∼1010 M⊙ ≤ M
* ≤ ∼1011 M⊙ with the Mitchell Integral-Field Spectrograph, reaching approximately three half-light ...radii in most cases. We extracted line-of-sight velocity distributions for the stellar and gaseous components. We find little evidence of transitions in the stellar kinematics of the galaxies in our sample beyond the central effective radius, with centrally fast-rotating galaxies remaining fast-rotating and centrally slow-rotating galaxies likewise remaining slow-rotating. This is consistent with these galaxies having not experienced late dry major mergers; however, several of our objects have ionized gas that is misaligned with respect to their stars, suggesting some kind of past interaction. We extract Lick index measurements of the commonly used H β, Fe5015, Mg b, Fe5270 and Fe5335 absorption features, and we find most galaxies to have flat H β gradients and negative Mg b gradients. We measure gradients of age, metallicity and abundance ratio for our galaxies using spectral fitting, and for the majority of our galaxies find negative age and metallicity gradients. We also find the stellar mass-to-light ratios to decrease with radius for most of the galaxies in our sample. Our results are consistent with a view in which intermediate-mass ETGs experience mostly quiet evolutionary histories, but in which many have experienced some kind of gaseous interaction in recent times.
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.