Relativistic Images in Strong Gravity

• Relativistic images are multiply-lensed images produced as photons wind around compact objects, serving as diagnostics for strong-field gravitational lensing.
• They are modeled using null geodesics with logarithmically diverging deflection angles and exponential decay in magnification near the photon sphere.
• Observations of angular positions, time delays, and spectral distortions in these images provide practical tests for general relativity and alternative gravity theories.

Relativistic images are multiply-lensed images produced when light rays from a background source traverse null geodesics in a strong gravitational field, winding around a compact object—typically a black hole—by large deflection angles. These images emerge through strong-field gravitational lensing as sequences accumulating near the photon sphere (unstable photon circular orbit), manifesting as faint, highly demagnified images close to the boundary of the black hole shadow or observed in spectral lines via distinct spectral substructures. Relativistic images probe the fundamental structure of spacetime in the immediate vicinity of compact objects, encode the spacetime parameters (e.g., mass, spin, charge), and are key diagnostics for testing general relativity in the strong-field regime as well as alternative theories of gravity.

1. Formation of Relativistic Images: Lens Theory and Geodesic Structure

Relativistic images originate when photons, emitted by a source near a compact lens (black hole or similar object), traverse null geodesics with deflection angles $\hat\alpha > 2\pi$. As a photon approaches the photon sphere (e.g., $r_{ph}=3M$ for Schwarzschild), the bending angle diverges logarithmically in the impact parameter $u$, according to

$$\hat\alpha(u) = -\bar{a}\,\ln\!\left(\frac{u}{u_m}-1\right) + \bar{b} + O(u-u_m),$$

where $u_m$ is the critical impact parameter for the photon sphere, and $\bar{a}$, $\bar{b}$ are metric-dependent strong-deflection coefficients (0810.2109, Ghazale et al., 28 Nov 2025, Sahu et al., 2013). The corresponding observables (image positions $\theta_n$, magnifications $\mu_n$, time delays $\Delta t_{n,n+1}$) are therefore encoded in the geometry near this critical region.

For a Schwarzschild black hole, the $r_{ph}=3M$0-th relativistic image is found at

$r_{ph}=3M$1

with $r_{ph}=3M$2, $r_{ph}=3M$3 the observer-lens distance (0810.2109). Magnifications decrease exponentially with $r_{ph}=3M$4, $r_{ph}=3M$5, rendering these images intrinsically faint (Ghazale et al., 28 Nov 2025).

In axisymmetric (Kerr) or electrically charged (Reissner–Nordström) backgrounds, the photon sphere splits into prograde and retrograde orbits with order-dependent $r_{ph}=3M$6, $r_{ph}=3M$7, and deflection coefficients $r_{ph}=3M$8; relativistic images shift and their properties depend sensitively on the spin $r_{ph}=3M$9, charge $u$0, and observer inclination $u$1 (Ghazale et al., 28 Nov 2025, Aratore et al., 2024). For black hole binaries with multiple photon spheres, distinct image families appear both inside and outside the critical angle (Patil et al., 2016).

2. Observables: Positions, Magnifications, and Time Delays

The core observables of relativistic images include:

• Angular Position: For each winding order $u$2, the image position is exponentially close to the shadow rim: $u$3, where the sequence accumulates towards $u$4 as $u$5 (0810.2109, Ghazale et al., 28 Nov 2025).
• Magnification: $u$6, strongly suppressed relative to primary/secondary images, with exact scaling set by the metric coefficients at the photon sphere (Aratore et al., 2021, Ghazale et al., 28 Nov 2025).
• Time Delay: The interval between adjacent image arrival times is nearly constant for spacetimes with a photon sphere (e.g., Schwarzschild/weakly-naked singularities):

$u$7

where $u$8 is the critical impact parameter (Sahu et al., 2013, 0810.2109). For spacetimes without a photon sphere, the delays decrease monotonically with $u$9 (Sahu et al., 2013).

• Deflection Signs: Effective deflection angles on the primary vs. secondary image side can change sign, and their detailed dependencies are critical for accurate modeling (0810.2109).

These properties are robust under wide variations in source alignment and lens-observer geometry; the universality of the position/delay formulas makes relativistic images powerful tools for precision black hole mass and distance measurements (0810.2109, Sahu et al., 2013).

3. Metric Sensitivity and Astrophysical Applications

Relativistic images serve as diagnostic probes of spacetime structure:

• Testing General Relativity and Beyond: The stepwise pattern in interferometric visibility functions reflects the sequence of images and is controlled by the metric's strong-deflection coefficients. The ratio of successive image fluxes or visibility step heights, $$\hat\alpha(u) = -\bar{a}\,\ln\!\left(\frac{u}{u_m}-1\right) + \bar{b} + O(u-u_m),$$0, provides a direct measurement of the Lyapunov exponent of the photon sphere orbit (Aratore et al., 2021).
• Spin and Inclination Diagnostics: In the Kerr metric, the relative positions and time delays of the primary and secondary relativistic images encode both spin $$\hat\alpha(u) = -\bar{a}\,\ln\!\left(\frac{u}{u_m}-1\right) + \bar{b} + O(u-u_m),$$1 and inclination $$\hat\alpha(u) = -\bar{a}\,\ln\!\left(\frac{u}{u_m}-1\right) + \bar{b} + O(u-u_m),$$2. Systematic inversion of measured $$\hat\alpha(u) = -\bar{a}\,\ln\!\left(\frac{u}{u_m}-1\right) + \bar{b} + O(u-u_m),$$3 yields unique $$\hat\alpha(u) = -\bar{a}\,\ln\!\left(\frac{u}{u_m}-1\right) + \bar{b} + O(u-u_m),$$4 pairs, strongly constraining the black hole parameters (Wang et al., 2022). Spin introduces order-dependent shifts in both the sky position and time-of-flight, reflecting frame-dragging (Aratore et al., 2024).
• Binary Lensing and Caustic Structure: Binary black holes with multiple photon spheres produce two infinite sequences of relativistic images—one outside and one inside the outer spherical radius. As the inter-black-hole separation increases, photon spheres can merge and vanish, switching the image morphology from infinite to finite with the emergence of novel radial caustics (Patil et al., 2016).

Astrophysical imaging/spectroscopy of the relativistic Fe K$$\hat\alpha(u) = -\bar{a}\,\ln\!\left(\frac{u}{u_m}-1\right) + \bar{b} + O(u-u_m),$$5 line profile harnesses these higher-order images. Although their integrated flux is $$\hat\alpha(u) = -\bar{a}\,\ln\!\left(\frac{u}{u_m}-1\right) + \bar{b} + O(u-u_m),$$6 of the total line, neglecting their contribution introduces systematic bias in inferred black hole spins (up to $$\hat\alpha(u) = -\bar{a}\,\ln\!\left(\frac{u}{u_m}-1\right) + \bar{b} + O(u-u_m),$$7; $$\hat\alpha(u) = -\bar{a}\,\ln\!\left(\frac{u}{u_m}-1\right) + \bar{b} + O(u-u_m),$$8 being the dimensionless spin) in spectroscopic analyses with next-generation X-ray missions (Falanga et al., 2021).

4. Relativistic Imaging in Alternate Contexts

Plasma-Modified Relativistic Images

The presence of plasma introduces frequency-dependent refractive effects, shifting the photon sphere outward, enlarging the angular radii of relativistic rings, and increasing their magnifications as the photon frequency approaches the plasma cutoff. The position and scaling formulas proceed with $$\hat\alpha(u) = -\bar{a}\,\ln\!\left(\frac{u}{u_m}-1\right) + \bar{b} + O(u-u_m),$$9, and all standard vacuum results are smoothly recovered as $u_m$0 (Tsupko et al., 2013).

Relativistic Image Doubling

In atmospheric Cherenkov imaging, relativistic air showers viewed at angles where the radial velocity component crosses subluminal values produce instantaneous image doubling: two images of the same shower, moving in opposite directions, appear in the focal plane with a sharp transition at a calculable altitude $u_m$1 (Nemiroff et al., 2019).

Quantum Field Theoretic Imaging

Quantum imaging protocols, such as ghost imaging with Unruh–DeWitt detectors in non-inertial frames, reveal that information about a quantum system can be reconstructed through relativistic detectors coupled to fields even in the absence of entanglement, with advantage over classical guessing quantified via signal contrast (Bornman et al., 2019).

5. Ray Tracing, Simulation Techniques, and Experimental Visualization

Numerical computation of relativistic images proceeds via backward ray tracing through the spacetime geometry, integrating the geodesic equations—either in 3+1 split from numerical simulations or in analytical spacetimes (typically Kerr or approximations to it) (Vincent et al., 2011, Vincent et al., 2012). Modern codes implement the full transfer of specific intensity, frequency shift, and radiative transfer to model emission from surfaces or accretion flows.

Recent laboratory experiments visualize the special-relativistic Terrell–Penrose effect in real time, producing images of Lorentz-contracted objects that reveal the "rotation" effect predicted by relativity for moving objects, distinct from Lorentz contraction itself. By time-gating fs-laser pulses, effective slow-light conditions enable experimental access to apparent shapes at relativistic velocities, confirming theoretical predictions (Hornof et al., 2024, Zhu, 12 Nov 2025). The mathematical framework for these visualizations employs the light-cone condition for apparent positions, Lorentz/Poincaré transformations, and explicit emission-time integrals to recover the observed "rotated" (not just contracted) shapes (Bajaj, 2021, Zhu, 12 Nov 2025).

6. Summary Table: Relativistic Image Properties in Key Geometries

Metric/Context	Deflection Scaling	Angular Position	Magnification	Time Delay	Key Dependence
Schwarzschild	$u_m$2	$u_m$3	$u_m$4	$u_m$5	$u_m$6, $u_m$7
Reissner–Nordström	$u_m$8, $u_m$9	similar, $\bar{a}$0 shifts	$\bar{a}$1 as $\bar{a}$2	$\bar{a}$3 as $\bar{a}$4 increases	$\bar{a}$5, $\bar{a}$6
Kerr (equatorial)	Separate pro/retrograde, $\bar{a}$7	$\bar{a}$8/$\bar{a}$9	$\bar{b}$0 with $\bar{b}$1	Shortened (pro), lengthened (retro)	$\bar{b}$2, $\bar{b}$3
Plasma-modified	$\bar{b}$4 as $\bar{b}$5	$\bar{b}$6	$\bar{b}$7 drastically	similar to vacuum	$\bar{b}$8, $\bar{b}$9
Binary BH (two photon spheres)	Logarithmic near each $\theta_n$0	Two sequences ($\theta_n$1)	Exponential fall by $\theta_n$2	Each $\theta_n$3	$\theta_n$4, inter-BH separation, $\theta_n$5

7. Observational Prospects and Theory Discrimination

Detection of relativistic images remains an extreme challenge: magnifications are typically sub-$\theta_n$6 as, and fluxes fall off exponentially with winding number. However, their signature in high-fidelity spectral line profiles (Fe K$\theta_n$7), interferometric visibility structures, and time-domain echoes (lag structures of variable sources) will become accessible with next-generation X-ray, mm-VLBI, and extended baseline observations (Falanga et al., 2021, Aratore et al., 2021, Wang et al., 2022).

The character of time-delay sequences (constant for photon-sphere spacetimes; decreasing for horizonless states), the scaling of magnification and positions, and the secular drift of image positions in Kerr or binary backgrounds provide direct and discriminating diagnostics for fundamental tests of the Kerr hypothesis, cosmic censorship, and the possible presence of alternative compact objects (Sahu et al., 2013, Patil et al., 2016). Precision modeling necessarily requires inclusion of these relativistic images, both as direct signatures and as sources of systematic bias in parameter inference (Falanga et al., 2021).