A leading-order comparison between fluid-gravity and membrane-gravity dualities

In this note, we have compared two different perturbation techniques that are used to generate dynamical black-brane solutions to Einstein equation in presence of negative cosmological constant. One is the `derivative expansion', where the gravity solutions are in one-to-one correspondence with the solutions of relativistic Navier-Stokes equation. The second is the expansion in terms of inverse power of space-time dimensions and here the gravity solutions are dual to a co-dimension one dynamical membrane, embedded in AdS space and coupled to a velocity field. We have shown that in large number of space-time dimensions, there exists an overlap regime between these two perturbation techniques and we matched the two gravity solutions along with their dual systems upto the first non-trivial order in the expansion parameter on both sides. In the process, we established a one-to-one map between dynamical black-brane geometry and the AdS space, which exists even when the number of dimensions is finite.


Introduction:
It is very hard to solve Einstein's equations -the key equation governing the dynamics of space-time, even at classical level. Only few exact solutions are known, mostly being static or stationary. To handle non-trivial dynamics, we have to take recourse of perturbation. Two such important perturbation schemes, which can handle dynamical fluctuations around static solutions even at non-linear level, are 'derivative expansion' [1][2][3][4][5][6] and expansion in inverse powers of dimension [7][8][9][10][11][12][13][14][15][16]. The first one generates 'black-hole' type solutions (i.e., space-time with singularity shielded behind the horizon) that are in one to one correspondence with the solutions of relativistic Navier Stokes equation 1 whereas the second one generates similar 'black hole type' solutions, but dual to the dynamics of a codimension one membrane embedded in the asymptotic geometry 2 . It is natural to ask whether it is possible to apply both the perturbation techniques simultaneously in any regime(s) of the parameter-space of the solutions, and if so, how the two solutions compare in those regimes. In this note we would like to answer these two questions. In a nutshell our final result is only what is expected.
• It is possible to apply both the perturbation techniques simultaneously. Further, in the regime where both D is large and derivatives are small in an appropriate sense, we could treat 1 D and ∂ µ (with respect to some length scale) as two independent small parameters, with no constraint on their ratio.
• In other words, if the metric dual to hydrodynamics is further expanded in inverse powers of dimension, it matches with the metric dual to membranedynamics, again expanded in terms of derivatives.
However, this matching is not at all manifest. We could see it only after some appropriate gauge or coordinate transformation of one solution to the other. The whole subtlety of our computation lies in finding the appropriate coordinate transformation.
The 'large-D' expansion technique, as described in [7], generates the dynamical black brane geometry in a 'split form' where the full metric could always be written as a sum of pure AdS metric and something else. In other words, the black-brane space-time, constructed through 'large-D' approximation would always admit a very particular point-wise map to pure AdS geometry. On the other hand, the space-time dual to fluid dynamics does not require any such map for its perturbative construction and apparently there is no guarantee that the particular map used in 'large-D' technique, would also exist for the dynamical blackbrane geometries, constructed in 'derivative expansion'. In this note, we have shown that the 'hydrodynamic metric' 3 indeed could be 'split' as required through an explicit computation upto first order in derivative expansion. This map could be constructed in any number of dimension and is independent of the 'large -D' approximation. After determining this map, we have matched these two different gravity solutions upto the first subleading order on both sides. We believe it would get more non-trivial at next order but we leave that for future. One interesting outcome of this exercise, is the matching of the dual theories of both sides. It essentially reduces to a rewriting of hydrodynamics in large number of dimensions, in terms of the dynamics of the membrane. After implementing the correct gauge transformation, we finally get a field redefinition of the fluid variables (i.e., fluid velocity and the temperature) in terms of membrane velocity and its shape 4 . We hope such a rewriting would lead to some new ways to view fluid and membrane dynamics and more ambitiously to a new duality between fluid and membrane dynamics in large number of dimensions, where gravity has no role to play (See [19], [16] for a similar discussion on such field redefinition and rewriting of fluid equations though in [19] the authors have taken the large D limit in a little different way than ours).
The organization of this note is as follows. In section - (2) we first discussed about the overlap regime of these two perturbation schemes. Next in section - (3) we discussed the map between the bulk of the 'blackhole' space-time and the pure AdS, mentioned above and described an algorithm to construct the map, whenever it exists. In section -(4) we compared the two metric and the two sets of dual equations (controlling the fluid-dynamics and the membrane dynamics respectively) within the overlap regime, upto the first subleading order on both sides. This section contains the main calculation of the paper. We worked out the map between these two sets of dual variables, leading to a map between large D relativistic hydrodynamics and the membrane dynamics. Finally in section - (6) we concluded and discussed the future directions.

The overlap regime
In this section we shall discuss whether we could apply both 'derivative expansion' and 1 D expansion simultaneously. We shall first define the perturbation parameters for each of these two techniques in a precise way and also fix the range of their validity. We shall see that these two parameters are completely independent of each other and therefore their ratio could be tuned to any value, large or small.
Next we shall compare the forms of the two metrics, determined using these two techniques, assuming the ratio (between the two perturbation parameters) to have any arbitrary value.

Perturbation parameter in 'derivative expansion'
Here we shall very briefly describe the method of 'derivative expansion'. See [17] for a more elaborate discussion.
The technique of 'derivative expansion' could be applied to construct a certain class of solutions to Einstein's equations in presence of negative cosmological constant in arbitrary dimension D.
The key gravity equation: λ is the inverse of AdS radius. From now on, we shall choose units such that λ is set to one. These gravity solutions are of black hole' type, meaning they would necessarily have a singularity shielded by some horizon [3]. They are in one-to-one correspondence with the solutions of relativistic Navier-Stokes equation in (D − 1) dimensional flat space-time (without any restriction on the value of D). In fact, we could use the hydrodynamic variables themselves to label the different gravity solutions, constructed using this technique of 'derivative expansion'. The labeling hydrodynamic variables are At the moment r H is just some arbitrary length scale, which would eventually be related to the horizon scale of the dual black-brane metric. {x µ }, µ = {0, 1, · · · , D −2} are the coordinates on the flat space-time whose metric is simply given by the Minkowski metric, η µν = Diag{−1, 1, 1, 1 · · · }.
'Derivative expansion' enters right into the definition of the hydrodynamic limit. The velocity and the temperature of a fluid are functions of space-time but the functional dependence must be slow with respect to the length scale r H (x). For a generic fluid flow at a generic point, it implies the following.
Choose an arbitrary point x µ 0 ; scale the coordinates (or set the units) such that in the transformed coordinate r H (x 0 ) = 1. Now the technique of derivative expansion would be applicable provided in this scaled coordinate system In other words, the number of ∂ α derivatives in a given term determines how suppressed the term is 5 . In terms of original x µ coordinates, each derivative ∂ µ cor- 5 The conditions as described in (2.2) are for a generic situation. For a particular fluid profile, it could happen that at a given point in space-time some nth order term is comparable to or even smaller than some (n + 1)th order term. One might have to rearrange the fluid expansion around such anomalous points if they exist, but they do not imply a 'breakdown' of hydrodynamic approximation. As long as all derivatives in appropriate dimensionless coordinates are suppressed compared to one, 'derivative expansion' could be applied. responds to r H∂µ . Therefore if we work in x µ (which, unlikex µ , are not defined around any given point) coordinates, the parameter that controls the perturbation is schematically ∼ r −1 H ∂ µ . 6 . The starting point of this perturbation is a boosted black-brane in asymptotically AdS space. The metric has the following form (in coordinates denoted as {r, x µ }, µ = {0, 1, · · · , D − 2}. Units are chosen so that dimensionful constant, λ, appearing in equation (2.1) is set to one) 7 .
is an exact solution to equation (2.1) provided u µ and r H are constants. Now the algorithm for 'derivative expansion' runs as follows. Suppose, u µ and r H are not constants but are functions of {x µ } . Equation (2.3) will no longer be a solution. If we evaluate the gravity equation E AB on (2.3), the RHS will certainly be proportional to the derivatives of u µ and r H . But u µ and r H being the hydrodynamic variables, their derivatives are 'small' at every point in the sense described in (2.2). Therefore a 'small' correction in the leading ansatz could solve the equation. The r dependence of these 'small corrections ' could be determined exactly while the {x µ } dependence would be treated in perturbation in terms of the labeling data u µ (x) and r H (x) and their derivatives. u µ (x) and r H (x) themselves would be constrained to satisfy the hydrodynamic equation, order by order in derivative expansion. While dealing with the full set of gravity equations (2.1), these equations on the hydrodynamic variables or the labeling data would emerge as the 'constraint equations' of the theory of classical gravity.

Perturbation parameter in 1
D expansion This is a perturbation technique, which is applicable only in a large number of spacetime dimension (denoted as D), as a series expansion in powers of 1 D . Clearly 1 D is the perturbation parameter (a dimensionless number to begin with) here, which 6 For a conformal fluid in finite dimension, there is only one length scale, set by the local temperature which also sets the scale of derivative expansion. But if we take D → ∞, T (x) and r H ∼ T (x) D are two parametrically separated scales and it becomes important to know which one among these two scales controls the derivative expansion. In the condition (2.2) we have chosen r H to be the relevant scale and set it to order O(1). Indeed the results in [5] seem to indicate that terms of different derivative orders in hydrodynamic stress tensor, dual to gravity are weighted by factors of r H ∼ T (x) D , and not T alone. Note that here the temperature of the fluid would scale as D, which is different from the D scaling of the temperature, imposed in [19]. 7 Note that the scaling of λ with D is upto us. At finite D it is of no relevance, but it matters while taking the large D limit. Here λ would be fixed to one as we would take D to ∞. must satisfy Unlike the derivative expansion, the 1 D expansion does not necessarily need the presence of cosmological constant, but we could also apply it if the cosmological constant is present provided we keep λ, the AdS radius (see equation (2.1) in subsection -2.1) fixed as we take D large. Note that the choice λ = 1, as we have done in previous subsection, is consistent with this 'D-scaling'.
The starting point here is the following metric.
whereḠ AB , ψ and O A are defined as follows.
1.Ḡ AB is a smooth metric of pure AdS geometry which we shall refer to as 'background'.
We could choose any coordinate as along as the metric is smooth and all components of the Riemann curvature tensors are of order O(1) or smaller in terms of large D -order counting.
2. ψ −D is a harmonic function with respect to the metricḠ AB .
3. O A is a null geodesic in the background satisfying O A n BḠ AB = 1 where n A is the unit normal on the constant ψ hypersurfaces (viewed as hypersurfaces embedded in the background).
The metric (2.4) would solve the Einstein's equations (2.1) at leading order (which turns out to be of order O(D 2 )) provided the divergence of the O(1) vector field, U A ≡ n A − O A with respect to the background metric is also of order O(1).
Naively equation (2.5) does not seem to constrain the vector field U A since each of its components along with their derivatives in every direction are of order O(1) (this is what we mean by an 'order O(1) vector field'). However, it is indeed a constraint within the validity-regime of 1 D expansion. We could apply large D techniques provided for a generic O(1) vector field V A ∂ A , its divergence is of order O(D) 8 . One easy way to ensure such scaling would be to assume that the dynamics is confined within a finite number of dimensions and the rest of the geometry is protected by some large symmetry [7]. From now on, we shall assume such symmetry to be present in all the dynamics we discuss, including the dual hydrodynamics, labeling the different geometries constructed in 'derivative expansion'. For example, we shall assume that the divergence of the fluid velocity u µ , which we shall denote by Θ(≡ ∂ µ u µ ), is always of order O(D), whereas the velocity vector itself is of order O(1). Now we shall briefly describe some general features of this leading geometry in (2.9). See [7] for a detailed discussion. Firstly note that with the above conditions, the hypersurface ψ = 1 becomes null and we could identify this surface with the event horizon of the full space-time. Also, if one is finitely away from the ψ = 1 hypersurface, the factor ψ −D vanishes for large D and the metric reduces to its asymptotic formḠ AB . Next consider the region of thickness of the order of O 1 D around ψ = 1 hypersurface. This is the region 9 , where 1 D expansion would lead to a nontrivial correction to the leading geometry. To see why, let us do the following coordinate transformation.
where {X A 0 } is an arbitrary point on the ψ = 1 hypersurface. In these new coordinates Now, ifx A is not as large as D, it is possible to expand ψ −D , O A andḠ AB around a coordinate system {z, y µ } for the background.
Here clearly the first term is of order O(1). The second term could potentially be of order O(D) since large number of indices are summed over. Still to precisely cancel against the last term, which certainly is of order O(D) as long as Vz z is not very small, it requires some fine tuning. Equation (2.5) says that U A ∂ A is such a fine-tuned vector field. 9 Following [7] , we shall refer to this region as 'membrane region' Note that from the second condition (see the discussion below equation (2.4)) it follows that Extrinsic curvature of (ψ = 1) surface = K| ψ=1 = D N A N BḠ AB + O(1) Substituting equation (2.8) in equation (2.7) we find Clearly at the very leading order, the metric will have non-trivial variation only along the direction of N A -the normal to the ψ = 1 hypersurface at point X A 0 . Variations along all other directions are suppressed by factors of 1 D . This is very similar to the metric in equation (2.3) where at leading order the non-trivial variation is only along a single direction -r. Therefore, within this 'membrane region', 1 D expansion would almost reduce to derivative expansion along directions other than N A provided the metric (2.9) solves equation(2.1) at very leading order. The conditions, listed below equation (2.4) along with equation (2.5) ensure that this is the case. Once the leading solution is found, the same algorithm, described in the previous subsection, would work and we could find the subleading corrections handling the variations of N A and O A along the constant ψ hypersurface. All such variations would be suppressed as long as none of the components of N A , O A and their derivatives (in the unscaled X A coordinates) are as large as D. In other words, we should be able choose a coordinate system, along the horizon (or the hypersurface ψ = 1) such that It is enough to impose this inequality only on the ψ = 1 hypersurface; the conditions listed below equation (2.4) will ensure that they are true on all constant ψ surfaces.
These conditions also specify the defining data (analogue of fluid-velocity and temperature in case of 'derivative expansion') for the class of metrics, generated by 1 D expansion. Here, the gravity solutions are expressed in terms of the auxiliary function ψ and the one-form O A dX A . These two auxiliary fields satisfy the second and the third conditions, listed below equation (2.4). However, the above mentioned conditions, being differential equations, could not fix the fields completely unless some boundary conditions are specified along any fixed surface. The most natural choice for this hypersurface is the surface given by ψ = 1, which, by construction, is the horizon of the full space-time geometry. Different metric solutions are classified by the shape of this surface and the components of O A projected along the surface. Just as in 'derivative expansion', we could solve for the metric correction only if these defining data (the projected O A field and the shape of the surface, encoded in its extrinsic curvature) satisfy the constraint equation, which we shall refer to as the 'membrane equation'.

Comparison between two perturbation schemes
In subsection-(2.2), we have seen that within the membrane region, O 1 D expansion is almost like 'derivative expansion' as described in subsection-(2.1). Still it is also clear that they are not quite the same. The leading ansatz itself looks quite different for the two schemes, and there is no question of overlap if these two techniques compute perturbations around two entirely different geometries. So, to find an 'overlap regime', the first step would be to see where in the parameter-space and in what sense, equation (2.3) and (2.7) describe the same leading geometry. Note that though the leading geometries look different algebraically, they both have similar geometric properties -namely the existence of a curvature singularity. In metric (2.3) it is located at r = 0 and the metric (2.7) is singular at ψ = 0. Also the singularity is shielded by some event-horizon 10 . To see the similarities more explicitly, let us first choose a coordinate system X A ≡ {ρ, X µ }, such that the background metric-Ḡ AB in equation (2.8) takes the form In this coordinate system, the following metric is an exact solution of equation (2.1) This is simply the Schwarsczchild black-brane solution, written in Kerr-Schild form. Now let us note the following features of this metric [7].
• The function ρ is harmonic with respect to the background upto So far, the way both the techniques of 'large-D expansion' and 'derivative expansion' are developed, the existence of a horizon is a must. It would be interesting to know whether we could depart from this condition and still apply either of these two techniques to construct 'horizon free' or non-singular smooth geometries. • The one form dρ ρ − ρ dt is null and satisfies the geodesic equation. Further, contraction of this one-form with the unit normal to constant ρ hypersurfaces is one. Hence this one form could be identified with the null one form O A dX A Hence it follows that the metric in (2.12), which is an exact solution of (2.1), could be cast in the form of our leading ansatz upto corrections subleading in 1 D expansion. We could also expand the metric in equation (2.12) around a given point on the horizon ρ = r H , the same way we have done (see equation (2.9)) in the previous subsection with the following set of identifications.
The very leading term in this expansion, once written in terms of N A and O A would have exactly the same form as that of the metric in equation (2.7). The main difference between our leading ansatz, equation (2.4) and equation (2.12) is that in the later N A and O A satisfy equation (2.13) everywhere along the horizon, in the same {ρ, y µ } coordinates. For our leading ansatz (2.4) also, it is true that we could always choose a local {ρ, t} coordinates by reversing the equations in (2.13). But for a generic ψ and O A , this could not be done globally and this is the reason why our leading ansatz is not an exact solution of (2.1). However, the deviation from the exact solution would clearly be proportional to the derivatives of N A and O A and therefore subleading. So finally we conclude that locally around a point on the horizon, the leading ansatz for 1 D expansion looks like a Schwarzschild black-brane written in a Kerr-Schild form with the local ρ and t coordinates, respectively oriented along the direction of the normal N A and the direction O A projected along the membrane ψ = 1. Now let us come to the leading ansatz for the metric in derivative expansion. As it is explained in detail in [2], the leading ansatz in derivative expansion, equation (2.3) , reduces to Schwarzschild black-brane in Eddington-Finkelstein coordinates if we choose r H = constant and u µ = {1, 0, 0, · · · }. Also locally at any point {x µ 0 }, we could always choose a coordinate system such that u µ (x 0 ) = {1, 0, 0, · · · }, or in other words by appropriate choice of coordinates locally the metric described in (2.3) could always be made to look like a Schwarzschild black-brane, though in a different gauge than in equation (2.4). Clearly the starting point of these different expansions are 'locally 'same and it is possible to have an overlap regime.
But the difference lies in the concept of 'locality' and also in the space of defining data. In case of 'large-D' expansion, the classifying data of the metric is specified on the horizon whereas for 'derivative expansion' it is defined on the boundary of AdS space. The range of validity for 'large-D' expansion is given in equation (2.10). If we replace ∂ A ψ −D | horizon by (−DN A ) the condition (2.10) reduces to the existence of coordinate system such that which looks very similar to the validity regime for 'derivative expansion' , as already mentioned in subsection (2.1) If we could somehow map each point on the boundary to a point on the horizon (viewed as a hypersurface embedded in the background), the same {x µ } coordinates could be used as coordinates along the horizon. In that case, whenever r H is of order O(1) in terms of 'large-D' order counting, the inequality (2.15) would imply equation (2.14). In other words, as D → ∞, all solutions of 'derivative expansion' could be legitimately expanded further in 1 D , though the reverse may not be true. Now we know that ∂ A and ∂ µ are simply related (without any extra factor of D) for the case of exact Schwarzschild black-brane solutions. This is just the well-known coordinate transformation one should use to go from Kerr-Schild to Eddington-Finkelstein form of the black-brane metric. This transformation also gives the required map from the horizon to boundary coordinates. Once perturbations are introduced on both sides, we expect the relation between these two sets of coordinate systems would get corrected, but in a controlled and perturbative manner, thus maintaining the above argument for the existence of overlap.
So in summary, there does exist a region of overlap between these two perturbative techniques. In this note, our goal is to match them in the regime of overlap. As it is clear from the above discussion, the key step involves determining the map between ∂ A and ∂ µ , which we are going to elaborate in the next section.

Transforming to 'large-D' gauge
From the discussion of section - (2) it follows that if the space-time dimension D is very large, we could always apply ' 1 D expansion' whenever 'derivative expansion' is applicable. Therefore a metric, corrected in derivative expansion in arbitrary dimension, when further expanded in 1 D , should reproduce the metric generated independently using the method of ' 1 D expansion'. More precisely if we take the metric of equation (4.1) from [5] and expand it in 1 D , it should match with the metric given in equation (8.1) of [7] after appropriate transformation. In this section our goal is to understand what these 'appropriate transformations' are.
Let us explain it in little more detail. As we have mentioned before, both of these two perturbative techniques generate black brane geometries, in terms of a set of 'dynamical data' , confined to a co dimension one hypersurface. In the first case it is the boundary of the Asymptotic AdS space and in the second case it is the event horizon viewed as a hypersurface embedded in pure AdS. So both the techniques require a map from the full spacetime geometry to a co dimension-one membrane. The details of this map are quite clear for the case of 'derivative expansion'. The data-set that distinguishes between different dynamical geometries, here is the profile of a relativistic conformal fluid (its velocity and temperature). In other words, given a unit normalized velocity field and temperature, defined on a (D − 1) dimensional flat space-time and satisfying the relativistic Navier-Stokes equation, we should be able to uniquely construct a D dimensional space-time with a dynamical event horizon such that its metric is a solution to (2.1). The (D − 1) dimensional space is identified with the conformal boundary of this D dimensional black-brane geometry, which we shall refer to as bulk. This construction [5] uses a very specific coordinate system, that encodes how a point in the bulk could be associated with a point in the boundary. In [20], the authors have also explained how to reverse the construction of [2], [5]. They have given an algorithm to read off the dual fluid variables starting from any black-brane geometry that admits derivative expansion, but written in arbitray coordinates. This explicitly proves the claim of one-to-one correspondence between the dynamical black-brane geometry, admitting derivative expansion and the fluid profile, satisfying relativistic Navier-Stokes equation. This algorithm has been heavily used to cast the rotating black-holes in the 'hydrodynamic form' [5].
Similarly according to [7], there exists a one-to-one correspondence between dynamical black-brane geometries in 1 D expansion and a codimension-one 'membrane dynamics' in pure AdS space, though [7] shows the correspondence in only one direction. It starts from a valid membrane data and integrate it outward towards infinity to construct the corresponding black-brane geometry. But to explicitly show this correspondence, we also need to know the reverse. In other words, we should know how to associate a point on the membrane to a point on the bulk and how to read off the membrane data, starting from a dynamical black-brane geometry that admits an expansion in 1 D , but written in some arbitrary coordinates. In the next subsection we shall formulate an algorithm to determine this 'membranebulk map', analogous to the discussion of [20] in the context of transforming the rotating black holes to the hydrodynamic form. 11 The 'large-D expansion' technique, as developed in [7], would always generate the dynamical black-brane metric G AB in a 'split' form. This 'split' is specified in terms of an auxiliary function ψ and an auxiliary vector field O A ∂ A . In terms of equation,

Bulk-Membrane map
whereḠ AB is the background and G (rest) AB is such that there exists a null geodesic vector field O A ∂ A in the background, satisfying The normalization of this null geodesic vector is determined in terms of the function ψ, defined as follows.
1. ψ −D is a harmonic function with respect to the metricḠ AB .
2. ψ = 1 hypersurface, when viewed as an embedded surface in full space-time, becomes the dynamical event horizon. This is how the boundary condition on ψ is specified.
After fixing ψ, the normalization of O A is fixed through the following condition.
where n A is the unit normal on the constant ψ hypersurfaces (viewed as hypersurfaces embedded in the background). The equations (3.1) and (3.2) together specify a map between two entirely different geometries, with metricḠ AB and G AB respectively, both satisfying equation (2.1). So if we want to recast an arbitrary dynamical black-brane metric, which admits 1 D expansion, in the form as described in (3.1), the first step would be to figure out this map or the 'split' of the space-time between 'background' and the 'rest', so that the equation (3.2) is obeyed. Now from the discussion of the previous subsection, we see that this 'map' is crucially dependent on the vector field O A ∂ A and the function ψ. But both of them are defined using the 'background' geometry and we immediately face a problem, since given an arbitrary black-brane metric, it is the 'background' that we are after. For example, given a black-brane metric we could always determine the location of the event horizon, but we would never know its embedding in the background, unless we know the 'split' and therefore we would not be able to construct the ψ function, by exploiting the harmonicity condition on ψ −D . If we do not know ψ we would not be able to orient or normalize O A , as required.
So we must have some equivalent formulation of this 'split' just in terms of the full space-time metric. The following observation allows us to do it. We could show that if G AB admits a split betweenḠ AB and G (rest) which is a null geodesic with respect toḠ AB , is also a null geodesic with respect to G AB .

Proof:
We know that where ∇ denotes the covariant derivative with respect toḠ AB and κ is the proportionality factor. We would like to show that SupposeΓ A BC denotes the Christoffel symbol corresponding to∇ A and Γ A BC denotes the Christoffel symbol corresponding to ∇ A . These two would be related as follows [7].Γ Here all raising and lowering of indices have been done usingḠ AB . Note that What we want to show simply follows from equation (3.4) So we could determine O A by solving the null geodesic equation with respect to the full space-time metric G AB . But to determine it fully, we also need to know κ, fixed by the normalization of O A . As mentioned before, the normalization used previously in the application of 'large-D' technique is not suitable for our purpose, since it requires the knowledge of the 'background' beforehand. But luckily the form of the 'split', which is defined by the condition O A G (rest) AB = 0 is independent of the normalization of O A . So we shall first determine another null geodesic field (let us denote it byŌ A to remind ourselves of the difference in normalization) which is affinely parametrized and whose inner-product with the normal to event horizon (which, upto normalization, could again be determined without any knowledge of the 'split') is one. Now we are at a stage to define the map the between the 'background' and the full space-time geometry.
Suppose {Y A } denote the coordinates in the background geometry (in our case pure AdS, the metric is denoted byḠ AB ) and {X A } are the coordinates of the full space-time (the dynamical black-brane, the metric is denoted by G AB ). Let us denote the invertible functions that give a one to one correspondence between these two spaces as {f A }.
The equations that will determine f A s are the followinḡ 12 HereŌ A are affinely parametrized the null geodesics in the full space-time geometries i.e.,Ō ·∇Ō A = 0 (3.8) Equation (3.8) would fixŌ A completely once we specify the angles it would make with the tangents of the horizons, which is effectively a set of (D − 1) numbers. Now what we are actually interested in is notŌ A but O A which is related toŌ A with a normalization. Therefore we are free to choose the normalization ofŌ A , since anyway we have to re-normalize it again. This will fix one of the (D − 1) initial conditions. Rest we shall keep arbitrary. We shall assumē A | horizon = some arbitrary functions of horizon cordinates (3.9) where N A is the null normal to the event horizon (with some arbitrary normalization) and l A (i) ∂ A s are the unit normalized space-like tangent vectors to the horizon. It turns out that the hydrodynamic metric could be split for a very specific choice of these spatial initial conditions and we shall fix them order by order in derivative expansion by matching the hydrodynamic and the 'large-D' metric. OnceŌ A is fixed (in terms of these arbitrary angles), we could determine f A s upto some integration constants by solving equation (3.7).
Equation (3.7) further says that if we apply the map (3.6) as a coordinate transformation on the 'background', then in the new {X A } coordinates the map would just be an 'identity' map and the full space-time metric G AB would admit the split as given in equation (3.1) satisfying (3.2) 13 . Once we have figured out how to split the full space-time metric into 'background' and the 'rest', we know how to view the event horizon as a surface embedded in the 'background' and therefore the auxiliary function ψ (by solving the harmonicity of ψ −D w.r.t the background) everywhere. Now we can normalizeŌ A as it has been done in [7]. Using these ψ and O A (appropriately normalized) one should be able to recast any arbitrary metric, that admits large-D expansion, exactly in the form of [7].

Bulk-Membrane map in metric dual to Hydrodynamics
In this subsection, we shall implement the above algorithm, described in the previous subsection, for the metric dual to hydrodynamics. For convenience we are summarizing the steps again.
• Determine the equation for the event horizon.
• Determine the null normal to the horizon.
• Solve equation (3.8) to determineŌ A everywhere. We need the normal, derived in previous step, to impose the boundary condition.
• Choose any arbitrary coordinate system {Y A }, where the 'background' has a smooth metric G AB .
• Now solve the equation (3.7) to determine the mapping functions f A 's.
For a generic dynamical metric, it is not easy to implement all these steps. But in this case what would help us is the 'derivative expansion' and the fact that f A 's are exactly known at zero derivative order; it is simply the coordinate transformation between Eddington-Finkelstein and Kerr-Schild form of a static black brane metric. 13 We would also like to emphasize that what we are describing here is not just a gauge or coordinate transformation. The 'split' mentioned in equation (3.1) is a genuine point-wise map between two entirely different geometries. Once we have figured out the 'map', we are free to transform the coordinates further; both G AB andḠ AB would change, but the 'map' will still be there.
Though the zeroth order transformation is already known, as a 'warm-up' exercise we shall re-derive it using the above algorithm. The condition of 'staticity' and translational symmetry of the metric allow us to solve relevant equations exactly in this case.

Zeroth order in 'derivative expansion':
At zeroth order in derivative expansion the metric dual to hydrodynamics has the following form We could read off the components of the metric and its inverse.
At zero derivative order both r H and u µ could be treated as constants, The event horizon and the null normal to it are given by Now we shall figure out the 'map' that will lead to the desired 'split' between 'background' and 'rest'. We have already determined the event horizon. Next we have to solve forŌ A , satisfying the conditions At zero derivative order, G AB has translational symmetry in all the x µ . The conditions onŌ A does not break this symmetry. HenceŌ A must have the form Now we shall process the condition that O A is a null vector field.
So finallyŌ A ∂ A = h 1 (r)∂ r 14 . Substituting this form ofŌ A in the geodesic equation we could see that h 1 (r) has to be a constant and then boundary condition simply says that h 1 (r) = 1 Now let us choose a coordinate system Y A = {ρ, y µ } for the 'background' where the metric takes the following form Again the symmetries motivate us to take the following form for the mapping, which gives the one to one correspondence between the background coordinates {Y A } = {ρ, y µ } and black-brane coordinates Let us apply the map (4.8) as a coordinate transformation on the background. In the new coordinates (where the map is just an 'identity') the background metric takes the following form Here we have suppressed the r dependence and derivative w.r.t r is denoted by prime ( ′ ). In this coordinates equation (3.7) takes the form These two equation could be solved very simply. The general solution where c 1 and c 2 are two arbitrary constants. We shall choose the plus sign in h(r) to make sure that whenever r increases, ρ also increases.
14 Actually there are two solution to (4.5). If we assume h 2 (r) = 0 and finite everywhere, then This implies that h 1 (r) will vanish at the horizon r = r H (which is a zero of the function f (r/r H )), contradicting the boundary condition onŌ r . Now we have to fix the integration constants. Note that once we know the map, we know the form of G (rest) AB , satisfying equation (3.2) by construction.
Now we further want that if D → ∞, the metric should reduce to its asymptotic form at any finite distance from the event horizon or in other words, G (rest) µν must vanish outside the 'membrane region' (a region with 'thickness' of the order of O 1 D around the 'membrane', see section (2.2)). This condition will force us to set c 1 = 0. The other constant c 2 is not appearing in the final form of the metric at all, so this ambiguity will remain here at this order and it is simply a consequence of the translational symmetry in x µ and y µ directions. For simplicity we shall also choose c 2 = 0. So the final form of the map at zeroth order ρ = r, y µ = x µ + u µ r (4.13)

First order in derivative expansion
In this subsection we shall extend the computation of the previous subsection upto the first order in derivative expansion. Here u µ and r H depends on x µ but any term that has more than one derivatives of u µ and r H has been neglected. All calculations presented in this subsection generically will have corrections at order O(∂ 2 ). At first order in derivative expansion the metric dual to hydrodynamics has the following form [5] Where, We shall often refer to this metric, described in equation (4.14), as 'hydrodynamic metric'. Here both r H and u µ are functions of x µ s; but they are not completely arbitrary. the hydrodynamic metric will solve the Einstein's equations (upto corrections of order O(∂ 2 )) provided the derivatives of r H and u µ satisfies the following equations 16 .
We read off the components of the metric and its inverse The horizon is still given by the surface (no correction at first order in derivative, though the normal gets corrected since ∂ µ r H is not negligible now.) We need the Christoffel symbols to compute the geodesic equation.
Γ r rr = 0, Γ µ rr = 0 Γ r αr = rf (r/r H ) + At first order in derivative expansion, the most general correction that could be added toŌ A , maintaining it as a null vector with respect to the first order corrected metric: We shall fix w 1 (r) and w 2 (r) using the geodesic equation. The r component of the geodesic equation gives the following.
From the µ component of the geodesic equation we find We could partially fix the integration constants using the boundary conditions. At horizonŌ Hence it follows that .
Next we have to solve for the 'mapping functions'. Let us choose the same coordinates {Y A }, as in the previous subsection so that background takes the form of equation (4.7). We expect that the mapping functions (4.13) will get corrected by first order terms in derivative expansion.
As before, we shall apply the map (4.25) as a coordinate transformation on the background. In the new coordinates (where the map is just an 'identity') the background metric takes the following form Substituting equation (4.26) in equation (3.7) we find The general solution for equation (4.27): In the new X A = {r, x µ } coordinates the metric of the background takes the following form In the last step we have rewritten G µν using the following identity Once we know the background, we could determineḠ rest AB . In this section we would like to expand the 'hydrodynamic metric' (already split into 'background' and 'rest' in the previous secion) in an expansion in 1 D and compare it against the metric described in [7].
This comparison involves two steps. The first one is of course an exact match of the two metric upto the required order. The second step involves the mapping of the evolution of the data. Let us explain it in a little more detail. As we have mentioned before, both 'hydrodynamic metric' and 'large -D' metric are determined in terms of data, defined on a co dimension one hypersurfaces -in the first case it is the velocity and temperature of a relativistic fluid living on the boundary of asymptotic AdS and in the second case it is the horizon viewed as a membrane embedded in the background with fluctuating shape and velocity. However we cannot choose the data arbitrarily. The hydrodynamic metric or the large D metric will solve the Einstein's equations only if the corresponding data satisfy certain evolution equation. For matching of these two metrics, the evolution of the data also should match. More precisely , we should be able to re express the membrane velocity and shape in terms of fluid velocity and temperature and further we have to show that once hydrodynamic equations are satisfied, the membrane equation is also true upto the required order.
Below we shall first compare the two metrics and in the next subsection we shall prove the equivalence of the evolution of these two sets of defining data.

Comparison between the two metrics
If the hydrodynamic metric has to match with the final metric described in [7], the first requirement is thatḠ rest µν must vanish as one goes finitely away from the horizon. This is possible providedC 3 is zero and also the function [F (r/r H ) − 1] has a certain type of fall-off behavior at large r. NowC 3 being an integration constant we could easily set it to zero. In appendix (A) we have analyzed the integral (4.15) and therefore the function [F (r/r H ) − 1]. It turns out that at large D this integral could be approximated as follows.
After substituting equation (5.1) and the value for the integration constantC 3 , the black-brane metric dual to hydrodynamics takes the following form 17 Also note that the vanishing has appropriate fall-off behavior (exponential decay in the scaled Z variable) as required by large D corrections where dS 2 background is given by equation (4.29) As we have mentioned before, the metric in [7] is described in terms of one auxiliary function ψ and one auxiliary null one-form O A dX A . For convenience we are quoting the metric here again.
Here ψ −D is harmonic with respect to the background with ψ = 1 being the event horizon of the full space-time and O A is simply proportional toŌ A determined in the previous subsection. The proportionality factor (let us denote it by the scalar function Φ(X)) is fixed using the condition that the component of O A along the unit normal of ψ = constant hypersurfaces is one everywhere. In terms of equations, the above conditions could be expressed as The metric in (5.5) will match exactly with the metric in (5.2) provided we set A 2 to zero and identify Φ 2 r 2 r H with the boundary condition that ψ = 1 should reduce to the horizon, i.e., the hypersurface given by r = r H , in an expansion in 1 D . Now we shall first determine ψ and then Φ. Note that both ψ and the norm of ∂ A ψ are scalar functions and it is much easier to compute them in a coordinate system where the background metric has a simple form. Therefore we shall solve the equation in the {ρ, y µ } coordinate system and then transform the answer to the {r, x µ } coordinates for final matching. First we need to know the position of the horizon in {Y A } coordinates since that will provide the required boundary condition for ψ. We know that in {X A } = {r, x µ } coordinates the horizon is at r = r H (x) + O(∂ 2 ). Now {X A } and {Y A } coordinates are related as follows.
The inverse transformation: Therefore in terms of {Y A } coordinates the horizon is at Here, for any term that is of first order in derivative to begin with, this coordinate transformation will generate change of order O(∂ 2 ) and therefore negligible in our computation. In the last line we have used equation (4.16).
Once we know the position of the horizon, we could solve for ψ. In {ρ, y µ } coordinates the expressions for ψ and its norm are as follows (see appendix (B for derivation).
Clearly this solution satisfies the boundary condition that ψ = 1 ⇒ ρ = r H (y) + O(∂ 2 ). Now we have to transform these quantities in {X A } coordinates. We shall first transform the quantity ρ r H (y) .
From equation (5.12) it follows that Substituting this solution in equation ( This is exactly what is required to have a match between the 'hydrodynamic metric' and the 'large-D' metric upto the expected order.

Comparison between the evolution of two sets of data
As mentioned before, the 'hydrodynamic metric' is defined in terms of the velocity and the temperature 18 of the relativistic conformal fluid moving in a flat Minkowski 18 The temperature and the horizon radius are related by the following relation space-time of dimension (D − 1). In case of large -D expansion, the metric is given in terms of a (D − 1) dimensional time-like fluctuating membrane embedded in pure AdS space-time with a dynamical velocity field on it. Both of these two sets of data are controlled by separate equations. For 'derivative expansion' , the governing equation of data is given in (4.16). In 'large-D' technique, the relevant equation is the following [7] ∇ · U = 0, Here the equation is written as an intrinsic equation on the membrane world-volume. All raising, lowering and contraction of the indices are done with respect to the induced metric on the dynamical membrane. U α is the velocity of the membrane, expressed in terms of its intrinsic coordinates. K βα is the extrinsic curvature of the membrane, expressed as a symmetric tensor on the membrane world-volume. K denotes its trace. P α γ is the projector perpendicular to U α . In this subsection, our goal is to show that equation Our first job would be to express the U α and K αβ in terms of velocity u µ and temperature (or r H ) of the relativistic fluid. Remember that though both u µ and U α are unit normalized velocity vector, they are defined on completely different spaces, one being a flat Minkowski metric and the other is the curved (both intrinsic and extrinsic curvature, being nonzero) membrane world volume. For convenience, we shall work in {Y A } = {ρ, y µ } coordinates where the background metric is simple. We shall first compute the unit normal to the membrane and different components of its extrinsic curvature, to begin with in terms of background coordinates and then we shall re-express it as an intrinsic symmetric tensor on the membrane. The unit normal to the membrane is given by The extrinsic curvature is defined as follows. Now let us choose {y µ } as the intrinsic coordinate on the membrane world volume. In this choice of coordinates, the extrinsic curvature K αβ will have the following structure.
Note that the first term in the RHS of equation (5.18) does not contribute at first order derivative expansion. After using equation (5.17) and (5.18), at this order the final expression for K µν turns out to be very simple (see appendix (C) for the details of the computation).
The induced metric on the membrane is given by Once we know K αβ , U α and the induced metric on the membrane, we could compute each term in the equation (5.15).
As it is clear from the notation, in the LHS of each equation the relevant metric is the induced metric on the membrane whereas in RHS it is the flat Minkowski metric η αβ . Substituting equations (5.24) in equation (5.15) we could easily show that membrane equation follows as a consequence of fluid equation.
In this context let us mention the work in [16]. Here the authors have computed the boundary stress tensor dual to a slowly varying membrane embedded in AdS. They have found the dual fluid velocity in terms of the membrane velocity. It could be easily checked that equation (5.23) is indeed the inverse of what they have found upto correction of order O(∂ 2 ).

Conclusion
In this note we have compared dynamical black-brane solutions to Einstein's equations (in presence of negative cosmological constant) generated by two different perturbative schemes, namely 'derivative expansion' and Large-dimension expansion. In both the cases, the space-time necessarily have an event horizon. We have shown that in large number of dimensions whenever 'derivative expansion' is applicable, we can expand the metric further in 1 D , (though the reverse may not be true always). We have found perfect match in this overlap regime of these two perturbative techniques upto first subleading order on both sides. One immediate interesting project would be to extend this calculation to the next order on both sides, since we already know both the 'hydrodynamic metric' and the 'large D metric' upto the second subleading order [2,18]. It would also be interesting to generalize this calculation to Einstein-Maxwell system in presence of negative cosmological constant, where also we know the metric on both sides upto the first subleading order [21][22][23][24].
In some sense, our analysis serves as a consistency test for these two methods. But this comparison could teach us something more. This is about the dual systems of these two gravity solutions. The dynamical black-brane metric generated by 'derivative expansion' in D dimension is dual to the relativistic conformal hydrodynamics living in (D −1) dimensional flat space-time. The variables of hydrodynamics are fluid velocity and temperature, which are the data that label different black-brane solutions in derivative expansion. On the other hand the metric generated in 'large D expansion' is dual to a codimension one dynamical membrane embedded in pure AdS and coupled with a velocity field. Here also the labeling data of the metric live on a (D − 1) dimensional hypersurface and they consist of a scalar function -the shape of the membrane and a unit normalized velocity field. This is very similar to hydrodynamics in terms of counting, though the governing equations and the physical significance of the variables are entirely different.
However, we have already seen that these two systems of equations are approximately equivalent after an appropriate field redefinition. In this note, we have verified it at the very leading order and we expect that the project of comparing the two metric upto second subleading order would extend this equivalence to the next order on both sides.
In fact it is expected that this equivalence is valid to all orders [16]. In other words, in the overlap regime, these two equations must be exactly equivalent to each other if we consider all orders on both sides [16], though to see this equivalence we need to re-express the variables of one side in terms of the other [16,19,25]. This equivalence actually involves some interesting resum of one series into the other. Even the leading term in derivative expansion can encode many terms of 1 D expansion and on the other hand the leading membrane equation might have information about many higher order transport coefficientas. At linearlized level, this has been nicely captured in the analysis in [26]. The frequencies of Quasi normal modes do exhibit such resum. In [16], the authors have proposed a resummed stress tensor that could exactly reproduce the fluid stress tensor exactly upto the first order in derivative expansion. It would be very interesting to understand this structure in full detail, at non linear level. This might lead to a fluid-membrane duality in large number of dimensions where gravity does not have any role to play.

A Analysis of F (r/r H )
In this section we shall evaluate the integral (4.15) in large D limit. For convenience we are quoting the equation here.
We would like to evaluate this integral systematically for large D. Let us first evaluate the integral for y ≥ 2. In this case, since D is very large, x D >> 1 throughout the range of integration. So we shall expand the integrand in the following way. Clearly the sums in the RHS of (A.3) are convergent for y ≥ 2. Let us denote the RHS as k(y). However, the expansion in (A.2) is not valid inside the 'membrane region', i.e., when y − 1 ∼ O 1 D and naively k(y) is not the answer for the integral. But consider the functionk(y) = F (y)−k(y). This function vanishes for all y ≥ 2 and also by construction it is a smooth function at y = 2 (none of the derivatives diverge). Hencek(y) must vanish for every y. So we conclude, for every allowed y (i.e., y ≥ 1) F (y) = 1 + Note that F (y) reduces to 1 as y → ∞ as required in section (4.2). Now we would like to expand F (y) in a series in 1 D , where y is in the membrane regime.
In this regime F (y) takes the following form In this note we considering only the first subleading correction in 1 D expansion. Therefore F (y) could be set to 1 for our purpose.

B Derivation of ψ in {Y
In this section we shall give the derivation of ψ as mentioned in eq (5.7). We want to solve ψ such that ∇ 2 ψ −D = 0. Where ∇ is the covariant derivative with respect to the background metric = Π ρ µ ∂ ρ n ν − Γ ρ ρν n ρ + Π α µ ∂ α n ν − Γ ρ αν n ρ = −δ α µ Γ ρ αν n ρ = ρ 2 η µν (C.3) Now, as mentioned in (5.18) in terms of the intrinsic coordinates on the membrane the extrinsic curvature will have the structure The trace of the extrinsic curvature For the calculation of only the extrinsic curvature we need background metric, where for the rest of the calculation we require induced metric on the horizon. The induced metric on the horizon is given by