Surface transport in plasma-balls

We study the surface transport properties of stationary localized configurations of relativistic fluids to the first two non-trivial orders in a derivative expansion. By demanding that these finite lumps of relativistic fluid are described by a thermal partition function with arbitrary stationary background metric and gauge fields, we are able to find several constraints among surface transport coefficients. At leading order, besides recovering the surface thermodynamics, we obtain a generalization of the Young-Laplace equation for relativistic fluid surfaces, by considering a temperature dependence in the surface tension, which is further generalized in the context of superfluids. At the next order, for uncharged fluids in 3+1 dimensions, we show that besides the 3 independent bulk transport coefficients previously known, a generic localized configuration is characterized by 3 additional surface transport coefficients, one of which may be identified with the surface modulus of rigidity. Finally, as an application, we study the effect of temperature dependence of surface tension on some explicit examples of localized fluid configurations, which are dual to certain non-trivial black hole solutions via the AdS/CFT correspondence.


Introduction
The theory of hydrodynamics provides us with a tractable effective framework to study the low-energy near-equilibrium states in any finite temperature system with a well behaved microscopic description. Although the description of these states in terms of the microscopic degrees of freedom may be very complicated, hydrodynamics allows us to describe JHEP06(2016)015 them with a few effective degrees of freedom -the fluid fields. This effective theory is constructed purely on the basis of symmetries inherent to the microscopic theory, in addition to certain empirical assumptions like the second law of thermodynamics (see [1] 1 ). The information of the underlying field theory is encoded, in a phenomenological way, in the transport coefficients that characterize the effective macroscopic description. 2 Although hydrodynamics is a very old and well studied subject, recently there has been a renewed interest in it, particularly after the discovery of its connections with black hole dynamics in the context of the AdS/CFT correspondence [7]. These recent studies have valuably contributed to the improved understanding of the structural features of the subject and has led to the unraveling of new transport phenomenon [2,3,[8][9][10][11]. Most of these latest developments mainly focus on the description of states in which the same phase of the fluid fills the entire spacetime, which is taken to be a non-compact pseudo-Riemanian manifold. In other words, the hydrodynamics that has been explored in most of these recent developments is the effective theory of a class of states, which does not involve any fluid surface or a phase boundary.
In this paper we proceed to analyze the situation where the class of states described by the effective theory of hydrodynamics is extended to incorporate the states that include a fluid surface, which separates two phases. We will mainly focus on finite lumps of relativistic fluids, which occupies only a finite subspace of an otherwise non-compact spacetime. One of the concrete examples of our set up is described in [12], where metastable finite lumps of the deconfined phase of N = 4 SYM is separated from the confined phase by a phase boundary. In the large N limit, such a situation can be described by a metastable fireball of plasmafluid separated from the vacuum by a fluid surface. 3 Although we have a set up similar to [12] at the back of our mind, we wish to clarify that in this paper we have taken a purely field theory perspective and we make no use of the AdS/CFT correspondence in any way. 4 We would like to highlight the fact that the behaviour of the microscopic field theory, at or near the surface, can in general be quite different from that in the bulk. This difference would be captured by new surface transport coefficients in the effective theory. Some of these new surface transport coefficients would simply encode the way in which the bulk transport coefficients are modified at the surface, while others would represent entirely new transport properties, particular to the existence of the fluid surface.
A very well known example of such a surface phenomenon, at the leading order in derivative expansion, is surface tension. In this paper, we study the surface transport coefficients, at the subleading order, and investigate the relations that may exists between them and the bulk transport coefficients. We would like to emphasize, that the surface transport coefficients carry entirely new information about the microscopics and modify the 1 See also [2][3][4][5][6] for a more recent use and discussion of the second law of thermodynamics in the context of hydrodynamics. 2 These transport coefficients are often expressible in terms of correlators of symmetry currents (Kubo formulae), which may be evaluated directly from the microscopic quantum field theory. 3 Although we have this specific set up in mind, our constructions can be straightforwardly generalized to describe the surface transport properties of any phase boundary. 4 On several occasions, in this paper, we use the word 'bulk' which would always mean the bulk of the fluid in contrast to its surface. It should never be taken to imply the holographic dual.

JHEP06(2016)015
fluid equations at the surface very non-trivially. For instance, knowledge of the equation of state 5 for a fluid, tells us nothing about how the surface tension depends on temperature.
The study of surface transport has been carried out to some extent in the context of fluids which are confined to a thin submanifold of spacetime [13][14][15]. This is the case in which the bulk fluid is not present or, alternatively, its pressure and higher-order transport coefficients vanish at the surface. In such situations, a systematic analysis based on an effective action approach [13], underlying symmetries and positivity of the entropy current [14,15] have been used to constrain the form of the equations of motion up to second order in derivatives. One of the particular features inherent to the study of dynamical surfaces is that, due to possible deformations along directions transverse to it, new transport coefficients appear encoding the response of the surface to bending. One such transport coefficient is the surface modulus of rigidity [13]. 6 The novelty in this paper resides on the fact that the system we study is a more intricate one, in which both the fluid residing in the bulk and the fluid living on the surface constitute the same system.
One of the central simplifying assumptions that we shall make in this paper is that we will only consider stationary fluid configuration. For the case of stationary space-filling relativistic fluids without any surface, the constitutive relations and hence the equations of motion could be significantly constrained from a simple physical criterion. This criterion is that we demand the symmetry currents, including the stress tensor, to follow from an equilibrium partition function [19][20][21][22][23].
The stationary fluid configurations are considered in the presence of non-trivial background fields, like the metric and the gauge fields corresponding to other conserved charges. These background fields are considered to be slowly varying, with respect to the length scale associated with the radius of time-circle, in this finite temperature description. These slowly varying background fields serve as sources in the relativistic fluid equations. On solving these equations, the fluid variables are expressed in terms of these background sources. Now, if we substitute this solution of the fluid variables, back into a putative action for stationary configurations, we obtain the equilibrium partition function expressed as a functional of background fields. Since, the fluid equations, and hence the solutions of fluids fields, are constructed in a derivative expansion, therefore the partition function could be expanded in a derivative expansion, in terms of these background fields and their derivatives.

JHEP06(2016)015
It is straightforward to compute the symmetry currents from the partition function by varying it with respect to the background fields. Then, we proceed to compare the symmetry currents so obtained, with that which is expressed in terms of the fluid variables through the constitutive relations. This comparison not only yields the expressions for the fluid variables in terms of the background fields, specific to the stationary configurations under consideration, but also provides non-trivial relations between the transport coefficients (see section 1.1 for more details). Although derived by analyzing stationary fluid configurations, these relations between transport coefficients are expected to hold even away from equilibrium. 8 In this paper, one of our principal goals would be to adopt a such a method to constrain transport properties at the surface of relativistic fluids. (1.1) Here, the metric functions, depends only on the spatial coordinates x, and ∂ t is our timelike Killing vector. Here g ij is the metric on spatial manifold obtained by reducing on the time-circle, 10 which we shall denote by N s . In some of our discussions, we will also include a conserved global U(1). The background gauge fields for this U(1) take the form Since none of the functions depend on time, all the quantities of interest, including the conserved currents, can be dimensionally reduced on the time-circle, whose radius we take to be T 0 . It is possible to restrict to this reduced language, focusing only on N s , for most of the discussions relating to the partition function. Among the reduced quantities, time-translation invariance survives as a gauge invariance corresponding to the Kaluza-Klein gauge field a i (x). All our constructions starting from the partition function, must be manifestly invariant under this Kaluza-Klein gauge transformation.
Since the U(1) gauge fields in (1.2) transform non-trivially under the Kaluza-Klein gauge transformation, it is convenient to define a new set of shifted gauge fields which are 8 These constraints were found to be identical to the equalities among transport coefficients that follow from the considerations of the second law of thermodynamics. See [24,25] for more details on this connection. 9 The discussion in this subsection is generally applicable in all dimensions, but while performing the analysis, particularly in section 3, we shall specialize to four dimensions. Also we shall choose the Levi-Civita connection to define the covarinat derivative on N . 10 This can be done by identifying all the points on the orbits generated by the time-like Killing vector.

JHEP06(2016)015
manifestly invariant under it 11 , T 0 } constitutes the set of background data in terms of which the partition function is to be expressed. Now, since we wish to describe a finite lump of relativistic fluid, we will assume that the fluid is confined to a sub-manifold of N of the same dimensionality as the spacetime, which we shall denote by M. The fluid surface is considered to be a co-dimension one hypersurface. We shall denote the fluid surface by f (x) = 0, where f (x) is taken to be independent of time, following our stationary assumption. In fact, f (x) can be taken to be a spatial coordinate itself, which is positive inside the fluid and negative outside. The region inside, f (x) = 0 is M, which can again be reduced on the time-circle to obtain M s . Here M s is also a sub-manifold of N s , with a compact boundary. We furthermore assume that the boundary of M does not have boundaries itself.
The normal vector orthogonal to the fluid surface is a spatial vector, with a vanishing inner product with the time-like Killing vector. The partition function of interest, after performing the trivial time integral, can be schematically written as where θ(f ) is a distribution functional of the surface function f (x) and is introduced to encode the variation of the fluid fields at the surface. Here, θ(f ) contains all the information of the surface. In particular, it has a dependence on the dimensionless parameterτ = T 0 τ , with τ being the length scale associated with the surface thickness. All such non-universal dependence of θ(f ), which are sensitive to the microscopics, are left implicit throughout our analysis. Realistically θ(f ) is a distribution as shown schematically in figure 1. The notation θ(f ) is purposely used to indicate the fact that, in the limiting case where the parameterτ is small, this distribution may be well approximated by a Heaviside step function. 12 We will also assume that the size of the fluid configuration, i.e. the average length scale associated with M s , is much greater than τ as well as that associated with the temperature. We would like to expand S in a derivative expansion. Keeping in mind the reparameterization invariance of the surface, this derivative expansion can be schematically performed in the following way (1.6) 11 Note that our definition of A0 here differs from that in [19], in that we do not include the shift with respect to the chemical potential, which may be thought to have been absorbed in the constant part of A0. 12 Besides the θ(f ) andδ(f ), the surface transport coefficients may also be dependent onτ . Here we shall also leave that implicit (see section 1.2 for more details).

JHEP06(2016)015
where S j (i) is a collection of terms containing i number of derivatives on the background fields B. S 0 (i) are the terms in the bulk of the fluid and they are exactly the ones that have been considered for space filling fluids, in the earlier constructions of stationary fluid partition functions [19,[21][22][23]26]. The dots at the end, in (1.6), denote terms where more than one derivatives act on θ(f ). On all such terms, an integration by parts can be performed and they can be cast into the same form as the second term in (1.6). On performing such an integration by parts, the modified S 1 (i) now must contain terms involving derivatives of the normal vector n. These new kind of terms specific to the existence of the surface are simply the ones involving the extrinsic curvature of the surface and its derivatives. Indeed, these are the type of terms considered in the analysis of effective actions for fluids confined to a thin surface [13][14][15]. We may write, without any loss of generality 13 Here, δ(f ) denotes the derivative of the distribution θ(f ). This notation is again purposely chosen, so that, in the limit where θ(f ) approximates to Heaviside step function, δ(f ) approximates to the Dirac delta function. Finally, if we can reliably approximate θ(f ) and δ(f ) to the Heaviside step function and the Dirac delta functions respectively, we may write (1.7) as 14 The second term in (1.9) is the main focus of this paper, and in particular cases, we shall provide the explicit forms of the surface partition function, up to the first non-trivial orders in derivatives. Upon variation of the partition function (1.7) with respect to the background metric, the stress tensor that we get has the form Note that although we were able to remove the derivatives of delta function in the partition function by an integration by parts, such derivatives are still present in the expression for the local stress tensor. The remaining symmetry currents also have a structure similar to (1.10).
There is a very important and interesting role played by the function f (x) in the partition function (1.7). We may derive an equation of motion for f (x) by extremizing the partition function with respect to it. We can think of this equation as the one which 13 Note that reparameterization invariance of the surface fixes the dependence on √ ∂f · ∂f . Therefore any other additional dependence on this quantity has not been considered.
14 Here we have to use the fact, √ ∂f · ∂f = γ g , where γ is the determinant of the induced metric on ∂Ms, γij = gij − ninj.

Fluid variables and choice of frames
Now let us consider the description in terms of the original fluid variable {u µ , T, . . . }. We can write down a stress tensor in terms of these fluid fields and their derivatives purely based on symmetry considerations, which has the same structure as (1.10), namely The transport coefficients are the functions of scalar fluid fields, that multiply specific symmetry structures whenT µν (i) is further expanded in a derivative expansion. The number of scalar structures that goes into the partition function is, in general, much less than the allowed linearly independent symmetry structures which arises in the stress tensor (1.11). Therefore comparing T C µν with T P F µν gives us very non-trivial constraints on the transport coefficients. This exercise, when executed for the fluid configurations with a surface, not only gives relations among surfaces transport coefficients but also relates them with some of the bulk transport coefficients.
While writing (1.11), there is a crucial issue of the choice of fluid frames. Since we do not want the fluid to spill out of the surface, we must require 12) to be true at all orders in derivative expansion. Also, we want the surface to be moving in the same way as the bulk of the fluid. Therefore a suitable frame choice should ensure that the fluid fields does not jump discontinuously at the surface. We should point out

JHEP06(2016)015
that one of the most popular frame choice -the Landau frame, in which u µ is an eigenvector of the full stress tensor, is not a suitable frame choice in this respect. This is because since the stress tensor has new corrections at the surface, the Landau frame choice introduces discontinuities in the fluid variables at the surface. Also, trying to impose the condition (1.12) in addition to the Landau frame condition may turn out to be more constraining than necessary. A suggestion for a suitable frame choice may be to work with a Landau frame condition only for the bulk stress tensor. 16 That is, we may impose both in the bulk and at the surface of the fluid. Here E is only the bulk energy density. This, in particular, would imply that the direction of energy flow at the surface of the fluid is not along the fluid flow, which is not a problem at all. However, although this cures one of the problems (that there would be no separate corrections to the fluid fields at the surface), imposing (1.12) in addition may still be over constraining, particularly at the second and higher orders in derivatives. This problem may be easier to visualize, if we remember the results of [19] for 3+1 dimensional uncharged fluids. In that paper, after comparison with the partition function, in the Landau frame, it was found that the velocity u µ , which is identical to the Killing vector of (1.1) at the leading order, receives nontrivial corrections in terms of derivatives of the background fields at second order in the derivative expansion. So if we now wish to impose (1.12) on that result, it would be automatically satisfied at the leading and first order. But at second order, it would imply non-trivial constrains involving the background field and the normal vector n µ , which may be too restrictive.
Hence, the most suitable choice of frame for this problem is to choose a frame where (1.12) is a part of the frame fixing condition. This is possible to implement only because, at leading order the fluid velocity u µ is proportional to the time-like Killing vector of (1.1) therefore (1.12) is automatically satisfied. The higher order corrections can always be manipulated by a frame choice. In fact, we may foliate the bulk of the fluid M, with constant f (x) surfaces, thus extending n µ throughout M. Then u µ n µ = 0 can be chosen to be part of frame fixing condition throughout the bulk of the fluid. 17 The remaining part of the frame choice can be implemented by imposing a condition similar to (1.13) but projected orthogonal to n µ . We will refer to this frame as the orthogonal-Landau-frame. 18 While performing the partition function analysis in section 3.3, we shall make this frame choice.
There is another very interesting point of view while describing the fluid surface in terms of the fluid variables. We may consider two separate sets of fluids variables one in the bulk {u µ , T } blk and the other at the surface {u a , T } sur . The surface has one less fluid variable because it is one dimension less than the bulk. Now we can regard the surface JHEP06(2016)015 equation of motion, 19 as being the dynamical equations for the variables {u a , T } sur . In this equation, the bulk variables only act as sources and we can solve them to obtain specific solutions {u s a , T s } sur . Subsequently we should solve the bulk fluid equations with the boundary condition 20 where e µ a is the projector on to the tangent space of the surface. Thus, in this point of view, the fluid equations are solved with dynamical 21 boundary conditions. The dynamics of the boundary conditions are given by the surface equations of motion. In this paper, we strive to constrain the form of this surface equation using the framework of the equilibrium partition function.

A brief summary of results
At first, in section 2 we consider perfect fluids in arbitrary dimensions. The partition function for perfect fluids with a surface can be written as where just like P can be identified with the pressure P in the bulk of the perfect fluid, C is identified with the surface tension χ. Comparison of the stress tensor constructed through symmetry arguments, with that following from (1.15) yields the expected surface thermodynamics. The component of the stress tensor conservation equation normal to the surface, at the surface, reads where K is the extrinsic curvature of the surface and a µ is the fluid acceleration. This is a modified version of the Young-Laplace equation where the term proportional to the acceleration is new compared to its original form. This additional term is non-zero only if the surface entropy χ S is non-zero, i.e., if there is a non-trivial temperature dependence of surface tension. This term has a very simple physical interpretation. If the surface entropy is non-zero, it implies that there are non-trivial degrees of freedom localized at the surface. Then the additional term accounts for the centripetal acceleration of these degrees of freedom in the force balance equation that (1.16) represents.
Since the acceleration term in (1. 16) has not been widely considered in the literature before, we analyzed the consequence of this term on some simple fluid configurations in 19 By this we mean the part of the equation of motion proportional toδ(f ) and its derivatives. Note that the number of fluid equations at the boundary is equal to the number of dimensions M rather than ∂M, which is one higher than the number of fluid variables {ua, T }sur. The extra equation, may be thought of as the equation of motion for the function f (x). 20 This boundary condition should be implemented at all orders in derivatives. 21 Here the word 'dynamical' is to interpreted only in a restricted sense, since we only consider a stationary equilibrium ansatz. Some dynamics is still present in our stationary assumption, in contrast to a static one.

JHEP06(2016)015
section 4. For this purpose, as a sample system, we choose to revisit the localized configurations of deconfined plasma of large N , strongly coupled N = 4 Yang-Mills theory, compatified down to 2 + 1 dimensions on a Scherk-Schwarz circle, that were constructed in [27]. These configurations are dual to exotic black holes in Scherk-Schwarz compactified AdS 5 . However, due to the unavailability of the exact dependence of surface tension on temperature for this system, the surface tension was taken to be a constant, 22 in the analysis of [27]. In section 4, we suitably parameterized this ignorance and studied the change in the phase diagram for the configurations, as we varied this parameter. We found that turning on this parameter introduced an upper bound on the surface velocity. This arises from the fact that the surface temperature dips below the phase transition temperature, when the bound is overshot. This results in the termination of the phase curve at a specific point, see figure 2. This has important consequences for the existence of a phase transition between the ball and ring configurations. We find that the phase transition may not exist for large values of this parameter. 23 Moving on to the case of finite lumps of superfluids, at zeroth order in derivatives, P and C in the partition function (1.15), now would also depend on A 0 and the norm of the superfluid velocity ξ µ [23]. We find that (1.16) is further modified in the case of superfluids to become (1.17) In section 3 we consider the case of uncharged fluids in 3+1 dimensions, where the first corrections to the perfect fluid partition function occurs at second order in the bulk and at first order on the surface. The full partition function upto this order, including the parity odd sector, takes the form Here P i are the three independent coefficients that were considered in [19], while B i are the new three surface transport coefficients. The terms proportional to B 1 and B 2 can also be viewed as bulk total derivative terms, while the term B 3 , a term which was studied in [13], eventually contributes to the modulus of rigidity. As pointed out before, the surface transport coefficients in (1.18) may also depend oñ τ (the dimensionless ratio of surface thickness and T 0 ), which we leave implicit here. If any particular limit is taken on this parameterτ , it may directly influence the surface transport coefficients in (1.18), particularly B 3 .
For the stress tensor T µν C , which follows from symmetry considerations, there are 31 surface terms that can be written down, which are linearly independent for stationary 22 The value of surface tension at the phase-transition temperature which was previously computed in [12] was used for this constant value. 23 More specifically, if χ ∼ T −ζ , then we found for ζ 1/3 the phase transition would cease to exist between the ball and the ring.

JHEP06(2016)015
configurations. Now, taking into consideration the fact that P i correspond to three independent transport coefficient in the bulk fluid, we are able to derive 28 relations between the 31 surface transport coefficients and the 3 independent bulk transport coefficients, as it has been explicated in section 3.

Perfect fluids
In this section we will study fluids at zeroth order in derivatives. At this order, the only surface effect is encoded in the surface tension, which is extremely well studied. However, it is very instructive to re-derive the known physics in the language of partition functions described in section 1. We will also get the occasions to discuss a few effects related to the temperature dependence of surface tension and surface tension in superfluids which has not been widely discussed in the literature.

Ordinary uncharged perfect fluids in arbitrary dimensions
At first, let us briefly review the partition function for space filling ordinary perfect fluids as discussed in [19]. The partition function in terms of the metric sources can be written as, The functional form of P is to be determined from microscopics. Let us now evaluate the stress tensor from the above partition function by using [19] T Evaluating these formulae explicitly for (2.1) we get where T = T 0 e −σ . By comparing (2.3) with the zeroth order form of the stress tensor that follows from symmetry considerations we get while the fluids fields are found to be Note that (2.5) is identical to the condition on pressure and energy density that follows from thermodynamics. In this way, we are able to derive the thermodynamic properties of the fluid by comparison with the partition function. In some sense, pressure and energy

JHEP06(2016)015
density can be thought of as zeroth order transport coefficients, which are related by thermodynamic relations which follow from the partition function analysis. Following the above procedure, we wish to write down a partition function for perfect fluids in equilibrium, confined within a surface (which itself is dynamically determined by minimization of the free energy). Respecting the principles of KK-gauge invariance for writing down the partition function and reparameterization invariance of the surface, the partition function is given in terms of two unknown functions 24 In order to obtain the stress tensor, we have to vary the partition function (2.8) with respect to the background metric fields. 25 In fact, using (2.2), explicitly we find Now we have to compare (2.9) with the stress tensor that may be written from symmetry arguments (1.11) to this particular order, namely, The bulk stress tensor components T (0) µν are given by (2.4), while the components of the surface stress tensor also have a similar form Here, P µν = G µν + u µ u ν − n µ n ν is the projector orthogonal to both the velocity vector and the normal to the surface. In (2.11) we also introduced χ E and χ, which are, respectively, the surface energy and the surface 'pressure', also known as surface tension. Comparing (2.9) with (2.11), we can therefore identify while the fluids fields are again given by (2.6). 26 This identification (2.12) had also been done in [13]. Just like the bulk perfect fluid, (2.12) implies the thermodynamic identity

13)
24 Note that, as described in section 1, under suitable assumptions on θ(f ) and δ(f ), (2.8) may also be written as, (2.7) 25 The functional W is to be taken to be a functional of the metric functions and the function f , defining the surface. All these functions are independent functions of the coordinates and must be treated independently. 26 Note that with this the continuity of the fluids fields, as we move from the bulk to the boundary is maintained. Also, (2.6) implies u · n = 0, is automatically satisfied, at the order of perfect fluids.

JHEP06(2016)015
where χ S = −∂χ/∂T is the surface entropy. Therefore, if the surface tension depends non-trivially on the fluid temperature, it means that the surface entropy is non-zero and hence that there are active degrees of freedom living on the surface of the fluid. When the surface entropy vanishes (that is the surface tension is constant), the surface tension is equal to the negative of surface energy. The conservation of the stress tensor (2.10), implies where we have defined∇ µ (. . . . In the bulk (2.14) will give rise to the usual equation bulk conservation equation ∇ µ T µν (0) , while at the surface, it gives rise to the condition This equation is in fact a Carter equation with a force term [28], and is extensively used in the context of (mem)-brane hydrodynamics [13][14][15]29]. As with any Carter equation, (2.15) gives rise to two physically different sets of equations, obtained by projecting both orthogonally and tangentially to the fluid surface. Before doing so, let us note that that is, for a perfect fluid, the bulk contribution to (2.15) is only present in the normal component of (2.15). Thus, if we project (2.15) along the fluid surface with the projector e a ν such that n ν e a ν = 0, we obtain e a ν∇ µ T µν Here the index a = 0, 1, 2 labels the directions along the surface. Equation (2.17) expresses the conservation of the surface stress tensor along the surface. Note that if we consider higher derivative corrections, this equation will also receive a contribution from the bulk stress tensor, which would signify energy and momentum transport from the bulk to the surface. For the perfect fluid, however, such transport does not take place. The component of (2.15) normal to the surface, describing the elastic degrees of freedom of the surface, is more interesting and provides us with the condition that determines the position of the surface. For perfect fluids, it reduces to Now, given that the surface stress tensor (2.11) is orthogonal to the normal vector n µ , we can rewrite the second term in (2.18) in the following way

JHEP06(2016)015
where a µ = u ν ∇ ν u µ is the fluid acceleration, and K = G µν K µν is the mean extrinsic curvature. If there are no active degrees of freedom in the boundary and the surface tension is constant implying that χ E = χ, then (2.20) immediately reduces to which is the Young-Laplace equation. Thus (2.20) is a generalization of the Young-Laplace equation, when there are non-trivial degrees of freedom living on the surface of the fluid. Such generalization had not been previously considered in the works of [27,30,31]. In section 4, we shall examine the consequence of the presence of this additional term in (2.20) for some simple fluid configurations. It is noteworthy that (2.20) can also be directly obtained from the partition function (2.8). In terms of the partition function, this is simply given by the extremization of the partition function with respect to the surface function f , which is the equation of motion f . This is intuitively expected since the location of the fluid surface is obtained by the minimization of the free energy. If we vary (2.8) with respect to f , at the leading order we find where K = g ij∇ i n j , with∇ being the spatial covariant derivative defined with respect to the reduced metric. 28 Given the thermodynamical relations (2.5), (2.12) and remembering that in terms of the background fields the fluid acceleration is given by (2.20) reduces to (2.22).

Zeroth order superfluids with a surface
For the case of superfluids, the zeroth order stress tensor is modified in order to include the superfluid velocity ξ µ [32,33]. The bulk stress tensor has the form and it is accompanied by a conserved current Just like in the case of ordinary fluids, in the presence of a surface there are surface stress tensor and current contributions, which read It can be explicitly checked that this form of the boundary stress tensor and current follows from the following partition function

26)
28 K is related to the full extrinsic curvature by K = K + n · ∂σ.

JHEP06(2016)015
where ξ is the norm of the superfluid velocity. The bulk term of this partition function was first derived in [23]. In the presence of the surface we must also assume that the superfluid velocity and the normal to the surface are mutually orthogonal n µ ξ µ | f =0 = 0. As in the previous sections, the stress tensor is conserved and the normal component of the conservation of the boundary stress tensor gives the generalized Young-Laplace equation.
To leading order, the bulk and surface currents are conserved separately, with no current flowing from the bulk to the surface. The generalized Laplace-Young equation (2.20), is further modified in the case of superfluids. It takes the form Note that the new term is present even if there is no temperature dependence in the surface tension, as long as the goldstone boson also constitutes an active degree of freedom on the surface. Also this modified equation is also applicable to the case when there is an emergent goldstone boson only at the surface of a fluid, a situation which is reminiscent of topological insulators in the context of fluids.

Next to leading order corrections for uncharged fluids
In this section we shall consider the next to leading order corrections for uncharged fluids with a surface. The principal goal of this section is to demonstrate that there are only three new equilibrium transport coefficients on the surface of the uncharged fluid, at the next to leading order, two of which are parity even and the other one being parity odd. Two of these new boundary terms in the partition function, also precisely coincide with two possible bulk total derivative terms. Here, we work out the interplay between these new surface transport coefficients and the bulk second order transport coefficients.

Partition function at next to leading order
In order to write down the first corrections to the partition function (2.8), we need to write down all KK-gauge invariant scalar terms at higher order in derivatives. As it was observed in [19], the bulk of the fluid does not receive any corrections at first order. In other words, there are no KK-gauge invariant scalar bulk terms at first order, which can be written in terms of the sources of an uncharged fluid. However, at the surface of the fluid, there is an additional geometrical structure, which is the vector normal to the fluid surface. This allows us to write down two possible parity even scalar terms which may constitute the partition function. These are Here K is the trace of the extrinsic curvature reduced along the time direction Now for parity odd terms, there is only one possible parity odd scalar

JHEP06(2016)015
which must be taken into account while writing down the partition function. For spacefilling fluids, it was not possible to construct any parity odd term in the partition function, at the second order [19]. Therefore the existence of this term suggests that even uncharged fluids may have parity odd transport when surface effects are considered. Since, the fluid surface is co-dimension one, this term is particularly reminiscent of a parity odd transport that can exist in 2 + 1 dimensions [19,20]. Including these new surface terms, the partition function takes the following form It is important to note that if we include these first order surface terms in the partition function then we must also include the second order bulk terms for consistency. For instance, the surface stress tensor following from the partition function in (3.3) may get contributions from total derivative terms at second order in the bulk. In fact, we must point out that two of the new terms that we have added in (3.3) may also be written as a bulk total derivative terms at second order. Thus, including the bulk second order terms, which was written down in [19], 29 we have where P i denote the three independent transport coefficients at second order for a fluid without surfaces. Now, as pointed out before we can write the new surface terms as a bulk term in the following way where we have defined Now, as it was shown in [19], P 1 , P 2 and P 3 were determined in terms of the bulk second order equilibrium coefficients. Also eliminating the P 1 , P 2 and P 3 from those relations, gave rise to 5 relations among the eight possible second order equilibrium transport coefficients.

JHEP06(2016)015
It is clear that the terms proportional to P 4 and P 5 (or any total derivative term) will not enter the bulk stress tensor. 30 But they definitely contribute non-trivially to the surface stress tensor. It is important to note that both the form of the partition function (3.4) and (3.5) are equivalent and describe the same system. Hence, everything physical that is evaluated from them, such as the surface stress tensor or the Young-Laplace equation, must be identical.

Corrections to the stress tensor
Once we have written down the partition function (3.5), it is immediate to evaluate the stress tensor by varying with respect to the background fields using (2.2). It is convenient to split the bulk and surface contributions, up to second order, in the followind way The bulk stress tensor remains the same as that computed in [19] and involves only the the coefficients P 1 , P 2 and P 3 while the surface contribution is obtained as terms proportional toδ(f ) and its derivatives 31 This stress tensor must satisfy the conservation equation which gives rise to two separate sets of equations as in section 2, one determining the position of the surface and the other the conservation of the surface stress tensor along surface directions. Indeed, by explicitly using (3.9), one can verify that the tangential projection of (3.10) is automatically verified -a trivial consequence of diffeomorphism invariance along the surface. 30 This is because, a bulk total derivative can always be written as

7)
for any Vi. Since the variation of such terms, with respect to the metric, always lies within the derivative, hence such terms can only contribute to the surface stress tensor and never to the bulk stress tensor. 31 Here we have kept the term proportional to ∇inj − ∇jni. This term depends on how ni is extended away from the fluid surface. We may choose to perform this extension so that this anti-symmetric derivative of ni is zero. However, we perform our analysis without such an assumption so that, if our equations is applied in some generalized circumstance where a more natural extension of ni away from the surface demands this term to be non-zero.  Table 2. The first order vectors on the surface projected on the surface and orthogonal to the velocities. Here P µν = G µν +u µ u ν −n µ n ν , is the projector orthogonal to both u µ and n µ . The spatial components of this projector, projects orthogonal to n i in a given time slice, that is P ij = g ij −n i n j . Also, note that in the reduced language V 0 (j) = 0.

Constraints on surface transport coefficients from the equilibrium partition function
In this section we write down the next to leading order surface stress tensor in equilibrium by classifying all the terms allowed by symmetries. We then reduce this stress tensor along the time circle and compare it with that which follows from the partition function. This allows us to see a rich interplay between the surface transport coefficients and the bulk second order coefficients. Under the assumption of time independence, at first order on the surface, the non-zero linearly independent terms have been classified in tables 1, 2 and 3. The presence of the vectors u µ and n µ breaks down the local Lorentz symmetry at the surface to a smaller subgroup. The classification is based on transformation properties of the surface quantities under this preserved subgroup. We refer to the objects as scalars, vectors and tensors, depending on their transformation properties under this subgroup. Note that we have defined K µν = ∇ (µ n ν) .
We would like to point out that in table 2, we have not included the term P µα u β ∇ β n α because upon reduction it evaluates to the same result as P µα n ν ω να . Hence, in the stationary equilibrium case under consideration, these two terms are not independent. Also, 1 P µα P νβ K αβ − 1 2 P µν P αβ K αβ P ik P jm K km − 1 2 g km K Table 3. Symmetric traceless tensor. Upon reduction, T i 0 and T 00 components of the tensor vanishes.
owing to the identity 6) . Also note that in table 3, P αβ K αβ is distinct from K, since u α u β K αβ is non-zero. In the stationary case, P αβ K αβ reduces to K, while K = K + n · ∂σ, as shown in table 1.
Since we would like to have, the velocity at the surface, to be equal to the bulk velocity evaluated at the surface, there is no freedom in choosing a frame at the surface once the bulk frame has been chosen, as discussed in section 1.1.2. In order to respect the continuity of fluid variables and to naturally impose the condition (1.12), we shall proceed with the frame choice as described in section 1.1.2. This frame choice only constrains the form of the bulk stress tensor and leaves the surface stress tensor unconstrained. Therefore, while constructing the surface stress tensor at first order, we have to write down all possible terms that are allowed by symmetry without imposing any restrictions. We have where {S (i) , V ν (i) , T µν } are specified in the second column of tables 1, 2 and 3, respectively. The corresponding surface transport coefficients are denoted by As we may already foresee, among these 4 × 4 + 2 × 7 + 1 = 31 transport coefficients, only 3 are independent. The rest are determined in terms one another or bulk second order transport coefficients. We will now work out these relations.
If we consider the reduction of (3.11) along the time direction, we obtain the following reduced stress tensor where {S red (a) , V i (a) , T ij } are specified in the third column of tables 1, 2 and 3, respectively.

JHEP06(2016)015
Now comparing (3.12) with (3.9) we get (3.13) Let us recall from [19] that P 1 , P 2 and P 3 may be expressed in terms of the bulk transport coefficients in the following way 32 Finally, eliminating the P i variables from (3.13), we can summarize the following relations involving surface first order coefficients and second order bulk transport coefficients (3.16) These relations (3.16) are one of the central results of the paper. Let us now highlight some of the most interesting aspects of these relations (3.16). The last three relations in (3.16) relate bulk transport coefficients to those in the boundary. This shows that the linear response to particular deformation of the surface is intimately related to some, otherwise unrelated, transport coefficient in the bulk. Particularly interesting is the fact that t and κ 1 are proportional to each other. This physically implies that the linear response to a longitudinal stretch of the surface is entirely determined by how the fluid reacts to a change in background curvature. 32 Here we pick up only three specific relations; they can be expressed in several other ways using the relations between bulk second order coefficients, as was obtained in [19]. The bulk second order coefficients in (3.15) appeared in the Landau-frame stress tensor in the following way Tµν = T κ1R µν + κ2 K µν + λ3 ω α µ ω αν + λ4 a µ a ν + Pµν (ζ2R + ζ3R αβ u α u β + ξ3ω 2 + ξ4a 2 ) . (3.14) For further details of the conventions, we refer the reader to [19]. In the orthogonal-Landau-frame as defined in section 1.1.2, which is the most suitable bulk frame for describing the fluid configurations with a surface, the stress tensor takes the form (A.5). Note that the bulk transport coefficients appearing in (3.15) have very similar physical meanings in both the frames. See appendix A for further details.

JHEP06(2016)015
Another noteworthy fact is that the parity odd term introduced in (3.3), is reminiscent of the possible parity odd term in 2 + 1 dimensional space-filling fluids discussed in [19,20]. It leads to two non-zero parity odd coefficients, namely s2 3 and v1 5 . The scalar S 2 is proportional to n · , while V µ 5 is non-zero only when the acceleration at the surface a µ has a component parallel to the surface (see tables 1 and 2). It is interesting to note that although, space-filling uncharged fluids do not have any parity odd stationary transport at next to leading order, such a transport may exist at the surface of a finite lump of the same fluid.
Some of these constraints can be anticipated from the structural aspects of the conservation equation (3.10) on an arbitrary surface stress tensor, as explained in appendix B. In appendix C, the remaining constraints are also obtained through an entropy current argument, particularly adapted to deal with the stationary transport coefficients.

Description in terms of original fluid variables
In this section we lift the partition function of stationary neutral fluids (3.4) to a fourdimensional covariant action. 33 This action assumes that existence of a spacetime Killing vector field k µ with modulus k = −G µν k µ k ν along which the fluid flows are aligned, i.e., u µ = k µ /k as in [13,20]. We also assume that the surface is characterized by the same bulk Killing vector field restricted to the surface such that k| f =0 = −H ab k a k b , where H ab is the induced metric on the surface. Generically, we may write the effective action as the sum of a bulk and surface parts, This effective action, for neutral fluids up to second order, as in the case of the partition function, is described in terms of six transport coefficients, 18) whereR is the Ricci scalar of the spacetime, ω 2 = ω µν ω µν and a 2 = a µ a µ . Here, the pressure, the surface tension and all transport coefficients are functions of the local fluid temperature T which is given in terms of the global temperature T 0 and the modulus k via the relation T 0 = kT .
The bulk part of this action has been written down in [20], and the coefficientsP i measure the response of the fluid to background curvature, vorticity and acceleration. The surface part of this action was analyzed in [13] in arbitrary spacetime dimensions. However, there, as explained in appendix B, since the bulk pressure vanished at the surface, the scalars a µ n µ and K were not independent. Here the response coefficientB 3 is the surface modulus of rigidity of the surface fluid [13] whileB 1 encodes the response to centrifugal acceleration on the surface. Furthermore, dimension-dependent scalars were not analyzed 33 Note that we are using the terminology action in a slightly different way than [34]. This is because in the presence of a surface we can view (3.18) as an action functional for the surface f . Indeed, the surface part of (3.18), to leading order, is equivalent to the DBI action for co-dimension one branes when χ is constant and no worldvolume or background fields are present.

JHEP06(2016)015
in [13]. In this case, the scalar µ n µ is well known from the study of 2 + 1 parity odd fluids [20] and encodes the response due to vorticity at the surface. Despite being written in four spacetime dimensions, the action (3.18) generalises to arbitrary spacetime dimension withB 2 = 0.
The equations of motion can be derived from the action (3.18) by performing a general diffeomorphism of the form δG µν = 2∇ (µ ξ ν) and decomposing ξ µ into tangential and normal components to the surface such that ξ µ = e µ a ξ a + n µ n ν ξ ν as in [13]. The surface part of the variation of the action (3.18) yields where δ ξ denote variations along the co-vector field ξ µ and where we have defined This variation leads to two sets of equations of motion [13], which can equivalently be obtained from (3.10). One expresses conservation of the surface stress tensor in directions tangential to the surface, and is automatically satisfied for the stress tensor obtained for each contribution in (3.18). Indeed, the stress tensor (3.11) with the coefficients (3.13) satisfies (3.21). The other equation is the modified Young-Laplace equation, due to the presence of ∂ ρδ (f ) corrections in the surface stress tensor, T ab sur K ab = n ρ ∇ a ∇ b T acρ (2) + n ρ T abρ (2) n λ n µRλ aµb + T µν (0) n µ n ν , (3.22) where the stress tensor T ab sur , obtained by directly varying (3.18) with respect to the induced metric on the surface, is given in terms of the components (1.11) as 34 e a µ e a ν + n ρ T abρ whereR µνλρ is the Riemann tensor of the spacetime and where we have defined v c = u c ρ n λ ∇ λ n ρ . Note that if we do not consider higher order corrections, then T µνρ (2) vanishes and T ab (1) takes the perfect fluid form. In this case, equations (3.21) and (3.22) reduce to (2.15) and (2.17).
In order to understand the relation between these transport coefficients and the ones appearing in the partition function of section 3.3, we reduce (3.18) over the time circle and obtain Comparison with (3.4) leads to the identification of the pressure and surface tension as P = P and C = χ as well as to the relations between higher and lower dimensional transport coefficients. More precisely, we find and, if written in the form (3.5), then we can readily identify

Fluid configurations in 2+1 dimensions
In this section, we shall construct a few simple stationary fluid configurations to demonstrate the relevance of a non-zero surface entropy. The modifications of the fluid equations, in particular the Young-Laplace equation (2.20), when there is a non-zero surface entropy or equivalently a non-trivial dependence of surface tension on temperature, was discussed in section 2. Here, we will work out some particular solutions to these equations and explore the consequences of a non-zero surface entropy on the phase diagram of these fluid configurations.
We will keep our focus only on perfect fluids in 2+1 dimensions, ignoring possible higher derivative corrections. For working out explicit configurations we need the knowledge of the equation of state, for which we have to consider a specific system. For this purpose, we will consider the system of localized deconfined plasma of N = 4 Yang-Mills theory, compactified down to 2+1 on a Scherk-Schwarz circle, dual to rotating black holes and black rings in Scherk-Schwarz compactified AdS 5 , via the AdS/CFT correspondence [12,27].
Here, we shall revisit the analysis of [27] in order to find out how the results there are modified in the presence of non-zero surface entropy. In this section, we would like to briefly present our main results and therefore refer the reader to [27] and the references therein (also see [30,36]), for the background material.
One of the key assumptions in [27] was that the surface tension was constant, which we would like to relax here. This was assumed mainly because the value of surface tension for the interface between the confined and deconfined phase of N = 4 Yang-Mills was only known at the critical temperature [12]. The full dependence of surface tension on temperature for this interface is still unknown, but we would like to parameterize this ignorance in a suitable way and study its consequences.

Equation of state and thermodynamic quantities
The configurations that we shall deal with here (these are the ones that were already found in [27] in 2 + 1 dimensions), have the feature that the temperature is constant throughout the surface of the configuration. Also there is an empirical fact that the surface tension must decrease with the increase in temperature. This expectation follows from the fact that otherwise the surface entropy would become negative. Taking these two observations into consideration, and assuming that the surface temperature is very near to (slightly above) the phase transition temperature T c , we can consider the following dependence of the surface tension on temperature where χ 0 is the value of surface tension at T c . 35 We would like to emphasize that for the configurations of interest, the value of temperature at the fluid surface is in general different from the critical temperature. But, as we shall demonstrate later in this section, for all the configurations that we consider here, the surface temperature remains very close to T c , thus justifying our assumption. Also for the metastable plasma configurations to exist we must have that T s ≥ T c . As it will turn out, this condition will play a very important role in our analysis. In (4.1), κ is the parameter that parameterizes our ignorance about the exact temperature dependence of surface tension. Positivity of surface entropy implies that κ must be positive. The results in [27] were obtained with κ = 0 and the main goal in this section is to study how those results are modified as we turn on κ.
Using the thermodynamical relations from section 2, we find that the surface energy density χ E and the surface entropy χ S are respectively given by In addition, we also need the explicit form of the equation of state of the bulk fluid. Following [27], we take this to be where ρ is the energy density and ρ 0 is the shift in the free energy due to the Scherk-Schwarz compactification (see [27], for more details). The other thermodynamic quantities simply follow from (4.3) via the thermodynamic relations. In particular, the bulk entropy density and temperature can be expressed in terms of the energy density as follows (4.4) 35 In [12], it was found that at the critical temperature χ0 = 2.0 π 2 N 2 2 T 2 c for the plasma-balls in N = 4 SYM.

Spinning ball and ring
Before proceeding and state our results, we will briefly mention our conventions, while the rest of the details can be checked in [27]. The fluid configurations are in 2 + 1 dimensional flat space with metric ds 2 = −dt 2 + dr 2 + r 2 dφ 2 . (4.5) We seek rigidly rotating fluid configurations with the velocity vector u µ = γ{1, 0, ω} and γ = 1 − r 2 ω 2 1/2 . As an ansatz for rigidly rotating stationary configurations, the surfaces are are taken to be constant r slices in the spacetime (4.5). 36 Since the bulk fluid equations are not affected by κ, the solution in the bulk remains identical to [27]. The energy density in the bulk of the fluid has the form C being a constant of integration. At the inner and outer surfaces (denoted by the subscript − and + respectively), the Young-Laplace equation enforces (4.7) For the rotating balls, there is no inner surface and therefore no condition associated with it. The second term in (4.7) is precisely the acceleration term in the modified Young-Laplace equation (2.20). This additional term, in this boundary condition, is one of the two new modifications in our analysis compared to that in [27]. Now we proceed and obtain the phase diagram for the rotating balls and rings. We wish to plot the total entropy versus the total angular momentum at fixed total energy, for these configurations. The total energy E and the total angular momentum L is simply obtained by integrating the T tt and r 2 T tφ components of the stress tensor (2.10).
The total entropy S, is obtained by integrating the time component of the entropy current J µ s = s u µ . This is equivalent to integrating γs, where s is the total entropy density, including contributions from the surface where s is given by (4.4) while χ s is given by (4.2). This inclusion of the surface contribution to (4.8) is the second important modification in our analysis. Following [27], we introduce dimensionless quantities 36 The function f in the previous sections, defining the surface, is taken to be f = r+ − r for the outer surface and f = r − r− for the inner surface.

JHEP06(2016)015
where we have also defined a velocity v =rω. For the rotating ball we havẽ where v + is the velocity at the surface of the ball. In turn, for the rotating ring we havẽ where v + and v − are respectively the velocities at the outer and inner surface of the ring.
We must point out that v + , v − andω are not independent parameters for the rings. In fact, they must be related by the condition that the following two functions, must be identical (4.12) As expected, these expressions reduce to their counterparts in [27] when we set κ = 0.

Phase diagram for spinning balls and rings
The phase diagram that emerges out of (4.10) and (4.11) has been plotted in figure 2. We have plotted the total entropyS versus total angular momentumL for a fixed total energỹ E. We have used the same fixed value of energyẼ = 40, as in [27], so as to facilitate easy comparison. In fact, in both the plots in figure 2, we have displayed the κ = 0 phase diagram with light gray lines. We have displayed the phase diagram for two values of κ = 0.1, 0.5. The dark line represents the rotating plasma-ball solution while the blue and the green lines represent the rotating fat and thin plasma-ring respectively.
The main qualitative difference that we find here, as compared to [27], is that neither the spinning ball nor the rings reach zero entropy. This is because, at a given energyẼ, there is an upper bound for the velocity at the outer surface v + , which lies below 1, for non-zero κ. This upper bound on velocity is the point, where the curves terminate, while the curves for κ = 0 continue to zero entropy as v + approaches 1. This bound on v + arises from the fact that the temperature at the outer surface T + , reaches the phase transition temperature T c at the upper bound for v + . At higher values of v + , even if it remains below 1, the temperature at the surface would drop below T c and the configuration would cease to exist. We have demonstrated this behaviour of the surface temperature T + , in figure 3.
In figure 3, the temperature at the outer surface of the rotating ball and the fat ring have been plotted as a function of the velocity at the outer surface v + at fixed energỹ E = 40. For the thin ring, the behaviour of temperature is identical to that of the fat ring. The various lines represent values of κ ranging from 0.1 to 0.9, where the darkest line corresponds to 0.9. As it is apparent from figure 3, the value of v + for which the temperature dips below the dotted blue line, representing the phase transition temperature, decreases with the increase in κ.
For the rings, there is also a lower bound on v + , below which the solutions ceases to exist. This was also present for κ = 0. Also, as it is apparent from figure 3, the surface temperature for all the configurations remains very close to T c . This justifies our initial assumption that, in this analysis, we have taken the value of the surface tension and surface entropy evaluated at T c in (4.2).
The important consequence of this qualitative difference is that, for sufficiently large values of κ the phase transition between the ball and the ring configurations may disappear. As we can see from figure 2, such a phase transition does not exist for κ = 0.5. The critical value of κ at which this phase transition ceases to exist is approximately 0.33 ≈ 1/3. Thus, we see that the temperature dependence of the surface tension can crucially affect the existence of the phase transitions between fluid configurations. In the dual gravity picture, this would have important consequence for the phase transition between black holes of different horizon topologies. This calls for a future investigation, along the lines of [12] from the gravity side to ascertain the value of κ.
Finally, we would like to observe that the parameters determining the validity of our analysis are the same as in [27]. The first of such parameters is given by the change in the fluid temperature over the scale of the mean free path ∆u ≡ωv + /(1 − v 2 + ), which must be small for the fluid approximation to be valid. The other parameter is v + /ω for the ball and {v − , v + − v − }/ω for the rings, which must be large for δ(f ) to be well approximated by the Dirac delta. These parameters have a κ dependence throughω for a fixed value of energy. Both these parameters are not significantly affected by the value of κ (for the values of it that we have used), in the range of parameters that we consider. Therefore, we expect the validity of our result to be as good as that in [27].

Discussions
In this paper, after performing a systematic analysis of the nature of surface transport in relativistic fluids, we were able to significantly constrain the structural form of the fluid equations at the surface. We have focused on some particular cases during our analysis, namely perfect fluids in arbitrary dimensions and the next to leading order corrections to 3+1 dimensional relativistic normal fluids. Although we have a specific set up at the back of our minds, as indicated in section 1, our construction may be useful in more general settings, like when boundaries between different fluid phases are present. Since we do not use any particular form for the distributions θ(f ) and δ(f ), they can be suitably chosen to model a wide variety of situations. In order to serve a more general purpose, it may be particularly useful to study the nonrelativistic surface effects. This may be achieved by taking a non-relativistic limit of our set up following [37][38][39].
There are some immediate extensions of our work that are worth investigating. For instance, it would be very interesting to work out the next to leading order surface effects in superfluids. Due to the interplay between the vector which is normal to the surface, and the superfluid velocity, there may be a very rich, but yet unexplored surface transport properties in this case. In fact, while analyzing the zeroth order superfluids in section 2.2, we noticed a new term, in the modified Young-Laplace equation (2.27), which has not been widely considered in the literature.
It would be very interesting to understand the implication of this term on the thermodynamics of finite lumps of superfluids. This may be accomplished by undertaking an analysis of various possible superfluid configurations along the lines discussed in section 4. In fact, such an analysis may also provide direct hints towards the existence of hairy JHEP06(2016)015 black-rings or hairy black holes with other exotic horizon topologies in Scherk-Schwarz compactified AdS spacetimes via the AdS/CFT correspondence.
In section 4, we have analyzed the effects of temperature dependence of surface tension on phase diagram of some simple fluid configurations. We found that this effect can be very significant, especially while drawing conclusions about the existence of a phase transition between the ball and ring type configurations. Since, we have considered our sample system to be the same as in [27], our observation may have direct relevance for the existence of phase transition between spinning black holes and black-rings in Scherk-Schwarz compactified AdS 5 . In particular, if the surface tension for these configurations scales as the inverse of temperature, then our analysis suggests that the existence of such a phase transition cannot be reliably predicted by a fluid dynamical analysis. This observation calls for generalization of [12], to deduce the exact dependence of surface tension on temperature.
Another interesting extension of our work is the possible generalization to embedded fluids with surfaces of higher co-dimension and its application to the description of asymptotically flat and AdS black holes. As it is well known, both Myers-Perry black holes and the higher-dimensional Kerr-AdS black holes admit ultraspinning regimes [40,41]. Moreover, it was shown in [41,42] that these regimes can be described by a rotating fluid disc with a boundary, where the fluid is moving at the speed of light. These analytic solutions, therefore, could allow us to extract some of the new surface transport coefficients that we have found in this work and hence study the physical and stability properties of these black holes using the description of fluid dynamics with surfaces.
We would like point out one curious feature related to anomaly induced transport properties (see [43][44][45] for the most recent discussions on this). In our constructions here, we have treated the bulk of the fluid by multiplying the partition function of space filling fluids with a θ(f ) function. This procedure was justified (see section 1) by noting that θ(f ) denoted the change in the bulk transport coefficients at or near the surface. Since the usual transport coefficients are macroscopic parameters, representing the microscopic UV theory only in an effective way, this procedure of introducing the θ(f ) function is perfectly well-defined. However, there may exist certain terms in the partition function whose coefficients must be a constant as a consequence of gauge invariance [19]. The terms representing transport due to anomalies (for instance the term W anom in [19]), also falls within a similar category, since their form is fixed by the criterion that they reproduce the right anomaly coefficient locally everywhere in spacetime, including at the surface of the fluid. Such an anomaly coefficient is not an effective macroscopic parameter but a parameter of the microscopic theory. Therefore, those terms cannot be straightforwardly handled in an effective way by multiplying the, with θ(f ). 37 We postpone further analysis of these terms, in the context of fluid surfaces, to future work.
It would be interesting to see how the transport coefficients discussed in this paper fit into the classification of [5,46]. Further, recently there has been significant progress in formulating dissipative fluid dynamics in terms of an action [47,48]. It would be very

JHEP06(2016)015
interesting to understand how the presence of surfaces generalizes these constructions. In fact, it would be particularly interesting to understand the time evolution of fluid surfaces, involving dissipation. If we are able to incorporate time dependence, in a controlled fashion within our set up, it may have some relevance to situations concerning dynamical formation of surfaces, in the interface of two phases, which are described by the Cahn-Hilliard equations [49].

A Frame transformation in the bulk
In this appendix, we shall perform a frame transformation from the Landau-frame in the bulk of the fluid, to the orthogonal-Landau-frame which was defined in section 1.1.2.
In the presence of the surface at f (x) = 0, we can choose our coordinates so that one of the spatial coordinates vanishes at the surface. Let us refer to such a coordinate by f . For sufficiently well behaved spacetimes, the constant f surfaces would foliate the entire spacetime, including the bulk of the fluid. Every point on the constant f surfaces would admit a well defined, outward pointing normal vector, which we refer to as n µ . This provides us with an extension of the normal vector on the surface throughout the spacetime. 38 As discussed in section 1.1.2, instead of imposing the Landau-frame condition in the bulk, we make a slightly different frame choice, which is given by to all orders, everywhere in the bulk of the fluid. Here H µν = G µν − n µ n ν , is the projector orthogonal n µ , which is defined throughout the bulk of the fluid. It possible to impose this

JHEP06(2016)015
condition everywhere in the bulk, since we now have a definition of n µ extended throughout the bulk of the fluid. This immediately ensures that the fluid velocity is orthogonal to the normal vector at the surface of the fluid, where n µ is unambiguously defined. This frame transformation may be achieved by simply redefining Now, at the leading order, for stationary configurations, the fluid velocity can be oriented along the time-like killing vector, preserving the Landau-frame condition. Since, the surface has a trivial time evolution for stationarity configurations, this immediately implies n · u (0) must be zero. This is no longer true at higher orders and we need to perform a frame transformation by higher derivative terms in order to achieve (A.2).
In the partition function construction presented in [19], the fluid velocity and temperature in the bulk, are solved in terms of the background fields. The first non-trivial corrections to the fluid velocity occurs at the second order in derivatives. This implies that we have to perform a frame transformations with second order terms, and the transformation should have a form like (A.3). In fact, the exact form of the required transformation can be read off from the second order corrections to the fluid velocity in [19] 39 where P µν = G µν + u µ u ν , is the projector orthogonal to the fluid velocity and the ellipsis denote the higher order order corrections that may be necessary to keep (A.2) intact. We should take the coefficientsṽ 1 andṽ 2 to be the same as the ones appearing in the second order velocity corrections, worked out in [19]. This frame transformations directly impacts the form of the second order stress tensor. In the Landau-frame the second order bulk stress tensor was given by (3.14), and is now modified to T µν = T κ 1R µν + κ 2 K µν + λ 3 ω α µ ω αν + λ 4 a µ a ν + P µν (ζ 2R + ζ 3Rαβ u α u β + ξ 3 ω 2 + ξ 4 a 2 ) + ũ (µ n ν) v 1 a α ω αρ n ρ + v 2 P ρ α ∇ β ω βα n ρ .

(A.5)
Here v m = 2(E + P )ṽ m . All the second order fluid quantities in (A.5) are to be expressed in terms of the transformed velocityũ µ , and the vector field n µ is the extension of the normal vector at the surface throughout the bulk, as explained above. Now, we know from the partition function analysis that there are only 3 independent bulk transport coefficients and therefore v 1 and v 2 must be related to the rest of the transport coefficients, through two new relations. The necessary frame transformation to ensure (A.2) fixes them to be (A.6) 39 Note that the combination of the second order terms are chosen such that in the stationary situation, when expressed in terms of the background data, it reduces to the velocity corrections obtained in [19]. This choice may not be unique, specially when applying this trick to arbitrary orders, but as long as we focus only on stationary configurations, all such frames would be equivalent.

JHEP06(2016)015
These relations must hold in addition to the five relation between the rest of the transport coefficients in (A.5) as explicated in [19]. The fluid velocityũ µ that will now be obtained in terms of the background data, when we compare (A.5) with the bulk partition function in (1.18), will be automatically projected orthogonal to n µ , which has been ensured due to the frame choice (A.2).

B General constraints on the stress tensor
In this appendix we discuss generic constrains and symmetries of the surface stress tensor (3.11). The full spacetime stress tensor, including the bulk contribution, to second order in derivatives, can be decomposed as in (1.11), where the surface stress tensor (3.11) is written in the form Here, the structure T µν (1) denotes the contribution to the surface stress tensor of a monopole source of stress while T µνρ (2) denotes the contribution of a dipole source of stress. When applying this decomposition to (3.11) we easily read off T µνρ (2) = s1 4 P µν + s2 4 u µ u ν + s3 4 n µ n ν + s4 4 u (µ n ν) n ρ + v1 7 u (µ P ν)ρ + v2 7 n (µ P ν)ρ , (B.2) while T µν (1) includes all the other surface stress tensor components. However, we can impose additional constraints which follows from the fact that the the stress tensor (B.1) enjoys a symmetry, for which its components transform as (see [35] for more details) δT µν (1) = − µνa v a , δT µνρ (2) = µνa e ρ a , (B.3) for some coefficients µνa and where we recall that we have defined v a = e a ρ n µ ∇ µ n ρ . This transformation arises due to the freedom of introducing (D − 1) redundant delta functions in (B.1), so that (B.1) could have been written as . Therefore, the tangential derivatives of the distributionδ(f ) are integrated out and the coefficients T µνa (2) can be removed [35]. This implies that the terms involving v1 7 and v2 7 in (B.2) can be set to zero. 40 However, the stress tensor that follows from the partition function is obtained in a fixed gauge, as far as the transformations (B.3) are concerned. In that case, although we do get v1 7 = 0, v2 7 however, is related to other transport coefficients (see (3.16)).
Furthermore, for a stress tensor of the form (1.11), there is a perturbative symmetry that allows to displace the surface located f = 0 by a small amount ε such that f → f + ε. This symmetry expresses the freedom of defining the surface theory on a specific 40 In this paper we assume that the surface does not have boundaries. However, if there were boundaries then this symmetry would not present at the surface boundary and there we would need to impose µνa ηa = 0| ∂M for a normal covector ηa to the surface boundary.

JHEP06(2016)015
infinitely thin slice of a surface with finite thickness. Looking at figure 1, this means slightly displacing the dashed vertical line into another location within the distributionδ(f ). Under this infinitesimal displacement, the form of the stress tensor (1.11) is unchanged but its components have varied according to From (B.5) we see that this transformation induces a contribution proportional to n ρ ∂ ρδ (f ). If we take T µν (1) to have the perfect fluid form at leading order, then by appropriately choosing ε we could work with a surface for which either s1 4 or s2 4 vanish. However, since the transformation (B.5) induces a term proportional to T µν (0)δ (f ), such choice of surface would require to introduce the bulk pressure P | f =0 as an independent scalar in the surface part of the partition function. For this reason we have decided to work with the scalars K and a µ n µ instead. 41 There are structural consistency conditions on the components of the surface stress tensor (B.1), which, therefore, are not all independent and thus cannot be freely chosen. These consistency conditions arise due to the fact that we are working with the expansion (B.1) to a particular order and they can be derived by carefully analysing the conservation equation (3.10) to this particular order. The resulting conditions must hold in any physical situation, including time-dependent settings. One of these conditions constrains the dipole source of stress such that [35] n µ n ν n ρ T µνρ (2) = 0 , (B.6) and also, for codimension-1 surfaces, it follows that we must have e a µ n ν n ρ T µνρ which is a trivial consequence of there being no transverse two-plane on which the surface can rotate [13]. In turn, both conditions imply the constraints The remaining conditions determine the normal components of the monopole source T µν (1) in terms of the dipole source of stress T µνρ (2) . In particular we must have that If we now evaluate the normal components of T µν (1) in (3.11) we obtain n µ n ν T µν (1) = (s3 1 a ρ n ρ + s3 2 K + s3 3 ρ n ρ ) , (B.10) which upon comparison with (B.9) leads to the constraints When the bulk pressure vanishes at the surface, which is the case studied in [14], then K and a µ nµ are not independent.

JHEP06(2016)015
Moreover, the remaining normal components of T µν (0)surf must respect the following condition e a µ n ν T µν (B.12) Performing the same operation in (3.11) yields e a µ n ν T µν which, again, upon comparison with (B.12) leads to the constraints As it can be quickly verified, the constraints (B.8), (B.11) and (B.14) are a subset of the contraints captured by the partition function analysis (3.16). From these considerations, we see that the only components which are left unconstrained due to stress tensor conservation are the surface components of the stress tensor T ab (1) and the dipole source components T abρ (2) , also known as the bending moment. The remaining constraints in (3.16) in the case of relativistic fluids can be obtained by demanding positivity of the entropy current, an analysis which is carried in appendix C.

C Entropy current constraints
In this appendix we analyze the constraints on the transport coefficients that arise from the positivity of the entropy current and show that both the partition function and the effective action capture these constraints. This can be done by analyzing the divergence of the entropy current for a membrane subjected to external forces. The equations of motion were given in (3.21) and (3.22). The fact that the equations of motion only involve T ab sur , T abρ (2) and T µν (0) , signifies that only these three structures are required in order specific the dynamics of the membrane. However, we need to specify what the conservation equation for the surface entropy current is. The full entropy current can be expanded analogously to the stress tensor, where the surface part, up to first order, is given by For stationary configurations, we require (C.1) to be divergence free. In the bulk, this simply results in the bulk conservation equation ∇ µ J µ s(0) = 0 while in the surface this results in Here we have assumed that the entropy current can be obtained from the partition function/action in a similar way as a U(1) charge current, in the spirit of [5]. From the effective action (3.18), the surface entropy current can be obtained via the variation [14,15] J a s = ∂I surf ∂T u a , (C.5) though, depending on the type of corrections that the bulk action receives, there may be contributions to J a s due to bulk terms. Note that we have assumed that it is always possible to write the entropy current that follows from an action in the form (C.5), which does not include all terms allowed by symmetry. The reason for this is that for stationary configurations we are always free to add total derivative terms to the entropy current, which are also divergence free, such that it takes the form (C.5). This is true for uncharged fluids up to second order in derivatives, both for the bulk entropy current and for the surface entropy current. Since (C.3) depends only on J a s and J µ s(0) it is only necessary to classify these terms in order to obtain the constraints on the fluid transport.
Before proceeding and classifying possible terms that can appear in the different relevant structures, it is important to properly define fluid frames both in the bulk and in the boundary. Because, in principle, we can be placing a completely different fluid on the surface of another bulk fluid, we should consider two fluids, each described by their own fluid variables. In this appendix, we use tilde quantities to describe the bulk fluid, which is characterized by the set of bulk variables (T ,ũ µ ), while the surface fluid is characterized by the set of surface variables (T, u µ ). However, in order to fully specify the system, we need to impose boundary conditions on the bulk fluid variables. These boundary conditions were described in (1.14) and a natural consequence of them is that (C.6) These boundary conditions are dynamical, in the sense that the evolution of (T ,ũ µ ), described by the bulk equations to leading order whereΘ = ∇ µũ µ , must be subjected to the boundary conditions (1.14), which are dynamically determined by the surface evolution equations to leading order Note that (C.6) states that derivatives of the temperature along the fluid flows are equal in the bulk and in the surface. Derivatives of the temperature tangentially to the surface, but perpendicular to the fluid flows, are also guaranteed by (1.14) to be the same in both the bulk and in the surface. This can be seen by tangentially projecting the second equation in (C.7) and comparing it with the second equation in (C.8). In order to more clearly JHEP06(2016)015 present our results, we will impose (1.14) from the get go while carefully keeping track of the derivatives ofũ µ andT using equations (C.7). We now make a few comments regarding fluid frame transformations of the bulk and surface fluids. Frame transformations of the bulk fluid variables (T ,ũ µ ) allows us to set the bulk stress tensor in the orthogonal-Landau frame (A.2), that is, where T µν (0) are the higher derivative corrections to T µν (0) . However, due to the boundary condition (1.14), the restriction of such frame transformations to the surface at f = 0 will induce a frame transformation of the surface fluid variables (T, u µ ), if they are defined in the Landau frame. However, since we are working with stationary fluids up to second order in bulk derivatives, then such bulk frame transformations are second order in derivatives. Since the surface fluid quantities are only expanded to first order, then bulk frame transformations do not affect the surface stress tensor neither the surface entropy current. On the other hand, we can perform a first order frame transformation of the surface fluid variables (T, u µ ) and set the surface stress tensor in the Landau frame where again, T ‫)ג(‬ ab sur are the higher derivative corrections to T ab sur . This is not an elegant choice, in the sense that if we impose (C.10) for the surface stress tensor then we cannot simultaneously impose (C.9) for the bulk stress tensor evaluated at the surface. However, this is still a convenient choice because, as it will be explained below, the analysis of the divergence of the entropy current is insensitive to such frame transformations of the bulk stress tensor at the surface besides the fact that it reduces the number of structures appearing in the stress tensor to 4. We note, furthermore, that frame transformations do not modify the components T abρ (2) of the stress tensor [14]. We now proceed and write the relevant terms that enter the several structures involved. We note that, since we are working in the orthogonal-Landau frame (C.9)-(C.10) and since the entropy conservation equation (C.3) only involves the projection of eq. (3.21) along u µ , we find, using (A.5), that T µν (0) u ν n µ | f =0 = − 1 2 v 1ãαω αν n ν + v 2 P ν α ∇ βω βα n ν . (C.11) Moreover, the fact that only the contraction matters T µν (0) u ν n µ | f =0 for the divergence of the entropy current implies that surface fluid frame transformations do not affect it. On the other hand, the classification of J µ s(0) has already been done in [4]. As we will see later in this section, when comparing with the partition function and effective action, the contraction J µ s(0) n µ | f =0 for stationary configurations can be written as J µ s(0) n µ | f =0 = π 1 Θa µ n µ + π 2ω µb n µ a b + π 3 u a ∇ a (a µ n µ ) + π 4 u a a b K ab + π 5 ∇ a (ω µa n µ ) + π 6 u a ∇ a K + π 7 u a ∇ b K b a +s v 1ãαω αν n ν + v 2 P ν α ∇ βω βα n ν , (C.12)

JHEP06(2016)015
where the last line above is due to the frame change (A.4). Here the transport coefficients π i are only functions of T since we have restricted it to the surface. As the analysis of [4] shown, in a generic situation only 5 of the π i coefficients are independent but when stationarity is imposed, only 3 are independent [19]. Furthermore, the analysis of the remaining structures has largely been done in [14] but because of the presence of the bulk fluid and parity odd transport, we have T ab sur = χH ab +(χ+χ E )u a u b +α 1 KP ab +α 2 a µ n µ P ab +α 3 µ n µ P ab +α 4 P a c P b d K cd , (C.13) T abρ (2) = ϑ 1 H ab n ρ + ϑ 2 u a u b n ρ , (C.14) J a s = su a + γ 1 Ku a + γ 2 a µ n µ u a + γ 3 µ n µ u a + γ 4 u b K b a + γ 5ω µa n µ + γ 6 abc u b a c . (C. 15) It is worthwhile keeping in mind that all the 12 surface transport coefficients are only functions of T . 42 Comparing with the work in [14] for a free membrane, we note that the terms α 2 , α 3 , ϑ 2 , γ 3 , γ 5 , γ 6 were not present. As mentioned in section B, when the bulk pressure P vanishes at f = 0, then according to (2.20), the scalars K and a µ n µ are not independent and hence we need to include α 2 and γ 2 . Consequently, we must also include the coefficient ϑ 2 . Removing it would require, as discussed in section B, considering an extra term in T ab sur of the form P | f =0 P ab . The term γ 6ω µa n µ is a consequence of the presence of the bulk degrees of freedom. For stationary configurations, it may be replaced by a term of the form P a c u b K b c at the leading order; however, we must include it as a separate terms since they differ at higher orders. 43 Finally, the coefficients α 3 , γ 3 , γ 6 are well known in the context of parity odd fluids in 2+1 dimensions [26] and, as noticed in [26], γ 3 and γ 6 are not independent. Due to the freedom of adding to (C.15) a total derivative term of the form ∇ b α abc u c for an arbitraryα, shifting the coefficients γ 3 → γ 3 −α and γ 6 → γ 6 + T ∂ Tα −α, then only the linear combination γ 7 = T ∂γ 3 ∂T + γ 6 − γ 3 , (C. 16) is invariant under this shift. The surprisingly simple form of (C.13) can be obtained from (3.11) by bringing it to the Landau frame. Given these structures, we now impose the entropy conservation equation (C.3) and, using (3.21), obtain an expression of the following form ∇ a J a s − J µ s(0) n µ | f =0 = β 1 ΘK + β 2 Θa µ n µ + β 3 u b ∇ a K a b + β 4 u a a b K ab + β 5 u a ∇ a K + β 6 u a ∇ a (a µ n µ ) + β 7 σ ab K ab + β 8ω µb n µ a b + β 9 ∇ b ω µb n µ + β 10 Θ µ n µ . (C.17) We note here that the effect of the last line of (C.12) is cancelled by the force term (C.11). Close inspection of (C.17) leads us to conclude that all terms involved are linear in fluid 42 We could have added other terms to J a s such as terms proportional to a a and v a , however, these terms would be required to vanish at the end and hence, for simplicity, we have not considered them. 43 We do not consider terms which vanish for stationary configurations, upto the order we are keep tract of.
In principle, such terms should be considered as well but it is possible to show that they do not contribute to the analysis of the entropy current of stationary configurations.

JHEP06(2016)015
data. Since we must set the r.h.s. of (C.17) to zero, according to (C.3), then all β i coefficients must individually vanish, i.e., we must require that β i = 0 , ∀ i = 1, . . . , 10. We can solve these constraints in terms of ϑ i , γ i and the external force coefficients π i leading to α 4 = γ 4 T , γ 5 = π 5 , γ 2 T = π 3 + ϑ 2 , γ 1 T = π 6 T − ϑ 1 , γ 4 T = π 7 T + 2ϑ 1 , ∂γ 5 ∂T T = −π 2 , α 3 T −1 = −γ 7 T ∂T ∂s , (C. 18) In the case of no external forces π i = 0, no parity odd terms and the bulk pressure at the surface vanishing P | f =0 = 0, these constraints reduce to those in [14], while, instead, if we have parity odd terms and require no bending corrections ϑ 1 = ϑ 2 = 0, hence only α 3 , γ 3 , γ 6 remain, then this reduces to the result of [20,26]. There are a total of 10 relations in (C.18) relating the 6 surface transport coefficients α 1 , α 2 , α 3 , α 4 , ϑ 1 , ϑ 2 to the 7 external coefficients π i and the 6 entropy current coefficients γ i . Of these relations, 3 of them recover relations between bulk transport coefficients which were already known from a bulk analysis [4]. Specifically, these are the second relation in the first line, the second relation in the second line and the first relation in the third line of (C.18). From the 7 remaining relations, 6 of them determine the transport coefficients α 1 , α 3 , α 4 in terms of the transport coefficients ϑ 1 , ϑ 2 and the external coefficients π i while the remaining relation, namely the second relation in the third line of (C.18), relates the transport coefficient α 3 with the linear combination γ 7 defined in (C.16) in terms of the entropy current coefficients γ 3 , γ 6 as observed in [20,26]. Therefore, all surface coefficients appearing in (C.13)-(C.15) are determined in terms of the 3 transport coefficients α 3 , ϑ 1 , ϑ 2 and 3 independent π i coefficients, as expected, since it is indeed the number of independent scalars in both the partition function (3.4) and the action (3.18). Furthermore, from the partition function analysis, a total of 28 constraints were obtained in (3.16). From the analysis of appendix B, we have obtained a total of 14 constraints from (B.8), (B.11) and (B.14). However, as explained in appendix B the coefficients v1 7 and v2 7 may be removed by the transformation (B.3) and constitutes 2 of the 28 constraints in (3.16). The fact that (3.11) is not in the Landau gauge adds 9 extra constraints. Therefore, we have that 28-14-2-9=3, which is exactly the number of constraints that we have obtained from an entropy current analysis.

C.1 Comparison with the action and the partition function
We now compare the constraints obtained in (C.18) with the results obtained from the action (3.18) and the partition function (3.4). We begin by comparing (C.12) with the general form of the entropy current introduced in [4]. Using the notation of [19], the entropy current up to second order, and ignoring the first order corrections which vanish in equilibrium, can be written as J µ s(0) =sũ µ + ∇ ν A 1 ũ µ ∇ νT −ũ ν ∇ µT + ∇ ν A 2Tω µν + A 3 R µν − 1 2 g µνR ũ ν

(D.6)
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.