Dark matter content and internal dynamics of NGC 4697: nmagic particle models from slit data and planetary nebula velocities

De Lorenzi, Flavio ; Gerhard, Ortwin ; Saglia, Roberto P. ; Sambhus, Niranjan ; Debattista, Victor P. ; Pannella, Maurilio ; Méndez, Roberto H.

In: Monthly Notices of the Royal Astronomical Society, 2008, vol. 385, no. 4, p. 1729-1748

Add to personal list
    Summary
    We present a dynamical study of NGC 4697, an almost edge-on, intermediate-luminosity, E4 elliptical galaxy, combining new surface brightness photometry, new as well as published long-slit absorption-line kinematic data, and published planetary nebula (PN) velocity data. The combined kinematic data set extends out to ≃5 arcmin (≃4.5Re) and allows us to probe the galaxy's outer halo. For the first time, we model such a data set with the new and flexible χ2-made-to-measure particle code nmagic. We extend nmagic to include seeing effects, introduce an efficient scheme to estimate the mass-to-light ratio, and incorporate a maximum-likelihood technique to account for discrete velocity measurements. For modelling the PN kinematics, we use line-of-sight velocities and velocity dispersions computed on two different spatial grids, and we also use the individual velocity measurements with the likelihood method, in order to make sure that our results are not biased by the way we treat the PN measurements. We generate axisymmetric self-consistent models as well models including various dark matter haloes. These models fit all the mean velocity and velocity dispersion data with χ2/N < 1, both in the case with only luminous matter and in potentials including quite massive haloes. The likelihood analysis together with the velocity histograms suggest that models with low-density haloes such that the circular velocity vc≲ 200 km s−1 at 5Re are not consistent with the data. A range of massive haloes with vc≃ 250 km s−1 at 5Re fit the PN data best. To derive stronger results would require PN velocities at even larger radii. The best-fitting models are slightly radially anisotropic; the anisotropy parameter β≃ 0.3 at the centre, increasing to β≃ 0.5 at radii ≳2Re