Searching the Parameters of Dark Matter Halos on the Basis of Dwarf Galaxies’ Dynamics

L. M. Chechin^1,2, T. K. Konysbayev²

¹V.G. Fessenkov Astrophysical Institute, Almaty, Kazakhstan.
²Al-Farabi Kazakh National University, Almaty, Kazakhstan.

DOI: 10.4236/jmp.2016.79089                

Abstract

Article devoted to searching the parameters of dark matter halos on the base of dwarf galaxies’ dynamics (Messier 32 and Leo I). For doing this, we propose the new approach founded on construction the coupled elliptical trajectory for a probe body in the gravitational fields of Newtonian potential and potential of dark matter’s halo. This allows more accuracy estimate its central density for the Navarro-Frenk-White profile and free parameter for the Einasto profile . Our result is in good correlation with results of other authors that are got by different numerical methods.

Keywords

Dwarf Galaxies, Dark Matter Halo, Central Density of Dark Matter Profile, Parameters of Dark Matter’s Halo

Received 11 March 2016; accepted 24 May 2016; published 27 May 2016

1. Introduction

One of the urgent problems of modern cosmology is searching the dark matter properties. Dark matter―an unusual type of cosmic substance, which is in the overall energy balance of the Universe approximately equals 27% [1] . Astronomical observations suggest that dark matter is concentrated around the large space objects such as galaxies and their clusters. Dark matter forms a halo which mass enriches about 90% of the galaxy total mass [2] . (Note, that dark matter may form different types of clumps without existence of any baryonic objects [1] ).

It should be noted that distribution of dark matter in the halo of galaxies is not uniform―it concentrates in their centers and decreases to peripheries. The corresponding distribution function of dark matter (profile) is usually founded on numerical methods that are modeling the dynamics of stars in galaxies. Today a number of profiles are known [3] - [7] , which include the unknown dark matter density in the centers of galaxies, as well as a number of free parameters.

Finding numerical values of these quantities is one of unsolved cosmological problems that from our viewpoint, can be partially solved by examining the dynamics of dwarf galaxies. Note that influence of dark energy and dark matter on galaxies’ dynamics successfully described even in the framework of Newtonian approach [2] [8] [9] . Therefore, we’ll be used it in our searching for estimation the magnitudes of the dark matter halo’s parameters.

2. On the Dynamics of Test Bodies in the Dark Matter Halo’s of Massive Galaxies

Now consider the following dynamic model―in the gravitational field of a massive galaxy, surrounded by a halo of dark matter, the probe body moves. The prototype of this model, for example, is the model of dwarf galaxy Messier 32 motion in the Andromeda galaxy or the motion of dwarf galaxy Leo I in the Milky Way. As for the distribution function of dark matter we’ll choose two most known of them―Navarro-Frenk-White profile and Einasto profile.

2.1 Navarro-Frenk-White Profile of Dark Matter

The Navarro-Frenk-White profile looks as follows [10]

(1)

Here and hereinafter―the density of dark matter in the center of galaxy,―the size of its halo,― the current radius. Assume that following relationship exists between them. Therefore, the profile function can be expanded in Taylor series with respect to above introduced small parameter. With the first order on accuracy we get

(2)

Then the potential energy of the dark matter of the field can be written down as

(3)

where coefficients in standard designations are

(4)

Later on we’ll use first term in (3) only because such potential energy leads to a closed trajectory that corresponds to the real observations of dwarf galaxies’ movement. Consequently, based on the conservation laws of energy and momentum for the test body dynamics, its trajectory can be write down in the standard form

(5)

(6)

(7)

From expressions (5)-(7) sees that it describes the elliptical trajectory that is similar to standard trajectory for one body in Newtonian mechanics [11]

(8)

(9)

where―constant of Newtonian potential for a gravitating mass. Marking and it’s clear that expressions (5) and (8) can be written down in the form of coupled elliptical trajectory

(10)

Now require that trajectories (5) and (8) coincide each other at the real movement of probe body. To substantiate this declaration we use the graphs of rotation curves for the Newtonian movement and for the movement of body in the gravitation field of dark matter. It is known that for the first case such curve is the hyperbola, while for the second case―the quasi-logarithmic line (for example, see [12] [13] ). Observations have shown that these curves cross at the critical point and the corresponding mean velocities at it approximately

equals. (Really, for the dwarf galaxy Messier 32 the observable velocities interval is [14] ; for the dwarf galaxy Leo I― [15] ). Because of that we’ll use these values for estimation of unknown dark matter halo’s parameters.

Here it’s necessary to point out that later on we’ll consider distances no larger than. At these distances curved line in the gravitation field of dark matter approximately looks like the straight line type of Note that such relation also follows from the equality centrifugal force and gravitational force of dark matter that relates to the first term in potential energy (3).

Put that angles of trajectories are coincide also, i.e. Hence, the following equality takes place

(11)

Since the total energy is larger than kinetic one and it, in its turn, is greater than potential energy , we can assume that Then expression (11) simplifies―

(12)

Now from (12) can be found the value of an unknown factor in the first term of expression (3)

(13)

Comparison (4) and (13) allows get the expression of central density of dark matter’s halos

(14)

For numerical estimations assume that the probe mass m in (14) equals to one. Thus roughly put,. After substitution the above mentioned numerical values for the position of probe body and its velocity, we obtain the following estimate of the central density of dark matter’s halos

(15)

It is interesting to compare this result with the previously obtained similar values. For example, in [16] . However, in [17] it was shown that the central part of dark matter density should not be larger than Furthermore, in [3] [18] [19] it is shown that the central density of dark matter . Analysis of these estimates gives that the most acceptable central density of the dark matter’s halo, probably, lays within the interval.

2.2. Einasto Profile for Dark Matter

Note that Navarro-Frank-White profile is completely identified. Nevertheless, in literature are known profiles of dark matter that contain from one [20] up to several [4] free parameters. In our article we’ll consider Einasto profile [20]

(16)

with one unknown parameter. Here as usual―the density of dark matter in the center of galaxy,―the size of its halo,―the current radius. Based on the previously adopted relationship between size of the halo and the current radius, the exponential function of Einasto profile may be expanded in Taylor series. Restricting ourselves by two terms only we get

(17)

In doing this the new additional condition we have used here―the smallness of the power parameter. Then Einasto profile in the approximate form is written down as

(18)

where

(19)

The potential energy of dark matter field in this case is as follows

(20)

where

(21)

Now for further calculations we use only the second term in (22), because it leads to a closed trajectory as before. Therefore, using the conservation laws of energy and momentum, its trajectory can be written as

(22)

(23)

(24)

From the expression (24) sees that it describes an elliptical trajectory that is similar to the standard trajectory in Newtonian mechanics. Repeating the previous arguments about the procedure of analyzing trajectories and their shapes, we find the relation

(25)

Hence, the expression of an unknown coefficient is

(26)

From the comparison of expressions (21), (26) and usage (19) we get

. (27)

Important to notice that expression (27), besides dynamic characteristics of a particle E and M, contains the current radius r and unknown parameter. Basing on our previous article [21] , it is possible to estimate and this parameter.

We write down the expression of central density of dark matter’s halo in the form

(28)

Here the coefficient X has various numerical values which are known in literature. In fact, from [3] it follows that from [18] ―from [19] ―from [21] ―(see also (15)). As for the central density of dark matter’s halo in Einasto profile, we write down it as

(29)

where Y is unknown coefficient. Substituting the necessary numerical values into (27) we find―for the dwarf galaxy Messier 32 for the dwarf galaxy Leo I―

From comparison of (27), (28) and (29) follows the expression of unknown parameter

(30)

Using the found values X and Y we can obtain this parameter. As for the galaxy Messier 32 and above given values of X, we find: And finally, for the Galaxy Leo I when we have: (At calculations we took into account that for such galaxy).

From the saying above we conclude―free parameter in Einasto profile must satisfy the following interval

(31)

Note that in article [1] the corresponding interval is closed to (31) and equals Also important to emphasize that in [22] the next interval of parameter have been argued

3. Conclusions

In the article, the new method for search characteristics of dark matter’s halo is proposed. Searching the Navarro-Frenk-White profile shows that the central part of density of dark matter’s halos must satisfy the following

magnitude This estimation allowed find the interval for free parameter of Einasto profile

Remarkable that it is close to similar intervals given in recent articles [22] - [24] , especially in [1] . That is why expedient to compare few methods of the Einasto parameter finding.

Authors of [22] used the semidegenerate thermal self-gravitating general relativistic fermionic gas as the dark matter model. As the result they got the dark matter profile in the form of broken line of four regions. First of them describes by rotation curve type of (the correctness of which on small distances we point out before), fourth of them―by the quasi-logarithmic line (that is correct for the large space scales). The peculiarity of such profile consists in existence of two intermediate regions for the moderate distances.

For describing this broken line in terms of Einasto profile (and not only) they used results of the survey THINGS were have been obtained highest quality rotational curves for 34 nearby spiral and irregular galaxies. Naturally that these curves depend not only the dark matter, but on some baryonic parts (baryonic substrate, relativistic gas, etc.) and their quantitative relation also. Therefore authors got the cited above range of Einasto parameter

In the framework of N-body simulations―“Aquarius” project with number of particles is 4.4 billions―au- thors [23] used Einasto profile in the form. Approximately their result is and looks like the quasi-permanent parameter. Nevertheless they pointed out that the magnitude of Einasto parameter slightly varies from halo to halo because different halos cannot, in general, be identical. Therefore power parameter must have the more width numerical interval.

Really, in [24] in the framework of “Aquarius” project also, but for one order smaller number of particles, was shown that Einasto parameter lays within the interval. In doing this authors emphasize that this parameter changes together with growth of galaxy’s mass. They also point out that if the shape of Einasto profile tends to isotherml one, while for Gaussian profile. So, the sought-for interval in reality must be.

From this examination, it’s clear the difference between ours models―in our case model based on the dwarf galaxies Newtonian dynamics and allowed get the improved estimate of dark matter core. Moreover, it’s very simple and allowed get the very plausible interval for the Einasto parameter. In contrary, the model in [22] is relativistic one, it’s based on the precisian curve lines observations but have the essentially width interval for Einasto parameter.

Two above mentioned other articles [23] [24] , as sees, are based on the N-body simulations method and give the closed to our power parameter. Together with that these articles devoted to study other unsolved problems―mass, velocity dispersion and anisotropy profiles of cold dark matter halos and its evolution, while they are not been under our consideration.

Therefore all of these models must be regard as additional each other which give rather satisfactory description of dark matter halo on different cosmological scales and for different space objects.

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	Berezinsky, B.S., Dokuchaev, B.I. and Eroshenko, Y.N. (2014) Advances in Physical Sciences, 184, 3-42. (In Russian)
[2]	Dolgachev, V.P., Domozhilova, L.M. and Chernin, A.D. (2003) Astronomy Reports, 47, 728-732.
[3]	Iocco, F., Pato, M., Bertone, G. and Jetzer, P. (2011) Journal of Cosmology and Astroparticle Physics, 11, 29.
[4]	Avila-Reese, V., Firmani, C., Klypin, A. and Kravtsov, A.V. (1999) Monthly Notices of the Royal Astronomical Society, 310, 527-539.
[5]	Burkert, A. (1995) The Astrophysical Journal Letters, 447, 171-175.
[6]	Catena, R. and Ullio, P. (2010) Journal of Cosmology and Astroparticle Physics, 8, 1-21.
[7]	Evans, N.W. and An, J. (2006) Physical Review D, 73, Article ID: 023524.
[8]	Kahn, F.D. and Woljter, L. (1959) The Astrophysical Journal, 130, 705-717.
[9]	Chechin, L.M. and Ibraimova, A.T. (2015) Astronomical and Astrophysical Transactions, in Press.
[10]	Navarro, J.F., Frenk, C.S. and White, S.D.M. (1996) The Astrophysical Journal, 462, 563.
[11]	Landau, L.D. and Lifschitz, E.M. (1988) Mechanics. M. Nauka.
[12]	Kirilov, A.A. and Turaev, D. (2006) Monthly Notices of the Royal Astronomical Society, 371, L31-L35.
[13]	Kravtsov, A.V. and Klypin, A.A. (1998) The Astrophysical Journal, 502, 48.
[14]	Howley, K.M., Guhathakurta, P., van der Marel, R., Geha, M., Kalirai, J., Yniguez, B., Kirby, E., Cuillandre, J.-C. and Gilbert, K. (2013) The Astrophysical Journal, 765, 22.
[15]	Lokas, E.L., Klimentowski, J., Kazantzidis, S. and Mayer, L. (2008) Monthly Notices of the Royal Astronomical Society, 390, 625-634.
[16]	Hideyoshi, A. (2010) Advances in Space Research, 45, 1007-1014.
[17]	Chechin, L.M. (2013) International Journal of Astronomy and Astrophysics, 3, 285-290.
[18]	Nesti, F. and Salucci, P. (2012) The Local Dark Matter Density. VIII International Workshop on the Dark Side of the Universe, Rio de Janeiro, 10-15 June 2012.
[19]	Zacek, V. (2007) Dark Matter. Fundamental Interactions: Proceedings of the 22nd Lake Louise Winter Institute, World Scientific Publishing Co. Pte. Ltd., 170-206.
[20]	Einasto, J. (2000) Astronomical Soc. Pacific Conference Series, 252, 85.
[21]	Chechin, L.M. and Konysbayev, T.K. (2016) Reports of the National Academy of Sciences of Kazakhstan. (In Russian, in Press)
[22]	Siutsou, I., Argüelles, C.R. and Ruffini, R. (2015) Astronomy Reports, 59, 656-666.
[23]	Navarro, J.F., Ludlow, A., Springel, V., et al. (2009) MNRAS, 402, 21-24.
[24]	Dutton, A.A. and Andrea Macciò, V. (2014) MNRAS, 441, 3359-3374.