Magnetized particle motion around magnetized Schwarzschild-MOG black hole

In this paper, we have presented the studies of the motion of magnetized particles and energetic processes around Schwarzschild black holes in modified gravity (MOG). The study of circular stable orbits shows that orbits of magnetized particles can not be stable for the values of magnetic coupling parameter β≥1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta \ge 1$$\end{document}. It was also shown that the range of stable circular orbits increases with the increase of both MOG and magnetic coupling parameters, while the effects of magnetic interaction stronger than the gravity. It was obtained that the increase of the MOG parameter causes the increase of center-of-mass energy collision of magnetized particles. Moreover, we have analyzed how to mimic the magnetic interaction with the spin of Kerr and Schwarzschild-MOG black holes. We have obtained that the magnetic coupling parameter can mimic the spin parameter a≤0.15\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a \le 0.15$$\end{document} (a≤0.28\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a \le 0.28$$\end{document}) giving the same radius of innermost contour(co)-rotating orbits at the values of the parameter β∈(-1,1)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta \in (-1,1)$$\end{document} and the MOG parameter in the range α∈(-0.17,0.28)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha \in (-0.17,0.28)$$\end{document} while the MOG parameter α∈(-0.7,0.9)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha \in (-0.7, 0.9)$$\end{document} mimics spin parameter of the black hole with the range |a|∈(0,1)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|a| \in (0,1)$$\end{document}.


Introduction
Detection of the dark energy and the dark matter in the Universe at the end of the 20th century changed our imaginations a e-mail: haydarovk378@gmail.com b e-mail: javlon@astrin.uz (corresponding author) c e-mail: ahmadjon@astrin.uz d e-mail: satimbay@yandex.ru e e-mail: dilfuza.fizic@mail.ru about Nature. The gravitational field described by the general relativity can not fully explain the nature of the dark energy and the dark matter. Among the different ways of explaining their nature the modified gravity (MOG) proposed in [1] became one of the most promising model. The theory proposed in [1] considers the scalar and massive vector field and usually called scalar-tensor-vector gravity (STVG). The motivation of introducing MOG is that being the classical field theory the general relativity breaks down at short length scales. In order to consider the quantum effects one needs the modification of the theory and MOG is one of the way to modify the general relativity. Introduced massive vector field with the source charge Q = √ αG M, where G is the gravitational constant, M is the mass of the central object, and α is the new coupling parameter, causes repulsive force, and becomes significant at the quantum level.
The non-rotating and rotating black hole solutions within MOG theory have been obtained in [2] and called Schwarzschil-MOG and Kerr-MOG black hles, respectively. The test particle motion around Schwarzschild-MOG black hole has been studied in [3], and it was shown the orbits become more stable due to the presence of a vector field in STVG theory. The several interpretations of MOG theory through solar system tests [1], galaxy rotation curve [4,5], through X-ray observations [6], black hole shadow [2,7], the study of thermodynamics [8], supernovae [9], gravitational lensing [10], quasinormal modes [11]. The thermodynamic properties of MOG theory have been studied in [12]. Epicyclic frequencies in Kerr-MOG black hole discussed in [13]. Quasi-periodic oscillations around Kerr-MOG black holes have been studied in [14]. Test particle dynamics near Kerr-MOG black hole have been considered in [15]. The gravitational instability in the Newtonian limit of MOG has been discussed in [16].
In the presence of electromagnetic field around black hole, one may study the motion of particles with non-zero spin and magnetic dipole momentum. In Ref. [42] it was shown that around Schwarzschild black hole immersed in an external asymptotically uniform magnetic field the magnetized particles can move along stable non-geodesic, spatially circular equatorial orbits with the radius smaller than Innermost stable circular orbits (ISCO). The study has been extended to the case of Kerr black hole in Ref. [43].
One of the authors of this paper has studied the magnetized particle motion around non-Schwarzschild black hole in the presence of a magnetic field [44]. Other our study is devoted to acceleration of magnetized particle around a rotating black hole in quintessence [27]. One of the autors of this paper has been involved to study the high energy collision of magnetized particles around a Hořava-Lifshitz black hole [28]. Magnetized particle acceleration around a Schwarzschild black hole in a magnetic field has been analyzed in [30]. One may find the analysis of magnetized particle motion around braneworld black hole in [45,46]. The magnetized particle motion in conformal gravity has been study in [47].
This work is devoted to studying the effect MOG to magnetized particle motion around a black hole and acceleration process of this type of particles. The paper is organized as follows: Sect. 2 is devoted to study the electromagnetic field and magnetized particle motion around Schwarzschild-MOG black hole. The magnetized particle acceleration near the Schwarzschil-MOG black hole has been studied in Sect. 3. We consider some astrophysical applications of our results in Sect. 4. We summarize our results in Sect. 5.
Throughout this work we use signature (−, +, +, +) for the space-time and geometrized unit system G N = c = 1 (however, for an astrophysical application we have written the speed of light explicitly in our expressions). The Latin indices run from 1 to 3 and the Greek ones from 0 to 3.

Magnetized particle motion around Schwarzschild-MOG black holes in magnetic field
The line element of Schwarzschild black holes in modified gravity can be written as [2] where and α is coupling parameter of MOG gravity. Consider the Schwarzschild-MOG black hole immersed in an asymptotically uniform magnetic field. The electromagnetic four-potential can be found using the Wald method [17] and expressed as: where B 0 is external uniform magnetic field. The non zero components of the electromagnetic tensor can be easily calculated using the definition F μν = A ν,μ − A μ,ν and have the following form The nonzero orthonormal components of magnetic field in the rest frame of the comoving observer have the following form Now we construct the equations of motion of magnetized particles around Schwarzschild black hole immersed in the external asymptotically uniform magnetic field in MOG theory. The Hamilton-Jacobi equation for magnetized particles can be expressed in the following form [42] where m is mass of the particle, S is the actin for magnetized particle in the curved spacetime background, the product of polarization and electromagnetic field tensors D μν F μν is responsible of the interaction between the external magnetic field and magnetized particles. The expression for the polarization tensor D μν corresponding to the magnetic moment of magnetized particles has the following form [42]: where μ ν is the four-vector of magnetic dipole moment and u ν is four-velocity of the particles in the rest frame of the fiducial comoving observer, being orthogonal to the magnetic moment. The electromagnetic field tensor F αβ can be expressed through electric E α and magnetic B α field components as Using the condition given in Eq. (8) one can easily calculate the interaction quantity D · F in the following form where μ = |μ| = μˆi μˆi = μ is the module of the dipole magnetic moment of the magnetized particles and L[λα] is the function of the coordinates and the other parameters of the spacetime around the black hole as well as magnetic field defining the tetrad λα attached to the comoving fiducial observer.
Here we study the orbital motion of magnetized particles around the Schwarzschild-MOG black hole in the weak magnetic interaction approximation, in other word we dropped out higher orders of D μν F μν 2 → 0. The conserved quantities in the equatorial plane are the angular momentum p φ = L and energy p t = −E of the particle. Now one can express the action of the magnetized particle in the form which can be used to seperate the variables in Hamilton-Jacobi equation. The equation of radial motion of the magnetized particles can be found in the following forṁ where newly introduced effective potential of radial motion has the form: where β = 2μB 0 /m is the magnetic coupling parameter and l = L/m is the specific angular momentum of the magnetized particle. The condition for the circular orbits can be expressed as follows: The first condition in (14) allows one to find the possible values of the magnetic coupling parameter β for circular orbits The second condition in (14) gives us We consider the particle at the equatorial plane with magnetic dipole moment perpendicular to the equatorial plane, thus the components of the external magnetic field measured by the observer in a frame of comoving with the particle take the following form Inserting Eq. (17) into (10) we can find the interaction part of the Eq. (7) where e = f − 2 r 2 − 1 2 [42] and is angular momentum and has following form Comparing Eq. (18) with Eq. (10) we get Finally, the magnetic coupling parameter β(r ; l, E, α) for stable circular orbits has the following form Equation (21) implies that a magnetized particle with magnetic coupling parameter β corresponds to circular stable orbit r with the energy E and angular momentum l. Figure 1 shows the radial dependence of the magnetic coupling parameter β for the different values of MOG parameter α for the fixed values of the specific energy and angular momentum. One can see from the top panel of Fig. 1 (when l 2 = 10) that the increase of the MOG parameter, α, causes the decrease of the maximum value of the magnetized parameter. One may also see that with the increase of the parameter, α, the loci where the parameter β = 0 shifts to the observer at infinity. However, when l 2 > 12 (bottom panel of Fig. 1), the maximum value of the magnetic coupling parameter β does not depend on the MOG parameter, α. Moreover, the local minimum of the β increases with the increase of α parameter and disappears for higher values of the parameter α. Now we analyze the values of the magnetic coupling parameter corresponding to the stable orbits of the magnetized particles. It can be found using following set of equations [42,44]: One can see (22) contains two equations with five parameters of the particle (β, r, l, E) and spacetime (α), so its solution can be parametrized in terms of any two of five independent variables. In order to solve the system of equations, it is better to use the magnetic coupling parameter and radius of the stable orbits r as free parameters. First, we will find the specific energy E and the angular momentum l of the magnetized particle as functions of radial coordinates and α parameter as: The expression (23) corresponds to the possible values of the specific energy of the magnetized particle at stable circular orbits. Figure 2 demonstrates radial dependence of the minimum values of the specific energy of magnetized particles at stable orbits for the different values of α parameter. One can see from Fig. 2 with the increase of α parameter the value of the specific energy decreases and the distance where the particle energy is zero increases.
Substituting (23) into (21) one may easily calculate the minimum value of the magnetic coupling parameter of magnetic particles for the given value of the specific energy in   Figure 3 illustrates the radial dependence of magnetic coupling parameter for different values of α parameter for the fixed value of the specific angular momentum l = 2 √ 2. One can see that with the increase the value of the parameter, α, the maximum value of the magnetic coupling parameter, β, and the distance where the magnetic coupling parameter, β, is zero increase.
Consider the upper limit for the angular momentum that the particle can be in stable circular orbits. The minimum value of the specific angular momentum corresponding to the minimum value of the magnetic coupling parameter can be found through the solution of the following condition ∂β min /∂r = 0 with respect to l 2 : Fig. 4 The radial dependence of minimal value of specific angular momentum of the magnetized particle for the different values of the α parameter  Figure 4 shows the radial dependence of the minimum value of the specific angular momentum for the different values of the MOG parameter, α. One can see from Fig. 4 that the value of the specific angular momentum corresponding to a stable circular orbit and the radius where the angular momentum is maximum increase with the increase the parameter α.
The extreme value of the magnetic coupling parameter β can be found by omitting Eq. (25) into Eq. (24) in the following form Figure 5 show the range of the magnetic coupling parameter corresponds to stable circular orbits. In the figure dashed lines correspond to the minimum value of the β parameter at l = 0 (freely falling magnetized particle) and solid ones correspond to the extreme value of the parameter β. Gray, light-blue and light-red colored areas correspond to the values of MOG parameter α = 0, α = 0.1 and α = 0.2, respectively. One can see from Fig. 5 that the distance where both the minimum and extreme values of the magnetic coupling parameter β are zero shifts to the observer at infinity. However, the width of the area corresponding to the fixed value of the magnetic coupling parameter, β, does not depend on the parameter α.
Thus, the extreme value of the parameter β corresponds to maximum value of the critical stable circular orbits r max and it can be found through the solution of the following equation with respect to r The minimum value for the circular stable orbits can be found solving the following equation with respect to r , The distance between maximum and minimum radius of circular stable orbits r = r max − r min give us the allowed area for the stable orbits for a magnetized particle. That means the circular stable orbits of a magnetized particle with given β are confined in the range r max (β; α) > r > r min (β; α). However, one can see from Eqs. (24) and (26) it is quite complicated to obtain the analytic solutions of Eqs. (27) and (28). We solve the equation numerically and present the results in a table form.
The area of stable circular orbits of magnetized particles for the different values of MOG parameter is presented in Table 1 corresponding to the angular moment from 0 to l min . One can see from the table the range r increases as the increase of both parameters α and β, however, the effect of the β parameter is stronger than effect of the parameter α.
One may express the dependence of minimum values of the specific angular momentum on magnetic coupling parameter β solving by the equation β = β min with respect to the specific angular momentum l

Particles collisions near the Schwarzschild-MOG black hole immersed in magnetic field
In this section, we investigate particle acceleration near the black hole of two colliding particles in the presence of the external magnetic field. Here we will study the effect of MOG parameter and external magnetic field to the center-of-mass energy of the colliding particles in-falling from infinity with specific (normalized to its mass) energies E 1 and E 2 . The center of mass energy of the two particles with the same mass m can be found using the expression [54]: where v α i (i = 1, 2) are the 4-velocities of the colliding particles.

Collision of two magnetized particles
The four-velocity of the magnetized particle at equatorial plane (θ = 0) has the following components: Consequently, the expression for center of mass energy of the two magnetized particles gets the following form: Figure 6 illustrates the radial dependence of center-ofmass energy of magnetized particles for the different values of the MOG parameter, α. The plots are taken for the values of magnetized parameter β 1 = β 2 = 1/2 and it was considered head-on collision of particles with the specific angular momentum l 1 = 2, l 2 = −2. One can see from Fig. 6 that the center-of-mass energy increases with the increase of the parameter α. The distance where the center-of-mass energy reaches its maximum increases with the increase of the parameter α.

Collision of two magnetized and charged particles
Now we will study the collisions of magnetized and charged particles. The four-velocity of a charged particle can be found using the Lagrangian for the charged particle in curved space in the presence of electromagnetic field: where e is electric charge of the particle. The energy and the angular momentum have the following form: and the four-velocity of the charged particle at equatorial plane has the following components: where ω B = eB/(2mc) is the cyclotron frequency responsible for the interaction between magnetic field and electric charge. One can easily find the expression for center-of-mass energy of magnetized and charged particles using Eqs. (37), (32), and (31): Figure 7 shows the radial dependence of center-of-mass energy of the colliding magnetized and charged particles near the Schwarzschild-MOG black hole. Plots are taken for the values of specific angular momentum of the particles l 1 = 2, l 2 = −2 and the magnetic coupling and cyclotron parameters β = ω = 0.5 for the different values of the parameter α. One can see from both panels of Fig. 7 that at large distances the center-of-mass energy disappears due to repulsive Lorentz forces which means that at large distance, the collision does not happen. However, in the case of collision of the magnetized and negative charged particle the distance where the center-of-mass energy disappears smaller Fig. 7 The radial dependence of the center-of-mass energy of collision of charged and magnetized particles with the same initial energy E 1 = E 2 = 1, around Schwarzschild-MOG black hole for the different values of α parameter. Top and bottom panels correspond to the cases of collisions of magnetized particle with negatively and positively charged particles, respectively than the case of the collision of magnetized particle with positively charged particles due to the orientation of magnetic field: in the first case the Lorentz force has repulsive nature and in the second one attractive nature. Moreover, in both cases the value of the center of mass energy increases with the increase of MOG parameter.

Collision of two magnetized and neutral particles
Now one can easily write the equations of motion for neutral particles around Schwarzschild-MOG black hole in the following form: Using the Eqs. (32), (39), and (31) one can easily find the expression of center-of-mass energy of collision of neutral Fig. 8 The radial dependence of the center-of-mass energy of collision of magnetized and neutral particles with the same initial energy E 1 = E 2 = 1, around Schwarzschild-MOG black hole for the different values of α parameter and magnetized particles in the following form: Figure 8 presents the radial dependence of center-of-mass energy of collision of neutral and magnetized particles with magnetic coupling parameter β = 0.5 for the different values of MOG parameter α. In this case one may also see that the increase of the parameter α causes the increase of the centerof-mass energy of the collision. However, the energy does not disappear due to the absence of Lorentz forces in the case of neutral particle. It implies that neutral particles can collide with other particles at any distance.

Collision of two charged particles
Here we will consider energetic process from the collisions of two charged particles. The expression for center of mass energy of the two charged particle can be obtained inserting Eq. (37) into (31) and we have The radial dependence of the center-of-mass energy of collisions of charged particles near Schwarzschild-MOG black hole for the different values of the parameter α is presented in Fig. 9. One can see from the figure in all cases the center-of-mass energy increases with the increase of the MOG parameter, α. In the case of negative and positive charged particles, the distance where the energy disappears increases with the increase of the value of the parameter α (bottom panel). Moreover, collisions of charged particles with the same sign shows that the center-of-mass energy decreases again at the larger distances due to repulsive Coulomb interaction between the colliding particles and the collision does not take place at higher energies (top panels).

Astrophysical applications
As an astrophysical applications of the studies of magnetized particles around Schwarzchild-MOG black holes, we consider analysis of ISCO radius for the magnetized particle around the black hole in an external asymptotically uniform magnetic field and rotating Kerr black holes. In other words, we look for the answer to the question: can magnetic interaction mimic the MOG and/or rotation parameters in the spectral fitting method implying necessity of additional methods in order to distinguish the Kerr black hole, the Schwarzschild-MOG black hole and presence of the external magnetic field. The expression for the radius of ISCO of the test particles around Kerr BH is given by the relations [58] where Our aim is to study and compare ISCO radius of a magnetized particle around: (1) Kerr black hole, (2) Schwarzschild-MOG black hole and Schwarzschild black hole immersed in an external asymptotically uniform magnetic field, corresponding the dimensionless values of rotation parameter a, α parameter and magnetic coupling parameter β in the range −1 to 1.
We show ISCO profiles of the magnetized particle around Schwarzschild black hole in the magnetic field (red-dashed line), Schwarzschild-MOG (dot-dashed blue line) and Kerr black holes (gray solid line) in Fig. 10. In this plot, the positive (negative) values of the magnetic coupling parameter correspond to the same (opposite) direction of the external magnetic field with dipole momentum of particles and we consider negative spin parameter as co-rotation of the particle. One can see from the profiles that the parameters can give the same ISCO radius corresponding to their different values. This implies that the parameters can mimic each others in the observations of ISCO radius of magnetized particles.
First, we consider the particle motion around Kerr and Schwarzschild-MOG black holes for the same ISCO radius.   Figure 11 shows the set of values of spin parameter of Kerr black hole and MOG parameter corresponding to the same value of the ISCO radius. One can see that positive (negative) values of the MOG parameter, α can mimic spin of Kerr black hole for the case of co-rotation (contour-rotation) of particles around the Kerr black hole. One may also see that the MOG Fig. 11 Relations between MOG and rotation parameters giving the same ISCO of the magnetized particles parameter mimics at α ∈ (−0.7, 0.9) spin parameter of the black hole at the range |a| ∈ (0, 1). Now we consider magnetized particle motion around Scwarzchild-MOG and Kerr black hole. Figure 12 shows the ISCO profiles of magnetized and neutral particles around Schwarzschild-MOG black hole in an external asymptotically uniform magnetic field and Kerr black hole. One can see from the figure that ISCO radius of particles around Kerr black hole r isco = 3M(7M) corresponds to the spin parameter a = 0.728533 (a = −0.395583) which fits the  The similar plot as Fig. 11, but for rotation and MOG parameters value of MOG parameter at α ∈ (−0.730477, −0.562357) (α ∈ (−0.081132, 0.424942)) for magnetized particles with the magnetic coupling parameters β ∈ (−0.5, 0.5). Now we will compare the ISCO radius of magnetized particles around the Kerr black hole and Schwarzschild black hole in the magnetic field and show how the magnetic field around Schwarzschild black hole can mimic the spin of Kerr black hole.
In Fig. 13 we show the relations between the magnetic coupling parameter, β, and spin parameter, a, of Kerr black giving the same ISCO radius. One can see from Fig. 13 that positive (negative) values of β parameter can mimic innermost co(contour)-rotating orbits of the particles around the Kerr black hole, giving the same radius of ISCO. One can see from Fig. 13 that the magnetic coupling parameter can mimic spin of Kerr parameter a ≤ 0.15 (a ≤ 0.28) giving the same radius of innermost contour(co)-rotating orbits at the values of the magnetic coupling parameter β ∈ (−1, 1).
Finally, we will investigate the magnetized particle motion around Swarzschild-MOG black hole and Schwarzschild black hole in the magnetic field.

Summary and discussions
In this work, we have studied the motion of magnetized particles around Schwarzschil-MOG black hole immersed in an external asymptotically magnetic field. Analysis of circular orbits show that the maximum value for magnetic coupling parameter corresponding to the specific energy E = √ 0.9 and angular momentum l = √ 10 decreases with the increase of MOG parameter. However, the maximum of the parameter β does not change with the change of the parameter α for the values l 2 > 12. The minimum value of the magnetic coupling parameter near the Schwarzschild-MOG black hole increases with increase of the α parameter and the minimum value disappears at α ≥ 0.1.
We have also studied the specific energy responsible for stable circular orbits and obtained that the energy decreases as the increase of parameter α. The studies of minimum and extreme values of β magnetic coupling parameter show that orbits of the magnetized particles can not be stable at β ≥ 1.
Numerical calculations of the range between maximum and minimum stable orbits expand with increasing both α and β parameters.
The investigations of collisions of magnetized particles have shown that the center-of-mass energy of the collisions increase with the increase of MOG-parameter, α.
Moreover, we have considered the similar studies of magnetized particles around Schwarzschild-MOG black hole immersed in an external asymptotically uniform magnetic field and Kerr black hole. Obtained that the magnetic coupling parameters β ∈ (−0.5, 0.5) and the MOG parameter at α ∈ (−0.730477, −0.562357) (α ∈ (−0.081132, 0.424942)) can mimic the values of the spin parameter at a = 0.728533 (a = −0.395583).
In our future works, in order to deeply understand the nature of modified gravity near the rotating black holes, we plan to extend the studies of magnetized particles motion around rotating Kerr-MOG black hole immersed in an external asymptotically uniform magnetic field. 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 .