Reheating consistency condition on the classically conformal U (1) B − L Higgs inﬂation model

,


I. INTRODUCTION
The Standard Model of particle physics may be regarded as an outcome of gauging the global symmetries that were initially introduced for classification of particles.The U (1) Y hypercharge originates from the work of Nakano, Nishijima [1,2] and Gell-Mann [3], while the SU (2) L symmetry dates back to the work of Heisenberg [4] who introduced the concept of isospin.The SU (3) c quantum numbers were introduced in the 1960s for the analysis of hyperons [5,6] (see also [7]).This view may be useful for investigating a theory beyond the Standard Model.Indeed, there exists a U (1) B−L (baryon number minus lepton number) global symmetry in the Standard Model, which is usually considered accidental.Gauging the U (1) B−L , one obtains a theory beyond the Standard Model, that is the U (1) B−L gauge extended Standard Model.It is endowed with a new U (1) B−L gauge boson Z .Breaking of this gauge symmetry at low energy is accomplished by a new complex scalar field Φ, which plays the role of the Higgs boson for the U (1) B−L symmetry.Furthermore, theoretical consistency requires three chiral fermions (right-handed neutrinos) for anomaly cancellation.The minimal matter contents of the U (1) B−L extended Standard Model are thus the Standard Model particles, plus three right-handed neutrinos N R and the U (1) B−L Higgs boson Φ, as listed in Table I.
The Standard Model is known to have several issues, and interestingly, many of them find natural solutions in the U (1) B−L extention.The small but nonvanishing (left-handed) neutrino masses indicated by neutrino oscillations, for example, are naturally generated through the seesaw mechanism as the right-handed neutrinos acquire Majorana masses when the U (1) B−L symmetry is spontaneously broken.Lepton asymmetry can also be generated by the decay of the right-handed neutrinos, which may later be converted into the baryon asymmetry of the Universe in the so-called baryogenesis via leptogenesis scenario.Cosmic inflation may also be explained in the framework of the U (1) B−L -extended Standard Model, as the U (1) B−L Higgs field Φ can play the role of the inflaton, the field responsible for the dynamics of inflation.A simple, observationally viable and phenomenologically well-motivated model of cosmic inflation is constructed by allowing the Φ field to nonminimally couple to gravity.The model has been a subject of much attention and has been studied actively from various aspects [8][9][10][11][12][13].
In this article, we examine the U (1) B−L Higgs inflation model focusing on the reheating process after inflation.In inflationary cosmology, it is common to use the number of e-folds, N k , as a parameter that quantifies the expansion of the Universe during inflation, or more specifically, between the horizon exit of the scale of the cosmic microwave background (CMB) and the end of inflation.The typical range of N k is between 50 to 70; some uncertainly is usually assumed due to the modeldependent specifics of the reheating process.Once the scenario of particle cosmology is specified, however, this N k is calculable in principle.The purpose of the paper is to carry out the computation of N k in the case of the U (1) B−L Higgs inflation model.The structure of the model is that the inflationary dynamics is controlled by two parameters, the breaking scale of the U (1) B−L symmetry v BL and the U (1) B−L gauge coupling g I at low energy.We will find the relation between the e-folding number and those two parameters of the U (1) B−L -extended Standard Model, and show, by solving the renormalization group equations and eliminating the uncertainties associated with the reheating process, that the prediction for the CMB spectrum is determined by v BL and g I .At present, those parameters are largely unconstrained, neither by collider experiments or by CMB observations.We argue that these parameters will be severely constrained by near future precision measurements of the CMB spectrum, or the U (1) B−L Higgs inflation scenario will be ruled out entirely.The rest of the paper is organized as follows.We review the U (1) B−L Higgs inflation model in next section and examine the reheating process of this model in Sec.III.In Sec.IV we discuss the ranges of the model parameters v BL and g I that are of interest to us.We solve the cosmological evolution together with the renormalization group (RG) equations in Sec.V to find the prediction of the inflationary scenario.We conclude in Sec.VI with brief comments.The appendices contain supplementary mathematical details and technical notes.I.We consider the classically conformal model, that is, the scalar potential is given by where λ Φ , λ H and λ are dimensionless couplings.We assume the mixing is small, 0 < λ 1, so that the dynamics of Φ and H are separate during inflation.We may decompose the complex scalar Φ into two real scalars φ and χ as The real component φ is assumed to have a large initial value and drive inflation, that is, it plays the role of the inflaton.The χ field does not play any significant role below.Including the nonminimal coupling of Φ to gravity, the Jordan frame action for the inflaton sector is where M P = 2.435 × 10 18 GeV is the reduced Planck mass, ξ is a dimensionless parameter and is the RG-improved effective action [14,15].The running quartic coupling λ(µ) is the coupling λ Φ in (1) evaluated at the RG scale µ, and the second term is a constant that ensures the potential vanishes at the symmetry breaking global minimum φ = v BL (discussed more below).
The model is analyzed conveniently in the Einstein frame where the scalar field is minimally coupled to gravity, upon rescaling of the metric g µν → Ω(φ) g µν with The canonically normalized scalar field σ in the Einstein frame is related to φ by The scalar potential in the Einstein frame is in terms of which the slow roll parameters are defined as Under the slow roll approximation, the amplitude of the curvature perturbation at comoving scale k is which is to be compared with the measurement value1 A s at the pivot scale k.The scalar spectral index and the tensor-to-scalar ratio are expressed using the slow roll parameters as, The coupling λ(µ) is subject to the RG flow.We focus on the regime where the effects of the Yukawa couplings y i M and the running of the nonminimal coupling ξ are negligible.Then the RG equations for the self coupling λ and the U (1) B−L gauge coupling g are, at 1-loop order, We interpret the quantum corrections in the presence of nonminimal coupling as follows [17].The renormalization scale of the Jordan frame, in which the theory is defined, is given by the field φ.The renormalization scale (of mass dimension one) is then appropriately rescaled in the Einstein frame, in which measurements are made.Thus the renormalization scale µ that appears in the effective potential in the Einstein frame (8) takes the form [18] This is also the renormalization scale µ used in the RG equations ( 13) and (14).Note that the scale µ asymptote to a constant value µ → M P / √ ξ at large φ, thus the RG running slows down and stops at high energy.This behavior is in accord with the presumed UV finiteness of the theory near the Planck scale; there the metric and hence the length scale is blurred by the quantum gravity effects and above certain energy the concept of scale loses its meaning.
At low energy, the U (1) B−L gauge symmetry is broken by the Coleman-Weinberg mechanism.The symmetry breaking vacuum φ = v BL (where we live) satisfies the stationarity condition We suppose that the symmetry breaking scale is much lower than the inflationary scale (v BL M P / √ ξ).Then the renormalization scale is µ ≈ φ near φ = v BL , see (15), thus the distinction between the Einstein frame and the Jordan frame is unimportant at low energy.The condition (16) gives a relation between λ and g, where we have used the fact that in the perturbative regime the 96 g 4 term dominates the right hand side of (13).The subscript I (for IR) denotes values at the potential minimum φ = v BL .Note that λ I is negative, as it should in the symmetry breaking minimum.The mass of the Z boson and that of the inflaton are The masses of the right-handed neutrinos are given by the Majorana Yukawa coupling y i M as The offset term of the potential ( 5) is now written as The symmetry breaking scale v BL , the gauge coupling g I and the Yukawa coupling y i M at the potential minimum φ = v BL are treated as input parameters of the model.In particular, v BL and g I control the inflationary dynamics.Let the renormalization scale at the potential minimum From there the RG equation ( 14) for the gauge coupling is solved as, up to the scale µ relevant for the inflationary dynamics.Using (22), the RG equation for the self coupling (13) can be numerically integrated so that the effective potential of the inflaton ( 8) can be evaluated.The slow roll parameters are then given by ( 9), (10) as functions of φ.
To find the field value φ = φ e at which inflation ends, we use the condition that one of the slow roll parameters becomes unity, V (φ e ) = 1.The horizon exit of the CMB scale takes place at a larger value of the inflaton field φ = φ k , and there the amplitude of the curvature perturbation (11) at the pivot scale k is normalized by the observational value [16].This normalization fixes the nonminimal coupling ξ.The number of e-folds for the cosmic expansion between the horizon exit of the CMB scale and the end of inflation is In the standard slow roll paradigm of inflationary cosmology, it is a common practice to consider this e-folding number N k as a free parameter reflecting the uncertainty of the reheating process.In the next section we examine the concrete reheating process of the inflationary model based on the U (1) B−L -extended Standard Model and evaluate the e-folding number.

III. REHEATING AFTER CLASSICALLY CONFORMAL U (1)B−L HIGGS INFLATION
A salient feature of this cosmological model based on the classically conformal potential (1) is that the quartic term dominates the potential at high energy, as the mass term is generated by the Coleman-Weinberg mechanism only at the scale where the U (1) B−L symmetry is broken.Thus, at the end of inflation when the amplitude of φ is still large, the potential is essentially quartic.The symmetry breaking mass term becomes important as the oscillating amplitude of φ becomes small due to redshift.The reheating process thus proceeds stepwise: after inflation, the inflaton oscillates in the potential which is approximately quartic, and as the amplitude of the oscillations is damped by the redshift the inflaton starts to feel the presence of the mass term (19), and then starts to oscillate in the approximately quadratic potential about the symmetry breaking minimum φ = v BL .Eventually, as the Hubble expansion rate H becomes comparable to the decay rate Γ of the inflaton, the energy deposited in the inflaton is converted into the radiation of relativistic Standard Model particles and the Universe becomes thermalized.
The transition from the oscillations in the quartic-like potential to the oscillations in the quadratic-like potential is important, since the expansion rate of the Universe changes there and the prediction of the inflationary model is affected.At the transition, the inflaton that was swinging with a large amplitude fails to go over the central maximum of the double well potential.This situation is characterized by the condition that the kinetic term of the inflaton becomes comparable to the potential hight at the central maximum V E (φ = 0) = V 0 .Thus the inflaton energy density at this moment is approximately We assume that the decay of the inflaton and the ensuing thermalization of the Universe takes place after this quartic-quadratic transition.This condition is written with H the Hubble expansion rate at the transition from the quartic oscillation regime to the quadratic oscillation regime.If the decay rate Γ is larger than H , the inflaton will decay immediately after the transition and thus corresponds to the case when the condition (25) is saturated.One may also consider possible decay of the inflaton condensate into radiation during the oscillations in the quartic potential [19].We discuss this effect in Appendix B. It is found that this effect is negligible if a condition slightly weaker than (25) is satisfied.
Using the Friedman equation and ( 24), the condition ( 25) is rewritten, up to a factor of O(1), as We will see how this condition constrains the model parameters in Sec.IV.

A. The number of e-folds
We now evaluate the number of e-folds based on this picture, assuming otherwise the standard thermal history of the Universe.We denote the comoving wave number of the CMB scale by k.Then the scale factor a k and the Hubble parameter H k at the horizon exit of the CMB scale are related by k = a k H k .We write the scale factor at the end of inflation as a e , at the quartic-quadratic transition as a , at the thermalization of the Universe (end of reheating) as a th , at the matter-radiation equality as a eq , and the scale factor today as a 0 .Then one obtains an obvious relation a e a a a th a th a eq a eq a 0 where H 0 = 100 h km s −1 Mpc −1 with h = 0.674 [20] is the Hubble parameter today.The logarithm of the first factor N k ≡ ln(a k /a e ) is the e-folding number of inflation that we wish to evaluate.From the end of inflation to the quartic-quadratic transition, we may write a e /a = (ρ /ρ e ) 1/4 , where ρ is (24) and ρ e is the energy density at the end of inflation, which is roughly twice the potential energy, ρ e 2V e .We used the fact that when a scalar field oscillates in a quartic potential the Universe undergoes a radiation-dominant like expansion.Likewise, from the quartic-quadratic transition to the thermalization of the Universe we may write a /a th = (ρ th /ρ ) 1/3 , where ρ th is the energy density at thermalization and we have used the fact that when a scalar field oscillates in a quadratic potential the Universe undergoes a matter-dominant like expansion.The evaluation of the remaining factors is standard, e.g.[21,22].From the thermalization to the matter-radiation equality, entropy conservation and the Stefan-Boltzmann law give a th /a eq = (ρ eq /ρ th ) 1/4 (g eq * /g th * ) 1/12 , where ρ eq is the energy density at the matter-radiation equality and g th * , g eq * are the numbers of relativistic degrees of freedom at the thermalization and the matter-radiation equality, respectively.The factor a eq /a 0 = 1/(1 + z eq ) is the redshift of the matter-radiation equality.We use the slow roll Friedman equation to write the Hubble parameter at the time of the horizon exit of the wave number k in terms of the the potential V k , as ). Assembling all those pieces we find the e-folding number N k between the horizon exit of the comoving wave number k and the end of inflation, Apart from the uncertainly of the reheating temperature T R hidden in ρ th = π 2 g th * T 4 R /30, the e-folding number is determined by the potential (8) and can be evaluated once the dynamics of the U (1) B−L Higgs field is known 2 .In order to evaluate the reheating temperature we need to consider the decay modes of the inflaton.

B. Decay of the inflaton
Eq. (19) shows that the Z mass is heavier than the inflaton mass in the perturbative regime (g 1).Thus the decay of the inflaton into the Z boson is kinematically forbidden 3 .Also, the inflaton is a Standard Model singlet and it cannot decay through the Standard Model gauge interactions.Thus the dominant decay channel of the inflaton is through the Standard Model Higgs field.
Let us use the unitary gauge and rewrite the scalar potential (1) as We may neglect 4 quantum corrections for λ H and λ.The stationarity conditions ∂V /∂h = 0 and ∂V /∂φ = 0 at the Using (31), the last term of ( 32) is shown to be negligible, justifying the relation ( 17) that we used as the boundary conditions for the inflationary model.We also find 2 Evaluation of V k requires the value of N k (23) which needs to match (28).This can be done consistently in numerics. 3It has been pointed out in [23] that for λ 3 × 10 −4 , violent preheating into the longitudinal mode of the gauge boson may take place in the first few oscillations of the inflaton, due to spikelike features of the conformal factor Ω. This potentially leads to an issue of unitarity as the decay products have extremely high momenta ∼ √ λM P .In our model this unitarity bound corresponds to g I 0.13, which is somewhat stronger than the bound from perturbativity (see (45) and Fig. 1 below). 4For example, λ H 0.1 and its quantum corrections are O(g 4  2 /16π 2 ), which is negligible.
The Higgs mass is m h = 125.25 GeV [24].Now we may think of two separate cases, (i) when the inflaton mass is heavier than twice the Higgs mass m φ > 2m h , and (ii) when the inflaton mass is lighter than twice the Higgs mass m φ < 2m h .Let us discuss those two cases in turn.
Below in this section we consider the fields shifted about the minimum h In this case the inflaton may decay into two Higgs through the direct coupling in (30), The decay rate is where the factor of 4 is to take into account the effects of mass mixing [25][26][27][28].According to the standard perturbative picture of reheating5 , the inflaton starts to decay when the Hubble parameter becomes smaller than the decay rate Γ φ .Assuming that the thermalization is instantaneous6 , we have Γ φ H th and ρ th in ( 28) is evaluated as The reheating temperature is then found to be 2. m φ < 2m h In the second case, when the inflaton is lighter than 2m h , the process φ → hh is kinematically forbidden.If the mass range is m h < m φ < 2m h , the process φφ → hh is possible, but reheating through this process is not possible as the decay rate Γ(φφ → hh) redshifts faster than the Hubble expansion rate.When m φ < 2m h , the inflaton may instead decay through the mixing with the Higgs boson.The mass matrix of the scalars is diagonalized by rotating the fields The rotation angle is where we have used m b = 4.2 GeV, m c = 1.3 GeV and m τ = 1.8 GeV.Thus, using 2 sin 2 θ = 1 − 1/ √ 1 + tan 2 2θ and (42) the decay rate is determined by m φ and v BL .The condition for the decay is Γ φ H th .Thus the Friedman equation gives ρ th 3M 2 P Γ 2 φ and the reheating temperature is similar to (39), with Γ φ now replaced by Γ φ .

IV. CONSTRAINTS ON THE PARAMETERS
Before discussing the cosmological prediction of the inflationary model in the next section, let us summarize the constraints on the two parameters v BL and g I .
First of all, we assume that perturbative quantum field theory is valid up to the scale of inflation.As the condition of perturbativity we demand that the gauge coupling is perturbative up to the Planck scale 7 g(µ = M P ) < 1. (44) Using the solution (22) of the RG equation this condition is written, with Other conditions concern the decay of the inflaton, so let us consider the two separate cases as we did in the previous section.
1. m φ > 2m h 7 It may be somewhat more natural to consider α ≡ g 2 /4π < 1 at µ → M P / √ ξ as the criterion of perturbativity.We however use the slightly tighter condition (44) for the sake of practical convenience, as it is ξ-independent and leaves some margin from the singular regions that are numerically difficult to handle.
Using (19), the condition m φ > 2m h is written The condition (25) that the inflaton decay after the quartic-quadratic transition is written using (37) as It can be checked that ( 46) is a stronger constraint than (47) when v BL > 1.23 × 10 7 GeV.

m φ < 2m h
In this case the condition m φ < 2m h is written and the condition that the inflaton decays after the transition (25) reads where Γ φ is evaluated using (42) and (43).
Fig. 1 shows the constraints on the symmetry breaking scale v BL and the gauge coupling at low energy g I as described above.The green region is excluded by the perturbativity condition (45) and the blue region is excluded by the requirement that the decay of the inflaton takes place after the quartic-quadratic transition, (47) for the left panel and (49) for the right panel.The orange region is excluded by the condition on the inflaton mass, m φ > 2 m h for the left panel and m φ < 2 m h for the right panel.Light inflaton m φ 10 6 TeV is excluded in both cases, and there are both upper and lower bounds for g I in the m φ > 2 m h case, whereas in the m φ < 2 m h case g I is only bounded from above.
Let us also comment on the bounds that come from the collider experiments.Fig. 2 shows the bounds on v BL and g I obtained from the search for high-mass dilepton resonance by the ATLAS detector in the Large Hadron Collier [33] (139 fb −1 proton-proton collisions at a center of mass energy √ s = 13 TeV).The lower left region to the red curve is excluded.Also shown are the green and orange lines, that are the perturbativity bound and the m φ = 2m h line as in Fig. 1.Comparing Fig. 2 and Fig. 1, the bounds from the ATLAS experiments are seen to provide no further constraints on the parameter space of the inflationary model as the region is already excluded by the condition H < Γ φ (in the case of m φ > 2 m h ) or H < Γ φ (in the case of m φ < 2 m h ).

V. CMB SPECTRUM OF THE U (1)B−L HIGGS INFLATION MODEL
Let us now discuss the prediction of the cosmological model.

A. Numerical method
In order to determine the set of parameters that meet the consistency requirements and to calculate the resulting spectrum of the CMB, we employ the following procedure to solve the slow roll equation of motion and the RG equations.For a specified value of the symmetry breaking scale v BL , we choose a set of parameters N (test) k and g I so that (v BL , g I ) is within the allowed parameter region discussed in Sec.IV.Then the slow roll equation and the RG equations can be numerically integrated, using as the e-folding number of ( 23) to identify the field value φ k at the horizon exit of the CMB scale.The value of the nonminimal coupling ξ is adjusted so the the amplitude of the curvature perturbation P R matches the Planck normalization value [16] at the pivot scale.Then the cosmological evolution is determined by the set of three parameters (v BL , N (test) k , g I ), and we may evaluate the e-folding number defined by the formula (28), which, in general, differs from the value of N (test) k .We then make a scan of the parameter g I (but v BL and N (test) k kept fixed) to see if N k of (28) can be adjusted to be GeV GeV.The left panel shows the scalar spectral index ns and the tensor-to-scalar ratio r.The right panel shows ns and the U (1)B−L coupling g at the symmetry breaking minimum φ = vBL.The end points marked with • correspond to the lower bound of gI , limited by the m φ > 2m h condition or the Γ < H condition.The end points mark with • correspond to the upper bound of gI given by the perturbativity condition.The background contours are the 68% and 95% confidence level Planck+BICEP/Keck 2018 results [34] (blue), and the LiteBIRD [35] (green) and CMB-S4 [36] (red) 1-and 2-σ prospects for a fiducial model with r = 0.
the same value as is found in the range of Sec.IV, then the solution meets all consistency requirements.If this procedure fails, then that means there is no cosmological solution compatible with the reheating consistency requirement.
We carried out the parameter scan within the allowed regions of Fig. 1, and have found solutions satisfying these requirements.In the case of m φ > 2m h , for v BL 10 6 GeV there exist consistent cosmological solutions between the upper and lower bounds of g I .In contrast, when m φ < 2m h we only found consistent solutions in narrow vicinities of m φ = m h , where the mixing angle between φ and h becomes θ = π/4.While this situation may be of some phenomenological interest, it is outside of our initial assumption that the inflaton dynamics is independent of the Standard Model Higgs field during inflation, and we thus will not examine this case further.

B. Numerical results
We thus discuss the results for the m φ > 2 m h case below.Fig. 3 shows the solutions.The left panel is the prediction for the CMB spectrum, the primordial tilt n s against the tensor-to-scalar ratio r for the consistent solutions as described above.The curves indicate solutions for fixed values of v BL = 10 6 GeV to 10 12  GeV and the background contours shaded in blue are the 68% and 95% confidence level constraints of the recent Planck+BICEP/Keck 2018 constraints [34].It is seen that the prediction of the model comfortably sits inside the 68% contour, up to v BL 10 12 GeV.The right panel shows the same set of solutions on the n s -g I plane.In both panels, the endpoints marked with a filled/blank circle correspond to the lower/upper bound of g I shown on the left panel of Fig. 1.The v BL = 10 6 GeV solution is seen to be trimmed by the Γ φ < H constraint, as one can see by comparing with Fig. 1, left panel.
In Fig. 3, the prospect constraint contours by the Lite-BIRD and CMB-S4, for a r = 0 fiducial model are shown in green and red.The prediction of the cosmological model studied here is clearly outside the 2-σ contours, and thus would be strongly disfavored if those projects bring null results.If, on the other hand, the tensor mode is detected, the measurements of the CMB spectrum would give significant constraints on the parameter space of the U (1) B−L Higgs inflation model.

VI. FINAL REMARKS
We have examined the reheating process of the inflationary scenario based on the U (1) B−L extension of the Standard Model, and formulated the condition of consistency in terms of the number of e-folds.We then solved the equation of the inflationary dynamics along with the RG equations to identify solutions that meet these requirements.The results show that the predictions of the CMB spectrum are in excellent agreement with current observational constraints.It is also suggested that the proposed model could be tested by future experiments, such as LiteBIRD and CMB-S4.Our aim was to address the previously overlooked aspects of model construction and to provide a clearer prediction for cosmological observables by incorporating the consistency condition from the reheating process.
The primary focus of this paper has been the analysis of a simple inflationary model, which is characterized by two key parameters: the U (1) B−L breaking scale (v BL ) and the U (1) B−L gauge coupling (g I ) at low energy.The U (1) B−L extension of the Standard Model is a well-motivated theory beyond the Standard Model, and this example may be considered as one of the best candidate cosmological scenarios based on particle phenomenology.Clearly, our analysis can be extended to more involved cosmological models.For instance, the U (1) B−L model can be extended to the U (1) X model that allows for the mixing of the U (1) B−L and U (1) Y gauge symmetries without violating the anomaly cancellation condition, as described for example in [13].Additionally, inflationary models based on supersymmetric extensions of the Standard Model, such as those discussed in [37,38], may also be worthy of exploration.As upcoming observational cosmology projects are poised to bring new results in the near future, particularly with regards to the CMB B-model polarization, it is a promising time to re-evaluate the reheating dynamics of these inflationary scenarios.represents the energy transmission into radiation.Thus the depletion of the inflaton energy by the decay into radiation is negligible if where t is the time when the inflaton starts to oscillate in the quadratic potential.Using (B7) and evaluating φ o ∼ v BL at t = t , the condition (B13) is equivalent to 1 4π Now using (31), (33)  .

(B16)
This is seen to be a slightly weaker condition than (47).

Scattering picture
Instead of the decay φ → H † H, one may alternatively consider the scattering process φφ → H † H.The initial φ's are assumed to be condensed, at rest with energy ω φ .The number density of the inflaton quanta φ must then obey the Boltzmann equation with The energy of the inflaton quanta ω φ may be evaluated as We shall show that (B17) is equivalent to (B9), up to a numerical factor.We first note that n φ = ρ φ /ω φ , and that ω φ redshifts as ω φ ∝ φ o ∝ 1/a, so that coinciding with (B18) up to a factor of ∼ 8, that may be attributed to the approximations we have used.
II. HIGGS INFLATION IN THEU (1)B−L-EXTENDED STANDARD MODEL We consider the minimal U (1) B−L extension of the Standard Model, with the gauge group SU (3) c × SU (2) L ×U (1) Y ×U (1) B−L .The particle contents are the Standard Model particles supplemented by three generations of singlet leptons (the right-handed neutrinos) and a complex scalar with the U (1) B−L charge 2, as listed in Table

2 , y c sin θ/ √ 2 ,
which is small in general, apart from the accidental narrow region of m φ m h .Thus the field h is almost h, and φ is almost φ in generic cases.This almost-inflaton φ couple to the bb, cc, τ τ of the Standard Model with the Yukawa couplings y b sin θ/ √ y τ sin θ/ √ 2 respectively, where y b , y c , y τ are the Standard Model Yukawa couplings for b, c and τ .The decay rate of φ into the Standard Model particles is then

FIG. 2 .
FIG. 2. Constraints from the ATLAS experiments[33] at Run 2 of the Large Hadron Collider (the center-of-mass energy √ s = 13 TeV and integrated luminosity 139 fb −1 ), recast into bounds on the parameters vBL and gI of the U (1)B−Lextended Standard Model.The region lower left to the red curve is excluded.The green and orange lines are respectively the perturbativity limit and the m φ = 2m h line as in Fig.1.

FIG. 3 .
FIG.3.The prediction of the U (1)B−L Higgs inflation model, with the requirement of the reheating consistency taken into account.This is the m φ > 2 m h case of Fig.1.The value of vBL is varied as 10 6 , 10 7 , 10 8 , 10 9 , 10 10 , 10 11 , 10 12 GeV.The left panel shows the scalar spectral index ns and the tensor-to-scalar ratio r.The right panel shows ns and the U (1)B−L coupling g at the symmetry breaking minimum φ = vBL.The end points marked with • correspond to the lower bound of gI , limited by the m φ > 2m h condition or the Γ < H condition.The end points mark with • correspond to the upper bound of gI given by the perturbativity condition.The background contours are the 68% and 95% confidence level Planck+BICEP/Keck 2018 results[34] (blue), and the LiteBIRD[35] (green) and CMB-S4[36] (red) 1-and 2-σ prospects for a fiducial model with r = 0.

TABLE I .
Representations and charges of the particle con- tents in the U (1)B−L-extended Standard Model.The subscripts L/R are the chiralities and the index i = 1, 2, 3 indicates the generations of the fermions. ) I m ϕ <2m h case H ★ <Γ ϕ ∼ m ϕ >2m h Nonperturbative FIG.1.The constraints on the parameters (vBL, gI ), by the conditions that (i) the inflaton decays after it starts to oscillate about the minimum φ = vBL, (ii) the inflaton is heavier (left panel) or lighter (right panel) than 2m h , and (iii) the U (1)B−L gauge coupling is perturbative up until the scale of inflation.The blank regions are unconstrained.
The decay amplitude for φ → H † H is As the universe expands like radiation dominated a ∝ √ t in this regime, the inflaton redshifts as φ o ∼ φ e t e /t, where t e and φ e are the cosmic time and the background inflaton value at the end of inflation.Using H(t = t e ) = 1/2t e and the slow roll equation of motion, we findThe energy density of the inflaton ρ φ and that of the radiation ρ rad evolve according to The total energy density ρ total = ρ φ + ρ rad thus evolves as ρ total ∝ a −4 .The equation (B9) is solved asρ φ (t) = ρ φ (t e ) t e tThe factor (t e /t) 2 is due to the dilution by the cosmic expansion and the exponential factor with and evaluating λ in our model of the Coleman-Weinberg symmetry breaking as