Kerr-MOG-(A)dS black hole and its shadow in scalar-tensor-vector gravity theory

Abstract

The scalar-tensor-vector gravity (STVG) theory has attracted significant interest due to its ability to effectively address the issue of galaxy rotation curves and clusters of galaxies without considering the influence of dark matter. In this paper, we construct rotating black hole solutions with a cosmological constant in the STVG theory (i.e., Kerr-MOG-(A)dS black hole solutions), where the import of a gravitational charge as a source modifies the gravitational constant, determined by $G=G_{\text{N}}(1+\alpha)$. For Kerr-MOG-dS spacetime, the observer is situated at a specific location within the domain of outer communication, rather than being located infinitely far away. Since black hole shadows are shaped by light propagation in spacetime, the interaction between the MOG parameter and the cosmological constant is expected to produce novel effects on these shadows. As the cosmological constant $\Lambda$ increases, the apparent size of the black hole shadow on the screen decreases. Additionally, the shadow expands with an increase in the MOG parameter $\alpha$, reaching a maximum at a certain value, and its shape becomes more rounded under an arbitrary rotation parameter, which leads to degeneracy between different black hole parameters. However, by employing numerical backward ray-tracing techniques, we have found that gravitational lensing and the frame-dragging effect effectively distinguish this degeneracy. Our work contributes to a deeper understanding of black holes in modified gravity, their observational signatures, and constraints.

I Introduction

Although the no-hair theorem suggests that black holes are characterized by their mass, spin, and electric charge, however, the electric charge is usually negligible for astrophysical black holes. The reason is that a black hole with mass $M$ and electric charge $Q$ will not gravitationally attract particles with mass $m$ and electric charge $e$ if $eQ>Mm$. Given the $m/e$ ratio for electrons is approximately $10^{-21}$, it is highly unlikely for black holes to maintain a significant electric charge Gibbons (1975). Moreover, processes including vacuum polarization, pair production, and neutralization ensure that stellar-mass black holes cannot sustain significant electric charge. Any accumulated charge dissipates quickly and escapes observation by gravitational wave detectors Gibbons (1975); Blandford and Znajek (1977); Cardoso et al. (2016). However, by performing general-relativistic simulations of charged black holes, Bozzola et al. have proposed that the gravitational wave event GW150914 could have a charge-to-mass ratio as high as $0.3$ Bozzola and Paschalidis (2021); Gupta et al. (2021); Carullo et al. (2022). Consequently, the term “charge” could be more appropriately interpreted as dark charge or gravitational charge, sparking renewed interest in the scalar-tensor-vector gravity (STVG) theory, also known as modified gravity (MOG) theory, which was first proposed nearly two decades ago Moffat (2006). The STVG theory introduces additional massive vector and scalar fields to the metric tensor field, positing that the associated “charge” is a gravitational charge denoted by $Q_{5}\equiv\kappa M$, which is determined by a body’s mass $M$ and its interaction with the vector field.

On the observational side, despite numerous attempts to detect dark matter, all experimental efforts up to now have yet to achieve recognized success An et al. (2023); Ning et al. (2023). Therefore, on the theoretical side, STVG theory has attracted a great deal of attention because it not only successfully explains the rotation curves of galaxies and the properties of galaxy clusters, but also does so without considering dark matter Moffat and Rahvar (2013, 2014); Moffat and Toth (2015); Moffat (2014); Rahvar (2022). In addition, observational data from the Event Horizon Telescope (EHT) on the shadows of Sgr A* and M87* Akiyama et al. (2019a, b, 2022a) support the predictions of STVG theory and there are criteria to distinguish this theory from general relativity Moffat and Toth (2020). The axially symmetric solution in STVG theory provided by Moffat et al. is called the Kerr-MOG black hole Moffat (2015a). To date, several new black hole solutions related to this theory have been obtained Moffat (2021); Moffat and Toth (2021); Liu et al. (2023a). Building on the established solutions in STVG theory, extensive research has been conducted on various aspects of black holes, including their thermodynamic properties Mureika et al. (2016), geodesic motion Lee and Han (2017), shadows Moffat (2015b); Wang et al. (2019); Wei and Liu (2018); Guo et al. (2018); Sheoran et al. (2018); Zhang et al. (2024); Qiao et al. (2020); Konoplya and Zhidenko (2021); Younsi et al. (2016), and quasi-normal modes Manfredi et al. (2018); Jiang et al. (2024), etc.

Despite much progress in the last few years in constructing black hole solutions in STVG theory and investigating their physical properties, rotating (A)dS black holes in STVG theory remain to be the virgin territory and thus need to be explored deeply. On the other hand, our current universe is in a state of accelerated expansion, and the simplest explanation for such accelerated expansion is a positive cosmological constant Riess et al. (1998); Perlmutter et al. (1999). Thus, it is significant to construct rotating dS black hole solutions in STVG theory. Additionally, the study of rotating AdS black holes has already shed light on the nature of gravity through gauge-gravity dualities Maldacena (1998); Gubser et al. (1998); Witten (1998), so it is very important and remarkable to seek the rotating AdS black hole solution in STVG theory. These three aspects motivate us to conduct the present work.

In this work, we aim to extend the Kerr-MOG solution Moffat (2015a) to more general cases with a cosmological constant, thereby obtaining Kerr-MOG-(A)dS black holes. Given the recent groundbreaking observations by the EHT, theoretical research on the properties of black holes frequently focuses on the analysis of black hole images Chen et al. (2023). Hence, we have also investigated the effects of the MOG parameter and a positive cosmological constant on black hole shadows, as well as the numerical estimation of the angular radius of the supermassive black holes Sgr A* and M87* within this theoretical framework. The remaining part of this paper is organized as follows. In Sec. II, we provide a concise overview of the STVG theory and give the Kerr-like solution of the field equations with a cosmological constant, i.e., Kerr-MOG-(A)dS black holes. In Sec. III, utilizing the geodesic equation, we derive the orbital equations for photons in the context of Kerr-MOG-dS spacetime. In Sec. IV, we examine the shadows of black holes, with a primary focus on the apparent shapes and distortions of the shadow. Also, we numerically calculate these two observables, $R_{s}$ and $\delta_{s}$, to characterize the shadow. Sec. V focuses on providing an overview of our conclusions while also contemplating the possible development and exploration of related research topics in the future.

II STVG theory and black hole solution

In this section, we first briefly review the STVG theory and then proceed to construct new rotating black hole solutions using the Kerr-Schild-(A)dS metric ansatz as a seed solution.

II.1 STVG theory

The action corresponding to STVG theory is given by Moffat (2006)

$$\mathcal{S}=\mathcal{S}_{\text{M}}+\mathcal{S}_{\text{G}}+\mathcal{S}_{\phi}+\mathcal{S}_{\text{S}},$$ (1)

with

$$\mathcal{S}_{\text{G}}= \frac{1}{16\pi}\int d^{4}x\sqrt{-g}\left[\frac{1}{G}\left(R-2\Lambda\right)\right],$$ (2)

$$\mathcal{S}_{\phi}= \int d^{4}x\sqrt{-g}\left(-\frac{1}{4}B^{ab}B_{ab}+\frac{1}{2}\mu^{2}\phi^{a}\phi_{a}\right),$$

$$\mathcal{S}_{\text{S}}= \int d^{4}x\sqrt{-g}\left[\frac{1}{G^{3}}\left(\frac{1}{2}g^{ab}\nabla_{a}G\nabla_{b}G-V(G)\right)\right.$$

$$\left.+\frac{1}{\mu^{2}G}\left(\frac{1}{2}g^{ab}\nabla_{a}\mu\nabla_{b}\mu-V(\mu)\right)\right],$$

where the symbol $\mathcal{S}_{\text{M}}$ represents the matter action, $\mathcal{S}_{\text{G}}$ represents the action for Einstein gravity, and $R$ and $\Lambda$ are the Ricci scalar and cosmological constant, respectively. Additionally, the action $\mathcal{S}_{\phi}$ describes a Proca field $\phi^{a}$ with mass $\mu$. The field strength, represented by $B_{ab}=\partial_{a}\phi_{b}-\partial_{b}\phi_{a}$, satisfies the following equations:

$$\nabla_{b}B^{ab}=\frac{1}{\sqrt{-g}}\partial_{b}\left(\sqrt{-g}B^{ab}\right)=0,$$ (3)

$$\partial_{c}B_{ab}+\partial_{a}B_{bc}+\partial_{b}B_{ca}=0.$$ (4)

The potentials $V(G)$ and $V(\mu)$ are linked to the scalar fields $G(x)$ and $\mu(x)$ within the action $\mathcal{S}_{\mathrm{S}}$, respectively.

Under the assumption of vacuum solution, the action $\mathcal{S}$ can be simplified to:

$$\mathcal{S}= \int d^{4}x\sqrt{-g}\Big{[}\frac{1}{16\pi G}\left(R-2\Lambda\right)$$ (5)

$$-\frac{1}{4}B^{ab}B_{ab}+\frac{1}{2}\mu^{2}\phi^{a}\phi_{a}\Big{]}.$$

The field equation of motion derived from varying the action in Eq. (5) with respect to the metric is given by

$$\mathcal{G}_{ab}+\Lambda g_{ab}=-8\pi GT_{ab},$$ (6)

where $\mathcal{G}_{ab}=R_{ab}-\frac{1}{2}Rg_{ab}$ is Einstein tensor, and the energy-momentum tensor is

$$\begin{aligned} T_{ab}=&-\frac{1}{4\pi}\Big{(}B_{a}\hskip 0.28453pt^{c}B_{bc}-\frac{1}{4}g_{ab}B^{cd}B_{cd}\Big{)}\\ &+\frac{\mu^{2}}{4\pi}\Big{(}\phi_{a}\phi_{b}-\frac{1}{2}g_{ab}\phi^{c}\phi_{c}\Big{)}.\end{aligned}$$ (7)

The first and second terms of the energy-momentum tensor are of the order $\sim(\partial\phi)^{2}$ and $\sim\mu^{2}\phi^{2}$, respectively. To solve the equation (6), one can typically use the weak field approximation and neglect the higher-order terms including $g^{2}$, $g\phi$, and $\phi^{2}$ Moffat and Rahvar (2013); Rahvar (2022). Concurrently, according to Refs. Moffat and Rahvar (2014), the particle mass of the $\phi^{a}$ field in the present universe can be fitted as $m_{\phi}\sim 10^{-28}eV$, making it negligible for a black hole solution. In STVG theory, it is important to note that the fifth force charge $Q_{5}$ is tied to a particle’s inertial mass as $Q_{5}=\kappa M$. This parameter $\kappa$ affects the gravitational constant, given by $G=G_{\text{N}}+\kappa^{2}$, with the convention $\kappa^{2}=\alpha G_{\text{N}}$. Then, we can obtain the relationship as $G=(1+\alpha)G_{\text{N}}$, which can also be thought of as a constant unaffected by the spacetime coordinates. By setting $\alpha$ to zero, the theory reduces to General Relativity (GR). Thus, $\alpha$ can be considered as a parameter that measures the deviation of STVG from GR.

In the static case, assuming the vector field $\phi_{a}$ takes the form Moffat (2015a)

$$\phi_{a}=\left(-\sqrt{\alpha G_{\text{N}}}\frac{M}{r},0,0,0\right),$$ (8)

and setting $G_{\text{N}}=1$ and $\beta=\alpha/(1+\alpha)$, the corresponding asymptotic (A)dS metric solution is given by Liu et al. (2023a)

$$ds^{2}= -\left(1-\frac{2M_{\text{D}}}{r}+\frac{\beta M^{2}_{\text{D}}}{r^{2}}-\frac{\Lambda}{3}r^{2}\right)dt^{2}$$ (9)

$$+\left(1-\frac{2M_{\text{D}}}{r}+\frac{\beta M^{2}_{\text{D}}}{r^{2}}-\frac{\Lambda}{3}r^{2}\right)^{-1}dr^{2}$$

$$+r^{2}d\vartheta^{2}+r^{2}\sin^{2}\vartheta d\varphi^{2}.$$

The symbol $M_{\text{D}}$ stands for the ADM mass, which can be related to the Newtonian mass $M$ through the equation $M_{\text{D}}=(1+\alpha)M$ Sheoran et al. (2018).

II.2 Rotating-(A)dS solutions

Rotating black hole solutions play a crucial role in astrophysics. Considering the total spin angular momentum represented by $J=Ma$, adjustments to the vector potential and metric, equations (8) to (9), are necessary to incorporate the spin parameter. In order to obtain the axisymmetric solutions for asymptotically (A)dS spacetimes, we adopt the Kerr-Schild-(A)dS ansatz Gibbons et al. (2005); Malek and Pravda (2011)

$$g_{ab}=\bar{g}_{ab}-2\mathcal{H}k_{a}k_{b},$$ (10)

where $\mathcal{H}$ is an arbitrary scalar function. The vector $\mathbf{k}$ is null with respect to both the full metric $g_{ab}$ and the Kerr-(A)dS metric $\bar{g}_{ab}$, with the latter serving as the background. Alternatively, by introducing the coordinate $(\bar{t},\bar{r},\bar{\vartheta},\bar{\varphi})$, the line elements can be expressed as follows:

$$d\bar{s}^{2}= -\left(1-\frac{\Lambda}{3}\bar{r}^{2}\right)\frac{\Delta_{\bar{\vartheta}}}{\chi}d\bar{t}^{2}+\frac{\Sigma}{\Delta_{\bar{r}}}d\bar{r}^{2}$$ (11)

$$+\frac{\Sigma}{\Delta_{\bar{\vartheta}}}d{\bar{\vartheta}}^{2}+\frac{(\bar{r}^{2}+a^{2})\sin^{2}{\bar{\vartheta}}}{\chi}d{\bar{\varphi}}^{2}$$

$$-2\mathcal{H}\left(\frac{\Delta_{\bar{\vartheta}}}{\chi}d\bar{t}+\frac{\Sigma}{\Delta_{\bar{r}}}d\bar{r}-\frac{a\sin^{2}{\bar{\vartheta}}}{\chi}d{{\bar{\varphi}}}\right)^{2},$$

with

$$\Sigma= \bar{r}^{2}+a^{2}\cos^{2}{\bar{\vartheta}}, \Delta_{\bar{\vartheta}}= 1+\frac{\Lambda}{3}a^{2}\cos^{2}{\bar{\vartheta}},$$ (12)

$$\chi= 1+\frac{\Lambda}{3}a^{2}, \Delta_{\bar{r}}= (\bar{r}^{2}+a^{2})\left(1-\frac{\Lambda}{3}\bar{r}^{2}\right).$$

For the metric ansatz (11), the Ricci scalar is

$$R=g^{ab}R_{ab}=4\Lambda-2\left(\frac{\partial^{2}}{\partial\bar{r}^{2}}+\frac{4\bar{r}}{\Sigma}\frac{\partial}{\partial\bar{r}}+\frac{2}{\Sigma}\right)\mathcal{H},$$ (13)

and the vacuum gravitational field equations (with cosmological constant) are $R=4\Lambda$. Therefore, one can obtain the following partial differential equation for the $\mathcal{H}$ function:

$$\left(\frac{\partial^{2}}{\partial\bar{r}^{2}}+\frac{4\bar{r}}{\Sigma}\frac{\partial}{\partial\bar{r}}+\frac{2}{\Sigma}\right)\mathcal{H}=0.$$ (14)

The above equation has the general solution,

$$\mathcal{H}=\frac{GQ_{5}^{2}}{2\Sigma}-\frac{GM\bar{r}}{\Sigma},$$ (15)

where $GM$ and $GQ_{5}^{2}$ are constants of integration, as in Ref. Moffat (2015a). In STVG theory, the charge $Q_{5}>0$ represents a fifth force with a gravitational origin and is proportional to the mass of the source particle, given by $Q_{5}^{2}=\alpha G_{\text{N}}M^{2}$. Together with the relation $G=G_{\text{N}}(1+\alpha)$, $\mathcal{H}$ can be rewritten as follows:

$$\mathcal{H}=\frac{G^{2}_{\text{N}}M^{2}\alpha\left(1+\alpha\right)}{2\Sigma}-\frac{G_{\text{N}}M\bar{r}(1+\alpha)}{\Sigma}.$$ (16)

Then, the seed metric (11) with scalar function (16) can be expressed as

$$ds^{2}= -\frac{\Delta_{\bar{\vartheta}}}{\chi}\left(1-\frac{\Lambda\bar{r}^{2}}{3}+\frac{h\Delta_{\bar{\vartheta}}}{\chi\Sigma}\right)d\bar{t}^{2}+\frac{\Sigma}{\Delta_{\bar{r}}}\left(1-\frac{h}{\Delta_{\bar{r}}}\right)d\bar{r}^{2}$$

$$+\frac{\Sigma}{\Delta_{\bar{\vartheta}}}d\bar{\vartheta}^{2}+\left(\frac{\bar{r}^{2}+a^{2}}{\chi}-\frac{ha^{2}\sin^{2}\bar{\vartheta}}{\chi^{2}\Sigma}\right)\sin^{2}\bar{\vartheta}d\bar{\varphi}^{2}$$

$$-\frac{2h\Delta_{\bar{\vartheta}}}{\chi\Delta_{\bar{r}}}d\bar{t}d\bar{r}+\left(\frac{\Delta_{\bar{\vartheta}}d\bar{t}}{\chi\Sigma}+\frac{d\bar{r}}{\Delta_{\bar{r}}}\right)\frac{2ha\sin^{2}\bar{\vartheta}}{\chi}d\bar{\varphi},$$ (17)

where

$$h=G_{\text{N}}M(1+\alpha)(G_{\text{N}}M\alpha-2\bar{r}).$$ (18)

The non-zero components of the effective tensor, defined as $E_{ab}=R_{ab}-\frac{1}{2}Rg_{ab}+\Lambda g_{ab}+8\pi GT_{ab}$, associated with the metric, are

$$E_{00}=G_{\alpha}\Delta_{\bar{\vartheta}}\left[(1-\bar{r}^{2}\Lambda/3)\Sigma_{2}+h\Delta_{\bar{\vartheta}}\right]+8\pi GT_{00},$$ (19)

$$E_{01}=G_{\alpha}h\chi\Delta_{\bar{r}}^{-1}\Delta_{\bar{{\vartheta}}}+8\pi GT_{01},$$

$$E_{03}=-G_{\alpha}\Delta_{\bar{\vartheta}}\left(h+2\Delta_{\bar{r}}\right)a\sin^{2}\bar{\vartheta}+8\pi GT_{03},$$

$$E_{11}=G_{\alpha}\left(h-\Delta_{\bar{r}}\right)\chi^{2}\Delta_{\bar{r}}^{-2}\Sigma^{2}+8\pi GT_{11},$$

$$E_{13}=-G_{\alpha}h\chi\Delta_{\bar{\vartheta}}^{-1}\Sigma a\sin^{2}\bar{\vartheta}+8\pi GT_{13},$$

$$E_{22}=G_{\alpha}\chi^{2}\Delta_{\bar{\vartheta}}^{-1}\Sigma^{2}+8\pi GT_{22},$$

$$E_{33}=G_{\alpha}\sin^{2}\bar{\vartheta}\left[(\bar{r}^{2}+a^{2})\Sigma_{2}+ha^{2}\sin^{2}\bar{\vartheta}\right]+8\pi GT_{33},$$

where

$$G_{\alpha} =G_{\text{N}}^{2}M^{2}\alpha(1+\alpha)\chi^{-2}\Sigma^{-3},$$ (20)

$$\Sigma_{2} =\Sigma+a^{2}\Delta_{\bar{\vartheta}}-a^{2}\left(1-r^{2}\Lambda/3\right)(\cos^{2}\bar{\vartheta}-\sin^{2}\bar{\vartheta}).$$

The physical significance of these parameters becomes clearer in the Boyer-Lindquist coordinates $(t,r,\vartheta,\varphi)$, which offer a more appropriate asymptotic form for the metric. Applying the Boyer-Lindquist transformation,

$$d\bar{r}= dr,~{}~{}~{}d{\bar{\vartheta}}=d\vartheta,$$ (21)

$$d\bar{t}= dt-\frac{h}{\left(1-\frac{\Lambda}{3}r^{2}\right)\left(\Delta_{\bar{r}}+h\right)}dr,$$

$$d\bar{\varphi}= d\varphi-\frac{\Lambda}{3}adt-\frac{ah}{\left(r^{2}+a^{2}\right)\left(\Delta_{\bar{r}}+h\right)}dr,$$

the metric (II.2) can be rewritten in the following form,

$$ds^{2}= -\frac{\Delta_{r}}{\Sigma}\left(dt-\frac{a\sin^{2}\vartheta}{\chi}d\varphi\right)^{2}+\frac{\Sigma}{\Delta_{r}}dr^{2}+\frac{\Sigma}{\Delta_{\vartheta}}d\vartheta^{2}$$ (22)

$$+\frac{\Delta_{\vartheta}\sin^{2}\vartheta}{\Sigma}\left(adt-\frac{(r^{2}+a^{2})}{\chi}d\varphi\right)^{2},$$

with

$$\Delta_{\vartheta}=\Delta_{\bar{\vartheta}},~{}~{}~{}~{}~{}\Delta_{r}=\Delta_{\bar{r}}-2G_{\text{N}}M_{\text{D}}r+\frac{\alpha G_{\text{N}}^{2}M^{2}_{\text{D}}}{(1+\alpha)}.$$

To ensure the effective tensor $E_{ab}=0$ is satisfied, we can simply obtain the MOG vector potential as:

$$\phi=\phi_{a}dx^{a}=\frac{M\sqrt{G_{\text{N}}\alpha}r}{\Sigma}\left(-dt+\frac{a\sin^{2}\vartheta}{\chi}d\varphi\right).$$ (23)

In Appendix A, we provide the specific forms of all the energy-momentum tensors contained in equations (LABEL:E00), which are sourced from the MOG vector potential (23). Equations (22) and (23) describe the Kerr-MOG-(A)dS black hole solutions. When $\Lambda=0$, the solutions recovers the usual Kerr-MOG solution Moffat (2015a); when $a=0$, it reduces to the static (A)dS cases Liu et al. (2023a). This metric (22) has the similar form as the Kerr-Newman-(A)dS metric. However, in this case, the ADM mass $M_{\text{D}}$ is associated with the coupling parameter $\alpha$, and the source of the vector field is gravitational. This distinction is crucial because astrophysical objects, including black holes, are electrically neutral.

The equations of the horizons are given by:

$$\Delta_{r}=\frac{\Lambda}{3}\left(r-r_{m}\right)\left(r-r_{h}\right)\left(r_{c}-r\right)\left(r+r_{b}\right)=0,$$ (24)

where $r_{m}$, $r_{h}$, and $r_{c}$ correspond to the Cauchy horizon, event horizon, and dS cosmological horizon, respectively Liu et al. (2023a). The parameters $r_{b}$, $r_{m}$, $r_{h}$, and $r_{c}$ satisfy the condition $r_{b}+r_{m}+r_{h}+r_{c}=0$. Moreover, $r_{c}\geq r_{h}>r_{m}>0$ and $r_{b}$ is negative.

III Photon orbits

In this section, we provide a brief derivation of photon trajectories in the gravitational field of Kerr-MOG-(A)dS black holes. When a black hole is located between the observer and the extended background source, not all photons emitted by the source can reach the observer after being deflected by the black hole gravitational field. The apparent shape of the black hole is represented by its shadow contour, which is closely related to the geodesics of photons within the black hole spacetime.

The geodesics corresponding to photons geometry are defined by the Hamilton-Jacobi equation:

$$\frac{\partial\mathcal{S}}{\partial\tau}=-\frac{1}{2}g^{ab}\frac{\partial\mathcal{S}}{\partial x^{a}}\frac{\partial\mathcal{S}}{\partial x^{b}},$$ (25)

where $\tau$ is an affine parameter along the null geodesics, and $\mathcal{S}$ is the Jacobi action. We can set up the Killing fields $\xi^{a}=(\partial/\partial t)^{a}$ and $\psi^{a}=(\partial/\partial\varphi)^{a}$, then the corresponding energy $\mathcal{E}$ and the angular momentum in the direction of the axis of symmetry $L$ can be defined as Wald (1984); Zeng et al. (2022); Liu et al. (2023b):

$$\mathcal{E}= -g_{ab}\xi^{a}\dot{x}^{b}=-g_{00}\dot{t}-g_{03}\dot{\varphi},$$ (26)

$$L= g_{ab}\psi^{a}\dot{x}^{b}=g_{03}\dot{t}+g_{33}\dot{\varphi}.$$ (27)

When $\mathcal{S}$ is separable, the Jacobi action for the photons can be expressed in a straightforward and general form:

$$\mathcal{S}=-\mathcal{E}t+L\varphi+\mathcal{S}_{r}(r)+\mathcal{S}_{\vartheta}(\vartheta).$$ (28)

Combining Eqs. (25) and (28) and setting $G_{\text{N}}=1$ , we obtain two equations of motion for the propagating photons

$$\dot{t}= \frac{\left(r^{2}+a^{2}\right)k}{\Sigma\Delta_{r}}+\frac{a\left(L\chi-a\mathcal{E}\sin^{2}\vartheta\right)}{\Sigma\Delta_{\vartheta}},$$ (29)

$$\dot{\varphi}= \frac{a\chi k}{\Sigma\Delta_{r}}+\frac{\chi\left(L\chi\csc^{2}\vartheta-a\mathcal{E}\right)}{\Sigma\Delta_{\vartheta}},$$ (30)

with

$$k=\left(r^{2}+a^{2}\right)\mathcal{E}-aL\chi.$$ (31)

Substituting Eqs. (29) and (30) into the null geodesics equation $0=-g_{ab}\dot{x}^{a}\dot{x}^{b}$, one can obtain

$$\frac{k^{2}-\rho^{4}\dot{r}^{2}}{\Delta_{r}}-\frac{\rho^{4}\dot{\vartheta}^{2}+\left(L\chi\csc\vartheta-a\mathcal{E}\sin\vartheta\right)^{2}}{\Delta_{\vartheta}}=0.$$ (32)

By using the Carter constant $\mathcal{K}$ to separate variables in the above equation, we can obtain two additional equations of motion for photons Carter (1968)

$$\dot{r}= \frac{1}{\Sigma}\sqrt{\mathcal{R}(r)},$$ (33)

$$\dot{\vartheta}= \frac{1}{\Sigma}\sqrt{\Theta(\vartheta)}.$$ (34)

The explicit forms of the radial and angular equations are as follows:

$$\mathcal{R}= k^{2}-\left[\mathcal{K}+\left(L\chi-a\mathcal{E}\right)^{2}\right]\Delta_{r},$$ (35)

$$\Theta= \left[\mathcal{K}+\left(L\chi-a\mathcal{E}\right)^{2}\right]\Delta_{\vartheta}-\left(a\mathcal{E}\sin\vartheta-L\chi\csc\vartheta\right)^{2}.$$ (36)

These equations determine the behavior of light propagation within the spacetime influenced by rotating black holes.

The radial motion equation of the photon is of particular interest, as the silhouette of the black hole shadow can be derived from the orbits with constant $r=r_{p}$ that fulfill the conditions

$$\mathcal{R}(r)\big{|}_{r=r_{p}}=0, \frac{d}{dr}\mathcal{R}(r)\big{|}_{r=r_{p}}=0.$$ (37)

To determine the contour of black hole shadow, one can introduce the impact parameters $\xi=L/\mathcal{E}$ and $\eta=\mathcal{K}/\mathcal{E}^{2}$, and solve the system of equations accordingly Chandrasekhar (1984). The resulting impact parameters will provide the necessary information:

$$\xi= \frac{M_{\mathrm{D}}\left(a^{2}+2M\alpha r_{p}-3r^{2}_{p}\right)+\chi r_{p}\left(a^{2}+r^{2}_{p}\right)}{a\chi\left[M-r_{p}+M\alpha+\frac{1}{3}r_{p}\Lambda\left(a^{2}+2r_{p}^{2}\right)\right]},$$ (38)

$$\eta= \frac{r_{p}^{2}\mathcal{W}-M_{\mathrm{D}}\Lambda r_{p}^{5}\left(2a^{2}+\frac{4}{3}M\alpha r_{p}\right)+r_{p}^{6}-r_{p}^{6}\chi^{2}}{a^{2}\left[M-r_{p}+M\alpha+\frac{1}{3}r_{p}\Lambda\left(a^{2}+2r_{p}^{2}\right)\right]^{2}},$$ (39)

where

$$\mathcal{W}= 4Mr_{p}(1+\alpha)\left[a^{2}+M^{2}\alpha(1+\alpha)\right]$$ (40)

$$-4M^{2}\alpha(1+\alpha)\Delta_{r}-r_{p}^{2}\left[r_{p}-3M(1+\alpha)\right]^{2}.$$

These two impact parameters are equivalent to the results provided in Ref. Wang et al. (2019) when $\Lambda=0$. For simplicity, we set $M=1$ in subsequent calculations.

IV SHADOWS OF BLACK HOLES

While AdS solutions play a crucial role in AdS/CFT correspondences, they lack empirical support from specific astronomical observations. Therefore, in this section, we primarily focus on the black hole shadows in Kerr-MOG-dS spacetime. For this, we need to adopt the analysis approach of A. Grenzebach Grenzebach et al. (2014) and the numerical backward ray-tracing method Wang et al. (2018a, 2021); Hu et al. (2021); Zhong et al. (2021); Chen and Jing (2024) to investigate the shadow of Kerr-MOG black holes in the presence of a positive cosmological constant.

IV.1 Apparent shape

In the Kerr-MOG-dS spacetime, there exist two unstable light orbits within the domain of outer communication, which extends from the outer horizon $r_{h}$ to the cosmological horizon $r_{c}$: a prograde orbit and a retrograde orbit. Both orbits serve as limit curves for lightlike geodesics that spiral towards them. Similarly to the circular lightlike geodesics in the Schwarzschild spacetime Perlick and Tsupko (2022), all spherical lightlike geodesics in the domain of outer communication are also unstable and can serve as limit curves for past-oriented light rays from an observer position that spiral towards them. Hence, the photon region determines the boundary curve of the black hole shadow. Assume that the observer is positioned at $(r_{0},\vartheta_{0})$, where $r_{h}<r_{0}<r_{c}$. We define an orthonormal tetrad as Eiroa and Sendra (2018):

$$e_{0}= \left.\frac{r^{2}+a^{2}}{\sqrt{\Sigma\Delta_{r}}}\partial_{t}+\frac{a\chi}{\sqrt{\Sigma\Delta_{r}}}\partial_{\varphi}\right|_{(r_{0},\vartheta_{0})},$$ (41)

$$e_{1}= \left.\sqrt{\frac{\Delta_{\vartheta}}{\Sigma}}\partial_{\vartheta}\right|_{(r_{0},\vartheta_{0})},$$

$$e_{2}= \left.-\frac{\chi\csc\vartheta}{\sqrt{\Sigma\Delta_{\vartheta}}}\partial_{\varphi}-\frac{a\sin\vartheta}{\sqrt{\Sigma\Delta_{\vartheta}}}\partial_{t}\right|_{(r_{0},\vartheta_{0})},$$

$$e_{3}= \left.-\sqrt{\frac{\Delta_{r}}{\Sigma}}\partial_{r}\right|_{(r_{0},\vartheta_{0})}.$$

It is easy to verify that $e_{0}$, $e_{1}$, $e_{2}$, and $e_{3}$ are orthonormal. In this chosen tetrad, $e_{0}$ represents the observer four-velocity, $e_{3}$ corresponds to the spatial direction towards the center of the black hole, and the combination $e_{0}\pm e_{3}$ is tangential to the principal null congruences of the metric. For each light ray $\lambda(s)$ with coordinate representation $t(s),r(s),\vartheta(s),\varphi(s)$, the general form of tangent vector is

$$\dot{\lambda}=\dot{t}\partial_{t}+\dot{r}\partial_{r}+\dot{\vartheta}\partial_{\vartheta}+\dot{\varphi}\partial_{\varphi}.$$ (42)

Following the approach in Ref. Grenzebach et al. (2014), the tangent vector at the observer can be expressed as

$$\dot{\lambda}=\gamma(\sin\theta\cos\psi e_{1}+\sin\theta\sin\psi e_{2}+\cos\theta e_{3}-e_{0}),$$ (43)

where $\gamma$ is a scalar factor and can be solved by substituting Eqs. (29), (30), (33), and (34), and Eq. (43) into Eq. (42). The result gives

$$\gamma=-\frac{k}{\sqrt{\Sigma\Delta_{r}}},$$ (44)

and then yield

$$\sin\theta= \left.\frac{\sqrt{\eta\Delta_{r}+\left(a-\chi\xi\right)^{2}\Delta_{r}}}{r^{2}+a^{2}-a\chi\xi}\right|_{r=r_{0}},$$ (45)

$$\sin\psi= \left.\frac{\csc\vartheta\left(a\cos^{2}\vartheta-a+\chi\xi\right)}{\sqrt{\eta\Delta_{\vartheta}+\left(a-\chi\xi\right)^{2}\Delta_{\vartheta}}}\right|_{\vartheta=\vartheta_{0}}.$$ (46)

In order to use analytical parameter to roperesent the boundary curve of the shadow Grenzebach et al. (2014), we utilized the stereographic projection from the celestial sphere onto a plane, where standard Cartesian coordinates are employed,

$$x(r_{p})= -2\tan\left[\frac{\theta(r_{p})}{2}\right]\sin\left[\psi(r_{p})\right],$$ (47)

$$y(r_{p})= -2\tan\left[\frac{\theta(r_{p})}{2}\right]\cos\left[\psi(r_{p})\right].$$ (48)

This is a set of parameter equations concerning the radius of the photon sphere $r_{p}$ Perlick and Tsupko (2022).

For the non-rotating black hole case with $a=0$, the expression for the impact parameter has a singularity, necessitating separate consideration of the static case. Using the circular orbit condition (37) and setting $\eta+\xi^{2}\chi^{2}:=1/X$ Grenzebach et al. (2014), the Eqs. (45)-(46) have the form:

$$\sin\theta^{(a=0)}= \left.\frac{\sqrt{\Delta_{r}|_{a=0}}}{r^{2}\sqrt{X}}\right|_{r=r_{0}},$$ (49)

$$\sin\psi^{(a=0)}= \left.\sqrt{X}\chi\xi\csc\vartheta\right|_{\vartheta=\vartheta_{0}},$$ (50)

where

$$X=\frac{1}{2r_{p}^{2}}-\frac{(1+\alpha)}{2r_{p}^{3}}-\frac{\Lambda}{3},$$ (51)

and the radius of the photon sphere is

$$r_{p}=\frac{3}{2}(1+\alpha)\left(1+\sqrt{1-\frac{8\alpha}{9(1+\alpha)}}\right).$$ (52)

In this case, the system of Eqs. (47)-(48) will become parameter equations only concerning the impact parameter $\xi$. Note that the expression of the photon sphere radius is independent of $\Lambda$. From this, one can deduce that the shape of a static black hole is a standard circle, with its radius dependent on the MOG parameter $\alpha$.

Figure 1: Shadows cast by black holes with a cosmological constant, as seen by an observer at $r_{0}=50$ and $\vartheta_{0}=\pi/2$. The (a) and (c) are static cases $a/M=0$, corresponding to $\Lambda=10^{-4}$ and $\Lambda=10^{-3}$, respectively. The (b) and (d) are rotating case $\alpha=0.02$, corresponding to $\Lambda=10^{-4}$ and $\Lambda=10^{-3}$, respectively.

Figure 2: Shadows cast by black holes with a spin $a/M=0.99$, as seen by an observer at $r_{0}=50$ and $\vartheta_{0}=\pi/2$. The (a), (b), (c) and (d) corresponding to $\alpha=0$, $\alpha=1$, $\alpha=2$ and $\alpha=3$, respectively.

Figure 3: Shadows cast by black holes with a cosmological constant, as seen by an observer at $r_{0}=50$ and $\vartheta_{0}=\pi/2$. The (a) and (c) are spin $a/M=0.5$, corresponding to $\Lambda=10^{-4}$ and $\Lambda=10^{-3}$, respectively. The (b) and (d) are spin $a/M=0.99$, corresponding to $\Lambda=10^{-4}$ and $\Lambda=10^{-3}$, respectively.

For simplicity, we consider fixing the position of the observer at $r_{0}=50$ and $\vartheta_{0}=\pi/2$. Where the angle $\vartheta_{0}$ is defined as the inclination angle between the observer’s line of sight and the rotation axis of the black hole. Form Eqs. (47) and (48) one can depict the contour of the black hole shadow under varying MOG parameter $\alpha$, spin $a$, or cosmological constant $\Lambda$. In Fig. 1, the shadows of static black holes and rotating black holes are shown under $\Lambda=10^{-4}$ or $\Lambda=10^{-3}$. When the cosmological constant is relatively small, the shadow scale exhibits an increase with the growth of the MOG parameter. However, as the cosmological constant increases, this trend changes. This shift can be attributed to the observer’s position $r_{0}$ approaching the cosmological horizon $r_{c}$, causing the observer’s position to have a more significant effect on the shadow size than the effects of the MOG modification. In Fig. 2, we present the shadow contour of a near-extreme black hole with different $\alpha$ and $\Lambda$, i.e., $a/M=0.99$. The results show that, as the MOG parameter increases, the range of the photon sphere radius narrows, and its contour increasingly resembles a standard circle. When the cosmological constant is relatively small, the shadow size increases with the growth of the MOG parameter. For a fixed $a$ and $\alpha$, the apparent shape of the shadow varies with $\Lambda$. As $\Lambda$ decreases, the shadow’s size expands. In Fig. 3, the same phenomenon that is clearly visible in Fig. 1-(c) can also be observed in rotating black holes. Furthermore, when the MOG parameter increases significantly, the rotation-induced distortion of the photon orbits becomes less noticeable. It is worth noting that in Figs. 1-3, the purple solid line represents the shadow of the Schwarzschild black hole, serving as a scale reference.

The contours of static versus rotating black holes will be extremely similar when the MOG parameter is sufficiently large Camilloni et al. (2024), making it difficult to infer black hole parameters solely from the black hole shadow. However, does the frame-dragging effect caused by rotation in black hole spacetimes also affect the MOG parameter? For this, we can use the numerical backward ray-tracing method, which adopts the same orthonormal tetrad as the analytical methods in this paper and sets the number of pixels to $n=2048$. Figs. 4 and 5 effectively demonstrate the distortion of space by MOG black holes and the gravitational lensing effects these objects produce. Specifically, the gravitational lensing observed can be utilized to examine theories of gravity Liu et al. (2024); Chowdhuri et al. (2023); Ghosh and Bhattacharyya (2022). We can observe that although increasing the MOG parameter makes the contours of rotating black holes approach those of static black holes, it does not affect the frame-dragging effect induced by rotation. Moreover, Figs. 4 and 5 serve as double checks for Figs. 1-(a) and 2-(b), respectively, with each displaying the same black hole contours.

Figure 4: Shadows cast by static black holes ($a/M=0$) with a cosmological constant, as seen by an observer at $r_{0}=50$ and $\vartheta_{0}=\pi/2$.

Figure 5: Shadows cast by extremely rotating black holes ($a/M=0.99$) with a cosmological constant, as seen by an observer at $r_{0}=50$ and $\vartheta_{0}=\pi/2$.

IV.2 Observable measurements

Up to this point, we have determined that parameters such as spin $a$, MOG parameter $\alpha$, and cosmological constant $\Lambda$ will significantly affect the apparent shape of the shadow in Kerr-MOG-dS spacetime. To characterize the shadow, it is essential to have easily measurable and reliable properties of observable measurements. We adopt two such observable measurements: the on-screen radius $R_{s}$ and distortion parameter $\delta_{s}$, following an approach analogous to the one taken in Ref. Hioki and Maeda (2009).

Figure 6: Shadows cast by black holes, as seen by an observer at $r_{0}=50$ and $\vartheta_{0}=\pi/2$. The observable measurements for the apparent shape of the black hole are the radius of the reference circle $R_{s}$ and the distortion parameter $\delta_{s}=\mathrm{D}/R_{s}$, where $\mathrm{D}$ is the difference between the left endpoints of the reference circle and of the shadow.

The parameter $R_{s}$ characterizes the on-screen size of the shadow, corresponding to the radius of the reference circle depicted as a blue dashed line in Fig. 6. The circle passes through three points: the top position $P_{A}=(x_{t},y_{t})$, the bottom position $P_{B}=(x_{b},y_{b})$, and the point $P_{C}=(x_{r},0)$, which corresponds to the unstable retrograde circular orbit as observed from the equatorial plane by an observer. Additionally, the point $P_{O}=(x_{o},0)$ represents the center of the reference circle, and $x_{o}$ can be obtained from the coordinates of $P_{A},P_{B}$, and $P_{C}$:

$$x_{o}=\frac{x^{2}_{r}-x^{2}_{t}-y^{2}_{t}}{2(x_{r}-x_{t})}.$$ (53)

Afterwards, the difference between the shaded left-hand point $P_{D}=(x_{d},0)$ and the reference circle left-hand point $P_{\mathcal{D}}=(\tilde{x}_{r},0)$ needs to be taken into account, with the size of this difference being evaluated by $\mathrm{D}=|x_{d}-\tilde{x}_{r}|$. Further, the two observable measurements take the form

$$R_{s}= \frac{(x_{t}-x_{r})^{2}+y_{t}^{2}}{2|x_{r}-x_{t}|},$$ (54)

$$\delta_{s}= \frac{\mathrm{D}}{R_{s}}.$$ (55)

Figure 7: observable measurements $R_{s}$ (left) and $\delta_{s}$ (right) as functions of $\alpha$ are presented, respectively, for $\Lambda=10^{-4}$ (top) and $\Lambda=10^{-3}$ (bottom). In all cases, the observer is placed at $r_{0}=50$ and $\vartheta_{0}=\pi/2.$

We numerically calculated these two observable measurements under the condition of fixing the observer position at $r_{0}=50$ with an inclination angle of $\vartheta_{0}=\pi/2$. The results are presented in Fig. 7. Note that the static case is not included in the results because the contour of the shadow coincides with the reference circle, resulting in no distortion. Figs. 7-(b) and 7-(d) show that the cosmological constant $\Lambda$ has little effect on $\delta_{s}$. But $\delta_{s}$ decreases as the parameter $\alpha$ increases, which suggests that the Kerr-MOG-dS black hole gets less deformed compared to the Kerr-dS black hole. Figs. 7-(a) and 7-(c) clearly show the radius $R_{s}$ of the shadow. From these results, it becomes evident that when the cosmological constant $\Lambda=10^{-4}$, $R_{s}$ increases as the parameter $\alpha$ increases with $\alpha\in[0,3]$. However, when the cosmological constant $\Lambda=10^{-3}$, there is a maximum value for $R_{s}$ , corresponding to $\alpha\approx 1.815$. This is consistent with the peculiar phenomena observed in Figs. 1-(c), 3-(c), and 3-(d) describing the shadow cast by MOG-dS black holes in Sec. IV.

For a approximatively estimation, utilizing the Kerr-MOG-dS metric (22), we calculate the angular radius of a black hole shadow, defined as $\theta_{\text{BH}}=\bar{R}_{s}\frac{\mathcal{M}}{D_{O}}$, with $D_{O}$ representing the distance from the observer to the black hole. The radius $\bar{R}_{s}$, which is proportional to the on-screen radius $R_{s}$ as described by the geometric relationships in equations (47) and (48) (see detail in footnote 1). Specifically, for a black hole possessing a mass of $\mathcal{M}$ and situated $D_{O}$ away from the observer, we express the angular radius $\theta_{\text{BH}}$ quantitatively as $\theta_{\text{BH}}=9.87098\times 10^{-6}\bar{R}_{s}\left(\frac{\mathcal{M}}{M_{\odot}}\right)\left(\frac{1\text{kpc}}{D_{O}}\right)\mu\text{as}$ Amarilla and Eiroa (2012).

$\theta_{\text{BH}}(\mu arcsec)$	$\alpha=0$	$\alpha=0.1$	$\alpha=0.2$	$\alpha=0.3$	$\alpha=0.4$
$\Lambda=0.0000$	24.737	26.747	28.738	30.712	32.670
$\Lambda=0.0001$	23.644	25.562	27.460	29.342	31.208
$\Lambda=0.0002$	22.497	24.317	26.118	27.902	29.671
$\Lambda=0.0003$	21.288	23.004	24.701	26.382	28.047
$\Lambda=0.0000$	24.751	26.760	28.750	30.722	32.679
$\Lambda=0.0001$	23.658	25.573	27.470	29.351	31.215
$\Lambda=0.0002$	22.510	24.328	26.127	27.909	29.677
$\Lambda=0.0003$	21.298	23.013	24.709	26.388	28.053

Table 1: The numerical estimation of the angular radius of the supermassive black hole Sgr A* in our galaxy using the metric of a Kerr-MOG-dS black hole. The top four rows correspond to $a/M=0.1$; the bottom four rows to $a/M=0.99$.

$\theta_{\text{BH}}(\mu arcsec)$	$\alpha=0$	$\alpha=0.1$	$\alpha=0.2$	$\alpha=0.3$	$\alpha=0.4$
$\Lambda=0.0000$	19.859	21.473	23.081	24.656	26.229
$\Lambda=0.0001$	18.982	20.522	22.046	23.556	25.055
$\Lambda=0.0002$	18.062	19.522	20.968	22.401	23.821
$\Lambda=0.0003$	17.090	18.457	19.816	21.164	22.517
$\Lambda=0.0000$	19.871	21.484	23.064	24.664	26.236
$\Lambda=0.0001$	18.993	20.531	22.054	23.563	25.061
$\Lambda=0.0002$	18.071	19.531	20.975	22.406	23.826
$\Lambda=0.0003$	17.099	18.475	19.837	21.185	22.521

Table 2: The numerical estimation of the angular radius of the supermassive black hole M87* using the metric of a Kerr-MOG-dS black hole. The top four rows correspond to $a/M=0.1$; the bottom four rows to $a/M=0.99$.

In Tables 1 and 2, we present the calculated angular radius of the black holes at the Galactic center, Sgr A*, and M87*, respectively, using the Kerr-MOG-dS metric. Here, we utilize the latest observations which indicate that the mass of the black hole Sgr A* is $\mathcal{M}=4.0\times 10^{6}M_{\odot}$ with an observer distance of $D_{O}=8.3kpc$ Akiyama et al. (2022b), and for the black hole M87*, the mass is $\mathcal{M}=6.5\times 10^{9}M_{\odot}$ with an observer distance of $D_{O}=16.8Mpc$ Akiyama et al. (2019b). It should be noted that we have highlighted the data corresponding to the observed range of angular diameters for the Sgr A* and M87* black holes from recent observations. By combining it with the data in Tables 1 and 2, one can find that there is room for the theoretical model of such a Kerr-MOG-dS black hole. Furthermore, the inclusion of the cosmological constant provides a larger parameter space for the MOG parameters.

V CONCLUSIONS AND OUTLOOKS

In this work, we have extend the Kerr-MOG black hole solution Moffat (2015a) to more general cases with a cosmological constant (namely, the Kerr-MOG-dS and Kerr-MOG-AdS black hole solutions) and investigate the black hole shadows of the Kerr-MOG-dS case. Due to the presence of the positive constant $\Lambda$, the observer should be positioned between the black hole’s event horizon and the cosmological horizon, rather than at infinity. The apparent shape of the shadow depends on the observer’s position.

A summary of three novel consequences of the shadows of the Kerr-MOG-dS black hole obtained in this paper are listed in the following order:

• Size: As the cosmological constant increases, the size of the black hole shadow decreases. But the size of the shadow does not appear to be monotonic as $\alpha$ changes. For example, it would reach a maximum at $\alpha\approx 1.815$ with $\Lambda=10^{-3}$.

• Shape: The shape of the shadow becomes more rounded as the MOG parameter $\alpha$ increases, causing the contours of rotating black holes to approach those of static black holes. However, this change does not impact the frame-dragging effect induced by rotation. The distortion of space is consistent with that observed in GR.

• Observable measurements: For observable measurements, we consider the radius $R_{s}$, which is associated with the apparent size, and the distortions $\delta_{s}$, which relate to the deformation of the shadow. The results show that, with the position of the observer held constant and the cosmological constant sufficiently large, the radius $R_{s}$ reaches a maximum value as $\alpha$ increases. This increase in $\alpha$ concurrently reduces the distortion $\delta_{s}$.

Our results indicate that the interplay between the cosmological constant and MOG parameter leads to intriguing phenomena in the black hole shadow contours. If the Earth were situated at the edge of the universe in this black hole spacetime, the black hole spin parameter and MOG parameter would become degenerate. Additionally, we used the latest observational data from M87* and Sgr A* to impose certain parameter constraints on Kerr-MOG-dS black holes, confirming the validity of STVG theory.

There are two promising further topics to be pursued in the future. An intriguing topic is to explore other observable measurements of the shadow of the Kerr-MOG-dS black hole. To prevent certain shadow degeneracy caused by black hole parameters, we can explore other distortion parameters, such as the oblateness parameter $K_{s}$, and the ”thickness” parameter $T$ of the shadow Wang et al. (2018b). In order to visualize the shadow in a real astronomical setting, one should rely on the emissions originating from the accretion disk surrounding the black hole Zhang et al. (2021, 2022). And the other is to investigate the properties of the Kerr-MOG-AdS black hole, such as black hole thermodynamics Christodoulou (1970); Bardeen et al. (1973); Wu and Wu (2019); Cong et al. (2021); Ahmed et al. (2023a), thermodynamic topological classification Wei et al. (2022); Wu (2023a); Wu and Wu (2023); Wu (2023b), null hypersurface caustics Al Balushi and Mann (2019); Imseis et al. (2021); Wu et al. (2021), phase transition criticality Kubiznak and Mann (2012); Ahmed et al. (2023b); Wu et al. (2024), Joule-Thomson expansion Ökcü and Aydıner (2018), etc.

Appendix A Energy-momentum tensor

In this appendix, we provide the energy-momentum tensor corresponding to the metric (II.2). The MOG vector potential (23) in the coordinates $(\bar{t},\bar{r},\bar{\vartheta},\bar{\varphi})$ is given as follows:

$$\bar{\phi}= -\frac{\sqrt{G_{\text{N}}\alpha}M\bar{r}}{\Sigma}\left(\frac{\sin^{2}\bar{\vartheta}}{\chi}+\cos^{2}\bar{\vartheta}\right)d\bar{t}-\frac{h\sqrt{G_{\text{N}}\alpha}M\bar{r}}{\Delta_{r}\Delta_{\bar{r}}}d\bar{r}+\frac{\sqrt{G_{\text{N}}\alpha}M\bar{r}}{\chi\Sigma}a\sin^{2}\bar{\vartheta}d\bar{\varphi},$$ (56)

All non-zero components of $T_{ab}$ are then given as follows:

$$T_{00}= -\frac{1}{24\pi\chi^{2}\Sigma^{3}}G_{\text{N}}M^{2}\alpha\Delta_{\bar{\vartheta}}\left[(3-\bar{r}^{2}\Lambda)\Sigma_{2}+3h\Delta_{\bar{\vartheta}}\right], T_{01}=-\frac{1}{8\pi\chi\Delta_{\bar{r}}\Sigma^{2}}G_{\text{N}}M^{2}\alpha h\Delta_{\bar{\vartheta}},$$ (57)

$$T_{11}= -\frac{1}{8\pi\Delta_{\bar{r}}^{2}\Sigma}G_{\text{N}}M^{2}\alpha\left(\Delta_{r}-2\Delta_{\bar{r}}\right), T_{03}=\frac{1}{8\pi\chi^{2}\Sigma^{3}}G_{\text{N}}M^{2}\alpha\Delta_{\bar{\vartheta}}\left(\Delta_{r}+\Delta_{\bar{r}}\right)a\sin^{2}\bar{{\vartheta}},$$

$$T_{22}= -\frac{1}{8\pi\Delta_{\bar{\vartheta}}\Sigma}G_{\text{N}}M^{2}\alpha, T_{13}=\frac{1}{8\pi\chi\Delta_{\bar{r}}\Sigma^{2}}G_{\text{N}}M^{2}\alpha ha\sin^{2}\bar{\vartheta},$$

$$T_{33}= -\frac{\sin^{2}\bar{{\vartheta}}}{8\pi\chi^{2}\Sigma^{3}}G_{\text{N}}M^{2}\alpha\left[(\bar{r}^{2}+a^{2})\Sigma_{2}+ha^{2}\sin^{2}\bar{\vartheta}\right].$$

This set of equations corresponds one-to-one with the effective tensor component equations (LABEL:E00). It is clear that this will lead to the effective tensor $E_{ab}=0$, ensuring that the solutions satisfy the field equations.