Quasi-Conformal Models and the Early Universe

Extensions of the Standard Model and general relativity featuring a UV fixed point can leave observable implications at accessible energies. Although mass parameters such as the Planck scale can appear through dimensional transmutation, all fundamental dimension-4 operators can (at least approximately) respect Weyl invariance at finite energy. An example is the Weyl-squared term, whose consistency and observational consequences are studied. This quasi-conformal scenario emerges from the UV complete quadratic gravity and is a possible framework for inflation. We find two realizations. In the first one the inflaton is a fundamental scalar with a quasi-conformal non-minimal coupling to the Ricci scalar. In this case the field excursion must not exceed the Planck mass by far. An example discussed in detail is hilltop inflation. In the second realization the inflaton is a pseudo-Goldstone boson (natural inflation). In this case we show how to obtain an elegant UV completion within an asymptotically free QCD-like theory, in which the inflaton is a composite scalar due to new strong dynamics. We also show how efficient reheating can occur. Unlike the natural inflation based on Einstein gravity, the tensor-to-scalar ratio is well below the current bound set by Planck. In both realizations mentioned above, the basic inflationary formulae are computed analytically and, therefore, these possibilities can be used as simple benchmark models.


Introduction
In 2015 (and with a recent update in 2018) the Planck collaboration [1,2] was able to exclude (or set stringent bounds on) several inflationary models. Indeed, better determinations of many observables related to the cosmic microwave background (CMB) (such as the tensor-to-scalar ratio r, the scalar spectral index n s and the curvature power spectrum P R ) were provided. Further improvements are expected in the next future: CMB Stage 4 (S4) will be active soon 1 . Therefore, early universe cosmology continues to be an exciting research area.
Despite this big progress, several models of the early universe are still allowed and we have only a partial understanding of how the universe initially expanded and evolved. Such large ambiguity can be reduced by looking for theoretical reasons to discriminate among the various possibilities. One way to do so is to focus on UV-complete theories that are applicable to the early universe.
In Ref. [3] it was shown how a relativistic field theory of all interactions (gravity included) can reach infinite energy. In order to achieve this goal, terms quadratic in the curvature are added to the Einstein-Hilbert action, so this theory is also known as quadratic gravity (QG) (see [4] for a review). The UV completion is obtained by demanding that the Weyl symmetry breaking terms vanish in the infinite energy limit and all couplings enjoy a UV fixed point. Therefore, this is a way of implementing asymptotic safety. An important difference between this possibility and the original proposal made by Weinberg in [5], though, is the presence of only a finite number of terms in the fundamental action, which guarantees the theory to be predictive and calculable. Another interesting property of QG is the possibility to solve the hierarchy problem [3,6], in the sense that the Higgs mass can be made technically natural [7]. When one lowers the energy from the UV-conformal fixed point mentioned above, the dimensionless Weyl-symmetry breaking operators are generated, but since they are sourced only by multiloop diagrams [8][9][10], they can remain small down to the inflationary energies. This corresponds to a scenario for the early universe, which we call "quasi-conformal" 2 .
One of the quadratic terms in the action is the squared of the Weyl tensor, which is necessary to keep the theory UV complete. Featuring higher time derivatives, this Weyl-squared term brought a number of issues, which, however, have been addressed in the recent years. One of the purposes of this paper is to review and extend these arguments in favour of the viability of such term.
Another aspect that has not been investigated so far is the observational predictions of the quasi-conformal scenario. The natural arena to study these effects is the early universe, which typically involves extremely high energies and field scales, sometimes reaching the Planck mass M P . This is the main topic of the present paper 3 .
The "quasi" of "quasi-conformal" is important. It reminds us that the Einstein-Hilbert term and other Weyl-breaking parameters can be present in the Lagrangian, generated, for example, through dimensional transmutation and various quantum effects [3,6,[15][16][17][18][19][20][21]. Note that an exactly Weyl invariant model would not be physically different from a Weyl-breaking one; indeed, any exactly conformal model can be written as a non-conformal one by fixing a gauge for Weyl symmetry and, conversely, any model of gravity coupled to an arbitrary matter sector can always be made exactly Weyl invariant 4 . What is proposed here is quasi-conformality, in which Weyl symmetry is approximate at high energy, but not exact apart from the strict UV limit. This is also consistent with the presence of the Weyl anomaly, which would anyhow render exact conformality impossible (when the theory is not at the fixed point, for instance at accessible energies).
A first question one can ask is whether the Standard Model (SM) Higgs can play the role of the inflaton in this scenario. The answer is negative because the requirement of a quasi-conformal non-minimal coupling (to the Ricci scalar) leads essentially to the same problem one has for a vanishing non-minimal coupling [22]: the SM Higgs cannot be responsible for inflation and the generation of the observed perturbations at the same time. We have verified the validity of this statement for a quasi-conformal non-minimal coupling by using the precise 2-loop effective potential as in [23][24][25][26][27]. A way to circumvent this problem is to demand a large [28,29] (or order 10 [30][31][32][33]) non-minimal coupling, which, however, would completely break Weyl symmetry in the UV.
Another apparently natural candidate for the inflaton is the effective spinless field corresponding to the square of the Ricci scalar, as proposed by Starobinsky [34]. However, in the quasiconformal scenario the coefficient of this term has to be very small because is not Weyl-invariant, while the observed value of the curvature power spectrum requires a very large coefficient.
In this paper we discuss two implementations of the quasi-conformal scenario. The first one is the class of models where the inflaton features a moderate field excursion: with "moderate" we mean that it does not exceed the Planck scale by far. This is because a quasi-conformal nonminimal coupling prohibits, as we will discuss in Sec. 3, field excursions bigger than √ 6M P . A specific realization that we study in detail here is hilltop inflation [35].
Another interesting implementation that we investigate in this paper is inflation triggered by a pseudo-Goldstone boson, what is known as "natural inflation" in the context of the non-conformal Einstein gravity [36]. Indeed, in this case the inflaton does not need to a have a (quasi-)Weyl invariant non-minimal coupling being morally the phase of some field representation of a spontaneously broken global group. An attractive feature of natural inflation is the natural flatness of the potential, which is protected by Goldstone theorem from large quantum effects. This is preserved in (quasi-)conformal realizations and we will study possible ways of distinguishing this setup from the one based on Einstein gravity.
The paper is organized as follows. In the next section we describe in full generality the quasiconformal scenario. As anticipated before, the conformal Weyl-squared term is introduced and we review and extend arguments which indicate the viability of this term. In Sec. 3 we study the hilltop quasi-conformal model, while in Sec. 4 we focus on the natural-inflation option, investigating how the universe inflated and subsequently evolved in this case. Finally, Sec. 5 offers our conclusions and some outlook of this work.

The quasi-conformal scenario and the Weyl-squared term
The quasi-conformal scenario studied in this work is a subset of the theories with scale invariance in the UV. Therefore, we start here by reviewing this class of theories. In the most general case, the field content features gauge fields with field strength F A µν , Weyl fermions ψ j and real scalars φ a . These fields should be considered as fundamental. Other composite states can appear below some confinement scales as usual.
The full action in the so-called Jordan frame is where g is the determinant of the metric g µν . We describe in turn the two pieces L SI and L SI . L SI is the scale invariant part. Since we are interested in theories with an approximate Weyl symmetry in the UV we split the first term as L SI = L Weyl + L Weyl , where the first part is invariant under Weyl transformations and the second part is not invariant. If one rewrites the coefficients ζ ab as ζ ab = ξ ab + δ ab /6 one recovers the usual notation for the non-minimal couplings ξ ab . Here W 2 ≡ W µνρσ W µνρσ , where W µνρσ is the Weyl tensor. L SI contains dimensionless parameters only: the gravitational ones f 0 and f 2 , the nonminimal couplings ξ ab and the Yukawa Y a ij , quartic couplings λ abcd and the gauge couplings in the covariant derivatives D µ .

L
SI contains instead the massive parameters and, although irrelevant in the UV, they play an important role in the IR: whereM P is the reduced Planck scale, Λ is the cosmological constant and we have introduced generic dimensionful parameters, m 2 ab , A abc and M ij , in the scalar and fermion sectors. This scenario has the remarkable property of renormalizability [37,38]. Moreover, for small enough values of f 2 , namely the hierarchy problem is solved [3,6]: gravity does not lead to an excessive radiative contribution to the Higgs mass, which can, therefore, obtain naturally [7] its observed value. Indeed, in this case the Higgs field acquires a shift symmetry, which protects its mass from large quantum effects. However, this theory was thought to suffer from two problems:

I
The clash between asymptotic freedom (or, more generally, asymptotic safety) and stability (understood as the absence of tachyons) [3,4,39,40]: whenever the parameters are chosen to ensure stability, perturbation theory features a Landau pole in the RG flow of f 0 .

II
The presence of a field, with spin-2 and mass but with an unusual sign of the kinetic term [37] (i.e. a ghost-like field). This is a generic issue of theories with more than two time-derivatives in the Lagrangian, which have a classical Hamiltonian that is unbounded from below, as shown by Ostrogradsky [41]. The theory above has indeed more than two derivatives because of the Weyl-squared term, W 2 .
However, in the following we will review (and extend in the case of Problem II), recent works that have shown how to address these issues. The solution of problem I is particularly relevant for this work because it leads to the quasi-conformal scenario, as we now describe. In Ref. [3] it was pointed out that Problem I is an artefact of perturbation theory in f 0 : when this parameter grows, perturbativity in f 0 is lost, but one can develop an expansion in 1/f 0 showing that no Landau poles are present, provided that all fundamental scalars have asymptotically Weyl-invariant couplings (1/f 0 → 0, ζ ab → 0) and the remaining couplings approach fixed points in the UV limit. In other words, the theory flows to a conformal behavior in the infinite energy limit and, in this way, can be UV complete. Although flowing to conformal gravity at infinite energy is consistent, at finite energy conformal invariance is broken by the scale anomaly and the R 2 term and non-vanishing values of ζ ab are generated. However, as discussed in [3], this is a multiloop effect (see [8][9][10] and references therein). Therefore, 1/f 0 and ζ ab can remain tiny down to the inflationary scales and the subsequent stages of the evolution of the universe. This leads to the quasi-conformal scenario, which is investigated in this work 5 .
Let us now turn to Problem II. Although higher derivatives can potentially be dangerous for a theory, it can be shown that this does not need to be so for the theory we consider here. In order to present a quantitative argument it is necessary to proceed perturbatively: we split g µν as where g cl µν is a classical background that solves the classical equations of motion (EOMs) andĥ µν is a quantum fluctuation. Although a non-perturbative approach would be desirable, that is not currently within the reach of our understanding of quantum field theory 6 . At this point it is clear that the W 2 term poses two problems: one that concerns the classical evolution (related to g cl µν ) and another one that regards the quantum dynamics (encoded inĥ µν ).
The fundamental question at the classical level is the following: can we avoid the possible instabilities due to the Ostrogradsky theorem? One obvious way to avoid runaways is to take f 2 1, such that the Weyl-squared term would have a negligible effect on any experimentally observable quantities. However, one would like to keep f 2 small both to see the effects of this term and to solve the hierarchy problem (see (2.5)). Moreover, when the massive spin-2 field is not tachyonic, M 2 2 > 0, it is also asymptotically free: f 2 goes to zero in the infinite energy limit [40,43,44]. In Ref. [45] it was shown that the answer to the question above is positive for the theory considered here even for a small f 2 (see also Ref. [46] for a related calculation). The key properties used in [45] was the fact that the massive spin-2 field decouples as f 2 → 0 and is not tachyonic. We here highlight the main steps of the argument in [45].
• First, one observes that in the free-field limit the Hamiltonian of the massive spin-2 field is where Q α and P α are the associated canonical variables and conjugate momenta and the helicity sum is over α = ±2, ±1, 0 because this massive particle has spin 2. The manifestation of the Ostrogradsky theorem in this case is the overall minus sign. However, despite that sign there are no instabilities in the free-field limit (that sign cancels in the EOM).
• The effective field theory (EFT) approach tells us that at energies below M 2 we should not find runaways even if the massive spin-2 field has an order-one coupling, f 2 ∼ 1.
• The intermediate case 0 < f 2 < 1 must have intermediate energy thresholds (above which the runaways might be activated).
• The weak coupling case f 2 1 (compatible with Higgs naturalness) must have an energy thresholds much larger than M 2 : we could see the effect of the massive spin-2 field without runaways.
The complete argument of Ref. [45] leads to the following energy thresholds. When the energies E associated with the derivatives of the spin-2 fields (both the massless graviton and the massive one) satisfy (for the derivatives of the spin-2 fields) (2.7) and the energies E associated with the matter sector (due to derivatives of matter fields or matterfield values times coupling constants) fulfill the runaways are avoided. Even for f 2 as small as 10 −8 this threshold is high enough to accommodate the whole cosmology 7 [45]. We will keep in mind these bounds when analysing the quasi-conformal models of the early universe in Secs. 3 and 4. Although the argument above gives a satisfactory solution of the classical Ostrogradsky problem it is still needed to address the issues raised by the quantization. One reason is the fact that quantum effects might lead to tunneling above the energy threshold even if the classical fields satisfy the bounds in (2.7) and (2.8). However, we know that renormalizability implies that the quantum Hamiltonian governingĥ µν is bounded from below [4,37] and, therefore, this dangerous tunneling should not occur. On the other hand, renormalizability also implies that the space of states is endowed with an indefinite metric (with respect to which the quantum mechanical "position" q and momentum p operators are self-adjoint) [4,37,47,48] . The presence of an indefinite metric, therefore, leads to the question: how can we define probabilities consistently?
The crucial observation here is that such metric is not the one that should be used in the Born rule to compute probabilities. To find the correct norm we start by a definition of observables that generalizes the one usually given in quantum mechanics to avoid any reference to a specific norm. We define "observable" any operator A with a complete set of eigenstates 8 {|a } [4]: for any state |ψ there is a decomposition |ψ = a c a |a for some coefficients c a . Moreover, we interpret |a as the state where A assumes certainly the value corresponding to a (the deterministic part of the Born rule). The correct way to approach this problem is to recall how experiments are performed. Experimentalists prepare a large number of times N the same state, so one is naturally led to consider the direct product of |ψ with itself N times: The next step is to define a frequency operator F a that counts the number of times N a there is the value a in the state |a 1 ...|a N : 7 Furthermore, the fatal runaways above such threshold give an (anthropic) rationale for a homogeneous and isotropic universe [45]. 8 One can show that the basic operators q and p as well as the Hamiltonian in the theory considered here have complete eigenstates at any order in perturbation theory.
Refs. [49,50] showed that in the sense that all coefficients in the basis |a 1 ...|a N of both F a |Ψ N and p a |Ψ N converge to the same quantities. This gives us a strong physical motivation for assuming the complete Born rule to compute probabilities: the probability of observing |a if the system is prepared in the state |ψ is given by p a . The born rule can be expressed in terms of a norm on the space of states, but this is not the indefinite metric mentioned above: Indeed, for any pair of states |ψ 1 , |ψ 2 , we denote the indefinite metric mentioned above with ψ 2 |ψ 1 and the positive norm which gives the Born rule in (2.9) is different and defined by The second equation above gives the definition of the "norm" operator P A . Eq. (2.9) clearly shows that all probabilities are positive and that they sum up to one at any time (the theory is unitary) 9 .
A remaining issue that has been pointed out is that higher-derivative theories may generate a violation of micro causality in its decay processes [51]. However, whenever M 2 < H, where H is the Hubble rate during inflation 10 , the width of the massive spin-2 field, Γ 2 , (the only one that could lead to micro acausality) is always much smaller than H as Γ 2 M 2 in this case [54][55][56]. Therefore, these potential acausal processes are actually quickly diluted by the expansion of the universe. In the opposite case, M 2 > H, the Weyl-squared term does not contribute to the experimentally observable quantities as it essentially decouples [42]. Therefore, the Weyl-squared term effectively respects the causality principle Having addressed the potential issues raised by the W 2 , we now turn to the analysis of the quasi-conformal scenario.

Hilltop-inflation realization
We describe here a simple realisation of this scenario. The inflaton-gravity (IG) sector features a single scalar field φ non-minimally coupled to gravity with Lagrangian (3.1) 9 A different proposal for unitarizing this theory was found in [52,53]. In this work we do not adopt this alternative approach because it is not (currently) developed enough to provide us with predictions for the CMB observables. On the other hand, for the approach discussed here, the predictions for such observables have been found in [42]. 10 Note that H satisfies the observational bound H < 2.7 · 10 −5M P (95% CL) [1,2] so H M P and M 2 < H implies f 2 1 (see Eq. (2.6)).
where the non-minimal coupling is taken close to the conformal value ξ ≈ −1/6 (or ζ = ξ + 1/6 ≈ 0). The potential will be chosen of the form We neglect the present-day cosmological constant Λ here as completely negligible compared to the early-universe scales. Such a potential is the simplest one that guarantees the presence of a maximum at a moderate field value (zero in this case), which in turn can lead to hilltop inflation [35,57], and a minimum where φ can lie after inflation. This is what one needs to have inflation for ξ ≈ −1/6, indeed, in order for the effective Planck mass M 2 P + ξφ 2 ≈ M 2 P − φ 2 /6 to remain real, inflation should occur at moderate field values (not exceeding √ 6M P ). The non-minimal coupling −ξφ 2 R/2 can be eliminated through the conformal transformation Note that, in order for this field redefinition to be well-defined, we have to require |φ| <M P / |ξ| when ξ < 0, so in particular this constraint applies for ξ = −1/6. This is nothing but the requirement of having a real effective Planck mass discussed above. In the frame obtained with the transformation above (called the Einstein frame) gravity is canonically normalized. Indeed, the action (after the conformal transformation) becomes The non-canonical kinetic term can be made canonical through the field redefinition φ = φ(χ) defined by with the conventional condition 11 φ(χ = 0) = 0. One can find a closed expression of χ as a function of φ [58]: (3.10) 11 For a generic number N s of scalar fields φ a (a = 1..., N s ) Eq. (3.1) gets replaced by and going to the Einstein frame one finds where The function χ(φ) can be inverted to obtain φ(χ) because Eq. (3.5) implies dχ/dφ > 0. Thus, χ feels a potential (3.11) The label E recalls that this is the Einstein frame potential, as opposed to the so-called Jordan frame one, Eq. (3.2). It is difficult to compute analytically φ(χ) for a generic value of ξ. But for ξ = −1/6 we have simply and therefore This is the reason why an inflationary scenario with a quasi-conformally coupled inflaton can be analytically studied. In the following, as an example, we will present the analytic inflationary formulae in the case of hilltop inflation, with Jordan frame potential (3.2), but an analytic treatment can always be carried out whenever one has a closed form for the Jordan frame potential. By inserting (3.13) inside (3.11) one finds (3.14) Note that this potential does not have a singularity for any value of χ although V (φ)/Ω 4 (φ) diverges at φ = √ 6M P (having set ξ = −1/6). The reason why this happens is because the function φ(χ) in (3.13) reaches the value √ 6M P only for χ → ∞, which is due to the fact that the singularity of V (φ)/Ω 4 (φ) at φ = √ 6M P is compensated by a singularity of K(φ) at the same point. There is an analogous effect for generic values of ξ < 0 [59]. Also note that V E has two stationary points for v < √ 6M P : a maximum at χ = 0 and a vanishing minimum at χ = √ 6M P arctanh v/ √ 6M P ; for v > √ 6M P the latter minimum disappears and the only stationary point is χ = 0, which becomes a minimum. However, as commented before, for v > √ 6M P the effective Planck mass M 2 P − v 2 /6 becomes imaginary and the conformal transformation in (3.3) is ill-defined, thus one must exclude this case. Fig. 1 gives V as a function of φ and V E as a function of χ for two values of v. V and V E are very similar for small enough field values. However, in the opposite limit they differ.

Computing inflationary observables
Let us take as a starting point the Friedmann-Robertson-Walker (FRW) metric and (3.9) Generically, the Ricci scalar of the field metric K ab is not zero so K ab is not flat. (k = 0, ±1). Correspondingly, we assume the scalar field χ to be a function of t only. Then the Einstein equations and the scalar equations imply the following equations for a(t) and the spatially homogeneous field χ(t)χ where a dot denotes a derivative with respect to t and a prime is a derivative with respect to χ. From Eqs. (3.16) and (3.17) one can derive (3.18), which is, therefore, dependent. During inflation we assume that the energy density is dominated by the scalar field χ and set k = 0; moreover, the following slow-roll parameters should be small For ξ = −1/6 one obtains the analytic expressions In the slow-roll approximation the EOMs (3.16)-(3.18) arė The scalar spectral index n s and the tensor-to-scalar ratio r are given by [42] n s = 1 − 6 + 2η, where r E is the tensor-to-scalar ratio when the pure gravitational Lagrangian is the Einstein-Hilbert one [42]. The difference between r and r E is due to W 2 . The Weyl-squared term (if M 2 < H) also leads to an isocurvature mode with power spectrum P B , such that the ratio is a factor 3/16 smaller than r E [42]. Eqs. (3.27) and (3.28) are valid in a generic inflationary model [42] and so they apply to the natural-inflation model that will be studied in Sec. 4 too. A general observational bound on r have been set in Ref. [45] (see Fig. 7 in that article), which translates into a bound on r E because, as mentioned above, r = 3r E /16. Moreover, we find the following analytic expressions n s (χ) = 12 sech 2 and, since in the slow-roll approximation H is given by .

(3.30)
The number of e-folds N e as a function of the field χ is given by where, in the slow-roll approximation, the function N is and χ end is the field value at the end of inflation, determined by (χ end ) = 1: we find an analytic expression for χ end too, The integration constant in the expression of N (χ) above plays no role because what matters is the difference between N computed in different points, Eq. (3.31). Also, one obtains a closed form for N (χ): where "const" represents an arbitrary constant with respect to χ. By requiring that the measured power spectrum [2], is reproduced for a field value corresponding to an appropriate value of e-folds one fixes λ φ . In Fig. 2 it is shown that there are initial field values χ for which we can have N e ∼ 60 e-folds and obtain n s and r E (and, therefore, r) in good agreement with the bounds in [1,2], preserving the slow-roll conditions, 1, η 1. Since r E < 0.1 [45] the isocurvature power spectrum parameterized by r also satisfies the bounds in [1,2]. By choosing λ φ ≈ 7 × 10 −12 one obtains the measured P R in Eq. (3.35). A numerical analysis shows that, generically, one has a good agreement with the inflationary observables if v is close enough to √ 6M P . Finally, we have explicitly checked in this hilltop implementation of the quasi-conformal scenario that the bounds in (2.7) and (2.8) are satisfied during the whole cosmology and, therefore, the possible runaways due to the Ostrogradsky theorem are avoided.

Higgs naturalness
In the previous section we have seen that successful inflation can occur if the vacuum expectation value (VEV) v in the inflaton sector is around (but not exceeding much) the Planck scale.
The fluctuation of φ around v acquires a mass M φ ≡ 2λ φ v, which, for the typical orders of magnitude v ∼M P and λ φ ∼ 10 −11 (obtained from the amplitude of the curvature fluctuations) is of order 10 13 GeV. Given that M φ lead to a quantum correction to the Higgs mass squared 36) Figure 2: The slow-roll parameters and η defined in Eq. (3.19), the number of e-folds N e given in (3.31) and n s and r (computed with Eqs. (3.25) and (3.26)) as a function of χ. We have set v = 2.43M P (the quantities plotted do not depend on λ φ ).
where λ hφ is the quartic portal coupling between φ and the Higgs (appearing in the potential as +λ hφ φ 2 h 2 , where h is the physical Higgs field), one obtains an upper bound on λ hφ to ensure the naturalness of the Higgs mass, which requires δM 2 h M 2 h [60]. For the orders of magnitude mentioned above we obtain λ hφ 10 −20 .
This is a tiny number, but one should bear in mind that λ hφ is the only coupling between the inflaton and the observable sector (besides gravity) and, therefore, a small λ hφ is not destabilized by quantum corrections. Gravity can take part in the running of λ hφ , but the agravity scenario studied in [6] softens gravity above a certain energy scale and its contribution turns out to be negligibly small when the condition in (2.5) is used. The conditions in (2.5) and (3.37) ensure that the shift symmetry of the Higgs field is good enough to protect its mass from gravity and inflaton quantum corrections, respectively, leading to a technically natural [7] Higgs mass.

Complexity of UV-completions
The hilltop inflationary model provides us with an existing proof of quasi-conformal models of inflation. Given that such scenario is motivated by possible UV-completions of Einstein gravity one eventually want to embed hilltop inflation in a field theory without a high-momentum cutoff. In Refs. [61][62][63] it was shown that renormalizable field theories can be made asymptotically free with respect to all couplings (including the quartic couplings λ abcd ) if embedded in a theory with a (semi-)simple gauge group and with a complex field content. The field content can be made simpler in asymptotically safe extensions of the SM, which, however, present some challenges [64,65]. The scenario described in Sec. 2 is general enough to embed this SM extensions. Such an embedding for the specific case of hilltop inflation, although interesting, goes beyond the scope of the present work because, as we will see in Sec. (4), there are simple and elegant asymptotically free theories in which the inflaton is identified with a composite scalar. In these theories all fundamental scalars feature a quasi-conformal non-minimal coupling, but, as will be discussed in Sec. 4, the inflaton, being a pseudo-Goldstone boson, has a negligible (and thus non-conformal) non-minimal coupling.

Natural-inflation realizations
Another way of implementing the quasi-conformal scenario for the early universe is identifying the inflaton φ with the pseudo-Goldstone boson associated with the breaking of a global symmetry. In the context of the non-conformal Einstein-gravity case, this scenario has been proposed in [36] and is known as "natural inflation" because the potential of a pseudo-Goldstone boson is naturally flat thanks to Goldstone theorem. This good feature is preserved when one constructs a naturalinflation model in the quasi-conformal scenario, but, as we will see, some predictions differ with respect to the Einstein-gravity case, such that one can distinguish between these two proposals.
First note that (pseudo-)Goldstone bosons do not need to have a Weyl-invariant non-minimal coupling to gravity in the (even exact) conformal scenario; their ξ-couplings can be different from −1/6. Indeed, if Φ is a scalar field in some representation of the spontaneously broken group G S , its non-minimal coupling is proportional to ξ Φ |Φ| 2 R, where |Φ| 2 is invariant under G S , and the Goldstone boson (being morally some phase of Φ) does not feature non-minimal couplings.
The inflaton potential of natural inflation V N (φ) can be written as [36] V N (φ) = Λ 4 1 + cos N φ f , (4.1) where Λ and f are two energy scales and N is a natural number. Given that V N is even and periodic with period 2πf /N we restrict ourselves to the interval 12 φ ∈ [0, πf /N ]. Like in Sec. 3.1, 12 As far as inflation is concerned, it would be equivalent to consider the potential with a different sign, , and to restrict ourselves to the interval φ ∈ [πf /N, 2πf /N ]. Therefore, the two signs (considered in [36]) are both covered here.
we take an FRW metric with k = 0 to describe inflation, when the energy density of the universe was dominated by the inflaton contribution.
The slow-roll quantities are given by In the slow-roll approximation the spectral index n s and the tensor-to-scalar ratio in Einstein gravity are given by As discussed before, adding the Weyl-squared term leaves n s unchanged, but gives the following tensor-to-scalar ratio [42] where H is the inflationary Hubble rate, which, in the slow-roll approximation, is given by H 2 (φ) = V N (φ)/3M 2 P . Therefore, for natural inflation . (4.7) Note that when H M 2 the tensor-to-scalar ratio is much smaller than the one predicted within Einstein gravity, as clear from (4.6). This is interesting because natural inflation is on the verge of being excluded when implemented in Einstein gravity [2]. Therefore, the Weyl-squared term can rescue this well-motivated model when H > M 2 .
The number of e-folds N e as a function of the field φ is here given by where, in the slow-roll approximation, the function N is and φ end is defined by (φ end ) = 1, which gives (4.10) Figure 3: The slow-roll parameters, and η, the spectral index n s , the tensor-to-scalar ratio r and the curvature power spectrum P R in natural inflation as a function of f . We have set Λ ≈ 6 × 10 −3M P to fit the observed value of the curvature power spectrum P R and chosen N = 1 and f 2 = 10 −8 (the value of f 2 influences the plot of r only). The bounds from the latest Planck analysis from 2018 [2] are also shown.
One also obtains an analytic expression for N (φ): Finally the curvature power spectrum (at horizon exit) is and its expression for natural inflation is (4. 13) We see that the main inflationary predictions of the model can be determined analytically, like for the fundamental conformally coupled inflatons discussed before.
In Fig. 3 we show the predictions of the model for the observables n s , r and P R together with the behavior of and η to show the validity of the slow-roll approximation. In the bottom plot on the left one can clearly see that the Weyl-squared term allows perfect compatibility between natural inflation and the latest CMB data, including the value of r. We also see that for N = 1 the scale f should be around f ≈ 6.6M P while Λ ≈ 6 × 10 −3M P . Since r E < 0.1 the isocurvature power spectrum parameterized by r (see Eq. (3.28)) also satisfies the most recent bounds (see Fig. 7 in [45]).
Like for the hilltop case of Sec. 3, we have also explicitly checked for the natural-inflation implementation that the bounds in (2.7) and (2.8) are satisfied during the whole cosmology and, therefore, the possible runaways due to the Ostrogradsky theorem are avoided.
Note that the non-minimal coupling between the Ricci scalar and φ might be generated by quantum corrections. However, it should be small, of order of the scale N Λ 2 /f at which the shift symmetry of the φ field is broken divided by the Planck mass: one has to divide byM P because the non-minimal coupling is dimensionless and because the non-minimal coupling should disappears when gravity is decoupled. This gives a negligibly small 13 non-minimal coupling of the inflaton, that is ξ φ ∼ Λ 2 /(fM P ) ∼ 10 −5 for an order one N , Λ ≈ 6 × 10 −3M P and f ≈ 6.6M P .

Asymptotically free UV-completion
An asymptotically free UV-completion of this scenario can be obtained as follows. Take a version of QCD with a confinement scale around the Planck scale. Introduce three flavors of tilde-quarks 14 : q = {ũ,d,s}. By using the analogy with QCD, we identify f with the decay constant ofπ 0 . We take f ≈ 6.6M P , as suggested by the results in Fig. 3. Like in QCD the strong dynamics forms condensates with a typical scale 15 κ ∼ q q . We will denote the gauge group with SU(3) I to emphasize that it is the group related to the inflaton sector. From the knowledge of QCD we obtain κ ≈ 3f 3 . (4.14) We add now, just like in QCD, quark masses such that the axial part of the global SU(3) f flavor group (which rotates {ũ,d,s}) is explicitly broken. Then we have the following spectrum of pseudo-Goldstone bosons which can be found with EFT methods [66]. Now, by choosing an inverted hierarchy mũ md, ms we obtain that the lightest pseudo-Goldstone boson isK 0 ; this gives just one inflaton field,K 0 , and avoids the isocurvature bounds of Planck. Note that from (4.1) the mass of the inflaton is m φ = N Λ 2 /f and, using the numerical values in Fig. 3 (Λ ≈ 6 × 10 −3M P and N = 1 ) as well as f ≈ 6.6M P , one obtains the typical value m φ ≈ 5.4 × 10 −6M P . The hierarchy Λ f corresponds to choosing the Lagrangian masses of the lightest tilde-quarks 16 mq to be really tiny in units of M P : indeed, from m 2 φ ∼ 4mqκ/f 2 (see (4.15)-(4.18)), f ≈ 6.6M P , m φ ≈ 5.4 × 10 −6M P and (4.14) one finds mq ∼ m φ (m φ /12f ) ∼ 10 −13M P .

Higgs naturalness
The naturalness of the Higgs mass is elegantly achieved in these realizations. One can consider an asymptotically free theory with a gauge group where G h is the group factor under which the Higgs is charged. For example, in [63] it was shown that the SM can be embedded in an asymptotically free theory with the trinification gauge group G 333 ≡ SU(3) L ×SU(3) C ×SU(3) R . In this case the total gauge group will be Other examples can be found by looking at the other asymptotically free matter sectors found in [61,62]. Even more examples can be generated by considering identical copies of these asymptotically free matter contents. As long as the couplings of these visible sectors (charged under G h ) with the inflaton sector are small the Higgs mass is protected from large radiative corrections by an approximate shift symmetry acting on the Higgs field. As discussed in Sec. 2, this approximate symmetry is respected by gravitational interactions if (2.5) is enforced.

Reheating
Given that there is an asymptotically free embedding of natural inflation we study further the early-universe implications of this scenario. Here we describe how reheating can occur.
The Lagrangian in the Einstein frame can be computed through the quantity Ω 2 defined in (3.8). Note that the fundamental scalars φ a undergo quantum fluctuations [68] of order H/(2π) M P during inflation, therefore, assuming all of them of this order of magnitude one finds (4.21) Taking now H 2π ≈ 10 13 GeV (4.22) (having used the setup of Fig. 3) gives Ω 2 ≈ 1 − 3 × 10 −12 N s . So Ω 2 is very close to one whenever N s 10 12 . We assume this is satisfied. Then, as far as reheating is concerned, the Einstein-frame Lagrangian is basically equal to the Jordan frame one, but with negligibly small non-minimal couplings (see Eqs. (3.7), (3.8) and (3.9) and, for the matter fields with non-zero spin, see the Einstein-frame Lagrangian in Ref. [4]).
During inflation the inflaton energy density ρ I is very high; in particular it is much higher than the energy density of the fundamental scalars: However, ρ I decreases with time. When m φ ȧ(t)/a(t) the inflaton oscillates. Taking m φ ≈ 5.4 × 10 −6M P as in the example of Fig. 3, the inflaton oscillations occur when the Hubble rate is say a couple of orders of magnitude below, H osc ∼ 5 × 10 −8M P , and the corresponding inflaton energy density is ρ osc To reheat the visible universe we need eventually During inflation, as mentioned above, φ a undergo quantum fluctuations of order H/(2π). If λ abcd 1, which we assume and is typically true in asymptotically free theories, this value is kept for a long time, while = 1 at the end of inflation, so the Hubble rate decreases more rapidly and ρ I quickly reaches ρ osc I . Assuming, for example, the λ abcd to be all of the same order,λ, (4.25) requiresλ For example, forλ ∼ 10 −2 this means 17 N s 10 2 . One should keep in mind, however, that some of the λ abcd can be much smaller or vanishing.
When (4.26) holds, assuming the instantaneous reheating (equating the scalar to the thermal energy density), one obtains the reheating temperature , (4.28) where the number of relativistic degrees of freedom g * has to be large because of (4.27). For the numerical values used in Fig. 3  which is typically large enough and in particular much larger than the electroweak scale. Note that what allows us to heat the universe is the presence of many weakly coupled scalars, which have sizable couplings to the observed particles and undergo large quantum fluctuations.

Conclusions and outlook
In this paper we explored the observable predictions of the quasi-conformal scenario, in which the dimensionless couplings in the action respect (at least approximately) Weyl symmetry. This behavior is suggested by the UV complete QG, which can hold up to infinite energy if all couplings reach a conformal fixed point at infinite energies [3]. QG can also render the Higgs mass (technically) natural by allowing a good enough shift symmetry of the Higgs field even in the presence of gravity.
We have included in the discussion the Weyl-squared term and reviewed and extended the arguments in favour of its viability in Sec. 2. Basically the discussion can be addressed at two levels. At the classical level, it was shown [45] that the possible Ostrogradsky runaways are avoided if the typical energies are below a threshold (given in (2.7) and (2.8)), which is high enough to describe the whole cosmology. At quantum level unitarity is restored by computing probabilities by an appropriate positive norm, which is singled out by an experimental approach to the definition of probabilities. Moreover, the possible acausal effects are extremely diluted by the expansion of the universe.
The predictions of the Weyl-squared term and the quasi-conformal non-minimal coupling of fundamental scalars have been worked out and compared with the available data of the early universe. We considered two implementations of the quasi-conformal option.
In the first one, discussed in Sec. 3, the inflaton is identified with a fundamental scalar with a moderate (< √ 6M P ) field excursion, such that the effective Planck mass never becomes imaginary. An example is hilltop inflation, which agrees well with the most recent data provided by the Planck collaboration. We pointed out, however, the complexity required to UV complete models of this sort.
This brought us to the second implementation, which we investigated in Sec. 4: a pseudo-Goldstone boson as the inflaton (natural inflation), which provides a very good rationale for the flatness of the potential. In this case the non-minimal coupling can be nearly zero, rather than conformal, because a pseudo-Goldstone boson has very small non-derivative interactions. In Einstein-gravity implementations of natural inflation, the prediction of r is rather large and on the verge of being excluded. On the other hand, in the presence of the Weyl-squared term, r is drastically reduced for a natural Higgs mass. This can revive natural inflation. Another advantage of this realization is the possibility to be elegantly UV-completed. We have shown that an asymptotically free QCD-like sector can provide an inflaton of this type: it can be identified with the composite scalar analogous to the K 0 meson in ordinary strong interactions. Furthermore, we have also shown that a satisfactory reheating can occur in this case.
We have checked, both in the hilltop and in the natural-inflation case, that the bounds (2.7) and (2.8) on the energies to avoid the Ostrogradsky instabilities are satisfied.
Furthermore, in both cases the main inflationary formulae can be analytically computed, potentially providing simple benchmark models of inflation.
We take advantage of the concluding section to point out that this scenario, appropriately implemented to include a UV-complete sector containing all SM particles, as discussed in Sec. 4.2, may solve all phenomenological problems of the SM. For example, dark matter, baryogenesis and neutrino masses might be due to right-handed neutrinos [69,70] (which in some models with UV fixed points are necessary [63]) and the electroweak vacuum stability may be achieved thanks to the extra scalars generically present in the theory (see Refs. [71][72][73][74] for the stabilization through extra scalars). Therefore, we hope that the results of this paper will lead to scientific activities along these lines.