Does the detection of primordial gravitational waves exclude low energy inflation?

We show that a detectable tensor-to-scalar ratio $(r\ge 10^{-3})$ on the CMB scale can be generated even during extremely low energy inflation which saturates the BBN bound $\rho_{\rm inf}\approx (30 {\rm MeV})^4$. The source of the gravitational waves is not quantum fluctuations of graviton but those of $SU(2)$ gauge fields, energetically supported by coupled axion fields. The curvature perturbation, the backreaction effect and the validity of perturbative treatment are carefully checked. Our result indicates that measuring $r$ alone does not immediately fix the inflationary energy scale.


I. INTRODUCTION
The inflationary paradigm has been successful over the past few decades to serve as a mechanism to produce the observed inhomogeneities in the universe such as the cosmic microwave background (CMB) anisotropies and large-scale structure (LSS), while resolving the conceptual difficulties in the hot big bang scenario. An important prediction in the framework is generation of the B-mode polarization in the CMB [1], whose signal is conventionally quantified by the tensor-to-scalar ratio r ≡ P h /P ζ | k=kCMB . The current bound is r < 0.07 at k CMB = 0.05 Mpc −1 with 95% confidence [2], and a number of proposed missions are expected to improve the bound to O(10 −3 ) (see e.g. [3]). The conventional relationship between the tensor-to-scalar ratio and the Hubble parameter during inflation is where H inf is the Hubble parameter during inflation and P ζ ≈ 2.2 × 10 −9 has been used [4]. An immediate implication of (1) is that detection of r would fix the inflationary scale at such high energy levels as beyond our current experimental reach. Considering the ongoing and upcoming experimental efforts for B-mode detection, it is right time to test the validity of the conventional prediction (1). In general, the value of r at cosmological scales can be estimated as the spectrum of the energy fraction of gravitational wave (GW) at the horizon crossing divided by P ζ r P −1 where ρ inf ≡ 3M 2 Pl H 2 inf and dρ GW /d ln k H 2 M 2 Pl P h at the horizon crossing. The energy density of GW from the vacuum fluctuations produced during the quasi de Sitter expansion must be characterized by the Hubble scale dρ vac GW /d ln k H 4 inf , leading to the conventional relation r vac ∝ H 2 inf . On the other hand, if GW is induced by another energy source, the conventional relation (1) may be altered. Provided that an energy source ρ s generates GWs with efficiency γ, one generally expects which can be significant even if ρ s ρ inf and γ 1 thanks to the smallness of P ζ . Conventionally, however, an efficient energy transfer from a source to GW has been assumed to be rather difficult. The reasoning is rooted in the decomposition theorem in cosmology, which states that perturbations around a homogeneous and isotropic background can be decomposed into scalar, vector and tensor sectors that are mutually decoupled at the linearized order. Since GW is the only tensor degree of freedom in the Einstein gravity, we have no choice but use the source term from scalar δS or vector perturbation δV i which is schematically written as ijk are operators traceless and transverse in the indices ij that depend on time and spatial derivatives. However, the decomposition theorem bans the existence of such operators at the linear order. Although the second order effects (e.g. ∂ i δS∂ j δS, δV i δV j ) arXiv:1705.01533v1 [astro-ph.CO] 3 May 2017 are allowed to generate GW, the efficiency of the energy transfer is suppressed, because the coefficients of the source term effectively becomes the order of perturbation, O There is a loophole in this argument. If O (V ) ijk in (4) consists of the background vector fieldV i (t), GW can be sourced at linear order byV i δV j . It is known that SU (2) gauge fields can achieve this without disrupting background isotropy by taking a particular configuration. 1 Moreover this isotropic configuration is realized as an attractor solution, if SU (2) gauge fields are coupled to a rolling pseudo-scalar field [6]. Therefore SU (2) gauge fields can source the GW through the termsV i δV j without violating the isotropy of the universe at the linear order, thus with a high efficiency of the energy transfer.
As we shall see later, the energy source ρ s to generate GW is the (linear) perturbation of an SU (2) gauge field. It is produced as quantum fluctuations and thus acquires the amplitude O(H inf ) around the horizon crossing. In addition, however, it experiences a transient instability around horizon crossing and is amplified by an exponential factor. As a result, the energy fraction of the source and the efficiency factor of energy transfer in (3) are given by where s now denotes the perturbation of SU (2) gauge field, ρ A is its background energy density, and m Q is the SU (2) mass parameter in the units of H inf . For values of m Q with H inf √ Ω A e 2m Q O(10 12 ) GeV, one can realize a detectable r even in the case of low-energy inflation.

II. SPECTATOR AXION-SU (2) MODEL
In our consideration of GW production, we leave the gravity sector as the standard Einstein-Hilbert and the inflation model unspecified, which is also responsible for generating the observed curvature perturbation. We then consider the axion-SU (2) sector with the action [7] (see also [8]): where χ is a pseudo-scalar field (axion) with a cosine-type potential are the field strength of SU (2) gauge field and its dual, respectively, and λ is a dimensionless coupling constant. At the background level, it is shown that the isotropic configuration of the SU (2) gauge fields, A a 0 = 0 and A a i = δ a i a(t)Q(t), is an attractor solution while the vev of χ(t) slowly rolls down its potential [6,7]. At the perturbation level, δA a µ contains two scalar δQ, M , two vector M i and two tensor t ij polarizations as dynamical degrees of freedom [6,7]. Interestingly, t ij is coupled to the metric tensor modes h ij already at the linear order, and only one circular polarization mode of t ij is substantially amplified due to a transient instability around the horizon crossing. It then efficiently sources one polarization of GW h ij at the linear order, if m Q ≡ gQ/H > √ 2 [9]. Therefore we focus on t ij among the perturbations of A a µ . The Einstein equation at the background yields Pl H 2 , and dot denotes the cosmic time derivative. The inflaton part ρ φ and φ ≡ −ρ φ /6M 2 Pl H 3 depend on the inflation model, and ρ t and t ≡ −ρ t /6M 2 Pl H 2 denote the contributions from the perturbation t ij on the background dynamics, which will be discussed later. The equations of motion for χ(t) and Q(t) arë where we include the backreaction terms, T Q BR and T χ BR , from t ij . Without the backreaction, one can show that the effective potential of Q uplifted by the coupling to χ acquires a non-zero minimum at Q min µ 4 sin(χ/f )/3gλH 1/3 , if χ slowly rolls and the coupling is sufficiently strong [6,7]. The tensor perturbations consist of t ij and h ij , and each of them can be decomposed into the circular polarization modes t R/L and h R/L , respectively. At the linearized order, one finds their equations of motion coupled together among the same polarizations, written in the Fourier space as [7], where x ≡ k/aH and ψ R,L (t, k) are the mode functions of the canonical gravitational wave, ψ ij ≡ aM Pl h ij /2.
While t R/L are sourced by ψ R/L in principle, the former is always parametrically larger than the latter for our concern, and thus ignoring the right-hand side of (11) is a justified approximation. We have also neglected slowroll suppressed and subdominant terms in (11) and (12).
Here, ξ(t) ≡ λχ/2f H is well approximated by m Q + m −1 Q in the slow-roll regime. Without loss of generality m Q is assumed to be positive, and then t R becomes unstable for x max > x > x min , with x max,min ≡ m Q +ξ ±(m 2 Q +ξ 2 ) 1/2 . Assuming m Q = const., we obtain the homogeneous solution to (11) as where W β,α (z) is the Whittaker function with α ≡ −i 2m Q ξ − 1/4 and β ≡ −i(m Q + ξ). We have used the WKB solution in the sub-horizon limit, t R (k/aH → ∞) = (2k) −1/2 (2x) β e ix , as the initial condition. Then t R is amplified around the horizon crossing by the factor of e π 2 (m Q +ξ) W β,α (−2ix min ) ≈ e 1.85m Q , while it decays as matter, ρ t ∝ a −3 i.e. t R ∝ a −1/2 , on super-horizon scales. The source term for ψ R/L reads and the generated t R sources ψ R , producing additional GW. Using (13), one can obtain the sourced ψ R by using Green's function method, giving the GW power spectrum where F 2 ≈ 2e 3.62m Q and its full expression can be found in [7]. Note that (13) and (15) assume constant B , m Q and ξ, while to determine their values and time variations one needs to solve the background dynamics, (7)-(10).

III. CHECKLIST
In order to settle the final allowed strength of GW signals from this model, we need to ensure some computational and observational consistencies. We list them and show the resulting parameter region in the following subsections. Since this constraints are proportional to e 1.85m Q as mentioned in the main text, the coupling constant g shown in the plot is rescaled by this factor. In the upper yellow shaded region, the backreaction is expected to be strong and disrupts the background evolution. In the lower blue shaded region, the energy fraction of the gauge field is significant enough to make the scalar spectral index becomes too red beyond the 2σ region of Planck constraints for r = 10 −3 . The black dotted contours for the values of H inf are superimposed in the case with r = 10 −3 .

A. Backreaction
The produced t R (13) backreacts on the background dynamics through eqs. (7)-(10) with the terms where we ignore the sub-leading backreaction from t L or ψ R,L . We first estimate these contributions analytically. Using (13) and background relation ξ ∼ = m Q + m −1 Q and changing variables into x = k/aH with the integration domain 0 < x < x max , one can write |ρ t | = H 4 I ρ (m Q ), |T χ BR | = λH 4 I χ (m Q )/f and T Q BR = gH 3 I Q (m Q ), where all the I's approximately follow I ρ,χ,Q ∝ e 3.7m Q . For a given value of g, these terms would easily dominate (7), (9) and (10) for large m Q , if one took m Q as a free parameter. However, this would infer that strong backreaction prevents the system from reaching such a parameter region. The conditions to ensure that each of ρ t and T χ,Q BR is subdominant in (7), (9) and (10) are translated into upper bounds on g, where G ρ,χ,Q ∝ I −1/2 ρ,χ,Q . In Fig. 1, we show the strongest constraints coming from G χ , though they are almost degenerate.
For large m Q , the backreaction is not completely negligible even in the allowed region shown in Fig. 1. In those cases, one has to resort to full numerical calculations simultaneously solving all equations of motion for background fields, (8)-(10) and for perturbations, (11) and (12) with full source terms included. Fig. 2 shows our numerical result for the following parameters: where the corresponding maximum of m Q is around 44.
The tensor-to-scalar ratio r R = P h R /P ζ (k * ) where k * is the pivot scale for CMB observations indeed exceeds the detectable limit 10 −3 even with such a extremely low inflationary energy scale ∼ 36 MeV.

B. Curvature Perturbation
Previous attempts to generate GW from scalar or vector fields are tightly constrained by the CMB observation on the curvature perturbation ζ [5,10]. In our model, the inflaton fluctuation δφ is assumed to be responsible for generating ζ compatible with the CMB observation. Contributions from the other scalar modes δχ, δQ and M to ζ are negligible, unless χ becomes a curvaton [7].
In addition, we investigate another channel in which the second order effect of t R produces the inflaton perturbation, t R t R → δφ, through the gravitational interaction. This effect arises only at the second order due to the absence of linear couplings between δφ and t R , while the sourcing of t R to the GW is first-order, thus ζ (s) = −Hδφ (s) /φ ∝ (t R /M Pl ) 2 is expected to be negligible for the parameter range of our interest. We will address this effect in detail in the upcoming work.
Even though the part of δφ sourced by the second order of t R or the linear order of the scalar perturbations in the axion-SU (2) sector only has negligible effects, that of δφ originated from its own vacuum fluctuations can be influenced by the background fields χ and Q, due to their contribution toḢ. As a result, the spectral index in our model reads, where in the last step we have used A B φ , χ , true with the parameters of our interest. The Planck measures n s = 0.9645 ± 0.0049 [4], and without assuming an accidental cancellation between B and η φ , we require a bound on B as where t * denotes the time of the horizon crossing of the CMB modes. Note that this constraint can be relaxed if η φ is positive. When (22) saturates, in our model, B can explain the red-tilted curvature perturbations without a huge hierarchy of slow-roll parameters η φ φ . It is a quite intriguing possibility for small-field inflationary models since all slow-roll parameters are naively expected to be equivalently small in that class of inflation. We numerically checked that n S (k * ) within 2σ of Planck constraints is realized solely by B for the parameters (20).
The bound (22) is translated into a lower bound on We plot this as the light-blue shaded region in Fig. 1.

C. Perturbativity
Since the amplitude of t R is substantially amplified due to the instability in our model, we need to ensure that it does not invalidate our perturbative calculation. We thus impose that the 1-loop contribution to the two-point function t R t R should be negligible to that of the tree level. The terms −F a µν F aµν /4 + λχF a µνF aµν /(4f ) lead to three-and four-point vertices, and it can be shown that their one-loop diagrams give contributions of the same order [10]. We here focus on the latter and demonstrate that the perturbativity condition gives no additional bounds on the model parameters. The four-point interaction Hamiltonian readŝ H (4) giving rise to, using the in-in formalism, where we ignored the left-handed mode. Defining R t as the ratio of (25) divided by the tree-level contribution, (2π) 3 δ(k + k )|t k (τ )| 2 , evaluated at the time when t R reaches its maximum value, we ensure R t 1 to safely ignore the higher-order loops and to justify the perturbative approach. Evaluating (25) with (13), we verify that R t 1 is satisfied up to m Q = 50 for r = 10 −3 .

IV. CONCLUSION
The main message of this Letter is that the detection of primordial gravitational waves does not necessarily exclude low-energy inflation. Once an SU (2) gauge field The energy densities of the inflaton ρ φ (blue), the axion ρχ (yellow), the SU (2) gauge field ρA (green) and the backreaction from the amplified perturbation ρt R (red) are shown. The inflationary energy scale is as low as ρ φ ≈ (30MeV) 4 . The horizontal axis represents the backward e-folds, and N * corresponds to the scale at which the tensor perturbation is maximal. This scale can be identified as the CMB pivot scale since the inflaton sector is independent of the gauge field sector. Although ρt R is negative while the instability is getting stronger, the total energy of the SU (2) gauge field is always positive. (Right panel) The tensor-to-scalar ratio r due to the sourced GW only with the right-handed polarization. On the CMB scale k * it exceeds the threshold value r = 10 −3 (yellow dashed line) and thus it is detectable for the upcoming CMB missions.
has a background configuration that respects the spatial rotation, its perturbations are coupled to the GW at the linear order. The former is amplified by instabilities around the horizon crossing, whose power is then linearly transferred to the latter. We have demonstrated that the GW power spectrum produced from this mechanism can be as significant as at detectable levels respecting all the consistency conditions, even if the inflationary energy scale is close to the BBN bound. Having a possible alternative source of GW, it is crucial to discriminate the generation mechanism of primordial GW to reveal the true energy scale of inflation. Fortunately, our model has the following distinct predictions to be distinguished from the conventional vacuum GW. (i) The fully parity-violating GW may be detected through CMB temperature and B-mode (TB) or E-mode and Bmode polarization (EB) cross-correlation by the upcoming satellite mission such as LiteBIRD [11]. (ii) Our model produces a sizable tensor non-Gaussianity with a particular shape [12]. (iii) The conventional consistency relation, n T = −r vac /8, is broken, where n T is the tensor spectral index. With the future observation, these signatures will carry important information for rigorous determination of inflationary energy scale.