Constraints on Inflaton Higgs Field Couplings

According to the best-fit parameters of the Standard Model, the Higgs field's potential reaches a maximum at a field value $h \sim 10^{10-11}$ GeV and then turns over to negative values. During reheating after inflation, resonance between the inflaton and the Higgs can cause the Higgs to fluctuate past this maximum and run down the dangerous side of the potential if these fields couple too strongly. In this paper, we place constraints on the inflaton-Higgs couplings such that the probability of the Higgs entering the unstable regime during reheating is small. To do so, the equations of motion are approximately solved semi-analytically, then solved fully numerically. Next the growth in variance is used to determine the parameter space for $\kappa$ and $\alpha$, the coupling coefficients for inflaton-Higgs cubic and quartic interactions, respectively. We find the upper bounds of $\kappa<1.6 \times 10^{-5} m_\phi \sim 2.2 \times 10^8$ GeV and $\alpha<10^{-8}$ to allow the Higgs to remain stable in most Hubble patches during reheating, and we also find the full two parameter joint constraints. We find a corresponding bound on the reheat temperature of $T_\text{reh} \lesssim 9.2 \times 10^9$ GeV. Additionally, de Sitter temperature fluctuations during inflation put a lower bound on inflaton-Higgs coupling by providing an effective mass for the Higgs, pushing back its hilltop during inflation. These additional constraints provide a lower bound on $\alpha$, while $\kappa$ must also be non-zero for the inflaton to decay efficiently.


I. INTRODUCTION
The confirmation of the Higgs particle at the LHC means that for the first time, we have a unitary theory of particle physics -the Standard Model (SM).Combined with the graviton, the theory appears to have internal consistency down to the Planck scale.However, the existence of the Higgs, the first discovered scalar in nature, opens up a new type of stability problem; as we now recap.The dimension 4 Higgs potential is given by where H is the Higgs doublet.This work will focus on the magnitude of H which will be denoted as the Higgs field, h, while its angular components are reorganized into the longitudinal modes of the W ± and Z bosons, and will not play a direct role for us here.As is well known, this potential has a non-zero local minimum at 246 GeV, which is classically stable and denotes our vacuum.At high energies, the Higgs interactions with various SM particles causes the value of the self coupling λ to evolve under renormalization group equations.For the central values of SM parameters, this causes λ to turn negative at large values of the Higgs field [1][2][3][4][5].This means the potential reaches a maximum at some large value of the Higgs field, h max , and then quickly drop to negative values.Figure 1 provides a plot of the effective potential from solving the two-loop renormalization group equations.This shows that the peak occurs at energies ∼ 10 10−11 GeV, depending sensitively on the top quark mass.Here we plot in black the central top mass, while the standard deviation values are given by the dashed curves [6,7].(Note that the initial dip associated with the usual electroweak vacuum is not directly relevant as the µ h h 2 /2 term in Eq. 1 is negligible at such high energies).The huge fluctuation required for the FIG. 1.
The two loop renormalized effective potential of the Higgs within the SM.The turnover at high scales around hmax (V ′ h (hmax) = 0) is focussed on.Blue-dashed curve is mtop = 172.5 GeV, black-solid curve is mtop = 172.9GeV, and purple-dashed curve is mtop = 173.3GeV.After the turnover, the potential runs negative.
Higgs field to go over, or tunnel through, this hilltop implies that it is not an issue in the late universe where the Higgs is centered at the electroweak scale 246 GeV.
However, in the very early universe, it is anticipated that the universe was at such extremely high energies that the Higgs is in danger of fluctuation to the dangerous side of the potential.Therefore, one should consider the serious possibility that the Higgs field could have fluctuated past h max during an early era; of particular focus here was an early era of inflation.If this had occurred, the Higgs could have run down the potential towards infinity, and the universe would presumably undergo a catastrophic crunch.Important earlier work in this area, includes Refs.[8][9][10][11][12][13][14][15][16][17][18].
One may potentially avoid this problem by significantly lowering the scale of inflation or utilizing a large effective mass of the Higgs during inflation, which we shall discuss later.But there is still a possible disaster that can take place after inflation during p/reheating.In particular, one anticipates a direct coupling between the Higgs and the inflaton, and this can cause (parametric or tachyonic) resonance during reheating that exponentially increases in the value of the Higgs.In fact a coupling between the inflaton and Higgs is required for efficient reheating of the SM as it is the only renormalizable way a (gauge singlet) inflaton can couple to the SM.Potential ways to handle this problem has been discussed in the literature, including Refs.[19][20][21][22][23][24][25][26].
Altogether, constraints on inflaton-Higgs couplings are needed to determine if a simple model in which the SM is taken seriously to high scales is possible; this is the focus of this work.We will examine all dimension 4 couplings between the Higgs and the inflaton, and derive combined bounds on this parameter space from demanding stability both during and after inflation -see ahead to Figure 12 for a final result.We use this to place a bound on the reheating temperature and discuss possible consequences.

II. INFLATON-HIGGS COUPLING
In conventional models of inflation, the exponential expansion of the universe is driven by a heavy scalar field ϕ -the inflaton.During inflation, the inflaton's potential was relatively constant with the field at very large values.As a result, the inflaton dominated the universe, and so the energy density of the early universe can be approximated as simply the sum of the inflaton's kinetic and potential energy As the universe expanded, most fields fell down their potential towards zero.However, since the inflaton potential was roughly constant at this time, its field, and in turn, ρ(t), would have also remained roughly constant.
The Friedman equation where H(t) is the Hubble parameter and a(t) is the scale factor, says that if ϕ undergoes sufficient Hubble friction and moves slowly, ρ(t) is nearly constant, so the solution for a(t) is approximately exponential a(t) ∝ e Ht .As the inflaton continues to roll, inflation will eventually end.Our primary interest here is in the era immediately following the end of inflation, (although we shall also consider further bounds from the inflationary phase itself).During this subsequent era, the inflaton must have decayed into the SM in order to reheat the universe.Plausibly, this was dominated by the inflaton's decay into the Higgs as it is the only possible renormalizable channel (for a gauge singlet inflaton).All other particles in the SM must couple to the inflaton via higher order operators and therefore may be highly suppressed.

A. Action and Inflaton Evolution
The action for gravity, the Higgs, and the inflaton is (units The −κ ϕ h 2 /2 term provides perturbative inflaton decay into two Higgs particles as illustrated in the upper part of Figure 2.Moreover, it is critical for understanding how the Higgs will grow during reheating, which can be resonant.The coefficient κ (which has units of mass) specifies its strength.The only other allowed renormalizable interaction is ∼ −α ϕ 2 h 2 /2 and will be discussed in Section VIII (see lower part of Figure 2).Other terms in the action, including terms for all other SM particles and their interactions, may be small during reheating because they can only couple to the (gauge singlet) inflaton through higher dimension operators, and will not be included in our simplified analysis.
With the above action, the following is the Heisenberg equation of motion for the inflaton To first approximation, the inflaton is a classical field with small quantum fluctuations Since the quantum fluctuations relative to the classical background are very small in standard inflationary models, the inflaton can be approximated as simply its classical, background field ϕ(t) when we study the impact on the Higgs.During reheating, the inflaton's potential can be approximated by a mass term around its minimum Using this potential in the (classical version of) Eq. 5 gives So the equation of motion is that of a damped harmonic oscillator with Hubble acting as a friction term.

B. Higgs Evolution Equations
The Heisenberg equation of motion for the Higgs is For simplicity, the partial derivative of the Higgs potential can be dropped from the equation of motion.It turns out that this term is around five orders of magnitude smaller than κϕ ĥ when we expand around the Higgs vacuum.It therefore has a negligible effect during the first stages of evolution.Of course, once the inflaton has completely decayed into the Higgs, the κ ϕ ĥ term goes away and the shape of the Higgs potential becomes critical.If the Higgs is on the other side of the hilltop at that point, then it would run down the potential and cause a disaster.Therefore, while the Higgs potential can be neglected during the first stages of reheating, it will be critical in our understanding of the final analysis.Therefore we shall incorporate its effects qualitatively by noting that this problem can occur.Ultimately, a more precise analysis would involve the inclusion of this term throughout the entire analysis.This renders the equations nonlinear, with a ∼ λ ĥ3 term.This could be handled with full lattice simulations.However, this is beyond the scope of the current work.Our results are therefore approximate and can be improved upon with more detailed simulations.However, the qualitative and semiquantitative results found here are anticipated to provide a reasonable estimate.
In this absence of the regular Higgs potential term, the resulting equation is linear.So it is useful to take a Fourier transform and the equation of motion is then given in terms of the wavenumber, k, as Since the equation of motion is now linear with respect to the Higgs, and the background metric and inflaton is approximated as carrying spatial translation invariance, the Higgs field can be solved in terms of mode functions v k (t), as follows Since the annihilation and creation operators are constant in time, they can be factored out, and thus the mode functions must satisfy the same equation of motion Thus need to solve this set of ordinary differential equations (ODEs) in order to solve the theory at the linear level.
In the free theory, where the Higgs and inflaton are not coupled and κ = 0, the solution to this equation of motion is where the particular solution is chosen to maintain the Higgs' commutation relations.In the coupled theory, this solution is an initial condition which the Higgs evolves away from as the inflaton evolves.We are interested in computing the evolution, as the coupling to the inflaton may cause a radical growth in its variance.This will imply the Higgs would have a significant probability of being on the dangerous side of the effective potential V h .The variance is readily calculated from Eq. 11 to be In the initial vacuum state and in absence of coupling, the variance in the Higgs field can be calculated by inserting Eq. 13 into Eq.14 or, after integrating out the angular components, While this integral is formally UV divergent, we will be only interested in a smoothed field with a cutoff set by a bubble size (to be discussed in Section V), so in fact the relevant fluctuations are finite (the two-point correlation function is finite on the scales of interest).Introducing non-zero values of κ corresponds to solutions that evolve away from Eq. 13.By dividing out the variance of the interacting theory by the variance of the free theory, the growth in variance can be determined.

III. SIMPLIFIED ANALYSIS
In the next Section, we shall solve these equations numerically.However, for now it is useful to gain intuition and develop approximate semi-analytical results.
The solution to Eq. 12 is different for each value of k.We will choose a discrete values for k and then sum the results.Eq. 12 can first be approximated by treating the time dependent prefactors as varying adiabatically slowly.This means it becomes a type of Mathieu equation.Solving this gives a better understanding of the best range of k values to sum over.The Mathieu equation has the canonical form Its solutions are of the form where f 1,2 (τ ) are periodic.If the so-called Floquet exponent µ is real the growth is exponential, otherwise the growth is absent.The Mathieu equation requires a harmonic driving term which can come from the inflaton field.As in Eq. 7, the inflaton potential resembles that of a harmonic oscillator during reheating, and so the field can be approximated as The inflaton's equation of motion 8, includes a Hubble friction term.This can be incorporated in this approximation by allowing the envelope ϕ amp to decrease with time as the universe expands.In particular, since reheating is approximately a matter dominated era (as the inflaton's oscillations lead to the pressure averaging to zero), the Hubble parameter during this time can be ap-proximated as so that This approximation is not precise right at the start of reheating, the end of inflation, but becomes more accurate over time as the universe transitions from the accelerating phase to an effective matter dominated phase.Under these approximations ϕ amp (t) decreases as To ensure that these are reasonable approximations, they are compared to numerical solutions for H(t) and ϕ(t) from solving Eqs. 2, 3, and 8.This is shown in Figure 3.As expected, these approximations become increasingly accurate at later times (the late time off-set in ϕ is just a phase shift).
To generate these and other plots, we take the inflaton mass to be m ϕ ≈ 1.4 × 10 13 GeV (a value expected in chaotic inflation [27]).Calculations are done in units of m ϕ , such that the reduced Planck mass is GeV ≈ 1.7 × 10 5 m ϕ , and Pl is the field value at the end of inflation under the quadratic potential approximation.
In addition to approximating the inflaton field as a classically oscillating function, the expansion terms that depend on time must be removed from Eq. 12 in order to approximate it as a Mathieu equation.The adiabtic approximation involves replacing in the inflaton field by decreasing ϕ amp with time according to Eq. 22 so as to capture its effects.Similarly, a(t) can be removed from the equation of motion with the following substitution Here, k phys (t) is the physical wavenumber at a particular point in time and decreases as the universe expands.By allowing k phys (t) and ϕ amp (t) to decrease with time, expansion will be approximately captured in the solution, and so the Hubble term can be dropped in this first analysis.Substituting these approximations into Eq.12 gives Comparing Eq. 24 to Eq. 17, the following change of variables must be made to reproduce a Mathieu equation The resulting Floquet exponents µ are plotted as contours in Figure 4 as a function of k phys (t) and ϕ amp (t) for κ = 4 × 10 −5 m ϕ .The red lines in Figure 4 represent a discrete set of values of comoving wavenumber k, which redshift to the left according to k phys (t) ∝ 1/a(t) and downwards according to ϕ amp (t) ∝ 1/a(t) 3/2 as the universe expands.From this approximation and further numerical analysis, the range 0.2 ≤ k/m ϕ ≤ 3 and resolution ∆k/m ϕ = 0.2, was found to be sufficient as all these lines pass through the yellow bands of significant resonance.

IV. FULL NUMERICAL ANALYSIS
With a range of k values chosen, v k (t) can be numerically calculated according to Eq. 12 and the Higgs variance can be calculated according to Eq. 14.The growth in the Higgs field is measured as the ratio between the variance in the coupled theory and the variance in the free theory For two values of κ, this growth is shown in the upper and lower panels of Figure 5.These plots illustrate how sensitive the Higgs field is to a small change in κ as a factor of two causes a difference in growth of 10 5 .To determine if the Higgs field will grow past h max , only the endpoint of the variance plots need to be considered.Figure 6 shows the growth of the Higgs field at late times as a function of κ.There appears to be a dip in Higgs growth near κ ≈ 10 −5.7 m ϕ where the Higgs field apparently decreases at late times.This is possible is because the inflaton is oscillating, and so resonance may just coincidentally occur when the inflaton is at a minimum for this value of κ.But more typically, larger κ tends to lead to much larger growth.In the following section, the amount of growth required for the Higgs to go over its hilltop will be discussed.

V. CREATING A BUBBLE
The probability distribution of any ground state in quantum mechanics with a quadratic potential is Gaussian.The Higgs field's probability distribution can therefore be approximated as a Gaussian, since its potential is primarily quadratic with a quartic term that only contributes a small correction.So, the probability that the field would grow past h max where the potential V h (h)

peaks is approximately
where we use the fact that we will be typically in the tail of the distribution, so an integral past the peak provides a small correction.
To compute this value, it is important to first note that an individual point in the Higgs field cannot grow independently of the points around it due to field gradients.Rather, a type of bubble will form with a characteristic size, R bub , and grow.The size of the bubble depends on its height.Numerically, this can be explained by energy conservation.The energy density of the Higgs field is given by where V h (h) is defined by Eq. 1.As the field fluctuates up the potential, the field around it will tend to grow as well so that the gradient and time derivatives balance this increase in the potential and the energy density does not change too quickly.Given that the field's gradient is approximately equal to its time derivative (for such relativistic bubbles), this balancing act requires that Noting that the slope of the bubble is approximately its height divided by its characteristic radius, the gradient can be approximated as Therefore, balancing the gradient with the potential means the height of a bubble is approximately proportional to the inverse of the bubble's size 1 2 which implies Furthermore, the variance ⟨h 2 ⟩ is also proportional to 1/R bub .This is because when calculating the variance for a bubble according to Eq. 14, the integral only needs to be computed up to wavelengths of about the size of the bubble.Shorter wavelengths would not affect the overall shape of the bubble and can therefore be neglected.For the vacuum, this gives a variance of So, it is clear that the size of the bubble will drop out from the final probability calculations.
A more precise analysis has determined that the ratio of field height to vacuum fluctuations is indeed of this order, but with a corrected O(1) prefactor [28] Generally, λ is a function of the Higgs field h, due to the renormalization group flow of the effective potential.It is convenient to have an approximate fitting function for it.A useful form is [8] For simpler calculations, λ can be approximated as a constant λ ≈ −0.008.This gives a potential that is somewhat accurate in the regime V h (h) < 0, as seen in Figure 7, which is the region of interest due to the turn-over in the Higgs potential.
In the vacuum, Eq. 33 can therefore be directly plugged into Eq.27 with and λ = −0.008 to find the probability of the Higgs going over its hilltop today, when the inflaton coupling has no influence.The result is that P ∼ e −3300 which recovers the well known result that the Higgs in the ordinary vacuum today is very stable.

VI. LIMITS ON CUBIC COUPLING
To find the probability of the Higgs reaching h max after resonance after inflation, the vacuum variance in Eq. 33 for the fitting function λ eff (h) versus using a simple constant λ = −0.008.The latter function requires some translations in field value in order to show this approximate matching, but it is only the curvature determined by λ that is of direct importance here.
(with h bub → h max ) needs to be rescaled by the increased variance after resonance, which we earlier computed.For a reasonable approximation, the probability can therefore be given as A very generous limit would be to demand that which would put the Higgs field at ∼ 40% chance of growing past h max and becoming unstable.Solving for the growth in variance for this limit gives Meanwhile, the most conservative limits would have no more than one Hubble patch go over the hilltop.This is a more reasonable requirement, because even if one patch goes over it could expand to cause the rest of the universe into a crunch.Since there were about ∼ e 3Ne Hubble patches, where N e ≈ 55, this would correspond to which leads to If one Hubble patch were to fall down the Higgs potential, however, the patches around it could pull it back to the safe side.Therefore, a more reasonable estimate may allow for a little more growth than Eq.39 Each of these limits is indicated in Figure 6, where the horizontal dotted lines represent where the Higgs field growth is 20 (cyan), 100 (orange), and 3300 (red).

VII. BOUND ON REHEAT TEMPERATURE
Limits on inflaton-Higgs couplings can also put constraints on the reheat temperature of the universe.First, the inflaton decay rate can be calculated under the reasonable assumption that a gauge-singlet inflaton primarily decays into the Higgs.This is because it is only to the Higgs that one can form renormalizable couplings.While other interactions are dimension 5 and above, which are plausibly suppressed by a very high scale, such as M Pl .At late times, after any resonance has occurred, perturbative tree-level decay into the Higgs from the cubic coupling can be most important where g h = 4 is the number of components in the Higgs.
Using the bound on κ in Eq. 42 and assuming an inflaton mass of ∼ 1.4 × 10 13 GeV, this corresponds to a bound on the decay rate of Γ(ϕ → hh) ≲ 2.2 × 10 26 sec −1 .
As the inflaton decays, the large amount of energy stored in its mass is converted into kinetic energy in the significantly lighter Higgs and then the rest of the SM, producing a bath of relativistic particles.These particles will rapidly thermalize at time t ∼ 1/Γ.Since Hubble at this time is around H ∼ 1/t as the universe becomes matter dominated during reheating, and so it is true that H ∼ Γ at this time.Given an approximate value for H, the energy density of the thermal bath can be calculated from the Friedman equation (eq.3).The temperature of this thermalized sea of radiation can then be calculated from its energy density where g = 106.75represents the degrees of freedom in the SM (if other new degrees of freedom are present, we assume it does not increase the total g too much).These calculations give a reheat temperature of Therefore, the calculated constraints on κ from Eq. 42 gives a reheat temperature bound of where the uncertainty arises from the uncertainty in inflaton mass.(This can be compared to the work in Ref. [14] that also placed a bound on the reheat temperature.)Interestingly, this temperature is well below the temperature required in GUT scale models of baryogenesis that may help to explain the matter to anti-matter asymmetry; so this is potentially quite important.

VIII. INCLUDING QUARTIC COUPLING
While the resonance due to the −κ ϕ h 2 /2 term in Eq. 4 may dominate the growth in the Higgs field, an inflaton annihilation term −α ϕ 2 h 2 /2 in the action can also be considered The Feynman diagram for this annihilation is given in the lower part of Figure 2.This changes the equation of motion for the mode functions to Following the same steps as before, the Mathieu approximations can be used again to analyze the effects of varying α by first setting κ = 0 and noticing that Thus, the change of variables for the pure α theory analogous to Eq. 25 are (50) The Floquet exponent for α = 10 −8 (with κ = 0) with respect to ϕ amp and k phys is shown in Figure 8.
Figure 9 shows the numerically calculated growth for a pure α theory as well as the previously discussed constraints.Unlike in the pure κ theory, the growth in the Higgs fields fluctuates even as α increases steadily.The reason for this may be due to the fact that the bands of resonance in Figure 8 are much more narrow than in Figure 4.As result, the universe may spend very little time in the resonance bands when κ = 0.Meanwhile, ϕ is consistently oscillating as a cosine, and so ϕ can accidentally pass through a maximum or minimum as it redshifts through the instability bands.As a result, a small shift in the bands can result in an altered amount of resonance either up or down.But the overall trend for large α is growth.
Limiting the Higgs growth to be a factor of around 100, as discussed earlier, gives the bound α < 10 −8 .

IX. COMBINING CUBIC AND QUARTIC
Finally, the growth in the Higgs field with both cubic κ and quartic α couplings included can be calculated numerically from Eq. 48.The resulting constraints are given in Figure 10.The plot recovers the limits as were previously discussed in this paper, if one sets one of the parameters towards zero.The contour lines again appear jagged, which is due to the same reason the pure α (and to some extent pure κ) results in Figure 9 were not smooth.Figure 10 is one of our primary results, however we shall constrain the allowed parameter space further in the next section.

X. INSTABILITY DURING INFLATION
Apart from resonance caused by inflaton couplings during reheating, the Higgs is also at risk of going over its hilltop during inflation due to the large energy scales at this time.This is because even though the inflaton is slowly rolling during inflation and there is no resonance into the Higgs, there can be important quantum fluctuations.Specifically, the Gibbons-Hawking de Sitter temperature of inflation is [29] As a result, any light scalar field fluctuates as ⟨h 2 ⟩ ≈ H inf /2π per Hubble time.Since inflation lasted about N e ≈ 55 Hubble times, the total fluctuation after random walks becomes where H end represents Hubble at the end of inflation and may be close to ∼ 10 13 GeV (it cannot be significantly larger, as this would overproduce gravitational waves).Such fluctuations already exceed the h max values seen in Figure 1.Interestingly, coupling to the inflaton can actually help this issue.During inflation, the inflaton field is very large and roughly constant.As a result, the coupling terms in the potential ∆V = 1 2 (κϕ inf +αϕ 2 inf )h 2 can be interpreted as an additional effective Higgs mass This effective mass can push h max to higher energies during inflation and help to avoid causing an instability during inflation.Combining this new effective contribution to the Higgs potential, with the renormalized potential, we obtained the total Higgs potential in Figure 11.
In order to save the Higgs from creating an instability, it is necessary that h (inf) max ≫ ⟨h 2 ⟩ ∼ 10 13 GeV.The same probability arguments from earlier can be applied to find the required h (inf) max .For a generous constraint, would put the probability the Higgs growing too large during inflation at around 40%.This corresponds to For a very conservative estimate, the probability of reaching h (inf) max can be set again to e −3Ne so that only one Hubble patch may become unstable FIG. 11.Higgs potential including additional effective masses of 0 (solid black), 1.6 × 10 12 GeV (dashed red), 10 13 GeV (dashed orange), and 2.8 × 10 13 GeV (dashed cyan).This is for the central top mass value of mtop = 172.9GeV.The red, orange, and cyan correspond to the generous, intermediate, and conservative limits on the effective mass, respectively.
This puts the hilltop at By adding the effective mass to the Higgs potential, the turnover h max can be determined as a function κ and α.The generous constraint corresponds to m eff > m crit = 1.6 × 10 12 GeV while the conservative constraint corresponds to the bound m eff > m crit = 2.8 × 10 13 GeV.An intermediate value can be taken to be m eff > m crit = 10 13 GeV.Note that towards the end of inflation ϕ decreases to the value ϕ inf , decreasing m eff and therefore this provides the tightest bound.By estimating the end of inflation as ϕ inf = √ 2 M Pl , the following constraint can be put on the coupling coefficients These constraints can be added to Figure 10 to produce Figure 12, which is our primary result.This figure extends beyond existing results and is useful to summarize the space of allowed couplings.
In this Figure, the white (blank) region is allowed, as the Higgs is sufficiently stable.While the colored regions are ruled out, due to either an instability during reheating (upper right) or an instability during inflation (lower left).Here we focus on κ > 0; for κ < 0, the constraints from inflation are more severe as this coupling worsens the problem (we are adopting the convention that ϕ inf > 0).
Figure 12 implies that α must be non-zero in order to avoid instabilities during inflation.At the same time, it must be true that κ is non-zero in order for inflaton decay to properly reheat the universe.While κ = 0 would not cause an instability during reheating or during inflation, it would mean that the inflaton could only produce Higgs particles via annihilation.This process requires two inflaton particles to find each other which becomes less and less probable with time.

XI. CONCLUSION
In this work, we have taken the SM seriously to high energies.Under this assumption, we have shown that (for weak couplings) the inflaton cubic coupling to the Higgs must have a coupling strength of at most κ < 1.6 × 10 −5 m ϕ ∼ 2.2 × 10 8 GeV, while the coefficient for the quartic coupling is constrained by α < 10 −8 .This arises from demanding that the probability of the Higgs field going over its hilltop remains low despite resonance effects.The upper bound on κ also places a bound on the the reheat temperature of T reh ≲ 9.2 × 10 9 GeV, since in this framework inflaton decays into Higgs would dominate.This may have important implications for classes of models of baryogenesis, which often appeal to extremely high temperatures.Further constraints arise during inflation from de Sitter fluctuations, implying α must be non-zero, while κ must also be non-zero for the inflaton to decay efficiently.Further work includes considering other possible couplings between the inflaton and SM particles.Although these would be higher dimension operators, it would be important to determine under what conditions this may impact the results found here.
Other directions include the work in Ref. [26], where it was found that if the inflaton-Higgs cubic coupling κ is significantly larger than that focussed on here (namely κ ∼ 0.5 m ϕ ), and if the inflaton is moderately light (m ϕ ≲ 10 12 GeV), one can cure the Higgs potential entirely (this would be in a region far to the right of that displayed in Fig. 12).It would be useful to explore a larger portion of parameter space to consider all these possibilities.

FIG. 3 .
FIG. 3. Upper: Comparison of exact and approximate solutions of the inflaton field.Lower: Comparison of exact and approximate solutions of the Hubble parameter.

FIG. 4 .
FIG.4.Floquet exponent for Higgs resonance due to cubic coupling decay for κ = 4 × 10 −5 m ϕ .Brighter yellow regions represent regions of greater resonance and a larger characteristic exponent as specified by the color legend.The red lines represent how a particular wavenumber k evolves as k phys (t) and ϕamp(t) decrease with the expansion of the universe.

FIG. 8 .
FIG.8.Floquet exponent for Higgs resonance due to quartic coupling for α = 10 −8 (and κ = 0).Brighter yellow regions represent regions of greater resonance and a larger characteristic exponent as specified by the color legend.The red lines represent how a particular wavenumber k evolves as k phys (t) and ϕamp(t) decrease with the expansion of the universe.

FIG. 12 .
FIG.12.Parameter space due to both reheating and inflation.The (dashed) cyan area corresponds to violating the most conservative constraints m eff > 2.8 × 10 13 GeV.The red area corresponding to violating the most generous constraint m eff > 1.6 × 10 12 GeV is not contained within the plot shown.The orange area corresponds to violating an intermediate value of m eff > 10 13 GeV.