In this Letter, we investigate the claim of Valluri et al. (2003), namely that the use of Schwarzschild dynamical (orbit-based) models leads to an indeterminacy regarding the estimation of the free parameters like the central black hole mass and the stellar mass-to-light ratio of the galaxy under study. We examine this issue with semi-analytic two-integral models, which are not affected by the intrinsic degeneracy of three-integral systems. We however confirm the Valluri et al. result, and observe the so-called widening of the chi2 contours as the orbit library is expanded. We also show that, although two-dimensional data coverage help in constraining the orbital structure of the modelled galaxy, it does not in principle solve the indeterminacy issue, which mostly originates from the discretization of such an inverse problem. We show that adding regularization constraints stabilizes the confidence level contours on which the estimation of black hole mass and stellar mass-to-light ratio are based. We therefore propose to systematically use regularization as a tool to prevent the solution to depend on the orbit library. Regularization, however, introduces an unavoidable bias on the derived solutions. We hope that the present Letter will trigger some more research directed at a better understanding of the issues addressed here.

We present a unique set of nested stellar kinematical maps of NGC 3377 obtained with the integral-field spectrographs OASIS and SAURON. We then construct general axisymmetric dynamical models for this galaxy, based on the Schwarzschild numerical orbit superposition technique applied to these complementary measurements. We show how these two datasets constrain the mass of the central massive object and the overall mass-to-light ratio of the galaxy by probing the inner and outer regions respectively. The simultaneous use of both datasets leads us to confirm the presence of a massive black hole with a mass of $M_{BH} = 7_{-5}^{+4} 10^{7} M_\sun$ (99.7% confidence level), with a best-fit stellar mass-to-light ratio $\Upsilon_I = 2.1 \pm 0.2$ (for an assumed edge-on inclination).

We test the formation of low-luminosity elliptical galaxies through collisionless mergers of unequal-mass disk galaxies. The kinematic properties of a small survey of simulated merger end-products with initial disk mass ratios of 3:1 is compared to a sample of seven low-luminosity galaxies observed by Rix et al. that were chosen photometrically to be "ellipticals". In this paper, we go beyond a comparison in terms of global properties (using one number to characterize a model or a galaxy e.g. <a4>, ellipticity at some fixed radius or central velocity dispersion) and examine the detailed kinematics as function of galactocentric distance. The merger remnants are ``observed'' through a slit along the major and minor axis, using a pixel binning and slit width similar to the one used during the spectroscopic observations. We compare the rotational support (v/sigma) of the merger remnants, the normalized rotation on the minor axis and the major axis Gauss-Hermite moments h3, to that of the observed galaxies. The N-body remnants are very flattened when viewed edge-on (<epsilon> = 0.6) and should be inclined before making a fair comparison with the Rix et al. data set (which has <epsilon> = 0.3). When the merger remnants are appropriately inclined, their v/sigma profiles rise slower than the observed ones. Note that even when the remnants are viewed edge-on, the v/sigma profiles do not match the observations. The detailed comparison of the observations with purely collisionless 3:1-merger remnants suggests that this kind of dissipationless merger is not likely to be the dominant formation channel for low-luminosity elliptical galaxies. (abridged)

We present direct observational constraints on the orbital distribution of the stars in the giant elliptical NGC 2320. Long-slit spectra along multiple position angles are used to derive the stellar line-of-sight velocity distribution within one effective radius. In addition, the rotation curve and dispersion profile of an ionized gas disk are measured from the [OIII] emission lines. After correcting for the asymmetric drift, we derive the circular velocity of the gas, which provides an independent constraint on the gravitational potential. To interpret the stellar motions, we build axisymmetric three-integral dynamical models based on an extension of the Schwarzschild orbit- superposition technique. We consider two families of gravitational potential, one in which the mass follows the light (i.e. no dark matter) and one with a logarithmic gravitational potential. Using chi^2-statistics, we compare our models to both the stellar and gas data to constrain the value of the V-band mass-to-light ratio Upsilon-V. We find Upsilon-V = 15.0 \pm 0.6 h75 for the mass-follows-light models and Upsilon-V = 17.0 \pm 0.7 h75 for the logarithmic models. For the latter, Upsilon-V is defined within a sphere of 15'' radius. Models with radially constant Upsilon-V and logarithmic models with dark matter provide comparably good fits to the data and possess similar dynamical structure. Across the full range of Upsilon-V permitted by the observational constraints, the models are radially anisotropic in the equatorial plane over the radial range of our kinematical data (1'' < r < 40''). Along the true minor axis, they are more nearly isotropic. (abridged)

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)

We have constructed axisymmetric dynamical models of the edge-on S0 galaxy NGC 4342: simple two-integral Jeans models as well as fully general, three-integral models using a modified version of Schwarzschild's orbit superposition technique. The two-integral models suggest a black hole (BH) of 3 or $6\times 10^8 M_\odot$, depending on the data set. The three-integral models can fit all ground-based and HST data simultaneously, but only when a central BH is included. Models without BH are ruled out at better than 99.73% confidence level. We determine a BH mass of $3.0^{+1.7}_{-1.0} \times 10^8 M_\odot$. This corresponds to 2.6% of the bulge mass, making NGC 4342 one of the galaxies with the highest BH mass to bulge mass ratio currently known.

We present axisymmetric dynamical models of the edge-on S0 galaxy NGC 4342. A combination of observations from the ground and with the Hubble Space Telescope (HST) has shown that NGC 4342 rotates rapidly and has a strong central increase in velocity dispersion. We construct simple two-integral Jeans models as well as fully general, three-integral models. The latter are built using a modified version of Schwarzschild's orbit-superposition technique. The two-integral Jeans models suggest a black hole (BH) mass between 3 and 6x10^8 Msun, depending on the data set used to constrain the model, but they fail to fit the details of the observed kinematics. The three-integral models can fit all ground-based and HST data simultaneously, but only when a central BH is included. Models without BH are ruled out to a confidence level better than 99.73 per cent. We determine a BH mass of 3x10^8 Msun. This corresponds to 2.6 per cent of the total mass of the bulge, making NGC 4342 one of the galaxies with the highest BH mass to bulge mass ratio currently known. The models that best fit the data do not have a two-integral phase-space distribution function. They have rather complex dynamical structures: the velocity anisotropies are strong functions of radius reflecting the multi-component structure of this galaxy. The best fit model without BH tries to fit the high central velocity dispersion by placing stars on radial orbits. The measured rotation velocities, however, restrict the amount of radial anisotropy such that the central velocity dispersion measured with the HST can only be fit when a massive BH is included in the models.

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 describe a powerful technique to model and interpret the stellar line-of-sight velocity profiles of galaxies. Following Schwarzschild's modeling approach, a representative library of orbits is calculated in a given potential; then the non-negative superposition of these orbits is determined to fit best a given set of observational constraints. Our implementation incorporates several new features: (i) we calculate velocity profiles and represent them by a Gauss-Hermite series. This allows us to constrain the orbital anisotropy in the fit. (ii) we take into account the error on each observational constraint to obtain an objective chi2 measure for the quality-of-fit. Only projected, observable quantities are included in the fit, and aperture binning and seeing convolution of the data are properly taken into account. This scheme is valid for any geometry, but here we focus on spherical geometry and the issue of dark halos around elliptical galaxies. We model radially extended velocity profiles of the E0 galaxy NGC 2434, and find that constant M/L models are clearly ruled out, regardless of the orbital anisotropy. To study how much dark matter is needed, we considered a sequence of cosmologically motivated `star+halo' potentials, which are specified by the stellar mass-to-light ratio Gamma and the characteristic halo velocity, V_200 (from Navarro et al. 1996). The star+halo models provide an excellent fit to the data, with Gamma=3.35+-0.25 (in B-band solar units) and V_200=450+-100km/s. The best-fitting potential has a circular velocity Vc that is constant (at ~300km/s) to within 10% between 0.2--3 effective radii. In NGC 2434 roughly half of the mass within an effective radius appears to be dark.