1. Introduction
The Standard Model of quarks and leptons is enormously successful, it has passed many precision tests, and is here to stay. However, if the Standard Model were complete, the universe would have no matter: no dark matter, little baryonic matter, and no neutrino masses. “The New Minimal Standard Model” [1] is an extension that aims to “include the minimal number of new degrees of freedom to accommodate convincing (e.g., >5σ) evidence for physics beyond the Minimal Standard Model”. But this aim has a moving target: as new data becomes available, the model may need to be amended accordingly. The inclusion of a “minimal number of new degrees of freedom” is in accordance with the absence of new particles at the LHC. The purpose of the present study is to see if the New Minimal Standard Model is consistent with the new data on dark matter that has recently become available, and, if necessary, update the model accordingly.
Let us briefly describe the New Minimal Standard Model [1]. First, the Standard Model Lagrangian is extended to include classical gravity. Next, a gauge singlet real scalar Klein-Gordon field with
parity is added for dark matter. Dark energy is described by the cosmological constant
. Two gauge singlet Majorana neutrinos are added to account for neutrino masses and mixing (leaving one neutrino massless until data requires otherwise), and also to obtain baryogenesis via leptogenesis. Finally, a real gauge singlet scalar field is included to implement inflation.
The outline of this article is as follows. Measurements of dark matter properties are presented in Section 2. Scalar, vector and sterile neutrino dark matter models are studied in Sections 3 to 5. We close with conclusions.
2. Measured Properties of Dark Matter
Fits to spiral galaxy rotation curves [2] [3] [4], and studies of galaxy stellar mass distributions [5] [6] [7], independently obtain the following dark matter scenario. Dark matter is in thermal and diffusive equilibrium with the Standard Model sector in the early universe, i.e. no freeze-in, and decouples (from the Standard Model sector and from self-annihilation) while still ultra-relativistic, i.e. no freeze-out. The decoupling occurs at a temperature
to not upset Big Bang Nucleosynthesis. Dark matter has zero chemical potential. The root-mean-square velocity of non-relativistic dark matter particles, at expansion parameter
, is
, where
. Dark matter becomes non-relativistic at an expansion parameter
. Dark matter is warm with a free-streaming cut-off wavenumber
. The corresponding free-streaming transition mass is
, comparable with the mass of the Milky Way. The dark matter particle mass is
(
) for scalar (vector) dark matter. There is evidence in favor of boson dark matter with a significance of 3.5σ [7]. The number of boson degrees of freedom is limited to
or 2, i.e. to scalar or vector dark matter. The ultra-relativistic dark matter temperature, relative to the photon temperature, after e+e− annihilation, is measured to be
(
) for scalar (vector) dark matter. All uncertainties have 68% confidence. These numbers are obtained from Table 4 of [7] for the boson scenario that assumes that non-relativistic dark matter particles reach non-relativistic thermal equilibrium (NRTE) (i.e. the non-relativistic Bose-Einstein momentum distribution) due to their dark matter-dark matter elastic scatterings. The relations between
and
and
for zero chemical potential are [7]:
(1)
(2)
In the case of negligible dark matter elastic scattering, the non-relativistic dark matter retains its ultra-relativistic thermal equilibrium (URTE), i.e. the ultra-relativistic Bose-Einstein momentum distribution, and the measurements are
,
,
,
, and
(
) and
(
) for scalar (vector) dark matter. The relations between
and
and
for zero chemical potential are [7]:
(3)
(4)
For an overview of these measurements see [8]. To make this article self-contained, Figure 1 presents forty-six independent measurements of
from fits to spiral galaxy rotation curves [4]. From
we calculate the warm dark matter free-streaming cut-off wavenumber
[7]. This cut-off wavenumber is also obtained from galaxy stellar mass distributions as shown in Figure 2 [7]. These independent measurements are consistent!
The current limit on dark matter self interaction cross-section is
with 95% confidence [14] [15]. A tentative measurement obtains
[16]. If this measurement is confirmed, dark matter retains URTE.
The current limits on dark matter particle mass are
for fermions, and
for bosons [14].
In the present study we will assume this specific dark matter scenario, and ask the following questions. What dark matter interactions lead to this scenario? How is dark matter created? How does dark matter and the Standard Model sector come into thermal and diffusive equilibrium? How do they decouple?
Figure 1. Forty-six independent measurements of the expansion parameter
at which dark matter particles become non-relativistic (uncorrected for dark matter halo rotation). Each measurement was obtained by fitting the rotation curves of a spiral galaxy in the Spitzer Photometry and Accurate Rotation Curves (SPARC) sample [9] with the indicated total luminosity at 3.6 μm. Full details of each fit are presented in [4].
Figure 2. Distribution of stellar masses of galaxies at redshift
compared with predictions. From this data, and similar distributions corresponding to
, and 8, we obtain the power spectrum cut-off wavenumber
. Figure from [7]. The data are from [10] [11] [12] [13].
How does dark matter acquire mass? Why is dark matter stable (relative to the age of the universe)? And, why is the measured dark matter particle mass
so tiny compared to the Higgs boson mass
?
Notes: For a discussion of tensions between measurements of, and limits on, thermal relic dark matter mass see [7] [8]. We should mention that the observed galaxy mass distribution presented in Figure 2 is in tension with Lyman-α forest studies [17]. The 3.5σ confidence in favor of boson dark matter mentioned above, based on spiral galaxy rotation curves and galaxy stellar mass distributions, does not include the Tremaine-Gunn limit on fermion dark matter mass [18] [19]. Including this limit would strengthen the confidence. However, the Tremaine-Gunn limit needs to be revised in view of resent observations on dwarf spheroidal “satellites” of the Milky Way [20] [21] [22] [23].
3. Scalar Dark Matter
The measured dark matter properties allow scalar or vector dark matter, with fermion dark matter disfavored but not ruled out. We begin with the real scalar field
of [1]. To attain thermal and diffusive equilibrium between dark matter and the Standard Model sector we need to add a coupling between the two. The simplest renormalizable coupling is proportional to
since
is the only Standard Model gauge singlet scalar with mass dimension ≤ 2.
is the Higgs boson field. The interaction rates
, relative to the universe expansion rate, scale as 1/T, so equilibrium is approached towards the future, and statistical equilibrium needs to be achieved by
to avoid freeze-in. Decoupling occurs when the Higgs boson
becomes non-relativistic at
. Thereafter the reaction rates become exponentially suppressed because the Higgs bosons annihilate, and only the tail of the
particle momentum distribution is above threshold. With
there is no freeze-out if
is stable. A super-renormalizable interaction proportional to
needs to be avoided because it leads to a ratio of number densities
that depends on T. For this reason, and to obtain a stable
, and to avoid extra parameters in the potential
, we impose a
symmetry
.
Therefore, we consider a gauge singlet real Klein-Gordon scalar dark matter field
, with
symmetry
, and portal coupling to the Higgs boson [1]. Here we present a brief review of the model to see if it can describe the observed properties of dark matter. To the Standard Model Lagrangian we add
(5)
and a contact coupling to the Higgs field
:
(6)
(We are omitting the metric factor
.) After electroweak symmetry breaking (EWSB) the Higgs doublet, in the unitary gauge, has the form
(7)
with real
, the interaction Lagrangian becomes
(8)
and dark matter particles acquire a mass squared
(9)
assumed to be >0. We note that
is absolutely stable since there is no interaction term with a single
.
The running of coupling parameters to 1-loop or 2-loop order can be found in [24] [25] [26] [27]. Some center of mass cross-sections are
(10)
(11)
where
is the Mandelstam variable. The reaction rates are exponentially suppressed at
or
. These interactions bring dark matter into thermal and diffusive equilibrium with the Standard Model sector at
if
and 10−6, respectively. The Higgs boson invisible decay rate for
is
(12)
Requiring this decay rate to be less than the limit on the invisible width of the Higgs boson (≈0.013 GeV [14] ) implies
. In summary, we require
.
As an example, take
and
, so
. Then there is fine tuning in (9):
. Note that to achieve
as low as 73 eV starting from
requires fine tuning between two unrelated input parameters with dimensions of mass.
Let us now check whether non-relativistic dark matter acquires the non-relativistic Bose-Einstein momentum distribution due to elastic scattering. The cross-section at
(neglecting interference with (14)),
(13)
implies that the mean time between collisions of dark matter particles at
is less than the age of the universe even for
, so, in this model, non-relativistic dark matter has non-relativistic thermal equilibrium. The cross-section (neglecting interference with (13)),
(14)
also corresponds to collisional dark matter if
.
Dark matter decouples from the Standard Model sector at
when the Higgs bosons become non-relativistic. As the universe expands and cools, particles and antiparticles that become non-relativistic annihilate heating the Standard Model sector without heating dark matter, or neutrinos if they have already decoupled. For decoupling at
we expect the temperature of ultra-relativistic dark matter, relative to the photon temperature, after e+e− annihilation, to be
[14], which can be compared with the measured ratio
[7].
The cross-section limit
[14] at
, and (13), implies
, so the present model is ruled out. If we lower
to this value,
and the Standard Model sector do not achieve statistical equilibrium at
.
4. Vector Dark Matter
To reduce the dark matter-dark matter elastic scattering cross-section, and to relieve the fine tuning in the model of Section 3, we attempt reaching the small
in two steps.
To the Standard Model Lagrangian we add a complex scalar field
that is invariant with respect to the local
transformation
. The corresponding vector gauge boson
acquires mass due to the breaking of the
symmetry of the ground state. In the present model,
is the dark matter candidate, and
decays to
. The dark matter sector is known in the literature as the “Abelian Higgs model”.
The relevant part of the Lagrangian is
(15)
(16)
(17)
(18)
(19)
and the Standard Model sector have no charges in common.
For
,
,
, and
, there is symmetry breaking, and the fields
and
acquire vacuum expectation values [28]
(20)
if
and
. In unitary gauge, the real amplitudes A and
become the longitudinal components of
and
, respectively, and the complex amplitude
becomes the longitudinal components of
and
. The mass eigenstates are
(21)
(22)
and the mixing angle is
(23)
To bring
into thermal and diffusive equilibrium with the Standard Model sector without exceeding the limit on the invisible width of the Higgs boson
requires
as in Section 3. Some reactions of interest are
(24)
(25)
(26)
The couplings of
,
, and
(up to order 4) are proportional to
,
,
,
,
,
,
,
, and
. The kinematics allow
to decay only to γ’s or ν’s. However, we note that there is no coupling with a single
, so
is absolutely stable.
Our challenge is to choose parameters so that
attains statistical equilibrium with the Standard Model sector at
, but
does not; and we need the decay
to occur after
has decoupled from the Standard Model sector, and while
is still ultra-relativistic, i.e. within the temperature range
; and that
.
Case
: Let us assign the high mass eigenstate to
, and the low mass eigenstate to
(the opposite case will be considered below). A particular solution of interest has
, so
as in the Standard Model. A bench-mark scenario with
is
,
,
,
,
, and
. To meet all requirements, there is fine tuning of
to lower
: the relative difference of the two terms in the numerator of (20) is ≈10−6.
The reaction rate of
, relative to the expansion rate of the universe
, is
at
, so this coupling is strong. For
,
at
, so this coupling is also strong. For
,
(200) at
(
), so
does not attain statistical equilibrium with
, or with the Standard Model sector, at
. The decay rate of
, relative to the expansion rate of the universe
, is
(3 × 104) at
(
), so indeed we have arranged that the decay occurs after
has decoupled, and while
is still ultra-relativistic, i.e. in the temperature range
.
The cross-section for
at
is
(27)
This cross-section implies that the mean dark matter particle interaction rate is much less than the expansion rate of the universe
at all temperatures, so, in this model, non-relativistic dark matter retains the ultra-relativistic Bose-Einstein momentum distribution.
The two
’s in the decay
have correlated polarizations, so the average number of boson degrees of freedom, needed to calculate the dark matter density (see (21) of [7] ) is
. Then, from (3) and (4), the measured values for this scenario are
, and
.
For zero chemical potential, the number of
per unit volume, given by the ultra-relativistic Bose-Einstein distribution, is
(28)
where the number of boson degrees of freedom of
is
. After the decay
(29)
Each
in 8 orbitals of momentum 2p decays to two
’s corresponding to one orbital with momentum p, so
(30)
Integrating, we obtain
. So, the predicted ratio is
, to be compared with the measured value
.
The cross-section limit
[14] at
, and (27), implies
, in agreement with the benchmark solution. The tentative measurement
[16], if confirmed, would imply
, which is in agreement with the benchmark solution within uncertainties!
In summary, the vector model with
is consistent with all currently measured properties of dark matter. There is fine tuning to obtain the small required symmetry breaking of the ground state of
.
Case
: Let us now assign the high mass eigenstate to
, and the low mass eigenstate to
. Again, as an example, we consider the case
, so
as in the Standard Model, and
. A benchmark solution with
is
,
,
,
,
, and
. When particles
become non-relativistic at
, they decay mostly to the Standard Model sector: reactions
are much faster than the universe expansion rate, while
and
are much slower, so the universe is left with no dark matter.
Assigning charges QS to Standard Model particles, to enhance or replace the contact interaction between
and
, does not lead to compelling alternative models.
5. Sterile Neutrino Dark Matter
Observations of spiral galaxy rotation curves and of galaxy stellar mass distributions favor boson over fermion dark matter with a significance of 3.5σ [7], so we should not yet rule out fermion dark matter. Sterile neutrinos have been studied extensively as dark matter candidates [29] [30] [31]. In this section we briefly review sterile neutrinos and see if they are consistent with the measured properties of dark matter presented in Section 2.
We extend the Standard Model with a gauge singlet neutrino
with a Majorana mass
. This is the measured mass for the case of fermion dark matter retaining ultra-relativistic thermal equilibrium (URTE), see Table 4 of [7]. We will refer to the two irreducible representations of the proper Lorentz group of dimension 2 as “Weyl_L” and “Weyl_R”. For simplicity we focus on one generation.
and
are two-component Weyl_L and Weyl_R fields, respectively.
and
transform as Weyl_R and Weyl_L fields, respectively, where
is a Pauli matrix.
,
,
, and
are scalars with respect to the proper Lorentz group.
To include Weyl spinors into the Standard Model, it is convenient to use 4-component Dirac spinor notation. Our metric is
. The matrices A and C are defined, in any basis, as
, and
[31], with
,
, and
. We define
, and the charge conjugate field
. Then
, and
. A Dirac spinor that satisfies
is a Majorana spinor (
is an arbitrary phase).
In a Weyl basis [30],
,
,
(31)
(32)
(33)
Note that
, and
, so these are Majorana fields. With this notation the Majorana fields
and
can mix. Note however that
and
are distinct:
has weak interactions while
does not.
,
,
,
, and
are scalars with respect to the proper Lorentz group. The neutrino mass term after electroweak symmetry breaking has the form [31]
(34)
where
is a Dirac mass (
is a Yukawa coupling), and
is a Majorana mass. We consider the case
. The mass eigenstates are [31]
(35)
where
,
, and
.
Let us now consider dark matter production. We are interested in the reactions
, or
. First, we verify that the produced
is a coherent superposition of
and
. The coherence factor is [32]
(36)
with
. Since we are interested in energy E of
of order
, we take its uncertainty to be
, so the first factor is 1.
is the mean time between
interactions. The propagation time of
and
over which their wave packets cease to overlap is the decoherence time [32]
(37)
Taking the wave packet duration
, we estimate
for the small value of
being considered. In conclusion,
, and
and
do not become decoherent between
interactions, so we must take into account their oscillations.
Consider initial conditions for
production to be
and
. Then, from (35), we obtain the probabilities
to observe a Weyl_L neutrino, and
to create a sterile neutrino,
(38)
with
, and
. Note that
for
, but this probability is suppressed by a factor
for small
. Equation (38) describes the oscillation between the active and sterile neutrinos. Similar phenomenology has been confirmed in neutrino flavor oscillation experiments.
The cross-section
is given by Eq. (50.25) of [14]. Multiplying by
we obtain the cross-section for sterile neutrino production
. We find that the production mechanism
to bring
into statistical equilibrium with the Standard Model sector at
, and decouple at
, fails because of the interference factor
. (Note: In Figure 11 of [2] I did not include this factor so that figure is wrong.)
The production channel
is negligible.
6. Conclusions
Accurate, detailed and redundant measurements of dark matter properties have recently become available [7]. We have studied scalar, vector and sterile neutrino dark matter models in the light of these measurements. The vector dark matter model presented in Section 4 is (arguably) the renormalizable model with the least number of new degrees of freedom that is consistent with all current observations, and replaces the scalar dark matter model of Section 3 [1] that is ruled out. The sterile neutrino dark matter production mechanism studied in Section 5 did not meet experimental constraints.
New insights pose new questions. If nature has chosen the vector dark matter of Section 4, why do the two terms in the numerator of (20) cancel to 1 part in 106? Similar questions can be made regarding the cosmological constant
, or the strong CP phase
. Do the scalars
and/or
participate/cause inflation? Baryogenesis via leptogenesis (arguably) requires sterile Majorana neutrinos. How are they produced? What is the origin, if any, of their masses?
How can we move forward? A signal in direct dark matter searches would rule out the vector model. Indirect searches may find an excess of photons (or neutrinos!) with energy ≈36 eV, ≈53 eV, or ≈62 eV, if dark matter is unstable and decays. Such a signal would also rule out the vector dark matter model. Collider experiments may discover an invisible Higgs decay width. Further progress will come from the cosmos: more studies of disk galaxy rotation curves, and galaxy stellar mass distributions (these studies can enhance the boson/fermion discrimination, and perhaps can observe the predicted tail of the boson warm dark matter power spectrum cut-off factor
[7] ), galaxy formation simulations, the “small scale crisis” (missing satellites, too big to fail, galaxy core vs. cusp, large voids), super massive black holes at galaxy centers (Einstein condensation may occur at the galaxy center), revised constraints on fermion dark matter mass from the Tremaine-Gunn limit, and tighter constraints on dark matter self-interactions. It is necessary to understand the tensions between the Lyman-α forest studies and the observed galaxy stellar mass distributions, see Figure 2. Studies of dark matter halo rotation in disk galaxies are also needed.