Horndeski Rotating Black Hole

• Rotating Horndeski black holes are exact solutions in scalar-tensor gravity, extending the Kerr metric by incorporating a scalar hair parameter.
• The additional scalar hair alters the horizon structure, shadow geometry, and quasinormal mode spectra, offering observable differences from general relativity.
• Observational probes such as EHT imaging, QPO analysis, and polarimetric studies provide practical constraints on the scalar hair and test alternative gravity models.

A rotating Horndeski black hole is an exact or numerically constructed solution to the field equations of Horndeski gravity—the most general 4D scalar-tensor theory with a single scalar degree of freedom and second-order equations of motion. These objects generalize the Kerr black hole by including an additional parameter ("scalar hair") characterizing deviations from general relativity (GR) due to nonminimal scalar–curvature couplings or scalar field profiles. Rotating Horndeski black holes have unique phenomenology in their horizon structure, spacetime geometry, electromagnetic observables (e.g., shadow, accretion-image morphology), energetics, and linear dynamics (quasinormal modes, superradiance, instabilities) that can be directly linked to the underlying scalar–tensor modifications (Heydari-Fard et al., 24 Oct 2025, Lei et al., 2023, Afrin et al., 2021).

1. Metric Structure and Hair Parameterization

The canonical rotating Horndeski black hole metric stems from a Newman–Janis complexification of a static, spherically symmetric solution with scalar hair. In Boyer–Lindquist type coordinates $(t,r,\theta,\varphi)$, the general line element is:

$$ds^2 = -\frac{\Delta - a^2\sin^2\theta}{\Sigma}\,dt^2 + \frac{\Sigma}{\Delta}\,dr^2 + \Sigma\,d\theta^2 + \frac{2a\sin^2\theta[\Delta-(r^2+a^2)]}{\Sigma}dt\,d\varphi + \frac{\sin^2\theta[(r^2+a^2)^2-\Delta a^2\sin^2\theta]}{\Sigma}\,d\varphi^2$$

where

$$\Sigma = r^2 + a^2\cos^2\theta, \qquad \Delta = r^2 + a^2 - 2M r + h r \ln\left(\frac{r}{2M}\right)$$

$M$ is the mass, $a$ the spin, and $h$ the dimensionful scalar hair parameter. For $h = 0$, this metric reduces precisely to the Kerr solution. The function $\Delta$ encodes the scalar hair modification, derived from the backreaction of a nontrivial scalar field allowed by the shift-symmetric quartic sector of Horndeski theory. The hair affects the horizon structure, the ergosphere, and frame-dragging properties, with parameters constrained by regularity and physical requirements (Heydari-Fard et al., 24 Oct 2025, Afrin et al., 2021, Heydari-Fard et al., 2023).

A minority of constructions (e.g., trace-anomaly/stealth solutions, disformal Kerr) introduce further generalizations: noncircularity, explicit $\theta$-dependence of horizon location, or algebraically general Petrov type I, but the aforementioned class contains all current phenomenologically relevant axisymmetric solutions (Fernandes, 2023, Achour et al., 22 Dec 2025).

2. Horizon Structure and Geometric Properties

The event horizon is the largest positive root of $\Delta(r_H) = 0$: $$ds^2 = -\frac{\Delta - a^2\sin^2\theta}{\Sigma}\,dt^2 + \frac{\Sigma}{\Delta}\,dr^2 + \Sigma\,d\theta^2 + \frac{2a\sin^2\theta[\Delta-(r^2+a^2)]}{\Sigma}dt\,d\varphi + \frac{\sin^2\theta[(r^2+a^2)^2-\Delta a^2\sin^2\theta]}{\Sigma}\,d\varphi^2$$0 For $$ds^2 = -\frac{\Delta - a^2\sin^2\theta}{\Sigma}\,dt^2 + \frac{\Sigma}{\Delta}\,dr^2 + \Sigma\,d\theta^2 + \frac{2a\sin^2\theta[\Delta-(r^2+a^2)]}{\Sigma}dt\,d\varphi + \frac{\sin^2\theta[(r^2+a^2)^2-\Delta a^2\sin^2\theta]}{\Sigma}\,d\varphi^2$$1 and $$ds^2 = -\frac{\Delta - a^2\sin^2\theta}{\Sigma}\,dt^2 + \frac{\Sigma}{\Delta}\,dr^2 + \Sigma\,d\theta^2 + \frac{2a\sin^2\theta[\Delta-(r^2+a^2)]}{\Sigma}dt\,d\varphi + \frac{\sin^2\theta[(r^2+a^2)^2-\Delta a^2\sin^2\theta]}{\Sigma}\,d\varphi^2$$2, two horizons exist; for $$ds^2 = -\frac{\Delta - a^2\sin^2\theta}{\Sigma}\,dt^2 + \frac{\Sigma}{\Delta}\,dr^2 + \Sigma\,d\theta^2 + \frac{2a\sin^2\theta[\Delta-(r^2+a^2)]}{\Sigma}dt\,d\varphi + \frac{\sin^2\theta[(r^2+a^2)^2-\Delta a^2\sin^2\theta]}{\Sigma}\,d\varphi^2$$3, only a single horizon appears at $$ds^2 = -\frac{\Delta - a^2\sin^2\theta}{\Sigma}\,dt^2 + \frac{\Sigma}{\Delta}\,dr^2 + \Sigma\,d\theta^2 + \frac{2a\sin^2\theta[\Delta-(r^2+a^2)]}{\Sigma}dt\,d\varphi + \frac{\sin^2\theta[(r^2+a^2)^2-\Delta a^2\sin^2\theta]}{\Sigma}\,d\varphi^2$$4.

The static limit (ergosurface) satisfies $$ds^2 = -\frac{\Delta - a^2\sin^2\theta}{\Sigma}\,dt^2 + \frac{\Sigma}{\Delta}\,dr^2 + \Sigma\,d\theta^2 + \frac{2a\sin^2\theta[\Delta-(r^2+a^2)]}{\Sigma}dt\,d\varphi + \frac{\sin^2\theta[(r^2+a^2)^2-\Delta a^2\sin^2\theta]}{\Sigma}\,d\varphi^2$$5, leading to $$ds^2 = -\frac{\Delta - a^2\sin^2\theta}{\Sigma}\,dt^2 + \frac{\Sigma}{\Delta}\,dr^2 + \Sigma\,d\theta^2 + \frac{2a\sin^2\theta[\Delta-(r^2+a^2)]}{\Sigma}dt\,d\varphi + \frac{\sin^2\theta[(r^2+a^2)^2-\Delta a^2\sin^2\theta]}{\Sigma}\,d\varphi^2$$6. The horizon's angular velocity is: $$ds^2 = -\frac{\Delta - a^2\sin^2\theta}{\Sigma}\,dt^2 + \frac{\Sigma}{\Delta}\,dr^2 + \Sigma\,d\theta^2 + \frac{2a\sin^2\theta[\Delta-(r^2+a^2)]}{\Sigma}dt\,d\varphi + \frac{\sin^2\theta[(r^2+a^2)^2-\Delta a^2\sin^2\theta]}{\Sigma}\,d\varphi^2$$7 and the surface gravity is given by: $$ds^2 = -\frac{\Delta - a^2\sin^2\theta}{\Sigma}\,dt^2 + \frac{\Sigma}{\Delta}\,dr^2 + \Sigma\,d\theta^2 + \frac{2a\sin^2\theta[\Delta-(r^2+a^2)]}{\Sigma}dt\,d\varphi + \frac{\sin^2\theta[(r^2+a^2)^2-\Delta a^2\sin^2\theta]}{\Sigma}\,d\varphi^2$$8 demonstrating explicit dependence on the scalar hair (Heydari-Fard et al., 2023, Donmez, 2024).

In more exotic cases (e.g., disformal Kerr), the horizon and ergosphere structures remain identical to Kerr, and causal pathologies (closed timelike curves) are avoided, but the Petrov classification shifts to type I, and off-equatorial geodesics acquire corrections (Achour et al., 22 Dec 2025).

3. Photon Dynamics, Black Hole Shadows, and Plasma Effects

Null geodesics in the rotating Horndeski background can be separated in the Hamilton–Jacobi framework, allowing direct calculation of photon regions and black hole shadow shapes. For the spacetime with (possibly nonuniform) plasma, the full Hamiltonian becomes: $$ds^2 = -\frac{\Delta - a^2\sin^2\theta}{\Sigma}\,dt^2 + \frac{\Sigma}{\Delta}\,dr^2 + \Sigma\,d\theta^2 + \frac{2a\sin^2\theta[\Delta-(r^2+a^2)]}{\Sigma}dt\,d\varphi + \frac{\sin^2\theta[(r^2+a^2)^2-\Delta a^2\sin^2\theta]}{\Sigma}\,d\varphi^2$$9 with separability for physically motivated plasma profiles (uniform, radial, latitudinal). The celestial coordinates of the shadow boundary depend on the normalized impact parameters $$\Sigma = r^2 + a^2\cos^2\theta, \qquad \Delta = r^2 + a^2 - 2M r + h r \ln\left(\frac{r}{2M}\right)$$0, evaluated at the photon region, and the shadow outline satisfies (Heydari-Fard et al., 24 Oct 2025): $$\Sigma = r^2 + a^2\cos^2\theta, \qquad \Delta = r^2 + a^2 - 2M r + h r \ln\left(\frac{r}{2M}\right)$$1 for a uniform plasma of density parameter $$\Sigma = r^2 + a^2\cos^2\theta, \qquad \Delta = r^2 + a^2 - 2M r + h r \ln\left(\frac{r}{2M}\right)$$2.

Increasing $$\Sigma = r^2 + a^2\cos^2\theta, \qquad \Delta = r^2 + a^2 - 2M r + h r \ln\left(\frac{r}{2M}\right)$$3 at fixed $$\Sigma = r^2 + a^2\cos^2\theta, \qquad \Delta = r^2 + a^2 - 2M r + h r \ln\left(\frac{r}{2M}\right)$$4 increases the shadow size and distortion (for instance, with $$\Sigma = r^2 + a^2\cos^2\theta, \qquad \Delta = r^2 + a^2 - 2M r + h r \ln\left(\frac{r}{2M}\right)$$5, $$\Sigma = r^2 + a^2\cos^2\theta, \qquad \Delta = r^2 + a^2 - 2M r + h r \ln\left(\frac{r}{2M}\right)$$6 changing from 0 to $$\Sigma = r^2 + a^2\cos^2\theta, \qquad \Delta = r^2 + a^2 - 2M r + h r \ln\left(\frac{r}{2M}\right)$$7 increases the horizontal shadow diameter by $$\Sigma = r^2 + a^2\cos^2\theta, \qquad \Delta = r^2 + a^2 - 2M r + h r \ln\left(\frac{r}{2M}\right)$$8). Inhomogeneous plasma profiles (radial, latitudinal) shrink the shadow until it vanishes beyond a critical plasma density. The shadow position is shifted by spin-induced frame dragging, and the shadow is further distorted with increasing $$\Sigma = r^2 + a^2\cos^2\theta, \qquad \Delta = r^2 + a^2 - 2M r + h r \ln\left(\frac{r}{2M}\right)$$9 and/or $M$0 (Heydari-Fard et al., 24 Oct 2025, Afrin et al., 2021).

The Event Horizon Telescope observations of M87* shadow ($M$1as; $M$2; $M$3 Mpc) set direct observational upper limits on the allowed plasma densities and scalar hair parameters (Heydari-Fard et al., 24 Oct 2025, Afrin et al., 2021). Bias analysis demonstrates that for $M$4 even next-generation VLBI arrays will not distinguish shadows from Kerr, but larger hair ($M$5) will be testable (Afrin et al., 2021).

4. Quasinormal Modes, Superradiance, and Instabilities

The linear dynamics of massive scalar perturbations in the rotating Horndeski background obey a separable wave equation, yielding both quasinormal modes (QNMs) and quasibound states (QBSs). The effective radial equation (in tortoise coordinates) features a potential well deepened by the hair parameter $M$6, and the standard boundary conditions (ingoing at the horizon, outgoing or exponentially decaying at infinity) lead to:

• For massless scalar perturbations, QNM frequencies have real parts suppressed and imaginary parts reduced (longer lifetimes) as $M$7 becomes more negative. No instability is found.
• For massive perturbations, scalar clouds and superradiant instability can develop (via the standard critical condition $M$8), but the onset threshold, and width of the instability region in $M$9 space, are modified relative to Kerr. Negative $a$0 enables superradiant instability at lower spins and for slightly higher scalar masses.

The matrix method—a numerical eigensolver applied to the ODEs—yields spectra, revealing that increasing negative hair relaxes the ultra-spinning threshold for superradiant instability compared to pure Kerr (Lei et al., 2023, Jha et al., 2022). However, the ultimate instability window for massive scalar fields remains controlled primarily by the horizon angular velocity, so the qualitative criteria are preserved (Jha et al., 2022).

5. Accretion Disk Imaging and X-ray/Polarimetric Probes

Physical observables beyond the shadow include the electromagnetic and polarimetric imaging of thin and thick accretion disks, and timing features from orbiting plasma.

• Thin disk models: In Novikov–Thorne disks on rotating Horndeski backgrounds, the ISCO radius decreases monotonically with increasing $a$1, supporting higher disk efficiency and hotter, more luminous spectra. The central shadow shrinks, but the lensed photon ring brightens, and relativistic Doppler and lensing effects are modified by the hair (Heydari-Fard et al., 2023).
• Thick disk/radiative transfer: For magnetized, conical thick disks, the strongest hair-dependence is found in the photon ring diameter ($a$2), which can be shifted by $a$3as for $a$4 varying from 0 to the critical value at extremality (for $a$5), exceeding the angular resolution of next-generation space-VLBI (BHEX). The total observed flux and brightness of the ring decrease as $a$6 increases, primarily due to strengthened gravitational redshift. The ring diameter is relatively insensitive to accretion-flow microphysics, so it serves as a robust hair diagnostic (Wan et al., 30 Nov 2025).
• Near-horizon polarization: The electric vector position angle (EVPA) of the photon ring is governed dominantly by the geometric properties of the spacetime and shows small but potentially measurable dependence on the hair. For slowly rotating black holes, increasing $a$7 raises the near-horizon EVPA, while for high spins, the trend reverses. This effect is more pronounced at low observer inclinations and specific azimuths (Chen et al., 16 Sep 2025).
• Quasiperiodic oscillations (QPOs): Epicyclic frequencies (orbital, radial, vertical) in the equatorial plane are shifted by the hair, yielding tight empirical constraints when fit to observed QPO triplets from X-ray binaries. Analyses combining MCMC techniques and QPO frequencies produce bounds of $a$8 to $a$9, much tighter than those from previous accretion simulations (Wu et al., 19 Aug 2025).
• Bondi-Hoyle accretion: In Bondi–Hoyle accretion, large negative hair moves the stagnation point of the accretion shock cone towards the horizon and can suppress or destroy the cone entirely for $h$0. QPO modes in the cone are also modified as hair increases (Donmez, 2024).

6. Thermodynamics and Conserved Quantities

Thermodynamic quantities retain their Kerr-like form but depend on the hair parameter:

$h$1

Komar charges and the first law acquire additional scalar-hair terms. The Smarr formula and most thermodynamic relations survive up to the replacement of $h$2 evaluated at the outer horizon, with corrections proportional to $h$3 (Walia et al., 2021).

7. Observational Signatures and Constraints

Astrophysical observations, notably the EHT shadow of M87*, have been used to directly constrain the allowed parameter space of rotating Horndeski black holes. For the observed shadow diameter $h$4as and ring circularity deviation $h$5, viable solutions must satisfy

$h$6

for nearly edge-on inclination (Afrin et al., 2021). Systematic bias analyses show that for $h$7, current EHT sensitivity cannot distinguish Horndeski from Kerr, but higher-precision VLBI or polarimetric measurements could probe larger hair. Weak-field lensing (e.g., Sgr A*-S2, $h$8) offers independent constraints, implying $h$9 for most SMBH systems (Walia et al., 2021).

Future X-ray timing (QPOs), high-frequency polarimetry, and horizon-scale imaging (with BHEX/ngEHT or space-VLBI) are projected to test the existence of nonzero scalar hair at the $h = 0$0 level, potentially ruling out large classes of Horndeski-like deviations from general relativity (Wan et al., 30 Nov 2025, Chen et al., 16 Sep 2025, Wu et al., 19 Aug 2025).

---

References:

• (Heydari-Fard et al., 24 Oct 2025) Influence of plasma on the observational appearance of rotating black holes in Horndeski gravity
• (Lei et al., 2023) Scalar field perturbation around a rotating hairy black hole: quasinormal modes, quasibound states and superradiant instability
• (Afrin et al., 2021) Testing Horndeski Gravity from EHT Observational Results for Rotating Black Holes
• (Heydari-Fard et al., 2023) Thin accretion disk images of rotating hairy Horndeski black holes
• (Jha et al., 2022) Superradiant energy extraction from rotating hairy Horndeski black holes
• (Chen et al., 16 Sep 2025) Near-horizon polarized images of rotating hairy Horndeski black holes
• (Wu et al., 19 Aug 2025) Parameter constraints on Horndeski rotating black hole through quasiperiodic oscillations
• (Donmez, 2024) Bondi-Hoyle-Lyttleton Accretion around the Rotating Hairy Horndeski Black Hole
• (Wan et al., 30 Nov 2025) Probing the Scalar Hair of Rotating Horndeski Black Holes through Thick Disk Images
• (Achour et al., 22 Dec 2025) Circular Disformal Kerr: An Exact Rotating Black Hole Beyond GR
• (Walia et al., 2021) Rotating Black Holes in Horndeski Gravity: Thermodynamic and Gravitational Lensing