Clouds are ubiquitous in exoplanet atmospheres and represent a challenge for the model interpretation of their spectra. Complex cloud models are too numerically costly for generating a large number of spectra, while more efficient models may be too strongly simplified. We aim to constrain the atmospheric properties of the directly imaged planet HR 8799e with a free retrieval approach. We use our radiative transfer code petitRADTRANS for generating spectra, which we couple to the PyMultiNest tool. We added the effect of multiple scattering which is important for treating clouds. Two cloud model parameterizations are tested: the first incorporates the mixing and settling of condensates, the second simply parameterizes the functional form of the opacity. In mock retrievals, using an inadequate cloud model may result in atmospheres that are more isothermal and less cloudy than the input. Applying our framework on observations of HR 8799e made with the GPI, SPHERE and GRAVITY, we find a cloudy atmosphere governed by disequilibrium chemistry, confirming previous analyses. We retrieve that ${\rm C/O}=0.60_{-0.08}^{+0.07}$. Other models have not yet produced a well constrained C/O value for this planet. The retrieved C/O values of both cloud models are consistent, while leading to different atmospheric structures: cloudy, or more isothermal and less cloudy. Fitting the observations with the self-consistent Exo-REM model leads to comparable results, while not constraining C/O. With data from the most sensitive instruments, retrieval analyses of directly imaged planets are possible. The inferred C/O ratio of HR 8799e is independent of the cloud model and thus appears to be a robust. This C/O is consistent with stellar, which could indicate that the HR 8799e formed outside the CO$_2$ or CO iceline. As it is the innermost planet of the system, this constraint could apply to all HR 8799 planets.

We combine data from ALMA and MUSE to study the resolved (~300 pc scale) star formation relation (star formation rate vs. molecular gas surface density) in cluster galaxies. Our sample consists of 9 Fornax cluster galaxies, including spirals, ellipticals, and dwarfs, covering a stellar mass range of ~10^8.8 - 10^11 M_Sun. CO(1-0) and extinction corrected Halpha were used as tracers for the molecular gas mass and star formation rate, respectively. We compare our results with Kennicutt (1998) and Bigiel et al. (2008). Furthermore, we create depletion time maps to reveal small-scale variations in individual galaxies. We explore these further in FCC290, using the 'uncertainty principle for star formation' (Kruijssen & Longmore, 2014a) to estimate molecular cloud lifetimes, which we find to be short (<10 Myr) in this galaxy. Galaxy-averaged depletion times are compared with other parameters such as stellar mass and cluster-centric distance. We find that the star formation relation in the Fornax cluster is close to those from Kennicutt (1998) and Bigiel et al. (2008}), but overlaps mostly with the shortest depletion times predicted by Bigiel et al. (2008). This slight decrease in depletion time is mostly driven by dwarf galaxies with disturbed molecular gas reservoirs close to the virial radius. In FCC90, a dwarf galaxy with a molecular gas tail, we find that depletion times are a factor >~10 higher in its tail than in its stellar body.

Large-amplitude Sgr A* near-infrared flares result from energy injection into electrons near the black hole event horizon. Astrometry data show continuous rotation of the emission region during bright flares, and corresponding rotation of the linear polarization angle. One broad class of physical flare models invokes magnetic reconnection. Here we show that such a scenario can arise in a general relativistic magnetohydrodynamic simulation of a magnetically arrested disc. Saturation of magnetic flux triggers eruption events, where magnetically dominated plasma is expelled from near the horizon and forms a rotating, spiral structure. Dissipation occurs via reconnection at the interface of the magnetically dominated plasma and surrounding fluid. This dissipation is associated with large increases in near-infrared emission in models of Sgr A*, with durations and amplitudes consistent with the observed flares. Such events occur at roughly the timescale to re-accumulate the magnetic flux from the inner accretion disc, 10h for Sgr A*. We study near-infrared observables from one sample event to show that the emission morphology tracks the boundary of the magnetically dominated region. As the region rotates around the black hole, the near-infrared centroid and linear polarization angle both undergo continuous rotation, similar to the behavior seen in Sgr A* flares.

Globular clusters (GCs) are found ubiquitously in massive galaxies and due to their old ages, they are regarded as fossil records of galaxy evolution. Spectroscopic studies of GC systems are often limited to the outskirts of galaxies, where GCs stand out against the galaxy background and serve as bright tracers of galaxy assembly. In this work, we use the capabilities of the Multi Unit Explorer Spectrograph (MUSE) to extract a spectroscopic sample of 722 GCs in the inner regions ($\lesssim 3 R_\text{eff}$) of 32 galaxies in the Fornax cluster. These galaxies were observed as part of the Fornax 3D project, a MUSE survey that targets early and late-type galaxies within the virial radius of Fornax. After accounting for the galaxy background in the GC spectra, we extracted line-of-sight velocities and determined metallicities of a sub-sample of 238 GCs. We found signatures of rotation within GC systems, and comparing the GC kinematics and that of the stellar body shows that the GCs trace the spheroid of the galaxies. While the red GCs prove to closely follow the metallicity profile of the host galaxy, the blue GCs show a large spread of metallicities but they are generally more metal-poor than the host.

Globular cluster (GC) systems of massive galaxies often show a bimodal colour distribution. This has been interpreted as a metallicity bimodality, created by a two-stage galaxy formation where the red, metal-rich GCs were formed in the parent halo and the blue metal-poor GCs were accreted. This interpretation, however, crucially depends on the assumption that GCs are exclusively old stellar systems with a linear colour-metallicity relation (CZR). The shape of the CZR and range of GC ages are currently under debate, because their study requires high quality spectra to derive reliable stellar population properties. We determined metallicities with full spectral fitting from a sample of 187 GCs with high spectral signal-to-noise ratio in 23 galaxies of the Fornax cluster that were observed as part of the Fornax 3D project. The derived CZR from this sample is non-linear and can be described by a piecewise linear function with a break point at ($g - z$) $\sim$ 1.1 mag. The less massive galaxies in our sample ($M_\ast < 10^{10} M_\odot$) appear to have slightly younger GCs, but the shape of the CZR is insensitive to the GC ages. Although the least massive galaxies lack red, metal-rich GCs, a non-linear CZR is found irrespective of the galaxy mass, even in the most massive galaxies ($M_\ast \geq 10^{11} M_\odot$). Our CZR predicts narrow unimodal GC metallicity distributions for low mass and broad unimodal distributions for very massive galaxies, dominated by a metal-poor and metal-rich peak, respectively, and bimodal distributions for galaxies with intermediate masses (10$^{10}$ $\leq$ $M_\ast < 10^{11} M_\odot$) as a consequence of the relative fraction of red and blue GCs. The diverse metallicity distributions challenge the simple differentiation of GC populations solely based on their colour.

Extragalactic planetary nebulae (PNe) are detectable through relatively strong nebulous [OIII] emission and act as direct probes into the local stellar population. Because they have an apparently universal invariant magnitude cut-off, PNe are also considered to be a remarkable standard candle for distance estimation. Through detecting PNe within the galaxies, we aim to connect the relative abundances of PNe to the properties of their host galaxy stellar population. By removing the stellar background components from FCC 167 and FCC 219, we aim to produce PN luminosity functions (PNLF) of these galaxies, and thereby also estimate the distance modulus to these two systems. Finally, we test the reliability and robustness of our novel detection and analysis method. It detects unresolved point sources by their [OIII] 5007{\AA} emission within regions that have previously been unexplored. We model the [OIII] emissions in the spatial and spectral dimensions together, as afforded to us by the Multi Unit Spectroscopic Explorer (MUSE), and we draw on data gathered as part of the Fornax3D survey. For each source, we inspect the properties of the nebular emission lines to remove other sources that might hinder the safe construction of the PNLF, such as supernova remnants and HII regions. As a further step, we characterise any potential limitations and draw conclusions about the reliability of our modelling approach through a set of simulations. By applying this novel detection and modelling approach to integral field unit observations, we report for the distance estimates and luminosity-specific PNe frequency values for the two galaxies. Furthermore, we include an overview of source contamination, galaxy differences, and possible effects on the PNe populations in the dense stellar environments.

The Galactic Center black hole Sagittarius A* is a variable NIR source that exhibits bright flux excursions called flares. The low-flux density turnover of the flux distribution is below the sensitivity of current single-aperture telescopes. We use the unprecedented resolution of the GRAVITY instrument at the VLTI. Our light curves are unconfused, overcoming the confusion limit of previous photometric studies. We analyze the light curves using standard statistical methods and obtain the flux distribution. We find that the flux distribution of SgrA* turns over at a median flux density of (1.1\pm0.3)mJy. We measure the percentiles of the flux distribution and use them to constrain the NIR K-band SED. Furthermore, we find that the flux distribution is intrinsically right-skewed to higher flux density in log space. Flux densities below 0.1mJy are hardly ever observed. In consequence, a single powerlaw or lognormal distribution does not suffice to describe the observed flux distribution in its entirety. However, if one takes into account a power law component at high flux densities, a lognormal distribution can describe the lower end of the observed flux distribution. We confirm the RMS-flux relation for Sgr~A* and find it to be linear for all flux densities in our observation. We conclude that Sgr~A* has two states: the bulk of the emission is generated in a lognormal process with a well-defined median flux density and this quiescent emission is supplemented by sporadic flares that create the observed power law extension of the flux distribution.

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

The star S2 orbiting the compact radio source Sgr A* is a precision probe of the gravitational field around the closest massive black hole (candidate). Over the last 2.7 decades we have monitored the star's radial velocity and motion on the sky, mainly with the SINFONI and NACO adaptive optics (AO) instruments on the ESO VLT, and since 2017, with the four-telescope interferometric beam combiner instrument GRAVITY. In this paper we report the first detection of the General Relativity (GR) Schwarzschild Precession (SP) in S2's orbit. Owing to its highly elliptical orbit (e = 0.88), S2's SP is mainly a kink between the pre-and post-pericentre directions of motion ~ +- 1 year around pericentre passage, relative to the corresponding Kepler orbit. The superb 2017-2019 astrometry of GRAVITY defines the pericentre passage and outgoing direction. The incoming direction is anchored by 118 NACO-AO measurements of S2's position in the infrared reference frame, with an additional 75 direct measurements of the S2-Sgr A* separation during bright states ('flares') of Sgr A*. Our 14-parameter model fits for the distance, central mass, the position and motion of the reference frame of the AO astrometry relative to the mass, the six parameters of the orbit, as well as a dimensionless parameter f_SP for the SP (f_SP = 0 for Newton and 1 for GR). From data up to the end of 2019 we robustly detect the SP of S2, del phi = 12' per orbital period. From posterior fitting and MCMC Bayesian analysis with different weighting schemes and bootstrapping we find f_SP = 1.10 +- 0.19. The S2 data are fully consistent with GR. Any extended mass inside S2's orbit cannot exceed ~ 0.1% of the central mass. Any compact third mass inside the central arcsecond must be less than about 1000 M_sun.

The Galactic center black hole candidate Sgr A* is the best target for studies of low-luminosity accretion physics, including with near-infrared and submillimeter wavelength long baseline interferometry experiments. Here we compare images and spectra generated from a parameter survey of general relativistic MHD simulations to a set of radio to near-infrared observations of Sgr A*. Our models span the limits of weak and strong magnetization and use a range of sub-grid prescriptions for electron heating. We find two classes of scenarios can explain the broad shape of the submillimeter spectral peak and the highly variable near-infrared flaring emission. Weakly magnetized "disk-jet" models where most of the emission is produced near the jet wall, consistent with past work, as well as strongly magnetized (magnetically arrested disk) models where hot electrons are present everywhere. Disk-jet models are strongly depolarized at submillimeter wavelengths as a result of strong Faraday rotation, inconsistent with observations of Sgr A*. We instead favor the strongly magnetized models, which provide a good description of the median and highly variable linear polarization signal. The same models can also explain the observed mean Faraday rotation measure and potentially the polarization signals seen recently in Sgr A* near-infrared flares.

We report on the determination of electron densities, and their impact on the outflow masses and rates, measured in the central few hundred parsecs of 11 local luminous active galaxies. We show that the peak of the integrated line emission in the AGN is significantly offset from the systemic velocity as traced by the stellar absorption features, indicating that the profiles are dominated by outflow. In contrast, matched inactive galaxies are characterised by a systemic peak and weaker outflow wing. We present three independent estimates of the electron density in these AGN, discussing the merits of the different methods. The electron density derived from the [SII] doublet is significantly lower than than that found with a method developed in the last decade using auroral and transauroral lines, as well as a recently introduced method based on the ionization parameter. The reason is that, for gas photoionized by an AGN, much of the [SII] emission arises in an extended partially ionized zone where the implicit assumption that the electron density traces the hydrogen density is invalid. We propose ways to deal with this situation and we derive the associated outflow rates for ionized gas, which are in the range 0.001--0.5 M$_{\odot}$yr$^{-1}$ for our AGN sample. We compare these outflow rates to the relation between $\dot{M}_{out}$ and $L_{AGN}$ in the literature, and argue that it may need to be modified and rescaled towards lower mass outflow rates.

We present stellar metallicity measurements of more than 600 late-type stars in the central 10 pc of the Galactic centre. Together with our previously published KMOS data, this data set allows us to investigate, for the first time, spatial variations of the nuclear star cluster's metallicity distribution. Using the integral-field spectrograph KMOS (VLT) we observed almost half of the area enclosed by the nuclear star cluster's effective radius. We extract spectra at medium spectral resolution, and apply full spectral fitting utilising the PHOENIX library of synthetic stellar spectra. The stellar metallicities range from [M/H]=-1.25 dex to [M/H]> +0.3 dex, with most of the stars having super-solar metallicity. We are able to measure an anisotropy of the stellar metallicity distribution. In the Galactic North, the portion of sub-solar metallicity stars with [M/H]<0.0 dex is more than twice as high as in the Galactic South. One possible explanation for different fractions of sub-solar metallicity stars in different parts of the cluster is a recent merger event. We propose to test this hypothesis with high-resolution spectroscopy, and by combining the metallicity information with kinematic data.

The inner regions of the discs of high-mass young stellar objects (HMYSOs) are still poorly known due to the small angular scales and the high visual extinction involved. We deploy near-infrared (NIR) spectro-interferometry to probe the inner gaseous disc in HMYSOs and investigate the origin and physical characteristics of the CO bandhead emission (2.3-2.4 $\mu$m). We present the first GRAVITY/VLTI observations at high spectral (R=4000) and spatial (mas) resolution of the CO overtone transitions in NGC 2024 IRS2. The continuum emission is resolved in all baselines and is slightly asymmetric, displaying small closure phases ($\leq$8$^{\circ}$). Our best ellipsoid model provides a disc inclination of 34$^{\circ}$$\pm$1$^{\circ}$, a disc major axis position angle of 166$^{\circ}$$\pm$1$^{\circ}$, and a disc diameter of 3.99$\pm$0.09 mas (or 1.69$\pm$0.04 au, at a distance of 423 pc). The small closure phase signals in the continuum are modelled with a skewed rim, originating from a pure inclination effect. For the first time, our observations spatially and spectrally resolve the first four CO bandheads. Changes in visibility, as well as differential and closure phases across the bandheads are detected. Both the size and geometry of the CO-emitting region are determined by fitting a bidimensional Gaussian to the continuum-compensated CO bandhead visibilities. The CO-emitting region has a diameter of 2.74$\pm^{0.08}_{0.07}$ mas (1.16$\pm$0.03 au), and is located in the inner gaseous disc, well within the dusty rim, with inclination and $PA$ matching the dusty disc geometry, which indicates that both dusty and gaseous discs are coplanar. Physical and dynamical gas conditions are inferred by modelling the CO spectrum. Finally, we derive a direct measurement of the stellar mass of $M_*\sim$14.7$^{+2}_{-3.6}$ M$_{\odot}$ by combining our interferometric and spectral modelling results.

Infrared observations of Sgr A* probe the region close to the event horizon of the black hole at the Galactic center. These observations can constrain the properties of low-luminosity accretion as well as that of the black hole itself. The GRAVITY instrument at the ESO VLTI has recently detected continuous circular relativistic motion during infrared flares which has been interpreted as orbital motion near the event horizon. Here we analyze the astrometric data from these flares, taking into account the effects of out-of-plane motion and orbital shear of material near the event horizon of the black hole. We have developed a new code to predict astrometric motion and flux variability from compact emission regions following particle orbits. Our code combines semi-analytic calculations of timelike geodesics that allow for out-of-plane or elliptical motions with ray tracing of photon trajectories to compute time-dependent images and light curves. We apply our code to the three flares observed with GRAVITY in 2018. We show that all flares are consistent with a hotspot orbiting at R$\sim$9 gravitational radii with an inclination of $i\sim140^\circ$. The emitting region must be compact and less than $\sim5$ gravitational radii in diameter. We place a further limit on the out-of-plane motion during the flare.

We study the evidence for a diversity of formation processes in early-type galaxies by presenting the first complete volume-limited sample of slow rotators with both integral-field kinematics from the ATLAS3D Project and high spatial resolution photometry from the Hubble Space Telescope. Analysing the nuclear surface brightness profiles of 12 newly imaged slow rotators, we classify their light profiles as core-less, and place an upper limit to the core size of ~10 pc. Considering the full magnitude and volume-limited ATLAS3D sample, we correlate the presence or lack of cores with stellar kinematics, including the proxy for the stellar angular momentum and the velocity dispersion within one half-light radius, stellar mass, stellar age, $\alpha$-element abundance, and age and metallicity gradients. More than half of the slow rotators have core-less light profiles, and they are all less massive than $10^{11}$ Msun. Core-less slow rotators show evidence for counter-rotating flattened structures, have steeper metallicity gradients, and a larger dispersion of gradient values than core slow rotators. Our results suggest that core and core-less slow rotators have different assembly processes, where the former are the relics of massive dissipation-less merging in the presence of central supermassive black holes. Formation processes of core-less slow rotators are consistent with accretion of counter-rotating gas or gas-rich mergers of special orbital configurations, which lower the final net angular momentum of stars, but support star formation. We also highlight core fast rotators as galaxies that share properties of core slow rotators and core-less slow rotators. Formation processes similar to those for core-less slow rotators can be invoked to explain the assembly of core fast rotators, with the distinction that these processes form or preserve cores.[Abridged]

Our objective is to estimate the C/O ratio in the atmosphere of beta Pictoris b and obtain an estimate of the dynamical mass of the planet, as well as to refine its orbital parameters using high-precision astrometry. We used the GRAVITY instrument with the four 8.2 m telescopes of the Very Large Telescope Interferometer to obtain K-band spectro-interferometric data on $\beta$ Pic b. We extracted a medium resolution (R=500) K-band spectrum of the planet and a high-precision astrometric position. We estimated the planetary C/O ratio using two different approaches (forward modeling and free retrieval) from two different codes (ExoREM and petitRADTRANS, respectively). Finally, we used a simplified model of two formation scenarios (gravitational collapse and core-accretion) to determine which can best explain the measured C/O ratio. Our new astrometry disfavors a circular orbit for $\beta$ Pic b ($e=0.15^{+0.05}_{-0.04}$). Combined with previous results and with Hipparcos/GAIA measurements, this astrometry points to a planet mass of $M = 12.7\pm{}2.2\,M_\mathrm{Jup}$. This value is compatible with the mass derived with the free-retrieval code petitRADTRANS using spectral data only. The forward modeling and free-retrieval approches yield very similar results regarding the atmosphere of beta Pic b. In particular, the C/O ratios derived with the two codes are identical ($0.43\pm{}0.05$ vs $0.43^{+0.04}_{-0.03}$). We argue that if the stellar C/O in $\beta$ Pic is Solar, then this combination of a very high mass and a low C/O ratio for the planet suggests a formation through core-accretion, with strong planetesimal enrichment.

We present near-infrared interferometric data on the Seyfert 2 galaxy NGC 1068, obtained with the GRAVITY instrument on the European Southern Observatory Very Large Telescope Interferometer. The extensive baseline coverage from 5 to 60 M\lambda allowed us to reconstruct a continuum image of the nucleus with an unrivaled 0.2 pc resolution in the K-band. We find a thin ring-like structure of emission with a radius r = 0.24+/-0.03 pc, inclination i = 70+/-5 deg, position angle PA = -50+/-4 deg, and h/r < 0.14, which we associate with the dust sublimation region. The observed morphology is inconsistent with the expected signatures of a geometrically and optically thick torus. Instead, the infrared emission shows a striking resemblance to the 22 GHz maser disc, which suggests they share a common region of origin. The near-infrared spectral energy distribution indicates a bolometric luminosity of (0.4-4.7) x 10^45 erg/s, behind a large A_K ~ 5.5 (A_V ~ 90) screen of extinction that also appears to contribute significantly to obscuring the broad line region.

The formation and the evolution of protoplanetary disks are important stages in the lifetime of stars. The processes of disk evolution and planet formation are intrinsically linked. We spatially resolve with GRAVITY/VLTI in the K-band the sub au-scale region of 27 stars to gain statistical understanding of their properties. We look for correlations with stellar parameters, such as luminosity, mass, temperature and age. Our sample also cover a range of various properties in terms of reprocessed flux, flared or flat morphology, and gaps. We developed semi-physical geometrical models to fit our interferometric data. Our best models correspond to smooth and wide rings, implying that wedge-shaped rims at the dust sublimation edge are favored, as found in the H-band. The closure phases are generally non-null with a median value of ~10 deg, indicating spatial asymmetries of the intensity distributions. Multi-size grain populations could explain the closure phase ranges below 20-25 deg but other scenarios should be invoked to explain the largest ones. Our measurements extend the Radius-Luminosity relation to ~1e4 Lsun and confirm the significant spread around the mean relation observed in the H-band. Gapped sources exhibit a large N-to-K band size ratio and large values of this ratio are only observed for the members of our sample that would be older than 1 Ma, less massive, and with lower luminosity. In the 2 Ms mass range, we observe a correlation in the increase of the relative age with the transition from group II to group I, and an increase of the N-to-K size ratio. However, the size of the current sample does not yet permit us to invoke a clear universal evolution mechanism across the HAeBe mass range. The measured locations of the K-band emission suggest that these disks might be structured by forming young planets, rather than by depletion due to EUV, FUV, and X-ray photo-evaporation.

We use VLTI/GRAVITY near-infrared interferometry measurements of 8 bright, Type 1 AGN to study the size and structure of hot dust heated by the central engine. We partially resolve each source, and report Gaussian FWHM sizes in the range 0.3-0.8 milliarcseconds. In all but one object, we find no evidence for significant elongation or asymmetry (closure phases < 1 deg). The effective physical radius increases with bolometric luminosity as found from past reverberation and interferometry measurements. The measured sizes for Seyfert galaxies are systematically larger than for the two quasars in our sample when measured relative to the previously reported R ~ L^1/2 relationship explained by emission at the sublimation radius. This could be evidence of evolving near-infrared emission region structure as a function of central luminosity.

The velocity distribution of stars is a sensitive probe of the gravitational potential of the Galaxy, and hence of its dark matter distribution. In particular, the shape of the dark halo (e.g. spherical, oblate, or prolate) determines velocity correlations, and different halo geometries are expected to result in measurable differences. Here we explore and interpret the correlations in the $(v_R, v_z)$-velocity distribution as a function of position in the Milky Way. We selected a high-quality sample of stars from the Gaia DR2 catalogue and characterised the orientation of the velocity distribution or tilt angle over a radial distance range of $[4-13]~$kpc and up to $3.5~$kpc away from the Galactic plane while taking into account the effects of the measurement errors. We find that the tilt angles change from spherical alignment in the inner Galaxy ($R\sim4~$kpc) towards more cylindrical alignments in the outer Galaxy ($R\sim11~$kpc) when using distances that take a global zero-point offset in the parallax of $-29~\mu$as. However, if the amplitude of this offset is underestimated, then the inferred tilt angles in the outer Galaxy only appear shallower and are intrinsically more consistent with spherical alignment for an offset as large as $-54~\mu$as. We further find that the tilt angles do not seem to strongly vary with Galactic azimuth and that different stellar populations depict similar tilt angles. Therefore we introduce a simple analytic function that describes the trends found over the full radial range. Since the systematic parallax errors in Gaia DR2 depend on celestial position, magnitude, and colour in complex ways, it is not possible to fully correct for them. Therefore it will be particularly important for dynamical modelling of the Milky Way to thoroughly characterise the systematics in astrometry in future Gaia data releases.

The 31 brightest galaxies (m_B < 15 mag) inside the virial radius of the Fornax cluster were observed from the centres to the outskirts with the Multi Unit Spectroscopic Explorer on the Very Large Telescope. These observations provide detailed high-resolution maps of the line-of-sight kinematics and line strengths of the stars and ionised gas reaching 2-3 Re for 21 early-type galaxies and 1-2 Re for 10 late-type galaxies. The majority of the galaxies are regular rotators, with eight hosting a kinematically distinct core. Only two galaxies are slow rotators. The mean age, total metallicity, and [Mg/Fe] abundance ratio in the bright central region inside 0.5 Re and in the galaxy outskirts are presented. Extended emission-line gas is detected in 13 galaxies, most of them are late-type objects with wide-spread star formation. The measured structural properties are analysed in relation to the galaxies' position in the projected phase space of the cluster. This shows that the Fornax cluster appears to consist of three main groups of galaxies inside the virial radius: the old core; a clump of galaxies, which is aligned with the local large-scale structure and was accreted soon after the formation of the core; and a group of galaxies that fell in more recently.

We present a study of the three dimensional structure, kinematics, and age distribution of the Orion OB association, based on the second data release of the Gaia satellite (Gaia DR2). Our goal is to obtain a complete picture of the star formation history of the Orion complex and to relate our findings to theories of sequential and triggered star formation. We select the Orion population with simple photometric criteria, and we construct a three dimensional map in galactic Cartesian coordinates to study the physical arrangement of the stellar clusters in the Orion region. The map shows structures that extend for roughly $150 \, \mathrm{pc}$ along the line of sight, divided in multiple sub-clusters. We separate different groups by using the density based clustering algorithm DBSCAN. We study the kinematic properties of all the groups found by DBSCAN first by inspecting their proper motion distribution, and then by applying a kinematic modelling code based on an iterative maximum likelihood approach, which we use to derive their mean velocity, velocity dispersion and isotropic expansion. By using an isochrone fitting procedure we provide ages and extinction values for all the groups. We confirm the presence of an old population ($\sim 15$ Myr) towards the 25 Ori region, and we find that groups with ages of $12-15 \, \mathrm{Myr}$ are present also towards the Belt region. We notice the presence of a population of $\sim 10$ Myr also in front of the Orion A molecular cloud. Our findings suggest that star formation in Orion does not follow a simple sequential scenario, but instead consists of multiple events, which caused kinematic and physical sub-structure. To fully explain the detailed sequence of events, specific simulations and further radial velocity data are needed.

We present a 0.16% precise and 0.27% accurate determination of R0, the distance to the Galactic Center. Our measurement uses the star S2 on its 16-year orbit around the massive black hole Sgr A* that we followed astrometrically and spectroscopically for 27 years. Since 2017, we added near-infrared interferometry with the VLTI beam combiner GRAVITY, yielding a direct measurement of the separation vector between S2 and Sgr A* with an accuracy as good as 20 micro-arcsec in the best cases. S2 passed the pericenter of its highly eccentric orbit in May 2018, and we followed the passage with dense sampling throughout the year. Together with our spectroscopy, in the best cases with an error of 7 km/s, this yields a geometric distance estimate: R0 = 8178 +- 13(stat.) +- 22(sys.) pc. This work updates our previous publication in which we reported the first detection of the gravitational redshift in the S2 data. The redshift term is now detected with a significance level of 20 sigma with f_redshift = 1.04 +- 0.05.

We used deep MUSE observations to perform a stellar-kinematic and population analysis of FCC 153 and FCC 177, two edge-on S0 galaxies in the Fornax cluster. The geometrical definition of the different structural components of these two galaxies allows us to describe the nature of their thick disks. These are both old, relatively metal poor and [Mg/Fe]-enhanced, and their star formation history (SFH) reveals a minor younger component whose chemical properties suggest its later accretion. Moreover, the outer regions of these geometrically defined thick disks show higher values of metallicity and lower values of [Mg/Fe]. These stars probably formed in the thin-disk region and they were dynamically heated to form the flares present in these two galaxies. We propose different formation scenarios for the three populations of these thick disks: in-situ formation, accretion and disk heating. A clear distinction in age is found between the metal poor and [Mg/Fe]-enhanced thick disks (old, $\sim 12-13$ Gyr), and the metal rich and less [Mg/Fe]-enhanced thin disks (young, $\sim 4-5$ Gyr). These two galaxies show signs of relatively recent star formation in their thin disks and nuclear regions. While the thin disks show more continuous SFHs, the nuclei display a rather bursty SFH. These two galaxies are located outside of the densest region of the Fornax cluster where FCC 170 resides. This other edge-on S0 galaxy was studied by \citet{Pinna2019}. We compare and discuss our results with this previous study. The differences between these three galaxies, at different distances from the cluster center, suggest that the environment can have a strong effect on the galaxy evolutionary path.

To date, infrared interferometry at best achieved contrast ratios of a few times $10^{-4}$ on bright targets. GRAVITY, with its dual-field mode, is now capable of high contrast observations, enabling the direct observation of exoplanets. We demonstrate the technique on HR8799, a young planetary system composed of four known giant exoplanets. We used the GRAVITY fringe tracker to lock the fringes on the central star, and integrated off-axis on the HR8799e planet situated at 390 mas from the star. Data reduction included post-processing to remove the flux leaking from the central star and to extract the coherent flux of the planet. The inferred K band spectrum of the planet has a spectral resolution of 500. We also derive the astrometric position of the planet relative to the star with a precision on the order of 100$\,\mu$as. The GRAVITY astrometric measurement disfavors perfectly coplanar stable orbital solutions. A small adjustment of a few degrees to the orbital inclination of HR 8799 e can resolve the tension, implying that the orbits are close to, but not strictly coplanar. The spectrum, with a signal-to-noise ratio of $\approx 5$ per spectral channel, is compatible with a late-type L brown dwarf. Using Exo-REM synthetic spectra, we derive a temperature of $1150\pm50$\,K and a surface gravity of $10^{4.3\pm0.3}\,$cm/s$^{2}$. This corresponds to a radius of $1.17^{+0.13}_{-0.11}\,R_{\rm Jup}$ and a mass of $10^{+7}_{-4}\,M_{\rm Jup}$, which is an independent confirmation of mass estimates from evolutionary models. Our results demonstrate the power of interferometry for the direct detection and spectroscopic study of exoplanets at close angular separations from their stars.

The stellar initial mass function (IMF) regulates the baryonic cycle within galaxies, and is a key ingredient to translate observations into physical quantities. Although for decades it was assumed to be universal, there is now growing observational evidence showing that the center of massive early-type galaxies host an enhanced population of low-mass stars compared to the expectations from the Milky Way. Moreover, these variations in the IMF have been found to be related to the radial metallicity variations in massive galaxies. We present here a two-dimensional stellar population analysis of the massive lenticular galaxy FCC 167 (NGC 1380) as part of the Fornax3D project. Using a newly developed stellar population fitting scheme, we derive a full two-dimensional IMF map of an early-type galaxy. This two-dimensional analysis allows us go further than a radial analysis, showing how the metallicity changes along a disc-like structure while the IMF follows a distinct, less disky distribution. Thus, our findings indicate that metallicity cannot be the sole driver of the observed radial IMF variations. In addition, a comparison with the orbital decomposition shows suggestive evidence of a coupling between stellar population properties and the internal dynamical structure of FCC 167, where metallicity and IMF maps seem to track the distribution of cold and warm orbits, respectively.

During its orbit around the four million solar mass black hole Sagittarius A* the star S2 experiences significant changes in gravitational potential. We use this change of potential to test one part of the Einstein equivalence principle: the local position invariance (LPI). We study the dependency of different atomic transitions on the gravitational potential to give an upper limit on violations of the LPI. This is done by separately measuring the redshift from hydrogen and helium absorption lines in the stellar spectrum during its closest approach to the black hole. For this measurement we use radial velocity data from 2015 to 2018 and combine it with the gravitational potential at the position of S2, which is calculated from the precisely known orbit of S2 around the black hole. This results in a limit on a violation of the LPI of $|\beta_{He}-\beta_{H}| = (2.4 \pm 5.1) \cdot 10^{-2}$. The variation in potential that we probe with this measurement is six magnitudes larger than possible for measurements on Earth, and a factor ten larger than in experiments using white dwarfs. We are therefore testing the LPI in a regime where it has not been tested before.

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

In a dynamically relaxed cluster around a massive black hole a dense stellar cusp of old stars is expected to form. Previous observations showed a relative paucity of red giant stars within the central 0.5 pc in the Galactic Center. By co-adding spectroscopic observations taken over a decade, we identify new late-type stars, including the first five warm giants (G2-G8III), within the central 1 arcsec 2 (0.04 {\times} 0.04 pc^2) of the Galaxy. Our findings increase the number of late-type stars to 21, of which we present deep spectra for 16. The updated star count, based on individual spectral classification, is used to reconstruct the surface density profile of giant stars. Our study, for the first time, finds a cusp in the surface number density of the spectroscopically identified old (>3 Gyr) giants population (m K<17) within 0.02-0.4 pc described by a single power law with an exponent {\Gamma}= 0.34 {\pm} 0.04.

We present and discuss the stellar kinematics and populations of the S0 galaxy FCC 170 (NGC 1381) in the Fornax cluster, using deep MUSE data from the Fornax 3D survey. We show the maps of the first four moments of the stellar line-of-sight velocity distribution and of the mass-weighted mean stellar age, metallicity and [Mg/Fe] abundance ratio. The high-quality MUSE stellar kinematic measurements unveil the structure of this massive galaxy: a nuclear disk, a bar seen as a boxy bulge with a clear higher-velocity-dispersion X shape, a fast-rotating and flaring thin disk and a slower rotating thick disk. Whereas their overall old age makes it difficult to discuss differences in the formation epoch between these components, we find a clear-cut distinction between metal-rich and less [Mg/Fe]-enhanced populations in the thin-disk, boxy-bulge and nuclear disk, and more metal-poor and [Mg/Fe]-enhanced stars in the thick disk. Located in the densest region of the Fornax cluster, where signs of tidal stripping have been recently found, the evolution of FCC 170 might have been seriously affected by its environment. We discuss the possibility of its "pre-processing" in a subgroup before falling into the present-day cluster, which would have shaped this galaxy a long time ago. The thick disk displays a composite star formation history, as a significant fraction of younger stars co-exist with the main older thick-disk population. The former sub-population is characterized by even lower-metallicity and higher-[Mg/Fe] values, suggesting that these stars formed later and faster in a less chemically evolved satellite, which was subsequently accreted. Finally, we discuss evidence that metal-rich and less [Mg/Fe]-enhanced stars were brought in the outer parts of the thick disk by the flaring of the thin disk.

Galaxies continuously reprocess their interstellar material. One can therefore expect changing dust grain properties in galaxies which have followed different evolutionary pathways. Determining the intrinsic dust grain mix of a galaxy helps in reconstructing its evolutionary history. Early-type galaxies occasionally display regular dust lanes in their central regions. Due to the relatively simple geometry and composition of their stellar bodies, these galaxies are ideal to disentangle dust mix variations from geometric effects. We therefore model the various components of such a galaxy (FCC 167). We reconstruct its recent history, and investigate the possible fate of the dust lane. MUSE and ALMA observations reveal a nested ISM structure. An ionised-gas disk pervades the central regions of FCC 167, including those occupied by the main dust lane. Inward of the dust lane, we also find a disk/ring of cold molecular gas where stars are forming and HII regions contribute to the ionised-gas emission. Further in, the gas ionisation points towards an active galactic nucleus and the fuelling of a central supermassive black hole from its surrounding ionised and molecular reservoir. Observational constraints and radiative transfer models suggest the dust and gas are distributed in a ring-like geometry and the dust mix lacks small grains. The derived dust destruction timescales from sputtering in hot gas are short and we conclude that the dust must be strongly self-shielding and clumpy, or will quickly be eroded and disappear. Our findings show how detailed analysis of individual systems can complement statistical studies of dust-lane ETGs.

The Galactic Center black hole Sgr A* is the archetypical example of an underfed massive black hole. The extremely low accretion rate can be understood in radiatively inefficient accretion flow models. Testing those models has proven to be difficult due to the lack of suitable probes. Radio and submm polarization measurements constrain the flow very close to the event horizon. X-ray observations resolving the Bondi radius yield an estimate roughly four orders of magnitude further out. Here, we present a new, indirect measurement of the accretion flow density at intermediate radii. We use the dynamics of the gas cloud G2 to probe the ambient density. We detect the presence of a drag force slowing down G2 with a statistical significance of approx 9 {\sigma}. This probes the accretion flow density at around 1000 Schwarzschild radii and yields a number density of approx. 4 x 10^3 cm^-3. Self-similar accretion models where the density follows a power law radial profile between the inner zone and the Bondi radius have predicted similar values.

We study the three dimensional arrangement of young stars in the solar neighbourhood using the second release of the Gaia mission (Gaia DR2) and we provide a new, original view of the spatial configuration of the star forming regions within 500 pc from the Sun. By smoothing the star distribution through a gaussian filter, we construct three dimensional density maps for early-type stars (upper-main sequence, UMS) and pre-main sequence (PMS) sources. The PMS and the UMS samples are selected through a combination of photometric and astrometric criteria. A side product of the analysis is a three dimensional, G-band extinction map, which we use to correct our colour-magnitude diagram for extinction and reddening. Both density maps show three prominent structures, Scorpius-Centaurus, Orion, and Vela. The PMS map shows a plethora of lower mass star forming regions, such as Taurus, Perseus, Cepheus, Cassiopeia, and Lacerta, which are less visible in the UMS map, due to the lack of large numbers of bright, early-type stars. We report the finding of a candidate new open cluster towards $l, b \sim 218.5^{\circ}, -2^{\circ}$, which could be related to the Orion star forming complex. We estimate ages for the PMS sample and we study the distribution of PMS stars as a function of their age. We find that younger stars cluster in dense, compact clumps, and are surrounded by older sources, whose distribution is instead more diffuse. The youngest groups that we find are mainly located in Scorpius-Centaurus, Orion, Vela, and Taurus. Cepheus, Cassiopeia, and Lacerta are instead more evolved and less numerous. Finally, we find that the three dimensional density maps show no evidence for the existence of the ring-like structure which is usually referred to as the Gould Belt.

We report the detection of continuous positional and polarization changes of the compact source SgrA* in high states ('flares') of its variable near- infrared emission with the near-infrared GRAVITY-Very Large Telescope Interferometer (VLTI) beam-combining instrument. In three prominent bright flares, the position centroids exhibit clockwise looped motion on the sky, on scales of typically 150 micro-arcseconds over a few tens of minutes, corresponding to about 30% the speed of light. At the same time, the flares exhibit continuous rotation of the polarization angle, with about the same 45(+/-15)-minute period as that of the centroid motions. Modelling with relativistic ray tracing shows that these findings are all consistent with a near face-on, circular orbit of a compact polarized 'hot spot' of infrared synchrotron emission at approximately six to ten times the gravitational radius of a black hole of 4 million solar masses. This corresponds to the region just outside the innermost, stable, prograde circular orbit (ISCO) of a Schwarzschild-Kerr black hole, or near the retrograde ISCO of a highly spun-up Kerr hole. The polarization signature is consistent with orbital motion in a strong poloidal magnetic field.

The highly elliptical, 16-year-period orbit of the star S2 around the massive black hole candidate Sgr A* is a sensitive probe of the gravitational field in the Galactic centre. Near pericentre at 120 AU, ~1400 Schwarzschild radii, the star has an orbital speed of ~7650 km/s, such that the first-order effects of Special and General Relativity have now become detectable with current capabilities. Over the past 26 years, we have monitored the radial velocity and motion on the sky of S2, mainly with the SINFONI and NACO adaptive optics instruments on the ESO Very Large Telescope, and since 2016 and leading up to the pericentre approach in May 2018, with the four-telescope interferometric beam-combiner instrument GRAVITY. From data up to and including pericentre, we robustly detect the combined gravitational redshift and relativistic transverse Doppler effect for S2 of z ~ 200 km/s / c with different statistical analysis methods. When parameterising the post-Newtonian contribution from these effects by a factor f, with f = 0 and f = 1 corresponding to the Newtonian and general relativistic limits, respectively, we find from posterior fitting with different weighting schemes f = 0.90 +/- 0.09 (stat) +\- 0.15 (sys). The S2 data are inconsistent with pure Newtonian dynamics.

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

The Fornax cluster provides a uniquely compact laboratory to study the detailed history of early-type galaxies and the role played by environment in driving their evolution and their transformation from late-type galaxies. Using the superb capabilities of the Multi Unit Spectroscopic Explorer on the Very Large Telescope, high-quality integral-field spectroscopic data were obtained for the inner regions of all the bright ($m_B\leq15$) galaxies within the virial radius of Fornax. The stellar haloes of early-type galaxies are also covered out to about four effective radii. State-of-the-art stellar dynamical and population modelling allows to aim in particular at better characterising the disc components of fast-rotating early-type galaxies, constraining radial variations in the stellar initial-mass functions and measuring the stellar age, metallicity, and $\alpha$-element abundance of stellar haloes in cluster galaxies. This paper describes the sample selection, observations, and overall goals of the survey, and provides initial results based on the spectroscopic data, including the detailed characterisation of stellar kinematics and populations to large radii; decomposition of galaxy components directly via their orbital structure; the ability to identify globular clusters and planetary nebulae, and derivation of high-quality emission-line diagnostics in the presence of complex ionised gas.

SECCO1 is an extremely dark, low-mass (M_star=10^5 M_sun), star-forming stellar system lying in the Low Velocity Cloud (LVC) substructure of the Virgo cluster of galaxies, and hosting several HII regions. Here we review our knowledge of this remarkable system, and present the results of (a) additional analysis of our panoramic spectroscopic observations with MUSE, (b) the combined analysis of Hubble Space Telescope and MUSE data, and (c) new narrow-band observations obtained with OSIRIS@GTC to search for additional HII regions in the surroundings of the system. We provide new evidence supporting an age as young as 4 Myr for the stars that are currently ionising the gas in SECCO1. We identify only one new promising candidate HII region possibly associated with SECCO1, thus confirming the extreme isolation of the system. We also identify three additional candidate pressure-supported dark clouds in Virgo among the targets of the SECCO survey. Various possible hypotheses for the nature and origin of SECCO1 are considered and discussed, also with the help of dedicated hydrodynamical simulations showing that a hydrogen cloud with the characteristics of SECCO1 can likely survive for >1 Gyr while traveling within the LVC Intra Cluster Medium.

We use the first data release of the Gaia mission to explore the three dimensional arrangement and the age ordering of the many stellar groups towards the Orion OB association, aiming at a new classification and characterization of the stellar population. We make use of the parallaxes and proper motions provided in the Tycho Gaia Astrometric Solution (TGAS) sub-set of the Gaia catalogue, and of the combination of Gaia and 2MASS photometry. In TGAS we find evidence for the presence of a young population, at a parallax $\varpi \sim 2.65 \, \mathrm{mas}$, loosely distributed around some known clusters: 25 Ori, $\epsilon$ Ori and $\sigma$ Ori, and NGC 1980 ($\iota$ Ori). The low mass counterpart of this population is visible in the color-magnitude diagrams constructed by combining Gaia and 2MASS photometry. We study the density distribution of the young sources in the sky. We find the same groups as in TGAS, and also some other density enhancements that might be related to the recently discovered Orion X group, the Orion dust ring, and to the $\lambda$ Ori complex. We estimate the ages of this population and we infer the presence of an age gradient going from 25 Ori (13-15 Myr) to the ONC (1-2 Myr). We confirm this age ordering by repeating the Bayesian fit using the Pan-STARRS1 data. The estimated ages towards the NGC 1980 cluster span a broad range of values. This can either be due to the presence of two populations coming from two different episodes of star formation or to a large spread along the line of sight of the same population. Our results form the first step towards using the Gaia data to unravel the complex star formation history of the Orion region in terms of the different star formation episodes, their duration, and their effects on the surrounding interstellar medium.

The Galactic Centre is of fundamental astrophysical interest, but existing near-infrared surveys fall short to cover it adequately. Here we introduce the GALACTICNUCLEUS survey, a JHKs imaging survey of the centre of the Milky Way with a 0.2" angular resolution. We present the observations of Field 1 of our survey, centred approximately on SgrA* with an approximate size of 7.95 ' x 3.43 '. We describe the observational set-up and data reduction pipeline and discuss the quality of the data. Finally, we present some preliminary analysis of the data. The data were acquired with the near-infrared camera HAWK-I at the ESO VLT. Short readout times in combination with the speckle holography algorithm allowed us to produce final images with a stable, Gaussian PSF of 0.2" FWHM. Astrometric calibration is achieved via the VVV survey and photometric calibration is based on the SIRIUS/IRSF survey. The quality of the data is assessed by comparison between observations of the same field with different detectors of HAWK-I and at different times. We reach 5 sigma detection limits of approximately J = 22, H = 21, and Ks = 20. The photometric uncertainties are less than 0.05 at J < 20, H < 17 and Ks < 16. We distinguish five stellar populations in the colour-magnitude diagrams; three of them appear to belong to foreground spiral arms, and the other two correspond to a high- and a low-extinction star groups at the Galactic Centre. We use our data to analyse the near-infrared extinction curve and conclude that it can be described very well by a power-law with an index of JHKs = 2.31 +- 0.03. We do not find any evidence that this index depends on the position along the line-of-sight, or on the absolute value of the extinction. We produce extinction maps that show the clumpiness of the ISM at the Galactic Centre. Finally, we estimate that the majority of the stars have solar or super-solar metallicity.

We constrain the possible presence of a central black hole (BH) in the center of the Large Magellanic Cloud (LMC). This requires spectroscopic measurements over an area of order a square degree, due to the poorly known position of the kinematic center. Such measurements are now possible with the impressive field of view of the Multi Unit Spectroscopic Explorer (MUSE) on the ESO Very Large Telescope. We used the Calcium Triplet (~850nm) spectral lines in many short-exposure MUSE pointings to create a two-dimensional integrated-light line-of-sight velocity map from the ~$10^8$ individual spectra, taking care to identify and remove Galactic foreground populations. The data reveal a clear velocity gradient at an unprecedented spatial resolution of 1 arcmin$^{2}$. We fit kinematic models to arrive at a $3\sigma$ upper-mass-limit of $10^{7.1}$ M$_{Sun}$ for any central BH - consistent with the known scaling relations for supermassive black holes and their host systems. This adds to the growing body of knowledge on the presence of BHs in low-mass and dwarf galaxies, and their scaling relations with host-galaxy properties, which can shed light on theories of BH growth and host system interaction.

GRAVITY is a new instrument to coherently combine the light of the European Southern Observatory Very Large Telescope Interferometer to form a telescope with an equivalent 130 m diameter angular resolution and a collecting area of 200 m$^2$. The instrument comprises fiber fed integrated optics beam combination, high resolution spectroscopy, built-in beam analysis and control, near-infrared wavefront sensing, phase-tracking, dual beam operation and laser metrology [...]. This article gives an overview of GRAVITY and reports on the performance and the first astronomical observations during commissioning in 2015/16. We demonstrate phase tracking on stars as faint as m$_K$ ~ 10 mag, phase-referenced interferometry of objects fainter than m$_K$ ~ 15 mag with a limiting magnitude of m$_K$ ~ 17 mag, minute long coherent integrations, a visibility accuracy of better than 0.25 %, and spectro-differential phase and closure phase accuracy better than 0.5{\deg}, corresponding to a differential astrometric precision of better than 10 microarcseconds ({\mu}as). The dual-beam astrometry, measuring the phase difference of two objects with laser metrology, is still under commissioning. First observations show residuals as low as 50 {\mu}as when following objects over several months. We illustrate the instrument performance with the observations of archetypical objects for the different instrument modes. Examples include the Galactic Center supermassive black hole and its fast orbiting star S2 for phase referenced dual beam observations and infrared wavefront sensing, the High Mass X-Ray Binary BP Cru and the Active Galactic Nucleus of PDS 456 for few {\mu}as spectro-differential astrometry, the T Tauri star S CrA for a spectro-differential visibility analysis, {\xi} Tel and 24 Cap for high accuracy visibility observations, and {\eta} Car for interferometric imaging with GRAVITY.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

We construct triaxial dynamical models for the Milky Way nuclear star cluster using Schwarzschild's orbit superposition technique. We fit the stellar kinematic maps presented in Feldmeier et al. (2014). The models are used to constrain the supermassive black hole mass M_BH, dynamical mass-to-light ratio M/L, and the intrinsic shape of the cluster. Our best-fitting model has M_BH = (3.0 +1.1 -1.3)x10^6 M_sun, M/L = (0.90 +0.76 -0.08) M_sun/L_{sun,4.5micron}, and a compression of the cluster along the line-of-sight. Our results are in agreement with the direct measurement of the supermassive black hole mass using the motion of stars on Keplerian orbits. The mass-to-light ratio is consistent with stellar population studies of other galaxies in the mid-infrared. It is possible that we underestimate M_BH and overestimate the cluster's triaxiality due to observational effects. The spatially semi-resolved kinematic data and extinction within the nuclear star cluster bias the observations to the near side of the cluster, and may appear as a compression of the nuclear star cluster along the line-of-sight. We derive a total dynamical mass for the Milky Way nuclear star cluster of M_MWNSC = (2.1 +-0.7)x10^7 M_sun within a sphere with radius r = 2 x r_eff = 8.4 pc. The best-fitting model is tangentially anisotropic in the central r = 0.5-2 pc of the nuclear star cluster, but close to isotropic at larger radii. Our triaxial models are able to recover complex kinematic substructures in the velocity map.

We present the results of VLT-MUSE integral field spectroscopy of SECCO1, a faint, star-forming stellar system recently discovered as the stellar counterpart of an Ultra Compact High Velocity Cloud (HVC274.68+74.0), very likely residing within a substructure of the Virgo cluster of galaxies. We have obtained the radial velocity of a total of 38 individual compact sources identified as HII regions in the main and secondary body of the system, and derived the metallicity for 18 of them. We provide the first direct demonstration that the two stellar bodies of SECCO1 are physically associated and that their velocities match the HI velocities. The metallicity is quite uniform over the whole system, with a dispersion sigma_12+log(O/H/)=0.08, lower than the uncertainty on individual metallicity estimates. The mean abundance, 12+log(O/H)=8.44, is much higher than the typical values for local dwarf galaxies of similar stellar mass. This strongly suggests that the SECCO~1 stars were born from a pre-enriched gas cloud, possibly stripped from a larger galaxy. Using archival HST images we derive a total stellar mass of ~1.6 X 10^5 M_sun for SECCO1, confirming that it has a very high HI to stellar mass ratio for a dwarf galaxy, M_HI/M_*~ 100. The star formation rate, derived from the H_alpha flux is a factor of more than 10 higher than in typical dwarf galaxies of similar luminosity.

Knowing the metallicity distribution of stars in the Galactic Centre has important implications for the formation history of the Milky Way nuclear star cluster. However, this distribution is not well known, and is currently based on a small sample of fewer than 100 stars. We obtained near-infrared K-band spectra of more than 700 late-type stars in the central 4 pc^2 of the Milky Way nuclear star cluster with the integral-field spectrograph KMOS (VLT). We analyse the medium-resolution spectra using a full-spectral fitting method employing the G\"ottingen Spectral library of synthetic PHOENIX spectra. The derived stellar metallicities range from metal-rich [M/H]>+0.3 dex to metal-poor [M/H]<-1.0 dex, with a fraction of 5.2(^{+6.0}+{-3.1}) per cent metal-poor ([M/H]<-0.5 dex) stars. The metal-poor stars are distributed over the entire observed field. The origin of metal-poor stars remains unclear. They could originate from infalling globular clusters. For the metal-rich stellar population ([M/H]>0 dex) a globular cluster origin can be ruled out. As there is only a very low fraction of metal-poor stars in the central 4 pc^2 of the Galactic Centre, we believe that our data can discard a scenario in which the Milky Way nuclear star cluster is purely formed from infalling globular clusters.

We present a 1.4 GHz Karl G. Jansky Very Large Array (VLA) study of a sample of early-type galaxies (ETGs) from the volume- and magnitude-limited ATLAS-3D survey. The radio morphologies of these ETGs at a resolution of 5" are diverse and include sources that are compact on sub-kpc scales, resolved structures similar to those seen in star-forming spiral galaxies, and kpc-scale radio jets/lobes associated with active nuclei. We compare the 1.4 GHz, molecular gas, and infrared (IR) properties of these ETGs. The most CO-rich ATLAS-3D ETGs have radio luminosities consistent with extrapolations from H_2-mass-derived star formation rates from studies of late-type galaxies. These ETGs also follow the radio-IR correlation. However, ETGs with lower molecular gas masses tend to have less radio emission relative to their CO and IR emission compared to spirals. The fraction of galaxies in our sample with high IR-radio ratios is much higher than in previous studies, and cannot be explained by a systematic underestimation of the radio luminosity due to the presence extended, low-surface-brightness emission that was resolved-out in our VLA observations. In addition, we find that the high IR-radio ratios tend to occur at low IR luminosities, but are not associated with low dynamical mass or metallicity. Thus, we have identified a population of ETGs that have a genuine shortfall of radio emission relative to both their IR and molecular gas emission. A number of mechanisms may conspire to cause this radio deficiency, including a bottom-heavy stellar initial mass function, weak magnetic fields, a higher prevalence of environmental effects compared to spirals and enhanced cosmic ray losses.

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

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

The Galactic centre hosts a crowded, dense nuclear star cluster with a half-light radius of 4 pc. Most of the stars in the Galactic centre are cool late-type stars, but there are also >100 hot early-type stars in the central parsec of the Milky Way. These stars are only 3-8 Myr old. Our knowledge of the number and distribution of early-type stars in the Galactic centre is incomplete. Only a few spectroscopic observations have been made beyond a projected distance of 0.5 pc of the Galactic centre. The distribution and kinematics of early-type stars are essential to understand the formation and growth of the nuclear star cluster. We cover the central >4pc^2 of the Galactic centre using the integral-field spectrograph KMOS. We extracted more than 1,000 spectra from individual stars and identified early-type stars based on their spectra. Our data set contains 114 bright early-type stars: 6 have narrow emission lines, 23 are Wolf-Rayet stars, 9 stars have featureless spectra, and 76 are O/B type stars. Our wide-field spectroscopic data confirm that the distribution of young stars is compact, with 90% of the young stars identified within 0.5 pc of the nucleus. We identify 24 new O/B stars primarily at large radii. We estimate photometric masses of the O/B stars and show that the total mass in the young population is >12,000M_sun. The O/B stars all appear to be bound to the Milky Way nuclear star cluster, while less than 30% belong to the clockwise rotating disk. The central concentration of the early-type stars is a strong argument that they have formed in situ. A large part of the young O/B stars is not on the disk, which either means that the early-type stars did not all form on the same disk or that the disk is dissolving rapidly. [abridged]

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

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

NGC1266 is a nearby lenticular galaxy that harbors a massive outflow of molecular gas powered by the mechanical energy of an active galactic nucleus (AGN). It has been speculated that such outflows hinder star formation (SF) in their host galaxies, providing a form of feedback to the process of galaxy formation. Previous studies, however, indicated that only jets from extremely rare, high power quasars or radio galaxies could impart significant feedback on their hosts. Here we present detailed observations of the gas and dust continuum of NGC1266 at millimeter wavelengths. Our observations show that molecular gas is being driven out of the nuclear region at $\dot{M}_{\rm out} \approx 110 M_\odot$ yr$^{-1}$, of which the vast majority cannot escape the nucleus. Only 2 $M_\odot$ yr$^{-1}$ is actually capable of escaping the galaxy. Most of the molecular gas that remains is very inefficient at forming stars. The far-infrared emission is dominated by an ultra-compact ($\lesssim50$pc) source that could either be powered by an AGN or by an ultra-compact starburst. The ratio of the SF surface density ($\Sigma_{\rm SFR}$) to the gas surface density ($\Sigma_{\rm H_2}$) indicates that SF is suppressed by a factor of $\approx 50$ compared to normal star-forming galaxies if all gas is forming stars, and $\approx$150 for the outskirt (98%) dense molecular gas if the central region is is powered by an ultra-compact starburst. The AGN-driven bulk outflow could account for this extreme suppression by hindering the fragmentation and gravitational collapse necessary to form stars through a process of turbulent injection. This result suggests that even relatively common, low-power AGNs are able to alter the evolution of their host galaxies as their black holes grow onto the M-$\sigma$ relation.

Galactic archeology based on star counts is instrumental to reconstruct the past mass assembly of Local Group galaxies. The development of new observing techniques and data-reduction, coupled with the use of sensitive large field of view cameras, now allows us to pursue this technique in more distant galaxies exploiting their diffuse low surface brightness (LSB) light. As part of the Atlas3D project, we have obtained with the MegaCam camera at the Canada-France Hawaii Telescope extremely deep, multi--band, images of nearby early-type galaxies. We present here a catalog of 92 galaxies from the Atlas3D sample, that are located in low to medium density environments. The observing strategy and data reduction pipeline, that achieve a gain of several magnitudes in the limiting surface brightness with respect to classical imaging surveys, are presented. The size and depth of the survey is compared to other recent deep imaging projects. The paper highlights the capability of LSB--optimized surveys at detecting new prominent structures that change the apparent morphology of galaxies. The intrinsic limitations of deep imaging observations are also discussed, among those, the contamination of the stellar halos of galaxies by extended ghost reflections, and the cirrus emission from Galactic dust. The detection and systematic census of fine structures that trace the present and past mass assembly of ETGs is one of the prime goals of the project. We provide specific examples of each type of observed structures -- tidal tails, stellar streams and shells --, and explain how they were identified and classified. We give an overview of the initial results. The detailed statistical analysis will be presented in future papers.

Within the central 10pc of our Galaxy lies a dense nuclear star cluster (NSC), and similar NSCs are found in most nearby galaxies. Studying the structure and kinematics of NSCs reveals the history of mass accretion of galaxy nuclei. Because the Milky Way (MW) NSC is at a distance of only 8kpc, we can spatially resolve the MWNSC on sub-pc scales. This makes the MWNSC a reference object for understanding the formation of all NSCs. We have used the NIR long-slit spectrograph ISAAC (VLT) in a drift-scan to construct an integral-field spectroscopic map of the central 9.5 x 8pc of our Galaxy. We use this data set to extract stellar kinematics both of individual stars and from the unresolved integrated light spectrum. We present a velocity and dispersion map from the integrated light and model these kinematics using kinemetry and axisymmetric Jeans models. We also measure CO bandhead strengths of 1,375 spectra from individual stars. We find kinematic complexity in the NSCs radial velocity map including a misalignment of the kinematic position angle by 9 degree counterclockwise relative to the Galactic plane, and indications for a rotating substructure perpendicular to the Galactic plane at a radius of 20" or 0.8pc. We determine the mass of the NSC within r = 4.2pc to 1.4 x 10^7 Msun. We also show that our kinematic data results in a significant underestimation of the supermassive black hole (SMBH) mass. The kinematic substructure and position angle misalignment may hint at distinct accretion events. This indicates that the MWNSC grew at least partly by the mergers of massive star clusters. Compared to other NSCs, the MWNSC is on the compact side of the r_eff - M_NSC relation. The underestimation of the SMBH mass might be caused by the kinematic misalignment and a stellar population gradient. But it is also possible that there is a bias in SMBH mass measurements obtained with integrated light.

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

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

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

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

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

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

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

We present new high spatial resolution Karl G. Jansky Very Large Array (VLA) HI absorption and Very Long Baseline Array (VLBA) continuum observations of the Active Galactic Nucleus (AGN)-driven molecular outflow candidate NGC 1266. Although other well-known systems with molecular outflows may be driven by star formation in a central molecular disk, the molecular mass outflow rate reported in Alatalo et al. (2011) in NGC 1266 of 13 M$_{\odot}$ year$^{-1}$ exceeds star formation rate estimates from a variety of tracers. This suggests that an additional energy source, such as an AGN, may play a significant role in powering the outflow. Our high spatial resolution HI absorption data reveal compact absorption against the radio continuum core co-located with the putative AGN, and the presence of a blueshifted spectral component re-affirms that gas is indeed flowing out of the system. Our VLBA observations at 1.65 GHz reveal one continuum source within the densest portion of the molecular gas, with a diameter d < 8 mas (1.2 pc), a radio power $P_{\mathrm{rad}}$ = 1.48 $\times$ 10$^{20}$ W Hz$^{-1}$, and a brightness temperature $T_{\mathrm{b}}$ > 1.5 x 10$^7$ K that is most consistent with an AGN origin. The radio continuum energetics implied by the compact VLBA source, as well as archival VLA continuum observations at lower spatial resolution, further support the possibility that the AGN in NGC 1266 could be driving the molecular outflow. These findings suggest that even low-level AGNs may be able to launch massive outflows in their host galaxies.

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

For galaxies hosting supermassive black holes (SMBHs), it has been observed that the mass of the central black hole (M_BH) tightly correlates with the effective or central velocity dispersion (sigma) of the host galaxy. The origin of this M_BH - sigma scaling relation is assumed to lie in the merging history of the galaxies but many open questions about its origin and the behavior in different mass ranges still need to be addressed. The goal of this work is to study the black-hole scaling relations for low black-hole masses, where the regime of intermediate-mass black holes (IMBHs) in globular clusters (GCs) is entered. We collect all existing reports of dynamical black-hole measurements in globular clusters, providing black-hole masses or upper limits for 14 candidates. We plot the black-hole masses versus different cluster parameters including total mass, velocity dispersion, concentration and half-mass radius. We search for trends and test the correlations in order to quantify their significance using a set of different statistical approaches. For correlations showing a large significance we perform a linear fit, accounting for uncertainties and upper limits. We find a clear correlation between the mass of the IMBH and the velocity dispersion of the globular cluster. As expected, the total mass of the globular cluster then also correlates with the mass of the IMBH. While the slope of the M_BH - sigma correlation differs strongly from the one observed for SMBHs, the other scaling relations M_BH - M_TOT, and M_BH - L are similar to the correlations in galaxies. Significant correlations of black-hole mass with other cluster properties were not found in the present sample.

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

Intermediate-mass black holes (IMBHs, 10^2-10^5 M_sun) fill the gap between stellar-mass black holes and supermassive black holes (SMBHs). Simulations have shown that IMBHs may form in dense star clusters, and therefore may still be present in these smaller stellar systems. We investigate the Galactic globular cluster NGC 5286 for indications of a central IMBH using spectroscopic data from VLT/FLAMES, velocity measurements from the Rutgers Fabry Perot at CTIO, and photometric data from HST. We run analytic spherical and axisymmetric Jeans models with different central black-hole masses, anisotropy, mass-to-light ratio, and inclination. Further, we compare the data to a grid of N-body simulations without tidal field. Additionally, we use one N-body simulation to check the results of the spherical Jeans models for the total cluster mass. Both the Jeans models and the N-body simulations favor the presence of a central black hole in NGC 5286 and our detection is at the 1- to 1.5-sigma level. From the spherical Jeans models we obtain a best fit with black-hole mass M_BH=(1.5+-1.0)x10^3 M_sun. The error is the 68% confidence limit from Monte Carlo simulations. Axisymmetric models give a consistent result. The best fitting N-body model is found with a black hole of 0.9% of the total cluster mass (4.38+-0.18)x10^5 M_sun, which results in an IMBH mass of M_BH=(3.9+-2.0)x10^3 M_sun. Jeans models give lower values for the total cluster mass. Our test of the Jeans models with N-body simulation data shows that this discrepancy has two reasons: The influence of a radially varying M/L profile, and underestimation of the velocity dispersion as the measurements are limited to bright stars. We conclude that detection of IMBHs in Galactic globular clusters remains a challenging task unless their mass fractions are above those found for SMBHs in nearby galaxies. [abridged]

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

For a sample of nine Galactic globular clusters we measured the inner kinematic profiles with integral-field spectroscopy that we combined with existing outer kinematic measurements and HST luminosity profiles. With this information we are able to detect the crucial rise in the velocity-dispersion profile which indicates the presence of a central black hole. In addition, N-body simulations compared to our data will give us a deeper insight in the properties of clusters with black holes and stronger selection criteria for further studies. For the first time, we obtain a homogeneous sample of globular cluster integral- field spectroscopy which allows a direct comparison between clusters with and without an intermediate-mass black hole.

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

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

The formation of supermassive black holes at high redshift still remains a puzzle to astronomers. Their growth becomes reasonable only when starting from a massive seed black hole with mass of the order of 10^2 - 10^5 M_SUN. Intermediate-mass black holes (IMBHs) are therefore an important field of research. Especially the possibility of finding them in the centers of globular clusters has recently drawn attention. The search for IMBHs in the centers of globular clusters could therefore shed light on the process of black-hole formation and cluster evolution. We are investigating six galactic globular clusters for the presence of an IMBH at their centers. Based on their kinematic and photometric properties, we selected the globular clusters NGC 1851, NGC 1904 (M79), NGC 5694, NGC 5824, NGC 6093 (M80) and NGC 6266 (M62). We use integral field spectroscopy in order to obtain the central velocity-dispersion profile of each cluster. We compute the cluster photometric center and the surface brightness profile using HST data. After combining these datasets we compare them to analytic Jeans models. We use varying M/L_V profiles for clusters with enough data points in order to reproduce their kinematic profiles in an optimal way. Finally, we vary the mass of the central black hole and test whether the cluster is better fitted with or without an IMBH. We present the statistical significance, including upper limits, of the black-hole mass for each cluster. NGC 1904 and NGC 6266 provide the highest significance for a black hole. Jeans models in combination with a M/L_V profile obtained from N-body simulations (in the case of NGC 6266) predict a central black hole of M_BH = (3 +- 1) x 10^3 M_SUN for NGC 1904 and M_BH = (2 +- 1) x 10^3 M_SUN for NGC 6266. Furthermore, we discuss the possible influence of dark remnants and mass segregation at the center of the cluster on the detection of an IMBH.

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

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

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

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

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

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

We use interferometric CO observations to compare the extent, surface brightness profiles and kinematics of the molecular gas in CO-rich Atlas3D early-type galaxies (ETGs) and spiral galaxies. We find that the molecular gas extent is smaller in absolute terms in ETGs than in late-type galaxies, but that the size distributions are similar once scaled by the galaxies optical/stellar characteristic scale-lengths. Virgo cluster ETGs have less extended molecular gas reservoirs than field counterparts. Approximately half of ETGs have molecular gas surface brightness profiles that follow the stellar light profile. These systems often have relaxed gas out to large radii, suggesting they are unlikely to have had recent merger/accretion events. A third of the sample galaxies show molecular gas surface brightness profiles that fall off slower than the light, and sometimes show a truncation. We suggest that ram pressure stripping and/or the presence of hot gas has compressed/truncated the gas in these systems. The remaining galaxies have rings, or composite profiles, that we argue can be caused by the effects of bars. We investigated the kinematics of the molecular gas using position-velocity diagrams, and compared the observed kinematics with dynamical model predictions, and the observed stellar and ionised gas velocities. We confirm that the molecular gas reaches beyond the turnover of the circular velocity curve in ~70% of our CO-rich Atlas3D ETGs. In general we find that in most galaxies the molecular gas is relaxed and dynamically cold. Molecular gas is a better direct tracer of the circular velocity than the ionised gas, justifying its use as a kinematic tracer for Tully-Fisher and similar analyses. (abridged)

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

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

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

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

We report the detection of five high-velocity stars in the core of the globular cluster NGC 2808. The stars lie on the the red giant branch and show total velocities between 40 and 45 km/s. For a core velocity dispersion sigma_c = 13.4 km/s, this corresponds to up to 3.4 sigma_c. These velocities are close to the estimated escape velocity (~ 50 km/s) and suggest an ejection from the core. Two of these stars have been confirmed in our recent integral field spectroscopy data and we will discuss them in more detail here. These two red giants are located at a projected distance of ~ 0.3 pc from the center. According to their positions on the color magnitude diagram, both stars are cluster members. We investigate several possible origins for the high velocities of the stars and conceivable ejection mechanisms. Since the velocities are close to the escape velocity, it is not obvious whether the stars are bound or unbound to the cluster. We therefore consider both cases in our analysis. We perform numerical simulations of three-body dynamical encounters between binaries and single stars and compare the resulting velocity distributions of escapers with the velocities of our stars. We compare the predictions for a single dynamical encounter with a compact object with those of a sequence of two-body encounters due to relaxation. If the stars are unbound, the encounter must have taken place recently, when the stars were already in the giant phase. After including binary fractions and black-hole retention fractions, projection effects, and detection probabilities from Monte-Carlo simulations, we estimate the expected numbers of detections for all the different scenarios. Based on these numbers, we conclude that the most likely scenario is that the stars are bound and were accelerated by a single encounter between a binary of main-sequence stars and a ~ 10 M_sun black hole.

Globular clusters are an excellent laboratory for stellar population and dynamical research. Recent studies have shown that these stellar systems are not as simple as previously assumed. With multiple stellar populations as well as outer rotation and mass segregation they turn out to exhibit high complexity. This includes intermediate-mass black holes which are proposed to sit at the centers of some massive globular clusters. Today's high angular resolution ground based spectrographs allow velocity-dispersion measurements at a spatial resolution comparable to the radius of influence for plausible IMBH masses, and to detect changes in the inner velocity-dispersion profile. Together with high quality photometric data from HST, it is possible to constrain black-hole masses by their kinematic signatures. We determine the central velocity-dispersion profile of the globular cluster NGC 2808 using VLT/FLAMES spectroscopy. In combination with HST/ACS data our goal is to probe whether this massive cluster hosts an intermediate-mass black hole at its center and constrain the cluster mass to light ratio as well as its total mass. We derive a velocity-dispersion profile from integral field spectroscopy in the center and Fabry Perot data for larger radii. High resolution HST data are used to obtain the surface brightness profile. Together, these data sets are compared to dynamical models with varying parameters such as mass to light ratio profiles and black-hole masses. Using analytical Jeans models in combination with variable M/L profiles from N-body simulations we find that the best fit model is a no black hole solution. After applying various Monte Carlo simulations to estimate the uncertainties, we derive an upper limit of the back hole mass of M_BH < 1 x 10^4 M_SUN (with 95 % confidence limits) and a global mass-to-light ratio of M/L_V = (2.1 +- 0.2) M_SUN/L_SUN.

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

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

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

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

Supermassive black holes (SMBHs) are fundamental keys to understand the formation and evolution of their host galaxies. However, the formation and growth of SMBHs are not yet well understood. One of the proposed formation scenarios is the growth of SMBHs from seed intermediate-mass black holes (IMBHs, 10^2 to 10^5 M_{\odot}) formed in star clusters. In this context, and also with respect to the low mass end of the M-sigma relation for galaxies, globular clusters are in a mass range that make them ideal systems to look for IMBHs. Among Galactic star clusters, the massive cluster $\omega$ Centauri is a special target due to its central high velocity dispersion and also its multiple stellar populations. We study the central structure and dynamics of the star cluster $\omega$ Centauri to examine whether an IMBH is necessary to explain the observed velocity dispersion and surface brightness profiles. We perform direct N-body simulations to follow the dynamical evolution of $\omega$ Centauri. The simulations are compared to the most recent data-sets in order to explain the present-day conditions of the cluster and to constrain the initial conditions leading to the observed profiles. We find that starting from isotropic spherical multi-mass King models and within our canonical assumptions, a model with a central IMBH mass of 2% of the cluster stellar mass, i.e. a 5x10^4 M_{\odot} IMBH, provides a satisfactory fit to both the observed shallow cusp in surface brightness and the continuous rise towards the center of the radial velocity dispersion profile. In our isotropic spherical models, the predicted proper motion dispersion for the best-fit model is the same as the radial velocity dispersion one. (abridged)

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

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

Grigori Kuzmin was a very gifted dynamicist and one of the towering figures in the distinguished history of the Tartu Observatory. He obtained a number of important results in relative isolation which were later rediscovered in the West. This work laid the foundation for further advances in the theory of stellar systems in dynamical equilibrium, thereby substantially increasing our understanding of galaxy dynamics.

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

Intermediate-mass black holes (IMBHs) are of interest in a wide range of astrophysical fields. In particular, the possibility of finding them at the centers of globular clusters has recently drawn attention. IMBHs became detectable since the quality of observational data sets, particularly those obtained with HST and with high resolution ground based spectrographs, advanced to the point where it is possible to measure velocity dispersions at a spatial resolution comparable to the size of the gravitational sphere of influence for plausible IMBH masses. We present results from ground based VLT/FLAMES spectroscopy in combination with HST data for the globular cluster NGC 6388. The aim of this work is to probe whether this massive cluster hosts an intermediate-mass black hole at its center and to compare the results with the expected value predicted by the $M_{\bullet} - \sigma$ scaling relation. The spectroscopic data, containing integral field unit measurements, provide kinematic signatures in the center of the cluster while the photometric data give information of the stellar density. Together, these data sets are compared to dynamical models and present evidence of an additional compact dark mass at the center: a black hole. Using analytical Jeans models in combination with various Monte Carlo simulations to estimate the errors, we derive (with 68% confidence limits) a best fit black-hole mass of $ (17 \pm 9) \times 10^3 M_{\odot}$ and a global mass-to-light ratio of $M/L_V = (1.6 \pm 0.3) \ M_{\odot}/L_{\odot}$.

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

We make use of interferometric CO and HI observations, and optical integral-field spectroscopy to probe the origin of the molecular and ionised interstellar medium (ISM) in local early-type galaxies (ETGs). We find that 36\pm5% of our sample of fast rotating ETGs have their ionised gas kinematically misaligned with respect to the stars, setting a strong lower limit on the importance of externally acquired gas (e.g. from mergers and cold accretion). Slow rotators have a flat distribution of misalignments, indicating that the dominant source of gas is external. The molecular, ionised and atomic gas in all the detected galaxies are always kinematically aligned, even when they are misaligned from the stars, suggesting that all these three phases of the ISM share a common origin. In addition, we find that the origin of the cold and warm gas in fast-rotating ETGs is strongly affected by environment, despite the molecular gas detection rate and mass fractions being fairly independent of group/cluster membership. Galaxies in dense groups and the Virgo cluster nearly always have their molecular gas kinematically aligned with the stellar kinematics, consistent with a purely internal origin. In the field, however, kinematic misalignments between the stellar and gaseous components indicate that >46% of local fast-rotating ETGs have their gas supplied from external sources. We discuss several scenarios which could explain the environmental dichotomy, but find it difficult to simultaneously explain the kinematic misalignment difference and the constant detection rate. Furthermore, our results suggest that galaxy mass may be an important independent factor associated with the origin of the gas, with the most massive fast-rotating galaxies in our sample (M_K<-24 mag; stellar mass of >8x10^10 Msun) always having kinematically aligned gas. (abridged)

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

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

We look at the morphology of fast and slow rotator early-type galaxies. Edge-on fast rotators are lenticular galaxies. They appear like spiral galaxies with the gas and dust removed, and in some cases are flat ellipticals with disky isophotes. Fast rotators are often barred and span the same full range of bulge fractions as spiral galaxies. The slow rotators are rounder and are generally consistent with being genuine, namely spheroidal-like, elliptical galaxies. We propose a revision to the tuning-fork diagram by Hubble as it gives a misleading description of ETGs. We study for the first time the kinematic morphology-density T-Sigma relation using fast and slow rotators to replace lenticulars and ellipticals. We find that our relation is cleaner than using classic morphology. Slow rotators are nearly absent at the lowest density environments [f(SR)<2%] and generally constitute a small fraction [f(SR)~4%] of the total galaxy population in the environments explored by our survey, with the exception of the densest core of the Virgo cluster [f(SR)~20%]. We find a clean log-linear relation between the fraction f(Sp) of spiral galaxies and the local galaxy surface density. The existence of a smooth kinematic T-Sigma relation in the field excludes processes related to the cluster environment as main contributors to the apparent conversion of spirals into fast-rotators in low-density environments. It shows that the segregation is driven by local effects at the small-group scale. Only at the largest densities in the Virgo core does the f(Sp) relation break down and steepens sharply, while the fraction of slow-rotators starts to significantly increase. This suggests that a different mechanism is at work there. (Abridged)

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

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

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

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

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

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

Boxy and peanut-shaped bulges are seen in about half of edge-on disc galaxies. Comparisons of the photometry and major-axis gas and stellar kinematics of these bulges to simulations of bar formation and evolution indicate that they are bars viewed in projection. If the properties of boxy bulges can be entirely explained by assuming they are bars, then this may imply that their hosts are pure disc galaxies with no classical bulge. A handful of these bulges, including that of the Milky Way, have been observed to rotate cylindrically, i.e. with a mean stellar velocity independent of height above the disc. In order to assess whether such behaviour is ubiquitous in boxy bulges, and whether a pure disc interpretation is consistent with their stellar populations, we have analysed the stellar kinematics and populations of the boxy or peanut-shaped bulges in a sample of five edge-on galaxies. We placed slits along the major axis of each galaxy and at three offset but parallel positions to build up spatial coverage. The boxy bulge of NGC3390 rotates perfectly cylindrically within the spatial extent and uncertainties of the data. This is consistent with the metallicity and alpha-element enhancement of the bulge, which are the same as in the disk. This galaxy is thus a pure disc galaxy. The boxy bulge of ESO311-G012 also rotates very close to cylindrically. The boxy bulge of NGC1381 is neither clearly cylindrically nor non-cylindrically rotating, but it has a negative vertical metallicity gradient and is alpha-enhanced with respect to its disc, suggesting a composite bulge comprised of a classical bulge and bar (and possibly a discy pseudobulge) [abridged] Even this relatively small sample is sufficient to demonstrate that boxy bulges display a range of rotational and population properties, indicating that they do not form a homogeneous class of object.

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

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

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

The Atlas3D project is a multi-wavelength survey combined with a theoretical modeling effort. The observations span from the radio to the millimeter and optical, and provide multi-colour imaging, two-dimensional kinematics of the atomic (HI), molecular (CO) and ionized gas (Hbeta, [OIII] and [NI]), together with the kinematics and population of the stars (Hbeta, Fe5015 and Mgb), for a carefully selected, volume-limited (1.16*10^5 Mpc^3) sample of 260 early-type (elliptical E and lenticular S0) galaxies (ETGs). The models include semi-analytic, N-body binary mergers and cosmological simulations of galaxy formation. Here we present the science goals for the project and introduce the galaxy sample and the selection criteria. The sample consists of nearby (D<42 Mpc) morphologically-selected ETGs extracted from a parent sample of 871 galaxies (8% E, 22% S0 and 70% spirals) brighter than M_K<-21.5 mag (stellar mass M_Star>6*10^9 M_Sun). We analyze possible selection biases and we conclude that the parent sample is essentially complete and statistically representative of the nearby galaxy population. We present the size-luminosity relation for the spirals and ETGs and show that the ETGs in the Atlas3D sample define a tight red sequence in a colour-magnitude diagram, with few objects in the transition from the blue cloud. We describe the strategy of the SAURON integral-field observations and the extraction of the stellar kinematics with the pPXF method. We give an overview of the characteristics of the other main datasets already available for our sample and of the ongoing modelling projects.

The Galactic globular cluster omega Centauri is a prime candidate for hosting an intermediate mass black hole. Recent measurements lead to contradictory conclusions on this issue. We use VLT-FLAMES to obtain new integrated spectra for the central region of omega Centauri. We combine these data with existing measurements of the radial velocity dispersion profile taking into account a new derived center from kinematics and two different centers from the literature. The data support previous measurements performed for a smaller field of view and show a discrepancy with the results from a large proper motion data set. We see a rise in the radial velocity dispersion in the central region to 22.8+-1.2 km/s, which provides a strong sign for a central black hole. Isotropic dynamical models for omega Centauri imply black hole masses ranging from 3.0 to 5.2x10^4 solar masses depending on the center. The best-fitted mass is 4.7+-1.0x10^4 solar masses.

(Abridged) We present a stellar population analysis of the absorption line strength maps for 48 early-type galaxies from the SAURON sample. Using the line strength index maps of Hbeta, Fe5015, and Mgb, measured in the Lick/IDS system and spatially binned to a constant signal-to-noise, together with predictions from up-to-date stellar population models, we estimate the simple stellar population-equivalent (SSP-equivalent) age, metallicity and abundance ratio [alpha/Fe] over a two-dimensional field extending up to approximately one effective radius. We find a large range of SSP-equivalent ages in our sample, of which ~40% of the galaxies show signs of a contribution from a young stellar population. The most extreme cases of post-starburst galaxies, with SSP-equivalent ages of <=3 Gyr observed over the full field-of-view, and sometimes even showing signs of residual star-formation, are restricted to low mass systems(sigma_e <= 100 k/ms or ~2x10^10 M_sol). Spatially restricted cases of young stellar populations in circumnuclear regions can almost exclusively be linked to the presence of star-formation in a thin, dusty disk/ring, also seen in the near-UV or mid-IR. The flattened components with disk-like kinematics previously identified in all fast rotators (Krajnovi\'c et al.) are shown to be connected to regions of distinct stellar populations. These range from the young, still star-forming circumnuclear disks and rings with increased metallicity preferentially found in intermediate-mass fast rotators, to apparently old structures with extended disk-like kinematics, which are observed to have an increased metallicity and mildly depressed [alpha/Fe] ratio compared to the main body of the galaxy. The slow rotators generally show no stellar population signatures over and above the well known metallicity gradients and are largely consistent with old (>=10 Gyr) stellar populations.

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

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

We report observations of Lyman Alpha Blob 1 (LAB1) in the SSA 22 protocluster region (z=3.09) with the integral-field spectrograph SAURON. We increased the signal-to-noise in the spectra by more than a factor three compared to our previous observations. This allows us to probe the structure of the LAB system in detail, examining its structure in the spatial and wavelength dimensions. We find that the emission from the system comes largely from five distinct blobs. Two of the emission regions are associated with Lyman Break Galaxies, while a third appears to be associated with a heavily obscured submillimeter galaxy. The fourth and fifth components do not appear to be associated with any galaxy despite the deep imaging that is available in this field. If we interpret wavelength shifts in the line centroid as velocity structure in the underlying gas, many of these emission systems show evidence of velocity shear. It remains difficult to distinguish between an underlying rotation of the gas and an outflow driven by the central object. We have examined all of the line profiles for evidence of strong absorption features. While several systems are better fitted by the inclusion of a weak absorption component, we do not see evidence for a large-scale coherent absorption feature such as that seen in LAB2.

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

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

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

We present HI observations of the elliptical galaxy NGC 2974, obtained with the Very Large Array. These observations reveal that the previously detected HI disc in this galaxy (Kim et al. 1988) is in fact a ring. By studying the harmonic expansion of the velocity field along the ring, we constrain the elongation of the halo and find that the underlying gravitational potential is consistent with an axisymmetric shape. We construct mass models of NGC 2974 by combining the HI rotation curve with the central kinematics of the ionised gas, obtained with the integral-field spectrograph SAURON. We introduce a new way of correcting the observed velocities of the ionised gas for asymmetric drift, and hereby disentangle the random motions of the gas caused by gravitational interaction from those caused by turbulence. To reproduce the observed flat rotation curve of the HI gas, we need to include a dark halo in our mass models. A pseudo-isothermal sphere provides the best model to fit our data, but we also tested an NFW halo and Modified Newtonian Dynamics (MOND), which fit the data marginally worse. The mass-to-light ratio M/L_I increases in NGC 2974 from 4.3 (M/L_I)sun at one effective radius to 8.5 (M/L_I)sun at 5 Re. This increase of M/L already suggests the presence of dark matter: we find that within 5 Re at least 55 per cent of the total mass is dark.

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

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

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

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

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

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

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

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

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

As part of an effort to enlarge the number of well-studied Magellanic-type galaxies, we obtained broadband optical imaging and neutral hydrogen radio synthesis observations of the dwarf irregular galaxy ESO 364-G 029. The optical morphology characteristically shows a bar-like main body with a one-sided spiral arm, an approximately exponential light distribution, and offset photometric and kinematic centers. The HI distribution is mildly asymmetric and, although slightly offset from the photometric center, roughly follows the optical brightness distribution, extending to over 1.2 Holmberg radii (where mu_B = 26.5 mag/arcsec^2). In particular, the highest HI column densities closely follow the bar, one-arm spiral, and a third optical extension. The rotation is solid-body in the inner parts but flattens outside of the optical extent. The total HI flux F_HI = 23.1 pm 1.2 Jy km/s, yielding a total HI mass M_HI= (6.4 pm 1.7) x 10^8 Msun (for a distance D = 10.8 pm 1.4 Mpc) and a total HI mass-to-blue-luminosity ratio M_HI/L_B = (0.96 pm 0.14) Msun / Lsun,B (distance independent). The HI data suggest a very complex small-scale HI structure, with evidence of large shells and/or holes, but deeper observations are required for a detailed study. Follow-up observations are also desirable for a proper comparison with the Large Magellanic Cloud, where despite an optical morphology very similar to ESO 364-G 029 the HI bears little resemblance to the optical.

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

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

We are conducting the first systematic 3D spectroscopic imaging survey to quantify the properties of the atomic gas (HI) in a distance-limited sample of 28 Seyfert galaxies and a sample of 28 inactive control galaxies with well-matched optical properties (the VHIKINGS survey). This study aims to address the role of the host galaxy in nuclear activity and confront outstanding controversies in optical/IR imaging surveys. Early results show possible relationships between Seyfert activity and HI extent, content and the prevalence of small, nearby gas-rich dwarf galaxies (M(HI)~10^7 Msun); results will be tested via rigorous comparison with control galaxies. Initial results from our optical followup study of 15 of our galaxies using the SAURON integral field unit on the WHT suggest a possible difference between Seyfert and inactive stellar and gaseous kinematics that support the conclusion that internal kinematics of galaxies are the key to nuclear activity.

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

The Multi Unit Spectroscopic Explorer (MUSE) is a second-generation VLT panoramic integral-field spectrograph under preliminary design study. MUSE has a field of 1x1 arcmin**2 sampled at 0.2x0.2 arcsec**2 and is assisted by the VLT ground layer adaptive optics ESO facility using four laser guide stars. The simultaneous spectral range is 465-930 nm, at a resolution of R~3000. MUSE couples the discovery potential of a large imaging device to the measuring capabilities of a high-quality spectrograph, while taking advantage of the increased spatial resolution provided by adaptive optics. This makes MUSE a unique and tremendously powerful instrument for discovering and characterizing objects that lie beyond the reach of even the deepest imaging surveys. MUSE has also a high spatial resolution mode with 7.5x7.5 arcsec**2 field of view sampled at 25 milli-arcsec. In this mode MUSE should be able to obtain diffraction limited data-cubes in the 600-930 nm wavelength range. Although the MUSE design has been optimized for the study of galaxy formation and evolution, it has a wide range of possible applications; e.g. monitoring of outer planets atmosphere, environment of young stellar objects, super massive black holes and active nuclei in nearby galaxies or massive spectroscopic surveys of stellar fields in the Milky Way and nearby galaxies.

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

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

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

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

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

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

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

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

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

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

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

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

High-velocity galactic outflows, driven by intense bursts of star formation and black hole accretion, are invoked by current theories of galaxy formation to terminate star formation in the most massive galaxies and to deposit heavy elements in the intergalactic medium. From existing observational evidence on high-redshift galaxies, it is unclear whether such outflows are localized to regions of intense star formation just a few kiloparsecs in extent, or whether they instead have a significant impact on the entire galaxy and its surroundings. Here we present two-dimensional spectroscopy of a star-forming galaxy at redshift z=3.09 (seen 11.5 Gyr ago, when the Universe was 20 per cent of its current age): its spatially extended Ly-alpha emission appears to be absorbed by HI in a foreground screen covering the entire galaxy, with a lateral extent of at least 100 kpc and remarkable velocity coherence. It was plausibly ejected from the galaxy during a starburst several 1E8 yr earlier and has subsequently swept up gas from the surrounding intergalactic medium and cooled. This demonstrates the galaxy-wide impact of high-redshift superwinds.

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

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

We have used the SAURON panoramic integral field spectrograph to study the structure of the Ly-alpha emission-line halo, LAB1, surrounding the sub-millimeter galaxy SMM J221726+0013. This emission-line halo was discovered during a narrow-band imaging survey of the z=3.1 large-scale structure in the SSA22 region. Our observations trace the emission halo out to almost 100 kpc from the sub-millimeter source and identify two distinct Ly-alpha `mini-haloes' around the nearby Lyman-break galaxies. The main emission region has a broad line profile, with variations in the line profile seeming chaotic and lacking evidence for a coherent velocity structure. The data also suggests that Ly-alpha emission is suppressed around the sub-mm source. Interpretation of the line structure needs care because Ly-alpha may be resonantly scattered, leading to complex radiative transfer effects, and we suggest that the suppression in this region arises because of such effects. We compare the structure of the central emission-line halo with local counter-parts, and find that the emission line halo around NGC 1275 in the Perseus cluster may be a good local analogue, although the high redshift halo is factor of ~100 more luminous and appears to have higher velocity broadening. Around the Lyman-break galaxy C15, the emission line is narrower, and a clear shear in the emission wavelength is seen. A plausible explanation for the line profile is that the emission gas is expelled from C15 in a bipolar outflow, similar to that seen in M82.

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

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

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

We present the results of the analysis of a set of medium resolution spectra, obtained by the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope, of the emission line gas present in the nuclei of a complete sample of 21 nearby, early-type galaxies with radio jets (the UGC FR-I Sample). For each galaxy nucleus we present spectroscopic data in the region of H-alpha and the dervived kinematics. We find that in 67% of the nuclei the gas appears to be rotating and, with one exception, the cases where rotation is not seen are either face on or have complex central morphologies. We find that in 62% of the nuclei the fit to the central spectrum is improved by the inclusion of a broad component. The broad components have a mean velocity dispersion of 1349 +/- 345 km\s and are redshifted from the narrow line components (assuming an origin in H-alpha) by 486 +/- 443 km\s.

We present radial velocities for 29 B- and A-type stars in the field of the nearby association Perseus OB2. The velocities are derived from spectra obtained with AURELIE, via cross correlation with radial velocity standards matched as closely as possible in spectral type. The resulting accuracy is ~2 - 3 km s$^{-1}$. We use these measurements, together with published values for a few other early-type stars, to study membership of the association. The mean radial velocity (and measured velocity dispersion) of Per OB2 is 23.5 \pm 3.9 km s$^{-1}$, and lies ~15 km s$^{-1}$ away from the mean velocity of the local disk field stars. We identify a number of interlopers in the list of possible late-B- and A-type members which was based on Hipparcos parallaxes and proper motions, and discuss the colour-magnitude diagram of the association.

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.

The Jeans equations relate the second-order velocity moments to the density and potential of a stellar system. For general three-dimensional stellar systems, there are three equations and six independent moments. By assuming that the potential is triaxial and of separable Staeckel form, the mixed moments vanish in confocal ellipsoidal coordinates. Consequently, the three Jeans equations and three remaining non-vanishing moments form a closed system of three highly-symmetric coupled first-order partial differential equations in three variables. These equations were first derived by Lynden-Bell, over 40 years ago, but have resisted solution by standard methods. We present the general solution here. We consider the two-dimensional limiting cases first. We solve their Jeans equations by a new method which superposes singular solutions. The singular solutions, which are new, are standard Riemann-Green functions. The two-dimensional solutions are applied to non-axisymmetric discs, oblate and prolate spheroids, and also to the scale-free triaxial limit. We then extend the method of singular solutions to the triaxial case, and obtain a full solution. The general solution can be expressed in terms of complete (hyper)elliptic integrals which can be evaluated in a straightforward way, and provides the full set of second moments which can support a triaxial density distribution in a separable triaxial potential. (abridged)

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

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

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

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

We modify the curvature condition for the existence of self-consistent scale-free discs, introduced by Zhao, Carollo & de Zeeuw. We survey the parameter space of the power-law discs, and show that the modified curvature condition is in harmony with the results of Schwarzschild's numerical orbit superposition method. We study the orbital structure of the power-law discs, and find a correlation between the population of centrophobic banana orbits and the non-self-consistency index. We apply the curvature condition to other families of scale-free elongated discs and find that it rules out a large range of power-law slopes and axis ratios. We generalize the condition, and apply it, to three-dimensional scale-free axisymmetric galaxy models.

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

The GAIA space observatory was recently approved as Cornerstone 6 of ESA's science program, to be launched no later than mid-2012. It will provide a stereoscopic and kinematic census of about 10^9 stars throughout our Galaxy (and into the Local Group) complete to V=20 mag, amounting to about 1% of the Galactic stellar population. Combined with astrophysical information for each star, provided by on-board multi-colour photometry and (limited) spectroscopy, the positional and radial velocity measurements will have the precision necessary to quantify the early formation, and subsequent dynamical, chemical and star formation evolution of our Galaxy: GAIA will establish when the stars in our Galaxy formed, when and how the Galaxy was assembled, and how the dark matter is distributed. The GAIA data will also allow detection and orbital classification of ~5 x 10^4 extra-Solar planetary systems, provide a comprehensive survey of 10^5-10^6 minor bodies in our Solar System, of galaxies in the nearby Universe, of some 5 x 10^5 quasars, and will test general relativity and cosmology.

We describe a practical method for constructing axisymmetric two-integral galaxy models (with distribution functions of the form f(E,L_z), in which E is the orbital energy, and L_z is the vertical component of the angular momentum), based on Schwarzschild's orbit superposition method. Other f(E,L_z)-methods are mostly based on solving the Jeans equations or on finding the distribution function directly from the density, which often places restrictions on the shape of the galaxy. Here, no assumptions are made and any axisymmetric density distribution is possible. The observables are calculated (semi-)analytically, so that our method is faster than most previous, fully numerical implementations. Various aspects are tested extensively, the results of which apply directly to three-integral Schwarzschild methods. We show that a given distribution function can be reproduced with high accuracy and investigate the behaviour of the parameter that is used to measure the goodness-of-fit. Furthermore, we show that the method correctly identifies the range of cusp clopes for which axisymmetric two-integral models with a central black hole do not exist.

We present surface brightness profiles for 56 of the 78 spiral galaxies observed in the HST/{\tt NICMOS/Camera-2} F160W snapshot survey introduced in Carollo et al. (2001; paper I), as well as surface brightness profiles for 23 objects out of the 41 that were also observed in the F110W filter. We fit these surface brightness profiles with the Nuker law of Lauer et al. (1995), and use the smooth analytical descriptions of the data to compute the average nuclear stellar cusp slopes $<\gamma>$ in the 0.1$''-0.5''$ radial range. Our main result is the startling similarity between the nuclear stellar cusp slopes $<\gamma>$ in the near-infared compared to those derived in the visual passband. This similarity has several implications: (1.) There are typically little or no optical-NIR {\it global} color gradients, and thus no global stellar population variations, inside $\sim 50$-100 pc from the nucleus in nearby spirals. (2.) The large observed range of the strength of the nuclear stellar cusps seen in the HST optical study of spiral galaxies reflects a physical difference between galaxies, and is not an artifact caused by nuclear dust and/or recent star formation. (3.) The dichotomy between $R^{1/4}$ bulges, with steep nuclear stellar cusps $<\gamma > \sim 1$, and exponential bulges, with shallow nuclear stellar cusps $<\gamma > < 0.3$, is also not an artifact of the effects of dust or recent star formation. (4) A surrounding massive disk appears to have no effect on the rise of the stellar density distribution within the innermost hundred pc of the

This is the first of two papers where we present the analysis of an HST NICMOS-Cam2 near-infared (NIR) snapshot survey in the F160W (H) filter for a sample of 78 spiral galaxies selected from the UGC and ESOLV catalogs. For 69 of these objects we provide nuclear color information derived by combining the H data either with additional NICMOS F110W (J) images or with V WFPC2/HST data. Here we present the NIR images and the optical-NIR color maps. We focus our attention on the properties of the photometrically-distinct `nuclei' which are found embedded in most of the galaxies, and provide measurements of their half-light radii and magnitudes in the H (and when available, in the J) band. We find that: (i) In the NIR, the nuclei embedded in the bright early- to intermediate-type galaxies span a much larger range in brightness than the nuclei which are typically found embedded in bulgeless late-type disks; (ii) Nuclei are found in both non-barred and barred hosts, in large-scale ($\gta$ 1 kpc) as well as in nuclear (up to a few 100pc) bars; (iii) There is a significant increase in half-light radius with increasing luminosity of the nucleus in the early/intermediate types; (iv) The nuclei of early/intermediate-type spirals cover a large range of optical-NIR colors, from $V-H\approx$ -0.5 to 3. Some nuclei are bluer and others redder than the surrounding galaxy, indicating the presence of activity or reddening by dust in many of these systems; (v) Some early/intermediate nuclei are elongated and/or slightly offset from the isophotal center of the host galaxy. On average, however, these nuclei appear as centered, star-cluster-like structures similar to those which are found in the late-type disks. [abridged]

Early results are reported from the SAURON survey of the kinematics and stellar populations of a representative sample of nearby E, S0 and Sa galaxies. The survey is aimed at determining the intrinsic shape of the galaxies, their orbital structure, the mass-to-light ratio as a function of radius, the age and metallicity of the stellar populations, and the frequency of kinematically decoupled cores and nuclear black holes. The construction of the representative sample is described, and its properties are illustrated. A comparison with long-slit spectroscopic data establishes that the SAURON measurements are comparable to, or better than, the highest-quality determinations. Comparisons are presented for NGC 3384 and NGC 4365 where stellar velocities and velocity dispersions are determined to a precision of 6 km/s, and the h3 and h4 parameters of the line-of-sight velocity distribution to a precision of better than 0.02. Extraction of accurate gas emission-line intensities, velocities and line widths from the datacubes is illustrated for NGC 5813. Comparisons with published line-strengths for NGC 3384 and NGC 5813 reveal uncertainties of < 0.1 A on the measurements of the Hbeta, Mgb and Fe5270 indices. Integral-field mapping uniquely connects measurements of the kinematics and stellar populations to the galaxy morphology. The maps presented here illustrate the rich stellar kinematics, gaseous kinematics, and line-strength distributions of early-type galaxies. The results include the discovery of a thin, edge-on, disk in NGC 3623, confirm the axisymmetric shape of the central region of M32, illustrate the LINER nucleus and surrounding counter-rotating star-forming ring in NGC 7742, and suggest a uniform stellar population in the decoupled core galaxy NGC 5813.

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

We describe a new technique, kinemetry, to quantify kinematic maps of early-type galaxies in an efficient way. We present the first applications to velocity fields obtained with the integral-field spectrograph SAURON.

We describe MIFS, a second generation integral-field spectrograph for the VLT, operating in the visible wavelength range. It combines a 1'x1' field of view with the improved spatial resolution provided by multi-conjugate adaptive optics and covers a large simultaneous spectral range (0.6-1.0 microns). A separate mode exploits the highest spatial resolution provided by adaptive optics. With this unique combination of capabilities, MIFS has a wide domain of application and a large discovery potential. The MIFS low-spatial resolution mode (sampled at 0.2") combined with the initial MCAO capabilities planned for the VLT will provide ultra-deep fields with a limiting magnitude for spectroscopy of R~28. MIFS will improve the present day detection limit of Lyman alpha emitters by a factor of 100, and will detect low-mass star-forming galaxies to z~7. The MIFS high-spatial resolution mode (3"x3" field sampled at 0.01") is optimized for the next step in (MC)AO. It will probe, e.g., the relationship between supermassive central black holes and their host galaxy and the physics of winds from accretion disks in young stellar objects at unprecedented spatial resolution. MIFS will extend Europe's lead in integral-field spectroscopy. It capitalizes on new developments in adaptive optics, and is a key step towards instrumentation for OWL.

We construct phase-space distribution functions for the oblate, cuspy mass models of Sridhar & Touma, which may contain a central point mass (black hole) and have potentials of St\"ackel form in parabolic coordinates. The density in the ST models is proportional to a power $r^{-\gamma}$ of the radius, with $0<\gamma<1$. We derive distribution functions $f(E, L_z)$ for the scale-free ST models (no black hole) using a power series of the energy $E$ and the component $L_z$ of the angular momentum parallel to the symmetry axis. We use the contour integral method of Hunter & Qian to construct $f(E, L_z)$ for ST models with central black holes, and employ the scheme introduced by Dejonghe & de Zeeuw to derive more general distribution functions which depend on $E$, $L_z$ and the exact third integral $I_3$. We find that self-consistent two- and three-integral distribution functions exist for all values $0 < \gamma < 1$.

Early results are reported from the SAURON survey of the kinematics and stellar populations of nearby cluster and field E, S0, and Sa galaxies. We present maps of the Hbeta and [OIII] emission-line distribution and kinematics for NGC 5813 and NGC 7742.

We present initial results of our analysis of line emission produced in gas disks found at the centers of a sample of nearby, radio galaxies with radio jets. We obtained data using STIS (The Space Telescope Imaging Spectrograph) at three parallel slit positions on the nucleus of each galaxy. This allows us to map the H$\alpha$ + [NII] flux, the gas radial velocity and the velocity dispersion. We find evidence of rotating disks in 11 of the sample galaxies and we can not currently rule out a rotating disk model for the remaining eight. For rotating systems, we find that the minimum central enclosed mass is greater than or similar to the predicted black hole mass based on ground-based stellar velocity dispersions. By modeling the gas dynamics we will go on to constrain the masses of the black holes. We will also investigate the properties of the gas disks themselves, giving us an insight into fueling, ionization mechanisms and the structure of the central regions.

A new integral-field spectrograph, SAURON, is described. It is based on the Tiger principle, and uses a lenslet array. SAURON has a large field of view and high throughput, and allows simultaneous sky subtraction. Its design is optimized for studies of the stellar kinematics, gas kinematics, and line-strength distributions of nearby early-type galaxies. The instrument design and specifications are described, as well as the extensive analysis software which was developed to obtain fully calibrated spectra, and the associated kinematic and line-strength measurements. A companion paper reports on the first results obtained with SAURON on the William Herschel Telescope.

The GAIA astrometric mission has recently been approved as one of the next two `cornerstones' of ESA's science programme, with a launch date target of not later than mid-2012. GAIA will provide positional and radial velocity measurements with the accuracies needed to produce a stereoscopic and kinematic census of about one billion stars throughout our Galaxy (and into the Local Group), amounting to about 1 per cent of the Galactic stellar population. GAIA's main scientific goal is to clarify the origin and history of our Galaxy, from a quantitative census of the stellar populations. It will advance questions such as when the stars in our Galaxy formed, when and how it was assembled, and its distribution of dark matter. The survey aims for completeness to V=20 mag, with accuracies of about 10 microarcsec at 15 mag. Combined with astrophysical information for each star, provided by on-board multi-colour photometry and (limited) spectroscopy, these data will have the precision necessary to quantify the early formation, and subsequent dynamical, chemical and star formation evolution of our Galaxy. Additional products include detection and orbital classification of tens of thousands of extra-Solar planetary systems, and a comprehensive survey of some 10^5-10^6 minor bodies in our Solar System, through galaxies in the nearby Universe, to some 500,000 distant quasars. It will provide a number of stringent new tests of general relativity and cosmology. The complete satellite system was evaluated as part of a detailed technology study, including a detailed payload design, corresponding accuracy assesments, and results from a prototype data reduction development.

We have observed a complete sample of 21 nearby (D < 70h^{-1} Mpc) Fanaroff & Riley Type I galaxies with HST/WFPC2 and detected dust disks and lanes in 19 of them. The radio jets are roughly perpendicular to the dust which is used to constrain the Doppler boosting factors of the radio jet and cores. The VLBA core flux correlates with the central Halpha+[NII] flux which might indicate that the VLBA core is dominated by an isotropic component. Twelve galaxies show nuclear optical sources. We discuss various possible origins for this emission.

Hipparcos parallaxes fix distances to individual stars in the Hyades cluster with an accuracy of 6%. We use the Hipparcos (and Tycho-2) proper motions, which have a larger relative precision than the trigonometric parallaxes, to derive ~3 times more precise distance estimates, by assuming that all members share the same space motion. The improved parallaxes as a set are statistically consistent with the Hipparcos parallaxes. The new parallaxes confirm that the Hipparcos measurements are correlated on small angular scales, consistent with the limits specified in the Hipparcos Catalogue, though with significantly smaller `amplitudes' than claimed by Narayanan & Gould. The colour-absolute magnitude diagram of the cluster based on the new paral- laxes shows a well-defined main sequence with two gaps/turn-offs. These features provide the first direct observational support of Boehm-Vitense's prediction that (the onset of) surface convection in stars affects their B-V colours. We present and discuss the theoretical HRD for an objectively defined set of 88 high-fidelity members of the cluster as well as the delta Scuti star theta^2 Tau, the giants delta^1, theta^1, epsilon, and gamma Tau, and the white dwarfs V471 Tau and HD 27483 (all of which are also members). The precision with which the new parallaxes place individual Hyades members in the Hertz- sprung-Russell diagram is limited by (systematic) uncertainties related to the transformations from observed colours and absolute magnitudes to effective temperatures and luminosities. The new parallaxes provide stringent constraints on the calibration of such transformations when combined with theoretical stellar evolutionary modelling, tailored to the chemical composition and age of the Hyades, over the large stellar mass range probed by Hipparcos.

We report the first wide-field mapping of the kinematics and stellar populations in the E3 galaxy NGC 4365. The velocity maps extend previous long-slit work. They show two independent kinematic subsystems: the central 300 x 700 pc rotates about the projected minor axis, and the main body of the galaxy, 3 x 4 kpc, rotates almost at right angles to this. The line-strength maps show that the metallicity of the stellar population decreases from a central value greater than solar, to one-half solar at a radius of 2 kpc. The decoupled core and main body of the galaxy have the same luminosity-weighted age, of ~14 Gyr, and the same elevated magnesium-to-iron ratio. The two kinematically distinct components have thus shared a common star formation history. We infer that the galaxy underwent a sequence of mergers associated with dissipative star formation that ended >12 Gyr ago. The misalignment between the photometric and kinematic axes of the main body is unambiguous evidence of triaxiality. The similarity of the stellar populations in the two components suggests that the observed kinematic structure has not changed substantially in 12 Gyr.

Hipparcos proper motions and trigonometric parallaxes allow the derivation of secular parallaxes which fix the distances to individual stars in the Hyades cluster to an accuracy of \sim 2 percent. The resulting color-absolute magnitude diagram for 92 high-fidelity single members of the cluster displays a very narrow main sequence, with two turn-offs and associated gaps. These occur at the locations where the onset of surface convection affects the B-V colors, as predicted by Boehm-Vitense thirty years ago. The new distances provide stringent constraints on the transformations of colors and absolute magnitudes to effective temperatures and luminosities, and on models of stellar interiors.

We use milli-arcsecond accuracy astrometry (proper motions and parallaxes) from Hipparcos and from radio observations to retrace the orbits of 56 runaway stars and nine compact objects with distances less than 700 pc, to identify the parent stellar group. It is possible to deduce the specific formation scenario with near certainty for two cases. (i) We find that the runaway star zeta Ophiuchi and the pulsar PSR J1932+1059 originated about 1 Myr ago in a supernova explosion in a binary in the Upper Scorpius subgroup of the Sco OB2 association. The pulsar received a kick velocity of about 350 km/s in this event, which dissociated the binary, and gave zeta Oph its large space velocity. (ii) Blaauw & Morgan and Gies & Bolton already postulated a common origin for the runaway-pair AE Aur and mu Col, possibly involving the massive highly-eccentric binary iota Ori, based on their equal and opposite velocities. We demonstrate that these three objects indeed occupied a very small volume \sim 2.5Myr ago, and show that they were ejected from the nascent Trapezium cluster. We identify the parent group for two more pulsars: both likely originate in the 50 Myr old association Per OB3, which contains the open cluster alpha Persei. At least 21 of the 56 runaway stars in our sample can be linked to the nearby associations and young open clusters. These include the classical runaways 53 Arietis (Ori OB1), xi Persei (Per OB2), and lambda Cephei (Cep OB3), and fifteen new identifications, amongst which a pair of stars running away in opposite directions from the region containing the lambda Ori cluster. Other currently nearby runaways and pulsars originated beyond 700 pc, where our knowledge of the parent groups is very incomplete.

We briefly describe the SAURON project, aimed at determining the intrinsic shape and internal dynamics of spheroids. We focus here on the ability of SAURON to identify gaseous and stellar disks and to measure their morphology and kinematics. We illustrate some of our results with complete maps of NGC3377, NGC3623, and NGC4365.

We have analysed $V$, $H$ and $J$ HST images for a sample of early- to late-type spiral galaxies, and reported elsewhere the statistical frequency of $R^{1/4}$-law and exponential bulges in our sample as a function of Hubble type, and the frequency of occurrence and structural properties of the resolved central nuclei hosted by intermediate- to late-type bulges and disks (see references in the text). Here we use these data to show that: {\it (i)} The $V-H$ color distribution of the $R^{1/4}$ bulges peaks around $<V-H> \sim 1.3$, with a sigma $\Delta (V-H) \sim 0.1$ magnitudes. Assuming a solar metallicity, these values correspond to stellar ages of $\approx 6 \pm 3$ Gyrs. In contrast, the $V-H$ color distribution of the exponential bulges peaks at $<V-H> \sim 0.9$ and has a sigma $\Delta (V-H) \sim 0.4$ mags. This likely implies significantly smaller ages and/or lower metallicities for (a significant fraction of the stars in) the exponential bulges compared to the $R^{1/4}$-law spheroids. {\it (ii)} Most of the central nuclei hosted by the exponential bulges have $V-H$ and $J-H$ colors which are compatible with relatively unobscured stellar populations. Assuming no or little dust effects, ages $\gta 1$ Gyrs are suggested for these nuclei, which in turn imply masses of about a few 10$^6$ to a few 10$^7$ M$_\odot$, i.e., sufficient to dissolve progenitor bars with masses consistent with those inferred for the exponential bulges by their luminosities. {\it (iii)} While different bulge-nucleus pairs cover a large range of $V-H$ colors, each bulge-nucleus pair has quite similar $V-H$ colors, and thus possibly similar stellar populations. (abridged)

This paper presents the properties of a family of scale-free triaxial haloes. We adduce arguments to suggest that the velocity ellipsoids of such models are aligned in conical coordinates. We provide an algorithm to find the set of conically aligned velocity second moments that support a given density against the gravity field of the halo. The case of the logarithmic ellipsoidal model -- the simplest triaxial generalisation of the familiar isothermal sphere -- is examined in detail. The velocity dispersions required to hold up the self-consistent model are analytic. The velocity distribution of the dark matter can be approximated as a triaxial Gaussian with semiaxes equal to the velocity dispersions. There are roughly twenty experiments worldwide that are searching for evidence of scarce interactions between weakly-interacting massive-particle dark matter (WIMPs) and detector nuclei. The annual modulation signal, caused by the Earth's rotation around the Sun, is a crucial discriminant between WIMP events and the background. The greatest rate is in June, the least in December. We compute the differential detection rate for energy deposited by the rare WIMP-nucleus interactions in our logarithmic ellipsoidal halo models. Triaxiality and velocity anisotropy change the total rate by up to 40 %, and have a substantial effect on the amplitude of the annual modulation signal. The overall rate is greatest, but the amplitude of the modulation is weakest, in our radially anisotropic halo models. Even the sign of the signal can be changed. Restricting attention to low energy events, the models predict that the maximum rate occurs in December, and not in June.

Milli-arcsecond astrometry provided by Hipparcos and by radio observations makes it possible to retrace the orbits of some of the nearest runaway stars and pulsars to determine their site of origin. The orbits of the runaways AE Aurigae and mu Columbae and of the eccentric binary iota Orionis intersect each other about 2.5 Myr ago in the nascent Trapezium cluster, confirming that these runaways were formed in a binary-binary encounter. The path of the runaway star zeta Ophiuchi intersects that of the nearby pulsar PSR J1932+1059, about 1 Myr ago, in the young stellar group Upper Scorpius. We propose that this neutron star is the remnant of a supernova that occurred in a binary system which also contained zeta Oph, and deduce that the pulsar received a kick velocity of about 350 km/s in the explosion. These two cases provide the first specific kinematic evidence that both mechanisms proposed for the production of runaway stars, the dynamical ejection scenario and the binary-supernova scenario, operate in nature.

We run numerical simulations of the disruption of satellite galaxies in a Galactic potential to build up the entire stellar halo, in order to investigate what the next generation of astrometric satellites will reveal by observing the halo of the Milky Way. We generate artificial DIVA, FAME and GAIA halo catalogues, in which we look for the signatures left by the accreted satellites. We develop a method based on the standard Friends-of-Friends algorithm applied to the space of integrals of motion. We find this simple method can recover about 50% of the different accretion events, when the observational uncertainties expected for GAIA are taken into account, even when the exact form of the Galactic potential is unknown. The recovery rate for DIVA and FAME is much smaller, but these missions, like GAIA, should be able to test the hierarchical formation paradigm on our Galaxy by measuring the amount of halo substructure in the form of nearby kinematically cold streams with for example, a two-point correlation function in velocity space.

The peculiar elliptical galaxy IC 1459 (M_V = -21.19, D = 16.5 Mpc) has a fast counterrotating stellar core, stellar shells and ripples, a blue nuclear point source and strong radio core emission. We present results of a detailed HST study of IC 1459, and in particular its central gas disk, aimed a constraining the central mass distribution. We obtained WFPC2 narrow-band imaging centered on the Halpha+[NII] emission lines to determine the flux distribution of the gas emission at small radii, and we obtained FOS spectra at six aperture positions along the major axis to sample the gas kinematics. We construct different dynamical models for the Halpha+[NII] and Hbeta kinematics that include a supermassive black hole, and in which the stellar mass distribution is constrained by the observed surface brightness distribution and ground-based stellar kinematics. All models are consistent with a black hole mass in the range Mbh=1-4 x 10^8 Msun, and models without a black hole are always ruled out at high confidence.

We have extended the spectral dynamics formalism introduced by Binney & Spergel, and have implemented a semi-analytic method to represent regular orbits in any potential, making full use of their regularity. We use the spectral analysis code of Carpintero & Aguilar to determine the nature of an orbit (irregular, regular, resonant, periodic) from a short-time numerical integration. If the orbit is regular, we approximate it by a truncated Fourier time series of a few tens of terms per coordinate. Switching to a description in action-angle variables, this corresponds to a reconstruction of the underlying invariant torus. We then relate the uniform distribution of a regular orbit on its torus to the non-uniform distribution in the space of observables by a simple Jacobian transformation between the two sets of coordinates. This allows us to compute, in a cell-independent way, all the physical quantities needed in the study of the orbit, including the density and in the line-of-sight velocity distribution, with much increased accuracy. The resulting flexibility in the determination of the orbital properties, and the drastic reduction of storage space for the orbit library, provide a significant improvement in the practical application of Schwarzschild's orbit superposition method for constructing galaxy models. We test and apply our method to two-dimensional orbits in elongated discs, and to the meridional motion in axisymmetric potentials, and show that for a given accuracy, the spectral dynamics formalism requires an order of magnitude fewer computations than the more traditional approaches.

It is now generally believed that galaxies were built up through gravitational amplification of primordial fluctuations and the subsequent merging of smaller precursor structures. The stars of the structures that assembled to form the Milky Way should now make up much or all of its bulge and halo, in which case one hopes to find "fossil" evidence for those precursor structures in the present distribution of halo stars. Confirmation that this process is continuing came with the discovery of the Sagittarius dwarf galaxy, which is being disrupted by the Milky Way, but direct evidence that this process provided the bulk of the Milky Way's population of old stars has so far been lacking. Here we show that about ten per cent of the metal--poor stars in the halo of the Milky Way, outside the radius of the Sun's orbit, come from a single coherent structure that was disrupted during or soon after the Galaxy's formation. This object had a highly inclined orbit about the Milky Way at a maximum distance of $\sim$ 16 kpc, and it probably resembled the Fornax and Sagittarius dwarf spheroidal galaxies.

We present and analyse HST/WFPC2 broad- and narrow-band observations of the central regions of 19 nearby radio-loud early-type galaxies. Together with two more galaxies they form a complete sample of Fanaroff & Riley Type I galaxies. We obtained V- and I-band images and narrow-band images centered on the \HalphaNII emission lines. We use archival data for six galaxies. We describe the data reduction, give isophotal fits and analyse the central emission-line gas and dust distributions. Our main conclusions are the following. Although obscuration by dust inhibits a direct determination of central cusp slopes, the data suggest most but not all galaxies have shallow cores. Dust is detected in all but two galaxies. There is a wide variety of central dust morphologies, ranging from central disks to lanes and irregular distributions. The analysis suggests the difference between disks and lanes is intrinsic and not due to different viewing angles. Central emission-line gas is detected in all galaxies. Extended low-surface-brightness emission is always associated with the dust features. In a subsequent paper we will present a detailed analysis of the relation between these central properties and the nuclear activity.

We describe an improved, practical method for constructing galaxy models that match an arbitrary set of observational constraints, without prior assumptions about the phase-space distribution function (DF). Our method is an extension of Schwarzschild's orbit superposition technique. As in Schwarzschild's original implementation, we compute a representative library of orbits in a given potential. We then project each orbit onto the space of observables, consisting of position on the sky and line-of-sight velocity, while properly taking into account seeing convolution and pixel binning. We find the combination of orbits that produces a dynamical model that best fits the observed photometry and kinematics of the galaxy. A key new element of this work is the ability to predict and match to the data the full line-of-sight velocity profile shapes. A dark component (such as a black hole and/or a dark halo) can easily be included in the models. We have tested our method, by using it to reconstruct the properties of a two-integral model built with independent software. The test model is reproduced satisfactorily, either with the regular orbits, or with the two-integral components. This paper mainly deals with the technical aspects of the method, while applications to the galaxies M32 and NGC 4342 are described elsewhere (van der Marel et al., Cretton & van den Bosch). (abridged)

Since the previous (1990) edition of this meeting enormous progress in the field of OB associations has been made. Data from X-ray satellites have greatly advanced the study of the low-mass stellar content of associations, while astrometric data from the Hipparcos satellite allow for a characterization of the higher-mass content of associations with unprecedented accuracy. We review recent work on the OB associations located within 1.5 kpc from the Sun, discuss the Hipparcos results at length, and point out directions for future research.

We investigate what the proposed ESA astrometric satellite GAIA will reveal by observing the halo of the Milky Way. Specifically, we look for halo streams which are the signatures left by the merging/accretion events experienced by a typical galaxy like the Milky Way. We run numerical simulations of the disruption of satellite galaxies in a Galactic potential to generate artificial GAIA halo catalogs. We recover the streams by searching for peaks in angular momentum space.

A comprehensive census of the stellar content of the nearby OB associations is presented, based on Hipparcos positions, proper motions, and parallaxes. Moving groups are identified by combining de Bruijne's refurbished convergent point method with the `Spaghetti method' of Hoogerwerf & Aguilar. Monte Carlo simulations are used to estimate the expected number of interloper field stars. Astrometric members are listed for 12 young stellar groups, out to a distance of ~650 pc. These are the 3 subgroups Upper Scorpius, Upper Centaurus Lupus and Lower Centaurus Crux of Sco OB2, as well as Vel OB2, Tr 10, Col 121, Per OB2, alpha Persei (Per OB3), Cas-Tau, Lac OB1, Cep OB2, and a new group designated as Cep OB6. The selection procedure corrects the list of previously known astrometric and photometric B- and A-type members, and identifies many new members, including a significant number of F stars, as well as evolved stars, e.g., the Wolf-Rayet stars gamma^2 Vel (Vel OB2) and EZ CMa (Col 121), and the classical Cepheid delta Cep in Cep OB6. In the nearest associations the later-type members include T Tauri objects and other pre-main sequence stars. Astrometric evidence for moving groups in the fields of R CrA, CMa OB1, Mon OB1, Ori OB1, Cam OB1, Cep OB3, Cep OB4, Cyg OB4, Cyg OB7, and Sct OB2, is inconclusive, due to their large distance or unfavorable kinematics. The mean distances of the well-established groups are systematically smaller than previous estimates. The mean motions display a systematic pattern, which is discussed in relation to the Gould Belt. Six of the 12 detected moving groups do not appear in the classical list of nearby OB associations. The number of unbound young stellar groups in the Solar neighbourhood may be significantly larger than thought previously.

With the steady increase of the sample size of observed microlenses towards the central regions of the Galaxy, the main source of the uncertainty in the lens mass will shift from the simple Poisson noise to the intrinsic non-uniqueness of our dynamical models of the inner Galaxy, particularly, the Galactic bar. We use a set of simple self-consistent bar models to investigate how the microlensing event rate varies as a function of axis ratio, bar angle and velocity distribution. The non-uniqueness of the velocity distribution of the bar model adds a significant uncertainty (by about a factor of 1.5) to any prediction of the lens mass. Kinematic data and self-consistent models are critical to lift the non-uniqueness. We discuss the implications of these results for the interpretation of microlensing observations of the Galactic bulge. In particular we show that Freeman bar models scaled to the mass of the Galactic bulge/bar imply a typical lens mass of around $0.8M_\odot$, a factor of 3-5 times larger than the value from other models.

We present a method for measuring small deviations from axisymmetry of the potential of a filled gas disk. The method is based on a higher order harmonic expansion of the full velocity field of the disk. This expansion is made by first fitting a tilted-ring model to the velocity field of the gas disk and subsequently expanding the velocity field along each ring into its harmonic terms. We use epicycle theory to derive equations for the harmonic terms in a distorted potential. The phase of each component of the distortion can vary with radius. We show that if the potential has a distortion of harmonic number m, the velocity field as seen on the sky exhibits an m-1 and m+1 distortion. As is to be expected, the effects of a global elongation of the halo are similar to an m=2 spiral arm. The main difference is that the phase of the spiral arm can vary with radius. Our method allows a measurement of epsilon_pot sin(2 phi_2), where epsilon_pot is the elongation of the potential and phi_2 is one of the viewing angles. Using \hi data, one can probe the potential at radii beyond the stellar disk, into the regime where dark matter is thought to be the dominant dynamical component. The method is applied the spiral galaxies NGC 2403 and NGC 3198 and the harmonic terms are measured up to ninth order. We find epsilon_pot sin(2 phi_2) to be 0.064 +/- 0.003 for NGC 2403 and 0.019 +/- 0.003 for NGC 3198. More galaxies should be examined to separate viewing angle from elongation in a statistical way.

We present a new procedure for the identification of moving groups. It is a combination of two independent member selection methods. One is a modern implementation of a classical convergent point method for proper motion data. The other objectively identifies moving groups in velocity space using proper motions and parallaxes. We briefly describe both methods, and illustrate their power by applying them to Hipparcos measurements of a field containing the Upper Scorpius subgroup of the Scorpio-Centaurus association (Sco OB2). We show how our membership selection procedure not only improves the list of previously known B and A-type members, but also identifies many new members, including a significant number of F stars. We apply our procedure to other nearby OB associations elsewhere in these proceedings (Hoogerwerf et al.; de Zeeuw et al.).

We present the first results of a comprehensive census of the stellar content of the nearby OB associations based on Hipparcos positions, proper motions and parallaxes for 12842 candidate member stars distributed over 21 fields on the sky. We use a new method to identify moving groups in these fields (see de Bruijne et al., these proceedings). Previously, astrometric membership in nearly all the nearby OB associations was known only for stars with spectral types earlier than B5. The Hipparcos measurements now allow us to identify members down to late F. This census provides a firm basis for studies of galactic and extragalactic star forming regions.

We have developed a new procedure to identify moving groups, and have applied it to carry out a census of the nearby OB associations based on Hipparcos parallaxes and proper motions (see de Bruijne et al. and de Zeeuw et al. elsewhere in this volume). Here we present three illustrative cases. For $\alpha$ Persei our method allows us to refine the bright end of the membership list, while for Collinder 121 the Hipparcos data change the whole appearance of the association. Finally, we report the discovery of a new association in the field of Cepheus OB2.

Axisymmetric dynamical models are constructed for the E3 galaxy M32 to interpret high spatial resolution stellar kinematical HST data. Models are studied with two-integral phase-space distribution functions, and with fully general three- integral distribution functions. The latter are built using a new extension of Schwarzschild's orbit superposition approach. Models are constructed for inclinations of 90 and 55 degrees. No model without a nuclear dark object can fit the combined ground-based and HST data, independent of the dynamical structure of M32. Models with a nuclear dark object of 3.4 x 10^6 solar masses do provide an excellent fit. The inclined models provide the best fit, but the inferred dark mass does not depend sensitively on the assumed inclination. The models that best fit the data are not two-integral models, but like two-integral models they are azimuthally anisotropic. An extended dark object can fit the data only if its half-mass radius is r_h < 0.08 arcsec (=0.26 pc), implying a central dark matter density exceeding 1 x 10^8 solar masses / pc^3. This density is high enough to rule out most plausible alternatives to a massive black hole. The dynamically inferred dark mass is identical to that suggested by existing adiabatic black hole growth models for HST photometry of M32. The low activity of M32 implies either that only a very small fraction of the gas that is shed by evolving stars is accreted onto the black hole, or that accretion proceeds at very low efficiency, e.g. in an advection-dominated mode. (shortened version)

We have obtained spectra through small apertures centered on the nuclear region and major axis of M32, with the HST/FOS. A detailed analysis and reduction of the data is presented, including: (i) new calibrations and modeling of the FOS aperture sizes, PSF and line-spread-functions; (ii) determination of the aperture positioning for each observation; and (iii) accurate wavelength calibration, template matching, and kinematical analysis of the spectra. This yields measurements of the stellar rotation velocities and velocity dispersions with five times higher spatial resolution than the best available ground-based data. The inferred velocities provide the highest angular resolution stellar kinematical data obtained to date for any stellar system. The observations show a steeper rotation curve and higher central velocity dispersion than ground-based data. The rotation velocity is observed to be 30 km/s at 0.1" from the nucleus; roughly twice the value measured from the ground at this distance. The nuclear dispersion measured through the smallest FOS aperture (0.068" square) is 156 +/- 10 km/s. The average of four independent dispersion measurements at various positions inside the central 0.1" is 126 km/s, with a RMS scatter of 21 km/s. The nuclear dispersion measured from the ground is only 85-95 km/s, whereas the dispersion outside the central arcsec is only 45-55 km/s. These results significantly strengthen previous arguments for the presence of a massive nuclear black hole in M32. Detailed dynamical models are presented in a series of companion papers.

Massive black holes are thought to reside at the centres of many galaxies, where they power quasars and active galactic nuclei. But most galaxies are quiescent, indicating that any central massive black hole present will be starved of fuel and therefore detectable only through its gravitational influence on the motions of the surrounding stars. M32 is a nearby, quiescent elliptical galaxy in which the presence of a black hole has been suspected; however, the limited resolution of the observational data and the restricted classes of models used to interpret this data have made it difficult to rule out alternative explanations, such as models with an anisotropic stellar velocity distribution and no dark mass or models with a central concentration of dark objects (for example, stellar remnants or brown dwarfs). Here we present high-resolution optical HST spectra of M32, which show that the stellar velocities near the centre of this galaxy exceed those inferred from previous ground-based observations. We use a range of general dynamical models to determine a central dark mass concentration of (3.4 +/- 1.6) x 10^6 solar masses, contained within a region only 0.3 pc across. This leaves a massive black hole as the most plausible explanation of the data, thereby strengthening the view that such black holes exist even in quiescent galaxies.

The Jeans equations give the second moments or stresses required to support a stellar population against the gravity field. A general solution of the Jeans equations for arbitrary axisymmetric scale-free densities in flattened scale-free potentials is given. A two-parameter subset of the solution for the second moments for the self-consistent density of the power-law models, which have exactly spheroidal equipotentials, is examined in detail. In the spherical limit, the potential of these models reduces to that of the singular power-law spheres. We build the physical three-integral distribution functions that correspond to the flattened stellar components. Next, we attack the problem of finding distribution functions associated with the Jeans solutions in flattened scale-free potentials. The third or partial integral introduced by de Zeeuw, Evans and Schwarzschild for Binney's model is generalised to thin and near-thin orbits moving in arbitrary axisymmetric scale-free potentials. The partial integral is a modification of the total angular momentum. For the self-consistent power-law models, we show how this enables the construction of simple three-integral distribution functions. The connexion between these approximate distribution functions and the Jeans solutions is discussed in some detail.

Recently, observations with the Hubble Space Telescope have revealed small stellar discs embedded in the nuclei of a number of ellipticals and S0s. In this paper we construct two-integral axisymmetric models for such systems. We calculate the even part of the phase-space distribution function, and specify the odd part by means of a simple parameterization. We investigate the photometric as well as the kinematic signatures of nuclear discs, including their velocity profiles (VPs), and study the influence of seeing convolution. The rotation curve of a nuclear disc gives an excellent measure of the central mass-to-light ratio whenever the VPs clearly reveal the narrow, rapidly rotating component associated with the nuclear disc. Steep cusps and seeing convolution both result in central VPs that are dominated by the bulge light, and these VPs barely show the presence of the nuclear disc, impeding measurements of the central rotation velocities of the disc stars. However, if a massive BH is present, the disc component of the VP can be seen in the wing of the bulge part, and measurements of its mean rotation provide a clear signature of the presence of the BH. This signature is insensitive to the uncertainties in the velocity anisotropy, which often lead to ambiguity in the interpretation of a central rise in velocity dispersion as due to a central BH.

This paper presents two families of phase-space distribution functions (DFs) that generate scale-free spheroidal mass densities in scale-free spherical potentials. The `case I' DFs are anisotropic generalizations of the flattened f(E,L_z) model, which they include as a special case. The `case II' DFs generate flattened constant-anisotropy models, in which the constant ratio of rms tangential to radial motion is characterized by Binney's parameter beta. The models can describe the outer parts of galaxies and the density cusp structure near a central black hole, but also provide general insight into the dynamical properties of flattened systems. The dependence of the intrinsic and projected properties on the model parameters and the inclination is described. The observed ratio of the rms projected line-of-sight velocities on the projected major and minor axes of elliptical galaxies is best fit by the case II models with beta > 0. These models also predict non-Gaussian velocity profile shapes consistent with existing observations. The distribution functions are used to model the galaxies NGC 2434 (E1) and NGC 3706 (E4), for which stellar kinematical measurements out to two effective radii indicate the presence of dark halos (Carollo et al.). The velocity profile shapes of both galaxies can be well fit by radially anisotropic case II models with a spherical logarithmic potential. This contrasts with the f(E,L_z) models studied previously, which require flattened dark halos to fit the data.

A general solution of the Jeans equations for oblate scale-free logarithmic potentials is given. This provides all possible second velocity moments that can hold up a stellar population of flattened scale-free density against the gravity field. A two-parameter subset of second moments for the self-consistent density of Binney's model is examined in detail. These solutions have the desirable property that the observable dispersions in the radial and proper motions can be given explicitly. In the spherical limit, the potential of these models reduces to that of the singular isothermal sphere. The Jeans solutions for scale-free densities of arbitrary flattening that can correspond to physical three-integral distribution functions are identified. The problem of finding distribution functions associated with the Jeans solutions in flattened scale-free logarithmic potentials is then investigated for Binney's model. An approximate solution of the collisionless Boltzmann equation is found which provides a third (partial) integral of good accuracy for thin and near-thin tube orbits. It is a modification of the total angular momentum. This enables the construction of many simple three-integral distribution functions. The kinematic properties of these approximate DFs are shown to agree with a subset of the Jeans solutions --- which are thereby confirmed as good approximations to physical solutions.

Luminous galaxies and their dark halos are likely to have triaxial shapes. The construction of distribution functions for triaxial systems is a hard problem, due to the existence of only one exact integral of motion, the energy $E$. Even so, progress has been made by use of analytic methods for special potentials, numerical combination of orbit densities (Schwarzschild's method), and N-body simulations. Despite these efforts, we are a long way from understanding -- let alone constructing -- the full variety of models needed to interpret the photometric and kinematic data on galaxies that is now available. Some of the main unsolved problems include the range of intrinsic shapes and figure rotation speeds for which stable dynamical equilibria can be built, and the dynamical role of cusps and central black holes.

Recent work on the construction of spherical, Axisymmetric and triaxial dynamical models for elliptical galaxies is reviewed briefly, including their role in providing evidence for dark halos and central black holes. The different orbital structures and shapes of low-mass and giant elliptical galaxies provide essential constraints on scenarios of galaxy formation.

A simple numerical scheme is presented for the construction of three-integral phase-space distribution functions for oblate galaxy models with a gravitational potential of St\"{a}ckel form, and an arbitrary axisymmetric luminous density distribution. The intrinsic velocity moments can be obtained simultaneously with little extra effort. The distribution of the inner and outer turning points of the short-axis tube orbits that are populated can be specified freely, and is chosen in advance. The entire distribution function is then derived from the density by an iterative scheme that starts from the explicitly known distribution function of the thin-orbit (maximum streaming) model, in which only the tubes with equal inner and outer turning points are populated. The versatility and limitations of this scheme are illustrated by the construction of a number of self-consistent three-integral flattened isochrone models of Kuzmin--Kutuzov type, and by investigation of special cases where the scheme is tractable analytically. This includes the behaviour of the distribution functions in the outer regions of the models. The scheme converges rapidly for models containing orbits with ratios of the outer to inner turning point as large as ten, and is particularly suited for the construction of tangentially anisotropic flattened models, self-consistent as well as non-consistent. The algorithm simplifies in the disk and spherical limit, and can be generalized to triaxial models.

A simple algorithm is presented for the calculation of the projected line-of-sight velocity profiles (VPs) of non-rotating anisotropic spherical stellar dynamical models with a phase-space distribution function of the Osipkov-Merritt type. The velocity distribution in these models is isotropic inside the anisotropy radius r_a and becomes increasingly radially biased at larger radii. VP shape parameters are presented for a family of models in which the luminous mass density has an r^{-gamma} power-law cusp at small radii and an r^{-4} power-law fall-off at large radii. Self-consistent models and models in which the luminous matter is embedded in a dark halo are discussed. The effects of changes in the cusp slope gamma and in the anisotropy radius r_a are documented, and the area in the (gamma,r_a)-plane that contains physical models is delineated. The shapes of the VPs of the models show a considerable (and observable) variation with projected galactocentric radius. These models will be useful for interpreting the data on the VP shapes of elliptical galaxies that are now becoming available.

The contour integral method of Hunter & Qian is applied to axisymmetric galaxy models in which the distribution function (DF) is of the form f=f(E,L_z), where E and L_z are the classical integrals of motion in an axisymmetric potential. A practical way to construct the unique even part of the DF for such systems is presented. It is applied to models, both oblate and prolate, in which the mass density is stratified on similar concentric spheroids. The spheroids with scale-free densities are discussed in detail. These provide useful approximations to the behaviour of more realistic models in the limit of small and large radii. The self-consistent case is treated, as well as the case in which there are additional contributions to the potential from a central black hole or dark halo. The range of axis ratios and density profile slopes is determined for which spheroidal power--law cusps with a central black hole have a physical two--integral DF. More generally, the f(E,L_z) DFs are discussed for a set of spheroidal `(alpha,beta)-models', characterized by a power-law density cusp with slope alpha at small radii, and a power-law density fall-off with slope [alpha + 2 beta] at large radii. As an application, the DF is constructed for the (alpha,beta) model with a 1.8 x 10^6 solar mass black hole used by van der Marel et al. to interpret their high spatial resolution spectroscopic data for M32. The results confirm that the model fits the observed line-of-sight velocity profiles remarkably well. The model is used to calculate predictions for future spectroscopic observations with the HST.

Walraven photometry of established and probable members of the Orion OB1 association is presented. Effective temperature, surface gravity, luminosity and mass are derived for all stars, using atmosphere models by Kurucz (1979). Absolute magnitudes are calculated using the Straizys and Kuriliene (1981) tables. Distance moduli and visual extinctions are determined. A comparison of the visual extinctions to IRAS $100\um$ data shows that the near edge of the Orion A and B clouds lies at a distance of $\sim 320\pc$, while the far edge is at $\sim 500\pc$. A method for deriving the ages of the subgroups by comparing theoretical isochrones to the observations in the log g, log T plane is presented. The derived ages suggest, contrary to earlier studies, that subgroup 1b is younger than 1c, which can possibly be explained by past geometries of the system of stars and gas. The initial mass function for Orion OB1 is derived with the aid of the Kolmogorov-Smirnoff test. Through extensive simulations, we show that it is very difficult to derive accurately the IMF from the available data. To within somewhat weak limits the IMF is found to be of the form $\xi (\log M)=AM^{-1.7\pm 0.2}$ for all subgroups. The energy output of the subgroups in the form of stellar winds and supernovae is calculated and compared to the observed size and expansion velocity of the Orion-Eridanus bubble. It is shown that the energy output of the association can account for the morphology and kinematics of the ISM.