Quintessential effects on quasiperiodic oscillations in 4D Einstein–Gauss–Bonnet gravity

Dynamics of test particles in the space-time of a quintessential black hole (BH) in a novel four-dimensional Einstein–Gauss–Bonnet (EGB) theory is studied. First, we analyse the properties of the space-time of the BH and possible values of the coupling α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document} and quintessential parameters q that allows the existence of the event horizon. The scalar invariants of the space-time and effects of the parameters on innermost stable circular orbits (ISCOs) radius are also studied. Moreover, fundamental frequencies such as Keplerian and harmonic oscillations are investigated together with applications to quasi-periodic oscillations (QPOs) where the powerful way to test gravity theories using astrophysical observations of twin-peak QPOs is discussed.


Introduction
Over the last years, we got several observational tests of general relativity [1,2,11,12]. The theory is overall well tested in both weak and strong regimes, and it meets several fundamental problems. Among other issues, it is incompatible with quantum theory, it requires the inclusion of additional concepts (such as dark matter and dark energy) and solutions of its field equations have non-removable singularities. However, theory allows for modifications and alterations to fix those issues and such modified theories do exist in numbers and studying them is important to be able to compare them against observation data and constrain and estimates parameters of the theories.
One of such modifications is Einstein-Gauss-Bonnet (EGB) theory that is defined in higher dimensional space-time (D > 4). According to Lovelock theorem [26] the only theory of gravity in four-dimensional space-time is the general relativity with non-zero cosmological constant. But recently [13] proposed a new four-dimensional Einstein-Gauss-Bonnet solution obtained using a new approach to the Lovelock theorem. To achieve that, the authors rescaled the Gauss-Bonnet term by the factor 1/(D −4) and took the four-dimensional limit. The solution can be used to study the particle motion in this new model. Background space-time in this case represents a non-standard black hole. based on the study, new tests can be devised which will be used to obtain constraints on EGB parameters.
In Ref. [20], gravitational stability was studied by means of photon motion and quasi-normal modes. It was shown that there is a limiting value for coupling constant α (0 < α < 0.15) appearing in the solutions, above which the black hole (BH) becomes gravitationally unstable.
In Ref. [21], authors obtained electrically charged 4D EGB BH solution. Quantum effects of 4D EGB BH have been explored in Ref. [23], while the solution in 3D EGB has been studied in Ref. [22]. Analysing the gravitational collapse of spherical homogeneous dust in 4D EGB gravity authors of [27] had shown its qualitative similarity to the case of Einstein's gravity. Quasi-normal modes and superradiant process in 4D EGB gravity have been explored in Ref. [41]. The stability of the solution in 4D EGB gravity through the odd parity perturbation method has been studied in Authors of Ref. [5] investigated the stability of 4D EGB solution using odd parity perturbation method. Exploring the spinning particle motion around 4D EGB BH Authors of [40] studied motion of spinning particles around 4D EGB BH and had shown the existence of two separate orbits for each value pair of spin and orbital angular momentum. The geodesic motion of test particles has been explored in Ref. [14].
Meanwhile, there are objections to the validity of aforementioned new approach and solutions produced by it. Authors of the original research [13] claimed that 4D solution is not well defined in the comments to the paper. There are arguments, that it is only possible to get solutions for D > 4 when Gauss-Bonnet term is present [15]. Other arguments against the solution of 4D EGB gravity and corresponding approach are presented in Refs. [7,8,37].

Quintessence in 4D Einstein-Gauss-Bonnet black holes
The action for the EGB gravitational theory is described by the following expression: where R is the Ricci scalar, α is the GB constant, L m is the quintessential (anisotropic) matter field whose non-vanishing stress-energy tensor components, L G B is Lagrangian of the EGB theory given as where R μνγ σ is Riemann tensor and R μν is Ricci tensor. The field equation in EGB gravity has the following form: where the Einstein tensor, and the Lanczos tensor, Fig. 1 The dependence of radius of the event horizon on the quintessential parameter q for different values of the parameter ω q and the GB coupling parameter T μν defines the stress tensor of the quintessential field which has the following components in the Kiselev black hole solution [19]: The solution for the space-time around the quintessential black holes in 4D EGB gravitational field theory have been found in Ref. [34] and in spherical coordinates (x α = {t, r, θ, φ}) it has the following form: with the radial metric function where q and ω q are parameters of quintessential field, M is the total mass of the BH. In general, the event horizon of the space-time around a black hole can be found by the condition g tt (r h ) = g rr (r h ) = 0. Due to the complicated form of the solution of the equation for the event horizon, we only provide graphics for the different values of the GB and ω q parameters. Figure 1 presents the behaviour of event horizon radii of quintessential black holes in 4D EGB gravity, with respect to the variation of the quintessential parameter q for different values of parameters α and ω q . It is seen that as q = α = 0 the black lines cross The same relation as shown in Fig.2, but for the event horizon radius, the same as Schwarzschild black hole r h = 2M radius at a corresponding critical value of the quintessential parameter, and the critical value increases as the α parameter is increased as well as the minimum of the event horizon radius.
Next we are interested to determine values of the parameters ω q and q for which the metric (6) can be a black hole solution. One can get the set values of these parameters solving a system of equations f (r ) = f (r ) = 0. However, it is impossible to have analytical solutions for the q and α parameters, due to the complexity of the equation, so, one may analyse this problem graphically and numerically.
Our detailed analysis in Fig.2 show that the maximum values of both α and q parameters and their configurations increase as ω q increases. It is interesting to study conditions on the parameters q and α that allow the black hole solution (6) to have event horizon r h = 2M, in other words, the solution that reflects the Schwarzschild black hole one and does not allow distinguishing the quintessential black hole solution from the Schwarzschild one. Moreover, we have found out that the limit of the metric function given in Eq.(6) at the large values of radial coordinate turns to one for the Schwarzschild case in 4D EGB gravitational theory. It does mean that there are no quintessential field effects at large distances.
In Fig.3, we show similar analysis on the relations provided in Fig.3 for the value of event horizon radius r h = 2M, the same as Schwarzschild black hole event horizon. One can see from the figure that the relation is linear, i.e. as the α parameter increases the quintessential field parameter increases linearly, and the slope is larger at higher values of ω q parameter.

Scalar invariants
In this subsection, we show detailed analysis of the space-time properties around the quintessential black holes. We aimed to study the scalar invariants: Ricci scalar, square of Ricci tensor and Kretschmann scalar.

Ricci scalar
Ricci scalar R = g μν R μν of the space-time around the quintessential black hole has the form where In the GR limit, the pure quintessential black hole's Ricci tensor takes the form Above limit shows that in the absence of the quintessence space-time becomes Ricci flat R = 0. The square of Ricci tensor Next we study the square of Ricci tensor R = R ab R ab of the space-time around the black hole, and it takes the form The GR limit has the form, in Schwarzschild limit it tends to zero, R = 0.

Kretschmann scalar
Finally, we calculate the Kretschmann scalar of the space-time around the black hole, and it has the following form: The GR limit (α = 0) of the Kretschmann scalar of the space-time around the black hole, It has the form K = 48M 2 /r 6 in the Schwarzschild case.

Particle motion
In this section, we study the dynamics of neutral test particles in the space-time of quintessential black hole in 4D EGB theory. The Lagrangian for a neutral particle with mass m The conserved quantities in the dynamics of the particles read where L and E are the special angular momentum and energy of test particles, respectively. Equations of motion for a test particle are governed by the normalization condition where is 0 and −1 for massless and massive particles, respectively. The equation of motion of the test particles can be governed by timelike geodesics of the space-time around the black hole, and it can be found using Eq. (15). The equation for the radial motion of test particles in the equatorial plane where θ = π/2,θ = 0 reads as, where the effective potential Figure 4 demonstrates the radial dependence of the effective potential for radial motion of test particles in the equatorial plane of the space-time around the quintessential black hole in 4D EGB theory for different values of the GB and quintessential field parameters. It is seen from the figure that the presence of quintessential field decreases the maximum in the values of effective potential, while positive (negative) values of the α parameter cause to increase (decrease) it and the position where the effective potential takes maximum value slightly shifts to the black hole. Moreover, the increase of ω q parameter also decreases the effective potential and the particles can have orbits closer to the central black hole.
The circular motion of test particles around a black hole can be described by the conditions: Conditions given in Eq. (18) help to find the expression for angular momentum of the particles at circular orbits, as well as corresponding to specific energy.
The radial profiles of the specific energy (in the left panel) and angular momentum (in the right panel) are shown in Fig. 5 for different values of the parameters α, q and ω q . It is seen that the presence of quintessential field decreases both the energy and angular momentum, and the effect of the GB space is essential near the compact object. Moreover, the asymptotic values of the energy and angular momentum decrease as ω = −0.5 compared to their values at ω q = −0.35. Fig. 6 The dependence of ISCO radius around the quintessential black holes in 4D EGB gravity from the quintessential parameter Fig. 7 The same relations shown in Fig.2, but for the ISCO radius being the same as for the Schwarzschild black hole

Innermost stable circular orbits
The circular orbits of test particles around compact gravitational objects can be described in standard way using the condition V eff ≥ 0. In fact, the innermost stable circular orbit (ISCO) radius is a solution of equation It is also hard, and sometimes even impossible, to solve Eq. (20) with respect to the radial coordinate. So, we can solve it numerically and analyse graphically. Below, we will show dependence of ISCO radius of test particles on the quintessential parameter q, by varying the GB and ω q parameters.
According to Fig.6, ISCO expands under the effects of quintessential field. The expansion is linear at ω q = −0.35, while at ω q = −0.5 we have two possible orbits for the same values of parameters and there is an upper limit for parameter q where two solutions of Eq. (20) coincide. Moreover, negative (positive) values of the GB parameter cause weak increase (decrease) of the ISCO radius.
We are also interested in the possibility of a quintessential black hole producing the same gravitational field the as the Schwarzschild black hole. We will analyse relations between q and α parameters providing r ISCO = 6M. In Fig.7, we have shown relations between the GB and quintessential parameters providing the radius of ISCO equal to 6M, mimicking the Schwarzschild black hole case. It is seen that the relations have linear character, which is similar to the results obtained in Fig.3.

Keplerian frequency
The angular momentum of Keplerian orbits is described by so-called Keplerian frequency, which has the form and in the space-time around a static black hole, it has the following form Now we calculate the Keplerian frequency in space-time around the quintessential black hole in 4D EGB theory using expression (22) and we get Furthermore, for the qualitative analysis to understand the values for frequencies, we will convert the geometrized unit of the frequency to Hz using the following relation: One can see from Fig. 8 that the presence of quintessential field and positive values of the GB parameter cause the Keplerian frequency near the photonsphere to be slightly lower, while negatives of α and values of ω q near −1 increase it. This can be useful to propose new tests of the proposed gravity and solution using the QPO data.

Harmonic oscillations
Here we will investigate the fundamental frequencies of the radial and vertical oscillations of test particles around the BH which are useful to explain the origin mechanisms of QPOs. The frequencies can be calculated, assuming there are radial r → r 0 + δr and vertical θ → θ 0 + δθ perturbations to the circular orbits. Using where radial and vertical angular frequencies are noted Ω r and Ω θ , respectively, and it is assumed that they are measured by a distant observer. The expressions for these frequencies have the form, After a little algebra we can have the following analytic form of the radial and the vertical frequencies in space-time around the quintessential black hole in 4D EGB gravity: where Now, we analyse the effects of parameters q, ω q and α on the radial oscillations graphically. Radial dependencies of radial frequencies are plotted on Fig.9 for the values of q = 0.01, ω q = −0.35 and −0.6, and α ± 0.5. By comparing the curves, one can easily see that at q = 0.01, ω q = −0.35 in GB frame, the frequency is less than in Schwarzschild case. In 4D EGB gravity, when α = 0.5 at the same values of the quintessential parameters, the frequency is more than in GR. Frequency decreases at negative values of α as ω q decreases.

QPO models in 4D quintessential field
In this section we study the values of upper and lower frequencies of twin peak QPOs around the quintessential black holes in 4D gravity comparing results with the Schwarzschild and Kerr black hole cases [31,33,36]. In this way, according to QPO models, we obtain radial dependencies of the frequencies, as well as the GB and quintessential parameters. Then we will provide the values of the frequencies by fixing the black hole parameters, and plot all possible values of the lower and upper frequencies at the distances from the corresponding ISCO radius to larger distances. In this work, as a simple example, we use the following wellknown model called relativistic precession (RP) model [10,32,35]. According to this model, the upper and lower frequencies are defined by the radial and orbital frequencies in the form ν U = ν φ and ν L = ν φ − ν r , respectively. But in modified RP1 and RP2 models they are ν U = ν θ , ν L = ν φ −ν r and ν U = ν φ , ν L = ν θ −ν r , respectively.
In Fig.10, we have shown the set of possible relations (values) of upper and lower frequencies of twin peaked QPOs in the RP model. Here we compare the values in the frames of GR (around the Schwarzschild and Kerr black holes) and 4D EGB gravity in the presence of quintessential field effects. In analysis of a twin peak QPO object, we can find the position of the object on the diagram and if it lies on the black line, then we can say the central black hole of the QPO is the pure Schwarzschild black hole. In the case when it's positioned over the black line (corresponding to the Schwarzschild black hole case) then the central black hole can not exactly be the Schwarzschild one. It can be a rotating Kerr black hole, or it is a 4D quintessential black hole with corresponding spin and quintessential parameters. In this case, there is a problem (not possible) to distinguish the nature of central black hole: it can be either Kerr or quintessential black hole in 4D EGB theory. Let us see another scenario where the position of another twin peak QPO lying under the black line, the source can only be a black hole in 4D EGB gravity.

Conclusions
In this work, we have analysed the properties of the space-time of the quintessential BH in 4D EGB gravitational theory and possible values for the coupling α and quintessential parameters q that provide the existence of event horizon. The case when the quintessential black hole event horizon lies at 2M, the same as the Schwarzschild black hole horizon radius with corresponding parameters, is also discussed. The scalar invariants of the spacetime have also been studied and their GR limits are shown. We also have studied a test particle dynamics in the space-time of a quintessential black hole in a novel 4D Einstein-Gauss-Bonnet theory. The effects of the quintessential and GB parameters on the effective potential for radial motion, specific energy and angular momentum of the particles corresponding to their circular orbits have been investigated. The effects of the parameters on ISCO radius have also been studied, and the relations between the quintessential and GB parameters providing 6M ISCO radius are shown, reflecting the same gravitational field effects on ISCO as the Schwarzschild black hole. Moreover, fundamental frequencies such as Keplerian and harmonic oscillations are also studied together with applications to quasi-periodic oscillations (QPOs) where we discussed the way to test gravity theories through astrophysical observations of twin-peak QPOs.