Relaxation of the Cosmological Constant
RRelaxation of the Cosmological Constant
Peter W. Graham, David E. Kaplan, and Surjeet Rajendran
2, 3 Stanford Institute for Theoretical Physics, Department of Physics, Stanford University, Stanford, CA 94305 Department of Physics & Astronomy, The Johns Hopkins University, Baltimore, MD 21218 Berkeley Center for Theoretical Physics, Department of Physics, University of California, Berkeley, CA 94720
We present a model that naturally tunes a large positive cosmological constant to a small cos-mological constant. A slowly rolling scalar field decreases the cosmological constant to a smallnegative value, causing the universe to contract, thus reheating it. An expanding universe with asmall positive cosmological constant can be obtained, respectively, by coupling this solution to anymodel of a cosmological bounce and coupling the scalar field to a sector that undergoes a technicallynatural phase transition at the meV scale. A robust prediction of this model is a rolling scalar fieldtoday with some coupling to the standard model. This can potentially be experimentally probed ina variety of cosmological and terrestrial experiments, such as probes of the equation of state of darkenergy, birefringence in the cosmic microwave background and terrestrial tests of Lorentz violation.
I. INTRODUCTION
The observed accelerated expansion of the Universe is well described by the existence of a small cosmologicalconstant. However, quantum corrections to this quantity are much larger than the observed value. One mighthope that the mysteries of quantum gravity hold the solution, but dangerous contributions come from very well-known physics at scales where spacetime curvature is weak (for example, finite corrections to vacuum energy fromthe electron mass). In this regard, the problem can be seen as one of fine-tuning, where contributions, known andunknown, conspire to cancel to generate the small value detected today.One approach to this puzzle is the introduction of an exponentially large number of universes, in which the vacuumenergy appears to be a random variable taking on different values in each. Anthropic selection then determines whichuniverse we are likely to appear in, based on the existence of structure or other arguments, and the assumption thata number of other parameters of our universe (such as the baryon-to-photon ratio, the dark matter abundance andthe value of the primordial density fluctuations) are the same over this exponentially large number of universes.A more natural solution could come from the dynamical relaxation of the cosmological constant in the early universevia a slowly rolling field in a potential. Indeed, such a model was attempted by Abbott [1] and others (e.g. [2]), anda similar model was successfully implemented in a solution to the gauge hierarchy problem [3]. Attempts to solvethe cosmological constant (CC) problem fell short as they invariably result in an empty universe. This is a robustproblem with relaxation: the CC can only be sensed through gravity, but gravity is universal and thus couples to thetotal space-time curvature. One way to solve this empty universe problem is to make the universe undergo a bounceafter the relaxation of the CC. This framework also solves many other thorny issues that confront solutions of the CCproblem such as the problem of cosmological phase transitions affecting the CC after relaxation [14]. Motivated bythese considerations, in [14], we discussed ways to obtain a bouncing cosmology but did not address the relaxation ofthe CC.In this paper, we present a dynamical relaxation model for the CC problem. Here, a rolling scalar field takes theuniverse from a natural, large positive cosmological constant (CC) to a small negative one. At this point the universewill begin to contract. The contraction increases the energy in the field(s) responsible for the tuning of the CC. Atsome density, this increased energy, through a small coupling, reheats other matter. The energy density in this matterblue-shifts as the universe continues to contract. We take as an assumption that this sector can trigger dynamicsthat causes the scale factor to bounce at short distances, allowing the universe to expand and produce our observedcosmological history. This bounce could occur through vorticity as in [14], but any other possible bounce model(e.g. through NEC violating fluids) would work as well. The relaxation mechanism is independent of the bounce andcosmology that comes after. The simplest model (section II) naturally tunes a CC scale as large at 10 MeV to anegative CC of scale 1 meV. We show how a few additional fields (section III) and stages of rolling allow one to scana CC from scales much higher than 1 TeV, and end with a positive
CC of order the current critical density. Thesemodels are experimentally testable (section IV), through astrophysical, cosmological, and laboratory probes.This existence proof factorizes the solution of the CC problem into an infra-red (IR) part that accomplishes thesensitive tuning of the CC and an ultra-violet (UV) sector whose purpose is to accomplish a cosmological bounce athigh densities. Importantly, the UV dynamics is decoupled from the IR tuning. This existence proof highlights theimportance of short-distance descriptions of a cosmological bounce, and presents the opportunity to re-imagine thesource of the initial perturbations often credited to inflation. a r X i v : . [ h e p - ph ] F e b II. SIMPLE MODEL
We now show a simple model that naturally tunes the CC from up to ∼ (10MeV) to ∼ − (1meV) and subsequentlyreheats the universe when the universe contracts. In the following section, we will show how one can increase theinitial cosmological constant and end on a small positive one.In this model, a rolling scalar field starts from a point with large vacuum energy (abiding eternal inflation bounds).As it rolls down, the CC decreases, eventually going through zero. At this point, the universe begins to contract ata parametrically smaller negative CC. The universe contracts to a large energy density, and assuming a cosmologicalbounce, re-expands until today, while the field value does not evolve significantly, thus keeping the vacuum energysmall. Finally, reheating is shown to be trivially accomplished by extracting energy from the rolling condensate anddumping it into a thermal bath through derivative couplings.Remarkably, all of this can be accomplished with the dynamics of the following example model: L = 12 ( ∂φ ) + 14 F (cid:48) F (cid:48) + ¯ ψ ( i /D − m ψ ) ψ − m A (cid:48) A (cid:48) + g φ − φf F (cid:48) ˜ F (cid:48) , (1)Here, φ is a scalar field with a softly broken shift symmetry, A (cid:48) µ is a massive photon whose gauge field strength is F (cid:48) and ψ is a charged massive Dirac fermion. We use a mostly negative metric and we have defined the value φ = 0to be the point of vanishing cosmological constant for convenience. The rolling of φ decreases the CC and its kineticenergy is eventually converted to the gauge bosons A (cid:48) . The energy in this radiation reheats the universe, producingboth the standard model and the degrees of freedom necessary to cause the universe to bounce. A. Rolling to -meV The dynamics of interest (see Figure 1) start with an initial condition of an expanding universe, and a large negativevalue for φ , namely φ = φ < | φ | (cid:29) M p (the reduced Planck scale). Here, there is positive vacuum energy,Λ ≡ g φ , and assuming initially ˙ φ (cid:28) Λ , φ slow rolls down its potential in a vacuum dominated universe.We would like the universe to evolve to a small cosmological constant and avoid eternal inflation. This puts theconstraint g > (cid:112) / πH (cid:39) (cid:112) / π Λ / (3 M p ) (2)or a limit on the highest CC that can be relaxed of Λ (cid:46) (cid:112) gM p . (3)The rolling continues until φ reaches a value φ = φ ∼ − M p , where the kinetic energy surpasses the potential energyand φ is less than a Hubble time away from the origin. Now, with the kinetic energy increasing, the potential energybecomes negative and the Hubble scale decreases at an increasing rate, as one can clearly see from the FriedmannEquations: H = 8 π G N (cid:18)
12 ˙ φ − g φ (cid:19) (4)˙ H = − πG N ˙ φ (5)The Hubble rate H vanishes in a finite time when a value φ ∼ M p is reached. To see this analytically, take Eq. (5)and integrate from point φ to where Hubble vanishes:0 − H = − (cid:90) πG N ˙ φ dt = − πG N (cid:90) ˙ φdφ (6) H > πG N ˙ φ ∆ φ (7)where the 2 subscript indicates the values at point φ where kinetic and potential energy are equal. The inequalitycomes from the fact that H is monotonically decreasing from (5) making ˙ φ monotonically increasing due to its equationof motion. This allows us to replace ˙ φ in the integral with its minimum (initial) value to generate the inequality.Thus, φ traverses a finite distance (of order M p ) in a finite time (as can also be shown numerically). Because ˙ H < H continues to decrease below zero and the universe begins to contract. The potential energy at thispoint is ∼ − Λ ≡ − g M p . I l ~ - dog land FIG. 1: Evolution of φ . The field starts rolling from point 1 with vacuum energy Λ . At point 2, φ ’s kinetic energy equals itspotential energy, Λ , a distance ∼ M p away from the origin. At point 3, the Hubble scale passes through zero and the vacuumenergy is ∼ − Λ . The point 4 represents the position of φ after a period of kinetic-energy dominated contraction, where thevacuum energy as decreased to ∼ − Λ log (1 /a rh ) and a rh is the scale factor at reheating. At point 5, the reheated universehas expanded until today and φ has moved a negligible amount from point 4. B. Kinetic Energy during Contraction
As the universe contracts, the kinetic energy of φ quickly dominates the potential energy and blue-shifts as ˙ φ =˙ φ a − , where a is the scale-factor of the Friedmann-Robertson-Walker metric and ˙ φ is the velocity of φ at the pointwhere H = 0 (taking a = 1 at that point). Taking the kinetic energy to be dominant, it is simple to compute thedistance φ travels down its potential while the universe contracts:∆ φ = (cid:90) ˙ φ dt ≈ − (cid:90) ˙ φ a √ (cid:112) πG N / φ a da ≈ √ M p log (1 /a ) (8)where we used the first Friedmann equation (4), the fact that H <
0, and the definition of the reduced Planck mass.We see that the distance traveled by φ (and thus the change in the potential energy) is only logarithmically sensitiveto the scale factor. For example, if the universe contracts to a scale where the energy density is ˙ φ = Λ , where a = (Λ / Λ hot ) / , then ∆ φ = √ M p (2 /
3) log (Λ hot / Λ ), or about 60 M p for Λ hot = 1 TeV and Λ = 1 meV. Thisdefines the point 4 in Figure 1. And so in the time that φ rolls parametrically not far beyond point 3, the energydensity in the universe (kinetic energy in φ ) increases all the way up to an essentially arbitrarily high scale Λ .The final limit on the highest CC that can be naturally scanned down to ∼ meV then arises as follows. In orderto avoid tuning we want the negative CC reached after contraction to be ∼ meV in magnitude. The value of theCC at point 5 in Figure 1 is very close to the value at point 4 (as we will see below). So we need the value at point 4 g φ (cid:46) meV . Combining this with equation (8) and equation (3) givesΛ (cid:46) (cid:18) √ hot / Λ ) (cid:19) meV M p (9)Or Λ (cid:46)
10 MeV for Λ hot = 1 TeV and Λ = 1 meV, with only a very weak dependence on Λ hot . So this model cannaturally reduce a CC of 10 MeV to the observed value, reducing fine-tuning by roughly 40 orders of magnitude. C. Reheating A (cid:48) Now we utilize φ ’s coupling to the massive vector to convert the kinetic energy of φ into a thermal bath of A (cid:48) µ and ψ . We describe this process in two stages:(1st Stage) A population of vectors will be produced when ˙ φ/f > m A (cid:48) due to an instability in the mode equationfor the vectors [15, 16, 22]: ¨ A (cid:48)± + ( k ± k ˙ φf + m A (cid:48) ) A (cid:48)± = 0 (10)where A (cid:48)± are the spatial Fourier transforms of circularly polarized modes of the vector A (cid:48) . The A (cid:48)− modes with wavenumbers k < ( ˙ φ/f ) and k > ( f m A (cid:48) / ˙ φ ) will be exponentially growing modes. Assuming an initial fluctuation of order A (cid:48)− ∼ k in each mode (the minimum set by quantum mechanics), the largest energy density comes from modes k ∼ ˙ φf in the massive vector and so grows up to ρ A (cid:48) ∼ ( ˙ φ/f ) e
2( ˙ φ/f ) t th , (11)where t th is the thermalization time scale. This initial density sources the thermal destruction of the condensate aswe see below.In order to be sure that the thermal calculations we rely on below (in (14)) are valid, we would like to have the A (cid:48) thermalize when their energy density is above m A (cid:48) so that the temperature they thermalize at is above m A (cid:48) . Evenif this is violated, the thermal friction may well still stop the rolling of φ as we want, but it is outside the regime ofvalidity of the thermal field theory calculations that have been done, so we will choose to avoid this region. As wewill see this is easy to do.We compute a process that will start thermalizing the A (cid:48) . Assuming the fermion ψ ’s mass is of order m A (cid:48) , the crosssection for scattering A (cid:48) A (cid:48) → ψψ is σv ∼ α m A (cid:48) (12)where α = e π and e is the charge of ψ , and we have assumed that ˙ φf ∼ m A (cid:48) so the A (cid:48) produced by the rolling φ aresemi-relativistic. We want the vector to thermalize only after the energy density in A (cid:48) has reached m A (cid:48) . The numberdensity of A (cid:48) at this point is n ∼ m A (cid:48) . Then the scattering rate at this point isΓ A (cid:48) A (cid:48) → ψψ ∼ nσv ∼ α m A (cid:48) (13)We need this to take longer than the time it takes φ to produce this energy density which is, using (11), t th ∼ f / ˙ φ ) log ( m A (cid:48) f / ˙ φ ). Again, assuming ˙ φf ∼ m A (cid:48) , this simply requires that α < / √
2. Satisfying this bound allowsthe energy density in the vectors to grow to at least m A (cid:48) before thermalization begins.(2nd Stage) After there are some A (cid:48) particles around, they will produce ψ particles, and they can then thermalizerapidly. This then leads rapidly to a thermal bath which will then cause thermal friction of φ (see e.g. [17]). We willtake the thermal friction coefficient to be: Γ ∼ π α T f (14)where T is the temperature of the thermal bath, and the coupling α is the renormalized coupling at T . This is theresult for the analogous friction in the case of pure Yang-Mills, and the numerical factor is a result of normalizationdifference with [17]. There may be numerical differences between the computed Yang-Mills case and the yet to becomputed Abelian case with charged fermions, but it is unlikely to be parametrically different (at least with respectto f and T dependence), and thus we will use this value as a rough estimate of the friction. This coefficient will damp and ultimately suppress the φ rolling such that its energy density is negligible relativeto that of the radiation. We want to see how long this process takes. Initially while this damping is happening, thebath is at some temperature T which is less than the kinetic energy in φ . Then the rate at which energy density isbeing taken out of the φ rolling and put into the thermal bath is dρdt ∼ π α T f ˙ φ (15)Setting the energy density in the thermal bath to ρ ∼ T we find dTdt ∼ π α ˙ φ f (16)This means that it takes the longest time to get it up to the highest temperature (the time is dominated by the UV).We can see how long it takes to remove an O(1) fraction of the kinetic energy in φ by setting ˙ φ to be a constant equalto its initial value when we want to reheat: ˙ φ ∼ T . Then we see how long it takes T to get up to this value, callthis time ∆ t reheat . We can see that ∆ t reheat ∼ f π α T (17)Requiring that this happen within a Hubble time then means that we set ∆ t reheat (cid:46) H ∼ M p T . This then leads tothe requirement that f (cid:46) π α T reheat M p . (18)Plugging in ˙ φ ∼ T ∼ f m A (cid:48) , we also have the constraint f (cid:46) π / ( m A (cid:48) M p ) / . From equation (18) we can seethat f has a wide range of possible values even for a fairly low reheating temperature. And f can even be all the wayup near M p if we take a high reheating temperature. D. Bouncing
The contracting universe needs to bounce (evolve to an expanding universe) so that it can re-expand and reproduceour cosmic history. The dynamics responsible for the bounce can be decoupled from the tuning of the cosmologicalconstant. This is easily accomplished - after all the kinetic energy of φ is dumped by thermal friction into A (cid:48) , the A (cid:48) can reheat the degrees of freedom responsible for the bounce through weak couplings. Around the time of thebounce, these degrees of freedom must effectively violate the null energy condition or be able to trigger vorticity inextra dimensions as in [14]. In order to bounce, these degrees of freedom need to blue-shift faster than the othermatter content in the universe so that they are relevant at the short distances where the bounce occurs. Further,when the universe re-expands after the bounce, this matter must return to its original state so that the tuning ofthe CC is not affected. This can likely be guaranteed if the behavior of this sector is determined by thermodynamics(such as a temperature), wherein the re-expansion of the universe would cool this sector, returning it to its originalstate. However, we leave an explicit model of the way to trigger the bounce starting with our CC relaxation modelfor future work.Before the bounce φ rolls only a short distance as seen above. And all the initial kinetic energy of φ from contractionis dumped into the ( A (cid:48) µ , ψ ) sector by thermal friction, so φ is then moving very slowly. Once the universe bounces andre-expands, it is dominated by radiation, specifically in the ( A (cid:48) µ , ψ ) sector, and the kinetic energy in φ is never aboveits terminal velocity value ˙ φ ∼ g / Γ , which (as can be easily shown) keeps φ from rolling a significant amount. At One can also build a model where the sector being heated is non-Abelian and the mass scale m A (cid:48) corresponds to the confinement scaleof the strong group. In order to avoid the generation of larger barriers in the φ potential, on can add an additional massless quark. Andalso, this model is close enough to that studied in [17] to suggest the damping rate Γ is parametrically the same. any time, even if thermal friction becomes smaller than Hubble friction, φ will not roll more than ∆ φ ∼ (Λ /T ) M p in a Hubble time. So φ does not roll significantly during the entire contaction, bounce, and subsequent expansion ofthe universe. Thus the potential energy of φ is not changed significantly and so the dynamical relaxation solution forthe CC is not spoiled. E. Reheating the Standard Model
The last step would be to reheat the rest of the universe (namely the Standard Model sector). This can beaccomplished by coupling the massive vector to the normal matter through mixing with the photon or through higherdimensional operators, and allowing the Standard Model to thermalize at some point when the temperature is higherthat the scale of big bang nucleosynthesis.A kinetic mixing with the hypercharge gauge boson, (cid:15)F (cid:48) µν F µνY would allow the vector A (cid:48) decay into standard modelparticles with a rate Γ decay ∼ α Y (cid:15) m A (cid:48) . Equating this with the Hubble scale, this gives the temperature of theUniverse at the time of decay: T d ∼ α / Y (cid:15) (cid:112) m A (cid:48) M p . On the other hand, the mixing of the vectors produces aneffective coupling to φ of the form (cid:15) ( φ/f ) F Yµν ˜ F µνY . If we require – though it may not be necessary – the rate ofthis instability (from the analogous version of equation (10) for photons), to be less than Hubble, ˙ φ(cid:15) /f < H , thennone of the dynamics described above change. This constraint is most sensitive at the lowest values of Hubble, where˙ φ ∼ g /H ∼ H M p , thus requiring f > (cid:15) M p . One can show that these constraints are trivial to satisfy. III. EPICYCLES
The model presented above naturally takes a large cosmological constant and relaxes it to a parametrically smaller(albeit negative) one, converts the energy from this sector to a hot standard model, and (after a bounce) producesnormal big bang cosmology with a tiny cosmological constant. In this section, we will show how a few additionaldegrees of freedom will allow: (a) dynamics that produce a small positive cosmological constant and (b) naturalrelaxation from much larger cosmological constants, while maintaining technical naturalness.
A. Positive CC
Suppose the CC has been reduced to ∼ − meV , with the energy in the rolling field dumped into other forms ofmatter. At this stage, the universe starts crunching and the energy density in these matter fields will blue-shift. Thisenergy can be used to trigger a technically natural phase transition at the ∼ meV scale, resulting in an addition tothe vacuum energy and the CC changing from ∼ − meV to ∼ +meV . This transition is not fine-tuned so long asthe CC is changed to a positive value of roughly the same size as, or greater than, the small negative value it hadafter relaxation. Once the universe has already started to crunch, changing the CC by ∼ +meV does not change thedynamics of the universe as its energy density is dominated by the rapidly blue-shifting matter or radiation density.Thus, the rest of the cosmic evolution necessary to implement our framework such as the bounce and the subsequentre-expansion of the universe are unaffected by the transition necessary to achieve a positive CC. There are likely manyways to accomplish this goal, we present one such example in Appendix A. B. Larger Cutoff
The principal difficulty in achieving a larger cutoff in the model presented in section II is that the slope of φ needs tobe sufficiently large when the vacuum energy is big in order to avoid eternal inflation. This large slope induces a kineticenergy for φ that causes it to roll well beyond ∼ meV . To achieve a larger cutoff, we can introduce an additionalrolling scalar field (Φ) that has a steeper slope. We start the universe with both φ and Φ rolling. When the CC islarge, the rolling of Φ provides the clock necessary to avoid eternal inflation. As the CC approaches ∼ (10 MeV) ,we need to create barriers that stop the rolling of Φ. Once Φ is stopped at ∼ (10 MeV) , the rolling of φ will furtherrelax the CC down to ∼ − meV .How can we naturally trigger barriers for Φ? The key idea is to observe that when the CC is large, the large Hubblefriction results in a low terminal velocity for Φ. As the CC drops, Hubble friction decreases, resulting in a largerterminal velocity. We use this to trigger the barriers. The increased velocity of Φ can trigger instabilities in gaugefields to which Φ is derivatively coupled, such as in the models discussed in section II. The energy released in thisprocess can be used to raise barriers for Φ. There are many ways to accomplish this goal - we present a proof ofconcept model in Appendix B. Note, this initial stage of relaxation does not require a bounce since the relaxationends at a relatively high, positive value of the CC. Models of this kind could potentially also be used to simply relaxthe value of the CC in inflationary relaxion models where the CC could be reduced from the cut-off to the weak scale. IV. POTENTIAL SIGNATURES
There are four generic elements of our construction: a rolling scalar field φ that cancels the bare cosmologicalconstant, a bounce in our immediate past to reheat the universe, a phase transition that should occur at scales ∼ meV in order to push the cosmological constant to slightly positive values after it becomes negative and strongdynamics at various scales (for example, ∼
10 MeV) that enable the cutoff of the theory to be above the TeV scale.Each of these elements can be separately tested.Any dynamical relaxation model requires a field that scans the CC. So a relatively model-independent signature ofthis framework is that the kinetic energy of the rolling field φ gives rise to a non trivial equation of state for dark energy.Current bounds on the equation of state of dark energy imply that the velocity ˙ φ of the field is (cid:47) . [18]. Inthe simplest of our models, we expect ˙ φ ∼ g M pl / meV with g (cid:47) meV /M pl . Thus, current and near future probesof the equation of state of dark energy are constraining the simplest models that can solve the cosmological constantproblem. In addition to cosmological probes, the kinetic energy ˙ φ can also be probed in laboratory experiments if φ has couplings to the standard model. Of course, this coupling is necessary at some level since the kinetic energy of φ has to reheat the universe just before the bounce, ultimately resulting in our existence. Radiative stability of φ andefficient reheating implies that φ must couple derivatively to the standard model, much like an axion. There are twoleading interactions that can be experimentally probed: the coupling of φ to electromagnetism and nucleon/electronspins via the operators φf a F ˜ F and ∂ µ φf a ¯Ψ γ µ γ Ψ respectively. The electromagnetic coupling is already constrained– if ˙ φ ∼ . , current constraints on B modes in the CMB require f a (cid:39) M pl as this coupling will cause thepolarization of CMB photons to rotate as they propagate through the evolving dark energy [19]. Interactions withnucleon/electron spins can potentially be probed through tests of Lorentz symmetry since the evolving dark energyprovides a cosmic background that is being searched for in these experiments [20]. These signatures are relativelymodel-independent signatures of a dynamical relaxation solution to the CC problem.In our model, a cosmic bounce is required. This could be detected through a cosmological background of stochasticgravitational waves. The hubble scale during a bounce is not constant - thus, the gravitational wave spectrum wouldexhibit a sharp feature corresponding to the minimum of the bounce, unlike inflationary cosmology that produces anearly scale invariant spectrum. The detection of stochastic gravitational waves at different frequency bands wouldenable experimental discrimination between these two possible cosmological scenarios in our immediate past.One of the simplest ways to obtain a slightly positive cosmological constant after the rolling of φ makes it slightlynegative is to reheat a hidden sector that undergoes a phase transition at the ∼ meV scale. This suggests that theuniverse could contain a hidden sector of dark radiation around the meV scale, with a phase transition likely to occurin this sector. Such a transition would also indicate an evolving equation of state of dark energy. Moreover, it wouldalso be interesting to directly search for dark radiation in laboratory experiments, building on the work that hasoccurred in recent years on searching for ultra-light dark matter.Finally, we expect the existence of confining sectors at scales such as ∼ MeV in order to push the cutoff of thetheory to scales above ∼ TeV. It would be interesting to develop techniques to search for such confining sectors - forexample, this sector might contain degrees of freedom such as glueballs which can interact with the standard model.A generic operator analysis suggests that these interactions are suppressed. However, since these particles are light,they could conceivably be probed in high statistics intensity frontier experiments [23].
V. DISCUSSION AND CONCLUSIONS
We have shown a technically natural way to solve the CC problem. Our framework takes a large positive CC andreduces it to a small negative CC through dynamical relaxation. This causes the universe to crunch. At the sametime, the relaxation process also naturally dumps energy into a new sector. The energy densities in this sector blue-shifts during contraction thus reheating the universe to high temperatures. By using the energy in the new sector,we are also able to naturally push the CC to positive values after the universe begins to crunch. This new sector hasto ultimately be responsible for instigating a cosmological bounce so that the universe can re-expand, giving rise tothe present day universe that had a hot big bang but with a small CC, though we do not model the bounce here.This framework overcomes all the obstacles faced by dynamical relaxation methods to solve the CC problem suchas the empty universe problem in Abbott’s model, Weinberg’s “no-go” theorem [12] and the problem of standardmodel phase transitions. Many aspects of this framework lead to testable consequences, some of which require thedevelopment of new experimental probes to target these specific signatures.An important fact about this construction is that the UV dynamics of the bounce are decoupled from the IRrelaxation process. Moreover, the CC itself does not change significantly during the bouncing phase. There is thusconsiderable freedom to attach the IR relaxation phase to any UV dynamics permitting a cosmological bounce [4–11, 13, 14]. To be considered a complete solution to the CC problem, we need to identify the specific mechanism thatwould allow the hot matter in the crunching universe to trigger a bounce. This requires a better understanding ofthe matter sources necessary to create a bounce.To reproduce observational facts about our universe, it is important to identify mechanisms that would give riseto the scale invariant spectrum of perturbations that have been observed in the CMB. Since our model is largelyjust inflation, but with reheating accomplished through a bounce, there are elements of scale invariance built intothe mechanism. For example, the rolling scalar field will have nearly scale invariant fluctuations until the CC goesthrough zero. It would be interesting to see if these fluctuations could seed the observed spectrum of perturbations.At the very least, a period of inflation could follow the bounce.Another direction worthy of exploration is to see if the relaxion paradigm that solves the hierarchy problem canbe successfully incorporated into this CC relaxation mechanism. Cosmological relaxation appears to be the onlydynamical mechanism that has the potential to solve the naturalness problems associated with both the CC and theweak scale. These phenomena find a natural home in a universe that is much older than conventionally assumed —something that is observationally possible and theoretically interesting.
Acknowledgments
We thank Peter Arnold, Savas Dimopoulos, Michael Fedderke, Nemanja Kaloper, Mikko Laine, Guy Moore andRaman Sundrum for discussions. SR was supported in part by the NSF under grants PHY-1638509 and PHY-1507160, and the Simons Foundation Award 378243. PWG acknowledges the support of NSF grant PHY-1720397,DOE Early Career Award de-sc0012012. DEK acknowledges the support of NSF grant PHY-1214000. This work wassupported in part by Heising-Simons Foundation grants 2015-037, 2015-038, and 2018-0765, DOE HEP QuantISEDaward
Appendix A: Axion Model for Positive CC
As discussed in the main article, a simple mechanism for generating a positive CC is by allowing the thermal bathto generate a phase transition to a vacuum with a higher vacuum energy (by an amount meV ). Here we present anexplicit model, though many others are possible.Take the following low-energy potential for an axion-like field, χ : V ( χ ) = Λ cos nχf − ˜Λ cos χf (A1)where n is any small integer bigger than 2 and λ > meV and ˜Λ ∼ meV. This potential is shown in Figure 2. Wetake this potential to be periodic with period ∼ f . Take χ to be in some random minimum after the CC relaxationdiscussed above. During the crunching universe, it is possible for the Λ sector to thermalize while the ˜Λ sector doesnot. This can happen if the ˜Λ confinement scale is much higher, but the sector has a small quark mass allowing˜Λ (cid:28) Λ. During the contracting phase, once the temperature of the Λ sector rises above its confinement scale, thebarriers would disappear. Then χ will begin to roll toward the minimum of the ˜Λ potential. But then this velocitywill rapidly blue-shift because of the contraction. It is easy to check that there is a large parameter space where thislarge velocity will cause χ to go around the entire ∼ f period of the periodic potential many times. As the universecools, it will ultimately end up in other random mimimum, in general different from the one it started in. And thereis thus an O (1) chance it will be a higher minimum with a net positive CC. This method is not tuned so long as theoriginal relaxation mechanism tuned the CC down to a negative value of ∼ meV or lower. Then the amount addedto the CC by this axion field is what determines the CC today. χ - - V FIG. 2: A sketch of the potential in equation (A1). The blue line is the potential at low temperatures, the orange line is thepotential at high temperatures where only the second term in equation (A1) is on. This axion begins in some random minimumin its potential and then ends up in a different one. The vacuum energies of the minima are split by meV so this can easilyraise the negative CC left after relaxation to a positive CC. Appendix B: Higher Cutoff
In this section, we show a simple, proof of principle, model with one extra rolling field (beyond the one that wasalready used in the model of Section II) that allows us to push the cutoff well above ∼
10 MeV. More generally, thissetup describes how the CC can be relaxed from ∼ Λ a to a lower, positive value ∼ Λ d . This stage can precede themodel presented in section II, thus instrumenting the full tuning. In fact it could also replace the model of SectionII as the lowest stage of relaxation, except that since it leaves a positive CC the universe does not naturally startcrunching and heat up, so there would have to be some other way to avoid the ‘empty universe’ problem.The model of this stage is also one of a rolling field tuning the CC. During the rolling, however, the field is coupledto a non-Abelian group (which can generate barriers) plus an Abelian group (to generate additional friction): L = 12 ( ∂ Φ) + 14 G (cid:48) G (cid:48) + 14 F (cid:48) F (cid:48) + κ Φ − Φ f th G (cid:48) ˜ G (cid:48) − Φ f A F (cid:48) ˜ F (cid:48) , (B1)where G (cid:48) and F (cid:48) are the field strengths of the non-Abelian and Abelian groups (with indices suppressed), and theremaining parameters are couplings. During the rolling, the coupling of the rolling field, Φ, to non-Abelian gaugebosons produces a non-trivial background temperature for the non-Abelian group (akin to what happens in warminflation [21]). The temperature remains high enough for a long period during which the instanton-generated potentialbarriers do not form. Once the Hubble scale becomes low enough, the Abelian group begins to extract energy fromthe rolling during which the background temperature decreases. When the temperature become low enough, thebarriers form and stop the rolling of Φ, fixing the CC to a value parametrically smaller than its initial value. Duringthis entire period – as we will show at the end of this sub-section – the model of the previous section does not evolvemuch until Hubble becomes small enough that it can begin to roll consequently.Thus, while the model in (B1) is quite simple, the dynamics associated with rolling are quite non-trivial, and wedescribe them in the following subsections chronologically for a rolling field. A summary of the motion is presentedin Figure 3.
1. Early Friction: Hubble
The field Φ is taken to be slowly rolling (as an initial condition) at point (a) with vacuum energy Λ a ) and Hubblescale H ( a ) ∼ Λ a ) /M p . In addition, there is assumed to be a small background temperature in the non-Abelian group’s0 I I : n a b C d a b C d FIG. 3: Evolution of Φ. In the upper figure, the field rolls from (a) where Hubble friction dominates and the backgroundtemperature grows quickly, to (b) where thermal friction dominates and the temperature grow very slowly, to (c) where Φ’srolling is dominantly suppressed by coherent production of Abelian fields and the temperature drops. At point (d), thetemperature drops to the confinement scale of the non-Abelian group and barriers form. degrees of freedom, a temperature much smaller than Λ ( a ) , but larger than the strong group’s confinement scale, whichwe will define as Λ th . The latter will be the scale of the non-perturbative, low-energy potential for Φ.During to the rolling, the background temperature is maintained in the non-Abelian sector due to thermal frictionΓ( T ) as described in the last section via Eq. (14). One can compute an approximate steady-state temperature duringthis rolling period by requiring the red-shifting/cooling of the plasma bath is counteracted by the heating from thecondensate: 0 (cid:39) dT dt ≈ − HT + (cid:18) π α (cid:48) ( T ) T f th (cid:19)
12 ˙Φ (B2)where in front of T there is really a factor that includes the number of degrees of freedom in the non-Abelian sectorwhich we take to be ∼
1. When Hubble friction dominates, we can take ˙Φ (cid:39) ( κ /H ), its terminal velocity duringslow roll. Thus, for temperatures well above the confinement scale, and taking a nominal value of the coupling α (cid:48) (and ignoring its temperature dependence), one finds a quasi-steady-state temperature of T (cid:39) κ / (9 H ¯ f ) , (B3)where we defined ¯ f ≡ f th / (16 πα (cid:48) / ).1
2. Late Friction: Thermal
As Φ rolls, the Hubble scale will slowly drop and the temperature will rise, eventually reaching a point when thedominant friction is due to the thermal bath, i.e. , Γ( T ) (cid:39) T / ¯ f ∼ H . This again produces a terminal velocity forΦ, namely ˙Φ (cid:39) κ / Γ (cid:39) κ ¯ f / ( T ) , and a temperature, using (B2), of T (cid:39) (cid:18) κ ¯ f H (cid:19) . (B4)As Hubble slowly decreases, the temperature stays nearly constant. Without the Abelian sector, the temperaturewould eventually dominate Hubble and the universe would become radiation dominated (as in warm inflation [15]),and Φ would roll to negative values of the cosmological constant. We instead would like the temperature to dropbefore this happens.
3. Mode Instability: Cooling
The coupling of Φ to the Abelian sector produces an instability in some modes of the gauge fields. An effectivenegative mass term appears in the equation of motion for the A (cid:48) + modes with k < ˙Φ /f A as in equation (10), but with m A (cid:48) set to zero. Following the analysis of [15], one can show that a quasi-steady-state is reached when a rolling fieldis coupled in this way to an Abelian group and vacuum energy dominates the energy density of the universe. In thisregime, ˙Φ (cid:39) ξf A H (B5)where ξ = 12 π log (cid:32) · π α A M p f A κ V (Φ) (cid:33) ∼ −
100 (B6)where α A is the fine-structure constant for the Abelian group. The velocity in (B5) becomes the terminal velocitywhen this instability becomes the dominant energy-loss mechanism, which is equivalent to when (B5) is smaller thanthe terminal velocity due to thermal friction, κ / Γ which happens when:
H < κ ¯ f ¯ f A T (B7)where for simplicity, we define ¯ f A ≡ ξf A . With this velocity, thermal friction still extracts energy from the rollinggenerating a quasi-steady state temperature: T (cid:39) ¯ f A H f (B8)which thus decreases as Φ rolls. The rolling thus eventually stops as the temperature drops to the confinement scaleof the non-Abelian group and barriers in the Φ potential (from instanton effects) begin to form. This must occurwhile vacuum energy is still dominating, both over the temperature bath and over the energy density in the Abelianfields, which is estimated to be [15] κ ¯ f A .
4. Constraints on the Initial Cosmological Constant
Now we can use the above behavior to put constraints on parameters, including the initial value of the CC. Theconstraints are as follows:Λ th (cid:46) T ( a ) Temp is higher than confinement scale at start.Λ th ∼ T ( d ) Temp lowers to confinement scale at end.Λ th > κ f th Barrier’s slope beats underlying slope.Λ d ) > κ ¯ f A CC dominates over Abelian mode growth.2Additional constraints, such as the requirement that the final vacuum energy is greater than the energy density inthe thermal bath and that thermal fluctuations do not lead to eternal inflation at any time during the scan, can beshown to be weaker than those above for the ranges of parameters here.Combining these constraints with the steady-state temperatures at the beginning and end of Φ’s roll (Equations(B3) and (B8), respectively), one can derive the following bound:Λ ( a ) < Λ / d ) M / p ¯ f / (B9)where we have taken α (cid:48) = O (1). A strict constraint on Λ ( a ) could come from requiring f th = 16 πα (cid:48) / ¯ f to be greaterthan Λ ( a ) (the CC at top), producing: Λ ( a ) < Λ / d ) M / p (B10)which, for Λ ( d ) = 10 MeV, produces a bound Λ ( a ) <
300 GeV. If instead, one requires f th to be only greater than thehighest temperature achieved, which occurs when thermal friction gives way to mode instability – when Eq.s (B4)and (B8) are equal – then one finds: Λ ( a ) < Λ / d ) M / p (B11)or for Λ ( d ) = 10 MeV, Λ ( a ) <
300 TeV.One can easily show that the lower stage field φ , described in the previous section, does not evolve appreciably in itspotential. This in fact is not a phenomenological constraint, as the upper stage described in this section will alwaysstop at a positive CC, allowing the lower stage to work (as long as φ is not at the bottom of its potential). Nevertheless,one can check that during the upper stage, both classical rolling and evolution due to quantum fluctuations of φ arenegligible.In principle there could be several axion fields with varying scales in their potentials. This could give rise to severalstages of the relaxation model described here which could raise the initial CC far beyond the 300 GeV - 300 TeVscales found here. We leave the exploration of multiple stages for future work. [1] L. F. Abbott, Phys. Lett. , 427 (1985). doi:10.1016/0370-2693(85)90459-9[2] T. Banks, Phys. Rev. Lett. , 1461 (1984). doi:10.1103/PhysRevLett.52.1461[3] P. W. Graham, D. E. Kaplan and S. Rajendran, “Cosmological Relaxation of the Electroweak Scale,” Phys. Rev. Lett. , no. 22, 221801 (2015) doi:10.1103/PhysRevLett.115.221801 [arXiv:1504.07551 [hep-ph]].[4] P. J. Steinhardt and N. Turok, Science , 1436 (2002) doi:10.1126/science.1070462 [hep-th/0111030].[5] P. Creminelli and L. Senatore, JCAP , 010 (2007) doi:10.1088/1475-7516/2007/11/010 [hep-th/0702165].[6] D. Battefeld and P. Peter, Phys. Rept. , 1 (2015) doi:10.1016/j.physrep.2014.12.004 [arXiv:1406.2790 [astro-ph.CO]].[7] P. Creminelli, M. A. Luty, A. Nicolis and L. Senatore, JHEP , 080 (2006) doi:10.1088/1126-6708/2006/12/080 [hep-th/0606090].[8] V. A. Rubakov, Theor. Math. Phys. , 1651 (2006) [Teor. Mat. Fiz. , 409 (2006)] doi:10.1007/s11232-006-0149-7[hep-th/0604153].[9] P. Creminelli, A. Nicolis and E. Trincherini, JCAP , 021 (2010) doi:10.1088/1475-7516/2010/11/021 [arXiv:1007.0027[hep-th]].[10] V. A. Rubakov, Phys. Usp. , 128 (2014) [Usp. Fiz. Nauk , no. 2, 137 (2014)] doi:10.3367/UFNe.0184.201402b.0137[arXiv:1401.4024 [hep-th]].[11] K. Balasubramanian and S. P. Dabholkar, arXiv:1401.7015 [hep-th].[12] S. Weinberg, Rev. Mod. Phys. , 1 (1989). doi:10.1103/RevModPhys.61.1[13] L. Alberte, P. Creminelli, A. Khmelnitsky, D. Pirtskhalava and E. Trincherini, JHEP , 022 (2016)doi:10.1007/JHEP12(2016)022 [arXiv:1608.05715 [hep-th]].[14] P. W. Graham, D. E. Kaplan and S. Rajendran, Phys. Rev. D , no. 4, 044003 (2018) doi:10.1103/PhysRevD.97.044003[arXiv:1709.01999 [hep-th]].[15] M. M. Anber and L. Sorbo, Phys. Rev. D , 043534 (2010) doi:10.1103/PhysRevD.81.043534 [arXiv:0908.4089 [hep-th]].[16] A. Hook and G. Marques-Tavares, JHEP , 101 (2016) doi:10.1007/JHEP12(2016)101 [arXiv:1607.01786 [hep-ph]].[17] M. Laine and A. Vuorinen, Lect. Notes Phys. , pp.1 (2016) doi:10.1007/978-3-319-31933-9 [arXiv:1701.01554 [hep-ph]].[18] N. Aghanim et al. [Planck Collaboration], arXiv:1807.06209 [astro-ph.CO].[19] Y. Akrami et al. [Planck Collaboration], arXiv:1807.06211 [astro-ph.CO].[20] M. Pospelov and M. Romalis, Phys. Today , 40 (2004). doi:10.1063/1.1784301[21] R. Rangarajan, arXiv:1801.02648 [astro-ph.CO]. [22] B. A. Campbell, N. Kaloper, R. Madden and K. A. Olive, Nucl. Phys. B399