Quasiperiodic oscillation around charged black holes in Einstein–Maxwell-scalar theory

In the present paper, first, we study the event horizon properties of charged black holes (BHs) in Einstein Maxwell-scalar (EMS) gravity. Then, we investigate the circular motion of test particles’ around the BH in the EMS gravity. We also analyze the effects of the EMS parameters on the position of innermost circular orbits (ISCOs), energy, and angular momentum of the test particles corresponding to circular orbits. We provide detailed studies of the efficiency of energy release from EMS BHs based on the Hartle–Thorne model and fundamental frequencies of oscillations of particles along their circular stable orbits. Moreover, we have explored possible values of upper and lower frequencies of twin-peak quasiperiodic oscillations (QPOs) around the BHs. Finally, we obtain relationships between the BH charge and the EMS parameters using observational data from the QPOs detected in the microquasars: GRS 1905+105, GRO J 1655-40, H 1745+322, and XTE 1550-564.


Introduction
The standard theory of gravitational field known as general relativity (GR) was introduced by Albert Einstein in 1915. GR as gravitational field theory has been successfully tested in both weak and strong field regimes [1][2][3][4]. However, GR meets several fundamental problems which cannot be resolved within the basic principles of the theory. One of the above-mentioned problems is the existence of the singularity at the origin of most general relativistic solutions. On the other hand, the current resolutions of experiments and observations leave open the window for the development of alternative/modified theories of gravity, which can be further used to resolve the fundamental open issues of gravitational interaction.
One of the possible ways of extension of GR is to consider the low energy limit of string theory, where one may introduce the dilaton scalar field. This field will be included in the action as an additional supplement term to the Einstein-Hilbert action and has the form of the axion, gauge fields, and another non-trivial coupling of dilaton to fields. Particularly, some aspects of such combination like the causal structures and thermodynamic properties of black hole (BH) solution with dilaton have been investigated in Refs. [5][6][7][8][9][10][11][12][13]. Further, one may extend the black hole solution using the cosmological constant. The models with the presence of negative cosmological constants are prominent candidates for grav-ity theories defined in higher dimensions. Detailed analysis of these models and corresponding solutions can be found in Refs. [14][15][16][17][18][19]. One may also consider the heterotic string theory, where the scalar dilaton field is coupled to the electromagnetic field tensor [6]. Here, we plan to study the properties of the black hole solution within the Einstein-Maxwellscalar (EMS) field theory using the dynamics of test particles.
The tests of the field theories based on the observations of compact objects allow us to get the constraints on the field parameters. Particularly, the modified and alternative theories of gravity and corresponding solutions describing the compact objects have been successfully tested using X-ray data from astrophysical objects by the Authors of Refs. [20][21][22]. One may also develop the tests of the general relativity and other models of gravitational interaction based on the test particle motion [23,24].
One of the interesting features of the dynamics of test particles appears with the analysis of circular/bound orbits around compact objects. Particularly, circular orbits of test particles may lead to the phenomenon observed as quasiperiodic oscillations (QPOs) in astrophysics. QPOs are objects observed in the X-ray band from microquasars representing BH surrounded by matter flowing and a companion star. QPO has been first detected by analyzing the power spectra of the flux from X-ray binary pulsars [25]. After that pioneering discovery, a big number of works have been dedicated to studying of QPOs and corresponding models. Among a big number of models proposed to describe the behavior of the origin and nature of QPOs those models based on particle motion are the most promising. According to these models, the innermost stable circular orbits of test particles determine the frequencies of the QPOs. Thus, the investigations of QPOs using models based on circular orbits of test particles may lead to the test of gravity theories in the strong field regime. Particularly, the current high resolution of the observations of QPOs frequencies allows us to test the models of gravity and corresponding solutions and get the constraints on field parameters. Depending on the value of the frequencies, one may distinguish two types of QPOs: low frequency (up to 30 Hz, LF QPOs) and high frequency (up to 500 Hz, HF QPOs). A big number of works are devoted to studying the models where QPOs are considered as a result of the collective motion of particles in the accretion disk of different type [26][27][28][29][30][31][32][33][34][35][36][37][38][39][40][41][42]. These models can be further used to develop a new test of EMS theory, and here we plan to explore this possibility.
Here we plan to study the spacetime properties as well as neutral particle dynamics around BH within EMS gravity. The paper has organized as follows: Sect. 2 is devoted to studying the spacetime properties of curvature scalars. Then in Sect. 3 we investigate the massive test particle motion around the compact object. We discuss the fundamental frequencies and their application to QPOs in Sect. 4. In Sect. 5 we study the quasi-periodic oscillations of circular orbits of the particles around BH in EMS gravity. Using the observation data on QPO objects, we obtained constraints on parameters of EMS gravity in Sect. 6. We conclude and summarize the obtained results of the paper in Sect. 7. Our study in this paper uses a space-like signature (−, +, +, +) for the spacetime and geometrized unit of the system where G = c = 1. Latin (Greek) indices run from 1 (0) to 3.

Black holes in Einstein-Maxwell-scalar theory
We start with a brief review of BH solutions in EMS theory described by the action [5,6,43] where ∇ α is the covariant derivative, g is the determinant of the metric tensor g μν , R is the Ricci scalar of the curvature, φ is the massless scalar field, F αβ is the electromagnetic field tensor, K (φ) is the coupling function between the dilaton and the electromagnetic fields. The exact analytical solutions for spacetime around the static black hole in general form have been found in Ref. [43] in the following form, where U (r ) and f (r ) are radial functions. In Ref. [43] Authors have used the following special forms for the function K (φ) and obtained the solution in the form Mr , It is worth noting that it has been assumed that the vector and the dilaton fields depend on the radial coordinate only as [43] The spacetime metric (1) with the metric functions (3) turns to Schwarzschild solution when β = 0 and/or Q = 0, also it covers the Reissner-Nordström (RN) one when γ = 0 and β = 1. Now, we will analyze the event horizon structure of the spacetime (1) defined by the metric functions (3) using the One may set the following simple conditions with the aim of detailed analyses of the minimum radius of the outer horizon and extreme values of BH charge. In Fig. 2

The motion of test particles around BH in EMS theory
The equation of motion of test particles with the rest mass m around a BH can be found using the Lagrangian of from The constants of motion read, where E and L are the specific energy and angular momentum of the particle, respectively. Equations of motion for a test particle are governed by the normalization condition where takes the values 0 and −1 for massless and massive particles, respectively. For test neutral particles with non-zero rest mass, the equation of motion can be governed by time-like geodesics of spacetime. Using the condition (9) together with the Eq. (8) one can easily obtain the equation of motion in the following form [44][45][46] where K is the Carter constant corresponding to the total angular momentum of the particle.
Restricting the motion of the particle to a constant plane with θ = const, andθ = 0, and the Carter constant reads as K = L 2 / sin 2 θ . Consequently, one can easily obtain the equation of the radial motion in the standard forṁ where the effective potential of the motion of test particles has the following form The radial dependence of the effective potential for the radial motion of test particles in the spacetime around charged EMS BHs is shown in Fig. 3 for the different values of the parameters of EMS theory β and γ . In this figure, we fix the BH charge and angular momentum of the particles as Q = M and L = 4.5M, respectively, and consider the equatorial plane where θ = π/2. One can see from Fig. 3 that the maximum value of the effective potential decreases with the increase of γ at β = 1 case. While the increase of β causes the increase of the maximum in the effective potential. This leads to the fact that with the increase of β the radius of circular orbits where the effective potential takes maximum shifts towards the BH. Now we consider the circular motion of test particles at the equatorial plane (θ = π/2). Circular orbits of the particles can be determined by the following conditions: Using the conditions (16) one may easily obtain the expressions for specific angular momentum and energy of the particle corresponding to the circular obits around the charged BHs in EMS gravity in the following form where The radial dependence of specific energy and specific angular momentum of the test particles at circular orbits around charged BHs in EMS theory are shown in Fig. 4. Here, we have fixed the BH charge as Q/M = 1/2. From Fig. 4 one can see that the presence of the BH charge decreases the minimum value of both energy and angular momentum. However, the minimum energy and angular momentum increase with the increase of β. On the other hand, the increase of the parameter γ causes an increase in the minimum values of energy and angular momentum. Figure 5 shows allowed values of the energy and angular momentum for bounded orbits around EMS BHs for the different values of the parameters β and γ and for the fixed value of the BH charge Q = M. One may see from the figure that the particle energy may lie between its maximum and minimum (E min ≤ E ≤ E max ) for a fixed value of the angular momentum at L > L cr , while at L = L cr the maximum and minimum energies equals to each other and takes its critic value which is shown as black dots. The critical values of both energy and angular momentum decrease with the increase of β for the fixed values of γ parameters, while it increases with respect to the increase of γ for fixed β.

Innermost stable circular orbits
The radius of innermost stable circular orbits (ISCO) of test particles can be defined using an additional condition V eff ≥ 0. Using this condition and Eqs. (17) and (18) one may obtain the equation for ISCOs of the form One may get an analytical solution corresponding to the case when γ = 0 as where We have obtained the numerical solution of Eq. (19) for the case of when γ = 0. The dependence of the BH charge on corresponding values of ISCO radius is shown in Fig. 6 for different values of the EMS parameters β and γ . One can easily see from this dependence that at Q = 0 i.e. Schwarzschild limit the ISCO radius equals to 6M and the increase of the BH Fig. 6 Dependence of ISCO radius from the BH charge for different values of EMS parameters β and γ charge forces to decrease the ISCO radius, and it decreases rapidly with increasing β. While the presence of the parameter γ causes the increase of the radius for the fixed values of Q (Q < Q extr ).
In Fig. 7 we have shown relationships between energy (right panel) and angular momentum (left panel) of test particles at ISCO around EMS charged BHs for the different values of EMS parameters β and γ . One can see from the figure that the increase of γ causes the increase of both energy and angular momentum corresponding to ISCOs, while β affects oppositely.
Relationships between energy and angular momentum of test particles at ISCO around charged BHs in EMS gravity, shown in Fig. 8 for the various values of β and γ . Similar effects can be observed as it is seen in Fig. 7.

Energy efficiency
The energy of test particles along Keplerian orbits in the accretion disk around BHs takes its minimum value. According to the thin accretion disk mode, the energy depends on  Relationships between E ISCO and L ISCO for different values of parameters β and γ at the range of the BH charge from 0 to corresponding Q extr the properties of the spacetime circular geodesics [47]. The energy efficiency of the spacetime around BHs indicates the maximum extracted energy of the radiation of infalling particles into the BH, and it can be calculated as the difference of their rest energy and orbital energy at corresponding ISCOs (E ISCO ) using the following expression where E ISCO characterizes the ratio of the binding energy. The dependence of efficiency of energy release from EMS BHs from the BH charge Q is shown in Fig. 9 for various values of the EMS parameters β and γ . Here, we have compared the results obtained in EMS gravity with the results for the RN BH case. One can see from the figure that the efficiency increases as the BH charge increases, up to about 16%.

Fundamental frequencies
In this section, we study the fundamental frequencies of test particles orbiting the charged BH in EMS gravity. Particularly, we analyze the frequency of Keplerian orbits and frequencies of oscillations of the particles along radial and vertical directions with respect to the circular orbits.

Keplerian frequency
The angular velocity Ω K =φ/ṫ of test particles in circular orbits or so-called Keplerian orbits around static BHs is defined as After mathematical calculations, the Keplerian frequency in spacetime (1) with metric functions (3) takes the following In order to estimate the value of the fundamental frequencies one may express them in the unit of Hz multiplying the frequencies by the factor c 3 /(2π G M). Here, we use the values of the speed of light at vacuum as c = 3 · 10 8 m/s as well as the gravitational constant as G = 6.67 · 10 −11 m 3 /(kg 2 · s).
The radial profiles of the Keplerian frequencies of test particles around charges EMS BHs are presented in Fig. 10 for different values β and γ parameters and fixed value of Q/M = 1. One can see from the figure that the Keplerian frequency decreases with the increase of the parameter γ when β = 1. It is observed from the comparison of bluedashed and red large-dashed lines that an increase of the parameter β causes an increase in the frequency.

Harmonic oscillations
The small perturbation of radial (r → r 0 + δr ) and vertical (θ → θ 0 + δθ) coordinates of the test particle, in their stable circular orbits at the equatorial plane around the BH, cause to oscillate the particles along these axes.
One can find the equation of the oscillations by expanding the effective potential in terms of the coordinates, r and θ , and using the condition for the extremes of the effective potential V eff (r 0 , θ 0 ) = 0 and ∂ r (θ) V eff = 0 as where are the square of radial and vertical angular frequencies of the particles around the BH measured by a distant observer, respectively. After small algebraic calculations, one easily gets expressions for the radial and vertical frequencies as, where Ω φ is the angular velocity of the particle measured by an observer at infinity. Using the Eq. (26) and the metric function (3) one can easily find the corresponding fundamental frequencies. Figure 11 shows frequencies of the radial oscillations of test particles around charged BH in EMS theory for various values of EMS parameters β and γ for the fixed value of the BH charge as Q/M = 1. One can see from the figure that the frequency increases in the presence of the BH charge. Moreover, for the cases, when β = 1/2 and γ = 1 the maximum in the radial frequencies increases, while at β > 1 and γ = 1 the frequency increases sufficiently. The increase of both γ and β parameters causes a shift in the position of orbits, where the frequency of radial oscillations reaches its maximum, towards the central BH.

QPOs around BHs in EMS theory
In this section, we focus our attention on studies of twin peak QPOs around charged BHs in EMS gravity in the frame of relativistic precession (RP) model [48]. The upper and lower frequencies of the QPOs using combinations of the frequencies of the radial, vertical and orbital oscillations of the test particles along the stable circular orbits around the BH as ν U = ν φ and ν L = ν φ − ν r , respectively.
In Fig. 12 we plot relationships between the upper and lower frequencies of twin-peak QPOs from charged BHs in EMS gravity in the RP model. Here, we fix the BH mass as M = 15M and the charge as Q = M. It is observed from this figure that the QPO frequencies grow in the presence of the BH charge when β = 1, while in cases of β < 1 the QPO can be observed near the LF QPOs. Similarly, the increase of γ also causes decreasing frequencies at β = 1. However, in β > 1 cases the frequencies take value near kHz QPOs. We have shown black dots where the black solid line and other colored lines cross, in zoom on the top of the figure zoomed. It does mean that at a fixed value of the BH charge and EMS parameters, the upper and lower frequencies can be the same as it is in the pure Schwarzschild case, but the QPO may occur at different distances from the central BH. However, the colored lines do not cross each other, implying there are no degeneracy values in the EMS parameters. Figure 13 demonstrates the dependencies of the radius of QPO orbits with the frequency ratio 3:2 and 5:4 for various values of the EMS parameters. It is observed from the figure that the QPO orbits located out of ISCO and QPOs with a ratio near 1 shine at the orbits close to ISCOs. It is known that ISCO radius is one of the most critical properties of BHs to be measured. It is still a problematic issue in observations of stellar mass and intermittent BHs to solve. From this point of view, the studies may play a role in being a solution to this  issue. It implies that one may calculate the ISCO radius using frequency data from twin peak QPOs with a frequency ratio less than 5:4. Moreover, the QPO radius comes closer to the ISCO at higher values of the BH charge. On the other hand, it is possible to get constraints and relations between the EMS parameters and BH charge using the results of studies of QPO orbits.

Constraints on EMS gravity parameters using QPO frequencies
Here we try to get constraints on the EMS parameters β and γ and the BH charge around EMS BHs using data from QPOs in the following microquasars:  On the other hand, the BHs' masses are not also too much different from each other. Our expectation is that their charge and the EMS parameters for these microquasars must take values that are close to each other.
Since there are two independent parameters in EMS theory including the BH charge, it is not possible to get exact values for the parameters using the above-mentioned observational data from the objects. However, one may get relationships between the BH charge and the parameter γ by keeping the parameter β as a constant with the allowed values. Now, we try to get the relationships on the EMS parameter γ and the BH charge Q for the fixed values of the parameter β solving the following equations, in the frame of RP model [40,54,55], where ν ob U and ν ob L are observational values of the upper and lower frequencies of the QPOs. In order to get relationships between the possible values of the EMS parameters for the above-mentioned objects, we use their observational parameters. Below, we provide the relationships for the abovementioned objects.
In Fig. 15 we provide dependence of charge of central BHs in the microquasars GRS 1905+105, GRO J 1655-40, H 1745+322, and XTE 1550-564 from the EMS parameter γ for the fixed values of β. One can see from the figure that the BH charge increases linearly with the increase of the γ parameter and the increase of β causes to decrease in the charge. It is seen that the dependencies for each value of β are close to each other because the frequency ratios of the QPO objects are close to 3:2.
In order to show a clear picture of the difference in the relationships, we provide the sectional view of the dependence in the range γ ∈ (0.1 − 0.15) in Fig. 16, while, in order to see the charge of the RN BHs, we plot the top left panel from γ = 0 to 0.1. One can see from the panel that the charge of the BHs in the microquasars take the values between ∼ 1.2 M -1.215M , and the BH charges increase with the increase of the parameters γ and β. Moreover, it is observed that the difference between the lines, corresponding to BH charge profiles, increases with the increase of β.

Conclusion
In the present work, we have studied profiles of event horizon radius with respect to the variation of the EMS gravity parameters and the BH charge. It is found that when γ → −∞ the extreme value of the BH charge tends to zero and the spacetime (3) turns to the pure Schwarzschild one. It has been shown that the γ parameter has an upper value that depends on the parameter β (see Fig. 2). Relationships between extreme charge, maximum in the values of γ and β parameters have obtained as γ ∼ β and Q extr ∼ β − 1 2 . We have also investigated the dynamics of test particles in the spacetime of charged BHs in EMS theory. Analytic expressions for the specific angular momentum and energy of test particles at circular orbits are found, and the effects of the parameters β and γ on them are discussed. It has been shown that an increase of the BH charge causes decreasing the minimum of both energy and angular momentum, however, the presence of γ and β parameters increases.
We have analyzed EMS field effects on ISCO radius, energy, and angular momentum of particles at ISCOs, as well as the energy efficiency. It has been obtained that ISCO radius (energy efficiency) rapidly decreases (increases) with the increase of the BH charge Q with compare to RN BH case (and reaches up to about 16%).
Moreover, we have explored fundamental frequencies, such as Keplerian and radial frequencies of oscillations of particles orbiting the EMS BHs. As an astrophysical application of the studies, we have developed the RP model for twin peak QPOs.
It is well known that one of the problematic issues in observations of stellar mass and intermittent BHs is ISCO radius measurements, being an essential property of BHs. From this point of view, the study plays an important role in solving the problem. We have shown that ISCO radius can be measured by using frequency data from twin peak QPOs with a frequency ratio less than 5:4 with high accuracy. Moreover, the QPO radius comes closer to the ISCO at higher values of the BH charge. The frequency ratio becomes 1 when it is generated at ISCO.
Finally, we have got constraints on the BH charge and EMS parameters for observed QPOs in the microquasars GRS 1905+105, GRO J 1655-40, H 1745+322, and XTE 1550-564 which have frequency ratios close to each other and 3:2 with the stellar mass central BH.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: The present work is pure theoretical, and we have not used any observational data.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 . SCOAP 3 supports the goals of the International Year of Basic Sciences for Sustainable Development.