Source author record

Eric Emsellem

Eric Emsellem appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

106works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

106 published item(s)

preprint2026arXiv

Distance measurements from the internal dynamics of globular clusters: Application to the Sombrero galaxy (M104)

Globular clusters (GCs) are dense star clusters found in all massive galaxies. Recent work has established that they follow a tight relation between their internal stellar velocity dispersion $σ$ and luminosity, enabling accurate distance measurements. In this work, we aim to apply this GC velocity dispersion (GCVD) distance method to measure the distance to M104 (NGC 4594, the Sombrero galaxy). We have measured internal stellar velocity dispersions for 85 globular clusters (GCs) and one ultra-compact dwarf galaxy around M104 using high-resolution multi-object integrated-light spectroscopy with FLAMES/GIRAFFE on the Very Large Telescope. The measured velocity dispersions range from $σ= 4 - 30$ km s$^{-1}$, with a mean uncertainty of $Δσ= 2.5$ km s$^{-1}$. For a subset of 77 GCs with $V$-band magnitudes and reliable velocity dispersion measurements above $σ> 4$ km s$^{-1}$, we constructed the $M_V$-$σ$ relation to measure the distance to M104, finding $D=9.00\pm0.29$ (stat.)~$\pm0.26$ (sys.) Mpc. The GCs follow the Milky Way and M31 $M_V-σ$ relation closely, with the exception of the luminous ultra-compact dwarf SUCD1, which is nearly one magnitude brighter than the mean relation. 29 GCs in the sample have sizes determined from Hubble Space Telescope imaging which allowed us to determine their masses and $V$-band dynamical mass-to-light ratios (M/L$_V$). We find a mean $<M/L_V> = 2.6 M_{\odot}/L_{\odot}$ for the luminous ($M_V < -8$ mag) M 104 GCs, which is higher than the Milky Way GCs, but is reminiscent of the brightest GCs in Centaurus A. With the exception of SUCD1, the GCs of M104 follow the GCVD relation irrespective of their mass-to-light ratio.

preprint2026arXiv

The impact of ram pressure on the radio spectral index and magnetic field of NGC 4522: A high-resolution VLA continuum study

We present high-resolution Very Large Array (VLA) continuum observations at S-band ($3$ GHz, $560$ pc scale) and X-band ($10$ GHz, $200$ pc scale) of the ram-pressure-stripped Virgo galaxy NGC 4522, to investigate the characteristics of its radio continuum, spectral index, and magnetic field under the influence of the intracluster medium (ICM). The total radio continuum shows an asymmetry that extends northwest, mirroring the HI gas distribution, but showing distinct features in the extraplanar regions. The spectral index steepens systematically from $α\sim-0.6$ in the main disk to $α\sim-1.1$ in the outer disk. We find that the spectral index behavior of the outer disk is mainly due to an ICM shock that can re-accelerate electrons and a significant reduction of thermal emission. Intriguingly, extraplanar clouds exhibit exceptionally flat spectral indices ($α\sim-0.2$ to $0$), resulting from a combination of significantly enhanced thermal emission and pronounced spectral aging of the non-thermal component. Although some of these regions correlate with H$α$, others do not. We propose that the mixing between the ICM and interstellar medium (ISM) is an alternative mechanism that enhances thermal emission independently of star formation. Polarized continuum emissions are highly asymmetric, preferentially distributed along the ICM wind side, and the polarization fraction increases radially outward from the galactic midplane, indicating that the polarized emission is strongly influenced by the ICM wind. Our results show how and where the ICM substantially affects the ISM, and also demonstrate that high-frequency observations are crucial for analyzing the radio continuum of ram pressure stripping galaxies.

preprint2022arXiv

A 2-3 mm high-resolution molecular line survey towards the centre of the nearby spiral galaxy NGC 6946

The complex physical, kinematic, and chemical properties of galaxy centres make them interesting environments to examine with molecular line emission. We present new $2-4$" (${\sim}75{-}150$ pc at $7.7$ Mpc) observations at 2 and 3 mm covering the central $50$" (${\sim}1.9$ kpc) of the nearby double-barred spiral galaxy NGC 6946 obtained with the IRAM Plateau de Bure Interferometer. We detect spectral lines from ten molecules: CO, HCN, HCO$^+$, HNC, CS, HC$_3$N, N$_2$H$^+$, C$_2$H, CH$_3$OH, and H$_2$CO. We complemented these with published 1mm CO observations and 33 GHz continuum observations to explore the star formation rate surface density ${Σ_{\mathrm{SFR}}}$ on 150 pc scales. In this paper, we analyse regions associated with the inner bar of NGC 6946 $-$ the nuclear region (NUC), the northern (NBE), and southern inner bar end (SBE) and we focus on short-spacing corrected bulk (CO) and dense gas tracers (HCN, HCO$^+$, and HNC). We find that HCO$^+$ correlates best with ${Σ_{\mathrm{SFR}}}$, but the dense gas fraction ($f_{\mathrm{dense}}$) and star formation efficiency of the dense gas (${\mathrm{SFE_{dense}}}$) fits show different behaviours than expected from large-scale disc observations.The SBE has a higher ${Σ_{\mathrm{SFR}}}$, $f_{\mathrm{dense}}$, and shocked gas fraction than the NBE. We examine line ratio diagnostics and find a higher CO(2-1)/CO(1-0) ratio towards NBE than for the NUC. Moreover, comparison with existing extragalactic datasets suggests that using the HCN/HNC ratio to probe kinetic temperatures is not suitable on kiloparsec and sub-kiloparsec scales in extragalactic regions. Lastly, our study shows that the HCO$^+$/HCN ratio might not be a unique indicator to diagnose AGN activity in galaxies.

preprint2022arXiv

A tale of two DIGs: The relative role of HII regions and low-mass hot evolved stars in powering the diffuse ionised gas (DIG) in PHANGS-MUSE galaxies

We use integral field spectroscopy from the PHANGS-MUSE survey, which resolves the ionised interstellar medium at ${\sim}50$ pc resolution in 19 nearby spiral galaxies, to study the origin of the diffuse ionised gas (DIG). We examine the physical conditions of the diffuse gas by first removing morphologically-defined HII regions and then binning the low-surface-brightness areas to achieve significant detections of the key nebular lines. A simple model for the leakage and propagation of ionising radiation from HII regions is able to reproduce the observed distribution of H$α$ in the DIG. Leaking radiation from HII regions also explains the observed decrease in line ratios of low-ionisation species ([SII]/H$α$, [NII]/H$α$ and [OI]/H$α$) with increasing H$α$ surface brightness ($Σ_{Hα}$). Emission from hot low-mass evolved stars, however, is required to explain: (1) the enhanced low-ionisation line ratios observed in the central regions of some galaxies; (2) the observed trends of a flat or decreasing [OIII]/H$β$ with $Σ_{Hα}$; and (3) the offset of some DIG regions from the locus of HII regions in the Baldwin-Phillips-Terlevich (BPT) diagram, extending into the area of low-ionisation (nuclear) emission-line regions (LI[N]ERs). Hot low-mass evolved stars make a small contribution to the energy budget of the DIG (2% of the galaxy-integrated H$α$ emission), but their harder spectra make them fundamental contributors to [OIII] emission. The DIG might result from a superposition of two components, an energetically dominant contribution from young stars and a more diffuse background of harder ionising photons from old stars. This unified framework bridges observations of the Milky Way DIG with LI(N)ER-like emission observed in nearby galaxy bulges.

preprint2022arXiv

Stellar feedback in M83 as observed with MUSE -- I. Overview, an unprecedented view of the stellar and gas kinematics and evidence of outflowing gas

We present a large VLT/MUSE mosaic (3.8 x 3.8 kpc) of the nearby spiral galaxy M83, with a spatial resolution ~20 pc. We obtained the kinematics of the stars and ionised gas, and compared them with molecular gas kinematics from ALMA CO(2-1). We separated the ionised gas into HII regions and diffuse ionised gas (DIG) and determined the fraction of Ha luminosity originating from the DIG (f_DIG). We observe that both stars and gas trace the galactic disk rotation, as well as a fast-rotating nuclear component, likely connected to secular processes driven by the galactic bar. In the gas kinematics, we observe a stream east of the nucleus, redshifted with respect to the disk. The stream is surrounded by an extended ionised gas region with enhanced velocity dispersion and a high ionisation state, which is largely consistent with being ionised by slow shocks. We interpret this feature as either the superposition of the disk and an extraplanar layer of DIG, or as a bar-driven inflow of shocked gas. A double Gaussian component fit to the Ha line also reveals the presence of a nuclear biconic structure whose axis of symmetry is perpendicular to the bar. The two cones appear blue- and redshifted along the line of sight and stand out for having an Ha emission separated by up to 200 km s-1 from that of the disk, and a high velocity dispersion ~80-200 km s-1. At the far end of the cones, we observe that the gas is consistent with being ionised by shocks. These features had never been observed before in M83; we postulate that they are tracing a starburst-driven outflow shocking into the surrounding ISM. Finally, we obtain f_DIG ~ 13% in our field of view. We inspect the emission of the HII regions and DIG in `BPT' diagrams, finding that in HII regions photoionisation accounts for 99.8% of the Ha flux, whereas the DIG has a mixed contribution from photoionisation (94.9%) and shocks (5.1%). [abridged]

preprint2022arXiv

The Gas-Star Formation Cycle in Nearby Star-forming Galaxies II. Resolved Distributions of CO and H$α$ Emission for 49 PHANGS Galaxies

The relative distribution of molecular gas and star formation in galaxies gives insight into the physical processes and timescales of the cycle between gas and stars. In this work, we track the relative spatial configuration of CO and H$α$ emission at high resolution in each of our galaxy targets, and use these measurements to quantify the distributions of regions in different evolutionary stages of star formation: from molecular gas without star formation traced by H$α$ to star-forming gas, and to HII regions. The large sample, drawn from the Physics at High Angular resolution in Nearby GalaxieS ALMA and narrowband H$α$ (PHANGS-ALMA and PHANGS-H$α$) surveys, spans a wide range of stellar mass and morphological types, allowing us to investigate the dependencies of the gas-star formation cycle on global galaxy properties. At a resolution of 150 pc, the incidence of regions in different stages shows a dependence on stellar mass and Hubble type of galaxies over the radial range probed. Massive and/or earlier-type galaxies exhibit a significant reservoir of molecular gas without star formation traced by H$α$, while lower-mass galaxies harbor substantial HII regions that may have dispersed their birth clouds or formed from low-mass, more isolated clouds. Galactic structures add a further layer of complexity to relative distribution of CO and H$α$ emission. Trends between galaxy properties and distributions of gas traced by CO and H$α$ are visible only when the observed spatial scale is $\ll$ 500 pc, reflecting the critical resolution requirement to distinguish stages of star formation process.

preprint2022arXiv

The PHANGS-MUSE survey -- Probing the chemo-dynamical evolution of disc galaxies

We present the PHANGS-MUSE survey, a programme using the MUSE IFS at the ESO VLT to map 19 massive $(9.4 < \log(M_{*}/M_\odot) < 11.0)$ nearby (D < 20 Mpc) star-forming disc galaxies. The survey consists of 168 MUSE pointings (1'x1' each), a total of nearly 15 Million spectra, covering ~1.5 Million independent spectra. PHANGS-MUSE provides the first IFS view of star formation across different local environments (including galaxy centres, bars, spiral arms) in external galaxies at a median resolution of 50~pc, better than the mean inter-cloud distance in the ionised interstellar medium. This `cloud-scale' resolution allows detailed demographics and characterisations of HII regions and other ionised nebulae. PHANGS-MUSE further delivers a unique view on the associated gas and stellar kinematics, and provides constraints on the star formation history. The PHANGS-MUSE survey is complemented by dedicated ALMA CO(2-1) and multi-band HST observations, therefore allowing us to probe the key stages of the star formation process from molecular clouds to HII regions and star clusters. This paper describes the scientific motivation, sample selection, observational strategy, data reduction and analysis process of the PHANGS-MUSE survey. We present our bespoke automated data-reduction framework, which is built on the reduction recipes provided by ESO, but additionally allows for mosaicking and homogenisation of the point spread function. We further present a detailed quality assessment and a brief illustration of the potential scientific applications of the large set of PHANGS-MUSE data products generated by our data analysis framework. The data cubes and analysis data products described in this paper represent the basis for the first PHANGS-MUSE public data release and are available in the ESO archive and via the Canadian Astronomy Data Centre.

preprint2021arXiv

PHANGS-HST: Star Cluster Spectral Energy Distribution Fitting with CIGALE

The sensitivity and angular resolution of photometric surveys executed by the Hubble Space Telescope (HST) enable studies of individual star clusters in galaxies out to a few tens of megaparsecs. The fitting of spectral energy distributions (SEDs) of star clusters is essential for measuring their physical properties and studying their evolution. We report on the use of the publicly available Code Investigating GALaxy Emission (CIGALE) SED fitting package to derive ages, stellar masses, and reddenings for star clusters identified in the Physics at High Angular resolution in Nearby GalaxieS-HST (PHANGS-HST) survey. Using samples of star clusters in the galaxy NGC 3351, we present results of benchmark analyses performed to validate the code and a comparison to SED fitting results from the Legacy ExtraGalactic Ultraviolet Survey (LEGUS). We consider procedures for the PHANGS-HST SED fitting pipeline, e.g., the choice of single stellar population models, the treatment of nebular emission and dust, and the use of fluxes versus magnitudes for the SED fitting. We report on the properties of clusters in NGC 3351 and find, on average, the clusters residing in the inner star-forming ring of NGC 3351 are young ($< 10$ Myr) and massive ($10^{5} M_{\odot}$) while clusters in the stellar bulge are significantly older. Cluster mass function fits yield $β$ values around -2, consistent with prior results with a tendency to be shallower at the youngest ages. Finally, we explore a Bayesian analysis with additional physically-motivated priors for the distribution of ages and masses and analyze the resulting cluster distributions.

preprint2021arXiv

PHANGS-MUSE: the HII region luminosity function of local star-forming galaxies

We use an unprecedented sample of about 23,000 HII regions detected at an average physical resolution of 67pc in the PHANGS-MUSE sample to study the extragalactic HII region Ha luminosity function (LF). Our observations probe the star-forming disk of 19 nearby spiral galaxies with low inclination and located close to the star formation main sequence at z=0. The mean LF slope $α$ in our sample is =1.73 with a $σ$ of 0.15. We find that $α$ decreases with the galaxy's star formation rate surface density and argue that this is driven by an enhanced clustering of young stars at high gas surface densities. Looking at the HII regions within single galaxies we find that no significant variations occur between the LF of the inner and outer part of the star-forming disk, whereas the LF in the spiral arm areas is shallower than in the inter-arm areas for six out of the 13 galaxies with clearly visible spiral arms. We attribute these variations to the spiral arms increasing the molecular clouds arm--inter-arm mass contrast and find suggestive evidence that they are more evident for galaxies with stronger spiral arms. Furthermore, we find systematic variations in $α$ between samples of HII regions with high and low ionization parameter q and argue that they are driven by the aging of HII regions.

preprint2021arXiv

The MUSE Extremely Deep Field: Evidence for SFR-induced cores in dark-matter dominated galaxies at z=1

Disc-halo decomposition on rotationally supported star-forming galaxies (SFGs) at $z>1$ are often limited to massive galaxies ($M_\star>10^{10}~M_\odot$) and rely on either deep Integral Field Spectroscopy data or stacking analyses. We present a study of the dark matter (DM) content of nine $z\approx1$ SFGs selected Using the brightest [OII] emitters in the deepest Multi-Unit Spectrograph Explorer (MUSE) field to date, namely the 140hr MUSE Extremely Deep Field, we perform disk-halo decompositions on 9 low-mass SFGs (with $10^{8.5}<M_\star<10^{10.5}~M_\odot$) using a novel 3D modeling approach, which together with the exquisite S/N allows us to measure individual rotation curves to $3\times R_e$. The disk-halo decomposition includes a stellar, DM, gas, and occasionally a bulge component. The DM component primarily uses the generalized $α,β,γ$ profile or a Navarro-Frenk-White (NFW) profile. The disk stellar masses $M_\star$ obtained from the [OII] disk-halo decomposition agree with the values inferred from the spectral energy distributions. While the rotation curves show diverse shapes, ranging from rising to declining at large radii, the DM fractions within the half-light radius $f_{\rm DM}(<R_e)$ are found to be 60\% to 95\%, extending to lower masses (densities) recent results on massive SFGs with $M_\star>10^{10}~M_\odot$. The DM halos show constant surface densities of $\sim100~M_\odot$ pc$^{-2}$. Half of the sample shows a strong preference for cored over cuspy DM profiles. The presence of DM cores appears to be related to galaxies with stellar-to-halo mass $\log M_\star/M_{\rm vir}\approx-2.5$. In addition, the cuspiness of the DM profiles is found to be a strong function of the recent star-formation activity. Both of these results are interpreted as evidence for feedback-induced core formation in the Cold Dark Matter context.

preprint2021arXiv

The PHANGS-HST Survey: Physics at High Angular resolution in Nearby GalaxieS with the Hubble Space Telescope

The PHANGS program is building the first dataset to enable the multi-phase, multi-scale study of star formation across the nearby spiral galaxy population. This effort is enabled by large survey programs with ALMA, VLT/MUSE, and HST, with which we have obtained CO(2-1) imaging, optical spectroscopic mapping, and high resolution UV-optical imaging, respectively. Here, we present PHANGS-HST, which has obtained five band NUV-U-B-V-I imaging of the disks of 38 spiral galaxies at distances of 4-23 Mpc, and parallel V and I band imaging of their halos, to provide a census of tens of thousands of compact star clusters, and multi-scale stellar associations. The combination of HST, ALMA, and VLT/MUSE observations will yield an unprecedented joint catalog of the observed and physical properties of ~100,000 star clusters, associations, HII regions, and molecular clouds. With these basic units of star formation, PHANGS will systematically chart the evolutionary cycling between gas and stars, across a diversity of galactic environments found in nearby galaxies. We discuss the design of the PHANGS-HST survey, and provide an overview of the HST data processing pipeline and first results, highlighting new methods for selecting star cluster candidates, morphological classification of candidates with convolutional neural networks, and identification of stellar associations over a range of physical scales with a watershed algorithm. We describe the cross-observatory imaging, catalogs, and software products to be released. These high-level science products will seed a broad range of investigations, in particular, the study of embedded stellar populations and dust with JWST, for which a PHANGS Cycle 1 Treasury program to obtain eight band 2-21 $μ$m imaging has been approved.

preprint2020arXiv

A spectroscopic study of MATLAS-2019 with MUSE: an ultra-diffuse galaxy with an excess of old globular clusters

The MATLAS deep imaging survey has uncovered a plethora of dwarf galaxies in the low density environment it has mapped. A fraction of them are unusually extended and have a low-surface brightness. Among these so-called ultra-diffuse galaxies, a few seem to host an excess of globular clusters. With the integral-field unit spectrograph MUSE we have observed one of these galaxies - MATLAS J15052031+0148447 (MATLAS-2019) - located towards the nearby group NGC 5846 and measured its systemic velocity,age, and metallicity, and that of its globular clusters candidates. For the stellar body of MATLAS-2019 we derive a metallicity of -1.33+0.19-0.01 dex and an age of 11.2+1.8-0.8 Gyr. For some of the individual GCs and the stacked GC population, we derive consistent ages and metallicities. From the 11 confirmed globular clusters and using a Markov Chain Monte Carlo approach we derived a dynamical mass-to-light ratio of 4.2+8.6-3.4M/L. This is at the lower end of the luminosity-mass scaling relation defined by the Local Group dwarf galaxies. Furthermore, we couldn't confirm nor reject the possibility of a rotational component of the GC system. If present, this would further modify the inferred mass. Follow-up observations of the globular cluster population and of the stellar body of the galaxy are needed to assess whether this galaxy is lacking dark matter like it was suggested for the pair of dwarf galaxies in the field of NGC 1052, or if this is a miss-interpretation arising from systematic uncertainties of the method commonly used for these systems and the large uncertainties of the individual globular cluster velocities.

preprint2020arXiv

Distances to PHANGS Galaxies: New Tip of the Red Giant Branch Measurements and Adopted Distances

PHANGS-HST is an ultraviolet-optical imaging survey of 38 spiral galaxies within ~20 Mpc. Combined with the PHANGS-ALMA, PHANGS-MUSE surveys and other multiwavelength data, the dataset will provide an unprecedented look into the connections between young stars, HII regions, and cold molecular gas in these nearby star-forming galaxies. Accurate distances are needed to transform measured observables into physical parameters (e.g., brightness to luminosity, angular to physical sizes of molecular clouds, star clusters and associations). PHANGS-HST has obtained parallel ACS imaging of the galaxy halos in the F606W and F814W bands. Where possible, we use these parallel fields to derive tip of the red giant branch (TRGB) distances to these galaxies. In this paper, we present TRGB distances for 11 galaxies from ~4 to ~15 Mpc, based on the first year of PHANGS-HST observations. Five of these represent the first published TRGB distance measurements (IC 5332, NGC 2835, NGC 4298, NGC 4321, and NGC 4328), and eight of which are the best available distances to these targets. We also provide a compilation of distances for the 118 galaxies in the full PHANGS sample, which have been adopted for the first PHANGS-ALMA public data release.

preprint2020arXiv

Dynamical Equilibrium in the Molecular ISM in 28 Nearby Star-Forming Galaxies

We compare the observed turbulent pressure in molecular gas, $P_\mathrm{turb}$, to the required pressure for the interstellar gas to stay in equilibrium in the gravitational potential of a galaxy, $P_\mathrm{DE}$. To do this, we combine arcsecond resolution CO data from PHANGS-ALMA with multi-wavelength data that traces the atomic gas, stellar structure, and star formation rate (SFR) for 28 nearby star-forming galaxies. We find that $P_\mathrm{turb}$ correlates with, but almost always exceeds the estimated $P_\mathrm{DE}$ on kiloparsec scales. This indicates that the molecular gas is over-pressurized relative to the large-scale environment. We show that this over-pressurization can be explained by the clumpy nature of molecular gas; a revised estimate of $P_\mathrm{DE}$ on cloud scales, which accounts for molecular gas self-gravity, external gravity, and ambient pressure, agrees well with the observed $P_\mathrm{turb}$ in galaxy disks. We also find that molecular gas with cloud-scale ${P_\mathrm{turb}}\approx{P_\mathrm{DE}}\gtrsim{10^5\,k_\mathrm{B}\,\mathrm{K\,cm^{-3}}}$ in our sample is more likely to be self-gravitating, whereas gas at lower pressure appears more influenced by ambient pressure and/or external gravity. Furthermore, we show that the ratio between $P_\mathrm{turb}$ and the observed SFR surface density, $Σ_\mathrm{SFR}$, is compatible with stellar feedback-driven momentum injection in most cases, while a subset of the regions may show evidence of turbulence driven by additional sources. The correlation between $Σ_\mathrm{SFR}$ and kpc-scale $P_\mathrm{DE}$ in galaxy disks is consistent with the expectation from self-regulated star formation models. Finally, we confirm the empirical correlation between molecular-to-atomic gas ratio and kpc-scale $P_\mathrm{DE}$ reported in previous works.

preprint2020arXiv

Formation channels of slowly rotating early-type galaxies

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 ATLAS3D 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 ~10 pc. Considering the full magnitude and volume-limited ATLAS3D sample, we correlate the presence or lack of cores with stellar kinematics, including the proxy for the stellar angular momentum and the velocity dispersion within one half-light radius, 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}$ Msun. Core-less slow rotators show evidence for counter-rotating flattened structures, have steeper metallicity gradients, and a larger dispersion of gradient values than core slow rotators. Our results suggest that core and core-less slow rotators have different assembly processes, where the former 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 and core-less slow rotators. 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.[Abridged]

preprint2020arXiv

On the accretion of a new group of galaxies onto Virgo: I. Internal kinematics of nine in-falling dEs

Galaxy environment has been shown to play an important role in transforming late-type, star-forming galaxies to quiescent spheroids. This transformation is expected to be more severe for low-mass galaxies (M < $10^{10}$ \(M_\odot\)) in dense galaxy groups and clusters, mostly due to the influence of their past host halos (also known as pre-processing) and their present-day environments. For the first time, in this study, we investigate a sample of nine early-type dwarf galaxies (dEs) that were accreted as a likely bound group onto the Virgo galaxy cluster about 2-3 Gyr ago. Considering this special condition, these nine dEs may provide a test bed for distinguishing between the influence of the Virgo galaxy cluster and the effects of the previous host halo on their current properties. Specifically, we use VLT/MUSE integral-field unit spectra to derive their kinematics and specific angular momentum ($λ_{R}$) profiles. We observe a spread in the $λ_{R}$ profiles of our sample dEs, finding that the $λ_{R}$ profiles of half of them are as high as those of low-mass field galaxies. The remaining dEs exhibit $λ_{R}$ profiles as low as those of Virgo dEs that were likely accreted longer ago. Moreover, we detect nebular emission in one dE with a gas velocity offset suggesting ongoing gas stripping in Virgo. We suggest that the low-$λ_{R}$ dEs in our sample were processed by their previous host halo, prior to their infall to Virgo, and that the high-$λ_{R}$ dEs may be experiencing ram pressure stripping in Virgo.

preprint2020arXiv

PHANGS CO kinematics: disk orientations and rotation curves at 150 pc resolution

We present kinematic orientations and high resolution (150 pc) rotation curves for 67 main sequence star-forming galaxies surveyed in CO (2-1) emission by PHANGS-ALMA. Our measurements are based on the application of a new fitting method tailored to CO velocity fields. Our approach identifies an optimal global orientation as a way to reduce the impact of non-axisymmetric (bar and spiral) features and the uneven spatial sampling characteristic of CO emission in the inner regions of nearby galaxies. The method performs especially well when applied to the large number of independent lines-of-sight contained in the PHANGS CO velocity fields mapped at 1'' resolution. The high resolution rotation curves fitted to these data are sensitive probes of mass distribution in the inner regions of these galaxies. We use the inner slope as well as the amplitude of our fitted rotation curves to demonstrate that CO is a reliable global dynamical mass tracer. From the consistency between photometric orientations from the literature and kinematic orientations determined with our method, we infer that the shapes of stellar disks in the mass range of log($\rm M_{\star}(M_{\odot})$)=9.0-10.9 probed by our sample are very close to circular and have uniform thickness.

preprint2020arXiv

SDSS-IV MaNGA: Spatially resolved star formation in barred galaxies

Bars inhabit the majority of local-Universe disk galaxies and may be important drivers of galaxy evolution through the redistribution of gas and angular momentum within disks. We investigate the star formation and gas properties of bars in galaxies spanning a wide range of masses, environments, and star formation rates using the MaNGA galaxy survey. Using a robustly-defined sample of 684 barred galaxies, we find that fractional (or scaled) bar length correlates with the host's offset from the star-formation main sequence. Considering the morphology of the H$α$ emission we separate barred galaxies into different categories, including barred, ringed, and central configurations, together with H$α$ detected at the ends of a bar. We find that only low-mass galaxies host star formation along their bars, and that this is located predominantly at the leading edge of the bar itself. Our results are supported by recent simulations of massive galaxies, which show that the position of star formation within a bar is regulated by a combination of shear forces, turbulence and gas flows. We conclude that the physical properties of a bar are mostly governed by the existing stellar mass of the host galaxy, but that they also play an important role in the galaxy's ongoing star formation.

preprint2020arXiv

When Gas Dynamics Decouples from Galactic Rotation: Characterizing ISM Circulation in Disk Galaxies

In galactic disks, galactic rotation sets the bulk motion of gas, and its energy and momentum can be transferred toward small scales. Additionally, in the interstellar medium, random and noncircular motions arise from stellar feedback, cloud-cloud interactions, and instabilities, among other processes. Our aim is to comprehend to which extent small-scale gas dynamics is decoupled from galactic rotation. We study the relative contributions of galactic rotation and local noncircular motions to the circulation of gas, $Γ$, a macroscopic measure of local rotation, defined as the line integral of the velocity field around a closed path. We measure the circulation distribution as a function of spatial scale in a set of simulated disk galaxies and we model the velocity field as the sum of galactic rotation and a Gaussian random field. The random field is parameterized by a broken power law in Fourier space, with a break at the scale $λ_c$. We define the spatial scale $λ_{\rm eq}$ at which galactic rotation and non-circular motions contribute equally to $Γ$. For our simulated galaxies, the gas dynamics at the scale of molecular clouds is usually dominated by noncircular motions, but in the center of galactic disks galactic rotation is still relevant. Our model shows that the transfer of rotation from large scales breaks at the scale $λ_c$ and this transition is necessary to reproduce the circulation distribution. We find that $λ_{\rm eq}$, and therefore the structure of the gas velocity field, is set by the local conditions of gravitational stability and stellar feedback.

preprint2019arXiv

Direct Detection of Black Hole-Driven Turbulence in the Centers of Galaxy Clusters

Supermassive black holes (SMBHs) are thought to provide energy that prevents catastrophic cooling in the centers of massive galaxies and galaxy clusters. However, it remains unclear how this "feedback" process operates. We use high-resolution optical data to study the kinematics of multi-phase filamentary structures by measuring the velocity structure function (VSF) of the filaments over a wide range of scales in the centers of three nearby galaxy clusters: Perseus, Abell 2597 and Virgo. We find that the motions of the filaments are turbulent in all three clusters studied. There is a clear correlation between features of the VSFs and the sizes of bubbles inflated by SMBH driven jets. Our study demonstrates that SMBHs are the main driver of turbulent gas motions in the centers of galaxy clusters and suggests that this turbulence is an important channel for coupling feedback to the environment. Our measured amplitude of turbulence is in good agreement with Hitomi Doppler line broadening measurement and X-ray surface brightness fluctuation analysis, suggesting that the motion of the cold filaments is well-coupled to that of the hot gas. The smallest scales we probe are comparable to the mean free path in the intracluster medium (ICM). Our direct detection of turbulence on these scales provides the clearest evidence to date that isotropic viscosity is suppressed in the weakly-collisional, magnetized intracluster plasma.

preprint2019arXiv

Mapping electron temperature variations across a spiral arm in NGC 1672

We report one of the first extragalactic observations of electron temperature variations across a spiral arm. Using MUSE mosaic observations of the nearby galaxy NGC 1672, we measure the [N II]5755 auroral line in a sample of 80 HII regions in the eastern spiral arm of NGC1672. We discover systematic temperature variations as a function of distance perpendicular to the spiral arm. The electron temperature is lowest on the spiral arm itself and highest on the downstream side. Photoionization models of different metallicity, pressure, and age of the ionizing source are explored to understand what properties of the interstellar medium drive the observed temperature variations. An azimuthally varying metallicity appears to be the most likely cause of the temperature variations. The electron temperature measurements solidify recent discoveries of azimuthal variations of oxygen abundance based on strong lines, and rule out the possibility that the abundance variations are artefacts of the strong-line calibrations.

preprint2019arXiv

The Gas-Star Formation Cycle in Nearby Star-Forming Galaxies I. Assessment of Multi-scale Variations

The processes regulating star formation in galaxies are thought to act across a hierarchy of spatial scales. To connect extragalactic star formation relations from global and kpc-scale measurements to recent cloud-scale resolution studies, we have developed a simple, robust method that quantifies the scale dependence of the relative spatial distributions of molecular gas and recent star formation. In this paper, we apply this method to eight galaxies with roughly 1 arcsec resolution molecular gas imaging from the PHANGS-ALMA and PAWS surveys that have matched resolution, high quality narrowband Halpha imaging. At a common scale of 140pc, our massive (log(Mstar/Msun)=9.3-10.7), normally star-forming (SFR/Msun/yr=0.3-5.9) galaxies exhibit a significant reservoir of quiescent molecular gas not associated with star formation as traced by Halpha emission. Galactic structures act as backbones for both molecular and HII region distributions. As we degrade the spatial resolution, the quiescent molecular gas disappears, with the most rapid changes occurring for resolutions up to about 0.5kpc. As the resolution becomes poorer, the morphological features become indistinct for spatial scales larger than about 1kpc. The method is a promising tool to search for relationships between the quiescent or star-forming molecular reservoir and galaxy properties, but requires a larger sample size to identify robust correlations between the star-forming molecular gas fraction and global galaxy parameters.

preprint2019arXiv

The lifecycle of molecular clouds in nearby star-forming disc galaxies

It remains a major challenge to derive a theory of cloud-scale ($\lesssim100$ pc) star formation and feedback, describing how galaxies convert gas into stars as a function of the galactic environment. Progress has been hampered by a lack of robust empirical constraints on the giant molecular cloud (GMC) lifecycle. We address this problem by systematically applying a new statistical method for measuring the evolutionary timeline of the GMC lifecycle, star formation, and feedback to a sample of nine nearby disc galaxies, observed as part of the PHANGS-ALMA survey. We measure the spatially-resolved ($\sim100$ pc) CO-to-H$α$ flux ratio and find a universal de-correlation between molecular gas and young stars on GMC scales, allowing us to quantify the underlying evolutionary timeline. GMC lifetimes are short, typically 10-30 Myr, and exhibit environmental variation, between and within galaxies. At kpc-scale molecular gas surface densities $Σ_{\rm H_2}\geqslant8$M$_{\odot}$pc$^{-2}$, the GMC lifetime correlates with time-scales for galactic dynamical processes, whereas at $Σ_{\rm H_2}\leqslant8$M$_{\odot}$pc$^{-2}$ GMCs decouple from galactic dynamics and live for an internal dynamical time-scale. After a long inert phase without massive star formation traced by H$α$ (75-90% of the cloud lifetime), GMCs disperse within just 1-5 Myr once massive stars emerge. The dispersal is most likely due to early stellar feedback, causing GMCs to achieve integrated star formation efficiencies of 4-10% These results show that galactic star formation is governed by cloud-scale, environmentally-dependent, dynamical processes driving rapid evolutionary cycling. GMCs and HII regions are the fundamental units undergoing these lifecycles, with mean separations of 100-300 pc in star-forming discs. Future work should characterise the multi-scale physics and mass flows driving these lifecycles.

preprint2019arXiv

The MUSE Atlas of Disks (MAD): Ionized gas kinematic maps and an application to Diffuse Ionized Gas

We have obtained data for 41 star forming galaxies in the MUSE Atlas of Disks (MAD) survey with VLT/MUSE. These data allow us, at high resolution of a few 100 pc, to extract ionized gas kinematics ($V, σ$) of the centers of nearby star forming galaxies spanning 3 dex in stellar mass. This paper outlines the methodology for measuring the ionized gas kinematics, which we will use in subsequent papers of this survey. We also show how the maps can be used to study the kinematics of diffuse ionized gas for galaxies of various inclinations and masses. Using two different methods to identify the diffuse ionized gas, we measure rotation velocities of this gas for a subsample of 6 galaxies. We find that the diffuse ionized gas rotates on average slower than the star forming gas with lags of 0-10 km/s while also having higher velocity dispersion. The magnitude of these lags is on average 5 km/s lower than observed velocity lags between ionized and molecular gas. Using Jeans models to interpret the lags in rotation velocity and the increase in velocity dispersion we show that most of the diffuse ionized gas kinematics are consistent with its emission originating from a somewhat thicker layer than the star forming gas, with a scale height that is lower than that of the stellar disk.

preprint2016arXiv

A low upper mass limit for the central black hole in the late-type galaxy NGC 4414

We present our mass estimate of the central black hole in the isolated spiral galaxy NGC 4414. Using natural guide star adaptive optics assisted observations with the Gemini Near-Infrared Integral Field Spectrometer (NIFS) and the natural seeing Gemini Multi-Object Spectrographs-North (GMOS), we derived two-dimensional stellar kinematic maps of NGC 4414 covering the central 1.5 arcsec and 10 arcsec, respectively, at a NIFS spatial resolution of 0.13 arcsec. The kinematic maps reveal a regular rotation pattern and a central velocity dispersion dip down to around 105 km/s. We constructed dynamical methods using two different methods: Jeans anisotropic dynamical modeling and axisymmetric Schwarzschild modeling. Both modeling methods give consistent results, but we cannot constrain the lower mass limit and only measure an upper limit for the black hole mass of Mbh= 1.56 x 10^6 Msun(at 3 sigma level) which is at least 1 sigma below the recent Mbh-sigma_e relations. Further tests with dark matter, mass-to-light ratio variation and different light models confirm that our results are not dominated by uncertainties. The derived upper limit mass is not only below the Mbh-sigma_e relation, but is also five times lower than the lower limit black hole mass anticipated from the resolution limit of the sphere of influence. This proves that via high quality integral field data we are now able to push black hole measurements down to at least five times less than the resolution limit.

preprint2016arXiv

Exploring the mass assembly of the early-type disc galaxy NGC3115 with MUSE

We present MUSE integral field spectroscopic data of the S0 galaxy NGC 3115 obtained during the instrument commissioning at the ESO Very Large Telescope (VLT). We analyse the galaxy stellar kinematics and stellar populations and present two dimensional maps of their associated quantities. We thus illustrate the capacity of MUSE to map extra-galactic sources to large radii in an efficient manner, i.e., ~4 Re, and provide relevant constraints on its mass assembly. We probe the well known set of substructures of NGC 3115 (its nuclear disc, stellar rings, outer kpc-scale stellar disc and spheroid) and show their individual associated signatures in the MUSE stellar kinematics and stellar populations maps. In particular, we confirm that NGC 3115 has a thin fast rotating stellar disc embedded in a fast rotating spheroid, and that these two structures show clear differences in their stellar age and metallicity properties. We emphasise an observed correlation between the radial stellar velocity, V, and the Gauss-Hermite moment, h3, creating a "butterfly" shape in the central 15" of the h3 map. We further detect the previously reported weak spiral and ring-like structures, and find evidence that these features can be associated with regions of younger mean stellar ages. We provide tentative evidence for the presence of a bar, despite the fact that the V-h3 correlation can be reproduced by a simple axisymmetric dynamical model. Finally, we present a reconstruction of the two dimensional star formation history of NGC 3115 and find that most of its current stellar mass was formed at early epochs (>12Gyr ago), while star formation continued in the outer (kpc-scale) stellar disc until recently. Since z~2, and within ~4 Re, we suggest that NGC 3115 has been mainly shaped by secular processes.

preprint2016arXiv

MUSE crowded field 3D spectroscopy of over 12,000 stars in the globular cluster NGC 6397 - II. Probing the internal dynamics and the presence of a central black hole

We present a detailed analysis of the kinematics of the Galactic globular cluster NGC 6397 based on more than ~18,000 spectra obtained with the novel integral field spectrograph MUSE. While NGC 6397 is often considered a core collapse cluster, our analysis suggests a flattening of the surface brightness profile at the smallest radii. Although it is among the nearest globular clusters, the low velocity dispersion of NGC 6397 of <5km/s imposes heavy demands on the quality of the kinematical data. We show that despite its limited spectral resolution, MUSE reaches an accuracy of 1km/s in the analysis of stellar spectra. We find slight evidence for a rotational component in the cluster and the velocity dispersion profile that we obtain shows a mild central cusp. To investigate the nature of this feature, we calculate spherical Jeans models and compare these models to our kinematical data. This comparison shows that if a constant mass-to-light ratio is assumed, the addition of an intermediate-mass black hole with a mass of 600M_sun brings the model predictions into agreement with our data, and therefore could be at the origin of the velocity dispersion profile. We further investigate cases with varying mass-to-light ratios and find that a compact dark stellar component can also explain our observations. However, such a component would closely resemble the black hole from the constant mass-to-light ratio models as this component must be confined to the central ~5arcsec of the cluster and must have a similar mass. Independent constraints on the distribution of stellar remnants in the cluster or kinematic measurements at the highest possible spatial resolution should be able to distinguish the two alternatives.

preprint2016arXiv

New insights on the formation of nuclear star clusters

Nuclear Clusters (NCs) are common stellar systems in the centres of galaxies. Yet, the physical mechanisms involved in their formation are still debated. Using a parsec-resolution hydrodynamical simulation of a dwarf galaxy, we propose an updated formation scenario for NCs. In this 'wet migration scenario', a massive star cluster forms in the gas-rich disc, keeping a gas reservoir, and growing further while it migrates to the centre via a combination of interactions with other substructures and dynamical friction. A wet merger with another dense cluster and its own gas reservoir can occur, although this is not a pre-requisite for the actual formation of the NC. The merging process does significantly alter the properties of the NC (mass, morphology, star formation history), also quenching the on-going local star formation activity, thus leading to interesting observational diagnostics for the physical origin of NCs. A population of lower mass clusters co-exist during the simulation, but these are either destroyed via tidal forces, or have high angular momentum preventing them to interact with the NC and contribute to its growth. The proposed updated scenario emphasises the role of gas reservoirs associated with the densest star clusters formed in a gas-rich low-mass galaxy.

preprint2016arXiv

SDSS IV MaNGA - Spatially resolved diagnostic diagrams: A proof that many galaxies are LIERs

We study the spatially resolved excitation properties of the ionised gas in a sample of 646 galaxies using integral field spectroscopy data from SDSS-IV MaNGA. Making use of Baldwin-Philips-Terlevich diagnostic diagrams we demonstrate the ubiquitous presence of extended (kpc scale) low ionisation emission-line regions (LIERs) in both star forming and quiescent galaxies. In star forming galaxies LIER emission can be associated with diffuse ionised gas, most evident as extra-planar emission in edge-on systems. In addition, we identify two main classes of galaxies displaying LIER emission: `central LIER' (cLIER) galaxies, where central LIER emission is spatially extended, but accompanied by star formation at larger galactocentric distances, and `extended LIER' (eLIER) galaxies, where LIER emission is extended throughout the whole galaxy. In eLIER and cLIER galaxies, LIER emission is associated with radially flat, low H$α$ equivalent width of line emission ($<$ 3 Å) and stellar population indices demonstrating the lack of young stellar populations, implying that line emission follows tightly the continuum due to the underlying old stellar population. The H$α$ surface brightness radial profiles are always shallower than $\rm 1/r^{2}$ and the line ratio [OIII]$λ$5007/[OII]$λ$3727,29 (a tracer of the ionisation parameter of the gas) shows a flat gradient. This combined evidence strongly supports the scenario in which LIER emission is not due to a central point source but to diffuse stellar sources, the most likely candidates being hot, evolved (post-asymptotic giant branch) stars. Shocks are observed to play a significant role in the ionisation of the gas only in rare merging and interacting systems.

preprint2016arXiv

SDSS-IV MaNGA IFS Galaxy Survey --- Survey Design, Execution, and Initial Data Quality

The MaNGA Survey (Mapping Nearby Galaxies at Apache Point Observatory) is one of three core programs in the Sloan Digital Sky Survey IV. It is obtaining integral field spectroscopy (IFS) for 10K nearby galaxies at a spectral resolution of R~2000 from 3,622-10,354A. The design of the survey is driven by a set of science requirements on the precision of estimates of the following properties: star formation rate surface density, gas metallicity, stellar population age, metallicity, and abundance ratio, and their gradients; stellar and gas kinematics; and enclosed gravitational mass as a function of radius. We describe how these science requirements set the depth of the observations and dictate sample selection. The majority of targeted galaxies are selected to ensure uniform spatial coverage in units of effective radius (Re) while maximizing spatial resolution. About 2/3 of the sample is covered out to 1.5Re (Primary sample), and 1/3 of the sample is covered to 2.5Re (Secondary sample). We describe the survey execution with details that would be useful in the design of similar future surveys. We also present statistics on the achieved data quality, specifically, the point spread function, sampling uniformity, spectral resolution, sky subtraction, and flux calibration. For our Primary sample, the median r-band signal-to-noise ratio is ~73 per 1.4A pixel for spectra stacked between 1-1.5 Re. Measurements of various galaxy properties from the first year data show that we are meeting or exceeding the defined requirements for the majority of our science goals.

preprint2016arXiv

The Atlas3D project -- XXXI. Nuclear radio emission in nearby early-type galaxies

We present the results of a high-resolution, 5 GHz, Karl G. Jansky Very Large Array study of the nuclear radio emission in a representative subset of the Atlas3D survey of early-type galaxies (ETGs). We find that 51 +/- 4% of the ETGs in our sample contain nuclear radio emission with luminosities as low as 10^18 W/Hz. Most of the nuclear radio sources have compact (< 25-110 pc) morphologies, although < 10% display multi-component core+jet or extended jet/lobe structures. Based on the radio continuum properties, as well as optical emission line diagnostics and the nuclear X-ray properties, we conclude that the majority of the central 5 GHz sources detected in the Atlas3D galaxies are associated with the presence of an active galactic nucleus (AGN). However, even at sub-arcsecond spatial resolution, the nuclear radio emission in some cases appears to arise from low-level nuclear star formation rather than an AGN, particularly when molecular gas and a young central stellar population is present. This is in contrast to popular assumptions in the literature that the presence of a compact, unresolved, nuclear radio continuum source universally signifies the presence of an AGN. Additionally, we examine the relationships between the 5 GHz luminosity and various galaxy properties including the molecular gas mass and - for the first time - the global kinematic state. We discuss implications for the growth, triggering, and fueling of radio AGNs, as well as AGN-driven feedback in the continued evolution of nearby ETGs.

preprint2016arXiv

The growth of the central region by acquisition of counter-rotating gas in star-forming galaxies

Galaxies grow through both internal and external processes. In about 10% of nearby red galaxies with little star formation, gas and stars are counter-rotating, demonstrating the importance of external gas acquisition in these galaxies. However, systematic studies of such phenomena in blue, star-forming galaxies are rare, leaving uncertain the role of external gas acquisition in driving evolution of blue galaxies. Based on new measurements with integral field spectroscopy of a large representative galaxy sample, we find an appreciable fraction of counter-rotators among blue galaxies (9 out of 489 galaxies). The central regions of blue counter-rotators show younger stellar populations and more intense, ongoing star formation than their outer parts, indicating ongoing growth of the central regions. The result offers observational evidence that the acquisition of external gas in blue galaxies is possible; the interaction with pre-existing gas funnels the gas into nuclear regions (< 1 kpc) to form new stars.

preprint2016arXiv

The Next Generation Virgo Cluster Survey (NGVS). XIII. The Luminosity and Mass Function of Galaxies in the Core of the Virgo Cluster and the Contribution from Disrupted Satellites

We present measurements of the galaxy luminosity and stellar mass function in a 3.71 deg$^2$ (0.3 Mpc$^2$) area in the core of the Virgo cluster, based on $ugriz$ data from the Next Generation Virgo Cluster Survey (NGVS). The galaxy sample consists of 352 objects brighter than $M_g=-9.13$ mag, the 50% completeness limit of the survey. Using a Bayesian analysis, we find a best-fit faint end slope of $α=-1.33 \pm 0.02$ for the g-band luminosity function; consistent results are found for the stellar mass function as well as the luminosity function in the other four NGVS bandpasses. We discuss the implications for the faint-end slope of adding 92 ultra compact dwarfs galaxies (UCDs) -- previously compiled by the NGVS in this region -- to the galaxy sample, assuming that UCDs are the stripped remnants of nucleated dwarf galaxies. Under this assumption, the slope of the luminosity function (down to the UCD faint magnitude limit, $M_g = -9.6$ mag) increases dramatically, up to $α= -1.60 \pm 0.06$ when correcting for the expected number of disrupted non-nucleated galaxies. We also calculate the total number of UCDs and globular clusters that may have been deposited in the core of Virgo due to the disruption of satellites, both nucleated and non-nucleated. We estimate that ~150 objects with $M_g\lesssim-9.6$ mag and that are currently classified as globular clusters, might, in fact, be the nuclei of disrupted galaxies. We further estimate that as many as 40% of the (mostly blue) globular clusters in the core of Virgo might once have belonged to such satellites; these same disrupted satellites might have contributed ~40% of the total luminosity in galaxies observed in the core region today. Finally, we use an updated Local Group galaxy catalog to provide a new measurement of the luminosity function of Local Group satellites, $α=-1.21\pm0.05$.

preprint2016arXiv

The Next Generation Virgo Cluster Survey XVI. The Angular Momentum of Dwarf Early-Type Galaxies from Globular Cluster Satellites

We analyze the kinematics of six Virgo cluster dwarf early-type galaxies (dEs) from their globular cluster (GC) systems. We present new Keck/DEIMOS spectroscopy for three of them and reanalyze the data found in the literature for the remaining three. We use two independent methods to estimate the rotation amplitude (Vmax) and velocity dispersion (sigma_GC) of the GC systems and evaluate their statistical significance by simulating non-rotating GC systems with the same number of GC satellites and velocity uncertainties. Our measured kinematics agree with the published values for the three galaxies from the literature and, in all cases, some rotation is measured. However, our simulations show that the null hypothesis of being non-rotating GC systems cannot be ruled out. In the case of VCC1861, the measured Vmax and the simulations indicate that it is not rotating. In the case of VCC1528, the null hypothesis can be marginally ruled out, thus, it might be rotating although further confirmation is needed. In our analysis, we find that, in general, the measured Vmax tends to be overestimated and the measured sigma_GC tends to be underestimated by amounts that depend on the intrinsic Vmax/sigma_GC, the number of observed GCs (N_GC), and the velocity uncertainties. The bias is negligible when N_GC>~20. In those cases where a large N_GC is not available, it is imperative to obtain data with small velocity uncertainties. For instance, errors of <2km/s lead to Vmax<10km/s for a system that is intrinsically not rotating.

preprint2015arXiv

Assessing the Jeans Anisotropic Multi-Gaussian Expansion method with the Illustris simulation

We assess the effectiveness of the Jeans-Anisotropic-MGE (JAM) technique with a state-of-the-art cosmological hydrodynamic simulation, the Illustris project. We perform JAM modelling on 1413 simulated galaxies with stellar mass M^* > 10^{10}M_{sun}, and construct an axisymmetric dynamical model for each galaxy. Combined with a Markov Chain Monte Carlo (MCMC) simulation, we recover the projected root-mean-square velocity (V_rms) field of the stellar component, and investigate constraints on the stellar mass-to-light ratio, M^*/L, and the fraction of dark matter f_{DM} within 2.5 effective radii (R_e). We find that the enclosed total mass within 2.5 R_e is well constrained to within 10%. However, there is a degeneracy between the dark matter and stellar components with correspondingly larger individual errors. The 1 sigma scatter in the recovered M^*/L is 30-40% of the true value. The accuracy of the recovery of M^*/L depends on the triaxial shape of a galaxy. There is no significant bias for oblate galaxies, while for prolate galaxies the JAM-recovered stellar mass is on average 18% higher than the input values. We also find that higher image resolutions alleviate the dark matter and stellar mass degeneracy and yield systematically better parameter recovery.

preprint2015arXiv

P-MaNGA: Full spectral fitting and stellar population maps from prototype observations

MaNGA (Mapping Nearby Galaxies at Apache Point Observatory) is a 6-year SDSS-IV survey that will obtain resolved spectroscopy from 3600 $Å$ to 10300 $Å$ for a representative sample of over 10,000 nearby galaxies. In this paper, we derive spatially resolved stellar population properties and radial gradients by performing full spectral fitting of observed galaxy spectra from P-MaNGA, a prototype of the MaNGA instrument. These data include spectra for eighteen galaxies, covering a large range of morphological type. We derive age, metallicity, dust and stellar mass maps, and their radial gradients, using high spectral-resolution stellar population models, and assess the impact of varying the stellar library input to the models. We introduce a method to determine dust extinction which is able to give smooth stellar mass maps even in cases of high and spatially non-uniform dust attenuation. With the spectral fitting we produce detailed maps of stellar population properties which allow us to identify galactic features among this diverse sample such as spiral structure, smooth radial profiles with little azimuthal structure in spheroidal galaxies, and spatially distinct galaxy sub-components. In agreement with the literature, we find the gradients for galaxies identified as early-type to be on average flat in age, and negative (- 0.15 dex / R$_e$ ) in metallicity, whereas the gradients for late-type galaxies are on average negative in age (- 0.39 dex / R$_e$ ) and flat in metallicity. We demonstrate how different levels of data quality change the precision with which radial gradients can be measured. We show how this analysis, extended to the large numbers of MaNGA galaxies, will have the potential to shed light on galaxy structure and evolution.

preprint2015arXiv

P-MaNGA: Gradients in Recent Star Formation Histories as Diagnostics for Galaxy Growth and Death

We present an analysis of the data produced by the MaNGA prototype run (P-MaNGA), aiming to test how the radial gradients in recent star formation histories, as indicated by the 4000AA-break (D4000), Hdelta absorption (EW(Hd_A)) and Halpha emission (EW(Ha)) indices, can be useful for understanding disk growth and star formation cessation in local galaxies. We classify 12 galaxies observed on two P-MaNGA plates as either centrally quiescent (CQ) or centrally star-forming (CSF), according to whether D4000 measured in the central spaxel of each datacube exceeds 1.6. For each galaxy we generate both 2D maps and radial profiles of D4000, EW(Hd_A) and EW(Ha). We find that CSF galaxies generally show very weak or no radial variation in these diagnostics. In contrast, CQ galaxies present significant radial gradients, in the sense that D4000 decreases, while both EW(Hd_A) and EW(Ha) increase from the galactic center outward. The outer regions of the galaxies show greater scatter on diagrams relating the three parameters than their central parts. In particular, the clear separation between centrally-measured quiescent and star-forming galaxies in these diagnostic planes is largely filled in by the outer parts of galaxies whose global colors place them in the green valley, supporting the idea that the green valley represents a transition between blue-cloud and red-sequence phases, at least in our small sample. These results are consistent with a picture in which the cessation of star formation propagates from the center of a galaxy outwards as it moves to the red sequence.

preprint2015arXiv

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 Atlas3D 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 timescales. Using homogeneously derived effective radii and dynamically-determined galaxy masses, we present the distribution of stellar population parameters on the Mass Plane (M_JAM, Sigma_e, R_maj), 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 today's local Universe, approximately 50% 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_sun), which themselves formed 90% 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.

preprint2015arXiv

The HI Tully-Fisher Relation of Early-Type Galaxies

We study the HI K-band Tully-Fisher relation and the baryonic Tully-Fisher relation for a sample of 16 early-type galaxies, taken from the ATLAS3D sample, which all have very regular HI disks extending well beyond the optical body (> 5 R_eff). We use the kinematics of these disks to estimate the circular velocity at large radii for these galaxies. We find that the Tully-Fisher relation for our early-type galaxies is offset by about 0.5-0.7 magnitudes from the relation for spiral galaxies. The residuals with respect to the spiral Tully-Fisher relation correlate with estimates of the stellar mass-to-light ratio, suggesting that the offset between the relations is mainly driven by differences in stellar populations. We also observe a small offset between our Tully-Fisher relation with the relation derived for the ATLAS3D sample based on CO data representing the galaxies' inner regions (< 1 R_eff). This indicates that the circular velocities at large radii are systematically 10% lower than those near 0.5-1 R_eff, in line with recent determinations of the shape of the mass profile of early-type galaxies. The baryonic Tully-Fisher relation of our sample is distinctly tighter than the standard one, in particular when using mass-to-light ratios based on dynamical models of the stellar kinematics. We find that the early-type galaxies fall on the spiral baryonic Tully-Fisher relation if one assumes M/L_K = 0.54 M_sun/L_sun for the stellar populations of the spirals, a value similar to that found by recent studies of the dynamics of spiral galaxies. Such a mass-to-light ratio for spiral galaxies would imply that their disks are 60-70% of maximal. Our analysis increases the range of galaxy morphologies for which the baryonic Tully-Fisher relations holds, strengthening previous claims that it is a more fundamental scaling relation than the classical Tully-Fisher relation.

preprint2015arXiv

The Next Generation Virgo Cluster Survey. VI. The Kinematics of Ultra-compact Dwarfs and Globular Clusters in M87

The origin of ultra-compact dwarfs (UCDs)--objects larger and more massive than typical globular clusters (GCs), but more compact than typical dwarf galaxies--has been hotly debated in the 15 years since their discovery. Even whether UCDs should be considered galactic in origin, or simply the most extreme GCs, is not yet settled. We present the dynamical properties of 97 spectroscopically confirmed UCDs (rh >~10 pc) and 911 GCs associated with central cD galaxy of the Virgo cluster, M87. Our UCDs, of which 89% have M_star > ~2X10^6 M_sun and 92% are as blue as the classic blue GCs, nearly triple the sample of previous confirmed Virgo UCDs, providing by far the best opportunity for studying the global dynamics of a UCD system. We found that (1) UCDs have a surface number density profile that is shallower than that of the blue GCs in the inner ~ 70 kpc and as steep as that of the red GCs at larger radii; (2) UCDs exhibit a significantly stronger rotation than the GCs, and the blue GCs seem to have a velocity field that is more consistent with that of the surrounding dwarf ellipticals than with that of UCDs; (3) UCDs have a radially increasing orbital anisotropy profile, and are tangentially-biased at radii < ~ 40 kpc and radially-biased further out. In contrast, the blue GCs become more tangentially-biased at larger radii beyond ~ 40 kpc; (4) GCs with M_star > 2X10^6 M_sun have rotational properties indistinguishable from the less massive ones, suggesting that it is the size, instead of mass, that differentiates UCDs from GCs as kinematically distinct populations. We conclude that most UCDs in M87 are not consistent with being merely the most luminous and extended examples of otherwise normal GCs. The radially-biased orbital structure of UCDs at large radii is in general agreement with the "tidally threshed dwarf galaxy" scenario.

preprint2015arXiv

The Next Generation Virgo Cluster Survey. XII. Stellar Populations and Kinematics of Compact, Low-Mass Early-Type Galaxies from Gemini GMOS-IFU Spectroscopy

We present Gemini GMOS-IFU data of eight compact low-mass early-type galaxies (ETGs) in the Virgo cluster. We analyse their stellar kinematics, stellar population, and present two-dimensional maps of these properties covering the central 5"x 7" region. We find a large variety of kinematics: from non- to highly-rotating objects, often associated with underlying disky isophotes revealed by deep images from the Next Generation Virgo Cluster Survey. In half of our objects, we find a centrally-concentrated younger and more metal-rich stellar population. We analyze the specific stellar angular momentum through the lambdaR parameter and find six fast-rotators and two slow-rotators, one having a thin counter-rotating disk. We compare the local galaxy density and stellar populations of our objects with those of 39 more extended low-mass Virgo ETGs from the SMAKCED survey and 260 massive ($M>10^{10}$\Msun) ETGs from the A3D sample. The compact low-mass ETGs in our sample are located in high density regions, often close to a massive galaxy and have, on average, older and more metal-rich stellar populations than less compact low-mass galaxies. We find that the stellar population parameters follow lines of constant velocity dispersion in the mass-size plane, smoothly extending the comparable trends found for massive ETGs. Our study supports a scenario where low-mass compact ETGs have experienced long-lived interactions with their environment, including ram-pressure stripping and gravitational tidal forces, that may be responsible for their compact nature.

preprint2015arXiv

Unveiling the counter-rotating nature of the kinematically distinct core in NGC5813 with MUSE

MUSE observations of NGC5813 reveal a complex structure in the velocity dispersion map, previously hinted by SAURON observations. The structure is reminiscent of velocity dispersion maps of galaxies comprising two counter-rotating discs, and may explain the existence of the kinematically distinct core (KDC). Further evidence for two counter-rotating components comes from the analysis of the higher moments of the stellar line-of-sight velocity distributions and fitting MUSE spectra with two separate Gaussian line-of-sight velocity distributions. The emission-line kinematics show evidence of being linked to the present cooling flows and the buoyant cavities seen in X-rays. We detect ionised gas in a nuclear disc-like structure, oriented like the KDC, which is, however, not directly related to the KDC. We build an axisymmetric Schwarzschild dynamical model, which shows that the MUSE kinematics can be reproduced well with two counter-rotating orbit families, characterised by relatively low angular momentum components, but clearly separated in integral phase space and with radially varying contributions. The model indicates that the counter-rotating components in NGC5813 are not thin discs, but dynamically hot structures. Our findings give further evidence that KDCs in massive galaxies should not necessarily be considered as structurally or dynamically decoupled regions, but as the outcomes of the mixing of different orbital families, where the balance in the distribution of mass of the orbital families is crucial. We discuss the formation of the KDC in NGC5813 within the framework of gas accretion, binary mergers and formation of turbulent thick discs from cold streams at high redshift.

preprint2014arXiv

A kinematically distinct core and minor-axis rotation: the MUSE perspective on M87

We present evidence for the presence of a low-amplitude kinematically distinct component in the giant early-type galaxy M87, via datasets obtained with the SAURON and MUSE integral-field spectroscopic units. The MUSE velocity field reveals a strong twist of ~140 deg within the central 30 arcsec connecting outwards such a kinematically distinct core to a prolate-like rotation around the large-scale photometric major-axis of the galaxy. The existence of these kinematic features within the apparently round central regions of M87 implies a non-axisymmetric and complex shape for this galaxy, which could be further constrained using the presented kinematics. The associated orbital structure should be interpreted together with other tracers of the gravitational potential probed at larger scales (e.g., Globular Clusters, Ultra Compact Dwarfs, Planetary Nebulae): it would offer an insight in the assembly history of one of the brightest galaxies in the Virgo Cluster. These data also demonstrate the potential of the MUSE spectrograph to uncover low-amplitude spectral signatures.

preprint2014arXiv

Connection between dynamically derived IMF normalisation and stellar population parameters

We report on empirical trends between the dynamically determined stellar initial mass function (IMF) and stellar population properties for a complete, volume-limited sample of 260 early-type galaxies from the Atlas3D project. We study trends between our dynamically-derived IMF normalisation and absorption line strengths, and interpret these via single stellar population- (SSP-) equivalent ages, abundance ratios (measured as [alpha/Fe]), and total metallicity, [Z/H]. We find that old and alpha-enhanced galaxies tend to have on average heavier (Salpeter-like) mass normalisation of the IMF, but stellar population does not appear to be a good predictor of the IMF, with a large range of normalisation at a given population parameter. As a result, we find weak IMF-[alpha/Fe] and IMF-age correlations, and no significant IMF-[Z/H] correlation. The observed trends appear significantly weaker than those reported in studies that measure the IMF normalisation via low-mass star demographics inferred through stellar spectral analysis.

preprint2014arXiv

Identification of old tidal dwarfs near early-type galaxies from deep imaging and HI observations

It has recently been proposed that the dwarf spheroidal galaxies located in the Local Group disks of satellites (DoSs) may be tidal dwarf galaxies (TDGs) born in a major merger at least 5 Gyr ago. Whether TDGs can live that long is still poorly constrained by observations. As part of deep optical and HI surveys with the CFHT MegaCam camera and Westerbork Synthesis Radio Telescope made within the ATLAS3D project, and follow-up spectroscopic observations with the Gemini-North telescope, we have discovered old TDG candidates around several early-type galaxies. At least one of them has an oxygen abundance close to solar, as expected for a tidal origin. This confirmed pre-enriched object is located within the gigantic, but very low surface brightness, tidal tail that emanates from the elliptical galaxy, NGC 5557. An age of 4 Gyr estimated from its SED fitting makes it the oldest securely identified TDG ever found so far. We investigated the structural and gaseous properties of the TDG and of a companion located in the same collisional debris, and thus most likely of tidal origin as well. Despite several Gyr of evolution close to their parent galaxies, they kept a large gas reservoir. Their central surface brightness is low and their effective radius much larger than that of typical dwarf galaxies of the same mass. This possibly provides us with criteria to identify tidal objects which can be more easily checked than the traditional ones requiring deep spectroscopic observations. In view of the above, we discuss the survival time of TDGs and question the tidal origin of the DoSs.

preprint2014arXiv

Overview of the SDSS-IV MaNGA Survey: Mapping Nearby Galaxies at Apache Point Observatory

We present an overview of a new integral field spectroscopic survey called MaNGA (Mapping Nearby Galaxies at Apache Point Observatory), one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) that began on 2014 July 1. MaNGA will investigate the internal kinematic structure and composition of gas and stars in an unprecedented sample of 10,000 nearby galaxies. We summarize essential characteristics of the instrument and survey design in the context of MaNGA's key science goals and present prototype observations to demonstrate MaNGA's scientific potential. MaNGA employs dithered observations with 17 fiber-bundle integral field units that vary in diameter from 12" (19 fibers) to 32" (127 fibers). Two dual-channel spectrographs provide simultaneous wavelength coverage over 3600-10300 A at R~2000. With a typical integration time of 3 hr, MaNGA reaches a target r-band signal-to-noise ratio of 4-8 (per A, per 2" fiber) at 23 AB mag per sq. arcsec, which is typical for the outskirts of MaNGA galaxies. Targets are selected with stellar mass greater than 1e9 Msun using SDSS-I redshifts and i-band luminosity to achieve uniform radial coverage in terms of the effective radius, an approximately flat distribution in stellar mass, and a sample spanning a wide range of environments. Analysis of our prototype observations demonstrates MaNGA's ability to probe gas ionization, shed light on recent star formation and quenching, enable dynamical modeling, decompose constituent components, and map the composition of stellar populations. MaNGA's spatially resolved spectra will enable an unprecedented study of the astrophysics of nearby galaxies in the coming 6 yr.

preprint2014arXiv

The Atlas3D project - XXIV. The intrinsic shape distribution of early-type galaxies

We use the Atlas3D sample to perform a study of the intrinsic shapes of early-type galaxies, taking advantage of the available combined photometric and kinematic data. Based on our ellipticity measurements from the Sloan Digital Sky Survey Data Release 7, and additional imaging from the Isaac Newton Telescope, we first invert the shape distribution of fast and slow rotators under the assumption of axisymmetry. The so-obtained intrinsic shape distribution for the fast rotators can be described with a Gaussian with a mean flattening of q=0.25 and standard deviation sigma_q = 0.14, and an additional tail towards rounder shapes. The slow rotators are much rounder, and are well described with a Gaussian with mean q = 0.63 and sigma_q =0.09. We then checked that our results were consistent when applying a different and independent method to obtain intrinsic shape distributions, by fitting the observed ellipticity distributions directly using Gaussian parametrisations for the intrinsic axis ratios. Although both fast and slow rotators are identified as early-type galaxies in morphological studies, and in many previous shape studies are therefore grouped together, their shape distributions are significantly different, hinting at different formation scenarios. The intrinsic shape distribution of the fast rotators shows similarities with the spiral galaxy population. Including the observed kinematic misalignment in our intrinsic shape study shows that the fast rotators are predominantly axisymmetric, with only very little room for triaxiality. For the slow rotators though there are very strong indications that they are (mildly) triaxial.

preprint2014arXiv

The Atlas3D project - XXVI. HI discs in real and simulated fast and slow rotators

One quarter of all nearby early-type galaxies (ETGs) outside Virgo host a disc/ring of HI with size from a few to tens of kpc and mass up to ~1e+9 solar masses. Here we investigate whether this HI is related to the presence of a stellar disc within the host making use of the classification of ETGs in fast and slow rotators (FR/SR). We find a large diversity of HI masses and morphologies within both families. Surprisingly, SRs are detected as often, host as much HI and have a similar rate of HI discs/rings as FRs. Accretion of HI is therefore not always linked to the growth of an inner stellar disc. The weak relation between HI and stellar disc is confirmed by their frequent kinematical misalignment in FRs, including cases of polar and counterrotating gas. In SRs the HI is usually polar. This complex picture highlights a diversity of ETG formation histories which may be lost in the relative simplicity of their inner structure and emerges when studying their outer regions. We find that LCDM hydrodynamical simulations have difficulties reproducing the HI properties of ETGs. The gas discs formed in simulations are either too massive or too small depending on the star formation feedback implementation. Kinematical misalignments match the observations only qualitatively. The main point of conflict is that nearly all simulated FRs and a large fraction of all simulated SRs host corotating HI. This establishes the HI properties of ETGs as a novel challenge to simulations.

preprint2014arXiv

The ATLAS3D Project - XXVIII. Dynamically-driven star formation suppression in early-type galaxies

We present measurements of the star formation rate (SFR) in the early-type galaxies (ETGs) of the ATLAS3D sample, based on Wide-field Infrared Survey Explorer (WISE) 22um and Galaxy Evolution Explorer (GALEX) far-ultraviolet emission. We combine these with gas masses estimated from 12CO and HI data in order to investigate the star formation efficiency (SFE) in a larger sample of ETGs than previously available. We first recalibrate (based on WISE data) the relation between old stellar populations (traced at Ks-band) and 22um luminosity, allowing us to remove the contribution of 22um emission from circumstellar dust. We then go on to investigate the position of ETGs on the Kennicutt-Schmidt (KS) relation. Molecular gas-rich ETGs have comparable star formation surface densities to normal spiral galaxy centres, but they lie systematically offset from the KS relation, having lower star formation efficiencies by a factor of ~2.5 (in agreement with other authors). This effect is driven by galaxies where a substantial fraction of the molecular material is in the rising part of the rotation curve, and shear is high. We show here for the first time that although the number of stars formed per unit gas mass per unit time is lower in ETGs, it seems that the amount of stars formed per free-fall time is approximately constant. The scatter around this dynamical relation still correlates with galaxy properties such as the shape of the potential in the inner regions. This leads us to suggest that dynamical properties (such as shear or the global stability of the gas) may be important second parameters that regulate star formation and cause much of the scatter around star-formation relations.

preprint2014arXiv

The Atlas3D project -- XXIX. The new look of early-type galaxies and surrounding fields disclosed by extremely deep optical images

Galactic archeology 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, multi--band, images of nearby early-type galaxies. We present here a catalog of 92 galaxies from the Atlas3D sample, that are located in low to medium density environments. The observing strategy and data reduction pipeline, that 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 is 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 halos 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 is 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.

preprint2014arXiv

The interplay between a galactic bar and a supermassive black hole: nuclear fueling in a sub-parsec resolution galaxy simulation

We study the connection between the large-scale dynamics and the gas fueling toward a central black hole via the analysis of a Milky Way-like simulation at sub-parsec resolution. This allows us to follow a set of processes at various scales (e.g., the triggering of inward gas motion towards inner resonances via the large-scale bar, the connection to the central black hole via mini spirals) in a self-consistent manner. This simulation provides further insights on the role of shear for the inhibition of star formation within the bar in regions with significant amount of gas. We also witness the decoupling of the central gas and nuclear cluster from the large-scale disc, via interactions with the black hole. This break of symmetry in the mass distribution triggers the formation of gas clumps organised in a time-varying 250 pc ring-like structure, the black hole being offset by about 70 pc from its centre. Some clumps form stars, while most get disrupted or merge. Supernovae feedback further creates bubbles and filaments, some of the gas being expelled to 100 pc or higher above the galaxy plane. This helps remove angular momentum from the gas, which gets closer to the central dark mass. Part of the gas raining down is being accreted, forming a 10~pc polar disc-like structure around the black hole, leading to an episode of star formation. This gives rise to multiple stellar populations with significantly different angular momentum vectors, and may lead to a natural intermittence in the fueling of the black hole.

preprint2014arXiv

The Next Generation Virgo cluster Survey. V. modelling the dynamics of M87 with the Made-to-Measure method

We study the dynamics of the giant elliptical galaxy M87 from the central to the outermost regions with the made-to-measure (M2M) method. We use a new catalogue of 922 globular cluster line-of- sight velocities extending to a projected radius of 180 kpc (equivalent to 25 M87 effective radii), and SAURON integral field unit data within the central 2.4 kpc. 263 globular clusters, mainly located beyond 40 kpc, are newly observed by the Next Generation Virgo Survey (NGVS). For the M2M modelling, the gravitational potential is taken as a combination of a luminous matter potential with a constant stellar mass-to-light ratio and a dark matter potential modelled as a logarithmic potential. Our best dynamical model returns a stellar mass-to-light ratio in the I band of M/LI = 6.0(+ -0.3) M_sun/L_sun with a dark matter potential scale velocity of 591(+ -50) km/s and scale radius of 42(+ -10) kpc. We determine the total mass of M87 within 180 kpc to be (1.5 + - 0.2) 10^13 M_sun. The mass within 40 kpc is smaller than previous estimates determined using globular cluster kinematics that did not extend beyond 45 kpc. With our new globular cluster velocities at much larger radii, we see that globular clusters around 40 kpc show an anomalously large velocity dispersion which affected previous results. The mass we derived is in good agreement with that inferred from ROSAT X-ray observation out to 180 kpc. Within 30 kpc our mass is also consistent with that inferred from Chandra and XMM-Newton X-ray observations, while within 120 kpc it is about 20% smaller. The model velocity dispersion anisotropy beta parameter for the globular clusters in M87 is small, varying from -0.2 at the centre to 0.2 at 40 kpc, and gradually decreasing to zero at 120 kpc.

preprint2014arXiv

The Next Generation Virgo Cluster Survey. VIII. The Spatial Distribution of Globular Clusters in the Virgo Cluster

We report on a large-scale study of the distribution of globular clusters (GCs) throughout the Virgo cluster, based on photometry from the Next Generation Virgo Cluster Survey, a large imaging survey covering Virgo's primary subclusters to their virial radii. Using the g', (g'-i') color-magnitude diagram of unresolved and marginally-resolved sources, we constructed 2-D maps of the GC distribution. We present the clearest evidence to date showing the difference in concentration between red and blue GCs over the extent of the cluster, where the red (metal-rich) GCs are largely located around the massive early-type galaxies, whilst the blue (metal-poor) GCs have a more extended spatial distribution, with significant populations present beyond 83' (215 kpc) along the major axes of M49 and M87. The GC distribution around M87 and M49 shows remarkable agreement with the shape, ellipticity and boxiness of the diffuse light surrounding both galaxies. We find evidence for spatial enhancements of GCs surrounding M87 that may be indicative of recent interactions or an ongoing merger history. We compare the GC map to the locations of Virgo galaxies and the intracluster X-ray gas, and find good agreement between these baryonic structures. The Virgo cluster contains a total population of 67300$\pm$14400 GCs, of which 35% are located in M87 and M49 alone. We compute a cluster-wide specific frequency S_N,CL=$2.8\pm0.7$, including Virgo's diffuse light. The GC-to-baryonic mass fraction is e_b=$5.7\pm1.1\times10^{-4} $and the GC-to-total cluster mass formation efficiency is e_t=$2.9\pm0.5\times10^{-5}$, values slightly lower than, but consistent with, those derived for individual galactic halos. Our results show that the production of the complex structures in the unrelaxed Virgo cluster core (including the diffuse intracluster light) is an ongoing process.(abridged)

preprint2014arXiv

The role of turbulence in star formation laws and thresholds

The Schmidt-Kennicutt relation links the surface densities of gas to the star formation rate in galaxies. The physical origin of this relation, and in particular its break, i.e. the transition between an inefficient regime at low gas surface densities and a main regime at higher densities, remains debated. Here, we study the physical origin of the star formation relations and breaks in several low-redshift galaxies, from dwarf irregulars to massive spirals. We use numerical simulations representative of the Milky Way, the Large and the Small Magellanic Clouds with parsec up to subparsec resolution, and which reproduce the observed star formation relations and the relative variations of the star formation thresholds. We analyze the role of interstellar turbulence, gas cooling, and geometry in drawing these relations, at 100 pc scale. We suggest in particular that the existence of a break in the Schmidt- Kennicutt relation could be linked to the transition from subsonic to supersonic turbulence and is independent of self-shielding effects. This transition being connected to the gas thermal properties and thus to the metallicity, the break is shifted toward high surface densities in metal-poor galaxies, as observed in dwarf galaxies. Our results suggest that together with the collapse of clouds under self-gravity, turbulence (injected at galactic scale) can induce the compression of gas and regulate star formation.

preprint2013arXiv

Fueling AGN-I: How the Global Characteristics of the Central Kiloparsec of Seyferts differ from Quiescent Galaxies

In a matched sample of Seyfert and quiescent galaxies we simultaneously probe the stellar and molecular gas kinematics from 1 kpc down to 50 pc with the aim of identifying the dynamical processes dictating black hole accretion rates. This first paper compares the global characteristics of a sample of ten galaxies. We find several differences within a radius of 500 pc that are correlated with Active Galactic Nucleus (AGN) activity. The Seyferts have: (1) a more centrally concentrated nuclear stellar surface brightness with lower stellar luminosities beyond a radius of 100 pc, (2) a lower stellar velocity dispersion within a radius of 200 pc, (3) elevated H2 luminosity out to a radius of at least 250 pc, and (4) more centrally concentrated H2 surface brightness profiles. These observed differences can be interpreted as evidence for Seyfert galaxies having a dynamically cold (in comparison to the bulge) nuclear structure composed of a significant gas reservoir and a relatively young stellar population. This structure is undetected (and possibly does not exist) in quiescent galaxies. The presence of such a nuclear structure in Seyfert galaxies provides evidence for inflow of the surrounding interstellar medium since the nuclear stellar population requires a supply of gas from which to form. The fueling of a Seyfert AGN is thus associated with the formation of a dynamically cold component of gas and stars on scales of hundreds of parsecs.

preprint2013arXiv

Is the black hole in NGC1277 really over-massive?

A claim has been made by van den Bosch et al. (2012) that NGC1277 hosts an over-massive BH with a mass larger than half its spheroid mass. We revisit this claim by examining the predictions from dynamical realisations based on new MGE models of NGC1277. We present realisations which fit the observed photometry. M/L is fixed following scaling relations which predict a Salpeter-like IMF. A model without a BH provides a surprisingly good fit of the observed kinematics outside the unresolved central region, but not, as expected, of the central dispersion and h4 values. A model with a MBH of 5 10^9 Msun allows to fit the dispersion profile, consistently with models of the same mass and M/L in vdB+12. It departs from the central h4 values by only about twice the given uncertainty. A slightly varying M/L or the addition of high velocity stars would further lower the need for a very massive BH. These results do not rule out the presence of an over-massive BH at the centre of NGC1277. However, they lead us to advocate the use of 3-sigma confidence intervals for derived MBH as better, more conservative, guidelines for such studies. We caution for the use of ill-defined spheroidal components as an input for scaling relations, and emphasise the fact that a MBH in the range 2-5 10^9 Msun would represent less than 5% of the spheroid mass of our models. This would make the BH in NGC1277 consistent or just twice as large as what a recent version of the MBH-sigma relation predicts. We examine the impact of the presence of a bar by running simulations from the same MGE model but with extreme anisotropies. An inner small bar forms, and an end-on view gets closer to fitting the central dispersion profile without the need for a central BH, while adding a black hole of 2.5 10^9 Msun, in line with the prediction from scaling relations, allows to fit the dispersion peak and h3 profiles. EDITED and ABRIDGED.

preprint2013arXiv

NGC1266 as a Local Candidate for Rapid Cessation of Star Formation

We present new Spectrographic Areal Unit for Research on Optical Nebulae (SAURON) integral-field spectroscopy and Swift Ultraviolet Optical Telescope (UVOT) observations of molecular outflow host galaxy NGC 1266 that indicate NGC 1266 has experienced a rapid cessation of star formation. Both the SAURON maps of stellar population age and the Swift UVOT observations demonstrate the presence of young ($< 1$ Gyr) stellar populations within the central 1 kpc, while existing Combined Array for Research in Millimeter-wave Astronomy (CARMA) CO(1--0) maps indicate that the sites of current star formation are constrained to the inner few hundred parsecs of the galaxy only. The optical spectrum of NGC 1266 from Moustakas & Kennicutt (2006) reveal a characteristic post-starburst (K+A) stellar population and Davis et al. (2012) confirm that ionized gas emission in the system originate from a shock. Galaxies with K+A spectra and shock-like ionized gas line ratios may comprise an important, overlooked segment of the post-starburst population, containing exactly those objects in which the AGN is actively expelling the star-forming material. While AGN activity is not the likely driver of the post-starburst event that occurred 500 Myr ago, the faint spiral structure seen in the Hubble Space Telescope (HST) Wide-field Camera 3 (WFC3) Y-, J- and H-band imaging seems to point to the possibility of gravitational torques being the culprit. If the molecular gas were driven into the center at the same time as the larger scale galaxy disk underwent quenching, the AGN might be able to sustain the presence of molecular gas for $\gtrsim 1$ Gyr by cyclically injecting turbulent kinetic energy into the dense molecular gas via a radio jet, inhibiting star formation.

preprint2013arXiv

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>6*10^9 Msun) Atlas3D sample of 260 early-type galaxies. We construct detailed axisymmetric dynamical models (JAM), 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. We derive accurate total M/L and dark matter fractions f_DM, within a sphere of radius r=Re. We also measure the stellar M/L and derive a median dark matter fraction f_DM=13%. We find that the thin two-dimensional subset spanned by galaxies in the (M_JAM,sigma_e,R_e) coordinates system, which we call the Mass Plane (MP) has an observed rms scatter of 19% and an intrinsic one of 11%. The MP satisfies the scalar virial relation M_JAM sigma_e^2 R_e within our tight errors. However, the details of how both Re and sigma_e are determined are critical in defining the precise deviation from the virial exponents. We revisit the (M/L)-sigma_e relation, which describes most of the deviations between the MP and the FP. The best-fitting relation is (M/L) sigma_e^0.72 (r-band). It provides an upper limit to any systematic increase of the IMF mass normalization with sigma_e. We study of the link between sigma_e and the galaxies circular velocity V_circ within 1Re (where stars dominate) and find the relation max(V_circ)~1.76*sigma_e, which has an observed scatter of 7%. The accurate parameters described in this paper are used in the companion Paper XX of this series to explore the variation of global galaxy properties, including the IMF, on the projections of the MP. [Abridged]

preprint2013arXiv

The Atlas3D project - XX. Mass-size and mass-sigma distributions of early-type galaxies: bulge fraction drives kinematics, mass-to-light ratio, molecular gas fraction and stellar initial mass function

We study the (M,sigma) and (M,Re) projections of the thin Mass Plane (MP) (M,sigma,Re) which describes the distribution of the galaxy population. The distribution of galaxy properties on the MP is characterized by: (i) a zone of exclusion described by two power-laws joined by a break at M 3e10 Msun. This results in a break in the mean M-sigma relation with Msigma^2.3 and M sigma^4.7 at small and large sigma respectively; (ii) a mass M 2e11 Msun which separates a population dominated by flat fast rotator with disks and spiral galaxies at lower masses, from one dominated by quite round slow rotators at larger masses; (iii) below that mass the distribution of ETGs properties tends to be constant along lines of constant sigma; (iv) it forms a parallel sequence with the distribution of spiral galaxies; (v) at even lower masses, the distribution of fast rotator ETGs and late spirals naturally extends to that of dwarf ETGs (Sph) and dwarf irregulars (Im) respectively. We show via dynamical models that sigma traces the bulge fraction, which drives the observed trends in M/L, Hbeta, colour, IMF and molecular gas fraction. We interpret this as due to a combination of two main effects: (i) an increase of the bulge fraction which increases sigma and greatly enhances the likelihood for a galaxy to have its star formation quenched, and (ii) dry merging, increasing galaxy mass along lines of nearly constant sigma, while leaving the population unchanged. [Abriged]

preprint2013arXiv

The ATLAS3D Project -- XXIII. Angular momentum and nuclear surface brightness profiles

[Abridged and Edited] We investigate nuclear light profiles in 135 ATLAS3D galaxies for which the Hubble Space Telescope (HST) imaging is available and compare them to the large scale kinematics obtained with the SAURON integral-field spectrograph. Specific angular momentum, lambda_R, correlates with the shape of nuclear light profiles, where cores are typically found in slow rotators and core-less galaxies are fast rotators. Cores are found only in massive galaxies and only in systems with the stellar mass M>8x10^10 Msun. Based on our sample, we, however, see no evidence for a bimodal distribution of nuclear slopes. The best predictor for finding a core is based on the stellar velocity dispersion within an effective radius, sigma_e, and specific angular momentum, where cores are found for lambda_R<0.25 and sigma_e>160 km/s. We estimate that only about 10% of nearby early-type galaxies contain cores. Furthermore, we show that there is a genuine population of fast rotators with cores. We also show that core fast rotators are morphologically, kinematically and dynamically different from core slow rotators. The cores of fast rotators could harbour black holes of similar masses to those in core slow rotators, but typically more massive than those found in core-less fast rotators. Core-less galaxies, and especially core-less fast rotators, are under-luminous in the diffuse X-ray emission, but the presence of a core does not imply high X-ray luminosities. We postulate a possible population of core-less galaxies among slow rotators, which can not be explained as face-on discs, but comprise a genuine sub-population of slow rotators. These galaxies are typically less massive and flatter than core slow rotators, and show evidence for dynamical cold structures and exponential photometric components. We discuss possible processes for formation of cores and their subsequent preservation.

preprint2013arXiv

The Next Generation Virgo Cluster Survey - Infrared (NGVS-IR): I. A new Near-UV/Optical/Near-IR Globular Cluster selection tool

The NGVS-IR project (Next Generation Virgo Survey - Infrared) is a contiguous near-infrared imaging survey of the Virgo cluster of galaxies. It complements the optical wide-field survey of Virgo (NGVS). The current state of NGVS-IR consists of Ks-band imaging of 4 deg^2 centered on M87, and J and Ks-band imaging of 16 deg^2 covering the region between M49 and M87. In this paper, we present the observations of the central 4 deg^2 centered on Virgo's core region. The data were acquired with WIRCam on the Canada-France-Hawaii Telescope and the total integration time was 41 hours distributed in 34 contiguous tiles. A survey-specific strategy was designed to account for extended galaxies while still measuring accurate sky brightness within the survey area. The average 5σlimiting magnitude is Ks=24.4 AB mag and the 50% completeness limit is Ks=23.75 AB mag for point source detections, when using only images with better than 0.7" seeing (median seeing 0.54"). Star clusters are marginally resolved in these image stacks, and Virgo galaxies with μ_Ks=24.4 AB mag arcsec^-2 are detected. Combining the Ks data with optical and ultraviolet data, we build the uiK color-color diagram which allows a very clean color-based selection of globular clusters in Virgo. This diagnostic plot will provide reliable globular cluster candidates for spectroscopic follow-up campaigns needed to continue the exploration of Virgo's photometric and kinematic sub-structures, and will help the design of future searches for globular clusters in extragalactic systems. Equipped with this powerful new tool, future NGVS-IR investigations based on the uiK diagram will address the mapping and analysis of extended structures and compact stellar systems in and around Virgo galaxies.

preprint2013arXiv

The Next Generation Virgo Cluster Survey. IV. NGC 4216: A Bombarded Spiral in the Virgo Cluster

We present an investigation into the origins of a series of interlaced narrow filamentary stellar structures, loops and plumes in the vicinity of the Virgo Cluster, edge-on spiral galaxy, NGC 4216 that were previously identified by the Blackbird Telescope. Using the deeper, higher-resolution and precisely calibrated optical CFHT/MegaCam images obtained as part of the Next Generation Virgo Cluster Survey (NGVS), we confirm the previously identified features and identify a few additional structures. The NGVS data allowed us to make a physical study of these low-surface brightness features and investigate their origin. The likely progenitors of the structures were identified as either already catalogued VCC dwarfs or newly discovered satellites caught in the act of being destroyed. They have the same g-i color index and likely contain similar stellar populations. The alignment of three dwarfs along an apparently single stream is intriguing, and we cannot totally exclude that these are second-generation dwarf galaxies being born inside the filament from the debris of an original dwarf. The observed complex structures, including in particular a stream apparently emanating from a satellite of a satellite, point to a high rate of ongoing dwarf destruction/accretion in the region of the Virgo Cluster where NGC 4216 is located. We discuss the age of the interactions and whether they occurred in a group that is just falling into the cluster and shows signs of so-called "pre-processing" before it gets affected by the cluster environment, or in a group which already ventured towards the central regions of Virgo Cluster.

preprint2012arXiv

A 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. It depends on the stellar initial mass function (IMF) describing the distribution of stellar masses when the population formed. Consequently knowledge of the IMF is critical to virtually every aspect of galaxy evolution. More than half a century after the first IMF determination, no consensus has emerged on whether it is universal in different galaxies. Previous studies indicated that the IMF and the dark matter fraction in galaxy centres cannot be both universal, but they could not break the degeneracy between the two effects. Only recently indications were found that massive elliptical galaxies may not have the same IMF as our Milky Way. Here we report unambiguous evidence for a strong systematic variation of the IMF in early-type galaxies as a function of their stellar mass-to-light ratio, producing differences up to a factor of three in mass. This was inferred from detailed dynamical models of the two-dimensional stellar kinematics for the large Atlas3D representative sample of nearby early-type galaxies spanning two orders of magnitude in stellar mass. Our finding indicates that the IMF depends intimately on a galaxy's formation history.

preprint2012arXiv

Discovery of a giant HI tail in the galaxy group HCG 44

We report the discovery of a giant HI tail in the intra-group medium of HCG 44 as part of the Atlas3D survey. The tail is ~300 kpc long in projection and contains ~5x10^8 M_sun of HI. We detect no diffuse stellar light at the location of the tail down to ~28.5 mag/arcsec^2 in g band. We speculate that the tail might have formed as gas was stripped from the outer regions of NGC 3187 (a member of HCG 44) by the group tidal field. In this case, a simple model indicates that about 1/3 of the galaxy's HI was stripped during a time interval of <1 Gyr. Alternatively, the tail may be the remnant of an interaction between HCG 44 and NGC 3162, a spiral galaxy now ~650 kpc away from the group. Regardless of the precise formation mechanism, the detected HI tail shows for the first time direct evidence of gas stripping in HCG 44. It also highlights that deep HI observations over a large field are needed to gather a complete census of this kind of events in the local Universe.

preprint2012arXiv

Gemini GMOS and WHT SAURON integral-field spectrograph observations of the AGN driven outflow in NGC 1266

We use the SAURON and GMOS integral field spectrographs to observe the active galactic nucleus (AGN) powered outflow in NGC 1266. This unusual galaxy is relatively nearby (D=30 Mpc), allowing us to investigate the process of AGN feedback in action. We present maps of the kinematics and line strengths of the ionised gas emission lines Halpha, Hbeta, [OIII], [OI], [NII] and [SII], and report on the detection of Sodium D absorption. We use these tracers to explore the structure of the source, derive the ionised and atomic gas kinematics and investigate the gas excitation and physical conditions. NGC 1266 contains two ionised gas components along most lines of sight, tracing the ongoing outflow and a component closer to the galaxy systemic, the origin of which is unclear. This gas appears to be disturbed by a nascent AGN jet. We confirm that the outflow in NGC 1266 is truly multiphase, containing radio plasma, atomic, molecular and ionised gas and X-ray emitting plasma. The outflow has velocities up to \pm900 km/s away from the systemic velocity, and is very likely to be removing significant amounts of cold gas from the galaxy. The LINER-like line-emission in NGC 1266 is extended, and likely arises from fast shocks caused by the interaction of the radio jet with the ISM. These shocks have velocities of up to 800 km/s, which match well with the observed velocity of the outflow. Sodium D equivalent width profiles are used to set constraints on the size and orientation of the outflow. The ionised gas morphology correlates with the nascent radio jets observed in 1.4 GHz and 5 GHz continuum emission, supporting the suggestion that an AGN jet is providing the energy required to drive the outflow.

preprint2012arXiv

How to recover both velocity components in discs of barred galaxies with integral-field spectroscopy

We present a new method that derives both velocity components in the equatorial plane of a barred stellar disc from the observed line-of-sight velocity, assuming geometry of a thin disc. The method can be applied to large departures from circular motion, and does not require multipole decomposition. It is based on assumptions that the bar is close to steady-state (i.e. does not evolve fast), and that both morphology and kinematics are symmetrical with respect to the major axis of the bar. We derive the equations used in the method, and analyze the effect of observational errors on the inferred velocity fields. We show that this method produces meaningful results via a simple toy model. We also apply the method on integral-field data of NGC 936, for which we recover both velocity components in the disc. Knowing both velocity components in the disc, i.e. the non-observable transverse velocity in addition to the line-of-sight velocity, puts additional constraints on dynamical models and allows for new ways of determining parameters that are crucial in characterizing galaxies.

preprint2012arXiv

Simulations of binary galaxy mergers and the link with Fast Rotators, Slow Rotators, and Kinematically Distinct Cores

We study the formation of early-type galaxies (ETGs) through mergers with a sample of 70 high-resolution numerical simulations of binary mergers of disc galaxies. These simulations encompass various mass ratios, initial conditions and orbital parameters. We find that binary mergers of disc galaxies with mass ratios of 3:1 and 6:1 are nearly always classified as Fast Rotators according to the Atlas3D criterion: they preserve the structure of the input fast rotating spiral progenitors. Major disc mergers (mass ratios of 2:1 and 1:1) lead to both Fast and Slow Rotators. Most of the Slow Rotators hold a stellar Kinematically Distinct Core (KDC) in their 1-3 central kilo-parsec: these KDCs are built from the stellar components of the progenitors. The mass ratio of the progenitors is a fundamental parameter for the formation of Slow Rotators in binary mergers, but it also requires a retrograde spin for the progenitor galaxies with respect to the orbital angular momentum. The importance of the initial spin of the progenitors is also investigated in the library of galaxy mergers of the GalMer project.

preprint2012arXiv

The Atlas3D project - XII. Recovery of the mass-to-light ratio of simulated early-type barred galaxies with axisymmetric dynamical models

We investigate the accuracy in the recovery of the stellar dynamics of barred galaxies when using axisymmetric dynamical models. We do this by trying to recover the mass-to-light ratio (M/L) and the anisotropy of realistic galaxy simulations using the Jeans Anisotropic Multi-Gaussian Expansion (JAM) method. However, given that the biases we find are mostly due to an application of an axisymmetric modeling algorithm to a non-axisymmetric system and in particular to inaccuracies in the de-projected mass model, our results are relevant for general axisymmetric modelling methods. We run N-body collisionless simulations to build a library with various luminosity distribution, constructed to mimic real individual galaxies, with realistic anisotropy. The final result of our evolved library of simulations contains both barred and unbarred galaxies. The JAM method assumes an axisymmetric mass distribution, and we adopt a spatially constant M/L and anisotropy beta_z=1-sigma_z^2/sigma_R^2 distributions. The models are fitted to two-dimensional maps of the second velocity moments V_rms=sqrt(V^2+sigma^2) of the simulations for various viewing angles (position angle of the bar and inclination of the galaxy). We find that the inclination is generally well recovered by the JAM models, for both barred and unbarred simulations. For unbarred simulations the M/L is also accurately recovered, with negligible median bias and with a maximum one of just Delta(M/L)<1.5% when the galaxy is not too close to face on. At very low inclinations (i<30 deg) the M/L can be significantly overestimated (9% in our tests, but errors can be larger for very face-on views). For barred simulations the M/L is on average (when PA=45 deg) essentially unbiased, but we measure an over/under estimation of up to Delta(M/L)=15% in our tests. The sign of the M/L bias depends on the position angle of the bar as expected. [Abridged]

preprint2012arXiv

The Atlas3D project - XIV. The extent and kinematics of molecular gas in early-type galaxies

We use interferometric CO observations to compare 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 scale-lengths. Virgo cluster ETGs have less extended molecular gas reservoirs than field counterparts. 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. We suggest that ram pressure stripping and/or the presence of hot gas has compressed/truncated the gas in these systems. 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 ionised gas velocities. We confirm that the molecular gas reaches beyond the turnover of the circular velocity curve in ~70% of our CO-rich Atlas3D ETGs. In general we find that in most galaxies the molecular gas is relaxed and dynamically cold. Molecular gas is a better direct tracer of the circular velocity than the ionised gas, justifying its use as a kinematic tracer for Tully-Fisher and similar analyses. (abridged)

preprint2012arXiv

The ATLAS3D project - XVII. Linking photometric and kinematic signatures of stellar discs in early-type galaxies

[Abridged] We analyse the morphological structures in galaxies of the ATLAS3D sample by fitting a single Sersic profile and decomposing all non-barred objects (180 of 260 objects) in two components parameterised by an exponential and a general Sersic function. The aim of this analysis is to look for signatures of discs in light distributions of nearby early-type galaxies and compare them to kinematic properties. Using Sersic index from single component fits for a distinction between slow and fast rotators, or even late- and early-type galaxies, is not recommended. Assuming that objects with n>3 are slow rotators (or ellipticals), there is only a 22 per cent probability to correctly classify objects as slow rotators (or 37 per cent of previously classified as ellipticals). We show that exponential sub-components, as well as light profiles fitted with only a single component of a low Sersic index, can be linked with the kinematic evidence for discs in early-type galaxies. The median disk-to-total light ratio for fast and slow rotators is 0.41 and 0.0, respectively. Similarly, the median Sersic indices of the bulge (general Sersic component) are 1.7 and 4.8 for fast and slow rotators, respectively. Overall, discs or disc-like structures, are present in 83 per cent of early-type galaxies which do not have bars, and they show a full range of disk-to-total light ratios. Discs in early-type galaxies contribute with about 40 per cent to the total mass of the analysed (non-barred) objects. The decomposition into discs and bulges can be used as a rough approximation for the separation between fast and slow rotators, but it is not a substitute, as there is only a 59 per cent probability to correctly recognise slow rotators. Kinematics (i.e. projected angular momentum) remains the best approach to mitigate the influence of the inclination effects.

preprint2012arXiv

The ATLAS3D project - XVIII. CARMA CO imaging survey of early-type galaxies

We present the Combined Array for Research in Millimeter Astronomy (CARMA) ATLAS3D molecular gas imaging survey, a systematic study of the distribution and kinematics of molecular gas in CO-rich early-type galaxies. Our full sample of 40 galaxies (30 newly mapped and 10 taken from the literature) is complete to a 12CO(1-0) integrated flux of 18.5 Jy km/s, and it represents the largest, best-studied sample of its type to date. A comparison of the CO distribution of each galaxy to the g-r color image (representing dust) shows that the molecular gas and dust distributions are in good agreement and trace the same underlying interstellar medium. The galaxies exhibit a variety of CO morphologies, including discs (50%), rings (15%), bars+rings (10%), spiral arms (5%), and mildly (12.5%) and strongly (7.5%) disrupted morphologies. There appear to be weak trends between galaxy mass and CO morphology, whereby the most massive galaxies in the sample tend to have molecular gas in a disc morphology. We derive a lower limit to the total accreted molecular gas mass across the sample of 2.48x10^10 Msuns, or approximately 8.3x10^8 Msuns per minor merger within the sample, consistent with minor merger stellar mass ratios.

preprint2012arXiv

The ATLAS3D Project - XXI. Correlations between gradients of local escape velocity and stellar populations in early-type galaxies

We explore the connection between the local escape velocity, V_esc, and the stellar population properties in the ATLAS3D survey, a complete, volume-limited sample of nearby early-type galaxies. We make use of ugriz photometry to construct Multi-Gaussian Expansion models of the surface brightnesses of our galaxies. We are able to fit the full range of surface brightness profiles found in our sample, and in addition we reproduce the results of state-of-the-art photometry in the literature with residuals of 0.04 mags. We utilise these photometric models and SAURON integral-field spectroscopy, combined with Jeans dynamical modelling, to determine the local V_esc derived from the surface brightness. We find that the local V_esc is tightly correlated with the Mgb and Fe5015 linestrengths and optical colours, and anti-correlated with the Hbeta linestrength. In the case of the Mgb and Colour - V_esc relations we find that the relation within individual galaxies follows the global relation between different galaxies. We intentionally ignored any uncertain contribution due to dark matter since we are seeking an empirical description of stellar population gradients in early-type galaxies that is ideal for quantitative comparison with model predictions. We identify a population of galaxies that occur only at low V_esc that exhibit negative gradients in the Mgb - and Colour -V_esc relations. These galaxies typically have young central stellar populations and contain significant amounts of molecular gas and dust. Combining these results with N-body simulations of binary mergers we use the Mgb-V_esc relation to constrain the possible number of dry mergers experienced by the local early-type galaxy population - a typical massive ETG can have experienced only ~1.5 major mergers before becoming a significant outlier in the Mgb-V_esc relation. [Abridged]

preprint2012arXiv

The Planetary Nebulae Population in the Nuclear Regions of M31: the SAURON view

Following a first study of the central regions of M32 that illustrated the power of integral-field spectroscopy (IFS) in detecting and measuring the [O III]λ5007 emission of PNe against a strong stellar background, we turn to the very nuclear PN population of M31, within 80 pc of its centre. We show that PNe can also be found in the presence of emission from diffuse gas and further illustrate the excellent sensitivity of IFS in detecting extragalactic PNe through a comparison with narrowband images obtained with the Hubble Space Telescope. We find the nuclear PNe population of M31 is only marginally consistent with the generally adopted form of the PNe luminosity function (PNLF). In particular, this is due to a lack of PNe with absolute magnitude M5007 brighter than -3, which would only result from a rather unfortunate draw from such a model PNLF. We suggest that the observed lack of bright PNe in the nuclear regions of M31 is due to a horizontal-branch population that is more tilted toward less massive and hotter He-burning stars, so that its progeny consists mostly of UV-bright stars that fail to climb back up the asymptotic giant branch (AGB) and only of few, if any, bright PNe powered by central post-AGB stars. These results are also consistent with recent reports on a dearth of bright post-AGB stars towards the nucleus of M31, and lend further support to the idea that the metallicity of a stellar population has an impact on the way the horizontal branch is populated and to the loose anticorrelation between the strength of the UV-upturn and the specific number of PNe that is observed in early-type galaxies. Finally, our investigation also serves to stress the importance of considering the same spatial scales when comparing the PNe population of galaxies with the properties of their stellar populations.

preprint2012arXiv

The SAURON project-XXI. The spatially-resolved UV-line strength relations of early-type galaxies

The unexpected rising flux of early-type galaxies at decreasing ultraviolet (UV) wavelengths is a long-standing mystery. One important observational constraint is the correlation between UV-optical colours and Mg2 line strengths found by Burstein et al. (1988). The simplest interpretation of this phenomenon is that the UV strength is related to the Mg line strength. Under this assumption, we expect galaxies with larger Mg gradients to have larger UV colour gradients. By combining UV imaging from GALEX, optical imaging from MDM and SAURON integral-field spectroscopy, we investigate the spatially-resolved relationships between UV colours and stellar population properties of 34 early-type galaxies from the SAURON survey sample. We find that galaxies with old stellar populations show tight correlations between the FUV colours (FUV-V and FUV-NUV) and the Mgb index, Hβ index and metallicity [Z/H]. The equivalent correlations for the Fe5015 index, α-enhancement [α/Fe] and age are present but weaker. We have also derived logarithmic internal radial colour, measured line strength and derived stellar population gradients for each galaxy and again found a strong dependence of the FUV-V and FUV-NUV colour gradients on both the Mg b line strength and the metallicity gradients for galaxies with old stellar populations. In particular, global gradients of Mg b and [Z/H] with respect to the UV colour across galaxies are consistent with their local gradients within galaxies, suggesting that the global correlations also hold locally. From a simple model based on multi-band colour fits of UV upturn and UV-weak galaxies, we have identified a plausible range of parameters that reproduces the observed radial colour profiles. In these models, the centers of elliptical galaxies, where the UV flux is strong, are enhanced in metals by roughly 60% compared to UV-weak regions.

preprint2011arXiv

A Magellan-IMACS-IFU Search for Dynamical Drivers of Nuclear Activity. I. Reduction Pipeline and Galaxy Catalog

Using the Inamori Magellan Areal Camera and Spectrograph (IMACS) integral-field unit (IFU) on the 6.5m Magellan telescope, we have designed the first statistically significant investigation of the two-dimensional distribution and kinematics of ionized gas and stars in the central kiloparsec regions of a well-matched sample of Seyfert and inactive control galaxies selected from the Sloan Digital Sky Survey. The goals of the project are to use the fine spatial sampling (0.2 arcsec/pixel) and large wavelength coverage (4000-7000A) of the IMACS-IFU to search for dynamical triggers of nuclear activity in the central region where active galactic nucleus (AGN) activity and dynamical timescales become comparable, to identify and assess the impact of AGN-driven outflows on the host galaxy and to provide a definitive sample of local galaxy kinematics for comparison with future three-dimensional kinematic studies of high-redshift systems. In this paper, we provide the first detailed description of the procedure to reduce and calibrate data from the IMACS-IFU in `long mode' to obtain two-dimensional maps of the distribution and kinematics of ionized gas and stars. The sample selection criteria are presented, observing strategy described and resulting maps of the sample galaxies presented along with a description of the observed properties of each galaxy and the overall observed properties of the sample.

preprint2011arXiv

Discovery of an AGN-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 and requires a nested Gaussian fit to be properly described. Interferometric observations reveal a massive, centrally- concentrated molecular component with a mass of 1.1x10^9 Msuns and a molecular outflow with a molecular mass of 2.4x10^7 Msuns. The molecular gas close to the systemic velocity consists of a rotating, compact nucleus with a mass of about 4.1x108 Msuns within a radius of approximately 60 pc. This compact molecular nucleus has a surface density of \approx 2.7 \times 10^4 Msuns/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 Msuns/yr leads to a depletion timescale of about 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 (AGN). 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.

preprint2011arXiv

The Atlas-3D project - IX. The merger origin of a fast and a slow rotating Early-Type Galaxy revealed with deep optical imaging: first results

The mass assembly of galaxies leaves imprints in their outskirts, such as shells and tidal tails. The frequency and properties of such fine structures depend on the main acting mechanisms - secular evolution, minor or major mergers - and on the age of the last substantial accretion event. We use this to constrain the mass assembly history of two apparently relaxed nearby Early-Type Galaxies (ETGs) selected from the Atlas-3D sample, NGC 680 and NGC 5557. Our ultra deep optical images obtained with MegaCam on the Canada-France-Hawaii Telescope reach 29 mag/arcsec^2 in the g-band. They reveal very low-surface brightness (LSB) filamentary structures around these ellipticals. Among them, a gigantic 160 kpc long tail East of NGC 5557 hosts gas-rich star-forming objects. NGC 680 exhibits two major diffuse plumes apparently connected to extended HI tails, as well as a series of arcs and shells. Comparing the outer stellar and gaseous morphology of the two ellipticals with that predicted from models of colliding galaxies, we argue that the LSB features are tidal debris, the star-forming objects near NGC 5557, long lived Tidal Dwarf Galaxies and that each of the two ETGs was assembled during a relatively recent, major wet merger, which likely occurred at a redshift below z = 0.5. The inner kinematics of NGC 680 is typical for fast rotators which make the bulk of nearby ETGs in the Atlas-3D sample. NGC 5557 belongs to the poorly populated class of massive, round, slow rotators that are predicted by semi-analytic models and cosmological simulations to be the end-product of a complex mass accretion history, involving ancient major mergers and more recent minor mergers. Our observations suggest that under specific circumstances a single binary merger may dominate the formation history of such objects and thus that at least some massive ETGs may form at relatively low redshift (abridged).

preprint2011arXiv

The ATLAS3D project - II. Morphologies, kinemetric features and alignment between photometric and kinematic axes of early-type galaxies

[Abridged] We use the ATLAS3D sample of 260 early-type galaxies (ETGs) to study the apparent kinematic misalignment angle, Psi, defined as the angle between the photometric and kinematic major axis. We find that 71% of nearby ETGs are strictly aligned systems (Psi > 5 deg), an additional 14% have 5 < Psi < 10 deg and 90% of galaxies have Psi < 15 deg. Taking into account measurement uncertainties, 90% of galaxies can be considered aligned to better than 5 deg, suggesting that only a small fraction of early-type galaxies (~10%) are not consistent with axisymmetry within the projected half-light radius. We identify morphological features such as bars and rings (30%), dust structures (16%), blue nuclear colours (6%) and evidence of interactions (8%). We use kinemetry to analyse the mean velocity maps and separate galaxies in two broad types of regular and non-regular rotators. We find 82% of regular rotators and 17% non-regular rotators, with 2 galaxies that we were not able to classify due to data quality. The non-regular rotators are typically found in dense regions and are massive. The majority of galaxies do not have any specific kinemetric features, but we highlight here the frequency of the kinematically distinct cores (7%) and the aligned double peaks in the velocity dispersion maps (4%). Most of the galaxies that are misaligned have complex kinematics and are non-regular rotators. While the trends are weak, there is a tendency that large values of Psi are found in galaxies at intermediate environmental densities and among the most massive galaxies in the sample. We suggest that the most common formation mechanism for ETGs preserves the axisymmetry of the disk progenitors. Less commonly, the formation process results in a triaxial galaxy with much lower net angular momentum.

preprint2011arXiv

The ATLAS3D project - III. A census of the stellar angular momentum within the effective radius of early-type galaxies: unveiling the distribution of Fast and Slow Rotators

We provide a census of the apparent stellar angular momentum within 1 Re of a volume-limited sample of 260 early-type galaxies (ETGs) in the nearby Universe, using integral-field spectroscopy obtained in the course of the ATLAS3D project. We exploit the LambdaR parameter to characterise the existence of two families of ETGs: Slow Rotators which exhibit complex stellar velocity fields and often include stellar kinematically Distinct Cores (KDCs), and Fast Rotators which have regular velocity fields. Our complete sample of 260 ETGs leads to a new criterion to disentangle Fast and Slow Rotators which now includes a dependency on the apparent ellipticity (Epsilon). It separates the two classes significantly better than the previous prescription, and than a criterion based on V/Sigma: Slow Rotators and Fast Rotators have LambdaR lower and larger than kFSxSQRT(Epsilon), respectively, where kFS=0.31 for measurements made within 1 Re. We show that the vast majority of early-type galaxies are Fast Rotators: these have regular stellar rotation, with aligned photometric and kinematic axes (Paper II, Krajnovic et al. 2011}, include discs and often bars and represent 86% (224/260) of all early-type galaxies in the volume-limited ATLAS3D sample. Fast Rotators span the full range of apparent ellipticities from 0 to 0.85, and we suggest that they cover intrinsic ellipticities from about 0.35 to 0.85, the most flattened having morphologies consistent with spiral galaxies. Only a small fraction of ETGs are Slow Rotators representing 14% (36/260) of the ATLAS3D sample of ETGs. Of all Slow Rotators, 11% (4/36) exhibit two counter-rotating stellar disc-like components and are rather low mass objects (Mdyn<10^10.5 M_Sun). All other Slow Rotators (32/36) appear relatively round on the sky (Epsilon_e<0.4), tend to be massive (Mdyn>10^10.5 M_Sun), and often (17/32) exhibit KDCs.

preprint2011arXiv

The Atlas3D project - V. The CO Tully-Fisher relation of early-type galaxies

We demonstrate here using both single-dish and interferometric observations that CO molecules are an excellent kinematic tracer, even in high-mass galaxies, allowing us to investigate for the first time the CO Tully-Fisher relation of early-type galaxies. We compare the Tully-Fisher relations produced using both single-dish and interferometric data and various inclination estimation methods, and evaluate the use of the velocity profile shape as a criterion for selecting galaxies in which the molecular gas extends beyond the peak of the rotation curve. We show that the gradient and zero-point of the best-fit relations are robust, independent of the velocity measure and inclination used, and agree with those of relations derived using stellar kinematics. We also show that the early-type CO Tully-Fisher (CO-TF) relation is offset from the CO-TF of spirals by 0.98 \pm 0.22 magnitude at Ks-band, in line with other results. The intrinsic scatter of the relation is found to be \approx 0.4 magnitudes, similar to the level found in the spiral galaxy population. Next generation facilities such as the Large Millimeter Telescope (LMT) and the Atacama Large Millimeter/Sub-millimeter Array (ALMA) should allow this technique to be used in higher-redshift systems, providing a simple new tool to trace the mass-to-light ratio evolution of the most massive galaxies over cosmic time.

preprint2011arXiv

The Atlas3D Project - VI. Simulations of binary galaxy mergers and the link with Fast Rotators, Slow Rotators, and Kinematically Distinct Cores

We study the formation of early-type galaxies through mergers with a sample of 70 high-resolution (softening length < 60 pc and 12*10^6 particles) numerical simulations of binary mergers of disc galaxies and 16 simulations of ETG remergers. These simulations, designed to accompany observations and models conducted within the Atlas3D project, encompass various mass ratios (from 1:1 to 6:1), initial conditions and orbital parameters. The progenitor disc galaxies are spiral-like with bulge to disc ratios typical of Sb and Sc galaxies. We find that binary mergers of disc galaxies with mass ratios of 3:1 and 6:1 are nearly always classified as Fast Rotators according to the Atlas3D criterion (based on the lambda_R parameter): they preserve the structure of the input fast rotating spiral progenitors. They have intrinsic ellipticities larger than 0.5, cover intrinsic lambda_R values between 0.2 and 0.6, within the range of observed Fast Rotators. Major disc mergers (mass ratios of 2:1 and 1:1) lead to both Fast and Slow Rotators. Most of the Fast Rotators produced in major mergers have intermediate flattening, with ellipticities between 0.4 and 0.6. Most Slow Rotators formed in these binary disc mergers hold a stellar Kinematically Distinct Core (KDC) in their 1-3 central kilo-parsec: these KDCs are built from the stellar components of the progenitors. Besides a handful of specific observed systems -- the counter-rotating discs (2-sigma galaxies) -- these therefore cannot reproduce the observed population of Slow Rotators in the nearby Universe. The mass ratio of the progenitors is a fundamental parameter for the formation of Slow Rotators in these binary mergers, but it also requires a retrograde spin for the earlier-type (Sb) progenitor galaxy with respect to the orbital angular momentum. (Abridged)

preprint2011arXiv

The Atlas3D project - VII. A new look at the morphology of nearby galaxies: the kinematic morphology-density relation

We look at the morphology of fast and slow rotator early-type galaxies. Edge-on fast rotators are lenticular galaxies. They appear like spiral galaxies with the gas and dust removed, and in some cases are flat ellipticals with disky isophotes. Fast rotators are often barred and span the same full range of bulge fractions as spiral galaxies. The slow rotators are rounder 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. We study for the first time the kinematic morphology-density T-Sigma 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 [f(SR)<2%] and generally constitute a small fraction [f(SR)~4%] of the total galaxy population in the environments explored by our survey, with the exception of the densest core of the Virgo cluster [f(SR)~20%]. We find a clean log-linear relation between the fraction f(Sp) of spiral galaxies and the local galaxy surface density. The existence of a smooth kinematic T-Sigma relation in the field excludes processes related to the cluster environment 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. Only at the largest densities in the Virgo core does the f(Sp) relation break down and steepens sharply, while the fraction of slow-rotators starts to significantly increase. This suggests that a different mechanism is at work there. (Abridged)

preprint2011arXiv

The Atlas3D project - X. On the origin of the molecular and ionised gas in early-type galaxies

We make use of interferometric CO and HI observations, and optical integral-field spectroscopy to probe the origin of the molecular and ionised interstellar medium (ISM) in local early-type galaxies (ETGs). We find that 36\pm5% of our sample of fast rotating ETGs have their ionised 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, 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 ISM share a common origin. In addition, we find that the origin of the cold and warm gas in fast-rotating ETGs 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. In the field, however, kinematic misalignments between the stellar and gaseous components indicate that >46% of local fast-rotating ETGs have their gas supplied from external sources. We discuss several scenarios which could explain the environmental dichotomy, but 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 >8x10^10 Msun) always having kinematically aligned gas. (abridged)

preprint2011arXiv

The Atlas3D project -- IV: the molecular gas content of early-type galaxies

We have carried out a survey for 12CO J=1-0 and J=2-1 emission in the 260 early-type galaxies of the volume-limited Atlas3D sample, with the goal of connecting their star formation and assembly histories to their cold gas content. This is the largest volume-limited CO survey of its kind and is the first to include many Virgo Cluster members. Sample members are dynamically hot galaxies with a median stellar mass 3\times 10^{10} Msun; they are selected by morphology rather than colour, and the bulk of them lie on the red sequence. The overall CO detection rate is 56/259 = 0.22 \error 0.03, with no dependence on K luminosity and only a modest dependence on dynamical mass. There are a dozen CO detections among the Virgo Cluster members; statistical analysis of their H_2 mass distributions and their dynamical status within the cluster shows that the cluster's influence on their molecular masses is subtle at best, even though (unlike spirals) they seem to be virialized within the cluster. We suggest that the cluster members have retained their molecular gas through several Gyr residences in the cluster. There are also a few extremely CO-rich early-type galaxies with H_2 masses >= 10^9 Msun, and these are in low density environments. We do find a significant trend between molecular content and the stellar specific angular momentum. The galaxies of low angular momentum also have low CO detection rates, suggesting that their formation processes were more effective at destroying molecular gas or preventing its re-accretion. We speculate on the implications of these data for the formation of various sub-classes of early-type galaxies.

preprint2011arXiv

The Atlas3D Project -- XI. Dense molecular gas properties of CO-luminous early-type galaxies

Surveying eighteen 12CO-bright galaxies from the ATLAS3D early-type galaxy sample with the Institut de Radio Astronomie Millimétrique (IRAM) 30m telescope, we detect 13CO(1-0) and 13CO(2-1) in all eighteen galaxies, HCN(1-0) in 12/18 and HCO+(1-0) in 10/18. We find that the line ratios 12CO(1-0)/13CO(1-0) and 12CO(1-0)/HCN(1-0) are clearly correlated with several galaxy properties: total stellar mass, luminosity-weighted mean stellar age, molecular to atomic gas ratio, dust temperature and dust morphology. We suggest that these correlations are primarily governed by the optical depth in the 12CO lines; interacting, accreting and/or starbursting early-type galaxies have more optically thin molecular gas while those with settled dust and gas discs host optically thick molecular gas. The ranges of the integrated line intensity ratios generally overlap with those of spirals, although we note some outliers in the 12CO(1- 0)/13CO(1-0), 12CO(2-1)/13CO(2-1) and HCN/HCO+(1-0) ratios. In particular, three galaxies are found to have very low 12CO(1-0)/13CO(1-0) and 12CO(2-1)/13CO(2-1) ratios. Such low ratios may signal particularly stable molecular gas which creates stars less efficiently than 'normal' (i.e. below Schmidt-Kennicutt prediction), consistent with the low dust temperatures seen in these galaxies.

preprint2011arXiv

The Atlas3D project -- XIII. Mass and morphology of HI in early-type galaxies as a function of environment

We present the Atlas3D HI survey of 166 nearby early-type galaxies (ETGs) down to M(HI)~10^7 M_sun. We detect HI in ~40% of all ETGs outside the Virgo cluster and in ~10% of all ETGs inside it. This demonstrates that it is common for non-cluster ETGs to host HI. The HI morphology varies from regular discs/rings (the majority of the detections) to unsettled gas distributions. The former are either small discs (M(HI)<10^8 M_sun) confined within the stellar body and sharing the same kinematics of the stars, or large discs/rings (M(HI) up to 5x10^9 M_sun) extending to tens of kpc from the host galaxy and frequently kinematically decoupled from the stars. Neutral hydrogen provides material for star formation in ETGs. Galaxies with central HI exhibit signatures of star formation in ~70% of the cases, ~5 times more frequently than galaxies without central HI. The central ISM is dominated by molecular gas. In ETGs with a small gas disc the conversion of HI into H_2 is as efficient as in spirals. The ETG HI mass function has M*~2x10^9 M_sun and slope=-0.7. ETGs host much less HI than spirals as a family. However, a significant fraction of them is as HI-rich as spirals. The main difference between ETGs and spirals is that the former lack the high-column-density HI typical of the bright stellar disc of the latter. We find an envelope of decreasing M(HI) with increasing environment density. The gas-richest ETGs live in the poorest environments (where star-formation is more common), galaxies in the centre of Virgo have the lowest HI content, and the cluster outskirts are a transition region. We find an HI morphology-density relation. At low environment density HI is mostly distributed on large discs/rings. More disturbed HI morphologies dominate environment densities typical of rich groups, confirming the importance of processes occurring on a galaxy-group scale for the evolution of ETGs.

preprint2011arXiv

The ATLAS3D Project-- VIII: Modelling the Formation and Evolution of Fast and Slow Rotator Early-Type Galaxies within $Λ$CDM

We propose a simple model for the origin of fast and slow rotator early-type galaxies (ETG) within the hierarchical $Λ$CDM scenario, that is based on the assumption that the mass fraction of stellar discs in ETGs is a proxy for the specific angular momentum expressed via $λ_R$. Within our model we reproduce the fraction of fast and slow rotators as a function of magnitude in the \atl survey, assuming that fast rotating ETGs have at least 10% of their total stellar mass in a disc component. In agreement with \atl observations we find that slow rotators are predominantly galaxies with $ M_* > 10^{10.5}$ M$_{\odot}$ contributing $\sim 20%$ to the overall ETG population. We show in detail that the growth histories of fast and slow rotators are different, supporting the classification of ETGs into these two categories. Slow rotators accrete between $\sim 50% -90%$ of their stellar mass from satellites and their most massive progenitors have on average up to 3 major mergers during their evolution. Fast rotators in contrast, accrete less than 50% and have on average less than one major merger in their past. We find that the underlying physical reason for the different growth histories is the slowing down and ultimately complete shut-down of gas cooling in massive galaxies. Once cooling and associated star formation in disc stops, galaxies grow via infall from satellites. Frequent minor mergers thereby, destroy existing stellar discs via violent relaxation and also tend to lower the specific angular momentum of the main stellar body, lowering $λ_R$ into the slow rotator regime. Abridged...

preprint2011arXiv

The Planetary Nebulae Population in the Central Regions of M32: the SAURON view

Extragalactic Planetary Nebulae (PNe) are not only useful as distance signposts or as tracers of the dark-matter content of their host galaxies, but constitute also good indicators of the main properties of their parent stellar populations. Yet, so far, the properties of PNe in the optical regions of galaxies where stellar population gradients can be more extreme have remained largely unexplored, mainly because the detection of PNe with narrow-band imaging or slit-less spectroscopy is considerably hampered by a strong stellar background. Integral-field spectroscopy (IFS) can overcome this limitation, and here we present a study of the PN population in the nearby compact elliptical M32. Using SAURON data taken with just two 10-minutes-long pointings we have doubled the number of known PNe within the effective radius of M32, detecting PNe five times fainter than previously found in narrow-band images that collected nearly the same number of photons. Furthermore, by carefully accounting for the incompleteness of our survey we could conclude, despite having only 15 sources, that the central PNe population of M32 is consistent with the generally adopted shape for the PNe Luminosity Function and its typical normalization observed in early-type galaxies. Finally, owing to the proximity of M32 and to UV images taken with HST, we could identify the most likely candidates for the central star of a subset of our detected PNe and conclude that these stars are affected by substantial amounts of circumstellar dust extinction, a finding that could reconcile the intriguing discrepancy previously reported in M32 between model predictions and observations for the later stages of stellar evolution. Considering the modest time investment on a 4m-class telescope that delivered these results, this work illustrates the potential of future IFS studies for the central PNe population of early-type galaxies.

preprint2011arXiv

The SAURON project - XVIII. The integrated UV-linestrength relations of early-type galaxies

Using far (FUV) and near (NUV) ultraviolet photometry from guest investigator programmes on the Galaxy Evolution Explorer (GALEX) satellite, optical photometry from the MDM Observatory and optical integral-field spectroscopy from SAURON, we explore the UV-linestrength relations of the 48 nearby early-type galaxies in the SAURON sample. Identical apertures are used for all quantities, avoiding aperture mismatch. We show that galaxies with purely old stellar populations show well-defined correlations of the integrated FUV-V and FUV-NUV colours with the integrated Mgb and Hbeta absorption linestrength indices, strongest for FUV-NUV. Correlations with the NUV-V colour, Fe5015 index and stellar velocity dispersion are much weaker. These correlations put stringent constraints on the origin of the UV-upturn phenomenon in early-type galaxies, and highlight its dependence on age and metallicity. In particular, despite recent debate, we recover the negative correlation between FUV-V colour and Mg linestrength originally publicised by Burstein et al. (1988), which we refer to as the "Burstein relation", suggesting a positive dependence of the UV-upturn on metallicity. We argue that the scatter in the correlations is real, and present mild evidence that a strong UV excess is preferentially present in slow-rotating galaxies. We also demonstrate that most outliers in the correlations are galaxies with current or recent star formation, some at very low levels. We believe that this sensitivity to weak star formation, afforded by the deep and varied data available for the SAURON sample, explains why our results are occasionally at odds with other recent but shallower surveys. This is supported by the analysis of a large, carefully-crafted sample of more distant early-type galaxies from the Sloan Digital Sky Survey (SDSS), more easily comparable with current and future large surveys.

preprint2011arXiv

The SAURON Project - XX. The Spitzer [3.6] - [4.5] colour in early-type galaxies: colours, colour gradients and inverted scaling relations

We investigate the [3.6] - [4.5] Spitzer-IRAC colour behaviour of the early-type galaxies of the SAURON survey, a representative sample of 48 nearby ellipticals and lenticulars. We investigate how this colour, which is unaffected by dust extinction, can be used to constrain the stellar populations in these galaxies. We find a tight relation between the [3.6]-[4.5] colour and effective velocity dispersion, a good mass-indicator in early-type galaxies. Contrary to other colours in the optical and near-infrared, we find that the colours become bluer for larger galaxies. The relations are tighter when using the colour inside r_e, rather than the much smaller r_e/8 aperture, due to the presence of young populations in the central regions. We also obtain strong correlations between the [3.6]-[4.5] colour and 3 strong absorption lines (H beta, Mg b and Fe 5015). Comparing our data with the models of Marigo et al., which show that more metal rich galaxies are bluer, we can explain our results in a way consistent with results from the optical, by stating that larger galaxies are more metal rich. The blueing is caused by a strong CO absorption band, whose line strength increases strongly with decreasing temperature and which covers a considerable fraction of the 4.5 micron filter. In galaxies that contain a compact radio source, the [3.6]-[4.5] colour is generally slightly redder than in the other galaxies, indicating small amounts of either hot dust, non-thermal emission, or young stars near the center. We find that the large majority of the galaxies show redder colours with increasing radius. abbreviated...

preprint2010arXiv

A collisional origin for the Leo ring

Extended HI structures around galaxies are of prime importance to probe galaxy formation scenarios. The giant HI ring in the Leo group is one of the largest and most intriguing HI structures in the nearby Universe. Whether it consists of primordial gas, as suggested by the apparent absence of any optical counterpart and the absence of an obvious physical connection to nearby galaxies, or of gas expelled from a galaxy in a collision is actively debated. We present deep wide field-of-view optical images of the ring region obtained with MegaCam on the CFHT. They reveal optical counterparts to several HI and UV condensations along the ring, in the g', r', and i' bands, which likely correspond to stellar associations formed within the gaseous ring. Analyzing the spectral energy distribution of one of these star-forming regions, we found it to be typical for a star-forming region in pre-enriched tidal debris. We then use simulations to test the hypothesis that the Leo ring results from a head-on collision between Leo group members NGC 3384 and M96. According to our model which is able to explain, at least qualitatively, the main observational properties of the system, the Leo ring is consistent with being a collisional ring. It is thus likely another example of extended intergalactic gas made-up of pre-enriched collisional debris.

preprint2010arXiv

The Atlas3D project -- I. A volume-limited sample of 260 nearby early-type galaxies: science goals and selection criteria

The Atlas3D project is a multi-wavelength survey combined with a theoretical modeling effort. The observations span from the radio to the millimeter and optical, and provide multi-colour imaging, two-dimensional kinematics of the atomic (HI), molecular (CO) and ionized gas (Hbeta, [OIII] and [NI]), together with the kinematics and population of the stars (Hbeta, Fe5015 and Mgb), for a carefully selected, volume-limited (1.16*10^5 Mpc^3) 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) morphologically-selected ETGs extracted from a parent sample of 871 galaxies (8% E, 22% S0 and 70% spirals) brighter than M_K<-21.5 mag (stellar mass M_Star>6*10^9 M_Sun). We analyze 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 give an overview of the characteristics of the other main datasets already available for our sample and of the ongoing modelling projects.

preprint2010arXiv

The SAURON project - XVII. Stellar population analysis of the absorption line strength maps of 48 early-type galaxies

(Abridged) 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 Hbeta, Fe5015, and Mgb, measured in the Lick/IDS system and spatially binned to a constant signal-to-noise, together with predictions from up-to-date stellar population models, we estimate the simple stellar population-equivalent (SSP-equivalent) age, metallicity and abundance ratio [alpha/Fe] over a two-dimensional field extending up to approximately one effective radius. We find a large range of SSP-equivalent ages in our sample, of which ~40% 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(sigma_e <= 100 k/ms or ~2x10^10 M_sol). 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 disk/ring, also seen in the near-UV or mid-IR. The flattened components with disk-like kinematics previously identified in all fast rotators (Krajnović et al.) are shown to be connected to regions of distinct stellar populations. These range from the young, still star-forming circumnuclear disks and rings with increased metallicity preferentially found in intermediate-mass fast rotators, to apparently old structures with extended disk-like kinematics, which are observed to have an increased metallicity and mildly depressed [alpha/Fe] ratio compared to the main body of the galaxy. The slow rotators generally show no stellar population signatures over and above the well known metallicity gradients and are largely consistent with old (>=10 Gyr) stellar populations.

preprint2009arXiv

OASIS integral-field spectroscopy of the central kpc in 11 Seyfert 2 galaxies

We map narrow-line regions (NLRs) of 11 nearby Seyfert 2 galaxies with the optical integral-field spectrograph OASIS mounted at CFHT. We model emission-line profiles of 5 forbidden-line doublets and 2 Balmer lines, correcting for the underlying stellar absorption by reconstructing stellar spectra with synthetic evolutionary stellar population models. For each of the 11 targets, we present 2D maps of surface brightness in the observed emission lines, diagnostic line intensity ratios, gas kinematics (mean line-of-sight velocity and velocity dispersion), electron density, and interstellar reddening, and we plot spatially resolved spectral-diagnostic diagrams. The stellar data are represented by maps of mean line-of-sight (LOS) velocities and of the relative mass fractions of the young stellar populations. The gas velocity fields in 80% of the sample exhibit twisted S-shaped isovelocity contours, which are signatures of non-circular orbits and indicate non-axisymmetric gravitational potentials, gas motions out of the galactic plane, or possible outflows and inflows. Based on the kinematic measurements, we identified a possible nuclear ring or radial gas flow in NGC 262 (Mrk 348), not reported before. Eight of the eleven observed objects exhibit strongly asymmetric or multi-component emission-line profiles, in most cases confined to an elongated region passing through the galactic centre, perpendicular to the major axis of emission.

preprint2009arXiv

The SAURON project - XIII. SAURON-GALEX study of early-type galaxies: the ultraviolet colour-magnitude relations and Fundamental Planes

We present GALEX far (FUV) and near (NUV) ultraviolet imaging of 34 nearby early-type galaxies from the SAURON representative sample of 48 E/S0 galaxies, all of which have ground-based optical imaging from the MDM Observatory. The surface brightness profiles of nine galaxies show regions with blue UV-optical colours suggesting recent star formation. Five of these show blue integrated UV-optical colours that set them aside in the NUV integrated colour-magnitude relation. They also have other properties confirming they have had recent star formation, in particular H_beta absorption higher than expected for a quiescent population and a higher CO detection rate. NUV-blue galaxies are generally drawn from the lower stellar velocity dispersion and thus lower dynamical mass part of the sample. We have also constructed the first UV Fundamental Planes and show that NUV blue galaxies bias the slopes and increase the scatters. If they are eliminated the fits get closer to expectations from the virial theorem. Although our analysis is based on a limited sample, it seems that a dominant fraction of the tilt and scatter of the UV Fundamental Planes is due to the presence of young stars in preferentially low-mass early-type galaxies. Interestingly, the UV-optical radial colour profiles reveal a variety of behaviours, with many galaxies showing signs of recent star formation, a central UV-upturn phenomenon, smooth but large-scale age and metallicity gradients, and in many cases a combination of these. In addition, FUV-NUV and FUV-V colours even bluer than those normally associated with UV-upturn galaxies are observed at the centre of some quiescent galaxies. Four out of the five UV-upturn galaxies are slow rotators. These objects should thus pose interesting challenges to stellar evolutionary models of the UV-upturn.

preprint2009arXiv

The SAURON project XVI: On the Sources of Ionisation for the Gas in Elliptical and Lenticular Galaxies

Following our study on the incidence, morphology and kinematics of the ionised 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, 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 [OIII]/Hb and [NI]/Hb line ratios. We find a tight correlation between the stellar surface brightness and the flux of the Hb recombination line across our sample, which points to a diffuse and old stellar source as the main contributor of ionising photons in early-type galaxies, with post-asymptotic giant branch (pAGB) stars being still the best candidate based on ionising-balance arguments. Other ionising sources such as a central AGN, OB-stars, shocks or the interaction between the hot and warm phases of the interestellar medium are found to play only a limited or localised role in powering the diffuse nebular emission observed in our sample galaxies. These results lead us to investigate the relative importance of stellar and AGN photoionisation in explaining the ionised-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 early-type galaxies that display modest values for the equivalent width of the [OIII] line (less than ~2.4Å) and LINER-like [OIII]/Hb values, the nebular emission is truly powered by an AGN.

preprint2008arXiv

Dynamical Evolution of AGN Host Galaxies - Gas In/Out-Flow Rates in 7 NUGA Galaxies

To examine the role of the host galaxy structure in fueling nuclear activity, we estimated gas flow rates from several kpc down to the inner few 10 pc for seven nearby spiral galaxies, selected from the NUGA sample (NUclei of GAlaxies). We calculated gravitational torques from near-IR images and determined gas in/out-flow rates as a function of radius and location within the galactic disks, based on high angular resolution interferometric observations of molecular (CO using PdBI) and atomic (HI using the VLA) gas. The results are compared with kinematic evidence for radial gas flows and the dynamical state of the galaxies (via resonances) derived from several different methods. We show that gravitational torques are very efficient at transporting gas from the outer disk all the way into the galaxies centers at ~100 pc; previously assumed dynamical barriers to gas transport, such as the Corotation Resonance of stellar bars, seem to be overcome by gravitational torque induced gas flows from other non-axisymmmetric structures. The resulting rates of gas mass inflow range from 0.01 to 50 solar masses per year and are larger for the galaxy center than for the outer disk. Our gas flow maps show the action of nested bars within larger bars for 3 galaxies. Non-circular streaming motions found in the kinematic maps are larger in the center than in the outer disk and appear to correlate only loosely with the in/out-flow rates as a function of radius. We demonstrate that spiral gas disks are very dynamic systems that undergo strong radial evolution on timescales of a few rotation periods (e.g. 5 times 10^8 yrs at a radius of 5 kpc), due to the effectiveness of gravitational torques in redistributing the cold galactic gas.

preprint2008arXiv

Specific Angular Momentum of Disc Merger Remnants and the $λ_R$-Parameter

We use two-dimensional kinematic maps of simulated binary disc mergers to investigate the $λ_R$-parameter, which is a luminosity weighted measure of projected angular momentum per unit mass. This parameter was introduced to subdivide the SAURON sample of early type galaxies in so called fast $λ_R > 0.1$ and slow rotators $λ_R < 0.1$. Tests on merger remnants reveal that $λ_R$ is a robust indicator of the true angular momentum content in elliptical galaxies. We find the same range of $λ_R$ values in our merger remnants as in the SAURON galaxies. The merger mass ratio is decisive in creating a slow or a fast rotator in a single binary merger, the former being created mostly in an equal mass merger. Slow rotators have a $λ_R$ which does not vary with projection. The confusion rate with face-on fast rotators is very small. Merger with low gas fractions form slow rotators with smaller ellipticities and are in much better agreement with the SAURON slow rotators. Remergers of merger remnants are slow rotators but tend to have too high ellipticities. Fast rotators maintain the angular momentum content from the progenitor disc galaxy if merger mass ratio is high. Some SAURON galaxies have values of $λ_R$ as high as our progenitor disc galaxies.

preprint2007arXiv

Integral Field Spectrographs: a user's view

We easily tend to think of Integral-Field Spectrographs (IFS) along two opposing trends: as either the beautiful combination between photometry and spectroscopy, or as our worst nightmare including the dark side of both worlds. I favour a view where each IFS is considered individually, as one instrument with specific performances which can be used optimally for a certain range of scientific programs. It is indeed true that data-wise, IFS do sometime merge the characteristics of classic (e.g., long-slit) spectrographs with annoying issues associated with Imagers. This is in fact the price to pay to access a drastically different perspective of our favourite targets. The challenge is then to provide the necessary tools to properly handle the corresponding data. However, this should certainly not be thought as something specific to IFS: such a challenge should be accepted for any instrument, and most importantly solved prior to its delivery at the telescope.

preprint2007arXiv

Scalable N-body code for the modelling of early-type galaxies

Early-type galaxies exhibit a wealth of photometric and dynamical structures. These signatures are fossil records of their formation and evolution processes. In order to examine these structures in detail, we build models aimed at reproducing the observed photometry and kinematics. The developed method is a generalization of the one introduced by Syer and Tremaine (1996), consisting in an N-body representation, in which the weights of the particles are changing with time. Our code is adapted for integral-field spectroscopic data, and is able to reproduce the photometric as well as stellar kinematic data of observed galaxies. We apply this technique on SAURON data of early-type galaxies, and present preliminary results on NGC 3377.

preprint2006arXiv

Young Kinematically Decoupled Components in Early-Type Galaxies

We present results from a series of follow-up observations of a sub-sample of the representative SAURON survey elliptical (E) and lenticular (S0) galaxies using the OASIS integral-field spectrograph. These observations focus on the central 10" x 10", with roughly double the spatial resolution of the SAURON observations. This increased spatial resolution reveals a number of interesting and previously unresolved features in the measured stellar kinematics and absorption-line strengths. We find that galaxies exhibiting the youngest global stellar populations (as measured with SAURON) often contain a distinctly young central region (on scales of a few hundred parsec or less) compared to the rest of the galaxy. Moreover, these compact, young components are found to be mostly counter-rotating with respect to the rest of the galaxy. Given that there is no well-established reason for such young components to `prefer' counter- over co-rotation, this finding raises the following questions: How common are these small KDCs as a function of age? Why are there more young than old compact KDCs? Where are the equivalent co-rotating components? We explore these questions using simple simulated velocity fields and stellar population models, and find that the fading of the young component as it evolves, coupled with the fact that counter-rotating components are more easily detected in the velocity field, may help explain the observed trends.

preprint2005arXiv

A Bar Signature and Central Disk in the Gaseous and Stellar Velocity Fields of NGC5448

We analyse SAURON kinematic maps of the inner kpc of the early-type (Sa) barred spiral galaxy NGC5448. The observed morphology and kinematics of the emission-line gas is patchy and perturbed, indicating clear departures from circular motion. The kinematics of the stars is more regular, and display a small inner disk-like system embedded in a large-scale rotating structure. We focus on the [OIII] gas, and use a harmonic decomposition formalism to analyse the gas velocity field. The higher-order harmonic terms and the main kinematic features of the observed data are consistent with an analytically constructed simple bar model. The bar model is derived using linear theory, considering an m=2 perturbation mode, and with bar parameters which are consistent with the large-scale bar detected via imaging. We also study optical and near infra-red images to reveal the asymmetric extinction in NGC5448, and we recognise that some of the deviations between the data and the analytical bar model may be due to these complex dust features. Our study illustrates how the harmonic decomposition formalism can be used as a powerful tool to quantify non-circular motions in observed gas velocity fields.

preprint2003arXiv

Difficulties with Recovering The Masses of Supermassive Black Holes from Stellar Kinematical Data

We investigate the ability of three-integral, axisymmetric, orbit-based modeling algorithms to recover the parameters defining the gravitational potential (M/L ratio and black hole mass Mh) in spheroidal stellar systems using stellar kinematical data. We show that the potential estimation problem is generically under-determined when applied to long-slit kinematical data of the kind used in most black hole mass determinations to date. A range of parameters (M/L, Mh) can provide equally good fits to the data, making it impossible to assign best-fit values. We illustrate the indeterminacy using a variety of data sets derived from realistic models as well as published observations of the galaxy M32. In the case of M32, our reanalysis demonstrates that data published prior to 2000 are equally consistent with Mh in the range 1.5x10^6-5x10^6 solar masses, with no preferred value in that range. While the HST/STIS data for this galaxy may overcome the degeneracy in Mh, HST data for most galaxies do not resolve the black hole's sphere of influence and in these galaxies the degree of degeneracy allowed by the data may be substantial. We investigate the effect on the degeneracy of enforcing smoothness (regularization) constraints. However we find no indication that the true potential can be recovered simply by enforcing smoothness. For a given smoothing level, all solutions in the minimum-chisquare valley exhibit similar levels of noise. These experiments affirm that the indeterminacy is real and not an artifact associated with non-smooth solutions. (Abridged)