Observable Optical Signatures, Particle Dynamics and Epicyclic Frequencies of Mod(A)Max Black Holes

In this work, we investigate the observable optical signatures of the Mod(A)Max black hole spacetime. We analyze key optical features, including the photon sphere, black hole shadow, and photon trajectories, and examine how these observables depend on the underlying geometric parameters, such as the electric charge and the Mod(A)Max coupling parameter. We further study the dynamics of neutral test particles in the vicinity of the black hole by deriving the effective potential within the Hamiltonian formalism. Using this potential, we obtain the specific energy and specific angular momentum for test particles on circular orbits of fixed radius, as well as the innermost stable circular orbit (ISCO), and explore how the geometric parameters influence these quantities and the ISCO radius. Finally, we derive the epicyclic (azimuthal, radial, and vertical) frequencies to analyze quasi-periodic oscillations (QPOs) exploring how the geometric parameters influences these and discuss their physical implications.

💡 Research Summary

This paper investigates the observable optical signatures and particle dynamics of black holes arising from the Modified Maxwell (ModMax) and its phantom counterpart, the Mod(A)Max, which are nonlinear extensions of classical electrodynamics. The authors start from the spherically symmetric line element
ds² = –f(r) dt² + dr²/f(r) + r²(dθ² + sin²θ dφ²)
with the lapse function f(r) = 1 – 2M/r + η e^{–γ}Q²/r², where M is the mass, Q the electric charge, γ the ModMax coupling parameter, and η = +1 for the ordinary ModMax solution and η = –1 for the phantom Mod(A)Max case.

Photon sphere and shadow: By constructing the Lagrangian for null geodesics, the conserved energy E and angular momentum L lead to an effective potential V_eff = L²/(r² f(r)). Circular photon orbits satisfy E² = V_eff and dV_eff/dr = 0, yielding the photon‑sphere radius
r_ph = (3M ± √(9M² – 8η e^{–γ}Q²))/2.
For η = +1 the charge term reduces the photon‑sphere size, while for η = –1 it enhances the gravitational pull, expanding the photon sphere. The shadow radius observed at infinity is r_sh = r_ph/√{f(r_ph)}; consequently, increasing Q shrinks the shadow in the normal ModMax case but enlarges it in the phantom Mod(A)Max case. The coupling γ modulates these trends through the exponential factor e^{–γ}, suppressing the charge contribution as γ grows. These dependencies provide potentially observable signatures for Event Horizon Telescope (EHT) measurements.

Neutral test‑particle dynamics: Using the Hamiltonian formalism, the effective potential for massive particles is V_eff^p = (1 + L²/r²) f(r). Imposing dV_eff^p/dr = 0 gives analytic expressions for the specific energy E(r) and specific angular momentum L(r) of circular equatorial orbits. The innermost stable circular orbit (ISCO) is found by solving d²V_eff^p/dr² = 0. Numerical analysis shows that in the normal ModMax case the ISCO radius decreases with larger Q and γ, allowing particles to orbit closer to the horizon. In contrast, the phantom Mod(A)Max case exhibits an outward shift of the ISCO as Q grows, reflecting the attractive nature of the phantom charge term. This behavior directly impacts accretion‑disk inner edges and the emitted X‑ray spectra.

Epicyclic frequencies and QPOs: The azimuthal (Keplerian) frequency is
Ω_φ = √{M/r³ – η e^{–γ}Q²/r⁴}.
Radial and vertical epicyclic frequencies follow from second‑order perturbations around circular orbits:
Ω_r² = Ω_φ²(1 – 6M/r + 9η e^{–γ}Q²/r²),
Ω_θ² = Ω_φ²(1 – 3M/r + 2η e^{–γ}Q²/r²).
In the normal ModMax scenario, increasing Q or γ lowers Ω_r and Ω_θ, shifting high‑frequency quasi‑periodic oscillations (HF‑QPOs) to lower values. Conversely, in the phantom Mod(A)Max case the same increase raises the epicyclic frequencies, moving HF‑QPOs to higher frequencies. This distinct signature offers a way to discriminate between the two theories using timing observations of X‑ray binaries.

Thermodynamics: The Hawking temperature is T = f’(r_h)/(4π) = (1/(4π r_h))