Accretion onto a static spherically symmetric regular MOG dark compact object

In astrophysics, the process of a massive body acquiring matter is referred to as accretion. The extraction of gravitational energy occurs as a result of the infall. Since it converts gravitational energy into radiation, accretion onto dark compact objects, e.g. black holes, neutron stars, and white dwarfs is an extremely significant process in the astrophysical context. Accretion process is a fruitful way to explore the features of modified gravity (MOG) theories by testing the behavior of their solutions associated with dark compact objects. In this paper, we study the motion of electrically neutral and charged particles moving in around a regular spherically symmetric MOG dark compact object to explore their related innermost stable circular orbit (ISCO) and energy flux. Then, we turn to investigate the accretion of perfect fluid onto the regular spherically symmetric MOG dark compact object. We obtain analytical expressions for four-velocity and proper energy density of the accreting fluid. We see that the MOG parameter increases the ISCO radius of either electrically neutral or charged test particles while it decreases the corresponding energy flux. Moreover, the energy density and the radial component of the four-velocity of the infalling fluid decrease by increasing the MOG parameter near the central source.

Phenomenologically, dark compact objects are an extensive family of astrophysical objects, which include black holes, neutron stars, white dwarfs, etc.In theoretical point of view, dark compact objects could be predicted in the context of extended gravity theories as well as in scenarios of the beyond the standard model of particle physics [1].Recently, observations of LIGO/Virgo proved the existence of binary black holes mergers through detection of gravitational waves [2][3][4][5], and additionally the Event Horizon Telescope (EHT) revealed the existence of supermassive black holes in center of galaxy M87 [6][7][8][9][10][11][12][13] and Milky Way [14][15][16][17][18][19].Therefore, it can be naturally anticipated that future advances in the field of gravitational wave astronomy and very long baseline interferometry will reveal new species of compact objects.On the other hand, it is fascinating to understand how and at what limits a dark compact object tends to be a black hole by increasing its compactness, which makes interesting the study of dark compact objects from a mathematical viewpoint.
General Theory of Relativity (GR) designed by Albert Einstein, besides a lot of achievements in explaining observations and predicting astonishing phenomena, is not yet the complete theory to describe gravitational interaction and corresponding events in the Universe.Reproduction of the rotation curves of nearby galaxies [20,21], mass profiles of galaxy clusters [22,23], intrinsic singularities at the center of black holes, etc are some examples of the failures of GR.Additionally, GR requires the cosmological constant term Λ to explain the positively accelerated expansion of the Universe at late-time [24,25].One interesting way to reform GR is to restructure the geometric part of GR through different approaches that can e.g.result in the so-called MOdified Gravity (MOG), which is a Scalar-Tensor-Vector (STVG) theory to describe gravitational interaction [26], proposed and developed by John W. Moffat.A massive vector field ϕ in addition to three scalar field as the mass of the vector field μ, the effective gravitational constant G, and the vector field coupling ξ are responsible for expressing the gravitational effects of spacetime in MOG setup.The MOG theory has several achievements in describing astrophysical observations, such as clarifying the rotation curves of many galaxies and the dynamics of galactic clusters without dark matter [27][28][29][30][31][32], in addition to compatibility with Planck 2018 data [33].Moreover, several black hole solutions including non-rotating and rotating ones [34] even with extra dimensions [35], cosmological solutions [36][37][38] and also, non-stationary solutions for inhomogeneity distributions of mass-energy in spacetime [39] are released within the framework of MOG theory in recent years.Also, many theoretical and observational efforts have been done to understand the MOG theory features and how it work in different situations [40][41][42][43][44][45][46][47][48][49][50][51][52][53].Interestingly, the solution describing the regular rotating and non-rotating MOG dark compact object has been recently explored in Ref. [54].The shadow behaviour of the regular rotating and non-rotating MOG dark compact object is investigated in Ref. [55].
Accretion is a process of particles being dragged onto a dark compact object.This process releases extra energy into surroundings, which is a source of some astronomical phenomena [56,57]; for instance the production of powerful jets, high-energy radiation, and quasars.A flattened structure made by rotating gaseous materials that slowly spiral into a massive central body is called an accretion disk.Accretion disks typically form around compact objects when interstellar matter exists.Accretion disks of compact objects are results of rotating gaseous materials in unstable bounded orbits [56,57].Under some conditions, the gas particles fall into gravitational potential of the compact objects, which causes gravitational energy in the form of heat.The inner portion of the accretion disk cools down as a result of the conversion of some heat into radiation [56,57].The electromagnetic spectrum of the emitted radiation can be analyzed when it reaches radio, optical, or X-ray telescopes.The motion of the gas particles, which may also be related to the structure and nature of the central mass, determines the properties of this radiation.As a result, studying accretion disk emission spectra can provide fruitful astrophysical data.Hence the accretion disks of compact objects drawn a lot of attention and have been studied in several cases in the literature .
The regular non-rotating (spherically symmetric) MOG dark compact object [54], which can be formed from the collapse of stellar object, can be tested in astrophysical phenomena, e.g.accretion process.It is the reason we found it interesting to study accretion disk onto the regular MOG dark compact object.In this regard, we also aim to study the motion of electrically neutral and charged test particles moving in this spacetime, and explore their corresponding energy flux.The rest of the paper is organized as follows.In Section II, we review the MOG field equations, and then, we introduce the regular MOG dark compact object spacetime and its features.Next, we study the motion of electrically neutral and charged test particles travelling in the regular MOG dark compact object spacetime in Section III.Then, we investigate the static spherically symmetric accretion in Section IV.Finally, we end with some conclusions in Section V.

II. ACTION AND FIELD EQUATIONS OF STVG THEORY
The total action in the theory of STVG is in the form of [26] where S GR is the Einstein-Hilbert action, S M is the action of all possible matter sources, S ϕ is the action of the (spin 1 graviton) vector field ϕ µ possessing the mass μ as one of the scalar fields in the theory, and S S is the action of three scalar fields, which can be expressed as follows in which g µν is the background metric tensor and g is the corresponding determinant, R is the Ricci scalar constructed by contracting R µν as the Ricci tensor, G is a scalar field in the setup, which known as the enhanced Newtonian parameter, ξ is third scalar field in the setup as the vector field coupling, V 1 (ϕ), V 2 (G), V 3 (μ), and V 4 (ξ) are the corresponding potentials of the vector field ϕ µ , and three scalar field G, μ, and ξ, respectively, and and also ∇ µ stands for the covariant derivative in the spacetime.
In the STVG theory, T µν = (M ) T µν + (ϕ) T µν + (S) T µν is the total stress-energy tensor, in which the stress-energy tensor of matter sources is (M ) T µν , the stress-energy tensor of the scalar fields is (S) T µν , and the stress-energy tensor of the vector field is for which V 1 (ϕ) = 0.One can find the full field equations of the STVG framework by variation of the action S concerning the inverse of the metric tensor, which yields [26] G in which the Einstein tensor is defied as G µν = R µν − 1 2 g µν R.
A. Regular MOG static spherically symmetric dark compact object The line element of the regular MOG static spherically symmetric dark compact object were found under the following assumptions [54] • The vector field is massless, i.e., μ = 0, since one can prove that for MOG compact objects, e.g., black holes possessing horizons, the mass of the vector field in the setup is zero.
• The enhanced Newtonian parameter G is defined as a constant depending on the free dimensionless parameter α so that G = G N (1 + α) where G N is the Newtonian constant.Furthermore, the gravitational source charge of the vector field is Q g = √ αG N M where M is the source mass.Here, we set G N = 1.
• The matter-free field equations of STVG setup is considered since the MOG dark compact object is a vacuum solution of the framework.
The above assumptions result in S M = S S = 0 and consequently, we have (M ) T µν = (S) T µν = 0. Thus, the field equations ( 6) now reduce to the following form Solving the last equation by following the procedure introduced in Ref. [54] leads to the line element of the regular MOG static spherically symmetric dark compact object as follows where dΩ 2 = dθ 2 + sin 2 θdφ 2 is the line element of the unit 2-sphere, and also we have defined which satisfies the weak energy condition [80,81].The MOG dark compact object possesses a critical value for α as α crit = 0.674 [54], so that for α ≤ α crit it has two horizons.It is worth mentioning that the (spin 1 graviton) vector field produces a repulsive gravitational force, which prevents the collapse of the MOG dark compact object to a MOG black hole with horizon.
Setting α = 0 in the line element (8) recovers the Schwarzschild black hole in GR.Moreover, the asymptotic behavior of the MOG compact object in the limit of r → ∞ is deduced as follows When, α ≤ α crit , the two horizons of the regular MOG static spherically symmetric dark compact object in the limit of r → ∞ can be found as When, α > α crit , there is a naked regular MOG static spherically symmetric dark compact object with no horizon.On the other hand, approaching the source, i.e., r → 0, the MOG dark compact object behaves to the form Therefore, the spacetime metric of the MOG dark compact object is regular so that f (0) = 1.Additionally, one can verify that the Kretschmann scalar R µνλσ R µνλσ in addition to the Ricci scalar R in the spacetime metric are regular at r = 0.
For the static spherically symmetric system, the gravitational redshift z at the asymptotic distance r to an observer is gathered as follows where the radius of the MOG dark compact object is R.For α < α crit , in the limit of r → ∞, the gravitational redshift of the compact object becomes infinite on the horizon r + and for α > α crit it has a finite value.Based on the observational data, however, one anticipates that the regular MOG dark compact object is adequately dark to be compatible with binary X-ray observations, so that α ∼ α crit [54].

III. MOTION OF TEST PARTICLE IN MOG DARK COMPACT OBJECT SPACETIME
The geodesic structure of the spacetime of the MOG compact object governs the trajectory of a test particle.In this section, we investigate the time-like geodesics around the regular MOG static spherically symmetric dark compact object through Lagrangian formalism [47,48,78,82,83].
Under temporal translation and rotation around the axes of symmetry, the line element (8) of the MOG dark compact object associated with the metric coefficient ( 9) is invariant since this spacetime is static and spherically symmetric.Therefore, the spacetime of the regular MOG dark compact object possesses two Killing vectors as follows These Killing vectors imply two conserved (constants) quantities for the motion of the test particle in the spacetime, which we aim to find them in the following.We plan to investigate the trajectory of both electrically neutral and charged test particles motion around the regular MOG dark compact object.

A. Motion of electrically neutral test particle
The Lagrangian of a test particle moving in the spacetime of the regular MOG dark compact object is expressed as where over-dot stands for derivative with respect to the affine parameter τ .The four-velocity of the test particle is defined as ẋµ ≡ u µ = (u t , u r , u θ , u φ ).We interested in the planar motion of the particle on the equatorial plane with θ = π 2 .Thus, utilizing the Euler-Lagrange equation one can find two conserved quantities of the particle motion corresponding with two Killing vectors as follows where E and L as two conserved quantities are the total energy and the total angular momentum per unit mass of the particle, respectively.Moreover, using the Euler-Lagrange equation, we can find dθ dτ = θ ≡ u θ = 0 in addition to Based on the normalization condition for the four-velocity of the test particle, i.e., u µ u µ = 1 and utilizing Eqs. ( 9), (17), and ( 18) one can find where V ef f is the effective potential of the test particle, which is defined as Effective potential analysis is significant in studying geodesic structure.The location of the circular orbits, for example, is determined by the local extremum of the effective potential.Figure 1 illustrates the behavior of the effective potential V ef f for the MOG dark compact object in comparison with the Schwarzschild case in GR.From Fig. 1 we see that increasing the value of the parameter α leads to increment of the effective potential.

FIG. 1:
The illustration of V ef f of the regular static spherically symmetric MOG dark compact object versus r for different values of α.The black solid line is for the case of Schwarzschild solution in GR.

Stable circular orbits around regular MOG dark compact object
The main characteristic of circular orbits is ṙ = r = 0 or equivalently u r = ur = 0. Hence, from Eqs. ( 17)-( 19) one can verify that for circular orbits, E 2 = V ef f and consequently, dV ef f dr = 0 must be satisfied.Solving these two equations simultaneously by using Eqs.( 9), (17), and (18) results in the following relations for the total (specific) energy E, total (specific) angular momentum L, and the angular velocity Ω φ ≡ dφ dt = u φ u t for the test particle in MOG dark compact object background where a prime stands for differentiation with respect to radial coordinate r and also we have defined According to Eqs. ( 22)-( 24) one can see that the condition 2f (r) − rf ′ (r) > 0 for existence of the circular orbits is required in order the total energy and the total angular momentum to be real.Figure 2 demonstrates the behavior of the E 2 versus r from which we can see that growing the parameter α leads to amplify the specific energy of the test particle in the spacetime of the regular MOG dark compact object while far from the source, the energy becomes almost constant.The corresponding curve of the Schwarzschild solution in GR is also shown in Fig. 2 which has always smaller values than the regular MOG dark compact object case.Figure 3 is the illustration of L 2 versus r associated with the regular MOG dark compact object in comparison with Schwarzschild solution in GR for different values of α, so that again increasing it results in growing the value of L 2 .All of these figures have smaller values of E 2 and L 2 than the corresponding ones in the case of the Schwarzschild solution in GR.In Fig. 4 we see the curves of Ω 2 φ versus r for the regular MOG dark compact object in comparison with the Schwarzschild case in which we see that increasing α leads to reduction of the value of Ω 2 φ so that the curve of the Schwarzschild case contains higher values of Ω 2 φ than corresponding ones in the regular MOG dark compact object.The location of the stable circular orbits correspond to the local minimum of the effective potential.Accordingly, an innermost (marginally) stable circular orbit (ISCO) needs the conditions to be satisfied.The existence of ISCO, r ISCO is purely a relativistic phenomenon.Instead of classical mechanics in which the effective potential possesses just one minimum; in GR however, the effective potential can generate either a local maximum and minimum or no extremum, relying on the choice of L in the effective potential.A stable outer and an unstable inner circular orbit for the test particle is related to this extremum.ISCO is where the stable and unstable circular orbits coincide for a specific value of L. Due to the complexity of the metric coefficient function (9) the explicit analytical form of ISCO associated with the regular MOG dark compact object is not available.Hence, solving equation set (27)

Radiant energy flux
In accretion process, the falling particles at infinity from rest will accrete onto the source mass.During the process, the gravitational energy of these falling particles will release and then convert into the electromagnetic radiation [56,84].One can express the radiation flux of the accretion disc around the central mass in the following form, which depends on the specific angular momentum, the specific energy, and the angular velocity of the falling test particle [56,84] where Ṁ is the accretion rate and g = det(g µν ) = −r 4 sin 2 θ is the determinant of the background metric tensor associated with the line element (8), so that on the equatorial plane, we have g = −r 4 .Inserting Eqs. ( 22)-( 24) into Eq.( 28) and also, using numerical data in Table I, one can find an approximate expression for radiation flux as follows where we have defined Figure 5 is the illustration of energy flux F(r) over r associated with the regular MOG dark compact object so that Fig. 5a is related to α = 0.09 and Fig. 5b is related to α = 0.2.From Fig. 5 we see that the energy flux in the setup is zero at r < r ISCO and then at r = r ISCO it grows from zero to infinity at r > r ISCO and after that, it again becomes zero at far from the source.Comparing Figs.5a and 5b demonstrates that increasing the value of α leads to decrease the energy in the setup.It also should be noted that the thermodynamical equilibrium is basic demand for the model describing the steadystate accretion disk.As a result, the radiation emitted from the accretion disk surface is equivalently as the black body spectrum [46,56,84].This means that the energy flux and the effective temperature of the accretion disk can be related with the well-known Stefan-Boltzman law F(r) = σ SB T 4 in which σ SB is Stefan-Boltzman constant.Therefore, from this law, one can find the effective temperature T of the accretion disk.Furthermore, at the distance d with the inclination angle γ to the central mass, the luminosity of accretion disk can be found as [46,56,84] where I(υ) is the thermal energy flux as function of frequency υ, while υ e = υ(1 + z) is the emitted frequency at the redshift z.Calculating the luminosity of accretion disk from the above equation is not possible, analytically due to the complexity of the relations.

B. Motion of electrically charged test particle
Assuming magnetic coupling process [85][86][87][88][89] in the vicinity of regular MOG dark compact object, the energy and angular momentum can be transferred from the dark compact object to the accretion disk.Therefore, on the horizon of the dark compact object, the strength of the magnetic field is expressed as [48] where the index (h) stands for horizon, c is the speed of light, and m p is the magnetization parameter so that m p = 1 means the equipartition state for the accretion and magnetic coupling process.Theoretical and experimental evidence demonstrate that the magnetic field can be exist in the surroundings of black holes and other compact objects [90][91][92].
Here, we suppose a weak magnetic field whose energy cannot influence the background geometry [93].Accordingly, this type of regular MOG dark compact object is called weakly magnetized.
Following the procedure introduced in Refs.[48,89,91], we aim to calculate the magnetic field in the surroundings of the regular MOG dark compact object.The Killing vectors introduced in Eq. ( 14) satisfy the following Killing vector equation [94] where 2 = ∂ µ ∂ µ is the d'Almmbert operator.In Lorentz gauge, the above equation is equivalent with the Maxwell equation for four-potential in which A µ is the four-potential and " ; " shows covariant derivative.The expression is related to a weak magnetic field, which is homogeneous at the spatial infinity with the strength B.Moreover, the magnetic field four-vector can be defined as follows where ϵ µνλσ is the Levi-Cività symbol, F λσ = A ν;µ − A µ;ν is the Maxwell tensor, andw ν is the four-velocity of a local observer at rest, which can be written as Utilizing Eqs. ( 34)-( 37) results in the magnetic field four-vector expression, which on the equatorial plane is as follows It is assumed that the magnetic field is directed upward along the z-axis at spatial infinity [95].Figure 6 demonstrates the illustration of B θ in vicinity of the regular MOG dark compact object versus r for different values of α.We see from Fig. 6 that far from the regular MOG dark compact object, the magnetic field is almost vanishing.Also, the effect of α on B θ is to reduce its strength.The Lagrangian of an electrically charged test particle with rest mass m and electric charge q travelling in the spacetime of the regular MOG dark compact object is expressed as Similar to the previous section, using Euler-Lagrange equation ( 16), on the equatorial plane we can find where ũµ is the four-velocity of the electrically charged test particle, Ẽ and L are the specific energy and specific angular momentum of the electrically charged test particle respectively.Again, on the equatorial plane, one can employ the normalization condition ũµ ũµ = 1 to gain where Ṽeff is the effective potential of the electrically charged test particle, which is Figure 7 is the illustration of Ṽeff versus r for the electrically charged test particle moving in the spacetime of the regular MOG dark compact object for different values of α.From Fig. 3 we see that increasing α firstly leads to increase the effective potential while in far from the source decreases it.

Stable circular orbits around regular MOG dark compact object for electrically charged particles
Similar to previous section, the conditions E 2 = V ef f and dV ef f dr = 0 must be satisfied for circular orbits.Therefore, one can solve these equations simultaneously by using Eqs.( 9), (40), and ( 41) to obtain the following relations Eqs. ( 44)-( 46) in the limit of q → 0 reduce to Eqs. ( 22)-( 24).The location of ISCO for the massive electrically charged test particle moving in regular MOG dark compact object satisfies the conditions (27).As we previously mentioned, the explicit analytical form of ISCO for the electrically charged test particle associated with the regular MOG dark compact object is not available due to complexity of the metric coefficient (9).Thus, one can numerically solve the equations set (27) by using, for instance, Wolfram Mathematica (v13.1) to obtain the numerical values of the ISCO for the electrically charged test particle moving in the spacetime of the MOG regular dark compact object.To do this, in Table II we again set M = 1 and for different values of α we collect the numerical values of rISCO , LISCO , and ẼISCO corresponding with the electrically charged test particle for the regular static spherically symmetric MOG dark compact object.Table II demonstrates that increasing the value of α causes to increase the ISCO radius of electrically charged test particle associated with the regular MOG dark compact object.Comparing Tables I and II shows us that the values of ISCO related to the weakly magnetised regular MOG dark compact object are smaller than the corresponding ones related to the regular MOG dark compact object.Therefore, the electric charge of the test particle and the magnetic field in the vicinity of the source affect the ISCO radius by reducing it.Similar to previous section, one can find the energy flux associated with the massive electrically charged particle moving in the regular MOG dark compact object spacetime by inserting Eqs. ( 44)-( 46) into Eq.(28).Also, the corresponding luminosity of the accretion disk can be found by Eq. (31).However, due to the lengthy and complexity of the related equations, it cannot be solved analytically.

IV. ACCRETION ONTO REGULAR MOG DARK COMPACT OBJECT
In this section, we aim to find the basic dynamical equations and parameters associated with the accretion onto the regular MOG dark compact object following the procedure performed in Refs.[76,78].To do this, we take into account the spherically symmetric accretion in the equatorial plane with θ = π 2 .Additionally, we assume that the accreting matter is inflowing perfect fluid onto the regular MOG dark compact object.

A. Dynamical equations
The perfect fluid stress-energy tensor is expressed as where p, ρ, and v µ are pressure, energy density, and four-velocity of the perfect fluid, respectively.On the equatorial plane, the only non-vanishing four-velocity components are v µ = (v t , v r , 0, 0).To be precise, the four-velocity of the perfect fluid v µ and the four-velocity of the test particle u µ in previous section are equivalent since the inflowing fluid, in fact, travels on the time-like geodesics creating the accretion disk around the MOG compact object.Therefore, the trajectory of inflowing fluid and the test particle in previous section are identical.In other words, the test particle in the previous section is assumed here as perfect fluid.On the other hand, according to the normalization condition for the four-velocity of the perfect fluid v µ v µ = 1 one can find where the condition v r < 0 must be satisfied since the accretion is an inward flow of matter while the assumption v t > 0 is taken into account because we interested in forward flow in time.
From the conservation of the stress-energy tensor, i.e., T µν ;ν = 0 in which (;) stands for covariant derivative, we can find the following relation where A 0 is a constant of integration.Additionally, we can project the stress-energy tensor conservation law onto the perfect fluid four-velocity to the form of which results in the following relation By integrating, the last equation yields where A 1 is a constant of integration.Since the condition u r < 0 holds, one can deduce where A 2 is an integration constant.Equation of mass flux in the setup is given by where on the equatorial plane results in the following relation where A 3 is an integration constant.

B. Dynamical parameters
Isothermal fluids with the equation of state p = ωρ where ω is the equation of state parameter are taken into account.During the motion of these fluids, the temperature remains constant.Consequently, Eqs. ( 52), (53), and (55) yields where A 4 is an integration constant.Inserting p = ωρ into the lest equation yields the v r as follows Figure 8 is the graph of v r versus r for the regular static spherically symmetric MOG dark compact object in comparison with Schwarzschild black hole in GR.From Fig. 8, we see that the fluid with the radial element of its four-velocity corresponding with each curve (associated with each value of α) begins to move towards the regular MOG dark compact object from rest at large r, as previously mentioned.Then, it approaches the regular static spherically symmetric MOG dark compact object to again reach the rest state.Furthermore, from Fig. 4, we see that decreasing the parameter α leads to increase the value of v r far from the source.Moreover, the curve of Schwarzschild case goes to infinity by approaching the regular static spherically symmetric MOG dark compact object.(ω + 1) Figure 9 demonstrates the illustration of ρ versus r for the regular static spherically symmetric MOG dark compact object in comparison with Schwarzschild black hole in GR.

C. Mass evolution
The central source mass of a black holes as well as a dark compact object is a dynamic quantity over time.Accretion process, for example, leads to grow their mass by accreting the surrounding matter onto them.The mass change measure or accretion rate of the regular MOG dark compact object can be obtained through Ṁ ≡ dM dt = − T r t dS in which the surface element of the object is dS = ( √ −g) dθdφ and also T r t = (p + ρ)v t v r .Consequently, the accretion rate Ṁ can be obtained as where the definitions A 0 ≡ −A 1 A 2 and A 2 ≡ (p ∞ + ρ ∞ ) f (r ∞ ) are assumed.Finally, we gain One can use Eq. ( 60) to obtain a relation between the initial mass M i and the mass in arbitrary time t as follows where the critical accretion time is defined as . At t = t cr , the mass of the regular MOG dark compact object approaches infinity in a finite time.

D. Critical Accretion
In accretion process, the inward flow of the fluid from rest at far from the source (regular MOG dark compact object) begins to move and continues to accelerate due to the gravitational field of the central source.During the inward flow motion of the fluid towards the source, it reaches sonic (critical) point, where the four-velocity of the fluid coincides the local speed of sound c s .From this critical point to the central source, the inward flow accelerated motion has supersonic velocities.A radial velocity gradient is needed to find the critical point.

FIG. 3 :FIG. 4 :
FIG. 2: The plot of E 2 of the regular static spherically symmetric MOG dark compact object versus r for different values of α.The black line is for the case of Schwarzschild solution in GR.

FIG. 5 :
FIG. 5:The behavior of F(r) versus r for the regular static spherically symmetric MOG dark compact object.

FIG. 6 :
FIG. 6: The illustration of B θ around the regular static spherically symmetric MOG dark compact object versus r for different values of α, where we have set M = 1.The black solid line is for the case of Schwarzschild solution in GR.

FIG. 7 :
FIG. 7: The illustration of Ṽeff of the massive electrically charged test particle moving in the spacetime of the regular static spherically symmetric MOG dark compact object versus r for different values of α.The black solid line is for the case of Schwarzschild solution in GR.

FIG. 8 :
FIG. 8: The illustration of v r of the regular static spherically symmetric MOG dark compact object versus r for different values of α.The black solid line is for the case of Schwarzschild solution in GR.

FIG. 9 :
FIG. 9: The illustration of ρ of the regular static spherically symmetric MOG dark compact object versus r for different values of α.The black solid line is for the case of Schwarzschild solution in GR.

TABLE I :
The numerical values of r ISCO , L ISCO , and E ISCO for a test particle moving in the regular static spherically symmetric MOG dark compact object spacetime associated with different values of α. α r ISCO L ISCO E ISCO

TABLE II :
The numerical values of rISCO , LISCO , and ẼISCO for an electrically charged test particle moving in the regular static spherically symmetric MOG dark compact object spacetime associated with different values of α.