lsection
[5.8em]
\contentslabel2.3em
\contentspage
Higgs inflation
Javier Rubio
Institut für Theoretische Physik, RuprechtKarlsUniversität Heidelberg,
Philosophenweg 16, 69120 Heidelberg, Germany
——————————————————————————————————————————
Abstract
The properties of the recently discovered Higgs boson together with the absence of new physics at collider experiments allows us to speculate about consistently extending the Standard Model of particle physics all the way up to the Planck scale. In this context, the Standard Model Higgs nonminimally coupled to gravity could be responsible for the symmetry properties of the Universe at large scales and for the generation of the primordial spectrum of curvature perturbations seeding structure formation. We overview the minimalistic Higgs inflation scenario, its predictions, open issues and extensions and discuss its interplay with the possible metastability of the Standard Model vacuum.
——————————————————————————————————————————
Email: j.rubio@thphys.uniheidelberg.de
Contents
1 Introduction and summary
Inflation is nowadays a wellestablished paradigm [Starobinsky:1980te, Guth:1980zm, Mukhanov:1981xt, Linde:1981mu, Albrecht:1982wi, Linde:1983gd] able to explain the flatness, homogeneity and isotropy of the Universe and the generation of the primordial density fluctuations seeding structure formation [Hawking:1982cz, Starobinsky:1982ee, Sasaki:1986hm, Mukhanov:1988jd]. In spite of the phenomenological success, the inflaton’s nature remains unknown and its role could be played by any particle physics candidate able to imitate a slowlymoving scalar field in the very early Universe.
In spite of dedicated searches, the only outcome of the Large Hadron Collider (LHC) experiments till date is a scalar particle with mass [Aad:2012tfa, Chatrchyan:2012xdj, Tanabashi:2018oca]
(1.1) 
and properties similar to those of the Standard Model (SM) Higgs. The mass value (1.1) is certainly particular since it allows to extend the SM up to the Planck scale without leaving the perturbative regime [Shaposhnikov:2009pv]. The main limitation to this appealing scenario is the potential instability of the electroweak vacuum at high energies. Roughly speaking, the value of the Higgs selfcoupling following from the SM renormalization group equations decreases with energy till a given scale and starts increasing thereafter. Whether it stays positive all the way up to the Planck scale, or turns negative at some intermediate scale depends, mainly, on the interplay between the Higgs mass and the top quark Yukawa coupling extracted from the reconstructed MonteCarlo top mass in collider experiments [Butenschoen:2016lpz], cf. Fig. 1. Neglecting the effect of gravitational corrections, the critical value separating the region of absolute stability from the metastability/instability^{1}^{1}1The metastability region is defined as the parameter space leading to vacuum instability at energies below the Planck scale but with an electroweak vacuum lifetime longer than the age of the Universe. regions is given by [Bezrukov:2014ina]
(1.2) 
with the strong coupling constant at the boson mass. Within the present experimental and theoretical uncertainties [Bezrukov:2014ina, Butenschoen:2016lpz, Espinosa:2016nld], the SM is compatible with vacuum instability, metastability and absolute stability [Bezrukov:2012sa, Degrassi:2012ry, Buttazzo:2013uya, Espinosa:2015qea, Espinosa:2015kwx], with the separation among the different cases strongly depending on the ultraviolet completion of gravity [Branchina:2013jra, Branchina:2014usa, Branchina:2014rva], cf. Fig. 2.
In the absence of physics beyond the SM, it is certainly tempting to identify the recently discovered Higgs boson with the inflaton condensate. Unfortunately, the Higgs selfinteraction significantly exceeds the value leading to a sufficiently long inflationary period without generating an excessively large spectrum of primordial density perturbations [Linde:1983gd]. The situation is unchanged if one considers the renormalization group enhanced potential following from the extrapolation of the SM renormalization group equations to the inflationary scale [Isidori:2007vm, Hamada:2013mya, Fairbairn:2014nxa]. The simplest way out is to modify the Higgs field kinetic term in the largefield regime. In Higgs inflation^{2}^{2}2An alternative possibility involving a derivative coupling among the Einstein tensor and the Higgs kinetic term was considered in Refs. [Germani:2010gm, Germani:2010ux, Fumagalli:2017cdo]. this is done by including a direct coupling to the Ricci scalar [Bezrukov:2007ep], namely^{3}^{3}3 Prior to Ref. [Bezrukov:2007ep], the effect of nonminimal couplings had been extensively studied in the literature (see for instance Refs. [Minkowski:1977aj, Zee:1978wi, Smolin:1979uz, Spokoiny:1984bd, Futamase:1987ua, Salopek:1988qh, Fakir:1990iu, Fakir:1990eg, Makino:1991sg, Fakir:1992cg, Kaiser:1994vs, CervantesCota:1994zf, CervantesCota:1995tz, Komatsu:1997hv]), but never with the SM Higgs as the inflaton.
(1.3) 
with the Higgs doublet and a dimensionless constant to be fixed by observations. The inclusion of the nonminimal coupling (1.3) can be understood as an inevitable consequence of the quantization of the SM in a gravitational background, where this term is indeed required for the renormalization of the energymomentum tensor [Callan:1970ze, Birrell:1982ix].
When written in the Einstein frame, the Higgs inflation scenario displays two distinct regimes. At low energies, it approximately coincides with the SM minimally coupled to gravity. At high energies, it becomes a chiral SM with no radial Higgs component [Bezrukov:2009db, Dutta:2007st]. In this latter regime, the effective Higgs potential becomes exponentially flat, allowing for inflation with the usual chaotic initial conditions. The associated predictions depend only on the number of efolds of inflation, which is itself related to the duration of the heating stage. As the type and strength of the interactions among the Higgs field and other SM particles is experimentally known, the duration of this entropy production era can be computed in detail [GarciaBellido:2008ab, Bezrukov:2008ut, Bezrukov:2011sz, Repond:2016sol, Ema:2016dny], leading to precise inflationary predictions in perfect agreement with observations [Akrami:2018odb].
The situation becomes more complicated when quantum corrections are included. The shape of the inflationary potential depends then on the values of the Higgs mass and top Yukawa coupling at the inflationary scale. In addition to the plateau already existing at treelevel [Bezrukov:2014bra, Fumagalli:2016lls, Fumagalli:2016sof, Enckell:2016xse], the renormalization group enhanced potential can develop a quasiinflection point along the inflationary trajectory [Allison:2013uaa, Bezrukov:2014bra, Hamada:2014iga, Bezrukov:2014ipa, Rubio:2015zia, Fumagalli:2016lls, Enckell:2016xse, Bezrukov:2017dyv, Rasanen:2017ivk] or a hilltop [Fumagalli:2016lls, Rasanen:2017ivk, Enckell:2018kkc], with different cases giving rise to different predictions.
Although a precise measurement of the inflationary observables could be understood as an interesting consistency check between cosmological observations and particle physics experiments [Espinosa:2007qp, Barvinsky:2008ia, DeSimone:2008ei, Bezrukov:2008ej, Bezrukov:2009db, Barvinsky:2009ii, Popa:2010xc, Bezrukov:2012sa, Salvio:2013rja], the low to highenergy connection is subject to unavoidable ambiguities related to the nonrenormalizability of the model [Barbon:2009ya, Burgess:2009ea, Burgess:2010zq, Bezrukov:2010jz, George:2013iia, Bezrukov:2014ipa, Rubio:2015zia, George:2015nza, Fumagalli:2016lls, Enckell:2016xse, Bezrukov:2017dyv]. In particular, the finite parts of the counterterms needed to renormalize the treelevel action lead to localized jumps in the SM renormalization group equations when connected to the chiral phase of Higgs inflation. The strength of these jumps encodes the remnants of the ultraviolet completion and cannot be determined within effective field theory approach [Bezrukov:2010jz, Bezrukov:2014ipa, Burgess:2014lza]. If the finite parts are significantly smaller than the associated coupling constants, Higgs inflation leads to a direct connection among the SM parameters measured at collider experiments and the large scale properties of the Universe, provided that the former do not give rise to vacuum instability. On the contrary, if the jumps in the coupling constants are large, the relation between high and lowenergy physics is lost, but Higgs inflation can surprisingly occur even if the electroweak vacuum is not completely stable [Bezrukov:2014ipa].
In this article we review the minimalistic Higgs inflation scenario, its predictions, open issues and extensions, and discuss its interplay with the potential metastability of the SM vacuum. The paper is organized as follows:

The general framework is introduced in Section 2. To illustrate the effect of nonminimal couplings, we consider an induced gravity scenario in which the effective Newton constant is completely determined by the Higgs vacuum expectation value. Having this toy model in mind, we reformulate Higgs inflation in the socalled Einstein frame in which the coupling to gravity is minimal and all nonlinearities appear in the scalar sector of the theory. After emphasizing the pole structure of the Einsteinframe kinetic term and its role in the asymptotic flatness of the Higgs inflation potential, we compute the treelevel inflationary observables and discuss the decoupling properties of the SM degrees of freedom.

The limitations of Higgs inflation as a fundamental theory are reviewed in Section 3. In particular, we present a detailed derivation of the cutoff scales signaling the violation of perturbative unitarity in different scattering processes and advocate the interpretation of Higgs inflation as an effective field theory to be supplemented by an infinite set of higher dimensional operators. Afterwards, we adopt a selfconsistent approach to Higgs inflation and formulate the set of assumptions leading to a controllable relation between low and highenergy observables. Based on the resulting framework, we analyze the contribution of quantum corrections to the renormalization group enhanced potential and their impact on the inflationary observables. We finish this section by discussing the potential issues of Higgs inflation with the metastability of the SM vacuum.

Some extensions and alternatives to the simplest Higgs inflation scenario are considered in Section 4. In particular, we address the difference between the metric and Palatini formulations of the theory and its extension to a fully scale invariant framework [Shaposhnikov:2008xb, GarciaBellido:2011de, Blas:2011ac, Bezrukov:2012hx, GarciaBellido:2012zu, Rubio:2014wta, Trashorras:2016azl, Karananas:2016kyt, Ferreira:2016vsc, Ferreira:2016wem, Ferreira:2016kxi, Casas:2017wjh, Ferreira:2018qss, Ferreira:2018itt, Casas:2018fum]. The inflationary predictions in these models are put in one to one correspondence with the pole structure of the Einsteinframe kinetic term, allowing for an easy comparison with the results of the standard Higgs inflation scenario.
Overall, we intend to complement the existing monographs in the literature [Bezrukov:2013fka, Bezrukov:2015dty, Moss:2015fma] by i) providing a further insight on the classical formulation of Higgs inflation and by ii) focusing on the uncertainties associated with the nonrenormalizability of the theory and their impact on model building.
2 General framework
The inflationary paradigm is usually formulated in terms of conditions on the local flatness on an arbitrary potential, which can in principle contain a large number of extrema and slopes [Artymowski:2016pjz]. This flatness is usually related to the existence of some approximate shiftsymmetry, which, for the purposes of Higgs inflation, is convenient to reformulate as a nonlinear realization of approximate scaleinvariance.
2.1 Induced gravity
Let us start by considering an induced gravity scenario
(2.1) 
involving a scalar field , a vector field and a fermion field , with interactions similar to those appearing in the SM of particle physics when written in the unitary gauge . The quantity stands for the standard kinetic term, which for simplicity we take to be Abelian. In this toy model, the effective Newton constant is induced by the scalar field expectation value,
(2.2) 
In order for to be wellbehaved, the nonminimal coupling is restricted to take positive values. This condition is equivalent to require the semipositive definiteness of the scalar field kinetic term, as can be easily seen by performing a field redefinition .
An important property of the induced gravity action (2.1) is its invariance under scale transformations
(2.3) 
Here is an arbitrary constant, compactly denotes the various fields in the model and ’s are their corresponding scaling dimensions. The consequences of this dilatation symmetry are more easily understood in a minimallycoupled frame displaying the standard EinsteinHilbert term. This Einstein frame is achieved by a Weyl redefinition of the metric^{4}^{4}4In spite of its extensive use in the literature, we avoid referring to pointwise rescalings of the metric as “conformal transformations.” For a comprehensive discussion on the differences between Weyl and conformal symmetries, see for instance Refs. [Karananas:2015ioa, Karananas:2015eha].
(2.4) 
together with a Weyl rescaling of the vector and fermion fields,
(2.5) 
After some trivial algebra we obtain an Einsteinframe action^{5}^{5}5In order to keep the notation as simple as possible, we will not introduce different notations for the quantities defined in different Weylrelated frames. In particular, the implicit Lorentz contractions in this article should be understood to be performed with the metric of the frame under consideration.
(2.6) 
containing a noncanonical term for the field. The coefficient of this kinetic term,
(2.7) 
involves a quadratic pole at and a constant
(2.8) 
varying between zero at and when . The field kinetic term can be made canonical by performing an additional field redefinition,
(2.9) 
mapping the vicinity of the pole at to . The resulting action
(2.10) 
is invariant under shift transformations , with a constant. The exponential mapping in Eq. (2.9) indicates that such translational symmetry is nothing else than the nonlinear realization of the original scale invariance we started with in Eq. (2.1) [Csaki:2014bua]. The Einsteinframe transition in Eq. (2.4) is indeed equivalent to the spontaneous breaking of dilatations, since we implicitly required the field to acquire a nonzero expectation value. The canonically normalized scalar field is the associated Goldstone boson and as such it is completely decoupled from the matter fields and . The nonminimal coupling to gravity effectively replaces by in all dimension4 interactions involving conformal degrees of freedom. Note, however, that this decoupling statement does not apply to scaleinvariant extensions including additional scalar fields [Kaiser:2010ps, Kaiser:2012ak, Kaiser:2013sna, GarciaBellido:2011de, Bezrukov:2012hx, Karananas:2016kyt] or other nonconformal interactions such as terms [Starobinsky:1980te, Gorbunov:2010bn, Gorbunov:2012ij, Gorbunov:2012ns].
2.2 Higgs inflation from approximate scale invariance
Although the toy model presented above contains many of the key ingredients of Higgs inflation, it is not phenomenologically viable. In particular, the Einsteinframe potential is completely shiftsymmetric and does not allow for inflation to end. On top of this limitation, the scalar field is completely decoupled from all conformal fields, excluding the possibility of entropy production and the eventual onset of a radiationdominated era. All these phenomenological limitations are intrinsically related to the exact realization of scale invariance and as such they should be expected to disappear once a (sizable) dimensionfull parameter is included into the action. This is precisely what happens in Higgs inflation. The total Higgs inflation action [Bezrukov:2007ep]
(2.11) 
contains two dimensionfull parameters: the reduced Planck GeV and the Higgs vacuum expectation value GeV responsible for the masses within the SM Lagrangian density . Among these two scales, the Planck mass is the most important one at the large field values relevant for inflation. To illustrate how the inclusion of modifies the results of the previous section, let us consider the graviscalar part of Eq. (2.11) in the unitary gauge , namely
(2.12) 
with
(2.13) 
the usual SM symmetrybreaking potential. As in the induced gravity scenario, the inclusion of the nonminimal coupling to gravity changes the strength of the gravitational interaction and makes it dependent on the Higgs field,
(2.14) 
In order for the graviton propagator to be welldefined at all values, the nonminimal coupling must be positive.^{6}^{6}6Models with negative have been considered in the literature [Herranen:2014cua, Kamada:2014ufa, Figueroa:2017slm]. In this type of scenarios the gravitational instability at large field values can be avoided by replacing the quadratic coupling by a designed function remaining smaller than during the whole field regime. If , this requirement translates into a weakening of the effective Newton constant at increasing Higgs values. For nonminimal couplings in the range , this effect is important in the largefield regime , but completely negligible otherwise.
As we did in Section 2.1, it is convenient to reformulate Eq. (2.12) in the Einstein frame by performing a Weyl transformation with
(2.15) 
In the new frame, all the nonlinearities associated with the nonminimal Higgsgravity interaction are moved to the scalar sector of the theory,
(2.16) 
which contains now a nonexactly flat potential
(2.17) 
and a noncanonical kinetic sector resulting from the rescaling of the metric determinant and the nonhomogeneous part of the Ricci scalar transformation. The kinetic function
(2.18) 
shares some similarities with that in Eq. (2.7). In particular, it contains two poles located respectively at and . The first pole is an inflationary pole, like the one appearing in the induced gravity scenario. This pole leads to an enhanced friction for the field around and allows for inflation to happen even if the potential is not sufficiently flat. The second pole is a Minkowski pole around which the Weyl transformation equals one and the usual SM action is approximately recovered. To see this explicitly, we carry out an additional field redefinition ,^{7}^{7}7Note that all equations till this point hold even if the nonminimal coupling is fielddependent [Ezquiaga:2017fvi, Masina:2018ejw].
(2.19) 
to recast Eq. (2.16) in terms of a canonically normalized scalar field . This differential equation admits an exact solution [GarciaBellido:2008ab]
(2.20) 
In terms of the original field , we can distinguish two asymptotic regimes
(2.21) 
separated by a critical value
(2.22) 
The comparison between these approximate expressions and the exact field redefinition in Eq. (2.20) is shown in Fig. 3.
The large hierarchy between the transition scale and the electroweak scale allows us to identify in practice the vacuum expectation value with . In this limit, the Einsteinframe potential (2.17) can be rewriten as
(2.23) 
with
(2.24) 
At we recover the usual Higgs potential (up to highly suppressed corrections, cf. Section 3.1). At the Einsteinframe potential becomes exponentially stretched and approaches the asymptotic value at . The presence of in Eq. (2.11) modifies also the decoupling properties of the Higgs field as compared to those in the induced gravity scenario. In particular, the masses of the intermediate gauge bosons and fermions in the Einsteinframe, ,^{8}^{8}8Here we use a compact notation for the gauge boson couplings, namely and for the and bosons respectively, with and the gauge couplings of the and SM groups and the weak mixing angle. The coupling denotes a generic Yukawa coupling.
(2.25) 
coincide with the SM masses in the small field regime () and evolve towards constant values proportional to in the largefield regime (). The transition to the Einsteinframe effectively replaces by in all (nonderivative) SM interactions. This behavior allows us to describe the Einsteinframe matter sector in terms of a chiral SM with vacuum expectation value [Bezrukov:2009db, Dutta:2007st].
2.3 Treelevel inflationary predictions
The flattening of the Einsteinframe potential (2.23) due to the pole allows for inflation with the usual slowroll conditions even if the potential is not sufficiently flat. Let us compute the inflationary observables in the corresponding region , where
(2.26) 
The statistical information of the primordial curvature fluctuations generated by a singlefield model like the one under consideration is mainly encoded in the twopoint correlation functions of scalar and tensor perturbations, or equivalently in their Fourier transform, the power spectra. Following the standard approach [Mukhanov:1990me], we parametrize these spectra in an almost scaleinvariant form,
(2.27) 
and compute the inflationary observables
(2.28) 
with
(2.29) 
the first and second slowroll parameters and the primes denoting derivatives with respect to . The quantities in (2.28) should be understood as evaluated at a field value , with
(2.30) 
the efold number at which the reference scale in Eq. (2.27) exits the horizon, i.e. . Here,
(2.31) 
stands for the field value at the end of inflation, which is defined, as usual, by the condition . Equation (2.30) admits an exact inversion,
(2.32) 
with the lower branch of the Lambert function and
(2.33) 
a rescaled number of efolds. Inserting Eq. (2.32) into (2.28) we get the following analytical expressions for the primordial scalar amplitude,
(2.34) 
its spectral tilt,
(2.35) 
and the tensortoscalar ratio
(2.36) 
At large , these predictions display an interesting attractor behavior, very similar to that appearing in attractor scenarios [Ferrara:2013rsa, Kallosh:2013yoa, Galante:2014ifa] (see also Ref. [Artymowski:2016pjz]). Indeed, by taking into account the lower bound on the Lambert function [DBLP:journals/corr/Chatzigeorgiou16],
(2.37) 
we can obtain the approximate expressions^{9}^{9}9Note that the expressions contain rather than .
(2.38) 
at . The free parameter (or equivalently the nonminimal coupling ) can be fixed by combining Eq. (2.34) with the normalization of the primordial spectrum at large scales [Akrami:2018odb],
(2.39) 
Doing this, we get a relation
(2.40) 
among the nonminimal coupling , the number of efolds and the Higgs selfcoupling .
The precise value of the number of efolds in Eqs. (2.38) and (2.40) depends on the whole postinflationary expansion and, in particular, on the duration of the heating stage. As the strength of the interactions among the Higgs field and the SM particles is experimentally known, the entropy production following the end of inflation can be computed in detail [GarciaBellido:2008ab, Bezrukov:2008ut, Repond:2016sol].^{10}^{10}10This allows, for instance, to distinguish Higgs inflation from Starobinsky inflation [Starobinsky:1980te, Bezrukov:2011gp]. The depletion of the Higgscondensate is dominated by the nonperturbative production of massive intermediate gauge bosons, which, contrary to the SM fermions, can experience bosonic amplification. Once created, the and bosons can decay into lighter SM fermions with a decay probability proportional to the instantaneous expectation value of the Higgs field . The onset of the radiationdomination era is determined either by i) the time at which the Higgs amplitude approaches the critical value where the effective potential becomes quartic or by ii) the moment at which the energy density into relativistic fermions approaches that of the Higgs condensate; whatever happens first. The estimates in Refs. [GarciaBellido:2008ab, Bezrukov:2008ut, Repond:2016sol] provide a range
(2.41) 
with the lower and upper bounds associated respectively with the cases i) and ii) above. For the upper limit of this narrow window, we have and we can rewrite Eq. (2.40) as a relation between and ,
(2.42) 
Note that a variation of the Higgs selfcoupling in this equation can be compensated by a change of the a priori unknown nonminimal coupling to gravity. For the treelevel value , the nonminimal coupling must be significantly larger than one, but still much smaller than the value leading to sizable modifications of the effective Newton constant at low energies. In this regime, the parameter is very close to its maximum value . This effective limit simplifies considerably the expression for the critical scale separating the low and highenergy regimes,
(2.43) 
and collapses the inflationary predictions to the attractor values [Bezrukov:2007ep]
(2.44) 
in very good agreement with the latest results of the Planck collaboration [Akrami:2018odb]. Note that, although computed in the Einstein frame, these predictions could have been alternatively obtained in the nonminimally coupled frame (2.12), provided a suitable redefinition of the slowroll parameters in order to account for the Weyl factor relating the two frames [Makino:1991sg, Fakir:1992cg, Komatsu:1999mt, Tsujikawa:2004my, Koh:2005qp, Chiba:2008ia, Weenink:2010rr, Flanagan:2004bz, Chiba:2013mha, Postma:2014vaa, Jarv:2014hma, Ren:2014sya, Jarv:2015kga, Kuusk:2015dda, Jarv:2016sow, Burns:2016ric, Karam:2017zno, Karamitsos:2017elm, Karamitsos:2018lur].
3 Effective field theory interpretation
The presence of gravity makes Higgs inflation perturbatively nonrenormalizable [Barbon:2009ya, Burgess:2009ea, Burgess:2010zq, Bezrukov:2010jz] and forbids its interpretation as an ultraviolet complete theory. The model should be therefore understood as an effective description valid up to a given cutoff scale [Bezrukov:2010jz, George:2015nza]. This cutoff could either indicate the onset of a strongly coupled regime to be studied within the model by nonperturbative techniques (such as resummations, lattice simulations or functional renormalization studies) [Aydemir:2012nz, Calmet:2013hia, Saltas:2015vsc, Escriva:2016cwl] or the appearance of new degrees of freedom beyond the initiallyassumed SM content [Giudice:2010ka, Barbon:2015fla].
3.1 The cutoff scale
A priori, the cutoff scale of Higgs inflation could coincide with the Planck scale, where gravitational effects should definitely taken into account. Although quite natural, the identification of these two energy scales may not be theoretically consistent, since other interactions could lead to violations of treelevel unitarity at a lower energy scale. An estimate^{11}^{11}11This procedure does not take into account possible cancellations among scattering diagrams, as those taking place, for instance, in models involving a singlet scalar field not minimally coupled to gravity [Hertzberg:2010dc]. of the cutoff scale can be obtained by expanding the fields around their background values, such that all kind of higher dimensional operators appear in the resulting action [Bezrukov:2010jz, Ferrara:2010in]. The computation is technically simpler in the original frame (2.11). In order to illustrate the procedure let us consider the graviscalar sector in Eq. (2.12). Expanding the fields around their background values and ,
(3.1) 
we obtain the following quadratic Lagrangian density for the perturbations and
(3.2)  
with denoting the trace of the metric excitations. For nonvanishing , the last term in this equation mixes the trace of the metric perturbation with the scalar perturbation [Barvinsky:2008ia, DeSimone:2008ei, Barvinsky:2009fy, Barvinsky:2009ii]. To identify the different cutoff scales one must first diagonalize the kinetic terms. This can be done by performing a redefinition of the perturbations with
(3.3)  
(3.4) 
Once Eq. (3.2) has been reduced to a diagonal form, we can proceed to read the cutoff scales. The easiest one to identify is that associated with purely gravitational interactions,
(3.5) 
which coincides with the effective Planck scale in Eq. (2.12). For scalargraviton interactions, the leadingorder higherdimensional operator is , where
(3.6) 
Although we have focused on the graviscalar sector of the theory, the lack of renormalizability associated with the nonminimal coupling to gravity permeates all SM sectors involving the Higgs field. One could study, for instance, the scattering of intermediate and bosons. Since we are working in the unitary gauge, it is sufficient to consider the longitudinal polarization. The modification of the Higgs kinetic term at large field values changes the delicate pattern of cancellations in the SM and leads to a treelevel unitarity violation at a scale
(3.7) 
Note that the above scales depend on the background field . For small field values (), the cutoffs (3.5), (3.6) and (3.7) coincide with those obtained by naively expanding the theory around the electroweak scale, namely , [Burgess:2009ea, Barbon:2009ya, Burgess:2010zq, Hertzberg:2010dc, Atkins:2010yg]. At large field values, (), the suppression scale depends on the particular process under consideration. For the graviscalar cutoff grows quadratically till , where it becomes linear in and traces the dynamical Planck mass in that regime, . On the other hand, the gauge cutoff smoothly interpolates between at and at . Note that all cutoffs scales become linear in at . This means that any operator constructed out of them, the Higgs field and some Wilson coefficients approaches a scaleinvariant form at large field values, namely
(3.8) 
3.2 Relation between high and lowenergy parameters
In what follows we will assume that the ultraviolet completion of the theory respects the original symmetries of the treelevel action, and in particular the approximate scale invariance of Eq. (2.12) in the largefield regime and the associated shiftsymmetry of its Einsteinframe formulation. This strong assumption forbids the generation of dangerous higherdimensional operators that would completely spoil the predictivity of the model. In some sense, this requirement is not very different from the one implicitly assumed in other inflationary models involving transPlanckian field displacements.
The minimal set of higherdimensional operators to be included on top of the treelevel action is the one generated by the theory itself via radiative corrections [Bezrukov:2010jz, Bezrukov:2014ipa]. The cancellation of the loop divergences stemming from the original action requires the inclusion of an infinite set of counterterms with a very specific structure. As in any other nonrenormalizable theory, the outcome of this subtraction procedure depends on the renormalization scheme, with different choices corresponding to different assumptions about the ultraviolet completion of the theory. Among the different subtractions setups, a dimensional regularization scheme involving a fielddependent subtraction point [Bezrukov:2009db]
(3.9) 
fits pretty well with the approximate scalesymmetry of Eq. (2.11) at largefield values.^{12}^{12}12The use of other schemes such as PauliVillars regularization or standard dimensional regularization with fieldindependent subtraction point leads to dilatationsymmetry breaking and the consequent bending of the Higgs inflation plateau due to radiative corrections, see for instance Refs. [Barvinsky:2008ia, DeSimone:2008ei, Barvinsky:2009fy, Barvinsky:2009ii]. Given this frame and scheme, the minimal set of higherdimensional operators generated by the theory can be computed in any Weylrelated frame provided that all fields and dimensionfull parameters are appropriately rescaled. The computation becomes particularly simple in the Einsteinframe, where the Weylrescaled renormalization point coincides with the standard fieldindependent prescription of renormalizable field theories, . A general counterterm in dimensional regularization contains a finite part and a divergent part in the form of a pole in , with the dimension of spacetime. The coefficient of the pole is chosen to cancel the loop divergences stemming from the original action. Once this divergent part is removed, we are left with the finite contribution . The strength of this term encodes the remnants of a particular ultraviolet completion and cannot be determined within the effective field theory approach [Bezrukov:2010jz, Bezrukov:2014ipa, Burgess:2014lza]. From a quantitative point of view, the most relevant terms are related to the Higgs and topquark interactions. In the Einsteinframe at one loop, they take the form [Bezrukov:2014ipa]
(3.10) 
where the primes denote again derivatives with respect to . Note that these operators differ, as expected, from those appearing in the treelevel action. This means that, while the contribution can be removed by a selfcoupling redefinition, the finite parts and should be promoted to new couplings constants. Once the associated operators are added to the treelevel action, the reevaluation of radiative corrections will generate additional contributions beyond the original oneloop result. These contributions come together with new finite parts that must be again promoted to novel couplings with their own renormalization group equations. The iteration of this scheme leads to a renormalized action including an infinite set of higherdimensional operators constructed out of the function and its derivatives. For small field values, the function becomes approximately linear (, ) and one recovers the SM nonminimally coupled to gravity up to highly suppressed interactions. In this limit, the coefficients of the infinite set of counterterms can be eliminated by a redefinition of the low energy couplings, as happens in a renormalizable theory. When evolving towards the inflationary region, the function becomes approximately constant (, ) and some of the previously absorbed finite parts are dynamically subtracted. The unknown finite parts modify therefore the running of the SM couplings at the transition region , such that the SM masses at the electroweak scale cannot be unambiguously related to their inflationary counterparts without a precise knowledge of the ultraviolet completion [Bezrukov:2014ipa, Burgess:2014lza, Hertzberg:2011rc].
If the finite contributions are of the same order as the loops generating them, the tower of higher dimensional operators generated by radiative corrections can be truncated [Bezrukov:2014ipa]. In this case, the effect of the 1loop threshold corrections can be imitated by an effective change^{13}^{13}13This replacement implicitly neglects the running of the finite parts and in the transition region . [Bezrukov:2014ipa]
(3.11) 
with and given by the SM renormalization group equations. We emphasize, however, that the truncation of the renormalization group equations is not essential for most of the results presented below, since, within the selfconsistent approach to Higgs inflation, the functional form of the effective action is almost insensitive to it [Bezrukov:2010jz, Bezrukov:2014ipa, Bezrukov:2017dyv].
3.3 Potential scenarios and inflationary predictions
To describe the impact of radiative corrections on the inflationary predictions, we will make use of the renormalization group enhanced potential. This is given by the one in Eq. (2.26) but with the Higgs selfcoupling replaced by its corresponding running value ,
(3.12) 
Note that we are not promoting the nonminimal coupling within to a running coupling —as done, for instance, in Ref. [Ezquiaga:2017fvi]—but rather assuming it to be constant during inflation. This is indeed a reasonable approximation since the oneloop beta function determining the running of [Bezrukov:2009db, Yoon:1996yr],
(3.13) 
is rather small for realistic values of the couplings constant at the inflationary scale, [Bezrukov:2017dyv, Masina:2018ejw] (see also Ref. [Salvio:2017oyf]).
Although, strictly speaking, the renormalization group enhanced potential is not gauge invariant, the gauge dependence is small during slowroll inflation, especially in the presence of extrema [Espinosa:2015qea, Espinosa:2016nld, Cook:2014dga]. In the vicinity of the minimum of , we can use the approximation [Bezrukov:2014bra]
(3.14) 
with the parameters , and depending on the inflationary values of the Einsteinframe Higgs and top quark masses, according to the fitting formulas [Bezrukov:2014bra]
(3.15)  
with and in GeV. As seen in the last expression, the parameter , standing for the derivative of the beta function for at the scale of inflation, is rather insensitive to the Higgs and top quark mass values at that scale and can be wellapproximated by . The choice
(3.16) 
with optimizes the convergence of perturbation theory [Bezrukov:2009db, Bezrukov:2008ej], while respecting the asymptotic symmetry of the treelevel action (2.12) and its nonlinear shiftsymmetric Einsteinframe realization. In the second equality, we have introduced an effective parameter to facilitate the numerical computation of the inflationary observables.
A simple inspection of Eqs. (3.12) and (3.14) allows us to distinguish three regimes:

Noncritical regime/Universal: If , the effective potential (3.12) is almost independent of the radiative logarithmic correction and can be well approximated by its treelevel form (2.26). Consequently, the inflationary observables retain their treelevel values [Fumagalli:2016lls, Enckell:2016xse, Bezrukov:2014bra, Bezrukov:2017dyv], cf. Fig. 4.

Critical regime: If , the first two derivatives of the potential are approximately zero, , leading to the appearance of a quasiinflection point at
(3.17) Qualitatively, the vast majority of inflationary efolds in this scenario takes place in the vicinity of the inflection point , while the inflationary observables depend on the form of the potential as some value .
Given the small value of the Higgs selfcoupling in this scenario, , the nonminimal coupling can be significantly smaller than in the universal regime, , while still satisfying the normalization condition (2.39) [Allison:2013uaa, Bezrukov:2014bra, Hamada:2014iga, Hamada:2014wna]. This drastic decrease of the nonminimal coupling alleviates the treelevel unitary problems discussed in Section 3.1 by raising the cutoff scale.
For small values, the tensortoscalar ratio can be rather large, [Allison:2013uaa, Bezrukov:2014bra, Hamada:2014iga, Hamada:2014wna] (see also Ref. [Masina:2018ejw]). Note, however, that although CMB data seems to be consistent with the primordial power spectrum at large scales, the simple expansion in (2.27) cannot accurately describe its global behavior since the running of the spectral tilt and its scale dependence also become considerably large, cf. Fig. 5.
The nonmonotonic evolution of the slowroll parameter in the vicinity of the inflection point leads to the enhancement of the spectrum of primordial density fluctuations at small and intermediate scales. It is important to notice at this point that the standard slowroll condition may break down if the potential becomes extremely flat and the inertial contribution in the equation of motion for the inflation field is not negligible as compared with the Hubble friction [Kannike:2017bxn, Germani:2017bcs, GarciaBellido:2017mdw]. In this regime, even the classical treatment is compromised since stochastic effects can no longer be ignored [Starobinsky:1994bd, Vennin:2015hra, Pattison:2017mbe, Ezquiaga:2018gbw].
Figure 5: (Left) Spectral tilt running in critical Higgs inflation as a function of the tensortoscalar ratio and the spectraltilt [Rasanen:2017ivk]. (Right) Scale dependence of the spectraltilt in the same case [Rasanen:2017ivk]. The purple dots indicate the universal/noncritical Higgs inflation regime. The boundaries on the righthand side of the figures correspond to the constraint on the number of efolds following the heating estimates in Refs. [GarciaBellido:2008ab, Bezrukov:2008ut, Repond:2016sol]. For lower heating efficiencies, the boundaries move to the left, decreasing the spectral tilt but not significantly affecting the tensortoscalar ratio [Rasanen:2017ivk]. If we restrict ourselves to situation in which the slowroll approximation is satisfied during the whole inflationary trajectory,[Bezrukov:2017dyv]^{14}^{14}14The onset of the slowroll regime prior to the arrival of the field to the inflection point and its dependence on preinflationary conditions was studied in Ref. [Salvio:2017oyf], where a robust inflationary attractor was shown to exist. the height and width of the generated bump at fixed spectral tilt are correlated with the tensortoscalar ratio , cf. Fig. 4. Contrary to some claims in the literature [Ezquiaga:2017fvi], the maximum amplitude of the powerspectrum compatible with the 95% C.L Planck contours [Bezrukov:2017dyv] is well below the critical threshold needed for primordial black hole formation [Carr:2016drx, Bird:2016dcv, Carr:2017jsz] (see, however, Refs. [Ezquiaga:2018gbw] and [Rasanen:2018fom]). This conclusion is unchanged if one considers the effect of noninstantaneous threshold corrections [Bezrukov:2017dyv], which could potentially affect the results given the numerical proximity of the inflection point to the upper boundary of the transition region, .

Hilltop regime: If the potential develops a new minimum at large field values [Fumagalli:2016sof, Rasanen:2017ivk]. This minimum is separated from the electroweak minimum by a local maximum where hilltop inflation can take place [Boubekeur:2005zm, Barenboim:2016mmw]. This scenario is highly sensitive to the initial conditions since the inflaton field must start on the electroweak vacuum side and close enough to the local maximum in order to support an extended inflationary epoch. On top of that, the fitting formulas in (3.3) may not be accurate enough for this case, since they are based on an optimization procedure around the minimum. The tensortoscalar ratio in this scenario differs also from the universal/noncritical Higgs inflation regime, but contrary to the critical case, it is decreased to , rather than increased [Rasanen:2017ivk, Fumagalli:2016lls].
3.4 Vacuum metastability and hightemperature effects
The qualitative classification of scenarios and predictions presented in the previous section depends on the inflationary values of the Higgs and top quark masses and holds independently of the value of their electroweak counterparts. In particular, any pair of couplings following from the SM renormalization group equations can be connected to a wellbehaved pair of couplings in the chiral phase by a proper choice of the unknown threshold corrections. This applies also if the SM vacuum is not completely stable. Some examples of the 1loop threshold correction needed to restore the universal/noncritical Higgs inflation scenario beyond and GeV are shown in Fig. 6. For a detailed scan of the parameter space see Refs. [Fumagalli:2016lls, Enckell:2016xse].
The nontrivial interplay between vacuum stability and threshold corrections generates an additional minimum at large field values. Provided the usual chaotic initial conditions, the Higgs field will evolve in the transPlanckian field regime, inflating the Universe while moving towards smaller field values. Since the new minimum is significantly wider and deeper than the electroweak one, it seems likely that the Higgs field will finish its postinflationary evolution there. Note, however, that this conclusion is strongly dependent on the ratio of the Higgs energy density to the second minimum depth. If this ratio is large, the entropy production at the end of inflation may significantly modify the shape of the potential, triggering its stabilization and allowing the Higgs field to evolve towards the desired electroweak vacuum [Bezrukov:2008ut].
The oneloop finite temperature corrections to be added on top of the Einsteinframe renormalization group enhanced potential take the form [Linde:1978px]
(3.18) 
with the plus and minus signs corresponding respectively to fermions and bosons and standing for the Einsteinframe masses in Eq. (2.25). The most important contributions in Eq. (3.18) are associated with the top quark and the electroweak bosons, with the corresponding coupling constants and evaluated at and , in order to minimize the radiative corrections [Kajantie:1995dw].
A detailed analysis of the universal/non critical Higgs inflation scenario reveals that the temperature of the decay products generated during the heating stage exceeds generically the temperature at which the unwanted secondary vacuum at large field values disappears [GarciaBellido:2008ab, Bezrukov:2008ut], see Fig 6.^{15}^{15}15A detailed scan of the parameter space assuming instantaneous conversion of the inflaton energy density into a thermal bath was performed in Ref. [Enckell:2016xse]. The stabilization becomes favored for increasing values^{16}^{16}16The larger is, the shallower and narrower the “wrong” minimum becomes, cf. Fig. 6. and holds even if this scale is as low as GeV [Bezrukov:2008ut]. The thermallycorrected potential enables the Higgs field to relax to the SM vacuum. After the heating stage, the temperature decreases as the Universe expands and the secondary minimum reappears, first as a local minimum and eventually as the global one. When that happens, the Higgs field is already trapped in the electroweak vacuum. Although the barrier separating the two minima prevents a direct decay, the Higgs field could still tunnel to the global minimum. The probability for this to happen is, however, very small and the lifetime of SM vacuum significantly exceeds the life of the Universe [Anderson:1990aa, Arnold:1991cv, Espinosa:1995se, Espinosa:2007qp]. Universal/noncritical Higgs inflation with a graceful exit can therefore take place for electroweak SM pole masses leading to vacuum metastability at energies below the inflationary scale. [Bezrukov:2008ut].
The situation changes completely if one considers the critical Higgs inflation scenario. In this case, the energy of the Higgs condensate is comparable to the depth of the secondary minimum and symmetry restoration does not take place. Unless the initial conditions are extremely finetuned, the Higgs field will relax to the minimum of the potential at Planckian values, leading with it to the inevitably collapse of the Universe [Felder:2002jk]. The success of critical Higgs inflation requires therefore the absolute stability of the electroweak vacuum [Bezrukov:2008ut].
4 Variations and extensions
Many variations and extensions of Higgs inflation have been considered in the literature, see for instance Refs. [Lerner:2010mq, Giudice:2010ka, BenDayan:2010yz, Kamada:2010qe, Lerner:2011ge, Arai:2011nq, Steinwachs:2013tr, Kamada:2012se, Einhorn:2012ih, Greenwood:2012aj, Kanemura:2012ha, Choudhury:2013zna, Oda:2014rpa, Xianyu:2014eba, Oda:2014mza, Ellis:2014dxa, He:2014ora, Kamada:2015sca, Okada:2015zfa, Cai:2015gla, Lazarides:2015xia, vandeBruck:2015gjd, Takahashi:2016uyq, Calmet:2016fsr, Ellis:2016spb, Okada:2016ssd, Ge:2016xcq, Marian:2017bvz, Ema:2017ckf, Ema:2017rqn, He:2018gyf, Chen:2018ucf]. In what follows we will restrict ourselves to those proposals that are more closely related to the minimalistic spirit of the original scenario. In particular, we will address a Palatini formulation of Higgs inflation and the embedding of the model to a fully scale invariant framework.
4.1 Palatini Higgs inflation
In the usual formulation of Higgs inflation presented in Section 2.2, the action is minimized with respect to the metric. This procedure implicitly assumes the existence of a LeviCivita connection depending on the metric tensor and the inclusion of a YorkHawkingGibbons term ensuring the cancellation of a total derivative term with novanishing variation at the boundary [York:1972sj, Gibbons:1976ue]. One could alternatively consider a Palatini formulation of gravity in which the metric tensor and the connection are treated as independent variables and no additional boundary term is required to obtain the equations of motion [Ferraris1982]. Roughly speaking, this formulation corresponds to assuming an ultraviolet completion involving different gravitational degrees of freedom.
Although the metric and Palatini formulations of General Relativity give rise to the same equations of motion [Ferraris1982], this is not true for scalartensor theories as Higgs inflation. To see this explicitly let us consider the Higgs inflation action in Eq. (2.12) with and a nonLeviCivita connection.^{17}^{17}17For a recent generalization built from the Higgs, the metric and the connection and involving only up to two derivatives see Ref. [Rasanen:2018ihz]. Performing a Weyl rescaling of the metric with given by Eq. (2.15) we obtain an Einsteinframe action
(4.1) 
containing a potential (2.17) and a noncanonical kinetic term with
(4.2) 
and
(4.3) 
Note that the kinetic function (4.2) differs from that obtained in the metric formulation, see Eq. (2.18). In particular, it does not contain the part associated with the nonhomogeneous transformation of the Ricci scalar, since is now invariant under Weyl rescalings. For the purposes of inflation, this translates into a modification of the residue of the inflationary pole at with respect to the metric case. While the metric value of in Eq. (2.18) is bounded from above (cf. Eq. (2.8)), it can take positive arbitrary values in the Palatini formulation (cf. Eq. (4.3)). Performing a field redefinition
(4.4) 
to canonically normalize the field kinetic term, we can rewrite the graviscalar action (4.1) at as
(4.5) 
with
(4.6) 
The comparison of the latest expression with Eq. (2.24) reveals some important differences between the metric and Palatini formulations. In both cases, the effective Einsteinframe potential smoothly interpolates between a lowenergy quartic potential and an asymptotically flat potential at large field values. Note, however, that the transition in the Palatini case is rather direct and does not involve the quadratic piece appearing in the metric formulation. On top of that, the flatness of the asymptotic plateau is different in the two cases, due to the effective change in . The Palatini dependence has a strong impact on the inflationary observables. In the largefield regime they read
(4.7) 
with
(4.8) 
a rescaled number of efolds and
(4.9) 
the inflaton value at the end of inflation (), with corresponding to the limit and to the end of inflation in a minimally coupled theory. A relation between the nonminimal coupling , the selfcoupling and the number of efolds can be obtained by taking into account the amplitude of the observed power spectrum in Eq. (2.39),
(4.10) 
A simple inspection of Eq. (4.7) reveals that the predicted tensortoscalar ratio in Palatini Higgs inflation is within the reach of current or future experiments [Matsumura:2016sri] only if , which, assuming , requires a very small coupling . For a discussion of unitarity violations in the Palatini formulation see Ref. [Bauer:2010jg].
4.2 HiggsDilaton model
The existence of robust predictions in (noncritical) Higgs inflation is intimately related to the emerging dilatation symmetry of its treelevel action at large field values. The uplifting of Higgs inflation to a completely scaleinvariant setting was considered in Refs. [Shaposhnikov:2008xb, GarciaBellido:2011de, Blas:2011ac, Bezrukov:2012hx, GarciaBellido:2012zu, Rubio:2014wta, Trashorras:2016azl, Karananas:2016kyt, Ferreira:2016vsc, Ferreira:2016wem, Ferreira:2016kxi, Ferreira:2018itt, Casas:2017wjh, Ferreira:2018qss, Casas:2018fum]. In the unitary gauge , the graviscalar sector of the HiggsDilaton model considered in these papers takes the form
(4.11) 
with
(4.12) 
a scaleinvariant version of the SM symmetrybreaking potential and positive dimensionless parameters. The existence of a welldefined gravitational interactions at all field values requires the nonminimal gravitational couplings to be positivedefinite, i.e. . In the absence of gravity, the ground state of Eq. (4.11) is determined by the scaleinvariant potential (4.12). For and , the vacuum manifold extends along the flat directions . Any solution with breaks scale symmetry spontaneously and induces nonzero values for the effective Planck mass and the electroweak scale.^{18}^{18}18 Among the possible values of in the presence of gravity, the case seems also preferred [GarciaBellido:2011de, Allen:1987tz, Jalmuzna:2011qw], see also Refs. [Antoniadis:1985pj, Tsamis:1992sx, Tsamis:1994ca, Antoniadis:2006wq, Polyakov:2009nq, Wetterich:2017ixo]. The relation between these highly hierarchical scales is set by finetuning . For this small value, the flat valleys in the potential are essentially aligned and we can safely approximate for all inflationary purposes.
To compare the inflationary predictions of this model with those of the standard Higgsinflation scenario, let us perform a Weyl rescaling followed by a field redefinition [Casas:2017wjh]