Fluid/gravity correspondence and the CFM black brane solutions

We consider the lower bound for the shear viscosity-to-entropy density ratio, obtained from the fluid/gravity correspondence, in order to constrain the post-Newtonian parameter of brane-world metrics. In particular, we analyse the Casadio-Fabbri-Mazzacurati (CFM) effective solutions for the gravity side of the correspondence and argue that including higher order terms in the hydrodynamic expansion can lead to a full agreement with the experimental bounds, for the Eddington-Robertson-Schiff post-Newtonian parameter in the CFM metrics. This lends further support to the physical relevance of the viscosity-to-entropy ratio lower bound and fluid/gravity correspondence. Hence we show that CFM black branes are, effectively, Schwarzschild black branes.


I. INTRODUCTION
Black hole solutions of the Einstein equations in general relativity (GR) are useful tools to investigate the space-time structure and quantum effects in any possible theory of gravity. In particular, models with extra spatial dimensions are prominent candidates as extensions of Einstein GR [1], leading to important consequences, not only for black hole physics [2][3][4], but also for particle physics, cosmology, as well as for the astrophysics of supermassive objects [5]. In this context, black strings play a notable role in the quest for physically viable and realistic models [6][7][8].
Fluid dynamics is an effective description of interacting quantum field theories (QFTs) at long wavelengths [9]. As a low energy effective field theory, fluid dynamics is determined by a derivative expansion of the local fluid variables, and is used to describe near equilibrium systems. The derivative expansion is determined by the symmetries of the system and its thermodynamical features.
Transport coefficients, such as viscosities and conductivities, assess how perturbations propagate under relaxation and, besides being theoretically predictable, they can be measured experimentally.
The AdS/CFT correspondence, which relates the coupled dynamics of SU (N ) gauge theories to gravity in anti-de Sitter (AdS) spaces [10][11][12][13], has been used to connect those two realms of modern physics, namely gravity and QFT. It holographically relates a strongly interacting QFT at large N with the dynamics provided by gravity in asymptotically AdS space-time. In the long wavelength limit, the fluid/gravity correspondence has been recently proposed as a framework related to the AdS/CFT correspondence [14][15][16], mapping black holes in asymptotically AdS space-times to the fluid dynamics of a strongly coupled boundary field theory. One of the most remarkable predictions of AdS/CFT and fluid/gravity correspondence is the ratio of the shear viscosity-to-entropy density, which is universal for a large class of strongly coupled (gauge theory) isotropic plasmas [9]. The universality of this ratio plays a prominent role for gauge theories that are dual to gravitational backgrounds [14][15][16][17][18][19][20]. Moreover, the same results were obtained by employing the Green-Kubo formula in Ref. [21].
On the other hand, black strings present hydrodynamic features such as viscosity, diffusion rates, diffusion constants and other transport coefficients [22], besides temperature and entropy.
From the perspective of the holographic principle [10], a black string corresponds to a certain finitetemperature QFT in lower dimensions. In this context, the hydrodynamic features of the horizon of a black string can be identified with the hydrodynamic behaviour of a dual theory. For such QFTs, the ratio of the shear viscosity η to entropy density s is bounded by the universal value [17] η s 4 π k B 6.08 × 10 −13 K s , for an extended range of thermal QFTs. The above inequality holds for all known substances, including, for example, liquid helium, water, and the quark-gluon plasma, produced at the Relativistic Heavy Ion Collider (RHIC). Therefore, it has been conjectured to represent a universal lower bound for all materials 1 , which is called the Kovtun-Starinets-Son (KSS) bound [16,17]. In the case of Einstein gravity, the ratio η/s equals 1/4π. This bound is obtained in bulk supergravity from a recent calculation in N = 4 supersymmetric SU (N c ) Yang-Mills theories, in the regime of infinite N c and large 't Hooft where ζ(3) is the Apéry constant, and the next-to-leading term represents the first string theory correction to GR. The KSS bound is a powerful tool for studying strongly interacting systems, like the quark-gluon plasma, and is further employed to analyse trapped atomic gases. Strongly interacting Fermi gases of atoms have been noticed to be ruled by hydrodynamics [23], presenting finite shear viscosity at finite temperature. The bound (2) holds in relativistic QFTs at finite temperature and chemical potential [16,17], for at least a single component non-relativistic gas of particles with either spin zero or spin 1/2 [24]. The entropy of the dual QFT equals the entropy of a black string, which is proportional to the area A (volume) of its event horizon in Planck units, According to the AdS/CFT duality, any given asymptotically AdS bulk geometry should be equivalent to specific states of the gauge theory on the boundary. In fact, the AdS bulk geometry can be mapped into vacuum states of the gauge theory, whereas AdS-Schwarzschild black holes are related to thermal states of the corresponding gauge theory on the boundary. Hence, a fluid can be associated with the bulk black hole solution. Moreover, the bulk dynamics is governed by 1 We shall mostly use units with = c = kB = 1, four-dimensional indices will be represented by Greek letters, µ = 0, . . . , 3 and five-dimensional indices by capital latin letters, A = 0, . . . , 4.
solution is equivalent to a thermal state, corresponding to black strings in AdS [26]. The dynamical framework can be realised by a system in local thermodynamical equilibrium, corresponding to a non-uniform black string in AdS, which evolves according to effective hydrodynamical equations corresponding to the Einstein's equations [25]. Such a system will subsequently relax towards a global equilibrium state, corresponding to a uniform black string in AdS [26].
Motivated by the above picture, we shall here explore the fluid/gravity duality for the fivedimensional black string [5,27,28] associated with the effective four-dimensional Casadio-Fabbri-Mazzacurati (CFM) solutions [29]. CFM metrics generalise the Schwarzschild solution and have a parametrised post-Newtonian (PPN) form with parameter β, which is the usual Eddington-Robertson-Schiff parameter used to describe the classical tests of GR 2 . This parameter allows for a direct comparison with experimental data in the weak-field limit, which is accurate enough to include most solar system tests [30]. In general, one can construct a dual fluid stress tensor by solving the bulk Einstein equations in asymptotically AdS spacetime and taking the long wavelength limit, but there are only a few black hole solutions that admit such a hydrodynamic dual on the boundary, even asymptotically. The prototype is the Schwarzschild metric, however the CFM solutions can also be obtained as four-dimensional vacuum brane solutions of Einstein gravity with a cosmological constant. In a second-order derivative theory of five-dimensional gravity interacting with other fields in AdS, for instance type II B supergravity in AdS 5 × S 5 , the CFM solutions can be realised as universal sub-sectors, regarding pure gravity with negative cosmological constant, in the long-wavelength limit.
On the gravity side, we shall thus consider a five-dimensional CFM black string, associated with a four-dimensional CFM black hole, and the four-dimensional dual fluid description will necessarily carry an imprint of the parameter β. The fluid/gravity correspondence will then be shown to provide a bound for β, corresponding to the universal KSS bound (2). In particular, both the shear viscosity and black string entropy will depend on β, and Eq. (2) will read which is precisely what gives rise to the bound on β. We shall also argue that this bound can be made to agree with observations, provided higher order terms are included in the hydrodynamic expansion of the stress tensor.
This paper is organised as follows: in Section II the CFM solutions are reviewed; in Section III, the Green-Kubo formula is used to bound the PPN parameter through the shear viscosity-toentropy density ratio in the linear regime, and we show how the estimates improve by including higher-order hydrodynamical terms in Section IV; finally we discuss our results in Section V.

II. CASADIO-FABBRI-MAZZACURATI BRANE-WORLD SOLUTIONS
Solutions of the Einstein field equations on the brane are not uniquely determined by the matter energy density and pressure, since gravity can propagate into the bulk and generates a Weyl term onto the brane itself. Taking that into account, the CFM metrics [27] are vacuum brane solutions, like the tidally charge metric in Ref. [31], and contain a PPN parameter β, which can be measured on the brane. The case β = 1 corresponds to the exact Schwarzschild solution on the brane, and extends into a homogeneous black string in the bulk. Furthermore, it was observed in Ref. [29] that β ≈ 1, in solar system measurements. More precisely, the deflection of light in the classical tests of GR provides the bound |β − 1| 0.003 [30]. The parameter β also measures the difference between the inertial mass and the gravitational mass of a test body, besides affecting the perihelion shift and describing the Nordtvedt effect [29]. Finally, measuring β provides information regarding the vacuum energy of the brane-world or, equivalently, the cosmological constant [1,29].
We recall the effective four-dimensional Einstein equations on the brane can be expressed as where κ 2 5 = 8 π G 5 and G N = κ 2 5 σ/48 π (σ being the brane tension), T µν denotes the energymomentum tensor of brane matter, T ≡ T µ µ , and E µν is the Weyl tensor term. For brane vacuum, T µν = 0, and absorbing the bulk cosmological constant Λ 4 into the bulk warp factor, the above field equations reduce to We are in particular interested in static spherically symmetric systems on the brane, whose general metric can be written as The CFM metrics are obtained by relaxing the condition N (r) = B −1 (r), valid, for example, for the Schwarzschild and Reissner-Nordström metrics. In fact, for black strings, such a condition will result in a central singularity extending all along the extra dimension and a singular bulk horizon [29], a configuration which moreover suffers of the well-known Gregory-Laflamme instability [32]. The CFM solution I is obtained, by fixing N (r) equal to the Schwarzschild form and then determining B(r) from the field equations (6), whereas the CFM solution II follows from the same procedure but starting from a metric coefficient N (r) of the Reissner-Nordström form [31].

A. CFM solution I
For the first case, the metric coefficients in Eq. (7) are given by The solution (8) depends upon just one parameter β and the Minkowski vacuum is recovered for M → 0. The horizon radius r = R on the brane is then determined by the algebraic equation 1/B(R) = 0, and this black hole is either hotter or colder than the Schwarzschild black hole of equal mass M depending upon the sign of (β − 1) [29]. For example, assuming β = 5/4, one finds two solutions equal to the Schwarzschild radius R S = 2 G N M . Note also that the four-dimensional Kretschmann scalar K (I) = R µνρσ R µνρσ diverges for r = 0 and r = 3 G N M /2 < R S .

B. CFM solution II
The second solution for the metric coefficients reads In this case the classical radius R for the black hole horizon is given by R = R S and R = R S (β − 1/4). Moreover the Kretschmann scalar diverges at each black hole horizon [5], hence indicating a physical singularity.
A comprehensive analysis of the causal structure and further features on both CFM solutions can be found in Refs. [5,27,29].

III. BOUNDING THE PPN PARAMETER IN THE LINEAR REGIME
As previously mentioned, hydrodynamics can be viewed as a low-energy effective description for the low-momentum regime of correlation functions. We are particularly interested in the Green-Kubo formula in the linear response theory, where an important fluid transport coefficient, the shear viscosity, arises in response to a perturbation in the fluid stress-energy tensor. The Kubo formula can further relate the shear viscosity to the absorption cross section of low-energy gravitons [19] 3 .
Let us consider an action S, and introduce a set of sources J a coupled to a set of operators O a , where Using the interaction picture in QFT [13], one can immediately see that the variation where G R (q) denotes the retarded Green function in momentum space with q µ = (ω, k). This simple linear formalism can be generalised to curved space-times and non-vanishing VEV in order to describe a fluid in the conditions of interest here.
A common way to derive the Kubo formula is, in fact, to couple four-dimensional gravity g µν to the fluid and determine how the fluid energy-momentum tensor responds to gravitational perturbations 4 . In particular, for small metric perturbations which vary slowly in space and time, g µν =ḡ µν +h µν with h µν 1, the stress tensor can be doubly expanded in the metric fluctuation h µν and in gradients [33]. This procedure yields a generalisation of Eq. (10) to higher orders for the operator T µν , namely where G µν|... R are retarded n-point correlators, with the measurement point z having the largest time, and we also assume the metric fluctuations h ρσ vanish in the far past. The fluid response is then obtained from the stress tensor conservation law ∇ µ T µν = 0, together with the condition that the fluid describes a conformal theory, namely, T µ µ = 0. We shall here analyse the first-order formalism, and illustrate how our results can be further refined by considering second-order terms in the next section.
At zero order in derivatives, the energy momentum tensor reads where u µ is the fluid four-velocity, denotes its energy density, P the pressure andḡ µν is again the unperturbed metric on the four-dimensional boundary where the fluid lives 5 .
Since in the linear response theory the metric can be seen as the source of the stress-energy tensor, the response of the two-point function of the associated perturbed metric can then be written as the first order term T µν in the expansion (13) [34], where the retarded Green function is G µν|αβ R = T µν (x) T αβ (y) . In the gauge/gravity duality, the metric is dual to the stress-energy tensor [12]. Hence, from the gravity point of view of the correspondence, the above perturbation should arise from metric fluctuations h µν of the appropriate black brane metric.
The first-order expansion (in derivatives) of the hydrodynamical stress-energy tensor, also known as the constitutive equation, includes dissipative terms, as well as shear and bulk viscosities. In curved space-times, the first order contribution is given by [35] T µν where η is the shear viscosity, ζ the bulk viscosity, ∇ µ represents the covariant derivative in the (generally curved) space-time metricḡ µν . In addition, P µν =ḡ µν + u µ u ν is the projection tensor along spatial directions, that enables us to express the constitutive equation in a covariant manner.
We shall next restrict ourselves to a particular type of perturbations which is simpler to treat using the AdS/CFT correspondence. In fact, we suppose, as usual, that the contribution to the shear viscosity η comes only from the component h 12 = h 12 (t) of h µν , with all other h µν = 0 [13,34].
Since we assumed such fluctuations around thermal equilibrium are small, we can say the fluid has uniform temperature T (x µ ) = T 0 and is at rest in the chosen frame, that is u µ = (1, u i = 0).
Moreover, this perturbation has spin-2 with respect to the spatial SO(3) group and it cannot excite linear order fluctuations of the velocity (which is a vector) or of the temperature (a scalar).
Therefore, u i = 0 and T = T 0 remain valid up to linear order.
The shear viscosity of the dual theory can be computed from gravity in a number of equivalent approaches [13,14,16]. We shall employ the Kubo formula, which relates the viscosity to equilibrium correlation functions, in order to compute the T 12 component of the stress-energy tensor, in a large distance and long time scale regime. The non-vanishing contribution in the covariant derivative arises from the Christoffel symbol The component ∇ 2 u 1 is obtained analogously and the other components vanish. Hence, only the first two terms inside the brackets in Eq. (16) contribute to T µν (1) , and we have By taking the Fourier transform of Eq. (18), one obtains [34] δ T (1)12 (ω, k = 0) = i ω η h 12 .
A perturbed fluid Lagrangian is correspondingly given by for which Eq. (12) reads where (19) represents the same linear response relation as Eq. (20), and the Green-Kubo formula can therefore be obtained as [13,34] where denotes the imaginary part.
Our first concern here is to provide analytic expressions for asymptotically AdS black strings which solve the full five-dimensional Einstein equations. The CFM black strings are "tubular", codimension one, black branes in AdS [5], corresponding to a late time behaviour of their dual fluid. Two families of such solutions are obtained from the same assumption that leads to the CFM brane metrics reviewed in Section II, namely by relaxing the condition that the metric coefficient g tt is (minus) the inverse of the metric coefficient g rr . These two families of black strings are likewise expressed in terms of the PPN parameter β [29], and reduce to the five-dimensional Schwarzschild black string solution, when β → 1. CFM black branes in AdS can be also obtained as a perturbative approach of CFM black holes on the brane [5]. Our five-dimensional metrics are, in particular, given by where , the length is related to the AdS curvature [26,36] and the metric coefficients N and B equal the CFM expressions (8) or (9). The entropy of these black strings can be computed by using the Hawking-Bekenstein formula [16], where R is the horizon radius determined by 1/B(R) = 0 and V 3 = dΩ. Moreover, the volume density of entropy s = S/V 3 must be finite, which requires that β < 5/4 in Eqs. (8) and (9). It is also convenient to trade the radial coordinate r for the dimensionless u = R/r, and define the lengthM = G N M , so that and the metric finally reads Our main aim is now to predict a value for β exclusively from the KSS bound (2), implied by the fluid/gravity correspondence. The background metric (25) can be perturbed as g AB → g AB +h AB [14,17]. As mentioned before, the contribution to the shear viscosity comes from a single component of the metric fluctuations. Denoting this term by φ = φ(t, u, x), the corresponding wave equation reads [16,34] By taking the Fourier transform along the time direction and disregarding the spatial momentum, φ e i ω t Φ(u), one finds which can be specialised for both CFM metrics I and II, respectively provided by Eqs. (8) and (9), that is where V = V I,II , with for the CFM case I and for the CFM case II. Following the gauge/gravity prescription, we can then proceed to calculate the retarded Green function [14][15][16], for both cases in Eqs. (29) and (30), and then derive the shear viscosity from the Green-Kubo formula (21). A bound on the PPN parameter β can finally be obtained from the KSS lower bound (2). From the propagator (31), by inserting Eq. (21) in Eq. (4), one in fact finds where we remarked that both s and η do depend on β. However, the solutions to Eq. (31) are cumbersome linear combinations of integral exponentials and Hermite functions, and we shall therefore just display the results graphically, in the near horizon approximation.
The resulting bound for β is displayed graphically in Fig. 1 [30]. Both experimental bounds fall outside the bands shown in Fig. 1, although they are close in value. Nevertheless, we shall argue in the next Section how the situation can be significantly improved.

IV. SECOND ORDER IMPROVEMENTS
The results obtained in Section III are not far from the observational data, nevertheless they still appear to be ruled out by the experimental bounds. This could signal that the fluid/gravity correspondence fails for these systems, or that the approximations we employed are too strong. In particular, Eq. (32) was obtained in the linear hydrodynamics regime, corresponding to the first two terms in Eq. (13), that is Eqs. (14) and (16). We can instead tackle the non-linear dynamics by keeping higher-order terms in the expansion (13) of the stress-energy tensor in velocity and gradients multiplied by the corresponding hydrodynamic coefficients [25]. In particular, we shall here include the second-order regime of a dual hydrodynamic prescription [9,25,33,37,38], in order to refine the bounds on the PPN parameter β.
The second-order hydrodynamic analysis is rather involved [9,25,33,37,39,40], so let us first introduce some compact notation. For any second rank tensor O αβ , we denote O αβ ≡ 1 2 P αµ P βν O (µν) − 1 3 P αβ P µν O µν , which is transverse to the fluid motion, u µ O µν = 0, and traceless g µν O µν = 0. In addition, σ αβ = P αµ P βν ∇ (µ u ν) − 2 3 P αβ P µν ∇ µ u ν is the shear tensor and Ω αβ = P αµ P βν ∇ [µ u ν] denotes the vorticity tensor, as they appear in the T µν (1) of Eq. (16). We next express derivatives along the fluid 4-velocity in terms of transverse derivatives by means of the zero-order equations of motion [33], where ∇ µ ⊥ ≡ P µα ∇ α and ∇ ⊥ · u = ∇ · u ≡ ∇ µ u µ . All possible contributions to the stress-energy tensor are σ µν (∇·u), and R µν , where R αβγδ is the Riemann tensor. But the only ones that transform homogeneously under Weyl transformations are [33] A µν Eq. (34) can be recast as [33] u ρ ∇ ρ σ µν + 1 2 σ µν (∇·u), which by using (33), reduces to the linear combination A µν The second order dissipative part in the expansion of the stress-energy tensor (13) must therefore read where the κ term contributes to the two-point Green's function of the stress-energy tensor. The λ i terms are non-linear in velocity and the parameter τ Π can be thought of as the relaxation time [33].
It is worth mentioning that the second-order hydrodynamic transport coefficients κ and λ i have been calculated from the leading order AdS/CFT correlators at finite chemical potential and in the presence of fundamental matter [33,38,39]. They correspond to the black brane horizon response properties. Moreover, these coefficients have been derived for strongly coupled plasma described holographically by an Einstein-scalar action in the bulk, as functions of the temperature in lattice QCD thermodynamics [41]. Nonetheless, there is not a final word concerning the values of these parameters. For conformal theories at finite temperature, and zero chemical potential, they scale with suitable powers of the temperature and may depend on coupling constants and eventually on the rank of the gauge group [33,42]. In fact, there is an intrinsic ambiguity in second order viscosity parameters in relativistic hydrodynamics [42]. For example, one of the coefficients can be calculated as λ 3 = −4 lim k 1 →0 G xy|0x|0y . The expressions thus obtained are very similar to those given in Refs. [33,39], namely where γ = (g 2 N c ) −3/2 ζ(3)/8 and T again denotes the system temperature. We naturally choose the latter to equal the black string Hawking temperature, given by Eqs. (13) and (21) in Ref. [29] for the CFM solutions I and II respectively, and depending on M and β. The strong coupling limit γ → 0 can be also obtained [14,25]. We emphasise that our results below are independent of changes in the above transport coefficients, at least for the parameters derived in Refs. [14,25,42].
Similar steps to the ones yielding Eq. (17) and the Kubo formula (21) can be repeated. If one now takes, as usual in the literature, h 12 = h 12 (t, z), Eq. (36) yields where dots denote time derivatives. The retarded Green function then reads which is the analogue of the Kubo formula with second-order derivative terms taken into account.
We then see the kinetic parameters τ Π and κ are the coefficients of the ω 2 and k 2 terms in the low-momentum expansion of the retarded propagator.
Upon employing the above expressions and repeating the analysis of Sect. III, one finds new bounds for the PPN parameter β, and preliminary results are shown in Fig. 2. All terms in Eq. (36) are included, with transport coefficients provided by Eqs. (37a)-(37e) correspond to the second term for M = 10 5 M • . The numerical stability of these results has been checked by considering several different orders of magnitude for the t'Hooft coupling, and the differences were always found to be negligibly small. Moreover, refined expressions involving more transport coefficients, like in Ref. [41], could be used. However, such contributions in the viscosity-to-entropy ratio would amount at most to order 10 −5 corrections in the bounds for β. Similarly, we tested other possible choices for the second-order viscosity parameters (37a)-(37e) and the differences with respect to our results were always found to be negligibly small. A comment is due about the relaxation time τ Π in Eq. (37b), which turns out to affect the viscosity-to-entropy ratio by a negligible order of 10 −6 . This transport coefficient, which appears in the first term on the r. T µν (n) are taken into account in Eq. (13). This is in agreement with the conjecture of Ref. [51], that the Schwarzschild black brane is the unique static black hole solution of vacuum gravity with a negative cosmological constant.

V. DISCUSSION AND OUTLOOK
In this work, we derived the shear viscosity-to-entropy ratio η/s as a function of the PPN parameter β of candidate CFM black strings from the Green-Kubo formula in the context of the fluid/gravity correspondence. On assuming the KKS universal bound (2), we then determined the allowed range of β for different black string mass. Our results show a possible slight dependence of β on the mass in the range M • < M < 10 4 M • in the linear order hydrodynamic expansion (see Fig. 1). The non-linear hydrodynamical regime was just started to be investigated including second order terms. In particular, when the terms displayed in Eq. (36) are included, with the transport coefficients (37a)-(37e), the theoretical bounds are modified and appear in agreement the four-dimensional observational bounds (see Fig. 2). The discrepancy with experimental data shown in Fig. 1 could thus be simply related to the approximations employed in order to derive Eq. (32) from the Green-Kubo formula, for instance, the near horizon limit and the linearised hydrodynamics regime. Of course, this will be the subject of further investigations. We showed that the KSS result implies that the CFM black branes are, effectively, Schwarzschild black branes, when sufficiently high-order terms in the boundary energy-momentum tensor (13) are taken into account, in agreement with the conjectured unicity of the Schwarzschild black brane as a static black-hole solution of vacuum gravity with a negative cosmological constant [51].
As more possible developments, we recall that higher derivative gravity corrections have been considered [33,45]. The ratio of shear viscosity-to-entropy density indicates that the KSS bound (2) could then be violated, and the new lower bound η/s ≥ 4/25 π was proposed [21], based on causality of the dual field theory. In this context, it seems that such higher derivative gravity corrections make the CFM solutions inadequate as a black brane prescription. Indeed, from the bound η/s 4/25 π, the PPN parameter of these metrics would be β 0.7432, which contradicts the theoretical bound β ≥ 1. A correction for η/s was also given in Ref. [45] in the form of a linear combination of the coefficients α i of the scalars R, R µν R µν and R µνρσ R µνρσ in the Gauss-Bonnet contribution to the action. More precisely, in the Gauss-Bonnet context, η/s (1 − 8 α 3 )/4 π and the KSS bound is violated for α 3 > 0. However, since the CFM solutions are obtained with no Gauss-Bonnet terms, one can set α 3 = 0 and there is no anisotropy, so that the KSS bound (32) is the only one that does hold for CFM metrics. The quark-gluon plasma produced in the ultra-relativistic collision of heavy ions at RHIC and LHC [46] is a strongly coupled fluid that is represented by the anisotropic evolution of the fluid in the initial stages before isotropisation. Recently, there has been some interest in modelling this anisotropy at strong coupling [47,48] and in studying how observables may be affected by it. Some of the studies include the computation of the shear viscosity-to-entropy density ratio [49,50]. Anisotropic scenarios can be further considered to include such possibilities.