Baryogenesis via leptogenesis in multi-field inflation

In multi-field reheating after modular $j$-inflation we investigate the conditions under which baryogenesis via non-thermal leptogenesis can be successfully realized. We introduce three heavy right-handed neutrinos to the non-supersymmetric Standard Model of particle physics, assuming hierarchical neutrino masses. Considering a typical mass for the first right-handed neutrino of the order of $10^{11}~GeV$, suggested from the seesaw mechanism and also from concrete $SO(10)$ grand unification models, we obtain the allowed parameter space for viable baryogenesis. An upper bound for the inflaton mass as well as a lower bound for its branching ratio to the pair of lightest right-handed neutrinos are found and reported.


Introduction
Inflation [1,2] is widely accepted as the standard paradigm of the early Universe. The first reason is due to the fact that several long-standing puzzles of the Hot Big-Bang model, such as the horizon, flatness, and monopole problems, find a natural explanation in the framework of inflationary Universe. In addition, and perhaps the most intriguing feature of inflation, is that it gives us a causal interpretation of the origin of the Cosmic Microwave Background (CMB) temperature anisotropies [3], while at the same time it provides us with a mechanism to explain the Large-Scale Structure (LSS) of the Universe, since quantum fluctuations during the inflationary era may give rise to the primordial density perturbations [4].
Although single-field slow-roll inflation provides us with the best fit to the data, considering multi-field a grigorios.panotopoulos@tecnico.ulisboa.pt b nelson.videla@pucv.cl inflation offers a wide range of new features, which go beyond the predictions of single-field scenarios, and which could be detected. As a few examples, we may mention isocurvature perturbations, correlated or anticorrelated with the curvature perturbation, and detectable level of non-Gaussianity. In particular, for the current observational constraints on isocurvature perturbations, see Ref. [5]. More fundamentally, theories where the interactions of the Standard Model particles are unified with gravity, such as supergravity [6] and Superstring theory [7,8], give rise to multiple fields instead of one. For a comprehensive review on multi-field inflation, see e.g. [9,10].
The Standard Model (SM) of particle physics has been extremely successful describing very accurately a vast amount of observational data at energies that span many orders of magnitude. Despite its success, however, it is widely accepted that the SM is the low energy limit of some underlying fundamental theory. Perhaps the most straightforward evidence for physics beyond the SM is the tiny neutrino masses in the sub-eV range seen in solar, atmospheric, reactor and accelerator neutrino experiments [11,12,13,14,15,16,17]. Right-handed neutrinos are very well motivated hypothetical particles postulated to exist due to their appealing properties, which can be summarized as follows: they can explain small neutrino masses via the seesaw mechanism [18], they fit very nicely in the spinorial 16-dimensional representation of SO(10) GUT group [19], and finally they can explain the baryon asymmetry in the Universe via leptogenesis, see the next paragraph.
One of the goals of successful inflation must be the explanation of the baryon asymmetry in the Universe, which comprises on of the biggest challenges of modern cosmology. Both primordial Big Bang nucleosynthesis [20] and data from the CMB temperature anisotropies [21,22,23] show that the baryon-to-photon ratio is a tiny number, η B = 6.19 × 10 −10 [24]. This number should be calculable within the framework of known particle physics. Although several mechanisms exist, perhaps the most elegant one is leptogenesis [25]. A lepton asymmetry via the out-of-equilibrium decays of right-handed neutrinos is generated first, and then this lepton asymmetry is partially converted into baryon asymmetry via non-perturbative "sphaleron" effects [26].
Of particular interest is the non-thermal leptogenesis scenario [27,28,29,30,31,32], since the lepton asymmetry within the framework of non-thermal leptogenesis is proportional to the reheating temperature after inflation. The two key parameters of the big bang, namely the baryon asymmetry and the reheating temperature, are therefore linked together. Recently in [33] the author studied multifield reheating in interacting theories where the inflaton trajectory is weakly curved, and he realized this scenario in a particular example of modular inflation. It is the goal of this article to investigate the conditions under which baryogenesis via non-thermal leptogenesis in the scenario of multifield reheating considered in [33] can be successfully realized. Our work is organized as follows: In the next section we briefly present the theoretical framework, while the numerical results are discussed in section three. Finally we conclude summarizing our main findings in the last section.
2 Theoretical framework

Modular j-inflation
We briefly present modular j-inflation. For more details the interested reader may consult [33,34,35]. In the multifield scenario there are a bunch of inflatontype fields φ I , and another bunch of decay products χ A . The total Lagrangian density takes the usual form where the non-trivial kinetic term L kin is given by with G IJ , G AB being the metric tensors in the field space for the φ fields and the χ fields, respectively, while the interaction term is given by with a potential V for the φ fields only, another potential U for the χ fields only, and an interaction term between the two types of scalar fields with a coupling constant g I,A for the triscalar interactions (with dimensions of mass), and another dimensionless coupling constant h I,A for the bi-quadratic interactions, leading to decay and scattering processes of the form Modular inflation, not to be confused with moduli inflation, where inflation is driven by moduli fields coming from Superstring compactifications on Calabi-Yau manifolds, is a two-field inflationary model with modular functions for the field target space. In modular j-inflation, a particular example of modular inflation, the field space has a non-trivial metric geometry characterized by the Poincaré metric with the two inflaton fields being φ 1 = µτ 1 and φ 2 = µτ 2 , where µ is some mass scale. The inflaton potential is given by where Λ is another mass scale, and F (x) is a dimensionless function. The potential can be Taylor expanded around its minimum φ 0 , and assuming that the difference φ − φ 0 is small, the inflaton potential may be approximated by a monomial of the form with φ = φ 1 being the dominant component along the inflationary trajectory, which is assumed to be only weakly curved. Therefore, during reheating the second inflaton field will be irrelevant, and the inflaton potential looks like the one of a usual single-inflationary model. Despite the similarity, however, the phenomenology of the scenario discussed here is different compared to genuine single-field inflationaty models, such as Natural Inflation [36,37,38] or chaotic inflation based on a mass term for the inflaton, see a couple of comments in the end of the next section. 3 In the multi-field reheating realized in the modular scenario of [33] the inflaton potential is computed to be leading to an inflaton mass In addition, the inflaton decays into bosonic decay products χ via a triscalar interaction term in the Lagrangian density with a coupling constant g

Non-thermal leptogenesis
In the scenario of non-thermal leptogenesis, the lepton asymmetry is given by where m φ is the mass of the inflaton, ǫ is the CPviolation asymmetry factor, T reh is the reheating temperature after inflation, and BR(φ → N 1 N 1 ) ≡ BR is the branching ratio of the inflaton decay channel into a pair of right-handed neutrinos φ → N 1 N 1 . The CP-violation asymmetry factor is defined by where Γ = Γ (N → lH) andΓ = Γ (N →lH † ) and it can be written down as where the maximum CP-asymmetry factor (assuming hierarchical neutrino masses, M 1 ≪ M 2 , M 3 ) has been computed to be [40] ǫ max = 3 8π with v = 246 GeV being the Higgs vacuum expectation value, M 1 being the mass of the right-handed neutrino of the first family N 1 , and ∆m 2 atm = 2.5 × 10 −3 eV 2 being the atmospheric neutrino mass difference from neutrino oscillation data [24].
During reheating [41,42,43] the particle production is taken into account via a phenomenological approach in which an additional term is added into the Klein-Gordon equation where Γ φ is the total inflaton decay. The reheating temperature after inflation is given by [33] T reh = 90 where M pl = 2.4 × 10 18 GeV is the reduced Planck mass, and g * counts for the relativistic degrees of freedom. In the SM g * = 106.75 [39]. The inflaton decays into right-handed neutrinos, and in the scenario of modular inflation realized in [33] into bosonic decay products χ as well. Therefore the total inflaton decay width Γ φ = Γ f + Γ b has two contributions, a fermionic one where y i is the Yukawa coupling of the inflaton to the right-handed neutrino N i , and a bosonic one [33] By definition and we may easily solve for the Yukawa coupling to express y 1 in terms of g, m φ for a given branching ratio as follows Finally, the initial lepton asymmetry Y L = n L /s is converted into baryon asymmetry Y B = n B /s via sphaleron effects [26] where n is the number density of leptons or baryons, s is the entropy density of radiation, s = (2π 2 h * T 3 )/45, and the conversion factor a is computed to be a = (24 + 4N H )/(66 + 13N H ) [44], with N H being the number of Higgs doublets in the model. In the SM without supersymmetry with only one Higgs doublet N H = 1 and a = 28/79. We thus obtain the final expression for the baryon asymmetry assuming that m φ > 2M 1 , sinδ = 1 (maximum CP asymmetry factor), and that the other two channels φ → N 2 and φ → N 3 are kinematically closed. The 4 baryon asymmetry is related to the baryon-to-photon ratio, and therefore Y B takes the observational value Using the expressions for the total inflaton decay rate as well as the reheating temperature after inflation we may express g as a function of the inflaton mass for a given branching ratio and accordingly for y 1 We see that g and y 1 scale differently with the inflaton mass, namely g ∼ m

Numerical results
In the following we shall assume a typical value for the mass of the first right-handed neutrino M 1 = 10 11 GeV [45]. Then the model is characterized by three free parameters, namely the Yukawa coupling y 1 , the triscalar coupling constant g and the inflaton mass m φ . It is more convenient, however, to work with the branching ratio BR(φ → N 1 N 1 ) since it is a dimensionless number taking values in the range [0,1].
Imposing the observational constraint Y B = 7.89 × 10 −11 we obtain the coupling constant g as a function of the inflaton mass for different values of the branching ratio shown in Fig. 1. Note that in the usual baryogenesis via non-thermal leptogenesis scenario, at least in non-supersymmetric models, the inflaton decays into right-handed neutrinos only, there are no scalar decay products, and the branching ratio BR ≃ 1 (see e.g. [32]). Therefore we have considered here low branching rations up to 0.1. We see that m φ ≫ g as anticipated in [33]. For each point of the curves shown in Fig. 1 we can compute the Yukawa coupling and the reheating temperature after inflation using the formulas presented in the previous section. We find that both y 1 and T reh increase with the inflaton mass, as shown in Fig. 2 and 3. For a given inflaton mass the reheating temperature after inflation decreases with the branching ratio. Therefore, the scenario studied here predicts a higher T reh compared to the standard discussion where BR ≃ 1.
Since the model is non-supersymmetric, there are no bounds on T reh due to the gravitino problem [46,47].
Non-thermal leptogenesis, however, works if T reh ≪ M 1 , and therefore we impose the condition T reh ≤ M 1 /100, which implies an upper bound for the inflaton mass proportional to the BR, and therefore for a viable baryobenesis via non-thermal leptogenesis the inflaton mass must take values in the range (29) and this in turn implies a lower bound on the branching ratio, BR ≥ 0.003. In figures 1-3 the upper bound of the inflaton mass is shown. We see that the specific scenario with a low branching ratio discussed in this work requires a relatively light inflaton, m φ < 10 13 GeV , while in the chaotic inflationary model based on a mass term for the inflaton, (1/2)m 2 φ 2 , it is well-known that the COBE normalization requires an inflaton mass m > 10 13 GeV . Therefore, this single-field model cannot work here.
Finally, we may now show the mass scale µ as a function of the mass scale Λ for a given value of the inflaton mass. But before that, since any viable inflationary model first should be compatible with the spectral index n RR and tensor-to-scalar ratio r bounds, we briefly summarize here the main results obtained in [48]. According to that work, within the framework of modular j-inflation the spectral index was found to be n RR = 0.96, while r was found to take values in the range 10 −8 ≤ r ≤ 0.08, compatible with the values reported by the Planck collaboration [22,23] as well as the BICEP2/Keck/Planck collaboration [49]. In addition, the scale Λ is allowed to take values in the range 10 −6 ≤ Λ/M p ≤ 10 −4 , and it is lower than the mass scale in Natural Inflation, where it is of the order of the GUT scale [36,37,38].
In Fig. 4 we show the mass scale µ as a function of the mass scale Λ for three different values of m φ . The scale Λ takes values in its allowed range mentioned before, while the scale µ is of the order of the GUT scale.

Conclusions
To summarize, in this article we have studied baryogenesis via non-thermal leptogenesis in multi-field reheating realized in a particular example of modular inflation. We have assumed hierarchical neutrino masses in the seesaw mechanism scenario introducing three heavy right-handed neutrinos N i without supersymmetry, and we have investigated under what conditions the  model is viable. The inflaton φ decays into heavy righthanded neutrinos, and into bosonic decay products too.
We have focused to the case of a small branching ratio φ → N 1 N 1 , and we have obtained the allowed parameter space corresponding to successful baryogenesis. An upper bound for the inflaton mass as well as a lower bound for its branching ratio into the pair of lightest right-handed neutrinos are obtained and reported. As a final remark, we have assumed that perturbative reheating applies. Non-perturbative preheating effects af-   ter inflation deserve a more detailed analysis. We hope to be able to address this point in a future work.