[Abridged] We present the apparent stellar angular momentum of 300 galaxies across the Hubble sequence, using integral-field spectroscopic data from the CALIFA survey. Adopting the same $\lambda_\mathrm{R}$ parameter previously used to distinguish between slow and fast rotating early-type (elliptical and lenticular) galaxies, we show that spiral galaxies as expected are almost all fast rotators. Given the extent of our data, we provide relations for $\lambda_\mathrm{R}$ measured in different apertures, including conversions to long-slit 1D apertures. Our sample displays a wide range of $\lambda_\mathrm{Re}$ values, consistent with previous IFS studies. The fastest rotators are dominated by relatively massive and highly star-forming Sb galaxies, which preferentially reside in the main star-forming sequence. These galaxies reach $\lambda_\mathrm{Re}$ values of $\sim$0.85, are the largest galaxies at a given mass, and display some of the strongest stellar population gradients. Compared to the population of S0 galaxies, our findings suggest that fading may not be the dominant mechanism transforming spirals into lenticulars. Interestingly, we find that $\lambda_\mathrm{Re}$ decreases for late-type Sc and Sd spiral galaxies, with values than in occasions puts them in the slow-rotator regime. While for some of them this can be explained by their irregular morphologies and/or face-on configurations, others are edge-on systems with no signs of significant dust obscuration. The latter are typically at the low-mass end, but this does not explain their location in the classical ($V/\sigma$,$\varepsilon$) and ($\lambda_\mathrm{Re}$,$\varepsilon$) diagrams. Our initial investigations, based on dynamical models, suggest that these are dynamically hot disks, probably influenced by the observed important fraction of dark matter within R$_\mathrm{e}$.

Einstein's theory of gravity, General Relativity, has been precisely tested on Solar System scales, but the long-range nature of gravity is still poorly constrained. The nearby strong gravitational lens, ESO 325-G004, provides a laboratory to probe the weak-field regime of gravity and measure the spatial curvature generated per unit mass, $\gamma$. By reconstructing the observed light profile of the lensed arcs and the observed spatially resolved stellar kinematics with a single self-consistent model, we conclude that $\gamma = 0.97 \pm 0.09$ at 68% confidence. Our result is consistent with the prediction of 1 from General Relativity and provides a strong extragalactic constraint on the weak-field metric of gravity.

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.

Galaxy formation entails the hierarchical assembly of mass, along with the condensation of baryons and the ensuing, self-regulating star formation. The stars form a collisionless system whose orbit distribution retains dynamical memory that can constrain a galaxy's formation history. The ordered-rotation dominated orbits with near maximum circularity $\lambda_z \simeq1$ and the random-motion dominated orbits with low circularity $\lambda_z \simeq0$ are called kinematically cold and kinematically hot, respectively. The fraction of stars on `cold' orbits, compared to the fraction of stars on `hot' orbits, speaks directly to the quiescence or violence of the galaxies' formation histories. Here we present such orbit distributions, derived from stellar kinematic maps via orbit-based modelling for a well defined, large sample of 300 nearby galaxies. The sample, drawn from the CALIFA survey, includes the main morphological galaxy types and spans the total stellar mass range from $10^{8.7}$ to $10^{11.9}$ solar masses. Our analysis derives the orbit-circularity distribution as a function of galaxy mass, $p(\lambda_z~|~M_\star)$, and its volume-averaged total distribution, $p(\lambda_z)$. We find that across most of the considered mass range and across morphological types, there are more stars on `warm' orbits defined as $0.25\le \lambda_z \le 0.8$ than on either `cold' or `hot' orbits. This orbit-based "Hubble diagram" provides a benchmark for galaxy formation simulations in a cosmological context.

Based on MUSE integral-field data we present evidence for a radial variation at the low-mass end of the stellar initial-mass function (IMF) in the central regions of the giant early-type galaxy NGC4486 (M87). We used state-of-the-art stellar population models and the observed strength of various IMF-sensitive absorption-line features to solve for the best low-mass tapered "bimodal" form of the IMF, while accounting also for radial variations in stellar metallicity, the overall $\alpha$-elements abundance, and the abundance of individual elements such as Ti, O, Na and Ca. Our analysis reveals a strong IMF gradient in M87, corresponding to an exceeding fraction of low-mass stars compared to the case of the Milky Way toward the center of M87 that drops to nearly Milky-way levels by 0.4 $R_e$. This IMF gradient is found to correlate well with both the radial profile for stellar metallicity and for $\alpha$-elements abundance but not with stellar velocity dispersion. Such IMF variations correspond to over a factor two increase in stellar mass-to-light M/L ratio compared to the case of a Milky-way like IMF, consistent with other investigations into IMF gradients in early-type galaxies, including recent dynamical constraints on M/L radial variations in M87 by Oldham & Auger. In addition to constraining the IMF in M87 we also looked into the abundance of Sodium, which turned up to be super-Solar over the entire radial range of our MUSE observations and to exhibit a considerable negative gradient. These findings suggest an additional role of metallicity in boosting the Na-yields in the central, metal-rich regions of M87 during its early and brief star-formation history. Our work adds the case of M87 to the few objects that as of today have radial constraints on their IMF or [Na/Fe] abundance, while also illustrating the accuracy that MUSE could bring to this kind of investigations.

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.

Schwarzschild orbit-based dynamical models are widely used to uncover the internal dynamics of early-type galaxies and globular clusters. Here we present for the first time the Schwarzschild models of late-type galaxies: an SBb galaxy NGC 4210 and an S0 galaxy NGC 6278 from the CALIFA survey. The mass profiles within $2\,R_e$ are constrained well with $1\sigma$ statistical error of $\sim 10\%$. The luminous and dark mass can be disentangled with uncertainties of $\sim 20\%$ and $\sim 50\%$ respectively. From $R_e$ to $2\,R_e$, the dark matter fraction increases from $14\pm10\%$ to $18\pm10\%$ for NGC 4210 and from $15\pm10\%$ to $30\pm20\%$ for NGC 6278. The velocity anisotropy profiles of both $\sigma_r/\sigma_t$ and $\sigma_z/\sigma_R$ are well constrained. The inferred internal orbital distributions reveal clear substructures. The orbits are naturally separated into three components: a cold component with near circular orbits, a hot component with near radial orbits, and a warm component in between. The photometrically-identified exponential disks are predominantly made up of cold orbits only beyond $\sim 1\,R_e$, while they are constructed mainly with the warm orbits inside. Our dynamical hot components are concentrated in the inner regions, similar to the photometrically-identified bulges. The reliability of the results, especially the orbit distribution, are verified by applying the model to mock data.

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 the detection of supermassive black holes (BHs) in two Virgo ultracompact dwarf galaxies (UCDs), VUCD3 and M59cO. We use adaptive optics assisted data from the Gemini/NIFS instrument to derive radial velocity dispersion profiles for both objects. Mass models for the two UCDs are created using multi-band Hubble Space Telescope (HST) imaging, including the modeling of mild color gradients seen in both objects. We then find a best-fit stellar mass-to-light ratio ($M/L$) and BH mass by combining the kinematic data and the deprojected stellar mass profile using Jeans Anisotropic Models (JAM). Assuming axisymmetric isotropic Jeans models, we detect BHs in both objects with masses of $4.4^{+2.5}_{-3.0} \times 10^6$ $M_{\odot}$ in VUCD3 and $5.8^{+2.5}_{-2.8} \times 10^6$ $M_{\odot}$ in M59cO (3$\sigma$ uncertainties). The BH mass is degenerate with the anisotropy parameter, $\beta_z$; for the data to be consistent with no BH requires $\beta_z = 0.4$ and $\beta_z = 0.6$ for VUCD3 and M59cO, respectively. Comparing these values with nuclear star clusters shows that while it is possible that these UCDs are highly radially anisotropic, it seems unlikely. These detections constitute the second and third UCDs known to host supermassive BHs. They both have a high fraction of their total mass in their BH; $\sim$13% for VUCD3 and $\sim$18% for M59cO. They also have low best-fit stellar $M/L$s, supporting the proposed scenario that most massive UCDs host high mass fraction BHs. The properties of the BHs and UCDs are consistent with both objects being the tidally stripped remnants of $\sim$10$^9$ M$_\odot$ galaxies.

Dedicated photometric and spectroscopic surveys have provided unambiguous evidence for a strong stellar mass-size evolution of galaxies within the last 10 Gyr. The likely progenitors of today's most massive galaxies are remarkably small, disky, passive and have already assembled much of their stellar mass at redshift z=2. An in-depth analysis of these objects, however, is currently not feasible due to the lack of high-quality, spatially-resolved photometric and spectroscopic data. In this paper, we present a sample of nearby compact elliptical galaxies (CEGs), which bear resemblance to the massive and quiescent galaxy population at earlier times. Hubble Space Telescope (HST) and wide-field integral field unit (IFU) data have been obtained, and are used to constrain orbit-based dynamical models and stellar population synthesis (SPS) fits, to unravel their structural and dynamical properties. We first show that our galaxies are outliers in the present-day stellar mass-size relation. They are, however, consistent with the mass-size relation of compact, massive and quiescent galaxies at redshift z=2. The compact sizes of our nearby galaxies imply high central stellar mass surface densities, which are also in agreement with the massive galaxy population at higher redshift, hinting at strong dissipational processes during their formation. Corroborating evidence for a largely passive evolution within the last 10 Gyr is provided by their orbital distribution as well as their stellar populations, which are difficult to reconcile with a very active (major) merging history. This all supports that we can use nearby CEGs as local analogues of the high-redshift, massive and quiescent galaxy population, thus providing additional constraints for models of galaxy formation and evolution.

Mrk 1216 is a nearby, early-type galaxy with a small effective radius of 2.8 kpc and a large stellar velocity dispersion of 308 km/s for its K-band luminosity of 1.4x10^11 L_sun. Using integral-field spectroscopy assisted by adaptive optics from Gemini North, we measure spatially resolved stellar kinematics within ~450 pc of the galaxy nucleus. The galaxy exhibits regular rotation with velocities of \pm 180 km/s and a sharply peaked velocity dispersion profile that reaches 425 km/s at the center. We fit axisymmetric, orbit-based dynamical models to the combination of these high angular resolution kinematics, large-scale kinematics extending to roughly three effective radii, and Hubble Space Telescope imaging, resulting in a constraint of the mass of the central black hole in Mrk 1216. After exploring several possible sources of systematics that commonly affect stellar-dynamical black hole mass measurements, we find a black hole mass of (4.9\pm1.7)x10^9 M_sun and a H-band stellar mass-to-light ratio of 1.3\pm0.4 M_sun/L_sun (1-sigma uncertainties). Mrk 1216 is consistent with the local black hole mass - stellar velocity dispersion relation, but is a factor of ~5-10 larger than expectations from the black hole mass - bulge luminosity and black hole mass - bulge mass correlations when conservatively using the galaxy's total luminosity or stellar mass. This behavior is quite similar to the extensively studied compact galaxy NGC 1277. Resembling the z~2 quiescent galaxies, Mrk 1216 may be a passively evolved descendant, and perhaps reflects a previous era when galaxies contained over-massive black holes relative to their bulge luminosities/masses, and the growth of host galaxies had yet to catch up.

We present the stellar kinematic maps of a large sample of galaxies from the integral-field spectroscopic survey CALIFA. The sample comprises 300 galaxies displaying a wide range of morphologies across the Hubble sequence, from ellipticals to late-type spirals. This dataset allows us to homogeneously extract stellar kinematics up to several effective radii. In this paper, we describe the level of completeness of this subset of galaxies with respect to the full CALIFA sample, as well as the virtues and limitations of the kinematic extraction compared to other well-known integral-field surveys. In addition, we provide averaged integrated velocity dispersion radial profiles for different galaxy types, which are particularly useful to apply aperture corrections for single aperture measurements or poorly resolved stellar kinematics of high-redshift sources. The work presented in this paper sets the basis for the study of more general properties of galaxies that will be explored in subsequent papers of the survey.

Coevolution between supermassive black holes (BHs) and their host galaxies is universally adopted in models for galaxy formation. In the absence of feedback from active galactic nuclei, simulated massive galaxies keep forming stars in the local Universe. From an observational point of view, however, such coevolution remains unclear. We present a stellar population analysis of galaxies with direct BH mass measurements and the BH mass-{\sigma} relation as a working framework. We find that over-massive BH galaxies, i.e., galaxies lying above the best-fitting BH mass-{\sigma} line, tend to be older and more {\alpha}-element enhanced than under-massive BH galaxies. The scatter in the BH mass-{\sigma}-[{\alpha}/Fe] plane is significantly lower than in the standard BH mass-{\sigma} relation. We interpret this trend as an imprint of active galactic nucleus feedback on the star formation histories of massive galaxies.

We present a study of the excitation conditions and metallicity of ionized gas (Z_gas) in eight nearby barred and unbarred spiral galaxies from the VIRUS-P Exploration of Nearby Galaxies (VENGA) survey, which provides high spatial sampling and resolution (median ~ 387 pc), large coverage from the bulge to outer disc, broad wavelength range (3600-6800 Ang.), and medium spectral resolution (~ 120 km s^-1 at 5000 Ang.). Our results are: (1) We present high resolution gas excitation maps to differentiate between regions with excitation typical of Seyfert, LINER, or recent star formation. We find LINER-type excitation at large distances (3-10 kpc) from the centre, and associate this excitation with diffuse ionized gas (DIG). (2) After excluding spaxels dominated by Seyfert, LINER, or DIG, we produce maps with the best spatial resolution and sampling to date of the ionization parameter q, star formation rate, and Z_gas using common strong line diagnostics. We find that isolated barred and unbarred spirals exhibit similarly shallow Z_gas profiles from the inner kpc out to large radii (7-10 kpc or 0.5-1.0 R_25). This implies that if profiles had steeper gradients at earlier epochs, then the present-day bar is not the primary driver flattening gradients over time. This result contradicts earlier claims, but agrees with recent IFU studies. (3) The Z_gas gradients in our z ~ 0 massive spirals are markedly shallower, by ~ 0.2 dex kpc^-1, than published gradients for lensed lower mass galaxies at z ~ 1.5-2.0. Cosmologically-motivated hydrodynamical simulations best match this inferred evolution, but the match is sensitive to adopted stellar feedback prescriptions.

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.

In this letter we describe how we use stellar dynamics information to constrain the shape of the stellar IMF in a sample of 27 early-type galaxies from the CALIFA survey. We obtain dynamical and stellar mass-to-light ratios, $\Upsilon_\mathrm{dyn}$ and $\Upsilon_{\ast}$, over a homogenous aperture of 0.5~$R_{e}$. We use the constraint $\Upsilon_\mathrm{dyn} \ge \Upsilon_{\ast}$ to test two IMF shapes within the framework of the extended MILES stellar population models. We rule out a single power law IMF shape for 75% of the galaxies in our sample. Conversely, we find that a double power law IMF shape with a varying high-mass end slope is compatible (within 1$\sigma$) with 95% of the galaxies. We also show that dynamical and stellar IMF mismatch factors give consistent results for the systematic variation of the IMF in these galaxies.

According to the Virial Theorem, all gravitational systems in equilibrium sit on a plane in the 3D parameter space defined by their mass, size and second moment of the velocity tensor. While these quantities cannot be directly observed, there are suitable proxies: the luminosity Lk, half-light radius Re and dispersion sigma_e. These proxies indeed lie on a very tight Fundamental Plane (FP). How do the black holes in the centers of galaxies relate to the FP? Their masses are known to exhibit no strong correlation with total galaxy mass, but they do correlate weakly with bulge mass (when present), and extremely well with the velocity dispersion through the Mbh = sigma_e^5.4 relation. These facts together imply that a tight plane must also exist defined by black hole mass, total galaxy mass and size. Here I show that this is indeed the case using a heterogeneous set of 230 black holes. The sample includes BHs from zero to 10 billion solar masses and host galaxies ranging from low surface brightness dwarfs, through bulge-less disks, to brightest cluster galaxies. The resulting BH-size-luminosity Mbh=(Lk/Re)^3.8 has the same amount of scatter as the M-sigma relation and is aligned with the galaxy FP, such that it is just a re-projection of sigma. The inferred BH-size-mass relation is Mbh=(M_star/Re)^2.9. These relationships are universal and extend to galaxies without bulges. This implies that the black hole is primarily correlated with its global velocity dispersion and not with the properties of the bulge. I show that the classical bulge--mass relation is a projection of the M-sigma relation. When the velocity dispersion cannot be measured (at high-z or low dispersions), the BH-size-mass relation should be used as a proxy for black hole mass in favor of just galaxy or bulge mass.

This paper describes the Third Public Data Release (DR3) of the Calar Alto Legacy Integral Field Area (CALIFA) survey. Science-grade quality data for 667 galaxies are made public, including the 200 galaxies of the Second Public Data Release (DR2). Data were obtained with the integral-field spectrograph PMAS/PPak mounted on the 3.5m telescope at the Calar Alto Observatory. Three different spectral setups are available, i) a low-resolution V500 setup covering the wavelength range 3749-7500 AA (4240-7140 AA unvignetted) with a spectral resolution of 6.0 AA (FWHM), for 646 galaxies, ii) a medium-resolution V1200 setup covering the wavelength range 3650-4840 AA (3650-4620 AA unvignetted) with a spectral resolution of 2.3 AA (FWHM), for 484 galaxies, and iii) the combination of the cubes from both setups (called COMBO), with a spectral resolution of 6.0 AA and a wavelength range between 3700-7500 AA (3700-7140 AA unvignetted), for 446 galaxies. The Main Sample, selected and observed according to the CALIFA survey strategy covers a redshift range between 0.005 and 0.03, spans the color-magnitude diagram and probes a wide range of stellar mass, ionization conditions, and morphological types. The Extension Sample covers several types of galaxies that are rare in the overall galaxy population and therefore not numerous or absent in the CALIFA Main Sample. All the cubes in the data release were processed using the latest pipeline, which includes improved versions of the calibration frames and an even further improved im- age reconstruction quality. In total, the third data release contains 1576 datacubes, including ~1.5 million independent spectra. It is available at http://califa.caha.es/DR3.

Megamaser disks provide the most precise and accurate extragalactic supermassive black hole masses. Here we describe a search for megamasers in nearby galaxies using the Green Bank Telescope (GBT). We focus on galaxies where we believe that we can resolve the gravitational sphere of influence of the black hole and derive a stellar or gas dynamical measurement with optical or NIR observations. Since there are only a handful of super massive black holes (SMBH) that have direct black hole mass measurements from more than one method, even a single galaxy with a megamaser disk and a stellar dynamical black hole mass would provide necessary checks on the stellar dynamical methods. We targeted 87 objects from the Hobby-Eberly Telescope Massive Galaxy Survey, and detected no new maser disks. Most of the targeted objects are elliptical galaxies with typical stellar velocity dispersions of 250 km/s and distances within 130 Mpc. We discuss the implications of our non-detections, whether they imply a threshold X-ray luminosity required for masing, or possibly reflect the difficulty of maintaining a masing disk around much more massive (>10^8 Msun) black holes at low Eddington ratio. Given the power of maser disks at probing black hole accretion and demographics, we suggest that future maser searches should endeavour to remove remaining sample biases, in order to sort out the importance of these covariant effects.

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.

NGC 1042 is a late type bulgeless disk galaxy which hosts a low luminosity Active Galactic Nuclei (AGN) coincident with a massive nuclear star cluster. In this paper, we present the integral-field-spectroscopy (IFS) studies of this galaxy, based on the data obtained with the Mitchell spectrograph on the 2.7 meter Harlan J. Smith telescope. In the central $100\textrm{-}300\ \mathrm{pc}$ region of NGC 1042, we find a circumnuclear ring structure of gas with enhanced ionization, which we suggest is mainly induced by shocks. Combining with the harmonic decomposition analysis of the velocity field of the ionized gas, we propose that the shocked gas is the result of gas inflow driven by the inner spiral arms. The inflow velocity is $\sim 32\pm10\ \mathrm{km}\ \mathrm{s}^{-1}$ and the estimated mass inflow rate is $\sim 1.1\pm0.3 \times 10^{-3}\ \mathrm{M}_{\odot}\ \mathrm{yr}^{-1}$. The mass inflow rate is about one hundred times the blackhole's mass accretion rate ($\sim 1.4 \times 10^{-5}\ \mathrm{M}_{\odot}\ \mathrm{ yr}^{-1}$), and slightly larger than the star formation rate in the nuclear star cluster ($7.94 \times 10^{-4}\ \mathrm{M}_{\odot}\ \mathrm{yr}^{-1}$), implying that the inflow material is enough to feed both the AGN activity and the star formation in the nuclear star cluster. Our study highlights that secular evolution can be important in late-type unbarred galaxies like NGC 1042.

The nearby lenticular galaxy NGC 1277 is thought to host one of the largest black holes known, however the black hole mass measurement is based on low spatial resolution spectroscopy. In this paper, we present Gemini Near-infrared Integral Field Spectrometer observations assisted by adaptive optics. We map out the galaxy's stellar kinematics within ~440 pc of the nucleus with an angular resolution that allows us to probe well within the region where the potential from the black hole dominates. We find that the stellar velocity dispersion rises dramatically, reaching ~550 km/s at the center. Through orbit-based, stellar-dynamical models we obtain a black hole mass of (4.9 \pm 1.6) x 10^9 Msun (1-sigma uncertainties). Although the black hole mass measurement is smaller by a factor of ~3 compared to previous claims based on large-scale kinematics, NGC 1277 does indeed contain one of the most massive black holes detected to date, and the black hole mass is an order of magnitude larger than expectations from the empirical relation between black hole mass and galaxy luminosity. Given the galaxy's similarities to the higher redshift (z~2) massive quiescent galaxies, NGC 1277 could be a relic, passively evolving since that period. A population of local analogs to the higher redshift quiescent galaxies that also contain over-massive black holes may suggest that black hole growth precedes that of the host galaxy.

We investigate the compact, early-type galaxy NGC 1281 with integral field unit observations to map the stellar LOSVD out to 5 effective radii and construct orbit-based dynamical models to constrain its dark and luminous matter content. Under the assumption of mass-follows-light, the H-band stellar mass-to-light ratio (M/L) is {\Upsilon} = 2.7(+-0.1) {\Upsilon}_{sun}, higher than expected from our stellar population synthesis fits with either a canonical Kroupa ({\Upsilon} = 1.3 {\Upsilon}_{sun}) or Salpeter ({\Upsilon} = 1.7 {\Upsilon}_{sun}) stellar initial mass function. Such models also cannot reproduce the details of the LOSVD. Models with a dark halo recover the kinematics well and indicate that NGC 1281 is dark matter dominated, making up ~ 90 per cent of the total enclosed mass within the kinematic bounds. Parameterised as a spherical NFW profile, the dark halo mass is 11.5 < log(M_{DM}/M_{sun}) < 11.8 and the stellar M/L is 0.6 < {\Upsilon} < 1.1. However, this stellar M/L is lower than predicted by its old stellar population. Moreover, the halo mass within the kinematic extent is ten times larger than expected based on {\Lambda}CDM predictions, and an extrapolation yields cluster sized dark halo masses. Adopting {\Upsilon} = 1.7 {\Upsilon}_{sun} yields more moderate dark halo virial masses, but these models fit the kinematics worse. A non-NFW model might solve the discrepancy between the unphysical consequences of the best-fitting dynamical models and models based on more reasonable assumptions for the dark halo and stellar mass-to-light ratio, which are disfavoured according to our parameter estimation.

We present a dynamical analysis to infer the structural parameters and properties of the two nearby, compact, high velocity dispersion galaxies MRK1216 & NGC1277. Combining deep HST imaging, wide-field IFU stellar kinematics, and complementary long-slit spectroscopic data out to 3 R_e, we construct orbit-based models to constrain their black hole masses, dark matter content and stellar mass-to-light ratios. We obtain a black hole mass of log(Mbh/Msun) = 10.1(+0.1/-0.2) for NGC1277 and an upper limit of log(Mbh/Msun) = 10.0 for MRK1216, within 99.7 per cent confidence. The stellar mass-to-light ratios span a range of Upsilon_V = 6.5(+1.5/-1.5) in NGC1277 and Upsilon_H = 1.8(+0.5/-0.8) in MRK1216 and are in good agreement with SSP models of a single power-law Salpeter IMF. Even though our models do not place strong constraints on the dark halo parameters, they suggest that dark matter is a necessary ingredient in MRK1216, with a dark matter contribution of 22(+30/-20) per cent to the total mass budget within 1 R_e. NGC1277, on the other hand, can be reproduced without the need for a dark halo, and a maximal dark matter fraction of 13 per cent within the same radial extent. In addition, we investigate the orbital structures of both galaxies, which are rotationally supported and consistent with photometric multi-S\'ersic decompositions, indicating that these compact objects do not host classical, non-rotating bulges formed during recent (z <= 2) dissipative events or through violent relaxation. Finally, both MRK 1216 and NGC 1277 are anisotropic, with a global anisotropy parameter delta of 0.33 and 0.58, respectively. While MRK 1216 follows the trend of fast-rotating, oblate galaxies with a flattened velocity dispersion tensor in the meridional plane of the order of beta_z = delta, NGC 1277 is highly tangentially anisotropic and seems to belong kinematically to a distinct class of objects.

Located in the Perseus cluster, NGC 1271 is an early-type galaxy with a small effective radius of 2.2 kpc and a large stellar velocity dispersion of 276 km/s for its K-band luminosity of 8.9x10^{10} L_sun. We present a mass measurement for the black hole in this compact, high-dispersion galaxy using observations from the integral field spectrograph NIFS on the Gemini North telescope assisted by laser guide star adaptive optics, large-scale integral field unit observations with PPAK at the Calar Alto Observatory, and Hubble Space Telescope WFC3 imaging observations. We are able to map out the stellar kinematics on small spatial scales, within the black hole sphere of influence, and on large scales that extend out to four times the galaxy's effective radius. We find that the galaxy is rapidly rotating and exhibits a sharp rise in the velocity dispersion. Through the use of orbit-based stellar dynamical models, we determine that the black hole has a mass of (3.0^{+1.0}_{-1.1}) x 10^9 M_sun and the H-band stellar mass-to-light ratio is 1.40^{+0.13}_{-0.11} M_sun/L_sun (1-sigma uncertainties). NGC 1271 occupies the sparsely-populated upper end of the black hole mass distribution, but is very different from the Brightest Cluster Galaxies (BCGs) and giant elliptical galaxies that are expected to host the most massive black holes. Interestingly, the black hole mass is an order of magnitude larger than expectations based on the galaxy's bulge luminosity, but is consistent with the mass predicted using the galaxy's bulge stellar velocity dispersion. More compact, high-dispersion galaxies need to be studied using high spatial resolution observations to securely determine black hole masses, as there could be systematic differences in the black hole scaling relations between these types of galaxies and the BCGs/giant ellipticals, thereby implying different pathways for black hole and galaxy growth.

We study a sample of eight massive galaxies that are extreme outliers (3-5$\sigma$) in the M$_{\bullet}$-M$_\mathrm{bulge}$ local scaling relation. Two of these galaxies are confirmed to host extremely large super massive black holes (SMBHs), whereas the virial mass estimates for the other six are also consistent with having abnormally large SMBHs. From the analysis of their star formation histories and their structural properties we find that all these extreme outliers can be considered as relic galaxies from the early (z$\sim$2) Universe: i.e. they are compact (R$_{\mathrm{e}}$$<$2 kpc) and have purely old stellar populations (t$\gtrsim$10 Gyr). In order to explain the nature of such deviations from the local relations, we propose a scenario in which the hosts of these \"uber-massive SMBHs are galaxies that have followed a different evolutionary path than the two-phase growth channel assumed for massive galaxies. Once the SMBH and the core of the galaxy are formed at z$\sim$2, the galaxy skips the second phase, remaining structurally untouched and without further mass and size increase. We show that if the outliers had followed the normal evolutionary path by growing in size via merger activity, the expected (mild) growth in mass would place them closer to the observed local relations. Our results suggest that the SMBH growth epoch for the most massive galaxies stopped $\sim$10Gyr ago.

There has been significant controversy over the mechanisms responsible for forming compact stellar systems like ultra compact dwarfs (UCDs), with suggestions that UCDs are simply the high mass extension of the globular cluster (GC) population, or alternatively, the liberated nuclei of galaxies tidally stripped by larger companions. Definitive examples of UCDs formed by either route have been difficult to find, with only a handful of persuasive examples of stripped-nucleus type UCDs being known. In this paper we present very deep Gemini/GMOS spectroscopic observations of the suspected stripped nucleus UCD NGC 4546-UCD1 taken in good seeing conditions (< 0.7"). With these data we examine the spatially resolved kinematics and star formation history of this unusual object. We find no evidence of a rise in the central velocity dispersion of the UCD, suggesting that this UCD lacks a massive central black hole like those found in some other compact stellar systems, a conclusion confirmed by detailed dynamical modelling. Finally we are able to use our extremely high signal to noise spectrum to detect a temporally extended star formation history for this UCD. We find that the UCD was forming stars since the earliest epochs until at least 1-2 Gyr ago. Taken together these observations confirm that NGC 4546-UCD1 is the remnant nucleus of a nucleated dwarf galaxy that was tidally destroyed by NGC 4546 within the last 1-2 Gyr.

Aims. We aim at finding candidates of potential survivors of high-redshift compact galaxies in SDSS, as targets for more detailed follow-up observations. Methods. From the virial theorem it is expected that for a given mass, compact galaxies have stellar velocity dispersion higher than the mean due to their smaller sizes. Therefore velocity dispersion coupled with size (or mass) is an appropriate method to select relics, independent of the stellar population properties. Based on these consideration we design a set of criteria using distribution of early-type galaxies from SDSS on the log$_{10}$(R$_{0}$)-log$_{10}$($\sigma_{0}$) plane to find the most extreme objects on it. Results. We find 76 galaxies at 0.05 < z < 0.2, which have properties similar to the typical quiescent galaxies at high redshift. We study how well these galaxies fit on well-known local universe relations of early-type galaxies such as the fundamental plane, the red sequence or mass-size relations. As expected from the selection criteria, the candidates are located in an extreme corner of mass-size plane. However, they do not extend as deeply into the so-called zone of exclusion as some of the high-redshift compact galaxies ('red nuggets') found at high redshift, being a factor 2-3 less massive at a given intrinsic scale size. Our candidates are systematically offset from scaling relations of average early-type galaxies, while being in the mass-size range expected for passive evolution of the red nuggets from their high redshift to the present. Conclusions. The 76 selected candidates form a well suited set of objects for further follow-up observations. We argue that selecting a high velocity dispersion is the best way to find analogues of compact high redshift galaxies in the local universe.

We have conducted an optical long-slit spectroscopic survey of 1022 galaxies using the 10m Hobby-Eberly Telescope (HET) at McDonald Observatory. The main goal of the HET Massive Galaxy Survey (HETMGS) is to find nearby galaxies that are suitable for black hole mass measurements. In order to measure accurately the black hole mass, one should kinematically resolve the region where the black hole dominates the gravitational potential. For most galaxies, this region is much less than an arcsecond. Thus, black hole masses are best measured in nearby galaxies with telescopes that obtain high-spatial resolution. The HETMGS focuses on those galaxies predicted to have the largest sphere-of-influence, based on published stellar velocity dispersions or the galaxy fundamental plane. To ensure coverage over galaxy types, the survey targets those galaxies across a face-on projection of the fundamental plane. We present the sample selection and resulting data products from the long-slit observations, including central stellar kinematics and emission line ratios. The full dataset, including spectra and resolved kinematics, is available online. Additionally, we show that the current crop of black hole masses are highly biased towards dense galaxies and that especially large disks and low dispersion galaxies are under-represented. This survey provides the necessary groundwork for future systematic black hole mass measurement campaigns.

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.

Ultracompact dwarf galaxies (UCDs) are among the densest stellar systems in the universe. These systems have masses up to 200 million solar masses, but half light radii of just 3-50 parsecs. Dynamical mass estimates show that many UCDs are more massive than expected from their luminosity. It remains unclear whether these high dynamical mass estimates are due to the presence of supermassive black holes or result from a non-standard stellar initial mass function that causes the average stellar mass to be higher than expected. Here we present the detection of a supermassive black hole in a massive UCD. Adaptive optics kinematic data of M60-UCD1 show a central velocity dispersion peak above 100 km/s and modest rotation. Dynamical modeling of these data reveals the presence of a supermassive black hole with mass of 21 million solar masses. This is 15% of the object's total mass. The high black hole mass and mass fraction suggest that M60-UCD1 is the stripped nucleus of a galaxy. Our analysis also shows that M60-UCD1's stellar mass is consistent with its luminosity, implying many other UCDs may also host supermassive black holes. This suggests a substantial population of previously unnoticed supermassive black holes.

We present new constraints on the density profiles of dark matter (DM) halos in seven nearby dwarf galaxies from measurements of their integrated stellar light and gas kinematics. The gas kinematics of low mass galaxies frequently suggest that they contain constant density DM cores, while N-body simulations instead predict a cuspy profile. We present a data set of high resolution integral field spectroscopy on seven galaxies and measure the stellar and gas kinematics simultaneously. Using Jeans modeling on our full sample, we examine whether gas kinematics in general produce shallower density profiles than are derived from the stars. Although 2/7 galaxies show some localized differences in their rotation curves between the two tracers, estimates of the central logarithmic slope of the DM density profile, gamma, are generally robust. The mean and standard deviation of the logarithmic slope for the population are gamma=0.67+/-0.10 when measured in the stars and gamma=0.58+/-0.24 when measured in the gas. We also find that the halos are not under concentrated at the radii of half their maximum velocities. Finally, we search for correlations of the DM density profile with stellar velocity anisotropy and other baryonic properties. Two popular mechanisms to explain cored DM halos are an exotic DM component or feedback models that strongly couple the energy of supernovae into repeatedly driving out gas and dynamically heating the DM halos. We investigate correlations that may eventually be used to test models. We do not find a secondary parameter that strongly correlates with the central DM density slope, but we do find some weak correlations. Determining the importance of these correlations will require further model developments and larger observational samples. (Abridged)

We analyze line-of-sight velocity and proper motion data of stars in the Galactic globular cluster M15 using a new method to fit dynamical models to discrete kinematic data. Our fitting method maximizes the likelihood for individual stars and, as such, does not suffer the same loss of spatial and velocity information incurred when spatially binning data or measuring velocity moments. In this paper, we show that the radial variation in M15 of the mass-to-light ratio is consistent with previous estimates and theoretical predictions, which verifies our method. Our best-fitting axisymmetric Jeans models do include a central dark mass of $\sim2 \pm 1\cdot 10^3M_\odot$, which can be explained by a high concentration of stellar remnants at the cluster center. This paper shows that, from a technical point of view and with current computing power, spatial binning of data is no longer necessary. This not only leads to more accurate fits, but also avoids biased mass estimates due to the loss of resolution. Furthermore, we find that the mass concentration in M15 is significantly higher than previously measured, and is in close agreement with theoretical predictions for core-collapsed globular clusters without a central intermediate-mass black hole.

We present a new framework for modelling discrete kinematic data. Current techniques typically involve binning. Our approach works directly with the discrete data and uses maximum-likelihood methods to assess the probability of the dataset given model predictions. We avoid making hard cuts on the datasets by allowing for a contaminating population in our models. We apply our models to discrete proper-motion and line-of-sight-velocity data of Galactic globular cluster omega Centauri and find a mildly radial velocity anisotropy beta = 0.10 +/- 0.02, an inclination angle i = 50 +/- 1 deg, a V-band mass-to-light ratio Upsilon = 2.71 +/- 0.05 Msun/Lsun and a distance d = 4.59 +/- 0.08 kpc. All parameters are in agreement with previous studies, demonstrating the feasibility of our methods. We find that the models return lower distances and higher mass-to-light ratios than expected when we include proper motion stars with high errors or for which there is some blending. We believe this not a fault of our models but is instead due to underestimates or missing systematic uncertainties in the provided errors.

We have developed spherically symmetric dynamical models of dwarf spheroidal galaxies using Schwarzschild's orbit superposition method. This type of modelling yields constraints both on the total mass distribution (e.g. enclosed mass and scale radius) as well as on the orbital structure of the system (e.g. velocity anisotropy). This method is thus less prone to biases introduced by assumptions in comparison to the more commonly used Jeans modelling, and it allows us to derive the dark matter content in a robust way. Here we present our results for the Sculptor dwarf spheroidal galaxy, after testing our methods on mock data sets. We fit both the second and fourth velocity moment profile to break the mass-anisotropy degeneracy. For an NFW dark matter halo profile, we find that the mass of Sculptor within 1 kpc is M_1kpc = 1.03 \pm 0.07 x 10^8 Msol, and that its velocity anisotropy profile is tangentially biased and nearly constant with radius. The preferred concentration (c ~ 15) is low for its dark matter mass but consistent within the scatter found in N-body cosmological simulations. When we let the value of the central logarithmic slope \alpha{} vary, we find that the best-fit model has \alpha{} = 0, although an NFW cusp or shallower is consistent at 1 \sigma{} confidence level. On the other hand, very cuspy density profiles with logarithmic central slopes \alpha{} < -1.5 are strongly disfavoured for Sculptor.

We present orbit-based dynamical models and stellar population analysis of galaxy SDSS J151741.75-004217.6, a low-redshift (z=0.116) early-type galaxy (ETG) which, for its moderate luminosity, has an exceptionally high velocity dispersion. We aim to determine the central black hole mass (M_bh), the i-band stellar mass-to-light ratio, and the low-mass slope of the initial mass function (IMF). Combining constraints from HST imaging and longslit kinematic data with those from fitting the SDSS spectrum with stellar populations models of varying IMF, we show that this galaxy has a large fraction of low-mass stars, significantly higher than implied even by a Salpeter IMF. We exclude a Chabrier/Kroupa as well as a unimodal (i.e. single-segment) IMF, while a bimodal (low-mass tapered) shape is consistent with the dynamical constraints. Thereby, our study demonstrates that a very bottom-heavy IMF can exist even in an L* ETG. We place an upper limit of ~10^{10.5}M_sun on M_bh, which still leaves open the possibility of an extremely massive BH.

Recent observations have shown that compact nuclear star clusters (NSCs) are present in up to 80% of galaxies. However, detailed studies of their dynamical and chemical properties are confined mainly to spiral galaxy hosts, where they are more easily observed. In this paper we present our study of the NSC in FCC 277, a nucleated elliptical galaxy in the Fornax cluster. We use a combination of adaptive optics assisted near-infrared integral field spectroscopy, Hubble Space Telescope imaging, and literature long slit data. We show that while the NSC does not appear to rotate within our detection limit of ~6 km/s, rotation is detected at larger radii, where the isophotes appear to be disky, suggesting the presence of a nuclear disk. We also observe a distinct central velocity dispersion drop that is indicative of a dynamically cold rotating sub-system. Following the results of orbit-based dynamical modelling, co-rotating as well as counter-rotating stellar orbits are simultaneously needed to reproduce the observed kinematics. We find evidence for varying stellar populations, with the NSC and nuclear disk hosting younger and more metal rich stars than the main body of the galaxy. We argue that gas dissipation and some level of merging have likely played an important role in the formation of the nucleus of this intermediate-mass galaxy. This is in contrast to NSCs in low-mass early- type galaxies, which may have been formed primarily through the infall of star clusters.

We present the survey design, data reduction, and spectral fitting pipeline for the VIRUS-P Exploration of Nearby Galaxies (VENGA). VENGA is an integral field spectroscopic survey, which maps the disks of 30 nearby spiral galaxies. Targets span a wide range in Hubble type, star formation activity, morphology, and inclination. The VENGA data-cubes have 5.6'' FWHM spatial resolution, ~5A FWHM spectral resolution, sample the 3600A-6800A range, and cover large areas typically sampling galaxies out to ~0.7 R_25. These data-cubes can be used to produce 2D maps of the star formation rate, dust extinction, electron density, stellar population parameters, the kinematics and chemical abundances of both stars and ionized gas, and other physical quantities derived from the fitting of the stellar spectrum and the measurement of nebular emission lines. To exemplify our methods and the quality of the data, we present the VENGA data-cube on the face-on Sc galaxy NGC 628 (a.k.a. M 74). The VENGA observations of NGC 628 are described, as well as the construction of the data-cube, our spectral fitting method, and the fitting of the stellar and ionized gas velocity fields. We also propose a new method to measure the inclination of nearly face-on systems based on the matching of the stellar and gas rotation curves using asymmetric drift corrections. VENGA will measure relevant physical parameters across different environments within these galaxies, allowing a series of studies on star formation, structure assembly, stellar populations, chemical evolution, galactic feedback, nuclear activity, and the properties of the interstellar medium in massive disk galaxies.

We measure the radial profile of the 12CO(1-0) to H_2 conversion factor (Xco) in NGC 628. The H\alpha emission from the VENGA integral field spectroscopy is used to map the star formation rate surface density (\Sigma_{SFR}). We estimate the molecular gas surface density (\Sigma_{H2}) from \Sigma_{SFR} by inverting the molecular star formation law (SFL), and compare it to the CO intensity to measure Xco. We study the impact of systematic uncertainties by changing the slope of the SFL, using different SFR tracers (H\alpha vs. far-UV plus 24\mu m), and CO maps from different telescopes (single-dish and interferometers). The observed Xco profile is robust against these systematics, drops by a factor of 2 from R~7 kpc to the center of the galaxy, and is well fit by a gradient \Delta log(Xco)=0.06\pm0.02 dex kpc^-1. We study how changes in Xco follow changes in metallicity, gas density, and ionization parameter. Theoretical models show that the gradient in Xco can be explained by a combination of decreasing metallicity, and decreasing \Sigma_{H2} with radius. Photoelectric heating from the local UV radiation field appears to contribute to the decrease of Xco in higher density regions. Our results show that galactic environment plays an important role at setting the physical conditions in star forming regions, in particular the chemistry of carbon in molecular complexes, and the radiative transfer of CO emission. We caution against adopting a single Xco value when large changes in gas surface density or metallicity are present.

All massive galaxies likely have supermassive black holes at their centers, and the masses of the black holes are known to correlate with properties of the host galaxy bulge component. Several explanations have been proposed for the existence of these locally-established empirical relationships; they include the non-causal, statistical process of galaxy-galaxy merging, direct feedback between the black hole and its host galaxy, or galaxy-galaxy merging and the subsequent violent relaxation and dissipation. The empirical scaling relations are thus important for distinguishing between various theoretical models of galaxy evolution, and they further form the basis for all black hole mass measurements at large distances. In particular, observations have shown that the mass of the black hole is typically 0.1% of the stellar bulge mass of the galaxy. The small galaxy NGC4486B currently has the largest published fraction of its mass in a black hole at 11%. Here we report observations of the stellar kinematics of NGC 1277, which is a compact, disky galaxy with a mass of 1.2 x 10^11 Msun. From the data, we determine that the mass of the central black hole is 1.7 x 10^10 Msun, or 59% its bulge mass. Five other compact galaxies have properties similar to NGC 1277 and therefore may also contain over-sized black holes. It is not yet known if these galaxies represent a tail of a distribution, or if disk-dominated galaxies fail to follow the normal black hole mass scaling relations.

We present a new stellar dynamical mass measurement of the black hole in the nearby, S0 galaxy NGC 3998. By combining laser guide star adaptive optics observations obtained with the OH-Suppressing Infrared Imaging Spectrograph on the Keck II telescope with long-slit spectroscopy from the Hubble Space Telescope and the Keck I telescope, we map out the stellar kinematics on both small spatial scales, well within the black hole sphere of influence, and on large scales. We find that the galaxy is rapidly rotating and exhibits a sharp central peak in the velocity dispersion. Using the kinematics and the stellar luminosity density derived from imaging observations, we construct three-integral, orbit-based, triaxial stellar dynamical models. We find the black hole has a mass of M_BH = (8.1_{-1.9}^{+2.0}) x 10^8 M_sun, with an I-band stellar mass-to-light ratio of M/L = 5.0_{-0.4}^{+0.3} M_sun/L_sun (3-sigma uncertainties), and that the intrinsic shape of the galaxy is very round, but oblate. With the work presented here, NGC 3998 is now one of a very small number of galaxies for which both stellar and gas dynamical modeling have been used to measure the mass of the black hole. The stellar dynamical mass is nearly a factor of four larger than the previous gas dynamical black hole mass measurement. Given that this cross-check has so far only been attempted on a few galaxies with mixed results, carrying out similar studies in other objects is essential for quantifying the magnitude and distribution of the cosmic scatter in the black hole mass - host galaxy relations.

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.

We study the mass distribution in the late-type dwarf galaxy NGC 2976 through stellar kinematics obtained with the VIRUS-P integral-field spectrograph and anisotropic Jeans models as a test of cosmological simulations and baryonic processes that putatively alter small-scale structure. Previous measurements of the H-alpha emission-line kinematics have determined that the dark matter halo of NGC 2976 is most consistent with a cored density profile. We find that the stellar kinematics are best fit with a cuspy halo. Cored dark matter halo fits are only consistent with the stellar kinematics if the stellar mass-to-light ratio is significantly larger than that derived from stellar population synthesis, while the best-fitting cuspy model has no such conflict. The inferred mass distribution from a harmonic decomposition of the gaseous kinematics is inconsistent with that of the stellar kinematics. This difference is likely due to the gas disk not meeting the assumptions that underlie the analysis such as no pressure support, a constant kinematic axis, and planar orbits. By relaxing some of these assumptions, in particular the form of the kinematic axis with radius, the gas-derived solution can be made consistent with the stellar kinematic models. A strong kinematic twist in the gas of NGC 2976's center suggests caution, and we advance the mass model based on the stellar kinematics as more reliable. The analysis of this first galaxy shows promising evidence that dark matter halos in late-type dwarfs may in fact be more consistent with cuspy dark matter distributions than earlier work has claimed.

We have developed a spherically symmetric dynamical model of a dwarf spheroidal galaxy using the Schwarzschild method. This type of modelling yields constraints both on the total mass distribution (e.g. enclosed mass and scale radius) as well as on the orbital structure of the system modelled (e.g. velocity anisotropy). Therefore not only can we derive the dark matter content of these systems, but also explore possible formation scenarios. Here we present preliminary results for the Sculptor dSph. We find that the mass of Sculptor within 1kpc is 8.5\times10^(7\pm0.05) M\odot, its anisotropy profile is tangentially biased and slightly more isotropic near the center. For an NFW profile, the preferred concentration (~15) is compatible with cosmological models. Very cuspy density profiles (steeper than NFW) are strongly disfavoured for Sculptor.

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

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

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.

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

VENGA is a large-scale extragalactic IFU survey, which maps the bulges, bars and large parts of the outer disks of 32 nearby normal spiral galaxies. The targets are chosen to span a wide range in Hubble types, star formation activities, morphologies, and inclinations, at the same time of having vast available multi-wavelength coverage from the far-UV to the mid-IR, and available CO and 21cm mapping. The VENGA dataset will provide 2D maps of the SFR, stellar and gas kinematics, chemical abundances, ISM density and ionization states, dust extinction and stellar populations for these 32 galaxies. The uniqueness of the VIRUS-P large field of view permits these large-scale mappings to be performed. VENGA will allow us to correlate all these important quantities throughout the different environments present in galactic disks, allowing the conduction of a large number of studies in star formation, structure assembly, galactic feedback and ISM in galaxies.

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

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.

Most of the super massive black hole mass estimates based on stellar kinematics use the assumption that galaxies are axisymmetric oblate spheroids or spherical. Here we use fully general triaxial orbit-based models to explore the effect of relaxing the axisymmetric assumption on the previously studied galaxies M32 and NGC 3379. We find that M32 can only be modeled accurately using an axisymmetric shape viewed nearly edge-on and our black hole mass estimate is identical to previous studies. When the observed 5 degrees kinematical twist is included in our model of NGC 3379, the best shape is mildly triaxial and we find that our best-fitting black hole mass estimate doubles with respect to the axisymmetric model. This particular black hole mass estimate is still within the errors of that of the axisymmetric model and consistent with the M-sigma relationship. However, this effect may have a pronounced impact on black hole demography, since roughly a third of the most massive galaxies are strongly triaxial.

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

We investigate how well the intrinsic shape of early-type galaxies can be recovered when both photometric and two-dimensional stellar kinematic observations are available. We simulate these observations with galaxy models that are representative of observed oblate fast-rotator to triaxial slow-rotator early-type galaxies. By fitting realistic triaxial dynamical models to these simulated observations, we recover the intrinsic shape (and mass-to-light ratio), without making additional (ad-hoc) assumptions on the orientation. For (near) axisymmetric galaxies the dynamical modelling can strongly exclude triaxiality, but the regular kinematics do not further tighten the constraint on the intrinsic flattening significantly, so that the inclination is nearly unconstrained above the photometric lower limit even with two-dimensional stellar kinematics. Triaxial galaxies can have additional complexity in both the observed photometry and kinematics, such as twists and (central) kinematically decoupled components, which allows the intrinsic shape to be accurately recovered. For galaxies that are very round or show no significant rotation, recovery of the shape is degenerate, unless additional constraints such as from a thin disk are available.

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

[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 construct axisymmetric and triaxial galaxy models with a phase-space distribution function that depends on linear combinations of the three exact integrals of motion for a separable potential. These Abel models, first introduced by Dejonghe & Laurent and subsequently extended by Mathieu & Dejonghe, are the axisymmetric and triaxial generalisations of the well-known spherical Osipkov-Merritt models. We show that the density and higher order velocity moments, as well as the line-of-sight velocity distribution (LOSVD) of these models can be calculated efficiently and that they capture much of the rich internal dynamics of early-type galaxies. We build a triaxial and oblate axisymmetric galaxy model with projected kinematics that mimic the two-dimensional kinematic observations that are obtained with integral-field spectrographs such as SAURON. We fit the simulated observations with axisymmetric and triaxial dynamical models constructed with our numerical implementation of Schwarzschild's orbit-superposition method. We find that Schwarzschild's method is able to recover the internal dynamics and three-integral distribution function of realistic models of early-type galaxies.

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.

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

We construct orbit-based axisymmetric dynamical models for the globular cluster M15 which fit groundbased line-of-sight velocities and Hubble Space Telescope line-of-sight velocities and proper motions. This allows us to constrain the variation of the mass-to-light ratio M/L as a function of radius in the cluster, and to measure the distance and inclination of the cluster. We obtain a best-fitting inclination of 60+/-15 degrees, a dynamical distance of 10.3+/-0.4 kpc and an M/L profile with a central peak. The inferred mass in the central 0.05 parsec is 3400 Msun, implying a central density of at least 7.4x10^6 Msun pc^-3. We cannot distinguish the nature of the central mass concentration. It could be an IMBH or it could be large number of compact objects, or it could be a combination. The central 4 arcsec of M15 appears to contain a rapidly spinning core, and we speculate on its origin.

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 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 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 determine the dynamical distance D, inclination i, mass-to-light ratio M/L and the intrinsic orbital structure of the globular cluster omega Cen, by fitting axisymmetric dynamical models to the ground-based proper motions of van Leeuwen et al. and line-of-sight velocities from four independent data-sets. We correct the observed velocities for perspective rotation caused by the space motion of the cluster, and show that the residual solid-body rotation component in the proper motions can be taken out without any modelling other than assuming axisymmetry. This also provides a tight constraint on D tan i. Application of our axisymmetric implementation of Schwarzschild's orbit superposition method to omega Cen reveals no dynamical evidence for a significant radial dependence of M/L. The best-fit dynamical model has a stellar V-band mass-to-light ratio M/L_V = 2.5 +/- 0.1 M_sun/L_sun and an inclination i = 50 +/- 4 degrees, which corresponds to an average intrinsic axial ratio of 0.78 +/- 0.03. The best-fit dynamical distance D = 4.8 +/- 0.3 kpc (distance modulus 13.75 +/- 0.13 mag) is significantly larger than obtained by means of simple spherical or constant-anisotropy axisymmetric dynamical models, and is consistent with the canonical value 5.0 +/- 0.2 kpc obtained by photometric methods. The total mass of the cluster is (2.5 +/- 0.3) x 10^6 M_sun. The best-fit model is close to isotropic inside a radius of about 10 arcmin and becomes increasingly tangentially anisotropic in the outer region, which displays significant mean rotation. This phase-space structure may well be caused by the effects of the tidal field of the Milky Way. The cluster contains a separate disk-like component in the radial range between 1 and 3 arcmin, contributing about 4% to the total mass.

We have used the ARGUS multi-object spectrometer at the CTIO 4m Blanco telescope to obtain 2756 radial velocity measurements for 1966 individual stars in the globular cluster omega Centauri brighter than blue photographic magnitude of about 16.5. Of these, 1589 stars are cluster members. A comparison with two independent radial velocity studies, carried out by Suntzeff & Kraft and by Mayor et al., demonstrates that the median error of our measurements is below 2 km/s for the stars brighter than B-magnitude 15, which constitute the bulk of the sample. The observed velocity dispersion decreases from about 15 km/s in the inner few arcmin to about 6 km/s at a radius of 25 arcmin. The cluster shows significant rotation, with a maximum amplitude of about 6 km/s in the radial zone between 6 and 10 arcmin. In a companion paper by van de Ven et al., we correct these radial velocities for the perspective rotation caused by the space motion of the cluster, and combine them with the internal proper motions of nearly 8000 cluster members measured by van Leeuwen et al., to construct a detailed dynamical model of omega Centauri and to measure its distance.

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.