## MeV scale model of SIMP dark matter, neutrino mass and leptogenesis

## Subhendra Mohanty, a Ayon Patra, a Tripurari Srivastava a

a Theoretical Physics Division, Physical Research Laboratory, Ahmedabad-380009, India

E-mail: mohanty@prl.res.in , ayon@okstate.edu , tripurari@prl.res.in

Abstract: We consider a simple extension of the Standard Model with two singlet scalar fields and three heavy right-handed neutrinos. One of the scalar fields serves as an MeV scale dark matter and its stability is ensured by the introduction of an extra Z 2 symmetry. The second scalar (which is even under the Z 2 symmetry) generates the mass term of the scalar, contributes to the 3 → 2 annihilation process required for the correct relic density of the dark matter and it also contributes to the leptogenesis. The right-handed neutrinos are responsible for the generation of light neutrino masses through Type-I seesaw mechanism. The decay of the heavy right-handed neutrino can generate the lepton asymmetry which can then be converted to baryon asymmetry through sphaleron transitions.

|   Contents | Contents                       |    |
|------------|--------------------------------|----|
|          1 | Introduction                   |  2 |
|          2 | Model                          |  3 |
|          3 | Relic abundance                |  5 |
|          4 | Self interaction               |  9 |
|          5 | Neutrino mass generation       | 12 |
|          6 | Leptogenesis                   | 14 |
|          7 | Neutrinoless Double Beta Decay | 16 |
|          8 | Conclusion                     | 17 |

## 1 Introduction

The Weakly interacting massive particle (WIMP) [1, 2] paradigm of dark matter (DM) is severely constrained due to the non-observation of any signal in DM direct detection experiments [3-5]. This calls for an extension of the DM parameter space to a lower mass range which is not very well probed in direct detection experiments. The amount of dark matter in the universe is also well measured from various observations, e.g., CMB [6] and large scale structure surveys [7]. At the scale of galaxies, the cold dark matter model is in good agreement with the observation. However there are certain problems like lesser observation of dwarf galaxies (missing satellite problem) [8, 9] and the absence of very bright galaxies (too big too fail problem) [10, 11] which are both predicted from N-body simulations. Simulations also predict that in the galactic distribution of dark matter, there is a large concentration at the core which is again not observed in rotation curves of galaxies (core cusp problem) [12, 13]. These problems of ΛCDM model can be solved by introducing large self-interactions in the dark matter [14-16].

The observed relic density, which in the WIMP paradigm, could be explained by decoupling of the 2 → 2 annihilation no longer works for MeV scale DM. On the other hand, the abundance of an MeV scale DM naturally satisfies the relic density constraints if there exists a 3 → 2 mechanism for DM annihilation [17-32]. This requires a large self-interaction among the DM particles and is popularly known as the strongly interacting massive particle (SIMP) DM scenario [33-37]. In the usual 3 → 2 mechanism, three DM particles can annihilate into two DM particles in the final state. In our model, we have invoked a Z 2 symmetry to ensure the stability of DM ( φ ), which however prevents this 3 φ to 2 φ annihilation process. We therefore introduce a second scalar ( δ ) which assists in obtaining the correct relic density via a φφφ → φδ annihilation process. We write the most general potential involving the scalars φ and δ which respects the Z 2 -odd nature of φ . In this potential, by giving δ a vacuum expectation value, we generate a mass of the φ and the δ in the ∼ 10 MeV scale. Although the scalar δ does not contribute to the relic density, it can interact with the right-handed neutrinos to generate the observed baryon asymmetry of the universe through low scale leptogenesis [38-42]. In this model we have introduced three right-handed neutrinos (two of them with masses around 10 TeV and one much lighter with mass around 10 MeV) to obtain the light neutrino masses by Type-I seesaw mechanism [43-48] while also being able to generate the required lepton asymmetry. The decay process of the heavy right-handed neutrino can produce the lepton asymmetry through the interference of the tree-level and the one-loop diagram involving the lighter right-handed neutrino and the scalar δ .

The rest of the paper is organized as follows. In Sec. 2 we introduce the model and present a detailed analysis of the scalar sector. The calculation of the DM relic density is performed in Sec. 3 where we scan over the allowed parameters which provide the correct experimentally observed relic abundance. In Sec. 4 we look at the constraints on the selfinteraction cross-section for a SIMP DM arising from various astrophysical observations and study their effect on our allowed parameter space. A detailed analysis of the neutrino sector is performed in Sec. 5. Our scan results show that a clear hierarchy among the

Table 1 : Particle spectrum for the Model. Here i is the generation index.

| Field                   |   SU (3) C |   SU (2) L | U (1) Y   | Z 2   |
|-------------------------|------------|------------|-----------|-------|
| Q L i =   u d   L   |          3 |          2 | 1 3       | +     |
| u R i                   |          3 |          1 | 4 3       | +     |
| d R i                 |          3 |          1 | - 2 3     | +     |
| l L i =   ν e   L   |          1 |          2 | -1        | +     |
| e R i                   |          1 |          1 | -2        | +     |
| N i                     |          1 |          1 | 0         | +     |
| H =    H + H 0    |          1 |          2 | 1         | +     |
| φ                       |          1 |          1 | 0         | -     |
| δ                       |          1 |          1 | 0         | +     |

neutrino Dirac masses is required to satisfy the experimentally observed neutrino masssquared differences and mixings. In Sec. 6 we study the mechanism to achieve low scale leptogenesis in our framework. Sec. 7 includes a discussion on the possible constraints arising from neutrino-less double beta decay processes. We present our conclusions in Sec. 8.

## 2 Model

We extend the SM framework with two singlet scalars and three right-handed neutrinos. The gauge group of the SM remains unchanged but we need to introduce a discrete Z 2 symmetry to ensure the stability of the dark matter. The particle spectrum of the model is given in Table 1. The quarks and leptons remain the same as in the SM but we introduce three right-handed neutrinos ( N i ) which help in generating the light neutrino masses and mixings via Type-I seesaw mechanism. Similar to the SM, the Higgs doublet H is responsible for the electroweak symmetry breaking and generation of quark and lepton masses. As can be seen from Table 1, the scalar singlet φ is odd under the Z 2 symmetry while all other particles are even. The scalar φ thus cannot decay and becomes the DM candidate in this model. The other scalar singlet δ can acquire a non zero vacuum expectation value (VEV) and generate the mass for the DM. The scalar δ actually has multiple roles in this

model - in addition to contributing to the DM mass, it helps in obtaining the correct DM relic density by assisting in the 3 → 2 DM annihilation process and also plays an important role in generating the lepton asymmetry. One of the heavy right-handed neutrinos can decay into a light fermion and a scalar boson directly through a tree level process but this cannot produce any Charge-Parity (CP) asymmetry. The one-loop corrections to this process includes vertex correction and self-energy corrections, and their interference with the tree-level decay process give rise to the required CP asymmetry. This will be discussed in more details in section 6.

The non-zero vacuum expectation values for the scalar fields are

$$\langle H \rangle = v _ { H }, \ \langle \delta \rangle = v _ { \delta },$$

while the scalar potential in the model is given as:

$$\mathcal { L } & = \frac { \lambda _ { 1 1 } } { 4! } \phi ^ { 4 } + \frac { \mu _ { \phi } ^ { 2 } } { 2 } \phi ^ { 2 } + \frac { \lambda _ { 1 2 } } { 4 } \phi ^ { 2 } \delta ^ { 2 } + \frac { \mu _ { \delta } ^ { 2 } } { 2 } \delta ^ { 2 } + \frac { \mu _ { 2 2 } } { 3! } \delta ^ { 3 } + \frac { \lambda _ { 2 2 } } { 4! } \delta ^ { 4 } + \frac { \mu _ { 2 1 } } { 2 } \phi ^ { 2 } \delta } \\ & + \frac { \lambda _ { 1 3 } } { 2 } \phi ^ { 2 } H ^ { \dagger } H + \mu _ { 2 3 } \delta H ^ { \dagger } H + \frac { \lambda _ { 2 3 } } { 2 } \delta ^ { 2 } H ^ { \dagger } H + \mu _ { H } ^ { 2 } H ^ { \dagger } H + \lambda _ { 3 3 } \left [ H ^ { \dagger } H \right ] ^ { 2 }. \\ \cdots & \quad \cdots \, \dots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \end{$$

Measurement of the Higgs invisible decay branching ratio puts stringent constraints on some of these couplings. Recent results from the LHC gives an upper limit on the invisible decay of the Higgs boson as BR( H → invisible) &lt; 0 19 [49]. . In our model, the process contributing to it comes from the H → φφ channel. Calculating the decay width of this channel gives:

$$\Gamma ( H \rightarrow \phi \phi ) = \frac { \lambda _ { 1 3 } ^ { 2 } v _ { H } ^ { 2 } } { 3 2 \pi m _ { H } } \left ( 1 - \frac { 4 m _ { \phi } ^ { 2 } } { m _ { H } ^ { 2 } } \right ) ^ { \frac { 1 } { 2 } }.$$

The numerical value of λ 13 &lt; O (10 -2 ) ensures the Higgs invisible decay constraint are satisfied.

The scalar field φ , being odd under Z 2 , does not mix with the other fields and becomes the DM candidate with its mass given as

$$M _ { \phi } ^ { 2 } = \left ( \frac { \lambda _ { 1 2 } } { 2 } v _ { \delta } ^ { 2 } + 2 \lambda _ { 1 3 } v _ { H } ^ { 2 } + \mu _ { 1 2 } v _ { \delta } + 2 \mu _ { \phi } ^ { 2 } \right ).$$

In order to obtain the correct relic density, M φ must be in the MeV range with λ 12 ∼ 1. Hence we choose v δ to be of order MeV as well. The value of v H , on the other hand, is required to be v H ≡ v EW = 174 GeV. This forces λ 13 to be extremely small ( glyph[lessorsimilar] 10 -6 ) as otherwise a large fine-tuning will be required to obtain the low DM mass.

The tiny λ 13 coupling also plays an important role in evading the DM direct detection experimental constraints. Current dark matter direct detection experiments are based on nucleon recoil energy measurement from their scattering with the DM [50]. However, for sub GeV dark matter, these bounds are quite weak although recent experiments, e.g., CRESST-III [51], Xenon1T [52] and CDEX-1B [53] etc. provide limit for some range in sub GeV DM. The small λ 13 coupling results in the DM scattering with the quarks (via a Higgs boson) to be highly suppressed compared to the limit from direct detection experiments. There are other experiments, e.g., Xenon10 [54], super CDMS [55] and SENSEI [56] that

use DM electron scattering to also provide constraints over cross-section. Since the Higgs couplings to the DM particle and to the electrons are both extremely small, we can easily satisfy the experimental bounds in this case as well.

The neutral component of the doublet scalar H 0 and the singlet δ fields mix to form a 2 × 2 mass-squared matrix given as

$$M _ { H } ^ { 2 } = \left [ \begin{array} { c } \frac { 1 } { 2 } \lambda _ { 2 2 } v _ { \delta } ^ { 2 } + 2 \lambda _ { 2 3 } v _ { H } ^ { 2 } + 6 \mu _ { 2 2 } v _ { \delta } + 2 \mu _ { \delta } ^ { 2 } & - \frac { v _ { \delta } } { 3 \sqrt { 2 } v _ { H } } \left ( \lambda _ { 2 2 } v _ { \delta } ^ { 2 } + 1 8 \mu _ { 2 2 } v _ { \delta } + 1 2 \mu _ { \delta } ^ { 2 } \right ) \\ - \frac { v _ { \delta } } { 3 \sqrt { 2 } v _ { H } } \left ( \lambda _ { 2 2 } v _ { \delta } ^ { 2 } + 1 8 \mu _ { 2 2 } v _ { \delta } + 1 2 \mu _ { \delta } ^ { 2 } \right ) & 4 \lambda _ { 3 3 } v _ { H } ^ { 2 } \end{array} \right ].$$

As can clearly be seen, the mixing between the H 0 and the δ field is proportional to v δ v H and hence is extremely small.

## 3 Relic abundance

The SM Higgs mixings with the singlet scalars in our model remain very small as was discussed in Section. 2. The additional scalars φ and δ thus remain largely secluded. The relic abundance of dark matter can be computed via annihilation of dark matter φ through 3 → 2 process. Choosing m φ &lt; m δ &lt; 2 m φ and imposing the discrete Z 2 symmetry, the dominant contribution for the annihilation of φ is through φφφ → φδ . In addition to this, the δδ → φφ process also needs to be included to compute the evolution of number density of φ . The Feynman diagrams for the relevant processes contributing to the relic abundance are given in the Fig. 2.

Figure 1 : Feynman diagram for processes contributing to the calculation of relic density.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00909v1-figure-07ca8eac-5e20-413a-885f-d080fc1357b8.jpg)

Evolution of number density of dark matter is governed by the set of Boltzmann's equations given by:

$$\frac { d n _ { \phi } } { d t } + 3 H n _ { \phi } = - \frac { 1 } { 3 } \left \langle \sigma v ^ { 2 } \right \rangle _ { 3 \rightarrow 2 } \left ( n _ { \phi } ^ { 3 } - n _ { \phi } n _ { \delta } n _ { \phi } ^ { e q } \right ) + \left \langle \sigma v \right \rangle _ { 2 \rightarrow 2 } \left ( n _ { \phi } ^ { 2 } - n _ { \delta } ^ { 2 } \right ) \text{ \quad \ } \text{(3.1)}$$

$$\widetilde { \frac { d n _ { \delta } } { d t } } + 3 H n _ { \delta } = + \frac { \breve { 1 } } { 6 } \left \langle \sigma v ^ { 2 } \right \rangle _ { 3 \rightarrow 2 } ( n _ { \phi } ^ { 3 } - n _ { \phi } n _ { \delta } n _ { \phi } ^ { e q } ) - \left \langle \sigma v \right \rangle _ { 2 \rightarrow 2 } ( n _ { \phi } ^ { 2 } - n _ { \delta } ^ { 2 } ). \text{ \quad } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{\quad } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \end{$$

Here the number densities are defined as:

$$n _ { i } = g _ { i } \int \frac { d ^ { 3 } p _ { i } } { ( 2 \pi ) ^ { 3 } } f _ { i },$$

where f i is the distribution function for i-th particle. The factors in front of thermal average cross-sections are to avoid double counting (see [57]). We can perform a change of variables to solve the equations by assuming Y i = n /s i and x = m /T φ where s is the entropy density given as (2 π g T 2 ∗ 3 / 45) and T is the temperature. The number density evolution equation can be written as:

$$\frac { d Y _ { \phi } } { d x } = - \frac { x s ^ { 2 } } { 3 H ( m _ { \phi } ) } \left \langle \sigma v ^ { 2 } \right \rangle _ { 3 \to 2 } ( Y _ { \phi } ^ { 3 } - Y _ { \phi } Y _ { \delta } Y _ { \phi } ^ { e q } ) + \frac { x s } { H ( m _ { \phi } ) } \left \langle \sigma v \right \rangle _ { 2 \to 2 } ( Y _ { \phi } ^ { 2 } - Y _ { \delta } ^ { 2 } ) \quad ( 3. 4 )$$

$$\frac { d Y _ { \delta } } { d x } = + \frac { \stackrel { \cdot } { x s ^ { 2 } } } { 6 H ( m _ { \phi } ) } \left \langle \sigma v ^ { 2 } \right \rangle _ { 3 \to 2 } ( Y _ { \phi } ^ { 3 } - Y _ { \phi } Y _ { \delta } Y _ { \phi } ^ { e q } ) - \frac { \cdot \cdot \cdot } { H ( m _ { \phi } ) } \left \langle \sigma v \right \rangle _ { 2 \to 2 } ( Y _ { \phi } ^ { 2 } - Y _ { \delta } ^ { 2 } ), \quad ( 3. 5 )$$

where H m ( φ ) = √ π g 2 ∗ 90 m 2 φ M pl .

Thermal average cross-section can be computed from the model as follows:

$$\langle \sigma _ { 3 \to 2 } v ^ { 2 } \rangle = \frac { 1 } { n _ { \phi } ^ { 3 } } \int d \Pi _ { \phi } d \Pi _ { \phi } d \Pi _ { \phi } d \Pi _ { \delta } d \Pi _ { \phi } ( 2 \pi ) ^ { 4 } \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

$$\delta ^ { 4 } ( p _ { 1 } ( \phi ) + p _ { 2 } ( \phi ) + p _ { 3 } ( \phi ) - p _ { 4 } ( p _ { \delta } ) - p _ { 5 } ( \phi ) ) f _ { \phi } ^ { 3 } | \mathcal { M } _ { \phi \phi \phi \rightarrow \delta \delta } | ^ { 2 }, \quad ( 3. 7 )$$

where Π i = d p 3 i (2 π ) 3 2 E i .

For the 3 → 2 process, the thermal average cross-section can be written as [58]:

$$\langle \sigma _ { 3 \rightarrow 2 } v ^ { 2 } \rangle = \frac { \lambda _ { 1 1 } ^ { 2 } \mu _ { \text{eff} } ^ { 2 } } { 6 4 \pi m _ { \phi } ^ { 3 } } \sqrt { \left ( 1 - \frac { ( m _ { \phi } + m _ { \delta } ) ^ { 2 } } { 9 m _ { \phi } ^ { 2 } } \right ) \left ( 1 - \frac { ( m _ { \delta } - m _ { \phi } ) ^ { 2 } } { 9 m _ { \phi } ^ { 2 } } \right ) } \frac { 1 } { 6 4 m _ { \phi } ^ { 4 } },$$

where µ eff = µ 12 + λ 12 2 v δ . Here, we have assumed φ to be non-relativistic. We solve the coupled Boltzmann equations numerically to evaluate the yield Y φ at freeze out temperature. The relic abundance can be obtained from it and is given by:

$$\Omega _ { \phi } h ^ { 2 } \approx 2. 7 5 5 \times 1 0 ^ { 8 } m _ { \phi } Y _ { \phi } ( T _ { f } )$$

We show the freezeout of φ in Fig. 2 by plotting the evolution of Y φ with x (= m /T φ ), where T is temperature. Here we take the masses m φ = 40 MeV, m δ = 60 MeV, and coupling λ 11 = 4. The black solid line in the figure represents the value of Y φ while the dotted red line is the Maxwell-Boltzmann equilibrium distribution function.

We analyze the relic abundance of DM to obtain the parameters that satisfy the PLANCK data. The recent PLANCK observation limit on DM relic abundance is given as [6]:

$$\Omega h ^ { 2 } = 0. 1 2 0 6 \pm 0. 0 0 2 1.$$

We analyze the model in the light of the constraints from the Planck data for DM relic density. To have annihilation via mentioned channel, we choose m φ &lt; m δ &lt; 2 m φ and randomly vary the values of λ 11 and λ 12 in their perturbative range of [0 , 4 π ]. In Fig. 3, we show the region of parameter space that satisfy the 1 σ range of Planck data in the M φ -λ 11 plane. The plot reveals that as the DM mass increases, a larger value of λ 11 is required to satisfy the relic density constraint. This is quite expected as the annihilation

Figure 2 : Temperature evolution of number density of φ with v δ = 50 MeV. Black solid line represents the yield Y φ evolution with dimension less variable x . Red dashed line shows the variation of Y eq φ .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00909v1-figure-0ac83946-02cc-4bd6-b0b1-b02dd32865b4.jpg)

Figure 3 : Correlation of coupling and mass of DM that satisfy the relic abundance from PLANCK.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00909v1-figure-0b4b98bd-fdfd-4657-aec8-0ee22133dbc6.jpg)

cross-section decreases with an increase in the DM mass and it increases with an increase in the coupling. Thus, simply increasing the DM mass itself would result in an overabundance of DM density and the value of λ 11 should also increase to satisfy the observed bound.

Next, we study the effect of changing the mass difference between the DM and the singlet scalar field. In Fig. 4 we plot the scalar boson mass with respect to the DM

Figure 4 : Allowed values of scalar mass which satisfy the relic density constraint for DM mass varying from 10 MeV to 100 MeV.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00909v1-figure-38d9450d-c678-4276-8c9d-8c0d26296490.jpg)

mass. For a given dark matter mass m φ , the scalar boson mass M δ is varied between M φ &lt; M δ &lt; 2 M φ . The plot also represents the magnitude of the self-coupling ( λ 11 ) required for each point through a color coding. Similar to the previous plot, an increase in the DM mass must be accompanied by an increase in the self coupling in order to satisfy the relic density constraints. It is also evident from the plot that for a fixed DM mass, the required value of λ 11 increases as M δ increases. This is again quite natural since an increase in the scalar boson mass results in a reduction of the phase space factor as the mass difference between the initial and final states get smaller. Hence a larger value of the coupling would be required to satisfy the observed relic density.

Finally we focus on the effect of varying the VEV of the scalar field on the DM relic density. Fig. 5 shows five different plots corresponding to five different values of v δ and their corresponding allowed parameter points in the M φ -λ 11 plane. Here we have chosen the singlet scalar mass to be degenerate with the DM mass and λ 11 = λ 12 . The purple points corresponding to the v δ = 0 1 MeV can only satisfy the relic density constraints . for large values of the coupling. Since the δφφ coupling is proportional to the v δ , a small value of the VEV requires a large value of the quartic coupling in order to produce enough annihilation of the DM particle. Similarly as the VEV increases one can satisfy the relic density constraints for smaller values of the quartic couplings as well. Another interesting feature of Fig. 5 is that the allowed range of DM masses increases as the value of the VEV increases. A quick glance at the 3 → 2 scattering process reveals that its cross-section is

Figure 5 : Points satisfying the DM relic density for different values of v δ . The purple points correspond to v δ = 0 1 MeV, green for . v δ = 1 MeV, blue for v δ = 10 MeV, orange points represent v δ = 50 MeV and yellow points are for v δ = 100 MeV.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00909v1-figure-eaf919d7-e021-4d82-8d9c-bedda915294b.jpg)

proportional to v 2 δ and hence a larger value of v δ will naturally allow m φ to be large as well.

## 4 Self interaction

Self-interacting DM can solve many of the discrepancies between experimental observations and N-body simulations [14, 15]. As an example, numerical simulations in ΛCDM framework predict the DM density radially diverges at the core of galaxy clusters while actual observations suggest a flat DM distribution [12, 13]. Self-interaction among the DM particles allow them to be scattered by each other and can result in a uniform DM density at the core as is expected from experimental observations. Additionally, observations from galaxy clusters can also put severe constraints on the self-scattering cross-section of DM particles. The mean free path for scattering of DM can be given as:

$$\lambda _ { \text{scatt} } & = \frac { \sigma _ { \text{self} } } { m _ { \phi } } \rho \, v,$$

where σ self , ρ , v and m φ are the self-scattering cross-section, density, velocity and mass of dark matter respectively. Observations from galaxy clusters provide the density at the core and a typical velocity ∼ 50 kms -1 , while the value of λ scatt is around the size of the galaxy cluster. Knowing all these values from astrophysical observations, one can easily provide constraints on the ratio of the self-scattering cross-section to the mass of the DM particle. The constraints on σ self /m φ from different observations are as follows:

Table 2 : Constraints on DM self-scattering cross-section from various observations.

| σ self /m φ (cm 2 / g)     | Observations                                                  |
|----------------------------|---------------------------------------------------------------|
| ∼ (1 . 7 ± 0 . 7) × 10 - 4 | Bright cluster galaxies in the 10 kpc core of Abell 3827 [59] |
| ∼ 0 . 1                    | Cores in clusters [60, 61]                                    |
| ∼ 1 . 5                    | Abell 3827 subhalos [62]                                      |
| ∼ 1                        | Abell 520 cluster [63-65]                                     |
| glyph[lessorsimilar] 1     | Halo shapes and Bullet cluster [66, 67]                       |

Figure 6 : Feynman diagrams for the processes contributing to the self-interaction of dark matter.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00909v1-figure-923aec1b-0f85-4a28-9f22-b7a65f0e6257.jpg)

In our model, the diagrams contributing to self-scattering are presented in the Fig. 6. The amplitude can be written as

$$i \mathcal { M } = i \left ( \lambda _ { 1 1 } + \mu _ { \text{eff} } ^ { 2 } \frac { 1 } { ( s - m _ { \delta } ^ { 2 } ) } + \mu _ { \text{eff} } ^ { 2 } \frac { 1 } { ( t - m _ { \delta } ^ { 2 } ) } + \mu _ { \text{eff} } ^ { 2 } \frac { 1 } { ( u - m _ { \delta } ^ { 2 } ) } \right ),$$

where s , t , u are the Mandelstam variables and µ eff = µ 12 + λ 12 2 v δ . The self-scattering cross-section thus becomes

$$\sigma _ { \text{self} } = \frac { 1 } { 6 4 \pi m _ { \phi } ^ { 2 } } | \mathcal { M } | ^ { 2 }.$$

We compute the self scattering cross-section and analyze the bounds mentioned above. The first term in Eq. 4.2 usually gives the dominant contribution to the self-scattering process and hence we study the allowed range of λ 11 as a function of the DM mass. Table 2 clearly shows that different observations suggest different limits on σ self /m φ so it is not possible to choose a common range satisfying all of them. We have thus chosen an upper bound of σ self /m φ glyph[lessorsimilar] 10 cm 2 / g in order to incorporate all the constraints in a single plot. In the Figure 7 we have shown the allowed parameter space that satisfy the self-scattering bounds. We have divided the plots into various regions, i.e., (i) σ self /m φ glyph[lessorsimilar] 0 1 cm . 2 / g (blue), (ii) 0 1 cm . 2 / g glyph[lessorsimilar] σ self / m φ glyph[lessorsimilar] 1 cm 2 / g (red), (iii) 1 cm 2 / g glyph[lessorsimilar] σ self / m φ glyph[lessorsimilar] 1 5 cm . 2 / g (green), (iv) 1 5 cm . 2 / g glyph[lessorsimilar] σ self / m φ glyph[lessorsimilar] 3 cm 2 / g (yellow) and (v) 3 cm 2 / g glyph[lessorsimilar] σ self / m φ glyph[lessorsimilar] 10 cm 2 / g (cyan). The left panel of Fig. 7 corresponds to v δ = 10 MeV while the right panel is for v δ = 100 MeV. As can be seen from the plots, an increase in the DM mass requires

a larger value of λ 11 to satisfy the σ self /m φ constraints (except for the blue region with σ self /m φ glyph[lessorsimilar] 0 1 cm . 2 / g). This is quite expected as the DM self-scattering cross-section is proportional to its coupling strength. A larger value of m φ requires a larger value of σ self and hence a larger value of λ 11 in order to satisfy the lower bound on σ self /m φ . The absence of this feature in the blue region is because there is no lower limit on σ self /m φ here. This dependence on λ 11 is more prominent for v δ = 10 MeV (left panel) compared to v δ = 100 MeV (right panel). For a low value of v δ , the dominant contribution to the self-scattering cross-section is from the contact interaction (four scalar interaction) diagram as the other process is suppressed by the propagator of δ . As we increase v δ , the value of µ eff increases and the second term starts contributing and hence the dependence on λ 11 is flattened out a bit. This is clearly visible on the right panel of Fig. 7 where v δ = 100 MeV. Another distinction between the two plots is that the allowed range of λ 11 for each region is much larger for the left panel plot compared to the right panel. As the second term in Eq. 4.2 contributes significantly for the v δ = 100 MeV case, a smaller value of λ 11 is required here. On the other hand, for v δ = 10 MeV case, the entire self-interaction cross-section is due to the first term which is directly proportional to λ 2 11 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00909v1-figure-2b0e50ae-e77b-4362-a118-9128ad7d29f5.jpg)

0.02

0.04

0.06

0.08

0.10

m

ϕ

(

GeV

)

Figure 7 : Regions showing the parameter space that satisfy σ self /m φ glyph[lessorsimilar] 10 cm 2 / g. Color coding for various divisions in σ self /m φ limit: σ self /m φ glyph[lessorsimilar] 0 1 . cm 2 / g (blue), 0 1 . cm 2 / g glyph[lessorsimilar] σ self / m φ glyph[lessorsimilar] 1 cm 2 / g (red), 1 cm 2 / g glyph[lessorsimilar] σ self / m φ glyph[lessorsimilar] 1 5 . cm 2 / g (green), 1 5 cm . 2 / g glyph[lessorsimilar] σ self / m φ glyph[lessorsimilar] 3 cm 2 / g (yellow) and 3 cm 2 / g glyph[lessorsimilar] σ self / m φ glyph[lessorsimilar] 10 cm 2 / g (cyan).

In Fig. 8 we consider the σ self /m φ constraints arising from the recently observed Bullet and Abell galaxy clusters. The Bullet cluster bound suggests σ self /m φ glyph[lessorsimilar] 1 cm 2 / g [66] while Abell cluster suggest 1 cm 2 / g glyph[lessorsimilar] σ self / m φ glyph[lessorsimilar] 3 cm 2 / g [68]. The left panel in Fig. 8 represents the Bullet cluster bound and the right panel is for the Abell cluster bound for v δ = 10 MeV. As the two observations do not agree with each other, their preferred regions of allowed parameters are also completely disjointed in our model. It is thus not possible for us to choose an allowed range of parameters based on these two observations.

12

10

8

6

4

2

0

11

λ

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00909v1-figure-26b9866f-1d7f-413a-a87a-5054e7b805bd.jpg)

0.02

0.04

0.06

0.08

0.10

m

ϕ

(

GeV

)

Figure 8 : Shaded regions in λ 11 and m φ plane that satisfy the Bullet cluster (left) and Abell cluster data (right). VEV of δ is assumed to be 10 MeV.

## 5 Neutrino mass generation

The inclusion of three right-handed neutrinos ( N i ) in the model helps generate the light neutrino masses through type-I seesaw mechanism. The neutrino specific terms in the Lagrangian are given as

$$\mathcal { L } _ { Y } \supset [ Y _ { D _ { i j } } L _ { i } ^ { T } i \sigma _ { 2 } H N _ { j } + H. C. ] + M _ { N _ { i j } } N _ { i } ^ { c } N _ { j } + f _ { i j } N _ { i } ^ { c } N _ { j } \delta.$$

The neutrino mass matrix is thus a 6 × 6 matrix which can be written as

$$\begin{pmatrix} 0 & M _ { D } \\ M _ { D } ^ { T } & M _ { R } \end{pmatrix},$$

where M D ij = Y D ij v H is the Dirac mass term while M R ij = M N ij + f ij v δ is the Majorana mass term for the right-handed neutrino. Without loss of generality we can choose one of the right-handed neutrinos to be much lighter with mass around a few MeV. This will help in generating the required baryon asymmetry through leptogenesis as will be discussed in the next section.

The experimentally observed 3 σ ranges [69, 70] for the neutrino mass-squared differences and the corresponding Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix elements for normal hierarchy case are given in Table 3. In our scenario, we choose the two heavy right-handed neutrinos to be almost degenerate with a mass of around 10 TeV while the lighter right-handed neutrino mass is chosen to be around 10 MeV. For this mass spectrum of the singlet neutrinos, we perform a scan over the neutrino Dirac masses ( M D ij ) which can satisfy the experimental constraints. Fig. 9 shows a scatter plot of the Dirac masses which satisfy the neutrino mass-squared differences and mixings as given in Table 3.

It can be seen from Fig. 9 that there is a hierarchy in the neutrino Dirac masses. As we have chosen the lightest singlet neutrino ( N 1 ) to be only around 10 MeV, the mass terms M D 11 , M D 21 and M D 31 are much smaller than the other terms in order to satisfy the seesaw structure of the mass generation mechanism. As is clearly visible in the left panel

11

λ

12

10

8

6

4

2

0

Table 3 : Experimental 3 σ ranges for light neutrino parameters from NuFIT [69, 70].

| 6.79 × 10 - 5 eV 2 < ∆ m 2 21 < 8.01 × 10 - 5 eV 2   | 6.79 × 10 - 5 eV 2 < ∆ m 2 21 < 8.01 × 10 - 5 eV 2                                                                                  |
|------------------------------------------------------|-------------------------------------------------------------------------------------------------------------------------------------|
| 2.432 × 10 - 3 eV 2 < ∆ m 2 31 < 2.618 × 10 - 3 eV 2 | 2.432 × 10 - 3 eV 2 < ∆ m 2 31 < 2.618 × 10 - 3 eV 2                                                                                |
| 0 . 275 < sin 2 θ 12 < 0 . 350                       | 0 . 275 < sin 2 θ 12 < 0 . 350                                                                                                      |
| 0 . 427 < sin 2 θ 23 < 0 . 609                       | 0 . 427 < sin 2 θ 23 < 0 . 609                                                                                                      |
| 0 . 02046 < sin 2 θ 13 < 0 . 02440                   | 0 . 02046 < sin 2 θ 13 < 0 . 02440                                                                                                  |
| U PMNS                                               |       0 . 797 → 0 . 842 0 . 518 → 0 . 585 0 . 143 → 0 . 156 0 . 244 → 0 . 496 0 . 467 → 0 . 678 0 . 646 → 0 . 772       |
| U PMNS                                               | 0 . 287 → 0 . 525 0 . 488 → 0 . 693 0 . 618 → 0 . 749                                                                               |

Figure 9 : Neutrino Dirac masses satisfying the experimental mass-squared differences and mixings with heavy neutrino masses M N 1 = 10 MeV, M N 2 ≈ 10 TeV, M N 3 ≈ 10 TeV.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00909v1-figure-ea779560-380e-4784-a1bb-d2b027d815f9.jpg)

of the figure, these three small masses appear almost degenerate compared to the other terms. To show the hierarchy between them, we have included the right panel which only shows these three smaller masses. For both the figures, we have plotted the Dirac masses with respect to the M D 11 as it shows how the other mass terms are varying in magnitude as M D 11 increases.

The left panel of Fig. 9 gives a clear indication that M D 33 is the largest in each case. This can be understood as an effect of our choice of normal hierarchy for the neutrino mass spectra. Since ν 3 is the heaviest, it implicitly requires the Dirac mass term M D 33 to be the largest in each case. The plot also shows that if we neglect the much smaller M D 11 , M D 21 and M D 31 mass terms, the value of M D 13 is the smallest among the rest. A similar feature can be observed in the right panel of Fig. 9 where M D 31 is the smallest. This is quite expected as the two terms M D 13 and M D 31 are responsible for generation of the θ 13 mixing between ν 1 and ν 3 . The smallness of the experimentally measured value of θ 13 thus

requires these two Dirac mass terms to be quite small as well.

Another interesting feature in this figure is that the magnitude of M D 12 and M D 21 increases as M D 11 increases. Since the M D 12 and M D 21 are responsible for generation of the mixing between ν 1 and ν 2 , it is quite obvious that these term have to increase as the diagonal terms increase. This feature can also be seen in M D 13 and M D 31 mass terms though it is not as pronounced as the other case, maybe due to the smallness of θ 13 .

## 6 Leptogenesis

The observed baryon asymmetry of the universe can be explained in this model through leptogenesis. The CP asymmetry is generated through the interference of the tree-level and one-loop diagrams for the heavy neutrino decay. Let us consider the decay of one of the heavy right-handed neutrinos into a final state lepton and a scalar boson. There are two possible decays here. The heavy neutrino ( N 2 or N 3 ) can decay into an active neutrino and Higgs boson or it can decay into the lightest right-handed neutrino ( N 1 ) and singlet scalar boson ( δ ). Assuming the final state leptons and scalar boson masses are much smaller

Figure 10 : Diagrams contributing to leptogenesis. The required CP asymmetry is generated through the interference of the tree-level process with the one-loop diagrams.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00909v1-figure-b5094954-a3e0-4e99-a4de-9190d7a2db1a.jpg)

compared to the right-handed neutrino, the required lepton number asymmetry will be given as

$$\epsilon _ { 1 } = - \sum _ { i } \left [ \frac { \Gamma ( N _ { 2 } \to \bar { l } _ { i } H ^ { * } ) - \Gamma ( N _ { 2 } \to l _ { i } H ) } { \Gamma _ { t o t } ( N _ { 2 } ) } \right ].$$

Here N 2 represents the heavy right-handed neutrino and Γ tot ( N 2 ) is its total decay width given as

$$\Gamma _ { t o t } ( N _ { 2 } ) = \frac { ( Y _ { D _ { 2 i } } ^ { \dagger } Y _ { D _ { 2 i } } ) + | f _ { 1 2 } | ^ { 2 } } { 4 \pi } M _ { N _ { 2 } }.$$

The Y D term corresponds to the decay into a SM Higgs boson and a light neutrino while the f coupling is due to the decay into a light singlet neutrino and singlet scalar boson. The one-loop diagrams for these decay processes will involve the generic self-energy and vertex corrections due to the Y D term in the Yukawa Lagrangian. In addition, there will be an extra self-energy correction due to the fN N δ c 1 j term as well. Fig. 10 shows the tree level and one-loop Feynman diagrams contributing to the process. As a result, the CP asymmetry in this case will be given as

$$\epsilon _ { 1 } = \frac { 1 } { 8 \pi } \left ( [ g _ { V } ( x ) + g _ { S } ( x ) ] \mathcal { T } _ { 2 3 } + g _ { S } ( x ) \mathcal { S } _ { 2 3 } \right ),$$

where g V ( x ) = √ x { 1 -(1 + x )ln[(1 + x /x ) ] } , g S ( x ) = √ x/ (1 -x ) with x = M 2 N 3 /M 2 N 2 ,

$$\mathcal { T } _ { 2 3 } = \frac { \text{Im} [ ( Y _ { D _ { 2 i } } Y _ { D _ { 3 i } } ^ { \dagger } ) ^ { 2 } ] } { ( Y _ { D _ { 2 i } } ^ { \dagger } Y _ { D _ { 2 i } } ) + | f _ { 2 1 } | ^ { 2 } }$$

and

$$\mathcal { S } _ { k 1 } = \frac { \text{Im} [ ( Y _ { D _ { 2 i } } Y _ { D _ { 3 i } } ^ { \dagger } ) ( f _ { 2 1 } f _ { 3 1 } ^ { \dagger } ) ] } { ( Y _ { D _ { 2 i } } ^ { \dagger } Y _ { D _ { 2 i } } ) + | f _ { 2 1 } | ^ { 2 } }.$$

Here the g V ( x ) term arrives from the vertex correction while the g S ( x ) terms are due to the self energy corrections. The S k 1 term is purely due to the new contribution from the self energy diagram involving the singlet scalar field as discussed before.

It is quite interesting to see that when M N 2 glyph[similarequal] M N 3 (i.e. x glyph[similarequal] 1), the correction from the self energy term is dominant and can significantly enhance the asymmetry factor given in Eq. 6.3. The lepton asymmetry in this case can be approximately given by

$$\epsilon _ { 1 } \simeq - \frac { 1 } { 1 6 \pi } \left [ \frac { M _ { N _ { 3 } } } { v ^ { 2 } } \frac { I m [ ( Y _ { D } ^ { * } m _ { \nu } Y _ { D } ^ { \dagger } ) _ { 2 2 } ] } { ( Y _ { D } ^ { \dagger } Y _ { D } ) _ { 2 2 } + | f _ { 2 1 } | ^ { 2 } } + \frac { I m [ ( Y _ { D } Y _ { D } ^ { \dagger } ) _ { 2 3 } ( f _ { 2 1 } f _ { 3 1 } ^ { \dagger } ) ] } { ( Y _ { D } ^ { \dagger } Y _ { D } ) _ { 2 2 } + | f _ { 2 1 } | ^ { 2 } } \right ] R$$

where R may be recognized as a resonance factor defined as R ≡ | M N 2 | / ( | M N 3 |-| M N 2 | ). In the absence of the f ij couplings, only the first term of Eq. 6.6 will contribute and a value of R ∼ 10 6 -7 is required to get the correct amount of lepton number asymmetry. This leads to almost degenerate masses for the heavy right-handed neutrinos with M N 3 -M N 2 ≈ 10 -2 GeV resulting in an unnaturally large fine-tuned scenario. The presence of the extra couplings between the singlet scalar and singlet neutrinos can alleviate this problem by enhancing the second term in Eq. 6.6. The f 21 coupling though is constrained from the the out-of-equilibrium condition which requires the total decay width of N 2 to be less than the Hubble constant. Thus we get

$$\Gamma _ { N _ { 2 } } < H | _ { T = M _ { N _ { 2 } } }$$

where H is the Hubble expansion rate. This produces an upper bound on f 21 with

$$\sqrt { | f _ { 2 1 } | ^ { 2 } } < 3 \times 1 0 ^ { - 4 } \sqrt { M _ { N _ { 2 } } / 1 0 ^ { 9 } ( \text{GeV} ) }.$$

This constraint however is not applicable for the f 31 coupling and it can be taken to be much larger. The B -L asymmetry generated in this scenario is given as

$$Y _ { B - L } = - \eta \epsilon _ { 1 } Y _ { N _ { 2 } } ^ { e q },$$

where η is the efficiency factor and Y eq N 2 glyph[similarequal] 45 π 4 ζ (3) g ∗ k B 3 4 is the number density of N 2 in thermal equilibrium at T &gt; &gt;M N 1 . Here k B is the Boltzmann constant and g ∗ is the effective number of degrees of freedom. This asymmetry can then be transferred into a baryon asymmetry through sphaleron induced processes at the electroweak symmetry breaking scale.

If one chooses f 21 ∼ Y D 21 with f 31 ∼ χY D 31 , then the second term in Eq. 6.6 increases by a factor of χ and the required asymmetry can be obtained even with a much smaller resonance factor R . If we choose χ ∼ 10 , the value of 3 R required is only around 10 3 -4 and even a mass difference of a few GeV between the two heavy neutrinos could be allowed in order to generate the B -L asymmetry. Thus we arrive at a scenario which is free from the extreme fine-tuning between the heavy neutrino masses and is instead realized through a simple hierarchy in the f ij couplings which is much more natural.

## 7 Neutrinoless Double Beta Decay

Figure 11 : Feynman diagram for dominant contribution to the neutrinoless double beta decay process mediated by the light right-handed neutrino in this model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00909v1-figure-f9a0b4b1-5b9f-4759-9bb5-eb1ebf7aec41.jpg)

The low mass of N 1 can results in significant bound from neutrinoless double beta decay (0 ν 2 β ) processes [71-73]. The Feynman diagram in Fig. 11 shows the dominant contribution to 0 ν 2 β signal in this model. As is quite evident from the diagram, the 0 ν 2 β decay

process will be directly proportional to the mass of the light right-handed neutrino and its mixing with the electron neutrino. This process can thus produce significant constraints on the allowed mixing between the active and sterile neutrinos. A detailed calculation of this diagram and its effect on the neutrino spectrum will be discussed in an upcoming paper.

## 8 Conclusion

In this paper we have constructed a unified model of a singlet scalar self-interacting DM with a mass of ∼ 10 MeV. We have given a 3 → 2 annihilation mechanism for generating the required relic density of the DM with a second singlet scalar in the final state. This second singlet acquires a non-zero VEV which helps in generating the DM mass. Solving the coupled Boltzmann equations, we identify the allowed region of parameter space that satisfy the experimentally observed bounds on DM relic abundance. Our results show that a larger DM mass requires a larger value of the self-interaction coupling λ 11 to satisfy the DM relic density. This is because the 3 → 2 annihilation process is mediated by the DM particle and hence is inversely proportional to its mass. One thus needs to increase the DM self-interaction coupling. A similar result can be obtained by simply increasing the scalar VEV while keeping λ 11 constant as this also helps in enhancing the DM annihilation cross-section. Self-interacting DM can have stringent bounds from several astrophysical observations. We analyze the model for bounds from many such observations focusing specifically on the recently observed Bullet and Abell galaxy clusters. Since these two bounds are mutually exclusive, we have presented our allowed parameter regions which can satisfy the constraints from either of them.

We have also included three right-handed neutrinos which help generate the light neutrino masses via Type-I seesaw mechanism. We choose two of the right-handed neutrinos to be heavy around 10 TeV while one of them remains light at around 10 MeV. Detailed scans of the neutrino parameters were performed to identify a distinct pattern in their Dirac mass terms which is required to satisfy the experimental observations. This model can also explain the observed baryon asymmetry of the universe. The interaction term of the right-handed neutrinos with the scalar singlet can help generate a large lepton asymmetry which can then be converted into baryon asymmetry through sphaleron transitions. One of the possible way to test this model is through neutrino-less double beta decay as the light right-handed neutrino can have significant contribution to the process. The absence of any such signals will put stringent bounds on the model.

## References

- [1] E. W. Kolb and M. S. Turner, The Early Universe , Front. Phys. 69 (1990) 1-547.
- [2] L. Roszkowski, E. M. Sessolo, and S. Trojanowski, WIMP dark matter candidates and searchescurrent status and future prospects , Rept. Prog. Phys. 81 (2018), no. 6 066201, [ arXiv:1707.06277 ].
- [3] S. C. Kim et al., New Limits on Interactions between Weakly Interacting Massive Particles and Nucleons Obtained with CsI(Tl) Crystal Detectors , Phys. Rev. Lett. 108 (2012) 181301, [ arXiv:1204.2646 ].

- [4] XENON Collaboration, E. Aprile et al., Physics reach of the XENON1T dark matter experiment , JCAP 1604 (2016), no. 04 027, [ arXiv:1512.07501 ].
- [5] XENON Collaboration, E. Aprile et al., Search for WIMP Inelastic Scattering off Xenon Nuclei with XENON100 , Phys. Rev. D96 (2017), no. 2 022008, [ arXiv:1705.05830 ].
- [6] Planck Collaboration, N. Aghanim et al., Planck 2018 results. VI. Cosmological parameters , arXiv:1807.06209 .
- [7] A. Refregier, Weak gravitational lensing by large scale structure , Ann. Rev. Astron. Astrophys. 41 (2003) 645-668, [ astro-ph/0307212 ].
- [8] B. Moore, S. Ghigna, F. Governato, G. Lake, T. R. Quinn, J. Stadel, and P. Tozzi, Dark matter substructure within galactic halos , Astrophys. J. 524 (1999) L19-L22, [ astro-ph/9907411 ].
- [9] A. A. Klypin, A. V. Kravtsov, O. Valenzuela, and F. Prada, Where are the missing Galactic satellites? , Astrophys. J. 522 (1999) 82-92, [ astro-ph/9901240 ].
- [10] M. Boylan-Kolchin, J. S. Bullock, and M. Kaplinghat, Too big to fail? The puzzling darkness of massive Milky Way subhaloes , Monthly Notices of the Royal Astronomical Society: Letters 415 (07, 2011) L40-L44,
- [ http://oup.prod.sis.lan/mnrasl/article-pdf/415/1/L40/2881682/415-1-L40.pdf ].
- [11] M. Boylan-Kolchin, J. S. Bullock, and M. Kaplinghat, The Milky Ways bright satellites as an apparent failure of CDM , Monthly Notices of the Royal Astronomical Society 422 (04, 2012) 1203-1218,
- [ http://oup.prod.sis.lan/mnras/article-pdf/422/2/1203/3464467/mnras0422-1203.pdf ].
- [12] A. Burkert, The Structure of dark matter halos in dwarf galaxies , IAU Symp. 171 (1996) 175, [ astro-ph/9504041 ]. [Astrophys. J.447,L25(1995)].
- [13] P. Salucci, A. Lapi, C. Tonini, G. Gentile, I. Yegorova, and U. Klein, The Universal Rotation Curve of Spiral Galaxies. 2. The Dark Matter Distribution out to the Virial Radius , Mon. Not. Roy. Astron. Soc. 378 (2007) 41-47, [ astro-ph/0703115 ].
- [14] D. N. Spergel and P. J. Steinhardt, Observational evidence for self-interacting cold dark matter , Phys. Rev. Lett. 84 (Apr, 2000) 3760-3763.
- [15] S. Balberg, S. L. Shapiro, and S. Inagaki, Selfinteracting dark matter halos and the gravothermal catastrophe , Astrophys. J. 568 (2002) 475-487, [ astro-ph/0110561 ].
- [16] S. Tulin and H.-B. Yu, Dark Matter Self-interactions and Small Scale Structure , Phys. Rept. 730 (2018) 1-57, [ arXiv:1705.02358 ].
- [17] 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), no. 2 021301, [ arXiv:1411.3727 ].
- [18] Y. Hochberg, E. Kuflik, and H. Murayama, SIMP Spectroscopy , JHEP 05 (2016) 090, [ arXiv:1512.07917 ].
- [19] N. Bernal, C. Garcia-Cely, and R. Rosenfeld, WIMP and SIMP Dark Matter from the Spontaneous Breaking of a Global Group , JCAP 1504 (2015), no. 04 012, [ arXiv:1501.01973 ].
- [20] S.-M. Choi and H. M. Lee, SIMP dark matter with gauged Z 3 symmetry , JHEP 09 (2015) 063, [ arXiv:1505.00960 ].

- [21] U. K. Dey, T. N. Maity, and T. S. Ray, Light Dark Matter through Assisted Annihilation , JCAP 1703 (2017), no. 03 045, [ arXiv:1612.09074 ].
- [22] J. M. Cline, H. Liu, T. Slatyer, and W. Xue, Enabling Forbidden Dark Matter , Phys. Rev. D96 (2017), no. 8 083521, [ arXiv:1702.07716 ].
- [23] A. Kamada, H. Kim, and T. Sekiguchi, Axionlike particle assisted strongly interacting massive particle , Phys. Rev. D96 (2017), no. 1 016007, [ arXiv:1704.04505 ].
- [24] S.-M. Choi, Y. Hochberg, E. Kuflik, H. M. Lee, Y. Mambrini, H. Murayama, and M. Pierre, Vector SIMP dark matter , JHEP 10 (2017) 162, [ arXiv:1707.01434 ].
- [25] S.-Y. Ho, T. Toma, and K. Tsumura, A Radiative Neutrino Mass Model with SIMP Dark Matter , JHEP 07 (2017) 101, [ arXiv:1705.00592 ].
- [26] M. Heikinheimo, K. Tuominen, and K. Langble, Hidden strongly interacting massive particles , Phys. Rev. D97 (2018), no. 9 095040, [ arXiv:1803.07518 ].
- [27] Y. Hochberg, E. Kuflik, and H. Murayama, Twin Higgs model with strongly interacting massive particle dark matter , Phys. Rev. D99 (2019), no. 1 015005, [ arXiv:1805.09345 ].
- [28] Y. Hochberg, E. Kuflik, R. Mcgehee, H. Murayama, and K. Schutz, Strongly interacting massive particles through the axion portal , Phys. Rev. D98 (2018), no. 11 115031, [ arXiv:1806.10139 ].
- [29] M. H. Namjoo, T. R. Slatyer, and C.-L. Wu, Enhanced n-body annihilation of dark matter and its indirect signatures , JHEP 03 (2019) 077, [ arXiv:1810.09455 ].
- [30] U. K. Dey, T. N. Maity, and T. S. Ray, Boosting Assisted Annihilation for a Cosmologically Safe MeV Scale Dark Matter , Phys. Rev. D99 (2019), no. 9 095025, [ arXiv:1812.11418 ].
- [31] S.-M. Choi, H. M. Lee, Y. Mambrini, and M. Pierre, Vector SIMP dark matter with approximate custodial symmetry , arXiv:1904.04109 .
- [32] S. Bhattacharya, P. Ghosh, and S. Verma, SIMPler realisation of Scalar Dark Matter , arXiv:1904.07562 .
- [33] N. Bernal and X. Chu, Z 2 SIMP Dark Matter , JCAP 1601 (2016) 006, [ arXiv:1510.08527 ].
- [34] S.-M. Choi and H. M. Lee, Resonant SIMP dark matter , Phys. Lett. B758 (2016) 47-53, [ arXiv:1601.03566 ].
- [35] S.-M. Choi, Y.-J. Kang, and H. M. Lee, On thermal production of self-interacting dark matter , JHEP 12 (2016) 099, [ arXiv:1610.04748 ].
- [36] N. Bernal, X. Chu, and J. Pradler, Simply split strongly interacting massive particles , Phys. Rev. D95 (2017), no. 11 115023, [ arXiv:1702.04906 ].
- [37] B. Chauhan, Sub-MeV Self Interacting Dark Matter , Phys. Rev. D97 (2018), no. 12 123017, [ arXiv:1711.02970 ].
- [38] M. Fukugita and T. Yanagida, Baryogenesis Without Grand Unification , Phys. Lett. B174 (1986) 45-47.
- [39] A. Pilaftsis, CP violation and baryogenesis due to heavy Majorana neutrinos , Phys. Rev. D56 (1997) 5431-5451, [ hep-ph/9707235 ].
- [40] A. Pilaftsis and T. E. J. Underwood, Resonant leptogenesis , Nucl. Phys. B692 (2004) 303-345, [ hep-ph/0309342 ].

- [41] S. K. Kang and C. S. Kim, Extended double seesaw model for neutrino mass spectrum and low scale leptogenesis , Phys. Lett. B646 (2007) 248-252, [ hep-ph/0607072 ].
- [42] S. K. Kang and A. Patra, keV Sterile Neutrino Dark Matter and Low Scale Leptogenesis , J. Korean Phys. Soc. 69 (2016), no. 8 1375-1382, [ arXiv:1412.4899 ].
- [43] P. Minkowski, µ → eγ at a Rate of One Out of 10 9 Muon Decays? , Phys. Lett. 67B (1977) 421-428.
- [44] T. Yanagida, HORIZONTAL SYMMETRY AND MASSES OF NEUTRINOS Conf. Proc. , C7902131 (1979) 95-99.
- [45] O. Sawada and A. Sugamoto, eds., Proceedings: Workshop on the Unified Theories and the Baryon Number in the Universe , (Tsukuba, Japan), Natl.Lab.High Energy Phys., Natl.Lab.High Energy Phys., 1979.
- [46] M. Levy, J. L. Basdevant, D. Speiser, J. Weyers, R. Gastmans, and M. Jacob, eds., QUARKS AND LEPTONS. PROCEEDINGS, SUMMER INSTITUTE, CARGESE, FRANCE, JULY 9-29, 1979 , vol. 61, 1980.
- [47] P. Van Nieuwenhuizen and D. Z. Freedman, eds., SUPERGRAVITY. PROCEEDINGS, WORKSHOP AT STONY BROOK, 27-29 SEPTEMBER 1979 , 1979.
- [48] R. N. Mohapatra and G. Senjanovic, Neutrino Mass and Spontaneous Parity Violation , Phys. Rev. Lett. 44 (1980) 912.
- [49] CMS Collaboration, A. M. Sirunyan et al., Search for invisible decays of a Higgs boson produced through vector boson fusion in proton-proton collisions at √ s = 13 TeV , Phys. Lett. B793 (2019) 520-551, [ arXiv:1809.05937 ].
- [50] S. Kang, S. Scopel, G. Tomar, and J.-H. Yoon, Present and projected sensitivities of Dark Matter direct detection experiments to effective WIMP-nucleus couplings , Astropart. Phys. 109 (2019) 50-68, [ arXiv:1805.06113 ].
- [51] CRESST Collaboration, F. Petricca et al., First results on low-mass dark matter from the CRESST-III experiment , in 15th International Conference on Topics in Astroparticle and Underground Physics (TAUP 2017) Sudbury, Ontario, Canada, July 24-28, 2017 , 2017. arXiv:1711.07692 .
- [52] XENON Collaboration, E. Aprile et al., Dark Matter Search Results from a One Ton-Year Exposure of XENON1T , Phys. Rev. Lett. 121 (2018), no. 11 111302, [ arXiv:1805.12562 ].
- [53] CDEX Collaboration, Z. Z. Liu et al., Constraints on spin-independent nucleus scattering with sub-GeV WIMP dark matter from the CDEX-1B Experiment at CJPL , arXiv:1905.00354 .
- [54] R. Essig, T. Volansky, and T.-T. Yu, New Constraints and Prospects for sub-GeV Dark Matter Scattering off Electrons in Xenon , Phys. Rev. D96 (2017), no. 4 043017, [ arXiv:1703.00910 ].
- [55] SuperCDMS Collaboration, R. Agnese et al., Search for Low-Mass Dark Matter with CDMSlite Using a Profile Likelihood Fit , Phys. Rev. D99 (2019), no. 6 062001, [ arXiv:1808.09098 ].
- [56] SENSEI Collaboration, M. Crisler, R. Essig, J. Estrada, G. Fernandez, J. Tiffenberg, M. Sofo haro, T. Volansky, and T.-T. Yu, SENSEI: First Direct-Detection Constraints on sub-GeV Dark Matter from a Surface Run , Phys. Rev. Lett. 121 (2018), no. 6 061803, [ arXiv:1804.00088 ].

- [57] K.-C. Yang, Thermodynamic Evolution of Secluded Vector Dark Matter: Conventional WIMPs and Nonconventional WIMPs , arXiv:1905.09582 .
- [58] A. Berlin, D. Hooper, and G. Krnjaic, Thermal Dark Matter From A Highly Decoupled Sector , Phys. Rev. D94 (2016), no. 9 095019, [ arXiv:1609.02555 ].
- [59] R. Massey et al., The behaviour of dark matter associated with four bright cluster galaxies in the 10 kpc core of Abell 3827 , Mon. Not. Roy. Astron. Soc. 449 (2015), no. 4 3393-3406, [ arXiv:1504.03388 ].
- [60] O. D. Elbert, J. S. Bullock, M. Kaplinghat, S. Garrison-Kimmel, A. S. Graus, and M. Rocha, A testable conspiracy: Simulating baryonic effects on self-interacting dark matter halos , The Astrophysical Journal 853 (jan, 2018) 109.
- [61] M. Kaplinghat, S. Tulin, and H.-B. Yu, Dark matter halos as particle colliders: Unified solution to small-scale structure puzzles from dwarfs to clusters , Phys. Rev. Lett. 116 (Jan, 2016) 041302.
- [62] F. Kahlhoefer, K. Schmidt-Hoberg, J. Kummer, and S. Sarkar, On the interpretation of dark matter self-interactions in Abell 3827 , Mon. Not. Roy. Astron. Soc. 452 (2015), no. 1 L54-L58, [ arXiv:1504.06576 ].
- [63] Spethmann, Christian, Veerm¨e, Hardi, Sepp, Tiit, Heikinheimo, Matti, Deshev, Boris, a Hektor, Andi, and Raidal, Martti, Simulations of galaxy cluster collisions with a dark plasma component , A&amp;A 608 (2017) A125.
- [64] M. J. Jee, H. Hoekstra, A. Mahdavi, and A. Babul, HUBBLE SPACE TELESCOPE/ADVANCED CAMERA FOR SURVEYS CONFIRMATION OF THE DARK SUBSTRUCTURE IN a520 , The Astrophysical Journal 783 (feb, 2014) 78.
- [65] F. Kahlhoefer, K. Schmidt-Hoberg, M. T. Frandsen, and S. Sarkar, Colliding clusters and dark matter self-interactions , Mon. Not. Roy. Astron. Soc. 437 (2014), no. 3 2865-2881, [ arXiv:1308.3419 ].
- [66] 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, [ arXiv:0704.0261 ].
- [67] A. H. G. Peter, M. Rocha, J. S. Bullock, and M. Kaplinghat, Cosmological Simulations with Self-Interacting Dark Matter II: Halo Shapes vs. Observations , Mon. Not. Roy. Astron. Soc. 430 (2013) 105, [ arXiv:1208.3026 ].
- [68] F. Kahlhoefer, K. Schmidt-Hoberg, J. Kummer, and S. Sarkar, On the interpretation of dark matter self-interactions in Abell 3827 , Monthly Notices of the Royal Astronomical Society: Letters 452 (07, 2015) L54-L58,

[ http://oup.prod.sis.lan/mnrasl/article-pdf/452/1/L54/3241655/slv088.pdf ].

- [69] I. Esteban, M. C. Gonzalez-Garcia, A. Hernandez-Cabezudo, M. Maltoni, and T. Schwetz, Global analysis of three-flavour neutrino oscillations: synergies and tensions in the determination of θ 2 3 , δ C P , and the mass ordering , JHEP 01 (2019) 106, [ arXiv:1811.05487 ].
- [70] 'Nufit 4.1 (2019).' www.nu-fit.org.
- [71] W. Rodejohann, Neutrino-less Double Beta Decay and Particle Physics , Int. J. Mod. Phys. E20 (2011) 1833-1930, [ arXiv:1106.1334 ].

- [72] F. F. Deppisch, M. Hirsch, and H. Pas, Neutrinoless Double Beta Decay and Physics Beyond the Standard Model , J. Phys. G39 (2012) 124007, [ arXiv:1208.0727 ].
- [73] L. Graf, F. F. Deppisch, F. Iachello, and J. Kotila, Short-Range Neutrinoless Double Beta Decay Mechanisms , Phys. Rev. D98 (2018), no. 9 095023, [ arXiv:1806.06058 ].
Subhendra Mohanty Ayon Patra Tripurari Srivastava
MeV scale model of SIMP dark matter, neutrino mass and leptogenesis
We consider a simple extension of the Standard Model with two singlet scalar
fields and three heavy right-handed neutrinos. One of the scalar fields serves
as an MeV scale dark matter and its stability is ensured by the introduction of
an extra $Z_2$ symmetry. The second scalar (which is even under the $Z_2$
symmetry) generates the mass term of the scalar, contributes to the $3
\rightarrow 2$ annihilation process required for the correct relic density of
the dark matter and it also contributes to the leptogenesis. The right-handed
neutrinos are responsible for the generation of light neutrino masses through
Type-I seesaw mechanism. The decay of the heavy right-handed neutrino can
generate the lepton asymmetry which can then be converted to baryon asymmetry
through sphaleron transitions.