Running non-minimal inflation with stabilized inflaton potential

In the context of the Higgs model involving gauge and Yukawa interactions with the spontaneous gauge symmetry breaking, we consider λϕ4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda \phi ^4$$\end{document} inflation with non-minimal gravitational coupling, where the Higgs field is identified as the inflaton. Since the inflaton quartic coupling is very small, once quantum corrections through the gauge and Yukawa interactions are taken into account, the inflaton effective potential most likely becomes unstable. In order to avoid this problem, we need to impose stability conditions on the effective inflaton potential, which lead to not only non-trivial relations amongst the particle mass spectrum of the model, but also correlations between the inflationary predictions and the mass spectrum. For concrete discussion, we investigate the minimal B-L\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B-L$$\end{document} extension of the standard model with identification of the B-L\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B-L$$\end{document} Higgs field as the inflaton. The stability conditions for the inflaton effective potential fix the mass ratio amongst the B-L\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B-L$$\end{document} gauge boson, the right-handed neutrinos and the inflaton. This mass ratio also correlates with the inflationary predictions. In other words, if the B-L\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B-L$$\end{document} gauge boson and the right-handed neutrinos are discovered in the future, their observed mass ratio provides constraints on the inflationary predictions.


Introduction
Current understanding about the origin of our universe is that, for a very brief moment at the beginning, our universe went through a period of rapid accelerated expansion known as inflation. Inflation scenario [1][2][3][4] was originally proposed to solve serious problems in the Standard Big-Bang Cosmology, namely, the horizon, flatness and monopole problems. In addition and more importantly from the view point of the current cosmological observations, inflation provides a mechanism to create primordial density fluctuations of the a e-mail: okadan@ua.edu b e-mail: draut@crimson.ua.edu early universe which seed the formation of large scale structure of the universe that we see today. In a simple inflation scenario, inflation is driven by a single scalar field (inflaton) that slowly rolls down to its potential minimum (slow-roll inflation). During the slow-roll era, the inflaton energy is dominated by a slowly varying potential, which causes the universe to undergo a phase of an accelerated expansion. Quantum fluctuations of the inflaton field are stretched to macroscopic scales by inflation to yield the primordial density fluctuations. After inflation, the inflaton decays to the Standard Model (SM) particles and the decay products heat up the universe (reheating). The success of big bang nucleosynthesis scenario requires the reheating temperature to be T R 1 MeV.
Recently the Planck 2015 results [5] have set an upper bound on the tensor-to-scalar ratio as r 0.11 while the best fit value for the spectral index (n s ) is 0.9655 ± 0.0062 at 68% CL. Hence, the simple chaotic inflationary scenario with the inflaton potentials V ∝ φ 4 and V ∝ φ 2 are disfavored because their predictions for r are too large. Among many inflation models, λφ 4 inflation with non-minimal gravitational coupling (ξφ 2 R, where φ is inflaton, R is the scalar curvature, and ξ is a dimensionless coupling) is a very simple model, which can satisfy the constraints by the Planck 2015 with ξ 0.001 [6][7][8].
Given that we need interactions between SM particles and inflaton for a successful reheating of the universe, a more compelling inflation scenario would be where the inflaton field plays another important role in particle physics. As an example of such a scenario, we may consider the (general) Higgs model, where a scalar (Higgs) field plays the crucial role to spontaneously break the gauge symmetry of the model, and we identify the Higgs field as inflation. The SM Higgs inflation [9][10][11][12][13][14][15][16][17][18] is nothing but this scenario, where the SM Higgs boson plays the role of inflaton with non-minimal gravitational coupling. For the observed Higgs boson mass of around 125 GeV, the stability of SM Higgs potential is very sensitive to the initial value of top quark pole mass.
Hence, we need more precise measurements for top quark pole mass [19] to obtain a conclusion about the SM Higgs potential stability. 1 For a large top quark pole mass, for example, M t = 173.34 GeV [24], the SM Higgs potential turns out to be unstable, and some extension of the SM is necessary to realize the SM Higgs inflation [25,26]. 2 It is interesting to apply the same idea as the SM Higgs inflation to the general Higgs model and identify the Higgs field of the model (not the SM Higgs field) as the inflaton in the presence of non-minimal gravitational coupling. For a simple example, see [37].
As in the SM, the Higgs field in the general Higgs model has the gauge, Yukawa and quartic Higgs interactions. For complete analysis of inflation scenario in the Higgs model, we consider the effective inflaton/Higgs potential by taking quantum corrections into account. In fact, we see that quantum corrections most likely cause an instability of the effective inflaton potential. Note that unless the non-minimal coupling (ξ ) is very large, the quartic inflaton coupling is very small [6][7][8]. Hence, quantum corrections to the effective potential are dominated by the gauge and Yukawa interactions. We consider the renormalization group (RG) improved effective potential described as where φ denotes inflaton, and λ(φ) is the running quartic coupling satisfying the (one-loop) RG equation of the form, Here, g and Y are the gauge and Yukawa couplings, respectively, and C g and C Y are positive coefficients whose actual values are calculable once the particle contents of the model are defined. Since the quartic coupling is very small, we have neglected terms proportional to λ (λ 2 term and the anomalous dimension term). The solution to the RG equation is controlled by g and Y , which are much larger than λ and independent of λ. Therefore, we expect that unless the beta function is extremely small, the running inflaton quartic coupling λ is driven to be negative in the vicinity of the inflation initial value, in other words, the effective potential has true minimum (far) away from the vacuum set by the Higgs potential at the tree level. A simple way to avoid this problem is to require the beta function to vanish at the initial inflaton value (the stationary condition of λ with respect to φ), namely, C g g 4 −C Y Y 4 = 0. 3 This condition leads to a relation between g and Y , equivalently, a mass relation between the gauge boson and fermion in the Higgs model. Since the Higgs quartic coupling at low energy is evaluated by solving the RG equation, in which the gauge and Yukawa couplings dominate, the resultant Higgs mass also has a relation to the gauge and fermion masses. The stability of the effective potential also requires the positivity of the second derivative of the potential, which leads to another constraint on the gauge and Yukawa couplings. In the slow-roll inflation, the inflationary predictions are determined by the slow-roll parameters defined with the potential and its derivatives, and therefore, the inflationary predictions have a correlation with the mass spectrum of the Higgs model.
In order to explicitly show the mass relation and the correlation between the particle mass spectrum and inflationary predictions, we take the minimal B − L model as an example. This model is a very simple, well-motivated extension of the SM, where the global B − L (baryon number minus lepton number) in the SM is gauged. Three right-handed neutrinos and the B − L Higgs field (which is identified as the inflaton) are introduced for the cancellation of the gauge and gravitational anomaly and the B − L gauge symmetry breaking, respectively. Associated with the B − L gauge symmetry breaking, the B − L gauge boson and the right-handed neutrinos acquire their masses. With the generation of the Majorana right-handed neutrino masses, the seesaw mechanism [40][41][42][43][44][45] for the light neutrino mass generation is automatically implemented in this model. Analyzing the RG evolutions of the B − L sector and the effective inflaton (B − L Higgs) potential, we show the particle mass spectrum and its correlation to the inflationary predictions. Through the correlation, the Planck 2015 results provide us with constraints on the particle mass spectrum. This paper is organized as follows. In the next section, we briefly review the λφ 4 inflation with non-minimal gravitational coupling at the tree level, and discuss the inflationary predictions in the light of the Planck 2015 results. In Sect. 3, we introduce the minimal B − L extension of the SM and calculate the RG improved effective Higgs potential. We show the particle mass spectrum derived from the stability conditions of the effective potential and its correlation to the inflationary predictions. We then compare our results for various values of the non-minimal gravitational coupling ξ with the Planck 2015 results. In Sect. 5, we discuss reheating scenario in the B − L Higgs inflation for the completion of our inflationary scenario. Section 6 is devoted to conclusions.

Non-minimal λφ 4 inflation at tree level
In the Jordan frame, the action of our inflation model is given by (hereafter we always work in the Planck unit, 19 GeV is the Planck mass) where f (φ) = (1 + ξφ 2 ) with ξ being a positive, dimensionless parameter, and the inflaton potential is Using the conformal transformation, g Eμν = f (φ)g μν , the action in the Einstein frame is descried as In the Einstein frame with a canonical gravity sector, we describe the theory with a new inflaton field (σ ) which has a canonical kinetic term. The relation between σ and the original inflaton field φ is given by where a prime denotes the derivative with respect to φ. The inflationary slow-roll parameters in terms of the original scalar field (φ) are expressed as where V E is the potential in the Einstein frame in terms of original field φ given by (2.6) The amplitude of the curvature perturbation R is given by which should satisfy 2 R = 2.195 × 10 −9 from the Planck 2015 results [5] with the pivot scale chosen at k 0 = 0.002 Mpc −1 . The number of e-folds is given by where φ I is the inflaton value at horizon exit corresponding to the scale k 0 , and φ e is the inflaton value at the end of inflation, which is defined by max[ (φ e ), |η(φ e )|] = 1. The value of N depends logarithmically on the energy scale during inflation as well as on the reheating temperature, and is typically around 50-60. 4 The slow-roll approximation is valid as long as the conditions 1, |η| 1 and ζ 1 hold. In this case, the inflationary predictions, the scalar spectral index n s , the tensor-to-scalar ratio r , and the running of the spectral index α = dn s d ln k , are given by (2.9) Here the slow-roll parameters are evaluated at φ = φ I . Figure 1 shows inflationary predictions for the nonminimal λφ 4 inflation at the tree level for N = 50 (dasheddotted) and N = 60 (solid). Top panels show r vs. ξ (left) and n s vs. ξ (right). Both r and n s show asymptotic behavior for both small and large ξ values. In the minimal λφ 4 inflation limit with ξ = 0, we obtain r 0.31 (0.26) and n s 0.942 (0.951) for N = 50 (60). The plots also show that, for a larger e-holding number, we obtain a larger n s while a smaller r . The bottom-left panel shows the tree-level quartic coupling λ as a function of ξ . Note that λ is very small unless ξ 1. The inflationary predictions for n s and r for various values of ξ are depicted in the bottom-right panel along with the results from the measurements by Planck 2015 and Planck+BICEP2/Keck Array [5]. We see that the inflationary predictions for ξ 0.001 are consistent with the observations.

Running B − L Higgs inflation and stability of inflaton potential
In order to investigate the Higgs inflation with the stabilized inflaton potential, in this paper we take the minimal B − L extension of the SM as an example, where the anomalyfree U(1) B−L gauge symmetry is introduced along with a scalar field ϕ and three right-handed neutrinos NR i . The particle contents of our model are listed in Table 1. This model requires three generations of right-handed neutrinos to cancel all the gauge and gravitational anomalies. The B − L gauge symmetry is broken by the vacuum expectation value (VEV) of ϕ in its Higgs potential of  Table 1 Particle contents of the minimal B − L model. In addition to the SM particle contents, the right-handed neutrino NR i (i = 1, 2, 3 denotes the generation index) and a complex scalar ϕ are introduced Associated with the gauge symmetry breaking, the righthanded neutrinos acquire their Majorana masses through the Yukawa interaction, where we have taken the degenerate mass spectrum for the right-handed neutrinos, for simplicity. After the B − L symmetry breaking with the Higgs VEV ϕ = v B L / √ 2, the particle masses are given by Let us now consider the B − L Higgs inflation scenario. The action in the Jordan frame is given by 2 in the unitary gauge with the physical Higgs field φ identified as the inflaton. In the Einstein frame, the RG improved effective inflaton potential at the one-loop level is given by [47][48][49] where ≡ φ/ 1 + ξφ 2 , and we have neglected v B L , which is much smaller than the Planck mass. 5 The RG equations of the couplings at the one-loop level are given by [37] where a factor s defined as is assigned to each term in the RG equations associated with only the physical Higgs boson loop corrections [9][10][11][12][13][14][15][16]. 6 We find that this s-factor has no effect in our numerical analysis for a parameter region satisfying g 2 , Y 2 λ, and one may fix s = 1 as a good approximation. In our RG analysis, we have neglected the RG evolution of ξ , whose effect is found to be negligible. It has been pointed out in the last paper in [9][10][11][12][13] and emphasized recently in [39,50], some undetermined parameters are involved in the presence of the non-minimal gravitational coupling, when we connect couplings at high energies to those at low energies in the RG equations. However, as we will see in the following, all couplings in our model are very small, and for such small couplings, we can expect the effect of the undetermined parameters is negligibly small.
Let us now investigate the stability of the effective inflaton potential. In our analysis throughout this paper, we set the initial values of λ to be the one obtained in the tree-level analysis at the initial inflaton value φ = φ I , equivalently, I . Then, we consider the RG improved effective inflaton potential by taking into account the RG evolution of the quartic coupling with the initial condition at I . Since RGE evolves logarithmically, to a good approximation we can analytically solve the RGE of λ around I such that β λ is taken to be constant. The solution is given by 5 In this RG improved effective potential, we have identified the renormalization scale with according to the prescription proposed in [47][48][49], where it has been shown that the effective potential is frameindependent with this choice. In fact, we are especially interested in the parameter region of ξ 1 in this paper, so that φ and the identification reduces to the usual one. 6 There are a few different prescriptions for computing quantum corrections in the presence of the non-minimal gravitational coupling [17,18]. For recent, detailed computations of quantum corrections, see [47][48][49] and their results of 1-loop beta functions with the s-factor.
β λ ( I ) , (3.8) where I is the initial inflation scale and 0 is the scale such that λ( 0 ) = 0. So we require β λ ( I ) 16π 2 × λ( I ) to ensure that the potential is stable. As we have seen in the previous section (see Fig. 1), the quartic coupling is very small unless ξ 1. Consider ξ = 1, such that λ( I ) 10 −10 , we thus require β λ ( I ) 10 −8 . Hence, imposing β λ ( I ) = 0 is a good approximation for stability of the potential.
Unless ξ 1, the quartic coupling is very small, hence the beta function of the quartic coupling is approximately given by 7 The RG evolution is controlled by g and Y , which are independent of the initial value of the inflaton quartic coupling. Figure 2 shows the RG evolution of the inflaton quartic coupling in the vicinity of the initial inflaton value for various values of g and Y with a fixed ξ = 1.
In the left panel, the solid, the dashed and the dot-dashed lines denote the RG evolutions for g = 0.01, 0.011 and 0.009, respectively, with the fixed value of Y = 0.0237. In the right panel, the solid, the dashed and the dot-dashed lines denote the RG evolutions for Y = 0.0237, 0.0214 and 0.0261, respectively, with the fixed value of g = 0.01. Note that β λ ( I ) = 0 is satisfied with the parameter choice for the solid lines, g = 0.01 and Y = 0.0237. We can see from Fig. 2 that if the condition of β λ = 0 is violated even with ±10 % deviations for the values of g or Y , the running quartic coupling quickly becomes negative in the vicinity of I (see the dashed and the dot-dashed lines). This fact indicates that the B−L gauge symmetry breaking vacuum at φ = v B L is unstable and the effective potential develops a true vacuum with a negative cosmological constant. The quantum corrections through the gauge and Yukawa coupling completely change our inflationary scenario from the one at the tree level.
In order to avoid this instability, we impose not only the condition of β λ = 0 but also dβ λ /d > 0 at I . 8 From the 7 In this paper, we are interested in this case, otherwise the beta function is so small that the inflaton quartic coupling is almost RG invariant. Although the tree-level analysis is valid in this case, the gauge and Yukawa couplings are too small to yield any impacts in the experimental point of view. 8 To be precise, β λ is not necessary to be exactly zero to stabilize the inflaton potential. We may consider a more general condition 96g 4 − 3Y 4 λ 2 . Since λ 1, β λ = 0 is a good approximation. In addition, although g, Y λ can approximately satisfy this condition, we are not interested in such very small couplings for which RG evolutions are negligible. first condition, the Yukawa coupling Y is determined by the gauge coupling, which we take as a free parameter in our analysis, along with the others, ξ and v B L . The second condition ensures that the effective potential is monotonically increasing in the vicinity of I , and yields a lower bound on g > g min . 9 When we analyze the global structure of the effective potential, we can notice that there exists an upper bound on g < g max . For a large g > g max , the effective potential develops a local minimum at < I , so that the inflaton field will be trapped in this minimum after inflation. 10 A second inflation then takes place until the vacuum transition 9 From Eq. (3.9) we find g 2 min ∼ λ. The discussion here is applicable only for g 2 λ, and g min is not the general lower bound on the gauge coupling. For g 2 , Y 2 λ, the running effect on all the couplings is negligible and hence the inflaton potential is stable. 10 There is an interesting possibility that such a local minimum can be lifted up once thermal effect on the scalar potential is taken into account [50]. In our case, as we will see in Sect. 5, reheating temperature is not high enough to remove this local minimum. from this local minimum to the true B − L symmetry breaking vacuum. To avoid this problem, the parameter region is restricted to be in the range of g min < g < g max . Figure 3 shows the effective potential for various values of the gauge coupling (g) for fixed ξ = 0.1. In the left panel, the solid line depicts the effective potential for g = 0.041, while the dot-dashed line is for g = 0.046. In this example, we find the upper bound as g max = 0.0425. We can see that the effective potential develops a local minimum for g = 0.046 > g max . For various values of g < g max , the effective potentials are shown in the right panel.

Inflationary predictions and low energy observables
Under the stability conditions, g min < g < g max with various value of ξ , we now calculate the inflationary predictions with the effective inflaton potential. Since we refer the results in the tree-level analysis for λ( I ) for a fixed ξ and impose the stability condition β λ ( I ) = 0, our prediction for the tensor-to-scalar ratio r is the same as the one obtained in the tree-level analysis. However, the RG evolution of the inflaton quartic coupling alters the other inflationary predictions, n s and α, from those obtained in the tree-level analysis, because they are calculated by the second and third derivatives of the effective potential (see Eqs. (2.5) and (2.9)).
Let us first derive an analytic formula for the deviations of n s from the tree-level prediction. In the effective inflaton potential, the inflation quartic coupling is not a constant, but a function of = φ/ 1 + ξφ 2 . With the stability condition β λ ( I ) = 0, we calculate the derivatives of the effective potential as where a prime denotes the derivative with respect to φ, V E tree and V E tree are evaluated to be the same as those in the treelevel analysis, and we have used λ (φ I ) ∝ β λ ( I ) = 0 under the stability condition. Rewriting d/dφ in terms of d/d , we have where tree and η tree are the slow-roll parameters evaluated in the tree-level analysis, and Here, a dot is the derivative with respect to . Sinceβ λ is non-zero, η is deviated from its tree-level value η tree . Using Eqs. (3.6) and (3.9), we obtain where in the last expression we have used the stability condition Y 4 32g 4 at I . Finally, we arrive at an expression for the spectral index as (see Eq. (2.9)) 11 11 In the same way, we can express the running of the spectral index, α, in terms of ξ , I , λ and g. However, as we will see in the following numeral analysis, the predicted α values are found to be very small and always consistent with the Planck 2015 results. Thus, we omit the expression for α. n s = n tree s + 2 η n tree s + For a fixed ξ , I and λ( I ) are determined by the tree-level analysis, so that the inflationary prediction is controlled by the gauge coupling g, which is the free parameter in our analysis. In our model, there are only three free parameters, ξ , v B L and g, with a fixed e-folding number N = 50/60. Once we fix ξ and N , I and λ( I ) are fixed by the tree-level analysis for inflation, and the inflationary predictions except for r are controlled by the gauge coupling g with its relation to the Yukawa coupling Y led by the stability condition β λ ( I ) = 0. In the B − L model, the particle mass spectrum is determined by the gauge, Yukawa and inflation quartic couplings at the scale v B L (see Eq. (3.3)), which are obtained by solving the RG equations in Eq. (3.6) from μ = I to μ = v B L . Since λ( I ) is very small, its RG evolution is determined by g and Y as we can see from its RG equation. Considering all of these facts, we expect that there exists a non-trivial mass relation in the particle spectrum and a non-trivial correlation between the inflationary predictions and the particle mass spectrum. In the following, the results of our numerical analysis will show such non-trivial relations. For simplicity, we fix v B L so as to yield m Z = 2g(v B L )v B L = 3 TeV, to be consistent with the current results from the search for Z boson resonance at the Large Hadron Collider [51][52][53]. Figure 4 shows the resultant inflationary predictions for a variety of ξ values with the input values of g in the range of g min < g < g max at I (N = 60), along with the contours given by the Planck 2015 results. As we have discussed above, the prediction for the tensor-to-scalar ratio is the same as the one in the tree-level analysis, while the predicted spec-   Fig. 4, we can see that, for ξ 1, the results show sizable deviations for g ∼ g max from those at the tree-level analysis depicted as the diagonal dashed line. Interestingly, the Planck 2015 results provide upper bounds on g, which are more severe than g max for ξ 0.001. In our numerical analysis, we can see that I approaches 1/ √ ξ from a smaller value as we increase ξ 1, and hence the deviation of the predicted n s value from the one in the tree-level analysis becomes smaller as we can see from Eq. (4.6) with the limit 1 − ξ 2 I → 0. In order to see the inflationary predictions as a function of g, we show our results in Fig. 5 for ξ = 0.0029, 0.01 and 0.1, with N = 60. The top-left and top-right panels show the inflationary predictions of n s and α = dn s /d ln k as a function of g in the range of g min < g < g max . The prediction for r is the same in the tree-level analysis, and r 0.108, 0.045 and 0.008, respectively, for ξ = 0.0029, 0.01 and 0.1. For a larger value of ξ , g max becomes larger. As g is lowered, the predicted n s value approaches the treelevel prediction. By numerically solving the RG equations for the couplings in Eq. (3.6) for a fixed g, we obtain the particle mass spectrum with m Z = 3 TeV. The mass ratio m N /m Z is shown in the bottom-left panel, while the bottomright panel shows the inflaton mass. 12 In both panels, the solid lines from left to right correspond to the results for ξ = 0.0029, 0.01 and 0.1, respectively. The resultant mass ratio is almost independent of ξ , but shows a splitting for g 5 × 10 −4 . For such a very small g, its corresponding Y determined by the stability condition is also very small, and hence both g and Y are almost RG invariant and the mass ratio is determined by Y ( I )/g( I ). However, in this case, the condition g 2 , Y 2 λ is no longer valid, and Y determined by β λ ( I ) = 0 depends on the input λ values. This is the reason why the bottom-left panel shows the splitting among three solid lines for g 5 × 10 −4 . For g 0.005, the RG evolution of λ is mainly determined by g and Y in its beta function, since λ( I ) is extremely small. For a very small g value 0.005, the effect of g and Y on the RG evolution of λ becomes negligible, and λ(v B L ) λ( I ).  Fig. 6. The dashed lines denote the results for N = 50, while the solid lines for N = 60. The inflationary predictions show a sizable deference for the two different N values, as shown in the tree-level analysis in Sect. 2. On the other hand, the particle mass spectrum weakly depends on N .
We show the results for the mass spectrum for large ξ values and N = 60 in Fig. 7. The solid lines from left to right corresponds to the results for ξ = 1, 10, 50, 150, 500 and 1000. For g 0.05, the resultant mass ratio in the left panel shows ξ -dependence. This is because g 2 , Y 2 λ is no longer valid for such a small g value, and the Y value determined by β λ ( I ) depends on λ( I ). We also show the

Reheating after inflation
Any successful inflation scenario requires the transition to the Standard Big Bang Cosmology after inflation. This happens via the decay of the inflaton into the SM particles; during the era the inflaton is oscillating around it potential minimum, and the decay products then reheat the universe. 13 We estimate the reheating temperature after inflation T R by using = H = T 2 R π 2 90 g * , where is the inflaton decay width, and g * is the effective degrees of freedom for relativistic SM particles when the reheating occurs, so that T R 0.55 (100/g * ) 1/4 M P .
(5.1) 13 In general, reheating can occur through the parametric resonance more effectively than the perturbative decay of the inflaton. For a detailed discussion of "preheating" in the Higgs inflation, see [54,55]. In our analysis we do not consider the preheating and estimate the reheating temperature using Eq. (5.1). The true reheating temperature could be much higher than our result.
From the success of big bang nucleosynthesis, we take a model-independent lower bound on the reheating temperature as T R 1 MeV. As we discussed in the previous section, the inflaton is much lighter than the Z boson and right-handed neutrinos. Thus, its decay width to the SM particles through off-shell process mediated by the heavy particles is too small to satisfy the lower bound T R 1 MeV. An efficient reheating process is possible when the inflaton in general has a coupling with the SM Higgs doublet (H ) such as (see, for example, Ref. [56,57] for phenomenology of inflaton through the coupling with the SM Higgs boson) Although this coupling is crucial for the reheating process, we assume λ 1 not to change our analysis for the B − L model in the previous sections. Since the inflaton quartic coupling is very small and hence the inflaton is light, it is most likely that the inflaton can decay to the SM particles only through its mixing with the SM Higgs boson through the λ coupling. After the B − L and electroweak symmetry breakings, we diagonalize the scalar mass matrix of the form where h is the SM Higgs boson, and φ 1 and φ 2 are the mass eigenstates. The relations among the mass parameters and the mixing angle are the following: where v SM = 246 GeV is the Higgs doublet VEV, and m h and m φ are the masses for h and φ, respectively.
The left panel in Fig. 9 shows the mass eigenvalue m φ 2 as a function of sin θ for various m φ values with m h = 125 GeV. In the θ = 0 limit, m φ 2 = m φ , while m φ = m h for θ = π/2. Although in this plot we show the results by using Eq. (5.5), we only consider the case with θ 1 as mentioned above, otherwise our results obtained in the previous sections are changed in the presence of a sizable λ . For θ 1, the mass eigenstate φ 2 (φ 1 ) is almost identical to φ (h).
The inflaton can decay to the SM particles through the mixing with the SM Higgs boson. We calculate the inflaton decay width as  The sharp drop for each solid line corresponds to the fact that m φ 2 becomes zero for a certain sin θ value as shown in the left panel. We find that the universe is sufficiently heated up for a mixing angle in the suitable range shown in Fig. 8. Finally, we check the theoretical consistency of our analysis. When we introduce the coupling in Eq. (5.2), the beta function of the inflaton quartic coupling is modified to 16π 2 β λ → 16π 2 β λ + 2λ 2 . (5.6) In order not to change our results in the previous sections by the introduction of λ , λ 2 should be negligibly small in the beta function. We then impose a condition For m 2 φ m 2 h , we obtain from Eq. (5.5) From Figs. 5, 7 and 9, we can see that this condition is satisfied for a large potion of the parameters space.

Conclusions
The inflationary universe is the standard paradigm in modern cosmology, which not only solves the problems in the Stan-dard Big Bang Cosmology, but also provide the primordial density fluctuations necessary for generating the large scale structure of the present universe. As a simple and successful inflationary scenario, we have considered the λφ 4 inflation with non-minimal gravitational coupling. With a suitable strength of the non-minimal coupling, the inflationary predictions of this scenario becomes perfectly consistent with the Planck 2015 results. It is more interesting if the inflaton can also play some crucial role in particle physics. We have considered the general Higgs model with the gauge and Yukawa interactions with the spontaneous gauge symmetry breaking. In the presence of the non-minimal gravitational coupling, the Higgs field can also play the role of inflaton. The analysis with the Higgs potential at the tree-level leads to the inflationary predictions consistent with the cosmological observations. However, once we take quantum corrections, the effective inflaton potential most likely becomes unstable. This is because the inflaton quartic coupling is extremely small in a large portion of the parameters space and the effective potential is controlled by the gauge and Yukawa couplings independently of the quartic coupling. In the renormalization group improved effective potential, we see that the running quartic coupling becomes negative in the vicinity of the initial inflaton value, indicating the instability of the effective potential. In order to avoid this problem, we have imposed the stability condition of vanishing the beta function of the inflation quartic coupling. This condition leads to a non-trivial relation between the gauge and fermion masses. Since the renormalization group evolution of the inflaton quartic coupling is mainly controlled by the gauge and Yukawa coupling, the inflation mass at low energy is determined by the couplings. Therefore, the mass spectrum of the gauge boson, fermion and inflation shows a non-trivial relation.
Since the inflaton potential is modified from the treelevel one, the inflationary predictions are altered from those obtained by the tree-level analysis. Although the prediction of the tensor-to-scalar ratio remains the same under the condition of the vanishing beta function, the predictions for the scalar spectral index and the running of the spectral index can be significantly altered. The fact that the effective potential is controlled by the gauge and Yukawa couplings implies a correlation between the inflationary predictions and the particle mass spectrum. Therefore, the observables at the gauge symmetry breaking scale correlate with the inflationary predictions which determined by physics at an extremely high energy compared to the gauge symmetry breaking scale.
By taking the minimal B − L extension of the Standard Model as a simple example, we have shown such a nontrivial relation in the particle mass spectrum driven by the stability condition of the effective inflaton potential. We also have calculated the inflationary predictions from the effective potential and found their dependence of the B − L gauge coupling. Therefore, the new particle mass spectrum of the B − L model, once observed, has an implication to the inflationary predictions. On the other hand, more precise measurements of the inflationary predictions yield a constraint on the B − L particle mass spectrum.
For completeness, we have also investigated reheating after inflation. Since the inflation is lighter than the Z boson and the right-handed neutrinos, its reheating process through the heavy particles are not efficient, and the resultant reheating temperature is too low to be consistent with the bound from big bang nucleosynthesis. We then introduce a coupling between the inflaton and the Standard Model Higgs doublet. Through the mixing with the Standard Model Higgs boson, the inflaton can decay into the Standard Model particles and the universe can be heated up with a sufficiently high reheating temperature. We have found that this happens with a sufficiently small coupling between the inflaton and the Higgs doublet and such a small has essentially no effect on our analysis for the particle mass spectrum and the inflationary predictions.