Alison F. Crocker

Systematic variation of the stellar initial mass function in early-type galaxies

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 ATLAS3D 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 galaxys formation history.

The ATLAS3D project XV: benchmark for early-type galaxies scaling relations from 260 dynamical models: mass-to-light ratio, dark matter, fundamental plane and mass plane

We study the volume-limited and nearly mass-selected (stellar mass M-stars greater than or similar to 6 x 10(9) M circle dot) ATLAS(3D) 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 equivalent to L x (M/L)(e) approximate to 2 x 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 (M-JAM, sigma(e), R-e(maj)) 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 Sigma(e) is measured within that isophote. The MP satisfies the scalar virial relation M-JAM proportional to sigma R-2(e)e(maj) within our tight errors. This show that the larger scatter in the Fundamental Plane (FP) (L, Sigma(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 Sigma(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)-Sigma(e) relation, which describes most of the deviations between the MP and the FP. The best-fitting relation is (M/L)(e) sigma(0.72)(e) (r band). It provides an upper limit to any systematic increase of the IMF mass normalization with Sigma(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 Sigma(e) and the galaxies circular velocity V-circ within 1R(e) (where stars dominate) and find the relation max (V-circ) approximate to 1.76 x Sigma(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.

The ATLAS3D project - XIII. Mass and morphology of H I in early-type galaxies as a function of environment

We present the ATLAS3D H i survey of a volume-limited, complete sample of 166 nearby early-type galaxies (ETGs) brighter than MK=-21.5. The survey is mostly based on data taken with the Westerbork Synthesis Radio Telescope, which enables us to detect H i down to 5 x 1065 x 107 M? within the survey volume. We detect similar to 40 per cent of all ETGs outside the Virgo galaxy cluster and similar to 10 per cent of all ETGs inside it. This demonstrates that it is common for non-cluster ETGs to host H i. The morphology of the detected gas varies in a continuous way from regular, settled H i discs and rings to unsettled gas distributions (including tidal or accretion tails) and systems of clouds scattered around the galaxy. The majority of the detections consist of H i discs or rings (1/4 of all ETGs outside Virgo) so that if H i is detected in an ETG it is most likely distributed on a settled configuration. These systems come in two main types: small discs [ M?], which are confined within the stellar body and share the same kinematics of the stars; and large discs/rings [M(H i) up to 5 x 109 M?], which extend to tens of kpc from the host galaxy and are in half of the cases kinematically decoupled from the stars. Neutral hydrogen seems to provide material for star formation in ETGs. Galaxies containing H i within similar to 1Re exhibit signatures of on-going star formation in similar to 70 per cent of the cases, approximately five times more frequently than galaxies without central H i. The interstellar medium (ISM) in the centre of these galaxies is dominated by molecular gas, and in ETGs with a small gas disc the conversion of H i into H2 is as efficient as in spirals. The ETG H i mass function is characterized by M*similar to 2 x 109 M? and by a slope a similar to-0.7. Compared to spirals, ETGs host much less H i as a family. However, a significant fraction of all ETGs are as H i-rich as spiral galaxies. The main difference between ETGs and spirals is that the former lack the high-column-density H i typical of the bright stellar disc of the latter. The ETG H i properties vary with environment density in a more continuous way than suggested by the known Virgo versus non-Virgo dichotomy. We find an envelope of decreasing M(H i) and M(H i)/LK with increasing environment density. The gas-richest galaxies live in the poorest environments (as found also with CO observations), where the detection rate of star formation signatures is higher. Galaxies in the centre of Virgo have the lowest H i content, while galaxies at the outskirts of Virgo represent a transition region and can contain significant amounts of H i, indicating that at least a fraction of them has joined the cluster only recently after pre-processing in groups. Finally, we find an H i morphologydensity relation such that at low environment density (measured on a local scale) the detected H i is mostly distributed on large, regular discs and rings, while more disturbed H i morphologies dominate environment densities typical of rich groups. This confirms the importance of processes occurring on a galaxy-group scale for the evolution of ETGs.

The CO-to-H2 Conversion Factor and Dust-to-gas Ratio on Kiloparsec Scales in Nearby Galaxies

We present ~kiloparsec spatial resolution maps of the CO-to-H_2 conversion factor (α_(CO)) and dust-to-gas ratio (DGR) in 26 nearby, star-forming galaxies. We have simultaneously solved for α_(CO) and the DGR by assuming that the DGR is approximately constant on kiloparsec scales. With this assumption, we can combine maps of dust mass surface density, CO-integrated intensity, and H I column density to solve for both αCO and the DGR with no assumptions about their value or dependence on metallicity or other parameters. Such a study has just become possible with the availability of high-resolution far-IR maps from the Herschel key program KINGFISH, ^(12)CO J = (2-1) maps from the IRAM 30 m large program HERACLES, and H I 21 cm line maps from THINGS. We use a fixed ratio between the (2-1) and (1-0) lines to present our α_(CO) results on the more typically used ^(12)CO J = (1-0) scale and show using literature measurements that variations in the line ratio do not affect our results. In total, we derive 782 individual solutions for α_(CO) and the DGR. On average, α_(CO) = 3.1 M_☉ pc^(–2) (K km s^(–1))^(–1) for our sample with a standard deviation of 0.3 dex. Within galaxies, we observe a generally flat profile of α_(CO) as a function of galactocentric radius. However, most galaxies exhibit a lower α_(CO) value in the central kiloparsec—a factor of ~2 below the galaxy mean, on average. In some cases, the central α_(CO) value can be factors of 5-10 below the standard Milky Way (MW) value of α_(CO,MW) = 4.4 M_☉ pc^(–2) (K km s^(–1))^(–1). While for α_(CO) we find only weak correlations with metallicity, the DGR is well-correlated with metallicity, with an approximately linear slope. Finally, we present several recommendations for choosing an appropriate α_(CO) for studies of nearby galaxies.

Discovery of an active galactic nucleus driven molecular outflow in the local early-type galaxy NGC 1266

We report the discovery of a powerful molecular wind from the nucleus of the non-interacting nearby S0 field galaxy NGC 1266. The single-dish CO profile exhibits emission to ?400?km?s?1 and requires a nested Gaussian fit to be properly described. Interferometric observations reveal a massive, centrally concentrated molecular component with a mass of 1.1 ? 109 M ? and a molecular outflow with a molecular mass of 2.4 ? 107 M ?. The molecular gas close to the systemic velocity consists of a rotating, compact nucleus with a mass of about 4.1 ? 108 M ? within a radius of 60?pc. This compact molecular nucleus has a surface density of 2.7 ? 104 M ??pc?2, more than two orders of magnitude larger than that of giant molecular clouds in the disk of the Milky Way, and it appears to sit on the Kennicutt-Schmidt relation despite its extreme kinematics and energetic activity. We interpret this nucleus as a disk that confines the outflowing wind. A mass outflow rate of 13 M ? yr?1 leads to a depletion timescale of 85 Myr. The star formation in NGC 1266 is insufficient to drive the outflow, and thus it is likely driven by the active galactic nucleus. The concentration of the majority of the molecular gas in the central 100?pc requires an extraordinary loss of angular momentum, but no obvious companion or interacting galaxy is present to enable the transfer. NGC 1266 is the first known outflowing molecular system that does not show any evidence of a recent interaction.

KINGFISH—Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel: Survey Description and Image Atlas

The KINGFISH project (Key Insights on Nearby Galaxies: a Far-Infrared Survey with Herschel) is an imaging and spectroscopic survey of 61 nearby (d < 30 Mpc) galaxies, chosen to cover a wide range of galaxy properties and local interstellar medium (ISM) environments found in the nearby universe. Its broad goals are to characterize the ISM of present-day galaxies, the heating and cooling of their gaseous and dust components, and to better understand the physical processes linking star formation and the ISM. KINGFISH is a direct descendant of the Spitzer Infrared Nearby Galaxies Survey (SINGS), which produced complete Spitzer imaging and spectroscopic mapping and a comprehensive set of multiwavelength ancillary observations for the sample. The Herschel imaging consists of complete maps for the galaxies at 70, 100, 160, 250, 350, and 500 μm. The spectral line imaging of the principal atomic ISM cooling lines ([O I] 63 μm, [O III] 88 μm, [N II] 122,205 μm, and [C II] 158 μm) covers the subregions in the centers and disks that already have been mapped in the mid-infrared with Spitzer. The KINGFISH and SINGS multiwavelength data sets combined provide panchromatic mapping of the galaxies sufficient to resolve individual star-forming regions, and tracing the important heating and cooling channels of the ISM, across a wide range of local extragalactic ISM environments. This article summarizes the scientific strategy for KINGFISH, the properties of the galaxy sample, the observing strategy, and data processing and products. It also presents a combined Spitzer and Herschel image atlas for the KINGFISH galaxies, covering the wavelength range 3.6–500 μm. All imaging and spectroscopy data products will be released to the Herschel user-generated product archives.

HERSCHEL FAR-INFRARED AND SUBMILLIMETER PHOTOMETRY FOR THE KINGFISH SAMPLE OF NEARBY GALAXIES

New far-infrared and submillimeter photometry from the Herschel Space Observatory is presented for 61 nearby galaxies from the Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel (KINGFISH) sample. The spatially integrated fluxes are largely consistent with expectations based on Spitzer far-infrared photometry and extrapolations to longer wavelengths using popular dust emission models. Dwarf irregular galaxies are notable exceptions, as already noted by other authors, as their 500 μm emission shows evidence for a submillimeter excess. In addition, the fraction of dust heating attributed to intense radiation fields associated with photodissociation regions is found to be (21 ± 4)% larger when Herschel data are included in the analysis. Dust masses obtained from the dust emission models of Draine & Li are found to be on average nearly a factor of two higher than those based on single-temperature modified blackbodies, as single blackbody curves do not capture the full range of dust temperatures inherent to any galaxy. The discrepancy is largest for galaxies exhibiting the coolest far-infrared colors.

The ATLAS3D project – X. On the origin of the molecular and ionized gas in early-type galaxies

We make use of interferometric CO and HI observations, and optical integral-field spectroscopy from the ATLAS 3D survey to probe the origin of the molecular and ionised interstellar medium (ISM) in local early-type galaxies. We find th at 36±5% of our sample of fast rotating early-type galaxies have their ionised gas kinematically misaligned with respect to the stars, setting a strong lower limit on the importance of e xternally acquired gas (e.g. from mergers and cold accretion). Slow rotators have a flat distri bution of misalignments, indicating that the dominant source of gas is external. The molecular, ionised 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 interst ellar medium share a common origin. In addition, we find that the origin of the cold and warm ga s 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 (pres umably stellar mass loss). In the field, however, kinematic misalignments between the stellar and gaseous components indicate that at least 42±5% of local fast-rotating early-type galaxies have their ga s supplied from external sources. When one also considers evidence of accretion present in the galaxies’ atomic gas distributions, > 46% of fast-rotating field ETGs are likely to have acquired a d etectable 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 ga laxies in our sample (MK < 24 mag; stellar mass of �8×10 10 M⊙) always having kinematically aligned gas. This mass dependence appears to be independent of environment, suggesting it is caused by a separate physical mechanism.

The ATLAS3D project - XXV : Two-dimensional kinematic analysis of simulated galaxies and the cosmological origin of fast and slow rotators

We present a detailed two-dimensional stellar dynamical analysis of a sample of 44 cosmological hydrodynamical simulations of individual central galaxies with stellar masses of 2 × 1010 M⊙ ≲ M* ≲ 6 × 1011 M⊙. Kinematic maps of the stellar line-of-sight velocity, velocity dispersion and higher order Gauss-Hermite moments h3 and h4 are constructed for each central galaxy and for the most massive satellites. The amount of rotation is quantified using the λR-parameter. The velocity, velocity dispersion, h3 and h4 fields of the simulated galaxies show a diversity similar to observed kinematic maps of early-type galaxies in the ATLAS3D survey. This includes fast (regular), slow and misaligned rotation, hot spheroids with embedded cold disc components as well as galaxies with counter-rotating cores or central depressions in the velocity dispersion. We link the present-day kinematic properties to the individual cosmological formation histories of the galaxies. In general, major galaxy mergers have a significant influence on the rotation properties resulting in both a spin-down as well as a spin-up of the merger remnant. Lower mass galaxies with significant (≳18 per cent) in situ formation of stars since z ≈ 2, or with additional gas-rich major mergers - resulting in a spin-up - in their formation history, form elongated (ɛ ˜ 0.45) fast rotators (λR ˜ 0.46) with a clear anticorrelation of h3 and v/σ. An additional formation path for fast rotators includes gas-poor major mergers leading to a spin-up of the remnants (λR ˜ 0.43). This formation path does not result in anticorrelated h3 and v/σ. The formation histories of slow rotators can include late major mergers. If the merger is gas rich, the remnant typically is a less flattened slow rotator with a central dip in the velocity dispersion. If the merger is gas poor, the remnant is very elongated (ɛ ˜ 0.43) and slowly rotating (λR ˜ 0.11). The galaxies most consistent with the rare class of non-rotating round early-type galaxies grow by gas-poor minor mergers alone. In general, more massive galaxies have less in situ star formation since z ˜ 2, rotate slower and have older stellar populations. We discuss general implications for the formation of fast and slowly rotating galaxies as well as the weaknesses and strengths of the underlying models.

The ATLAS3D Project – XXX. Star formation histories and stellar population scaling relations of early-type galaxies

We present the stellar population content of early-type galaxies from the ATLAS^(3D) survey. Using spectra integrated within apertures covering up to one effective radius, we apply two methods: one based on measuring line-strength indices and applying single stellar population (SSP) models to derive SSP-equivalent values of stellar age, metallicity, and alpha enhancement; and one based on spectral fitting to derive non-parametric star formation histories, mass-weighted average values of age, metallicity, and half-mass formation time-scales. Using homogeneously derived effective radii and dynamically determined galaxy masses, we present the distribution of stellar population parameters on the Mass Plane (M_(JAM), σ_e, R^(maj)_e), showing that at fixed mass, compact early-type galaxies are on average older, more metal-rich, and more alpha-enhanced than their larger counterparts. From non-parametric star formation histories, we find that the duration of star formation is systematically more extended in lower mass objects. Assuming that our sample represents most of the stellar content of todays local Universe, approximately 50 per cent of all stars formed within the first 2 Gyr following the big bang. Most of these stars reside today in the most massive galaxies (>10^(10.5) M⊙), which themselves formed 90 per cent of their stars by z ∼ 2. The lower mass objects, in contrast, have formed barely half their stars in this time interval. Stellar population properties are independent of environment over two orders of magnitude in local density, varying only with galaxy mass. In the highest density regions of our volume (dominated by the Virgo cluster), galaxies are older, alpha-enhanced, and have shorter star formation histories with respect to lower density regions.