Analytically Approximation Solution to R 2 Gravity

.


Introduction
More than a century ago, Newtonian gravity failed to explain the anomalies in the orbits of Mercury.In 1915, Einstein introduced his theory of General Relativity (GR), which not only described the answer to the physics problems of that time, including the issue of Mercury but also its ability in its interesting predictions like the gravitational lensing, gravitational waves, black holes, etc.However, despite all of these successes, in recent years this theory, failed to be consistent with observation and was not able to describe the dark energy and dark matter that are confirmed by observations.Also, GR was not renormalizable and, therefore, could not be quantized by using the conventional methods in quantum field theory.It was shown that if one demands the renormalization to be satisfied at one-loop, the Einstein-Hilbert action should be supplemented by higher-order curvature terms [1].Later on, Einstein's theory of GR had many forms of modifications such as: f (R) gravity with R as the Ricci scalar [2], [3], [4], f (R, T ) gravity, where T is the trace of energy-momentum tensor [5], [6], f (G) gravity with G the Gauss-Bonnet invariant [7] etc.All of these modified theories have received much attention to investigating the accelerated expansion of our universe, flat rotation curves of galaxies, wormhole behavior [8], [9], [10], [11].One of the familiar extensions of GR is to include terms such as where α i and β i are constants, and α i = β i = 0 corresponds to the Einstein-Hilbert action.The special case α 2 = 0, α i = 0(i = 2), and β i = 0 is R 2 gravity and known as the Starobinsky model [12].In f (R) gravity, the interaction of spacetime and matter is different from GR.In GR, gravity appears as the curvature of spacetime, where the source of this curvature is all forms of matter.
In the absence of any mass or energy, spacetime can become completely flat.What f (R) gravity does is allow spacetime to act as a source of its curvature so that there can still be some curvature even if spacetime is empty of matter.Therefore, as the universe expands and the matter density decreases, some curvature may remain that is capable of driving cosmic acceleration.Many applications have been developed in the framework of the f (R) gravity such as earlytime inflation [12], [13], [14], cosmological phases [15], [16], [17], [18], [19], [20], [21], gravitational wave detection [22], [23], the stability of the solutions [24], [25], [26] and other different issues [27].Interestingly, the duality between gravity and quantum field theory that has been introduced in [28] and its applications considered in different theories and f (R) gravity such as [29]- [45].In addition, the different solutions of f (R) gravity have been studied.Among these solutions, the authors in [46] and [47] deduced the static spherically symmetric solutions by involving a perfect fluid.The black hole solution with/without electric charge is presented in [48].Many analytic spherically symmetric solutions are derived in [50], [51], [52].
In this work, following the papers of quadratic [53] and cubic gravity [54], we concentrated on the analytic approximate black hole solutions of f (R) gravity using a continued-fraction expansion and their stability through Quasi-Normal Modes (QNMs) and thermodynamic.
The paper is organized as follows: In Section 2, first, we give a brief review of f (R) gravity and apply the field equations of the theory to static-spherically symmetric spacetime.Then, to obtain the full black hole solution, we combine the asymptotic and near horizon solutions using a continued-fraction expansion.In Section 3, using the QNMs we study the dynamical stability of the black hole solutions.Finally, in Section 4 we give our concluding remarks.

Basic equations
The action of f (R) gravity can be written as The variation to action gives the vacuum field equations as where is the d'Alembertian operator and f R = df /dR.The trace of the field equations takes the form Inserting Eq. (4) in Eq. (3) we get In this paper, we consider f (R) as follows Inserting into the field equation ( 5), one can get For this model α > 0 from the Dolgov and Kawasaki stability point of view [55].We consider the following spherically symmetric and static line element for describing the geometry of spacetime As we know, generic static, spherically, and symmetric metrics do not need to obey g tt g rr = −1 necessarily, but field equation ( 7) admits solutions with this property [56], [57], to which case we shall restrict our consideration here.By inserting the metric into the field equations (7), the differential equations for h(r) become We would like to emphasize that, utilizing the other components of field equations gives the same result.Expanding the function h(r) around the event horizon r + we obtain and then inserting these expressions into equations ( 9), we find where r + , h 1 and h 2 are undetermined constants of integration.The other near horizon constants are provided in the Appendix B. In the large r limit, we linearize the field equations about the flat background where H(r) is to be determined by the field equations, and we linearize the differential equation by keeping terms only to order .The resulting differential equation for H(r) takes the form In the large r limit we obtain that can be solved as here, the pF q(a, b, z) is the generalized hypergeometric function [58], [59].For c 1 = c 2 = c 3 = H 1 and in the large r, H(r) becomes Finally, the solution is We wish to obtain an approximate analytic solution that is valid near the horizon and at large r.To reach this goal, we employ a continued-fraction expansion and write with Even at lowest order, the continued fraction approximation does a good job of approximating the solution everywhere outside the horizon.This only gets better as more terms are included.
Here, we truncated the continued-fraction expansion at order 4.This method give an accurate analytic expression approximating for the metric for the whole space outside the event horizon, and not only near the black hole or far from it (r + ≤ r < ∞).By expanding (18) near the horizon (x → 0) and the asymptotic region (x → 1) we obtain for the lowest order expansion coefficients, with the remaining a i given in terms of (r + , h 1 , h 2 ).Also, we provided these expressions in the Appendix A. The result is an approximate analytic solution for metric functions everywhere outside the horizon.For a static space time we have a timelike Killing vector ξ = ∂ t everywhere outside the horizon.Therefore, one can obtain We computed the entropy as follows [70], [71] We now consider the thermodynamics of these black hole solutions, whose basic equations are the first law and Smarr formula dM = T dS , where there are no pressure/volume terms since we have set Λ = 0. From Eq. ( 24) we have yielding the mass parameter as a function of the horizon radius.Regarding the asymptotic flat spacetimes, we can define the total mass at spatial infinity as the Arnowitt-Deser-Misner (ADM) mass.Besides, the Komar mass makes sense for calculating total mass of a black hole solutions enjoying a timelike Killing vector.In addition to these two approaches, one can use alternative methods like Misner-Sharp prescription.Although the Misner-Sharp mass is just the Schwarzschild-like mass in the Schwarzschild space, it can be extended to 4-dimensional f(R) gravity (see for example: [65], [66], [67] for higher dimensions).In the Appendix C, we showed that the asymptotic behavior of the mass ( 25) is the same as the mass of the Schwarzschild black hole.Therefore, one can interpret the mass (25) as ADM mass.We now impose the first law (23), which becomes where differential equation that must be satisfied by h 1 (r + ) and h 2 (r + ).Therefore, we achieve the prominent equation as follows In order to solve the differential equation (27), we need to define h 2 .To do so, we assume h 2 = g(α).Recall that, near the horizon, the metric function is expanded as where the constants with i > 2 are determined by the field equations in terms of other parameters.We will demand that this expansion has a smooth α → 0 limit.Therefore, the expansions for the first four terms are Demanding h 3 to have a smooth behavior at α → 0 limit, we must take Then, for h 4 we have Now, by inserting (33) and (34) into the h 5 and h 6 , one can find h 5 = 0 and g = 0.By continuing the procedure in this way, one can find that the other coefficients of expansion should be zero.
Therefore, h 2 can be obtained as By inserting h 2 into the equation ( 27), one can obtain a differential equation only for h 1 .By solving it we obtained three solutions for h 1 h (1) (2) where For α > 0 and c 1 > 0 we have just one real solution and also for α > 0 and c 1 ≤ 0 we have three real solutions.Inserting ( 37)-( 41) into the thermodynamical quantities and plotting them, one can obtain Fig. (1).In this figure, we have illustrated the temperature, mass, and entropy for different values of parameters in the left, middle, and right panels, respectively.As can be seen, there are two kinds of solutions with four branches in such a way that the black dashed lines are the Schwarzschild-like solution and the solid lines non-Schwarzschild-like solution.The black solid line is a non-physical black hole solution due to having a negative temperatures and mass.In figure (2) we have shown the behavior of heat capacity and free energy in the left and right panels, respectively.As can be seen, the heat capacity is negative for all values of r + which shows the black holes locally are unstable.The free energy for a non-Schwarzschild-like solution is decreasing function which shows the solutions globally are stable.In Figure (3) we present the solutions for h(r) with different values of horizons in the left, middle, and right panels that depicting the full continued-fraction solution (18) along with its comparison to the near-horizon solution with blue dashed lines and large-r series expansions solution with black dashed lines.We see that the continued-fraction expansion converges to both of these two approximations.Similar to quadratic gravity, we find two groups of solutions [68], [69], [53].For the first group, the metric function is increasing functions in r ≥ r + .These solutions are generalizations of the Schwarzschild black hole solution and reduce to the Schwarzschild as α → 0. The non-Schwarzschild solutions are physically distinct from the first group, and as α → 0, it does not go back to the Schwarzschild metric.Here, we would like to solve the field equation ( 8) and the metric near the center of a black hole.Therefore, we consider the following expansion around the origin as Similar to previous cases for near horizons and asymptotes, one can obtain the metric as with c 0 = 1 and c 0 = −2/9.The metric is fully determined and there is no free parameter.The Kretschmann scalar and Ricci scalars near the origin behaves as As can be seen, for the case c 0 = 1, the behavior of the metric near the origin corresponds to an Anti-de Sitter space-time and the Kretschmann and Ricci scalar have a finite values (K = 3/200α2 , R = −2/3α) (with a nonzero mass).Regular solutions with an Anti-de Sitter core have been studied with a nonlinear electromagnetic source in [72].But, for the case c 0 = −2/9, the Kretschmann and Ricci scalar do not become zero and there is a singularity.In Figure 4, we have shown the full solutions of the field equations from center to infinity by starting from the near origin metric (43) and asymptotic solution (17), for c 0 = 1 (a solution without horizon 2 ) and for c 0 = −2/9 (black hole case).

Dynamical Stability
In this section, we are interested in calculating the Quasi-Normal Modes (QNMs) of constructed black hole solutions to investigate the dynamical stability of obtained black hole solutions undergoing scalar perturbations.To do so, we used the Klein-Gordon equation for the massless real scalar probe φ(r) as follows φ(r) = 0 .
It is important to note that, since we assume that the metric ( 8) is the vacuum solution of the field equation ( 9) and also the perturbation of the scalar field is second order then can not contribute at the linear perturbation to the field equation [73], [74].After the separation of variables and changing the radial coordinate to the tortoise coordinate, one can find the secondorder Schrodinger-like Ordinary Differential Equation (ODE) for a radial coordinate as follows where here, h(r) has been obtained in (18).The boundary conditions for asymptotically flat black holes are ϕ ∼ e ±iωx , x → ±∞ .
Now, in order to solve the differential equation (47) with the conditions (49), we used the Mashhoon's Method [75], [76].According to this method the QNM's of a potential barrier are related to the bound states of the inverted potential.The frequency of waves must be complex, ω = ω r +iω i .The imaginary and real parts are related to the damping time scale (τ i = 1/ω i ) and oscillation time scale (τ r = 1/ω r ), respectively.The proper QNM's using this method can be obtained as where Our results are presented in tables ( 1), ( 2) and ( 3) and plotted in figures (5).A couple of points should be mentioned.First, from the figures and the tables, one can conclude that the quasinormal frequencies have a positive imaginary part, which shows that both the Schwarzschild and non-Schwarzschild black holes are stable under the scalar perturbations.Second, it is observed for both black holes, with increasing j, the values of the real part of the frequencies increase while the imaginary part of the frequencies decreases, which is a sign of more stability of black holes.Also, as can be seen from the figures and tables with increasing n, the imaginary part of the frequencies increase which shows that the black hole becomes more unstable.

Conclusion
f (R) gravity is a natural generalization of GR and has been received extensive investigation in past decades.In particular, the black hole solutions of this gravity are interesting.To do so, in this paper, we obtained analytical approximate black hole solutions in the framework of R 2 gravity.To achieve this goal, first, we obtained the near horizon and asymptotic solutions and then use these to obtain a complete solution using a continued-fraction expansion.It is important to note that the continued fraction expansion is applicable just from horizon to infinity.Two general types of solutions numerically have been obtained, i.e. black hole solutions and the solutions without horizon and singularity.In addition, the black hole solutions are categorized in two types, i.e.Schwarzschild and non-Schwarzschild-like solutions.Our results showed that there are obvious differences between Schwarzschild and non-Schwarzschild-like solutions.In the non-Schwarzschild-like solutions unlike the Schwarzschild case by increasing the radius of horizon mass and entropy decrease while temperature increases.We studied the stability of the solutions and we found that the solutions thermodynamically and dynamically are stable.We would like to emphasize that in this work similar to our previous works [53], [54], we assumed that the near horizon constant h 1 is a function of r + .This assumption is eligible, because h 1 is proportional to temperature according to Equation (21).
Another interesting direction to extend our work will be to investigate the non-vacuum, rotating black hole, and other solutions of the theory by using the continued-fraction expansion.
If we assume H 1 = −2M , the total energy of black hole becomes which is the same as ADM mass.While the mass (25) in the case of small α 1 becomes for three solutions of h 1 in Eqs. ( 36)- (41).It is straightforward to show that for c 1 = −1/4, the asymptotic value of the mass (59) is the same as Schwarzschild mass.

Figure 1 :
Figure 1: Plots of T , M and S in terms of r + for α = 0.5, c 1 = −0.1406.In all panels the dashed line curves indicate Schwarzschild-like behavior and solid line curves non-Schwarzschildlike behavior.

Figure 2 :
Figure 2: The behavior of heat capacity (left) and free energy (right) in terms of r + for α = 0.5, c 1 = −0.1406.

55 Figure 3 :
Figure 3: The plots of metric for α = 0.5, c 1 = −0.1406.The red solid line is the full continuedfraction solution, the blue dashed line is the near-horizon solution and black dashed line is the large-r solution.

Table 1 :
The values of the quasinormal frequencies in the (small) non-Schwarzschild geometry for the values of parameters r + = 0.5, M = 0.55, α = 0.5.

Table 2 :
The values of the quasi-normal frequencies in the (large) Schwarzschild geometry for the values of parameters r + = 1.5, M = 0.57, α = 0.5.

Table 3 :
The values of the quasinormal frequencies in the (large) non-Schwarzschild geometry for the values of parameters r + = 2.0, M = 0.37, α = 0.5.