# Single field inflation with modulated potential in light of the Planck and BICEP2

## Abstract

The recently released BICEP2 data detected the primordial B-mode polarization in the Cosmic Microwave Background (CMB) map which strongly supports for a large tensor-to-scalar ratio, and thus, is found to be in tension with the Planck experiment with no evidence of primordial gravitational waves. Such an observational tension, if confirmed by forthcoming measurements, would bring a theoretical challenge for the very early universe models. To address this issue, we in the present paper revisit a single field inflation model proposed in (1); (2) which includes a modulated potential. We show that this inflation model can give rise to a sizable negative running behavior for the spectral index of primordial curvature perturbation and a large tensor-to-scalar ratio. Applying these properties, our model can nicely explain the combined Planck and BICEP2 observations. To examine the validity of analytic calculations, we numerically confront the predicted temperature and B-mode power spectra with the latest CMB observations and explicitly show that our model is consistent with the current data.

###### pacs:

98.80.Cq###### pacs:

98.80.-k, 98.80.Cq## I Introduction

The inflationary hypothesis of the very early universe, since was proposed in the early 1980s (3); (4); (5) (see also (6); (7); (8) for early works), has become the dominant paradigm for understanding the initial conditions for the hot big bang cosmology. Within this context, a well established picture of causally generating cosmological perturbations in the primordial epoch has been greatly developed theoretically and observationally in the past decades. In particular, a significant prediction of nearly scale-invariant power spectra of primordial density perturbations based on the inflationary paradigm has been verified to high precision by the CMB observations in recent years (9); (10). Inflationary cosmology also predicted a nearly scale-invariant power spectrum of primordial tensor perturbations (11), which can give rise to the CMB B-mode polarization as detected by the BICEP2 collaboration (12). Assuming that all polarization signals were contributed by inflationary gravitational waves, the BICEP2 experiment implies that a nonzero value of the ratio between the spectra of tensor and scalar modes, dubbed as the tensor-to-scalar ratio , has been discovered at more than confidence level (CL) with a tight constraint as: at CL ( with foreground subtracted). This observation, if eventually verified by other ongoing experiments, implies a large amplitude of primordial gravitational waves and hence has significant theoretical implications on various early universe models.

Such a large amplitude of primordial tensor spectrum as indicated by the BICEP2, however, is in certain tension with another CMB experiment, the Planck result with at CL. Thus, the combination of these two data sets leads to a critical challenge for theoretical models of the very early universe. In the literature, there are some discussions of this issue from either the perspective of new physics beyond the inflationary CDM paradigm in which more degrees of freedom are introduced, see e.g. (13); (14); (15); (16); (17); (18); (19); (20), or from the propagations of photons after decoupling, see e.g. (21). However, it remains to be interesting to investigate the possibility of resolving such an observational challenge within the framework of single field inflation.

In order to address the issue of observational tension, we take a close look at prior assumptions made by the Planck and BICEP2. The Planck’s result of was obtained by, the assuming a constant spectral index of primordial curvature perturbations. As was pointed out in (22); (9); (12), however, if a nonzero running of the spectral index is allowed in the data analysis, there exist reasonable degeneracies among the spectral index , the running of the spectral index , and the tensor-to-scalar ratio . For example, an allowance of the running spectral index can at most enhance the upper bound of the tensor-to-scalar ratio of the Planck data to at CL (23)(9)(12) (see also (24) for theoretical discussions). Accordingly, a possibly existing running of the spectral index can efficiently circumvent the tension issue between the Planck and BICEP2 data. This phenomenological scenario, however, is not easy to be achieved in usual single field slow-roll inflation models. Therefore, we in the present paper study a notable mechanism of generating a negative running behavior for the spectral index of primordial curvature perturbations. In particular, we analyze a model of single field inflation as proposed in (1); (2), of which the inflaton field has a modulated potential. We will show that this model gives rise to a negative value of and also yields a relatively large value of . Therefore, our model can provide a theoretical interpretation in reconciling the tension issue existing between the Planck and BICEP2 measurements.

The present paper is organized as follows. In Sec. II we briefly review the single field inflation model with modulated potential and discuss its theoretical motivation from the physics of extra dimensions. Then, in Sec. III we perform analytic and numerical calculations on the background dynamics as well as theoretical predictions on primordial power spectra, respectively. We show that this model generally produces a large amplitude of tensor-to-scalar ratio and a negative running spectral index. Afterwards, we confront the theoretical predictions of this model with the combined Planck and BICEP2 data under a class of fixed parameter values in Sec. IV. Our numerical computation nicely demonstrate the model can explain these observations consistently. We conclude in Sec. V with a general discussion.

## Ii Single field inflation model with modulated potential

The single field inflation model with modulated potential was first proposed in Refs. (1); (2). In these earliest papers, the modulation acts as a rapid oscillating term added to the so-called natural inflation potential (25), which can help generating large running spectral index. In this section, we will firstly take a brief review of the development of this kind of model.

Natural inflation model was motivated by the idea to connect cosmic inflation and particle physics, where the inflaton is considered as a Pseudo-Nambu-Goldstone Boson (PNGB) from the spontaneous breaking of a global symmetry (25). In this model, shift symmetry has been introduced to make the inflaton potential flat and stable, and to maintain sufficient e-folding numbers (26). Similar to the axion (27), the inflaton potential has the form

(1) |

where is the scale of spontaneous symmetry breaking and the cosine term is thought to be produced by some non-perturbative effects which break the symmetry explicitly at a relative low scale . So the inflaton has the mass at the order of . However, the flatness condition of the inflaton potential requires the scale of spontaneous symmetry breaking to be larger than the Planck scale and it is expected that at such a high scale the global symmetries are violated explicitly by the quantum gravity effects (28); (29). These effects introduces a modulation to the potential, which was considered in (1) from the viewpoint of effective field theory. By considering the higher dimensional operators without derivatives (due to the global symmetry breaking at ), the model proposed in (1) has the potential

(2) |

where is a small number and is large, the phase is physically unrelevant and can be set to zero. It is also possible to add a constant to the above potential to make its minimum vanish. A large will modulate the potential with rapid oscillations and superimpose a series bumps into the otherwise featureless potential. But for sufficiently small , the amplitude of the oscillations can be controlled to be small to protect the overall picture of inflation. The slow-roll conditions are violated mildly and the predicted scalar power spectrum has strong oscillations with sizeable running but still allowed by the observations, as shown in (1); (30). Another prediction of this model is the enhanced wiggles in the matter power spectrum. These features are possible to be detected by future experiments (31).

In 2003, Arkani-Hamed et al. proposed the extra-dimensional version of natural inflation(extranatural inflation) (32), where a five-dimensional Abelian gauge field is considered and the fifth dimension is compactified on a circle of radius . The extra component propagating in the bulk is considered as a scalar field from the four-dimensional view and the inflaton is identified as the gauge-invariant Wilson loop , with a 5D gauge coupling constant . The 4D effective Lagrangian at energies below can be written as:

(3) |

with the four-dimensional effective gauge coupling constant. The non-local potential is generated in the presence of fields charged under the Abelian symmetry in the bulk (34); (33), for a massless field with charge the potential for the Wilson loop is

(4) |

where the “” and “” represent the bosonic and fermionic fields respectively and we have neglected higher power terms. By defining and adding a constant term, the potential has almost the same form with that natural inflation (1) and the effective decay constant is (we have assumed is of order unity) . There are some advantages of this model compared with the old natural inflation based on the 4D PNGB, such as that the effective decay constant can be naturally greater than for a sufficiently small coupling constant , and gravity-induced higher-dimensional operators are generally exponentially suppressed as long as the extra dimension is larger than the Planck length due to the extra dimension nature. However, as same as the natural inflation, this model predicts a scalar spectrum with negligible running of the spectral index. In order to have large running index, in Ref.(2) the authors generalized this model to the case including multiple charged fields under the Abelian symmetry. If we simply consider one massless and one massive fields coupled to gauge field, i.e., and , as studied in detail in (2), the effective potential for the inflaton then becomes:

(5) |

where again we have neglected higher power terms, are the charges of these two fields, and is related to as

(6) |

where the constant for the bosonic and fermionic fields respectively. For a large mass , has a small value. If the ratio of these two charges , we get the same potential as the model (2) mentioned above. It was found in (2) that this model can produce a scalar spectrum with negative running as large as .

In all, the model with modulations (2) or (5) are well-motivated, and different from the simplest slow-roll inflation models, it produces an oscillating scalar spectrum with significant running. However, the tensor-to-scalar ratio produced in (2) is very small which cannot be consistent with the BICEP2 data. In this paper we will investigate whether this model can produce a large suggested by BICEP2 and at the same time a sizeable running to alleviate the tension between the Planck and BICEP2 data. Note that some related studies has been done in Ref. (35), and here we will revisit this problem in more detail.

## Iii Inflationary dynamics with non-vanishing running spectral index

In this section we perform the analytical and numerical analyses of the inflationary solution described by this model. In particular, we analyze the dynamics of the slow roll parameters during inflation. Consider a canonical scalar field with the potential given by (5).

The inflationary dynamics can be characterized by a series of slow roll parameters, of which the expressions are given by,

(7) |

where we have introduced

(8) |

Note that, inflation requires the above slow roll parameters to be much less than unity. Accordingly, the inflationary e-folding number follows:

(9) |

where we have applied the approximation . Following the standard procedure of inflationary perturbation theory (36); (37), the power spectra of primordial curvature and tensor perturbations of this model can be expressed as:

(10) |

and correspondingly, the tensor-to-scalar ratio is defined as

(11) |

Moreover, one can define the spectral index of primordial curvature perturbations and the associated running spectral index as follows,

(12) | |||

(13) |

In regular inflation models, the slow roll parameters scale as: and during inflation. Moreover, the spectral index is of order and is of order . Thus, a large running behavior of the spectral index is difficult to be achieved due to the suppression effect by . However, it is interesting to notice that, in the model under consideration, the slow roll approximations can be slightly broken for a short while due to the inclusion of the rapid oscillating term in the potential.

In our model, inflation ceases when reaches with , and one can have the initial value for the inflaton to be the value at the moment of Hubble-crossing. Near the Hubble-crossing, one can get the expressions for the slow roll parameters approximately,

(14) |

When , and , from the above approximation one can get the value of in the same order of and , and hence a relatively large running behavior can be obtained.

In the following we perform the numerically calculation of our model. In Fig.1 we plot the evolutions of the slow-roll parameters , and with respect to the e-folding number . In our figure plot, inflation begins from the right side where the slow-roll parameters are small, and ends at the left side where they present some oscillatory behavior with their amplitude approaching 1. The pivot scale, which corresponds to , crosses the Hubble radius at the time when , marked with black dotted line. The numerical results depends only on three parameters in the model, namely , and . In the numerical calculation, we take three groups of parameter choices (see the caption), and in order to have a comparison, we also plot the cases of natural inflation model. One could see that at the pivot scale is almost of the same order as and (To help see more clearly, we also plot the zoomed-in figures around the pivot scale, with the vertical coordinates of the same range.) in our model, while is negligible in natural inflation model. Therefore, as has been analyzed above, one can observe a considerable running behavior of the spectral index around this point.

One can directly relate the slow-roll parameters with the perturbation variables of a canonical single field inflation model. In Fig.2 we plot the evolution of the spectral index of scalar perturbation , the running of the spectral index , and the tensor-to-scalar ratio . In the plot, we take the range from to which is able to cover the range () used in Planck and BICEP2 paper. We also marked with a vertical dotted line the pivot scale, which reenters the Hubble radius and eventually can be observed by today’s experiments. From the plot we can see that both our model and natural inflation can give a large , as needed by the BICEP2’s data. However, one obvious difference between the two models is that our model is able to yield a negative running spectral index for the power spectrum of scalar perturbations, roughly of the order that can be applied to reconcile the Planck and BICEP2 data (10); (12), while the running behavior from the model of natural inflation is negligible. Such a considerable running behavior obtained in our model arises from a large-valued parameter , which bring to the same order of and so that it has dominant contribution to the expression (13) of in comparison with other two terms.

A sizable negative running has the possibility to make the spectral index vary efficiently with scales, e.g., from blue tilt to red tilt. This could lead to some observable features on the power spectrum , namely, a bump would appear on the plot, or the amplitude on small region might get suppressed, which can be useful in the explanation of small anomaly. In Fig.3 we plot the amplitudes of scalar perturbations under various parameter choices. We can see that, although at the pivot scale the amplitudes in these cases are almost the same, which are consistent with the data, they can be very different at small regions.

## Iv Fitting the cosmological data

With the analyses performed in the above section, we have shown that our model can indeed have a sizable negative running as well as a large . In this section, we directly confront our model to the observational data to see how it reconciles the Planck and BICEP2 data.

In Fig. 4 we present our results in plot. We plot our model with two groups of parameter choices (blue and red), which can both fit the Planck+BICEP2 data very well. For each choice, we considers two cases which inflation continues for 60 (solid lines) and 50 e-foldings (dashed lines). For a comparison, we also plot natural inflation models with and , presented with magenta lines. The lines grows as grows, making its prediction of and more and more close to our model, and also more and more close to the allowed space by the contour. In this plot, we have chosen the pivot scale as .

We have also showed the TT and BB spectrum in Fig. 5 and Fig. 6, with all the color lines have the same parameter-correspondence as in Fig.2. We see that on the large region all the lines glues together indicating a degeneracy of the parameters, and fit the data very well. On small regions, the lines deviate from each other, but since the error bars on this region are quite large, the lines are still in consistency with the data. However, one might also notice that when the line best fits the Planck temperature spectrum, i.e. the cyan one, gives smaller auto-correlation compared with BICEP2’s data. Conversely, the line which fits BICEP2’s data much better (the green one) gives larger temperature power spectrum which fails to explain the small anomaly. This phenomenon can be easily understood: the scalar and tensor spectrum are linked by tensor-scalar ratio , which won’t change too much with in our model, as can be seen in Fig.2. Therefore, a raising/lowering of scalar spectrum at large scales corresponds to the same behavior of tensor spectrum. Furthermore, we also plot the red line as an intermediate case, which will not deviate too much from the data points in either TT or BB spectrum. We hope the global fitting of full parameter space can provide us a better parameter choice for both Planck TT spectrum and BICEP2 BB spectrum, which we will leave for future investigations.

## V Conclusions

The BICEP2 group has released the results of the CMB polarization measurement, which strongly hints to an existence of a large amplitude of primordial gravitational waves. This result, however, is in tension with the Planck data released last year when they are interpreted by the standard -parameter (without ) model. One simple approach to alleviate this tension is to take into account the running of the spectral index of the curvature perturbations, which in our paper is characterized by . The inclusion of this parameter can greatly relax the observational constraint on the tensor-to-scalar ratio due to their degeneracy (22). Recently, a numerical global simulation of the model reveals that the combined Planck and BICEP2 data favor a negatively valued running spectral index of (23).

This observational implication, while puts forward a challenge to slow roll inflation models, can be nicely implemented by a single field inflation with modulated potential (1); (2) as demonstrated in the present paper, where we treat the modulation as a rapid oscillating term. This is because, although the whole inflationary dynamics is dominated by the regular slow roll part of the potential, this rapid oscillating term can relatively violate the slow roll approximation during some local evolutions. In the specific model considered, we explicitly show that the parameter which is associated with the running and a higher order slow roll parameter in the normal slow roll inflation model, can be enhanced to the value as large as the first order slow roll parameters and . Therefore this model can give rise to a considerable and negative running spectral index. In this paper, we performed the numerical calculation of the model in detail by solving the dynamics of slow-roll parameters and perturbation variables, and then fitted them to the combined Planck and BICEP2 data. From the numerical results, one can easily see that, with a large-valued parameter near the pivot scale, the spectrum index can be changed from value larger than to value smaller than performing a negative running feature. Due to this running behavior, the power spectrum can be suppressed at small region. At the same time this model can also produce gravitational waves of large amplitudes as long as the effective decay constant is large enough. Hence the inflation model under consideration provides a consistent interpretation of the combined Planck and BICEP2 data.

The present model has further implications for the observations. As we have mentioned the modulation of the rapid oscillating term could amplify the wiggles of the CMB temperature spectrum (30) and the matter power spectrum and thus is of observable interest to the future experiments (31). Although this model is motivated by observational phenomena, it deserves mentioning that this model has interesting connections with other inflation models. For example, we have mentioned that our model can reduce to a natural inflation model when the modulation term is small enough. The extranatural inflation including higher power terms has been recently investigated in (38). The axion-monodromy inflation with modulations was studied in (39); (40), while its supergravity version was discussed in (41). In the literature, there are other studies on deriving a large running of the spectral index from various approaches, for instances, see (42); (43); (44); (45); (46); (47); (48); (49); (50); (51); (52); (53); (54); (55); (56); (57); (58); (59).

Additionally, the wiggles in power spectrum can induce features on non-Gaussianities, especially of the squeezed shape, since from the consistency relation we roughly have (60), where wiggles in may affect (61). The features on non-Gaussianities are expected to be detected by the future observations. We will discuss these issues as a sequence of this work in a future project.

As a final remark, note that there may be other approaches of reconciling the tension between the Planck and BICEP2 data within the framework of inflationary cosmology, such as to suppress the scalar spectrum at large scales by a double field inflation model (44) or using a step-like process (58). The existence of nontrivial tensor spectral index may also work, which needs to be accompanied by simulation of model. To address this issue, we would like to numerically scan the full parameter space and check all available regions allowed by observations, which will be the future topic.

###### Acknowledgements.

We are grateful to Junqing Xia and Hong Li for helpful discussions. YW, SL and XZ are supported by NSFC under grants Nos. 11121092, 11033005, 11375220 and also by the CAS pilotB program. ML is supported by Program for New Century Excellent Talents in University and by NSFC under Grants No. 11075074. CYF is supported in part by physics department at McGill university.### References

- X. Wang, B. Feng, M. Li, X. -L. Chen and X. Zhang, Int. J. Mod. Phys. D 14, 1347 (2005) [astro-ph/0209242].
- B. Feng, M. -Z. Li, R. -J. Zhang and X. -m. Zhang, Phys. Rev. D 68, 103511 (2003) [astro-ph/0302479].
- A. H. Guth, Phys. Rev. D 23, 347 (1981).
- A. D. Linde, Phys. Lett. B 108, 389 (1982).
- A. Albrecht and P. J. Steinhardt, Phys. Rev. Lett. 48, 1220 (1982).
- A. A. Starobinsky, Phys. Lett. B 91, 99 (1980).
- L. Z. Fang, Phys. Lett. B 95, 154 (1980).
- K. Sato, Mon. Not. Roy. Astron. Soc. 195, 467 (1981).
- P. A. R. Ade et al. [Planck Collaboration], arXiv:1303.5076 [astro-ph.CO].
- P. A. R. Ade et al. [Planck Collaboration], arXiv:1303.5082 [astro-ph.CO].
- A. A. Starobinsky, JETP Lett. 30, 682 (1979).
- P. A. R. Ade et al. [BICEP2 Collaboration], arXiv:1403.3985 [astro-ph.CO].
- V. Miranda, W. Hu and P. Adshead, arXiv:1403.5231 [astro-ph.CO].
- J. -F. Zhang, Y. -H. Li and X. Zhang, arXiv:1403.7028 [astro-ph.CO].
- J. -Q. Xia, Y. -F. Cai, H. Li and X. Zhang, arXiv:1403.7623 [astro-ph.CO].
- D. K. Hazra, A. Shafieloo, G. F. Smoot and A. A. Starobinsky, arXiv:1404.0360 [astro-ph.CO].
- K. M. Smith, C. Dvorkin, L. Boyle, N. Turok, M. Halpern, G. Hinshaw and B. Gold, arXiv:1404.0373 [astro-ph.CO].
- Y. -F. Cai and Y. Wang, arXiv:1404.6672 [astro-ph.CO].
- Z. -G. Liu, H. Li and Y. -S. Piao, arXiv:1405.1188 [astro-ph.CO].
- Y. -F. Cai, arXiv:1405.1369 [hep-th].
- B. Feng, M. Li, J. -Q. Xia, X. Chen and X. Zhang, Phys. Rev. Lett. 96, 221302 (2006) [astro-ph/0601095]; M. Li and X. Zhang, Phys. Rev. D 78, 103516 (2008) [arXiv:0810.0403 [astro-ph]]; W. Zhao and M. Li, arXiv:1402.4324 [astro-ph.CO]; W. Zhao and M. Li, arXiv:1403.3997 [astro-ph.CO]. S. Lee et al, arXiv: 1403.5585.
- H. Li and J. -Q. Xia, JCAP 1211, 039 (2012) [arXiv:1210.2037 [astro-ph.CO]].
- Hong Li, Jun-Qing Xia, and Xinmin Zhang, [arXiv:1404.0238 [astro-ph.CO]].
- Y. Gong, arXiv:1403.5716 [gr-qc].
- K. Freese, J. A. Frieman and A. V. Olinto, Phys. Rev. Lett. 65, 3233 (1990); F. C. Adams, J. R. Bond, K. Freese, J. A. Frieman and A. V. Olinto, Phys. Rev. D 47, 426 (1993) [hep-ph/9207245].
- E. D. Stewart, Phys. Lett. B 391, 34 (1997) [hep-ph/9606241]; Phys. Rev. D 56, 2019 (1997) [hep-ph/9703232]; S. Dodelson and E. Stewart, Phys. Rev. D 65, 101301 (2002) [astro-ph/0109354].
- R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38, 1440 (1977); R. D. Peccei and H. R. Quinn, Phys. Rev. D 16, 1791 (1977).
- M. Kamionkowski and J. March-Russell, Phys. Lett. B 282, 137 (1992); R. Holman, et al., Phys. Lett. B 282, 132 91992); S. M. Barr and D. Seckel, Phys. Rev. D 46, 539 (1992).
- R. Kallosh, A. D. Linde, D. A. Linde and L. Susskind, Phys. Rev. D 52, 912 (1995) [hep-th/9502069].
- J. Liu, H. Li, J. Xia and X. Zhang, JCAP 0907, 017 (2009) [arXiv:0901.2033 [astro-ph.CO]].
- C. Pahud, M. Kamionkowski and A. R.Liddle, Phys. Rev. D 79, 083503 (2009) [arXiv:0807.0322 [astro-ph]].
- N. Arkani-Hamed, H. -C. Cheng, P. Creminelli and L. Randall, Phys. Rev. Lett. 90, 221302 (2003) [hep-th/0301218]; N. Arkani-Hamed, H. -C. Cheng, P. Creminelli and L. Randall, JCAP 0307, 003 (2003) [hep-th/0302034].
- A. Delgado, A. Pomarol and M. Quiros, Phys. Rev. D 60, 095008 (1999) [hep-ph/9812489].
- Y. Hosotani, Phys. Lett. B 126, 309 (1983); 129, 193 (1983).
- M. Czerny and F. Takahashi, arXiv:1401.5212 [hep-ph]; M. Czerny, T. Higaki and F. Takahashi, arXiv:1403.0410 [hep-ph].
- E. D. Stewart and D. H. Lyth, Phys. Lett. B 302, 171 (1993) [gr-qc/9302019].
- V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, Phys. Rept. 215, 203 (1992).
- K. Kohri, C. S. Lim and C. -M. Lin, arXiv:1405.0772 [hep-ph].
- T. Higaki, T. Kobayashi, O. Seto and Y. Yamaguchi, arXiv:1405.0775 [hep-ph].
- R. Flauger, L. McAllister, E. Pajer, A. Westphal and G. Xu, JCAP 1006, 009 (2010) [arXiv:0907.2916 [hep-th]]. R. Easther and R. Flauger, JCAP 1402, 037 (2014) [arXiv:1308.3736 [astro-ph.CO]]. T. Kobayashi, O. Seto and Y. Yamaguchi, arXiv:1404.5518 [hep-ph].
- R. Kallosh, A. Linde and B. Vercnocke, arXiv:1404.6244 [hep-th].
- J. E. Lidsey and R. Tavakol, Phys. Lett. B 575, 157 (2003) [astro-ph/0304113].
- M. Kawasaki, M. Yamaguchi and J. ’i. Yokoyama, Phys. Rev. D 68, 023508 (2003) [hep-ph/0304161].
- B. Feng and X. Zhang, Phys. Lett. B 570, 145 (2003) [astro-ph/0305020].
- D. J. H. Chung, G. Shiu and M. Trodden, Phys. Rev. D 68, 063501 (2003) [astro-ph/0305193].
- M. Bastero-Gil, K. Freese and L. Mersini-Houghton, Phys. Rev. D 68, 123514 (2003) [hep-ph/0306289].
- M. Yamaguchi and J. ’i. Yokoyama, Phys. Rev. D 68, 123520 (2003) [hep-ph/0307373].
- D. -S. Lee, L. -Z. Fang, W. -L. Lee and Y. -Y. Charng, J. Korean Phys. Soc. 45, S203 (2004) [astro-ph/0403055].
- N. Kogo, M. Sasaki and J. ’i. Yokoyama, Phys. Rev. D 70, 103001 (2004) [astro-ph/0409052].
- F. Paccetti Correia, M. G. Schmidt, and Z. Tavartkiladze, Nucl. Phys. B 739, 156 (2006) [hep-th/0504083].
- G. Ballesteros, J. A. Casas and J. R. Espinosa, JCAP 0603, 001 (2006) [hep-ph/0601134].
- Q. -G. Huang and M. Li, Nucl. Phys. B 755, 286 (2006) [astro-ph/0603782].
- R. Gonzalez Felipe and N. M. C. Santos, Phys. Rev. D 78, 023519 (2008) [arXiv:0711.0022 [hep-ph]].
- T. Matsuda, JHEP 0802, 099 (2008) [arXiv:0802.3573 [hep-th]].
- T. Kobayashi and F. Takahashi, JCAP 1101, 026 (2011) [arXiv:1011.3988 [astro-ph.CO]].
- M. Peloso, L. Sorbo and G. Tasinato, arXiv:1401.7136 [astro-ph.CO].
- M. Czerny, T. Kobayashi and F. Takahashi, arXiv:1403.4589 [astro-ph.CO].
- C. R. Contaldi, M. Peloso and L. Sorbo, arXiv:1403.4596 [astro-ph.CO].
- A. Ashoorioon, K. Dimopoulos, M. M. Sheikh-Jabbari and G. Shiu, arXiv:1403.6099 [hep-th].
- J. M. Maldacena, JHEP 0305, 013 (2003) [astro-ph/0210603]; C. Cheung, A. L. Fitzpatrick, J. Kaplan and L. Senatore, JCAP 0802, 021 (2008) [arXiv:0709.0295 [hep-th]].
- J. -O. Gong, K. Schalm and G. Shiu, Phys. Rev. D 89, 063540 (2014) [arXiv:1401.4402 [astro-ph.CO]].