We estimate the mass of the intermediate-mass black hole at the heart of the dwarf elliptical galaxy NGC 404 using Atacama Large Millimeter/submillimeter Array (ALMA) observations of the molecular interstellar medium at an unprecedented linear resolution of ~0.5 pc, in combination with existing stellar kinematic information. These ALMA observations reveal a central disc/torus of molecular gas clearly rotating around the black hole. This disc is surrounded by a morphologically and kinematically complex flocculent distribution of molecular clouds, that we resolve in detail. Continuum emission is detected from the central parts of NGC 404, likely arising from the Rayleigh-Jeans tail of emission from dust around the nucleus, and potentially from dusty massive star-forming clumps at discrete locations in the disc. Several dynamical measurements of the black hole mass in this system have been made in the past, but they do not agree. We show here that both the observed molecular gas and stellar kinematics independently require a ~5x10$^5$ Msun black hole once we include the contribution of the molecular gas to the potential. Our best estimate comes from the high-resolution molecular gas kinematics, suggesting the black hole mass of this system is 5.5$^{+4.1}_{-3.8}\times$10$^5$ Msun (at the 99% confidence level), in good agreement with our revised stellar kinematic measurement and broadly consistent with extrapolations from the black hole mass - velocity dispersion and black hole mass - bulge mass relations. This highlights the need to accurately determine the mass and distribution of each dynamically important component around intermediate-mass black holes when attempting to estimate their masses.

Galaxy clusters have long been theorised to quench the star-formation of their members. This study uses integral-field unit observations from the $K$-band Multi-Object Spectrograph (KMOS) - Cluster Lensing And Supernova survey with Hubble (CLASH) survey (K-CLASH) to search for evidence of quenching in massive galaxy clusters at redshifts $0.3<z<0.6$. We first construct mass-matched samples of exclusively star-forming cluster and field galaxies, then investigate the spatial extent of their H$\alpha$ emission and study their interstellar medium conditions using emission line ratios. The average ratio of H$\alpha$ half-light radius to optical half-light radius ($r_{\rm{e},\rm{H}\alpha}/r_{\rm{e},R_c}$) for all galaxies is $1.14\pm0.06$, showing that star formation is taking place throughout stellar discs at these redshifts. However, on average, cluster galaxies have a smaller $r_{\rm{e},\rm{H}\alpha}/r_{\rm{e},R_c}$ ratio than field galaxies: $\langle r_{\rm{e},\rm{H}\alpha}/r_{\rm{e},R_c}\rangle = 0.96\pm0.09$ compared to $1.22\pm0.08$ (smaller at a 98\% credibility level). These values are uncorrected for the wavelength difference between H$\alpha$ emission and $R_c$-band stellar light, but implementing such a correction only reinforces our results. We also show that whilst the cluster and field samples follow indistinguishable mass-metallicity (MZ) relations, the residuals around the MZ relation of cluster members correlate with cluster-centric distance; galaxies residing closer to the cluster centre tend to have enhanced metallicities (significant at the 2.6$\sigma$ level). Finally, in contrast to previous studies, we find no significant differences in electron number density between the cluster and field galaxies. We use simple chemical evolution models to conclude that the effects of disc strangulation and ram-pressure stripping can quantitatively explain our observations.

We present the first asymmetric drift (AD) measurements for unresolved stellar populations of different characteristic ages above and below 1.5 Gyr. These measurements sample the age-velocity relation (AVR) in galaxy disks. In this first paper we develop two efficient algorithms to extract AD on a spaxel-by-spaxel basis from optical integral-field spectroscopic (IFS) data-cubes. The algorithms apply different spectral templates, one using simple stellar populations and the other a stellar library; their comparison allows us to assess systematic errors in derived multi-component velocities, such as template-mismatch. We test algorithm reliability using mock spectra and Monte Carlo Markov Chains on real data from the MaNGA survey in SDSS-IV. We quantify random and systematic errors in AD as a function of signal-to-noise and stellar population properties with the aim of applying this technique to large subsets of the MaNGA galaxy sample. As a demonstration of our methods, we apply them to an initial sample of seven galaxies with comparable stellar mass and color to the Milky Way. We find a wide range of distinct AD radial profiles for young and old stellar populations.

We present the KMOS-CLASH (K-CLASH) survey, a K-band Multi-Object Spectrograph (KMOS) survey of the spatially-resolved gas properties and kinematics of 191 (predominantly blue) H$\alpha$-detected galaxies at $0.2 \lesssim z \lesssim 0.6$ in field and cluster environments. K-CLASH targets galaxies in four Cluster Lensing And Supernova survey with Hubble (CLASH) fields in the KMOS $IZ$-band, over $7'$ radius ($\approx2$-$3$ Mpc) fields-of-view. K-CLASH aims to study the transition of star-forming galaxies from turbulent, highly star-forming disc-like and peculiar systems at $z\approx1$-$3$, to the comparatively quiescent, ordered late-type galaxies at $z\approx0$, and to examine the role of clusters in the build-up of the red sequence since $z\approx1$. In this paper, we describe the K-CLASH survey, present the sample, and provide an overview of the K-CLASH galaxy properties. We demonstrate that our sample comprises star-forming galaxies typical of their stellar masses and epochs, residing both in field and cluster environments. We conclude K-CLASH provides an ideal sample to bridge the gap between existing large integral-field spectroscopy surveys at higher and lower redshifts. We find that star-forming K-CLASH cluster galaxies at intermediate redshifts have systematically lower stellar masses than their star-forming counterparts in the field, hinting at possible "downsizing" scenarios of galaxy growth in clusters at these epochs. We measure no difference between the star-formation rates of H$\alpha$-detected, star-forming galaxies in either environment after accounting for stellar mass, suggesting that cluster quenching occurs very rapidly during the epochs probed by K-CLASH, or that star-forming K-CLASH galaxies in clusters have only recently arrived there, with insufficient time elapsed for quenching to have occured.

Galaxy properties are known to correlate most tightly with the galaxy effective stellar velocity dispersion $\sigma_{\rm e}$. Here we look for {\em additional} trends at fixed $\sigma_{\rm e}$ using 1339 galaxies ($M_\ast \gtrsim 6\times10^9$ M$_\odot$) with different morphologies in the MaNGA (DR14) sample with integral-field spectroscopy data. We focus on the gradients ($\gamma_{\rm rms} \equiv \sigma(R_{\rm e}/4)/\sigma_{\rm e}$) of the stellar root-mean-square velocity ($V_{\rm rms} \equiv \sqrt{V^2 + \sigma^2}$), which we show traces the total mass density gradient $\gamma_{\rm tot}$ derived from dynamical models and, more weakly, the bulge fraction. We confirm that $\gamma_{\rm rms}$ increases with $\sigma_{\rm e}$, age and metallicity. We additionally find that these correlations still exist at fixed $\sigma_{\rm e}$, where galaxies with larger $\gamma_{\rm rms}$ are found to be older and more metal-rich. It means that mass density gradients contain information of the stellar population which is not fully accounted for by $\sigma_{\rm e}$. This result puts an extra constraint on our understanding of galaxy quenching. We compare our results with galaxies in the IllustrisTNG hydrodynamical simulations and find that, at fixed $\sigma_{\rm e}$, similar trends exist with age, the bulge fraction, and the total mass density slope but, unlike observations, no correlation with metallicity can be detected in the simulations.

We study the link between the kinematic-morphology of galaxies, as inferred from integral-field stellar kinematics, and their relation between mass and star formation rate (SFR). Our sample consists of $\sim 3200$ galaxies with integral-field spectroscopic data from the MaNGA survey with available determinations of their effective stellar angular momentum within the half-light radius $\lambda_{R_e}$. We find that for star-forming galaxies, namely along the star formation main sequence (SFMS), the $\lambda_{R_e}$ values remain large and almost unchanged over about two orders of magnitude in stellar mass, with the exception of the lowest masses $\mathcal{M}_{\star}\lesssim2\times10^{9} \mathcal{M}_{\odot}$, where $\lambda_{R_e}$ slightly decreases. The SFMS is dominated by spiral galaxies with small bulges. Below the SFMS, but above the characteristic stellar mass $\mathcal{M}_{\rm crit}\approx2\times10^{11} \mathcal{M}_{\odot}$, there is a sharp decrease in $\lambda_{R_e}$ with decreasing star formation rate: massive galaxies well below the SFMS are mainly slow-rotator early-type galaxies, namely genuinely spheroidal galaxies without disks. Below the SFMS and below $\mathcal{M}_{\rm crit}$ the decrease of $\lambda_{R_e}$ with decreasing SFR becomes modest or nearly absent: low-mass galaxies well below the SFMS, are fast-rotator early-type galaxies, and contain fast-rotating stellar disks like their star-forming counterparts. We also find a small but clear environmental dependence for the massive galaxies: in the mass range $10^{10.9}-10^{11.5} \mathcal{M}_{\odot}$, galaxies in rich groups or denser regions or classified as central galaxies have lower values of $\lambda_{R_e}$. While no environmental dependence is found for galaxies of lower mass. We discuss how our results can be understood as due to the different star formation and mass assembly histories of galaxies with varying mass.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

I present a flexible solution for the axisymmetric Jeans equations of stellar hydrodynamics under the assumption of an anisotropic (three-integral) velocity ellipsoid aligned with the spherical polar coordinate system. I describe and test a robust and efficient algorithm for its numerical computation. I outline the evaluation of the intrinsic velocity moments and the projection of all first and second velocity moments, including both the line-of-sight velocities and the proper motions. This spherically-aligned Jeans Anisotropic Modelling (JAM_sph) method can describe in detail the photometry and kinematics of real galaxies. It allows for a spatially-varying anisotropy, or stellar mass-to-light ratios gradients, as well as for the inclusion of general dark matter distributions and supermassive black holes. The JAM_sph method complements my previously derived cylindrically-aligned JAM_cyl and spherical Jeans solutions, which I also summarize in this paper. Comparisons between results obtained with either JAM_sph or JAM_cyl can be used to asses the robustness of inferred dynamical quantities. As an illustration, I modelled the Atlas3D sample of 260 early-type galaxies with high-quality integral-field spectroscopy, using both methods. I found that they provide statistically indistinguishable total-density logarithmic slopes. This may explain the previously-reported success of the JAM method in recovering density profiles of real or simulated galaxies. A reference software implementation of JAM_sph is included in the publicly-available JAM software package.

We construct the first comprehensive dynamical model for the high-quality subset of stellar kinematics of the Milky Way disc, with full 6D phase-space coordinates, provided by the Gaia Data Release 2. We adopt an axisymmetric approximation and use an updated Jeans Anisotropic Modelling (JAM) method, which allows for a generic shape and radial orientation of the velocity ellipsoid, as indicated by the Gaia data, to fit the mean velocities and all three components of the intrinsic velocity dispersion tensor. The Milky Way is the first galaxy for which all intrinsic phase space coordinates are available, and the kinematics are superior to the best integral-field kinematics of external galaxies. This situation removes the long-standing dynamical degeneracies and makes this the first dynamical model highly over-constrained by the kinematics. For these reasons, our ability to fit the data provides a fundamental test for both galaxy dynamics and the mass distribution in the Milky Way disc. We tightly constrain the volume average total density logarithmic slope, in the radial range 3.6--12 kpc, to be $\alpha_{\rm tot}=-2.149\pm 0.055$ and find that the dark halo slope must be significantly steeper than $\alpha_{\rm DM}=-1$ (NFW). The dark halo shape is close to spherical and its density is $\rho_{\rm DM}(R_\odot)=0.0115\pm0.0020$ M$_\odot$ pc$^{-3}$ ($0.437\pm0.076$ GeV cm$^{-3}$), in agreement with previous estimates. The circular velocity at the solar position $v_{\rm circ}(R_{\odot}) = 236.5\pm 3.1$ km s$^{-1}$ (including systematics) and its gently declining radial trends are also consistent with recent determinations.

We study a sample of 148 early-type galaxies in the Coma cluster using SDSS photometry and spectra, and calibrate our results using detailed dynamical models for a subset of these galaxies, to create a precise benchmark for dynamical scaling relations in high-density environments. For these galaxies, we successfully measured global galaxy properties, modeled stellar populations, and created dynamical models, and support the results using detailed dynamical models of 16 galaxies, including the two most massive cluster galaxies, using data taken with the SAURON IFU. By design, the study provides minimal scatter in derived scaling relations due to the small uncertainty in the relative distances of galaxies compared to the cluster distance. Our results demonstrate low ($\leq$55% for 90$^{th}$ percentile) dark matter fractions in the inner 1$R_{\rm e}$ ~of galaxies. Owing to the study design, we produce the tightest, to our knowledge, IMF-$\sigma_e$ relation of galaxies, with a slope consistent with that seen in local galaxies. Leveraging our dynamical models, we transform the classical Fundamental Plane of the galaxies to the Mass Plane. We find that the coefficients of the mass plane are close to predictions from the virial theorem, and have significantly lower scatter compared to the Fundamental plane. We show that Coma galaxies occupy similar locations in the (M$_*$ - $R_{\rm e}$) and (M$_*$ - $\sigma_e$) relations as local field galaxies but are older. This, and the fact we find only three slow rotators in the cluster, is consistent with the scenario of hierarchical galaxy formation and expectations of the kinematic morphology-density relation.

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, $\alpha$-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]

We present near-infrared observations of 42 gravitationally lensed galaxies obtained in the framework of the KMOS Lensed Emission Lines and VElocity Review (KLEVER) Survey, a program aimed at investigating the spatially resolved properties of the ionised gas in 1.2<z<2.5 galaxies by means of a full coverage of the YJ, H and K near-infrared bands. Detailed metallicity maps and gradients are derived for a sub-sample of 28 galaxies from reconstructed source plane emission line maps, exploiting the variety of different emission line diagnostics provided by the broad wavelength coverage of the survey. About 85% of these galaxies are characterised by metallicity gradients shallower than 0.05 dex/kpc and 89% are consistent with flat slope within 3$\sigma$ (67% within 1$\sigma$), suggesting a mild evolution with cosmic time. In the context of cosmological simulations and chemical evolution models, the presence of efficient feedback mechanisms and/or extended star formation profiles on top of the classical "inside-out" scenario of mass assembly is generally required to reproduce the observed flatness of the metallicity gradients beyond z$\sim$1 . Three galaxies with significantly (> 3$\sigma$) "inverted" gradients are also found, showing an anti-correlation between metallicity and star formation rate density on local scales, possibly suggesting recent episodes of pristine gas accretion or strong radial flows in place. Nevertheless, the individual metallicity maps are characterised by a variety of different morphologies, with flat radial gradients sometimes hiding non-axisymmetric variations on kpc scales which are washed out by azimuthal averages, especially in interacting systems or in those undergoing local episodes of recent star formation.

Galaxy interactions leave imprints in the motions of their stars, and so observing the two-dimensional stellar kinematics allows us to uncover their formation process. Slow rotators, which have stellar orbits dominated by random motions, are thought to be the fossil relics of a sequence of multiple gas-poor mergers, in an environment where the cold gas required to form new stars is nearly absent. Indeed, observations of a handful of nearby galaxy clusters have indicated that slow rotators are preferentially found in the gas-poor, dense cores of clusters, which themselves must form by merging of smaller groups. However, the generality of this result and connection between kinematics and environment is currently unclear, as recent studies have suggested that, at given stellar mass, the environment does not influence the formation of slow rotators. Here we address this issue by combining a careful quality-assessed sample selection with two-dimensional stellar kinematics from a large galaxy survey and a novel photometric classification approach where kinematics are unavailable. We obtain the first complete census of the location of massive slow rotators in ten large clusters: in all cases, slow rotators are extremely rare and generally trace the clusters density peaks. This result unambiguously establishes that massive slow rotators are the relics of violent hierarchical cluster formation.

We investigate the nature of star-forming galaxies with reduced specific star formation rate (sSFR) and high stellar masses, those `green valley' objects that seemingly cause a reported bending, or flattening, of the star-forming main sequence. The fact that such objects host large bulges recently led some to suggest that the internal formation of bulges was a late event that induced the sSFRs of massive galaxies to drop in a slow downfall, and thus the main sequence to bend. We have studied in detail a sample of 10 galaxies at $0.45<z<1$ with secure SFR from Herschel, deep Keck optical spectroscopy, and HST imaging from CANDELS allowing us to perform multi-wavelength bulge to disc decomposition, and to derive star formation histories for the separated bulge and disc components. We find that the bulges hosted in these systems below main sequence are virtually all maximally old, with ages approaching the age of the Universe at the time of observation, while discs are young ($\langle$ T$_{50}\rangle \sim 1.5$ Gyr). We conclude that, at least based on our sample, the bending of the main sequence is, for a major part, due to rejuvenation, and we disfavour mechanisms that postulate the internal formation of bulges at late times. The very old stellar ages of our bulges suggest a number density of Early Type Galaxies at $z=1-3$ higher than actually observed. If confirmed, this might represent one of the first direct validations of hierarchical assembly of bulges at high redshifts.

Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) is acquiring integral-field spectroscopy for the largest sample of galaxies to date. By 2020, the MaNGA Survey --- one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) --- will have observed a statistically representative sample of 10$^4$ galaxies in the local Universe ($z\lesssim0.15$). In addition to a robust data-reduction pipeline (DRP), MaNGA has developed a data-analysis pipeline (DAP) that provides higher-level data products. To accompany the first public release of its code base and data products, we provide an overview of the MaNGA DAP, including its software design, workflow, measurement procedures and algorithms, performance, and output data model. In conjunction with our companion paper Belfiore et al., we also assess the DAP output provided for 4718 observations of 4648 unique galaxies in the recent SDSS Data Release 15 (DR15). These analysis products focus on measurements that are close to the data and require minimal model-based assumptions. Namely, we provide stellar kinematics (velocity and velocity dispersion), emission-line properties (kinematics, fluxes, and equivalent widths), and spectral indices (e.g., D4000 and the Lick indices). We find that the DAP provides robust measurements and errors for the vast majority ($>$99%) of analyzed spectra. We summarize assessments of the precision and accuracy of our measurements as a function of signal-to-noise, and provide specific guidance to users regarding the limitations of the data. The MaNGA DAP software is publicly available and we encourage community involvement in its development.

The majority of galaxy group catalogues available in the literature use the popular friends-of-friends algorithm which links galaxies using a linking length. One potential drawback to this approach is that clusters of point can be link with thin bridges which may not be desirable. Furthermore, these algorithms are designed with large-scales galaxy surveys in mind rather than small-scale, local galaxy environments, where attention to detail is important. Here we present a new simple group finder algorithm, TD-ENCLOSER, that finds the group that encloses a target galaxy of interest. TD-ENCLOSER is based on the kernel density estimation method which treats each galaxy, represented by a zero-dimensional particle, as a two-dimensional circular Gaussian. The algorithm assigns galaxies to peaks in the density field in order of density in descending order ("Top Down") so that galaxy groups "grow" around the density peaks. Outliers in under-dense regions are prevented from joining groups by a specified hard threshold, while outliers at the group edges are clipped below a soft (blurred) interior density level. The group assignments are largely insensitive to all free parameters apart from the hard density threshold and the kernel standard deviation, although this is a known feature of density-based group finder algorithms, and operates with a computing speed that increases linearly with the size of the input sample. In preparation for a companion paper, we also present a simple algorithm to select unique representative groups when duplicates occur. TD-ENCLOSER produces results comparable to those from a widely used catalogue, as shown in a companion paper. A smoothing scale of 0.3 Mpc provides the most realistic group structure.

Galaxy morphology is inextricably linked to environment. The morphology-density relation quantifies this relationship. However, optical morphology is only loosely related to the kinematic structure of galaxies, and about two thirds of elliptical ("spheroidal") galaxies are actually misclassified face-on disks, and would appear flattened if view edge on. A more robust classification is the slow/fast rotator classification which describes the kinematic structure. Slow and fast rotators form a bimodality in galaxy properties and are thought to follow distinct evolutionary paths, and so a kinematic morphology-density (kT-$\Sigma$) relation is more meaningful. To date the kT-$\Sigma$ relation has only been studied for a handful of nearby clusters,or across large numbers of clusters but with incomplete coverage. Here, we combine stellar kinematics obtained with the Sloan Digital Sky Survey's (SDSS) Mapping Nearby Galaxies at Apache Point Observatory survey with classifications using a novel visual method to obtain the largest complete census of slow and fast rotators in nearby galaxy groups and clusters. To account for incompleteness in the SDSS spectroscopic catalogue, we combine the catalogue with the photometric catalogue, which we clean using empirical criteria. We test our visual classification method and find our false negative rate to be only $\sim$6% (0% for $M \gtrsim 10^{11.7}M_{\odot}$). In contrast, our false positive rate is about 50% implying half of all slow rotator candidates will be confirmed as fast rotators if stellar kinematics become available. Hence, our misclassification is essentially random with only a weak dependence on stellar mass, and as slow rotators are intrinsically rare, the absolute number of misclassifications will be small compared to the sample size. The result is the largest complete census of massive slow rotator candidates and their environments to date.

It has been observed that low redshift early-type galaxies can be separated into slow and fast rotators according to a proxy of specific stellar angular momentum, $\lambda_{R_e}$. Detailed studies of a handful of nearby clusters have shown that slow rotators are generally found at the centres of clusters where the number density is highest, whereas the fast rotators trace the trend followed by early-type galaxies of increasing in number with local density. In this paper, we study the environmental distribution of slow and fast rotators using the stellar kinematics of about 3900 galaxies from the Sloan Digital Sky Survey's Mapping Nearby Galaxies at Apache Point Observatory survey. For galaxies in groups closer than $z=0.08$ that are not observed with MaNGA but satisfy the necessary conditions for slow rotators, we visually assign slow/fast rotator classifications to obtain a complete sample. Our final catalogue contains about 900 groups of five or more members. We observe the kinematic morphology-density (kT-$\Sigma$) relation for each group and find an increasing fraction of massive slow rotators with increasing number density. We provide evidence suggesting that the observed lack of trends in angular momentum with environment at fixed stellar mass is partly because the maximum density varies between clusters, and that the locations of massive slow rotators are strongly correlated with peak densities in galaxy groups and clusters. We conclude that the (projected) number density relative to the cluster peak is more fundamental than the absolute number density in influencing the abundance of slow rotators. We find that the kT-$\Sigma$ relation does exist at fixed stellar mass, and we rule out the hypothesis that the kT-$\Sigma$ relation is a result of dynamical friction alone, instead arguing that massive slow rotators grow hierarchically in tandem with their host clusters.

As part of the mm-Wave Interferometric Survey of Dark Object Masses (WISDOM), we present a measurement of the mass of the supermassive black hole (SMBH) in the nearby early-type galaxy NGC 0383 (radio source 3C 031). This measurement is based on Atacama Large Millimeter/sub-milimeter Array (ALMA) cycle 4 and 5 observations of the 12CO(2-1) emission line with a spatial resolution of 58x32pc2 (0."18x0."1). This resolution, combined with a channel width of 10 km/s, allows us to well resolve the radius of the black hole sphere of influence (measured as R_SOI = 316pc = 0."98), where we detect a clear Keplerian increase of the rotation velocities. NGC 0383 has a kinematically-relaxed, smooth nuclear molecular gas disc with weak ring/spiral features. We forward-model the ALMA data cube with the Kinematic Molecular Simulation (KinMS) tool and a Bayesian Markov Chain Monte Carlo method to measure a SMBH mass of (4.2+/-0.7)x10^9 Msun, a F160W-band stellar mass-to-light ratio that varies from 2.8+/-0.6 Msun/Lsun in the centre to 2.4+/-0.3 Msun/Lsun at the outer edge of the disc and a molecular gas velocity dispersion of 8.3+/-2.1 km/s (all 3-sigma uncertainties). We also detect unresolved continuum emission across the full bandwidth, consistent with synchrotron emission from an active galactic nucleus. This work demonstrates that low-J CO emission can resolve gas very close to the SMBH (~140,000 Schwarzschild radii) and hence that the molecular gas method is highly complimentary to megamaser observations as it can probe the same emitting material.

SDSS-IV MaNGA (Mapping Nearby Galaxies at Apache Point Observatory) is the largest integral-field spectroscopy survey to date, aiming to observe a statistically representative sample of 10,000 low-redshift galaxies. In this paper we study the reliability of the emission-line fluxes and kinematic properties derived by the MaNGA Data Analysis Pipeline (DAP). We describe the algorithmic choices made in the DAP with regards to measuring emission-line properties, and the effect of our adopted strategy of simultaneously fitting the continuum and line emission. The effect of random errors are quantified by studying various fit-quality metrics, idealized recovery simulations and repeat observations. This analysis demonstrates that the emission lines are well-fit in the vast majority of the MaNGA dataset and the derived fluxes and errors are statistically robust. The systematic uncertainty on emission-line properties introduced by the choice of continuum templates is also discussed. In particular, we test the effect of using different stellar libraries and simple stellar-population models on the derived emission-line fluxes and the effect of introducing different tying prescriptions for the emission-line kinematics. We show that these effects can generate large ($>$ 0.2 dex) discrepancies at low signal-to-noise and for lines with low equivalent width (EW); however, the combined effect is noticeable even for H$\alpha$ EW $>$ 6~\AA. We provide suggestions for optimal use of the data provided by SDSS data release 15 and propose refinements on the \DAP\ for future MaNGA data releases.

Different massive black hole mass - host galaxy scaling relations suggest that the growth of massive black holes is entangled with the evolution of their host galaxies. The number of measured black hole masses is still limited, and additional measurements are necessary to understand the underlying physics of this apparent co-evolution. We add six new black hole mass (MBH) measurements of nearby fast rotating early-type galaxies to the known black hole mass sample, namely NGC 584, NGC 2784, NGC 3640, NGC 4570, NGC 4281 and NGC 7049. Our target galaxies have effective velocity dispersions ({\sigma}e) between 170 and 245 km s^(-1), and thus this work provides additional insight into the black hole properties of intermediate-mass early-type galaxies. We combine high-resolution adaptive-optics SINFONI data with large-scale MUSE, VIMOS and SAURON data from ATLAS3D to derive two-dimensional stellar kinematics maps. We then build both Jeans Anisotropic Models and axisymmetric Schwarzschild models to measure the central black hole masses. Our Schwarzschild models provide black hole masses which are consistent with recent MBH-{\sigma}e scaling relations. NGC 3640 has a velocity dispersion dip and NGC 7049 a constant velocity dispersion in the center, but we can clearly constrain their lower black hole mass limit. We conclude our analysis with a test on NGC 4570 taking into account a variable mass-to-light ratio (M/L) when constructing dynamical models. When considering M/L variations linked mostly to radial changes in the stellar metallicity, we find that the dynamically determined black hole mass from NGC 4570 decreases by 30%. Further investigations are needed in the future to account for the impact of radial M/L gradients on dynamical modeling.

We present high angular resolution (0.3" or 37 pc) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of the CO(2-1) line emission from a central disc in the early-type galaxy NGC 524. This disc is shown to be dynamically relaxed, exhibiting ordered rotation about a compact 1.3mm continuum source, which we identify as emission from an active supermassive black hole (SMBH). There is a hole at the centre of the disc slightly larger than the SMBH sphere of influence. An azimuthal distortion of the observed velocity field is found to be due to either a position angle warp or radial gas flow over the inner 2.5". By forward-modelling the observations, we obtain an estimate of the SMBH mass of $4.0^{+3.5}_{-2.0}\times10^8\,\mathrm{M_\odot}$, where the uncertainties are at the $3\sigma$ level. The uncertainties are dominated by the poorly constrained inclination and the stellar mass-to-light ratio of this galaxy, and our measurement is consistent with the established correlation between SMBH mass and stellar velocity dispersion. Our result is roughly half that of the previous stellar dynamical measurement, but is consistent within the uncertainties of both. We also present and apply a new tool for modelling complex molecular gas distributions.

We present the first measurement of the mass of a supermassive black hole (SMBH) in the nearby double-barred spiral galaxy NGC 3504 as part of the Measuring Black Holes Below the Milky Way ($M_{\star}$) mass galaxies (MBHBM$_{\star}$) Project. Our analysis is based on Atacama Large Millimeter/submillimeter Array (ALMA) Cycle-5 observations of the ${\rm ^{12}CO(2-1)}$ emission line. NGC 3504 has a circumnuclear gas disk (CND), which has a relatively high-velocity dispersion of 30 \kms. Our dynamical models of the CND yield a \Mbh~of $M_{\rm BH}=1.02^{+0.18}_{-0.15}\times10^7$\Msun and a mass-to-light ratio in $H$-band of \ml$_{\rm H}=0.66^{+1.44}_{-0.65}$ (\Msun/\Lsun). This black hole (BH) mass is consistent with BH--galaxy scaling relations. We also detect a central deficiency in the ${\rm ^{12}CO(2-1)}$ integrated intensity map with a diameter of 2.7 pc at the putative position of the SMBH. However, this hole is filled by a dense gas tracer ${\rm CS(5-4)}$ that peaks at the galaxy center found in one of the three low-velocity-resolution continuum spectral correlators. The ${\rm CS(5-4)}$ line has the same kinematics with the ${\rm ^{12}CO(2-1)}$ line within the CND, suggesting that it is also an alternative transition for measuring the central \Mbh~in NGC 3504 probably more accurately than the current commonly used of ${\rm ^{12}CO(2-1)}$ due to its centralization.

Three basic ingredients are required to generate a simple stellar population (SSP) library, i.e., an initial mass function (IMF), a stellar evolution model/isochrones, and an empirical/theoretical stellar spectral library. However, there are still some uncertainties to the determination and understanding of these ingredients. We perform the spectral fitting to test the relative parameter offsets between these uncertainties using two different stellar population models, two different empirical stellar libraries, two different isochrones, and the Salpeter and Chabrier IMFs. Based on these setups, we select five SSP libraries generated with the Galaxev/STELIB and Vazdekis/MILES models, and apply them to the pPXF full-spectrum fitting of both MaNGA and mock spectra. We find that: 1) Compared to the Galaxev/STELIB model, spectral fitting qualities with the Vazdekis/MILES model have significant improvements for those metal-rich (especially over-solar) spectra, which cause better reduced $\chi^2$ distributions and more precisely fitted absorption lines. This might due to the lack of metal rich stars in the empirical STELIB library, or code improvement of the Vazdekis model. 2) When applying the Vazdekis/MILES model for spectral fitting, the IMF variation will lead to not only a systematic offset in $M_*/L_r$, but also offsets in age and metallicity, and these offsets increase with increasing stellar population ages. However, the IMF-variation caused metallicity offsets disappear in the case of Galaxev/STELIB based libraries. 3) The Padova2000 model provides a better match to the MaNGA galaxy spectra at [M/H]$_L<-1.0$, while the BaSTI model match the local galaxy spectra better at [M/H]$_L>-1.0$. Current tests suggest that spectral fitting with the Vazdekis/MILES+BaSTI combination would be a better choice for local galaxies.

We improve the dynamical black hole (BH) mass estimates in three nearby low-mass early-type galaxies--NGC 205, NGC 5102, and NGC 5206. We use new \hst/STIS spectroscopy to fit the star formation histories of the nuclei in these galaxies, and use these measurements to create local color--mass-to-light ratio (\ml) relations. We then create new mass models from \hst~imaging and combined with adaptive optics kinematics, we use Jeans dynamical models to constrain their BH masses. The masses of the central BHs in NGC 5102 and NGC 5206 are both below one million solar masses and are consistent with our previous estimates, $9.12_{-1.53}^{+1.84}\times10^5$\Msun~and $6.31_{-2.74}^{+1.06}\times10^5$\Msun~(3$\sigma$ errors), respectively. However, for NGC 205, the improved models suggest the presence of a BH for the first time, with a best-fit mass of $6.8_{-6.7}^{+95.6}\times10^3$\Msun~(3$\sigma$ errors). This is the least massive central BH mass in a galaxy detected using any method. We discuss the possible systematic errors of this measurement in detail. Using this BH mass, the existing upper limits of both X-ray, and radio emissions in the nucleus of NGC 205 suggest an accretion rate $\lesssim$$10^{-5}$ of the Eddington rate. We also discuss the color--\mleff~relations in our nuclei and find that the slopes of these vary significantly between nuclei. Nuclei with significant young stellar populations have steeper color--\mleff~relations than some previously published galaxy color--\mleff~relations.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

We present a series of 20 cosmological zoom simulations of the formation of massive galaxies with and without a model for AGN feedback. Differences in stellar population and kinematic properties are evaluated by constructing mock integral field unit (IFU) maps. The impact of the AGN is weak at high redshift when all systems are mostly fast-rotating and disc-like. After $z \sim 1$ the AGN simulations result in lower mass, older, less metal rich and slower rotating systems with less disky isophotes - in general agreement with observations. Two-dimensional kinematic maps of in-situ and accreted stars show that these differences result from reduced in-situ star formation due to AGN feedback. A full analysis of stellar orbits indicates that galaxies simulated with AGN are typically more triaxial and have higher fractions of x-tubes and box orbits and lower fractions of z-tubes. This trend can also be explained by reduced late in-situ star formation. We introduce a global parameter, $\xi_3$ , to characterise the anti-correlation between the third-order kinematic moment $h_3$ and the line-of-sight velocity ($v_{los}/{\sigma}$), and compare to ATLAS$^{3D}$ observations. The kinematic asymmetry parameter $\xi_3$ might be a useful diagnostic for large integral field surveys as it is a kinematic indicator for intrinsic shape and orbital content.

A new class of quiescent galaxies harboring possible AGN-driven winds has been discovered using spatially resolved optical spectroscopy from the ongoing SDSS-IV MaNGA survey. These galaxies, termed "red geysers", constitute $5-10\%$ of the local quiescent population and are characterized by narrow bisymmetric patterns in ionized gas emission features. Cheung et al. argued that these galaxies host large-scale AGN-driven winds that may play a role in suppressing star formation at late times. In this work, we test the hypothesis that AGN activity is ultimately responsible for the red geyser phenomenon. We compare the nuclear radio activity of the red geysers to a matched control sample with similar stellar mass, redshift, rest frame $NUV-r$ color, axis ratio and presence of ionized gas. We have used the 1.4 GHz radio continuum data from VLA FIRST survey to stack the radio flux from the red geyser and control samples. In addition to a 3 times higher FIRST detection rate, we find that red geysers have a 5$\sigma$ higher level of average radio flux than control galaxies. After restricting to rest-frame $NUV-r$ color $>$ 5 and checking mid-IR WISE photometry, we rule out star formation contamination and conclude that red geysers are associated with more active AGN. Red geysers and a possibly-related class with disturbed H$\alpha$ emission account for 40\% of all radio-detected red galaxies with $\rm log~(M_\star/M_\odot) < 11$. Our results support a picture in which episodic AGN activity drives large-scale-relatively weak ionized winds that may provide a feedback mechanism for many early-type galaxies.

By inverting the distributions of galaxies' apparent ellipticities and misalignment angles (measured around the projected half-light radius $R_{\rm e}$) between their photometric and kinematic axes, we study the intrinsic shape distribution of 189 slow rotator early-type galaxies with stellar masses $2\times 10^{11} M_{\odot}<M_\ast<2\times 10^{12} M_{\odot}$, extracted from a sample of about 2200 galaxies with integral-field stellar kinematics from the DR14 of the SDSS-IV MaNGA IFU survey. Thanks to the large sample of slow rotators, Graham+18 showed that there is clear structure in the misalignment angle distribution, with two peaks at both $0^{\circ}$ and $90^{\circ}$ misalignment (characteristic of oblate and prolate rotation respectively). Here we invert the observed distribution from Graham+18. The large sample allows us to go beyond the known fact that slow rotators are weakly triaxial and to place useful constraints on their intrinsic triaxiality distribution (around $1R_{\rm e}$) for the first time. The shape inversion is generally non-unique. However, we find that, for a wide set of model assumptions, the observed distribution clearly requires a dominant triaxial-oblate population. For some of our models, the data suggest a hint for a minor triaxial-prolate population, but a dominant prolate population is ruled out.

The majority of spiral and elliptical galaxies in the Universe host very dense and compact stellar systems at their centres known as nuclear star clusters (NSCs). In this work we study the stellar populations and star formation histories (SFH) of the NSCs of six nearby galaxies with stellar masses ranging between $2$ and $8\times10^9~{\rm M_{\odot}}$ (four late-type spirals and two early-types) with high resolution spectroscopy. Our observations are taken with the X-Shooter spectrograph at the VLT. We make use of an empirical simple stellar population (SSP) model grid to fit composite stellar populations to the data and recover the SFHs of the nuclei. We find that the nuclei of all late-type galaxies experienced a prolonged SFH, while the NSCs of the two early-types are consistent with SSPs. The NSCs in the late-type galaxies sample appear to have formed a significant fraction of their stellar mass already more than $10$ Gyr ago, while the NSCs in the two early-type galaxies are surprisingly younger. Stars younger than $100$ Myr are present in at least two nuclei: NGC 247 and NGC 7793, with some evidence for young star formation in NGC 300's NSC. The NSCs of the spirals NGC 247 and NGC 300 are consistent with prolonged \in situ star formation with a gradual metallicity enrichment from $\sim-1.5$ dex more than $10$ Gyr ago, reaching super-Solar values few hundred Myr ago. NGC 3621 appears to be very metal rich already in the early Universe and NGC 7793 presents us with a very complex SFH, likely dominated by merging of various massive star clusters coming from different environments.

Using mock spectra based on Vazdekis/MILES library fitted within the wavelength region 3600-7350\AA, we analyze the bias and scatter on the resulting physical parameters induced by the choice of fitting algorithms and observational uncertainties, but avoid effects of those model uncertainties. We consider two full-spectrum fitting codes: pPXF and STARLIGHT, in fitting for stellar population age, metallicity, mass-to-light ratio, and dust extinction. With pPXF we find that both the bias in the population parameters and the scatter in the recovered logarithmic values follows the expected trend. The bias increases for younger ages and systematically makes recovered ages older, $M_*/L_r$ larger and metallicities lower than the true values. For reference, at S/N=30, and for the worst case ($t=10^8$yr), the bias is 0.06 dex in $M_*/L_r$, 0.03 dex in both age and [M/H]. There is no significant dependence on either E(B-V) or the shape of the error spectrum. Moreover, the results are consistent for both our 1-SSP and 2-SSP tests. With the STARLIGHT algorithm, we find trends similar to pPXF, when the input E(B-V)<0.2 mag. However, with larger input E(B-V), the biases of the output parameter do not converge to zero even at the highest S/N and are strongly affected by the shape of the error spectra. This effect is particularly dramatic for youngest age, for which all population parameters can be strongly different from the input values, with significantly underestimated dust extinction and [M/H], and larger ages and $M_*/L_r$. Results degrade when moving from our 1-SSP to the 2-SSP tests. The STARLIGHT convergence to the true values can be improved by increasing Markov Chains and annealing loops to the "slow mode". For the same input spectrum, pPXF is about two order of magnitudes faster than STARLIGHT's "default mode" and about three order of magnitude faster than STARLIGHT's "slow mode".

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

We present mass estimates of supermassive black holes in six nearby fast rotating early-type galaxies (NGC4339, NGC4434, NGC4474, NGC4551, NGC4578 and NGC4762) with effective stellar velocity dispersion around 100 km/s. We use near-infrared laser-guide adaptive optics observations with the GEMINI/NIFS to derive stellar kinematics in the galactic nuclei, and SAURON observations from the ATLAS3D Survey for large-scale kinematics. We build axisymmetric Jeans anisotropic models and axisymmetric Schwarzschild dynamical models. Both modelling approaches recover consistent orbital anisotropies and black hole masses within 1-2sigma confidence level, except for one galaxy for which the difference is just above the 3sigma level. Two black holes (NGC4339 and NGC4434) are amongst the largest outliers from the current black hole mass - velocity dispersion relation, with masses of $(4.3^{+4.8}_{-2.3})\times10^7$ and $(7.0^{+2.0}_{-2.8})\times10^7$ M$_\odot$, respectively ($3\sigma$ confidence level). The black holes in NGC4578 and NGC4762 lie on the scaling relation with masses of $(1.9^{+0.6}_{-1.4})\times10^7$ and $(2.3^{+0.9}_{-0.6})\times10^7$ M$_\odot$, respectively (3sigma confidence level). For two galaxies (NGC4474 and NGC4551) we are able to place upper limits on their black holes masses ($<7\times10^6$ and $<5\times10^6$ M$_\odot$, respectively, $3\sigma$ confidence level). The kinematics for these galaxies clearly indicate central velocity dispersion drops within a radius of 35 pc and 80 pc, respectively. These drops cannot be associated with cold stellar structures and our data do not have the resolution to exclude black holes with masses an order of magnitude smaller than the predictions. Parametrizing the orbital distribution in spherical coordinates, the vicinity of the black holes is characterized by isotropic or mildly tangential anisotropy.

We present a detailed study of the nuclear star clusters (NSCs) and massive black holes (BHs) of four of the nearest low-mass early-type galaxies: M32, NGC205, NGC5012, and NGC5206. We measure dynamical masses of both the BHs and NSCs in these galaxies using Gemini/NIFS or VLT/SINFONI stellar kinematics, Hubble Space Telescope (HST) imaging, and Jeans Anisotropic Models. We detect massive BHs in M32, NGC5102, and NGC5206, while in NGC205, we find only an upper limit. These BH mass estimates are consistent with previous measurements in M32 and NGC205, while those in NGC5102 and NGC5206 are estimated for the first time, and both found to be $<$$10^6~M_{\odot}$. This adds to just a handful of galaxies with dynamically measured sub-million $M_{\odot}$ central BHs. Combining these BH detections with our recent work on NGC404's BH, we find that 80\% (4/5) of nearby, low-mass ($10^9-10^{10}~M_{\odot}$; $\sigma_{\star}\sim20-70$ km/s) early-type galaxies host BHs. Such a high occupation fraction suggests the BH seeds formed in the early epoch of cosmic assembly likely resulted in abundant seeds, favoring a low-mass seed mechanism of the remnants, most likely from the first generation of massive stars. We find dynamical masses of the NSCs ranging from $2-80\times10^6~M_{\odot}$ and compare these masses to scaling relations for NSCs based primarily on photometric mass estimates. Color gradients suggest younger stellar populations lie at the centers of the NSCs in three of the four galaxies (NGC205, NGC5102, and NGC5206), while the morphology of two are complex and are best-fit with multiple morphological components (NGC5102 and NGC5206). The NSC kinematics show they are rotating, especially in M32 and NGC5102 ($V/\sigma_{\star}\sim0.7$).

We examine the internal properties of the most massive ultracompact dwarf galaxy (UCD), M59-UCD3, by combining adaptive optics assisted near-IR integral field spectroscopy from Gemini/NIFS, and Hubble Space Telescope (HST) imaging. We use the multi-band HST imaging to create a mass model that suggests and accounts for the presence of multiple stellar populations and structural components. We combine these mass models with kinematics measurements from Gemini/NIFS to find a best-fit stellar mass-to-light ratio ($M/L$) and black hole (BH) mass using Jeans Anisotropic Models (JAM), axisymmetric Schwarzschild models, and triaxial Schwarzschild models. The best fit parameters in the JAM and axisymmetric Schwarzschild models have black holes between 2.5 and 5.9 million solar masses. The triaxial Schwarzschild models point toward a similar BH mass, but show a minimum $\chi^2$ at a BH mass of $\sim 0$. Models with a BH in all three techniques provide better fits to the central $V_{rms}$ profiles, and thus we estimate the BH mass to be $4.2^{+2.1}_{-1.7} \times 10^{6}$ M$_\odot$ (estimated 1$\sigma$ uncertainties). We also present deep radio imaging of M59-UCD3 and two other UCDs in Virgo with dynamical BH mass measurements, and compare these to X-ray measurements to check for consistency with the fundamental plane of BH accretion. We detect faint radio emission in M59cO, but find only upper limits for M60-UCD1 and M59-UCD3 despite X-ray detections in both these sources. The BH mass and nuclear light profile of M59-UCD3 suggests it is the tidally stripped remnant of a $\sim$10$^{9-10}$ M$_\odot$ galaxy.

We measure $\lambda_{R_e}$, a proxy for galaxy specific stellar angular momentum within one effective radius, and the ellipticity, $\epsilon$, for about 2300 galaxies of all morphological types observed with integral field spectroscopy as part of the MaNGA survey, the largest such sample to date. We use the $(\lambda_{R_e}, \epsilon)$ diagram to separate early-type galaxies into fast and slow rotators. We also visually classify each galaxy according to its optical morphology and two-dimensional stellar velocity field. Comparing these classifications to quantitative $\lambda_{R_e}$ measurements reveals tight relationships between angular momentum and galaxy structure. In order to account for atmospheric seeing, we use realistic models of galaxy kinematics to derive a general approximate analytic correction for $\lambda_{R_e}$. Thanks to the size of the sample and the large number of massive galaxies, we unambiguously detect a clear bimodality in the $(\lambda_{R_e}, \epsilon)$ diagram which may result from fundamental differences in galaxy assembly history. There is a sharp secondary density peak inside the region of the diagram with low $\lambda_{R_e}$ and $\epsilon < 0.4$, previously suggested as the definition for slow rotators. Most of these galaxies are visually classified as non-regular rotators and have high velocity dispersion. The intrinsic bimodality must be stronger, as it tends to be smoothed by noise and inclination. The large sample of slow rotators allows us for the first time to unveil a secondary peak at +/-90 degrees in their distribution of the misalignments between the photometric and kinematic position angles. We confirm that genuine slow rotators start appearing above a stellar mass of 2\times10^{11} M_{\odot}$ where a significant number of high-mass fast rotators also exist.

We perform full spectrum fitting stellar population analysis and Jeans Anisotropic modelling (JAM) of the stellar kinematics for about 2000 early-type galaxies (ETGs) and spiral galaxies from the MaNGA DR14 sample. Galaxies with different morphologies are found to be located on a remarkably tight mass plane which is close to the prediction of the virial theorem, extending previous results for ETGs. By examining an inclined projection (`the mass-size' plane), we find that spiral and early-type galaxies occupy different regions on the plane, and their stellar population properties (i.e. age, metallicity and stellar mass-to-light ratio) vary systematically along roughly the direction of velocity dispersion, which is a proxy for the bulge fraction. Galaxies with higher velocity dispersions have typically older ages, larger stellar mass-to-light ratios and are more metal rich, which indicates that galaxies increase their bulge fractions as their stellar populations age and become enriched chemically. The age and stellar mass-to-light ratio gradients for low-mass galaxies in our sample tend to be positive ($\rm centre<outer$), while the gradients for most massive galaxies are negative. The metallicity gradients show a clear peak around velocity dispersion $\log_{10} \sigma_{\rm e}\approx 2.0$, which corresponds to the critical mass $\sim 3\times 10^{10}M_{\odot}$ of the break in the mass-size relation. Spiral galaxies with large mass and size have the steepest gradients, while the most massive ETGs, especially above the critical mass $M_{\rm crit}\ge 2\times 10^{11} M_{\odot}$, where slow rotator ETGs start dominating, have much flatter gradients. This may be due to differences in their evolution histories, e.g. mergers.

We present results on the structural properties of massive passive galaxies in three clusters at $1.39<z<1.61$ from the KMOS Cluster Survey. We measure light-weighted and mass-weighted sizes from optical and near-infrared Hubble Space Telescope imaging and spatially resolved stellar mass maps. The rest-frame $R$-band sizes of these galaxies are a factor of $\sim2-3$ smaller than their local counterparts. The slopes of the relation between the stellar mass and the light-weighted size are consistent with recent studies in clusters and the field. Their mass-weighted sizes are smaller than the rest frame $R$-band sizes, with an average mass-weighted to light-weighted size ratio that varies between $\sim0.45$ and $0.8$ among the clusters. We find that the median light-weighted size of the passive galaxies in the two more evolved clusters is $\sim24\%$ larger than for field galaxies, independent of the use of circularized effective radii or semi-major axes. These two clusters also show a smaller size ratio than the less evolved cluster, which we investigate using color gradients to probe the underlying $M_{*}/L_{H_{160}}$ gradients. The median color gradients are $\nabla{z-H} \sim-0.4$ mag dex$^{-1}$, twice the local value. Using stellar populations models, these gradients are best reproduced by a combination of age and metallicity gradients. Our results favor the minor merger scenario as the dominant process responsible for the observed galaxy properties and the environmental differences at this redshift. The environmental differences support that clusters experience accelerated structural evolution compared to the field, likely via an epoch of enhanced minor merger activity during cluster assembly.

We present data for 16 galaxies in the overdensity JKCS 041 at $z \simeq 1.80$ as part of the K-band Multi-Object Spectrograph (KMOS) Cluster Survey (KCS). With 20-hour integrations, we have obtained deep absorption-line spectra from which we derived velocity dispersions for seven quiescent galaxies. We combined photometric parameters derived from Hubble Space Telescope images with the dispersions to construct a fundamental plane (FP) for quiescent galaxies in JKCS 041. From the zero-point evolution of the FP, we derived a formation redshift for the galaxies of $z_{form} = 3.0\pm0.3$, corresponding to a mean age of $1.4\pm0.2$ Gyrs. We tested the effect of structural and velocity dispersion evolution on our FP zero point and found a negligible contribution when using dynamical mass-normalized parameters ($\sim 3\%$), but a significant contribution from stellar-mass-normalized parameters ($\sim 42 \%$). From the relative velocities of the galaxies, we probed the three-dimensional structure of these 16 confirmed members of JKCS 041, and found that a group of galaxies in the south west of the overdensity had systematically higher velocities. We derived ages for the galaxies in the different groups from the FP. We found the east-extending group had typically older galaxies ($2.1\substack{+0.3\\-0.2}$ Gyrs), than those in the south-west group ($0.3\pm0.2$ Gyrs). Although based on small numbers, the overdensity dynamics, morphology, and age results could indicate that JKCS 041 is in formation and may comprise two merging groups of galaxies. The result could link large-scale structure to ages of galaxies for the first time at this redshift.

We investigate the variation of black hole masses (Mbh) as a function of their host galaxy stellar mass (Mstar) and half-light radius (Re). We confirm that the scatter in Mbh within this plane is essentially the same as that in the Mbh - sigma relation, as expected from the negligible scatter reported in the virial mass estimator sigma_v^2=GxMstar/(5xRe). All variation in Mbh happens along lines of constant sigma_v on the (Mstar, Re) plane, or Mstar $\propto$ Re for Mstar <2x10^11 Msun. This trend is qualitatively the same as those previously reported for galaxy properties related to stellar populations, like age, metallicity, alpha enhancement, mass-to-light ratio and gas content. We find evidence for a change in the Mbh variation above the critical mass of Mcrit ~ 2x10^11 Msun. This behaviour can be explained assuming that Mbh in galaxies less massive than Mcrit can be predicted by the Mbh - sigma relation, while Mbh in more massive galaxies follow a modified relation which is also dependent on Mstar once Mstar >Mcrit. This is consistent with the scenario where the majority of galaxies grow through star formation, while the most massive galaxies undergo a sequence of dissipation-less mergers. In both channels black holes and galaxies grow synchronously, giving rise to the black hole - host galaxy scaling relations, but there is no underlying single relation that is universal across the full range of galaxy masses.

As part of the mm-Wave Interferometric Survey of Dark Object Masses (WISDOM) project we present an estimate of the mass of the supermassive black hole (SMBH) in the nearby fast-rotating early-type galaxy NGC4429, that is barred and has a boxy/peanut-shaped bulge. This estimate is based on Atacama Large Millimeter/submillimeter Array (ALMA) cycle-2 observations of the 12CO(3-2) emission line with a linear resolution of ~13 pc (0"18x0"14). NGC4429 has a relaxed, flocculent nuclear disc of molecular gas that is truncated at small radii, likely due to the combined effects of gas stability and tidal shear. The warm/dense 12CO(3-2) emitting gas is confined to the inner parts of this disc, likely again because the gas becomes more stable at larger radii, preventing star formation. The gas disc has a low velocity dispersion of 2.2$^{+0.68}_{-0.65}$ km/s. Despite the inner truncation of the gas disc, we are able to model the kinematics of the gas and estimate a mass of (1.5$\pm0.1^{+0.15}_{-0.35}$)$\times$10$^8$ Msun for the SMBH in NGC4429 (where the quoted uncertainties reflect the random and systematic uncertainties, respectively), consistent with a previous upper limit set using ionised gas kinematics. We confirm that the V-band mass-to-light ratio changes by ~30% within the inner 400 pc of NGC4429, as suggested by other authors. This SMBH mass measurement based on molecular gas kinematics, the sixth presented in the literature, once again demonstrates the power of ALMA to constrain SMBH masses.

We present dynamical measurements of the central mass-to-light ratio ($M/L$) of a sample of 27 low-mass early-type ATLAS$^{3D}$ galaxies. We consider all ATLAS$^{3D}$ galaxies with 9.7$<$log(M$_\star/$M$_\odot)$$<$10.5 in our analysis, selecting out galaxies with available high-resolution Hubble Space Telescope (HST) data, and eliminating galaxies with significant central color gradients or obvious dust features. We use the HST images to derive mass models for these galaxies and combine these with the central velocity dispersion values from ATLAS$^{3D}$ data to obtain a central dynamical $M/L$ estimate. These central dynamical $M/L$s are higher than dynamical $M/L$s derived at larger radii and stellar population estimates of the galaxy centers in $\sim$80\% of galaxies, with a median enhancement of $\sim$14\% and a statistical significance of 3.3$\sigma$. We show that the enhancement in the central $M/L$ is best described either by the presence of black holes in these galaxies or by radial IMF variations. Assuming a black hole model, we derive black hole masses for the sample of galaxies. In two galaxies, NGC 4458 and NGC 4660, the data suggests significantly over-massive BHs, while in most others only upper limits are obtained. We also show that the level of $M/L$ enhancements we see in these early-type galaxy nuclei are consistent with the larger enhancements seen in ultracompact dwarf galaxies (UCDs), supporting the scenario where massive UCDs are created by stripping galaxies of these masses.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

We present the analysis of the fundamental plane (FP) for a sample of 19 massive red-sequence galaxies ($M_{\star} >4\times10^{10} M_{\odot}$) in 3 known overdensities at $1.39<z<1.61$ from the KMOS Cluster Survey, a guaranteed time program with spectroscopy from the K-band Multi-Object Spectrograph (KMOS) at the VLT and imaging from the Hubble Space Telescope. As expected, we find that the FP zero-point in $B$ band evolves with redshift, from the value 0.443 of Coma to $-0.10\pm0.09$, $-0.19\pm0.05$, $-0.29\pm0.12$ for our clusters at $z=1.39$, $z=1.46$, and $z=1.61$, respectively. For the most massive galaxies ($\log M_{\star}/M_{\odot}>11$) in our sample, we translate the FP zero-point evolution into a mass-to-light-ratio $M/L$ evolution finding $\Delta \log M/L_{B}=(-0.46\pm0.10)z$, $\Delta \log M/L_{B}=(-0.52\pm0.07)z$, to $\Delta \log M/L_{B}=(-0.55\pm0.10)z$, respectively. We assess the potential contribution of the galaxies structural and stellar velocity dispersion evolution to the evolution of the FP zero-point and find it to be $\sim$6-35 % of the FP zero-point evolution. The rate of $M/L$ evolution is consistent with galaxies evolving passively. By using single stellar population models, we find an average age of $2.33^{+0.86}_{-0.51}$ Gyr for the $\log M_{\star}/M_{\odot}>11$ galaxies in our massive and virialized cluster at $z=1.39$, $1.59^{+1.40}_{-0.62}$ Gyr in a massive but not virialized cluster at $z=1.46$, and $1.20^{+1.03}_{-0.47}$ Gyr in a protocluster at $z=1.61$. After accounting for the difference in the age of the Universe between redshifts, the ages of the galaxies in the three overdensities are consistent within the errors, with possibly a weak suggestion that galaxies in the most evolved structure are older.

We observed twelve nearby HI -detected early-type galaxies (ETGs) of stellar mass $\sim 10^{10}M\odot \leq M_* \leq \sim 10^{11}M\odot$ with the Mitchell Integral-Field Spectrograph, reaching approximately three half-light radii in most cases. We extracted line-of-sight velocity distributions for the stellar and gaseous components. We find little evidence of transitions in the stellar kinematics of the galaxies in our sample beyond the central effective radius, with centrally fast-rotating galaxies remaining fast-rotating and centrally slow-rotating galaxies likewise remaining slow-rotating. This is consistent with these galaxies having not experienced late dry major mergers; however, several of our objects have ionised gas that is misaligned with respect to their stars, suggesting some kind of past interaction. We extract Lick index measurements of the commonly-used H$\beta$, Fe5015, Mg\, b, Fe5270 and Fe5335 absorption features, and we find most galaxies to have flat H$\beta$ gradients and negative Mg, b gradients. We measure gradients of age, metallicity and abundance ratio for our galaxies using spectral fitting, and for the majority of our galaxies find negative age and metallicity gradients. We also find the stellar mass-to-light ratios to decrease with radius for most of the galaxies in our sample. Our results are consistent with a view in which intermediate-mass ETGs experience mostly quiet evolutionary histories, but in which many have experienced some kind of gaseous interaction in recent times.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

We present here parameter distributions and fundamental scaling relations for 190 galaxies as part of the Spectroscopy and H-bang Imaging of Virgo cluster galaxies (SHIVir) survey. We find the distribution of galaxy velocities to be bimodal about $V_{\rm circ} \sim 125$ km ${\rm s^{-1}}$, hinting at the existence of dynamically unstable modes in the inner regions of galaxies. An analysis of the Tully-Fisher relation (TFR) of late-type galaxies (LTGs) and fundamental plane (FP) of early-type galaxies (ETGs) is also presented, yielding a compendium of galaxy scaling relations. The slope and zero-point of the Virgo TFR match those of field galaxies, while scatter differences likely reflect distinct evolutionary histories. The velocities minimizing scatter for the TFR and FP are measured at large apertures where the baryonic fraction becomes subdominant. While TFR residuals remain independent of any galaxy parameters, FP residuals (i.e. the FP "tilt") correlate strongly with the dynamical-to-stellar mass ratio, yielding stringent galaxy formation constraints. Furthermore, we construct a stellar-to-total mass relation (STMR) for ETGs and LTGs and find linear but distinct trends over the range $M_{*} = 10^{8-11} M_{\odot}$. Stellar-to-halo mass relations (SHMRs), which probe the extended dark matter halo, can be scaled down to masses estimated within the optical radius, showing a tight match with the Virgo STMR at low masses; however, possibly inadequate halo abundance matching prescriptions and broad radial scalings complicate this comparison at all masses. While ETGs appear to be more compact than LTGs of the same stellar mass in projected space, their mass-size relations in physical space are identical. The trends reported here call for validation through well-resolved numerical simulations.

We study trends in the slope of the total mass profiles and dark matter fractions within the central half-light radius of 258 early-type galaxies, using data from the volume-limited ATLAS$^{\mathrm{3D}}$ survey. We use three distinct sets of dynamical models, which vary in their assumptions and also allow for spatial variations in the stellar mass-to-light ratio, to test the robustness of our results. We confirm that the slopes of the total mass profiles are approximately isothermal, and investigate how the total-mass slope depends on various galactic properties. The most statistically-significant correlations we find are a function of either surface density, \(\Sigma_e\), or velocity dispersion, \(\sigma_e\). However there is evidence for a break in the latter relation, with a nearly universal logarithmic slope above \(\log_{10}[\sigma_e/(\si{km~s^{-1}})]\sim 2.1\) and a steeper trend below this value. For the 142 galaxies above that critical \(\sigma_e\) value, the total mass-density logarithmic slopes have a mean value \(\left\langle\gamma^\prime\right\rangle = -2.192 \pm 0.016\) (\(1\sigma\) error) with an observed rms scatter of only \(\sigma_{\gamma^\prime}=0.167 \pm 0.016\). Considering the observational errors, we estimate an intrinsic scatter of \(\sigma_{\gamma^\prime}^\mathrm{intr} \approx 0.15\). These values are broadly consistent with those found by strong lensing studies at similar radii and agree, within the tight errors, with values recently found at much larger radii via stellar dynamics or HI rotation curves (using significantly smaller samples than this work).

As a part of the mm-Wave Interferometric Survey of Dark Object Masses (WISDOM) project, we present an estimate of the mass of the supermassive black hole (SMBH) in the nearby fast-rotator early-type galaxy NGC 3665. We obtained Combined Array for Research in Millimeter Astronomy (CARMA) B and C array observations of the $^{12}$CO$(J=2-1)$ emission line with a combined angular resolution of $0".59$. We analysed and modelled the three-dimensional molecular gas kinematics, obtaining a best-fit SMBH mass $M_{\rm BH}=5.75^{+1.49}_{-1.18} \times 10^{8}$ $M_{\odot}$, a mass-to-light ratio at $H$-band $(M/L)_{H}=1.45\pm0.04$ $(M/L)_{\odot, H}$, and other parameters describing the geometry of the molecular gas disc (statistical errors, all at $3\sigma$ confidence). We estimate the systematic uncertainties on the stellar $M/L$ to be $\approx0.2$ $(M/L)_{\odot, H}$, and on the SMBH mass to be $\approx0.4\times10^{8}$ $M_{\odot}$. The measured SMBH mass is consistent with that estimated from the latest correlations with galaxy properties. Following our older works, we also analysed and modelled the kinematics using only the major-axis position-velocity diagram, and conclude that the two methods are consistent.

As part of the mm-Wave Interferometric Survey of Dark Object Masses (WISDOM) project, we present an estimate of the mass of the supermassive black hole (SMBH) in the nearby fast-rotating early-type galaxy NGC4697. This estimate is based on Atacama Large Millimeter/submillimeter Array (ALMA) cycle-3 observations of the 12CO(2-1) emission line with a linear resolution of 29 pc (0.53"). We find that NGC4697 hosts a small relaxed central molecular gas disc with a mass of 1.6x10^7 Msun, co-spatial with the obscuring dust disc visible in optical Hubble Space Telescope (HST) imaging. We also resolve thermal 1mm continuum emission from the dust in this disc. NGC4697 is found to have a very low molecular gas velocity dispersion, $\sigma_{gas}=1.65^{+0.68}_{-0.65}$ km/s. This seems to be partially because the giant molecular cloud mass function is not fully sampled, but other mechanisms such as chemical differentiation in a hard radiation field or morphological quenching also seem to be required. We detect a Keplerian increase of the rotation of the molecular gas in the very centre of NGC4697, and use forward modelling of the ALMA data cube in a Bayesian framework with the KINematic Molecular Simulation (KinMS) code to estimate a SMBH mass of ($1.3_{-0.17}^{+0.18})\times10^8$ Msun and an i-band mass-to-light ratio of $2.14_{-0.05}^{+0.04}$ Msun/Lsun (at the 99% confidence level). Our estimate of the SMBH mass is entirely consistent with previous measurements from stellar kinematics. This increases confidence in the growing number of SMBH mass estimates being obtained in the ALMA era.

We perform Jeans anisotropic modeling (JAM) on elliptical and spiral galaxies from the MaNGA DR13 sample. By comparing the stellar mass-to-light ratios estimated from stellar population synthesis (SPS) and from JAM, we find a similar systematic variation of the initial mass function (IMF) as in the earlier $\rm ATLAS^{3D}$ results. Early type galaxies (elliptical and lenticular) with lower velocity dispersions within one effective radius are consistent with a Chabrier-like IMF while galaxies with higher velocity dispersions are consistent with a more bottom heavy IMF such as the Salpeter IMF. Spiral galaxies have similar systematic IMF variations, but with slightly different slopes and larger scatters, due to the uncertainties caused by higher gas fractions and extinctions for these galaxies. Furthermore, we examine the effects of stellar mass-to-light ratio gradients on our JAM modeling, and find that the trends from our results becomes stronger after considering the gradients.

We present a sample of 34 spectroscopically confirmed BzK-selected ~1e11 Msun quiescent galaxies (pBzK) in the COSMOS field. The targets were initially observed with VIMOS on the VLT to facilitate the calibration of the photometric redshifts of massive galaxies at z >~ 1.5. Here we describe the reduction and analysis of the data, and the spectrophotometric properties of these pBzK galaxies. In particular, using a spatially resolved median 2D spectrum, we find that the fraction of stellar populations with ages <1 Gyr is at least 3 times higher in the outer regions of the pBzK galaxies than in their cores. This results in a mild age gradient of ~<0.4 Gyr over ~6 kpc and suggests either the occurrence of widespread rejuvenation episodes or that inside-out quenching played a role in the passivization of this galaxy population. We also report on low-level star formation rates derived from the [OII]3727A emission line, with SFR_OII = 3.7-4.5 Msun/yr. This estimate is confirmed by an independent measurement on a separate sample of similarly-selected quiescent galaxies in the COSMOS field, using stacked far-infrared data (SFR_FIR = 2-4 Msun/yr). This second, photometric sample also displays significant excess at 1.4 GHz, suggestive of the presence of radio-mode AGN activity.

To ascertain whether photometric decompositions of galaxies into bulges and disks are astrophysically meaningful, we have developed a new technique to decompose spectral data cubes into separate bulge and disk components, subject only to the constraint that they reproduce the conventional photometric decomposition. These decompositions allow us to study the kinematic and stellar population properties of the individual components and how they vary with position, in order to assess their plausibility as discrete elements, and to start to reconstruct their distinct formation histories. An initial application of this method to CALIFA integral field unit observations of three isolated S0 galaxies confirms that in regions where both bulge and disc contribute significantly to the flux they can be physically and robustly decomposed into a rotating dispersion-dominated bulge component, and a rotating low-dispersion disc component. Analysis of the resulting stellar populations shows that the bulges of these galaxies have a range of ages relative to their discs, indicating that a variety of processes are necessary to describe their evolution. This simple test case indicates the broad potential for extracting from spectral data cubes the full spectral data of a wide variety of individual galaxy components, and for using such decompositions to understand the interplay between these various structures, and hence how such systems formed.

I start by providing an updated summary of the penalized pixel-fitting (pPXF) method, which is used to extract the stellar and gas kinematics, as well as the stellar population of galaxies, via full spectrum fitting. I then focus on the problem of extracting the kinematic when the velocity dispersion $\sigma$ is smaller than the velocity sampling $\Delta V$, which is generally, by design, close to the instrumental dispersion $\sigma_{\rm inst}$. The standard approach consists of convolving templates with a discretized kernel, while fitting for its parameters. This is obviously very inaccurate when $\sigma<\Delta V/2$, due to undersampling. Oversampling can prevent this, but it has drawbacks. Here I present a more accurate and efficient alternative. It avoids the evaluation of the under-sampled kernel, and instead directly computes its well-sampled analytic Fourier transform, for use with the convolution theorem. A simple analytic transform exists when the kernel is described by the popular Gauss-Hermite parametrization (which includes the Gaussian as special case) for the line-of-sight velocity distribution. I describe how this idea was implemented in a significant upgrade to the publicly available pPXF software. The key advantage of the new approach is that it provides accurate velocities regardless of $\sigma$. This is important e.g. for spectroscopic surveys targeting galaxies with $\sigma\ll\sigma_{\rm inst}$, for galaxy redshift determinations, or for measuring line-of-sight velocities of individual stars. The proposed method could also be used to fix Gaussian convolution algorithms used in today's popular software packages.

We explore the nucleus of the nearby 10$^9M_{\odot}$~early-type galaxy (ETGs), NGC~404, using Hubble Space Telescope (HST)/STIS spectroscopy and WFC3 imaging. We first present evidence for nuclear variability in UV, optical, and infrared filters over a time period of 15~years. This variability adds to the already substantial evidence for an accreting black hole at the center of NGC~404. We then redetermine the dynamical black hole mass in NGC~404 including modeling of the nuclear stellar populations. We combine HST/STIS spectroscopy with WFC3 images to create a local color-$M/L$~relation derived from stellar population modeling of the STIS data. We then use this to create a mass model for the nuclear region. We use Jeans modeling to fit this mass model to adaptive optics (AO) stellar kinematic observations from Gemini/NIFS. From our stellar dynamical modeling, we find a 3$\sigma$ upper limit on the black hole mass of $1.5\times10^5M_{\odot}$. Given the accretion evidence for a black hole, this upper limit makes NGC~404 the lowest mass central black hole with dynamical mass constraints. We find that the kinematics of H$_2$ emission line gas show evidence for non-gravitational motions preventing the use of gas dynamical modeling to constrain the black hole mass. Our stellar population modeling also reveals that the central, counter-rotating region of the nuclear cluster is dominated by $\sim$1~Gyr old populations.

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.

The kinematics and stellar populations of the low-mass nearby S0 galaxy NGC 5102 are studied from integral field spectra taken with the Multi-Unit Spectroscopic Explorer (MUSE). The kinematic maps reveal for the first time that NGC 5102 has the characteristic 2 sigma peaks indicative of galaxies with counter-rotating discs. This interpretation is quantitatively confirmed by fitting two kinematic components to the observed spectra. Through stellar population analysis we confirm the known young stellar population in the centre and find steep age and metallicity gradients. We construct axisymmetric Jeans anisotropic models of the stellar dynamics to investigate the initial mass function (IMF) and the dark matter halo of the galaxy. The models show that this galaxy is quite different from all galaxies previously studied with a similar approach: even within the half-light radius, it cannot be approximated with the self-consistent mass-follows-light assumption. Including an NFW dark matter halo, we need a heavy IMF and a dark matter fraction of 0.37+-0.04 within a sphere of one effective radius to describe the stellar kinematics. The more general model with a free slope of the dark matter halo shows that slope and IMF are degenerate, but indicates that a light weight IMF (Chabrier-like) and a higher dark matter fraction, with a steeper (contracted) halo, fit the data better. Regardless of the assumptions about the halo profile, we measure the slope of the total mass density to be -1.75+-0.04. This is shallower than the slope of -2 of an isothermal halo and shallower than published slopes for more massive early type galaxies.

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.

Integral field unit spectrographs allow the 2D exploration of the kinematics and stellar populations of galaxies, although they are generally restricted to small fields-of-view. Using the large field-of-view of the DEIMOS multislit spectrograph on Keck and our Stellar Kinematics using Multiple Slits (SKiMS) technique, we are able to extract sky-subtracted stellar light spectra to large galactocentric radii. Here we present a new DEIMOS mask design named SuperSKiMS that explores large spatial scales without sacrificing high spatial sampling. We simulate a set of observations with such a mask design on the nearby galaxy NGC 1023, showing that the kinematic and metallicity measurements can reach radii where the galaxy surface brightness is several orders of magnitude fainter than the sky. Such a technique is also able to reproduce the kinematic and metallicity 2D distributions obtained from literature integral field spectroscopy in the innermost galaxy regions. In particular, we use the simulated NGC 1023 kinematics to model its total mass distribution to large radii, obtaining comparable results with those from published integral field unit observation. Finally, from new spectra of NGC 1023 we obtain stellar 2D kinematics and metallicity distributions that show good agreement with integral field spectroscopy results in the overlapping regions. In particular, we do not find a significant offset between our SKiMS and the ATLAS$^\rm{3D}$ stellar velocity dispersion at the same spatial locations.

Using high resolution $N$-body simulations, we recently reported that a dynamically cool inner disk embedded in a hotter outer disk can naturally generate a steady double-barred (S2B) structure. Here we study the kinematics of these S2B simulations, and compare them to integral-field observations from ATLAS$^{3D}$ and SAURON. We show that S2B galaxies exhibit several distinct kinematic features, namely: (1) significantly distorted isovelocity contours at the transition region between the two bars, (2) peaks in $\sigma_\mathrm{LOS}$ along the minor axis of inner bars, which we term "$\sigma$-humps", that are often accompanied by ring/spiral-like features of increased $\sigma_\mathrm{LOS}$, (3) $h_3-\bar{v}$ anti-correlations in the region of the inner bar for certain orientations, and (4) rings of positive $h_4$ when viewed at low inclinations. The most impressive of these features are the $\sigma$-humps, these evolve with the inner bar, oscillating in strength just as the inner bar does as it rotates relative to the outer bar. We show that, in cylindrical coordinates, the inner bar has similar streaming motions and velocity dispersion properties as normal large-scale bars, except for $\sigma_z$, which exhibits peaks on the minor axis, i.e., humps. These $\sigma_z$ humps are responsible for producing the $\sigma$-humps. For three well-resolved early-type S2Bs (NGC 2859, NGC 2950, and NGC 3941) and a potential S2B candidate (NGC 3384), the S2B model qualitatively matches the integral-field data well, including the "$\sigma$-hollows" previously identified. We also discuss the kinematic effect of a nuclear disk in S2Bs.

The dark matter (DM) haloes around spiral galaxies appear to conspire with their baryonic content: empirically, significant amounts of DM are inferred only below a universal characteristic acceleration scale. Moreover, the discrepancy between the baryonic and dynamical mass, which is usually interpreted as the presence of DM, follows a very tight mass discrepancy acceleration (MDA) relation. Its universality, and its tightness in spiral galaxies, poses a challenge for the DM interpretation and was used to argue in favour of MOdified Newtonian Dynamics (MOND). Here, we test whether or not this applies to early-type galaxies. We use the dynamical models of fast-rotator early-type galaxies by Cappellari et al. based on ATLAS$^{3D}$ and SLUGGS data, which was the first homogenous study of this kind, reaching ~4 $R_e$, where DM begins to dominate the total mass budget. We find the early-type galaxies to follow an MDA relation similar to spiral galaxies, but systematically offset. Also, while the slopes of the mass density profiles inferred from galaxy dynamics show consistency with those expected from their stellar content assuming MOND, some profiles of individual galaxies show discrepancies.

Quiescent galaxies with little or no ongoing star formation dominate the galaxy population above $M_{*}\sim 2 \times 10^{10}~M_{\odot}$, where their numbers have increased by a factor of $\sim25$ since $z\sim2$. Once star formation is initially shut down, perhaps during the quasar phase of rapid accretion onto a supermassive black hole, an unknown mechanism must remove or heat subsequently accreted gas from stellar mass loss or mergers that would otherwise cool to form stars. Energy output from a black hole accreting at a low rate has been proposed, but observational evidence for this in the form of expanding hot gas shells is indirect and limited to radio galaxies at the centers of clusters, which are too rare to explain the vast majority of the quiescent population. Here we report bisymmetric emission features co-aligned with strong ionized gas velocity gradients from which we infer the presence of centrally-driven winds in typical quiescent galaxies that host low-luminosity active nuclei. These galaxies are surprisingly common, accounting for as much as $10\%$ of the population at $M_* \sim 2 \times 10^{10}~ M_{\odot}$. In a prototypical example, we calculate that the energy input from the galaxy's low-level active nucleus is capable of driving the observed wind, which contains sufficient mechanical energy to heat ambient, cooler gas (also detected) and thereby suppress star formation.

We observed the lenticular galaxy NGC 3998 with the Mitchell Integral-Field Spectrograph and extracted line-of-sight velocity distributions out to 3 half-light radii. We constructed collisionless orbit models in order to constrain NGC 3998's dark and visible structure, using kinematics from both the Mitchell and SAURON instruments. We find NGC 3998 to be almost axisymmetric, seen nearly face on with a flattened intrinsic shape - i.e., a face-on fast-rotator. We find an I-band mass-to-light ratio of $4.7_{-0.45}^{+0.32}$ in good agreement with previous spectral fitting results for this galaxy. Our best-fit orbit model shows a both a bulge and a disc component, with a non-negligible counter-rotating component also evident. We find that relatively little dark matter is needed to model this galaxy, with an inferred dark mass fraction of just $(7.1^{+8.1}_{-7.1})\%$ within one half-light radius.

We report a tight linear relation between the HI circular velocity measured at 6 $R_{\rm e}$ and the stellar velocity dispersion measured within 1 $R_{\rm e}$ for a sample of 16 early-type galaxies with stellar mass between $10^{10}$ and $10^{11}$ $\mathrm{M}_\odot$. The key difference from previous studies is that we only use spatially resolved $v_\mathrm{circ}$(HI) measurements obtained at large radius for a sizeable sample of objects. We can therefore link a kinematical tracer of the gravitational potential in the dark-matter dominated outer regions of galaxies with one in the inner regions, where baryons control the distribution of mass. We find that $v_\mathrm{circ}$(HI) = 1.33 $\sigma_\mathrm{e}$ with an observed scatter of just 12 percent. This indicates a strong coupling between luminous and dark matter from the inner- to the outer regions of early-type galaxies, analogous to the situation in spirals and dwarf irregulars. The $v_\mathrm{circ}$(HI)-$\sigma_\mathrm{e}$ relation is shallower than those based on $v_\mathrm{circ}$ measurements obtained from stellar kinematics and modelling at smaller radius, implying that \vcirc\ declines with radius -- as in bulge-dominated spirals. Indeed, the value of $v_\mathrm{circ}$(HI) is typically 25 percent lower than the maximum $v_\mathrm{circ}$ derived at $\sim0.2\ R_\mathrm{e}$ from dynamical models. Under the assumption of power-law total density profiles $\rho \propto r^{-\gamma}$, our data imply an average logarithmic slope $\langle\gamma\rangle=2.18\pm0.03$ across the sample, with a scatter of 0.11 around this value. The average slope and scatter agree with recent results obtained from stellar kinematics alone for a different sample of early-type galaxies.

We analyse the sizes, colour gradients, and resolved stellar mass distributions for 36 massive and passive galaxies in the cluster XMMUJ2235-2557 at z=1.39 using optical and near-infrared Hubble Space Telescope imaging. We derive light-weighted S\'ersic fits in five HST bands ($i_{775},z_{850},Y_{105},J_{125},H_{160}$), and find that the size decreases by ~20% going from $i_{775}$ to $H_{160}$ band, consistent with recent studies. We then generate spatially resolved stellar mass maps using an empirical relationship between $M_{*}/L_{H_{160}}$ and $(z_{850}-H_{160})$ and use these to derive mass-weighted S\'ersic fits: the mass-weighted sizes are ~41% smaller than their rest-frame $r$-band counterparts compared with an average of ~12% at z~0. We attribute this evolution to the evolution in the $M_{*}/L_{H_{160}}$ and colour gradient. Indeed, as expected, the ratio of mass-weighted to light-weighted size is correlated with the $M_{*}/L$ gradient, but is also mildly correlated with the mass surface density and mass-weighted size. The colour gradients $(\nabla_{z-H})$ are mostly negative, with a median value of $\sim0.45$ mag dex$^{-1}$, twice the local value. The evolution is caused by an evolution in age gradients along the semi-major axis ($a$), with $\nabla_{age} = d \log(age) / d \log(a)$ $\sim-0.33$, while the survival of weaker colour gradients in old, local galaxies implies that metallicity gradients are also required, with $\nabla_{Z} = d \log(Z) / d \log(a)$ $\sim-0.2$. This is consistent with recent observational evidence for the inside-out growth of passive galaxies at high redshift, and favours a gradual mass growth mechanism, such as minor mergers.

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.

Observations of galaxy isophotes, longs-slit kinematics and high-resolution photometry suggested a possible dichotomy between two distinct classes of E galaxies. But these methods are expensive for large galaxy samples. Instead, integral-field spectroscopic can efficiently recognize the shape, dynamics and stellar population of complete samples of early-type galaxies (ETGs). These studies showed that the two main classes, the fast and slow rotators, can be separated using stellar kinematics. We showed there is a dichotomy in the dynamics of the two classes. The slow rotators are weakly triaxial and dominate above $M_{\rm crit}\approx2\times10^{11} M_\odot$. Below $M_{\rm crit}$, the structure of fast rotators parallels that of spiral galaxies. There is a smooth sequence along which, the metals content, the enhancement in $\alpha$-elements, and the "weight" of the stellar initial mass function, all increase with the CENTRAL mass density slope, or bulge mass fraction, while the molecular gas fraction correspondingly decreases. The properties of ETGs on galaxy scaling relations, and in particular the $(M_{\ast}, R_{\rm e})$ diagram, and their dependence on environment, indicate two main independent channels for galaxy evolution. Fast rotators ETGs start as star forming disks and evolve trough a channel dominated by gas accretion, bulge growth and quenching. While slow rotators assemble near the center of massive halos via intense star formation at high redshift, and remain as such for the rest of their evolution via a channel dominated by gas poor mergers. This is consistent with independent studies of the galaxies redshift evolution.

At present, the main challenge to the interpretation of variations in gravity-sensitive line strengths as driven by a non-universal initial mass function (IMF), lies in understanding the effect of other parameters describing unresolved stellar populations, such as elemental abundance ratios. We combine various TiO-based, IMF-sensitive indicators in the optical and NIR spectral windows, along with the FeH-based Wing-Ford band to break this degeneracy. We obtain a significant radial trend of the IMF slope in XSG1, a massive early-type galaxy (ETG), with velocity dispersion sigma~300km/s, observed with the VLT/X-SHOOTER instrument. In addition, we constrain both the shape and normalization of the IMF based only on a stellar population analysis. We robustly rule out a single power-law to describe the IMF, whereas a power law tapered off to a constant value at low masses (defined as a bimodal IMF) is consistent with all the observational spectroscopic data and with the stellar M/L constraints based on the Jeans Anisotropic Modelling method. The IMF in XSG1 is bottom-heavy in the central regions (corresponding to a bimodal IMF slope Gb~3, or a mass normalization mismatch parameter alpha~2), changing towards a standard Milky-Way like IMF (Gb~1.3; alpha~1) around half of the effective radius. This result, combined with previous observations of local IMF variations in massive ETGs, reflects the varying processes underlying the formation of this type of galaxies, between the central core and the outer regions.

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.

We study the star formation history for a sample of 154 galaxies with stellar mass $10^{10}\lesssim M_{\ast}\lesssim 10^{12} M_{\odot}$ in the redshift range $0.7 < z < 0.9$. We do this using stellar population models combined with full-spectrum fitting of good quality spectra and high resolution photometry. For a subset of 68 galaxies ($M_{\ast}\gtrsim 10^{11} M_{\odot}$) we additionally construct dynamical models. These use an axisymmetric solution to the Jeans equations, which allows for velocity anisotropy, and adopts results from abundance matching techniques to account for the dark matter content. We find that: (i) The trends in star formation history observed in the local universe are already in place by $z\sim1$: the most massive galaxies are already passive, while lower mass ones have a more extended star formation histories, and the lowest mass galaxies are actively forming stars; (ii) we place an upper limit of a factor 1.5 to the size growth of the massive galaxy population; (iii) we present strong evidence for low dark matter fractions within $1R_{\rm e}$ (median of 9 per cent and 90th percentile of 21 per cent) for galaxies with $M_{\ast} \gtrsim 10^{11} M_{\odot}$ at these redshifts; and (iv) we confirm that these galaxies have, on average, a Salpeter normalisation of the stellar initial mass function.

NGC 4395 is a bulgeless spiral galaxy, harboring one of the nearest known type 1 Seyfert nuclei. Although there is no consensus on the mass of its central engine, several estimates suggest it to be one of the lightest massive black holes (MBHs) known. We present the first direct dynamical measurement of the mass of this MBH from a combination of two-dimensional gas kinematic data, obtained with the adaptive optics assisted near infrared integral field spectrograph Gemini/NIFS, and high-resolution multiband photometric data from Hubble Space Telescope's Wide Field Camera 3 (HST/WFC3). We use the photometric data to model the shape and stellar mass-to-light ratio (M/L) of the nuclear star cluster. From the Gemini/NIFS observations, we derive the kinematics of warm molecular hydrogen gas as traced by emission through the H$_2$ 1--0 S(1) transition. These kinematics show a clear rotational signal, with a position angle orthogonal to NGC 4395's radio jet. Our best fitting tilted ring models of the kinematics of the molecular hydrogen gas contain a black hole with mass $M=4_{-3}^{+8}\times 10^5$ M$_\odot$ (3$\sigma$ uncertainties) embedded in a nuclear star cluster of mass $M=2 \times 10^6$ M$_\odot$. Our black hole mass measurement is in excellent agreement with the reverberation mapping mass estimate of Peterson et al. (2005), but shows some tension with other mass measurement methods based on accretion signals.

We investigate the stellar population properties of a sample of 24 massive quenched galaxies at $1.25<z_\mathrm{spec}<2.09$ identified in the COSMOS field with our Subaru/MOIRCS near-IR spectroscopic observations. Tracing the stellar population properties as close to their major formation epoch as possible, we try to put constraints on the star formation history, post-quenching evolution, and possible progenitor star-forming populations for such massive quenched galaxies. By using a set of Lick absorption line indices on a rest-frame optical composite spectrum, the average age, metallicity [Z/H], and $\alpha$-to-iron element abundance ratio [$\alpha$/Fe] are derived as $\log(\mathrm{age}/\mathrm{Gyr})=0.04_{-0.08}^{+0.10}$, $\mathrm{[Z/H]}=0.24_{-0.14}^{+0.20}$, and $[\alpha/\mathrm{Fe}]=0.31_{-0.12}^{+0.12}$, respectively. If our sample of quenched galaxies at $\langle z \rangle = 1.6$ is evolved passively to $z=0$, their stellar population properties will align in excellent agreement with local counterparts at similar stellar velocity dispersions, which qualifies them as progenitors of local massive early-type galaxies. Redshift evolution of stellar population ages in quenched galaxies combined with low redshift measurements from the literature suggests a formation redshift of $z_\mathrm{f} \sim 2.3$ around which the bulk of stars in these galaxies have been formed. The measured [$\alpha$/Fe] value indicates a star formation timescale of $\lesssim 1$ Gyr, which can be translated into a specific star formation rate of $\simeq 1\,\mathrm{Gyr}^{-1}$ prior to quenching. Based on these findings, we discuss identifying possible progenitor star-forming galaxies at $z \simeq 2.3$. We identify normal star-forming galaxies, i.e, those on the star-forming main sequence, followed by a rapid quenching event, as likely precursors of the quenched galaxies at $\langle z \rangle = 1.6$ presented here.

Cappellari (2008) presented a flexible and efficient method to model the stellar kinematics of anisotropic axisymmetric and spherical stellar systems. The spherical formalism could be used to model the line-of-sight velocity second moments allowing for essentially arbitrary radial variation in the anisotropy and general luminous and total density profiles. Here we generalize the spherical formalism by providing the expressions for all three components of the projected second moments, including the two proper motion components. A reference implementation is now included in the public JAM package available at http://purl.org/cappellari/software

We study the total mass-density profile for a sample of 14 fast-rotator early-type galaxies (stellar masses $10.2<\log M_\ast/M_\odot<11.7$). We combine observations from the SLUGGS and Atlas3D surveys to map out the stellar kinematics in two-dimensions, out to a median radius for the sample of four half-light radii $R_e$ (or 10 kpc), and a maximum radius of 2.0-6.2 $R_e$ (or 4-21 kpc). We use axisymmetric dynamical models based on the Jeans equations, which allow for a spatially varying anisotropy, and employ quite general profiles for the dark halos, and in particular do not place any restriction on the profile slope. This is made possible by the availability of spatially extended two-dimensional kinematics. We find that our relatively simple models provide a remarkably good description of the observed kinematics. The resulting total density profiles are well described by a nearly-isothermal power law $\rho_{\rm tot}(r)\propto r^{-\gamma}$ from $R_e$/10 to at least 4$R_e$, the largest average deviation being 11%. The average logarithmic slope is $\langle\gamma\rangle=2.19\pm0.03$ with observed rms scatter of just $\sigma_\gamma=0.11$. This scatter out to large radii, where dark matter dominates, is as small as previously reported by lensing studies around $r\approx R_e/2$, where the stars dominate. Our bulge-halo conspiracy places much tighter constraints on galaxy formation models. It illustrates the power of two-dimensional stellar kinematics observations at large radii. It would now be important to test the generality of our results for different galaxy types and larger samples.

NGC 4258 is the galaxy with the most accurate (maser-based) determination for the mass of the supermassive black hole (SMBH) in its nucleus. In this work we present a two-dimensional mapping of the stellar kinematics in the inner 3.0 x 3.0 arcsec = 100 x 100 pc of NGC 4258 using adaptative-optics observations obtained with the Near-Infrared Integral Field Spectrograph of the GEMINI North telescope at a 0.11 arcsec (4 pc) angular resolution. The observations resolve the radius of influence of the SMBH, revealing an abrupt increase in the stellar velocity dispersion within 10 pc from the nucleus, consistent with the presence of a SMBH there. Assuming that the galaxy nucleus is in a steady state and that the velocity dispersion ellipsoid is aligned with a cylindrical coordinate system, we constructed a Jeans anisotropic dynamical model to fit the observed kinematics distribution. Our dynamical model assumes that the galaxy has axial symmetry and is constructed using the multi-gaussian expansion method to parametrize the observed surface brightness distribution. The Jeans dynamical model has three free parameters: the mass of the central SMBH, the mass-luminosity ratio of the galaxy and the anisotropy of the velocity distribution. We test two types of models: one with constant velocity anisotropy, and another with variable anisotropy. The model that best reproduces the observed kinematics was obtained considering that the galaxy has radially varying anisotropy, being the best-fitting parameters with 3$\sigma$ significance $M_\bullet=4.8^{+0.8}_{-0.9}\times 10^7\,{\rm M_\odot}$ and $\Gamma_k = 4.1^{+0.4}_{-0.5}$. This value for the mass of the SMBH is just 25 per cent larger than that of the maser determination and 50 per cent larger that a previous stellar dynamical determination obtained via Schwarzschild models.

MaNGA (Mapping Nearby Galaxies at Apache Point Observatory) is a 6-year SDSS-IV survey that will obtain resolved spectroscopy from 3600 $\AA$ to 10300 $\AA$ 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.

We present a high spatial resolution ($\approx 20$ pc) of $^{12}$CO($2-1$) observations of the lenticular galaxy NGC4526. We identify 103 resolved Giant Molecular Clouds (GMCs) and measure their properties: size $R$, velocity dispersion $\sigma_v$, and luminosity $L$. This is the first GMC catalog of an early-type galaxy. We find that the GMC population in NGC4526 is gravitationally bound, with a virial parameter $\alpha \sim 1$. The mass distribution, $dN/dM \propto M^{-2.39 \pm 0.03}$, is steeper than that for GMCs in the inner Milky Way, but comparable to that found in some late-type galaxies. We find no size-linewidth correlation for the NGC4526 clouds, in contradiction to the expectation from Larson's relation. In general, the GMCs in NGC4526 are more luminous, denser, and have a higher velocity dispersion than equal size GMCs in the Milky Way and other galaxies in the Local Group. These may be due to higher interstellar radiation field than in the Milky Way disk and weaker external pressure than in the Galactic center. In addition, a kinematic measurement of cloud rotation shows that the rotation is driven by the galactic shear. For the vast majority of the clouds, the rotational energy is less than the turbulent and gravitational energy, while the four innermost clouds are unbound and will likely be torn apart by the strong shear at the galactic center. We combine our data with the archival data of other galaxies to show that the surface density $\Sigma$ of GMCs is not approximately constant as previously believed, but varies by $\sim 3$ orders of magnitude. We also show that the size and velocity dispersion of GMC population across galaxies are related to the surface density, as expected from the gravitational and pressure equilibrium, i.e. $\sigma_v R^{-1/2} \propto \Sigma^{1/2}$.

We obtained Hubble Space Telescope/Wide Field Camera 3 imaging of a sample of ten of the nearest and brightest nuclear clusters residing in late-type spiral galaxies, in seven bands that span the near-ultraviolet to the near-infrared. Structural properties of the clusters were measured by fitting two-dimensional surface brightness profiles to the images using GALFIT. The clusters exhibit a wide range of structural properties. For six of the ten clusters in our sample, we find changes in the effective radius with wavelength, suggesting radially varying stellar populations. In four of the objects, the effective radius increases with wavelength, indicating the presence of a younger population which is more concentrated than the bulk of the stars in the cluster. However, we find a general decrease in effective radius with wavelength in two of the objects in our sample, which may indicate extended, circumnuclear star formation. We also find a general trend of increasing roundness of the clusters at longer wavelengths, as well as a correlation between the axis ratios of the NCs and their host galaxies. These observations indicate that blue disks aligned with the host galaxy plane are a common feature of nuclear clusters in late-type galaxies, but are difficult to detect in galaxies that are close to face-on. In color-color diagrams spanning the near-UV through the near-IR, most of the clusters lie far from single-burst evolutionary tracks, showing evidence for multi-age populations. Most of the clusters have integrated colors consistent with a mix of an old population (> 1 Gyr) and a young population (~100-300 Myr). The wide wavelength coverage of our data provides a sensitivity to populations with a mix of ages that would not be possible to achieve with imaging in optical bands only.

Even 10 billion years ago, the cores of the first galaxy clusters are often found to host a characteristic population of massive galaxies with already suppressed star formation. Here we search for distant cluster candidates at z~2 using massive passive galaxies as tracers. With a sample of ~40 spectroscopically confirmed passive galaxies at 1.3<z<2.1, we tune photometric redshifts of several thousands passive sources in the full 2 sq.deg. COSMOS field. This allows us to map their density in redshift slices, probing the large scale structure in the COSMOS field as traced by passive sources. We report here on the three strongest passive galaxy overdensities that we identify in the redshift range 1.5<z<2.5. While the actual nature of these concentrations is still to be confirmed, we discuss their identification procedure, and the arguments supporting them as candidate galaxy clusters (likely mid-10^13 M_sun range). Although this search approach is likely biased towards more evolved structures, it has the potential to select still rare, cluster-like environments close to their epoch of first appearance, enabling new investigations of the evolution of galaxies in the context of structure growth.

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.

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.

I review our understanding of classic dynamical scaling relations, relating luminosity, size and kinematics of early-type galaxies. Using unbiased determinations of galaxy mass profiles from stellar dynamical models, a simple picture has emerged in which scaling relations are driven by virial equilibrium, accompanied by a trend in the stellar mass-to-light ratio (M/L). This picture confirms the earliest insights. The trend is mainly due to the combined variation of age, metallicity and the stellar initial mass function (IMF). The systematic variations best correlate with the galaxy velocity dispersion, which traces the bulge mass fraction. This indicates a link between bulge growth and quenching of star formation. Dark matter is unimportant within the half-light radius, where the total mass profile is close to isothermal ($\rho\propto r^{-2}$).

We present the results on the stellar Initial Mass Function (IMF) normalisation of 68 massive ($M_\ast=10^{11}-10^{12}$$M_\odot$) Early-Type Galaxies (ETGs) at redshift of $\sim$1. This was achieved by deriving the stellar Mass-to-Light ratio (M/L) of the galaxies through axis-symmetric dynamical modelling and comparing it to the same derived via stellar population modelling through full spectrum fitting. The study also employs an Abundance Matching technique to account for the dark matter within the galaxies. The results demonstrate that massive ETGs at high redshifts on average have a Salpeter-like IMF normalisation, while providing observational evidence supporting previous predictions of low dark matter fraction in the inner regions ($\lesssim$ 1R$_{\rm e}$) of galaxies at higher redshift.

We present an investigation about the shape of the initial mass function (IMF) of early-type galaxies (ETGs), based on a joint lensing and dynamical analysis, and on stellar population synthesis models, for a sample of 55 lens ETGs identified by the Sloan Lens ACS (SLACS) Survey. We construct axisymmetric dynamical models based on the Jeans equations which allow for orbital anisotropy and include a dark matter halo. The models reproduce in detail the observed \textit{HST} photometry and are constrained by the total projected mass within the Einstein radius and the stellar velocity dispersion ($\sigma$) within the SDSS fibers. Comparing the dynamically-derived stellar mass-to-light ratios $(M_*/L)_{\rm dyn}$, obtained for an assumed halo slope $\rho_{\rm h}\propto r^{-1}$, to the stellar population ones $(M_*/L)_{\rm pop}$, derived from full-spectrum fitting and assuming a Salpeter IMF, we infer the mass normalization of the IMF. Our results confirm the previous analysis by the SLACS team that the mass normalization of the IMF of high $\sigma$ galaxies is consistent on average with a Salpeter slope. Our study allows for a fully consistent study of the trend between IMF and $\sigma$ for both the SLACS and \ATLAS samples, which explore quite different $\sigma$ ranges. The two samples are highly complementary, the first being essentially $\sigma$ selected, and the latter volume-limited and nearly mass selected. We find that the two samples merge smoothly into a single trend of the form $\log\alpha =(0.38\pm0.04)\times\log(\sigma_{\rm e}/200\,\mathrm{km~s}^{-1})+(-0.06\pm0.01)$, where $\alpha=(M_*/L)_{\rm dyn}/(M_*/L)_{\rm pop}$ and $\sigma_{\rm e}$ is the luminosity averaged $\sigma$ within one effective radius $R_{\rm e}$. This is consistent with a systematic variation of the IMF normalization from Kroupa to Salpeter in the interval $\sigma_{\rm e}\approx90-270\,\mathrm{km~s}^{-1}$.

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.

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.

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.

Galaxy masses play a fundamental role in our understanding of structure formation models. This review addresses the variety and reliability of mass estimators that pertain to stars, gas, and dark matter. The different sections on masses from stellar populations, dynamical masses of gas-rich and gas-poor galaxies, with some attention paid to our Milky Way, and masses from weak and strong lensing methods, all provide review material on galaxy masses in a self-consistent manner.

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.

We present integral field spectroscopy of 10 early-type galaxies in the nearby, low-mass, Fornax cluster, from which we derive spatially resolved stellar kinematics. Based on the morphologies of their stellar velocity maps we classify 2/10 galaxies as slow rotators, with the remaining 8 galaxies fast rotators. Supplementing our integral field observations with morphological and kinematic data from the literature, we analyse the `kinematic' type of all 30 galaxies in the Fornax cluster brighter than M_K = -21.5 mag (M_* ~ 6 x 10^9 M_sun). Our sample's slow rotator fraction within one virial radius is 7(^+4_-6) per cent. 13(^+8_-6} per cent of the early-type galaxies are slow rotators, consistent with the observed fraction in other galaxy aggregates. The fraction of slow rotators in Fornax varies with cluster-centric radius, rising to 16(^+11_-8) per cent of all kinematic types within the central 0.2 virial radii, from 0 per cent in the cluster outskirts. We find that, even in mass-matched samples of slow and fast rotators, slow rotators are found preferentially at higher projected environmental density than fast rotators. This demonstrates that dynamical friction alone cannot be responsible for the differing distributions of slow and fast rotators. For dynamical friction to play a significant role, slow rotators must reside in higher mass sub-halos than fast rotators and/or form in the centres of groups before being accreted on to the cluster.

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.

The stellar initial mass function (IMF) is a key parameter to study galaxy evolution. Here we measure the IMF mass normalization for a sample of 68 field galaxies in the redshift range 0.7 to 0.9 within the Extended Groth Strip. To do this we derive total (stellar + dark matter) mass-to-light [$(M/L)_{\rm dyn}$] using axisymmetric dynamical models. Within the region where we have kinematics (about one half-light radius), the models assume: (i) that mass-follows-light, implying negligible differences between the stellar and total density profiles; (ii) constant velocity anisotropy ($\beta_{\rm z}\equiv1-\sigma_z^2/\sigma_R^2=0.2$); (iii) that galaxies are seen at the average inclination for random orientations (i.e. $i=60^\circ$, where $i=90^\circ$ represents edge-on). The dynamical models are based on anisotropic Jeans equations, constrained by HST/ACS imaging and the central velocity dispersion of the galaxies, extracted from good-quality spectra taken by the DEEP2 survey. The population $(M/L)_{\rm pop}$ are derived from full-spectrum fitting of the same spectra with a grid of simple stellar population models. Recent dynamical modelling results by the $ATLAS^{3D}$ project and numerical simulations of galaxy evolution indicate that the dark matter fraction within the central regions of our galaxies should be small. This suggest that our derived total $(M/L)_{\rm dyn}$ should closely approximate the stellar $M/L$. Our comparison of $(M/L)_{\rm dyn}$ and $(M/L)_{\rm pop}$ then imply that for galaxies with stellar mass $M_\odot \geq 10^{11}$ $M_{\odot}$, the $average$ normalization of the IMF is consistent with a Salpeter slope, with a substantial scatter. This is similar to what is found within a similar mass range for nearby galaxies.

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.

We present a study of the cold gas contents of the Atlas3D early-type galaxies, in the context of their optical colours, near-UV colours, and H\beta\ absorption line strengths. Early-type (elliptical and lenticular) galaxies are not as gas-poor as previously thought, and at least 40% of local early-type galaxies are now known to contain molecular and/or atomic gas. This cold gas offers the opportunity to study recent galaxy evolution through the processes of cold gas acquisition, consumption (star formation), and removal. Molecular and atomic gas detection rates range from 10% to 34% in red sequence early-type galaxies, depending on how the red sequence is defined, and from 50% to 70% in blue early-type galaxies. Notably, massive red sequence early-type galaxies (stellar masses > 5e10 Msun, derived from dynamical models) are found to have HI masses up to M(HI)/Mstar ~ 0.06 and H_2 masses up to M(H$_2$)/Mstar ~ 0.01. Some 20% of all massive early-type galaxies may have retained atomic and/or molecular gas through their transition to the red sequence. However, kinematic and metallicity signatures of external gas accretion (either from satellite galaxies or the intergalactic medium) are also common, particularly at stellar masses <= 5e10 Msun, where such signatures are found in ~ 50% of H$_2$-rich early-type galaxies. Our data are thus consistent with a scenario in which fast rotator early-type galaxies are quenched former spiral galaxies which have undergone some bulge growth processes, and in addition, some of them also experience cold gas accretion which can initiate a period of modest star formation activity. We discuss implications for the interpretation of colour-magnitude diagrams.

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.

We present a detailed two-dimensional stellar dynamical analysis of a sample of 44 cosmological hydrodynamical simulations of individual central galaxies and their satellites. Kinematic maps of the stellar line-of-sight velocity, velocity dispersion, and higher-order Gauss-Hermite moments $h_3$ and $h_4$ are constructed for each central galaxy and for the most massive satellites. The amount of rotation is quantified using the $\lambda_{\mathrm{R}}$-parameter. The velocity, velocity dispersion, $h_3$, and $h_4$ fields of the simulated galaxies show a diversity similar to observed kinematic maps of early-type galaxies in the ATLAS$^{\rm{3D}}$ survey. This includes fast (regular), slow, and misaligned rotation, hot spheroids with embedded cold disk components as well as galaxies with counter-rotating cores or central depressions in the velocity dispersion. We link the present day kinematic properties to the individual cosmological formation histories of the galaxies. In general, major galaxy mergers have a significant influence on the rotation properties resulting in both a spin-down as well as a spin-up of the merger remnant. Lower mass galaxies with significant in-situ formation of stars, or with additional gas-rich major mergers - resulting in a spin-up - in their formation history, form elongated fast rotators with a clear anti-correlation of $h_3$ and $v/\sigma$. An additional formation path for fast rotators includes gas poor major mergers leading to a spin-up of the remnants. This formation path does not result in anti-correlated $h_3$ and $v/\sigma$. The galaxies most consistent with the rare class of non-rotating round early-type galaxies grow by gas-poor minor mergers alone. In general, more massive galaxies have less in-situ star formation since $z \sim 2$, rotate slower and have older stellar populations. (shortened)

The distribution of galaxies on the mass-size plane as a function of redshift or environment is a powerful test for galaxy formation models. Here we use integral-field stellar kinematics to interpret the variation of the mass-size distribution in two galaxy samples spanning extreme environmental densities. The samples are both identically and nearly mass-selected (stellar mass M*>6e9 Msun) and volume-limited. The first consists of nearby field galaxies from the Atlas3D parent sample. The second consists of galaxies in the Coma Cluster (Abell 1656), one of densest environments for which good resolved spectroscopy can be obtained. The mass-size distribution in the dense environment differs from the field one in two ways: (i) spiral galaxies are replaced by bulge-dominated disk-like fast-rotator early-type galaxies (ETGs), which follow the SAME mass-size relation and have the SAME mass distribution as in the field sample; (ii) the slow rotator ETGs are segregated in mass from the fast rotators, with their size increasing proportionally to their mass. A transition between the two processes appears around the stellar mass M_crit=2e11 Msun. We interpret this as evidence for bulge growth (outside-in evolution) and bulge-related environmental quenching dominating at low masses, with little influence from merging, while significant dry mergers (inside-out evolution) and halo-related quenching driving the mass and size growth at the high-mass end. The existence of these two processes naturally explains the diverse size evolution of galaxies of different masses and the separability of mass and environmental quenching.

Chronos is our response to ESA's call for white papers to define the science for the future L2, L3 missions. Chronos targets the formation and evolution of galaxies, by collecting the deepest NIR spectroscopic data, from the formation of the first galaxies at z~10 to the peak of formation activity at z~1-3. The strong emission from the atmospheric background makes this type of survey impossible from a ground-based observatory. The spectra of galaxies represent the equivalent of a DNA fingerprint, containing information about the past history of star formation and chemical enrichment. The proposed survey will allow us to dissect the formation process of galaxies including the timescales of quenching triggered by star formation or AGN activity, the effect of environment, the role of infall/outflow processes, or the connection between the galaxies and their underlying dark matter haloes. To provide these data, the mission requires a 2.5m space telescope optimised for a campaign of very deep NIR spectroscopy. A combination of a high multiplex and very long integration times will result in the deepest, largest, high-quality spectroscopic dataset of galaxies from z=1 to 12, spanning the history of the Universe, from 400 million to 6 billion years after the big bang, i.e. covering the most active half of cosmic history.

[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.

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]

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]

The masses of the supermassive black-holes found in galaxy bulges are correlated with a multitude of galaxy properties, leading to suggestions that galaxies and black-holes may evolve together. The number of reliably measured black-hole masses is small, and the number of methods for measuring them is limited, holding back attempts to understand this co-evolution. Directly measuring black-hole masses is currently possible with stellar kinematics (in early-type galaxies), ionised-gas kinematics (in some spiral and early-type galaxies) and in rare objects which have central maser emission. Here we report that by modelling the effect of a black-hole on the kinematics of molecular gas it is possible to fit interferometric observations of CO emission and thereby accurately estimate black hole masses. We study the dynamics of the gas in the early-type galaxy NGC4526, and obtain a best fit which requires the presence of a central dark-object of 4.5(+4.2-3.0)x10^8 Msun (3 sigma confidence limit). With next generation mm-interferometers (e.g. ALMA) these observations could be reproduced in galaxies out to 75 megaparsecs in less the 5 hours of observing time. The use of molecular gas as a kinematic tracer should thus allow one to estimate black-hole masses in hundreds of galaxies in the local universe, many more than accessible with current techniques.

For early-type galaxies, the ability to sustain a corona of hot, X-ray emitting gas could have played a key role in quenching their star-formation history. Yet, it is still unclear what drives the precise amount of hot gas around these galaxies. By combining photometric and spectroscopic measurements for the early-type galaxies observed during the Atlas3D integral-field survey with measurements of their X-ray luminosity based on X-ray data of both low and high spatial resolution we conclude that the hot-gas content of early-type galaxies can depend on their dynamical structure. Specifically, whereas slow rotators generally have X-ray halos with luminosity L_X,gas and temperature T values that are in line with what is expected if the hot-gas emission is sustained by the thermalisaton of the kinetic energy carried by the stellar-mass loss material, fast rotators tend to display L_X,gas values that fall consistently below the prediction of this model, with similar T values that do not scale with the stellar kinetic energy as observed in the case of slow rotators. Considering that fast rotators are likely to be intrinsically flatter than slow rotators, and that the few L_X,gas-deficient slow rotators also happen to be relatively flat, the observed L_X,gas deficiency in these objects would support the hypothesis whereby flatter galaxies have a harder time in retaining their hot gas. We discuss the implications that a different hot-gas content could have on the fate of both acquired and internally-produced gaseous material, considering in particular how the L_X,gas deficiency of fast rotators would make them more capable to recycle the stellar-mass loss material into new stars than slow rotators. This is consistent with the finding that molecular gas and young stars are detected only in fast rotators in the Atlas3D sample, and that fast rotators tend to dustier than slow rotators. [Abridged]

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]{\lambda}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.

[Abridged] We present a detailed study of the physical properties of the molecular gas in a sample of 18 molecular gas-rich early-type galaxies (ETGs) from the ATLAS$ 3D sample. Our goal is to better understand the star formation processes occurring in those galaxies, starting here with the dense star-forming gas. We use existing integrated $^{12}$CO(1-0, 2-1), $^{13}$CO(1-0, 2-1), HCN(1-0) and HCO$^{+}$(1-0) observations and present new $^{12}$CO(3-2) single-dish data. From these, we derive for the first time the average kinetic temperature, H$_{2}$ volume density and column density of the emitting gas, this using a non-LTE theoretical model. Since the CO lines trace different physical conditions than of those the HCN and HCO$^{+}$ lines, the two sets of lines are treated separately. We also compare for the first time the predicted CO spectral line energy distributions (SLEDs) and gas properties of our molecular gas-rich ETGs with those of a sample of nearby well-studied disc galaxies. The gas excitation conditions in 13 of our 18 ETGs appear analogous to those in the centre of the Milky Way. Such results have never been obtained before for ETGs and open a new window to explore further star-formation processes in the Universe. The conclusions drawn should nevertheless be considered carefully, as they are based on a limited number of observations and on a simple model. In the near future, with higher CO transition observations, it should be possible to better identify the various gas components present in ETGs, as well as more precisely determine their associated physical conditions. To achieve these goals, we show here from our theoretical study, that mid-J CO lines (such as the $^{12}$CO(6-5) line) are particularly useful.

We study the global efficiency of star formation in high resolution hydrodynamical simulations of gas discs embedded in isolated early-type and spiral galaxies. Despite using a universal local law to form stars in the simulations, we find that the early-type galaxies are offset from the spirals on the large-scale Kennicutt relation, and form stars 2 to 5 times less efficiently. This offset is in agreement with previous results on morphological quenching: gas discs are more stable against star formation when embedded in early-type galaxies due to the lower disc self-gravity and increased shear. As a result, these gas discs do not fragment into dense clumps and do not reach as high densities as in the spiral galaxies. Even if some molecular gas is present, the fraction of very dense gas (above 10^4 cm-3) is significantly reduced, which explains the overall lower star formation efficiency. We also analyse a sample of local early-type and spiral galaxies, measuring their CO and HI surface densities and their star formation rates as determined by their non-stellar 8um emission. As predicted by the simulations, we find that the early-type galaxies are offset from the Kennicutt relation compared to the spirals, with a twice lower efficiency. Finally, we validate our approach by performing a direct comparison between models and observations. We run a simulation designed to mimic the stellar and gaseous properties of NGC524, a lenticular galaxy, and find a gas disc structure and global star formation rate in good agreement with the observations. Morphological quenching thus seems to be a robust mechanism, and is also consistent with other observations of a reduced star formation efficiency in early-type galaxies in the COLD GASS survey. This lower efficiency of star formation is not enough to explain the formation of the whole Red Sequence, but can contribute to the reddening of some galaxies.

This is an addendum to the paper by Cappellari (2008, MNRAS, 390, 71), which presented a simple and efficient method to model the stellar kinematics of axisymmetric stellar systems. The technique reproduces well the integral-field kinematics of real galaxies. It allows for orbital anisotropy (three-integral distribution function), multiple kinematic components, supermassive black holes and dark matter. The paper described the derivation of the projected second moments and we provided a reference software implementation. However only the line-of-sight component was given in the paper. For completeness we provide here all the six projected second moments, including radial velocities and proper motions. We present a test against realistic N-body galaxy simulations.

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]

Since Edwin Hubble introduced his famous tuning fork diagram more than 70 years ago, spiral galaxies and early-type galaxies (ETGs) have been regarded as two distinct families. The spirals are characterized by the presence of disks of stars and gas in rapid rotation, while the early-types are gas poor and described as spheroidal systems, with less rotation and often non-axisymmetric shapes. The separation is physically relevant as it implies a distinct path of formation for the two classes of objects. I will give an overview of recent findings, from independent teams, that motivated a radical revision to Hubble's classic view of ETGs. These results imply a much closer link between spiral galaxies and ETGs than generally assumed.

NGC 1266 is a nearby field galaxy observed as part of the ATLAS3D survey (Cappellari et al. 2011). NGC 1266 has been shown to host a compact (< 200 pc) molecular disk and a mass-loaded molecular outflow driven by the AGN (Alatalo et al. 2011). Very Long Basline Array (VLBA) observations at 1.65 GHz revealed a compact (diameter < 1.2 pc), high bright- ness temperature continuum source most consistent with a low-level AGN origin. The VLBA continuum source is positioned at the center of the molecular disk and may be responsible for the expulsion of molecular gas in NGC 1266. Thus, the candidate AGN-driven molecular outflow in NGC 1266 supports the picture in which AGNs do play a significant role in the quenching of star formation and ultimately the evolution of the red sequence of galaxies.

We detail the rich molecular story of NGC 1266, its serendipitous discovery within the ATLAS3D survey (Cappellari et al. 2011) and how it plays host to an AGN-driven molecular outflow, potentially quenching all of its star formation (SF) within the next 100 Myr. While major mergers appear to play a role in instigating outflows in other systems, deep imaging of NGC 1266 as well as stellar kinematic observations from SAURON, have failed to provide evidence that NGC 1266 has recently been involved in a major interaction. The molecular gas and the instantaneous SF tracers indicate that the current sites of star formation are located in a hypercompact disk within 200 pc of the nucleus (Fig. 1; SF rate ~ 2 Msuns/yr). On the other hand, tracers of recent star formation, such as the H{\beta} absorption map from SAURON and stellar population analysis show that the young stars are distributed throughout a larger area of the galaxy than current star formation. As the AGN at the center of NGC 1266 continues to drive cold gas out of the galaxy, we expect star formation rates to decline as the star formation is ultimately quenched. Thus, NGC 1266 is in the midst of a key portion of its evolution and continued studies of this unique galaxy may help improve our understanding of how galaxies transition from the blue to the red sequence (Alatalo et al. 2011).

[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.

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)

The evolution of masses and sizes of passive (early-type) galaxies with redshift provides ideal constraints to galaxy formation models. These parameters can in principle be obtained for large galaxy samples from multi-band photometry alone. However the accuracy of photometric masses is limited by the non-universality of the IMF. Galaxy sizes can be biased at high redshift due to the inferior quality of the imaging data. Both problems can be avoided using galaxy dynamics, and in particular by measuring the galaxies stellar velocity dispersion. Here we provide an overview of the efforts in this direction.

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.

In order to try and understand its origins, we present high-quality long-slit spectral observations of the counter-rotating stellar discs in the strange S0 galaxy NGC 4550. We kinematically decompose the spectra into two counter-rotating stellar components (plus a gaseous component), in order to study both their kinematics and their populations. The derived kinematics largely confirm what was known previously about the stellar discs, but trace them to larger radii with smaller errors; the fitted gaseous component allows us to trace the hydrogen emission lines for the first time, which are found to follow the same rather strange kinematics previously seen in the [OIII] line. Analysis of the populations of the two separate stellar components shows that the secondary disc has a significantly younger mean age than the primary disc, consistent with later star formation from the associated gaseous material. In addition, the secondary disc is somewhat brighter, also consistent with such additional star formation. However, these measurements cannot be self-consistently modelled by a scenario in which extra stars have been added to initially-identical counter-rotating stellar discs, which rules out Evans & Collett's (1994) elegant "separatrix-crossing" model for the formation of such massive counter-rotating discs from a single galaxy, leaving some form of unusual gas accretion history as the most likely formation mechanism.

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.

As a demonstration of the capabilities of the new Oxford SWIFT integral field spectrograph, we present first observations for a set of 14 early-type galaxies in the core of the Coma cluster. Our data consist of I- and z-band spatially resolved spectroscopy obtained with the Oxford SWIFT spectrograph, combined with r-band photometry from the SDSS archive for 14 early- type galaxies. We derive spatially resolved kinematics for all objects from observations of the calcium triplet absorption features at \sim 8500 {AA} . Using this kinematic information we classify galaxies as either Fast Rotators or Slow Rotators. We compare the fraction of fast and slow rotators in our sample, representing the densest environment in the nearby Universe, to results from the ATLAS3D survey, finding the slow rotator fraction is \sim 50 per cent larger in the core of the Coma cluster than in the Virgo cluster or field, a 1.2 {\sigma} increase given our selection criteria. Comparing our sample to the Virgo cluster core only (which is 24 times less dense than the Coma core) we find no evidence of an increase in the slow rotator fraction. Combining measurements of the effective velocity dispersion {\sigma_e} with the photometric data we determine the Fundamental Plane for our sample of galaxies. We find the use of the average velocity dispersion within 1 effective radius, {\sigma_e}, reduces the residuals by 13 per cent with respect to comparable studies using central velocity dispersions, consistent with other recent integral field Fundamental Plane determinations.

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.

[Abridged] We present the results of new near-IR spectroscopic observations of passive galaxies at z>1.4 in a concentration of BzK-selected galaxies in the COSMOS field. The observations have been conducted with Subaru/MOIRCS, and have resulted in absorption lines and/or continuum detection for 18 out of 34 objects. This allows us to measure spectroscopic redshifts for a sample almost complete to K(AB)=21. COSMOS photometric redshifts are found in fair agreement overall with the spectroscopic redshifts, with a standard deviation of ~0.05; however, ~30% of objects have photometric redshifts systematically underestimated by up to ~25%. We show that these systematic offsets in photometric redshifts can be removed by using these objects as a training set. All galaxies fall in four distinct redshift spikes at z=1.43, 1.53, 1.67 and 1.82, with this latter one including 7 galaxies. SED fits to broad-band fluxes indicate stellar masses in the range of ~4-40x10^10Msun and that star formation was quenched ~1 Gyr before the cosmic epoch at which they are observed. The spectra of several individual galaxies have allowed us to measure their Hdelta_F and Dn4000 indices, which confirms their identification as passive galaxies, as does a composite spectrum resulting from the coaddition of 17 individual spectra. The effective radii of the galaxies have been measured on the HST/ACS F814W image, confirming the coexistence at these redshifts of passive galaxies which are substantially more compact than their local counterparts with others that follow the local size-stellar mass relation. For the galaxy with best S/N spectrum we were able to measure a velocity dispersion of 270+/-105 km/s, indicating that this galaxy lies closely on the virial relation given its stellar mass and effective radius.

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]

A combination of ground-based and spacecraft observations has uncovered two black holes of 10 billion solar masses in the nearby Universe. The finding sheds light on how these cosmic monsters co-evolve with galaxies.

The use of integral field units is now commonplace at all major observatories offering efficient means of obtaining spectral as well as imaging information at the same time. IFU instrument designs are complex and spectral images have typically highly condensed formats, therefore presenting challenges for the IFU data reduction pipelines. In the case of the VLT VIMOS-IFU, a fringe-like pattern affecting the spectra well into the optical and blue wavelength regime as well as artificial intensity variations, require additional reduction steps beyond standard pipeline processing. In this research note we propose an empirical method for the removal of the fringe-like pattern in the spectral domain and the intensity variations in the imaging domain. We also demonstrate the potential consequences for data analysis if the effects are not corrected. Here we use the example of deriving stellar velocity, velocity dispersion and absorption line-strength maps for early-type galaxies. We derive for each spectrum, reduced by the ESO standard VIMOS pipeline, a correction-spectrum by using the median of the eight surrounding spectra as a proxy for the unaffected, underlying spectrum. This method relies on the fact that our science targets (nearby ETGs) cover the complete FoV of the VIMOS-IFU with slowly varying spectral properties and that the exact shape of the fringe-like pattern is nearly independent and highly variable between neighboring spatial positions. We find that the proposed correction methods for the removal of the fringe-like pattern and the intensity variations in VIMOS-IFU data-cubes are suitable to allow for meaningful data analysis in our sample of nearby early-type galaxies. Since the method relies on the scientific target properties it is not suitable for general implementation in the pipeline software for VIMOS.

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{\beta} index and metallicity [Z/H]. The equivalent correlations for the Fe5015 index, {\alpha}-enhancement [{\alpha}/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.

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.

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.

Surveying eighteen 12CO-bright galaxies from the ATLAS3D early-type galaxy sample with the Institut de Radio Astronomie Millim\'etrique (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.

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.

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...

We propose a simple model for the origin of fast and slow rotator early-type galaxies (ETG) within the hierarchical $\Lambda$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 $\lambda_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 $\lambda_R$ into the slow rotator regime. Abridged...

[Abridged] We present ground-based MDM V-band and Spitzer/IRAC 3.6um-band photometric observations of the 72 representative galaxies of the SAURON Survey. In combination with the SAURON stellar velocity dispersion measured within an effective radius (se), this allows us to explore the location of our galaxies in the main scaling relations. We investigate the dependence of these relations on our recent kinematical classification of early-type galaxies (i.e. Slow/Fast Rotators) and the stellar populations. Slow Rotator and Fast Rotator E/S0 galaxies do not populate distinct locations in the scaling relations, although Slow Rotators display a smaller intrinsic scatter. Surprisingly, extremely young objects do not display the bluest (V-[3.6]) colours in our sample, as is usually the case in optical colours. This can be understood in the context of the large contribution of TP-AGB stars to the infrared, even for young populations, resulting in a very tight (V-[3.6]) - se relation that in turn allows us to define a strong correlation between metallicity and velocity dispersion. Many Sa galaxies appear to follow the Fundamental Plane defined by E/S0 galaxies. Galaxies that appear offset from the relations correspond mostly to objects with extremely young populations, with signs of on-going, extended star formation. We correct for this effect in the Fundamental Plane, by replacing luminosity with stellar mass using an estimate of the stellar mass-to-light ratio, so that all galaxies are part of a tight, single relation. The new estimated coefficients are consistent in both photometric bands and suggest that differences in stellar populations account for about half of the observed tilt with respect to the virial prediction. After these corrections, the Slow Rotator family shows almost no intrinsic scatter around the best-fit Fundamental Plane.

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)

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).

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)

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)

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.

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.

We present observations and dynamical models of the stellar nuclear clusters (NCs) at the centres of NGC 4244 and M33. We then compare these to an extensive set of simulations testing the importance of purely stellar dynamical mergers on the formation and growth of NCs. Mergers of star clusters are able to produce a wide variety of observed properties, including densities, structural scaling relations, shapes (including the presence of young discs) and even rapid rotation. Nonetheless, difficulties remain, most notably that the second order kinematic moment V_rms = (V^2 + sigma^2)^(1/2) of the models is too centrally peaked to match observations. This can be remedied by the merger of star clusters onto a pre-existing nuclear disc, but the line-of-sight velocity V is still more slowly rising than in NGC 4244. Our results therefore suggest that purely stellar dynamical mergers cannot form NCs, and that gas dissipation is a necessary ingredient for at least ~50% of a NC's mass. The negative vertical anisotropy found in NGC 4244 however requires at least 10% of the mass to be accreted as stars, since gas dissipation and in situ star formation leads to positive vertical anisotropy.

[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.

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.

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.

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.

The molecular gas content of local early-type galaxies is constrained and discussed in relation to their evolution. First, as part of the Atlas3D survey, we present the first complete, large (260 objects), volume-limited single-dish survey of CO in normal local early-type galaxies. We find a surprisingly high detection rate of 22%, independent of luminosity and at best weakly dependent on environment. Second, the extent of the molecular gas is constrained with CO synthesis imaging, and a variety of morphologies is revealed. The kinematics of the molecular gas and stars are often misaligned, implying an external gas origin in over a third of the systems, although this behaviour is drastically diffferent between field and cluster environments. Third, many objects appear to be in the process of forming regular kpc-size decoupled disks, and a star formation sequence can be sketched by piecing together multi-wavelength information on the molecular gas, current star formation, and young stars. Last, early-type galaxies do not seem to systematically obey all our usual prejudices regarding star formation, following the standard Schmidt-Kennicutt law but not the far infrared-radio correlation. This may suggest a greater diversity in star formation processes than observed in disk galaxies. Using multiple molecular tracers, we are thus starting to probe the physical conditions of the cold gas in early-types.

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.

The mass assembly of galaxies leaves various imprints on their surroundings, such as shells, streams and tidal tails. The frequency and properties of these fine structures depend on the mechanism driving the mass assembly: e.g. a monolithic collapse, rapid cold-gas accretion followed by violent disk instabilities, minor mergers or major dry / wet mergers. Therefore, by studying the outskirts of galaxies, one can learn about their main formation mechanism. I present here our on-going work to characterize the outskirts of Early-Type Galaxies (ETGs), which are powerful probes at low redshift of the hierarchical mass assembly of galaxies. This work relies on ultra-deep optical images obtained at CFHT with the wide-field of view MegaCam camera of field and cluster ETGs obtained as part of the Atlas-3D and NGVS projects. State of the art numerical simulations are used to interpret the data. The images reveal a wealth of unknown faint structures at levels as faint as 29 mag arcsec-2 in the g-band. Initial results for two galaxies are presented here.

[Abridged] We present two-dimensional line-of-sight stellar kinematics of the lens galaxy in the Einstein Cross, obtained with the GEMINI 8m telescope, using the GMOS integral-field spectrograph. The velocity map shows regular rotation up to ~100 km/s around the minor axis of the bulge, consistent with axisymmetry. The velocity dispersion map shows a weak gradient increasing towards a central (R<1") value of sigma_0=170+/-9 km/s. We deproject the observed surface brightness from HST imaging to obtain a realistic luminosity density of the lens galaxy, which in turn is used to build axisymmetric dynamical models that fit the observed kinematic maps. We also construct a gravitational lens model that accurately fits the positions and relative fluxes of the four quasar images. We find that the resulting luminous and total mass distribution are nearly identical around the Einstein radius R_E = 0.89", with a slope that is close to isothermal, but which becomes shallower towards the center if indeed mass follows light. The dynamical model fits to the observed kinematic maps result in a total mass-to-light ratio (M/L)_dyn=3.7+/-0.5 M_sun/L_sun,I (in the I-band). This is consistent with the Einstein mass M_E = 1.54 x 10^10 M_sun divided by the (projected) luminosity within R_E, which yields a total mass-to-light ratio of (M/L)_E=3.4 M_sun/L_sun,I, with an error of at most a few per cent. We estimate from stellar populations model fits to colors of the lens galaxy a stellar mass-to-light ratio (M/L)_* from 2.8 to 4.1 M_sun/L_sun,I. Although a constant dark matter fraction of 20 per cent is not excluded, dark matter may play no significant role in the bulge of this ~L* early-type spiral galaxy.

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.

We demonstrate that the comparison of Tully-Fisher relations (TFRs) derived from global HI line widths to TFRs derived from the circular velocity profiles of dynamical models (or stellar kinematic observations corrected for asymmetric drift) is vulnerable to systematic and uncertain biases introduced by the different measures of rotation used. We therefore argue that to constrain the relative locations of the TFRs of spiral and S0 galaxies, the same tracer and measure must be used for both samples. Using detailed near-infrared imaging and the circular velocities of axisymmetric Jeans models of 14 nearby edge-on Sa-Sb spirals and 14 nearby edge-on S0s drawn from a range of environments, we find that S0s lie on a TFR with the same slope as the spirals, but are on average 0.53+/-0.15 mag fainter at Ks-band at a given rotational velocity. This is a significantly smaller offset than that measured in earlier studies of the S0 TFR, which we attribute to our elimination of the bias associated with using different rotation measures and our use of earlier type spirals as a reference. Since our measurement of the offset avoids systematic biases, it should be preferred to previous estimates. A spiral stellar population in which star formation is truncated would take ~1 Gyr to fade by 0.53 mag at Ks-band. If S0s are the products of a simple truncation of star formation in spirals, then this finding is difficult to reconcile with the observed evolution of the spiral/S0 fraction with redshift. Recent star formation could explain the observed lack of fading in S0s, but the offset of the S0 TFR persists as a function of both stellar and dynamical mass. We show that the offset of the S0 TFR could therefore be explained by a systematic difference between the total mass distributions of S0s and spirals, in the sense that S0s need to be smaller or more concentrated than spirals.

We present an analysis of deep WSRT observations of the HI in 33 nearby early-type galaxies selected from a sample studied earlier at optical wavelengths with the SAURON integral-field spectrograph. The sample covers both field environments and the Virgo cluster. Our analysis shows that gas accretion plays a role in the evolution of field early-type galaxies, but less so for those in clusters. For detection limits of a few times 10^6 Msun, HI is detected in about 2/3 of the field galaxies, while <10% of the Virgo objects are detected. In about half of the detections, the HI forms a regularly rotating disc or ring. All HI discs have counterparts of ionised gas and inner HI discs are also detected in molecular gas. The cold ISM is dominated by molecular gas (M_H2/M_HI ~ 10). We conclude that accretion of HI is common for field early-type galaxies, but the amount of material involved is usually small. Cluster galaxies appear not to accrete HI. The few galaxies with a significant young sub-population all have inner gas discs, but for the remaining galaxies there is no trend between stellar population and HI. Some early-type galaxies are very gas rich, but only have an old population. The stellar populations of field galaxies are typically younger than those in Virgo. This is likely related to differences in accretion history. In about 50% of the galaxies we detect a central continuum source. In many objects this emission is from a low-luminosity AGN, in some it is consistent with the observed star formation. Galaxies with HI in the central regions are more likely detected in continuum. This is due to a higher probability for star formation to occur in such galaxies and not to HI-related AGN fuelling. (Abridged)

(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\'c 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.

We study resolution effects in numerical simulations of gas-rich and gas-poor major mergers, and show that the formation of slowly-rotating elliptical galaxies often requires a resolution that is beyond the present-day standards to be properly modelled. Our sample of equal-mass merger models encompasses various masses and spatial resolutions, ranging from about 200pc and 10^5 particles per component, typical of some recently published major merger simulations, to up to 32pc and 10^3 M_sun in simulations using 2.4 x 10^7 collisionless particles and 1.2 x 10^7 gas particles, among the highest resolutions reached so far for gas-rich major merger of massive disc galaxies. We find that the formation of fast-rotating early-type galaxies, that are flattened by a significant residual rotation, is overall correctly reproduced at all such resolutions. However, the formation of slow-rotating early-type galaxies, which have a low residual angular momentum and are supported mostly by anisotropic velocity dispersions, is strongly resolution-dependent. The evacuation of angular momentum from the main stellar body is largely missed at standard resolution, and systems that should be slow rotators are then found to be fast rotators. The effect is most important for gas-rich mergers, but is also witnessed in mergers with an absent or modest gas component. The effect is robust with respect to our initial conditions and interaction orbits, and originates in the physical treatment of the relaxation process during the coalescence of the galaxies. Our findings show that a high-enough resolution is required to accurately model the global properties of merger remnants and the evolution of their angular momentum. The role of gas-rich mergers of spiral galaxies in the formation of slow-rotating ellipticals may therefore have been underestimated.

We present observations of a very massive galaxy at z=1.82 which show that its morphology, size, velocity dispersion and stellar population properties that are fully consistent with those expected for passively evolving progenitors of today's giant ellipticals. These findings are based on a deep optical rest-frame spectrum obtained with the Multi-Object InfraRed Camera and Spectrograph (MOIRCS) on the Subaru telescope of a high-z passive galaxy candidate (pBzK) from the COSMOS field, for which we accurately measure its redshift of z=1.8230 and obtain an upper limit on its velocity dispersion sigma_star<326 km/s. By detailed stellar population modeling of both the galaxy broad-band SED and the rest-frame optical spectrum we derive a star-formation-weighted age and formation redshift of t_sf~1-2 Gyr and z_form~2.5-4, and a stellar mass of M_star~(3-4)x10^{11} M_sun. This is in agreement with a virial mass limit of M_vir<7x10^{11}M_sun, derived from the measured sigma_star value and stellar half-light radius, as well as with the dynamical mass limit based on the Jeans equations. In contrast with previously reported super-dense passive galaxies at z~2, the present galaxy at z=1.82 appears to have both size and velocity dispersion similar to early-type galaxies in the local Universe with similar stellar mass. This suggests that z~2 massive and passive galaxies may exhibit a wide range of properties, then possibly following quite different evolutionary histories from z~2 to z=0.

We examine the nuclear morphology, kinematics, and stellar populations in nearby S0 galaxy NGC 404 using a combination of adaptive optics assisted near-IR integral-field spectroscopy, optical spectroscopy, and HST imaging. These observations enable study of the NGC 404 nucleus at a level of detail possible only in the nearest galaxies. The surface brightness profile suggests the presence of three components, a bulge, a nuclear star cluster, and a central light excess within the cluster at radii <3 pc. These components have distinct kinematics with modest rotation seen in the nuclear star cluster and counter-rotation seen in the central excess. Molecular hydrogen emission traces a disk with rotation nearly orthogonal to that of the stars. The stellar populations of the three components are also distinct, with half of the mass of the nuclear star cluster having ages of ~1 Gyr (perhaps resulting from a galaxy merger), while the bulge is dominated by much older stars. Dynamical modeling of the stellar kinematics gives a total nuclear star cluster mass of 1.1x10^7 Msol. Dynamical detection of a possible intermediate mass black hole is hindered by uncertainties in the central stellar mass profile. Assuming a constant mass-to-light ratio, the stellar dynamical modeling suggests a black hole mass of <1x10^5 Msol, while the molecular hydrogen gas kinematics are best fit by a black hole with mass of 4.5x10^5 Msol. Unresolved and possibly variable dust emission in the near-infrared and AGN-like molecular hydrogen emission line ratios do suggest the presence of an accreting black hole in this nearby LINER galaxy.

We summarize the recent results of our survey of the nearest nuclear star clusters. The purpose of the survey is to understand nuclear star cluster formation mechanisms and constrain the presence of black holes using adaptive optics assisted integral field spectroscopy, optical spectroscopy, and HST imaging in 13 galaxies within 5 Mpc. We discuss the formation history of the nuclear star cluster and possible detection of an intermediate mass BH in NGC 404, the nearest S0 galaxy.

We show that high quality laser guide star (LGS) adaptive optics (AO) observations of nearby early-type galaxies are possible when the tip-tilt correction is done by guiding on nuclei while the focus compensation due to the changing distance to the sodium layer is made 'open loop'. We achieve corrections such that 40% of flux comes from R<0.2 arcsec. To measure a black hole mass (Mbh) one needs integral field observations of both high spatial resolution and large field of view. With these data it is possible to determine the lower limit to Mbh even if the spatial resolution of the observations are up to a few times larger than the sphere of influence of the black hole.

We investigate the accuracy of mass determinations M_BH of supermassive black holes in galaxies using dynamical models of the stellar kinematics. We compare 10 of our M_BH measurements, using integral-field OASIS kinematics, to published values. For a sample of 25 galaxies we confront our new M_BH derived using two modeling methods on the same OASIS data.

We construct mass models of 28 S0-Sb galaxies. The models have an axisymmetric stellar component and a NFW dark halo and are constrained by observed Ks-band photometry and stellar kinematics. The median dark halo virial mass is 10^12.8 Msun, and the median dark/total mass fraction is 20% within a sphere of radius r_1/2, the intrinsic half-light radius, and 50% within R_25. We compare the Tully-Fisher relations of the spirals and S0s in the sample and find that S0s are 0.5 mag fainter than spirals at Ks-band and 0.2 dex less massive for a given rotational velocity. We use this result to rule out scenarios in which spirals are transformed into S0s by processes which truncate star formation without affecting galaxy dynamics or structure, and raise the possibility of a break in homology between spirals and S0s.

We combine SAURON integral field data of a representative sample of local early-type, red sequence galaxies with Spitzer/IRAC imaging in order to investigate the presence of trace star formation in these systems. With the Spitzer data, we identify galaxies hosting low-level star formation, as traced by PAH emission, with measured star formation rates that compare well to those estimated from other tracers. This star formation proceeds according to established scaling relations with molecular gas content, in surface density regimes characteristic of disk galaxies and circumnuclear starbursts. We find that star formation in early-type galaxies happens exclusively in fast-rotating systems and occurs in two distinct modes. In the first, star formation is a diffuse process, corresponding to widespread young stellar populations and high molecular gas content. The equal presence of co- and counter-rotating components in these systems strongly implies an external origin for the star-forming gas, and we argue that these star formation events may be the final stages of (mostly minor) mergers that build up the bulges of red sequence lenticulars. In the second mode of star formation, the process is concentrated into well-defined disk or ring morphologies, outside of which the host galaxies exhibit uniformly evolved stellar populations. This implies that these star formation events represent rejuvenations within previously quiescent stellar systems. Evidence for earlier star formation events similar to these in all fast rotating early-type galaxies suggests that this mode of star formation may be common to all such galaxies, with a duty cycle of roughly 1/10, and likely contributes to the embedded, co-rotating inner stellar disks ubiquitous in this population.

We present a preliminary analysis of the HI properties of early-type galaxies in the ATLAS3D sample. Using WSRT data for ~100 galaxies outside the Virgo cluster and data from the Alfalfa project for galaxies inside Virgo, we discuss the dependence of HI properties on environment. We detect HI in about half of the galaxies outside Virgo. For these systems, the HI morphology and kinematics change as a function of environment, going from regular, rotating systems around isolated galaxies to progressively more disturbed structures for galaxies with neighbours or in groups. In denser environment, inside Virgo, nearly none of the galaxies contains HI. We discuss future work in this field which will be enabled by next-generation, pre-SKA radio instruments. We present a simulated Apertif HI observation of an ATLAS3D early-type galaxy, showing how its appearance and detection level vary as a function of redshift.

We review the concepts of the Voronoi binning technique (Cappellari & Copin 2003), which optimally solves the problem of preserving the maximum spatial resolution of general two-dimensional data, given a constraint on the minimum signal-to-noise ratio (S/N). This is achieved by partitioning the data in an adaptive fashion using a Voronoi tessellation with nearly hexagonal lattice. We review astrophysical applications of the method to X-ray data, integral-field spectroscopy, Fabry-Perot interferometry, N-body simulations, standard images and other regularly or irregularly sampled data. Voronoi binning, unlike adaptive smoothing, produces maps where the noise in the data can be visually assessed and spurious artifacts can be recognized. The method can be used to bin data according to any general criterion and not just S/N. It can be applied to higher dimensions and it can be used to generate optimal adaptive meshes for numerical simulations.

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\AA) and LINER-like [OIII]/Hb values, the nebular emission is truly powered by an AGN.

(Abridged) We present mass models of a sample of 14 spiral and 14 S0 galaxies that constrain their stellar and dark matter content. For each galaxy we derive the stellar mass distribution from near-infrared photometry under the assumptions of axisymmetry and a constant Ks-band stellar mass-to-light ratio, (M/L)_Ks. To this we add a dark halo assumed to follow a spherically symmetric NFW profile and a correlation between concentration and dark mass within the virial radius, M_DM. We solve the Jeans equations for the corresponding potential under the assumption of constant anisotropy in the meridional plane, beta_z. By comparing the predicted second velocity moment to observed long-slit stellar kinematics, we determine the three best-fitting parameters of the model: (M/L)_Ks, M_DM and beta_z. These simple axisymmetric Jeans models are able to accurately reproduce the wide range of observed stellar kinematics, which typically extend to ~2-3 Re or, equivalently, ~0.5-1 R_25. We find a median stellar mass-to-light ratio at Ks-band of 1.09 (solar units) with an rms scatter of 0.31. We present preliminary comparisons between this large sample of dynamically determined stellar mass-to-light ratios and the predictions of stellar population models. The stellar population models predict slightly lower mass-to-light ratios than we measure. The mass models contain a median of 15 per cent dark matter by mass within an effective radius Re, and 49 per cent within the optical radius R_25. Dark and stellar matter contribute equally to the mass within a sphere of radius 4.1 Re or 1.0 R_25. There is no evidence of any significant difference in the dark matter content of the spirals and S0s in our sample.

We measured stellar velocity dispersions sigma and dynamical masses of 9 massive (M~10^11 Msun) early-type galaxies (ETG) from the GMASS sample at redshift 1.4<z<2.0. The sigma are based on individual spectra for two galaxies at z~1.4 and on a stacked spectrum for 7 galaxies with 1.6<z<2.0, with 202-h of exposure at the ESO Very Large Telescope. We constructed detailed axisymmetric dynamical models for the objects, based on the Jeans equations, taking the observed surface brightness (from deep HST/ACS observations), PSF and slit effects into account. Our dynamical masses M_Jeans agree within ~30% with virial estimates M_vir=5*Re*sigma^2/G, although the latter tend to be smaller. This suggests that sizes are not underestimated by more than a similar fraction. Our M_Jeans also agrees within a factor <2 with the M_pop previously derived using stellar population models and 11 bands photometry. This confirms that the galaxies are intrinsically massive. The inferred mass-to-light ratios M/L_U in the very age-sensitive rest frame U-band are consistent with passive evolution in the past ~1 Gyr (formation redshift z_f~3). A 'bottom-light' stellar Initial Mass Function (IMF) appears to be required to ensure close agreement between M_Jeans and M_pop at z~2, as it does at z~0. The GMASS ETGs are on average more dense than their local counterpart. However a few percent of local ETGs of similar dynamical masses also have comparable sigma and mass surface density Sigma_50 inside Re.

We study the origin of the scaling relations of early-type galaxies (ETGs) by constructing detailed models of the stellar dynamics for the K-band selected, volume-limited ATLAS3D sample of 263 nearby ETGs, spanning a large range of masses and stellar velocity dispersions (60 < sigma < 350 km/s).

[abridged] We present observations of NGC524 and NGC2549 with LGS AO obtained at GEMINI North telescope using the NIFS IFU in the K band. The purpose of these observations, together with previously obtained observations with the SAURON IFU, is to determine the masses (Mbh) of the supermassive black holes (SMBH). The targeted galaxies were chosen to have central light profiles showing a core (NGC524) and a cusp (NGC2549), to probe the feasibility of using the galaxy centre as the NGS required for LGS AO. We employ an innovative `open loop' technique. The data have spatial resolution of 0.23" and 0.17" FWHM, showing that high quality LGS AO observations of these objects are possible. We construct axisymmetric three-integral dynamical models which are constrained with both the NIFS and SAURON data. The best fitting models yield Mbh=(8.3 +2.7 -1.3) x 10^8 Msun for NGC524 and Mbh=(1.4 +0.2 -1.3) x 10^7 Msun for NGC2549 (all errors are at the 3 sigma CL). We demonstrate that the wide-field SAURON data play a crucial role in the M/L determination increasing the accuracy of M/L by a factor of at least 5, and constraining the upper limits on Mbh. The NIFS data are crucial in constraining the lower limits of Mbh and in combination with the large scale data reducing the uncertainty by a factor of 2 or more. We find that the orbital structure of NGC524 shows significant tangential anisotropy, while at larger radii both galaxies are consistent with having almost perfectly oblate velocity ellipsoids. Tangential anisotropy in NGC524 coincides with the size of SMBH sphere of influence and the core region in the light profile. We test the accuracy to which Mbh can be measured using seeings obtained from typical LGS AO observations, and conclude that for a typical conditions and Mbh the expected uncertainty is of the order of 50%.

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.

We present the results of an investigation of the local escape velocity (Vesc) - line strength index relationship for 48 early type galaxies from the SAURON sample, the first such study based on a large sample of galaxies with both detailed integral field observations and extensive dynamical modelling. Values of Vesc are computed using Multi Gaussian Expansion (MGE) photometric fitting and axisymmetric, anisotropic Jeans' dynamical modelling simultaneously on HST and ground-based images. We determine line strengths and escape velocities at multiple radii within each galaxy, allowing an investigation of the correlation within individual galaxies as well as amongst galaxies. We find a tight correlation between Vesc and the line-strength indices. For Mgb we find that this correlation exists not only between different galaxies but also inside individual galaxies - it is both a local and global correlation. The relation within individual galaxies has the same slope and offset as the global relation to a good level of agreement, though there is significant intrinsic scatter in the local gradients. We transform our line strength index measurements to the single stellar population (SSP) equivalent age (t), metallicity ([Z/H]) and enhancement ([$\alpha$/Fe]) and carry out a principal component analysis of our SSP and Vesc data. We find that in this four-dimensional parameter space the galaxies in our sample are to a good approximation confined to a plane, given by $\log \mathrm({V}_{\mathrm{esc}}/500\mathrm{km/s}) = 0.85 \mathrm{[Z/H]} + 0.43 \log (\mathrm{t}/\mathrm{Gyrs})$ - 0.20. It is surprising that it is a combination of age and metallicity that is conserved; this may indicate a 'conspiracy' between age and metallicity or a weakness in the SSP models.

We use the integral-field spectrograph SAURON to measure the stellar line-of-sight velocity distribution and absorption line strengths out to four effective radii (Re) in the early-type galaxies NGC 3379 and NGC 821. With our newly developed observing technique we can now probe these faint regions in galaxies that were previously not accessible with traditional long-slit spectroscopy. We make optimal use of the large field-of-view and high throughput of the spectrograph: by adding the signal of all ~1400 lenslets into one spectrum, we obtain sufficient signal-to-noise in a few hours of observing time to reliably measure the absorption line kinematics and line strengths out to large radius. We find that the line strength gradients previously observed within 1 Re remain constant out to at least 4 Re, which puts constraints on the merger histories of these galaxies. The stellar halo populations are old and metal-poor. By constructing orbit-based Schwarzschild dynamical models we find that dark matter is necessary to explain the observed kinematics in NGC 3379 and NGC 821, with 30 - 50 per cent of the total matter being dark within 4 Re. The radial anisotropy in our best-fit halo models is less than in our models without halo, due to differences in orbital structure. The halo also has an effect on the Mgb - Vesc relation: its slope is steeper when a dark matter halo is added to the model.

We measure the Tully-Fisher relations of 14 lenticular galaxies (S0s) and 14 spirals. We use two measures of rotational velocity. One is derived directly from observed spatially-resolved stellar kinematics and the other from the circular velocities of mass models that include a dark halo and whose parameters are constrained by detailed kinematic modelling. Contrary to the naive expectations of theories of S0 formation, we find no significant difference between the Tully-Fisher relations of the two samples when plotted as functions of both brightness and stellar mass.

We examine a sample of 30 edge-on spiral and S0 galaxies that have boxy and peanut-shaped bulges. We compute model stellar kinematics by solving the Jeans equations for axisymmetric mass distributions derived from K-band images. These simple models have only one free parameter: the dynamical mass-to-light ratio, which we assume is independent of radius. Given the simplicity of the modelling procedure, the model second velocity moments are strikingly good fits to the observed stellar kinematics within the extent of our kinematic data, which typically reach ~ 0.5-1 R25 (where R25 is the optical radius), or equivalently 2-3 Re (where Re is the effective or half-light radius). We therefore find no evidence for a dominant dark matter component within the optical disk of spiral galaxies. This is equally true of the S0s in our sample, which significantly extends previous observational constraints on dark matter in these galaxies. The predicted kinematics do deviate slightly but systematically from the observations in the bulge region of most galaxies, but we argue that this is consistent with the claim that boxy and peanut-shaped bulges are bars viewed edge-on.

We present a simple and efficient anisotropic generalization of the semi-isotropic (two-integral) axisymmetric Jeans formalism which is used to model the stellar kinematics of galaxies. The following is assumed: (i) a constant mass-to-light ratio M/L and (ii) a velocity ellipsoid that is aligned with cylindrical coordinates (R,z) and characterized by the classic anisotropy parameter beta_z=1-sigma_z^2/sigma_R^2. Our simple models are fit to SAURON integral-field observations of the stellar kinematics for a set of fast-rotator early-type galaxies. With only two free parameters (beta_z and the inclination) the models generally provide remarkably good descriptions of the shape of the first (V) and second (V_rms=sqrt{V^2+sigma^2}) velocity moments, once a detailed description of the surface brightness is given. This is consistent with previous findings on the simple dynamical structure of these objects. With the observationally-motivated assumption that beta_z>0, the method is able to recover the inclination. The technique can be used to determine the dynamical mass-to-light ratios and angular momenta of early-type fast-rotators and spiral galaxies, especially when the quality of the data does not justify more sophisticated modeling approaches. This formalism allows for the inclusion of dark matter, supermassive black holes, spatially varying anisotropy, and multiple kinematic components.

We present a determination of the mass of the supermassive black hole (BH) and the nuclear stellar orbital distribution of the elliptical galaxy Centaurus A (NGC5128) using high-resolution integral-field observations of the stellar kinematics. The observations were obtained with SINFONI at the ESO Very Large Telescope in the near-infrared (K-band), using adaptive optics to correct for the blurring effect of the earth atmosphere. The data have a spatial resolution of 0.17" FWHM and high S/N>80 per spectral pixel so that the shape of the stellar line-of-sight velocity-distribution can be reliably extracted. We detect clear low-level stellar rotation, which is counter-rotating with respect to the gas. We fit axisymmetric three-integral dynamical models to the data to determine the best fitting values for the BH mass M_BH=(5.5+/-3.0)*10^7 Msun (3sigma errors) and (M/L)_K=(0.65+/-0.15) in solar units. These values are in excellent agreement with previous determinations from the gas kinematics, and in particular with our own published values, extracted from the same data. This provides one of the cleanest gas versus stars comparisons of BH determination, due to the use of integral-field data for both dynamical tracers and due to a very well resolved BH sphere of influence R_BH~0.70". We derive an accurate profile of the orbital anisotropy and we carefully test its reliability using spherical Jeans models with radially varying anisotropy. We find an increase in the tangential anisotropy close to the BH, but the spatial extent of this effect seems restricted to the size of R_BH instead of that R_b~3.9" of the core in the surface brightness profile, contrary to detailed predictions of current simulations of the binary BH scouring mechanism. More realistic simulations would be required to draw conclusions from this observation.

We describe a new technique to measure stellar kinematics and line-strengths at large radii in nearby galaxies. Using the integral-field spectrograph SAURON as a 'photon-collector', we obtain spectra out to four effective radii (Re) in the early-type galaxy NGC 3379. By fitting orbit-based models to the extracted stellar velocity profile, we find that ~40% of the total mass within 5 Re is dark. The measured absorption line-strengths reveal a radial gradient with constant slope out to 4 Re.

We use two-dimensional kinematic maps of simulated binary disc mergers to investigate the $\lambda_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 $\lambda_R > 0.1$ and slow rotators $\lambda_R < 0.1$. Tests on merger remnants reveal that $\lambda_R$ is a robust indicator of the true angular momentum content in elliptical galaxies. We find the same range of $\lambda_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 $\lambda_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 $\lambda_R$ as high as our progenitor disc galaxies.

[Abridged] We analysed two-dimensional maps of 48 early-type galaxies obtained with the SAURON and OASIS integral-field spectrographs using kinemetry, a generalisation of surface photometry to the higher order moments of the line-of-sight velocity distribution (LOSVD). In the SAURON sample, we find that 31% of early-type galaxies are single component systems. 91% of the multi-components systems have two kinematic subcomponents, the rest having three. In addition, 29% of galaxies have kinematically decoupled components, nuclear components with significant kinematic twists. We find that the velocity maps of fast rotators closely resemble those of inclined disks, except in the transition regions between kinematic subcomponents. In terms of E/S0 classification, this means that 74% of Es and 92% of S0s have components with disk-like kinematics. For the majority of fast rotators, the kinematic axial ratios are equal to or less than their photometric axial ratios, contrary to what is predicted with isotropic Jeans models viewed at different inclinations. The position angles of fast rotators are constant, while they vary abruptly in slow rotators. Velocity dispersion maps of face-on galaxies have shapes similar to the distribution of light. We constructed local (bin-by-bin) h3 - V/sigma and h4 - V/sigma diagrams from SAURON observations. We confirm the classical anti-correlation of h3 and V\sigma, but we also find that h3 is almost zero in some objects or even weakly correlated with V/sigma. The distribution of h4 for fast and slow rotators is mildly positive on average. The difference between slow and fast rotators is traceable throughout all moments of the LOSVD, with evidence for different intrinsic shapes and orbital contents and, hence, likely different evolutionary paths.

We present interferometric observations resolving the CO emission in the four gas-rich lenticular galaxies NGC 3032, NGC 4150, NGC 4459, and NGC 4526, and we compare the CO distribution and kinematics to those of the stars and ionized gas. Counterrotation documents an external origin for the gas in at least one case (NGC 3032), and the comparisons to stellar and ionized gas substructures in all four galaxies offer insights into their formation histories. The molecular gas is found in kpc-scale disks with mostly regular kinematics and average surface densities of 100 to 200 \msunsqpc. The disks are well aligned with the stellar photometric and kinematic axes. In the two more luminous Virgo Cluster members NGC 4459 and NGC 4526 the molecular gas shows excellent agreement with circular velocities derived independently from detailed modeling of stellar kinematic data. There are also two puzzling instances of disagreements between stellar kinematics and gas kinematics on sub-kpc scales. In the inner arcseconds of NGC 3032 the CO velocities are significantly lower than the inferred circular velocities, and the reasons may possibly be related to the external origin of the gas but are not well understood. In addition, the very young population of stars in the core of NGC 4150 appears to have the opposite sense of rotation from the molecular gas.

We present a flexible and efficient method to construct triaxial dynamical models of galaxies with a central black hole, using Schwarzschild's orbital superposition approach. Our method is general and can deal with realistic luminosity distributions, which project to surface brightness distributions that may show position angle twists and ellipticity variations. The models are fit to measurements of the full line-of-sight velocity distribution (wherever available). We verify that our method is able to reproduce theoretical predictions of a three-integral triaxial Abel model. In a companion paper (van de Ven, de Zeeuw & van den Bosch), we demonstrate that the method recovers the phase-space distribution function. We apply our method to two-dimensional observations of the E3 galaxy NGC 4365, obtained with the integral-field spectrograph SAURON, and study its internal structure, showing that the observed kinematically decoupled core is not physically distinct from the main body and the inner region is close to oblate axisymmetric.

We discuss some recent integral field spectroscopy using the SAURON instrument of a sample consisting of 24 early-type spirals, part of the SAURON Survey, and 18 late-type spirals. Using 2-dimensional maps of their stellar radial velocity, velocity dispersion, and absorption line strength, it is now much easier to understand the nature of nearby galactic bulges. We discuss a few highlights of this work, and point out some new ideas about the formation of galactic bulges.

Using the unique dataset obtained within the course of the SAURON project, a radically new view of the structure, dynamics and stellar populations of early-type galaxies has emerged. We show that galaxies come in two broad flavours (slow and fast rotators), depending on whether or not they exhibit clear large-scale rotation, as indicated via a robust measure of the specific angular momentum of baryons. This property is also linked with other physical characteristics of early-type galaxies, such as: the presence of dynamically decoupled cores, orbital structure and anisotropy, stellar populations and dark matter content. I here report on the observed link between this baryonic angular momentum and a mass sequence, and how this uniquely relates to the building of the red sequence via dissipative/dissipationless mergers and secular evolution.

Supermassive black holes are a key element in our understanding of how galaxies form. Most of the progress in this very active field of research is based on just ~30 determinations of black hole mass, accumulated over the past decade. We illustrate how integral-field spectroscopy, and in particular our OASIS modeling effort, can help improve the current situation.

We present two-dimensional gas-kinematic maps of the central region in Centaurus A. The adaptive optics (AO) assisted SINFONI data from the VLT have a resolution of 0.12" in K-band. The ionized gas species (Br_gamma, [FeII], [SiVI]) show a rotational pattern that is increasingly overlaid by non-rotational motion for higher excitation lines in direction of Cen A's radio jet. The emission lines of molecular hydrogen (H_2) show regular rotation and no distortion due to the jet. The molecular gas seems to be well settled in the gravitational potential of the stars and the central supermassive black hole and we thus use it as a tracer to model the mass in the central +/-1.5". These are the first AO integral-field observations on the nucleus of Cen A, enabling us to study the regularity of the rotation around the black hole, well inside the radius of influence, and to determine the inclination angle of the gas disk in a robust way. The gas kinematics are best modeled through a tilted-ring model that describes the warped gas disk; its mean inclination angle is ~34deg and the mean position angle of the major axis is ~155deg. The best-fit black hole mass is M_BH~4.5x10^7 Msolar, based on a "kinematically hot" disk model where the velocity dispersion is included through the Jeans equation. This black hole mass estimate is somewhat lower than, but consistent with the mass values previously derived from ionized gas kinematics. It is also consistent with the stellar dynamical measurement from the same AO observations, which we present in a separate paper. It brings Cen A in agreement with the M_BH-sigma relation.

Motivated by recent progress in the study of early-type galaxies owing to technological advances, the launch of new space telescopes and large ground-based surveys, we attempt a short review of our current understanding of the recent star-formation activity in such intriguing galactic systems.

In this proceeding we make use of the two-dimensional stellar kinematics of a representative sample of E and S0 galaxies obtained with the SAURON integral-field spectrograph to reveal that early-type galaxies appear in two broad flavours, depending on whether they exhibit clear large-scale rotation or not. We measure the level of rotation via a new parameter LambdaR and use it as a basis for a new kinematic classification that separates early-type galaxies into slow and fast rotators. With the aid of broad-band imaging we will reinforce this finding by comparing our kinematic results to the photometric properties of these two classes.

In this proceeding we look at the relationship between the photometric nuclear properties of early-type galaxies from Hubble Space Telescope imaging and their overall kinematics as observed with the SAURON integral-field spectrograph. We compare the inner slope of their photometric profiles and the Slow/Fast rotator classes, defined by the amplitude of a newly defined LambdaR parameter, to show that slow rotators tend to be more massive systems and display shallower inner profiles and fast rotators steper ones. It is important to remark, however, that there is not a one-to-one relationship between the two photometric and kinematic groups.

2D stellar kinematics of 48 representative E and S0 galaxies obtained with the SAURON IFS reveal that early-type galaxies appear in two broad flavours, depending on whether they exhibit clear large-scale rotation or not. We define a new parameter LambdaR which involves luminosity weighted averages over the full 2D kinematic field, as a proxy to quantify the observed projected stellar angular momentum per unit mass. We use it as a basis for a new kinematic classification: early-type galaxies are separated into slow and fast rotators (SRs, FRs), depending on whether they have LambdaR values within their effective radius Re below or above 0.1, respectively. SRs and FRs are shown to be physically distinct classes of galaxies, a result which cannot simply be the consequence of a biased viewing angle. FRs tend to be relatively low luminosity galaxies. SRs tend to be brighter and more massive galaxies, but are still spread over a wide range of absolute magnitude. 3 slow rotators of our sample, among the most massive ones, are consistent with zero rotation. Remarkably, all other SRs contain a large kpc-scale KDC. All FRs show well aligned photometric and kinemetric axes, and small velocity twists, in contrast with most SRs which exhibit significant misalignments and velocity twists. In a companion paper (Paper X), we also show that FRs and SRs are distinct classes in terms of their orbital distribution. We suggest that gas is a key ingredient in the formation and evolution of FRs, and that the slowest rotators are the extreme evolutionary end point reached deep in gravitational potential wells where dissipationless mergers had a major role in the evolution, and for which most of the baryonic angular momentum was expelled outwards. (abridged)

We present absorption line-strength maps for a sample of 18 Sb-Sd galaxies observed using the integral-field spectrograph SAURON. The SAURON spectral range allows the measurement of the Lick/IDS indices Hbeta, Fe5015 and Mgb, which can be used to estimate the stellar population parameters. We present here the two-dimensional line-strength maps for each galaxy. From the maps, we learn that late-type spiral galaxies tend to have high Hbeta and low Fe5015 and Mgb values, and that the Hbeta index has often a positive gradient over the field, while the metal indices peak in the central region. We investigate the relations between the central line-strength indices and their correlations with morphological type and central velocity dispersion, and compare the observed behaviour with that for ellipticals, lenticulars and early-type spirals from the SAURON survey. We find that our galaxies lie below the Mg - sigma relation determined for elliptical galaxies and that the indices show a clear trend with morphological type. From the line-strength maps we calculate age, metallicity and abundance ratio maps and discuss the results from a one-SSP approach and from a two-SSP approach. Late-type galaxies are generally younger and more metal poor than ellipticals and have abundance ratios closer to solar values. We also explore a continuous star formation scenario, and try to recover the star formation history using the evolutionary models of Bruzual & Charlot (2003), assuming constant or exponentially declining star formation rate (SFR). We find a correlation between the e-folding time-scale tau of the starburst and the central velocity dispersion: more massive galaxies tend to have shorter tau, suggesting that the star formation happened long ago and has now basically ended, while for smaller objects with larger values of tau it is still active.

We present absorption line strength maps of a sample of 24 representative early-type spiral galaxies, mostly of type Sa, obtained as part of the SAURON survey of nearby galaxies using our custom-built integral-field spectrograph. Using high-quality spectra, spatially binned to a constant signal-to-noise, we measure several key age, metallicity and abundance ratio sensitive indices from the Lick/IDS system over a contiguous two-dimensional field including bulge and inner disc. We present maps of H beta, Fe 5015 and Mg b, for each galaxy The absorption line maps show that many galaxies contain some younger populations (<= 1 Gyr), distributed in small or large inner discs, or in circumnuclear star forming rings. In many cases these young stars are formed in circumnuclear mini-starbursts, which are dominating the light in the centres of some of the early-type spirals. These mini-starburst cause a considerable scatter in index-index diagrams such as Mg b- H beta and Mg b -Fe 5015, more than is measured for early-type galaxies. We find that the central regions of Sa galaxies display a wide range in ages, even within the galaxies. 50% of the sample show velocity dispersion drops in their centres. All of the galaxies of our sample lie on or below the Mg b- $\sigma$ relation for elliptical galaxies in the Coma cluster, and above the H beta absorption line - $\sigma$ relation for elliptical galaxies. If those relations are considered to be relations for the oldest local galaxies we see that our sample of spirals has a considerable scatter in age, with the largest scatter at the lowest $\sigma$. This is in disagreement with highly inclined samples, in which generally only old stellar populations are found in the central regions. All this can be understood if ... (see paper for rest of abstract)

We analyse the orbital distribution of E and S0 galaxies using SAURON integral-field stellar kinematics. We construct the anisotropy (V/sigma,epsilon) diagram, for the 48 E/S0 galaxies from the SAURON survey. For a subsample of 24 galaxies consistent with axisymmetry, we use three-integral axisymmetric Schwarzschild dynamical models to recover the detailed orbital distribution and we find good agreement with the anisotropy derived from the (V/sigma,epsilon) diagram. In the companion Paper IX we show that the early-type galaxies can be subdivided into two classes of systems with or without a significant amount of specific stellar angular momentum. Here we show that the two classes have different distributions on the (V/sigma,epsilon) diagram. The slow rotators are more common among the most massive systems and are generally classified as E from photometry alone. Those in our sample tend to be fairly round (epsilon<0.3), but can have significant kinematical misalignments, indicating that as a class they are moderately triaxial, and span a range of anisotropies (delta<0.3). The fast rotators are generally fainter and are classified either E or S0. They can appear quite flattened (epsilon<0.7), do not show significant kinematical misalignments (unless barred or interacting), indicating they are nearly axisymmetric, and span an even larger range of anisotropies (delta<0.5). These results are confirmed when we extend our analysis to 18 additional E/S0 galaxies observed with SAURON. All these results support the idea that fast rotators are nearly oblate and contain disk-like components. The role of gas must have been important for their formation. The slow rotators are weakly triaxial. Current collisionless merger models seem unable to explain their detailed observed properties. (Abridged)

We discuss SAURON absorption line strength maps of a sample of 24 early-type spirals, mostly Sa. From the Lick indices H$\beta$, Mg b and Fe 5015 we derive SSP-ages and metallicities. By comparing the scaling relations of Mg b and H beta and central velocity dispersion with the same relation for the edge-on sample of Falcon-Barroso et al. (2002) we derive a picture in which the central regions of Sa galaxies contain at least 2 components: one (or more) thin, disc-like component, often containing recent star formation, and another, elliptical-like component, consisting of old stars and rotating more slowly, dominating the light above the plane. If one defines a bulge to be the component responsible for the light in excess of the outer exponential disc, then many Sa-bulges are dominated by a thin, disc-like component containing recent star formation.

By detecting ionised-gas emission in 75% of the cases, the SAURON integral-field spectroscopic survey has further demonstrated that early-type galaxies often display nebular emission. Furthermore, the SAURON data have shown that such emission comes with an intriguing variety of morphologies, kinematic behaviours and line ratios. Perhaps most puzzling was the finding that round and slowly rotating objects generally display uncorrelated stellar and gaseous angular momenta, consistent with an external origin for the gas, whereas flatter and fast rotating galaxies host preferentially co-rotating gas and stars, suggesting internal production of gas. Alternatively, a bias against the internal production of ionised gas and against the acquisition of retrograde material may be present in these two kinds of objects, respectively. In light of the different content of hot gas in these systems, with slowly rotating objects being the only systems capable of hosting massive X-ray halos, we suggest that a varying importance of evaporation of warm gas in the hot interstellar medium can contribute to explain the difference in the relative behaviour of gas and stars in these two kinds of objects. Namely, whereas in X-ray bright and slowly rotating galaxies stellar-loss material would quickly evaporate in the hot medium, in X-ray faint and fast rotating objects such material would be allowed to lose angular momentum and settle in a disk, which could also obstruct the subsequent acquisition of retrograde gas. Evidence for a connection between warm and hot gas phases, presumably driven by heat conduction, is presented for four slowly rotating galaxies with CHANDRA observations.

We give two examples of spiral galaxies that show non-circular gas motions in the inner kiloparsecs, from SAURON integral field spectroscopy. We use harmonic decomposition of the velocity field of the ionized gas to study the underlying mass distribution, employing linear theory. The higher order harmonic terms and the main kinematic features of the observed data are consistent with an analytically constructed simple bar model. We also present maps of a number of strong absorption lines in M 100, derive simple stellar populations and correlate them with features in the gas kinematics.

We present high spatial resolution integral-field spectroscopy of 28 elliptical (E) and lenticular (S0) galaxies from the SAURON representative survey obtained with the OASIS spectrograph during its operation at the CFHT. These seeing-limited observations explore the central 8"x10" (typically one kiloparsec diameter) regions of these galaxies using a spatial sampling four times higher than SAURON (0.27" vs. 0.94" spatial elements), resulting in almost a factor of two improvement in the median PSF. These data allow accurate study of the central regions to complement the large-scale view provided by SAURON. Here we present the stellar and gas kinematics, stellar absorption-line strengths and nebular emission-line strengths for this sample. We also characterise the stellar velocity maps using the "kinemetry" technique, and derive maps of the luminosity-weighted stellar age, metallicity and abundance ratio via stellar population models. We give a brief review of the structures found in our maps, linking also to larger-scale structures measured with SAURON. We present two previously unreported kinematically-decoupled components (KDCs) in the centres of NGC3032 and NGC4382. We compare the intrinsic size and luminosity-weighted stellar age of all the visible KDCs in the full SAURON sample, and find two types of components: kiloparsec-scale KDCs, which are older than 8 Gyr, and are found in galaxies with little net rotation; and compact KDCs, which have intrinsic diameters of less than a few hundred parsec, show a range of stellar ages from 0.5 - 15 Gyr (with 5/6 younger than 5 Gyr), are found exclusively in fast-rotating galaxies, and are close to counter-rotating around the same axis as their host. (abridged)

(Abridged) We present deep WSRT observations of neutral hydrogen in 12 nearby elliptical and lenticular galaxies. The selected objects come from a representative sample of nearby galaxies earlier studied at optical wavelengths with the integral-field spectrograph Sauron. We detect HI - both in regular discs as well as in clouds and tails off-set from the host galaxy - in 70% of the galaxies. This detection rate is much higher than in previous, shallower single-dish surveys, and is similar to that for the ionised gas. The results suggest that at faint detection levels the presence of HI is a relatively common characteristic of field early-type galaxies. The observed total HI masses range between a few times 10^6 Msun to just over 10^9 Msun. The presence of regular disc-like structures is a situation as common as HI in offset clouds and tails around early-type galaxies. All galaxies where HI is detected also contain ionised gas, whereas no HI is found around galaxies without ionised gas. Galaxies with regular HI discs tend to have strong emission from ionised gas. In these cases, the similar kinematics of the neutral hydrogen and ionised gas suggest that they form one structure. We do not find a clear trend between the presence of HI and the global age of the stellar population or the global dynamical characteristics of the galaxies. Our observations support the idea that gas accretion is common and does not happen exclusively in peculiar early-type galaxies. For the majority of the cases the gas is likely acquired through merging, but the lack of correlation with the stellar population age suggests that smooth, cold accretion could be an alternative scenario, at least in some galaxies. In either cases, the data suggest that early-type galaxies continue to build their mass up to the present.

We combine Hubble Space Telescope spectroscopy and ground-based integral-field data from the SAURON and OASIS instruments to study the central black hole in the nearby elliptical galaxy NGC 3379. From these data, we obtain kinematics of both the stars and the nuclear gaseous component. Axisymmetric three-integral models of the stellar kinematics find a black hole of mass 1.4 (+2.6 / -1.0) x 10^8 M_sun (3 sigma errors). These models also probe the velocity distribution in the immediate vicinity of the black hole and reveal a nearly isotropic velocity distribution throughout the galaxy and down to the black hole sphere of influence R_BH. The morphology of the nuclear gas disc suggests that it is not in the equatorial plane; however the core of NGC 3379 is nearly spherical. Inclined thin-disc models of the gas find a nominal black hole of mass 2.0 (+/- 0.1) x 10^8 M_sun (3 sigma errors), but the model is a poor fit to the kinematics. The data are better fit by introducing a twist in the gas kinematics (with the black hole mass assumed to be 2.0 x 10^8 M_sun), although the constraints on the nature and shape of this perturbation are insufficient for more detailed modelling. Given the apparent regularity of the gas disc's appearance, the presence of such strong non-circular motion indicates that caution must be used when measuring black hole masses with gas dynamical methods alone.

(Abridged) We present observations of the stellar and gas kinematics for a representative sample of 24 Sa galaxies obtained with our custom-built integral-field spectrograph SAURON operating on the William Herschel Telescope. Our maps typically cover the bulge dominated region. We find a significant fraction of kinematically decoupled components (12/24), many of them displaying central velocity dispersion minima. They are mostly aligned and co-rotating with the main body of the galaxies, and are usually associated with dust discs and rings detected in unsharp-masked images. Almost all the galaxies in the sample (22/24) contain significant amounts of ionised gas which, in general, is accompanied by the presence of dust. The kinematics of the ionised gas is consistent with circular rotation in a disc co-rotating with respect to the stars. The distribution of mean misalignments between the stellar and gaseous angular momenta in the sample suggest that the gas has an internal origin. The [OIII]/Hbeta ratio is usually very low, indicative of current star formation, and shows various morphologies (ring-like structures, alignments with dust lanes or amorphous shapes). The star formation rates in the sample are comparable with that of normal disc galaxies. Low gas velocity dispersion values appear to be linked to regions of intense star formation activity. We interpret this result as stars being formed from dynamically cold gas in those regions. In the case of NGC5953, the data suggest that we are witnessing the formation of a kinematically decoupled component from cold gas being acquired during the ongoing interaction with NGC5954.

We present results from our ongoing effort to understand the nature and evolution of nearby galaxies using the SAURON integral-field spectrograph. In this proceeding we focus on the study of the particular case formed by the interacting galaxies NGC5953 and NGC5954. We present stellar and gas kinematics of the central regions of NGC5953. We use a simple procedure to determine the age of the stellar populations in the central regions and argue that we may be witnessing the formation of a kinematically decoupled component from cold gas being acquired during the ongoing interaction with NGC5954.

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.

We present absorption line strength maps of 48 representative elliptical and lenticular galaxies obtained as part of the SAURON survey of nearby galaxies. Using high-quality spectra, spatially binned to a constant signal-to-noise, we measure four key age, metallicity and abundance ratio sensitive indices from the Lick/IDS system over a two-dimensional field extending up to approximately one effective radius. We modify the classical Fe5270 index to define a new index, Fe5270S, which maximizes the usable spatial coverage of SAURON. Maps of Hbeta, Fe5015, Mgb and Fe5270S are presented for each galaxy. We use the maps to compute average line strengths integrated over circular apertures of one-eighth effective radius, and compare the resulting relations of index versus velocity dispersion with previous long-slit work. The metal line strength maps show generally negative gradients with increasing radius roughly consistent with the morphology of the light profiles. Remarkable deviations from this general trend exist, particularly the Mgb isoindex contours appear to be flatter than the isophotes of the surface brightness for about 40% of our galaxies without significant dust features. Generally these galaxies exhibit significant rotation. We infer from this that the fast-rotating component features a higher metallicity and/or an increased Mg/Fe ratio as compared to the galaxy as a whole. The Hbeta maps are typically flat or show a mild positive outwards radial gradient, while a few galaxies show strong central peaks and/or elevated overall Hbeta-strength likely connected to recent star-formation activity. For the most prominent post-starburst galaxies even the metal line strength maps show a reversed gradient. (abridged)

We present selected results from integral-field spectroscopy of 48 early-type galaxies observed as part of the SAURON survey. Maps of the Hbeta, Fe5015, Mgb and Fe5270 indices in the Lick/IDS system were derived for each of the survey galaxies. The metal line strength maps show generally negative gradients with increasing radius roughly consistent with the morphology of the light profiles. Remarkable deviations from this general trend exist, particularly the Mgb isoindex contours appear to be flatter than the isophotes of the surface brightness for about 40% of our galaxies without significant dust features. Generally these galaxies exhibit significant rotation. We infer from this that the fast-rotating component features a higher metallicity and/or an increased Mg/Fe ratio as compared to the galaxy as a whole. We also use the line strengths maps to compute average values integrated over circular apertures of one effective radius, and derive luminosity weighted ages and metallicities. The lenticular galaxies show a wide range in age and metallicity estimates, while elliptical galaxies tend to occupy regions of older stellar populations.

We present high spatial resolution near-infrared spectra and images of the nucleus of Centaurus A (NGC 5128) obtained with NAOS-CONICA at the VLT. The adaptive optics corrected data have a spatial resolution of 0.06" (FWHM) in K- and 0.11" in H-band, four times higher than previous studies. The observed gas motions suggest a kinematically hot disk which is orbiting a central object and is oriented nearly perpendicular to the nuclear jet. We model the central rotation and velocity dispersion curves of the [FeII] gas orbiting in the combined potential of the stellar mass and the (dominant) black hole. Our physically most plausible model, a dynamically hot and geometrically thin gas disk, yields a black hole mass of M_bh = (6.1 +0.6/-0.8) 10^7 M_sun. As the physical state of the gas is not well understood, we also consider two limiting cases: first a cold disk model, which completely neglects the velocity dispersion; it yields an M_bh estimate that is almost two times lower. The other extreme case is to model a spherical gas distribution in hydrostatic equilibrium through Jeans equation. Compared to the hot disk model the best-fit black hole mass increases by a factor of 1.5. This wide mass range spanned by the limiting cases shows how important the gas physics is even for high resolution data. Our overall best-fitting black hole mass is a factor of 2-4 lower than previous measurements. With our revised M_bh estimate, Cen A's offset from the M_bh-sigma relation is significantly reduced; it falls above this relation by a factor of ~2, which is close to the intrinsic scatter of this relation. (Abridged)

We present the stellar and gas kinematics of a sample of 18 nearby late-type spiral galaxies (Hubble types ranging from Sb to Sd), observed with the integral-field spectrograph SAURON at the 4.2-m William Herschel Telescope. SAURON covers the spectral range 4800-5380 A, allowing us to measure the Hbeta, Fe, Mgb absorption features and the emission in the Hbeta line and the [OIII], and [NI] doublets over a 33x41 arcsec field of view. The maps cover the nuclear region of these late-type galaxies and in all cases include the entire bulge. In many cases the stellar kinematics suggests the presence of a cold inner region, as visible from a central drop in the stellar velocity dispersion. The ionised gas is almost ubiquitous and behaves in a complicated fashion: the gas velocity fields often display more features than the stellar ones, including wiggles in the zero-velocity lines, irregular distributions, ring-like structures. The line ratio [OIII]/Hbeta often takes on low values over most of the field, probably indicating a wide-spread star formation.

We present a generalisation of surface photometry to the higher-order moments of the line-of-sight velocity distribution of galaxies observed with integral-field spectrographs. The generalisation follows the approach of surface photometry by determining the best fitting ellipses along which the profiles of the moments can be extracted and analysed by means of harmonic expansion. The assumption for the odd moments (e.g. mean velocity) is that the profile along an ellipse satisfies a simple cosine law. The assumption for the even moments (e.g velocity dispersion) is that the profile is constant, as it is used in surface photometry. We find that velocity profiles extracted along ellipses of early-type galaxies are well represented by the simple cosine law (with 2% accuracy), while possible deviations are carried in the fifth harmonic term which is sensitive to the existence of multiple kinematic components, and has some analogy to the shape parameter of photometry. We compare the properties of the kinematic and photometric ellipses and find that they are often very similar. Finally, we offer a characterisation of the main velocity structures based only on the kinemetric parameters which can be used to quantify the features in velocity maps (abridged).

We have developed an orbit-based method for constructing triaxial models of elliptical galaxies, which fit their observed surface brightness and kinematics. We have tested the method against analytical models with general distribution functions. Here, we present models that fit integral-field SAURON observations of NGC3379 and NGC4365.

We investigate with unprecedented accuracy the correlations between the dynamical mass-to-light ratio M/L and other global observables of E and S0 galaxies. We construct two-integral Jeans and three-integral Schwarzschild dynamical models for a sample of 25 E/S0 galaxies with SAURON integral-field stellar kinematics. We find a tight correlation of the form (M/L)=(3.80+/-0.14)*(sigma/200 km/s)^(0.84+/-0.07) between the dynamical M/L (in the I-band) and the luminosity-weighted second moment (sigma) of the line-of-sight velocity-distribution within Re. The observed rms scatter in M/L for our sample is 18%, while the inferred intrinsic scatter is ~13%. The (M/L)-sigma relation can be included in the remarkable series of tight correlations between sigma and other galaxy global observables. The comparison of the observed correlations with the predictions of the Fundamental Plane (FP), and with simple virial estimates, shows that the `tilt' of the FP of early-type galaxies, is due to a real M/L variation, while structural and orbital non-homology have a negligible effect. The virial mass is a reliable estimator of the mass in the central regions of galaxies. The best-fitting virial relation has the form (M/L)_vir=(5.0+/-0.1)*Re*sigma^2/(L*G). The comparison of the dynamical M/L with the (M/L)_pop inferred from the analysis of the stellar population, indicates a median dark matter fraction in early-type galaxies of ~30% of the total mass inside one Re. (Abridged)

We present results from our ongoing effort to understand the morphological and kinematical properties of early-type galaxies using the integral-field spectrograph SAURON. We discuss the relation between the stellar and gas morphology and kinematics in our sub-sample of 24 representative Sa spiral bulges. We focus on the frequency of kinematically decoupled components and on the presence of star formation in circumnuclear rings.

We investigate the well-known correlations between the dynamical mass-to-light ratio M/L and other global observables of elliptical (E) and lenticular (S0) galaxies. We construct two-integral Jeans and three-integral Schwarzschild dynamical models for a sample of 25 E/S0 galaxies with SAURON integral-field stellar kinematics to about one effective (half-light) radius R_e. The comparison of the dynamical M/L with the (M/L)_pop inferred from the analysis of the stellar population, indicates that dark matter in early-type galaxies contributes ~30% of the total mass inside one R_e, in agreement with previous studies, with significant variations from galaxy to galaxy. Our results suggest a variation in M/L at constant (M/L)_pop, which seems to be linked to the galaxy dynamics. We speculate that fast rotating galaxies have lower dark matter fractions than the slow rotating and generally more massive ones.

We present the emission-line fluxes and kinematics of 48 representative elliptical and lenticular galaxies obtained with our custom-built integral-field spectrograph SAURON. Hb, [OIII], and [NI] emission lines were measured using a new procedure that simultaneously fits both the stellar spectrum and the emission lines. Using this technique we can detect emission lines down to an equivalent width of 0.1A set by the current limitations in describing galaxy spectra with synthetic and real stellar templates, rather than by the quality of our spectra. Emission is detected in 75% of our sample galaxies, and comes in a variety of resolved spatial distributions and kinematic behaviours. The ionised-gas kinematics is rarely consistent with simple coplanar circular motions. However, the gas almost never displays completely irregular kinematics, generally showing coherent motions with smooth variations in angular momentum. In the majority of the cases the gas kinematics is decoupled from the stellar kinematics, and in half of the objects this decoupling implies a recent acquisition of gaseous material. Over the entire sample however, the distribution of the mean misalignment values between stellar and gaseous angular momenta is inconsistent with a purely external origin. Consistent with previous studies, the presence of dust features is always accompanied by gas emission while the converse is not always true. A considerable range of values for the [OIII]/Hb ratio is found both across the sample and within single galaxies. Despite the limitations of this ratio as an emission-line diagnostic, this finding suggests either that a variety of mechanisms is responsible for the gas excitation in E and S0 galaxies or that the metallicity of the interstellar material is quite heterogeneous.

In the past decade, most effort in the study of supermassive black holes (BHs) has been devoted to measuring their masses. This lead to the finding of the tight M_BH-sigma relation, which indicates the existence of strong links between the formation of the BH and of their host spheroids. Many scenarios have been proposed to explain this relation, and all agree on the key role of BHs' growth and feedback in shaping their host galaxies. However the currently available observational constraints, essentially BH masses and galaxy photometry, are not sufficient to conclusively select among the alternatives. A crucial piece of information of the black hole formation is recorded in the orbital distribution of the stars, which can only be extracted from high-resolution integral-field (IF) stellar kinematics. The introduction of IF spectrographs with adaptive optics on large telescopes opens a new era in the study of BHs by finally allowing this key element to be uncovered. This information will be complementary to what will be provided by the LISA gravitational wave satellite, which can directly detect coalescing BHs. Here an example is presented for the recovery of the orbital distribution in the center of the giant elliptical galaxy M87, which has a well resolved BH sphere of influence, using SAURON IF kinematics.

Morphological disturbances and gas kinematics of the SB0 galaxy NGC 7679=Arp 216 are investigated to get clues to the history of this highly composite object, where AGN and starburst signatures dominate each other in the X-ray and optical/IR regime, respectively. Perturbations of the ionized gas velocity field appear quite mild within 15'' (~5 kpc) from the center, so as it can be straightforwardly modeled as a circularly rotating disk. On the contrary, outside that radius, significant disturbances show up. In particular, the eastern distorted arm as well as the huge neutral hydrogen bridge connecting NGC 7679 with the nearby Seyfert spiral NGC 7682 unambiguously represent the vestige of a close encounter of the two objects dating back ~500 Myr ago. The relationship of such past event with the much more recent, centrally located starburst (not older than 20 Myr) cannot be easily established. Altogether, the classification of NGC 7679, turns out to be less extreme than that proposed in the past, being simply a (disturbed) galaxy where starburst and AGN activity cohexist with a starburst dominating the bolometric luminosity.

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.

We use integral-field observation of the stellar kinematics obtained with SAURON in combination with Schwarzschild dynamical models to revisit our understanding of the classic (V/sigma)-epsilon anisotropy diagram of early-type galaxies.

We summarise the results and achievements of integral-field spectroscopy of early-type galaxies, observed as part of a survey using both the SAURON and OASIS spectrographs. From the perspective of integral-field spectroscopy, these otherwise smooth and featureless objects show a wealth of structure, both in their stellar kinematics and populations. We focus on the stellar content, and examine properties on both kiloparsec scales with SAURON, and scales of 100's of parsecs with OASIS. These complementary studies reveal two types of kinematically distinct components (KDCs), differing primarily in their intrinsic sizes. In previous studies, KDCs and their host galaxies have generally been found to be unremarkable in other aspects. We show that large KDCs, typical of the well-studied cases, indeed show little or no age differences with their host galaxy. The KDCs detected with the higher spatial-resolution of OASIS are intrinsically smaller and include, in contrast, a significant fraction of young stars. We speculate on the relationship between KDCs and their host galaxies, and the implications for young populations in early-type galaxies.

We present results of our ongoing study of the morphology and kinematics of the ionised gas in 48 representative nearby elliptical and lenticular galaxies using the SAURON integral-field spectrograph on the 4.2m William Herschel Telescope. Making use of a recently developed technique, emission is detected in 75% of the galaxies. The ionised-gas distributions display varied morphologies, ranging from regular gas disks to filamentary structures. Additionally, the emission-line kinematic maps show, in general, regular motions with smooth variations in kinematic position angle. In most of the galaxies, the ionised-gas kinematics is decoupled from the stellar counterpart, but only some of them present signatures of recent accretion of gaseous material. The presence of dust is very common in our sample and is usually accompanied by gas emission. Our analysis of the [OIII]/Hbeta emission-line ratios, both across the whole sample as well as within the individual galaxies, suggests that there is no unique mechanism triggering the ionisation of the gas.

We analyse SAURON kinematic maps of the inner kpc of the early-type (Sa) barred spiral galaxy NGC 5448. 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 displays 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, which is derived using linear theory. 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.

We study the large-scale stellar and gaseous kinematics of the E4 galaxy NGC2974, based on panoramic integral-field data obtained with SAURON. We quantify the velocity fields with Fourier methods (kinemetry), and show that the large-scale kinematics is largely consistent with axisymmetry. We construct general axisymmetric dynamical models for the stellar motions using Schwarzschild's orbit-superposition method, and compare the inferred inclination and mass-to-light ratio with the values obtained by modelling the gas kinematics. Both approaches give consistent results. However we find that the stellar models provide fairly weak constraints on the inclination. The intrinsic orbital distribution of NGC2974, which we infer from our model, is characterised by a large-scale stellar component of high angular momentum. We create semi-analytic test models, resembling NGC2974, to study the ability of Schwarzschild's modelling technique to recover the given input parameters (mass-to-light ratio and inclination) and the distribution function. We also test the influence of a limited spatial coverage on the recovery of the distribution function (i.e. the orbital structure). We find that the models can accurately recover the input mass-to-light ratio, but we confirm that even with perfect input kinematics the inclination is only marginally constrained. This suggests a possible degeneracy in the determination of the inclination, but further investigations are needed to clarify this issue. For a given potential, we find that the analytic distribution function of our test model is well recovered by the three-integral model within the spatial region constrained by integral-field kinematics.

NGC5018, one of the weakest UV emitters among giant ellipticals (gE) observed with IUE, appears to consist of an optical stellar population very similar to that of the compact, dwarf elliptical M32, which is several magnitudes fainter in luminosity than NGC5018 and whose stellar population is know to be ~3 Gyr old. Here we show that the mid-UV spectra of these two galaxies are also very similar down to an angular scale hundreds times smaller than the IUE large aperture (as probed by HST/FOS UV spectra obtained through 0.86 arcsec apertures). This implies a reasonably close match of the populations dominating their mid-UV light (namely, their main-sequence turnoff stars). These data indicate that NGC5018 has, in its inner regions, a rather uniform dominance of a ~3 Gyr-old stellar population, probably a bit different in metallicity from M32. Combined with the various structures that indicate that NGC5018 is the result of a recent major merger, it appears that almost all of stars we see in its center regions were formed about 3 Gyr ago, in that merger event. NGC5018 is likely the older brother of NGC7252, the canonical gE-in-formation merger. As such, NGC5018 is perhaps the best galaxy which can tell us how a merger works, after the fireworks subside, to form a gE galaxy today. For this reason alone, the stellar populations in NGC5018 at all radii are worth studying in detail.

We present the stellar kinematics of 48 representative elliptical and lenticular galaxies obtained with our custom-built integral-field spectrograph SAURON operating on the William Herschel Telescope. The data were homogeneously processed through a dedicated reduction and analysis pipeline. All resulting SAURON datacubes were spatially binned to a constant minimum signal-to-noise. We have measured the stellar kinematics with an optimized (penalized pixel-fitting) routine which fits the spectra in pixel space, via the use of optimal templates, and prevents the presence of emission lines to affect the measurements. We have thus generated maps of the mean stellar velocity, the velocity dispersion, and the Gauss-Hermite moments h3 and h4 of the line-of-sight velocity distributions. The maps extend to approximately one effective radius. Many objects display kinematic twists, kinematically decoupled components, central stellar disks, and other peculiarities, the nature of which will be discussed in future papers of this series.

We present SAURON integral-field observations of the S0 galaxy NGC7332. Existing broad-band ground-based and HST photometry reveals a double disk structure and a boxy bulge interpreted as a bar viewed close to edge-on. The SAURON two-dimensional stellar kinematic maps confirm the existence of the bar and inner disk but also uncover the presence of a cold counter-rotating stellar component within the central 250 pc. The Hbeta and [OIII] emission line maps show that the ionised gas has a complex morphology and kinematics, including both a component counter-rotating with respect to the stars and a fainter co-rotating one. Analysis of the absorption line-strength maps show that NGC7332 is young everywhere. The presence of a large-scale bar can explain most of those properties, but the fact that we see a significant amount of unsettled gas, together with a few peculiar features in the maps, suggest that NGC7332 is still evolving. Interactions as well as bar-driven processes must thus have played an important role in the formation and evolution of NGC7332, and presumably of S0 galaxies in general.

We investigate the accuracy of the parametric recovery of the line-of-sight velocity distribution (LOSVD) of the stars in a galaxy, while working in pixel space. Problems appear when the data have a low signal-to-noise ratio, or the observed LOSVD is not well sampled by the data. We propose a simple solution based on maximum penalized likelihood and we apply it to the common situation in which the LOSVD is described by a Gauss-Hermite series. We compare different techniques by extracting the stellar kinematics from observations of the barred lenticular galaxy NGC 3384 obtained with the SAURON integral-field spectrograph.

We present a summary of high-spatial resolution follow-up observations of the elliptical (E) and lenticular (S0) galaxies in the SAURON survey using the OASIS integral field spectrograph. The OASIS observations explore the central 8x10" regions of these galaxies using a spatial sampling four times higher than SAURON, often revealing previously undiscovered features. Around 75% (31/48) of the SAURON E/S0s with central velocity dispersion >= 120 km/s were observed with OASIS, covering well the original SAURON representative sample. We present here an overview of this follow-up survey, and some preliminary results on individual objects, including a previously unreported counter-rotating core in NGC 4382; the decoupled stellar and gas velocity fields of NGC 2768; and the strong age gradient towards the centre of NGC 3489.

Kinematic and population studies show that bulges are generally rotationally flattened systems similar to low-luminosity ellipticals. However, observations with state-of-the-art integral field spectrographs, such as SAURON, indicate that the situation is much more complex, and allow us to investigate phenomena such as triaxiality, kinematic decoupling and population substructure, and to study their connection to current formation and evolution scenarios for bulges of early-type galaxies. We present the examples of two S0 bulges from galaxies in our sample of nearby galaxies: one that shows all the properties expected from classical bulges (NGC5866), and another case that presents kinematic features appropriate for barred disk galaxies (NGC7332).

We present dynamical modeling of nearby E4 galaxy NGC 2974. The galaxy was observed with the integral-field spectrograph SAURON. An axisymmetric stellar dynamical model is constructed using a three-integral orbit superposition Schwarschild modeling technique. We explore the dynamical structure and investigate the distribution and dynamics of the gas.

We have developed a method to construct realistic triaxial dynamical models for elliptical galaxies, allowing us to derive best-fitting parameters, such as the mass-to-light ratio and the black hole mass, and to study the orbital structure. We use triaxial theoretical Abel models to investigate the robustness of the method.

We describe our program for the dynamical modeling of early-type galaxies observed with the panoramic integral-field spectrograph SAURON. We are using Schwarzschild's numerical orbit superposition method to reproduce in detail all kinematical and photometric observables, and recover the intrinsic orbital structure of the galaxies. Since catastrophes are the most prominent features in the orbital observables, two-dimensional kinematical coverage is essential to constrain the dynamical models.

We present new techniques to perform adaptive spatial binning of Integral-Field Spectroscopic (IFS) data to reach a chosen constant signal-to-noise ratio per bin. These methods are required for the proper analysis of IFS observations, but can also be used for standard photometric imagery or any other two-dimensional data. Various schemes are tested and compared by binning and extracting the stellar kinematics of the Sa galaxy NGC2273 from spectra obtained with the panoramic IFS SAURON.

We have measured the two-dimensional kinematics and line-strength distributions of 72 representative nearby early-type galaxies, out to approximately one effective radius, with our panoramic integral-field spectrograph SAURON. The resulting maps reveal a rich variety in kinematical structures and linestrength distributions, indicating that early-type galaxies are more complex systems than often assumed. We are building detailed dynamical models for these galaxies, to derive their intrinsic shape and dynamical structure, and to determine the mass of the supermassive central black hole. Here we focus on two examples, the compact elliptical M32 and the E3 galaxy NGC4365. These objects represent two extreme cases: M32 has very regular kinematics which can be represented accurately by an axisymmetric model in which all stars rotate around the short axis, while NGC4365 is a triaxial galaxy with a prominent kinematically decoupled core, with an inner core that rotates about an axis that is nearly perpendicular to the rotation axis of the main body of the galaxy. Our dynamical models for these objects demonstrate that two-dimensional observations are essential for deriving the intrinsic orbital structure and dark matter content of galaxies.

The E3 giant elliptical galaxy IC1459 is the prototypical galaxy with a fast counterrotating stellar core. We obtained one HST/STIS long-slit spectrum along the major axis of this galaxy and CTIO spectra along five position angles. We present self-consistent three-integral axisymmetric models of the stellar kinematics, obtained with Schwarzschild's numerical orbit superposition method. We study the dynamics of the kinematically decoupled core (KDC) in IC1459 and we find it consists of stars that are well-separated from the rest of the galaxy in phase space. The stars in the KDC counterrotate in a disk on orbits that are close to circular. We estimate that the KDC mass is ~0.5% of the total galaxy mass or ~3*10^9 Msun. We estimate the central black hole mass M_BH of IC1459 independently from both its stellar and its gaseous kinematics. Some complications probably explain why we find rather discrepant BH masses with the different methods. The stellar kinematics suggest that M_BH = (2.6 +/- 1.1)*10^9 Msun (3 sigma error). The gas kinematics suggests that M_BH ~ 3.5*10^8 Msun if the gas is assumed to rotate at the circular velocity in a thin disk. If the observed velocity dispersion of the gas is assumed to be gravitational, then M_BH could be as high as ~1.0*10^9 Msun. These different estimates bracket the value M_BH = (1.1 +/- 0.3)*10^9 Msun predicted by the M_BH-sigma relation. It will be an important goal for future studies to assess the reliability of black hole mass determinations with either technique. This is essential if one wants to interpret the correlation between the BH mass and other global galaxy parameters (e.g. velocity dispersion) and in particular the scatter in these correlations (believed to be only ~0.3 dex). [Abridged]

We present results from a new and unique integral-field spectrograph, SAURON. It has a large field of view and high throughput and is primarily built for the study of stellar & gaseous kinematics and stellar populations in galaxies. Its aim is to carry out a systematic survey of the velocity fields, velocity dispersions, and line-strength distributions of nearby ellipticals, lenticular galaxies and spiral bulges. Its wide field is especially useful for the study of complicated velocity structures. Together with other spectroscopic data, images, and dynamical modelling, SAURON will help to constrain the intrinsic shapes, mass-to-light ratios, and stellar populations of early-type galaxies and spiral bulges.

We present dynamical models of the nearby compact elliptical galaxy M32, using high quality kinematical measurements, obtained with the integral-field spectrograph SAURON mounted on the William Herschel Telescope on La Palma. We also include STIS data obtained by Joseph et al. We find a best-fit black hole mass of (2.5 +/- 0.5) million solar masses and a stellar I-band mass-to-light ratio of (1.85 +/- 0.15) in solar units. For the first time, we are also able to constrain the inclination along which M32 is observed to (70 +/- 5) degrees. Combined with an averaged observed flattening of 0.73, this corresponds to an intrinsic flattening of approximately 0.68 +/- 0.03. These tight constraints are mainly caused by the use of integral-field data. We show this quantitatively by comparing with models that are constrained by multiple slits only. We show the phase-space distribution and intrinsic velocity structure of the best-fit model and investigate the effect of regularisation on the orbit distribution.

We present the SAURON project, which is aimed at studying the morphology, two-dimensional kinematics and stellar populations of a representative sample of elliptical galaxies and spiral bulges. SAURON, a dedicated integral-field spectrograph that is optimized for wide-field observations and has high throughput, was built in Lyon and is now operated at the WHT 4.2m telescope. At present, we have observed approximately two thirds of the seventy-two sample galaxies with SAURON. A comparison with published long-slit measurements demonstrates that the SAURON-data is of equal or better quality, and provides full two-dimensional coverage. The velocity and velocity dispersion fields exhibit a large variety of morphologies: from simple rotating systems to cylindrical, disky and triaxial velocity fields, bars and decoupled cores. Most of these kinematical signatures do not have counterparts in the light distribution. While some galaxies are consistent with axisymmetry, most are more complex systems than assumed previously. This suggests that the kinematical properties of nearby E/S0 galaxies do not agree with the often assumed simplistic two-family model, in which the giant non-rotating triaxial ellipticals are opposed to the fast-rotating axisymmetric faint ellipticals and S0s.

The panoramic integral-field spectrograph SAURON is currently being used to map the stellar kinematics, gaseous kinematics, and stellar populations of a large number of early-type galaxies and bulges. Here, we describe SAURON observations of cold stellar disks embedded in spheroids (NGC3384, NGC4459, NGC4526), we illustrate the kinematics and ionization state of large-scale gaseous disks (NGC4278, NGC7742), and we show preliminary comparisons of SAURON data with barred galaxy N-body simulations (NGC3623).

We use long-slit spectroscopy along the major axis of a sample of 23 nearby disk galaxies to study the kinematic properties of the ionized-gas component in their inner regions. For each galaxy, we derive the position-velocity diagram of the ionized gas from its emission lines. We discuss the variety of shapes observed in such position-velocity diagrams by comparing the gas velocity gradient, velocity dispersion and integrated flux measured in the inner (r +/-1'') and outer regions (r +/-4''). This kind of analysis allows the identification of galaxies which are good candidates to host a circumnuclear Keplerian gaseous disk rotating around a central mass concentration, and to follow up with Hubble Space Telescope observations.

The SAURON project will deliver two-dimensional spectroscopic data of a sample of nearby early-type galaxies with unprecedented quality. In this paper, we focus on the mapping of their stellar populations using the SAURON data, and present some preliminary results on a few prototypical cases.

We present new techniques to perform adaptive spatial binning of two-dimensional (2D) data to reach a chosen constant signal-to-noise ratio per bin. These methods are required particularly for the proper analysis of Integral Field Spectrograph (IFS) observations, but can also be used for standard photometric imagery. Various schemes are tested and compared using data obtained with the panoramic IFS SAURON.

We describe a simple, efficient, robust and fully automatic algorithm for the determination of a Multi-Gaussian Expansion (MGE) fit to galaxy images, to be used as a parametrization for the galaxy stellar surface brightness. In most cases the least-squares solution found by this method essentially corresponds to the minimax, constant relative error, MGE approximation of the galaxy surface brightness, with the chosen number of Gaussians. The algorithm is well suited to be used with multiple resolution images (e.g., Hubble Space Telescope [HST] and ground-based). It works orders of magnitude faster and is more accurate than currently available methods. An alternative, more computing intensive, fully linear algorithm, that is guaranteed to converge to the smallest Chi^2 solution, is also discussed. Examples of MGE fits are presented for objects with HST or ground-based photometry, including galaxies with significant isophote twist.

A summary of the SAURON project and its current status is presented. SAURON is a panoramic integral-field spectrograph designed to study the stellar kinematics, gaseous kinematics, and stellar populations of spheroids. Here, the sample of galaxies and its properties are described. The instrument is detailed and its capabilities illustrated through observational examples. These includes results on the structure of central stellar disks, the kinematics and ionization state of gaseous disks, and the stellar populations of galaxies with decoupled cores.

The nucleus of the bulge-dominated, multiply-barred S0/a galaxy NGC 2681 is studied in detail, using high resolution Hubble Space Telescope FOC and NICMOS imaging and FOS spectroscopy. The ionised gas central velocity dispersion is found to increase by a factor ~2 when narrowing the aperture from R~1.5" (ground) to R~0.1" (FOS). Dynamical modeling of these velocity dispersions suggests that NGC 2681 does host a supermassive black hole (BH) for which one can estimate a firm mass upper limit M_BH < 6*10^7 Solar Masses. This upper limit is consistent with the relation between the central BH mass and velocity dispersion M_BH - sigma known for other galaxies. The emission line ratios place the nucleus of NGC 2681 among LINERs. It is likely that the emission line region comes from a rather mild, but steady, feeding of gas to the central BH in this galaxy. The inner stellar population lacks any measurable color gradient (to a radius of 0.6 kpc) from the infrared to the ultraviolet, consistently with FOC, FOS and IUE data, all indicating that this system underwent a starburst ~1 Gyr ago that encompassed its whole interior, down to its very center. The most likely source of such a widely-distributed starburst is the dumping of tidally-extruded gas from a galaxy neighbor. If so, then NGC 2681 can be considered as the older brother of M82, seen face-on as opposed to the edge-on view we have for M82.

Newly obtained HST/WFPC2 images of the disturbed elliptical galaxy NGC 5018 show that the average amount of internal reddening due to the its complex ``dust web'' is as low as E(B-V)~0.02 within the IUE aperture, thus implying that its observed and intrinsic energy distributions do not differ significantly down to UV wavelengths. This, in turn, is quite relevant to the current debate on the age of its dominant stellar population.

We report the case of the early-type disk galaxy NGC 4672 as a new example of a galaxy characterized by the orthogonal geometrical decoupling between bulge and disk. The morphological features of this galaxy are discussed as well as the velocity curves and velocity dispersion profiles of stars and ionized gas along both its major and minor axis. We conclude that NGC 4672 has structural (i.e., a bulge elongated perpendicularly to the disk) and kinematical (i.e., a stellar core rotating perpendicularly to the disk) properties similar to those of the Sa NGC 4698. The presence of the isolated core suggests that the disk component is the end result of the acquisition of external material in polar orbits around a pre-existing oblate spheroid as in the case of the ring component of AM 2020-504, the prototype of polar ring ellipticals.

The first UV F300W and F218W WFPC2 observations of the rich galaxy cluster CL0939+4713 at z=0.41 are presented and discussed. UV/optical two-color and c-m diagrams of the sources detected in the F300W waveband are constructed. Thanks to pre-existing HST optical images of the same field a morphological classification for the majority of these objects is also provided. Moreover, taking advantage of recent redshift surveys along CL0939+4713 line of sight, separate diagrams comparing the properties of galaxies belonging to the cluster and to its close projected field are presented. Possible evolutionary effects in the UV from z~0.4 to present time are investigated by comparing the restframe [mid-UV-Optical] colors of galaxies in CL0939+4713 with balloon-borne data of the Coma cluster, as well as by resorting to suitable galaxy evolution models. Finally, current attempts to constrain the epoch of the UV-upturn onset in evolved populations by means of HST UV observations are discussed.

The R-band isophotal map of the Sa galaxy NGC 4698 shows that the inner region of the bulge structure is elongated perpendicularly to the major axis of the disk, this is also true for the outer parts of the bulge if a parametric photometric decomposition is adopted. At the same time the stellar component is characterized by an inner velocity gradient and a central zero-velocity plateau along the minor and major axis of the disk respectively. This remarkable geometric and kinematic decoupling suggests that a second event occurred in the formation history of this galaxy.

We report the case of the geometrical and kinematical decoupling between the bulge and the disk of the Sa galaxy NGC 4698. The R-band isophotal map of this spiral shows that the bulge structure is elongated perpendicularly to the major axis of the disk. At the same time a central stellar velocity gradient is found along the major axis of the bulge. We also present the Sa NGC 4672 as good candidate of a spiral hosting a bulge and a disk orthogonally decoupled with respect to one other. This decoupling of the two fundamental components of a visible galaxy suggests that the disk could represent a second event in the history of early-type spirals.

The complex phenomenology shown by the UV-bright, variable spike first detected with the Hubble Space Telescope (HST) at the center of the otherwise normal galaxy NGC 4552 is further investigated with both HST imaging (FOC) and spectroscopy (FOS). HST/FOC images taken in 1991, 1993, and 1996 in the near UV have been analyzed in a homogeneous fashion, showing that the central spike has brightened by a factor ~4.5 between 1991 and 1993, and has decreased its luminosity by a factor ~2.0 between 1993 and 1996. FOS spectroscopy extending from the near UV to the red side of the optical spectrum reveals a strong UV continuum over the spectrum of the underlying galaxy, along with several emission lines in both the UV and the optical ranges. In spite of the low luminosity of the UV continuum of the spike (~3*10^5 Lsolar), the spike is definitely placed among AGNs by current diagnostics based on the emission line intensity ratios, being just on the borderline between Seyferts and LINERs. Line profiles are very broad, and both permitted and forbidden lines are best modelled with a combination of broad and narrow components, with FWHM of ~3000 km s^-1 and ~700 km s^-1, respectively. This evidence argues for the variable central spike being produced by a modest accretion event onto a central massive black hole (BH), with the accreted material having possibly being stripped from a a star in a close fly by with the BH. The 1996 broad Halpha luminosity of this mini-AGN is ~5.6*10^37 erg s^-1, about a factor of two less than that of the nucleus of NGC 4395, heretofore considered to be the faintest known AGN.

FOS spectra and FOC photometry of two centrally located, UV-bright spikes in the elliptical galaxy NGC 4552 and the bulge-dominated early spiral NGC 2681, are presented. These spectra reveal that such point-like UV sources detected by means of HST within a relatively large fraction ~15% of bulges can be related to radically different phenomena. While the UV unresolved emission in NGC 4552 represents a transient event likely induced by an accretion event onto a supermassive black hole, the spike seen at the center of NGC 2681 is not variable and it is stellar in nature.

We have obtained high resolution UV images with the HST+WFPC2 of the central region of the dwarf elliptical galaxy NGC 205. Our images reveal that many of the hot UV stars previously detected and studied from the ground are actually multiple systems, open clusters and star associations. We have performed photometry of two such clusters and we find our data are consistent with stellar ages of 50 and 100 Myr, respectively. From the number of massive stars in NGC 205 we estimate that the star formation episode in this galaxy has turned ~1,000 Msolar of gas into stars over the last 100 Myr.

In this paper we demonstrate the possibility of inferring the presence of Keplerian gaseous disks using optical ground-based telescopes properly equipped. We have modeled the peculiar bidimensional shape of the emission lines in a sample of five S0-Sa galaxies as due to the motion of a gaseous disk rotating in the combined potential of a central point-like mass and of an extended stellar disk. The value of the central mass concentration estimated for four galaxies of the sample (NGC 2179, NGC 4343, NGC 4435 and NGC 4459) is ~10^9 Msolar. For the remaining galaxy NGC 5064 an upper limit of 5*10^7 Msolar is estimated.