Dark matter, neutrino mass, cutoff for cosmic-ray neutrino, and the Higgs boson invisible decay from a neutrino portal interaction

We study an effective theory beyond the standard model (SM) where either of the two additional gauge singlets, a Majorana fermion and a real scalar, constitutes all or some fraction of dark matter. In particular, we focus on the masses of the two singlets in the range of (10) MeV- (10) GeV with a neutrino portal interaction, which plays an important role not only in particle physics but also in cosmology and astronomy. We point out that the thermal dark matter abundance can be explained by (co-)annihilation, where the dark matter with a mass greater than 2 GeV can be tested in future lepton colliders, CEPC, ILC, FCC-ee and CLIC, in the light of the Higgs boson invisible decay. When the gauge singlets are lighter than (100) MeV, the interaction can affect the neutrino propagation in the universe due to its annihilation with cosmic background neutrino into the gauge singlets. Although in this case it can not be the dominant dark matter, the singlets are produced by the invisible decay of the Higgs boson at such a rate which is fully within reach of future lepton colliders. In particular, a high energy cutoff of cosmic-ray neutrino, which may account for the non-detection of Greisen-Zatsepin-Kuzmin (GZK) neutrino or the non-observation of the Glashow resonance, can be set. Interestingly, given the cutoff and the mass (range) of WIMPs, a neutrino mass can be " measured” kinematically.


Introduction
Weakly Interacting Massive Particles (WIMPs) are promising candidates for dark matter [1][2][3][4]. However, WIMPs with mass 6 GeV-(10 2 ) TeV have been severely constrained by the XENON, LUX and PandaX experiments [5][6][7][8][9][10][11][12]. This situation gives the motivation to investigate WIMPs lighter than a GeV. Such a WIMP should be a singlet of the Standard Model (SM) gauge group to avoid the LEP constraints [13]. If a gauge singlet dark matter is stabilized by a hidden symmetry, its possible interaction with the SM particles is represented by a portal coupling , where ( ) is an SM gauge singlet operator composed only of SM fields (only of hidden fields including WIMP).
It is interesting to study the neutrino portal interaction, i.e.
, where L is a Weyl spinor for a left-ϕ H handed lepton, and we will take Weyl representation hereafter; is the Higgs doublet field; the dot denotes the contraction of the SU(2) gauge indices, while the Lorentz indices are omitted. This is because this interaction can be not only a window from the SM to a dark sector but also affects neutrino and the Higgs boson physics.
Neutrino portal dark matter has been studied in several contexts: asymmetric dark matter [14,15], decaying dark matter [16], and WIMP dark matter [17][18][19][20]. The first part of this paper can be classified as the last one. In particular, we will focus on the dark matter mass range between ~ 10 MeV and ~ 10 GeV, which differs from the previous studies where the mass is greater than several GeV. More concretely, we take an effective field theory approach based on the strategy of simplicity, and focus on the simplest neutrino portal operator of dimension five, where ( ) is a Majorana fermion (a real scalar) carrying a hidden charge, and is a dimensionless coupling. Therefore, the lighter of and is stable.
We point out that and are restricted to be nearly degenerate to satisfy the neutrino mass constraints from the observations of cosmic microwave background (CMB) and baryon acoustic oscillations (BAO) [21]. Otherwise, a sizable neutrino mass would be produced radiatively from the neutrino portal interaction. This allows the Higgs boson to decay into and neutrino kinematically. Such an invisible decay rate could be measured in several future lepton colliders, such as the Circular Electron Positron Collider (CEPC), International Linear Collider (ILC), FCC-ee, and Compact Linear Collider (CLIC) [22][23][24][25][26], and thus could be a probe of the dark matter or neutrino physics.
We show that the observed thermal dark matter abundance can be explained by (co-)annihilation of WIMPs through the neutrino portal interaction. Furthermore, this dark matter, if heavier than around 2 GeV, can be tested in future lepton colliders.
In the second part, we study the neutrino propagation in the universe with the neutrino portal interaction. We show that neutrino propagation is affected only when the invisible decay of the Higgs boson is at such a rate which is fully within the sensitivity reach of future lepton colliders. This possibility is interesting because in the IceCube neutrino observatory [27,28] the cosmic-ray neutrino event above PeV has not yet been detected, especially of the Greisen-Zatsepin-Kuzmin (GZK) neutrinos [29][30][31]. Also, the Glashow resonance [32] has not been observed. We point out that the absence of high energy cosmic-ray neutrinos can be explained if annihilation of the neutrino-(anti)neutrino into WIMPs takes place before the neutrino arrives to Earth. Namely, a cutoff can be set from the neutrino portal interaction. Moreover, neutrino mass is constrained kinematically from the mass range of WIMPs with a given cutoff, e.g. for a cutoff of a few PeV, which could explain the nonobservation of the Glashow resonance, one of the neutrino masses is within 0.01-0.2 eV. On the other hand, for a cutoff around 10 PeV, which may explain the non-detection of GZK neutrino, one of the neutrino masses is within 0.008-0.1 eV. Therefore, the neutrino mass can be "measured" kinematically through the neutrino portal interaction.
A UV model is built to justify the setup and to study the experimental constraints for heavy particles relevant for generation of the higher dimensional terms. In this model, the neutrino mass can be dominantly obtained from the neutrino portal interaction.

ϕ ψ
The paper is organized as follows. In Section 2, we explain the model with several constraints and show that or can explain the thermal dark matter abundance. In Section 3, the propagation of cosmic-ray neutrinos with the neutrino portal interaction is discussed. In Section 4, the UV model is discussed. The last section is devoted to discussion and conclusions.

A simple effective theory for WIMP
To simplify the discussion, suppose that the Lagrangian additional to the SM, , has only one generation of neutrino, where the total Lagrangian is given by ; ( ) is the mass of ( ); is the potential of the scalar fields which is supposed to give a vanishing vacuum expectation value (VEV), , and additional mass squared, , to . We will neglect the Higgs portal term, , in , because the scalar mass is lighter than 10 GeV, and is sufficiently small if , 1) where v=174 GeV is the VEV of the Higgs field. A small portal coupling larger than the order of is stable under quantum corrections, where is the cutoff scale of the model, which could be smaller than M . The other dimension-five operators,  ,  , , and are suppressed due to approximate lepton number symmetry under which L and are 1 while the others are 0. 2) In particular, we suppose that a tree-level term induces a neutrino mass that is smaller, or of the same order, as the physical one. Note that the tree-level term is not generated in a UV model if all heavy particle masses and interactions preserve the lepton number (see Section 4).

Constraint from the neutrino mass
For the broken phase of electroweak symmetry, one obtains an interaction . It was pointed out that the neutrino mass is generated at the 1-loop level in this broken phase interaction [33][34][35][36]: 1) The neutrino portal models with an efficient Higgs portal interaction are studied in [17][18][19][20].
2) The coefficient of these terms are stable under quantum corrections if they are greater than . The quantum correction for will be discussed in the following.
We have taken the renormalization scale so that this is an on-shell renormalization. Since it is constrained by the CMB and BAO observati- while it is also constrained by the double beta decay experiment for an electron neutrino [37], the neutrino mass crucially restricts the mass range of the two WIMPs. In Fig. 1, the contour plot of the generated neutrino mass and its constraint (gray shaded region) are represented in the plane with =12 MeV. One sees that is restricted to be around , and smaller the M, smaller is the difference . Since one of and is stable, M has an upper bound for sufficient (co-) annihilation of or so as not to over-close the universe. Thus, and are constrained to be nearly degenerate, ϕ Note that this constraint disappears when is replaced by a complex scalar field with only a bilinear mass term because the lepton number symmetry recovers. How-m ϕ ≃ M ψ ever, let us pursue with the simple real scalar case with , The following discussion will be qualitatively the same in a specific parameter region but with complex extension of the scalar field.

Heavy boson decays in colliders
Since 20 GeV in our consideration, the anomalous decays of the Higgs, W and Z bosons into and a lepton are possible. Thus, in colliders this scenario is constrained and tested. In particular, the Higgs boson invisible decay is represented as The decay width of the process is obtained as is the Higgs boson mass and the decay products are approximated to be massless. Given the total decay width of the Higgs boson 4 MeV, the branching ratio of this process is estimated as
On the other hand, the decay rate of boson to a charged lepton and missing energy (Z boson to missing energy) can be estimated as at the leading order of the anomalous decay 1) , where is the decay rate of the subscript at the treelevel in the SM. The branching ratio of W boson to lepton + missing energy (Z boson to missing energy) differs from the SM by .

O
One finds that M can be as small as (100) GeV to be consistent with the current experiments. To be conservative, let us set a bound 2) M ≳ 400 GeV. (9) This is represented as the horizontal black band in MeV. The purple region may be excluded due to the neutrino effective number. In the gray region, the universe is over-closed. On the orange band, the thermal abundance of the lighter WIMP explains the dark matter. The pink region may be tested in future CMB/BAO observations. ϕ, ψ M ≲ 1 1) The processes with virtual emission and absorption are also included in the decay width.

2) For
TeV, one may care for the constraint for a heavy field in some UV models. The constraint in a UV model, which will be discussed Sec.4, is represented by a lower bound (34) similar to (9).

Thermal relic abundance of WIMP ϕ ψ
Let us now discuss the thermal relic abundance of the lighter of or , which annihilates into (anti-)neutrinos through the t(u)-channel, In the first row, one does not have or , because, by integrating out , the corresponding effective vertex vanishes in the equation-of-motion for external neutrinos. The total annihilation cross-sections times the relative velocity at the tree-level are given as, for annihilations of and , respectively. The terms are calculated by FeynRules and FormCalc [40,41], which are also used to confirm all amplitude calculations in this paper. The dark matter abundance is es-timated as ⟨σ eff v⟩ where is the thermally averaged annihilation cross-section given by is the modified Bessel function of the j-th kind) [42], where the co-annihilation effect is included; is the degree of freedom of the energy (entropy) density of the radiation, which is typically around 10-100 for 10 GeV; is the freeze-out temperature in units of , which is around 15-20; h=0.678. The region satisfying is represented by the orange band in Fig. 1. In Fig. 2, the contours of are shown (orange bands) at the limit . The width of the orange bands denotes the ambiguity of our calculation 1) . The gray regions in both figures denote the over-closure of the universe, . From the figures, one finds that the thermal dark matter can be tested in future lepton colliders with mass ϕ ψ If the lighter of or is part of the dark matter, the mass range where it can be tested increases.
In particular, masses greater than 6 GeV are now being tested in the Xenon1T, LUX, and PandaX experiments [43][44][45] . The above boundary is represented as the black dotted line in Fig. 2. It is interesting to note that we can have a cross-check if the dark matter is detected in the direct-detection experiments.

N eff and BBN
The mass of or should be larger than MeV, otherwise the created neutrinos from their annihilations will change by and could spoil the BBN [46][47][48][49]. According to [47], N eff has a slightly larger and may be tested by several future CMB observations such as the PIXIE and CMB-S4 experiments, as well as the BAO observation [50][51][52].   We now focus on the region where the lighter WIMP composes a fraction of the dark matter, , i.e. the region with sufficiently strong neutrino portal interaction. The observed dark matter abundance can be explained by the other dark matter components: WIMP with neutrino portal interaction of different generation (see Section 5 and footnote 4), superpartner 1) , inflaton [76][77][78][79][80][81][82][83][84][85], etc. This region is interesting because it would affect the neutrino propagation in the universe.
Although more statistics is needed, up to now no cosmic-ray neutrino event above several PeV has been detected in the IceCube experiment [27,28], and the Glashow resonance around 6 PeV has also not been observed [32]. Despite several detections of cosmic-ray events of other kinds of particles up to ~ 10 2 EeV, this fact implies that there may be a special cutoff for cosmic-ray neutrino. In particular, if some of the observed cosmic-rays around 10 2 EeV are protons, cosmic-ray neutrinos of (EeV) should also be observed. Protons of energy larger than (10 2 ) EeV interact with CMB photons and produce pions via the resonance, and hence lose energy before cosmic-ray neutrinos reach Earth. This scattering sets a GZK cutoff at an energy (10 2 ) EeV for protons [86,87], which explains the observed cutoff for high energy cosmic-ray events. In the GZK cutoff scenario, GZK neutrinos of energy (EeV) are produced from the decay of these pions [88] and should be detected at events/year in the IceCube neutrino observatory [29][30][31].
The absence of such energetic neutrino events can be explained from a viewpoint of particle physics 2) . Thanks to the neutrino portal interaction, annihilation between cosmic-rays and cosmic background neutrinos, is enhanced for sufficiently small M, so that before the neutrinos reach Earth they turn into WIMPs 3) . Namely, we propose that the neutrino portal interaction can set a cutoff for cosmic-ray neutrinos.
To set a cutoff, there are two conditions. First, the annihilation channel should be turned on at , and hence the center-of-mass energy of the neutrino-(anti) , should be greater than the threshold, or , at , as Here, is the typical energy of cosmic background neutrinos with temperature 2×10 −4 eV, and is the angle between the momenta of two neutrinos.
Secondly, the mean free path, , imposed by annihilation, should be shorter than the distance to the neutrino source. To discuss this, let us neglect for simplicity the neutrino oscillation 4) . Following [94], one obtains the mean free path of a neutrino given by Here, is the neutrino distribution function for cosmic background neutrinos, and is the helicity averaged neutrino-(anti)neutrino annihilation cross-section, which is, with , approximated as The neutrino flux from the source at a distance L is then reduced by a factor Here, we have neglected the effect of the redshift of due to the expansion of the universe, which would reduce the observed in IceCube by (10)% with Gpc. For instance, the predicted neutrino flux is represented in Fig. 3 by assuming a neutrino flux before annihilation of O ∼ O(1) 1) There are several typical lightest superpartners (LSPs) which might be the dominant dark matter, depending on SUSY breaking scenarios: gravitino LSP in gauge mediation [53], bino-like LSP with SUSY breaking in a gauge unified manner, wino-like LSP in anomaly mediation and simple SUSY breaking scenarios based on the anomaly mediation [54][55][56][57][58][59][60][61][62], N=2 superpartners in N=2 partial breaking [63][64][65][66][67][68][69][70][71][72][73][74][75], etc.
2)There are also astronomical explanations for the absence of neutrino events above several PeVs. For example, if the neutrino is originated from the galaxy clusters or starburst galaxies the non-observation of the Glashow resonance can be accounted for [89]. If the (10 2 ) EeV cosmic-rays observed are composed of heavy nuclei such as iron, the absence of the GZK neutrino event can be explained [90].
The first term is the best-fit power law in [28] , while the second term represents a toy GZK neutrino flux, for GeV, otherwise it is 0. (Realistic ones for GZK neutrino are given in [29][30][31].) The neutrino flux in our scenario is approximated as In Fig. 3, one finds that the neutrino flux does get a cutoff or an absorption band through the t-channel annihilation. Note that the cutoff is less efficient in a model where there is a significant s-channel annihilation/scatter-ing process. In fact, the s-channel process itself does not contribute like a "cutoff", but is an absorption line due to the quick decrease of the cross-section when the centerof-mass energy exceeds the threshold. Moreover, the scattering process between neutrino and (anti)neutrino is at tree-level if the s-channel process exists. This is constrained by the CMB observation [100] and the efficiency of the cutoff is bounded. In our case, the scattering process is 1-loop suppressed and this bound is much looser than the heavy boson decay.

Relation with heavy boson decay
The contour plot of is represented in Fig. 4. From the left panel, one finds that the neutrino flux originating from a distance L > O(10) Mpc (26) can be affected by the neutrino portal interaction.
As shown in the left panel of Fig. 4, we have checked that to obtain smaller than the scale of particle horizon size ~ 10 Gpc, i.e. where neutrino propagation in the universe could be affected, M should be smaller than ~ 2 TeV. From the right panel, where M is around the lower bound (9), one reads the upper bound of to be around (100) MeV. This upper bound decreases with larger M due to the scaling of the cross-section. Hence, one obtains the parameter range where the neutrino propagation in the universe is affected, Since the upper bound of M satisfies (10), the following is predicted: if the high energy neutrino flux in IceCube is affected by the neutrino portal interaction, the Higgs invisible decay is fully within reach of future lepton colliders. Let us emphasize again that M required here is   14), and we can not provide dominant dark matter whose interaction affects cosmic-ray neutrinos. However, to explain the neutrino oscillation, an extension with several flavors of or is needed. (See conclusions and discussion.) In this case, some of the flavors could be the dominant dark matter, while some could affect the spectra of cosmic-ray neutrinos. We note that in this case the dark matter should be lighter than the particles relevant for the cutoff, and thus the dark matter mass is lighter than (100) MeV.

Measuring the physical neutrino mass range
Neutrino flux carries information on annihilation during its propagation. In particular, as in a collider, one can measure kinematically once the cutoff and the masses of and are given. This implies that a mass scale can be obtained for one of the neutrinos if its mass is greater than .
Even if just the mass range of is given, one can predict the neutrino mass range. Since , with the cutoff scale fixed, one obtains which implies that ( ) has a lower bound corresponding to Eq.(16) ( ). If at the lower bound of Eq.(16), the lower bound of one of neutrino masses is predicted.
For instance, for a cutoff at 6 PeV, one obtains a lower bound of neutrino mass eV. With a cutoff 6 PeV, which may explain the non-observation of the Glashow resonance, the lower bound for neutrino mass increases. Thus the neutrino mass range is predicted as The neutrino flux around the lower limit is illustrated by the purple dotted line in Fig. 3. This mass range covers the atmospheric neutrino scale of 0.05 eV.
If the GZK neutrino source originates from Gpc away from Earth, and should be smaller than . This can be found in the right panel of Fig. 4, because for any this is almost satisfied. Then, for a cutoff around 10 PeV, which may explain the non-detection of GZK neutrino (See the blue dashed line in Fig. 3), the neutrino mass range can be estimated as The lower bound, where annihilation is most efficient, is close to the solar neutrino scale of 0.009 eV.

UV model
In the previous sections, we studied a dimension 5 op-M ≲ TeV erator with two additional gauge singlets. It is questionable whether there is a UV model, and if there is, whether the constraints for heavy particles in the UV model restrict our scenario, especially for .
To suppress the term in order to satisfy the constraint (4) at the tree-level, the UV model should also have an approximate lepton number conservation. One of such UV models is given by where S and N are gauge singlet Weyl fermions with lepton number 1 and -1, respectively, and we have omitted the kinetic terms. For later convenience, we introduce the Yukawa coupling y, the decay constant f, the order parameter , and the mass parameter satisfying due to the approximate lepton number symmetry. We have forbidden the term at the tree-level by imposing symmetry, under which and the spurion have a charge 1/4, 1/2, 1/2 and 1/4, respectively. Leptons can carry a charge 1/2, so that Yukawa couplings are allowed. Thus, this symmetry is identified to be spontaneously broken down due to the VEV of some scalar field.
By making a shift , one finds that the neutrino portal interaction appears as The neutrino portal term (1st term) is decoupled from the heavy fields, S and N. Moreover, does not appear after integrating out the heavy fields up to 1-loop level, because does not directly couple to the heavy field. This fact can be also checked by integrating out the heavy fields in the terms of Eq. (31) after diagonalizing the fermion mass matrix. Interestingly, with this UV model, the neutrino mass is generated purely radiatively through the neutrino portal interaction. The Higgs portal term, , could be suppressed if is a pseudo-Nambu-Goldstone boson with breaking scale f (See discussion for a concrete non-linear sigma model) 1) . There are several constraints for N, because it be- haves as a right-handed neutrino with Yukawa coupling y and mass [101,102]. If we adopt the constraint in [102] for a heavy right-handed neutrino, which dominantly mixes with neutrino, is required. On the other hand, this effective theory should have from the viewpoint of perturbative unitarity, and when is identified as a pion-like field, which, respectively, turn out to be In this paper, a simplest neutrino portal interaction for WIMPs was investigated, especially for the lightest WIMP mass in the range of , where it is not yet severely constrained by direct detections. Neutrino portal interaction is interesting because it can affect not only collider physics for the Higgs boson but also neutrino physics.
We pointed out that the constraint for radiatively generated neutrino mass seriously restricts the parameter space so that the two WIMPs are nearly degenerate. Due to this restriction, the Higgs boson can decay into WIMPs plus a neutrino, and this invisible decay can be searched for in future lepton colliders, CEPC, ILC, FCC-ee and CLIC.

≳
We showed that the neutrino portal interaction can successfully (co-)annihilate the lightest WIMP, and that the WIMP relic abundance can explain the observed abundance of dark matter. Such a neutrino portal dark matter could be tested in future lepton colliders for mass 2 GeV. In the case that WIMPs explain a small fraction of the dark matter abundance, the neutrino propagation in the universe can be significantly affected. We pointed out that this region can be fully tested in future lepton colliders. In particular, this region can set a cutoff for cosmic-ray neutrino and can explain the non-detection of the GZK neutrino events, and the non-observation of the Glashow resonance in IceCube. Moreover, a neutrino mass can be "measured" kinematically from the scale of the cutoff and WIMP mass.
Using a UV model, we have justified our set up and showed that neutrino mass can be dominantly generated by the neutrino portal interaction.
Since there are generations in the SM, it is natural to make an extension of the neutrino portal interaction (1) to the case of 3 generations, e.g. , where i, j denotes the generation and is the dimensionless coupling in the mass basis of ( ). The neutrino mass matrix is generated with where is K from Eq. (3) but with ( ) replaced by the mass of ( ). In this extension, several parameter regions previously discussed can be simultaneously realized by the neutrino portal interactions of different generations.
Let us now provide a natural realization of the UV model Eq. (31) for our relevant parameter ranges where and M are sufficiently small. A light scalar suggests a naturalness problem. One of the solutions of this problem 1) is to identify as a pseudo-Nambu-Goldstone boson. Consider now the spontaneous breaking of an approximate SU(2)×U(1) global symmetry to U(1) by some non-perturbative effect, in analogy with the chiral symmetry breaking in QCD. If all explicit breaking terms of SU(2)×U(1) can be identified as spurions with even charges under the residual U(1), this residual symmetry contains an exact symmetry. The U(1) charged pion, say , is odd and contains as . This possibility not only explains the smallness of , but also allows a rather small decay constant, f, for the composite scalar , like the pion decay constant in QCD.
To be concrete, let us consider the following non-linear realized Lagrangian for pions,  Here, is SU(2)×U(1) symmetric Lagrangian, while the terms that explicitly break SU(2)×U(1) are collected in ; is a matter doublet with U(1) charge -1/2 and lepton number -1, while the fermion S carries lepton number 1; is the VEV of an SU(2) doublet operator with U(1) charge -1/2, and the second term of Eq. (36) turns out to be the second and third terms in Eq. (31). The mass parameters are the explicit breaking terms of the SU(2)×U(1) symmetry, and can be smaller than and f naturally. In particular, the un-Z 2 1) Alternatively, this may indicate that a SUSY extension of the SM has a SUSY breaking soft scale around MeVs in the odd sector, while that in the SM sector is above TeV to survive the experimental constraints. A candidate is a gauge mediation scenario [53], where sparticles charged under the SM gauge group gain weight via gauge interaction, while a singlet scalar acquires a highly suppressed mass either from higher order correction or the gravity effects. In this model, with additional light particles which are assumed to be lighter than the Higgs boson, the possibility of testing in future lepton colliders is even increased. Although the neutrino mass constraint is alleviated, and can deviate from , for given values of the mass and the cross-section for dark matter-dark matter max (M ψ , mφ ≃ m ϕ ) (neutrino-(anti)neutrino), the increase of leads to the increase of the neutrino portal coupling 1/M. Thus, the Higgs invisible decay rate is even enhanced for the regions of thermal dark matter and affects the propagation of cosmic-ray neutrino.
I would like to thank Adam Falkowski, Fapeng Huang, and Hao Zhang for collaboration at an early stage of this work. I also thank Hiroyuki Ishida and Yingnan Mao for useful discussions and thank Hiromasa Takaura for carefully reading the manuscript. Moreover, I thank the referee for carefully checking the calculations as well as pointing out typos and confusing statements.