Constraining the general reheating phase in the attractor inflationary cosmology
Abstract
In this paper we constrain some aspects of the general postinflationary phase in the context of superconformal attractor models of inflation. In particular, we provide constraints on the duration of the reheating process, , and on the reheating temperature, , simulating possible and future results given by the nextgeneration of cosmological missions. Moreover, we stress what kinds of equationofstate parameter, , are favored for different scenarios. The analysis does not depend on the details of the reheating phase and it is performed assuming different measurements of the tensortoscalar ratio .
pacs:
Valid PACS appear hereI Introduction
The theory of inflation is the most promising process for the description of the early Universe 1 (). The detection of the inflationary background of gravitational waves (GW) could provide a further confirmation of the goodness of the inflationary paradigm and represents one of the most important goals of the current research activity. The Planck mission constraints at C.L. on a pivot scale of Mpc 2 (); 3 (). In the next future, we might detect inflationary gravitational waves by forthcoming polarization missions of the Cosmic Macrowave Background (CMB) (like LITEBIRD or COrE, see 4 (); 5 (); 6 (); 7 (); 8 () for details) or by gravitational waves experiments (ALIA and BBO missions, see 9 () for a review). However, the knowledge of the inflationary process could shed light on several aspects of the postinflationary phase, the reheating phase, responsible for the production of the observable entropy in the Universe 10 (); 11 (); 12 (); 13 (); 14 (); 15 (); 16 (); 17 (); 18 (); 19 (); 20 (); 21 (). In this work, we study what kinds of constraints we may have on some aspects of the postaccelerated phase assuming future detections of the tensortoscalar ratio in the range of . In particular, we assume the attractor models of inflation (see 22 (); 23 (); 24 (); 25 (); 26 (); 27 () for details) as the model for the early Universe and we provide constraints on the macro reheating variables, i.e., the number of foldings during the reheating phase, , and the reheating temperature, . In doing this, we outline the possible mean values of the equationofstate parameter (EoSp), , compatible with each given observation. Furthermore, we show how the constraints on these parameters are sensible to the uncertainty on the scalar spectral index, . Currently, PLANCK mission fixes (see 2 (); 3 ()) which provides a very large uncertainty on . The paper is organized in the following way. In Sec. II, we recall the basic concepts of inflation and discuss the reheating phase introducing the fundamental parameters, , and . In Sec. III, we reconstruct the main parameters of the reheating stage with respect to the mean value of the EoSp of the reheating fluid. In Sec. IV, we discuss the results and the implications of the work. In this paper we use the natural units of particle and cosmology .
Ii Inflation, Reheating and Potential reconstruction
Inflation and quantum fluctuations
The inflationary paradigm provides a fast accelerated expansion of the early Universe, on energy scale close to the one of the simmetry breaking of the grand unified interaction ( GeV). Mathematically, we have , where is the cosmic scale factor and is the number of foldings, the number of exponential growth before the end of inflation. The most popular scenario for inflation involves the simplest version of a scalar field (neutral, homogeneous, minimally coupled to gravity and canonically normalized) with an effective autointeraction potential sufficiently flat and characterized by a global minimum. Then, the earlycosmological action assumes the following generic form:
(1) 
where is the Ricci scalar, is the metric tensor, is its determinant and is the reduced Planck mass. The scalar field, is often called “inflaton”, while the function is the inflationary potential. The nature of the inflaton field is an open issue. However, it is common belief that could arise from the low energy limit of some more fundamental cosmological theory. The inflationary process starts with the scalar field that explores the flat ragion or “false vacuum” with a very slow dynamics: . Useful quantities to describe this phase are the so called slow rollparameters,
(2) 
During this stage we have the production of the cosmological perturbations related to the large scale structure and to the temperature fluctuations of the cosmic background radiation 28 (); 29 (); 30 (); 31 (); 32 (); 33 (); 34 (); 35 (); 36 (). The power spectrum and the spectral index for the generated scalar perturbations are respectively,
(3) 
and
(4) 
while the corresponding quantities for the tensor sector are
(5) 
and
(6) 
Finally, the tensortoscalar ratio amplitude is given by
(7) 
Note that all of these quantities are functions of the scalar field and are calculated at first order in the slow roll parameters at horizon crossing of the pivot scale . The observables and are the quantities useful to test the “degrees of compatibility” of a given model with the cosmological data. As broadly discussed in 37 (); 38 (), we can rewrite the observables in terms of the number of foldings , by using the relation
(8) 
where is the variation of the field between the horizon crossing of quantum mode and the end of inflation, while
are the parameters of the potential function.
In conclusion, we have relations of the form , .
The end of inflation: the reheating phase
The inflationary process stops when the field approaches the value . In general, one has . It is straightforward to show that the energy density at this epoch is given by
(9) 
In the following, the field quickly relaxes toward the real minimum of the potential and here it takes mass, as well as starts to oscillate and decay, producing the radiation characterizing the radiation dominated era of standard Big Bang Model. This phase is called the reheating phase and can be described by a single reheating fluid, with a “mean” EoSp, . The physics of reheating is extremely complicated and depends on many factors such as the dominant operator terms of the potential about the minimum and the details of the inflaton physics. The first and simplest proposed model for the description of the reheating era is sometimes called “elementary theory of reheating” 9 (); 10 (); 11 () although we can have more complicated effects like preheating and turbulent stages 13 (); 14 (); 15 (); 16 (); 17 (); 18 (); 19 () (a complete review can be found in 20 () while the first reheating constraints are given in 21 ()). The simplest model provides a first stage of coherent oscillations of the inflaton field about the minimum of a simple quadratic potential. The value of EoSp is fixed to be zero, . The oscillations produce a cold gas of inflatonparticles that decay to relativistic particles. Afterwards, a thermalization stage takes place in which the particles strongly interact with each other and reach the thermal equilibrium. The duration of the reheating phase is of the order of , where is the inflaton decay rate. However, independent on the details of the reheating mechanism we can introduce two macro quantities: the number of foldings during the reheating stage, and the temperature at the end of the process, . Indeed, exploiting the fact that the number of foldings before the end of inflation depends on the complete history of the Universe after the inflationary phase (and so on ), we can show that for any model holds (see references 39 (); 40 (); 41 (); 42 ()),
(10) 
Here, are the parameters of the model and are the known cosmological quantities (CMB photon temperature , Hubble rate ,…). The function is defined as
(11) 
where is the effective EoSp of the conserved “reheating fluid”, is the pivot scale ( Mpc, in our case), is the remaining number of foldings before the end of inflation when crosses the horizon, GeV is the current Hubble rate, GeV is the CMB temperature, is the energy density when leaves the horizon, is the energy density at the end of inflation and is the number of relativistic degrees of freedom at the onset of the radiation dominance. Now, the above equations allow to derive the duration of the reheating phase, once a model is chosen. Note that, when the function equals zero, the reheating occurs instantaneously providing the maximum value for . Eq.(10) tells us that the duration of the reheating era depends on the mean value of the EoSp, i.e, on the nature of the reheating fluid. In particular, we should note that if then if . Alternatively, if we need . Actually, we need for to exit from the inflationary phase. On the other side, we need to avoid close to the minimum of . Then, this second case appears unnatural in the context of quantum field theory although remain a possibility 39 (). In literature one can find different numerical studies on this question, such as 21 (). The parameter allows to calculate the “general reheating temperature” (GRT) that results in
(12) 
In the limit of , the instantaneous scenario for the reheating mechanism is recovered. In this case, the “instantaneous reheating temperature” (IRT) comes out to be
(13) 
Then, Eq.(12) shows an exponential modulation of the instantaneous case. As we can expect, the energy scale of the reheating gets lower as its duration becomes larger. Moreover, the final reheating temperature also depends on . Now, the inflationary observables depend on the and on the value of the EoSp by the inverting Eq.(10). In this sense (once a model is chosen and for fixed values of the cosmological parameters in Eq.(10)), different values of the EoSp provide different pairs . On the other hand, measurements of both and allow to reconstruct a given model (see 37 () for details), and from it and using the equations Eq.(10) and Eq.(12). In the following, we apply such a procedure in the context of the attractor inflationary models.
Iii The attractor models of inflation: an overview
The attractor models of inflation represent a very interesting class of inflationary models. This class can rise from different scenarios, but the most advanced version is actually realized in the supergravity landscape (22 (); 23 (); 24 (); 25 (); 26 (); 27 ()). In particular we have two classes of attractor potentials. The first one is the socalled Tmodel class for which the prototype potential is given by
(14) 
The second is the Emodel class defined by
(15) 
These classes of models are particular interesting especially by virtue of its capability to interpolate a broad range of predictions in the plane. Herein, we want to focus on the simplest version of the Eclass, where
(16) 
In general, is the model parameter defined as:
(17) 
In such a way, within supergravity context, the parameter is related to the curvature of the Kälher geometry associated with the inflaton field
(18) 
In the limit of the false vacuum is recovered as shown in Fig.(1).
The relation between the vacuum and the number of foldings is given by solving the equation of motion Eq.(8):
(19) 
Note that, when we apply the condition , we can simplify the solution as
(20) 
Using the condition , we can show the vacuum expectation value of the scalar field at the end of the accelerated phase, results in
(21) 
This value is crucial to determine once is fixed, inverting Eq.(19) by the Lambert function. The predicted scalar spectral index and the tensortoscalar ratio of the attractor models (in the limit of small ) are given by the following relations:
(22) 
The dependency of and with respect to the number of foldings is quite common: it is shared by a lot of important models such as the Tmodels themselves, the Starobinsky scenario (for ) 43 (); 44 (), the GoncharovLinde model 45 (), the Higgs inflation 46 () or the so called string moduli inflation 47 (); 48 (); 49 (); 50 (), because the shape of the potential during the inflationary stage provides . The plateau of the potential function is replaced by a pure quadratic shape if . However, we studied the local potential reconstruction (see 51 () for a review) of the Emodels with in 37 () which provides
(23) 
with
(24) 
Starting from these definitions we can now reconstruct the postaccelerated epoch.
Iv Constraints on the reheating phase in early attractor Universe
The reconstruction of the postinflationary phase is possible when a robust statistical information is available on the input parameters and . For this goal, we simulate values of and , randomly extracted from a multivariate Gaussian distribution of the form
(25) 
Here
(26) 
where
(27) 
The parameters and are the mean and rms values of the scalar index while and are the corresponding values for the tensortoscalar ratio. The parameter is the correlation coefficient. Let us start, considering and , compatibly with the current Planck bounds, and three different values for () with the same uncertainty . We also assume to take into account our ignorance about possible correlations among and . The vacuum expectation value of the scalar field at the end of inflation is reconstructed through (or ) by the Eq.(21). In Tab.I we report the results for the three chosen values of . Note that, these results represent the extremal values that assumes in the considered cases: in principle inflation can stops before, in general when and this induces larger values for . In Fig.(2) we show the distribution of the values of associated with a detection . The variable is useful to derive the energy density of the Universe when inflation stops, , by the relation Eq.(9). The number of foldings is reconstructed once and are known by Eq.(19). Alternatively, one can also use Eq.(20) if the condition is taken into account. In Tab.II we report the related Monte Carlo results. These kinds of  values on implies a great uncertainty on the number of foldings and especially, on the reheating temperature due to Eq.(10) and Eq.(12).
mean value   value  

0.001  0.781  
0.002  0.889  
0.003  0.950 
mean value   value  

0.001  63.7  
0.002  63.0  
0.003  62.5 
This piece of evidence suggests that improvements of the estimation on could be strongly reduce the uncertainty on the inflationary number folds . Nowadays, foreseen cosmological experiments aim to have an uncertainty or even better. In this paper we try to be optimistic and we set . As we will discuss in the last section, it is hard to get such a sensitivity from a single experiment but rather by combining future several experiments. We assume two possible scenarios for the scalar spectral index

,

,
and for both cases we consider two possible realizations for the tensortoscalarratio
with and . In the analysis we stress what kinds of EoSp are allowed for each cases ( larger or less then ). As we can see from Tab.III, in the first case there is an independence on the order of magnitude that the mean value of takes on. In all considered cases of , the function results to be positive and low values of the EoS, , are allowed. Once the parameter called results to be setted then, the mean value of increases as gets larger. Consequently, the final reheating temperature gets lower. When we increase the value of the scalar spectral index up to , we deal with two different situations. In the range of , we have so the reheating is characterized by a large mean value of the EoS, (Tab.IV). In this case, once results are fixed, gets lower as increases. Therefore, the reheating temperature becomes larger. In the range of , instead, we have and so the required values of the EoS turns to be . Consequently, we recover the same behavior seen when (see Tab.V).
V Discussion and conclusions
In this paper, we reconstructed some aspects of the postinflationary phase in the context of superconformal attractor (supergravity) models of inflation. In particular, we focused our study on constraining the duration of the reheating era, , and the reheat scale of the Universe, i.e. the reheating temperature, . The reconstruction of these two parameters is strongly related to the mean value of the EoSp of the reheating fluid, . In our study we have outlined the existence of transitions. As we pass from a mean value of the scalar index to in the range of , then we have a change in the sign of the function : from to . This corresponds to a change of the allowed EoS from to . A this point, moving toward larger values of the function turns to be progressively positive. The case of is not favored from the standpoint of QFT but it still remain a possibility and in this sense it represents a compelling challenge in the context of early fundamental physics, as outlined also in 40 (). Scenarios in which are characterized by a specific property. Once is fixed, the duration of the reheating gets larger as increases and the reheating temperature gets lower. The same happen if we fix with a varying . In the case in which we have the opposite behavior as summarized in the previous section. Anyway, we want to stress that these results are overall in agreement with the analytic studies done by Eshaghi et all as well as Ueno and Yamamoto 42 (). For example, a detecion of and () with a standard equation of state , implies a very low reheating temperature of the order of GeV. At the same time a detection of and again with a standard equation of state, implies a reheating temperature GeV. Therefore we need , in order to move toward prolonged reheating with low temperatures, tipically as outlined especially by Eshaghi et all 42 (). Low values of the tensor to scalar ratio with are not consistent with standard equation of state and show an opposite behaviour: we have low reheating temperature for low values of . Another important aspects of this work is in the resulting distributions for . Providing robust constraints on the duration of the reheating is fundamental for understanding how much important is the hypothetical preheating phase. To have an efficient preheating phase we need or more. In some cases, (see Tab.V for , for instance) the is tiny, so the possibility for a preheating seems to be not so significantly even considering the relative errors. On the other hand, the evidence of a prolonged reheating phases implies very low reheating temperatures. As well summarized by Dai, Kamionkowsi and Wang in 40 (), there could be tensor fluctuations produced during inflation which reenter in the microphysical or Hubble horizon during the cosmic reheating history. These fluctuations therefore were stretched on very small scales compared to the current Hubble radius . In particular, these GW should be locally detectable on scales of the order of the solar system, about folds below the CMB scales as suggested in 52 (); 53 (), for instance. Therefore, they could bring fundamental cosmological and particle physics information about the reheating phase. The robustness of the presented analysis is related to the power of forthcoming future experiments to improve the current datasets. The first important point is in the capability of efficiently exploring ranges of the tensortoscalar ratio below the PLANCK mission sensitivity. The second important point is in the capability of reducing the uncertainty on the scalar spectral index . As we saw in the previous section, our forecasting are quite sensitive to some optimal constraint on . In fact, in our simulations we assume that we will be able to constrain up to the order of whereas the actual limit given by PLANCK are of the order of . This improvement, looking at the nextgeneration experiments could be a reachable goal. Indeed, missions such as EPIC 4 () or CoRE (LiteCore) 5 () combined with results coming from experiments like PRISM 8 (), EUCLID 54 () or by the 21cm surveys 55 (); 56 () could able to reduce the total uncertainty on at the required level of our analysis.
Acknowledgements.
This work has been supported in part by the STaI “Uncovering Excellence” grant of the University of Rome “Tor Vergata”, CUP E82L15000300005. A.D.M would like to thank Nicola Bartolo and Sabino Matarrese for the useful discussion on the early version of the paper. P.C thanks “L’isola che non c’è S.r.l” for the support.References
Footnotes
 preprint: APS/123QED
References
 A. H. Guth, “The Inflationary Universe: A Possible Solution to the Horizon and Flatness Problems,” Phys. Rev. D 23 (1981) 347; A. D. Linde, “A New Inflationary Universe Scenario: A Possible Solution of the Horizon, Flatness, Homogeneity, Isotropy and Primordial Monopole Problems,” Phys. Lett. 108B (1982) 389; A. Albrecht and P. J. Steinhardt, “Cosmology for Grand Unified Theories with Radiatively Induced Symmetry Breaking,” Phys. Rev. Lett. 48 (1982) 1220; S. W. Hawking and I. G. Moss, “Supercooled Phase Transitions in the Very Early Universe,” Phys. Lett. 110B (1982) 35; A. D. Linde, “Chaotic Inflation,” Phys. Lett. 129B (1983) 177.
 P. A. R. Ade et al. [Planck Collaboration], “Planck 2015 results. XIII. Cosmological parameters,” Astron. Astrophys. 594 (2016) A13 [arXiv:1502.01589 [astroph.CO]].
 P. A. R. Ade et al. [Planck Collaboration], “Planck 2015 results. XX. Constraints on inflation,” Astron. Astrophys. 594 (2016) A20 [arXiv:1502.02114 [astroph.CO]].
 J. Bock et al. [EPIC Collaboration], “Study of the Experimental Probe of Inflationary Cosmology (EPIC)Intemediate Mission for NASA’s Einstein Inflation Probe,” arXiv:0906.1188 [astroph.CO].
 F. R. Bouchet et al. [COrE Collaboration], “COrE (Cosmic Origins Explorer) A White Paper,” arXiv:1102.2181 [astroph.CO]; F. Finelli et al. [CORE Collaboration], “Exploring Cosmic Origins with CORE: Inflation,” arXiv:1612.08270 [astroph.CO].
 A. Kogut et al., “The Primordial Inflation Explorer (PIXIE): A Nulling Polarimeter for Cosmic Microwave Background Observations,” JCAP 1107 (2011) 025 [arXiv:1105.2044 [astroph.CO]].
 T. Matsumura et al., “Mission design of LiteBIRD,” J. Low. Temp. Phys. 176 (2014) 733 [arXiv:1311.2847 [astroph.IM]].
 P. Andre et al. [PRISM Collaboration], “PRISM (Polarized Radiation Imaging and Spectroscopy Mission): A White Paper on the Ultimate Polarimetric SpectroImaging of the Microwave and FarInfrared Sky,” arXiv:1306.2259 [astroph.CO].
 J. Crowder and N. J. Cornish, “Beyond LISA: Exploring future gravitational wave missions,” Phys. Rev. D 72 (2005) 083005 [grqc/0506015].
 L. F. Abbott, E. Farhi and M. B. Wise, “Particle Production in the New Inflationary Cosmology,” Phys. Lett. 117B (1982) 29.
 A. D. Dolgov and A. D. Linde, “Baryon Asymmetry in Inflationary Universe,” Phys. Lett. 116B (1982) 329.
 A. Albrecht, P. J. Steinhardt, M. S. Turner and F. Wilczek, “Reheating an Inflationary Universe,” Phys. Rev. Lett. 48 (1982) 1437;
 L. Kofman, A. D. Linde and A. A. Starobinsky, “Reheating after inflation,” Phys. Rev. Lett. 73 (1994) 3195 [hepth/9405187];
 L. Kofman, “The Origin of Matter in the Universe: Reheating after Inflation”, (1996) arXiv:astroph/9605155
 L. Kofman, A. D. Linde and A. A. Starobinsky, “Towards the theory of reheating after inflation,” Phys. Rev. D 56 (1997) 3258 [hepph/9704452];
 B. R. Greene, T. Prokopec and T. G. Roos, “Inflaton decay and heavy particle production with negative coupling,” Phys. Rev. D 56 (1997) 6484 [hepph/9705357];
 R. Micha and I. I. Tkachev, “Relativistic turbulence: A Long way from preheating to equilibrium,” Phys. Rev. Lett. 90 (2003) 121301 [hepph/0210202].
 R. Micha and I. I. Tkachev, “Turbulent thermalization,” Phys. Rev. D 70 (2004) 043538 [hepph/0403101].
 D. I. Podolsky, G. N. Felder, L. Kofman and M. Peloso, “Equation of state and beginning of thermalization after preheating,” Phys. Rev. D 73 (2006) 023501 [hepph/0507096].
 R. Allahverdi, R. Brandenberger, F. Y. CyrRacine and A. Mazumdar, “Reheating in Inflationary Cosmology: Theory and Applications,” Ann. Rev. Nucl. Part. Sci. 60 (2010) 27 [arXiv:1001.2600 [hepth]].
 J. Martin and C. Ringeval, Phys. Rev. D 82 (2010) 023511 [arXiv:1004.5525 [astroph.CO]].
 R. Kallosh and A. Linde, “Universality Class in Conformal Inflation,” JCAP 1307 (2013) 002 [arXiv:1306.5220 [hepth]]; R. Kallosh and A. Linde, “Multifield Conformal Cosmological Attractors,” JCAP 1312 (2013) 006 [arXiv:1309.2015 [hepth]].
 S. Ferrara, R. Kallosh, A. Linde and M. Porrati, “Minimal Supergravity Models of Inflation,” Phys. Rev. D 88 (2013) no.8, 085038 [arXiv:1307.7696 [hepth]].
 R. Kallosh and A. Linde, “Superconformal generalizations of the Starobinsky model,” JCAP 1306 (2013) 028 [arXiv:1306.3214 [hepth]].
 R. Kallosh, A. Linde and D. Roest, “Superconformal Inflationary Attractors,” JHEP 1311 (2013) 198 [arXiv:1311.0472 [hepth]].
 M. Galante, R. Kallosh, A. Linde and D. Roest, “Unity of Cosmological Inflation Attractors,” Phys. Rev. Lett. 114 (2015) no.14, 141302 [arXiv:1412.3797 [hepth]].
 R. Kallosh and A. Linde, “Cosmological Attractors and Asymptotic Freedom of the Inflaton Field,” JCAP 1606 (2016) no.06, 047 [arXiv:1604.00444 [hepth]].
 A. A. Starobinsky, “Spectrum of relict gravitational radiation and the early state of the universe,” JETP Lett. 30 (1979) 682 [Pisma Zh. Eksp. Teor. Fiz. 30 (1979) 719]. V. F. Mukhanov and G. V. Chibisov, “Quantum Fluctuations and a Nonsingular Universe,” JETP Lett. 33 (1981) 532 [Pisma Zh. Eksp. Teor. Fiz. 33 (1981) 549].
 L. F. Abbott and M. B. Wise, “Constraints on Generalized Inflationary Cosmologies,” Nucl. Phys. B 244 (1984) 541.
 S. W. Hawking, “The Development of Irregularities in a Single Bubble Inflationary Universe,” Phys. Lett. 115B (1982) 295
 A. A. Starobinsky, “Dynamics of Phase Transition in the New Inflationary Universe Scenario and Generation of Perturbations,” Phys. Lett. 117B (1982) 175; A. D. Linde, “Scalar Field Fluctuations in Expanding Universe and the New Inflationary Universe Scenario,” Phys. Lett. 116B (1982) 335;
 S. W. Hawking and I. G. Moss, “Fluctuations in the Inflationary Universe,” Nucl. Phys. B 224 (1983) 180.
 J. M. Bardeen, P. J. Steinhardt and M. S. Turner, “Spontaneous Creation of Almost Scale  Free Density Perturbations in an Inflationary Universe,” Phys. Rev. D 28 (1983) 679;
 V. F. Mukhanov, “Gravitational Instability of the Universe Filled with a Scalar Field,” JETP Lett. 41 (1985) 493 [Pisma Zh. Eksp. Teor. Fiz. 41 (1985) 402].
 F. Lucchin and S. Matarrese, Phys. Rev. D 32 (1985) 1316.
 E. D. Stewart and D. H. Lyth, “A More accurate analytic calculation of the spectrum of cosmological perturbations produced during inflation,” Phys. Lett. B 302 (1993) 171 [grqc/9302019].
 A. Di Marco, P. Cabella and N. Vittorio, “Reconstruction of attractor supergravity models of inflation,” Phys. Rev. D 95 (2017) no.2, 023516 [arXiv:1703.06472 [astroph.CO]].
 J. Martin, C. Ringeval and V. Vennin, “EncyclopÃ¦dia Inflationaris,” Phys. Dark Univ. 56 (2014) 75 [arXiv:1303.3787 [astroph.CO]].
 A. R. Liddle and S. M. Leach, “How long before the end of inflation were observable perturbations produced?,” Phys. Rev. D 68 (2003) 103503 [astroph/0305263].
 L. Dai, M. Kamionkowski and J. Wang, “Reheating constraints to inflationary models,” Phys. Rev. Lett. 113 (2014) 041302 [arXiv:1404.6704 [astroph.CO]].
 J. B. Munoz and M. Kamionkowski, “EquationofState Parameter for Reheating,” Phys. Rev. D 91 (2015) no.4, 043521 [arXiv:1412.0656 [astroph.CO]].
 M. Eshaghi, M. Zarei, N. Riazi and A. Kiasatpour, “CMB and reheating constraints to attractor inflationary models,” Phys. Rev. D 93 (2016) no.12, 123517 [arXiv:1602.07914 [astroph.CO]]; Y. Ueno and K. Yamamoto, “Constraints on attractor inflation and reheating,” Phys. Rev. D 93 (2016) no.8, 083524 [arXiv:1602.07427 [astroph.CO]].
 A. A. Starobinsky, “A New Type of Isotropic Cosmological Models Without Singularity,” Phys. Lett. 91B (1980) 99.
 M. B. Mijic, M. S. Morris and W. M. Suen, “The Cosmology: Inflation Without a Phase Transition,” Phys. Rev. D 34 (1986) 2934.
 A. S. Goncharov and A. D. Linde, “Chaotic Inflation Of The Universe In Supergravity,” Sov. Phys. JETP 59 (1984) 930 [Zh. Eksp. Teor. Fiz. 86 (1984) 1594]; A. B. Goncharov and A. D. Linde, “Chaotic Inflation in Supergravity,” Phys. Lett. 139B (1984) 27.
 F. L. Bezrukov and M. Shaposhnikov, “The Standard Model Higgs boson as the inflaton,” Phys. Lett. B 659 (2008) 703 [arXiv:0710.3755 [hepth]].
 J. P. Conlon and F. Quevedo, “Kahler moduli inflation,” JHEP 0601 (2006) 146 [hepth/0509012].
 M. Cicoli, C. P. Burgess and F. Quevedo, “Fibre Inflation: Observable Gravity Waves from IIB String Compactifications,” JCAP 0903 (2009) 013 [arXiv:0808.0691 [hepth]].
 M. Cicoli, F. G. Pedro and G. Tasinato, “Polyinstanton Inflation,” JCAP 1112 (2011) 022 [arXiv:1110.6182 [hepth]].
 C. P. Burgess, M. Cicoli, S. de Alwis and F. Quevedo, “Robust Inflation from Fibrous Strings,” JCAP 1605 (2016) no.05, 032 [arXiv:1603.06789 [hepth]].
 J. E. Lidsey, A. R. Liddle, E. W. Kolb, E. J. Copeland, T. Barreiro and M. Abney, “Reconstructing the inflation potential : An overview,” Rev. Mod. Phys. 69 (1997) 373 [astroph/9508078].
 A. R. Liddle, “Can the gravitational wave background from inflation be detected locally?,” Phys. Rev. D 49 (1994) 3805 Erratum: [Phys. Rev. D 51 (1995) 4603] [grqc/9307036].
 R. BarKana, “Limits on direct detection of gravitational waves,” Phys. Rev. D 50 (1994) 1157 [astroph/9401050].
 L. Amendola et al., “Cosmology and Fundamental Physics with the Euclid Satellite,” arXiv:1606.00180 [astroph.CO].
 Y. Mao, M. Tegmark, M. McQuinn, M. Zaldarriaga and O. Zahn, “How accurately can 21 cm tomography constrain cosmology?,” Phys. Rev. D 78 (2008) 023529 [arXiv:0802.1710 [astroph]].
 A. Kosowsky, “Can we test inflationary expansion of the early universe?,” Physics 3, 103 (2010)