1. Introduction
The formation of galaxies is qualitatively different if dark matter is warm instead of cold. The cold dark matter limit is singular. It turns out that understanding galaxies, and the formation of galaxies, is less difficult if dark matter is warm. So we add to the cold dark matter ΛCDM cosmology one more parameter, namely the temperature-to-mass ratio of dark matter, and let observations decide whether dark matter is warm or cold. These lecture notes do bring new understanding of dark matter, and constrain the properties of the dark matter particles.
There is an extensive literature on warm dark matter, e.g. [1] and references therein. Here, our main focus of attention is the understanding of a relation between observables and an adiabatic invariant in the early universe, which is either an unlikely coincidence, or the solution to the long-standing dark matter problem.
2. Warm Dark Matter
Let us consider warm dark matter as a non-relativistic classical noble gas of particles of mass m, density ρ and temperature T. By “classical” we mean that the velocity distribution of the non-relativistic particles of dark matter in the early universe is assumed to be the Maxwell distribution, i.e. is not degenerate. For a discussion on how dark matter might have acquired the Maxwell distribution of velocities, see [2] . By “noble” we mean that collisions, if any, do not excite internal
states of the particles. Recall that for a non-relativistic gas,
.
is the root-mean-square thermal velocity of the dark matter particles. Recall that for adiabatic expansion of a noble gas,
(here V is the volume), so
(1)
is an adiabatic invariant. We will review these concepts in more detail in the next lecture. We want to stress that Equation (1) is valid for a collisionless or collisional gas, and is valid whether, or not, the particles are bouncing off the walls of an expanding box of volume V.
To verify these statements, let us now consider a free particle in a homogeneous universe with expansion parameter
, normalized to
at the present time
. The particle has velocity
at the space point
. In time
the particle advances
and arrives at the space point
. Due to the expansion of the universe,
moves away from
with velocity
(2)
So the velocity of the free particle relative to
, when it coincides with
, is reduced by
, so
(3)
Note that this expression is in agreement with Equation (1). As the universe expands, velocities decrease in proportion to
, the temperature decreases in proportion to
, and the density decreases in proportion to
.
It is not necessary to invoke General Relativity. The preceding arguments are valid also in the non-relativistic physics of Newton. During galaxy formation, the dark matter particles have a velocity field in addition to the thermal velocity. Whenever the velocity field diverges, dark matter cools as described in the preceding paragraph.
We define the “adiabatic invariant” of warm dark matter as the comoving root-mean-square thermal velocity of the dark matter particles:
(4)
is the mean dark matter density of the universe at the present time (throughout we use the notation and the values of parameters of [3] ). The sub-index h stands for the dark matter halo. We will use the sub-index b for “baryons”, mostly hydrogen and helium. These sub-indices will be used only when needed.
Consider a free observer in a density peak in the early universe. This observer feels no gravity and “sees” warm dark matter expand adiabatically, reach maximum expansion, and then contract adiabatically into the core of a galaxy. The adiabatic invariant (4) in the early universe remains constant in the core of the galaxy throughout the formation of the galaxy. This statement is non-trivial, so we will study it in detail in the following lectures, and finally will validate it with observations.
To the cold dark matter cosmology ΛCDM, that has six parameters [3] , we add one more parameter, namely the adiabatic invariant
, and obtain the warm dark matter cosmology ΛWDM. As we shall learn in the following lectures, we are able to measure
and
in the core of spiral galaxies, and will therefore be able to obtain
, and then decide whether dark matter is warm or cold.
3. The Exponential Isothermal Atmosphere
This lecture is included to remind the reader of results we will be using later on. (For more background, I recommend studying the awe inspiring Feynman Lectures on “The exponential atmosphere” [4] .)
James Clerk Maxwell presented several clever arguments to obtain the distribution of velocities of the particles in a gas in thermal equilibrium. The number of particles per unit phase space volume
is proportional to
:
(5)
where
is the particle momentum. The momenta
are assumed isotropic (later on we will lift this assumption when needed). This Maxwell distribution has been validated experimentally (which settles the issues of the clever arguments), and is our point of departure in these lectures. Using the definite integrals
(6)
the following well known results are readily obtained:
(7)
where the pressure P is defined as twice the momenta
of particles with
traversing unit area in unit time.
Let N/V be the number of particles per unit volume of the gas. Then the normalized Maxwell distribution is
(8)
The number of particles with z-momenta in the interval
to
, traversing unit area in unit time, obtained from (8), is
(9)
Now consider gas in equilibrium in a column in a uniform gravitational field
with
. The potential energy of a particle at altitude
is
. We wish to obtain the distribution of momenta at altitude h. We consider particles at altitude 0 with
and
, with z-momenta in the interval
to
. These particles reach altitude h with z-momenta in the interval
to
, where
(10)
Subtracting these two equations, and keeping first order terms, obtains
(11)
Multiplying by dt obtains
(12)
This equation is Liouville’s theorem in one dimension: the volume of phase space occupied by a set of particles moving in a conservative field is a constant of the motion.
In the steady state, the flux
of particles at altitude h with z-momenta in the interval
to
is the same as the flux
of particles at altitude 0 with z-momenta in the interval
to
, since these are the same particles. Therefore, expressing
in terms of variables at altitude h, we obtain:
(13)
where
(14)
So, the density varies with altitude as
(15)
Comparing (13) with (9), we obtain the distribution of particles at any altitude:
(16)
where the energy of a particle is
(17)
Equation (16) is the Boltzmann distribution. The proportionality constant in (16) is independent of altitude. Note that T in (16) is the same T as in (8), so the column of gas is isothermal! These equations are valid for general potential energies
not explicitly dependent on t.
An amazing result of these calculations is that the root-mean-square velocity
of the particles of the gas is independent of altitude! This is because only the more energetic particles with the Maxwell distribution reach a higher altitude. Another amazing result is that the column of gas is isothermal because of the Maxwell distribution of momenta, not because of thermal contact, or not, with the walls (if any) of the gas column.
We have not considered particle collisions. It turns out that results are unchanged because collisions conserve energy (another non-trivial statement).
If the particles are collisionless, the thermal velocities may not be isotropic, in which case we will consider separately each component of the thermal velocity.
One more equation that we will be using in the sequel is
(18)
in agreement with (15). This equation expresses that the difference of pressure
supports the weight of the gas between z and
. This result may again be surprising since there are no membranes at z and
on which the particles can bounce off. Equation (18) expresses conservation of momentum.
4. The Cored Isothermal Sphere
Let us repeat the arguments of the preceding lecture, but this time we consider a self-gravitating gas with spherical symmetry. The equations to be solved are
(19)
The rotation velocity
of a test particle in a circular orbit of radius r is given by
(20)
First we seek a particular solution of the form
, with
independent of r, i.e. we limit the scope of the present lectures to the isothermal case. There is a single well known solution of this form (with
):
(21)
The potential energy difference between particles at r and a reference rc is
, so (15) is valid. Note that to obtain a solution it is necessary to provide
.
The total mass of the halo can be defined as
with
. Then
, which is the Tully-Fisher relation if
is proportional to the absolute luminosity.
Let us review Liouville’s theorem for the problem at hand. Let dn be the number of particles occupying the phase-space volume
. If the force field is conservative, i.e. if (17) is satisfied with
not explicitly dependent on t, then, in the course of the motion of the dn particles, the corresponding volume of phase-space moves in phase-space, and may change shape, but its volume turns out to be a constant of the motion. Also, no particles cross the boundary of this phase-space volume, so dn is also a constant of the motion. The phase-space density
at time t is equal to the phase-space density
at time
, and so can only be a function of constants of the motion, e.g. the energy E. If this function of E is exponential the system is isothermal, i.e.
is independent of
. In the present case of the isothermal sphere, E does not depend on the directions of
or
, so the phase-space volume can be written as
, with
, and, in particular,
. If convenient, the phase-space volume can also be written as
, where
is the momentum transverse to
. Then the flux at r per unit area and time of particles with
in
, and similarly at
, are
(22)
(23)
Note that
, so
.
Equations (19) can be solved numerically in radial steps dr, starting from
. To start the numerical integration we need to provide two boundary conditions, e.g.
and
, in addition to
. For solutions with no black hole at
we can take
. For these isothermal spheres with a core the solution for
is (21), while the solution for
is
(24)
The two asymptotes meet at a “core radius”
(25)
We note that the adiabatic invariant (1) is a minimum at
where
is a minimum. Since the cored isothermal sphere has formed from a density perturbation in the early universe, its two parameters
and
are related by the adiabatic invariant (4):
(26)
Further justification of (26), and observational validation, will be given in following lectures. Therefore, in the present work, the cored isothermal sphere, and its Boltzmann distribution of dark matter particles, is defined by a single independent parameter,
or
. Equation (26) assumes isotropic velocities (valid in the core due to spherical symmetry, and also valid if the dark matter halo is isothermal).
We note that a measurement of
obtains
, and a measurement of
obtains
, and together they obtain the adiabatic invariant
.
5. An Example
As an extreme example, with a very large
, let us consider the spiral galaxy UGC11914 (note that our final measurements of
are obtained from rotation curves of dwarf galaxies [6] ). The rotation curves of UGC11914, observed by the SPARC collaboration [5] , are presented in Figure 1. We fit these rotation curves by solving the equations:
Figure 1. Top: Observed rotation curve
(dots) and the baryon contribution
(triangles) of galaxy UGC11914 [5] . The solid lines are obtained by numerical integration as explained in the text. Bottom: Mass densities of baryons and dark matter obtained by the numerical integration.
(27)
(28)
(29)
(30)
and
are the contributions of baryons and dark matter to the gravitation field
. We have included
and
to account for rotation. These parameters are quite uncertain, though a rough estimate is
and
[7] or less [8] . We can eliminate
and
from the numerical integration by replacing
, and
. We start the numerical integration at the first measured point at
, and end at the last measured point, in steps dr. To start the numerical integration we need six boundary conditions:
,
,
,
,
and
, where
(31)
(32)
to include the possibility of a black hole at the center. Good fits are obtained assuming
and
are independent of r. We vary the six boundary conditions to minimize the
between the observed and calculated rotation velocities, and obtain:
(33)
(34)
(35)
Uncertainties are statistical from the fit. The core radius is 0.70 kpc. The core dark matter density
is 2.6 × 108 times the mean dark matter density of the universe
. We note that this galaxy indeed has a black hole at its center. From
and
we obtain the adiabatic invariant (26):
(36)
where the uncertainty is statistical only.
Let us mention that similar results are obtained from galaxies spanning 3.5 orders of magnitude in absolute luminosity [6] , validating that the adiabatic invariant in the core of galaxies is conserved, confirming that
is of cosmological origin, and that dark matter is indeed warm!
If we add an r-dependence to
we generally obtain a higher
of the fit, so indeed, within uncertainties, the dark matter halo is isothermal, at least out to the observed rotation curves, indicating that the particles have the Maxwell-Boltzmann momentum distribution. We note that
is of the same order of magnitude as
, because dark matter and baryons fall into the same potential well. Since thermal equilibrium implies
, we conclude that, unless
, dark
matter is not in thermal equilibrium with baryons (mostly hydrogen atoms). So, on galactic scales, we can neglect dark matter-baryon interactions. If dark matter feels only the gravitational interaction, it can be shown that deviations of its trajectory due to particle-particle encounters, or interchange of energy with baryons, can be neglected on galactic scales. We also note that the ratio of dark matter to baryon densities in the galaxy is generally of the order of the universe average.
The galaxy density (21) reaches
at
(37)
The age of the universe at, say, redshift
, is 3 × 1016 s. In this time a particle with constant velocity 198 km/s reaches 0.2 Mpc, less than
, and much farther than the last observed rotation velocity.
From (36), neglecting
, we estimate that dark matter becomes non-relativistic at expansion parameter
(38)
i.e. after
annihilation, and while the universe is still dominated by radiation [3] .
Now let us do the following back-of-the-envelope approximate calculations. An ultra-relativistic gas with zero chemical potential has a number density of particles
[3] . At expansion parameter
the temperature is
. The present number density is
, if dark matter particles do not decay or annihilate at
, i.e. if there is no “freeze-out”. From these equations we obtain
, (39)
or
. (40)
For scalar dark matter, i.e.
and
, and neglecting
, we estimate the mass of the dark matter particles from (36):
. At expansion parameter
the photon temperature is
. From the preceding equations we obtain
(41)
or, for our example,
. That this ratio is of order 1, given that the dark matter particle mass
is uncertain over 89 orders of magnitude [3] , is surely telling us something! Furthermore, note that dark matter is sufficiently cooler than photons to (marginally) evade the “thermal relic” mass limits obtained from the Lyman-α forest [9] , and sufficiently cool to not spoil the success of Big-Bang Nucleosynthesis [7] . Finally, within experimental uncertainties, dark matter is in thermal and diffusive equilibrium with the Standard Model sector, and decouples from this sector, at T somewhere between the top quark mass
and the temperature
of the deconfinement-confinement transition from quarks to hadrons (decoupling at a lower temperature compromises the agreement with Big Bang Nucleosynthesis [7] ). For a proper treatment of the preceding estimates see [2] [10] . It turns out that Equations (40) and (41) change by less than 1%. For a summary of measurements and their interpretation, see [11] . For details of each measurement see the citations in [9] .
Excersize to test the Boltzmann factor
: From the measured dark matter density in the vecinity of the sun [3] , and other Milky Way data that you can find in the literature, estimate
. Is your answer consistent with (36)?
Exersizes: The isothermal assumption can be lifted by using the observed galaxy density runs
and
[12] [13] . Study these articles and obtain the adiabatic invariant
from each of them. Are your answers consistent with (36)?
6. Galaxy Formation
Let us solve, by numerical integration, the following equations that describe the formation of a galaxy: Newton’s equation
(42)
the continuity equations
(43)
(44)
and Euler’s equations
(45)
(46)
is the gravitation field,
is the proper (not comoving) coordinate vector,
and
are the mass densities,
and
are
the velocity fields, and
and
are the radial (1-dimensional)
velocity dispersions, i.e. thermal velocities, of dark matter and baryons (mostly hydrogen), respectively. Equations (45) and (46) express the conservation of momentum. The static limits of Equations (42) - (46) are Equations (27) - (30). These equations need to be supplemented by the equations of state of dark matter and baryons at
in order to obtain a solution. The Maxwell distributions of baryons and dark matter have different temperatures
and
, respectively. The dark matter velocity field converging onto the core has two consequences: it increases both the core density
and temperature
, maintaining constant the dark matter adiabatic invariant (26) (as explained in lecture 2), and similarly for baryons. When the dark matter halo reaches a final stationary state its cored isothermal halo is defined by a single parameter
.
An example, discussed in [15] , is presented in Figure 2. For baryons we take
, corresponding to hydrogen decoupling from photons at
[16] . We set the initial
and
as shown in Figure 2. We integrate the equations in steps dt, and for each t, in steps dr, starting at
, to calculate the new
and
. The above equations need to be supplemented by the adiabatic conditions, so, for each step of t, we set, only once, i.e. for all r:
, (47)
and
. (48)
Equation (47) is justified out to at least the last observed rotation velocity by the nearly flat rotation curves. The isothermal prescription (47) cannot be correct beyond
where the universe is expanding homogeneously. However, this is not a problem since, beyond
,
and
are zero. While the hydrogen and helium gas remains adiabatic, i.e. until excitations, ionization, radiation, shocks and star formation become significant, we require (48). Figure 2 illustrates the formation of a core, and how the galaxy reaches a stationary isothermal state starting from
and progressing outwards. For additional examples, see [14] and [15] .
Note that, independently of the details of how the core reaches its final density
, including spherical or filament collapse and galaxy mergers, this density fully specifies the cored isothermal sphere with its temperature, core radius, Boltzmann distribution, and adiabatic invariant
(up to a rotation, and perhaps also a relaxation, correction ≤ 1, that observationally is found to be in the approximate range 1 to 0.5 [6] ).
Finally, let us mention that in the cold dark matter scenario, i.e. in the limit
, we are left with no hydrostatic Equations (27) - (30), or hydrodynamical Equations (42) - (46).
Figure 2. Shown is the formation of a warm dark matter plus baryon galaxy. The densities
and
, and velocity fields
and
, are presented at time-steps that increase by factors 1.4086 (or
for the dot-dashed lines). The initial density perturbations are Gaussian, with parameters listed in the figure. In this example, dark matter is warm with
.
.
is increased from 10−4 to 1 keeping constant baryon angular momentum. Note the initial velocity field corresponding to the expansion of the universe. Note that both warm dark matter and baryons acquire cores. The baryon core is due to the baryon angular momentum conservation. The warm dark matter core is due to the adiabatic invariant feedback (47). Note the effect of this feedback on
at the final times shown (studied in [14] ). An isothermal steady state is reached for both warm dark matter and baryons, starting at
and then extending out to larger r, at least out to the observed rotation curves. This Figure is taken from [15] .
7. Conclusions
We have studied galaxies, and galaxy formation, assuming that dark matter is warm. This scenario requires the addition of one parameter to the cold dark matter ΛCDM cosmology, namely the adiabatic invariant
defined in (4). We find that the formation of galaxies, all the way from linear perturbations in the early universe, until the galaxies reach an isothermal stationary state, conserves the adiabatic invariant
in the core to a good approximation. The arguments in favor of this view are not without objections, so observational validation is necessary. This adiabatic invariant feedback and the initial conditions determine the final core density
, and hence the entire cored isothermal sphere including its core radius and temperature (at least out to the observed spiral galaxy rotation curves). Note that it is not necessary to invoke dark matter self interactions, baryonic feedback, or fermionic dark matter quantum degeneracy, to justify a core (though these phenomena may influence the core). The observed spiral galaxy rotation curves allow measurements of
(up to a rotation, and perhaps also a relaxation, correction ≤ 1, that observationally is found to be in the approximate range 1 to 0.5 [6] ). These measurements are consistent for galaxies with absolute luminosities spanning 3.5 orders of magnitude [6] , so the analysis is validated by observations, and the interpretation that
is of cosmological origin is confirmed. Independent measurements of
are obtained by studying the consequences of the warm dark matter free-streaming suppression factor
of the cold dark matter power spectrum of density perturbations, i.e. galaxy stellar mass distributions, galaxy rest-frame ultra-violet luminosity distributions, first galaxies, and their effect on the reionization optical depth [17] . All of these measurements of
are consistent, solve the small scale problems of ΛCDM [18] , and constrain the warm dark matter particle properties [9] [11] [19] .