# BICEP2 constrains Composite Inflation

## Abstract

In light of BICEP2, we reexamine single field inflationary models in which the inflaton is a composite state stemming from various four-dimensional strongly coupled theories. We study in the Einstein frame a set of cosmological parameters, the primordial spectral index and tensor-to-scalar ratio , predicted by such models. We confront the predicted results with the joint Planck data, and with the recent BICEP2 data. We constrain the number of e-foldings for composite models of inflation in order to obtain a successful inflation. We find the minimal composite inflationary model is fully consistent with the Planck data. However it is in tension with the recent BICEP2 data. The observables predicted by the glueball inflationary model can be consistent with both Planck and BICEP2 contours if a suitable number of e-foldings are chosen. Surprisingly, the super Yang-Mills inflationary prediction is significantly consistent with the Planck and BICEP2 observations.

PACS: 1.10.-z, 11.30.Rd, 98.80.-k

Keywords: Composite Inflation, Non-minimal Coupling, Strongly Interacting Field Theories,

Planck and Bicep2 Constraints.

July 14, 2014

BICEP2 constrains Composite Inflation

Phongpichit Channuie

Department of Physics, School of Science, Walailak University,

Thasala District, Nakhon Si Thammarat, 80160, Thailand

Email: phongpichit.ch@wu.ac.th

###### Contents:

## 1 Introduction

Nowadays, inflationary models gain a lot of interest. The inflationary paradigms [1, 2, 3, 4, 5] tend to solve important issues, e.g. the magnetic monopoles, the flatness, and the horizon problems, and subsequently provide the mechanism for generation of density perturbations as seed for the formation of large scale structure in the universe. Most models of inflation so far were formulated by introducing new scalar fields (called inflaton) with a nearly flat potential appearing in many paradigms, e.g. [6, 7, 8, 9, 10], and even in superstring [11, 12] and supergravity theories [13, 14, 15, 16].

However, the theories featuring elementary scalar fields are unnatural meaning that quantum corrections generate unprotected quadratic divergences which must be fine-tuned away if the models must be true till the Planck energy scale. Hence, this is the main reason why the authors in [17, 18, 19] investigated inflation in which the inflaton need not be an elementary degree of freedom. Recent investigations show that it is possible to construct models in which the inflaton emerges as a composite state of a four-dimensional strongly coupled theory. We called these alternative paradigms composite inflation. There were other models of super or holographic composite inflation [20, 21, 22, 23, 24, 25, 26, 27]. In order to constrain the inflationary theory, recent investigations contain hundreds of different scenarios [28].

Most frequently, we added a term in which the inflaton field, , non-minimally coupled to gravity to the action, say . This term has purely phenomenological origin. The reason resides from the fact that one want to relax the unacceptable large amplitude of primordial power spectrum generated if one takes or small. However, it is instructive here to provide some insight on how a non-minimal coupling can be naturally formed. An instructive analysis of the generated coupling of a composite scalar field to gravity has been initiated in the Nambu-Jona-Lasinio (NJL) model [29]. In this case, a composite field is a chiral condensate, , from the NJL model which couples to gravity via a coupling term very similar to that we deploy in this work. The presence of the induced non-minimal coupling of the boundstate object to gravity, with a coupling constant, and the (Ricci) scalar curvature, can be implemented from the symmetric phase of the theory, i.e. a massless case, or from the broken phase when providing the non-minimal term, where with the vacuum expectation value of . However, the results of are the same. In principle, we can transform such a term into another form by applying the conformal transformation (see detailed discussion in [30]). Among many frames, the Jordan frame and the Einstein frame are those discussed in the community. Roughly speaking, physics looks different in two different conformal frames. However, physical conclusions remains the same in a weak gravitational field limit [30].

In this work, we study observables, i.e. the scalar spectral index and tensor-to-scalar ratio , in the Einstein frame predicted by the composite inflationary models recently proposed. We will briefly present the concept of conformal transformation. We then compute the power spectrum for the primordial curvature perturbations and the tensor-to-scalar ratio for various composite inflationary models in the Einstein frame. Next we constrain the models by placing our results into the plane implemented by using the observational bound for and from the Planck data and also confront the results with the recent BICEP2 data. Finally, some comments about our results are made in the last section.

## 2 Composite Setup and Conformal Transformation

We consider a generic strongly-interacting theory non-minimally coupled to gravity. For the model of inflation, we identify the inflaton with one of the lightest composite states of the theory and denote it with . This state has mass dimension . The action for the composite models studied below can be written in the general form

(1) |

where and can in general be an arbitrary function of the composite field . gives rise to the non-minimal coupling . To recover the ordinary Einstein theory of gravity at low energy, we deduce where is the vacuum expectation value of the field at the end of inflation and is the reduced Planck energy scale. According to the models of single-field inflation [17, 18, 19] in which the inflaton is a composite state stemming from various strongly-interacting theories well-known in particle physics, in this work, we write and in the form as follows:

(2) |

where the non-minimal coupling to gravity is controlled by the dimensionless coupling . With the non-minimal coupling term phenomenologically introduced, one need to diagonalize into another form by applying a conformal transformation. The conformal transformation transforms a metric tensor into another metric according to the rule

(3) |

such that

(4) |

Hereafter we drop tildes to express variables in the Einstein frames. Applying the conformal transformation leads to the Einstein frame and the action reads

(5) | |||||

where prime denotes derivative over the field . We arrived at an involved kinetic term for the inflaton. It is convenient to introduce a canonically normalized field related to via

(6) |

In terms of the canonically normalized field we come up with the standard fashion

(7) |

In this work, we will examine the dynamics in the Einstein frame, and therefore express the slow-roll parameters in terms of and :

(8) | |||||

(9) |

where “ ” denotes derivative with respect to the field . Slow-roll inflation ends when , with the corresponding field value . The horizon exist when the field value equals which is determined by the number of e-foldings . This parameter is given by

(10) |

Here we express the number of e-foldings for the change of the field (or equivalently ) from to and will use the above mathematical framework to examine the composite-field models in Sec.(4).

## 3 Inflationary Observables

In this section, we will briefly review the study of the scalar and tensor perturbations. This leads to obtaining the required inflationary predictions, e.g. the spectral index and tensor-to-scalar ratio. Our starting point here is the action given in Eq.(7). We will derive the parameters by closely following the work presented in [37] for “Generalised G-inflation” using the unitary gauge. In the unitary gauge , it is kind of tradition to begin with the perturbed metric as [31]

(11) |

where is a lab function, and is the shift vector which are respectively defined as:

(12) |

where and scalar perturbations and is a tensor perturbation satisfying the condition . Physically, the lapse function represents the rate of flow of proper time with respect to the coordinate time. Substituting the perturbed line-element given in Eq.(11) into the action (7) and expanding it to the second order, we obtain [37]

(13) |

Note that with the definition in [37] in this case the sound speed equals unity, . The tensor perturbations can be canonically renormalised via the following re-definitions:

(14) |

Plugging these new parameters, the above quadratic action becomes

(15) |

where a prime represents derivative with respective to the conformal time . It is common to transform the action and re-express the resulting one in terms of the Fourier modes, and one can basically show that the evolution equation from the action in Eq. (15) is given by

(16) |

The normalized solution to this perturbation equation can be conventionally written in terms of the Hankel function so that we obtain the power spectrum of the primordial tensor perturbation at the horizon exits [37] as

(17) |

where

(18) |

and the tensor spectral tilt is given by . Now we turn to considering the scalar fluctuations. This can be achieved by setting in the perturbed line-element mentioned in Eq. (12). Here we proceed in the same manner as for the tensor perturbations. For the scalar perturbations, we arrive with the second order action [37]

(19) |

where is already given in Eq. (8), and the sound speed is also unity, , in this case. The above action can be canonically renormalised via the following re-definitions:

(20) |

Plugging these new parameters, the above quadratic action becomes

(21) |

where a prime represents derivative with respective to the conformal time . It is common to transform the action and re-express the resulting one in terms of the Fourier modes, and one can basically show that the evolution equation from the action (21) is given by

(22) |

The normalized solution to this perturbation equation can be conventionally written in terms of the Hankel function so that we obtain the power spectrum of the primordial scalar perturbation at the horizon exits [37] as

(23) |

where

(24) |

and the spectral index is given by . Moreover, the tensor-to-scalar ratio is defined by

(25) |

Recently, the Planck satellite data showed that the spectral index of curvature perturbations is constrained to be ( CL) and ruled out the exact scale-invariance () at more than confident level (CL), whilst the tensor-to-scalar ratio is bounded to be ( CL). Most recently, the observed -mode power spectrum provides the tensor-to-scalar ratio with disfavoured at CL [43]. These constraints are used to falsify the most popular and simple inflationary models.

## 4 Predictions from Composite Inflation

According to the composite paradigms considered in this work, it was potentially shown that the models of composite inflation nicely respect tree-level unitary for the scattering of the inflation field during inflation all the way to the Planck scale. In this section, we consider composite inflationary models that can be described by the primordial power spectrum observables consisting of the spectral index, , and the tensor-to-scalar ratio, , derived in the Einstein frame. In the following computations, we assume a large non-minimal coupling and a large field approximation, i.e. , is also implemented. For later convenience, we express the inflationary parameters in terms of the number of e-foldings..

Minimal Composite Inflation (MCI)

The first model of composite inflation we will examine is recently investigated in [17]. According to this paradigm, they engaged the simplest models of technicolor passing precision tests well known as the minimal walking technicolor (MWT) theory [33, 34, 35, 36] with the standard (slow-roll) inflationary paradigm as a template for composite inflation. The inflaton field in this case is the lightest composite state emerging from a bilinear condensate of techni-quarks with . With the large field approximation, the action in the Jordan frame reads [17]

(26) |

where is a self coupling which is constrained by the underlying theory to be of the order of unity. For this model, we have

(27) |

Here we can diagonalize the inflaton-gravity sector by performing the conformal transformation already discussed in the previous section. After taking the conformal transformation, the resulting action in the Einstein frame reads

(28) |

where for a large field approximation

(29) |

This leads to the potential in the Einstein frame as

(30) |

With the large field approximation, we can derive the following slow-roll parameter in terms of the field as

(31) |

Inflation ends when such that

(32) |

In the large field limits, the number of e-foldings reads

(33) |

Using eqs.(8)-(10), it is straightforward to express the inflationary predictions in terms of the number of e-foldings, and find the Einstein frame parameters:

(34) |

Here we kept to the lowest order in to derive the spectral index and the tensor-to-scalar ratio. It is worthy to note here that the higher corrections of and are suppressed to these expressions. This model provides and for . Unfortunately, from Eq. (34), the predictions are in tension with the recent BICEP2 observations since a large value of cannot be produced in this model, see detailed discussions in the last section.

Glueball Inflation (GI)

In this section, we consider the work presented in [18] as the model motivated by a pure Yang-Mills theory. In this case, the inflaton is the lightest glueball field with . The reason why we want to continue analysing composite inflation by following such model resides from the fact that this theory features the archetype of any composite paradigm in flat space and consequently of models of composite inflation. The low-energy effective Lagrangian of the lightest glueball state can be found in [38, 39, 40].

It is worthy to note here that the theory we are using describes the ground state of pure Yang-Mills theory, and of course is not the simple theory. Another difference from the theory is that the form of the effective potential, before coupling to gravity, is completely fixed by the underlying gauge theory. The action in the Jordan frame in this case reads [18]

(35) |

which yields

(36) |

Imposing the conformal transformation, the resulting action in the Einstein frame reads

(37) |

which leads to the potential in the Einstein frame for a large field approximation as

(38) |

Here we have left the explicit dependence on the field possessing unity canonical dimension instead of using the canonically normalized new scalar field . However, we can express the resulting potential in terms of the field by substituting into the potential given in Eq. (38). With the large field approximation, we can derive the following slow-roll parameter

(39) |

Inflation ends when such that

(40) |

In the large field limits, the number of e-foldings reads

(41) |

Using eqs.(8)-(10), the inflationary predictions in terms of the number of e-foldings in the Einstein frame parameters read

(42) |

where we have assumed . We also have kept to the lowest order in to derive the spectral index and the tensor-to-scalar ratio and neglected the higher corrections of and suppressed to these expressions. The predictions of this model from Eq. (42) provide and for . Concretely, more detailed discussions will be added in the last section.

Super Yang-Mills Inflation (SYMI)

The application of the supersymmetric version of a pure Yang-Mills (SYM) has been investigated in various cornerstones. Here we consider for inflation the bosonic part of the Veneziano-Yankielowicz effective theory. For inflationary model building, the gluino-ball state in the super Yang-Mills theory is assigned as the inflaton with . The effective Lagrangian in supersymmetric gluodynamics was constructed in [41]. Having already stamped as the model of composite inflation, the authors of [19] took the scalar component part of the super-glueball and coupled it non-minimally to gravity. The action of the theory in term of the canonical-dimension field reads [19]

(43) |

with a number of colours and a -independent quantity. With the action given above, we find

(44) |

However, the gravity-scalar coupled sector can basically be diagonalised by imposing the conformal transformation. We find the resulting action in the Einstein frame as

(45) |

where

(46) |

where we have also left the explicit dependence on the field instead of using the canonically normalized new field . With the large field approximation, we can derive the following slow-roll parameter

(47) |

Inflation ends when such that

(48) |

Performing the similar approximations to those of the previous subsection, the number of e-foldings for this model in the large limit reads

(49) |

Using eqs.(8)-(10), the inflationary predictions, i.e. , in terms of the number of e-foldings in the Einstein frame read

(50) |

In order to derive these parameters, we have kept to the lowest order in and neglect the higher corrections of and . We find for this model from Eq. (50) that and for . In the last section, we will summarize our findings for this model.

Orientifold Inflation (OI)

The authors of [19] examined the supersymmetric low-energy effective action to study inflation driven by the gauge dynamics of SU(N) gauge theories adding one Dirac fermion in either the two-index antisymmetric or symmetric representation of the gauge group. Such theories are known as orientifold theories. Here the gluino field of supersymmetric gluodynamics is replaced by two Weyl fields which can be formed as one Dirac spinor. For investigating the inflationary scenario, the orientifold sector non-minimally coupled to gravity in the Jordan frame action is implemented in [19]. For this investigation, we write the action by using the real part of the field as [19]

(51) |

where , and hereafter we will keep only leading order in . With the action given above, we find

(52) |

As usually proceeded in the standard fashion, we impose the conformal transformation and then find the resulting action in the Einstein frame

(53) |

where

(54) |

with being a number of colours. Note that at large- the theory features certain super Yang-Mills properties, i.e. . With this limit, the transformed potential reduces to that of (46). With the large field limit, we can derive the following slow-roll parameter

(55) |

Inflation ends when such that

(56) |

where can be directly obtained from Eq. (48). According to this model, the number of e-foldings in the large limit is approximately given by

(57) |

Using eqs.(8)-(10), the predictions in the Einstein frame for and of this model can consequently written in terms of the number of e-foldings as

(58) |

Notice that for large the observables given above features the Super Yang-Mills inflation since .

## 5 Summary

In this last section we summarise our findings by which we confront the predictions of composite inflationary models with Planck and also very recent BICEP2 data and then plot our results on the constraint contours illustrated in Fig. (1). Moreover, recent study in the Jordan frame the primordial spectral index and tensor-to-scalar ratio predicted by composite models has been discussed in [32]. For this work, the following is our remarks.

The first model of composite inflation (MCI) is inspired by models of dynamical electroweak symmetry breaking with the Higgs-like potential. Here the Higgs sector of the SM is replaced by a new underlying four-dimensional gauge dynamics free from fundamental scalars and the Higgs field is therefore composite state. We discover for model that the predictions can satisfy the CL of the Planck data for . Moreover, the model predictions can be consistent with the Planck data up to the CL if we choose . However, we find that the MCI model is in tension with the recent BICEP2 data. This is so since the model predictions yield quite small values of both and if we consider large limit. For a small case, there have been discussed in [32]. In order to satisfy up to the CL of Planck data, we discover the number of e-foldings should not be greater than .

The underlying theory of the second paradigm (GI) features the archetype of any composite paradigm in flat space and consequently of models of composite inflation. The observables predicted by GI model lie inside the CL of the Planck for , whilst for this model can be consistent with the region of the Planck contours. Interestingly, we discover for this model that the predictions can satisfy the CL of the BICEP2 data for , and up to CL if we choose . In order to satisfy upto the CL of Planck and BICEP2 data, we discover the number of e-foldings should not be greater than and , respectively.

Finally, the SYMI model features the supersymmetric version of a pure Yang-Mills properties. The model provide an interesting class of inflationary model. Apparently, the SYMI predictions are significantly consistent with the Planck and BICEP2 constraints. We find that in order to satisfy at the CL of the BICEP2 data, the model prefers . The model can be consistent with the 68 CL for both Planck and simultaneously BICEP2 data for , and with CL for . For the OI model, the present of parameter yields small change of the predictions. However, the results from the OI model coincide with those of the SYMI for a large number of colours. In order to satisfy the and CL of the BICEP2 data, we discover the number of e-foldings should not be greater than and , respectively.

Another crucial consequence for the model of inflation is the (p)reheating mechanism. We anticipate to investigate this mechanism by following closely references [44, 45]. Rather interestingly, for instance, the author of [46] studied the consequent of inflation as seed of the present intergalactic magnetic field. However, the author claimed that the results after making a number of simplifying approximations should be considered to be preliminary. Therefore, it is very interesting to study the mechanism for generating an intergalactic magnetic field based on the composite inflationary manners.

Acknowledgment : The author thanks Arunee Shunava and Worakrit Thida for assisted graphic works. Remote conversation with the CP-Origins’s Team, F. Sannino, J. Joergensen and O. Svendsen, is also acknowledged. This work is fully supported by the “Research Fund for DPST Graduate with First Placement” under Grant No. 033/2557.

### References

- A. A. Starobinsky, “A New Type Of Isotropic Cosmological Models Without Singularity,” Phys. Lett. B 91, 99 (1980).
- K. Sato, Mon. Not. R. Astron. Soc. 195, 467 (1981).
- K. Sato, ”Cosmological baryon-number domain structure and the first order phase transition of a vacuum”, Phys. Lett. 99B, 66 (1981).
- D. Kazanas, “Dynamics Of The Universe And Spontaneous Symmetry Breaking,”Astrophys. J. 241 L59 (1980).
- A. H. Guth, “The Inflationary Universe: A Possible Solution To The Horizon And Flatness Problems,” Phys. Rev. D 23, 347 (1981).
- A. D. Linde,“A New Inflationary Universe Scenario: A Possible Solution of the Horizon, Flatness, Homogeneity, Isotropy and Primordial Monopole Problems,” Phys. Lett. B 108, 389 (1982)
- A. Albrecht and P. J. Steinhardt,“Cosmology for Grand Unified Theories with Radiatively Induced Symmetry Breaking,” Phys. Rev. Lett. 48, 1220 (1982).
- A. D. Linde,“Chaotic Inflation,” Phys. Lett. B 129, 177 (1983).
- K. Freese, J. A. Frieman and A. V. Olinto,“Natural inflation with pseudo - Nambu-Goldstone bosons,” Phys. Rev. Lett. 65, 3233 (1990).
- F. C. Adams, J. R. Bond, K. Freese, J. A. Frieman and A. V. Olinto,“Natural inflation: Particle physics models, power law spectra for large scale structure, and constraints from COBE,” Phys. Rev. D 47, 426 (1993) [hep-ph/9207245].
- J. R. Espinosa, A. Riotto and G. G. Ross, “D - term inflation in superstring theories,” Nucl. Phys. B 531, 461 (1998) [hep-ph/9804214].
- J. A. Casas, G. B. Gelmini and A. Riotto, “F term inflation in superstring theories,” Phys. Lett. B 459, 91 (1999) [hep-ph/9903492].
- M. Yamaguchi, “Supergravity based inflation models: a review,” Class. Quant. Grav. 28, 103001 (2011) [arXiv:1101.2488 [astro-ph.CO]].
- F. Farakos, A. Kehagias and A. Riotto, “On the Starobinsky Model of Inflation from Supergravity,” Nucl. Phys. B 876, 187 (2013) [arXiv:1307.1137].
- A. Linde, “Chaotic inflation in supergravity and cosmic string production,” Phys. Rev. D 88, 123503 (2013) [arXiv:1303.4435 [hep-th]].
- S. Ferrara, R. Kallosh, A. Linde and M. Porrati, “Minimal Supergravity Models of Inflation,” arXiv:1307.7696 [hep-th].
- P. Channuie, J. J. Joergensen and F. Sannino, “Minimal Composite Inflation,” JCAP 1105, 007 (2011) [arXiv:1102.2898 [hep-ph]].
- F. Bezrukov, P. Channuie, J. J. Joergensen and F. Sannino, “Composite Inflation Setup and Glueball Inflation,” Phys. Rev. D 86, 063513 (2012) [arXiv:1112.4054 [hep-ph]].
- P. Channuie, J. J. Jorgensen and F. Sannino, “Composite Inflation from Super Yang-Mills, Orientifold and One-Flavor QCD,” Phys. Rev. D 86, 125035 (2012) [arXiv:1209.6362 [hep-ph]].
- M. Cvetic, T. Hubsch, J. C. Pati and H. Stremnitzer,“A Natural Origin Of Inflation Within A Class Of Supersymmetric Preon Models,” Phys. Rev. D 40, 1311 (1989).
- S. D. Thomas, “Moduli inflation from dynamical supersymmetry breaking,” Phys. Lett. B 351, 424 (1995) [hep-th/9503113].
- J. Garcia-Bellido, “Dual inflation,” Phys. Lett. B 418, 252 (1998) [hep-th/9707059].
- R. Allahverdi, K. Enqvist, J. Garcia-Bellido and A. Mazumdar, “Gauge invariant MSSM inflaton,” Phys. Rev. Lett. 97, 191304 (2006) [hep-ph/0605035].
- K. Hamaguchi, K. -I. Izawa and H. Nakajima, “Supersymmetric Inflation of Dynamical Origin,” Phys. Lett. B 662, 208 (2008) [arXiv:0801.2204 [hep-ph]].
- N. Evans, J. French and K. -y. Kim, “Holography of a Composite Inflaton,” JHEP 1011, 145 (2010) [arXiv:1009.5678 [hep-th]].
- G. L. Alberghi and R. Casadio, “Non-Gaussianity from Compositeness,” arXiv:1010.4395 [astro-ph.CO].
- G. L. Alberghi, “CMB: A Look Inside the Inflaton,” arXiv:0902.4787 [gr-qc].
- J. Martin, C. Ringeval and V. Vennin, “Encyclopaedia Inflationaris,” arXiv:1303.3787 [astro-ph.CO].
- C. T. Hill and D. S. Salopek, “Calculable nonminimal coupling of composite scalar bosons to gravity,” Annals Phys. 213, 21 (1992).
- Y. Fujii and K. Maeda, “The Scalar-Tensor Theory of Gravitation”, Cambridge Monographs on Mathematical Physics, 2004 (first published in printed format in 2003).
- R. L. Arnowitt, S. Deser and C. W. Misner, “The Dynamics of general relativity,” Gen. Rel. Grav. 40, 1997 (2008) [gr-qc/0405109].
- P. Channuie and K. Karwan, “Composite Inflation confronts BICEP2 and PLANCK,” arXiv:1307.2880 [hep-ph].
- F. Sannino and K. Tuominen, “Orientifold theory dynamics and symmetry breaking,” Phys. Rev. D 71, 051901 (2005) [hep-ph/0405209].
- D. K. Hong, S. D. H. Hsu and F. Sannino, “Composite Higgs from higher representations,” Phys. Lett. B 597, 89 (2004) [hep-ph/0406200].
- D. D. Dietrich, F. Sannino and K. Tuominen, “Light composite Higgs and precision electroweak measurements on the Z resonance: An Update,” Phys. Rev. D 73, 037701 (2006) [hep-ph/0510217].
- D. D. Dietrich, F. Sannino and K. Tuominen, “Light composite Higgs from higher representations versus electroweak precision measurements: Predictions for CERN LHC,” Phys. Rev. D 72, 055001 (2005) [hep-ph/0505059].
- T. Kobayashi, M. Yamaguchi and J. ’i. Yokoyama, “Generalized G-inflation: Inflation with the most general second-order field equations,” Prog. Theor. Phys. 126, 511 (2011) [arXiv:1105.5723 [hep-th]].
- J. Schechter, “Effective Lagrangian with Two Color Singlet Gluon Fields,” Phys. Rev. D 21, 3393 (1980).
- A. A. Migdal and M. A. Shifman, “Dilaton Effective Lagrangian in Gluodynamics,” Phys. Lett. B 114, 445 (1982).
- J. M. Cornwall and A. Soni, “Couplings Of Low Lying Glueballs To Light Quarks, Gluons, And Hadrons,” Phys. Rev. D 29, 1424 (1984).
- G. Veneziano and S. Yankielowicz, “An Effective Lagrangian for the Pure N=1 Supersymmetric Yang-Mills Theory,” Phys. Lett. B 113, 231 (1982).
- P. A. R. Ade et al. [Planck Collaboration], “Planck 2013 results. XXII. Constraints on inflation,” arXiv:1303.5082 [astro-ph.CO].
- P. A. R. Ade et al. [BICEP2 Collaboration], “BICEP2 I: Detection Of B-mode Polarization at Degree Angular Scales,” arXiv:1403.3985 [astro-ph.CO].
- S. Tsujikawa, K. -i. Maeda and T. Torii, “Preheating with nonminimally coupled scalar fields in higher curvature inflation models,” Phys. Rev. D 60, 123505 (1999) [hep-ph/9906501].
- Y. Watanabe and E. Komatsu, “Reheating of the universe after inflation with f(phi)R gravity,” Phys. Rev. D 75, 061301 (2007) [gr-qc/0612120].
- B. Ratra, “Cosmological ’seed’ magnetic field from inflation,” Astrophys. J. 391, L1 (1992).