# Black hole photon ring beyond General Relativity: an integrable parametrization

Jibril Ben Achour, Eric Gourgoulhon, Hugo Roussille

► To cite this version:

Jibril Ben Achour, Eric Gourgoulhon, Hugo Roussille. Black hole photon ring beyond General Relativity: an integrable parametrization. Journal of Cosmology and Astroparticle Physics, 2025, 10, pp.012. 10.1088/1475-7516/2025/10/012. hal-05163702

HAL Id: hal-05163702

https://hal.science/hal-05163702v1

Submitted on 2 Nov 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
PAPER • OPEN ACCESS

# Black hole photon ring beyond General Relativity: an integrable parametrization

To cite this article: Jibril Ben Achour et al JCAP10(2025)012

View the [article online](#) for updates and enhancements.

You may also like

- *Machine and Deep Learning-driven Angular Momentum Inference from BHEX Observations of the n = 1 Photon Ring*
  Joseph Farah, Jordy Davelaar, Daniel Palumbo et al.

- *Multifrequency Models of Black Hole Photon Rings from Low-luminosity Accretion Disks*
  Tejahni Desire, Alejandro Cárdenas-Avendaño and Andrew Chael

- *Measuring Spin from Relative Photon-ring Sizes*
  Avery E. Broderick, Paul Tiede, Dominic W. Pesce et al.
# Black hole photon ring beyond General Relativity: an integrable parametrization

Jibril Ben Achour<sup>ID,a,b,c</sup> Éric Gourgoulhon<sup>ID,d,e</sup> and Hugo Roussille<sup>ID,c</sup>

<sup>a</sup>Arnold Sommerfeld Center for Theoretical Physics, Munich, Germany

<sup>b</sup>Munich Center for Quantum Sciences and Technology, Munich, Germany

<sup>c</sup>Université de Lyon, ENS de Lyon, Laboratoire de Physique, CNRS UMR 5672, Lyon 69007, France

<sup>d</sup>LUX, CNRS UMR 8262, Observatoire de Paris — PSL, Sorbonne Université Paris, 5 place Jules Janssen, 91190 Meudon, France

<sup>e</sup>Laboratoire de Mathématiques de Bretagne Atlantique, CNRS UMR 6205, Université de Bretagne Occidental, 6 avenue Victor Le Gorgeu, 29200 Brest, France

E-mail: jibril.ben-achour@ens-lyon.fr, eric.gourgoulhon@obspm.fr, hugo.roussille@ens-lyon.fr

**ABSTRACT:** In recent years, the shape of the photon ring in black holes images has been argued to provide a sharp test of the Kerr hypothesis for future black hole imaging missions. In this work, we confront this proposal to beyond Kerr geometries and investigate the degeneracy in the estimations of the black hole parameters using the circlipse shape proposed by Gralla and Lupsasca. To that end, we consider a model-independent parametrization of the deviations to the Kerr black hole geometry, dubbed Kerr off shell (KOS), which preserves the fundamental symmetry structure of Kerr known as the Killing tower. Besides exhibiting a Killing tensor and thus a Carter-like constant, all the representants of this family also possess a Killing-Yano tensor and are of Petrov type D. The allowed deviations to Kerr, selected by the symmetry, are encoded in two free functions which depend respectively on the radial and polar angle coordinates. Using the symmetries, we provide an analytic study of the radial and polar motion of photon trajectories generating the critical curve, to which the subrings composing the photon ring converge. This allows us to derive a ready-to-use closed formula for the parametric critical curve in term of the free functions parametrizing the deviations to Kerr. Using this result, we confront the circlipse fitting function to four examples of Kerr-like objects and we show that it admits a high degree of degeneracy. At a given inclination, the same circlipse can fit both a Kerr black hole of a given mass and spin (M, a) or a modified rotating black hole with different mass and spin parameters (M, a) and a new parameter α. Therefore, future tests of the Kerr hypothesis could be achieved only provided one can measure independently the mass and spin of the black hole to break this degeneracy.

**KEYWORDS:** Exact solutions, black holes and black hole thermodynamics in GR and beyond, gravity, GR black holes, astrophysical black holes

ARXIV ePRINT: 2506.09882

<figure><img src="image_2.png" /></figure>
Contents

<table><tbody><tr><td>1</td><td>The Kerr off shell family</td><td>4</td></tr><tr><td>1.1</td><td>Preserving Kerr hidden symmetries</td><td>4</td></tr><tr><td>1.2</td><td>Comparison with previous parametrizations</td><td>7</td></tr><tr><td>2</td><td>Geodesic motion and photon ring</td><td>8</td></tr><tr><td>2.1</td><td>Generalized Carter constant and geodesic equations</td><td>8</td></tr><tr><td>2.2</td><td>Photon ring of the Kerr off-shell metrics</td><td>11</td></tr><tr><td>2.3</td><td>Screen angular coordinates and the critical curve</td><td>15</td></tr><tr><td>3</td><td>Parameter estimations beyond Kerr: concrete examples</td><td>18</td></tr><tr><td>3.1</td><td>The Kerr black hole</td><td>18</td></tr><tr><td>3.2</td><td>Radial deformations to Kerr</td><td>20</td></tr><tr><td>3.3</td><td>Polar deformations to Kerr</td><td>24</td></tr><tr><td>4</td><td>Discussion</td><td>25</td></tr></tbody></table>

The first interferometric observations of the supermassive black hole at the center of the M87 galaxy by the Event Horizon Telescope (EHT) in 2019 [1–3], followed by the first image of the Sagittarius A* black hole at the center of our galaxy in 2022 [4, 5], has opened a new era for horizon scale electromagnetic astronomy and a new window to test the strong field regime of gravity. These observations have revealed the existence of an unresolved and asymmetric bright ring, as expected from the theoretical predictions and simulations. Since then, important efforts have been devoted to clarify how these interferometric signatures of ultra-compact objects could be used to test General Relativity (GR) and its different extensions<fnref n="1" /> [32–44].

The gravitational lensing induced by a Kerr black hole<fnref n="2" /> on photons emitted from a generic source possesses universal properties predicted by GR. In particular, photons can be momentarily captured and loop arbitrary many times around the black hole before escaping to infinity. These photons generically produce on the screen detector an infinite sequence of photon subrings corresponding to the arbitrary number of loops made around the hole [47]. These subrings converge towards a curve called the *critical curve* which stands as a key prediction of GR first derived in [48]. In this work, the *photon ring* will refer to the limit of the photon subrings and therefore be confounded with the critical curve.<fnref n="3" /> It depends on the black hole mass *M*, spin *a* and the observer inclination $y_O$ (w.r.t the black hole rotation axis). This critical curve, being the limit of the sequence of photon subrings, is

<footnote n="1">¹Important efforts have focused on investigating the shape of the photon ring of exotic compact objects, among which hairy black holes [6–20], wormholes [21–24], boson stars [25, 26], parity-violating black holes [27], regular black holes [28] and quantum corrected black holes [29]. To a large extent, these study use numerical ray-tracing methods to determine the shadow of the ultra-compact object. Theory-agnostic results relating the size of black hole hair to the size of the photon ring have been presented in [30, 31].</footnote>

<footnote n="2">²See [45, 46] for the study of the lensing induced by a Kerr-Newman black hole.</footnote>

<footnote n="3">³In the recent literature, the term “photon ring” is rather used for the set of all photon subrings [49, 50].</footnote>
not an observable. Nevertheless, as the sequence converges quickly, the critical curve is the fundamental theoretical object upon which an observable can be built (see ref. [50] for a discussion and limitations). In particular, it is expected to be contained in the so far unresolved bright ring observed by the EHT. This challenges the possibility to directly test the shape of the photon ring, and thus the GR prediction, with current data. Observation of the primary photon ring and subsequent subrings will require a space-based interferometric mission [50, 51].

In the last five years, important efforts have been done to i) characterize universal observables free from astrophysical bias allowing for a sharp test of GR [49, 52–58], and ii) to clarify the feasibility of a future space-based mission aiming at measuring these observables [50, 59–62]. In a series of works, it was shown how the intensity profile of the photon (sub)ring image (the amplitude and phase of the signal) depends on its geometrical shape (the diameter and width of the ring) [54, 55]. This relation allows one to reconstruct the shape of the photon ring starting from possibly partial information, for instance if one has only access to the amplitude of the signal. Prospects for the measurement of this universal signatures using a space-based interferometric mission was further explored with encouraging results, motivating the recent Black Hole Explorer (BHEX) proposal [51, 63]. With these exciting recent developments allowing to test GR and its extensions using photon ring science in the near future, it becomes crucial to provide tractable ready-to-use parametrization of beyond Kerr geometries. The key goal is to provide a sufficiently general parametrization encompassing a large set of possible and physically reasonable deviations, while at the same time allowing one to analytically derive the photon ring observables mentioned above, i.e. the critical curve and its intensity profile.

To that end, it is crucial that such new parametrization preserves the explicit and hidden symmetries of the Kerr geometry which are responsible for the integrability of the geodesic motion, and thus to obtain close analytical expression for the critical curve. A first initial agnostic parametrization allowing for a Carter-like constant was presented by Johanssen and Psaltis (JP) in [64, 65], followed by a series of works generalizing the JP proposal [66–72]. So far all these new proposals have focused on including deformation to the radial potential experienced by test particles while still preserving the existence of a non-trivial rank-2 Killing tensor which generates a Carter-like constant. However, this is not the most fundamental hidden symmetry of the Kerr black hole. Indeed, the Carter charge and the associated Killing tensor of the Kerr geometry, respectively identified by Carter in [73] and Walker and Penrose in [74], actually descend from a hierarchy of Killing objects known as the Killing tower. This structure was first identified in the context of higher dimensional rotating black holes in [75–77]. Concretely, the Killing tower is built from the existence of a non-degenerate conformal Killing-Yano rank-2 tensor from which all other Killing objects can be canonically derived. Interestingly, the family of geometries which admits such Killing tower (and not only the existence of its Killing tensor) has been fully characterized and is dubbed the *Kerr off shell* (KOS) family (see [78] for a review). The associated metric is given by (1.1) below. The goal of this work is to use this new KOS family of geometries as a new agnostic parametrization for Kerr-like objects and show the new interesting features it allows. There are two main outcomes coming from this new parametrization.
First, the KOS family is parametrized by two free functions denoted hereafter ($\Delta_r(r), \Delta_y(y)$) where $y \propto \cos\theta$ refers to the polar angle. These free functions modify respectively the radial and polar potentials (2.10b) experienced by test particles. Comparing to the previous parametrization [67–71], the free function $\Delta_r(r)$ encompassed a large set of radial deformations treated so far and coincides with the deformation introduced in [71] (denoted $A_5(r)$ there). However, to our knowledge, the polar deformation encoded in $\Delta_y(y)$ is new. Therefore, the KOS parametrization allows for a new deformation of the Kerr geometry while still preserving its most fundamental hidden symmetries. As we shall see, the new polar deformation allows to parametrize new subtle modifications affecting the critical angles $\theta_{\pm}$ characterizing the polar motion of the test particles which plays a non-trivial role when computing the image of source lensed by a Kerr-like object. Notice that general polar deformations to Kerr have already been considered in [79], but the parametrization introduced in that work does not preserve the rank-2 Killing tensor.

Second, since the KOS family is defined by the Killing tower structure, the geodesic motion is integrable and one can derive a generalized Carter constant given by (2.6). Quite remarkably, the form of the KOS metric and the simple expression for its Killing objects allows one to obtain the closed formula (2.47) describing the parametric shape of the critical curve in terms of the free functions ($\Delta_r(r), \Delta_y(y)$). This close formula is the main result of this work. It provides a ready-to-use parametric expression for the critical curve which allows us to test any given deformations encompassed in the KOS family.

In the second part of this work, we use this new general parametric expression for the KOS critical curve to explore the robustness of testing GR using the shape of the photon ring as proposed by Gralla and Lupsasca in [55]. The strategy adopted in [55, 59] relies on fitting the parametric expression for the critical curve by a *phoval curve* (or equivalently a circlinepse) on the observer’s screen. The even contribution to this fitting curve encoding the projected diameter $d_\phi$, is characterized by three real parameters ($R_0, R_1, R_2$) while its odd contribution, i.e. the centroid $C_\phi$, involves two additional parameters. The idea is that focusing on the even contribution $d_\phi$, one can fit to a great accuracy the Kerr critical curve and relate the three parameters ($R_0, R_1, R_2$) to the black hole mass and spin ($M, a$) and to the observer inclination $y_O$. In turn, the amplitude of the visibility (the Fourier transform) of the photon ring image is uniquely related to the profile $d_\phi$ [54]. While this suggests that the measurement of the shape of the photon ring can provide a sharp test of the Kerr hypothesis, and thus of GR, the procedure which consists in introducing a fitting *phoval* function could a priori remains blind to other Kerr-like geometries and thus admits some degeneracy in the parameter estimations. To understand these possible degeneracies, we confront the *phoval fitting function* proposed in [59] to four examples of Kerr-like geometries which involve a new parameter (additionally to the mass and spin) and which belong to our KOS family. Besides two well known examples given by the Kerr-MOG geometry and the regular Simpson-Visser rotating black hole, we consider two new modifications of Kerr which have not been treated so far and which are allowed by the KOS parametrization. For each of the four examples, we present the following results:

* First, we provide analytical expressions for the radial and polar extend of the photon ring and for the angular momentum $l_c$ and Carter constant $k_c$ labeling critical spherical
photon orbits. Notice that these quantities are usually not accessible due to either the
high complexity of the model or to its non-integrability.

* Second, for each examples, we test the degeneracy in the parameter estimation using the circlipse fitting function. We show that at maximal inclination, i.e. $y_O = 0$ which corresponds to the equatorial plane, an accurate fit of the Kerr black hole for a given $(M, a)$ can always be degenerate with another equally accurate fit for the Kerr-like object for another set of parameters $(M, a, \alpha)$ where $\alpha$ is the new parameter encoding the deviation to Kerr. Finally, we provide a plot of the shadow for each examples for different values of the new parameter $\alpha$.

This article is organized as follows. In section 1, we present the *Kerr off shell* (KOS) family involving the two free functions and we review the Killing tower structure. We further compare this new parametrization of Kerr-like objects with previous works. Section 2 is devoted to the geodesic motion in the generic KOS geometry. We derive the Killing tensor charge and discuss the radial and polar motion and their associated potentials (2.11a)–(2.11b). Finally, we present a detail discussion on the photon ring structure in terms of the two free functions. In section 2.3, we derive the parametric form of the critical curve, given in (2.47), for the KOS family. The last section 3 is devoted to the study of several concrete examples of beyond Kerr geometries which are used to test the limitations of the *phoval fitting function*. We conclude by a discussion of our results and the perspectives it opens.

# 1 The Kerr off shell family

In this section, we introduce the new parametrization for Kerr-like objects which we will use
throughout this work: the Kerr off-shell family of geometries. We present briefly its ruling
structure, namely the Killing tower and we compare this new beyond-Kerr-parametrization
to the existing ones.

## 1.1 Preserving Kerr hidden symmetries

By virtue of the no-hair theorem of four dimensional GR (see [80] for a review), the Kerr
geometry is the unique asymptotically flat stationary and axi-symmetric vacuum solution
of GR. Therefore, it stands as our best model for the observed rotating ultracompact
astrophysical objects, in particular supermassive black holes sitting at the center of galaxies.
Nevertheless, even if no observational evidence point towards any deviation to this model
so far, it is crucial to provide efficient parametrization of possible deviations that could
be captured in the future.

As is well known, the Kerr geometry enjoys a set of explicit and hidden symmetries which
play a key role in the derivations of its properties. Besides the time and azimuthal isometries,
it also admits a Killing-Yano two-form which canonically equips the geometry with a Killing
tensor. This later object generates the well-known Carter constant which makes the geodesic
motion integrable [73, 81, 82], allowing us to classify the different allowed trajectories of
massive and massless particles on the Kerr geometry [83, 84]. In turn, the presence of a Killing-
Yano tensor generates a set of conserved charges which allows one to integrate the motion of
any test spinning particle<fnref n="4" /> (up to linear order in the spin) on the Kerr background [87, 88], allowing us to study properties of extreme-mass-ratio inspirals [89, 90]. Finally, the same hidden Killing symmetries are also responsible for the well-known separability of spin-s wave equations allowing one to study the scattering problem of various test fields with the Kerr black hole [91]. Therefore, testing the Kerr hypothesis can be understood as testing the fundamental symmetries of the Kerr black hole. From this point of view, it appears interesting to search for a parametrization of Kerr-like geometries which preserves the Kerr symmetries. Consider the following stationary and axi-symmetric metric

$$ds^2 = -\frac{\Delta_r}{\Sigma} (d\tau + y^2 d\varphi)^2 + \frac{\Delta_y}{\Sigma} (d\tau - r^2 d\varphi)^2 + \frac{\Sigma}{\Delta_r} dr^2 + \frac{\Sigma}{\Delta_y} dy^2 \quad (1.1)$$

where

$$\Sigma = r^2 + y^2 \qquad \Delta_r := \Delta_r(r) \qquad \Delta_y := \Delta_y(y) \quad (1.2)$$

while the inverse metric reads

$$g^{\mu\nu} \partial_{\mu} \partial_{\nu} = \frac{1}{\Sigma} \left[ -\frac{1}{\Delta_r} (r^2 \partial_{\tau} + \partial_{\varphi})^2 + \frac{1}{\Delta_y} (y^2 \partial_{\tau} - \partial_{\varphi})^2 + \Delta_r (\partial_r)^2 + \Delta_y (\partial_y)^2 \right] \quad (1.3)$$

The metric is dubbed off shell because it does not have to be a solution of Einstein's equations: it is parametrized by the two free functions $\Delta_r(r)$ and $\Delta_y(y)$. The metric (1.1) is a subset of a general class of metrics considered by Carter [73] to ensure the separability of the Hamilton-Jacobi and Schrödinger equations. In this subset, the function $\Sigma$ is not modified w.r.t the Kerr case. Contrary to the geometries of the general class considered by Carter which are not algebraically special, this restriction ensures that the spacetime geometry is of Petrov type D.

While the functions $\Delta_r$ and $\Delta_y$ are free, they must satisfy several constraints in order for the spacetime (1.1) to describe our object of interest, i.e. an asymptotically flat black hole. First, asymptotic flatness constrains the asymptotic behavior of $\Delta_r$: one must have $\Delta_r \sim r^2$ as $r \to \infty$. Second, in order for the metric (1.1) to have an event horizon, $\Delta_r$ must be such that $\Delta_r(r) > 0$ for $r > r_h$ and $\Delta_r(r_h) = 0$; $r_h$ corresponds to the radius of the outermost horizon. Third, $\Delta_y$ must be positive for all possible values of $y$ to conserve the Lorentzian signature. Finally, since the metric exhibits a curvature singularity for $\Sigma = 0$, we assume $r_h > 0$ in order to avoid the presence of a naked singularity. Furthermore, a metric satisfying these constraints might still exhibit a conical or Misner string singularity, depending on the exact shape of the functions $\Delta_r$ and $\Delta_y$. Requiring these singularities to be removable yields additional constraints on these functions [92, 93] which would need to be further explored.

The Kerr metric belongs to this family. To see this, one can perform the following change of coordinates<fnref n="5" />

$$y = a \cos \theta, \qquad \varphi = \frac{\phi}{a}, \qquad \tau = t - a\phi, \quad (1.4)$$

<footnote n="4">⁴A set of conserved charges for an extended quadrupolar test object on the Kerr background was found in [85, 86] assuming that the test object is a black hole. It was further shown in [87] that a sufficient number of these conserved charges are in involution, ensuring the integrability of this system.</footnote>

<footnote n="5">⁵We can see that the change of variables (1.4) is singular in the Schwarzschild limit $a \to 0$. Nevertheless, one does indeed recover the Schwarzschild black hole by performing the change of variables (1.4), the choice (1.5) and then taking the $a \to 0$ limit. One can note that such parameter singularities are common when one recovers common solutions from the Plebanski-Demianski family of spacetimes [94].</footnote>
which relates the coordinates $(\tau, r, y, \varphi)$ for the Boyer-Lindquist coordinates $(t, r, \theta, \phi)$. Then, the Kerr metric is recovered by setting

$$ \Delta_r^{\text{Kerr}} = r^2 + a^2 - 2Mr, \quad \Delta_y^{\text{Kerr}} = a^2 - y^2. \qquad (1.5) $$

Within this family, one can also find all the other representatives of the Plebianski-Demianski family which are labeled by the NUT charge, the acceleration and the cosmological constant [73]. We shall consider these extended family in the last section. The key interest in working with this family of geometries is that it is the most general set of geometries which possess the same hidden symmetries as the Kerr metric. These hidden symmetries are encoded in the so called Killing tower structure.

To understand the structure, it is crucial to understand that all the symmetries of Kerr descend (can be canonically derived) from the existence of a principal tensor, i.e. a closed non-degenerate conformal Killing-Yano 2-form which is defined by the following equations:

$$ \nabla_{\mu} p_{\nu\alpha} = 2g_{\mu[\nu}h_{\alpha]}, \quad \nabla_{[\mu}p_{\nu\alpha]} = 0, \quad h_{\alpha} = \frac{1}{3}\nabla_{\mu}p^{\mu}_{\alpha}. \qquad (1.6) $$

For all the members of the off-shell Kerr family, this tensor takes the form<fnref n="6" />

$$ p_{\mu\nu}dx^{\mu} \wedge dx^{\nu} = ydy \wedge (d\tau - r^2 d\varphi) - rdr \wedge (d\tau + y^2 d\varphi). \qquad (1.7) $$

Using this object, one can construct the associated (dual) KY tensor, i.e. $f_{\mu\nu} = \frac{1}{2}\epsilon_{\mu\nu\alpha\beta}p^{\alpha\beta}$, which satisfies $\nabla_{(\mu}f_{\nu)}_{\alpha} = 0$ and which reads

$$ f_{\mu\nu}dx^{\mu} \wedge dx^{\nu} = rdy \wedge (d\tau - r^2 d\varphi) + ydr \wedge (d\tau + y^2 d\varphi). \qquad (1.8) $$

These two anti-symmetric tensors can then be combined to construct a conformal Killing tensor and a Killing tensor. The Killing tensor which will generate the Carter-like constant is given by

$$ K_{\mu\nu}dx^{\mu}dx^{\nu} = f_{\mu\alpha}f^{\alpha}_{\nu}dx^{\mu}dx^{\nu}, \\ = \frac{1}{\Sigma} \left[ y^2 \Delta_r (d\tau + y^2 d\varphi)^2 + r^2 \Delta_y (d\tau - r^2 d\varphi)^2 \right] + \Sigma \left[ \frac{r^2 dy^2}{\Delta_y} - \frac{y^2 dr^2}{\Delta_r} \right], \qquad (1.9) $$

and satisfies $\nabla_{(\mu}K_{\nu\alpha)} = 0$. Finally, one can identify the killing vectors (KV) generating the time and azimuthal isometries from the above objects. Indeed, the so-called primary vector and the second vector can be derived from the principal tensor $p_{\mu\nu}$ and the Killing tensor $K_{\mu\nu}$ via the relations

$$ \xi^{\alpha} \partial_{\alpha} = \frac{1}{3} \nabla_{\mu} p^{\mu \alpha} \partial_{\alpha} = \partial_{\tau}, \qquad (1.10) $$

$$ \chi^{\alpha} \partial_{\alpha} = -K^{\alpha}_{\mu} \xi^{\mu} \partial_{\alpha} = \partial_{\varphi}. \qquad (1.11) $$

One recovers in this way the two commuting KVs of this stationary and axisymmetric spacetime. This elegant structure was initially identified and explored in [75–77] when constructing higher dimensional rotating black hole geometries preserving the Kerr symmetries

<footnote n="6">⁶Computations performed in this section have been checked by means of a SageMath notebook posted at https://nbviewer.org/url/relativite.obspm.fr/notebooks/Kerr_off_shell.ipynb.</footnote>
(see the review [78], page 44, section 3.2, for more details). This KOS family of four-dimensional geometries is uniquely defined by the existence of this Killing tower [75]. Therefore, it naturally provides a new parametrization for Kerr-like objects which is selected solely on symmetry arguments. Let us now discuss the main novelty it provides compared to the parametrizations proposed so far.

## 1.2 Comparison with previous parametrizations

The search for a suitable parametrization of Kerr-like objects was initiated by Johannssen and Psaltis in [95]. Modifications to the initial JP proposal were studied in [66–72], providing new parametrizations preserving some of the Kerr symmetries, in particular allowing for a Carter-like constant. In the following, we compare the KOS parametrization with the latest parametrizations proposed so far.

A general parametrization of Kerr-like objects was explored in [71]. Consider the metric

$$ds^2 = g_{tt}dt^2 + 2g_{t\phi}d\phi + g_{rr}dr^2 + g_{\theta\theta}d\theta^2 + g_{\phi\phi}d\phi^2, \qquad (1.12)$$

where the metric functions read explicitly

$$g_{tt} = -\frac{\tilde{\Sigma}(A_5 - a^2 A_2^2 \sin^2 \theta)}{\rho^4}, \qquad g_{t\phi} = \frac{a A_5 (A_5 - A_0) \tilde{\Sigma} \sin^2 \theta}{\rho^4}, \qquad g_{\theta\theta} = \tilde{\Sigma}, \qquad (1.13a)$$

$$g_{\phi\phi} = \frac{\tilde{\Sigma} \sin^2 \theta A_5 (A_1^2 - a^2 A_5 \sin^2 \theta)}{\rho^4}, \qquad g_{rr} = \frac{\tilde{\Sigma}}{A_5}, \qquad (1.13b)$$

with

$$\tilde{\Sigma} = r^2 + a^2 \cos^2 \theta + f(r) + g(\theta), \qquad (1.14)$$

$$\rho^4 = a^4 A_2^2 A_5 \sin^4 \theta + a^2 \sin^2 \theta (A_0^2 - 2A_0 A_5 - A_1^2 A_2^2) + A_1^2 A_5. \qquad (1.15)$$

The functions $A_i(r)$, $f(r)$ and $g(\theta)$ are arbitrary.<fnref n="7" /> The Petrov type of this general geometry was studied in detail in [71], showing that, without further conditions, such geometry is not algebraically special, i.e. is of Petrov type I. Within this general family, the authors identified two restricted classes which admit a Carter-like constant.

* The first class (hereafter Class I) is given by imposing

$$A_0 = A_1 = r^2 + a^2, \qquad A_2 = 1, \qquad f = g = 0 \qquad (1.16)$$

leading to

$$g_{tt} = -\frac{(A_5(r) - a^2 \sin^2 \theta)}{\Sigma}, \qquad g_{t\phi} = \frac{a(r^2 + a^2 - A_5(r)) \sin^2 \theta}{\Sigma}, \qquad g_{\theta\theta} = \Sigma, \qquad (1.17a)$$

$$g_{\phi\phi} = \frac{(r^2 + a^2)^2 - a^2 A_5(r) \sin^2 \theta}{\Sigma}, \qquad g_{rr} = \frac{\Sigma}{A_5(r)}, \qquad (1.17b)$$

where $\Sigma = r^2 + a^2 \cos^2 \theta$. A direct comparison shows that this class is a sub-case of the KOS family such that

$$\Delta_r(r) = A_5(r), \qquad \Delta_y(y) = \Delta_y^{\text{Kerr}}(y) = a^2 - y^2 = a^2 \sin^2 \theta. \qquad (1.18)$$

<footnote n="7">⁷The initial JP parametrization is recovered for $A_0 = A_1 A_2$.</footnote>
This family of geometries only contains radial deformations while remaining blind to any polar corrections. Nevertheless, it contains a large class of corrections studied so far, in particular all the regular black holes built by deforming the function $\Delta_r^{\text{Kerr}}(r) = r^2 + a^2 - 2Mr$ to $\Delta_r(r) = r^2 + a^2 - 2M(r)r$ with a suitable function of $M(r)$. As shown in [71], all members of this class are of Petrov type D. The fact that all these geometries belong to the KOS family implies that each of them admits not only a Killing tensor, but a killing tensor which descends from a Killing-Yano tensor. This is coherent with the known theorems relating the existence of KY two-form and the Petrov type of a given geometry. Indeed, only Petrov type D, N, and O geometries are equipped with a KY two form [96].

* The second restricted Class II corresponds to

$$A_0 = A_1 = r^2 + a^2 + f(r), \quad A_2 = 1, \quad g = 0 \qquad (1.19)$$

while all the functions $A_5(r)$ and $f(r)$ are left free. As shown in [71], these geometries are of Petrov type I. As an example, one can obtain the Kerr-Sen black hole by setting all the free functions to zero while imposing $f(r) = \alpha r$ and $A_5(r) = r^2 + a^2 - 2Mr + \alpha r$. Since this class of geometries allows for corrections to the function $\Sigma = r^2 + a^2 \cos^2 \theta$, and is of Petrov type I, it cannot be covered by the KOS family which is by construction Petrov type D. Therefore, representatives of this second class provide interesting examples of non-algebraically special geometry which admits a Killing tensor but no Killing-Yano twoform.

Therefore, compared to the parametrization studied in [71], the KOS family allows us to single out the Petrov type D geometries belonging to the restricted class I, but also a whole new set of geometries admitting polar deformations while still equipped with a Killing-Yano tensor.

To conclude, we point that a parametrization of Kerr-like objects with general polar and radial deformations have been studied in [79]. However, the allowed freedom is so general that the considered geometries fail to be algebraically special and thus preserve the Kerr symmetries tied to the Petrov type D. The polar deformation $\Delta_y(y)$ of the KOS family can be understood as the only allowed polar deformation selected by the Killing tower structure. We are now ready to study the geodesic motion and the photon ring structure for this KOS family.

## 2 Geodesic motion and photon ring

In this section, we study the geodesic motion on a general KOS geometry in the Hamiltonian formalism. We derive the Carter-like constant and discuss how the deformations $(\Delta_r(r), \Delta_y(y))$ affect the radial and polar motion of test particles. Then, we study the unstable radial and polar positions for the geodesic motion and provide close formula for the allowed range of the radial and polar coordinates for such unstable equilibriums.

### 2.1 Generalized Carter constant and geodesic equations

Consider the general KOS metric (1.1). The phase space of the geodesic motion on this background is described by the coordinates $\tau, r, y, \varphi$ and the impulsions $p_\tau, p_r, p_y$ and $p_\varphi$.
The canonical commutation relations read

$$ \{\tau, p_{\tau}\} = 1, \quad \{r, p_{r}\} = 1, \quad \{y, p_{y}\} = 1, \quad \{\varphi, p_{\varphi}\} = 1, \quad (2.1) $$

while all other Poisson brackets are zero. Notice that

$$ p_r = \frac{\Sigma}{\Delta_y} p^y \quad \text{and} \quad p_r = \frac{\Sigma}{\Delta_r} p^r. \qquad (2.2) $$

The geodesic Hamiltonian generating the dynamics of test particles is given by

$$ H = \frac{1}{2} g^{\mu\nu} p_{\mu} p_{\nu} = \frac{1}{2\Sigma} \left[ -\frac{1}{\Delta_r} (r^2 p_{\tau} + p_{\varphi})^2 + \frac{1}{\Delta_y} (y^2 p_{\tau} - p_{\varphi})^2 + \Delta_r p_r^2 + \Delta_y p_y^2 \right]. \quad (2.3) $$

As expected for a stationary and axi-symmetric metric, the variables $(\tau, \phi)$ are cyclic, such that $p_{\tau}$ and $p_{\varphi}$ are constants of motion, i.e. $\{p_{\tau}, H\} = 0$ and $\{p_{\varphi}, H\} = 0$. They correspond to the conserved charges induced by the time and azimuthal Killing vectors $\xi^{\mu}\partial_{\mu} = \partial_{\tau}$ and $\chi^{\mu}\partial_{\mu} = \partial_{\varphi}$ given by

$$ E = -\xi^{\mu} p_{\mu} = -p_{\tau} \quad L = \chi^{\mu} p_{\mu} = p_{\varphi}. \qquad (2.4) $$

Notice that since our coordinates $(\tau, \varphi)$ are not the standard BL coordinates $(t, \phi)$, the above charges $(E, L)$ do not coincide with the standard energy and angular momentum derived in the BL coordinates.<fnref n="8" /> Then, the existence of Killing tensor (descending from the Killing-Yano twoform) provides another conserved charge quadratic in the momenta. This generalized Carter constant is given by

$$ K = K^{\mu\nu} p_{\mu} p_{\nu} = \frac{1}{\Sigma} \left[ \frac{y^2}{\Delta_r} (r^2 p_{\tau} + p_{\varphi})^2 + \frac{r^2}{\Delta_y} (y^2 p_{\tau} - p_{\varphi})^2 - y^2 \Delta_r p_r^2 + r^2 \Delta_y p_y^2 \right]. \quad (2.6) $$

One can verify that $\{K, H\} = 0$, which means that the quantity $K$ is indeed conserved along a geodesic. It is useful for the following to simplify this constant of motion. To that end, we use the trivial constant of motion given by the mass $m$ of the particle, i.e. its hamiltonian, such that $2H = p_{\mu}p^{\mu} = -m^2$. Using (2.2) and (2.3), one can simplify as:

$$ K = \frac{1}{\Delta_y} \left[ (y^2 p_{\tau} - p_{\varphi})^2 + \Sigma^2 (p^y)^2 \right] + y^2 m^2 = \frac{1}{\Delta_y} \left[ (y^2 E + L)^2 + \Sigma^2 (p^y)^2 \right] + y^2 m^2, \quad (2.7) $$

or equivalently as

$$ K = \frac{1}{\Delta_r} \left[ (r^2 p_{\tau} + p_{\varphi})^2 - \Sigma^2 (p^r)^2 \right] + r^2 m^2 = \frac{1}{\Delta_r} \left[ (r^2 E - L)^2 - \Sigma^2 (p^r)^2 \right] + r^2 m^2. \quad (2.8) $$

These four constants of motion $(p_{\tau}, p_{\varphi}, K, H)$ commute with each other, which imply that only 4 out of the initial 8 degrees of freedom are left in the phase space, thus making the

$$ E = E^{\text{BL}} \quad L = a (L^{\text{BL}} - aE^{\text{BL}}) \qquad (2.5) $$

which will be useful when comparing to the Kerr case.

<footnote n="8">⁸Using the change of coordinates (1.4), one finds that the two Killing vectors are given by $\xi^{\mu}\partial_{\mu} = \partial_{t}$ and $\chi^{\mu}\partial_{\mu} = a(\partial_{\phi} + a\partial_{t})$ in terms of BL coordinates. Introducing the momenta $(p_{t}, p_{\phi})$ and the energy $E^{\text{BL}} = -p_{t}$ and angular momentum $L^{\text{BL}} = p_{\phi}$ associated to the BL coordinates, one obtains therefore</footnote>
geodesic motion integrable. With the phase space at hand, we can now derive the dynamics of a test particle.

Let us introduce an affine parameter denoted $\lambda$. Hamilton's equations for $\tau$, $r$, $y$ and $\varphi$ are given by

$$ \Sigma \frac{d\tau}{d\lambda} = -\frac{r^2(r^2 p_{\tau} + p_{\varphi})}{\Delta_r} + \frac{y^2(y^2 p_{\tau} - p_{\varphi})}{\Delta_y}, \qquad (2.9a) $$

$$ \Sigma \frac{dr}{d\lambda} = \Delta_r p_r, \qquad (2.9b) $$

$$ \Sigma \frac{dy}{d\lambda} = \Delta_y p_y, \qquad (2.9c) $$

$$ \Sigma \frac{d\varphi}{d\lambda} = -\frac{r^2 p_{\tau} + p_{\varphi}}{\Delta_r} - \frac{y^2 p_{\tau} - p_{\varphi}}{\Delta_y}, \qquad (2.9d) $$

These equations can be further simplified by expressing them in terms of the conserved charges $(E, L, K, m)$. One can extract the expression for $p_r$ and $p_y$ from the expressions for $m^2$ and $K$, yielding

$$ \Sigma \frac{d\tau}{d\lambda} = \frac{r^2(r^2 E - L)}{\Delta_r} - \frac{y^2(y^2 E + L)}{\Delta_y}, \qquad (2.10a) $$

$$ \Sigma \frac{dr}{d\lambda} = \pm \sqrt{V_r(r)}, \qquad (2.10b) $$

$$ \Sigma \frac{dy}{d\lambda} = \pm \sqrt{V_y(y)}, \qquad (2.10c) $$

$$ \Sigma \frac{d\varphi}{d\lambda} = \frac{r^2 E - L}{\Delta_r} + \frac{y^2 E + L}{\Delta_y}, \qquad (2.10d) $$

where the two potentials $V_r(r)$ and $V_y(y)$ are given by

$$ V_r(r) = (r^2 E - L)^2 - \Delta_r (K + r^2 m^2), \qquad (2.11a) $$

$$ V_y(y) = \Delta_y (K - m^2 y^2) - (y^2 E + L)^2. \qquad (2.11b) $$

This provides the general set of first order equations describing the geodesic motion on an arbitrary representant of the KOS family. The two potentials (2.11a) and (2.11b) govern respectively the radial and polar motion of the test particle.<fnref n="9" /> The classification of the geodesic

$$ \left(\frac{d\theta}{d\lambda}\right)^2 = \frac{1}{a^2 - y^2} \left(\frac{dy}{d\lambda}\right)^2. $$

This leads to the following equation:

$$ \Sigma^2 \left(\frac{d\theta}{d\lambda}\right)^2 = \Theta(y), \quad \Theta(y) = \frac{\Delta_y (K - m^2 y^2) - (y^2 E + L)^2}{a^2 - y^2}. $$

One recovers the Kerr potentials using (2.5) and (1.18), i.e.

$$ V_r(r) = \mathcal{R}(r) = [(r^2 + a^2)E^{\text{BL}} - aL^{\text{BL}}]^2 - \Delta(r^2 m^2 + K), $$

$$ V_y(y) = \Theta(\theta) = K - \left( \frac{L^{\text{BL}}}{\sin\theta} - aE^{\text{BL}} \sin\theta \right)^2 - m^2 a^2 \cos^2\theta. $$

<footnote n="9">⁹In order to compare (2.11a) and (2.11b) to the standard potentials involved in the radial and polar motion around Kerr potentials, it is useful to write the equation for $\theta(\lambda)$ instead of $y(\lambda)$, using that $\theta = \arcsin(y/a)$ such that</footnote>
motion on the Kerr black hole has a long story [48, 97–103], a full characterization of the set of allowed trajectories being still an active research field [83, 84]. See [104] for a detailed discussion. The different Kerr geodesics can be classified in terms of the range of the constants of motion (E, L, K, m) and the spin a.<fnref n="10" /> A first condition for causal geodesics is that

$$V_r(r) \ge 0 \quad \text{and} \quad V_y(y) \ge 0 \qquad (2.14)$$

at any point (r, y) of the geodesic flow. Assuming that $\Delta_y(y) \ge 0$, one sees that the positivity of $V_y(y)$ implies that the generalized Carter constant is always positive or vanishing, i.e.

$$K \ge 0. \qquad (2.15)$$

In the Kerr case, where $\Delta_y^{\text{Kerr}}(y) = a^2 - y^2 = a^2 \sin^2 \theta \ge 0$, this condition is satisfied and one recovers the positivity of the Carter constant. For the Kerr black hole, one can split the study between the trajectories with $L = 0$ and the ones with $L \ne 0$. While the first are allowed to reach the rotation axis and cross it, the later family never reach it and are either confined between two critical polar angles or constrained to stay at a given angle $\theta$. One can then distinguish several families of trajectories with specific qualitative behavior depending on the range of the constants (E, L, K, a). See [83, 84, 105, 106] for details.

Obviously, such classification cannot be done in our case without selecting concrete functions $\Delta_r(r)$ and $\Delta_y(y)$. Nevertheless, one can derive the general expressions for the equilibrium positions of the radial and polar motion for a general KOS geometry and study the set of critical spherical photon trajectories which characterize the photon ring.

## 2.2 Photon ring of the Kerr off-shell metrics

The photon ring is generated by photon trajectories that asymptote (unstable) bound photon orbits — the so-called “photon spherical orbits” to be discussed below. Thus, let us fix the particle mass to zero, i.e. $m = 0$ and let us focus on the null geodesics with $E \ne 0$. This allows one to introduce the reduced constants of motion given by

$$\ell = \frac{L}{E} \quad \text{and} \quad k = \frac{K}{E^2}. \qquad (2.16)$$

$$Q = \frac{V_y(0)}{a^2} = \frac{K\Delta_y(0) - L^2}{a^2} = \frac{K\Delta_y(0) - a^2(L^{\text{BL}} - aE^{\text{BL}})^2}{a^2} \qquad (2.12)$$

which for Kerr reduces indeed to

$$Q = K - (L^{\text{BL}} - aE^{\text{BL}})^2. \qquad (2.13)$$

<footnote n="10">¹⁰Notice that it is standard to introduce another constant $Q$ which is well suited to study the equatorial plane, i.e. $y = 0$ or $\theta = \pi/2$. In our case, this alternative Carter constant can be written</footnote>
Imposing $m = 0$, and using the new constants, the reduced geodesic equations read<fnref n="11" />

$$ \frac{\Sigma}{E} \frac{d\tau}{d\lambda} = \frac{r^2(r^2 - \ell)}{\Delta_r} - \frac{y^2(y^2 + \ell)}{\Delta_y}, \qquad (2.17a) $$

$$ \frac{\Sigma}{E} \frac{dr}{d\lambda} = \pm \sqrt{\mathcal{V}_r(r)}, \qquad (2.17b) $$

$$ \frac{\Sigma}{E} \frac{dy}{d\lambda} = \pm \sqrt{\mathcal{V}_y(y)}, \qquad (2.17c) $$

$$ \frac{\Sigma}{E} \frac{d\varphi}{d\lambda} = \frac{r^2 - \ell}{\Delta_r} + \frac{y^2 + \ell}{\Delta_y}, \qquad (2.17d) $$

where the two potentials $\mathcal{V}_r(r)$ and $\mathcal{V}_y(y)$ are given by

$$ \mathcal{V}_r(r) = \frac{V_r(r)}{E} = (r^2 - \ell)^2 - \Delta_r k, \qquad (2.18a) $$

$$ \mathcal{V}_y(y) = \frac{V_y(y)}{E} = \Delta_y k - (y^2 + \ell)^2. \qquad (2.18b) $$

In the following, we identify the spherical photon orbits which play the role of the critical orbits beyond which photon can escape to infinity. The photon orbits which are infinitely close to the critical orbits are the ones generating the image of the photon ring. They are characterized by a set of parameters $(\ell_c, k_c)$ which we now determine.

### 2.2.1 Radial motion and spherical photon orbits

A photon spherical orbit<fnref n="12" /> is a null geodesic at fixed value, say $r_0$, of the coordinate $r$ [107].

Spherical photon orbits obey $dr/d\lambda = 0$ for all values of $\lambda$. According to (2.17b), this implies

$$ \mathcal{V}_r(r_0) = 0, \qquad \mathcal{V}'_r(r_0) = 0. \qquad (2.19) $$

These conditions translate into two equations for three unknowns $(r_0, \ell_c, k_c)$ given by

$$ (r_0^2 - \ell_c)^2 - \Delta_r(r_0)k_c = 0, \qquad (2.20a) $$

$$ 4r_0(r_0^2 - \ell_c) - \Delta'_r(r_0)k_c = 0. \qquad (2.20b) $$

Given some $r_0$, we look for the solutions $(\ell_c, k_c)$ labeling to allowed photon trajectories. Let us assume that $\Delta'_r(r_0) \neq 0$. Then, one can solve for $k$:

$$ k_c = \frac{4r_0(r_0^2 - \ell_c)}{\Delta'_r(r_0)}. \qquad (2.21) $$

Injecting this solution in the last equation, one obtains

$$ (\ell_c - r_0^2) \left( \ell_c - r_0^2 + \frac{4r_0 \Delta_r(r_0)}{\Delta'_r(r_0)} \right) = 0, \qquad (2.22) $$

<footnote n="11">¹¹Since we are focusing on photons that are reaching a distant observer, we may set $|E| = E$ in the following equations, given that one has necessarily $E > 0$ for such photons (only photons confined to the ergoregion can have $E \le 0$).</footnote>

<footnote n="12">¹²Here, the denomination “spherical” should be understood in the topological sense: indeed, the metric induced by $(\tau = \text{cst}, r = \text{cst})$ on (1.1) is not the metric of a 2D sphere in general.</footnote>
which has two solutions. The first one is given by $(\ell_c, k_c) = (r_0^2, 0)$. Now injecting $k=0$ in the polar potential, one finds $V_y(y) = -(y^2 + \ell)^2$ which has to be positive for all $y$. Thus, for $y=0$, one has $\ell_c = 0$ which implies $r_0 = 0$ and one finds that this solution is not admissible because it is clearly under the outermost horizon. The second solution $(\ell_c, k_c)$ corresponds to

$$ \ell_c = r_0 \left( r_0 - \frac{4\Delta_r(r_0)}{\Delta'_r(r_0)} \right), \quad k_c = \frac{16r_0^2\Delta_r(r_0)}{\Delta'_r(r_0)^2}. \qquad (2.23) $$

Notice that as expected, $(\ell_c, k_c)$ do not depend on the free function $\Delta_y$. Now, in order to fully characterize the locus of points allowing for spherical photon orbits, we have to further impose that $V_y(y) \ge 0$.

At the equator, i.e. at $y=0$, one has

$$ \mathcal{V}_y(0) \ge 0 \quad \rightarrow \quad \Delta_y(0)k_c \ge \ell_c^2. \qquad (2.24) $$

The above condition translates into the following inequality

$$ \frac{r_0^2}{\Delta'_r(r_0)^2} \left[ 16\Delta_r(r_0)\Delta_y(0) - (r_0\Delta'_r(r_0) - 4\Delta_r(r_0))^2 \right] \le 0 \qquad (2.25) $$

At this stage, one needs to pick up a given function $\Delta_r(r)$ to find the associated positions $r_0$ satisfying this inequality. In the case of Kerr, one finds that the condition reduces to a cubic polynomial equation from which one can extract the range of $r_0$. In the following, we shall assume that the function $\Delta_r(r)$ is such that the largest set of values of $r_0$ allowed by this inequality belong to an interval

$$ r_0 \in [r_{\text{ph}}^{-}, r_{\text{ph}}^{+}], \qquad (2.26) $$

where we have

$$ r_{\text{ph}}^{-} < r_{\text{ph}}^{+}, \quad \Delta_{r}(r_{\text{ph}}^{-}) > 0 \quad \text{and} \quad \Delta_{r}(r_{\text{ph}}^{+}) > 0. \qquad (2.27) $$

This concludes the study of the spherical photon orbits in the KOS geometry. Given one representant of the family, one can determine the couple $(\ell_c, k_c)$ labeling these critical orbits and the radial positions $r_{ph}^{\pm}$ for which spherical orbits exist. At this stage, we are left with studying the polar motion.

### 2.2.2 Polar motion and stable equilibrium positions

In order to capture the different allowed polar motions, it is instructive to recast the polar equation of motion into an equation for a 1d particle in a potential. By setting $d\lambda' = |E|\sqrt{\Delta_y}d\lambda/\Sigma$, using the fact that $\Delta_y \ge 0$, this equation can be recast into

$$ \left(\frac{dy}{d\lambda'}\right)^2 + \Theta(y) = k \quad \text{with} \quad \Theta(y) = \frac{(y^2 + \ell)^2}{\Delta_y}. \qquad (2.28) $$

This can be interpreted as a one-dimensional motion in a potential $\Theta(y)$ and total energy $k$. The allowed trajectories therefore have to satisfy

$$ 0 \le \Theta(y) \le k. \qquad (2.29) $$
The derivatives of the potential are given by

$$ \Theta'(y) = \frac{y^2 + \ell}{\Delta_y^2} (4y\Delta_y - (y^2 + \ell)\Delta'_y) \quad (2.30a) $$

$$ \Theta''(y) = \frac{8y^2}{\Delta_y} + \frac{4(y^2 + \ell)}{\Delta_y^2}(\Delta_y - 2y\Delta'_y) + \frac{(y^2 + \ell)^2}{\Delta_y^3}(2\Delta_y'^2 - \Delta_y\Delta_y'') \quad (2.30b) $$

Let us denote the different possible roots of $\Delta_y(y)$ by $y_*^\pm$. These positions bound the available domain for $y$. In the case of Kerr, they correspond to the axis of rotation.

In order to understand the polar motion we need to identify the equilibriums. A polar equilibrium at $y = y_e$ satisfies $\Theta(y_e) = 0$ and $\Theta'(y_e) = 0$. Depending on the stability of the various equilibriums, several cases arise. For simplicity, we assume that $\Delta_y$ is concave at $y = 0$: this implies $\Delta_y''(0) < 0$.

* First, we see that $y_e = 0$ is always a solution: indeed, since $\Delta_y$ is even, one has $\Delta'_y(0) = 0$. One then obtains

$$ \Theta_y''(0) = \frac{\ell}{\Delta_y(0)} \left( 4 - \ell \frac{\Delta_y''(0)}{\Delta_y(0)} \right) . \quad (2.31) $$

The position is stable only if $\Theta_y''(0) > 0$. This corresponds to

$$ \ell > 0 \quad \text{or} \quad \ell < \frac{4\Delta_y(0)}{\Delta_y''(0)}. \quad (2.32) $$

If $\ell = 0$, stability depends on the third derivative of $\Theta$.

* Second, we see that when $\ell < 0$, $y_e = \pm\sqrt{-\ell}$ are two solutions as well. When they exist, these solutions are always stable, since $\Theta''(\pm\sqrt{-\ell}) > 0$.

Depending on the exact expression of $\Delta_y$, there might be other positions for which $\Theta' = 0$. For simplicity, we assume that this is not the case. The possible movements then depend on the value of $\ell$:

* if $\ell > 0$, there are two positions $\pm y_m$ such that $\Theta(\pm y_m) = k$ and the polar motion corresponds to oscillations between $\pm y_m$;

* if $4\Delta_y(0)/\Delta_y''(0) < \ell < 0$, the movement depends on the value of $k$: either the equation $\Theta(y) = k$ has two solutions $\pm y_m$ and the motion corresponds to oscillations around the equatorial plane between these two values, or it has four solutions $\pm y_m$ and $\pm y_v$ (with $0 < y_v < y_m$) and we can have oscillations either in $[-y_m, -y_v]$ and $[y_v, y_m]$;

* if $\ell < 4\Delta_y(0)/\Delta_y''(0)$, the movement depends on the relative values of $y_*^\pm$ and the roots of $\Theta(\pm y_m) = k$ — one should note that this case does not exist in the Kerr limit because $y_*^\pm = \pm a$ in this case;

* if $\ell = 0$, the motion depends on the exact formulation of $\Delta_y$ through the third derivative of $\Theta$.

Let us recall that for the critical orbits relevant for the photon ring image, the parameters $(\ell_c, k_c)$ are already fixed in term of the radial position $r_0$ by the solution (2.23) for the spherical photon orbits. Having discussed the radial and polar motion of spherical photon orbits in the KOS family, we now turn to the black hole image and the derivation of the critical curve.
## 2.3 Screen angular coordinates and the critical curve

We now derive the projected image of the photon ring as seen on the screen of an asymptotic observer. To that end, we consider an observer at fixed spatial coordinates given by $(r_O, y_O, \varphi_O)$. We assume that $r_O \gg r_h$ while $y_O$ encodes the inclination of the observer w.r.t the black rotation axis, i.e. $y_O = 0$ corresponds to an asymptotic observer in the equatorial plane of the black hole.

We further assume that the observer has zero angular momentum. The orthogonal frame associated to this observer is given by

$$e_{(\tau)} = \sqrt{\frac{r^4 \Delta_y - y^4 \Delta_r}{\Sigma \Delta_r \Delta_y}} \left( \partial_{\tau} + \frac{y^2 \Delta_r + r^2 \Delta_y}{r^4 \Delta_y - y^4 \Delta_r} \partial_{\varphi} \right), \qquad (2.33a)$$

$$e_{(r)} = \sqrt{\frac{\Delta_r}{\Sigma}} \partial_r, \qquad (2.33b)$$

$$e_{(y)} = \sqrt{\frac{\Delta_y}{\Sigma}} \partial_y, \qquad (2.33c)$$

$$e_{(\varphi)} = \sqrt{\frac{\Sigma}{r^4 \Delta_y - y^4 \Delta_r}} \partial_{\varphi}. \qquad (2.33d)$$

which satisfies the properties that the observer's wordline is normal to the hypersurfaces of constant time $\tau$ and whose 4-velocity is orthogonal to $\partial_{\varphi}$. Now, we assume that this observer is equipped with a screen centered on the direction to the black hole such that $e_{(r)}$ is normal to the screen. In the large $r$ limit, this frame reduces to

$$e_{(\tau)} = \partial_{\tau}, \quad e_{(r)} = \partial_{r}, \quad e_{(y)} = \frac{\sqrt{\Delta_y(y_O)}}{r_O} \partial_y, \quad e_{(\varphi)} = \frac{1}{r_O \sqrt{\Delta_y(y_O)}} \partial_{\varphi}. \qquad (2.34)$$

such that $e_{(\tau)}$ is the unit vector orthogonal to constant-$\tau$ hypersurfaces.

Now let us consider the photon trajectories generating the image of the black hole on the observer's screen. The observer will collect photons whose trajectories satisfy

$$\mathcal{V}_r(r_O) = (r_O^2 - \ell)^2 - \Delta_r(r_O)k \ge 0, \qquad (2.35a)$$

$$\mathcal{V}_y(y_O) = \Delta_y(y_O)k - (y_O^2 + \ell)^2 \ge 0. \qquad (2.35b)$$

At large $r$, one has always $\mathcal{V}_r(r_O) \sim r_O^4$, which implies that (2.35a) is always satisfied. On the other hand, (2.35b) translates into a condition on the parameters $\ell$ and $k$ at a given observer's inclination $y_O$. A photon that follows a geodesic $\Gamma$ has the following 4-momentum at the location of $O$:

$$\mathbf{p} = p^{\tau} e_{(\tau)} + p^{r} e_{(r)} + \frac{p^{y} r_{O}}{\sqrt{\Delta_{y}(y_{O})}} e_{(y)} + p^{\varphi} r_{O} \sqrt{\Delta_{y}(y_{O})} e_{(\varphi)} = E(e_{(\tau)} + \mathbf{n}), \qquad (2.36)$$

with $\mathbf{n}$ a unit spacelike vector that describes the photon's direction of propagation, seen from the point of view of $O$, i.e.

$$\mathbf{n} = p^r e_{(r)} + \frac{p^y r_O}{\sqrt{\Delta_y(y_O)}} e_{(y)} + p^\varphi r_O \sqrt{\Delta_y(y_O)} e_{(\varphi)}. \qquad (2.37)$$
As seen from the observer, the projected momenta is given by $-\mathbf{n}$. The relevant information for the photon ring image is carried by the projected momenta of the photon on the observer's screen. Splitting the vector $-\mathbf{n}$ as $-\mathbf{n} = -\mathbf{n}_{\perp} - \mathbf{n}_{\parallel}$, one has

$$-\mathbf{n}_{\parallel} = -\frac{p^{\varphi} r_{\mathcal{O}}}{E \sqrt{\Delta_{y}(y_{\mathcal{O}})}} e_{(y)} - \frac{p^{\varphi} r_{\mathcal{O}} \sqrt{\Delta_{y}(y_{\mathcal{O}})}}{E} e_{(\varphi)}. \qquad (2.38)$$

Now, let us introduce dimensionless Cartesian coordinates $(\alpha, \beta)$ on the screen such that the black hole's rotation axis appears as the $\beta$ axis. We then have

$$-\mathbf{n}_{\parallel} = \alpha \mathbf{e}_{(\alpha)} + \beta \mathbf{e}_{(\beta)} \quad \text{with} \quad \mathbf{e}_{(\alpha)} = \mathbf{e}_{(\varphi)} \quad \text{and} \quad \mathbf{e}_{(\beta)} = -\mathbf{e}_{(y)}, \qquad (2.39)$$

such that comparing to (2.38), one finds the following relation between the couple $(\alpha, \beta)$ and the components of the projected photon momenta:

$$\alpha = -\frac{p^{\varphi} r_{\mathcal{O}} \sqrt{\Delta_{y}(y_{\mathcal{O}})}}{E} \quad \text{and} \quad \beta = \frac{p^{y} r_{\mathcal{O}}}{E \sqrt{\Delta_{y}(y_{\mathcal{O}})}}. \qquad (2.40)$$

Now, using (2.17d) and (2.17c) to express $p^y$ and $p^\varphi$, one obtains

$$\alpha = -\frac{r_{\mathcal{O}} \sqrt{\Delta_{y}(y_{\mathcal{O}})}}{E(r_{\mathcal{O}}^{2} + y_{\mathcal{O}}^{2})} \left( \frac{r_{\mathcal{O}}^{2} E - L}{\Delta_{r}(r_{\mathcal{O}})} + \frac{y_{\mathcal{O}}^{2} E + L}{\Delta_{y}(y_{\mathcal{O}})} \right), \quad \beta = \pm \frac{r_{\mathcal{O}}}{r_{\mathcal{O}}^{2} + y_{\mathcal{O}}^{2}} \left( k - \frac{(y_{\mathcal{O}}^{2} + \ell)^{2}}{\Delta_{y}(y_{\mathcal{O}})} \right)^{1/2}. \qquad (2.41)$$

Finally, taking the large $r_{\mathcal{O}}$ limit, the above relations reduce to

$$\alpha = -\frac{\sqrt{\Delta_{y}(y_{\mathcal{O}})}}{r_{\mathcal{O}}} \left( 1 + \frac{y_{\mathcal{O}}^{2} + \ell}{\Delta_{y}(y_{\mathcal{O}})} \right), \qquad (2.42a)$$

$$\beta = \pm \frac{1}{r_{\mathcal{O}}} \left( k - \frac{(y_{\mathcal{O}}^{2} + \ell)^{2}}{\Delta_{y}(y_{\mathcal{O}})} \right)^{1/2}. \qquad (2.42b)$$

These equations relates the coordinates $(\alpha, \beta)$ of the impact of the photon on the screen with its parameters $(\ell, k)$ and the position and inclination of the observer $(r_{\mathcal{O}}, y_{\mathcal{O}})$. Now, the photon ring is generated by critical photons orbits which have parameters $(\ell, k)$ infinitely close to the critical parameters $(\ell_c, k_c)$ characterizing spherical photon orbits.<fnref n="13" /> Therefore, these geodesics remain trapped in the region $r_0 \in [r_{ph}^+, r_{ph}^-]$ before escaping to infinity and eventually reaching the observer's screen. The black hole shadow is obtained as the image of these geodesics on a screen far away from the black hole. Therefore, using the expressions (2.23)

$$\mathcal{V}_r(r_p) = \mathcal{V}_r(r_a) = 0 \quad \text{and} \quad \mathcal{V}_r(r) \ge 0 \quad \text{when} \quad r \in [r_p, r_a].$$

For the Kerr metric, this situation in investigated in [108]: due to the roots system of the potential, no such nonspherical orbit is found. For the Kerr off-shell family, despite careful investigation, we found no expression of $\Delta_r$ that allows for such orbits. This is the case in particular for all examples studied in this work.

<footnote n="13">¹³Nonspherical closed photon orbits around the black hole could exist and also contribute to the photon ring. A nonspherical photon orbit exists between a perihelion $r_p$ and aphelion $r_a$ if and only if</footnote>
for the critical parameters ($\ell_c, k_c$) derived in section 2.2.1, the cartesian coordinates ($\alpha_c, \beta_c$) take the following form:

$$ \alpha_c = -\frac{1}{r_{\mathcal{O}}\sqrt{\Delta_y(y_{\mathcal{O}})}} \left[ \Sigma(r_0, y_{\mathcal{O}}) + \Delta_y(y_{\mathcal{O}}) - \frac{4r_0\Delta_r(r_0)}{\Delta'_r(r_0)} \right], \quad (2.43a) $$

$$ \beta_c = \pm \frac{1}{r_{\mathcal{O}}\sqrt{\Delta_y(y_{\mathcal{O}})}} \left[ -16r_0^2 \frac{\Delta_r(r_0)^2}{\Delta'_r(r_0)^2} + 8r_0 \frac{\Delta_r(r_0)}{\Delta'_r(r_0)^2} [2r_0\Delta_y(y_{\mathcal{O}}) + \Sigma(r_0, y_{\mathcal{O}})\Delta'_r(r_0)] - \Sigma^2(r_0, y_{\mathcal{O}}) \right]^{1/2} . \quad (2.43b) $$

With these expressions, one has access to the positions of the impact of the nearly critical photons building the critical curve. Notice that at this stage, the couple ($\alpha_c, \beta_c$) depends on the parameters of the black hole, the position of the observer's screen ($r_{\mathcal{O}}, y_{\mathcal{O}}$) which is fixed, and finally on the position $r_0$ at which the photon is emitted within the photon shell, i.e. the region spanned by spherical photon orbits. Contrary to the other parameters, one has $r_0 \in [r_{\text{ph}}^{-}, r_{\text{ph}}^{+}]$ such that this parameter is not fixed. Therefore, for now on, we write the couple of Cartesian coordinates as $\alpha_c(r_0)$ and $\beta_c(r_0)$ to make the dependency on this parameter explicit.

It is useful to introduce alternative coordinates on the screen of the observer. Let us decompose the dimensionless Cartesian coordinates ($\alpha(r_0), \beta(r_0)$) into polar ones as

$$ \alpha(r_0) := \mathcal{R}(r_0) \cos \phi(r_0), \quad \beta(r_0) := \mathcal{R}(r_0) \sin \phi(r_0), \quad (2.44) $$

where $\mathcal{R}(r_0)$ is a dimensionless radius and $\phi(r_0)$ the angle w.r.t the $\alpha$-axis. Notice that here, $r_0$ plays the role of a curvilinear abscissa for the critical curve. Let us introduce a basis of orthogonal unit vectors on the plane ($\alpha, \beta$) denoted ($\vec{e}_{\alpha}, \vec{e}_{\beta}$). The position on the curve is given by

$$ \vec{r}_c = \mathcal{R}_c(r_0)\vec{u} \quad \text{with} \quad \vec{u} = \cos \phi_c(r_0)\vec{e}_{\alpha} + \sin \phi_c(r_0)\vec{e}_{\beta}. \quad (2.45) $$

A direct computation shows that the coordinates ($\mathcal{R}_c(r_0), \phi_c(r_0)$) are related to the black hole parameters by

$$ \mathcal{R}_c(r_0) = \sqrt{\alpha_c^2(r_0) + \beta_c^2(r_0)} = \frac{1}{r_{\mathcal{O}}} \left[ \Delta_y(y_{\mathcal{O}}) + 2\Sigma(r_0, y_{\mathcal{O}}) + 8[2r_0 - \Delta'_r(r_0)] \frac{\Delta_r(r_0)}{\Delta'_r(r_0)^2} \right]^{1/2}, \quad (2.46a) $$

$$ \tan[\phi_c(r_0)] = \frac{\beta_c(r_0)}{\alpha_c(r_0)} = \mp \frac{\sqrt{16r_0^2 \Delta_r(r_0) \Delta_y(y_{\mathcal{O}}) - [\Sigma(r_0, y_{\mathcal{O}}) \Delta'_r(r_0) - 4r_0 \Delta_r(r_0)]^2}}{\Delta'_r(r_0)(\Sigma(r_0, y_{\mathcal{O}}) + \Delta_y(y_{\mathcal{O}})) - 4r_0 \Delta_r(r_0)}. \quad (2.46b) $$

These relations are not yet the ones relevant to interpret the interferometry signal. As discussed in detail in [54, 59], the relevant information is contained in the angle $\phi(r_0)$ which encodes the ratio between the rate of change of the curve w.r.t the $\alpha$-axis and the rate of change of the curve w.r.t the $\beta$-axis. Concretely, this angle is given by

$$ \tan[\phi_c(r_0)] = -\frac{\partial_{r_0} \alpha_c(r_0)}{\partial_{r_0} \beta_c(r_0)} = \pm \frac{\sqrt{16r_0^2 \Delta_r(r_0) \Delta_y(y_{\mathcal{O}}) - [\Delta'_r(r_0) \Sigma(r_0, y_{\mathcal{O}}) - 4r_0 \Delta_r(r_0)]^2}}{4r_0(\Delta_r(r_0) - \Delta_y(y_{\mathcal{O}})) - \Sigma(r_0, y_{\mathcal{O}}) \Delta'_r(r_0)}. \quad (2.47) $$

This is the key relation introduced in [59] to reconstruct the image of the critical curve.

This concludes the derivation of the critical curve and its properties for the general KOS family. We now use these general results to discuss the status of the test of the Kerr hypothesis using the shape of the critical curve.
## 3 Parameter estimations beyond Kerr: concrete examples

In this section, we use the previous results to address two key points.

First, we demonstrate that the formula obtained in the previous section allows us to efficiently capture the position and shape of the photon ring for different examples of Kerr-like objects belonging to the KOS family. Two of these examples are known models but the remaining two exhibit new types of corrections which have not appeared in the literature so far, illustrating the power of the new parametrization. Contrary to most investigation which only focus on the position of the equatorial radius of the photon ring, we derive the positions of the horizon, the radial and polar positions of the photon ring and the critical parameters ($l_c, k_c$) for spherical photon orbits for each examples.

Second, we use these concrete examples to discuss to what extend the shape of the photon ring can provide a sharp test of the Kerr hypothesis, and therefore of GR. In [59], it was proposed to fit the critical curve (2.47), denoted $f(\varphi)$, by a phoval, described by

$$f(\phi) = R_0 + \sqrt{R_1^2 \sin^2 \phi + R_2^2 \cos^2 \phi + (X - \chi) \cos \phi + \arcsin(\chi \cos \phi)}, \quad (3.1)$$

where $R_0, R_1, R_2, X$ and $\chi$ are scalar parameters. The link of each of these parameters with $M, a$ and the inclination $y_O$ can only be found numerically. They key point is that this phoval (or circplise) function fit to a great accuracy the Kerr critical curve. In the following, we confront this phoval fitting function introduced in [59] to the four examples and study the degeneracy in the estimation of the black hole parameters.

### 3.1 The Kerr black hole

As a first test of our general close formula, let us treat the special case of the Kerr black hole. Consider again the following change of coordinates

$$y = a \cos \theta, \quad \varphi = \frac{\phi}{a}, \quad \tau = t - a\phi, \quad (1.4)$$

which relates the coordinates $(\tau, r, y, \varphi)$ for the Boyer-Lindquist coordinates $(t, r, \theta, \phi)$. Then, the Kerr metric is recovered by choosing

$$\Delta_r^{\text{Kerr}} = r^2 + a^2 - 2Mr, \quad \Delta_y^{\text{Kerr}} = a^2 - y^2 = a^2 \sin^2 \theta. \quad (1.5)$$

The Kerr outer horizon is located at

$$r_h = M + \sqrt{M^2 - a^2}. \quad (3.2)$$

The parameters ($l_c, k_c$) for the spherical photon orbits of Kerr satisfy (2.23) which reads

$$k_c = \frac{4r_0^2(r_0^2 - 2Mr_0 + a^2)}{(r_0 - M)^2}, \quad l_c = -\frac{r_0(r_0^2 - 3Mr_0 + 2a^2)}{r_0 - M}, \quad (3.3)$$

where $r_0$ is a radius satisfying $r_0 > M$ and belonging to the photon ring. The range of $r_0$ is constrained by the inequality (2.25). Plugging the Kerr functions (1.5), this equation reduces to the standard cubic polynomial equation

$$r^3 - 6Mr^2 + 9M^2r - 4a^2M \le 0, \quad (3.4)$$
<figure><img src="image_3.png" /><figcaption>Figure 1. Critical curve of the Kerr black hole for $M = 1$, $y_O = 0$ and various values of $a$.</figcaption></figure>

which determines the minimal and maximal equatorial radii of the Kerr photon ring given by

$$r_{\text{ph}}^{\pm} = 2M \left[ 1 + \cos\left(\frac{2}{3}\xi^{\pm}\right) \right] \quad \text{with} \quad \xi^{\pm} = \arccos\left(\pm\frac{a}{M}\right). \qquad (3.5)$$

such that $r_0 \in [r_{\text{ph}}^{+}, r_{\text{ph}}^{-}]$. As we have seen, the polar motion depends on the value of the couple $(\ell, k)$. Besides the equilibrium positions $y_e^{\pm} = \{0, \pm\ell\}$, one has also other positions described in section (2.2.2).

Therefore, the closed formula obtained in the previous section consistently reproduce the known results for the positions of the photon ring and the parameters of the spherical photon orbits for the Kerr black hole. Using the above results, we can plot the parametric form of the critical curve (2.47) for different values of the spin $a$. It is shown in figure 1.

Having tested our formula on the standard Kerr black hole example, our goal now is to consider different examples of Kerr-like objects and reproduce the very same study for each of them. However, before presenting these Kerr-like objects, let us already make an important remark. It is standard to work with dimensionless screen coordinates $(r_O\alpha/r_*, r_O\beta/r_*)$ where $r_O$ is the radial position of the observer and $r_*$ is usually taken to be the mass of the black hole (see figure 1). Therefore, the mass of the black hole will directly affect the overall scaling of the image. However, while the mass of the Kerr black hole admits a well-defined interpretation as an asymptotic charge of the geometry, this is no longer true in modified gravity. Indeed the standard Komar charge one uses in GR is generically modified when one considers alternatives (and in particular higher order) lagrangians. While the modifications to the Komar charge can be computed in principle using the covariant phase space method, this point has not been studied systematically such that given an ad hoc Kerr-like geometry, there is no way to know what one should interprete what is the quasi-local asymptotic mass of the object as this notion depends on the theory which admits this geometry as a solution. In order to bypass this difficulty when dealing with ad hoc models for Kerr-like object, one way is to use the position of the horizon of the geometry which can be extracted solely form the metric in order to rescale the screen coordinates. This is the choice we adopt in the following when plotting the shadow of Kerr like objects.
## 3.2 Radial deformations to Kerr

In this section, we consider three examples of deviations to Kerr which are purely radial, i.e. they correspond to different choices of $\Delta_r(r)$. The first two are known examples of Kerr-like objects that have been studied in the literature: the Kerr-MOG black hole derived in [109] and the regular rotating black hole introduced by Simpson and Visser in [110]. The third example is a new type of correction that does not appear elsewhere.

### 3.2.1 The Kerr-MOG black hole

Consider first The Kerr-MOG black hole discussed in [109]. This rotating black hole geometry is a solution of the so called Scalar-Vector-Tensor (SVT) theory introduced in [111]. Its shadow was previously studied in [112]. The geometry is characterized by a new parameter $\alpha$ which can be understood as a coupling constant of the SVT lagrangian which changes the Newton constant $G_N$ to effective Newton constant of the form $G = (1 + \alpha)G_N$. It can be shown that the new Kerr-like solution of this theory enjoys a modified radial function $\Delta_r(r)$ given by:

$$ \Delta_r(r) = r^2 - 2GMr + a^2 + \alpha G_N GM^2 = \Delta^{\text{Kerr}} + \alpha G_N GM^2, \quad (3.6) $$

where $G = G_N(1 + \alpha) = 1 + \alpha$ when $G_N = 1$. The outer horizon is located at

$$ r_h = (1 + \alpha)M \left( 1 + \sqrt{1 - \frac{a^2}{(1 + \alpha)^2 M^2} - \frac{\alpha}{1 + \alpha}} \right) \quad (3.7) $$

By construction, this Kerr-like geometry belongs to the KOS family and therefore, one can use the previous results to obtain analytical expressions for the different locus of points of its photon rings and the parameters of its spherical photon orbits. The couple $(\ell_c, k_c)$ is given by

$$ \ell_c = - \frac{r(2a^2 + 2\alpha(\alpha + 1)M^2 - 3(\alpha + 1)Mr + r^2)}{r - (\alpha + 1)M}, \quad (3.8a) $$

$$ k_c = \frac{4r^2 (a^2 + \alpha(\alpha + 1)M^2 - 2(\alpha + 1)Mr + r^2)}{(\alpha M + M - r)^2}. \quad (3.8b) $$

The parameter $r_0$ ranges between the minimal and maximal positions (written at first order in $\alpha$ for clarity here)

$$ r_{\text{ph}}^{\pm} = 2M \left[ 1 + \cos\left(\frac{2}{3}\xi^{\pm}\right) \right] + \frac{1}{3}\alpha M \left[ 7 + \frac{2a^2 \cos\left(\frac{2}{3}\xi^{\pm}\right)}{a^2 + M^2\left(\cos\left(\frac{2}{3}\xi^{\pm}\right) + \cos\left(\frac{4}{3}\xi^{\pm}\right)\right)} \right] + \mathcal{O}(\alpha^2), \quad (3.9) $$

where $\xi^{\pm}$ is the quantity defined for Kerr in (3.5). Notice that since the geometry does not have any polar deformation, i.e. $\Delta_y(y) = \Delta_y^{\text{Kerr}}(y)$, the expressions for the equilibrium polar positions $y_e^{\pm}$ have the same expressions as for Kerr.

With these results, one can use the general formula (2.47) to plot the critical curve of the Kerr-MOG black hole for i) four different value of the new parameter $\alpha$ and ii) two different values of the inclination $y_o$. It is represented in figure 2.

As a first remark, one can see that the deviations w.r.t the Kerr shadow (depicted in green) are qualitatively more important at high inclination, i.e. for $y_o = 0$, in agreement
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Critical curve of the Kerr-MOG black hole for $M = 1$, $a = 0.9$ and different values of $y_O$.</figcaption></figure>

<table><caption>**Table 1.** Example of degeneracy between the Kerr black hole and the Kerr-MOG model.</caption><thead><tr><th>Parameters</th><th><i>M</i></th><th><i>a</i></th><th>&alpha;</th><th><i>y</i><sub>O</sub></th><th><i>R</i><sub>0</sub></th><th><i>R</i><sub>1</sub></th><th><i>R</i><sub>2</sub></th><th>Residuals</th></tr></thead><tbody><tr><td>Kerr</td><td>1</td><td>0.2</td><td></td><td>0</td><td>10.0</td><td>0.394</td><td>0.370</td><td>1.88 &times; 10<sup>-5</sup></td></tr><tr><td>MOG</td><td>0.922</td><td>0.193</td><td>0.101</td><td>0</td><td>10.0</td><td>0.395</td><td>0.371</td><td>1.86 &times; 10<sup>-5</sup></td></tr></tbody></table>

with previous studies [41]. However, the shape of the critical curve represented here does not provide in itself a precise way to compare the Kerr and Kerr-MOG black holes. To confront the two models, one needs to perform a parameter estimations of the black holes parameters.

To that end, one can use the fitting function introduced by Gralla and Lupsasca and which consists in the phoval (or circlipse) defined by the equation (3.1) to test the two models. This is summarized in the table 1. As a result, we find that this fitting function admits a non-trivial degeneracy in the estimation of the parameters ($R_0, R_1, R_2$). Concretely, for this critical curve, the circlipse fitting test cannot distinguish between a Kerr black hole geometry with mass and spin ($M, a$) = (1, 0.2) and a Kerr-MOG black hole with parameters ($M, a$) = (0.922, 1.923) and $\alpha$ = 0.101. This example provides the first explicit and concrete proof that the circlipse test does not, by itself, provide a sharp test of the Kerr hypothesis. It suggests that in order to properly test the Kerr hypothesis, one has to be able to test the relation between the mass and angular momentum which characterize a given geometry. Obtaining independent measures of the mass and the spin of the black holes appear therefore as a mandatory task to break the degeneracy presented here. We stress that this example (and the analytical derivation of the critical curve used here) could not have been obtained without the use of the KOS parametrization introduced in this work. Let us now turn to other examples of Kerr-like objects.

### 3.2.2 The regular Simpson-Visser rotating model

We now consider a model of rotating regular black hole. Such models have been studied intensively in the literature [113–115]. To a large extent, the regularization of the interior singularity of spherically symmetric black holes have been performed by introducing an
effective de Sitter and thus repulsive core. Yet, other type of regularization have been studied and it has been shown that one could obtain interesting ad hoc regular black hole geometries by replace the deep interior geometry by a Minkowski core [116]. This approach was extended to rotating black hole and an explicit regular Kerr-like black hole dubbed the “Eyes of the Storm” (EOS) was presented in [110]. The modification amounts at the a change of the mass which becomes $M(r) = Me^{-\ell M/r}$ where $\ell$ is the new UV cut-off encoding the scale at which the effective quantum effects are triggered. For a given $(M, a)$, $\ell = 0$ yields the Kerr black hole. Furthermore, possible values of $\ell$ are bounded from above: $0 \le \ell \le l_{extr}$. Image of this regular black hole have already been studied in [117]

The modification amounts at a change of $\Delta(r)$ of the form

$$ \Delta_{\text{EOS}} = r^2 + a^2 - 2rM(r) = r^2 + a^2 - 2rMe^{-\ell M/r}, \qquad (3.10) $$

belonging thus to the KOS family. Unfortunately, the choice of the exponential function in the function $M(r)$ does not allows one to find closed analytical formulae for the positions of the horizons. At best, one can provide an expression for the correction to the Kerr horizon when $\ell$ is considered small, which gives

$$ r_h^{\text{EOS}} = M + \sqrt{M^2 - a^2 - 2M\ell + \mathcal{O}(\ell^2)} < r_h^{\text{Kerr}} \qquad (3.11) $$

Using the previous relations, we now derive the parameters $(\ell_c, k_c)$ for the critical spherical photon orbits which read

$$ \ell_c = \frac{r^2 \left( M(3r - M\ell) - (2a^2 + r^2) e^{\frac{M\ell}{r}} \right)}{r^2 e^{\frac{M\ell}{r}} - M(M\ell + r)} \qquad (3.12) $$

$$ k_c = \frac{16r^2 \left( a^2 + r \left( r - 2Me^{-\frac{M\ell}{r}} \right) \right)}{\left( 2r - \frac{2Me^{-\frac{M\ell}{r}}(M\ell+r)}{r} \right)^2} \qquad (3.13) $$

where the position of the photon ring is given (at first order in $\ell$ for clarity) by

$$ r_{\text{ph}}^{\pm} = 2M \left[ 1 + \cos\left(\frac{2}{3}\xi^{\pm}\right) \right] + \frac{2}{3}\ell M \left[ -2 + \frac{a^2}{a^2 + M^2\left(\cos\left(\frac{2}{3}\xi^{\pm}\right) + \cos\left(\frac{4}{3}\xi^{\pm}\right)\right)} \right] + \mathcal{O}(\ell^2), \quad (3.14) $$

where $\xi^{\pm}$ is the quantity defined for Kerr in (3.5). The critical curve of the EOS black hole for a few different values of the parameter $\ell$ and at two different inclinations are shown in figure 3. While the green curve represents the Kerr black hole critical curve, we vary the parameter $\ell$ up to its nearly maximal value corresponding to the extremality. Following the same procedure as above, we also find that the circplispe fitting function suffers from a degeneracy as shown by the concrete example of parameters estimation presented in table 2. This second example provides thus a confirmation of the result obtained for the study of the Kerr-MOG model. We can now turn to yet new examples of Kerr-like objects which have not been treated so far.

### 3.2.3 Log corrections

An interesting outcome of the KOS parametrization used in this work is to allow us to consider new ad hoc modification to the Kerr geometry which will not break the Killing tower.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Critical curve of the Simpson-Visser black hole for $M = 1$, $a = 0.9$ and different values of $y_O$.</figcaption></figure>

<table><caption>**Table 2.** Example of degeneracy between the Kerr black hole and the Simpson-Visser model. The observer is at $y_O = 0$ in both cases.</caption><thead><tr><th>Parameters</th><th><i>M</i></th><th><i>a</i></th><th><i>ℓ</i></th><th><i>y</i><sub>O</sub></th><th><i>R</i><sub>0</sub></th><th><i>R</i><sub>1</sub></th><th><i>R</i><sub>2</sub></th><th>Residuals</th></tr></thead><tbody><tr><td>Kerr</td><td>1</td><td>0.2</td><td>0</td><td>0</td><td>10.0</td><td>0.394</td><td>0.370</td><td>1.88 × 10<sup>-5</sup></td></tr><tr><td>SV</td><td>1.32</td><td>0.1</td><td>0.82</td><td>0</td><td>10.0</td><td>0.402</td><td>0.380</td><td>1.53 × 10<sup>-5</sup></td></tr></tbody></table>

As a concrete illustrative example, we now consider a one-parameter family of deviations from Kerr defined by

$$ \Delta_r = r^2 - 2Mr + a^2 + qM^2 \log\left(\frac{r}{M}\right), \qquad (3.15) $$

This log-modification still preserves the asymptotical flatness condition as $\Sigma \sim r^2$ at large $r$. However, to our knowledge, it has not been considered in the literature. We write analytical formulas at first order in the dimensionless variable $q$ for clarity; however, all results are exact in $q$, in particular the plots shown on figure 4. The horizon $r_h$ of the Kerr black hole is changed at first order as follows:

$$ r_h = M + \sqrt{M^2 - a^2} - q \frac{M}{2} \frac{\log\left(1 + \sqrt{1 - (a/M)^2}\right)}{\sqrt{1 - (a/M)^2}} + \mathcal{O}(q^2). \qquad (3.16) $$

A direct computation shows that the couple $(\ell_c, k_c)$ is given by

$$ \ell_c = r \left( r - \frac{4r \left(a^2 + M^2 q \log\left(\frac{r}{M}\right) + r(r - 2M)\right)}{M^2 q - 2Mr + 2r^2} \right), \qquad (3.17) $$

$$ k_c = \frac{16r^4 \left(a^2 + M^2 q \log\left(\frac{r}{M}\right) + r(r - 2M)\right)}{\left(M^2 q - 2Mr + 2r^2\right)^2}. \qquad (3.18) $$

The parameter $r_0$ ranges between the minimal and maximal positions given by

$$ r_{\text{ph}}^{\pm} = 2M \left[ 1 + \cos\left(\frac{2}{3}\xi^{\pm}\right) \right] + qd_{\pm}^{(1)} + \mathcal{O}(q^2), \quad \xi^{\pm} = \arccos\left(\pm\frac{a}{M}\right). \qquad (3.19) $$
<figure><img src="image_6.png" /><figcaption>Figure 4. Critical curve of the log-corrected black hole for $M = 1$, $a = 0.9$ and different values of $y_{\mathcal{O}}$.</figcaption></figure>

The corrections $d_{\pm}^{(1)}$ are such that

$$d_{\pm}^{(1)} = -\frac{M(4\log[4\cos^2(\xi^{\pm}/3)]-1)(\cos(4\xi^{\pm}/3)+\cos(2\xi^{\pm}))(2\log[4\cos^2(\xi^{\pm}/3)]-1)(a/M)^2}{2+2\cos(2\xi^{\pm})+3\cos(4\xi^{\pm}/3)+3\cos(2\xi^{\pm}/3)-(a/M)^2} \quad (3.20)$$

This first example shows that one can use the symmetries of the Kerr off shell system to analytically capture the modifications of the photon ring for some new deviations w.r.t. Kerr. Using the results of 2.3, we can numerically trace the shadow of the black hole. The result is shown on figure 4.

### 3.3 Polar deformations to Kerr

We now turn to a situation in which $\Delta_y$ deviates from the Kerr solution while $\Delta_r$ is unchanged. In order to maintain the reflection symmetry, we take the parametrization

$$\Delta_y = a^2 - y^2 + py^4. \quad (3.21)$$

Similarly to section 3.2.3, we present results at first order in $p$ for clarity but our results are exact in this variable. The horizon $r_h$ is unchanged since $\Delta_r$ is the same, as well as $\ell_c$, $k_c$ and $r_{ph}^{\pm}$ which are the same as in the Kerr case (see section 3.1). However the value of $y_m$, which corresponds to the maximal value of $y$ for a null geodesic with $\ell > 0$, is modified: one has

$$y_m = \sqrt{\frac{1}{2} \left( -(k+2\ell) + \sqrt{k} \sqrt{k+4(a^2+\ell)} \right) + qy_m^{(1)} + \mathcal{O}(q^2)}, \quad (3.22)$$

with

$$y_m^{(1)} = - \frac{(k - \sqrt{k} \sqrt{4a^2 + k + 4\ell})(\sqrt{k} \sqrt{4a^2 + k + 4\ell} - k - 2\ell)^{3/2}}{4\sqrt{2}a^2(-\sqrt{k} \sqrt{4a^2 + k + 4\ell} + 4a^2 + k + 4\ell)}. \quad (3.23)$$

Using the results of section 2.3, we can numerically trace the shadow of the black hole. The result is shown on figure 5.
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Critical curve of the polar-corrected black hole for $M = 1$, $a = 0.9$ and different values of $y_{\odot}$. When $y_{\odot} = 0$, the black hole is seen from its equatorial plane, and the polar corrections do not come into play: the critical curves for all values of $p$ are the same.</figcaption></figure>

## 4 Discussion

In this work, we have shown how the Kerr off shell (KOS) family of geometries introduced in [75–77] allows one to efficiently parametrize beyond Kerr geometries. Using this new parametrization, we have studied the deformations of the shape of the black hole photon ring image and confronted the ability to distinguish between Kerr and beyond Kerr geometries using the circlipse fitting function proposed by Gralla and Lupsasca. Let us summarize the main points of our study.

**The KOS parametrization.** The KOS family of metrics stands as the most general set of geometries preserving the Killing tower. This structure uncovered in [75] encodes the fundamental explicit and hidden symmetries of the Kerr geometry. In particular, any KOS representant is of Petrov type D and is equipped with a rank-2 closed conformal Killing-Yano tensor and its dual. In turn, this ensures that the geometry also admits a Killing tensor which makes the geodesic motion integrable. Remarkably, despite preserving these hidden symmetries, the KOS parametrization still allows for a large set of deformations to the Kerr geometry encoded in the two free functions $\Delta_r(r)$ and $\Delta_y(y)$. On one hand, $\Delta_r(r)$ allows for general radial deformations and coincides with the restricted class I family of beyond Kerr geometries studied in [71] where $\Delta_r(r) = A_5(r)$. This radial deformation allows one to reproduce a large class of *ad hoc* models studied so far in the literature. On the other hand, the KOS parametrization provides another type of polar deformations encoded in $\Delta_y(y)$ which have not been studied so far. By construction, such polar deformations preserve the symmetries of Kerr, allowing for the analytic exploration of new Kerr-like objects such as the example studied in section 3.3.

**Analytic characterization of the photon shell.** Since the image of the black hole photon ring is generated by the geodesics which are nearly critical, i.e. geodesics which are infinitely close to the unstable critical spherical photon orbits, we have provided a detailed study of
this family of geodesics. One of the remarkable outcomes of the KOS family is to allow us to obtain an analytical characterization without fixing the free functions of the metric. Concretely, for an arbitrary KOS geometry, the radial and polar motions are determined by the two potentials (2.18a)–(2.18b) which depend respectively on the radial and polar deformations $\Delta_r(r)$ and $\Delta_y(y)$. The critical spherical photon orbits which satisfy (2.2.1) are characterized by a reduced angular momentum and Carter constant ($\ell_c, k_c$) given by (2.23). These expressions depend solely on the function $\Delta_r(r)$ and the radial equilibrium position $r_0$ which can be obtained only once a given choice of $\Delta_r(r)$ is made. The extent of the region in which the spherical orbits are allowed, i.e. the photon shell, is constrained by the inequality (2.25). For all the examples studied in this work, the allowed position $r_0$ is restricted to a range $[r_{ph}, r_{ph}^+]$ with $r_{ph}^- > r_h$ where $r_h > 0$ stands for the position of the outer horizon. One the other hand, we have provided a detailed classification of the possible polar equilibrium positions for such critical spherical photon orbits. Therefore, for a given KOS representant, the above formula allows one to fully characterize the photon shell geometry.

**Non-spherical unstable critical photon orbits.** Notice that a Kerr-like geometry could in principle allow for non-spherical orbits within the photon shell, i.e. photon orbits with $V'_r(r) \neq 0$ which admit radial turning points $r_*$ (at which $V_r(r_*) = 0$). The existence of such orbits depends on the form of the function $\Delta_r(r)$. For the Kerr geometry, one can show that such non-spherical critical orbits do not exist. For the general KOS geometry, we have not been able to provide a definite proof that they are forbidden. However, for all the examples considered here, one can show that such non-spherical critical orbit are not allowed.

**Analytic parametrization for the critical curve.** The main outcome of this study is the analytical derivation of the critical curve for an arbitrary KOS geometry given by eq. (2.47). This result provides a ready-to-use parametrization of the critical curve of beyond Kerr geometries which depends explicitly on the two free functions $\Delta_r(r)$ and $\Delta_y(y)$. To our knowledge, such analytic expressions has not been obtained yet for such a general class of rotating compact objects. At this stage, it is important to underline that the critical curve is not an observable *per se*. The interest in having an analytic prediction for this object lies in that fact that the photon subrings converge quickly towards this theoretical curve, making it the key object to parametrize. Let us emphasize that such an analytic result is made possible by the integrability of the geodesic motion and thus by the preservation of the Kerr symmetries.

**The circlipse fitting function and its degeneracy.** Equipped with this analytical parametrization of the critical curve, we have confronted the circlipse fitting functions on four different examples of Kerr-like geometries. We show that at maximal inclination, i.e. $y_O = 0$, an accurate fit of the Kerr black hole for a given $(M, a)$ can always be degenerate with another equally accurate fit for the Kerr-like object for another set of parameters $(M, a, \alpha)$ where $\alpha$ is the new parameter encoding the deviation to Kerr. Therefore, alone, the circlipse does not allow for a test the Kerr hypothesis, confirming recent results [41]. In order to break this degeneracy, an independent measurement of the mass and spin of the compact object is needed (see [118, 119]). Notice that this degeneracy between different models of
compact objects add up to other type of degeneracy, in particular degeneracy coming from the astrophysical modeling of the disk accretion [50]. It follows that in order to use black hole imaging as a sharp test of GR, important challenges are ahead of us to systematically alleviate the different degeneracies.

**Kinematical versus dynamical tests.** Another interesting outcome of working out the critical curve and the resulting image of beyond Kerr geometries is to underline some conceptual difficulties when dealing with the notion of mass and spin. Indeed, for the Kerr black hole which solve the Einstein equations, the notions of mass and spin are well-understood in term of asymptotic Komar charges whose expressions are derived from the Einstein-Hilbert action. However, when considering an *ad hoc* geometry which is not solution of any given theory of gravity, such as the Simpson-Visser model studied here, we have no recipe to identify a well defined notion of quasi-local mass or angular momentum. For this reason, while considering *ad hoc* geometries provide a useful first step to study beyond Kerr phenomenology, it is mandatory to associate such models to given theories to interpret the different parameters of the geometry as well-defined asymptotic quasi-local charges.<fnref n="14" /> This task represents another challenge in order to use black hole imaging to test alternatives theories of gravity such as scalar-tensor theories. (See [121] for review of rotating black holes and possible solution generating method in such modified theories). An interesting example treated here is the Kerr-MOG black hole, which is a *bona fide* solution of a scalar-vector-tensor theory where the additional parameter $\alpha$ plays the role of a coupling constant in the lagrangian of the theory. Thus, constraining the parameter $\alpha$ allows one to confront two theories and not only two geometries. From this perspective, the natural next step will be to identify suitable family of theories admitting the KOS geometries as exact solutions. This will be presented in the near future.

**Acknowledgments**

We thank Frédéric Vincent for useful discussions. EG acknowledges funding by l'Agence Nationale de la Recherche, projects StronG ANR-22-CE31-0015-01 and Einstein PPF ANR-23-CE40-0010-02. JBA and HR thank the Perimeter Institute for Theoretical Physics for its hospitality at the beginning of this project.

**References**

[1] EVENT HORIZON TELESCOPE collaboration, *First M87 Event Horizon Telescope results. I. The shadow of the supermassive black hole*, Astrophys. J. Lett. 875 (2019) L1 [arXiv:1906.11238] [INSPIRE].

[2] EVENT HORIZON TELESCOPE collaboration, *First M87 Event Horizon Telescope results. IV. Imaging the central supermassive black hole*, Astrophys. J. Lett. 875 (2019) L4 [arXiv:1906.11241] [INSPIRE].

[3] EVENT HORIZON TELESCOPE collaboration, *First M87 Event Horizon Telescope results. V. Physical origin of the asymmetric ring*, Astrophys. J. Lett. 875 (2019) L5 [arXiv:1906.11242] [INSPIRE].

<footnote n="14">¹⁴See [120] for a study of quasi-local energy using the covariant phase space in Horndeski.</footnote>
[4] EVENT HORIZON TELESCOPE collaboration, *First Sagittarius A* Event Horizon Telescope results. I. The shadow of the supermassive black hole in the center of the Milky Way*, *Astrophys. J. Lett.* **930** (2022) L12 [arXiv:2311.08680] [INSPIRE].

[5] EVENT HORIZON TELESCOPE collaboration, *First Sagittarius A* Event Horizon Telescope results. III. Imaging of the galactic center supermassive black hole, *Astrophys. J. Lett.* **930** (2022) L14 [arXiv:2311.09479] [INSPIRE].

[6] P.V.P. Cunha, C.A.R. Herdeiro, E. Radu and H.F. Runarsson, *Shadows of Kerr black holes with scalar hair*, *Phys. Rev. Lett.* **115** (2015) 211102 [arXiv:1509.00021] [INSPIRE].

[7] J.W. Moffat, *Modified gravity black holes and their observable shadows*, *Eur. Phys. J. C* **75** (2015) 130 [arXiv:1502.01677] [INSPIRE].

[8] P.V.P. Cunha, C.A.R. Herdeiro, E. Radu and H.F. Runarsson, *Shadows of Kerr black holes with and without scalar hair*, *Int. J. Mod. Phys. D* **25** (2016) 1641021 [arXiv:1605.08293] [INSPIRE].

[9] F.H. Vincent, E. Gourgoulhon, C. Herdeiro and E. Radu, *Astrophysical imaging of Kerr black holes with scalar hair*, *Phys. Rev. D* **94** (2016) 084045 [arXiv:1606.04246] [INSPIRE].

[10] P.V.P. Cunha and C.A.R. Herdeiro, *Shadows and strong gravitational lensing: a brief review*, *Gen. Rel. Grav.* **50** (2018) 42 [arXiv:1801.00860] [INSPIRE].

[11] P.V.P. Cunha, C.A.R. Herdeiro and E. Radu, *EHT constraint on the ultralight scalar hair of the M87 supermassive black hole*, *Universe* **5** (2019) 220 [arXiv:1909.08039] [INSPIRE].

[12] M. Khodadi, A. Allahyari, S. Vagnozzi and D.F. Mota, *Black holes with scalar hair in light of the Event Horizon Telescope*, *JCAP* **09** (2020) 026 [arXiv:2005.05992] [INSPIRE].

[13] J. Ben Achour et al., *On rotating black holes in DHOST theories*, *JCAP* **11** (2020) 001 [arXiv:2006.07245] [INSPIRE].

[14] F. Long, S. Chen, M. Wang and J. Jing, *Shadow of a disformal Kerr black hole in quadratic degenerate higher-order scalar-tensor theories*, *Eur. Phys. J. C* **80** (2020) 1180 [arXiv:2009.07508] [INSPIRE].

[15] Q. Gan, P. Wang, H. Wu and H. Yang, *Photon ring and observational appearance of a hairy black hole*, *Phys. Rev. D* **104** (2021) 044049 [arXiv:2105.11770] [INSPIRE].

[16] J. Kumar, S.U. Islam and S.G. Ghosh, *Investigating strong gravitational lensing effects by supermassive black holes with Horndeski gravity*, *Eur. Phys. J. C* **82** (2022) 443 [arXiv:2109.04450] [INSPIRE].

[17] M. Afrin, R. Kumar and S.G. Ghosh, *Parameter estimation of hairy Kerr black holes from its shadow and constraints from M87*, *Mon. Not. Roy. Astron. Soc.* **504** (2021) 5927 [arXiv:2103.11417] [INSPIRE].

[18] G. Antoniou, A. Papageorgiou and P. Kanti, *Probing modified gravity theories with scalar fields using black-hole images*, *Universe* **9** (2023) 147 [arXiv:2210.17533] [INSPIRE].

[19] X.-J. Wang et al., *Rings and images of Horndeski hairy black hole illuminated by various thin accretions*, *Phys. Rev. D* **107** (2023) 124052 [arXiv:2304.10015] [INSPIRE].

[20] B. Rodríguez, J. Chagoya and C. Ortiz, *Shadows of black holes in dynamical Chern-Simons modified gravity*, *arXiv:2403.13062* [INSPIRE].

[21] T. Ohgami and N. Sakai, *Wormhole shadows*, *Phys. Rev. D* **91** (2015) 124020 [arXiv:1704.07065] [INSPIRE].

[22] R. Shaikh, *Shadows of rotating wormholes*, *Phys. Rev. D* **98** (2018) 024044 [arXiv:1803.11422] [INSPIRE].
[23] F. Lamy, E. Gourgoulhon, T. Paumard and F.H. Vincent, *Imaging a non-singular rotating black hole at the center of the galaxy*, Class. Quant. Grav. **35** (2018) 115009 [arXiv:1802.01635] [INSPIRE].

[24] F. Rahaman et al., *Shadows of Lorentzian traversable wormholes*, Class. Quant. Grav. **38** (2021) 215007 [arXiv:2108.09930] [INSPIRE].

[25] F.H. Vincent et al., *Imaging a boson star at the galactic center*, Class. Quant. Grav. **33** (2016) 105015 [arXiv:1510.04170] [INSPIRE].

[26] P. Grandclément, *Light rings and light points of boson stars*, Phys. Rev. D **95** (2017) 084011 [arXiv:1612.07507] [INSPIRE].

[27] C.-Y. Chen, *Rotating black holes without $Z_2$ symmetry and their shadow images*, JCAP **05** (2020) 040 [arXiv:2004.01440] [INSPIRE].

[28] N. Tsukamoto, *Black hole shadow in an asymptotically-flat, stationary, and axisymmetric spacetime: the Kerr-Newman and rotating regular black holes*, Phys. Rev. D **97** (2018) 064021 [arXiv:1708.07427] [INSPIRE].

[29] Y.-H. Shu and J.-H. Huang, *Circular orbits and thin accretion disk around a quantum corrected black hole*, Phys. Lett. B **864** (2025) 139411 [arXiv:2412.05670] [INSPIRE].

[30] R. Ghosh, S. Sk and S. Sarkar, *Hairy black holes: nonexistence of short hairs and a bound on the light ring size*, Phys. Rev. D **108** (2023) L041501 [arXiv:2306.14193] [INSPIRE].

[31] R. Ghosh and C. Singha, *Can rotating black holes have short hairs?, Phys. Rev. D **111** (2025) 044008 [arXiv:2501.08767] [INSPIRE].

[32] N. Tsukamoto, Z. Li and C. Bambi, *Constraining the spin and the deformation parameters from the black hole shadow*, JCAP **06** (2014) 043 [arXiv:1403.0371] [INSPIRE].

[33] M. Khodadi and E.N. Saridakis, *Einstein-Aether gravity in the light of event horizon telescope observations of M87*, Phys. Dark Univ. **32** (2021) 100835 [arXiv:2012.05186] [INSPIRE].

[34] S.E. Gralla, *Can the EHT M87 results be used to test general relativity?, Phys. Rev. D **103** (2021) 024023 [arXiv:2010.08557] [INSPIRE].

[35] F.H. Vincent et al., *Geometric modeling of M87 as a Kerr black hole or a non-Kerr compact object*, Astron. Astrophys. **646** (2021) A37 [arXiv:2002.09226] [INSPIRE].

[36] M. Khodadi, G. Lambiase and D.F. Mota, *No-hair theorem in the wake of Event Horizon Telescope*, JCAP **09** (2021) 028 [arXiv:2107.00834] [INSPIRE].

[37] K. Glampedakis and G. Pappas, *Can supermassive black hole shadows test the Kerr metric?, Phys. Rev. D **104** (2021) L081503 [arXiv:2102.13573] [INSPIRE].

[38] A. Stepanian, S. Khlghatyan and V.G. Gurzadyan, *Black hole shadow to probe modified gravity*, Eur. Phys. J. Plus **136** (2021) 127 [arXiv:2101.08261] [INSPIRE].

[39] S. Vagnozzi et al., *Horizon-scale tests of gravity theories and fundamental physics from the Event Horizon Telescope image of Sagittarius A*, Class. Quant. Grav. **40** (2023) 165007 [arXiv:2205.07787] [INSPIRE].

[40] A.E. Broderick, K. Salehi and B. Georgiev, *Shadow implications: what does measuring the photon ring imply for gravity?, Astrophys. J. **958** (2023) 114 [arXiv:2307.15120] [INSPIRE].

[41] S. Staelens et al., *Black hole photon rings beyond general relativity*, Phys. Rev. D **107** (2023) 124026 [arXiv:2303.02111] [INSPIRE].

[42] M. Khodadi, S. Vagnozzi and J.T. Firouzjacee, *Event Horizon Telescope observations exclude compact objects in baseline mimetic gravity*, Sci. Rep. **14** (2024) 26932 [arXiv:2408.03241] [INSPIRE].
[43] Q. Yue, Z. Xu and M. Tang, Testing extended theories of gravity via black hole photon rings, Eur. Phys. J. C 85 (2025) 906 [arXiv:2505.21314] [INSPIRE].

[44] R. Carballo-Rubio, H. Delaporte, A. Eichhorn and P.G.S. Fernandes, Non-minimal light-curvature couplings and black-hole imaging, arXiv:2505.21431 [INSPIRE].

[45] C.-Y. Wang, D.-S. Lee and C.-Y. Lin, Null and timelike geodesics in the Kerr-Newman black hole exterior, Phys. Rev. D 106 (2022) 084048 [arXiv:2208.11906] [INSPIRE].

[46] B.-R. Chen, T. Hsieh and D.-S. Lee, Null geodesics in extremal Kerr-Newman black holes, Phys. Rev. D 111 (2025) 024058 [arXiv:2411.00458] [INSPIRE].

[47] S.E. Gralla, D.E. Holz and R.M. Wald, Black hole shadows, photon rings, and lensing rings, Phys. Rev. D 100 (2019) 024018 [arXiv:1906.00873] [INSPIRE].

[48] J.M. Bardeen, Timelike and null geodesics in the Kerr metric, in Proceedings, Ecole d'Eté de Physique Théorique: Les Astres Occlus: Les Houches, France, August, 1972, 215–240, (1973), p. 215 [INSPIRE].

[49] M.D. Johnson et al., Universal interferometric signatures of a black hole’s photon ring, Sci. Adv. 6 (2020) eaaz1310 [arXiv:1907.04329] [INSPIRE].

[50] H. Paugnat, A. Lupsasca, F. Vincent and M. Wielgus, Photon ring test of the Kerr hypothesis: variation in the ring shape, Astron. Astrophys. 668 (2022) A11 [arXiv:2206.02781] [INSPIRE].

[51] M.D. Johnson et al., The black hole explorer: motivation and vision, Proc. SPIE Int. Soc. Opt. Eng. 13092 (2024) 130922D [arXiv:2406.12917] [INSPIRE].

[52] S.E. Gralla, A. Lupsasca and A. Strominger, Observational signature of high spin at the Event Horizon Telescope, Mon. Not. Roy. Astron. Soc. 475 (2018) 3829 [arXiv:1710.11112] [INSPIRE].

[53] E. Himwich, M.D. Johnson, A. Lupsasca and A. Strominger, Universal polarimetric signatures of the black hole photon ring, Phys. Rev. D 101 (2020) 084020 [arXiv:2001.08750] [INSPIRE].

[54] S.E. Gralla, Measuring the shape of a black hole photon ring, Phys. Rev. D 102 (2020) 044017 [arXiv:2005.03856] [INSPIRE].

[55] S.E. Gralla, A. Lupsasca and D.P. Marrone, The shape of the black hole photon ring: a precise test of strong-field general relativity, Phys. Rev. D 102 (2020) 124004 [arXiv:2008.03879] [INSPIRE].

[56] S. Hadar, D. Kapec, A. Lupsasca and A. Strominger, Holography of the photon ring, Class. Quant. Grav. 39 (2022) 215001 [arXiv:2205.05064] [INSPIRE].

[57] A. Cárdenas-Avendaño and A. Lupsasca, Prediction for the interferometric shape of the first black hole photon ring, Phys. Rev. D 108 (2023) 064043 [arXiv:2305.12956] [INSPIRE].

[58] H. Jia, E. Quataert, A. Lupsasca and G.N. Wong, Photon ring interferometric signatures beyond the universal regime, Phys. Rev. D 110 (2024) 083044 [arXiv:2405.08804] [INSPIRE].

[59] S.E. Gralla and A. Lupsasca, Observable shape of black hole photon rings, Phys. Rev. D 102 (2020) 124003 [arXiv:2007.10336] [INSPIRE].

[60] F.H. Vincent, S.E. Gralla, A. Lupsasca and M. Wielgus, Images and photon ring signatures of thick disks around black holes, Astron. Astrophys. 667 (2022) A170 [arXiv:2206.12066] [INSPIRE].

[61] K. Salehi, R. Kumar Walia, D.O. Chang and P. Kocherlakota, Influence of observer’s inclination and spacetime structure on photon ring observables, Phys. Rev. D 111 (2025) 104057 [arXiv:2411.15310] [INSPIRE].
[62] R. Kumar Walia, P. Kocherlakota, D.O. Chang and K. Salehi, *Spacetime measurements with the photon ring*, Phys. Rev. D **111** (2025) 104074 [arXiv:2411.15119] [INSPIRE].

[63] A. Lupsasca et al., *The black hole explorer: photon ring science, detection, and shape measurement*, Proc. SPIE Int. Soc. Opt. Eng. **13092** (2024) 130926Q [arXiv:2406.09498] [INSPIRE].

[64] T. Johannsen, *Regular black hole metric with three constants of motion*, Phys. Rev. D **88** (2013) 044002 [arXiv:1501.02809] [INSPIRE].

[65] T. Johannsen, *Photon rings around Kerr and Kerr-like black holes*, Astrophys. J. **777** (2013) 170 [arXiv:1501.02814] [INSPIRE].

[66] V. Cardoso, P. Pani and J. Rico, *On generic parametrizations of spinning black-hole geometries*, Phys. Rev. D **89** (2014) 064007 [arXiv:1401.0528] [INSPIRE].

[67] G.O. Papadopoulos and K.D. Kokkotas, *Preserving Kerr symmetries in deformed spacetimes*, Class. Quant. Grav. **35** (2018) 185014 [arXiv:1807.08594] [INSPIRE].

[68] R.A. Konoplya, Z. Stuchlík and A. Zhidenko, *Axisymmetric black holes allowing for separation of variables in the Klein-Gordon and Hamilton-Jacobi equations*, Phys. Rev. D **97** (2018) 084044 [arXiv:1801.07195] [INSPIRE].

[69] Z. Carson and K. Yagi, *Asymptotically flat, parameterized black hole metric preserving Kerr symmetries*, Phys. Rev. D **101** (2020) 084030 [arXiv:2002.01028] [INSPIRE].

[70] Z. Carson and K. Yagi, *Probing beyond-Kerr spacetimes with inspiral-ringdown corrections to gravitational waves*, Phys. Rev. D **101** (2020) 084050 [arXiv:2003.02374] [INSPIRE].

[71] K. Yagi, S. Lomuscio, T. Lowrey and Z. Carson, *Regularizing parametrized black hole spacetimes with Kerr symmetries*, Phys. Rev. D **109** (2024) 044017 [arXiv:2311.08659] [INSPIRE].

[72] K. Salehi, A.E. Broderick and B. Georgiev, *Photon rings and shadow size for general axisymmetric and stationary integrable spacetimes*, Astrophys. J. **966** (2024) 143 [arXiv:2311.01495] [INSPIRE].

[73] B. Carter, *Hamilton-Jacobi and Schrödinger separable solutions of Einstein's equations*, Commun. Math. Phys. **10** (1968) 280 [INSPIRE].

[74] M. Walker and R. Penrose, *On quadratic first integrals of the geodesic equations for type [22] spacetimes*, Commun. Math. Phys. **18** (1970) 265 [INSPIRE].

[75] P. Krtous, D. Kubiznak, D.N. Page and V.P. Frolov, *Killing-Yano tensors, rank-2 Killing tensors, and conserved quantities in higher dimensions*, JHEP **02** (2007) 004 [hep-th/0612029] [INSPIRE].

[76] D. Kubiznak and V.P. Frolov, *Hidden symmetry of higher dimensional Kerr-NUT-AdS spacetimes*, Class. Quant. Grav. **24** (2007) F1 [gr-qc/0610144] [INSPIRE].

[77] P. Krtous, V.P. Frolov and D. Kubiznak, *Hidden symmetries of higher dimensional black holes and uniqueness of the Kerr-NUT-(A)dS spacetime*, Phys. Rev. D **78** (2008) 064022 [arXiv:0804.4705] [INSPIRE].

[78] V.P. Frolov, P. Krtous and D. Kubiznak, *Black holes, hidden symmetries, and complete integrability*, Living Rev. Rel. **20** (2017) 6 [arXiv:1705.05482] [INSPIRE].

[79] R. Konoplya, L. Rezzolla and A. Zhidenko, *General parametrization of axisymmetric black holes in metric theories of gravity*, Phys. Rev. D **93** (2016) 064015 [arXiv:1602.02378] [INSPIRE].

[80] P.T. Chrusciel, J. Lopes Costa and M. Heusler, *Stationary black holes: uniqueness and beyond*, Living Rev. Rel. **15** (2012) 7 [arXiv:1205.6112] [INSPIRE].
[81] M. Walker and R. Penrose, On quadratic first integrals of the geodesic equations for type [22] spacetimes, Commun. Math. Phys. 18 (1970) 265 [INSPIRE].

[82] B. Carter, Killing tensor quantum numbers and conserved currents in curved space, Phys. Rev. D 16 (1977) 3395 [INSPIRE].

[83] S.E. Gralla and A. Lupsasca, Null geodesics of the Kerr exterior, Phys. Rev. D 101 (2020) 044032 [arXiv:1910.12881] [INSPIRE].

[84] G. Compère, Y. Liu and J. Long, Classification of radial Kerr geodesic motion, Phys. Rev. D 105 (2022) 024075 [arXiv:2106.03141] [INSPIRE].

[85] G. Compère and A. Druart, Complete set of quasi-conserved quantities for spinning particles around Kerr, SciPost Phys. 12 (2022) 012 [arXiv:2105.12454] [INSPIRE].

[86] G. Compère, A. Druart and J. Vines, Generalized Carter constant for quadrupolar test bodies in Kerr spacetime, SciPost Phys. 15 (2023) 226 [arXiv:2302.14549] [INSPIRE].

[87] P. Ramond, On the integrability of extended test body dynamics around black holes, Class. Quant. Grav. 42 (2025) 065019 [arXiv:2402.02670] [INSPIRE].

[88] P. Ramond, Symplectic mechanics of relativistic spinning compact bodies I: covariant foundations and integrability around black holes, arXiv:2210.03866 [INSPIRE].

[89] V. Witzany, Hamilton-Jacobi equation for spinning particles near black holes, Phys. Rev. D 100 (2019) 104030 [arXiv:1903.03651] [INSPIRE].

[90] V. Skoupý and V. Witzany, Analytic solution for the motion of spinning particles in Kerr spacetime, Phys. Rev. Lett. 134 (2025) 171401 [arXiv:2411.16855] [INSPIRE].

[91] S.A. Teukolsky, Rotating black holes — separable wave equations for gravitational and electromagnetic perturbations, Phys. Rev. Lett. 29 (1972) 1114 [INSPIRE].

[92] F.J. Ernst, Removal of the nodal singularity of the C-metric, J. Math. Phys. 17 (1976) 515.

[93] M. Astorino, Removal of conical singularities from rotating C-metrics and dual CFT entropy, JHEP 10 (2022) 074 [arXiv:2207.14305] [INSPIRE].

[94] J.F. Plebanski and M. Demianski, Rotating, charged, and uniformly accelerating mass in general relativity, Annals Phys. 98 (1976) 98 [INSPIRE].

[95] T. Johannsen and D. Psaltis, Testing the no-hair theorem with observations in the electromagnetic spectrum: I. Properties of a quasi-Kerr spacetime, Astrophys. J. 716 (2010) 187 [arXiv:1003.3415] [INSPIRE].

[96] C. Batista and E.B. Santos, Conserved quantities for the free motion of particles with spin, Phys. Rev. D 101 (2020) 104049 [arXiv:2004.02966] [INSPIRE].

[97] D.C. Wilkins, Bound geodesics in the Kerr metric, Phys. Rev. D 5 (1972) 814 [INSPIRE].

[98] W. Schmidt, Celestial mechanics in Kerr space-time, Class. Quant. Grav. 19 (2002) 2743 [gr-qc/0202090] [INSPIRE].

[99] S. Drasco and S.A. Hughes, Rotating black hole orbit functionals in the frequency domain, Phys. Rev. D 69 (2004) 044015 [astro-ph/0308479] [INSPIRE].

[100] R. Fujita and W. Hikida, Analytical solutions of bound timelike geodesic orbits in Kerr spacetime, Class. Quant. Grav. 26 (2009) 135002 [arXiv:0906.1420] [INSPIRE].

[101] R.W. O'Shaughnessy, Transition from inspiral to plunge for eccentric equatorial Kerr orbits, Phys. Rev. D 67 (2003) 044004 [gr-qc/0211023] [INSPIRE].

[102] Y. Mino, Perturbative approach to an orbital evolution around a supermassive black hole, Phys. Rev. D 67 (2003) 084027 [gr-qc/0302075] [INSPIRE].
[103] J. Dexter and E. Agol, A fast new public code for computing photon orbits in a Kerr spacetime, Astrophys. J. 696 (2009) 1616 [arXiv:0903.0620] [INSPIRE].

[104] A. Druart, The motion of test bodies around Kerr black holes, Ph.D. thesis, U. Brussels, Brussels, Belgium (2023) [arXiv:2307.02589] [INSPIRE].

[105] S. Chandrasekhar, The mathematical theory of black holes, Oxford University Press, Oxford, U.K. (1998) [DOI:10.1093/oso/9780198503705.001.0001].

[106] B. O'Neill, The geometry of Kerr black holes, Dover, Mineloa, NY, U.S.A. (2014) [INSPIRE].

[107] E. Teo, Spherical photon orbits around a Kerr black hole, Gen. Rel. Grav. 35 (2003) 1909 [INSPIRE].

[108] É. Gourgoulhon, Geometry and physics of black holes, https://luth.obspm.fr/~luthier/gourgoulhon/bh16/bholes.pdf (2025).

[109] J.W. Moffat, Black holes in modified gravity (MOG), Eur. Phys. J. C 75 (2015) 175 [arXiv:1412.5424] [INSPIRE].

[110] A. Simpson and M. Visser, The eye of the storm: a regular Kerr black hole, JCAP 03 (2022) 011 [arXiv:2111.12329] [INSPIRE].

[111] J.W. Moffat, Scalar-tensor-vector gravity theory, JCAP 03 (2006) 004 [gr-qc/0506021] [INSPIRE].

[112] W. Liu et al., Kerr-MOG-(A)dS black hole and its shadow in scalar-tensor-vector gravity theory, JCAP 08 (2024) 035 [arXiv:2406.00579] [INSPIRE].

[113] R. Carballo-Rubio, F. Di Filippo, S. Liberati and M. Visser, Phenomenological aspects of black holes beyond general relativity, Phys. Rev. D 98 (2018) 124009 [arXiv:1809.08238] [INSPIRE].

[114] C. Lan, H. Yang, Y. Guo and Y.-G. Miao, Regular black holes: a short topic review, Int. J. Theor. Phys. 62 (2023) 202 [arXiv:2303.11696] [INSPIRE].

[115] R. Torres, Regular rotating black holes: a review, arXiv:2208.12713 [INSPIRE].

[116] A. Simpson and M. Visser, Regular black holes with asymptotically Minkowski cores, Universe 6 (2019) 8 [arXiv:1911.01020] [INSPIRE].

[117] C.A. Benavides-Gallego, S. Shashank and H. Xu, Observing the eye of the storm I: testing regular black holes with LVK and EHT observations, arXiv:2411.13897 [INSPIRE].

[118] M. Grould, F.H. Vincent, T. Paumard and G. Perrin, General relativistic effects on the orbit of the S2 star with GRAVITY, Astron. Astrophys. 608 (2017) A60 [arXiv:1709.04492] [INSPIRE].

[119] C. Bambi et al., Towards precision measurements of accreting black holes using X-ray reflection spectroscopy, Space Sci. Rev. 217 (2021) 65 [arXiv:2011.04792] [INSPIRE].

[120] M. Minamitsuji and K.-I. Maeda, Black hole thermodynamics in Horndeski theories, Phys. Rev. D 108 (2023) 084061 [arXiv:2308.01082] [INSPIRE].

[121] J. Ben Achour, Dost theories as disformal gravity: from black holes to radiative spacetimes, Eur. Phys. J. C 85 (2025) 424 [arXiv:2412.04135] [INSPIRE].