Horizon feedback inflation

We consider the effect of the Gibbons–Hawking radiation on the inflaton in the situation where it is coupled to a large number of spectator fields. We argue that this will lead to two important effects – a thermal contribution to the potential and a gradual change in parameters in the Lagrangian which results from thermodynamic and energy conservation arguments. We present a scenario of hilltop inflation where the field starts trapped at the origin before slowly experiencing a phase transition during which the field extremely slowly moves towards its zero temperature expectation value. We show that it is possible to obtain enough e-folds of expansion as well as the correct spectrum of perturbations without hugely fine-tuned parameters in the potential (albeit with many spectator fields). We also comment on how initial conditions for inflation can arise naturally in this situation.


Introduction
Cosmological inflation is the leading paradigm which explains the horizon, flatness and defect problem of the a e-mail: malcolm.fairbairn@kcl.ac.uk b e-mail: t.markkanen@imperial.ac.uk c e-mail: david.rodriguez@kcl.ac.uk extremely successful FLRW hot big bang model as well as explaining the source of the initial density perturbations observed to exist in the CMB [1]. Furthermore the exponential expansion seems to have an elegant explanation in field theory as being sourced by the potential energy of a field which rolls only slowly down to its minimum, its kinetic energy being red-shifted by this rapid expansion [2][3][4]. It is also very difficult to think up alternatives to inflation which are natural, even with significant modifications of general relativity and those which do exist often create the wrong spectrum of perturbations [5].
Unfortunately, there are several problems with the standard inflationary scenario. The most difficult problem is (arguably) the fact that in order for inflation to start in the first place, one needs to find a Hubble patch in the early Universe across the entirety of which the potential energy dominates the kinetic energy of the field [6]. Another way of putting this is that in order to solve the horizon problem, one creates another one at earlier times. Another problem is that in order for the kinetic energy to be redshifted by the expansion for long enough to explain the horizon and flatness problems, either the inflaton field has to be transplanckian during inflation or the shape of the potential below the Planck scale has to be extremely flat, in other words the extent in the scalar field direction has to be much larger in GeV than its height in the energy direction.
Cosmological inflation usually requires a period of de-Sitter or at least quasi de-Sitter expansion in order to obtain the many e-folds required to solve the horizon problem. The Gibbons-Hawking temperature associated with the cosmological horizon in de Sitter space [7] T H = H 2π , (1.1) plays a central role in inflation as it acts as the source of quantum fluctuations in the inflaton field which source density perturbations. This thermal radiation is analogous to the thermal population observed by an accelerating observer i.e.
Unruh radiation [8] and closely resembles the Hawking radiation which surrounds a black hole [9]. In the case of a black hole, the thermal Hawking radiation escapes to infinity and energy is conserved only by postulating that the mass of the black hole is correspondingly released -a hypothesis which cannot be proved in semi-classical quantum gravity without including back-reaction but which does fit coherently into the theory of black hole thermodynamics [10]. Black hole evaporation has an interpretation of resulting from a flux of particles with negative mass going into the black hole. A similar physical interpretation for de Sitter radiation is possible, where the Cosmological Constant is reduced from the addition of negative vacuum (zero-point) energy to the overall energy density of de Sitter [11]. The author of [8] has famously said 'you could cook your steak by accelerating it' in regards to the physicalness of Unruh radiation [12]. In light of the analogy between the Unruh and de Sitter radiation, assigning similarly physical characteristics to the radiation associated with the de Sitter horizon motivates its inclusion in the Einstein equations [13,14]. In a first principle approach this can be achieved in a prescription where the energy-momentum tensor sourcing gravity is coarse grained to include only the observable degrees of freedom, i.e. those inside the cosmological horizon [11,15,16].
Including a source of continuous particle production in de Sitter space is well-known to bring about a qualitative change in the system's behaviour: when given enough time exponential expansion will cease implying that de Sitter space with particle production is not stable. In addition to thermal radiation from the horizon a de Sitter instability has been speculated to result from a variety other underlying mechanisms ranging from quantum gravity and infrared divergences in propagators to environment induced decoherence and the second law of thermodynamics [11,[36][37][38][39].
A de Sitter instability unavoidably leads to a dynamical cosmological constant and multiple studies have looked at the possibility that particle production could gradually reduce the value of the cosmological constant from a phenomenological point of view [40][41][42][43][44][45][46], see [47,48] for observational investigations. It has been argued that since particle production from gravitational fields seems to be an irreversible process, this reduction is inevitable from a thermodynamical perspective [33,[49][50][51][52]. Normally this effect is extremely small and would not really be a practical way of ensuring a small cosmological constant. In particular the rate of decrease would not get rid of the cosmological constant fast enough in the late Universe to explain today's cosmic acceleration.
It is interesting to consider the situation where the cosmological constant is not identified as a geometric term in the Einstein field equation, but is rather the energy density of a field which is located at a stable point in its potential where dV /dϕ = 0. In this situation the same logic would imply that terms in the Lagrangian that set the scale of the potential would also decay over time due to the gravitational production of particles. For example for the potential and for a field resting at the metastable origin ϕ = 0, the particle production associated with the space-time curvature arising from the non-zero potential energy at the origin should have the effect of making ϕ 0 decrease with time.
Another implication of the Gibbons-Hawking temperature would be the possibility that it might create nonnegligible finite temperature contributions to the potential energy of scalar fields which are evolving in the quasi de Sitter background.
In thermal inflation, a thermal sector with the equation of state of radiation and coupled to the inflaton keeps the field trapped at the origin for a few e-folds until the thermal radiation has been redshifted away, decreasing its temperature [53]. If however the origin of the thermal radiation is the Gibbons-Hawking temperature associated with the horizon then its temperature will be constant and will not decrease over time if there is no corresponding change in the vacuum energy. It is clear however that in this situation only a radiation bath with a large number of degrees of freedom, all of which are coupled to the scalar field, would be able to keep the field trapped at the origin. Also once located and stablised there due to this horizon temperature, the field would essentially be stuck at the origin for all time, leading to a graceful exit problem. A constant T H would also violate the continuity equation [14].
If however ϕ 0 changes over time due to the production of radiation, one can imagine a situation where ϕ 0 and consequently the energy density at the origin, the rate of expansion of the Universe and the temperature of the thermal radiation all decrease with time. The stable point at the origin then eventually becomes tachyonic, and field becomes free to roll away from ϕ = 0, setting the initial conditions for inflation (see also [54] for a similar idea in a different context). We argue that this could happen and that thermal effects would subsequently slow the phase transition sufficiently such that successful inflation can take place.
In Sect. 2 we will go through the equations of this scenario in more detail and study the dynamics of the field and how it might produce enough e-folds of expansion. Then in Sect. 3 we will study the inflationary predictions, namely the perturbations and the spectral tilt, as well as comparing our analytic estimates to a numerical treatment.

Phase transition with decaying vacuum energy
Let us examine the situation with the potential (1.2) in the presence of a thermal sector characterized by the Gibbons-Haking temperature (1.1) in more detail. We start by writing down the Friedmann equations Our model will consist of a scalar field ϕ and importantly N conformal fields that couple to ϕ and are in thermal equilibrium with the horizon as described by (1.1). This will lead to and additional temperature component for the energy and pressure densities ρ =φ and in addition to (1.2) the potential contains a contribution from the thermalised conformal fields We will throughout work in the approximation where the energy density of the thermal component is subdominant to that of the potential and in particular to the vacuum energy piece This condition will turn out to be easily satisfied for a large parameter range i.e. 1 from which it is quite apparent that a thermal sector with the Gibbons-Hawking temperature is inconsistent with strict de Sitter space, but will lead toḢ < 0. Furthermore, it is easy to see that the continuity equatioṅ can only be satisfied if ϕ 0 is not strictly constant but still provides the source for the continuous particle creation required for maintaining T H in the conformal fields despite the dilution from the expansion of space. This one may understand from the situation when ϕ = 0 giving ρ ∼ ρ , but ρ + p ∼ T 4 H . However as we will see we will not be relying on this feature of the current scenario in the current work other than to set initial conditions, the dynamics of ϕ 0 will be irrelevant by the time we come to calculate observables. Also we have deemed more natural to keep the usual interpretation of a strict coupling constant for the other constants of our theory, λ and g, and only allow ϕ 0 change with time. Since the continuity Eq. (2.6) provides only one constraint in principle one could consistently allow also the other parameters to vary, at least from a purely phenomenological point of view.
From (2.5) we may conclude that the change in H due to the thermal sector is very gradual: the first Hubble slow roll parameter for ϕ = 0 reads and if initially N (720π 2 M P /H init ) 2 then 1 is clearly satisfied. For the parameters in the example situation that we will present later 10 −7 . The evolution of the Hubble expansion rate as a function of the number of e-folds is then given by where H SB is the value of the Hubble expansion at symmetry breaking. (The number of e-folds N e are obtained simply by integrating d N e = Hdt.) However, the Friedman equations also show for ϕ = 0 that since ρ = 3H 2 M 2 P ∼ (λ/4)ϕ 4 0 , ϕ 0 will decrease slower than T H . Therefore, the effective mass squared for ϕ will eventually cross over to negative values, even if initially Simply put, at some point the system will undergo a phase transition. This phase transition in contrast to [53] is extremely gradual. As we will show, it can take several thousands of e-folds to complete. This is due to the special nature of the thermal radiation as given by T H : the thermal bath is continuously replenished by the decay of ϕ 0 and hence does not dilute in the usual fashion.
After the phase transition ϕ will start rolling to the new minimum and plays the role of the inflaton in the usual sense. The evolution of ϕ can be characterized with (2.9) to consist of three regions, m 2 0, m 2 ∼ 0 and m 2 0.
2.1 m 2 0: gradual decay of ρ When the effective mass of the field (2.9) is very large and positive the minimum of the potential is at ϕ = 0 and supposing the field is in the minimum (2.6) leads (in the approximation that neglectsḢ ) to Fig. 1 In the unbroken phase the field is at rest at ϕ = 0 with practically no quantum fluctuations, φ2 ∼ 0. The suppression of the quantum fluctuations is caused by the large effective mass making the field heavy with respect to the background curvature which is consistent with the requiremenṫ as given by the continuity equation (2.6). As discussed after (2.5) for most cases the change in the Hubble is very gradual which justifies neglecting the derivatives of H from the right side of the Einstein equation. Also interestingly the late time limit of (2.10) is independent of the initial Hubble rate indicating that the case m 2 0 exhibits late time attractor evolution which is independent of initial conditions: after a sufficiently long time with m 2 0 the system will always relax to a configuration with H ∼ (M 2 P /t) 1/3 and ϕ =φ = 0. For the remaining analysis we will choose the attractor configuration as our initial condition.
What is apparent from this section is that when m 0 the system has fairly unremarkable behaviour: the field sits put in its vacuum state and the vacuum energy ρ gradually decays. The quantum fluctuations around the mean we will denote asφ. When the field is heavy m 0 they are similarly suppressed The case m 0 is illustrated in Fig. 1. As discussed after (2.9) eventually the effective mass will vanish and the system will undergo a phase transition leading to interesting dynamics for ϕ.

m 2 ∼ 0: large quantum fluctuations
The phase transition happens when the effective mass of the field (2.9) vanishes at the origin. The value of the Hubble parameter at and soon after the symmetry breaking transition is approximately given by the first Friedmann equation from (2.1) as Close to symmetry breaking the classical (mean) field ϕ is almost stationary, but the quantum fluctuations φ 2 very large since the field is effectively light m ∼ 0 However, importantly, T H and hence m change very gradually and therefore immediately after symmetry breaking we expect a long epoch of very large quantum fluctuations before the rolling becomes dominant as illustrated in Fig. 2. The behaviour and magnitude of these fluctuations may be analytically solved via the stochastic formalism [55,56]. This epoch of large quantum fluctuations is quite important in our model since as we will further discuss in Sect. 2.3, when the T H has dropped enough the classical rolling will take over for which the initial condition will be dynamically set by the period of quantum jumping. We note that during this epoch where m 2 ∼ 0 we would expect the large density perturbations to lead to the formation of domain walls, but since we expect to obtain many more than 60 e-folds of inflation after this period, we expect these evils to be swept outside the horizon.
In order to obtain a qualitative picture of the system's behavior as the representative value for the classical dynamics one may take the square root of the variance ϕ 2 ≡ φ 2 .
Here we settle for solving this expectation value by using the Langevin equation in the Hartree approximation relegating a more complete discussion to Sect. 3.1. The relevant equation (2.14) During this epoch there is little dynamics in φ 2 and the two relevant limiting cases from the above are 15) which are separated by the threshold i.e. for −m 2 m 2 t (−m 2 m 2 t ) the results coincide with those on the left (right) side of (2.15).
The solutions (2.15) are approached only at the saturated limit when the system has been given enough time to equilibrate. The exact time this takes depends on the parameters of the potential. It is known that for a quartic theory the equilibration time scale in terms of e-folds is given by 1/ √ λ [57,58]. If the time scale for the change in the Hubble rate in terms of e-folds or −1 is much longer than the equilibration time our approximation of using the results at the saturated limit is valid. With (2.4), (2.9) at m ≈ 0 and (2.7) we can write the condition −1 λ −1/2 as which again is easy to satisfy for large N . Also, in our model the scale of symmetry breaking (2.13) can be tuned by choosing λ,N and g in our potential (2.3) with smaller scales corresponding to a slower dynamics as is evident from (2.7). Hence the condition in (2.17) can also be understood to imply the freedom to choose H SB low enough such that the saturated expressions in (2.15) are a good approximation.
In the absence of the horizon entropy thermalising a large number of fields coupled to the scalar field, this kind of symmetry breaking would not lead to good inflationary initial conditions since the Kibble mechanism would lead to large fluctuations from horizon to horizon with large field gradients that would prevent inflation from starting in the first place. In this scenario, the thermal corrections to the potential prevent the kinetic energy dominated regions from running straight down to |ϕ| = ϕ 0 . Eventually therefore as the thermal dissipation continues to gently facilitate the slow phase transition and the field's classical position and motion starts to dominate over its quantum fluctuations, a region should emerge somewhere where the field is coherent enough across several horizons such that the well known difficulties obtaining inflationary initial conditions are overcome. Fig. 3 The onset of classical rolling occurs when the steepness of the potential has increased enough. The initial value for the field is determined by the preceding random-walking epoch when m ∼ 0 2.3 m 2 0: classical rolling Strictly speaking one can only talk about a classical "rolling" once the minimum is further away from the origin than the stochastic vacuum expectation value that the field ϕ 2 ≡ φ 2 would obtain from the quantum fluctuations. As a first approximation one can say that the field starts rolling when in one e-fold the size of a single quantum jump (H/2π ) is smaller that the distance traveled due to the classical rolling which is the opposite to the usual condition for eternal inflation [59]. More accurately we can use (2.14) for the dynamics of φ 2 . As shown in the previous section, at the time of symmetry breaking (m = 0), the variance of the field is approximately 2 λ H φ 2 2 = H 3 4π 2 . At times soon after symmetry breaking the mass term remains negligible m 2 φ 2 3λ φ 2 2 and the field will remain constant at the value acquired at symmetry breaking until the mass is relevant, i.e. −m 2 − 3λ φ 2 = 0. At this point, the potential will be steep enough for classical slow roll to start and this classical motion will dominate over the quantum fluctuations (see Fig. 3). Note that this condition agrees with the estimate (2.18) since d dt φ 2 = H 3 4π 2 . Once the classical rolling starts the mass is given by and the value of the field and the speed are Once classical rolling has been triggered, our model gives rise to the usual slowly rolling inflation. We emphasize that the initial conditions for inflation are set dynamically by the large quantum fluctuations prior to classical rolling and hence are not free parameters. Similarly, the start of slow roll is triggered dynamically once the potential has acquired sufficient steepness and occurs for a wide range of values for ϕ 0 , in particular also for ϕ 0 /H 1. Finally, we remind the reader that the neat attractor behaviour of the solutions prior to symmetry breaking (2.10) also exhibit an independence from initial conditions. For these reasons we can conclude that the model presented here successfully evades all the usual fine-tuning issues of inflationary models.

Inflationary predictions
To have a successful inflationary scenario, we need enough e-folds (around 50-60) and we need to obtain the right perturbations P R and spectral tilt n s 50-60 e-folds before the end of inflation. As shown in the previous section, more than 60 e-folds of expansion can be obtained very easily in this scenario -the scale of inflation H can be set to be much smaller than the Planck mass so 1 continues for many e-folds before H differs significantly from the value at symmetry breaking H SB , meaning that we get enough inflation. Here we show how the perturbations are generated once the field is classically rolling to the minimum. The perturbations in the spatially flat gauge ( = 0) are defined as [60] which we note will in our model lead to the usual expression encountered in models of warm inflation [61]. The power spectrum takes the form where from now on we set ϕ 2 ≡ φ2 . In contrast to single field slow roll inflation, in the current setting the perturbations and the tilt are given by where we used the slow roll conditionφ = −V /3H , and dropped the subdominant kinetic pieceφ 2 in (2.5). Furthermore, we note that is negligible in the calculation of the tilt.
The difference with the usual single field slow roll inflation comes from the contribution of the thermal sector to the dynamical Friedman equation (2.5), which both changes the vacuum energy decay and modifies the usual calculation for the spectrum becauseḢ ∝ H 4 . Because of this, the perturbations follow the inverse of the usual P R ∝ H 4 /φ 2 behaviour. This means that the usual expression for the tilt n s − 1 ∝ η is not true and we get n s − 1 ∝ −η, where η = M 2 P V V is the usual slow roll parameter. So in order to obtain a red spectrum the inflaton needs to be evolving in the regime where the potential is convex (V > 0) i.e. beyond the inflection point at V = 0, contrary to the usual small field case.
There are a variety of different combinations of parameters that can give rise to the correct inflationary perturbations, however in what follows we will present a situation where the 50-60 e-folds of inflation we are interested in starts very soon after the symmetry breaking occurs. In this situation, the value of η is naturally small enough to give us the right tilt because soon after symmetry breaking V = 0 by definition. There is therefore no need to fine tune our tree-level parameters to give us a flat potential near the origin effectively bypassing the issue that the natural tree-level values are argued to give rise to η 1 [62].
Our system has 4 free parameters (N , λ, ϕ 0 and g). However for the sake of clarity we choose here the final e-folds of inflation to occur soon after symmetry breaking and hence we can approximate the value of ϕ 0 during inflation by which can be derived by making use of m SB = 0 in (2.9) and (2.13). With the above we effectively reduce the degrees of freedom from four to three. Also, if the observable part of the inflationary spectrum is going to take place soon after symmetry breaking then as soon as the classical rolling starts we want the field to be in the red tilt regime. For this to occur we need to ensure that d dt λ at the inflection point, otherwise the minimum would move faster away from the origin than the field. If this were to occur, the spectral tilt would be blue (after symmetry breaking). By making use of very similar steps that led to (3.4) and the Hartree approximation for φ 2 from (2.14) one may show this condition to lead to the order of magnitude constraint 1 g 2 < 15 4π 2 . So from now on, we will set for simplicity g = 1. However we also emphasize that a red tilt can be obtained a long time after symmetry breaking, so g = 1 is only a choice for the forthcoming calculation for the estimates of the perturbations and the tilt, and none of the above derivations relies on a specific value of g.
After this choice, we effectively have 2 degrees of freedom left (N , λ) that will determine the perturbations and the spectral tilt soon after the classical roll to the minimum starts. It proves convenient to define a new parameter, α = ϕ/ −m 2 /3λ, where α = 1 means that the field is at the inflection point and α = √ 3 means that the field is at the minimum of the potential. We can then approximate the perturbations and the tilt from (3.3) and with the help of equations from Sect. 2.3 by It is interesting to note that the order of magnitude of the spectral tilt does not depend on N while the magnitude of the perturbations depends on both (λ and N ). A smaller value of λ will make our spectrum become more scale invariant and will allow us to reduce the number of conformal fields that need to be present to give us the right spectrum. Hence we can shift the fine tuning between a small value of λ and a large number of spectator fields.
Since in this model inflation takes place between the inflection point and the new minimum emerging due to broken symmetry, the field excursion of ϕ scales as ϕ ∼ −m 2 /λ. Furthermore, since for our choice g = 1 this takes place soon after symmetry breaking where the mass parameter can be estimated as −m 2 ∼ √ λH 2 , which with the Hubble rate from (2.13) and (3.4) indicates that for a large number of spectator fields the field excursions are sub-Planckian. For the parameters agreeing with observations this turns out to be true. This indicates that despite begin a small field model of inflation, the initial conditions are not fine-tuned but arise naturally as the favoured attractor solutions unlike what is usually encountered [63].
The approximations and estimates set out here agree with the numerical solution that we will turn to now.

Numerical solution
In this section, we look at the numerical solution corresponding to the parameter choices we have made above.
Before symmetry breaking the field lies at the origin ϕ = 0 until ϕ 0 drops sufficiently for the potential at the origin to become tachyonic. Classically of course, the field would then not move anywhere because dV /dϕ = 0 at the origin, and if we were to introduce a small perturbation away from ϕ = 0 as an initial condition, then classically our final solution would depend strongly on this perturbation. To get around this problem we evolve the field using the Langevin equation which takes into account the stochastic quantum fluctuations the field receives where ξ is a Gaussian white noise with zero mean and unit variance.
For each combination of the remaining two free parameters (N , λ) we perform 10 4 numerical realisations of the We know that to obtain a red spectrum (n s < 1) the field must lie between the inflection point and the minimum. Figure 4 shows how the variance of the field φ 2 evolves with respect to the number of e-folds of expansion for the parameters λ = 10 −2 , N = 10 5.5 , and g = 1 (making ϕ 0 = 10 −1.41 M P ). Note that the initial epoch where we set the horizontal axis equal to zero is chosen arbitrarily. We also plot in the same figure the variance of the field that one would expect using the Starobinsky-Yokoyama prescription for the variance in a perfect de Sitter spacetime corresponding to a vacuum energy equal to our evolving vacuum energy [55], assuming instant equilibration for the probability distribution. The field expectation value lags behind this estimator, showing the importance of solving the Langevin equation. The field approximately remains constant at this value until the minimum has dropped enough to make the potential steep and for it to possess an inflection point. We also plot the position of the minimum as it moves out towards its zero temperature value and the position of the inflection point, showing that we remain on the good side to obtain a red spectrum. Note that in a period of time during inflation corresponding to 60 e-folds of expansion, the variation of ϕ 0 is negligible relative to the variation in ϕ so we assume that it is constant for calculational simplicity. Figure 5 shows the average evolution of the magnitude of perturbations for the same parameters as a function of the spectral tilt. We can see that these parameters can give rise to the correct combination of amplitude and spectrum for  [64,65].
The power spectrum of tensor fluctuations is given by the energy scale of inflation ( for example see equation 216 in [60]), therefore for the parameters chosen, the tensor to scalar ratio that we get is r = 0.17. This value is still within the 95% CL from the Planck data [1] although lower values can be obtained by choosing a different set of parameters, since the scale of symmetry breaking can be easily tuned, Eq. (2.13).
In order to end inflation, we assume the temperature corresponding to the expansion rate during inflation falls below the mass of the particles which form the thermal bath affecting the potential for the scalar field. We set this by hand to give us 60 e-folds after the epoch corresponding to the good values of P R and n s .

Discussion/Conclusion
In this work we have presented a possible new mechanism for inflation the early Universe. We have argued that if one takes the Gibbons-Hawking temperature associated with the horizon of de Sitter space seriously, one is lead to a couple of conclusions that can affect the evolution of quantum fields in the early Universe significantly. In particular we have argued that if there are enough fields coupled to the scalar field which takes the role of the inflaton, their Horizon induced temperature leads to thermal corrections to the potential which can affect its expectation value. Since these thermal fluctuations do not redshift as rapidly as normal radiation, this effect can last for a significant number of e-folds of expansion.
To summarise the mechanism -We consider a real scalar field with a Z 2 symmetric potential with minima at ±ϕ 0 . If the field starts at the origin, the non-zero energy density leads to de Sitter expansion and a cosmological horizon with an entropy. The field is trapped at the origin due to its coupling to a large number of conformal spectator fields which have a non-zero temperature associated with this horizon entropy. The height of the potential at the origin decays slowly as a result of the back-reaction of the thermal radiation on the parameter ϕ 0 in the Lagrangian. The expansion rate therefore decreases as does the temperature T H of the thermal radiation until the mass at the origin becomes tachyonic, at which point a phase transition occurs and the field rolls away from the origin towards its zero temperature minimum at |ϕ| = ϕ 0 . This phase transition however occurs extremely slowly due to the same finite temperature corrections to the potential. During this period of rolling we are able to obtain not only the correct number of e-folds but also the correct perturbations and spectral tilt as measured in the CMB. This setup has two attractive features:-• Normally, for a field to act as an inflaton, it must have a super Planckian expectation value and/or very finely tuned parameters. The mechanism we outline in this paper enables a potential which would otherwise not be flat enough to give rise to enough e-folds of inflation to do so due to the thermal effects of multiple spectator fields • The initial conditions for inflation may arise naturally as attractor solutions due to a slowly occurring phase transition. This happens despite the fact that inflation occurs with sub-Planckian values for the inflaton. This is the result of the parameters in the Lagrangian changing due to the back-reaction of the thermal radiation.
It is clear that this scenario is not necessarily a panacea for all of the problems of inflation. We need to assume that parameters in the Lagrangian decay over time due to the back reaction of Hawking radiation at the horizon. While there are many respected physicists who believe that such behaviour is probable and perhaps necessary, it is clear that further theoretical investigation and debate is required to put such speculation on a stronger footing. We also require quite a large number of fields which are coupled to the inflaton in order to obtain the thermal braking required to obtain enough e-folds of inflation, in the example we have put forward here about 10 5.5 . The masses of the fields need to be chosen such that inflation ends 50-60 e-folds after the epoch where the good perturbations and spectral index are set. In return for this cost, we obtain a theory of inflation which doesn't require transplackian field excursions, does not require extremely small parameters in the Lagrangian (we use a value of λ ∼ 0.01) and which may naturally explain the initial conditions for inflation. Finally, it seems quite challenging to understand how the difference between this scenario and normal inflation could be distinguished experimentally.