To constrain the properties of dark matter, we study spiral galaxy rotation curves measured by the THINGS collaboration. A model that describes a mixture of two self-gravitating non-relativistic ideal gases, “baryons” and “dark matter”, reproduces the measured rotation curves within observational uncertainties. The model has four parameters that are obtained by minimizing a x 2 between the measured and calculated rotation curves. From these four parameters, we calculate derived galaxy parameters. We find that dark matter satisfies the Boltzmann distribution. The onset of Fermi-Dirac or Bose-Einstein degeneracy obtains disagreement with observations and we determine, with 99% confidence, that the mass of dark matter particles is m h> 16 eV if fermions, or m h > 45 eV if bosons. We measure the root-mean-square velocity of dark matter particles in the spiral galaxies. This observable is of cosmological origin and allows us to obtain the root-mean-square velocity of dark matter particles in the early universe when perturbations were still linear. Extrapolating to the past we obtain the expansion parameter at which dark matter particles become non-relativistic: ahNR=[4.17±0.34(STAT)±2.50(SYST)]×10−6. Knowing we then obtain the dark matter particle mass mh=69.0±4.2(stat)±31.0(syst)eV, and the ratio of dark matter-to-photon temperature Th/T=0.389±0.008(stat)±0.058(syst) after e+e− annihilation while dark matter remains ultra-relativistic. We repeat these measurements with ten galaxies with masses that span three orders of magnitude, and angular momenta that span five orders of magnitude, and obtain fairly consistent results. We conclude that dark matter was once in thermal equilibrium with the (pre?) Standard Model particles (hence the observed Boltzmann distribution) and then decoupled from the Standard Model and from self-annihilation at temperatures above m μ. These results disfavor models with freeze-out or freeze-in. We also measure the primordial amplitude of vector modes, and constrain the baryon-dark matter cross-section: . Finally, we consider sterile Majorana neutrinos as a dark matter candidate.
The dark matter density in the core of spiral galaxies can exceed 107 times the mean dark matter density of the Universe. To learn about the properties of dark matter, we study the rotation curves of spiral galaxies measured by the THINGS collaboration [
We use the standard notation for cosmology as defined in [
v obs 2 ≡ v 2 = v b 2 + v h 2 , v b 2 = v disk 2 + v bulge 2 + v gas 2 . (1)
These rotation velocities correspond, by definition, to test particles in circular orbits of radius r in the plane of the galaxy. There is good agreement between the THINGS analysis and previous measurements by K. Begeman [
From the measured rotation curves v ( r ) and v b ( r ) we can obtain directly g b = − v b 2 / r and g h = − ( v 2 − v b 2 ) / r . From g h = − M h ( r ) G / r 2 (assuming spherical symmetry) we can obtain
ρ h = 1 4 π G r 2 d ( r v h 2 ) d r , (2)
and similarly for ρ b . We use these equations as cross-checks. Note that ρ h (r)
and ρ b ( r ) are insensitive to whether or not the dark matter is rotating.
In the present analysis we would like to go further, i.e. constrain the equation of state of dark matter. To this end we compare the measured rotation curves with the following simplified model of a mixture of two self-gravitating non-relativistic ideal gases: “baryons” and “dark matter”. We assume that the interactions between baryons and dark matter can be neglected.
Galaxy | 〈 v r h 2 〉 ′ 1 / 2 [km/s] | 〈 v r b 2 〉 ′ 1 / 2 [km/s] | ρ h ( r min ) [ 10 − 2 M ⊙ ⋅ pc − 3 ] | ρ b ( r min ) [ 10 − 2 M ⊙ ⋅ pc − 3 ] | χ 2 / d .f . |
---|---|---|---|---|---|
NGC 2403 | 101 ± 3 | 63.4 ± 1.8 | 7.5 ± 1.4 | 30.4 ± 6.2 | 7.2/20 |
NGC 2841 | 220 ± 3 | 168 ± 3 | 9.3 ± 0.7 | 9.7 ± 1.2 | 22.1/22 |
NGC 2903 | 142 ± 3 | 117 ± 3 | 14.6 ± 2.1 | 25.2 ± 3.1 | 29.0/24 |
NGC 2976 | 129 ± 177 | 49 ± 6 | 4.0 ± 2.7 | 17.0 ± 3.7 | 9.8/16 |
NGC 3198 | 104 ± 3 | 90 ± 3 | 4.5 ± 0.8 | 4.9 ± 1.2 | 8.3/22 |
NGC 3521 | 153 ± 10 | 136 ± 5 | 22.9 ± 8.6 | 43.4 ± 11.8 | 6.8/20 |
NGC 3621 | 126 ± 5 | 78.8 ± 1.0 | 2.6 ± 0.5 | 26.6 ± 1.5 | 13.2/22 |
DDO 154 | 36.5 ± 3.7 | 22.3 ± 1.7 | 1.3 ± 0.3 | 1.9 ± 0.8 | 7.5/18 |
NGC 5055 | 144 ± 4 | 130 ± 3 | 28.2 ± 6.8 | 41.4 ± 10.1 | 19.2/24 |
NGC 7793 | 85.5 ± 5.0 | 52.2 ± 1.5 | 8.0 ± 1.6 | 29.7 ± 4.4 | 36.3/16 |
∇ ⋅ g b = − 4 π G ρ b , ∇ ⋅ g h = − 4 π G ρ h , (3)
g = g b + g h , g b ≡ − v b 2 r , g h ≡ − v h 2 r , v 2 ≡ v b 2 + v h 2 , (4)
∇ P b = ρ b ( g + κ b v 2 r e ^ r ) , ∇ P h = ρ h ( g + κ h v 2 r e ^ r ) , (5)
P b = 〈 v r b 2 〉 ρ b and P h = 〈 v r h 2 〉 ρ h ( 1 ± ρ h π 3 / 2 ℏ 3 2 〈 v r h 2 〉 3 / 2 N f , b m h 4 + ⋯ ) . (6)
〈 v r b 2 〉 and 〈 v r h 2 〉 are the mean-square of the velocity components in the radial direction r of baryons and dark matter respectively. Newton’s Equation (3) determine the gravitational field g ( r ) = g b ( r ) + g h ( r ) due to the densities ρ b ( r ) of baryons and ρ h ( r ) of dark matter. Equations (5) express momentum conservation and are valid even for collisionless gases. The centrifugal acceleration terms proportional to κ b and κ h are inserted to study the rotation of baryons and dark matter.
We consider dark matter to be a mixture of interacting or non-interacting particles of masses m h i . Let n h i ( r , v r ) d v r be the number density of dark matter particles of mass m h i with radial component of velocity between v r and v r + d v r . Then
ρ h = ∑ i m h i ∫ 0 ∞ n h i ( r , v r ) d v r , P h ≡ 〈 v r h 2 〉 ρ h , (7)
where
〈 v r h 2 〉 ≡ ∑ i m h i ∫ 0 ∞ v r 2 ⋅ n h i ( r , v r ) d v r ∑ i m h i ∫ 0 ∞ n h i ( r , v r ) d v r , (8)
and similarly for baryons. The pressure P h ( r ) is the momentum component in the radial direction in the galactic plane per unit time traversing unit area at r with v r > 0 .
Equations (6) are the equations of state of the gasses. For the dark matter halo we have included a term due to the onset of degeneracy of fermions (upper signs) or bosons (lower signs). N f ( N b ) is the number of fermion (boson) degrees of freedom. We take N f = 2 for one flavor of spin-up and spin-down sterile Majorana neutrinos. Results can be amended for other cases. For bosons we take N b = 1 . Fermi-Dirac and Bose-Einstein degeneracy will be considered in Section 7. For now we set m h to some large value, e.g. m h = 500 eV, so the last term in Equation (6) is negligible.
We integrate numerically Equations (3) to (6) from r min to r max along a radial direction in the galactic plane. Hence ∇ P h = e ^ r d P h / d r and ∇ ⋅ g h = ( 1 / r 2 ) d ( r 2 g h ) / d r , and similarly for baryons.
Variables κ b and 〈 v r b 2 〉 , and also κ h and 〈 v r h 2 〉 , occur in the combinations
〈 v r b 2 〉 ′ ≡ 〈 v r b 2 〉 1 − κ b , 〈 v r h 2 〉 ′ ≡ 〈 v r h 2 〉 1 − κ h , (9)
if the last term in (6) is negligible. To lift this degeneracy and obtain these variables separately, we need information in addition to the galaxy rotation curves.
We fit four parameters to minimize the χ 2 between the measured and calculated rotation velocities v ( r ) and v b ( r ) : the starting densities ρ h ( r min ) and ρ b ( r min ) , and the “reduced” root-mean-square radial velocities 〈 v r h 2 〉 ′ , and 〈 v r b 2 〉 ′ . These parameters are boundary conditions at r → 0 and r → ∞ . Assuming 〈 v r h 2 〉 ′ and 〈 v r b 2 〉 ′ are independent of r we obtain fits to the galaxy rotation curves that are in agreement with observations, within the experimental uncertainties, as shown in Figures 1-10, and in
We define M b ( r ) as the baryon mass contained within r, and similarly for M h ( r ) . Equations (3) to (6) allow the definition of several galaxy parameters: the “equal density radius” r ed with ρ b ( r ed ) = ρ h ( r ed ) , the “galaxy proper radius”
〈 v r h 2 〉 ′ 1 / 2 | 〈 v r b 2 〉 ′ 1 / 2 | ρ h ( r min ) | ρ b ( r min ) | |
---|---|---|---|---|
〈 v r h 2 〉 ′ 1 / 2 | 1.000 | −0.833 | −0.802 | 0.860 |
〈 v r b 2 〉 ′ 1 / 2 | −0.833 | 1.000 | 0.775 | −0.831 |
ρ h ( r min ) | −0.802 | 0.775 | 1.000 | −0.639 |
ρ b ( r min ) | 0.860 | −0.831 | −0.639 | 1.000 |
r g at which M b ( r g ) / M h ( r g ) = Ω b / Ω c , and the corresponding “primordial comoving radius” r c defined by M b ≡ 4 3 π r c 3 Ω b ρ crit . The “baryon mass” M b
converges, so we define M b ≡ M b ( r g ) . In the early Universe, when the perturbation that formed the galaxy was still linear, the matter of the galaxy was contained in the expanding sphere of proper radius r g , and comoving radius r c at
time t g . At time t g the expansion parameter is r g / r c . The results, obtained with the four fitted parameters listed in
As an example, let us consider the spiral galaxy NGC 2403. A distribution of velocities along a particular line of sight is presented in figure 3 of [
Consider an expanding sphere of proper radius r g and comoving radius r c at time t g containing the matter that will become spiral galaxy NGC 2403. We assume adiabatic primordial perturbations, so the velocity fields of baryons and
Galaxy | M b [ 10 10 M ⊙ ] | r ed [kpc] | r g [kpc] | r c [Mpc] | t g [Gyr] |
---|---|---|---|---|---|
NGC 2403 | 0.9 × e ± 0.22 | 2.2 × e ± 0.28 | 13.6 × e ± 0.24 | 0.7 × e ± 0.07 | 0.047 × e ± 0.25 |
NGC 2841 | 15 × e ± 0.13 | 4.1 × e ± 0.12 | 48.5 × e ± 0.18 | 1.8 × e ± 0.04 | 0.077 × e ± 0.20 |
NGC 2903 | 10 × e ± 0.46 | 3.5 × e ± 0.20 | 74 × e ± 0.47 | 1.6 × e ± 0.15 | 0.18 × e ± 0.48 |
NGC 2976 | 0.4 × e ± 0.58 | 1.8 × e ± 0.41 | 6.3 × e ± 0.47 | 0.5 × e ± 0.19 | 0.022 × e ± 0.45 |
NGC 3198 | 9 × e ± 1.0 | 1.9 × e ± 1.0 | 128 × e ± 1.0 | 1.5 × e ± 0.34 | 0.42 × e ± 1.0 |
NGC 3521 | 69 × e ± 1.8 | 5.0 × e ± 0.87 | 403 × e ± 2.0 | 3.0 × e ± 0.60 | 0.86 × e ± 1.5 |
NGC 3621 | 3.9 × e ± 0.18 | 5.8 × e ± 0.15 | 38 × e ± 0.20 | 1.1 × e ± 0.06 | 0.11 × e ± 0.22 |
DDO 154 | 0.07 × e ± 0.48 | 1.1 × e ± 0.65 | 7.0 × e ± 0.46 | 0.3 × e ± 0.16 | 0.062 × e ± 0.46 |
NGC 5055 | 62 × e ± 1.3 | 2.3 × e ± 0.77 | 416 × e ± 1.5 | 2.9 × e ± 0.42 | 0.95 × e ± 1.5 |
NGC 7793 | 0.5 × e ± 0.23 | 1.6 × e ± 0.19 | 10 × e ± 0.27 | 0.6 × e ± 0.08 | 0.04 × e ± 0.30 |
Galaxy | log 10 { L b / [ M ⊙ km / s ⋅ pc ] } | Fermion mh> | Boson mh> |
---|---|---|---|
NGC 2403 | 12.63 ± 0.17 | 18 eV | no limit |
NGC 2841 | 14.65 ± 0.15 | 12 eV | 35 eV |
NGC 2903 | 14.61 ± 0.44 | 18 eV | 40 eV |
NGC 2976 | 11.93 ± 0.43 | no limit | no limit |
NGC 3198 | 14.76 ± 0.94 | 16 eV | 45 eV |
NGC 3521 | 16.37 ± 1.58 | 10 eV | no limit |
NGC 3621 | 13.73 ± 0.15 | no limit | no limit |
DDO 154 | 10.86 ± 0.36 | 30 eV | no limit |
NGC 5055 | 16.32 ± 1.07 | no limit | no limit |
NGC 7793 | 12.13 ± 0.19 | 25 eV | no limit |
of dark matter are the same. Then the angular momenta of baryons and dark matter in the sphere of proper radius r g are in the ratio L b / L h = M b / M h ( r g ) = Ω b / Ω h . According to the model of hierarchical formation of galaxies [
to 1 / r 2 . If dark matter rotates with velocity κ h v the angular momentum of dark matter of the galaxy of radius r g is L h ≈ ( π / 8 ) M h κ h v r g . From the numerical integration for galaxy NGC 2403, L b = 0.15 M b κ b v r g with κ b = 0.98 . Solving for κ h we obtain κ h ≈ 0.143 .
For galaxy NGC 3198 we obtain L b = 0.16 M b κ b v r g with κ b = 0.98 , and κ h ≈ 0.163 .
Excellent fits to the rotation curves are obtained from Equations (3) to (6) with 〈 v r b 2 〉 ′ and 〈 v r h 2 〉 ′ independent of r, and non-degenerate dark matter. Since for the galaxies under consideration κ b ≈ 1 and κ h ≈ 0 , we will neglect the possible dependence of κ b and κ h on r. From (5) and (6) we obtain
ρ h ∝ exp [ − ϕ ( r ) 〈 v r h 2 〉 ′ ] , n h i ( r , v r ) ∝ exp [ − v r 2 / 2 − ϕ ( r ) 〈 v r h 2 〉 ′ ] , (10)
where ϕ ( r ) ≡ − ∫ g ( r ) d r is the gravitational potential. Note that v r 2 / 2 + ϕ ( r ) is conserved. So, within the observational uncertainties of the galaxy rotation curves, dark matter satisfies the Boltzmann distribution (10). This is a non-trivial
result. The derivation of the Boltzmann distribution in statistical mechanics assumes that the “system” under study is in “thermal equilibrium” with a “reservoir”,
and implies the equipartition theorem, e.g. 1 2 m i v i 2 = 1 2 m j v j 2 , which does not
hold between dark matter and baryons. There is no obvious “reservoir” interacting with the dark matter, so the well known derivation of the Boltzmann
distribution may not apply to dark matter, unless dark matter was once in thermal equilibrium with “something”. Does the observed Boltzmann distribution suggest that dark matter was once in statistical equilibrium with the primordial (pre?) Standard Model cosmological soup?
So far we have not considered the last term in Equation (6), i.e. we have set m h
to some high value. Note that we obtain excellent fits to the galaxy rotation curves assuming P h ( r ) = 〈 v r h 2 〉 ρ h ( r ) with 〈 v r h 2 〉 independent of r. By lowering the value of m h in the fits we can probe non-linearities between P h ( r ) and ρ h ( r ) . In particular, we study the onset of Fermi-Dirac or Bose-Einstein degeneracy of dark matter composed of particles of mass m h . The dark matter density and pressure are [
ρ h = 〈 v r h 2 〉 3 / 2 N f , b m h 4 ( 2 π ) 3 / 2 ℏ 3 [ e μ ′ 1 3 / 2 ∓ e 2 μ ′ 2 3 / 2 + e 3 μ ′ 3 3 / 2 ∓ e 4 μ ′ 4 3 / 2 + ⋯ ] , (11)
P h = 〈 v r h 2 〉 5 / 2 N f , b m h 4 ( 2 π ) 3 / 2 ℏ 3 [ e μ ′ 1 5 / 2 ∓ e 2 μ ′ 2 5 / 2 + e 3 μ ′ 3 5 / 2 ∓ e 4 μ ′ 4 5 / 2 + ⋯ ] , (12)
where μ ′ ≡ μ / ( k T h ) , and upper signs are for fermions, and lower signs are for bosons. μ is the chemical potential. In the course of integrating Equations (3) to (6) numerically we need to obtain ρ h ( r ) given P h ( r ) and 〈 v r h 2 〉 . We do this by solving (12) for e μ ′ and substituting the result in (11). We consider up to third order terms in the series of (11) and (12) until the onset of full degeneracy for fermions:
P h = ρ h 5 / 3 6 2 / 3 π 4 / 3 ℏ 2 5 N f 2 / 3 m h 8 / 3 , (13)
or the onset of Einstein condensation for bosons:
P h = 〈 v r h 2 〉 ρ h ⋅ 0.5136, (14)
where the numerical factor is the ratio of the two series with μ = 0 .
For bosons 0 < e μ ′ < 1 . Einstein condensation sets in when e μ ′ reaches 1. As m h is lowered we find that e μ ′ at r min may exceed 1 indicating the onset of Einstein condensation. However the χ 2 of the fit to the galaxy rotation curves already exclude this value of m h even when we increase the starting value of g ( r min ) = − v 2 / r min within experimental bounds. In conclusion, the galaxy rotation curves exclude Einstein condensation in the studied galaxies.
Suppose a galaxy (not in
Fits to rotation curves of galaxy NGC 3198 for several values of m h are summarized in
Fermions | |||||||||
---|---|---|---|---|---|---|---|---|---|
mh [eV]: | 500 | 100 | 30 | 20 | 17 | 16 | 15 | 14 | 13 |
χ2: | 8.2 | 8.1 | 10.1 | 9.9 | 12.2 | 19.3 | 30.5 | 44.9 | 61.4 |
Bosons | |||||||||
mh [eV]: | 500 | 100 | 60 | 50 | 45 | 40 | 37 | 35 | 33 |
χ2: | 8.2 | 8.3 | 10.0 | 13.6 | 17.7 | 24.9 | 31.4 | 36.3 | 41.2 |
m h > 16 eV for fermions, and m h > 45 eV for bosons, if they dominate the dark matter density. Einstein condensation sets in at m h = 35 eV.
As another example let us consider galaxy NGC 2403. For fermions, degeneracy sets in at m h ≈ 40 eV, full degeneracy begins at m h ≈ 30 eV, and the limit we can set based on the increase of the χ 2 of the fit by 9 units is m h = 18 eV. For bosons, degeneracy sets in at m h ≈ 50 eV, Einstein condensation begins at m h = 47 eV, and for NGC 2403 no limit on m h can be set based on the increase of χ 2 . Additional examples are presented in
In conclusion, the observed spiral galaxy rotation curves disfavor dark matter dominated by axion-like bosons with mass less than 45 eV.
To be specific let us consider galaxy NGC 2403. The observed rotation curves of this galaxy, presented in
A model of the hierarchical formation of galaxies is described in References [
v h rms ≈ 3 〈 v r h 2 〉 1 / 2 ( ρ ¯ h ρ h ( r min ) ) 1 / 3 ≈ 63 km / s . (15)
The absolute velocity v h rms has contributions from three spatial components, hence the factor 3 . The last factor is a correction corresponding to adiabatic compression. According to the model of the hierarchical formation of galaxies, Equation (15) should be valid even with merging of galaxies. A correction due to mixing, expected to be of order O(1), may be needed. Extrapolating to the past, we find that the expansion parameter at which dark matter particles become non-relativistic is
a h NR ≈ r g r c v h rms c ≈ 4.1 × 10 − 6 , (16)
and the corresponding photon temperature is T 0 / a h NR ≈ 57 eV. Note that v h rms ∝ 1 / a for non-relativistic particles, and the expansion parameter at time t g is a = r g / r c . Equation (16) assumes that dark matter particles decoupled while still ultra-relativistic. To be specific, we assume that particles of mass m h dominate dark matter, and that these particles have N b = 0 boson degrees of freedom, and N f = 2 fermion degrees of freedom. The number density of dark matter particles at expansion parameter a h NR can be calculated in two ways:
n h NR = Ω c ρ crit m h a h NR 3 = 1.20205 π 2 ( k T h NR ℏ c ) 3 { N b + 3 4 N f } , (17)
where the ultra-relativistic expression to the right assumes zero chemical potential. Substituting m h c 2 ≈ 3.15 k T h NR , we obtain m h ≈ 70 eV. The ratio of the temperatures of dark matter and photons is T h / T ≈ m h / ( 3.15 × 57 eV ) ≈ 0.39 after e + e − annihilation while dark matter remains ultra-relativistic. Repeating these calculations for other galaxies, and including statistical uncertainties derived from the uncertainties of the measured galaxy rotation curves, we obtain
Note that M b and r g drop out of the expression for a h NR :
a h NR = 3 〈 v r h 2 〉 1 / 2 c [ Ω c ρ crit ρ h ( r min ) ] 1 / 3 , (18)
with r min → 0 . This equation can be understood directly: it expresses 〈 v r h 2 〉 1 / 2 ∝ 1 / a , were the last factor in (18) is the “expansion parameter” in the core of the galaxy.
Galaxy | 10 6 × a h NR | mh [eV] | T h / T |
---|---|---|---|
NGC 2403 | 4.09 ± 0.38 | 70.0 ± 4.9 | 0.387 ± 0.009 |
NGC 2841* | 8.29 ± 0.33 | 41.2 ± 1.3 | 0.462 ± 0.006 |
NGC 2903 | 4.60 ± 0.32 | 64.1 ± 3.4 | 0.399 ± 0.008 |
NGC 2976 | 6.44 ± 2.33 | 49.8 ± 13.5 | 0.433 ± 0.039 |
NGC 3198 | 4.99 ± 0.44 | 60.3 ± 4.0 | 0.407 ± 0.010 |
NGC 3521 | 4.27 ± 0.82 | 67.8 ± 9.7 | 0.391 ± 0.019 |
NGC 3621 | 7.26 ± 0.76 | 45.5 ± 3.6 | 0.447 ± 0.012 |
DDO 154 | 2.65 ± 0.47 | 96.9 ± 13.0 | 0.347 ± 0.016 |
NGC 5055 | 3.75 ± 0.41 | 74.7 ± 6.1 | 0.379 ± 0.011 |
NGC 7793 | 3.39 ± 0.43 | 80.6 ± 7.6 | 0.369 ± 0.012 |
Average | 4.17 ± 0.34 | 69.0 ± 4.2 | 0.389 ± 0.008 |
Alternatively, (18) expresses adiabatic expansion of a non-relativistic ideal “noble” gas: T h V γ − 1 = constant with γ = 5 / 3 . By “noble” we mean that collisions (if any) between dark matter particles do not excite internal degrees of freedom (if any) of these particles.
The independent measurements of m h presented in
a h NR = [ 4.17 ± 0.34 ( stat ) ] × 10 − 6 , m h = 69.0 ± 4.2 ( stat ) eV , (19)
and the average T h / T , after e + e − annihilation while dark matter is still ultra-relativistic, is
T h / T = 0.389 ± 0.008 ( stat ) , (20)
where T is the photon temperature. Note that m h ∝ a h NR − 3 / 4 and T h / T ∝ a h NR 1 / 4 , so the correlated uncertainties of m h and T h / T can be derived directly from the uncertainty of a h NR .
The near coincidence of the two temperatures T h and T is evidence that dark matter was once in thermal equilibrium with the (pre?) Standard Model particles, which explains the observed Boltzmann distribution. From the preceding analysis we also conclude that dark matter decoupled while still ultra-relativistic, so scenarios with freeze-out, or freeze-in are disfavored.
Let us summarize. We assume that ultra-relativistic dark matter has zero chemical potential, and obtain, from the spiral galaxy rotation curves, a value of T h / T that is consistent with thermal equilibrium of dark matter with the Standard Model sector at some time in the early Universe. This is a highly significant result since T depends on T 0 while T h does not. If we would have obtained a different value of the “adiabatic invariant” 〈 v r h 2 〉 1 / 2 ρ h ( r min ) − 1 / 3 , we would have concluded that the chemical potential of dark matter is different from zero, and/or dark matter was never in thermal equilibrium with the Standard Model sector. So the observed adiabatic invariant is strong evidence that dark matter has zero chemical potential and was once in thermal equilibrium with the Standard Model particles.
The statistical uncertainties (at 68% confidence) presented so far are derived from the observational uncertainties of the galaxy rotation curves, and contributions from h and Ω c .
The main systematic uncertainty comes from Equation (15). This Equation is justified by the model of the hierarchical formation of galaxies [
Another correction of O(1) is due to the non-spherical symmetry of the galaxies. Yet another correction of O(1) may come from the details of the transition from ultra-relativistic to non-relativistic dark matter.
Pending detailed studies, our preliminary estimates of systematic uncertainties are:
a h NR = [ 4.17 ± 0.34 ( stat ) ± 2.50 ( syst ) ] × 10 − 6 , (21)
m h = 69.0 ± 4.2 ( stat ) ± 31.0 ( syst ) eV , (22)
T h / T = 0.389 ± 0.008 ( stat ) ± 0.058 ( syst ) . (23)
Note that (22) and (23) follow directly from (21): their uncertainties are correlated. In comparison, the current uncertainty of m h spans 70 orders of magnitude for fermions, and 90 orders of magnitude for bosons [
We consider dark matter that at some time in the history of the Universe was in thermal equilibrium with the (pre?) Standard Model particles, and decoupled from these particles and from self-annihilation while still ultra-relativistic. To be specific, we consider dark matter to be dominated by a single family of sterile Majorana neutrinos, i.e. N b = 0 boson degrees of freedom, and N f = 2 fermion degrees of freedom (for spin-up and spin-down). We also assume three families of active Majorana neutrinos with N f = 2 each.
As the universe expands and cools, Standard Model particles and antiparticles that become non-relativistic annihilate heating the Standard Model sector conserving entropy, without heating sterile and active neutrinos if they have already decoupled.
For example, assume the decoupling temperature is in the range m b to m W . Then T h / T is calculated as follows:
T h T = ( 4 × 43 11 × 345 ) 1 / 3 = 0.357. (24)
Decoupling temperature range | T h / T in range m h to m e | m h |
---|---|---|
m H to m t | 0.344 | 99.9 ± 3.1 eV |
m W to m H | 0.345 | 98.9 ± 3.1 eV |
m b to m W | 0.357 | 89.5 ± 2.8 eV |
m τ to m b | 0.372 | 78.6 ± 2.5 eV |
m c to m τ | 0.378 | 75.0 ± 2.4 eV |
m s to m c | 0.399 | 64.1 ± 2.0 eV |
m e to m μ | 0.714 | 11.2 ± 0.4 eV |
These numbers can be found in the
2 + 7 8 6 ( 4 11 ) 4 / 3 + 7 8 2 ( 4 × 43 11 × 345 ) 4 / 3 2 + 7 8 6 ( 4 11 ) 4 / 3 = 3.391 3.363 = 1.0084 , (25)
due to one family of sterile Majorana neutrinos. Equivalently, the effective number of light neutrinos increases from N ν = 3 to 3.062, well within BBN experimental bounds, i.e. N ν ≲ 4 [
An example of a sterile neutrino dark matter candidate is presented in
Consider a hydrogen atom, or a proton, of mass m p , in a circular orbit of velocity V b within the core of a galaxy. The velocity V b will decay with a relative rate τ b − 1 ≡ d V b / ( V b ⋅ d t ) ≈ ρ h σ h b V b / m p due to collisions with dark matter particles with cross-section σ h b . Taking τ b equal to the age of the universe, V b ≈ 150 km/s, and ρ h ≈ 0.2 M ⊙ / pc 3 (see
σ h b ≲ 2 × 10 − 26 cm 2 , (26)
independently of m h .
The mean time between collisions of a dark matter particle with baryons is τ h ≈ m p / ( ρ b σ h b 3 〈 v r h 2 〉 1 / 2 ) . Dark matter is not in thermal equilibrium with baryons so τ h is larger than the age of the universe. With data of galaxy NGC 2841 we obtain
σ h b ≲ 2 × 10 − 26 cm 2 . (27)
These limits are useful for dark matter masses m h below the reach of direct or indirect searches, e.g. m h < 50 keV. Planed direct low mass WIMP searches with the NEWS-SNO project will reach dark matter masses ≈ 0.1 GeV [
The mean baryon density of the universe is 6.2 × 10 − 9 M ⊙ / pc 3 , while the baryon density in the core of galaxy NGC 2403 is 0.3 M ⊙ / pc 3 , so matter has collapsed by a factor 5 × 10 7 in this case. What prevented ρ h ( r ) or ρ b ( r ) from collapsing all the way to ∞ as r → 0 ? In Figures 1-10 we see no increase of ρ h ( r ) or ρ b ( r ) as r → 0 down to the first measured point at r min . Why a core instead of a cusp? The four parameters ρ h ( r min ) , ρ b ( r min ) , 〈 v r h 2 〉 ′ , and 〈 v r b 2 〉 ′ can fit any core radius. So, what determines these four boundary conditions?
Possibilities to consider are Bose-Einstein of Fermi-Dirac degeneracy? (the observed galaxy rotation curves obtain m h at the edge of degeneracy as shown in Section 7), dark matter or baryon crossing time?, non-spherical initial density perturbations?, angular momentum of the primordial rotation field?, dominance of baryons in the core?, ongoing collapse of the core? (however baryons are already supported by the centrifugal force as shown in Section 4), etc.
In Figures 1-10 we observe that the baryon core radius, the radius r ed of equal densities, and the dark matter core radius are similar, and increase in this order. We also observe that the baryon density dominates at small radii. These observations suggest that baryons radiate energy conserving angular momentum, fall to the galactic plane, acquire increasingly circular orbits, until baryons are supported by the centrifugal force with κ b → 1 . Thereafter baryons can fall no further. For example, a proton or hydrogen atom falling into the galaxy collides with other baryons, falls to the disk and acquires circular velocity κ b v . Once the central gravitational potential is determined mostly by baryons, dark matter rests in this given potential.
Note that as r → ∞ , ρ h → v 2 / ( 4 π G r 2 ) with v 2 = 2 〈 v r h 2 〉 ′ . Note that (18) obtains ρ h ( 0 ) given 〈 v r h 2 〉 and a h NR , and therefore constrains the dark matter core radius.
The suggestion is that the primordial angular momentum in the sphere of proper radius r g and comoving radius r c during the expansion phase at time t g is the same as the angular momentum L = L b ( Ω c + Ω b ) / Ω b of the galaxy with L b listed in
Perturbations of density, gravitational potential, and velocity in the non-relativistic universe dominated by dark matter have four independent modes [
Most dark matter should have been non-relativistic at the time when there was a galactic mass inside the horizon [
At an expansion parameter a h NR ≈ 4.2 × 10 − 6 , when dark matter with mass m h ≈ 69 eV becomes non-relativistic, the universe is still dominated by radiation, see
Note that protons are non-relativistic when the smallest observed galaxies and globular clusters enter the horizon. Compared to the cold dark matter scenario with m h ≫ 5 keV, for m h = 69 eV the power spectrum of density fluctuations on scales smaller than M b = 20 × 10 10 M ⊙ becomes suppressed. If the dark matter inside the horizon becomes nearly homogeneous while ultra-relativistic (as expected), then the suppression factor is Ω b / ( Ω b + Ω c ) ≈ 0.16 assuming adiabatic initial conditions.
The hierarchical formation of galaxies is shown in figure 6 of Reference [
Dark matter with mass ≈ 69 eV also has tension with studies of the power spectra of density perturbations on scales less than 10 h−1 Mpc with the “Lyman-alpha forest” [
We conclude that the rotation curves of regular spiral galaxies are described, within observational uncertainties, by Equations (3) to (6). These equations describe the stationary state of two self-gravitating non-relativistic ideal gases, “baryons” and “dark matter”, that do not interact significantly with each other except for gravity. These equations can be integrated numerically given four parameters (that are boundary conditions). These four parameters are obtained by minimizing a χ 2 between the measured and calculated spiral galaxy rotation curves.
From these studies, we obtain a wide range of quantitative results on the physics of galaxies, dark matter, and cosmology. In particular, we have presented precision measurements of the dark matter mass m h and temperature T h ( a ) .
The most significant result is the ratio of the temperatures of dark matter and photons after e + e − annihilation while dark matter is still ultra-relativistic:
T h / T = 0.389 ± 0.008 ( stat ) ± 0.058 ( syst ) . (28)
The plausibility that this result is correct lies in the consistency of nine complete and independent measurements of the mass m h of particles of dark matter with spiral galaxies that span three orders of magnitude in mass, and five orders of magnitude in angular momentum. The ratio (28) could have been orders of magnitude larger than, or smaller than, unity, because T depends on T 0 , while T h is independent of T 0 . The ratio T h / T depends on the measured adiabatic invariant 〈 v r h 2 〉 1 / 2 ρ h ( r min ) − 1 / 3 , so this analysis provides evidence that dark matter was once in thermal equilibrium with the (pre?) Standard Model particles. It is therefore likely that dark matter particles will have one of the precise masses listed in
The present studies pose new questions. Why is the spiral galaxy core radius similar in galaxies that span three orders of magnitude in mass? It seems likely that this core radius is of cosmological origin. In fact Equation (18) is already a constraint. It also seems likely that the angular momentum of spiral galaxies is of cosmological origin. Studies of angular momentum correlation between galaxies should settle this issue. If the angular momenta are cosmological, we have a handle on the amplitude of the primordial vector modes. The importance of this avenue of research lies in the fact that the velocities of these vector modes grow towards the past conserving angular momentum and hence could dominate the origin of the Universe.
The author declares no conflicts of interest regarding the publication of this paper.
Hoeneisen, B. (2019) A Study of Dark Matter with Spiral Galaxy Rotation Curves. International Journal of Astronomy and Astrophysics, 9, 71-96. https://doi.org/10.4236/ijaa.2019.92007