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

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

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

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

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 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, but these are not very helpful, as they contain six independent moments. By assuming that the potential is triaxial and of separable Staeckel form, the mixed moments vanish in confocal ellipsoidal coordinates. 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. They were first derived by Lynden-Bell in 1960, but have resisted solution by standard methods. Here we present the general solution by superposition of singular solutions.

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

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

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

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

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

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

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

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.

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.

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.

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.

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.