MGD-decoupled black holes, anisotropic fluids and holographic entanglement entropy

The holographic entanglement entropy (HEE) is investigated for a black hole under the minimal geometric deformation (MGD) procedure, created by gravitational decoupling via an anisotropic fluid, in an AdS/CFT on the brane setup. The respective HEE corrections are computed and confronted to the corresponding corrections for both the standard MGD black holes and the Schwarzschild ones.


I. INTRODUCTION
corrections, for MGD-decoupled spacetimes generated by anisotropic fluids, are derived and discussed in Sect. III in two important regimes, considering the boundary far from the event horizon and almost on it. More analysis and discussions about these results are scrutinized in Sect. IV.

II. MGD-DECOUPLING AND ANISOTROPIC BLACK HOLES
The MGD procedure can be realized as a mechanism that is usually employed to derive high energy corrections to GR. Denoting by σ the brane tension, systems with energy E σ neither feel the self-gravity effects nor the bulk effects, which then allows the recovery of GR in such energies. An infinitely rigid brane scenario, representing the 4-dimensional GR brane pseudo-Riemannian manifold, can be implemented in the σ → ∞ limit. The most strict brane tension bound, σ 2.81 × 10 −6 GeV 4 , was derived in the extended MGD-decoupling context [16].
Starting from any straightforward static, spherically symmetric, source of gravity, corresponding to some stressenergy tensor T µν , one can iteratively introduce more intricate gravitational sources. This procedure is described by [2,3] where α (m) are constant parameters that drive the weakly-coupled effects of the sources T (m) µν . This procedure holds whenever the backreaction among the sources is disregarded, namely ∇ µ T (m) µν = 0, for any natural number m. Hereon natural units will be used.
The MGD decoupling for a perfect fluid can be introduced by considering Einstein's equations where κ 2 denotes the Newton coupling constant. One assumes that the total stress-energy tensor is endowed with two contributions,T where T µν = (ρ + p) u µ u ν − p g µν is the 4-dimensional stress-energy tensor for a perfect fluid with 4-velocity field u µ , density ρ and pressure p. The term Θ µν in Eq. (3) describes an additional source whose weakly-coupling to gravity is driven by α, inducing anisotropy. As the Einstein tensor does satisfy the Bianchi identity, the source in Eq. (3) also satisfies the Eq. (1). Given a static, spherically symmetric, metric ds 2 = −e ν(r) dt 2 + e λ(r) dr 2 + r 2 dθ 2 + r 2 sin 2 θ dφ 2 , the fluid 4-velocity reads u µ = e −ν/2 δ µ 0 , for r ∈ [0, R], where R = ∞ 0 ρ(r )r 3 dr / ∞ 0 ρ(r )r 2 dr emulates the compact star surface radius. The metric (4) must satisfy Einstein's equations (2), yielding where f = df /dr. Besides, the equality θ 3 3 = θ 2 2 comes from spherical symmetry. Eq. (1) yields [2,3] p + ν 2 (ρ + p) + α ν 2 θ 0 0 − θ 1 The limit α → 0, corresponding to the perfect fluid case, is then formally redeemed. One can define the effective density, and the effective radial and tangential pressures, respectively, by [2,3] It leads the Θ µν tensor to induce a coefficient of anisotropy of a compact stellar configuration, The MGD-decoupling can now be applied to the case at hand by simply noting that the stress-energy tensor (3) is precisely of the form (1), with T µν as the one of a perfect fluid, α (1) = α and T (1) µν ≡ 0, for all m ≥ 2. The components of the diagonal metric g µν that solve the complete Einstein equations (2) and satisfy the MGD readg µν = g µν = g (1) µν , for µ = ν = 1, andg 11 = g 11 + α g (1)11 [2,3]. Hence, solely the radial metric component carries signatures of Θ µν . One can solve the Einstein's equations for a perfect fluid T µν , and then the remaining quasi-Einstein equations for the source Θ µν , where the divergence-free quasi-Einstein tensorG µν = G µν + Γ µν , with Γ µν = Γ µν (g ρσ ) denotes a metric dependent tensor that is divergence free [2,54]. A feasible solution of the coupled system in (11), for a perfect fluid, reads where m(r) denotes the Misner-Sharp mass, r-dependent, function, which by measuring the amount of energy within a sphere of areal radius r, provides a coherent quasilocal definition to the curvature-producing energy in black holes. The source Θ µν effects on the perfect fluid solution {ξ, µ, ρ, p} can be then encoded in the MGD into the radial component of the perfect fluid geometry in (13). Namely, the general solution is given by Eq. (4) with ν(r) = ξ(r) and where f = f (r) is the MGD function to be determined from Eqs. (12), given by [54] Black hole solutions are derived from an EoS when one determines the vacuum MGD function f . The MGD metric will therefore read For the Schwarzschild solution, the surface r H = 2M is both a Killing horizon and an exterior marginally trapped surface [54]. For the MGD-decoupled Schwarzschild metric (19), the component g tt (r) = e ν(r) equals the Schwarzschild standard form, having coordinate singularity at r = r H , being also a Killing horizon. However, the causal horizon at r = r h is such that g rr (r h ) = e −λ(r h ) = 0, or equivalently, One then demands that r h ≥ 2M , in such a way that the surface r = r H is either concealed behind, or at most coincides with, the causal horizon. Eq. (19), corresponding to the MGD-decoupled metric, can represent a black hole only if the causal and the Killing horizons coincide, namely, r h = 2M = r H . This condition will be assumed in what follows. Now let one considers an anisotropic configuration, governed by a EoS of type [54] with a and b constants. Conformal configurations are obtained by setting a = −1 and b = −2, whereas barotropic configurations require the choice a = −1/K and b = 0 [54]. Eqs. (16 -18) then yield the ODE for the MGD function whose general solution for r > r H = 2M is given by where > 0 represents a length scale, and with obviously b = 2. The range A > 0 is required for deriving an asymptotic flat solution. Hence, The effective density, and the radial and tangential pressures are respectively given bẙ They diverge at r c = BM, that is a physical singularity at the range 0 < r c < r H for 0 < B < 2. For B > 2, r c represents a physical singularity outer the Killing horizon, r H , being this case physically forbidden, to preclude naked singularities. Moreover, the effective radial and tangential pressures are related bẙ Since A > 0, one concludes that both the radial and tangential pressures do have contrary signs. In addition, the MGD effective density and radial pressure satisfẙ Since A > 0, the asymptotic behaviour in Eq. (30) demands the range A ∈ (0, 1], for the density does not change sign in the sector r ∈ (2M, +∞). For α negative, the effective density is positive. For b ∈ (2, 4), the coordinate singularity r c attains negative values. Hence, no extra singularity, besides the usual Schwarzschild one at r = 0, exist, if there are no solutions, r 0 > 0, of the equation e −λ(r) = 0 [2]. This is satisfied whenever α > 0, for any A > 0, that is, for a > 1. There is also a second coordinate singularity, solution of e −λ(r) = 0, consisting of when α < 0. To produce a physically viable black hole solution, this solution r 0 have to attain lower values than the existing singularity. If b ∈ (2, 4), then the existing singularity r = 0 yields r 0 ≤ 0,. Equivalently, if and |α| satisfy Otherwise, if b < 0 or b > 4, the relevant singularity occurs at 0 < r c < r H , but r 0 > r c makes this case inappropriate. The final conclusion is thus that the linear EoS (21) always produces black holes, with a Schwarzschild-like physical singularity at r = 0, if b ∈ (2, 4) and a > 1.
Not being so generalist, and to further explore the physical content of this model, one adopts hereon a = 2 and b = 3, implying that A = 2 and B = −1. This choice allows no other additional physical singularity than the well known r = 0 one. In fact, Refs. [2,3] analyzed the particular case b = 3 with a > 1, to ensure that solutions are asymptotically flat.

III. HEE FOR A BLACK HOLE FROM ANISOTROPIC FLUID UNDER MGD
The entanglement entropy (EE), S A , of some manifold A, with boundary ∂A, is an important quantity in 4dimensional QFTs. Indeed, it represents the von Neumann entropy of the (reduced) density matrix, in the case where degrees of freedom, into a 3-dimensional space-like submanifold, B, are stretched out [64]. The EE quantifies how A is correlated to B, measuring the amount of entropy in A, by an observer isolated from B. Therefore, there is a part of the AdS 5 bulk from which one can compute S A in the gauge/gravity duality. At a zero temperature regime, any system in QFT is characterized by a pure state |Ψ , where the respective density matrix reads ρ = |Ψ Ψ|. When one splits the quantum system into A and B, an observer isolated from B describes the quantum system by a reduced density matrix This equality is violated by finite temperature QFTs. The inequality S A+B ≤ S A + S B is always valid. Considering a QFT on a 4-dimensional spacetime splitting, R × Σ 3 , into some timelike direction and a 3-dimensional spacelike manifold, Σ 3 , a 3-dimensional submanifold B ⊂ Σ 3 represents a geometric complement of A. An ultraviolet cut-off, b, prevents the EE to diverge. The coefficient that drives this divergence at the continuum limit depends on the area of ∂A, The conformal Poincaré metric, that characterizes AdS 5 bulk geometry, reads The CFT 4 that is dual to gravity in the bulk lives on the R 1,3 boundary, at z = 0, with coordinates (x 0 , x i ). The AdS 5 conformal coordinate, z, is the energy scale.
Besides the pure AdS 5 spacetime (34), AdS 5 black holes can be also regarded. Indeed, the boundary ∂A of A can be extended to some manifold γ A , whose boundary equals ∂A. Consequently, the EE S A in the CFT 4 reads [64][65][66].
where κ 2 5 is the 5-dimensional Newton coupling constant. A. HEE corrections

Far from the horizon
For computing HEE corrections, the first procedure consists to settle the boundary manifold in a sector that is very far away from the event horizon, named r = r L , which is still at a finite large distance. Now, let a circle in polar coordinates, be defined by the angle θ = θ 0 , carrying the role to confine the entangling surface. Additionally, let r = r(θ) be the areal radius coordinate describing the minimal surface, whose boundary depicts the the entanglement manifold.
One needs to derive the minimum of the area function, where the boundary condition, r(θ 0 ) = r L , which is required to find the minimal surface. Eq. (36) can be rewritten, when one substitutes y = cos θ, reading where y 0 = cos θ 0 , the derivative with respect to y was denoted by a dot and G = G (r(y)) ≡ exp [λ (r(y))]. After finding the global minimum of Eq. (37), one obtains the HEE by directly employing Eq. (35). Varying Eq. (37) with respect to r(y) yields which is a ordinary differential equation (ODE) that is highly nonlinear, whose most simple solution is derived by attributing G ≡ G (r(y)) = 1, i.e., r(y) = k 0 /y. According to [68], it is possible to attain nontrivial solutions of Eq. (38) through series expansions for both the G (r(y)) and r(y) functions, respectively expressed by considering as a dimensionless small parameter encoding black hole mass M -to-r L ratio, that is, = M r L . Just the O( ) terms in (39a) carry corrections when one considers a collapsing black hole case.
To determine the G function, up to the 2 nd -order, one has the following metric of this specific spacetime permeated by anisotropic fluid discussed in Sect. II: where dΩ 2 = dθ 2 + sin 2 θdϕ 2 and one denotes α 2 = −ω 2 . Therefore, one can achieve the g-functions, displayed in Eq. (39a), i.e., It is important to stress out that a dimensional analysis was necessary, to indicate the MGD parameter linked to the expansion parameter as ω = ψM . Clearly, it is always possible to construct higher orders terms in Eq. (39a), whenever necessary. Just for the current goal, one stopped at the second order. As it was demonstrated in [68,69], Eq. (38) is solved -order by order -employing g 1 (y) and g 2 (y) in Eq. (41).
To determine the modifications to the HEE up to the 2 nd -order, one needs to compute the r-functions as follows. Henceforth, one gets 1 st -order terms in after expanding Eq. (38), and using the g-functions listed in Eq. (41). Hence, the 1 st -order ODE can be expressed as The resolution of Eq. (42) brings the constant of integrations A 1 and A 2 , which have their values determined by precluding any type of singularity at y = 1, i.e., y = cos θ ∈ [cos θ 0 , 1]. Thus, it is imperative to establish A 2 = 2r L . Next, the boundary condition r 1 (y) = 0 implies that A 1 = − (y 0 + 2 log[y 0 /(1 + y 0 )]) r L . So, the first r-function reads r 1 (y) = r L 2 y − y 0 − 2 log 1 + y 1 + y 0 + 2 log y y 0 .
Importantly, there is a subtle restriction due to limitations in the perturbative expansion, as aforementioned in Ref. [68]. In this sense, one enforces that, once y = 0 cannot be ever achieved, then r 1 (y) is well defined in the range θ 0 < π/2 or, in an equivalent way, y ∈ (0, 1). Now, one goes to the 2 nd -order in , proceeding as early and employing the r 1 (y) function, displayed in Eq. (43), to determiner Proceeding analogously as it has been made to obtain r 1 (y), one has Again, the computation of the constant of integration, which is important to reach the r-function above, was realized by eliminating the divergences at y = 1 and the boundary condition r 2 (y 0 ) = 0. Finally, one realizes the expansion I = I 0 + I 1 + 2 I 2 , within the area formula showed in Eq. (37). Henceforth, the r-functions Eq. (43,45) need to be engaged to calculate each order of contribution for the HEE, i.e., S Anis = S Anis 0 + S Anis 1 + S Anis 2 + · · · , where we stand ourselves up to the 2 nd -order term. The expression for the 0 th -order is written as Besides, the expression for the 1 st -order reads whereas the 2 nd -order is given by Emulating some of the useful results in Refs. [68,69], one notices no difference confronting the 0 th -order terms of HEE for the MGD, the Schwarzschild or the spacetime permeated by an anisotropic fluid, thus Ξ 0 =Ξ 0 = 1. In addition, one recalls that where ξ is a parameter characterizing a pure MGD spacetime which carries the brane tension signature [69][70][71].
A comparative analysis can be established for the HEE corrections of those three spacetimes mentioned in the last paragraph. Henceforth, one defines the ratio, when likening the n th -order corrections of HEE for the spacetime described by an anisotropic fluid under the MGD approach and the pure MGD one, where S Anis n and S MGD n are the n th -order corrections of HEE for those spacetimes early mentioned. On the other hand, when ξ = 0, then S MGD n = S Schw n , which is the n th -order correction of HEE for the Schwarzschild spacetime. For this specific case, one definesΞ n = S Anis n /S Schw n . The next step can be accomplished, through the calculation of the quantities displayed in Eq. (50). Thereupon, in a straightforward way, one gets Ξ 0 = 1,Ξ 0 = 1 and whileΞ 1 = 1. Since ξ < 0, a higher value of the brane tension induces a lower Ξ 1 .
To continue the comparative analysis with the MGD spacetime and, consequently, with Schwarzschild ones, let one resets ω = ψM to determine As an upper bound, one has the saturation value −α 2 = ω 2 = M 2 , according to Eq. (32), which is equivalent to make ψ = 1 . Another relevant point regards the brane tension. Higher negative values of ξ do not bring meaningful influence to the ratio. Thus, one restricts the analysis dealing with the interval 0 ≤ ξ ≤ −1.
In Eq. (52), the ratio Ξ 2 is positive, whereas y 0 is close to zero. It is very important to emphasize that the limit interferes with the sign of Ξ 2 . Then, one assumes the lower limit of integration -a practical small numerical value that can be used -to investigate the behavior of such ratio. Another interesting point regards the brane tension itself. Bigger negative values of ξ do not bring meaningful influence to the ratio Ξ 2 . From now on, the respective analysis of Ξ 2 is restricted to the interval −1 ≤ ξ ≤ 0. Fig. 2 shows the profile of Ξ 2 under the influence of ξ, with ψ assuming different values.

FIG. 2.
On the left, the behavior of the ratio Ξ2 considering the brane tension varying in the range −1 ≤ ξ ≤ 0. On the right, the low influence of bigger negative values for the brane tension. We adopted y0 = 10 −10 6 . Fig. 3 clarifies how the ratio of the 2 nd -order correction is affected by the anisotropic fluid, regarding a fixed brane tension in MGD spacetimes. Of course, with ξ = 0 one recovers such correction for a Schwarzschild spacetime, which is also demonstrated.
FIG. 3. The thick line (the curve in black) shows the ratio of HEE 2 nd -order corrections between a spacetime permeated by an anisotropic fluid and a Schwarzschild one. The dashed line (the curve in green) displays the ratio with the brane tension in an intermediary value, whereas the dot-dashed line goes with the extreme value of the brane tension that has been considered here. In both fixed ξ-curves, one has the anisotropic indicator ψ varying from the null to the saturation value. Again, y0 = 10 −10 6 was adopted.
For smaller values of the anisotropy parameter ψ, the ratios Ξ 2 are bigger ones. Even though one considers such aspect, the maximum ratio,Ξ 2 = 1, occurs when the brane tension is infinite and ψ = 0. Within the chosen interval for ξ, the brane under the weakest tension (ξ = −1) leads to a ratio of around 37%, since ψ is closest to 0. Besides, S Anis 2 over S Anis 1 is negligible due to stands at -order, thus the first law of HEE still holds, since δS Anis = S Anis − S Anis 0 ∝ M even up to the 2 nd -order in the calculations. Holding the first law of HEE is a crucial feature of the HEE paradigm [68,69]. This low impact of the anisotropic fluid could be caused by the localization of the boundary. Therefore, the next section treats the opposite situation to localize the boundary, which means settling it very near to the horizon to observe whether the influence this source is meaningful.

Almost on the horizon
A definitive test must occur while the boundary is localized close the event horizon [68,69] for the background under the influence of an anisotropic fluid, which is characterized by the metric (40).
where r H = 2M and r U = ω − r H /2 . Apparently, there is an extra horizon in metric (53). However, sustained by Eq. (32), one has no extra singularity when r U < 0, which is equivalent to ω < M . One sets a boundary almost on the horizon adopting a small displacement like ρ 0 = √ r H , once 1, which models the entangling surface as a circumference, while θ = θ 0 . In this sense, such setup induces the following metric on the t-constant manifold, that is, which considers p ≡ p(ρ) = ρ 2 + r H and q ≡ q(ρ) = 4(ρ 2 + 3 2 rH) 2 (ρ 2 +2rH+rU)(ρ 2 +rH+rU) , while ρ ≡ ρ(θ). To determine the HEE, at first, it is necessary to compute ρ, which leads to the minimization of the surface area, i.e., whereĨ = pq(1 − y 2 )ρ 2 + p 2 1/2 . Once again, the change of variable y = cos θ, has been employed to get ρ ≡ ρ(y). Therefore, the minimization of Eq. (55) with respect to y gives The endeavor of looking for an analytical solution of Eq. (56) drops off, in front of the strong nonlinearity of such differential equation. Hence, one must implement a perturbative method, adopting the expansion since ρ 1 (y 0 ) = √ r H and ρ 2 (y 0 ) = 0, considering the boundary condition ρ(y 0 ) = 0.
The 0 th -order term in Eq. (57) is absent just to avoid an area having a value greater than one valued at the point (ρ = ρ 0 , θ = θ 0 ). Thus, searching for the ρ-functions up to 2 nd -order, one inserts Eq. (57) into Eq. (56). Thus, at 1 st -order in , one gets In this manner, one can notice the influence of the main parameters of the gravitational decoupling within the minimal geometric deformation under an anisotropic fluid as a source. One determines the solution for Eq. (58), which is with τ = 1 2 −1 + 1 − 4 κ ) and P τ (y) denotes the Legendre polynomial of 1 st kind. This solution is regular at y = 1, with boundary condition ρ 1 (y 0 ) = √ r H .
Moving to the next order, that is, towards the 2 nd -order in , one arrives at a Legendre equation that is similar to Eq. (58), reading which has a solution like ρ 2 (y) = BP τ (y). By the other side, a boundary condition ρ 2 (y 0 ) = 0 implies B = 0. Hence, ρ 2 (y) = 0, which simplifies the ρ-function uniquely to the 1 st -order in . As a final development, with the ρ-function in hands, one calculates the area of the entangling surface. Firstly, expanding the integrand of Eq. (55) in , which yields Secondly, one substitutes the Eq. (62) into Eq. (55), executes the expansion of Area, which is Area = Area 0 + Area 1 + Area 2 + · · · . Restricting to the 2 nd -order, the correspondent HEE has the following contributions: Without an efficient method to calculate analytical solutions ofS Anis 2 , a clever procedure to deal with it consists of using numerical computation techniques. Indeed, numerically, Fig. 4 shows two specific values of κ. In fact, the saturation value (κ s = 1.125) and the value κ 0 = 1 that recovers the Schwarzschild solution. A general scenario based on this order of correction is displayed in Fig. 5, considering a valid range for κ and the complete size of the subsystem y 0 . One can realize that theS Anis 2 maximum value occurs when one takes into account the lower limit of integration at y 0 ≈ −0.53 and the saturation value κ s = 1.125. The shape of the 3D surface showed in Fig. 5 has a similar behavior, for all κ inside the range of validity.
Regarding the ratio of 2 nd -order corrections of HEE, taking the entangling surface almost on the horizon, by a calibrated brane tension (represented by particular ζ values 2 ), whereasS Anis Sat 2 stands for the saturated 2 nd -order correction of HEE for the spacetime permeated by an anisotropic fluid. Besides, taking the saturation value κ s = 1.125, one gets ζ s = −1/9 ≈ − 0.1111, which means that one has a saturated correspondence between a background with an extra anisotropic fluid and a generic MGD spacetime. Each curve adopts a percentual attenuation of ζ s . One might see clearly the influence of the anisotropy while one looks for the saturation value κ s = 1.125 and comparing with the Schwarzschild case. Another important feature shown in Fig. 6 is that higher ratio occurs with the totality of the interval of integration instead of using the limit y 0 ≈ −0.53, where the maximum 2 nd -order contribution happens for the saturation value κ s = 1.125. Finally, Fig. 7 displays the complete spectrum of ζ values, regarding all sizes of the subsystem to widely clarify the behavior of the ratioΞ 2 . Clearly, one notices that the most accentuated growing of the ratioΞ 2 occurs in the sector with the largest size of the subsystem and values of ζ close to zero, which is profoundly related to high values of the brane tension.

IV. FINAL REMARKS AND CONCLUSIONS
The MGD procedure has been applied to derive new solutions, acting as a solid and reliable mechanism to calculate high energy corrections to the GR. Therefore, the addition of an extra anisotropic source into the stress-energy tensor can originate a family of anisotropic black holes. This important line is combined with another one, which is delineated by the HEE conjecture for asymptotically flat spacetimes. Hence, the HEE was calculated for an anisotropic black hole considering boundaries -the situs to settle the entangling surfaces -far from the horizon and almost on it, permitting to define corresponding dual field theories. In the absence of a complete analytical solution, an expansion was used to determine the HEE up to the 2 nd -order in both scenarios that those boundaries were settled up. In addition, comparative investigations were implemented, considering order-by-order corrections of the HEE for MGD and Schwarzschild spacetimes with respect to anisotropic black hole.
For a boundary localized far from the horizon, the HEE of an anisotropic black hole is different of the HEE for an original MGD black hole, after one looks their respective corrections beyond the 0 th -order. The 1 st -order of correction does not have any relationship with the MGD parameter driven by the brane tension, shown by Eq. (47). By the way, it is the same value for the same order of correction to the HEE for a Schwarzschild black hole. The parameter takes a large radial distance to set the boundary, then can be the reason for such equality. Even though, comparing the 1 st -order corrections of HEE between anisotropic black holes and MGD ones, there is an important contribution caused by the brane tension carried by the parameter ξ in Eq. (47), and has featured in the ratio Ξ 1 with its sole dependence on such parameter. Thus, smaller brane tension values rise the impact of an anisotropic source, in front of a standard MGD formulation, clearly demonstrated by Eq. (51).
Towards the next order, a wider scenario rises. In fact, two parameters enter in the respective ratio Ξ 2 , displayed in Eq. (52). First, the full scenario about the influence of the brane tension in the anisotropy parameter was studied, restricting to the interval −1 ξ 0, where it is possible to see the major contribution of the brane tension aggregated with the anisotropy parameter ψ. Fig. 3 illustrates the influence of ψ and ξ. Immediately, it is possible to notice that the brane tension is rather than the weighting of the source with an anisotropic fluid, while ψ is small. On the other hand, very close to the saturation value of ψ, all ratios go to 0, which implies the high weight of the anisotropic fluid in the attenuation of the 2 nd -order correction of HEE. Besides, considering the values at the extrema of the ξ-interval, the ratio was obtained, indicating that the lower the brane tension, the bigger the ratio is. The choice of the largest distance to settle the entangling surface was imperative to determine 2 nd -order correction of HEE for the anisotropic black hole, even less than the correspondent order of correction of HEE for pure MGD spacetimes under low brane tension values and also for the Schwarzschild black hole, prompting the calculations of HEE closer to the event horizon. At the end, the first law of HEE holds -since δS Anis ∝ M -which is a crucial point based on the HEE paradigm.
Almost on the horizon, the boundary there establishes a HEE up to 2 nd -order correction in an adapted -expansion, employed to compute the minimal area of a respective entangling surface. The 0 th -order of HEE matches with the same order of correction of the HEE for the MGD spacetime as well as for the Schwarzschild one. The vanishing of the 1 st -order is also determined in the present case. The 2 nd -order correction of HEE has been computed numerically, considering the saturation value κ s = 1.125 that corresponds to the upper bound of −α 2 M 2 , which has been imposed by the requirement of no extra singularity. It heads to an important restriction, when comparing with the MGD spacetime. In fact, additional anisotropic sources bring a restrictive range for brane tensions, that is, the saturation value κ s = 1.125 gives a lower bound ζ s = −1/9 working as a limit for the brane tension. Obeying this physical constraint, the curve of values for such order of correction is displaced when compared with that correspondent order of correction of HEE for a Schwarzschild spacetime, as shown in Fig. 4 and Fig. 5. This last one emphasizes the influence of the size of the subsystem with the lower limit y 0 . Concerning the maximum value of the 2 nd -order correction, it has occurred with an inferior limit at y 0 ≈ −0.53 (while y = 1 works as the superior one) to the size of the subsystem -see Fig. 4.
Adding on, the ratioΞ 2 is plotted in Fig. 6 to demonstrate the impact of the anisotropy parameter on the 2 nd -order correction of HEE, when compared to the same order of correction of HEE for both the Schwarzschild and the MGD solutions. The full size of the dual subsystem leads to the ratioΞ 2 as higher as the saturation value gets attenuation. If one looks the 2 nd -order corrections of HEE for a saturated anisotropic black hole and the Schwarzschild one, there is an 11% increment of the ratioΞ 2 with the largest possible size of the subsystem, against a 4.4% increment with the lower extremum at y 0 = −0.53, which means a 247% increment from the minor percentage. Besides, Fig. 7 reveals the accentuated growing ofΞ 2 after tackling the largest size of the subsystem and lower values of the anisotropy parameter κ, which is linked with the brane tension parameter ζ 3 .