Kerr-MOG Black Holes in Modified Gravity

• Kerr-MOG black holes are defined by three parameters (M, a, α) that modify the traditional Kerr metric by introducing an effective gravitational charge.
• They exhibit unique shadow morphology, thermodynamics, and geodesic structures that enhance gravitational lensing and energy extraction compared to Kerr black holes.
• Observational methods such as EHT and polarimetric imaging offer actionable tests for constraining the modified gravity parameter α in strong-field regimes.

A Kerr-MOG black hole is the stationary, axisymmetric, vacuum solution of Scalar-Tensor-Vector Gravity (STVG, commonly called MOG) representing a rotating, uncharged compact object with three key parameters: ADM mass $M$, spin parameter $a$, and the dimensionless modified gravity parameter $\alpha$. This solution generalizes the classic Kerr spacetime of general relativity by replacing the gravitational coupling $G_N$ with $G = G_N (1+\alpha)$ and mimicking the inclusion of an effective gravitational charge. The Kerr-MOG metric has profound implications for black hole thermodynamics, causal structure, geodesics, black hole shadows, electromagnetic and gravitational lensing, energy extraction mechanisms, and strong-field tests of gravity.

1. Kerr-MOG Metric Structure and Global Properties

The Kerr-MOG metric in Boyer–Lindquist coordinates $(t, r, \theta, \phi)$ is

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

where

$$\rho^2 = r^2 + a^2\cos^2\theta, \qquad \Delta(r) = r^2 - 2 G_N (1+\alpha) M r + a^2 + \alpha G_N^2 (1+\alpha) M^2.$$

The ADM mass and angular momentum are $M_{\rm ADM} = M(1+\alpha)$ and $J = M(1+\alpha)a$ (Sheoran et al., 2017). The metric reduces to Kerr for $a$0.

Horizons and extremality. The roots of $a$1 yield the Cauchy ($a$2) and event ($a$3) horizons: $a$4 A regular black hole requires $a$5, with extremality at equality.

Ergosurface and singularities. The ergosphere is the locus $a$6, which yields a deformed surface compared to Kerr due to $a$7. The metric possesses a ring singularity at $a$8 and possible closed timelike curve regions for "superspinar" over-spinning cases (Pradhan et al., 2024, Moffat, 2014).

Causal and topological features. In the zero-mass limit, the spacetime becomes a traversable wormhole with a ring throat, paralleling Gibbons–Volkov for Kerr (Pradhan et al., 2024).

2. Geodesic Structure, Photon Orbits, and Shadows

Geodesics and integrability. The Hamilton–Jacobi equation separates as in Kerr, enabling definition of conserved energy $a$9, axial angular momentum $\alpha$0, and Carter constant $\alpha$1 (Moffat, 2014, Li et al., 2024). The effect of $\alpha$2 appears in the form of the effective gravitational charge and in the horizon function.

Photon spheres and spherical photon orbits. Spherical photon orbits satisfy $\alpha$3 and $\alpha$4, where $\alpha$5 is the radial geodesic potential. In Kerr-MOG, all spherical photon orbits outside the horizon are radially unstable. The location and number of photon orbits depend sensitively on $\alpha$6 and $\alpha$7. The allowed spin range is reduced: $\alpha$8 (Li et al., 2024).

Black-hole shadow and observational signatures. The boundary of the shadow corresponds to the projection of critical photon orbits onto a distant observer’s sky. For equatorial observers,

$\alpha$9

where $G_N$0 and $G_N$1 for spherical photon orbits. The shadow diameter (for near-extremal spin) is (Guo et al., 2018)

$G_N$2

The shadow shrinks monotonically as $G_N$3 increases, in contrast with the static (Schwarzschild-like) MOG black hole, for which the shadow grows with $G_N$4 (Moffat, 2015, Guo et al., 2018). This nonmonotonic effect in the rotating case is a distinctive feature of the Kerr-MOG geometry.

Energy emission and lensing in plasma. In plasma, the shadow is modified; inhomogeneous plasma models tighter constrain $G_N$5 with current EHT observations (Yasmin et al., 25 Feb 2025). The deflection angle in vacuum acquires an overall factor $G_N$6, enhancing lensing relative to Kerr (Övgün et al., 2018).

3. Dynamical and Thermodynamical Properties

Surface gravity, temperature, and entropy.

The angular velocity of the event horizon is

$G_N$7

The surface gravity and Hawking temperature are

$G_N$8

The Bekenstein–Hawking entropy generalizes as

$G_N$9

All thermodynamical quantities smoothly reduce to Kerr (for $G = G_N (1+\alpha)$0), with novel $G = G_N (1+\alpha)$1-dependent corrections (Mureika et al., 2015, Pradhan, 2017).

Stability and remnants. Kerr-MOG black holes feature new thermodynamic structure: a small, cold extremal remnant at $G = G_N (1+\alpha)$2 is guaranteed for any $G = G_N (1+\alpha)$3, and there exists a branch of locally stable black holes with positive heat capacity ($G = G_N (1+\alpha)$4), unlike Kerr (Mureika et al., 2015).

First law and Smarr formula. The first law and Smarr relation incorporate a gravitational "charge" term, with the identifications

$G = G_N (1+\alpha)$5

where $G = G_N (1+\alpha)$6 plays an algebraic role analogous to electric charge (Mureika et al., 2015, Pradhan, 2017).

Kerr-MOG/CFT correspondence. The extremal Kerr-MOG black hole satisfies a holographic duality with central charge $G = G_N (1+\alpha)$7 and left-moving temperature $G = G_N (1+\alpha)$8, with Cardy entropy matching the Bekenstein–Hawking value (Pradhan, 2017).

4. Astrophysical and Observational Signatures

Shadows and EHT constraints. The "shadow diameter" measured by the Event Horizon Telescope places bounds on the pair $G = G_N (1+\alpha)$9 for sources such as M87* and Sgr A*. Current data favor $(t, r, \theta, \phi)$0 for moderate spin in M87*, but acceptable spin intervals are broadened relative to Kerr (Zhang et al., 2024, Wang et al., 12 Nov 2025).

Polarimetric imaging. The orientation and spatial pattern of linear polarization near the photon ring carry clear imprints of both $(t, r, \theta, \phi)$1 and $(t, r, \theta, \phi)$2. The radius at which the characteristic polarization pattern "twists" by $(t, r, \theta, \phi)$3 increases with $(t, r, \theta, \phi)$4 and can be measured with next-generation mm-VLBI (Wang et al., 12 Nov 2025).

Strong-field geodesic observables.

• The ISCO and ISSO radii increase with $(t, r, \theta, \phi)$5.
• Nodal and periastron precession frequencies are amplified; nodal precession (Lense–Thirring) increases monotonically with both $(t, r, \theta, \phi)$6 and $(t, r, \theta, \phi)$7, periastron advance grows with $(t, r, \theta, \phi)$8 but can be suppressed by spin (Wang et al., 4 Jul 2025).
• Redshift–blueshift spectroscopy of stars on close orbits provides a direct method to distinguish Kerr-MOG from Kerr: the range of allowed spectral shifts increases with $(t, r, \theta, \phi)$9, providing future tests of strong gravity near Sgr A* (Sheoran et al., 2017).

Gravitational lensing. Weak-field lensing is enhanced by a scaling factor $$ds^2 = -\frac{\Delta}{\rho^2}(dt - a \sin^2\theta d\phi)^2 + \frac{\sin^2\theta}{\rho^2}\big[(r^2 + a^2)d\phi - a dt\big]^2 + \frac{\rho^2}{\Delta}dr^2 + \rho^2 d\theta^2,$$0 relative to Kerr, affecting angular separations and time delays of lensed images (Övgün et al., 2018).

Causal and frame-dragging signatures, lensing-induced ring deformations, and time-variability (period of hot-spot recurrence) all offer independent measures sensitive to $$ds^2 = -\frac{\Delta}{\rho^2}(dt - a \sin^2\theta d\phi)^2 + \frac{\sin^2\theta}{\rho^2}\big[(r^2 + a^2)d\phi - a dt\big]^2 + \frac{\rho^2}{\Delta}dr^2 + \rho^2 d\theta^2,$$1 (Guo et al., 2018, Moffat, 2015).

5. Particle, Field, and Energy Extraction Dynamics

Charged-particle and plasma dynamics. The dynamics of charged particles around Kerr-MOG black holes, especially in the presence of magnetic fields, is influenced by the combination of mechanical parameters (energy $$ds^2 = -\frac{\Delta}{\rho^2}(dt - a \sin^2\theta d\phi)^2 + \frac{\sin^2\theta}{\rho^2}\big[(r^2 + a^2)d\phi - a dt\big]^2 + \frac{\rho^2}{\Delta}dr^2 + \rho^2 d\theta^2,$$2, angular momentum $$ds^2 = -\frac{\Delta}{\rho^2}(dt - a \sin^2\theta d\phi)^2 + \frac{\sin^2\theta}{\rho^2}\big[(r^2 + a^2)d\phi - a dt\big]^2 + \frac{\rho^2}{\Delta}dr^2 + \rho^2 d\theta^2,$$3, spin $$ds^2 = -\frac{\Delta}{\rho^2}(dt - a \sin^2\theta d\phi)^2 + \frac{\sin^2\theta}{\rho^2}\big[(r^2 + a^2)d\phi - a dt\big]^2 + \frac{\rho^2}{\Delta}dr^2 + \rho^2 d\theta^2,$$4), the magnetic field coupling, and $$ds^2 = -\frac{\Delta}{\rho^2}(dt - a \sin^2\theta d\phi)^2 + \frac{\sin^2\theta}{\rho^2}\big[(r^2 + a^2)d\phi - a dt\big]^2 + \frac{\rho^2}{\Delta}dr^2 + \rho^2 d\theta^2,$$5. Increasing $$ds^2 = -\frac{\Delta}{\rho^2}(dt - a \sin^2\theta d\phi)^2 + \frac{\sin^2\theta}{\rho^2}\big[(r^2 + a^2)d\phi - a dt\big]^2 + \frac{\rho^2}{\Delta}dr^2 + \rho^2 d\theta^2,$$6 tends to enlarge chaotic regions in phase space but is generally subdominant compared to $$ds^2 = -\frac{\Delta}{\rho^2}(dt - a \sin^2\theta d\phi)^2 + \frac{\sin^2\theta}{\rho^2}\big[(r^2 + a^2)d\phi - a dt\big]^2 + \frac{\rho^2}{\Delta}dr^2 + \rho^2 d\theta^2,$$7 or $$ds^2 = -\frac{\Delta}{\rho^2}(dt - a \sin^2\theta d\phi)^2 + \frac{\sin^2\theta}{\rho^2}\big[(r^2 + a^2)d\phi - a dt\big]^2 + \frac{\rho^2}{\Delta}dr^2 + \rho^2 d\theta^2,$$8 (Xu et al., 2024). Magnetic Penrose process efficiencies and the structure of stable circular orbits are both modulated by $$ds^2 = -\frac{\Delta}{\rho^2}(dt - a \sin^2\theta d\phi)^2 + \frac{\sin^2\theta}{\rho^2}\big[(r^2 + a^2)d\phi - a dt\big]^2 + \frac{\rho^2}{\Delta}dr^2 + \rho^2 d\theta^2,$$9 and can lead to particle trapping phenomena not present in pure Kerr (Khan et al., 2023).

Superradiance. Scalar field superradiant amplification is suppressed in Kerr-MOG relative to Kerr: both the critical frequency for amplification and the maximum gain are lower. The entire amplification band is "redshifted" for fixed spin, and growth rates for superradiant instabilities are reduced. This has implications for constraints based on black hole spin-down rates and GW echoes (Wondrak et al., 2018).

Energy extraction: Comisso–Asenjo Process (CAP) and Blandford–Znajek. The MOG deformation amplifies magnetic reconnection-driven energy extraction power and efficiency, making the CAP more effective relative to Blandford–Znajek for fixed spin. Even a small $$\rho^2 = r^2 + a^2\cos^2\theta, \qquad \Delta(r) = r^2 - 2 G_N (1+\alpha) M r + a^2 + \alpha G_N^2 (1+\alpha) M^2.$$0 (current EHT upper bound for Sgr A*) can boost CAP modestly relative to GR (Khodadi et al., 2023).

Cosmic censorship, overspinning, and stability.

• Test particle and test-field analyses show that for extremal and near-extremal Kerr-MOG black holes, finely tuned absorption events can drive the solution beyond the extremal bound, violating weak cosmic censorship (WCCC) (Liang et al., 2018, Düztaş, 2019).
• However, when absorption is modeled as a series of infinitesimal, adiabatic steps, or when the horizon’s response is properly taken into account, the horizon is preserved and WCCC is restored at $$\rho^2 = r^2 + a^2\cos^2\theta, \qquad \Delta(r) = r^2 - 2 G_N (1+\alpha) M r + a^2 + \alpha G_N^2 (1+\alpha) M^2.$$1 (Liang et al., 2018, Düztaş, 2019).

6. Methodologies and Numerical Techniques

Advanced symplectic integration algorithms (e.g., PRK$$\rho^2 = r^2 + a^2\cos^2\theta, \qquad \Delta(r) = r^2 - 2 G_N (1+\alpha) M r + a^2 + \alpha G_N^2 (1+\alpha) M^2.$$2) are essential for accurately simulating chaotic charged-particle motion, ensuring conservation of Hamiltonian constraints and distinguishing physical chaos from numerical artifacts (Xu et al., 2024). Poincaré section analyses, Fast Lyapunov Indicators (FLI), and systematic parameter scans over $$\rho^2 = r^2 + a^2\cos^2\theta, \qquad \Delta(r) = r^2 - 2 G_N (1+\alpha) M r + a^2 + \alpha G_N^2 (1+\alpha) M^2.$$3 clarify dynamical transitions in orbit structure. Ray-tracing based on the full Kerr-MOG metric and consistent radiative transfer codes (e.g., symphony, ARCMANCER) are employed to calculate synthetic images and polarization morphologies (Zhang et al., 2024, Wang et al., 12 Nov 2025).

7. Outlook: Distinguishability, Constraints, and Theoretical Implications

Kerr-MOG black holes inherit much of the integrable, algebraically-special structure of Kerr but are characterized by reduced extremal spin, modified precessional dynamics, distinctive shadow morphology, and altered superradiant and lensing properties. Observational constraints from EHT, magnetized disk modeling, polarization, and orbit tracking all favor $$\rho^2 = r^2 + a^2\cos^2\theta, \qquad \Delta(r) = r^2 - 2 G_N (1+\alpha) M r + a^2 + \alpha G_N^2 (1+\alpha) M^2.$$4–$$\rho^2 = r^2 + a^2\cos^2\theta, \qquad \Delta(r) = r^2 - 2 G_N (1+\alpha) M r + a^2 + \alpha G_N^2 (1+\alpha) M^2.$$5 for stellar and supermassive black holes. The combination of polarimetric imaging, high-resolution timing, and spectroscopic techniques offers a multi-modal route to empirically distinguish Kerr-MOG from Kerr in the next generation of strong-field experiments (Zhang et al., 2024, Wang et al., 12 Nov 2025, Sheoran et al., 2017).

In sum, Kerr-MOG spacetimes provide a working test-bed for strong-field extensions of GR and for the phenomenology of black holes in modified gravity, preserving separability and horizon structure while encoding distinctive, parameterizable departures in both mathematical and observational domains.