Effects of Two Quantum Correction Parameters on Chaotic Dynamics of Particles Near Renormalized Group Improved Schwarzschild Black Holes

A renormalized group improved Schwarzschild black hole spacetime contains two quantum correction parameters. One parameter $\gamma$ represents the identification of cutoff of the distance scale, and another parameter $\Omega$ stems from nonperturbative renormalization group theory. The two parameters are constrained by the data from the shadow of M87* central black hole. The dynamics of electrically charged test particles around the black hole are integrable. However, when the black hole is immersed in an external asymptotically uniform magnetic field, the dynamics are not integrable and may allow for the occurrence of chaos. % English Editor: %Please verify that the intended meaning has been retained. Employing an explicit symplectic integrator, we survey the contributions of the two parameters to the chaotic dynamical behavior. It is found that a small change of the parameter $\gamma$ constrained by the shadow of M87* black hole has an almost negligible effect on the dynamical transition of particles from order to chaos. However, a small decrease in the parameter $\Omega$ leads to an enhancement in the strength of chaos from the global phase space structure. A theoretical interpretation is given to the different contributions. The term with the parameter $\Omega$ dominates the term with the parameter $\gamma$, even if the two parameters have same values. % English Editor: Please verify %that the intended meaning has been retained. In particular, the parameter $\Omega$ acts as a repulsive force, and its decrease means a weakening of the repulsive force or equivalently enhancing the attractive force from the black hole. On the other hand, there is a positive Lyapunov exponent that is universally given by the surface gravity of the black hole when $\Omega\geq 0$ is small and the external magnetic field vanishes.


I. INTRODUCTION
The images of the shadows of the supermassive black holes M87* at the center of the giant elliptical galaxy and Sgr A* in the Galactic Center [1,2] directly confirm the existence of black holes.They have not only been successful for testing Einstein's theory of general relativity in the strong field regime but have also been useful for testing other theories of gravity.Although the theory of general relativity has been proven to be a quite effective theory of gravity, it cannot completely describe various situations like the accelerating expansion of the universe and galaxy rotation curves [3].In this sense, general relativity is not the final theory of gravity and is necessarily extended to cure its shortcomings.Modified theories of gravity can be considered as generalizations of general relativity; in fact, there have been a number of modified gravity theories, such as f (R) gravity [4], scalar-tensor theory, tensor-vector-scalar theory, Hořava-Lifshitz gravity, Gauss-Bonnet theory, Einstein-AEther theory and quantum gravity (see, e.g., review papers [5,6] for more details on these modified gravity theories).
In addition to the mentioned-above limitations in ap- * Electronic address: 21200006@sues.edu.cn;wuxin˙1134@sina.complication, the existence of curvature singularities is one of the shortcomings of the theory of general relativity.Because these singularities are hidden behind event horizons, the pathologies at the singular region have no causal connection with the physics in the exterior region.The motion of physical objects outside the event horizons can be accurately described by general relativity black hole solutions.However, how physical objects fall inside is unknown due to the singular behavior of general relativity black hole solutions.Thus, the whole spacetime that contains both the exterior of a black hole and the interior of black hole is desired.Some regular black hole solutions can avoid the black hole singularity problem.Bardeen black holes [7] are a class of regular black hole models that are not exact solutions of the Einstein field equations.Of course, there are also singularity-free black hole solutions that are exact solutions of the Einstein field equations coupled to suitable nonlinear electrodynamics [8].The quantum gravity effects from loop quantum gravity and the renormalization group improvement (RGI) technique can also eliminate the emergence of black hole singularities.Bonanno and Reuter [9] gave an RGI Schwarzschild black hole obtained from the renormalization group improvement of the Einstein-Hilbert action with a running Newton constant.This black hole solution is based on the quantum effects in Schwarzschild black hole geometry.In practice, this solution is simply the Schwarzschild black hole, where Newton's constant is modified as a vari-able quantity depending on distance or time.There are two quantum correction parameters, including one parameter γ regarding the identification of cutoff of the distance scale and another parameter Ω stemming from nonperturbative renormalization group theory.When Ω is replaced with −Ω, the RGI Schwarzschild black holes become scale-dependent Planck stars [10].
Taking the lenses of the supermassive black holes Sgr A* and M87*, the authors of [11] investigated weak and strong deflection gravitational lensing of the RGI Schwarzschild black hole.They also provided a constraint in the quantum effects on the two quantum correction parameters γ and Ω.The authors of [12] discussed the circular motion of electrically charged particles around the RGI Schwarzschild black hole in the presence of an external asymptotically uniform magnetic field.They found that the radius of innermost stable circular orbit (ISCO) increases with an increase in the parameter γ but decreases with an increase in the parameter Ω or the magnetic field.RGI Schwarzschild black holes and the scale-dependent Planck stars can be distinguished in terms of the particles' periodic motions based on the precession of the S2 star from the GRAVITY observations, the shadow result of Sgr A* from EHT, and the gravitational waveforms [10].The authors of [13] showed that the variation in the parameters of the RGI Schwarzschild black hole causes a transition from periodic motions of time-like particles around the black hole to quasi-periodic motions.With the aid of the 3:2 resonance in X-ray binaries, quasi-periodic oscillation frequencies for the epicyclic motions of charged particles around the black hole immersed in an external asymptotically uniform magnetic field were studied.Such high frequency quasi-periodic oscillations in black hole systems arise from resonant phenomena between orbital and epicyclic motion of accreting matter [14][15][16].Periodic orbits and the epicyclic frequencies of quasi-circular motions around other quantum-corrected black holes were also considered in Refs.[17][18][19].
Although the external asymptotically uniform magnetic field is too small to affect the black hole geometry, it can exert an important influence on the trajectories of charged particles moving near the black hole.The trajectories are nonintegrable, and even allow for the chaotic behavior.Chaos is one kind of nonlinear phenomenon of a dynamical system, which exhibits sensitive dependence on initial conditions.The existence of chaos in some general relativity or modified gravity black holes has been reported in the literature, e.g., [20][21][22][23][24][25].Chaotic behavior may be useful to explain some astrophysical phenomena.Chaotic scattering of charged particles arises from the black hole gravitational field combined with the external uniform magnetic field.It provides a mechanism for charged particle acceleration along the magnetic field lines and is helpful for the formation of relativistic jets [26].In addition, self-similar fractal structures for characterizing the nature of chaotic photon motions can be observed in the shadow of nonin-tegrable black hole spacetimes [27,28].
As claimed above, the effect of varying the RGI black hole parameters on the ISCO dynamics of the charged particles was shown in Ref. [12], and the transition from periodic motions of time-like particles around the RGI black hole to quasi-periodic motions was studied in Ref. [13].Unlike the two works, the present paper mainly considers how varying the RGI black hole parameters causes a transition from regular dynamics to chaotic dynamics according to the constraint of the two quantum correction parameters γ and Ω in Ref. [11].For this purpose, we briefly introduce the dynamical equations of charged particles around the RGI black hole surrounded by an external asymptotically uniform magnetic field in Section 2.Then, we numerically investigate the regular and chaotic dynamics of charged particles in Section 3. Finally, the main results are summarized in Section 4.

II. DYNAMICAL MODEL OF CHARGED PARTICLES
An RGI Schwarzschild black hole metric [9] is briefly introduced.Then, an external asymptotically uniform magnetic field surrounding the black hole is described in terms of electromagnetic four-potential.Dynamical equations of charged particles around the black hole immersed in such an external magnetic field are obtained from a Hamiltonian system.

A. RGI Schwarzschild black hole metric
Suppose that the Newton's gravitational constant is G 0 .Bonanno and Reuter [9] modified the gravitational constant as a position-dependent quality where M is the mass of Schwarzschild black hole.γ and Ω are two new spatio-temporal quantum correction parameters without dimensions.The former parameter corresponds to the identification of cutoff of the distance scale, and the latter parameter corresponds to nonperturbative renormalization group theory.In this case, the function f (r) = 1 − 2G(r)M/c 2 is written as where the original Newton's gravitational constant G 0 and the speed of light c are taken as geometric units, G 0 = c = 1.On the other hand, quantum effects are included in the geometry of Schwarzschild black hole by means of the techniques of renormalization group.Considering the two points, Bonanno and Reuter obtained an RGI Schwarzschild black hole metric The two quantum correction parameters are constrained by the data from the shadow of the M87 central black hole.The data from the shadow of the M87 central black hole give the constraints of the two parameters in the ranges of 0.02 ≤ γ ≤ 0.22 and 0.165 ≤ Ω ≤ 9.804 [11].
Unlike Schwarzschild spacetime with γ = Ω = 0, RGI Schwarzschild spacetime with γ ̸ = 0 and Ω ̸ = 0 has finite curvature values of the Ricci scalar, square of the Ricci tensor and the Kretschmann scalar t r = 0, and γ ̸ = 0 [12].That is to say, the quantum gravity effects with the loop quantum gravity and the RGI approach completely rule out the formation of the black hole singularity.However, the existence of one or two horizons is still possible in the RGI Schwarzschild spacetime.In fact, there are two horizons for 0 < ε < 1, a single horizon for ε = 1, and no horizon for ε > 1, where Ω = εΩ + .Ω + is a root of the equation f (r) = 0 in the following form B. Electromagnetic four-potential An external asymptotically uniform magnetic field surrounding the RGI Schwarzschild black hole is assumed to be perpendicular to the equatorial plane θ = π/2.The magnetic field has an asymptotic value B. The spacetime of an RGI Schwarzschild black hole is asymptotically flat, and has two conserved Killing vectors along the time-like and space-like directions, ξ µ t = (1, 0, 0, 0) and ξ µ φ = (0, 0, 0, 1).According to the Wald method [29], an electromagnetic four-potential is a linear combination of the two Killing vectors in the form where c t and c φ are constant coefficients.A µ is a solution of the source-less Maxwell field equations.Taking c t = 0 and c φ = B/2, the authors of [12] gave an electromagnetic four-potential with only one nonzero covariant component This electromagnetic four-potential is used to describe the asymptotically uniform magnetic field.
Because the spacetime of an RGI Schwarschild black hole deals with the quantum gravity correction, it belongs to a kind of modified gravity theory.In general, such a modified gravity black hole is a nonvacuum solution or is non-Ricci-flat.In this case, the two constant coefficients c t and c φ from the Wald's formulas should be generalized as functions of the coordinates so that the obtained fourpotential is an exact solution to the source-less Maxwell equations, as was claimed by Azreg-Aïnou [30].See also the latest paper by Cao et al. [31].Does the covariant potential ( 6) not satisfy the source-less Maxwell equations?To answer this question, we rewrite Equation (2) as f Considering 1/r as a small value, we expand the function f (r) as This expression is an asymptotical expansion function with f 1 = M and f 2 = 0 mentioned in Ref. [30].Noticing the neutral black hole with charge Q = 0 in Equation ( 23) of [30], we can easily know that the covariant potential component ( 6) is what we want.In other words, the potential A φ satisfies the Maxwell equations.

C. Hamiltonian system
The magnetic field given by the four-potential ( 6) is so weak that it has a negligible effect on the spacetime geometry.In spite of this, it strongly affects the motion of a charged test particle with mass m and charge q.The motion is governed by the following Hamiltonian where p µ stands for a generalized covariant momentum In practice, it stems from one set of Hamilton canonical equations of the Hamiltonian system ( 9) ẋµ = ( ṫ, ṙ, θ, φ) represents a four-velocity, which is a derivative of coordinate x µ with respect to the proper time τ .Another set of Hamilton canonical equations of the Hamiltonian system (9) where g αν ,µ = ∂g αν /∂x µ .They show the presence of two constants of motion E and L correspond to the particle's specific energy E and the z-direction angular momentum, respectively.
For simplicity, dimensionless operations are given to the related variables and parameters in terms of scale transformations: , and H → mH.Thus, m and M in the above equations are two mass factors that can be eliminated.In this way, the Hamiltonian (9) becomes a system of two degrees of freedom where β = qB.For a time-like geodesic, the rest mass or the four-velocity ẋµ satisfying the relation g µν ẋµ ẋν = −1 corresponds to the conserved Hamiltonian This is a third motion constant of the system (9).If β ̸ = 0, there is no separation in the variables in the Hamilton-Jacobi equation of the Hamiltonian system (15).Therefore, the system ( 9) or ( 15) does not accept a fourth motion constant and is nonintegrable.If β = 0, the fourth motion constant is admitted and the RGI Schwarzschild spacetime is dynamically integrable.These facts show that the external asymptotically uniform magnetic field in Equation ( 6) does not change the spacetime geometry (3) but induces the nonintegrability of charged particle motions.

III. NUMERICAL INVESTIGATIONS
A time-transformed explicit symplectic integration algorithm is designed for the system (15).Then, it is used to provide some insight into contributions of the dynamical parameters to chaotic dynamics.

A. Setup of an explicit symplectic integrator
As aforementioned, the Hamiltonian system (15) is nonintegrable.Numerical integration methods are very convenient to solve such a nonintegrable system.Unlike traditional numerical methods such as Runge-Kutta integrators, symplectic integrators exhibit good numerical performance during a long-term integration of a Hamiltonian system.They conserve the symplectic geometric structure and give no secular drift to errors of motion constants during a long-term integration of a Hamiltonian system.In general, explicit symplectic integrators are superior to implicit ones at same orders in computational efficiency.Most of the Hamiltonians from curved spacetimes have no separation of variables or are not split into two explicitly integrable parts, and thus do not seem to allow for the application of explicit symplectic methods.However, more than two explicitly integrable splitting terms can exist in some curved spacetimes like the Schwarzschild black hole spacetime, and explicit symplectic methods are available without doubt [32].Unfortunately, the Hamiltonian for the Kerr black hole metric are not directly split into more than two explicitly integrable terms.Nevertheless, this problem can be solved with the help of an appropriate time transformation to the Hamiltonian and explicit symplectic methods can still work well [33,34].In what follows, we consider construction of an explicit symplectic integrator for the Hamiltonian system (15) along the idea of time transformation.
Let the proper time τ be regarded as a new coordinate q 0 = τ and its corresponding momentum be p 0 with p 0 = −H = 1/2 ̸ = p t .We have an extended phasespace Hamiltonian in the extended phase-space made of (p r , p θ , p 0 ; r, θ, q 0 ) Clearly, J is always identical to zero, i.e., J = 0, whenever the Hamiltonian H is conservative or not.Choosing a time transformation we have a new time transformation Hamiltonian The time-transformed Hamiltonian K has the following splitting pieces: where all sub-Hamiltonian systems are expressed as The six splitting parts have analytical solutions that are explicit functions of the new time w.Their solvers correspond to K 1 , K 2 , K 3 , K 4 , K 5 , and K 6 in turn.Setting h as a time step, we have a second-order explicit symplectic algorithm for the Hamiltonian (20) Because the Hamiltonian ( 20) is equivalent to the Hamiltonian (15), the symplectic method (28) designed for the Hamiltonian ( 20) is also suitable for the Hamiltonian (15).The solutions at the new time w should be changed into those at the proper time τ .
The integration time step is h = 1.The parameters are E = 0.995, L = 4.1, β = 4 × 10 −4 , γ = 0.22, and Ω = 0.166.The initial conditions are p r = 0 and θ = π/2.The initial separation are r = 15 for Orbit 1 and r = 92 for Orbit 2. The initial value of p θ (> 0) is given by Equations ( 15) and ( 16).The method S 2 yields no secular growth in the Hamiltonian errors ∆H = H + 1/2 of the two orbits in Figure 1a.This result is fit for a property of a symplectic integrator.When the integration lasts 10 7 steps, the proper time τ and the new time w are almost the same, as shown in Figure 1b.
In fact, the two tested orbits have different dynamical behaviors.Orbit 1 exhibits a closed torus on the Poincaré section θ = π/2 with p θ > 0 in Figure 1c.Such a torus is a characteristic of regular dynamics.Orbit 2 behaves as a thin area with many random distribution points and therefore is chaotic.Why is the Hamiltonian error for the regular orbit 1 larger than that for the chaotic orbit 2? It is because the regular orbit 1 has a smaller average period.Although no period exists in the chaotic orbit 2, an average period T is possible.When T increases, the Hamiltonian error with an order of (h/T ) 2 decreases.

B. Contributions of the parameters to chaotic dynamics
As mentioned above, the technique of Poincaré sections is a good description of phase space structures in the system (15) or (20).Taking the parameters E = 0.995, L = 4.5, γ = 0.02 and Ω = 8.166, we plot Poincaré sections for three values of the magnetic parameter β =0.0004, 0.0008, 0.0012 in Figure 2a-c.The plotted five orbits of charged particles around the RGI Schwarzschild black hole are regular for β =0.0004 in Figure 2a.Only one of the five orbits is chaotic for β =0.0008 in Figure 2b.However, the five orbits are chaotic for β =0.0012 in Figure 2c.These facts show that the occurrence of chaos becomes easier as the magnetic parameter increases from a statistical viewpoint of the global phase space structure.The result on a larger value of the magnetic parameter easily inducing chaos is also suitable for the Schwarzschild case with γ = Ω = 0 in Figure 2d-f.Comparing the two cases of an RGI Schwarzschild black hole and a Schwarzschild black hole, we find that the five orbits are regular for β =0.0004 in panels a and d, but chaotic for β =0.0012 in panels c and f.Chaos seems to be stronger in panel f than in panel c.In fact, an explicit difference between the two cases is that for β =0.0008,only one of the five orbits is chaotic in the RGI Schwarzschild case in panel b, whereas the five orbits are chaotic in the Schwarzschild case in panel e.Thus, for a given magnetic field parameter under some circumstances, chaos seems to be much easily induced in the Schwarzschild black hole spacetime than in the RGI Schwarzschild black hole spacetime.
Letting the parameters be L = 4.5, γ = 0.089, Ω = 3.166 and β = 8 × 10 −4 , we continue to consider the motion of charged particles around the RGI Schwarzschild black hole.When the energies are given three values of E = 0.999, E = 0.996, and E = 0.991 in Figure 3ac, the decrease in the energy typically decreases the extent of chaos and even rules out the occurrence of chaos.However, the extent of chaos is typically strengthened when the particle's angular momentum L decreases from L = 6.5 in Figure 3d to L = 5.5 in Figure 3e and to L = 4 in Figure 3f, where the parameters are Ω = 1.166,E = 0.995, γ = 0.089 and β = 8 × 10 −4 .Similar results exist as the renormalization group parameter Ω decreases from Ω = 9.166 in Figure 3g to Ω = 2.166 in Figure 3h and to Ω = 0.166 in Figure 3i, where the parameters are β = 8 × 10 −4 , E = 0.995, L = 4.5 and γ = 0.089.
The above-mentioned orbits can be checked by the maximal Lyapunov exponent, which was defined in Ref. [35] as where d(w) is a proper distance of two nearby orbits at the proper time w and d(0), one at the starting time.λ tending to zero shows regular dynamics for β =0.0004 in Figure 4a, but two positive stabilizing values of λ correspond to chaotic dynamics for β =0.0008, 0.0012.In particular, β =0.0012 with a larger Lyapunov exponent leads to stronger chaos than β =0.0008.A quicker method to distinguish between order and chaos is a fast Lyapunov indicator (FLI) of Froeschlé and Lega [36].This indicator was redefined in terms of the two nearby-orbit method [37] as The FLIs that exponentially increase with time log 10 w describe the onset of chaos for β =0.0008, 0.0012 in Figure 4b, while the FLI that algebraically grows with time indicates the presence of order for β =0.0004.The methods of λ and FLIs in Figure 4c-h display that a larger energy or renormalization group parameter has larger Lyapunov exponent and FLI for a given orbit, but a larger angular momentum corresponds to smaller Lyapunov exponent and FLI.Thus, the effects of the parameters E, L, and Ω described by the methods of λ and FLIs are consistent with those given by the Poincaré map method in Figures 2 and 3.
Taking the orbit with the initial radius r = 40 in Figure 2 as a tested orbit, we trace the FLI depending on the magnetic parameter β.The FLI for a given value β is obtained after integration time w = 5 × 10 6 .The FLIs ≥ 18.6 indicate disorder dynamics, while the FLIs < 18.6 show regular dynamics.In this way, the regular dynamics and chaotic dynamics in the RGI Schwarzschild black hole spacetime are clearly shown through the dependence of the FLI on the magnetic parameter β in Figure 5a.The dynamical features in the Schwarzschild black hole spacetime are also drawn in Figure 5b.The relations for the dependence of the FLI on the magnetic parameter are nearly the same in the two spacetimes.When the magnetic parameters are larger than some values, chaos can be allowed in the two spacetimes.The increase in the magnetic parameter makes the chaoticity become easy and strong under appropriate circumstances.This result does not mean that a larger magnetic parameter must induce stronger chaos, but is based on most of the values of the magnetic parameter larger than some values from the statistical viewpoint of the global phase space structure.However, there are some small differences between the two spacetimes.For β ∈ [0.000713, 0.001046], regular dynamics occurs in the RGI Schwarzschild case, but chaotic dynamics does in the Schwarzschild case.For β ∈ [0.001847, 0.001928], the RGI Schwarzschild case allows for chaos, while the Schwarzschild case does not.
Using the FLIs corresponding to two parameter spaces, we detect chaos from order in the two parameter spaces.For example, the initial radius is r = 40 and the FLIs are obtained from the two-parameter space (Ω, β), where 0.165 ≤ Ω ≤ 9.804 and 0.0002 ≤ β ≤ 0.0018.It can be seen clearly from Figure 6a that an increase in β enhances the strength of chaos, while an increase in Ω weakens the strength of chaos.Such an effect of varying the renormalization group parameter Ω on the dynamical transition of particles from order to chaos also exists in the twoparameter spaces (Ω, E) and (Ω, L) in Figure 6b,c.In addition, chaos occurs easily with the energy E increasing, but does not with the angular momentum L increasing.These results are also similar to those in Figures 2  and 3. Now, let us consider the impact of the distance scale γ ∈ [0.02, 0.22] on the dynamical transition.When the starting separation r = 40 and the parameters E = 0.995, L = 5.157, β = 8 × 10 −4 and Ω = 0.166 are given, three values γ = 0.08, 0.12, and 0.19 allow for the existence of chaos, as shown via the methods of Poincaré sections, Lyapunov exponents, and FLIs in Figure 7a-c.The dynamical transition is insensitive dependence on a variation of γ ∈ [0.02, 0.22] in the parameter spaces (γ, β), (γ, E), and (γ, L) of Figure 6d-f.The dynamical transi-tion insensitive dependence is based on the parameter γ constrained in Ref. [11].However, the ISCO radius increasing with the increase in the parameter γ in Ref. [12] is due to the parameter γ dissatisfying this constraint.
The main results can be concluded from Figures 2-7.Under some circumstances, chaos becomes easier as the energy E and the magnetic parameter β increase but does not when the angular momentum L and the renormalization group parameter Ω increase.A small variation in the distance scale γ seems to exert a negligible influence on the dynamical transition.In fact, these results are obtained from the statistical viewpoint of the global phase space structure.In order to explain them, we expand the third and fourth terms of Equation ( 15) outside the horizons as The second term presents that the black hole gives an attractive force contribution to the particle motion.The third term acts as a magnetic field attractive force and causes chaos to easily occur.The sixth term and the fifth term yield a repulsive force effect, which can somewhat suppress chaos.If γ ≪ 2Ω and γ ∈ [0.02, 0.22], the contribution of γ to the seventh term seems is negligible; if γ > 2Ω and γ ∈ [0.02, 0.22], the seventh term is still denominated by the sixth term.Why would a smaller value of the renormalization group parameter much easily induce chaos than a larger value?The question can be answered with the help of the universal relations between the chaotic test particle motions around the black holes and the surface gravities of such black hole horizons [38,39].Considering the motion of a particle near a horizon of a spherically symmetric black hole, the authors of [38] obtained a universal Lyapunov exponent provided by the surface gravity of the black hole.The authors of [39] estimated the radial distance of a massless and chargeless particle very near to the event horizon having exponential growing nature.For simplicity, we take γ = 0 in the present black hole geometry ( 2) and ( 3).The surface gravity of the black hole is expressed in [38] as Given Ω = 0, the spacetime (3) is the Schwarzschild black hole with the horizon r H = 2.In the case, the surface gravity near the horizon is When Ω = 1, the spacetime (3) also has a horizon r H = 1, and the surface gravity near the horizon is If 0 < Ω < 1, the black hole geometry (3) has two horizons r H = r ± H = 1± √ 1 − Ω, and the surface gravity near the horizon r = r + H is That is, this inequality is always satisfied for the case of 0 < Ω < 1.However, κ < 0 near the horizon r = r − H .In fact, the surface gravity of the black hole near the horizon was defined as a Lyapunov exponent λ H = κ [38].Thus, the Lyapunov exponent λ H near the horizon r = r H or r = r + H is positive for the case of 0 ≤ Ω < 1.The horizon would influence chaotic behavior in the motion of charged particles in the Hamiltonian system (15).Although this analysis is based on the case of γ = 0, it is also suitable for the case of γ ̸ = 0 because the small distance scale γ almost has a negligible influence on the dynamics.

IV. SUMMARY
An RGI Schwarzschild black hole metric is a nonrotating black hole solution in the quantum gravity as a kind of modified gravity theory.In this metric, the Newton's gravitational constant is modified as a positiondependent quality.There are also two spatio-temporal quantum correction parameters, including the identification of cutoff of the distance scale and a nonperturbative renormalization group theory.No black hole singularity exists due to the quantum gravity effects from the loop quantum gravity and the RGI approach.The RGI Schwarzschild black hole spacetime is similar to the Schwarzschild black hole spacetime that is integrable.In other words, the motion equations of test particles around the RGI Schwarzschild black hole have formally analytical solutions, and the particles' motion is regular.
When an external asymptotically uniform magnetic field surrounds the RGI Schwarzschild black hole, it is so weak that it does not change the spacetime geometry.However, it strongly affects the motion of electrically charged test particles near the RGI Schwarzschild black hole when the ration of charge to mass is large.Even the motion of charged particles is nonintegrable and may allow for the occurrence of chaos.A time-transformed explicit second-order symplectic integrator exhibits good long-term performance and effectively simulates the dynamical evolution of the nonintegrable system.
Using the methods of Poincaré sections, Lyapunov exponents and the FLIs, we investigate a varying parameter how to affect the dynamical transition from regular dynamics to chaotic dynamics.From a statistical viewpoint, chaos occurs easily under some circumstances as the specific energy and the magnetic parameter increase, but does not when the angular specific momentum and the renormalization group parameter Ω increase.A small variation in the distance scale has a negligible effect on the dynamical transition.An interpretation is given to the dependence of the dynamical transition on one varying parameter.In particular, a smaller value of the renormalization group parameter would much easily induce chaos than a larger value because of a positive Lyapunov exponent that is universally given by the surface gravity of the black hole when Ω ≥ 0 is small and the external magnetic field vanishes.In this case, the horizon would influence chaotic behavior in the motion of charged particles around the black hole surrounded by the external magnetic field.
Author Contributions: J.L. made contributions to the software, and writing-original draft.X.W. contributed to the supervision, conceptualization, writingreview and editing, and funding acquisition.All authors have read and agreed to the published version of the manuscript.
Funding: This research was supported by the National Natural Science Foundation of China (Grant No. 11973020).
Data Availability Statement: All of the data are shown as the figures and formula.No other associated data.

Fig. 1 :
Fig. 1: (a) The Hamiltonian errors ∆H = H + 1/2 for two orbits solved by the explicit symplectic integrator with step size h = 1.The parameters are E = 0.995, L = 4.1, β = 4 × 10 −4 , γ = 0.22, and Ω = 0.166.The initial conditions are pr = 0 and θ = π/2, and the initial separations are r = 15 for Orbit 1 colored Black and r = 92 for Orbit 2 colored Red.The errors show no secular drift.(b) The relation between proper time τ and the new time w.The two times are nearly identical.(c) Poincaré sections of the two orbits at the plane θ = π/2 with p θ > 0. Orbit 1 is regular, while Orbit 2 is chaotic.

Fig. 2 :
Fig. 2: Poincaré sections.The parameters are E = 0.995 and L = 4.5.In the RGI Schwarzschild spacetime with γ = 0.02 and Ω = 8.166, the magnetic field parameters are (a) β = 0.0004, (b) β = 0.0008, and (c) β = 0.0012.In the Schwarzschild spacetime with γ = Ω = 0, the magnetic field parameters are (d) β = 0.0004, (e) β = 0.0008, and (f) β = 0.0012.Clearly, chaos occurs easily in the two spacetimes as the magnetic field parameter increases.For β = 0.0004, the dynamical structure without the onset of chaos in the RGI Schwarzschild spacetime in panel (a) is similar to that in the Schwarzschild spacetime in panel (d).For β = 0.0008 in panels (b) and (e), more orbits can be chaotic in the Schwarzschild spacetime.For β = 0.0012, chaos becomes stronger in the Schwarzschild spacetime in panel (f) than in the RGI Schwarzschild spacetime in panel (c).

Fig. 4 :
Fig. 4: The maximal Lyapunov exponents λ and the Fast Lyapunov Indicators (FLIs).(a) λ and (b) FLIs for the tested orbit with the initial radius r = 33, and the parameters E = 0.995, L = 4.5, γ = 0.089 and Ω = 2.166, but three different values of the parameter β.(c) λ and (b) FLIs for the tested orbit having the initial radius r = 25, and the parameters same as those of Figure 3a-c.(e) λ and (f) FLIs for the tested orbit with the initial radius r = 80 and the parameters corresponding to those of Figure 3d-f.(g) λ and (h) FLIs for the tested orbit having the initial radius r = 40 and the parameters same as those of Figure 3g-i.

Fig. 5 :
Fig. 5: (a) Dependence of the FLI on the magnetic parameter β in Figure 2a-c.The initial radius is r = 40.The FLI for a given value β is obtained after integration time w = 5 × 10 6 .The FLIs ≥ 18.6 indicate disorder dynamics, while the FLIs < 18.6 show regular dynamics.(b) Dependence of the FLI on the magnetic parameter β in Figure 2d-f.

Fig. 6 :
Fig. 6: Distributions of two parameters corresponding to order and chaos via the FLIs for the tested orbit with the initial radius r = 40.(a) The two-parameter space is (Ω, β) and the other parameters are E = 0.995, L = 4.5 andγ = 0.089.(b) The two-parameter space is (Ω, E) and the other parameters are the same as those in Figure 3a.(c) The two-parameter space is (Ω, L) and the other parameters are the same as those in Figure 3d.