Innermost stable circular orbit near dirty black holes in magnetic field and ultra-high-energy particle collisions

We consider the behavior of the innermost stable circular orbit (ISCO) in the magnetic field near “dirty” (surrounded by matter) axially symmetric black holes. The cases of near-extremal, extremal, and nonextremal black holes are analyzed. For nonrotating black holes, in the strong magnetic field ISCO approaches the horizon (when backreaction of the field on the geometry is neglected). Rotation destroys this phenomenon. The angular momentum and radius of ISCO look model-independent in the main approximation. We also study the collisions between two particles that results in the ultra-high energy Ec.m.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{\text {c.m.}}$$\end{document} in the center-of-mass frame. Two scenarios are considered—when one particle moves on the near-horizon ISCO or when collision occurs on the horizon, one particle having the energy and angular momentum typical of ISCO. If the magnetic field is strong enough and a black hole is slowly rotating, Ec.m.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{\text {c.m.}}$$\end{document} can become arbitrarily large. The kinematics of the high-energy collision is discussed. As an example, we consider the magnetized Schwarzschild black hole for an arbitrary strength of the field (the Ernst solution). It is shown that backreaction of the magnetic field on the geometry can bound the growth of Ec.m.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{\text {c.m.}}$$\end{document}.


Introduction
The motion of particles in the vicinity of black holes is a subject that has been continuing to attract interest until recently. In doing so, a special role is played by circular orbits-see, e.g., Refs. [1,2] and references therein. Especially, this concerns an innermost stable circular orbit (ISCO). It is important in phenomena connected with accretion disc and properties of cosmic plasma [3,4]. Apart from astrophysics, such a kind of orbits possesses a number of nontrivial features and, therefore, is interesting from the theoretical viewpoint. In a classic paper [5] it was shown that in the extremal limit ISCO a e-mail: zaslav@ukr.net approaches the horizon. As a result, some subtleties arise here since the horizon is a lightlike surface, so a massive particle cannot lie within it exactly. Nowadays, near-horizon circular orbits for near-extremal and extremal rotating black holes are still a subject of debates [6][7][8][9].
Quite recently, a new circumstance came into play that makes the properties of ISCO important in a new context. Namely, it is the ISCO that turns out to be a natural venue for the realization of the so-called BSW effect. Several years ago, it was shown by Bañados, Silk, and West that if two particles collide near the black hole horizon of the extremal Kerr metric, their energy E c.m. in the center-of-mass (CM) frame can grow unboundedly [10]. These findings stimulated further study of high-energy collisions near black holes. The validity of the BSW effect was extended to extremal and nonextremal more general black holes. It was also found that there exists a version of this effect near nonrotating electrically charged black holes [11]. Another version of ultra-high-energy collisions reveals itself in the magnetic field, even if a black hole is neutral, vacuum, and nonrotating, so it is described by the Schwarzschild metric [12]. Generalization to the case when the background is described by the Kerr metric was done in [13].
In the BSW effect, one of the colliding particles should be so-called critical. It means that the energy and the angular momentum (or electric charge) of this particle should be finetuned. In particular, the corresponding critical condition is realized with good accuracy if a particle moves on a circular orbit close to the horizon. Therefore, an ISCO can play a special role in ultra-high-energy collisions in astrophysical conditions. Without the magnetic field, this was considered in [14] for the Kerr black hole and in [15] for more general rotating black holes. Kinematically, the effect is achieved due to collision of a rapid typical so-called usual particle (without fine-tuning) and the slow fine-tuned particle on the ISCO [16] (see also below).
In [12,13] collisions were studied just near the ISCO in the magnetic field. In both cases, a black hole was taken to be a vacuum one. Meanwhile, in astrophysical conditions, black holes are surrounded by matter. By definition, such black holes are called "dirty", according to the terminology suggested in Ref. [17]. (We would like to stress that it is matter but not the electromagnetic field that makes a black hole dirty.) The aim of our work is two-fold, since two different issues overlap here. The first one concerns the properties of ISCO near dirty black holes in a magnetic field, so both matter and the magnetic field are present. The second issue concerns the scenarios of high-energy particle collisions near such orbits. We derive general asymptotic formulas for the position of the ISCO in the magnetic field, which are used further for the evaluation of E c.m. and examining of two scenarios of the BSW effect near ISCO.
In Refs. [12,13], it was assumed that the magnetic field is weak in the sense that backreaction of the magnetic field on the metric is negligible but, at the same time, it is strong in the sense that it affects he motion of test particles. Such a combination is self-consistent, since the dimensionless parameter b that controls the magnetic field strength contains a large factor q/m, relevant for the motion of particles. Our approach is model-independent and is not restricted by some explicit background metric. Therefore, most of the formulas apply also to the metrics which are affected by the magnetic field. On the other hand, if the magnetic field is too strong, its backreaction on the metric can change the properties of E c.m. itself, as will be seen below. Thus we discuss two new features absent from previous works in the sense that both matter and magnetic field are taken into account in a modelindependent way.
It is worth noting that high-energy collisions in the magnetic field were studied also in another context, including scenarios not connected with ISCO-see [18][19][20][21].
In general, it is hard to find and analyze ISCO even in the Kerr or Kerr-Newman cases [1,2]. However, it is the proximity to the horizon that enables us to describe some properties of ISCO, even not specifying the metric (so we work in a model-independent way) and even with a magnetic field. This can be considered as one manifestation of universality of black hole physics.
The paper is organized as follows. In Sect. 2, the metric and equations of motion are presented. In Sect. 3, we give basic equations that determine ISCO. In Sect. 4, we consider ISCO in a magnetic field for near-extremal black holes and analyze the cases of small and large fields. In Sect. 5 the cases of nonrotating (but dirty) and slowly rotating black holes are discussed. As we have two small parameters (slow rotation and inverse field strength), we consider different relations between them separately. In Sect. 6 we show that if a black hole rotates, even in the limit of a strong magnetic field ISCO does not tend to the horizon radius. In Sect. 7, it is shown that for extremal nonrotating black holes, for large b, ISCO approaches the horizon radius. In Sect. 8, it is shown that this property is destroyed by rotation. In Sect. 9, general formulas for E c.m. for particle collisions in the magnetic field are given. In Sect. 10, we find the velocity of a particle on ISCO and argue that a kinematic explanation of high-energy collisions is similar to that for the BSW effect [16]. In Sect. 11, we apply general formulas of collision to different black hole configurations and different scenarios. In Sect. 12, the exact solution of the Einstein-Maxwell equations (static Ernst black hole) is chosen as a background for collisions. This enables us to evaluate the role of backreaction of the magnetic field on E c.m. . In Sect. 13, the main results are summarized. Some technical points connected with cumbersome formulas are put in the appendix.
Throughout the paper we use units in which the fundamental constants are G = c = 1.

Metric and equations of motion
Let us consider the metric of the form where the metric coefficients do not depend on t and φ. The horizon corresponds to N = 0. We also assume that there is an electromagnetic field with the four-vector A μ where the only nonvanishing component equals In vacuum, this is an exact solution with B = const [22]. We consider configuration with matter (in this sense a black hole is "dirty"), so in general B may depend on r and θ .
Let us consider motion of test particles in this background. The kinematic momentum p μ = mu μ , where m is the particle's mass, the four-velocity u μ = dx μ dτ , where τ is the proper time, x μ are coordinates. Then the generalized momentum is equal to q is the particle's electric charge. Due to the symmetry of the metric, P 0 = −E and P φ = L are conserved, where E is the energy, and L is the angular momentum. We consider motion constrained within the equatorial plane, so θ = π 2 . Redefining the radial coordinate r → ρ, we can always achieve that within this plane. Then the equations of motion givė ( 1 0 ) A dot denotes differentiation with respect to the proper time τ . As usual, we assume the forward in time conditionṫ > 0, so X ≥ 0. Hereafter, we use the notations The subscripts "+", "0" denote quantities calculated on the horizon and ISCO, respectively.
In what follows, we will use the Taylor expansion of the quantity ω near the horizon. We denote x = ρ − ρ + , where ρ + is the horizon radius. Then
In general, it is impossible to find exact solutions of Eqs. (17) and (18). Therefore, in the next sections we analyze separately different particular situations, with main emphasis made on the near-horizon region. In doing so, we develop different versions of the perturbation theory that generalize the ones of [13]. The radius of ISCO, its energy, and angular momentum are represented as some series with respect to the corresponding small parameter, truncated at the leading or subleading terms, similarly to [13].

Near-extremal black holes
Let us consider a nonextremal black hole. In what follows, we are interested in the immediate vicinity of the horizon and use the Taylor series for corresponding quantities. Then near the horizon we have the expansion where κ has the meaning of the surface gravity. By definition, we call a black hole near-extremal if where x 0 = ρ 0 − ρ + . Then for the lapse function we have an expansion near ISCO, Taking into account (18), after straightforward (but somewhat cumbersome) calculations, one can find that

)
where P and d P dx are to be taken at x = x 0 or, with the same accuracy, at x = 0 (i.e., on the horizon). Then From (16) and (21) we have Using (23), (24), and (9) we derive an equation for the value of the angular momentum L 0 on ISCO: where To have a well- We are interested in the positive root according to (23). Then and, in a given approximation, where we neglected the difference between dβ dx + and dβ dx 0 . Equations (33)-(37) give the expression for H after substitution into (28). To avoid cumbersome expressions, we leave it in the implicit form. Now, two different limiting cases can be considered.

Small magnetic field
which agrees with Eq. (44) of Ref. [15]. It follows from (16) and (30) that Let us consider a small but nonzero b. We can find from (33) that

Large magnetic field
According to (27), there exists a finite lim B→∞ H = H ∞ . In doing so, we find from (33), (9), and (30) Thus according to (27), in general the radius of ISCO depends on the value of the magnetic field via the coefficient H . However, there is an exception. Let Then so the dependence on b drops out from the quantity H and, correspondingly, from the ISCO radius (27). One can check easily that the conditions (49) are satisfied for the nearextremal Kerr metric in the magnetic field. This agrees with Eq. (38) of [13] where the observation was made that in the main corrections of the order κ 2/3 the magnetic field does not show up. Thus this is the point where dirty black holes behave qualitatively differently from the Kerr metric in that the dependence of the ISCO radius on b is much stronger than in the Kerr case. It is instructive to evaluate the relation between H (0) and H (∞) for vanishing and large magnetic fields, which results, according to (27), in different values of the corresponding radii x 0 . The dependence on the magnetic field is due to the term 1 P d P dx in the denominator, One can find that However, in general they can differ significantly. Say, for As a result, the ISCO radius (27) also may vary over wide range.

Slowly rotating black hole
Now, we assume that κ is not small, so the first term in (19) dominates. Here, we will consider different cases separately.

Nonrotating black hole
Here, we generalize the results known for the Schwarzschild black hole [25,26] to a more general metric of a dirty static black hole. In Eqs. (14) and (15) we should put a 1 = 0 = a 2 . For a finite value of the magnetic field parameter b, ISCO lies at some finite distance from the horizon. However, now we will show that in the limit b → ∞, the radius of ISCO tends to that of the horizon with x 0 ∼ b −1 .
We will show that this indeed happens, provided the term with L in (9) is large and compensates the second one with b. Correspondingly, we write where For what follows, we introduce the quantity (1). Then, near the horizon, where x is small, we can use the Taylor expansion where Now, β 0 1 but, by assumption, β is finite. In terms of the variable it can be rewritten as It is clear from the above formulas that the expansion with respect to the coordinate x is equivalent to the expansion with respect to inverse powers of the magnetic field b −1 , so for b 1 this procedure is reasonable. Then, after substitution of (62), we can represent (14) and (15) in the form of an expansion with respect to β −1 0 : Here, the coefficients at leading powers are equal to Then, in the main approximation, we have the equations C 0 = 0 and S 1 = 0 which give us Equations (69) and (71) agree with [12,13]. It is interesting that in terms of the variables u, L 0 R + , and b the result (68) looks model-independent in the main approximation. This can be thought of as a manifestation of the universality of black hole physics near the horizon. The dependence on a model reveals itself in higher-order corrections.

Extremely slow rotation
Now, we consider rotation as the perturbation. Here, the angular velocity of rotation is the smallest parameter. Correspondingly, in Eqs. (18) and (15) we neglect the term L 2 , since it contains ω 2 . More precisely, we assume so from (184) and (185) we have In the particular case of the slowly rotating Kerr metric, There are two kinds of corrections-due to the magnetic field and due to rotation. One can check that the presence of rotation leads to the appearance in the series (64) and (65) of half-integer inverse powers of β 0 , in addition to integer ones. In the main approximation, we consider both kinds of corrections as additive contributions. Omitting details, we list the results: It follows from (8) and (16) that For the slowly rotating Kerr metric, R + ≈ 2M, In the main approximation the difference between the Boyer-Lindquist coordinate r and the quasiglobal one ρ has the same order a 2 and can be neglected. Then They agree with the results of Sect. 3 B 2 of [13]. It is seen from (82)-(84) that the fractional corrections have the order a * b 3/2 and are small in accordance with (74). In the more general case, the small parameter of the expansion corresponds to (73), so it is the quantity

Modestly slow rotation
Let now, instead of (72) and (73), the opposite inequalities hold: It turns out (see the details in the appendix) that It follows from (190), (191), and (199) that

Kerr metric
In the case of the slowly rotating Kerr black hole, Eq. (81) entails where we used (96). One should compare this result to that in [13]. Now, R + = 2M, the horizon radius of the Kerr metric r + ≈ 2M 1− a * 2 4 . Equation (53) of [13] gives us whence x 0 = r 0 −r + ≈ R + 16 a * 2 , which coincides with (99). It is seen from (91) and (98) that the angular momentum takes the value which coincides with Eq. (55) of [13]. Also, one finds that In Eq. (90) one should take into account that ω + depends on r + , which itself can be expressed in terms of a * and M.

ISCO for rotating nonextremal black holes in a strong magnetic field
In the previous section we saw that in the limit b → ∞ the ISCO radius does not coincide with that of the horizon, which generalizes the corresponding observation made in Sect. III B 3 of [13]. Now, we will see that this is a general result which is valid for an arbitrary degree of rotation and finite κ (so, for generic nonextremal black holes). It is worth noting that for b = 0 it was noticed that the near-horizon ISCO are absent [15,27]. However, for b 1 the corresponding reasonings do not apply, so we must consider this issue anew.
We have to analyze Eqs. (17) and (18) in which (16) is taken into account.
Neglecting higher-order corrections, we can rewrite them in the form (1) Let us suppose that β is finite or, at least, β b. Then it follows from (57), (9), and (55) that dβ dx ∼ b and L ∼ b. Also, x ∼ b −1 according to (61), N ∼ √ x ∼ b −1/2 . However, it is impossible to compensate the term with L 2 in (105) having the order b 2 .
(2) Let β ∼ L ∼ b. Then in (104) the first term has the order b 2 and cannot be compensated. Thus we see that, indeed, in the limit b → ∞ the assumption about x → 0 leads to contradictions, so the ISCO radius does not approach the horizon.

Extremal nonrotating black hole
Up to now, we considered the case of a nonextremal black hole, so the surface gravity κ was arbitrary or a small quantity but it was nonzero anyway. Let us discuss now the case of the extremal black hole, so κ = 0 exactly. We pose the question: is it possible to get an ISCO such that for b → ∞ the ISCO radius tends to that of the horizon? Now, we will see that this is indeed possible for a nonrotating black hole (ω = 0). We assume that the electric charge that can affect the metric is negligible. The extremal horizon appears due to the properties of the matter that surrounds the horizon, which is possible even in the absence of the electric charge, provided the equation of state obeys some special conditions [23,24].
For ISCO close to the horizon we can use the expansion in which we drop the terms of the order x 3 and higher. Now we show that the case under discussion does exist with a finite quantity β. We can use now (62) in which only the first term is retained, so where u is given by Eq. (61). Then Eq. (10) reads where Equations (14) and (15) reduce to They have the solution Correspondingly, Eqs. (107) and (13) give us We can also find the angular momentum on ISCO Thus for big b there is ISCO outside the horizon that tends to it in the limit b → ∞, so that the quantity x 0 → 0.

Extremal rotating black hole
Now, we consider the same question but now for rotating black holes: is it possible to have ISCO in the near-horizon region (as closely as we like) for the extremal BH, when κ = 0? Mathematically, it would mean that Then Eqs. (17) and (21) with κ = 0 give for small x Equation (18) with terms of the order x 2 0 and higher neglected, gives rise to Then the main terms in (119) and (120) entail For b 1, assuming for definiteness that d > D (d is defined according to (32)), one finds from (9) and (121) that The terms x 2 0 in (119) and x 0 in (120) give us, with ( 121) taken into account The system is overdetermined; Eq. (127) cannot be satisfied in general. In principle, one can consider (127) as a restriction on the black hole parameters. This is similar to the situation for the extremal Kerr-Newman metric (b = 0), where ISCO near the horizon exists only for the selected value of the angular momentum, approximately equal to a M ≈ 1 √ 2 [8,9]. However, we will not discuss such exceptional cases further. Generically, the answer to our question is negative, so the ISCO radius does not approach the horizon in the limit b → ∞.

Particle collisions: general formulas
Let two particles collide. We label their characteristics by indices 1 and 2. Then, in the point of collision, one can define the energy in the center of mass (CM) frame as Here, is the total momentum, and is the Lorentz factor of their relative motion. For motion in the equatorial plane in the external magnetic field (2), one finds from the equations of motion (6) and (7) that Here, ε = +1, if the particle moves away from the horizon, and ε = −1 if it moves toward it. Now, there are two scenarios relevant in our context. We call them the O-scenario and H-scenario according to the terminology of [15]. Correspondingly, we will use the superscripts "O" and "H".

O-scenario
Particle 1 moves on ISCO. As V 1 (ρ 0 ) = 0 on ISCO, the formula simplifies to As we are interested in the possibility to get γ as large as one likes, we will consider the case when the ISCO is close to the horizon, so N is small. In doing so, we will assume that (X 2 ) = 0, so particle 2 is usual according to the terminology of [15]. We also must take into account Eq. (16), whence For ISCO close to the horizon, the first term dominates and we have 9.2 H-scenario Now, particle 1 leaves ISCO (say, due to additional collision) with the corresponding energy E = E(x 0 ) and angular momentum L = L(x 0 ), which corresponds just to ISCO. This particle moves toward the horizon where it collides with particle 2.
Mathematically, it means that we should take the horizon limit N → 0 first in Eq. (131). We assume that both particles move toward the horizon, so ε 1 ε 2 = +1. Then where all quantities are to be calculated on the horizon. For small X 1 , when we see from (135) that Now, where X 0 = E 0 − ω 0 L corresponds to ISCO. With (12) and (16) taken into account, in the main approximation Now, we apply these formulas to the different cases considered above.

Kinematics of motion on ISCO
It is instructive to recall that the general explanation of high E c.m. consists in the simple fact that a rapid usual particle, having a velocity close to speed of light, hits the slow particle that has parameters approximately equal to the critical values. This was explained in detail in [16] for the standard BSW effect (without considering collision near ISCO). Does this explanation retain its validity in the present case? One particle that participates in a collision is usual, so it would cross the horizon with a velocity approaching the speed of light in an appropriate stationary frame (see below). We consider the near-horizon ISCO, so the velocity of a usual particle is close to the speed of light. Now, we must check what happens to the velocity of a particle on ISCO.
To describe the kinematic properties, it is convenient to introduce the tetrads that in the local tangent space enable us to use formulas similar to those of special relativity. A natural and simple choice is the tetrad of a so-called zeroangular observer (ZAMO) [5]. It reads N (1, 0, 0, 0), Here, It is also convenient to define the local three-velocity [5] according to From the equations of motion (13)- (22) and the formulas for the tetrad components, we obtain the component v (2) = 0 for equatorial motion. Then, introducing also the absolute value of the velocity v according to one can find that Equation (150) was derived in [16] for the case when the magnetic field is absent. We see that its general form does not depend on the presence of such a field.
For a circle orbit, Eq. (16) should hold. Comparing it with (150), we find that which has the same form as for the static case [12]. Now we can consider different cases depending on the value of the magnetic field and the kind of black hole.
For large b, accoding to (59), the quantity β is proportional to b,and grows, v → 1. However, this case is not very interesting since the individual energy (48) diverges itself.

Nonrotating or slowly rotating nonextremal black holes
According to Eq. (68), β ≈ 1 √ 3 . Slow rotation adds only small corrections to this value. Thus, rather unexpectedly, we again see that on ISCO This value coincides for the near-extremal Kerr without a magnetic field and a nonrotating or slowly rotating dirty black hole in the strong magnetic field.

Modestly rotating nonextremal black hole
It follows from (93) and (88) that |β| 1. However, as the energy of a particle on ISCO (96) tends to inifnity, this case is also not so interesting.
To summarize, in all cases of interest (when an individual energy is finite), β remains finite even in the strong magnetic field. Correspondingly, v < 1 on ISCO and the previous explanation of the high E c.m. [16] applies. For less interesting cases, when an individual energy diverges, we have collision between two rapid particles but their velocities are not parallel and this also gives rise to high γ 0 (see Eq. 20 in [16]).

O-scenario
Using (134), (21) and (29) one obtains In a strong magnetic field, with b 1, using the expression (45) for β, we obtain In the near-extremal Kerr case, As a result, which coincides with Eq. (61) of [13], in which the limit b → ∞ should be taken.

H-scenario
Now, due to (121), Eq. (139) gives X 1 = 0. It means that in the expansion (21) we must retain the first correction in the expression for N , when it is substituted into (139). As a result, we have There are also terms of the order x 2 0 ∼ κ 4/3 but they are negligible as compared to κ. Correspondingly, Eq. (135) gives In a strong magnetic field, with b 1, using (45) again we obtain Thus in both versions, for b 1 the effect is enhanced due to the factor b. For b = 0 we return to [15].
In the Kerr case, which corresponds to Eq. (59) of [13], in which b 1.
11.2 Extremely slowly rotating or nonrotating black hole

H-scenario
Using (78) and neglecting in (139) the second term (rotational part), we get 11.3 Modestly rotating black holes in strong magnetic field

H-scenario
In a similar manner, one can obtain from (137), (139), and (95) that (E H c.m. ) 2 ∼ b with a somewhat cumbersome coefficient, which we omit here.
Both these scenarios are less interesting, since according to (96), the individual energy E 0 ∼ b diverges itself in the limit b → ∞.

Backreaction of magnetic field: Ernst static black hole
Now, we illustrate the obtained results using the metric of a static magnetized black hole [28] that can be considered as the generalization of the Schwarzschild solution. This will also allow us to elucidate the role of backreaction due to the magnetic field on the behavior of E c.m. , which bounds the BSW effect. The metric reads r + = 2M is the horizon radius, B is a constant parameter. It follows from (11) (with B replaced withB) and (167) that Many important details of the particle's motion in this background can be found in Ref. [29].
Calculating the corresponding coefficients according to (58)-(60) and substituting them into (182)-(187), we obtain It follows from (186) that For the energy of collision we have from (160) and (162) where When ξ 1, there is agreement with the results for the Schwarzschild metric [12,13] since Eq. (170) turns into (69) and (172) turns into (70). It is interesting that for any ξ , the velocity of a particle on ISCO is equal to 1/2 like this happens for ξ 1. The approach under discussion works well also for ξ 1, provided that the ISCO lies close to the horizon to ensure large E c.m. , i,e, N 2 1.According to (167) which bounds E c.m. , which begins to decrease when B increases. One should also bear in mind that it is impossible to take the limit ξ → ∞ literally since the geometry becomes singular. In particular, the component of the curvature tensor R θφ θφ grows like ξ 2 . The maximum possible E c.m. is achieved when ξ ∼ 1, then z ∼ q m . The example with the Ernst metric shows that strong backreaction of the magnetic field on the geometry may restrict the growth of E c.m. to such extent that even in spite of large b, the effect disappears because of the factor ξ that enters the metric. It is of interest to consider the exact rotating magnetized black hole [30] that generalizes the Kerr metric but this problem certainly needs separate treatment.

Summary and conclusion
We obtained characteristics of ISCO and the energy in the CM frame in two different situations. For the near-extremal case, we considered the BSW effect. Previous results applied to the weakly magnetized Kerr metric or dirty black holes without the magnetic field. Now, we took into account both factors, so generalized the previous results for the case when both matter and magnetic field are present. In doing so, there is qualitative difference between dirty rotating black holes and the Kerr one. Namely, the radius of ISCO depends on the magnetic field strength b already in the main approximation with respect to small surface gravity κ in contrast to the case of the vacuum metric [13], where this dependence reveals itself in the small corrections only.
For extremal black holes, we showed that, due to the strong magnetic field, there exists the near-horizon ISCO that does not have a counterpart in the absence of this field. Correspondingly, we described the effect of high-energy collisions near these ISCO.
We demonstrated that rotation destroys near-horizon ISCO both for the nonextremal and extremal horizons, so lim b→∞ r 0 (b) = r + . However, if the parameter responsible for rotation is small, E c.m. is large in this limit.
For slowly rotating black holes we analyzed two different regimes of rotation thus having generalized previous results on the Kerr metric [13]. The parameters of expansion used in calculations and the results agree with the Kerr case. In particuar, for modestly slow rotation the individal energy of the particle on ISCO is unbound.
In the main approximation, the expressions for the ISCO radius and angular momentum in dimensionless variables are model-independent, so here one can see universality of black hole physics.
We also found the three-velocity of a particle on ISCO in the ZAMO frame. It turned out that for slowly rotating dirty black holes in the magnetic field it coincides with the value typical of the Kerr metric without a magnetic field, v ≈ 1 2 . Correspondingly, the previous explanation of the high E c.m. as the result of collisions of very fast and slow particles [16] retains its validity in the scenarios under discussion as well.
In previous studies of the BSW effect in the magnetic field [12,13,21], some fixed background was chosen. In this sense, the magnetic field was supposed to be weak in that it did not affect the metric significantly (although it influenced strongly motion of charged particles). Meanwhile, the most part of the formulas obtained in the present work applies to generic background and only asymptotic behavior of the metric near the horizon was used. Therefore, they apply to the backgrounds in which the magnetic field enters the metric itself, with reservation that the surface gravity κ = κ(b), etc. In particular, we considered the static magnetized Ernst black hole and showed that strong backreaction of the magnetic field on the geometry bounds the growth of E c.m. .
Thus we embedded previous scenarios of high-energy collisions in the magnetic field near ISCO in the vicinity of black holes [12,13] and took into account the influence of the magnetic field on the metric.
Throughout the paper, it was assumed that the effect of the electric charge on the metric is negligible. It is of interest to extend the approach of the present work to the case of charged black holes.

Equation (18) reads
where now, with a given accuracy, Substituting L R + ≈ b into (193) and (194 ) and assuming h < 0, one finds the system of two equations . This system can be solved easily. There are two roots here but only one of them satisfies the condition h = 0: (For β > 0, one can obtain the equation 6y − h 1 = −3δ √ y but in combination with (198) it would give y < 0, which is unacceptable since outside the horizon we should have y > 0. Thus this case should be rejected.) Then, using (188), (192), and (12), we find the results (93)-(97).