Thermodynamic Topology of Kiselev-AdS Black Holes within f (R, T) gravity

Black hole thermodynamics is a fascinating field that explores the intersection of thermodynamics, quantum mechanics, and general relativity. To obtain insights into quantum gravity, researchers investigate various approaches, including the study of black holes and their thermodynamic properties, which serve as a promising tool [1, 2]. This approach combines classical thermodynamics with statistical mechanics to interpret the entropy of black holes as the number of quantum degrees of freedom near the event horizon. The thermodynamics of black holes have been explored within different frameworks, with results documented in [3−19]. A key concept in this field is the Hawking-Page phase transition [20], which highlights a specific relationship between gravity and thermodynamics. This concept suggests a phase transition between thermal radiation in AdS space and black holes, known as confinement/deconfinement phase transitions in the context of the AdS/CFT correspondence. Additionally, in the thermodynamics of black holes, the study of phase transitions and critical points consistently garners significant attention[21−25]. Notably, the structural similarities observed between the liquid-to-gas phase transition in van der Waals fluids and the first-order phase transition between small and large black holes are particularly noteworthy.

Recently, innovative methods have emerged to analyze and compute critical points and phase transitions in black hole thermodynamics. One notable approach is the topological method, as detailed in [26−29]. To adopt a topological viewpoint in thermodynamics, Duan’s topological current ϕ mapping theory is highly recommended. SW Wei et al. proposed two distinct methods for studying topological thermodynamics, focusing on temperature and generalized free energy functions[26, 27]. The temperature method involves analyzing the temperature function by removing pressure and utilizing the auxiliary and topological parameter $ 1/\sin \theta $. Based on these assumptions, the potential is then constructed. The generalized free energy function method treats black holes as defects in the thermodynamic parameter space. Their solutions can be explored using the generalized off-shell free energy. In this framework, the stability and instability of the black hole solutions are indicated by positive and negative winding numbers, respectively. For more study, see [26−56].

Researchers, by utilizing black hole actions and extracting their effective potential, have discovered that the maximum and minimum of the potential function, which indicate changes in the field structure, can be examined by transforming the potential into a vector field and mapping it onto a two-dimensional space in ($ r, \Theta $) coordinates. In this scenario, these maximum and minimum, similar to winding numbers, manifest as rotations around zero points or poles. Since the mathematical structure of this method was designed for scalars, SW Wei extended it to thermodynamics and used it to study temperature, meaning that the coordinates of the location of temperature changes could also play the role of the pole or zero point for a vector field made of temperature[26, 27, 57−60]. In the T method, the goal is to obtain critical points and determine their position using the topological charge of the system. In this method, we do not know much about the phase transition and how to do it, which makes it a simpler method. Then SW Wei and his colleagues extended this method to a more comprehensive potential function, namely the internal energy or, more specifically, the Helmholtz free energy. In the F method, we use the Helmholtz free energy the structure of the τ function, and the topological charge to determine which phase transition the structure undergoes. If the process has certain maxima and minima in a radial interval, then the system will undergo a first-order phase transition with the creation of an intermediate black hole, and we will observe the process of creation (maximum) and annihilation (minimum) during the process. However, if the function does not have a maximum or minimum and undergoes a uniform temperature trend during the radial changes, the small black hole will directly and gradually turn into a large black hole, and we will witness a second-order phase transition. Both scenarios can be identified using the F method[26, 27, 57−60]. In other words, the F method and the T method are two distinct approaches used to study the thermodynamic topology of black holes. The generalized free energy method involves analyzing the free energy landscape of black holes. By examining the free energy, researchers can identify critical points and phase transitions. The stability of black hole solutions is determined by the sign of the second derivative of the free energy. Positive values indicate stability, while negative values indicate instability. The method uses the concept of topological charge, which is related to the winding numbers of the free energy function. These winding numbers help classify black holes based on their topological properties[26, 27, 57−60]. T Method focuses on the temperature function of black holes. It involves studying the temperature function by removing pressure and introducing an auxiliary topological parameter. Similar to the F method, the T method also considers topological aspects, but it does so through the lens of temperature variations and their impact on the system. Both methods provide a topological perspective on black hole thermodynamics, allowing researchers to classify black holes based on their topological properties. Both methods aim to determine the stability and instability of black hole solutions, although they use different functions (free energy vs. temperature) to achieve this. While the F method focuses on free energy and the T method on temperature, they are complementary approaches that together offer a more comprehensive understanding of black hole thermodynamics. By using both the F method and the T method, researchers can gain a deeper and more nuanced understanding of the thermodynamic properties and stability of black holes. Each method provides unique insights that, when combined, enhance our overall knowledge of black hole thermodynamics[26, 27, 57−60].

The concept of the topological photon sphere in black holes is a fascinating area of study in theoretical physics. The photon sphere is a region where gravity is so strong that photons (light particles) are forced to travel in orbits. This sphere is crucial for understanding the properties of black holes, including their shadow and the bending of light around them[61, 62]. In the context of topology, researchers explore the properties and behaviors of these photon spheres using topological methods. These methods can provide insights into the stability and structure of the photon sphere, as well as its interactions with other physical phenomena in the black hole's vicinity.

The $ f(R, T) $ gravity model presents an intriguing extension of Einstein's equations of motion by incorporating terms related to the Ricci scalar and the trace of the energy-momentum tensor. This theoretical framework has recently attracted considerable attention due to its potential to explain various cosmological phenomena. Additionally, it offers an alternative gravity theory that can account for dark matter and dark energy. Numerous studies have been conducted in the context of $ f(R, T) $ gravity, as noted in the literature[63−75]. According to general relativity (GR), black holes possess gravitational forces so intense that nothing, not even light, can escape their pull. It is fascinating to consider that the effects of $ f(R, T) $ gravity might be more significant for compact objects like black holes. Therefore, exploring black holes within the $ f(R, T) $ gravity framework is particularly compelling. Recognizing this, Santos et al. formulated a black hole solution similar to Kiselev's within the $ f(R, T) $ gravity framework[76]. Subsequently, the rotating version of this black hole solution has been further elaborated upon[77].

Kiselev-AdS black holes within the $ f(R, T) $ gravity framework are an intriguing area of study in theoretical physics. These black holes are solutions to the modified Einstein field equations, incorporating both the Ricci scalar R and the trace of the energy-momentum tensor T. This framework allows for the exploration of various cosmological and astrophysical phenomena, including the effects of dark energy and dark matter[77−79]. Recent research has focused on constructing rotating Kiselev black holes within this framework using the revised Newman-Janis algorithm. These solutions include special cases like the Kerr and Kerr-Newman black holes and exhibit unique spacetime structures influenced by the parameters of the $ f(R, T) $ function. Additionally, studies have shown that Kiselev-AdS black holes in $ f(R, T) $ gravity can exhibit phase transitions similar to those seen in thermodynamic systems, such as van der Waals-like and Hawking-Page-like transitions. These phase transitions are influenced by the state parameter ω and the $ f(R, T) $ gravity parameter γ, providing insights into the nature of dark energy and the thermodynamic properties of black holes[77−79]

These concepts motivated us to study the Kiselev-AdS black holes from $ f(R, T) $ gravity based on thermodynamic topology. We organize the paper as follows: In Section 2, we briefly introduce the mentioned black hole. Also, we study the curvature singularities and energy conditions of Kiselev-AdS black holes from $ f(R, T) $ gravity in Section 3. We investigate the thermodynamic topology from the perspective of topology using the generalized off-shell Helmholtz free energy and temperature methods and study the photon sphere of Kiselev-AdS black holes from $ f(R, T) $ gravity in Section 4. Finally, we describe the results of our work in detail in Section 5.

In this study, we aim to obtain a static spherically symmetric AdS black hole solution as proposed by Kiselev within the framework of $ f(R, T) $ gravity. In our analysis, we consider the specific form $ f(R, T) = R + 2 f(T) $[78, 104]. One can assume a static geometry in a spherically symmetric space-time, described by the following line element,

With respect to [78] and using the form $ f(T) = \gamma T $, the unknown functions is obtained as follows,

Here, γ is the $ f(R, T) $ gravity parameter, M is the mass of the black hole, and k is a constant of integration. By setting $\gamma \rightarrow 0$, the well-known Kiselev-AdS black hole solution is recovered as follows in General Relativity [105],

The black hole mass can be determined by solving $ B(r_+) = 0 $ as follows,

The Hawking temperature defined as $ T = \kappa / 2\pi $. So, it can be expressed as follows,

We use the well-known Hawking-Bekenstein area law to express the entropy $ S_+ $ in terms of the event horizon radius $ r_+ $ as follows:

To provide a suitable proof of the singularity and uniqueness of the black hole solution, we have to perform some analyses in terms of scalar invariants such as the Ricci scalar, the Ricci squared scalar, and the Kretschmann scalar. The Ricci scalar with respect to the relevant metric can be given as,

The Ricci squared function, on the other hand, can be expressed in relation to the metric function as,

The Kretschmann scalar can be represented as,

Upon deeper inspection of the expressions $ (7) $, $ (8) $ and $ (9) $, it can be seen that the black hole metric is singular for all admissible values of the parameters ω and γ. Indeed, the existence of the singularity results from the mass term and the relevant term of the constant k in the black hole metric. By constraining $ \omega <-1/3 $ and $ \gamma<0 $, the singularity will disappear due to these imposing conditions. Even so, to rule out the singularity made from the mass and charge terms, it may be worth outlining a routine with a non-linear charge distribution function analogous to Ref[80]. Along this work, we will not consider such a situation and stick to the metric function (37) for the rest of the analysis. A few remarks about the scalars mentioned show that the Ricci scalar and the Ricci squared are not a function of the black hole mass M. On the other hand, the Kretschmann scalars are functions of the black hole mass M, so any variation in the black hole mass, Kiselev $ f (R, T) $ parameters can bring about substantial variations in these same scalars. Hence, the scalars under consideration show that the black hole solution is indeed unique and that the AdS background, as well as the Kislev structure, significantly modify the black hole spacetime. So, an analysis of scalar invariant limits at $ r = 0 $ might provide as follows,

On the other hand, analyzing at a large distance, the relevant behavior of the scalar invariants can be provided as follows,

This in turn shows that the Ricci scalar, the Ricci squared scalar, and the Kretschmann scalar comprise a finite term at large distances. In a nutshell, these scalars are proof that our black hole solution is completely consistent and that the parameters of Kiselev, as well as the $ f(R, T) $ parameter theory, do not affect the asymptotic behavior of spacetime. The above analysis can be demonstrated based on the three scalar invariants graph. For that reason, Figure 1 displays the radial variations of scalar invariants in spacetime (1), taking into account the Ricci scalar (Fig. 1(a)), the square of the Ricci tensor (Fig. 1(b)), and the Kretschmann scalar (Fig. 1(c)) for various values of the parameter γ with $ \omega=1/3 $. The figure reveals that the scalar invariants show a declining trend along the radial coordinates. Furthermore, as the value of ω and γ decrease, the Ricci squared scalar and the Kretschmann scalar increase. On the other hand, a proportional relation is observed between the variation of the Ricci scalar and the parameter γ. Similarly, it is found that all three scalars drop to zero when r is outside the event horizon. We now look at the weak, strong, and dominant energy conditions for the source fluid. The energy-momentum tensor for the anisotropic source fluid is given by [81−84]

Figure 1.  (color online) Ricci, Ricci squared and Kretschmann scalars representation for various values of black hole system.

where ρ is the energy density, $ P_i\quad (i = 1, 2, 3) $ are the pressures for the source fluid, and $ e_i^\mu $ represents the vielbein components. Eq. 11 refers to the source of the Einstein field equations that might be needed to derive a black hole solution. We consider the static and spherically symmetric line element proposed in Eq. 1, and we require that this metric satisfies the Einstein equations according to the auxiliary term in Eq. 2. To be more specific, because we are looking for an Anti-de Sitter-type contribution with a negative cosmological constant, we obtain that the Einstein field equations give[81−84]

Notice that we employ the convention $ G = 1/(8\pi) $, whereas [81] uses $ G = 1 $, G being Newton's constant. These results are broadly in agreement with the ones given in [81, 81−84], where $ P_r $ is the longitudinal pressure, while $ P_\theta $ and $ P_\phi $ are the transverse pressures. Moreover, the term $ 8\pi P_\Lambda $, which is included in equations (12) and (13), is essential to ensure the presence of the negative cosmological constant in the Einstein field equations and according to current studies [85, 86], one can interpret it as a thermodynamic pressure, such that $ P_\Lambda=-\Lambda/8\pi= 3/8\pi \ell^2 $. We will now examine the energy conditions in question for the black hole solution [87−90].

Figure 10.  (color online) (a): The topological potential H(r) for Kiselev-AdS black holes in f(R, T) gravity, (b): The normal vector field n in the $ (r-\theta) $ plane. The photon sphere are located at $ (r,\theta)=(1.615119922842,1.57) $ with respect to $ (\gamma=-2.68, m=1,l=1,k=1) $

Figure 11.  (color online) (a): The topological potential H(r) for Kiselev-AdS Black Holes, (b): The normal vector field n in the $ (r-\theta) $ plane. The photon sphere are located at $ (r,\theta)=(0.573037724709,1.57),(r,\theta)=(1.658519337581,1.57) $ with respect to $ (\gamma=-3, m=1,l=1,k=1) $

Figure 13.  (color online) (a): Metric function with different γ for Kiselev-AdS black holes in f(R, T) gravity, (b): The normal vector field n in the $ (r-\theta) $ plane. The photon sphere are located at $ (r,\theta)=(1.591046,1.57) $ with respect to $ (\gamma=-6.605, m=1,l=1,k=1) $. ??: The normal vector field n in the $ (r-\theta) $ plane. The photon sphere are located at $ (r,\theta)=(1.107899985366,1.57),(r,\theta)=(2.282979884819,1.57) $ with respect to $ (\gamma=-6.45, m=1,l=1,k=1) $

The elements of the stress-energy tensor $ T^{\mu\nu} $ for Kiselev-AdS black hole within $ f(R, T) $ gravity are as follows,

● For each time vector $ t^\mu $, the weak energy condition (WEC) demands that $ T_{\mu\nu}\, t^\mu t^\nu\geqslant0 $ everywhere. This is equivalent to [91]

and so $ \rho+P_r=0 $ and

Given that $ \omega>1/3 $ and that $ \gamma>0 $, the WEC is satisfied. It is evident from Fig. 2a that, given a span of horizon radii, the energy density ρ is positive in accordance with the value of fixed ω and for a spectral value of the parameter γ.

Figure 2.  (color online) Energy conditions using $ \omega=1/3 $ and $ k=1 $.

● For each null vector $ t^\mu $ in the entire spacetime, the null energy condition (NEC) requires that $ T_{\mu\nu}\, t^\mu t^\nu\geqslant0 $. When $ \omega>1/3 $, the NEC predicts $ \rho+P_r \geqslant 0 $, which is identically zero, and $ \rho+P_\theta=\rho+P_\phi \geqslant 0 $, which is satisfied for equation (17).

● The strong energy condition (SEC) states that, for every time vector $ t^\mu $, $ T_{\mu\nu}\, t^\mu t^\nu\geqslant 1/2 \,T_{\mu\nu} t^\nu t_\nu $ globally, assuming that [91]

Upon careful inspection, the following restriction is the only one that satisfies the SEC:

Graphical analysis shows that energy conditions such as WEC, NEC, and SEC vary as a function of the variable r (see Fig. 2). Thus, we observe that for each choice of values of ω with a range of values of the parameter γ, it is observed that the WEC is satisfied for the considered case such as $ \omega=1/3 $ with $ 16\leqslant\gamma\leqslant84.5 $ (Fig. 2(a)). Similarly, the NEC seems to be satisfied also for the case such that $ \omega=1/3 $ with $ 16\leqslant\gamma\leqslant84.5 $ (Fig. 2(b)). The Kiselev predominantly satisfies the SEC, as shown in Fig. 2(c). This is akin in contrast to the quintessence of dark energy within GR. Indeed, a violation of the SEC is understood as a violation of the attractive pattern of gravity, as exemplified by dark energy accelerating the expansion of the universe in cosmological studies, together with the matter content of the background of a regular black hole, whose singularity has been superseded by a Sitter core. For brevity, it is worth noting that the general mapping between SEC and gravity behavior is the focus of an ongoing search in the literature. It is well known that the validity of the SEC in a gravity environment is typically related to the attractive aspect of gravity. This is true in GR, where the SEC must be assumed to secure the attractive nature of the theory - focus theorem [92, 93]. Nevertheless, such a connection need not have entirely general validity in extended gravity. A paradigmatic illustration is provided by $ f(R) $ gravity, shown in [94], involving the Raychaudhuri equation. It has already been revealed that even supposing the standard SEC, the Raychaudhuri equation can present positive inputs from spacetime geometry, generally considered a possible sign of repulsive gravity. As a result, the SEC/attractive gravity paradigm appears invalid in this context. A piece of stringent and comprehensive evidence that the attractiveness of gravity is no longer guaranteed by the SEC in extended gravity is beyond the present analysis. This task will be left to the future.

The photon sphere is a region surrounding a black hole where the gravitational pull is so intense that it forces photons into circular orbits. This sphere is crucial because it defines the boundary within which light can be trapped in orbit around the black hole, leading to intriguing phenomena such as gravitational lensing and the black hole's shadow[95−103]. Thermodynamic topology examines the topological characteristics of black holes about their thermodynamic behavior. This method aids in understanding phase transitions and the stability of black holes. By studying the topological charge and critical points, researchers can categorize different phases of black holes and forecast their behavior under various conditions. For example, in hyperscaling-violating black holes, the topological charge can reveal whether the black hole is stable or if it will transition to a different phase. These concepts are essential for comprehending the complex nature of black holes and their interactions with light and matter. In the following section, we will discuss these concepts in the context of Kiselev-AdS black holes within $f(R, T)$ gravity.

To investigate the thermodynamic properties of black holes, various quantities are utilized. For instance, mass and temperature can describe the generalized free energy. Considering the relationship between mass and energy in black holes, we can represent our generalized free energy function as a standard thermodynamic function in the following form [27].

where τ denotes the Euclidean time period, while T (the inverse of τ) represents the temperature of the ensemble. The generalized free energy is on-shell only when $ \tau = \tau_{H} = \dfrac{1}{T_{H}} $. As stated in [27], a vector ϕ is constructed as follows,

Where $ \phi^{\Theta} $ diverges, the vector direction points outward at $ \Theta = 0 $ and $ \Theta = \pi $. The ranges for $ r_{H} $ and Θ are $ 0 \leq r_{H} \leq \infty $ and $ 0 \leq \Theta \leq \pi $, respectively. Using Duan's ϕ-mapping topological current theory, a topological current can be defined as follows,

Given $ n = (n^1, n^2) $, where $ n^1 = \dfrac{\phi^r}{\|\phi\|} $ and $ n^2 = \dfrac{\phi^\Theta}{\|\phi\|} $, Noether's theorem ensures that the resulting topological currents are conserved.

To determine the topological number, we reformulate the topological current [27],

The Jacobi tensor determine as,

The Jacobi vector reduces to the standard Jacobi when $ \mu=0 $, as demonstrated by $ J^{0}\left(\dfrac{\phi}{x}\right)=\dfrac{\partial(\phi^1,\phi^2)}{\partial(x^1,x^2)} $. Eq.(23) shows that $ j^{\mu} $ is only non-zero when $ \phi=0 $. Through some calculations, we can express the topological number or total charge W as follows,

Here, $ \beta_i $ denotes the positive Hopf index, which counts the loops of the vector $ \phi^a $ in the ϕ space when $ x^\mu $ is near the zero point $ z_i $. Meanwhile, $ \eta_i=\text{sign}(j^0(\phi/x)_{z_i})=\pm 1 $. The quantity $ \widetilde{\omega}_i $ represents the winding number for the i-th zero point of ϕ in Σ. Note that the winding number is independent of the shape of the region where the calculation occurs. The value of the winding number is directly related to black hole stability, with a positive (negative) winding number corresponding to a stable (unstable) black hole state. Using Eqs. (4), (6), and (21), we derive the generalized Helmholtz free energy for Kiselev-AdS black holes within the framework of $ f(R, T) $ gravity.

Based on the discussion in the previous section, the form of the function $ (\phi^{r}, \phi^\Theta) $ is determined as follows,

The unit vectors $ \mathbf{n}_1 $ and $ \mathbf{n}_2 $ are computed using Eq.(28). Next, we find the zero points of the $ \phi^{r} $ component by solving $ \phi^{r} = 0 $ and derive an expression for τ as follows,

We face a single zero point in Figs. 3(b), 4(b), 4(d), 5(b), 6(b), indicating one topological charge determined by the free parameters mentioned in the study. This charge corresponds to the winding number and is located within the blue contour loops at coordinates $(r, \theta)$. The sequence of the illustrations is determined by the parameters ω and γ; for example, in Fig. (3), we have ($\omega = 0$ with $\gamma = -1.9, -0.4$), and for Figs. (4) through (6), ω takes the values (ω=$\dfrac{1}{3}$ with $\gamma = 16, 84.5$), (ω=$-\dfrac{2}{3}$ with $\gamma = -5.7, -5.61$), and (ω=$-\dfrac{4}{3}$ with $\gamma = -6.63, -6.685$), respectively. In all Figs, we set $k = \ell = 1$. The findings from all Figs show that the distinctive feature of total topological charge $ +1 $ is the zero point enclosed within the contour.

Figure 3.  (color online) The curve corresponding to Eq.(30) is depicted in Figs. 3(a)and 3(c). In Figs. 3(b) and 3(d), the colorful arrows illustrate the vector field (n) on a segment of the $ (r-\theta) $ plane for Kiselev-AdS black holes within f(R, T) gravity, with parameters $ (\ell=k=1, \omega=0, \gamma=-1.9) $ and $ (\ell=k=1, \omega=0, \gamma=-0.4) $, respectively. The ZPs are positioned at $ (r, \theta) $ on the circular loops.

Figure 4.  (color online) The curve corresponding to Eq.(30) is depicted in Figs. 4(a) and 4(c). In Figs. 4(b) and 4(d), the ZPs are located at $ (r, \theta) $ on the circular loops with parameters $ (\ell=k=1, \omega=\frac{1}{3}, \gamma=16) $ and $ (\ell=k=1, \omega=\frac{1}{3}, \gamma=-84.5) $.

Figure 5.  (color online) The curve corresponding to Eq.(30) is depicted in Figs. 5(a) and 5(c). In Figs. 5(b) and 5(d), the ZPs are located at $ (r, \theta) $ on the circular loops with parameters $ (\ell=k=1, \omega=-\frac{2}{3}, \gamma=-5.7) $ and $ (\ell=k=1, \omega=-\frac{2}{3}, \gamma=-5.61) $.

Figure 6.  (color online) The curve corresponding to Eq.(30) is depicted in Figs. 6(a) and 6(c). In Figs. 6(b) and 6(d), the ZPs are located at $ (r, \theta) $ on the circular loops with parameters $ (\ell=k=1, \omega=-\frac{4}{3}, \gamma=-6.63) $ and $ (\ell=k=1, \omega=-\frac{4}{3}, \gamma=-6.685) $.

Additionally, our findings indicate that by increasing the parameter γ for $(\omega = 0, -\dfrac{2}{3})$, the number of topological charges increases, as shown in Figs. 3(b), 3(d), 5(b), 5(d), but the total topological charge remains $W = +1$. Conversely, for $(\omega=-\dfrac{4}{3})$, the number of topological charges increases by decreasing the value of the parameter γ. Generally, for $(\dfrac{1}{3})$, changing the parameter γ does not alter the number of topological charges. However the the parameters γ and ω affect the number of topological charges in the entire system, while the total topological charges for all modes are equal $W = +1$. Also, Figs. 3(d), 5(d), 6(d) depict three topological charges $(\widetilde{\omega} = +1, -1, +1)$, resulting in a total topological charge of $W = +1$. In Figs. 3(a), 3(c), 4(a), 4(c), 5(a), 5(c), 6(a), 6(c), we plotted the trajectory corresponding to Eq.(30) across various free parameter values. The results are summarized in Table 1.

Free parameters	$\tilde{\omega}$	W
$(\ell=k=1, \omega=0, \gamma=-1.9)$	$+1$	$+1$
$(\ell=k=1, \omega=0, \gamma=-0.4)$	$+1,-1,+1$	$+1$
$(\ell=k=1, \omega=\frac{1}{3}, \gamma=16)$	$+1$	$+1$
$(\ell=k=1, \omega=\frac{1}{3}, \gamma=84.5)$	$+1$	$+1$
$(\ell=k=1, \omega=-\frac{2}{3}, \gamma=-5.7)$	$+1$	$+1$
$(\ell=k=1, \omega=-\frac{2}{3}, \gamma=-5.61)$	$+1,-1,+1$	$+1$
$(\ell=k=1, \omega=-\frac{4}{3}, \gamma=-6.63)$	$+1$	$+1$
$(\ell=k=1, \omega=-\frac{4}{3}, \gamma=-6.685)$	$+1,-1,+1$	$+1$

Table 1.  Summary of the results for F-Model

In the extended thermodynamic framework, the temperature T of a system like a black hole is expressed as a function of pressure P, and an additional parameter ξ. Critical points are identified using the conditions,

Recent research indicates that each critical point can have a topological charge, which can be conventional or novel. To eliminate the parameter P, the relation:$ \left(\dfrac{\partial T}{\partial r_{H}}\right)_{P, \xi} = 0 $ is used. To study the topological charge, the thermodynamic function Φ is defined as,

where $\dfrac{1}{\sin \theta}$ simplifies the calculations. A new vector field $\phi = (\phi^{r_H}, \phi^\theta)$ is defined using Duan’s ϕ-mapping theory,

The vector field ϕ is zero at $\theta = \dfrac{\pi}{2}$, helping to identify critical points. The points $\theta = 0$ and $\theta = \pi$ serve as boundaries in the parameter space. When the vector field $\phi^a$ is zero, its topological current $j^\mu$ becomes non-zero. Conversely, contour C can be parameterized by $ \vartheta \in (0, 2\pi) $ with $ r = a \cos \vartheta + r_0 $ and $ \theta = b \sin \vartheta + \dfrac{\pi}{2} $, where $ (r_0, \dfrac{\pi}{2}) $ represents the center of the contour. We introduce a new quantity to measure the deflection of the vector field along the contour,

where $ i, j = S, \theta $. Consequently, the topological charge $ Q_t $ is given by,

The temperature is reformulated by eliminating the pressure parameter,

Subsequently, by substituting equation (35) into equation (31), the thermodynamic function Φ for the black hole can be derived as follows,

Topological charges can be identified due to the presence of critical points within acceptable regions for free parameters. By applying Eqs. (30), (31), (32), (35), and (36) for $\omega = 0$ and $\omega = -4/3$ with various values of γ as shown in Figs. 7(a), 7(c), 8(a) and 8(c), the topological charges of each critical point can be determined. Consequently, from Figs. (7) and (8), we find that $Q_{CP1} = -1$. Additionally, we conclude that $Q_{total} = -1$.

Figure 7.  (color online) The arrows represent the vector field n on the $ r_H-\theta $ plane for the Kiselev-AdS black holes within f (R, T) gravity with $ (k =\ell= 1, \omega = 0, \gamma = -0.4) $ and $ (k =\ell= 1, \omega = 0, \gamma = -1.9) $ in Figs. 7(a) and 7(c), respectively. Isobaric and spinodal curves (blue lines) for the Kiselev-AdS Black Holes within f (R, T) gravity in Figs. 7(b) and 7(d)

Figure 8.  (color online) The arrows represent the vector field n on the $ r_H-\theta $ plane for the mentioned black hole with $ (k =\ell= 1, \omega = -4/3, \gamma = -6.63) $ and $ (k =\ell= 1, \omega = -4/3, \gamma = -6.685) $ in Figs. 8(a) and 8(c), respectively. Isobaric and spinodal curves (blue lines) in Figs. 8(b) and 8(d)

As illustrated in Figs. 7(b), 7(d), 8(b) and 8(d), the conventional critical points represent the local minima and maxima of the spinodal curve. The range of various parameters for the black hole in question is constrained by the conditions of the event horizon, which has been thoroughly examined in [78]. Therefore, for other values of the parameter ω, specifically $1/3$ and $-2/3$, the topological charges were not obtained within this range. Consequently, the discussed topics apply only to (dust) and (phantom) fields.

Referring to [106−109], we start with a regular potential,

By examining this potential, we can determine the radius of the photon sphere, which is located at,

We then introduce a vector field $ \phi = (\phi^r, \phi^\theta) $, defined as follows,

Consequently, the total charge is given by,

In conclusion, the presence of a zero point within a closed curve indicates that the charge Q is precisely equal to the winding number. For more study please see [108]. Now with respect to B(r) in Eq. 2 the above functions will be as follows,

Since in this article, our main goal is to investigate the effect of gravitational corrections on the model, for this aim, we will study the structure of the photon sphere and its parameter range in three parts. It should be noted that for the sake of simplicity, we consider the values of $ M=1, \ell=1, k=1 $ for all cases in this section.

We have drawn the metric function for $ \omega=0 $ and different γ. As can be seen in Fig. 9, the metric function for $ \gamma>-2.6911 $ always has a root. For this case, the general form of the Eqs.(37) and (39) will be as follows,

Figure 9.  (color online) Metric function with different γ for Kiselev-AdS black holes in $f(R, T) $ gravity

Based on Table 2, as evident from Figs. 9, 10 and 11, in the first region, the potential structure possesses only a single local maximum or an unstable photon sphere with a negative unit charge. However, in the second region, with the disappearance of the event horizon, a local minimum (stable photon sphere) emerges beyond the horizon. This results in the total topological charge becoming zero, rendering the space in the form of a naked singularity, which is also confirmed by the metric function.

Kiselev-AdS Black Holes	Fix parametes	Conditions	TTC	$(R_{PLPS})$
unstable photon sphere	$k=1,m=1,l=1$	$-2.6911\geq\gamma $	$-1$	$1.616399193583$
naked singularity	$k=1,m=1,l=1$	$-8.37<\gamma< -2.6911$	$ 0 $	$-$

Table 2.  $R_{PLPS}$: The minimum or maximum possible radius for the appearance of an unstable photon sphere. TTC: Total Topological Charge

An interesting aspect of this model, in comparison to previous works, is the elimination of the forbidden region. In this model, beyond the singularity region, contrary to previous models where a forbidden zone existed, the metric function exhibits a single root at $ \gamma = -8.38 $, and a photon sphere appears. However, whether these γ values are permissible within Einstein's equations and whether they satisfy the necessary energy conditions requires further scrutiny.

As can be seen in Fig. 12a, with respect to $ \omega=-2/3 $ the metric function for $ \gamma>-5.7903 $ always has a root. Now with respect to Eqs. (37) and (38) we have,

Figure 12.  (color online) (a): Metric function with different γ for Kiselev-AdS black holes in f(R, T) gravity, (b): The normal vector field n in the $ (r-\theta) $ plane. The photon sphere are located at $ (r,\theta)=(1.53024,1.57) $ with respect to $ (\gamma=-5.68, m=1,l=1,k=1) $. (c): The normal vector field n in the $ (r-\theta) $ plane. The photon sphere are located at $ (r,\theta)=(1.164619828617,1.57),(r,\theta)=(2.495766445262,1.57) $ with respect to $ (\gamma=-6.1, m=1,l=1,k=1) $

The behavior of the model based on the ranges obtained for the γ parameter is fully displayed in Table 3 we have drawn a sample from each area in Fig. 12.

Kiselev-AdS Black Holes	Fix parametes	Conditions	TTC	$(R_{PLPS})$
unstable photon sphere	$k=1,m=1,l=1$	$-5.7902\geq\gamma $	$-1$	$1.616366040729$
naked singularity	$k=1,m=1,l=1$	$-6.8544<\gamma< -5.7902$	$ 0 $	$-$

Table 3.  $R_{PLPS}$: The minimum or maximum possible radius for the appearance of an unstable photon sphere. TTC: Total Topological Charge

In this scenario, the impact of corrections operates in the exact opposite manner compared to previous cases. Unlike the previous instances where corrections influenced the system by increasing the γ, in this case, it is the reduction of γ that has a more significant effect. With respect to $ \omega=-4/3 $ the metric function for $ \gamma<-6.5935 $ always has a root. For this case, the general form of the equations due to Eqs. (37) and (39) will be as follows,

According to the ranges obtained for the γ parameter, the behavior of the model and the total topological charge can be seen in the Table 4, which in Fig. 13 we have drawn a sample from each area.

Kiselev-AdS Black Holes	Fix parametes	Conditions	TTC	$(R_{PLPS})$
unstable photon sphere	$k=1,m=1,l=1$	$-6.5936\leq\gamma $	$-1$	$1.593135505582$
naked singularity	$k=1,m=1,l=1$	$-6.5936<\gamma< -5.8$	$ 0 $	$-$

Table 4.  $R_{PLPS}$: The minimum or maximum possible radius for the appearance of an unstable photon sphere. TTC: Total Topological Charge

Before concluding this section, it is important to note that the studies and calculations for the radiation region ($ \omega = 1/3 $) did not present any new insights that would lead to different conclusions compared to the three plotted regions, except for maintaining the trend of changes in the γ parameter range. Therefore, to avoid redundancy, this region has been omitted.

In this paper, we have thoroughly investigated the topological charge and the conditions for the existence of the photon sphere (PS) in Kiselev-AdS black holes within $f(R, T)$ gravity. We established their topological classifications using two distinct methods based on Duan’s topological current ϕ-mapping theory viz temperature and the generalized Helmholtz free energy method. By analyzing the critical and zero points (topological charges and topological numbers) for different parameters, we revealed that the Kiselev parameter ω and the $f(R, T)$ gravity parameter γ significantly influence the number of topological charges of black holes, providing novel insights into their topological classifications. Our findings indicate that for given values of the free parameters, there exist total topological charges ($Q_{total} = -1$) for the T-method and total topological numbers ($W = +1$) for the generalized Helmholtz free energy method. Notably, in contrast to the scenario where $\omega = 1/3$, increasing the parameter γ in other cases results in an increase in the number of total topological charges for the black hole. Interestingly, for the phantom field ($\omega = -4/3$), decreasing the parameter γ leads to an increase in the number of topological charges.

Additionally, our study of the photon sphere reveals that the simultaneous presence of γ and ω effectively expands the permissible range for γ, allowing the model to exhibit black hole behavior over a larger domain. It is evident that with the stepwise reduction of ω, the region covered by singularity diminishes and becomes more restricted. An intriguing aspect of our findings, compared to previously studied models [106−109], is the elimination of the forbidden region in this model. In other words, in the investigated areas, there is no region where both the ϕ function and the metric function simultaneously lack solutions. Also, at the end, we fully checked the curvatures singularities, and energy conditions for the mentioned black hole In conclusion, our research provides a comprehensive understanding of the topological properties of Kiselev-AdS black holes within $f(R, T)$ gravity. On the other hand, the influence of the parameters ω and γ on the topological charges and the conditions for the existence of the photon sphere offer valuable insights into the complex behavior of these black holes. These findings pave the way for further exploration and deeper understanding of the topological aspects of black holes in modified gravity theories.

While our study has shed light on several important aspects of Kiselev-AdS black holes within $f(R, T)$ gravity, it also opens up several intriguing questions for future research:

1. How do these findings translate to other modified gravity theories? Can similar topological classifications be observed in other contexts, such as $f(R)$ or $f(T)$ gravity alone?

2. What are the dynamic properties of these black holes under perturbations? How do the topological charges evolve, and what implications does this have for black hole stability?

3. How do quantum corrections affect the topological properties of these black holes? Can we extend our classical findings to a quantum regime, and what new phenomena might emerge?

4. How do these topological properties manifest in higher-dimensional black holes? Does the dimensionality of spacetime introduce new complexities or simplifications in the topological classifications?

Exploring these questions will not only deepen our understanding of black holes in modified gravity theories but also potentially uncover new and exciting aspects of gravitational physics.