Origin of inflation in CFT driven cosmology: $R^2$-gravity and non-minimally coupled inflaton models

We present a detailed derivation of the recently suggested new type of hill-top inflation [arXiv:1509.07270] originating from the microcanonical density matrix initial conditions in cosmology driven by conformal field theory (CFT). The cosmological instantons of topology $S^1\times S^3$, which set up these initial conditions, have the shape of a garland with multiple periodic oscillations of the scale factor of the spatial $S^3$-section. They describe underbarrier oscillations of the inflaton and scale factor in the vicinity of the inflaton potential maximum, which gives a sufficient amount of inflation required by the known CMB data. We build the approximation of two coupled harmonic oscillators for these garland instantons and show that they can generate inflation consistent with the parameters of the CMB primordial power spectrum in the non-minimal Higgs inflation model and in $R^2$ gravity. In particular, the instanton solutions provide smallness of inflationary slow-roll parameters $\epsilon$ and $\eta<0$ and their relation $\epsilon\sim\eta^2$ characteristic of these two models. We present the mechanism of formation of hill-like inflaton potentials, which is based on logarithmic loop corrections to the asymptotically shift-invariant tree level potentials of these models in the Einstein frame. We also discuss the role of $R^2$-gravity as an indispensable finite renormalization tool in the CFT driven cosmology, which guarantees the non-dynamical (ghost free) nature of its scale factor and special properties of its cosmological garland type instantons. Finally, as a solution to the problem of hierarchy between the Planckian scale and the inflation scale we discuss the concept of a hidden sector of conformal higher spin fields.

Abstract We present a detailed derivation of the recently suggested new type of hill-top inflation originating from the microcanonical density matrix initial conditions in cosmology driven by conformal field theory (CFT) [arXiv:1509.07270]. The cosmological instantons of topology S 1 × S 3 , which set up these initial conditions, have the shape of a garland with multiple periodic oscillations of the scale factor of the spatial S 3 -section. They describe underbarrier oscillations of the inflaton and scale factor in the vicinity of the inflaton potential maximum, which gives a sufficient amount of inflation required by the known CMB data. We build the approximation of two coupled harmonic oscillators for these garland instantons and show that they can generate inflation consistent with the parameters of the CMB primordial power spectrum in the non-minimal Higgs inflation model and in R 2 gravity. In particular, the instanton solutions provide smallness of inflationary slow-roll parameters ǫ and η < 0 and their relation ǫ ∼ η 2 characteristic of these two models. We present the mechanism of formation of hill-like inflaton potentials, which is based on logarithmic loop corrections to the asymptotically shift-invariant tree level potentials of these models in the Einstein frame. We also discuss the role of R 2 -gravity as an indispensable finite renormalization tool in the CFT driven cosmology, which guarantees the non-dynamical (ghost free) nature of its scale factor and special properties of its cosmological garland type instantons. Finally, as a solution to the problem of hierarchy between the Planckian scale and the inflation scale we discuss the concept of a hidden sector of conformal higher spin fields.

Introduction
The problem of initial conditions in cosmology [1,2] as a source of the inflationary scenario starts attracting attention again. A considerable success in explaining the CMB data [3,4] by amplification of primordial quantum cosmological pertubations [5] in the model of R 2 gravity [6] and the non-minimal Higgs inflation model [7,8] is called in question by objections against the origin of inflation [9]. These objections are based on the statement that distributions of initial conditions violate a widely accepted naturalness assumption that all forms of inflaton energy (kinetic, gradient and potential) should initially have the same Planckian scale magnitude [10,11]. Three known sources of these distributions -pure no-boundary [1] and "tunneling" [2] quantum states of the Universe and the Fokker-Planck equation for coarse-grained cosmological evolution [12] -in their turn suffer from intrinsic difficulties associated with a missing clear canonical quantization ground, insufficient amount of generated inflation, anthropic (observer dependence) problems [13] and rather contrived multiverse measure problem [14], etc.
The goal of this work is to give a detailed derivation of a recently suggested new model of hill-top inflation [15] in the theory of CFT driven cosmology [16,17] which is likely to circumvent the difficulties of the above type. This theory represents the synthesis of two main ideas -new concept of the cosmological microcanonical density matrix as the initial state of the Universe and application of this concept to the system with a large number of quantum fields conformally coupled to gravity. It plays an important role within the cosmological constant and dark energy problems. In particular, its statistical ensemble is bounded to a finite range of values of the effective cosmological constant, it incorporates inflationary stage and is potentially capable of generating the cosmological acceleration phenomenon within the so-called Big Boost scenario [16,18]. Moreover, as was noticed in [19], the CFT driven cosmology provides perhaps the first example of the initial quantum state of the inflationary Universe, which has a thermal nature of the primordial power spectrum of cosmological perturbations [20]. This suggests a new mechanism for the red tilt of the CMB anisotropy, complementary to the conventional mechanism which is based on a small deviation of the inflationary expansion from the exact de Sitter evolution [5].
This setup has a clear origin in terms of operator quantization of gravity theory in the Lorentzian signature spacetime and is based on a natural notion of the microcanonical density matrix as a projector on the space of solutions of the quantum gravitational Dirac constraints -the system of the Wheeler-DeWitt equations [17,21]. Its statistical sum has a representation of the Euclidean quantum gravity (EQG) path integral [16,17] over metric g µν and matter fields Φ which are periodic on the Euclidean spacetime with a time compactified to a circle S 1 . As shown in [16,17], this statistical sum is approximately calculable and has a good predictive power in the gravitational model with the primordial cosmological constant Λ and the matter sector which mainly consists of a large number N of free (linear) fields φ conformally coupled to gravity -conformal field theory (CFT) with the action S CF T [ g µν , Φ ], Important point, which allows one to overstep the limits of the usual semiclassical expansion, consists here in the possibility to omit the integration over conformally noninvariant matter fields and spatially-inhomogeneous metric modes on top of a dominant contribution of numerous conformal species. Integrating them out one obtains the effective gravitational action S eff [ g µν ] which differs from (2) by S CF T [ g µν , Φ ] replaced with Γ CF T [ g µν ]the effective action of Φ on the background of g µν , On Friedmann-Robertson-Walker (FRW) background this action is exactly calculable by using the local conformal transformation to the static Einstein universe and well-known local trace anomaly. The resulting Γ CF T [ g µν ] turns out to be the sum of the anomaly contribution and free energy of conformal matter fields on the sphere S 3 at the temperature determined by the circumference of the compactified time dimension S 1 . 1 This is the main calculational advantage provided by the local Weyl invariance of Φ conformally coupled to g µν .
Physics of the CFT driven cosmology is entirely determined by this effective action. Solutions of its equations of motion, which give a dominant contribution to the statistical sum, are the cosmological instantons of S 1 ×S 3 topology, which have the Friedmann-Robertson-Walker metric with a periodic lapse function N (τ ) and scale factor a(τ ) -functions of the Euclidean time belonging to the circle S 1 [16]. These instantons serve as initial conditions for the cosmological evolution a L (t) in the physical Lorentzian spacetime. The latter follows from a(τ ) by analytic continuation a L (t) = a(τ * + it) at the point of the maximum value of the Euclidean scale factor a + = a(τ * ). The fact that these instantons exist only in the finite range of Λ implies the restriction of the microcanonical ensemble of universes to this range, which from the viewpoint of string theory can be interpreted as the solution of the landscape problem for stringy vacua.
As was originally mentioned in [18] this scenario can incorporate a finite inflationary stage if the model (2) is generalized to the case when Λ is replaced by a composite operator Λ(φ) = V (φ)/M 2 P -the potential of the inflaton φ slowly varying during the Euclidean and inflationary stages and decaying in the end of inflation by a usual exit scenario. The goal of this paper is to develop such a generalization which starts with the replacement of (2), S[ g µν , Φ ] → S[ g µν , φ, Φ ], by the action with the 1 Note that in previous works on trace anomaly applications in cosmology the contribution of the static Einstein universe was either restricted to the case of the vacuum state (Casimir energy) [22] or recovered for a particular value of integration constant in the stress tensor conservation law [23,24]. Here it follows directly from the density matrix prescription for the initial state of the Universe, which endows the Einstein universe with the radiation gas of CFT particles and makes our setting of the problem essentially different from previous studies. 3 inflaton field, whose potential V (φ) simulates the effect of the primordial cosmological constant. Under this replacement the restriction of the primordial cosmological constant range of the above type, Λ(φ) = V (φ)/M 2 P , becomes a selection of the range of φ or fixation of the initial conditions for inflation. These initial conditions is a principal goal of this paper.
As we will see CFT driven cosmology realizes these initial conditions in the form of the new type of hill-top inflation originating from the underbarrier oscillations of the inflaton φ and the scale factor a in the vicinity of local maxima of V (φ). Thus, these initial conditions are not the Planckian scale parameters introduced by hands from some ad hoc arguments of naturalness [10,9], but rather become derivable from microcanonical partition function as calculable parameters of the corresponding saddle point configurations.
Dynamical nature of the physical mode simulating the effect of the cosmological constant has also another important implication. The properties of instantons in the CFT driven cosmology critically depend on the vanishing of one of the conformal anomaly parameters (the coefficient of R). It is renormalization ambiguous, because it can be arbitrarily changed by the addition of the local R 2 term, and therefore can be renormalized to zero to provide needed properties of CFT initial conditions [16,15]. Quite interestingly, this finite renormalization goes beyond an ad hoc assumption of [16] and can be a result of inclusion of the Starobinsky R 2 model, one part of which would provide this renormalization in the anomaly part of the action and the other part would give a scalar mode simulating the effect of the inflaton field -effective cosmological term decaying at the exit from inflation.
The plan of the paper is as follows. In Sect.2 we overview the model of CFT driven cosmology with a fundamental cosmological constant. Sect.3 is devoted to the generalization of this model to the case of the dynamical inflaton field minimally coupled to gravity. It suggests the approximation of coupled harmonic oscillators for cosmological instantons in the regime of the Euclidean "slow roll", which set up initial conditions for inflation in physical spacetime, starting in the vicinity of the top of the inflaton potential. Sect.4 demonstrates formation of the relevant hill-like effective potential in the model of Higgs inflation with a strong non-minimal coupling of the inflaton to the scalar curvature. In Sect. 5 a similar mechanism is discussed for the Starobinsky model of R 2 -gravity which is shown to play a double role in the CFT scenario: finite renormalization of the quantum effective action, which provides special properties of the cosmological instantons, and generation of the dynamical inflaton feeding the CFT scenario with the hill-like potential. Sect.6 contains an attempt to solve the hierarchy problem in the CFT cosmology by using a hidden sector of numerous conformal higher spin fields and the concluding Sect.7 briefly discusses observational prospects of the model. Two appendices contain technical details of the harmonic oscillator approximation and its Euclidean "slow-roll" regime.

Model with a fundamental cosmological constant
For cosmology with S 1 × S 3 topology and FRW metric (4) its effective action where a ′ ≡ da/N dτ . The first three terms in curly brackets of (6) represent the Einstein action with a fundamental cosmological constant Λ ≡ 3H 2 (H is the corresponding Hubble parameter). The constant B is a coefficient of the contributions of the conformal anomaly and vacuum (Casimir) energy (B/2a) on a conformally related static Einstein spacetime mentioned in Introduction. This constant, expresses via the coefficient β of the Gauss-Bonnet term E = R 2 µναγ − 4R 2 µν + R 2 in the trace anomaly of conformal matter fields It should be emphasised here that this effective action is independent of the anomaly coefficients α and γ, because it is assumed that α is renormalized to zero by a local counterterm, This guarantees absence of higher derivative terms in (6) -non-ghost nature of the scale factor -and simultaneously gives the renormalized Casimir energy a partiular value proportional to B/2 = β/16π 2 M 2 P [25]. Both of these properties are critically important for the instanton solutions of effective equations. The coefficient γ of the Weyl tensor term C 2 µναβ does not enter (6) because C µναβ identically vanishes for any FRW metric.
Finally, F (η) is the free energy of conformal fields also coming from this Einstein space -a typical boson or fermion sum over CFT field oscillators with energies ω on a unit 3-sphere, η playing the role of the inverse temperature -an overall circumference of S 1 in the S 1 × S 3 instanton, calculated in units of the conformal time (8).
The statistical sum (1) is dominated by the solutions of the effective equation, δS eff /δN (τ ) = 0, which in the gauge N = 1 reads This is the modification of the Euclidean Friedmann equation by the anomalous B-term and the radiation term C/a 4 . The constant C here characterizes the sum of the Casimir energy and the energy of thermally excited particles with the inverse temperature η given by (8). This quadratic equation (12) can be solved forȧ 2 , to give a periodic oscillation of a between its maximal and minimal values a ± , provided that at a − we have a turning point with a vanishingȧ, which means that a 2 − > B. This inequality immediately yields the first two restrictions on the range of H 2 and C whereas the third one follows from the requirement of reality of turning points a ± , As shown in [16,17,19] the solutions of this integrodifferential equation 2 give rise to the set of periodic S 3 × S 1 instantons with the oscillating scale factorgarlands that can be regarded as the thermal version of the Hartle-Hawking instantons. The scale factor oscillates m times (m = 1, 2, 3, ...) between the maximum and minimum values (16), a − ≤ a(τ ) ≤ a + , so that the full period of the conformal time (8) is the 2mmultiple of the integral between the two neighboring turning points of a(τ ),ȧ(τ ± ) = 0, This value of η is finite and determines effective temperature T = 1/η as a function of G = 1/8πM 2 P and Λ = 3H 2 . This is the artifact of a microcanonical ensemble in cosmology [17] with only two freely specifiable dimensional parameters -the gravitational and cosmological constants. According to (17) these garland-type instantons exist only in the limited range of the cosmological constant Λ = 3H 2 [16]. In view of (18) and (19) they belong to the curvilinear domain in the two-dimensional plane of the Hubble constant H 2 and the amount of radiation constant C (each instanton being represented by a point in this plane). In this domain they form a countable, m = 1, 2, ..., sequence of one-parameter families -curves interpolating between the lower straight line boundary C = B − B 2 H 2 and the upper hyperbolic boundary C = 1/4H 2 . Each curve corresponds to a respective m-folded instantons of the above type. Therefore, the spectrum of admissible values of Λ has a band structure, each band being a projection of the m'th curve to the H 2 axis. The sequence of bands of ever narrowing widths with m → ∞ accumulates at the upper bound of this range H 2 max = 1/2B. The lower bound H 2 min -the lowest point of m = 1 family -can be obtained numerically for any field content of the model.
For a large number of conformal fields N, and therefore a large β ∝ N, the both bounds are of the order m 2 P /N. Thus the restriction (17) suggests a kind of 1/N solution of the cosmological constant problem, because specifying a sufficiently high number of conformal fields one can achieve a primordial value of Λ well below the Planck scale where the effective theory applies, but high enough to generate a sufficiently long inflationary stage.

Minimally coupled scalar field inflaton in CFT cosmology
Generalization of the action (2) to the effective Λ generated by a scalar field implies the transition from the effective minisuperspace action (6) to which leads to the replacement in the effective Friedmann equation (12) and to the origin of dynamical equation for the inflaton φ, The possibility of existence of periodic solutions of this system of equations qualitatively follows from counting its integration constants. Because of the fact that Eq.(23) is of the first order in derivatives, their number is three -the initial value of the inflaton φ 0 at some τ = τ 0 , its initial time derivativeφ 0 and the initial vale of the scale factor a 0 . Thus the solution has the form Requirement of periodicity with some period T implies the following equations From these three equations one can determine three initial condition parameters as functions of the period, which are locally unique provided the nonvanishing Ja- The period T does not remain a free parameter, because of the additional requirement. Analytical continuation of the Euclidean solution to the real solution in the Lorentzian time τ = τ * + it requires that time derivatives of φ(τ ) and a(τ ) at τ * vanish. The existence of such τ * is guaranteed for one of the variables, say , a 0 (T ) ) = 0, in view of its periodic motion away from φ 0 and back to φ 0 during the period of oscillations. Now if we demand that at the same moment the time derivative of a(τ ) is also zero, then we get the equation for T , which gives a locally unique solution T = T (τ * ). This counting of constants of motion guarantees that, if some approximate periodic solution satisfying the condition of Euclidean-Lorentzian junction exists, then perturbation theory will guarantee the periodicity of perturbed solution. Let us construct such an approximate solution which will serve as a lowest order of perturbation theory for the exact periodic trajectory.
To begin with, note that periodicity imposes a restriction on the shape of the potential V (φ) -integration of Eq.(24) over the period gives This implies that the gradient of V changes sign at some point φ 0 , so that inflaton oscillations occur in the vicinity of the extremum of V . Critical point is that these oscillations can take place only in the vicinity of the potential maximum in the underbarrier regime, which is shown schematically in Fig.2. This follows from a simple fact that our equations of motion are in Euclidean time, and they admit periodic oscillations of the scale factor only in the Euclidean time. Therefore, inflaton oscillations take place also in the underbarrier regime, which implies the vicinity of a potential maximum. Critical difference from the tunneling prescription of [1,2] here is that selection of the inflaton potential maxima follows not from the properties of the probability distribution, but from the requirement of a periodic solution for a saddle point of the partition function path integral. This solution gives rise to a new type of hilltop inflation which starts at the turning point φ * on the slope of the potential close to its maximum.

Harmonic oscillator approximation
It is obvious that solutions described above for a constant primordial Λ remain approximately true for the system of equations (23)-(24) when the inflaton field and its effective cosmological "constant" (22) do not vary much and their variations are not rapid. This is consistent with sufficiently slow variations of the scale factor, because from the inflaton equation of motion (24) it follows that the relevant rate of change of (22) is also small for smallφ andȧ, If we disregard the friction term in (24) and expand the potential in the vicinity of its maximum at φ 0 , small oscillations of φ will be nearly harmonic, with some amplitude ∆ φ = φ 0 − φ * , where φ * is a turning point shown on Fig.2.
A related quasi-harmonic behaviour of the scale factor is analytically available for solutions close to the upper hyperbolic boundary of the domain (19) of Fig.1. In a small vicinity of this boundary smallness ofȧ is guaranteed, because at C = 1/4H 2 the scale factor is pinched between the coincident a − and a + , a(τ ) = a − = a + . If we use the notations introduced in [27] then the proximity to the upper boundary is determined by smallness of ∆, ∆ ≪ 1. In these variables the admissible domain of solutions looks as a unit triangle (the quadrangle in terms of the quantities ε and d ≡ ∆/ε of [27]) and the scale factor oscillating between the turning points a ± can be parameterized by the variable y running between −1 and 1 according to In terms of y the equation (15) forȧ 2 reads aṡ All these equations hold also for the inflaton model. In this case the parameters ε and ∆, depending on H 2 , become the functions of time. However, if they are sufficiently slowly varying, time derivativesε and∆ can be disregarded. Then the above equation forȧ reads as the equation forẏ the variable y describes a harmonic oscillator with a slowly varying frequency ω and a unit amplitude, so that Conditions of applicability of this approximation, that should guarantee smallness of the friction term in the inflaton oscillator and a small change of ε and H during the period of oscillation 2π/ω, are Sinceȧ/a ∼ ω∆/2, these bounds lead to where we have introduced the notation In addition to the above bounds we also need a small rate of change of ε and ∆ necessary for the transition from Eq. (38) to Eq. (39). Smallness ofε is guaranteed by the bound (49) derived above, whereas a small |∆| ≪ ω∆ requires a much stronger bound, because∆ ≃ −(1/∆)Wȧ/a ∼ W ω ≪ ω∆ or Interestingly, the last bound -smallness of∆ -can be replaced by the opposite limit W ≫ ∆. For largė ∆ Eq.(39) does not hold, and oscillations of y become strongly anharmonic. However, as shown in Appendix A, the period of these oscillations remains the same 2π/ω and y behaves as y ≃ sin ωτ 2 provided the function (50) is slowly varying compared to the oscillations of the scale factor,Ẇ ≪ ωW . This leads to the extra bound on µ, µ ≪ ω, so that together with (48) the frequency of inflaton oscillations µ belongs to a limited range which is nonempty for assumed values of ∆ ≪ 1. Thus, the inflaton field and the scale factor represent two coupled quasi-harmonic oscillators with frequencies µ and ω. Periodicity of their motion implies their frequencies to be commensurable where m and n are some integer numbers. Therefore, from the definition (43) of the frequency ω it follows that Thus the approximation of two coupled oscillators works for a wide range of instantons lying close to the upper hyperbolic boundary even when the parameter ∆ is rapidly varying in time -as fast as the the scale factor and violating the bound (52). The lower bound for µ in (53) implies fast oscillations of a with a slower motion of φ, m ≫ n ≥ 1. This conclusion is also confirmed by considering the slow roll smallness parameters for the inflation stage originating from the transition through the turning point φ * .

Inflation stage and its slow roll parameters
Hill-top inflation histories in Lorentzian time, φ L (t) and a L (t), originate by analytic continuation of the Euclidean solutions (34) and (45) to τ = 2mπ/ω + it, where m ≫ 1 is the number of oscillations of the scale factor in the garland instanton. For a small time t the linearized Lorentzian solutions read while at later times nonlinear effects start dominating, so that the Lorentzian version of the nonlinear equation (23) with Λ replaced by the dynamical inflaton energy density, enters the stage. This equation can be rewritten in the manifestly Friedmann form with the effective Planck mass M eff (ρ) depending on the full matter density ρ, and ρ together with the inflaton energy density ρ φ includes the primordial radiation of the CFT cosmology [18] 8 The further evolution for large t consists in the fast quasi-exponential expansion during which the primordial radiation gets diluted, the inflaton field decays by a conventional exit scenario and goes over into the quanta of conformally non-invariant fields produced from the vacuum. They get thermalized and reheated to give a new post-inflationary radiation. Thus the primordial radiation of numerous conformal species and the inflaton energy get replaced with the radiation of non-conformal particles, ρ → ρ rad , which give rise to the radiation dominated Universe. With ρ rad dropping down below the sub-Planckian energy scale, ρ rad ≪ M 4 P /β, the effective Planck mass (60) tends to its Planckian value, M eff (ρ) → M P , and one obtains a standard general relativistic inflationary scenario for which the initial conditions were prepared by our CFT garland instanton [18,19].
The parameters of this inflation scenario are determined by the properties of the inflaton potential at the nucleation point φ * = φ(τ * ). For the quadratic potential (32) at τ * we have φ * = φ 0 − ∆ φ and V * ≡ V (φ * ) = 3M 2 P H 2 in view ofφ * = 0. Therefore the expressions for inflationary slow roll parameters at this point read Note that they are related by the equation which implies for ∆ φ ∼ M P the typical relation η * ∼ ǫ 2 * for slow roll parameters in the Starobinsky model [6] or in the model with a non-minimally coupled inflaton [26]. With H given by (55) the second slow roll smallness parameter (62), fails to be small unless Bµ 2 ≪ 1 and m 2 ≫ n 2 . This gives additional ground for the scale factor oscillations to be much faster than those of the inflaton, m ≫ n, and we consider this limit below.

Fast oscillations of the scale factor
For m ≫ n ≥ 1 the number m of the scale factor oscillations during the full period means that we consider m-fold instanton garland for which the energy scale is very close to the upper bound of the range (17) [16] The Euclidean solutions of the CFT cosmology in this limit can be called slow-roll ones, because the rate of change of the scale factor is much higher than that of the inflaton field. The Euclidean version of the slowroll regime is, however, rather peculiar, because in contrast to the Lorentzian case with monotonically changing variables here the scale factor and inflaton are oscillating functions of time. The details of these solutions including, in particular, the derivation of this asymptotics for H 2 (first given in [16]) are presented in Appendix B. Here we give a simplified overview of these solutions and their relation to the conventional slow roll parameters of inflation in Lorentzian theory.
Comparison of (66) with (55) implies that Since n ≥ 1 the lower bound on m is exponentially high, and the ratio n/m in (65) becomes exponentially small, As a result the solution becomes very close to the upper quantum gravity scale -the cusp of the curvilinear triangle on Fig.1 and the corresponding slow roll smallness parameter (65) expresses in terms of the quantity Bµ 2 = B| V ′′ * |, In view of the known CMB data for n s = 1−6ǫ * +2η * ≃ 1 + 2η * ≃ 0.96, this quantity is thus supposed to be also very small, Bµ 2 ∼ 0.01. Now, bearing in mind that we have the admissible range of ∆, ∆ ≤ 2 Bµ 2 e −2π/ √ Bµ 2 . This range is, however, further reduced by the requirement of the harmonic oscillator approximation (40), ∆ ≪ ε, and finally reads This bound is stronger than the requirement of a valid harmonic oscillator approximation (48), µ ≫ ω∆, or Bµ 2 , because in the full range of values of Bµ 2 the last bound is higher than (72) by the factor of Bµ 2 2 e π/ √ Bµ 2 ≥ π 2 e 2 2 ≫ 1. This establishes the range of the amplitude of scale factor oscillations ∆. To match with the slow roll smallness parameter η * ∼ −0.01 it should be exponentially Below we will estimate the bound on the first smallness parameter ǫ * in the CFT cosmology modelling initial conditions for the non-minimal Higgs inflation. It requires the knowledge of the amplitude of the inflaton field oscillations ∆ φ . As we will see, this amplitude will have a typical sub-Planckian value, so that a typical relation ǫ * ∼ η 2 * characteristic of the Starobinsky model or the model of the non-minimal Higgs inflation will hold and signify that ǫ * adds a negligible contribution to the CMB spectral parameter and provides a very small tensor to scalar ratio r = 16ǫ * .

Non-minimal Higgs inflation: the mechanism of hill like inflaton potential
In what follows we want to advocate that the CFT cosmology can serve as a source of initial conditions for the non-minimally coupled Higgs inflation model [7,28] which, together with the Starobinsky model of R 2inflation [6], is considered as one of the most promising models fitting the CMB data [3,4]. There is a twofold reason for that because, firstly, the Higgs inflation model at the quantum level has a natural mechanism of forming a hill-top potential and, secondly, it provides a relation ε * ∼ η 2 * ≪ |η * |, η * < 0, which establishes a strong link between the observable value of the CMB spectral parameter n s = 1 − 6ǫ * + 2η * ≃ 1 + 2η * ≃ 0.96 (74) and the value of the Higgs mass discovered at LHC [8,29,30,31]. This relation, as we will see below, will be provided by the bound on the amplitude of the inflaton oscillations ∆ φ in the underbarrier regime. The inflationary model with a non-minimally coupled Higgs-inflaton H, ϕ 2 = H † H, as any other semiclassical model, has a low-derivative part of its effective action (appropriate for the inflationary slow-roll scenario), Its coefficient functions contain together with their treelevel part the logarithmic loop corrections and include the dependence on the UV normalization scale µ, Numerical coefficients A, C, F are determined by contributions of quantum loops of all particles and represent beta functions of the corresponding running coupling constants -quartic self-coupling λ, non-minimal coupling of the Higgs field to curvature ξ -and anomalous dimension of ϕ. In the one-loop approximation these logarithmic corrections comprise the Coleman-Weinberg potential for the Higgs inflaton and the relevant correction to the non-minimal curvature coupling. When ξ ≫ 1 the inflationary stage and the corresponding CMB parameters critically depend only on A and the part of C linear in ξ, C = 3ξλ [8,28]. Other coefficients, the normalization scale µ inclusive, are irrelevant in the leading order of the slow roll expansion. Inflation and its CMB are easy to analyze in the Einstein frame of fieldsĝ µν , φ, in terms of which the action (75) has a minimal coupling of the inflaton to curvature, U = M 2 P /2, a canonical normalization of the inflaton field,Ĝ = 1, and a new inflaton potential, These Einstein frame fields are related to the Jordan frame of Eq.(75) by the equationŝ Due to the presence of leading logarithmic terms the Einstein frame potential starts decreasing to zero for large values of ϕ, In this limit the original Jordan frame field ϕ and the Einstein frame field φ are related by ϕ ≃ MP √ ξ e φ/ √ 6MP , so that this asymptotics looks likeV (φ) ≃ π 2 √ 6 9λξ 2 φ AM 5 P → 0, φ → ∞. Of course, this behavior cannot be extrapolated to infinity, because semiclassical expansion fails at transplanckian energies, but the maximum of the potential, which is reached at someφ,V ′ (φ) = 0, ϕ 2 ≃ µ 2 exp 16π 2 /3λ , corresponds for ξ ≫ 1 to a small subplanckian value of energŷ Therefore the potential starts bending down in the domain where the semiclassical expansion is still applicable and where it acquires a shape of suitable for our hill-top inflation scenario.
In fact a similar qualitative behavior holds in any order of loop expansion, because in the leading logarithm approximation of l-th loop order both V (ϕ) and U (ϕ) grow like l-th power of the logarithm while their ratio V /U 2 in the Einstein frame potential (80) decreases like (ln(ϕ/µ)) −l . This property was confirmed numerically within RG resummation of leading logarithms in [28] in the model of non-minimal Higgs inflation, which is illustrated by the plot ofV on Fig.3. where includes the contributions of the Higgs field ϕ 2 ≡ H † H non-minimally coupled to the metric and of the other Standard model fields denoted by ellipses. Quantization of this theory, which we perform in the original Jordan frame 3 , results in the effective action of the gravitating Higgs model (75) and the effective action of the CFT sector (3) Now we rewrite the Higgs effective action in the Einstein frame of fieldsĝ µν and φ according to (79) and take it on the FRW background,Γ Higgs [ĝ µν , φ ]| F RW = Γ Higgs [â,N , φ ]. Then it reads as a classical part of the action (21) but in terms of the Einstein frame scale factorâ and lapseN and with quantum corrected potential V (φ) of the above hill-like shape, The CFT quantum effective action which was obtained on the FRW background in the original Jordan frame by the conformal transformation method, Γ CF T [ g µν ]| F RW = Γ CFT [ a, N ], has the following form in terms of the original a and N Now we have to rewrite this CFT action in terms of the Einstein frame variables, Γ CFT [ a, N ] =Γ CFT [â,N , φ ], so that the full effective action -the sum of (88) and (90) will be parameterized in the Einstein frame. Under the replacement (in which we disregard logarithmic corrections in U or absorb them in the running coupling constants) the conformal time remains unchanged in view of its local conformal invariance, η = dτ N/a = dτN /â, and only the derivative of the scale factor undergoes the transition This can be interpreted as as an additional contribution of the conformal anomaly due to the transition from the Jordan frame to the Einstein frame. For a slowly varying scalar field,  (21) in terms of hatted variables. One can apply to it the above analysis. The additional restriction |ȧ|/a ≫ |φ|/ √ 6M P (we omit hats from now on) implies the bound This is stronger than the bounds (51) and (52). On account of the first bound (48), µ ≫ ω∆ or ∆ 2 /ε ≪ µ 2 /H 2 , it takes the form This suppresses ǫ in (64) even below its value in the Starobinsky or Higgs inflation models and makes the estimate for the spectral parameter n s even less sensitive to the value of ǫ.

CFT driven cosmology and the Starobinsky R 2 model
Non-minimal Higgs inflation model is very similar to the Starobinsky R 2 -model [6] from the viewpoint of inflation theory predictions [34]. Here we will show that inclusion of this R 2 -model into the full action of the CFT cosmology is very important because it not only supplies a dynamical degree of freedom with the hill like inflaton potential, but also provides the theory with a necessary finite renormalization (11) of the trace anomaly coefficient α.
Note that the contribution of the conformal anomaly to the effective actions (6) and (21) originates from the Wess-Zumino procedure of integrating the trace anomaly (10) along the orbit of the conformal group g µν = e σḡ µν . The resulting Wess-Zumino action for σ is just the difference of effective actions calculated on two members of this orbit g µν andḡ µν . It reads [35] where all barred quantities are built in terms ofḡ µν andD is the barred version of the fourth-order Paneitz This expression has an important property -with α = 0 the only higher order (quartic) derivatives of σ, contained in the combinationḡ 1/2 σDσ − 1 9 g 1/2 R 2 in the second line above, completely cancel out, and the resulting Wess-Zumino action does not acquire extra higher-derivative degrees of freedom [16]. With a nonzero α the same property holds for the renormalized action Γ R CF T [ g ] in (11). The increment of this action along the orbit of the local conformal group becomes α-independent and acquires the following minimal form The increment of this action along the orbit of the local conformal group becomes α-independent and acquires the following minimal form where again all barred quantities are built in terms of the metricḡ µν . On the FRW background spacetime of S 1 × S 3 -topology with the metric (4) and conformally related metric of the Einstein static universe ds 2 ≡ ds 2 /a 2 = dη 2 + dΩ 2 (3) , dη = N dτ /a, this difference of effective actions -the conformal anomaly contribution to (6) -equals 12 Another part of the full CFT action in (6) is the effective action of CFT fields on the Einstein universe, which consists of the free energy F (η) and the linear in η contribution of the Casimir energy E 0 . Prior to the finite renormalization (11) it equals where a particular dependence of the Casimir energy on the coefficients of the conformal anomaly was observed in the cosmological context [22] and universally derived in the class of conformally-flat spacetimes from the normalization of this energy to zero in flat spacetime [25]. 4 Since (α/384π 2 ) d 4 xḡ 1/2R2 = 3αη/16, the renormalization (11) also leads to the renormalization of the Casimir energy which acquires a particular value ∼ B/2 corresponding to α = 0. This value of E R 0 contributes to the full action of the model (6), and it is essentially responsible for the particular properties of the garland instantons [16,36]. 5 This is, of course, equivalent to the well-known statement that the coefficient of R in the trace anomaly can always be renormalized to zero by the counterterm quadratic in Ricci scalar [37], which is admissible from the viewpoint of UV renormalization due to its locality. However, we want to emphasize here that this renormalization results in twofold consequences -CFT quantum corrections preserving the non-dynamical nature of the scale factor and a particular value of the Casimir energy ∼ Bm 2 P /2 = 3β/8 which universally expresses via the topological (Gauss-Bonnet) coefficient in the conformal anomaly. In this respect, Γ CF T in the formalism of Sects.1-4 above should be everywhere labeled by R, which implies that this is the renormalized action which incorporates these properties. Below we show that this finite renormalization can be enforced by the inclusion 4 Flat spacetime can also be connected to the Einstein universe by another special conformal transformation, the relevant vacuum stress tensors being related via the conformal factor of this transformation. This gives the dependence of E 0 on β and α [25]. 5 This value of the Casimir energy guarantees that the limiting case of the garland instantons with the radiation constant C → 0 corresponds to chains of "touching" exact spheres S 4 whose contribution is suppressed to zero by infinite positive effective action [36]. This suppression excludes the infrared catastrophe of the Hartle-Hawking no-boundary state leading, as is well known, to the anti-intuitive conclusion that the quantum origin of an infinitely large universe is infinitely more probable than that of the finite one. Other values of the Casimir energy would lead to conical singularities in this set of instantons and would leave the resolution of this issue ambiguous. of the Starobinsky model R 2 -term as a part of the full action.
Indeed, the renormalization (11) can be viewed as a replacement of the Einstein-Hilbert action by the action of the Starobinsky model [6] with the coupling constant ξ of the curvature squared term This allows one to rewrite the full action of the theory with a generic set of conformal fields as a combination of the Starobinsky model with a particular value of the coupling constant ξ = α/96π 2 and the renormalized CFT action of the non-dynamical dilaton a 2 = e σ , (103) In its turn the Starobinsky action, which contains an additional conformal degree of freedom, can be rewritten in terms of the scalar-tensor theory with a nonminimal curvature coupling of the auxiliary scalar field ϕ, On the solution of the equation of motion for ϕ this scalar field expresses in terms of the scalar curvature ϕ 2 = R and, thus, recovers the original purely metric representation of the Starobinsky model. This brings us to the equivalent formulation of the full effective action as Naively here the field ϕ does not have a kinetic term, though of course it is hidden in the non-minimal coupling of ϕ to gravity. Transition to the Einstein frame according to (81) with G = 0 and U = (M 2 P + ξϕ 2 )/2 gives the relation between the Jordan frame fields and those of the Einstein frame φ andĝ µν and the Einstein frame potentialV Here the modulus of φ is chosen in order to cover the range of negative ϕ by negative values of φ.
If the slow roll condition (93) holds, then this transition leaves the CFT part of the action unchanged similarly to the case of the Higgs inflation model discussed 13 above. Therefore, on the FRW background the action (105) takes the form of our original effective action (21)(where a and φ are understood as hatted Einstein frame variables) with the minimal dynamical inflaton having the potential (107). This potential is asymptotically shift invariant at large φ, but logarithmic radiative corrections render it a hill-top shape by the mechanism discussed in the previous section. Thus, we can apply the dynamical inflaton scenario considered above.
Moreover, the inclusion of the Starobinsky model becomes indispensable if we put forward as a guiding principle a necessity to preserve the garland instantons and their dynamical inflaton generalization. This is because its R 2 term is the only means to render a nondynamical nature of the dilaton (scale factor) mode and a particular value of the Casimir energy ∼ B/2. The restriction on the CFT model for this to hold is the positivity of the overall value of α. It is necessary for the positivity of ξ in (102) -an admissible range of ξ providing unitarity -absence of ghosts in the Starobinsky model (quite paradoxically corresponding to a negative definite Euclidean R 2 -action). This restriction is rather mild, because for low conformal spins α is dominated by a positive contribution of vector particles which should only be not outnumbered by scalar bosons and fermions, where N 0 , N 1/2 and N 1 are the numbers of scalar, spinor (Weyl) and vector particles respectively. Another restriction follows from the observation that the inflaton potential plateau in (80) has over-Planckian scale 24π 2 M 4 P /α, because, in contrast to the Higgs model, ξ = α/96π 2 is small for a typical overall value of α = O(1). This difficulty can be circumvented by demanding a large value of α ≫ 24π 2 , which requires either a huge number of vector bosons or higher spin conformal particles (see below). Another possibility is to start, from the very beginning, instead of (105) with the model consisting of the large ξ Starobinsky model and the CFT theory so that even with the negative total α the sub-Planckian bound onV , would imply the bound on ξ, 4ξ + α/24π 2 ≫ 1. Then Eq.(109) can be interpreted as the way R 2 -gravity plays its double role -part of it performs finite renormalization of α to zero, Γ CF T → Γ R CF T , while the rest of it generates a dynamical inflaton feeding the CFT scenario with the potential (110). Needless to say that a similar mechanism can be attained by using a pure R 2 term in Eqs.(102)-(105) without the Einstein-Hilbert term. This would correspond to a classically scale-invariant theory in which the Einstein-Hilbert term arises by dimensional transmutation due to quantum corrections -UV renormalizable models of this type (with additional R 2 µν terms in the action) were recently considered in [38,39].

Hierarchy problem and higher spin conformal fields
The major difficulty in the construction of a realistic inflationary model via SLIH scenario is the hierarchy problem. Its inflation scale H 2 ≃ 1/2B requires the parameter B to be exceedingly large in order to belong to the sub-Planckian domain compatible with the CMB data. When expressed in terms of the coefficient β of the overall trace anomaly (10) the corresponding inflaton potential read For the Higgs inflation with a large non-minimal coupling ξ ∼ 10 4 the (Einstein frame) energy density at the start of inflation is V * ∼ 10 −11 M 4 P [8], so that for the total beta we must have β ∼ 10 13 . (112) In order to reach this value with the conventional low spin particle phenomenology characteristic of the Standard Model one would need unrealistically high numbers of conformal invariant scalar bosons N 0 , Dirac fermions N 1/2 and vector bosons N 1 in the expression for β β = 1 180 Hidden sector of so numerous low spin and very weakly interacting particles does not seem to be realistic. However, unification of interactions inspired by the ideas of string theory, holographic duality [40] and higher spin gauge theory [41] suggests that this hidden sector might contain conformal higher spin (CHS) fields [42,43,44], so that the total value of β consists of the additive sum of their partial contributions Recently there was essential progress in the theory of these fields. In particular, it was advocated in [42,43] that the values of β s can be explicitly calculated for conformal fields of an arbitrary spin s, described by totally symmetric tensors and (Dirac) spin-tensors, where ν s is their respective number of dynamical degrees of freedom -polarizations (negative for fermions). Though these fields serve now basically as a playground for holographic AdS/CFT duality issues and suffer from the problems of perturbative unitarity, which is anticipated to be restored only at the non-perturbative level, 6 it is worth trying to exploit them as a possible solution of hierarchy problem in the CFT driven cosmology. Strong motivation for this is that partial contributions of individual higher spins rapidly grow with the spin as s 6 , so that the tower of spins up to some large S generates the total value of β β ≃ 7 180 S 0 ds s 6 = S 7 180 (117) (for simplicity we consider only bosons and assume that every higher spin species is taken only once, N s = 1). Therefore, in order to provide the hierarchy bound (112) the maximal spin should be S ∼ 100, which corresponds to the following estimate of the total number of particle modes (polarizations) in the hidden sector of the theory

Conclusions
Thus, for a wide range of parameters satisfying the bounds (40) and (53) the CFT cosmology with a dynamical inflaton having a hill-like potential has instantons which are very close to the garland type solutions of the model with a primordial cosmological constant. These garland instantons are described by the approximation of two coupled oscillators and generate a new type of hill-top inflation depicted on Fig.2. Exponentially high number of the garland instanton folds m, corresponding to the upper bound on the effective cosmological constant in the range (17)-(18), guarantees smallness of slow roll parameters ǫ and η beginning with their values ǫ * and η * at the Euclidean-Lorentzian transition point. They turn out to satisfy the relations ǫ ∼ η 2 and η < 0, characteristic of the well-known non-minimal Higgs inflation or Starobinsky R 2 gravity models. Hill-like potential in these models can be generated by logarithmic loop corrections to their tree-level asymptotically shift-invariant potential, so that with the bounds (94) or (95) on the amplitude and frequency of oscillations of cosmological instantons, CFT cosmology can be regarded as a source of initial conditions for the Higgs inflation with a strong nonminimal coupling or the Starobinsky R 2 gravity. The bound (95) implies the relation ǫ ≪ η 2 at the onset of inflation, which even stronger bounds the tensor to scalar ratio in the CMB signal of these models.
It was shown that the R 2 -gravity becomes indispensable if we put forward as a guiding principle a necessity to preserve special properties of garland instantons in the CFT scenario. This is because the R 2 term is the only means to render a non-dynamical -and therefore stable non-ghost -nature of the scale factor mode and a particular value of the Casimir energy in (6) and (13). In fact it plays a double role -a part of it serves as a renormalization tool for the conformal anomaly effective action, while the rest of it contributes a dynamical inflaton feeding the CFT scenario with the asymptotically shift-invariant potential which is converted by radiative corrections to the hill-like shape.
A major difficulty with this scenario is the hierarchy problem between the subplanckian scale of inflation and the CFT scale M 4 P /β. Possible solution consists in invoking a hidden sector of conformal higher spin fields which can generate a large value of trace anomaly coefficient β (see Eq.(112)). Though this mechanism of a large β is very vulnerable to criticism regarding unitarity problem, problem of naturalness and fine tuning, it fits presently very popular ideas of string theory, higher spin gauge theory [41] and holographic duality [42], and it will be considered in much detail in the coming publication [45].
For the inverted quadratic potential (32) garland instantons establish a strong correlation between a small magnitude of the inflationary parameter |η| ≃ 0.02 and exponentially high bound (68) on the number of garland folds, m ≥ exp(π/ Bµ 2 ) ≃ exp(π √ 2/ 3|η|). This is another serious difficulty of the CFT model, which makes the admissible range of instantons exceedingly narrow (73) and too close to the upper bound of the cosmological constant range (17). In particular, this makes thermal effects of primordial CFT radiation absolutely negligible, because the thermal contribution to the CMB red tilt ∆n thermal s derived in [20] turns out to be exponentially suppressed by an exponentially high value of m. In the l'th CMB multipole it approximately reads [20] ∆n thermal whereβ is the average value of β per one effective conformal degree of freedomβ ≃ s β s N s / s ν s N s . Smallness of this contribution for large m would unfortunately destroy experimental verifiability of this model, because the thermal pattern of the power spectrumtemperature of the CMB temperature [20] -seems to be the only observable effect which distinguishes density matrix intial conditions of the Universe from its initial pure quantum state.
Fortunately, this problem can be solved by considering potentials more realistic than (32) with their convexity |V ′′ * | at the onset of inflation point being smaller than the one at their top µ 2 = |V ′′ 0 |. For such potentials correlation between small |η| and large m no longer holds, and the estimate for Bµ 2 can be raised to yield a moderate O(1) lower bound on m. This makes the thermal tilt (119) largely determined by the value ofβ following from the CFT particle content of the model.
Remarkably, for the tower of CHS fields up to S = 100, chosen above as a solution of the hierarchy problem, the estimate for the average value of β per one conformal degree of freedom amounts in view of (112) and (118) to the valueβ ∼ 10 7 (we assume that each higher spin species enters only once, N s = 1). With this value ofβ and m = O(1) the thermal correction (119) to the spectral index appears in its third decimal order, ∆n thermal s ∼ −0.001, which is within the reach of a coming technology. This means that a potential resolution of the hierarchy problem via CHS fields simultaneously makes measurable a thermal contribution to the CMB red tilt, which is complementary to the conventional tilt of the CMB spectrum [5].
The final remark regarding the CMB signal of the CFT model is that its precise power spectrum has not yet been found, except the thermal contribution (119) found in [20]. Conventional dependence of CMB parameters on ǫ and η [5] might be modified by a nontrivial speed of sound caused by nontrivial kinetic terms in the effective action (98) -the issue subject to current research [46].
Slow rate of change of W during the period of oscillation of a, τ − ≤ τ ≤ τ + , |Ẇ | ≪ ω|W |, allows us to integrate the above equation for ∆ and obtain in view of (37) Note that the quantities (16), a 2 ± ≡ (1 ± ∆)/2H 2 , via H = H(τ ) are now functions of time a ± = a ± (τ ), and a + (τ + ) is a real constant -the maximum value of the scale factor during a given period of oscillations. The above quadratic equation can be solved for the quantity in terms of the variable y running between the values ±1 during the individual period of the scale factor oscillations From this relation it follows that during this period of oscillation λ (and correspondingly ∆) change from the minimal value λ + to the maximal value λ + + 1, A fast change of ∆ means that the differentiated version of the equation (37), a 2 B = 1+∆y 1−ε , takes the form where we disregardε terms. Solving this equation with respect toȧ/a we havė If we use this relation in Eq.(38) whose right hand side in terms of y reads aṡ where ω is the same slowly varying in time frequency as (43). This equation is the generalization of Eq. (39) by the factor (1 − y/2λ(y)) 2 = 1 taking into account rapid variations of ∆. This factor which renders the oscillations of y very anharmonic does not, however, change the period of these oscillations (A.14) In the leading order approximation in ∆/ε ≪ 1 this period coincides with the period of harmonic oscillations of (44) for any positive value of λ + . In fact, with λ(y) given by (A.7) the equation for y(τ ), (A. 16) In the limit λ + → ∞ this solution reduces to the harmonic law (44), whereas for λ + → 0 (that is ∆ ≪ W ) it reads as the function not differentiable at sin ωτ 2 = 0, y = sin ωτ 2 . (A.17) The range of τ where −1 ≤ y ≤ 0 in this limit shrinks to these singular points with sin ωτ 2 = 0.
Appendix B.: "Slow roll" Euclidean solutions in CFT cosmology "Slow roll" solutions in the Euclidean regime are analogous to the Lorentzian solutions, except that the scale factor and the inflaton are oscillating instead of respectively monotonically growing and decreasing, and the scale factor oscillations are much faster than those of the inflaton field. Oscillations of φ with a constant frequency µ and rapid oscillations of a with a slowly varying frequency ω should be commensurable. Then the integer number n of full oscillations of φ within the time period T = 2πn/µ corresponds to m, m ≫ n, oscillations of a. The resulting total increase of the phase of these oscillations by 2πm is given by the integral where H 2 * = H 2 (τ * ) is the value of H at the beginning of the full period τ * = 0 -the turning point of the solution, one finds leading corrections in ∆ and ∆ φ to time dependent H 2 and ε = 1 − 2BH 2 , sin 2 (µ + ω * 2 )τ 2µ + ω * + sin 2 ω * τ T i which according to Appendix A is equal to 2π/ω i with As a result the leading corrections to the frequency commensurability equation (54) and η have the following form where all starred quantities are taken at τ * . Suppression of these corrections by the factor ∆∆ 2 φ /M 2 P is obvious from the dynamical equation (31) for H 2 , while the extra suppressing factor n 2 /m 2 ≪ 1 comes from integration of slow φ-harmonics modulated by rapid oscillations of the scale factor. As we will see below, these corrections of the magnitude ∆ ε × n 2 m 2 (B.24) will be negligibly small, so that they can be discarded in the leading order, and moreover the star labels of H, ε and ∆ can also be omitted in this approximation. Therefore, as in Sect. 3 where the first inequality follows from (B.29). One can show that the case of R ≤ 2/Bµ 2 m 2 leads to the contradiction, because the resulting value of R becomes too high, R ≥ 1 − 2∆ 2 , and violates the low temperature limit R(η) → 0 at η → ∞. On the contrary, the case of does not lead to new bounds except the positivity of ∆.
The further analyzes relies on the low temperature limit of the free energy F (η) for large η, F (η) ≃ −N 0 e −ω0η , where N 0 is the number of lowest spin conformal fields having the lowest value ω 0 of the energy of field-theoretical oscillator on S 3 (ω 0 = 1 for a spin zero conformal field). In this limit the function R(η) becomes The above bounds are consistent with the admissible range (53) of µ vs ω, ω ≫ µ ≫ ω∆, which can be rewritten as 1 ≫ µ 2 ω 2 ≃ n 2 m 2 ≫ ∆ 2 , (B.41)