Conformal GUT inflation, proton lifetime and non-thermal leptogenesis

In this paper, we generalize Coleman-Weinberg (CW) inflation in grand unified theories (GUTs) such as $\text{SU}(5)$ and $\text{SO}(10)$ by means of considering two complex singlet fields with conformal invariance. In this framework, inflation emerges from a spontaneously broken conformal symmetry. The GUT symmetry implies a potential with a CW form, as a consequence of radiative corrections. The conformal symmetry flattens the above VEV branch of the CW potential to a Starobinsky plateau. As a result, we obtain $n_{s}\sim 1-\frac{2}{N}$ and $r\sim \frac{12}{N^2}$ for $N\sim 50-60$ e-foldings. Furthermore, this framework allow us to estimate the proton lifetime as $\tau_{p}\lesssim 10^{40}$ years, whose decay is mediated by the superheavy gauge bosons. Moreover, we implement a type I seesaw mechanism by weakly coupling the complex singlet, which carries two units of lepton number, to the three generations of singlet right handed neutrinos (RHNs). The spontaneous symmetry breaking of global lepton number amounts to the generation of neutrino masses. We also consider non-thermal leptogenesis in which the inflaton dominantly decays into heavy RHNs that sources the observed baryon asymmetry. We constrain the couplings of the inflaton field to the RHNs, which gives the reheating temperature as $10^{6}\text{ GeV}\lesssim T_{R}<10^{9}$ GeV.


I. INTRODUCTION
Primordial inflation is a successful paradigm for the description of the early Universe and it is strongly supported by the current observational data [1][2][3][4]. Primordial perturbations, when the scales exiting the horizon (k ∼ aH), are eventually responsible for the structure formation in the Universe. From Planck 2015 [4,5], the key observables of inflation, namely, the scalar tilt and the ratio of tensor to scalar power spectra, are constrained as n s = 0.968 ± 0.006, r < 0.09 at 95% confidence level. The CMB power spectra is observed to be nearly adiabatic, scale invariant and Gaussian [4,6]. Although the physical nature of the inflaton is still uncertain [7,8], the models based on f (R) or a canonical scalar field with a flat potential are favoured with respect to the data. Since the inflationary scale is in general expected to be ∼ 10 16 GeV, it is natural consider the inflaton to be a scalar field associated with grand unified theory (GUT) groups, such as SU (5) and SO (10). The Shafi-Vilenkin (SV) model [9] is one of the first realistic model of inflation which was based on SU(5) GUT [10]. In this framework, inflation is a result of the spontaneous breaking of SU(5) → SU(3) c × SU(2) L × U(1) Y by a GUT field (24-plet adjoint Higgs) and an inflaton, which is a SU (5) singlet that rolls down to a vacuum expectation value (VEV). The success of the SV model is that it can lead to a successful baryogenesis after inflation and predicts a proton life time above the current lower bound [11,12]. In this model, the scalar field potential is of a Coleman-Weinberg (CW) form, according to which primordial gravitational waves are constrained by 0.02 ≤ r ≤ 0.1 [13]. Although the SV model is well within the current bounds of Planck 2015, several extensions of this model were studied to get smaller values of r. In [14][15][16][17][18], CW inflation was studied in the context of induced gravity, non-minimal coupling and brane-world scenario, where the tensor to scalar ratio was obtained to be r ∼ O 10 −2 − O 10 −3 . We note that all these modifications necessarily introduce an additional parameter whose value determines the shape of the inflaton potential in the Einstein frame 1 .
Moreover, extensions of the SV model within particle physics offer rich physics beyond the Standard Model (SM). Therefore, the SV model is embedded in a higher gauge group such as SO (10), which can be broken to the SM via an intermediate group G 422 = SU(4) c × SU (2) L × SU (2) R [21,22]. Obtaining successful inflation in SO (10) is more a Electronic address: sravan@ubi.pt, korumilli@sustc.edu.cn, sravan.korumilli@rug.nl b Electronic address: pmoniz@ubi.pt 1 For example, in the case of non-minimal coupling of inflaton to Ricci scalar (ξ), the value of ξ 1 in order to get shape of the potential to be like a Starobinsky plateau [19,20]. realistic with additional benefits to explain physics beyond SM, such as neutrino physics, matter anti-matter asymmetry through non-thermal leptogenesis, monopoles and dark matter (DM) [12]. For example, Ref. [23] considered a complex singlet scalar being coupled to right handed neutrinos (RHNs), followed by implementing type I seesaw mechanism. This approach unified inflation with Majorana DM together with the scheme of generating neutrino masses. In [24] an additional U(1) B−L symmetry was considered in the SM i.e., SU(3) c × SU(2) L × U(1) Y × U (1) B−L , where B − L symmetry can be spontaneously broken when the scalar field takes the VEV. In this setup, we can explain the baryon asymmetry of the Universe through non-thermal leptogenesis [22,[25][26][27]. Recently, CW inflation was studied in an extension with SO(10) and E 6 , pointing out the possibility of observing primordial monopoles [28].
Apart from models based on GUT theories, the Starobinsky model based on the R 2 gravity modification and the Higgs inflation [1,29,30] occupy a privileged position, with practically equal predictions in the (n s , r) plane where N is the number of e-foldings before the end of inflation. There has been a growing interest on embedding these models in string theory and supergravity (SUGRA) aiming for a UV completion [31,32]. Recently, a UV completion of the Starobinsky model was proposed in the context of non-local gravity inspired from string field theory [33,34]. Starobinsky like models were also developed in N = 1 SUGRA, namely, no scale [35] and α−attractor models [36] where an additional physical parameter leads to any value of r < 0.1. In [37] α−attractor models were studied in the non-slow-roll context where a new class of potentials were shown to give the same predictions. On the other hand, Higgs inflation is particularly interesting due to the fact that Higgs was the only scalar so far found at LHC. But for it to be an inflaton candidate compatible with CMB data, we require a very large non-minimal coupling (ξ 1) to Ricci scalar. It was known that a scalar field with large non-minimal coupling gives rise to a R 2 term considering 1-loop quantum corrections. Consequently, renormalization group (RG) analysis shows that Higgs inflation is less preferable compared to Starobinsky model [38,39]. This result not only applies to Higgs inflation but also to any arbitrary scalar with very large non-minimal coupling. Furthermore, in both R 2 and Higgs inflation the inflaton field rolls down to zero after inflation 2 . Differently, in GUT theories the inflaton field acquires a VEV due to its interaction with GUT fields.
The main goal of this paper is to generalize the SV model in order to achieve r ∼ O 10 −3 without introducing any additional parameters that affect the flatness of the inflaton potential (in Einstein frame), coasting towards a Starobinsky plateau 3 . In our construction, we introduce conformal symmetry (or local scale invariance) in a GUT model. It was shown by Wetterich [40] that scale symmetries play a crucial role in the construction of realistic cosmological models based on particle physics. Moreover, scale symmetries successfully explain the hierarchy of different scales such as the Planck and the Higgs mass [41][42][43][44]. Therefore, it is natural to consider scale invariance in constructing an inflationary scenario, through which we can obtain a dynamical generation of the Planck mass, inflationary scale and particle physics scales beyond SM. In this regard, we consider two complex singlet fields X , Φ of SU(5) or SO (10) and couple them to the Ricci scalar and adjoint Higgs field (Σ), such that the total action would be conformally invariant. We obtain inflation as a result of spontaneous breaking the conformal and GUT symmetries. The former occurs due to gauge fixing of one singlet field to a constant for all spacetime and the latter occurs due to Σ field taking its GUT VEV. Here the inflaton is identified with the real part of the second singlet (φ = √ 2Re [Φ]), whereas the imaginary part is the corresponding Nambu-Goldstone boson is assumed to pick up a mass due to the presence of small explicit soft lepton number violation terms in the scalar potential [23]. We also assume Φ carries two units of lepton number and it is coupled to the RHNs. Near the end of inflation, the inflaton is supposed to reach its VEV and also the global lepton number is violated. Thereafter, we study the dominant decay of inflaton into heavy RHNs producing non-thermal leptogenesis. We compute the corresponding reheating temperatures and also discuss the issue of producing the observed baryon asymmetry. Our study completes with an observationally viable inflationary scenario, predicting proton life time, neutrino masses and producing non-thermal leptogenesis from heavy RHNs.
The paper is briefly organized as follows. In Sec. II, we describe toy models with conformal and scale invariance. We identify the interesting aspects of spontaneous symmetry breaking of these symmetries leading to viable inflationary scenarios. In Sec. III, we briefly present the SV model and the computation of the proton life time. In Sec. IV we propose our generalization of the SV model by introducing an additional conformal symmetry. We report the inflationary predictions of the generalized model together with estimates of proton life time. In Sec. V we further explore the nature of inflaton couplings to the SM Higgs and singlet RHNs through type I seesaw mechanism. In the view of the dominant decay of the inflaton into heavy RHNs, we constrain the Yukawa couplings of the inflaton field compatible with the generation of light left handed neutrino masses. In Sec. VI we implement non-thermal leptogenesis and compute the reheating temperatures corresponding to the dominant decay of inflaton to heavy RHNs. We additionally comment on the necessary requirements for the production of observed baryon asymmetry through CP violation decays of RHNs. In Sec. VII we summarize our results pointing to future steps. We provide an Appendix. A summarizing the effects of geometric destabilization from fields space of inflaton and the presence of heavy fields in our model. In this paper we follow the units = 1, c = 1, m 2 P = 1 8πG .

II. CONFORMAL VS SCALE INVARIANCE
Models with global and local scale invariance (Weyl invariance (or) conformal invariance) are often very useful to address the issue of hierarchies in both particle physics and cosmology [41][42][43][45][46][47]. Models with these symmetries contains no input mass parameters. The spontaneous breaking of those symmetries induced by the VEV's of the scalar fields present in the theory, generates a hierarchy of mass scales e.g., Planck mass, GUT scale and neutrino masses 4 . Moreover, it is a generic feature that scale or conformal symmetry breaking induce a flat direction in the scalar field potential [40], which makes these models even more interesting in the context of inflation. Another motivation to consider scale invariance for inflationary model building comes from CMB power spectra which is found to be nearly scale invariant [4].
In this section, we present firstly a toy model (with two fields) that is (global) scale invariant and present the generic form of (scale invariant) potentials and their properties. We review the presence of a massless Goldstone boson that appears as a result of spontaneous breaking of global scale invariance. In the following, we discuss the two field conformally invariant model, in which case the presence of a massless Goldstone boson can be removed by appropriate gauge fixing. The resultant spontaneous breaking of conformal symmetry (SBCS) turns to be very useful to obtain a Starobinsky like inflation 56 . We will later explore the role of SBCS in a more realistic inflationary setting based on GUTs.

A. Scale invariance
Here we discuss a toy model with two scalar fields (in view of Refs. [40,51,60,61]) and point out interesting features that we later utilize in our construction.
A generic two field global scale invariant action can be written as where α, β are constants and ρ = φ χ , the generic function f φ χ here can be treated as quartic self coupling of the field φ [40,61]. The action (2) is scale invariant, i.e., invariant under global scale transformations g µν → e −2λ g µν , φ → e λ φ , χ → e λ χ for any constant λ (dilatation symmetry).
The interesting property here is that if φ ∝ χ there exists a flat direction for the field φ (see [40] for detailed analysis). This will be more useful in the context of local scale invariant model. Let us consider a scale invariant potential of the form where the couplings can in general depend on the ratio of the two fields i.e., φ/χ. If for example, we assume the couplings to be independent of the ratio of the two fields and consider the spontaneous breaking of scale symmetry i.e., the case with φ = 0, χ = 0, thus, as a result of minimizing the potential, we arrive at [61] with λ 2 m = λ φ λ χ and λ m < 0. In (5) we can re-define the coupling asλ then the potential (5) looks like a simple quartic potential We can alternatively have a potential of the form which also satisfies the constraint (3) and is different from (4). We will later see that the form of the potential in (8) gives a viable inflationary scenario. From (5) -(8) we can crucially learn how to define couplings as a function of the ratio of two fields in a scale invariant model. Of course, we only considered here simple toy models. However, we note that such field dependent couplings can be expected to arise in string theory and were applied in the context of early Universe [62]. The spontaneous breaking of scale symmetry occurs when one of the fields develops a VEV (let us take the field χ) which induces the emergence of a corresponding massless Goldstone boson (dilaton)χ = √ 6M ln that is associated with an arbitrary scale M ∝ m P [40]. By performing a Weyl rescaling of the metric g µν →g µν = χ √ 6M 2 g µν and φ →φ = M √ 6χ φ we indeed observe that the fieldχ is massless since the potential becomes independent of the fieldχ Although interesting cosmology and particle physics can be developed based on the scale invariant models, we need to constrain the implications of the massless dilaton present in the system [44]. It was shown that the dilaton can be gauged away if we consider a model with local scale symmetry [63].

B. Conformal invariance
A general action that is invariant under local scale transformations g µν → Ω −2 (x) g µν , φ → Ω(x)φ , χ → Ω(x)χ can be written as where the potential in the above action should also satisfy the condition (3). From the above action we can define an effective Planck mass m 2 ef f = χ 2 −φ 2 6 which evolves with time. In these theories, we would recover the standard Planck scale m P when the fields reach their VEV. Note that the field χ contains a wrong sign for the kinetic term but it is not a problem as we can gauge fix the field at χ = constant = √ 6M for all spacetime where M ∼ O (m P ). This particular gauge choice is so called c−gauge 7 which spontaneously breaks the conformal symmetry. It was argued that the theories in this gauge are of interest especially in cosmological models based on particle physics [44]. We will further see in this paper that fixing the scale M sources the hierarchy of mass scales related to inflation and particle physics (e.g., neutrino masses). In the inflationary models based on GUTs it is natural that the field φ takes a non-zero VEV, i.e., φ = 0 in which case it is useful to assume 6M 2 − φ 2 = 6m 2 P in order to generate Planck mass. Moreover, its also necessary to keep the evolution of the field φ √ 6M in order to avoid an anti-gravity regime.
in (10), SBCS via gauge fixing χ = √ 6m P leads to the Einstein frame action in and it is written as We can see that the above action leads to a Starobinsky like inflation as the potential acquires a plateau when ϕ m P (i.e., φ → √ 6m P ). In this case the inflaton rolls down to a zero VEV by the end of inflation because of the gauge fixing χ = √ 6m P and consequently Einstein gravity is recovered. In the next sections, we will study realistic GUT inflationary models where the inflaton rolls down to non-zero VEV and sources interesting implications in particle physics sector.

III. COLEMAN-WEINBERG GUT INFLATION
In this section we briefly review the Shafi-Vilenkin model [9,64]. It is one of the first realistic model of inflation which was based on SU(5) grand unified theory (GUT) . In this framework a new scalar field φ, a SU(5) singlet was considered and it weakly interacts with the GUT symmetry breaking field (adjoint) Σ and fundamental Higgs field H 5 . The tree level scalar potential is given by where the coefficients a, b, α and β are taken to be of the order of 8 g 2 , therefore the radiative corrections in (Σ, H 5 ) sector can be neglected. The coefficient γ takes a relatively smaller value and 0 < λ i g 2 and λ 1 max λ 2 2 , λ 2 3 . The GUT field Σ which is a 5 × 5 matrix can diagonalized as 7 Supergravity gauge was first realized in the context of 2T − physics based SUGRA models [44,63] where it was shown to be useful to obtain geodesic completeness of the theory. We follow this gauge choice in this paper as it allow us to explain hierarchy of scales in our model. 8 The field Σ interacts with vector boson X with a coupling constant g Various symmetry breaking patterns of SU(5) were studied in [65], among which the one with SU (5) where σ is scalar field that emerges from spontaneous breaking of SU (5). Substituting it in (12) the equations of motion for the σ field read where λ c = a + 7 15 b. Taking λ 2 λ c , the σ field quickly evolves to its local minimum of the potential given by Adding the radiative corrections due to the couplings − λ2 2 φ 2 TrΣ 2 and λ3 2 φ 2 H † 5 H 5 , the effective potential of φ gets to the CW form given by [9,64] where The (φ , σ) sector of effective potential is given by and µ = φ denotes the VEV of φ at the minimum, V 0 and C are dimensionfull and dimensionless constants respectively. Substituting (16) in (19) we obtain the effective potential for the field φ in the direction of σ ∝ φ.
which is the vacuum energy density i.e., V (φ = 0) and the constant C can be chosen such that V (φ = µ) = 0. Therefore, the potential (19) can be written as Following (16) the GUT field σ reaches its global minimum only when the inflaton field reaches its VEV by the end of inflation. The inflationary predictions of this model were reported in detail in [11,12]. This model was shown to be in good agreement with the spectral index n s = 0.96 − 0.967 and the tensor to scalar ratio 0.02 ≤ r ≤ 0.1, which is well consistent with the Planck 2015 data [4,13]. From the VEV of the singlet field φ we can compute the masses of superheavy gauge bosons as Taking A ∼ The key prediction of GUT models is proton decay p → π 0 + e + mediated by X, Y gauge bosons. The life time of proton can be computed using where m pr is proton mass and α G ∼ 1/40 is the GUT coupling constant. The current lower bound on proton life time is given by τ p > 1.6 × 10 34 years indicates M X ∼ 4 × 10 15 GeV [66,67].

IV. GUT INFLATION WITH CONFORMAL SYMMETRY
As discussed in Sec. II, conformal symmetry is useful to generate flat potentials and the hierarchy of mass scales. Therefore, embedding conformal symmetry in GUT inflation is more realistic and helpful to generate simultaneously a Planck scale m P along with the mass scale of X Bosons M X ∼ 10 15 GeV that sources proton decay. In this section, we extend the previously discussed CW inflation by means of introducing conformal symmetry in SU(5) GUT theory. We then obtain an interesting model of inflation by implementing spontaneous breaking of conformal symmetry together with GUT symmetry 9 . We start with two complex singlet fields 10 ) is identified as the inflaton. Gauge fixing the fieldX causes SBCS as discussed in Sec. II. It is worth to note here that the same framework we study here, based on SU(5) GUT, can be easily realized in the SO(10) GUT. Therefore, the two complex singlets of SU(5) considered here are also singlets of SO(10) [12,22].
The conformally invariant action with complex SU(5) singlet fields Φ,X can be written as where Here we assume the Φ field coupling to the Higgs field H 5 is negligible and not very relevant during inflation. We consider that the singlet field Φ is weakly coupled to the adjoint field Σ through the following tree level potential where the coefficients a ∼ b ∼ g 2 (gauge couplings g 2 ∼ 0.3). Following the discussion in section II we assume the coupling constants are field dependent, i.e., in (24) the coupling constants can be read asλ 2 = λ 2 f Φ X ,λ 1 = λ 1 f 2 Φ X which depend on the ratio of the fields Φ,X . We consider With the tree level potential in (24), the action (23) is conformally invariant under the following transformations The SBCS occurs with gauge fixingX =X * = √ 3M , where M ∼ O (m P ). We assume inflation to happen in a direction ImΦ = 0. Therefore, for the stability of inflaton trajectory we require the mass of ImΦ to be 11 m 2 ImΦ H 2 inf . To arrange this, we can add a new term to the potential (24) as such that the mass of the ImΦ in the inflationary direction is m 2 we can have m 2

ImΦ ImΦ=0
H 2 inf during inflation. In this way, we can successfully obtain the stability of the inflaton trajectory during inflation [71]. Note that when multiple non-minimally scalar fields are involved, it is in generally expected to induce geometrical destabilization effects due to the negative curvature of the fields space in the Einstein frame. This topic has been extensively studied in recent years [72][73][74][75][76][77]. In Appendix. A, we present details of fields space geometry and argue that these effects might be negligible in the model we study herein, deferring a detailed quantitative analysis for future investigations. Similarly to the SV model, we also consider here SU (5) Likewise to the SV model, we assume λ 1 λ 2 a, b and due to the coupling − λ2 2 φ 2 TrΣ 2 f φ √ 6M , the GUT field σ reaches to its local field dependent minimum given by 12 Note that the above local minimum of the GUT field remains the same even though there is a non-minimal coupling with the Ricci scalar. We can easily understand this by conformally transforming the action (23) into the Einstein frame.
After SU(5) symmetry breaking, the X gauge bosons become superheavy, whereas the field σ continues to follow the behavior of the field φ. The tree level potential for (φ, σ) sector is given by Substituting (29) in (23) and rescaling the field φ → 1 + λ2 λc φ, we obtain Since λ 1 λ 2 , the effective potential for the inflaton field φ due to the radiative corrections becomes where δV is the counter term, µ is the VEV of the field φ and V 0 is a constant. Using (29), choosing an appropriate , a normalization constant such that V ef f (φ = µ) = 0 and the vacuum energy density such that V (φ = 0) = V 0 = Aµ 4 4 , we obtain where A ∼λ 2 2 16π 2 . We note here that the CW potential we considered is the standard one obtained from 1-loop correction in Minkowski space-time. In the de Sitter background, 1-loop corrections are in principle different and their significance was discussed in literature [79][80][81]. Recently, in Ref. [82], it was argued that during slow-roll inflation we can neglect the contribution of 1-loop corrections in the gravity sector. In addition, the contributions from higher loops can also be neglected by the consideration of the slow-rolling scalar field. Refs. [83,84] provide quantum corrections calculated for the cases of non-minimally coupled scalar fields.
In order to get a Planck mass m P dynamically generated by the end of inflation, we should take the corresponding VEV of the inflaton field as We can see that M ≥ m P and µ → √ from (25) and by doing a conformal transformation of the action (31) with the effective potential (33) into Einstein frame, we obtain (expressing in the units of m P = 1) Under the conformal transformation, the mass scales in the Einstein frame must be redefined as µ 2 → µ 2 6M 2 − φ 2 −1 . This is very much an equivalent procedure to the 1-loop analysis of Higgs inflation. See Refs. [85][86][87][88] for a detailed discussion on the equivalence between Jordan and Einstein frames which exactly matches, if we redefine the mass scales accordingly by conformal factor. Subsequently, substituting (33) The kinetic term of (36) is similar the no-scale models [35]. Canonically normalizing the scalar field as φ = √ 6M tanh ϕ √ 6 yields the Einstein frame potential The corresponding VEV of the canonically normalized field is ϕ = √ 6 arctan µ √ 6M . The potential in (37) is a flattened version of CW potential (20). Concretely, due to SBCS, the shape of the potential above VEV φ > µ gets significantly flattened. In Fig. 1 we compare the CW potential of the SV model with the modified form (37) we obtained in our case. The shape of the potential reaches a plateau like in Starobinsky model when ϕ 1 i.e., φ → √ 6M . Inflation always starts near the plateau and continues to evolve as φ (25) is always positive and consequently that avoids an anti gravity regime. Note that the flat potential (37) is significantly different from the one of CW inflation, studied with positive non-minimal coupling in [15]. In the next subsection we show that the inflationary observables for the potential (37) exactly match those of Starobinsky and Higgs inflation.

A. Inflationary predictions and proton lifetime
We assume the standard Friedmann-Lemaître-Robertson-Walker (FLRW) background. Let us define the general definitions of slow-roll parameters as where H is the Hubble parameter and the prime denotes derivative with respect to e-folding number N = ln a (t) before the end of inflation. The scalar power spectrum is given by The scalar power spectrum amplitude at pivot scale k = 0.002 M pc −1 is measured to be P R * = 2.2 × 10 −9 [4]. The scalar spectral index up to the first orders in slow roll parameters is given by The running and running of running spectral index known as [89] The ratio of tensor to scalar power spectrum is The potential (37) when ϕ 1 (AV branch) can be approximated as The equation of motion of the canonically normalized field is which during the slow-roll regime reduces to where we took H inf ≈ V E (ϕ) 3 and N * is the 60 e− foldings before the end of inflation. Computing the slow-roll parameter using (45) we obtain Using (46) we can write the predictions for the scalar tilt (40) and tensor to scalar ratio (42) as which exactly match with the predictions of Starobinsky and Higgs inflation [1,30]. We emphasize that the predictions of our model in (47) are almost independent of the VEV of the inflaton field φ = µ.
In Table. I we present the inflationary predictions of the model together with the corresponding X bosons mass and proton life time using (21) and (22). We also show our results for the case when the inflaton field rolls from above VEV (AV) i.e., when φ > µ. The predictions of below VEV (BV) branch i.e., when φ < µ are not very interesting as those are nearly same in the original CW inflation without any conformal symmetry [12]. This is evident from Fig. 1 where we can see only the AV branch of the potential significantly different in our case, whereas the BV branch is nearly same as in the SV model. Therefore, our interest in this paper is restricted to AV branch. For this case, from Table. I we can see that the inflationary predictions of the model almost remains the same for any value of inflaton VEV. Note that even though the inflaton field values are trans-Planckian, the values of n s , r remain the same. This is due to the fact that when ϕ µ the shape of the potential is exponentially flat like in Starobinsky model. Therefore, inflationary predictions only depend on the potential plateau rather than the field values (shift symmetry).
In Fig. 2 we depict the evolution of field φ (also for the canonically normalized field ϕ) and slow-roll parameter for particular parameter values.   For both plots, we have taken µ = 1.12 mP.

V. TYPE I SEESAW MECHANISM AND NEUTRINO MASSES
In this section, we further extend our model through a type I seesaw mechanism with global lepton number symmetry, whose spontaneous breaking leads to the generation of neutrino masses. In this framework, we suppose the singlet field Φ carries two units of lepton number and is coupled to the three generation of singlet right handed Majorana neutrinos (RHNs), from [23] where l is the lepton doublet, τ 2 is the second Pauli matrix. Here Y D is the Yukawa coupling matrix of the SM Higgs coupling to the left handed neutrinos and Y N is the coupling matrix of the singlet field to the three generations of Majorana RHNs ν i R . In principle, we can also weakly couple the inflaton with the SM Higgs boson as We note that even with the new potential in (49), conformal symmetry in (23) can be preserved by the following additional transformations 13 Applying SBCS viaX =X * = √ 3M and computing 1-loop corrections due to the additional couplings to neutrinos (48) and SM Higgs, the effective potential of the field φ becomes where In (52) we assume the coupling constant Y i N to be at least O (10) smaller thanλ 2 and λ h Y i N , such that β f ∼ 20λ 2 2 and µ f ∼ µ. Therefore, during inflation the coupling of a singlet field to the adjoint scalar Σ dominates. Consequently, the inflationary predictions in Table. I are unaffected by these additional couplings to Higgs and singlet neutrinos. Moreover, since we impose λ h Y i N , the inflaton field dominantly decays to RHNs rather than to SM Higgs. Lets consider that the lepton number violation happens at a scale when φ = µ. Computing the mass matrix of singlet and doublet neutrinos in the basis of ν L , ν R , using the Einstein frame potential of (48), we have where v 2 = 246 GeV is the electroweak vacuum. The light left handed neutrino mass can be obtained from perturbative diagonalization of (53) as The mass of heavy RHNs is given by The essence of the seesaw mechanism is the generation of neutrino masses, resulting in light left handed neutrinos and heavy right handed neutrinos. Both are related here to the VEV of the inflaton field. The current Planck data indicates the sum of light neutrino masses constrained as m νi < 0.23 eV [90]. Therefore considering the light neutrino mass to be m ν L ∼ O(0.1) eV, (54) gives a relation Taking Y D ∼ O 10 −1 and from Table. I imposing µ ∼ 1.2 m P − 24.37 m P , we get 2.5 × 10 −6 Y i N 1.0 × 10 −5 . This supports our previous assumptions after (52), that the couplings to the RHNs have a negligible effect for inflation. Our generalization of the SV model successfully fits into explaining the origin of neutrino masses. We can also take Y D < O 10 −1 which results in smaller values for Y N < O 10 −6 . Taking Y N ∼ 10 −6 , the heavy RHN mass will be around m ν R ∼ 4 × 10 12 GeV. For Y N < O 10 −6 we can lower the masses of RHNs. In the next section we aim to study reheating in our inflationary scenario, taking into account the constraints we have derived so far.

VI. REHEATING AND NON-THERMAL LEPTOGENESIS
We consider reheating through a dominant decay of the inflaton into heavy RHNs 1415 which requires m ϕ 2m ν R . The mass of the canonically normalized field ϕ at the minimum of the potential is given by the second derivative of the potential (37) where we have taken a value for A ∼ 5 × 10 −12 from Table. I. We implement the scheme of non-thermal leptogenesis proposed in [22,96] which can give rise to baryogenesis through CP violating decays of RH Majorana neutrinos. In this section, we closely follow [25][26][27]. We consider: • Hierarchical masses for RHNs m ν 1 To arrange this we require the coupling constants to be Y N1 Y N2 ∼ Y N3 . We assume that the inflaton decays equally into the two heavy RHNs ν 2,3 R and the corresponding reheating temperature can be computed using [24,25] T R = 90 where g * = 105.6 is the number of relativistic degrees of freedom and the decay rate is given by 14 The inflaton could also decay into Higgs field but we have chosen the coupling of the Higgs field to the inflaton as λ h Y i N O 10 −6 . For these couplings, the decay rate of the inflaton to a pair of Higgs bosons is negligible [91,92]. However, there can be a period of parametric resonance in the phase of preheating right after the end of inflation, during which the number of Higgs particles can grow exponentially [91,92]. Around the VEV, the inflaton potential (37) can be approximated as Then we can apply the results of [91,92] to estimate the effect of parametric resonance. The inflaton field oscillates around the minimum asφ (t) ≈φ A (t) sin(mt) , is the amplitude of oscillations of the inflaton field. The regime of parametric resonance occurs as far aŝ h 8π ϕ and whenφ A drops to smaller values then standard perturbation theory dominates. To estimate the effect of parametric resonance in our case we compute the number of oscillations at the end of parametric resonance (N f ). Following estimates from [91] we 2π where t f is the instant when parametric resonance ends, by means of As a result, we can further obtain [91] N f ∼ λ h m P 6πmϕ 1 , since mϕ ∼ O 10 −6 from (57) and λ h 10 −6 in our case. Therefore, the effects of parametric resonance in our case is negligible for our chosen values of inflaton-Higgs couplings. 15 We ignore the effects of non-minimally coupled heavy fields Im[Φ], σ during preheating or reheating due to non-trivial fields space geometry in the Einstein frame [77,[93][94][95]. We defer these interesting studies for future investigation.
The masses of heavy RHNs are m ν 2,3 , which for Y 2,3 N ∼ 10 −8 − 10 −6 we have m ν 2,3 R ∼ 10 10 − 10 12 GeV. In Fig.  3 we plot the possible reheating temperatures 16 of our case taking c 1 ≈ 0 and c 2 = c 3 = 1. • The decays of RH Majorana neutrinos ν i R break the lepton number conservation and leads to CP violation. There are two decay channels where H and l denote the Higgs field and the lepton doublets of the SM. The (lepton asymmetry generated by the CP violation) decay of ν i R is measured by the following quantity CP asymmetry i can be computed for the dominant decays of ν 2,3 R using [26,[101][102][103] where f (y) = √ y −1 + (y + 1) ln 1 + 1 y , g (y) = √ y y − 1 .
Here, we only aim to constrain the range of values for i leaving for future studies the explicit computation of constraining Yukawa matrix Y ij D [25].
The lepton asymmetry is given by where n L is the difference between number of leptons and anti-leptons and s indicates the entropy density, Br i denotes the branching ratio • The production of RH Majorana neutrinos happens non-thermally and sufficiently late so that the produced lepton asymmetry sources the baryon asymmetry at a later stage. This essentially requires m ν 1 R T R so that the later decay of lightest RH Majorana neutrino ν 1 R does not wash away the produced lepton asymmetry by the heavy ones. We assume there is an accidental B − L conservation 17 such that sphaleron process is active which brings a part of the above lepton asymmetry into the baryon asymmetry (see Ref. [104][105][106] for details). As the reheating temperature in our case is T R ∼ 10 6 − 10 9 GeV (see Fig. 3), we take Y 1 N ∼ 10 −10 − 10 −9 such that m ν 1 R ∼ 10 8 − 10 9 GeV . Therefore, with values m ν 2,3 R ∼ 10 10 − 10 12 GeV , m ν 1 R ∼ 10 8 − 10 9 GeV and T R ∼ 10 6 − 10 9 GeV, we have met the conditions for successful leptogenesis which are m ν 2 Baryon asymmetry is proportional to the lepton asymmetry as n B s 28 79 The baryon asymmetry which is measured by the ratio of the difference between the number of baryons minus the anti-baryons n B to the entropy density in the present Universe is constrained [90] in the following form n B s = (6.05 ± 0.06) × 10 −10 .

VII. CONCLUSIONS
Coleman-Weinberg inflation [9] has been a successful and realistic model based on GUT and is consistent with the current Planck data with r 0.02 [13]. In this work, we have further generalized the framework of CW inflation with an additional conformal symmetry. Spontaneous breaking of conformal symmetry is useful to create a hierarchy of mass scales, therefore it is natural to realize this symmetry in GUT models. In this respect, two complex singlet fields of SU (5) or SO(10) were considered and are coupled to the GUT fields in a suitable manner. We have showed that this setup, upon spontaneous breaking of GUT and conformal symmetry, leads to an interesting inflationary scenario driven by the real part of the singlet field. In our model, the above VEV branch of CW potential gets flattened to a Starobinsky plateau, allowing for n s ∼ 1− 2 N and r ∼ 12 N 2 for N ∼ 50−60 number of e-foldings. Therefore, our model is observationally fits with the same predictions of the Starobinsky and Higgs inflation. Moreover, the VEV of the inflaton affects the masses of the superheavy gauge bosons that mediate the proton decay. We calculated the corresponding estimates for the proton life time above the current lower bound from Super-K data as τ p p → π 0 + e + > 1.6 × 10 34 . In the next step, we introduced a coupling between the complex singlet field with the generation of three singlet RHNs, where the singlet field is assumed to carry two units of lepton number. We implemented a type I seesaw mechanism, where spontaneous symmetry breaking of global lepton number results in generating neutrino masses. We put an upper bound to the inflaton couplings to RHNs, assuming inflation is dominated by inflaton couplings to GUT field. For the non-thermal leptogenesis to happen, we have considered a dominant decay of the inflaton into some of the RHNs and obtained the corresponding reheating temperatures as 10 6 GeV T R < 10 9 GeV. Furthermore, our proposed extension of CW inflation can be tested within future CMB and particle physics experiments [107].
In this work, we mainly restricted to a non-supersymmetric construction of GUT inflation with conformal symmetry. It would be interesting to consider this model in GUT based SUGRA framework with superconformal symmetries, which we defer for future investigations.
Let us consider f φ √ 6 = 1 − φ 2 6M 2 . As discussed in Sec. IV during the inflationary regime φ → √ 6M and since we have assumed λ 2 λ c , then we have σ 2 φ 2 . In this phase, action (A1) effectively reduces to We can re-parametrize the fields as φ = ρ sin θ τ = ρ cos θ , where θ = tan −1 φ τ . Note that inflaton trajectory considered in Sec. IV corresponds to θ = 0. Substituting (A6) in the action (A5) we get Conformally transforming the above action into the Einstein frame gives (in the units of m P = 1) The above action can be rewritten as the following, by introducing a field metric where . (A10) Here I, J = 1, 2 for ρ, θ respectively. Computing the Ricci tensor and Ricci scalar for the metric (A10) we obtain Notice that the Ricci scalar associated to the fields space is negative and unit. It is very similar to several Starobinsky like models of inflation and α− attractor models of SUGRA 19 , for which it was shown that geometrical destabilization could only occur towards the end of inflation [73]. The point to emphasize the following. With a suitable choice of potential for θ, fields space geometrical effects on inflationary epoch can be heavily suppressed (see e.g., [108] for more details). However, in recent studies, effects of heavy fields during preheating epoch have been explored in multifield non-canonical, non-minimal models of inflation [77,[93][94][95]109]. We opted to ignore such effects in our investigation and assume likewise that inflaton dominantly decays into the RHNs in Sec. VI. We leave for future an analysis of this interesting aspect.