_{1}

^{*}

It is well-known that equilibrium density profiles of gases and fluids in gravitational potentials have an
*r*
^{−1} dependence, where r is the radius. This is found in both astronomical observations and detailed simulations in spherically-symmetric geometries. It is also well-known that the standard equation for hydrostatic equilibrium does not produce such solutions. This paper utilizes a Lagrangian formulation that produces a closed-form
*r*
^{−1} solution and identifies the needed terms to be added to the standard equation for hydrostatic equilibrium to obtain such a solution. Variants of the
*r*
^{−1} solution avoid a density singularity at the origin and a total-enclosed mass singularity at infinity. The resulting solutions are shown to be in good agreement with well-established density profiles of ordinary matter in galaxies, dark matter haloes, and the atmosphere of earth. Comparisons are made to solutions based on the standard hydrostatic equation, including solutions based on the Lane-Emden equation. The origins of differences are explained.

Early numerical solutions for dark matter (DM) density profiles were found to have dependence inversely proportional to radius [^{−1} profile, where r is the radius from the nominal center of the distribution, in a region surrounding the origin [^{−1} profile is “ubiquitous” [

To analyze the spatial distribution of a gas of particles, N-body simulations [

( 1 / ρ ) ( d P / d r ) = − m v M e n c ( r ) G / r 2 . (1)

here ρ is the number density of DM, P is the pressure, r is the radius, m υ is the mass of a DM particle, M e n c ( r ) is the enclosed mass, and G is the gravitational constant. This equation can be solved using the well-known Lane-Emden formulation [

P = c γ ρ γ , (2)

where c γ is a constant for a given polytropic exponent γ . An inhomogeneous form of the Lane-Emden equation can be used when ordinary matter (OM) is present as shown later in the paper. As shown in Section 4, calculations using Equations (1) and (2) are inconsistent with an r^{−1} solution. Such solutions have no cusp at the origin and are a poor match to the standard Sersic or Einasto profiles based on observations [

To address this, a generalization of Equation (1) is derived using a Lagrangian formulation. This derivation is given in Section 2. A partial validation of this derivation is given in Section 3 using the properties of a well-known density profile of a gas in a gravitational potential, namely, the earth’s atmosphere. Section 4 presents behavior of solutions near the origin. Section 5 presents the behavior of solutions far from the origin and presents a sample calculation for a dark-matter halo. Section 6 compares the solutions to Einasto and de-projected Sersic profiles. Section 7 compares the solutions to results from the Lane Emden equation. Section 8 provides a brief summary.

The assumed kinetic and potential energy densities for a medium in a spherically symmetric geometry under the influence of gravity is given by

K . E . ( r ) = P ( r ) , P . E . ( r ) = − ρ ( r ) m v M e n c ( r ) G / r , (3)

where the variables are defined as in Equation (1). In the topic of this paper, M e n c ( r ) consists of both ordinary and dark matter. The ordinary (radiant) matter has mass density denoted ρ m , r a d ( r ) . Note that the pressure is treated as a kinetic energy for the obvious reason. This will allow a variational approach with a Lagrangian density. Explicit relativistic effects and the possible impact of angular momentum are neglected in this initial analysis

The pressure is assumed to be a function of density only. As an example, the number density for fermions in Equation (3) assuming thermodynamic equilibrium is

ρ = n s / ( 2 π 2 ℏ 3 ) ∫ p 2 d p / [ exp ( { [ ( p c ) 2 + ( m v c 2 ) 2 ] 1 / 2 − m v c 2 − μ v } / k T ) + 1 ] , (4)

where p is the fermion momentum, m v is its mass, n s is the number of spin states, μ v is the chemical potential, and ћ is Planck’s constant. The other variables have been defined earlier. In Equation (4), the chemical potential is set to zero for non-interacting particles. With Equation (4) the pressure for fermions has the usual form in the non-relativistic limit, and is given by

P = c 5 / 3 ρ 5 / 3 ， (5)

where c 5 / 3 = 1.914 ℏ 2 / m v assuming 2 spin states for a fermion. For the relativistic case, the exponent is 4/3, and c 4 / 3 = 5.536 ℏ c . For the general case of arbitrary polytropic exponent γ and assuming c γ is of the form α ћ a c b m ν c with α a dimensionless constant of order unity, one obtains c γ = α ћ 3 ( γ − 1 ) c 5 − 3 γ m ν 4 − 3 γ . With this relation and m v = 1 eV / c 2 , for example, one obtains c 5 / 3 = 7.42 × 10 − 32 J ⋅ m − 2 For γ = 1 and 2 one finds c 1 ≈ m ν c 2 = 16 × 10 − 20 J and c 2 ≈ ћ 3 / ( m ν 2 c ) = 1.22 × 10 − 39 J ⋅ m 3 , respectively.

This initial general development also considers independent variations of pressure with density and temperature, ( ∂ P / ∂ ρ ) T and ( ∂ P / ∂ T ) ρ Moving on to the second term of Equation (3), the gravitational potential energy term is of the usual form, involving the enclosed mass computed starting from the center of mass of the overall mass distribution. It also has units of energy density. One may use the kinetic and potential energy density expressions of Equation (3) to find the density of DM or OM versus radius in a spherically-symmetric geometry for a single species of matter. This is done by taking variations of the Lagrangian density, K.E.-P.E. with respect to ρ and T. The result for DM or OM is (with ordinary matter or dark matter held fixed, respectively)

[ ( ∂ P / ∂ ρ ) T + m v M e n c ( r ) G / r ] δ ρ ( r ) + 4 π ρ m v 2 [ ∫ 0 r r ′ 2 d r ′ δ ρ ( r ′ ) ] G / r + ( ∂ P / ∂ T ) ρ δ T ( r ) = 0. (6)

This equation should hold for arbitrary small variations δ ρ ( r ) and δ T ( r ) at each point r in mechanical and thermodynamic equilibrium at temperature T(r) and density ρ ( r ) . The two last terms look like they could pose a problem. To address the last term, assume Equation (2) applies, so that T is a function of ρ only, and δ T ( r ) = ( d T / d ρ ) δ ρ ( r ) . Then

[ ( ∂ P / ∂ ρ ) T + ( ∂ P / ∂ T ) ρ d T / d ρ + m v M e n c ( r ) G / r ] δ ρ ( r ) + 4 π ρ m v 2 [ ∫ 0 r r ′ 2 d r ′ δ ρ ( r ′ ) ] G / r . (7)

Note that the first two terms are d P / d ρ . To address the last term, differentiate Equation (7) with respect to r to obtain

[ d / d r ( ∂ P / ∂ ρ ) + 4 π r 2 m v ( ρ m v + ρ m , r a d ) G / r − m v M e n c ( r ) G / r 2 + 4 π ρ m v 2 r 2 G / r ] δ ρ ( r ) + [ d P / d ρ + m v M e n c ( r ) G / r ] d δ ρ ( r ) / d r + [ d / d r ( ρ / r ) ] m v 2 G [ 4 π ∫ 0 r r ′ 2 d r ′ δ ρ ( r ′ ) ] . (8)

The integral in the last term in Equation (8) can then be eliminated using Equation (7) to give

[ d / d r ( ∂ P / ∂ ρ ) + 4 π r 2 m v ( 2 ρ m v + ρ m , r a d ) G / r − m v M e n c ( r ) G / r 2 ] δ ρ ( r ) + 4 π ρ m v 2 r 2 G / r ] δ ρ ( r ) + [ d P / d ρ + m v M e n c ( r ) G / r ] d / d r ( δ ρ ( r ) ) − ( ρ / r ) − 1 [ d / d r ( ρ / r ) ] [ d P / d ρ + m v M e n c ( r ) G / r ] δ ρ ( r ) = 0. (9)

The second term in Equation (9) is the change in the potential energy of the system at radius r due to mass at radius r, and can be dropped if the enclosed mass is understood to be the enclosed mass strictly less than r.

To address the term involving d / d r [ δ ρ ( r ) ] the second line in Equation (9) can be integrated by parts, accounting for the volume integral of the Lagrangian density, as well as the r^{2} factor that is associated with it. After considerable algebra, the result is

{ [ d P / d ρ + m v M e n c ( r ) G / r ] ( 1 / ρ ) d ρ / d r + [ m v M e n c ( r ) G / r + d P / d ρ ] / r } δ ρ ( r ) = 0. (10)

This equation now has the required form so that the expression in curly brackets is zero:

( 1 / ρ ) d P / d r + ( 1 / ρ ) ( d ρ / d r ) [ m v M e n c ( r ) G / r ] = − m v M e n c ( r ) G / r 2 − ( d P / d ρ ) / r . (11)

Equation (11) is the result for strict equilibrium for a medium with a pressure that depends on density only in a spherically-symmetric geometry, and with a single species of matter. The limiting case of the standard hydrostatic equation is obtained by neglect of the second terms on both sides of Equation (11).

Based on the usual derivation of Equation (1), the standard hydrostatic equation should be valid when mechanical equilibrium applies but the stricter action equilibrium of Equation (11) does not apply. The term involving ( 1 / ρ ) ( d ρ / d r ) in Equation (11) arises because of the energy required to maintain a density gradient in a gravity potential, and the term involving − ( d P / d ρ ) / r arises because of the form of the divergence of the radial part of the bulk modulus in a spherical geometry.

It is also worth noting that Equation (11) can be written as

d ln ( ρ ) / d ln ( r ) = − 1 . (12)

This simplification indicates a solution for density of the form

ρ ~ ρ 0 ( r 0 / r ) . (13)

This result is both remarkable and perhaps expected. It arises because of the spherical geometry and the r^{−1} potential as discussed above. It matches the Navarro-Frenk-White (NFW) result [

{ [ 1 + r λ M / ( ρ m v 2 G ) ] − 1 [ ( 1 / ρ ) d ρ / d r ] [ d P / d ρ + m v M e n c ( r ) G / r + λ ρ ] } δ ρ ( r ) + ( { 2 − [ 1 + r λ M / ρ m v 2 G ] − 1 } [ d P / d ρ + m v M e n c ( r ) G / r + λ ρ ] / r ) δ ρ ( r ) + [ λ M ( 4 π m v 2 G ) d / d r ( ρ / r ) ] ∫ r ∞ d r ′ r ′ 2 δ ρ ( r ′ ) / [ λ M + ρ m v 2 G / r ] = 0. (14)

The last term arises from the variation of the total-mass constraint with respect to density. For r small (close to the origin) or large (far from the origin), one might expect this integral to be approximately zero, since the total mass constraint requires that the radial integral from zero to infinity of the density variations be identically zero. This version of the equation will be explored in subsequent sections of this paper.

As a basic check of the results of Section 2, one may consider the case of earth’s atmosphere. In the earth’s atmosphere, it is well-known that (a) the temperature is not constant with altitude in the troposphere [

First the version of the generalized equation for hydrostatic equilibrium given by Equation (11) is considered. For the earth’s atmosphere at distance z above sea level, one has

( 1 / ρ ) d P / d r + ( 1 / ρ ) ( d ρ / d z ) m a i r g z = − m a i r g − d P / d ρ / ( r e + z ) . (15)

here m a i r is the atomic weight of air, approximately (29.0) (1.67 × 10^{−27}) kg, and g is the gravitational constant at or near earth’s surface, 9.8 m·sec^{−2}. The radius r e is measured from earth center. Next consider an equation of state

P ( z ) = k ρ ( z ) T ( z ) , (16)

where T(z) is the temperature, which is treated as an input from [^{−28} J·m^{−1} and that m a i r g = 4.7 × 10 − 25 J ⋅ m − 1 . Hence the last term of Equation (15) can be neglected relative to m a i r g . Using expression (16) and Equation (15) one obtains

( k / ρ ) [ d ρ / d z T ( z ) + ρ ( z ) d T / d z ] + ( 1 / ρ ) ( d ρ / d z ) m a i r g z = − m a i r g . (17)

This equation simplifies to (neglecting terms of order z/r_{e}),

( 1 / ρ ) ( d ρ / d z ) = [ − m a i r g − k d T / d z ] / [ k T ( z ) + m a i r g z ] . (18)

This can be solved analytically:

ρ ( z ) = ρ 0 k T ( z = 0 ) / [ k T ( z ) + m a i r g z ] . (19)

here ρ 0 is the density at sea level, which is an input. The solution of Equation (19) gives poor agreement with measured data, with significantly larger density than measured at higher altitudes. This will be seen to arise from the underlying assumption of strict action (mechanical and thermal) equilibrium. Equation (19) may be compared to the isothermal approximation, which amounts to retention of only the first terms on both sides of Equation (17) and assuming constant temperature. One obtains the usual exponential atmosphere in this case,

ρ ( z ) = ρ 0 exp [ − m a i r g z / k T 0 ] , (20)

where T_{0} is the temperature of the atmosphere at sea level.

There is also the adiabatic equation for pressure and density. The above Equations (17)-(19) assumes strict action equilibrium. The adiabatic equation assumes that hot parcels of air near the surface of the earth rise without significant heat exchange with the surrounding air. Hence it corresponds to a process in quasi-static equilibrium that is not in strict mechanical and thermal equilibrium at all layers. With this insight, the adiabatic result can be derived from Equation (17) by neglect of the last terms on both sides of the equation, and assuming the gas satisfies a polytropic relation (2) with a polytropic exponent γ ~ 1.4 . With these assumptions, one obtains the result from the standard hydrostatic equation, which is the well-known formula for density [

ρ ( z ) = ρ 0 [ 1 − m a i r g z / ( c p T 0 ) ] 1 / ( γ − 1 ) . (21)

here c_{p} is the specific heat at constant pressure of the atmosphere at sea level. Note that c p T 0 = [ ( γ − 1 ) / γ ] P 0 / ρ 0 , where P_{0} is the pressure at sea level. For _{p} used is 1004.7 J·kg^{−1}·K^{−1}, the value of T_{0} used is 290.1 K. The latter is consistent with the ideal gas law and the sea-level values of pressure and density given in [

It should be noted that Equation (21) gives nonsensical results for m a i r g z / ( c p T 0 ) > 1 , which corresponds to an altitude z of about 30 km. However, the result for density is quite accurate for the lowest 15 km of the atmosphere. To remedy the errors above 15 km, the additional terms in Equation (15) can be restored gradually above 15 km, representing the return to true equilibrium at higher altitudes. The altitude dependent factor used to restore the full Equation (15) is given by Equation (22):

f ( z ) = H ( z − 15 ) { 1 − exp [ − ( z − 15 ) 2 / 20 2 ] } , (22)

where z is in km in this expression and H(z) is the Heaviside function. This factor gives a gaussian onset to equilibrium and multiplies the terms ( 1 / ρ ) ( d ρ / d z ) m a i r g z and d P / d ρ / ( r e + z ) in Equation (15). A numerical integration of (15) with this factor is shown in

Comparing the cyan and blue curves for both pressure and density in

As discussed at the end of Section 1, the unconstrained equilibrium solution gives infinite density at the origin. This defies the Fermi-Dirac equation, which implies a maximum density due to fermion degeneracy pressure. It also seems counter-intuitive, at least in the absence of a gravitational singularity. To remedy this, a density constraint (as well as a total particle constraint) is added to the Lagrangian density, with a result given by Equation (14). It is not difficult to see that Equation (14) can be written

{ d ln ( ρ ) / d ln ( r ) + 1 + 2 [ r λ M / ( ρ m v 2 G ) ] } δ ρ ( r ) + [ ( 4 π m v 2 G ) r d / d r ( ρ / r ) ] ∫ r ∞ d r ′ r ′ 2 δ ρ ( r ′ ) / [ 1 + ρ m v 2 G / ( r λ M ) ] × [ 1 + r λ M / ( ρ m v 2 G ) ] / [ d P / d ρ + m v M e n c ( r ) G / r + λ ρ ] = 0 (23)

The last term of Equation (23), comprising the last two lines, can be estimated in several ways. The simplest approach is to assume that the density variations in the integral average to zero, yielding zero for this term. However, a more careful approach is appropriate, and this can be done assuming the zeroth-order solution, Equation (13). It is further assumed r λ M / ( ρ m v 2 G ) is much less than 1, to investigate behavior near the origin. Then the various factors in the last term of Equation (23) are

r d / d r ( ρ / r ) = − 2 ρ 0 r 0 / r 2 , (24a)

1 + r λ M / ( ρ m v 2 G ) = 1 + O ( r λ M / ( ρ m v 2 G ) ) , (24b)

[ 1 + ρ m v 2 G / ( r λ M ) ] − 1 = r 2 λ M / ( ρ 0 r 0 m v 2 G ) + o ( r λ M / ( ρ m v 2 G ) ) , (24c)

And m v M e n c ( r ) G / r = 2 π m v 2 G ρ 0 r 0 r ≡ C 0 r . (24d)

The “big-oh” and “little-oh” notations are used in Equations (24b) and (24c). With Equations (24), the last term of Equation (23) can then be estimated as

− 8 π λ M ∫ r ∞ d r ′ r ′ 2 δ ρ ( r ′ ) / ( d P / d ρ + C 0 r + λ ρ ) . (25)

Note also that with constrained total particle number, the variations should satisfy

4 π ∫ 0 ∞ d r ′ r ′ 2 δ ρ ( r ′ ) = 0 , (26)

which allows Equation (25) to be re-written as

+ 8 π λ M ∫ 0 r d r ′ r ′ 2 δ ρ ( r ′ ) / ( d P / d ρ + C 0 r + λ ρ ) . (27)

With bounded variations δ ρ ( r ′ ) in 0 to r, this expression tends to zero as r tends to zero, when the denominator is not zero (as should be the case). Hence, however the last term is estimated, the summary result for Equation (23) near the origin is

{ d ln ( ρ ) / d ln ( r ) + 1 + 2 [ r λ M / ρ m v 2 G ] } δ ρ ( r ) = 0 . (28)

Note that this has two possible solutions for any given r. The first solution is obtained by the usual approach, setting the quantity inside curly brackets to zero. This gives the r^{−1} solution presented in Equation (13) near the origin, when r λ M / ( ρ m v 2 G ) is much less than 1. A second solution to Equation (28) is obtained by setting the variation δ ρ ( r ) to zero, i.e.,

ρ = constant in a spherical region. (29)

This latter solution is consistent with an incompressible medium, as occurs when the density is limited by Fermi-Dirac statistics. It is also consistent with the observed features of most celestial bodies, in which the density is approximately constant at their core, and also often in spherical shells. It is also approximately consistent with observations for dark matter in central regions that are less than a kpc in diameter [

Next consider the standard hydrostatic equation. One can assume a solution of the form ρ 0 r 0 r − 1 for the density and substitute into Equation (1) with Equation (2) and check for consistency. One finds that

( 1 / ρ ) ( d P / d r ) = − m ν M e n c G / r 2 → c γ γ ρ γ − 2 d ρ / d r = − 2 π m ν 2 G ρ 0 r 0 . (30)

One notes that the left-hand side varies as r − γ whereas the right-hand side is a constant under the assumption of ρ ~ r − 1 . Hence for any γ except possibly zero the standard hydrostatic equation is not consistent with an r^{−1} solution.

One can also analyze the standard equation for hydrostatic equilibrium in a gravitational field assuming constant number density ρ 0 in the vicinity of the origin of a spherical mass distribution. A constant density near the origin is self-consistent in Equation (1) as can be seen as follows. Assuming a polytropic relation P = c γ ρ γ in Equation (1) one obtains

c γ γ ρ γ − 2 d ρ / d r = − ( 4 π / 3 ) m v 2 G ρ 0 r . (31)

here ρ 0 is the sum of OM and DM mass density at the origin, normalized to the mass of the DM particles, with the assumption that OM dominates. This expression shows that the derivative of the dark or ordinary matter is zero at the origin, and gradually becomes more negative, consistent with constant density. A solution of Equation (31) for dark matter gives the following result for γ > 1 :

ρ ( r ) = [ ρ 01 γ − 1 − ρ 01 γ − 1 C 1 r 2 ] 1 / ( γ − 1 ) = ρ 01 [ 1 − C 1 r 2 ] 1 / ( γ − 1 ) = ρ 01 [ 1 − ( r / r a ) 2 ] 1 / ( γ − 1 ) , (32)

where ρ 01 is the (equivalent) number density at the origin, and C 1 = ( 2 π / 3 ) m v 2 G ρ 0 ( γ − 1 ) / ( c γ γ ρ 01 γ − 1 ) . Note that C 1 has units of ·m^{−2}, and so one may set C 1 ≡ r a − 2 , and this is used in the last equality of Equation (32). The solution is clearly approximate because for r greater than r a the density can be become complex and also because the density variation of DM in Equation (32) is not included in the enclosed mass in Equation (31). However, it does show that a constant density to order of ( r / r a ) 2 is consistent with the standard hydrostatic equation.

To summarize Section 4, it is found that setting the variations of the Lagrangian density to zero yields an alternative, constant-density solution. This solution is consistent with the observed density profile of many celestial bodies near their origin, including observationally-inferred dark matter profiles. The standard hydrostatic equation also shows solutions consistent with approximately constant density at the origin but is seen to be inconsistent with r^{−1} solutions. With these considerations in mind, the constant-density solution is used in a region of radius r_{0} near the origin in the following sections.

The mass-constrained Equation (14) is most relevant here. The last term in the equation is assumed negligible in this case, since the integral of any physical density variations from r to infinity should tend to zero for large r in the case of constrained total particle number. The result is

[ 1 + r λ M / ( ρ m v 2 G ) ] − 1 [ ( 1 / ρ ) d ρ / d r ] [ d P / d ρ + m v M e n c ( r ) G / r + λ ρ ] + { 2 − [ 1 + r λ M / ρ m v 2 G ] − 1 } [ d P / d ρ + m v M e n c ( r ) G / r + λ ρ ] / r = 0. (33)

Inspection shows that this equation can be simplified considerably, giving, as in Equation (28),

d ln ( ρ ) / d ln ( r ) = − 1 − 2 [ r λ M / ρ m 2 G ] . (34)

This result shows that the solution is independent of the density constraint λ ρ with the neglect of the last term of Equation (14). This equation may further be put into dimensionless form by defining λ M / ( m v 2 G ) ≡ ρ c / r c , where is ρ c the cutoff density and r c is the cutoff radius where the density drops due to the constraint on the total particle count. The resulting equation is

d ln ( ρ ) / d ln ( r ) = − 1 − 2 [ ( r / r c ) ( ρ c / ρ ) ] . (35)

This equation may be solved numerically or perturbatively. For a perturbative treatment, assume ( r / r c ) ( ρ c / ρ ) is much less than 1, corresponding to radii near the origin but not at the origin. Then the zeroth order solution is given by Equation (12),

ρ ( r ) = ρ 0 ( r 0 / r ) + O [ ( r / r c ) ( ρ c / ρ ) ] . (36)

where the “big-oh” notation is used to indicate the order of the approximation. Substituting this solution back into Equation (35) gives

d ln ( ρ ) / d ln ( r ) = − 1 − 2 [ ( r 2 / r c r 0 ) ( ρ c / ρ 0 ) ] . (37)

This equation is readily solved to give

ρ ( r ) = ρ 0 ( r 0 / r ) exp { − ( r 2 − r 0 2 ) [ ρ c / ( ρ 0 r 0 r c ) ] } + O { [ ( r / r c ) ( ρ c / ρ ) ] 2 } = ρ 0 ( r 0 / r ) { 1 − ( r 2 − r 0 2 ) [ ρ c / ( ρ 0 r 0 r c ) ] } + O { [ ( r / r c ) ( ρ c / ρ ) ] 2 } = ρ 0 r 0 r [ 1 + ( r 2 r 0 2 − 1 ) ( ρ c r 0 ρ 0 r c ) ] + O { [ ( r / r c ) ( ρ c / ρ ) ] 2 } = ρ 0 r 01 r ( 1 + r 2 / r 02 2 ) + O { [ ( r / r c ) ( ρ c / ρ ) ] 2 } (38)

In the last expression, r 01 = r 0 / ( 1 − μ ) , r 02 2 = r 0 2 ( 1 − μ ) / μ , and μ = ρ c r 0 / ( ρ 0 r c ) . With μ less than 1 and [ ( r / r c ) ( ρ c / ρ ) ] 2 < 1 , this expression is similar to the NFW result with a 1/r dependence near the origin and a 1/r^{3} dependence far from the origin. It is also similar to the double-power-law models and in particular the ( α , β , γ ) models [

To check the above results, Equation (35) may be integrated numerically from the origin for some relevant input values. In accord with Section 4, the density is assumed constant for r < r 0 . The resulting density is shown in plot (a) of _{c} and for Equation (38) for r less than about 0.7r_{c}. The plot on the right shows the normalized enclosed-mass curve for the numerical result in plot (a). The normalization of enclosed mass uses the results of Section 4, accounting for the region of radius r_{0} of constant density ρ 0 near the origin and a 1/r dependence outside this region to a radius of 0.7r_{c} (at which the 1/r dependence is seen to be no longer valid). The result is

M e n c , t o t ≡ 2 π m v ρ 0 r 0 [ ( 0.7 r c ) 2 − r 0 2 / 3 ] . (39)

It should be noted that the results of

In both this section and Section 4, the last term of Equation (14) is neglected with some justification. However, there is no particular justification for its neglect in regions neither near nor far from the origin. The retention of this last term could introduce perturbations into the nominal equation which are arbitrary, but which might average in some sense to zero. It should also be observed that in cases of interest, the parameter λ M is extremely small compared to the other energy densities involved, so in such cases the neglect of the last term of Equation (14) is justified. For example, with λ M = m v 2 G ρ c / r c , with r_{c} = 92 kpc for a halo like that of the Milky Way [^{2} and 5 keV/c^{2}, one obtains λ M = 2.13 × 10 − 90 and 1.89 × 10^{−77} J·m^{−3}, respectively, using Equation (4). These values of λ M should be compared with the pressure at the origin, which from Equations (4) and (5) is about 3.22 × 10^{−10} and 1.14 × 10^{−4} J·m^{−3} for the same two masses, respectively. Even with densities 200 times lower as occur might occur at the edge of the density profile, the value of λ M is negligible in comparison. These considerations provide some justification for neglect of the last term of Equation (14).

Taking all the above sections into consideration, one finds that r λ M / ( ρ m v 2 G ) = ( r / r c ) ( ρ c / ρ ) is always less than 1 where there is appreciable density, so the resulting solution is

ρ ( r ) = { ρ 0 inasphericalregionabouttheoriginofradius r 0 , or ρ ( r ) satisfies d ln ρ d ln r = − 1 − 2 ( r r c ) ( ρ c ρ ) (40)

This solution has a central core region that is qualitatively consistent with the size of inner profiles of low-spectral brightness galaxies [

For the region just outside the core, the logarithmic slope given by Equation (40) is −1 or less, which is consistent with most past results. A comparison of the solution to the ( α , β , γ ) = ( 2 , 3 , 1 ) model is given in the previous section. Equation (40) can also be compared to the well-known de-projected Sersic and Einasto models [

ρ r a d = ρ r a d 0 ( r / R e ) − p n exp [ − b n ( r / R e ) 1 / n ] . (41)

The parameter ρ r a d 0 is obtained by setting the volume integral of Equation (41) equal to the measured or inferred ordinary mass of the galaxy. The variable R e denotes the radius which encloses 1/2 the total light of the galaxy. The dpS profile is also used with some success for characterization of DM density profiles versus radius. The other two parameters in Equation (41) are given conveniently and approximately from Equations (19) and (27) of [

p n = 1.0 − 0.6097 / n + 0.05463 / n 2 , and (42)

b n = 2 n − 1 / 3 + 0.009876 / n . (43)

The parameter n ranges from roughly 2 to 4 for DM haloes of galaxies and clusters, from the same reference. With these values of n, the leading exponent p n in Equation (42) is −0.7 to −0.85, which is not too different from −1, and [

d ln ( ρ d p S ) / d ln ( r ) = − p n − ( b n / n ) ( r / R e ) 1 / n . (44)

Another density distribution used for characterization of DM is the Einasto distribution, which for the purposes of this paper is given by

ρ E i n = ρ 0 exp [ − d n ( r / r e ) 1 / n ] , (45)

where r e is the exp ( − d n ) radius of the DM galactic halo, n is typically 4 to 7 in this case, and d n is given approximately by Equation (24) of [

d n ≈ 3 n − 1 / 3 + 0.0079 / n , for n > 0.5 . (46)

For the Einasto distribution, the logarithmic derivative is

d ln ( ρ E i n ) / d ln ( r ) = − ( d n / n ) ( r / r e ) 1 / n . (47)

Inspection of Equation (44) finds a total logarithmic slope of about − 3 + 1 / ( 3 n ) + 0.6097 / n at r = R e for the dpS model. The leading term matches the logarithmic slope of Equation (40) at ( r / r c ) ( ρ c / ρ ) = 1 . Note that the onset of steeper logarithmic slopes occurs for smaller r with the dpS model than for Equation (40). The primary potential disadvantage of the dpS model relative to Equation (40) is the infinite density at the origin (but if there is a gravitational singularity this may be an advantage).

The Einasto model does not have infinite density at the origin. The logarithmic slope of the Einasto model is also zero at the origin, matching the core solution of Equation (40). These may be viewed as positive features if the true density distribution is not infinite at the origin. However, the derivative d ρ E i n / d r is infinite at the origin for n > 1. The logarithmic slope of the Einasto model also gives a value of about −3 + 1/(3n) at r = r e using Equation (46). This is an approximate match of the logarithmic slope of Equation (40) at ( r / r c ) ( ρ c / ρ ) = 1 . However, as with the dpS model, the onset of steeper logarithmic slopes, slopes less than −1 but greater than −3, occurs for smaller r than for Equation (40).

This different behavior of the model solutions can be addressed with the possible solutions or extensions of Equation (14). One approach is to add angular momentum. The net effect of angular momentum would be to limit the extent of the galactic halo and would result in a more rapid decline in density than given in Equation (40). However, faster declines are found in simulations without appeal to angular momentum. As a second approach, one may consider the results of Section 3, which indicate that if the medium is not in a strict equilibrium, then solutions with faster declines occur. The outer perimeters of a large halo are not expected to be in strict equilibrium because of the process of aggregation of subhaloes [

( 1 / ρ ) d P / d r = − m v M e n c ( r ) G / r 2 [ 1 + 2 ( r / r c ) ( ρ c / ρ ) ] . (48)

Consider the region r which is significantly greater than r 0 but also significantly less than r c ( ρ / ρ c ) , to address the issue identified in the previous paragraph. Using the techniques of Equations (30) to (32), one may assume an initial ρ 0 ( r 0 / r ) solution to evaluate the enclosed mass, but this time without the expectation for consistency of the resulting density profile. One obtains an equation similar to Equation (21) in this range of r:

ρ ( r ) = ρ 0 [ 1 − r / r γ ] 1 / ( γ − 1 ) = ρ 0 [ 1 + r / r γ ] − 1 / γ − 1 + o ( r / r γ ) + o [ ( r / r c ) ( ρ c / ρ ) ] (49)

where

r γ = c γ γ / [ 2 π ( γ − 1 ) G m υ 2 r 0 ρ 0 2 − γ ] . (50)

For related states of matter, such as neutron stars, a polytropic exponent γ ranging from 3/2 to 2 is often used [

One may also compare the logarithmic derivative of Equation (48) to the dpS and Einasto forms. From Equation (48), the logarithmic derivative is given by

d ln ( ρ ) / d ln ( r ) = − m v M e n c ( r ) G / r [ 1 + 2 ( r / r c ) ( ρ c / ρ ) ] / c γ γ ρ γ − 1 . (51)

Substituting ρ ~ 1 / r s , s > 0, gives a right-hand side proportional to r 2 + s ( γ − 2 ) for ( r / r c ) ( ρ c / ρ ) < 0.5 . This is always negative (with the overall minus sign on the right of Equation (51)). When the right-hand side is less than −1, the density will decline faster than the r^{−1} solution. One can match the r 1 / n dependence of the dpS and Einasto logarithmic derivatives of Equations (44) and (47) with exponent n of 2 to 7, by setting 2 + s ( γ − 2 ) equal to n^{−1}. The requisite values of s are in the range of ( 13 / 7 ) / ( 2 − γ ) to ( 3 / 2 ) / ( 2 − γ ) , for γ < 2 and in the stated range of values of r. Hence, the overall radial dependence of the Einasto and de-projected Sersic models can be obtained with the assumption of deviation from equilibrium, using insights from Section 3. One also observes that including the mass constraint term will invariably result in convergent mass for decreasing density with increasing radius.

It is also instructive to consider the case in which constituent particles are non-relativistic and use the standard hydrostatic relation Equation (1) and the polytropic relation Equation (2). This leads to a familiar equation for the density and provides an analytic estimate for the scale size of cosmic structures. When Equation (2) for non-relativistic fermions is inserted in the standard equation for hydrostatic equilibrium, Equation (1), one obtains

( 5 / 3 ) c 5 / 3 ρ − 1 / 3 d ρ / d r = − m ν M e n c ( r ) G / r 2 . (52)

Multiplying Equation (52) by r^{2}, and then including the ρ − 1 / 3 factor in the differential, one obtains

( 5 / 2 ) r 2 c 5 / 3 d ρ 2 / 3 / d r = − m ν M e n c ( r ) G . (53)

Dividing by m υ 2 G and applying the operator r − 2 d / d r to both sides of Equation (53) gives

[ 5 c 5 / 3 / ( 2 m υ 2 G ) ] r − 2 d / d r ( r 2 d ρ 5 / 3 / d r ) = C 5 / 3 ∇ 2 ρ 2 / 3 = − r − 2 d / d r M e n c ( r ) / m υ (54)

where C 5 / 3 = 5 c 5 / 3 / ( 2 m υ 2 G ) . Note that the first equality in Equation (54) assumes spherical symmetry. Account for ordinary radiant matter as part of the enclosed mass. Then

M e n c ( r ) = 4 π ∫ 0 r d r ′ r ′ 2 [ m ν ρ ( r ′ ) + ρ m , r a d ( r ′ ) ] , (55)

where ρ m , r a d ( r ) is the density of ordinary matter, as defined in Section 2. Also normalize the number density by its value ρ 0 at the origin. With the notatio ρ 1 = ρ / ρ 0 the resulting dimensionless equation is

[ 5 c 5 / 3 / ( 8 π m υ 2 G ρ 0 1 / 3 ) ] ∇ 2 ρ 1 2 / 3 = − ρ 1 − ρ m , r a d / ( m ν ρ 0 ) . (56)

This equation is a form of the well-known Lane-Emden equation for polytropic spheres [

L 5 / 3 = [ 5 c 5 / 3 / ( 8 π m υ 2 G ρ 0 1 / 3 ) ] 1 / 2 (non-relativistic). (57)

has units of length and sets the scale of the decay of the density distribution. For relativistic fermions with polytropic exponent 4/3 the corresponding length scale is

L 4 / 3 = [ c 4 / 3 / ( π m υ 2 G ρ 0 2 / 3 ) ] 1 / 2 (relativistic). (58)

The parameter L 5 / 3 is shown in ^{3}^{/2} for non-relativistic fermions [^{2}, and lower temperature, 2 to 8 K, yield values of the scale parameter L 5 / 3 which correspond to the peak of the large-scale mass power spectrum of about 300 Mpc [^{2} and temperatures of 2 to 8 K correspond to scales sizes L 5 / 3 of superclusters of the order of 30 Mpc [^{2} or more for these temperatures give L 5 / 3 of 3 kpc or less, corresponding to the diameter of dwarf galaxies. One may note that the temperatures and masses shown correspond to non-relativistic conditions.

L_{5/3} | Mass (eV/c^{2}) | ||||||
---|---|---|---|---|---|---|---|

T (K) | 0.025 | 0.1 | 0.4 | 1.6 | 5.0 | 10 | 20 |

2 | 391 | 34.7 | 3.07 | 0.27 | 0.037 | 0.011 | 0.0033 |

8 | 275 | 24.5 | 2.17 | 0.192 | 0.026 | 0.0078 | 0.0023 |

100 | 133 | 12.7 | 1.14 | 0.102 | 0.0139 | 0.0041 | 0.0006 |

r_{5/3} | Mass (eV/c^{2}) | ||||||

T (K) | 0.025 | 0.1 | 0.4 | 1.6 | 5.0 | 10 | 20 |

2 | 3.1 × 10^{8} | 2.4 × 10^{6} | 1.9 × 10^{4} | 147 | 2.72 | 0.240 | 0.021 |

8 | 1.5 × 10^{8} | 1.2 × 10^{6} | 9390 | 73.4 | 1.36 | 0.120 | 0.011 |

100 | 3.6 × 10^{7} | 3.2 × 10^{5} | 2620 | 20.7 | 0.385 | 0.033 | 0.003 |

_{0} = 1 kpc. The table highlights the rapid variability of this scale parameter with mass. Masses of 1.2 to 1.6 eV/c^{2} give scale sizes corresponding to the peak of the large-scale mass power spectrum, at a nominal background temperature of 2 K, and masses of 10 to 50 eV/c^{2} yield scale sizes that correspond to radii of large to dwarf galaxies. Note that the masses shown in this table are not consistent with recent estimates of dark-matter particle mass [_{0}, the region in which the density is approximately constant and equal to ρ 0 . This region is a corollary of the assumption of fermionic matter in the absence of a gravitational singularity.

^{10} solar masses [^{8} m^{−3} out to a radius r_{c} of 300 Mpc. This choice of number density floor is based on the estimated neutrino density in standard cosmology [^{−1} solution. For this solution, the scale parameter r γ of Equation (50) and

Input Parameter | Values | Comment |
---|---|---|

Temperature at origin, T | 2, 25 K | Sets density at origin using Equation (4), μ υ = 0 . |

Masses of particles, m υ | 0.025 - 0.3 eV/c^{2} | |

Source mass at origin | 9 × 10^{10} solar masses | Similar to Milky Way OM mass |

Inner scale r_{0} | 77 kpc | Input to Equation (40) |

Temperature T_{out} at outer radius r_{c} | 0.1 K | Sets density ρ c at outer radius in Equation (40) using Equation (4). |

Polytropic exponent | 5/3 | Input to Lane-Emden equation |

In summary, this Section compares the results of the Lane-Emden equation with that of the generalized hydrostatic equation. It is seen that the former produces a different dependence of scale size on temperature and particle mass that does the generalized hydrostatic equation, and the shapes of the solutions are significantly different near the origin; the former is flat, the latter has a cusp at or near the origin. However, in both cases, a narrow range of masses and temperatures yield scale sizes that are consistent with observed galactic and large-scale structure as seen from

In summary, Section 2 shows a derivation of a generalized hydrostatic equation in spherically symmetric geometries with a single species of matter satisfying a polytropic relation. Section 3 provides a partial validation of the generalized hydrostatic Equation (10) using the earth’s atmosphere. Section 4 shows that there is a constant-density solution for the Lagrangian formulation in a spherical region about the origin (or in spherical regions about the origin), and that the r^{−1} solution near the origin is not self-consistent with the standard hydrostatic equation, as expected, when the added terms of the generalized equation are not included. Section 5 explains how Equation (14) gives a convergent enclosed mass far from the origin and shows an example solution of the differential equation. Section 6 demonstrates that the solutions can match the radial behavior of accepted models of OM and DM density versus radius away from the origin. Section 7 presents solutions of the generalized equation for the density and compares it to the standard hydrostatic equation with a polytropic exponent in a gravitational potential. The overall conclusion is that a generalized equation for hydrostatic equilibrium, Equation (40), gives a cuspy solution for DM that is more consistent with observations and simulations (e.g., [^{−1} density profile is largely independent of the properties of the constituent material, this profile should indeed be “ubiquitous,” as discussed by other authors, wherever there is strict equilibrium.

Portions of this work were presented in Paper APR19-000356 at the 2019 April Meeting of the American Physical Society.

The author declares no conflicts of interest regarding the publication of this paper.

Holmes, R.B. (2020) Density Profiles of Gases and Fluids in Gravitational Potentials from a Generalization of Hydrostatic Equilibrium. Journal of Modern Physics, 11, 648-667. https://doi.org/10.4236/jmp.2020.115042