Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Isospin QCD as a Laboratory for Dense QCD
Next Article in Special Issue
Coupled Quintessence Inspired by Warm Inflation
Previous Article in Journal
Wormhole Restrictions from Quantum Energy Inequalities
Previous Article in Special Issue
Radial Oscillations of Dark Matter Stars Admixed with Dark Energy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stochastic Tunneling in de Sitter Spacetime

1
Department of Physics, Kobe University, Kobe 657-8501, Japan
2
International Center for Quantum-Field Measurement Systems for Studies of the Universe and Particles (QUP), KEK, Tsukuba 305-0801, Japan
3
Particle Theory and Cosmology Group, Center for Theoretical Physics of the Universe, Institute for Basic Science (IBS), Daejeon 34126, Republic of Korea
*
Author to whom correspondence should be addressed.
Universe 2024, 10(7), 292; https://doi.org/10.3390/universe10070292
Submission received: 2 June 2024 / Revised: 8 July 2024 / Accepted: 10 July 2024 / Published: 11 July 2024
(This article belongs to the Special Issue Cosmological Models of the Universe)

Abstract

:
Tunneling processes in de Sitter spacetime are studied by using the stochastic approach. We evaluate the Martin–Siggia–Rose–Janssen–de Dominicis (MSRJD) functional integral by using the saddle-point approximation to obtain the tunneling rate. The applicability conditions of this method are clarified using the Schwinger–Keldysh formalism. In the case of a shallow potential barrier, we reproduce the Hawking–Moss (HM) tunneling rate. Remarkably, in contrast to the HM picture, the configuration derived from the MSRJD functional integral satisfies physically natural boundary conditions. We also discuss the case of a steep potential barrier and find an interesting Coleman–de Luccia (CDL) bubblelike configuration. Since the starting point of our analysis is the Schwinger–Keldysh path integral, which can be formulated in a more generic setup and incorporates quantum effects, our formalism sheds light on further studies of tunneling phenomena from a real-time perspective.

1. Introduction

Tunneling is an important nonperturbative phenomenon in cosmology. Hence, it is crucial to understand tunneling processes in curved spacetime. In particular, tunneling in de Sitter spacetime is worth studying in detail. The reason is that the de Sitter spacetime is the simplest non-trivial curved spacetime. Moreover, the results obtained there can be applicable to tunneling phenomena during inflation. Since both Minkowski and de Sitter spacetime have the maximal symmetry, there exist a similarity. Indeed, the Euclidean instanton method for evaluating the tunneling rate in Minkowski spacetime [1,2] is applicable to tunneling processes in de Sitter spacetime [3]. However, the extension of the Euclidean method to the de Sitter spacetime has a difficulty in interpretation. In the case of the de Sitter spacetime, there are two saddle solutions: Coleman–De Luccia (CDL) [3] and Hawking–Moss (HM) instantons [4]. The latter instanton is a homogeneous solution. Therefore, it is not straightforward to interpret the tunneling rate as that for the tunneling process from a false vacuum to a true vacuum [5,6].
To circumvent the difficulty, a real-time formalism might be useful. In the case of the de Sitter spacetime in the flat chart, the stochastic formalism has been used for investigating time evolution of fluctuations during inflation [7,8]. In this formalism, Langevin equations describe the dynamics of long-wavelength modes. Stochastic noises stem from short-wavelength quantum modes. The stochastic formalism has also been used to analyze the tunneling processes in de Sitter spacetime [9,10,11,12,13,14,15]. However, most of the previous works treated only the HM transition by utilizing the Fokker–Plank (FP) equation for a homogeneous scalar field, which can be derived from the Langevin equation. Obviously, the HM transition cannot describe the bubble nucleation in contrast to CDL tunneling. To describe the bubble, we have to include the spatial gradient term, and the FP equation becomes a functional differential equation instead of a partial differential equation. Hence, at first glance, it seems difficult to treat the CDL bubble by the stochastic approach. However, we make an observation that there exists a path integral representation to solutions of the FP equation. In condensed matter physics, it is called the Martin–Siggia–Rose–Janssen–de Dominicis (MSRJD) functional integral [16,17,18,19].
In this paper, we use the MSRJD path integral formula to study the tunneling processes on the fixed de Sitter background. The applicability conditions of this method are clarified based on the Schwinger–Keldysh formalism. We consider a potential V ( ϕ ) which has a local minimum at ϕ = ϕ false , a true minimum at ϕ true , and a local maximum at ϕ top as shown in Figure 1. We use the saddle-point method and evaluate the tunneling rate by computing the action for the tunneling configurations [20]. We investigate both HM and CDL tunneling processes. In the case of a shallow potential barrier, the conventional HM tunneling rate is reproduced. Remarkably, the tunneling configuration in the MSRJD functional integral method has physically natural boundary conditions in contrast to the Euclidean method. We then clarify that the HM tunneling rate is the tunneling rate of a (roughly Hubble sized) coarse-grained patch transitioning from ϕ false to ϕ top . In the case of a steep potential barrier, we find an interesting CDL bubblelike configuration which is obtained by solving the solution of a scalar field equation in Euclidean anti-de Sitter space even though we do not work in the imaginary time. Our results show how the bubble nucleation process could be described in the stochastic approach.
This paper is organized as follows. In Section 2, we present the setup and review the stochastic approach. In Section 3, the MSRJD functional integral is constructed from the Langevin equations. We also outline the derivation of the MSRJD functional integral starting from the Schwinger–Keldysh formalism. In Section 4, we apply the MSRJD functional integral to the case of the shallow potential barrier and study the HM tunneling process. In Section 5, we consider the case of the steep potential barrier and find the CDL bubblelike configuration. We derive the tunneling rate using the configuration and compare it with the one obtained in Section 4. We also compare our results with those in the Euclidean method. Section 6 is devoted to the conclusion. Some technical details are presented in appendices.

2. Stochastic Approach

In this section, we review the stochastic approach [7,8].

2.1. Setup

We consider a real scalar field ϕ in a (3 + 1)-dimensional de Sitter background. The background geometry is described by the metric
d s 2 = d t 2 + a ( t ) 2 d x 2 = a ( η 2 ) ( d η 2 + d x 2 ) ,
with the scale factor
a ( t ) = e H t = 1 H η ,
where H is the Hubble constant. The action is given by
S = d t d 3 x a 3 1 2 ϕ ˙ 2 1 2 a 2 ( ϕ ) 2 V ( ϕ ) ,
where a dot denotes a derivative with respect to t and V ( ϕ ) is a potential function. From this action, the canonical conjugate momentum Π is defined as Π : = a 3 ϕ ˙ and the Hamilton‘s equations of motion are given by
ϕ ˙ = a 3 Π , Π ˙ = a 2 ϕ a 3 V ( ϕ ) .

2.2. Langevin Equation

First, we divide the quantum fields ϕ and Π into an infrared (IR) part and an ultraviolet (UV) part, respectively, as
ϕ = ϕ IR + ϕ UV , Π = Π IR + Π UV .
The UV part is defined as
ϕ UV : = d 3 k ( 2 π ) 3 θ ( k k c ( t ) ) a k u k ( t ) e i k · x + a k u k ( t ) e i k · x , Π UV : = d 3 k ( 2 π ) 3 θ ( k k c ( t ) ) a ( t ) 3 a k u ˙ k ( t ) e i k · x + a k u ˙ k ( t ) e i k · x ,
where u k is a mode function that is specified later, and a k and a k are annihilation and creation operators, respectively. Here, θ ( k k c ( t ) ) is a step function, and k c ( t ) is the cut-off scale defined as
k c ( t ) : = ε a ( t ) H ,
where ε is a constant. The value of ε is discussed later. Substituting Equation (6) into Hamilton‘s Equation (4), we obtain the following equations;
ϕ ˙ IR = a 3 Π IR + ξ ϕ , Π ˙ IR = a 2 ϕ IR a 3 V ( ϕ IR ) + ξ Π ,
where ξ ϕ and ξ Π are defined as
ξ ϕ : = k ˙ c ( t ) d 3 k ( 2 π ) 3 δ ( k k c ( t ) ) a k u k ( t ) e i k · x + a k u k ( t ) e i k · x , ξ Π : = k ˙ c ( t ) d 3 k ( 2 π ) 3 δ ( k k c ( t ) ) a k a ( t ) 3 u ˙ k ( t ) e i k · x + a k a ( t ) 3 u ˙ k ( t ) e i k · x .
To derive Equations (8) and (9), we assumed that the time evolution of UV modes could be well approximated by the free theory that is defined around the false vacuum. The mode function u k ( t ) introduced in (9) then satisfies the following equation
u ¨ k + 3 H u ˙ k + a 2 k 2 u k + V ( ϕ false ) u k = 0 .
From this equation, we see the mode function only depends on k : = | k | . Hence, we simply denote the mode function as u k . Note that this mode function is used for UV modes k > k c ( t ) .
Let us consider statistical properties of ξ ϕ and ξ Π . The annihilation and creation operators a k and a k satisfy the following commutation relations:
[ a k , a k ] = ( 2 π ) 3 δ ( 3 ) ( k k ) , [ a k , a k ] = [ a k , a k ] = 0 .
Defining the vacuum state | 0 as
a k | 0 = 0 , k ,
we can calculate the one-point and two-point correlation functions of ξ ϕ and ξ Π as
0 | ξ α ( t , x ) | 0 = 0 , ( α , β = ϕ , Π ) , 0 | ξ α ( t , x ) ξ β ( t , x ) | 0 = 1 2 π 2 k ˙ c ( t ) k c ( t ) 2 sin ( k c ( t ) r ) k c ( t ) r g α β ( t ) δ ( t t ) ,
where r : = | x x | , and
g ϕ ϕ ( t ) : = | u k c ( t ) | 2 g Π Π ( t ) : = a ( t ) 6 | u ˙ k c ( t ) | 2 g ϕ Π ( t ) = ( g Π ϕ ) : = a ( t ) 3 u k c ( t ) u ˙ k c ( t ) . .
Because of the Wick theorem, all higher correlation functions can be decomposed into two-point functions. In the stochastic approach, we replace the operator ξ α by the real random field which has the following statistical properties
ξ α ( t , x ) ξ = 0 , ( α , β = ϕ , Π ) , ξ α ( t , x ) ξ β ( t , x ) ξ = 1 2 π 2 k ˙ c ( t ) k c ( t ) 2 sin ( k c ( t ) r ) k c ( t ) r Re [ g α β ( t ) ] δ ( t t ) ,
where ξ represents an expectation value with the distribution function of ξ α . With this prescription, Equation (8) can be interpreted as the Langevin equations with the noise ξ α stemming from quantum fluctuations. Since ξ α are now regarded as classical variables, we ignore the imaginary parts of g ϕ Π and g Π ϕ which come from non-commutativity. Accordingly, IR variables ϕ IR and Π IR are also regarded as classical stochastic variables in this prescription. The validity of this classical approximation is discussed in the next section.

2.3. Bunch–Davies Vacuum

Now we define the vacuum state | 0 by specifying the mode function u k . Equation (10) is the equation for the mode k > k c . Let us assume that we can ignore the term V in solving (10). This can be justified when | V | is sufficiently small to satisfy max ( k c / a ) 2 , H 2 | V | . Under this approximation, (10) can be solved as
u k = ( k η ) 3 2 C 1 H 3 2 ( 1 ) ( k η ) + C 2 H 3 2 ( 2 ) ( k η ) ,
where C 1 , C 2 are constants of integration, and H ν ( 1 ) ( z ) , H ν ( 2 ) ( z ) are Hankel functions. We choose the Bunch–Davies vacuum, i.e., the vacuum that coincides with Minkowski‘s vacuum at past infinity ( η ) so that
u k H η 2 k e i k η , ( k η ) .
This determines C 1 , C 2 as
C 1 = π 4 k H k , C 2 = 0 .
Thus, we obtain the mode function
u k = H η 2 k 1 i k η 1 e i k η .
To check (17), we used the following formulae
H 3 2 ( 1 ) ( z ) = 2 π z 1 i z 1 e i z , H 3 2 ( 2 ) ( z ) = 2 π z 1 i z 1 e i z .
Substituting (19) into (14) and using (7), we obtain
g ϕ ϕ = H 2 η 2 2 k c 1 k c 2 η 2 + 1 = 1 2 H a 3 ε 3 + ε 1 g Π Π = k c 2 H 2 η 2 = H a 3 2 ε g ϕ Π = g Π ϕ = 1 2 1 k c η + i = 1 2 ε 1 + i 2 .
From (8), we have to scale g α β as g ϕ ϕ g ϕ ϕ , g Π Π g Π Π / ( H 2 a 6 ) and g ϕ Π g ϕ Π / ( H a 3 ) in order to compare these quantities. Then, g ϕ ϕ and g Π Π become dominant for the case ε 1 and ε 1 , respectively.

2.4. Reduced Langevin Equation

Let us focus on physics at the super-horizon scales by taking ε 1 and assume a shallow potential so that we can use the massless approximation to evaluate the noise correlations. In this setup, we can use (21), from which we find that g ϕ ϕ becomes dominant. We then ignore the noise for momentum ξ Π . We also ignore the gradient term in (8) because we focus on the physics at super-horizon scales. Since the potential is assumed to be shallow, the time variation in ϕ IR ( t ) is then very small. Therefore, we may integrate the second equation of (8) as
Π IR a 3 3 H V ( ϕ IR ) .
Substituting this relation into the first equation of (8), we obtain
ϕ ˙ IR = V ( ϕ IR ) 3 H + ξ ϕ ,
where the two-point correlation function of ξ ϕ is given at the leading order by
ξ ϕ ( t , x ) ξ ϕ ( t , x ) ξ = H 3 4 π 2 j 0 ( k c r ) δ ( t t ) , j 0 ( k c r ) : = sin ( k c r ) k c r .
Equation (23) is the Langevin equation with the noise correlation (24) derived in [8].

3. Path Integral Formalism for the Stochastic Approach

In this section, we introduce a path integral approach for the Langevin equation following [21]. From now on, we replace the IR field and its conjugate momentum as ( ϕ IR , Π IR ) ( ϕ c , Π c ) . As we see in Section 3.2, this notation corresponds to the Keldysh basis. In Section 3.2, we provide a first-principle derivation of the stochastic approach based on the Schwinger–Keldysh formalism following [22,23]1. We briefly show the sketch of the derivation in the main text, and details of derivations are given in Appendix A.

3.1. MSRJD Functional Integral Representation

3.1.1. Zero-Plus-One-Dimensional Theories

Equation (23) can be regarded as the equation determining the position of a particle ϕ ( t ) with the white noise ξ ϕ defined by Equation (24). For k c r 1 , we have sin ( k c r ) / ( k c r ) 1 as a good approximation. Let us discretize the field and noise as ϕ c ( t ) ϕ i , ξ ϕ ( t ) ξ i ϕ , ( i Z ) . Then, Equation (23) reads
X i : = ϕ i ϕ i 1 + Δ t V ( ϕ i ) 3 H ξ i ϕ = 0 ,
where Δ t is the temporal discretization interval. Note that the discretization (25) is called Ito discretization2. Denoting a solution of (23) by ϕ c [ ξ ] , the expectation value of an observable O ξ can be formally expressed as
O [ ϕ c ( t ) ] ξ = D ϕ c O [ ϕ c ] δ ( ϕ c ϕ c [ ξ ] ) ξ = D ϕ c O [ ϕ c ] | δ X δ ϕ c | δ ( X ) ξ .
where D ϕ c : = i d ϕ i is the functional measure, δ ( ϕ c ϕ c [ ξ ] ) : = i δ ( ϕ i ϕ i [ ξ ] ) , and | δ X / δ ϕ c | is the determinant of { δ X i / δ ϕ j } . When we choose the above discretization (25), δ X / δ ϕ c becomes a triangular matrix with a unit in the diagonal components, i.e., the functional determinant is unity. Thus, substituting (25) into (26), we obtain
O [ ϕ c ( t ) ] ξ = D ϕ c O [ ϕ c ] i δ ϕ i ϕ i 1 + Δ t V ( ϕ i ) 3 H ξ i ϕ ξ .
Representing the delta functions in terms of a Fourier integral and taking the continuum limit, we obtain the following path integral representation
O [ ϕ c ( t ) ] ξ = D ( ϕ c , ϕ ˜ ) O [ ϕ c ] exp i d t ϕ ˜ ϕ ˙ c + V ( ϕ c ) 3 H ξ ϕ ξ .
After averaging over the noise assuming Gaussian statistics, we obtain the Martin–Siggia–Rose–Janssen–de Dominicis (MSRJD) functional integral
O [ ϕ c ( t ) ] ξ = D ( ϕ c , ϕ ˜ ) O [ ϕ c ] exp d t i ϕ ˜ ϕ ˙ c + V ( ϕ c ) 3 H H 3 8 π 2 ϕ ˜ 2 .
Putting O [ ϕ c ] = δ ( ϕ c ϕ c ( t ) ) and taking the boundary conditions ϕ c ( t ) = ϕ c and ϕ c ( t ) = ϕ c , we obtain the transition probability p ( ϕ c , t | ϕ c , t )
p ( ϕ c , t | ϕ c , t ) = ϕ c ( t ) = ϕ c ϕ c ( t ) = ϕ c D ( ϕ c , ϕ ˜ ) exp d t i ϕ ˜ ϕ ˙ c + V ( ϕ c ) 3 H H 3 8 π 2 ϕ ˜ 2 .
Finally, introducing new variables Π Δ : = i ϕ ˜ , we obtain the “phase space” path integral
p ( ϕ c , t | ϕ c , t ) = ϕ c ( t ) = ϕ c ϕ c ( t ) = ϕ c D ( ϕ c , Π Δ ) exp d t ( Π Δ ϕ ˙ c H ( ϕ c , Π Δ ) ) ,
where we defined the Hamiltonian as
H ( ϕ c , Π Δ ) : = V ( ϕ c ) 3 H Π Δ H 3 8 π 2 Π Δ 2 .
We emphasize that this is just a change of integration variable, not the Wick rotation. A set of saddle-point equations is defined as the Hamilton equations for the Hamiltonian defined above.

3.1.2. Three-Plus-One-Dimensional Theories

Let us consider the case of the full Langevin Equations (8), (13) and (21). When the space is discretized, (8) can be seen as the dynamics of many body particles. Thus, the procedure in the previous subsection can be utilized. Taking the continuum limit, we obtain the following expression which is analogous to (30),
p ( ϕ c ( x ) , t | ϕ c ( x ) , t ) = ϕ c ( t , x ) = ϕ c ( x ) ϕ c ( t , x ) = ϕ c ( x ) D ( ϕ c , Π c , ϕ Δ , Π Δ ) × exp i d 4 x Π Δ ϕ ˙ c ϕ Δ Π ˙ c H ( ϕ c , Π c , ϕ Δ , Π Δ ) .
We define the Hamiltonian H ( ϕ c , Π c , ϕ Δ , Π Δ ) as
H ( ϕ c , Π c , ϕ Δ , Π Δ ) = Π c Π Δ a 3 ( a 2 ϕ c a 3 V ( ϕ c ) ) ϕ Δ i 2 α , β d 4 x X α ( x ) G α β ( x , x ) X β ( x ) ,
where the Greek indices α , β = ( ϕ , Π ) label the Δ fields as ( X ϕ , X Π ) = ( Π Δ , ϕ Δ ) , and G α β denotes the correlations of noises:
G α β ( x , x ) : = ξ α ( x ) ξ β ( x ) ξ ,
where the RHS is given in (15). We use the formula (33) in Section 4.2 and Section 5. Note that we have not changed the integration variables as ( ϕ Δ , Π Δ ) ( i ϕ Δ , i Π Δ ) .
In Section 4 and Section 5, we take the initial state and the final state to be the false vacuum and the field configurations after tunneling (true vacuum or bubble), respectively.

3.2. Transition Probability from the Schwinger–Keldysh Formalism

In the above prescription, we treat IR fields ( ϕ IR , Π IR ) as classical stochastic variables. To quantify the validity of this approximation, we derive the stochastic approach from the first principles based on the Schwinger–Keldysh formalism. Note that we use the label “IR” for IR fields in this subsection following the notation adopted in Appendix A.
Suppose that the system is in the false vacuum | Ψ false at the past infinity t = , and it evolves into a certain quantum state | { ϕ IR ( T , x ) } x D at the final time t = T , where | { ϕ IR ( T , x ) } x D is an eigenstate of an IR field ϕ ^ IR ( x ) at spatial points x in the domain D with an eigenvalue ϕ IR ( T , x ) . Here, ϕ ^ IR ( x ) is an IR field in the Schrödinger picture3. The transition probability p for this process is given by
p ( { ϕ IR ( T , x ) } x D ) = { ϕ IR ( T , x ) } x D | U ^ ( t , ) | Ψ false 2 ,
where U ^ ( t , t ) describes the unitary time evolution from t to t. The RHS is the Fourier component of the generating functional Z [ J IR ] for IR fields, and we have
p ( { ϕ IR ( T , x ) } x D ) = x D d J IR ( T , x ) Z [ J IR ( T ) ] e i J IR ( T , x ) ϕ IR ( T , x ) { J IR ( T , y ) = 0 } y D ,
where the generating functional Z [ J IR ( T ) ] is defined by
Z [ J IR ( T ) ] : = Tr U ^ ( T , ) | Ψ false Ψ false | U ^ ( T , ) x e i J IR ( T , x ) ϕ IR ( T , x ) .
We can evaluate Z [ J IR ( T ) ] nonperturbatively in the IR sector based on the method developed in [22,23]. Our strategy is that we integrate out UV sector k > k c ( t ) perturbatively in nonlinear couplings to obtain Z [ J IR ( T ) ] . This is analogous to tracing out environmental degrees of freedom to evaluate the reduced density matrix for the system under considerations. Each step we take can be summarized as follows:
  • First, we derive a path integral representation of Z [ J IR ( T ) ] . We then split the integration variables into UV fields and IR fields such that the integration contour of UV variables of ϕ is closed; see discussions around Equation (A4) for more details4.
  • We perform the integration over UV variables and evaluate an IR effective action, the so-called Feynman–Vernon influence functional [25] perturbatively.
  • We substitute the obtained expression of Z [ J IR ( T ) ] into (37), giving rise to the path integral expression for the transition probability p.
Technical details are shown in Appendix A. At step 2, we assume that the quantum state for nonzero modes is given by the Bunch–Davies vacuum state for a free field defined around the false vacuum, and that the exact zero mode provides a non-fluctuating classical background field configuration ϕ = ϕ false .
In this way, we can obtain the path integral expression for the probability p as
p ( { ϕ IR ( T , x ) } x D ) = D ( ϕ c IR , ϕ Δ IR , Π c IR , Π Δ IR ) e i S δ C fin δ C ini .
Here, δ C fin and δ C ini fix the boundary conditions of the path integral at the final time and the initial time, respectively:
δ C fin : = x D δ ( ϕ c IR ( T , x ) ϕ IR ( T , x ) ) y δ ( ϕ Δ IR ( T , y ) ) ,
δ C ini : = x δ ϕ c IR ( , x ) ϕ false δ Π c IR ( , x ) δ [ C detail ] ,
where δ [ C detail ] is not relevant for the main points here, and we define it in Appendix A. A term S may be understood as an IR effective action,
S : = d 4 x Π Δ IR ϕ ˙ c IR ϕ Δ IR Π ˙ c IR H ( ϕ c IR , Π c IR , ϕ Δ IR , Π Δ IR ) δ H ,
δ H : = a 3 ( t ) V ( ϕ c IR + ( ϕ Δ IR / 2 ) ) V ( ϕ c IR ( ϕ Δ IR / 2 ) ) V ( ϕ c IR ) ϕ Δ IR + δ H higher ,
where H is given in (34). δ H higher denotes the higher-order terms in the coupling constant. Note that Equation (39) allows us to understand the transition probability p in terms of the stochastic dynamics: this point is discussed in Appendix A.3.

Validity of Equation (33)

Now, we find that S coincides with the term in the exponent of (33) up to the term δ H . Hence, (39) shows that the result (33) in the previous section can be justified from the first principles when the term δ H is negligible. Suppose that the potential V ( ϕ ) is sufficiently flat in the regime of interest such that the perturbation theory works when integrating out UV modes. In this case, we can systematically calculate δ H higher perturbatively, and we have δ H higher = 0 at the leading order. We set δ H higher = 0 in the present analysis. We should keep in mind however that the nonperturbative physics at UV scales k > k c ( t ) = ε a ( t ) H is lost by setting δ H higher = 0 . From this viewpoint, we should choose ε as large as possible.
Even after setting δ H higher = 0 , we have nonzero δ H . In the classical approximation, we may regard the Δ -variables ( ϕ Δ IR , Π Δ IR ) as tiny quantities and keep only the terms which are linear in the Δ -variables in S , leading to δ H 0 . In the case ε 1 , this approximation would work thanks to the squeezing of quantum fluctuations at super-horizon scales. However, it is more subtle if we can set δ H = 0 in the case ε 1 . We revisit this issue in Section 5.

4. Hawking–Moss Tunneling

In this section, we discuss the Hawking–Moss (HM) tunneling from the perspective of the MSRJD functional integral. We calculate the tunneling rate defined by (31) and (33) in Section 4.1 and Section 4.2, respectively.
Hereafter, for the numerical calculations, we consider the potential [26] in Figure 2;
V ( ϕ ) = g 2 4 ( ϕ 2 α 2 ) 2 + 4 β 3 ( α ϕ 3 3 α 3 ϕ 2 α 4 ) ,
for which the derivative with respect to ϕ is given by a simple form
V ( ϕ ) = g 2 ( ϕ α ) ( ϕ + α ) ( ϕ + α β ) .
For this potential, we have ϕ false = α , ϕ true = α and ϕ top = β α . The dimensionless parameter β must be 0 < β < 1 , otherwise ϕ = α is not a false vacuum. The potential at each point takes the following values:
V ( α ) = 0 , V ( α ) = 4 β g 2 α 4 3 , V ( β α ) = g 2 α 4 12 ( 1 β ) 3 ( 3 + β ) .
The height of potential barrier Δ V between the false vacuum and the true vacuum is
Δ V : = V ( β α ) V ( α ) = g 2 α 4 12 ( 1 β ) 3 ( 3 + β ) .
The mass of fluctuations around a given point ϕ is
V ( ϕ ) = α 2 g 2 3 ( ϕ / α ) 2 + 2 β ( ϕ / α ) 1 .
In our analysis, we work in the quantum field theory on the fixed de Sitter background. This assumption is valid when the condition max | V ( α ) | , | V ( β α ) | , | V ( α ) | 3 M pl 2 H 2 is satisfied. This is always satisfied when the following condition is imposed
H 3 g M pl 2 H g α 2 2.9 × 10 12 GeV g 1.12 × 10 4 140 g α / H 2 ,
where we use 0 < β < 1 and choose parameters ( g , g α / H ) = ( 1.12 × 10 4 , 140 ) as a benchmark point for later convenience. This shows that we can ignore the backreaction consistently with the weak coupling g 1 , provided that H is well below the Planck scale M pl 2.44 × 10 18 GeV. Note that under this condition, our discussion below is also applicable in the case of a quasi de Sitter background. However, it is non-trivial to generalize our analysis to the case when a non-negligible backreaction on H exists. This is because it is necessary to extend the first-principle derivation of the stochastic approach to account for the backreaction. This is in itself an interesting issue which we leave for future work.

4.1. The Case for the Reduced Langevin Equation

Let us start from the simpler formula (31). Hamilton‘s equations for the Hamiltonian H ( ϕ c , Π Δ ) defined in (32) are given by
ϕ ˙ c = V ( ϕ c ) 3 H H 3 4 π 2 Π Δ Π ˙ Δ = V ( ϕ c ) 3 H Π Δ .
There exist many Hamiltonian flow lines in phase space ( ϕ c , Π Δ ) that represent the solution of Hamilton‘s equations as shown in Figure 3. Each flow line can be specified by the value of the Hamiltonian since it is conserved on the given flow line. In the Figure 3, there are two important flow lines corresponding to H ( ϕ c , Π Δ ) = 0 . Since the Hamiltonian is the quadratic polynomial of Π Δ , there are two Π Δ satisfying H ( ϕ c , Π Δ ) = 0 :
Π Δ = 0 , 8 π 2 3 H 4 V .
In the case of Π Δ = 0 , Equation (48) reduces to
ϕ ˙ c = V ( ϕ c ) 3 H .
Similarly, in the case of Π Δ = 8 π 2 3 H 4 V ( ϕ c ) , Equation (48) reduces to
ϕ ˙ c = V ( ϕ c ) 3 H .
These two flow lines have three intersections which correspond to the false vacuum, the top of the potential hill and the true vacuum from left to right. On the flow line with Π Δ = 0 , the false and true vacuums are stable and the top of the potential is unstable, while on the flow line with Π Δ 0 , the stability is reversed. The different signs of the potential forces account for the reverse of the stability at the intersections. We see that the tunneling configuration exists; starting from the false vacuum ( ϕ c , Π Δ ) = ( α , 0 ) , going to the top of the potential ( ϕ c , Π Δ ) = ( β α , 0 ) through the flow line with Π Δ 0 and finally going to the true vacuum through the flow line with Π Δ = 0 .
We numerically solve the second equation with the initial conditions ϕ c ( t ) = α and switch to the first equation when the field reaches the top of potential ϕ c ( t ) = β α . Thus, we obtain the tunneling configuration in Figure 4.
Now, we evaluate the transition probability (31) by substituting the above tunneling configurations. The action in (31) can be calculated without concrete solutions for ϕ c ;
I = t t d t Π Δ ϕ ˙ c H ( ϕ c , Π Δ ) = 8 π 2 3 H 4 t t d t ϕ ˙ c V ( ϕ c ) = 8 π 2 3 H 4 Δ V ,
where t is the time at which ϕ c reaches the top of the potential, and Δ V is the difference of energy density between the top of the potential and the false vacuum. The trajectory from the top of the potential to the true vacuum does not contribute to the transition probability because H ( ϕ c , Π Δ ) = 0 and Π Δ = 0 on its trajectory. Thus, the transition probability of the tunneling can be evaluated as
p ( α , t | α , t ) exp 8 π 2 3 H 4 Δ V .
This shows a complete agreement with the result of the HM instanton on the fixed background [27].
Comments on this result are in order. It is pointed out in [27] that when the potential V ( ϕ ) has several degenerate local maxima, all of them yield the same tunneling action in the Euclidean method, despite differing distances from the false vacuum. In our formalism, however, a transition to a distant maximum is described by a saddle-point solution that passes through all the intermediate maxima, and the transition probability is given by the product of the factors on the RHS of (53) for each maximum. The aforementioned subtlety is absent in our formalism.
Since the analysis in Section 4.1 is based on (23) describing the stochastic dynamics of the IR field ϕ c at a single spatial point, Equation (53) calculates the tunneling probability at a single spatial point, not on the whole universe5. Hence, physically speaking, (53) provides the tunneling probability of a coarse-grained patch with a physical radius ( ε H ) 1 . This point also becomes manifest in Section 4.2.

4.2. The Case for the Full Langevin Equation

Next, we evaluate the tunneling rate based on the formula (33). An advantage of (33) is that one can discuss various non-trivial dynamics in the whole universe covered by the flat chart in principle.
In the previous section, we focused on the dynamics at a single spatial point and found the HM tunneling process as a non-trivial solution of Hamilton‘s equations. In this section, we discuss the corresponding process in the global picture based on (33) by identifying a non-trivial configuration that satisfies Hamilton‘s equations in a good approximation. We then reproduce the result (52). We also estimate the characteristic time scale of the tunneling. In Section 4.2.1, we remark some technical complications that arise in the global picture.
We start with simplifying the expression (33). As discussed in Section 2, G ϕ ϕ becomes dominant, and the other terms can be neglected for ε 1 in (34). Under this approximation, the exponent of the integrand of (33) is linear in ϕ Δ . We can then perform the path integral D ϕ Δ in (33), yielding the product of delta functions x δ Π ˙ c a 2 ϕ c + a 3 V ( ϕ c ) . These delta functions give an equation of motion for Π c and eliminate the path integral D Π c . We solve this equation as Π c a 3 V ( ϕ c ) / 3 H , ignoring the spatial gradient and the time variation in ϕ c . Furthermore, we change the integration variable as Π Δ i Π Δ in (33). Consequently, (33) reduces to
p ( ϕ c ( x ) , t | ϕ c ( x ) , t ) ϕ c ( t , x ) = ϕ c ( x ) ϕ c ( t , x ) = ϕ c ( x ) D ( ϕ c , Π Δ ) exp d 4 x Π Δ ϕ ˙ c H HM ( ϕ c , Π Δ ) ,
where the Hamiltonian H HM is defined as
H HM ( ϕ c , Π Δ ) : = 1 3 H V ( ϕ c ) Π Δ H 3 8 π 2 Π Δ Π ¯ Δ ,
and Π ¯ Δ is given by
Π ¯ Δ ( t , x ) : = d 3 x j 0 k c ( t ) | x x | Π Δ ( t , x ) = d 3 k ( 2 π ) 3 4 π 2 H 3 H 2 2 k 3 δ ( t t k ) Π Δ ( t , k ) e i k · x .
Here, we define t k by the condition k c ( t k ) = k . Hamilton‘s equations are
ϕ ˙ c = 1 3 H V ( ϕ c ) H 3 4 π 2 Π ¯ Δ ,
Π ˙ Δ = 1 3 H V ( ϕ c ) Π Δ .
Here, the term Π ˙ Δ ( t , x ) on the LHS of (57b) is defined as the Fourier transform of Π ˙ Δ ( t , k ) i.e., Π ˙ Δ ( t , x ) : = ( 2 π ) 3 d 3 k Π ˙ Δ ( t , k ) e i k · x . The same rule applies to other IR variables in the real space as in Equation (A24). Intuitively, this is because the stochastic formalism is formulated in the momentum space, and its Fourier transform provides the stochastic formalism in the real space; see Appendix A for more details.
Equation (57) admits a trivial configuration Π Δ ( t , k ) = 0 under which we have H HM = 0 , and (57a) becomes ϕ ˙ c = V ( ϕ c ) / ( 3 H ) , describing the standard classical time evolution. Now, we are interested in the non-trivial configuration which satisfies (57). In the previous section, we focused on the non-trivial dynamics at the single spatial point x = x 0 and it was found that the tunneling configuration was obtained by Π Δ ( t , x 0 ) = 8 π 2 3 H 4 V ( ϕ c ( t , x 0 ) ) . This configuration can be naturally extended to the one in the global picture by considering the following configuration
Π Δ ( t , k ) = 8 π 2 3 H 4 V ( ϕ c ( t , x 0 ) ) e i k · x 0 .
From now on, we show that this configuration satisfies a set of Hamilton Equations (57) in a good approximation. For this configuration, we have
Π Δ ( t , x ) = 8 π 2 3 H 4 V ( ϕ c ( t , x 0 ) ) f ( r ; t ) ,
Π ˙ Δ ( t , x ) = 8 π 2 3 H 4 V ( ϕ c ( t , x 0 ) ) ϕ ˙ c ( t , x 0 ) f ( r ; t ) ,
Π ¯ Δ ( t , x ) = 8 π 2 3 H 4 V ( ϕ c ( t , x 0 ) ) j 0 k c ( t ) r ,
where r : = | x x 0 | , and f ( r ; t ) : = ( 2 π ) 3 d 3 k θ ( k c ( t ) k ) e i k · ( x x 0 ) . The functions j 0 ( k c ( t ) r ) and f ( r ; t ) rapidly oscillate for r k c 1 ( t ) . Hence, they exponentially decay to zero after they are smeared over the Hubble time: for instance, j 0 and f becomes the smooth function after the suitable coarse-graining as6
j 0 ( k c ( t ) r ) time c . g . W ( r ; t ) : = d t w ( t , t ) j 0 k c ( t ) r 1 ( r k c 1 ( t ) ) 0 ( r k c 1 ( t ) ) , f ( r ; t ) time c . g . W f ( r ; t ) : = d t w ( t , t ) f ( r ; t ) k c 3 ( t ) / ( 6 π 2 ) ( r k c 1 ( t ) ) 0 ( r k c 1 ( t ) ) .
Here, we perform the smearing by using a window function w ( t , t ) which is approximately constant in the domain | t t | H 1 while it exponentially decays to zero at | t t | H 1 . We also impose the normalization condition d t w ( t , t ) = 1 . We can choose w ( t , t ) = H π 1 / 2 exp [ H 2 ( t t ) 2 ] for example. Note that the property (60) is based on the behavior j 0 ( k c r ) 1 and f ( r ; t ) k c 3 / ( 6 π 2 ) for r k c 1 and the rapid oscillations of j 0 ( k c r ) and f 0 ( r ; t ) for r k c 1 . Hence, the property (60) is robust against the choice of w ( t , t ) .
Therefore, substituting Equations (59c) and (60) into Equation (57a) and performing the smearing in time, we obtain
ϕ ˙ c ( t , x ) = 1 3 H V ( ϕ c ( t , x ) ) 2 V ( ϕ c ( t , x 0 ) ) j 0 k c ( t ) r
time c . g . V ( ϕ c ( t , x ) ) / ( 3 H ) ( r k c 1 ( t ) ) V ( ϕ c ( t , x 0 ) ) / ( 3 H ) ( r k c 1 ( t ) ) .
This exhibits the non-trivial dynamics only in the domain r k c 1 , i.e., a coarse-grained patch centered at a point x 0 with a physical radius ( ε H ) 1 . Here, we used ϕ c ( t , x ) ϕ c ( t , x 0 ) for r k c 1 ( t ) . In the second line, it is assumed that V ( ϕ c ( t , x 0 ) ) does not cancel the rapid oscillations of j 0 ( k c ( t ) r ) so that their product V j 0 decays exponentially at r k c 1 after the suitable smearing. This assumption is valid thanks to the approximate constancy of super-horizon fluctuations ϕ c ( t , k ) over the Hubble time: ϕ c ( t , k ) ϕ c ( t , k ) for | t t | H 1 . Adopting the same assumption, from (59a) and (59b), we also obtain the smeared expressions of Π Δ and Π ˙ Δ which turn out to satisfy (57b) under the condition (57a). Hence, we conclude that (59) describes a non-trivial saddle-point solution in a good approximation. Furthermore, after the smearing, we can find H HM 0 for that configuration.
Now, we use the configuration (58) or (59) to obtain the tunneling configuration shown in Figure 4 again. We can calculate the action for this tunneling configuration similarly to (52). By substituting the exact expressions (59) and (61a), we obtain7
I = H 3 8 π 2 t t d t d 3 x Π Δ ( t , x ) Π ¯ Δ ( t , x ) = 8 π 2 3 H 4 t t d t V ( ϕ c ( t , x 0 ) ) V ( ϕ c ( t , x 0 ) ) 3 H d 3 x j 0 k c ( t ) r f ( r ; t ) = 8 π 2 3 H 4 t t d t ϕ ˙ c ( t , x 0 ) V ( ϕ c ( t , x 0 ) ) = 8 π 2 3 H 4 Δ V .
In the first line, we used (57a) to eliminate ϕ ˙ c from the integrand. In the second line, we substituted the concrete configurations (59). In the third line, we used (61a) and performed the spatial integral as d 3 x j 0 k c ( t ) r f ( r ; t ) = 1 . Thus, we reproduced the previous result (52). Our analysis confirms that (62) gives the tunneling probability of a coarse-grained patch with a physical radius ( ε H ) 1 .
It is worth mentioning the time scale of the HM tunneling. We may naively calculate the typical time scale of the HM tunneling as t t = ϕ false ϕ top d ϕ c 3 H V ( ϕ c ) . Here, the spacetime argument of the field is suppressed. However, this integral is in general divergent since V = 0 at both ends of the integral. In the vicinity of the domain where V = 0 , we expect that the quantum fluctuations of the field would play an important role. Since the typical size of quantum fluctuations accumulated over the Hubble time is H / ( 2 π ) , we may define the regulated quantity t HM as
t HM : = ϕ false + H 2 π ϕ top H 2 π d ϕ c 3 H V ( ϕ c ) ,
and we expect that t HM would correctly characterize the typical time scale of the HM tunneling. In particular, for our potential (42), we have
t HM = 3 H 2 g 2 α 2 1 + β 3 + β 1 β log 1 1 β H / 2 π α log 1 2 H / 2 π α 1 + 1 + β H / 2 π α 1 H · O ( H / g α ) 2 .
Since we have ( H / g α ) 2 ( H 2 / V ) ϕ false ϕ c ϕ top 1 for the shallow potential, we conclude that the time scale of HM tunneling is much longer than the Hubble time.

4.2.1. Remark

So far it has been found that starting from (33), the configuration (58) approximately solves Hamilton‘s Equation (57) and reproduces the previous result (52) when it is substituted into the action. There is actually an important reason behind why we should choose the configuration (58) to evaluate the HM tunneling process.
To see this, it is important to realize an important difference between the current Hamilton Equation (57) and the previous one (48); the equation for ϕ c now contains Π ¯ Δ rather than Π Δ . Due to this difference, the structures of Hamilton‘s Equation (57) are understood as follows;
  • The equation for ϕ c (57a) is determined once Π ¯ Δ is specified;
  • Π ¯ Δ ( t , x ) is understood as initial conditions for (57b) since Π ¯ Δ ( t , x ) contains only the boundary modes satisfying k = k c ( t ) . Hence, Π Δ ( t , x ) is obtained by solving (57b) for a given Π ¯ Δ under the condition (57a).
On top of that, substituting Hamilton‘s equation for ϕ c into the action, we can rewrite the first line of (62) as
I = 1 2 t t d t d 3 k ( 2 π ) 3 Π Δ ( t , k ) Π Δ ( t , k ) H 2 2 k 3 δ ( t t k ) .
The RHS depends only on the value of Π Δ ( t , k ) | k = k c ( t ) . Hence, we can evaluate the tunneling action once Π ¯ Δ is specified, provided that (65) or equivalently (33) is reliable.
We claim that (33) is valid when Π ¯ Δ is chosen so that Π Δ ( t > t k , k ) becomes a smooth function in time. Otherwise, it is not justified to ignore the higher-order corrections δ H higher which are neglected in (33). Intuitively, this states that the perturbation theory tends to be broken down around exotic configurations. When Π Δ is a smooth function, we can relate Π ¯ Δ to Π Δ after the suitable smearing in time:
Π ¯ Δ ( t , x ) time c . g . d 3 x Π Δ ( t , x ) W ( | x x | ; t ) V Π Δ ( t , x ) .
Here, we also used Π Δ ( t , x ) Π Δ ( t , x ) for | x x | k c 1 and defined V : = ( 4 π / 3 ) k c 3 ( t ) d 3 x W ( | x x | ; t ) . We then require that Hamilton‘s equation for Π Δ (57b) is satisfied in a good approximation to ensure the validity of (33). This constrains the choice of Π ¯ Δ . Indeed, the configuration (58) approximately solves (57b) and satisfies the smoothness and henceforth (66) after the smearing. This would be the reason why our analysis with (58) based on (33) can reproduce the correct result.

4.3. On the ε -Independence of the Hawking–Moss Tunneling

In Section 4.1 and Section 4.2, the tunneling probability of a coarse-grained patch with a physical radius ( ε H ) 1 is calculated. Our results coincide with the HM tunneling and are independent of ε .
The ε -independence of the results would be a consequence of the scale independence of the dynamics of light scalar fields at super-horizon scales. To see this, it is useful to notice that a coarse-grained patch of physical radius ( ε 1 H ) 1 at a time t = T expands to a patch of larger physical radius ( ε 2 H ) 1 at a later time t = T + δ t > T with ε 2 = ε 1 exp [ H δ t ] < ε 1 1 . The value of the IR field does not evolve from t = T to T + δ t because the fluctuations of a light scalar field at super-horizon scales are approximately time-independent,
d 3 k ( 2 π ) 3 ϕ c ( t , k ) e i k · x θ ( k c ( t ) k ) d 3 k ( 2 π ) 3 ϕ c ( t + δ t , k ) e i k · x θ ( k c ( t ) k ) ,
unless δ t is much longer than H 1 . Hence, we can relate the coarse-grained dynamics of different ε 1 by considering the time shift t t + δ t with the value of the IR field being kept fixed. This means that the tunneling probability from ϕ = ϕ false to ϕ true of a patch with physical radius ( ε H ) 1 should be independent of ε .

5. Coleman–de Luccia Tunneling

In the previous section, the result of the HM instanton was reproduced by using (33). We expect that the Formula (33) can describe not only the HM tunneling but also the CDL tunneling. In other words, we expect a flow line starting from the false vacuum to the bubble configurations as illustrated in Figure 5. In this section, we concretely show an interesting configuration by following the prescription in the previous section.

5.1. Coleman–de Luccia Bubble Solution

In the Euclidean method, the CDL instanton is dominant rather than the HM instanton when the potential barrier is steep [28,29,30,31]. By taking ε 1 , we can study nonperturbative physics at sub-horizon scales such as the formation of a bubble. In this case, G Π Π becomes dominant and the other noises can be ignored in (34). Under this approximation, the exponent of the integrand of (33) is linear in Π Δ . We can then perform the path integral D Π Δ in (33), yielding the product of delta functions x δ ϕ ˙ c a 3 Π c . These delta functions eliminate the path integral D Π c 8. Furthermore, we change the integration variable as ϕ Δ i ϕ Δ in (33), leading to
p D ( ϕ c , ϕ Δ ) exp d 4 x a 3 ϕ ˙ Δ ϕ ˙ c H CDL ( ϕ c , ϕ Δ )
with appropriate boundary conditions; we consider bubble configurations as boundary conditions later. Here, the Hamiltonian H CDL is defined by
H CDL ( ϕ c , ϕ Δ ) : = a 2 ϕ c a 3 V ( ϕ c ) ϕ Δ a 3 2 ϕ Δ ( x ) ϕ ¯ Δ ( x ) ,
with
ϕ ¯ Δ ( t , x ) : = H k c 4 ( t ) 4 π 2 a ( t ) d 3 x j 0 k c ( t ) | x x | ϕ Δ ( t , x ) .
We would like to find the appropriate configuration ϕ ¯ Δ which describes the tunneling process. Following the choice made in Section 4, let us suppose that the appropriate choice is given by the configuration for which the Hamiltonian (69) vanishes:
ϕ ¯ Δ ( x ) = 0 , 2 ( a 2 2 ϕ c V ( ϕ c ) ) .
Note that this condition is imposed up to the suitable smearing in time because ϕ ¯ Δ ( t , x ) includes only the boundary Fourier mode k = k c ( t ) . As discussed in Section 4.2.1, ϕ Δ becomes sufficiently smooth for the appropriate ϕ ¯ Δ , and hence we also suppose
ϕ ¯ Δ H 2 ε 3 π ϕ Δ ( x )
after the suitable averaging over the Hubble time, similarly to (66). We use the relation (72) to evaluate the action later. This would work at least for the purpose of estimating the action.
In principle, we do not need to use the estimate (72) in evaluating the action since it is determined once ϕ ¯ Δ is specified based on the analogous logic discussed around (65). For this purpose, we need to specify the configurations { ϕ Δ ( t , k ) | k = k c ( t ) } in momentum space which result in the sufficiently smooth ϕ Δ that solves the Hamilton equations. It is not easy to find such appropriate configurations precisely, however. Hence, we leave a more careful analysis on the choice of ϕ ¯ Δ for future work. In this study, instead, we focus on how we can proceed with the analysis for the given configuration (71) while adopting the estimate (72) and how the bubble nucleation process could be described in the stochastic approach.
Using the relation (72), the Hamiltonian (69) and the configuration (71) become
H CDL a 2 ϕ c a 3 V ( ϕ c ) ϕ Δ H 2 ε 6 π a 3 ϕ Δ 2 ,
ϕ Δ ( x ) = 0 , 6 π H 2 ε ( a 2 2 ϕ c V ( ϕ c ) ) .
Hamilton’s equations under the constraints (74) give the following equation of motions for ϕ c :
ϕ ¨ c + 3 H ϕ ˙ c = ( V ( ϕ c ) a 2 2 ϕ c ) ,   ( ϕ Δ = 0 ) , ϕ ¨ c + 3 H ϕ ˙ c = V ( ϕ c ) a 2 2 ϕ c , ( ϕ Δ = 6 π H 2 ε ( a 2 2 ϕ c V ( ϕ c ) ) ) .
Similar to the previous case (50), the signs of the potential term and the gradient term are flipped in the equation of motion with the non-trivial ϕ Δ 0 . Hence, we can expect the tunneling process to be realized. From (68), we find that only the solutions satisfying ϕ Δ 0 contribute to the action. Then, we focus on the second equation in (75).
Interestingly, this equation is the classical equation of motion in the Euclidean anti-de Sitter (AdS) space which is defined by the embedding equation
X 0 2 + X 1 2 + X 2 2 + X 3 2 + X 4 2 = H 2
in five-dimensional Minkowski spacetime
d s 2 = d X 0 2 + d X 1 2 + d X 2 2 + d X 3 2 + d X 4 2 .
In fact, the d’Alembert operator t 2 + 3 H t + a 2 2 can be obtained from the following induced metrics;
d s 2 = d t 2 + e 2 H t d x 2 = 1 H 2 η 2 ( d η 2 + d x 2 ) = H 2 ( d ρ 2 + sinh 2 ρ d Ω 2 ) ,
where d Ω 2 : = d θ 1 2 + sin 2 θ 1 d θ 2 2 + sin 2 θ 1 sin 2 θ 2 d θ 3 2 and 0 < ρ < , 0 θ 1 π , 0 θ 2 π and 0 θ 3 2 π . The second metric is that in the Poincaré coordinates and the third one is that in the global coordinates. The explicit coordinate transformations are given by
X 0 = 1 2 η 2 + x 2 + H 2 η = H 1 cosh ρ , X 1 = x 1 H η = H 1 sinh ρ sin θ 1 cos θ 2 , X 2 = x 2 H η = H 1 sinh ρ sin θ 1 sin θ 2 cos θ 3 , X 3 = x 3 H η = H 1 sinh ρ sin θ 1 sin θ 2 sin θ 3 , X 4 = 1 2 η 2 + x 2 H 2 η = H 1 sinh ρ cos θ 1 .
For later convenience, we derive the relation between ( η , r : = | x | ) and ( ρ , σ 1 ) as
H η = 1 cosh ρ sinh ρ cos θ 1 , r = H 1 sinh ρ sin θ 1 cosh ρ sinh ρ cos θ 1 .
From this, we can check that a point ( η , r ) = ( H 1 , 0 ) is mapped to ρ = 0 for any θ 1 [ 0 , π ] . Except for this, there is one-to-one correspondence between points in the ( η , r ) -plane with η 0 and r 0 and points in the ( ρ , θ 1 ) -plane with ρ > 0 and 0 θ 1 π . On a constant- η slice in the ( η , r ) -plane, the value of ρ for a given spatial point r is given by
ρ = arccosh 1 2 H 2 ( η 2 + r 2 ) + 1 H η .
This function increases monotonically in r for a given η in the region r > 0 .
From the global coordinate (78), we see the Euclidean AdS spacetime has O ( 4 ) symmetry. Hence, we assume that the field ϕ c depends only on ρ . Thus, the second equation of (75) reads
d 2 ϕ c d ρ 2 + 3 tanh ρ d ϕ c d ρ = V ( ϕ c ) H 2 .
Note that this is the same equation utilized in the Euclidean method [3,5]. Imposing the boundary conditions [1,3]
lim ρ ϕ c = α , d ϕ c d ρ | ρ = 0 = 0 ,
we obtain the bubble solutions as shown in Figure 6. Note that the latter condition ensures the continuity of η ϕ c ( η , r ) and η 2 ϕ c ( η , r ) at ( H η , H r ) = ( 1 , 0 ) .
It is useful to fit the bubble configurations by the following fitting function (Figure 6);
ϕ c ( ρ ) = α tanh μ ( ρ ρ ¯ ) ,
where μ and ρ ¯ represent the thickness and the position of the bubble wall, respectively. Note that for β 0.6 , the deviation from the true vacuum at ρ = 0 becomes significant [26], and the fitting (84) becomes bad.
We also have to examine the boundary conditions in the Hamilton flow. As in the case of HM tunneling, we expect several intersections where ϕ Δ = 0 (Figure 5). Since we have the concrete ϕ Δ (74) and the bubble configurations (84), it is possible to obtain the curves in the ( r , η ) -plane on which ϕ Δ = 0 . Using (74) and (75), such curves are given by
F ( η , r ) : = η 2 ϕ c 2 η η ϕ c = 0 .
Note that this is equivalent to solving a 2 2 ϕ c V ( ϕ c ) = 0 . Substituting (74) into (85) and using the fitting formula (84), we find four solutions which are expressed as four lines in the ( η , r ) -plane; two of them are placed at η = with fixed r and r = with fixed η . These correspond to ρ = where ϕ = α (Figure 6). Also, we numerically find other two non-trivial hypersurfaces (Figure 7a). It can be seen that one of the curves is totally spacelike but another is partially timelike. Thus, it seems natural to take the η = curve as an initial time slice on which ϕ is a false vacuum and a totally spacelike non-trivial curve as a final time slice on which the bubble is nucleated. We denote the region between the two spacelike curves as Σ .
As a matter of fact, the configurations on the final time slice has a bubble. It can be checked as follows; the value of ϕ c at the origin exceeds the top of the potential hill. In the case of Figure 7a where we set ( β , α g / H ) = ( 0.5 , 140 ) , we have ϕ c / α 0.1 > 0.5 = β at H r = 0 , H η 0.7 .
Now, we can evaluate the action for the bubble configurations for Σ (Figure 7b,d). We use the fitting function (84) for the numerical integration.
I Σ d 4 x a 3 ϕ Δ ϕ ¨ c + 3 H ϕ ˙ c H CDL = 6 π H 6 ε Σ d ρ d θ 1 d θ 2 d θ 3 sin 2 θ 1 sin θ 2 sinh 3 ρ ( cosh ρ sinh ρ cos θ 1 ) 4 η 2 ϕ c 2 η η ϕ c 2 = 24 π 2 H 2 ε Σ d ρ d θ 1 sin 2 θ 1 sinh 3 ρ [ ( sinh ρ cosh ρ cos θ 1 ) 2 ( cosh ρ sinh ρ cos θ 1 ) 2 2 ϕ c ρ 2 + 3 sinh ρ cosh ρ cos θ 1 cosh ρ sinh ρ cos θ 1 + sin 2 θ 1 tanh ρ ( cosh ρ sinh ρ cos θ 1 ) 2 ϕ c ρ ] 2 : = 24 π 2 α 2 H 2 ε I ˜ ( μ , ρ ¯ ) .
In the second line, we used H CDL = 0 , (74) and (75). Also, we performed a coordinate transformation ( t , r ) ( ρ , θ 1 ) . In the third line, the following relations were used
2 η ϕ c ( ρ ) η = 2 η ρ η ϕ c ( ρ ) ρ = 2 H 2 ( cosh ρ sinh ρ cos θ 1 ) ( sinh ρ cosh ρ cos θ 1 ) ϕ c ( ρ ) ρ ,
and
2 ϕ c ( ρ ) η 2 = H 2 [ ( sinh ρ cosh ρ cos θ 1 ) 2 2 ρ 2 + ( sinh ρ cosh ρ cos θ 1 ) ( cosh ρ sinh ρ cos θ 1 ) + sin 2 θ 1 tanh ρ ρ ] ϕ c ( ρ ) .
Here, the partial differentiation with respect to η was taken while keeping r constant.

5.2. Discussions

5.2.1. Appropriate Choice of ε

Since our results (86) depend on ε , we need to choose some specific value of ε to predict the tunneling rate. We emphasize that the ε dependence does not imply the pathology of our result. Rather our results should depend on ε because the nonperturbative effects from the UV modes are missed in our formalism, as discussed in Section 3.2. We then choose ε as large as possible to evaluate the tunneling rate.
Our analysis is valid only when the value of ε lies in a certain range. Below, we briefly discuss necessary conditions for ε to justify our analysis.

Lower Bound on ε

We need to choose a sufficiently large ε so that IR fields can describe the bubble configuration. This imposes that ( ε H ) 1 should be smaller than the typical physical size of the bubble or the bubble wall. This condition is satisfied by taking ε μ for our parameter choices.

Upper Bound on ε

Our analysis is based on Equation (33), which is invalid when the quantum fluctuations of IR fields are too large. This imposes an upper bound on ε because the size of the fluctuations in ϕ c at the scale k k c ( t ) is proportional to ε :
δ ϕ c | k k c : = d 3 k ( 2 π ) 3 δ ln ( k / k c ( t ) ) | ϕ ^ ( t , k ) | 2 1 / 2 k c / a 2 π = ε H 2 π ,
where we assumed ( k c / a ) 2 V ( ϕ c ) 9 and approximated ϕ as a massless scalar field.
As discussed in Section 3.2, (33) is justified when the term δ H is negligible. Then, we define the ratio Q of the term δ H to the final term on the RHS of (73) as10
Q : = δ H H 2 ε a 3 ϕ Δ 2 / ( 6 π ) ϕ Δ = 6 π H 2 ε ( a 2 2 ϕ c V ( ϕ c ) ) = 3 π 2 2 ε 2 H 4 V ( ϕ c ) a 2 2 ϕ c V ( ϕ c ) ,
and we impose Q < 1 . In the second equality, we used δ H = ( a 3 / 24 ) V ( ϕ c ) ϕ Δ 3 . Similarly to (87), we can estimate the size of the quantum fluctuations of a 2 2 ϕ . Then, we estimate Q at a large ε as Q g 2 α ε / H , where we also used | V ( ϕ c ) | g 2 α for | ϕ c | α . The condition ε < H / ( g 2 α ) follows from Q < 1 .

Summary

In summary, it is necessary to choose ε satisfying the following condition to justify our analysis,
μ ε ε max : = H / ( g 2 α ) 755 × 1.12 × 10 4 g 140 g α / H .
This suggests that, at least when g is sufficiently tiny, we can reliably choose a large ε . Note that we cannot make the action I given in (86) arbitrarily small even if we take ε max in the free-theory limit g 0 ; for ε = ε max , we have I | ε = ε max g 1 ( g α / H ) 3 I ˜ ( μ , ρ ¯ ) . Roughly speaking, the small g α / H corresponds to the bubble radius being larger than the Hubble radius. Since we are interested in the bubble being smaller than the Hubble radius, g α / H cannot be taken small arbitrarily.
A few more comments on the ε -dependence of our results are in order. Our result (86) shows that the value of I becomes larger for smaller choice of ε . This is because the size of the quantum fluctuations of modes k ε a ( t ) H is proportional to ε , while the nonperturbative physics caused by UV modes k > ε a ( t ) H is truncated in our formalism. However, this interpretation implicitly assumes that the modes k ε max a ( t ) H are relevant for the tunneling process. We expect that this is the case because the kinetic energy of the fluctuations in modes k ε max a ( t ) H is estimated as ( ε max H ) 4 , which is much smaller than the height of potential barrier Δ V g 2 α 4 for our parameter choices.
By contrast, very short scale physics which is insensitive to the detailed structure of V ( ϕ ) near its origin might be irrelevant for the tunneling process. Hence, we expect that the tunneling rate converges to some finite value in the limit ε . It would be interesting to see whether the tunneling rate converges to a finite value as ε increases, and if so, from which value of ε this convergence begins.

5.2.2. Hawking–Moss vs. Coleman–de Luccia

From the Formula (33), we derived the two configurations which describe the tunneling from the false vacuum to the true vacuum. One is the HM configuration, and the other is the CDL-like configuration. Using these configurations, we evaluated the probability of the tunneling. Here, we compare these probabilities. Denoting the action (86) and (52) as I bubble and I HM , the ratio of the two actions is given by
γ : = I bubble I HM = 9 α 2 H 2 Δ V ε I ˜ ( μ , ρ ¯ ) = 1 ε 108 ( 1 β ) 3 ( β + 3 ) H 2 α 2 g 2 I ˜ ( μ , ρ ¯ )
For example, if we choose β = 0.5 and α g / H = 140 (see Figure 6), the fitting parameters become μ = 7.32515 and ρ ¯ = 0.34397 , and the ratio becomes γ = 12.9840 / ε . For these parameters, the upper bound of ϵ is given by (88), and then γ can be smaller than one for an appropriately small g. This means that the CDL bubble configuration is dominant rather than the HM configuration. Results for other parameters are also shown in Table 1.

5.2.3. Comparison with the Euclidean Method

We also computed the CDL action in the Euclidean method (Appendix B) and compared the ratios of the CDL action to the HM action (Table 1). For ε > ε , our tunneling process becomes more probable than the one predicted by the Euclidean method. As discussed in Section 5.2.1, we can consistently choose ε > ε for certain choices of parameters ( g , α , β ) . Perhaps our results may indicate the presence of a tunneling process which is more probable than the Euclidean method.

5.2.4. Bubble Nucleation Hypersurface and the Subsequent Evolution

We found the non-trivial spacelike hypersurface on which the bubble was nucleated (Figure 7a), where the hypersurface was given by the condition F ( η , r ) = 0 or equivalently a 2 2 ϕ c V ( ϕ c ) = 0 . The subsequent evolution was then described by the one shown in the first line of (75), which is the standard classical equation of motion. The configuration on the hypersurface gave the initial data for the classical dynamics after the bubble nucleation.
The field value on the spacelike hypersurface does not reach the true vacuum but rather, it lies between the true vacuum and the top of the potential as mentioned in Section 5.1. This is possible because the location of the hypersurface is defined by the condition a 2 2 ϕ c V ( ϕ c ) = 0 , where the gradient force is balanced with the potential force. Though it is non-trivial to solve the evolution starting from the general spacelike hypersurfaces, it is desirable to solve it to fully understand the formation of the true-vacuum bubble in our scenario.
We conclude this section by pointing out that our formalism naturally predicts that the condition ϕ ¯ Δ = 0 (or ϕ Δ = 0 via (72)) defines the hypersurface on which the quantum dynamics is switched to the classical dynamics and the non-trivial field configuration is nucleated. This would be the generic prediction of our formalism that holds true once the appropriate configuration ϕ ¯ Δ is specified.

6. Conclusions

We studied the tunneling processes on the de Sitter background by using the stochastic approach. A novel point is that we applied the MSRJD functional integral to the problem. In this formalism, the tunneling rate was obtained by evaluating the functional integral with the saddle-point approximation. Using this method, we investigated both the HM and the CDL tunnelings.
For the HM case, we first analyzed (31), which is the MSRJD functional integral of the stochastic Equation (23) on a single spatial point. We succeeded in deriving the tunneling solution in the “phase space” (Figure 3) which represented the tunneling process from the false vacuum, through the top of the potential hill, and finally to the true vacuum (Figure 4). This solution had natural tunneling boundary conditions which could not be obtained from the Euclidean method. The tunneling rate for this configuration coincided with that of the HM instanton at the leading order. In Section 4.2, we succeeded in re-deriving this result starting from (33), which described the stochastic dynamics in the global region covered by the flat chart. We also estimated the time scale of the HM tunneling as (63), which became much longer than the Hubble time scale. Our analysis clarified the physical picture of the HM instanton, i.e., the HM transition probability represents the transition probability of a coarse-grained patch with a physical radius ( ε H ) 1 . Some technical complications that arose in dealing with (33) and the way of handling them were remarked on in Section 4.2.1. Physically, these complications were due to the spatial correlations of stochastic noises.
For the CDL case, we found the configurations which described the bubble nucleation process. The tunneling rates for the bubble configurations depended on the cut-off scale dividing IR and UV fields. We argued that this dependence came from the truncation of the nonperturbative effects from UV modes. We also discussed the valid choice of the cut-off scale for which the MSRJD functional integral would be reliable. This consideration was based on the first-principle derivation of the stochastic approach from the Schwinger–Keldysh formalism. With an appropriate cut-off scale, it turned out that the CDL tunneling rate was larger than the HM one for the steep potential barrier. However, as mentioned above (73), it has not yet been investigated if our configuration can really be used for estimating the tunneling action. This aspect would be important for evaluating the tunneling action, being left for future work. Nonetheless, we believe that our study clarifies how we can proceed with the analysis and how the bubble nucleation process could be described in the stochastic approach; for instance, our formalism can naturally define the location of the hypersurface on which the quantum dynamics is switched to the classical dynamics and the non-trivial field configuration is nucleated.
It is interesting to understand the relation between the stochastic approach and the Euclidean method. We compared the CDL tunneling rate in our method with that in the Euclidean method and found the former became larger than the latter for the potential considered with a certain cut-off scale. It would be worth investigating the meaning of this result. In particular, we need to know a precise relation between the CDL instanton and our configuration. Intriguingly, our method is also related to that discussed in [32,33,34,35,36], where the bubble accidentally appears from the initial quantum fluctuations. In our method, however, we assumed that we could neglect the nonperturbative effects from UV modes, the quantum component of the potential of the form V ( ϕ c ) ϕ Δ 3 , and that the configuration (74) could be used to estimate the tunneling rate. The first two were summarized into the term δ H . We also assumed a fixed-background spacetime. The relaxation of these assumptions would be important to seek the relations among the three methods. It would also be useful to apply our path integral method to the tunneling in flat space and make a comparison with the Euclidean method.
Once this relaxation is achieved, it would also be interesting to apply our method to tunneling phenomena where the backreaction to the background geometry is non-negligible and make comparisons with recently discussed methods. For instance, the formalism of the Wheeler–DeWitt equation [37,38,39] and the tunneling potential [40] must be related to ours because the HM exponent is reproduced. The tunneling in the black hole spacetime [41,42] is also interesting. In particular, [42] discussed the HM transition with a black hole from the viewpoint of the stochastic approach.
Our method using the saddle-point approximation effectively amounts to solving real-time quantum dynamics, which reduces to Starobinsky‘s stochastic approach [7,8] at leading order when the super-horizon dynamics in the de Sitter space is considered. However, since our starting point is the Schwinger–Keldysh path integral, which can be formulated in a more generic setup and incorporates all the quantum effects in principle, we believe that our formalism sheds light on further studies of tunneling phenomena from a real-time perspective.
One of the advantages of our method is that it is applicable to a dynamical setup. The application of our method to the inflation models such as the chain inflation [43] and the warm inflation [44] is also intriguing. Using our method, we can study the tunneling process in the inflationary background under dissipation and fluctuations coming from the circumference. We leave these issues for future work.

Author Contributions

Conceptualization, T.M., J.S. and J.T.; methodology, T.M., J.S. and J.T.; software, T.M.; investigation, T.M., J.S. and J.T.; writing—original draft, T.M., J.S. and J.T.; visualization, T.M. All authors have read and agreed to the published version of the manuscript.

Funding

T.M. was supported by JST SPRING, grant number JPMJSP2148 and JSPS KAKENHI grant number JP23KJ1543. J.S. was in part supported by JSPS KAKENHI grant numbers JP17H02894, JP17K18778, JP20H01902, JP22H01220. J.T. is supported by IBS under the project code IBS-R018-D1.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Stochastic Approach from the First Principle

In this section, we evaluate the generating functional for IR fields Z [ J IR ( T ) ] defined in Equation (38) to obtain (39). Note that each step we need to take for evaluating Z [ J IR ( T ) ] is summarized in Section 3.2.

Appendix A.1. Path Integral Representation

A path integral representation of Z [ J IR ( T ) ] takes the following standard form,
Z [ J IR ( T ) ] = D ( ϕ + , ϕ , Π + , Π ) e i J IR · ( ϕ + IR + ϕ IR ) / 2   × e i ( S H + S H ) x δ ( ϕ + ( T , x ) ϕ ( T , x ) ) Ψ 0 [ ϕ + ] Ψ 0 [ ϕ ] ,
where Ψ 0 [ ϕ ] is an initial wave functional. We specify it more concretely in Appendix A.2.1. S H ± = t = T d 4 x ( Π ± ϕ ˙ ± H [ ϕ ± , Π ± ] ) is a Hamiltonian action. After the rotation of the basis ( X + , X ) ( X c , X Δ ) : = ( X + + X 2 , X + X ) with X = ϕ , Π , (A1) is written as
Z [ J IR ( T ) ] = D ( ϕ c , ϕ Δ , Π c , Π Δ ) e i J IR · ϕ c IR e i S H [ ϕ c , ϕ Δ , Π c , Π Δ ] x δ ( ϕ Δ ( T , x ) ) ρ 0 [ ϕ c , ϕ Δ ] .
Here, we defined S H [ ϕ c , ϕ Δ , Π c , Π Δ ] : = ( S H + S H ) | ( X + , X ) ( X c , X Δ ) with X = ϕ , Π . We also defined ρ 0 [ ϕ c , ϕ Δ ] : = Ψ 0 [ ϕ c + ( ϕ Δ / 2 ) ] Ψ 0 [ ϕ c ( ϕ Δ / 2 ) ] .

Appendix A.2. Nonperturbative Generating Functional for IR Sector

Now, we explain how to calculate Z [ J IR ( T ) ] nonperturbatively in the IR sector. This allows us to capture the nonperturbative physics of the sector k k c ( t ) . Our strategy is to split the integration variables into UV modes k > k c ( t ) and IR modes k k c ( t ) for each time step t, and perform the integration over UV variables to obtain Z [ J IR ( T ) ] . One may simply split the integration variables X a ( t , k ) 11, where X = ( ϕ , Π ) and a = ( c , Δ ) , by the following replacement in (A2):
X a ( t , k ) X a IR ( t , k ) ( t t k ) X a UV ( t , k ) ( t < t k ) ,
where t k is defined by the condition k c ( t k ) = k . It is then tempting to perform the integration over UV variables X a UV by using the Schwinger–Keldysh (or closed-time-path) formalism. However, the expression obtained from (A2) after the replacement (A3) does not have the product of delta functions k δ ( ϕ Δ UV ( t k 0 + , k ) ) at the final time t = t k 0 + 12. That is, the time contour for each UV ϕ -variable with modes 0 < k k c ( T ) is not closed at the final time. Hence, the usual Schwinger–Keldysh formalism does not apply for evaluating the integration over UV variables. To resolve this issue, Refs. [22,23] proposed a new way of splitting integration variables. The splitting is given by (A3) for X = Π , while for ϕ -variable, it is defined by the following rule rather than (A3):
ϕ c ( t , k ) ϕ c IR ( t , k ) ( t > t k ) ϕ c UV ( t , k ) ( t t k ) , ϕ Δ ( t , k ) ϕ Δ IR ( t , k ) ( t t k ) ϕ Δ UV ( t , k ) ( t < t k ) .
In this replacement (A4), the variables ϕ Δ UV ( t k , k ) are absent while the variables ϕ c UV ( t k , k ) are present. Therefore, in the expression obtained from (A2) after the new replacement, the time contour for every UV ϕ -variable including those for modes k with 0 < k k c ( T ) is closed. This allows us to perform the integration over UV variables based on the Schwinger–Keldysh formalism as we see below. Note that there is no subtlety in the replacement of integration variables X a ( t , k ) for a zero mode k = 0 and those for deep UV modes k > k c ( T ) : in both replacement rules mentioned above, they are simply replaced by IR variables X a IR ( t , k ) and UV variables X a UV ( t , k ) , respectively. The time contours for them are closed thanks to the product of delta functions in the original expression (A2).
After the new replacement, the term S H [ ϕ c , ϕ Δ , Π c , Π Δ ] in the exponent of (A2) is decomposed into three pieces: purely IR terms S H IR : = S H [ ϕ c IR , ϕ Δ IR , Π c IR , Π Δ IR ] , purely UV terms S H UV : = S H [ ϕ c UV , ϕ Δ UV , Π c UV , Π Δ UV ] and the remaining IR–UV mixing terms S mixed :
S H [ ϕ c , ϕ Δ , Π c , Π Δ ] S H IR + S H UV + S mixed .
In terms of these quantities, Z [ J IR ( T ) ] can be formally written as [22,23]
Z [ J IR ( T ) ] = D ( ϕ c IR , ϕ Δ IR , Π c IR , Π Δ IR ) e i J IR · ϕ c IR e i S H IR e i Γ × k k c ( T ) δ ( ϕ Δ IR ( T , k ) ) 0 < k k c ( T ) δ ( ϕ c IR ( t k , k ) ) ,
where i Γ is the effective action for IR fields,
e i Γ : = D ( ϕ c UV , ϕ Δ UV , Π c UV , Π Δ UV ) e i ( S H UV + S mixed ) k > 0 δ ( ϕ Δ IR ( t f ( k ) , k ) ) ρ 0 [ ϕ c UV , ϕ Δ UV ; ϕ c IR , ϕ Δ IR ] ,
with t f ( k ) : = min [ T , t k ] . Here, ρ 0 [ ϕ c , ϕ Δ ] ρ 0 [ ϕ c UV , ϕ Δ UV ; ϕ c IR , ϕ Δ IR ] due to the replacement of variables. The IR variables ( ϕ c IR , ϕ Δ IR ) in the argument of ρ 0 contain only a zero mode. Since the time path for each UV ϕ -variable is closed in (A7), we can calculate i Γ perturbatively as usual by writing the connected diagrams with n external IR fields that are connected by UV propagators. The diagrammatic rules such as vertex factors and symmetry factors follow the standard rules of the Scheinger–Keldysh formalism for given vertexes in S H UV + S mixed .
Interestingly, due to the non-trivial replacement rule (A4), the Π ϕ ˙ terms in S H also contribute to S mixed [11,22,23,24]. Referring to such contributions and the remaining terms in S mixed as S tr and S mix - int , respectively, we have S mixed = S tr + S mix - int with 13
S tr : = T d t d 3 k ( 2 π ) 3 δ ( t t k ) ϕ Δ IR ( t , k ) Π c UV ( t , k ) Π Δ IR ( k , t ) ϕ c UV ( t , k ) .
The bi-linear vertexes exist only at the UV→IR transition time t k for given modes with 0 < k k c ( T ) . Physically, these vertexes set the initial fluctuations for such modes in (A6).

Appendix A.2.1. Integrate out Short-Wavelength Modes

Now, we calculate i Γ perturbatively by specifying the setup more concretely. We organize the perturbation theory around the false vacuum state ρ 0 for which the expectation value of the zero mode is ϕ false . We assume that the fluctuations of zero modes are tiny enough so that the couplings between zero-mode fluctuations and those of nonzero modes k > 0 are switched off in the past infinity t = . We then take the initial state of nonzero modes to be the Bunch–Davies vacuum state for a free field, where the time evolution of a free field is defined by the quadratic action expanded around the homogeneous background ϕ false . These assumptions can be implemented by writing the initial quantum state as
ρ 0 [ ϕ c UV , ϕ Δ UV ; ϕ c IR , ϕ Δ IR ] = ρ BD [ ϕ c UV , ϕ Δ UV ; ϕ false ] ρ 0 [ ϕ c IR ϕ false , ϕ Δ IR ]
which is properly normalized as
1 = k > 0 d ϕ c UV ( , k ) d ϕ Δ UV ( , k ) ρ BD [ ϕ c UV , ϕ Δ UV ; ϕ false ] = d ϕ c IR ( , 0 ) d ϕ Δ IR ( , 0 ) ρ 0 [ ϕ c IR ϕ false , ϕ Δ IR ] .
ρ BD [ ϕ c UV , ϕ Δ UV ; ϕ false ] is the density matrix for nonzero modes. An initial state for the zero mode is given by ρ 0 [ ϕ c IR ϕ false , ϕ Δ IR ] , for which we have
Ψ ( ) | ϕ ^ ( k = 0 ) | Ψ ( ) = d ϕ c IR ρ 0 [ ϕ c IR ϕ false , ϕ Δ IR ] ϕ c IR = ϕ false ( 2 π ) 3 δ ( 0 ) .
Combining this with Ψ ( ) | ϕ ^ ( k ) | Ψ ( ) = 0 for nonzero modes k > 0 , we have
Ψ ( ) | ϕ ^ ( k ) | Ψ ( ) = ϕ false ( 2 π ) 3 δ ( k ) .
Hence, the VEV of the IR field ϕ ^ IR ( t , x ) in the past infinity t = is given by ϕ false .
The interaction vertexes for calculating i Γ are the bi-linear vertexes (A8) and nonlinear vertexes coming from interacting potential V ( ϕ ) . We assume that the potential is sufficiently flat under the region we are interested in. It is then natural to separate i Γ into the leading-order pieces and higher-order terms in the coupling constant as follows,
e i Γ = ρ 0 [ ϕ c IR ϕ false , ϕ Δ IR ] e i Γ LO k > 0 e i d 4 x δ H higher , e i Γ LO k > 0 : = D ( ϕ c UV , ϕ Δ UV , Π c UV , Π Δ UV ) e i S H , free UV + S tr k > 0 δ ( ϕ Δ IR ( t f ( k ) , k ) ) ρ BD [ ϕ c UV , ϕ Δ UV ; ϕ false ] .
Here, S H , free UV [ ϕ c UV , ϕ Δ UV , Π c UV , Π Δ UV ; ϕ false ] is the free part of the Hamiltonian action defined around the false vacuum as explained above. Higher-order corrections to Γ in the coupling constant are represented by the term δ H higher , which can be evaluated perturbatively. We have δ H higher = 0 at the leading order.
To calculate i Γ LO k > 0 by writing down the connected diagrams with external IR fields, we only need to consider the bi-linear vertexes S tr given in (A8). The number of diagrams are only four, and i Γ LO k > 0 can be calculated as [22,23]
i Γ LO k > 0 = 1 2 d t d t d 3 k ( 2 π ) 3 X α ( t , k ) g α β ( t , t , k ) X β ( t , k ) = 1 2 d 4 x d 4 x X α ( x ) G α β ( x , x ) X β ( x ) .
Here, the Greek indices α , β = ( ϕ , Π ) label the IR- Δ fields: ( X ϕ , X Π ) = ( Π Δ IR , ϕ Δ IR ) . Substituting Equation (A13) into (A6), we find
Z [ J IR ( T ) ] = D ( ϕ c IR , ϕ Δ IR , Π c IR , Π Δ IR ) e i J IR · ϕ c IR ρ 0 [ ϕ c IR ϕ false , ϕ Δ IR ] e i ( S H IR + Γ LO k > 0 ) e i d 4 x δ H higher × k k c ( T ) δ ( ϕ Δ IR ( T , k ) ) 0 < k k c ( T ) δ ( ϕ c IR ( t k , k ) ) .
Now, we introduce i Γ k = 0 as
e i Γ k = 0 [ ϕ Δ IR , Π Δ IR ] : = d ϕ c IR ρ 0 [ ϕ c IR ϕ false , ϕ Δ IR ] e i Π Δ IR ( ϕ c IR ϕ false ( 2 π ) 3 δ ( 0 ) ) .
By definition, we have
d ϕ c IR ρ 0 [ ϕ c IR ϕ false , ϕ Δ IR ] e i Π Δ IR ϕ c IR ϕ false ( 2 π ) 3 δ ( 0 ) = d ϕ c IR e i Π Δ IR ϕ c IR e i Γ k = 0 [ ϕ Δ IR , Π Δ IR ] δ ϕ c IR ϕ false ( 2 π ) 3 δ ( 0 ) .
Substituting (A17) into (A15), we find
Z [ J IR ( T ) ] = D ( ϕ c IR , ϕ Δ IR , Π c IR , Π Δ IR ) e i J IR · ϕ c IR e i ( S H IR + Γ LO k > 0 + Γ k = 0 ) e i d 4 x δ H higher × k k c ( T ) δ ( ϕ Δ IR ( T , k ) ) δ ( ϕ c IR ( t k , k ) ϕ false ( 2 π ) 3 δ ( k ) ) ,
by defining t k | k = 0 = .

Appendix A.2.2. Long-Wavelength Sector

For later convenience, let us focus on the term d t Π c IR ( t , k ) ϕ ˙ Δ IR ( t , k ) in S H IR . We use the following trick to perform the integration by parts:
δ ( ϕ Δ IR ( T , k ) ) e i t k T d t Π c IR ( t , k ) ϕ ˙ Δ IR ( t , k ) = δ ( ϕ Δ IR ( T , k ) ) d Π c IR ( t k 0 + , k ) e i t k 0 + T d t Π c IR ( t , k ) ϕ ˙ Δ IR ( t , k ) δ ( Π c IR ( t k 0 + , k ) ) = δ ( ϕ Δ IR ( T , k ) ) d Π c IR ( t k 0 + , k ) e i t k 0 + T d t Π ˙ c IR ( t , k ) ϕ Δ IR ( t , k ) δ ( Π c IR ( t k 0 + , k ) ) .
Thanks to this trick, (A18) can be written as
Z [ J IR ( T ) ] = D ( ϕ c IR , ϕ Δ IR , Π c IR , Π Δ IR ) D ˜ Π c IR e i J IR · ϕ c IR e i ( S H , ( d ) IR + S H , ( s ) IR + Γ LO k > 0 + Γ k = 0 ) e i d 4 x δ H higher × k k c ( T ) δ ( ϕ Δ IR ( T , k ) ) δ ( ϕ c IR ( t k , k ) ϕ false ( 2 π ) 3 δ ( k ) ) δ ( Π c IR ( t k 0 + , k ) ) ,
where D ˜ Π c IR : = k k c ( T ) d Π c IR ( t k 0 + , k ) . Below, we simply omit this integration measure to simplify the notation. We also decompose S H IR as S H IR = S H , ( d ) IR + S H , ( s ) IR up to the boundary terms, which vanish thanks to the trick (A19) as
S H , ( d ) IR = d 4 x Π Δ IR ϕ ˙ c IR a 3 Π c IR ϕ Δ IR Π ˙ c IR a ( t ) 2 ϕ c IR + a 3 ( t ) V ( ϕ c IR ) ,
S H , ( s ) IR = d 4 x a 3 ( t ) V ( ϕ c IR + ( ϕ Δ IR / 2 ) ) V ( ϕ c IR ( ϕ Δ IR / 2 ) ) V ( ϕ c IR ) ϕ Δ IR .
Here, the spacetime argument of IR variables are omitted. We introduce the IR variables and their time derivatives as
X IR ( t , x ) : = d 3 k ( 2 π ) 3 θ ( k c ( t ) k ) X IR ( t , k ) e i k . x ,
X ˙ IR ( t , x ) : = d 3 k ( 2 π ) 3 θ ( k c ( t ) k ) X ˙ IR ( t , k ) e i k . x ,
where X IR represents the IR variables: X IR = ( ϕ c IR , ϕ Δ IR , Π c IR , Π Δ IR ) . Precisely speaking, for X IR = Π c IR , we should replace k c ( t ) by k c ( t + 0 + ) so that we have k c ( t ) = k at t = t k 0 + . One may be concerned about the subtlety of the time derivatives ϕ ˙ c IR ( t , x ) and Π ˙ c IR ( t , x ) in (A21): in general we have X ˙ IR ( t , x ) t ( X IR ( t , x ) ) , where the LHS is defined by (A24), due to the presence of a time-dependent step function. However, there is no subtlety thanks to the initial conditions for all IR modes k k c ( T ) of the form
ϕ c IR ( t k , k ) = ϕ false ( 2 π ) 3 δ ( k ) , Π c IR ( t k 0 + , k ) = 0 ,
which are imposed by the delta functions in the final line in (A20). The corresponding conditions in the real space for all x are
X ˙ IR ( t , x ) = t ( X IR ( t , x ) ) for X IR = ( ϕ c IR , Π c IR ) ,
ϕ c IR ( , x ) = ϕ false , Π c IR ( , x ) = 0 .
To implement these conditions in the path integral, we define δ C ini by (40b) with
δ C detail : = t T X IR = ( ϕ c IR , Π c IR ) x δ X ˙ IR ( t , x ) t ( X IR ( t , x ) ) .
We find from Equations (37) and (A20) that the transition probability can be written in terms of the path integral as
p ( { ϕ IR ( T , x ) } T , x D ) = D ( ϕ c IR , ϕ Δ IR , Π c IR , Π Δ IR ) e i S e i Γ k = 0 δ C fin δ C ini ,
where δ C fin is defined in (40a), and the term S in the exponent is given by
S = ( S H , ( d ) IR + S H , ( s ) IR + Γ LO k > 0 ) d 4 x δ H higher
which agrees with the definition of S given in (41). When the zero-mode fluctuations are negligible, we can set Γ k = 0 = 0 . This is the setup adopted in the main text. Then, (A28) coincides with (39).

Appendix A.3. Stochastic Interpretation of the Tunneling

It is useful to see how the stochastic equations emerge as a consequence of performing the path integral over IR variables nonperturbatively. This leads to the stochastic interpretation of the tunneling probability. A final result of this section is (A39).
For simplicity, let us ignore higher-order corrections δ H higher for a while. We start by introducing the auxiliary fields ( ξ ϕ , ξ Π ) as
e i Γ LO k > 0 + Γ k = 0 + S H , ( s ) IR = D ( ξ ϕ , ξ Π ) P [ ξ ϕ , ξ Π ] e i d 4 x ξ Π ( x ) ϕ Δ IR ( x ) ξ ϕ ( x ) Π Δ IR ( x ) .
We substitute (A30) into (A20). Then, the exponent becomes linear in ϕ Δ IR ( t , k ) or Π Δ IR ( t , k ) with t 0 t < T :
Z [ J IR ( T ) ] D ( ϕ c IR , ϕ Δ IR , Π c IR , Π Δ IR ) D ˜ Π c IR e i J IR · ϕ c IR D ( ξ ϕ , ξ Π ) P [ ξ ϕ , ξ Π ] × exp i T d t d 3 k ( 2 π ) 3 Π Δ IR ( t , k ) L ϕ ( t , k ) ϕ Δ IR ( t , k ) L Π ( t , k ) Θ ( k c ( t ) k ) × k k c ( T ) δ ( ϕ Δ IR ( T , k ) ) δ ( ϕ c IR ( t k , k ) ϕ false ( 2 π ) 3 δ ( k ) ) δ ( Π c IR ( t k 0 + , k ) ) ,
where
L ϕ ( t , k ) : = ϕ ˙ c IR ( t , k ) a 3 ( t ) Π c IR ( t , k ) ξ ϕ ( t , k ) ,
L π ( t , k ) : = Π ˙ c IR ( t , k ) + a ( t ) k 2 ϕ c IR ( t , k ) + a 3 ( t ) V ( ϕ c IR ) ( t , k ) ξ Π ( t , k ) .
Consequently, we can perform the path integral D ( ϕ Δ IR , Π Δ IR ) , leading to
Z [ J IR ( T ) ] k k c ( T ) d ϕ c IR ( T , k ) e i J IR · ϕ c IR D ( ξ ϕ , ξ Π ) P [ ξ ϕ , ξ Π ] L ϕ = L Π = 0 & A 25 .
Now, we regard auxiliary fields ( ξ ϕ , ξ Π ) as noise variables whose correlations are given by
ξ α ( t , k ) ξ β ( t , k ) : = D ( ξ ϕ , ξ Π ) P [ ξ ϕ , ξ Π ] ξ α ( t , k ) ξ β ( t , k ) = i δ δ X α ( t , k ) i δ δ X β ( t , k ) e i Γ LO k > 0 + Γ k = 0 + S H , ( s ) IR ϕ Δ IR = Π Δ IR = 0 ,
where α , β = ( ϕ , Π ) and ( X ϕ , X Π ) = ( Π Δ IR , ϕ Δ IR ) following the notation in (A14). We also have the corresponding expression in the momentum space. If we ignore S H , ( s ) and i Γ k = 0 , we have
ξ α ( t , k ) ξ β ( t , k ) = ( 2 π ) 3 δ k + k Re g α β ( t ) δ ( t t ) δ ( t t k ) .
Equation (A33) shows that the effective dynamics of IR fields are described by the set of Langevin equations L ϕ = L Π = 0 with the initial condition (A25) and the noise correlations (A35).
In the real space, the dynamics is described by the set of Langevin equations,
ϕ ˙ c IR ( x ) = a 3 ( t ) Π c IR ( x ) + ξ ϕ ( x ) ,
Π ˙ c IR ( x ) = a ( t ) 2 ϕ c IR ( x ) a 3 V ( ϕ c IR ( x ) ) + ξ Π ( x ) ,
with the initial conditions (A26) and the noise correlations
ξ α ( x ) ξ β ( y ) = G α β ( x , y ) .
Equation (A33) takes the following form in the real space:
Z [ J IR ( T ) ] x d ϕ c IR ( T , x ) e i J IR · ϕ c IR D ( ξ ϕ , ξ Π ) P [ ξ ϕ , ξ Π ] A 26 & A 36 .
Substituting (A38) into (37), we have
p ( { ϕ IR ( T , x ) } T , x D ) y d ϕ c IR ( T , y ) x D δ ( ϕ c IR ( T , x ) ϕ IR ( T , x ) )   × D ( ξ ϕ , ξ Π ) P [ ξ ϕ , ξ Π ] A 26 & A 36 .
This expression provides a stochastic interpretation of the transition probability; the dynamics of ϕ c IR are described by the set of Langevin equations, and the transition probability is simply understood as the probability to realize the field configuration ϕ c IR ( T , x ) = ϕ IR ( T , x ) in the domain x D at the final time T starting from the initial configuration ϕ c IR ( , x ) = ϕ false . In particular, the set of Langevin equations is given by (8) with the noise correlations (15) when we ignore the terms δ H higher , S H , ( s ) IR , and Γ k = 0 . Note that such terms simply correct the Langevin equations14 and do not invalidate the above stochastic interpretation.

Appendix B. Coleman–de Luccia Tunneling in Euclidean Method

In this appendix, we provide the CDL action in the de Sitter background with the Euclidean method.
Generally, the ways of Euclideanization for the de Sitter spacetime depend on the coordinates. In this paper, we discuss the field theory on the metric (1), and then we choose the following coordinate again15.
d s 2 = 1 H 2 η 2 ( d η 2 + d x 2 ) .
By taking the wick rotation as η i τ , the coordinate becomes
d s 2 = 1 H 2 τ 2 ( d τ 2 + d x 2 ) .
This is the Euclidean AdS space discussed in Section 5, and we can take the global coordinate
d s 2 = H 2 ( d ρ 2 + sinh 2 ρ d Ω 2 ) .
We also impose that the field ϕ only depends on ρ . Therefore, the Euclidean action I E and the equation of motion becomes
I E = 2 π 2 0 d ρ sinh 3 ρ 1 2 H 2 ( ρ ϕ ) 2 + V ( ϕ ) H 4 ,
d 2 ϕ d ρ 2 + 3 tanh ρ d ϕ d ρ = V ( ϕ ) H 2 .
The second equation is the same as (82), and we can obtain the same bounce solutions discussed in Section 5. For numerical calculations, we consider the potential (42) and fitting function (84). Also, we take the following non-dimensionalizations;
ϕ = α ϕ ˜ ,
V ( ϕ ) = g 2 α 4 V ˜ ( ϕ ˜ ) .
Therefore, the Euclidean action (A43) becomes
I E = 2 π 2 α 2 H 2 0 d ρ sinh 3 ρ 1 2 ( ρ ϕ ˜ ) 2 + g 2 α 2 H 2 V ˜ ( ϕ ˜ ) : = 2 π 2 α 2 H 2 I ˜ E .
The ratio of the Coleman–de Luccia Euclidean action I bounce E to the Hawking–Moss action I HM is given by
γ E = I bounce E I HM = H 2 g 2 α 2 9 I ˜ bounce E ( 1 β ) 3 ( β + 3 ) .

Notes

1
For earlier discussions of a stochastic approach based on the Schwinger–Keldysh formalism, see, e.g., [11,24].
2
Generally, we can choose other discretization like Stratonovich‘s discretization. This ambiguity does not affect the result since the amplitudes of the noise do not depend on fields within the range of our approximation.
3
Precisely, ϕ ^ IR ( x ) is defined in terms of the the Schrödinger picture field ϕ ^ ( k ) in momentum space as ϕ ^ IR ( x ) : = d 3 k ( 2 π ) 3 ϕ ^ ( k ) e i k · x θ ( k c ( t ) k ) . We have suppressed the trivial time dependence stemming from the step function.
4
There is a subtlety that modes satisfying k k c ( T ) are initially regarded as the “UV” degrees of freedom (DoFs) while they become “IR” DoFs due to the accelerating expansion of the spacetime. However, by adopting this splitting procedure, we can use the Schwinger–Keldysh (or closed-time-path) formalism to evaluate the integration over UV variables first as usual.
5
This is consistent with the observation made in [6] that an HM solution corresponds to the transition over a region of a Hubble horizon volume.
6
Note that this behavior may motivate one to replace the term j 0 ( k c ( t ) r ) in (55) by the step function
j 0 ( k c ( t ) r ) θ ( 1 k c ( t ) r ) .
This is the approximation adopted in [8].
7
Note that we can also estimate the action I by performing the coarse-graining in time suitably and substituting the smeared expression of Π Δ into the action:
I 8 π 2 3 H 4 t t d t V ( ϕ c ( t , x 0 ) ) d 3 x ϕ ˙ c ( t , x ) W f ( r ; t ) 8 π 2 3 H 4 t t d t ϕ ˙ c ( t , x 0 ) V ( ϕ c ( t , x 0 ) ) = 8 π 2 3 H 4 Δ V ,
reproducing the previous result (52). In the first approximate equality, we used H HM = 0 , which is well satisfied after the smearing. In the second one, we used ϕ c ( t , x ) ϕ c ( t , x 0 ) for r k c 1 ( t ) and the condition d 3 x W f ( r ; t ) = 1 which follows from the definition of W f . This estimate may be more analogous to the calculation of the action in Section 4.1 since the approximate locality of the dynamics is restored after the smearing.
8
This is the reason why we illustrate flow lines in the ( ϕ c , ϕ Δ ) -plane in Figure 5.
9
This is compatible with the condition ε < α / H when the weak coupling g 1 is considered.
10
We can also compare δ H with the second and the third term on the RHS of (73). Such considerations do not change our estimate of ε max in (88).
11
In the main text, we use the notation X k ( t ) to write variables in momentum space. However, we adopt the notation X ( t , k ) in this appendix since we have many subscripts such as c and Δ .
12
Here, 0 + is the infinitesimal time step which is introduced to obtain the path integral representation of the unitary time evolution as usual.
13
Precisely speaking, we need to perform the UV–IR splitting in the path integral with an infinitesimal discrete time step for deriving S tr correctly. After taking the continuum limit, we obtain (A8); see also [23]. Furthermore, the form of S mix - int depends on the model. In our case, we have
S mix - int = d 4 x [ V ( ϕ + UV + ϕ + IR ) V ( ϕ + UV ) V ( ϕ + IR ) V ( ϕ UV + ϕ IR ) + V ( ϕ UV ) + V ( ϕ IR ) ] ϕ ± UV = ϕ c UV ± ( ϕ Δ UV / 2 ) , ϕ ± IR = ϕ c IR ± ( ϕ Δ IR / 2 ) .
14
For instance, leading-order corrections to the stochastic dynamics at the super-horizon scales are calculated in [22].
15
The analytic continuation of this coordinate is examined in [5].

References

  1. Coleman, S.R. The Fate of the False Vacuum. 1. Semiclassical Theory. Phys. Rev. D 1977, 15, 2929–2936. [Google Scholar] [CrossRef]
  2. Callan, C.G., Jr.; Coleman, S.R. The Fate of the False Vacuum. 2. First Quantum Corrections. Phys. Rev. D 1977, 16, 1762–1768. [Google Scholar] [CrossRef]
  3. Coleman, S.R.; De Luccia, F. Gravitational Effects on and of Vacuum Decay. Phys. Rev. D 1980, 21, 3305. [Google Scholar] [CrossRef]
  4. Hawking, S.W.; Moss, I.G. Supercooled Phase Transitions in the Very Early Universe. Phys. Lett. B 1982, 110, 35–38. [Google Scholar] [CrossRef]
  5. Rubakov, V.A.; Sibiryakov, S.M. False vacuum decay in de Sitter space-time. Theor. Math. Phys. 1999, 120, 1194–1212. [Google Scholar] [CrossRef]
  6. Brown, A.R.; Weinberg, E.J. Thermal derivation of the Coleman-De Luccia tunneling prescription. Phys. Rev. D 2007, 76, 064003. [Google Scholar] [CrossRef]
  7. Starobinsky, A.A. Stochastic de Sitter (Inflationary) Stage in the Early Universe. Lect. Notes Phys. 1986, 246, 107–126. [Google Scholar] [CrossRef]
  8. Starobinsky, A.A.; Yokoyama, J. Equilibrium state of a selfinteracting scalar field in the De Sitter background. Phys. Rev. D 1994, 50, 6357–6368. [Google Scholar] [CrossRef] [PubMed]
  9. Goncharov, A.S.; Linde, A.D.; Mukhanov, V.F. The Global Structure of the Inflationary Universe. Int. J. Mod. Phys. A 1987, 2, 561–591. [Google Scholar] [CrossRef]
  10. Linde, A.D. Hard art of the universe creation (stochastic approach to tunneling and baby universe formation). Nucl. Phys. B 1992, 372, 421–442. [Google Scholar] [CrossRef]
  11. Tolley, A.J.; Wyman, M. Stochastic tunneling in DBI inflation. J. Cosmol. Astropart. Phys. 2009, 10, 006. [Google Scholar] [CrossRef]
  12. Noorbala, M.; Vennin, V.; Assadullahi, H.; Firouzjahi, H.; Wands, D. Tunneling in Stochastic Inflation. J. Cosmol. Astropart. Phys. 2018, 09, 032. [Google Scholar] [CrossRef]
  13. Hashiba, S.; Yamada, Y.; Yokoyama, J. Particle production induced by vacuum decay in real time dynamics. Phys. Rev. D 2021, 103, 045006. [Google Scholar] [CrossRef]
  14. Camargo-Molina, J.E.; Rajantie, A. Phase transitions in de Sitter: The stochastic formalism. arXiv 2022, arXiv:2204.02875. [Google Scholar] [CrossRef]
  15. Camargo-Molina, J.E.; González, M.C.; Rajantie, A. Phase Transitions in de Sitter: Quantum Corrections. arXiv 2022, arXiv:2204.03480. [Google Scholar] [CrossRef]
  16. Martin, P.C.; Siggia, E.D.; Rose, H.A. Statistical Dynamics of Classical Systems. Phys. Rev. A 1973, 8, 423–437. [Google Scholar] [CrossRef]
  17. De Dominicis, C.; Brezin, E.; Zinn-Justin, J. Field Theoretic Techniques and Critical Dynamics. 1. Ginzburg-Landau Stochastic Models Without Energy Conservation. Phys. Rev. B 1975, 12, 4945. [Google Scholar] [CrossRef]
  18. Janssen, H.K. On a Lagrangean for classical field dynamics and renormalization group calculations of dynamical critical properties. Z. Phys. B Condens. Matter 1976, 23, 377–380. [Google Scholar] [CrossRef]
  19. De Dominicis, C.; Peliti, L. Field Theory Renormalization and Critical Dynamics Above t(c): Helium, Antiferromagnets and Liquid Gas Systems. Phys. Rev. B 1978, 18, 353–376. [Google Scholar] [CrossRef]
  20. Elgart, V.; Kamenev, A. Rare event statistics in reaction-diffusion systems. Phys. Rev. E 2004, 70, 041106. [Google Scholar] [CrossRef]
  21. Altland, A.; Simons, B.D. Condensed Matter Field Theory, 2nd ed.; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar] [CrossRef]
  22. Tokuda, J.; Tanaka, T. Statistical nature of infrared dynamics on de Sitter background. J. Cosmol. Astropart. Phys. 2018, 02, 014. [Google Scholar] [CrossRef]
  23. Tokuda, J.; Tanaka, T. Can all the infrared secular growth really be understood as increase of classical statistical variance? J. Cosmol. Astropart. Phys. 2018, 11, 022. [Google Scholar] [CrossRef]
  24. Morikawa, M. Dissipation and Fluctuation of Quantum Fields in Expanding Universes. Phys. Rev. D 1990, 42, 1027–1034. [Google Scholar] [CrossRef] [PubMed]
  25. Feynman, R.P.; Vernon, F.L., Jr. The Theory of a general quantum system interacting with a linear dissipative system. Ann. Phys. 1963, 24, 118–173. [Google Scholar] [CrossRef]
  26. Weinberg, E.J. Classical Solutions in Quantum Field Theory: Solitons and Instantons in High Energy Physics; Cambridge Monographs on Mathematical Physics; Cambridge University Press: Cambridge, UK, 2012. [Google Scholar] [CrossRef]
  27. Weinberg, E.J. Hawking-Moss bounces and vacuum decay rates. Phys. Rev. Lett. 2007, 98, 251303. [Google Scholar] [CrossRef] [PubMed]
  28. Jensen, L.G.; Steinhardt, P.J. Bubble Nucleation and the Coleman-Weinberg Model. Nucl. Phys. B 1984, 237, 176–188. [Google Scholar] [CrossRef]
  29. Jensen, L.G.; Steinhardt, P.J. Bubble Nucleation for Flat Potential Barriers. Nucl. Phys. B 1989, 317, 693–705. [Google Scholar] [CrossRef]
  30. Hackworth, J.C.; Weinberg, E.J. Oscillating bounce solutions and vacuum tunneling in de Sitter spacetime. Phys. Rev. D 2005, 71, 044014. [Google Scholar] [CrossRef]
  31. Batra, P.; Kleban, M. Transitions Between de Sitter Minima. Phys. Rev. D 2007, 76, 103510. [Google Scholar] [CrossRef]
  32. Braden, J.; Johnson, M.C.; Peiris, H.V.; Pontzen, A.; Weinfurtner, S. New Semiclassical Picture of Vacuum Decay. Phys. Rev. Lett. 2022, 123, 031601. [Google Scholar] [CrossRef]
  33. Blanco-Pillado, J.J.; Deng, H.; Vilenkin, A. Flyover vacuum decay. J. Cosmol. Astropart. Phys. 2019, 12, 001. [Google Scholar] [CrossRef]
  34. Hertzberg, M.P.; Rompineve, F.; Shah, N. Quantitative Analysis of the Stochastic Approach to Quantum Tunneling. Phys. Rev. D 2020, 102, 076003. [Google Scholar] [CrossRef]
  35. Tranberg, A.; Ungersbäck, G. Bubble nucleation and quantum initial conditions in classical statistical simulations. J. High Energy Phys. 2022, 09, 206. [Google Scholar] [CrossRef]
  36. Hertzberg, M.P.; Yamada, M. Vacuum Decay in Real Time and Imaginary Time Formalisms. Phys. Rev. D 2019, 100, 016011. [Google Scholar] [CrossRef]
  37. Kristiano, J.; Lambaga, R.D.; Ramadhan, H.S. Coleman-de Luccia Tunneling Wave Function. Phys. Lett. B 2019, 796, 225–229. [Google Scholar] [CrossRef]
  38. Cespedes, S.; de Alwis, S.P.; Muia, F.; Quevedo, F. Lorentzian vacuum transitions: Open or closed universes? Phys. Rev. D 2021, 104, 026013. [Google Scholar] [CrossRef]
  39. Maniccia, G.; De Angelis, M.; Montani, G. WKB Approaches to Restore Time in Quantum Cosmology: Predictions and Shortcomings. Universe 2022, 8, 556. [Google Scholar] [CrossRef]
  40. Espinosa, J.R. Fresh look at the calculation of tunneling actions including gravitational effects. Phys. Rev. D 2019, 100, 104007. [Google Scholar] [CrossRef]
  41. Gregory, R.; Moss, I.G.; Oshita, N. Black Holes, Oscillating Instantons, and the Hawking-Moss transition. J. High Energy Phys. 2020, 07, 024. [Google Scholar] [CrossRef]
  42. Gregory, R.; Moss, I.G.; Oshita, N.; Patrick, S. Hawking-Moss transition with a black hole seed. J. High Energy Phys. 2020, 09, 135. [Google Scholar] [CrossRef]
  43. Freese, K.; Spolyar, D. Chain inflation: ‘Bubble bubble toil and trouble’. J. Cosmol. Astropart. Phys. 2005, 07, 007. [Google Scholar] [CrossRef]
  44. Berera, A. Warm inflation. Phys. Rev. Lett. 1995, 75, 3218–3221. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The schematic picture of the potential V ( ϕ ) which has a false vacuum at ϕ = ϕ false , true vacuum at ϕ = ϕ true and the top of the potential at ϕ = ϕ top .
Figure 1. The schematic picture of the potential V ( ϕ ) which has a false vacuum at ϕ = ϕ false , true vacuum at ϕ = ϕ true and the top of the potential at ϕ = ϕ top .
Universe 10 00292 g001
Figure 2. The potential of (42) for various β ’s. From top to bottom, the curves correspond to β = 0 , 0.25 , 0.5 , 0.75 , 1 .
Figure 2. The potential of (42) for various β ’s. From top to bottom, the curves correspond to β = 0 , 0.25 , 0.5 , 0.75 , 1 .
Universe 10 00292 g002
Figure 3. The Hamilton flow for (48) with α = H , g = 0.4 and β = 0.1 . The solid lines correspond to H ( ϕ c , Π Δ ) = 0 . From left to right, the intersections of the solid lines correspond to the false vacuum, the top of the potential hill and the true vacuum.
Figure 3. The Hamilton flow for (48) with α = H , g = 0.4 and β = 0.1 . The solid lines correspond to H ( ϕ c , Π Δ ) = 0 . From left to right, the intersections of the solid lines correspond to the false vacuum, the top of the potential hill and the true vacuum.
Universe 10 00292 g003
Figure 4. The tunneling configuration for the homogeneous fields is plotted for β = 0.1 and α g / H = 10 . We choose the initial condition as ϕ c ( 10 15 ) / α = 1 + 10 5 and solve the second equation of (50). When ϕ c / α reach the top of potential ϕ c / α = β , we switch to the first equation of (50) with the initial condition ϕ c / α = β + 10 5 .
Figure 4. The tunneling configuration for the homogeneous fields is plotted for β = 0.1 and α g / H = 10 . We choose the initial condition as ϕ c ( 10 15 ) / α = 1 + 10 5 and solve the second equation of (50). When ϕ c / α reach the top of potential ϕ c / α = β , we switch to the first equation of (50) with the initial condition ϕ c / α = β + 10 5 .
Universe 10 00292 g004
Figure 5. A schematic picture of Hamiltonian flows in the “field” phase space. We expect a flow line which starts from the false vacuum (left intersection) and reaches the bubble configuration (right intersection) through a nonzero ϕ Δ ( x ) path.
Figure 5. A schematic picture of Hamiltonian flows in the “field” phase space. We expect a flow line which starts from the false vacuum (left intersection) and reaches the bubble configuration (right intersection) through a nonzero ϕ Δ ( x ) path.
Universe 10 00292 g005
Figure 6. The tunneling configuration for the inhomogeneous field. The solid curve is the numerical solution for β = 0.5 and α g / H = 140 . We took the initial conditions as ϕ c ( 10 15 ) / α = 1 4.98284 × 10 2 and ϕ c ( 10 15 ) / α = 0 . The dashed curve is the fitting function (84) where μ = 7.32515 and ρ ¯ = 0.34397 .
Figure 6. The tunneling configuration for the inhomogeneous field. The solid curve is the numerical solution for β = 0.5 and α g / H = 140 . We took the initial conditions as ϕ c ( 10 15 ) / α = 1 4.98284 × 10 2 and ϕ c ( 10 15 ) / α = 0 . The dashed curve is the fitting function (84) where μ = 7.32515 and ρ ¯ = 0.34397 .
Universe 10 00292 g006
Figure 7. (a): The two non-trivial curves on which (85) is satisfied. We inserted the mesh to confirm whether the curves are spacelike or timelike. It can be seen that the upper line is spacelike but the lower one is partially timelike. (b): The integration region for (86). There is another spacelike curve at η = . (c): The two non-trivial curves in ( ρ , θ 1 ) plane. (d): The integration region for (86) in the ( ρ , θ 1 ) plane. All plots are with β = 0.5 and α g / H = 140 .
Figure 7. (a): The two non-trivial curves on which (85) is satisfied. We inserted the mesh to confirm whether the curves are spacelike or timelike. It can be seen that the upper line is spacelike but the lower one is partially timelike. (b): The integration region for (86). There is another spacelike curve at η = . (c): The two non-trivial curves in ( ρ , θ 1 ) plane. (d): The integration region for (86) in the ( ρ , θ 1 ) plane. All plots are with β = 0.5 and α g / H = 140 .
Universe 10 00292 g007
Table 1. The values of γ for several parameter sets. The first and second ones are potential parameters. The third one is the initial condition for the equation of motion (82). The fourth and fifth ones are fitting parameters for (84). The sixth and seventh ones are the numerical results of (89) and (A48). The last one is the cut-off parameter ϵ for which (89) coincides with (A48).
Table 1. The values of γ for several parameter sets. The first and second ones are potential parameters. The third one is the initial condition for the equation of motion (82). The fourth and fifth ones are fitting parameters for (84). The sixth and seventh ones are the numerical results of (89) and (A48). The last one is the cut-off parameter ϵ for which (89) coincides with (A48).
β g 2 α 2 / H 2 1 ϕ c ( 10 15 ) / α μ ρ ¯ ε γ γ E ϵ
0.3804.87210 × 10 5 6.225631.05197355.0890.2470041437.58
0.4805.20273 × 10 3 5.996450.6525764.11850.0909796704.758
0.5803.89106 × 10 2 5.634440.47524733.95260.0599470566.376
0.31403.50945 × 10 4 8.139900.67330552.25150.04866781073.64
0.41408.96003 × 10 3 7.824190.45965919.29120.0239903804.123
0.51404.98284 × 10 2 7.325150.3439712.98400.0172063754.606
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Miyachi, T.; Soda, J.; Tokuda, J. Stochastic Tunneling in de Sitter Spacetime. Universe 2024, 10, 292. https://doi.org/10.3390/universe10070292

AMA Style

Miyachi T, Soda J, Tokuda J. Stochastic Tunneling in de Sitter Spacetime. Universe. 2024; 10(7):292. https://doi.org/10.3390/universe10070292

Chicago/Turabian Style

Miyachi, Taiga, Jiro Soda, and Junsei Tokuda. 2024. "Stochastic Tunneling in de Sitter Spacetime" Universe 10, no. 7: 292. https://doi.org/10.3390/universe10070292

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop