Properties of Spherically Symmetric Black Hole in the Framework of GBD Modified Gravitational Theory

1. Introduction

Einstein established the famous theory of special relativity in 1905, which laid the foundation for the relativistic view of spacetime. Based on his deep thinking about special relativity, Einstein established the theory of general relativity (GR) in 1915. GR is considered one of the most elegant theories in the field of physics, and has received significant experimental support. For example: (1) General relativity passed rigorous tests within the solar system; (2) Its applications in cosmology have also been extremely successful; (3) The detection of gravitational waves by the advanced Laser Interferometer Gravitational-Wave Observatory (LIGO) in 2015 provided strong support for the theoretical predictions made by GR 100 years ago [1]; (4) The fact that the Event Horizon Telescope [2] discovered the existence of supermassive black holes in the center of the Milky Way galaxy in 2019 greatly enhanced interest and confidence for us in the research of general relativity theory. However, it is also evident that many problems and challenges remain in the study of gravitational physics, such as the problem of gravitational quantization, the nature of dark energy, the problem of inflationary universe and the problem of space-time singularity, etc. To address these issues, many modified gravity theories have been developed and discussed, such as f (R) theory [3,4], f (T) theory [5,6], f (R, T) theory [7], and f (G) theory [8], among others. The exploration of modified theories of gravity remain a hot topic in current research of gravitational physics. As one type of modified gravity theory, the Generalized Brans-Dicke (GBD) theory has been studied in cosmology, gravitational wave physics, and other fields [9,10,11]. In this paper we will explore some issues in the framework of GBD theory.

Black hole (BH) is one of the important prediction of general relativity. After years of observation and exploration, people finally captured the first image of a BH in 2019 through the Event Horizon Telescope. For many years, the research on BH physics has attracted the attention of lots of physicists and astronomers. Different types of black holes, such as static BHs [12], dynamic BHs [13], spherically symmetric BHs [14], axially symmetric BHs [15], and exotic BHs [16], have been intensively discussed. The study of the thermodynamic laws of BHs areas suggests that black holes, as special celestial bodies, seem to have thermal properties. It is well known that thermal properties of BHs under different types has been widely explored. The research on the corrected BH entropy is a hot topic in BH thermodynamics, and some studies on the correction of black hole entropy under modified gravity theory have been investigated, such as conformal field theory [17,18], string theory [19,20], and others. This article mainly explore the relevant properties of thermodynamic correction entropy of static spherically symmetric BHs within the framework of GBD modified gravity theory. Additionally, the stability analysis of the geodesic motion of particles around a BH is an important research topic in the field of gravitational physics, which plays a vital role in exploring the properties of black holes and gravity. Under different contexts of BH, people have conducted extensive research on the stability of black holes through the application of geodesic deviation equation, such as charged black holes in f (T) theory [21], rotating (anti-) de-Sitter black holes in f (R) theory [22], and non-trivial black holes [23], etc [24,25,26,27]. In this paper, we also discuss the stability of the geodesic motion of particles around black holes in the GBD gravity theory.

The structure of our paper is as follows. The first section is an introduction. In the second part of this paper, we focus on the correction of thermodynamic entropy of a static spherically symmetric BH within the framework of GBD modified gravity theory. We analyse the stability of geodesic motion of particles in black holes in the third section. The fourth part is the conclusion of this article.

2. $\mathrm{Entropy}\mathrm{Correction}\mathrm{of}\mathrm{Black}\mathrm{Hole}\mathrm{Thermodynamics}\mathrm{in}\mathrm{GBD}\mathrm{modified}\mathrm{gravity}\mathrm{theory}$

We first briefly introduce the GBD theory, whose action is written as [9]:

$$S=\int d^{4}x\sqrt{-g}[\varphi f\left(R\right)-\frac{\omega }{\varphi }{\partial }_{\mu }\varphi {\partial }^{\mu }\varphi +\frac{16\pi }{c^4}{L}_m(g_{\mu \nu },{\psi }_m)].$$

(1)

Here, $\varphi$ denotes the Brans-Dicke scalar field, $\omega$ is the coupling constant, $L_m$ represents Lagrange density of matter field. By applying the variational principle, the gravitational field equation and scalar field equation of GBD theory can be derived as follows:

$$\begin{array}{c}\hfill \varphi [f_R{R}_{\mu \nu }-\frac{1}{2}f\left(R\right)g_{\mu \nu }]-({\nabla }_{\mu }{\nabla }_{\nu }-g_{\mu \nu }\square )\left(\varphi f_R\right)+\frac{1}{2}\frac{\omega }{\varphi }{g}_{\mu \nu }{\partial }_{\sigma }{\partial }^{\sigma }\varphi -\frac{\omega }{\varphi }{\partial }_{\mu }\varphi {\partial }_{\nu }\varphi =8\pi T_{\mu \nu }.\end{array}$$

(2)

$$\begin{array}{c}\hfill f\left(R\right)+2\omega \frac{\square \varphi }{\varphi }-\frac{\omega }{{\varphi }^2}{\partial }_{\mu }\varphi {\partial }^{\mu }\varphi =0.\end{array}$$

(3)

We consider a specific parametric model[28,29]: $\varphi \left(r\right)={\varphi }_0{r}^{-a}$ and $f\left(R\right)=R+\beta R^{(-n)}$. Under background of the spherically symmetric space-time line elements:

$$\begin{array}{c}\hfill d{\tau }^2=-B\left(r\right)d{t}^2+A\left(r\right)d{r}^2+r^{2}d{\theta }^2+r^{2}sin{\theta }^{2}d{\phi }^2.\end{array}$$

(4)

with $A\left(r\right)=B{\left(r\right)}^{-1}$, the solution of the field equation can be calculated when $n\ge -1$

$$\begin{array}{c}\hfill B\left(r\right)=1+\frac{C_1}{r^2}+\frac{C_2}{r^2}.\end{array}$$

(5)

where $C_1$ and $C_2$ are two constant parameters. Obviously, expression (5) is similar in form to the Reissner-Nordstrom solution for charged celestial bodies in Einstein’s general relativity (if the integration constants are set to: $C_1=-2GM$, $C_2=\frac{G{Q}^2}{4\pi }$, then the two solutions are the same in form). Therefore, in the subsequent discussion,we consider that the values of these two parameters satisfy: $C_1\le 0$ and $C_2\le 0$.

The study of black hole thermodynamics and entropy correction has always been an important issue in the gravitational physics. Next, we investigate the corrected entropy of spherically symmetric BH in the framework of GBD theory. Consider the partition function in the form [30]:

$$\begin{array}{c}\hfill Z\left(\beta \right)={\int }_0^{\infty }\rho \left(E\right)e^{(-\beta E)}dE.\end{array}$$

(6)

where $T=\frac{1}{\beta }$ is the temperature in units of the Boltzmann constant ${\kappa }_B$. $\rho \left(E\right)$ is the density of state, which can be given from (6) by performing the inverse Laplace transform (keeping E fixed) [31,32]:

$$\begin{array}{c}\hfill \rho \left(E\right)=\frac{1}{2\pi i}{\int }_{c-i\infty }^{c+i\infty }Z\left(\beta \right)e^{\beta E}d\beta =\frac{1}{2\pi i}{\int }_{c-i\infty }^{c+i\infty }{e}^{S\beta }d\beta .\end{array}$$

(7)

In Eq. (7), the exact entropy as a function of temperature (not just at equilibrium) can be written as[33]

$$\begin{array}{c}\hfill S\left(\beta \right)=lnZ\left(\beta \right)+\beta E.\end{array}$$

(8)

It is formally defined as the sum of entropies of subsystems of the thermodynamic system, which are small enough to be themselves in equilibrium [33]. Considering the method of steepest descent around the saddle point: ${\beta }_0=\frac{1}{T_0}$ , we have $S_0^{{}^{\prime }}\equiv \frac{\partial S\left(\beta \right)}{\partial \beta }{|}_{\beta ={\beta }_0}=0$ with $T_0$ being the equilibrium temperature. On the other hand, we can receive the following form of entropy function by expanding it about $\beta ={\beta }_0$,

$$\begin{array}{c}\hfill S\left(\beta \right)=S\left({\beta }_0\right)+S^{{}^{\prime }}\left({\beta }_0\right)(\beta -{\beta }_0)+\frac{S^{{}^{\prime \prime }}\left({\beta }_0\right)}{2}{(\beta -{\beta }_0)}^2+....\end{array}$$

(9)

Then Eq.(9) can be simplified as

$$\begin{array}{c}\hfill S\left(\beta \right)=S_0+\frac{1}{2}{(\beta -{\beta }_0)}^2{S}^{{}^{\prime \prime }}\left({\beta }_0\right)+...\end{array}$$

(10)

with

$$\begin{array}{c}\hfill S_0\equiv S\left({\beta }_0\right),S_0^{{}^{\prime \prime }}\equiv \frac{{\partial }^{2}S\left(\beta \right)}{\partial {\beta }^2}{|}_{\beta ={\beta }_0}.\end{array}$$

(11)

where $S\left(\beta \right)$ represents the entropy at any temperature, and $S_0$ is the entropy calculated by the Beckenstein-Hawking area law. As can be seen from Equation(9), the correction of entropy is only described by the term of $S_0^{{}^{\prime \prime }}$. Using Equation(8), it can be deduced that

$$\begin{array}{c}\hfill S^{{}^{\prime }}\left(\beta \right)=\frac{1}{Z\left(\beta \right)}\frac{\partial Z\left(\beta \right)}{\partial \beta }+E.\end{array}$$

(12)

$$\begin{array}{c}\hfill S^{{}^{\prime \prime }}\left(\beta \right)=-\frac{1}{Z^2\left(\beta \right)}{\left(\frac{\partial Z\left(\beta \right)}{\partial \beta }\right)}^2+\frac{1}{Z\left(\beta \right)}\frac{{\partial }^{2}Z\left(\beta \right)}{\partial {\beta }^2}.\end{array}$$

(13)

Consider the energy of a canonical ensemble [33]:

$$\begin{array}{c}\hfill E=\overline{E}=-\frac{\partial }{\partial \beta }lnZ=-\frac{1}{Z}\frac{\partial Z}{\partial \beta }.\end{array}$$

(14)

$$\begin{array}{c}\hfill \frac{\partial \overline{E}}{\partial \beta }=\frac{1}{Z^2\left(\beta \right)}{\left(\frac{\partial Z\left(\beta \right)}{\partial \beta }\right)}^2-\frac{1}{Z\left(\beta \right)}\frac{{\partial }^{2}Z\left(\beta \right)}{\partial {\beta }^2}.\end{array}$$

(15)

Then we have:

$$\begin{array}{c}\hfill \frac{\partial \overline{E}}{\partial \beta }=-[<E^2>-<E{>}^2]=-{\kappa }_B{T}^{2}C.\end{array}$$

(16)

here C is the heat capacity. So, we finally get

$$\begin{array}{c}\hfill S_0^{{}^{\prime \prime }}=T^{2}C.\end{array}$$

(17)

Applying it to a thermodynamic system of BH, the saddle point is represented as ${\beta }_0=\frac{1}{T_H}$ by taking Hawking temperature $T_H$ instead of T. Following the method used in reference[34], the entropy can be rewritten by introducing a parameter $\gamma$ to track the correction term

$$\begin{array}{c}\hfill S=S_0-\gamma \left(C{T}^2\right).\end{array}$$

(18)

For GBD modified gravity theory, $S_0$ can be derived as follows

$$\begin{array}{c}\hfill S_0=\pi r_{+}^2{\phi }_0\frac{1}{r_{+}^{\alpha }}(1-2\frac{\beta }{R^{n+1}}).\end{array}$$

(19)

Comparing the spherically symmetric BH solution(5) under this theory with the Schwarzschild BH solution in the GR theory, the integration constant $C_1$ can be easily expressed as $C_1=-2M$, namely.

$$\begin{array}{c}\hfill B\left(r\right)=1-\frac{2M}{r}+\frac{C_2}{r^2}.\end{array}$$

(20)

Using equations(19) and (20), we can derive:

$$\begin{array}{c}\hfill S_0=\pi r_{+}^{2-\alpha }.\end{array}$$

(21)

Let ${\phi }_0=1$ , expression of the event horizon radius can be written as

$$\begin{array}{c}\hfill r_{+}={\pi }^{-\frac{1}{2-\alpha }}{S}_0^{\frac{1}{2-\alpha }}.\end{array}$$

(22)

Obviously, the relationship between the BH mass and the event horizon radius in GBD theory meets

$$\begin{array}{c}\hfill M=\frac{r_{+}}{2}+\frac{C_2}{2r_{+}}\end{array}$$

(23)

Using expression (23), we draw the picture of the black hole mass relative to the event horizon radius for different parameter values $C_2$, as shown in Figure 1. From Figure 1, we observe that for the case $C_2=15$, the black hole mass reaches the minimum value $M\approx 3.87$ when the event horizon radius is $r_{+}\approx 3.87$.

Combining Eqs. (21) and(23), we obtain

$$\begin{array}{c}\hfill M=\frac{1}{2}+{\pi }^{\frac{1}{2-\alpha }}{S}_0^{-\frac{1}{2-\alpha }}({\pi }^{\frac{1}{2-\alpha }}{S}_0^{\frac{2}{2-\alpha }}+C_2).\end{array}$$

(24)

Then we can get the concrete expressions of temperature and heat capacity for the GBD black hole as follows:

$$\begin{array}{c}\hfill T=\frac{\partial M}{\partial S_0}=\frac{{\pi }^{-\frac{1}{2-\alpha }}{S}_0^{\frac{-1+\alpha }{2-\alpha }}}{2-\alpha }-\frac{{\pi }^{\frac{1}{2-\alpha }}{S}_0^{\frac{\alpha -3}{2-\alpha }}({\pi }^{-\frac{1}{2-\alpha }}{S}_0^{\frac{2}{2-\alpha }}+C_2)}{2*(2-\alpha )}\end{array}$$

(25)

$$\begin{array}{c}\hfill C=\frac{\partial M}{\partial T}=T\frac{\partial S_0}{\partial T}=-\frac{S_0*(-2+\alpha )({\pi }^{\frac{2}{-2+\alpha }}(-1+\alpha )-S_0^{\frac{2}{-2+\alpha }}{C}_2)}{{\pi }^{-\frac{2}{2-\alpha }}(-1+\alpha )-S_0^{\frac{2}{-2+\alpha }}(-3+\alpha )C_2}.\end{array}$$

(26)

According to these expressions, we plot the variation of BH temperature (upper graph) and heat capacity (lower graph) with respect to the event horizon radius in Figure 2. From Figure 2 (upper), we observe that the BH temperature is positive only within a specific range of the event horizon radius. For instance, for the parameter values $\alpha =0.01$ and $C_2=15$, the BH temperature is positive when $r_{+}>3.87$; whereas for $r_{+}<3.87$, it corresponds to a physically meaningless negative temperature region. Moreover, the black hole temperature reaches its maximum value when the event horizon radius is $r_{+}\approx 6.73$. The lower graph of Figure 2 shows that the black hole undergoes a phase transition in the framework of GBD modified gravity theory, and the location of the phase transition depends on the model parameter values. For example, for $\alpha =0.01$ and $C_2=15$, the black hole is in an unstable phase with negative heat capacity in the range $0<r_{+}<3.87$, and the heat capacity of BH has its minimum value at the event horizon radius $r_{+}\approx 2.87$ ; On the other hand, for the range $r>3.87$, the black hole is in a stable phase with positive heat capacity. Clearly, the black hole undergoes a phase transition at $r_{+}\approx 3.87$, where the BH heat capacity $C=0$.

By combining equations(18),(25) and (26), we can derive the corrected entropy expression for black holes under GBD theory, given by:

$$\begin{array}{c}\hfill S=S_0+\frac{{\left(r_{+}^{2-\alpha }\right)}^{\frac{\alpha }{2-\alpha }}\gamma {[-1+{\left(r_{+}^{2-\alpha }\right)}^{\frac{2}{-2+\alpha }}{C}_2]}^3}{4\pi (-2+\alpha )[1-\alpha +{\left(r_{+}^{2-\alpha }\right)}^{\frac{2}{-2+\alpha }}(-3+\alpha )C_2]}.\end{array}$$

(27)

Taking the parameters $\gamma =1$, $\alpha =0.01$ and $C_2=15$ into account, we can use equation (27) to plot the corrected entropy as a function of the event horizon radius (as shown in Figure 3). From Figure 3, we observe that the entropy S with the correction term decreases rapidly with increasing event horizon radius and reaches its minimum value $S\approx 5.44$ at $r_{+}\approx 1.09$, then gradually increases. The plot of $S_0$ indicates that it increases monotonically with $r_{+}$ increasing. Comparing the two plots of S and $S_0$ , we find that the correction term has a significant effect mainly in the region of small values $r_{+}$, and the entropy-increasing effect is very prominent. As for large values $r_{+}$, the effect of the correction term can be almost ignored (the two curves almost overlap). In addition, we also calculate the influence of other model parameter values on the entropy. From Figure 3 (left), we can see that the variation of the model parameter $C_2$ has little effect on the corrected entropy when $r_{+}$ is large.

3. $\mathrm{Analysis}\mathrm{of}\mathrm{geodesic}\mathrm{motion}\mathrm{stability}\mathrm{of}\mathrm{particles}\mathrm{in}\mathrm{black}\mathrm{holes}\mathrm{under}\mathrm{GBD}\mathrm{theory}$

In this section, we analyze the motion of particles and their properties around black holes in the GBD theoretical framework. The motion of a free particle in a gravitational field is described by the geodesic equation:

$$\begin{array}{c}\hfill \frac{d^2{x}^{\mu }}{d{p}^2}+{\Gamma }_{\nu \lambda }^{\mu }\frac{d{x}^{\nu }}{dp}\frac{d{x}^{\lambda }}{dp}=0.\end{array}$$

(28)

Here, p represents the orbital parameter. The component forms of geodesic equation can be written as

$$\begin{array}{c}\hfill \frac{d^{2}t}{d{p}^2}+\frac{B^{{}^{\prime }}}{B}\frac{dt}{dp}\frac{dr}{dp}=0.\end{array}$$

(29)

$$\begin{array}{c}\hfill \frac{d^{2}r}{d{p}^2}+\frac{B{B}^{{}^{\prime }}}{2}{\left(\frac{dt}{dp}\right)}^2-\frac{1}{2B}{\left(\frac{dr}{dp}\right)}^2-rB{\left(\frac{d\theta }{dp}\right)}^2-rBsi{n}^2\theta {\left(\frac{d\phi }{dp}\right)}^2=0.\end{array}$$

(30)

$$\begin{array}{c}\hfill \frac{d^2\theta }{d{p}^2}+\frac{2}{r}\frac{d\theta }{dp}\frac{dr}{dp}-sin\theta cos\theta {\left(\frac{d\phi }{dp}\right)}^2=0.\end{array}$$

(31)

$$\begin{array}{c}\hfill \frac{d^2\phi }{d{p}^2}+\frac{2}{r}\frac{d\phi }{dp}\frac{dr}{dp}+2cot\theta \frac{d\phi }{dp}\frac{d\theta }{dp}=0.\end{array}$$

(32)

Where the prime represents the derivative with respect to the radial coordinate r. For a gravitational field with spherical symmetry, without loss of generality, we select the initial position and velocity of the particle to be on the equatorial plane, i.e., $\theta =\frac{\pi }{2}$ and $\frac{d\theta }{dp}=0$. Using equation (31), we obtain $\frac{d^2\theta }{d{p}^2}=0$, which indicates that the particle motion will always remain on the equatorial plane. For particles with non-zero rest mass, the orbital parameter is taken to be the proper time $\tau$. Further derivation leads to two equations for the motion of particles in a gravitational field:

$$\begin{array}{c}\hfill \frac{dt}{d\tau }=\frac{E}{B}.\end{array}$$

(33)

$$\begin{array}{c}\hfill r^2\frac{d\phi }{d\tau }=J.\end{array}$$

(34)

Here E and J are two constants of integration, representing the energy and angular momentum of a unit mass particle. Additionally, using the normalization condition of the four-velocity $g_{\mu \nu }{U}^{\mu }{U}^{\nu }=-1$ and equations (33) and (34), we obtain:

$$\begin{array}{c}\hfill {\left(\frac{dr}{d\tau }\right)}^2=E^2-B\left(r\right)(1+\frac{J^2}{r^2}).\end{array}$$

(35)

We define the effective potential as:

$$\begin{array}{c}\hfill V_{eff}\left(r\right)\equiv B\left(r\right)(1+\frac{J^2}{r^2}).\end{array}$$

(36)

Then equation (35) becomes:

$$\begin{array}{c}\hfill {\left(\frac{dr}{d\tau }\right)}^2=E^2-V_{eff}\left(r\right).\end{array}$$

(37)

Clearly, the properties of the gravitational potential in the GBD theory can be described by equation (36), which depends on the relative position and angular momentum of the particles. Combining equations (5) and (36), we see that when the radius is $r\to \infty$ , the effective potential is $V_{eff}=1$ . Furthermore, equation (37) indicates that when $E^2=V_{eff}$, i.e., $\frac{dr}{d\tau }=0$, the orbital radius r is constant, and the trajectory of particle is circular.

In black hole physics, it is meaningful to discuss the innermost stable circular orbit (ISCO). The ISCO is determined by the following expressions:

$$\begin{array}{c}\hfill \frac{d{V}_{eff}}{dr}{|}_{r=r_{ISCO}}=0.\end{array}$$

(38)

$$\begin{array}{c}\hfill \frac{d^2{V}_{eff}}{d{r}^2}{|}_{r=r_{ISCO}}=0.\end{array}$$

(39)

Therefore, for the GBD theory, we derive the relationship that the ISCO needs to satisfy:

$$\begin{array}{c}\hfill \frac{3r^2{C}_1^2+8C_2^2+C_1(r^3+9r{C}_2)}{r^2(r{C}_1+2C_2)}=0\end{array}$$

(40)

Solving the above equation numerically, we can plot the variation of the ISCO radius with respect to the parameters (as shown in Figure 4). Our calculations show that the radius of the ISCO in the GBD modified gravity theory framework will continuously increase with increasing values of the parameters $C_1$ and $C_2$.

We also investigate the properties of angular momentum of particles in the GBD theory. Considering the condition of circular orbits in the equatorial plane:

$$\begin{array}{c}\hfill \theta =\frac{\pi }{2},\frac{d\theta }{d\tau }=0,\frac{dr}{d\tau }=0\end{array}$$

(41)

we obtain the expressions for t and $\phi$ by calculating the components of the motion equation (28):

$$\begin{array}{c}\hfill {\left(\frac{d\phi }{d\tau }\right)}^2=\frac{B^{\prime }\left(r\right)}{r[2B\left(r\right)-r{B}^{\prime }\left(r\right)]}\end{array}$$

(42)

$$\begin{array}{c}\hfill {\left(\frac{dt}{d\tau }\right)}^2=\frac{2}{B\left(r\right)-r{B}^{\prime }\left(r\right)}\end{array}$$

(43)

From equations (42) and(43), we derive the expression for angular velocity as:

$$\begin{array}{c}\hfill \Omega =\frac{\dot{\phi }}{\dot{t}}=\sqrt{\frac{{\partial }_{r}B\left(r\right)}{2r}}\end{array}$$

(44)

Furthermore, substituting equation (5) into equation (44) yields

$$\begin{array}{c}\hfill \Omega =\sqrt{\frac{r^3-C_{1}r-2C_2}{2r^4}}\end{array}$$

(45)

Using equation (45), we plot the picture of angular velocity as a function of the radius r in Figure 5, where parameter values $C_1=-10$ and $C_2=15$ are selected (corresponding to the two thin-solid lines in Figure 5). For this case we can see that, for these parameter values, the angular velocity sharply increases in the range of radius r :2.09-3.05, reaches its maximum value $\Omega \approx 0.41$ at radius $r\approx 3.05$, and then $\Omega$ decreases to a steady state. In addition, to discuss the variation of angular velocity with respect to the parameters $C_1$ and $C_2$ , we also plot the curves for different values of these parameters in Figure 5. For example, in the left panel of Figure 5, we consider $C_2=15$, $C_1$ with values of -8, -10, and -12, respectively. In the right panel of Figure 5, we consider $C_1=-10$, $C_2$ with values of 10, 15, and 20, respectively. From Figure 5, we observe that the values of angular velocity decrease with increasing values of $C_1$ or $C_2$.

Also, we can derive to get the specific expressions of energy and angular momentum of unit mass particles in GBD modified theory as follows:

$$\begin{array}{c}\hfill E=\frac{r^2+C_{1}r+C_2}{r^2}\sqrt{\frac{2r^2}{2r^2+3C_{1}r+4C_2}}\end{array}$$

(46)

$$\begin{array}{c}\hfill J=\sqrt{-\frac{r^2(C_{1}r+2C_2)}{2r^2+3C_{1}r+4C_2}}\end{array}$$

(47)

Figure 6 shows the variation of energy and angular momentum of a unit mass particle with radius. From the figure, it can be seen that when the parameter values are taken as $C_1=-10$ and $C_2=15$ , the energy sharply decreases in the radius interval 12.62-24.86, and reaches its minimum value $E\approx 0.93$ at the radius $r\approx 24.86$. It then gradually increases. The variation trend of angular momentum is similar to that of energy, also reaching its minimum value $J\approx 15.67$ at the radius $r\approx 24.86$. The variation of E and J for different parameter values $C_1$ and $C_2$ can be seen in detail in Figure 6. In the plots on the left, takes the value of $C_2=15$, while in the plots on the right, takes the value of $C_1=-10$.

To study the issue of circular orbit stability of particles moving around a spherically symmetric black hole in GBD theory, we consider the application of the geodesic deviation equation. The geodesic deviation equation is given by:

$$\begin{array}{c}\hfill \frac{d^2{\xi }^{\alpha }}{d{\tau }^2}+2{\Gamma }_{\mu \nu }^{\alpha }\frac{d{x}^{\mu }}{d\tau }\frac{d{\xi }^{\nu }}{d\tau }+\frac{d{x}^{\mu }}{d\tau }\frac{d{x}^{\nu }}{d\tau }{\xi }^{\beta }{\partial }_{\beta }{\Gamma }_{\mu \nu }^{\beta }=0\end{array}$$

(48)

where ${\xi }^{\alpha }$is the deviation four-vector. By substituting the spherically symmetric line element (4) and using equations (41)-(43), we can derive the expressions for the geodesic deviation equation components:

$$\begin{array}{c}\hfill \frac{d^2{\xi }^0}{d{\phi }^2}+\frac{B^{\prime }}{B}\frac{dt}{d\phi }\frac{d{\xi }^1}{d\phi }=0\end{array}$$

(49)

$$\begin{array}{c}\hfill \frac{d^2{\xi }^1}{d{\phi }^2}+B{B}^{\prime }\frac{dt}{d\phi }\frac{d{\xi }^0}{d\phi }-2rB\frac{d{\xi }^3}{d\phi }+[\frac{1}{2}{\left(\frac{dt}{d\phi }\right)}^2(B^{{\prime }^2}+B{B}^{\prime \prime })-(B+r{B}^{\prime })]{\xi }^1=0\end{array}$$

(50)

$$\begin{array}{c}\hfill \frac{d^2{\xi }^2}{d{\phi }^2}+{\xi }^2=0\end{array}$$

(51)

$$\begin{array}{c}\hfill \frac{d^2{\xi }^3}{d{\phi }^2}+\frac{2}{r}\frac{d{\xi }^1}{d\phi }=0\end{array}$$

(52)

Obviously, the solution to equation(51) can be expressed as ${\xi }^2={\zeta }^2{e}^{i\phi }$. This indicates that the circular orbits of test particles initially in the equatorial plane will undergo harmonic vibration under perturbations. Therefore, the circular orbit of particle motion is stable. For other equations, assuming the solution takes the form:

$$\begin{array}{c}\hfill {\xi }^0={\zeta }^0{e}^{i\omega \phi }\end{array}$$

(53)

$$\begin{array}{c}\hfill {\xi }^1={\zeta }^1{e}^{i\omega \phi }\end{array}$$

(54)

$$\begin{array}{c}\hfill {\xi }^3={\zeta }^3{e}^{i\omega \phi }\end{array}$$

(55)

then substituting equations(53)-(55) into equations (49), (50), and(52), and considering the stability requirement for circular orbit motion, we derive the following constraints:

$$\begin{array}{c}\hfill {\omega }^2=3B-2r{B}^{\prime }+\frac{rB{B}^{\prime \prime }}{B^{\prime }}>0\end{array}$$

(56)

By substituting expression (5) into the above equation(56), we show in Figure 7 the dependence of the parameter ${\omega }^2$ on the circular orbit radius r. When the integral constants are taken as $C_1=-10$ and $C_2=15$, we find that the stable region of circular orbits ${\omega }^2>0$ is : $r>24.86$. The effects of other model parameter values on the parameter ${\omega }^2$ and the corresponding stable circular orbit regions are shown in Figure 6 (the left figure with $C_2$ taken as 15 and the right figure with $C_1$ taken as -10).

4. Conclusion

The emergence of challenging problems such as gravity quantization and the origin of dark matter and dark energy has provided motivation for finding gravity theories beyond Einstein’s general relativity. Researchers in the field of gravity have made many efforts and practices in exploring modified or extended theories of general relativity, and the study of applying modified gravity theories to astrophysics and cosmology has increasingly received people’s attention. This article focuses on the properties of a BH solution from the so-called GBD theory. We investigates the thermodynamic corrected entropy problem of static spherically symmetric BHs and the stability of particle geodesics in BHs under the framework of the GBD modified gravity theory. Since studying the structure of geodesics in strong gravitational fields plays an important role in exploring the physical characteristics of compact objects, we firstly quantitatively analyze the contribution of the corrected term to the entropy in the GBD gravity. It is found that the effect of the correction term on entropy is significant when the value of the event horizon radius $r_{+}$ is small, while its effect can be negligible when $r_{+}$ is large. Secondly, the article also analyzes and shows the dependence relationship of black hole mass, temperature, heat capacity, etc. on the event horizon radius. Based on the motion equations, we investigate the relevant properties of physical quantities such as energy, gravitational potential energy, angular velocity, and angular momentum of test particles under the GBD theory. Finally, the article applies the geodesic deviation equation to analyze the stable circular orbit of test particles in the GBD theory and finds that the particle undergoes harmonic motion in the equatorial plane, and the range of the restricted area that the radial coordinate needs to satisfy to achieve stable circular orbit motion is given.