The large $D$ membrane paradigm for general four-derivative theory of gravity with a cosmological constant

We find the membrane equations which describe the leading order in $1/D$ dynamics of black holes in the $D\rightarrow\infty$ limit for the most general four-derivative theory of gravity in the presence of a cosmological constant. We work up to linear order in the parameter determining the strength of the four-derivative corrections to the gravity action and hence there are no ghost modes in the theory. We find that the effective membrane equations we obtain are the covariant version of the membrane equations in absence of the cosmological constant. We also find the world-volume stress tensor for the membrane whose conservation gives the membrane equations. We apply the membrane equations to predict the light quasi-normal mode spectrum of black holes and black branes in the theory of gravity under consideration.

1 Introduction The equations of general theory of relativity governing the dynamics of spacetime are a set of coupled partial differential equations on the metric of spacetime. Since, there are no algorithmic way to solve these equations, exact analytic solutions to these equations are extremely difficult to find and there are very few such solutions known. Of them a class of interesting solutions are those of static or stationary black holes. Even though these solutions are exact, it is next to impossible to study any non-trivial dynamical processes involving these black holes analytically. For the vacuum Einstein equations there are no free parameter which can be tuned to suitable limits to get more analytic control over the problem, so that we can gain more insight into the structure of these equations and solutions.
In the last few years starting with the work of [1][2][3][4][5], the limit of taking the spacetime dimensionality to infinity has been found to be particularly useful to get more analytic control over the dynamics of black holes. In particular, it was shown in [6,7] that in the D → ∞ limit there is a decoupled sector of finite number of light quasi-normal modes of black holes, which completely determine the late time small amplitude dynamics. These modes live in a region of width of O(1/D) about the horizon of the black hole. In addition there also is an infinite tower of heavy quasi-normal modes which die off exponentially fast at late times. Since there is a parametric separation between the energies of the light and heavy modes, there must be an autonomous system of non-linear equations determining the complete late dynamics of the black holes and which in the small amplitude limit gives rise to these light quasi-normal modes.
The effective theory of these light modes were discovered independently in two different incarnations. In one of the approaches [8,9] the effective equations are obtained in terms of effective mass and momentum variable associated with the black hole. In the other approach [10] the theory of light modes is described by a system of equations determining the dynamics of a co-dimension one membrane moving in asymptotic spacetime of the black hole. The linearised spectrum about static black holes obtained from both these approaches matched with the light quasinormal modes mentioned above.
Since then both these approaches have been extended in diverse directions. Both these approaches have been generalised to black holes with charge and/or in presence of cosmological constant [11][12][13][14][15]. In particular the effective equations obtained in the mass and momentum variable were particularly useful in understanding the endpoint of IR instability of black strings/branes [16] in flat spacetime (see also [17]). The relation between the hydrodynamic limit and the large D limit have also been explicated in the mass-momentum approach in [18] and in the membrane picture in [19,20]. The large D formalism has been applied to obtain many interesting gravity solutions. It has been applied to find new interesting stationary black-bar solutions in [21,22] and also to find analytic bulk solutions in AdS dual to Bjorken flows [23].
The equivalence between the effective equations obtained in these two different approaches have not been established in complete generality but in many different situations the equivalence between the two approaches have been established: e.g the equivalence of the effective equations for stationary and static black holes have been established in [24] and the equivalence of the effective equations about black branes have been established in [25]. For the purposes of this paper, we will focus our attention on the membrane picture of large D black hole dynamics.
A particularly interesting observation was made in [26] where it was shown using the membrane picture that the entropy current (determining the entropy of the dual black hole) has a positive definite divergence. This information was encoded in the second order in 1/D membrane equations obtained in [27]. Hence, the membrane equations somehow know about the second law of black hole thermodynamics without it being used as an input anywhere in the formalism explicitly. The immediate natural question that one can ask is if the large D membrane picture can tell something about the status of the second law of black hole thermodynamics for the most general classical theory of gravity. The status of the second law of black hole thermodynamics is not yet understood in full generality. A candidate entropy for most general higher curvature theory of gravity was found in [28] which follows second law for linearised dynamics about a stationary solution. Another candidate entropy which obeys second law for higher curvature gravity in situations which preserve spherical symmetry was obtained in [29]. The membrane approach can be used to fill some of this gap by finding a candidate entropy preserving the second law for the most general theory of gravity in the large D limit.
Keeping this goal in mind, it is imperative that the effective equations for black hole dynamics in large D limit are derived in many different scenarios involving higher derivative theories of gravity to gain a lot of intuition about the structure of the membrane equations in these theories. The large D formalism to study black hole dynamics for the simplest consistent higher derivative theory of gravity namely Einstein-Gauss-Bonnet (EGB) gravity was first applied to study the quasi-normal modes of static black holes from the gravity picture in [30]. The authors of this paper and other collaborators worked out the effective equations for black holes and black branes dynamics in the large D limit in terms of the mass and momentum variables in [31][32][33]. The large D effective equations in terms of the membrane variables for the EGB gravity (to linear order in the Gauss-Bonnet (GB) parameter) was worked out in [34].
In this paper we take forward the investigation in this direction to find the membrane effective equations to leading order in 1/D in most general four-derivative theory of gravity in the presence of a negative cosmological constant. We will work perturbatively order-by-order in the parameter determining the strength of the fourderivative part of the gravity action and we will work with only those branches of solutions which have a nice limit to solutions of Einstein-Hilbert action. This makes sure that there are no ghost modes in the theory under consideration. We will provide all our results up to linear oder in this perturbation series.
We will try to make this paper as self-contained as possible. Since the rest of the paper will contain many technical details involving the computation, we provide a brief summary of the results obtained here.

Summary of results
We will work with the most general theory of four-derivative gravity whose action is given by where l is related to the cosmological constant Λ by .
κ is the parameter specifying the strength of the four-derivative term in the gravitational action. We will show that for our purposes a scaled parameterκ related to κ byκ = κ(D − 4)(D − 3), will be of particular interest. We show that to linear order inκ the membrane equation determining the leading order in large D dynamics of black holes is given by where all the derivatives and dot products in the above equation are taken w.r.t the AdS spacetime to which the black hole asymptotes to. The AdS radius of the asymptotic spacetime is not completely determined in terms of the cosmological constant 'l', but rather is modified in presence of the higher derivative terms in the action as The above scaling is true to leading order in large D. The membrane is defined by its shape and a unit normalised velocity vector field u M in its world-volume. K M N is the extrinsic curvature of the membrane and K is the trace of the extrinsic curvature. P M A is the projector orthogonal to the velocity field in the world-volume of the membrane.
We also find a stress tensor in the membrane world-volume whose conservation gives the membrane equations mentioned above. This stress tensor is given by where, As an application of the membrane equations, we compute the spectrum of linearised fluctuations about the static spherical membrane and the static flat membrane. These spectra give predictions for the light quasi-normal modes about static black hole and black branes in the gravity theory under consideration. The scalar and vector frequencies of the spectrum of quasi-normal mode about the static black hole are given by and where j represents the angular momentum number of spherical harmonics of the background SO(D − 2) isometry group. The frequency of scalar quasi-normal mode about the black brane are obtained for two different scalings in D for the momenta.
The quasi-normal mode for momenta of the order of O(1) are given by The spectra when the momentum of the modes are O( √ D) are given by In all the results for quasi-normal mode about black branes mentioned above k denotes the amplitude of the momenta along the planar directions.

Most general four-derivative theory of gravity
The most general action for a theory of four-derivative gravity in the presence of a negative cosmological constant can be written as Here, κ is a constant with the dimension of the square of length and a 1 · · · a 5 are the arbitrary dimensionless coefficients determining the relative strength of the five possible terms with a total of four derivatives acting on the metric. We take the point of view that the four-derivative terms are corrections over the Einstein-Hilbert action and we will be only looking for solutions that can be perturbatively connected to the solutions of the Einstein-Hilbert action. In this paper, we will provide all results up to linear order in κ only. We also take the length scale corresponding to the cosmological constant given by l −1/2 to be an O(1) number in any spacetime dimension.
Since the terms appearing with coefficients a 4 , a 5 can be recast as the divergence of a vector, they can be integrated by parts and converted to boundary terms. So, these terms do not contribute to the bulk gravity equation of motion and we will drop these terms from the action from now on. The action that we will work with in this paper is given by 1 (2.2) 1 The dual hydrodynamics of this action has been studied in [35,36].
This action reduces to the Einstein-Gauss-Bonnet action for the particular values of the coefficients a 1 = 1, a 2 = −4 and a 3 = 1. The gravity equation of motion corresponding to this action is given by It is easy to see from the above equation of motion that for the combination of the coefficients for which the action reduces to the Einstein-Gauss-Bonnet action, there are no terms in the equation with more than two derivatives acting on the metric 2 . Hence, for generic choice of coefficients there can be ghost modes (modes which are ill-behaved in the sense that they can travel faster than the speed of light) in the theory. Since, we will be working perturbatively in the parameterκ about solutions of two-derivative gravity and so we will not encounter these ghost modes.

The leading order ansatz from static black hole
The aim of this work is to find dynamical solutions with horizons of the gravity equations (2.3) in a series in 1/D. In the limit of D → ∞ each new term in the series adds to the accuracy of the solution in the D → ∞ limit 3 . To start this procedure we need to find a solution which to leading order in large D solves for the gravity equation and also has a "horizon". We will follow the algorithm mentioned in [10,11,13,14,27,34] of using the metric of static black holes in Kerr-Schild coordinates to find the starting ansatz solution.
The metric of the static spherically symmetric black hole which asymptotes to AdS spacetime and solves the gravity equation (2.3) (to linear order in κ) is given by where f (r) = 1 + lr 2 +κÃl 2 r 2 − g(r), We can do a coordinate transformation under which we replace the length element dt with dt = dv to arrive at the metric ds 2 = ds 2 AdS + g(r) 1 + lr 2 +kÃl 2 r 2 1 + lr 2 +kÃl 2 r 2 dt + dr and ds 2 AdS = −dt 2 (1 + lr 2 +kÃl 2 r 2 ) + is the AdS spacetime to which the black hole asymptotes to. The AdS radius of this spacetime is given by 1 √l where,l = l(1 +κÃl). The above metric is the Kerr-Schild form of the metric of the static black hole in asymptotic AdS spacetime. The trace of the extrinsic curvature of the surface r = r h embedded in the asymptotic AdS spacetime defined above is given by Using this, the Kerr-Schild metric (3.4) can be written in a more covariant form as where, and as already mentioned, ds 2 AdS is the metric of the AdS spacetime whose length scale is 1 √l .
We elevate ψ and u to the status of analytic functions of the spacetime coordinates with O(1) derivatives while keeping the following properties intact dψ.dψ so that n · n = 1, u · u = −1, u · n = 0. Of primary interest to us will be the near-horizon "membrane region" where, ψ − 1 ∼ R D . To leading order in large D then ψ ∼ e −R and following the arguments in [13,14] it can be shown that the metric in this region solves the equation of motion provided ∇ · u = O(1/D).
From discussions in subsequent sections it will be clear that the above constraint to leading order in large D will imply that the velocity field does not point in the direction of the radius of the isometry direction. Once this is satisfied the metric in the membrane region to leading order maps to the metric of static black hole to leading order and hence obviously solves for the gravity equations. The region of spacetime where ψ ≪ 1 is in the interior of the "horizon" and we will be able to maintain this at subsequent orders in the 1/D perturbation series. Since this region is causally disconnected from the perspective of the observer outside the horizon it plays no role in the physics in that region and hence we will not worry about this here.
Thus the metric (3.8) can be taken to be the leading order ansatz metric for our 1/D perturbation theory as it solves the gravity equations to leading order in the regions of interest.

The dual membrane, its data and choice of auxiliary conditions
The ansatz metric (3.8) can be thought of as being expressed in terms of data of a membrane propagating in the asymptotic AdS spacetime of the black hole. The data defining the membrane comprise of its shape defined by the surface ψ = 1 (which specifies the horizon of the spacetime to leading order in large D) and a unitnormalised time-like velocity field u µ defined in the world-volume of the membrane. The fact that the velocity data is defined only in the world-volume can be inferred from the orthonormality of the velocity and the normal vector to the surface ψ = 1 in (3.8). The large D limit will allow us to express further corrections to the metric in a systematic 1/D expansion, in terms of higher derivatives of the shape and velocity data of the membrane. This is the reason we call it the "large D membrane paradigm". As explained above, we will be interested in finding the metric corrections in the 'membrane region' of width of the order O(1/D) about the surface ψ = 1. To get local expressions of the metric in terms of the data of the membrane located at ψ = 1 we need to define the extension of the shape and velocity data off the membrane surface. The extension of the shape data of the surface is related to the choice of the family of surfaces of which the surface ψ = 1 is a particular member e.g. this unique surface can be part of two (of many) distinct family of surfaces say ψ 2 − ψ = 0 or ψ = 1 + f (ψ), where, f (ψ = 1) = 0. The choice of the extension of the velocity vector of the membrane surface may be thought of as the worldvolume velocity fields defined on each member of the family of surfaces that one chooses. There are many different choices of doing this and the expression of the metric corrections will depend on these choices, but it is clear from the procedure that the physics content of these different metric corrections in terms of the actual data on the membrane surface remains unchanged. Following earlier works, we will keep calling this choice the choice of auxiliary conditions of the membrane data and we will make the following choice used in [10,27,34]  The above auxiliary conditions can be physically thought of as the lifting of the shape data by demanding that the normal vector at each point follows its geodesic and the velocity vector is parallel transported along this geodesic.

The effective metric, dilaton field and the effective gravity equations
In D spacetime dimension, the number of independent components of the metric and the number of gravity equations are D(D+1)

2
. Clearly, in the large D limit, both the number of equations and variables are infinite. We solve the problem of correcting the metric at subleading order in 1/D using computers and so we will use the methods used in [10,11,27,34] to reduce the infinite number of gravity equations to a finite number effective equations using the fact that the spacetime under considerations preserves a large SO(D−p−3) isometry . A convenient way to parametrise spacetime that preserves SO(D − p − 3) isometry is to write the metric as Here g µν is an effective metric in the p + 3 directions where, non-trivial dynamics takes place and we also have an effective dilaton field φ determining the non-trivial warping of the isometry direction of the full spacetime. Both the effective metric and the dilaton fields depend only on the p + 3 dimensional coordinate x µ . Due to symmetry properties of the spacetime the gravity equations of motion split in two distinct classes: the first one having components only along the p + 3 x µ directions and the second class having components only along the d dimensional spherical directions. Also due to symmetry considerations all the i, j equations are proportional to the unit sphere metric Ω ij . The dimension of the isometry group d appears as a parameter in the equations.
The component of the gravity equations in the p + 3 directions are given by is the Laplacian w.r.t the full metric of spacetime. The quantities in this equation evaluate in terms of the effective metric of the reduced spacetime and the dilaton field to the following expressions Whereas the equations in spherical directions are given by where,

Effective equations in near horizon region coordinates
The ansatz metric has a fast direction dψ along which derivatives acting on the metric can be of order O(D). We find it convenient to go to a "scaled" coordinate system based around any arbitrary point in the membrane region so that the maximum order of derivatives along the fast direction acting on the metric becomes O(1). We use the procedure already mentioned in [10,11,34] for this purpose. We briefly review the procedure in this section. First of all we zoom into the membrane region about any arbitrary point x µ 0 using the following coordinate transformation, 5 Due to this coordinate change, the effective p + 3 dimensional metric and the derivative of the dilaton field in the 'y a ' coordinates transform as Here we have introduced a new field corresponding to the derivative of the dilaton field namely χ a . Doing the above-mentioned scalings of the coordinates and fields we now have a system where metrics measure O(1) distances for unit coordinate distances and also the maximum order of derivatives on the fields is O(1). We write down the effective gravity equations obtained in the last section in these scaled system 5 We will use 1 D−3 , rather than 1 D as the perturbation parameter throughout this paper as this is the natural variable that appears in the blackening factor of black holes in arbitrary dimensions.
now. The equations in the p + 3 directions of the reduced spacetime are given by where, The equations in the isometry directions which are proportional to the unit sphere metric Ω ij are given by (6.6)

A choice of coordinate system
Having specified the set of equations that we will work with in the last section we will now explain the choice of the coordinate system that we use for our computation in this section. To do that it will be useful to understand the reformulation of the asymptotic metric of the AdS black brane geometry which solves for the gravity equations (2.3) in terms of the effective spacetime. The metric written above has an ISO(D −2) isometry. This metric can be recast in a manner where only a SO(D − p − 3) isometry is manifested in the following way Here we have introduced an ad hoc spacial direction along dS along which the manifested ISO(D − 2) isometry is broken into SO(D − p − 3) × ISO(p). This way of writing the metric helps us to identify the effective metric and dilaton field for the spacetime to be given by 3) The shape and velocity functions in terms of which the ansatz metric is defined are functions of only the p + 3 directions of the reduced spacetime 6 . In addition the velocity field has components only along these p + 3 directions. Of the p + 3 directions 3 of the directions are distinct namely those along: 1) The direction of the velocity field 'u', 2) the direction of the normal to the membrane 'n ∝ dψ' and 3) the direction along which the isometry is broken from i.e. 'dS'. Before explicitly stating the choice of the coordinates let us evaluate the divergence of the velocity vector field in the coordinates (7.2). It is given by We need to make sure that the divergence of the vector vanishes at leading order in large D so that the leading order ansatz vanishes. We choose u r ∼ O(1/D) and u S ∼ O(1/D) for this purpose. This makes sure that the divergence of vector field is O(D 0 ) and we denote this as ∇ · u = U, (7.5) where U is an O(D 0 ) quantity. We now follow the convention of [10,11,34] to choose 3 of the coordinates to be along n, O = n − u, and dZ = dS − (n · dS)dS. We choose O as one of the coordinate directions so that we have a null (ingoing) coordinate which makes it easier to impose regularity conditions at the horizon. The dZ direction is chosen so that it is orthogonal to both n and u which is made sure by the choice of u r and u S mentioned above. We assign the following names to the coordinates in the membrane region about an arbitrary point with coordinates In terms of these coordinates the ansatz metric to leading order in large D is given by The ansatz metric is not a solution to the gravity equations at subleading orders in 1/D. We need to add corrections to the ansatz metric and the dilaton field in a 1/D expansion to get a solution to the gravity equations up to the desired order. The corrected metric and dilaton field can be parametrised as µν,κ , is the effective metric part and φ (0) is the dilaton field part of the ansatz metric.
The gravity equations and their effective counterpart equations are a set of nonlinear coupled partial differential equations and there are no well defined algorithmic procedure to extract analytic solutions out of them. The ansatz metric that we work with has a fast direction along R direction in our coordinate system and hence, when the gravity equations act on the ansatz metric, the piece that is left unsolved at the subleading order is a function of the coordinate R only with coefficients of the various R-dependent parts in the function being ultra-local functions of shape and velocity data. To be more precise, the derivatives transverse to the R direction acting on the metric are always suppressed by O(1/D) w.r.t leading order in R derivatives and hence, around a given point in the membrane region they only contain information about local data of the membrane like the extrinsic curvature of the membrane or the derivative of the velocity field about that point. These local data then appear as coefficient of the non-trivial R-dependent terms in the gravity equation at subleading order. For this reason the correction metric at this order can be thought of as function only of R coordinate and the gravity equation acting on them produce homogeneous pieces in R derivatives acting on the correction metric. Hence, the gravity equations that we need to solve become ordinary differential equations (ODEs) in R coordinate acting on the correction pieces with sources being determined by local membrane data like extrinsic curvatures etc. In particular the differential equations that we arrive at are particularly simple to solve.

Tensor structures of the equations
As explained above, in the large D limit the gravity equations that we need to solve are ODEs involving the metric corrections with sources determined by the membrane data. The corrections to the effective spacetime metric can be decomposed according to their tensor structures w.r.t the p equivalent directions orthogonal to u, n and dZ directions. The ODEs can also be classified according to this tensor structure in the sense that in a given tensor sector the homogeneous parts of the ODEs involve metric corrections of that particular sector only. Pushing the argument a bit further it is easy to see that the inhomogeneous parts of the ODEs comprising the shape and velocity functions also appear in the same tensor structure as the homogeneous parts. The different tensor sectors are completely decoupled. This is of particular use to us for solving the problem on a computer as each tensor structure can be solved separately.
Similar strategy has also been followed in [10,11,27,34] where the independent shape and velocity data which can appear as sources for the leading order ODEs have been classified in the tensor sectors in the following manner.
The above set of independent sources are obtained after imposing the orthonormality constraints n · n = 1, u · u = −1 and u · n = 0, and the auxiliary conditions n · ∇n = 0 and n · ∇u = 0.

The choice of gauge
The metric corrections are arbitrary up to coordinate redefinitions (diffeomorphism).
To get a well-defined corrected metric at subleading order, we need to make a choice of gauge so that the redundant diffeomorphism degrees of freedom are removed. We make the following gauge choice (found particularly useful in earlier works on the large D membrane paradigm) This gauge choice when implemented in our choice of patch coordinates, the most general form of metric correction in the different tensor sectors can be written as (9.1)

Boundary conditions and regularity of solutions
The equations of gravity that we work with (2.3) and its counterpart in effective spacetime, have terms up to two derivatives acting on the Riemann tensor (hence, four derivatives acting on the metric). Naively it looks like that the ODEs that we need to solve in the large D limit are of order four. But we work only to linear order in the variableκ and we work towards finding only those corrections to the metric that smoothly match to solutions of two derivative Einstein gravity in theκ → 0 limit. This is encoded in the way we have parametrised the corrections to the ansatz metric order-by-order inκ in the previous section. All of these together make sure that the ODEs we need to solve are of degree two for both h (1) ab and h (1) ab,κ . Since, the equations have order two, the solutions to them come with two undetermined constants. We explain below how to fix these constants.
We impose the regularity conditions on the metric corrections that they are analytic functions everywhere in the membrane region. This we do since we work with a coordinate system where to start with there are no singularities anywhere in the membrane region. We also need to impose boundary condition that the metric corrections vanish outside the membrane region, i.e. as R → ∞. This makes sure that the metric of the spacetime approaches asymptotic AdS spacetime exponentially fast. These boundary and regularity conditions fix most of the undetermined constants for us.
The boundary and regularity conditions though are not enough to take care of all of the undetermined constants. We find that we are left with two undetermined constant in the scalar sector and a one vector number worth of constant in the vector sector of the metric corrections. The origin of this can be traced back to ambiguities in the definition of the shape and velocity functions in the ansatz metric as noticed in [13-15, 27, 34]. The ansatz metric that we use in this paper is of the form 7 Here, the notation F (ψ) means that F is a function of spacetime coordinates via its dependence on the shape function ψ. The function F (ψ) must satisfy the additional constraints Under a generic small change in the shape and velocity function parametrised by ψ → ψ + δψ and δu → u + δu, the change in the ansatz metric is given by The form of the function F 8 is such that the leading order ansatz metric remains unaffected by changes in the shape and velocity function of the form where, δũ and δψ are O(1). The redefinition of the velocity field is not entirely free but is constrained by the requirement that the redefinition keeps the velocity field unit normalised. This requires that δũ · u = 0.
Due to the above redefinition, the direction of the normal vector remains along dψ up to errors of order O(1/D 2 ). In absence of any redefinition, the surface ψ = 1 is a null surface w.r.t the ansatz metric. After the redefinition, the norm of the normal to the redefined surface n ∝ dψ is given by This fixes one of the undetermined constants in the scalar sector.
Another interesting observation is that the velocity vector field raised w.r.t the asymptotic spacetime is the generator of the event horizon i.e. the surface ψ = 1 to leading order in 1/D. Under the redefinition of the velocity field that we mention, it can be shown that the generator of event horizon on the surface ψ = 1 is given by Imposing the boundary conditions that we mentioned here, in addition to the regularity and fall off conditions, the first order metric is completely determined in terms of the shape and velocity data of the dual membrane along with the auxiliary conditions.

The constraint Einstein equation and the membrane equations
The second order differential equations that we need to solve acquire the properties of the Einstein-Hilbert equations. This is easy to see because the ODEs at linear order inκ have homogeneous parts containing the terms in the metric corrections given by h ab,κ . Since, these metric corrections already come with aκ multiplied with them their homogeneous parts are completely determined by the action of the twoderivative part of the Einstein equations on them. Hence, not all of the equations are dynamical in nature 9 . Some of the equations are constraints on the metric data at any given slice of the evolution direction. These equations are given by where, E µν = 0 are the gravity equations, G να is the inverse of the metric of the spacetime and n α is the normal to the ψ = constant slices. All other components of the Einstein equation are dynamical in nature in the sense that given the data on a given ψ = constant slice, these equations determine the data at the subsequent slices. The interesting property of Einstein equations is that once we demand that the dynamical equations are solved everywhere and the constraint equations are solved on any ψ = constant slice then the constraint equations are solved everywhere in the spacetime. In the language of ODEs that we solve the dynamical equations are of order two and the constraint equations are of order one.
We separate the constraint equations according to their tensor structures into scalar and vector equations. Scalar constraints are obtained by picking the component of (11.1) along u, dz and O directions. The scalar constraint equations of interest to us will be the components along u and dZ. The structure of the equation along u direction is given by where, H 1 (R) and H 2 (R) are functions of R with membrane data dependent sources. Since, the constraint equations carry same information irrespective of the choice of slice, we evaluate the above equation at R = 0. Using the boundary conditions on S V z and S V z,κ and assuming that the metric corrections will evaluate to regular functions everywhere, it is clear that the homogeneous parts vanish at R = 0. Hence, the above equations are satisfied at R = 0 only if the sources of the homogeneous parts vanish at R = 0. This is not a content less statement as has been observed in [13,14,27,34]. In fact this condition puts the following constraint on the membrane data U = 0. (11.3) 9 Our ODEs do not determine the time evolution of the system but instead can be thought of determining the "dynamics" along the dψ direction.

The other scalar constraint equation along dZ direction is given by
The homogeneous parts of the above equation evaluate to zero at R = 0 and following similar arguments as before, the above equation gives the following constraint on the membrane data

The vector component of the constraint equation evaluates to
The homogeneous part of the above equation again vanishes at R = 0 and hence, puts the following constraint on the vector data of the membrane (11.7) The constraint on the membrane data obtained are in fact dynamical equations on the membrane data and hence we call them the "membrane equations". Our procedure to correct the ansatz metric to arrive at solutions to the gravity equations at subleading order in 1/D involves "patching" together many local solutions to form a global solution. The membrane equations constrain one-derivative data of the membrane and hence can be thought of as the leading order obstruction to patching together any arbitrary data to form a black hole solution. If and only if the membrane satisfies the membrane equations written above, one can construct a first order corrected solution from the membrane data.

Solution to the dynamical equations
Once the membrane equations are satisfied, the dynamical equations can be solved to arrive at regular corrections to the ansatz metric at first order in 1/D by demanding regularity and imposing the boundary conditions on the metric mentioned earlier. The metric correction in the different tensor sectors obtained are given by

Geometrized metric correction
All the metric corrections we have written down heavily use the presence of the dS direction. The dS direction is specific to the choice of what the value of p is. It has been beautifully explained in [11] that a solution which preserves a SO(D − p − 3) isometry also preserves a SO(D − q − 3) isometry for q ≫ p. This is possible only if the solutions obtained by above secretly possess a symmetry under redistribution of finite number of coordinates between isometry and non-isometry directions. The mechanism by which the solutions possess this symmetry has been explained in [11]. The idea is that the solution obtained can be written in terms of geometric quantities of a membrane propagating in full spacetime (not just the effective p + 3 dimensional spacetime) so that the dynamics preserves SO(D − p − 3) isometry for any arbitrary p. This mechanism was called "geometrisation" of the solutions. The full spacetime has only two distinct directions: a) along the velocity vector field and b) along the normal to the membrane. The metric corrections that need to be added to the ansatz metric can be decomposed into tensor structures with respect to the directions orthogonal to these two distinct directions. Together with our gauge choice, i.e H M N O N = 0, the most general metric correction in the tensor decomposition mentioned above can be written as where P AB is the projector orthogonal to n and u vectors and We need to put an extra factor of 1/D in front of the H T r to make sure that the corrections do not destroy the property that the leading order ansatz metric solves for the gravity equations at leading order in 1/D.
The map between the metric correction written in this manner and the metric correction written in the effective spacetime formalism is explicitly given in Appendix (C). This can be used to write the expressions of the metric correction in geometric form as given here

Geometrized membrane equation
The equations determining the "membrane equations" can also be written in a geometric way. In the geometric way of writing, the "vector" membrane equations can be written in terms of the geometric vector data along the direction orthogonal to the u and n vectors as where P M N = δ M N + u M u N is the projector orthogonal to the velocity vector in the membrane world-volume. In the above equation K M N is the extrinsic curvature of the membrane surface and K is the trace of the extrinsic curvature of the membrane. The covariant derivative and the dot products in the above equation are taken w.r.t asymptotic AdS spacetime of the black hole with effective radius L 2 AdS = The interpretation of the constraint on scalar membrane data coming from the constraint gravity equation projected along the velocity vector is already manifested in (11.3) and is given by Giving geometric interpretation to the other scalar membrane equation requires a bit more work. For this we need to calculate the divergence of the geometric vector membrane equation written above to leading order in 1/D. This is given by It can be shown easily that to leading order in large D where, R AB is the Ricci tensor of the induced metric on the world-volume of the membrane. Let the intrinsic curvature of the background spacetime be given by R αβγδ , R αβ , R. The relation between the intrinsic curvature of the induced metric on the membrane and the extrinsic curvature of the membrane are given by the Gauss-Codacci relations Using (14.4) and (14.5) in the divergence of the vector membrane equation one gets (14.6) It can easily be checked that the above equation when written in terms of the data of the membrane written in effective spacetime language matches with the second scalar constraint on the membrane data coming from the constraint gravity equation evaluated along dZ.

Comparison to Charged Membranes with cosmological constant
The leading order membrane equations dual to charged black hole dynamics in presence of a cosmological constant were derived in [15] and are given by It is particularly interesting to observe that when in the first and last equation written above we use the map Q → κa 3 K D we arrive at the geometric form of the scalar and vector membrane equations derived by us. Also under this map the equation determining the charge variable reduce to which to the relevant order is the divergence of vector equation written above. The map between the charge and shape variable mentioned above is precisely which converts the leading order in 1/D ansatz metric written in [15] to the ansatz metric written by us in (3.8). We find this map 10 very intriguing and we think it may be useful to understand its origin better.
Another interesting thing we notice in this context is that naively the charge equation written above seems to have a piece which is not covariantisation of the flat space answer obtained in [11], namely the piece u·R·u K . But it is worth noticing that for both the charge equation and the divergence of the vector equation mentioned above this factor always comes int he combination which can also be written as ∇ M (u·∇u M )

K
. Written this way this piece is still a covariantisation of the answer written in flat spacetime. We think that this is a more natural way to express the membrane equations from the context of the expression of the charge current derived in [26], the relevant part of which is given by The conservation of the above charge current gives the charge membrane equation and in the conservation equation the term that appears is ∇ µ (u · ∇)u µ . So, in summary we can conclude that all the membrane equations in asymptotic AdS spacetime that we know of at leading order in large D are covariantisation of the answer in asymptotic flat spacetime.

Stress tensor
The membrane dual to black holes in leading order in large D has a world-volume stress tensor in the sense that its conservation gives rise to the membrane equations.
This stress-tensor is given by where, The equation for conservation of the stress tensor along u m is given by Using the Gauss-Codacci relations given in the last section one can arrive at the following useful identities Combining the above equations we get (15.5) The right hand side is precisely the divergence of the vector membrane equation of motion and is O(1) and hence, we arrive at the scalar membrane equation Similarly, the conservation of the stress tensor orthogonal to the velocity vector in the membrane world-volume gives which upon using the Gauss-Codacci relation gives precisely the vector membrane equation of motion. Hence, the membrane dynamics can be stated in a very compact form as the conservation of the stated stress energy tensor in the world-volume of the membrane.

Light quasi-normal modes
It has been shown in [10,11,13,14,27] that the large D membrane equations can predict the spectrum of light quasi-normal modes. The membrane equations can be used to predict the spectrum of light modes for situations where an explicit answer is not known from gravity side. In [31] the spectrum of light quasi-normal modes for static black holes in Gauss-Bonnet gravity was predicted from the effective equations derived in terms of effective mass and momentum variables. The prediction for these modes from the effective membrane equations in [34] matched with the results derived in [31] in the common regime of validity. In this section we compute the spectrum of linearised fluctuations about a static spherical membrane and a static planar membrane to predict the spectrum of light quasi-normal modes for static black holes and static black branes in the theory of gravity under consideration. We start with the computation for spherical membrane in the next subsection.

Light quasi-normal modes of black hole
To study the spectrum of linearised fluctuations about static-round membranes we consider a coordinate system adapted to the spherical symmetry of the problem. The metric of the AdS spacetime to which the black holes asymptote to and in which the membrane propagates, in these coordinates is given by ds 2 = −dt 2 1 + lr 2 +κÃl 2 r 2 + dr 2 1 + lr 2 +κÃl 2 r 2 + r 2 dΩ 2 D−2 . (16.1) We will express all our answers in terms of the variablel given by where, the effective AdS length scale of the AdS spacetime is given by L 2 AdS . Let the static spherical surface be given by r = 1.
We will be concerned with small fluctuations about this configuration which can be parametrised as r = 1 + ǫδr(t, θ a ). (16.3) The unit normalised space-like normal to this surface is given by The extrinsic curvature is given by The superscript s indicates the fact that this is the extrinsic curvature computed as a spacetime quantity. In the membrane equation of motion the pull back of the extrinsic curvature on to the membrane world-volume enters and the relation between the two is given by where x A are the spacetime coordinates and x µ are the world-volume coordinates. The above equation is to be evaluated at the position of the membrane given by The non-zero components of the pulled back extrinsic curvature are given by 11 (16.6) The trace of the extrinsic curvature K can as well be evaluated in the bulk at the membrane surface and its expression is given by K = g tt ∇ t n t + r rr ∇ r n r + g ab ∇ a n b | r=1+ǫδr , = ǫ∂ 2 t δr (1 +l 3/2 ) +l 1 +l where, ∇ (s) denotes gradient w.r.t the unit sphere coordinates. The induced metric on the membrane world-volume is given by Since, the background configuration is static the background velocity field points in the dt direction. The unit normalised background velocity field is given by u = − 1 +ldt, 11 The details of the computation are provided in the appendix (B.1). and the fluctuations about it can be parametrised as u = − 1 +ldt + ǫδu t (t, θ a )dt + ǫδu a (t, θ a )dθ a . (16.8) Imposing the condition that the above velocity field is unit normalised w.r.t the induced metric up to linear order in the amplitude fluctuation parameter ǫ gives δu t (t, θ a ) = −l δr 1 +l .
The scalar membrane equation ∇ · u = 0 then evaluates to We can write the vector membrane equation as We evaluate the non-zero components of the projector orthogonal to the u vector in the world-volume of the membrane given by to be Collecting the expressions written above, the component of the linearised vector membrane equation along the angular direction is given by The vector equation above can be decomposed into a part which is the gradient of a scalar function and pure divergenceless vector part as given below Once can take a divergence of the above equation to get which in the configuration that we are working with evaluates to ǫ 1 +l ∇ (s) a ∇ 2 δu a D + 2 − (1 +κa 3 (1 +l))∂ t ∇ (s) a δu a + (1 +l)∇ (s) a δu a −l∇ (s) a δu a We can decompose the angular component of the velocity vector field also into a gradient of a scalar part and a divergenceless vector part as The background has an SO(D − 2) × R (1) isometry and all the functions are characterised by their decomposition into spherical harmonics and the frequency in each spherical harmonic sector e.g.
where the spherical harmonics obey Substituting the above decomposition of δr into the equation for δr gives the spectrum of the scalar fluctuations to be Substituting the decomposition of δv a into vector spherical harmonics sectors given by δv a = j,m d j,m Y j,m a e −iw j v t into the above equation gives the spectrum of the vector quasi-normal mode frequency as These results match with the scalar and vector quasi-normal mode frequencies of AdS black hole found in [7] and [13] in the limitκ → 0. For Gauss-Bonnet combination i.e. a 1 = 1, a 2 = −4 and a 3 = 1, the quasi-normal mode frequencies agree with [30,31] in the appropriate limit.
It is easy to see that the mode j = 0 is a zero mode for the scalar mode. This can be seen both from the expression of the frequencies of the spectrum given above and also from the fact that the decoupled equation for δr is trivial for j = 0 and the scalar membrane equation for j = 0 gives w = 0. This mode corresponds to a uniform scaling of the size of the black hole and is obviously a solution of the gravity equations. There are no zero modes corresponding to j = 1 scalar mode. The presence of this zero mode would have implied that translations and boost of the centre of mass of the black hole are solutions of the gravity equation. This is not the case here as AdS spacetime has a non-trivial warping.
Similarly, the vector sector has a zero mode at j = 1 and corresponds to uniform rotation of the black hole about an axis, which also is obviously a solution of the gravity equations.

Quasi-normal mode of a black brane
In this subsection, we compute the frequency of linearised fluctuations about a membrane dual to a AdS black brane in the four-derivative gravity theory under consideration. The asymptotic spacetime of the black brane in which the membrane propagates is given by where,l = l(1 +κÃl).
We consider a small amplitude fluctuation around the static membrane solutions dual to black brane r = 1 parametrised by r = 1 + ǫδr(t, x a ). (16.20) The unit normal to this membrane surface is given by Non-trivial components of the pullback of the spacetime extrinsic curvature on the membrane r = 1 + ǫδr are given by The trace of the extrinsic curvature on the membrane surface takes the form The induced metric on the membrane is We parametrise small fluctuation around static velocity configuration u = − l r 2 dt by u = − l dt + ǫδu t (t, a)dt + ǫδu a (t, a)dx a . Explicit expressions for the terms appearing in the vector membrane equation are this scaling limit the terms in the above equation that contribute at leading and subleading order in large D limit are We keep the first subleading parts to get a non-zero answer for the frequency. Substituting δr = δr 0 e −iwst e ikax a in (16.34), we arrive at the spectrum of the scalar fluctuations δr = δr 1 0 e −iw 1 s t e ikax a + δr 2 0 e −iw 2 s t e ikax a , (16.36) We again consider the most general form of δu a (16.37) and substitute it in the vector membrane equation to get (16.38) and the frequency of the vector fluctuations to be In the absence of higher derivative correction, these results agree with [13].

Outlook and future directions
The main result of this paper is the computation of the membrane equations dual to the dynamics of black holes in the most general four-derivative theory of gravity to leading order in a 1/D expansion. Like in [34] we have worked to linear order in the parameter specifying the relative strength of the four-derivative terms in the action w.r.t the two-derivative Einstein-Hilbert part of the action. The AdS radius of the asymptotic spacetime of the solutions is modified in presence of the higher derivative terms as L 2 AdS = 1 l(1+κa 1 l) . The effective equation of the membrane when expressed in terms of quantities expressed in this AdS spacetime is the covariantisation of the corresponding membrane equation dual to black holes dynamics in asymptotic flat spacetime for the four-derivative theory of gravity. Together with the observation made for Einstein-Hilbert gravity in [13] this can be seen as a general property of membranes dual to large D dynamics of black holes in absence of charge: To leading order in large D membranes dual to black holes in asymptotic AdS spacetime have equations which are covariant versions of the corresponding equations in flat spacetime. We do not see any obvious reason for this to be true and we think that this needs to be understood better.
Another interesting thing that we observed is the map between membranes dual to four-derivative theory of gravity up to linear order inκ and those dual to charged black holes [11,15], continue to persist even in the presence of the cosmological constant. Under the map the charged membrane equations map exactly to the membrane equations derived here. In fact the non-trivial appearance of the curvature of the background spacetime in the charge case appears non-trivially in the divergence of the vector membrane equation here via the Gauss-Codacci equations. Since, there is no known independent theory of membranes which can be dual to black hole dynamics in general, the above two observations can be thought of as data points to get better intuition about the "transport coefficients" of the membranes. It will be interesting to extend the analysis to find the dual membranes in more general backgrounds to get a better understanding of the structure of the membrane equations and stress tensor. A first step in this direction in the mass and momentum approach has been taken in [37]. The authors of this paper derived the effective equations in AdS in presence of arbitrary boundary deformation with particular large D scaling. Due to the choice of this scaling they observe the presence of the effect of curvature of the asymptotic metric even at leading order. It will be interesting to study similar phenomena in the membrane picture as this may help us in formulating the effect of background curvature on the membrane.
It will also be interesting to extend the analysis to subleading order both inκ and in 1/D. The analysis at subleading order inκ will help us to understand if the phenomenon of covariantisation of the membrane equation still persists or not. The analysis at subleading order in 1/D will be important from the point of view of collecting more data about the second order membrane equations which played crucial role in the derivation of the entropy current for two-derivative theory of gravity.
We again do a coordinate change of the form dV = lr 2 (1 +κÃl)dt + In the D → ∞, at r ≫ r h , the function g(r) vanishes, which implies that in the large D limit, blackening factor is non-trivial only around a small thickness around the horizon. It is also worth to notice that the background geometry(A.7) is modified in the presence of the higher derivative parameter. The above metric (A.6) can be written in a more covariant form given by Components of ∇ M n N are ∇ t n t = −ǫ ∂ 2 t δr l r 2 −l 3/2 r 2 , ∇ t n r = ǫ ∂ t δr l r 2 , ∇ t n a = −ǫ ∂ t ∂ a δr l r 2 , ∇ r n t = 2ǫ ∂ t δr l r 2 , ∇ r n r = 0, ∇ r n a = 2ǫ ∂ a δr l r 2 , ∇ a n t = −ǫ ∂ a ∂ t δr l r 2 , ∇ a n r = ǫ ∂ a δr l r 2 , ∇ a n b = −ǫ ∂ a ∂ b δr l r 2 + l l r 2 δ ab . (B.7) Components of the spacetime extrinsic curvature are K s tt = −ǫ ∂ 2 t δr l r 2 −l 3/2 r 2 , K s tr = ǫ ∂ t δr l r 2 , K s ta = −ǫ ∂ t ∂ a δr l r 2 , K s rr = 0, K s ra = ǫ ∂ a δr l r 2 and K s ab = −ǫ ∂ a ∂ b δr l r 2 + l l r 2 δ ab . ∇ t u t = 0, ∇ t u a = ǫ ∂ t δu a + l ∂ a δr , ∇ a u t = 0, C Map between metric corrections in geometric and effective spacetime This appendix can be thought of as a generalisations of the explicit map worked out in Appendix B of [11]. The projector in full spacetime orthogonal to the normal and velocity vector can also be expressed as where, p A B is the projector orthogonal to n, u and the z direction in the effective spacetime and Ω AB is the unit sphere metric along the ismometry directions and it has non-zero component only along the θ i directions specifying the isometry directions.
We will now find the map between the geometric form of the metric correction of (13.2) and the correction written in the effective spacetime of (9.1). Let us start with the metric correction H