An Introduction to Particle Dark Matter

We review the features of Dark Matter as a particle, presenting some old and new instructive models, and looking for their physical implications in the early universe and in the process of structure formation. We also present a schematic of Dark Matter searches and introduce the most promising candidates to the role of Dark Matter particle.


I. INTRODUCTION
In this review we employ units = c = k B = 1 units, except where otherwise stated.
A. Why do we need Dark Matter?
Cosmic Microwave Background (CMB) observations strongly suggest that the universe is, at large scales, nearly spatially flat [6], giving an observed value for the present energy density ρ 0 close to its critical value: corresponding to an average density of 10 protons per cubic meter; in Eq. (1) H 0 ≈ 70 km/s/Mpc is the value of the Hubble factor today (i.e. the Hubble constant) and G N ≈ 6.67×10 −11 N m 2 /kg 2 is Newton's gravitational constant.
This energy density is made up of various contributions [6].Mostly (about 70%) of Dark Energy (DE), which is responsible for the accelerated expansion of the universe, followed by Dark Matter (DM, about 26%), which is responsible for the gravitational collapse of ordinary matter (usually called "baryons" and which amounts to the remaining 4%) and eventually for the formation of the structures that we observe in the sky.
There are several probes able to weigh the Dark Energy (DE) contribution to this average density: for example, CMB data [6], the distribution of galaxies on large scales, also known as the Large Scale Structure (LSS) of the universe [3] and a byproduct of the latter, the so-called Baryon Acoustic Oscillations (BAO) [7].Moreover, we can also discriminate the contribution of ordinary baryonic matter and DM, via e.g.Big Bang Nucleosynthesis (BBN) and again CMB observations, see e.g.Refs.[30,67].In fact, we are able to determine the amount of DM in the universe, which is the following: where ρ DM,0 is the DM density evaluated today (usually a 0 subscript means, in cosmology, that the corresponding quantity has to be evaluated at the present time, i.e. at the age of the universe t 0 ).The number in Eq. (2) can be rewritten in units more suitable to astrophysics and particle physics as follows: where M ⊙ ≈ 2 × 10 30 kg is the mass of the Sun.
Note that ρ DM,0 is 10 5 times less dense than the average density in clusters of galaxies and 10 6 less dense than galaxies.Thus, we can conclude, from the point of view of perturbations about a homogeneous and isotropic background, that the universe is in a highly nonlinear regime.
Dark Matter is a crucial ingredient, necessary in order to produce the large inhomogeneities that form structures in the universe.Indeed, looking at the CMB sky, we measure very small temperature fluctuations, proportional to the baryonic density fluctuations, which we infer being of order: (at recombination, i.e. for z ≈ 1100) , where δ b is the baryonic density contrast, i.e. the ratio of density fluctuations in the baryonic density and its background value.This number is important because we also know that subhorizon pressureless matter (as baryons are) overdensities in the linear regime grow linearly with the scale factor in the matter-dominated epoch of the evolution of the universe.Therefore, from the recombination epoch at z rec ≈ 1100 baryon fluctuations, which are finally free to grow, have not enough time for entering the nonlinear regime, but they would be of order 10 −1 today, so still linear.This issue is solved if a different "dark" species, which decoupled from photons much earlier than baryons did, is present, giving time for its perturbations to grow enough in order to trigger, at the recombination epoch, a faster collapse of baryonic fluctuations and thus eventually making structure formation possible.
For this reason, it is in general difficult to explain consistently structure formation without recurring to DM. See e.g.[31] for the case against MOdified Newtonian Dynamics (MOND) [61] and its covariant version TeVeS [13].

B. Microscopic features of Dark Matter
We can get some clues about the microscopic nature of DM by trying to answer the following questions: 1. Is Dark Matter actually dark?
It is generally believed that DM is not bright, in the sense that DM particles do not interact electromagnetically and so are not able to scatter photons; in other words they do not possess an electric charge.However, observation does not exclude DM particles to have indeed a tiny electric charge [62].

Is Dark Matter collisionless?
In order for collisional effects in the DM sector to be relevant, the mean free path λ MFP of a DM particle must be larger than the typical size of a cluster of galaxies, i.e. ≈ 1 Mpc, which has a density of order ρ cl ≈ 1 GeV/cm 3 .Therefore, we get, for a DM particle with mass m = 1 GeV, the following result for the cross section σ: which is a cross section whose value is typical of the strong interaction.This suggests that DM might not be completely collisionless and makes models of Self-interacting Dark Matter (SIDM) appealing.It is worth also to mention scenarios of Partially Interacting DM, where only a subdominant component of the total DM contribution is self-interacting [35,36].

Is Dark Matter classical?
We know that DM must be gravitationally bounded on scales at least as large as the size of a dwarf spheroidal (dSph) galaxy.We can use this size in order to derive an upper bound on the de Broglie wavelength λ dB of a DM particle and, consequently, a lower bound on its mass.Indeed, consider a DM particle of mass m with velocity v ≈ 100 km/s.Using the Planck's constant value (which we recover for the moment) h Pl ≈ 4 × 10 −15 eV•s, and the fact that the particle is non-relativistic, the de Broglie wavelength λ dB = h Pl /(mv) is given by: This result means that in order to have λ dB < 1 kpc ≈ 3 × 10 21 cm we need a mass m > 10 −22 eV , which is not a very stringent lower bound.
On the other hand, we know that quantum effects can drastically change such a lower bound e.g. if DM is a fermion, because in this case its phase space density is bounded by the so-called Pauli blocking effect (i.e. two fermions cannot occupy the same quantum state, because of Pauli's exclusion principle).Using the simplest model of a singular isothermal sphere with ρ dSph = σ 2 dSph /(2πG N r 2 ), with r radius of the sphere and σ 2 dSph a constant velocity dispersion one obtains the so-called Tremaine-Gunn limit [75]: where g is the number of spin states of the DM particle.The difference of 23 orders of magnitude in the above two bounds in Eqs. ( 7) and ( 8) for the DM particle mass teaches us how important quantum effects might be in the description of DM.

Is Dark Matter a fluid?
It is possible to describe DM as a fluid?In this case, we can forget about the particle nature of DM and simply describe it by using its density and velocity flow.This is also called fluid approximation.We can get some constraints on the coarse graining of such a fluid, i.e. on the scale down to which the fluid approximation works, by requiring that it must not disrupt e.g.clusters of stars.In fact, if the grains constituting this fluid are too big, e.g.black holes, they must have effects on bound compact systems such as globular clusters.In Ref. [55] it is showed that this would be the case for grains with mass m ≥ 10 6 M ⊙ .
The estimate energy exchanged in collisions among globular clusters and black holes, using the impulse approximation according to which the interactions are instantaneous (as if the bodies were rigid spheres), demand that the energy must be smaller than the gravitational binding energy of the bound compact system (roughly proportional to G N M 2 /R, where M is the mass of the system and R its size), giving a mass of order m 10 3 M ⊙ ≈ 10 70 eV, see e.g.Ref. [39].
In conclusion, from the above phenomenological considerations we can list the following prop-

erties of particle DM:
Mass:: There are more than 90 orders of magnitude of mass range available for bosons; 70 for fermions.
Interactions:: A self-interaction is possible, if of the order at most of the strong interaction, so of the order of the MeV.It is in principle also possible to have interactions with ordinary Standard Model particles, as long as such interactions do not involve emission of photons, otherwise DM halos would shine and be visible.Particularly promising as DM candidates are massive particles that interact only via weak interactions, the so-called WIMPs (Weakly Interacting Massive Particle).
Abundance:: It has to be enough in order to satisfy the observational constraint Ω DM ≈ 0.26.
As we shall see later, this is achieved naturally for WIMPs, and such occurrence is usually referred to as the WIMP miracle.direct detection experiments at present constrain proton-WIMP interactions to be below the level of a cross section of 10 −47 cm 2 for a mass around the weak scale (i.e. about 100 GeV); the absence of emission of gamma-rays from DM-dominated (this DM-domination is inferred of course by gravitational effects) nearby structures such as dSph galaxies implies a pair-annihilation cross section times relative velocity smaller than the required one for thermal production in s-wave annihilation (which means that the thermal averaged cross section is velocity-independent) for masses between around 10 and 100 GeV, with weaker limits at larger masses.Finally, no conclusive result from colliders has been reported, with resulting limits which depend upon the specific WIMP model.
We recommend to the reader the recent review on WIMP searches given in Ref. [10] for further details.

II. THERMAL DECOUPLING
Thermal decoupling provides a successful framework able to explain the origin of species in the early universe.Built on a synergy among Statistical Mechanics, Particle Physics and General Relativity, it allows to make highly accurate testable predictions.
The key idea of thermal decoupling is based on the interplay between the rate Γ of a certain interaction process and the expansion of the universe, described by the Hubble factor H, both timedependent.When Γ ≫ H, the interaction is very efficient and, if we are dealing with annihilation and pair productions, they are in equilibrium (more precisely, chemical equilibrium).However, because of the expansion of the universe, in general Γ becomes at some time smaller than the Hubble factor H. At this point the expansion of the universe thwarts the interaction, making it inefficient.Taking again the example of the annihilation, when Γ ≈ H it stops to be efficient and the species which was previously annihilating is said to freeze-out.Species that were originally in thermal equilibrium with the primordial plasma and then freeze-out are said to be thermally produced.
The reaction rate Γ for a given interaction is specified by the relevant cross section σ, the number density n of the particle species and the relative velocity v: We know from Statistical Mechanics that the number density at thermal equilibrium at temperature T is given by (see Ref. [67] for a detailed derivation of these formulae): where the subscripts refer to the fact that particles of a given species in thermal equilibrium are mostly relativistic or non-relativistic depending whether the temperature of the thermal bath is much larger or much smaller than the mass of the particle species under consideration.
The Hubble rate H is given by the Friedmann equation: Using Bose-Einstein statistics, one can show that: where g is again the number of spin states of the particle boson species.For fermions, the above result still holds true if multiplied by a factor 7/8.We can write in a compact way the radiation density as follows: where g eff is an effective total number of spin states of the relativistic particle species, taking into account the factor 7/8 for relativistic fermions.Note that g eff varies with time because massive species becomes non-relativistic as long as the thermal bath cools down, but, after DM and the heaviest particles have decoupled from the primordial plasma, g eff remains basically of order 10.
Therefore, we can finally combine our knowledge of Cosmology (i.e. the Friedmann equation) and Statistical Mechanics and obtain: where M P ≈ 10 19 GeV is the Planck mass.

A. Thermal relics
Relics are particles left over when their annihilation reactions are no more efficient, i.e. when their abundance attains a constant (comoving, i.e. not taking into account the standard dilution due to the expansion of the universe) value.If a species freezes-out when relativistic then it is dubbed hot thermal relic, whereas it is a cold thermal relic if it freezes-out when non-relativistic.

Hot thermal relics
A simple calculation of a hot thermal relic abundance can be carried out e.g. for neutrinos, which constitute the cosmic neutrino background (CνB) once decoupled.To this purpose, let us consider the interactions: of neutrino annihilation and neutrino pair production through fermion-antifermion annihilation.
Weak interaction rule this process, hence the cross section is given by: where G F ≈ 10 −5 GeV −2 is the Fermi coupling constant and T ν is the neutrino thermal bath temperature.The freeze-out condition Γ = H gives us, using Eq. ( 9): where we have taken into account that neutrinos are relativistic (i.e.v = 1, since we are working in natural c = = 1 units).For the same reason, from Eq. (10) we have that n ν ∼ T 3 ν and thus combining the above relations we get: For the above estimate we have used Fermi effective Lagrangian of the weak interaction and the fact that the particle is relativistic.The first approximation is reasonable because m ν ≪ m W , with m W mass of the W-boson, so that the interaction can be considered pointlike, and T ν ≫ m ν is a condition which holds essentially through all the thermal history of the universe (maybe up to recombination, if we consider m ν = 0.1 eV).
So far we have obtained the freeze-out temperature of the relic neutrinos background, but how do we now evaluate their abundance?Let us introduce the quantity Y ν = n ν /s, where s is the entropy density.Since S = sa 3 is conserved then Y ν ∝ n ν a 3 is a constant.Therefore: where T ν,0 is the neutrino temperature today.In particular, we can write: and if we assume neutrinos to be non-relativistic today, so that ρ ν,0 = n ν,0 m ν , we obtain: and from this relation one determines the so-called Cowsik-McClelland bound [25]: where h is such that H 0 = 100 h km/s/Mpc.

Cold thermal relics
Let us see now how to proceed in the case of a cold relic of mass m χ , for which we must use the non-relativistic expression (11) for the number density: Combining the latter with the freeze-out condition: and defining the new variable x ≡ m χ /T χ , we easily obtain: For the case of WIMPs we have m χ ≈ 10 2 GeV, and σ ≃ G 2 F m 2 χ , and the solution of the above equation becomes: From this, we can calculate the thermal relic density of WIMPs.As first step, let us write: It is important to stress that in the above computation we are assuming a DM-symmetric universe,1 , i.e. there is roughly the same amount of DM and anti-DM.Assuming again conservation of the entropy density s, so that aT = constant, we can use the relation we finally obtain: In the above estimate we have neglected the relative velocity of the cold relics, that could be estimated by considering the equipartition theorem: from which we get: Now let us consider the cross section of a WIMP: where E EW ≈ 200 GeV is the electroweak energy scale.We have then: which is of the correct order of the observed value Ω DM ≈ 0.26.This is the so-called WIMP miracle.However, note that in general this is not a special feature of WIMPs.Indeed, in order to guarantee the result in Eq. ( 34) one has just to use the following values: from which we can read the pretty weak constraint: Now it is natural to ask ourselves which is the range of masses expected for cold relics.Since we want to preserve the unitarity of the interaction (see e.g.[78] about the upper bound on a cross section coming from the requirement of unitarity of the scattering matrix), the cross section is bounded as follows: so that from which we finally get the unitarity limit:2 On the other hand, for WIMPs the cross section is given by σ ∼ G 2 F m 2 χ .Therefore, the smaller the WIMP mass is the earlier the freeze-out takes place and the higher the today WIMP abundance would be, eventually even overclosing the universe.The bound for which Ω χ = 1 is called the Lee-Weinberg bound [56].Plugging instead σ ∼ G 2 F m 2 χ into Eq.( 30), we obtain: showing that a WIMP particle with mass m χ = 10 − 100 GeV provides the correct present DM abundance.

B. Boltzmann equation
So far we have used qualitative arguments in to characterise a relic density, but of course precise calculations are based on the solution of the Boltzmann equation, which can be written in the following schematic way: where L is the Liouville operator and Ĉ the collision term (also an operator).The former is given by: in its non-relativistic and covariant forms (boldface letters represent usual 3-vectors and the dot a scalar product).The latter simplifies in a FLRW background because in this case the probability density function f depends only on the energy and on the time only due to the cosmological principle of homogeneity and isotropy: and the covariant Liouville operator can be written in the following form: The number density of a given species is given by the integral over the reduced phase space of f , i.e. the integral over the proper momentum volume: where g is, as usual, the number of spin (or helicities, for a massless species) states.Using this definition into the expression for the Liouville operator integrated in the momentum space we get for the left hand side of the (integrated) Boltzmann equation: Assuming a 2-to-2 interaction among species of the type the collision factor integrated in the momentum space of particle 1 can be written in the following form (see e.g.Ref. [53]): In the above formula, the superscripts "eq" refers to quantities (in this case the number densities) in thermal equilibrium, i.e. the expressions in Eqs. ( 10) and ( 11) that we have already encountered.
The average • • • is a thermal one, to be specified in a moment, σ is the total cross section of the interaction between species 1 and 2: i.e. summed over all the possible final states, and where we have defined: where p 1 • p 2 is the product of the two four-momenta of the particles 1 and 2, m 1 and m 2 are their masses and E 1 and E 2 their energies.This v Møl is the so-called Møller velocity and allows one to write the interaction rate in a covariant form, similar to the non-covariant expression, i.e.Γ = n 1 n 2 σv Møl .The thermal average is defined as: If the particle 1 and 2 are the same, i.e. the interaction under investigation is an annihilation process, Boltzmann equation can be finally written as follows: This equation must be solved numerically.Typically one expands σv in powers of the relative velocity v and the zero order term, i.e. assuming a constant σv , is usually called the s-wave approximation.
The above Boltzmann equation is the standard tool by means of which we are able to compute relic abundances.However, there exist important exceptions, such as resonances, thresholds and co-annihilation, for which that standard tool fail.We do not discuss them here, but refer the interested reader to Ref. [44].

C. Modified expansion history and relic abundance: the Kination example
Until now we have considered what happens to the thermal history of a certain particle species by focusing on the left hand side of the freeze-out condition Γ = nσv ≃ H.We have shown that e.g. for WIMPs the resulting relic abundance Ω χ is in good agreement with the observed value, i.e. ρ χ,0 ≈ ρ 0 .On the other hand, it is worth to also study the implications of modifying instead the expansion history of the universe, e.g. the evolution of the Hubble factor H. To this extent, let us consider a quintessence dark energy model accounting for the dynamics of H, e.g.
an homogeneous real scalar field φ, whose energy density and pressure are defined as follows: For simplicity, we assume a constant equation of state ω φ ≡ P φ /ρ φ , which implies, solving the continuity equation, that ρ φ ∝ a −3(1+ω φ) .When the field kinetic energy dominates over the potential one can see that ω φ ≈ 1 and ρ φ ∝ a −6 .This phase is called Kination [69].
Let us introduce a scale of temperature T KRE , which represents the temperature of kination and radiation equality, i.e. when the energy density of the scalar field and that of the radiation component are equal.We can express the Hubble factor above this threshold and the freeze-out condition as follows: since for T > T KRE the scalar field dominates and thus H ∝ a −3 ∝ T 3 , where the last proportionality comes from thermal equilibrium.Now, suppose that e.g.WIMPs decouple from the primordial plasma during the kination epoch.How would their relic abundance be affected, with respect to the standard case that we have seen earlier?First of all, using Eqs.( 28) and ( 29), we can establish also in the present case that: Then, we can estimate the ratio between the relic abundance of WIMPs (or any other cold species, since the ratio we are computing cancels out the cross section) in presence of a kination epoch and without it: where we have used the approximate result T standard freeze−out ≈ m χ /20 and the bound T KRE T BBN , since the kination epoch has to take place before the Big Bang Nucleosynthesis, otherwise it would spoil it.
The above example of kination teaches us that from a modified expansion history, different from the usual radiation-dominated epoch, we can get a huge enhancement of the cold thermal relic density.

D. Dark Matter after chemical decoupling
It is possible for DM in the early universe to be in kinetic equilibrium still after chemical decoupling.This happens due to elastic scattering processes between DM and e.g.Standard Model fermions: Here the relative velocity is unity because the temperature of the thermal bath is considered to be much higher than the fermion mass, which is therefore relativistic.In the case of a prototypical WIMP we assume that: Assuming m χ ≫ T , hence nonrelativistic WIMPs, then from equipartition one has on average that: and since m χ is the dominant energy scale, collisions entail a momentum transfer of the order of the thermal bath temperature, i.e. δp ≈ T .This implies that (δp)/p ≃ T /m χ ≪ 1. Upon a collision, a WIMP can gain or lose momentum, thus, as for the case of a random walk, after N collisions a WIMP particle momentum changes on average of the relative amount √ N (δp)/p.
When this quantity is of order one, the WIMP has changed momentum significantly and this corresponds to the following number of collisions: from which we can estimate a typical kinetic decoupling temperature for a WIMP particle, through the usual formula Γ ≈ H: from which we obtain for the kinetic decoupling temperature T kd : We can see how this results affects structure formation.In fact, for temperatures larger than T kd WIMPs cannot form halos because they are coupled to the primordial plasma and thus any fluctuation in their density is washed away.We can associate to T kd the following mass scale: where Earth's mass is: The above cutoff mass is the total mass of WIMPs contained inside the Hubble radius at the epoch in the history of the universe where the temperature of the thermal bath was precisely T kd .
So the first structures to collapse might be these minihalos with mass proportional to Earth's mass, which subsequently merge to form bigger and bigger halos, in a so-called bottom-up scenario of structure formation.However, note that the kinetic decoupling scale differs significantly for different theories of DM.
The picture changes dramatically in the case of hot relics.Taking neutrinos as an example, they decouple, by definition, when T ≫ m ν , and they can collapse and form structures only on large scales because on smaller ones they stream away (free-streaming), washing out any inhomogeneities.We can quantify this claim by computing the distance d ν travelled by the neutrino up to the moment in which T ≈ m ν , i.e. when it becomes non-relativistic, i.e. no more hot: Assuming again H ≈ T2 /M P , we then have d ν ≈ M P /m 2 ν and thus the cutoff mass for hot relics becomes: where one uses How can we compare these results with observation?Observational constraints from Lyman-α experiments give a cutoff mass scale of order: and we can conclude that DM can be thermally produced at most at 10 keV scale.In summary, structure formation differs strikingly between hot and cold DM, the former implying a Top-Down scenario, whereas the latter a Bottom-Up one.In Ref. [26] it has been showed that numerical simulations of structure formation in a universe with cold DM match the observed large scale structure of the universe better than the hot DM scenario, which is thus disfavoured.

III. DETECTABILITY OF PARTICLE DARK MATTER
We know that DM particles interact gravitationally with ordinary matter, but it is important to understand if they also couple through other known or unknown interactions.This is crucial in order to establish if there is a chance to directly detect DM particles via experiments in colliders or detectors in orbit.

A. Direct detection
Detecting particles which interact weakly interaction has always been known to be difficult.
For example, Bethe and Peierls in 1934 [14] estimated the cross section for the weak process: to be of order: and concluding that it is therefore absolutely impossible to observe process of this kind.Anyway, Bethe and Peierls position was perhaps too pessimistic, and since 1953 neutrinos have been abundantly detected in inelastic process such as the one in Eq. ( 69), whereas elastic neutrino scattering was observed in 1973 in the bubble chamber Gargamelle [47].
Let us use WIMPs again as prototypical DM particles and look for the relevant range of energies and masses that we could expect on the basis of direct detection.A mass m χ WIMP with velocity v colliding and bouncing off a much heavier nucleus of mass m N transmits a maximal recoil momentum p = 2m χ v. Hence the nucleus acquires an energy: We can bound the maximal velocity that a DM particle could have with the escape velocity of the galaxy in which it is located, which is usually of order v max ≈ 500 km/s.This implies an energy of order of the keV for DM particles with mass of the order of the GeV.
We also have bounds coming from the event rate R = Kφσ, where is the number of targets per kg of material with atomic number A and the incident flux φ is: from which we get: R = 0.06 events kg day In order to achieve direct detection experimentalists have to deal with two crucial points: enough signal events and enough background suppression, which is not an easy task considering the weakness of the DM signal.In particular, sources of noise are: • Slowly decaying "primeval" nuclides such as Uranium, Thorium and Potassium-40, whose abundance is about 10 −4 and half-life is ≈ 10 9 years; • Rare, fast decaying trace elements like tritium and Carbon-14, whose abundance is about 10 −8 and half-lives of the order of 10 years.
In order to achieve shielding from these environmental noises experiments are often conducted underground, as in the Italian Gran Sasso National Laboratory, where 1400 tons of rock protect the experimental setup from cosmic rays.
There is also the possibility of detecting a DM signal in the radioactive background of our planet, for example through seasonal and diurnal modulations, see for example Refs.[24,33] or by directional information, see e.g.Ref. [73].
Let us consider the direct detection event rate: where N T is the number of target nuclei, n χ the DM particle density, v the DM particle velocity over which the average is performed and E R is the recoil energy, which has the following expression: where q is the transferred momentum, m T is the target mass, and µ is the reduced mass.Taking the differential, one has and thus we can write: where f (v) is the velocity probability distribution function of the DM halo.
Of course, we need to calculate the scattering cross section in order to evaluate the above quantity.In the nonrelativistic limit, the scattering matrix element is given by the Fourier transform of the WIMP-nucleus potential: where f | and |i are the initial and final quantum vector states.To lowest order in v, the potential can be approximated to be composed by axial and spin-independent contact interactions: where the G s,a represent, respectively, the effective scalar and axial couplings between DM and nucleons, and s χ and s n are the DM particle and nucleon spins, respectively.
Assuming that a given microscopic theory of DM (i.e.some DM Lagrangian) is given, how we can infer something about the DM-nucleus cross section?The proper way to deal with this kind of questions is via Effective Field Theory (EFT), i.e. assuming the mass of the mediator particle of the interaction between DM and Standard Model particles to be larger than the typical momentum transfer.As an example, in the case of milli-electric charged DM the cross section is of the form: where α is the fine-structure constant, Z is the nucleus atomic number and ǫ ≪ 1 is the DM electric charge.ones.In both cases the rates are: where ρ DM is of course the DM density for a DM particle of mass m χ (ρ DM /m χ is thus the number density, for a cold DM species), σ is the annihilation cross section and v the relative velocity of the annihilating DM particles, τ dec is the lifetime of the unstable, decaying DM particle and N SM are the numbers of Standard Model species (hence the subscript SM) into which DM particles can annihilate or decay.
We can infer informations about the product σv by DM thermal production history, but how can we obtain informations about the decay rate?If we suppose that such decay is mediated by some mass-scale M particle, then a combination of the type: would give: which is a too short-lived candidate.If we try another combination of powers, i.e.
then we would get which is more interesting, because the time factor 10 27 s is much larger than the age of the universe (10 17 s).As mass scale M, the Grand Unification Scale 10 16 GeV has been chosen above.
If we ask ourselves what we can say about the annihilation final state, our answer is highly model-dependent and different internal symmetries give very different results.For example: • If DM particles belong to an SU(2) multiplet, then we expect well-defined combinations of ZZ, W W final states; • In Universal Extra Dimension (UED) theories [49], DM is the Kaluza-Klein KK-1 mode • We expect a special selection rule, e.g.helicity suppression for Majorana fermion (analogous to charged pion decay): The decay/annihilation rate of DM can be constrained indirectly in several ways, a possible classification being the following: • Very indirect.Looking for DM effects induced in astrophysical objects or in cosmological observations.
• Pretty indirect.Using probes that do not trace back to the annihilation event, since their trajectories are bent as the particles propagate.For example, cosmic rays.
• Not-so-indirect.Using neutrinos and gamma rays, which have the great added advantage of traveling in straight lines.

Very indirect probes
Examples of very indirect probes include e.g.
• Solar Physics.It is possible that DM could affect Sun's core temperature, or the sound speed in its interior; • Neutron Star Capture.DM can lead to Neutron star capture that eventually leads to the formation of black holes (notably e.g. in the context of asymmetric DM); • Supernova and Stars, in which DM could be responsible for cooling processes; • Protostars, e.g.WIMP-fueled population-III stars; • Planets warming; • Cosmological observation, where the DM content has strong implication for e.g.Big Bang Nucleosynthesis, or the Cosmic Microwave Background spectrum, also affecting the time of the recombination, or the structure formation process.

Pretty indirect probes
A good idea is to use rare cosmic rays events where antiprotons and positrons are relatively abundant due to inelastic scattering off interstellar medium protons.
An interesting probe is given by antideuterons, or even anti-3 He, for which the relevant process is: and which has a large energy threshold, about 17 GeV, so typically large momentum, even if caused by very low momentum DM particles.
Positrons, and in part antiprotons, have attracted attention due to "anomalies" reported by the detector PAMELA (part of the AMS-02 project mounted on the International Space Station) [5].
The general scheme for treating galactic CR is with diffusion (the so-called leaky-box) models; defining the differential in the energy number density of cosmic rays as dn dE = ψ (x, E, t) we can write a master equation of the form: where D(E) represents the diffusion coefficient, b(E) describes energy losses and Q include all the possible sources.Note that the above master equation can be made arbitrarily more sophisticated, for example adding the effect of cosmic rays convection, or diffusive re-acceleration, or considering fragmentations and decay.However, we do not consider these possibilities here.
Let us assume assume as boundary conditions a cylindrical slab for our galaxy with typical dimensions of R ≈ 10 kpc for the radius and h ≈ 1 kpc for the height.Let us also assume an energy dependence for the diffusion coefficient of the following form: with E 0 ≈ 1 GeV, D 0 ≈ 10 28 cm 2 /s and δ ≈ 0.7.We can furthermore simplify the diffusion equation assuming a steady-state regime, i.e. assuming ψ to be constant.In this case the typical diffusion and energy loss time-scales look like: and the diffusion equation simplifies as follows: and has the simple solution: We should distinguish distinguish among different energy phases.For cosmic rays accelerated via a Fermi mechanism, i.e. through multiple reflections in a magnetic field and with no energy losses (so in this case we neglect τ loss ), then one has: which is a good estimate when considering cosmic rays protons, for which indeed energy losses are irrelevant.In the case of cosmic rays electrons instead energy losses are efficient only above a certain energy threshold, so that we can write: where u ph is the background photon energy density, B is the ambient magnetic field (measured in micro-Gauss µG), b 0 IC ≈ 0.76 and b 0 sync ≈ 0.025 × 10 −16 GeV/s represent the energy loss coefficients for inverse Compton scattering and synchrotron radiation emission, respectively.Therefore, we expect a broken power-law for ψ(E), since: Instead, the secondary-to-primary (secondary cosmic rays are those produced from the original particles which are accelerated to high energies and constitute the primary ones) ratios are generically: The broken power law in Eqs. ( 97) and ( 98) is in very good agreement with observation, but the prediction for the secondary-to-primary ratio is at odds with the observed rising positron fraction. Indeed: • There is no excess of antiprotons, so DM should be leptophilic, which is possible but not generic; • There is no observed secondary radiation due to bremsstrahlung or inverse Compton scattering; • A very large pair annihilation rate is required for thermal production, which leads to unseen gamma-ray or radio emission, i.e.
σv ≈ 10 −24 cm 3 s • m χ 100 GeV (100) A possible way out is to consider the presence of a nearby point sources, such as a pulsar, that injects a burst of positrons and that entails the following correction: Then, we can estimate the age and the distance from us of the putative pulsar as follows (using E = 100 GeV): An interesting possibility is to use anisotropy in their relative distribution in order to disentangle pulsars from DM contributions, but complications appear since the Larmor radius of a typical pulsar magnetic field is of the order of the solar system size, so a very high spatial resolution is needed.

Not-so-indirect probes
Even if very hard, it is possible to detect neutrinos from astrophysical sources.The idea is very interesting, since DM could be gravitationally captured in celestial bodies, here experiencing accretion and starting annihilation.If the latter process is in equilibrium, we expect that a large flux of neutrinos could escape, since they have a large mean free path in matter.
The best target is the Sun, since it is sufficiently large, very close to us and naturally emits low-energy neutrinos.Let us estimate the above-mentioned process quantitatively.First of all, Sun's capture rate C ⊙ can be written down as follows: where the capture cross-section (DM particle χ captured by e.g. a proton p) can be estimated as: and hence: The number of accreted DM particles N is: where being A ⊙ the annihilation rate and E ⊙ the evaporation rate, taking into account the effect of evaporation of DM particles due to their thermal velocity.The annihilation can be written as follows: where t ⊙ ≈ 4.5 Gyr, or 10 17 s, is the age of the Sun.Defining the equilibration time as: then thermal DM is in equilibration as long as WIMP-nucleon cross section is larger than 10 −41 cm 2 .With equilibration, the flux of neutrinos of flavour f only depends on the capture rate and is: where D ⊙ is the diffusion coefficient of the Sun and the inj subscript on the right hand side refers to the injected flux, i.e. the intrinsic flux which the Sun naturally produces.
The number of events of µ-neutrinos detection at the IceCube4 observatory is expected to be: The most probable final states upon DM annihilation are W W , ZZ or leptons-antileptons pairs.
So far no anomalous events from the Sun have been observed, and detection from Earth is even less promising.Instead, an opportunity comes from low-energy-threshold subdetectors, such as DeepCore and PINGU (both belonging to IceCube).
Another interesting idea is to look for photons coming from a possible coupling between DM and Standard Model particles.For example, a "smoking-gun" signal would be a monochromatic, spatially extended gamma-ray emission, because a monochromatic gamma-ray line with a diffuse morphology has no astrophysical counterparts.From data by the FermiLAT 5 Bringmann [18] and Weniger [79] identified a 130 GeV line from the galactic center that was initially explained by recurring to DM annihilation.Unfortunately, the signal turned out to be a statistical fluke, since it was too narrow and its significance did not increase with time.Finally, when the data were processed with the Pass 8 software no excess was detected.
Generally, we classify photon signals into two categories; in the primary one we have prompt emission of photons and bremsstrahlung radiation;6 in the secondary one, photons undergo inverse Compton scattering or constitute synchrotron radiation.
The average photon energies of the second kind of processes (inverse Compton scattering or synchrotron radiation) are of the order: where E ′ 0 is the average final energy of the scattered photon, γ e is the electron Lorentz factor and E 0 is the photon thermal energy which can be, for example: for the CMB thermal bath or for different instances of interstellar medium.The typical energy scales E e for electrons produced in DM decay or annihilation and the energy of scattered CMB photons can be estimated as follows: while the typical frequency of the synchrotron radiation is: The prompt emission of photons instead depends only on the annihilation final state and the target of choice; the flux of photons φ γ produced by DM from a given direction ψ within a solid angle ∆Ω is: where dl is an integration along the line of sight and f is the sum over all the possible final states.It is useful to define the socalled J-factor as the expression above between curly brackets: in order to characterise typical sources of radiation for an angular region that varies from 1 to 0.1 degrees.Some of them are the following: a. Dwarf Spheroidal Galaxies • Draco, J ≈ 10 19 GeV 2 /cm 5 ± a factor 1.5; • Ursa Minor, J ≈ 10 19 GeV 2 /cm 5 ± a factor 1.5; • Segue, J ≈ 10 29 GeV 2 /cm 5 ± a factor 3.
In order to have detection we need to collect enough photons (i.e.we need a high signal-to-noise ratio).The number of detected photons is given by: where A eff is the effective area of the detector and t obs is the typical observation times, and in the second approximation we considering typical values from the Fermi-LAT experiment.The energy range is of order: so that the photon flux can be written as: from which we obtain a J factor of order: that is a bit bigger than the typical values of J given above for dwarf spheroidal galaxies.Finally, we can compute the following rough constraint on DM pair annihilation rate depending on the DM mass: For a more detailed analysis see for example Ref. [4].In addition, we have also production of monochromatic photons from χ + χ → γ + γ, for which:

The Galactic center excess
In Ref. [40] the authors realised, after analysing data from the the Fermi gamma-ray space telescope, that there exists an excess of gamma-rays over the background in the center of our Galaxy.
Such excess was also subsequently reported independently and with a variety of different assumptions for the background in Refs.[2,27,42].If we ask ourselves what might produce such an excess, it is very easy to fit the data with a signal from DM annihilation; indeed both the morphology and the spectrum are compatible with some annihilation processe, and also compatible with constraints from CMB and observations of dwarf spheroidal galaxies.On the other hand, trying to explain such an excess with emissions of astrophysical origin, e.g.unresolved pulsars, appears problematic for a variety of reasons.
However, it is possible that the starting assumptions of these analyses are in fact wrong, e.g.
there is no excess at all and our models for the diffusion emission are simply inadequate for current data.Indeed, diffusion emission is a complicated physical process that involves numerous ingredients; cosmic rays are transported from their acceleration sources through the interstellar medium; most gamma-rays are produced by p-p inelastic collisions via the production of neutral pions that decay and release radiation, or by bremsstrahlung from charged cosmic rays that decelerate due to interstellar hydrogen.
Why do we expect that our model of emission might be not realistic?The main reason is that the following assumptions are usually used to evaluate the quantity of radiation that we expect from the galactic center: i) a bi-dimensional gas density distribution; ii) a bi-dimensional model of cosmic rays propagation; iii) the hypothesis of steady-state and iv) a simplistic source distribution for the cosmic rays.It turns out that each of these assumptions costs a systematic effect which is of the same order as the excess itself.How can we improve the overall picture towards the next generation of diffuse gamma-ray models?We should consider a 3-dimensional model for the gas density distribution and for the propagation of cosmic rays, and the effect of transient processes, such as the presence of cosmic ray bursts, together with a physically motivated choice for the cosmic rays source distributions.
In Ref. [21] the authors show how using physically-motivated models for the location of the acceleration sites of cosmic rays leads to a significant reduction of the Galactic center excess, as well as to drastic changes to its spectrum and morphology.Additionally, these novel diffuse emission models provide globally a better fit to the Galactic emission as a whole.
While significant progress is made in the direction of better and more physically motivated models for the diffuse emission, the discrimination between unresolved point sources and the diffuse emission is also part of the issue of extracting an excess and understanding its properties.Beyond and because of these considerations, we remain skeptic about firmly establishing a conclusive DM detection signal from the Galactic Center; while we still are optimistic about the possibility of detecting DM from gamma-rays.

Collider production of Dark Matter particles
In principle DM particles can be produced in collider experiments, but the expected production rate is very low if compared to Galactic DM fluxes.A possible idea is to look for anomalous events with missing energy and Standard Model particles, such as for example monojets or monophotons.
We can use two different approach: 1. Top-down: pick a model and scan the parameter space (e.g supersymmetry or unified theories).
2. Bottom-up: use some effective field theory (EFT) or simplified models to sketch of how DM could manifest itself at colliders.
In the EFT approach one comes up with a list of effective operators compatible with the constraints on DM coming from collider experiments [41] and then apply the following algorithm: compute the production cross section and simulate the events, and finally set limits on the effective operator scale.
The main issue of such approach is that EFT is meaningful only in a certain range of energy, whose scale say Λ is given by: where M is the interaction particle mediator mass, and the g i 's are the couplings to Standard Model particles and particles beyond the Standard Model.This entails the question of whether the constraints that we assumed make sense, i.e. if the typical energy of the reaction (say the momentum transfer q) is smaller then 4πΛ.
An alternative approach is to use simplified models for the interactions of DM with Standard Model particles; for the coupling with a scalar S and a vector V µ mediator, we have for example the following Lagrangians: where the y and the g's are coupling constants, the γ µ are Dirac's matrices, the symbols ⊃ means that the contributions on the right hand side are only a part of the total Lagrangian (there still is the Standard Model one) and h.c.means "hermitian conjugate".
The idea is to set meaningful constraints on combinations of mediator mass and couplings constants for given DM masses.Then, we can compare with direct detection results.
An additional probe of DM within colliders is given by the invisible decay of the Higgs boson into DM, with an interaction term of the form: where H is the Higgs scalar (not to be confused with the Hubble factor) and λ Hχ χ is the DM-Higgs coupling.
The invisible Higgs decay modes may stem from decays of the Higgs boson directly into DM particles, or into new particles that in turn decay into DM particle, see for example Ref. [1].

Axions and Axions searches
Axions are hypothetical particles postulated in order to solve the strong CP problem in Quantum Chromodynamics (QCD) [74].Several production mechanisms can lead to abundant production of axions in the early universe.If the universe is indeed filled with such primordial particles they can be, in principle, good DM candidates.The QCD Lagrangian is given by: where G a µν is the gauge-field (the gluon) tensor, the n q's are the quark fields, D µ is the gaugecovariant derivative.The last term, where θ is a coupling and ǫ µνρσ is the Levi-Civita symbol, is innocuous perturbatively, being a total derivative, but it enters phenomenologically via nonperturbative QCD effects, producing a large neutron-electron dipole moment: where the last bound is the experimental result [66].The Peccei-Quinn solution [65] consists of promoting θ to a dynamical variable which is driven to vanish by its classical potential.Then one postulates a global (quasi-)symmetry of the theory U(1) P Q (broken by perturbative effects); such a symmetry is spontaneously broken at a scale f a .In this way, the axion turns out to be the (pseudo) Nambu-Goldstone boson associated with the group U(1) P Q .Its mass is given by: where Λ QCD ≈ 200 MeV is the QCD energy scale.As shown in Ref. [77], the ground state of the axion potential solves the strong CP problem.These QCD effects produce effective couplings (slightly model-dependent) to fermions and photons: where a represents here the axion field, N is the number of fermion species considered, γ 5 = iγ 0 γ 1 γ 2 γ 3 , E and B are the electric and magnetic fields and, as usual, the g's are the coupling constants.
The setup is similar if we consider Axion Like Particles (ALP): we still have a new global U(1) symmetry which is spontaneously broken by a hidden Higgs-like mechanism at scale v h .We parametrize the Higgs field as : The potential for the ALP field a(x) is flat, and depending on the model realizations one generates couplings to SM particles: where F µν is the electromagnetic field tensor.Because of the coupling to Standard Model particles, in particular photons, axions decay into two photons: To have a sufficiently long-lived axion we must demand also: The axions can have a dramatic impact on stars due to their Compton-like and Bremsstrahlung-like interactions: which produce an axion luminosity, for example for the Sun, of order: Given the observed luminosity of the Sun, axion emission would require an enhanced nuclear energy production, which in turns would imply a larger neutrino flux.
Axions can also have impact on the physics of supernovae.If the axion mass is small enough, axions would free-stream out of supernovae, causing their cooling; the corresponding axion luminosity is of order: For consistency we should require that this luminosity does not exceed the value of the thermal neutrino luminosity L ν ≈ 10 53 erg/s , which would happen for: On the other hand, if axions are too massive, they get trapped and not contribute to supernovae luminosity efficiently.It turns out that axions would affect supernovae burst duration in the mass range: How can axions be produced?If they are produced thermally, this should happen via interaction among gluons and quarks of the form: a + g ↔ q + q or g + g , or a + q (q) ↔ g + q (q) .
This kind of process comes from the QCD Lagrangian term proportional to θ, from which we can also roughly estimate the following cross section: where α s is the coupling constant relative to the strong interaction.From the above formula we can obtain a rough estimate of the temperature at which the processes considered above go out of equilibrium and the axions freeze-out: T a,freeze−out ≈ 10 11 GeV f a 10 12 GeV .
However, at lower temperatures (below the QCD phase transition scale) the interactions with pions must be considered: Such process is in equilibrium when: Finally, we obtain the density of axions as a function of their mass: but we know already that this does not work since hot DM is not good for structure formation.
Also, other constraints on axion mass are also not good, because in order to have Ω a ≈ Ω DM ≈ 0.26 we would need m a ≈ 10 eV, and the resulting axion lifetime would be too short.
But then, how about non-thermal production?We can consider for example the so-called misalignment mechanism for producing cold axions, which results in the following cosmological density, see Ref. [72]: where θ 1 is a mixing angle.If we assume that the value for θ 1 is given by its root mean square we get: Ω a,mis, RMS h 2 ≈ 0.13 m a 10 µeV −1.18 . (153) On the other hand, if the distribution for θ 1 has a non trivial topology, a network of axion strings and domain walls can be generated; the resulting cosmological density is in general modeldependent and the results of different groups are in general in disagreement.For example, Refs.

Sterile neutrinos and the 3.5 keV line puzzle
One of the few still open problems of the Standard Model of particle physics is the fact that neutrinos are not massless.A simple way out is to consider a set of n gauge singlet fermions N a (the sterile neutrino field) so that: where / ∂ ≡ γ µ ∂ µ and Lα represent the ordinary neutrino field.The mass matrix can be computed to be the following: where the Higgs field vacuum expectation value H is usually written as iv/ √ 2 (do not confuse v with any velocity).Just for simplicity, assume only one sterile neutrino family and only one ordinary neutrino family, so that the above matrix become a 2 × 2 matrix: It is easy then to compute its eigenvalues, which are: hence, if the condition yv ≪ M is satisfied, then λ + ≫ λ − .In general, since λ + λ − = y 2 v 2 /2, if one eigenvalue grows the other decreases and viceversa.This is the so-called Seesaw mechanism.
So if yv ≪ M is satisfied, we have one large eigenvalue M corresponding to the sterile neutrino mass and one small y 2 v 2 /M eigenvalue corresponding to the ordinary neutrino mass.
Sterile neutrinos mix via explicit (but possibly very small) mixing with ordinary neutrinos; as such, they decay (into the 3 Standard Model neutrinos): where m νs is the sterile neutrino mass, θ is the mixing angle and where we have used θ −2 mν s keV −5 ≫ 1 because sterile neutrinos are fermions and thus m νs > keV [75].
Sterile neutrino are produced essentially via freeze-in [45].Briefly, the freeze-in mechanism relates to particles which are not in thermal equilibrium with the primordial plasma but still interact (and are produced) with some of the species which instead are in thermal equilibrium with the primordial plasma.This interaction stops being effective when the latter species freeze-out, so that the weakly interacting species are said to freeze-in.For sterile neutrinos we can estimate the production rate as follows: where the mixing angle θ is usually considered to be constant, or in some refined models has a temperature dependence: Abundance of sterile neutrinos can be estimated as the following value: [32]: An additional important effects that can be considered are the Mikheyev-Smirnov-Wolfenstein effect [60,81], and the Shi-Fuller resonant production effect [70].
What if sterile neutrinos are not produced thermally?A simple alternative is to consider nonthermal production from singlet scalar coupling with the field S: and (or) SH † H, S 2 H † H, that gives: Sterile neutrinos are also interesting from the standpoint of structure formation, since: Another interesting feature of sterile neutrinos is that they can provide an answer for the so called pulsar kick problem, see for example Ref. [15].
What about detection of sterile neutrino DM?The key mechanism is loop-mediated sterile neutrino decay via weak interactions and mass-mixing with active neutrinos, leading to a naive estimate for the decay width and the photon flux from decay of: elsewhere, or we can use some different probe rather than the spectrum.About the former option, there is no trace of signal from dSph galaxies [58], nor from stacked galaxies and groups with low plasma temperature [9], nor from Andromeda [50].Also, there is no signal from dedicated observation of Draco dSph [51].What about using something different than spectrum?We can use morphology, e.g.look at where the 3.5 keV radiation come from.However, morphology of Perseus and the Milky way shows that decaying DM is strongly disfavoured [20].
In Any competent theoretician can fit any given theory to any given sets of facts.
The 3.5 keV line is an exciting excuse for a new mechanism generating a DM signal.Consider for example an interaction of the form [28]: which models interactions between DM and interstellar plasma, resulting in a signal which is proportional to the product of the DM and gas densities, and it is also a good thermal relic.Why should we be excited by this model?
• It is a new indirect detection channel; • it has an unmistakable signature, free of background; • Is a good model, in the sense that it is economic, with a natural UV completion and a thermal relic DM.
Moreover, it is a highly falsifiable scenario for the following reasons: • Its line shape, e.g.geometric average of thermal DM velocities, can be resolved by e.g. the Hitomi/astro-h satellite; 7 William of Occam, c. 1286-1347. 8Quoted in Ref. [57], Sec.2.5.1,The psychology of astronomers and astrophysicists.
• It has unique morphology • It has unique target-dependence • Lines could appear anywhere in the spectrum, from eV, to UV, to X-ray.
Possible future tests include implications on structure formation and small scale structure.

IV. DARK MATTER BESTIARIUM
In this section we introduce a list of notable candidate models of DM particle.

A. Gravitinos
Gravitinos are prototypical DM with Planck-suppressed interactions [64].The prototype Lagrangian is: and the resulting cross sections are: where we used the symbol G to indicate a gravitino particle.Finally, the freeze-out temperature is: Thermal gravitinos are hot relics (in fact they were the first super-symmetric DM candidate ever proposed) and their abundance is calculated as: However, this calculation neglects single-gravitino processes that maintain gravitinos out of equilibrium, for example: where V is a gauge boson and λ its gauginos (its supersymmetric partner).Generically, there is a gravitino overproduction problem, so we need to dilute them away: where T RH is the reheating temperature.
Gravitinos are also produced from next-to-lightest supersymmetric particles (NLSP) decays [11]: Light gravitinos are already ruled out, whereas heavier gravitinos are still viable for certain combinations of mass and T RH .Moreover, they might also be NLSP with long lifetimes: An important question is of course: can we ever hope to detect gravitinos?There are basically two strategies; the first one is to use long-lived s-leptons produced at LHC trapping them into water tanks, and then wait their decay to gravitinos.The second one is to use neutrino telescopes to search for neutrino-induced long lived s-leptons:

B. WIMPzillas and super-heavy Dark Matter candidates
This interesting class of models grounds on the presence of non-thermally produced super heavy particles, the so-called WIMPzillas [52], with the following properties: • The DM particle is never in thermal equilibrium; • The particle mass is comparable to the inflaton mass say M φ ; • The particle lifetime is much longer than the age of the universe.
The following generic expression gives us the density of such heavy particles produced by gravitational interaction at the end of inflation [22,54]: where Ω R is the radiation density and and T RH is the reheating temperature.Such kind of particles could be detected in ultra high energy cosmic rays detectors.
Another interesting super heavy DM candidates are the strangelets [80], i.e. macroscopic clumps of quark matter, with radius between 10 mm to 10 cm and masses around 10 9 − 10 18 g.They are in principle detectable via cosmic rays collisions.

C. Self-interacting Dark Matter
In order to have at least one interaction for particle during the age of the universe this class of models must satisfy: which in turn implies a cross section of order: Events of this type could be object of study in DM halos, which can be considered as nonrelativistic colliders.For example, dwarfs galaxy are equivalent to the B-Factory collider, Milky way-like galaxies reproduce events like in the LEP, while big clusters of galaxies show events that on Earth are reproducible at the LHC.
Which model of particle physics model could work for self-interacting DM?We have a quite large landscape; consider for example a kind of Dark QCD with glueballs DM [16]: Another possibility is given by dark atoms, whose relic density computation is still a work in progress.Another one: light mediator [76] such that: Effects should turn off at large velocities but are significant at small scales or low velocities.We should assume that the mediator has a small coupling ǫ to Standard Model fermions.
Another possibility is given by a model with a complex scalar field and potential [59]: which gives us the following cross section: from which, using m s = 1 MeV and λ ≃ 0.5 we get

D. Asymmetric Dark Matter
The idea of asymmetric DM, formerly known as Technocosmology [63], is that there exists a new dark baryon quantum number B D in addition to the standard baryon number B V , with DM being the lightest charged particle under B D .We can define the following linear combination: which are conserved or broken at some energy scale.If the universe have a generalised global baryon symmetry the first combination is unbroken at all energy, while the second one is broken.
At the breaking scale we have a net production of ∆B broken such that If the universe does not have a generalised global baryon symmetry, instead some asymmetry is generated either on B D or B V .The mass of such asymmetric DM candidate is [46]: where the ratios η in the right hand side of the above equation are the charge-to-entropy ratios.For a baryon-symmetric universe we have: if the chemical decoupling of the two sectors happens when DM is non-relativistic.

E. Minimality
With minimality we mean a theory which contains a minimal number of ingredients, e.g.heavy neutrinos.SU(2) neutrinos are excluded by direct detection, unless we introduce a suppressed coupling to the Z boson.Another possibility is to consider a minimal number of new fields, e.g. a real scalar singlet.In this case [38] the Lagrangian becomes: and: The "original" minimal DM model (MDM) is given by an SU(2) L multiplet with a given spin and mass, see Ref. [23]: We have to require vanishing hypercharge and electric charge equal to: This implies that no operators could mediate decay (i.e.terms of type χeH, → χ = eH) and we are left with a very constrained set of possibilities: n = 5 M = 9.4T eV spin 1/2 fermion.
Such heavy SU(2) particles have important Sommerfeld enhancement effects, affecting indirect DM constraints.
Other minimal possibility is to consider an inert doublet model; e.g. a new SU(2) Higgs doublet with Z 2 symmetry [29].It turns out that lightest Z 2 -odd particle is stable and a good WIMP.

F. Dark Photons
Another interesting possibility is given by dark photons, see for example Ref. [34].Here we have gauge bosons associated with a new U(1) ′ symmetry, with a kinetic mixing term: m A ′ ≪ 2m e , can be long-lived enough (207) These are produced in oscillation or processes like: γ + e ± → A ′ + e ± e + + e − → A ′ + γ , or from dark photon condensates, (similar to the axion misalignment mechanism): If new dark photons are charged under a non-abelian symmetry group, kinetic mixing term is automatically prohibited, and the particle is stable; the Lagrangian is of type: and has interesting phenomenology: e.g. for SU(3) we have semi-annihilation and bright gamma ray lines (∼ α instead of ∼ α 2 , ψ i ψ j → ψ k γ).

V. CONCLUSIONS
The book of physics contains several blank pages when it comes to the fundamental nature of DM as a particle.The discovery of the particle nature of DM will likely offer a portal on new physics beyond the Standard Model, and opportunities to understand how the universe on large scales is shaped by the microscopic properties of DM.It is important for the scholar who embarks in the study of particle DM to be keenly aware of the scope and nature of work that others have produced thus far; especially useful it is to be familiar with the "bag of tricks" that scholars have accumulated so far and to be able to master the key tools of this trade.These notes are meant to help with this critical and exciting task.

C
. Observational constraints on Dark Matter interactions While only loose model-independent constraints exist on DM self-interaction and on DM interaction with Standard Model particles, a broad portfolio of experimental probes has been deployed over the years with the aim of searching for signals coming from specific DM particle models.The best-developed searches focus on WIMPs, and seek to capture the elusive signals from elastic scattering of WIMPs off nucleons (direct detection), signals of WIMP pair-annihilation or decay (indirect detection), and the pair-production of WIMPs at colliders in conjunction with detectable Standard Model particles.So far, while several signals have been claimed to be possibly associated with different DM candidates, including WIMPs, no conclusive detection has been agreed upon.Null results from of hypercharged gauge boson, thus the scattering matrix element M is proportional to the fermion hypercharge Y f , i.e. |M| 2 ∝ |Y f | 4 , and the preferred annihilation modes are the up quarks Y u L = 4/3 and charged leptons Y e R = 2.
Another interesting possibility of detection is the indirect one, namely to use debris (i.e.Standard Model particles) of DM pair annihilation, which are likely to be relevant if DM is a cold thermal relic, or to use decay of DM particles into Standard Model (hence in principle detectable) summary, we are left with the open question: Dark Matter or Potassium?We can look for a possible answer in the old literature: Entia non sunt multiplicanda praeter necessitatem, 7 .Occam's razor favours Potassium against sterile neutrinos because of the lack of signal from Draco and because of the morphology of the thermal lines.But what if it is another type of DM? Any challenge is also an opportunity and an interesting riddle for theorists; according to Redman's