Vacuum bubbles may nucleate and expand during the inflationary epoch in the early universe. After inflation ends the bubbles quickly dissipate their kinetic energy they come to rest with respect to the Hubble flow and eventually form black holes. The fate of the bubble itself depends on the resulting black hole mass. If the mass is smaller than a certain critical value the bubble collapses to a singularity. Otherwise the bubble interior inflates forming a baby universe which is connected to the exterior FRW region by a wormhole. A similar black hole formation mechanism operates for spherical domain walls nucleating during inflation. As an illustrative example we studied the black hole mass spectrum in the domain wall scenario assuming that domain walls interact with matter only gravitationally. Our results indicate that depending on the model parameters black holes produced in this scenario can have significant astrophysical effects and can even serve as dark matter or as seeds for supermassive black holes. The mechanism of black hole formation described in this paper is very generic and has important implications for the global structure of the universe. Baby universes inside supercritical black holes inflate eternally and nucleate bubbles of all vacua allowed by the underlying particle physics. The resulting multiverse has a very nontrivial spacetime structure with a multitude of eternally inflating regions connected by wormholes. If a black hole population with the predicted mass spectrum is discovered it could be regarded as evidence for inflation and for the existence of a multiverse.
['High Energy Physics - Theory', 'Astrophysics - Cosmology and Nongalactic Astrophysics', 'General Relativity and Quantum Cosmology']
Black holes and the multiverse
{'Black holes and the multiverse': 'Jaume Garriga a,b , Alexander Vilenkin b and Jun Zhang b  \na Departament de Fisica Fonamental i Institut de Ciencies del Cosmos, Universitat de Barcelona, Marti i Franques, 1, 08028, Barcelona, Spain b Institute of Cosmology, Department of Physics and Astronomy, Tufts University, Medford, MA 02155, USA', 'Abstract': 'Vacuum bubbles may nucleate and expand during the inflationary epoch in the early universe. After inflation ends, the bubbles quickly dissipate their kinetic energy; they come to rest with respect to the Hubble flow and eventually form black holes. The fate of the bubble itself depends on the resulting black hole mass. If the mass is smaller than a certain critical value, the bubble collapses to a singularity. Otherwise, the bubble interior inflates, forming a baby universe, which is connected to the exterior FRW region by a wormhole. A similar black hole formation mechanism operates for spherical domain walls nucleating during inflation. As an illustrative example, we studied the black hole mass spectrum in the domain wall scenario, assuming that domain walls interact with matter only gravitationally. Our results indicate that, depending on the model parameters, black holes produced in this scenario can have significant astrophysical effects and can even serve as dark matter or as seeds for supermassive black holes. The mechanism of black hole formation described in this paper is very generic and has important implications for the global structure of the universe. Baby universes inside super-critical black holes inflate eternally and nucleate bubbles of all vacua allowed by the underlying particle physics. The resulting multiverse has a very non-trivial spacetime structure, with a multitude of eternally inflating regions connected by wormholes. If a black hole population with the predicted mass spectrum is discovered, it could be regarded as evidence for inflation and for the existence of a multiverse.', 'I. INTRODUCTION': "A remarkable aspect of inflationary cosmology is that it attributes the origin of galaxies and large-scale structure to small quantum fluctuations in the early universe [1]. The fluctuations remain small through most of the cosmic history and become nonlinear only in relatively recent times. Here, we will explore the possibility that non-perturbative quantum effects during inflation could also lead to the formation of structure on astrophysical scales. Specifically, we will show that spontaneous nucleation of vacuum bubbles and spherical domain walls during the inflationary epoch can result in the formation of black holes with a wide spectrum of masses.  \nThe physical mechanism responsible for these phenomena is easy to understand. The inflationary expansion of the universe is driven by a false vacuum of energy density ρ i ≈ const. Bubble nucleation in this vacuum can occur [2] if the underlying particle physics model includes vacuum states of a lower energy density, ρ b < ρ i . 1 We will be interested in the case when ρ b > 0. Once a bubble is formed, it immediately starts to expand. The difference in vacuum tension on the two sides of the bubble wall results in a force F = ρ i -ρ b per unit area of the wall, so the bubble expands with acceleration. This continues until the end of inflation (or until the energy density of the inflating vacuum drops below ρ b during the slow roll). At later times, the bubble continues to expand but it is slowed down by momentum transfer to the surrounding matter, while it is also being pulled inwards by the negative pressure of vacuum in its interior. Eventually, this leads to gravitational collapse and the formation of a black hole.  \nThe nature of the collapse and the fate of the bubble interior depends on the bubble size. The positive-energy vacuum inside the bubble can support inflation at the rate H b = (8 πGρ b / 3) 1 / 2 , where G is Newton's constant. If the bubble expands to a radius R > ∼ H -1 b , its interior begins to inflate 2 . We will show that the black hole that is eventually formed contains a ballooning inflating region in its interior, which is connected to the exterior region by a wormhole. On the other hand, if the maximum expansion radius is R glyph[lessmuch] H -1 b , then internal inflation does not occur, and the bubble interior shrinks and collapses to a singularity. 3  \nBubbles formed at earlier times during inflation expand to a larger size, so at the end of inflation we expect to have a wide spectrum of bubble sizes. They will form black holes with a wide spectrum of masses. We will show that black holes with masses above a certain critical mass have inflating universes inside.  \nThe situation with domain walls is similar to that with vacuum bubbles. It has been shown in Refs. [4] that spherical domain walls can spontaneously nucleate in the inflating false vacuum. The walls are then stretched by the expansion of the universe and form black holes when they come within the horizon. Furthermore, the gravitational field of domain walls is known to be repulsive [5, 6]. This causes the walls to inflate at the rate H σ = 2 πGσ , where σ is the wall tension. We will show that domain walls having size R > ∼ H -1 σ develop a wormhole structure, while smaller walls collapse to a singularity soon after they enter the cosmological horizon. Once again, there is a critical mass above which black holes contain inflating domain walls connected to the exterior space by a wormhole.  \nWe now briefly comment on the earlier work on this subject. Inflating universes contained inside of black holes have been discussed by a number of authors [7-9]. Refs. [8, 9] focused on black holes in asymptotically flat or de Sitter spacetime, while Ref. [7] considered a different mechanism of cosmological wormhole formation. The possibility of wormhole formation in cosmological spacetimes has also been discussed in Ref. [14], but without suggesting a cosmological scenario where it can be realized. Cosmological black hole formation by vacuum bubbles was qualitatively discussed in Ref. [15], but no attempt was made to determine the resulting black hole masses. Black holes formed by collapsing domain walls have been discussed in Refs. [16-18], but the possibility of wormhole formation has been overlooked in these papers, so their estimate of black hole masses applies only to the subcritical case (when no wormhole is formed).  \nIn the present paper, we shall investigate cosmological black hole formation by domain walls and vacuum bubbles that nucleated during the inflationary epoch. We shall study the spacetime structure of such black holes and estimate their masses. We shall also find the black hole mass distribution in the present universe and derive observational constraints on the particle physics model parameters. We shall see that for some parameter values the black holes produced in this way can serve as dark matter or as seeds for supermassive black  \nholes observed in galactic centers.  \nThe paper is organized as follows. Sections II and III are devoted, respectively, to the gravitational collapse of vacuum bubbles and domain walls. Section IV deals with the mass distribution of black holes produced by the collapse of domain walls, and Section V with the observational bounds on such distribution. Our conclusions are summarized in Section VI. A numerical study of test domain walls is deferred to Appendix A, and some technical aspects of the spacetime structure describing the gravitational collapse of large domain walls are discussed in Appendix B.", 'II. GRAVITATIONAL COLLAPSE OF BUBBLES': 'In this Section we describe the gravitational collapse of bubbles after inflation, when they are embedded in the matter distribution of a FRW universe. Three different time scales will be relevant for the dynamics. These are the cosmological scale  \nt ∼ H -1 = (3 / 8 πGρ m ) 1 / 2 , (1)  \nwhere ρ m is the matter density, the scale associated with the vacuum energy inside the bubble,  \nt b ≡ H -1 b = (3 / 8 πGρ b ) 1 / 2 , (2)  \nand the acceleration time-scale due to the repulsive gravitational field of the domain wall  \nt σ ≡ H -1 σ = (2 πGσ ) -1 . (3)  \nIn what follows we assume that the the inflaton transfers its energy to matter almost instantaneously, and for definiteness we shall also assume a separation of scales,  \nt i glyph[lessmuch] t b , t σ , (4)  \nwhere t i ∼ H -1 i is the Hubble radius at the time when inflation ends. The relation (4) guarantees that the repulsive gravitational force due to the vacuum energy and wall tension are subdominant effects at the end of inflation, and can only become important much later.  \nThe interaction of bubbles with matter is highly model dependent. For the purposes of illustration, here we shall assume that matter is created only outside the bubble, and has reflecting boundary conditions at the bubble wall. For the case of domain walls, which will  \nbe discussed in the following Section, we shall consider the case where matter is on both sides of the wall. As we shall see, this leads to a rather different dynamics.  \nThe dynamics of spherically symmetric vacuum bubbles has been studied in the thin wall limit by Berezin, Kuzmin and Tkachev [8]. 4 The metric inside the bubble is the de Sitter (dS) metric with radius H -1 b . The metric outside the bubble is the Schwarzschild-de Sitter metric 5 (SdS), with a mass parameter which can be expressed as  \nM = 4 3 π ( ρ b -ρ i ) R 3 +4 πσR 2 [ ˙ R 2 +1 -H 2 b R 2 ] 1 / 2 -8 π 2 Gσ 2 R 3 . (5)  \nHere, ρ i is the vacuum energy outside the bubble, R is the radius of the bubble wall, and ˙ R ≡ dR/dτ , where τ is the proper time on the bubble wall worldsheet, at fixed angular position. Eq. (5) can be interpreted as an energy conservation equation. 6', 'A. Initial conditions': "A bubble nucleating in de Sitter space has zero mass, M dS = 0, so the right hand side of (5) vanishes during inflation. Assuming that R , ˙ R and the bubble wall tension are continuous in the transition from inflation to the matter dominated regime, we have  \n4 3 πρ m ( t i ) R 3 i ≈ 4 3 πρ b R 3 i +4 πσR 2 i [ ˙ R 2 i +1 -H 2 b R 2 i ] 1 / 2 -8 π 2 Gσ 2 R 3 i . (6)  \nHere, we have assumed that the transfer of energy from the inflaton to matter takes place almost instantaneously at the time t i , so that ρ m ( t i ) ≈ ρ i , and that the vacuum energy ρ 0 outside the bubble is completely negligible after inflation,  \nρ 0 glyph[lessmuch] ρ m . (7)  \nNote that the left hand side of Eq. (6),  \nM i ≡ 4 3 πρ m ( t i ) R 3 i (8)  \nis the 'excluded' mass of matter in a cavity of radius R i at time t i , which has been replaced by a bubble of a much lower vacuum energy density ρ b . The difference in energy density goes into the kinetic energy of the bubble wall and its self-gravity corrections, corresponding to the second and third terms in Eq. (6).  \nLet us now calculate the Lorentz factor of the bubble walls with respect to the Hubble flow at the time t i . Assuming the separation of scales (4), the right hand side of Eq. (6) is dominated by the kinetic term and we have  \n˙ R i ≈ 1 4 H 2 i R i t σ . (9)  \nWe wish to consider the motion of the bubble relative to matter. The metric outside the bubble is given by  \nds 2 = -dt 2 + a 2 ( t )( dr 2 + r 2 d Ω 2 ) . (10)  \nHere, a ( t ) ∝ t β is the scale factor, where β = 1 / 2 for radiation, or β = 2 / 3 for pressureless matter, as may be the case if the energy is in the form of a scalar field oscillating around the minimum of a quadratic potential. The proper time on the worldsheet is given by  \ndτ = √ 1 -a 2 r ' 2 dt, (11)  \nwhere a prime indicates derivative with respect to t , and we have  \ndR dτ = HR + ar ' √ 1 -a 2 r ' 2 , (12)  \nwith H = a ' /a = β/t . The typical size of bubbles produced during inflation is at least of order t i , although it can also be much larger. Hence, at the end of inflation we have H i R i > ∼ 1. On the other hand, the second term in the numerator of (12) is bounded by unity, ar ' < 1. Combining (9) and (12) we have  \nγ i ≡ 1 √ 1 -a 2 i r ' 2 i ∼ t σ t i , (13)  \nand using (4), we find that the motion of the wall relative to matter is highly relativistic, γ i glyph[greatermuch] 1.", 'B. Dissipation of kinetic energy': "The kinetic energy of the bubble walls will be dissipated by momentum transfer to the surrounding matter. Let us now estimate the time-scale for this to happen. The force acting  \non the wall due to momentum transfer per unit area in the radial direction can be estimated as  \nT ˆ 0ˆ r ∼ ρ m U ˆ 0 U ˆ r ∼ -ρ m γ 2 , (14)  \nwhere U ˆ µ is the four-velocity of matter in an orthonormal frame in which the wall is at rest, hatted indices indicate tensor components in that frame, and we have used that the motion of matter is highly relativistic in the radial direction 7 . The corresponding proper acceleration of the wall is given by  \nα ∼ -ρ m γ 2 σ . (15)  \nIn the rest frame of ambient matter, we have  \ndγ dt = γ 3 vv ' = vα ≈ α, (16)  \nwhere v = ar ' is the radial outward velocity of the wall. The second equality just uses the kinematical relation between acceleration and proper acceleration for the case when velocity and acceleration are parallel to each other. Combining (15) and (16) we have  \ndγ γ 2 ∼ -t σ dt t 2 . (17)  \nSolving (17) with the initial condition (13) at t = t i , it is straightforward to show that the wall loses most of its energy on a time scale  \n∆ t ∼ t 3 i t 2 σ glyph[lessmuch] t i , (18)  \nand that the wall slows down to γ ∼ 1 in a time of order ∆ t ' ∼ t 2 i /t σ glyph[lessmuch] t i , which is still much less than the Hubble time at the end of inflation.  \nSince the bubble motion will push matter into the unperturbed FRW region, this results in a highly overdense shell surounding the bubble, exploding at a highly relativistic speed. The total energy of the shell is E shell ≈ M i , with most of this energy contained in a thin layer of mass M shell ∼ ( t i /G )( R i /t σ ) 2 expanding with a Lorentz factor γ shell ∼ γ 2 i ∼ ( t σ /t i ) 2  \nrelative to the Hubble flow. Possible astrophysical effects of this propagating shell are left as a subject for further research.  \nAfter the momentum has been transferred to matter, the bubble starts receding relative to the Hubble flow. As we shall see, for large bubbles with R i glyph[greatermuch] t i the relative shrinking speed becomes relativistic in a time-scale which is much smaller than the Hubble scale. To illustrate this point, let us first consider the case of pressureless matter.", 'C. Bubbles surrounded by dust': 'Once the bubble wall has stopped with respect to the Hubble flow, it is surrounded by an infinitesimal layer devoid of matter, with a negligible vacuum energy ρ 0 which we shall take to be zero. Since the matter outside of this layer cannot affect the motion of the wall, from that time on we can use Eq. (5) with ρ i replaced by the vanishing vacuum energy after inflation, ρ 0 = 0,  \nM bh = 4 3 πρ b R 3 +4 πσR 2 [ ˙ R 2 +1 -H 2 b R 2 ] 1 / 2 -8 π 2 Gσ 2 R 3 . (19)  \nThe mass parameter M bh for the Schwarzschild metric within the empty layer can be estimated from the initial condition ˙ R i ≈ H i R i , at the time t = t i +∆ t ≈ t i when the bubble is at rest with respect to the Hubble flow. Using the separation of scales (4), we have  \nM bh ≈ ( 4 3 πρ b +4 πσH i ) R 3 i . (20)  \nThe motion of the bubble is such that it eventually forms a black hole of mass M bh , nested in an empty cavity of co-moving radius  \nr c ≈ R i a ( t i ) . (21)  \nA remarkable feature of (20) is that M bh is proportional to the initial volume occupied by the bubble, with a universal proportionality factor which depends only on microphysical parameters. Under the condition (4), this factor is much smaller than the initial matter density ρ m ( t i ), and we have  \nM bh ∼ ( t 2 i t 2 b + t i t σ ) M i glyph[lessmuch] M i , (22)  \nwhere M i is the excluded mass given in (8). On dimensional grounds, the first term in parenthesis in (20) is of order η 4 b , where η b is the energy scale in the vacuum inside the  \nbubble, whereas the second term is of order η 3 w η 2 i /M p , where η w is the energy scale of the bubble wall and η i is the inflationary energy scale. The second term is Planck suppressed, and so in the absence of a large hierarchy between η w , η i and η b , the first term will be dominant. In this case we have M bh ≈ ( t i /t b ) 2 M i . However, we may also consider the limit when the vacuum energy inside the bubble is completely negligible, and then we have M bh ∼ ( t i /t σ ) M i .  \nIt was shown by Blau, Guendelman and Guth [9] that the equation of motion (19) can be written in the form  \n( dz d ˜ τ ) 2 + V ( z ) = E, (23)  \nwhere the rescaled variables z and ˜ τ are defined by  \nz 3 = H 2 + 2 GM bh R 3 , (24)  \n˜ τ = H 2 + 4 H σ τ, (25)  \nwith τ the proper time on the bubble wall trajectory, and H + defined by  \nH 2 + = H 2 b +4 H 2 σ . (26)  \nThe potential and the conserved energy in Eq. (23) are given by  \nV ( z ) = -( 1 -z 3 z 2 ) 2 -g 2 z , (27)  \nE = -16 H 2 σ (2 GM bh ) 2 / 3 H 8 / 3 + , (28)  \nwhere  \ng = 4 H σ H + . (29)  \nThe shape of the potential is plotted in Fig. 1. The maximum of the potential is at  \nz 3 = z 3 m ≡ 1 2 [ √ 8 + (1 -g 2 / 2) 2 -(1 -g 2 / 2) ] . (30)  \nNote that 0 ≤ g ≤ 2, and 1 ≤ z m ∼ 1 in the whole parameter range.  \nThe qualitative behaviour of the bubble motion depends on whether the Schwarzshild mass M bh is larger or smaller than a critical mass defined by  \nM cr = ¯ M g 3 z 6 m (1 -g 2 / 4) 1 / 2 3 √ 3( z 6 m -1) 3 / 2 , (31)  \nwhere  \n1  \nR  \n= 2  \nρ  \n=  \nρ  \n0  \nt  \nR  \n= 2  \nGM  \nρ  \n0  \nt  \n= 0  \nt  \n=  \nt  \nV  \n(  \nz  \n)  \nz  \nConsider a spherical domain  \ndifferent time-scales are relevant to  \nand the other is the acceleration t  \n1  \n-  \nσ  \nt  \n,  \nthe repulsive field  \nconformally stretched by cosmologi  \n.  \nConsider a spherical domain wall embedded in a spatially flat matter dominated unive  \ndifferent time-scales are relevant to the dynamics of such a wall. One is the cosmological scale  \n/lessmuch  \nand the other is the acceleration time-scale due to the repulsive gravitational field of the do  \nt  \nσ  \n∼  \n<!-- image -->  \nquickly shrinks under its tension a  \nHere, we will be primarily inter  \nσ  \nFor  \nR  \nt  \nσ  \n,  \nthe repulsive field can be ignored.  \nIn this case, for  \nt  \nR  \nt  \n,  \nthe doma  \nconformally stretched by cosmological expansion. Eventually, when the wall falls within the h  \n/lessmuch  \n/lessmuch  \n/lessmuch  \nthe matter around it while its size  \nquickly shrinks under its tension and forms a black hole of radius R S = 2 GM ∼ t 2 /t σ /lessmuch t . Here, we will be primarily interested in the opposite limit, where R /greatermuch t σ . In this case, the the matter around it while its size grows faster than the ambient expansion rate. As we shal leads to the formation of a wormho a baby universe, and in the ambie FIG. 1: The potential in Eq. (27), for g = 0 (upper curve), g = 1 (middle curve) and g = 2 (lower curve).  \nleads to the formation of a wormhole. The dust which was originally in the interior of the wall  \nspherical region of of vacuum.  \na baby universe, and in the ambient FRW universe we are left with a black hole remnant c  \nspherical region of of vacuum.  \nand a dust cosmology. 1.1 Matching Schwarzschild to a dust cosmology 1.1 Matching Schwarzschil Although the expression for the critical mass is somewhat cumbersome, we can estimate it as  \nBefore we consider the effect of  \nBefore we consider the effect of the domain wall, let us first discuss the matching of Sch and a dust cosmology. ¯ M ≡ t b 2 G . (32)  \nConsider the Schwarzschild metric Consider the Schwarzschild metric GM cr ∼ Min { t σ , t b } . (33)  \nds 2 = -( 1 -2 GM R ) dT 2 + ( 1 -2 GM R ) -1 dR 2 + R 2 d Ω 2 . In Lemaitre coordinates, this takes the form ds 2 = -( 1 In Lemaitre coordinates, this takes On dimensional grounds t b ∼ M p /η 2 b , while t σ ∼ M 2 p /η 3 σ . In the absence of a hierarchy between η b and η σ we may expect GM cr ∼ t b , but the situation where GM cr ∼ t σ is of course also possible.  \nds  \n2  \n=  \n-  \nare defined by the relations', 'where τ and ρ 1. Small bubbles surrounded by dust': 'dτ = ± dT -√ 2 GM R ( 1 -2 GM R ) -1 dR, dρ = ∓ dT + √ R 2 GM ( 1 -2 GM R ) -1 dR. Subtracting (3) from (4), we have 3 2 / 3 1 / 3 where τ and ρ are defined by the r dτ dρ For M bh < M cr , we have V ( z m ) > E , and the bubble trajectory has a turning point. In this case it is easy to show that with the condition (4), the bubble trajectory is confined to the left of the potential barrier, z i < 1 ≤ z m , where z i is the value of z which corresponds to the initial radius R i . The bubbles grow from small size up to a maximum radius R max and then recollapse. The turning point is determined from (19) with ˙ R = 0. For M bh glyph[lessmuch] M cr we have  \nwhere an integration constant has been absorbed by a shift in the origin of the  \nρ  \ncoordin  \nexpression for  \nT  \nas a function of  \nτ  \nand  \nρ  \ncan be found from (3) as:  \n±  \nT  \n= 4  \nGM  \n1  \n3  \n(  \nR  \n)  \n3  \n/  \n2  \nR = [ 2 ( τ + ρ ) ] (2 GM ) , M bh ≈ 4 3 πρ b R 3 max +4 πσR 2 max . (34)  \nSubtracting (3) from (4), we have  \n√  \nR  \n-  \n√  \nR  \ndτ  \n2  \n+  \ndρ  \n2  \n+  \nR  \n2  \nd  \nΩ  \n2  \n.  \n2  \nGM  \nR  \n2  \nGM ρ  \n=  \nconst. r  \n= 0  \nr  \n2  \nGM  \n+  \n2  \nGM  \nS  \n-  \nR  \n= 0  \nR  \n= 2  \nGM  \ni  \nH  \n-  \n1  \n= 2  \nGM  \nbh  \nbh  \ntanh  \n1  \n2  \nGM  \n-  \nρ,  \n,T  \nbh  \nn  \na  \ni  \nn  \n,T  \nw  \na  \nt  \nl  \nh  \ne  \nw  \na  \nl  \nl  \nR  \ne  \n→∞  \nl  \ne  \n-  \ns  \ne  \nc  \ne  \nd  \na  \ni  \ns  \ny  \nl  \ng  \ni  \nm  \nn  \ne  \nn  \no  \nl  \na  \nb  \na  \nd  \no  \nn  \nr  \n.  \nc  \nk  \no  \ns  \nd  \nd  \ni  \nm  \nu  \ne  \ne  \nc  \ne  \nd  \nE  \nh  \ni  \nt  \nt  \nw  \n.  \nv  \no  \ne  \nh  \ne  \nl  \nI  \nn  \ne  \nl  \ni  \ne  \na  \ns  \nr  \ne  \ns  \ns  \nt  \no  \no  \n1  \ni  \nn  \ny  \ng  \n=  \no  \n2  \nl  \no  \nm  \nGM  \ns  \nco  \nd  \nt  \ndus  \nn  \ni  \nf  \nt  \nn  \nt  \no  \ni  \nm  \nm  \na  \no  \nr  \nl  \ni  \ne  \nt  \nt  \n,  \nb  \ng  \nf  \ni  \ni  \nt  \nd  \n2  \n.  \ni  \nn  \nw  \ns  \nc  \na  \ni  \nt  \nu  \ns  \ns  \nt  \n0  \nr  \n=  \nr  \n0  \n=  \nr  \n.  \nst  \ncon  \n=  \nρ  \nGM  \n2  \n-  \ndus  \n=  \nρ  \n∞  \nbh  \n→-  \n,T  \nbh  \n,T  \na  \ndo  \nGM  \nca  \nl  \n2  \n=  \nR  \nm  \n0  \n→-  \nt  \na  \ni  \na  \nn  \ns  \ni  \nr  \nu  \nc  \nSphe  \nl  \ns  \nh  \nl  \nh  \na  \ne  \nw  \nr  \nR  \np  \nR  \nS  \n=  \n→∞  \na  \n=  \ne  \nn  \np  \ni  \na  \n2  \nGM  \ni  \nt  \n=  \n0  \nn  \na  \nl  \nt  \ni  \nr  \ni  \nw  \n,T  \nbh  \na  \n=  \n0  \nC  \nl  \ns  \ns  \nl  \nu  \no  \ne  \nd  \nt  \nl  \nt  \nbh  \nr  \ni  \ne  \ne  \nGM  \n2  \n.  \nGM  \nO  \nv  \nff  \nn  \nh  \nt  \na  \nt  \ne  \n1  \n)  \no  \nz  \ni  \nr  \no  \nh  \ne  \nh  \nt  \nn  \ni  \nh  \nt  \ni  \nw  \ns  \nl  \nl  \nl  \na  \nσ  \n/t  \n2  \n.  \nt  \ne  \nr  \nl  \nl  \na  \nw  \ne  \nh  \nt  \n,  \ne  \ns  \n/lessmuch  \na  \nc  \n,  \ne  \ne  \ns  \nl  \nl  \na  \nh  \ns  \ne  \nw  \ns  \nA  \n.  \ne  \nt  \na  \nr  \nt  \no  \n.  \ne  \ns  \nr  \ne  \nv  \ni  \nn  \nu  \nd  \ne  \nt  \na  \nn  \ni  \nm  \no  \ns  \ne  \no  \ng  \nl  \nl  \na  \nw  \ne  \nh  \nt  \nf  \no  \nr  \no  \ni  \nr  \ne  \n=  \nt  \nρ  \no  \ns  \nz  \nr  \na  \nw  \nh  \nc  \nS  \nf  \no  \nρ  \ng  \nn  \ni  \nh  \nGM  \nc  \nt  \na  \nm  \n=  \nr  \ncon  \nst  \n.  \ni  \n=  \nd  \nu  \nn  \nr  \ni  \nv  \ne  \ns  \ne  \nt  \ne  \n∼  \no  \nc  \na  \nc  \nl  \ne  \nn  \nr  \nt  \ng  \ni  \na  \no  \nl  \no  \nl  \no  \ne  \nm  \ns  \nm  \nd  \ne  \nh  \na  \nt  \nt  \nm  \nff  \no  \no  \na  \nc  \nn  \nl  \ns  \na  \ne  \nh  \nt  \nd  \nm  \na  \nw  \ni  \nn  \nn  \nw  \ni  \ni  \na  \ns  \nl  \nl  \nt  \nn  \ni  \nh  \nm  \na  \no  \ne  \nh  \n.  \nt  \nr  \nh  \no  \ne  \nd  \nh  \n.  \n2  \nt  \n,  \nΩ  \n/lessmuch  \nm  \nd  \nt  \nσ  \nt  \ni  \ni  \nt  \n∼  \nt  \ne  \nl  \na  \nc  \ns  \nl  \na  \nc  \ni  \ng  \no  \nl  \no  \nm  \ns  \n/lessmuch  \no  \ng  \nt  \ny  \ns  \ni  \nh  \nt  \n∼  \nn  \nd  \ne  \nt  \ns  \ny  \nc  \nt  \nn  \na  \nn  \nm  \ne  \nr  \ne  \n2  \nl  \no  \n-  \nh  \n<!-- image -->  \n-  \no  \nσ  \nt  \nt  \n.  \nm  \no  \nd  \nR  \nσ  \n,  \nv  \ne  \nl  \ne  \na  \nh  \nr  \nc  \ne  \nt  \nl  \ne  \ne  \nt  \nu  \nw  \nn  \ni  \na  \nm  \no  \nd  \ne  \nh  \nt  \n,  \nσ  \nt  \nm  \nv  \ni  \nr  \ns  \nGM  \ne  \nl  \nd  \nr  \ni  \no  \ne  \nw  \ns  \nt  \nbh  \ne  \ne  \nr  \na  \ne  \nz  \no  \nw  \nl  \ne  \nl  \n,  \ni  \nl  \ne  \np  \n,  \ns  \na  \nl  \nl  \nl  \nl  \ne  \ns  \ni  \ns  \nt  \no  \nt  \nh  \nn  \ne  \ns  \ng  \na  \ns  \nh  \nw  \nn  \na  \ne  \nh  \n/lessmuch  \nt  \n,  \nw  \nt  \ne  \ns  \nf  \no  \n.  \nr  \n2  \nr  \nΩ  \ne  \nd  \no  \nl  \n2  \nh  \nR  \n+  \n2  \ne  \ne  \no  \nt  \nc  \n1  \ns  \ni  \na  \ni  \ne  \nd  \n,  \ndR  \nh  \nc  \n-  \nn  \nz  \nw  \nS  \n)  \no  \nn  \nf  \ng  \nGM  \n2  \n-  \n1  \n(  \nGM  \n2  \n√  \na  \nr  \ni  \ns  \nc  \nh  \n.  \ndR  \n1  \n-  \n)  \nGM  \n2  \n-  \n1  \n(  \nR  \n2  \nGM  \n2  \n√  \n-  \n2  \ndT  \nΩ  \nR  \n±  \n+  \n=  \ndR  \nm  \ny  \nt  \nn  \ni  \nt  \nh  \ns  \nn  \no  \ni  \nt  \na  \n1  \nl  \ne  \n-  \nr  \ne  \nh  \nt  \nh  \ne  \nc  \nm  \nρ  \nd  \nGM  \n2  \n+  \n2  \nτ  \nd  \n-  \n=  \n2  \nds  \nGM  \n)  \nt  \n→-∞  \n2  \n-  \ns  \nn  \nm  \ni  \nm  \ns  \no  \nc  \np  \nn  \nr  \nr  \nw  \ni  \nR  \ny  \nb  \n2  \nGM ρ  \n=  \nconst. r  \n= 0  \nr  \n=  \nr  \nτ  \nd  \n2  \nh  \n2  \n.  \n+  \ndT  \n∓  \n=  \nρ  \nd  \n.  \n.  \ne  \nt  \na  \nn  \ni  \nd  \nr  \no  \no  \nc  \nρ  \ne  \nh  \nt  \nf  \no  \nn  \ni  \ng  \ni  \nr  \no  \ne  \nh  \nt  \nn  \ni  \nt  \nf  \ni  \nh  \ns  \na  \ny  \nb  \nd  \ne  \nb  \nr  \no  \ns  \nb  \na  \nn  \ne  \ne  \nb  \ns  \na  \nh  \nt  \nn  \na  \nt  \ns  \nConsider a spherical domain wall embedded in a spatially flat matter dominated universe.  \nu  \np  \nm  \nw  \na  \nf  \ne  \no  \ne  \nl  \ny  \nl  \ni  \nr  \na  \na  \nm  \nl  \ni  \nh  \nt  \na  \nb  \ne  \nh  \nt  \nn  \n.  \nf  \nt  \ne  \ne  \nh  \no  \nt  \nm  \nu  \nu  \nc  \na  \nc  \nt  \nd  \ne  \nt  \ns  \ndifferent time-scales are relevant to the dynamics of such a wall. One is the cosmological scale  \nt  \nh  \nm  \nn  \nu  \ni  \nv  \nff  \ne  \nf  \ne  \ni  \no  \nzs  \nr  \na  \nw  \nh  \na  \nt  \ndus  \nl  \ni  \nh  \nc  \ns  \nz  \nd  \ne  \nn  \nfi  \ne  \nn  \ne  \ni  \no  \n+  \nΩ  \nR  \n1  \nw i t h R S i n c e t h e 2 GM R Subtracting (3) from (4), we have ∆ R R i ≈ σ ρ b H i glyph[lessmuch] 1 , (35)  \n-  \ne  \nv  \na  \nh  \ne  \nw  \n,  \n)  \n4  \n(  \n1  \nA  \nn  \no  \nc  \ncon st . r = 0 r = r 0 ρ = t e r d o m i n a t e d u n i v e r s e . T w o i s t h e c o s m o l o g i c a l s c a l e t ∼ H -1 , g r a v i t a t i o n a l fi e l d o f t h e d o m a i n w a l l c a s e , f o r t /lessmuch R /lessmuch t σ , t h e d o m a i n w a l l i s u a l l y , w h e n t h e w a l l f a l l s w i t h i n t h e h o r i z o n , i t f r a d i u s R S = 2 GM ∼ t 2 /t σ /lessmuch t . w h e r e R /greatermuch t σ . I n t h i s c a s e , t h e w a l l r e p e l s e n t e x p a n s i o n r a t e . A s w e s h a l l s e e , t h i s l l y i n t h e i n t e r i o r o f t h e w a l l g o e s i n t o h a b l a c k h o l e r e m n a n t c y s t e d i n a t h e m a t c h i n g o f S c h w a r z s c h i l d a n t σ ∼ ( F o r c o n f o r m a l l q u i c k l y s h r i n k H e r e , w e w i l l t h e m a t t e r a r o u n d l e a d s t o t h e f o r m a t i o n a b a b y u n i v e r s e , a n d i n s p h e r i c a l r e g i o n o f o f v a c B e f o r e w e c o n s i d e r t h e a n d a d u s t c o s m o l o g y . 1 . 1 M a t c hi ng Sc h w a r zs c hi l d C o n s i d e r t h e S c h w a r z s c h i l d m e t r i c ds 2 = -( 1 -2 GM R ) dT I n L e m a i t r e c o o r d i n a t e s , t h i s t a k e s t h e f o r m ds 2 = -d τ 2 + 2 GM R d ρ w h e r e τ a n d ρ a r e d e fi n e d b y t h e r e l a t i o n s d τ = ± dT -√ 2 GM R ( 1 -2 GM R ) -d ρ = ∓ dT + √ R 2 GM ( 1 -2 GM R ) -1 dR . S u b t r a c t i n g ( 3 ) f r o m ( 4 ) , w e h a v e R = [ 3 2 ( τ + ρ ) ] 2 / 3 ( 2 GM ) 1 / 3 w h e r a n i n t e g r a t i o n c o n s t a n t h a s b e e n a b s o r b e d b y a s h i e x p r e s s f o r T a s a f u n c t i o n o f τ a n d ρ c a n b e f o u n d f r ± T = 4 GM   1 3 ( R 2 GM ) 3 / 2 + g i v e n i n ( 5 ) . s e c o n d i n t e g r a t i o n m e t r i c ( 2 ) i s c h r o n o u s , t h e and the other is the acceleration time-scale due to the repulsive gravitational field of the domai t σ ∼ ( Gσ ) -1 . For R /lessmuch t σ , the repulsive field can be ignored. In this case, for t /lessmuch R /lessmuch t σ , the domain conformally stretched by cosmological expansion. Eventually, when the wall falls within the horiz quickly shrinks under its tension and forms a black hole of radius R S = 2 GM ∼ t 2 /t σ /lessmuch t . Here, we will be primarily interested in the opposite limit, where R /greatermuch t σ . In this case, the wall the matter around it while its size grows faster than the ambient expansion rate. As we shall se leads to the formation of a wormhole. The dust which was originally in the interior of the wall goe a baby universe, and in the ambient FRW universe we are left with a black hole remnant cyste spherical region of of vacuum. Before we consider the effect of the domain wall, let us first discuss the matching of Schwarz and a dust cosmology. 1.1 Matching Schwarzschild to a dust cosmology Consider the Schwarzschild metric ds 2 = -( 1 -2 GM R ) dT 2 + ( 1 -2 GM R ) -1 dR 2 + R 2 d Ω 2 . In Lemaitre coordinates, this takes the form ds 2 = -dτ 2 + 2 GM R dρ 2 + R 2 d Ω 2 . where τ and ρ are defined by the relations dτ = ± dT -√ 2 GM R ( 1 -2 GM R ) -dR, dρ = ∓ dT + √ R ( 1 -2 GM ) -1 dR. 1 Sphe r i ca l do R = 2 GM S R = 0 R = ρ 0 t = 0 t = t i H -1 = 2 V ( z ) z C o n s i d e r a s p h e r i c a l d i ff e r e n t t i m e -s c a l e s a r e a n d t h e o t h e r i s t h e a c c t σ ∼ ( G σ ) -1 . F o r R /lessmuch t σ , t h e r e c o n f o r m a l l y s t r e t c h e d b q u i c k l y s h r i n k s u n d e r i t H e r e , w e w i l l b e p r i t h e m a t t e r a r o u n d i t w l e a d s t o t h e f o r m a t i o a b a b y u n i v e r s e , a n d s p h e r i c a l r e g i o n o f o B e f o r e w e c o n s i d a n d a d u s t c o s m o l o g 1 . 1 M a t c hi ng Sc C o n s i d e r t h e S c h w a r I n L e m a i t r e c o o r d i n a w h e r e τ a n d ρ a r e d S u b t r a c t i n g ( 3 ) f r o m w h e r e a n i n t e g r a t i o n FIG. 2: Causal diagram showing the formation of a black hole by a small vacuum bubble, with mass M bh glyph[lessmuch] M cr , in a dust dominated spatially flat FRW. At the time t i when inflation ends, a small bubble (to the left of the diagram) with positive energy density ρ b > 0, initially expands with a large Lorentz factor relative to the Hubble flow. The motion is slowed down by momentum transfer to matter, turning around due to the internal negative pressure and the tension of the bubble wall, and eventually collapsing into a Schwarzschild singularity. At fixed angular coordinates, the timelike dashed line at the edge of the dust region represents a radial geodesic of the Schwarzschild metric (see Appendix B). Such geodesic approaches time-like infinity with unit slope. The reason is that a smaller slope would correspond to the trajectory of a particle at a finite radial coordinate R , while the geodesic of a particle escaping to infinity has unbounded R . The thin relativistic expanding shell of matter produced by the bubble is not represented in the figure. This shell disturbs the homogeneity of dust near the boundary represented by the dashed line. For σH i /ρ b glyph[lessmuch] 1, the first term in (20) dominates, and the turning point occurs after a small fractional increase in the bubble radius  \n2  \n/  \n3  \n± T = 4 GM   1 3 ( R 2 GM ) 3 / 2 + √ R 2 GM -tanh -1 √ R 2 GM   R max R i ≈ ( 3 σH i ρ b ) 1 / 3 . (36)  \nwith  \nR  \ngiven in (5). A second integration constant has been absorbed by a shift in the  \nT  \ncoord  \nSince the metric (2) is synchronous, the lines of constant spatial coordinate are geodesics, and  \nτ  \nR = [ 3 2 ( τ + ρ ) ] (2 GM ) 1 / 3 , where an integration constant has been absorbed by a shift in the origin of the expression for T as a function of τ and ρ can be found from (3) as: which happens after a time-scale much shorter than the expansion time, ∆ t ∼ ( σ/ρ b ) glyph[lessmuch] t i . For 1 glyph[lessmuch] σH i /ρ b glyph[lessmuch] ( R i H i ) 3 / 2 we have  \n1  \n)  \ndR  \n,  \nn  \ns  \ny  \nρ  \ncoordinate.  \no  \ns  \ni  \nh  \nt  \n,  \ns  \ne  \n2  \n1  \n+  \n(  \nt  \n2  \nt  \ni  \nr  \ne  \ne  \no  \n=  \nf  \nρ  \nu  \nr  \ne  \ny  \nt  \nw  \n=  \n2  \nds  \n∼  \ns  \nt  \nR  \nρ  \n=  \n-  \nd  \nR  \nR  \nR  \na  \nh  \nb  \na  \ns  \ny  \n/lessmuch  \ne  \nn  \nG  \nm  \nn  \ni  \na  \nm  \no  \nd  \ne  \nh  \nt  \nf  \no  \nd  \nl  \ne  \nfi  \nSchwarzschild"singularity" s d c o n a l R l l f . I  \nn  \nh  \nl  \nfi  \nf  \nd  \ne  \nd  \nd  \n(  \ni  \nt  \nρ  \nm  \ni  \ns  \nh  \n(  \n2  \n)  \nt  \n.  \ny  \n2  \nd  \nl  \nρ,  \n3  \n/  \n2  \n]  \n)  \nρ  \n+  \nτ  \n(  \n3 2  \n[  \n=  \nR  \n,  \n3  \n/  \n1  \n)  \n1  \nGM  \n(  \n2  \n(  \n)  \nFinally, for σH i /ρ b glyph[greatermuch] ( R i H i ) 3 / 2 we have  \nR max R i ≈ ( H i R i ) 1 / 2 . (37)  \nIn the last two cases, R max glyph[greatermuch] R i . From Eq. (19), we find that for R i < ∼ R glyph[lessmuch] R max the bubble radius behaves as R ∝ τ 1 / 3 . This is in contrast with the behaviour of the scale factor in the matter dominated universe, a ( t ) ∝ t 2 / 3 , and therefore the bubble wall decouples from the Hubble flow in a time-scale which is at most of order t i . Once the bubble radius reaches its turning point at R = R max , the bubble collapses into a Schwarzschild singularity. This is represented in Fig. 2. Denoting by t H the time it takes for the empty cavity of co-moving radius r c = R i /a ( t i ) to cross the horizon, we have  \nt H ∼ ( H i R i ) 2 R i . (38)  \nNote that t H glyph[greatermuch] ( H i R i ) 1 / 2 R i > ∼ R max , so the black hole forms on a time-scale much shorter than t H .', '2. Large bubbles surrounded by dust': "For M bh > M cr we have V ( z m ) > E , so the trajectory is unbounded and the bubble wall grows monotonically towards infinite size in the asymptotic future. The unbounded growth of the bubble is exponential, and starts at the time t ∼ Min { t b , t σ } . This leads to the formation of a wormhole which eventually 'pinches off', leading to a baby universe (see Figs. 3 and 4). The bubble inflates at the rate H b , and transitions to the higher energy inflatinary vacuum with expansion rate H i will eventually occur, leading to a multiverse structure [3, 12]. Initially, a geodesic observer at the edge of the matter dominated region can send signals that will end up in the baby universe. However, after a time t ∼ GM bh , the wormhole closes and any signals which are sent radially inwards end up at the Schwarzschild singularity.", 'D. Bubbles surrounded by radiation': "Due to the effects of pressure, the dynamics of bubbles surrounded by radiation is somewhat more involved than in the case of pressureless dust. Here, we will only present a  \na  \ne  \nn  \ns  \nd  \nt  \nl  \na  \np  \nh  \ne  \nh  \na  \nb  \ne  \nB  \na  \nM  \nr  \nt  \nd  \nc  \nq  \ne  \nd  \na  \nr  \ns  \nb  \ni  \ne  \no  \nu  \nt  \nn  \ni  \na  \nσ  \nF  \nd  \ni  \nn  \nρ  \nff  \nd  \n0  \nt  \nC  \nr  \ne  \nt  \n(  \n∼  \no  \nf  \nc  \nH  \nk  \no  \nr  \nl  \ne  \na  \nm  \nt  \ny  \nc  \nf  \nd  \na  \nh  \ni  \no  \nu  \nt  \ne  \nn  \na  \nb  \na  \no  \nu  \nl  \nr  \ns  \nc  \ne  \nt  \ny  \nr  \nt  \ne  \nr  \nm  \ns  \n,  \nt  \nt  \nn  \nr  \ne  \nh  \ni  \ne  \nw  \nc  \nhi  \nS  \nt  \ny  \nc  \ne  \ns  \nt  \n,  \nh  \nr  \ne  \nv  \ng  \ne  \nh  \nG  \nR  \nl  \na  \nh  \nw  \nl  \nr  \n=  \n=  \no  \nn  \ne  \nσ  \nn  \nt  \no  \n)  \n2  \n0  \ns  \nGM  \ni  \nt  \nt  \nd  \ni  \nt  \nh  \n1  \n-  \ns  \n/lessmuch  \ny  \ni  \ne  \na  \nf  \ne  \ni  \ne  \no  \ns  \nr  \no  \nc  \nm  \nng  \nh  \nw  \na  \nr  \nr  \no  \ns  \nn  \no  \no  \nn  \nl  \nn  \nw  \no  \nr  \ne  \n,  \no  \nk  \nt  \ns  \ni  \nl  \nu  \nm  \nl  \nn  \na  \na  \ne  \nm  \ne  \n.  \nt  \nr  \nσ  \ne  \nu  \nb  \nd  \nt  \ni  \nn  \no  \nf  \ns  \no  \nSc  \nz  \ns  \nc  \ni  \ng  \nh  \nh  \n2  \nd  \ny  \nf  \ne  \nr  \n.  \nw  \ni  \nl  \n=  \nt  \na  \nd  \nr  \n=  \nr  \ne  \nr  \n,  \nt  \nn  \ne  \ni  \no  \nd  \nv  \nt  \nn  \ni  \na  \nt  \nn  \nc  \nh  \nzs  \nm  \n-  \na  \nk  \ne  \ns  \n2  \nds  \nt  \nds  \nt  \ne  \n=  \nh  \nr  \ne  \nl  \n±  \ndT  \nρ  \n)  \na  \ni  \ndT  \n+  \n2  \no  \nn  \n=  \ns  \n-  \n√  \n2  \n/  \n√  \n2  \nR  \nR  \n(  \n1  \n1  \n-  \n/  \n,  \n3  \ne  \n)  \nt  \nn  \nGM  \n3  \n(  \n]  \na  \ns  \nm  \n2  \nh  \ni  \nGM  \nf  \n(  \nn  \nt  \n3  \nh  \ne  \nd  \nn  \nh  \ni  \na  \n)  \n1  \n-  \n2  \ns  \n:  \n√  \nb  \na  \ny  \nt  \ne  \nGM  \na  \na  \ns  \nr  \nR  \n-  \nh  \ne  \ni  \nf  \ng  \nt  \ne  \ni  \no  \nn  \nd  \ne  \ni  \ng  \ni  \nn  \n,  \nρ  \nt  \ns  \nh  \ni  \ne  \nc  \ns  \n,  \nT  \na  \no  \nf  \nc  \nd  \n-  \nS  \nt  \nR  \ni  \na  \ns  \ni  \nt  \nh  \nc  \ns  \nh  \ne  \ne  \np  \nr  \nr  \ns  \na  \nH  \np  \nl  \nt  \ne  \nd  \ni  \ni  \nw  \no  \nh  \nf  \nt  \nu  \ne  \nc  \nr  \n-  \nGM  \nt  \nρ  \n=  \nt  \nm  \na  \n.  \nl  \ndo  \nl  \nca  \nr  \n0  \n=  \nh  \ns  \ne  \nh  \nr  \ne  \ne  \nb  \nt  \nm  \ni  \ns  \na  \nl  \na  \nh  \nu  \ne  \nhi  \ni  \nc  \nt  \nc  \nn  \nh  \no  \nd  \ne  \no  \nr  \nτ  \ne  \nm  \nff  \nl  \n-  \ne  \n1  \nr  \ni  \na  \na  \np  \ny  \nt  \nr  \ne  \nw  \ni  \nR  \n=  \nc  \nr  \nc  \nu  \nl  \nc  \ne  \ni  \nl  \ni  \nc  \ns  \no  \nn  \ny  \nt  \no  \na  \n.  \ne  \nd  \n=  \nm  \nGM  \na  \ng  \no  \ns  \no  \nm  \np  \nρ  \n.  \ncon  \ne  \nn  \nd  \na  \nm  \nr  \nt  \ni  \no  \nr  \nt  \nd  \ne  \na  \ns  \nt  \nv  \nu  \ne  \nn  \ng  \no  \ns  \nc  \no  \nm  \nc  \nh  \no  \ne  \na  \ni  \nt  \na  \nn  \nt  \nr  \ni  \ne  \nc  \na  \nl  \ni  \nl  \nt  \nl  \ns  \na  \na  \no  \ne  \nt  \nn  \nm  \nd  \nh  \n,  \nw  \nt  \nt  \nσ  \ni  \na  \no  \ne  \nn  \nm  \nw  \ni  \ni  \na  \ns  \nl  \nl  \nt  \nn  \ni  \nh  \nh  \no  \ne  \n.  \nt  \nt  \ni  \ni  \ns  \ne  \nl  \no  \nh  \nσ  \n/lessmuch  \nt  \n,  \na  \ns  \nr  \nh  \ne  \nc  \nn  \nt  \n/t  \ni  \ns  \n.  \na  \nt  \ne  \nr  \na  \ne  \n,  \nz  \no  \nw  \nl  \ne  \nl  \ni  \nl  \ne  \nl  \nh  \ne  \ns  \nl  \ns  \nw  \nA  \ne  \ni  \nt  \ne  \ni  \nf  \nr  \nr  \nn  \no  \np  \ns  \ni  \ns  \nt  \nh  \n,  \ns  \na  \ne  \nn  \ne  \ns  \ng  \na  \ne  \nw  \nh  \no  \nt  \no  \nl  \nl  \nt  \nc  \nn  \nm  \na  \ni  \nt  \no  \ni  \nt  \nn  \na  \ni  \ne  \ny  \nt  \nn  \ne  \ne  \nh  \nc  \nb  \no  \nk  \nl  \na  \nu  \ni  \nt  \ns  \ns  \no  \nn  \nc  \nm  \nh  \ns  \nc  \ny  \nd  \ng  \no  \nm  \nFIG. 4: A space-like slice of an inflating bubble connected by a wormhole to a dust dominated flat FRW universe.  \n<!-- image -->  \ni  \ns  \ns  \ni  \ni  \nv  \nm  \no  \nn  \ns  \no  \nn  \nt  \ni  \ne  \nz  \nh  \ns  \nr  \nm  \no  \nc  \nt  \nt  \no  \ndR  \nρ  \nd  \ni  \ni  \ns  \n.  \nc  \nn  \no  \na  \nt  \nt  \nh  \ne  \n.  \ne  \no  \nm  \nb  \nf  \ni  \ne  \nt  \na  \nr  \ne  \no  \nl  \nn  \nh  \ne  \ns  \ng  \ne  \nt  \ne  \nt  \nr  \n.  \nF  \nd  \no  \nT  \nR  \no  \nw  \nh  \nW  \nm  \na  \ne  \ni  \nR  \nl  \no  \nl  \ns  \ny  \ni  \no  \nm  \nn  \ns  \nc  \nn  \nb  \na  \ns  \nfi  \nt  \n2  \n=  \n∞  \nt  \nr  \ni  \na  \nl  \na  \nw  \nt  \n,  \n=  \n0  \nb  \ni  \ne  \n0  \na  \nn  \nl  \ne  \nn  \n0  \ni  \nbh  \na  \n-  \ne  \ne  \nt  \ni  \nn  \nl  \no  \ng  \ni  \nh  \na  \nConsider a spherical domain wall embedded in a spatially flat matter dominated universe.  \nTwo  \nr  \ni  \ne  \ny  \na  \ne  \nl  \nd  \ni  \ni  \ns  \nt  \nl  \ni  \nz  \ns  \ne  \ne  \nd  \na  \no  \nl  \nl  \nm  \nw  \no  \nu  \ne  \nh  \nr  \ns  \nb  \na  \nt  \ni  \ne  \ni  \nm  \nF  \nn  \nv  \ne  \nr  \nt  \nd  \nh  \n=  \nf  \no  \ne  \n0  \nr  \nw  \nt  \nt  \nn  \nc  \ni  \ns  \nu  \ns  \nr  \nt  \ni  \nh  \ni  \nd  \n(  \ne  \n2  \n)  \nt  \nr  \nw  \na  \ni  \nn  \n,T  \nbh  \n=  \n2  \nl  \n2  \nGM  \nd  \nr  \n1  \ns  \ne  \na  \ne  \n2  \nw  \na  \ni  \ne  \nn  \ne  \nt  \n1  \nGM  \nh  \nw  \n(  \ne  \n4  \n)  \n=  \ns  \n(  \n3  \n=  \n)  \n<!-- image -->  \n2  \nGM  \n1  \n1  \n/lessmuch  \nl  \nl  \n2  \ns  \ns  \nr  \nh  \ne  \ne  \ny  \nm  \nr  \nl  \ne  \n-  \n2  \nGM ρ  \n=  \nconst. r  \n= 0  \nr  \n=  \nr  \nd  \ns  \nc  \nf  \ng  \nh  \na  \n-  \nGM  \n2  \nh  \nz  \nw  \nS  \n2  \n)  \n-  \n1  \nR  \nm  \nρ = r s e . T w o c a l e t ∼ H -1 , t h e d o m a i n w a l l R /lessmuch t σ , t h e d o m a i n w a l l i s w a l l f a l l s w i t h i n t h e h o r i z o n , i t = 2 GM ∼ t 2 /t σ /lessmuch t . r e R /greatermuch t σ . I n t h i s c a s e , t h e w a l l r e p e l s e n t e x p a n s i o n r a t e . A s w e s h a l l s e e , t h i s o r i g i n a l l y i n t h e i n t e r i o r o f t h e w a l l g o e s i n t o a r e l e f t w i t h a b l a c k h o l e r e m n a n t c y s t e d i n a a l l , l e t u s fi r s t d i s c u s s t h e m a t c h i n g o f S c h w a r z s c h i l d a dus t co s m o l o g y ( -2 GM R ) dT 2 + ( 1 -2 GM R ) -1 dR 2 + R 2 d Ω 2 . ( 1 ) t h e f o r m -d τ 2 +2 GM R d ρ 2 + R 2 d Ω 2 . GM ( 1 -2 GM R ) -1 dR , 2 GM R ) -1 R ρ 0 t C o n s d i ff e r e n t t i a n d t h e o t h e t σ ∼ ( G σ ) -1 . F o r R /lessmuch t σ , t c o n f o r m a l l y s t r e t c h e d q u i c k l y s h r i n k s u n d e r i t H e r e , w e w i l l b e p r t h e m a t t e r a r o u n d i t w l e a d s t o t h e f o r m a t i o n o f a a b a b y u n i v e r s e , a n d i n t h e s p h e r i c a l r e g i o n o f o f v a c u u m . B e f o r e w e c o n s i d e r t h e e ff e c t o a n d a d u s t c o s m o l o g y . 1 . 1 M a t c hi ng Sc h w a r zs c hi l d t o a C o n s i d e r t h e S c h w a r z s c h i l d m e t r i c ds 2 = -( 1 -2 GM R ) dT 2 + ( I n L e m a i t r e c o o r d i n a t e s , t h i s t a k e s t h e f o r m ds 2 = -d τ 2 + 2 GM R d ρ 2 + R h e r e τ a n d ρ a r e d e fi n e d b y t h e r e l a t i o n s d τ = ± dT -√ 2 GM R ( 1 -2 GM R ) -1 dR d ρ = ∓ dT + √ R 2 GM ( 1 -2 GM R ) -1 ( 4 ) , w e h a v e R = [ 3 2 ( τ + ρ ) ] 2 / 3 ( 2 h a s b e e n a b s o r b e d a n d ρ c a n b e R different time-scales are relevant to the dynamics of such a wall. One is the cosmological scale and the other is the acceleration time-scale due to the repulsive gravitational field of the domain wall t σ ∼ ( Gσ ) -1 . For R /lessmuch t σ , the repulsive field can be ignored. In this case, for t /lessmuch R /lessmuch t σ , the domain wall is conformally stretched by cosmological expansion. Eventually, when the wall falls within the horizon, it quickly shrinks under its tension and forms a black hole of radius R S = 2 GM ∼ t 2 /t σ /lessmuch t . Here, we will be primarily interested in the opposite limit, where R /greatermuch t σ . In this case, the wall repels the matter around it while its size grows faster than the ambient expansion rate. As we shall see, this leads to the formation of a wormhole. The dust which was originally in the interior of the wall goes into a baby universe, and in the ambient FRW universe we are left with a black hole remnant cysted in a spherical region of of vacuum. Before we consider the effect of the domain wall, let us first discuss the matching of Schwarzschild and a dust cosmology. 1.1 Matching Schwarzschild to a dust cosmology Consider the Schwarzschild metric ds 2 = -( 1 -2 GM R ) dT 2 + ( 1 -2 GM R ) -1 dR 2 + R 2 d Ω 2 . (1) FIG. 3: Formation of a black hole by a large vacuum bubble, with M bh glyph[greatermuch] M cr , in a dust dominated spatially flat FRW. In this case, the bubble does not collapse into the Schwarzschild singularity. Instead, at the time t cr ∼ Min { t b , t σ } the size of the bubble starts growing exponentially in a baby universe, which is connected by a wormhole to the parent dust dominated FRW universe. Initially, a geodesic observer at the edge of the dust region can send signals through the wormhole into the baby universe. This is represented by the blue arrow in the Figure. However, after a proper time t ∼ 2 GM bh , the wormhole 'closes' and any signals which are sent radially inwards end up at the Schwarzschild singularity.  \nIn Lemaitre coordinates, this takes the form  \nds  \n2  \n=  \n-  \ndτ  \nf  \no  \n)  \n(  \nn  \ng  \nt  \n+  \nR  \nr  \nu  \nco  \ndus  \nt  \ndρ  \n2  \nS  \ni  \n2  \na  \nb  \n-  \n-  \n(2  \nGM  \n)  \n2  \nGM  \nR  \n2  \nGM  \nR  \n1  \n/  \n3  \nd  \nΩ  \nr  \n3  \n2  \n)  \na  \ne  \ni  \na  \nr  \n.  \n(2)  \n-  \nn  \ne  \n1  \ne  \ns  \ndR,  \n(3)  \np  \nh  \nw  \n)  \n,  \n(5)  \nwhere an integration constant has been absorbed by a shift in the origin of the  \nρ  \ncoordinate.  \nThe  \ncan be found from (3) as:  \nR  \n-  \ntanh  \n-  \n√  \n2  \nGM  \n-  \nρ,  \n(6)  \ngiven in (5). A second integration constant has been absorbed by a shift in the  \nT  \ncoordinate.  \nSince the metric (2) is synchronous, the lines of constant spatial coordinate are geodesics, and  \nτ  \nis the  \n-  \n1  \ni  \n.  \nt  \nt  \no  \nn  \nt  \ng  \nn  \ni  \nr  \nt  \no  \nf  \no  \ne  \nx  \nn  \ns  \nr  \ne  \ndR.  \n(4)  \nΩ  \nn  \ns  \nc  \no  \ns  \nT  \na  \nh  \nn  \ni  \nR  \nt  \nn  \ni  \nτ  \nH  \n,  \n∼  \n1  \n-  \n1  \n,  \n.  \ndR  \nGM  \ny  \nb  \nu  \no  \n2  \nf  \n/  \n3  \nGM  \ng  \nn  \n(  \nGM  \n4  \n.  \n(  \ni  \ne  \ng  \nv  \ne  \n=  \n±  \no  \no  \nT  \n5  \nn  \n2  \n3  \nd  \nc  \no  \ns  \ne  \ni  \nn  \ne  \nA  \n)  \nt  \no  \nr  \nn  \nn  \nh  \ny  \ns  \no  \nc  \ns  \ni  \n2  \nc  \n(  \ni  \ne  \nr  \nm  \nh  \ne  \ni  \n)  \nt  \nt  \nc  \nS  \nr  \nw  \na  \nn  \ni  \nn  \nf  \nn  \nt  \nu  \nf  \nc  \na  \nc  \nt  \n1  \n2  \nd  \nw  \n(  \n5  \n)  \nm  \nm  \nS  \ne  \nn  \nc  \nf  \ni  \nh  \no  \nm  \nt  \nh  \nR  \n/lessmuch  \nl  \nρ  \n=  \nd  \nh  \nf  \nd  \nfi  \no  \nl  \nR  \no  \ni  \ns  \nr  \nr  \nd  \ni  \nn  \na  \n(  \nt  \n6  \ne  \n)  \n.  \nT  \nh  \ne  \nwhere  \nτ  \nand  \nρ  \nare defined by the relations  \nst  \n2  \n+  \n√  \n2  \nGM  \n2  \n/  \n3  \n2  \nGM  \nR  \n(  \n(  \n1  \nR  \ndR  \nh  \ni  \ns  \na  \nl  \nd  \nc  \nr  \n2  \nd  \nΩ  \nR  \n2  \n+  \n2  \n.  \n0  \n)  \n1  \n(  \nρ  \n=  \nqualitative analysis of this case, since a detailed description will require further numerical studies which are beyond the scope of the present work.  \nInitially, as mentioned in Subsection II A, the Lorentz factor of the bubble wall relative to the Hubble flow at the end of inflation is of order γ i ∼ t σ /t i glyph[greatermuch] 1. As the wall advances into ambient radiation with energy density ρ m , a shock wave forms. The shocked fluid is initially at rest in the frame of the wall, and its proper density is of order ρ s ∼ γ 2 ρ m , where γ is the Lorentz factor of the shocked fluid relative to the ambient [11]. In the frame of shocked radiation, the shock front moves outward at a subsonic speed β ≈ 1 / 3, and therefore the pressure of the shocked gas has time to equilibrate while the wall is pumping energy into it. Consider a small region of size  \nL glyph[lessmuch] t i (39)  \nnear the wall. At such scales, the wall looks approximately flat, and gravitational effects (including the expansion of the universe) can be neglected. In the frame of ambient radiation, energy conservation per unit surface of the bubble wall can be expressed as  \nγ i σ ∼ γ 2 ρ m ∆ t + γσ. (40)  \nThe left hand side of this relation represents the initial energy of the wall. The first term in the right hand side is the energy of the layer of shocked fluid. This is of order γρ s ∆ τ ∼ γ 2 ρ m ∆ t , where ∆ τ is the proper time ellapsed since the shock starts forming, and ∆ t = ( t -t i ) = γ ∆ τ is the corresponding interval in the ambient frame. The second term in the right hand side is just the energy of the wall at time t . Here, we have neglected the rest energy of the ambient fluid, which is of order γ -2 relative to the shocked fluid.  \nEq. (40) shows that the bubble wall looses most of its energy on a time scale ∆ t ∼ σ/ ( γ i ρ m ) ∼ ( t i /t σ ) 2 t i glyph[lessmuch] t i . This is the characteristic time which is required for the pressure of the shocked fluid to accelerate the wall away from it, and is actually in agreement with the rough estimate (18), which is based on a simplified fluid model. The recoil velocity of the wall cannot grow too large, though, because the backward pressure wave can only travel at the speed of sound. Hence, the relative Lorentz factor γ w between the wall and the shocked fluid will be at most of order one. The initial energy of the wall per unit surface is is γ i σ i ∼ G/t i . After time of order t i , this energy will be distributed within a layer of width t i near the position of the wall, and its density will be comparable to the ambient energy  \ndensity ρ m ∼ G/t 2 i . By that time, the shock will have dissipated, and the wall will have, at most, a mildly relativistic speed with respect to the ambient Hubble flow.  \nFor t > ∼ t i , the expansion of the universe, and of the bubble, start playing a role. The wall's motion becomes then dominated by the effect of its tension, and of the vacuum energy density in its interior. Again, the fate of gravitational collapse depends on whether the bubble is larger or smaller than a critical size.", '1. Small bubbles surrounded by radiation': "Shortly after t ∼ t i , the bubble wall will be recoiling from ambient radiation at mildly relativistic speed, due to the pressure wave which proceeds inward at the speed of sound, in the aftermath of the shock. For t > ∼ t i , we expect that the wall will only be in contact with rarified radiation dripping ahead of the pressure wave. Neglecting the possible effect of such radiation, the dynamics proceeds as if the bubble were in an empty cavity, along the lines of our discussion in Subsection II C. In that case, the bubble has a conserved energy M bh given approximately by Eq. (20). For M bh glyph[lessmuch] M cr , the bubble initially grows until it reaches a turning point, at some maximum radius R max , after which it collapses to form a black hole of mass M bh . As mentioned after Eq. (37), for R i < ∼ R glyph[lessmuch] R max , the radius grows as  \nR ∝ τ 1 / 3 , (41)  \nwhere τ is proper time on the wall. Note that the ambient radiation of the FRW universe expands faster than that, as a function of its own proper time, since a ( t ) ∝ t 1 / 2 . This means that even in the absence of external pressure, the bubble motion decouples from the Hubble flow on the timescale t i , with a tendency for the gap between the bubble and the ambient Hubble flow to grow in time. In what follows, and awaiting confirmation from numerical studies, we assume that the influence of radiation on the wall's motion is negligible thereafter.  \nEventually the wall will collapse, driven by the negative pressure in its interior and by the wall tension. We can estimate the mass of the resulting black hole as the sum of bulk and bubble wall contributions, as in Eq. (20):  \nM bh ∼ ( 4 3 πρ b +4 πσH i ) R 3 i . (42)  \nIn a radiation dominated universe, the co-moving region of initial size R i crosses the horizon at the time  \nt H ∼ R 2 i t i , (43)  \nand we can write (42) as  \nwhere in the last relation we have used that for M bh glyph[lessmuch] M cr , we have t i < ∼ R i glyph[lessmuch] t b , t σ . Hence, the black hole radius is much smaller than the co-moving region affected by the bubble at the time of horizon crossing.  \nGM bh ∼ ( t i R i t 2 b + R i t σ ) t H glyph[lessmuch] t H , (44)", '2. Large bubbles surrounded by radiation': 'For M bh glyph[greatermuch] M cr , the situation is more complicated. A wormhole starts developing at t ∼ GM cr , when the growth of the bubble becomes exponential. On the other hand, the wormhole stays open during a time-scale of order ∆ t ∼ GM bh glyph[greatermuch] GM cr . During that time, some radiation can flow together with the bubble into the baby universe, thus changing the estimate (42) for M bh . On general grounds, we expect that the radius of the black hole is at most of the size of the cosmological horizon at the time t H when the co-moving scale corresponding to R i crosses the horizon  \nGM bh < ∼ t H . (45)  \nThe study of modifications to Eq. (42) for M bh glyph[greatermuch] M cr requires numerical analysis and is left for further research.', 'III. GRAVITATIONAL COLLAPSE OF DOMAIN WALLS': 'In the previous sections we considered bubbles with ρ b glyph[lessmuch] ρ i . This condition led to a highly boosted bubble wall at the time when inflation ends. Here we shall concentrate instead on domain walls, which correspond to the limit ρ b = ρ i . Domain walls will be essentially at rest with respect to the Hubble flow at the end of inflation. Another difference is that for the case of bubbles we assumed that the inflaton field lives only outside the bubble, so its energy is transfered to matter that lives also outside. Here, we will explore the alternative possibility where matter is created both inside and outside the domain wall. This seems  \nto be a natural choice, since domain walls are supposed to be the result of spontaneous symmetry breaking, and we may expect the physics of both domains to be essentially the same. As we shall see, the presence of matter inside the domain wall has a dramatic effect on its dynamics.  \nFor t glyph[lessmuch] t σ , the graviational field of the walls is completely negligible. Walls of superhorizon size R glyph[greatermuch] t i are initially at rest with respect to the Hubble flow, and are conformally stretched by the expansion of the universe, R ∝ a . Depending on their initial size, their fate will be different. As we shall see, small walls enter the horizon at a time t H glyph[lessmuch] t σ and then they begin to decouple from the Hubble flow. Depending on the interaction of matter with the domain wall, they may either collapse to a black hole singularity, or they may develop long lived remnants in the form of pressure supported bags of matter, bounded by the wall. Larger walls with t H glyph[greatermuch] t σ start creating wormholes at the time t ∼ t σ , after which the walls undergo exponential expansion in baby universes.', 'A. Domain walls surrounded by dust': 'The case of walls surrounded by dust is particularly simple, since the effect of matter outside the wall can be completely ignored. As we did for bubbles, here we may distinguish between small walls, for which self-gravity is negligible, and large walls, for which it is important.', '1. Small walls': "Let us start by discussing the situation where the domain wall is 'impermeable'. Small walls with t H glyph[lessmuch] t σ continue to be dragged by the Hubble flow even after they have crossed the horizon. Unlike dust, radiation outside the wall may continue exerting pressure on it, and initially its effect cannot be neglected. The radius of the wall will stop growing once the wall tension balances the difference in pressure ∆ p between the inside and the outside:  \n∆ p ∼ σ R . (68)  \nAfter that, the pressure outside decays as ( t H /t ) 2 and becomes negligible compared with the pressure inside, which is given by p ∼ (1 /Gt 2 H )( t H /R ) 4 . From (68) we then find that the equilibrium radius is of order  \nR bag ∼ ( t σ /t H ) 1 / 3 t H glyph[greatermuch] t H . (69)  \nThe mass of the bag of radiation is of order  \nM bag ∼ t H G ( t H t σ ) 1 / 3 . (70)  \nTherefore the gravitational potential at the surface of the bag is very small Φ ∼ ( t H /t σ ) 2 / 3 glyph[lessmuch] 1.  \nhave γ w = 1 + O ( t 2 ∗ /t 2 σ ), or v ∼ ( t ∗ /t σ ) glyph[lessmuch] 1.  \nNext, we may consider the case of 'permeable' walls. Like in the case of dust, the dynamics of small walls is well described by the Nambu action in a FRW universe. The mass of the resulting black holes can be estimated as  \nM bh ≈ 16 πC r σt 2 H , (71)  \nwhere in a radiation dominated universe t H is determined from the condition  \nR ( t H ) = 2 t H . (72)  \nThe numerical coefficient is found by using the numerical approach described in Appendix A, and we find  \nC r ≈ 0 . 62 . (73)  \nNote that, although permeable walls do form black holes, both in the dust dominated and radiation dominated eras, their masses are much smaller than the masses of the corresponding bags of matter trapped by impermeable walls.", '2. Large domain walls surrounded by dust': "If a domain wall has not crossed the horizon by the time t σ , its radius starts growing exponentially and decouples from the Hubble flow, forming a wormhole.  \nTo illustrate this process, we may consider an exact solution where a spherical domain wall of initial radius  \nR ∗ > t σ , (58)  \nis surrounded by pressuereless dust. The solution is constructed as follows. At some chosen time t ∗ , we assume that the metric inside and outside the wall is initially a flat FRW metric with scale factor a ∝ t 2 / 3 . We shall also assume that the wall is separated from matter, on both sides, by empty layers of infinitesimal width. (As we shall see, these layers will grow with time.) By symmetry, the metric within these two empty layers takes the Schwarzschild form  \nThe mass parameter M -for the interior layer can in principle be different from the mass parameter M + for the exterior. Let us now find such mass parameters by matching the metric in the empty layers to the adjacent dust dominated FRW.  \nds 2 = -( 1 -2 GM ± R ) dT 2 + ( 1 -2 GM ± R ) -1 dR 2 + R 2 d Ω 2 . (59)  \nBy continuity, the particles of matter at the boundary of a layer must follow a geodesic of both Schwarzschild and FRW. Such geodesics originate at the white-hole/cosmological singularity, and they end at time-like infinity (see Fig. 6). The Schwarzschild metric is independent of T , and geodesics satisfiy the conservation of the corresponding canonical momentum:  \n˙ T = ( 1 -2 GM ± R ) -1 C ± . (60)  \nHere C ± is a constant, and a dot indicates derivative with respect to proper time. Initially, the geodesic is in the white hole part of the Schwarzschild solution (Region I in Fig. 6), and C ± is positive or negative depending on whether the particle moves to the right or to the left. Using (1 -2 GM ± /R ) ˙ T 2 -(1 -2 GM ± /R ) -1 ˙ R 2 = 1, we have  \n˙ R 2 = 2 GM ± R + C 2 ± -1 . (61)  \nFor geodesics with the escape velocity, we have 8  \nC ± = 1 . (62)  \nThe equation of motion (61) then takes the form  \n( ˙ R R ) 2 = 2 GM ± R 3 . (63)  \nThis has the solution R ∝ t 2 / 3 , where t is proper time along the geodesic of the dust particle, matching the behaviour of the flat FRW expansion factor. The geodesics at the boundaries of the two infinitesimal empty layers initially coincide at R = R ∗ . Assuming that the initial matter density is the same inside and outside the wall, comparison of (63) with the Friedmann equation leads to the conclusion that M + and M -have the same value:  \nM bh ≡ M ± = 4 π 3 ρ m ( t ∗ ) R 3 ∗ = 3 t H 4 G . (64)  \nHere, M bh is simply the total mass of matter inside the domain wall, which is also the mass which is excised from the exterior FRW, and t H is the time at which the boundary of such excised region crosses the horizon.  \nThe dynamics of the bubble wall in between the two empty layers is determined from Israel's matching conditions, which lead to the equations [9]  \n˙ T = ± ( 1 -2 GM bh R ) -1 R t σ , (65)  \nC  \n1  \no  \ni  \n.  \nn  \ne  \na  \n1  \ns  \nc  \nfi  \nn  \nd  \no  \nd  \ne  \no  \ni  \ns  \nd  \nM  \nr  \nr  \ne  \nb  \nd  \nt  \na  \nh  \nn  \ne  \na  \nt  \nh  \nt  \nt  \n,  \nh  \nco  \ni  \ns  \nt  \ne  \nn  \n=  \ndus  \ni  \nτ  \nl  \nl  \nw  \na  \n∓  \n+  \nh  \nt  \nl  \na  \np  \nh  \ne  \na  \nb  \ne  \nB  \nn  \na  \ne  \nd  \na  \nr  \ns  \nb  \ni  \ne  \nq  \nc  \nu  \no  \ni  \nt  \nn  \nc  \nH  \na  \nσ  \nf  \nk  \no  \nl  \ne  \na  \nm  \nt  \ny  \nc  \nf  \nd  \no  \nu  \nt  \ne  \nhi  \nc  \nS  \nng  \nw  \nl  \na  \nr  \ns  \ne  \nt  \nc  \no  \nu  \ny  \nr  \nt  \ne  \n,  \nt  \nt  \nn  \nr  \ne  \nh  \ni  \ne  \nw  \no  \nr  \ne  \nv  \ng  \ne  \ni  \ne  \ns  \nw  \na  \nf  \nr  \no  \nc  \nm  \nr  \no  \ns  \nn  \no  \no  \nn  \nl  \nSc  \nz  \nr  \nc  \n2  \nr  \ne  \n,  \no  \nf  \ns  \no  \ni  \ng  \nh  \n1  \nm  \na  \no  \nd  \ny  \nn  \nf  \ne  \n.  \nw  \na  \ns  \ni  \ny  \nh  \nl  \nr  \nzs  \nm  \no  \na  \nd  \ng  \nl  \n-  \nds  \no  \n=  \ns  \nm  \na  \ne  \nt  \n(  \ne  \nk  \ns  \n2  \nds  \nt  \ne  \na  \ndT  \n±  \ndT  \nρ  \nR  \n= 2  \nGM  \nρ  \n1  \nc  \nr  \nr  \nv  \nt  \nh  \nhi  \ni  \n-  \nt  \nt  \n= 0  \n2  \n2  \nf  \nc  \nGM  \nR  \no  \n=  \n0  \nh  \ne  \nt  \n-  \nρ  \ns  \np  \nGM  \na  \nn  \ni  \nr  \nd  \nσ  \nd  \n∼  \nm  \n(  \nGσ  \nFor  \nR  \ne  \n2  \n∞  \nd  \nd  \n-  \ne  \nSphe  \nR  \n1  \n2  \n=  \n=  \nR  \n0  \nρ  \nt  \nC  \nff  \nd  \ni  \nn  \nt  \nGM  \n0  \ns  \no  \ne  \nn  \nn  \nr  \ng  \no  \ns  \nρ  \n=  \n-  \n2  \nt  \nρ  \n=  \nr  \nt  \nt  \nm  \nGM  \na  \n.  \ncon  \n=  \n0  \nd  \na  \nm  \nst  \ni  \no  \no  \ne  \nn  \nr  \nt  \ng  \no  \nv  \nu  \nr  \nn  \n.  \nr  \ni  \ne  \ns  \ne  \nt  \ne  \n-  \nw  \nH  \nn  \nn  \na  \nl  \nl  \nw  \ni  \ni  \na  \ns  \nl  \nl  \nn  \ni  \nh  \nm  \ni  \na  \no  \ne  \nd  \nh  \nh  \nσ  \nm  \nh  \no  \ne  \n.  \nt  \nr  \nh  \nt  \ni  \nl  \ns  \nc  \na  \nc  \nl  \ni  \nl  \no  \nl  \no  \ns  \nm  \nd  \ne  \nh  \ne  \nc  \nn  \no  \nt  \na  \n∼  \na  \no  \ne  \nd  \nh  \nt  \nσ  \nt  \n,  \nw  \ni  \nn  \nt  \n/t  \n/lessmuch  \nt  \n,  \na  \ns  \ne  \nc  \n2  \nh  \ni  \ns  \ns  \nt  \nt  \nI  \n/lessmuch  \nl  \n/lessmuch  \nh  \ne  \nfi  \nf  \nd  \nR  \na  \nl  \nw  \nf  \na  \n2  \nl  \nl  \n∼  \n.  \nn  \nA  \ne  \nr  \no  \nn  \ne  \n.  \na  \nt  \nr  \ni  \nt  \ne  \ni  \nn  \no  \ns  \nw  \nz  \no  \na  \ne  \nw  \na  \ns  \nh  \ne  \ne  \nh  \no  \nf  \nr  \ne  \nt  \ne  \nm  \nt  \na  \nr  \nl  \nm  \nh  \ne  \no  \nk  \nl  \nGM  \nt  \nσ  \ni  \na  \n/greatermuch  \n=  \nR  \nS  \ne  \nr  \nR  \nx  \np  \nn  \ns  \nt  \ni  \ne  \nl  \ny  \nn  \nh  \nd  \na  \nh  \nc  \ns  \nc  \ns  \ns  \ny  \na  \ne  \nn  \ne  \nt  \ni  \na  \ng  \ni  \nn  \nl  \nw  \ni  \nh  \nb  \nu  \ni  \nt  \na  \nt  \ns  \nfi  \nρ  \n=  \nb  \ni  \nl  \no  \nr  \ne  \nr  \ne  \nl  \na  \nl  \nl  \no  \ng  \no  \n,  \n→-∞  \ns  \nco  \n2  \nGM  \nGM  \n2  \nSpherical domain wall in dust cosmology  \nm  \nl  \n-  \nConsider a spherical domain wall embedded in a spatially flat matter  \ndifferent time-scales are relevant to the dynamics of such a wall. One is the c  \nt  \nR  \nand the other is the acceleration time-scale due to the repulsive gravitation  \n→-∞  \nρ  \n=  \n2  \n-  \na  \n1  \n+  \n(  \nConsider a spherical domain wall embedded in a spatially flat matte  \n)  \nIn this case, for  \nt  \n(  \n4  \ndT  \ndifferent time-scales are relevant to the dynamics of such a wall. One is the  \ndus  \nGM  \n)  \n2  \n2  \nGM  \n-  \n/lessmuch  \nR  \nconformally stretched by cosmological expansion. Eventually, when the wall  \n/lessmuch  \nt  \na  \n<!-- image -->  \n)  \nf  \ns  \no  \nc  \ni  \n+2  \nm  \nt  \no  \ne  \ne  \n,  \nthe repulsive field can be ignored.  \nIn this case, for  \nt  \nI  \ni  \nn  \n1  \n.  \nd  \n(  \ne  \nσ  \n=  \nd  \nu  \nt  \ne  \nd  \nn  \nS  \nt  \n=  \nt  \no  \nf  \ny  \ne  \nv  \no  \nl  \ni  \ni  \nl  \nm  \nr  \ni  \nca  \nb  \na  \nm  \nl  \ndo  \ni  \ni  \ng  \no  \ne  \nf  \nm  \nu  \ns  \ne  \nr  \nt  \ns  \nl  \nd  \ni  \ns  \nc  \nk  \nt  \ns  \n/lessmuch  \nand the other is the acceleration time-scale due to the repulsive gravitati  \nt  \n.  \nΩ  \nm  \n2  \nd  \n∼  \nh  \n/lessmuch  \nFor  \nR  \nR  \nt  \nh  \nr  \nm  \nf  \n)  \nR  \no  \n5  \nS  \n(  \nt  \nt  \nn  \ne  \nn  \nt  \nu  \nσ  \n,  \nthe domain wall is  \n-  \nσ  \n,  \nthe repulsive field can be ignored.  \nIn this case, for  \nt  \n1  \ne  \na  \n,  \ny  \nl  \nw  \nd  \n2  \n-  \n1  \nHere, we will be primarily interested in the opposite limit, where  \nR  \nt  \nt  \na  \ni  \ne  \no  \nR  \nh  \na  \nc  \ns  \n2  \nd  \nτ  \n= 2  \nGM  \n=  \n+  \n-  \nt  \nσ  \nh  \ne  \nn  \n/lessmuch  \na  \na  \ns  \na  \n/lessmuch  \n/lessmuch  \ne  \nl  \ns  \ne  \no  \nA  \nr  \nt  \nh  \nl  \nb  \nr  \ni  \ne  \no  \nl  \ni  \nh  \ne  \nr  \ne  \nc  \nw  \nr  \nu  \ne  \nm  \nc  \no  \ns  \nconformally stretched by cosmological expansion. Eventually, when the wall falls within the horizon, it  \n,  \nh  \nw  \n.  \nr  \nt  \nh  \no  \nf  \nl  \ne  \ne  \n-  \nr  \no  \n+  \no  \nSc  \nh  \ni  \ns  \na  \nl  \n(  \ns  \nR  \ne  \nh  \nt  \n.  \nconformally stretched by cosmological expansion. Eventually, when the wa  \nw  \n2  \nτ  \nb  \nthe matter around it while its size grows faster than the ambient expansion  \nquickly shrinks under its tension and forms a black hole of radius  \nR  \ne  \n=  \ns  \no  \nn  \nT  \n→-  \n,  \na  \nm  \nl  \n2  \nl  \nw  \nl  \ni  \n-  \n2  \nGM  \n=  \n+  \n2  \n√  \nR  \ni  \nn  \nR  \nm  \na  \no  \nR  \n√  \nd  \nl  \nGM  \na  \ni  \n→∞  \n)  \n,  \nT  \nc  \ne  \nr  \n/  \nh  \nleads to the formation of a wormhole. The dust which was originally in the interior of the wall goes into  \nc  \nr  \nw  \nl  \no  \na  \ni  \nr  \n= 2  \nGM  \nt  \n2  \nS  \nb  \ni  \n-  \ns  \ni  \ne  \nt  \nh  \n2  \ne d y -s c a l e e i g n a n s i o n . E a b l a c k h o p p o s i t e r t h a n t h w h i c h w a s s e w e a r e l e t u s fi r g y -1 dR n a t e . ∼ /lessmuch Here, we will be primarily interested in the opposite limit, where R /greatermuch t σ . In this case, the wall repels the matter around it while its size grows faster than the ambient expansion rate. As we shall see, this Here, we will be primarily interested in the opposite limit, where R t σ the matter around it while its size grows faster than the ambient expansio i t l r e p e l s e s l l s e e , t t h e w a l l g o e r e m n a n t c y s t m a t c h i n g o f R ) -1 dR 2 GM R d ρ 2 a n d a M a t r t h e S c n a t e s , t t h e r e l a leads to the formation of a wormhole. The dust which was originally in the i a baby universe, and in the ambient FRW universe we are left with a black spherical region of of vacuum. FIG. 6: A large domain wall causing the birth of a baby universe.  \nt  \no  \nGM  \nt  \na  \nn  \nv  \nl  \ne  \nr  \ne  \ne  \n(  \nr  \ns  \nh  \nρ  \n2  \np  \ns  \n]  \na  \nr  \ny  \ne  \na  \ne  \n3  \nl  \nc  \na  \n(  \ns  \ns  \n2  \n-  \na  \nr  \nd  \ne  \ni  \ni  \nm  \nh  \nf  \n(  \nf  \nm  \nt  \nn  \na  \ne  \nt t i o n t i ( 1 -and  \nr  \nl  \ne  \n1  \ne  \nc  \nc  \n-  \na  \nh  \ne  \nt  \ni  \ns  \nr  \nGM  \ne  \nt  \nh  \nt  \n3  \no  \nt  \ne  \nn  \nr  \na  \no  \ni  \nn  \na  \n)  \nh  \ne  \nt  \nn  \nd  \na  \ne  \nb  \nd  \nn  \ny  \nt  \ne  \nn  \nh  \nt  \n)  \n1  \n/  \n.  \n1  \n)  \n-  \ns  \nσ  \nt  \n:  \n(  \nσ  \nG  \nR  \ne  \n,  \n/lessmuch  \nh  \no  \nr  \n3  \nt  \nl  \ni  \na  \nl  \ng  \no  \nr  \nr  \nm  \nf  \n∼  \n-  \n1  \n2  \nF  \nn  \no  \no  \ni  \nc  \nk  \nu  \nc  \nq  \n-  \ni  \n√  \nR  \nGM  \na  \na  \ns  \nr  \ne  \nh  \ni  \nf  \ng  \nt  \ne  \no  \nn  \nd  \nt  \ns  \nρ  \n,  \nh  \ni  \ne  \ns  \n,  \nT  \na  \nn  \nd  \nc  \ne  \na  \nm  \nd  \nu  \nb  \na  \nx  \nn  \nc  \np  \ne  \nR  \nGM  \ns  \n2  \nr  \nh  \ne  \nt  \nm  \ni  \ns  \nt  \ne  \na  \nf  \n+  \nu  \nt  \ne  \nr  \n2  \nn  \nR  \nΩ  \ns  \nd  \nv  \n2  \n,  \na  \nl  \nl  \nl  \no  \ne  \nf  \no  \nt  \nd  \nt  \n.  \ns  \nng  \nh  \nz  \nw  \n2  \n=  \na  \nk  \nquickly shrinks under its tension and forms a black hole of radius  \nR  \ne  \na  \nh  \ne  \na  \ns  \ni  \n+  \nl  \ni  \nds  \ny  \nd  \n)  \ni  \nr  \nd  \nh  \nt  \n2  \ns  \na  \ns  \np  \n.  \nleads to the formation of a wormhole. The dust which was originally in the  \nc  \nhi  \ne  \ns  \ni  \no  \ni  \no  \nz  \ns  \nn  \nd  \nS  \nw  \nz  \n2  \no  \nb  \ne  \no  \n1  \n.  \nc  \nh  \n1  \nd  \nds  \ni  \no  \nd  \na baby universe, and in the ambient FRW universe we are left with a black hole remnant cysted in a  \nR  \ngiven in (5). A second integration constant has been absorbed by a shift in the  \nT  \ncoordinate.  \nwith  \nR  \ngiven in (5). A second integration constant has been absorbed by a  \nr  \nd  \no  \na baby universe, and in the ambient FRW universe we are left with a bla  \na  \ns  \n/greatermuch  \no  \n= 2  \n/greatermuch  \n=  \nBefore we consider the effect of the domain wall, let us first discuss the  \nwith  \nR  \ngiven in (5). A second integration constant has been absorbed by  \nT  \nS  \nt  \ni  \nτ  \na σ , t h e r e p u l s i v e fi e l d y s t r e t c h e d b y c o s m o l o g i c a l y s h r i n k s u n d e r i t s t e n s i o n a n d f H e r e , w e w i l l b e p r i m a r i l y i n t e r e s t e d t h e m a t t e r a r o u n d i t w h i l e i t s s i z e g r o w s l e a d s t o t h e f o r m a t i o n o f a w o r m h o l e . T h e a b a b y u n i v e r s e , a n d i n t h e a m b i e n t F R W u s p h e r i c a l r e g i o n o f o f v a c u u m . B e f o r e w e c o n s i d e r t h e e ff e c t o f t h e d o m a i n a n d a d u s t c o s m o l o g y . 1 . 1 M a t c hi ng Sc h w a r zs c hi l d t o a dus t co s C o n s i d e r t h e S c h w a r z s c h i l d m e t r i c ds 2 = -( 1 -2 GM R ) dT 2 + ( 1 -2 R I n L e m a i t r e c o o r d i n a t e s , t h i s t a k e s t h e f o r m ds 2 = -d τ 2 + 2 GM R d ρ 2 + R 2 d w h e r e τ a n d ρ a r e d e fi n e d b y t h e r e l a t i o n s d τ = ± dT -√ 2 GM R ( 1 -2 GM R ) -1 dR d ρ = ∓ dT + √ R 2 GM ( 1 -2 GM R ) -1 dR . S u b t r a c t i n g ( 3 ) f r o m ( 4 ) , w e h a v e R = [ 3 2 ( τ + ρ ) ] 2 / 3 ( 2 GM ) 1 / 3 , w h e r e a n i n t e g r a t i o n c o n s t a n t h a s b e e n a b s o r b e d b y a s h i f t i n t h e o r i g i n o f t h e ρ x p r e s s i o n f o r T a s a f u n c t i o n o f τ a n d ρ c a n b e f o u n d f r o m ( 3 ) a s : ± T = 4 GM   1 3 ( R 2 GM ) 3 / 2 + √ R 2 GM -t a n h -1 √ R 2 GM   -ρ , i v e n i n ( 5 ) . A s e c o n d i n t e g r a t i o n c o n s t a n t h a s b e e n a b s o r b e d b y m e t r i c ( 2 ) i s s y n c h r o n o u s , t h e l i n e s o f c o n s t a n t s p a t i a l c o o r d spherical region of of vacuum. Before we consider the effect of the domain wall, let us first discuss the matching of Schwarzschild and a dust cosmology. 1.1 Matching Schwarzschild to a dust cosmology Consider the Schwarzschild metric ds 2 = -( 1 -2 GM R ) dT 2 + ( 1 -2 GM R ) -1 dR 2 + R 2 d Ω 2 . (1) In Lemaitre coordinates, this takes the form ds 2 = -dτ 2 + 2 GM R dρ 2 + R 2 d Ω 2 . (2) where τ and ρ are defined by the relations dτ = ± dT -√ 2 GM R ( 1 -2 GM R ) -1 dR, (3) dρ = ∓ dT + √ R 2 GM ( 1 -2 GM R ) -1 dR. (4) Subtracting (3) from (4), we have R = [ 3 2 ( τ + ρ ) ] 2 / 3 (2 GM ) 1 / 3 , (5) where an integration constant has been absorbed by a shift in the origin of the ρ coordinate. The expression for T as a function of τ and ρ can be found from (3) as: ± T = 4 GM   1 3 ( R 2 GM ) 3 / 2 + √ R 2 GM -tanh -1 √ R 2 GM   -ρ, (6) with spherical region of of vacuum. Before we consider the effect of the domain wall, let us first discuss th and a dust cosmology. 1.1 Matching Schwarzschild to a dust cosmology Consider the Schwarzschild metric ds 2 = -( 1 -2 GM R ) dT 2 + ( 1 -2 GM R ) -1 dR 2 + R In Lemaitre coordinates, this takes the form ds 2 = -dτ 2 + 2 GM R dρ 2 + R 2 d Ω 2 . where τ and ρ are defined by the relations dτ = ± dT -√ 2 GM R ( 1 -2 GM R ) -1 dR, dρ = ∓ dT + √ R 2 GM ( 1 -2 GM R ) -1 dR. Subtracting (3) from (4), we have R = [ 3 2 ( τ + ρ ) ] 2 / 3 (2 GM ) 1 / 3 , where an integration constant has been absorbed by a shift in the origi expression for T as a function of τ and ρ can be found from (3) as: ± T = 4 GM   1 3 ( R 2 GM ) 3 / 2 + √ R 2 GM -tanh -1 √ R 2 s c h i l d . ( 1 ) . ( 2 ) 2 GM R ) -1 dR , ( 3 ) 2 GM R ) -1 dR . ( 4 ) ( 5 ) i n o f t h e ρ c o o r d i n a t e . T h e ( 6 ) c o o r d i n a t e . τ i s t h e I n L e m a i t r e w h e r e τ a n d ρ a r e d e fi d S u b t r a c t i n g ( 3 ) f r o m ( 4 ) , w h e r e a n i n t e x p r e s s i o n and a dust cosmology. 1.1 Matching Schwarzschild to a dust cosmology Consider the Schwarzschild metric ds 2 = -( 1 -2 GM R ) dT 2 + ( 1 -2 GM R ) -1 dR 2 + R 2 In Lemaitre coordinates, this takes the form ds 2 = -dτ 2 + 2 GM R dρ 2 + R 2 d Ω 2 . where τ and ρ are defined by the relations dτ = ± dT -√ 2 GM R ( 1 -2 GM R ) -1 dR, dρ = ∓ dT + √ R 2 GM ( 1 -2 GM R ) -1 dR. Subtracting (3) from (4), we have R = [ 3 2 ( τ + ρ ) ] 2 / 3 (2 GM ) 1 / 3 , where an integration constant has been absorbed by a shift in the origin expression for T as a function of τ and ρ can be found from (3) as: ± T = 4 GM   1 3 ( R 2 GM ) 3 / 2 + √ R 2 GM -tanh -1 √ R 2 GM ˙ R 2 = 2 GM bh R + ( R t σ ) 2 -1 . (66) The double sign ± in (65) refers to the embedding of the domain wall as seen from the interior and exterior Schwarzshild solutions respectively. The embedding is different, since there has to be a jump in the extrinsic curvature of the worldsheet as we go from one side of the wall to the other. Comparing (65) and (66) with (60) and (61) with C ± = 1 and M ± = M bh , we see that, for R > t σ , the wall moves faster than the dust particle at the edge of the expanding interior ball of matter. As seen from the outside, the wall moves in the direction opposite to the exterior FRW Hubble flow. Consequently, the wall never runs into matter and its motion is well described by (66) thereafter. For R > t σ Eq. (66) has no turning points and the radius continues to expand forever. The expansion reaches exponential behaviour R ∝ e τ/t σ after the time when the first term in the right hand side of (66) can be neglected and the second one becomes dominant. Here τ is the proper time on the worldsheet. This expansion is of course much faster than that of the Hubble flow, and takes place in the baby universe (see Figs. 6 and 7) . Loosely speaking, a wormhole in the extended Schwarzschild solution forms and closes on a time-  \ne  \nSince the metric (2) is synchronous, the lines of constant spatial coordinate are geodesics, and  \nτ  \nis the  \nh  \nR  \ng  \nh  \nc  \ne  \nt  \nSince the metric (2) is synchronous, the lines of constant spatial coordinat  \nSince the metric (2) is synchronous, the lines of constant spatial coordinate  \nw  \ne  \nn  \nt  \nc  \ni  \ni  \nS  \nd  \n,  \nC  \nc  \nΩ  \nn  \nn  \nGM  \nh  \nm  \nw  \ni  \nd  \nn  \nΩ  \no  \n2  \nl  \nd  \no  \nn  \nn  \nl  \n,  \na  \nr  \n/t  \nσ  \nB  \ne  \nr  \ns  \nt  \nt  \ne  \nf  \nu  \ni  \nt  \ns  \ne  \nl  \n,  \nt  \nc  \nh  \nr  \nd  \ny  \ns  \ny  \no  \nm  \nl  \nt  \nFIG. 7: A space-like slice of a baby universe with an expanding ball of matter in it, connected by a wormhole to an asymptotic dust dominated FRW universe. The wormhole is created by the repulsive gravitational field of an inflating domain wall, depicted as a red line.  \n<!-- image -->  \nscale ∆ t ∼ GM bh . This coincides with the time t H when the empty cavity containing the black hole crosses the horizon. An observer at the edge of the cavity can initially send signals into the baby universe, but not after the time when H -1 = 2 GM bh .  \nIn the above discussion, we did not make any assumptions about the time t ∗ . The exact solution illustrating the formation of a wormhole is obtained by a matching procedure, and can be constructed for any initial value of t ∗ . As we have seen, the condition R ∗ > t σ ensures that the wall will not run into matter. The wall evolves in vacuum, and the Hubble flow of matter remains undisturbed after t ∗ . Moreover, it is easy to show that for t ∗ glyph[lessmuch] t σ and R ∗ > t σ , the wall will be approximately co-moving with the Hubble flow at the time t ∗ , with a relative velocity of order v ∼ ( t ∗ /t σ ) glyph[lessmuch] 1. 9  \nγ w = -g µν U µ W ν = ( 1 -9 t 2 ∗ 4 R 2 ∗ )   √ 1 + 9 t 2 ∗ 4 t 2 σ -( 3 t ∗ 2 R ∗ ) 4 ∓ 9 t 2 ∗ 4 R ∗ t σ   . (67)  \nwhere we have used (60), (61), (65) and (66), with 2 GM bh = (4 / 9 t 2 ∗ ) R 3 ∗ . Using t ∗ glyph[lessmuch] t σ and R ∗ > t σ , we  \nNonetheless, in order to make contact with the setup of our interest, we must consider a broader set of initial conditions such that the wall can be smaller than t σ at the time t i when inflation ends, R i < t σ . In this case, the wall may initially expand a bit slower than the matter inside. Consequently some matter may come into contact with the wall. However, since superhorizon walls are approximately co-moving for t glyph[lessmuch] t σ , the effect of this interaction will be small, and we expect the qualitative behaviour to be very similar to that of the exact solution discussed above provided that the size of the wall becomes larger than t σ at some time t ∗ < ∼ t σ .", '2. Large domain walls surrounded by radiation': 'For t H ∼ t σ the estimates for M bh and M bag which we found in the previous subsection converge to the value  \nM bh ∼ t H G . (74)  \nFor t H glyph[greatermuch] t σ a wormhole will start forming near the time t σ , and some of the radiation may follow the domain wall into the baby universe. Nonetheless, the size of the resulting black hole embedded in the FRW universe cannot be larger than the cosmological horizon, and therefore M bh ≤ t H / 2 G . Continuity suggests that the estimate (74) may be valid also for large bubbles in a radiation dominated universe. A determination of M bh for this case requires a numerical simulation and is left for further research.', 'IV. MASS DISTRIBUTION OF BLACK HOLES': "In the earlier sections we have outlined a number of scenarios whereby black holes can be formed by domain walls or by vacuum bubbles nucleated during inflation. The resulting mass distribution of black holes depends on the microphysics parameters characterizing the  \nwalls and bubbles and on their interaction with matter. We will not attempt to explore all the possibilities here and will focus instead on one specific case: domain walls which interact very weakly with matter (so that matter particles can freely pass through the walls). A specific particle physics example could be axionic domain walls, whose couplings to the Standard Model particles are suppressed by the large Peccei-Quinn energy scale. The walls could also originate from a 'shadow' sector of the theory, which couples to the Standard Model only gravitationally.", 'A. Size distribution of domain walls': 'Spherical domain walls nucleate during inflation having radius R ≈ H -1 i ; then they are stretched by the expansion of the universe. At t -t n glyph[greatermuch] H -1 i , where t n is the nucleation time, the radius of the wall is well approximated by 10  \nR ( t ) ≈ H -1 i exp[ H i ( t -t n )] . (75)  \nThe wall nucleation rate per Hubble spacetime volume H -4 i is  \nλ = Ae -S . (76)  \nIn the thin-wall regime, when the wall thickness is small compared to the Hubble radius H -1 i , the tunneling action S is given by [4]  \nS = 2 π 2 σ/H 3 i . (77)  \nThis estimate is valid as long as S glyph[greatermuch] 1, or  \nσ > ∼ H 3 i . (78)  \nWith this assumption, the prefactor has been estimated in [19] as  \nA ∼ ( σ/H 3 i ) 2 . (79)  \nThe number of walls that form in a coordinate interval d 3 x and in a time interval dt n is  \ndN = λH 4 i e 3 H i t n d 3 x dt n . (80)  \nUsing Eq. (75) to express t n in terms of R , we find the distribution of domain wall radii [4],  \ndn ≡ dN dV = λ dR R 4 , (81)  \nwhere  \ndV = e 3 H i t d 3 x (82)  \nis the physical volume element. At the end of inflation, the distribution (81) spans the range of scales from R ∼ H -1 i to R max ∼ H -1 i exp( N inf ), where N inf is the number of inflationary e-foldings. For models of inflation that solve the horizon and flatness problems, the comoving size of R max is far greater than the present horizon.  \nAfter inflation, the walls are initially stretched by the expansion of the universe,  \nR ( t ) = a ( t ) a ( t i ) R i , (83)  \nwhere a ( t ) is the scale factor, t i corresponds to the end of inflation and R i = R ( t i ). The size distribution of the walls during this period is still given by Eq. (81).', 'B. Black hole mass distribution': 'As we discussed in Section III, the ultimate fate of a given domain wall depends on the relative magnitude of two time parameters: the Hubble crossing time t H when R ( t H ) = H -1 , where H = ˙ a/a is the Hubble parameter, and the time t σ ∼ H -1 σ = (2 πGσ ) -1 , when the repulsive gravity of the wall becomes dynamically important. Walls that cross the Hubble radius at t H glyph[lessmuch] t σ have little effect on the surrounding matter and can be treated as test walls in the FLRW background. The mass of such walls at Hubble crossing (disregarding their kinetic energy) is  \nM H ≈ 4 πσH -2 ∼ t H t σ M ( t H ) glyph[lessmuch] M ( t H ) , (84)  \nwhere  \nM ( t ) = 1 2 GH (85)  \nis the mass within a sphere of Hubble radius at time t . Once the wall comes within the Hubble radius, it collapses to a black hole of mass M = CM H with C ∼ 1 in about a Hubble time.  \nIn the opposite limit, t H glyph[greatermuch] t σ , the wall starts expanding faster than the background at t ∼ t σ and develops a wormhole, which is seen as a black hole from the FRW region. Assuming that the universe is dominated by nonrelativistic matter, we found that the mass of this black hole is  \nM ∼ M ( t H ) (86)  \nand its Schwarzschild radius is r g ∼ t H . The boundary between the two regimes, t H ∼ t σ , corresponds to black holes of mass M ∼ 1 /G 2 σ ∼ M cr , where M cr is the critical mass introduced in Section III. We note that depending on the wall tension σ , the critical mass M cr can take a wide range of values, including values of astrophysical interest. If we define the energy scale of the wall η as σ ∼ η 3 , then, as η varies from ∼ 1 GeV to the GUT scale ∼ 10 16 GeV , the critical mass varies from 10 17 M glyph[circledot] to 10 4 g .  \nThe estimate (86) applies for t H > ∼ t σ > ∼ t eq , where t eq is the time of equal matter and radiation densities. For walls that start developing wormholes during the radiation era, t σ < t eq , we could only obtain an upper bound,  \nM < ∼ M ( t H ) ( t σ < t eq ) . (87)  \nWe shall start with the case t σ > t eq , which can be fully described analytically. Note, however, that this case requires the energy scale of the wall to be η < ∼ 1 GeV , which is rather small by particle physics standards. The critical mass in this case is M cr > 10 17 M glyph[circledot] .', '1. t σ > t eq': 'Let us first consider black holes resulting from domain walls that collapse during the radiation era, t H < t eq , and have radii between R and R + dR at Hubble crossing. For such domain walls, the Hubble crossing time is t H = R/ 2, and their density at t > t H is  \ndn ( t ) ∼ λ ( R t ) 3 / 2 dR R 4 . (88)  \nThe mass distribution of black holes can now be found by expressing R in terms of M , R ∼ ( M/σ ) 1 / 2 , and substituting in Eq. (88):  \ndn = λ ( σ t 2 ) 3 / 4 dM M 7 / 4 . (89)  \nA useful characteristic of this distribution is the mass density of black holes per logarithmic mass interval in units of the dark matter density ρ DM ,  \nf ( M ) ≡ M 2 ρ DM dn dM . (90)  \nSince black hole and matter densities are diluted in the same way, f ( M ) remains constant in time. We shall evaluate it at the time of equal matter and radiation densities, t eq , when Eq. (88) derived for the radiation era should still apply by order of magnitude. With ρ DM ( t eq ) ∼ B -1 Gt 2 eq , where B ∼ 10, we have  \nf ( M ) ∼ λ ( σ ) 3 / 4 Gt 1 / 2 eq M 1 / 4 ∼ Bλ M 1 / 2 eq M 1 / 4 M 3 / 4 cr . (91)  \nHere, M eq ∼ t eq / 2 G ∼ 10 17 M glyph[circledot] is the dark matter mass within a Hubble radius at t eq .  \nFor black holes forming in the matter era, but before t σ ( t σ > t H > t eq ), Eqs. (88) and (91) are replaced by  \nf ( M ) ∼ Bλ ( M M cr ) 1 / 2 . (93)  \ndn = λ ( R t ) 2 dR R 4 , (92)  \nFinally, for black holes formed at t > t σ , the mass is  \nM ∼ M ( t H ) ∼ R/G. (94)  \nSubstituting this in Eq. (92), we find  \nf ( M ) ∼ Bλ. (95)  \nThe resulting mass distribution function is plotted in Fig. 8.', '2. t σ < t eq': 'We now turn to the more interesting case of t σ < t eq . (We shall see that observational constraints on our model parameters come only from this regime.) In this case, domain walls start developing wormholes during the radiation era, and we have only an upper bound (87) for the mass of the resulting black holes.  \nDomain walls that have radii R < t σ at horizon crossing collapse to black holes with M < M cr . The size distribution of such walls and the mass distribution of the black holes  \nFIG. 8: Mass distribution function for t σ > t eq .  \n<!-- image -->  \nGLYPH<159>  \nare still given by Eqs. (89) and (91), respectively. For t σ < R < t eq , let us assume for the moment that the bound (87) is saturated. Then, using R ∼ GM in Eq. (88), we obtain the mass distribution  \nf ( M ) ∼ Bλ ( M eq M ) 1 / 2 , (96)  \nwhere B ∼ 10, as before. This distribution applies for M cr < M < M eq . With the same assumption, for walls with R > t eq we find  \nf ( M ) ∼ Bλ. (97)  \nThe resulting mass distribution function is plotted in Fig. 9, with the parts depending on the assumed saturation of the mass bound shown by dashed lines.  \nThe assumption that the bound (87) is saturated for M > M cr appears to be a reasonable guess. We know that it is indeed saturated in a matter-dominated universe, and it yields a mass distribution that joins smoothly with the distribution we found for M < M cr . A reliable calculation of f ( M ) in this case should await numerical simulations of supercritical domain walls in a radiation-dominated universe. For the time being, the distribution we found here provides an upper bound for the black hole mass function.  \nFIG. 9: Mass distribution function for t σ < t eq .  \n<!-- image -->  \nGLYPH<159>', 'V. OBSERVATIONAL BOUNDS': 'Observational bounds on the mass spectrum of primordial black holes have been extensively studied; for an up to date review, see, e.g., [20]. For small black holes, the most stringent bound comes from the γ -ray background resulting from black hole evaporation:  \nf ( M ∼ 10 15 g ) < ∼ 10 -8 . (98)  \nFor massive black holes with M > 10 3 M glyph[circledot] , the strongest bound is due to distortions of the CMB spectrum produced by the radiation emitted by gas accreted onto the black holes [21]:  \nf ( M > 10 3 M glyph[circledot] ) < ∼ 10 -6 . (99)  \nOf course, the total mass density of black holes cannot exceed the density of the dark matter. Since the mass distribution in Fig. 9 is peaked at M bh = M cr , this implies  \nf ( M cr ) < 1 . (100)  \nThese bounds can now be used to impose constraints on the domain wall model that we analyzed in Section IV. As before, we shall proceed under the assumption that the mass bound (87) is saturated.  \nThe model is fully characterized by the parameters ξ = σ/H 3 i and H i /M pl , where H i is the expansion rate during inflation. The nucleation rate of domain walls λ depends only on  \nFIG. 10: Observational constraints on the distribution of black holes produced by domain walls. Red and purple regions mark the parameter values excluded, respectively, by small black hole evaporation and by gas accretion onto large black holes. The green region indicates the parameter values allowing the formation of superheavy black hole seeds.The solid straight line marks the parameter values where f ( M cr ) = 1, so the parameter space below this line is excluded by Eq. (100).  \n<!-- image -->  \nξ ,  \nλ ∼ ξ 2 e -2 π 2 ξ (101)  \n(see Eqs. (76), (77), (79)). The parameter space { ξ, H i /M pl } is shown in Fig. 10, with red and purple regions indicating parameter values excluded by the constraints (98) and (99), respectively. We show only the range ξ > ∼ 1, where the semiclassical tunneling calculation is justified. Also, for ξ glyph[greatermuch] 1 the nucleation rate λ is too small to be interesting, so we only show the values ξ ∼ few.  \nThe dotted lines in the figure indicate the values M cr ∼ 10 15 g ≡ M evap and M cr ∼ 10 3 M glyph[circledot] ≡ M accr . These lines, which mark the transitions between the subcritical and supercritical regimes in the excluded red and purple regions, are nearly vertical. This is because ξ ∼ 1 in the entire range shown in the figure, and therefore M cr ∼ ξ -1 ( M pl /H i ) 3 M pl depends  \nessentially only on H i .  \nThe horizontal upper boundaries of the red and purple regions correspond to the regime of M > M cr , where the mass function (96) depends only on λ . The corresponding constraints on λ can be easily found: λ < ∼ 10 -15 for the red region (evaporation bound) and λ < ∼ 10 -26 for the purple region (accretion bound). As we have emphasized, the mass function (96) that we obtained for super-critical domain walls represents an upper bound on the black hole mass function, and thus the red and purple regions in Fig. 10 may overestimate the actual size of the excluded part of the parameter space. We have verified that the regime t σ > t eq does not yield any additional constraints on model parameters.  \nA very interesting possibility is that the primordial black holes could serve as seeds for supermassive black holes observed at the galactic centers. The mass of such seed black holes should be M > 10 3 M glyph[circledot] [22] and their number density n 0 ( M ) at present should be comparable to the density of large galaxies, n G ∼ 0 . 1 Mpc -3 . The region of the parameter space satisfying this condition is marked in green in Fig. 10. (As before, the horizontal upper boundary of this region may overestimate the black hole mass function.)  \nThe solid straight line in Fig. 10 marks the parameter values where f ( M cr ) = 1, so the parameter space below this line is excluded by Eq. (100). This line lies completely within the red and purple excluded regions, suggesting that the bound (100) does not impose any additional constraints on the model parameters. We note, however, that the observationally excluded regions are plotted in Fig. 10 assuming that the expansion rate H i remains nearly constant during inflation. Small domain walls that form black holes with masses M evap ∼ 10 15 M glyph[circledot] are produced towards the very end of inflation, when the expansion rate may be substantially reduced. The tunneling action (77) depends on H i exponentially, so any decrease in H i may result in a strong suppression of small wall nucleation rate. The evaporation bound on the model parameters can then be evaded, while larger black holes that could form dark matter or seed the supermassive black holes may still be produced. For example, with ξ ≈ 2, a change in H i by 20% would suppress the nucleation rate by 10 orders of magnitude.  \nThe black holes that form dark matter in our scenario can have masses in the range 10 15 g < M bh < 10 24 g . For M bh < 10 15 g the black holes would have evaporated by now, and for M bh > 10 24 g the line f ( M cr ) = 1 in Fig. 10 enters the purple region excluded by the gas  \naccretion constraint. 11 We note that Ref. [23] suggested an interesting possibility that dark matter could consist of intermediate-mass black holes with 10 2 M glyph[circledot] < ∼ M bh < ∼ 10 5 M glyph[circledot] . Most of this range, however, appears to be excluded by the gas accretion constraint (99).  \nWe note also the recent paper by Pani and Loeb [26] suggesing that small black holes could be captured by neutron stars and could cause the stars to implode. They argue that this mechanism excludes black holes with masses 10 17 g < ∼ M bh < ∼ 10 23 g from being the dominant dark matter constituents. This conclusion, however, was questioned in Ref. [27].', 'VI. SUMMARY AND DISCUSSION': "We have explored the cosmological consequences of a population of spherical domain walls and vacuum bubbles which may have nucleated during the last N ∼ 60 e-foldings of inflation. At the time t i when inflation ends, the sizes of these objects would be distributed in the range t i < ∼ R i < ∼ e N t i , and the number of objects within our observable universe would be of the order N ∼ λe 3 N , where λ is the dimensionless nucleation rate per unit Hubble volume. We have shown that such walls and bubbles may result in the formation of black holes with a wide spectrum of masses. Black holes having mass above a certain critical value would have a nontrivial spacetime structure, with a baby universe connected by a wormhole to the exterior FRW region.  \nThe evolution of vacuum bubbles after inflation depends on two parameters, t b ∼ ( Gρ b ) -1 / 2 and t σ ∼ ( Gσ ) -1 , where ρ b > 0 is the vacuum energy inside the bubble and σ is the bubble wall tension. For simplicity, throughout this paper we have assumed the separation of scales t i glyph[lessmuch] t b , t σ . At the end of inflation, the energy of a bubble is equivalent to the 'excluded' mass of matter which would fit the volume occupied by the bubble, M i = (4 π/ 3) ρ m ( t i ) R 3 i , where ρ m is the matter density. This energy is mostly in the form of kinetic energy of the bubble walls, which are expanding into matter with a high Lorentz factor. Assuming that particles of matter cannot penetrate the bubble and are reflected from the bubble wall, this kinetic energy is quickly dissipated by momentum transfer to the surrounding matter on a time-scale much shorter than the Hubble time, so the wall comes to rest with respect to the Hubble flow. In this process, a highly relativistic and dense  \nexploding shell of matter is released, while the energy of the bubble is significantly reduced.  \nIf the bubble is surrounded by pressureless matter, it subsequently decouples from the Hubble flow and collapses to form a black hole of mass  \nM bh ∼ ( t 2 i t 2 b + t i t σ ) M i glyph[lessmuch] M i , (102)  \nwhich is sitting in the middle of an empty cavity. The fate of the bubble depends on whether its mass is larger or smaller than a certain critical mass, which can be estimated as  \nM cr ∼ Min { t b , t σ } /G. (103)  \nFor M bh < M cr , the bubble collapses into a Schwarzschild singularity, as in the collapse of usual matter. However, for M bh > M cr the bubble avoids the singularity and starts growing exponentially fast inside of a baby universe, which is connected to the parent FRW by a wormhole. This process is represented in the causal diagram of Fig. 3, of which a spatial slice is depicted in Fig. 4. The exponential growth of the bubble size can be due to the internal vacuum energy of the bubble, or due to the repulsive gravitational field of the bubble wall. The dominant effect depends on whether t b is smaller than t σ or vice-versa.  \nIf the bubble is surrounded by radiation, rather than pressureless dust, there is one more step to consider. After the bubble wall has transfered its momentum to the surounding matter and comes to rest with respect to the Hubble flow, the effects of pressure may still change the mass of the resulting black hole. We have argued that for M bh glyph[lessmuch] M cr the black hole mass, given by Eq. (44), is not significantly affected by the pressure of ambient matter. However, for M bh glyph[greatermuch] M cr this may be important, since some radiation may follow the bubble into the baby universe, thus affecting the value of the resulting black hole mass. This process seems to require a numerical study which is beyond the scope of the present work. On general grounds, however, we pointed out that the mass of the black hole is bounded by M bh < ∼ t H /G , where t H is the time when the co-moving region corresponding to the initial size R i crosses the horizon.  \nThe case of domain walls is somewhat simpler than the case of bubbles, since only the scale t σ is relevant. If the wall is small, with a horizon crossing time t H < t σ , and assuming that the only interaction with matter is gravitational, the wall will collapse to a black hole of mass M bh ∼ σt 2 H . Alternatively, if the wall is impermeable to matter, a pressure supported bag of matter will form, with a substantially larger mass which is given in Eqs. (51) or (70).  \nThese estimates apply respectively to the case where matter is non-relativistic or relativistic at the time when the wall crosses the horizon. Larger walls with t H > t σ will begin to inflate at t ∼ t σ , developing a wormhole structure. If at t ∼ t σ the universe is dominated by nonrelativistic matter, the mass of the resulting black hole is M bh ∼ t H /G , where t H is now the time when the comoving region affected by wall nucleation comes within the horizon. For t σ in the radiation era, the estimate of the mass is more difficult to obtain, since some radiation may flow into the baby universe following the domain wall. This issue requires a numerical study and is left for further research. As in the case of bubbles, for this case we were able to find only an upper bound on the mass, M bh < ∼ t H /G .  \nAsystematic analysis of all possible scenarios would require numerical simulations, and we have not attempted to do it in this paper. Nonetheless, for illustration, we have considered the case of a distribution of domain walls which interact with matter only gravitationally. Also, awaiting confirmation from a more detailed numerical study, we have assumed that large walls entering the horizon in the radiation era, with t σ < ∼ t H < ∼ t eq , lead to black holes which saturate the upper bound on the mass, M bh ∼ t H /G . Under these assumptions, we have shown that black holes produced by nucleated domain walls can have a significant impact on cosmology. Hawking evaporation of small black holes can produce a γ -ray background, and radiation emitted by gas accreted onto large black holes can induce distortions in the CMB spectrum. A substantial portion of the parameter space of the model is already excluded by present observational constraints, but there is a range of parameter values that would yield large black holes in numbers sufficient to seed the supermassive black holes observed at the centers of galaxies. For certain parameter values the black holes can also play the role of dark matter.  \nOur preliminary analysis provides a strong incentive to consider the scenario of black hole formation by vacuum bubbles, which has a larger parameter space. Note that the black hole density resulting from domain wall nucleation is appreciable only if ξ = σ/H 3 i ∼ 1, where H i is the expansion rate during inflation. This calls for a rather special choice of model parameters, with 12 σ ∼ H 3 i . The narrowness of this range is related to the fact that the nucleation rate of domain walls is exponentially sensitive to the wall tension.  \nFIG. 11: The eternally inflating multiverse generally has a very nontrivial spacetime structure, with a multitude of eternally inflating regions connected by wormholes. This is reminiscent of a well known illustration of the inflationary multiverse by Andrei Linde, which we reproduce here by courtesy of the author (see also [28] and references therein). In the present context the links between the balloon shaped regions should be interpreted as wormholes which are created by supercritical bubbles or domain walls.  \n<!-- image -->  \nOn the other hand, the bubble nucleation rate is highly model-dependent, and the parameter space yielding an appreciable density of black holes is likely to be increased. We note also that string theory suggests the existence of a vast landscape of vacua [25], so one can expect a large number of bubble types, some of them with relatively high nucleation rates.  \nSuper-critical black holes produced by vacuum bubbles have inflating vacuum regions inside. These regions will inflate eternally and will inevitably nucleate bubbles of the highenergy parent vacuum with the expansion rate H i , as well as bubbles of all other vacua allowed by the underlying particle physics [2, 3, 12, 28]. In some of the bubbles inflation will end, and bubbles nucleated prior to the end of inflation will form black holes with inflating universes inside. A similar structure would arise within the supercritical black holes produced by domain walls. Indeed, inflating domain walls can excite any fields coupled to  \nthem, even if the coupling is only gravitational. In particular, they can excite the inflaton field, causing it to jump back to the inflationary plateau [29, 30], or they may cause transitions to any of the accessible neighboring vacua [31]. We thus conclude that the eternally inflating multiverse generally has a very nontrivial spacetime structure, with a multitude of eternally inflating regions connected by wormholes. 13 This is reminiscent of a well known illustration of the inflationary multiverse, by Andrei Linde which we reproduce in Fig. 11. In the present context the links between the balloon shaped regions should be interpreted as wormholes. We note that the mass distributions of black holes resulting from domain walls and from vacuum bubbles are expected to be different and can in principle be distinguished observationally.  \nBlack hole formation in the very early universe has been extensively studied in the literature, and a number of possible scenarios have been proposed (for a review see [18, 32]). Black holes could be formed at cosmological phase transitions or at the end of inflation. All these scenarios typically give a sharply peaked mass spectrum of black holes, with the characteristic mass comparable to the horizon mass at the relevant epoch. In contrast, our scenario predicts a black hole distribution with a very wide spectrum of masses, scanning many orders of magnitude. 14 If a black hole population with the predicted mass spectrum is discovered, it could be regarded as evidence for inflation and for the existence of a multiverse.  \nWe would like to conclude with two comments concerning the global structure of supercritical black holes, which may be relevant to the so-called black hole information paradox (see Ref. [35] for a recent review). First of all, given that the black hole has a finite mass, the amount of information carried by Hawking radiation as the black hole evaporates into the asymptotic FRW universe will be finite. For that reason, it cannot possibly encode the quantum state of the infinite multiverse which develops beyond the black hole horizon. Second, the supercritical collapse leads to topology change, and to the formation of a baby universe which survives arbitrarily far into the future. In this context, it appears that unitary evolution should not be expected in the parent universe, once we trace over the degrees of freedom in the baby universes. Finally, we should stress that the formation of wormholes,  \nleading to topology change and to formation of baby multiverses within black holes, does not require any exotic physics. It occurs by completely natural causes from regular initial conditions. 15", 'VII. ACKNOWLEDGEMENTS': 'This work is partially supported by MEC FPA2013-46570-C2-2-P, AGAUR 2014-SGR1474, CPAN CSD2007-00042 Consolider-Ingenio 2010 (JG) and by the National Science Foundation (AV and JZ). J.Z. was also supported by the Burlingame Fellowship at Tufts University. A.V. is grateful to Tsutomu Yanagida for a useful discussion, suggesting the possible relevance of our scenario for dark matter and to Abi Loeb for several useful discussions.', 'Appendix A: Evolution of test domain walls in FRW': "In this appendix, we solve the evolution of test domain walls in a FRW universe numerically and find the numerical coefficient C r and C m defined in Sec. III.  \nThe metric for a FRW universe is  \nds 2 = a 2 ( η )( -dη 2 + dr 2 + r 2 d Ω 2 ) , (A1)  \nwhere η is the conformal time. The scale factor is given by  \na ( η ) = ( η η i ) p , (A2)  \nwhere p = 2 for the dust dominated background and p = 1 for the radiation dominated background. We choose η i = p , so that H i = 1.  \nThe domain wall worldsheet is parametrized by η , and the worldsheet metric can be written as  \nds 2 3 = a 2 ( η ) [( -1 + r ' 2 ) dη 2 + r 2 d Ω 2 ] . (A3)  \nIn this Appendix, prime indicates a derivative with respect to the conformal time η . The domain wall action is proportional to the worldsheet area,  \nS = 4 πσ ∫ dτ √ 1 -r ' 2 a 3 r 2 , (A4)  \nand the equation of motion is  \nr '' +3 a ' a r ' (1 -r ' 2 ) + 2 r (1 -r ' 2 ) = 0 . (A5)  \nWe solve this equation with the initial conditions  \nr ( η i ) = R i , r ' ( η i ) = 0 . (A6)  \nThe solutions are plotted in Fig. 5 in terms of the cosmological time t .  \nWe regard a black hole as formed when the radial coordinate drops to r < 10 -4 t i and determine its mass from Eq. (53), which in terms of the conformal time η takes the form  \nM ( η ) = 4 πσ a 2 r 2 √ 1 -r ' 2 . (A7)  \nThis mass grows while the wall is being stretched by Hubble expansion, but remains nearly constant during the late stages of the collapse.  \nAccording to Sec. III, the numerical coefficients C m and C r are defined as  \nC m, r = p 2 (1 + p ) 2 M bh 4 πσt 2 H , (A8)  \nwhere the Hubble-crossing time t H satisfies  \nR i a ( t H ) = 1 + p p t H , (A9)  \nThe dependence of C m,r on H i R i = R i is shown in Fig. 12. We see that the coefficients approach the values C m glyph[similarequal] 0 . 15 and C r glyph[similarequal] 0 . 62 for H i R i > 10.  \nFIG. 12: The dependence of C m and C r on H i R i = R i . For H i R i > 10, we find C m glyph[similarequal] 0 . 15 and C r glyph[similarequal] 0 . 62.  \n<!-- image -->", 'Appendix B: Large domain wall in dust cosmology': 'In this Appendix we describe in more detail the construction of the solution represented in Fig. 6, which corresponds to a large domain wall embedded in a dust cosmology. The gravitational field of the wall creates a wormhole structure leading to a baby universe. We shall start by cosidering the matching of Schwarzschild to a dust cosmology, and then the matching of two Schwarzschild metrics accross a domain wall. Finally, the different pieces are put together.', '1. Matching Schwarzschild to a dust cosmology': 'Consider the Schwarzschild metric  \nds 2 = -( 1 -2 GM R ) dT 2 + ( 1 -2 GM R ) -1 dR 2 + R 2 d Ω 2 . (B1)  \nIn Lemaitre coordinates, this takes the form  \nds 2 = -d ˆ t 2 + 2 GM R dρ 2 + R 2 d Ω 2 . (B2)  \nwhere ˆ t and ρ are defined by the relations  \nd ˆ t = ± dT -√ 2 GM R ( 1 -2 GM R ) -1 dR, (B3)  \ndρ = ∓ dT + √ R 2 GM ( 1 -2 GM R ) -1 dR. (B4)  \nAdding (B3) and (B4), we have  \nR = [ 3 2 ( ˆ t + ρ ) ] 2 / 3 (2 GM ) 1 / 3 , (B5)  \nwhere an integration constant has been absorbed by a shift in the origin of the ρ coordinate.  \nThe expression for T as a function of ˆ t and ρ can be found from (B4) as:  \n± T = 4 GM ( 1 3 ( R 2 GM ) 3 / 2 + √ R 2 GM -tanh -1 √ R 2 GM ) -ρ ± , (B6)  \nwith R given in (B5). A second integration constant has been absorbed by a shift in the T coordinate. Since the metric (B2) is synchronous, the lines of constant spatial coordinate are geodesics, and ˆ t is the proper time along them. The double sign in front of T corresponds to two possible choices of geodesic coordinates ρ ± . We shall consider geodesics which start from the singularity at R = 0 at the bottom of the Kruskal diagram, see Fig. 13. In the region R < 2 GM the T coordinate may grow or decrease depending on which spatial direction we go, while ˆ t and R both grow towards the future. A sphere of test particles at ρ = ρ 0 has the physical radius given by  \nR ∝ ( ˆ t + ρ 0 ) 2 / 3 . (B7)  \nThus, co-moving spheres with different values of the ρ coordinate have a similar behaviour, with a shifted value of proper time.  \nThe metric (B2) can be matched to a flat FRW metric  \nds 2 = -dt 2 + a 2 ( t )( dr 2 + r 2 d Ω 2 ) . (B8)  \nThe co-moving sphere with ρ = ρ 0 will be glued to a co-moving sphere with r = r 0 . Continuity of the temporal and angular components of the metric accross that surface requires  \nR = a ( t ) r 0 , (B9)  \nwith t = ˆ t + ρ 0 , and  \na ( t ) ∝ t 2 / 3 , (B10)  \ncorresponding to a dust cosmology. Such cosmology has a matter density given by ρ m = 1 / (6 πGt 2 ). Using (B5) with ρ = ρ 0 , it is then straightforward to show that  \nM = 4 π 3 ρ m R 3 . (B11)  \ne  \nr  \ne  \n)  \n,  \nτ  \nn  \nρ  \nr  \ni  \nn  \ne  \na  \nm  \na  \nC  \nt  \ne  \ne  \nr  \ne  \nd  \nd  \nρ  \n=  \n∓  \nd  \ny  \n1  \no  \ni  \n.  \nn  \nc  \nfi  \nn  \nd  \ne  \nb  \nt  \no  \nd  \na  \ns  \nd  \nM  \nr  \nt  \na  \nh  \na  \np  \nh  \ne  \nB  \nn  \na  \ne  \nf  \nd  \no  \nu  \nt  \ne  \nt  \na  \n1  \ns  \ne  \nd  \no  \ni  \nr  \nn  \ne  \nh  \n=  \n±  \ndT  \nρ  \n3  \na  \nw  \nb  \ne  \nd  \ne  \nρ  \nn  \nc  \n3  \n/  \n2  \n+  \n2  \ne  \na  \nv  \nI  \nL  \nd  \nh  \ne  \nR  \n[  \na  \na  \nb  \nn  \ns  \no  \nb  \nr  \ne  \nb  \ne  \nf  \no  \n(  \nτ  \n+  \nd  \nu  \n√  \nh  \na  \ns  \nn  \nt  \ne  \nGM  \nb  \ns  \np  \ne  \na  \nn  \nt  \ni  \na  \na  \nb  \nl  \ne  \nr  \ndT  \n+  \n2  \n√  \n2  \n/  \n3  \n(  \n2  \ni  \n]  \nb  \nd  \nR  \n-  \ns  \nc  \no  \no  \nr  \no  \nb  \nr  \ne  \nd  \n)  \ny  \nn  \nf  \na  \nr  \no  \nt  \ns  \nm  \na  \nn  \nh  \nf  \n(  \nh  \nd  \ni  \nn  \nb  \na  \ny  \nt  \ne  \na  \na  \nGM  \nt  \n3  \n1  \n-  \n2  \ni  \nn  \na  \n)  \ns  \n:  \n√  \nR  \nGM  \ns  \nr  \nh  \ne  \ni  \nf  \ng  \nt  \ne  \no  \nd  \n=  \n2  \ns  \nc  \nm  \nl  \nl  \ni  \nn  \ndus  \nt  \nco  \no  \ns  \nm  \nl  \nn  \ni  \np  \na  \ns  \na  \nt  \na  \nl  \nl  \ny  \nfl  \nm  \na  \nt  \nn  \na  \ns  \nt  \nh  \ne  \ni  \nt  \na  \nt  \ni  \ni  \nv  \ne  \na  \nn  \na  \nl  \nT  \n.  \nr  \ni  \nv  \ne  \no  \ne  \ns  \nt  \n/lessmuch  \nR  \nt  \n-  \nt  \nσ  \n,  \n1  \n,  \ns  \n/lessmuch  \nH  \nl  \na  \nl  \nf  \na  \nw  \n∼  \nh  \nr  \nc  \nr  \nO  \ne  \n.  \nv  \na  \nl  \nl  \ni  \ns  \ng  \nu  \ne  \nn  \nd  \ns  \ne  \nc  \na  \ns  \n,  \nf  \no  \nl  \ne  \nn  \nt  \na  \ne  \ns  \nw  \nh  \nl  \n,  \ni  \ni  \nh  \na  \nc  \nh  \ne  \nw  \np  \nu  \nl  \nt  \nt  \ne  \nn  \nr  \ni  \na  \nn  \nI  \nm  \n.  \nc  \na  \nl  \ny  \ng  \nt  \no  \nv  \nu  \nl  \no  \na  \nn  \nu  \nh  \no  \ne  \nd  \no  \nf  \ns  \nt  \ni  \no  \nt  \no  \nr  \ng  \nn  \n2  \nGM  \nn  \nw  \na  \ng  \n=  \no  \nt  \n,  \nT  \n→-  \n∞  \nb  \ne  \nd  \nd  \ne  \nd  \nc  \ni  \ns  \ne  \ne  \ni  \nd  \nt  \na  \nu  \na  \nm  \ne  \nn  \nl  \nm  \nl  \ne  \ny  \ne  \nd  \ns  \nw  \na  \nl  \nt  \no  \na  \nh  \nc  \na  \nr  \nd  \nm  \nc  \ns  \n.  \nE  \ne  \no  \no  \no  \nn  \ne  \nh  \nh  \nl  \ne  \ne  \nl  \nk  \na  \nr  \nl  \nf  \nr  \no  \nl  \no  \nd  \nf  \nd  \ni  \na  \ni  \nu  \ns  \ne  \nd  \nh  \na  \no  \nm  \nS  \n=  \nl  \nl  \nn  \nw  \ni  \ni  \nt  \nh  \nl  \nt  \n2  \n/t  \nσ  \n/lessmuch  \nt  \n.  \nt  \nh  \ne  \nh  \ns  \na  \nl  \nl  \nl  \nl  \ng  \na  \nt  \nw  \nn  \no  \ne  \ns  \nd  \ni  \nn  \ns  \nz  \nw  \nS  \nr  \na  \ny  \nt  \nc  \no  \nn  \np  \ne  \nc  \na  \ns  \n,  \nw  \ns  \na  \nl  \ni  \ns  \ns  \ni  \nt  \nh  \ne  \n∼  \nw  \ni  \n2  \nGM  \nI  \nσ  \nt  \n.  \nn  \nr  \no  \nn  \n.  \na  \na  \nn  \ns  \ni  \na  \nR  \n/greatermuch  \nx  \nR  \ne  \nr  \ne  \nt  \nt  \nw  \nh  \ne  \nn  \ni  \ne  \nt  \nl  \nm  \ni  \nt  \n,  \nm  \ne  \ne  \n2  \na  \nb  \n,  \nu  \ne  \nn  \na  \np  \nd  \no  \ne  \ni  \ng  \ni  \nn  \nm  \nl  \nl  \ny  \nw  \ni  \nn  \nt  \ni  \ne  \nt  \nA  \nl  \nr  \nl  \ne  \ni  \nn  \nt  \ni  \no  \nr  \no  \nf  \ne  \ns  \ne  \nn  \nt  \nh  \nt  \ne  \nr  \ne  \nm  \na  \ni  \nh  \nl  \ne  \nl  \ni  \nl  \no  \ne  \na  \nl  \nc  \nm  \ny  \ns  \ne  \nh  \na  \nc  \nf  \ng  \ne  \nr  \na  \nt  \ne  \nz  \no  \n,  \nh  \nn  \ni  \nk  \nl  \na  \nc  \ns  \nc  \nh  \no  \nr  \na  \ne  \ni  \nh  \nh  \nt  \nh  \n.  \nw  \ns  \ny  \nt  \nfi  \nh  \no  \nr  \nT  \na  \nt  \ni  \nf  \nt  \nh  \nw  \nn  \nw  \nh  \nl  \nl  \nt  \n,  \ns  \ne  \ne  \ns  \ns  \nt  \nh  \ni  \nt  \no  \ni  \nn  \ne  \no  \ns  \ni  \n(  \n1  \n)  \nm  \no  \ni  \ne  \nn  \na  \nt  \nd  \nu  \nn  \nr  \ni  \nv  \ne  \ns  \ne  \nT  \n.  \nw  \no  \n1  \nT  \n,  \no  \n1  \nw  \nH  \ne  \ns  \n.  \nr  \ni  \ne  \ne  \nt  \n-  \nw  \nH  \nl  \n∼  \nn  \nl  \ns  \nd  \nl  \ni  \nt  \nc  \na  \na  \ng  \nv  \nu  \nc  \ne  \nn  \no  \nm  \nm  \ni  \nl  \ns  \n2  \no  \no  \nr  \na  \nT  \nw  \no  \nm  \na  \nl  \no  \ne  \ni  \nl  \nf  \ne  \n)  \nd  \ne  \n.  \nr  \n(  \ns  \ne  \ni  \nv  \nu  \nl  \ne  \nn  \nfi  \nh  \nn  \na  \ne  \nt  \nm  \nc  \no  \nl  \nd  \no  \ns  \no  \nt  \nl  \no  \ne  \nc  \na  \no  \nn  \na  \nh  \ng  \nd  \nc  \nc  \ni  \nn  \ns  \na  \n,  \na  \n1  \nf  \n-  \nd  \no  \nH  \nl  \n∼  \ne  \nl  \nfi  \nd  \ns  \nt  \no  \nt  \nd  \nd  \nc  \na  \ni  \nl  \ni  \nv  \ni  \nt  \na  \no  \ng  \na  \nfi  \nt  \nt  \na  \nt  \ne  \na  \nl  \n-  \nl  \n,  \nl  \na  \nl  \na  \nn  \nw  \ni  \n∼  \nt  \na  \no  \ne  \nd  \nh  \na  \nl  \nl  \no  \nw  \ni  \n,  \nn  \nt  \nm  \nσ  \nt  \no  \ne  \nR  \n/lessmuch  \nh  \nf  \nd  \no  \nl  \ne  \n/lessmuch  \nf  \na  \na  \nl  \nl  \nh  \ne  \nd  \nm  \na  \nm  \ni  \na  \ns  \nl  \nl  \nt  \nn  \nw  \ni  \nn  \no  \ne  \ni  \nh  \nr  \nz  \no  \ni  \na  \no  \n,  \nm  \ns  \ne  \nd  \nh  \nt  \ni  \ni  \nt  \nh  \nn  \nh  \ni  \n.  \nt  \nh  \nw  \ni  \nn  \no  \nr  \np  \nh  \nr  \nl  \ne  \ne  \nl  \nt  \nw  \na  \ne  \nn  \nt  \nσ  \n,  \nn  \nt  \na  \n)  \ns  \nw  \n2  \n3  \nl  \n/t  \nd  \nR  \n(  \nl  \no  \ne  \n∼  \nt  \n/lessmuch  \nm  \nl  \nw  \ni  \na  \ni  \nl  \nl  \nt  \nh  \n,  \nz  \nt  \n,  \ni  \nt  \nn  \ni  \nσ  \n/lessmuch  \nw  \nl  \nl  \ns  \ns  \nf  \na  \nl  \nI  \n/lessmuch  \nw  \no  \ne  \nt  \nh  \n)  \nn  \n,  \na  \n1  \ne  \ns  \nf  \nr  \nGM  \nt  \nt  \nh  \n,  \n2  \nσ  \nn  \nt  \nS  \nR  \nR  \ns  \nu  \nl  \n=  \nw  \na  \n.  \nw  \nt  \nt  \nh  \ne  \nt  \nt  \ni  \ni  \ns  \nh  \nh  \ne  \nc  \nh  \nn  \nσ  \nn  \nt  \n2  \ni  \nr  \no  \ns  \nt  \na  \nt  \n2  \n/t  \no  \ni  \ne  \nh  \nh  \ne  \n(  \nr  \nGM  \no  \ns  \ni  \n=  \na  \n/greatermuch  \ne  \nl  \nf  \nr  \nh  \n/lessmuch  \ne  \ne  \nR  \nl  \nw  \na  \nh  \ny  \n,  \ne  \nw  \na  \nl  \nl  \ns  \nd  \nn  \nt  \ni  \ni  \nu  \ne  \nr  \n/lessmuch  \nb  \na  \nw  \nf  \ne  \nS  \nx  \np  \nR  \nGM  \nn  \n2  \n/t  \nn  \n)  \nA  \n.  \n4  \n∼  \na  \nσ  \n.  \n/lessmuch  \nw  \ne  \ns  \n.  \nl  \nl  \ns  \na  \ne  \ne  \nt  \nl  \nl  \ng  \na  \ns  \nh  \ne  \n,  \np  \nσ  \ne  \nr  \nc  \na  \nl  \ns  \ns  \nt  \nI  \nt  \nn  \nt  \nr  \n.  \nt  \ne  \nσ  \ni  \nn  \nt  \n/lessmuch  \ne  \ni  \nh  \nt  \nt  \ni  \nn  \nr  \nl  \n∼  \na  \ne  \nl  \ny  \ne  \nR  \nn  \n.  \n/greatermuch  \ns  \nt  \nI  \nh  \ni  \nm  \na  \no  \nS  \n=  \n2  \nr  \nl  \ns  \n,  \no  \nm  \ne  \na  \nw  \ni  \ni  \ni  \ne  \nw  \ng  \nt  \nR  \ne  \nn  \nt  \ni  \no  \n,  \na  \nh  \nr  \ni  \nl  \ne  \nl  \no  \nf  \nt  \nh  \ne  \ns  \nw  \nh  \n,  \nn  \ns  \nt  \ne  \nA  \ne  \nw  \nn  \na  \ne  \ne  \nr  \nh  \ne  \nc  \nh  \n.  \ns  \nr  \ne  \nt  \nl  \nl  \ns  \nw  \nh  \ns  \nc  \na  \no  \nn  \nk  \ns  \nl  \ni  \na  \na  \nn  \n.  \nb  \np  \nt  \na  \nt  \n)  \nt  \ni  \n(  \nσ  \nf  \n2  \nt  \nw  \ni  \nh  \ne  \nx  \na  \nn  \ne  \nn  \no  \nm  \nl  \nh  \ne  \ns  \ne  \nf  \nt  \nA  \ne  \nr  \no  \nn  \na  \nl  \ni  \nn  \ni  \nl  \ny  \nt  \ne  \ni  \nb  \nt  \ni  \na  \nr  \ne  \n/greatermuch  \nl  \ne  \nm  \ne  \na  \nn  \nx  \na  \nn  \n<!-- image -->  \nng  \nt  \nSphe  \nc  \nk  \ne  \ne  \n-  \ns  \nb  \nu  \nh  \nt  \ne  \n=  \nl  \n2  \ne  \nc  \ns  \na  \na  \ndo  \ni  \nn  \na  \nSphe  \nl  \n1  \nh  \nm  \nb  \ni  \na  \nu  \n.  \nT  \ne  \no  \nl  \ng  \ne  \ne  \no  \nd  \nw  \ni  \nr  \nr  \nc  \ni  \na  \n→∞  \ne  \nn  \nm  \nl  \nt  \ne  \ne  \ns  \nr  \nn  \no  \ne  \nn  \ndifferent time-scales are relevant to the dynamics of such a wall. One is the cosmological scale  \nt  \nt  \nl  \ne  \nds  \n,  \nSphe  \ni  \nl  \nca  \nl  \nt  \n2  \n,  \ni  \nn  \nm  \nGM  \nt  \na  \ns  \n,  \nt  \nh  \ni  \ns  \na  \np  \na  \na  \ne  \nr  \nt  \ns  \nh  \ne  \nT  \n→∞  \n-  \nr  \ni  \nc  \nr  \nds  \ne  \nr  \nm  \nT  \ne  \nr  \na  \nt  \nv  \no  \n,  \na  \nc  \nl  \ns  \n,  \no  \nk  \no  \ne  \nn  \n=  \na  \nf  \nc  \ne  \ns  \nr  \na  \nl  \nc  \nd  \na  \nc  \n→-  \ne  \nr  \na  \nl  \nn  \nb  \ne  \ns  \nn  \ni  \nm  \nl  \ne  \ne  \ni  \nl  \na  \nt  \no  \nt  \no  \nn  \nm  \ne  \ne  \nfi  \nl  \nd  \ng  \ni  \nn  \ny  \na  \n(  \nc  \n2  \ne  \nl  \nR  \nGM  \no  \nd  \nu  \na  \ni  \nv  \nt  \ny  \ne  \nl  \ne  \ng  \no  \nc  \na  \ne  \nS  \nl  \ni  \nk  \n2  \ny  \nc  \ny  \nt  \nh  \ne  \ni  \nGM  \n2  \ne  \nr  \ne  \nw  \nt  \nw  \nr  \nr  \ne  \nr  \na  \nl  \ne  \nd  \nh  \np  \ne  \ne  \na  \nt  \nr  \nr  \na  \n,  \ne  \nt  \nd  \nl  \na  \ne  \no  \nn  \nu  \nn  \ns  \n1  \nl  \nm  \nc  \na  \n-  \ne  \nc  \no  \ndT  \ndR  \no  \ns  \nl  \ns  \nr  \nd  \ny  \nn  \nb  \ns  \ne  \ne  \nf  \nx  \nw  \ne  \nr  \nw  \ne  \nH  \nR  \nt  \nt  \ne  \na  \nu  \ns  \ni  \nn  \np  \nI  \nr  \n.  \nR  \n2  \nd  \nd  \nf  \ni  \nc  \nn  \nn  \no  \ne  \nd  \n1  \n-  \ne  \na  \nd  \n+  \na  \no  \n-  \n√  \n2  \nl  \ni  \ne  \ns  \nfi  \na  \nt  \ne  \ni  \nn  \nn  \nb  \na  \nr  \n,  \nr  \nf  \nr  \ne  \nc  \nl  \no  \ni  \nn  \nd  \ns  \nx  \nt  \nr  \na  \ne  \no  \nt  \no  \ni  \ng  \nf  \nf  \ne  \nl  \ni  \nr  \ng  \nn  \nw  \nr  \n2  \nT  \no  \no  \no  \nR  \nn  \nt  \nv  \ne  \no  \n(  \nw  \nE  \nn  \ne  \nc  \ns  \nm  \nf  \n.  \n2  \nGM  \ng  \nd  \no  \nt  \ne  \nc  \nr  \ni  \nt  \nr  \ns  \nT  \nt  \ne  \n.  \n∼  \ns  \na  \nf  \ne  \nf  \na  \nr  \no  \nl  \nt  \nd  \nt  \ni  \nT  \nh  \nr  \ns  \no  \nw  \nh  \nu  \ne  \no  \nk  \no  \nh  \nc  \nS  \np  \ni  \ni  \ne  \nl  \nt  \n,  \nr  \ni  \nv  \n1  \ns  \ne  \nr  \ns  \nu  \nn  \nt  \nm  \ne  \ne  \nh  \nt  \na  \nb  \n-  \nt  \na  \nh  \nW  \nd  \nn  \ns  \ni  \nt  \n1  \nf  \na  \nt  \ni  \nn  \ne  \nd  \nh  \nt  \ni  \n,  \nn  \nh  \ne  \nm  \ne  \nc  \nr  \nc  \no  \nn  \nh  \no  \n(  \nr  \nzs  \n2  \n/  \n3  \nw  \ne  \ns  \np  \no  \nr  \ni  \ng  \ni  \nh  \nr  \nr  \na  \na  \ne  \nt  \ns  \no  \ng  \no  \na  \ni  \nt  \ne  \ni  \nn  \nw  \nh  \no  \nt  \ne  \nr  \nr  \no  \nl  \nl  \ni  \nt  \n(  \n5  \n)  \na  \nh  \nc  \nn  \no  \nk  \nl  \nc  \nm  \nh  \na  \na  \ne  \ns  \nl  \nc  \nn  \nm  \ng  \ny  \nt  \nn  \ne  \nh  \ne  \nh  \ns  \no  \nt  \na  \nr  \nl  \ne  \nh  \ns  \nc  \nu  \nt  \nd  \ni  \nn  \nh  \ne  \nl  \ns  \nt  \ne  \nfi  \nr  \ns  \ny  \na  \nt  \ni  \nl  \nl  \ni  \nc  \nh  \nl  \ni  \nn  \ni  \na  \ne  \nl  \n,  \ns  \nn  \nt  \ne  \nw  \ns  \nl  \ng  \ne  \n1  \n,  \nw  \n)  \nu  \ns  \nr  \ns  \n,  \no  \nu  \ns  \nn  \nw  \na  \nt  \nr  \ni  \nv  \no  \ni  \ni  \ns  \ns  \ne  \n,  \nc  \n2  \nd  \nR = [ 2 ( τ + ρ ) ] (2 GM ) , (5) where an integration constant has been absorbed by a shift in the origin of the ρ coordinate. The w h e r e a n p r e s s i o n f o ( 5 ) . A s e c ) i s s y n c h r K αβ = 1 2 ∂ n g αβ (B12)  \nng  \nc  \nr  \ns  \nt  \nt  \nh  \nw  \na  \ne  \nr  \ns  \n2  \nl  \ne  \nt  \n2  \nr  \nh  \ne  \nt  \nu  \na  \ne  \ns  \ni  \nt  \no  \nk  \nl  \na  \nt  \ng  \ni  \ni  \ni  \nn  \nf  \nw  \nc  \nb  \na  \ne  \nd  \no  \nρ  \nh  \nf  \nt  \nc  \ni  \no  \nr  \ne  \nn  \nb  \nt  \nT  \n.  \nh  \ne  \nf  \ne  \ng  \na  \nt  \ni  \ns  \nt  \ns  \no  \ns  \nh  \nn  \nc  \nc  \nm  \nd  \ne  \n(  \nu  \ni  \nt  \nt  \nfi  \nr  \ns  \no  \nt  \nu  \ns  \ns  \nh  \ns  \nc  \nu  \ni  \nl  \ne  \nr  \nl  \nm  \no  \nl  \n,  \nt  \nd  \n1  \nh  \n-  \no  \ng  \ny  \nR  \n2  \nh  \na  \n3  \n)  \n,  \n(  \n4  \n)  \n1  \nn  \nt  \ni  \no  \n-  \nt  \ns  \nρ  \ng  \n-  \nGM  \no  \ns  \no  \nm  \n√  \ny  \nb  \nco  \ns  \nr  \nb  \n2  \nt  \no  \ng  \na  \nl  \no  \nn  \ndus  \no  \ne  \ne  \n2  \na  \nGM  \nl  \ny  \nd  \nb  \ne  \n1  \n+  \n(  \nl  \nt  \ni  \ndT  \np  \na  \n)  \nm  \nm  \ns  \nr  \nn  \nt  \no  \nf  \n2  \nGM  \n2  \nd  \n)  \n-  \nr  \ni  \no  \ny  \nf  \ni  \ne  \ne  \ns  \nh  \nr  \ne  \nn  \nR  \nt  \na  \n1  \ng  \n2  \n1  \n-  \n+  \n2  \nGM  \na  \n)  \ndT  \n2  \n(  \no  \n2  \nGM  \n-  \n1  \n+  \n(  \n2  \nds  \na  \n)  \n=  \n2  \nd  \nτ  \n+  \nR  \nd  \nR  \nρ  \nh  \nr  \nm  \n-  \nf  \no  \no  \no  \nn  \nR  \ne  \na  \nt  \ni  \nl  \nh  \ne  \ndT  \n-  \n-  \n1  \nR  \n(  \n)  \n3  \n,  \nand the other is the acceleration time-scale due to the repulsive gravitational field of the domain wall t σ ∼ ( Gσ ) -1 . For R /lessmuch t σ , the repulsive field can be ignored. In this case, for t /lessmuch R /lessmuch t σ , the domain wall is conformally stretched by cosmological expansion. Eventually, when the wall falls within the horizon, it quickly shrinks under its tension and forms a black hole of radius R S = 2 GM ∼ t 2 /t σ /lessmuch t . Here, we will be primarily interested in the opposite limit, where R /greatermuch t σ . In this case, the wall repels the matter around it while its size grows faster than the ambient expansion rate. As we shall see, this leads to the formation of a wormhole. The dust which was originally in the interior of the wall goes into a baby universe, and in the ambient FRW universe we are left with a black hole remnant cysted in a spherical region of of vacuum. Before we consider the effect of the domain wall, let us first discuss the matching of Schwarzschild and a dust cosmology. 1.1 Matching Schwarzschild to a dust cosmology Consider the Schwarzschild metric ds 2 = -( 1 -2 GM R ) dT 2 + ( 1 -2 GM R ) -1 dR 2 + R 2 d Ω 2 . (1) In Lemaitre coordinates, this takes the form ds 2 = -dτ 2 + 2 GM R dρ 2 + R 2 d Ω 2 . (2) where τ and ρ are defined by the relations dτ = ± dT -√ 2 GM R ( 1 -2 GM R ) -1 dR, (3) dρ = ∓ dT + √ R 2 GM ( 1 -2 GM R ) -1 dR. (4) Subtracting ( ?? ) from ( ?? ), we have 3 2 / 3 1 / 3 a n d 1 . 1 M a C o n s i d e r t h e S c I n L e m a i t r e c o o r d i w h e r e τ a n d ρ a r e d e fi n e d b y t h d τ = ± d ρ = ∓ dT + √ u b t r a c t i n g ( 3 ) f r o m ( 4 ) , w e h a v e R = [ 3 2 ( τ + ρ ) ] 2 / 3 ( 2 GM w h e r e g r a t i o n c o n s t a n t h a s b e e n a b s o r b e d b y a s h i f t i n e x p r e s s i o n a s a f u n c t i o n o f τ a n d ρ c a n b e f o u n d f r o m ( 3 ) a s : = 4 GM   1 3 ( R 2 GM ) 3 / 2 + √ R 2 GM -t a n w i t h R g i v e n i n ( 5 ) . A s n t e g r a t i o n c o n s t a n t h a s b S i n c e t h e m e t r i c ( 2 ) i s s y n c t h e l i n e s o f c o n s t a fl a t m a t t e r d o m i n a t e d u n i v e r s e . T w o a l l . O n e i s t h e c o s m o l o g i c a l s c a l e t ∼ H -1 , e p u l s i v e g r a v i t a t i o n a l fi e l d o f t h e d o m a i n w a l l . I n t h i s c a s e , f o r t /lessmuch R /lessmuch t σ , t h e d o m a i n w a l l i s . E v e n t u a l l y, w h e n t h e w a l l f a l l s w i t h i n t h e h o r i z o n , i t b l a c k h o l e o f r a d i u s R S = 2 GM ∼ t 2 /t σ /lessmuch t . t h e o p p o s i t e l i m i t , w h e r e R /greatermuch t σ . I n t h i s c a s e , t h e w a l l r e p e l s s f a s t e r t h a n t h e a m b i e n t e x p a n s i o n r a t e . A s w e s h a l l s e e , t h i s . T h e d u s t w h i c h w a s o r i g i n a l l y i n t h e i n t e r i o r o f t h e w a l l g o e s i n t o e n t F R W u n i v e r s e w e a r e l e f t w i t h a b l a c k h o l e r e m n a n t c y s t e d i n a e ff e c t o f t h e d o m a i n w a l l , l e t u s fi r s t d i s c u s s t h e m a t c h i n g o f S c h w a r z s c h i l d Sc h w a r zs c hi l d t o a dus t co s m o l o g y S c h w a r z s c h i l d m e t r i c ds 2 = -( 1 -2 GM R ) dT 2 + ( 1 -2 GM R ) -1 dR 2 + R 2 d Ω 2 . ( 1 ) s , t h i s t a k e s t h e f o r m ds 2 = -d τ 2 +2 GM R d ρ 2 + R 2 d Ω 2 . ( 2 ) e l a t i o n s -√ 2 GM R ( 1 -2 GM R ) -1 dR , ( 3 ) R GM ( 1 -2 GM R ) -1 dR . ( 4 ) ) 1 / 3 , ( 5 ) t h e o r i g i n o f t h e ρ c o o r d i n a t e . T h e  -ρ , i n t h e T e s i c 1 R = 0 R = 2 GM , C o n s i d e r a s p h e r d i ff e r e n t t i m e -s c a l e s a r e a n d t h e o t h e r i s t h e a c c e l t σ ∼ ( G σ ) -1 . F o r R /lessmuch t σ , t h e r e p u l s i v e fi c o n f o r m a l l y s t r e t c h e d b y c o s m o l o g q u i c k l y s h r i n k s u n d e r i t s t e n s i o n a n d H e r e , w e w i l l b e p r i m a r i l y i n t e r e s t e d i t h e m a t t e r a r o u n d i t w h i l e i t s s i z e g r o w s f l e a d s t o t h e f o r m a t i o n o f a w o r m h o l e . T h e d u s a b a b y u n i v e r s e , a n d i n t h e a m b i e n t F R W u n i v s p h e r i c a l r e g i o n o f o f v a c u u m . B e f o r e w e c o n s i d e r t h e e ff e c t o f t h e d o m a i n a n d a d u s t c o s m o l o g y . 1 . 1 M a t c hi ng Sc h w a r zs c hi l d t o a dus t co s m o l C o n s i d e r t h e S c h w a r z s c h i l d m e t r i c = -( 1 -2 GM R ) dT 2 + ( 1 -2 GM R ) -I n L e m a i t r e c o o r d i n a t e s , t h e s t h e f o r m = -d τ 2 + 2 GM R d ρ 2 + R 2 d Ω 2 . w h e r e τ a n d ρ a r e d e fi n e d b y t h e r e l a t i d τ = ± dT -GM ( 1 -2 GM R ) -1 dR , d ρ = ∓ dT + √ R 2 GM -2 GM R ) -1 dR . S u b t r a c t i n g ( 3 ) f r o m ( 4 ) , w e h a v e R = [ 3 2 ( τ + ρ ) ] 2 / 3 ( 2 GM ) 1 / i n t e g r a t i o n c o n s t a n t h a s b e e n a b s o r b e d b y a s h i f t i n r T a s a f u n c t i o n o f τ a n d ρ c a n b e f o u n d f r o m ( 3 ) ± T = 4 GM   1 3 ( R 2 GM ) 3 / 2 + √ R 2 o n d i n t e g r a t i o n o n o u s , t h e 1 R = 0 R C o n s i d d i ff e r e n t t i m e a n d t h e o t h e r i s t σ ∼ ( G σ ) -1 . F o r R /lessmuch t σ , t h e r c o n f o r m a l l y s t r e t c h e d b y q u i c k l y s h r i n k s u n d e r i t s t e n H e r e , w e w i l l b e p r i m a r i l y i n t h e m a t t e r a r o u n d i t w h i l e i t s s i z e l e a d s t o t h e f o r m a t i o n o f a w o r m h o l e . a b a b y u n i v e r s e , a n d i n t h e a m b i e n t F R s p h e r i c a l r e g i o n o f o f v a c u u m . B e f o r e w e c o n s i d e r t h e e ff e c t o f t h e d o m a i a d u s t c o s m o l o g y . 1 . 1 a t c hi ng Sc h w a r zs c hi l d t o a dus t co s C o n s i d e r c h w a r z s c h i l d m e t r i c ds 2 = -( 1 -2 GM R I n L e m a i t r e c o o r d i n a h i s t a k e s t w h e r e τ a n d ρ a r e d e fi n e d b y S u b 1 Sphe R = 0 R = 2 GM C o n s i d e r a s p h d i ff e r e n t t i m e -s c a l e s a a n d t h e o t h e r i s t h e a c c t σ ∼ ( G σ ) -1 . F o r R /lessmuch t σ , t h e r e p u l s i v e c o n f o r m a l l y s t r e t c h e d b y c o s m o l q u i c k l y s h r i n k s u n d e r i t s t e n s i o n a n H e r e , w e w i l l b e p r i m a r i l y i n t e r e s t e h e m a t t e r a r o u n d i t w h i l e i t s s i z e g r o w s s t o t h e f o r m a t i o n o f a w o r m h o l e . T h e d a u n i v e r s e , a n d i n t h e a m b i e n t F R W u n s p h e r e g i o n o f o f v a c u u m . B e f e c o n s i d e r t h e e ff e c t o f t h e d o m a i n w a l a n d a d u s m o l o g y . 1 . 1 M a t c hi ng Sc h w a r zs c hi l d t o a C o n s i d e r t h e S c h w a r l d m e t r i c ds I n L e m a i t r e c o o r d i n a t e s , t h i s w h e r e τ a n d ρ a r e 1 Sphe r i ca l do R = 0 R = 2 GM , T →∞ R = C o n s i d e r a s p h e r i c a l d o m a i n d i ff e r e n t t i m e -s c a l e s a r e r e l e v a n t t o a n d t h e o t h e r i s t h e a c c e l e r a t i o n t i m e t σ ∼ ( G σ ) -1 . F o r R /lessmuch t σ , t h e r e p u l s i v e fi e l d c a n b e r m a l l y s t r e t c h e d b y c o s m o l o g i c a l e x p a n s i q u s h r i n k s u n d e r i t s t e n s i o n a n d f o r m s a b l a H e w i l l b e p r i m a r i l y i n t e r e s t e d i n t h e o p p o t h e m a t o u n d i t w h i l e i t s s i z e g r o w s f a s t e r t h a n l e a d s t o t h m a t i o n o f a w o r m h o l e . T h e d u s t w h i c h w a a b a b y u n i v e n d i n t h e a m b i e n t F R W u n i v e r s e w e a r e s p h e r i c a l r e g i o n v a c u u m . B e f o r e w e c o n s h e e ff e c t o f t h e d o m a i n w a l l , l e t u s fi r s a n d a d u s t c o s m o l o g y 1 . 1 M a t c hi ng Sc h w a hi l d t o a dus t co s C o n s i d e r t h e S c h w a r z s c h i l d m e ds 2 = -( 1 -GM I n L e m a i t r e c o o r d i n a t e s , t h i s t a k e s t h e w h e r e τ a n ρ a r e d e fi n e d FIG. 13: Lemaitre coordinates in the extended Schwarzschild diagram. The left panel shows the curves ρ = 0 for the two possible choices of the double sign in Eq.(B6). In the left panel, we adopt the lower sign, indicated by -, which corresponds to radial geodesics traveling towards region II from the white hole singularty at R = 0, for different values of ρ . Note that such geodesic coordinates only cover half or the extended Schwarzschild solution. Therefore, the mass parameter M in the Schwarzschild solution corresponds to the total mass of dust that would be contained within the sphere of radius R , where the two solutions are matched. The relation (B11) can also be derived by imposing that the dust particle at the edge of the FRW metric should be a geodesic of both FRW and Schwarzschild. This is explained in Subsection III A 2. Since the normal n µ to the the matching hypersurface ρ = const. or r = const. has only radial component, and the metric is diagonal, the extrinsic curvature K µν = n µ ; ν takes a simple form given by  \na  \n,  \nn  \no  \nd  \no  \nr  \nτ  \nn  \ni  \nd  \ni  \na  \nt  \nt  \nexpression for T as a function of τ and ± T = 4 GM   1 3 ( R 2 GM ) 3 / 2 + √ R 2 GM -tanh -1 √ R 2 GM   -ρ, (6) with R given in ( ?? ). A second integration constant has been absorbed by a shift in the T coordinate. Since the metric ( ?? ) is synchronous, the lines of constant spatial coordinate are geodesics, and τ is the h e e w i t h R g i v e n i n S i n c e t h e m e t r i c ( where ∂ n = a -1 ∂ r = ( R/ 2 GM ) 1 / 2 ∂ ρ is the normal derivative, and α and β run over the temporal and angular coordinates. It is straightforward to check that K αβ has only angular components, and that these are the same on both sides of the matching hypersurface:  \nw  \nds  \nρ  \ncan be found from (3) as:  \ne  \n.  \ns  \nK ΩΩ \' = g ΩΩ \' /R. (B13)  \nHere Ω and Ω \' run over the angular coordinates. The continuity of the metric and of the extrinsic curvature means that there is no distributional source at the junction.  \nThe matching of Schwarzshild and FRW metrics can be done in two different ways, which are represented in Fig. 14. An expanding ball of matter may be embedded in a Schwarzschild  \n6  \n)  \nf  \ns  \nn  \ni  \n2  \nt  \nh  \ne  \nr  \nr  \nk  \ne  \ns  \ne  \no  \na  \no  \nr  \nm  \nf  \nt  \nh  \nt  \nt  \nh  \ne  \ni  \ng  \nt  \n2  \ne  \nGM  \n=  \nR  \n±  \ndT  \nt  \nds  \nR  \n2  \n+  \nτ  \nd  \nτ  \n-  \n=  \ns  \n1  \nt  \no  \nn  \nd  \ni  \nl  \n√  \nt  \ne  \nb  \n=  \na  \ne  \nl  \na  \n=  \nr  \ne  \n2  \nd  \na  \ne  \ns  \nh  \nd  \nρ  \n±  \nr  \nw  \nb  \n=  \nτ  \n)  \ne  \n-  \ny  \nd  \nds  \nr  \nh  \nn  \na  \ns  \n:  \n2  \ne  \ne  \nt  \nfi  \ne  \nt  \nb  \nd  \nGM  \nm  \nf  \na  \nd  \nn  \n4  \nb  \no  \na  \n,  \n(  \ni  \np  \nt  \nh  \nl  \na  \nρ  \na  \nd  \ni  \nr  \nt  \ne  \ns  \nn  \na  \ne  \ni  \nr  \no  \nh  \na  \n)  \nt  \nt  \nt  \nr  \ni  \no  \nn  \nn  \nn  \ny  \n3  \n(  \na  \ng  \ns  \no  \nc  \ns  \nr  \nn  \ns  \nn  \na  \nn  \ne  \na  \nr  \ng  \nw  \nh  \nt  \nt  \no  \ne  \n(  \nn  \no  \n)  \n,  \n(  \nt  \n)  \ng  \nw  \n4  \nf  \nr  \no  \nT  \nt  \nn  \ne  \ni  \na  \nr  \nw  \nw  \nh  \n)  \nf  \n(  \nm  \nf  \nc  \n3  \np  \ng  \ni  \nt  \nr  \nn  \ns  \ni  \ne  \nx  \na  \ne  \no  \n3  \nn  \nn  \ne  \nh  \ni  \ni  \no  \ne  \nx  \nw  \ns  \np  \nn  \ns  \nr  \ne  \ne  \nr  \ne  \nu  \ns  \nt  \nb  \nS  \ni  \na  \nb  \nt  \nr  \nu  \nc  \nt  \nf  \na  \ns  \nc  \nt  \nc  \nn  \ni  \ne  \nl  \nS  \nn  \nh  \n-  \nc  \no  \nw  \nf  \nt  \ne  \nh  \ne  \nr  \nw  \nr  \nh  \nt  \nt  \nz  \nc  \n±  \nh  \nR  \no  \nM  \n√  \nd  \ne  \ns  \na  \na  \no  \na  \nds  \nt  \ns  \nn  \nl  \na  \ni  \no  \nn  \nd  \ns  \ni  \na  \nh  \nh  \nr  \nl  \ni  \np  \nh  \ne  \nt  \ni  \nm  \nw  \nf  \nu  \na  \ny  \nh  \na  \nl  \nc  \nn  \nt  \nd  \nh  \nt  \nu  \ne  \nr  \ni  \n(  \nm  \nm  \nr  \nl  \ne  \ne  \na  \na  \nr  \no  \nl  \nt  \nhi  \nB  \ne  \nd  \nhi  \nτ  \ns  \ng  \ny.  \nn  \ni  \nd  \nl  \ni  \nh  \nc  \n=  \nC  \ne  \nn  \n1  \ns  \nc  \nn  \nd  \ne  \nb  \nτ  \no  \ni  \no  \nd  \nd  \no  \ni  \ne  \nr  \nd  \ny  \nM  \nr  \nt  \na  \nh  \nn  \ne  \na  \nt  \nt  \ns  \nh  \n=  \n±  \ndT  \n+  \nw  \ne  \no  \ns  \nc  \nc  \nm  \nn  \no  \no  \nn  \nl  \no  \nf  \ns  \no  \nSc  \ns  \nc  \nt  \ne  \nhi  \nS  \nng  \nw  \nz  \na  \nr  \n1  \nc  \n2  \n-  \ni  \ng  \nh  \na  \nn  \no  \nd  \ny.  \nf  \ne  \nSphe  \nh  \nr  \nw  \ns  \ni  \n=  \nt  \na  \nzs  \nl  \nm  \no  \nd  \nv  \nn  \ni  \na  \nl  \nt  \nn  \nc  \nh  \no  \nh  \nf  \nt  \nu  \nca  \ne  \nhi  \ni  \nR  \nc  \n=  \ne  \nt  \ni  \na  \nd  \nr  \nc  \n0  \n2  \n(  \nk  \ne  \ns  \n2  \n-  \n1  \nt  \nr  \nR  \n-  \nh  \n=  \ns  \nds  \nt  \n,  \ne  \nt  \nr  \ne  \na  \n1  \nSpherical domain wall in dust cosmology  \ndT  \no  \ne  \nR  \no  \nd  \nn  \ny  \no  \n-  \nR  \n= 0  \nR  \n= 2  \n√  \ng  \n2  \nl  \ny  \nl  \n2  \no  \nl  \no  \ns  \ni  \ni  \ny  \nr  \ne  \ni  \nl  \na  \nh  \nu  \ne  \ne  \nm  \nff  \nw  \ni  \nl  \ni  \nt  \no  \na  \n.  \ne  \nc  \nt  \n2  \ns  \nr  \nm  \nm  \nb  \nGM  \n,  \nT  \no  \nf  \nn  \ns  \nn  \nt  \ni  \ne  \nz  \nh  \ni  \na  \n=  \nr  \nR  \ne  \no  \nl  \nn  \ne  \n→∞  \ne  \ng  \ni  \nn  \ns  \ng  \ne  \nc  \nd  \nt  \nR  \n= 0  \nT  \nc  \na  \na  \nn  \nl  \na  \nf  \no  \nm  \ne  \nd  \nr  \na  \n.  \nr  \nl  \ni  \nt  \ns  \np  \nt  \nr  \n=  \nl  \na  \nl  \ni  \n= t  \nm  \na  \ne  \ne  \na"  \nT  \n=  \nt  \ne  \nc  \nh  \ni  \ne  \n∞ d  \ne  \nd  \ns  \nd  \nl  \ni  \no  \no  \na  \nt  \ne  \no  \nn  \ns  \nf  \ne  \nu  \nu  \ni  \nb  \ns  \na  \nn  \nt  \nd  \nb  \na  \nu  \n-∞  \na  \ny  \nx  \nb  \nl  \na  \nn  \ni  \nn  \ni  \no  \ni  \nr  \nl  \ns  \na  \nc  \na  \nr  \ng  \no  \nb  \ne  \nl  \ne  \nf  \nl  \nc  \ne  \nc  \no  \nn  \ng  \ni  \np  \nm  \no  \na  \no  \ns  \nw  \nh  \nw  \nt  \no  \ns  \ny  \na  \na  \nh  \nw  \nc  \nSc  \no  \nl  \np  \ng  \nc  \no  \nh  \nr  \ns  \nd  \ni  \nm  \nl  \n1  \n-  \ne  \nm  \nf  \nR  \nn  \nd  \nl  \nr  \nm  \nl  \n(  \nt  \no  \n2  \nd  \na  \nτ  \n+  \nl  \nt  \nt  \na  \ni  \nr  \no  \nd  \no  \nn  \ni  \ne  \nr  \no  \ni  \nd  \nw  \no  \nr  \nc  \nT  \ne  \ns  \nh  \na  \ne  \ne  \ns  \nr  \ne  \no  \nb  \nt  \nt  \np  \nm  \nn  \ni  \na  \nl  \nr  \ne  \ns  \nt  \nm  \ni  \na  \nn  \n2  \ne  \nv  \n=  \nr  \na  \nl  \ns  \ne  \ni  \nm  \nρ  \nn  \nn  \na  \nl  \nu  \nn  \n-  \nb  \nc  \ns  \na  \nn  \ne  \nr  \nh  \nf  \ni  \ne  \nt  \np  \ni  \na  \n2  \ns  \ni  \nb  \nR  \nl  \no  \na  \nn  \ne  \n2  \n1  \n-  \n.  \nx  \ne  \nt  \ne  \nh  \nw  \ni  \ni  \nr  \n=  \nc  \na  \ns  \na  \nc  \ni  \ng  \nn  \ns  \ne  \no  \n(  \np  \no  \nn  \ne  \ns  \n1  \ni  \np  \n-  \no  \ny  \n.  \no  \ni  \nh  \nn  \nd  \nh  \n,  \no  \nk  \nw  \nl  \ne  \nConsider a spherical domain wall embedded in a spatially flat matter dominated universe.  \nTwo  \nx  \nn  \ne  \nv  \ne  \ni  \ns  \np  \nh  \np  \nc  \ns  \nh  \ne  \no  \ns  \na  \ne  \na  \ns  \nt  \nd  \nz  \nw  \ns  \ns  \ne  \nt  \ndo  \nfi  \nT  \n-  \ns  \n,  \nca  \nl  \n=  \n2  \ne  \ns  \na  \nk  \n→∞  \ns  \n-  \ne  \na  \nt  \ne  \na  \no  \nn  \nt  \n∞  \nt  \no  \ne  \nd  \n2  \nGM  \n)  \ns  \nt  \nI  \nΩ  \np  \na  \na  \np  \nr  \ng  \ne  \na  \ni  \nn  \na  \nl  \nn  \ni  \na  \no  \nn  \ns  \ni  \nt  \nm  \ne  \ns  \nt  \nc  \ne  \no  \ne  \nng  \nt  \nr  \na  \nl  \nt  \nh  \nc  \ni  \nR  \nl  \nc  \ns  \no  \ni  \nm  \nb  \nf  \nm  \nd  \nb  \no  \nh  \nn  \na  \nl  \nl  \nt  \nt  \na  \nu  \ns  \ni  \nc  \nt  \no  \ne  \n2  \ne  \nt  \nb"  \ndifferent time-scales are relevant to the dynamics of such a wall. One is the cosmological scale  \nt  \nn  \nf  \ne  \nh  \nh  \nc  \nr  \nt  \na  \nf  \nR  \na  \nb  \nt  \na  \nn  \nh  \ns  \nR  \nh  \nw  \na  \ns  \nl  \ni  \nS  \nt  \ne  \ne  \nc  \nu  \nh  \ni  \ni  \np  \no  \ns  \ni  \nr  \na  \nd  \ne  \nh  \nr  \nw  \nw  \nh  \n∼  \ns  \nn  \nn  \np  \nt  \nv  \n1  \nH  \ns  \nt  \nd  \n-  \n-  \n1  \ns  \ne  \np  \ni  \no  \nr  \n,  \ni  \nu  \nf  \nR  \ne  \nl  \ns  \np  \na  \nn  \nx  \na  \nn  \nw  \n,  \nt  \ne  \nB  \nt  \ni  \nm  \nw  \nf  \no  \nr  \nu  \nu  \na  \na  \nl  \na  \ns  \nh  \na  \nh  \nh  \nc  \nhi  \nr  \nt  \nw  \nw  \nS  \ns  \nh  \nds  \ne  \ni  \nh  \nca  \nl  \ne  \nl  \nt  \nw  \nu  \ni  \nds  \ns  \nn  \nm  \ne  \nfi  \nl  \nd  \nd  \nGM  \no  \ng  \nR  \nd  \n2  \na  \nc  \nR  \ni  \nc  \nn  \nn  \nd  \na  \nd  \nl  \nl  \ni  \ne  \nt  \nc  \nr  \no  \nn  \no  \na  \nr  \ne  \ne  \n,  \n.  \no  \nb  \nl  \nh  \ns  \nd  \nE  \nl  \nI  \n.  \ne  \nl  \nt  \nt  \nf  \ns  \nr  \nr  \nc  \ne  \na  \ne  \nv  \na  \nt  \ne  \ni  \nh  \nn  \nl  \ni  \nl  \nr  \ns  \na  \nr  \ne  \n2  \ne  \nl  \no  \nm  \nc  \nd  \ni  \na  \na  \nl  \ne  \nc  \nl  \n,  \ne  \nx  \nr  \ne  \n.  \nt  \ne  \nr  \ne  \nh  \nt  \na  \nb  \nt  \nn  \n→-  \ne  \na  \nm  \nv  \nT  \nl  \nu  \no  \nv  \ni  \ne  \nl  \ne  \no  \nt  \nh  \ne  \na  \ne  \na  \ns  \no  \nr  \ni  \nt  \ni  \nl  \nw  \nc  \nh  \ne  \nm  \nv  \ns  \nt  \ni  \nn  \nh  \ni  \n,  \nand the other is the acceleration time-scale due to the repulsive gravitational field of the domain wall  \na  \nd  \na  \nu  \nt  \nc  \ne  \nc  \nt  \nu  \nW  \nR  \nF  \nh  \nh  \na  \ne  \nt  \ni  \nd  \nt  \nfi  \nl  \nv  \nl  \nw  \np  \no  \nT  \nt  \nσ  \ne  \nr  \nt  \ns  \ni  \nw  \ns  \ne  \n.  \n∼  \na  \nc  \n1  \np  \nh  \nh  \ns  \n(  \nGσ  \n)  \nk  \no  \n-  \nv  \nh  \ni  \ni  \ns  \nl  \nH  \nFor  \nR  \na  \ni  \nr  \nv  \nu  \nr  \nn  \ne  \np  \nc  \nn  \n1  \ne  \ne  \no  \n-  \nl  \ns  \nt  \nt  \nw  \ne  \ns  \na  \no  \ni  \ny  \n∼  \ne  \nw  \nd  \nn  \nh  \ni  \nl  \no  \nt  \n.  \nm  \no  \no  \nm  \n/lessmuch  \ne  \nc  \nh  \ng  \nt  \nu  \nt  \nf  \nn  \nl  \ni  \ne  \nh  \na  \nσ  \nn  \ni  \nl  \ne  \nn  \ns  \nl  \na  \nd  \ns  \nc  \no  \nr  \nn  \nt  \n,  \nthe repulsive field can be ignored.  \nIn this case, for  \nt  \ny,  \ns  \nt  \nl  \nfi  \n-  \n√  \n2  \n+  \nl  \ne  \ni  \nt  \nn  \no  \nn  \ns  \na  \nt  \nl  \nr  \ne  \na  \nl  \na  \nw  \ng  \ns  \nt  \n,  \n2  \n,  \nSphe  \nu  \na  \nr  \nw  \ni  \ndT  \ni  \nr  \nt  \nr  \ne  \ne  \ny  \np  \nGM  \nt  \na  \ne  \na  \ndo  \ny  \ne  \ne  \na  \n→∞  \na  \n-  \nh  \nT  \ne  \na  \nl  \ns  \nc  \ne  \ne  \na  \nt  \nh  \nr  \ni  \nc  \nl  \no  \nr  \no  \nd  \np  \n1  \nu  \nl  \nm  \nc  \nh  \ne  \ne  \nr  \ne  \nt  \ns  \na  \no  \nc  \ny  \na  \n-  \nl  \ns  \nb  \nc  \nGM  \nr  \n2  \ni  \ni  \ns  \n,  \ns  \ne  \nr  \ni  \nd  \n=  \ni  \nR  \nn  \nl  \n=  \n0  \nc  \nd  \nm  \nr  \nSphe  \nca  \ne  \ns  \nR  \n±  \n-  \nc  \nh  \nh  \na  \nb  \no  \np  \nh  \ne  \nt  \nm  \nn  \nn  \ns  \ni  \no  \no  \nn  \nl  \ne  \ni  \nn  \nd  \no  \nt  \na  \nr  \ne  \ni  \ndR  \ns  \no  \ng  \nr  \nn  \nl  \ni  \ni  \ne  \nv  \ne  \ni  \no  \nt  \nn  \nt  \ny  \ns  \n=  \na  \nm  \nh  \ns  \nGM  \n2  \ne  \nt  \nr  \ng  \nm  \nt  \nn  \ne  \nl  \nt  \nσ  \ne  \nR  \ne  \nd  \no  \n(  \ns  \nw  \nE  \nn  \ng  \nc  \n.  \n2  \nGM  \nr  \ns  \nm  \nf  \nn  \ne  \ni  \na  \n+  \n√  \ni  \no  \ne  \ni  \ns  \n=  \n,  \nthe domain wall is  \ny  \ni  \nt  \n/lessmuch  \nd  \nt  \ne  \na  \ne  \ne  \nd  \nl  \n1  \nfl  \nd  \nt  \na  \nn  \ny  \nn  \nl  \nv  \n.  \nd  \ne  \nl  \ne  \ns  \na  \ne  \ni  \ni  \nn  \ni  \nl  \nl  \ni  \na  \ni  \nO  \ne  \nr  \ng  \ns  \nt  \nr  \ne  \nt  \n,  \n,  \ni  \ne  \na  \nt  \ns  \ne  \n1  \nn  \nw  \nr  \n.  \no  \na  \n1  \ni  \ni  \nM  \na  \nt  \nu  \ne  \nm  \ns  \nT  \nt  \nv  \ng  \nl  \nh  \nd  \na  \nt  \nt  \n.  \nt  \na  \ne  \ns  \ni  \nt  \nd  \nh  \no  \nn  \ni  \ne  \nr  \no  \nk  \no  \nh  \nc  \nd  \no  \ns  \nb  \na  \n/lessmuch  \ne  \ndT  \nr  \ns  \nn  \nm  \nτ  \nn  \nl  \ne  \nh  \ne  \n.  \ne  \no  \nz  \na  \ns  \ns  \ne  \nn  \nn  \nn  \ni  \nR  \nT  \ns  \nl  \nl  \na  \nm  \no  \ne  \na  \nr  \nh  \nb  \ni  \nt  \ni  \nm  \ne  \nu  \nl  \nt  \nf  \nW  \nd  \ns  \ni  \np  \nt  \no  \nr  \nt  \ns  \na  \n3  \n/  \nF  \nw  \ng  \ne  \nn  \nh  \no  \np  \ni  \n2  \ne  \nr  \ne  \n/  \n3  \nw  \n)  \n1  \na  \nh  \nu  \ne  \n2  \nGM  \nt  \na  \nd  \nt  \na  \n(  \ne  \nr  \n1  \ne  \ns  \nc  \ne  \nn  \nW  \ne  \no  \nt  \nr  \nR  \nc  \nd  \ns  \ne  \no  \ns  \no  \nl  \ni  \ny  \nl  \na  \ne  \nl  \n,  \ne  \ni  \ni  \nt  \nd  \nl  \ns  \nh  \ne  \nfi  \nr  \ns  \nn  \nt  \ne  \nw  \ns  \nl  \ni  \ng  \nn  \na  \nr  \nr  \ns  \nm  \nn  \nt  \na  \n,  \no  \nw  \ni  \na  \nn  \nn  \ni  \nv  \ne  \nl  \nr  \ne  \nt  \nu  \nf  \ne  \nu  \na  \nr  \ne  \nt  \ni  \ne  \nt  \nk  \nl  \na  \nt  \nh  \ni  \nn  \nf  \nb  \na  \nc  \ne  \nρ  \nh  \nt  \ns  \no  \nf  \nu  \nfi  \nt  \ns  \ns  \nr  \ni  \ng  \n,  \nr  \nw  \nf  \nt  \no  \nh  \ne  \n:  \ns  \na  \ns  \ni  \nn  \n)  \n3  \na  \nw  \ni  \na  \ni  \ne  \no  \nh  \nh  \nh  \ne  \nw  \nc  \nh  \n)  \nh  \ne  \nT  \nρ  \nconformally stretched by cosmological expansion. Eventually, when the wall falls within the horizon, it  \na  \n,  \nd  \ni  \nt  \nr  \ne  \nn  \nm  \ni  \no  \nz  \nr  \nf  \ni  \nn  \nh  \nw  \ni  \nu  \ns  \nl  \ne  \nr  \nh  \nr  \ne  \nn  \ne  \ni  \nt  \ne  \ne  \nb  \ny  \na  \n,  \ns  \nt  \nh  \no  \nt  \no  \ni  \nm  \ns  \nt  \na  \ne  \nR  \nC  \nn  \no  \n.  \nσ  \n=  \nr  \ns  \nc  \nt  \nh  \no  \nm  \nR  \nd  \nfi  \n-  \n-  \na  \nn  \nh  \na  \ne  \ne  \nl  \ns  \nd  \n-  \nc  \ne  \nt  \na  \ns  \ni  \ni  \nn  \no  \ne  \nm  \nd  \nn  \nd  \nGM  \ni  \nt  \nv  \na  \no  \n1  \nquickly shrinks under its tension and forms a black hole of radius  \nw  \nh  \ns  \na  \nw  \nw  \ni  \ne  \nt  \no  \nm  \nh  \ni  \no  \ne  \nr  \nl  \ns  \nt  \nn  \nz  \nR  \ni  \ne  \nn  \ne  \nm  \nw  \nC  \ne  \nco  \ne  \nt  \n2  \ne  \no  \nR  \ns  \ni  \nm  \nHere, we will be primarily interested in the opposite limit, where  \nn  \ni  \nr  \n-  \nl  \nx  \nh  \nW  \na  \ns  \nh  \nw  \nF  \nn  \nh  \ni  \nu  \nt  \no  \nr  \ne  \ne  \nr  \nl  \nc  \nc  \nm  \ni  \ny  \ni  \nt  \np  \ne  \ni  \nr  \nd  \ne  \nR  \nd  \nm  \nf  \nn  \n=  \nt  \no  \nS  \nR  \nl  \np  \nw  \nh  \no  \n2  \nr  \no  \nm  \n(  \nn  \nt  \nt  \ns  \nc  \nu  \ni  \nl  \nl  \ne  \nl  \no  \nd  \ns  \nm  \ns  \n1  \nt  \nfi  \nu  \nt  \nt  \ndus  \nd  \no  \na  \nl  \nl  \ne  \nl  \n,  \nw  \no  \nt  \nh  \n,  \ne  \na  \nthe matter around it while its size grows faster than the ambient expansion rate. As we shall see, this  \ny  \n/greatermuch  \na  \na  \nn  \ni  \ne  \nw  \nn  \nn  \na  \nt  \nf  \no  \nt  \nσ  \nh  \no  \na  \nt  \na  \nleads to the formation of a wormhole. The dust which was originally in the interior of the wall goes into  \no  \nn  \ns  \nh  \ni  \nt  \nt  \nl  \ni  \ne  \na  \nd  \ni  \na  \nn  \n2  \n.  \nn  \ni  \nh  \n.  \nn  \ni  \nc  \nw  \nw  \na baby universe, and in the ambient FRW universe we are left with a black hole remnant cysted in a  \nb  \nl  \nt  \na  \nd  \nt  \nn  \n,  \nl  \nn  \ni  \ne  \nr  \nt  \nv  \nh  \na  \nf  \na  \nr  \ne  \nt  \nh  \nspherical region of of vacuum.  \ne  \ne  \n)  \ni  \no  \nk  \ns  \nf  \no  \n,  \nn  \nσ  \no  \n/lessmuch  \ne  \n+  \n√  \ni  \nr  \nc  \n1  \ne  \nt  \n-  \ndR  \ne  \ne  \ni  \ni  \nz  \no  \no  \nh  \nBefore we consider the effect of the domain wall, let us first discuss the matching of Schwarzschild  \nl  \ne  \nn  \nn  \nm  \n/t  \na  \nl  \ns  \na  \nc  \na  \nl  \ns  \ne  \ni  \nd  \nr  \ns  \ni  \nt  \ni  \no  \nr  \nw  \no  \nm  \nw  \nand a dust cosmology.  \ni  \ns  \nc  \nc  \no  \nt  \ni  \nh  \nA  \ne  \na  \no  \ni  \nn  \no  \ne  \nh  \nh  \nn  \nSpherical domain wall in dust cosmology  \ni  \ng  \nm  \nt  \ne  \nw  \nt  \no  \no  \nl  \no  \n,  \n.  \nz  \na  \no  \ns  \nr  \no  \nf  \ne  \ne  \n.  \nt  \nMatching Schwarzschild to a dust cosmology  \nh  \n<!-- image -->  \nr  \nf  \nm  \n2  \n+2  \nR  \n=  \nGM  \np  \no  \nt  \nn  \na  \nx  \ns  \ne  \ni  \n<!-- image -->  \nn  \nr  \ns  \nc  \nn  \nt  \nt  \nr  \nσ  \nd  \nu  \no  \ns  \nΩ  \n2  \n2  \ns  \ni  \nl  \nt  \nl  \no  \ng  \ny  \nR  \n2  \nGM  \ns  \nb  \nco  \n√  \nco  \nh  \nl  \no  \n-  \na  \nn  \n-  \no  \nr  \ny  \ns  \n+  \ng  \nt  \na  \ne  \nb  \n2  \nn  \ndus  \no  \ne  \nb  \ndT  \ni  \np  \na  \nt  \nm  \na  \nGM  \nR  \n2  \n.  \ns  \nhi  \nn  \n2  \nGM  \nm  \ns  \n)  \ns  \nh  \na  \no  \nt  \nco  \nn  \nt  \na  \nd  \nn  \nt  \nt  \nl  \no  \no  \nc  \nc  \nl  \no  \ni  \n=  \nr  \nzs  \np  \ns  \na  \nt  \n)  \n2  \ns  \n1  \no  \nf  \nc  \nt  \nR  \nn  \no  \nn  \ns  \nc  \n-  \ndus  \ns  \n(  \na  \nt  \n-  \nn  \ne  \ns  \ni  \ne  \na  \ne  \nt  \ns  \nr  \nf  \nSpherical domain wall in dust cosmology  \nh  \nl  \nr  \nt  \nds  \nd  \ns  \n,  \ne  \nr  \ns  \ni  \nk  \ne  \na  \ne  \nt  \nt  \na  \ni  \nt  \nm  \nt  \nh  \n,  \ni  \nl  \nd  \ne  \nn  \ne  \ni  \nc  \nh  \ns  \no  \ns  \ns  \ne  \nr  \nc  \nt  \ns  \nh  \no  \n2  \nGM  \n2  \n+  \n2  \nds  \nt  \n1  \n-  \ndT  \no  \ns  \nn  \nt  \ng  \ne  \nt  \nt  \ni  \ns  \no  \nl  \na  \n-  \n(  \nGM  \nh  \ni  \ni  \na  \n,  \n2  \nGM ρ  \n=  \nconst.  \ni  \n2  \nd  \ns  \nh  \ne  \ns  \ns  \ni  \no  \nh  \nv  \nr  \ne  \nt  \na  \na  \nu  \nt  \ne  \n2  \nt  \nl  \ni  \na  \nn  \nc  \nh  \ns  \na  \ns  \ny  \nw  \ns  \nm  \ne  \n→∞  \n→-∞  \n-  \ne  \na  \nl  \nt  \nh  \nH  \nst  \nr  \nt  \nq  \nr  \nh  \na  \nt  \nt  \nn  \nt  \nfi  \no  \nt  \ncon  \ne  \na  \nt  \ng  \nz  \nΩ  \nR  \nρ  \na  \nl  \ni  \n2  \nt  \ncon  \n.  \na  \nGM  \no  \ne  \nr  \ni  \nh  \nl  \ne  \nt  \nt  \n/lessmuch  \nn  \ni  \nd  \ne  \no  \nh  \n,  \nr  \na  \nn  \ns  \ne  \nl  \ni  \nConsider the Schwarzschild metric  \no  \nt  \ns  \np  \nh  \ne  \nρ  \n= 2  \nGM,T  \n.  \na  \ns  \n.  \nt  \nh  \n.  \na  \n.  \nl  \nl  \nr  \na  \nc  \nc  \ne  \ni  \nl  \nt  \nc  \ns  \nt  \na  \ns  \nd  \na  \ne  \ns  \n2  \nl  \nt  \nσ  \ne  \nConsider a spherical domain wall embedded in a spatially flat matter dominated universe.  \nTwo  \nv  \ne  \ns  \nr  \nt  \ni  \ne  \n→∞  \n→-∞  \ns  \n1  \nw  \nr  \ne  \nS  \na  \ne  \ne  \ni  \n-  \nk  \ne  \nc  \nw  \nw  \n/t  \nn  \nt  \n,  \n-  \n1  \n.  \n1  \ne  \no  \nn  \nt  \nh  \nn  \n)  \n3  \n/  \n=  \n1  \n,  \n=  \nR  \nt  \ni  \nl  \ns  \ne  \no  \nr  \nr  \ne  \nk  \ne  \nt  \na  \no  \nr  \nR  \ne  \nm  \nd  \nf  \ni  \ne  \nl  \nc  \nh  \nz  \nh  \n-  \nr  \nA  \nR  \ng  \ni  \nu  \nn  \ne  \no  \ns  \nd  \no  \ni  \ni  \np  \ns  \np  \nt  \nn  \ni  \no  \nConsider a spherical domain wall embedded in a spatially flat ma  \nh  \ns  \n)  \n.  \nσ  \n/lessmuch  \nhi  \ne  \n5  \nGM  \nt  \nc  \nh  \n,  \n(  \nd  \nw  \nn  \na  \nn  \nc  \nB  \ni  \ne  \ndT  \nr  \ne  \n(  \no  \ng  \no  \ns  \nConsider a spherical domain wall embedded in a spatially flat matter dominated universe.  \nTwo  \no  \n√  \ns  \nR  \nm  \nl  \nh  \ns  \nm  \nm  \nc  \n2  \nGM  \ne  \ni  \ns  \n)  \nh  \ne  \nr  \nd  \no  \no  \nc  \nw  \ne  \n,  \ns  \nh  \ny  \nh  \ne  \nl  \n)  \n(  \n2  \na  \n3  \ns  \no  \nt  \n∼  \n√  \nm  \no  \nl  \nl  \ns  \nh  \ng  \ni  \nf  \na  \n=  \n/lessmuch  \nt  \nl  \ni  \na  \ns  \n)  \na  \n2  \ne  \nl  \no  \nt  \nt  \nR  \na  \nl  \n→∞  \nf  \nn  \nn  \nt  \nu  \nfi  \n,  \nh  \ns  \na  \nn  \n/greatermuch  \nh  \no  \na  \ne  \nt  \nc  \nt  \na  \nM  \n→-∞  \n-  \ng  \no  \nl  \nn  \nt  \na  \n2  \nGM  \no  \nw  \nd  \nl  \ny  \nh  \ng  \ni  \nr  \nh  \nc  \no  \ni  \ni  \ng  \nst  \ne  \nu  \n2  \n1  \nc  \na  \nl  \n,  \nl  \nh  \ni  \ni  \nh  \nt  \nh  \nt  \nb  \ni  \ns  \nd  \nds  \ns  \ne  \nd  \nt  \ni  \ny  \nn  \ni  \ne  \ni  \nl  \nd  \n1  \nh  \na  \no  \nl  \nS  \nw  \ng  \nc  \nSc  \n.  \n2  \n/t  \nt  \ni  \nt  \n/lessmuch  \ne  \n,  \n1  \nH  \n-  \nn  \n+  \nσ  \ns  \nH  \n.  \nl  \nt  \na  \nh  \n,  \ndT  \n2  \n+  \n1  \n-  \ne  \nt  \ndR  \n2  \ne  \nd  \nh  \na  \ne  \nh  \nt  \nl  \nl  \nl  \nσ  \n,  \nn  \na  \nt  \nl  \nc  \nn  \nc  \nw  \nh  \n+  \nR  \n2  \n/lessmuch  \ni  \nt  \ne  \n±  \na  \nd  \nΩ  \nC  \nd  \nw  \nl  \nx  \ns  \na  \nt  \n2  \n.  \n(1)  \nn  \n,  \nd  \ndifferent time-scales are relevant to the dynamics of such a wall. One is  \n=  \n(  \nt  \nh  \n3  \nGM  \n(  \n(  \n)  \ne  \nc  \nh  \nh  \nc  \n1  \ni  \nt  \ne  \nM  \ne  \nr  \nng  \n∼  \n-  \ni  \ne  \n,  \ns  \nr  \ne  \n-  \nn  \n:  \nn  \nR  \nw  \ndifferent time-scales are relevant to the dynamics of such a wall. One is the cosmological scale  \nt  \ns  \n-  \ni  \nm  \nh  \nt  \nGM  \ny  \n1  \nn  \na  \na  \ne  \ns  \nl  \nt  \nh  \nz  \ne  \ni  \nc  \nh  \nA  \ne  \ns  \nb  \ny  \nt  \na  \na  \ne  \ne  \no  \ns  \ne  \na  \nh  \nz  \nn  \ns  \nt  \nk  \nh  \ni  \nf  \nd  \ne  \nw  \nn  \na  \nt  \ns  \ns  \na  \nc  \nh  \no  \n2  \nr  \nc  \n.  \na  \nn  \nt  \na  \nw  \ns  \na  \ne  \nr  \n.  \nt  \nn  \nt  \na  \nl  \na  \nl  \nd  \ndifferent time-scales are relevant to the dynamics of such a wall. One is the cosmological scale  \nt  \nco  \nc  \nGM  \n-  \n2  \n-  \nt  \nt  \ne  \ni  \ni  \nn  \na  \nh  \nm  \nt  \nco  \n)  \ns  \np  \ns  \nR  \nR  \np  \nb  \nfl  \n2  \n=  \ni  \nf  \n2  \nu  \ns  \nt  \n/lessmuch  \n=  \ne  \nb  \n)  \nv  \nn  \nt  \ne  \nt  \n/  \nn  \nt  \n2  \ns  \n2  \nst  \nρ  \n(  \nm  \n1  \nρ  \ny  \nO  \na  \nl  \n∼  \nf  \nn  \no  \ncon  \nm  \nt  \no  \nw  \n/lessmuch  \ne  \ns  \nt  \nc  \no  \nh  \ni  \nl  \nl  \ny  \nt  \ncon  \na  \n.  \nR  \no  \ns  \nf  \ne  \no  \nS  \nh  \nz  \nt  \nd  \np  \nl  \nu  \n2  \nr  \nd  \nc  \n-  \ns  \no  \n1  \n)  \nΩ  \n2  \nd  \nd  \nd  \ni  \n2  \nc  \n.  \no  \n2  \n2  \n+  \no  \nr  \nm  \ny  \nΩ  \nf  \nb  \ni  \nτ  \nand  \nρ  \ncan be found from (3) as:  \nd  \ne  \nn  \ns  \nt  \na  \no  \nn  \nn  \nd  \ns  \ne  \ni  \nt  \ns  \n(  \nt  \nr  \nn  \nr  \nd  \n+  \n2  \nGM  \n2  \ni  \ne  \ne  \ng  \nv  \nt  \nd  \ne  \ni  \nn  \nΩ  \nds  \no  \nn  \ne  \n-  \n.  \ns  \ni  \nand the other is the acceleration time-scale due to the repulsive gravitational field of the domain wall  \nt  \nn  \nb  \nO  \na  \nΩ  \na  \nand the other is the acceleration time-scale due to the repulsive gravit  \nr  \nt  \ne  \nf  \na baby universe, and in the ambient FRW universe we are left with a  \ni  \no  \nl  \ne  \nd  \ne  \n1  \ni  \nz  \ns  \n1  \ni  \nh  \na  \nR  \ny  \nl  \ng  \nr  \na  \nt  \ni  \nc  \nk  \nm  \n.  \n+  \nR  \nl  \nr  \nR  \n)  \nt  \ns  \n/lessmuch  \nr  \nf  \ni  \nc  \nr  \ne  \nf  \no  \nt  \na  \nl  \nl  \nn  \nt  \nh  \nb  \n2  \ng  \ni  \ne  \nd  \nh  \nF  \ne  \n,  \na  \na  \ns  \nt  \n/t  \n∼  \ni  \ns  \nw  \no  \ni  \nm  \na  \nr  \ne  \nConsider the Schwarzschild metric  \ny  \ns  \ne  \ne  \n-  \nt  \ny  \nd  \ne  \ni  \nl  \nn  \nGM  \na  \n1  \nd  \nτ  \nl  \nr  \nR  \ne  \na  \ni  \ni  \nr  \nR  \nw  \ns  \nt  \no  \ni  \n-  \ns  \n)  \ne  \ng  \nz  \n-  \n+  \np  \nm  \n2  \n2  \nb  \n=  \nr  \ns  \ni  \ns  \nd  \nρ  \n-  \ni  \ne  \nr  \ni  \n+  \n2  \nGM  \n3  \nds  \ne  \n)  \ne  \nS  \na  \nm  \n=  \nfl  \nt  \nds = -In Lemaitre coordinates, this takes the form ds 2 = -dτ 2 + 2 GM R dρ 2 + R 2 d Ω 2 . (2) where τ and ρ are defined by the relations dτ = ± dT -√ 2 GM R ( 1 -2 GM R ) -1 dR, (3) dρ = ∓ dT + √ R 2 GM ( 1 -2 GM R ) -dR. (4) Subtracting ( ?? ) from ( ?? ), we have R = [ 3 2 ( τ + ρ ) ] 2 / 3 (2 GM ) 1 / 3 , (5) where an integration constant has been absorbed by a shift in the origin of the expression for T as a function of ± T = 4 GM  ( )  with R given in ( ?? ). A second integration constant has been absorbed by a shift in the T coordinate. Since the metric ( ?? ) is synchronous, the lines of constant spatial coordinate are geodesics, and τ is the w i t h R S i n c e t h e m ( 1 ) ( 2 ) R ) -1 dR , ( 3 ) GM R ) -1 dR . ( 4 ) ) 1 / 3 , ( 5 ) t h e o r i g i n o f h e ρ c o o r d i n a t e . T h e  -ρ , ( 6 ) i n t h e T c o n a t e . e s i c s , h e 1 . 1 C o n s i d e r I n L e m a i t r e c o w h e r e τ a n d ρ a r e d e fi n e d d τ d ρ = ∓ S u b t r a c t i n g ( 3 ) o m ( 4 ) , w e h a v e R = [ 3 2 ( τ + ρ ) ] a n i n t e g r a t i o n c o n s t h a s b e e n a b s o r b e d b y a s n f o r T a s a f u n c t i o n τ a n d ρ c a n b e f o u n d f r o m ( ± T = 4 GM   1 3 ( R 2 GM ) 3 / 2 + √ R 2 GM w i t h s e c o n d i n t e g r a t i o n c t a n t h S i n c e t c h r o n o u s , t h e l i n e s o o I n L e w h e r e τ a n d ρ a r e S u b t i n g ( 3 ) w h e e 1 . 1 C o n s i d I n L e m a i t r e c o o r d i h e r e τ a n d ρ a r e d e fi n S u s p B a n d a d 1 . 1 M a t c hi C o n s i d e r t h e S c h I n L e m a i t r e c o o r d i n a t e s , t h i w h e r e n d ρ a r e d e fi n e d b y t h e r a c t i n g →-∞ ρ = -2 w a l l e m b e d d e d i n a s v a t o t h e d y n a m i c s o f s u c h a c c e l e r a t i o i m e -s c a l e d u e t o t h e r e p e r e p u l s i v e fi e l d c a n b e i g n o r e d . I n t h i s c a r d b y c o s m o l o g i c a l e x p a n s i o n . E v e n t u a l l y , w q u i c k l y n d e r i t s t e n s i o n a n d f o r m s a b l a c k h o l e o f r a d i u s H e r e , w i l l b e p r i m a r i l y i n t e r e s t e d i n t h e o p p o s i t e l i m i t , w h e r e t h e m a t t e r a r o u n d i t w h i l e i t s s i z e g r o w s f a s t e r t h a n t h e a m b i e n t e x l e a d s t o t h e f o r m a t i o n o f a w o r m h o l e . T h e d u s t w h i c h w a s o r i g i n a l l y i n a b a b y u n i v e r s e , a n d i n t h e a m b i e n t F R W u n i v e r s e w e a r e l e f t w i t h a s p h e r i c a l r e g i o n o f o f v a c u u m . B e f o r e w e c o n s i d e r t h e e ff e c t o f t h e d o m a i n w a l l , l e t u s fi r s t d i s c u s s t h e a n d a d u s t c o s m o l o g y . 1 . 1 M a t c hi ng Sc h w a r zs c hi l d t o a dus t co s m o l o g y C o n s i d e r t h e S c h w a r z s c h i l d m e t r i c ds 2 = -( 1 -2 GM R ) dT 2 + ( 1 -2 GM R ) -1 dR 2 + R 2 d Ω 2 . I n L e m a i t r e c o o r d i n a t e s , t h i s t a k e s t h 2 GM R d ρ 2 + R 2 d 2 . w h e r e τ a n d ρ a r d e fi n e d b y t h e r e l a d τ = ± dT ( 1 -2 GM R ) -1 dR d ρ = ∓ dT + 2 GM R ) -1 dR . S u b t r a c t i n g ( 3 ) f r o , w e h a v e R = [ 3 2 ( τ + ρ ) ] ( 5 ) w h e r e a n i n t e g r a t i o n c o n h a s b e e n a b s o r b e d b y o r i g i n o f t h e ρ o r d i n a t e . T h e e x p r e s s i o n f o r T a s a f u n c t τ a n d ρ c a n b e f o u n d f r ± T = 4 GM  ( R 2 GM ) 3 / 2 + √ R 2 GM -  -ρ , R g i v e n i n ( 5 ) . A s e c o n d i n t o n c o a n t h a s b e e n a b S i h e m e t r i c ( 2 ) i s s y n c h r o n o u s , o n s t a n t s p a t i a l c 1 Sphe r i ca l do m a i n w a l l i n dus t co R = 0 R = 2 GM , T →∞ R = 2 GM , T →-∞ ρ = -2 GM C o n s i d e r a s p h e r i c a l d o m a i n w a l l e m b e d d e d i n a s p a t i a d i ff e r e n t t i m e -s c a l e s a r e r e l e v a n t t o d y n a m i c s o f s u c h a w a l l a n d t h e o t h e r i s t h e a c c e l e r a t i o n t i m e a l e d u e t o t h e r e p u l s i v e t σ ∼ ( G σ ) -1 . F o r R /lessmuch t σ , t h e r e u l s i v e fi e l d c a n b e n o r e d . I n t h i s c a s e , f o r c o n f o r m a l l y s t r e t c h e s m o l o g i c a l e x p a n s . E v e n t u a l l y , w h e n t h e w q u i c k l y s h r i n k s u o n a n d f o r m s a b l a h o l e o f r a d i u s R S = 2 GM H e r e , w e w i l e r e s t e d i n t h e o p p o e l i m i t , w h e r e R /greatermuch t σ . I t h e m a t t e r a r o u n d g r o w s f a s t e r t h a n e a m b i e n t e x p a n s i o n r a t l e a d s t o t h e f o r m a t i o T h e d u s t w h i c h w o r i g i n a l l y i n t h e i n t e r i o r a b b y u n i v e r s e , a n d i R W u n i v e r s e w e a r e f t w i t h a b l a c k h o l e r e p h e r i c a l r e g i o n o f o f v a c u B e f o r e w e c o n s i d e r t h e e a i n w a l l , l e t u s fi r s d i s c u s s t h e m a t c h i n g o f S a n d a d u s t c o s m o l o g y . 1 M a t c hi ng Sc h w a r zs c hi l d m o l o g y C o n s h e S c h w a r z s c h i l d m e t r i c ds 2 = -( 1 -2 GM R ) GM ) -1 dR 2 + R 2 d Ω I n L e m a i t r e c o o t e s , t h i s t a k e s t h e f o r m ds 2 = -d τ 2 + 2 GM R w h e r e τ a n d ρ a r e d e fi n e d h e r e l a t i o n s d √ 2 GM R 1 S u i n g ( 3 ) f r o m ( 4 ) , w e h a ρ = -2 GM d d e d i n a s p a t i a l l y n a m i c s o f s u c h a w a l l . m s c a l e d u e t o t h e r e p u l s i v e g r e fi e l d c a n b e i g n o r e d . I n t h i s c a s e , f o r t c o s m o l o g i c a l e x p a n s i o n . E v e n t u a l l y , w h e n t h e w a n s i o n a n d f o r m s a b l a c k h o l e o f r a d i u s R S = 2 GM n t e r e s t e d i n t h e o p p o s i t e l i m i t , w h e r e R /greatermuch t σ . I n u i t s s i z e g r o w s f a s t e r t h a n t h e a m b i e n t e x p a n s i o n r a t e . e f o r m a t o f a w o r m h o l e . T h e d u s t w h i c h w a s o r i g i n a l l y i n t h e i n t e r i o r o u n i v e r s e , a n d i n t h e a m b i e n t F R W u n i v e r s e w e a r e l e f t w i t h a b l a c k h o l e r e m i c a l r e g i o n o f o f v a c u u m . r e w e c o n s i d e r t h e e ff e c t o f t h e d o m a i n w a l l , l e t u s fi r s t d i s c u s s t h e m a t c h i n g t c o s m o l o g y . 1 M a t c hi ng Sc h w a r zs c hi l d t o a dus t co s m o l o g y C o n s i d e r t h e S c h w a r z s c h i l d m e t r i c ds 2 = -( 1 -2 GM R ) dT 2 + ( 1 -2 GM R ) -1 dR 2 + R 2 . ( 1 I n L e m a i t r e c o o r d i n a t e s , t h i s t a k e s t h e f o r m ds 2 = -d τ 2 + 2 R 2 . ( 2 ) w h e r e τ a n d ρ a r e d e fi n e d b y t h e r e l a t i o n s d τ = ± dT -√ 2 GM R ( 1 -2 R ( 3 ) d ρ = dT + √ R 2 GM ( 1 -2 GM R ) ( 4 S u b t r a c t i n g ( 3 ) f r o m ( 4 ) , w e h a v e R = [ 2 ) ] 2 / 3 ( 2 GM ) 1 / 3 , ( 5 ) w h e r e a n i n t e g r a t i o n c o n s t a n t h a s b e e n a b s o r b e a s h i f t i n t h e o r i g i n o f t h e T h e p r e s s i o n f o r T a s a f u n c t i o n o f τ a n d ρ c a n b e f o u m ( 3 ) a s : ± T = 4 GM   1 3 ( R 2 GM ) 3 / 2 + √ R 2 GM n h -1 √ R 2 GM   -ρ , i n ( 5 ) . A s d i n t e g r a t i o n c o n s t a n t h a s b e e n a d h i f t i n t h e T c ( 2 ) i s s y n c o u s , t h e l i n e s o f c o n s t a n t s p a t i a l c e o d e s i c s , ρ = -2 GM ρ = d e d i n a s p a t i a l l y fl a t n a m i c s o f s u c h a w a l l . O n e i s c a l e d u e t o t h e r e p u l s i v e g r a v i t a t i fi e l d c a n b g n o r e d . I n t h i s c a s e , f o r t /lessmuch R /lessmuch m o l o g i c a l e x p a n s i o n . E v e n t u a l l y , w h e n t h e w a l l f a l l s e n s i o n a n d f o r m s a b l a c k h o l e o f r a d i u s R S = 2 GM ∼ t 2 /t σ y i n t e r e s t e d i n t h e o p p o s i t e l i m i t , w h e r e R /greatermuch t σ . I n t h i s c a s s i z e g r o w s f a s t e r t h a n t h e a m b i e n t e x p a n s i o n r a t e . A s w e o r m h o l e . T h e d u s t w h i c h w a s o r i g i n a l l y i n t h e i n t e r i o r o f t h e w a e r s e , h e a m b i e n t F R W u n i v e r s e w e a r e l e f t w i t h a b l a c k h o l e r e m n a n t c y s r e g i o n o f o a c u u m . e f o r e w e c o n s i d e r t h e e ff e c t o f t h e d o m a i n w a l l , l e t u s fi r s t d i s c u s s t h e m a t c h i n g o f S c h w a r a d u s t c o s m o l o g y . a t c hi ng Sc h w a r zs c hi l d t o a dus t co s m o l o g y i d e r t h e S c h w a r z s c h i l d m e t r i c ds 2 = -( 1 -2 GM R ) dT 2 + ( 1 -2 GM R ) -1 dR 2 + R 2 d Ω 2 . ( 1 ) I n L e m a i t r e c o o r d i n a t e s , t h i s t a k e s t h e f o r m ds 2 = -d τ 2 + 2 GM R d ρ 2 + R ( 2 ) w h e r e τ a n d ρ a r e d e fi n e d b y t h e r e l a t i o n s d τ = ± dT -√ 2 GM R ( 1 -2 R ( 3 ) d ρ = ∓ dT + √ R 2 GM ( 1 -2 GM R ) -( 4 ) S u b t r a c t i n g ( 3 ) f r o m ( 4 ) , w e h a v e R = [ + ρ ) ] 2 / 3 ( 2 GM ) 1 / 3 , ( 5 ) w h e r e a n i n t e g r a t i o n c o n s t a n t h a s b e e n a b s o y a s h i f t i n t h e o r i g i n o f t h e ρ e x p r e s s i o n f o r T a s a f u n c t i o n o f τ a n d ρ c a n b e o m ( 3 ) a s : ± T = 4 GM   1 3 ( R 2 GM ) 3 / 2 + √ 2 GM n h -1 √ R 2 GM   -ρ , t h R g i v e n i n ( 5 ) . A s e c o n d i n t e g r a t i o n c o n s t a n t h a s b e e n d b y a s h i f t i n t h e T c o o r d i n a t e t h e m e t r i c ( 2 ) i s s y n c h r o n o u s , t h e l i n e s o f c o n s t a n t s p a t i a l n a t e a r e g e o d e s i c s , a n d τ i s t h e 1 Sphe r i ca l do m a i n w a l l i n dus t co s m o l R = 0 R = 2 GM , T →∞ R = 2 GM , T →-∞ ρ = -2 GM ρ = C o n s i d e r a s p h e r i c a l d o m a i n w a l l e m b e d d e d i n a s p a t i a l l y fl d i ff e r e n t t i m e -s c a l e s a r e r e l e v a n t t o t h e d y n a m i c s o f s u c h a w a l l . O n a n d t h e o t h e r i s t h e a c c e l e r a t i o n t i m e -s c a l e e t o t h e r e p u l s i v e g r a v i t σ ∼ ( G σ ) -1 . F o r R /lessmuch t σ , t h e r e p u l s i v e fi e l d c a n b e i g n o r e I n t h i s c a s e , f o r t /lessmuch c o n f o r m a l l y s t r e t c h e d b y c o s m o g i c a l e x p a n s i o n . E n t u a l l y , w h e n t h e w a l l f q u i c k l y s h r i n k s u n d e r i t s t e n f o r m s a b l a c k h o o f r a d i u s R S = 2 GM ∼ H e r e , w e w i l l b e p r i m i n t h e o p p o s i t e l i t , w h e r e R /greatermuch t σ . I n t h i s t h e m a t t e r a r o u n d i t w h i l f a s t e r t h a n t h e a m i e n t e x p a n s i o n r a t e . A s l e a d s t o t h e f o r m a t i o n o f a w u s t w h i c h w a s o r i a l l y i n t h e i n t e r i o r o f t h a b a b y u n i v e r s e , a n d i n t h e a m i v e r s e w e a r e l e f t t h a b l a c k h o l e r e m n a n s p h e r i c a l r e i o n o f o f v a c u u m . B e f o e w e c o n s i d e r t h e e ff e c t o f l l , l e t u s fi r s t d i s c t h e m a t c h i n g o f S c h w a a n d a d u s t c o s m o l o g y . 1 . 1 M a ng Sc h w a r zs c hi l d t o a dus g y C o n s i d e r t h e S r z s c h i l d m e t r i c 2 = -( 1 -2 GM R ) dT 2 + ( dR 2 + R 2 d Ω 2 . I n L e m a i t r e c o o r d i n a t e s t a k e s t h e f o r m ds 2 = -d τ 2 + 2 GM R d ρ 2 + R w h e r e τ a n d ρ a r e d e fi n e d b y t h e n s d τ = ± ( 1 -2 GM d ρ = b t r a c t i n g ( m ( 4 ) , w h a v e R t σ ∼ ( Gσ ) -1 . For R /lessmuch t σ , the repulsive field can be ignored. In this case, for t /lessmuch R /lessmuch t σ , the domain wall is conformally stretched by cosmological expansion. Eventually, when the wall falls within the horizon, it quickly shrinks under its tension and forms a black hole of radius R S = 2 GM ∼ t 2 /t σ /lessmuch t . Here, we will be primarily interested in the opposite limit, where R /greatermuch t σ . In this case, the wall repels the matter around it while its size grows faster than the ambient expansion rate. As we shall see, this leads to the formation of a wormhole. The dust which was originally in the interior of the wall goes into a baby universe, and in the ambient FRW universe we are left with a black hole remnant cysted in a spherical region of of vacuum. Before we consider the effect of the domain wall, let us first discuss the matching of Schwarzschild and a dust cosmology. 1.1 Matching Schwarzschild to a dust cosmology Consider the Schwarzschild metric ds 2 = -( 1 -2 GM R ) dT 2 + ( 1 -2 In Lemaitre coordinates, this takes the form ds 2 = -dτ 2 + 2 GM R dρ 2 + R 2 d Ω 2 . (2) where τ and ρ are defined by the relations dτ = ± dT -√ 2 GM R ( 1 2 GM R ) -1 dR, (3) dρ = ∓ dT + √ R 2 GM ( 1 -2 GM R ) dR. (4) Subtracting (3) from (4), we have R = [ 3 2 ( τ + ρ ) ] 2 / 3 (2 GM ) 1 / 3 , (5) where an integration constant has been absorbed by a shift in the origin of the ρ coordinate. The expression for T as a function of τ and ρ can be found from (3) as: ± T = 4 GM   1 3 ( R 2 GM ) 3 / 2 + √ R 2 GM -tanh -1 √ R 2 GM   -ρ, given in (5). A second integration constant has been absorbed by a shift in the T coordinate. Since the metric (2) is synchronous, the lines of constant spatial coordinate are geodesics, and τ is the t σ ∼ ( Gσ ) -1 . For R /lessmuch t σ , the repulsive field can be ignored. In this case, for t conformally stretched by cosmological expansion. Eventually, when the quickly shrinks under its tension and forms a black hole of radius R S = Here, we will be primarily interested in the opposite limit, where R the matter around it while its size grows faster than the ambient expa leads to the formation of a wormhole. The dust which was originally in spherical region of of vacuum. Before we consider the effect of the domain wall, let us first discuss and a dust cosmology. 1.1 Matching Schwarzschild to a dust cosmology ds 2 = -dτ 2 + 2 GM R dρ 2 + R 2 d Ω 2 . where τ and ρ are defined by the relations dτ = ± dT -√ 2 GM R ( 2 GM ) -1 dρ = ∓ dT + √ R Subtracting (3) from (4), we have R = [ 3 2 ( τ + ρ ) ] 2 / where an integration constant has been absorbed by a shift in the or expression for T as a function of τ and ρ can be found from (3) as: ± T = 4 GM   1 3 ( R 2 GM ) 3 / 2 + √ 2 with R and the other is the acceleration time-scale due to the repulsive gravitational field of the domain wall t σ ∼ ( Gσ ) -1 . For R /lessmuch t σ , the repulsive field can be ignored. In this case, for t /lessmuch R /lessmuch t σ , the domain wall is conformally stretched by cosmological expansion. Eventually, when the wall falls within the horizon, it quickly shrinks under its tension and forms a black hole of radius R S = 2 GM ∼ t 2 /t σ /lessmuch t . Here, we will be primarily interested in the opposite limit, where R /greatermuch t σ . In this case, the wall repels the matter around it while its size grows faster than the ambient expansion rate. As we shall see, this leads to the formation of a wormhole. The dust which was originally in the interior of the wall goes into a baby universe, and in the ambient FRW universe we are left with a black hole remnant cysted in a spherical region of of vacuum. Before we consider the effect of the domain wall, let us first discuss the matching of Schwarzschild and a dust cosmology. 1.1 Matching Schwarzschild to a dust cosmology Consider the Schwarzschild metric ds 2 = -( 1 -2 GM R ) dT In Lemaitre coordinates, this takes the form ds 2 = -dτ 2 R where τ and ρ are defined by the relations dτ = ± dT -√ 2 GM R ( 1 -2 GM R ) -1 dR, (3) dρ = ∓ dT + √ R 2 GM ( 1 -2 GM R ) -1 dR. (4) Subtracting (3) from (4), we have R = [ 3 2 ( τ + ρ ) ] 2 / 3 (2 where an integration constant has been absorbed by a shift in the origin of the ρ coordinate. The expression for T as a function of τ and ρ can be found from (3) as: ± T = 4 GM   1 3 ( R 2 GM ) 3 / 2 + √ R 2 GM -tanh -1 √ R 2 GM   -ρ, (6) with R given in (5). A second integration constant has been absorbed by a shift in the T coordinate. Since the metric (2) is synchronous, the lines of constant spatial coordinate are geodesics, and τ is the FIG. 14: A Schwarzschild solution can be matched to a dust FRW solution at ta ρ = ρ 0 hypersurface. This can be done in two different ways. The left panel shows an expanding ball of dust matched to a Schwarzschild interior, while the right panel shows a black hole interior embedded in a dust FRW exterior. exterior solution, as indicated in the left panel, or a black hole interior can be embedded in a dust FRW exterior, as indicated in the right panel. In this case, the black hole connects through a wormhole to the asymptotically flat region IV . 2. Domain wall in Schwarzschild Let us now consider the motion of a domain wall in a spherically symmetric vacuum. In Lemaitre coordinates, a radial vector which is tangent to the worldsheet is given by u µ = ( ˙ ˆ t, ˙ ρ, 0 , 0), where a dot indicates derivative with respect to the wall\'s proper time τ : ˙ ρ = dρ dτ . (B14) The normal to the worldsheet is then given by n µ = √ 2 GM R ( -˙ ρ, ˙ ˆ t, 0 , 0 ) . (B15) Since the normal is in the radial and temporal directions, the angular components of the extrinsic curvature K µν = n µ ; ν are given by the simple expression K ΩΩ \' = 1 2 n µ ∂ µ g ΩΩ \' = ( 2 GM R ˙ ρ + ˙ ˆ t ) g ΩΩ \' R . (B16) The normalization of the radial tangent vector u µ requires that ˙ ˆ t = √ 1 + 2 GM R ˙ ρ 2 . (B17)  \nl  \nGM  \nR  \na  \ne  \ng  \ns  \nr  \n-  \ni  \nx  \ni  \ne  \nc  \nS  \nw  \nn  \nv  \nn  \ni  \nR  \nt  \ne  \ng  \nr  \ni  \nc  \ne  \ne  \nm  \nh  \ne  \nt  \nt  \n)  \nt  \n3  \nc  \nt  \nR  \no  \na  \ns  \nn  \na  \n2  \n.  \nt  \nr  \nb  \nu  \n2  \nGM  \nS  \n-  \nΩ  \n1  \nR  \ns  \nf  \nc  \n.  \n2  \nd  \ny  \na  \n2  \nb  \n(  \n6  \nd  \ndR  \ne  \nr  \nb  \n√  \nt  \ng  \ni  \nv  \nc  \ns  \ndT  \nl  \no  \ns  \nu  \no  \nt  \nco  \nl  \n1  \ndus  \nm  \n-  \nr  \ne  \ne  \nw  \nd  \nΩ  \n.  \n(1)  \n,  \nGM  \no  \ne  \nh  \nh  \n.  \nl  \ne  \n1  \nw  \n1  \nR  \ni  \nh  \nt  \nt  \n2  \nGM  \ni  \no  \nm  \nT  \nl  \n3  \n/  \n2  \n,  \ns  \ndR  \nr  \n2  \nt  \nu  \nn  \nn  \n2  \nb  \nR  \nW  \na  \nw  \no  \na  \n-  \nh  \n2  \ndR  \nn  \nf  \n2  \nm  \ns  \na  \ni  \nt  \nw  \na  \nr  \ne  \nw  \nf  \no  \n,  \n3  \n1  \ns  \nρ  \ncoordinate.  \nThe  \n-  \n3  \n(  \n(  \n/  \n/  \n)  \nT  \n.  \nh  \n(  \no  \nf  \nr  \nt  \ny  \nt  \n)  \ne  \n.  \n(  \ni  \ne  \n2  \nGM  \nt  \ni  \n2  \nGM  \nn  \na  \ni  \nv  \ne  \nn  \nR  \n5  \n2  \ni  \nA  \n)  \n1  \nh  \ne  \nt  \ni  \ns  \n2  \n-  \n)  \n2  \ni  \nn  \nt  \nc  \n(  \n1  \nds  \n=  \na  \na  \nn  \nw  \nτ  \nn  \ns  \ny  \nR  \n:  \ns  \n2  \n)  \na  \n-  \n)  \ndT  \nf  \nGM  \ni  \nt  \nr  \n)  \ns  \ne  \n2  \nR  \na  \nw  \nn  \nt  \no  \np  \ne  \ni  \nh  \nf  \ni  \ne  \nm  \nF  \nn  \n2  \nm  \ndR  \ndρ  \nt  \ni  \no  \nn  \no  \nGM  \n)  \no  \nr  \n+  \nR  \n3  \nb  \nh  \nt  \ni  \ne  \n√  \ne  \n2  \n2  \n√  \nd  \nΩ  \n.  \n(2)  \n2  \nGM  \n-  \nR  \n1  \n2  \no  \ndR  \nd  \no  \ne  \nm  \nt  \n.  \nt  \n-  \nu  \nm  \nff  \nGM  \n)  \ndR  \nh  \n=  \na  \n+  \n(  \n(  \ns  \ns  \n1  \nds  \n.  \nn  \nr  \n.  \n(1)  \nc  \nρ,  \n(6)  \nR  \n2  \nd  \no  \no  \nco  \n2  \nR  \nt  \no  \nx  \ne  \nR  \n1  \ntanh  \n-  \n√  \ne  \nIn Lemaitre coordinates, this takes the form  \n+  \nR  \nd  \nΩ  \nw  \nd  \nr  \nr  \na  \nb  \ne  \n)  \n(  \n4  \nt  \nd  \nu  \nm  \no  \nd  \n-  \n-  \n2  \nρ  \nds  \nc  \nr  \nf  \nn  \nf  \nτ  \n+  \n(  \n3  \nR  \n√  \nGM  \n2  \nGM  \nt  \ndR  \na  \ni  \nn  \ny  \nt  \nt  \ni  \nn  \na  \ns  \nt  \ns  \nr  \nd  \nd  \n-  \nb  \n,  \nt  \nh  \nn  \na  \n-  \nf  \no  \nb  \ns  \nn  \no  \nd  \nf  \nr  \nr  \nb  \ni  \no  \nn  \ne  \no  \na  \nn  \nh  \ng  \nh  \ne  \nm  \nh  \n2  \n1  \na  \ndT  \n-  \n1  \n+  \n3  \n(  \no  \nR  \nd  \n6  \no  \n(  \nρ  \nm  \no  \nt  \na  \ni  \nl  \nb  \nt  \nt  \ne  \na  \na  \ne  \nn  \nr  \nc  \n2  \nGM  \nt  \ni  \nr  \na  \nr  \nb  \n2  \nGM  \n1  \nh  \n-  \ni  \no  \n-  \nn  \nt  \na  \n.  \nt  \n1  \nf  \nn  \ns  \n(  \n2  \nR  \n2  \n2  \n.  \nd  \n+  \nd  \nρ  \n2  \nc  \ns  \no  \na  \n-  \n2  \ny  \nn  \nGM  \ne  \ns  \nb  \n1  \ni  \ne  \ne  \nl  \n√  \nh  \n=  \nb  \ns  \no  \ne  \ng  \nt  \nr  \nc  \ni  \nr  \nd  \no  \no  \n-  \ns  \n)  \no  \ng  \nr  \nΩ  \nR  \ne  \ndT  \nn  \n-  \na  \nt  \ni  \n±  \nr  \na  \nd  \no  \no  \ndT  \n-  \n3  \ne  \nR  \n2  \nGM  \nR  \n=  \n√  \nρ  \n(2  \nGM  \n)  \n-  \ndT  \ndT  \n+  \n∓  \nd  \ne  \nR  \nGM  \n1  \nGM  \n1  \no  \n+  \nR  \n2  \n∓  \n1  \n/  \n3  \nv  \ne  \n)  \n,  \ntanh  \ni  \ns  \n-  \nc  \no  \nr  \nt  \nn  \nt  \nr  \nc  \nd  \no  \ni  \no  \nt  \n(6)  \nn  \n√  \nt  \nt  \na  \nh  \n2  \nu  \nf  \nn  \ns  \ns  \na  \na  \nr  \no  \n-  \ngiven in (5). A second integration constant has been absorbed  \ne  \nh  \ni  \nw  \nn  \nc  \ne  \nu  \nS  \nt  \n)  \nf  \ne  \nSince the metric (2) is synchronous, the lines of constant spatial coordi  \nc  \ng  \nn  \nn  \ne  \ne  \ni  \ni  \no  \no  \nn  \nr  \nn  \na  \nr  \nr  \nt  \nh  \ne  \nn  \nt  \ng  \nn  \nw  \nn  \ne  \nh  \na  \ne  \ni  \na  \nr  \ns  \nr  \ne  \ni  \nf  \no  \ne  \nx  \ns  \np  \nf  \no  \ni  \ne  \nr  \nn  \ns  \ns  \no  \no  \nT  \nn  \nt  \na  \ni  \no  \nn  \nu  \nc  \nf  \ns  \nT  \na  \na  \n3  \nb  \n-  \nt  \nτ  \n3  \ny  \nr  \n1  \n)  \no  \nt  \nr  \no  \nc  \ne  \nb  \ni  \nt  \n+  \n1  \n/  \n3  \n1  \nh  \nc  \n,  \n(5)  \nn  \nτ  \n(  \n1  \n∓  \nGM  \n)  \n.  \n-  \n2  \nt  \nhi  \n1  \n1  \n.  \nc  \n(  \ns  \nm  \n2  \n/  \n3  \nco  \n(  \n2  \ni  \n]  \ndus  \nt  \nρ  \nt  \na  \no  \ns  \nm  \nh  \nf  \n(  \nT  \nGM  \nt  \n3  \ni  \nn  \n∞  \na  \n)  \nm  \n→-  \nn  \nT  \ns  \n,  \nT  \ne  \n∞  \n:  \nb  \ne  \nb  \ne  \ne  \n→-  \nh  \nl  \nGM  \nw  \nd  \nm  \nl  \nl  \nh  \n-  \n1  \n2  \nd  \nd  \ny  \ne  \na  \na  \nl  \ne  \nt  \n=  \nw  \nt  \nd  \ni  \nn  \nR  \n√  \nt  \nR  \nh  \n2  \na  \nt  \nt  \nm  \no  \no  \nt  \n-  \ne  \ne  \ni  \nm  \nn  \ns  \ni  \nGM  \no  \nt  \no  \nb  \nn  \nd  \nl  \ny  \nt  \nc  \na  \na  \ni  \nt  \ni  \nt  \ne  \nl  \nn  \nt  \ns  \no  \ne  \nr  \nn  \nt  \n-  \ne  \nr  \nh  \ne  \ne  \na  \ni  \nf  \nr  \ni  \ng  \na  \na  \ni  \na  \nr  \ns  \ns  \np  \na  \nl  \ne  \ne  \nc  \nt  \ne  \na  \nv  \nu  \nr  \ne  \no  \nl  \nh  \ne  \nd  \ni  \nv  \nl  \nc  \ns  \no  \ny  \nt  \ne  \ns  \np  \nu  \ni  \ne  \nb  \nr  \nr  \ne  \nb  \ns  \nt  \na  \ny  \ns  \nr  \ne  \nd  \nt  \ne  \nh  \nd  \nc  \ni  \n.  \nn  \ni  \nl  \n,  \nt  \ne  \ns  \np  \n1  \nr  \ne  \n-  \nh  \ne  \nd  \n)  \nu  \nn  \nσ  \nt  \nb  \ni  \nt  \nr  \nσ  \nm  \nr  \nn  \nl  \nl  \n(  \nG  \nt  \ni  \nt  \ni  \ne  \nu  \nk  \nR  \nd  \ni  \nn  \nd  \no  \ni  \nh  \nt  \ne  \nr  \nτ  \nl  \ny  \ns  \ne  \nc  \nh  \ne  \nt  \na  \nr  \ni  \nr  \ni  \nl  \nt  \ni  \nm  \nw  \nd  \ni  \ni  \nu  \nt  \ns  \nl  \ne  \nw  \ni  \ns  \nf  \nn  \n/lessmuch  \ns  \nr  \nt  \nl  \ny  \np  \ni  \ne  \nd  \nr  \no  \ns  \nw  \nu  \na  \nb  \ni  \nr  \nr  \nF  \nl  \nl  \na  \nl  \nm  \no  \ne  \ne  \nw  \no  \nh  \no  \no  \na  \nt  \nh  \ni  \nn  \nk  \ni  \nw  \ni  \nr  \nt  \nt  \ns  \nh  \nn  \nv  \nw  \ni  \no  \nn  \na  \nn  \nd  \no  \nm  \nn  \na  \nr  \nf  \no  \nf  \na  \nn  \nr  \nd  \no  \no  \nw  \nt  \nt  \ne  \ni  \nc  \nr  \no  \nt  \nt  \ne  \na  \nv  \nu  \nb  \nn  \ny  \nh  \nl  \ni  \no  \na  \ns  \nt  \nr  \nt  \np  \nd  \nc  \na  \ny  \nh  \ne  \nb  \nB  \nb  \na  \ns  \na  \na  \nh  \nr  \nd  \np  \nB  \ne  \nf  \ne  \nf  \ns  \nd  \nd  \nu  \nM  \na  \n1  \nn  \n.  \na  \n1  \nwith  \nR  \ns  \nC  \no  \nn  \n.  \n1  \ne  \nn  \ns  \na  \nm  \n4  \nc  \n(  \nfi  \nl  \nn  \nv  \nt  \ne  \n(  \nl  \nA  \ne  \n/lessmuch  \ne  \nh  \nn  \nw  \na  \nl  \ne  \no  \nc  \nτ  \nF  \nw  \nr  \nm  \nt  \nw  \n.  \nr  \ne  \nc  \n1  \n]  \ne  \nz  \nb  \nn  \nm  \ns  \no  \nm  \nd  \nl  \n∓  \nn  \nl  \nv  \n→∞  \nGM  \nn  \nd  \nt  \n)  \na  \nh  \ns  \nl  \ni  \nr  \nds  \nh  \ni  \ns  \nc  \nh  \n∓  \nt  \nh  \ni  \n-  \nc  \nIsrael\'s matching conditions for a thin domain wall of tension σ imply that the extrinsic curvature should change by the amount  \n[ K ΩΩ \' ] = -2 t σ g ΩΩ \' , (B18)  \nas we go accross the wall in the positive ρ direction. Here  \nt σ = (2 πGσ ) -1 . (B19)  \nBy continuity of the metric accross the wall, the proper radius of the wall as a function of proper time, R ( τ ), should be the same as calculated from both sides. This requires that the trajectory of the wall as seen from the \'inside\', ρ -( ˆ t -), should be glued to a mirror image of this trajectory, ρ + ( ˆ t + ). Here, the double sign ± refers to the choice of Lemaitre coordinate system in Eq. (B6), see Fig. 13, but ρ + ( ˆ t + ) and ρ -( ˆ t -) have the same functional form, which we shall simply denote by ρ ( ˆ t ). The extrinsic curvature of the domain wall will be the same from both sides, but with opposite sign. Combining (B16), (B17) and (B18), we have  \n2 GM R ˙ ρ + √ 1 + 2 GM R ˙ ρ 2 = R t σ , (B20)  \nwhich has the solution  \n˙ ρ = ( 1 -2 GM R ) -1 [ -R t σ + √ R 2 GM ( R 2 t 2 σ + 2 GM R -1 ) ] . (B21)  \nThe trajectory of a domain wall gluing two segments of the Schwarzschild solution is depicted in Fig. 15.  \nNote that for R = t σ we have ˙ ρ = 0, and it is straightforward to check that  \n˙ ρ > 0 , for R > t σ , (B22)  \n˙ ρ < 0 , for R < t σ , (B23)  \na property which will be important for the following discussion.  \nA similar calculation in the Schwarzschild coordinates gives a simpler form for the equations of motion [9],  \n˙ R 2 = -V ( R ) , (B24)  \n˙ T = ± ( 1 -2 GM R ) -1 R t σ . (B25)  \nFIG. 15: The trajectory of a domain wall between two segments of Schwarzschild solutions with the same mass M > M cr = t σ / (3 √ 3 G ).  \n<!-- image -->  \nwhere  \nV ( R ) = 1 -2 GM R -( R t σ ) 2 . (B26)  \nEq. (B24) is analogous to the equation of motion for a non-relativistic particle of zero energy in a static potential V ( R ). The maximum of this potential is at  \nR m = ( GMt 2 σ ) 1 / 3 . (B27)  \nThe value of V ( R m ) is negative (or positive) for M larger (or smaller) than a critical mass M cr given by:  \nM cr = t σ 3 √ 3 G . (B28)  \nFor M < M cr , a domain wall which starts its expansion at the white hole singularity will bounce at some R < R m and recollapse at a black hole singularity. Here we are primarily interested in the opposite case M > M cr , where the expansion of the wall radius is unbounded and continues past the maximum of the potential at R m , towards infinite size (see Fig. 15).', '3. Domain wall in dust': "As mentioned around Eq. (B22), for R > t σ the wall recedes from matter geodesics with the escape velocity, on both sides of the wall. Because of that, we can construct an exact solution where an expanding ball of matter, described by a segment of a flat matter dominated FRW solution is glued to a vacuum layer represented by a segment of the Schwarzschild solution. This in turn is glued to a second vacuum layer also described  \nby a Schwarzschild solution of the same mass on the other side of the wall, as described in Subsection B 2. Finally, the second layer is glued to an exterior flat matter dominated FRW solution. The result is represented in Fig. 6.  \n- [1] V. F. Mukhanov and G. V. Chibisov, 'Quantum Fluctuation and Nonsingular Universe. (In Russian),' JETP Lett. 33 , 532 (1981) [Pisma Zh. Eksp. Teor. Fiz. 33 , 549 (1981)].\n- [2] S. R. Coleman and F. De Luccia, 'Gravitational Effects on and of Vacuum Decay,' Phys. Rev. D 21 , 3305 (1980).\n- [3] K. M. Lee and E. J. Weinberg, 'Decay of the True Vacuum in Curved Space-time,' Phys. Rev. D 36 , 1088 (1987).\n- [4] R. Basu, A. H. Guth and A. Vilenkin, 'Quantum creation of topological defects during inflation,' Phys. Rev. D 44 , 340 (1991).\n- [5] A. Vilenkin, 'Gravitational Field of Vacuum Domain Walls,' Phys. Lett. B 133 , 177 (1983).\n- [6] J. Ipser and P. Sikivie, 'The Gravitationally Repulsive Domain Wall,' Phys. Rev. D 30 , 712 (1984).\n- [7] H. Kodama, M. Sasaki, K. Sato and K. i. Maeda, 'Fate of Wormholes Created by First Order Phase Transition in the Early Universe,' Prog. Theor. Phys. 66 (1981) 2052.\n- [8] V. A. Berezin, V. A. Kuzmin and I. I. Tkachev, 'Thin Wall Vacuum Domains Evolution,' Phys. Lett. B 120 , 91 (1983).\n- [9] S. K. Blau, E. I. Guendelman and A. H. Guth, 'The Dynamics of False Vacuum Bubbles,' Phys. Rev. D 35 , 1747 (1987).\n- [10] M. V. Medvedev and A. Loeb, 'Dynamics of Astrophysical Bubbles and Bubble-Driven Shocks: Basic Theory, Analytical Solutions and Observational Signatures,' Astrophys. J. 768 , 113 (2013) [arXiv:1212.0330 [astro-ph.HE]].\n- [11] R. D. Blandford and C. F. McKee, 'Fluid dynamics of relativistic blast waves,' Phys. Fluids 19 , 1130 (1976).\n- [12] J. Garriga and A. Vilenkin, 'Recycling universe,' Phys. Rev. D 57 , 2230 (1998) [astroph/9707292].\n- [13] A. E. Everett, 'Observational consequences of a 'domain' structure of the universe,' Phys. Rev. D 10 , 3161 (1974).  \n- [14] H. Maeda, T. Harada and B. J. Carr, 'Cosmological wormholes,' Phys. Rev. D 79 , 044034 (2009) doi:10.1103/PhysRevD.79.044034 [arXiv:0901.1153 [gr-qc]].\n- [15] J. Garriga and A. Vilenkin, 'Watchers of the multiverse,' JCAP 1305 , 037 (2013) [arXiv:1210.7540 [hep-th]].\n- [16] J. Garriga and A. Vilenkin, 'Black holes from nucleating strings,' Phys. Rev. D 47 , 3265 (1993) [hep-ph/9208212].\n- [17] N. Tanahashi and C. M. Yoo, 'Spherical Domain Wall Collapse in a Dust Universe, Class. Quant. Grav. 32 , no. 15, 155003 (2015) [arXiv:1411.7479 [gr-qc]].\n- [18] M. Y. Khlopov, 'Primordial Black Holes,' Res. Astron. Astrophys. 10 , 495 (2010) [arXiv:0801.0116 [astro-ph]].\n- [19] J. Garriga, 'Nucleation rates in flat and curved space,' Phys. Rev. D 49 , 6327 (1994) [hepph/9308280].\n- [20] B. J. Carr, K. Kohri, Y. Sendouda and J. Yokoyama, 'New cosmological constraints on primordial black holes,' Phys. Rev. D 81 , 104019 (2010) [arXiv:0912.5297 [astro-ph.CO]].\n- [21] M. Ricotti, J. P. Ostriker and K. J. Mack, 'Effect of Primordial Black Holes on the Cosmic Microwave Background and Cosmological Parameter Estimates,' Astrophys. J. 680 , 829 (2008) [arXiv:0709.0524 [astro-ph]].\n- [22] N. Duechting, 'Supermassive black holes from primordial black hole seeds,' Phys. Rev. D 70 , 064015 (2004) [astro-ph/0406260].\n- [23] P. H. Frampton, 'The Primordial Black Hole Mass Range,' arXiv:1511.08801 [gr-qc].\n- [24] K. Harigaya, M. Ibe, K. Schmitz and T. T. Yanagida, 'Cosmological Selection of MultiTeV Supersymmetry,' Phys. Lett. B 749 , 298 (2015) doi:10.1016/j.physletb.2015.07.073 [arXiv:1506.00426 [hep-ph]].\n- [25] L. Susskind, 'The Anthropic landscape of string theory,' In *Carr, Bernard (ed.): Universe or multiverse?* 247-266 [hep-th/0302219].\n- [26] P. Pani and A. Loeb, 'Tidal capture of a primordial black hole by a neutron star: implications for constraints on dark matter,' JCAP 1406 , 026 (2014) [arXiv:1401.3025 [astro-ph.CO]].\n- [27] F. Capela, M. Pshirkov and P. Tinyakov, 'A comment on 'Exclusion of the remaining mass window for primordial black holes ...', arXiv:1401.3025,' arXiv:1402.4671 [astro-ph.CO].\n- [28] A. Linde, 'A brief history of the multiverse,' arXiv:1512.01203 [hep-th].\n- [29] J. Garcia-Bellido, J. Garriga and X. Montes, 'Quasiopen inflation,' Phys. Rev. D 57 , 4669  \n- (1998) doi:10.1103/PhysRevD.57.4669 [hep-ph/9711214].\n- [30] X. Montes, 'Renormalized stress tensor in one bubble space-times,' Int. J. Theor. Phys. 38 , 3091 (1999) doi:10.1023/A:1026620502136 [gr-qc/9904022].\n- [31] B. Czech, 'A Novel Channel for Vacuum Decay,' Phys. Lett. B 713 , 331 (2012) doi:10.1016/j.physletb.2012.06.018 [arXiv:1112.1638 [hep-th]].\n- [32] B. J. Carr, 'Primordial black holes: Do they exist and are they useful?,' astro-ph/0511743.\n- [33] S. Clesse and J. Garca-Bellido, 'Massive Primordial Black Holes from Hybrid Inflation as Dark Matter and the seeds of Galaxies,' Phys. Rev. D 92 , no. 2, 023524 (2015) [arXiv:1501.07565 [astro-ph.CO]].\n- [34] J. Garcia-Bellido, A. D. Linde and D. Wands, 'Density perturbations and black hole formation in hybrid inflation,' Phys. Rev. D 54 , 6040 (1996) [astro-ph/9605094].\n- [35] D. Harlow, 'Jerusalem Lectures on Black Holes and Quantum Information,' arXiv:1409.1231 [hep-th].\n- [36] E. Farhi and A. H. Guth, 'An Obstacle to Creating a Universe in the Laboratory,' Phys. Lett. B 183 , 149 (1987).\n- [37] A. Borde, A. H. Guth and A. Vilenkin, 'Inflationary space-times are incompletein past directions,' Phys. Rev. Lett. 90 , 151301 (2003) [gr-qc/0110012].\n- [38] R. Penrose, 'Gravitational collapse and space-time singularities,' Phys. Rev. Lett. 14 , 57 (1965)."}