The large D Membrane Paradigm For Einstein-Gauss-Bonnet Gravity

We find the equations of motion of membranes dual to the black holes in Einstein-Gauss-Bonnet (EGB) gravity to leading order in 1/D in the large D regime. We also find the metric solutions to the EGB equations to first subleading order in 1/D in terms of membrane variables. We propose a world volume stress tensor for the membrane whose conservation equations are equivalent to the leading order membrane equations. We also work out the light quasi-normal mode spectrum of static black holes in EGB gravity from the linearised fluctuations of static, round membranes. Also, the effective equations for stationary black holes and the spectrum of linearised spectrum about black string configurations has been obtained using the membrane equation for EGB gravity.All our results are worked out to linear order in the Gauss-Bonnet parameter.


Introduction
In the limit of spacetime dimensions (D) going to infinity, the late time dynamics (at time scales of order O(1)) of black holes consists of only a finite number of degrees of freedom. The primary evidence for this was presented in [1], where it was shown that the spectrum of light quasi-normal modes (with frequency of order O(1)) consists of two scalar and one vector mode at each value of angular momentum number of the background SO(D − 1) isometry group. There are an infinitely large number of heavy quasi-normal modes of frequency of the order O(D) which die off in time intervals of the order of O(1/D). Motivated by this the authors of [2] derived a set of effective non-linear equations determining this late time dynamics in terms of variables of a membrane propagating in flat spacetime. This duality between black hole and membranes propagating in flat spacetimes is known as the large D membrane paradigm. Details of this duality in Einstein-Maxwel system [3] and in presence of a cosmological constant [4,5] have also been worked out. In [6], a systematic algorithm was proposed to work out the membrane dynamics to to any order in 1/D and explicit computations to find the first subleading correction to membrane equations and next to leading order corrections in metric were presented 1 .
Since the large D limit has been very useful in understanding many aspects of black hole physics in two-derivative Einstein gravity, it is natural to apply this methodology to study black hole dynamics in higher derivative theories of gravity. Some initial progress in this direction has already been made for the simplest higher derivative theories of gravity namely the Einstein-Gauss-Bonnet gravity. In [19] the quasi-normal mode spectrum of static black holes in EGB gravity in the large D limit has been worked out. It was found that when the Gauss-Bonnet (GB) terms can be treated as perturbations over the Einstein-Hilbert terms in the equations of motion, the same number of characteristic light modes mentioned above are still present and the frequencies of the light modes are corrected by the presence of the GB terms. In addition the effective equations for static and stationary (slowly rotating) black holes has also been obtained in [20,21]. Some specific dynamical situations involving black strings has been worked out in [22] and the Hydro-elastic complementarity study of [14] was extended to Einstein-Gauss-Bonnet (EGB) gravity in [23].
In this paper we work out the details of the large D membrane paradigm for black holes in EGB gravity to leading order in 1/D. We work to linear order in the GB parameter (β) only. The degrees of freedom of the dual membrane are same as those for the membranes dual to black holes in two derivative Einstein gravity, namely its shape and a time-like velocity vector field defined on the membrane. We find the equation of motion of the membrane dual to the black holes in EGB gravity to leading order in a 1/D expansion. It is given by Here u is unit normalised time like velocity vector field on the membrane and K is the trace of extrinsic curvature K AB of the membrane embedded in flat spacetime. We also work out the metric of the dynamical black hole which solves the EGB equation to first subleading order in 1/D in terms of membrane variables(see equation (9.1) and (9.2)). Also, we write down an expression for the world volume stress tensor of the dual membrane given by The equations of conservation of this stress tensor in the world volume of the membrane embedded in flat spacetime is equivalent to the membrane equations written above. We also find the spectrum of light quasi-normal modes of static black holes in EGB gravity in the large D limit from the linearised fluctuations of the membrane about a static spherical configuration. The frequency of the scalar modes is given by and the frequency of the velocity mode is given by Here, l denotes the angular momentum number of the modes for the spherical harmonics in the respective tensor sectors with respect to the SO(D − 1) isometry direction of the background.
2 EGB gravity: lagrangian, equations of motion and large D scaling of the GB parameter

The lagrangian
The Einstein-Gauss Bonnet (EGB) gravity is the simplest in a series of consistent 2 higher derivative terms that can be added to the Einstein-Hilbert term. They are generically called the Lovelock theories of gravity. The lagrangian for Lovelock theories is given by and R αrβr µr νr is the Riemann tensor. (2.1) The square bracket above denotes anti-symmetrisation in the indices without any numerical factor. The leading term in the series (i.e. when n = 1) is the Einstein-Hilbert term. We will be working with only the first two terms in this series which together form the EGB gravity lagrangian. The explicit form of the EGB lagrangian is given by 3 .
We will be thinking of α as a perturbatively small parameter and in this paper we will be working to linear order in α only. From the structure of the lagrangian (2.1) we see that the parameters α n count the number of derivatives on the metric for the particular term in the lagrangian. Hence, effects linear in α come from 4 derivative terms in the action and so on. If we want to consider effects which are quadratic or higher order in α (which are in principal 8 derivative effects), we feel that it is necessary to also take into account other terms in the lagrangian which explicitly contain higher number of derivatives on the metric. Hence to unambiguously keep track of four derivative terms in the lagrangian we only work to linear order in α.
2 By consistent we mean there are no -super-luminal modes for linearised fluctuations about flat spacetime. For fluctuations about other spacetimes, some constraint on the parameters can be imposed to make sure there are no ghost propagating modes 3 Here without loss of generality we can set α 1 = 1 and α 2 = α

Equation of motion and large D scaling of the GB parameter
The variation of the EGB action w.r.t. the metric gives rise to the following equation of motion Static black hole metrics in D dimensions contain blackening factors which go like 1 − C r D−3 + ..., where r measures radial distance in spherical polar coordinates and C is a constant. In the large D limit the maximum order in D that can appear in the expression of Riemann tensor comes from two r derivatives acting on the blackening factor and hence R A BCD = O(D 2 ). The Ricci tensor is obtained from contraction of two of the indices of Riemann tensor and if this contraction happens along the unit sphere part of the Schwarzschild metric then an additional factor of D can arise in Ricci tensor. But we know that to zeroth order in α the black hole metric satisfies R AB = 0 and hence R AB = O(D 2 ). Similarly, since the black hole metric is also a zero of the Ricci Scalar to zeroth order in α, we have R = O(D 2 ). Hence, to linear order in α the maximum order in D coming from the Gauss-Bonnet part of the equations is O(D 4 ), as it is quadratic in internal curvatures of the metric.
So, naively it looks like that the O(α 0 ) terms in the equation (2.3) are subdominant as compared to the O(α) terms in the large D limit. This makes it implausible to have solution to the EGB equation in the large D limit which are perturbatively close to the solution of pure Einstein's equation. But we expect such solutions to exist as we consider the Lovelock theories to be corrections of Einstein-Hilbert gravity. The only plausible solution to this apparent puzzle is that α must be O(D −2 ). More specifically we will choose α to be given by where β is an O(1) quantity. We will find in the next section that this specific choice helps us to write the metric of the black hole in a manner where the only D dependence in the metric is in the fall off behaviour of the blackening factor as the black hole metric approaches the asymptotically flat spacetime. Written in terms β the EGB equation becomes

The Static Black Hole in EGB gravity and the starting ansatz
The static spherically symmetric Schwarzschild like black-hole solution for EGB gravity in Kerr-Schild form is given by where the blackening factor f (r) is given by In the above r h is the position of the horizon. Since we are interested in effects only linear in β, the relevant form of the blackening factor is The blackening factor can be written in a more covariant looking form as where K is the trace of extrinsic curvature of the surface r = r h , embedded in flat spacetime and is given by K = D−2 r h . The horizon is now denoted by the surface ψ = 1 since ψ = r r h . Let us now define two new vectors: 1) a 'velocity' vector u = −dt and a vector normal to constant ψ slices given by n = dr. The vectors are orthonormal in the sense that u.u = −1, n.n = 1 and n.u = 0 (3.5) where the dot products are taken w.r.t the flat spacetime. Let us also define a unique combination of these two vectors O = n − u. Using these variables the static black hole metric in EGB gravity can be written as In the large D limit the above metric has a fast direction along dψ as it goes to flat spacetime exponentially fast outside the horizon in this direction. This can be seen by probing regions of width of the order of O(1/D) outside the horizon using the change of variable With this change of variable in the D → ∞ limit the blackening factor becomes So, there is non-trivial warping of the spacetime only up to distances of the order of O(1/D) outside the horizon where it decays exponentially fast to the asymptotically flat spacetime.
One can now promote the quantities ψ and u to be scalar and vector functions of the spacetime coordinates with order O(1) derivatives. The spacetime is still flat when ψ − 1 ∼ O(1) and hence it continues to solve the EGB equations there. In the non-trivial region where ψ − 1 ∼ O(1/D), the EGB equations are no longer satisfied and we will need to correct the metric (3.6) to make sure it satisfies the EGB equation. This correction can be found in a perturbative series in 1/D and most of the paper will be dedicated to finding the leading order correction to the metric. We call the metric (3.6) the 'starting ansatz metric' 4 as it solves the EGB equations to the leading order approximation of the functions ψ and u being constants.
The starting ansatz metric can be thought of as being dual to a membrane propagating in flat-spacetime with the shape of the membrane being given by the equation ψ = 1. u can be thought of as the time-like velocity of points on the membrane along its world volume. Given the shape and velocity of the membrane the dual black hole metric is given by (3.6). The corrections to the metric will also be written in terms of the shape and velocity data on the membrane.

Auxilliary conditions on the data of the membrane
The information about the shape of the membrane dual to the black hole is contained entirely in the zeroes of the function ψ −1. And similarly, the data about the velocity field is contained entirely on the function u defined on the surface ψ = 1. But in order to get a local expression for the black hole metric we need to define a way to lift these functions away from the membrane surface. In writing the ansatz metric (3.6), we made a particular choice for these extensions of the ψ = 1 surface by globally defining ψ = r h r and u = −dt. This choice though convenient is not 'covariant' and hence difficult to implement in our computations. Instead we make use of the choices used in [2,6] and uplift the ψ function such that n.∇n = 0.
(4.1) 4 As a side remark we notice that the starting ansatz metric looks exactly like the staring ansatz for charged black holes in [3] if we make the replacement where 'Q ′ is the charge parameter used in the paper [3]. These similarities will also persist in some aspects of the final results.
Similarly, the u function is uplifted using the equation The covariant derivatives in the above expressions are defined w.r.t. the flat spacetime. We will call the above constraints on the ψ and u functions as the 'auxiliary conditions'. These auxiliary conditions are not essential to defining the membrane and hence the above choice is in no way unique and there are other valid choices for these auxiliary conditions (see e.g. [3,4]). The difference between various choices of auxiliary conditions will appear in the final expression for the metric corrections, but the physical content remains the same. When the EGB equations are applied to the ansatz metric (3.6), the equations are not solved beyond leading order in 1/D. The 'amount' by which they are off from being a solution will depend on the particular choice of auxiliary conditions and hence the correction that needs to be added to make them a solution also depends on the choice of the auxiliary conditions. But the final metric at a given order in 1/D should be the same once they are expressed in terms of the data on the membrane since, they solve the same set of equations with same set of boundary conditions.

The isometry ansatz and the effective equations
Taking a large D limit makes sense only when we keep some quantity fixed while taking this limit. If we are to gain insight into some aspects of finite D black hole dynamics from studying the problem in the large D limit the natural thing to do will be to work with configurations which preserve a large isometry SO(D − p − 3) as D → ∞, keeping p fixed. In these configurations all the non-trivial dynamics takes place in the finite number of p + 3 directions. A very useful method to implement this is to reduce the full D dimensional Einstein's equation into a set of effective equations in p + 3 dimensional 'effective spacetime' [2,3,6].
To do this we choose the most general metric which preserves the isometry mentioned above, namely, The equation ( where the barred quantities are w.r.t the p + 3 dimensional effective metric. The second set of equations with both legs in the i, j... directions are proportional to the metric of the d dimensional sphere and hence is effectively just one more equation given by Total number of degrees of freedom in the metric g µν are (p+3)(p+4) 2 and one degree of freedom coming from the scalar φ. This matches precisely with the total number of effective equations, so naively it looks like that we have a set of well defined deterministic set of equations.
6 Zooming in to the near horizon region and the scaled coordinates Since the non-trivial physics of black holes is confined to a region of width O(1/D) outside the horizon, it will be convenient to use a coordinate system which is obtained by zooming into the a small region of size 1/D about a point x 0 in this non-trivial region. This will be particularly useful to us as we will be performing most of our computations on computer algebra and working with these coordinates will make it much easier to keep track of orders of D 5 . The coordinates adopted to this zoomed in region is given by 6 But the metric in the y a coordinates still measures order O(1/D 2 ) distance for order O(1) coordinates distance as shown below.
We rescale our metrics in this region so that they measure order O(1) distance for order O(1) coordinates distance in the following manner.
Similarly, the derivatives in the y a coordinates become order O(1/D) And we define new scaled variable χ a given by The above rescalings make sure that in addition to working in the membrane region, we also work with metrics which measure O(1) distances and the derivatives are also O(1). This will help us to keep track of only those orders of D which appear solely due effects from the presence of black holes and the large isometry.
The effective equations (5.2) and (5.3) expressed in terms of the scaled metric G ab and χ a are given by and,

Choice of patch coordinate
Here we will specify a particular choice of scaled coordinates 'y a ' [2,3], that we will use for our computations. To start with, we note that there are three distinct directions in the geometry under consideration, namely, 1) the unit spacelike normal n µ to the ψ = 1 surface, 2) the unit time-like one-form 'velocity' u µ and 3) the distinct radial direction normal to the isometry directions i.e. dS. We will use a particular set of linear combination of these three directions as three of our basis vectors for our coordinate system namely, n, O = n − u and Z = dS − (n.dS)n. In addition we will use an arbitrary set of orthonormal coordinates Y i for the rest of the p basis vectors.
These Y i coordinate directions are also orthogonal to the hypersurface defined by the n, O and the Z directions. There is a lot of ambiguity in choosing the Y i coordinates but our final result will be written in a manner independent of the details of this choice. The scaled coordinates y a used by us is given by In these scaled coordinates ansatz metric (3.6) in the the effective p dimensional spacetime is given by Our primary task will be to add corrections to the leading order ansatz metric (3.4) of the form so that the full metric (7.3) satisfies the EGB equations to desired order in 1 D−3 . Note that h µν,β is the leading order ansatz metric given in (3.6).

Solving for the metric corrections
The functions ψ and u present in the ansatz metric have order O(1) derivatives in the global coordinates, hence these are almost constants with respect to the patch coordinates y a in the membrane region. But the metric has order D derivatives along the dψ direction. Hence, to the leading non-trivial order in 1/D, when the EGB equations act on the ansatz metric, they reduce to functions of the coordinate R only. This can easily be seen from the leading order in 1/D piece of the ansatz metric written in (7.2). The coefficients of the various terms in these functions of R are composed of Taylor expansion coefficients of ψ and u about the particular point under consideration. So, we see that the corrections that we need to add to the ansatz metric need only be functions of R. Hence, the EGB equations acting on the ansatz metric plus the corrections become ordinary differential equations (ODE) in R for the corrections with inhomogeneous sources formed with taylor expansion coefficients of the shape and velocity functions. The schematic form of these equations look like where D 0 and D 1 are ordinary differential operators in R. This simplification would not have been possible in finite D as there would not have been a parametric seperation of scales between the derivatives along dψ and other directions.

Tensor Structure of the EGB equations
For the purpose of convenience of computation we will classify the equations according to their index structure w.r.t the directions orthogonal to the u, n and dS directions. We will classify the equations as scalars vectors and tensors depending upon the number of indices in the equations along the y i directions in the effective spacetime. Since the differential operators are linear, equations of a particular tensor structure can only contain metric corrections of that tensor structure. Similarly the source to the homogeneous parts of the differential equations must also transform under the same tensor structure. This simplifies the computation process as one can solve each sector independently, as there is no mixing between different tensor sectors.
Here we will list out all independent objects that can appear in the inhomogeneous part of the ODEs at the leading non-trivial order in 1/D. We will list out the independent objects in each tensor sector after imposing the following constraints on the shape and velocity functions.
In addition we also take into account the auxiliary conditions on the shape and velocity data, namely n.∇n = 0, n.∇u = 0 Our choice are the same as those of [2,3] and are given by

tensor sector
t µν 1 = P µα P νβ K αβ − P αβ p P θφ K θφ , and t µν 2 = P µα P νβ ∇ (α u β) − P αβ p P θφ ∇ θ u φ P µν is a projector orthogonal to the three distinct directions in the effective p + 3 dimensional geometry, namely u, n and dS and K µν is the extrinsic curvature of the constant ψ slice corresponding to the point of interest in the membrane region.

Gauge choice
We have to fix the gauge for the metric corrections to take care of the unphysical diffeomorphism degrees of freedom of the EGB equations. We will use the gauge choice of [3,6] given by The most general correction to the metric at leading order in the above mentioned coordinates system can be written as Here T ij is traceless in the sense that T ij δ ij = 0. The letters denoting the metric corrections: S, V and T indicate the tensor sector to which they belong.

Boundary conditions and definition of velocity and shape of the membrane
In our coordinate system the starting ansatz metric is regular everywhere except for the singularity inside the horizon. We would like to preserve this property and hence, we require that the metric corrections that we find are also regular everywhere in spacetime away from the above mentioned singularity. For regions far outside the membrane region the spacetime is flat and hence it is always regular. Since regions inside the horizon are causally disconnected from everything outside we will not be worried about the behaviour of metric corrections far inside. Hence, we need to make sure that the solutions are regular everywhere in the membrane region. Secondly, we will only be interested in those solutions which approach the flat space time outside the membrane region as the dψ direction is the fast direction for the metric. This can be implemented by making sure that the metric and scalar field corrections vanish as the coordinate R approaches infinity.
We find that even after using the above boundary conditions we are still left with some unfixed degrees of freedom in the solution. The ambiguity is worth one scalar plus one vector degree of freedom. The source of this ambiguity is the fact that the leading order ansatz solution is invariant if we redefine the shape function by ψ → ψ + 1 D δψ and the velocity function by u → u + 1 D δu. But these redefinitions add O(1/D) piece to the ansatz metric and hence they modify the corrections to the metric at O(1/D). Hence, to unambiguously specify the corrections we must give a precise definition of the shape and the velocity field. This is similar to choice of frame in fluid-gravity correspondence.
We adopt the geometric definitions of the shape and velocity fields, proposed in [6]. We demand that the position of the event horizon coincides with the membrane surface ψ = 1 7 . Hence, at ψ = 1 the normal to the surface must be a null vector w.r.t. the full metric. If G M N is the full metric of the spacetime at a desired order in 1/D then the condition that the ψ = 1 surface is event horizon is given by Similarly, the velocity field u is the generator of the event horizon for the starting ansatz metric and we demand that this property holds for all orders in 1/D, i.e The above two conditions can easily be shown (see [6]) to result in the following conditions on the metric corrections.
Once the regularity and boundary conditions are imposed, the final solution is unambiguously expressed in terms of the shape and velocity function of the membrane 8 .

The Constraint EGB Equations and the Equations of motion of the Membrane
As we mentioned earlier, the EGB equations are invariant under diffeomorphisms. To fix this ambiguity we chose a particular gauge for the metric corrections which reduce the number of independent components in the correction metric from (p+3)(p+4) 2 to (p+2)(p+3) 2 (in addition we have one scalar correction), though the number of components of the EGB equations remain (p+3)(p+4) 2 (plus one scalar equation). So, naively, it looks like that we have a over-determined system at hand. A Similar situation 7 The surface ψ = 1 is the event horizon for the starting ansatz metric to leading order in 1/D. Hence, we expect that the horizon of the metric with the subleading corrections to be in the vicinity of width 1/D about this surface. Since the original surface was a null hypersurface, we expect that the horizon of the corrected metric which is 1/D shifted from its zeroth order position to still be a null surface. We have not exactly found the horizon surface by extrapolating null rays from future infinity, but we expect that this procedure will also lead to the same surface. We thank Benjamin Withers for discussions on this matter 8 We also need to define the auxiliary conditions on the shape and velocity field which in our case is given by (4.1) arises in two derivative Einstein gravity and there the solution to this apparent problem is well known. There the solution lies in the fact that not all components of the Einstein's equations track the 'evolution' of a metric from a set of arbitrary initial data on a hypersurface. There are exactly p + 3 components of the equations which constrain the 'initial' data on the hypersurface itself. These equations put constraints on the possible initial data which one can use. Once a consistent set of initial data has been defined, rest of the (p+2)(p+3) 2 equations dynamically evolve this 'initial' data to give the metric correction at any point in spacetime. In addition once the constraint equations are solved on any given hypersurface (normal to the direction of evolution) and the dynamical equations are solved everywhere, then the constraint equations are automatically satisfied on any hypersurface orthogonal to the direction of evolution. It was shown in [28] that a similar structure of constraint and dynamical equations exist for EGB gravity also.
Since we will be interested in studying 'evolution' of the metric components using the EGB equations along the normal to the membrane direction, the constraint equations in this case are where n µ is the normal to the ψ = 1 hypersurface and G µν is the full metric at the desired order in 1/D. E µν are the EGB equations. In the classification of EGB equations according to their tensor structure mentioned above, there are two distinct type of constraint equation: scalar and vector. We will explicitly write down the constraint equations in both these sectors below and also demonstrate following arguments similar to [6] that these equations along with the imposition of the condition that the metric corrections are regular everywhere, gives rise to a set of constraints on the shape and velocity data on the membrane. We will call these the membrane equations of motion.

Scalar constraint equations
There are three scalar constraint equations formed by taking the components of the constraint equations C µ along u, dz and O directions. We explicitly write down the homogeneous parts of the constraint equations along u and dz directions below as these will determine the constraints on the scalar data coming from the shape and velocity of the membrane. The constraint equation along O direction gives no additional constraint on the data. The relevant scalar constraint equations are given by S u const (R) and S z const (R) schematically denote the sources for these ODEs which are formed from the independent scalar data listed in (8.1.1)

Vector constraint equations
The expression for the vector constraint equations is given by S V const,i are the sources for the vector constraint equations formed using the vector data listed in (8.1.1) Notice that in both the scalar and vector constraint equations the order of the equations is one. For the rest of the dynamical equations the order of the ODE's will be 2. This is in agreement with the discussion above about the constraint equations not being dynamical in nature. To completely determine the metric using the EGB equations we need to specify the metric components and their normal derivative as the initial data on a given hypersurface. And the constraint equation being only first order imposes restrictions on this initial data.

The Membrane equations of motion
Since the constraint equations can be solved on any constant R slice, we will chose the surface R = 0 for convenience. The regularity of the metric corrections everywhere in the membrane region makes the homogeneous parts of the equations go to zero at R = 0 and we arrive at the following two constraints on the scalar data The explicit form of the constraints on the shape and the velocity data is given by We call the above two equations the scalar membrane equations of motion. Similar arguments gives rise to the following constraint on the vector data S V const,i = 0 (8.10) Written explicitely this gives We call this the vector membrane equation of motion.

The dynamical equations and the leading order metric correction
In the above section we have used the constraint EGB equations to impose the condition of regularity on the metric corrections on the R = 0 slice and we found that to be able to impose this condition, the shape and velocity data on the membrane must satisfy certain constraints given in (8.9) and (8.11). Once the constraint equations are taken care of we can now proceed to solve for the metric corrections using the dynamical equations. We will not explicitly write down the dynamical EGB equations here and just write down the final solutions after all the boundary conditions have been imposed.
where P µν = g µν + u µ u ν − n µ n ν − z µ z ν 1 − n 2 S and n S = n.dS (8.12) 9 Geometric form of the metric and the membrane equation of motion We have written down the metric corrections and the membrane equations in the effective p+3 dimensional spacetime. The computation has been done using computer algebra. Following [2, 3, 6] we do the above computation for different values of p = 2, 3. The results can be recast in terms of shape and velocity data of the membrane in full spacetime (obviously preserving the huge isometry). The procedure for this has been explained in great details in [3] and we will not elaborate on that here. This way of writing the answer in terms of the shape ane velocity data of the full D dimensional spacetime was called as the geometric form of the solution. We write down below our answer in the geometrical form and it turns out to be independent of p.

Metric correction in geometric form
The most general form of the first order corrected metric in the geometric form (with our choice of gauge) is given by The metric corrections int he effective spacetime written in terms of the full spacetime metric mentioned above is given by K M N is the extrinsic curvature of the membrane surface ψ = 1 embedded in the full flat spacetime.

Membrane equations in the geometric form
Similarly, the membrane equations of motion of (8.9) and (8.11) can be recast in a geometrical form. In this language there is a single scalar equation obtained from the first line of (8.9) given by ∇.u = 0 (9.3) and a single vector equation obtained from (8.11) given by The second of (8.9) is simply the divergence of the above vector equation ( as was the case in [3]) and is given by This last equation contains no new information.

A world volume stress tensor
By comparing the equations of motion of the membranes dual to black holes in EGB gravity (9.4) and (9.5) with the corresponding quantities of [3] we see that they map to each other under the substitution We have already noticed that the starting ansatz metric for both EGB gravity and Einstein-Maxwel theory are the same under the above substitution. Taking cue from this we propose that the stress tensor for the membrane dual to EGB black holes at leading order in 1/D can be obtained from the stress tensor for membranes dual to charged black holes derived in [27], upon using the above mentioned substitution. The proposed leading order stress tensor for EGB gravity is then given by 9 This stress tensor is the world-volume stress tensor of the membrane and hence we will only be interested in the equation of conservation of this stress tensor in the world volume of the membrane moving in flat spacetime. The conservation equation has one free index and this can either be along the velocity field u (which is the only distinct direction in the world-volume), given by 2) or it could be in a direction orthogonal to the velocity field given by The explicit expression for E is To obtain the above expression we have used one of the Gauss-Codazzi equations for the surface embedded in flat spacetime (we call this the equation of the first type) given by ∇ M K M N = ∇ N K. 9 We have neglected terms in the stress tensor whose divergence gives O(1) quantities. This is so because we will be proving that the divergence of stress tensor which naively should have been atleast O(D 2 ) is actually only O(D) if equation of motion of the membrane is satisfied. We will not be interested in finding the explicit form of the O(1) quantity that will be left out in the conservation equation.
We massage the last two terms of the equation (10.4) to turn them into more familiar objects. First of all The first equality is obtained by neglecting the term with derivative on u as it is O(1). Then we commute the derivatives and use the fact that ∇ · u = O(1) for the leading order ansatz metric. The second line is obtained by using the second kind of Gauss-Codazzi equation for a flat spacetime background given by where R ABCD is the Riemann tensor of the induced metric on world-volume. Finally we neglect a term of O(1) to obtain the last line of the above equality. Also, The first equality is obtained from the fact that the commutator of second and third derivative turns out to be zero. The second equality is obtained because ∇ · u is O(1) and third line is obtained by using both of the Gauss-Codazzi equations or flat background. Now substituting the above two expressions in (10.4) we obtain If we now use the equation obtained by taking the divergence of the membrane equation of motion (9.5) we get But we already know that which is one order lower than the naive expectation. Hence, within the order of 1/D that we are interested in, E = 0 provided the membrane equation to leading order is satisfied.
Let us now shift our attention to E M which is given by In the above expression we have already neglected terms which are O(1) and have also used first of the Gauss-Codazzi equations. We massage the last term using similar arguments as above to get If we substitute the above expression (10.9) we find that it is of the order O(1) since the naive order O(D) terms is proportional to the equation of motion of the membrane. Hence, we conclude that to leading order in 1/D the proposed world volume stress tensor of the membrane is conserved provided the membrane equations are satisfied to leading order in 1/D or vice-versa.

Quasi-Normal Modes from Linearised Membrane Equation
In [2,6] the spectrum of linearised fluctuations of the membrane dual to black holes in Einstein-Hilbert gravity was found to leading and first subleading order in 1/D. It was demonstrated that the spectrum exactly reproduced the light quasi-normal mode spectrum of [1] to relevant order in 1/D. Similar analysis was undertaken in [3] and a prediction for light QNMs of static black holes in Einstein-Maxwel systems was given. Here we undertake a similar analysis and give a prediction for the light QNMs of static, spehrically symmetric black holes in EGB gravity. A static, spherical membrane in spehrical polar coordinates is given by r = 1 and u = −dt Let us parametrise small fluctuations about this static configurations by r = 1 + ǫδr(t, θ), and u = −dt + ǫ (δu t (t, θ)dt + δu a (t, θ)dθ a ) (11.1) We will now evaluate the membrane equations to first order in ǫ. The various terms in the equation evaluate to u.∇u a = ǫ∂ t δu a (11.2) A covariant derivative with a subscript or superscript (s) denote the covariant derivative w.r.t. the D − 2 dimensional unit sphere metric. Using the above expressions the linearised membrane equations become Taking a divergence of the first equation w.r.t the unit sphere metric and then using the second equation we get The shape fluctuation can be decomposed into a basis of spherical harmonics w.r.t the D − 2 dimensional unit sphere in the following manner This when substituted in the linearised membrane equations above gives a quadratic equation for the frequencies w s for the scalar mode perturbations at a given angular momentum number which when solved gives The velocity fluctuation δu a can be decomposed into a pure vector (divergence-less) part and a part that is derivative of a scalar function w.r.t. the unit sphere, Substituting this in the equation (11.3) and upon using the solution of δr from above gives the frequency for the vector mode perturbations to be The scalar modes described above are in agreement with the linearised fluctuations about a static spherically symmetric black hole in flat spacetime for EGB gravity obtained in equation (4.18) of [21] from an independent method of obtaining large D effective equations. We could not find a corresponding solution for the vector modes. There are no analytic expressions of light quasinormal modes in EGB gravity for the particular scaling of the GB parameter we worked with from the gravity side and hence, we could not check our solutions against anything else. But taking cue from the Einstein-Hilbert case we propose (11.6) and (11.8) to be the predictions for the light QNM spectrum of Static black holes in EGB gravity. For the scalar modes for l = 0 we get w s = 0 or 2i(1 − β), but only the w s = 0 is compatible with the second equation of (11.3). This zero mode corresponds to the time independent scaling of the spherical membrane which is trivially a solution of membrane equation. For l = 1 both the roots vanish and we seperately look at this mode in equation (11.4) which becomes −D∂ t δr + D∂ r δr + (1 + β)D∂ 2 t δr = 0 =⇒ δr l=1 = At + B (11.9) These two independent zero modes correspond to a time independent translation ( when A = 0) and an infinitesimal boost of the black hole when (B = 0). Similarly, there is a zero mode for the vector perturbations at l = 1 which simply correspond to imparting an infinitesimal rotation to the black hole . This line of discussion remains unaffected in the presence of a GB parameter as can be seen from similar observations in [3,6].

Summary and Discussions
In this paper we have worked out the large D membrane paradigm for Einstein-Gauss Bonnet Gravity to the leading non-trivial order in 1/D. We have worked out the leading corrections to the black hole metric of (3.6) which solves the EGB equations to first subleading order. In the process we have also worked out the leading order constraints that needs to be put on the shape and velocity data of the membrane to get regular solutions to the EGB equations. We call these the membrane equations (8.9), (8.11). Our expressions for metric correction and the membrane equations are valid to linear order in the Gauss-Bonnet parameter only.
We have also written down an expression for a world-volume stress tensor for the membrane based on similarity of the current system with the large D black holes of Einstein-Maxwel system. The equation of conservation of this stress tensor (10.1) correctly reproduces the equation of motion of the membrane. We have also worked out the spectra of light quasi-normal modes of static spherically symmetric black holes in EGB gravity in the leading large D limit.
In this paper we lay the groundwork towards the grander goal of understanding the status of second law of black hole thermodynamics for a general consistent classical theory of gravity with an action containing more than two derivatives. The membrane equations in the two derivative Einstein gravity at the first subleading order in 1/D was shown to be equivalent to the second law of black hole thermodynamics [6,27]. We expect that the correct interpretation of the membrane equations for EGB gravity at first subleading order in !/D may give us important insight on the status of the second law of black hole thermodynamics for these theories (atleast in the large D limit). The results of the current paper can be used as the base to compute the corrections to the membrane equations and the metric corrections to the first subleading order in 1/D.
One can also work out the membrane paradigm in some other higher derivative terms in the Lovelock theories of gravity, so that some general understanding of the membrane equation of motion and stress tensor can be obtained for these theories. This will also be important to in gaining some insight on the mechanism by which second law of black hole thermodynamics may be satisfied in general consistent theories of classical gravity.
In addition adapting the algorithm spelled out in [27] for higer derivative theories of gravity, one can figure out the space-time stress tensor localised on the membrane which is compatible with the linearised gravitational field at infinity from the black holes of EGB gravity. This stress tensor should also have the additional property that its conservation equation should correctly reproduce the equation of motions of the membrane. The relevant part of this stress tensor when evaluated on the world volume of the membrane should be equivalent to the expression of world-volume stress tensor given in (10.1). Also, a proposal for consistent membrane stress tensor at finite D which reduces to the correct stress tensor in large D as in [13] for EGB gravity will also be interesting.