Prepared for submission to JHEP
de Sitter Decays to Infinity
Patrick Draper,a Isabel Garcia Garcia,b and Benjamin Lillarda aDepartment of Physics, University of Illinois, Urbana, IL 61801 bKavli Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106, USA E-mail: [email protected], [email protected], [email protected]
Abstract: Bubbles of nothing are a class of vacuum decay processes present in some theories with compactified extra dimensions. We investigate the existence and properties of bubbles of nothing in models where the scalar pseudomoduli controlling the size of the extra dimensions are stabilized at positive vacuum energy, which is a necessary feature of any realistic model. We map the construction of bubbles of nothing to a four-dimensional Coleman-De Luccia problem and establish necessary conditions on the asymptotic behavior of the scalar potential for the existence of suitable solutions. We perform detailed analyses in the context of five-dimensional theories with metastable dS S1 vacua, using analytic 4 × approximations and numerical methods to calculate the decay rate. We find that bubbles of nothing sometimes exist in potentials with no ordinary Coleman-De Luccia decay process, and that in the examples we study, when both processes exist, the bubble of nothing decay rate is faster. Our methods can be generalized to other stabilizing potentials and internal manifolds. arXiv:2105.10507v1 [hep-th] 21 May 2021 Contents
1 Introduction1
2 Vacuum Decay in the Presence of Gravity3 2.1 Coleman-De Luccia formalism3 2.2 Bubble of nothing as a Coleman-De Luccia problem6 2.3 Sources of moduli potentials 10 2.4 Classical vs effective potentials 11
3 Bounce of Nothing 12 3.1 Bubble existence conditions 12 3.2 Structure of bubble solutions with dS S1 vacua 14 4 × 3.3 A toy model 16 3.4 A bounce of nothing 17 3.5 de Sitter decay rate 20 3.6 Application to more generic potentials 22
4 Other Exotic Bounces 22
5 Numerical Construction of Bounce Solutions 26 5.1 Preliminaries 26 5.2 Exponentially falling potentials 27 5.2.1 Parameter space of BON solutions 27 5.2.2 Instanton solutions and action 30 5.3 Exponentially growing potential: five-dimensional CC 32
6 Conclusions 34
A Solutions to Equations of Motion 37
B Analytic Continuation 43
C Specific Potential 44
1 Introduction
Our Universe is entering a phase of dark energy domination. To date, cosmological obser- vations are consistent with an epoch of late-time acceleration driven by a small and positive cosmological constant [1–3]. Although an epoch of accelerated expansion can be easily accommodated within gen- eral relativity, building full-fledged de Sitter vacua in the context of string theory — while seemingly not impossible [4,5] — has proven surprisingly hard [6–20]. This has led to con- jecture and debate about whether metastable de Sitter vacua can exist in consistent theories
– 1 – of quantum gravity [21–23], and if so, whether there is an upper bound on their lifetimes [24]. For example, it is well-known that in theories with compactified extra dimensions, if the effective potential for the geometric moduli admits minima of positive energy density, these minima are at best local. In turn, a four-dimensional de Sitter vacuum may ‘sponta- neously decompactify,’ decaying into a Universe with zero vacuum energy and spatial extra dimensions of infinite size [25–27]. Depending on the details of the higher-dimensional construction, additional decay channels can also be present and are generally described by either the Coleman-De Luccia (CDL) [28] or Hawking-Moss (HM) [29] instantons. In this paper, we focus on an alternative instability: the generalization of Witten’s bubble of nothing (BON) [30] to spacetimes with four-dimensional de Sitter factors and compact internal manifolds with stabilized moduli. Witten’s bubble describes an insta- bility of the Kaluza-Klein (KK) vacuum M4 S1. The gravitational instanton is the × five-dimensional Euclidean Schwarzschild solution, and the endpoint of the decay is the destruction of the spacetime. Via dimensional reduction, Witten’s bubble can be recast as a solution to a CDL problem of a real scalar field minimally coupled to gravity in four dimensions. The scalar field plays the role of the radial modulus setting the size of the S1, and Witten’s solution is recovered in the case of a vanishing scalar potential. This four- dimensional description, however, makes it natural to consider generalizations of the BON in the presence of a potential for the moduli, as first suggested in [31]. In this way, non- perturbative instabilities associated with moduli stabilization, including generalizations of the bubble of nothing, can be studied from the bottom-up. A number of studies have explored relatives of Witten’s bubble from a top-down per- spective, including in string compactifications [32–36], in higher-dimensional models [37], and in flux compactifications with stabilized moduli [38–43]. In particular, Refs. [40, 41] consider BON instabilities in 6D models with explicit mechanisms for moduli stabilization at a positive value of the vacuum energy. Our approach is complementary to previous work, attempting to factorize as much as possible the specific details of moduli stabilization from the existence and construction of bubbles of nothing. Our goals are to describe, from the four-dimensional perspective, existence criteria for BON solutions in theories with stabilized moduli; to develop analytic approximations for the BON action in the presence of a scalar potential; to compare the corresponding decay rate to those of the ordinary CDL and HM channels; and to study numerically those examples that are analytically intractable. This paper is a longer, more detailed companion to [44], where some of the main ideas and techniques were presented and illustrated with an eye toward minimal, streamlined presentation. Here we review and extend the results of [44], deriving some of the assertions, expanding the range of models studied, elaborating on the broader context of the work, and undertaking a more extensive numerical survey. We perform our most detailed analysis in the context of a five-dimensional model with a range of moduli potentials. Depending on the features of the potential, the BON may be the only accessible instability, or it may exist alongside more conventional CDL or HM bounces. In the latter case, the BON decay rate is always the fastest in the examples we study. Unlike the familiar CDL and HM solutions, the CDL solution corresponding to Wit- ten’s bubble exhibits singular boundary conditions at the center of the bounce, with di- vergent values for the scalar and its derivatives. In the ordinary CDL problem, the scalar field satisfies equations of motion analogous to those of a particle moving on an inverted
– 2 – potential. The solution can then be found via an overshoot-undershoot method, with the ‘particle’ starting at rest at the center of the bounce, and asymptotically rolling over to the (inverted) false vacuum. This picture can be extended to bounces that satisfy BON bound- ary conditions. In this case, the particle does not begin at rest, but rather it starts moving ‘from infinite distance with infinite velocity,’ corresponding to the singular boundary con- ditions. Implementing an overshoot-undershoot algorithm appropriate to accommodate these boundary conditions allows us to find BON instabilities of four-dimensional de Sitter vacua numerically, verifying and going beyond our analytic results. This paper is organized as follows. In Section2, we review the CDL formalism, the properties of Witten’s bubble most relevant to our analysis, and some of the familiar sources of moduli potentials that can be present in higher-dimensional models. In Section3 we study BON decays of four-dimensional de Sitter vacua. We find necessary conditions on the asymptotic form of the scalar potential for the CDL equations to admit a BON bounce solution, and, under some assumptions, we obtain an approximate analytic solution for the bounce. We evaluate the action of these instantons, showing that it remains finite in the limit of vanishing false vacuum energy density, and that for a large class of models it is mostly independent of the behavior of the scalar potential in the compactification regime. We discuss a qualitatively different type of instanton that also satisfy BON boundary conditions in Section4, and argue that these may be thought of as superpositions of a BON together with an ordinary CDL or HM bounce. In Section5 we present a detailed numerical study of BON bounce solutions and their Euclidean actions for a wide range of parameter space. We conclude in Section6 and outline directions for future work. A series of Appendices collect various additional results of a more detailed or tangential nature and are referred to throughout the text.
2 Vacuum Decay in the Presence of Gravity
In Section 2.1, we summarize the CDL formalism describing vacuum decay in the presence of gravity. We review the properties of Witten’s BON instability that will be most relevant to our work in Section 2.2, with special focus on how it can be recast as the solution to a four-dimensional CDL problem, as first done in [31]. In Section 2.3 we discuss several sources of moduli potentials that are naturally present in theories with extra dimensions. We conclude with some qualitative observations regarding the interpretation of moduli potentials in Section 2.4. As emphasized in the Introduction, our focus will be on instabilities of theories with compactified dimensions, stabilizing potentials for the geometric moduli, and local de Sitter vacua. Fig.1 shows the qualitative picture that we will have in mind throughout our work.
2.1 Coleman-De Luccia formalism The semiclassical treatment of vacuum decay in quantum field theory was developed in [45, 46]. To leading order, the decay rate of the false vacuum is given by the action of the O(4)-symmetric solution to the Euclidean field equations known as the “bounce.” In Euclidean signature the bounce interpolates between true and false vacuum, and via analytic continuation, it also provides the initial condition and real time evolution of true vacuum bubbles after nucleation. The decay rate per unit volume of the false vacuum is
– 3 –
compactification decompactification
+ 1 1 Figure 1: Theories featuring compactified extra dimensions, and moduli spaces that are stabilized at a positive value of the vacuum energy, are typically susceptible to decay into a higher-dimensional, decompactified theory, [25–27], as indicated by the φ 0 region in ≥ this figure. On the other hand, the behavior of the scalar potential in the compactification limit may vary widely depending on the details of the higher-dimensional construction. Our work will be concerned with instabilities of four-dimensional de Sitter vacua into the compactification regime, where φ 0. The location of the de Sitter false vacuum is ≤ arbitrary, and we choose φ 0 for convenience. fv ≡ given by 4 ∆S Γ v e− , (2.1) ∼ where v is a typical energy scale, and ∆S is the action of the bounce minus the action of the Euclidean false vacuum. Coleman and De Luccia further developed the formalism of [45, 46] to include the effects of gravitation [28]. In four dimensions, the O(4) symmetry of the bounce restricts the form of the relevant metrics to those compatible with the form:
2 2 ds4 = dξ + ρ(ξ) dΩ3, (2.2) V ( ) V ( ) U( ) where ξ is a radial coordinate and ρ sets the curvature radius of the transverse 3-sphere. For a real scalar field minimally coupled to gravity, the corresponding field equations read
3ρ0 ∂U(φ) φ00 + φ0 = , R( ) 0 (2.3) R( ) R5( ) 0 R5( ) ρ ∂φ ! !1 ! !1 2 2 ρ 2 ρ0 = 1 + 2 φ0 2U(φ) . (2.4) 6Mp − + + 1 1 1 1 The relevant boundary conditions accompanying Eqs. (2.3) and (2.4) depend on the topol- ogy of the solution. For a de Sitter false vacuum, and a scalar potential that is 0, the ≥ solution has the topology of an S4. In this case, ρ vanishes at two values of ξ, corresponding to a minimum and a maximum of the radial coordinate [47]. Without loss of generality, we can choose ξ 0. Then, ξ [0, ξ ], and ξ is defined via ρ(ξ ) ρ(0) = 0. min ≡ ∈ max maxcompactificationmax ≡ Any non-singular solution to the equations of motion must then also obey the boundary conditions φ0(0) = φ0(ξmax) = 0. decompactification
R !1 R 0 – 4 – ! For example, the Euclidean solution describing the de Sitter false vacuum is given by
ξ φ (ξ) 0, and ρ (ξ) = Λ sin , (2.5) dS ≡ dS Λ with Λ 3M 2/U . The traditional CDL bounce that describes false vacuum decay ≡ p fv is characterizedq by an internal region — the inside of the bubble — where φ lies in the vicinity of the true vacuum, and a transition to an exterior region where φ is near its value in the metastable phase. If false and true vacua are nearly degenerate, the bounce is well described within the “thin-wall” approximation, where the bubble wall thickness (the region over which the field transitions from false to true vacuum) is much smaller than its overall radius. For a Minkowski true vacuum, the curvature radius of a thin-wall bounce is given by ρ0 3σ ρ¯ 2 , with ρ0 = , (2.6) ' 1 + (ρ0/2Λ) Ufv
φfv 1 and σ is the bubble wall tension given by σ dφ 2(U(φ) U(φfv)) [28]. ' φtv − The tunneling rate is controlled by the Euclidean´ action of the bounce. The action of p a scalar field minimally coupled to Einstein gravity in four dimensions is given by
m2 1 S = d4x√g Pl + gµν∂ φ∂ φ + U(φ) (2.7) E ˆ − 2 R 2 µ ν ξmax = 2π2 dξ ρ3U (on-shell). (2.8) − ˆ0 The second line only applies “on-shell,” that is, when evaluating the action on a solution to the Euclidean field equations consistent with the ansatz of Eq. (2.2). In the thin-wall limit, the tunneling exponent of the CDL bounce is given by [28]
B0 ∆SCDL SE SE 2 2 , (2.9) ≡ CDL − dS ' [1 + (ρ0/2Λ) ]
with B 27π2σ4/(2U 3 ). For a small false vacuum energy density, U σ2/M 2, 0 ≡ fv fv p 2 4 24π Mp ∆SCDL . (2.10) ' Ufv The CDL bounce is guaranteed to exist provided the scalar potential satisfies
∂2U(φ) 4 > , (2.11) ∂φ2 Λ2 φ=φtop top
where φ is the location of the potential barrier, and Λ 3M 2/U . If this condition top top ≡ p top is not satisfied, the CDL bounce may fail to be present [48q–51]. However, even in the absence of a CDL solution, there is a second type of instability, corresponding to the
1 The thin-wall approximation is applicable provided Ufv − Utv Utop − Ufv, with Utop the value of the p scalar potential at the top of the potential barrier. To leading order, σ ∝ Utop, and corrections involving the difference Ufv − Utv are dropped within the thin-wall approximation.
– 5 – homogeneous solution where φ sits at the top of the potential barrier. This instanton was first discussed by Hawking and Moss (HM) [29], and it is given by
ξ φ φ , and ρ (ξ) = Λ sin . (2.12) HM ≡ top HM top Λ top In this case, the corresponding tunneling exponent reads
2 4 1 1 ∆SHM SE SE = 24π Mp . (2.13) ≡ HM − dS Ufv − Utop
A common feature of the CDL and HM instabilities is that their tunneling exponents diverge in the limit U 0, as can be seen from Eqs. (2.10) and (2.13). As the false and fv → true vacua become degenerate, the CDL and HM decay channels become ineffective, and the metastable vacuum becomes exponentially long-lived. We finish this section noting that Eq. (2.3) can be rewritten in integral form as follows:
ξ=ξ ξmax 3 max 3 ∂U ρ φ0 = dξ ρ . (2.14) ξ=0 ˆ0 ∂φ
In most cases, the left-hand side of Eq. (2.14) is known. For example, it vanishes for the CDL bounce, on account of the boundary conditions φ0(0) = φ0(ξmax) = 0, as well as on the de Sitter and HM solutions where φ is constant. This feature permits the construction of exact solutions to the equations of motion in certain simple cases, as we will illustrate in Section 3.3, where Eq. (2.14) supplies the initial conditions for the bounce solution.
2.2 Bubble of nothing as a Coleman-De Luccia problem The bubble of nothing discovered by Witten is an unusual instability of the KK vacuum M4 S1 [30]. The decay can be described semiclassically, in terms of a gravitational × instanton equivalent to the five-dimensional Euclidean Schwarzschild solution:
1 R2 − R2 ds = dr2 1 5 + r2dΩ + 1 5 dy2. (2.15) 5 − r2 3 − r2 In these coordinates, r [R , ) and y y + 2πR . As usual, the relationship between ∈ 5 ∞ ∼ 5 the y-coordinate periodicity and the Schwarzschild radius ensures the absence of a conical singularity at r = R5. The near-horizon metric is of the form
ds dλ2 + λ2dy˜2 + R2dΩ , (2.16) 5 ' 5 3 where λ (2R )(r R ), andy ˜ y/R . Indeed, the geometry remains smooth near ≡ 5 − 5 ≡ 5 r = R , and the topology of the solution is R2 S3. As depicted in the left panel of Fig.2, 5 p × the bubble geometry ends on a smooth cap at r = R5, which we will typically refer to variously as the bubble wall or core. For this BON solution, the scale R5 characterizes three distinct physical scales: the size of the bubble (the area radius of the transverse S3); the magnitude of local five-dimensional curvatures near its wall; and the asymptotic radius of the KK circle. As described in [30], the initial conditions for the nucleated bubble, as well as its Lorentzian-time evolution, can be obtained by analytic continuation of Eq. (2.15). Static time slices asymptote to the false vacuum, but near the origin the geometry appears as
– 6 – U =0BON R ⇡ 2 L/
L(⇠)=2⇡R exp 2 /M 3 p ⇣ q ⌘
⇠/R
Figure 10: Just L( )
1.0 U =0BON vacuum 0.8
0.6 5 R ⇡ 2
instanton L/ 0.4
0.2 2 L(⇠)=2⇡R5 exp 3 /Mp 0.0 ⇣q ⌘ 0.0 0.5 1.0 1.5 2.0 2.5 3.0 antipodal points of S3 ⇠/R
Figure 3: Diagram R5 Figure 11: Just L( ) Figure 2: Left: Diagram showing the geometry of the KK false vacuum (“vacuum”) and BON (“instanton”) solutions. The r direction lies along the symmetry axis, with the periodic y coordinate along the azimuth. In the false vacuum, φ 0, and the KK radius ≡ is independent of r. For Witten’s bubble, the KK radius approaches5 zero near r & R5, whereas far away from the bubble wall the geometry asymptotically approaches that of the false vacuum. Right: After dimensional reduction, the scalar modulus φ describes the proper length of the KK circle, L. As a function of the CDL coordinate, defined in Eq. (2.20), L vanishes at ξ = 0, corresponding to the location of bubble wall, whereas L 2πR in the asymptotic region ξ R . ' 5 5 though a hole has been carved out of the original manifold: from the four-dimensional perspective, spacetime ends on the growing three-sphere x2 t2 = R2. Nevertheless, − 5 the presence of a fifth dimension ensures that the solution remains smooth: although asymptotically the size of the KK radius is 2πR5, it smoothly shrinks to zero at the bubble wall located at x2 t2 = R2. Unlike in the ordinary picture of tunneling in quantum field − 5 theory, the bubble that gets nucleated is not a bubble of true vacuum, but rather it has ‘nothing’ inside. At first sight,Figure Witten’s 4: Left: Right: bubble appears to be a rather different instability from CDL tunneling, taking place in a purely gravitational higher dimensional theory. However, it was observed in [31] that the BON can be rewritten as a solution to a four-dimensional CDL problem via dimensional reduction. In the dimensionally reduced theory, the radial modulus responsible for2 setting the size of the KK circle is mapped onto the CDL scalar φ. The Einstein frame, canonically-normalized reduction is obtained via the following parametrization of the five-dimensional metric:
q 2 φ q 2 φ 3 M 2 3 M 2 ds5 = e− p ds4 + e p dy . (2.17)
– 7 – -U 1.5 (⇠)/Mp U( )
bon(⇠) 1.0 (ϕ⇠ R ) bon ⌧ 5 (⇠ R ) bon 5 BON 0.5
CDL
0.0
-0.8 -0.6 -0.4 -0.2 0.0 0.2
⇠/R5
Figure 5: Cartoon
1.5 (⇠)/M p ⇢(⇠)/R5 0.0 1.0 -0.2 1.5
-0.4 0.5 -0.6 1.0
-0.8 ⇢bon(⇠) bon(⇠) 0.0 0.5 -1.0 (⇠ R ) ⇢bon(⇠ R5) bon ⌧ 5 ⌧ (⇠ R ) ⇢bon(⇠ R5) -0.2 0.0 0.2 0.4 0.6 0.8 -1.2 bon 5 0.0 -1.4
0 1 2 3 4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2
⇠/R5 ⇠/R5
Figure 3: The five-dimensional EuclideanFigure 6: Schwarzschild Cartoon gravitational instanton, or Wit- ten’s BON, expressed in terms of the degrees of freedom φbon(ξ)(left) and ρbon(ξ)(right) of the corresponding CDL problem. The behavior of the solutions in the near-horizon (ξ R ) and asymptotic (ξ R ) regimes is3 indicated in the figure, following Eq. (2.21) 5 5 and Eq. (2.22), respectively.
The Euclidean Schwarzschild solution (2.15) can be written in this form, with
M 3 R2 φ (r) = p log 1 5 , (2.18) bon 2 2 − r2 r 1/2 1/2 R2 − R2 ds = dr2 1 5 + 1 5 r2dΩ . (2.19) bon − r2 − r2 3 The 4D metric in Eq. (2.19) can be further rewritten as in Eq. (2.2), with the CDL radial coordinate related to the Schwarzschild radial coordinate as2 r drˆ ξ(r) . (2.20) ≡ ˆ (1 R2/rˆ2)1/4 R5 − 5 The corresponding CDL degrees of freedom satisfy Eq. (2.3) and (2.4) with U(φ) 0, and ≡ the profile of φbon(ξ) and ρbon(ξ) is shown in Fig.3. Near the bubble wall, at ξ = 0,
2 3ξ 3ξ 1/3 φ (ξ R ) M log , and ρ (ξ R ) R , (2.21) bon 5 ' p 3 2R bon 5 ' 5 2R r 5 5 whereas asymptotically far from the bubble,
M 3 R 2 φ (ξ R ) p 5 , and ρ (ξ R ) ξ + γR , (2.22) bon 5 ' − 2 2 ξ bon 5 ' 5 r 2An exact expression for ξ can be obtained in terms of a hypergeometric function, given in AppendixA.
– 8 – with numerical constant
1 γ √πΓ 3 Γ 1 − 0.6. (2.23) ≡ 4 4 ' Exact solutions for φbon and ρbon can be found in AppendixA. Eq. (2.21) reflects an important difference between Witten’s BON and an ordinary CDL bounce: the former appears singular at ξ = 0, while the latter satisfies φ0(0) = 0. Instead, the BON must be thought of as a CDL problem with singular boundary conditions, describing tunneling from φ = 0 to φ = in the absence of a potential barrier. −∞ Dimensionally reducing the five-dimensional Einstein-Hilbert action one finds
2 4 Mp 1 µν Mp ∆S = d x√g + g ∂µφ∂νφ φ (2.24) ˆ (− 2 R 2 − √6 )
Mp rr 3 + √g d x√h∂rφ + ∆SGHY, (2.25) √6 ˆ∂ M where ∂ refers to the boundary manifold defined as the hypersurface r r , and ∆S M ≡ 0 GHY is the (appropriately regularized) Gibbons-Hawking-York boundary term evaluated on the boundary. On the BON solution, the combination M 2 +(∂φ)2 vanishes, and ∆S − p R GHY → 0 in the limit r .3 Integrating by parts the φ term and combining it with the first 0 → ∞ term of Eq. (2.25), the BON action can be written as
2 2 3 ∆Sbon = π Mp ρbon(ξ)φbon0 (ξ) (2.26) 3 ξ=0 r = π2M 2R2. (2.27) p 5 Although we have chosen to evaluate the right-hand side of Eq. (2.26) at ξ = 0, it follows 3 from inspecting Eq. (2.3) that the combination ρ φ0 is a constant on the BON. In particular,
3 3 2 ρbon(ξ)φbon0 (ξ) = MpR5 ξ. (2.28) r2 ∀ Compared to the five-dimensional perspective, the dimensionally-reduced formulation of the BON appears rather pathological. The four-dimensional geometry of Eq. (2.19) exhibits a naked singularity at r = R5 with diverging curvatures and scalar gradients. However, the combination M 2 + (∂φ)2 that enters into the four-dimensional action − p R identically vanishes, since five-dimensional Euclidean Schwarzschild is a solution of Ein- stein’s equations in vacuum. All 5D curvatures remain sub-Planckian, provided R5 is larger than the Planck length. The role of φ as a geometric modulus restricts the terms in the reduced action, organizing them into finite five-dimensional quantities. On the other hand, as was noted in [31], the advantage of the four-dimensional formula- tion is that it makes it straightforward to look for instanton solutions with BON boundary conditions in the presence of a scalar potential that stabilizes the geometric modulus. In the absence of supersymmetry, a runaway scalar potential will be generated by Casimir energies, so modulus stabilization is a necessity in any realistic construction. It is then
3 2 µ The boundary term appropriate for the BON metric can be written as SGHY = −Mp n ∂µVol(∂M) = 2 2 2 2 −6π Mp (r0 −R5/2), whereas the boundary term of the flat metric with the same boundary geometry reads (flat) 2 2 2p 2 2 (flat) SGHY = −6π Mp r0 1 − R5/r0. In total, ∆SGHY ≡ SGHY − SGHY , which vanishes in the limit r0 → ∞.
– 9 – interesting to ask whether BON solutions still exist and whether the corresponding decay rate is significantly modified. In the presence of a non-zero U(φ), the region near the bub- ble wall will be sensitive to the asymptotic form of the potential in the compactification direction, while the region near the metastable vacuum will depend on the details of the mechanism that stabilizes the modulus. In Section3, this useful factorization will allow us to classify BON instabilities of de Sitter vacua in terms of the behavior of the scalar po- tential in the compactification limit, while treating the region near the metastable vacuum with toy models.
2.3 Sources of moduli potentials As described in the introduction, we wish to employ a bottom-up strategy, remaining as agnostic as possible about the details of the stabilization mechanism leading to a four- dimensional de Sitter vacuum, and the detailed features of the effective potential relevant for the moduli. Nevertheless, it will be useful to keep in mind some of the sources of moduli potentials that appear naturally in theories with extra dimensions, especially as it pertains to their behavior in the compactification limit. This will allow us to gain some intuition about the regime of validity of some of the conclusions we draw in later sections. Sources of potentials for geometric moduli in theories with extra dimensions have been widely discussed in the literature (see e.g. [27]). For simplicity, let us first consider 4 + n-dimensional gravity, with the n extra di- mensions compactified. Allowing for a non-zero cosmological constant (CC), the higher- dimensional Einstein-Hilbert action reads M 2+n S = d4+nx√g + M 2+nΛCC , (2.29) E ˆ 4+n − 2 R4+n 4+n where M is the fundamental Planck mass. The Einstein frame dimensional reduction is obtained through the following parametrization: q q 2n φ 2 2 φ 2 n+2 Mp 2 n(n+2) Mp 2 ds4+n = e− ds4 + e dsn, (2.30) and the reduced action reads 2 Mp 1 Mp 2n S = d4x√g + gµν∂ φ∂ φ φ + U(φ) , (2.31) E ˆ µ ν (− 2 R 2 − 2 rn + 2 ) with 2 q φ q φ Mp 2n+4 2n U(φ) = e− n Mp + M 2ΛCC e− n+2 Mp . (2.32) − 2 Rn p 4+n For example, if the internal manifold is an n-sphere with curvature radius R , then = n Rn n(n 1)/R2 . In the compactification limit, the contribution to the scalar potential from − n this internal curvature always dominates over the contribution from the higher-dimensional CC, and therefore U as φ . → −∞ → −∞ More generally, the scalar potential of geometric moduli might receive contributions from other non-trivial field configurations present in the higher-dimensional theory, such as branes and fluxes. For example, an n form flux field wrapped around the internal − manifold leads to a contribution that scales as [27] q 3 2n φ U (φ) e− n+2 Mp . (2.33) flux ∝
– 10 – In the compactification limit, a contribution of this kind would dominate over that coming from internal curvature. Some simple (and non-exhaustive) possibilities for the behavior of the scalar potential in the compactification limit are therefore of the form
• U(φ) 0 (no flux, higher-dimensional CC, or internal curvature); ∼ q 2n φ • U(φ) M 2ΛCC e− n+2 Mp (no flux or internal curvature; CC dominated); ∼ p 4+n q M 2 2n+4 φ • U(φ) p e− n Mp (no flux; internal curvature dominated). ∼ − 2 Rn In Section 3.1 we will discuss what sources of moduli potentials are compatible with the existence of a BON instability of four-dimensional de Sitter vacua.
2.4 Classical vs effective potentials We have primarily discussed classical sources of moduli potentials, such as higher dimen- sional flux, the cosmological constant, and geometric curvature. These sources can produce strange, asymptotically divergent potentials upon dimensional reduction. If the potential is singular in the compactification limit, we should view the singularity similarly to the way we view the singularities of the lower-dimensional Ricci scalar and modulus kinetic term in those limits: they are artificial, induced by the choice of field-space coordinates, and they are absent in the higher-dimensional theory. What of quantum corrections to the effective potential? These may play some impor- tant role in the stabilization of the modulus, as in [4]. However, it is difficult to include these corrections systematically in BON solutions. Near the bubble wall, in the higher dimensional description, curvatures and momenta are of order the asymptotic KK scale. There is no systematic truncation in the derivative expansion, and the full effective action is needed to get a consistent quantum correction. In lieu of performing such a complete computation, we model the stabilization part of the potential, and assume that the asymp- totic potential, relevant near the bubble wall, is dominated by classical effects. In fact, this is typically a reasonable approximation even if quantum effects are important for the stabilization. Near the bubble wall at ξ = 0, field gradients are of order the KK scale, and they dominate the equations of motion (along with the classical potential, if it grows rapidly enough). Our classical BON solutions could become invalid if the EFT breaks down. Two ways this can happen are if (higher-dimensional) curvatures or stress energy become large compared to the Planck or string scales, or if a large tower of heavy states becomes light. The first possibility generally does not occur in BON geometries if the asymptotic volume modulus is larger than the Planck or string scales. The second possibility, the appearance of a tower of light states, is at first sight more concerning. Bubbles of nothing probe an infinite distance in moduli space, and the swampland distance conjecture states that under such an excursion an infinite tower of states should descend below the cutoff [52]. In ordinary KK theory, these states correspond to wound string modes, which become light as the radius of the KK circle goes to zero. However, this pathology does not actually arise in BON geometries, because the excursion is highly localized. Ref. [52] studied wound strings in the vicinity of a static BON and found that, while they become less massive than asymptotic wound states, they do not become lighter than the string scale. This has
– 11 – a simple explanation: near the smooth cap marking the bubble wall, wound strings simply slip off, becoming ordinary unwound states. There is no localized tower of light states, and string scale effects are expected to be small.
3 Bounce of Nothing
In this section, we focus on the description of BON instabilities of four-dimensional de Sitter vacua in the presence of a potential for the moduli. Our starting point is the Euclidean action of Eq. (2.25), but now supplemented by a scalar potential U(φ):
2 4 Mp 1 µν Mp SE = d x√g + g ∂µφ∂νφ φ + U(φ) . (3.1) ˆ (− 2 R 2 − √6 ) Since our choice of false vacuum is the four-dimensional de Sitter solution, no boundary terms need to be included in Eq. (3.1), as the manifolds relevant to our discussion are always compact. Provided that the four-dimensional metric is compatible with the ansatz of Eq. (2.2), the corresponding equations of motion are Eq. (2.3) and (2.4). Therefore our goal is to find solutions to the CDL field equations that asymptote to the de Sitter false vacuum, but that feature BON boundary conditions at the center of the bounce. In this section, we will focus on obtaining analytic approximations for both the bounce solution and its Euclidean action that determines the decay rate of the false vacuum. Rather than being comprehensive, we aim to establish the circumstances under which BON solutions can exist in the presence of a moduli potential, and to study the properties of the decay rate in the limit of small vacuum energy. In Section 3.1, we establish the conditions that moduli potentials need to satisfy in the compactification limit to be compatible with a bubble of nothing solution. We discuss the structure of solutions describing the decay of a dS S1 false vacuum into a BON 4 × in Section 3.2. In Section 3.3, we introduce a toy model for the potential of the radial modulus that allows us to obtain an analytic description of the bounce in Section 3.4. Finally, in Section 3.5 we compute the leading contributions to the tunneling exponent, paying special attention to its behavior in the limit U 0. In Section 3.6 we briefly fv → describe the application of our analytic methods to more general potentials.
3.1 Bubble existence conditions In order for BON-like instantons to exist, certain conditions must be satisfied by the scalar potential U(φ) in the compactification limit. We can establish these conditions by exam- ining the behavior of the CDL equations, Eqs. (2.3) and (2.4), near ξ = 0. We concentrate on moduli potentials with the asymptotic behavior U(φ) U eaφ/Mp in the compactifi- ∼ 0 cation direction, corresponding to φ . We emphasize, however, that we make no → −∞ prior assumption regarding the sign of a: it will turn out that some potentials that grow exponentially in the compactification limit also admit BON solutions. Near ξ = 0, we consider a generalization of Witten’s bubble given by the following ansatz: ξ ξ β φ(ξ) M α log , and ρ(ξ) ηR , (3.2) ' p ` ' 5 ` with parameters α, β, `, and η to be determined self-consistently. The space of solutions to the CDL equations depends on the dimensionality of the internal manifold. We focus on
– 12 – -U 1.5 (⇠)/Mp U( )
bon(⇠) 1.0 ( ⇠ R ) bon ⌧ 5 (⇠ R ) bon 5 BON 0.5
CDL
0.0
-0.8 -0.6 -0.4 -0.2 0.0 0.2
⇠/R5
Figure 5: Cartoon -U 1.5 1.5
U( ) 1.0 1.0
BON 0.5 0.5
CDL
0.0 0.0
-0.8 -0.6 -0.4 -0.2 -0.00.2 0.00.2 0.2 0.4 0.6 0.8
Figure 4: A diagram comparingFigure the BON 6:Cartoon and CDL boundary conditions. The BON initial conditions φ(0) , φ + ensure that φ(ξ) quickly approaches the false → −∞ 0 → ∞ vacuum. Its approach to φ(ξ) φ is determined by the details of the potential, which → fv is approximately quadratic in the neighborhood of φ φ . This direct approach to φ ≈ fv fv is indicated by the straight blue line in the figure. The CDL solution, on the other hand, begins at finite φ(0) = φCDL with φ0(0) =3 0. It subsequently “rolls” in the inverted 1.5 potential, as indicated by the dashed red line, until it approaches φfv. In the dashed blue line we show a second BON example, one in which the potential U(φ ) diverges → −∞ (with U(φ) shown in dark gray). The φ (0) BON initial conditions can compensate − 0 → ∞ for the divergence in U(φ), so that there is a BON solution even if there is no corresponding 1.0 CDL solution for that potential.
cases where the internal manifold contains an n-sphere submanifold that shrinks to zero size at the bubble wall. It is convenient to discuss the cases n = 1 and n 2 separately. ≥ 0.5 n = 1 Provided a > √6, it follows from inspecting Eqs. (2.3) and (2.4) that the ansatz Eq. (3.2) − provides a solution to the field equations near ξ = 0 so long as α = 2/3 and β = 1/3. ± 0.0 Further demanding that the corresponding five-dimensional geometry is free of conical p singularities imposes the restriction ` = 2R5/3 and selects the positive sign for α. In total, the small-ξ behavior of the solutions reads -0.2 0.0 0.2 0.4 0.6 0.8
2 3ξ 3ξ 1/3 φ(ξ) M log , and ρ(ξ) ηR . (3.3) ' p 3 2R ' 5 2R r 5 5 Eq. (3.3) is identical to the near-horizon geometry of the original Witten bubble, given in Eq. (2.21), except we now have an additional parameter η > 0. The physical significance of η will become apparent in Sections 3.2 and 3.4. Notice that a negative range of a is compatible with BON boundary conditions, despite the fact that U in the compact- | | → ∞ ification limit when a < 0. For example, this would be the case for a potential dominated
– 13 – by a five-dimensional CC, with a = 2/3, as discussed in Section 2.3. On the other CC − hand, no solutions compatible with our ansatz are possible when a √6. For example, p ≤ − axion flux wrapping the S1 yields a = √6, and the flux must be screened in order flux − for the S1 to shrink to zero size at the bubble wall. An example of this type was studied in [39]. n ≥ 2 The discussion above can be easily generalized to cases where a higher dimensional sphere shrinks to zero size at the center of the bounce. (The appropriate parametrization of the higher-dimensional metric was given in Eq. (2.30).) There are no solutions to the field equations when a < 2(n + 2)/n, and if a > 2(n + 2)/n the absence of a conical − − singularity is only possible in the case n = 1, corresponding to our previous discussion. p p However, when n 2, the marginal value a = 2(n + 2)/n admits a different class of ≥ − BON solutions. In this case, the ansatz of Eq. (3.2) corresponds to a smooth solution to p the Euclidean field equations provided 2n n 2R α = , β = , and ` = n , (3.4) rn + 2 n + 2 n + 2 where Rn is the radius of curvature of the n-sphere, and only if the overall scale of the potential is given by 2 Mp U0 = n(n 1) 2 . (3.5) − − 2Rn Inspecting the near-horizon metric reveals that the topology of these solutions is of the form R1+n S3. Overall, for this class of solutions to exist, the behavior of the scalar × potential in the compactification limit must therefore be of the form
2 n(n+2) φ Mp − n Mp U(φ) n(n 1) 2 e . (3.6) ' − − 2Rn Interestingly, this is precisely the behavior induced by the internal curvature of the n- sphere, as reviewed in Section 2.3. Examples in 6D Einstein-Maxwell theory with a shrink- ing S2 were studied in [41, 53]. Notice that the contribution to the scalar potential from n-form flux wrapping the internal manifold scales with an exponent a = 3 2n/(n + 2), as discussed in Sec- flux − tion 2.3. When n 2, a < 2(n + 2)/n, and therefore flux must either be absent, ≥ flux − p screened, or otherwise disappear near the bubble wall in order for BON solutions to exist, p in keeping with Gauss’ law. It would be interesting to generalize our analysis further to admit, for example, pointlike screening charges along the lines of [53].
In the remainder of this manuscript, we restrict our attention to potentials that satisfy a > √6 in the compactification direction, in which case the internal manifold that resolves − the apparent singularity at ξ = 0 must have dimension n = 1.
3.2 Structure of bubble solutions with dS S1 vacua 4 × As discussed in the Introduction, our focus on this paper is on the existence of BON instabilities of four-dimensional de Sitter vacua. Let us assume the existence of a false vacuum of the form dS S1, with Euclidean metric 4 × 2 2 2 ds5 = dξ + ρdS(ξ) dΩ3 + dy , (3.7)
– 14 – with y y + 2πR , and ρ (ξ) as given in Eq. (2.5). The parameter R is the radius of the ∼ 5 dS 5 KK circle in the false vacuum, and it is an input parameter defining the UV theory. Our goal is to construct instanton solutions that asymptotically approach this false vacuum, but that satisfy boundary conditions at ξ = 0 as given in Eq. (3.3). Unlike the original BON corresponding to an instability of M4 S1, the Euclidean geometry describing this × class of instantons will be compact, with ξ [0, ξ ] and ρ(ξ ) = 0. ∈ max max Eq. (3.3) corresponds to the small-ξ limit of the one-parameter family of BON solutions introduced in [44], given by
3 3/2 3/2 3/2 φη(ξ) Mp log η + φbon(ξ η− ), and ρη(ξ) η ρbon(ξ η− ). (3.8) ≡ r2 · ≡ · When ξ R , φ and ρ are indeed as in Eq. (3.3), whereas for ξ R , we find 5 η η 5 3 η3 R 2 φ (ξ) M log η 5 , and ρ (ξ) ξ + γη3/2R . (3.9) η ' 2 p − 2 ξ η ' 5 r ( ) Eq. (3.8) provide exact solutions to the Euclidean equations of motion when U(φ) 0. In ≡ the presence of a scalar potential, φη and ρη remain approximate solutions for a limited range of ξ, as we will soon discuss. As described in [44], it is only in this latter case that the parameter η acquires physical significance. Near ξ = 0, the five-dimensional metric corresponding to this solution can be written as
ds dλ2 + λ2y˜2 + η2R2dΩ , (3.10) 5 ' 5 3 with λ R (3ξ/2R )2/3 andy ˜ y/R . This expression clarifies the meaning of η: it ≡ 5 5 ≡ 5 parametrizes the curvature radius of the S3 of the near-horizon geometry, R ηR , and 3 ≡ 5 corresponds to the radius of the ‘hole’ that nucleates in spacetime. Having established that Eq. (3.3) provides acceptable boundary conditions at ξ = 0 in the presence of a scalar potential, we now turn to estimate the range of ξ over which the solution behaves approximately like the BON. Comparing the right-hand side of Eq. (2.3) to, say, the first term on the left side, we find (ignoring irrelevant (1) factors) O q 2+ 2 a ∂U/∂φ U R2 ξ 3 a 0 5 for ξ R , (3.11) φ ∼ M 2 R . 5 00 p 5 where we have used U(φ) U eaφ/Mp for the asymptotic form of the scalar potential in ∼ 0 the compactification limit. Provided that the overall scale of the potential satisfies U 0 M 2/R2, the above ratio remains small in the regime ξ R , since a < √6 for potentials p 5 . 5 − that admit BON solutions, as discussed in Section 3.1. The transition away from BON behavior therefore happens in the asymptotic regime ξ & R5. At that point, φ/Mp . 1, and the modulus is already approaching its value in the false vacuum. Expanding the scalar potential in this region as U(φ) U + 1 m2φ2, we have ' fv 2
∂U/∂φ 2 2 m ξ for ξ & R5. (3.12) φ00 ∼ The presence of a scalar potential becomes significant once ρ(ξ) ξ m 1. Self- ∼ ∼ − consistency of our analysis requires that the mass of the modulus in the false vacuum
– 15 – - U
⇢/R ⇢/R ϕ ⇤ ⇤ • •
=12 =12 ⇢/R • • ⇢ ⇢ R R instanton instanton 1 1 vacuum vacuum & & ⇤ ⇢ ⇢ iue5 iga for Diagram 5: Figure R R • 2 2
⇢/R =12 ⌧ ⌧ U U
⇤ =0BONSolution ⇤ =0BONSolution • iue3 iga for Diagram 3: Figure for Diagram 3: Figure ⇤ ⇢ R • nioa onsof points antipodal of points antipodal 1
U U =12 iue9 icws o potential toy Piecewise 9: Figure & iue4 iga R5 Diagram 4: Figure R5 Diagram 4: Figure = = iue8 icws o potential toy Piecewise 8: Figure 3.1 Sections in discussion the will by which by motivated minimum class given is wider Sitter is potential much de and toy a 3.2 for local This solution and a models. bounce with BON realistic potential analytic more approximate scalar of an toy construct a to us introduce allow to useful is model It toy A 3.3 atfrtedsrpino h O ntblt.O h te ad ihsaepotentials high-scale hand, other the On irrele- instability. of largely BON is the region regime, intermediate of satisfying (compactification) the description asymptotic in the the potential for and the vant vacuum of false shape the the between assumption, these Under treatment. satisfies quadratic. approximately remains For (3.41)). Eq. around discussion the (see ti ntbet ulaino BON. at a (3.13) vacuum of Eq. the although nucleation — 3.4, to Section decay unstable in CDL is see ordinary it will an we As admit vacuum. not Sitter does de a admits and quadratic with potentials smooth “realistic” the of for three solutions asymptotic All similar model. have toy quadratic will piecewise the potential of the utility if even solution, BON Witten potential the and 5 Figure ✏ • ✏ ⇢ ⇢ R R
instanton de Sitter Solution de Sitter Solution 1 2 iue4 iga for Diagram 4: Figure vacuum U & ⌧ ( U U ⇢ φ iue7 iga for Diagram 7: Figure
R U
U( )/U 2 fv 0 fv ilhv nipc ntebuc solution. bounce the on impact an have will ) U( )/U0 ⇠ ⇠ ⇤ ⌧ 2 2 =0BONSolution /R /R 0.0 0.5 1.0 1.5 +
m U iue3 iga for Diagram 3: Figure 0 = - : ⇤ 0.8 U =0BONSolution 1 2 iue3 iga for Diagram 3: Figure m 0 φ Left:
. U nioa onsof points antipodal U S S 2 & 2
⇢ ⇢ U ⇢/R ⇢/R U φ −∞ → 3 3 fv 1 1 1 U ⇢ = 1 iue4 iga R5 Diagram 4: Figure - = 1 2 /R M ( and and 1 0 =0 0.6 =0 φ • ✏ ✏ = . p eist aeantcal ffc nteeuto fmto for motion of equation the on effect noticeable a have to begins ) and 2 de Sitter Solution 5
h icws potential piecewise The de Sitter Solution /R . . n ewl sueti irrh ntermidro u analytic our of remainder the in hierarchy this assume will we and , Right: U ⇢ ⇢ 2 2 2 2
0 ⇢ ⇢ U M - • 5 2 . 0 0.4
iis sln as long as limits, U U
' • ' p ⇢ ⇠ U 1
ilntsaeti rpry n,i eea,tedtie features detailed the general, in and, property, this share not will ⇠ ⇠ /M
2 =0BON =0BON /R fv 2 + + /M toy U ⇠ ihpeeietypotential toy piecewise with , R R 2 p /R toy ⇠ ⇠ ( p - h quantity The /R /R fv 3 ⇢ 3 φ • 0.2 4 4 (
? 2 2 1 S 1
Near 5. Fig. of left the on depicted is ) ⇢ φ ⇢/R
< / / φ 3 1 ✏
1 p ⇢ p ) and ( =0 and =0
2 1 2 ≡ ξ ⇡ ⇡ 1 and ) 0.0 and ⇢
• ⇢/R ( . . 2 1 2 ⇢ ⇢ ⇢ ⇢ • • U U U 2 1 1 U M ⇢ U 2 and 0 fv 0 (
with , U
p ' ⇢ ( =0BON⇠ 0.2 ρ + 2 6– 16 – + mR ( φ U ihpeeietypotential toy piecewise with ,
,tgte ihmr elsi examples realistic more with together ), R ξ 4 U ⇢ 1 2 ? ,tgte ihmr elsi examples realistic more with together ), ⇠
( ,dsietesgicn ieecsi their in differences significant the despite ), m 2 /R toy φ
5 3
φ ⇢ = 4
slre hsfauegetyehne the enhances greatly feature This large. is ) 2
2 U( )/U U( )/U0 0 U( )/U0 L L/2⇡R ( 5 U φ 1.0 1.5 0.0 0.5 / p ' φ − -
( ⇠ 2 2 =0BON 0.8 ,with (3.13), Eq. from ) ⇡ + m 0.0 0.5 1.0 1.5 ssaleog that enough small is 1 φ )
2 R for for 2 ? U R - ⇢ , ⇠
• 5 2 L 0.2 1 =0BON /R - φ
q 3 ( φ φ 0.6 4 ( ⇠ 2 | 1
ξ )=2 • φ ≤ 2 3 swl prxmtdb the by approximated well is ) ≤
? / M p | per aıeysal — na¨ıvely stable appears 0 = 0.0 φ U ⌧ φ, p - 2 1 •
0.4 ⇡ ⇡ ( niae h on where point the indicates |
U( )/U
R 0 • • /M 1 /M /M ⇠ 1 5 ) fv /R | 0.2 1 exp ⇢ • p - φ ⇡ p p 1 • 0.2 ? ,teptnilis potential the 0, = ⇣ ⌧ HM • 0.4 ✏ 1 q =0 =0 HM • 0.0 φ U 3 2 1 . . ( U 2 2 /M | = 0.6 φ U M • ( | ≤ | 0 • 0.2 p p − 1 /M ⌘ ) fv 0 (3.13) . 2 p 0.8 φ M ? ? | φ p ) < ✏ 1 =0 =0 1 . . 2 2 U M 0 p The constant branch that extends toward the compactification limit (large, negative values of φ) is motivated by a property of Witten’s bubble discussed in Section 2.2, namely 3 that the quantity ρ φ0 remains constant (following Eq. (2.14)). This property also holds for a constant potential. Since the solutions we are looking for behave like Witten’s bubble near the center of the bounce, and later transition towards the false vacuum, this toy potential will allow us to construct analytic expressions for both the bounce and its Euclidean action that will be representative of more realistic models — following the discussion in Section 3.2, this toy model will provide approximate solutions for potentials that behave as U(φ) ∼ U eaφ/Mp in the compactification regime, so long as a > √6 and U M 2/R2. 0 − 0 . p 5 More intricate models can be devised that are also analytically tractable, e.g. by con- catenating other quadratic functions to construct a smooth local maximum in U(φ), and we discuss some of these possibilities and their implications in Section4. In this section, however, we focus on Eq. (3.13) because it successfully captures the features of more generic potentials most relevant for the description of BON solutions.
3.4 A bounce of nothing We now proceed to find an analytic solution to the equations of motion in the presence of the potential given in Eq. (3.13), with the boundary conditions at ξ = 0 specified in Eq. (3.3). It will be convenient to distinguish three spatial regions, according to the different behaviors of the bounce:
• A core region near the bubble wall, 0 ξ ξ , where the effects of the scalar potential ≤ ≤ 1 are subdominant and the solution behaves approximately like Witten’s BON;
• a transition region ξ ξ ξ that interpolates between the BON and false vacuum 1 ≤ ≤ 2 solutions;
• and an asymptotic region ξ ξ, far away from the bubble wall, where φ(ξ) is 2 ≤ exponentially close to the metastable de Sitter vacuum.
For the toy potential, the ξ1 that separates the core and transition regions is implicitly defined by φ(ξ1) = φ1. Within the core region, the potential U(φ φ1) = U0 is constant, 3 ≤4 and the quantity (ρ φ0) is exactly conserved as in Witten’s bubble. In the transition region, (ρ3φ ) decreases from its initial value towards ρ3φ 0. Unlike 0 0 → ξ1, the ξ2 outer boundary of the transition region is not sharply defined. Eventually φ0 becomes small enough that φ(ξ) φ is approximately constant, and the spacetime is ≈ fv approximately de Sitter. Later in this section we demonstrate that ξ (few) m 1, 2 ∼ O × − and that neither the parameter η nor the bounce action depend on the precise value of ξ2. The picture outlined in the bullet points above is sketched in Fig.6, for some potential with de Sitter radius Λ R . In the core region, 0 ξ ξ , ρ(ξ) is well approximated by 5 ≤ ≤ 1 the BON solution, which becomes approximately linear for ρ R . This ρ 1 behavior & 5 0 ' continues throughout the transition region, ξ ξ ξ , as long as ρ M / U(φ). For 1 ≤ . 2 p ξ ξ , the modulus φ(ξ) φ is nearly constant, and ρ(ξ) is well approximated by a & 2 ≈ fv p de Sitter solution.
4 2 2 Both ρ(ξ) and φ(ξ) include small perturbations, of order R5U0/Mp 1 compared to a bona-fide BON solution, as shown in Appendix A.2. In the regime ξ . R5 the corrections are negligibly small.
– 17 – 12
⇤ =12R5, ⌘ =1 10 1.5 • R 5 ⇢ 8 ⇠ + 1 de Sitter Solution 1.0 ⇢ U fv 5 '
5 =0BONSolution= 6 U U /R ⇢ /R =0BON ⇢ 0.5 U 4
⇢ ⇤ 0.0 • 2 ⌧ 2 ⇢ R -0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 • 1 & 5 ⇠/R5 0 0 5 10 15 20
⇠/R5
Figure 5: Diagram for ⇢ and ⇢ Figure 6: We compare the ρ(ξ) profile for the U = 01 BON2 (black) to a de Sitter vacuum 2 with positive U(φfv) (blue). The de Sitter length is Λ = Mp 3Mp /Ufv, and for this example we have taken η 1. The analytic model described inq Section 3.4 stitches these ' solutions together by constructing a transition region in which ρ 1 for both of these 0 ≈ geometries, corresponding to the region between ρ1 and ρ2 on the plot. Such a region exists if the de Sitter radius is large, Λ R . The precise value of ρ is determined 5 1 by solving the equations of motion for a particular potential. As we show in the inset diagram, as long as ρ1 & R5 (ξ1 & 0.4 R5), the U = 0 BON solution is well approximated ⇤ =12R by its ξ R expansion with ρ 1. The point ρ indicates whereU =0BONφ(ρ ρ ) φ 5 0 ' 2 & 2 ≈ fv 2⇡ • approaches the false vacuum. Its definition is2 somewhatp/ arbitrary, but it should satisfy 3 4 ⇢1 φ(ρ2) φfv φ1. R + − ✏ de Sitter Solution ⇠ = =0BONSolution ⇢ ' 5 U /R ⇢ U R /R ⇡
⇢ =0BON Core region U 2 L/ 3 In the region⇢2φ ⇤ φ , our choice of potential ensures that the quantity ρ φ remains • ⌧ ≤ 1 0 constant, and⇢ R set by the boundary conditions of Eq. (3.3). In this region, the CDL equations • 1 & Eqs. (2.3) and (2.4) can be rewritten as ⇠/R L(⇠)=2⇡R exp 2 /M 5 3 p ⇣ q ⌘ ⇠/R 3 R4 ρ2 ρ3φ = M R2η3, and ρ = 1 + 5 , ⇠/R (3.14) Figure0 3: Diagramp for5⇢1 and ⇢2 0 4 2 r2 s 4ρ − Λ0 Figure 6: Diagram for ⇢1 and ⇢2, with L( ) with Λ 3M 2/U . The last term inside the square root satisfies ρ Λ . This is not 0 ≡ p 0 0 an additionalq assumption on the properties of our potential, but rather a consequence of those already discussed in Section 3.3. Namely, the requirement that the overall scale of 2 2 the scalar potentialvacuum away from the asymptotic compactification regime remains . Mp /R5 allows us to neglect this last term, including3 into the regime where ξ R and ρ 1. & 5 0 ' In this core region, not far from the bubble wall, ρ(ξ) and φ(ξ) are given by Eq. (3.3) and Eq. (3.9) in the regimes where ξ R and ξ R , respectively. Appendix A.2 5 5 provides the complete solutions for ρ and φ, including the leading perturbations from the instanton
– 18 – antipodal points of S3
Figure 4: Diagram R5
2 constant U = U0 term in the equations of motion, and Appendix A.3 lists the subleading terms in the ξ R series expansions. 5 Transition region Once φ enters the quadratic region of the potential, φ φ , Eq. (2.3) can be conveniently ≥ 1 written as d2φ 3 dφ + m2φ. (3.15) dρ2 ρ dρ ' where we have taken ρ 1, which is a self-consistent assumption. The relevant solution 0 ' is given by5 K (mρ) φ(ρ) C 1 . (3.16) ' mρ
Demanding continuity of φ0 and φ across the transition allows us to obtain an expression for both C and η: 3 2 3 η R5 C Mp 2 , (3.17) ' − r2 ρ1K2(mρ1) m2R2 2e γE log η η3 5 log − , (3.18) ' 4 mρ 1 where ρ ρ(ξ ) ξ . Provided mR 1, log η 1, and we may expand the previous 1 ≡ 1 ' 1 5 expression around η = 1:
R m2R2 2e γE m2R2 1 η = 3 1 + 5 log − = 1 + 5 log + (1) . (3.19) R ' 4 mρ 4 mR O 5 1 5 In the regime mρ 1, Eq. (3.16) is of the form 2 3 R5 φ(ρ) Mp 2 (mρ 1), (3.20) ' − r2 ρ which coincides with the leading behavior of the BON solution, up to corrections suppressed by further powers of R /ρ. On the other hand, in the asymptotic regime mρ 1, we find 5 e mρ φ(ρ) − (mρ 1), (3.21) ∝ −(mρ)3/2 and the solution approaches the false vacuum exponentially fast. As anticipated in Sec- tion 3.2, ρ m 1 signals the beginning of the asymptotic regime where the bounce solution ∼ − quickly approaches the false vacuum. Motivated by the ρ m 1 functional form of φ(ρ), − we define
1 ρ (few) m− (3.22) 2 ≡ O × as the end of the transition region, where (ρ3φ ) R2M is much smaller than its initial 0 5 p value at ξ = 0.
5 K1(mρ) I1(mρ) In full generality, the solution to Eq. (3.15) is of the form φ(ρ) = c1 mρ +c2 mρ . The coefficient c2 is exponentially suppressed compared to c1, and the approximation c2 ' 0 is appropriate for our purposes.
– 19 – False vacuum region
Once ρ & ρ2, φ is exponentially close to its value in the false vacuum, and ρ is given by
ξ + γη3/2R ρ(ξ) Λ sin 5 , (3.23) ' Λ ! with γ 0.6, as defined below Eq. (2.22). Indeed, this is just the de Sitter solution of ' Eq. (2.5), shifted by (R ) towards smaller ξ. O 5 Comment on exponentially growing potentials The reader may be concerned that the approximate solution constructed above will break down in the presence of a potential that grows exponentially in the compactification limit. We now show that this is not the case, provided the rate of growth is not too close to the limiting value of the exponent, a = √6, from Section 3.1. − Recall that Eq. (2.14), using the relevant ξ = 0 boundary conditions given in Eq. (3.3), can be written as 3 ξmax ∂U M R2η3 = dξ ρ3 . (3.24) p 5 ˆ r2 − 0 ∂φ This expression allows us to estimate the change in η due to the region near the center of the bounce, where the potential is growing exponentially fast. Expanding around η = 1, and ignoring irrelevant (1) factors, the shift in η can be approximated using the small ξ O expansions of Appendix A.3, with the result
R5 1 3 ∂U δη 2 dξ ρ . (3.25) ∼ −MpR5 ˆ0 ∂φ
Combining Eq. (3.3) with the asymptotic form of the scalar potential, U U eaφ/Mp , ∼ 0 ' we find 2 U0R a δη 2 − . (3.26) ∼ Mp 1 + a/√6 In the approach to the critical value of the exponent, a √6, the small ξ expansions of → − Appendix A.3 break down, and our approximation η 1 is no longer valid. However, for ' negative values of a that are not too close to this limiting case, δη remains small provided U M 2/R2. Moreover, if U M 2m2 M 2/R2, η remains well approximated by 0 . p 5 0 . p p 5 Eq. (3.19).
3.5 de Sitter decay rate We now turn to the evaluation of the Euclidean action for the bounce solution obtained in Section 3.4. Evaluating Eq. (3.1) on a solution to the Euclidean field equations, and subtracting the action of the de Sitter false vacuum, the tunneling exponent governing the decay rate is given by
ξmax πΛ 2 2 3 2 3 2 3 ∆S = π Mp ρ(ξ) φ0(ξ) 2π dξ ρ U + 2π dξ ρdSUfv. (3.27) 3 ξ=0 − ˆ0 ˆ0 r
The first term arises after integrating by parts the φ term of Eq. (3.1), together with the restriction that ρ(ξmax) = 0, since the geometry is compact. The second term is obtained
– 20 – by evaluating the bulk terms in Eq. (3.1) on a solution to the Euclidean field equations, and the last term corresponds to minus the Euclidean action of the metastable de Sitter vacuum. The first term in Eq. (3.27) is given by
2 2 3 ∆S π Mp ρ(ξ) φ0(ξ) (3.28) ξ=0 ≡ 3 ξ=0 r = π2M 2R2η3 (3.29) p 5 3 1 = π2M 2R2 1 + m2R2 log + (m2R2) , (3.30) p 5 4 5 mR O 5 5 where in the last step we have used Eq. (3.19) to expand around η = 1. By breaking the last two terms in Eq. (3.27) into the different contributions coming from the core, transition, and false vacuum regions identified at the beginning of Section 3.4, the bounce action can be written as
2 2 3 2 ∆S = π Mp η R5 + Score + Strans + Sout. (3.31)
Let us discuss these three different terms separately. The contribution to the bounce action from the core region is given by
ξ1 ξ1+γR5 2 3 2 3 Score 2π dξ ρ(ξ) U + 2π dξ ρdS(ξ) Ufv, (3.32) ≡ − ˆ0 ˆ0 2 2 3π Mp 4 2 4 2 R5 (U0 Ufv) + (π U0R5) (3.33) ≈ − 16 φ1 − O = (π2M 2R4m2) + (π2U R4), (3.34) O p 5 O 0 5 where in the last step we have used the fact that (U U )/φ2 m2 (cf. Eq. (3.13)). 0 − fv 1 ∼ From the transition region, we find
ξ2 ξ2+γR5 2 3 2 3 Strans 2π dξ ρ U + 2π dξ ρdSUfv, (3.35) ≡ − ˆξ1 ˆξ1+γR5 3π2 1 M 2m2R4 log + (1) , (3.36) ≈ − 8 p 5 mR O 5 where the integral is dominated by the region near ρ m 1 where the potential is approx- ∼ − imately quadratic. Finally, similar to the ordinary CDL bounce solution, the contribution to the vacuum-subtracted action from the asymptotic region where the bounce is close to the false vacuum is negligibly small,
ξmax πΛ 2 3 2 3 Sout 2π dξ ρ(ξ) Ufv + 2π dξ ρdS(ξ) Ufv 0. (3.37) ≡ − ˆξ2 ˆξ2+γR5 '
Provided U M 2m2, and that m 1/R, the leading contributions to the bounce 0 . p action come from Eq. (3.30) and the log-enhanced terms in Eq. (3.36). In total,
3 1 ∆S π2M 2R2 1 + m2R2 log + (m2R2) . (3.38) ' p 5 8 5 mR O 5 5
– 21 – Eq. (3.38) provides the leading correction to the BON tunneling exponent for the nucleation of a BON within a four-dimensional de Sitter vacuum in the limit of small energy density. Unlike the CDL and HM instantons of Section 2.1, the BON action remains finite as U 0. fv → Finally, let us comment on the contribution to the bounce action from potentials that grow exponentially in the compactification limit. Using Eq. (3.24), we can write the change in the bounce action coming from the first and second terms in Eq. (3.27) as follows:
R5 M ∂U δ(∆S) 2π2 dξ ρ3 U + p . (3.39) ≈ − ˆ √6 ∂φ 0 With U U eaφ/Mp in the compactification regime, we find ∼ 0 δ(∆S) π2U R4. (3.40) ∼ − 0 5 This δ(∆S) is subdominant compared to the leading correction to the BON action in 2 2 Eq. (3.38), provided U0 . Mp m . This estimate of the change in ∆S is verified by the more careful treatment of the asymptotically exponential potential of Appendix A.3.
3.6 Application to more generic potentials 1 For ρ . m− , both φ(ξ) and ρ(ξ) match the Witten BON solutions to leading order, 1 2 2 whether U(φ) = U0 or U(φ) = Ufv + 2 m φ . This indifference to the shape of the potential in the ξ R limit extends even to exponentially growing potentials of the form U(φ) = 5 U exp(aφ/M ) for √6 < a < 0, as we discuss in Section 3.1 and Appendix A.3. Once 0 p − ρ m 1, on the other hand, the shape of the potential has a direct impact on the equations ∼ − of motion, affecting how quickly φ(ξ) approaches the false vacuum. Given two potentials with matching quadratic expansions about the false vacuum, but different φ asymptotic behaviors, there is a small-R √U /M limit in which the → −∞ 5 0 p solutions for φ(ξ) are expected to match. The value of φ at ρ m 1, ∼ − 2 2 m R5 3 φ? Mp, (3.41) ≡ − 2 r2 provides a helpful diagnostic. If U(φ) is approximately quadratic throughout the range φ? < φ < φfv, then the solutions for ρ(ξ) and φ(ξ) should be approximated by the piecewise 2 2 quadratic model, and the leading m R5 log(1/(mR5)) fractional corrections to η and the action are expected to match, even if the φ φ part of the potential is substantially ? different. Possibilities of this nature are shown in the right panel of Figure5.
4 Other Exotic Bounces
The bounce discussed in Section3 effectively concatenates Witten’s BON solution near the center of the bounce (where the effects of the scalar potential are subdominant) to the de Sitter false vacuum in the asymptotic regime. This class of bounces may be present even when the potential is monotonically increasing to the left of the false vacuum, and therefore no other instabilities are present in the compactification limit. Depending on the detailed features of the potential in the compactification regime, there may be additional saddle point solutions. For example, if U(φ) exhibits a local
– 22 – maximum to the left of φ = 0, an ordinary HM solution will also be present. If a second, deeper, minimum is present to the left of the false vacuum, then an ordinary CDL solution might also exist. For instance, this would be the case if the scalar potential features similar behavior in the compactification and decompactification regime (e.g. as indicated by the black curve in Fig.1), as might be the case in constructions where the higher-dimensional theory exhibits some equivalence between the limits of small and large radius [54]. It is not hard to see how these additional CDL or HM instantons (if present) will also have a BON counterpart. Just like the instability discussed in Section3 “stitches” together the BON and false vacuum solutions, other instantons can be built by stitching the BON to the CDL/HM solutions. The interpretation of this class of “hybrid” solutions is, however, not immediately clear: they have larger action than the other solutions, and they probably admit two negative fluctuation modes. We discuss their physical meaning further below. In this section, we briefly describe the properties of this hybrid class of instantons. For simplicity, we will restrict our discussion to a scalar potential with the qualitative features of the black curve in Fig.1, and given by the following toy model
1 1 + β Uˆ(φ) = U β e3aφ/Mp 1 eaφ/Mp + β eaφ/Mp with a > 0. (4.1) 0 0 3 − 2 1
The parameters β0 and β1 are determined by demanding that Uˆ(φ) exhibits a local mini- mum at φmin = 0, as well as a local maximum such that Uˆ(φtop) = U0. Exact expressions for β0 and β1 involve the solution of a cubic equation, with the results given in AppendixC. As an expansion in δ U /U for δ 1, the coefficients are approximated by ≡ fv 0 81 1 8δ β = + (δ), and β = + + (δ2), (4.2) 0 4 O 1 3 81 O and the position of the barrier lies at φ /M = a 1(log 3 + (δ)). The mass of the top p − − O scalar is approximately
2 ˆ 2 2 ∂ U U0 27a 19 2 m 2 2 1 + δ + (δ ) . (4.3) ≡ ∂φ ≈ Mp 2 27 O φfv
Figure7 shows Uˆ(φ) for several values of δ. It also shows Uˆ(φ) for several values of a on ˆ 2 U(L), as a function of L(φ) = 2πR5 exp 3 φ/Mp , demonstrating the impact of a on the potential barrier. q 2 2 As in Section3, we will make the assumption that U0 . Mp /R5, and we will often restrict our analysis to the thin-wall limit in which δ 1, with the goal of deriving approximate analytical expressions. Deviations from the thin-wall limit are hard to treat analytically, but they can be explored numerically, as we do in Section5.
BON-CDL Instanton When δ = U /U 1, the potential in Eq. (4.1) admits a thin-wall CDL bounce, with the fv 0 properties summarized in Section 2.1. Up to order one numbers, and for δ not too small, the value of φ inside the CDL bubble approximately satisfies
1 27 U0 1 U0 φ /M a− log a− log , (4.4) CDL p ∼ − 4 U ∼ − U fv fv
– 23 – 1.2 1.2 a =1 a =2 a =5 1.0 1.0 δ=0.9 =0.2 δ=0.7 0.8 0.8 δ=0.5 0 0 δ=0.3 0.6 /U /U 0.6 ) ) δ=0.1 ( ( 0.4 ˆ ˆ U U 0.4
0.2 0.2 • fv L =2⇡R 0.0 0.0 5 • -5 -4 -3 -2 -1 0 0.0 0.2 0.4 0.6 0.8 1.0 1.2
a /Mp L/2⇡R5
Figure 20: Potentialˆ U( ) Figure 7: Left: The scalar potential U(φ), as given in Eq. (4.1), for several values of δ U /U , as indicated in the figure. A potential of this kind, exhibiting a de Sitter ≡ fv 0 vacuum that is separated by a finite barrier from the compactification limit, admits hybrid BON-CDL and BON-HM instanton solutions. Right: To illustrate the physical meaning of a, we show Uˆ with fixed δ = 0.2 as a function of the proper length of the KK direction, 2 L = 2πR5 exp 3 φ/Mp . Larger values of a create a steeper potential near the false vacuum at L = 2qπR5, while also flattening U(L) at L πR5. a =2 ˆ ˆ so that U(φCDL) . U(0). Following [28], the surface tension is givenUfv =0 by.1 U0
0 σ = dφ 2 (U(φ) U(0)) . (4.5) 0 ˆ − δ=0 φCDL r /U )
Because of the( exponential falloff of the potential at large and negative φ, the lower limit ˆ of integrationU may be replaced by , and the integral done exactly. One finds −∞ 1 1/2 1 1/2 σ 2√6a− M U + (δ) a− M U . (4.6) ' p 0 O ∼ p 0 The bounce radius and corresponding tunneling exponent are as given in Eqs. (2.6) and 2 2 2 (2.9). For instance, in the limit of small energy density, Ufv σ /M ,ρ ¯ 4M /σ, and /Mp p ' p ∆S is given by Eq. (2.9). By analogy with the CDLFigure instanton, 21: Potential let usU define( , ⇤) φ 3 0 log η. (4.7) Mp ≡ r2 For the bounce of Section3, φ 0 up to small corrections, approximately matching the 0 ≈ 10 value of φfv of the false vacuum. A hybrid BON-CDL solution may be obtained by instead choosing φ φ . With this choice, φ approaches φ in the regime ξ R , at which 0 ' CDL CDL & 5 point the solution transitions into the familiar CDL instanton. The ratio of the BON radius to that of the KK circle in the false vacuum is now given by
√2/3 q φ q φ R3 2 0 2 CDL Ufv a η = = e 3 Mp e 3 Mp (4.8) R ' ∼ U 5 0
– 24 – where the final inequality again uses the estimate Eq. (4.4) which is only valid if δ is not too small; otherwise, it must be determined by solving the shooting problem numerically. The Euclidean action of this solution is approximately given by that of Eq. (2.9), on top of the contribution from the BON branch at ξ = 0, that is
∆S π2M 2R2η3 + ∆S (4.9) BON-CDL ' p 5 CDL √ 2 4 6φCDL 24π M 2 2 2 Mp p π Mp R5e + . (4.10) ' Ufv where in the last step we have used the expression for ∆SCDL appropriate in the limit U 0. fv → BON-HM Instanton
The ordinary HM instanton corresponds to the homogeneous solution φHM = φtop. As before, choosing BON boundary conditions at ξ = 0 with φ φ allows us to stitch 0 ' HM together the BON and HM solutions. In this case, we have
q q 2 φ0 2 φtop √2/3 3 M 3 M η = e p e p 3− a , (4.11) ' ≈ where we have used φ /M a 1 log 3, as appropriate for the δ 1 limit of our toy top p ' − − model. The Euclidean action of this solution is approximately given by
2 2 2 3 ∆SBON-HM π Mp R5η + ∆SHM (4.12) ' √ 2 2 2 6 2 4 1 1 π M R 3− a + 24π M . (4.13) ≈ p 5 p U − U fv top
Physical Interpretation of the Hybrid Solutions As can be seen from Eqs. (4.10) and (4.13), the action of the hybrid class of instantons is always larger than that of the pure BON solution discussed in Section3, as well as that of the ordinary CDL or HM bounce. Thus, these saddle points never provide the leading con- tribution to the semiclassical expansion. Furthermore, on physical grounds, it would seem that a better interpretation of these solutions is that they are already counted by (coinci- dence regions of) the product of the dilute gas sums for the BON and ordinary CDL/HM bubbles. They almost certainly possess two negative fluctuation modes, corresponding to the fact that they are best thought of as two superposed bounces rather than one. The fact that they provide an exact solution to the Euclidean equations of motion is, in this sense, something of a curiosity, a consequence of the CDL ansatz of a spherically symmet- ric instanton. On the other hand, when we analyze the space of solutions numerically, we will find that there is a limit in which the distinction between the BON and CDL/HM parts of the solutions dissolves, and the hybrid and pure BON branches of solutions match smoothly onto each other. Thus, somewhere in this regime we expect a transition in the number of negative modes.
– 25 – 5 Numerical Construction of Bounce Solutions
In this section, we perform a numerical exploration of BON bounces of four-dimensional de Sitter vacua. Our numerical analysis complements and extends the discussion of Sections 3 and4. In Section 5.1 we introduce and discuss the main features of the overshoot- undershoot algorithm that we implement to numerically explore these bounces. In Sec- tion 5.2 we present numerical results, for the bounce and its action, for potentials that fall off exponentially in the compactification regime. For concreteness, we focus on the toy model of Eq. (4.1). As discussed in Section4, this example features the hybrid BON-CDL and BON-HM solutions, that can also be resolved numerically. Finally, in Section 5.3 we focus instead on potentials that grow in the compactification regime, taking as an example the type of exponential growth that would arise if the five-dimensional theory featured a positive CC.
5.1 Preliminaries The analytic methods of Section3 treat the scalar potential as a small perturbation to the U(φ) 0 BON solution, an approach that is valid in the limit U M 2/R2. In ≡ 0 p 5 the mR 1 regime, where the solution is well-approximated by Witten’s bubble until 5 φ reaches the quadratic basin near the false vacuum, the bounce and its action can be estimated using the results of the piecewise quadratic model of Eq. (3.13). This estimate is independent of the asymptotic behavior of U(φ). For larger values of mR5, where U(φ) ceases to be approximately quadratic near φ φ , our piecewise toy model is no longer ∼ ? a valid tool for predicting the bounce action ∆S, or the value of η = R3/R5. We must instead follow a numerical approach.6 We construct numeric bounce profiles for φ and ρ using the shooting method. The initial data for ρ, φ, and their derivatives are defined at some ξ = ξinit close to zero (ξ R ). In our case the BON initial conditions are fully determined by R and η. init 5 5 Solutions for ρ(ξ) and φ(ξ) are then obtained from an iterative Runge–Kutta method on some finite grid of points (ξ , ξ + δξ, ξ + 2δξ, . . .) with some grid spacing δξ R . init init init 5 For generic values of (η, R5), the solution for φ(ξ) can be classified into one of three possible categories:
• undershoot: there is some turnaround point where φ (ξ ) = 0, where φ(ξ ) > 0 turn | turn | φ , and φ(ξ ξ ) diverges towards . There is no point at which φ(ξ) = φ ; | fv| turn −∞ fv • overshoot: the solution for φ(ξ) crosses the false vacuum at some φ(ξcross) = φfv, and diverges to + at ξ ξ ; ∞ cross • bounce: bounce solutions lie on the boundaries between undershoot and overshoot solutions on the (η, R ) plane. They approach φ (ξ) 0 as ξ ξ , where for 5 0 → → max de Sitter vacua ξmax > 0 is defined as the point where ρ(ξmax) = 0.
ξ ξmax Only the bounces are acceptable solutions, as any physical solution must satisfy φ (ξ) → 0 −−−−−→ 0. Bounce solutions can be found numerically using a recursive bisection approach, starting
6 2 2 0 If it is still the case that U0 Mp /R5, then it may be possible to apply the ρ ' 1 approximation to solve the equations of motion for a different piecewise model, one that more closely resembles U(φ)— we investigate another building block for such an approach in Appendix A.4. However, even these more 2 2 complicated analytic models cease to be reliable in regions of parameter space where U0 & Mp /R5.
– 26 – with an undershoot and an overshoot solution and testing points in the (η, R5) plane lying between the two solutions, until the interval between (η, R ) and (η, R ) becomes 5 |under 5 |over sufficiently small. Given that R5 is a property of the vacuum, we keep it fixed, and vary η to find η = R3/R5 for the bounce solution. Thus we may think of the bubble size R3 as the shooting parameter for BON geometries. As discussed in Section3 and Appendix A.3, even a potential U U exp(aφ/M ) ∼ 0 p with exponential growth in the φ limit can have a limited, even negligible effect → −∞ on the solutions for φ and ρ in the region ξ R . From Eq. (A.17) and Eq. (A.18), we 5 see that the leading potential-dependent corrections to ρ and φ are suppressed by factors of (U R2/M 2)(ξ /R )p, where p = 2 1 + a/√6 . For a > 2/3 this is small even 0 5 p init 5 − compared to the subleading (ξ/R )4/3 term in the series expansion of the U(φ) 0 Witten 5 p ≡ bubble. As a √6 the potential-dependent terms become more important, especially if 2 2 → − U0R5/Mp is not small. As discussed in Section 3.1, no solutions consistent with the ansatz of Eq. (3.2) exist when a √6. As long as we keep comfortable distance away from ≤ − a √6, however, an exponentially growing potential is no impediment to numerically ≈ − solving the equations of motion. In the remainder of this section, we explore the full set of bounce solutions in the 2 2 2 2 (η, R5) plane, from the U0 ≪ Mp /R5 limit to larger values U0 & Mp /R5. As it is this combination of dimensionful quantities that dictates whether the results are compatible with the analytic model of Section3, we define it as a single parameter
2 U0R5 u0 2 . (5.1) ≡ Mp 5.2 Exponentially falling potentials 5.2.1 Parameter space of BON solutions In this section, we focus on potentials of the form of Eq. (4.1) (depicted in Fig.7), with a single de Sitter vacuum at φfv = 0; a local maximum at φHM = φtop; and a global Minkowski vacuum at φ . This class of potentials admit CDL or HM solutions in → −∞ addition to the BON solutions we construct. When u0 & 1 we find that the BON and CDL/HM branches of solutions merge together, with ∆S approaching the action of the CDL solution, ∆S ∆SCDL. In these cases we find that there is a maximum value of ≈ (max) u0 such that potentials with u0 > u0 always lead to overshoot solutions. This implies that the corresponding KK theories are stable with respect to BON formation, unless the 7 potential hosts additional branches of bounce solutions at larger values of u0. Particularly for solutions on the BON-CDL/HM branches, it can be more convenient to express the shooting parameter η in terms of the combination φ0 = Mp 3/2 log η introduced in Eq. (4.7). This is the asymptotic value of φ(ξ) for the U(φ) 0 BON in ≡p the ξ R limit, Eq. (3.9). We can equally well think of it as our shooting parameter 5 for fixed R . In the limit u 1, the ξ R behavior of φ and ρ nearly match that of 5 0 ≪ 5 the traditional CDL bounce, with φ(ξ R ) φ and φ (ξ R ) 0. Similarly, the 5 ' 0 0 5 ≈ ξ R solution for ρ(ξ) is approximately linear, with ρ ξ + (R ). For this reason, we 5 ' O 5 expect to find a bounce solution when the shooting parameter φ0 is similar to the shooting 7For example, Ref. [55] explores potentials where the CDL thin-wall limit does not apply, which lead to multiple CDL solutions. These models would have multiple branches of BON-CDL hybrid solutions, with max various values of u0 .
– 27 – usd ahcrele h eino vrho ouin o ahptnil h regions varying The ∆ satisfying action potential. solutions each subtracted undershoot for solutions typical overshoot the + of include region “undershoots” the labeled lies the curve case each this valueOutside In Hawking–Moss the solution. on instead CDL converge no is there values smaller and For markers. diamond and with the branches In BON-CDL branch. right-hand have smaller the the on solutions case, parameter CDL latter the the vary we in and solutions; two or of zero are of there value Typically solution. Moss sstse,teBNCLhbi one approach bounces hybrid BON-CDL at the terminate satisfied, branches is BON-HM the limit this In ∆ variable shooting the for solutions the shows abscissa The cin ∆ actions input the by determined solutions bounce BON 8 Figure h rfie aeaBNlk oeisd faCLH-iergo.I diin hr salso is there addition, In region. CDL/HM-like a solutions; of of inside branch core the this BON-like labeling a “BON-CDL/HM” have nomenclature profiles the the for reason the is This solution: HM) (or CDL ordinary the for needed parameter h O eairtastosdrcl oflevcu eairfrti ls fsolutions.) of class this conditions for initial behavior vacuum false to directly transitions behavior BON the erfrt the to refer we O-D n O-Mnmnltr o yrdbuc ouin with solutions bounce hybrid for nomenclature BON-HM and BON-CDL ∞ S
φ 2 R5 U0/M 0 p swl saysltosta silt about oscillate that solutions any as well as , ≈ φ ae nteetodt ons ecnifrta nthe in that infer can we points, data two these on Based 0.0 0.5 1.0 1.5 2.0 2.5 0 orsod oteBNH/O-D rnh nti ltw fix we plot this In branch. BON-HM/BON-CDL the to corresponds - ∆
0.5 q ∼ a S φ eut o h hoigparameters. shooting the for Results : h small The . S HM 0 M orsod ote“ue O,o BNF”bac,wietelre value larger the while branch, “BON-FV” or BON, “pure” the to corresponds ≈ p 0.0 ( ntesmall the in
mR 2 ∆ φ R U /M undershoots overshoots 5 0 p 0 S ≈ CDL φ 5
0.5 q ) CDL(HM) HM FV 2 pr”BNsltosa h O-Vbac nti eto (because section this in branch BON-FV the as solutions BON “pure” 0 S log(1 S < a safnto fsotn parameter shooting of function a as HM 1.0 examples u HM / 0 0 ( φ < a / mR ii,wihfor which limit, ahBNF rnhtriae at terminates branch BON-FV Each . overshoots undershoots φ odmntaetetasto ewe O-MadBON- and BON-HM between transition the demonstrate to iue16: Figure 1.5 HM 0 iue1:Potential 17: Figure 0 5 =0 φ ( .T ac the match To (3.19). Eq. by predicted solution BON )) < u 0 0 a . → edt nesotsltos n usd hsrange this outside and solutions, undershoot to lead 0 2.0 8 0 = φ HM 1) R CDL . a a a 5 6 2.5 8– 28 – a ≈ = = =
, 1.2 1.6 2 vs 8 overshoots 0 δ ≤ 0 hw ntepo ihsldcrl,square, circle, solid with plot the on shown , φ / . φ 0 = 8 CDL 1 0 HM , 0 Left: . φ R 3.0 1 a a a sntsatisfied, not is (2.11) Eq. condition the 0, = for CdL u = = = 0 . . φ 5 0.6 0.8 1 φ ihBNH ouin aeactions have solutions BON-HM with 0 sapoiaey∆ approximately is 8 0 =0 U (or φ HM rnhsof branches 0 6= 0 n h o potential toy the and → HM ( → =0 . ahcrecrepnst aiyof family a to corresponds curve Each . 8 φ ii the limit 0 φ hw iha pncrl marker. circle open an with shown , Right: φ HM ) 0 . HM 8
φ 4
→ S U0/Mp ) 0 10 20 30 40 50 60 70 .
0 · For . nuisof units in φ φ
CDL
0 u eew hwtevacuum- the show we Here overshoots 0 ihfixed with , 0 a a → u ≥ ≥ tsmall at φ 0 1 CdL 0 ii,terneof range the limit, 0 1 1 S → φ → . . onesolutions bounce 1 hr q (2.11) Eq. where 1 HM 0 HM U / ˆ δ O solutions BON 0 2 with 0 φ ( ≡ =0 HM h Hawking– the , φ ' 0 φ fE.(4.1). Eq. of ) u δ U ( ≈ 59 . ξ 0 8 3 fv n have and , 0 = φ . ∞ → /U 2 φ CDL(HM) M 0 0 4 . < and 8 a a a a a a a a p 0 = 4 ======(5.2) HM 5 2 1.6 1.2 1.1 1 0.8 0.6 /U ) 0. → . 5 0 8 . , 14
12 δ=0.001 δ=0.2 δ=0.005 δ=0.5 10 δ=0.02 δ=0.8 2 p 8
/M a=5 0
U 6 q 5
R 4
2
0 0.5 1 2 5 R3/R5
Figure 9: The trajectoriesFigure of the 14: BON Parameter bounce solutions space R3/R in5 the (η, R5) plane are shown for a = 5 and various 0 < δ < 1. Each value of δ has a BON-FV branch of solutions, with R R (φ 0) in the u 1 limit where the potential becomes unimportant. 3 ≈ 5 0 ≈ 0 The BON-CDL branch of solutions has R < R (φ φ ) in the u 1 limit, and an 3 5 0 ≈ cdl 0 action ∆S S that is large compared to π2M 2R2. When u 1, the potential U(φ) ≈ CDL p 5 0 & is no longer a small perturbation to the Witten bubble, and we find that the BON-FV and BON-CDL branches merge together, with u0 for the bounce solution bounded from (max) (max) above by some u0 . For u0 > u0 , all solutions to the equations of motion are of the overshoot type. As δ 0 the BON-CDL action diverges as δ 1, and the unification of the → − BON-FV and BON-CDL branches is pushed to ever-larger values of u 1. 0 the (φ ,R 0) initial conditions produce overshoots. Whether the undershoot region 0 5 ≈ extends to arbitrary large u0 or whether it is capped at some finite value remains to be determined. 0 /U In Fig.8 ) we show the space of bounce instantons for Uˆ(φ) with δ Ufv/U0 = 0.8 and
( ≡ with varying values of a. In the a 1.1 cases the existence of a CDL solution is guaranteed ˆ U ≥ by Eq. (2.11), and indeed we find that the u 0 limit leads us to φ φ , shown as 0 → 0 ≈ CDL solid black markers on the left panel of the figure. With a 1.0, on the other hand, there ≤ is no CDL solution, and instead we find that φ smoothly approaches φ when u 0. 0 HM 0 → As we increase u0, we find that the existence of undershoot solutions is bounded from above, where the BON-CDL (or BON-HM) branch merges with the matching BON-FV set (max) of bounce solutions. Thus for a given potential there is an upper limit on u0, u0 u0 , a /Mp ≤ which can also be understood as an upper bound on R , R R(max), for a fixed barrier 5 5 ≤ 5 height. Figure 15: Potential U( ) The right panel of Fig.8 shows the action ∆ S as a function of φ0. As predicted in Eq. (4.9), the u 0 limit of the BON-CDL solution has an action that approaches 0 → ∆S S , or ∆S S in the BON-HM cases. Between the φ φ and φ ' CDL ' HM 0 ≈ fv 0 ≈ φCDL(HM) endpoints of the curves in Fig.8, there are bounce solutions with ∆ S larger than ∆S by some (1) fraction. CDL(HM) O In Fig.9 we show another example, with7 a = 5 kept constant and 0.001 δ 0.8 ≤ ≤ varied, and in this case we interpret the shooting parameter in terms of η = R3/R5.
– 29 – Every value of δ corresponds to a Uˆ(φ) that has a CDL solution in this case, thanks to the larger value of a. As δ 0, the BON-FV and BON-CDL branches become more → distinct, and solutions exist for larger values of u0. This is a consequence of the increasingly large difference between ∆S π2M 4/U and ∆S π2M 2R2 in the U 0 CDL ∼ p fv BON ∼ p 5 fv → limit. As a first approximation one could estimate the upper bound on u0 by setting ∆S π2(M R(max))2. This approach works well if u(max) 1, but is less accurate for CDL ∼ p 5 0 . u(max) 1, where the BON action receives large corrections. 0 5.2.2 Instanton solutions and action Field profiles
Now that we have established which values of η and R5 correspond to bounce solutions, we can inspect the solutions for ρ and φ and calculate the bounce action ∆S. 2 In Fig. 10, we show the φ(ξ) profiles for fixed u0 = 10− on the BON-FV and BON- CDL branches, for fixed a = 5 and a variety of δ values, 0.001 δ 0.999. On the ≤ ≤ BON-FV branch, in the ξ R region, φ(ξ) matches the Witten BON profile with φ 0 . 5 0 ≈ independently of the value of δ. For ξ & R5, the potential becomes important: for example, at larger δ, the solutions for φ(ξ) take longer to approach the false vacuum. For smaller values of δ . 0.5 on the BON-CDL branch, the profile for φ(ξ) has two distinct regions: an inner core, ξ R, where φ(ξ) closely matches the U = 0 Witten solution with φ φ ; . 0 ≈ CDL and a CDL-like region ξ R , with φ(ξ) φ until ξ R . In the BON-CDL plot & 5 ≈ CDL 5 we also show the CDL solution as a dashed line. For ξ & 3R5 the BON-CDL and CDL solutions match nearly exactly, particularly for δ 1. In Fig. 11, we show the solutions for ρ(ξ) for the same set of potentials U(a, δ). On the BON-FV branch the potential barrier has very little effect on ρ(ξ): it is well approximated by ρ(ξ) ρ (ξ) for ξ (few) R at small ξ, and ρ(ξ) ρ (ξ + (R )) for R ' bon . O × 5 ' dS O 5 5 . ξ ξ . This is the kind of solution shown in Fig.6, where the transition across the ≤ max potential barrier occurs in a region where ρ 1 to good approximation. The behavior on 0 ' the BON-CDL branch is more varied. Near δ 0, the solutions are well approximated by ≈ the CDL thin-wall limit. We can see from Fig. 10 that φ(ξ) φ remains constant for ≈ CDL R ξ (10)R , and that the transition from φ(ξ) φ to φ(ξ) φ occurs over a 5 . . O 5 ≈ CDL ≈ fv relatively small range of ξ. The corresponding δ 0.02 solutions for ρ are very similar to ≤ each other, each having an approximately triangular profile. In the thin-wall approximation the value ofρ ¯ 4.6 M /√U at the transition is estimated by Eq. (2.6). At larger δ 0.5, ≈ p 0 & the thin-wall assumptions are no longer as accurate. In particular, the thickness of the bubble wall increases, while the bubble itself becomes smaller. As δ 1 the solution for ρ → approaches a different well-understood limit: the sinusoidal Hawking-Moss solution defined in Eq. (2.12),
3/2 ρδ 1(ξ & η R5) ρHM(ξ + (R5)). (5.3) → ≈ O For ξ . R5, the solution for ρ at any value of δ is well approximated by the U = 0 BON ρη(ξ), Eq. (3.3), where in this case the deviation of η from 1 can be large (see Fig.9).
Action Given the numerical profiles φ(ξ) and ρ(ξ), the de Sitter-subtracted Euclidean action ∆S can be obtained from Eq.(3.27). The action obtained from the shooting profiles are shown in Fig. 12 for the Uˆ(φ) potential with a = 5. This plot shows the main results of our work.
– 30 – 2 2 BON-FV U0 =0.01 Mp /R5
a =5 p /M ) ⇠ (
BON-FV U =0.01 M 2/R2 1 0 p 5 δ=0.999 a =5 δ=0.995 ⇠/R 0.100 5 δ=0.98 p δ=0.8
/M Figure 22: BON-FV (⇠), function of ⇠/R5 ) δ=0.5 ⇠
( 0.010
δ=0.2 δ=0.02 0.001 δ=0.005 δ=0.001
10-4 0 2 4 6 8 10 12 14
⇠/R5 2.0 Figure 20: BON-FV (⇠), function of ⇠/R U =0.01 M 2/R2 5 0 p 5 δ=0.001 δ=0.02 1.5 a =5 δ=0.2
p δ=0.5
/M δ=0.8
) 1.0 ⇠
( δ=0.98 δ=0.999 0.5 BON-CDL 2 2 BON-CDL U0 =0.01 Mp /R5 CDL 0.0 0.1 0.5 1 5 10 a =550 100 ⇠/R5 p 2 2 Figure 10: Profiles of φ(Figureξ) with 23:U0 = BON-CdL 0.01 M /R (⇠on), functionthe BON-FV of ⇠/R (top) and BON-CDL /M p 5 5 )
(bottom) branches,⇠ for a = 5 and various 0 < δ < 1. For the BON-FV solutions the (
0 < ξ . R5 profiles are nearly identical: only for ξ & R5 does the shape of the potential have much impact on the profiles for φ(ξ). On the BON-CDL branch with small δ . 0.5, φ(ξ) quickly approaches φ(ξ) φ , where it lingers for a while before passing through ' CDL the barrier. For δ & 0.5 the approach to the false vacuum is more immediate. The BON- CDL φ(ξ) is shown on a log-linear plot, so that11 the BON-like (ξ . R5) and CDL-like (ξ R ) parts of the solution are simultaneously visible. For δ 0.5 the transition from 5 ≤ φ(ξ) φCDL to φ(ξ) φfv is rapid, implying the validity of the thin-wall approximation. As ' ≈ ⇠/R5 δ 1 there is no intermediate region of constant φ(ξ) φ , so the thin-wall description → ≈ CDL does not apply. Figure 21: BON-CdL (⇠), function of ⇠/R5
At small u 1, the BON-FV action is well approximated by ∆S π2M 2R2. This is the 0 . ' p 5 lowest-action BON solution for a given input R5, and is therefore the solution relevant for determining the decay rate. On the BON-CDL branch, the u 0 action is approximately 0 → 10
– 31 – δ=0.001 δ=0.005 δ=0.02 δ=0.2 δ=0.5 δ=0.8 δ=0.98 δ=0.995 δ=0.999 60 5 2 2 2 2 BON-FV U0 =0.01Mp /R5 U0 =0.01Mp /R5 BON-CDL 50 a =5 4 CDL thin a =5 wall limit 40 p p 3 1.5 /M /M 30 0 0 5 U 1.0 U 2 p p /R ⇢ ⇢ ⇢ 20 0.5
0.0 0.0 0.2 0.4 0.6 0.8 1 HM 10 limit ⇠/R5 ⇡ ⇡ 1 0 0 0 0 50 100 150 0 2 4 6 8 10
⇠pU0/Mp ⇠pU0/Mp
Figure 20: ⇢(⇠)withinset2 Left:2 BON-FV Right: BON-CdL Figure 11: Profiles for ρ(ξ) for U0R5/Mp = 0.01 on the BON-FV (left) and BON-CDL (right) branches. For the BON-FV branch the solutions for ρ(ξ) are nearly identical to the false vacuum solution ρfv(ξ) for ξ & R5. This encompasses the majority of the domain 0 ξ ξ πΛ, which is why we only show δ = 0.001 ... 0.2. For the BON- ≤ ≤ max ≈ CDL solutions, the profiles for ρ(ξ) follow the ξ R linear growth of Eq. (2.22) out to 5 ξ 10R or larger, with a relatively sharp transition from ρ 1 to ρ 1 in the case of & 5 0 ≈ 0 ≈ − small δ. As δ 1, the BON-CDL solutions for ρ(ξ) more closely match the Hawking–Moss → behavior ρ(ξ) ρHM(ξ) for all ξ & R5, even though the a = 5 CDL and HM solutions for ≈ 2 2 BON-FV φ(ξ) are distinct. U0 =0.01 Mp /R5
a =5 constant and equal to the ordinary CDL action, ∆S ∆SCDL. In the δ & 0.5 examples, p ≈ the maximum value of R5 for which there exist solutions is reasonably well approximated by
2 /M 2 2 (max) ) π M R ⇠ ∆S . For sufficiently small values of u , the field profiles and action on p 5 ( CDL 0
∼ the BON-FV branch are well approximated by the piecewise quadratic model of Section 3.3, with the shooting parameter and the action given in Eqs. (3.19) and (3.38) respectively. As demonstrated in Section 5.2.2, the BON-CDL branch solutions for ρ closely match the CDL solution, except in the ξ . R5 region where ρ(ξ) and φ(ξ) are Witten-like. As a result, the small R5 BON-CDL action is well described by ⇠/R ∆S (u 1) 5 ∆S + π2M 2R2η3, (5.4) BON-CDL 0 ≈ CDL p 5 Figure 21: BON-FV (⇠), function of ⇠/R5 where ∆SCDL is the dS-subtracted action for the CDL solution. It is clear from Fig.9 that R3 and R5 may be quite different from each other for the BON-CDL hybrid solutions: in the δ 0 thin wall limit of the a = 5 potentials, for example, R 0.5 R . → 3 ≈ 5 5.3 Exponentially growing potential: five-dimensional CC 10 In Section 3.1, we showed that the existence of BON solutions depends on the asymptotic form of the classical scalar potential. For potentials U(φ) U eaφ/Mp , we found that ∼ 0 a > √6 is a necessary constraint. In the previous section we have focused on potentials − with a > 0, and therefore asymptote to zero in the compactification limit. But the bound a > √6 is weaker, permitting even exponential growth of the potential as long as the rate −
– 32 – 105 a =5 104
10003 =0.001 4 p 10 =0.005
/M 2 =0.02 0 10
U =0.2 · 1 S 10 =0.5
=0.8 0 101 2 =0.98 2 R 5 2 p 1 M =0.995 10 ⇡ Sbon ⇡ 2 10 3 100.001 0.01 0.05 0.10 0.50 1 5 10
2 R5 U0/Mp q 4 Figure 12: The de Sitter-subtractedFigure 24: Action Euclidean as function actionof ∆SR5is shown in units of Mp /U0 as a function of R5, for a = 5 and various 0 < δ < 1. The dashed gray line shows the ∆S π2M 2R2 expectation from the small u 1 limit of the BON-FV solutions. It ≈ p 5 0 tracks the numerically-computed action closely on this branch for all u 1. This branch 0 is always of lowest action and determines the tunneling rate. On the BON-CDL branch the action approaches a constant, the ordinary CDL action ∆S SCDL, which is large 4 2 4 ≈ compared to Mp /U0 for δ . 0.5. For δ 1, SCDL 24π Mp /(δ U0), grows as 1/δ. At 2 2 2 ≈ · larger values of u0, where π Mp R5 ∆SCDL, the BON-FV and BON-CDL branches of ∼ 2 2 bounce solutions meet. (This method of estimating the value of u0 = U0R5/Mp at which 2 2 the two branches meet becomes imprecise once u0 & 1, where π Mp R5 is no longer a good approximation to the BON-FV action. is not too large.8 Now we consider cases where the potential satisfies a > √6 yet grows − exponentially in the compactification limit. The simplest examples of this type arise in the dimensional reduction of theories with positive higher-dimensional CC, c.f. Eq. (2.32). Here the asymptotic growth of the four-dimensional12 potential is exponential, but slow enough that the appropriate bubble of nothing initial conditions as ξ 0 still match Eq. (3.3). → Several example potentials of this type are shown in Fig. 13. Specifically, they corre- spond to the following modification of the Uˆ of Eq. (4.1):
CC 2 Uˆ(φ) 2 φ Λ5 Mp Uˆλ(φ) U0 + λ exp , λ . (5.5) ≡ U0 −r3 Mp !! ≡ U0 8Furthermore, a broad class of even faster-growing potentials probably admit a different type of BON solution under some plausible assumptions about the spectrum. Fluxes, for example, violate the inequality above, but can be screened by suitable dynamical sources, leading to a hybrid Schwinger production–BON decay process. In this case the BON core lies inside the region where the flux has been screened, and the modulus moves in a softer potential. A particular example of this type was studied in [40]. A general treatment of such cases should be possible along the lines of our study above, but we defer it for future work.
– 33 – 4 a =2 Ufv =0.1 U0 3 =1.00
0 =0.75
/U =0.50 ) 2
( =0.25 ˆ
U =0.10 1 =0.02
0 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 /Mp
Figure 13: The scalar potentialFigureUˆ ( 16:φ) of Potential Eq. (5.5U),( for , several⇤) values of λ ΛCCM 2/U . λ ≡ 5 p 0
The coefficient for the exponential growth, a = 2/3, is motivated by an n = 1 scenario − in which U(φ) is dominated by the five-dimensional CC, ΛCC, as in Eq. (2.32). This is small p 5 enough compared to the limiting value a = √6 that the field solutions for ξ R are − 5 closely approximated by the U = 0 Witten solutions. Here U0 is the height of the barrier in the λ 0 limit: however, when λ 1 the exponential overwhelms the Uˆ contribution, → and there ceases to be any local maximum in Uˆ (φ). Then it is clear that there is no λ=0.8 ordinary CDL solution that tunnels in the compactification direction.9 Nonetheless BON solutions can still exist. This is intuitively clear from the shooting problem in the inverted potential – there may be noovershoots point φ0 < 0 where a rolling ball can start from rest and reach the false vacuum at the origin, but it can still reach the false vacuum by shooting in from 2 p overshoots at high velocity. undershoots
−∞ /M We solve0 these cases numerically with similar methods as described in Section 5.1. Sample shootingU parameter solutions, field profiles, and actions are shown in Figs. 14, 15, q and 16, respectively.5 In Fig. 14 we provide a comparison with the analytic solution for the R shooting parameter derived in Section3. Agreement is good for small u0, as expected: the potential, despite its exponential growth, does not grow fast enough to significantly impact the shooting solution. Similarly, in Fig. 15 we see that the modulus profile continues to HM CdL match Witten’s bubble well at small ξ . R5, and Fig. 16 shows that the action continues 2 2 2 to be dominated by the π R5Mp term in the small u0 regime. 0/ HM 6 Conclusions Figure 17: R5 vs 0 for =0.8 In this work we have studied bubble of nothing instabilities of four-dimensional de Sitter vacua in theories with compactified extra dimensions. By mapping the problem of finding the relevant gravitational instanton to solving a four-dimensional CDL problem, we have been able to treat the moduli potential in a general, model-independent way, separating the 8 9Although there is probably still one in the decompactification direction; our potentials are not meant to be realistic on this side, since it does not influence the BON solutions in which we are interested.
– 34 – httedcyrt susprse ntelmto ml aumeeg s oga this as long potential (so the If energy vacuum Universe). our small in of likely seems limit and as the supersymmetry, constant, restore in cosmological also unsuppressed positive not is does with limit rate stabilization decay modulus the survive that can they that are corrections calculable small, 1 provides scale potential exponent. KK the tunneling the problem, for the When the even to irrelevant. in channels finite scale become decay energy instantons remains other largest HM which it the in or and this limit CDL a finite, either on vacuum, by is false based described the solution solutions in the density analytic to energy of (corresponding approximate vacuum of action vanishing bounce constructed Euclidean presence the have the of The The we center where observation. and the cases limit. near wall), irrelevant certain compactification bubble typically in the fact the even in solutions in is BON modulus fast potential that radial exponentially the shown the grows have of potential We stabilization scalar [30]. primarily with have in an compatible We is identified are manifold theory instabilities. internal Kaluza-Klein shrinking BON wall, of the bubble of where instability the existence case at the volume the zero on determining to focused data shrinks primary that sphere the internal are the of dimensionality the as qualitative and existence bubbles. the these determine of that criteria properties the from stabilization moduli of physics bounces. hybrid any of or function solution, a as conditions, the initial of on the of effect formulation little alternative has As potential growing limit regime. the exponentially this the the in in solution that (3.19), shooting indicating Eq. close, model, quite quadratic is piecewise the from mR expectation theoretical the parameter dimensionless of values five-dimensional a from contributions growing CC. exponentially by limit compactification the in 14 Figure ⌘ 1 ⌘ 1 λ 0.001 0.010 0.100 5 The . O hs o hnmnlgclproe,tems neetn rpriso hs bubbles these of properties interesting most the purposes, phenomenological for Thus, well as limit, compactification the in potential scalar the of behavior asymptotic The Left: 0.005 ( 1 m 2 ,i em fthe of terms in 1, R λ hoigprmtrsltoso h O-Vbac o oetasdominated potentials for branch BON-FV the on solutions parameter Shooting : 5 2 0.010 h au ftesotn variable shooting the of value The 1 = em htaentlgrtmclyenhanced. logarithmically not are that terms ) λ λ = = 1. 2. . and 0 R R ⌘ ⌘ 5 5 λ q q U U iue3:Buc solutions Bounce 35: Figure solutions Bounce 34: Figure fv fv 2 = a a 1 1 U U 0.050 =0 =0 ⇡ ⇡ =2 =2 0 0 /M /M m . oetas(elwadgen onthv n D rHM or CDL any have not do green) and (yellow potentials 0 . . 2 1 1 p p 2 2 4 1 4 1 0.100 m m U U ∼ 0 0 2 2 R R a λ 2 5 5 2 2 U mR log( log( Λ = 0 /M mR mR 5 5– 35 – 5 CC → p 2 17 5 5 M .For (4.3). Eq. of ) ) 0.500 h orcin to corrections the 1 p 2 η /U ⌘ ⌘ ssona ucinof function a as shown is ( ( 0 R R 1 ≤ 5 5 )for )for 2 .
.Tetikdse ieshows line dashed thick The 0. 0/Mp 0/Mp 2 3 4 0 1 0.0 U U Right: ( ( 0.5 √ S ) ) 1 √ u η si h rgnlBON original the in as , R R 0 u 1.0 ≈ 5 5 0 . eew show we Here q q o h aevalues same the for , r oiae by dominated are 1 U U 0 . 1.5 0 0 ,teagreement the 1, /M /M p p 2 2 2.0 √ u 0 2.5 o two for φ /R 0 the , 3.0 5 is U oiae oetas h oi n ahdlnscrepn to from correspond lines dashed and solid The 0 potentials. dominated 15 Figure praht h as aumi at curn eeat here occurring fast, is vacuum with false the to approach oe-eedn ngnrl u tlati oecsstecnrbto fcagdsources charged of is contribution rate the of the decay cases the in types some on literature the flux in the the of least of at in one impact but appeared that The general, have such 43]. in types [38–41, softened model-dependent both models manifold. be top-down of internal can explicit Examples the of potential context of allowed. the parts screened is of other either solutions growth onto Here, be most BON the off” must The compactifications. “slip process sphere or flux this collapsing exist. sources of In the dynamical still context over of may the wrapped nucleation in flux solutions active by that is become bubble to requires freedom arise potential, law of quickly could the degrees Gauss’ this too other tame where na¨ıvely if grow to example interest: that wall familiar of bubble potentials be the also even directly could near be Furthermore, solution can in BON paper naturally a this cases. curvature admit in arises from such described rather comes techniques to but The potential applied fine-tuning, in asymptotic manifold. arising internal terms a the collapsing other conditions to of the the contribution reflective boundary on dominant to not BON the important is of where equally This class cases is second potential equations. a the CDL is of the behavior there asymptotic 3.1, the Section when in seen have we limits, decompactifiaction faster. or be compactification can the process either BON in the decay CDL ordinary an admits limit. law power by . 5 epciey tsmall At respectively. 25, ( φ φ .A xetdfo h small the from expected As ). u ehd a eetne oohrcasso O ouin.Freape as example, For solutions. BON of classes other to extended be may methods Our bon u φ 0 ( ξ 10 = dse,ga)hstesoetapoc oteflevcu,a ti ie ythe by given is it as vacuum, false the to approach slowest the has gray) (dashed, ) ≈ φ h ouu rfieo h O-Vbac o xoetal rwn CC- growing exponentially for branch BON-FV the on profile modulus The : ( − ξ φ 4 ) bon , ∝ φ ( ( 1 ξ ξ /ξ cusas occurs ) nyapproaches only ) 2
iue4:Buc solutions Bounce 41: Figure (⇠)/M solutions Bounce 40: Figure (⇠)/M ahrta h xoetal damped exponentially the than rather p p 0.001 0.010 0.100 10 10 10 - - - 0.01 1 6 5 4 ξ h rfiermisBNlk,dsieterpdgot in growth rapid the despite BON-like, remains profile the ξ U U fv fv prahs2 approaches a a λ λ λ λ =0 =0 =2 =2 = = = = ξ 2.0 1.0 0.50 0.25 0.10 φ h departure the A.3, Appendix of expansions series . . fv 1 1 U U 6– 36 – nthe in
0 0
20
U
0
5
R
u
0
/M
p
⇠ =0 ⇠ R
2 =0
. . 25 25 /R /R 2 1 ξ 5 /
⌘ ⌘
5 5 .I h xmlswt larger with examples the In 3. ( (
R R U
R 0
R
5 5 5 u
0
)for )for /M 5
2
=10
ξ p
ii.The limit.
=10
10 2
O ∼
4 bon bon U U 4 ( ( ( ( ( φ ⇠ ⇠ R ) ) ( 5 ξ ,wiei h examples the in while ), ) ) ) 100 ∝ u u 0 0 e → − 10 = mξ ii indicated limit 0 nthe in − 4 and ξ u u 0 0 the R = 5 10 a =2 Ufv =0.1 U0 5 2 5 R 2 p M 2 λ=2. ⇡ λ=1. S/ 1 S ⇡2M 2R2 ⇡ p 5
0.5 0.01 0.05 0.10 0.50 1 2 R5 U0/Mp q 2 2 Figure 16: The actionFigure ∆S 38:on Action the BON-FV S/⇡ R branchas function for of exponentiallyR5 growing CC- 2 2 2 dominated potentials, shown here normalized by a factor of π Mp R5 to highlight the 2 2 2 leading deviations. At small u0 . 1, the action remains well approximated by π R5Mp , with corrections to ∆S that become large as mR5 & 1. is subleading [38]. Bottom-up methods like those used in this work could also be adapted to these cases, most easily in the point charge approximation [40], allowing the study of BON in flux compactifications while remaining agnostic about the precise form and origin of the potential in the vicinity of the false vacuum. This is an intriguing direction for future (⇠) work. bon What can we conclude about our universe? From the decay rate estimates and the
age of the universe we canp infer, for example, an upper bound on the size of compactified 1
dimensions [44]: for a collapsing/M internal S , a decay rate similar to that of Witten’s bubble )
⇠ 1 implies a lower bound( on the size of the KK circle, 2πR5 & 50Mp− . More generally, the presence of rapid BON instabilities can potentially shrink the region of the string landscape that is capable of giving rise to a vacuuma =2 compatible with our universe. Ufv =0.1 U0 R2U =10 4M 2 Acknowledgments 5 0 p
We thank Nathaniel Craig, Michael Dine, Paddy Fox, Arthur Hebecker, and Matt Reece for comments and discussions. PD and BL acknowledge⇠/R5 support from the US Department of Energy under Grant No. DE-SC0015655. The research of IGG is funded by the Gordon and Betty Moore FoundationFigure 39: through Bounce Grant solutions GBMF7392,⌘(R5)forU and ( by) the National Science Foundation through Grant No. NSF PHY-1748958.
A Solutions to Equations of Motion 19 In this Appendix we provide several exact and semi-exact analytic solutions to the equations of motion. Section A.1 reproduces the exact solutions for ρ(ξ) and φ(ξ) for the U = 0 Witten bubble [30]. In Section A.2 we show how the solutions for ρ(ξ) and φ(ξ) are modified in the presence of a constant potential, U = U0. This situation arises in the
– 37 – φ φ region of the toy potential in Section 3.3, shown in Figure5. Lastly, in Section A.3 ≤ 1 we provide the small ξ R solutions for ρ and φ in the presence of a potential with the 5 asymptotic form U(φ) U exp(aφ/M ) in the compactification limit. These corrections ∼ 0 p are particularly relevant to the numeric calculation of the bounce solutions, which depend sensitively on the initial conditions defined at ξ R . init 5 A.1 Witten Bubble of Nothing 3 In the absence of a potential, the quantity (ρ φ0) is exactly conserved, and the equations of motion for φ and ρ can be expressed as
3 2 2 3 3 3 2 2 η R5 ρ φ0 = M η R , ρ0 = 1 + , (A.1) p 2 5 2ρ2 r 3 2 where from the perspective of the four dimensional equations of motion, (η R5) is simply 3 an integration constant associated with (ρ φ0). Both equations can be integrated, using dρ dξ = , (A.2) 3 2 2 η R5 + 1 + 2ρ2 r where the + sign is appropriate for initial conditions with ρ0 > 0. By inspecting the ρ0 2 2 equation of motion from Eq. (A.1), it is clear that ρ0 is bounded from below by ρ0 = 1, meaning that ρ cannot change sign for 0 ρ < . 0 ≤ ∞ Eq. (A.2) can be integrated to produce an exact solution for ξ(ρ). Imposing ρ(0) = 0 so that there is no conical singularity at ξ = 0,
2ρ3 1 , 3 4ρ4 ξ(ρ) = F 2 4 . (A.3) 3η3R2 2 1 7 − η6R4 5 4 5
In the ρ R or ρ R limits, the series expansions for ξ(ρ) can be inverted to produce 5 5 ρ(ξ), which is often more useful. Rather than finding φ(ξ), the simpler analytic expression is φ(ρ), which we derive by integrating
M dρ 3 η3R2 2 φ φ η3R2 dφ = p 5 , − 0 = arcsinh 5 , (A.4) 6 4 3 2 η R5 2 ρ 3 Mp − 2ρ 1 + 4ρ4 r r q where we introduce a second integration constant,
3 φ0 Mp log η. (A.5) ≡ r2 Series Expansions: Near the center of the bubble core, ξ R, the function ξ(ρ) can be inverted through its series expansion,
3ξ 1/3 2η 2 3ξ 4/3 30η 4 3ξ 8/3 ρ(ξ R ) = ηR 1 + − − + (ξ4/R4) , 5 5 2R 7 2R − 539 2R O 5 5 5 5 ! (A.6)
– 38 – and φ(ρ) can be expressed as a function of ξ,
2 3ξ 9η 2 3ξ 4/3 531η 4 3ξ 8/3 φ(ξ R ) = M log − + − + (ξ4/R4) . 5 p 3 2R − 14 2R 2156 2R O 5 r 5 5 5 ! (A.7)
4/3 As a consequence of the singular ξ = 0 initial conditions, the (ξ/R5) corrections produce (1) terms in the equations of motion even in the ξ 0 limit. O → Well away from the bubble center, ξ R , ρ(ξ) and φ(ξ) are approximated by 5 3 2 6 4 3/2 Γ 4 η R5 5 4 ρ(ξ R5) = ξ + η R5 + 3 + (R /ξ ) (A.8) √ 2π 8ξ O 2 3 η3R2 Γ 3 η9/2R3 φ(ξ R ) = φ M 5 4 5 + (R4/ξ4) . (A.9) 5 0 p 2 3 5 − r2 2ξ − √2πξ O ! As ξ , ρ ξ grows linearly, while φ φ . → ∞ ≈ → 0 A.2 Constant Potential 3 As established in Eq. (2.14), (ρ φ0) is conserved if the potential is constant, even if U(φ) is nonzero,
U(φ) = U0. (A.10) The equations of motion can be solved with BON initial conditions, this time with
3 2 2 2 2 η R5 ρ dρ ρ0 = 1 + 2 2 , dξ = , (A.11) 2ρ − Λ 3 2 2 2 0 η R5 ρ + 1 + 2ρ2 2 − Λ0 r where in analogy with Λ for the de Sitter vacuum we have defined
2 3Mp Λ0 . (A.12) ≡ s U0 The exact solution to the equations of motion can be expressed in terms of elliptic functions. It can also be expanded in powers of R5/Λ0, which is helpful for capturing the leading perturbations to ξ(ρ) and φ(ρ) in the limit ρ `. To first order in (R2/Λ2) and (ρ2/Λ2), 5 0 0 8ρ4 3 1 3 4 3 2 1 + 6 4 1 1 4 2ρ 2 4 4ρ 5η R5 5η R5 4 2 4ρ ξ 2F1 7 + ρ 2F1 5 3 2 6 4 2 4ρ4 6 4 ' 3η R5 4 − η R5 24Λ0 − 4 − η R5 1 + 6 4 η R5 q (A.13) 1/2 2 φ φ η3R2 η3R2 2ρ2 ρ2 4ρ4 − − 0 arcsinh 5 + 5 arcsinh 1 + . (A.14) 3 M ' − 2ρ2 4Λ2 η3R2 − 2Λ2 η6R4 r p 0 5 0 5 The series expansions in ρ/Λ 1 can be extended to arbitrarily higher order by inte- 0 grating the series expansion of Eq. (A.11). For R Λ , the ρ R behavior of the 5 0 5 solutions is approximately de Sitter, with ρ 1 ρ2/`2. Only if R Λ is it necessary 0 ≈ − 5 ∼ 0 to abandon the series expansions in favor of the elliptic function solution for ξ(ρ). p
– 39 – A.3 Small ξ Limits Given some nonzero potential, the small ξ expansions for ρ(ξ) and φ(ξ) generically include potential-dependent terms not present in Eq. (A.6) and Eq. (A.7). This is clear from Eq. (A.2). More generically, for a potential with φ asymptotic behavior → −∞ U(φ) U exp(aφ/M ) (A.15) ≈ 0 p with a > √6 (here we consider only the case of an internal shrinking S1, or n = 1 in the − notation of Section 3.1), approximate solutions for ρ and φ can be derived by linearizing the equations of motion in the ξ R5 limit. The leading U0 dependent perturbations p enter the series expansion at order (ξ/R5) with a p 2 1 + , (A.16) ≡ √6 with the solutions 3ξ 1/3 2η 2 3ξ 4/3 30η 4 3ξ 8/3 ρ(ξ R ) ηR 1 + − − 5 ' 5 2R 7 2R − 539 2R 5 5 5 4(1 p) R2U 3ξ p + − 5 0 , (A.17) 9p(p + 1) M 2 2R p 5 ! 2 3ξ 9η 2 3ξ 4/3 531η 4 3ξ 8/3 φ(ξ R ) M log − + − 5 ' p 3 2R − 14 2R 2156 2R r 5 5 5 2(3 p) R2U 3ξ p + − 5 0 . (A.18) 3p(p + 1) M 2 2R p 5 ! The effect of the potential can be treated as a small perturbation in the ξ R limit 5 as long as a > √6 (p > 0). In the small R limit, R2U M 2, the impact of the − 5 5 0 p potential is further suppressed, which is why we have terminated the series expansion at 2 2 p 8/3 R5U0Mp− (ξ/R5) while keeping the (ξ/R5) contribution from the small ξ expansion of the Witten solution. As the point-and-shoot numeric solution is highly sensitive to the 8/3 initial conditions at ξ = ξinit, the (ξ/R5) level of precision is needed in some cases even when ξinit R5. 2 2 Working at linear order in R5U0/Mp , the small ξ solutions for ρ and φ are easily extended to potentials of the form U(φ) U exp(a φ/M ), (A.19) ≈ j j p Xj if the asymptotic form of the potential includes multiple terms with similar aj. In this case, the potential-dependent corrections to the Witten solutions for ρ(ξ) and φ(ξ) can be collected into
1/3 2 pj 3ξ 4(1 pj) R5Uj 3ξ δρ(ξ R5) ηR5 − 2 , (A.20) ' 2R5 9pj(pj + 1) Mp 2R5 Xj 2 2(3 p ) R2U 3ξ pj δφ(ξ R ) M − j 5 j (A.21) 5 ' p 3 3p (p + 1) M 2 2R r j j j p 5 X
– 40 – where
aj pj 2 1 + (A.22) ≡ √6 for each term in the sum.
Perturbation to η
At the same level of precision, it is possible to write the leading U0 dependence as a multiplicative correction to the BON solutions for ρ and φ0: 4(1 p) R2U 3ξ p ρ(ξ) ρ(η) (ξ) 1 + − 5 0 (A.23) ≈ bon 9p(p + 1) M 2 2R p 5 2 p 2(3 p) R5U0 3ξ φ0(ξ) φ0 (ξ) 1 + − (A.24) ≈ bon 3(p + 1) M 2 2R p 5 2 p 3 3 3 2 4 2p R5U0 3ξ ρ φ0 M η R 1 + − . (A.25) ≈ 2 p 5 3p M 2 2R r p 5 This is a convenient way to approximate the correction to η caused by the presence of the exponentially growing potential in the ξ R limit. In the analytic model of Fig.5 with 5 U = U for φ φ , the constant potential guarantees that (ρ3φ ) should be conserved. 0 ≤ 1 0 Incidentally, U = U0 corresponds to the p = 2 special case where the leading small ξ 3 correction to ρ φ0 vanishes, as it must. 3 Given a particular value of (ρ φ0) at φ = φ1, and imposing continuity in φ(ξ) and φ0(ξ), the solution for φ(ξ ξ ) is given in general by the linear combination of Bessel functions ≥ 1 K (mρ) I (mρ) φ(ρ ρ ) c 1 + c 1 . (A.26) ≥ 1 ' 1 mρ 2 mρ
If c2 is nonzero, then whether this solution is an overshoot or undershoot can be read off from the sign of c2. For the bounce solution, c2 must be exponentially suppressed. Imposing c 0 as a requirement on the initial conditions, the value of η(R ) for a given 2 ' 5 R5 must satisfy
2 (ρ3φ ) m2 2e γE log η 0 log − (A.27) ' 3 Mp ! 4 mρ1 r ξ=ξ1 2 p 2 4 2p R U0 3ξ1 m 1 η 1 η3R2 1 + − 5 log + (1) (A.28) − ≈ 5 3p M 2 2R 4 mR O p 5 5 3 as in Eq. (3.18), but where we have left the solution for η in terms of the value of (ρ φ0) at ξ . The η 1 expansion is still appropriate, with the leading terms 1 ' 2 2 m R5 1 2 2 η 1 + log + (m R5) (A.29) ≈ 4 mR5 O 2 p unchanged, assuming the R5U0(ξ1/R5) term remains small. This is achieved when ξ1 . 2R5/3, which in turn requires φ1 & 0.36 Mp. In this limit the transition from the BON-like initial condition to the approximately quadratic φ φ regime occurs promptly. If instead ∼ fv ξ1 & 2R5/3, then there may be an intermediate region where φ(ξ) is neither described by the Witten BON nor the Bessel function solution of Eq. (A.26), and the prediction for (η 1) is only reliable if R2U /M 2 is very small. − 5 0 p
– 41 – A.4 Piecewise Linear Model The analytic model of Section 3.4 is not the only potential with tractable equations of motion. If the potential is approximately linear in the region where ρ 1, analytic 0 ' solutions for ρ, φ and η can be derived in a similar manner as in the piecewise quadratic 2 model of Figure5. Models of this form are relevant especially at larger values of R5U0, where the potential becomes important before φ(ξ) is close enough to the false vacuum that U(φ) is approximately quadratic. Labeling the point where the solution for φ deviates from the Witten solution as φ?, in analogy with Eq. (3.41), we can expect a linear model to be particularly useful when φ? happens to fall in the neighborhood of an inflection point in U(φ). A characteristic example of the linear model interpolates between U(φ1) = U0 and U(φfv) = Ufv with constant slope: for example, U(φ) = U exp(a(φ φ )/M ) φ φ 0 − 1 p ≤ 1 U(φ) = U + (∆U/∆φ)(φ φ ) φ φ φ , (A.30) 0 − 1 1 ≤ ≤ fv U(φ) = U φ φ fv fv ≤ where ∆U U U , ∆φ φ φ . (A.31) ≡ 0 − fv ≡ fv − 1 As long as a √6 and R2U M 2, Eq. (A.25) indicates that (ρ3φ ) will remain approx- & − 5 0 . p 0 imately constant for 0 ξ R . For simplicity we specialize to the a = 0 case, where ≤ . 5 U(φ φ ) = U remains exactly constant, and (ρ3φ ) is conserved, ≤ 1 0 0 6 4 2 3 3 3 2 η R5 ρ ρ φ0 = Mp η R5, ρ0 = 1 + 4 2 , (A.32) r2 s 4ρ − Λ0 where for convenience we define 2 2 3Mp 3Mp Λ0 , Λ . (A.33) ≡ s U0 ≡ s Ufv 2 2 Up to corrections of order R5/Λ0, ρ and φ are approximately equal to the U = 0 Witten solutions. The values of ξ and ρ = ρ(ξ ) defined by φ(ξ ) φ can be obtained from 1 1 1 1 ≡ 1 Eq. (A.13) If (φ φ ) M , then ρ R , and ξ can be determined more simply from fv − 1 . p 1 & 5 1 the large ξ expansions of Eq. (A.13). The primary differences between the linear and quadratic models appear, of course, in the transition region ξ ξ ξ , where U(φ) and (ρ3φ ) change in value. If ρ R and 1 ≤ ≤ 2 0 1 & 5 ρ Λ , then the solution for ρ throughout the entire transition is well described by 2 0 3/2 ρ0 1, ρ(ξ ξ ξ ) ξ + γη R, (A.34) ' 1 ≤ ≤ 2 ≈ for the γ 0.6 defined in Eq. (2.23). Here ξ refers to the end of the transition, where ≈ 2 φ(ξ ) = φ . Unlike the quadratic model, the interval ξ ξ is precisely defined, because 2 fv 2 − 1 φ(ξ) approaches the false vacuum in finite time. The value of ξ2 is found by integrating φ0: (ξ ξ )(ξ ξ + 2ρ ) 3 φ(ξ ξ ξ ) φ + − 1 − 1 1 4 M η3R2 1 ≤ ≤ 2 ' 1 8ρ2(ξ ξ + ρ )2 2 p 5 1 − 1 1 r ∆U ρ2(ξ ξ )(ξ ξ + 2ρ ) , (A.35) − ∆φ 1 − 1 − 1 1 φ(ξ ξ ξ ) = φ . (A.36) 2 ≤ ≤ max fv
– 42 – For the linear model, the bounce condition Eq. (2.14) becomes quite simple:
ξ=ξ ξ2 3 2 ∆U 3 ρ φ0 = dξ ρ (A.37) ξ=ξ1 ∆φ ˆξ 1 3 ∆U ρ2 ∆U ρ4 ρ4 M η 3R2 − dρ ρ3 = − 2 − 1 . (A.38) 2 p 5 ' ∆φ ˆ ∆φ 4 r ρ1 Assuming a nonzero value of U , the value of (ρ3φ ) must vanish as ξ ξ . Otherwise, φ(ξ) fv 0 → 2 would continue moving in the positive direction with constant U(φ), while ρ Λ. Once → ρ decreases from its maximum, the value of φ0 would begin to increase again, ultimately approaching φ as ρ 0. 0 → ∞ → In the U R2 M 2 limit (R Λ ), the solution for η is approximately 0 5 p 5 0 1/2 η3/2R 1 3 M U log η 5 p 1 fv + (R /Λ ) (A.39) Λ 2 2 ∆φ U 5 0 ≈ 0 r ! r − 0 O ! 1/2 R2 U U η 1 + 5 0 − fv . (A.40) ≈ 2√6M ∆φ p ! Unlike Eq. (3.19) from the quadratic model, where (η 1) m2R2 log(mR ), this correction − ∼ 5 5 to η 1 is linear in R5. Calculating the bounce action ∆S, we find that the leading correction − 2 2 2 to ∆S = π Mp R5 is also quadratic in R5: ∆φ ∆S π2M 2η3R2 1 √6 + U R2/M 2 (A.41) ≈ p 5 − M O 0 5 p p ∆φ 3 3/4 R U U ∆S π2M 2R2 1 √6 + 5 0 − fv . (A.42) ≈ p 5 − M 2 M ∆φ/M p p s p !
So, the corrections to η and ∆S in the small R5 limit are both large compared to the analogous terms of the piecewise quadratic model. The linear model may be useful in cases where φ(ξ) diverges from φbon at values of φ? far away from the false vacuum. If U(φ φ?) is approximately linear in this region, then 2 2 2≈ we anticipate corrections to the π Mp R5 action that are quadratic in R5, as in Eq. (A.42). The fact that the U (linear) equations of motion can be integrated when ρ 1 also ∼ 0 ' permits the assembly of complicated piecewise models, composed of linear and quadratic sections, that can more accurately reconstruct a given U(φ). As long as ρ 1 throughout 0 ' the transition across the potential barrier, solutions for φ and ρ can be found by stitching together consecutive instances of Eq. (A.26) and Eq. (A.35).
B Analytic Continuation
2 2 Expanding the line element of the unit three-sphere as dΩ3 = dψ + sin(ψ) dΩ2, the O(4)- symmetric ansatz of Eq.(2.2) can be written as 2 2 2 2 dsE = dξ + ρ(ξ) (dψ + sin(ψ) dΩ2). (B.1) Performing the analytic continuation ψ π/2 + iχ, with χ R, we obtain the Lorentzian- → ∈ signature metric ds = ρ(ξ)2dχ2 + dξ2 + ρ(ξ)2 cosh(χ)2dΩ . (B.2) − 2
– 43 – The plane of symmetry of the metric given by ψ = π/2 plays the role of tE = 0, and the continuation ψ π/2 + iχ is equivalent to t it. → E → For Witten’s BON, in the regime r R , we have ρ(ξ) r ξ + γR , and the 5 ' ' 5 previous equations reads
ds = r2dχ2 + dr2 + r2 cosh(χ)2dΩ (B.3) − 2 = dt˜2 + dr˜2 +r ˜2dΩ , (B.4) − 2 where in the second step we have used the new coordinates:
r˜ r cosh(χ), ≡ (B.5) (t˜ r sinh(χ). ≡ As discussed in [30], the BON can be regarded as a kind of distorted Minkowski spacetime, subject to the restriction r2 =r ˜2 t˜2 R2. − ≥ 5 Our BON solution, on the other hand, asymptotes to the 4D de Sitter vacuum, as given in Eq.(3.23). The appropriate change of coordinates in this asymptotic regime is of the form r˜ Λ sin ξ+γR5 cosh(χ), ≡ Λ (B.6) t˜ Λ arctanh tan ξ+γR5 sinh(χ) . ≡ Λ h i Eq.(B.2) then reads 1 r˜2 r˜2 − ds 1 dt˜2 + 1 dr˜2 +r ˜2dΩ , (B.7) ' − − Λ2 − Λ2 2 which indeed corresponds to the static description of de Sitter spacetime.
C Specific Potential
For the numeric evaluation of the equations of motion, we specialize to potentials of the form
j aφ/Mp U(φ) = U0 bje , (C.1) Xj where the parameter a controls the overall scaling of the potential with φ, and the coeffi- cients bj are such that there is a local minimum U = Ufv > 0 at φfv, and a local maximum of U = UHM at some φHM < φfv. The simplest potentials in this family require only three nonzero bj, for example j = 1, 2, 3. The action is invariant with respect to shifts in φ, and we use this freedom to set φfv = 0. The equations of motion also do not depend on the overall normalization of U(φ): a factor of U0 can be absorbed into a redefinition of the coordinates ξ and ρ, which rescales the action by an overall constant:
ξ U /M 2 ξ, ρ U /M 2 ρ, S (M 4/U ) S. (C.2) → 0 p → 0 p → p 0 q q
– 44 – In this way the four independent parameters (a, b1, b2, b3) can be reduced to two parametric degrees of freedom: the original a, and a parameter δ defined as U δ fv . (C.3) ≡ U0
The coefficients bj in this scheme are given by
U(φ) 1 + β 1 = β β eaφ 1 e2aφ + e3aφ , (C.4) U 0 1 − 2 3 0 where β0 and β1 satisfy 3(9 5δ) 9 β = − + (3 + √δ)(1 √δ)2/3(1 + √δ)1/3 + (3 √δ)(1 + √δ)2/3(1 √δ)1/3 , 0 4 4 − − − 2δ 1 β1 = + , (C.5) β0 3 for 0 δ 1. The solutions for β and β can in principle be extended simply to δ < 0 ≤ ≤ 0 1 via √δ i√ δ. → − The location of the center of the barrier, φtop, is given by
a φtop = log(β1), (C.6) Mp and the Hawking–Moss solution is simply φ(ξ) = φtop.
References
[1] Supernova Cosmology Project collaboration, S. Perlmutter et al., Measurements of Ω and Λ from 42 high redshift supernovae, Astrophys. J. 517 (1999) 565–586, [astro-ph/9812133]. [2] SDSS collaboration, W. J. Percival et al., Baryon Acoustic Oscillations in the Sloan Digital Sky Survey Data Release 7 Galaxy Sample, Mon. Not. Roy. Astron. Soc. 401 (2010) 2148–2168,[0907.1660]. [3] Planck collaboration, N. Aghanim et al., Planck 2018 results. VI. Cosmological parameters, Astron. Astrophys. 641 (2020) A6,[1807.06209]. [4] S. Kachru, R. Kallosh, A. D. Linde and S. P. Trivedi, De Sitter vacua in string theory, Phys. Rev. D 68 (2003) 046005,[hep-th/0301240]. [5] V. Balasubramanian, P. Berglund, J. P. Conlon and F. Quevedo, Systematics of moduli stabilisation in Calabi-Yau flux compactifications, JHEP 03 (2005) 007,[hep-th/0502058]. [6] J. M. Maldacena and C. Nunez, Supergravity description of field theories on curved manifolds and a no go theorem, Int. J. Mod. Phys. A 16 (2001) 822–855,[hep-th/0007018]. [7] P. K. Townsend, Cosmic acceleration and M theory, in 14th International Congress on Mathematical Physics, 8, 2003. hep-th/0308149. [8] M. P. Hertzberg, S. Kachru, W. Taylor and M. Tegmark, Inflationary Constraints on Type IIA String Theory, JHEP 12 (2007) 095,[0711.2512]. [9] T. Wrase and M. Zagermann, On Classical de Sitter Vacua in String Theory, Fortsch. Phys. 58 (2010) 906–910,[1003.0029].
– 45 – [10] G. Shiu and Y. Sumitomo, Stability Constraints on Classical de Sitter Vacua, JHEP 09 (2011) 052,[1107.2925]. [11] I. Bena, M. Grana, S. Kuperstein and S. Massai, Giant Tachyons in the Landscape, JHEP 02 (2015) 146,[1410.7776]. [12] D. Kutasov, T. Maxfield, I. Melnikov and S. Sethi, Constraining de Sitter Space in String Theory, Phys. Rev. Lett. 115 (2015) 071305,[1504.00056]. [13] D. Andriot and J. Bl˚ab¨ack, Refining the boundaries of the classical de Sitter landscape, JHEP 03 (2017) 102,[1609.00385]. [14] J. Moritz, A. Retolaza and A. Westphal, Toward de Sitter space from ten dimensions, Phys. Rev. D 97 (2018) 046010,[1707.08678]. [15] S. Sethi, Supersymmetry Breaking by Fluxes, JHEP 10 (2018) 022,[1709.03554]. [16] U. H. Danielsson and T. Van Riet, What if string theory has no de Sitter vacua?, Int. J. Mod. Phys. D 27 (2018) 1830007,[1804.01120]. [17] Y. Hamada, A. Hebecker, G. Shiu and P. Soler, On brane gaugino condensates in 10d, JHEP 04 (2019) 008,[1812.06097]. [18] Y. Hamada, A. Hebecker, G. Shiu and P. Soler, Understanding KKLT from a 10d perspective, JHEP 06 (2019) 019,[1902.01410]. [19] F. Carta, J. Moritz and A. Westphal, Gaugino condensation and small uplifts in KKLT, JHEP 08 (2019) 141,[1902.01412]. [20] X. Gao, A. Hebecker and D. Junghans, Control issues of KKLT, 2009.03914. [21] G. Obied, H. Ooguri, L. Spodyneiko and C. Vafa, De Sitter Space and the Swampland, 1806.08362. [22] S. K. Garg and C. Krishnan, Bounds on Slow Roll and the de Sitter Swampland, JHEP 11 (2019) 075,[1807.05193]. [23] H. Ooguri, E. Palti, G. Shiu and C. Vafa, Distance and de Sitter Conjectures on the Swampland, Phys. Lett. B 788 (2019) 180–184,[1810.05506]. [24] A. Bedroya and C. Vafa 1909.11063. [25] M. Dine and N. Seiberg, Is the Superstring Weakly Coupled?, Phys. Lett. B 162 (1985) 299–302. [26] S. B. Giddings, The Fate of four-dimensions, Phys. Rev. D 68 (2003) 026006, [hep-th/0303031]. [27] S. B. Giddings and R. C. Myers, Spontaneous decompactification, Phys. Rev. D 70 (2004) 046005,[hep-th/0404220]. [28] S. R. Coleman and F. De Luccia, Gravitational Effects on and of Vacuum Decay, Phys. Rev. D 21 (1980) 3305. [29] S. Hawking and I. Moss, Supercooled Phase Transitions in the Very Early Universe, Adv. Ser. Astrophys. Cosmol. 3 (1987) 154–157. [30] E. Witten, Instability of the Kaluza-Klein Vacuum, Nucl. Phys. B 195 (1982) 481–492. [31] M. Dine, P. J. Fox and E. Gorbatov, Catastrophic decays of compactified space-times, JHEP 09 (2004) 037,[hep-th/0405190]. [32] M. Fabinger and P. Horava, Casimir effect between world branes in heterotic M theory, Nucl. Phys. B 580 (2000) 243–263,[hep-th/0002073]. [33] D. Brill and G. T. Horowitz, Negative energy in string theory, Phys. Lett. B 262 (1991) 437–443.
– 46 – [34] S. P. De Alwis and A. T. Flournoy, Closed string tachyons and semiclassical instabilities, Phys. Rev. D 66 (2002) 026005,[hep-th/0201185]. [35] B. S. Acharya, Supersymmetry, Ricci Flat Manifolds and the String Landscape, JHEP 08 (2020) 128,[1906.06886]. [36] I. Garcia Etxebarria, M. Montero, K. Sousa and I. Valenzuela, Nothing is certain in string compactifications, JHEP 12 (2020) 032,[2005.06494]. [37] G. Dibitetto, N. Petri and M. Schillo, Nothing really matters, JHEP 08 (2020) 040, [2002.01764]. [38] G. T. Horowitz, J. Orgera and J. Polchinski, Nonperturbative Instability of AdS(5) x S**5/Z(k), Phys. Rev. D 77 (2008) 024004,[0709.4262]. [39] J. J. Blanco-Pillado and B. Shlaer, Bubbles of Nothing in Flux Compactifications, Phys. Rev. D 82 (2010) 086015,[1002.4408]. [40] A. R. Brown and A. Dahlen, Bubbles of Nothing and the Fastest Decay in the Landscape, Phys. Rev. D 84 (2011) 043518,[1010.5240]. [41] J. J. Blanco-Pillado, H. S. Ramadhan and B. Shlaer, Decay of flux vacua to nothing, JCAP 10 (2010) 029,[1009.0753]. [42] J. J. Blanco-Pillado, B. Shlaer, K. Sousa and J. Urrestilla, Bubbles of Nothing and Supersymmetric Compactifications, JCAP 10 (2016) 002,[1606.03095]. [43] H. Ooguri and L. Spodyneiko, New Kaluza-Klein instantons and the decay of AdS vacua, Phys. Rev. D 96 (2017) 026016,[1703.03105]. [44] P. Draper, I. Garcia Garcia and B. Lillard, Bounce Of Nothing, 2105.08068. [45] S. R. Coleman, The Fate of the False Vacuum. 1. Semiclassical Theory, Phys. Rev. D 15 (1977) 2929–2936. [46] C. G. Callan, Jr. and S. R. Coleman, The Fate of the False Vacuum. 2. First Quantum Corrections, Phys. Rev. D 16 (1977) 1762–1768. [47] A. H. Guth and E. J. Weinberg, Could the Universe Have Recovered from a Slow First Order Phase Transition?, Nucl. Phys. B 212 (1983) 321–364. [48] L. G. Jensen and P. J. Steinhardt, Bubble Nucleation and the Coleman-Weinberg Model, Nucl. Phys. B 237 (1984) 176–188. [49] L. G. Jensen and P. J. Steinhardt, Bubble Nucleation for Flat Potential Barriers, Nucl. Phys. B 317 (1989) 693–705. [50] J. C. Hackworth and E. J. Weinberg, Oscillating bounce solutions and vacuum tunneling in de Sitter spacetime, Phys. Rev. D 71 (2005) 044014,[hep-th/0410142]. [51] P. Batra and M. Kleban, Transitions Between de Sitter Minima, Phys. Rev. D 76 (2007) 103510,[hep-th/0612083]. [52] P. Draper and S. Farkas, Transplanckian Censorship and the Local Swampland Distance Conjecture, JHEP 01 (2020) 133,[1910.04804]. [53] A. R. Brown and A. Dahlen, On ’nothing’ as an infinitely negatively curved spacetime, Phys. Rev. D85 (2012) 104026,[1111.0301]. [54] C. Angelantonj, M. Cardella and N. Irges, An Alternative for Moduli Stabilisation, Phys. Lett. B 641 (2006) 474–480,[hep-th/0608022]. [55] A. Rajantie and S. Stopyra, Standard Model vacuum decay in a de Sitter Background, Phys. Rev. D 97 (2018) 025012,[1707.09175].
– 47 –