Quasinormal modes and tidal responses of black holes in generic anisotropic matter environments

Abstract: We develop a perturbative framework for a black hole embedded in a generic, possibly anisotropic, matter environment under spherical symmetry. Our approach extends previous analyses restricted to vanishing radial pressure or to perturbative matter configurations. Within this framework, we derive an analytical generalization of the Einstein cluster that incorporates a polytropic radial pressure, and we investigate the properties of this solution. We show that both the geodesic structure and the axial quasinormal-mode spectrum remain predominantly governed by an overall gravitational redshift effect, while the radial pressure systematically enhances the environmental corrections. In contrast, the tidal Love numbers are substantially more sensitive, and can exhibit order-unity deviations, including vanishing and negative strictly static magnetic Love numbers for sufficiently large anisotropy. We present the full linearized equations, which can be applied to various extensions, including ringdown analysis and extreme-mass-ratio inspirals.

• The paper develops a relativistic perturbation framework for black holes in anisotropic fluids, including independent density, radial pressure, and tangential pressure and two axial-sector closure prescriptions.
• Radial pressure modifies axial quasinormal-mode frequencies through the leading redshift [1-(1+K)C, while realistic environmental corrections to geodesics and shadows remain highly suppressed.
• Magnetic tidal Love numbers provide the strongest environmental diagnostic, scaling as O(C) and showing prescription-dependent behavior, including a strictly static value that vanishes at K=5/17 and becomes negative above it.

Overview

This paper by Zhao and Pani develops a fully relativistic perturbative framework for a Schwarzschild-like black hole (BH) embedded in a generic, spherically symmetric, anisotropic effective fluid environment. The central motivation is that prior treatments of environmental effects on gravitational-wave (GW) observables — most notably the Einstein cluster (EC) construction — assume vanishing radial pressure, an idealization that fails for realistic dark matter (DM) halos with velocity dispersion, self-interactions, or superfluid behavior. The paper closes this gap in two ways: it derives the complete set of linearized axial and polar perturbation equations for a BH surrounded by a fluid with independent density $\rho$, radial pressure $p_r$, and tangential pressure $p_t$; and it constructs closed-form background solutions generalizing the EC to include a polytropic radial pressure $p_r = \mathcal{K}\rho^\gamma$. The analysis then quantifies how the radial sound speed parameter $\mathcal{K}$ modifies geodesic structure, axial quasinormal modes (QNMs), and magnetic tidal Love numbers.

Perturbative framework

The background is described by a static, spherically symmetric metric $$ds^2 = -e^{\nu}dt^2 + e^{\lambda}dr^2 + r^2d\Omega^2$$, sourced by an anisotropic stress-energy tensor with density, radial pressure, and tangential pressure. Introducing the mass function $M(r)$ via $e^{-\lambda} = 1 - 2M(r)/r$, the background reduces to the standard anisotropic Tolman-Oppenheimer-Volkoff (TOV) system, with the anisotropy parameter $\sigma = p_r - p_t$ entering explicitly. Since three equations govern five unknown functions, two constitutive relations must be supplied to close the system.

In the Regge-Wheeler gauge, perturbations split into axial and polar sectors. A key structural result is that the axial sector decouples from matter fluctuations only up to a source term proportional to $\sigma X$, where $p_r$0 is one of the extra fluid variables introduced by anisotropy. Consequently, unlike the isotropic case, the axial master equation is not automatically closed: a closure prescription for the fluid variables is required. The authors consider two physically motivated choices — the strictly static case ($p_r$1) and the irrotational case ($p_r$2, $p_r$3), the latter imposing vanishing vorticity on both $p_r$4 and $p_r$5. Each yields a Schrödinger-like master equation with effective potential

$p_r$6

$p_r$7

The full polar-sector equations are provided in supplementary material but not analyzed numerically. An appendix verifies that the formalism smoothly recovers the vacuum (Regge-Wheeler), isotropic (Kojima–Pani-type), EC, and weak-environment limits of previous work, including reproducing the multi-parameter expansion of Datta et al. at leading orders.

Background solutions and energy conditions

Adopting the EC mass profile $p_r$8 together with a polytropic radial equation of state, the authors obtain analytic metric functions for $p_r$9. For $p_t$0, $p_t$1 admits a direct interpretation as the squared radial sound speed, $p_t$2, while simultaneously controlling the anisotropy profile; $p_t$3 recovers the EC and $p_t$4 gives local isotropy.

A notable theoretical result concerns energy conditions. While the weak, null, and strong energy conditions hold everywhere, the dominant energy condition (DEC) is violated near the horizon, with the violation onset depending only weakly on halo parameters in the realistic regime of compactness $p_t$5. Requiring DEC validity up to the light ring imposes $p_t$6; requiring it only up to the innermost stable circular orbit (ISCO) yields $p_t$7, i.e., the causality bound. The authors argue these violations are not fatal because DM spike models predict matter depletion within roughly $p_t$8, precisely where the effective-fluid description is expected to break down anyway. This is nonetheless an assumption-laden argument: the model's physical consistency inside this radius rests on the coincidence between DEC violation and the expected breakdown region, which is plausible but not demonstrated from first principles.

Geodesic structure

Expanding to third order in the compactness $p_t$9 (realistic environments have $p_r = \mathcal{K}\rho^\gamma$0), the light-ring radius, orbital frequency, Lyapunov exponent, critical impact parameter, and ISCO all exhibit the same hierarchy: the leading correction at order $p_r = \mathcal{K}\rho^\gamma$1 is a pure global gravitational redshift, with genuine modifications to the orbital structure appearing only at $p_r = \mathcal{K}\rho^\gamma$2. Radial pressure systematically amplifies both effects through the replacement $p_r = \mathcal{K}\rho^\gamma$3. Two concrete consequences follow. First, the ISCO receives no $p_r = \mathcal{K}\rho^\gamma$4 correction for any $p_r = \mathcal{K}\rho^\gamma$5 — its first shift appears at second order. Second, the shadow deformation scales as $p_r = \mathcal{K}\rho^\gamma$6, rendering environmental effects on current shadow observations negligible.

Axial quasinormal modes

Because the effective potentials satisfy $p_r = \mathcal{K}\rho^\gamma$7, the QNM frequencies redshift as

$p_r = \mathcal{K}\rho^\gamma$8

consistent with the eikonal relation linking QNMs to light-ring quantities. Sixth-order WKB computations of the fundamental $p_r = \mathcal{K}\rho^\gamma$9 mode confirm this analytical prediction accurately for both closure prescriptions. Because $\mathcal{K}$0 is small in realistic configurations, the strictly static and irrotational cases give identical leading corrections, indicating a degree of universality: the axial ringdown response is governed predominantly by the background geometry rather than by the detailed dynamics of matter perturbations. Within the allowed parameter range ($\mathcal{K}$1), radial pressure can enhance the deviation from the vacuum spectrum by an order-unity factor relative to the EC case — a meaningful amplification, though the absolute effect remains suppressed by $\mathcal{K}$2.

Tidal Love numbers

The tidal response emerges as the most sensitive probe. Computing the static ($\mathcal{K}$3) axial response under each closure prescription, and adopting the normalization scale $\mathcal{K}$4 advocated as more natural in recent literature, the $\mathcal{K}$5 magnetic Love numbers become

$\mathcal{K}$6

Three results stand out. The strictly static Love number vanishes at $\mathcal{K}$7 and becomes negative for larger $\mathcal{K}$8, whereas the irrotational value remains positive throughout the allowed range. The relative difference between the two prescriptions, $\mathcal{K}$9, grows with $$ds^2 = -e^{\nu}dt^2 + e^{\lambda}dr^2 + r^2d\Omega^2$$0 even though the underlying equations reduce continuously to the isotropic limit as $$ds^2 = -e^{\nu}dt^2 + e^{\lambda}dr^2 + r^2d\Omega^2$$1 — confirming that the prescription dependence of relativistic magnetic Love numbers persists in anisotropic environments. Most importantly, since these deviations are $$ds^2 = -e^{\nu}dt^2 + e^{\lambda}dr^2 + r^2d\Omega^2$$2 with order-unity coefficients sensitive to $$ds^2 = -e^{\nu}dt^2 + e^{\lambda}dr^2 + r^2d\Omega^2$$3, tidal observables can carry environmental imprints comparable to or larger than those in the geodesic and ringdown sectors, making them the strongest candidate for detecting environmental effects in future GW observations.

Limitations and open questions

Several caveats bound the scope of these conclusions. The numerical analysis is restricted to the axial sector; polar perturbations, where matter fluctuations couple dynamically to the geometry, are left unexplored despite being expected to yield richer phenomenology relevant to EMRIs. The choice of closure prescription (strictly static versus irrotational) is shown not to affect leading-order QNMs, but its impact on the tidal sector — and more broadly the physical interpretation and microscopic justification of either prescription — remains open. The DEC violations near the horizon are tolerated on the assumption that they coincide with the breakdown region of the effective-fluid approximation; self-consistent models incorporating explicit DM spike profiles would be needed to verify this. Finally, all quantitative results assume nonrotating BHs and small compactness $$ds^2 = -e^{\nu}dt^2 + e^{\lambda}dr^2 + r^2d\Omega^2$$4; rotation, time-domain evolutions, and realistic spike profiles are identified as necessary extensions before direct observational application.

Conclusion

This work provides a unified, self-consistent relativistic treatment of BHs embedded in anisotropic fluids with nonzero radial pressure, generalizing the Einstein cluster paradigm and recovering all previously studied limits. Its principal findings are a robust universality of the axial QNM response — governed by a redshift $$ds^2 = -e^{\nu}dt^2 + e^{\lambda}dr^2 + r^2d\Omega^2$$5 regardless of closure prescription — and a markedly stronger sensitivity of magnetic tidal Love numbers, which can vanish or change sign at moderate anisotropy. The results establish tidal observables as the most promising channel for probing environmental effects around black holes, while delineating precisely which questions — polar dynamics, rotation, and fluid microphysics — remain to be addressed.