Gregory-Laflamme analysis of MGD black strings

The minimal geometric deformation (MGD), associated with the 4D Schwarzschild solution of the Einstein equations, is shown to be a solution of the pure 4D Ricci quadratic gravity theory, whose linear perturbations are then implemented by the Gregory-Laflamme eigentensors of the Lichnerowicz operator. The stability of MGD black strings is hence studied, through the correspondence between their Lichnerowicz eigenmodes and the ones associated with the 4D MGD solutions. Its is shown that there exists a critical mass driving the MGD black strings stability, above which the MGD black string is precluded from any Gregory-Laflamme instability. The general relativistic limit leads the MGD black string to be unstable, as expected, corresponding to the standard Gregory-Laflamme black string instability.


I. INTRODUCTION
Recent enterprises to establish theories of gravity, in a scope that is beyond the general relativity (GR), have been reasonably successful attempts to describe our universe. The prototypical Einstein-Hilbert action is well known to govern the GR Einstein equations. Adding terms to such an action, that are beyond the scalar curvature, does not immediately imply that the derived equations of motion are reduced to the GR Einstein equations, even at low energy regimes. A recently proposed setup, known as the minimal geometric deformation (MGD) [1][2][3][4], has been shown to be a successful framework procedure to generate realistic analytical solutions of the effective Einstein field equations, by deforming GR solutions [2,3]. The MGD procedure can provide new compact stellar distributions [4][5][6][7], with solid crusts due to a bulk Weyl fluid [8], that might collapse into black holes. This procedure is usually implemented by considering more intricate forms of the stress-energy tensor, that are beyond the ones describing perfect fluids in GR [6,9]. The MGD setup can be established for anisotropic spacetime fluids and other interesting examples already paved in the literature [1].
The key point of the MGD scenario was originally derived on brane-world models, wherein the 4D universe consists of the brane itself, that is a solution of an effective string action, placed into a 5D AdS bulk [10][11][12]. The brane has finite self-energy described as the brane tension σ, which can fluctuate. At low energy regimes, the brane tension emulates a chemical potential characterizing the energy expense for, at least hypothetically, creating any unit brane volume [13]. An infinite value of the brane tension dictates a perfectly rigid brane, corresponding to the GR limit. The brane tension can be approximated by a constant parameter related to the brane self-energy or an effective 4D brane cosmological constant, in an equilibrium brane configuration. However, more realistic scenarios, that further incorporate brane-world cosmology and includes the universe expansion and the reheating era as well, are described by variable tension branes [14][15][16]. Moreover, the brane tension has, still, zero mean local fluctuations, that are position dependent.
The MGD method is essentially based upon how 5D bulk constituents, encrypted in the 5D bulk Weyl tensor describing a Weyl fluid, can induce physical effects on the 4D brane. In fact, the 5D bulk Weyl fluid percolating the brane implements the MGD setup, responsible to deform the radial component of any static, spherically symmetric, metric, incorporating the inverse of the brane tension σ −1 . MGD solutions are led to the 4D Schwarzschild solution at low energy regimes, when σ → ∞.
One of our aims here is to provide an alternative description that does not make necessary the existence of a 5D bulk Weyl fluid, whose effects can be, for our purposes, emulated by a 4D modified Ricci quadratic gravity setup. In fact, the Sasaki-Shiromizu-Maeda procedure makes the 5D bulk pure gravity, described by the 5D Einstein equations with cosmological constant, to be projected onto the 4D brane [17]. The Gauss-Codazzi method then splits the 5D stress-energy tensor into a 4D energy-momentum tensor, that consists of terms describing, respectively, the energy and matter on the brane, Kaluza-Klein signatures from the bulk, and high energy corrections due to the 5D Weyl fluid [17,18]. The paradigm throughout the present paper relies on the perturbing the Schwarzschild-like solutions in 4D Ricci quadratic gravity [19], alternatively dismissing the necessity of a 5D Weyl fluid to describe the MGD solutions. Hence, this setup provides the study of MGD black holes and the stability of MGD black strings.
Standard black strings were shown to be unstable to long wavelength perturbations [20][21][22]. The first attempt to localise the black string close to the brane appeared in Refs. [23,24] . Gregory-Laflamme instabilities in the black string are related to the Rayleigh-Plateau instability, in fluid mechanics [25]. However, we shall show that there is a critical black string mass, above which the MGD black strings can be stable under long wavelength perturbations. Indeed, besides placing the MGD procedure, complementarily, as a legitimate 4D solution of the Ricci quadratic gravity theory, this paper is also devoted to study the stability of the MGD black strings under Gregory-Laflamme perturbations. In a previous work, the MGD black string was studied, showing that the brane tension robustly influences its event horizon [26]. The MGD scenario is suitable for the study of black strings, as brane 4D black holes can be produced by gravitational collapse [27]. The limitations of the approach in Ref. [26] regard the near-brane approximation of the black string [15,28], whose global form for arbitrary distances from the brane into the bulk, cannot be entirely provided. Although globally limited by the range of validity of a Taylor expansion into the bulk, the MGD black string revealed in Ref. [26] prominent features, beyond the standard (Schwarzschild) black string [18,29]. In fact, for some eras of the universe, the MGD black string event horizon was shown to have a finite range, since it has a collapsing throat in the bulk that collapses to a point along the extra dimension [26].
More precisely, generalized black strings features have been widely investigated, mainly in realistic brane-worlds with variable tension [15,30], also in the context of the MGD extensions [26,31]. However, all these approaches are valid near the brane, for a Taylor expansion of the brane metric onto the bulk is a good approximation for the purposes therein. For particular cases, the perturbative setup provides very interesting and relevant solutions for bulk regularity [15,30,32]. The Plateau-Rayleigh instability was also emulated in the MGD black string in Ref. [26]. Near enough to the brane, the Taylor expansion of the bulk metric -as a function of the variable brane tension, the 5D bulk Weyl tensor and the 5D cosmological constant as well -is an excellent approximation that provides the black string dripping into higher entropy 5D black holes states, after the formation of a black string throat along the extra dimension at finite cosmological times, also providing hints for its final state. Nonetheless, this phenomenon can just occur in a region of the bulk that is very near the brane to make the Taylor expansion to have a controlled error. Ref. [26] was the first to study the MGD black string, which was shown to have a finite extent along the extra dimension. Due to its intrinsic features, the original MGD procedure originated a MGD black string, whose analysis near the brane was an exact approach for its global behaviour [26,31], owing to a black string throat that is constricted by the variable brane tension as the 4D universe expands, reaching a event horizon with zero radius at finite cosmological time. Although valid for specific eras of the 4D universe, a final and general solution for the MGD black string stability problem lacks, still.
Here a complementary paradigm for performing the Gregory-Laflamme procedure shall be used, to derive the MGD black string stability analysis. Besides, the MGD solution is here shown to be a solution in modified gravity. The pivotal point to study the MGD black string stability is the Gregory-Laflamme procedure itself [20,33,34]. The method to be used is a correspondence between fluctuations about a 5D black string solution, which was shown to be related to black hole instabilities in 4D massive gravity theory and to black hole instabilities in 4D higher-derivative gravity theories [35][36][37]. Besides, black strings can be further be placed into the membrane paradigm, where the fluid/gravity correspondence plays a prominent role [39,40].
In order to study the MGD setup as a solution of 4D actions for the gravity, more intricate actions, beyond the Einstein-Hilbert one, can involve higher order derivatives, being obstructed by the Ostrogradsky theorem, proving the instability of some higher derivative theories of gravity. The ones that circumvent this theorem can play prominent roles as theories of gravity, due to their renormalizability, as for example the 4D Gauss-Bonnet Lagrangian, whose spacetime integral is a topological invariant. Ref. [19] introduced invariants constructed upon quadratic terms involving the Ricci tensor and scalar curvature.
Negative-eigenvalue eigentensors of the Lichnerowicz operator play a central role in the analysis both of static 4D non-Schwarzschild black hole solutions [41,42] and also in the analysis of time-dependent black-hole instabilities. Gregory-Laflamme instabilities [20] regarding fluctuations around MGD black string solutions are going to be studied in the context of black hole instabilities in 4D dimensional higher-derivative gravity theory [37]. However, in this work we shall not study the Gregory-Laflamme instabilities via bulk effects, but rather from their black hole effects on the brane using Ricci quadratic gravity theory. Linear perturbations around black hole solutions of massive gravity theories exhibit an unstable mode analogous to the Gregory-Laflamme instability of Schwarzschild black strings [35].
This paper is organized as follows: in Sect. II, the MGD is briefly reviewed and Sect. III is devoted to the study of the eigenmodes of the Lichnerowicz operator, in a Gregory-Laflamme analysis of the MGD black holes, constructed upon perturbing the Schwarzschild solution, using the (inverse of the) brane tension as the parameter driving the perturbation, originating the MGD solutions in the 4D Ricci quadratic gravity. The associated Lichnerowicz eigenvalue problem is split into time-dependent transverse traceless modes, being lead to a 1D radial wave equation with a potential, in tortoise coordinates. The results in a variable brane tension scenario are then compared to the GR limit. Sect. IV describe the Gregory-Laflamme analysis of the MGD black strings. Sect. V is dedicated to the concluding remarks and perspectives.

II. MINIMAL GEOMETRIC DEFORMATION (MGD)
The MGD procedure shall be briefly reviewed, in the original setup that induces high energy corrections to the GR on a 4D brane, percolated by a bulk Weyl fluid [1][2][3][4]. Greek indexes run through spacetime (brane) quantities, whereas capital roman indexes denote the bulk ones. The 4D Einstein equations encode an effective approach, when the Gauss-Codazzi method is applied to the 5D bulk Einstein equations with cosmological constant, R AB − 1 2 Rg AB = −Λ 5 g AB , where Λ 5 is the 5D cosmological constant, the g AB denote the bulk metric components, R AB is the Ricci tensor and R = g AB R AB . The brane tension shall be denoted by σ, hereon. The projected equations onto the brane, after the Shiromizu-Sasaki-Maeda procedure [17], can be expressed in the standard form where κ 2 4 is the 4D coupling constant and G µν = R µν − 1 2 Rg µν is the Einstein tensor, and the effective stress tensor reads, The brane tension of fluid membranes can explain the tiny value of the cosmological constant, based upon its small brane tension fluctuations [13]. The analogy is based upon a correspondence between statistical mechanics and quantum field theory. The description of a brane as a spacetime corresponds to the geometric portrait of spacetime itself as a 4D manifold, whose quantum fluctuations do correspond to the brane thermal fluctuations [14,16]. More precisely, in the membrane paradigm, the brane itself emulates the 4D universe wherein we live, having its self-gravity manifesting as a brane tension. The general-relativistic framework is based on an infinitely rigid brane, where σ → ∞. However the phenomenological bound for the brane tension σ 3.184 × 10 6 MeV 4 must be satisfied, as shown in Ref. [43] in the context of the MGD setup. The model of a fluid brane is better approached by considering an Eötvös fluid brane, for which the fluid dynamical brane tension depends on the brane temperature [44][45][46][47][48]. In fact, the universe has cooled down as it expanded, and the gradient of temperature is related by the time elapsing, also agreeing with the cosmological standard model.
The first component T µν in Eq. (1) denotes the brane matter stress tensor accounting a perfect fluid, T µν = ( + p)u µ u ν − pg µν , with as the density, p the pressure scalar field, u µ is the 4-velocity and g µν is the brane metric. The second term of (2), the tensor controls the local dominant high-energy regime, when σ. Besides, E µν represents the bulk Weyl tensor projection onto the brane, being related to eigenmodes of bulk gravitons [18]. For spherically symmetric, static, metrics, it reads where U is the bulk Weyl scalar, h µν is the projection tensor along the fluid velocity, and P is the stress tensor. Let us consider a static and spherically symmetric matter distribution, where dΩ 2 is the areal element of the 2-sphere. The spherical coordinates x 0 = ct, x 1 = r, x 2 = θ, x 3 = φ shall be used. Solving Eqs. (1), using the metric (5) and the tensor energy momentum (2), relates the density to the pressure of a fluid. Besides, the radial and tangential components of the effective pressure of the Weyl fluid read [2][3][4] p a = p + σ −1 where a = 0, 1. Effects of the 5D Weyl fluid (6) produce anisotropy in the projected fluid on the brane, which arises from the difference between the components of the pressure of the fluid, proportional to the inverse of the brane tension σ −1 . Hence, the isotropy of the perfect fluid is restored in the GR limit σ −1 → 0, namelẙ p r −p t = 0, as can be forthwith read off Eqs. (2) and (4), yielding Eq. (2) to read T µν = T µν . In order to obtain the deformation term in the metric, the temporal component is fixed, to deform the radial component g 11 (r)(= g rr (r)) = B(r), for r > R where M is the stellar distribution Misner-Sharp mass, R is average radius of the distribution and the parameter ς encodes a bulk-induced deformation of the vacuum, at the compact distribution surface r = R [2,4]. The MGD parameters are observationally constrained by the classical tests of GR [49] and here the particular case where shall be assumed, for b ∼ 1.348 a constant parameter [4,8,49]. The negativeness of ς compels the MGD star gravitational field to be weakened, as an effect of a finite brane tension, σ, when compared to the GR regime σ → ∞. Since ς ∼ σ −1 by Eq. (10), Eq. (9) can be written as for some radial function f (r), as indeed expected by the MGD procedure, since at the low energy regime σ −1 → 0, the Schwarzschild solution of the GR is recovered. Eq. (9) shows that the B(r) function admits two coordinate singularities, Just R 1 is a physical singularity. The scalar curvature, evaluated at R 2 , is finite. Moreover, since the solution is an analytical one up to the order O(σ −2 ), the horizon r 1 shall not be a physical singularity when higher order terms in the brane tension are taken into account.
Higher order Kretschmann scalars are also finite for R 1 and do not vanish, evaluated at R 1 . According to Ref. [2], it yields R 1 > R 2 , thus making R 2 unaccessible for observation, in the case of MGD black holes. Meanwhile, many applications were implemented in the MGD setup, mainly in stellar systems, also regarding MGD glueball stars and their improved range of observability at the LIGO and eLISA experiments [50], studies of acoustic perturbations in gas fluid flows [40] and their observability by lensing effects [51] as well.
Introducing the parameter the scalar curvature for the MGD metric (9) reads The next section is, in particular, devoted to place the MGD solution as a 4D one, in the Ricci quadratic gravity, besides studying their stability.

III. LICHNEROWICZ MODES AND STABILITY ANALYSIS OF MGD BLACK HOLES IN RICCI QUADRATIC GRAVITY
The 4D Ricci quadratic gravity is described by the Lagrangian [19], yielding its derived equations of motion to read Letting R µν = 0 after accomplishing the variation implies that [19] δR µν + ( whose trace yields δR−2[3α 2 +α 1 ] δR = 0. These equations have black solutions that can be lead, in particular to the Schwarzschild black hole for α 1 = 0 = α 2 . Integrating in the range between the black hole event horizon and the infinity yields δR = 0, whenever one chooses the correct sign of the term δR that corresponds to m 2 0 > 0 [19], for the for scalar mode [41,52]. Replacing δR = 0 into Eq. (17) yields where is the Lichnerowicz operator acting on any tensor field f αβ . The tensor δR µν was shown in Ref. [19] to be a transverse traceless tensor. Eq. (19) provides a transverse traceless eigentensor ψ µν of the Lichnerowicz operator, with eigenvalue λ = 1/α 1 . It provides a perturbation of the Schwarzschild solution, that here shall be identified with the MGD solution (8,9), whenever the perturbation parameter is identified to the MGD parameter in Eq. (10). The constraint α 1 < 0 must be imposed for massive spin-2 modes that are not tachyonic modes [19].
The Gregory-Laflamme analysis of the MGD black strings can be implemented by starting with the derivation of the negative eigenmodes of the Lichnerowicz operator. Ref. [19] made explicit the GR limit, σ −1 → 0, when B = A in Eq. (5), letting where γ ij = r 2 0 0 r 2 sin 2 θ and the involved functions, except the γ ij , are r-dependent one-variable functions. The requirement g µν ψ µν = 0, yields −ξ + χ + 2 ψ = 0. Besides, the transversality condition ∇ µ ψ µν = 0 implies that [19] The Lichnerowicz operator action on the ψ µν modes, and the Schwarzschild solution limit σ −1 → 0 with A(r) = B(r) in Eq. (5), together with Eq. (22), yield [19] L ψ 00 = −A 2 d 2 Taking into account λ = 1/α 1 in Eq. (21), and solving the tracefree and transversality conditions for ξ and ψ yields a radial ODE, The transverse traceless Lichnerowicz modes associated with the solutions of Eq. (28) of the form (22) can be obtained, for large distances, beyond the horizon event. In this limit, Eq. (28) reads χ + 4 r + λχ = 0, yielding the perturbation function χ = a 1 exp( √ −λr) + a 2 exp(− √ −λr), for a 1 , a 2 constants. It is related to one renormalizable mode, for λ < 0. Ref. [19] numerically showed that in full compliance to Ref. [53]. A transverse and traceless Lichnerowicz mode, with eigenvalue λ < 0 was shown to yield black hole solutions of type (5), ramificating from the Schwarzschild black hole ones when α 1 = λ −1 [19]. Indeed, transverse traceless perturbations of the Schwarzschild metric, g µν → g µν + h µν yielding δR µν = 1 2 L h µν . The MGD metric (5) can be implemented in the Ricci quadratic gravity setup by a linear perturbation of the standard GR solution (5), with respect to the inverse of the brane tension, whereÅ andB are functions to be determined a posteriori, and terms O(σ −2 ) are disregarded, due to the most recent bounds for the brane tension [43]. In fact, the underlying paradigm of the MGD procedure is based upon this consideration, as terms of order O(σ −2 ) are shown to be phenomenologically irrelevant [43]. Hence, the solution (30,31) is as analytical as the MGD itself. The perturbation (30,31) emulates the one in Ref. [19] with arbitrary parameter, but here the inverse of the brane tension is such a function that drives the perturbation. Eqs. (30,31) are shown in what follows to be equivalent to the MGD metric (5,8). In fact, Eqs. (21) taken together with δR = 0 yield the following equations: 2r 2 AÅ +r A+3 Å +r A+1 (ς B +ςB )+4ςB = 0 , (32) and 4r 2 rAÅ +B + α 1 2r 2 A 3A − 1 (ς B + 2ς B + ςB ) Our aim is to show that the MGD solution (8,9) satisfies the ansätzen (30,31). The parameter ς ∝ σ −1 in Eq. (10), regarding the brane tension, is well known to vary and cannot be taken as truly being a constant. Successful phenomenological models are led by Eötvös fluid branes, which takes into account the cosmological evolution and considers a timedependent variable brane tension that avoids any global anisotropy [14,16]. However, Eötvös fluid branes are based on FRW cosmologies, being ab initio isotropic. In order to encode CMB (cosmic microwave background) anisotropies, the brane tension can be more generally assumed to be r-dependent, without loss of generality, dismissing at a first analysis any angular dependence. In this approach, the brane tension can be made to depend upon the position of the black hole on the brane, impacting the Lichnerowicz modes. In fact, Refs. [54,55] show the relationship between bulk graviton massive modes stresses the brane shear anisotropy, using CMB to impose bounds on the anisotropy originated from these stresses. CMB anisotropies can be, thus, formulated in brane-world models [18,54,56], where bulk effects can determine the CMB anisotropy. Since inhomogeneity and anisotropy are observationally small, cosmological effects of anisotropy can be encrypted into the variable brane tension, as it describes the brane self-energy [56,57]. This setup implements the brane with inflation [58].
Eqs. (32) and (33) are satisfied for the MGD metric (5,8), when and the bound of the brane tension σ 3.184×10 6 MeV 4 is satisfied. It makes the ς-family of solutions (30,31) to be given by the MGD solutions (8) and (9), respectively. The classical tests of GR were used in the context of the MGD setup, to show that the parameter in Eq. (13) satisfies the bound (2.80 ± 3.45) × 10 −11 [49]. Besides, the bound of the brane tension σ 3.184 × 10 6 MeV 4 holds [43]. Ref. [59] proved that these effects are more feasible for mini-black holes, working as the MGD of the Schwarzschild solution. In fact, taking into account Eqs. (10,13), for a black hole of mass 10 9 M , where M denotes the solar mass M ≈ 1.989 × 10 30 kg, for (2.80±3.45)×10 −11 and σ 3.184×10 6 MeV 4 , it implies that the MGD of the radial Schwarzschild metric component reads B(r) = 1 − 2GM where α ≈ 10 −13 , for supermassive black holes. That is the reason why we can perform our perturbed solution. The effects of accretion onto the black hole, due to the MGD, might be substantial for mini-black holes, whose analysis was studied in Ref. [59].
The MGD black hole solutions (8,9) can be identified to the new perturbative black hole solutions (30,31) of Ricci quadratic gravity, whenever the perturbation parameter in Eqs. (30,31) is regarded as the inverse of the brane tension. Since the brane tension is a fundamental scalar field playing the role of the brane self-energy, it can be seen as a bath emulating either the cosmological constant or a minimum of a scalar potential that permeates the spacetime [58].
To plot the figure below, natural units are adopted where G = c = 1 = M , without loss of generality, following the computational routine in the Mathematica given in Ref. [61]. Following the analysis of the Schwarzschild perturbed black hole in Refs. [19,62], a similar qualitative profile is seen in Fig. 1. Although the finite value of the brane tension provides a slightly different plot, with respect to the GR limit σ −1 → 0, the asymptotic limit r → ±∞ is still similar to the GR limit. In fact, when r → −∞ both plots in Fig. 1 approach to 0, whereas the asymptotic limit when r → ∞ is a constant different of zero. The phenomenological finite upper bound for the variable brane tension then alters the exact value of the asymptotic constant for the plot in Fig. 1, when r → ∞. Whilst the qualitative profile of both finite and infinite (Schwarzschild) values of the brane tension are similar, the finite value of the brane tension imposes a different amplitude for the asymptotic outgoing waves e iω(t∓r ) .

IV. GREGORY-LAFLAMME ANALYSIS OF MGD BLACK STRINGS
The instability of the 4D Schwarzschild solution in Ricci quadratic gravity is directly related to the 5D black string instability through the Gregory-Laflamme procedure [20]. Transverse traceless linear perturbations are governed by a negative-eigenvalue mode of the Lichnerowicz operator action on Schwarzschild-type backgrounds [19,53], which was previously analyzed in Eq. (29). Stability properties of the MGD black hole solutions can be derived from the Lichnerowicz negative eigenvalue, that fixes the fringe of the MGD solutions (8,9), seen as a ς-family of Schwarzschild black holes (30,31), minimally geometric deformed. This is consistent with a variable brane tension, as previously discussed.
This was accomplished for the Schwarzschild case in the context of the 4D Ricci quadratic gravity. Now, taking into account the 4D static spherically symmetric metric ds 2 in Eq. (5), with the MGD coefficients (8,9), the MGD black string metric reads ds 2 = ds 2 + dy 2 , where one denotes by y the coordinate into the bulk, usually taken as the Gaussian coordinate [15,18]. In fact, a black string is, essentially, a higher-dimensional generalization of a black hole, whose event horizon has topology homeomorphic to S 2 × S 1 , also having asymptotical topology M d−1 × S 1 . A black string is a physical solution in lowenergy string theory, described by an extended gravitational object with an event horizon. [63]. The prototypical black string metric is a solution to the Einstein equations in 5D (or even higher-dimensional) gravity, being split into a Tangherlini-like black hole and an extra flat compact dimension. In the context of the braneworld scenarios, the matter trapped onto a brane might undergo gravitational collapse, where a black hole with the extended horizon into the transverse extra direction shall be produced. An observer on the brane shall realize a black hole on the brane, corresponding to a black string in the context of D-dimensional gravity, with metric where the g µν corresponds to the black hole geometry on the brane, and the coordinate y denotes the extra dimension.
The MGD black string analysis straightforwardly provides the determination of MGD black string instabilities, in a setup based upon massive gravity theories. In fact, a bulk massive gravity theory has transverse perturbations where the ψ µν are given by Eqs. (36,37) and driven by the Lichnerowicz condition ( L + ω 2 + m 2 ) ψ µν = 0 , where m denotes the mass of the 5D spin-2 bulk gravitons [35,36]. It shows the identity for time-dependent fluctuations between the MGD black string and the quadratic Ricci gravity, whenever ω 2 = −λ = − 1 α1 . The results assume the form (46), where the transverse traceless modes are given by Eqs. (36), (37), and (46).
Taking the metric (5,8,9) of the MGD black hole, with event horizon 2GM c 2 , the time dependent eigenmodes of the Lichnerowicz operator exponentially grow, if Since λ < 0 by Eq. (29) [19], the criterion for stability reads yielding complete stability, for all ω to be determined by M 2 − 4c 4 λ G 2 m 2 . This condition cannot be satisfied for massless gravity (m = 0), yielding the standard Gregory-Laflamme instability [20]. Nevertheless, if MGD black strings in 5D massive gravity theory, then for a MGD black hole event horizon 2GM c 2 , satisfying (48), the MGD black string is precluded from a Gregory-Laflamme instability.

V. CONCLUSIONS
We showed that the stability of MGD black strings can be analyzed from the deformation process that is inherent from the MGD construction itself. Different of the Gregory-Laflamme instability of the classical black string, there is a classical MGD black string mass above which the MGD black string is stable, under long wavelength perturbations. The variability of the brane tension is crucial in this point, for it to satisfy Eqs. (32,33), that cannot hold for constant values of the brane tension.
The MGD black hole solutions were shown to be solutions of 4D Ricci quadratic gravity. The Lichnerowicz operator negative eigenvalue rules the stability of the black hole, under perturbations, when its associated wave equations are identified to the Schrödinger-like 1D problem (40), with the potential (41). The results are plot in Fig. 1, that shows this potential for the Schwarzschild deformation in Ref. [19] and for the MGD, brane tensiondependent, solutions, for the maximum value of the brane tension in the phenomenological bound σ 3.184 × 10 6 MeV 4 [43]. Although the derived profile of both finite and asymptotic values of the brane tension are similar in the two plots in Fig. 1, the finite value of the brane tension was shown to yield different amplitudes to the asymptotic outgoing waves, when compared to an infinitely rigid brane.
Eq. (12) determines the event horizons of the MGD standard black hole. As the outer event horizon equals the Schwarzschild one, the study of the stability of the MGD black holes cannot be emulated from the Schwarzschild procedure [19], since it yields the same result. This intricate question may be solved and can be implemented through the MGD extensions [31], that is the most straightforward candidate, beyond the standard MGD procedure itself, to play the role of a different event horizon in the ansätzen analogous to Eqs. (30,31), for the extended MGD setup. However, this approach has been already revealed to involve hard computational tasks that are being still managed, with relative success [64]. Hence, a more speculative question regards the Bardeen-Cooper-Schrieffer (BCS) and its Cooper instability. In fact, Ref. [65] points to strange metals in condensed matter, driven by their duality to 4D Reissner-Nordström black holes. Holographic superconductors were shown to be dual to the Reissner-Nordström black hole, whose hair is implemented by a Higgs field that composes an atmosphere girdling the black hole, that has an instability shown to be dual to the BCS instability. We then want to answer whether the extended MGD instability can be dual to any condensed matter system.
The next steps also naturally consist of analyzing dilatonic brane-world models and the deformation of their black hole solutions, where the variable brane tension is a scalar field depending upon bulk scalar fields. In this scenario, the AdS/CFT was shown in Ref. [58] to be placed into brane-world models, where the effective stress tensor has some terms beyond Eq. (2) that involve the bulk dilaton scalar field and its derivatives. This extended tensor can be, hence, interpreted as the stress-energy tensor of the CFT matter, whereas the brane tension can be derived, via the Israel-Darmois junction conditions, as the Lie derivative of the bulk scalar field.