1 Introduction


September, 2017

Dark matter stability and one-loop neutrino mass generation based on Peccei-Quinn symmetry

Daijiro Suematsu111e-mail:  suematsu@hep.s.kanazawa-u.ac.jp

Institute for Theoretical Physics, Kanazawa University, Kanazawa 920-1192, Japan

We propose a model which is a simple extension of the KSVZ invisible axion model with an inert doublet scalar. Peccei-Quinn symmetry forbids tree-level neutrino mass generation and its remnant symmetry guarantees dark matter stability. The neutrino masses are generated by one-loop effects as a result of the breaking of Peccei-Quinn symmetry through a nonrenormalizable interaction. Although the low energy effective model coincides with an original scotogenic model which contains right-handed neutrinos with large masses, it is free from the strong problem.

1 Introduction

The standard model (SM) has been confirmed by the discovery of the Higgs scalar [1]. However, it is now considered to be extended to explain several experimental and observational data such as neutrino masses and mixings [2, 3], and dark matter (DM) [4]. Strong problem is also one of such problems suggested by an experimental bound of the electric dipole moment of a neutron [5]. Invisible axion models are known to give a simple and interesting solution to it [6, 7]. The KSVZ model, which is one of such realizations, is an extension of the SM by a complex singlet scalar and a pair of colored fermions. It has a global symmetry, which is violated only by the QCD anomaly and plays a role of Peccei-Quinn (PQ) symmetry [8]. If the spontaneous breaking of this symmetry occurs, a pseudo Nambu-Goldstone boson associated to this breaking called axion appears to solve the strong problem [9]. If the axion decay constant is large enough such as due to a vacuum expectation value (VEV) of the singlet scalar, the axion mass is very small and its coupling is extremely weak so as not to cause any contradiction with experiments and astrophysical observations [10].

On the other hand, the breaking is known to cause degenerate minima for the axion potential due to the QCD anomaly depending on both the field contents and the PQ charge assignment for them. As a result, the model is generally annoyed by the dangerous production of topologically stable domain walls [11]. It can be escapable only for unless one consider the domain wall free universe brought about by inflation. If a certain subgroup of remains as a discrete symmetry broken only by the QCD anomaly in a model with , it could present an interesting scenario in relation to the DM physics at the low energy regions.aaaThe similar idea has been discussed in several articles, recently [12]. However, the present model is different from them.

In this paper, we consider such a possibility in an extension of the KSVZ model, in which an inert doublet scalar and three right-handed neutrinos are added. The low energy effective model obtained from it after the breakdown of the symmetry is reduced to the original scotogenic neutrino mass model with an effective symmetry [13]. This symmetry could guarantee the stability of a lightest neutral component of the inert doublet scalar to give a DM candidate. The neutrino masses are generated through a one-loop effect as a result of the breaking. The relevant diagram is caused by both right-handed neutrinos and a nonrenormalizable interaction between the inert doublet scalar and the ordinary Higgs doublet. The model might be recognized as a well motivated simple framework at high energy regions for the original scotogenic model.

The remaining parts are organized as follows. In the next section, we introduce a model by fixing charge assignment of to the field contents. We discuss basic features of the model such as remnant effective symmetry, scalar mass spectrum, vacuum stability and so on. In section 3, phenomenological features such as neutrino mass generation, leptogenesis and DM abundance in this model are discussed. The consistency of the scenario is also studied from a viewpoint of the vacuum stability and a cut-off scale of the model. We summarize the paper in section 4.

2 An extension of the KSVZ model

The KSVZ model is constructed by introducing a singlet complex scalar and a vector-like colored fermions to the SM [6]. We assume as triplets of the color . Although they are singlets, they could have a suitable weak hypercharge , in general. This point is crucial for phenomenological consistency of the model as discussed below. The model has a global symmetry and its charge is assigned to and , but it is not assigned to the SM contents. We assume the existence of a gauge invariant Yukawa coupling so that the PQ mechanism could work to solve the strong problem. This requires that the PQ charge of these new ingredients should satisfy . On the other hand, this symmetry should be chiral to have the QCD anomaly and is satisfied. Thus, this is spontaneously broken through the VEV of .

The transformation for the colored fermions shifts the QCD parameter through the anomaly as [5, 11]


Since has a period , the model is invariant for where is an integer and . This means that the model could have a discrete symmetry after taking account of the QCD anomaly.bbbThe axion decay constant is related with the PQ symmetry breaking scale as by using this . If we assign the charge as , the model has and no degenerate minima in the axion potential. Thus, the model has no domain wall problem as is well known.ccc Although the model has domain walls bounded by the string caused from the spontaneous breaking, it is not topologically stable and then it can shrink and decay. As a result, no cosmological difficulty appears [14]. Here, we note that an effective symmetry could remain after the symmetry breaking due to although it is violated by the QCD anomaly. Since the SM contents are supposed to have no PQ charge, it could play an important role in the leptonic sector of the model to guarantee the stability of the lightest odd field in that sector, which could be DM.

If both and cannot couple with quarks, which occurs in case for example, they are stable and then its relic abundance has to be smaller than the DM abundance [15]. Even if its relic abundance satisfies such a condition, the existence of the fractionally charged hadrons is generally forbidden by the present bound obtained from the search of fractionally charged states. On the other hand, if we assign or to , all the hadrons can have integer charge. In that case, the relic abundance will restrict the mass into a narrow range such as TeV [15]. Moreover, they are allowed to couple with quarks through a renormalizable Yukawa interaction as long as their PQ charge is zero. For example, using the left handed quark doublet and the Higgs doublet or , the coupling is allowed for with and and also for with and . In these cases, decays to the SM fields through these couplings. can also decay via the mass mixing with induced by the coupling through . As a result, the mass has no constraint other than the bound obtained through the accelerator experiments. Anyway, in the model where the PQ charge is assigned as discussed above, the strong problem could be solved without inducing any cosmological and astrophysical difficulty, as long as the symmetry breaking scale satisfies .

Now, we consider a modification of this model by introducing an inert doublet scalar and three right-handed neutrinos . The PQ charge assignment of the fields contained in the model is shown in Table 1. Invariant terms under the assumed symmetry for the Yukawa couplings and the scalar potential of the relevant fields are summarized as


where is taken to be real and is a cut-off scale of the model. The quark generation index is abbreviated in the Yukawa coupling . We find that given in eq. (2) is the most general scalar potential up to the dimension 5.

0 0
2 0 2 1

Table 1   The hypercharge and the charge of new fields in the model. The SM contents are assumed to have no PQ charge. Parity for the effective symmetry which remains after the breaking is also listed.

After the symmetry breaking due to , , and are found to get masses such as , and , respectively. Since can decay to the SM fields through the second term in as discussed above, there is no thermal relic of in the present Universe. The effective model at the scale below could be obtained by integrating out [16]. This can be done by using the equation of motion for . As its result, we obtain the corresponding effective model whose scalar potential of the light scalars can be written as


where we use the shifted parameters which are defined as


We note that the model contains the neutrino Yukawa couplings between heavy right-handed neutrinos and the inert doublet scalar as shown in the above .

Vacuum stability condition for the scalar potential in eq. (3) is known to be given as [17]


and these should be satisfied at the energy region . On the other hand, at , both the same conditions for as eq. (5) except for the last one and new conditions


should be satisfied. The couplings in both regions should be connected through eq. (4). We can examine whether these conditions could be satisfied or not by using one-loop renormalization group equations (RGEs). This is the subject studied later.

This effective model obtained after the spontaneous breaking of is just the original scotogenic model [13].dddIn the case of , and then its subgroup could be broken by the electroweak anomaly also. However, since this breaking does not induce the decay of the lightest odd field, this can be considered to be a good symmetry in the effective model. This model connects the neutrino mass generation with the DM existence. It has been extensively studied from various phenomenological view points [18, 19, 20, 21, 22]. In the present case, the right-handed neutrinos do not have their masses in a TeV region but they are considered to be much heavier. The coupling which is crucial for the one-loop neutrino mass generation is derived from a nonrenormalizable term as a result of the PQ symmetry breaking. The model contains the inert doublet scalar which has odd parity of the remnant effective . It has charged components and two neutral components . Their mass eigenvalues can be expressed as


We suppose  TeV although it requires fine tuning because of . As a result of the effective symmetry, the lightest one among the components of is stable to be a DM candidate if it is neutral. If it is supposed to be , we find that this requires and as long as is satisfied. On the other hand, since is satisfied in eq. (7), the mass eigenvalues of the components are found to be degenerate enough so that the coannihilation processes among them are expected to be effective. This observation suggests that the abundance of could be suitably suppressed and then it could be a good DM candidate as the ordinary inert doublet model [23, 24]. The charged states with the mass of  TeV are also expected to be detected in the accelerator experiments.

3 Phenomenological features

3.1 Neutrino mass, leptogenesis and DM relic abundance

In this model, neutrino masses are forbidden at tree-level. However, since both the right-handed neutrino masses and the mass difference between and are induced after the breaking, the small neutrino masses can be generated radiatively through one-loop diagrams in the same way as the original scotogenic model. Since is satisfied, the neutrino mass formula can be approximately written as


where . In order to take account of the constraints from the neutrino oscillation data in the analysis, we may fix the flavor structure of neutrino Yukawa couplings at the one which induces the tri-bimaximal mixing [19]eeeAlthough a certain modification is required to reproduce the observed mixing in the lepton sector, this simplified example could give a rather good approximation for the present purpose as found from [21].


where the charged lepton mass matrix is assumed to be diagonal. In that case, the mass eigenvalues are estimated as


where .

As is known generally and found also from this mass formula, neutrino masses could be determined only by two right-handed neutrinos. It means that the mass and neutrino Yukawa couplings of a remaining right-handed neutrino could be free from the neutrino oscillation data as long as its contribution to the neutrino mass is negligible. In eq. (10), such a situation can be realized for . This is good for the thermal leptogenesis [25] since a sufficiently small neutrino Yukawa coupling makes the out-of-equilibrium decay of the right-handed neutrino possible.fffIf we consider the TeV scale right-handed neutrinos, leptogenesis requires fine degeneracy among the right-handed neutrinos for the resonance [26]. We need not consider such a possibility in the present case. We find that the squared mass differences required by the neutrino oscillation data could be explained if we fix the parameters relevant to the neutrino masses, for example, as


for  TeV. Using these values, we can estimate the expected baryon number asymmetry through the out-of-equilibrium decay of the thermal by solving the Boltzmann equation as done in [21]. The numerical analysis shows that the required baryon number asymmetry could be generated for  GeV, which is somewhat smaller than the Davidson-Ibarra bound [27] in the ordinary thermal leptogenesis. In case of the parameter set given in (11), we find if we assume and a maximal phase in the violation parameter . In Fig. 1, we plot as a function of . Its feature can be easily understood by taking account of eq. (11). If takes larger values, the neutrino Yukawa couplings become smaller to make the violation in the decay smaller but also the washout of the generated lepton number asymmetry smaller. On the other hand, if takes smaller values, the neutrino Yukawa couplings become larger to induce the reverse effects. This makes the required baryon number asymmetry be generated only for the in the limited regions as found in this figure.

Fig. 1  Baryon number asymmetry generated through the out-of-equilibrium decay of . is plotted as a function of for the parameter set shown in eq. (11), which can explain the neutrino mass differences required by the neutrino oscillation data. Horizontal dotted lines show the required value for .

The relic abundance of is tuned to the observed value if the couplings and take suitable values. In fact, since is assumed to be of  TeV in this scenario, the mass of each component of could be degenerate enough for wide range values of and as remarked at eq. (7). This makes the coannihilation among them effective enough to reduce the abundance [21]. We search the region of and , which realizes the required DM abundance as the relic abundance by taking the values of and as the ones given below eq. (11). They are suitable for the explanation of the neutrino oscillation data and the cosmological baryon number asymmetry. In the estimation of the DM relic abundance, we follow the procedure given in [28] where the coannihilation effects are taken into account.

We present a brief review of the procedure adopted here. The relic abundance is estimated as


where is the relativistic degrees of freedom. The freeze-out temperature of and are defined as


In these formulas, the effective annihilation cross section and the effective degrees of freedom are expressed asgggIn this part, we label as .


where is the thermally averaged (co)annihilation cross section and is the thermal equilibrium number density of . If the former is expanded by the thermally averaged relative velocity as , it could be approximated only by since is satisfied for the cold DM. Final states of the relevant (co)annihilation are composed only of the SM contents. The corresponding can be approximately calculated as [21, 24]


where and is defined by using given in eq. (7) as


Fig. 2   Points plotted by a red solid line in the plane can realize the required DM relic abundance as the relic abundance. The last condition in eq.(5) is satisfied at a region above the straight line which represents for a fixed .

We use this procedure to find the points in the plane, where the required DM abundance is realized by . In Fig. 2, we plot such points by a red solid line for  TeV and which are used in the previous part. In this figure, we take account of the condition which has been already discussed in relation to eq. (7). Moreover, if we use the Higgs mass formula , we find for  GeV and then the last condition in eq. (5) can be also plotted for a fixed in the same plane.hhhWe note that the second condition in eq. (5) is automatically satisfied if the last one is fulfilled. An allowed points are contained in the region above a straight line which is fixed by an assumed value of . We give two examples here. Although the DM abundance can be satisfied for the negative value of , we find that such cases contradict with the vacuum stability condition for given in eq. (5). The figure shows that and/or are required to take rather large values for realization of the DM abundance. This suggests that the RG evolution of the scalar quartic couplings could be largely affected if they are used as initial values at the weak scale. In that case, vacuum stability and perturbativity of the model could give constraints on the model. In the next part, we focus our study on this point.

Before proceeding to this subject, we comment on the contribution of the axion to the DM abundance and also a possible violation of by the quantum gravity effect. In this model, the axion could also contribute to the DM abundance through the misalignment mechanism. If the initial misalignment of the axion is written as , the axion contribution to the present energy density is estimated as [5]


The axion contribution to the DM abundance crucially depends on the scale of and . This estimation shows that it could be too small to give the required value for  GeV even if we assume .iiiThe estimation of the relic axion abundance has to take account of the contribution from the decay of string and domain walls. Depending on it, the upper bound on the PQ breaking scale seems to be somewhat ambiguous. While one group finds that the axion production is more efficient than the misalignment case [29], the other group finds that it is less efficient than the misalignment case [30]. Thus, the axion contribution to the DM abundance is sub-dominant or negligible for  GeV. In this region of , the result obtained for through the above study can be still applicable even if the axion contribution to the DM abundance is taken into account.

Although we assume that is exact in this study, continuous global symmetry is suggested to be violated by the quantum gravity. This possible effect on the PQ mechanism has been studied [31]. If the symmetry is violated by the gravity induced effective interaction which is suppressed by the Planck scale such as


it has been shown that should be satisfied for the PQ mechanism to give a solution to the strong problem in case that is of . If accidental appearance of global happens due to some discrete or continuous gauge symmetry [32], it might protect the PQ symmetry up to sufficiently higher order operators. The same breaking effect could also affect the axion CDM abundance [31]. If the contribution to the axion mass due to the quantum gravity is small compared to the one due to the QCD anomaly,  GeV is required for saturating by the axion contribution. Even if its contribution to the axion mass is larger than the one from the QCD anomaly within the bound which is required so as not to disturb the PQ mechanism,  GeV is required again for saturating . Thus, could play a dominant role in the DM abundance as long as is smaller than  GeV.

The stability for could be also violated through the same effect. The most effective processes for the decay are induced by nonrenormalizable Yukawa couplings such as


If the allowed dimension for these kind of operators is the same as the one which guarantees the PQ mechanism to work, the lifetime of could be longer than the age of our universe in case  TeV and as long as we take  GeV. If the lower dimension operators such as are allowed, its lifetime cannot be long enough to be the DM at the present universe.

3.2 Consistency of the scenario with a cut-off scale of the model

It is crucial to check what kind of values of the right-handed neutrino mass and could be consistent with a value of which is restricted by the axion physics. In this model, DM is identified with whose mass is of  TeV. In such a mass region, we find that its abundance is determined by the values of the scalar quartic couplings and . On the other hand, these couplings could affect the vacuum stability and also the perturbativity of the model through the radiative effects on the scalar quartic couplings . Here, we examine the consistency of the values of and required to realize of the DM abundance with these issues.jjjThe constraint due to the vacuum stability and the perturbativity is taken into account in the DM study of the inert doublet model on the basis of a different viewpoint from the present one [23, 24]. The consistency between fermionic DM and the vacuum stability is also studied in the scotogenic model [20, 33]. Since the breaking of the perturbativity is considered to be relevant to a scale for the applicability of the model, we could obtain an information for the cut-off scale . It allows us to judge whether the required value for by the neutrino masses and the leptogenesis could be induced through the VEV of .

The one-loop -functions for the scalar quartic couplings in the effective model at energy regions below are given as follows [34],


where is defined as . In these equations, we can expect that the positive contributions of and to the -functions of tend to save the model from violating the first two vacuum stability conditions in eq. (5). On the other hand, the same contributions of and could induce the breaking of the perturbativity of the model at a rather low energy scale since they could give large positive contributions to , and . Here, we identify a cut-off scale of the model with a scale where any of the perturbativity conditions and is violated.kkkSince the Landau pole appearing scale is expected to be near to this , it seems to be natural to identify with a cut-off scale of the model. In this case, should be satisfied. If is smaller than , the consistency of the scenario is lost.

We analyze this issue by solving the above one-loop RGEs at and also the ones at , which are given in Appendix. The quartic couplings in the tree-level potential at the energy scale are connected with the ones at through eq. (4). Since the masses of the right-handed neutrinos are considered to be heavy in the present model, they decouple at the scale to be irrelevant to the RGEs there. On the other hand, the mass of the colored fields can take any values larger than 1 TeV as discussed before, they can contribute to the RGEs at larger scales than their mass. In the present study, we assume that is light of  TeV but its Yukawa coupling with the ordinary quarks is small enough.lllIn the light case, study of the bound for this Yukawa coupling is an interesting subject related to the search of mixing with the ordinary quarks. However, it is beyond the scope of the present study and we do not discuss it here. Thus, they are considered to contribute substantially only to the -functions of the gauge couplings. In this study, we take its hypercharge as as shown in Table 1.

The free parameters in the scalar potential of the effective model (3) are and at as long as we assume  TeV.mmmQuartic couplings for are fixed as and at in the present study. As easily found from RGEs, larger values of make smaller. Among them, we should fix at a value used in the discussion of the neutrino mass and the leptogenesis. Both and are fixed at values determined through the DM relic abundance as shown in Fig. 2. We also have from the Higgs mass. From this point of view, is an only remaining parameter. Thus, if we solve the RGEs varying the value of for other fixed parameters, we can find checking the vacuum stability for each .

In the left panel of Fig. 3, as an example, we present the running of the scalar quartic couplings for the initial values , and at by assuming the breaking scale as  GeV. In the same panel, we also plot the value of as , which corresponds to the last one in eq. (5). In this example, we can see that the vacuum stability is kept until the cut-off scale  GeV. These values of and can naturally realize the assumed value for through the relation given in eq.(4) just by taking as a value of . This feature can be verified for other allowed values of and . Here, we note that the axion contribution to the DM abundance can be neglected for a value such as  GeV. In the right panel of Fig. 3, we plot as a function of for four sets of which are shown by black bulbs in Fig. 2. End points found in the two lines represent the value of for which the vacuum stability is violated before reaching . This figure shows that which is restricted to a rather narrow region can make appropriate values in order to realize a required value of for  GeV. This study suggests that the scenario could work well without strict tuning of the relevant parameters.

As found from the above study, the simultaneous explanation of the neutrino masses and the DM abundance could be preserved in this extended model in the same way as in the original scotogenic model. We should stress that no other additional constraint from the DM physics and the neutrino physics is brought about by taking the present scenario. The cosmological baryon number asymmetry is expected to be explained through the out-of-equilibrium decay of the lightest right-handed neutrino. The required right-handed neutrino mass could be smaller compared with the Davidson-Ibarra bound in the ordinary thermal leptogenesis [27]. This is consistent with the result in [21] where the mass bound of the right-handed neutrino for the successful leptogenesis is shown to be relaxed in the radiative neutrino mass model in comparison with the ordinary seesaw model .

Finally, we give brief comments on possible experimental signatures of the model. The present model might be examined through (i) the search of the DM and the charged scalars through the DM direct detection experiments and the accelerator experiments, (ii) the search of the mixing of with the ordinary quarks although it could be observed only in the light case, and (iii) the search of the axion whose coupling with photon is characterized by , where is the hypercharge of [15].

4 Summary

We have proposed an extension of the KSVZ invisible axion model so as to include a DM candidate and explain the small neutrino masses. An extra inert doublet scalar and three right-handed neutrinos are introduced as new ingredients. After the symmetry breaking, its subgroup could remain as a remnant effective symmetry, which is violated through the QCD anomaly but it can play the same role as the in the scotogenic neutrino mass model. Since only the new ones and have its odd parity, the model reduces to the scotogenic model which has in the leptonic sector. The neutrino masses are generated at one-loop level and the DM abundance can be explained by the thermal relics of the neutral component of . The cosmological baryon number asymmetry could be generated through the out-of-equilibrium decay of a right-handed neutrino in the same way as the ordinary thermal leptogenesis in the tree-level seesaw model. However, the bound for the right-handed neutrino mass can be relaxed in this model. Since this simple extension can relate the strong problem to the origin of neutrino masses and DM, it may be a promising extension of both the KSVZ model and the scotogenic model.


The -function for the scalar quartic couplings at are given as