Brane Cosmology and the self-tuning of the cosmological constant in the presence of bulk black holes

Motivated by the holographic self-tuning proposal of the cosmological constant, we generalize and study the cosmology of brane-worlds embedded in a higher-dimensional bulk black hole geometry. We describe the equations and matching conditions in the case of flat, spherical and hyperbolic slicing of the bulk geometry and find the conditions for the existence of a static solution. We solve the equations that govern dynamical geometries in the probe brane limit and we describe in detail the resulting brane-world cosmologies. Of particular interest are the properties of solutions when the brane-world approaches the black hole horizon. In this case the geometry induced on the brane is that of de Sitter, whose entropy and temperature is related to those of the higher dimensional bulk black hole.


Introduction, Results and Outlook
One of the biggest and yet unsettled problem in theoretical physics is the explanation of the observed very small positive cosmological constant of our expanding universe or equivalently the explanation of an un-natural small vacuum energy. This is the so called cosmological constant problem [1,2,3]. There are two ways to approach it from a low energy point of view, one is using the properties of effective quantum field theory (EQFT) and the other is using the framework of general relativity (GR). In more detail using the EQFT paradigm, one can describe with much success the standard model (SM) of observed particles in the absence of gravity. If we use EQFT to compute the energy density of the vacuum, it receives large renormalisations originating from UV modes. Therefore unless an extreme fine-tuning of parameters takes place, one expects a huge energy density from such calculations. Simultaneously the vacuum energy density sources the gravitational fields that are responsible for the curvature of the universe thus it affects the macroscopic or classical (IR) physics and the evolution of the observed universe today. This is in stark conflict with the currently observed very small value of the cosmological constant.
There are various proposals for the resolution of the cosmological constant problem in the literature. In this paper we will focus on the self-tuning 1 mechanism proposed in [4] that involves the use of the brane-world scenario and holography 2 . In such a framework, the simplest model that encapsulates the relevant physics, is that of a five dimensional bulk space time that is a solution to an Einstein scalar theory. The geometry in the UV asymptotes to that of an AdS 5 space time. In this geometry we have embedded a four dimensional brane that corresponds to the observed universe with the standard model fields. Employing the holographic principle, the five dimensional bulk gravity theory is dual to a strongly coupled four dimensional quantum field theory that interacts with the four dimensional weakly coupled theory of the standard model brane [5,6,7,8,9]. The self-tuning of the cosmological constant in this framework is achieved by the coupling of the two sectors: the one that lives on the brane (which we will colloquially call the SM) and the one that is dual to the five dimensional bulk gravitational theory (which we will denote as the "hidden sector " [11]). This interaction when we work in mass scales below that of any messenger fields that mediate interactions between the two sectors, can be described by introducing effective brane potentials in front of the kinetic terms of the fields that live on the standard model brane. In the present simple example, these potentials are functions of the scalar field of the five dimensional gravitational theory.
In the setup of [4] one should impose Israel junction conditions at the radial point where the four dimensional brane is inserted and solve them together with the five dimensional bulk equations of motion. In [4] the authors searched for stable solutions with zero cosmological constant induced on the standard model brane and they found that such solutions exist and are stable under small perturbations given some assumptions about the scalar potential induced on the brane. These solutions correspond to the Poincaré invariant vacua of the theory. Moreover, the brane stabilizes at some point r 0 in the bulk of the five dimensional space-time. The main novelty of the construction of [4] is demanding regularity in the IR of the bulk geometry as a basic principle 3 and searching for a consistent solution where the brane stabilizes somewhere in the bulk space time. This is done assuming general effective potentials that couple the bulk theory with the one that lives on the brane and the construction is then found to be able to give solutions that evade the issues (singular solution in the IR) present at [6,10,12] 4 . Finally a mechanism analogous to that introduced in [13] results in an effective localisation (for some regime of scales) of the gravitational force in four dimensions as perceived by observers that can perform experiments on the brane.
In [14] the authors studied time dependent curved solutions that are not static but rather move in the radial bulk direction. In contrast to [15], the boundary of the five dimensional bulk theory is the one of [4]. Thus, these time-dependent solutions correspond to time-dependent states of the theory whose vacuum (static) solutions were found in [4]. The time dependence of the position of the brane induces a time dependence on the brane metric which leads to a cosmological evolution as perceived from the observers on the brane. Although the authors found the appropriate Israel junction conditions and bulk equations of motion, they could not solve them analytically and thus relied on a probe limit analysis where the brane does not backreact on the bulk geometry (the bulk geometry remains static). In this limit, the time dependence enters only in the radial direction r(τ ). Even taking such a limit, it is still possible to find cosmological solutions in agreement with the general principles in [16].
In the moving brane solutions, if the brane is moving towards the IR, the brane universe has a contracting FRW metric. On the other hand, moving towards the UV one finds an expanding FRW universe. Moreover, taking the UV limit where the bulk geometry asymptotes to AdS one finds the induced geometry on the brane to be de Sitter with a Hubble parameter that depends on the UV values of the potentials induced on the brane. In the IR on the other hand, there are several possibilities depending on how the bulk geometry ends. For example we can have acceptable forms of singularities or horizons.
In this paper, our main focus will be to extend these previous studies in the case where the bulk solution corresponds to a thermal state of the dual holographic QFT. In such a case there will be a regular bulk horizon cloacking the bulk singularity. We believe this is one of the most interesting options physically, since it corresponds to geometries that describe the typical high energy states of the dual strongly coupled holographic theory according to the eigenstate thermalisation hypothesis [17,18]. Thus, we expect our results to be quite robust and universal, especially in the limit where the brane approaches the regular horizon. In addition, having as a motivation a solution to the cosmological constant problem (for observers living on the brane), we would like to understand the properties of states where the total vacuum energy of the combined system is quite large. If the induced curvature on the brane can remain naturally low even in such a case, this is a strong indication that the mechanism proposed in [4] is quite robust.
In practice, we follow the steps of [14], but this time we include regular black hole geometries in the five dimensional bulk. We first derive the full backreacting bulk equations of motion together with the Israel junction conditions across the brane, for three possible slicings of the bulk black hole geometry, namely flat, spherical and hyperbolic. These equations are coupled nonlinear ODE's and one can solve them only numerically. Nevertheless existence of solutions for some parameter space can be established without having the precise explicit solutions.
In addition if we assume that the bulk is dual to a large-N strongly coupled gauge theory, it is also natural to expect that the backreaction of the small number of SM degrees of freedom can be naturally kept low and hence there is a good motivation to further analyse a probe brane limit. In this limit since the only time dependence comes from the radial position of the brane, the action that governs the dynamics on the brane reduces to that of a quantum mechanical system. Moreover, the induced metric on the brane is of the FRW type and is controlled by the blackening factor, the scale factor of the bulk geometry and the position of the brane in the radial direction.
The next step is to study the induced geometries on the brane in the IR and the UV limit. The UV asymptotics in the case of flat slicing are the same as in [14] and in the case of spherical slicing the only difference is that instead of the Poincaré patch of AdS the bulk geometry asymptotes to the global AdS space-time. In the UV limit, the induced geometry will be again de Sitter with the corresponding Hubble parameter to be governed by the UV limits of the coupling potentials induced on the brane and thus we will not repeat the analysis here. We are mainly interested in the IR analysis, which is revealed to differ substantially from that of [14], since in our case the IR geometry is that of a black hole horizon. In particular we find that the induced geometry on the brane in the IR corresponds to the Poincaré patch of dS for all the possible slicings. Moreover the Hubble parameter and the scale factor of the induced metric depend on the horizon values of the induced potentials on the brane as well as on the first derivative of the blackening factor which is related to the temperature of the bulk black hole. The most interesting result of our analysis is a relation between the induced cosmology horizon entropy and temperature with the same quantities in the bulk black hole side. We believe that this is a point that certainly deserves more study, since it is usually quite hard to assign microstates to the dS entropy, while our analysis indicates that such an understanding could arise from realising dS in a purely QFT theoretical setup (once we understand in more detail the microscopic theory that is dual to such braneworld geometries 5 ).
The structure of the paper is the following: In section 2 we provide the general formulae of the self tuning theory and in section 3 we describe the matching conditions on the brane in the case of flat, spherical and hyperbolic slicing. In section 4 we solve these equations analytically in the probe limit where the brane does not backreact on the bulk geometry. We conclude with section 5 where we present the brane-world cosmology in the IR limit which in our case is captured by the near horizon region of the higher dimensional bulk black hole. This cosmology is that of de-Sitter space (dS), and this allows a natural relation between the entropy and temperature S dS and T dS , with those of the higher dimensional bulk black hole S BH , T BH . Our paper is complemented with various appendices where supplementary material as well as more detailed calculations are presented.

The Self-tuning theory
In this section we review the self-tuning mechanism of [4] and [14]. Our bulk theory is an Einstein -Dilaton theory in d + 1-dimensions. The bulk space-time coordinates are x a ≡ (r, x µ ). Moreover, we embed in the d + 1 bulk a d-dimensional brane that is parametrized by x µ . For such a system the most general 2-derivative action is Such a setup can be based on the ideas presented in [9,11]. with, 3) S GH is the Gibbons-Hawking term at the space-time boundary (e.g. the UV boundary if the bulk is asymptotically AdS), M d−1 ≡ (16πG 5 ) −1 is the bulk Planck scale, g ab is the bulk metric, R is its associated Ricci scalar, V (ϕ) is some bulk scalar potential, γ µν is the induced metric on the brane, R (γ) is the intrinsic curvature of the brane.
The ellipsis in (2.3) represent higher derivative terms of the gravitational sector fields (ϕ, γ µν ) as well as the action of the brane-localized fields (such as the "Standard Model" (SM)). W B (ϕ), Z B (ϕ) and U B (ϕ) are scalar potentials that are localized on the brane since they are generated by the quantum corrections of the branelocalized fields [9]. For example, W B (ϕ) contains the brane vacuum energy, which takes contributions from the brane matter fields. All of W B (ϕ), Z B (ϕ) and U B (ϕ) are cutoff dependent and they scale as, where Λ is the UV cutoff of the brane physics as described here. Its origin was motivated and described in [4].

Field equations and matching conditions
The bulk equations of motion are notice that they depend only on V (ϕ). The d-dimensional brane separates the bulk into two regions. The UV region is the one that extends from the point r 0 that the d dimensional brane is located and ends on the UV boundary of the d + 1 bulk spacetime (where the volume form becomes infinite). The IR region is the one that extends from the point r 0 that the d dimensional brane is located and ends at the interior of the bulk space time where the volume form becomes zero. For the case of bulk black hole geometries that we study in this paper the IR region includes the black hole horizon and ends at the black hole singularity. We symbolise by g U V ab , g IR ab and ϕ U V , ϕ IR the solutions for the metric and scalar field in the UV and IR regions of the brane respectively. The X IR U V symbolise the "jump" of a quantity X across the brane. Then, we can express Israel's junction conditions as 1. A continuity equation for the metric and scalar field: 2. The extrinsic curvature and the normal derivative of ϕ should satisfy the following discontinuity conditions: here γ is the induced metric on the brane, K µν is the extrinsic curvature of the brane, K = γ µν K µν its trace, and n a a unit normal vector to the brane, pointing into the IR.
The discontinuity equations (2.7) for an action of the form (2.3) are given explicitly by where ϕ 0 (x µ ) is the scalar field on the brane.

Embeddings of Black Holes
In this section we present the specific form that the relations (2.8), (2.9) take in the presence of bulk black holes with non-trivial scalar field profile. In the following subsections we examine in detail the case of flat 3.1 and spherical 3.2 slicing and in Appendix D we derive the equations in the case of hyperbolic slicing.

Embedding of a Black hole with flat slicing
In this subsection we examine what are the conditions (2.8), (2.9) in the case that the bulk geometry is a black hole with flat slicing with a metric of the form

1)
f (r) is the blackening factor and e 2A(r) is the scale factor of the metric. In order to compute the extrinsic curvature K µν and its trace K it is useful to bring the metric (3.1) into an ADM form where we use as "time" the radial bulk direction and decompose the metric as with N, N µ the lapse and shift functions and γ µν the induced metric on the hypersurfaces Σ r .
Using the unit normal vector on the surface n M = (1/N, −N µ /N) 6 we define the extrinsic curvature and its trace as For a metric of the form (3.1), we have that The extrinsic curvature components using (3.3) and (3.4) are and the trace of the extrinsic curvature is .
By substituting (3.5), (3.6), (3.7) at (2.8) we find For stabilization of the brane to occur at some point r 0 in the bulk five dimensional spacetime, the conditions (3.9), (3.11) and (3.12) should be satisfied simultaneously. That is the case only if which means that the derivative of the blackening factor should be continuous before and after the location of the brane r 0 . Using the superpotential formalism described in Appendix A.1 we can rewrite the equations (3.9) -(3.12) as where ′ symbolise derivatives with respect to ϕ. Both W IR , W U V are solutions to the superpotential equation Holographic self-tuning Following [4], we use the rules of holography to fix the integration constant C IR , of W IR by demanding regularity in the IR 7 . Then the matching conditions (3.14)-(3.16) fix the integration constant C U V for the UV superpotential W U V , the position of the brane ϕ 0 in the field space and how the blackening factor f (ϕ(r)) behaves around ϕ 0 .
The integration constant C U V through the holographic dictionary is related to the VeV of the operator dual to ϕ, thus the effect of the insertion of the brane is to change the VeV comparatively to the case without brane 8 .
In order to have self-tuning of the CC and not a fine tuning, one should be able to find the 4d Minkowski space geometry on the brane for generic values of the parameters. In the case at hand, the parameters of the model are the bulk and brane potentials, which contain the 4-d vacuum energy 9 .
In more detail the process is the following, for arbitrary values of the potentials V and W B in the bulk and on the brane, the UV side of the geometry adjusts itself dynamically, for given C U V and ϕ 0 , so that the induced geometry on the brane can be that of 4-d Minkowski. We observe that for arbitrary initial conditions for W U V at ϕ 0 , the space-time at ϕ 0 connects to the same UV AdS region. Consequently, we conclude that any value of W U V gives rise to a regular geometry that satisfies the same boundary conditions.
From the boundary field theory point of view, these geometries are distinct only due to the different VeV of the operator dual to ϕ. This VeV is related to the integration constant C U V that fixes W U V . We conclude that the UV geometry selfadjusts such that it can be pasted to the regular IR solution on the brane at ϕ 0 for any value of the parameters there.

Embedding of a Black hole with spherical slicing
In this subsection we study the case of the bulk geometry to be a black hole with spherical slicing using the following ansatz for the metric where f (r) is the blackening factor, e 2A(r) is the scale factor of the metric, dΩ d−1 is the metric of the unit transverse sphere and R is the radius of the transverse sphere. This ansatz is appropriate when the dual QFT to the bulk theory is defined on an R × S d−1 geometry. In the case of five bulk dimensions the explicit metric for the unit sphere is In order to compute the extrinsic curvature and its trace it is convenient to use the ADM form of the metric (3.2) with the following identifications where χ ij are the metric components of the d − 1 dimensional unit sphere 10 .
Substituting (3.20) in (3.3) the components of the extrinsic curvature are and the trace is .
The Ricci scalar induced on the brane is The Einstein tensor induced on the brane can be written in the compact form 11 Moreover, (2.9) becomes 11 In the case of a four dimensional brane the explicit form of the Einstein tensor is For stabilization of the brane to occur at some point r 0 in the bulk five dimensional spacetime, the conditions (3.27), (3.28) and (3.29) should be satisfied simultaneously. That is the case only if Additionally from (2.3), we observe that in order to have a positive Planck scale on the brane, U B should be positive 12 and thuṡ which means that the derivative of the blackening factor should be discontinuous at the point r = r 0 where the brane is located and decreasing from the UV (r 0 +ǫ, ǫ > 0) to the IR (r 0 − ǫ, ǫ > 0). Using the superpotential formalism presented in Appendix B.1 we can rewrite the equations (3.27) -(3.29) as where the ′ are derivatives with respect to ϕ and Moreover, both W IR , W U V are solutions to the superpotential equation (3.36) The discussion presented in the subsection 3.1 about the holographic self-tuning, follows similarly here. Additionally, in this case we have one more parameter U B on the brane that is related through (3.30) to the discontinuity of the blackening factor at the point where the brane is inserted.

Embedding of a Black Hole with hyperbolic slicing
The case of a bulk black hole with hyperbolic slicing is presented in detail in the appendix D. We briefly state here, that one follows the exact same steps as for the case of spherical slicing, the only difference in the present case is that the curvature (3.24) has an overall minus sign. Due to this sign flip in the curvature, the equivalent equation to (3.29) in the hyperbolic case is The equations (3.27), (3.28) and (3.30), (3.31), as well as the holographic self-tuning procedure remain the same.
To conclude this section we comment on the possibility of addressing the full backreacting problem. The bulk equations of motion (2.4) and (2.5) supplemented by the junction conditions (2.8) and (2.9) can be solved analytically only in the static case as was done in [4]. In the present note we are interested in the induced cosmology on the brane thus the metric on the brane should be time dependent. Unfortunately, in such a case one cannnot solve the equations analytically, only numerically. Since we want to have analytic control on our results, in the following we constrain our analysis in the probe limit where the brane does not backreact on the bulk geometry.

The probe brane limit
In this section we study the brane dynamics in the probe limit. In this limit one has to solve the bulk equations of motion (2.4) and (2.5) without taking into account the backreaction of the brane in the bulk geometry 13 . This translates into accepting that the bulk geometry is smooth across the brane and we simply have to set the (2.7) equal to zero. Consequently, the bulk geometries in our examples are given by (3.1) and (3.18) respectively.
The induced action on the brane before gauge fixing, in this regime is where ξ µ are world-volume coordinates and the hat indicates induced quantities. 13 As we saw in section 2 for a self-tuning mechanism to operate, the brane should backreact on the bulk geometry. In order to take the probe limit we need a regime of parameters such that the induced action of the brane is much smaller than the bulk one, hence the brane cosmological term cannot be big. Despite this shortcoming our analysis can still offer a qualitative understanding of self-tuning cosmologies as was discussed in [14].

Flat Slicing
In this subsection we study the probe brane limit in the case where the bulk black hole geometry has a flat slicing and is described by the metric (3.1). Moreover, we work in the static gauge, ξ µ = x µ , and hence the only dynamical variable is r(x µ ) that we allow it only to depend on time. Thus the brane in the probe limit has only one dynamical degree of freedom namely r(t). The induced metric on the brane is given by where in the rest of this section a dot denotes:˙≡ d dt . The induced Ricci scalar on the brane times the square root of the determinant of the induced metric is found to be (4.3) We then foliate the induced metric with spacelike surfaces of constant time. The extrinsic curvature and its trace of these surfaces can be computed to yield Substituting then (4.4) into (4.3) we have whereγ ij = δ ij e 2A , i, j = 1, ..., 3.
One can then express the action of the brane (4.1) as where we cancelled the total derivative term with the Gibbons-Hawking term at the boundaries of the time integral. The action is then reduced to the one of a one dimensional Lagrangian mechanics problem. Using the definition for the superpotential (A.9) we can write the action (4.6) as 8) where V 3 is the spatial volume and (4.9)

The general solution for flat slicing
The Hamiltonian corresponding to the action (4.8) is conserved since the action (4.8) is explicitly time-independent. The momentum conjugate to r is and the conserved Hamiltonian is where we have reabsorbed the overall M 3 V 3 into a new one, namely E. Given E, we can solve (4.11) forṙ From the solution of (4.12) we obtainṙ 2 as a function of e 2A(r) , f (r) and ϕ(r). Moreover, by performing the following coordinate transformation in (4.2) we obtain Then, we can write the induced metric on the brane as where τ is the proper time on the brane. Finally, by introducing the new variable y we can express (4.12) in this simpler form This is a cubic equation for y, that can be solved according to the analysis of Appendix C. Consequently, we have translated the problem of finding the trajectory of the probe brane in the bulk geometry into solving (4.16) for y(r) for fixed E 14 taking into account that y(r) ≥ 1.

Spherical Slicing
In this subsection we study the probe brane limit in the case where the bulk black hole geometry has a spherical symmetry and is described by the metric (3.18). Again, we work in the static gauge, ξ µ = x µ , and hence the only dynamical variable is r(x µ ) that we allow it only to depend on time. Thus the brane in the probe limit has only one dynamical degree of freedom namely r(t). The induced metric on the brane is given by where dΩ d−1 is the metric of the transverse unit sphere. In the case of five bulk dimensions the explicit metric is The induced Ricci scalar on the brane times the square root of the determinant of the induced metric is found to be We then foliate the induced metric with surfaces of constant time. The extrinsic curvature and its trace are (4.20) Substituting (4.20) in (4.19) we have 14 E acts as an initial condition for the brane's trajectory. The second initial condition is just a shift of the initial time point.
whereγ ij = χ ij R 2 e 2A , i, j = 1, ..., 3 and the action (4.1) can be written as where we cancelled the total derivative term with the Gibbons-Hawking term of the action.
Using the definition for the superpotential (A.9) we havė and the action (4.22) can be written in a more compact form as where V 3 is the spatial volume and F is the same as in the flat slicing case, namely (4.25)

The general solution for spherical slicing
We will again use the fact that the Hamiltonian to the action (4.24) is conserved since the action (4.24) is explicitly time-independent. The momentum conjugate to r is and the conserved Hamiltonian is 27) where we have reabsorb the overall M 3 V 3 into a new one, namely E.
Given E, we can solve (4.27) forṙ From the solution of (4.28) we obtainṙ 2 as a function of e 2A(r) , f (r) and ϕ(r). Moreover, by performing the following coordinate transformation in (4.17) we obtain Then we can write the induced metric on the brane as where τ is the proper time on the brane. Finally, by introducing the new variable y we can writeṙ and the equation (4.28) simplifies to (4.33) (4.33) is a cubic equation for y, and can be solved as shown in appendix C. Consequently, the problem of finding the trajectory of the probe brane in the bulk geometry translates into solving (4.33) for y(r) for fixed E and imposing y(r) ≥ 1.

Hyperbolic Slicing
For the bulk black hole with hyperbolic slicing the procedure is exactly the same as in the case of the spherical slicing described in the subsection 4.2. The induced metric on the brane is now where dH d−1 is the metric of the transverse hyperbolic space. In the case of five bulk dimensions the explicit metric is The induced action on the brane is where now V 3 is the volume of the hyperbolic space instead of the volume of the 3-sphere. The induced metric on the brane after the coordinate change (4.29) is and the third degree algebraic equation associated with the trajectory of the brane is Essentially, the differences in the equations in the case of spherical and hyperbolic slicings come from the relative minus sign of the curvature.

The Mirage Cosmology
After we have solved (4.16) or (4.33) we can determine the location of the brane r(τ ) as a function of τ by integrating then we invert (4.40) for r(τ ) in order to write the induced scale factor on the brane a(τ ) = e A(r(τ )) , (4.41) with the corresponding induced Hubble parameter being determined by The cosmological scale factor e A(r) is monotonically increasing with r from the IR to the UV. Moreover, f (y 2 − 1) ≥ 0 always and thus the contraction or expansion of the universe is determined only by the sign of (4.42). We conclude that as the brane moves further inside the bulk (IR) the universe contracts. Conversely as the brane heads towards the boundary of the bulk spacetime (UV) the universe expands [16].

Asymptotic cosmologies
In this section we present the induced geometry on the probe brane when it is located close to the IR and the UV limit of the bulk geometry. This is the cosmology as seen by the observer on the brane. The parameter that defines the cosmology on the brane is the scale factor a(τ ).
The UV region is reached as r → −∞ and the IR region when r → r h where the black hole horizon is located.
As the brane moves towards the IR the observer perceives a contracting geometry in contrast as the brane moves in the opposite direction to the UV he/she sees an expanding universe. Since r decreases monotonically as one moves from the IR to the UV, we can use the velocityṙ to detect the direction of the motion of the brane. Ifṙ < 0 it moves towards the boundary of the bulk geometry (UV) whereas ifṙ > 0 it moves towards the horizon of the bulk black hole (IR).
Since we work in the IR and the UV limits our analysis is valid only approximately and thus in the expressions below we use " ∼ " as opposed to " = " when necessary.
In this note we omit the UV analysis since it is identical to the one in [14]. We just state their result for completeness. The brane geometry is asymptotically a dS universe, with the scale factor and the Hubble parameter given by where η = ±, a 0 is set by initial conditions, ℓ and h W,U are constants governed by ϕ, f defined in the UV region.
In the next subsections we will proceed to discuss the induced cosmology on the brane close the horizon of the bulk black hole for the cases of flat and spherical slicing.

Expansion near a flat horizon
To study the cosmology induced on the brane close to the black hole horizon we have to first expand the blackening factor f (r), the scale factor of the bulk metric A(r) and the scalar field ϕ(r). For our purposes it is enough to keep up to order one with respect to r − r h , where r h is the location of the black hole horizon. A more detail analysis can be found in appendix A.2.
where f i , A i , ϕ i are dimensionless and the relations of the constants to the physical parameters of the solution can be found in the Appendix A.2.
Substituting the expansions in (4.41), (A.9) and (4.25) and keeping either zero order or up to order one in r − r h expansion we get Close to the horizon (4.16) becomes from the appendices C and C.1 we find the solution of (5.6) to be where C > 0 is a positive constant that depends on A h , A 1 , F h , E, ℓ, f 1 . Then using (4.29) we can compute the derivative of the radial direction with respect to the cosmological time τ on the brane and by integrating it once we obtain the radial position of the brane from the black hole horizon r − r h as a function of the cosmological time τ The corresponding Hubble parameter is From (5.5) and (5.10) we conclude that the induced geometry on the braneas the brane approaches the horizon of the black hole-is that of the Poincaré patch of de Sitter with scaling factor (5.5) and the Hubble parameter (5.10). Depending on whether the brane moves towards the r h or away from r h the observer that lives on the brane sees a contracting H < 0 or an expanding H > 0 universe respectively.
Moreover, close to the black hole horizon the dS radius is l dS ∼ |H| −1 and thus the temperature and the entropy associated to the dS static patch close to the bulk black hole horizon are: where from appendix A.2 we find that A 1 = − ℓ 2 V h 2f 1 and f 1 is related to the bulk black hole temperature as T BH = f 1 e A h 4π . Thus, the de Sitter temperature (determined by H) is directly related to the bulk black hole temperature and the same holds for the dS entropy with the bulk black hole entropy. Unfortunately the relation is not simple since the constant C is a complicated function of the black hole temperature that comes after solving the cubic equation presented in Appendix C.

Expansion near a spherical horizon
To study the cosmology induced on the brane close to the black hole horizon we have to first expand the blackening factor f (r), the scale factor of the bulk metric A(r) and the scalar field ϕ(r). For our purposes it is enough to keep up to order one with respect to r − r h , where r h is the location of the black hole horizon. A more detail analysis can be found in appendix B.2.
where f i , A i , ϕ i are dimensionless and the relations of the constants to the physical parameters of the solution can be found at the appendix B.2.
Substituting the expansions in (4.41), (A.9) and (4.25) and keeping either zero order or up to order one in r − r h expansion we get Close to the horizon (4.33) can be written as from the appendices C and C.2 we find the solution of (4.38) to be whereC > 0 is a positive constant that depends on A h , A 1 , F h , E, ℓ, f 1 , R. Then using (4.29) we can compute the derivative of the radial direction with respect to the cosmological time τ on the brane and by integrating it once we obtain the radial position of the brane from the black hole horizon r − r h as a function of the cosmological time τ The corresponding Hubble parameter is From (5.15) and (5.27) we conclude that the induced geometry on the brane -as the brane approaches the horizon of the black hole-is that of a part of global de Sitter 15 with scaling factor (5.15) and the Hubble parameter (5.27). Depending if the brane moves towards the r h or moving away from r h the observer that lives on the brane sees a contracting H < 0 or an expanding H > 0 universe respectively. At late times we can go to a flat slicing. The relation becomes τ sph ∼ τ f lat and the metric takes the asymptotic form We therefore find again the relations for the temperature and entropy and R ≡ ℓ R e −A 0 . f 1 is related to the bulk black hole temperature as T BH = f 1 e A h 4π .

Expansion near a hyperbolic horizon
In the case of the hyperbolic slicing the procedure is exactly the same as in the spherical case described in subsection 5.2, the only difference stems from the relative minus sign in the curvature (3.24), thus the expansions close to the horizon (5.12) -(5.14) are the same as well as (5.15). The expansion close to the horizon of the equation (4.38) is from the appendices C and C.2 we find the solution of (4.38) to be whereĈ > 0 is a positive constant that depends on A h , A 1 , F h , E, ℓ, f 1 , R. Then using (4.29) we can compute the derivative of the radial direction with respect to the cosmological time τ on the brane and by integrating it once we obtain the radial position of the brane from the black hole horizon r − r h as a function of the cosmological time τ The corresponding Hubble parameter is From (5.15) and (5.27) we conclude that the induced geometry on the braneas the brane approaches the horizon of the black hole-is that of a part of de Sitter with hyperbolic slicing 16 with scaling factor (5.15) and the Hubble parameter (5.27). Depending if the brane moves towards the r h or moving away from r h the observer that lives on the brane sees a contracting H < 0 or an expanding H > 0 universe respectively. At late times we can go to a flat slicing. The relation becomes τ hyp ∼ τ f lat and the metric takes the asymptotic form We therefore find again the relations for the temperature and entropy and R ≡ ℓ R e −A 0 . f 1 is related to the bulk black hole temperature as T BH = f 1 e A h 4π . To conclude this section we summarize our results. Remarkably, for all the slicings (flat, spherical and hyperbolic) of the bulk black hole geometry, we find that the cosmologies induced on the probe brane both in the UV and in the IR are de Sitter with different scaling factors and Hubble parameters. Thus we conclude that the cosmological models constructed in our paper interpolate between two de Sitter geometries. This is in contrast with the cases studied in [14] where the probe brane acquires a big bang/crunch singularity when its location is close to the bulk's IR. In particular in our regular horizon bulk geometries the big bang/crunch singularity is acquired only when the brane passes through the horizon and reaches the bulk black hole singularity.

A. Black-hole ansatz with a flat slicing
In this appendix we present the equations of motion for the Einstein Dilaton action when the ansatz of the metric is a black hole with flat slicing as the one in (3.1). For completeness we rewrite here the ansatz for the metric The equations of motion for the action (2.2) are where˙is the derivative with respect to the radial coordinate. We integrate (A.2b) once to obtaiṅ where C is an integration constant. where ′ is a derivative with respect to the scalar field. For f (r) = 1, and A(r) = r/ℓ,

A.2 Expansion near a flat horizon
In this subsection we present the expansion of the blackening factor f (r), the scalar field ϕ(r) and of the scale factor A(r) close to the black hole horizon at r = r h , We have defined the coefficients f i , A i , ϕ i such that they are dimensionless. Substituting the expansions above into (A.2a)-(A.3) we obtain and with the following definitions In order for the horizon to be regular, f 1 = 0 which implies The above is a solution even in the case where the potential vanishes at the horizon V (ϕ h ) = 0. The AdS-Schwarzschild black hole is obtained in the case that the scalar potential has an extremum on the horizon V ′ (ϕ h ) = 0. In this case ϕ = constant and A is linear in r.
In the case that we expand around r = r h , the solutions to the equations of motion are governed only by three parameters. This can be seen from the equation (A.24)-(A.27) combined with (A.21) which relates f 1 with ϕ 1 . This has the consequence that when V h = 0 then from (A.24) A 1 ∼Ȧ = 0 at r = r h .
In the case that V ′ (ϕ h ) = 0 (A.24) -(A.27) become the blackening factor can be resummed as where ϕ is constant. From this equation we can observe that the location of the horizon can be at any point.

B. Black hole ansatz with a spherical slicing
In this appendix we present the equations of motion for the Einstein Dilaton action when the ansatz of the metric is a black hole with spherical slicing as the one in (3.18). For completeness we rewrite here the ansatz for the metric with dΩ 2 d−1 we denote the (d − 1)-dimensional sphere metric with radius one and R is the radius of the sphere. We are interested in the case d = 4, but our method applies for cases d > 2. The equations of motion are 6Ä(r) +φ 2 (r) = 0, (B.2a) We A more detailed analysis can be found in [19].

B.1 The superpotential formalism
As in all the previous cases, we introduce a superpotential by defining We notice that (B.2a) is satisfied automatically and (B.2b) becomes and We integrate (B.8) to get with F referring to flat slicing and S to the spherical one.

B.2 Solutions near a spherical horizon
In this subsection we present the expansion of the blackening factor f (r), the scalar field ϕ(r) and of the scale factor A(r) close to the black hole horizon at r = r h , We have defined the coefficients f i , A i , ϕ i such that they are dimensionless. Substituting the expansions above into (B.2)-(B.2c) we obtain

C. Solving the cubic equation
In this appendix, we present the possible solutions of (4.16), (4.33). For a more detailed analysis one can look at the appendix D of [14].

C.1 Solution close to the flat horizon
We can write equation (4.16) as Close to the horizon it becomes Depending on the sign of ∆ 3 we have different solution given in C, but all have the form y ∼ C (r − r h ) 1/2 (C.9) where C > 0 a positive constant that depends on A h , A 1 , F h , E, ℓ, f 1 .

C.2 Solution close to the spherical and to the hyperbolic horizon
We can write (4.33) and (4.38) as where the upper sign is for the spherical case and the lower for the hyperbolic one, with and Close to the horizon they become (C.14) Depending on the sign of ∆ 3 we have different solutions given in appendix C, but all of them have the form y s ∼C (r − r h ) 1/2 , y h ∼Ĉ (r − r h ) 1/2 (C. 15) whereC,Ĉ > 0 positive constants that depends on A h , A 1 , F h , E, ℓ, f 1 , R.

D. Embedding of a Black hole with hyperbolic slicing
We now examine the case of hyperbolic slicing. This corresponds to the ansatz where dH d−1 is the metric of the transverse sphere. In the case of five bulk dimensions the explicit metric is dH 2 3 = dψ 2 + sinh 2 ψ dθ 2 + sin 2 θ dφ 2 (D. 2) The Ricci scalar induced on the brane is The Einstein tensor on the brane is For (D.5) and (D.6) to be consistent the following equation should then be satisfied Equation (2.9) then becomes We notice that the only difference comparatively to the spherical slicing is in (D.8) where there is a minus sign in front of the second term due to the opposite sign in the curvature. The equation (D.7) remains the same since the explicit form of the transverse metric drops out from the equations.