Issue 
A&A
Volume 541, May 2012



Article Number  A138  
Number of page(s)  8  
Section  Extragalactic astronomy  
DOI  https://doi.org/10.1051/00046361/201118704  
Published online  15 May 2012 
Isotropic dynamical models of NGC 4472 (M 49)
Astronomical Observatory, Volgina 7, 11060 Belgrade, Serbia
email: srdjan@aob.bg.ac.rs
Received: 21 December 2011
Accepted: 19 March 2012
Context. The existence of dark matter in some earlytype galaxies is wellestablished. We address one such case and study dynamical models of NGC 4472 by means of using both Newtonian and MOdified Newtonian Dynamics (MOND) approaches with particular emphasis on the latter.
Aims. We use the measurements of the radial velocities of globular clusters in NGC 4472 to test the predictions of dynamical models to fit a high and approximately constant value the velocity dispersion with radius for this galaxy.
Methods. The globular clusters of NGC 4472 are used as a tracer of the galaxy’s gravitational potential. We calculate the Jeans equation for both the Newtonian (massfollowslight and darkmatter models) and the MOND approaches assuming spherical symmetry and isotropic orbits.
Results. We solve the isotropic Jeans equation for a spherical approximation and find, assuming a rotational velocity of 50 kms^{1}, that the isotropic Newtonian massfollowslight models without dark matter cannot provide successful fits at radii larger than ~2 R_{e} and that in the Newtonian approach dark matter is needed. Of the four different MOND models that we tested, two models can provide a successful fit to the velocity dispersion throughout the whole galaxy without assuming any anisotropy (again assuming a rotational velocity of 50 kms^{1}), thus showing that in the case of NGC 4472 MOND can describe the dynamics of a massive earlytype galaxy with high values of velocity dispersion.
Key words: gravitation / galaxies: elliptical and lenticular, cD / galaxies: kinematics and dynamics / galaxies: halos / galaxies: individual: NGC 4472 (M 49)
© ESO, 2012
1. Introduction
Although some earlytype galaxies appear to lack dark matter (see, for example, the case of our nearest elliptical NGC 5128, see Samurović 2010) for some giant ellipticals all evidence inevitably suggests that they are found in massive dark haloes. One such example is a bright elliptical NGC 4472 (M 49) in the Virgo cluster. Numerous detailed studies of this galaxy (for example, Trinchieri et al. 1986; Loewenstein 1992; Saglia et al. 1993; Kronawitter et al. 2000; Zepf et al. 2000; Côté et al. 2003, hereafter C03) have shown that there appears to be a significant amount of dark matter in NGC 4472.
While the dark matter hypothesis remains the most widely used approach in dynamical modeling of all types of galaxies, MOND (MOdified Newtonian Dynamics) (Milgrom 1983) theory has also been shown to successfully fit the rotation curves of spiral galaxies, although when applied to earlytype galaxies the results have been ambiguous. The basic assumption of MOND is that the acceleration due to gravitational force does not simply depend upon the mass m as in the Newtonian approach, but that it has a more complex form, m/μ(a/a_{0}), where a_{0} is a universal constant, taken to be cm s^{2} (Famaey et al. 2007) and μ(a/a_{0}) is the MOND interpolating function (see Sect. 3.2 for the details on various interpolating functions used in this paper). For a ≫ a_{0}, the Newtonian acceleration is restored, whereas for a ≪ a_{0} the interpolation function becomes μ = a/a_{0}. To model the dynamics of galaxies, it is wellknown that the Jeans equation (for both Newtonian and MOND approaches) is one of the most important tools: the declining trend of the velocity dispersion for some galaxy is interpreted as proof of a constant masstolight ratio in a given galaxy, hence a lack of dark matter. The Jeans equation is widely used when one models the dynamics of the earlytype galaxies (e.g. Romanowsky et al. 2003; Samurović & Danziger 2005) and it has been shown in the case of the Milky Way (Samurović & Lalović 2011) that it can be successfully applied using the MOND approach to spiral systems such as our Galaxy. For ellipticals, the Jeans equation was used by Tiret et al. (2007), who modeled NGC 3379 and found that MOND can reproduce the observations on all scales. For NGC 5128, Samurović (2010) showed that dark matter is not needed out to 6.4 R_{e}, where R_{e} is the effective radius and both Newtonian massfollowslight and MOND models provide good fits to the observed velocity dispersion. Interestingly, Salinas et al. (2011) found that MOND predictions overestimate the observed velocity dispersion (declining profile) for the field elliptical NGC 7507. On the other hand, Richtler et al. (2008) found that the central galaxy of the Fornax cluster, NGC 1399 needs a hypothetical dark halo. We note that for both NGC 3379 and NGC 5128 the velocity dispersion shows a declining trend, which starts at ~200 kms^{1} for NGC 3379 and ~150 kms^{1} for NGC 5128. In addition, Angus et al. (2008) showed that by using the measured velocity distribution of satellite galaxies, obtained for Sloan Digital Sky Survey (SDSS) data, one can successfully model the observed declining lineofsight velocity dispersions. For NGC 1399, one encounters much higher values of the approximately constant value of velocity dispersion: σ ≈ 300 kms^{1}.
Kinematics data for NGC 4472 for the sample of GCs. Upper part: each bin (except for the outermost one) contains approximately 60 objects. Lower part: each bin has a constant width of 2 arcmin (except for the outermost one for which the width is 4 arcmin).
In this paper, we study the brightest optical member of the Virgo cluster NGC 4472, which is a noncentral galaxy (the galaxy at the center of the Virgo cluster is another wellstudied object, NGC 4486 = M 87). We use the globular cluster (GC) data from C03, who estimate that NGC 4472 contains ~6000 GCs projected within 100 kpc from its center as tracers of the mass distribution: the velocity dispersion of this galaxy is approximately constant and equal to ~320 kms^{1} throughout the whole galaxy providing an interesting way of testing the MOND approach on an earlytype galaxy with a high fraction of dark matter in the Newtonian approach. One should note that when solving the Jeans equation, the problem of the massanisotropy (Tonry 1983) may occur: for example, unaccountedfor anisotropies in the orbits of the tracers may mimic the existence (or lack) of the dark component. The example of NGC 4472 provides an opportunity to avoid this problem since it is shown that the orbits of GCs in this galaxy are isotropic (see C03 and also below). Thus, in this paper we study the dynamical isotropic Newtonian and MOND models of NGC 4472, assuming spherical symmetry.
The plan of the paper is as follows: in Sect. 2, we present the observational data related to NGC 4472; in Sect. 3, we solve the isotropic Jeans equation using two Newtonian models (one that takes into account that massfollowslight and the other dark matter), and the MOND approach (in four different models) and study the cumulative mass and the masstolight ratio profiles of NGC 4472, and finally in Sect. 4 we present our conclusions.
2. Observational data
The galaxy NGC 4472 ( = M 49) is one of the most massive galaxies in the Virgo cluster. In this paper, we adopt the distance of 16.7 ± 0.6 Mpc (Blakeslee et al. 2009). At this distance, 1 arcsec corresponds to 81 pc and 1 arcmin corresponds to 4.86 kpc. The effective radius of NGC 4472 is R_{e} = 104′′ ( = 8.42 kpc) and is taken from the RC3 catalog (de Vaucouleurs et al. 1991). This galaxy has a low systemic velocity, v_{gal} = 996 ± 18 kms^{1} (HyperLeda database). The total B magnitude is B_{T} = 9.3 (according to Blakeslee et al. 2009, based on the NED database). Throughout the paper, the Hubble constant h_{0} = 0.70 is assumed.
2.1. Globular clusters of NGC 4472
We analyzed the GCs data set of C03, which contains 263 GCs extending out to ~570 arcsec. We present the kinematics of NGC 4472 based on the GCs in Table 1 and in Fig. 1. The radial velocities of GCs were used to determine the kinematics of NGC 4472: we calculated, using standard definitions, the velocity dispersion, the skewness and kurtosis parameters, s_{3} and s_{4}, which describe the asymmetric and symmetric departures from the Gaussian in each bin. Negative values of kurtosis imply that tangential anisotropies dominate and positive values of kurtosis indicate the dominance of radial orbits. We determined the kinematics of NGC 4472 using two different binnings: (a) in the first one, we placed approximately 60 objects in each bin, except for the outermost bin where we have 84 GCs; and (b) in the second one, we kept the width of the each bin constant (two arcmin), again with the exception of the outermost one, which has the width of four arcmin (we grouped here two outermost bins, because in the bin between eight and ten arcmin only 8 GCs are present). As one can see from Table 1 and Fig. 1, both binnings produce very similar results and in our analysis below we concentrate on the binning “(a)” for simplicity but we also address the binning “(b)” when necessary. From Table 1 and Fig. 1, one can see that the profile of the velocity dispersion is approximately constant with a high value (σ ≈ 320 kms^{1}). One can also see that the values of s_{4} are consistent with zero throughout the whole galaxy in all bins and in all our models isotropy of the orbits will always be assumed (s_{4} = 0, i.e. β_{∗} = 0, see below). We note that when all GCs are put into one bin the s_{4} parameter is again consistent with zero: s_{4} = 0.178 ± 0.302 (and s_{3} = 0.188 ± 0.151 when all GCs are in one bin). This is in agreement with the detailed analysis of C03, who found that they “...have been able to show by direct means that an isotropic velocity ellipsoid is a highly accurate description of GC kinematics”.
Fig. 1 Kinematics of NGC 4472 based on the sample of the GCs. From top to bottom: radial velocity of the GCs in kms^{1}; velocity dispersion in kms^{1} calculated in a given bin and the s_{3} and s_{4} parameters, which describe symmetric and asymmetric departures from the Gaussian, respectively. Two different binnings were used: a) the binning with approximately 60 GCs per bin (except for the outermost one) (full circles) and b) the binning with a fixed width of the bin (2 armin, except for the outermost one where width is 4 arcmin) (open circles). See text for details. 
The ellipticity of the entire system of GCs is approximately 0.16 (Lee et al. 1998), and in our analysis described below we assumed spherical symmetry. This assumption had been previously applied in models of NGC 4472, such as in Loewenstein (1992), Kronawitter et al. (2000), Zepf et al. (2000), and C03, although as noted by Zepf et al. (2000) this “...is not strictly true for NGC 4472, but it not likely to be too far wrong given that the observed ellipticity of the globular cluster systems is 0.1–0.2 (e.g., Lee et al. 1998)”. More specifically, in the case of NGC 4472 we expect that the error in the masstolight ratio made using such an assumption is no greater than approximately 10% as shown in Kronawitter et al. (2000) (Sect. 5.1, where their estimates of uncertainties due to rotation and nonsphericity are presented).
For the Jeans modeling, we also needed the radial surface density profile of the GCs in NGC 4472. The profile that we used is based on the one presented in Rhode & Zepf (2001; see their Fig. 5 and Table 2): they fitted the radial surface density (from their Table 2) with the expression log N_{cl} = (1.91 ± 0.03) − (1.28 ± 0.04)log r, where N_{cl} is the surface density of GCs and r is in arcminutes. They note that this “single power law does not provide a good fit over the entire radial range”. This is why we performed a fitting procedure (on the same data) for the region for which we have GC data. The result is given in Fig. 2 and compared with the result of Rhode & Zepf (“RZ01” in Fig. 2). Our fit, given by log N_{cl} = (1.81 ± 0.05) − (1.10 ± 0.04)log r, resulted in a smaller χ^{2} value (χ^{2} = 0.153) with respect to the fit of Rhode & Zepf (χ^{2} = 0.765) and therefore α = dlnρ/dlnr = −2.10 was used in the models below. As can be seen from Fig. 2, our fit is not radically different from that obtained in Rhode & Zepf (2001) but in our models we wanted to include a fit that represents as closely as possible the distribution of GCs in the region that we analyze and from Fig. 2 one can see that the innermost region (inside ~1 R_{e}) is more accurately described by our fit. The quoted value of the fit by Rhode & Zepf is obtained over the whole radial range (out to ~23 arcmin), whereas our fit is performed for the data interior to ~10 arcmin.
Fig. 2 Radial distribution of GCs in NGC 4472. A power law is fitted to the radial surface density of GCs: N_{cl} ∝ r^{ − γ} (see text for details). The dashed line is the fit from the paper by Rhode & Zepf (2001). 
2.2. Xray data for NGC 4472
The Xrays emission of NGC 4472 has been wellstudied since this is one of the most Xray luminous earlytype galaxies in the local Universe (although it is fainter in Xrays than M 87). Irwin & Sarazin (1996) observed NGC 4472 using the ROSAT HighResolution Imager and determined the temperature profile of this galaxy, obtaining a fit T(r) = 0.93(r/arcmin)^{0.127} keV. Kraft et al. (2011) presented their result for a 100 ks XMMNewton observations of the hot gas in NGC 4472: from their Fig. 9, one can see that the temperature beyond ~5 kpc remains approximately constant, and T ≈ 1 keV. Therefore, in our estimates of the mass (and the masstolight ratio) based on the Xrays we used T = 1 keV.
The total mass of NGC 4472 interior to the radius r based on the Xray observations was calculated (as described in Sect. 3.3) using the equation of Kim & Fabbiano (1995) (1)where the exponent a is related to the temperature of the Xray halo (T ~ r^{ − a}) and in our calculations was taken to be zero for simplicity – this assumption was made because the weak dependence of the temperature on the radius is not expected to influence the estimate of the total mass significantly – and an error in the parameter a of at most 10% is expected; β is the slope used in the analytic King approximation model and is taken to be β = 0.43 (Trinchieri et al. 1986). The masstolight ratio in the Bband based on the Xrays was calculated as a function of radius r as (2)where B is the B magnitude of the galaxy inside radius r (Kim & Fabbiano 1995).
3. Dynamical models
For both approaches, we solve the Jeans equation (e.g. Binney & Tremaine Binney 2008) for a spherical approximation (3)The parameter a_{N;M} is an acceleration term that is different for each approach. In the Newtonian (“N”) approach, it is equal to , and for MOND (“M”) models, a_{M}μ(a_{M}/a_{0}) = a_{N} (in Sect. 3.2, we present the relevant details for the function μ). The parameter σ_{r} is the radial stellar velocity dispersion and α = dlnρ/dlnr is the slope of the tracer surface density ρ as given above; in the models below, we use α = −2.10. The rotation speed v_{rot} of NGC 4472 is taken to be kms^{1} (see C03, Table 3).
One can use a wellknown relation to describe the nonspherical nature of the stellar velocity dispersion, where we write β_{∗} to avoid confusion with the parameter used in Eqs. (1) and (2), given by (4)where and 0 < β_{∗} < 1 means that the orbits are predominantly radial (equivalent to s_{4} > 0), whereas for − ∞ ≤ β_{∗} < 0 the orbits are mostly tangential (equivalent to s_{4} < 0) (Gerhard 1993). As stressed above, the galaxy NGC 4472 represents a case for which isotropic distribution of the orbits is wellestablished (C03) and our estimate of the s_{4} parameter strongly suggests the absence of departures from a Gaussian, i.e. s_{4} ≈ 0 (Fig. 1).
In all models calculated below, we used the projected lineofsight velocity dispersion (e.g. Binney & Mamon 1982) given as (5)where the truncation radius, r_{t}, extends beyond the observed kinematical point of the largest galactocentric radius.
Cappellari et al. (2006) studied the correlation between masstolight ratio and velocity dispersion in the Iband in the inner parts where baryonic matter is expected to dominate for the sample of earlytype galaxies (NGC 4472 was not included in their sample). They found, using twointegral Jeans and threeintegral Schwarzschild dynamical models, that the expression that holds is (6)where the value of the velocity dispersion is assumed to that evaluated at the effective radius. When we insert σ_{e} = 310 kms^{1}, we obtain M/L_{I} = 5.49 ± 0.38. Since for NGC 4472, R − I = 0.55 (from the HyperLeda database), the masstolight ratio in the Rband becomes M/L_{R} = 6.66 ± 0.46. This value is somewhat higher than that obtained using the stellar population synthesis (SPS) models of Bell & de Jong (2001); using Kodama & Arimoto (1997) models based on the Salpeter initial mass function (IMF), one finds M/L_{R} = 5.42. Finally, to establish the masstolight ratio in the Bband, we use B − R = 1.64 (again from the HyperLeda database) to find M/L_{B} = 11.36 ± 0.78. This value is higher than that found by van der Marel (1991), M/L_{B} = 8.31 ± 0.17 (converted to the distance used in this paper). We note that from the SPS analysis of Bell & de Jong (see their Fig. 3) again using Kodama & Arimoto (1997) models based on the Salpeter IMF, one obtains M/L_{B} = 9.77, which is very close to the value used in our Jeans models. We find that the Salpeter IMF best explains the modeled M/L ratio (see also the discussion below in Sect. 3.2). In our models below, we always use M/L_{B} = 10.00, assuming that it includes only a stellar component, unless noted otherwise. In the discussion below (Sect. 3.3), we present several estimates of the masstolight ratio of NGC 4472 based on different methodologies.
3.1. Newtonian models
First, we develop a constant masstolight ratio model and consider that the stellar mass is distributed in the form of the standard Hernquist (1990) profile (7)which has two parameters of the total mass M_{T} and scale length, a, where R_{e} = 1.8153a. As stated above, we assume isotropy (β_{∗} = 0) and a constant masstolight ratio, M/L_{B} = 10.00, and the result of this model (using 28 < v_{rot} < 105 kms^{1}) is given in Fig. 3: the thick solid line is for v_{rot} = 50 kms^{1} and the hatched region encompassing this value is determined by v_{rot} = 28 kms^{1} (lower limit) and v_{rot} = 105 kms^{1} (upper limit). One can see that for this massfollowslight approximation, the observed velocity dispersion can be fitted interior to ~1 R_{e} and beyond that region dark matter begins to play an important dynamical role.
Fig. 3 Isotropic Jeans modeling of the projected velocity dispersion of NGC 4472 using a Newtonian approach: massfollowslight (M/L_{B} = 10) and NFW models are shown. One effective radius is indicated by a vertical dotted line. Solid and dashed lines are plotted using v_{rot} = 50 kms^{1}. A scaledup version of the Newtonian massfollowslight with M/L_{B} = 15 is plotted with a thin solid line. Two hatched regions (the upper one is for the NFW model and the lower one is the Newtonian constant masstolight ratio model with M/L_{B} = 10) for each approach take into account various values of the rotational velocity: the lower limit in both cases is for v_{rot} = 28 kms^{1} and the upper limit in both cases is for v_{rot} = 105 kms^{1}. 
Next, we create a Newtonian model with dark matter (again assuming 28 < v_{rot} < 105 kms^{1}). Numerous models are possible (see, for example, Bullock et al. 2001, for various profiles) but here we chose to present the model (again assuming isotropy), which includes the NFW (Navarro et al. 1997) halo with the following fitted parameters: ρ_{s} = 0.13 ± 0.02 M_{⊙} pc^{3} and r_{s} = 144 ± 5 arcsec ( ≈ 11.7 ± 0.4 kpc). This halo is added to a Hernquist spheroid that includes the stellar matter. The de Vaucouleurs profile is a fair representation of the observed stellar distribution in NGC 4472 between 7 arcsec and ~4 arcmin, where it begins to diverge although not dramatically (see Kim et al. 2000, Fig. 8), hence the use of the Hernquist profile to describe the de Vaucouleurs law is a useful approximation.This model is presented in Fig. 3: the dashed line represents v_{rot} = 50 kms^{1} and the hatched region encompassing it is determined as in the constant masstolight ratio model above. This model obviously provides a satisfactory fit to the observed velocity dispersion. To fit the outermost point (binning “a”, at ~4.2 R_{e}), one would need M/L_{B} ≈ 20, which means that interior to that point the contribution of dark matter is equal to that of baryonic matter. The quality of the fits is given in Table 2. The obtained fitted parameters imply the virial mass of ~2 × 10^{12} M_{⊙} and the high concentration parameter c = 15.5^{1} of the NFW model, which agrees with the results of the highresolution Nbody simulations of a ΛCDM cosmology (Bullock et al. 2001) for haloes at z = 0.
We also made a scaledup version of the Newtonian massfollowslight model with M/L_{B} = 15, which is plotted with a thin solid line in Fig. 3 and one can see that the quality of this fit is poorer than that of the model with the NFW halo ( and for the binnings “a” and “b” as described earlier). Although high, this value of the masstolight ratio is not inconsistent with the PEGASE2 models^{2} with x = −1.85 (where x is the logarithmic slope used in the parametrizations of the IMF, see also Sect. 3.2), for which for B − R ≈ 1.6 one obtains M/L_{B} ≈ 15.8 (see Bell & De Jong 2001, Fig. 4).
Taking both Newtonian models into account (and assuming the rotational velocity 28 < v_{rot} < 105 kms^{1}), one may reach the conclusion that a significant amount of dark matter is indeed needed in NGC 4472, which agrees with the results found in the literature.
3.2. MOND models
As mentioned in the introduction, we tested several MOND models using the Jeans equation assuming isotropy, a spherical approximation, and a constant masstolight ratio, M/L_{B} = 10.
The Newtonian acceleration is given as a_{N} = aμ(a/a_{0}), where a is the MOND acceleration, a_{0} is the universal constant given in the introduction, and μ(x) is the MOND interpolating function that is different for each MOND model. The following MOND models were tested in order to determine whether an additional dark matter component is needed or MOND alone can correctly describe the dynamics of NGC 4472: (i) the “simple” MOND formula from Famaey & Binney (2005), (ii) the “standard” formula (Sanders & McGaugh 2002), (iii) the Bekenstein’s “toy” formula (Bekenstein 2004), and (iv) the Zhao (2007) formula. The expressions for the circular velocities for the first three functions are given in Samurović & Ćirković (2008):
 1.
 2.
 3.

4.
Finally, we analyze the MOND model conceived by Zhao (2007); this model was created to explain the exact amount of dark energy and it belongs to the same family of MOND models as the “simple” model: (11)While the “simple” model is obtained for n = 1, the Zhao (2007) model is obtained for n = 3. Calculation of the circular velocity for n = 3 is much more complicated and is therefore done numerically.
Fig. 4 Isotropic MOND models of the projected velocity dispersion of NGC 4472. One effective radius is indicated by a vertical dotted line. In this figure for all the models different types of lines indicate that the rotational velocity v_{rot} = 50 kms^{1} is assumed. The hatched regions are determined by the rotational velocity of 28 kms^{1} (lower limit) and 105 kms^{1} (upper limit). To avoid confusion, the Zhao model has been labeled with “Z” on the right hand side of the figure. 
Fits for the Newtonian (constant masstolight ratio and NFW models) and the MOND isotropic models (assuming v_{rot} = 50 kms^{1}) for two different binnings (“a” and “b” as described in Sect. 2.1).
We again present the models for which the rotational velocity is 28 < v_{rot} < 105 kms^{1}. The quality of the fits is given in Table 2 for v_{rot} = 50 kms^{1}. From Fig. 4, one can see that all MOND models can provide a successful fit interior to ~ 2 R_{e}. Beyond ~2 R_{e}, only the “toy” model can provide a fit (χ^{2} = 1.53)^{3} to the observed velocity dispersion points, whereas the “simple” MOND model marginally fails to fit the outermost point (for the binning “a”, χ^{2} = 4.06). This is why we decided to present an additional fit of the “simple” model for which we tested a higher value of the masstolight ratio while still taking into account Eq. (6) (see also below) of M/L_{B} = 12. One can see that this higher value of the masstolight ratio (assuming no contribution of dark matter) manages to provide a fit to all the observed velocity dispersion points; the obtained fit is, together with the fit obtained using the “toy” MOND model, the closest of all tested MOND models, with χ^{2} = 1.58 for the binning “a”. The remaining two MOND models would need to have an anisotropic distribution of the orbits (tangential anisotropies with β_{∗} ~ −0.2, which are however inconsistent with either the available kinematics of NGC 4472 or a higher value of the masstolight ratio (implying that there is an additional dark component or the IMF based on the PEGASE2 models with x = −1.85, see Bell & de Jong 2001). In all our models, we set the masstolight ratio to 10 in the Bband (except for one case, where M/L_{B} = 12 was taken) and this was the only free parameter in our model. All other quantities entering the Jeans equation were fixed: β_{ ∗ } = 0 (isotropy was always assumed), α = −2.10, and 28 < v_{rot} < 105 kms^{1}. The resulting χ^{2} values (for v_{rot} = 50 kms^{1}) can be found in Table 2.
We note that in our MOND models we did not take into account the “external field effect”, which has implications for nonisolated systems. We used the results of Richtler et al. (2011), who estimated the external gravitational fields acting on Virgo and Coma galaxies. The galaxy NGC 4472 was on their list and they calculated the value of the external field at the position of this galaxy to be ~0.1a_{0} (see their Fig. 3) hence neglected the influence of this effect because this value was smaller than the values estimated in the inner parts of their galaxies. In a similar manner, we calculated the acceleration in the inner parts of NGC 4472, which we analyzed finding that indeed , where is the acceleration in the inner parts (interior to ~10 arcmin) and is the acceleration of the external gravitational field acting on NGC 4472.
3.3. The cumulative mass and the masstolight ratio of NGC 4472
We now present the estimates of the cumulative mass and the masstolight ratio based on various methodologies. The results are given in Fig. 5, where we provide the results based on two Newtonian estimates, namely one for Xray data and the socalled “tracers mass estimator” (TME; Evans et al. 2003), and four MOND formulas are given. The estimates based on Xrays are derived from Eqs. (1) and (2). The TME methodology enables one to estimate the enclosed mass from the projected positions and the lineofsight velocities of a given tracer population (such as GCs in the case of NGC 4472). The usual assumptions are that the tracer population is spherically symmetric and has a number density that obeys a power law, ρ(r) = ρ_{0}(a/r)^{γ}, where a is constant and γ is a parameter determined from the surface density of the tracer population (see details and comparisons with other techniques in e.g. Samurović 2007, Chap. 5), and these assumptions were adopted in the present work. The advantage of this method is independent method can be used to establish the total galaxy mass for various galaxies at large radii, and the disadvantage is the small number of systems that have been studied in this way so far. Additional observations of various tracers for a larger sample of earlytype galaxies will provide the opportunity to compare their independent estimates of cumulative masses.At ~7 arcmin (the outermost velocity dispersion point), one can see that the masstolight ratio is M/L_{B} ~ 20 in the Newtonian TME model, which is consistent with the result found in the Jeans modeling and shown in Fig. 3 and Table 2 (assuming rotational velocity 28 < v_{rot} < 105 kms^{1}).
In Fig. 5, one can see that interior to ~6 arcmin there is a discrepancy between the Xray and TME masses, which however decreases in the outermost regions (10 arcmin). The reason for this discrepancy might be a lack of hydrostatic equilibrium in NGC 4472: for an example of this we refer by Diehl & Statler (2007), who tabulate the position angles in the optical and Xray domains, PA_{opt} = 160 ± 0.5 and PA_{X} = 136 ± 13 degrees, respectively. This discrepancy warrants a further study. The Xray mass estimate at ~6 arcmin is much lower than that based on the TME methodology and does not suggest the existence of dark matter within this radius (inferred using the Jeans equation). Thus, from our Newtonian approach one may conclude that the TME methodology provides a more accurate estimate of the cumulative mass of NGC 4472.
Using both Newtonian estimates of the cumulative mass (TME and Xrays), the overall conclusion is that the estimated mass (and the masstolight ratio) clearly indicates that dark matter exists in the outermost regions. On the other hand, the predictions of various MOND formulas indicate that the mass and the masstolight ratio are consistent with the stellar population only throughout the whole galaxy. In Fig. 5, the estimates of the mass (and the masstolight ratio) based on the MOND models used in this paper are also given. The masses were calculated using the circular velocities of each model. All masstolight ratios of the MOND models tested in this paper tend to group, in the outermost region, around M/L_{B} ~ 10, except for the “toy” model, which has a somewhat higher value, M/L_{B} ~ 13. The results of the Jeans modeling presented in Fig. 4 imply that the “simple” MOND model with 10 < M/L_{B} < 12 together with the “toy” model with M/L_{B} = 10 provide the best fits of all tested MOND models to the observed velocity dispersion (assuming the rotational velocity of 50 kms^{1}). These values of the masstolight ratio are consistent with (a) the SPS analysis shown above and (b) the estimate obtained by calculating the cumulative masstolight ratio of NGC 4472 using the MOND formulas (shown in Fig. 5). A detailed analysis of the stellar populations in NGC 4472 will provide an insight into the stellar masstolight ratio of this galaxy and help us to disentangle the contribution of dark matter and thus to more tightly constrain the parameters necessary for the analysis of various MOND models (Samurović et al., in prep.).
Fig. 5 Cumulative Newtonian and MOND mass (upper panel) and Newtonian and MOND masstolight ratio (lower panel). In the Newtonian approach we used the Xrays with temperature T = 1 keV and the “tracer mass estimator” (TME). In the MOND approach we tested four different models: “simple”, “standard”, Bekenstein’s “toy” and Zhao model (see text for details). 
4. Conclusions
We have studied the dynamics of the earlytype galaxy NGC 4472 using the spherical Jeans equation assuming isotropy: the velocity dispersion out to ~7 arcmin (=4.2 R_{e}) was modeled in order to show whether the MOND theory can provide a successful fit to the observed data. Our conclusions are as follows:

1.
The kinematics of NGC 4472 were studied andthe velocity dispersion was found to remain approximatelyconstant at 320 kms^{1} throughout the whole galaxy. Wecalculated both skewness and kurtosis, which are both consistentwith zero implying that the orbits of GCs inNGC 4472 are close to isotropic. In all of ourdynamical models, we assumed isotropy to avoid themassanisotropy degeneracy which may lead to incorrectconclusions about the existence of dark matter (and the validity ofthe MOND theory) in earlytype galaxies.

2.
We solved the Jeans equation for two isotropic Newtonian models: massfollowslight and NFW models, both in the spherical approximation. We found that NGC 4472 indeed contains a significant amount of dark matter: interior to ~4.2 R_{e}, the contribution of dark matter is approximately equal to that of baryonic matter. To fit the outermost radial point, a high masstolight ratio of M/L_{B} ≈ 20 is needed. Using M/L_{B} = 10 based only on the baryonic matter (and in agreement with the SPS models), one can fit only the inner region of NGC 4472 (~1 R_{e}).

3.
We also solved the Jeans equation for a spherical approximation for four different MOND models, again assuming isotropy. We found (assuming a rotational velocity 28 < v_{rot} < 105 kms^{1}) that the “simple” (with 10 < M/L_{B} < 12) and “toy” (with M/L_{B} = 10) models can provide successful fits to the observed velocity dispersion, and we conclude that these models are consistent with a no darkmatter hypothesis. Two other tested MOND models, the “standard” and Zhao (2007) models could provide a fit to the observed velocity dispersion interior to ~2 R_{e}, whereas in the outer parts one would need to invoke anisotropies and/or additional dark matter. The future observations of tracers such as GCs in NGC 4472 will hopefully determine more precisely and to larger radii the dispersion and the rotational velocity in this galaxy thus providing a possibility of more tightly constraining dynamical models of this earlytype galaxy.
This value is high compared to the most recent simulations which use the latest Wilkinson Microwave Anisotropy Probe (WMAP) parameters (e.g. Duffy et al. 2008).
PEGASE2 (http://www2.iap.fr/users/fioc/PEGASE.html) is the code used to compute the spectral evolution of galaxies: using this code, one follows the evolution of the stars, gas, and metals for a law of star formation and a stellar IMF.
Acknowledgments
This work was supported by the Ministry of Education and Science of the Republic of Serbia through the project No. 176021, “Visible and Invisible Matter in Nearby Galaxies: Theory and Observations”. This research made use of the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. We acknowledge the usage of the HyperLeda database (http://leda.univlyon1.fr). The author gratefully acknowledges the valuable comments of the referee, which helped to significantly improve the quality of the manuscript and correct the errors present in its initial version.
References
 Angus, G. W., Famaey, B., Tiret, O., Combes, F., & Zhao, H.S. 2008, MNRAS, 383, L1 [NASA ADS] [Google Scholar]
 Bell, E. F., & de Jong, R. S. 2001, ApJ, 550, 212 [NASA ADS] [CrossRef] [Google Scholar]
 Bekenstein, J. 2004, Phys. Rev. D, 70, 083509 [NASA ADS] [CrossRef] [Google Scholar]
 Binney, J. J., & Mamon, G. 1982, MNRAS, 200, 361 [NASA ADS] [CrossRef] [Google Scholar]
 Binney, J. J., & Tremaine, S. 2008, Galactic Dynamics, second Edition (Princeton: Princeton Univ. Press) [Google Scholar]
 Blakeslee, J. P., Jordán, A., Mei, S., et al. 2009, ApJ, 694, 556 [NASA ADS] [CrossRef] [Google Scholar]
 Bullock, J. S., Kolatt, T. S., Sigad, Y., et al. 2001, MNRAS, 321, 559 [Google Scholar]
 Cappellari, M., Bacon, R., Bureau, M., et al. 2006, MNRAS, 366, 1126 [NASA ADS] [CrossRef] [Google Scholar]
 Côté, P., McLaughlin, D. E., Cohen, J. G., & Blakeslee, J. P. 2003, ApJ, 591, 850 (C03) [NASA ADS] [CrossRef] [Google Scholar]
 de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H. G. Jr., et al. 1991, Third Reference Catalogue of Bright Galaxies (New York: SpringerVerlag) [Google Scholar]
 Diehl, S., & Statler, T. S. 2007, ApJ, 668, 150 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Duffy, A. R., Schaye, J., Kay, S. T., & Dalla Vechia, C. 2008, MNRAS, 390, L64 [NASA ADS] [CrossRef] [Google Scholar]
 Evans, N. W., Wilkinson, M. I., Perrett, K. M., & Bridges, T. J. 2003, ApJ, 583, 752 [NASA ADS] [CrossRef] [Google Scholar]
 Famaey, B., & Binney, J. 2005, MNRAS, 363, 603 [NASA ADS] [CrossRef] [Google Scholar]
 Famaey, B., Gentile, G., Bruneton, J.P., & Zhao, H. S. 2007, Phys. Rev. D, 75, 063002 [NASA ADS] [CrossRef] [Google Scholar]
 Gerhard, O. 1993, MNRAS, 265, 213 [NASA ADS] [CrossRef] [Google Scholar]
 Hernquist, L. 1990, ApJ, 356, 359 [NASA ADS] [CrossRef] [Google Scholar]
 Irwin, J. A., & Sarazin, C. L. 1996, ApJ, 471, 683 [NASA ADS] [CrossRef] [Google Scholar]
 Kim, D.W., & Fabbiano, G. 1995, ApJ, 441, 182 [NASA ADS] [CrossRef] [Google Scholar]
 Kim, E., Lee, M. G., & Geisler, D. 2000, MNRAS, 314, 307 [NASA ADS] [CrossRef] [Google Scholar]
 Kodama, T., & Arimoto, N. 1997, A&A, 320, 41 [NASA ADS] [Google Scholar]
 Kraft, R. P., Forman, W. R., Jones, C., et al. 2011, ApJ, 727, 41 [NASA ADS] [CrossRef] [Google Scholar]
 Kronawitter, A., Saglia, R. P., Gerhard, O., & Bender, R. 2000, A&AS, 144, 53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lee, M. G., Kim, E., & Geisler, D. 1998, AJ, 115, 947 [NASA ADS] [CrossRef] [Google Scholar]
 Loewenstein, M. 1992, ApJ, 384, 474 [NASA ADS] [CrossRef] [Google Scholar]
 Milgrom, M. 1983, ApJ, 270, 365 [NASA ADS] [CrossRef] [Google Scholar]
 Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 [NASA ADS] [CrossRef] [Google Scholar]
 Rhode, K. L., & Zepf, S. E. 2001, AJ, 121, 210 [NASA ADS] [CrossRef] [Google Scholar]
 Richtler, T., Schuberth, Y., Hilker, M., et al. 2008, A&A, 478, L23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Richtler, T., Famaey, B., Gentile, G., & Schuberth, Y. 2011, A&A, 531, A100 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Romanowsky, A. J., Douglas, N. G., Arnaboldi, M., et al. 2003, Science, 5640, 1696 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Saglia, R. P., Bertin, G., Bertola, F., et al. 1993, ApJ, 403, 567 [NASA ADS] [CrossRef] [Google Scholar]
 Salinas, R., Richtler, T., Bassino, L. P., Romanowsky, A. J., & Schuberth, Y. 2012, A&A, 538, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Samurović S. 2007, Dark Matter in Elliptical Galaxies, Publications of the Astronomical Observatory of Belgrade, 81 [Google Scholar]
 Samurović, S. 2010, A&A, 514, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Samurović, S., & Ćirković, M. M. 2008, A&A, 488, 873 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Samurović, S., & Danziger, I. J. 2005, MNRAS, 363, 769 [NASA ADS] [CrossRef] [Google Scholar]
 Samurović, S., & Lalović, A. 2011, A&A, 531, A82 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sanders, R. H., & McGaugh, S. 2002, ARA&A, 40, 263 [NASA ADS] [CrossRef] [Google Scholar]
 Tiret, O., Combes, F., Angus, G. W., Famaey, B., & Zhao, H. S. 2007, A&A, 476, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tonry, J. L. 1983, ApJ, 266, 58 [NASA ADS] [CrossRef] [Google Scholar]
 Trinchieri, G., Fabbiano, G., & Canizares, C. R. 1986, ApJ, 310, 637 [NASA ADS] [CrossRef] [Google Scholar]
 van der Marel, R. P. 1991, MNRAS, 253, 710 [NASA ADS] [CrossRef] [Google Scholar]
 Zepf, S. E., Beasley, M. A., Bridges, T. J., et al. 2000, AJ, 120, 2928 [NASA ADS] [CrossRef] [Google Scholar]
 Zhao, H. S. 2007, ApJ, 671, L1 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Kinematics data for NGC 4472 for the sample of GCs. Upper part: each bin (except for the outermost one) contains approximately 60 objects. Lower part: each bin has a constant width of 2 arcmin (except for the outermost one for which the width is 4 arcmin).
Fits for the Newtonian (constant masstolight ratio and NFW models) and the MOND isotropic models (assuming v_{rot} = 50 kms^{1}) for two different binnings (“a” and “b” as described in Sect. 2.1).
All Figures
Fig. 1 Kinematics of NGC 4472 based on the sample of the GCs. From top to bottom: radial velocity of the GCs in kms^{1}; velocity dispersion in kms^{1} calculated in a given bin and the s_{3} and s_{4} parameters, which describe symmetric and asymmetric departures from the Gaussian, respectively. Two different binnings were used: a) the binning with approximately 60 GCs per bin (except for the outermost one) (full circles) and b) the binning with a fixed width of the bin (2 armin, except for the outermost one where width is 4 arcmin) (open circles). See text for details. 

In the text 
Fig. 2 Radial distribution of GCs in NGC 4472. A power law is fitted to the radial surface density of GCs: N_{cl} ∝ r^{ − γ} (see text for details). The dashed line is the fit from the paper by Rhode & Zepf (2001). 

In the text 
Fig. 3 Isotropic Jeans modeling of the projected velocity dispersion of NGC 4472 using a Newtonian approach: massfollowslight (M/L_{B} = 10) and NFW models are shown. One effective radius is indicated by a vertical dotted line. Solid and dashed lines are plotted using v_{rot} = 50 kms^{1}. A scaledup version of the Newtonian massfollowslight with M/L_{B} = 15 is plotted with a thin solid line. Two hatched regions (the upper one is for the NFW model and the lower one is the Newtonian constant masstolight ratio model with M/L_{B} = 10) for each approach take into account various values of the rotational velocity: the lower limit in both cases is for v_{rot} = 28 kms^{1} and the upper limit in both cases is for v_{rot} = 105 kms^{1}. 

In the text 
Fig. 4 Isotropic MOND models of the projected velocity dispersion of NGC 4472. One effective radius is indicated by a vertical dotted line. In this figure for all the models different types of lines indicate that the rotational velocity v_{rot} = 50 kms^{1} is assumed. The hatched regions are determined by the rotational velocity of 28 kms^{1} (lower limit) and 105 kms^{1} (upper limit). To avoid confusion, the Zhao model has been labeled with “Z” on the right hand side of the figure. 

In the text 
Fig. 5 Cumulative Newtonian and MOND mass (upper panel) and Newtonian and MOND masstolight ratio (lower panel). In the Newtonian approach we used the Xrays with temperature T = 1 keV and the “tracer mass estimator” (TME). In the MOND approach we tested four different models: “simple”, “standard”, Bekenstein’s “toy” and Zhao model (see text for details). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.