arxiv: v1 [hep-ph] 12 Sep 2018

Transcription

1 WU-HE arxiv: v1 [hep-ph] 12 Sep 2018 Affleck-Dine baryogenesis in the SUSY DFSZ axion model without R-parity Kensuke Akita 1 and Hajime Otsuka 2 1 Department of hysics, Tokyo Institute of Technology, Tokyo , Japan 2 Department of hysics, Waseda University, Tokyo , Japan Abstract We investigate the baryon asymmetry in the SUSY DFSZ axion model without R-parity. It turns out that the R-parity violating bilinear terms economically explain the atmospheric mass-squared difference of neutrinos and the appropriate amount of baryon asymmetry through the Affleck-Dine mechanism. In this model, the axion is a promising candidate for the dark matter and the axion isocurvature perturbation is mildly suppressed due to the large field values of eccei-quinn fields. Remarkably, in some parameter regions explaining the baryon asymmetry and the axion dark matter abundance, the proton decay will be explored in future experiments. address: kensuke@th.phys.titech.ac.jp address: h.otsuka@aoni.waseda.jp

2 Contents 1 Introduction 1 2 Model and constraints from experiments Setup Constraints on the R-parity violating couplings roton decay Baryon washout Neutrino masses Affleck-Dine baryogenesis in the SUSY DFSZ axion model with the R-parity violating terms Affleck-Dine mechanism otential for the AD/Q fields Dynamics of the AD/Q fields during inflation Dynamics of the AD/Q fields after inflation: H > m φ Dynamics of the AD/Q fields at H m φ and baryon asymmetry Saxion decay in the SUSY DFSZ model Axion isocurvature perturbation Conclusion and discussion 21 A General coupling between AD/Q fields and the minimum during inflation 22 2

3 1 Introduction The Standard Model SM) of elementary particle physics is well consistent with the experimental data up to TeV scale. However, neutrino oscillations reported in the SuperKamiokande [1], strong C problem in the SM, the excess of baryon over anti-baryon in the current Universe, and dark matter experiments indicate new physics beyond the SM. To explain these phenomenological problems, we focus on the supersymmetry SUSY) as an extension of the SM. So far, several models are proposed to explain the non-vanishing neutrino masses as represented by the see-saw mechanism, introducing the heavy right-handed Majorana neutrinos in the SM [2, 3, 4, 5, 6]. Among supersymmetric models, the R-parity violating SUSY scenario is a simplest approach to explain the tiny neutrino masses [7]. 1 Because of the R-parity violation, the lepton number violating interactions generate the neutrino masses without introducing a new particle in the framework of SUSY. In addition, the R-parity violating interactions help us to avoid the cosmological gravitino and moduli problems such as the overabundance of the lightest supersymmetric particle generating from the gravitino and moduli decays [9, 10, 11]. However, the sizable R-parity violating interactions cause several cosmological and phenomenological problems such as the unobservable proton decay, undetectable collider signatures of supersymmetric particles and washing out the primordial baryon asymmetry. Thus it is required to explain the smallness of R-parity violating interactions. To explain the smallness of R-parity violating interactions, we focus on the SUSY DFSZ axion model [12]. In this model, the strong C problem in the quantum chromodynamicsqcd) can be solved using the eccei-quinn Q) mechanism [13]. The Nambu-Goldstone boson called axion [14, 15] appears through the spontaneous symmetry breaking of global U1) Q and dynamically cancels the C phases in QCD. Furthermore this axion is the promising candidate for the dark matter and its coherent oscillation explains the current dark matter abundance [16]. Note that the SUSY DFSZ axion model controls the size of the µ-term and the R-parity violating couplings by U1) Q symmetry. It is then important to explore whether or not enough baryon asymmetry is realized in such an extension of the SM. In an inflationary era, the primordial baryon asymmetry is diluted away by the accelerated expansion of the Universe [17, 18]. The baryon asymmetry, in particular the baryons to photons ratio η required by the big-bang nucleosynthesis and cosmic microwave background [19], should be created after inflation. To obtain the sizable baryon asymmetry after inflation, we study the Affleck-Dine AD) mechanism [20, 21] in the SUSY DFSZ axion model without R-parity. Since Q fields couple to the baryon/lepton number violating terms to control the size of the R-parity violating couplings in this setup, it is nontrivial that the AD mechanism produces the appropriate amount of baryon asymmetry. The aim of this work is to reveal the origin of baryon asymmetry in the SUSY DFSZ axion model without R-parity in comparison with usual supersymmetric and R-symmetric axion models. In contrast to the R-parity violating AD mechanism proposed in Ref. [22], the R- parity violating terms in our model couple to Q fields. The dynamics of Q fields affect the magnitude of R-parity violation and the amount of baryon asymmetry. In this regard, our model is different from Ref. [22]. We find that there exist some parameter regions, explaining the 1 For more details, see, e.g., Ref. [8]. 1

4 current baryon asymmetry, axion dark matter abundance and the smallness of µ- and R-parity violating terms. It is also consistent with the unobservable proton decay, undetectable collider signatures of supersymmetric particles, the atmospheric mass-squared difference of neutrinos and washout of the baryon asymmetry. Remarkably, in this model, the axion isocurvature perturbation is mildly suppressed because of the enhancement of the Q breaking scale in the early Universe. Furthermore, in some parameter regions explaining baryon asymmetry and axion dark matter abundance, the proton decay will be explored in future experiments such as the HyperKamiokande [23]. In particular, this model is economical in the point of view of explaining the atmospheric neutrino mass data, the baryon asymmetry and the strong C problem in supersymmetric models. This is because we do not introduce a new field to explain the atmospheric neutrino mass data and the baryon asymmetry in the SUSY DFSZ axion model and we do not impose the discrete symmetry called R-parity which is usually imposed in supersymmetric models. This paper is organized as follows. In Sec. 2, we consider the SUSY DFSZ axion model without R-parity and show constraints on the R-parity violating couplings. In Sec. 3 we investigate the AD leptogenesis in more details. Then we study the dynamics of the AD and Q fields and estimate the baryon asymmetry, taking account of the dilution of saxion decay. The axion isocurvature perturbation is then mildly suppressed. Finally, we conclude in Sec Model and constraints from experiments 2.1 Setup In this paper, we consider the SUSY DFSZ axion model where the global U1) Q symmetry and the Q fields S 0, S 1, S 2 are introduced adding to the Minimal Supersymmetric SM MSSM) 2. The gauge symmetry of the MSSM does not prohibit both of the baryon and lepton number violating couplings in the superpotential. The so-called R-parity is often assumed in the MSSM to control these couplings, since the proton lifetime easily becomes shorter than the observational bounds. In the SUSY DFSZ model, the U1) Q symmetry plays this role instead of the conventional R-parity. It will be shown that the charge assignment of the U1) Q explains not only the smallness of the baryon and lepton number violating couplings, but also the size of the µ-parameter. We assume the following superpotential, W = W MSSM + W Rp / + W Q, 1) W MSSM = y u iju i Q j H u y d ijd i Q j H d y e ije i L j H d + y 0S 2 1 M H u H d, 2) / = y is1 3 L i H u + γ ijks 1 L i L j e k + γ ijk S 1 L i Q j d k + γ M W Rp M 2 2 See, as a review, Ref. [24]. M ijk S3 1 M 3 u i d j d k, 3) W Q = κs 0 S 1 S 2 f 2 ), 4) 2

5 S 0 S 1 S 2 H u H d u i d i Q i e i L i Q L Table 1: Q charges and lepton numbers of fields. where i, j, k = 1, 2, 3 are family indices. The Higgs fields, quarks and leptons are the usual representations of the MSSM gauge group, and the Q fields are singlet under the MSSM gauge group. y 0, y i, γ ijk, γ ijk, γ ijk, κ are dimensionless parameters which are typically of O1). The Q charges and also lepton numbers of each field are listed in Table 1. f in W Q is the present Q breaking scale and M GeV is the reduced lanck mass. When we arrange the Q fields to S 1 = Se A/f, S 2 = Se A/f, A behaves as the axion superfield. W Q induces a spontaneous symmetry breaking of U1) Q at the minimum S f. The present effective superpotential after the symmetry breaking is given by W eff =W Yukawa + µ 0 e 2A/f H u H d + µ i e 3A/f L i H u + λ ijk e A/f L i L j e k + λ ijke A/f L i Q j d k + λ ijke 3A/f u i d j d k, 5) where W Yukawa denotes the Yukawa coupling terms in Eq. 2) and µ 0 = y 0f 2, µ i = y if 3, M λ ijk = γ ijk f M M 2 ), λ ijk = γ ijk f M ) ) 3 f, λ ijk = γ ijk 6) M are defined. For M = GeV, f = GeV, y 0 = y i = 1, γ ijk = γ ijk = 1 2 we obtain and γ ijk = 1, µ 0 = GeV, µ i = GeV, λ ijk = , λ ijk = , λ ijk = ) Thus the smallness of the µ-parameter is explained by the Kim-Nilles mechanism [25], and at the same time the R-parity violating couplings are suppressed by the Froggatt-Nielsen mechanism [26] by the U1) Q assignment. The baryon number violating coupling constants λ ijk are highly suppressed, whereas the lepton number is sizably violated. These facts lead the stability under the single nucleon decay [27] and suppress effects from the washing out the baryon asymmetry [28, 29, 30, 31]. It may be possible that one imposes flavor dependent assignment of the U1) Q, which is related to the Flaxion models proposed by Refs. [32, 33, 34]. 2.2 Constraints on the R-parity violating couplings We shortly review experimental constraints on the R-parity violating terms. In the following, we only enumerate some severe bounds on the R-parity violating couplings. A more comprehensive review is written in Ref. [8]. 3

6 2.2.1 roton decay The observations of single nucleon decay give an important constraint on the trilinear R-parity violating terms [27]. A co-existence of the baryon and the lepton number violating operators induces a decay p π 0 l + mediated by d R squark in s-channel which has not been detected so far, thus giving the upper bounds on the trilinear R-parity violating couplings [35, 36], ) 2 m λ imkλ 11k < O1) d, 8) 5 TeV where i, k = 1, 2, 3, m = 1, 2 and m d is the typical down-type squark mass. The upper bounds are very severe, but the U1) Q -charge assignment leads λ λ 10 26, so that the proton lifetime is longer than the bound Baryon washout Other upper bounds on the trilinear R-parity violating terms come from the observation of the baryon asymmetry. If the baryon asymmetry is produced before the electroweak phase transition and the sphaleron process [37] is in the thermal equilibrium, the R-parity violating terms would erase the existing baryon asymmetry [28, 29, 30, 31]. To avoid this effects, the trilinear R-parity violating couplings are upper bounded by where m f is the typical mass of sfermions Neutrino masses ) m 1/2 λ, λ, λ < f, 9) 1 TeV We next focus on the bilinear R-parity violating couplings µ i. Constraints on the R-parity violating couplings come from the cosmological observations on neutrino masses because the bilinear R-parity violation gives rise to the neutrino mass. When we choose the basis in which the vacuum expectation values VEVs) of sneutrinos vanish and the Yukawa couplings of charged lepton are diagonal, the effective neutrino mass matrix at tree level is generated by the bilinear R-parity violating terms [38], M ν tree m νtree Σ 3 i=1 µ2 i µ 2 1 µ 1 µ 2 µ 1 µ 3 µ 1 µ 2 µ 2 2 µ 2 µ 3 µ 1 µ 3 µ 3 µ 2 µ ) The size of neutrino mass m νtree is given by m νtree m2 zcos 2 θ W M 1 + sin 2 θ W M 2 ) M 1 M tan 2 β) tan 2 ξ, 11) 4

7 where m Z, M 1, M 2 are the Z boson mass, gaugino masses respectively, and θ W is the Weinberg angle. tan β = v u /v d is the ratio between VEVs of two Higgs fields v u, v d and tan 2 ξ µ2 1 + µ µ ) µ 2 0 Since the rank of mass matrix is one, one of neutrinos acquires a mass at tree level. Although the other two neutrinos acquire masses by quantum corrections [39, 40], the detailed analysis of the other neutrino masses are beyond the scope of this paper. Indeed, one of the neutrino mass-squared differences can be explained by the bilinear R-parity violating terms, observed by the atmospheric neutrino observation [41] m 2 A ev. 13) If we explain the atmospheric neutrino observation without tuning of dimensionless parameters, the size of neutrino mass m νtree is constrained as m νtree ev which leads to the constraint on the bilinear R-parity violating couplings, ) µ 2 i ev1 + tan 2 M2 β) µ 2 1 TeV 0. 14) i In addition, the cosmological bound on neutrino masses i m ν i 0.11 ev m νtree 0.11 ev) [42] leads to another constraint on the bilinear R-parity violating couplings, ) µ 2 i tan 2 M2 β) µ 2 1 TeV 0. 15) i In our model, the R-parity violating couplings 7) avoid the constraints of Eqs. 8), 9), 14) and 15). Furthermore, the bilinear R-parity violating couplings can explain atmospheric mass-squared difference of neutrinos in some region of dimensionless parameters. Furthermore, in the low-scale SUSY scenario, we may observe the single nucleon decay for f GeV in future experiments such as HyperKamiokande [23]. 5

8 3 Affleck-Dine baryogenesis in the SUSY DFSZ axion model with the R-parity violating terms In this section, we study the Affleck-Dine mechanism exploiting the R-parity violating terms in the SUSY DFSZ axion model. A notable point of our scenario is that the AD field couples to the Q field and their dynamics are affected with each other. This behavior makes a different amount of the baryon asymmetry from the conventional one. In the following, we will take parameters which are consistent with the allowed region to explain the atmospheric neutrino mass-squared difference in Sec Affleck-Dine mechanism Here and in the following, we consider the Affleck-Dine leptogenesis exploiting one of the L i H u - directions of the scalar potential 3. In the conventional AD mechanism, the quadratic L i H u operators can not produce the baryon asymmetry [21]. However, in this model, the L i H u - direction coupled to the Q field is not quadratic and will produce the baryon asymmetry. We consider the LH u D-flat direction as the so-called AD field, namely L = 1 2 φ 0 ), H u = φ ). 16) In the next subsection, the potential for the AD field and Q fields is derived. The dynamics of the AD/Q fields during inflation is studied in Sec After inflation, we show their dynamics in an era H > m φ H m φ ) in Sec ). The resultant baryon asymmetry is extracted in Sec otential for the AD/Q fields We assume that the other D-flat directions have positive mass terms, so that we ignore effects from these directions in the following discussion. Note that some other D-flat directions coupled with the Q fields will be available for the Affleck-Dine mechanism. For more details, see Appendix A, where we discuss the AD mechanism with general couplings between Q field S 1 and AD field φ. The generalized setup can be applied to the R-parity conserving case. The potential of the AD field and Q fields depends on SUSY breaking scenarios and it is important for the AD mechanism how these fields couple to the inflaton. In this paper, we assume the gravity-or anomaly-)mediated SUSY breaking scenario, together with the F -term inflation. In the gauge-mediated SUSY breaking scenario, it is nontrivial that the AD mechanism works in this model. This is because the scalar potential is modified by the effect of the gauge-mediated supersymmetry breaking [43]. 3 If we impose the lepton number of S 1 as 1 3, the AD leptogenesis via L ih u cannot work. However, it may be possible that a certain baryon asymmetry is produced by other directions of the potential, e.g. LLe, based on a charge assignment different from the one in Table 1. The detailed calculation is one of the future works. 6

9 Let us consider the following superpotential, W Rp W = W inf I) + W Rp / S 1, φ) + W Q + W mix, / S 1, φ) = ys3 1φ 2, W Q = κs 0 S 1 S 2 f 2 ), 17) 2M 2 where W inf I) is the superpotential for the inflaton I, W Rp / φ) is originated from the first term in Eq.3) and W mix stands for possible mixing terms between the inflaton and the AD/Q fields. Here κ is a coupling constant. LH u -direction has a flavor dependence. Here and in what follows, we assume that the lepton asymmetry is produced by the flattest LH u -direction because the flattest direction produces the lepton asymmetry most efficiently. Then, y corresponds the smallest y i. The supergravity scalar potential is given by [44], [ V = e K/M 2 D a W )K ab D b W ) 3 W ], 2 18) M 2 with D a W = W Φ + K ) W, K ab 2 1 K =, 19) a Φ a M 2 Φ a Φ b where Φ a,b = I, φ, S 0, S 1, S 2 and K is the Kähler potential. We assume the following Kähler potential with non-minimal couplings, K = a Φ a Φ a + α φ φi I + M 2 β M 2 S 1S 1 I I + OM 3 ), 20) where coupling constants α, β are introduced. Here S 0 and S 2 have only the minimal Kähler potential. The inflaton scalar potential is given by V inf e K/M 2 F I K II F 3 W I M 2 inf I) ). 2 21) If α, β 1, the Hubble-induced mass terms are provided by the F -term of the inflaton, V Hubble c 0 H 2 S 0 2 c 1 H 2 S c 2 H 2 S 2 2 c 3 H 2 φ 2, 22) where c 0, c 1, c 2, c 3 are positive constants and H is the Hubble parameter. The signs of the mass terms indicate that S 0 and S 2 acquire large positive masses during inflation. There are also the soft SUSY-breaking mass terms, V soft = m 2 0 S m 2 1 S m 2 2 S m 2 φ φ 2. 23) 7

10 Next, let us consider the contributions from W Rp / S 1, φ), W Q and W mix. The F -term scalar potential is given by V F = κ 2 S 1 S 2 f κs 0S 2 3yS2 1φ 2 2+ κ 2 S 0 S y 2 S 1 6 φ 2. 24) 2M 2 In our study, we assume I M and the dynamics of the inflaton is basically separated from those of the AD/Q fields. Even in this case, the F -term of the inflaton F I may give significant effects to the dynamics of the AD/Q fields. Let us consider the following mixing superpotential W mix = α I W Rp / S 1, φ) + β I W Q, 25) M M where α, β are coupling constants. This superpotential induces the potential, ) F α W Rp / S 1, φ) + β W I Q + h.c., 26) M which are not suppressed even for I M. Therefore we expect that there are Hubbleinduced couplings in the scalar potential. In addition, there are soft SUSY-breaking terms, generating from the mixing superpotential, V A =a H H + a m m 3/2 )W Rp / S 1, φ) + b H H + b m m 3/2 )W Q + h.c. = a H H + a m m 3/2 ) ys3 1φ 2 + b H H + b m m 3/2 )κs 0 S 1 S 2 f 2 ) + h.c., 27) 2M 2 where a H, b H stand for the Hubble-induced couplings and a m, b m are soft SUSY breaking couplings. The size of the soft SUSY-breaking terms is represented by the gravitino mass m 3/2 as usual in the gravity mediation. Finally, the whole scalar potential is given by collecting the above contributions, Dynamics of the AD/Q fields during inflation M 4 V = V Hubble + V soft + V F + V A. 28) In this section, we show the dynamics of the AD/Q fields during inflation. We find a minimum of the scalar potential and consider the dynamics of the AD/Q fields around the minimum. The details of the calculation are shown in Appendix A. During inflation, H m 3/2, the soft SUSY-breaking terms can be neglected. We firstly focus on the phase dependent part of the potential. The fields can be decomposed to φ i = ˆφ i e iθ φ i, φi = S 0, S 1, S 2, φ, 29) 8

11 where ˆφ i and θ φi are real fields. The phase-dependent scalar potential is given by V phase = 2ˆκ 2 f 2 Ŝ 1 Ŝ 2 cos θ S1 + θ S2 ) 3ˆκŷ Ŝ0Ŝ2 ˆφ 1Ŝ2 2 cos θ S0 + θ S1 2θ S2 2θ φ + ζ) 2ˆκˆb H HŜ0)Ŝ1Ŝ2 cos θ S0 + θ S1 + θ S2 + η) + 2ˆκˆb H f 2 HŜ0) cos θ S0 + η) ŷâ H H Ŝ3 1 ˆφ 2 M 2 M 2 cos 3θ S1 + 2θ φ + ξ), 30) where the coupling constants with a hat stand for their absolute values and ζ Argκ y), η Argκb H ) and ξ Argya H ) are defined. The first line comes from the F -term potential V F and the last two lines come from the Hubble-induced A-terms in V A. If H Ŝ0 f 2 as shown in later, the second line can be neglected and the minimum of the phases are placed at θ S1 + θ S2 0, θ S0 + θ S2 2θ S1 2θ φ ζ, 3θ S1 + 2θ φ ξ. 31) Under this condition, the extremal condition for the radial directions V/ ˆφ i = 0 has a solution: Ŝ0 = ˆκ Ŝ2 3ŷ Ŝ1 2 ˆφ 2, ˆκ 2 Ŝ1 2 + ˆκ 2 Ŝ2 2 + c 0 H 2 2M 2 Ŝ2 f 2 Ŝ1, Ŝ1 k ˆφ, ˆφ 3kâH + 9k 6 â 2 H + 8c 19k 2 + 6k 4 ) 18k k 4 HM 2 ŷ ) 1 3, 32) where k is an O1) constant which depends on â H, c 1, c 3. At this minimum, our assumption prior to the estimation: H Ŝ0 f 2 H/ ˆφ f 2, 33) is satisfied and the estimation is self-consistent as long as H M. Since the masses of AD/Q fields m φi φ i = S 0, S 1, S 2, φ) at the extrema are estimated by using Eqs. 28) and 32), m S0 m S2 Ŝ1 m S1, m φ H, 34) S 0 and S 2 are expected to have large positive mass terms and to be fixed at the minimum during inflation. The curvature along the Ŝ1, ˆφ directions is determined by the mass matrix, 2 V 2 V 1 Ŝ1 Ŝ1 Ŝ1 ˆφ 2 2 V 2 V. 35) ˆφ Ŝ1 ˆφ ˆφ 9

12 Approximately, if a H c 1, c 3, the curvature is positive which is confirmed by numerical calculation. As long as this condition is satisfied, the AD field φ and the Q field S 1 have masses of Oâ H H). Therefore, all of the AD/Q fields will settle at the minimum during inflation. The axionic isocurvature perturbation [47, 48, 49, 50] and the baryonic isocurvature perturbation [51, 52, 53, 54, 55] will be so suppressed by the large VEVs that the model is consistent with the lanck observations [42]. The details are discussed in Sec Dynamics of the AD/Q fields after inflation: H > m φ After the end of inflation, the inflaton starts to oscillate around I min and the effect of Hubbleinduced A-term a H H on the dynamics of the AD/Q fields turns off. In this era, the Universe is dominated by the oscillation energy and the Hubble parameter evolves as H = 2/3t). The Kähler potential and superpotential are approximately K inf = I 2 + = I minδi + I min δi + δi 2 + OδI 3 ), W = 1 2 M infi I min ) = 1 2 M infδi 2 + OδI 3 ), 36) where δi = I I min and M inf is the inflaton mass. Since F I = M inf δi and the Hubble-induced A-term a H H is originated from F I, a H diminishes rapidly after the F -term inflation [55, 56] as long as the inflaton oscillates in the period M 1 inf which is shorter than the Hubble time H 1. The suppression factor is estimated to [56] H H inf, 37) where H inf is the Hubble parameter during inflation. The Hubble-induced A-term with a H H after inflation is given by ya H H 2 S 3 1φ 2 2H inf M 2 + h.c.. 38) Eventually, this Hubble-induced A-term diminishes as time evolves. The point Eqs. 32) is no longer the minimum of the potential and the value of this extrema changes. Then, the AD/Q fields start to roll down. Since Eqs. 28),32) and 34) show S 0 and S 2 are heavier than S 1 and φ until Ŝ1 is larger than f, we assume that S 0 and S 2 are fixed at the minimum in this epoch. The dynamics of S 1, φ obey the following equations of motion, d 2 S 1 dt t ds 1 dt + V S 1 = 0, d 2 φ dt + 2 dφ 2 t dt + V = 0. 39) φ The soft SUSY-breaking terms can also be neglected until H m φ m 3/2. 10

13 as We numerically solve the equations. Let us introduce new parameters z, s 1 and χ defined z = log H inf t, S 1 = s 1 2Hinf M 2 3ŷ ) 1 ) e z 3 2Hinf M 2 1, φ = χ e z 3, 40) 3ŷ from which Eqs. 39) become 2 s 1 z + 1 s z 4c χ z + 1 χ 2 3 z 4c s s 1 2 χ 4 s s 1 4 χ 2 s 1 4a H 9 e z s 2 1 χ 2 = 0, χ + 2 s 1 4 χ 2 χ s 1 6 χ 8a H 27 e z s 3 1 χ = 0. 41) The time variations of s 1 with s 1R = Res 1 ), s 1I = Ims 1 )) and χ with χ R = Reχ), χ I = Imχ)) are shown in the left panel and the right panel of Fig. 1, respectively. The redblue dash) lines represent the realimaginary) parts of scalar fields. For concreteness, ŷ = 1, â H = 3 5, c 1 = 1/4, c 3 = 1/5, ξ = 0 are assumed. Then, k is estimated as k for â 2 H c 1, c 3. Since the variations of field values depend on the Hubble parameter during inflation H inf, we assume that the Hubble parameter during inflation is fixed at H inf = 10 9 GeV which is consistent with the isocurvature perturbation discussed in Sec The initial conditions at inflation end t = t end = H 1 inf z = log2/3)) are θ s1 z=log2/3) = π 9, θ φ z=log2/3) = π 6 42) and s 1 z = 0, z=log2/3) χ z = 0. 43) z=log2/3) In addition, the initial conditions of s 1, χ are extracted from Eqs. 32), s 1 z=log2/3) = k χ z=log2/3), 3kâH + 9k χ z=log2/3) = â 2 H + 8c ) 19k 2 + 6k 4 1 ) 3. 18k k 4 44) 3 Hereafter we take k = 2 The trajectories start from t = H 1 inf z = log2/3)) and evolve to t = m 1 φ 1 TeV 1 z 14). Finally, s 1, χ at H m φ 1 TeVz 14) have the following values, s 1R 14.9, s 1I 5.42, χ R , χ I , z s 1R 5.54, z s 1I 2.02, z χ R 1.92, z χ I ) These values give initial values for the dynamics at H m φ which is analyzed in next subsection. The numerical calculation indicates that s 1 increases by O10) while χ diminishes by O10 3 ) as the time evolution for c 1 > c 3. This result will not depend on the parameters and 11

14 the initial values except for c 1, c 3 sensitively. However, the variations due to the dynamics will decrease for smaller H inf and c 1, c 3. Note that S 1 does not dominate the energy density of the Universe. Although s 1 increases during the time m φ H, the field value of Ŝ1 is still smaller than the reduced lanck mass, Thus, the inflaton dominates the Universe during m φ H, Ŝ 1 M. 46) ρ S1 H 2 Ŝ 2 1 ρ inf 3H 2 M 2. 47) In addition, the dynamics of s 1 and χ are characteristic for the c 1 = c 3 case. Although the point Eqs. 32) becomes a saddle point after inflation, s 1 and χ oscillate around the ridges of the saddle point. Due to the large field values of Ŝ1 and φ, the quantum fluctuations of Ŝ1 and φ will not affect the dynamics of s 1 and χ at the time m φ H. Fig. 2 shows the time variations of s 1 and χ for c 1 = c 3 = 1. The values of the other parameters and the initial conditions are the same as in Fig 1. From this perspective, we confirm that both s 1 and χ stay around O1) values. In the c 1 < c 3 case, we find that s 1 decreases while χ increases in typical parameter spaces. Since we must consider the dynamics of S 2 for small s 1, the dynamics of AD/Q fields after inflation are complicated and we leave it as a future work. In the next section, we therefore concentrate on two cases c 1 > c 3 and c 1 = c 3 by estimating the dynamics of S 1, φ for H m φ s 1 R, s 1 I 5 χ R, χ I z z Figure 1: The time variations of s 1, χ for c 1 > c 3 during m φ < H. In the left panel, we draw the time variations of s 1R, s 1I ), whereas, in the right panel, we draw the time variations of χ R, χ I ). The redblue dash) lines represent the realimaginary) parts of scalar fields. We set the parameters as ŷ = 1, â H = 5, c 1 = 1/4, c 3 = 1/5 and ξ = 0 and the initial conditions as θ s1 t=tend = π, θ 9 φ t=tend = π and ds 1 6 dt t=t end = 0, dχ dt t=t end = 0. We draw these time variations from H inf 10 9 GeVz log2/3)) to H m φ 1 TeVz 14). 12

15 s 1 R, s 1 I χ R, χ I z z Figure 2: The time variations of s 1, χ for c 1 = c 3 during m φ < H. In the left panel, we draw the time variations of s 1R, s 1I ), whereas, in the right panel, we draw the time variations of χ R, χ I ). The redblue dash) lines represent the realimaginary) parts of scalar fields. We set c 1 = c 3 = 1 and ξ = 0. The values of the other parameters and the initial conditions are the same as in Fig. 1. We draw these time variations from H inf 10 9 GeVz log2/3)) to H m φ 1 TeVz 14) Dynamics of the AD/Q fields at H m φ and baryon asymmetry In this subsection, we consider the dynamics of the AD/Q fields at H m 3/2 and finally estimate the amount of baryon asymmetry. At this epoch, the soft supersymmetry breaking effect becomes important. The AD/Q fields eventually start to oscillate around the minimum: Ŝ0 = ˆκ Ŝ2 ˆκ 2 Ŝ1 2 + ˆκ 2 Ŝ2 2 + c 0 H 2 3ŷ Ŝ1 2 ˆφ 2 2M 2 = 0, Ŝ1 f, Ŝ2 f, ˆφ = 0. 48) Then, S 0 and the scalar field orthogonal to the flat direction S 1 S 2 = f 2 have masses of order of f. These masses are heavier than the scale of soft mass Om 3/2 ) which is the same order of m φ. Thus, following the previous section, we can set S 1 S 2 = f 2 and Ŝ0 as Eqs. 32). In the following, we consider only the dynamics of S 1 and φ in this epoch. Since the masses of S 1 and φ are positive, S 1 and φ will go to the minimum Eqs. 48)). Because S 1 and φ will be enough damped, the potential of the AD field will be approximated as a quadratic one which does not depend on the Q fields. Thus the amount of baryon asymmetry will be conserved after S 1 f because we neglect C-violating terms. The above statement is numerically confirmed in the following analysis. First of all, let us analytically consider the dynamics of S 1, φ and estimate the amount of baryon asymmetry. The lepton number density is given by n L = i 2 dφ dt φ φ dφ dt ) = φ 2 dθ φ dt. 49) 13

16 It obeys the following equation of motion using the one of φ, ) dn L V dt + 3Hn L = Im φ φ, 50) and in the integral form, we obtain t ) Rt) 3 n L t) = dt Rt ) 3 yam m 3/2 S Im 1φ 3 2 t inf M 2 tosc ) = dt Rt ) 3 yam m 3/2 S 3 Im 1φ 2 t ) + dt Rt ) 3 yam m 3/2 S Im 1φ 3 2, 51) t osc t inf M 2 where R is the scale factor of the Universe and t inf H 1 inf is the time at the end of inflation. φ starts to oscillate after the time t osc defined as H osc = 2/3t osc ) m φ. The C-violating factor is defined as δ eff = sinξ + 3θ S1 + 2θ φ ). As explicitly checked in the numerical calculation, the second integration in the second line of Eq. 51) gives small effects to the baryon asymmetry since the sign phase factor δ eff changes rapidly after the AD/Q fields φ, S 1 start to oscillate. As a result, the lepton number Eq. 51) will be fixed at t = t osc and the lepton number density at t = t osc is estimated as where ɛ is defined as n L t osc ) 1 2 ɛâ mm 3/2 δ eff m3/2 M 2 ŷ ɛ = Ŝ1t osc ) 3 ˆφtosc ) 2 m3/2 M 2 ŷ M 2 ) 2 3, 52) ) ) To check the above statements, we numerically solve the equations of motion of S 1, φ and estimate the trajectories and the lepton asymmetry for H m φ. Here, we assume that the inflaton still dominates the Universe, namely matter-dominated Universe. Reparametrizing again as m3/2 M 2 S 1 S 1 ŷ m3/2 M 2 φ φ ŷ ) 1 3, ) 1 3, 54) 14

17 Eqs. 39) approximately become d 2 S 1 dt t ds 1 dt m2 3/2Ŝ1 ˆφ4 S 1 + 3m 2 3/2Ŝ4 ˆφ 1 2 S 1 3a m m 2 3/2S 2 1 φ 2 3/2 f 4 ˆφ4 S1 5 9m2/3 3/2 f 4 ˆφ4 + 9m2/3 2M 8/3 Ŝ1 6 d 2 φ dt t 2M 8/3 Ŝ 4 1 S m 2 φ 4c ) 1 S 9t 2 1 m 2 φ + 4c 2 9t 2 dφ dt m23/2ŝ4 ˆφ 1 2 φ + m 2 6 3/2Ŝ1 φ am m 2 3/2S 3 1 φ + 3/2 f 4 ˆφ2 S 4 φ 9m2/3 3/2 f 4 ˆφ2 S1φ 4 + 2M 8/3 Ŝ1 4 9m2/3 2M 8/3 Ŝ1 4 f 8 ˆφ4 2m 2/3 3/2 M S 16/3 1 Ŝ1 6 ) m 4/3 3/2 f 4 M 8/3 Ŝ1 4 9f 8 ˆφ2 2m 2/3 3/2 M 16/3 Ŝ1 4 m 2 φ 4c 3 9t 2 S 1 = 0, φ ) φ = 0, 55) from which the oscillation time is extracted as t osc = 2c 1/2 3 /3m φ ). Here we set κ = y = 1 for simplicity. In the following analysis, we numerically solve Eqs. 55) from the time t osc. First of all, let us consider the case with c 1 > c 3. We solve Eqs. 55) and draw the trajectories of S 1R = ReS 1 ), S 1I = ImS 1 )) in the right panel of Fig.3 and the trajectories of φ R = Reφ), φ I = Reφ)) in the left panel of Fig. 3. Here we set the parameters as m φ = 1, ˆλ = 1, â m = c 2 = 1, c 1 = 1/4, c 3 = 1/5, argλa m ) = π/3, f = 10 9 in units of m 3/2 = 10 3 GeV = 1 and the initial conditions as in Eqs. 45). We take the initial time as t osc = 2c 1/2 3 /3m φ ). Fig. 3 shows that the AD field φ and the Q field S 1 go to the minimum 48). Then, we can estimate ɛ numerically and the ratio of the numerical value in Eq. 51) to analytical value in Eq. 52) is shown in Fig. 4. It turns out that the numerical value coincides with the analytical one in Eq. 52) with δ eff = 1 after the AD/Q fields start to oscillate. However, we find the numerical value is a little larger than the analytical estimation and the ratio of numerical value to analytical one becomes O10) in other parameter regions. This is because the rotation of φ is complicated as in Fig 3 and then the lepton number Eq. 51) is not exactly fixed at t osc. The mass of θ φ is heavier than the mass of θ S1 because in Fig. 1, the field value of ˆφ decreases while Ŝ 1 increases and the field value of ˆφ is smaller than the one of S 1. For this reason, θ φ moves whereas θ S1 does not as shown in Fig 3. Similarly, we numerically solve Eqs. 55) and draw the trajectories of S 1, φ in Fig. 5 for the c 1 = c 3 case. Here we set c 1 = c 3 = 1 and the values of the other parameters are the same as in Fig. 3. The initial conditions of S 1, φ are set by the values of s 1, χ at z = 14. Then, we can estimate ɛ numerically and the ratio of numerical value in Eq. 51) to analytical one in Eq. 52) is shown in Fig. 6. We also find that the numerical value coincides with the analytical one in Eq. 52) with δ eff = 1 after the AD/Q fields start to oscillate. In this case, the mass of θ φ is comparable to the mass of θ S1 because the field value of ˆφ is the same order as the one of Ŝ1. As the result, the phases of S 1 and φ rotate at the same time, drawn in Fig 5. 15

18 S 1 I S 1 R Φ 1 I Φ 1 R Figure 3: The dynamics of S 1, φ for c 1 > c 3 during H < m φ. In the left right) panel, we draw the trajectories of S 1R, S 1I ) φ R, φ I )) as a function of z. We set the parameters, ˆλ = 1, â m = c 2 = 1, c 1 = 1/4, c 3 = 1/5, f = 10 9 in units of m 3/2 = 10 3 GeV = 1, argλa m ) = π/3 and the initial conditions as in Eqs. 45). 1.5 R 3 n L /Rt osc ) 3 n L t osc ) m Φ t Figure 4: The time evolution of the ratio between the numerical value Rt) 3 nt) 51) and R 3 t osc )nt osc ) 52) for c 1 > c 3. The horizontal axis corresponds m φ t. Here we set ɛ = and δ eff = 1. The parameters are the same as in Fig

19 S 1 I S 1 R Φ 1 I Φ 1 R Figure 5: The dynamics of S 1, φ for c 1 = c 3 during H < m φ. In the left right) panel, we draw the trajectories of S 1R, S 1I ) φ R, φ I )) as a function of z. We set the parameters, ˆλ = â m = c 1 = c 2 = c 3 = 1, f = 10 9 in units of m 3/2 = 10 3 GeV = 1, argλa m ) = π/3 and the initial conditions are the values of s 1, χ at z = 14 in Fig. 2. R 3 n L /Rt osc ) 3 n L t osc ) m Φ t Figure 6: The time evolution of the ratio between the numerical value Rt) 3 nt) 51) and R 3 t osc )nt osc ) 52) for c 1 = c 3 = 1. The horizontal axis corresponds m φ t. Here we set ɛ = and δ eff = 1. The parameters are the same as in Fig

20 The ratio of the lepton number density n L to the entropy density s after the reheating t = t reh ) is n L s = 1 ) 3 Rtosc ) nt osc) st reh ) Rt reh ) = ɛâ ) mm 3/2 δ eff T reh mφ M 2 2 3, 56) 8M 2 m2 φ ŷ where T reh is the reheating temperature. Then, the AD field decays and its energy density converts into radiation [57]. Consequently, the sphaleron effect relates the lepton asymmetry to the baryon asymmetry [37, 58, 59], n B s 8 n L 23 s. 57) Thus, the baryon asymmetry is estimated as n B s = ɛâ mm 3/2 δ eff T reh 23M 2 m2 φ mφ M 2 ŷ ) 2 3 ) Treh GeV ŷ ) Treh GeV ŷ ) 1 ) 1 3 m 3/2 ) mφ m φ ) mφ 3 m 3/2 m φ 1 TeV 1 TeV ) 1 3 ) 1 3 for c 1 = 1, c 4 3 = 1, ɛ = , for c 1 = c 3 = 1, ɛ = where we set â m = 1 and δ eff 1. It is remarkable that the obtained baryon asymmetry is well consistent with its current observed value and the tiny neutrino mass, simultaneously. Indeed, the baryon asymmetry 58) is not so sensitive to ŷ which is severely constrained by the neutrino mass discussed in Sec Finally, we comment on the Q-ball problem [60]. If the potential of the AD field φ is flatter than the quadratic one, the AD fields form the nontopological solitons called Q-balls [61, 62]. We have investigated the AD leptogenesis via LH u flat direction in the gravity-mediated SUSY breaking scenario. Then, at one-loop level, the soft mass term of the AD field φ has no contribution from gluino loops which is mainly responsible to make the potential of the AD field φ flatter [63, 64]. In addition, a contribution from the top Yukawa coupling changes the potential steeper. Thus, in our discussing gravity-mediation scenario, Q-balls are not formed in the LH u direction [64] and the baryon asymmetry is described by Eq. 58). 3.2 Saxion decay in the SUSY DFSZ model In this section, we discuss the dynamics of the saxion S 1 in more details. During inflation, the energy density of the saxion field should be smaller than that of the inflaton field, which constrains the Hubble scale during inflation, 58) ρ S1, ρ Φ H 2 infφ 2 < ρ inf 3H 2 infm 2. 59) 18

21 Now we have used the field values 32) during inflation and the masses of S 1, Φ are of OH inf ) as mentioned below Eq. 35). Thus, from Eqs. 32) and 59), the field values of S 1, φ during inflation have to satisfy S 1, φ < M. After inflation, the inflaton and saxion fields oscillate around their minimum and their energy densities decay in proportion to R 3. Consequently, the inflaton decays at the time t reh 1 ) 1/2 90 M Γ inf π 2 g T reh ) GeV 1 ) 1/2 ) GeV, 60) where g denotes the effective degrees of freedom and Γ inf is the total decay width of the inflaton. On the other hand, the total decay width of saxion field depends on the sparticle spectrum. For more details, see e.g. Ref. [65].) When m S1 > 2µ, the saxion mainly decays into Higgsino through the µ-term in Eq. 2). If such a decay is kinematically disallowed, the total decay width of saxion is dominated by the saxion decay into the C-even Higgs bosons h, H and gauge bosons W ± and Z. Note that saxion decays into axions are suppressed in our setup S 1 S 2 f, taking into account two Q fields S 1,2 and singlet field S 0 with U1) Q charges in Table 1 [66]. 4 It then allows us to avoid the dark radiation problem from the saxion decay. As a result, the total decay width of the saxion 5 ) 2 Γ S 1 µ 1) ms1 m 4π f S1 > 2µ) tot 61) 7 2π µ 4 f 2 m S1 m S1 < 2µ) g T reh determines the decay temperature of the saxion ) 1/ µ GeV T S 1) g dec 10 4 GeV ) 1/4 ) µ 2 GeV g 10 4 GeV ) ) GeV m S1 f GeV f 1/ GeV) m 4 S1 > 2µ) ) ms1 1/2 10 GeV) m 3 S1 < 2µ), 62) and the decay time of the saxion t S 1) dec 1 Γ S 1) tot 90 π 2 g ) 1/2 M T S 1) dec ) GeV 1 ) ) 1/ GeV. 63) g T S 1) dec Let us examine whether or not the saxion decay dilutes the lepton asymmetry via the entropy production from the saxion decay. The entropy dilution factor is determined by the ratio of the saxion decay temperature T S 1) dec and the saxion-radiation equality temperature T eq, that is T eq /T S 1) dec. When the energy density of saxion is equal to that of radiation, T eq is given by T eq 1 6 T reh Ŝ0 1 M ) T reh Hinf M ) 2/3, 64) 4 For the dark radiation constraints in models with multiple Q multiplets, we refer to e.g., Refs. [65, 67]. 5 Now we consider the µ > m A where m A is the mass of C-odd Higgs boson. 19

22 where Ŝ0 1 is the VEV of Ŝ1 during inflation. In our setup discussed so far, the amplitude of saxion Ŝ0 1 H inf M 2 )1/3 is of O10 3 M ) with H inf = 10 9 GeV and then T eq is not larger than T S 1) dec unless T reh > GeV. Hence, there is no entropy dilution. Finally, we comment on the lepton asymmetry generating from the saxion decay. The saxion field S 1 would be identified with the right-handed Majorana neutrino as seen in the superpotential 3). The leptogenesis scenario decaying from the Majorana neutrino have been discussed in the thermal [70] and non-thermal epoch [71]. Since in our case the saxion oscillates around its minimum soon after inflation, it is possible to generate the lepton asymmetry through the coupling Such a lepton asymmetry is determined by ) 2 S1 W 3 S 1 LH u. 65) n S 1) L s M ɛ n S 1 s, 66) where ɛ involves the C asymmetry factor determined by the saxion decay into leptons and Higgs bosons. At the reheating era, the number density of saxion n S1 becomes n S1 s n S 1 t S1) t=treh st reh ) osc ) ρt reh ) ρt S 1) osc ) 3 4 ) T reh m S1 t S 2 1) osc ) Ŝ Treh Hinf 2 M 10 7 GeV ) ) /3 GeV, where we use m S1 t S 1) osc ) H inf and Ŝ M. Since the lepton asymmetry n S 1) L /s is further suppressed by the factor ɛ which is proportional to the effective Yukawa coupling S 1 /M ) 2 in Eq. 65), the saxion produced lepton asymmetry is negligible even when T reh > 10 7 GeV. 3.3 Axion isocurvature perturbation H inf 67) In this model, the massless Nambu-Goldstone boson called axion exists through the spontaneous symmetry breaking of the U1) Q. This axion is the candidate for the dark matter and the present axion energy density is given by [69], ) 1.19 Ω a h θa 2 f a, 68) GeV where θ a is the misalignment angle of the axion. f a is the axion decay constant depending on the domain wall number N DW which is 6 for the DFSZ axion model [12], 2f f a =. 69) N DW 20

23 For f = GeV and θ a = 1.9, the axion energy density is coincident with the dark matter energy density Ω CDM However, such a massless boson is problematic in the early universe. In our model, the U1) Q symmetry is spontaneously broken during inflation and is not recovered after inflation. Because of U1) Q symmetry breaking, domain wall problem [45, 46] is not occured. The Q field S 1 gets the large VEV S 1 φ H inf M 2 )1/3 during inflation. It can suppress the axion isocurvature perturbation [47, 48, 49, 50]. The axion almost consists of the linear combination of θ S1 and θ φ for S 1, φ f during inflation, so the axion a takes the form, a 1 4 Ŝ ˆφ 2 2 Ŝ1 2 θ S1 3 ˆφ 2 θ φ ). 70) The Q breaking scale v during inflation is v max[ Ŝ1, ˆφ ] H inf M 2 ) 1/3. 71) The power spectrum of CDM isocurvature perturbation iso is ) 2 iso r 2 Hinf, 72) πvθ a where r is the ratio of the present axion energy density to the matter energy density, r = Ω a h 2 /Ω m h 2. The lanck constraint on the uncorrelated isocurvature perturbation [42] becomes H inf GeVθ 1 a ) GeV. 73) f a Thus the axion isocurvature perturbation is mildly suppressed due to the large v. analysis, we consider H inf 10 9 GeV to avoid this constraint. In our 4 Conclusion and discussion In this paper we have investigated the baryon asymmetry in the SUSY DFSZ axion model without R-parity. Such R-parity violating interactions are motivated not only by explaining the tiny neutrino masses, but also by avoiding the cosmological gravitino and moduli problems. In this model, the Affleck-Dine mechanism can work out via the coupling among Q fields and sleptons squarks). We reveal that the bilinear R-parity violating terms produce the appropriate amount of the baryon asymmetry in the parameter region, explaining the axion dark matter abundance, the smallness of µ- and R-parity violating interactions and atmospheric masssquared difference of neutrinos. Furthermore, in this model, the constraint for the Hubble parameter during inflation is mildly relaxed because the Q breaking scale is enhanced during inflation. Q-ball problem is also avoided if we use the LH u D-flat direction in the gravitymediated SUSY breaking scenario. 21

24 Although, in this paper, we have focused on the atmospheric mass-squared difference of neutrinos, it is interesting to discuss in more detail the neutrino masses and flavor mixings which will be the subject of future work. The SUSY DFSZ axion model without R-parity may explain the structure of the neutrino masses and flavor mixings without severe tunings of parameters. Acknowledgments We thank J. Kawamura for useful comments and improving the draft of this paper. We also thank T. Moroi for useful comments. K. A. is grateful to M. Yamaguchi for valuable discussion and also thank K. Fujikura and M. Ibe for useful comments. H. O. was supported in part by Grant-in-Aid for Young Scientists B) No. 17K14303) from Japan Society for the romotion of Science. A General coupling between AD/Q fields and the minimum during inflation In Sec.3, we consider the Q charges as Table 1 and LH u D-flat direction. However, we can consider other Q charges and other D-flat directions in the AD mechanism. Therefore, in this section, we consider the following general superpotential of the AD field and Q fields W AD S 1, Φ), W AD S 1, Φ) = y S m 1 φ n nm n+m 3, 74) where m, n are integers and y is a dimensionless parameter, giving rise to the following superpotential, W = W inf I) + W AD S 1, Φ) + W Q + W, W Q = κs 0 S 1 S 2 f 2 ), ) F W = α W AD S 1, φ) + β W I Q + h.c., 75) M where α, β are coupling constants. Let us investigate the minimum during inflation in this setup. After inflation, we must numerically study the dynamics of the AD/Q fields, and the results depend on the detail of parameters. In this respect, we will postpone the investigation of their dynamics for a future work, but the results will be similar to Secs and Here we consider n, m 2. In our model presented in Sec. 3, namely n = 2, m = 3, the following calculation is simplified. We can apply the following calculation to the R-parity conserving case. First, we consider the scalar potential for the AD/Q fields. We assume that the Kähler potential is the same as in Eq. 20). Then the scalar potential for the AD/Q field is described 22

25 as V = V Hubble + V soft + V A + V F, V A = a HH + a m m 3/2 ) y S m 1 φ n nm n+m 3 V F = κ 2 S 1 S 2 f κs 0S 2 my S m 1 + b HH + b mm 3/2 )κs 0 S 1 S 2 f 2 ) + h.c., 1 φ n nm n+m 3 2+ κ 2 S 0 S y 2 S 1 2m φ 2n 1), 76) M 2n+2m 6 where V Hubble and V soft are the same as in Eqs. 22) and 23). Let us investigate a minimum of the potential. Ignoring the soft supersymmetry breaking effect, the phase dependent term in this potential is m 1 V phase = 2ˆκ 2 f 2 Ŝ 1 Ŝ 2 cosθ S1 + θ S2 ) 2ˆκŷ mŝ1 ˆφn Ŝ 0 Ŝ 2 nm n+m 3 cosθ S0 + θ S2 m 1)θ S1 nθ φ + ζ ) 2ˆκ b ˆ H HŜ0)Ŝ1Ŝ2 cosθ S0 + θ S1 + θ S2 + η ) + 2ˆκ b ˆ H f 2 HŜ0) cosθ S0 + η ) m ˆφn 2ŷ â H H Ŝ1 nm n+m 3 cosmθ S1 + nθ φ + ξ ), 77) where ζ, η, ξ are some numerical constants. Then minimum of the phases are θ S1 + θ S2 0, θ S0 + θ S2 m 1)θ S1 nθ φ ζ, mθ S1 + nθ φ ξ. 78) We assume H S 0 < f 2 which is confirmed in the same way of Sec We set all dimensionless parameters as O1). Ignoring the soft mass terms, the first derivatives in each radial direction of field at the minimum of phase directions are V Ŝ0 V Ŝ1 V Ŝ2 2ˆκŜ2 ˆκŜ0Ŝ2 mŷ Ŝ 1 2ˆκ 2 Ŝ 2 Ŝ 1 Ŝ 2 f 2 ) m 1 ˆφn nm n+m 3 2mm 1)ŷ Ŝ 1 ) +2ˆκ 2 Ŝ 0 Ŝ c0 H 2 Ŝ 0, nm n+m 3 + 2ˆκ 2 Ŝ 02Ŝ 1 + 2mŷ ) 2 2m 1 Ŝ ˆφ2n 1) 1 M 2n+2m 6 2ˆκ 2 Ŝ 1 Ŝ 1 Ŝ 2 f 2 )+2ˆκŜ0 V ˆφ Ŝ 1 2mŷ m 1 ˆφn 1 M n+m 3 ˆκŜ0Ŝ2 mŷ Ŝ 1 m 2 ˆφn ˆκŜ0Ŝ2 mŷ Ŝ 1 ˆκŜ0Ŝ2 mŷ Ŝ 1 m 1 2mŷ â H HŜ1 ˆφn nm n+m 3 m 1 ˆφn nm n+m 3 m 1 ˆφn nm n+m 3 m 1 ˆφn nm n+m 3 ) +2c 2 H 2 Ŝ 2, ) + 2n 1)ŷ ) 2 Ŝ 1 M 2n+2m 6 ) 2c 1 H 2 Ŝ 1, 2m ˆφ2n 3 2ŷ â H HŜ1 M n+m 3 m ˆφn 1 2c 3 H 2 ˆφ. 79) 23

26 From the extremal conditions V ˆφ i Ŝ0 = Ŝ2 f 2 Ŝ1, ˆφ = 0, one of the extrema is given by ˆκ Ŝ2 mŷ Ŝ1 m 1 ˆφ n ˆκ 2 Ŝ1 2 + ˆκ 2 Ŝ2 2 + c 0 H 2 nm n+m 3, k m 2 m n âh + Ŝ1 k ˆφ, k 2m m 2 n 2 â 2 H + 4c 1k 2m 4 m2 m 1) 2k 2m 4 m2 m 1) n 2 + k 2m 2 m) n 2 + k 2m 2 m) HM n+m 3 ) 1 n+m 2, 80) ŷ where k is some numerical constant which depends on â H, c 1, c 2, n, and m. To get this extrema, we must satisfy the following condition during inflation, Ŝ1 f. 81) Note that in Sec.3, Eq. 81) is satisfied under H > m φ and f GeV. From the potential and Eq.81), S 0 and S 2 obtain large positive masses of Ŝ1. Thus we assume Ŝ0, Ŝ2 and all phase directions are fixed at the minimum. Then the mass matrix of Ŝ1 and ˆφ is 2 V 2 V where 1 2 V 2 Ŝ1 Ŝ1 1 2 Ŝ1 Ŝ1 2 V ˆφ Ŝ1 m2 m 1)2m 3)ŷ ) 2 2m 4 Ŝ ˆφ2n 1 n 2 M 2n+2m 6 mm m 2 1)ŷ â H HŜ1 ˆφn nm n+m V 2 Ŝ1 ˆφ 2m2 m 1)ŷ ) 2 Ŝ 1 nm 2n+2m 6 m 1 mŷ â H HŜ1 ˆφn 1 M n+m V 2 ˆφ ˆφ m2 2n 1)ŷ ) 2 Ŝ 1 nm 2n+2m 6 2m 3 ˆφ2n 1 Ŝ1 ˆφ 2 V, 82) ˆφ ˆφ + m2m 1)ŷ ) 2 2m 1) Ŝ ˆφ2n 1) 1 M 2n+2m 6 c 1 H 2, 83) + 2n 1)mŷ ) 2 2m 1 Ŝ ˆφ2n 3 1 M 2n+2m 6, 84) 2m 1) ˆφ2n 1) + 2n 3)n 1)ŷ ) 2 2m Ŝ ˆφ2n 4 1 M 2n+2m 6 n 1)ŷ â m H HŜ1 ˆφn 2 M n+m 3 c 3 H 2. 85) One can realize the positive eigenvalues of Eq. 82) under a H c 1, c 3. 24

27 References [1] Y. Fukuda et al. Super-Kamiokande Collaboration), hys. Rev. Lett. 81, ). [2]. Minkowski, hys. Lett. 67B 1977) [3] T. Yanagida, Conf. roc. C ) [4] M. Gell-Mann,. Ramond, and R. Slansky, Conf. roc. C ) , arxiv: [hep-th]. [5] R. N. Mohapatra and G. Senjanovic, hys. Rev. Lett ) 912. [6] J. Schechter and J. W. F. Valle, hys. Rev. D ) [7] L. J. Hall and M. Suzuki, Nucl. hys. B 231, ). [8] R. Barbier et al., hys. Rept ) 1 [hep-ph/ ]. [9] M. Endo, K. Hamaguchi and F. Takahashi, hys. Rev. Lett ) [hepph/ ]. [10] S. Nakamura and M. Yamaguchi, hys. Lett. B ) 389 [hep-ph/ ]. [11] M. Dine, R. Kitano, A. Morisse and Y. Shirman, hys. Rev. D ) [hepph/ ]. [12] M. Dine, W. Fischler and M. Srednicki, hys. Lett. B 104, ); A. R. Zhitnitsky, Sov. J. Nucl. hys. 31, ) [Yad. Fiz. 31, )]. [13] R. D. eccei and H. R. Quinn, hys. Rev. Lett. 38, ). [14] S. Weinberg, hys. Rev. Lett. 40, ). [15] F. Wilczek, hys. Rev. Lett. 40, ). [16] J. reskill, M. B. Wise and F. Wilczek, hys. Lett. 120B, ). ; L. F. Abbott and. Sikivie, hys. Lett. 120B, ). ; M. Dine and W. Fischler, hys. Lett. 120B, ). [17] A. D. Linde, Contemp. Concepts hys. 5, ), hep-th/ [18] D. H. Lyth and A. Riotto, hys. Rept. 314, ), hep-ph/ [19] S. Weinberg, Cosmology Oxford University ress, 2008). [20] I. Affleck and M. Dine, Nucl. hys. B ) 361. [21] M. Dine, L. Randall and S. D. Thomas, Nucl. hys. B ) 291 [hep-ph/ ]. 25