arxiv:astro-ph/ v2 17 Sep 2001

Transcription

1 Mon. Not. R. Aston. Soc. 000, (0000) Pinted 9 August 018 (MN LATEX style file v1.4) Velocity dispesions of dwaf spheoidal galaxies: dak matte vesus MOND Ewa L. Lokas Nicolaus Copenicus Astonomical Cente, Batycka 18, Wasaw, Poland axiv:asto-ph/ v 17 Sep August INTRODUCTION ABSTRACT We pesent pedictions fo the line-of-sight velocity dispesion pofiles of dwaf spheoidal galaxies and compae them to obsevations in the case of the Fonax dwaf. The pedictions ae made in the famewok of standad dynamical theoy of spheical systems with diffeent velocity distibutions. The stas ae assumed to be distibuted accoding to Sésic laws with paametes fitted to obsevations. We compae pedictions obtained assuming the pesence of dak matte haloes (with density pofiles adopted fom N-body simulations) vesus those esulting fom Modified Newtonian Dynamics (MOND). If the anisotopy of velocity distibution is teated as a fee paamete obsevational data fo Fonax ae epoduced equally well by models with dak matte and with MOND. If stella mass-to-light atio of 1M /L is assumed, the equied mass of the dak halo is M, two odes of magnitude bigge than the mass in stas.the deivedmondacceleationscale is a 0 = cm/s. In both cases a cetain amount of tangential anisotopy in the velocity distibution is needed to epoduce the shape of the velocity dispesion pofile in Fonax. Key wods: methods: analytical galaxies: dwaf galaxies: fundamental paametes galaxies: kinematics and dynamics Dwaf spheoidal (dsph) galaxies povide a unique testing field fo the theoies of stuctue fomation. Due to thei lage velocity dispesions they ae believed to be dominated by dak matte (heeafte DM, fo a eview see Mateo 1997). Because of thei small masses they also lie in the egime of small acceleations, whee, accoding to an altenative theoy of Milgom (1983a), the Modified Newtonian Dynamics (MOND) may opeate causing a simila incease of velocity dispesions. Until ecently the measuements of velocity dispesion in dwafs wee esticted to a single value, the cental velocity dispesion, on which all of the dynamical modelling was based. In ode to estimate the mass-to-light atio in such case a simplifying assumption of isotopic velocity distibution was usually made. Using this appoach many authos (see Mateo 1998 and efeences theein) estimated mass-tolight atios in dwafs to be much lage than stella values and intepeted this as a signatue of the pesence of DM. Howeve, Milgom (1995) analysed the data concening seven well known dsph galaxies and concluded that when MOND is applied thei mass-to-light atios ae consistent with stella values if all the obsevational eos ae popely taken into account. Given the velocity dispesion data as a function of distance, which ae becoming available now, we ae able to extact much moe infomation fom the obsevations. Ou conclusions do not have to ely on the assumed velocity anisotopy model, but the best-fitting model can actually be found. On the theoetical side some pogess has also been achieved. The DM distibution in haloes can now be eliably modelled using the esults of N-body simulations (e.g. Navao, Fenk & White 1997, heeafte NFW). The pupose of this wok is to exploe the two hypotheses, based on DM and on MOND, by calculating the pedicted velocity dispesion pofiles and compaing them to available data. In the case of DM ou analysis is somewhat simila to an ealie study of Pyo & Komendy (1990) who also consideed two-component models involving stas and DM. In the case of MOND, ou appoach hee is elated to that of Milgom (1984) who applied the Jeans equation to study density pofiles of equilibium systems with isotopic and isothemal velocity distibutions. Hee we assume instead the obseved density distibution of stas to obtain the velocity dispesions in MOND fo diffeent velocity anisotopy models. VELOCITY DISPERSION PROFILE The adial velocity dispesion σ () of stas can be obtained by solving the Jeans equation (Binney & Temaine 1987) c 0000 RAS

2 Ewa L. Lokas d d (νσ )+ β νσ = νg, (1) whee ν() is a 3D density of stas, g() is the gavitational acceleation associated with the potential Φ(), g = dφ/d, and β = 1 σ θ()/σ () is a measue of the anisotopy in the velocity distibution. In the absence of diect measuements of velocity anisotopy in dsph galaxies we have to conside diffeent models fo β. Ou choice hee is guided by the esults of N-body simulations of DM haloes and obsevations of elliptical galaxies. Fo DM haloes Thomas et al. (1998) find that, in a vaiety of cosmological models, the atio σ θ /σ is not fa fom unity and deceases slowly with distance fom the cente to each 0.8 at the viial adius. Obsevations of elliptical galaxies ae also consistent with the obits being isotopic nea the cente and somewhat adially anisotopic fathe away, although cases with tangential anisotopy ae also obseved (Gehad et al. 001). Ou fiducial anisotopy model will theefoe be that of isotopic obits: σ θ () = σ () and β = 0. A moe ealistic appoximation is povided by a model poposed by Osipkov (1979) and Meitt (1985) with β dependent on distance fom the cente of the object β OM =, () +a whee a is the anisotopy adius detemining the tansition fom isotopic obits inside to adial obits outside. This model coves a wide ange of possibilities fom isotopic obits ( a ) to adial obits ( a 0). Anothe possibility is that of tangential obits with the fiducial case of cicula obits when β. To cove intemediate cases between the isotopic and cicula obits we will conside a simple model of β =const< 0. The solution of the Jeans equation (1) with the bounday condition σ 0 at fo β OM anisotopy is νσ (β = β OM) = 1 + a while fo β=const we find ( + a)νg d, (3) νσ(β = const) = β β νg d. (4) Fom the obsevational point of view, an inteesting, measuable quantity is the line-of-sight velocity dispesion obtained fom the 3D velocity dispesion by integating along the line of sight (Binney & Mamon 198) σ los(r) = I(R) R ( ) 1 β R νσ d, (5) R whee I(R) is the suface bightness. Fo cicula obits, σ = 0, and one has σlos(r) = 1 ( ) R νv d, (6) I(R) R R whee V = (g) 1/ is the cicula velocity associated with the gavitational acceleation g(). Intoducing esults (3) o (4) into equation (5) and inveting the ode of integation the calculations of σ los can be educed to one-dimensional numeical integation of a fomula involving elementay functions in the case of β OM (see Mamon & Lokas, in pepaation) and special functions fo abitay β =const. We will now discuss diffeent factos enteing in the calculations of σ los, namely the density distibution of stas ν, the mass distibutions contibuted by stas and DM to gavitational acceleation g and how this acceleation is changed in MOND..1 Stas The distibution of stas is well known in a numbe of dsph galaxies fom suface bightness measuements and is usually well fitted by the exponential law o its genealization i.e. the Sésic pofile (Sésic 1968, see also Ciotti 1991) I(R) = I 0exp[ (R/R S) 1/m ], (7) whee I 0 is the cental suface bightness and R S is the chaacteistic pojected adius of the Sésic pofile. The Sésic pofile with paamete m vaying in the ange 1 m 10 has been found (Caon, Capaccioli & D Onofio 1993) to fit suface bightness pofiles of elliptical galaxies in lage mass ange bette than the standad de Vaucouleus law (coesponding to m = 4). Fo dsph systems the fits ae usually pefomed with m = 1, although is some cases othe values of m ae found to povide bette fits (e.g. Caldwell 1999). The 3D luminosity density ν() is obtained fom I(R) by depojection ν() = 1 π di dr (8) dr R. In the case of m = 1 we get ν(,m = 1) = I 0K 0(/R S)/(πR S), whee K 0(x) is the modified Bessel function of the second kind. Fo othe values of m in the ange 1/ m 10 an excellent appoximation fo ν() is povided by (Lima Neto, Gebal & Máquez 1999) ( ) pexp ( ) ν() = ν 0 [ ] 1/m (9) R S R S ν 0 = I 0Γ(m) R SΓ[(3 p)m] p = /m /m. The mass distibution of stas following fom (9) is γ[(3 p)m,(/r S) 1/m ] M () = ΥL tot, (10) Γ[(3 p)m] whee Υ=constis the mass-to-light atio fo stas, L tot is the total luminosity of the galaxy and γ(α,x) = x 0 e t t α 1 dt is the incomplete gamma function.. Dak matte We assume that the density distibution of DM is given by the so-called univesal pofile (NFW, see also Lokas & Mamon 001) ρ() δ cha = (11) ρ cit,0 (/ s)(1+/ s) with a single fitting paamete δ cha, the chaacteistic density. The scale adius s is defined by s = v/c whee v is the viial adius (the distance fom the cente of the halo within which the mean density is 00 times the pesent citical density, ρ cit,0) and c is the concentation paamete

3 Velocity dispesions of dwaf spheoidals 3 elated to the chaacteistic density by δ cha = 00c 3 g(c)/3 with g(c) = 1/[ln(1+c) c/(1+c)]. The value of c is known to depend on mass and will be extapolated hee to small masses chaacteistic of dwaf galaxies assuming the following fomula fitted to N-body esults fo ΛCDM cosmology of Jing & Suto (000) ( ) hmv c = 10.3, (1) 10 1 M whee M v = 800πvρ 3 cit,0/3 is the mass within viial adius and h is the Hubble constant in units of 100 km s 1 Mpc 1. In the following we will use h = 0.7. The mass distibution associated with the NFW pofile (11) is M NFW() = g(c)m v [ ln(1+cs) cs 1+cs ], (13) whee s = / v and fo h = 0.7 the viial adius scales with the viial mass so that v = 06[M v/(10 1 M )] 1/3 kpc..3 MOND In the case of MOND the gavitational acceleation is changed so that (Milgom 1984) g M = g N µ(g N/a 0), (14) whee g M and g N ae the modified and Newtonian gavitational acceleations, espectively. a 0 is the chaacteistic acceleation scale of MOND estimated to be of the ode of 10 8 cm/s. The function µ(x) has to epoduce Newtonian dynamics at lage acceleations, i.e. µ(x) 1 fo x, and is supposed to be µ(x) x fo x 0. A good appoximation is then povided by e.g. µ(x) = x(1+x ) 1/ (Milgom 1983b). Using this appoximation we find fom equation (14) [ ] g g M = N +g N(4a 0 +gn) 1/ 1/. (15) 3 RESULTS FOR THE FORNAX DWARF The stuctual paametes of the Fonax dwaf needed in the following computation wee taken fom Iwin & Hatzidimitiou(1995, heeafte IH). They fitted the suface bightness pofile of Fonaxwith the exponential distibution(m = 1in equation [7]) finding the Sésic adius ( e in thei notation) R S = 9.9 acmin = 0.35 kpc assuming the distance modulus m M = 0.4. The bightness of Fonax in V-band is M V = 13.0 mag so the total luminosity needed in equation (10)is L tot,v = L.Wealsoadoptthemass-to-light atio fo stas in this band to be Υ V 1M /L (Mateo et al. 1991). Figue 1 shows the pedicted line-of-sight velocity dispesion pofiles fo the Fonax dwaf togethe with the data fom Mateo (1997). The two panels pesent esults fo diffeent velocity anisotopy models. The uppe panel gives the esults in the case of isotopic obits obtained fom equations (3) and (5) in the limit of lage a. The lowe panel is fo cicula obits showing esults fom (6). The thee types of cuves in each panel coespond to the thee possibilities of the velocity dispesion pofile being geneated 1) by stas Figue 1. Uppe panel: pedicted line-of-sight velocity dispesion pofiles fo the Fonax dwaf in the case of isotopic obits. The dotted line shows the pediction obtained by solving Jeans equation fo stas only. The thee solid lines coespond to stas + DM with the dak halo masses of 10 8, 10 9 and M fom bottom to top. The thee dashed lines show esults fo stas in MOND with a 0 = 1, 3, cm/s fom bottom to top. The data fo Fonax ae fom Mateo (1997). Lowe panel: same as uppe panel, but fo cicula obits. moving unde Newtonian gavitational acceleation geneated by stas only (dotted cuves), ) by stas moving unde Newtonian gavitational acceleation geneated by stas and DM (solid lines), 3) by stas moving unde modified (MONDian) acceleation geneated by stas (dashed lines). Fo all cases the 3D distibution of stas ν() is given by (9). The gavitational acceleation fo the Newtonian cases 1) and ) is g = g N = GM()/ with M() = M () fo 1) and M() = M ()+M NFW() fo ), whee M () and M NFW() ae given by equations (10) and (13), espectively. In the case of MOND we have g = g M fom equation (15) withg N = GM ()/.Thetheesolidcuvesineachpanel coespond fom bottom to top to the viial masses of the

4 4 Ewa L. Lokas Figue. Best-fitting models fo Fonax based on stas+dm (solid line) and MOND (dashed line). The model involving DM has β = 1.4 and M v = M. The model with MOND equies β = 1 and a 0 = cm/s. dak halo M v = 10 8, 10 9 and M. The thee dashed lines show fom bottom to top esults obtained with diffeent chaacteistic MOND acceleation scales a 0 = 1, 3, cm/s. Figue 1 clealy shows that the stas alone cannot epoduce the obseved velocity dispesion pofile of Fonax, unless the stella mass-to-light atio is much lage than assumed. In the following we will fathe analyse the cases ) and 3) tying to find the best-fitting models. Compaing the pedictions shown in Figue 1 with the data we find that a cetain amount of tangential anisotopy is needed to epoduce the shape of the obseved velocity dispesion pofile. Indeed, tying moe adial anisotopy than in the case of isotopic obits by using small values of a (compaable to R S o v) in (3) and (5) makes the cuves in the uppe panel of Figue 1 decease even moe steeply with distance, contay to the tend obseved in the data. Theefoe in what follows we adopt the velocity anisotopy model of β =const< 0 and teat β as a fee paamete. To find the best models we pefomed a two-paamete fitting fo each case. The fitting paametes wee (M v,β) fo stas+dm case and (a 0,β) fo MOND. The esults of minimizing χ weighted by the vaiances of the data points ae shown in Figue. It tuns out that the DM model fits the data best fo β = 1.4 and M v = M. This viial mass is moe than two odes of magnitude lage than the stella mass ( M fo Υ V = 1M /L ). In the case of MOND the best-fitting paametes ae β = 1 and a 0 = cm/s. This value of a 0 is consistent with othe estimates (Milgom 1983b). Figue 3 shows the combined mass-to-light atio in V- band fo Fonax calculated fom M/L V = MNFW()+M (), (16) L V() whee M NFW() and M () ae given by equations (13) and Figue 3. The dependence of the mass-to-light atio of Fonax on distance. The mass comes fom stas and DM with the best fitting M v = M. (10) espectively, and L V() = M ()/Υ V is the luminosity distibution. The shape of the cuve eflects the fact that DM dominates both at small and lage distances. The divegence of M/L at the cente of the galaxy is due to the fact that the luminosity density (9) inceases towads the cente moe slowly than the DM cusp of 1 in (11). The gowth of M/L at lage distances is due to the behaviou of M NFW() which diveges logaithmically while the mass in stas is finite. The minimum of M/L V 11 appeas at a scale of the ode of the Sésic adius and is consistent with the values of the cental mass-to-light atio estimated by Mateo et al. (1991). 4 DISCUSSION We have found that the pedictions based on MOND with a 0 = cm/s fit the Fonax data equally well as models based on DM with viial mass of M v = M. In both cases a cetain amount of tangential anisotopy is needed to epoduce the shape of the obseved velocity dispesion pofile. It tuned out that the MOND pedictions in the case of Fonax dwaf need less tangential anisotopy than DM pedictions to epoduce the shape of the velocity dispesion pofile. It emains to be checked how the DM pedictions depend on the paticula shape of DM distibution. The paticula fitted values of a 0 and M v ae of couse sensitive to the obsevational paametes assumed. The measued luminosity, distance and stuctual paametes of the dwafs ae all known with limited accuacy, but the detailed analysis of the eos in them is beyond the scope of this pape. In the case of MOND the most doubtful input paamete is the stella mass-to-light atio in V-band. We have assumed hee Υ V 1M /L chaacteistic fo a elatively young stella population of Fonax, howeve this value may well be diffeent by a facto of two (Mateo et al. 1991). To estimate its effect we pefomed a simila MOND fitting as descibed in the pevious section with

5 Velocity dispesions of dwaf spheoidals 5 a seious poblem fo MOND (we agee on this point with Gehad & Spegel 199) since a 0 is supposed tobe aunivesal constant. Binging the value down by an ode of magnitude would equie vey high stella mass-to-light atio Υ V 10M /L o we would have to believe that the measuements seiously oveestimate the velocity dispesion. ACKNOWLEDGEMENTS Figue 4. Best-fitting models fo Daco (left panel) and Usa Mino (ight panel) based on stas+dm (solid lines) and MOND (dashed lines). The data ae fom Amandoff et al. (1997) Υ V = 0.5 and M /L, obtaining espectively a 0 = 4.3 and cm/s with the same best-fitting β. If Υ V = 0.5M /L is indeed the lowe limit of mass-to-light atio, then a 0 = cm/s is the uppe limit fo a 0. In the case of DM the stella mass-to-light is not likely to influence much the esults since the input fom stas to the velocity dispesion is small. Hee the most doubtful paamete is the concentation of DM haloes, which we adopted to follow the same dependence on mass as lage haloes actually studied in N-body simulations. Given, howeve, that this dependence is weak in the ΛCDM Univese, the possible eo associated with this facto is athe small. One moe eason to woy ae the tidal inteactions that may affect the dynamics of some dwafs and ou intepetations of thei velocity dispesions. Piatek & Pyo (1995) have pefomed numeical simulations of such effects and concluded that in Newtonian dynamics they should not affect much the infeed mass-to-light atios (but see Koupa 1997). Howeve, Bada & Milgom (000) found that in MOND, due to its non-lineaity, a galaxy system can be moe affected by an extenal field. Onemayaskifsimila esults holdfo othedsphgalaxies as well. Of those studied by IH the velocity dispesion pofiles with a few data points ae available only fo Sextans, Usa Mino and Daco dwafs (Hageaves et al. 1994a, 1994b, 1996; Amandoff, Olszewski & Pyo 1995; Amandoff, Pyo & Olszewski 1997). Applying the same analysis as we did fo Fonax to the 5 data points fo Daco and Usa Mino fom Amandoff et al. (1997) (with Υ V = 1M /L ) we find that again the MOND and DM models can fit the data equally well. Fo Daco the best-fitting models have tangential anisotopy, we get β = 0.9, M v = M ( 10 4 the mass in stas) fo DM and β = 1.5, a 0 = cm/s fo MOND. Fo Usa Mino the best fits equie a cetain amount of adial anisotopy which we model with (). We get a = v = 0.86R S = 0.16 kpc, M v = M (10 3 the mass in stas) fo DM and a = 0.01 v = 0.76R S = 0.14 kpc, a 0 = cm/s fo MOND. The pedictions of the best-fitting models togethe with the data ae shown in Figue 4. The case of Usa Mino is doubtful since it has been ecently obseved to possess tidal tails (Matínez-Delgado et al. 001) and theefoe may not be in viial equilibium. The high value of a 0 estimated fo Daco can, howeve, pose IwishtothankA.Gaham, P.Koupa,G.Mamon, C. Pyo, P. Salucci and the efeee, M. Mateo, fo useful comments and suggestions which helped to impove the pape. I am also gateful to Y. P. Jing fo poviding the esults fo the concentation of DM haloes in numeical fom. This eseach was patially suppoted by the Polish State Committee fo Scientific Reseach gant No. P03D0319. REFERENCES Amandoff T. E., Olszewski E. W., Pyo C., 1995, AJ, 110, 131 Amandoff T. E., Pyo C., Olszewski E. W., 1997, in Sofue Y., ed., Poc. IAU Symp. 184, The Cental Regions of the Galaxy and Galaxies. Kluwe, Dodecht, p. 35 Binney J., Mamon G. A., 198, MNRAS, 00, 361 Binney J., Temaine S., 1987, Galactic Dynamics. Pinceton Univ. Pess, Pinceton, chap Bada R., Milgom M., 000, ApJ, 541, 556 Caon N., Capaccioli M., D Onofio M., 1993, MNRAS, 65, 1013 Caldwell N., 1999, AJ, 118, 130 Ciotti L., 1991, A&A, 49, 99 Gehad O., Konawitte A., Saglia R. P., Bende R., 001, AJ, 11, 1936 Gehad O. E., Spegel D. N., 199, ApJ, 397, 38 Hageaves J. C., Gilmoe G., Iwin M. J., Cate D., 1994a, MN- RAS, 69, 957 Hageaves J. C., Gilmoe G., Iwin M. J., Cate D., 1994b, MNRAS, 71, 693 Hageaves J. C., Gilmoe G., Iwin M. J., Cate D., 1996, MN- RAS, 8, 305 Iwin M., Hatzidimitiou D., 1995, MNRAS, 77, 1354 Jing Y. P., Suto Y., 000, ApJ, 59, L69 Koupa P., 1997, NewA,, 139 Lima Neto G. B., Gebal D., Máquez I., 1999, MNRAS, 309, 481 Lokas E. L., Mamon G. A., 001, MNRAS, 31, 155 Matínez-Delgado D, Alonso-Gacía J., Apaicio A., Gómez- Flechoso M. A., 001, ApJ, 549, L63 Mateo M., 1997, in Anaboldi M. et al., eds, ASP Conf. Se. Vol. 116, The Natue of Elliptical Galaxies. Aston. Soc. Pac., San Fancisco, p. 59 Mateo M., 1998, in Richtle T., Baun J. M., eds, Poc. Bonn/Bochum-Gaduietenkolleg Wokshop, The Magellanic Clouds and Othe Dwaf Galaxies. Shake Velag, Aachen, p. 53 Mateo M., Olszewski E. W., Welch D. L., Fische P., Kunkel W., 1991, AJ, 10, 914 Meitt D., 1985, AJ, 90, 107 Milgom M., 1983a, ApJ, 70, 365 Milgom M., 1983b, ApJ, 70, 371 Milgom M., 1984, ApJ, 87, 571 Milgom M., 1995, ApJ, 455, 439 Navao J. F., Fenk C. S., White S. D. M., 1997, ApJ, 490, 493 Osipkov L. P., 1979, PAZh, 5, 77 Piatek S., Pyo C., 1995, AJ, 109, 1071 Pyo C., Komendy J., 1990, AJ, 100, 17 Sésic J. L., 1968, Atlas de Galaxies Austales, Obsevatoio Astonomico, Codoba Thomas P. A. et al., 1998, MNRAS, 96, 1061