Inflection point inflation and reheating

We revisit the inflection point inflation with an extended discussion to large field values and consider the reheating effects on the inflationary predictions. Parametrizing the reheating dynamics in terms of the reheating temperature and the equation of state during reheating, we show how the observationally favored parameter space of inflection point inflation is affected by reheating dynamics. Consequently, we apply the general results to the inflation models with non-minimal coupling, such as the SM Higgs inflation and 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 inflation.


Introduction
Cosmic inflation is the main paradigm for the early Universe that solves horizon, homogeneity, isotropy and relic problems in Standard Big Bang Cosmology and it predicts the anisotropies of Cosmic Microwave Background (CMB) and seeds the structure formation in the early Universe. Inflation is based on the initial stage of an exponential expansion in the presence of a slow-rolling scalar field with nonzero vacuum energy. Inflaton has been assumed to be a SM singlet and have an almost flat potential during inflation. Since the temperature of the Universe became almost zero after inflation ended, it is necessary to reheat the Universe by introducing the inflaton couplings to the SM particles. As the inflaton is a SM singlet, the reheating mechanism [1][2][3][4][5] is mostly unknown and it could have been more complicated than we naively think. In most cases, taking a range of number of efoldings required for inflation between N * = 50 and N * = 60 is attributed to the uncertainties in reheating dynamics. The recent results from Planck satellite constrain the spectral index with unprecedented precision [6] and the tensor-to-scalar ratio is best constrained by Planck + BICEP2 + Keck [7]. Therefore, it becomes more important to pin a e-mail: hminlee@cau.ac.kr down the uncertainties due to reheating dynamics in each class of inflation models [8][9][10][11].
We consider the general predictions for CMB anisotropies and the effects of reheating dynamics on them in the context of inflection point inflation models [12][13][14][15][16][17][18][19][20][21]. Unlike small or large inflation models, the inflection point inflation is based on specific values of the inflaton for which the inflaton potential has a zero curvature and a small slope while it is dominated by a constant vacuum energy during inflation. We keep up to the cubic terms near the inflection point by assuming that the higher order terms vanish due to small field excursion during inflation or for a symmetry reason of the full inflaton potential. It turns out that the effective potential for inflection point inflation is a good approximation in most cases, as far as the full potential shows such a specific point of the field space.
We divide our discussion into two cases, one with small field excursion and one with large field excursion. The effects of reheating dynamics is translated into the change of the number of efoldings required for solving the horizon problem. Depending on the reheating temperature and the equation of state during reheating, we constrain the parameter space of the effective inflaton couplings from the spectral index and the tensor-to-scalar ratio. Then, we also apply the general results to concrete models for inflection point inflation such as SM Higgs inflation and B − L Higgs inflation with non-minimal coupling where the renormalization group (RG) running effects become important for inflation. Both inflation models belong to the case with a large field excursion of the canonical inflaton. As a consequence, we identify the effective couplings of the inflaton potential with the fundamental parameters of the given model and show the parameter space where the tensor-to-scalar ratio is sizable. Depending on the models, the reheating effects become important for inflationary predictions, in particular, for the B − L Higgs inflation with a relatively small inflaton mass and a small Higgs-portal coupling.
The paper is organized as follows. We first present the general predictions of inflation observables in inflection point

Inflection point inflation
In this section, we summarize the general results for inflection point inflation. Some results on inflection point inflation have been obtained in the literature in a particular limit with small fields, for instance, in the context of warped D-brane inflation models [12,13]. We extend those results to the case with a large field excursion for which the inflationary predictions are different from the case with a small field excursion. There might be a concern on the initial condition in inflection point inflation because the inflaton should not have a large kinetic energy near the inflection point for generic initial conditions given away from the inflection point [22,23]. However, there might be dissipative effects of the inflaton alleviating the initial condition problem, as they could provide friction terms to set the initial condition for inflection point inflation [24].
We take the inflaton potential near the inflection point φ = φ 0 as where V 0 , λ 1 , λ 3 are constant parameters depending on the models. A sketch of the inflaton potential in the above form is given in Fig. 1. Henceforth, we work in units of the Planck mass set to M P = 1 and we take all the constant parameters above to be positive, without loss of generality. 1 We assumed that higher order terms are suppressed due to a small field excursion during inflation or for some symmetry reason protecting the potential against them. 1 For λ 1 < 0, the inflection point is near a false vacuum so there is a problem with graceful exit.
Then, the slow-roll parameters are computed as , the number of efoldings from φ to the end of inflation is Choosing λ 1 /λ 3 1, we can use the approximate form for the number of efoldings as In this case, it is necessary to take N max 1 for a large number of efoldings. By taking φ = φ * at horizon exit with |φ * − φ 0 | 1, the number of efoldings from horizon exit to the end of inflation is where η * ≡ η(φ * ) and η end ≡ η(φ end ). We note that for η * < 2π/N max and η end ∼ −1, the number of efoldings during inflation is N * N max /2. Consequently, from η * ≈ λ 3 V 0 (φ * − φ 0 ), we can solve Eq. (6) for φ * as follows, Plugging Eq. (7) into Eqs. (2) and (3) during inflation, we obtain the slow-roll parameters at horizon exit as follows, * = 1 2 On the other hand, inflation ends when end = 1. Thus, from end = 1 with Eq. (2), we can determine the inflaton field value at the end of inflation as For λ 1 /λ 3 1, the inflaton field excursion in inflection point inflation is determined mainly by the value of λ 3 /V 0 . From the inflation field value at the end of inflation in Eq. (11), the inflaton makes a trans-Planckian excursion during inflation for λ 3 /V 0 1 while becoming sub-Planckian for λ 3 /V 0 1. Therefore, from n s = 1 − 6 * + 2η * and r = 16 * , the spectral index and the tensor-to-scalar ratio becomes By using the above results, we find that the spectral index is correlated to the tensor-to-scalar ratio in this model as follows, Thus, for a fixed value of tan in Eq. (7), we can correlate between the spectral index and the tensor-to-scalar ratio. We note that the measured spectral index and the bound on the tensor-to-scalar ratio, are n s = 0.9652 ± 0.0047 and r < 0.10 at 95 % C.L., respectively, from Planck TT, TE, EE + low P [6]. Moreover, from the CMB normalization, A s = 1 Then, the sufficient number of efoldings and the spectral index constrain the following combination of the model parameters, Therefore, we can express the model parameters in terms of the tensor-to-scalar ratio r , N max and cos as Furthermore, using Eqs. (19) and (20), the useful discriminator between small and large field inflations, λ 3 /V 0 , is given by Therefore, for sub-Planckian inflation field excursions with λ 3 /V 0 1, the tensor-to-scalar ratio becomes about 10 −5 or less. On the other hand, for trans-Planckian field excursions, the tensor-to-scalar ratio as large as r = 0.1 can be achieved for λ 3 /V 0 ∼ 10 −2 or |φ end − φ * | ∼ 10. From Eqs. (20) and (21), we also find the ratio of λ 1 and λ 3 couplings as Then, for N max ∼ 100 and r cos 4 0.01, we get λ 1 λ 3 . In Fig. 2, the parameter space for λ 3 /V 0 vs N max is shown, under the condition that the spectral index is obtained within 1σ and 2σ of the Planck data, depending on the number of efoldings, N = 50 or 60. For given ratios of parameters, λ 3 /V 0 and λ 1 /V 0 , the CMB normalization (17) is achieved along a line in the parameter space in Fig. 2. Representative values of λ 1 , such as λ 1 = 10 −10 , 10 −13 , 10 −16 , 10 −19 , are overlaid in solid, dashed, dot-dashed and dotted red lines in the same figure. The boundary dividing the sub-and trans-Planckian inflaton excursions appears around r 10 −5 or 10 −6 . There is an isolated region for small λ 3 /V 0 and a relatively small N max in the left lower corner. In this region, however, our assumption that |φ * − φ 0 | < 1, namely, the fact that the inflation starts taking place near the inflection point is violated, so we don't consider the region any more.

Small-field case
When λ 3 /V 0 1, we have discussed that the field excursion during inflation is sub-Planckian. In this case, we can trust , satisfying the spectral index measured by Planck within 1σ (green) and 2σ (yellow) for the number of efoldings, N = 50 on left and N = 60 on right. For a given value of λ 1 , the CMB normalization is satisfied only along the red lines. Some values of tensor-to-scalar ratio, r = 0.1, 10 −3 , 10 −5 , 10 −7 , are shown in dotted, dot-dashed, dashed and solid lines (from bottom to top), respectively. The region with inflation field value at horizon exit given by |φ * − φ 0 | > 1 is shown in gray the perturbative expansion of the inflaton potential near the inflection point with higher order terms of (φ − φ 0 ). In this case, the model parameters have to satisfy √ λ 1 λ 3 V 0 λ 3 , which corresponds to the hierarchy, λ 1 V 0 λ 3 . Moreover, from Eq. (22) with N max ∼ 100 or as shown in the region with |φ * − φ 0 | < 1 in Fig. 2, we find that r 10 −5 or 10 −6 in the small-field case. In this case, since * |η * |, we can ignore the contribution from * in the spectral index. Then, from Eq. (14) with η end = −1, the spectral index becomes [12,13] Therefore, for N * = 50(60), we have derived the bounds on N max as 128 < N max < 141(162 < N max < 187) within 1σ from the Planck central value. This result can be seen clearly in Fig. 2 for the region with λ 3 /V 0 1.

Large-field case
When λ 3 /V 0 1, the inflaton can make a field excursion towards trans-Planckian values during inflation. Then, there is a concern about the higher order terms in the perturbative expansion near the inflection point. However, if the full inflaton potential is known to have higher order terms suppressed near the inflection point, the case with trans-Planckian field values remains an interesting possibility. In this case, the spectral index is given by Eq. (14). Moreover, from Eq. (22), a sizable tensor-to-scalar ratio can be obtained. For instance, for r ∼ 0.10 and N max ∼ 100, the model parameters have to satisfy a different hierarchy, λ 3 λ 1 V 0 . As shown in Fig. 2, for λ 3 /V 0 1 for which the inflaton field makes a trans-Planckian field excursion from the inflection point during inflation, the tensor-to-scalar ratio should be greater than r ∼ 10 −5 and the linear term should be larger than λ 1 ∼ 10 −15 in the region where the observed spectral index is obtained.

Reheating and number of efoldings
We parametrize the reheating dynamics in terms of the reheating temperature, at which the inflaton energy is completely transferred into the radiation of the SM particles, and the equation of state during reheating, w = p/ρ [25], which is restricted to − 1 3 ≤ w ≤ 1 for no more acceleration,ä ≤ 0. When the form of the inflaton potential is known after inflation, the equation of state might be in principle determined from the averaged energy density and pressure of the inflaton during reheating. However, it is known that a small inflaton coupling to spectator fields changes the equation of state by a large amount in a sufficiently short time scale during the process of preheating [25]. Moreover, in some cases, the precise form of the inflaton potential near the oscillation does not affect the inflection point inflation, e.g., in the case of inflation models with non-minimal coupling as will be discussed in a later section. Therefore, we treat the equation of state to be a free parameter.
Assuming that there is no entropy change between the reheating and the present, we obtain the number of efoldings between the horizon exit and the end of inflation required to solve the horizon problem up to the scale of k, as follows, where g * (T rh ) is the number of relativisitic species in thermal bath at reheating temperature T rh , and N rh is the number of efoldings during reheating, given by Taking V end V * and choosing g * (T rh ) 100 for T rh 100 GeV, the formula for the number of efoldings becomes We note that the last term depends on the inflaton potential and the second term depends on both inflation models and reheating dynamics. In Fig. 3, the dependence of the number of efoldings on the reheating temperature and the equation of state during inflation is given.
In the case of an instantaneous reheating for which the reheating temperature is given by T rh = (45V end /π 2 g * ) 1/4 (8.5 × 10 14 GeV)(r/10 −5 ) 1/4 , we would get N rh = 0, so the number of efoldings depends only on the inflaton potential. In this case, the number of efoldings for the inflection point inflation is given by N * = 53.7 + 1 4 ln(r/10 −5 ). If the reheating is delayed due to a late-time decay of the inflaton, depending on whether w > 1 3 or not, the contribution coming from reheating dynamics can increase or decrease the number of efoldings during inflation. The reason is the following. For w > 1 3 , the horizon radius of the Universe at recombination would have been larger than in the case of instantaneous reheating. Therefore, the largest scale of cosmological interest must have left the horizon earlier so that the number of efoldings during inflation gets larger. On the other hand, for w < 1

Examples of inflection point inflation
In this section, we illustrate some examples for inflection point inflation in the context of the SM Higgs inflation and the B − L Higgs inflation. Using the general results for the form of the effective potential near the inflection point in the previous section, we obtain the approximate formulas for the coefficients in each case and discuss the relation between the effective couplings and the fundamental parameters of the model.

SM Higgs inflation
First we employ the general results for inflection point inflation to the case of the SM Higgs inflation minimally coupled to gravity and show that the number of efoldings in this case is insufficient. Then, we continue our discussion for the SM Higgs inflation with a non-minimal coupling to gravity and comment on the effects of reheating in this case.

The case without non-minimal coupling
We discuss the SM Higgs inflation at the inflection point by using the general discussion in the previous section. The RG-improved Higgs potential takes the following form, where λ(h) is the running quartic coupling evaluated at the Higgs field value [26,27]. Here, near the criticality where the beta function for the effective Higgs quartic coupling vanishes, namely, β = dλ/d ln h = 0, the effective Higgs quartic coupling is given by where b is the two-loop coefficient given by b = 0.4/(4π) 4 and the Higgs field value at criticality is h cr ∼ 1. Now we expand the Higgs potential near the inflection point h = h 0 as follows, with where use is made of the inflection point condition Thus, the inflection point is located at h 0 < h cr . Requiring λ cr < 25 144 b for real λ 1 , λ 3 , and λ cr > b 16 for positive λ 1 , we can show that there is no false vacuum. We also note that V 0 is always positive. Taking λ cr ≈ b 16 ≈ 10 −6 for λ 1 V 0 , we obtain the other parameters of the potential as where a depends on the tuning between λ cr and b 16 and is matched to a ≈ cos by the CMB normalization in Eq. (17), and the inflection point as In this case, however, as the Planck bound on the tensorto-scalar ratio leads to h 0 (b/1.0 × 10 −5 ) 1/4 0.32, the number of efoldings during inflation is bounded by Taking b = 1.0 × 10 −5 in the SM, the number of efoldings is too small.

The case with non-minimal coupling
We introduce a non-minimal coupling of the SM Higgs to gravity. The tree-level Jordan-frame action for the inflaton h with non-minimal coupling is given by where ξ is the non-minimal gravity coupling. Then, by making a Weyl transformation of the metric with g μν = g E μν / 2 with = 1 + ξ h 2 , we obtain the Einstein-frame action as where the canonical inflaton field is related to the Jordan frame Higgs by and the inflaton potential in Einstein frame is given by In the original Higgs inflation with tree level λ, the COBE normalization requires a large non-minimal coupling, ξ ∼ 10 4 [28]. Therefore, in this case, the unitarity scale in the vacuum is not far from the Hubble scale during inflation, making the semi-classical approximation of the Higgs inflation questionable [29][30][31][32]. Nonetheless, it is known that the unitarity scale becomes background-field dependent so it becomes larger than the Hubble scale during inflation [33][34][35]. The Higgs inflation with non-minimal coupling is shown to be identical to the R 2 inflation or its scalar-tensor dual at the classical level [36,37], although there are distinctions at the loop level due to the would-be Goldstone bosons in the Higgs inflation case [29][30][31][32].
At the loop-level, the effective inflaton quartic coupling can be identified from the one-loop Higgs potential with rescaled inflaton-dependent masses of particles and it turns out to be a function of ϕ [26,[38][39][40][41]. Near the inflection point where the one-loop beta function vanishes, the loop corrections become more important and the effective inflaton quartic coupling is dominated by two-loop beta functions as follows, We note that for the SM Higgs boson, the two-loop beta function is given by b = 1.0 × 10 −5 ≡ b SM [26] whereas λ cr and ϕ cr depend on the precise values of top quark and gauge couplings in the low energy. The Higgs inflation near the inflection point was discussed in the literature [42][43][44][45] and it was generalized to the singlet scalar inflation with an extra U (1) gauge theory near the inflection point [46]. In this work, we revisit the Higgs inflation near criticality (meaning a zero quartic coupling and its vanishing beta function) from the point of view of inflection point inflation. Now we are making an expansion of the effective inflaton potential around the inflection point χ = χ 0 for the canonical Higgs (equivalently to ϕ = ϕ 0 ) as follows, where Here, the inflection point is determined from Henceforth, we assume that ξ 1. Then, from Eq. (42), taking ξ h 2 1 during inflation, the canonical field is related to the inflaton field in Jordan frame by χ . In this case, we also get ϕ 1 √ ξ 1 − 1 2ξ h 2 , resulting in dϕ dh 1 ξ 3/2 h 3 . Therefore, we get the derivatives of the potential as follows, Consequently, from Eqs. (46)-(48), the model parameters can be computed at the inflection point χ 0 or ϕ 0 1/ √ ξ as where we imposed the inflection point condition, to be small, we obtain λ cr b 16 and ϕ 0 e − 1 4 ϕ cr 1/ √ ξ , leading to the non-minimal coupling, ξ e 1 2 ϕ −2 cr . Then, for h 0 = c √ ξ with c 1, the effective couplings near the inflection point are determined as Here, we note that a depends on the tuning between λ cr and b 16 , determined by the low-energy parameters in the RG equation for the quartic coupling, and it is matched to a ≈ cos by the CMB normalization (17) for given ratios of inflaton potential parameters and N * from Eq. (6). As compared to the case with minimally coupled Higgs boson in Eqs. (35)-(37), the effective couplings of the inflection point inflation can be suppressed for large ξ and/or c, meaning that the range of the approximate inflection point is extended. We have checked that the higher order terms in the inflaton potential near the inflection point are suppressed by a sizable ξ or a large inflation field value in the relevant parameter space. Therefore, the number of efoldings during inflation is bounded by which is enhanced for a large non-minimal coupling and/or a large inflation field value at the inflection point, as compared to the case with minimally coupled Higgs boson in Eq. (39). From the Planck bound on the tensor-to-scalar ratio, we get the upper bound on the non-minimal coupling as ξ 9.6 (b/1.0 × 10 −5 ) 1/2 . Consequently, a relatively small non-minimal coupling is allowed in inflection point inflection as compared to the classical Higgs inflation. In this case, the field-dependent cutoff given by M P / √ ξ is much larger than the Hubble scale of order √ λ(ϕ I )M P /ξ , which is suppressed for a small λ(ϕ I ) for inflation field values ϕ I .
Imposing N max ∼ 120 in Eq. (60) with the bound on tensor-to-scalar ratio, ξ 9.6 (b/1.0 × 10 −5 ) 1/2 , we get the upper bound on the inflaton field value at the inflection point as c 2.6 a 1/3 , independent of the value of b.
On the other hand, using V 0 λ 3 0.2c 6 from Eqs. (53) and (55), we find a trans-Planckian field excursion of the canonical inflation field during inflation for c 1.3. Therefore, in the SM Higgs inflation with non-minimal coupling, the inflation field values near the inflection point are necessarily trans-Planckian for a = O(1) in terms of the canonical inflation field χ . Nonetheless, the inflaton field in the Jordan frame can be sub-Planckian, because h 0 = c/ √ ξ and/or a larger ξ from the bound on the tensor-to-scalar ratio.
In Fig. 5, employing the approximate potential up to cubic terms near the inflection point and taking the twoloop beta function coefficient to b = 1.0 × 10 −5 , we show the parameter space for ξ and c for a chosen tan = 0, −0.1, −0.4, −0.7 in red lines with N * = 50, 60 while the spectral index favored by Planck 1σ and 2σ bounds is shown in green and yellow regions, respectively. Furthermore, in the same figures, the lines with tensor-to-scalar ratio, r = 10 −1 , 10 −1.5 and r = 10 −1 , 10 −2 , 10 −3 are shown in blue lines, in upper and lower panels, respectively, and the region with inflation field value given by |φ * − φ 0 | > 1 at horizon exit is shown in gray too. Several values of h 0 are shown along the black lines in each plot, thus telling us that the Higgs field value at inflection point in Jordan frame is sub-Planckian in most of the parameter space. For a chosen value of tan , the λ 1 coupling is accordingly fixed by CMB normalization.
In the case with a small tan in the upper panel, the measured spectral index favors the non-minimal coupling ξ between 10 and 13 on left (13 and 17 on right) and the inflection point c between 2.5 and 2.9 in both cases. Moreover, in this case, the tensor-to-scalar ratio tends to be large, saturating the upper bound from Planck.
On the other hand, in the case with a large tan in the lower panel, where the region with |φ * − φ 0 | > 1 at horizon exit gets larger, the favored region is given by the nonminimal coupling ξ larger than 20 (30) and the inflection point c smaller than 2.8(2.6) on left (right). Moreover, in this case, the tensor-to-scalar ratio becomes smaller than the Planck bound in most of the parameter space favored by the measured spectral index. We note that | tan | is favored by the spectral index to be smaller than about 0.7.
In Fig. 6, we also depict the spectral index as a function of the non-minimal coupling in red lines, for several choices of inflection point h 0 and tan to (0.2, −0.4), (0.2, −0.7), (0.7, −0.1), (0.5, −0.4), in the plots clockwise. In each plot, the spectral index is bounded to blue line for the number of efoldings between N * = 50 and 60 and the spectral index favored by Planck 1σ and 2σ bounds is shown in green and yellow regions, respectively. Two-loop beta function coefficient is chosen to b = 1.0 × 10 −5 . We show the obtained tensor-to-scalar ratio in black line in the same figures as well. The region with inflation field value taking |φ * − φ 0 | > 1 at horizon exit is shown in gray line in the case with (0.5, −0.4).
We note that in the region with a relatively small nonminimal coupling in the lower panel of Fig. 6, satisfying the observed spectral index, the tensor-to-scalar ratio tends to be large and sensitive to the choice of the inflection point. The larger the value of h 0 , the smaller value of ξ the observed spectral index is obtained for, and the larger the tensor-toscalar ratio. Now we remark on the dependence of the inflationary predictions on the reheating dynamics. The SM Higgs couples to the rest of the SM particles strongly, so the reheating temperature is high and expected to be of similar order to the one for the classical SM Higgs inflation, which is 3 − 15 × 10 13 GeV [47,48]. Therefore, the number of efoldings during reheating could not be sizable, so the inflationary predictions are little dependent on the reheating dynamics. In this case, the infla-tionary predictions with instantaneous reheating or w = 1 3 are reliable.

B − L Higgs inflation
The U (1) B−L symmetry is an anomaly-free gauge theory that is introduced as a natural extension of the SM with three Right-Handed (RH) neutrinos for generating neutrino masses via see-saw mechanism [49,50]. In the U (1) B−L extension of the SM, we consider the effective potential for the B − L Higgs φ which carries +2 charge under the U (1) B−L and has a nonzero non-minimal coupling to gravity. Then, we identify the B − L Higgs field as the inflaton. The Coleman-Weinberg inflation in the context of the B − L Higgs inflation was discussed [51,52] and the B − L running inflation with non-minimal coupling was also recently studied [53][54][55][56].
Near the critical point φ = φ cr , we choose the B − L Higgs quartic coupling to λ φ (φ cr ) ≈ 0 and one-loop beta function to β (1) λ φ (φ cr ) ≈ 0. In this case, we need to choose the quartic coupling λ φ for the B − L Higgs to be small at low energies as for the SM Higgs. Then, ignoring λ φ in Eq. (A. 5) in Appendix A, we can obtain the relation between the RH neutrino Yukawa coupling y N and the B − L gauge coupling g B−L at the critical point as Since λ φ is small at low energies, the B − L Higgs tends to be light 2 so it decays dominantly into a pair of the SM Higgs bosons, reheating the SM sector. Then, the reheating temperature is given by 32π m φ for m φ m h with m h being the SM Higgs mass. Thus, imposing the bound on the mixing angle between the B − L Higgs and the SM Higgs from the Higgs signal strength, such as | sin θ | ≈ 2 The B − L symmetry can be broken with a tree-level negative mass squared for the B − L Higgs or by the Coleman-Weinberg mechanism [49,50]. In the latter case, the B − L Higgs without non-minimal coupling was discussed in view of Planck data as an inflaton candidate [51]. [58], we obtain the lower bound on the reheating temperature as Suppose that the Higgs-portal coupling λ H φ is very small in the UV so it is generated radiatively in the presence of the gauge kinetic mixing [49,50]. In this case, when the Higgs mass parameter is solely from the Higgs-portal coupling after the B−L breaking, we get |λ hφ | = [49,50], which satisfies the limit from the Higgs signal strength for v B−L m h . Therefore, the reheating temperature can be much smaller as follows, Therefore, depending on the equation of state during reheating, the reheating dynamics can affect much the inflationary predictions of the B − L Higgs inflation. In Fig. 7, we present the parameter space for the inflection point c = h 0 √ ξ vs the non-minimal coupling ξ in red lines, for a fixed tan = 0, −0.1, −0.4, from top to bottom panels. In each plot, we take the equation of state w during inflation between − 1 3 , 0, 2 3 , 1, in dotted, dotdashed, dashed and solid red lines, and the Planck 1σ and 2σ bounds on the spectral index are shown in green and yellow regions, respectively. We vary the reheating temperature to be T rh = 100, 10 7 , 10 15 GeV from left to right plots in each panel. The lower the reheating temperatures, the smaller the required number of efoldings at horizon exit becomes for w < 1 3 , opening the wider parameter space for ξ and c in the region with smaller values of c. Several values of inflection point h 0 are shown in each plot, showing that the inflaton field values in Jordan frame are sub-Planckian in most of parameter space. We also show the tensor-toscalar ratio with r = 10 −1 , 10 −1.5 , r = 10 −1 , 10 −2 , and r = 10 −1 , 10 −2 , 10 −3 , in blue lines, in the upper, middle and lower panels, respectively.
For a small tan as in the upper panel of Fig. 7, the tensorto-scalar ratio tends to be large, saturating the Planck bound, and the required non-minimal coupling is insensitive to the change of the equation of state, in the region favored by the measured spectral index. On the other hand, for a sizable tan as in the lower panel, the tensor-to-scalar ratio tends to be small, but the predictions for the spectral index and tensorto-scalar ratio are sensitive to the equation of state in the case of low reheating temperature. In this case, the required nonminimal coupling varies in a wide range of values.
In Fig. 8, we also show that the spectral index and the tensor-to-scalar ratio as a function of the non-minimal coupling, for a set of inflection point h 0 and tan to  The spectral index varies along the black line, depending on the equation of state, − 1 3 ≤ w < 0, 0 ≤ w < 2 3 , 2 3 ≤ w < 1, in red, blue and orange lines in each plot. The region with inflation field value taking |φ * − φ 0 | > 1 at horizon exit is shown in gray line in the lower panel. As the equation of state during reheating varies from w = − 1 3 to w = 1, the spectral index increases and falls into the region favored by Planck while the tensor-to-scalar ratio decreases. On the other hand, as the reheating temperature decreases, the parameter 3 ≤ w < 0, 0 ≤ w < 2 3 , 2 3 ≤ w < 1, in red, blue and orange lines, in each plot. Others are the same as in Fig. 6 space for small w widens, but it tends to be disfavored by the measured spectral index in Fig. 8. However, the effect of the decreasing equation of state at low reheating temperature can be compensated by taking a smaller inflection point, as shown in Fig. 7, such that the consistent spectral index is obtained.

Conclusions
We have considered the inflection point inflation as one of the slow-roll inflation models that might be insensitive to trans-Planckian problems usually encountered in large field inflation models. We showed that the inflection point inflation predicting a small tensor-to-scalar ratio is sensitive to the reheating temperature and the equation of state during inflation, whereas the case with a large tensor-to-scalar ratio tends to be robust and insensitive to the reheating dynamics. From the examples of the SM Higgs inflation and the B − L Higgs inflation where the reheating temperature is determined, we presented the inflationary predictions for the spectral index and the tensor-to-scalar ratio, depending on the reheating dynamics. In particular, in the B − L Higgs inflation, when the light inflaton has a small Higgs-portal coupling, the reheating temperature is small, so the reheating dynamics could affect the inflationary predictions much if the equation of state deviates from radiation.

Appendix: Two-loop renormalization group equations in the U(1) B−L extension
In the U (1) B−L extension of the SM, the quartic couplings in the scalar potential are given by V quartic = λ H |H | 4 + λ φ |φ| 4 + λ H φ |φ| 2 |H | 2 , (A.1) and the Yukawa couplings for the RH neutrinos are L Y = − 1 2 y N φ N c N +h.c. Then, when three RH Yukawa couplings are equal and the Higgs-portal coupling λ H φ and the gauge kinetic mixing between the SM hypercharge and U (1) B−L gauge bosons are small, the RG equations for the B − L Higgs quartic coupling, the RH neutrino Yukawa couplings and the B − L gauge coupling are given at two-loop order [49,50,57] by dλ φ d ln μ = β (1) λ φ + β (2) λ φ , (A.2) dy N d ln μ = β (1) y N + β (2) y N , where one-loop and two-loop beta functions for the B − L Higgs quartic coupling are, respectively, with g 1,2,3 being the SM gauge couplings and y t is the top Yukawa coupling.