Emergent Universe by Tunneling in a Jordan-Brans-Dicke Theory

In this work we study an alternative scheme for an Emergent Universe scenario in the context of a Jordan-Brans-Dicke theory, where the universe is initially in a truly static state supported by a scalar field located in a false vacuum. The model presents a classically stable past eternal static state which is broken when, by quantum tunneling, the scalar field decays into a state of true vacuum and the universe begins to evolve following the extended open inflationary scheme.


I. INTRODUCTION
The standard cosmological model (SCM) [1][2][3] and the inflationary paradigm [4][5][6][7] are shown as a satisfactory description of our universe [1][2][3]. However, despite its great success, there are still important open questions to be answered. One of these questions is whether the universe had a definite origin, characterized by an initial singularity or if, on the contrary, it did not have a beginning, that is, it extends infinitely to the past.
Theorems about spacetime singularities have been developed in the context of inflationary universes, proving that the universe necessarily has a beginning. In other words, according to these theorems, the existence of an initial singularity can not be avoided even if the inflationary period occurs, see Refs. [8][9][10][11][12]. In theses theorems it is demonstrated that null and time-like geodesics are generally incomplete in inflationary models, regardless of whether energy conditions are maintained, provided that the average expansion condition (H > 0) is maintained throughout of these geodesics directed towards the past, where H is the Hubble parameter.
In the EU scheme it is assumed that the universe emerged from a past eternal Einstein Static (ES) state to the inflationary phase and then evolves into a hot big bang era. These models do not satisfy the geometrical assumptions of the theorems [8][9][10][11][12] and they provide specific examples of nonsingular inflationary universes.
We can note that in these schemes of Emergent Universe are not all truly static during the static regime.
At this respect a new scheme for an EU model was proposed in Ref. [39], where the universe is initially in a truly static state supported by a scalar field located in a false vacuum, also see Refs. [40,41]. The universe begins to evolve when, by quantum tunneling, the scalar field decays into a state of true vacuum. For simplicity, in this first approach to this new scheme of EU, the model was developed in the context of General Relativity (GR).
In particular in Ref. [39] was concluded that this new mechanism for an Emergent Universe is plausible and could be an interesting alternative to the realization of the Emergent Universe scenario.
However, as this first model was developed in the context of General Relativity, the past eternal static period, suffer from classical instabilities associated with the instability of Einstein's static universe.
The ES solution is unstable to homogeneous perturbations, as was early discussed by Eddington in Ref. [42] and more recently studied in Refs. [43][44][45][46]. The instability of the ES solution ensures that any perturbations, no matter how small, rapidly force the universe away from the static state, thereby aborting the EU scenario.
This instability is possible to cure by going away from GR, for example, by consider a Jordan-Brans-Dicke (JBD) theory at the classical level, where it have been found that contrary to general relativity, a static universe could be classically stable, see Refs. [21,22,47].
In this work, we are interested in apply the scheme of Emergent Universe by Tunneling of Ref. [39] to EU models which present stable past eternal static regimes. In particular, we study this scheme in the context of a JBD theory, similar to the one studied in Refs. [21,22], but where the static solution is supported by a scalar field located in a false vacuum. In this case we are going to show that, contrary to what happens in Ref. [39], the ES solution is classically stable. The Jordan-Brans-Dicke [48] theory is a class of models in which the effective gravitational coupling evolves with time. The strength of this coupling is determined by a scalar field, the so-called Brans-Dicke field, which tends to the value G −1 , the inverse of the Newton's constant. The origin of Brans-Dicke theory is found in the Mach's principle according to which the property of inertia of material bodies arises from their interactions with the matter distributed in the universe. In modern context, Brans-Dicke theory appears naturally in supergravity models, Kaluza-Klein theories and in all known effective string actions [49][50][51][52][53][54][55].
In particular in this work we are going to consider that the universe is initially in a truly static state, which is supported by a scalar field ψ located in a false vacuum (ψ = ψ F ), see Fig. (2). The universe begins to evolve when, by quantum tunneling, the scalar field decays into a state of true vacuum. Then, a small bubble of a new phase of field value ψ W can be formed, and expands as it converts volume from high to low vacuum energy and feeds the liberated energy into the kinetic energy of the bubble wall. This process was first studied by Coleman & De Luccia in [56,57] in the context of General Relativity.
If the potential has a suitable form, inflation and reheating may occur inside the bubble as the field rolls from ψ W to the true minimum at ψ T , in a similar way to what happens in models of Open Inflationary Universes and Extended Open Inflationary Universes, see for example [58][59][60][61][62], where the interior of the bubble is modeled by an open Friedmann-Robertson-Walker universe.
The advantage of the EU by Tunneling scheme (and of the Emergent Universe in general), over the Eternal Inflation scheme is that it correspond to a realization of a singularity-free inflationary universe. As was discussed in Refs. [8][9][10][11][12], Eternal Inflation is usually future eternal but it is not past eternal, because in general space-time that allows for inflation to be future eternal, cannot be past null complete. On the other hand Emergent Universe are geodesically complete.
Notice that in the EU by Tunneling scheme, the metastable state which support the initial static universe could exist only a finite amount of time. Then, in this scheme of Emergent Universe, the principal point is not that the universe could have existed an infinite period of time, but that in theses models the universe is non-singular because the background where the bubble materializes is geodesically complete.
This implies that we have to consider the problem of the initial conditions for a static universe. Respect to this point, there are very interesting possibilities discussed for example in the early works on EU [14] and more recently in [39]. One of these options is to explore the possibility of an Emergent Universe scenario within a string cosmology context [63]. Other possibility is that the initial Einstein Static universe is created from "nothing" [64,65], see Refs. [66] for explicit examples. It is interesting to mention that the study of the Einstein Static solution as a preferred initial state for our universe have been considered in the past, where it has been proposed that entropy considerations favor the ES state as the initial state for our universe [44,45].
In this paper we consider a simplified version of this scheme, with the focus on studying the process of creation and evolution of a bubble of true vacuum in the background of an ES universe in the context of a JBD theory. This is motivated because we are mainly interested in the study of new ways of leaving the static period and begin the inflationary regime for Emergent Universe models which present a classically stable static state period.
In particular, in this paper we consider a JBD theory where one of the matter content of the model is a scalar field (inflaton) with a potential similar to Fig. (3) and we study the process of tunneling of the scalar field from the false vacuum U F to the true vacuum U T and the subsequent creation and evolution of a bubble of true vacuum in the background of an stable ES universe. The simplified model studied here contains the essential elements of the scheme we want to present, so we postpone the detailed study of the inflationary period, which occurs after the tunneling, for future work.
The paper is organized as follow. In Sect. II we study a Einstein static universe supported by a scalar field located in a false vacuum and its stability in the context of a JBD theory.
In Sect. III we study the tunneling process of the scalar field from the false vacuum to the true vacuum and the subsequent creation of a bubble of true vacuum in the background of the Einstein static universe for a JBD theory. In Sect. IV we study the evolution of the bubble after its materialization. In Sect. V we summarize our results.

II. FALSE VACUUM AND THE ES STATE IN JBD THEORIES
In this paper we consider a scheme for EU scenario where the universe is initially in a classically stable static state. This state is supported by a scalar field located in a false vacuum. The universe begins to evolve when, by quantum tunneling, the scalar field decays to a state of true vacuum. For this reason we will begin by studying the possibility of obtain a static and classically stable ES solution in this theory when the scalar field is in a false vacuum.
We consider the following JBD action for a self-interacting potential and matter, given by [67] where R is the Ricci scalar curvature, φ is the JBD field, ω is the JBD parameter, V (φ) is the potential associated to the JBD field, ψ is the scalar field (inflaton), U (ψ) is the scalar potential and L m denotes the Lagrangian density of a barotropic perfect fluid. In this theory 1/φ plays the role of the gravitational constant, which changes with time. This action also matches the low energy string action for ω = −1, see [55].
Following the EU scheme we consider a closed Friedmann-Robertson-Walker metric: where a(t) is the scale factor and t represents the cosmological time. The content of matter is modeled by a standard perfect fluid with an effective state equation given by P f = (γ − 1) ρ f , with γ constant, and a scalar field (inflaton) for which The scalar field potential U (ψ) es depicted in Fig.(2). The global minimum of U (ψ) is U T , but there is also a local minimum U F (the false vacuum).
We have considered that the early universe is dominated by two fluids because in our scheme of EU scenario, during the static regime, the scalar field ψ remains static at the false vacuum, in contrast to standard EU models where the scalar field rolls on the asymptotically flat part of the scalar potential [22]. Then, in order to obtain a static universe under these conditions, we have to included a standard perfect fluid. For simplicity we will consider that there are no interactions between the standard perfect fluid and the scalar field. The

Friedmann-Raychaudhuri equations become
The field equation for the JBD field is where V = dV (φ)/dφ, ρ = ρ f + ρ ψ , P = P f + P ψ and dots represent derivatives with respect to cosmological time.
The conservation equations for the scalar field and perfect standard fluid arë respectively.
The static universe is characterized by the conditions a = a 0 = Const.,ȧ 0 =ä 0 = 0 and (4), (5) and (6) the static solution for a universe dominated by a scalar field placed in a false vacuum and a standard perfect fluid is obtained if the following conditions are satisfied: and where V 0 = V (φ 0 ) and V 0 = (dV (φ)/dφ) φ=φ 0 and ρ f 0 is the energy density of the perfect fluid present in the static universe. These equations connect the equilibrium values of the scale factor and the JBD field with the energy density and the JBD potential at the equilibrium point. We now study the stability of this solution against small homogeneous and isotropic perturbations. In order to do this, we consider small perturbations around the static solutions for the scale factor, JBD field and inflaton field.
We set Then, we have where α(t) 1, β(t) 1 and λ(t) 1 are small perturbations. On the other hand, we have that at linear order in λ where δψ = ψ F λ(t) and δψ = ψ 0λ (t). As in our case (inflaton field in a local minimum of its potential)ψ F = 0 and U F = 0 then we obtain that δρ ψ = δP ψ = 0, at linear order. Then, we have that the evolution of λ is determined only by the following equation obtained from Eq. (7) and Eq. (14)λ where U F = (d 2 U (ψ)/dψ 2 ) ψ=ψ F . We can note that the evolution of λ is completely decoupled from the evolution of α and β and correspond to an oscillatory behavior given that U F > 0 in our case.
By introducing the expressions (12), (13) and (15) into equations (5) and (6), and keeping terms at the linear order in α, β and λ, we obtain the following set of coupled equations and From the system of equations (19) and (20) we can obtain the frequencies for small oscilla- Note that the static solution is stable if the inequality w 2 ± > 0 is satisfied. Then assuming that the parameter ω satisfies the constraint, (3 + 2ω) > 0, we find that the following inequalities must be fulfilled in order to have a stable static solution and From these inequalities we can conclude that for a universe dominated by a scalar field in a false vacuum and a perfect standard fluid, it is possible to find a solution where the universe is static and stable. In particular, we study numerical solutions to the equations (5), (6), (7) and (8), for initial conditions close to the static solution. We address the case in which the parameters take the following values that satisfy the stability conditions (22), (23) and (24) γ = 1; φ 0 = 0.5; a 0 = 1; where we have used units in which 8πG = 1. From the conditions of the static solution (9), (10), (11) and the stability condition (24) we obtain the values for V 0 , V 0 and V 0 .
In particular, for the numerical solution, we consider the following inflaton potential where we have used that around the local equilibrium point ψ ∼ ψ F the inflaton potential potential can be approximated to a quadratic function.
Also, we consider the following JBD potential which satisfies the static and stability conditions discussed above In Fig. (4) it is shown one numerical solution corresponding to a universe starting from an initial state not in the static solution but close to it. We can note that while the scalar field makes small oscillations around the false vacuum, the scale factor, JBD field and perfect fluid density present small oscillations around their equilibrium values. This tells us that the static solution is classically stable.
Then, in a purely classical field theory if the universe is static and supported by the scalar field located at the false vacuum U F , the universe remains static forever. Quantum mechanics makes things different because the scalar field ψ can tunnel through the barrier and by this process create a small bubble where the field value is ψ T . Depending on the background where the bubble materializes the bubble could expand or collapse, see Refs. [39,68,69].

III. BUBBLE NUCLEATION
In this section we study the tunneling process of the scalar field ψ from the false vacuum U F to the true vacuum U T , in the potential U (ψ) shown in Fig. (3) and the consequent creation of a bubble of true vacuum in the background of an Einstein static universe, in the context of a JBD Theory. In particular, we will consider the nucleation of a spherical bubble of true vacuum within the false vacuum. We will assume that the layer which separates the two phases (the wall) has a negligible thickness compared to the size of the bubble (the usual thin-wall approximation). The energy budget of the bubble consists of latent heat (the difference between the energy densities of the two phases) and surface tension.
In order to study the tunneling process in this model, we have to deal with the nonstandard gravitational interaction of the JBD theory. This was done in Ref. [70] and here we reproduce and adapt these results to the EU scheme. Regarding the study of bubble nucleation in JBD theories in other contexts see Refs. [71,72].
We begin with the JBD action (1) described in Sect. II and perform a weyl rescaling in order to remove the non-standard coupling of φ to R.
where φ −1 G = G, here G is the value of the Newton constant observed in the present. We choose the conformal factor Ω(x) to be Ω = √ 8πGφ. Then, the action in the rescaled theory can be expressed as follow where we have defined In the rescaled theory the probability of nucleation per unit of physical volume per time is given byΓ where P (t) is the bubble probability of nucleation andV (t) is the physical volume measured in the rescaled system at timet. Therefore, the equivalent rate in the original theory is given by where we have used that P (t) = P [t(t)], see Ref. [70].
In the rescaled formulation of the theory, the gravitational interaction has the standard form, and so we can expect gravitational effects similar to those of the standard theory.
Nevertheless as it is usually, in order to eliminate the problem of predicting the reaction of the geometry to an essentially a-causal quantum jump (the materialization of the bubble), we neglect during this computation the gravitational back-reaction of the bubble onto the space-time geometry. The gravitational back-reaction of the bubble will be consider in the next Section when we study the evolution of the bubble after its materialization. Thus, in this approach only the matter field ψ remains as a dynamical quantity in the action (30).
Then, by following Ref. [70] we determinate the rate of nucleation of a bubble of true vacuum in the remaining dynamical theory for the field ψ. The nucleation rateΓ is given bȳ where S E is the Euclidean action.
It is shown in Ref. [70] that the coefficient A has a net factor of b 2 with respect to a theory in which b = 1, this entails thatΓ(t) = b 2 Γ 0 , where Γ 0 is the nucleation rate for a normal scalar field theory with potential U (ψ). Since b is a function of the time-dependent JBD field, then we have a time-dependent nucleation rate in the rescaled theory. However, this time dependence disappears if we return to the original theory. The equations (28) and (29), together with the definition of φ G and Ω imply that b = Ω 2 . Using the equation (34) we find that the nucleation rate of the original theory is, see [70] The computation of Γ 0 in the context of the EU scheme, that is, the tunneling process of the scalar field ψ in a background of Einstein static universe in the context of General Relativity, was developed in Refs. [39][40][41], where it was found that where ε is the difference of energy density between the two phases (latent heat) and σ is the energy density of the wall.

IV. CLASSICAL EVOLUTION OF THE BUBBLE ON JORDAN-BRANS-DICKE THEORIES
In this section we study the evolution of true vacuum bubble after its nucleation via quantum tunnel effect. During this study we are going to consider the gravitational backreaction of the bubble.
We follow the approach and notation used in Ref. [73] where it is assumed that the bubble wall separates spacetime into two parts. The bubble wall is a timelike, spherically symmetric hypersurface, the interior of the bubble is our universe, according to the extended open inflation scheme [58][59][60][61][62] and the exterior correspond to the static universe discussed in Sec. II. In particular, we will use the scheme developed in Ref. [74] regarding the Darmois-Israel junction conditions [75,76] applied to a JBD Theory.

Let us start by summarize the formalism developed by Berezin, Kuzmin and Tkachev in
Ref. [73] for the analysis of a thin wall bubble in the context of General Relativity before to study this problem in the context of a Jordan-Brans-Dicke-Theory.
We begin by defining a time-like spherically symmetric hypersurface Σ, representing the world surface of the wall dividing space-time into two regions, V + (outside) and V − (inside).
Also, we define a space-like vector N µ , which is orthogonal to Σ and points from V − to V + .
It is convenient to introduce a normal Gaussian coordinate system (n, x i ) such that the hypersurface n = 0 corresponds to Σ. If we assume that the wall is infinitely thin, its surface energy-momentum tensor is Using the extrinsic curvature tensor defined by K ij ≡ N i;j and the Einstein's equations we can express the joint conditions on the wall as follows, see [75] [ where κ 2 ≡ 8πG, h ij is the 3-metric on Σ, and | denotes the three-dimensional covariant derivative. We use the brackets to represent the difference between the outer and inner values of any field variable, for example: [F ] ± ≡ F + − F − . By using Eq.(39) and Eq. (41) we eliminate K i − j obtaining: Since Eq. (40) and Eq.(42) do not contain K i − j , they may be useful when the geometry of interior of the bubble is unknown.
We are going to assume that the bubble wall Σ is spherically symmetric, then the intrinsic metric on the shell Σ is given by where R(τ ) is the radius of circumference of the wall and τ the proper time of the wall. The surface energy-momentum tensor S ij can be written as a perfect fluid: where σ,ω and u i = (1, 0, 0) are the surface energy density, surface pressure and a unit time-like vector tangent to Σ, respectively. Then we can rewrite Eq. (40) and Eq.(42) as As was discussed in Ref. [73], Eq.(45) and Eq.(46) determine the evolution of σ and R, in the context of General Relativity, as they give us information on how its radius evolves and about the energy density that accumulates on the outer side of the wall.
In our case, in order to derive the equations of motion of the bubble wall in the Jordan-Brans-Dicke theory we are going to follow a scheme similar to that discussed above, see Ref. [74].
From the JBD action Eq. (1) we obtained the following field equations whereT and primes denote derivatives with respect to φ. In this context it is useful to introduce the following surface energy-momentum tensor defined over Σ where the second equality is obtained from Eq.(48) and Eq. (49). With these definitions we can use the formalism described above (for General Relativity) in order to study the evolution of the bubble in the Jordan-Brans-Dicke context, we only have to replace T µν , S ij and κ 2 byT µν ,S ij y 1/φ, respectively. The junction condition for the Jordan-Brans-Dicke field φ is obtained from Eq(48), see Ref. [77], in particular we obtain: We note from this condition that in general V + and V − can not be both homogeneous.
By using the equations (49) and (50), we obtain As was mentioned, the exterior of the bubble (V + ) is described by the metric of a closed Friedmann-Robertson-Walker universe. We write this metric as follows: where r(χ + ) = sin(χ + ).
We can note from the conditions (51) that if one of the regions of space-time is homogeneous the other region is not. As in our case the outer region (V + ) is a homogeneous ES universe, then the region V − is generally inhomogeneous.
Substituting Eq. (44) into Eq. (56) and Eq. (57), we get The conditions of continuity for the metric on Σ are obtained from (43) and (58) as As usual, we consider that the content of matter in the outer region V + and in the inner region V − is a perfect fluid, with the following energy-momentum tensor where p, ρ and U µ are the pressure, energy density and 4-velocity of the perfect fluid, respectively. Then, by following Ref. [74], we can rewrite Eq.(59) and Eq.(60) as where So far we have summarized the calculations of Ref. [74] related to the evolution of a bubble in a Jordan-Brans-Dicke theory in a general background. Now we will concentrate on our particular case, where the bubble materializes and then evolves in a background corresponding to a closed and static universe, as the one described in Sec. II.
We assume that the outer region of the bubble is static, that is with a 0 and φ 0 constants, defined in Sec.II.
In the outer region, as content of matter we consider a scalar field ψ localized in a false vacuum, and a perfect fluid, see discussion in Sec. II. Then we have In the inner region, the matter content is described by the scalar field in the true vacuum (γ − = 0). Then, in the inner region we have Now, we develop the derivative of the left hand side of the equation (63) d Then the equation (63) takes the form where we have defined ε = (ρ + − ρ − )/σ. By using r(χ + ) = sin(χ + ) and R = a 0 r, we obtain If we introduce this into the Eq. (65), we have We are going to consider that the matter content of the bubble wall is the same as that of the outer region (V + ), then we haveγ = γ f . Therefore, we write Eq. (72) as follows Finally, Eq. (64) can be written as follows The evolution of the bubble wall is completely determined, in the outside coordinates, by the equations (74)- (76). We solved these equations numerically by consider different kind of matter content for the background, satisfying the stability conditions discussed in Sec. II. From these numerical solutions we found that once the bubble has materialized in the background of an ES universe in the context of a JBD theory, it grows filling the background space.
In particular in Figs The results for this case are shown in the Figure (5).
In the second example, we consider γ f = 1.3 and the following values for the radius and the JBD field The results for this case are shown in Figure (6).
In the third example, the matter content of the background is dust (γ f = 1) and we consider a 0 = 10 and φ 0 = 1.3. The results are shown in Figure (7).
In all these examples we have assumed the following initial conditions and units where 8πG = 1. From these examples we can note that the bubble of the new face, once materialized, grows to fill the background space without collapsing.

V. CONCLUSIONS
In this paper we study an alternative scheme for an Emergent Universe scenario called Emergent Universe by tunneling. In this scheme the universe is initially in a truly static state supported by a scalar field which is located in a false vacuum. The universe begins to evolve when, by quantum tunneling, the scalar field decays into a state of true vacuum.
The EU by tunneling scheme was originally developed in Ref. [40], in the context of General Relativity, where it was concluded that this mechanism is feasible as an EU scheme.
Nevertheless, this first model present the problem that the ES solution is classically instable.
The instability of the ES solution ensures that any perturbation, no matter how small, rapidly force the universe away from the static state, thereby aborting the EU scenario.
The present work is the natural extension of the idea presented in Ref. [40], but where the problem of the classical instability of the static solution is solved by going away from General Relativity and consider a JBD theory.
In particular, in this work we focus our study on the process of tunneling of a scalar field and the consequent creation and evolution of a bubble of true vacuum in the background of a classically stable Einstein Static universe. Our principal motivation is the study of new ways of leaving the static period and begin the inflationary regime in the context of Emergent Universe models.
In the first part of the paper, Sect. II, we study an Einstein static universe supported by a scalar field located in a false vacuum and its stability in the context of a JBD theory.
Contrary to General Relativity, we found that this static solution could be stable against isotropic perturbations if some general conditions are satisfied, see Eqs. (22)- (24). This modification of the stability behavior has important consequences for the emergent universe by tunneling scenario, since it ameliorates the fine-tuning that arises from the fact that the ES model is an unstable saddle in GR and it improves the preliminary model studied in Ref. [39]. In this study, for simplicity, we have not considered inhomogeneous or anisotropic perturbations. At this respect, the stability of the ES solution under anisotropic, tensor and inhomogeneous scalar perturbations have been studied in the context of JBD theories in Refs. [22,47]. It was found for theses JBD models that different from General Relativity [46] and others modified theories of gravity as f (R) [79] or modified Gauss-Bonnet gravity [80], that a static universe which is stable against homogeneous perturbations, could be also stable against anisotropic and inhomogeneous perturbations. We expect a similar behavior for our JBD model, were the static universe is supported by a scalar field located in a false vacuum. Then, we expect that in our case the inhomogeneous and anisotropic perturbations do not lead to additional instabilities. Nevertheless, we intend to return to these points in the near future by working an approach similar to that followed in Refs. [22,46,47,78].
In Sect. III we study the tunneling process of the scalar field from the false vacuum to the true vacuum and the consequent creation of a bubble of true vacuum in the background of Einstein static universe for a JBD theory. In particular we determinate the nucleation rate of the true vacuum bubble using the approaches developed in Ref. [70] and previous results obtained in Ref. [39].
The classical evolution of the bubble after its nucleation is studied in Sect. IV where we found that once the bubble has materialized in the background of an ES universe, it grows filling the background space. This demonstrates the viability of our EU model, since there is the possibility of having an open inflationary universe inside the bubble. During this study we consider the gravitational back-reaction of the bubble by using the formalism developed in Ref. [74] applied to a JBD theory. At this respect we found a system of coupled differential equations, which we solved numerically. Three specific examples of these solutions were shown in Sect. IV concerning to different background material contents.
It is worth to note that once the bubble has materialized, from conditions (51), it follows that if one of the regions of spacetime separated by the wall is homogeneous, then the other region is, in general, inhomogeneous [74]. Given that in our case the exterior of the bubble is a homogeneous universe, then the interior of the bubble will be, in general, inhomogeneous. However, since the degree of inhomogeneity depends on the difference in the energy density of the interior and the exterior of the bubble, it is possible in our case to decrease this inhomogeneity by adjusting the parameters of the static solution as was discussed in Ref. [74]. Then in our model, it is possible to study the feasibility of having an open inflationary universe inside the bubble. Nevertheless, given the similarities, we expect that the behavior inside the bubble of the EU by tunneling, will be similar to the models of single-field open and extended open inflation, as the ones studied in Refs. [58][59][60][61][62]. We expect to return to this point in the near future.