A statisticalmechanical explanation of dark matter halo properties
Key Laboratory of Frontiers in Theoretical Physics, Institute of
Theoretical Physics, Chinese Academy of Science,
Beijing
100190,
PR China
email: billkang@itp.ac.cn
Received:
27
May
2010
Accepted:
4
December
2010
Context. Cosmological Nbody simulations have revealed many empirical relationships of dark matter halos, yet the physical origin of these halo properties still remains unclear. On the other hand, the attempts to establish the statistical mechanics for selfgravitating systems have encountered many formal difficulties, and little progress has been made for about fifty years.
Aims. The aim of this work is to strengthen the validity of the statisticalmechanical approach we have proposed previously to explain the dark matter halo properties.
Methods. By introducing an effective pressure instead of the radial pressure to construct the specific entropy, we use the entropy principle and proceed in a similar way as previously to obtain an entropy stationary equation.
Results. An equation of state for equilibrated dark halos is derived from this entropy stationary equation, by which the dark halo density profiles with finite mass can be obtained. We also derive the anisotropy parameter and pseudophasespace density profile. All these predictions agree well with numerical simulations in the outer regions of dark halos.
Conclusions. Our work provides further support to the idea that statistical mechanics for selfgravitating systems is a viable tool for investigation.
Key words: methods: analytical / dark matter / galaxies: statistics / hydrodynamics
© ESO, 2011
1. Introduction
The investigation of structures and properties of dark matter halos is one of the important open questions in modern cosmology. Cosmological Nbody simulations have revealed many empirical relationships concerning dark halos, such as the density profile (Navarro et al. 1997; Moore et al. 1999; Einasto 1965), the velocity dispersion, the anisotropy parameter, and the pseudophasespace density profile (Navarro et al. 2010; Ludlow et al. 2010). Yet the physical origin of these halo properties still remains unclear.
For about fifty years, many authors have attempted statisticalmechanical approaches to investigate the properties of selfgravitating systems (to name a few, Antonov 1962; LyndenBell 1967; Shu 1978; Tremaine et al. 1986; Sridhar 1987; Stiavelli & Bertin 1987; White & Narayan 1987; Tsallis 1988; Spergel & Hernquist 1992; Soker 1996; Kull et al. 1997; Nakamura 2000; Chavanis 2002; Hjorth & Williams 2010). However, these attempts have encountered many formal difficulties, and little progress has been made (Mo et al. 2010).
In a recent work (He & Kang 2010), we employed a phenomenological entropy form of ideal gas, first proposed by White & Narayan (1987), to revisit this question. First, subject to the usual mass and energyconservation constraints, we calculated the firstorder variation of the entropy, and obtained an entropy stationary equation. Then, incorporated in the Jeans equation, and by specifying some functional form for the anisotropy parameter, we solved the two equations numerically and demonstrated that the anisotropy parameter plays an important role in attaining a density profile that is finite in mass, energy, and spatial extent. If incorporated again with some empirical density profile from simulations, the theoretical predictions of the anisotropy parameter and the radial pseudophasespace density in the outer regions of the dark halos agree very well with the simulation data. Finally, we calculated the secondorder variation, which reveals the seemingly paradoxical but actually complementary consequence that the equilibrium state of selfgravitating systems is the global minimum entropy state for the whole system, but simultaneously the local maximum entropy state for every and any finite volume element of the system. Our findings indicate that the statisticalmechanical approach should be a viable tool to account for these empirical relationships of dark halos, and may provide crucial clues to the development of the statistical mechanics of selfgravitating systems as well as other longrange interaction systems.
Despite its great success, we should point out that the specific entropy in that work is not defined in a consistent way. On the one hand, as addressed above, the anisotropy parameter is of great significance to get a finite density profile, whereas on the other hand, we had to use only the radial pressure (or velocity dispersion) for the entropy definition. Any combinations of other directions’ velocity dispersions that included the anisotropy parameter, failed to produce correct equations, and hence were rejected. Besides, another shortcoming of that work should be pointed out, which is the agreements between all predictions and the simulations are restricted only to a narrow radial range of r/r_{2} > 2.
The aim of the current work is to remedy these defects. With the analogy of the isotropic Jeans equation, we introduce the effective pressure instead of the radial pressure to define the entropy. In a similar procedure, we derive a new entropy stationary equation, which can be easily integrated to attain the equation of state of the equilibrated dark halos, which He & Kang (2010) failed to obtain with the radial pressure. It is this state equation that we employ to derive the halo density profile. Additionally, with this improved treatment, the radial range of the agreements between the predictions and simulations is increased by one order of magnitude.
This paper is organized as follows. In Sect. 2 we define the effective pressure and derive the equation of state for the selfgravitating dark matter halos. In Sect. 3 we solve this equation approximately, and compare all results with the simulations of dark halos. We present the summary and conclusion in Sect. 4.
2. Basic theory and formulae
2.1. Isotropic Jeans equation and effective pressure
After experiencing the violent relaxation process (LyndenBell 1967), the selfgravitating system settles into a virial equilibrium state, which is usually called the quasistationary state (Campa et al. 2009), which is described by the Jeans equation in spherical coordinates^{1} as (Binney & Tremaine 2008): (1)where , Φ is the gravitational potential, is the anisotropy parameter, and and are the tangential and radial velocity dispersion, respectively.
For β = 0, i.e. the isotropic case, the Jeans equation is (2)where is the radial pressure. Equation (2) resembles the hydrostatic equilibrium equation of an ordinary fluid, ∇p = −ρ∇Φ, in that the radial pressure gradient, the buoyancy, resists the gravitational force of the whole system. Enlightened by this similarity, we generalize Eq. (2) to the anisotropic case by defining the effective pressure P through: (3)Then, the Jeans equation, Eq. (1), can be rewritten in terms of the effective pressure P as (4)Below we will use the effective pressure to construct the specific entropy form.
2.2. Equation of state
As mentioned in the introduction, in He & Kang (2010) we employed the entropy principle and derived an entropy stationary equation by introducing the specific entropy, (5)This entropy form was first used by White & Narayan (1987) for the isotropic selfbounded gas sphere, but we assumed it can also be applied to the velocityanisotropic case. This is not a selfconsistent treatment, and also it is hard to generalize to include the anisotropy parameter.
Based on the above similarity between Eqs. (2) and (4), we replace p_{r} in Eq. (5) by the effective pressure P to make the entropy more appropriate for the anisotropic case, so that the total entropy is (6)where s′ = ln(P^{3/2}ρ^{−5/2}). Subject to the constraints of mass and energy conservation, we calculate the firstorder variation of the entropy in a similar process as the one in He & Kang (2010), to obtain the new entropy stationary equation as (7)where λ is the Lagrangian multiplier corresponding to the energy conservation. The lefthand side of this equation is formally the same as that of He & Kang (2010), except that p_{r} is replaced by the effective pressure P.
Equation (7) can be readily transformed into (8)and can be directly solved as (9)where μ is an integration constant. μ and λ are related to the total mass and energy of the dark halo, that is, they can be specified by the total mass M, and total energy E, as μ = μ(M,E) and λ = λ(M,E), and vice versa.
Equation (9) describes the relationship between ρ and P, which is exactly the equation of state of the equilibrated system, but P is also related to ρ through the differential equation, Eq. (4), thus Eq. (9) cannot directly provide us with the dark halo density profile. Below we will explore how to obtain the density profile from this state equation.
3. Results
From Eq. (9) we can see that ρ scales with P as ρ ~ λP at the center of the dark halo, but ρ ~ μP^{3/5} at the outskirts, since both ρ and P are large at the center, but small at the outskirts of the dark halos. This suggests that it would be advantageous to analyze the approximate solutions of these two cases, before we discuss the general solution of Eq. (9).
Fig. 1 Isothermal density profiles. The inner slope of the density profile is dependent on different λ. The Einasto profile, ln(ρ/ρ_{2}) = −2/α((r/r_{2})^{α} − 1), with α ≈ 0.17 (see Navarro et al. 2010), is also indicated for comparison. 

Open with DEXTER 
3.1. λ ≠ 0,μ = 0
As explained above, this case is the approximation of Eq. (9) at the halo center. If λ is renamed as (we will see λ must be positive in the following)^{2}, then with μ = 0, Eq. (9) is reexpressed as (10)Incorporating Eq. (4), and differentiating both sides of the equation with respect to r, we have (11)with the solution as , which resembles the equation that the isothermal gas satisfies (see Binney & Tremaine 2008, p. 303). Differentiating Eq. (11) again, we have (12)We show the solutions in Fig. 1, with the case of λ = 2 corresponding to the singular isothermal sphere^{3}, ρ ~ r^{2}. We can also see that λ with λ < 2 bends the singular isothermal solution toward a centrallycored density profile, and the size of the core depends on the value of λ, but the slope at the outskirts of the dark halos remains unchanged.
Fig. 2 Polytropic density profiles. Panels a) and b) correspond to the two approximate inner powerlaw solutions, ρ ∝ r^{n}, with the power index as n = −1.5 and n = 0, respectively. All density profiles are selftruncated, with the truncation radii varying with different μ. 

Open with DEXTER 
3.2. λ = 0,μ ≠ 0
In this case, the equation of state, Eq. (9), is reduced to (13)which resembles a polytropic process of a common gas. Differentiating this equation with respect to r, and incorporating Eq. (4), we get the following differential equation (14)As analyzed previously, this equation is just an approximation of Eq. (9) at the outskirts of the dark halos. As r reaches zero, we can see that the righthand side of the equation also reaches zero, so that we can approximate the solution with a powerlaw form, ρ ~ r^{n}, with the powerindex n = 0,−1.5. Then starting from the two powerlaw inner solutions at sufficiently small r, we numerically evaluate Eq. (14) with different μ, and show all results in Fig. 2. Since Eq. (14) is only valid at the outskirts of dark halos, we can see that regardless of the nonuniqueness of the inner slopes at small radii, all resulting density profiles behave as if they were truncated by the μterm of this equation, so that both the mass and the energy of the dark halos are not infinite.
Fig. 3 a) The density profiles and b) the logarithmic slopes corresponding to the density profiles in panel a). The profile and the slope derived from Eq. (15) are indicated as “solution” in the figure, with λ = 2.27,μ = 0.36. The major fitting formulae of the density profile, NFW, Moore, and Einasto profile, are also shown for comparison. 

Open with DEXTER 
3.3. λ ≠ 0,μ ≠ 0
With the heuristic results from the above two subsections, we explore the complete solution of this general case. It may not be easy to solve Eqs. (4) and (9) exactly, so we try to find the approximate solution to the density profile. We treat Eq. (10) as a first approximation and substitute it into the derivative of the state equation, Eq. (9), then we have (15)in which again is the renaming of λ, as . We numerically evaluate this equation, and the resulting density profile is exhibited in Fig. 3a. We also indicate in Fig. 3b the logarithmic slope of the density profile, i.e., γ = −dlnρ/dlnr. In the inner region ρ indeed behaves as the isothermal solution of Eq. (12). In the outer region, the density profile does not decline as fast as the polytropic solution of Eq. (14) because we only consider the first approximation, but it is still obvious that the whole density profile is roughly the superposition of the above two subsections’ solutions.
We can see that with an appropriate value of λ and μ, both the resulting density profile and its logslope agree very well with the fitting formulae at the radii r/r_{2} ≥ 0.2, the right part of the vertical dotted line in Fig. 3. At small radii with r/r_{2} < 0.2, the solution is qualitatively accepted, but the agreement is not as good as the case at large radii.
Fig. 4 The radial velocity dispersion , the anisotropy parameter β(r), and the pseudophasespace density Q(r) are shown in panels a), b), c), respectively. The results based on the resulting density profile in Sect. 3.3 are indicated as “solution” in the figure. The results based on the Einasto profile are also shown for comparison. The vertical scales of panels a) and c) are arbitrary. 

Open with DEXTER 
3.4. , β, and profiles
If the density profile ρ(r) is determined, we can obtain the effective pressure P by integrating Eq. (4). Numerical simulations indicate that the anisotropy parameter, β, is usually a small number (Navarro et al. 2010). From Eq. (3) we can see that when β is small, we can approximate the effective pressure P as (16)by which can be expressed as .
From Eqs. (1) and (15) we obtain (17)Then if , ρ, and hence the logarithmic slope of the density profile γ = −dlnρ/dlnr are given, β can be expressed as (18)where , .
We present a simple analysis of what Eq. (18) implies. In the inner region of dark halos ρ is large, the term including μ can be neglected, and from the simulation results (Navarro et al. 2010) we know both a(r) and b(r) vary slowly with r, so it seems that there is an approximately but not strictly linear relationship between β and γ. The existence of this linear relation has been observed in both Hansen & Moore (2006) and the inner region of Navarro et al. (2010). Because this analysis will not be valid for regions at larger radii, the linear relationship should not exist for the whole region. Indeed, the latest simulations by Navarro et al. (2010) and Ludlow et al. (2010) show that β(r) is not the monotonic function of r, that is, haloes are nearly isotropic near the center, radially biased to the maximum at some radius and approximately isotropic again in the outskirts, which is against the earlier monotonic result of Hansen & Moore (2006).
In the inner region, since a < 1 and b > 0, so from Eq. (18) we must have γ > 2β, and in the outer region, since γ > 2 and β < 1, we always have γ > 2β. So we have γ > 2β in the whole region, which is consistent with the result of An & Evans (2005) and Ciotti & Morganti (2010).
Next, we analyze the behavior of the pseudophasespace density (or ). From Eq. (18) we have (19)In the inner region we know from Navarro et al. (2010) that β, a, and b all vary slowly with r, so that dlnQ/dlnr roughly seems to be a constant. If we take Hansen & Moore (2006)’s result as a reference, a ≈ 0.4, b ≈ 0.32, and assume β ≈ 0.12 on the average, then dlnQ/dlnr ≈ −1.88, very close to the spherical secondaryinfall similarity solution (Bertschinger 1985).
With the density profile derived in the previous subsection (see Fig. 3a), and by Eq. (16), we calculate the approximate radial velocity dispersion, , which is shown in Fig. 4a. The velocity dispersion calculated with Einasto profile is also indicated in the figure for comparison. We can see that both the two results roughly agree with the simulation results of Navarro et al. (2010) when r/r_{2} > 0.2, as is consistent with the accuracy of the resulting density profile of Sect. 3.3. The discrepancy is significant when r/r_{2} < 0.2, which may suggest that there may be still some crucial physics that we have failed to capture in the current formulation.
Given the density profile ρ(r), and the velocity dispersion , and with the two parameters λ (or ) and μ being the same as we computed the density profile, we can also obtain β(r) via Eq. (18) as well as the pseudophasespace density Q(r). The results are exhibited in Figs. 4b and 4c, respectively. The corresponding results derived by using the Einasto profile are also indicated for comparison. From Fig. 4b, we notice that the nonmonotonicity of β(r) is once again reproduced in our current work, which agrees well with the simulation results of Navarro et al. (2010) and Ludlow et al. (2010). However, we find that both the position and the height of the peak of the β(r) profile deviate slightly from the simulations. We speculate that this discrepancy may be originating in the approximation we adopted for Eqs. (15)–(17). The result based on the Einasto profile cannot be accepted for r/r_{2} < 0.2, as is consistent with the previous results for both ρ and .
As expected, the resulting Q(r) profile closely follows a power law, r^{1.875}, near the center of the dark halo. However, the predictions do show a manifest curveup deviation from this power law at the outskirts of dark halos, which has been observed by the simulation result of Ludlow et al. (2010), and has already been reproduced by our previous study (He & Kang 2010).
4. Summary and conclusion
Cosmological Nbody simulations have revealed many empirical relationships concerning dark halos. Up to now, however, there are no robust physical origins for these empirical relationships. In this work, we aim to provide a unified statisticalmechanical explanation to the relationships concerning the matter density, anisotropy parameter and the pseudophasespace density profile. The main steps of this work are similar to those of He & Kang (2010), i.e., by using the entropy principle with mass and energyconservation as constraints to derive the entropy stationary equation, Eq. (7). Different from He & Kang (2010), however, based on the analogy with the isotropic Jeans equation, we phenomenologically introduce the effective pressure P instead of the radial pressure p_{r} to construct the specific entropy. The effective pressure is a key quantity in this work and proves to be more powerful than p_{r}.
The entropy stationary equation can be easily integrated to attain Eq. (9), the equation of state of the equilibrated dark halos, which He & Kang (2010) failed to attain with the radial pressure p_{r}. It is this equation of state that we employ to derive the halo density profile ρ(r). Notice that this density profile has nothing to do with β.
First, we study the approximate solutions of Eq. (9) in two asymptotic cases, i.e. the inner and outer region of the dark halos. The inner solution is the isothermal density profile, with the inner logarithmic slope varying with different parameter λ. The outer solution is the polytropic density profile, selftruncated at sufficiently large radii with a different parameter μ, so that the
density profile has finite mass, energy, and extent. The approximate complete solution can be roughly regarded as the superposition of the two asymptotic solutions. By choosing appropriate values of λ and μ, we can reproduce both the density profile and its logarithmic slope quite well in the region r/r_{2} > 0.2. In the inner region of r/r_{2} < 0.2, the solution seems to be qualitatively accepted, but the agreement between the solution and the fitting formulae is not quantitatively well. We speculate that there may be still some crucial physics that we have failed to capture in the current formulation.
Then, with the resulting density profile, we compute the radial velocity dispersion , and then β(r) and the pseudophasespace density Q(r) profile. We find that all these predictions agree quite well with the corresponding results of the latest simulations of Navarro et al. (2010) and Ludlow et al. (2010) in the region of r/r_{2} > 0.2. Despite the slight discrepancies from the simulation results in the inner regions, these agreements at the outer regions of dark halos strengthen the conclusions of He & Kang (2010), indicating once again the great success of the statisticalmechanical approaches for the selfgravitating systems.
Acknowledgments
D.B.K. is very grateful for the comments and suggestions of the anonymous referee. This work is supported by the National Basic Research Programm of China, No. 2010CB832805.
References
 An, J., & Evans, N. W. 2005, A&A, 444, 45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Antonov, V. A. 1962, Vestnik Leningrad Univ., 7, 135; English translation: Antonov, V. A. 1985, in Dynamics of Star Clusters, ed. J. Goodman, & P. Hut (Dordrecht: Reidel), 525 [Google Scholar]
 Bertschinger, E. 1985, ApJS, 58, 39 [NASA ADS] [CrossRef] [Google Scholar]
 Binney, J., & Tremaine, S. 2008, Galactic Dynamics, 2nd edn. (Princeton, New Jersay: Princeton University Press) [Google Scholar]
 Campa, A., Dauxois, T., & Ruffo, S. 2009, Phys. Rep., 480, 57 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Chavanis, P. H. 2002, in Proceedings of the Conference on Multiscale Problems in Science and Technology ed. N. Antonic, C. J. van Duijn, W. Jager, & A. Mikelic (Berlin: SpringerVerlag), 85 [Google Scholar]
 Ciotti, L., & Morganti, L. 2010, AIPC, 1242, 300 [NASA ADS] [Google Scholar]
 Einasto, J. 1965, Trudy Inst. Astrofizicheskogo AlmaAta, 51, 87 [NASA ADS] [Google Scholar]
 Hansen, S. H., & Moore, B. 2006, New A, 11, 333 [NASA ADS] [CrossRef] [Google Scholar]
 He, P., & Kang, D. B. 2010, MNRAS, 406, 2678 [NASA ADS] [CrossRef] [Google Scholar]
 Hjorth, J., & Williams, L. L. R. 2010, ApJ, 722, 851 [NASA ADS] [CrossRef] [Google Scholar]
 Kull, A., Treumann, R. A., & Böhringer, H. 1997, ApJ, 484, 58 [NASA ADS] [CrossRef] [Google Scholar]
 Ludlow, J. F., Springel, V., Vogelsberger, M., et al. 2010, MNRAS, 406, 137 [NASA ADS] [CrossRef] [Google Scholar]
 LyndenBell, D. 1967, MNRAS, 136, 101 [NASA ADS] [CrossRef] [Google Scholar]
 Mo, H. J., van den Bosch, F., & White, S. D. M. 2010, Galaxy Formation and Evolution (New York: Cambridge University Press) [Google Scholar]
 Moore, B., Quinn, T., Governato, F., Stadel, J., & Lake, G. 1999, MNRAS, 310, 1147 [NASA ADS] [CrossRef] [Google Scholar]
 Nakamura, T. K. 2000, ApJ, 531, 739 [NASA ADS] [CrossRef] [Google Scholar]
 Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 [NASA ADS] [CrossRef] [Google Scholar]
 Navarro, J. F., Ludlow, A., Springel, V., et al. 2010, MNRAS, 402, 21 [NASA ADS] [CrossRef] [Google Scholar]
 Shu, F. H. 1978, ApJ, 225, 83 [NASA ADS] [CrossRef] [Google Scholar]
 Soker, N. 1996, ApJ, 457, 287 [NASA ADS] [CrossRef] [Google Scholar]
 Spergel, D. N., & Hernquist, L. 1992, ApJ, 397, L75 [NASA ADS] [CrossRef] [Google Scholar]
 Sridhar, S. 1987, JA&A, 8, 257 [NASA ADS] [Google Scholar]
 Stiavelli, M., & Bertin, G. 1987, MNRAS, 229, 61 [NASA ADS] [Google Scholar]
 Tremaine, S., Hénon, M., & LyndenBell, D. 1986, MNRAS, 219, 285 [NASA ADS] [CrossRef] [Google Scholar]
 Tsallis, C. 1988, J. Stat. Phys., 52, 479 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 White, S. D. M., & Narayan, R. 1987, MNRAS, 229, 103 [NASA ADS] [Google Scholar]
All Figures
Fig. 1 Isothermal density profiles. The inner slope of the density profile is dependent on different λ. The Einasto profile, ln(ρ/ρ_{2}) = −2/α((r/r_{2})^{α} − 1), with α ≈ 0.17 (see Navarro et al. 2010), is also indicated for comparison. 

Open with DEXTER  
In the text 
Fig. 2 Polytropic density profiles. Panels a) and b) correspond to the two approximate inner powerlaw solutions, ρ ∝ r^{n}, with the power index as n = −1.5 and n = 0, respectively. All density profiles are selftruncated, with the truncation radii varying with different μ. 

Open with DEXTER  
In the text 
Fig. 3 a) The density profiles and b) the logarithmic slopes corresponding to the density profiles in panel a). The profile and the slope derived from Eq. (15) are indicated as “solution” in the figure, with λ = 2.27,μ = 0.36. The major fitting formulae of the density profile, NFW, Moore, and Einasto profile, are also shown for comparison. 

Open with DEXTER  
In the text 
Fig. 4 The radial velocity dispersion , the anisotropy parameter β(r), and the pseudophasespace density Q(r) are shown in panels a), b), c), respectively. The results based on the resulting density profile in Sect. 3.3 are indicated as “solution” in the figure. The results based on the Einasto profile are also shown for comparison. The vertical scales of panels a) and c) are arbitrary. 

Open with DEXTER  
In the text 