Testing alternative theories of gravity by fitting the hot-spot data of Sgr A*

We are fitting dynamics of electrically neutral hot-spot orbiting around Sagittarius A* (Sgr A*) source in Galactic center, represented by various modifications of the standard Kerr black hole (BH), to the three flares observed by the GRAVITY instrument on May 27, July 22, July 28, 2018. We consider stationary, axisymmetric, and asymptotically flat spacetimes describing charged BHs in general relativity (GR) combined with nonlinear electrodynamics, or reflecting the influence of dark matter (DM), or in so called parameterized dirty Kerr spacetimes, and test them using the hot-spot data. We show that the orbital frequencies as well as positions of the hot-spots orbiting the considered BHs fit the observed positions and periods of the flare orbits, and give relevant constraints on the parameters of the considered BH spacetimes and the gravity or other theories behind such modified spacetimes.


INTRODUCTION
General Relativity is an elegant theory of gravity that agrees with all observations at Solar System scale and beyond (Will 2001).The most famous results of GR are: the bending of light due to the gravitational field (Will 2015), the correction of precession of perihelion of Mercury (Will 2018), and the existence of gravitational waves (Abbott et al. 2016).Another important prediction of GR is the existence of astrophysical objects (BHs) with strong gravitational interaction from where nothing can escape, even the light.However, the nonlinear behaviour and strong-field structure of GR still remains elusive and difficult to test (Psaltis 2008).
Observational data based on the dynamics of whole Universe affirm that the major part of mass of the universe is invisible (in modern cosmology, this invisible mass is known as DM).Even larger part of invisible content of the Universe is related to the so called dark energy (DE) (Caldwell & Kamionkowski 2009) that can be well represented by the cosmological constant.For its relevance in astrophysical processes, see (Stuchlík & Hledík 1999;Stuchlík 2005;Slaný & Stuchlík 2005;Balaguera-Antolínez et al. 2007;Stuchlík et al. 2020).
Modern cosmological observations also reveal that our Universe is composed of 68.3% DE, 26.8% DM, and 4.9% ordinary matter (Planck Collaboration et al. 2014;Rezaei 2017).Dark matter surrounding the galaxies and clusters does not interact with baryonic matter but can be observed by its gravitational effects on visible matter.Babcock (Babcock 1939) examined the rotational speed of luminous objects in Andromeda galaxy and found that rotational speed continuously increases as one moves away from the center of these objects.This demonstrates that outer region of that luminous part is dominated by matter which does not shine.Zwicky (Zwicky 2009) found a large amount of unseen (nonluminous) matter in the Universe rather than the seen (luminous) and detected the non-luminous mass lying outside the luminous parts of the galaxies.Besides of these theoretical observations, there is no experimental success in detecting DM yet.
In addition to the need for DE and DM, in the study of BHs, a problem that appears in GR, is the presence of singularities that are points or set of points where the geodesic is interrupted and the physical quantities diverge (Hawking & Penrose 1970;Bronnikov & Rubin 2013).It is believed that the problem of singularity occurs because the theory is classical and that in the quantum theory of gravity this problem would be solved.This, together with some long-standing problems in GR (like difficulties in explaining the accelerated Universe and galaxy rotation curves, etc), has motivated the study of viable alternative theories of gravity.These theories, also known as modified theories of gravity, aim to reproduce GR in the weak-field regime, but they can differ substantially from it in the strong curvature regime, where non-linear effects become dominant.These modified theories of gravity are developed by modifying the matter or gravitational part of the Einstein-Hilbert action.
Astronomers have classified the astrophysical BH candidates into three major classes (depending on the mass of BHs): stellar-mass BHs with mass M ∼ 5M −20M situated in X-ray binary systems; super-massive BHs having M ∼ 10 5 M − 10 9 M found in galactic nuclei; and intermediate-mass BHs with M ∼ 10 2 M − 10 4 M (Narayan 2005).The third class of objects (intermediate-mass BHs) is still debatable because their observations are indirect and dynamical measurements of their masses are still lacking.
The Sagittarius A* (Sgr A*) is very compact and bright astronomical radio source situated at the Galactic Center of Milky Way, associated with the supermassive BH, and considered as a highly variable source across all the wavelengths (Eckart et al. 2005;Melia 2007).A precise measurement of its spin and mass is a long standing problem for astrophysicists.The mass of the supermassive BH Sgr A* estimated by the observed orbital parameters of the S0 star traveling in the BH's gravitational field with velocity of ∼ 10 3 km/s is M = (4.1 ± 0.4) × 10 6 M (Gillessen et al. 2009).The mass and spin of Sgr A* have been estimated using different approaches in several studies.(Dokuchaev 2014(Dokuchaev , 2016) ) has been estimated the values of mass M = (4.2± 0.2) × 10 6 M , and spin a = 0.65 ± 0.05 with the help of the observed quasi-periodic oscillations from the supermassive BH in the Galactic center in Xrays.However, the observed high-frequency QPOs could indicate a much higher spin on a ∼ 0.999 and relevance of the so called Aschenbach effect (Aschenbach 2004;Stuchlík et al. 2005) Near-infrared GRAVITY@ESO observations (Gravity Collaboration et al. 2018) have revealed the detection of two bright flares on July 22 and July 28, 2018 as well as a fainter flare on May 27, 2018, in the background of Galactic center massive BH.These flares are found with a peak approaching the flux of S2, and remain for 30-90 min.The GRAVITY observations of hot-spots nearby the ISCO of Sgr A* BH allow us to fit their orbital radii as well as orbital periods with the circular orbits of test particles orbiting Kerr BH with mass M ∼ 4 million M .
The dynamics of charged hot-spot around supermassive Kerr BH in the presence of magnetic field as well as the equatorial circular orbits fitting the observed periods and positions of three flares have been explored by Tursunov et al. (2020).Recently, (Matsumoto et al. 2020) investigated the July 22 flare in the background of several different models including geodesics, circular Keplerian, precessing pattern, and super-Keplerian pattern (a hot-spot moving faster than Keplerian along a circular trajectory), and fit the hot-spot trajectories.They proposed that a super-Keplerian circular orbits with orbital frequency Ω = 2.7 Ω k at r = 12.5M yields a better match to the data than the Keplerian orbits, where Ω k is the orbital frequency of Keplerian orbits.A model for the flares formulated on general relativistic magnetohydrodynamic simulations of magnetically arrested accretion disks which present the violent episodes of flux escape from BH magnetosphere has been discussed by Porth et al. (2021).
The nature and origin of hot-spots or flares still remains unclear.In addition, there is no observational information regarding shape of hot-spot.We assume that the hot-spot is a bound test, potentially neutral mass moving on a circular orbit in the background of Kerr spacetime, or in some of different BH spacetimes.We also assume that shape of the hot-spot as well as the surrounding environment does not change during one orbital timescale.
In this study, we consider the classical Kerr BHs, rotating regular BHs in GR (regular Bardeen, regular ABG), a well known Johannsen-Psaltis spacetime (non-Kerr BHs), and various other metrics in several theories of gravity, i.e., Kerr-Sen BHs in heterotic string theory, Born-Infeld BHs in Einstein-Born-Infeld theory, Kalb-Ramond BHs in heterotic string theory, Gauss-Bonnet BHs in Einstein-Gauss-Bonnet theory, Konoplya-Zhidenko BHs in an unknown alternative theory of gravity, rotating BHs in perfect fluid DM, as well as rotating regular BHs in Einstein-Yang-Mills theory, and explore the dynamics of neutral hot-spot orbiting the considering BHs.We calculate the corresponding orbital frequency, and examine the equatorial circular orbits of neutral hot-spot fitting the observed positions and periods of the three flares proposed by GRAVITY on July 22, May 27, and Jul 28, 2018.From dependence of the fitting lines on the parameters of various BH spacetimes we immediately indicate spacetime promising good fits for the data of the three flares.
Throughout the paper, we use the space-like signature (−, +, +, +) and the system of units in which c = 1 and G = 1.However, for the expressions with an astrophysical application and estimates, we use the units with the gravitational constant and the speed of light.Greek indices are taken to run from 0 to 3; Latin indices are related to the space components of the corresponding equations.

STATIONARY AND AXI-SYMMETRIC SPACETIMES
In four-dimensional GR, the no-hair theorem (Israel 1967;Carter 1971) states that the uncharged rotating BHs are uniquely characterized by only two parameters, the mass M , and spin a of the BH, and are governed by the Kerr metric.This metric is a unique axisymmetric, stationary, asymptotically flat, and vacuum solution of the Einstein field equations which possesses an event horizon but there is no closed timelike curves in an exterior domain.Due to the weak cosmic censorship conjecture (Penrose 1969), the central singularity is always behind the event horizon.However, the hypothesis that the astrophysical BH candidates are characterized by the Kerr spacetimes still lacks the direct evidence, furthermore, the GR has been tested only in the regime of weak gravity (Will 2014).For strong gravitational fields, the GR could be broken down and astrophysical BHs might not be the Kerr BHs as predicted by the no-hair theorem (Johannsen & Psaltis 2011).Several parametric deviations from the Kerr metric have been proposed to investigate the observational signatures in both the electromagnetic and gravitational-wave spectral that differ from the expected Kerr signals.
The line element of an arbitrary, stationary, axisymmetric, and asymptotically flat spacetime with refelction symmetry reads where metric components g αβ are functions of r, θ, and some additional parameters.In the following, we consider several stationary, axisymmetric, and asymptotically flat spacetimes both in GR and modified theories of gravity.

Classical BHs (Kerr) in GR
The nonzero components of the metric tensor g µν , describing the geometry of the well known classical neu-tral rotating Kerr BH, taking in the standard Boyer-Lindquist coordinates can be written in the form (Kerr 1963;Carter 1968) with where, M and a are the mass and rotation parameter of the BH, respectively.The spin parameter a is bounded by a ≤ M .The horizons for Kerr BH can be found by solving the condition ∆ = 0.

Charged BHs in GR
According to the no-hair theorem, BH solutions of the Einstein-Maxwell equations of GR (combining the field equations of gravity and electromagnetism) are fully characterised by their mass M , rotation parameter a, and electric charge Q.There are many kinds of charges such as electric, magnetic, tidal, and dyonic, etc.In the following we consider BH solutions with different charges.

Kerr-Newmann BHs
The non-zero components of metric co-effiecients of Kerr-Newman (KN) BH takes the form (Misner et al. 1973;Bicak et al. 1989) where Q is the electric charge of KN BH.For vanishing charge ( Q = 0), the KN BHs reduces to the Kerr BH solutions.

Braneworld BHs
Rotating charged BHs in the brany universe of the Randall-Sundrum type with infinite additional dimension are described by the Kerr geometry with an additional parameter, represented by the line element with following metric co-efficients (Aliev & Gümrükçüoǧlu 2005;Kotrlová et al. 2008;Stuchlík & Kotrlová 2009) with The tidal charge parameter β represents the interaction between brany BH and bulk spacetime, and can be both negative and positive.The negative tidal charge can provide a mechanism for spinning up the BH so that its rotation parameter exceeds its mass, which is not allowed in the framework of GR.

Dyonic charged BHs
A particle having both electric and magnetic charge is called dyon.The possibility of existence of dyonic BHs is either due to magnetic monopoles raised into grand unification theories or it may be primordial.The metric co-efficients for rotating charged dyonic BHs are given by (Kasuya 1982;Stuchlik 1983) with where Q e is electric charge, and Q m is the magnetic charge.

Kerr-Taub-Nut BHs
The Kerr-Taub-Nut (Newman-unti-tamburino) solution is an analytic type D vacuum solution of the Einstein equations and can be represented by (Demianski & Newman 1966;Miller 1973) where n is the gravitomagnetic or Nut charge, and for limiting value n = 0, Eq. ( 11) reduces to the usual Kerr BH solution.The Kerr-Taub-Nut spacetime is asymptotically non-flat due to the Nut charge, and there are string singularities on the symmetric axis.

KN-Taub-Nut BHs
The KN-Taub-Nut BH is stationary and axisymmetric non-vacuum object, completely described by mass, rotation, an electric charge, and nut parameter, written in the form (Demianski & Newman 1966;Miller 1973) with where Q n is the electric charge, and ñ is the Nut charge parameter of KN-Taubt-Nut BH.For Q n = 0, Eq. ( 13) reduces to the Kerr-Nut BH, and Q n = 0 = ñ leads to the Kerr BH.

Bumpy spacetimes in GR
It is possible that the spacetime around massive compact objects which are assumed to be BH is not described by the Kerr metric, but by a metric which can be considered as a perturbation of the Kerr metric, and are usually known as bumpy (non-Kerr) spacetimes (Collins & Hughes 2004).These spacetimes have multipoles and possesses some features that deviate slightly from the Kerr spacetime, reducing to the classical Kerr BH solutions when the deviation is zero.Here, we consider some bumpy spacetimes in GR.

Johannsen-Psaltis spacetime
In order to test the gravity in the region of strong gravitational field, Johannsen and Psaltis (Johannsen & Psaltis 2011) proposed a deformed Kerr-like metric which describes the geometry of a stationary, axisymmetric, and asymptotically flat vacuum spacetime, and the corresponding non-zero metric co-efficients can be written in the form where The deformation parameter determines the degree of variation that the BH is more oblate ( < 0), or prolate ( > 0) than the Kerr BH, and one can restore the Kerr metric for limiting case = 0.The Johannsen-Psaltis metric is a perturbation of the Kerr metric designed to avoid pathologies like naked singularities and closed timelike curves.

Hartle-Thorne spacetime
The other solution that deals with quadrupole to linear, and rotation to second order is the Hartle-Thorne metric.It is an approximate solution of GR equations and can describe an inner source for the compact object.The corresponding metric co-efficients read (Hartle & Thorne 1968) with and where, q 1 is the quadrupole parameter.The expressions for the Legendre functions , and P 2 (y) can be found in (Hartle & Thorne 1968).

Kerr-Q spacetime
The Kerr-Q metric is the simplest extension of Kerr metric that admit a quadruple and has been tested to be singularity free outside the horizon.It is derived from the rotating δ-metric using δ = (1+q 0 ), m = M/(1+q 0 ) and expanding to the first order in q 0 , second order in rotation parameter a, where M is the physical mass.The corresponding metric co-efficients read (Allahyari et al. 2020) where Here, q 0 is the quadrupole parameter that determines the deviations from the Kerr BH, and for q 0 = 0, Eq. ( 20) reduces to the Kerr BH.

Quasi-Kerr BHs
The general stationary axisymmetric neutral compact object can bee characterized by mass, multipole moments, and rotational parameter.The multipole moments are consist of a set of mass multipole moment M k and current multipole moment S k , here the subscript k of them is labeled by the angular inter eigenvalue k ≥ 0. The relation between the parameters of the multipole moments can be expressed as For classical Kerr BHs, the deviation δM k and δS k vanishes, while the BH solutions with the quadrupole moment takes the form where g Kerr αβ indicates the metric tensor for classical Kerr BHs, and the components of h αβ can be written as (Glampedakis & Babak 2006) with the functions F 1,2 (r) shown explicitly in Appendix A of (Glampedakis & Babak 2006).The deformation parameter ˜ indicates a small contribution to the quadrupole moment q of the compact object with the total mass M as and can take both negative and positive values.For vanishing ˜ = 0, Quasi-Kerr BHs reduce to the classical Kerr BH solutions.

Accelerating and rotating BHs
The accelerating and rotating BH solutions describe the gravitational field by a pair of uniformly accelerating Kerr-type BHs, which is a special case of the Plebański and Demiański spacetime that covered a large family of electro-vacuum type-D spacetimes including both the KN like solutions and the C-metric, the corresponding metric co-efficients read (Griffiths & Podolský 2005) with The parameter b determines the acceleration of the BH.The accelerating and rotating BHs have the same event and Cauchy horizons as the Kerr BH, but there also exits two other horizons which can be interpreted as the acceleration horizons, i.e., Thus, unlike in the usual Kerr BH spacetime, the physical region of this BH is situated in r + < r < r A , where ∆ A > 0 is satisfied.

Rotating regular BHs in GR
The regular BHs are non-singular exact solutions of the Einstein field equations minimally coupled to a nonlinear electrodynamics, satisfy the weak energy condition, and yield alteration to the classical BHs (Stuchlík & Schee 2015).The regular BHs are constructed to be regular everywhere, i.e., the Ricci scalar, and the components of the Riemann tensor are finite ∀ r ≥ 0. There are three well known regular BHs in GR, Bardeen rotating regular BHs, and Ayón-Beato-Garica (ABG) regular BHs, and Hayward regular BHs (Hayward 2006;Becerril et al. 2021), which we describe below.

Regular Bardeen BHs
The spacetime filled with a vacuum can give a proper discrimination at the final stage of gravitational collapse, replacing the future singularity.Based on this idea, Bardeen (Bardeen 1968) proposed the first regular BH solution named as Bardeen regular BH, according to whom there is no singularity but horizons can exist (Toshmatov et al. 2014;Stuchlík & Schee 2019).The matter field is a kind of magnetic field and the solution yields a modification of the classical Kerr BH solution.The non-zero metric co-efficients corresponding to the Bardeen regular BH read where Σ is defined by (4) and and the mass function m(r) takes the form The deviation parameter q can be recognized as a magnetic monopole charge of non-linear electrodynamics, which determines the deviation from the Kerr BH, and when we turn-off the non-linear electrodynamics (q = 0), one can recover the Kerr metric, and for vanishing spin (a = 0), we obtain the non-rotating Bardeen regular BH.

Regular ABG BHs
Another class of spherically symmetric static regular BH solutions was introduced by Ayón-Beato and García (Ayón-Beato & García 1998), and the rotating one is investigated by (Toshmatov et al. 2017a).The non-zero metric co-efficients of rotating ABG regular BH can be written in the form with where Q is the electric charge of the regular ABG BH, for vanishing charge (Q = 0), we obtain the Kerr metric, and for a = 0, Eq. ( 32) reduces to the case of nonrotating ABG BH.

Regular Hayward BHs
The regular Hayward BH solutions can be described by Eq. ( 29) with the mass function (Hayward 2006) where g is the deviation parameter, and for limiting case g → 0, one can recover the classical Kerr BH.

Rotating BHs in alternative theories of gravity
The late-time acceleration of the Universe is surely the most challenging problem in cosmology.Many cosmological observations indicate that the accelerated expansion of the Universe is due to the existence of mysterious form of energy known as DE.Modern astrophysical and cosmological models are also faced with two severe theoretical problems, that can be summarized as the DM (non or weakly interacting), and the DE problem.The two observations, namely, the mass discrepancy in galactic clusters, and the behavior of the galactic rotation curves, suggest the existence of a DM at galactic and extra-galactic scales.Recently, several modified theories of gravity have been proposed to address these two intriguing and exciting problems facing modern physics.These modified theories of gravity are constructed by modifying the gravitational or matter part of the Einstein-Hilbert action.In addition, both non-rotating as well as rotating BH solutions has been derived in these modified theories of gravity (Sen 1992;Moffat 2015;Shahzadi et al. 2019).In the following, we consider the several rotating BH solutions in many different modified theories of gravities.

Kerr-Sen BHs
Sen (Sen 1992) proposed a charged rotating BH solution to the equations of motion of the low-energy limit of the heterotic string theory, known as the Kerr-Sen solution.The metric coefficients of the corresponding line element can be written as with where Q K is the electric charge for Kerr-Sen BH, and the metric (35) reduces to the Kerr metric for the limiting value Q K → 0. The spacetime represented by the Kerr-Sen metric is not vacuum, and correspond to the KN case in the Einstein-Maxwell theory.There exists two horizons for a non-extremal Kerr-Sen BH, and can be determined by the condition ∆ K = 0 as where r − , and r + correspond to the inner and outer horizons of the BH respectively.The range of the parameter For extremal BH, both horizons coincide, and we have the condition

Einstein-Born-Infeld BHs
The gravitational field of a stationary and axisymmetric compact object with mass M , spin a and a non-linear electromagnetic source in the Einstein-Born-Infeld theory has been investigated by Julio Cirilo Lombardo (Julio Cirilo Lombardo 2004), and the metric coefficients of the spacetime read with where F denotes the Gauss hypergeometric function, Q BI shows the electric charge for Born-Infeld BH, η 2 (r) = Q 2 BI /r 4 β 2 , and β is the Born-Infeld parameter.For the limiting case β = 0, the metric (38) reduces to the Kerr metric while for β → ∞ (or Q(r) = Q = 0), we obtain KN metric.The Born-Infled parameter β can take any positive real value, and the charge Q BI is bounded by 0 < Q BI < 1.The metric Eq. ( 38) has curvature singularity at the points, where M = Q = 0, and Σ = 0.In an equatorial plane, it corresponds to a ring with radius a, and termed as a ring singularity.The properties of the rotating Einstein-Born-Infeld BH (38) are similar to that of the GR counterpart KN BH.Like the KN BH, Einstein-Born-Infeld BH is singular at ∆ BI = 0, and it admits thee static limit surface, two horizons like surfaces, and the event horizon.The horizons can be obtained by solving ∆ BI = 0, which are different from the KN BH, for details, see (Atamurotov et al. 2016).

Kalb-Ramond BHs
The Kalb-Ramond field is considered as a selfinteracting, second-rank antisymmetric tensor field in the heterotic string gravity.It can also be considered as a generalization of the electromagnetic potential with two indices, such that the gauge potential A α is replaced by the second-rank antisymmetric tensor field B αβ associated with the gauge-invariant rank-3 antisymmetric field strength et al. 2020).The stationay, axisymmetric, asymptotically solution of the modified field equations leads to the hairy BH solution The spontaneous Lorentz violating parameters s, and Γ are related to the vacuum expectation value of the Kalb-Ramond field and the non-minimal coupling parameter.
The free parameter s also known as Kalb-Ramond parameter, determines the potential deviation from the Kerr metric, and generalizes the KN metric.For the limiting case s = 0, Eq. ( 40) reverts to the Kerr BH, and to the KN BH for s = 1.The numerical solution of ∆ KR = 0 reveals that there exist only two positive real roots corresponding to the inner and outer horizons.Two distinct real positive roots infers the non-extremal BH, while no BH in the absence of real positive roots, i.e., no horizon exists.

Einstein-Gauss-Bonnet BHs
The uniqueness of the Einstein field equations is build on the Lovelock theorem which states that the GR with cosmological constant is the only theory of gravity in four-dimensional spacetime.However, the Einstein-Hilbert action is not unique in higher-dimensional spacetimes (d > 4, d represents the dimension of the spacetime), and the Einstein-Gauss-Bonnet theory of gravity is one of the interesting example in higher dimensions (Torii & Shinkai 2008).Recently, this theory has been proposed in lower dimensions by re-scaling the Gauss-Bonnet coupling parameter α → α/(d − 4), in the limit d → 4 (Glavan & Lin 2020).This theory has attracted much attention, and both non-rotating and rotating BHs have been proposed in this theory.The study of fourdimensional BH solutions presents a new approach to understand the Gauss-Bonnet gravity in low dimensions.The non-zero metric co-efficients of rotating Einstein-Gauss-Bonnet BH can the written in the form (Kumar & Ghosh 2020) with where α is the Einstein Gauss Bonnet coupling constant.

Konoplya-Zhidenko BHs
Konoplya and Zhidenko (Konoplya & Zhidenko 2016) proposed a rotating non-Kerr BH beyond GR and make an estimate for the possible deviations from the Kerr solution with the data of GW 150914, which can be regarded as a vacuum solution of an unknown alternative theory of gravity.The deformation changes the relation between the position of event horizon and BH mass, but preserves the asymptotic properties of Kerr spacetime.The non-zero metric co-efficients of the corresponding spacetime can be written in the form with The deformation parameter η of Konoplya-Zhidenko BH describes the deviations from the Kerr metric, and for vanishing η, one can obtain the usual Kerr spacetime.The presence of deformation parameter extends the allowed range of the spin parameter a and changes the geometry of BH in the strong field region.

Kerr-MOG BHs
The Kerr-MOG BHs are the stationary, axially symmetric and asymptotically flat solutions of field equations scalar-tensor-vector gravity which can be considered as another alternative to GR without DM in the present Universe, and can be described the line element with metric co-efficient (Moffat 2015;Kološ et al. 2020) where Here, G = G N (1 + α) is the enhanced gravitational constant, M is the mass of the BH, G N is Newton's gravitational constant.The dimensionless parameter α determines the gravitational field strength.For α = 0 and a = 0, the Kerr-MOG BH reduces to the Kerr and Schwarzschild-MOG BH, respectively.For details, see (Sharif & Shahzadi 2017, 2018).

Kaluza-Klein BHs
Kaluza-Klein BHs are the exact solutions of fivedimensional Kaluza-Klein (Einstein-Maxwell) theory.The simplest version of this theory is to study the GR in five dimensions.These theories are of great interest in string theory community because of their roles as lowenergy approximations to string theory.Larsen (Larsen 2000) proposed the most general BH solutions to the Einstein theory in five-dimensions and later dimensionally reducing the solution to four dimensions.The corresponding BH solution in an equatorial plane can be written using the follwing metric co-efficients (Larsen 2000) The unknown quantities takes the form where α ≡ a/M , x ≡ r/M , and b, γ are the dimensionless free parameters, for detail, see (Ghasemi-Nodehi et al. 2020).

BHs with Weyl corrections
The generalized Einstein-Maxwell theories have received a lot of attention recently because it contains higher derivative interactions and carries more information about the electromagnetic field.One of the simple generalized electromagnetic theories is the electrodynamics with Weyl corrections which involves a coupling between the Weyl tensor and Maxwell field.In this theory, the Lagrangian density of the electromagnetic field is modified as where F αβ is the electromagnetic tensor associated with the electromagnetic vector potential A α , and C αβγδ is the Weyl tensor.The coefficient α is a coupling constant with dimensions of length squared.The metric co-efficients describing the rotating BH with Weyl corrections can be written as (Chen & Jing 2014) with where q is the electric charge.For the limiting case q = α = 0, Eq. ( 51) reduces to the Kerr BH, and for α = 0, one can obtain the KN BH.

BHs in Rastall gravity
The Rastall theory of gravity is the modified GR, in which the usual energy-momentum conservation law (T αβ ;α = 0) is generalized to T αβ ;α = λR ,β , where T αβ is the energy-momentum tensor, λ is the Rastall parameter which represents the level of the energy-momentum conservation law in gravity theory.The metric co-efficients for rotating BHs in Rastall gravity reads (Xu et al. 2018a) with where N s is the surrounding fluid structure parameter, ψ is the Rastall coupling parameter, and ω s is the state parameter of surrounding fluid.For vanishing parameter N s , Eq. ( 53) goes over to the usual Kerr BH.For the limiting case ψ = 0, and −1 < ω s < −1/3, the metric shows the Kerr BHs surrounded by quintessence.

Charged Weyl BHs
In Weyl theory of gravity, the Einstein-Hilbert action is modified by a term proportional to the square of the Weyl tensor, and the action of the standard model of particle physics is modified to make it be conformally invariant.This theory has been appeared as one of the alternatives to compare the theoretical results with the cosmological as well as astrophysical observations, and also describe the cosmological parameters relevant to DE problem.The rotating BH solutions in this theory leads to the charged Weyl BH solutions (Fathi et al. 2021) with where m, q, and r is the mass, charge, and radius of source respectively, and ˜ is intended to recover the cosmological constituents of the spacetime, and has dimensions of m −2 .For λ > Q w , spacetime allows for two horizons; an event, and a cosmological horizon while there is a unique horizon for the extremal BH (λ = Q w ), and λ < Q w corresponds to the naked singularity.

Regular BHs in conformal massive gravity
The conformal massive gravity is an invariant theory under a conformal transformation of metric tensor as where Ω = Ω(x) is a nonsingular function of spacetime coordinates.It is noteworthy that the solutions of GR equations are a subset of the solutions of conformal gravity, and also it can describe the DM and DE scenarios.Recently, the regular BH solutions in this gravity has been proposed and the metric components takes the form (Jusufi et al. 2020) where Q c is the scalar charge, and λ 0 is the hair parameter.The metric is singular at the surface r = r sing , where Σ = 0.
2.5.12.Regular BHs in Einstein-Yang-Mills theory The dynamical interacting system of equations related to the non-abelian gauge theories defined on a curved spacetime is named as Einstein-Yang-Mills theory of gravity which describes the phenomenology of Yang-Mills fields interacting with the gravitational attraction, such as the electro-weak model or the strong nuclear force associated with quantum chromodynamics (Jusufi et al. 2021).The metric co-efficients of the regular, rotating, and magnetic charged BH solution with a Yang-Mills electromagnetic source in the non-minimal Einstein-Yang-Mills theory takes the form with where Q is the magnetic charge.This spacetime is free from the singularities and satisfies the energy conditions outside the outer horizon.For λ = 0, Eq. ( 60) reduces to the KN BH with a magnetic charge instead of an electric charge, and λ = Q = 0, we obtain the Kerr solution.

Hairy BHs
The modified Kerr BH solution so called rotating hairy BHs are surrounded by an axially symmetric "tensorvacuum" represented by a conserved energy-momentum tensor which could account for one or more fundamental fields (tensor, vector, or scalar fields representing any phenomenologically viable form of matter-energy, such as DM or DE).The energy-momentum tensor satisfies either the dominant energy condition or the strong energy condition in a region outside the event horizon.The metric co-efficients for hairy Kerr BHs can be written as (Contreras et al. 2021) with and the mass function m(r) takes the form The parameter α 1 determines the deviation from the Kerr BH, while the parameter α 0 measures the increase of entropy caused by the hair and must satisfy the condition α 0 ≤ 2M to ensure asymptotic flatness.For α 1 = 0, hairy BH solution reduces to the classical Kerr BH.

Rotating BHs modified by quintessence/matter field
Recently, with the help of Event Horizon Telescope's observations of BH shadows, it has been proposed that the existence of BHs in the universe is almost universally accepted (Event Horizon Telescope Collaboration et al. 2019).Inspired by this, many physicists have begun to study the interaction between DM (DM) and BHs (Kavanagh et al. 2020;Narzilloev et al. 2020;Xu et al. 2021a).Due to the existence of the supermassive BHs at the centers of galaxies, the strong gravitational potential of the BH concentrates a large amount of DM particles near the BH horizon (Gondolo & Silk 1999).The DM density increases by orders of magnitude due to the BH's gravitational field.Therefore, if DM particles can annihilate into gamma-ray radiation, the intensity of gamma-ray radiation near the BH will increase greatly, which provides a good opportunity for us to detect the DM annihilation signal.A series of DM models have been proposed in literature, some of them we consider here.

BHs in DM (dirty BHs)
The rotating BH solution surrounded by a spherical shell of DM can be expressed as (Pantig & Rodulfo 2020) with where the piecewise continuous mass function m(r), written in the form Here, ∆M < 0 and ∆M > 0 indicates the positive and negative energy density of matter, while r s , and ∆r s represent the inner radius, and thickness of the spherical shell of DM, respectively, for details see (Konoplya 2019).

BHs in perfect fluid DM
The non-zero metric co-efficients of rotating BH in perfect fluid DM can be written as (Hou et al. 2018) where The parameter k determines the perfect fluid DM intensity, and in the absence of perfect fluid DM, one can recover the Kerr metric.

BHs in cold DM halo
The strong gravity of a supermassive BH in the center of a galaxy could enhance the DM density significantly, producing a phenomenon known as "spike" (Gondolo & Silk 1999).But for the Navarro-Frenk-White density profile, a "cusp" problem occurs (de Blok 2010)a contradiction to the observations which show rather a flat density profile.However, for other DM models, i.e., scalar field dark DM, modified newtonian dynamics DM and warm DM, "cusp" is not produced in small scale.Motivated by these problems, Xu-et-al (Xu et al. 2018b) proposed the rotating BHs surrounded by DM halos, solution for the cold DM halos read with Here, ρ c denotes the density of the Universe at the moment when the halo collapsed and R c is the characteristic radius.For vanishing cold DM, Eq. ( 72) reduces to the Kerr BH with DM.

BHs in scalar field DM halo
The rotating BHs surrounded by scalar field DM halo can be described by (Xu et al. 2018b) with where, ρ s denotes the density of the Universe at the moment when the halo collapsed and R s is the characteristic radius.

Hayward BHs in perfect fluid DM
The solutions of the Einstein equations coupled to a nonlinear electromagnetic field in the presence of PFDM, representing the rotating and non-linear magnetic charged BHs surrounded by PFDM are given by (Ma et al. 2021) where where Q h is the magnetic charge of BH, and k denotes the intensity of PFDM.The value of k can be both positive and negative.In the absence of PFDM ( k = 0), one can obtain the rotating, non-linear magnetic charged BH, and for k = Q h = 0, Eq. ( 77) reduces to the Kerr BH.

BHs in DM spike
If the galactic center contains DM, then the existance of a supermassive BH in the galactic center would produce a cusp in the distribution of DM, known as DM spike (Gondolo & Silk 1999).In order to study the effects of DM spike on BHs, Nampalliwar et al. (Nampalliwar et al. 2021), proposed the rotating BH solutions immersed in DM spike, given by with where γ s = (9−2γ)/(4−γ), γ is the power-law index, ρ d denoted the density of DM spike, R sp shows the radius of DM spike, and R b is the inner edge of the DM spike.

Deformed BHs in DM spike
Recently, Xu et al. developed the deformed BH solutions immersed in DM spike, given by (Xu et al. 2021b) where Here, ρ R denotes the normalization of the DM density, α represents the power-law index, k 0 is the DM zero point parameter, R sp shows the radius of the DM spike, and R s is the Schwarzschild radius of BH.In the absence of DM spike (ρ R = 0), this BH reduces to the usual Kerr BH.

BHs in quintessence
The quintessence is dynamical and inhomogeneous scalar field having negative pressure, fully characterized by the equation ρ = wp, where ρ and p indicates the energy density and pressure, respectively.The non-zero metric co-efficients of rotating BH solution surrounded by quintessence can be expressed as (Xu & Wang 2017;Toshmatov et al. 2017b;Iftikhar & Shahzadi 2019) where c is the quintessential field parameter.There are three cases according to the value of state parameter ω, i.e., ω < −1, ω = 1, and −1 < ω < −1/3 , corresponds to the phantom energy, the cosmological constant and the quintessence, respectively (Sheoran et al. 2020).

ORBITAL PERIOD
The equations of motion for test particles in an alternative theories of gravity need not be geodesic.However, in the test-particle limit, equations of motion can be approximated as geodesics for a wide class of alternative theories, neglecting the spin of the small body (Vigeland et al. 2011).Here, we restrict our attention to theories where the modified equations of motion remain geodesic.
The orbital frequency describes the motion of test particles in the azimuthal direction, observed at radial infinity, defined by Ω φ = φ/ ṫ (where dot represents the derivative with respect to proper time τ ), and can be found with the help of the geodesic equation where Γ µ αβ are the Christoffel symbols which can be written in the form Due to the reflection and axi-symmetric properties of the spacetime, for the existence of equatorial circular orbits, we have Consequently, the equation that describes the particle motion in radial direction reduces to the following relation g tt,r ṫ2 + 2g tφ,r ṫ φ + g φφ,r φ2 = 0, ( Thus, the orbital frequency can be written the form where the upper and lower signs refer to the prograde, and retrograde orbits respectively.It is clear from Eq. ( 95), the orbital frequency Ω φ in independent of the metric co-efficients g rr and g θθ , while the partial derivatives of g tt , g φφ , and g tr with respect to the radial distance r are involved.The orbital frequency for those BHs having metric coefficients where some extra parameter is not multiplied with the radial distance r, is the same as for Kerr BH.We consider only those BHs which have different orbital frequency from that of Kerr BH.In the flare observations, the period of the orbital motion P is established, that is related to the orbital (Keplerian) frequency Ω φ by relation where period P is in minutes.This formula, hot-spot period as function of radius P (r), will be used for fitting observed flare period radius data in the following section.
Using the normalization condition p α p α = −µ 2 , evaluated in an equatorial plane, the effective potential V eff (r) can be expressed in the form where E = −p t , and L = p φ are interpreted as energy and axial angular momentum of a particle associated with Killing vector fields ξ µ (t) , and ξ µ (φ) , respectively.Effective potential is very important since it enables us to demonstrate the general properties of test particle dynamics, avoiding the necessity to solve the equations of motion.The circular equatorial orbits are governed by the condition (Kološ et al. 2015(Kološ et al. , 2017) ) V eff (r) = 0, dV eff (r) dr = 0. (98) The energy E and angular momentum L of circular orbits can be found by solving the Eq. ( 98), and the orbital frequency in terms of constants of motion can be written as Combining Eqs. ( 95), (99), along with the condition V eff (r) = 0, one can find the energy and angular momentum in terms of orbital frequency as where the upper and lower signs correspond to the prograde and retrograde orbits, respectively.The smallest stable equatorial circular orbits so called innermost stable circular orbits (ISCO) are governed by the Eq. ( 98) along with the condition d 2 V eff (r)/dr 2 = 0.The position of ISCO is one of the parameters that are very sensitive to the value of the BH spin.The location of ISCO for Schwarzschild BH (non-rotating) is situated at r = 6 from singularity.For rotating BHs, the ISCO of counter-rotating orbits move outwards the BH, while the position of ISCO for co-rotating orbits shift towards the BH.

FITTING TO THE FLARE DATA
Supermassive BH located at the dynamical center of our Galaxy has estimated mass of M ∼ 4 × 10 6 M from stellar dynamics around this gravitating center (Do et al. 2013), there is also more recent estimation to M = 4.30 × 10 6 M with a precision of about ±0.25% (GRAVITY Collaboration et al. 2021).The most prominent S2 star is orbiting Galactic center with pericenter at r p = 120 AU, which is bigger then 10 3 in geometrical units.The Galactic BH horizon diameter, which in geometrical units should be smaller then 4, is 0.3 AU only.Hence the mass M = 4.3 × 10 6 M measured be orbital dynamics of these S2 stars can be viewed as Newtonian far field limit on mass of relativistic compact object located at Galactic center.Hot-spot flares are observed from region much closer to the center, at r ∼ 0.7 AU (Venus orbit); with orbital period less then 60 minutes they are on fully relativistic trajectories, enabling us to explore true nature of our Galactic Center.
From formula (96), we notice that hot-spot period linearly depends on BH mass, see period-radius plots at Fig. 1.In the case of standard Kerr metric, all three BH masses M = (3, 4.3, 7) × 10 6 M fit the observed hot-spot flare data well, and the best fit can be seen for mass around M ∼ 4.5 × 10 6 M depending on BH spin.For large value of mass parameter, the circular orbits are situated above the center of error bars, while for small value of M , the circular orbits shift below the center.If BH mass M = 4.3 × 10 6 M will be used, we are able to fit observed data quite well, but still the hot-spot period appear to be little bit above center of error bars and one can conclude that the hot-spot is moving little bit faster then particle on geodesic in Kerr BH spacetime.Both co-rotating as well as counter-rotating circular orbits of neutral hot-spot in the background of Kerr BH with mass M = 4.3 × 10 6 M are depicted in Fig. 1.Orbiting periods around central BH are influenced by BH spin, but not strongly -all orbits fit the observed position and period of the three flares observed by GRAV-ITY.We compare the co-rotating (a = 1) as well as counter-rotating (a = −1) orbits with the Schwarzschild (a = 0) case and observe that the co-rotating orbits of hot-spot orbiting Kerr BH are situated above the orbits around Schwarzschild BH while the counter-rotating orbits lie below it.The dotted parts of the curves are plotted for regions below the ISCO position, while the solid part of curves show the behavior above the ISCO.If one would assume the hot-spot flare to be on stable circular orbit only, then counter-rotating orbits are not the option for higher BH spin (a < −0.4) since the ISCO is located too far away for counter-rotating orbit.
On the other hand, if one assumes the hot-spot could be also on unstable orbit below ISCO radius, then counterrotating orbits hot-spot orbits around mildly rotating BH a = −0.4 are giving better fits then co-rotating orbits.
Within this article, we consider some alternative BH spacetimes deviating from standard Kerr BH on observed radius-period data.We choose the BH mass M = 4.3 × 10 6 M , spin parameter a = 0.4, and for various different BH spacetimes, we fit the observed positions as well as periods of all three flares observed by GRAVITY in Fig. 2. Large variety (35 in total) of rotating BH spacetimes representing modifications of the standard Kerr geometry due to alternative gravity theories, or due to combinations of the standard GR vacuum Kerr spacetime with additional influences on the spacetime structure, we give the orbital periodradius relations for the BH with assumed fixed mass M = 4.3 × 10 6 M and dimensionless spin a = 0.4, with some representative values of the additional parameters representing the role of modifications of vacuum Kerr spacetime.
Total number of BH spacetimes used in this article is 35 so far.Only 31 plots of fits have been made, so there are 4 spacetimes, i.e, KN, braneworld BHs, dyonic  Hartle-Thorne q 1 =0.1 q 1 =0.5 q 1 =0.9 0.5 0.6 0.7 0.8 0.9 Kerr-Q charged BHs, and Kerr-MOG BHs, with the same orbital frequency as classical Kerr BH, hence they are not included in period-radius diagram (Fig. 2).The periodradius plots shown in Fig. 2 can give us some limitations on these alternative BH spacetime parameters where orbital period will be lowered down when such new spacetime parameters are introduced.We thus can immediately see the tendencies to increase or decrease the data fitting, being able to select the most favourable modified Kerr spacetimes.The Kaluza-Klein BHs, charged Weyl BHs, regular BHs and the BHs modified by DM field demonstrate clearly the best tendency to fit the data corresponding to observed three flares among all the considered BHs.Kerr-Sen Born-Infeld Weyl corrections Cold DM halo The detection of three bright "flares" in the neighbourhood of Galactic center supermassive BH that exhibited the orbital motion at a distance of about 6 − 11 gravitational radii from a 4.3 × 10 6 M BH has been declared by the near-infrared GRAVITY@ESO observations at 2.2 µm.The Ks-band observations also disclose that the flares are related with the orbiting luminous mass/hot-spots.We explore the dynamics of neutral hot-spot in the background of various stationary, axisymmetric, and asymptotically flat spacetimes, with the help of the three flares observed by GRAVITY on May 27, July 22, and July 28, 2018.We compared the co/counter-rotating equatorial hot-spot circular orbits around classical Kerr BH with different masses.
The equatorial circular orbits of hot-spot for most of the considering BHs are situated above the center of error bars of all three flares if we assume central object with mass M = 4.3 × 10 6 M and dimensionless spin a = 0.4, see Fig. 2. To obtain better fit of observed data, one can increase the central BH mass to M = 5×10 6 M or one can use some non-Kerr BH spacetime.In such alternative BH spacetimes new free parameters are introduced, which could be used to obtain better fit to observed GRAVITY data.The period-radius diagrams in Fig. 2 can give us some limitations on these alternative BH spacetime parameters where orbital period will be lowered down when such new spacetime parameters are introduced.Our results indicate that as favourable candidates can be considered the Kaluza-Klein BHs in Kaluza-Klein theory, charged Weyl BHs in Weyl gravity, regular BHs in conformal massive gravity and the BHs modified in DM field that are demonstrating the best tendency to fit the data corresponding to observed three flares among all the considered BHs.
All equatorial circular orbit period-radius relations of hot-spot orbiting the BHs, considering here both GR as well as modified theories of gravity, fit the observed position and period of all three flares observed by GRAVITY and it is hard to give conclusive answer on which alternative BH spacetime is correct -obviously more observations will be needed.However, it is instructive that between the best candidates belong all the considered cases of BHs influenced by DM in their vicinity.We also would like to note that the possible role of the electromagnetic interaction of slightly-charged hot spot (or slender torus) with magnetized Kerr BH has been considered in (Tursunov et al. 2020;Karas et al. 2021).Such non geodesic effects can also have strong influence on hot-spot dynamics and they are not considered in the article.

Figure 1 .
Figure1.Left: The relation between orbital period-radius, and position of three flares proposed by GRAVITY on Jully 22, May 27, and July 28, 2018 fitted with circular orbits of a neutral hot-spot in the background of Kerr BH with mass M = 4.3 × 10 6 M .The green curve is plotted for vanishing spin (a = 0) that corresponds to the Schwarzschild BH, while red and blue curves correspond to the co-rotating and counter-rotating orbits of hot-spot orbiting Kerr BH with spin a = ±1, fitting the observed positions and periods of the flares.The starting points of solid curves represent the ISCO positions, while dotted, and solid parts of the curves are plotted for below and above the ISCO position, respectively.Right: Orbits of neutral hot-spot moving around Kerr BH with three different BH masses M = (3, 4.3, 7) × 10 6 M , but having same spin a = 0.4, fitting the observed positions and periods of the flares.

Figure 2 .
Figure2.The relation between orbital period-radius, and position of three flares proposed by GRAVITY on July 22, May 27, and July 28, 2018 fitted with circular orbits of a neutral hot-spot for many different BHs both in GR as well as alternative theories of gravity.We choose the spin parameter a = 0.4, and mass M = 4.3 × 10 6 M for all BHs.Each plot represents three different orbits fitting the observed positions and periods of three flares.The black dotted curves indicate the Kerr limit.The red curves move towards the blue curves with the increase of the corresponding parameter.The text with each plot represents the designation of BH, and the values of parameters of corresponding BHs, used for the plots are also shown.We see that the circular equatorial orbits for most of the BHs are situated above the center of error bars of all three flares.