A minimal scale invariant axion solution to the strong CP-problem

We present a scale invariant extension of the Standard model allowing for the Kim-Shifman-Vainstein-Zakharov (KSVZ) axion solution of the strong CP problem in QCD. We add the minimal number of new particles and show that the Peccei-Quinn scalar might be identified with the complex dilaton field. Scale invariance, together with the Peccei-Quinn symmetry, is broken spontaneously near the Planck scale before inflation, which is driven by the Standard Model Higgs field. We present a set of general conditions which makes this scenario viable and an explicit example of an effective theory possessing spontaneous breaking of scale invariance. We show that this description works both for inflation and low-energy physics in the electroweak vacuum. This scenario can provide a self-consistent inflationary stage and, at the same time, successfully avoid the cosmological bounds on the axion. Our general predictions are the existence of coloured TeV mass fermion and the QCD axion. The latter has all the properties of the KSVZ axion but does not contribute to dark matter. This axion can be searched via it's mixing to a photon in an external magnetic field.

no new physics between Higgs vev and Planck scale [17][18][19]. However, the renormalizability of the theory is lost. This should not be thought as a serious drawback since the theory is non-renormalizable in the presence of gravity. Nevertheless, the model may be considered as a valid effective field theory with a cutoff scale 1 close to the Planck scale [20,21].
The SMASH model mentioned above includes an explicit intermediate mass scale (the axion decay constant). That provides with a sufficient level of fine tuning needed to obtain both the small Higgs mass and viable inflation (see the discussion in [9]). In this paper, we address a question whether it is possible to solve this hierarchy problem exploiting the idea of spontaneously broken scale invariance in the absence of heavy particles. In order to make SMASH scale invariant, one needs to add the dilaton field. However, the complex scalar field, responsible for the Peccei-Quinn symmetry breaking, is already present in the model.
The novel idea of this paper is to construct a model in which the Peccei-Quinn scalar also plays the role of the dilaton. In this setup, the QCD axion corresponds to a phase of the complex dilaton field. 2 Thus, following the idea of minimality, we do not add extra degrees of freedom, as compared to SMASH. At the same time, the scale invariance of the model might help to resolve the hierarchy problem. The purpose of this paper is to construct a model, stable with respect to the quantum corrections, which solves the strong CP problem and provides realistic cosmology and particle physics.
In the proposed framework, the scale symmetry must be broken together with the Peccei-Quinn symmetry. It had happened in the early Universe before inflation at the energies around the Planck scale. This means that the Peccei-Quinn symmetry has never been restored thermally. In order to avoid overproduction of the axion dark matter without finetuning the initial θ-angle, we require that during inflation QCD was in the strongly coupled phase. A way to obtain that is (following the idea of Ref. [27] 3 ) to let the QCD coupling constant to depend on the values of Higgs and dilaton fields. This would allow for normal QCD phenomenology in the electroweak vacuum, while at large values of the Higgs field (corresponding to inflation), the effective QCD coupling may well be larger than unity. In this regime, the axion can acquire a large mass. Since the fluctuations of the heavy (compared to the Hubble scale) axion during inflation are negligible, it can be left in the minimum of its potential. Therefore, axion dark matter is not produced in this case, avoiding this way all cosmological constraints.
Let us briefly describe the cosmological scenario, corresponding to the models under consideration. The scale invariance and Peccei-Quinn symmetry were broken by the dilaton vev before inflation driven by the Higgs field 4 . During inflation, QCD was strongly coupled and the axion was heavy, dynamically providing with zero θ-angle. Just after inflation, QCD becomes weakly coupled, with the same constant as it should be in the SM. The axion looses its mass and remains massless until the QCD phase transition. After reheating, the Universe evolves similarly to the νMSM scenario [2][3][4]. Namely, dark matter is provided by one of the sterile neutrinos. Two remaining sterile neutrinos generate the lepton asymmetry (which transfers to the baryon asymmetry) and SM neutrino masses. In this paper, we discuss several conditions which allow to obtain this scenario and satisfy all cosmological and laboratory constraints.
The paper is organized as follows. In Section 2 we discuss the problems arising in the simplest model for the Higgs inflation in the presence of a complex dilaton field. Then we consider the general scale invariant lagrangian with the complex dilaton being simultaneously the Peccei-Quinn scalar. In Section 3 we discuss several constraints on the action to fulfil all the phenomenological requirements consistently. We give an explicit example of a viable model which is stable against quantum corrections during the inflationary stage. In Section 4 we consider the general predictions for the scale invariant axion scenarios pointing out the distinguishing features for the experiments and observations. Section 5 is devoted to their distinctive observational signatures.

Axion as a phase of dilaton
The main idea behind this paper is to identify the complex Peccei-Quinn scalar with the dilaton field. In order to construct a model which is stable with respect to quantum corrections and, at the same time, can provide viable cosmology, solution for the strong CP problem and SM physics at low energy, several conditions has to be satisfied. First, QCD axion has to be massive during inflation and its mass must be larger than the Hubble scale. As it was mentioned, this can be reached if the effective QCD coupling is large in the early Universe. Second, the inflaton potential must provide the spectrum of scalar perturbations which is in agreement with the recent Planck data [33]. The next requirement concerns the extra fermion charged under SU(3) c . It must be heavy during inflation, in order not to suppress the axion mass [36], while in the SM vacuum its mass should be less than 10 TeV, to avoid large contributions to the Higgs mass [19]. Finally, since we want the model to be a predictive effective theory both during inflation and for low energy physics, we require that the (background-dependent) cutoff scale always must be much larger than all relevant scales.
In this section, we discuss the simplest extension of the Higgs dilaton model providing with a solution for the strong CP problem which is viable from the point of view of the low energy physics. However, this model cannot describe the inflationary stage. Therefore, we consider a general scale invariant lagrangian which may also allow for the viable cosmology.

The simplest model
Scale invariant model with the QCD axion. Let us start with the simplest model allowing for the simultaneous breaking of scale and Peccei-Quinn symmetry by the vev of the complex dilaton field charged under the Peccei-Quinn symmetry. We propose the scalar and gravitational sectors of the model to be of the same form as in the Higgs-dilaton model [26], except that the dilaton X is a complex field. Thus, the lagrangian looks where h stands for the radial component of the SM Higgs field in the unitary gauge and the Higgs vacuum expectation value (vev) is provided by the dilaton vev: α |X| = v. ξ, ξ h are dimensionless couplings of the fields to gravity. The term β|X| 4 actually provides a cosmological constant [13]. Therefore, β must be tiny, so we omit this term hereafter, since, in this paper, we do not discuss the dark energy problem. To implement the KSVZ axion scenario with the Peccei-Quinn scalar X, one has to add two extra Weyl fermions Q,Q in the fundamental and antifundamental representations of SU(3) c and with the SM hypercharges −1/3, 1/3, respectively. Their Yukawa couplings to the field X are The presence of the hypercharge for Q allows for their mixing with the SM quarks, making it possible for them to decay, avoiding the problem of their overabundance [31,32]. Furthermore, one needs an additional U(1) P Q symmetry which acts as a chiral rotation on Q and as a phase rotation of X. In this setup, the phase of the field X would be connected to the axion field, where F a is defined by the canonical normalization of the axion field. The field a acquires mass due to non-perturbative QCD effects and thus provides with a dynamical relaxation of the QCD θ-angle to zero value. This model is a viable solution to the strong CP problem from the point of view of the particle physics. It provides with the invisible axion which has suppressed couplings to all SM particles. Now we consider this particle in a cosmological context and discuss several problems arising in this way.
Cosmology of the axion.
The scalar sector of the described model coincides with the Higgs-dilaton model [26]. The difference is that the dilaton field is substituted by the modulus of the Peccei-Quinn field X. Thus, the inflationary stage can be obtained in the same way as in the original model. The values of ξ and ξ h are constrained by the tilt and amplitude of scalar perturbations, ξ < 0.01, ξ h ∼ 10 3 . The phase of field X (axion) do not affect the inflationary dynamics, however, it is be displaced from the QCD minimum. This misalignment can lead to the overproduction of the axion dark matter during QCD phase transition.
The standard cosmological scenario for the QCD axion dark matter involves the thermal restoration of the PQ-symmetry after inflation. This can not be done in our case because it would mean the restoration of scale invariance as well. However, the discussed scale invariant models can be the consistent effective description only in the broken phase, because they are defined only for large values of the dilaton field. Therefore, we deal with the scenarios without restoration of the PQ symmetry. The common problem of such models is that, in order to obtain the correct amount of axion dark matter, one needs to tune the initial value of the QCD θ-angle to be at most 10 −3 [22], though its natural value is of order unity. Another problem arising in large scale inflation is an amplitude of axion isocurvature perturbations which is incompatible with the Planck data [23]. To avoid both problems, we propose here that during inflation the axion was massive. This can be achieved if the QCD sector is modified during inflation in such a way that the effective QCD coupling constant is large at this stage. We write the non-minimal coupling between the Higgs, dilaton and gluons in the form, Here the scale invariant function f is such that for large Higgs field values during inflation it makes the coupling constant to be large. At the same time, we require f (0) = 1, such that we end up with 'standard' QCD in the electroweak vacuum.
In the strong coupling regime, the PQ symmetry (which is the axion shift symmetry) is broken nonperturbatively leading to the appearance of the axion potential. In the SM vacuum (at zero temperature) the potential reads, where Λ QCD is the characteristic scale at which the model enters the strong coupling regime. The value of Λ 0 QCD can be connected to the QCD coupling constant g 0 at some reference energy scale (here we take the Planck scale M P ). In the SM, the scale Λ 0 QCD at which the running of g 0 brings it to the value g 0 (Λ 0 QCD ) ∼ 1 can be found from renormalization group equation [24,25], In the inflationary background, the value of Λ QCD is modified compared to its usual value Λ 0 QCD . Namely, this scale is defined by the effective QCD coupling g 2 ef f = g 2 0 /f (h/|X|), taken at the reference scale M P . Solving the similar renormalization group equation one obtains, Therefore, the axion mass is From the lagrangian (1) one can obtain the effective axion field normalization during Higgs-dilaton inflation to be Considering Higgs-dilaton inflation, in order to obtain the tilt of scalar perturbations to be in agreement with the Planck data, the value of ξ should be smaller than 10 −2 [26]. Thus, to satisfy the condition m a > H, we need This consideration is valid only in the absence of charged under SU(3) c particles lighter than Λ QCD , otherwise, the axion mass gets suppressed [36]. Since the SM fermion masses are given by the large value of the Higgs field, they are heavy enough during inflation. However, we need to require that the mass of the additional quark Q is larger than Λ QCD , to obtain no suppression of the axion mass. In a simple model (2), this would provide an unnaturally large value o f the coupling y. The absence of the large impact to the Higgs mass requires m Q < 10 TeV [19] at low energies. The described tension is not possible to relax if the value of y is constant. However, we can assume that the value of y depends on the relation h/|X| in such a way that provides with the heavy quark during inflation and light quark in the low energy domain. Clearly, this idea calls for more general study which includes other scale invariant modifications of the lagrangian (1). In the next sections, we construct models which allow for obtaining viable inflation with a sufficiently massive axion.

The general scale invariant model
The most general scale invariant lagrangian for the Higgs field and complex dilaton field coupled to gravity with at most first derivatives reads Here the functions f 1,2,3 , V 1 and the parameter α are defined in such a way that, after the scale symmetry gets spontaneously broken, we recover the SM Higgs potential at low energies. The extra fermions in the model are represented by three sterile neutrinos and the fermion Q charged under SU(3) c introduced in the previous section. Thus, the lagrangian for the fermionic sector of the model can be written as Here L f,νM SM stands for the fermionic part of the νMSM lagrangian in which the masses of the sterile neutrinos are given by the dilaton vev X 5 and y(h/X) is the Yukava coupling which might depend on the values of the Higgs and dilaton fields. As it was discussed earlier, we also modify the QCD lagrangian as in (4), while the other sectors of the SM remain unaltered.

Einstein frame consideration
The lagrangian (11) is written in the explicitly scale invariant form. However, the presence of non-minimal couplings of scalar fields to gravity makes it not convenient to study the field dynamics during inflation in these variables. This lagrangian can be rewritten in the Einstein frame, in which the action for the gravity takes its canonical form and all scalars are minimally coupled. This can be done performing the Weyl transformation, and defining the new variables, we find that the scalar and gravitational sectors in the Einstein frame, read Here, the functions A, B, C, F a are related to the functions in the lagrangian (11), as The kinetic term of (16) can be diagonalized by the change of variables A straightforward calculation reveals that We can see that the two fields,r (hereafter we will refer to it as Einstein frame dilaton -EF dilaton) and the axion a, do not have a (classical) potential. This means that they possess shift symmetry. Therefore, the only field which can be responsible for inflation is θ, whose kinetic term reads Note that the field θ can be canonically normalized, something that is not possible for the EF dilaton and axion fields, since their kinetic terms depend on θ.

Constraints on the action
Let us list here several requirements that the lagrangian has to satisfy in order to obtain a viable scenario.
1. QCD has to be in the strong-coupling phase during inflation which constrains the function f (θ) in (4).
2. The axion mass during inflation is constrained from the both sides. The condition m > H provides with zero initial θ-angle and guarantees the absence of isocurvature perturbations. But the axion can not be too heavy, in order to obtain no overproduction of the dark radiation (see Appendix B).
3. We restrict ourselves to inflationary potentials which have only one minimum corresponding to the SM vacuum. Hence we require that at large values of θ, the potential converges to a constant of order ∆ 2 M 2 P with ∆ = 10 −5 , to be able to provide the scalar perturbations in agreement with the Planck data [33]. 4. To obtain the correct values for all the inflationary parameters, we require the following asymptotic form for the inflaton kinetic term at large θ, Although this is not unique choice 6 , we will use it since it is the simplest one. Note that this form corresponds to a class of α-attractor models [34] which are known to give predictions in a good agreement with the Planck data [33].

5.
Addressing the question about the quantum stability of the model during inflation, we have to require that the effective suppression scale for all non-renormalizable operators is much larger than the Hubble scale. For the class of models under discussion, we show in Appendix A that this cutoff is nothing else than Λ. Thus we need Λ ≫ H 7 .
6. The SU(3) c charged fermion Q must be heavy (m Q > H) during inflation, otherwise the axion mass is suppressed. On the other hand, a naturally small Higgs mass requires that at low energies its Yukawa coupling y is small enough. This, in turn, corresponds to the mass around 1 − 10 TeV. This allows to put constraints on the function y(θ).

General conditions for the functions in the lagrangian
We start with discussing the conditions that allow for a viable slow roll inflation with a heavy axion. From (17), we obtain For simplicity, let us assume that F a (θ) becomes a constant at large θ and F a < H. In this case, we have As for the inflaton kinetic term (20), we can simplify its asymptotic form for large θ within (23) and assuming γ/ξ ≪ 1. Besides that, let us also assume that the function f 2 grows not too fast. Namely, we require f ′ 2 (θ)/f 2 (θ) cf 2 /θ where c is a constant of order unity. 8 6 For example, in the usual Higgs dilaton model, the asymptotic form of the inflaton kinetic term reads (in our notifications) K θ ∝ 1/(θ 2 (θ 2 + ζ)) [35]. We leave the study of this case beyond the scope of this paper since it contains more free parameters than the case which we consider here. 7 Note that the here we discuss the cutoff scale calculated in the inflationary background. The cutoff scale in the vacuum can be smaller than the Hubble scale. In this case, one needs additional assumptions about the UV completion of the model. In the discussed constructions, we postulate the asymptotic shift symmetry for the inflaton field in the Einstein frame for large field values. This assumption protects the flatness of the inflaton potential from possible power law corrections. 8 The equality holds for polynomial functions.
Thus, the kinetic term for θ can be simplified an is given by To start with, let us consider the case in which the terms with f 2 dominate both in the numerator and denominator. Under this assumption, we obtain One can see that this model fails to describe inflation consistently, since its cutoff much smaller than the Planck mass. However, if the terms with f 3 dominate, then we get consistent inflation, since To conclude, the conditions for f 1 , f 2 , f 3 , V 1 at large θ leading to the consistent inflationary stage in the presence of a massive axion, read Although this class of solutions might not be unique, it is the one providing with a viable inflationary stage. Now let us turn to the modification of the kinetic term for the gluons (4). The function f has to make QCD strongly coupled during inflation, i.e. at large θ. The axion mass provided by the effective strong coupling scale Λ QCD , reads The axion mass must satisfy the inequality H < m < 10 15 GeV .
The upper limit corresponds to the production of the dark radiation, due to the disappearance of the axion mass after inflation. In the next section, this effect will be considered in details.
We focus on the case where the strongly-coupled phase of QCD finishes just after inflation, in order to avoid multiple transitions when crossing the critical value of θ before reheating. In the next section, we discuss a particular example for the various functions.
Finally, the mass term for the new heavy quark Q E in the Einstein frame is given by Here the fermion gets rescaled in order to make its kinetic term canonical. To have a quark with mass 1 TeV in the SM vacuum, we should set On the other hand, to obtain a heavy quark during inflation, we need

An explicit example
In this subsection we discuss the set of functions of the simplest form (that is, polynomials of the lowest degree) which satisfy all aforementioned conditions. The simplest choice of f 2 quadratic in θ, which corresponds to the model discussed in Sec. 2.1, does not work.
Therefore, we start with the set For this choice, the kinetic term for the inflaton (in the limit γ ≪ ξ) becomes The required potential correction V 1 (θ) is Using the above, we find that the axion constant becomes F a = M P γ/ξ. The inflaton potential can be expanded as The change of variables θ = e φ/Λ / √ µ, with makes it possible to canonically normalize the inflaton φ. The potential in terms of φ, reads Note that the inflation finishes around θ 2 = θ 2 e ∼ aξ/ξ h . In order to neglect the next terms in the 1/θ expansions for all functions in the lagrangian until the end of inflation, we must require θ e ≫ 1. To have the cutoff Λ ∼ M P , it would be enough to set d ≫ aξ. The next requirement is to verify that the next terms in the expansion of the kinetic term K θ can indeed be neglected. After a straightforward calculation, one obtains the extra condition These conditions do not contradict to each other, so they can be easily satisfied simultaneously.
Turning to the QCD part of the lagrangian, in order to fulfil the requirements of the previous subsection, we can set 9 This function interpolates between 1 for small θ and b for the inflationary region. Here we require that the QCD transition happens only once and right after inflation: θ tr ∼ θ e . The value of b is defined by the desired value of the axion mass m a : For example, if we consider γ = 10 −4 ξ than b lies in the range between 0.14 and 0.16, in order to provide the axion mass H < m a < 10 15 GeV. To satisfy the conditions for the Q quark, we need Note that large values of δ would not bring the theory to the strong coupling regime because there is an additional suppression of the operators by the cutoff scale. We discuss this in more detail in Appendix A. Summarising, we can set for example These values allow for consistent inflation, and also for an axion and quark Q which are heavy enough during inflation and light at low energies. All additional interactions, in the low energy limit, are suppressed by the scale of the dilaton vev. In Appendix B we prove that the outlined choice of the model parameters is stable with respect to radiative corrections during inflation as well as in the low energy domain. Although it is hard to describe the intermediate stage of preheating in this model, we do not expect any dangerous effects. The rigorous study of reheating is beyond the scope of this paper, we write here some comments on this issue. In order to leave the model in a predictive domain after inflation, the reheating temperature must not be high enough to restore the scale invariance: If this condition is satisfied (which should be the case for ξ 10 3 since we expect T reh ∼ 10 13 GeV, as in the Higgs-dilaton model [26]), then reheating finally provides the thermalized SM plasma with some admixture of coloured quarks Q. The latter are expected to decay via their mixing with the SM quarks. However, we should care about the production of light particles whose coupling to SM is suppressed because it may result in undesired relics in the Universe. We address this issue in the next section.

Phenomenology and predictions
We start with the discussion of the production of the Goldstone bosons present in the model, i. e. axion and dilaton. Once produced, they remain in the expanding Universe. Their energy density rescales as radiation. These particles do not decay because they have strongly suppressed interactions with the SM particles. For the dilaton production, we expect a significant suppression both for the perturbative and non-perturbative channels, similarly to the conclusions of [37]. In this model, we expect the reheating temperature to be comparable to the one in Higgs inflation, T ∼ 10 13 GeV [38], or even higher. Due to that, the inflaton has no time to decay into dilatons via the Planck suppressed couplings.
However, the latter is not true for the axion because, during the fast QCD transition after inflation, it's mass m a quickly changes from at least 10 13 GeV to zero. This may lead to an efficient production of relativistic axions. In Appendix B, we consider a simple analytical model providing an estimate for the energy density of axions created due to the non-adiabatic vanishing of the mass. If the duration of the transition, τ = 1/(βm a ), is smaller than the Hubble time and the period of inflaton oscillations, the axion density is given by We see that although the actual value of β (or τ ) is hard to extract from the model, the effect depends on this parameter only logarithmically. Thus, it seems reasonable to have log β 10. Hereafter we use ρ ∼ m 4 a as an order of magnitudes estimate. Let us compute the amount of axion dark radiation arising from the discussed effect. At reheating, the energy density reads Here we used the fact that the preheating stage was matter dominated, as in Higgs inflation [38]. The transition happened at the moment corresponding to H 0 ∼ 0.1H and scale factor a 0 . The index 'reh' here is related to the moment of reheating. The amount of axion dark radiation measured in the number of extra neutrino species is [37] ∆N ≈ 3ρ reh ρ SM , ρ SM = π 2 30 gT 4 reh .
Plugging in numbers, we obtain ∆N ∼ 10 −2 m 0 10 15 GeV We see that the amount of the dark radiation crucially depends on the axion mass. It might be significant and available for future CMB observations [39,40]. Given the present Planck bound [41], we are not able to increase m a further than 10 15 GeV. Thus, the axion mass lies between 10 13 and 10 15 GeV.
The discussed model provides a nearly invisible axion at low energies. In this domain, its effective decay constant is F a = M P / √ ξ, which can well be of the order of Planck scale or lower. The axion does not contribute to dark matter anymore. However, as we have seen, it can significantly contribute to dark radiation for some parameter choices. Dark matter can be provided by one of the sterile neutrinos, exactly as in the νMSM scenario. The other sterile neutrinos can be responsible for the lepton and baryon asymmetry as well as for the SM neutrino masses via the see-saw mechanism.
The possible experimental window for this kind of models may be connected to collider searches of strongly interacting quarks. The LHC limit on new quark masses is m Q 1 TeV [42], which is quite close to the naturalness bound m Q 10 TeV [19]. So the relatively light extra quarks may be probed at LHC and the remaining preferable window for their mass is narrow.
If the decay constant of the axion is significantly smaller than the Planck scale (which can be the case), then the axion can be probed via its mixing to the photon in a magnetic field.
All astrophysical bounds dealing with the axion-photon conversion (cooling of stars [43]) are relevant for the discussed model as well. Future experiments searching for solar axion emission [44,45], and for the light shining through the wall [46], are expected to put stronger constraints on the QCD axion parameters even if it does not contribute to the dark matter.
In this case, all experiments specially designed for searches of the dark matter axions would give null results.

Conclusions
In this paper, we examined the possibility to provide a solution to the strong CP problem in QCD in the framework of spontaneously broken scale invariance. This minimal extension of the Standard Model calls for the complex scalar field (SM singlet) and two Weyl fermions charged under SU(3) c and U(1) Y SM groups. The modulus of the scalar field plays the role of the dilaton, while its phase is the QCD axion. Thus, the model (in the SM vacuum) contains two light scalar fields-dilaton and axion-being the Nambu-Goldstone bosons of the broken scale invariance and Peccei-Quinn symmetry, respectively. We show that it is possible to construct a model satisfying all bounds for inflation and QCD axion by the price of considering non-linear lagrangians for the Higgs, gravity and QCD sectors. We prove that, under several conditions for the functions in the lagrangian, one can obtain an effective theory expanded over the classical inflationary background, as well as over the SM vacuum, that is stable with respect to the quantum corrections. Given that, we can conclude that the presented scale invariant models for the QCD axion might be a technically natural, self-consistent effective description for low energy physics as well as for inflation.
The distinct "signature" of the discussed scenarios is the presence of a QCD axion whose mass, decay constant and coupling to photons are connected as in the usual KSVZ model [7,8]. Unlike the usual case which is widely discussed in the literature, this particle does not contribute to dark matter, providing null results for the axion dark matter searches.
However, the axion mixes with the photon in the presence of a magnetic field, potentially allowing for its detection in laboratory experiments and astrophysical observations. The author is indebted to M. Shaposhnikov for bringing the attention to the problem and for the helpful discussions. A.T. thanks also G. Karananas and S. Troitsky for useful conversations. This work was partially supported by the Swiss National Science Foundation. The part of work related to the study of inflation and cosmological consequences of the discussed model was supported by Russian Science Foundation grant 14-22-00161.
If ζM l /Λ l < 1 one can obtain the lowestΛ from the term with n = l + 1, Therefore, we can take the cutoff scale of the model to beΛ ≃ Λ ∼ M P , given the assumption ζM l /Λ l < 1 is satisfied (which is usually the case). Let us check the validity of this conclusion for the mass term of the extra quark Q since we need the large dimensionless constant δ ∼ 10 10 (43). The corresponding operator looks, For large θ it can be expanded as We see that the cutoff scale is not lower than Λ ∼ M P within the discussed choice of parameters. Analogically, one can find the cutoff scale for other non-minimal operators in the lagrangian, including G µν G µν and the inflaton potential, to be Λ during inflation. As for the SM vacuum, all the additional interactions in the model are suppressed by the dilaton vacuum expectation value providing the cutoff scale to be of order M P / √ ξ. Therefore, the physics at low energies near the SM vacuum is the same as in the Standard Model.

B Radiative corrections to the model parameters
Let us first discuss the loop corrections that might affect the coupling constants in the low energy domain. All operators suppressed by the large dilaton vev (for example, corrections to the Higgs, dilaton and QCD kinetic terms) actually provide with negligible impact to the running of the dimensionless couplings. This fact, typical for the non-renormalizable models with large cutoff scale, was confirmed for the Higgs inflation [50] and Higgs-dilaton model [21]. Thus, we do not expect large loop contributions to the dimensionless couplings from the nonrenormalizable operators. From the other side, the small values of these couplings guarantee the stability of the choice of constants (43). The operators suppressed by the dilaton vev can receive corrections from renormalizable operators. However, these corrections are small provided that all the couplings are weak.
The separate discussion is needed for the dimensionless couplings ξ and ξ h . The large value ξ ∼ 10 3 of non-minimal coupling between the dilaton and gravity is stable because the dilaton is very weakly coupled to the SM sector. Similarly to the case of the Higgs inflation, this large value does not lead to significant contributions to other constants, since it translates to the SM sector only through the gravitons. Unlike the dilaton coupling ξ, the value of ξ h cannot be too small for the naturalness reasons. Namely, at one loop the SM couplings provide the running of ξ h [47,48] 10 , Here y t is top quark Yukawa coupling and g 1,2 are the SM gauge couplings. As it follows from this formula, the value ξ h 0.01 is stable under the loop corrections.
Let us now turn to the model expanded over the inflationary background. Technically, it is much simpler to consider the loop corrections in the Einstein frame, in terms of the diagonalized fields. Since the cutoff scale is larger than all the relevant mass parameters of the model (masses of the particles, effective Λ QCD scale), we expect all contributions to be suppressed by this cutoff scale. To illustrate this fact we consider the corrections to the kinetic term of gluons. For large θ ≫ θ e one can expand the QCD part of the lagrangian, Here we expand the field φ over its background valueφ. Before the end of inflation, the last term in the expansion of the function f is small due to the presence of the exponent.
One can see that all non-renormalizable interactions are suppressed by the scale Λ. Their contribution to the value of b is negligible compared to the usual running of the effective QCD coupling g 2 ef f = g 2 0 /b which is already accounted for in the rescaling of the Λ QCD (7). Therefore, the radiative corrections leave the value of b in the region 0.14 < b < 0.16 which is preferred for the axion phenomenology.
The similar expansions over the background inflaton field φ can be obtained for other non-renormalizable operators. These expansions contain the similar exponential factors e −φ/Λ which provide with an additional suppression of corrections to the bare parameter values that come from the interactions suppressed by the cutoff scale. Thus, the inflationary observables that depend on the inflaton potential are not significantly affected by the radiative corrections. This fact has been already established for the Higgs inflation [50], as well as for the Higgs-dilaton model [21] in which the inflaton has the exponentially flat potential. C Axion production due to the disappearance of mass Let us approximate the axion mass after inflation to be vanishing faster than all cosmologically relevant times (such as the Hubble time and period of the inflaton oscillations). In order to get an analytic solution we use the following mass dependence, m 2 (t) = m 2 0 (1 − th(βm 0 (t − t 0 )), β ≫ 1.
For such a function a Bogolubov coefficient and the energy density of axions can be found explicitly [49], Note the weak logarithmic dependence of the resulting energy density on the rate of mass vanishing. The described effect of axion production seems to be the main impact of the energy density of dark radiation. Other possible sources, such as the Universe expansion and the decay of the inflaton oscillations, can be neglected since they should provide with the result similar to the case of the dilaton. The latter is known to be negligible [37].