Soft leptogenesis in the NMSSM with a singlet right-handed neutrino superfield

In this work, we explore soft leptogenesis in the NMSSM framework extended by a right-handed neutrino superfield. We calculate the CP asymmetry, ε, and find it to be non-zero at tree level without using thermal effects for the final state particles. This is in contrast to soft leptogenesis in the MSSM extended by a right-handed neutrino superfield where thermal effects are essential. The difference arises due to the presence of a 3-body decay of the sneutrino in the NMSSM that violates lepton number at tree level. Apart from this, we also find that ε 6= 0 if the additional singlet scalar has a complex vacuum expectation value while all the other NMSSM parameters including the soft SUSY breaking ones relevant for CP asymmetry remain real. We estimate the order of magnitudes of these parameters to produce sufficient baryon asymmetry of the Universe. ∗ waleedabdallah@hri.res.in † abhass@prl.res.in ‡ aks@prl.res.in 1 ar X iv :1 91 1. 03 36 3v 1 [ he pph ] 8 N ov 2 01 9


I. INTRODUCTION
It is well known that the observable Universe has an asymmetry between baryons and anti-baryons [1], often called the problem of baryogenesis. Over the years, many mechanisms have been proposed to create this baryon asymmetry. More recently, leptogenesis [2,3] has become a highly favoured model for baryogenesis, specially because this mechanism is naturally linked to neutrino masses. Adding right-handed (RH) singlet heavy neutrinos to the standard model generates neutrino masses by the seesaw mechanism [4][5][6][7][8][9]. These RH neutrinos can also decay to produce a scalar and SM leptons and if the decay violates CP for example due to interference between tree-level and loop-level decays owing to complex couplings, a lepton asymmetry is generated. Since, in the standard model (SM), the B − L symmetry is exact while the B + L symmetry is broken by the electroweak (EW) sphaleron processes [10], these sphaleron processes can convert the generated lepton asymmetry to baryon asymmetry.
Soft leptogenesis [11][12][13] pertains to generating lepton asymmetry at the tree-level itself due to mixing between the particle and anti-particle states of the RH singlet sneutrino,Ñ , because of the presence of soft SUSY breaking terms 1 . In the most minimal soft leptogenesis setup using minimal supersymmetric standard model (MSSM) [25,26] extended by one RH neutrino superfield,N , a CP asymmetry in the RH sneutrino sector is created only when thermal masses for the final products are considered. The asymmetry is present because the thermal phase space factors are different for bosons and fermions. Another work featuring soft leptogenesis looks at CP violation not just due to mixing between particle and antiparticle initial states but in decays and in the interference of mixing and decay [27]. In Ref. [28], it is shown that considering most generic soft trilinear couplings and one loop self energy contributions for sneutrino decay it is possible to generate CP violation even without finite temperature effects within the same setup.
However, the MSSM suffers from the so-called µ-problem [29] -there is no explanation to why the SUSY scale preserving µ-term (a direct SUSY mass term for the Higgs fields) should be of the same order as the soft SUSY breaking terms. The most straightforward solution to the µ-problem comes by promoting the µ-parameter into a field whose vacuum expectation value (vev) is determined, like the other scalar field vevs, from the minimization of the 1 Soft leptogenesis in different types of SUSY frameowrk can be found in [14][15][16][17][18][19][20][21][22][23][24]. scalar potential along the new direction [30]. Naturally, it is expected to fall in the range of the other vevs, i.e., of order O(M SUSY ). The next-to-minimal supersymmetric standard model (NMSSM) (for review see [30,31]) is the most simple and elegant model to solve this problem, where a singlet superfieldŜ is introduced to the MSSM superfields which gets non-zero vev. The NMSSM can be extended by a set of RH neutrino superfields to generate masses for the SM light neutrinos by the type-I seesaw (see [32][33][34] for the MSSM extended by RH neutrino superfield). This has been explored earlier in Ref. [35]. This extension also keeps the R-parity conserved if the sneutrinos do not get vevs [35].
In this work, by using the NMSSM extended by the RH neutrino superfield, we present a soft leptogenesis scenario that creates a lepton asymmetry at the tree-level decay of the RH sneutrino without using thermal mass factors. The CP violation is achieved by the mixing between the particle and anti-particle states. We also show that it is possible to obtain non-zero CP asymmetry even when all the soft parameters are real. Since the soft terms are responsible for creating the CP asymmetry instead of needing flavour effects as in usual leptogenesis, using only one generation of the RH neutrino superfield is enough. Even so, the setup can be easily extended to get the experimentally observed SM neutrino mass hierarchies and their mixing angles pattern [36,37].
The paper is organised in the following manner. We setup the model and segregate the parts required for soft leptogenesis in the next section (Sec. II). In the one following that, i.e., Sec. III, we calculate the CP asymmetry produced by decays of the various particle present in the model that contribute to non-zero CP asymmetry parameter ε at the tree level. We talk about the decays ofÑ as well as the scalar S in the model. In Sec. IV, we discuss the most crucial and important constraints and give a simple expression for ε is given. In Sec. V, we give and discuss the results of our calculation. We find that for successfully generating the observed baryon asymmetry of the Universe, we need ε ≈ O(10 −6 ). We also discuss what this could mean for various parameters of the model including the soft ones.
We finally conclude in Sec. VI.

II. MODEL
In the NMSSM, an extra singlet superfieldŜ is added to the MSSM Higgs sector [30].
Assuming explicitly Z 3 symmetry, the superpotential for the NMSSM with a singlet RH neutrino superfieldN in terms of the new singlet superfieldŜ and the MSSM doublet superfieldsĤ u andĤ d will be as follows [30]: whereL i andQ i are the SU (2) doublet superfields of leptons and quarks;Ê i andD i (Û i ) denote singlet down (up)-type quark superfields, respectively, and Y 's, λ's and κ are dimensionless couplings with generation indices (i, j = 1, 2, 3). After the singlet S obtains a vacuum expectation value (vev) S , an effective µ-term is generated: µ eff = λ S , which solves the so-called µ-problem [29]. The soft SUSY-breaking Lagrangian is given by whereL i andÑ are the scalar components ofL i andN superfields, respectively. CP is spontaneously violated when the scalars H u , H d , S attain vevs with relative physical phases.
The vev of the singlet S is complex: Since leptogenesis occurs above the electroweak (EW) phase transition, we do not give vevs to the two Higgs doublets. In this case, spontaneous CP violation can occur only when sin δ = 0.

A. Terms relevant for soft leptogenesis
The terms from the superpotential required for leptogenesis via sneutrino decay are: Here we consider λ N , κ to be all real and positive. We also remove the i, j indices from the leptons and the u index from the Higgs superfield for brevity. The scalar potential is obtained using: with The fermionic part of the Lagrangian is given by: The soft SUSY-breaking Lagrangian terms that play a role in leptogenesis are: The superpotential and the soft breaking terms combine to give the following interactions forÑ and σ which could in principle contribute to soft leptogenesis due to mixing between the particle and anti-particle states through the soft terms: where σ = S − S .

III. CP ASYMMETRY
Because of the soft terms as well as the vev of S, there is a mixing between particle and anti-particle states of the sneutrino and the singlet scalar σ which is the dynamic part of S.
The squared mass matrices for the two of them are given by: where If A λ is real, the mass square eigenvalues of the sneutrino are: which have the following eigenstates: Similarly, one can write the mass square eigenvalues and eigenstates for the σ − σ * system.
Because of mixing between the particle and anti-particle states of sneutrino and singlet scalar, these systems are similar to K 0 −K 0 and B 0 −B 0 systems [38]. The evolution of these systems in the non-relativistic limit are driven by the Hamiltonian H defined as follows: where M is the mass matrix and Γ is the decay rate matrix of the corresponding system.
Finally, the decay rates of the time evolved particle and anti-particle states of theÑ −Ñ * and the σ − σ * system are calculated to get the final total CP asymmetry. Since bothÑ and σ are massive and can decay into each other as final products as can be seen in Eq. (11), there are two possibilities for decay. The first is when M 1 m σ . In this case,Ñ decays to produce the CP asymmetry. Decays of σ into a pair of right handed neutrino and sneutrinos is suppressed while if m σ M 1 , the 3-body decay ofÑ into σ, H * andL * is suppressed.
We consider these cases one-by-one.

A.Ñ decays
Irrespective of the mass of σ relative to the mass ofÑ , the sneutrino can decay into leptonic (sleptonic) and Higgs (Higgsino) final particles. The CP asymmetry generated from such a scenario was calculated in [11,12]. However, if M 1 m σ , the 3-body decay channel opens up and it leads to interesting consequences for soft leptogenesis as we show below.
For the snuetrino system, upto leading order in the off-diagonal terms 2 , the mass matrix can be calculated from the squared mass matrix in Eq. (12): The decay rate matrix can be written from the Eq. (11). It contains both diagonal as well as off-diagonal terms becauseÑ can decay into particle as well as anti-particle final states.
There are both 2-body and 3-body decays ofÑ assuming that the mass ofÑ is much larger than the mass of σ. with The solutions for the time evolution ofÑ andÑ * come from the Schrodinger like equation where ψ = {Ñ ,Ñ * } T . The solutions are obtained as, whereÑ 0 ,Ñ * 0 are the field values at t = 0 and a = (MÑ ) 11 If ∆ΓÑ ∆M as well, the argument of the trigonometric functions becomes ∆M t/2 such that we can write:Ñ where The CP asymmetry factor ε is defined as the ratio of the difference between the decay rates ofÑ andÑ * into final state particles with lepton number +1 and −1 to the sum of all the decay rates, i.e., where f,f are the final states with lepton number +1 and −1, respectively. This then gives us the following CP asymmetry parameter: The decay rates at the tree level itself for final states with lepton numbers ±1 are different because the factor p q = 1 as it is not a hermitian quantity. This requires non-zero offdiagonal terms to be present in the mass matrix as well as the decay rate matrix of the system with atleast one of them being complex.
We also note that while the usual soft leptogenesis done in the MSSM [11,12], to the leading order in soft terms, necessarily requires thermal phase space factors for the final state bosons (c B ) and fermions (c F ) to have ε ∝ ∆ BF = c B −c F c B +c F , we get an asymmetry even at zero temperature. This happens because of the presence of σ in the model which facilitates a 3-body decay which is not cancelled by the other terms. If we did not have this, at leading order in soft terms, 1 − The sum and difference of the ratios |q/p| 2 and |p/q| 2 can be written as: where Keeping terms upto the leading order in Γ 1 , we find that the sum of the ratios is 2 while the difference is twice the values of y/(x 2 − y 2 ) 1/2 with, such that the final CP asymmetry can be written as: From Eqs. (40) and (41) it is evident that there is non-zero CP asymmetry even with real A N , provided δ is sufficiently large.
Therefore it is possible to successfully generate non zero lepton asymmetry fromÑ −Ñ * system at the tree level without using thermal phase space factors for bosonic and fermionic final states. We discuss this more and give some numerical estimates in Sec. V for relevant parameters. For the moment, let's consider the decays of σ.

B. σ decays
Unlike theÑ decays where most of the decay products were massless, the final products of σ decay are massive. This creates two possible decay modes of σ according to the condition satisfied.
1. σ decays to N N andÑÑ . This happens when m 2 σ > 4M 2 1 , 2. σ decays only to N N . This happens when 16λ 2 The mass matrix and the Γ matrix of the σ − σ * system are respectively: where A non-relativistic Hamiltonian can be defined following Eq. (17). Immediately it can be seen that because of the absence of an off-diagonal term in the decay rate matrix of σ, the ratio corresponding to (p/q) 2 ofÑ decay, If we solve for the evolution of the σ − σ * system, the CP asymmetry parameter computed exactly analogously to theÑ −Ñ * system will be zero because of exact cancellation between the ratios |r/s| 2 and |s/r| 2 . Thus

IV. GENERAL CONSTRAINTS
The CP asymmetry depends in this model depends on a lot of parameters. However, we can constrain some of them by various considerations. In deriving the following general constraints, we take M 1 M ⇒ M 1 2λ N v S , M 1 m σ and A λ κv S for simplicity.
• The condition of out-of-equilibrium decays at T = M 1 is given by comparing the decay rate of theÑ with the Hubble parameter at T = M 1 : where Γ is the diagonal component of Eq. (20). Substituting it in Eq. (48) we get For M 1 m σ , α ≈ O(10 −2 ) and we may write Eq. (49) as • The way we derived the CP asymmetry requires well separated states [12], i.e., Γ ∆M as well as ∆Γ ∆M as stated before. This gives us two self-consistent limits: cos(3δ)ImA N + (sin(3δ) • Neutrino mass upper limits (m ν 0.1 eV [1,39]) put constraints on the Yukawa coupling strength Y N and the mass of the right handed (s)neutrino.
• Electric dipole moment calculations can constrain the CP violating phases that appear in the vevs of the two Higgs doublets and the scalar singlet S. In [40], they show that in principle δ u (the phase in the vev of H u should we go below the EW scale) and δ could be large as long as the relative phase is kept small. Since we keep vev only for S, we keep δ small. For more details about the EDM constraints on the NMSSM, see [41].

A. A simpler form for ε
We can write a simpler form for the CP asymmetry by using the approximations made and the general relationships between various parameters given above. The set of parameters governing ε are: To have enough asymmetry (ε O(10 −6 )), it is evident that κ 1 is required unless A N becomes very large. However large A N would violate the out-of-equlibrium decay condition in Eq. (50).

V. RESULTS & DISCUSSIONS
To obtain the baryon asymmetry of the Universe, η B we solve the simultaneous Boltzmann equations for theÑ number density, NÑ , and the B − L number density, N B−L which are as follows [42][43][44]: where KÑ = Γ H(z=1) is the Hubble parameter at z = 1 with z = M 1 T , N eq is the equilibrium number density ofÑ . They take the following forms: with m P l = 1.22×10 19 GeV being the Planck mass and g s is the number of relativistic degrees of freedom in NMSSM which we take ≈ 225 except for theÑ which is non-relativistic. In writing the Boltzmann equation for B − L number density, we neglect the ∆L = 2 scattering processes for washout and assume it is dominated mostly by inverse decays. The contribution to washout from the scattering processes is small because we are in the weak washout regime with KÑ 1. The final B − L number density thus created, N f B−L then converts to the baryon asymmetry by the sphaleron processes such that the ratio of the baryon number density to the photon number density, η B is: where g * g s 225, g 0 * is the effective number of relativistic degrees of freedom at recombination and a sph is the sphaleron conversion factor. Since we will solve the Boltzmann equations numerically, we use the complete form of ε given in Eq. (41) In Fig. 1, we show the typical value of the CP asymmetry that satisfies the observed We also assume A λ κv S for simplicity in finding the correct set of values for other parameters. Using these, we show the relation between A N and κ for different values of Y N in Figs. 2 and 3 as contour plots in ε.
In Fig. 2 we take a complex A N with ImA N = ReA N and vary κ and ReA N for Y N = 10 −4.8 (left panel) and Y N = 10 −4.7 (right panel). We find that κ 10 −5 and ReA N 10 7 GeV are needed to generate ε 10 −6 . The cyan colored region is ruled out by the out-of-equlibrium condition. The exculded region gets enhanced for larger value of Y N .
We keep A N real in Fig. 3 and vary κ and ReA N for similar values of Y N as before. It's clear that non zero asymmetry can be created even with real A N as long as δ = nπ. However compared to Fig. 2, we find that the bounds on A N and κ in Fig. 3 are stronger. Both Figs. 2 and 3 satisfy neutrino mass bounds and the conditions of Γ, ∆Γ ∆M .

VI. CONCLUSION
We have presented a new mechanism for soft leptogenesis in the context of the NMSSM with a singlet right handed neutrino superfield. Similar to soft leptogenesis in MSSM, we also generate CP asymmetry at the tree level owing to the CP violation occuring due to the difference between the mass and CP eigenstates similar to the K 0 −K 0 or the B 0 −B 0 systems. The difference lies in the fact that MSSM soft leptogenesis requires using thermal masses and phase space factors for boson and fermion final states without which there is no asymmetry. In the NMSSM where the singlet scalar S takes a vev, an asymmetry can be generated even with zero temperature fields. Further if there is spontaneous CP violation in the system with sin δ = 0, lepton asymmetry can be created without using any other complex parameter.
In the numerical analysis, we considered the mass scale of the RH sneutrino to be 10 7 GeV. This is well below the limits imposed by the cosmological gravitino problem [45]. We found that to generate sufficient asymmetry, one of the soft trilinear coupling A N needs to be 10 7 GeV and κ O(10 −5 ). This also tells us that A λ O(10 2 ) GeV. It would be interesting to explore the flavor effects in the present scenario that we leave for a future study.