arxiv: v2 [hep-ph] 7 Nov 2018

Transcription

1 DO-TH 18/12 CP3-Origins DRF90 Scaar Dark Matter, GUT baryogenesis and Radiative neutrino mass arxiv: v2 [hep-ph] 7 ov 2018 Wei-Chih Huang CP 3 -Origins, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark, Fakutät für Physik, Technische Universität Dortmund, Dortmund, Germany Heinrich Päs and Sinan Zeißner Fakutät für Physik, Technische Universität Dortmund, Dortmund, Germany Abstract We investigate an interesting correation among dark matter phenomenoogy, neutrino mass generation and GUT baryogenesis, based on the scotogenic mode. The mode contains additiona right-handed neutrinos and a second Higgs doubet Φ, both of which are odd under an imposed Z 2 symmetry. The neutra component of Φ, i.e. the ightest of the Z 2 -odd partices, is the dark matter candidate. Due to a Yukawa couping invoving Φ, and the Standard Mode eptons, the epton asymmetry is converted into the dark matter asymmetry so that a non-vanishing B L asymmetry can arise from (B L)-conserving GUT baryogenesis, eading to a nonzero baryon asymmetry after the sphaerons decoupe. On the other hand, Φ can aso generate neutrino masses radiativey. In other words, the existence of Φ as the dark matter candidate resuscitates GUT baryogenesis and reaizes neutrino masses. Eectronic address: huang@cp3.sdu.dk Eectronic address: heinrich.paes@tu-dortmund.de Eectronic address: sinan.zeissner@tu-dortmund.de 1

2 I. ITRODUCTIO The origin of the observed baryon asymmetry can not be accounted for within the Standard Mode (SM) and is one of the unresoved issues in partice physics and cosmoogy. The simpest Grand Unified Theory (GUT) based on the SU(5) mode, proposed by Georgi and Gashow in 1974 [1], features eptoquark gauge bosons which do mediate baryon number vioating processes, eading to proton decay. The mode, however, conserves the difference between the baryon and epton number B L. In other words, any generation of a baryon asymmetry from heavy gauge or Higgs boson decays, as discussed in Refs. [2 5], comes with an equa amount of epton asymmetry. These baryon and epton asymmetries, however, wi be washed out competey by non-perturbative sphaeron processes [6 8], which come into therma equiibrium when the temperature of the universe drops roughy beow GeV. The B L symmetry conservation aso exists in arger symmetry groups, such as SO(10), where the abeian U(1) B L is a subgroup. Therefore, as ong as U(1) B L is not broken when a baryon asymmetry is created, i.e., initiay B + L 0 but B L = 0, such a baryon asymmetry wi not survive the sphaeron processes. In principe, there are at east two ways to revive GUT baryogenesis. First, nonzero B L can sti be reaized in certain matter representations under SO(10) or arger groups as demonstrated, for instance, in Refs [9 14]. Second, Fukugita and Yanagida [15] (and a recent update, Ref. [16]) have proposed to incude right-handed neutrinos to resuscitate GUT baryogenesis, where the right-handed neutrino can be embedded into SU(5) as a singet or into the 16 of SO(10). A Majorana mass of, which can arise from the spontaneous symmetry breaking of L via the vacuum expectation vaue of a scaar or can simpy be imposed by hand, expicity vioates the origina B L symmetry. H H ν λ Φ 0 Φ 0 ν FIG. 1: eutrino mass generation in the scotogenic mode. Here, is the Yukawa couping and λ denotes the quartic couping between the Higgs (H) and the second doubet (Φ). In this paper, we revisit and extend the idea of Fukugita and Yanagida [15] in the context of the scotogenic mode [17]. In this mode, a second scaar SU(2) L doubet Φ is introduced which radiativey generates neutrino masses as shown in Fig. 1. At the same time, the neutra component of the doubet is a suitabe dark matter (DM) candidate because of an imposed Z 2 symmetry. 2

3 β Φ ν i Z ν i U iα W ± U iβ + Φ α β Φ Φ U iβ + + FIG. 2: DM (Φ) mediated epton number vioating processes which transfer a L asymmetry into that of + + f f In that both Φ and right-handed neutrinos are Z 2 -odd, ig the type-i seesaw ig Yukawa couping y H is forbidden (: SM epton doubet) but a new Yukawa couping Φ is aowed, which induces a washout of epton number. As iustrated h 0 in Fig. 2, the change of the h 0 epton number is accompanied by a change of the Φ number, L = Φ. In other words, the L asymmetry is transferred into a DM asymmetry. Moreover, part of the DM asymmetry f f further shifts to anigasymmetry ofigthe Higgs doubet (note that both doubets are equay charged under the SM gauge groups) because of the Φ H interactions: H ΦH and Φ ( ) Φ ( ) h 0 H ( ) H ( ). In this scenario, the L hasymmetry 0 can be maximay reduced down to one third of the initia vaue (instead of one-haf in the case without Φ1 where ony the y H couping exists [15, 16]) since Φ and H share the asymmetry. That is, the resuting fina B L asymmetry can be maximay one third of the initia B + L asymmetry generated by GUT baryogenesis. Taking into account the top (bottom) Yukawa couping, which is in therma equiibrium for temperatures T (10 12 ) GeV, the H asymmetry wi be transferred into quarks, eading to a arger epton 1 number washout. See Ref. [16] for more detais. ote that there exist many modes that reaize radiative neutrino masses and DM with discrete symmetries. It has been shown [18] in some of these modes incuding the scotogenic mode, the dark parity, used to protect the DM stabiity, is actuay reated to epton number L as ( 1) L+2j, where j is the partice spin. In other words, L and the DM parity are correated and that is the reason why the epton asymmetry is converted into the DM asymmetry in this work. If DM decoupes from the therma bath before the eectroweak phase transition (EWPT), the DM asymmetry wi be transformed back to H via the process ΦH H, which has ony a singe power of Botzmann suppression and is very efficient compared to the douby Botzmann-suppressed annihiation channes Φ HH and Φ ( ) Φ ( ) H ( ) H ( ). This conversion wi sighty increase the fina baryon asymmetry because the Yukawa coupings and the sphaerons wi redistribute the asymmetries among quarks, eptons and Higgs bosons. ote that after the EWPT, the asymmetry for the rea part of the neutra component H 0 wi vanish because of the Higgs vacuum expectation vaue [19], whereas the remaining degrees of freedom of H wi become the ongitudina component of W ± and Z. Simiary, the Φ 0 (neutra components of Φ) asymmetry wi aso vanish after the EWPT due to the efficient Φ H interactions 1 whie the Φ ± (charged components) asymmetry wi move to W ±. The fina 1 The interaction Φ 0 h Φ 0 h (h: SM Higgs boson after the EWPT) wi erase the Φ 0( ) asymmetry. 3

4 (B L)-conserving GUT baryogenesis L washouts and Φ H equiibrium After (B + L)-vioation sphaerons and EWPT L washouts (B L)-conservingΦ H equiibrium L washouts GUT baryogenesis After (B + L)-vioation sphaerons andφewpt H equiibrium B L H B L H B L H A B L FIG. 3: Pictoria iustration of asymmetry conversion in the presence of DM Φ and the righthanded neutrino. The L asymmetry generated from GUT baryogenesis is converted into a DM (Φ) asymmetry, and then is aso shared by H due to Φ H equiibrium. As a resut, the maxima B L asymmetry His one third of Bthe initia L B + ΦL asymmetry H from GUT B baryogenesisφas indicated H B L H B L H B L in the midde pane. If DM decoupes before the EWPT, the asymmetry wi be transferred back to the SM sector, increasing the fina B asymmetry as dispayed in the right pane, where the soid (dashed) Φ ( ) H ( ) Φ H ine corresponds to DM freeze-out before (after) the EWPT. See the text for more detais. L H λ 5 λ 3,4 B and L asymmetries wi stayφ ( ) unchanged since H ( ) the sphaerons become H ineffective after the EWPT. Fig. 3 eucidates the asymmetry transformation as a function of time. On the other hand, the DM reic abundance is mainy determined by the Higgs-DM coupings for TeV DM (such that DM freezes out prior to the EWPT) as shown in Fig. 4. Φ ( ) H ( ) Φ H Φ ( ) H ( ) Φ H Φ ( ) λ 5 λ 3,4 H ( ) H Φ ( ) λ 5 λ 3,4 H ( ) H FIG. 4: 1 Main DM annihiation processes which determine the DM density. The λ s denote the Higgs-DM quartic coupings, defined in Eq. (II.2). ote that the scotogenic mode aone can generate the baryon asymmetry, apart from reaizing neutrino masses and accommodating DM candidates, via heavy neutrino decays as first pointed out by Ref. [20] and foowed by more detaied studies [21 24]. The subject has been further deveoped recenty Ref. [25] which attains ow-scae eptogenesis without any degeneracy in the right-handed neutrino mass spectrum and Ref. [26] which features the KeV right-handed neutrino as a DM candidate. In this work, we instead focus on the L washout effects induced by the heavy neutrinos, and assume the new Yukawa couping is CP-conserving. In other words, the right-handed neutrino decays equay into eptons and 1 1 4

5 antieptons 2. This paper is organized as foows. In Section II, we briefy review the scotogenic mode and then deveop the formaism for epton number washout based on Botzmann equations in Section III. In Section IV, we expain how asymmetries are transferred between the DM and SM sectors and present our numerica resuts of the Botzmann equations. The reic abundance is cacuated in Section V where the DM direct search bounds from the XEO1T experiment are aso taken into account. Finay, we concude in Section VI. II. SCOTOGEIC MODEL The scotogenic mode has been proposed by E. Ma [17], where the neutrino mass is oopinduced by a second SU(2) L doubet scaar Φ and the right-handed neutrinos, both of which are odd under an imposed Z 2 symmetry. Thus, Yukawa coupings in the type-i seesaw, y ij i H j are forbidden and repaced by y ij i j. In principe both the Z 2 -odd and the neutra component of the Φ doubet coud be the DM candidates. However, the mass of being of interest for this work is above GeV, that is too heavy to thermay generate the correct reic density [28]. In the framework of SU(5), Φ can be embedded into the representation of 5, whie can be a singet. We here simpy assume that other partices, which are embedded in the same representation of SU(5) (or arger symmetry groups) as SM partices or Φ, are heavier than the scae of interest. Thus, ony the SM partices, Φ and are taken into account in the anaysis. In addition to the SM interactions, the Lagrangian reads with L ij i j + M k 2 c k k + V (Φ, H), (II.1) V (H, Φ) = µ 2 1 H 2 + µ 2 2 Φ 2 + λ 1 H 4 + λ 2 Φ 4 + λ 3 H 2 Φ 2 + λ 4 H Φ 2 + λ 5 2 ( (H Φ) 2 + h.c. ), (II.2) which is just the scaar potentia of the inert Higgs Doubet mode [29]. The radiative neutrino 2 ote that, whie this paper focuses on the asymmetry transfer between DM and eptons without considering the decay contribution, in the presence of CP vioation in the Yukawa coupings, the asymmetry from decays coud be sizabe. As pointed out in Ref. [27], tree-eve -mediated washout processes are out of equiibrium during the time of decays if m 10 7 GeV such that the L asymmetry from the decays can survive from washouts and account for the observed baryon asymmetry. The region of interest in this work, m GeV, fas into this region and hence the decay effect coud be important, depending on the vaues of the Yukawa coupings and the size of the CP phase(s). 5

6 mass matrix induced by oops of Φ and is [17] 3 where (m ν ) ij = k ( ) y ik y jk Mk 32π 2 ( m 2 R m 2 R M 2 k og m2 R M 2 k m 2 R = µ (λ 3 + λ 4 + λ 5 ) v 2, ) m 2 I m 2 I M og m2 I, (II.3) 2 k M 2 k m 2 I = µ (λ 3 + λ 4 λ 5 ) v 2, (II.4) with v = 246 GeV being the Higgs vacuum expectation vaue. ote that in order to obtain a non-vanishing neutrino mass, one must have m R m I, i.e., λ 5 0. We here are interested in the region of M GeV, m R m I TeV and m R m I m I m R. In this case, the neutrino mass matrix becomes (m ν ) ij = λ ( 5v 2 y ik ( ( jk) M 2 ) ) k og 1, (II.5) 32π 2 M k m 2 0 k where m 0 = m R+m I 2. To reproduce the observed neutrino mass squared difference responsibe for atmospheric neutrino osciations, the heaviest neutrino must be heavier than 0.05 ev or so, which corresponds to λ for M GeV and m 0 TeV, given of O(1). III. WASHOUT FORMALISM Due to the Hubbe expansion, a convenient quantity to describe the partice number density is Y n/s, which is the number density normaized to the entropy density s, i.e., the number per co-moving voume. The density Y is conserved in the absence of partice creation or annihiation. The Botzmann equation of a partice for an interaction a 1 a n f 1 f m is, zhs dy dz = a i,f j [a 1 a n f 1 f m ], (III.1) where H is the Hubbe parameter, z = M /T, and [a 1 a n f 1 f m ] = n n a1 n an n eq a 1 n eq γ eq (a 1 a n f 1 f m ) a n neq n f 1 n fm n eq f 1 n eq γ eq (f 1 f m a 1 a n ). f m (III.2) 3 ote that there is a factor of 1/2 missing in Ref. [17]; see, e.g., version 1 of Ref. [30] or Appendix C of Ref. [31]. 6

7 The therma rate γ eq is defined as γ eq (a 1 a n f 1 f m ) = [ d 3 p 2E (2π) 3 e E T ] f j [ d 3 p fj 2E fj (2π) 3 ] (2π) 4 δ 4 (p + a i [ d 3 p Ea ai i ] 2E ai (2π) 3 e T n p ai i=1 m p fj ) M 2, (III.3) where M 2 is the squared ampitude summing over initia and fina spins. To simpify the anaysis, we consider a 1+1 scenario, i.e., one generation of the SM eptons and one right-handed neutrino 4. Moreover, we assume that the scae of GUT baryogenesis is sighty beow the right-handed neutrino mass to avoid compications from finite-temperature effects (if, for exampe, the decay HL woud be kinematicay forbidden, the first processes in Fig. 2 woud not have resonance anymore, reducing the L washout effect) due to therma masses when T m [32]. For the L washout computation, we incude both L = 1 and L = 2 interactions. Foowing the notation of Ref. [32], the L = 2 washout processes incude Φ (with therma rate γ s ) and (γ t ) as dispayed in Fig. 2. The reevant L = 1 washout processes are Φ (γ D ), A (γ As ), Φ A (γ At1 ) and A (γ At2 ). We refer readers to our previous work [16] and references therein for more detais. ote that the previous work is based on the type-i seesaw mechanism whie in this work, it is another Yukawa couping that is responsibe for the washout processes. The formaism of washout computation is, however, simiar for the two cases. The resuting Botzmann equations incuding the epton washout and sphaerons [33, 34] processes read zhs dy B L dz zhs dy B+L dz zhs dy ( ) Y YB+L Y B L =2 2γ s + 4γ t + γ As Y eq + γ At1 + γ At2 2Y eq L [ ] Y YΦ 2b Φ γ s + 4γ t + γ As + γ At1 + γ At2 Y eq Y eq, (III.4) Φ ( ) Y YB+L Y B L = 2 2γ s + 4γ t + γ As Y eq + γ At1 + γ At2 2Y eq L [ ] Y YΦ + 2b Φ 2γ s + 4γ t + γ As + γ At1 + γ At2 Y eq Y eq M s Φ 2 α5 W Y B+L, (III.5) z ( dz = (γ Y D + 4γ As + 4γ At1 + 4γ At2 ) dy Φ dz dy B L dz, Y eq j=1 ) 1, (III.6) (III.7) where Y L(B) Y epton (baryon) Y anti-epton (anti-baryon) and Y eq is the equiibrium density of the corresponding (anti-)partice. The impact of the t- and b-yukawa coupings on the washout 4 For simpicity, we stick to the cases where the initia epton asymmetry is stored in the epton doubet. 7

8 processes can be characterized by the factor b Φ [16]: b Φ = { T GeV 1 5 T GeV. (III.8) Moreover, the chemica equiibrium for Φ ( ) Φ ( ) H ( ) H ( ) is reached if λ 2 5T 8π T 2 m P, (III.9) with m P = GeV. The Φ H chemica equiibrium is aways fufied for vaues of λ 5 of interest. The fina B and L asymmetries as functions of the fina B L asymmetry are Y fina B = c s Y fina B L, Y fina L = (c s 1)Y fina B L. (III.10) For non-supersymmetric modes the sphaeron conversion factor is c s = 28/79 [19, 35] if DM decoupes before the EWPT. On the other hand, if DM freezes out after the EWPT, the Φ 0 (Φ ± ) asymmetry wi just vanish (transfer into W ± ), and has no infuence on the fina B and L asymmetries as expained above. In this case, one has c s = 8/23 as we sha see beow. IV. ASYMMETRY TRASFER BETWEE DM AD SM SECTORS We now are in the position to expain how the washout processes can create a nonzero B L asymmetry and how asymmetries are transferred among different partices. For temperatures above GeV, the (B + L)-vioating 5 sphaerons are not in therma equiibrium and part of the epton asymmetry is moved to DM due to the washout processes induced by the Yukawa couping Φ. For both of the L = 1 and L = 2 interactions, the change in the epton number comes with an equa amount of the DM number change. ( The partia asymmetry of Φ is further converted into H through the interaction, λ 5 2 (H Φ) 2 + h.c. ). That is, after L washout one obtains a nonzero B L asymmetry: (B L) = ( Φ+ H). For the washout cacuation, m 0 = 5 TeV and λ 5 = 1 are assumed which guarantee that the Φ H interactions are aways in chemica equiibrium during the period of washout, i.e., (B L) = 2 Φ. The numerica resuts are presented in Fig. 5 with the initia B + L asymmetry from GUT baryogenesis injected at the scae of M /3 (eft pane) and M /10 (right pane). The contours represent the ratio of the fina B L to the initia B + L asymmetry, i.e., YB L fina initia /YB+L. A smaer B + L injection scae impies a shorter L washout period before the sphaerons kick in, and hence requires a arger Yukawa couping (a higher washout rate) to compensate. As a resut, the YB L fina initia /YB+L contours move upward in the right pane when compared to the eft one. 5 In the foowing, we wi use the shorthand notations B + L, B and L for (B + L)-, B- and L-vioating, respectivey. 8

9 Maxima L washout (the maxima fina B L asymmetry) denoted by dark red areas arises from the case in which the /L processes are very efficient before the sphaerons come into pay but become ineffective when the sphaerons are in therma equiibrium. Two minima B L asymmetry scenarios (white areas) correspond to situations where first L washouts have never been fast enough before the sphaerons destroy most of the initia B + L asymmetry, and second both /L and B + L processes are present and effective for a ong time, eading to vanishing B and L asymmetries 6. The back soid ine indicates the active neutrino mass m ν of 0.23 ev, the bound from Panck [36] on the sum of the active neutrino masses, whie the back dashed ine corresponds to m ν = m 2 atm 0.05 ev and the back dotted ine for m ν = m 2 so ev. If λ 5 is increased (decreased), according to Eq. (II.5) the back ines wi move downwards (upwards) accordingy. On the other hand, the active neutrino masses are not very sensitive to the precise vaue of m Φ due to the ogarithmic dependence. In order to obtain m ν = 0.05 ev and YB L fina initia /YB+L 10 2, M has to be roughy above GeV with In our previous work [16] with the type-i seesaw Yukawa couping H, one can achieve arger washout effects with YB L fina initia /YB+L 10 1 and at the same time reproduce m ν = The main difference in the presence of Φ is that the active neutrino mass is oop-induced and hence a arger Yukawa couping is needed. In this case, the washout processes ast for a onger time and coexist with the B+L sphaerons, eading to a smaer B L asymmetry. Depending on the initia B + L asymmetry, there exist regions of the parameter space capabe of reproducing YB L fina to account for the observed baryon asymmetry, YB fina = [36]. Assuming that, for exampe, the initia B + L asymmetry is of order 10 6 and the B + L injection scae is M /10, M can be as ow as GeV to reaize both the baryon asymmetry and the neutrino mass. In this case L washouts can sti be efficient enough to generate a non-vanishing B L asymmetry before the sphaerons destroy the entire B + L asymmetry. When the temperature drops beow GeV and becomes much smaer than M, L washouts are ineffective but the sphaeron processes start to destroy the B + L asymmetry. Later on, SM Yukawa coupings reach equiibrium to rearrange the asymmetry among eptons, quarks and the Higgs boson. One can repeat the anaysis of chemica equiibrium done in Refs. [19, 35, 37], incuding an extra constraint, µ Φ = µ H. To simpify the anaysis, we assume universa chemica potentias µ and µ er for the three eft-handed epton doubets and three right-handed eptons, respectivey, and a the Yukawa coupings are in therma equiibrium. This yieds µ q = 1 3 µ, µ ur = 1 6 µ, µ dr = 5 6 µ, µ er = 1 2 µ, µ H = 1 2 µ, µ Φ = 1 2 µ, (IV.1) (IV.2) 6 Again, the detaied anaysis can be found in our previous work, Ref. [16], with the different partice contents but with a simiar washout mechanism. 9

10 y' mν 0.23eV mν 0.050eV mν eV y' mν 0.23eV mν 0.050eV mν eV M GeV M GeV 10 7 FIG. 5: Contour pots of YB L fina /Y initia B+L with λ 5 set to unity and m 0 = 5 TeV. The eft (right) pane refers to the case of the initia B + L asymmetry being created at the scae of M /3 (M /10). The back soid, dashed and dotted ines present the neutrino mass of 0.23 ev, m 2 atm 0.05 ev and m 2 so ev. See the text for more detais. and thus the fina B and L asymmetries are B f = 8 (B L), 23 L f = 15 (B L), (IV.3) 23 which is different from the case in the absence of Φ with B f = 28 (B L) and L 79 f = 51 (B L) [19, 35, 37] for the SM. That is to say, Φ shares the asymmetry and sighty 79 reduces the baryon asymmetry for a given B L asymmetry. Finay, when the temperature fas beow m Φ, Φ begins to freeze out of the therma bath. The DM reic density wi be mainy determined by the quartic coupings λ 3,4,5 in Eq. (II.2), if they are arge compared to the gauge coupings and Yukawa coupings. In other words, the DM partice dominanty annihiates into the Higgs bosons. The interaction terms of λ 3 and λ 4 apparenty wi not change any asymmetries in Φ and H whie the λ 5 term, corresponding to Φ H H and Φ ( ) Φ ( ) H ( ) H ( ), wi shift the asymmetry from Φ to H. ote that the interaction Φ H H is aways much faster than the DM annihiation processes if λ 3,4,5 are of the same order. That is because the former interaction is singy Botzmannsuppressed but the atter ones are douby suppressed. The asymmetry conversion between Φ and H during freeze-out can be understood in the foowing simpe ways. Since Φ and H carry the same U(1) Y charge, the disappearance of Φ has to be compensated by the equa amount of H so that the tota U(1) Y is conserved. In the case where DM freezes out before the EWPT, the Φ asymmetry wi be transformed into that of H and further into those of the quarks and eptons. On the other hand, if DM freeze-out takes pace after the EWPT, due to the Φ H interactions the Φ 0 asymmetry wi 10

11 simpy vanish whie the Φ ± asymmetry wi transfer to that of W ±. Due to the fact that the sphaeron effects are not effective anymore beow the EWPT, both the L and B asymmetries are conserved quantities independent of the Φ asymmetry. The fina baryon number wi sighty increase by 2% if DM decoupes before the EWPT and hence the asymmetry conversion occurs. We woud ike to emphasize that regardess of the decouping time of DM, the fina DM abundance is not reated to the baryon asymmetry, even if the initia DM asymmetry is cosey connected to the initia B L (aso B) asymmetry. This is the price we have to pay in order to radiativey generate non-zero active neutrino masses via a non-zero λ 5. If Φ decoupes before the EWSB the interaction of λ 5 quicky shifts the Φ asymmetry into that of H as the density of Φ and decrease during freeze-out. Hence the fina density of Φ is ony determined by the annihiation of Φ and, simiar to symmetric DM scenarios. If Φ decoupes after the EWPT the asymmetry stored in Φ 0 and Φ 0 just vanishes due to the Φ H interactions as expained above. In addition, a non-zero λ 5 wi resut in a mass spitting between the two neutra components as indicated in Eq. (II.4). Thus the ightest neutra component is the DM partice, which is rea and is its own antipartice. The fina DM density wi ony be determined by the DM annihiation into two Higgs bosons. ote that a zero λ 5 woud yied correation between the fina DM abundance and the baryon asymmetry. In this case, the DM mass has to be around 5 GeV to reproduce the correct reic abundance. As Φ is a SU(2) doubet, the Z boson can decay into Φ, increasing the Z decay width. This, however, wi be excuded by the LEP bound. In other words, an asymmetric DM scenario cannot be reaized in this framework. V. DM RELIC DESITY AD DIRECT DETECTIO In this section, we compute the DM reic density and discuss direct search bounds. The study of DM phenomenoogy for inert Higgs doubet modes after eectroweak symmetry breaking has been studied, for instance, in Refs. [38 42], whie annihiation cross-sections in an unbroken phase have been computed in Ref. [43]. We here focus on the scenario in the atter case where TeV Φ decoupes before the EWPT and the main annihiation channes are Φ HH and Φ ( ) Φ ( ) H ( ) H as shown in Fig. 4. As we sha see ater, to achieve the correct DM density, the DM-Higgs coupings λ s have to be arger than unity and aso than the gauge and Yukawa coupings. Thus, we negect the gauge and fermion fina states in the computation. Since the Φ asymmetry is basicay zero during (and after) freeze-out, the computation of the DM reic density Ω Φ + Ω is essentiay the same as in the standard symmetric DM scenario and can be we approximated [44, 45] by Ω DM h cm 3 sec 1 σv Φ HH. (V.1) Here the thermay-averaged annihiation cross-section mutipied by the DM reative veocity 11

12 is σv Φ HH λ π m 2 Φ (V.2) where we assume λ 3 λ 5 and λ 4 = 0 for simpicity. ote that the mass degeneracy among components of Φ wi be ifted after eectroweak symmetry breaking. Heavy components of Φ wi decay into the ightest one but the tota reic density stays constant due to the unbroken Z 2 symmetry. In fact, for TeV DM the contribution from the Higgs vacuum expectation vaue to the DM mass is negigibe as can be seen from Eq. (II.4), i.e., m Φ m 0. λ λ = π Ω = Φ [ ] FIG. 6: Quartic couping λ 3 versus DM mass m Φ. The bue ine corresponds to the observed reic density whie the red dashed area is excuded by the XEO1T direct search resut [46]. The purpe ine represents the perturbativity imit 4π. On the other hand, Φ can interact with nuceons through the Higgs exchange and nu resuts from DM direct searches put constraints on the DM-Higgs couping λ 3. Again with the assumption of λ 3 λ 5 and λ 4 = 0, the DM-nuceon spin-independent cross-section is [47] σ SI = λ2 3f 2 4 π µ 2 m 2 n, (V.3) m 4 h m2 Φ where f = 0.3, µ = m n m Φ /(m n + m Φ ) and m n is the nuceon mass. In Fig. 6, we show the direct search bound from the XEO1T resut [46] denoted by the red ine 7 and the bue ine corresponds to the correct reic density, whie the purpe ine is the perturbativity imit. It is cear that XEO1T is unabe to probe the arge DM mass region as the DM-nuceon cross-section is inversey proportiona to the DM mass, eading to ow sensitivity. In addition, the DM annihiation cross-section is aso suppressed by the DM 7 The PandaX-II [48] and LUX [49] experiments yied simiar imits, whie the atest XEO1T resut [50] ony presents the bound for DM beow 1 TeV. 12

13 mass and λ 3 has to be arge in order to reproduce the correct DM density. Thus, for a arge DM mass m DM 19 TeV the theory is not perturbative anymore. This roughy agrees with the resut of Ref. [43], where 22.4 TeV is obtained by considering a contributions incuding the gauge bosons. VI. COCLUSIOS In this work, we have expored an interesting correation between DM, radiative neutrino masses and GUT baryogenesis, based on the scotogenic mode [17]. The mode contains a second Higgs doubet Φ together with right-handed neutrinos, both of which are odd under a Z 2 symmetry. The ightest one of the Z 2 -odd partices, Φ, is a DM candidate. Due to the Z 2 symmetry, the type-i seesaw Yukawa couping of the right-handed neutrinos to the Higgs boson is prohibited but a new couping Φ ( is the SM epton doubet) is aowed. Consequenty, the neutrino mass is radiativey induced by oops of Φ and. In the context of (B L)-preserving GUT baryogenesis, the additiona interaction Φ via -exchange shifts the L asymmetry into Φ such that a nonzero B L asymmetry can be generated. The net B L asymmetry wi be preserved by the (B + L)-vioating sphaeron effects and as a resut the observed baryon asymmetry can be obtained. Moreover, due to the interactions H ΦH and Φ ( ) Φ ( ) H ( ) H ( ), the asymmetry in Φ from L washouts wi be further transferred into H that heps to wash out more L, eading to a arger B L asymmetry. With two Higgs doubets, Φ and H, the induced B L asymmetry is at most one third (5/12 incuding the t- and b-yukawa couping effects) of the initia B + L asymmetry from GUT baryogenesis, which is arger than the asymmetry obtained in Ref. [15], where the type-i seesaw Yukawa couping is used to erase L and to produce a nonzero B L asymmetry. umericay, we have found that in order to generate a neutrino mass of m 2 atm(= 0.05 ev) and achieve YB L fina initia /YB+L O(10 2 ), the mass of the right-handed neutrino M has to be roughy arger than GeV for TeV DM. If the initia B + L asymmetry is sizabe ( 10 6 ) and the B + L injection scae is M /10, M can be as ow as GeV to accommodate both the baryon asymmetry and the neutrino mass, since L washouts can sti be efficient enough to create a non-vanishing B L asymmetry before the sphaerons competey destroy the B + L asymmetry. We have made sure that with Φ heavier than 3 TeV one can reproduce the observed reic abundance which requires an O(1) couping λ 3 and at the same time avoid the XEO1T direct search bounds. In this case Φ fas out of equiibrium before the eectroweak symmetry breaking and the asymmetry stored in Φ wi convert back into H. That sighty increases the fina baryon asymmetry. To summarize, we have estabished an intriguing correation among GUT baryogenesis, DM phenomenoogy and neutrino mass mechanism, where the existence of DM revives GUT baryogenesis and induces the radiative neutrino mass. 13

14 Acknowedgments The authors woud ike to thank Aveino Vicente for pointing out a missing factor of 1/2 in the neutrino mass matrix. This work is supported by Deutsche Forschungsgemeinschaft (DFG) Grant o. PA 803/10-1. S. Zeißner was supported by a feowship of Deutsche Studienstiftung. W. C. Huang is aso supported by Danish Counci for Independent Research Grant DFF The CP3-Origins centre is partiay funded by the Danish ationa Research Foundation, grant number DRF90. [1] H. Georgi and S. L. Gashow. Unity of A Eementary Partice Forces. Phys. Rev. Lett., 32: , [2] Motohiko Yoshimura. Unified Gauge Theories and the Baryon umber of the Universe. Phys. Rev. Lett., 41: , [Erratum: Phys. Rev. Lett.42,746(1979)]. [3] D. Toussaint, S. B. Treiman, Frank Wiczek, and A. Zee. Matter - Antimatter Accounting, Thermodynamics, and Back Hoe Radiation. Phys. Rev., D19: , [4] Steven Weinberg. Cosmoogica Production of Baryons. Phys. Rev. Lett., 42: , [5] Stephen M. Barr, Gino Segre, and H. Arthur Wedon. The Magnitude of the Cosmoogica Baryon Asymmetry. Phys. Rev., D20:2494, [6] Frans R. Kinkhamer and. S. Manton. A Sadde Point Soution in the Weinberg-Saam Theory. Phys. Rev., D30:2212, [7] Peter Brockway Arnod and Larry D. McLerran. Sphaerons, Sma Fuctuations and Baryon umber Vioation in Eectroweak Theory. Phys. Rev., D36:581, [8] Peter Brockway Arnod and Larry D. McLerran. The Sphaeron Strikes Back. Phys. Rev., D37:1020, [9] G. D. Coughan, Graham G. Ross, R. Homan, Pierre Ramond, M. Ruiz-Ataba, and J. W. F. Vae. Baryogenesis, Proton Decay and Fermion Masses in Supergravity GUTs. Phys. Lett., B158: , [10] K. S. Babu and R.. Mohapatra. Predictive neutrino spectrum in minima SO(10) grand unification. Phys. Rev. Lett., 70: , 1993, hep-ph/ [11] Bjorn Garbrecht, Tomisav Prokopec, and Michae G. Schmidt. SO(10)-GUT coherent baryogenesis. uc. Phys., B736: , 2006, hep-ph/ [12] Yoav Achiman. Spontaneous CP vioation in SUSY SO(10). Phys. Lett., B653: , 2007, hep-ph/ [13] K. S. Babu and R.. Mohapatra. B-L Vioating uceon Decay and GUT Scae Baryogenesis in SO(10). Phys. Rev., D86:035018, 2012, [14] K. S. Babu and R.. Mohapatra. Couping Unification, GUT-Scae Baryogenesis and eutron- Antineutron Osciation in SO(10). Phys. Lett., B715: , 2012, [15] M. Fukugita and T. Yanagida. Resurrection of grand unified theory baryogenesis. Phys. Rev. Lett., 89:131602, 2002, hep-ph/

15 [16] Wei-Chih Huang, Heinrich Päs, and Sinan Zeissner. eutrino assisted GUT baryogenesis - revisited. Phys. Rev., D97(5):055040, 2018, [17] Ernest Ma. Verifiabe radiative seesaw mechanism of neutrino mass and dark matter. Phys. Rev., D73:077301, 2006, hep-ph/ [18] Ernest Ma. Derivation of Dark Matter Parity from Lepton Parity. Phys. Rev. Lett., 115(1):011801, 2015, [19] Jeffrey A. Harvey and Michae S. Turner. Cosmoogica baryon and epton number in the presence of eectroweak fermion number vioation. Phys. Rev., D42: , [20] Ernest Ma. Common origin of neutrino mass, dark matter, and baryogenesis. Mod. Phys. Lett., A21: , 2006, hep-ph/ [21] Shoichi Kashiwase and Daijiro Suematsu. Baryon number asymmetry and dark matter in the neutrino mass mode with an inert doubet. Phys. Rev., D86:053001, 2012, [22] Shoichi Kashiwase and Daijiro Suematsu. Leptogenesis and dark matter detection in a TeV scae neutrino mass mode with inverted mass hierarchy. Eur. Phys. J., C73:2484, 2013, [23] J. Racker. Mass bounds for baryogenesis from partice decays and the inert doubet mode. JCAP, 1403:025, 2014, [24] Jackson D. Carke, Robert Foot, and Raymond R. Vokas. atura eptogenesis and neutrino masses with two Higgs doubets. Phys. Rev., D92(3):033006, 2015, [25] Thomas Huge, Moritz Patscher, and Kai Schmitz. Low-Scae Leptogenesis in the Scotogenic eutrino Mass Mode. 2018, [26] Sven Baumhozer, Vedran Brdar, and Pedro Schwaer. The ew νmsm : Radiative eutrino Masses, kev-scae Dark Matter and Viabe Leptogenesis with sub-tev ew Physics. 2018, [27] Diego Aristizaba Sierra, Chee Sheng Fong, Enrico ardi, and Eduardo Peinado. Coistered Baryogenesis. JCAP, 1402:013, 2014, [28] Kim Griest and Marc Kamionkowski. Unitarity Limits on the Mass and Radius of Dark Matter Partices. Phys. Rev. Lett., 64:615, [29] iendra G. Deshpande and Ernest Ma. Pattern of Symmetry Breaking with Two Higgs Doubets. Phys.Rev., D18:2574, [30] Aexander Mere and Moritz Patscher. Parity Probem of the Scotogenic eutrino Mode. Phys. Rev., D92(9):095002, 2015, [31] Aveino Vicente. Computer toos in partice physics. 2015, [32] G. F. Giudice, A. otari, M. Raida, A. Riotto, and A. Strumia. Towards a compete theory of therma eptogenesis in the SM and MSSM. uc. Phys., B685:89 149, 2004, hep-ph/ [33] Guy D. Moore. Do we understand the sphaeron rate? In Strong and eectroweak matter. Proceedings, Meeting, SEWM 2000, Marseie, France, June 13-17, 2000, pages 82 94, 2000, hep-ph/ [34] Michea D Onofrio, Kari Rummukainen, and Anders Tranberg. Sphaeron Rate in the Minima Standard Mode. Phys. Rev. Lett., 113(14):141602, 2014, [35] S. Yu. Khebnikov and M. E. Shaposhnikov. The Statistica Theory of Anomaous Fermion umber onconservation. uc. Phys., B308: ,

16 [36] P. A. R. Ade et a. Panck 2015 resuts. XIII. Cosmoogica parameters. Astron. Astrophys., 594:A13, 2016, [37] W. Buchmuer, R. D. Peccei, and T. Yanagida. Leptogenesis as the origin of matter. Ann. Rev. uc. Part. Sci., 55: , 2005, hep-ph/ [38] Abdessam Arhrib, Yue-Lin Sming Tsai, Qiang Yuan, and Tzu-Chiang Yuan. An Updated Anaysis of Inert Higgs Doubet Mode in ight of the Recent Resuts from LUX, PLACK, AMS-02 and LHC. JCAP, 1406:030, 2014, [39] Agnieszka Inicka, Maria Krawczyk, and Tania Robens. Inert Doubet Mode in ight of LHC Run I and astrophysica data. Phys. Rev., D93(5):055026, 2016, [40] ajimuddin Khan and Subhendu Rakshit. Constraints on inert dark matter from the metastabiity of the eectroweak vacuum. Phys. Rev., D92:055006, 2015, [41] Aexander Beyaev, Giacomo Cacciapagia, Igor P. Ivanov, Feipe Rojas-Abatte, and Marc Thomas. Anatomy of the Inert Two Higgs Doubet Mode in the ight of the LHC and non- LHC Dark Matter Searches. Phys. Rev., D97(3):035011, 2018, [42] Benedikt Eiteneuer, Andreas Goudeis, and Jan Heisig. The inert doubet mode in the ight of Fermi-LAT gamma-ray data: a goba fit anaysis. Eur. Phys. J., C77(9):624, 2017, [43] Camio Garcia-Cey, Michae Gustafsson, and Aejandro Ibarra. Probing the Inert Doubet Dark Matter Mode with Cherenkov Teescopes. JCAP, 1602(02):043, 2016, [44] Gerard Jungman, Marc Kamionkowski, and Kim Griest. Supersymmetric dark matter. Phys. Rept., 267: , 1996, hep-ph/ [45] Gianfranco Bertone, Dan Hooper, and Joseph Sik. Partice dark matter: Evidence, candidates and constraints. Phys. Rept., 405: , 2005, hep-ph/ [46] E. Aprie et a. First Dark Matter Search Resuts from the XEO1T Experiment. Phys. Rev. Lett., 119(18):181301, 2017, [47] James M. Cine, Kimmo Kainuainen, Pat Scott, and Christoph Weniger. Update on scaar singet dark matter. Phys. Rev., D88:055025, 2013, [Erratum: Phys. Rev.D92,no.3,039906(2015)]. [48] Xiangyi Cui et a. Dark Matter Resuts From 54-Ton-Day Exposure of PandaX-II Experiment. Phys. Rev. Lett., 119(18):181302, 2017, [49] D. S. Akerib et a. Resuts from a search for dark matter in the compete LUX exposure. Phys. Rev. Lett., 118(2):021303, 2017, [50] E. Aprie et a. Dark Matter Search Resuts from a One Tonne Year Exposure of XEO1T. 2018,