Charged particle dynamics in parabolic magnetosphere around Schwarzschild black hole

The study of charged particle dynamics in the combined gravitational and magnetic field can provide important theoretical insight into astrophysical processes around black holes. In this paper, we explore the charged particle dynamics in parabolic magnetic field configuration around Schwarzschild black hole, since the paraboloidal shapes of magnetic field lines around black holes are well motivated by the numerical simulations and supported by observations of relativistic jets. Analysing the stability of bounded orbits and using the effective potential approach, we show the possibility of existence of stable circular off-equatorial orbits around the symmetry axis. We also show the influence of radiation reaction force on the dynamics of charged particles, in particular on the chaoticity of the motion and Poincaré sections, oscillatory frequencies, and emitted electromagnetic spectrum. Applied to Keplerian accretion disks, we show that in parabolic magnetic field configuration, the thin accretion configurations can be either destroyed or transformed into a thick toroidal structure given the radiation reaction and electromagnetic-disk interactions included. Calculating the Fourier spectra for radiating charged particle trajectories, we find that the radiation reaction force does not affect the main frequency peaks, however, it lowers the higher harmonics making the spectrum more flat and diluted in high frequency range.


Introduction
Long-range forces provided by gravitational and electromagnetic (EM) interactions are of crucial importance for the proper understanding of high-energy processes around black holes (BHs).There are convincing observational evidences that magnetic fields (MFs) must be present in the vicinity of BHs [1,2].Orders of magnitude of MFs around BHs may vary from a few Gs up to 10 8 Gs and more, depending on the source generating the field.For stellar-mass BHs observed in X-ray binaries, the characteristic strength of MFs is of the order of 10 8 Gs, while for supermassive BHs (SMBHs), the characteristic strength is of the order of 10 1 − 10 4 Gs [1,3,4].Since the energy densities of MFs of such orders are not enough to make a sufficient contribution to the geometry of the background spacetime, in realistic astrophysical situations the spacetime metric around a BH can be fully described by the Kerr or Schwarzschild solution of the Einstein field equations.
The weakness of the MF around BH, in the sense that the MF does not affect the spacetime geometry, is compensated by the large value of the charge-to-mass ratio for elementary particles, whose motion will be essentially affected by MFs already of the order of few Gs.For a charged test particle with charge q and mass m moving in the vicinity of a BH with mass M immersed in MF of the strength B, one can introduce a dimensionless "magnetic parameter" B reflecting the ratio between the Lorentz force (LF) and gravitational force.For a relativistic electron orbiting a BH at the distance of the event horizon scale, one can estimate this ratio as follows that is very large due to large values of the specific charge q/m for elementary particles.Therefore, the influence of the EM force on the motion of charged particles can be dominant in realistic situations.Within the theory of general relativity, the gravitational influence of a central compact object (BH) can be well described by the Kerr or Schwarzschild spacetime [5], but unfortunately, we don't have such an elegant and well accepted solution for EM field influence -instead, there exist various models for BH magnetosphere depending on the complexity of included physical effects.As a zero approximation to the problem of BH magnetosphere one can use the vacuum solutions of Maxwell equations in the curved background as a starting point.In such approximation there is no source ρ = 0, or EM currents J α = 0 included, and the solution can be considered as an EM field of isolated BH in vacuum.However, electrovacuum solutions have limited astrophysical relevance, as they do not include the effects of material orbiting around BH in the form of a plasma.In the force-free approximation plasma is taken into account, but MF energy density is still dominating over matter energy density B 2 ρc 2 .In such models, the matter is used to generate EM currents inside the BH magnetosphere, but the matter is still locked in the MF.Plasma effect on the force-free magnetosphere of BH has been studied in the well-known work of Blandford & Znajek [6], where also the EM mechanism of rotational energy extraction from BH has been proposed.Several numerical techniques have been also employed, but the exact shape and intensity of the BH magnetosphere are not yet properly resolved, although a strong connection to the accretion processes is evident [7,8].One can also use the general relativistic magnetohydrodynamics (GRMHD) or general-relativistic particle-in-cell (GR-PIC) numerical simulations to test EM field configurations and try to construct a relevant BH magnetosphere model.A simple and still relevant BH magnetosphere model inspired by GRMHD and GRPIC simulations is the parabolic MF configuration [9][10][11][12].
MF effects on the astrophysical phenomena occurring in the vicinity of BHs cannot be neglected and the aim of this article is to find close connections between a simple realistic model for the BH magnetosphere with a radiative processes in the combined gravitational and EM fields.Simple charged particle model has limited validity in description of matter around BH -one should assume low particle density, i.e., that the mean time of two particle interaction is longer than orbital time around BH.One should also note that the typical Larmor radius of a particle in BH magnetosphere is much smaller than the Schwarzschild radius.Such low particle density assumption is correct only for the jet and for disks of low-luminosity active galactic nuclei (AGN), like Sgr A* or M87, while plasma inside thin disk around luminous AGN or XRBs may be quite collisional due to thin accretion disk's high density.
Our work can be seen as preliminary basic study of realistic BH magnetosphere, which should have according to the recent observation of synchrotron radiation and numerical modelling parabolic character [9].The article is structured as follows.We explore charged particle motion in parabolic BH magnetosphere model using effective potential approach in section 2, then we examine the influence of radiation reaction (RR) damping force on particle chaotic dynamics in section 3. We tested neutral Keplerian disk ionisation in section 4, estimate observable frequencies in section 5, and provide some astrophysical application in section 6.
Throughout the paper, we use the spacelike signature (−, +, +, +).Greek indices are taken to run from 0 to 3, while Latin indices run only over the spatial coordinates (r, θ, φ).However, for expressions having astrophysical relevance, we use the physical constants explicitly.

Charged particles in black hole magnetosphere
The line element describing the spacetime of nonrotating Schwarzschild BH with mass M is given as (2) where the function f (r) takes the form Hereafter, we put M = 1, i.e., we use dimensionless radial coordinate r (and time coordinate t).
In our magnetosphere model, we assume the axial symmetry and absence of electric field, hence the only non-zero component of EM four-vector potential A µ will be A φ , so that A µ = (0, 0, 0, A φ ).We will consider heuristic BH magnetosphere model given by the parabolic MF solution motivated by the GRMHD plasma simulations [12][13][14][15] where B ∈ (−∞, ∞) is MF intensity and w ∈ [0, 1.25] is declination of MF lines.Such parabolic field configuration should represent magnetic field inside BH jet region (close to z-axis) and more complex configuration is expected in accretion disk region (close to equatorial plane).One should note that ( 4) is a simplified version of the Blandford-Znajek (BZ) paraboloidal model [6] and it is not a solution of a vacuum Maxwell equation in Schwarzschild spacetime.The absolute value in A φ definition is making the solution to split in equatorial plane and preserving MF divergence, so MF lines go into the BH below the equatorial plane and outgoing above the equatorial plane, see Fig. 1.Such split magnetic configuration would require existence of conductive plasma current sheet in equatorial plane, thus a material orbiting around and falling into BH is supporting this model.In the case with w = 1 one gets a limit of BZ paraboloidal model and for w = 0 the BZ split monopole solution [6].The value, which is mostly used for BH magnetosphere in jet funnel is w = 3/4 [9].The orthonormal components of the MF B = (B r , B θ , B φ ), measured by zero-angular momentum ob-server (ZAMO) can be determined by [16] where η ijk is the Levi-Civita tensor in three dimensional space, and F ij = A j,i − A i,j is the antisymmetric tensor of the EM field.The non-vanishing orthonormal components of the parabolic MF take the form The magnitude of MF B as a function of r and θ takes the form The realistic accretion states from GRMHD simulations that get magnetically saturated do form such a split monopole state that lives transiently until accretion destroys it again (i.e. a thin current sheet at the base of a paraboloidal jet that looks like a split monopole nearby the BH) [17].The realistic or simulated BH magnetospheres [12,15] have a more complicated structure.Despite this, our simplified model (4) captures their most important features.One can distinguish three different regions around an accreting BH (see Fig. 1): the accretion disk/torus, where the matter density is high ρ ∼ ρ max and matter dominates over the turbulent MF; corona region, where the matter density is much lower ρ ≤ ρ max , but the matter still dominates over the still turbulent MF; jet funnel, where the matter component is very weak ρ ≤ 10 −6 ρ max and regular MF with parabolic shape dominates the region.The MF lines arising from the BH event horizon in equatorial plane are giving the boundary between the corona and jet funnel region.Corona-jet boundary condition is independent of magnetic intensity parameter B and reads In our simplified Schwarzschild BH magnetosphere model ( 4), the MF lines are not rotating and light cylinder will not be formed.Mostly due to the density differences, the mean free path for charged particles differs dramatically in all three (disk/corona/jet) regions.It is interesting to note that at least for M87 and Sgr A* parameters, also the disk is collisionless due to the small mass accretion rate and low plasma mass density.In this article, we will assume the description of collisionless charged test particle dynamics can be well applied at least in jet funnel and corona regions and hence we will focus on these two regions leaving the description of accretion disk to a more advanced GRMHD [9,15] or PIC [12,18] simulations.The charged test particle motion is described by the covariant Lorentz equation where u µ is the four-velocity of the particle with the mass m and charge q, normalized by the condition u µ u µ = −1, τ is the proper time of the particle, and is the antisymmetric tensor of the EM field.Using Hamiltonian formalism for the charged particle motion, we can write where the kinematical four-momentum p µ = mu µ is related to the generalized (canonical) four-momentum π µ by the relation that satisfies the Hamilton equations in the form The affine parameter ζ of the particle is related to its proper time τ by the relation ζ = τ /m.Due to the symmetries of the Schwarzschild spacetime (2) and the parabolic MF (4), one can easily find the conserved quantities that are the energy and the axial angular momentum of the particle and can be expressed as The dynamical equations for the charged particle motion in the Cartesian coordinates can be found by the transformation equations Introducing for convenience the specific parameters, energy E, axial angular momentum L, and magnetic parameter B, by the relations one can rewrite the Hamiltonian (11) in the form Fig. 3 For attracting LF (B < 0), we plot the radial positions of stable circular orbits for charged particles in the off-equatorial plane.Different MF configurations w = 0, 0.75, 1 have been plotted.The grey disk shows the BH horizon, the numbers on each solid curve indicate the values of the MF parameter B, while the points of these curves show the position of stable off-equatorial circular orbits.The dotted curve represents the offISCO for charged particles in the off-equatorial plane, while the dashed curve shows the corona/jet boundary.
where V eff (r, θ; L, B) denotes the effective potential given by the relation The terms in the parentheses correspond to the central force potential given by the specific angular momentum L, and EM potential energy given by the magnetic parameter B. The Hamiltonian (18) can be divided into dynamical part H D (first two terms containing dynamical momenta p r , p θ ) and potential part H P (last term).
The effective potential (19) shows clear symmetry (L, B) ↔ (−L, −B) that allows to distinguish the following two situations − minus configuration, here L > 0, B < 0 (equivalent to L < 0, B > 0) -MF and angular momentum parameters have opposite signs and the LF is attracting the charged particle to the z-axis, towards the BH.+ plus configuration, here L > 0, B > 0 (equivalent to L < 0, B < 0) -MF and angular momentum parameters have the same signs and the LF is repulsive, acting outward the BH.
In this article, without loss of generality, we use the positive angular momentum of a particle L > 0, while the magnetic parameter B can be both positive or negative.Fig. 4 For repulsive LF (B > 0), we plot the radial positions of the inISCOs for charged particles in equatorial plane, with different MF configurations w = 0, 0.1, 0.75, 1.The dotted vertical line shows the position of the BH horizon, while the dotdashed vertical line indicates the position of the inISCO for a neutral particle.

Off-and in-equatorial circular orbits
The charged particle motion is confined by the energetic boundary given by Let us properly investigate the features of the effective potential (19) represented in Fig. 2 that enables us to demonstrate the general properties of the charged particle dynamics, avoiding the necessity to solve the equations of motion.Since the V eff (r, θ) function is symmetrical above θ ∈ (0, π/2) and below θ ∈ (π/2, π) equato- rial plane, we can focus on the region θ ∈ (0, π/2) only, and we can substitute | cos θ| with cos θ function here.The effective potential as a function of coordinates x-z is plotted in Fig. 2.
The effective potential V eff (r, θ) function is zero at BH horizon and limits are for all values of magnetic parameters B and w ∈ (0, 1).Only for value w = 1, we have in the case of first limit The stationary points of the effective potential V eff (r, θ) function, where maxima or minima can exist, are given by the equations These equations for the region above the equatorial plane θ ∈ (0, π/2) are where we define the supplementary function The term containing | cos θ| in V eff (r, θ) function (19) will create the discontinuity in the first derivative of V eff (r, θ) with respect to θ in an equatorial plane θ = π/2, and such points must be also added when we would like to examine V eff (r, θ) function extrema.
There exist two classes of effective potential extrema solutions: off-equatorial plane extrema for attractive LF (B < 0), and in-equatorial plane extrema for repulsive LF (B > 0).
For MF parameter B < 0, we have minima located at off-equatorial plane while radial coordinate r can be calculated using angular momentum condition L = L off (r), where The new function Y (r) is given by Off-equatorial plane minima are limited by the minimal distance from the BH -so-called innermost stable circular orbits (offISCOs) given by the equation where the particle angular momentum L is given by Eq. ( 27).The radial positions of stable circular orbits as well as offISCO situated in the off-equatorial plane for different values of magnetic parameter B < 0 and declination parameter w are depicted in Fig. 3.Both offISCO and off-equatorial stable circular orbits get closer to the z-axis with the increase of the magnetic parameter B. The offISCOs are depending also on MF configuration, been close to the BH for pure magnetic monopole w = 0 than in the case of parabolic MF w = 3/4.Existence of the off-equatorial circular orbits in the background of magnetized BHs with external asymptotically uniform or dipole MFs [19][20][21] inspired the idea of investigation of the possible existence of charged off-equatorial toroidal structures orbiting around magnetized BH.One should note, that BHs can not have their own long-lived dipolar fields due to the no-hair theorem [22,23].Radiation belts filled with energetic charged particles exists around many different astrophysical objects (mostly with dipole magnetosphere), so also in the case of BH one can search for observational consequences.As it can be seen in Fig. 3, the off-equatorial plane minima are located closer to the z-axis (and within the jet funnel) with increasing MF magnitude B forming off-equatorial torodial structures in low density jet funnel filled with charged particle orbiting around BH on the off-equatorial circular orbits.For repulsive LF (B > 0), we have minima located in equatorial plane θ = π/2, while radial coordinate can be calculated using angular momentum condition L = L in (r), where Effective potential minima in equatorial plane are limited by in-equatorial innermost stable circular orbit (inISCO), given by the equation The radial position of inISCO for different magnitudes of the MF B and various values of parameter w, situated in an equatorial plane are depicted in Fig. 4. The inISCO radius is constant r inISCO = 6 only in the case of split magnetic monopole magnetosphere w = 0 [24].
For other values of declination parameter w of MF, the radii of inISCOs lie in the interval 2 ≤ r inISCO ≤6.For large MF magnitude, the inISCO radii to the BH horizon r inISCO = 2, however, when the MF is negligible, we obtain the limiting case r inISCO = 6.Similar behaviour can be observed for uniform MF [25,26].
The graphical behaviour of the particle energy E(r) as function of radial coordinate r for circular orbits located in both off-and in-equatorial planes is shown in Fig. 5. Before the offISCO (or inISCO) position, the energy of stable orbits decreases with the increase of radial coordinate, but after the offISCO (or inISCO) position, the energy increases monotonically as the radial distance increases.As it can be seen from Fig. 5, charged particle on the off-equatorial plane circular orbit has lower binding energy than particle on circular orbit in the equatorial plane.One can conclude that off-equatorial circular orbits are less stable against perturbation than the circular orbits located in equatorial plane.
Limiting cases, when the MF is negligible (|B| 1) or when the MF is dominating over gravity (|B| 1) are more astrophysically prevalent -examples of particle orbits for these two cases are plotted in Fig 6 .In order to see clearly the interplay between magnetic LF and gravitation BH attraction leading to strongly nonlinear and chaotic dynamics, one should use magnetic parameter |B| ∼ 1.The effect of attractive (B < 0) and repulsive (B > 0) LF is demonstrated on particle trajectories around off-and in-equatorial minima in Fig 7 .Effective potential minima can naturally serve as traps for particles and one could expect higher par-Fig.7 Charged particle trajectories (blue) around BH with parabolic magnetosphere on orbits around the off-and in-equatorial plane minima.Gray disk represents the BH, gray curves the MF lines, the dashed curve the energy boundary and the black dot the initial position of the particle.
ticle concentrations there, but the situation in/close to the equatorial plane must be described by a more complex approach than our single charged particle dynamics, for example by GRMHD of PIC simulations since Keplerian accretion disk formed by conductive plasma is expected to orbit BH in equatorial plane.

Radiating charged particle trajectories
In many astrophysically relevant scenarios, one cannot neglect the effects of RR due to the synchrotron radiation of charges in the vicinity of BHs, which are believed to be immersed in an external MF.Equation of motion for a charged particle undergoing RR force in curved spacetime is non-trivial [27].However, for elementary and subatomic particles, the equations can be simplified to the form [28] Du , where semicolon denotes the covariant coordinate derivative.The Eq. (32), which is a covariant form of the Landau-Lifshitz equation [29], is a habitual second order differential equation which satisfies the principle of inertia and does not contain runaway solutions.The exact equation of motion expressed for our particular parabolic magnetosphere are rather complex and hence we do not present exact formulas here, however we solve them numerically.The effect of RR is characterized by the radiation parameter Detailed discussion and derivation of (32) can be found in [28].One should note, that in our approach the charged particle is radiating due to acceleration given by the LF which has many orders of magnitude stronger effect than the so-called geodesic synchrotron radiation [30,31]; in our approach such radiation on geodesic given by the tail integral [28] can be neglected.Detailed analysis of the particular case of the motion of charged particles around Schwarzschild and Kerr BHs immersed into an external asymptotically uniform MF was presented in [28,32].
Trajectories of radiating charged particles around Schwarzschild BH immersed in parabolic MF are depicted in Fig. 8.For both attractive (B < 0) and repulsive (B > 0) LFs, the charged particle will lose energy and angular momenta due to the RR force, which will force it to spiral down to the BH.For radiating particle trajectories plotted in Fig. 8, we use the same initial conditions as for non radiating ones presented in Figs. 6 and 7.For repulsive (B > 0) LF, we can observe oscillation damping due to reaction leading to particle's temporary stabilization on circular orbit ("orbital parking").This "parking" effect can be observed only in Fig. 8 Influence of RR on charged particle dynamics in parabolic (w = 0.75) magnetosphere of BH.The initial conditions for particle dynamics are the same as in Figs. 6 and 7.The radiation parameter k has been chosen unrealistically large in order to clearly demonstrate the effect of RR force.Due to RR the particle spirals down to the BH.The damping of oscillations can be observed in the first row.Time scales related to RR depend on MF parameter B and specific charge, they are relevant for the strong MF, while negligible for weak one.
repulsive B > 0 case and similar effect for repulsive LF has been already observed in the case of BH in uniform MF [28].For parabolic MF, the situation is different than for uniform one, here the circular parking orbit is not stable and particle will be slowly in-spiraling down to the BH.The RR force acts as a damping force -the oscillatory character of the motion of charged particles changes because of the loss of energy and angular momentum and radiating charged particles will always be captured by the BH.The off-equatorial plane minima are unstable for radiating particles and charged particles can stay oscillating around off-equatorial plane minima only for a limited amount of time, which could still be quite long since RR is relevant only for relatively strong MFs.
It is well known that the neutral test particle motion in Schwarzschild background is regular and fully inte-grable.On the other hand, the charged particle dynamics around magnetized BHs is generally chaotic [21,33].The motion of charged particles is regular (integrable) only for magnetic monopole field configuration w = 1 due to monopole's spherical symmetry [24].Strongly non-linear and chaotic trajectories can be seen in the case of parabolic MF w = 3/4, see lower row in Fig. 6.Still, in parabolic MF not all orbits are chaotic: the trajectories corresponding to minima of effective potential or with small epicyclic oscillations around circular orbits are regular [25,34].Such almost circular orbits are very important from the astrophysical point of view as they govern the thin (Keplerian) accretion disks and the toroidal fluid configurations.
Charged particle chaotic behaviour can be analysed with the help of the Poincaré section (PS), presented in Fig. 9. From the distribution of the plotted points Fig. 9 Charged particle chaotic dynamics (blue curve) in parabolic MF (gray parabolic lines) around BH (gray circle) can be well demonstrated by the PS plots (last column).In the first row, we see complete regular motion represented by smooth the Kolmogorov-Arnold-Moser (KAM) torus section.In the second row, we observe that this KAM torus is already starting to break and chaos is emerging with this orbit.In the third row some special (unstable) periodic orbit is presented, with only single point in PS.In the fourth row one can observe completely chaotic orbit, with points fully covering its phase space, limited only by particle's energy (dashed black curve).We use MF parameters B = 2, w = 0.75 for all orbits.
in the PS, one can distinguish whether the motion is chaotic or not.For regular motion, the plotted points form a closed curve in the 2-dimensional r − p r plane, while in the case of the chaotic motion, PS points are distributed apparently randomly in the allowed region.When the energy of the particle is small, tori of small harmonic oscillation are not broken, see first row of Fig. 9.Such tori get broken when the energy of the particle is increased, see second row of Fig. 9, and the chaotic motion appears.The sea of chaos spreads in the PS with the increase of the particle energy, see fourth row of Fig. 9.The process of chaos increase with energy is sometimes interrupted by regular orbits, see example of unstable periodic orbit in third row of Fig. 9.
The comparison between PS for orbits of particles with and without the RR force is plotted in Fig. 10.The RR force acts as a damping force leading the radiating particle towards the BH attractor, hence trajectories at PS points are moving radially downwards in time, as one can see the visible accumulation of points even-Fig.10 PSs for chaotic particle orbit without (blue first row) and with (red second row) the influence of RR.The same trajectories as in Fig. 9 have been used.The RR is a non-linear damping force acting on a charged particle which can increase the chaoticity of a trajectory, but can also lead to the particle collapse to the BH.
tually leading towards r = 2, p r = −1.The radiation losses are highly non-linear during inspiral towards BH and particle may experience multiple resonances [35], which is demonstrated as structures in radiating particle PS.Unstable periodic orbit, third column in Fig. 10, will get perturbed by RR and become chaotic.Already chaotic, fourth column in Fig. 10, will remain chaotic.

Ionized Keplerian disks
Single charged particle dynamics is a toy model for plasma description, but even with this simple tool one can test combined gravitational and EM effects of BH magnetosphere keeping control on numerical scheme, which could be clouded in more complex GRMHD or PIC numerical simulations.To demonstrate the influence of the MF on processes in BH vicinity, we consider an ionization of Keplerian accretion disk consisting of electrically neutral test particles following circular geodesics and orbiting a Schwarzschild BH immersed in a parabolic MF, assuming the symmetry of Keplerian disk plane almost coinciding with the equatorial plane of the Schwarzschild geometry.
Realistic ionization model of the originally neutral matter has been applied e.g. for the description of the magnetic Penrose process [36,37], where the original neutral particle 1 splits into two oppositely charged particles -2nd and 3rd, with charges q 2 and q 3 , re-spectively.We thus assume conservation of the electric charge 0 = q 2 + q 3 , and the canonical momentum where π α(i) , i = {1, 2, 3} denotes the canonical momentum of 1st, 2nd and 3rd particles, respectively.Due to the charge conservation, the momentum conservation Eq. ( 34) takes the form In many realistic scenarios one of the newly created charged particles is much more massive than the other one, m 2 /m 3 1, such as in the case of a neutral atom ionization.Example of ionization mechanism for electrically neutral particle orbiting Schwarzschild BH immersed in parabolic MF is shown in Fig. 11.The magnitude of BH magnetosphere is set to B, but due to difference in masses of electrons and protons (ions), the LF (B parameter) and RR forces (k parameter) are very different.Therefore more massive charged particle (proton or ion) takes almost all the initial momentum of the original neutral particle, and the dynamical influence of the lighter charged particle (electron) can be neglected, so that In another realistic scenario, we can consider the Keplerian disk created by plasma modelled as a quasineutral soup of charged particles: electrons and ions, which are orbiting around the central object in circular orbits.If the disk is dense enough, the mean free path of charged particles is much shorter in comparison to the length of the orbit around the central object.This means that the charged particles orbit the BH together as a collective neutral body along circular geodesics as the influence of the MF is irrelevant to the motion, but at the edges of the disk its density decreases substantially and the charged particle motion starts to be significantly influenced by the MF as the mean free path becomes comparable to the orbital length.
In both scenarios neutral particles are ionized, while their mechanical momenta are conserved.This simple ionization model (36) has been studied previously in the field of rotating Kerr BH with uniform MF [38][39][40].Chaoticity of the motion of the ionized matter orbiting Schwarzschild BHs has been studied in detail in [33].
Parabolic MF around non-rotating BH is given by A φ component of the EM four-potential, hence particle's energy (p t ) is conserved during the ionization process (35).Particles from neutral Keplerian disk are on bounded orbits (−p t = E ≤ 1), hence they can not escape to infinity after ionization since this would require (E > 0), see (19).Charged particles are able to escape along open parabolic MF lines to infinity if they obtain initial kick in vertical direction -such a scenario has been described in the uniform MF case in [41][42][43], where the parameters of the initial kick in vertical direction were calculated.For Kerr BH, the situation is different since here BH rotation generates "electric" component A t even in the case of a static MF configuration.The presence of the A t component leads to the energy redis-tribution in the formula (34), which can lead the particle to escape from the bounded Keplerian circular orbit of BH with relatively high energy due to the so-called magnetic Penrose process [37,38,44,45].
Specific axial angular momentum L and specific energy E for neutral test particles on the inclined circular orbits from Keplerian accretion disk are given by [46] Due to a simple ionization condition (35), using the definitions of the specific energy and specific axial angular momentum (39), one can write for the specific axial angular momentum L and the specific energy E of the ionized test particle Only the specific angular momentum L is changed during the ionization, while the particle specific energy E remains constant.Charged particles from the ionized Keplerian disk can either be captured by BH, oscillate around circular orbit, or move chaotically along parabolic MF lines.The influence of the MF parameter B on the fate of the Keplerian disk orbiting a Schwarzschild BH in a position corresponding to the symmetry plane of the parabolic MF is demonstrated in Fig. 12.When the disk remains neutral or the MF is missing (B = 0 case), all the orbits remain circular -original inclined razor thin disk survives.For small values of the magnetic parameter (B = ±0.01cases), the charged particles will be unsettled from original circular orbits and ionized disk will start to oscillate in both radial and vertical directions; the ionized disk becomes slightly thick due to the epicyclic motion of the charged matter.When the MF is slightly increased (B = ±0.1 cases), complete destruction of the inner part of the Keplerian disk occurs in the attractive B = −0.1 case, as all the orbits end at the BH.In the repulsive B = 0.1 case, the regular epicyclic motion has a tendency to chaoticity and the ionized disk becomes thick.Further increasing the MF parameter (B = ±1 cases) the charged particle motion enters a fully chaotic regime causing a substantial thickening of the disk.For large values of the MF parameter (B = ±10 cases), the LF becomes the leading force of the particle motion.LF in both attractive and repulsive cases leads to similar results: the charged particles go up and down along MF lines realizing epicyclic motion around BH in clockwise (B > 0) or counter-clockwise (B < 0) directions.
The influence of the RR on the ionized Keplerian accretion disks orbiting Schwarzschild BH with parabolic MF can be observed for realistic values of the radiation parameter k only in a very long integration time; in Fig. 13 we have used k = 0.02 to make RR effect more visible.The RR force reduces the chaotic behaviour of the particle motion.Neutral disk ionization and subsequent radiation losses do not lead to effective filling of the off-equatorial plane minima (B < 0 case) with Fig. 13 PSs for particles from ionized Keplerian disk, see Fig. 12.Two upper plots correspond to trajectories without RR: the left one (B = 0.1) consists from regular trajectories forming closed loops, while the right one (B = 2) from chaotic trajectories, where the points are distributed apparently randomly; the grey dotted vertical lines show the inISCO positions.In the two lower plots, the RR is included and we can see that some trajectories are attracted to the BH.All radiating particle trajectories may eventually end up in BH, however for small B and k, or if the orbit is far away from the BH, the accretion takes much longer time.
charged particles -off-equatorial minima have higher energy than equatorial one, and off-equatorial minima are not attractors for radiating charged particles.Different processes similar to Earth's MF trapping charged particles from the solar wind in the Van Allen radiative belts can also be considered to form off-equatorial plane structures.It is worth mentioning that that the physics of pair cascades around BHs depend heavily on RR and Compton scattering [47][48][49].

Particle oscillations and observed frequencies
Ionization of particles from neutral accretion disk in the equatorial plane may lead to the disk oscillations for |B| < 1.If a test particle is slightly displaced from the equilibrium position, corresponding to a stable circular orbit, the particle starts to oscillate around the equilibrium position realizing thus epicyclic motion governed by the linear harmonic oscillations.In order to calculate the epicyclic frequencies of small perturbations to such orbits, we employ a linear stability analysis to compute the radial and vertical frequencies using a general method [25].Let us introduce a vector of canonical variables y K = (x i , p i ), so that we rewrite the equations of motion derived from the Hamiltonian (18) in the form ẏK = f K (y N ), (41) where dot denotes the derivative with respect to the proper time τ .We are interested in the solutions to the Eq. ( 41) near equilibrium position y K 0 corresponding to a circular orbit, where f K (y N 0 ) = 0, except N = 3. Linearizing y K = y K 0 + ζ K , where ζ K /y N 0 << 1 denotes a small deviation vector, one gets a set of linear differential equations with constant coefficients We decompose the solutions to this system into eigenvalues and eigenvectors of the Jacobian matrix ∂f K /∂y N .The eigenvalues λ specify the rate at which trajectories with small differences in initial conditions separate.Thus, we evaluate Eq. ( 42) using the Hamiltonian (18) with the circular orbit values y 0 and find a pair of eigenvalues that we interpret as radial and vertical angular frequencies given by where the unknown function χ takes the form ∂ 2 H p ∂θ 2 .Fig. 14 Off-equatorial radial profiles of the frequencies of small harmonic oscillations νr, ν θ , ν φ , and ν L of charged particle orbiting BH of the mass of M = 10M in external parabolic MF for different values of the parameters B and w. (45) There exists the third fundamental angular frequency of the epicyclic particle motion, namely the Keplerian (axial) frequency ω φ , given by The radial and vertical angular frequencies (43-44) measured by a local observer for the case of in-and offequatorial circular orbits have too complicated character when expressed explicitly.For w = 2 and in equatorial plane, Eqs.(43) and (46) reduce to the uniform MF configuration [25], while B = 0 leads to the neutral particle case [46].In contrast to the neutral case, there exists angular frequency so-called Larmor frequency for charged particles, which is associated with the MF only, and given by a relation where |B| can be found by Eq. ( 8).Due to the form of the Hamiltonian H P (18) in parabolic MF, the first and second derivatives of H P with respect to θ does not exist at the equatorial plane, hence the construction of ω θ is not possible at the equatorial plane.However, one can still perturb circular orbit in vertical directions around the equatorial plane, so that the perturbation will lead to the non-harmonic oscillations.Equations for the harmonic and non-harmonic oscillators are given, respectively, as follows where C θ is a constant.Large vertical perturbation from equatorial plane lead to chaotic dynamics, see the second row in Fig. 6.The locally measured angular frequencies ω r , ω θ , ω φ and ω L are connected to the angular frequencies measured by the static distant observers, Ω β , by the gravi- Fig. 15 In-equatorial radial profiles of the frequencies of small harmonic oscillations νr, ν φ , and ν L of charged particle orbiting BH of the mass of M = 10M in external parabolic MF for different values of the parameters B and w.
tational redshift transformation where (f (r)/E(r)) is the redshift coefficient, given by the function f (r) and the particle specific energy at the circular orbit E(r).If the fundamental frequencies of the small harmonic oscillations related to the distant observers, Ω β , are expressed in the physical units, their dimensionless form has to be extended by the factor c 3 /GM .Then the frequencies of the charged particle radial and latitudinal harmonic oscillations measured by a distant observer in physical units are given by The behaviour of the fundamental frequencies ω r , ω θ , ω φ and their ratios can help us to distinguish different shapes of charged particle epicyclic orbits in the vicinity of a stable circular orbit.
According to the Newtonian theory of gravitation, all the frequencies are equal, ω r = ω θ = ω φ , giving ellipse as the only possible bounded trajectory of a test particle around a gravitating spherically symmetric body.For uncharged particles moving around a Schwarzschild BH, the relation ω r < ω θ = ω φ holds, and there exists a periapsis shift for bounded elliptic-like trajectories implying the effect of relativistic precession that increases with decreasing radius of the orbit as the strong gravity region is entered [50].For charged particle frequencies around BH in presence of MF, one gets in general ω r = ω θ = ω φ , see [25,34,51].
The radial ν r , vertical ν θ , orbital ν φ and Larmor ν L frequencies measured by a distant observer in offequatorial plane (B < 0) for different MF configurations as a function of radial coordinate r are plotted in Fig. 14.We see that orbital ν φ and vertical ν θ frequencies coincide for all considered MF configurations.When the attractive LF is weak, the Larmor frequency is smaller than the orbital frequency, however, it coin-Fig.16 Power Spectral Density for three r(τ ) (blue), θ(τ ) (red) and φ(τ ) (green) coordinates along the particle's trajectory.Fundamental frequencies calculated analytically for circular orbit and presented in Figs. 14 and 15 are plotted as vertical dashed lines.The top row of figures corresponds to trajectories without RR k = 0, see, also Figs. 6 and 7. Cases 3 and 4 represent regular trajectories around off-and in-equatorial plane minima with distinguishable main peaks; The case 2 corresponds to the chaotic trajectory.Second row of figures is plotted for radiating particles, where the higher frequency peaks (overtones) are flattened due to radiation damping and the spectrum is more smooth (noisy) in higher frequencies.
cides with the orbital frequency as the magnetic parameter B is increased.For parabolic MF (w = 3/4), we can see that the radial frequency ν r is changing dramatically with the parameter B, while the vertical ν θ and orbital ν φ frequencies remain almost the same.For strong attractive LF (B < −1), the charged particle moves relatively slowly around BH at the off-equatorial plane minima, while oscillating relatively rapidly around MF lines.The radial frequency is too small for the case of paraboloidal MF with w = 1, so the radial profiles start to appear only far away from the BH, which is in concord with the results presented for offISCO in Fig. 3.
Fundamental frequencies for the case of inequatorial plane (B > 0) are plotted in Fig. 15.When the repulsive LF is weak, the Larmor frequency ν L is the smallest among all the frequencies, while the radial ν r and orbital ν φ frequencies coincide.However, when the repulsive LF is strong (B > 1), the radial ν r and Larmor ν L frequencies have different profiles in strong gravity regime close to the BH, but the frequencies ν r and ν L will come closer as the radial distance increases.
The effect of RR on Fourier spectra (Power Spectral Density) calculated from charged particle trajectories is shown in Fig. 16.Particles moving close to circular orbits have main peaks at fundamental frequencies given by Eqs. ( 43) -( 46).The higher harmonics (overtones) are also shown.The RR damping force does not change main frequency peaks, but it lowers the higher harmonics making the spectrum more flat and diluted in high frequency range (see cases 3 and 4 in Fig. 16 in comparison with the cases 3R and 4R, respectively).The RR makes the motion more irregular (chaotic) on short time scales.For both attractive (B < 0 case 3) and repulsive (B > 0 case 4) LFs, the orbital motion (φ, green) is considerably suppressed in high frequency range.We interpret this effect as follows: since the azimuthal motion (in φ direction) in quasi-circular orbits has the highest curvature, hence larger acceleration with respect to other directions of the motion, this component of the motion radiates away faster.Due to RR the trajectory becomes more smooth on short time scales, since the high curvature modes are dumped.Similar effect of RR damping can also be seen for chaotic trajectory (case 2)  the world constants in Eq. ( 17), the dimensionless parameter B, widely used in the present paper, takes the form reflecting the relative influence of the gravitational and MFs on the charged particle motion.As an example, let us estimate B for an electron around Galactic centre SMBH Sgr A*, which is currently one of the best known SMBH candidates.The equipartition value of the MF strength at the event horizon scale of Sgr A* is usually estimated to be around 10 − 100 Gs [1,55].Mass of the BH is measured to be around 4×10 6 M (see, e.g.[56]).Thus, we get for electrons at the Galactic centre For protons, the values of B SgrA * ( p + ) is of the order of 10 7 .In Table 2 we show the values of B for electron, proton, iron atom (without one electron), and a dust particle (without one electron) for various values of the MF strength and the three BH types (stellar mass, supermassive, and Galactic centre BH).Another robust way of estimating the MF in nearby the horizon region could be from observed relativistic jet, which is believed to be launched by the Blandford-Znajek mechanism [17].RR force acting on a charged particle is represented by the dimensionless parameter k which has the form The value of parameter k is much lower than that of B.
For protons, the values of k parameter is lower by the factor of m p /m e ≈ 1836, as in the case with B. Despite the weakness of parameter k as compared to B, it enters into the equations for ultrarelativistic particles as kB 2 , which can make the effect of the RR force considerably large.These values for RR have been estimated for Sgr A* and M87 in [17] and in more details for M87 in [57].
The motion of radiating charged particle is described by Eq. ( 32), containing several terms of different magnitudes.Let us estimate the contribution of each of the term in the equation separately and compare the individual terms with the gravitational force.The left hand side of ( 32) is purely gravitational term given by the Christoffel symbols, which can be taken in the dimensionless units to be ∼ 1.On the right hand side (32) one can distinguish the LF ∼ B and the contribution due to the RR consisting from the terms, which are proportional to ∼ kB and ∼ k 2 B 2 .Representative values for all four terms forces are given in Table 3.Here we assume the particles to be relativistic v ∼ 1, which allows us to use the magnitude estimates for distinguishing different forces of particle dynamics, although the exact values of individual forces at given point may differ.Large values of B in astrophysical settings suggest that the effects of EM interaction and RR can not be neglected since the LF and RR can be much stronger than gravity, but RR can overcome LF only for particles with large specific charge (electron) and for very strong MF of the order of B > 10 15 Gs.
The RR timescale can be estimated from the following equation [28] τ RR ≈ 3 2 where f (r) is a lapse function.Typical synchrotron RR timescales of particles in various MFs are given in Table 4.For the Galactic centre SMBH, the electron cooling timescale is of the order of 10 4 s, while orbiting timescale at inISCO is ∼ 10 3 s.In general, for ions Table 5 Magnitudes of some relevant astrophysical MFs and their related Larmor frequencies.
the cooling timescale is shorter than in case of electrons by the factor of (m p /m e ) 3 ∼ 10 10 .As one can see from the Table 4, the synchrotron radiation losses of ions and heavier particles are relevant in the presence of relatively strong MFs, i.e., much stronger than those corresponding to BHs, e.g., in the vicinity of neutron stars.
Fundamental frequencies for oscillating charged particle derived analytically in previous section can be relevant in astrophysical settings, e.g. in relation to the modulation of intensity flux (light curve) in quasiperiodic oscillations [25,51], observed in many BH systems or in the model of the Galactic centre flares components [58].Fundamental frequencies can also be closely related to frequencies of radiated EM radiation [30], hence they manifest themselves in produced EM spectrum and EM polarization.For strong MFs, the spectrum will be dominated by the Larmor frequency, estimated in Table 5 for some typical astrophysical MFs.As one can see that for typical AGNs (10 4 Gs) the Larmor frequency is of the order of hundreds of GHz, similar to the observational frequency of the Event Horizon Telescope (230 GHz).
The synchrotron radiation losses of charged particles can be relevant in a plasma surrounding BH, if the cooling timescale is shorter or comparable with the mean free time between collisions of particles in plasma.The later can be estimated from the equation [59] where lnΛ is a Coulomb logarithm, given by lnΛ = ln 3 2e 3 Note that τ ee represents the collisional timescale between electron-electron collisions in an electron-proton plasma.For proton-proton and electron-proton collisions at the same temperature, the timescales are typically longer.The representative estimates of collision timescales for given densities and temperatures of plasma are given in Table 6.Table 6 Collision timescale of electron-electron τee and ionion τee collisions in a plasma of given particle density and temperature.The plasma is considered quasi-neutral with equal temperatures of the components.The values are rounded to the orders of magnitude.
It is interesting to estimate the electron collision timescales for the environment of the Galactic centre SMBH.A particle density in a plasma surrounding Sgr A* can be obtained from the observed flares properties [58], to be of the order of n ∼ 10 7±1 cm −3 .Modeling the emission of simultaneous NIR and X-ray flares from Sgr A* the Lorentz factor of electrons are estimated to be of the order of γ e ∼ 10 3 , which corresponds to the plasma with electron temperature T e ∼ 10 9 K.Then, the electron-electron collision timescale is τ ee (Sgr A*)= 10 4±1 s.Note that this estimate is based on the flare observations, i.e., considered denser regions of the Galactic centre compared to its average surrounding.Such estimate depends on a number of assumptions -it may very well be that the density in a flaring region is lower, as long as the magnetization is high enough to power particles to Lorentz factors of ∼ 10 3 [60].At larger spatial scales, averaging the plasma environment, the collision timescale is less by a few orders of magnitude [61].Thus, the plasma surrounding Galactic centre can be considered as collisionless, so the considerations and estimates given in the current paper are fully applicable.In principle, there is more magnetic energy available in dilute (high sigma) regions, so flares could be powered from there.There are simulations showing that flares from Sgr A* can be powered from magnetized regions, see [17,60].This however does not change the main conclusion that the plasma is collisionless.

Conclusions
The effects of MFs around BHs cannot be neglected in many astrophysical settings.If the specific charge of a test particle q/m is large enough, as in the case of elementary particles, even a weak MF can significantly influence the particle dynamics, the position of the inner edge of an ionized Keplerian accretion disk, or the epicyclic frequencies for motion around BH.In this article we have tested a simple but realistic model for BH magnetosphere given by the parabolic MF and ex-plored influence of RR damping force on chaotic charge particle dynamics.We also briefly discussed the astrophysical relevance of the presented theoretical model.
We studied the parabolic MF configuration around BH in general form for different values of the inclination parameter w, limiting to the magnetic split-monopole solution for w = 0 as a special case.The parabolic configuration of MF is supported by observations and GRMHD numerical experiments [9].Generation of such MF geometry would require the presence of electric current sheets at the equatorial plane provided by ongoing accretion processes.The accretion disk is located in equatorial plane and outside of the equatorial plane, the matter density is expected to be much lower.Therefore, in discrete particle approach used in the present paper, the existence of charged particle off-equatorial orbits become especially important.
We have analysed stability of circular orbits of charged particles and found that the off-equatorial stable circular orbits may exist only in the case of attractive LF (B < 0), while the stable circular orbits at the equatorial plane can be observed only in the case of repulsive LF (B > 0).The effective potential minima situated at the off-equatorial plane get closer to the z-axis when the MF becomes stronger.In the case of the in-equatorial plane, the radius of inISCO of the charged particle decreases as the MF strength increases, approaching the BH horizon r inISCO = 2 at large values of the magnetic parameter B. In the case of the split magnetic monopole, the inISCO (or offISCO) remains at a constant radius r inISCO = 6.
We examined the regular and chaotic behaviour of the motion of charged particles with the help of PSs, which is a good tool illustrating the difference between regular and chaotic dynamics.We have found that charged particle motion around BH in parabolic MF can be chaotic under certain initial conditions.This result is in contrast with the neutral particle case or the case with the monopole MF configuration, where due to the symmetry of the EM and gravitational fields, chaos does not occur and the motion is always regular [24].We have also compared the change of the chaoticity of trajectories when RR force in included in the equations of motion and found that the RR acts as a damping force, creating new attractors in phase space, which are situated on circular orbit at the equatorial plane of the BH.We have found that the off-equatorial plane orbits are always unstable when the radiation is taken into account, since the RR force tends to shift the charged particle orbits to their guiding centers, which coincides with the position of the BH.The RR force has damping character of the influence on the charged particle dynamics and the final stage of the radiating particle motion is eventually a collapse into the BH if the integration time is long enough.The particle loses its energy and angular momentum due to RR cooling, which leads to spiraling down to the BH.For positive value of magnetic parameter B > 0, an oscillating charged particle experiencing the RR force radiates its oscillations away and settles down to the quasi-circular orbits, before eventually falling down into the BH.
We have also studied the influence of the RR force on the structure of a thin Keplerian accretion disk orbiting Schwarzschild BH in external parabolic MF.The effect of RR can be observed even for a relatively small values of B. For |B| 1, one can observe oscillations in the disk dynamics, for |B| ∼ 1, the disk trajectories enter into the turbulent chaotic regime and for |B| 1, the charged particle trajectories stick to MF lines.Due to RR force, disk loses energy and angular momentum and particle will end up into BH.The time scale of the destruction of the Keplerian disk by radiation strongly depends on the strength of MF and radiation parameter kB.
We have calculated the fundamental frequencies of a charged particle and produced EM spectrum, which can be useful in determining time scales of various astrophysical processes around BH and fitting observed data.Test particle frequencies are fundamental in the sense that they will appear in any more advanced and complicated model for BH neighborhood description, like GRMHD simulation as a zero approximation [62].It is worthwhile mentioning that our analytical solutions are very important to test the codes that evolve particles (e.g.GRPIC), for example in the codes of [63][64][65].One example where particle frequencies could be crucial is the quasi-periodic oscillations of the X-ray power density observed in quasars and microquasars.Fitting these quasi-periodic oscillation frequencies in timing spectra can help to verify BH magnetosphere model [51,66].Another observational window for BH accretion processes could be non-thermal synchrotron emission produced by charged particles moving in the BH magnetosphere [67].One charged particle dynamics and its oscillations around effective potential minima could be used to construct synthetic charged particle synchrotron spectra [57,68,69] We have also estimated the most relevant parameters of the discussed model in astrophysical situations and found the limits on the relevance of the effects of the LF and RR for different particle types and BH systems.For that, we gave the estimates of orbital, radiation, and collisional timescales of particles.We have shown that applied to the Galactic centre BH the plasma environment can be considered as collisionless and test particle approach can be used there.We also discussed the role of MF in fundamental frequencies, observable in relation to the so-called quasi-periodic oscillations.
Among the possible extensions of the presented parabolic magnetosphere model, which could be interesting to explore in future are the calculation of the radiated EM spectra and polarization, and inclusion of the BH rotation into the model.In relation to the existence of the off-equatorial orbits found in this paper, it would also be interesting to explore, whether the stable off-equatorial structures, existing in discrete particle approach can be found in numerical GRMHD or PIC simulations.This is a prediction that needs to be tested.

Fig. 1
Fig.1Parabolic MF configuration around Schwarzschild BH.Quasi neutral Keplerian accretion disk formed by particles on circular orbits close to the equatorial plane is restricted in its inner edge by the equatorial innermost stable circular orbit (inISCO).The dashed curve is the boundary between the accretion disk corona and jet funnel region close to the z-axis.

Fig. 2
Fig.2Effective potential V eff for charged particle dynamics as 2D function of Cartesian x, z coordinates.We can see the effective potential minima at the off-equatorial planes for B < 0 (top left) and the in-equatorial minima for B > 0 (top right).The lower row of figures shows the sections of the effective potential at given energy and the corresponding minima.Value of the parameter w has been set to w = 3/4, but similar behavior is observed for different values of w ∈ (0, 1).

Fig. 5
Fig.5Energy in effective potential minima, for off-equatorial (B < 0) and in-equatorial (B > 0) plane.We set the inclination parameter of the MF to w = 3/4.The black thick curve is plotted for neutral particles.Each point of the dashed curve represents the position of offISCO (or inISCO), while the numbers on solid curves indicate the values of magnetic parameter B.

Fig. 6
Fig. 6 Charged particles trajectories (blue) around BH with parabolic MF, representing two astrophysically relevant cases: when MF influence is negligible (|B| 1, first row), and when LF is dominating over gravity (|B| 1, second row).The gray circle represents the BH, gray curves MF lines, the dashed curve denotes the energy boundary, and the black dot the initial position of the particle.

Fig. 11
Fig. 11 Example of ionization mechanism for neutral particles in Keplerian accretion disk.The LF and RR force acting on electrons and protons (ions) are significantly different.In the presented case, the proton (blue) is oscillating around effective potential minima in the equatorial plane and influenced by LF only slightly (B = 0.01).On the other hand, the electron (red) is strongly influenced by the LF (B = −10) in the same MF, orbiting MF lines, while moving up, down, and around the BH.

Fig. 12
Fig. 12 Evolution of thin Keplerian accretion disks around Schwarzschild BH immersed into external parabolic MF.Initially, the particles from accretion disk follow circular geodesics with slight inclination from the equatorial plane (θ 0 = 1.5), while the MF is directed along the z-axis.We set the declination parameter of the parabolic MF lines to w = 3/4, radiation parameter to k = 10 −5 , and the values for the MF parameter to B = {±10, ±1, ±0.1, ±0.01, 0}.The black dashed closed curves are the ISCOs for neutral test particles, indicating the inner edge of the Keplerian accretion disk.In the absence of the MF (B = 0 case, middle figure), all the orbits remain in their circular shape and only the inclined razor-thin disk can be seen.If a slight EM interaction is switched-on (B = ±0.01cases, middle row), the charged particles creating a disk that is initially almost perpendicular to the MF lines start to oscillate with epicyclic frequencies around the circular orbit in both radial and vertical directions; the accretion disk becomes slightly thicker.Further, increasing the MF parameter (B = ±0.1,±1 cases), the motion of charged particles becomes chaotic and the accretion disk is either destroyed or modified into a thick toroidal structure.For B = −1 case, all the particles are captured by the BH, and complete destruction of the Keplerian disk can be observed.The LF dominates the particle motion when the parameter of the MF is large enough (|B| ≥ 10 cases).The charged particles spiral down and up along the MF lines, and move around the BH in the counter-clockwise (B < 0) and clockwise (B > 0) directions.

Table 3
Magnitudes of different forces acting on radiating charged particle in the dimensionless form of Eq. (32) for different values of magnetic field strengths, compared to the gravitational force of the order of ∼ 1.The estimates correspond to a relativistic electron in vicinity of a stellar mass BH with M = 10 M .

Table 4
Typical cooling times of electron τe, proton τp, and ionized iron atom (without one electron) for different MF strengths.For charged dust particles, the synchrotron cooling is irrelevant.