Thin shell dynamics in Lovelock gravity

We study matching conditions for a spherically symmetric thin shell in Lovelock gravity which can be read off from the variation of the corresponding first-order action. In point of fact, the addition of Myers' boundary terms to the gravitational action eliminates the dependence on the acceleration in this functional and such that the canonical momentum appears in the surface term in the variation of the total action. This procedure leads to junction conditions given by the discontinuity of the canonical momentum defined for an evolution normal to the boundary. In particular, we correct existing results in the literature for the thin shell collapse in generic Lovelock theories, which were mistakenly drawn from an inaccurate analysis of the total derivative terms in the system.


Introduction
Junction conditions describe the behavior of physical fields across surfaces where matter density is discontinuous. A simple example is found in electromagnetism: the integration of the Maxwell equations over a pillbox, which encloses a charge/current density, leads to a jump in the electric/magnetic field.
Junction conditions also play an essential role in gravitational collapse dynamics. In the case of Snyder-Oppenheimer collapse [1], one deals with a matter-vacuum interface where the density is given by a step function defined at the star boundary. By contrast, when treating thin shells, the source is localized by means of a stress tensor with a Dirac delta distribution. This allows to develop the analysis of spherically symmetric thin shells, studied in the theory of General Relativity (GR) by W. Israel [2]. The projection of the field equations to the shell frame, upon imposing continuity of the spacetime metric, enables to express the jump as a precise combination of the extrinsic curvature and its trace, sourced by the matter on the shell.
An important application of junction conditions is braneworld scenario [3], where the dynamics of the universe itself is described as a four-dimensional brane embedded in a five-dimensional bulk spacetime. The source is an energy-momentum tensor proportional to a delta function. In order to balance the presence of this function, a discontinuity appears on the other side of the equation. Thus, the dynamics of the brane is given by the Israel matching conditions. Junction conditions have been also obtained for gravity theories which represent modifications to GR. In higher-derivative gravity, junction conditions were first studied in Refs. [4,5]. Later, it was shown that different types of junction conditions can be obtained for a general quadratic theory of gravity, depending on how singular the metric is allowed to be [6]. The treatment is carried out in a similar way as in GR, by identifying the term with the highest normal derivative.
On the other hand, matching conditions can also be derived from the variational principle worked out in terms of the discontinuity of the canonical momentum of the theory. In order to make this manifest, the example of electromagnetism is revisited.
In Minkowski space the dynamics is governed by the Maxwell action where F µν = ∂ µ A ν − ∂ ν A µ is the Faraday tensor and A µ = A µ (x) is the electromagnetic four-potential minimally coupled to a four-current density J µ . Notice also that due to the dependence on quadratic terms in the first derivative of the fields, the variational principle for the action (1.1) is satisfied by Dirichlet boundary conditions. Indeed, the variation of the above equation gives Now, consider a localized source in a sheet ∂M which, in Cartesian coordinates x µ = (t, x, y, z), is placed at z = 0. The manifold is divided in two regions M + and M − , which share the boundary ∂M with the normal vector n µ = δ z µ and current density J µ = j µ δ(z). In this treatment, it is implicit the assumption that the transversal components of the gauge connection A i (where the Latin indices label the coordinates x i = (x, y, t)) are continuous across the shell. As the current density is given by a delta function, the integration in the normal direction can be performed, i.e., moving the source to the interface between the two regions. The variation of the action leads, on-shell, to a surface term of the form where [F ] = F + − F − refers to the difference of a quantity F between M + and M − .
Since δA i is arbitrary at the interface, the variational principle δI = 0 implies the matching condition It is straightforward to see that the last equation reproduce the known junction condi- for a generic normal vector to the boundary. Furthermore, for an evolution of the system along the normal coordinate z the associated canonical momentum is From this simple example, one can understand the analysis of junction conditions as coming from the discontinuity of the canonical momentum (associated to the normal evolution). 1 It is an appealing idea to extend this procedure to obtain junction conditions for GR.
However, in contrast with the Maxwell Lagrangian which is quadratic in first derivative terms, the Einstein-Hilbert (EH) Lagrangian depends on second order derivatives.
As a consequence, the variation of the EH action involves variations of the metric and its first derivative. In order to have a well-posed variational principle for Dirichlet boundary conditions, the term which is linear in the normal acceleration must be eliminated by the addition of the Gibbons-Hawking-York (GHY) boundary term. As a result, a first-order Lagrangian density is obtained and can be used to derive matching conditions for thin shells in gravity. This statement will be revised in the following sections. For these reasons, in this work we focus on Lovelock gravity: the most general covariant theory in a D-dimensional spacetime giving second order divergenceless field equations [7]. In presence of a boundary, akin to General Relativity, the variation the Lovelock actions produce boundary terms that depend on variations of the metric and its first derivatives. Consequently, a well-posed variational principle -for a Dirichlet condition on the boundary metric-requires the Lovelock action to be supplemented by Myers terms [9].
The above idea was used in Ref. [10] (see also Refs. [11,12,13,14]) to derive junction conditions for an Einstein-Gauss-Bonnet (EGB) braneworld. Knowing that Israel matching conditions can be derived by varying the EH action with the GHY term, they obtained generalized junction conditions for EGB theory from the variational princi- 1 In a more contemporary context, in axion electrodynamics, where the Maxwell Lagrangian is augmented by a pseudo-scalar field coupled to the Pontryagin term for U (1), the corresponding junction conditions can also be derived in this way [8].
ple with the corresponding Myers/Müller-Hoissen term [15]. Later, in Ref. [16], it was shown that the variation of a Lovelock density plus its corresponding Myers term reproduces the Hamiltonian variation made in Ref. [17]. In particular, the canonical momentum can be easily read off from the surface term of the variation of the Dirichlet action (see also Chapter 15 of Ref. [21]). It is then clear that the canonical momentum plays the role of a generalized Brown-York stress-tensor [18] in Lovelock gravity.
Therefore, it makes sense to study generalized junction conditions for a generic Lovelock action as the discontinuity in the canonical momenta where the brane itself is a boundary . This problem will be addressed in detail in this work.
As stressed in Ref. [22], the Dirichlet action is a first-order functional whose variation coincides with the one of the Hamiltonian action [17]. As a matter of fact, both actions produce a surface term of the form π ij δh ij . As the corresponding first-order Lagrangian L 1 is the one which enters in the definition of Hamiltonian, it is the proper function to which the canonical momentum is associated to.
As we shall see below, the variational principle and the use of adapted coordinates on the shell, allows to work out the shell dynamics in Lovelock gravity, as coming from the jump in the canonical momentum across it. In particular, we correct existing results in the literature, which are a consequence of a wrong analysis of junction conditions in Lovelock theory [30].
This work is organized as follows: In section 2, we briefly review the Dirichlet problem for Einstein-Hilbert, Einstein-Gauss-Bonnet and Lovelock theories in the Gaussnormal coordinate frame. In the process, the connection between first-order Lagrangians and the Dirichlet problem is revisited. With this in mind, we show in section 3 how the variational principle gives rise to junction conditions for a thin shell. Then, in section 4, we carry out the procedure described above to obtain explicit expressions for the junction conditions for Lovelock gravity. In the Conclusions, we summarize our results and explore future directions.

The Dirichlet action for Lovelock gravity
In what follows, we make extensive use of Gauss-normal coordinates where the spacetime is foliated as an infinite series of surfaces of constant w with normal n A = Nδ w A , each one with an induced metric h ab . In this frame, we consider the splitting of the spacetime indices as A = (w, a), such that the different components of the Riemann curvature tensor are given by the Gauss-Codazzi-Mainardi (GCM) relations in Appendix B. Furthermore, the normal derivative ∂ n , defined for the vector which produces the spacetime foliation (2.1), is expressed as

Einstein gravity
In General Relativity, the dynamics is described by the Einstein-Hilbert action in D = d + 1 dimensions where R is the Ricci scalar, g = det g AB is the metric determinant and Λ is the cosmological constant. A well-posed action principle is achieved if the Einstein-Hilbert action is supplemented by the Gibbons-Hawking-York (GHY) term [23,24] where K is the trace of the extrinsic curvature K ab = − 1 2N ∂ w h ab , where a prime stands for ∂ w . Arbitrary variations of this action give where the bulk expression is E AB = G AB +Λg AB -with G AB = R AB − 1 2 R g AB the Einstein tensor-and the boundary stress tensor, which can be read off from the surface term in the above equation, is The Dirichlet action (2.4) yields a well-defined variational problem, as it remains stationary when the metric is kept fixed at the boundary. As a matter of fact, it is manifest that its variation (2.5) vanishes when imposing the boundary condition δh ab = 0 at ∂M.
One can use Gaussian coordinates (2.1) and the GCM relations (B.3) in order to lift the GHY boundary term to the bulk, i.e., to express it as a part of the bulk Lagrangian.
In doing so, one can show that the Lagrangian density associated to the action (2.4) is equal to where the bar stands for a quantity computed with the boundary metric. This is the time-honored, first-order Arnowitt, Deser and Misner (ADM) Lagrangian of Einstein gravity [25], essential ingredient to construct a Hamiltonian in Einstein gravity. As a direct consequence of this fact, the tensor (2.6) can be identified with the canonical momentum, conjugate to the dynamic variable h ab , i.e., . (2.8) One may recover the full generality of the procedure, beyond the particular gauge choice (2.1). Indeed, in the general ADM decomposition, the extrinsic curvature K ab is extended to be Upon integrating by parts, a term proportional to N a∇ b π b a is added to the Lagrangian L 1 considered above. It can be explicitly checked that variations with respect to N, N a , and h ab of L 1 yield, respectively, the constraints E w w = 0 and E w a = 0, and the dynamical where ∂ n is given by Eq.(2.2).
For the above action, endowed with the Gibbons-Hawking boundary term, Einstein equation (2.11) implies the conservation of the canonical momentum (2.25). For a given set of Killing vectors {ξ i } the corresponding conserved quantities were derived in Ref. [18], which are now broadly known as Brown-York charges Here, u a is the normal to the co-dimension 2 surface Σ and σ, the determinant of its metric, typically expressed in Schwarzschild-like coordinates. As this is a backgrounddependent notion of energy and other conserved charges, the subscript [0] in the momentum denotes evaluation on the corresponding vacuum geometry, e.g., Minkowski or global de Sitter (dS) or anti-de Sitter (AdS) spacetimes.

Einstein-Gauss-Bonnet gravity
An arbitrary modification of General Relativity to include quadratic-curvature couplings in the action leads, in general, to fourth-order field equations. On the contrary, only for a precise combination of such terms -known as Gauss-Bonnet-the EOM are still of second order in derivatives of the metric. The action for Einstein-Gauss-Bonnet gravity is where α is the Gauss-Bonnet coupling. In order to turn this action compatible with Dirichlet boundary conditions, it is necessary to supplement it with a generalization of the Gibbons-Hawking term, i.e., .(2.14) For conventions on generalized Kronecker deltas of higher rank, see Appendix A. The variation of the above action takes the form where, in this case, the equations of motion are (2. 16) and the tensor π ab in the boundary term is (2.17) It is worthwhile to notice that the index structure of the part of π a b associated to the Gauss-Bonnet term makes apparent the fact the canonical momentum vanishes identically in D = 4. Therefore, it is clear that the Dirichlet problem for the metric h ab cannot be defined in the case of the addition of the Gauss-Bonnet term to the four-dimensional gravitational action 2 .
As in the GR case, the use of Gaussian coordinates and the GCM relations allows to lift the boundary term to the whole spacetime, such that it now appears as a bulk term.
In doing so, the second normal derivatives of the metric are eliminated from the bulk Lagrangian density [22], and the action (2.14) is then written as This is a first-order Lagrangian density, which depends on the first normal derivative of the induced metric and allows to identify the tensor (2.17) as the canonical momentum conjugate to h ab .
The components E w w and E w a of the equations of motion for Einstein-Gauss-Bonnet gravity take the form In addition, the component E a b is given by As in the case of General Relativity, the field equation E w a = 0 implies a conservation law written in terms of the tensor π ab . Notwithstanding the foregoing, this tensorial quantity is identified with the canonical momentum only if the corresponding Lagrangian for Einstein-Gauss-Bonnet gravity is of first-order in normal derivatives.
In adapted, Gaussian coordinates, the derivation of the canonical momentum in GR makes it equivalent to the Brown-York stress tensor [18]. The extension of the notion of Brown-York energy-momentum tensor to EGB gravity is then naturally realized by Eq.(2.17). However, the fact π ab correctly accounts for the the energy of, e.g., Boulware-Deser black holes [27], would depend on the value of the cosmological constant. For asymptotically flat solutions, the mass obtained is the correct one. In turn, for asymptotically AdS black holes, only a fraction of the mass is obtained, with a factor which depends both on the dimension and the GB coupling. The addition of local counterterms at the boundary happens to correct the factor and to remove infrared divergences at radial infinity [31,20].

Lovelock gravity
Lovelock gravity is the natural generalization of Einstein theory for D > 4. Evidence hinting at the appearance of such terms have been found in the low energy effective action of heterotic string theory [33,34,35] and six-dimensional Calabi-Yau compactifications of M-theory [36]. This fact brought a thirst for solutions of this model. Besides Boulware-Deser black hole solution, Wiltshire solved the EGB system with the Maxwell term included [37]. In turn, in Ref. [38] the Born-Infeld model was studied.
By dropping the spherically symmetric condition, topological black holes were found [39], featuring flat and hyperbolic transversal sections. Beyond the Gauss-Bonnet term, spherically symmetric black holes were discovered in maximally degenerate Lovelock gravity [41]. In Ref. [28], such solutions were extended to Lovelock Unique Vacuum (LUV) theories with intermediate multiplicity. Also, black holes were explored for Lagrangians containing a single Lovelock term plus cosmological constant [40].
The goal of this section is to review the construction of the first-order Lagrangian for Lovelock gravity.
Consider the Dirichlet action of a generic Lovelock theory where ⌊ · ⌋ is the floor function and α p is a set of arbitrary coupling constants. The term of degree p in the curvature is that has the property to be topological in D = 2p dimensions. Its corresponding Myers term it is such that, added on top of the bulk Lagrangian, guarantees a well-posed Dirichlet principle for the boundary metric h ab .
Thus, the variation of Eq.(2.19) reads where the equation of motion is the linear combination of the individual contributions coming from each term in the Lovelock series, that is, In doing so, the p-th term in the series produces the covariantly conserved tensor The surface term in Eq.(2.22) is proportional to the canonical momentum of the theory. This resembles standard derivations in Classical Mechanics, as the Myers term is responsible for turning the gravity action into a first-order functional [22], and therefore the momentum can be readily read off as the conjugate to the metric h ab as where, for the p-th Lovelock density, one obtains the associated piece for the canonical momentum π ab with L (p) 1 given by .
(2.28) Therefore, the Dirichlet action has been written explicitly as a first-order Lagrangian density. This allows to identify the tensor (2.25) as the canonical momentum conjugate to h ab .
The E w (p)w and E w (p)a components of the field equations for the Gaussian frame are In addition, the components E a (p)b are given by (2.30) One can also recognize in (2.29) that the second-derivative terms of the equations of motion are packed as derivatives of π ab .
The addition of surface terms and its role to describe properly the physical changes across an interface is specially relevant for the next sections. The presence of a shell which is itself a boundary to the spacetime geometry will provide the link between first-order Lagrangian for gravity and the collapsing shell dynamics.

Shell co-moving frame
In this section, we succinctly review the kinematic description of collapsing thin shells as it appears in, e.g., Refs. [26,29]. In turn, the collapse dynamics would appear from the junction conditions expressed in terms of the discontinuity of the canonical momentum π ab .
Consider a manifold M, separated in two regions by a spherical thin shell located at r = R s . As the collapse is radial, the inner (−) and outer (+) regions are described by the static spherically symmetric ansatz where t ± are the corresponding time coordinates and dΩ 2 = ω mn dx m dx n is the line element of a unit sphere S D−2 .
On the other hand, one may take a coordinate system x A = {λ, x a } of the type adapted to the shell frame, such that the normal direction λ is generated by the vector For an observer moving with the shell, the induced metric has the form with x a = {τ, θ 1 , ..., θ D−2 }.
One may parametrize the shell position using its proper time τ , in terms of the Schwarzschild-like coordinates r = R s (τ ) and t ± = t ± (τ ). Indeed, gluing the line elements Eq.(3.2) and Eq.(3.1), one obtains the relation where the dot stands for ∂ τ . In the following analysis, we drop the subscript ± in the metric function, as the treatment applies equally to either the interior or exterior regions.
Of particular usefulness is the definition of the factor which reduces to the standard relativistic factor when the spacetime is Minkowski.
The normal to the shell is a space-like unit vector which, as described in the coordinate system of the static black hole geometry n µ = ∂x A ∂y µ n A , takes the form By definition, the extrinsic curvature of the shell geometry is given by [pq] . (3.11) The energy-momentum tensor sourcing the discontinuity of the geometry is given by where L M is a generic matter Lagrangian density. The matter distribution is such that it appears localized on the shell.
The interface that divides the spacetime is given by a physical thin shell, whose stress tensor is the one of a perfect fluid In the co-moving frame, where the observer is at rest respect to the shell the relevant component of the stress tensor is This relation, combined with the corresponding equation of state for the shell matter will determine the collapse dynamics, once the junction conditions are imposed.

Discontinuity of the canonical momentum
As rendered manifest by the previous discussion, the use of Schwarzschild-like coordinates is suitable to describe both the interior and exterior regions in the present collapse setup. Of course, this metric can be put in a Gauss-normal form for a radial foliation of the spacetime. In turn, in the shell co-moving frame, the normal direction to the shell is λ which mixes up the radial and time Schwarzchild coordinates. Thus, junction conditions will be derived from a variational principle in adapted Gauss normal frame.
For any Lovelock theory written in first-order form (2.28) and taking the shell as a boundary, the variation of the action is given by 4 provided that the equations of motion hold. The integration in an infinitesimal interval (−ε, +ε) across the normal direction yields where one defines the limit by the left and by the right of h ab as Then, the (Riemannian) condition of g being smooth is relaxed and one considers instead a metric which is continuous at the shell location, that is, Localizing the energy-momentum tensor with a delta function at the shell, as given by Eq. (3.15), renders Eq.(3.17) into the following form Since variations of the metric δh ab are arbitrary, the variational principle requires π ab = 8πGS ab . (3.20) In this way, the condition δI +δI M = 0 is met by imposing the junction conditions at λ = λ(R s ), and the standard Dirichlet condition on the metric at radial infinity. In turn, the junction condition (3.20) describes the motion of the collapsing shell hypersurface.

General Relativity
As a warm up exercise, one may consider the gravitational collapse of thin shells in Einstein theory of gravity with negative cosmological constant.
As the matter density is given by a delta function the geometry 3.1 will jump across the shell characterized by the discontinuity in the metric function where Ω D−2 is the volume of the unit sphere S D−2 . The difference in the mass between the exterior and interior black holes is sourced by the (τ, τ ) component of the stress tensor, what leads to Thus, the junction condition (3.20) for General Relativity reduces to the relation i.e., in agreement with Refs. [2,29,14].
One can work out an alternative form, multiplying by γ + + γ − . This is particularly convenient in the case of incoherent dust (no pressure). Indeed, in that situation, one obtains the difference between the inner and outer mass

The Einstein-Gauss-Bonnet gravity
In the Lovelock series, the next gravity theory to be considered is Gauss-Bonnet. As one is interested in this term as a correction to GR appearing in the bulk action, its addition define Einstein-Gauss-Bonnet gravity. For the collapse of thin shells with spherical symmetry, the solution with the same symmetry is given by Boulware-Deser one [27] f 2 ± (r) = 1 + 4) and σ = ±1, which represents two branches of the theory.
For given values of the GB coupling α and the cosmological constant Λ, the theory may feature black hole solutions. As for the present treatment, one may assume the existence of such black holes. The results then are is given in terms of the function γ, defined in Eq.(3.5), and they equally apply for any function f 2 (r).
The component (τ, τ ) of the canonical momentum has the form such that the junction condition can then be written as This expression governs the shell dynamics in EGB gravity, and it consistently reproduces the results in the existing literature [13,14]. In the case the theory possesses a unique vacuum, which corresponds to the particular value for the GB couplingα = ℓ 2 4 , the above expression reduces to the relation The metric function for the black solutions around this critical point in the parametric space features a distinctive asymptotic behavior in the mass term, with a much slower falloff for large r. Thus, even though it is neither possible to solve explicitlyṘ 2 s nor to find a closed expression for ∆M, it is expected that the dynamics will be radically different for that case.

General Lovelock case
For a generic Lovelock theory of gravity, the canonical momentum is the linear combination of the corresponding term associated to every p−th term of the series with the contribution due to the p−th term in the Lovelock Lagrangian given by As a consequence, the junction condition for Lovelock gravity is expressed in terms of the following discontinuity A particular choice of the set of couplings in the Lovelock series leads to the socalled Lovelock Unique Vacum (LUV) theory [28]. This choice intends to free highercurvature gravity from undesirable instabilities that may trigger transitions between different vacua. The price to pay is that now global AdS space is a zero of the field equations with multiplicity k. Therefore, this class of gravity theories do not accept a linearization around AdS background, such that their black holes do not have the asymptotic behavior of the Schwarzschild solution.
As for the momentum tensor, the generic formula (4.9) turns into In this case, the sum can be factorized using the binomial expansion, such that the junction condition (3.20) for LUV theories can be cast in the form As the relativistic factor γ + and γ − depend onṘ 2 s , M + and M − , it turns a difficult task to properly isolate the shell velocity and, therefore, to determine the exact collapse dynamics from the initial conditions.

Comparison to existing literature
In the ref. [30] a different path to the derivation of the junction conditions in Lovelock gravity was taken. The procedure followed by these authors considers -as a starting point-the equations of motion (2.24), written down in the static black hole ansatz (3.1) in Schwarzschild-like coordinates (t, r), where the canonical momentum is given by The point is that, for the specific ansatz taken, the component (t, t) of the field equations can be further simplified as single derivative term (4.14) Once this total derivative was obtained, the idea in Ref. [30] was to perform a subsequent integration across the shell position, along the normal direction λ. This integration requires a projection between the coordinate systems (t, r) and (τ, λ) (a sort of rotation) which gives rise to additional factors depending on γ. Because, in this ansatz, the terms on top of ∂ n π t (p)t in Eq.(4.13) accidentally contribute to a total derivative, the matching conditions for the shell in Ref. [30] are mistakenly proportional to ∆M. Conversely, the canonical momentum cannot account -only by itself-for the difference in mass between the interior and the exterior regions, as the energy of the system in Lovelock gravity strongly depends on the multiplicity/degeneracy of the vacuum state [43,44].
An unambiguous way to rederive the junction conditions would require writing down the corresponding field equation in an schematic form Here, the bounded terms are contributions which go to zero after one performs the integration across the shell and takes the zero thickness limit, regardless the particular ansatz used. As a matter of fact, the above relation is an efficient form to pack the second-derivative terms (in the radial direction) as normal derivatives of the momentum. This is consistent with the picture developed by Deruelle et al. in Ref. [6], where junction conditions in higher-derivative gravity can be readily obtained from the highest-derivate term.

Conclusions
In this work, we have exploited the connection between junction conditions and the variational principle for a Dirichlet boundary condition on the metric. As a direct consequence of this treatment, we have studied the thin shell collapse in Lovelock gravity as coming from the discontinuity in the canonical momentum of the theory.
After reviewing, from the above standpoint, the derivation of shell dynamics in General Relativity and Einstein-Gauss-Bonnet gravity, we work out the corresponding expression for an arbitrary Lovelock theory. These derivations stress the link between the junction conditions and the discontinuity of the canonical momentum, which is properly identified once the Lagrangian adopts a first-order form.
In General Relativity, a spherically symmetry ansatz leads necessarily to the formation of a Schwarzschild solution, as dictated by the Birkhoff theorem. The only parameter accounting for global properties of the geometry is the black hole mass. Junction conditions in Einstein gravity sees the jump between the inner and outer black hole mass as proportional to the mass of the shell [2,29]. Later work derives a similar relation for LUV gravity theories [30]. In bold contrast to these results, the shell dynamics developed here involves nonlinear combinations of the γ factors. Therefore, the erroneous reasoning in ref. [30] is assuming that junction conditions can be readily read off from the total derivative term in the field equations for the static, spherically symmetric ansatz. Because the total derivative term acquires extra, accidental contributions on top of π ab , the picture of the junctions conditions as associated to the canonical momentum is lost.
As a prospect, it would be interesting to explore the physical implications of gravitational collapse of thin shells within a holographic framework, equipped with the tools presented here. In Ref. [32], the authors study thermalization in a boundary Conformal Field Theory, where the bulk spacetime is a solution to AdS gravity. In particular, this provides a dual gravitational setup in order to work out the time evolution of entanglement entropy of the boundary CFT. It is indeed a quite appealing idea to think of a similar model where higher curvature terms are included in the bulk gravity action [42]. As discussed in the previous section, the correct identification of the discontinuous quantities in the bulk geometry may turn of key importance when it comes to a proper holographic description of the system.

A Kronecker delta of rank p
The totally-antisymmetric Kronecker delta of rank p is defined as the determinant A contraction of k ≤ p indices in the Kronecker delta of rank p produces a delta of rank where N is the range of indices.

B Gauss-normal coordinates and Gauss-Codazzi-Mainardi relations
For the spacelike foliation (2.1), the components of the Christoffel symbol are where K ab is the extrinsic curvature The curvature tensors are given by In the notation used here, the prime stands for a partial derivative in w,∇ is the covariant derivative defined with the connection Γ a bc (h) associated to the boundary metric, andR abcd is the boundary Riemann tensor. Boundary indices are raised or lowered with the metric h ab . Then, the spacetime Ricci scalar and Ricci tensor can be expressed as In turn, for the static black hole ansatz (3.1), the nonvanishing components of the Christoffel symbol, in Schwarzschild coordinates, are

C Co-moving frame coordinates and continuity conditions
One can always introduce a set of adapted Gaussian coordinates x A = {λ, x a } = {λ, τ, θ (m) } and refer them to the Schwarzschild-like coordinates of the static interior/exterior geometry, as discussed in Ref. [32]. This is done by considering the change of coordinates in the bulk r = r (λ, τ ) and t = t (λ, τ ) . (C.1) In doing so, the metric (3.1) in the proper coordinates is given by Matching the above line element to the one in the adapted coordinates ds 2 = dλ 2 + h ab dx a dx b = dλ 2 − dτ 2 + R 2 s (λ, τ ) dΩ 2 , (C. 3) leads to the conditions It is straightforward to check that this treatment recovers the shell metric Eq.(3.2) for a constant λ. Therefore, without loss of generality, the shell is located at λ = 0.