A new mechanism of sterile neutrino dark matter production
Johannes Herms111johannes.firstname.lastname@example.org , Alejandro Ibarra222ibarra@tum.de , Takashi Toma333takashi.email@example.com
Physik-Department T30d, Technische Universität München,
James-Franck-Straße, D-85748 Garching, Germany
School of Physics, Korea Institute for Advanced Study, Seoul 02455, South Korea
One of the simplest extensions of the Standard Model consists in adding to the matter content a spin 1/2 Majorana fermion, singlet under the Standard Model gauge group, usually denominated sterile neutrino, right-handed neutrino or simply fermion singlet. The symmetries of the model allow a Majorana mass term for the sterile neutrino, with size which is a priori unrelated to the electroweak symmetry breaking scale and which can range between 0 and the Planck scale, as well as a Yukawa interaction of the sterile neutrino with the left-handed lepton doublets and the Higgs doublet. This model, in particular, leads to non-vanishing active neutrino masses after electroweak symmetry breaking through the renown seesaw mechanism [1, 2, 3, 4, 5].
The sterile neutrino has been advocated as a dark matter candidate. The Yukawa interaction leads to the production of sterile neutrinos in the very early Universe from interactions with the plasma of Standard Model particles . Nevertheless, current upper bounds on the strength of the Yukawa coupling from X-ray observations rule out this framework as the only mechanism of sterile neutrino dark matter production [7, 8]. On the other hand, the production rate could be enhanced in the presence of a lepton asymmetry , or by further extensions of the model e.g. by the decay of a scalar singlet which couples to the sterile neutrinos [10, 11, 12, 13]. Due to the smaller Yukawa couplings necessary to reproduce the correct dark matter abundance in these frameworks, the limits from X-ray observations are accordingly relaxed. However, limits from structure formation from the Lyman- forest and from Milky-Way satellite counts severely constrain the remaining parameter space of the model .
In this paper we present a mechanism of sterile neutrino production which is realized even when all portal interactions with the visible sector are negligibly small. Our framework requires the sterile neutrino to have sizable self-interactions, such that their initial population thermalized in a dark sector plasma at some very early stage of the history of the Universe. We show that the self-interactions induce the freeze-out of the sterile neutrino abundance through 4-to-2 annihilations, and we argue that, for appropriate parameters, the relic population of sterile neutrinos can account for all the dark matter of the Universe. Furthermore, the self-interactions could potentially alleviate the small scale problems of the cold dark matter paradigm (see e.g. ). Opening the portal interactions to the Standard Model, which in this framework are no longer directly constrained by the requirement of generating the observed dark matter abundance, allows for additional tests of the model from the decay of the sterile neutrino dark matter.
This framework is constructed along the lines of Strongly Interacting Massive Particles (SIMPs) as dark matter candidates. General considerations of the freeze-out of SIMPs via 3-to-2 and 4-to-2 annihilations were discussed in [16, 17]. Various models of SIMP dark matter have been proposed, where the dark matter candidate is a pseudoscalar [18, 19], a scalar [20, 21, 22, 23, 24], a spin 1/2 Dirac fermion , or a spin 2 boson . Here, we consider for the first time spin 1/2 Majorana fermions as SIMP dark matter. In this case, Lorentz invariance forbids the 3-to-2 annihilation, such that the freeze-out occurs via 4-to-2 annihilations. Furthermore, the Pauli exclusion principle requires the annihilation cross-section to be -wave suppressed, regardless of the nature of the self-interaction.
The paper is organized as follows. In Section 2 we introduce the framework where sterile neutrinos are SIMP dark matter candidates, with interaction described by a dimension six operator, and we discuss its phenomenology. In Section 3 we extend the model to include a scalar field which is responsible for the self interactions, and we discuss in detail the role of the scalar field in the cosmology, concretely in populating the dark sector plasma with sterile neutrinos. Finally, we present our conclusions in Section 4.
2 Sterile neutrinos as SIMP dark matter candidates
The sterile neutrino, that we denote by , is a Majorana fermion, singlet under the Standard Model gauge group. Then, the part of the Lagrangian involving only the sterile neutrino field, including operators up to dimension six, reads:
Besides, there exist portal interaction terms between the visible and the dark sectors, such as the Yukawa coupling , where is a lepton doublet and , with the Standard Model Higgs doublet. We will assume in what follows that all portal interactions with the Standard Model have negligible strength, such that they play no role in dark matter production and ensure its stability on cosmological time scales.
The dimension-six self-interaction term induces dark matter -to- scatterings, with a cross section
and which is constrained to be from Bullet Cluster observations ; this constraint translates into the lower limit . The -to- process has been invoked to alleviate the small scale problems of the cold dark matter paradigm (see e.g. ), however it does not change the particle number and does not play any role in dark matter production. On the other hand, the 4-to-2 process , induced by the diagrams shown in Fig. 1, does change the particle number and determines the epoch at which the sterile neutrinos can no longer be maintained in thermal equilibrium, thus setting their relic abundance.
We assume in what follows that dark matter particles were produced in the early stages of the Universe and that the strength of the dark matter self-interaction was sufficiently large to keep the dark plasma in thermal equilibrium, characterized by the temperature . Since we assume negligible portal interactions with the Standard Model plasma, the dark sector temperature, , is in general different from the visible sector temperature, .
The evolution of the dark matter number density, which we denote by , is dictated by the following Boltzmann equation:
where is the Hubble parameter, is the thermally averaged cross section for the 4-to-2 process and is the dark matter number density in equilibrium, all evaluated at the corresponding cosmic time . The Hubble parameter is given by:
with the gravitational constant, while and are respectively the energy densities in the visible and dark sectors:
with and the effective number of degrees of freedom in the corresponding sector at that cosmic epoch [27, 28]. Besides, the dark matter equilibrium number density reads, assuming Maxwell-Boltzmann statistics:
with the modified Bessel function of the second kind of integer order . Finally, the thermally averaged cross section for the 4-to-2 process is defined as 
where are the velocities of the initial state particles, while the -to- annihilation cross section can be obtained from the invariant amplitude from
where is the solid angle of any of the final state particles.
We have calculated the invariant amplitude from the Feynman diagrams depicted in Fig. 1, using FeynCalc [30, 31]. The final expression is rather lengthy and will not be shown here. Expanding in terms of the dark matter velocities (), one finds that the leading terms are proportional to fourth-order invariants of the velocities (e.g. , , ), namely the cross-section is -wave suppressed. This behavior can be easily understood from the Pauli exclusion principle. Each initial fermionic state is specified by its spin and orbital angular momentum . In a partial wave expansion, the state where all the particles have is incompatible with the Pauli exclusion principle, since this would require two of the fermions to have also the same spin quantum number. Then, the lowest order term in the partial wave expansion must contain two fermions with and two fermions with , thus giving the -wave suppression in the cross section.111The -wave suppression of the 4-to-2 cross section for Majorana fermions could be circumvented if the fermions have an additional quantum number (for instance, if they are charged under an unbroken global symmetry), or if they have higher spin. This argument can be generalized to other number changing processes involving self-interactions of spin 1/2 (or spin 3/2) Dirac or Majorana fermions.
Using the results from  for the thermal average of invariants involving velocities, we find for the thermally averaged cross section:
We note that and are functions of the temperature of the dark sector, whereas the Hubble parameter depends both on and . On the other hand, in the frameworks of relevance for our analysis and therefore the -dependence of the Hubble parameter can be neglected.
The dark matter relic abundance can be calculated using the instantaneous freeze-out approximation, which states that the dark matter yield, defined as the dark matter number density normalized to the entropy density in the visible sector, remains constant after freeze-out, defined by the time when the rate for the 4-to-2 process equals the Hubble parameter:
where the rate is given by
Here, we have used the fact that the thermally averaged cross section is -wave suppressed, and hence proportional to . We remark that by expanding the cross section in powers of the velocity we implicitly require that the dark matter freezes-out while already non-relativistic, which in turn requires a sufficiently low freeze-out temperature. To ensure the validity of our approximations we will conservatively require .
Finally, the dark matter abundance today can be calculated from
where and are the present time visible sector entropy density and critical density, while is the present time dark matter yield, which in the instantaneous freeze-out approximation is equal to the dark matter yield at the freeze-out epoch:
From Eqs. (10,11,12), and from requiring that the sterile neutrino relic abundance reproduces the dark matter abundance measured by Planck, , one obtains in the instantaneous freeze-out approximation the following approximate value of the dark sector freeze-out temperature:222This is an excellent approximation for freeze-out via 4-to-2 annihilations, due to the strong suppression of the reaction rate by a factor when the dark matter is non-relativistic. For comparison, the rate of 2-to-2 annihilations, relevant for the freeze-out of Weakly Interacting Massive Particles (WIMPs), is only suppressed by . We have confirmed the goodness of the approximation by solving exactly the Boltzmann equations for some selected values of the parameters, obtaining a result which differs by at most 10% from the one obtained with the instantaneous freeze-out approximation.
which depends on the dark matter mass and on the ratio between the temperatures of the visible and dark sectors at freeze-out. For and GeV , we find . Besides, we obtain for the thermally averaged cross section at the epoch :
The value of the suppression scale leading to the measured dark matter abundance, is shown in Fig. 2 for different values of the ratio of the freeze-out temperatures . In our analysis we conservatively disregard the region of the parameter space where , shown as an orange region, to ensure non-relativistic freeze-out, and the region where , shown as a gray region, to ensure the validity of the effective Lagrangian Eq. (1).
The remaining parameter space is constrained by astronomical observations, concretely by the effect of the self-interactions in the Bullet Cluster, which exclude the region of the parameter space where , shown as brown in the Figure. We find that the observed dark matter abundance can be reproduced in this framework, provided the dark matter mass is in the range GeV, which respectively correspond to a scale of the dimension-six operator, normalized to the dark matter mass, in the range and to a visible-to-hidden sector temperature ratio at freeze-out .
This model has observable consequences through the dark matter self-interactions. In particular, the self-interaction could provide a solution to the small scale structure problems if ; this region of the parameter space is shown as green in the Figure and lies in the region where the correct dark matter abundance is set by the 4-to-2 annihilations.
Additional signals of the model arise when the portal interactions to the Standard Model have non-vanishing strength. More specifically, the Yukawa coupling to the lepton doublet induces dark matter decay, which leads to a flux of cosmic gamma-rays, antimatter particles and neutrinos which, if sufficiently intense, could be disentangled from the background fluxes. Conversely, the non-observation of a statistically significant excess in the measured fluxes can be translated into upper limits on the fundamental parameters of the model. For instance, the non-observation of a significant monoenergetic signal in the gamma-ray sky (see e.g. ) can be translated into the following limit on the Yukawa coupling:
for . Such small values of the Yukawa coupling could be naturally accommodated in the model by extending its symmetry by an exact, or very mildly broken, discrete symmetry, under which the Standard Model fields are even while the sterile neutrino is odd. In the case that the symmetry is mildly broken, signals of the model could be observed with future instruments, such as the gamma-ray telescope e-ASTROGAM, which aims to increase the current sensitivity to the decay width for by more than one order of magnitude . The upper limit Eq. (18) also implies that the sterile neutrino dark matter candidate contributes negligibly to the active neutrino masses. Concretely, we find . Neutrino masses could on the other hand be generated by neutrino Yukawa interactions with other fermion singlets, or perhaps by other mechanisms. The observation of dark matter self interactions in astronomical objects, accompanied by the observation of a dark matter decay signal, notably a gamma-ray line in the energy range GeV, would provide support to the framework where the dark matter is constituted by sterile neutrinos as SIMPs.
An important question is whether the strength of the effective self-interaction which is necessary to reproduce the observed dark matter abundance can be naturally accommodated in an ultraviolet complete model. We will address this question in the next section with a simple toy model.
3 A toy model of dark matter self-interactions
We extend the particle content of the dark sector by one scalar singlet , with mass , that interacts with the sterile neutrino through a Yukawa coupling of the form:
Then, at the energy scales where the scalar can be integrated out, one recovers the effective interaction in Eq. (1), with . In this toy model, the Bullet Cluster constraint translates into the upper limit on the Yukawa coupling .
The process is induced in this model by the diagrams in Fig. 3. Away from resonances ( and ), the thermally averaged annihilation cross section reads:
where and the coefficients are given by
The value of the Yukawa coupling required to reproduce the observed dark matter abundance can be calculated for a given dark matter mass and mediator mass, and for a given value of the ratio of temperatures at freeze-out of the hidden and visible sectors, using Eqs. (20) and (16). The resulting values of the Yukawa coupling as a function of the dark matter mass are shown in Fig. 4, for the concrete cases (left panel) and 10 (right panel), for . The regions of the parameter space where our approach cannot be applied are indicated in Fig. 4 as an orange region and a gray region, and are bounded by the requirement that the dark matter freezes-out non-relativistically () and that perturbation theory remains valid (). Furthermore, the region indicated in brown implies a strength for the dark matter self-interaction in conflict with observations of the Bullet Cluster.
We find that the observed dark matter abundance can be reproduced by the thermal freeze-out of sterile neutrinos from the equilibrium density due to their self-interactions only, if the sterile neutrino mass is in the range (), the Yukawa coupling is (0.4) and the visible-to-hidden sector temperature ratio at freeze-out is for . We also find points in the parameter space where the small scale problems of the cold dark matter paradigm can be alleviated by dark matter self-interactions, which are shown in Fig. 4 as a green region.
We note that for large the effective theory limit is recovered, and therefore the right panel of Fig. 4 is a just a recast of Fig. 2 with . In contrast, for moderate the effective Lagrangian does not provide a good description of the model and accordingly the parameters required to reproduce the observed dark matter abundance differ. Concretely, the dark matter mass window widens both at low and high dark matter masses, and allows a wider range of values for , as illustrated in Fig. 2, left panel. We also note that for large the 4-to-2 cross section scales as . Therefore, in order to have a sufficiently large annihilation cross section, an increase in must be compensated by an increase in . As a result, for larger and larger the lower limit on correspondingly increases, eventually entering into conflict with our perturbativity requirement , thus setting the upper limit on the singlet scalar mass .
The presence of the singlet scalar in the spectrum opens new portals between the visible and the dark sector, through the couplings with the Standard Model Higgs and . These couplings, while not directly involved in the freeze-out of the 4-to-2 annihilation process, can be important to address the genesis of the dark sector plasma and hence provide a complete thermal history of the dark matter of our Universe, and in particular of the relation between the temperatures of the visible and dark sector plasmas.333This is in contrast to WIMP dark matter, where the same interaction that sets the final relic abundance is responsible for thermalizing the dark matter candidate with the Standard Model bath, erasing all previous history of the dark matter population.
More specifically, after electroweak symmetry breaking the particle content of the model contains a real scalar, , which is mainly composed of the Standard Model Higgs boson and has a small component of the scalar singlet , which we parametrize by the mixing angle . The Higgs boson is in general in thermal equilibrium with the visible sector plasma and populates the dark sector via freeze-in from the decays , induced by the interaction terms and . The temperature evolution of the dark sector energy density is described by the following Boltzmann equation :
which assumes that all dark sector particles are relativistic while freeze-in production is efficient, namely the epoch . In this Boltzmann equation, is the visible sector energy density, given in Eq. (5), is the Higgs number density, which at temperatures is well approximated by its equilibrium value,
and is the total decay rate into hidden sector particles:
Finally, one can calculate the ratio of dark-to-visible energy densities from solving Eq. (21), taking at high temperature as boundary condition. Correspondingly, the ratio of dark-to-visible sector temperatures can be calculated from Eq. (5).
The temperature ratio is time dependent. It is then convenient to work instead with the entropy ratio, which stays constant over time (as long as there is no entropy production nor energy transfer between the visible and dark sectors) and in particular stays constant during the freeze-out process [23, 16]. More specifically, we define the entropy ratio as , where is the entropy density of the visible sector, given in Eq. (14), and is the corresponding quantity for the dark sector. We calculate the value of the entropy ratio by evaluating at the temperature , when the freeze-in production is essentially completed. We obtain:
where was calculated using Eq. (21), and we have taken and at GeV.
Using the fact that the entropy density remains constant, it is straightforward to calculate from Eq. (13) the dark matter yield at the present time in the instantaneous freeze-out approximation. The result is
with .444The definition of follows from the effective number of degrees of freedom contributing to the hidden sector entropy density at freeze-out [27, 28]. Here, is determined by Eq. (10) and depends on the value of the thermally averaged cross-section for the 4-to-2 annihilation process, while is related to by the entropy ratio. Finally, the sterile neutrino relic abundance follows from Eq. (12), and can be readily confronted with the dark matter abundance measured by Planck.
We show in Fig. 5 the analog of Fig. 4, but plotting contours of the portal coupling or which is necessary to reproduce the observed dark matter abundance via freeze-out of the 4-to-2 annihilations (assuming that only one of them is non-vanishing), in a framework where the dark sector is populated via freeze-in from Higgs decays. In this manner, the viability of one point in the parameter space is dictated by fundamental parameters of the theory (namely or ), rather than the thermodynamical quantity . We restrict our analysis to values of the dark matter mass GeV to ensure that the hidden sector particles are relativistic during freeze-in, so that the Boltzmann equation Eq. (21) holds. Under these assumptions, and imposing the restrictions on the strength of the self-interaction from Bullet Cluster observations, we find that reproducing the observed dark matter abundance requires , or , for .
We have considered a scenario where the dark matter is constituted by Majorana sterile neutrinos with sizable self-interactions, which we describe by a dimension-six operator suppressed by the scale , and with negligible interactions with the Standard Model particles. We have assumed that at very early times the plasma of sterile neutrinos is in a thermalized state, with a temperature which is in general different to the temperature of the visible sector plasma, and we have calculated the sterile neutrino relic abundance from the freeze-out of the number changing 4-to-2 annihilation induced by the self-interaction.
We have shown that, for reasonable values of the model parameters, sterile neutrinos can account for all the dark matter of the Universe. The framework is further constrained by the effect of the self-interactions on astronomical objects, and is partly excluded by observations of the Bullet Cluster. Imposing that the sterile neutrinos freeze-out when already non-relativistic, our results indicate that their mass is GeV, that the suppression scale of the dimension-six operator is , and that the ratio of temperatures between the visible and dark sectors at freeze-out is . Furthermore, for some choices of the model parameters, the self-interactions can be strong enough to alleviate the small scale problems of the cold dark matter paradigm.
Opening the portal interactions to the Standard Model allows further tests of the model. Specifically, the feeble charged current interaction to the charged leptons would induce dark matter decay. Therefore, the observation of dark matter self interactions in astronomical objects, accompanied by the observation of a dark matter decay signal, notably a gamma-ray line in the energy range GeV, would provide support to our framework.
We have also analyzed the scenario where the self-interaction is induced by the exchange of a scalar singlet field. We have found that reproducing the observed dark matter relic abundance restricts the Yukawa coupling between the singlet scalar and the sterile neutrinos to be larger than 0.08 (0.4) and the sterile neutrino mass to be in the range () when the scalar mass is 3 (10) times larger than the sterile neutrino mass. This toy model also contains the necessary elements to construct a complete thermal history of the dark sector, as the portal interactions of the singlet scalar with the Standard Model Higgs can be responsible for populating the dark sector via freeze-in and for determining the temperature of the dark sector plasma.
This work has been partially supported by the DFG cluster of excellence EXC 153 “Origin and Structure of the Universe” and by the Collaborative Research Center SFB1258. T.T. acknowledges support from JSPS Fellowships for Research Abroad. Numerical computation in this work was carried out at the Yukawa Institute Computer Facility. T.T. thanks Camilo Garcia Cely for fruitful discussions.
-  P. Minkowski, at a Rate of One Out of Muon Decays?, Phys. Lett. 67B (1977) 421–428.
-  R. N. Mohapatra and G. Senjanovic, Neutrino Mass and Spontaneous Parity Violation, Phys. Rev. Lett. 44 (1980) 912.
-  T. Yanagida, HORIZONTAL SYMMETRY AND MASSES OF NEUTRINOS, Conf. Proc. C7902131 (1979) 95–99.
-  M. Gell-Mann, P. Ramond and R. Slansky, Complex Spinors and Unified Theories, Conf. Proc. C790927 (1979) 315–321, [1306.4669].
-  J. Schechter and J. W. F. Valle, Neutrino Masses in SU(2) x U(1) Theories, Phys. Rev. D22 (1980) 2227.
-  S. Dodelson and L. M. Widrow, Sterile-neutrinos as dark matter, Phys. Rev. Lett. 72 (1994) 17–20, [hep-ph/9303287].
-  S. Horiuchi, P. J. Humphrey, J. Onorbe, K. N. Abazajian, M. Kaplinghat and S. Garrison-Kimmel, Sterile neutrino dark matter bounds from galaxies of the Local Group, Phys. Rev. D89 (2014) 025017, [1311.0282].
-  D. Malyshev, A. Neronov and D. Eckert, Constraints on 3.55 keV line emission from stacked observations of dwarf spheroidal galaxies, Phys. Rev. D90 (2014) 103506, [1408.3531].
-  X.-D. Shi and G. M. Fuller, A New dark matter candidate: Nonthermal sterile neutrinos, Phys. Rev. Lett. 82 (1999) 2832–2835, [astro-ph/9810076].
-  A. Kusenko, Sterile neutrinos, dark matter, and the pulsar velocities in models with a Higgs singlet, Phys. Rev. Lett. 97 (2006) 241301, [hep-ph/0609081].
-  K. Petraki and A. Kusenko, Dark-matter sterile neutrinos in models with a gauge singlet in the Higgs sector, Phys. Rev. D77 (2008) 065014, [0711.4646].
-  A. Merle, V. Niro and D. Schmidt, New production mechanism for kev sterile neutrino dark matter by decays of frozen-in scalars, JCAP 1403 (2014) 028, [1306.3996].
-  J. König, A. Merle and M. Totzauer, keV Sterile Neutrino Dark Matter from Singlet Scalar Decays: The Most General Case, JCAP 1611 (2016) 038, [1609.01289].
-  A. Schneider, Astrophysical constraints on resonantly produced sterile neutrino dark matter, JCAP 1604 (2016) 059, [1601.07553].
-  S. Tulin and H.-B. Yu, Dark Matter Self-interactions and Small Scale Structure, 1705.02358.
-  E. D. Carlson, M. E. Machacek and L. J. Hall, Self-interacting dark matter, Astrophys. J. 398 (1992) 43–52.
-  Y. Hochberg, E. Kuflik, T. Volansky and J. G. Wacker, Mechanism for Thermal Relic Dark Matter of Strongly Interacting Massive Particles, Phys. Rev. Lett. 113 (2014) 171301, [1402.5143].
-  Y. Hochberg, E. Kuflik, H. Murayama, T. Volansky and J. G. Wacker, Model for Thermal Relic Dark Matter of Strongly Interacting Massive Particles, Phys. Rev. Lett. 115 (2015) 021301, [1411.3727].
-  Y. Hochberg, E. Kuflik and H. Murayama, SIMP Spectroscopy, JHEP 05 (2016) 090, [1512.07917].
-  N. Bernal, C. Garcia-Cely and R. Rosenfeld, WIMP and SIMP Dark Matter from the Spontaneous Breaking of a Global Group, JCAP 1504 (2015) 012, [1501.01973].
-  S.-M. Choi and H. M. Lee, SIMP dark matter with gauged Z symmetry, JHEP 09 (2015) 063, [1505.00960].
-  S.-Y. Ho, T. Toma and K. Tsumura, A Radiative Neutrino Mass Model with SIMP Dark Matter, JHEP 07 (2017) 101, [1705.00592].
-  N. Bernal and X. Chu, SIMP Dark Matter, JCAP 1601 (2016) 006, [1510.08527].
-  M. Heikinheimo, T. Tenkanen, K. Tuominen and V. Vaskonen, Observational Constraints on Decoupled Hidden Sectors, Phys. Rev. D94 (2016) 063506, [1604.02401].
-  X. Chu and C. Garcia-Cely, Self-interacting Spin-2 Dark Matter, Phys. Rev. D96 (2017) 103519, [1708.06764].
-  S. W. Randall, M. Markevitch, D. Clowe, A. H. Gonzalez and M. Bradac, Constraints on the self-interaction cross-section of dark matter from numerical simulations of the merging galaxy cluster 1e 0657-56, Astrophys. J. 679 (2008) 1173–1180, [0704.0261].
-  P. Gondolo and G. Gelmini, Cosmic abundances of stable particles: Improved analysis, Nucl. Phys. B360 (1991) 145–179.
-  L. Husdal, On Effective Degrees of Freedom in the Early Universe, Galaxies 4 (2016) 78, [1609.04979].
-  S.-M. Choi, H. M. Lee and M.-S. Seo, Cosmic abundances of SIMP dark matter, JHEP 04 (2017) 154, [1702.07860].
-  R. Mertig, M. Bohm and A. Denner, FEYN CALC: Computer algebraic calculation of Feynman amplitudes, Comput. Phys. Commun. 64 (1991) 345–359.
-  V. Shtabovenko, R. Mertig and F. Orellana, New Developments in FeynCalc 9.0, Comput. Phys. Commun. 207 (2016) 432–444, [1601.01167].
-  Planck collaboration, P. A. R. Ade et al., Planck 2015 results. XIII. Cosmological parameters, Astron. Astrophys. 594 (2016) A13, [1502.01589].
-  R. Essig, E. Kuflik, S. D. McDermott, T. Volansky and K. M. Zurek, Constraining Light Dark Matter with Diffuse X-Ray and Gamma-Ray Observations, JHEP 11 (2013) 193, [1309.4091].
-  e-ASTROGAM collaboration, A. De Angelis et al., The e-ASTROGAM mission, Exper. Astron. 44 (2017) 25–82, [1611.02232].
-  X. Chu, T. Hambye and M. H. G. Tytgat, The Four Basic Ways of Creating Dark Matter Through a Portal, JCAP 1205 (2012) 034, [1112.0493].