NGC 4494 is one of several intermediate-luminosity elliptical galaxies inferred to have an unusually diffuse dark matter halo. We use the chi^2-made-to-measure particle code NMAGIC to construct axisymmetric models of NGC 4494 from photometric and various kinematic data. The extended kinematics include light spectra in multiple slitlets out to 3.5 R_e, and hundreds of planetary nebulae velocities out to ~7 R_e, thus allowing us to probe the dark matter content and orbital structure in the halo. We use Monte Carlo simulations to estimate confidence boundaries for the halo parameters, given our data and modelling set-up. We find that the true potential of the dark matter halo is recovered within Delta G (merit function)<26 (Delta chi^2<59) at 70% confidence level (C.L.), and within Delta G<32 (Delta chi^2<70) at 90% C.L.. These numbers are much larger than the usually assumed Delta chi^2=2.3 (4.6) for 70% (90%) C.L. for two free parameters, perhaps case-dependent, but calling into question the general validity of the standard assumptions used for halo and black hole mass determinations. The best-fitting models for NGC 4494 have a dark matter fraction of about 0.6\pm0.1 at 5R_e (70% C.L.), and are embedded in a dark matter halo with circular velocity ~200 km/s. The total circular velocity curve (CVC) is approximately flat at v_c=220 km/s outside ~0.5R_e. The orbital anisotropy of the stars is moderately radial. These results are independent of the assumed inclination of the galaxy, and edge-on models are preferred. Comparing with the halos of NGC 3379 and NGC 4697, whose velocity dispersion profiles also decrease rapidly from the center outwards, the outer CVCs and dark matter halos are quite similar. NGC 4494 shows a particularly high dark matter fraction inside ~3R_e, and a strong concentration of baryons in the center.

The parallel code NMAGIC is an implementation of a particle-based method to create made-to-measure models in agreement with observations of galaxies. It works by slowly correcting the particle weights of an evolving N-body system, until a satisfactory compromise is achieved between the goodness of the fit to a given set of observational data, and some degree of smoothness (regularization) of the underlying particle model. We briefly describe the method together with a new regularization scheme in phase-space, which improves recovering the correct orbit structure in the models. We also mention some practical applications showing the power of the technique in investigating the dynamics of galaxies.

Made-to-measure methods such as the parallel code NMAGIC are powerful tools to build galaxy models reproducing observational data. They work by adapting the particle weights in an N-body system until the target observables are well matched. Here we introduce a moving prior regularization (MPR) method for such particle models. It is based on determining from the particles a distribution of priors in phase-space, which are updated in parallel with the weight adaptation. This method allows one to construct smooth models from noisy data without erasing global phase-space gradients. We first apply MPR to a spherical system for which the distribution function can in theory be uniquely recovered from idealized data. We show that NMAGIC with MPR indeed converges to the true solution with very good accuracy, independent of the initial particle model. Compared to the standard weight entropy regularization, biases in the anisotropy structure are removed and local fluctuations in the intrinsic distribution function are reduced. We then investigate how the uncertainties in the inferred dynamical structure increase with less complete and noisier kinematic data, and how the dependence on the initial particle model also increases. Finally, we apply the MPR technique to the two intermediate-luminosity elliptical galaxies NGC 4697 and NGC 3379, obtaining smoother dynamical models in luminous and dark matter potentials.

Following the seminal result of An & Evans, known as the central density slope-anisotropy theorem, successive investigations unexpectedly revealed that the density slope-anisotropy inequality holds not only at the center, but at all radii in a very large class of spherical systems whenever the phase-space distribution function is positive. In this paper we derive a criterion that holds for all spherical systems in which the augmented density is a separable function of radius and potential: this new finding allows to unify all the previous results in a very elegant way, and opens the way for more general investigations. As a first application, we prove that the global density slope-anisotropy inequality is also satisfied by all the explored additional families of multi-component stellar systems. The present results, and the absence of known counter-examples, lead us to conjecture that the global density slope-anisotropy inequality could actually be a universal property of spherical systems with positive distribution function.

Starting from the central density slope-anisotropy theorem of An and Evans (2006), recent investigations have shown that the involved density slope-anisotropy inequality holds not only at the center, but at all radii (i.e. globally) in a very large class of spherical systems with positive phase-space distribution function. Here we present some additional analytical cases that further extend the validity of the global density slope-anisotropy inequality. These new results, several numerical evidences, and the absence of known counter-examples, lead us to conjecture that the global density slope-anisotropy inequality could actually be a universal property of spherical systems with positive distribution function.

We present new optical long-slit data along 6 position angles of the bulge region of M31. We derive accurate stellar and gas kinematics reaching 5 arcmin from the center, where the disk light contribution is always less than 30%, and out to 8 arcmin along the major axis, where the disk makes 55% of the total light. We show that the velocity dispersions of McElroy (1983) are severely underestimated (by up to 50 km/s) and previous dynamical models have underestimated the stellar mass of M31's bulge by a factor 2. Moreover, the light-weighted velocity dispersion of the galaxy grows to 166 km/s, thus reducing the discrepancy between the predicted and measured mass of the black hole at the center of M31. The kinematic position angle varies with distance, pointing to triaxiality. We detect gas counterrotation near the bulge minor axis. We measure eight emission-corrected Lick indices. They are approximately constant on circles. We derive the age, metallicity and alpha-element overabundance profiles. Except for the region in the inner arcsecs of the galaxy, the bulge of M31 is uniformly old (>12 Gyr, with many best-fit ages at the model grid limit of 15 Gyr), slightly alpha-elements overabundant ([alpha/Fe]~0.2) and at solar metallicity, in agreement with studies of the resolved stellar components. The predicted u-g, g-r and r-i Sloan color profiles match reasonably well the dust-corrected observations. The stellar populations have approximately radially constant mass-to-light ratios (M/L_R ~ 4-4.5 for a Kroupa IMF), in agreement with stellar dynamical estimates based on our new velocity dispersions. In the inner arcsecs the luminosity-weighted age drops to 4-8 Gyr, while the metallicity increases to above 3 times the solar value.

General criteria to check the positivity of the distribution function (phase-space consistency) of stellar systems of assigned density and anisotropy profile are useful starting points in Jeans-based modeling. Here we substantially extend previous results, and we present the inversion formula and the analytical necessary and sufficient conditions for phase-space consistency of the family of multi-component Cuddeford spherical systems: the distribution function of each density component of these systems is defined as the sum of an arbitrary number of Cuddeford distribution functions with arbitrary values of the anisotropy radius, but identical angular momentum exponent. The radial trend of anisotropy that can be realized by these models is therefore very general. As a surprising by-product of our study, we found that the ``central cusp-anisotropy theorem'' (a necessary condition for consistency relating the values of the central density slope and of the anisotropy parameter) holds not only at the center, but at all radii in consistent multi-component generalized Cuddeford systems. This last result suggests that the so--called mass--anisotropy degeneracy could be less severe than what is sometimes feared.

It is well known that the density and anisotropy profile in the inner regions of a stellar system with positive phase-space distribution function are not fully independent. Here we study the interplay between density profile and orbital anisotropy at large radii in physically admissible (consistent) stellar systems. The analysis is carried out by using two-component (n-gamma,gamma_1) spherical self-consistent galaxy models, in which one density distribution follows a generalized gamma profile with external logarithmic slope n, and the other a standard gamma_1 profile (with external slope 4). The necessary and sufficient conditions for phase-space consistency are determined analytically, also in presence of a dominant massive central black hole, and the analytical phase-space distribution function of (n-gamma,1) models, and of n-gamma models with a central black hole, is derived for gamma=0,1,2. It is found that the density slope in the external regions of a stellar system can play an important role in determining the amount of admissible anisotropy: in particular, for fixed density slopes in the central regions, systems with a steeper external density profile can support more radial anisotropy than externally flatter models. This is quantified by an inequality formally identical to the ``cusp slope-central anisotropy'' theorem (An & Evans 2006), relating at all radii (and not just at the center) the density logarithmic slope and the anisotropy indicator in all Osipkov-Merritt systems [Abridged].

Dynamical properties of two-component galaxy models whose stellar density distribution is described by a gamma-model while the total density distribution has a pure r^(-2) profile, are presented. The orbital structure of the stellar component is described by Osipkov-Merritt anisotropy, while the dark matter halo is isotropic. After a description of minimum halo models, the positivity of the phase-space density (the model consistency) is investigated, and necessary and sufficient conditions for consistency are obtained analytically as a function of the stellar inner density slope gamma and anisotropy radius. The explicit phase-space distribution function is recovered for integer values of gamma, and it is shown that while models with gamma>4/17 are consistent when the anisotropy radius is larger than a critical value (dependent on gamma), the gamma=0 models are unphysical even in the fully isotropic case. The Jeans equations for the stellar component are then solved analytically; in addition, the projected velocity dispersion at the center and at large radii are also obtained analytically for generic values of the anisotropy radius, and it is found that they are given by remarkably simple expressions. The presented models, even though highly idealized, can be useful as starting point for more advanced modeling of the mass distribution of elliptical galaxies in studies combining stellar dynamics and gravitational lensing.