-
From the prediction of black holes by general relativity to the release of the first image of a black hole shadow in Virgo A* galaxy (M87) by the Event Horizon Telescope Collaboration in April 2019 [1], more than a hundred years have passed. Although we have made significant advancements in both observation and black hole theory, the unavoidable issue of singularities in black hole theory still remains. The singularity of a black hole has a series of peculiar properties, including infinite density and curvature, where currently known physical laws fail [2]. The black hole singularity problem originates from general relativity, but general relativity does not apply at the singularity, highlighting the incomplete development of our current theory of gravity. It is generally believed that solving this problem requires a complete theory of quantum gravity, through which the black hole singularity can be fully resolved. However, given that a complete theory of quantum gravity has not yet been discovered, using semiclassical theories to address the singularity problem remains extremely valuable.
In 1968, Bardeen presented the first solution for a regular black hole [3]. It avoids the singularity by incorporating a de Sitter core that generates negative pressure, thereby preventing a singular end-state of the gravitationally collapsed matter [4, 5]. Afterwards, Ayón-Beato and García considered the coupling between the gravitational field and a non-linear electromagnetic field [6]. Then, the Bardeen black hole solution could be recovered, indicating that the Bardeen model can be interpreted as a non-linear magnetic monopole in non-linear electrodynamics. Subsequently, more regular black holes, such as Ayón-Beato and García [7], Hayward [8], and Berej-Matyjasek-Trynieki-Wornowicz [9] black holes, were obtained by coupling non-linear electrodynamics with gravitational theory.
In astrophysical research, the spherically symmetric model is one of the simplest models. It can be divided into isotropic and anisotropic spherically symmetric celestial models based on whether the pressure is isotropic or anisotropic. Currently, most research focuses on isotropic fluids because astrophysical observations support isotropy [10, 11]. For example, the perfect Pascalian-fluid (isotropic-fluid) assumption is supported by theory and observation [12–14].
However, compared to isotropic models, studies on anisotropic models are relatively limited. Increasing theoretical evidence suggests that various interesting physical phenomena may occur within certain density ranges, leading to local anisotropy [15]. This has led to increased research on anisotropy. For instance, at the galactic level, dark force effects allow for an effective description in terms of general relativity sourced by an anisotropic fluid [16]. Local anisotropy in self-gravitating systems affects the physical properties, stability, and structure of stellar matter [17–19].
The electromagnetic field outside the Reissner-Nordström (R-N) black hole is the simplest anisotropic matter field, where the equation of state satisfies
$ p_{r}=-\varepsilon $ in the radial direction and$ p_{\theta}=p_{\phi}=\varepsilon $ in the angular directions within a local reference frame. Cho and Kim extended the R-N black hole solution to obtain a Schwarzschild black hole with an anisotropic matter field [20]. The anisotropic matter field in this black hole is described by the parameters ω and K, related as$ r_{0}^{2 \omega _{2}} \equiv \left(1-2 \omega _{2}\right) K \geq 0 $ . Subsequently, a rotating charged black hole with an anisotropic matter field [21] was proposed, and its geometry, thermodynamic properties, and energy extraction were studied. For further research related to anisotropic matter fields, see Refs. [22–27].We have reviewed research and developments on regular black holes and anisotropic matter fields. In this paper, we consider both to obtain the solution of a non-linear magnetic-charged black hole with an anisotropic matter field and the corresponding rotating black hole. We then further study the geometric and thermodynamic properties of this rotating black hole.
The structure of this paper is as follows. In Sec. II, we obtain the solution for a non-linear magnetic-charged black hole with an anisotropic matter field. In Sec. III, we extend this solution using the modified Newman-Janis algorithm to obtain a solution for a rotating and non-linear magnetic-charged black hole with an anisotropic matter field and study its energy-momentum tensor and weak energy condition. In Secs. IV and V, we explore the geometric and thermodynamic properties of the rotating black hole, respectively. In Sec. VI, we examine the efficiency of the Penrose process for extracting energy from the rotating black hole. Finally, in Sec. VII, we present the conclusion and discussion of this study.
-
We consider the coupling between Einstein gravity and a non-linear electromagnetic field in the presence of an anisotropic matter field, which is described by the following equations:
$ G_{\mu}^{\nu}=2\left[\frac{\partial {\cal{L}}(F)}{\partial F} F_{\mu \rho} F^{v \rho}-\delta_{\mu}^{\nu} {\cal{L}}(F)\right] +8\pi T_{\mu}^{\nu}, $
(1) $ \nabla_{\mu}\left(\frac{\partial {\cal{L}}(F)}{\partial F} F^{v \mu}\right)=0 , $
(2) $ \nabla_{\mu} * F^{\nu \mu}=0. $
(3) Here,
$ {\cal{L}} $ is a function of the invariant$ F \equiv \frac{1}{4} F_{\mu \nu} F^{\mu \nu} $ with$ F_{\mu \nu}=\partial_{\mu} A_{\nu}-\partial_{\nu} A_{\mu} $ , given by [5, 28]$ {\cal{L}}(F)=\frac{3 M}{|Q|^{3}} \frac{\left(2 Q^{2} F\right)^{3 / 2}}{\left[1+\left(2 Q^{2} F\right)^{3 / 4}\right]^{2}}, $
(4) where Q is the magnetic charge. In this study, we consider black holes with an anisotropic matter field. The energy-momentum tensor of an anisotropic matter field obtained from Ref. [20] is given by
$ T_{\nu}^{\mu}={\rm{diag}}\left(-\varepsilon, p_{r}, p_{\theta}, p_{\phi}\right), $
(5) where
$ p_{r}=-\varepsilon(r) $ and$ p_{\theta}=p_{\phi}=\omega \varepsilon(r) $ . The energy density is$ \varepsilon(r)=\frac{r_{0}^{2 \omega }}{8 \pi r^{2 \omega +2}}, $
(6) where
$ r_0 $ is a charge-like quantity of dimension of length defined by$ r_{0}^{2 \omega }=(1-2 \omega ) K $ .To obtain a static spherically symmetric black hole solution, we assume that the line element is
$ {\rm d} s^{2}=-f(r) {\rm d} t^{2}+f(r)^{-1} {\rm d} r^{2}+r^{2} {\rm d} \Omega^{2}, $
(7) $ f(r)=1-\frac{2 m(r)}{r} \text { and } {\rm d} \Omega^{2}={\rm d} \theta^{2}+\sin ^{2} \theta {\rm d} \phi^{2}, $
(8) and use the ansatz for Maxwell field,
$ F_{\mu \nu}=\left(\delta_{\mu}^{\theta} \delta_{\nu}^{\varphi}-\delta_{\nu}^{\theta} \delta_{\mu}^{\varphi}\right) B(r, \theta). $
(9) Using Eqs. (2) and (3), we can simplify Eq. (9) to obtain
$ F_{\mu \nu}=\left(\delta_{\mu}^{\theta} \delta_{\nu}^{\varphi}-\delta_{\nu}^{\theta} \delta_{\mu}^{\varphi}\right) Q \sin \theta, $
(10) where Q is the integration constant. Furthermore, we obtain
$ F=Q^{2} / 2 r^{4} $ . Note that the magnetic charge Q is defined in the following integral form:$ \frac{1}{4 \pi} \int_{S_{2}^{\infty}} {\boldsymbol{F}}=Q, $
(11) where
$ S_{2}^{\infty} $ is a two-sphere at the infinity. From the above results, it is possible to simplify the time component of Eq. (1) as$ \frac{2}{r^{2}} \frac{{\rm{d}} m(r)}{{\rm{d}} r}=\frac{6 M Q^{3}}{\left(r^{3}+Q^{3}\right)^{2}}+\frac{(1-2\omega )K}{2r^{2\omega+2 }}. $
(12) By integrating Eq. (12) from r to
$ \infty $ and using the relation$ M=\lim\limits_{r \rightarrow \infty}\left(m(r)-\dfrac{K}{2} r^{-2\omega +1}\right) $ , one finally obtains$ f(r)=1-\frac{2 M r^{2}}{r^{3}+Q^{3}}-\frac{K}{r^{2\omega }}. $
(13) For the above equation for
$ K=0 $ , we can derive the metric function$ f(r) $ of the Hayward-like black hole [8]. When$ Q=0 $ and$ \omega=1 $ , we can derive the metric function$ f(r) $ of the R-N black hole. It is found that K and q (q is the charge in the R-N black hole) have similar status in the metric function$ f(r) $ . Therefore, K is the global charge, and$ r_0 $ is a charge-like quantity of dimension of length defined by$ r_{0}^{2\omega} = (1 - 2\omega)K $ .In the introduction, we mentioned the non-linear magnetic-charged black hole, which is a regular black hole without a curvature singularity. Now, we consider whether an anisotropic matter field in the non-linear magnetic-charged black hole will alter this property. To investigate this, we calculate the curvature of the black hole as follows:
$ R=\frac{12 M Q^3 \left(2 Q^3-r^3\right)}{\left(Q^3+r^3\right)^3}+2 K (\omega -1) (2 \omega -1) r^{-2 (\omega +1)}, $
(14) $ \begin{aligned}[b] R_{\mu \nu} R^{\mu \nu}=\;&\frac{72 M^2 Q^6 \left(2 Q^6-2 Q^3 r^3+5 r^6\right)}{\left(Q^3+r^3\right)^6}\\ &+\frac{24 K M Q^3 (2 \omega - 1) r^{-2 (\omega +1)} \left(Q^3 (\omega - 1)-r^3 (2 \omega + 1)\right)}{\left(Q^3+r^3\right)^3}\\ &+2 \left(\omega ^2+1\right) (K-2 K \omega )^2 r^{-4 (\omega +1)}, \end{aligned} $
(15) $ \begin{aligned}[b] R_{\mu \nu \sigma \tau} R^{\mu \nu \sigma \tau} =\;& 4K^2 \left(4 \omega ^4+4 \omega ^3+5 \omega ^2+1\right) r^{-4 (\omega +1)} \\ & + \frac{16 K M }{\left(Q^3+r^3\right)^3 r^{2 (\omega +1)}} \left[Q^6 \left(2 \omega ^2-3 \omega +1\right) \right. \\ &\left. + Q^3 r^3 \left(-14 \omega ^2-9 \omega +2\right) + r^6 \left(2 \omega ^2+3 \omega +1\right)\right] \\ & + \frac{48 M^2 \left(2 Q^{12} - 2 Q^9 r^3 + 18 Q^6 r^6 - 4 Q^3 r^9+r^{12}\right)}{\left(Q^3+r^3\right)^6}. \end{aligned} $
(16) From the above results, we find that the curvature diverges at
$ r=0 $ owing to the influence of the anisotropic matter field, which indicates that the non-linear magnetic-charged black hole with an anisotropic matter field is not a regular black hole and has a singularity. -
In 1965, the Newman-Janis algorithm (NJA) [29] was first proposed by Newman and Janis. This algorithm obtains the Kerr metric by performing a complex coordinate transformation on the Schwarzschild line element. The NJA has been widely applied in the scientific community [21, 28, 30–38]. Later, Azreg-Aïnou modified the NJA [39, 40]; the only difference from the original one is that they skipped the complexification of coordinates.
In this section, we use the modified NJA to generalize the solution of the spherically symmetric non-linear magnetic-charged black hole with an anisotropic matter field to obtain the corresponding rotating black hole solution. First, let us review the modified NJA for deriving a general rotating black hole line element from a general static spherically symmetric black hole line element. The line element for a static spherically symmetric spacetime in Boyer-Lindquist (BL) coordinates is
$ {\rm d} s^{2}=-f(r) {\rm d} t^{2}+g^{-1}(r) {\rm d} r^{2}+h(r) {\rm d} \Omega^{2}, $
(17) where
$ {\rm d} \Omega^{2}={\rm d} \theta^{2}+\sin ^{2} \theta {\rm d} \phi^{2}. $
(18) Considering the coordinate transformation
$ {\rm d} u={\rm d} t-\frac{{\rm d} r}{\sqrt{f(r) g(r)}}, $
(19) the line element in Eddington-Finkelstein (EF) coordinates
$ (u, r, \theta, \phi) $ is obtained as$ {\rm d} s^{2}=-f(r) {\rm d} u^{2}-2 \sqrt{\frac{f(r)}{g(r)}} {\rm d} u {\rm d} r+h(r) {\rm d} \Omega^{2}. $
(20) Concerning the null tetrads that satisfy the relations
$ l_{\mu} l^{\mu}= n_{\mu} n^{\mu}=m_{\mu} m^{\mu}=l_{\mu} m^{\mu}= n_{\mu} m^{\mu}=0 $ and$ l_{\mu} n^{\mu}=-m_{\mu} \bar{m}^{\mu}= 1 $ , the contravariant metric tensor in Eq. (20) can be expressed in terms of the null tetrads as$ g^{\mu \nu}=-l^{\mu} n^{\nu}-l^{\nu} n^{\mu}+m^{\mu} \bar{m}^{\nu}+m^{\nu} \bar{m}^{\mu}, $
(21) where
$ \begin{aligned}[b] l^{\mu} & = \delta_{r}^{\mu}, \\ n^{\mu} & = \sqrt{\frac{g(r)}{f(r)}} \delta_{r}^{\mu} - \frac{f(r)}{2} \delta_{t}^{\mu}, \\ m^{\mu} & = \frac{1}{\sqrt{2 h(r)}} \delta_{\theta}^{\mu} + \frac{\rm i}{\sqrt{2 h(r)} \sin \theta} \delta_{\phi}^{\mu}, \end{aligned} $
$ \begin{aligned}[b] \bar{m}^{\mu} = \frac{1}{\sqrt{2 h(r)}} \delta_{\theta}^{\mu} - \frac{\rm i}{\sqrt{2 h(r)} \sin \theta} \delta_{\phi}^{\mu}. \end{aligned} $
(22) The key step of the modified NJA is the consideration of complex coordinate transformations in the
$ u-r $ plane,$ u \rightarrow u-{\rm i} a \cos \theta, \quad r \rightarrow r+{\rm i} a \cos \theta. $
(23) As a result of the transformation of Eq. (22), the metric function also takes a new form:
$ f(r) $ $ \rightarrow $ $ F(r, a, \theta) $ ,$ g(r) $ $ \rightarrow G(r, a, \theta) $ , and$ h(r) $ $ \rightarrow $ $ \Sigma=r^{2}+a^{2} \cos ^{2}\theta $ [39, 40]. Thus, the new form of null tetrads is given by$ \begin{aligned}[b] l^{\mu} & =\delta_{r}^{\mu}, \\ n^{\mu} & =\sqrt{\frac{G}{F}} \delta_{u}^{\mu}-\frac{1}{2} F \delta_{r}^{\mu}, \\ m^{\mu} & =\frac{1}{\sqrt{2 \Sigma}}\left[\delta_{\theta}^{\mu}+{\rm i} a\left(\delta_{u}^{\mu}-\delta_{r}^{\mu}\right) \sin \theta+\frac{\rm i}{\sin \theta} \delta_{\phi}^{\mu}\right], \\ \bar{m}^{\mu} & =\frac{1}{\sqrt{2 \Sigma}}\left[\delta_{\theta}^{\mu}-{\rm i} a\left(\delta_{u}^{\mu}-\delta_{r}^{\mu}\right) \sin \theta-\frac{\rm i}{\sin \theta} \delta_{\phi}^{\mu}\right]. \end{aligned} $
(24) Then, from Eq. (21), the contravariant components of the metric
$ g^{\mu\nu} $ can be obtained as$ \begin{aligned}[b] g^{u u} & =\frac{a^{2} \sin ^{2} \theta}{\Sigma}, \\ g^{u r} & =-\sqrt{\frac{G}{F}}-\frac{a^{2} \sin ^{2} \theta}{\Sigma}, \\ g^{u \phi} & =\frac{a}{\Sigma}, \\ g^{r r} & =G+\frac{a^{2} \sin ^{2} \theta}{\Sigma}, \\ g^{r \phi} & =-\frac{a}{\Sigma}, \\ g^{\theta \theta} & =\frac{1}{\Sigma}, \\ g^{\phi \phi} & =\frac{1}{\Sigma \sin ^{2} \theta}. \end{aligned} $
(25) The following transformation is used to convert from EF to BL coordinates:
$ {\rm d} u={\rm d} t+\lambda(r) {\rm d} r, \quad {\rm d} \varphi={\rm d} \phi+\chi(r) {\rm d} r, $
(26) where the functions
$ \lambda(r) $ and$ \chi(r) $ can be determined by ensuring that all the nondiagonal components of the metric tensor other than the coefficient$ g_{t \phi}\left(g_{\phi t}\right) $ are equal to zero. Therefore, we obtain that [39]$ \lambda(r)=-\frac{k(r)+a^{2}}{g(r) h(r)+a^{2}}, \quad \chi(r)=-\frac{a}{g(r) h(r)+a^{2}}, $
(27) with
$ k(r)=\sqrt{\frac{g(r)}{f(r)}} h(r), $
(28) and
$ \begin{aligned}[b] F(r, \theta)&=\frac{\left(g h+a^{2} \cos ^{2} \theta\right) \Sigma}{\left(k^{2}+a^{2} \cos ^{2} \theta\right)^{2}}, \\ G(r, \theta)&=\frac{g h+a^{2} \cos ^{2} \theta}{\Sigma}. \end{aligned} $
(29) Finally, the rotating solution corresponding to the spherically symmetric line element Eq. (17) is obtained as follows:
$ \begin{aligned}[b] {\rm d} s^{2}= & -\frac{\left(g h+a^{2} \cos ^{2} \theta\right) \Sigma}{\left(k+a^{2} \cos ^{2} \theta\right)^{2}} {\rm d} t^{2}+\frac{\Sigma}{g h+a^{2}} {\rm d} r^{2} \\ & -2 a \sin ^{2} \theta\left[\frac{k-g h}{\left(k+a^{2} \cos ^{2} \theta\right)^{2}}\right] \Sigma {\rm d} \phi {\rm d} t+\Sigma {\rm d} \theta^{2} \\ & +\Sigma \sin ^{2} \theta\left[1+a^{2} \sin ^{2} \theta \frac{2 k-g h+a^{2} \cos ^{2} \theta}{\left(k+a^{2} \cos ^{2} \theta\right)^{2}}\right] {\rm d} \phi^{2}. \end{aligned} $
(30) Above, we mentioned the modified NJA algorithm, which provides a general expression for a rotating black hole metric. Using the metric function from Eq. (8) in the previous section, we obtain
$ \begin{aligned}[b] g(r)&=f(r)=1-\frac{2 M r^{2}}{r^{3}+Q^{3}}-\frac{K}{r^{2\omega }}, \\ h(r)&=k(r)=r^{2}. \end{aligned} $
(31) Substituting Eq. (31) into Eq. (30), we obtain the line element of a rotating and non-linear magnetic-charged black hole with an anisotropic matter field as
$ \begin{aligned}[b] {\rm d} s^{2}= & -\left[1-\frac{2 \rho r}{\Sigma}\right] {\rm d} t^{2}+\frac{\Sigma}{\Delta} {\rm d} r^{2}-\frac{4 a \rho r \sin ^{2} \theta}{\Sigma} {\rm d} t {\rm d} \phi \\ & +\Sigma {\rm d} \theta^{2}+\sin ^{2} \theta\left[r^{2}+a^{2}+\frac{2 a^{2} \rho r \sin ^{2} \theta}{\Sigma}\right] {\rm d} \phi^{2}, \end{aligned} $
(32) where
$ \begin{aligned}[b] \Delta & =r^{2}-2 \rho r+a^{2}, \\ \Sigma & =r^{2}+a^{2} \cos ^{2} \theta, \\ 2 \rho & =\frac{2 M r^{3}}{r^{3}+Q^{3}}+K r^{-2\omega +1}. \end{aligned} $
(33) -
The nonvanishing components of the Einstein tensor
$ G_{\mu \nu} $ are given by$ \begin{aligned}[b] G_{t t} =\;&\frac{2\left(r^{4}+a^{2} r^{2}-2 r^{3} \rho-a^{4} \sin ^{2} \theta \cos ^{2} \theta\right) \rho^{\prime}}{\Sigma^{3}} \\ & -\frac{a^{2} r \sin ^{2} \theta \rho^{\prime \prime}}{\Sigma^{2}}, \\ G_{t \phi} =\;&\frac{2 a \sin ^{2} \theta\left[\left(a^{2}+r^{2}\right)\left(a^{2} \cos ^{2} \theta-r^{2}\right)+2 r^{3} \rho\right] \rho^{\prime}}{\Sigma^{3}} \\ & +\frac{a r \sin ^{2} \theta\left(a^{2}+r^{2}\right) \rho^{\prime \prime}}{\Sigma^{2}}, \\ G_{r r} =\;&-\frac{2 r^{2} \rho^{\prime}}{\Sigma \Delta}, \\ G_{\theta \theta} =\;&-\frac{2 a^{2} \cos ^{2} \theta \rho^{\prime}}{\Sigma}-r \rho^{\prime \prime}, \\ G_{\phi \phi} =\;&-\frac{a^{2} \sin ^{2} \theta\left(a^{2}+r^{2}\right)\left(a^{2}+\left(a^{2}+2 r^{2}\right) \cos 2 \theta\right) \rho^{\prime}}{\Sigma^{3}} \\ & -\frac{4 a^{2} r^{3} \sin ^{4} \theta \rho \rho^{\prime}}{\Sigma^{3}}-\frac{r \sin ^{2} \theta\left(a^{2}+r^{2}\right)^{2} \rho^{\prime \prime}}{\Sigma^{2}}, \end{aligned} $
(34) where the prime
$ \prime $ denotes the derivative with respect to r and$ \begin{aligned}[b] 2 \rho & =\frac{2 M r^{3}}{r^{3}+Q^{3}}+K r^{-2\omega +1} ,\\ 2 \rho^{\prime} & =\frac{6 M Q^3 r^2}{\left(Q^3+r^3\right)^2}+K (1-2 \omega ) r^{-2 \omega }, \\ 2 \rho^{\prime \prime}&=\frac{12 M Q^3 r \left(Q^3-2 r^3\right)}{\left(Q^3+r^3\right)^3}-2 K \omega (1-2\omega ) r^{-2 \omega -1}. \end{aligned} $
(35) To obtain the components of the energy-momentum tensor, we use the standard orthonormal basis [28, 41]:
$ \begin{aligned}[b] e_{(t)}^{\mu} & =\frac{1}{\sqrt{\Delta \Sigma}}\left(a^{2}+r^{2}, 0,0, a\right), \\ e_{(r)}^{\mu} & =\frac{\sqrt{\Delta}}{\sqrt{\Sigma}}(0,1,0,0), \\ e_{(\theta)}^{\mu} & =\frac{1}{\sqrt{\Sigma}}(0,0,1,0), \\ e_{(\phi)}^{\mu} & =-\frac{1}{\sqrt{\Sigma \sin ^{2} \theta}}\left(a \sin ^{2} \theta, 0,0,1\right). \end{aligned} $
(36) The components of the energy-momentum tensor derived from Eqs. (34) and (36) and the Einstein field equation
$ G_{\mu \nu} = 8 \pi T_{\mu \nu} $ are as follows:$ \begin{aligned}[b] T_{(t)(t)} & =\frac{1}{2} e_{(t)}^{\mu} e_{(t)}^{\nu} G_{\mu \nu}=\epsilon , \quad T_{(\theta)(\theta)} =\frac{1}{2} e_{(\theta)}^{\mu} e_{(\theta)}^{\nu} G_{\mu \nu}=p_{\theta}, \\ T_{(r)(r)} &=\frac{1}{2} e_{(r)}^{\mu} e_{(r)}^{\nu} G_{\mu \nu}=p_{r}, \quad T_{(\phi)(\phi)} =\frac{1}{2} e_{(\phi)}^{\mu} e_{(\phi)}^{\nu} G_{\mu \nu}=p_{\phi}, \end{aligned} $
(37) from which we obtain
$ \epsilon = -p_{r}=T_{(t)(t)}=\frac{1}{4 \pi} \frac{r^{2} \rho^{\prime}}{\Sigma^{2}}, $
(38) $ p_{\theta} = p_{\phi}=T_{(\theta)(\theta)}= \omega \left(\epsilon-\frac{2 \rho^{\prime}+r \rho^{\prime \prime}}{8 \pi \Sigma}\right). $
(39) When
$ Q=0 $ and$ K=\omega =0 $ , Eq. (39) corresponds to rotating black holes with an anisotropic matter field and non-linear magnetic-charged black holes, respectively. -
The weak energy condition ensures that the behavior of matter and energy is physically reasonable, preventing phenomena such as negative energy density. It states that, for any timelike vector
$ \xi^{\mu} $ at any point in spacetime, the stress-energy tensor$ T_{\mu \nu} $ satisfies$ T_{\mu \nu} \xi^{\mu} \xi^{\nu} \geq 0. $
(40) With the decomposition of the energy-momentum tensor
$ T_{\mu \nu} $ , the weak energy condition is equivalent to$ \epsilon \geq 0, \quad \epsilon+p_{i} \geq 0, $
(41) where
$ i=r $ , θ, φ. Substituting Eqs. (33), (38), and (39) into Eq. (41), we obtain$ \epsilon =\frac{r^2 \left(K (1-2 \omega ) r^{-2 \omega }+\frac{6 M Q^3 r^2}{\left(Q^3+r^3\right)^2}\right)}{8\pi\left(r^2+ a^2 \cos ^2\theta\right)^2} \geq 0, $
(42) $ \epsilon+p_{r}=0 , $
(43) $ \begin{aligned}[b] \epsilon+p_{\theta}=\;&\epsilon+p_{\phi}= \omega \epsilon+\frac{3 M Q^3 r^2 \omega \left(r^3-2 Q^3\right)}{4\pi\left(Q^3+r^3\right)^3 \left(a^2 \cos ^2\theta +r^2\right)}\\ &-\frac{K \omega \left(2 \omega ^2-3 \omega +1\right) r^{-2 \omega }}{8\pi(a^2 \cos ^2\theta +r^2)}\geq 0. \end{aligned} $
(44) From Eqs. (42), (43), and (44), we present the variations of
$ \epsilon $ and$ \epsilon+p_{\theta}/\epsilon+p_{\varphi} $ with r and$ \cos\theta $ under two sets of parameters in Figs. 1 and 2. It can be concluded from Figs. 1 and 2 that the weak energy condition is violated in a rotating, non-linear magnetic-charged black hole with an anisotropic matter field.Figure 1. (color online) Dependence of matter density
$ \epsilon $ and$ \epsilon+p_{\theta} $ on radius and angle for rotating and nonlinear magnetic-charged black hole with an anisotropic matter field for$ M=1 $ ,$ Q=0.5 $ ,$ a=0.2 $ ,$ K=1 $ , and$ \omega =0.2 $ .Figure 2. (color online) Dependence of matter density
$ \epsilon $ and$ \epsilon+p_{\theta} $ on radius and angle for a rotating and nonlinear magnetic-charged black hole with an anisotropic matter field for$ M=1 $ ,$ Q=0.5 $ ,$ a=0.9 $ ,$ K=1 $ , and$ \omega =0.2 $ .Based on the analysis of parameter ω in Ref. [21], when
$ 0 \leq \omega \leq \dfrac{1}{2} $ , the energy density is not sufficiently localized, causing the total mass to diverge as$ r \to \infty $ . As a result, subsequent studies primarily focused on systems with$ \omega > \dfrac{1}{2} $ to ensure that the spacetime asymptotically approaches flatness at infinity. -
The event horizon corresponds to the Killing horizon, which defines the "surface" of a black hole as a three-dimensional hypersurface. The line element Eq. (32) is singular at
$ \Delta= 0 $ , which corresponds to the event horizon of the black hole and satisfies$ \Delta =r^{2}+a^{2}-\frac{2 M r^{4}}{r^{3}+Q^{3}}-K r^{2\left(1-\omega \right)}=0 . $
(45) Based on the above equation, Fig. 3 shows the behavior of Δ as a function of r for different values of the parameters a, K, and ω with Q taking values of
$ 0 $ ,$ 0.3 $ ,$ 0.5 $ , and$ 0.7 $ , respectively. In Fig. 3, the first and second graphs correspond to the case of$ K<0 $ , while the third and fourth graphs correspond to the case of$ K>0 $ . According to the first and third graphs in Fig. 3, for$ \dfrac{1}{2} < \omega < 1 $ , the influence of Q on Δ is weakened as r increases, and this influence becomes negligible when r is large. This effect is noticeable within a small interval starting from$ 0 $ . According to the second and fourth graphs in Fig. 3, for$ \omega > 1 $ , the influence of Q on Δ is insignificant for both small and large values of r, with a noticeable effect only within a small interval of r. Therefore, the influence of Q on Δ mainly occurs within a small interval. The above analysis shows that we can analyze the possible roots of Δ by considering the case when$ Q = 0 $ [21]. By analyzing Eq. (45) and Fig. 3, we draw the following conclusions:Figure 3. (color online) Behavior of Δ as a function of r for different values of the parameters a, K, and ω when Q takes the values of 0, 0.3, 0.5, and 0.7.
● For
$ \dfrac{1}{2} < \omega < 1 $ , Δ is mainly dominated by the$ r^2 $ term for large r, while for smaller r, it is mainly influenced by$ r^2 $ and the non-linear magnetic-charged term. The influence of the nonlinear magnetic charge term causes Δ ($ Q \neq 0 $ ) to deviate from the case Δ ($ Q = 0 $ ) in a small interval starting from$ r = 0 $ . For both cases,$ K < 0 $ and$ K > 0 $ , the rotating black hole may have two horizons, namely the Cauchy and event horizons. For example, in the first and third graphs in Fig. 3, the Δ graphs correspond to the cases of$ K = -0.1 $ ,$ \omega = \dfrac{2}{3} $ and$ K = 0.1 $ ,$ \omega = \dfrac{2}{3} $ , respectively.● For
$ \omega>1 $ , Δ is mainly dominated by the$ r^2 $ term for large r, while for smaller r, it is primarily influenced by the anisotropic matter field term. The influence of the non-linear magnetic-charged term causes Δ ($ Q \neq 0 $ ) to deviate within a small range of r compared to the case Δ ($ Q = 0 $ ). For$ K < 0 $ , Δ starts from positive infinity at$ r \to 0 $ , decreases rapidly, and then increases again. Therefore, for$ \omega >1 $ and$ K < 0 $ , the black hole may have two horizons. For$ K > 0 $ , Δ starts from negative infinity at$ r \to 0 $ , increases rapidly, then decreases, and then increases again, or Δ starts from negative infinity at$ r \to 0 $ , increases rapidly, and then increases with the$ r^2 $ term. Therefore, for$ \omega >1 $ and$ K > 0 $ , the black hole may have three horizons. For example, in the second and fourth graphs in Fig. 3, they correspond to Δ images for$ K = -0.1 $ ,$ \omega = \dfrac{3}{2} $ and$ K = 0.018 $ ,$ \omega = \dfrac{3}{2} $ , respectively. -
The stationary limit surface, also known as the infinite redshift surface, is where the timelike Killing vector
$ K^\mu = \partial_t $ satisfies$ K^\mu K_\mu = 0 $ . Therefore, we obtain$ r^{2}+a^{2}\cos ^2\theta-\frac{2 M r^{4}}{r^{3}+Q^{3}}-K r^{2\left(1-\omega \right)}=0 . $
(46) The ergosphere is a region outside the black hole, bounded by the event horizon and the outer stationary limit surface. Energy extraction from a rotating black hole relies on the ergosphere, which corresponds to the area between the outermost red curve and the blue curve in each graph in Fig. 4. Using Eqs. (45) and (46), we present the behavior of the ergosphere in the xz-plane for different values of the parameters a, K, ω, and Q in Fig. 4. In the first row of ergosphere graphs in Fig. 4, we show that, when other parameters are kept constant and the parameter Q is varied, the area of the ergosphere increases with the parameter Q. In the second and third rows of ergosphere graphs in Fig. 4, we show that, when other parameters are kept constant and parameter a is varied, the area of the ergosphere increases with the parameter a. In the third row of ergosphere graphs in Fig. 4, we show that, for
$ K > 0 $ and$ \omega >1 $ , as the parameter a increases, the black hole transitions from having a single horizon to three horizons, and then to two horizons. -
We now understand that a black hole is not just a gravitational system but also a thermodynamic system. This insight originates from the proposal of black hole entropy [42, 43] and Hawking temperature [44, 45]. The Bekenstein-Hawking entropy is one-quarter of the event horizon area. Therefore, the entropy of a rotating and non-linear magnetic-charged black hole with an anisotropic matter field is given by
$ S=\frac{{\cal{A}}}{4}=\pi\left(r_{H}^{2}+a^{2}\right), $
(47) where
$ r_{H} $ is the event horizon, which satisfies Eq. (45). On the event horizon, the Hawking temperature is proportional to the surface gravity κ:$ \begin{aligned}[b] T=\frac{\kappa}{2 \pi}=\frac{\partial_{r_{H}} \Delta\left(r_{H}\right)}{4 \pi\left(r_{H}^{2}+a^{2}\right)} \end{aligned} $
$ \begin{aligned}[b] =\;&\frac{1}{4 \pi r_{H} \left(a^2+r_{H}^2\right) \left(Q^3+r_{H}^3\right)}(-a^2 \left(4 Q^3+r_{H}^3\right)+r_{H}^5\\ &+K r_{H}^{-2 \omega } \left(2 Q^3 r_{H}^2 (\omega +1)+r_{H}^5 (2 \omega -1)\right)-2 Q^3 r_{H}^2). \end{aligned} $
(48) The angular velocity of a rotating black hole at the event horizon is
$ \Omega=-\left.\frac{g_{t \phi}}{g_{\phi \phi}}\right|_{r=r_{H}}=\frac{a}{r_{H}^{2}+a^{2}} . $
(49) Above, we directly derived the entropy, temperature, and angular velocity of the rotating black hole through the relationship between thermodynamic and classical mechanical quantities. However, other thermodynamic quantities of the black hole, such as the electric potential, are not easily obtained using similar methods. Below, we use the first law of thermodynamics and the squared-mass formula of a rotating and non-linear magnetic-charged black hole with an anisotropic matter field to solve for various thermodynamic quantities of the rotating black hole.
For a rotating and non-linear magnetic-charged black hole with an anisotropic matter field, we consider S, Q, J, and
$ r_{0} $ as four independent thermodynamic quantities and express its first law of thermodynamics as [46]$ {\rm d} M=T {\rm d} S+\Phi {\rm d} Q+\Omega {\rm d} J+\Phi_{0} {\rm d} r_{0} ,$
(50) where Φ is the magnetic charge potential for Q, J is the angular momentum, and
$ \Phi_{0} $ is the electric potential for$ r_0 $ . From Eqs. (45) and (48) and$ J=Ma $ , we obtain the squared-mass formula for a rotating and non-linear magnetic-charged black hole with an anisotropic matter field as [47, 48]$ \begin{aligned}[b] M^2=\;&\frac{\pi }{S}J^2+\frac{\pi}{4S}\left(\frac{(\dfrac{S}{\pi}-a)^{\frac{3}{2}}+Q^{3}}{(\dfrac{S}{\pi}-a^{2})^{\frac{3}{2}}}\right)^{2}\\ &\times \left(\frac{S}{\pi} -\frac{r_{0}^{2\omega }}{(1-2\omega )}\left(\frac{S}{\pi}-a^{2}\right)^{\frac{1}{2}\left(1+\frac{r_{0}^{2\omega }}{K}\right)}\right)^{2}. \end{aligned} $
(51) From Eqs. (50) and (51), we obtain the following two thermodynamic quantities:
$ \begin{aligned} \Phi=\left(\frac{\partial M}{\partial Q}\right)_{S,J,r_{0}} =\frac{3Q^2}{2r_{H}^2(r_{H}^2+a^2)}\left(r_{H}^2+a^2-Kr_{H}^{2(1-\omega )}\right), \end{aligned} $
(52) and
$ \begin{aligned}[b] \Phi_{0}&=\left(\frac{\partial M}{\partial r_{0}}\right)_{S,Q,J}\\ &=\frac{\omega }{2\omega -1}\left(\frac{r_{H}^{3}+Q^{3}}{r_{H}^{3}}\right)\left(\frac{r_{H}r_{0}}{r_{H}^{2}+a^{2}}\right)\left(\frac{r_{0}}{r_{H}}\right)^{2(\omega -1)}. \end{aligned} $
(53) For Eq. (53), when
$ Q=a=0 $ and$ Q =0 $ ($ a\ne 0 $ ), it degenerates into a Schwarzschild black hole with an anisotropic matter field [20] and rotating black hole with an anisotropic matter field [21], respectively. For a thermodynamic black hole system, the thermodynamic stability can be studied through its heat capacity. Therefore, we obtain the heat capacity of a rotating and non-linear magnetic-charged black hole with an anisotropic matter field and analyze the stability of the rotating black hole. The heat capacity C, calculated at constant angular momentum and charge in the canonical ensemble, is given by$ \begin{aligned}[b] C =\;& T\frac{\partial S}{\partial T} = 2 \pi \left(a^2+r_{H}^2\right) \left(Q^3+r_{H}^3\right) r_{H}^{2 \omega +2} [-a^2 \left(4 Q^3+r_{H}^3\right) \\ &+K r_{H}^{-2 \omega } \left(2 Q^3 r_{H}^2 (\omega +1)+r_{H}^5 (2 \omega -1)\right)-2 Q^3 r_{H}^2+r_{H}^5] \\ &/[a^4 \left(4 Q^6+14 Q^3 r_{H}^3+r_{H}^6\right) r_{H}^{2 \omega } \\ &+a^2 r_{H}^2 [2 \left(5 Q^6+16 Q^3 r_{H}^3+2 r_{H}^6\right) r_{H}^{2 \omega } \\ &-K [2 Q^6 \left(2 \omega ^2+\omega -1\right)+2 Q^3 r_{H}^3 \left(4 \omega ^2-\omega +4\right) \\ &+r_{H}^6 (1-2 \omega )^2]]-r_{H}^4 [K [Q^6 \left(4 \omega ^2+6 \omega +2\right) \\ &+2 Q^3 r_{H}^3 \left(4 \omega ^2+3 \omega +5\right)+r_{H}^6 \left(4 \omega ^2-1\right)] \\ &+\left(-2 Q^6-10 Q^3 r_{H}^3+r_{H}^6\right) r_{H}^{2 \omega }]]. \end{aligned} $
(54) When
$ Q = 0 $ , the above equation degenerates into the heat capacity of a rotating black hole with an anisotropic matter field [21]. Based on Eq. (54), we plot the behavior of heat capacity as a function of$ r_{H} $ for different values of the parameters Q, a, K, and ω in Fig. 5. For a black hole thermodynamic system, the positive and negative values of the black hole heat capacity correspond to local thermodynamic stability and instability, respectively. In Fig. 5, the left and right graphs correspond to the heat capacity graphs for different values of the (K, ω) parameters and Q parameter, respectively. Note that, when the black hole event horizon radius satisfies$ 0 < r_{H} < r_{H}^{C} $ , the system is thermodynamically stable, while for$ r_{H} > r_{H}^{C} $ , it is thermodynamically unstable. At$ r_{H}=r_{H}^{C} $ , the black hole heat capacity diverges, which indicates a second-order phase transition. -
In 1969, British physicist Penrose proposed a theoretical process for extracting energy from a rotating black hole, known as the Penrose process. Here, we calculate the Penrose process for a rotating and non-linear magnetic-charged black hole with an anisotropic matter field [49–51] and obtain the energy gain and extraction efficiency of this process.
The basic idea of the Penrose process is that a particle A moving on the equatorial plane (
$ \theta =\pi/2 $ ,$ u^\theta =0 $ ) of the black hole enters the ergosphere, where it decays into two massless particles, B and C. The massless particle B falls into a negative energy orbit, while massless particle C escapes to infinity. According to energy conservation, the energy of massless particle C exceeds the initial energy of particle A. The four-momentum of a particle is$ P_{\mu} =g_{\mu \nu} u^{\nu} , $
(55) where
$ u^{\nu} $ is the four-velocity of a particle. Considering that particles moving on the equatorial plane of a black hole satisfy energy and momentum conservation, we obtain$ P_{t} =-E=g_{t t} u^{t}+g_{t \phi} u^{\phi } , $
(56) $ P_{\phi}=L=g_{\phi \phi} u^{\phi }+g_{t \phi} u^{t}, $
(57) where E and L are the energy and angular momentum of the test particle, respectively. The radial equation for the geodesic motion of a test particle, derived from Eqs. (56) and (57), under the condition
$ P_{\nu} P^{\nu}=-\delta $ , where$ \delta=-1 $ ,$ 0 $ ,$ 1 $ , corresponding to time-like, null, and space-like geodesics, respectively, is as follows:$ \begin{aligned}[b] r^{4} u^{r}=\;&E^{2}\left(r^{4}+a^{4}+a^{2}\left(2 r^{2}-\Delta\right)\right)\\ &+L^{2}\left(a^{2}-\Delta\right)-4 a L E \rho r+\delta r^{2} \Delta. \end{aligned} $
(58) Assuming that the decay of particle A occurs at the turning point of the geodesic (
$ u^{r}=0 $ ), we can derive from Eq. (58) and the condition$ u^{r} = 0 $ that$ E=\frac{2 a \rho r L \pm r \sqrt{\Delta} \sqrt{L^{2} r^{2}-\delta\left(r^{4}+a^{4}+a^{2}\left(2 r^{2}-\Delta\right)\right)}}{r^{4}+a^{4}+a^{2}\left(2 r^{2}-\Delta\right)}, $
(59) and
$ L=\frac{-2 a \rho r E \pm r \sqrt{\Delta} \sqrt{E^{2} r^{2}+\delta\left(r^{2}-2 \rho r\right)}}{\Delta-a^{2}} . $
(60) To have negative energy, the angular momentum must satisfy
$ L < 0 $ . Therefore, from Eq. (60), we obtain$ 1-\frac{2 \rho}{r}<\delta \frac{\Delta}{L^{2}} . $
(61) Given that
$ g_{tt}(\theta=\dfrac{\pi}{2})=-\left(1 - \dfrac{2\rho}{r}\right) $ , it is evident from the inequality that this phenomenon can occur within the ergosphere. We consider a particle A ($ \delta=-1 $ ) with energy$ E_A=1>0 $ and angular momentum$ L=L_A $ entering the ergosphere. It then decays into two massless particles B and C ($ \delta=0 $ ), with energies and angular momenta ($ E_B $ ,$ L_B $ ) and ($ E_C $ ,$ L_C $ ), respectively. The massless particle B that falls into the event horizon has negative energy, while the massless particle C that escapes to infinity has positive energy. From Eq. (60), we obtain the relationship between the angular momentum and energy for A, B, and C as$ \begin{aligned}[b] L_{A} & =\frac{-2 a \rho r+r \sqrt{\Delta} \sqrt{2 \rho r}}{\Delta-a^{2}}, \\ L_{B} & =\frac{-2 a \rho r-r^{2} \sqrt{\Delta}}{\Delta-a^{2}} E_{B}, \\ L_{C} & =\frac{-2 a \rho r+r^{2} \sqrt{\Delta}}{\Delta-a^{2}} E_{C}. \end{aligned} $
(62) The conservation of energy and angular momentum is expressed as follows:
$ E_{B}+E_{C}=E_{A}=1, $
(63) and
$ L_{B}+L_{C}=L_{A}. $
(64) The energies of two massless particles B and C obtained from Eqs. (62), (63), and (64) are given by
$ \begin{aligned}[b] E_{B}&=\frac{1}{2}\left(1-\sqrt{1+\frac{a^{2}-\Delta}{r^{2}}}\right)\\ &=\frac{1}{2}\left(1-\sqrt{\frac{2 M r^{2}}{r^{3}+Q^{3}}+K r^{-2\omega }}\right), \end{aligned} $
(65) and
$ \begin{aligned}[b] E_{C}=&\frac{1}{2}\left(1+\sqrt{1+\frac{a^{2}-\Delta}{r^{2}}}\right)\\ =&\frac{1}{2}\left(1+\sqrt{\frac{2 M r^{2}}{r^{3}+Q^{3}}+K r^{-2\omega }}\right) . \end{aligned} $
(66) We conclude from Eq. (66) that the energy of the massless particle C escaping to infinity is higher than the energy of the particle A that initially entered the ergosphere of the black hole. The energy gain
$ \Delta E $ can be obtained as$ \Delta E=\frac{1}{2}\left(\sqrt{\frac{2 M r^{2}}{r^{3}+Q^{3}}+K r^{-2\omega }}-1\right)=-E_{B} . $
(67) According to Eq. (65), the maximum energy gain occurs at the event horizon, and the maximum efficiency of the Penrose process is then given by
$ \begin{aligned}[b] E_{ff \max }&=\frac{E_{A}+\Delta E}{E_{A}}\\ &=\frac{1}{2}\left(1+\sqrt{\frac{2 M r_{H}^{2}}{r_{H}^{3}+Q^{3}}+K r_{H}^{-2\omega }}\right). \end{aligned} $
(68) The behavior of the maximum efficiency of the Penrose process as a function of a for different values of the parameters K, ω, and Q is shown in Fig. 6. Note that
$ E_{f f \max } $ increases with the spin parameter a and magnetic charge Q. The influence of the anisotropic matter field on$ E_{f f \max } $ is realized through its effect on the event horizon$ r_{H} $ . -
In this study, we first obtained the exact solution of the static, spherically symmetric non-linear magnetic-charged black hole with an anisotropic matter field and further calculated the curvature of the black hole. The results show that the curvature diverges at
$ r=0 $ , indicating that a non-linear magnetic-charged black hole with an anisotropic matter field has a singularity. We used a modified Newman-Janis algorithm to obtain the solution for a rotating and non-linear magnetic-charged black hole with an anisotropic matter field. Subsequently, we studied the energy-momentum tensor and weak energy condition for the rotating black hole, finding that this black hole violates the weak energy condition.We investigated the geometric and thermodynamic properties of a rotating and nonlinear magnetic-charged black hole with an anisotropic matter field. Regarding the geometric properties, we first analyzed the influence of the spin parameter a and magnetic charge Q on Δ in Fig. 3. We found that the influence of Q mainly occurs within a small interval. Further analysis of Fig. 3 and Δ (
$ Q = 0 $ ) [21] revealed that the rotating black hole may have two horizons for$ \dfrac{1}{2} < \omega < 1 $ , regardless of whether K is positive or negative. For$ \omega > 1 $ and$ K < 0 $ , Δ can have at most two roots, while for$ \omega > 1 $ and$ K > 0 $ , the black hole can have at most three roots, indicating the possible existence of three horizons. Further analysis of the influence of various parameters on the ergosphere of the black hole showed that the area of the ergosphere increases with the spin parameter a and magnetic charge Q.Regarding thermodynamic properties, we obtained the entropy, Hawking temperature, and angular velocity of the rotating and non-linear magnetic-charged black hole with an anisotropic matter field. To further derive other thermodynamic quantities, we presented the first law of thermodynamics for the black hole. Using
$ \Delta= 0 $ , we derived the squared mass formula and calculated other thermodynamic quantities. We further studied the stability of this black hole through its heat capacity and found that the black hole is thermodynamically stable when the event horizon radius satisfies$ 0 < r_{H} < r_{H}^{C} $ . However, when$ r_{H} > r_{H}^{C} $ , the black hole becomes thermodynamically unstable. At$ r_{H} = r_{H}^{C} $ , the heat capacity diverges, which indicates a second-order phase transition.Finally, we calculated the energy extraction efficiency of the Penrose process for a rotating and non-linear magnetic-charged black hole with an anisotropic matter field and found that the maximum efficiency increases with the spin parameter a and magnetic charge Q. The influence of the anisotropic matter field on the maximum efficiency of the Penrose process is realized through its impact on the event horizon.
In this study, we primarily analyzed black hole properties across various parameter values without imposing constraints on these parameters. The shape of the black hole shadow is sensitive to parameter changes. By combining observational data from the Event Horizon Telescope (EHT) and others, we can effectively eliminate parameter regions that do not match the measurements. In future research, we will constrain the parameter space of this black hole using observational features of the black hole shadow, such as the axis ratio and angular diameter.
Rotating and non-linear magnetic-charged black hole with an anisotropic matter field
- Received Date: 2025-04-29
- Available Online: 2025-10-15
Abstract: We present the solution for a non-linear magnetic-charged black hole with an anisotropic matter field and extend it to obtain the corresponding rotating black hole solution using the modified Newman-Janis algorithm. The event horizon and ergosphere of the rotating black hole are analyzed from a geometric perspective, revealing that the rotating black hole can have up to three horizons. The first law of thermodynamics and the squared-mass formula for the rotating black hole are derived from a thermodynamic perspective, based on which we obtain thermodynamic quantities and examine the thermodynamic stability of the rotating black hole. Additionally, we calculate the Penrose process for the rotating black hole, showing the influence of various black hole parameters on the maximum efficiency of the Penrose process.