Vacuum and nonvacuum black holes in a uniform magnetic field

We modify and generalize the known solution for the electromagnetic field when a vacuum, stationary, axisymmetric black hole is immersed in a uniform magnetic field to the case of nonvacuum black holes (of modified gravity) and determine all linear terms of the vector potential in powers of the magnetic field and the rotation parameter.

We modify and generalize the known solution for the electromagnetic field when a vacuum, stationary, axisymmetric black hole is immersed in a uniform magnetic field to the case of nonvacuum black holes (of modified gravity) and determine all linear terms of the vector potential in powers of the magnetic field and the rotation parameter.

I. THE MAGNETIC FIELD PROBLEM
A Killing vector ξ µ in vacuum (no stress-energy T µν ≡ 0) is endowed with the property of being parallel, that is proportional, to some vector potential A µ that solves the source-less (no currents J µ = F µν ;ν = ( |g|F µν ) ,ν / |g| ≡ 0) Maxwell field equations. So, ξ µ is itself a solution to the same source-less Maxwell field equations. In Ref. [1], this property was employed as an ansatz to determine the electromagnetic field of a vacuum, stationary, axisymmetric, asymptotically flat black hole placed in a uniform magnetic field that is asymptotically parallel to the axis of symmetry. The ansatz stipulates that the vector potential of the solution be in the plane spanned by the timelike Killing vector ξ µ t = (1, 0, 0, 0) and spacelike one ξ µ ϕ = (0, 0, 0, 1) of the stationary, axisymmetric black hole Since ξ µ t and ξ µ ϕ are pure geometric objects, they do not encode information on the applied magnetic filed B; such information is encoded in the coefficients (C t , C ϕ ). Here B is taken as a test field, so the metric of the stationary, axisymmetric black hole too does not encode any information on the applied magnetic field.
In this work, a spacetime metric has signature (+, −, −, −) and F µν = ∂ µ A ν − ∂ ν A µ . For neutral and charged black holes, Eq. (4.4) of Ref. [1] yields and respectively 1 , where B and Q are seen as perturbations, that is, if the metric of the background black hole is that of Kerr, then Qξ t µ /(2M) is, up to an additive constant, the one-form A µ dx µ = −Qr(dt − a sin 2 θdϕ)/ρ 2 of the Kerr-Newman black hole (with ρ 2 = r 2 + a 2 cos 2 θ). It is important to emphasize this point: The potential given by (1) and (3) is not an exact solution to the source-less 1 The sign "+" in (3) in front of Q is due to our metric-signature choice.

Maxwell equations
if the background metric is that of the charged black hole itself. Rather, it is a solution to (4) if the background metric is that of the corresponding uncharged black hole. For instance, in the Kerr background metric, the potential given by (1) and (3) is a solution to (4), but in the Kerr-Newman background metric the nonvanishing electric charge density J t and the ϕ current density J ϕ expand in powers of Q as which are zero to first order only. Even if rotation is suppressed (a = 0), J t is still nonzero: and its integral charge is also nonzero. Where does this electric charge density come from (the only existing electric charge is that of the black hole, which is confined inside the event horizon)? Because of the conservation of the total electric charge, the application of a uniform magnetic field does not generate current densities J µ outside the event horizon. Thus, as far as Q is considered as a perturbation, the potential given by (1) and (3) remains a good approximation for many astrophysical purposes. However, this fails to be the case if one is interested in the accretion phenomena that take place in the vicinity of the innermost stable circular orbit (ISCO) whose radius approaches that of the event horizon, for there the currents (J t , J ϕ ) cannot be neglected. One of the purposes of this paper is to provide an "exact" formula for the vector potential of a vacuum charged black hole immersed in a uniform magnetic filed parallel to its axis of symmetry. The purpose extends to include nonvacuum charged and uncharged black holes.

A. General considerations
The first thing we want to show in this section is that the expressions (2) are universal leading terms of more elaborate formulas for (C t , C ϕ ). The determination of these leading terms is purely geometrical and only depends on the asymptotic behavior of the metric of a (vacuum or nonvacuum), stationary, axisymmetric, asymptotically flat black hole. Here by asymptotical flatness we mean that, at spatial infinity, T µν → 0 ensuring Ricciflatness R µν → 0 and hence F µν ;ν → 0. The metric of such a black hole approaches that of a rotating star as the radial coordinate tends to infinity Using this, a direct integration of (4) yields the leading terms which were derived previously for the Kerr black hole (Eq. (6) of Ref. [2]). This behavior is universal and applies to all (vacuum or nonvacuum), stationary, axisymmetric, asymptotically flat black holes. The leading terms proportional to C t vanish as 1/r. Eqs. (6) are the expressions of the electromagnetic tensor F µν , expressed in spherical coordinates, of a uniform magnetic field, of strength B, parallel to the z axis provided we take There are three points to emphasize in the above derivation. First, notice that the derivation of (7) is valid whether the black hole is neutral or charged, for the presence of a charge, would certainly modify the expansion in (5), but would not modify the leading terms in (6). Second, we have made no assumption on the nature of the electrodynamics (linear or nonlinear) describing F µν . Thus, the value of C ϕ = B/2 applies equally to black holes with a linear electromagnetic source as well as to, generally regular (singularity-free), charged black holes with a nonlinear electromagnetic source provided they are asymptotically flat with a vanishing stress-energy at spatial infinity yielding R µν → 0 which, in turn, yields J µ → 0. Thirdly, we have made no link to general relativity nor to any of its modifications and extensions. Hence, no matter the theory of relativity describing the geometry and physics of spacetime and matter, the value of C ϕ = B/2 applies to all stationary, axisymmetric, asymptotically flat solutions (black holes, wormholes, etc ...) if they are placed in a uniform magnetic field B parallel to the symmetry axis. An instance of application of these arguments is the case of a non-Kerr black hole where it was shown that (7) applies [3].
As for C t , since we are dealing with asymptotically flat solutions, we are implicitly assuming that at spatial infinity linear electrodynamics (4) is sufficient for the description of the electromagnetic field. The charge of the black hole is given by the surface integral and this does not depend on the surface ∂Σ. This is the initial charge, if any, of the black hole. If ∂Σ is a surface of fixed radial coordinate r, then Σ µν = e µνθ ϕ |g|dθdϕ (where the totally antisymmetric symbol e µναβ is such that e trθ ϕ = +1). In the limit r → ∞, the metric of the black hole (5) is known yielding where we have used (7). In the limit r → ∞, we rederive the first equation in (2) [resp. in (3)] if the black hole is not charged (resp. charged). Corrections to (5) do not affect the leading term in (9).

B. Generalizing the ansatz
In more general situations where the Killing vector ξ µ is not a solution to the source-less Maxwell field equations, a linear combinations of all Killing vector with constant coefficient, which is also a Killing vector as in (1), may fail to be a solution to the sourceless Maxwell field equations (4). If the coefficients were taken as functions of the coordinates, it would be possible to determine them upon solving (4). This would allow to generalize Wald's formulas (2,3) and this is the purpose of this work. These generalizations are useful for a consistent analysis of (un)charged-particle dynamics around black holes.
It is worth mentioning that some specific generalizations of Wald's formulas (2, 3) to metrics not obeying the expansion (5) have been made but no general formulas were derived. For instance, it was shown [4] that (2) remains valid for the Kerr-Taub-NUT black hole but neither (2) nor (3) remains valid for black holes in Hořava-Lifshitz gravity and in braneworld where extensions to Wald's formulas (2, 3) have been performed in [5] and [6], respectively.
For the remaining part of this section, we rewrite the ansatz (1) in terms of (c t , c ϕ ) where C t = aB + c t and where the dependence on θ is not related to that on the rotation parameter a; that is, if rotation is suppressed (a = 0), the coefficients (c t , c ϕ ) may still depend on θ. Now, we need a metric formula to complete the task of integrating the source-less equations (4). There is no generic metric that describes all rotating black holes (and wormholes) nor a generic metric for static solutions. Nonetheless, the Gürses-Gürsey metric [18] where describes a variety of solutions including (A) Schwarzschild (a = 0), Reissner-Nordström (a = 0), Kerr, Kerr-Newman metrics, the Schwarzschild-MOG (a = 0) and Kerr-MOG black holes of the modified gravity (MOG) [7,8], and their trivial generalizations the Reissner-Nordström-MOG (a = 0) and Kerr-Newman-MOG, and some phantom Einstein-Maxwell-dilaton black holes [9,10]. It also includes (B) nonrotating regular black holes [11]- [17] and their rotating counterparts [19][20][21] as well as some nonrotating black holes of f(R) and f(T) gravities [22,23].
An important property of the metric (11), (12) is the following. In the most generic case, where f and (c t , c ϕ ) (10) are still arbitrary functions and for all values of the rotation parameter a , the first set of source-less equations (4) and two equations of the second set are satisfied: leaving only two differential equations to solve: Moreover, if rotation is suppressed (a = 0) then the coefficient (c t , c ϕ ) no longer depend on θ.
For the set (A) of singular solutions, the function f (r) is linear of the form where the constants ( f 1 , f 2 ) do not depend on the rotation parameter a. Table I gives the values of these constants for the set (A) of black hole solutions. The electromagnetic field of regular black holes does not satisfy the linear source-less Maxwell equations (4), except asymptotically, and their function f admits the expansion: f 1 r + f 2 + ∑ i=1 c i /r i . In this work we will not deal with metrics of regular black holes.

III. THE NONROTATING CASE: a = 0
We will first deal with the nonrotating case, where the black hole has the spherical symmetry, setting a = 0 the ansatz (10) reduces to:  (15) for the set (A) of singular nonrotating and rotating (vacuum and nonvacuum) black holes in terms of the mass M, the electric charge Q, and the universal ratio κ of the scalar charge Q s to the mass (of any particle and, particularly, of the MOG black hole). The Newtonian gravitational constant G N = 1. The black hole named "Reissner-Nordström-MOG" and "Kerr-Newman-MOG" were not derived in Ref. [7] since the author was not interested in astrophysical electrically charged solutions; these are trivial generalizations of the Schwarzschild-MOG and Kerr-MOG derived in Eqs. (11) and (35) of Ref. [7], respectively. NOMEN-CLATURE: "V" for "vacuum solution" (Ricci-flat: R µν = 0), "NV" for "nonvacuum solution" (non-Ricci-flat: R µν = 0), "N" for "neutral" (electrically uncharged), and "C" for "electrically charged".
Black hole for (c t , c ϕ ) do not depend on θ. With this last form of the ansatz, F tν ;ν = 0 implies where the prime denotes derivative with respect to r. This is solved by where (K 1 , K 2 ) are integration constants. Independently of the value of K 2 , the integral (8) performed on a sphere of radius r yields K 1 = −Q if the black hole is charged or K 1 = 0 if the black hole is neutral. We can take K 2 = 0, for this is an additive constant in the expression of the one-form A µ dx µ = −(Q/r)dt + K 2 dt+ terms proportional to dϕ. The result (18) was expected and constitutes a correction to Wald's term Q/(2M) in (3). Now, F ϕν ;ν = 0 reduces to the two equivalent differential equations where C ϕ = (B/2) + c ϕ , as defined earlier (10). In the generic case, where f is any function of r, a closed-form solution to (20) does not exist. For f linear (15), we ob-tain the solution where (L 1 , L 2 ) are integration constants. We immediately get rid of L 2 , for the physical electromagnetic tensor F µν would diverge at the horizons r ± = f 1 ± f 2 1 + 2 f 2 . At spatial infinity, C ϕ reduces to B/2 (7) yielding L 1 = B and For the linear case (15), we have obtained the solution for the vector potential For the Schwarzschild, normal or phantom Reissner-Nordström, Schwarzschild-MOG, and Reissner-Nordström-MOG black holes we obtain respectively Notice that it is the correction inside the parentheses, with respect to Wald's formula 3, that ensures the satisfaction of the source-less Maxwell equations (4). According to MOG theory [7,8], κ is the ratio of the scalar charge Q s of any particle to its mass m: κ ≡ Q s /m. This ratio is postulated to be universal and it is the same for all particles and massive bodies as black holes.
There is a couple of facts and conclusions to draw from (23).
1. We assume that the applied magnetic field is directed in the positive z axis (B > 0) and consider an uncharged black hole. On a charged particle, of electric charge q, the applied Lorentz magnetic force in the ϕ direction takes the form The third term is an extra term that would be missing had we applied Wald's formula (2). Recall that Wald's formulas (2) and (3) apply only to Schwarzschild and Kerr black holes. Notice that this extra force exists for neutral black holes of the MOG theory where f 2 = 0 (Table I). Since f 2 < 0, the extra force drives positively charged particles to accelerate (respectively, negatively charged particles to decelerate) in the increasing ϕ direction if they are approaching the axis of symmetry (dθ/ds < 0). This effect is new and cannot be neglected in the vicinity of the event horizon or the ISCO.
There is a similar extra magnetic term − 2qB f 2 sin θ cos θ in the applied Lorentz magnetic force in the θ direction.
These extra forces are attributable to the minimum coupling of the magnetic field, via the covariant derivative (4), to the stress-energy.
2. In fact, this effect, which was masked by the approximation made in (3), exists for all nonvacuum charged black holes no matter the nature of the stress-energy is. In our restriction (11), this effect has been derived in the linear case (15) but it should apply to neutral or charged generic configurations too where f expands as f 1 r + f 2 + power series in 1/r, as are the cases with the Ayón-Beato-García static black hole [12], the solutions derived in [17,24], and the black holes of the f (T) [23] and f (R) [25] gravities. (3) and (23) it was assumed that the magnetic field is a test field, thus neglecting its backreaction. The Einstein field equations are only approximately satisfied. If T µ ν (0) is the stress-energy corresponding to B = 0, then, for instance, in the nonrotating case all the new electromagnetic extra terms 2 , T µ ν EM , added to T µ ν (0) , when B is applied, are proportional to B 2 , if the background black hole is uncharged. The linear approximation (23) is valid if the end-behavior as r → ∞ of T µ ν EM is much smaller than that of T µ ν (0) yielding the constraints

In the derivation of both formulas
While the intergalactic magnetic field is supposed to be weak but these constraints show that (23) fails to provide a valid approximation as r → ∞. This very conclusion was stated in Ref. [26] concerning Wald's linear approximation (3): "... Wald's solution must break down as r → ∞, since the linearized solution is asymptotically flat, ...".
This conclusion extends also to charged solutions where, besides the constraint T t ϕ EM ∝ BQ/r ≪ 1, we have similar constraints to (30) 4. Ernst devised a procedure for generating Melvintype magnetic universes from Einstein-Maxwell solutions [27]. Being non-flat with a Melvin [28] asymptotic behavior, these solutions are useless for many astrophysical applications except in regions with strong magnetic field. The substitutes to Ernst's universes are the asymptotically flat solutions immersed in weak magnetic fields with linear vector potentials (3) and (23). All linear terms in powers of B have been determined in (23), however, we do not expect the next quadratic corrections to have simple mathematical structures even in the nonrotating case.

IV. THE ROTATING CASE: a = 0
This case is more involved. First of all, note that rotation mixes the electric field of the background black hole with the test magnetic field [30]. In this case it would not be possible to set the conditions constraining the test magnetic field, as we did in (30) and (31), unless the effects of rotation are weak, in which case the two constraints (30) and (31) remain valid to the linear approximation in the rotation parameter a. On the other hand, if the sources of the test magnetic field carry an electric charge density ρ, the dragging effects, which are proportional to a, cause the sources to accelerate in the geometry of the background black hole and thus enhance the magnetic field they generate 3 . Since in this work we are performing a general analysis in which we ignore the sources of the test magnetic field, it is safe to perform the analysis to the linear approximation in a to ensure that the generated B remains small, and the analysis remains valid, in all cases.
We seek a linear solution in a of the form (10) where we have written (c t , c ϕ ) and as sums of the new nonrotating contributions (23) and linear terms in a. It 3 It is well known that a charged rotating disk generates at its center a magnetic field linearly proportional to the disk's angular velocity. This extends to the case of an electromagnetic field around a Kerr black hole where the magnetic field, Eq. (29) of Ref. [30], is proportional to aq with q being the total charge of the electromagnetic source.
is understood from (32) that we restrict ourselves to the case f (r) = f 1 r + f 2 (15). The equations to solve (14) take the following forms in the a-linear approximation where, for instance, T (m,n) ≡ ∂ m+n T ∂r m ∂x n with x ≡ cos θ. Looking for solutions of the form where P 0 (x) = 1 and P 2 (x) = − 1 2 + 3x 2 2 are the Legendre polynomials, we were lead to the particular solutions 4 Finally, the expression of A µ takes the form To the linear approximation in a, the nonvanishing components of the electromagnetic tensor remain finite on the horizons , Using these expressions, it is straightforward to check that the r.h.s of (8) reduces to Q + O(a 2 ). The two invariants of the electromagnetic field, (where * F µν = 1 2 ǫ µναβ F αβ and ǫ µναβ = e µναβ |g| is the the totally antisymmetric tensor) too remain finite to the linear approximation in a.
As explained in the introduction section I, Wald's formulas (2) and (3) apply only to Schwarzschild and Kerr black holes where f (r) = f 1 r = Mr. However, had we tried to apply (3) to charged black holes with f (r) = Mr + f 2 (15) and metric of the form (11,12), we would obtain Now if we set f 1 = M in (38), we see that the terms proportional to M in both expressions (38) and (41)  Even if B is taken as a test field, this last term cannot be neglected near the axis of symmetry and in the vicinity of the horizon. As shown in (28), this produces the nonnegligible extra magnetic force on a charged particle This shows that the case of charged black holes lies beyond the realm of applicability of Wald's formula (3).
In the case of the Kerr-Newman black hole, besides what we just mentioned in the previous paragraph, Wald's formula 3 produces wrong and extra terms proportional to f 2 a = −Q 2 a/2 in F tr , F r ϕ , and F θ ϕ , and it yields a vanishing F tθ . Thus, the formulas (38) modify greatly the expressions of the forces acting on a charged particle and add extra terms to them. For instance, the force qF σ ϕ (dx σ /ds will still have the third extra term in (28) and the force qF σ θ (dx σ /ds will have two extra terms due to F tθ : Dropping the second term proportional to a 2 , there remains the term − qQ 2 Ba sin 2θ which, for all θ and g tt > 0, drives positively (resp. negatively) charged particles to (resp. away from) the axis of symmetry. Other more important applications will be given elsewhere [29].

Electromagnetic fields around rotating black holes
As mentioned in Sec. I, the electromagnetic field is taken as a test field in a given background metric; that is, the presence of the field does not modify the background metric. This approximation results in the two constraints (30) and (31).
To the best of our knowledge, there are no available exact, asymptotically flat, solutions to electromagnetic fields around (non)rotating black holes where the stressenergy of the electromagnetic field is taken into consideration except the well-known Kerr and Kerr-Newman solutions. Among existing approximate solutions, we find the stationary axisymmetric given in Refs. [30][31][32][33][34]. Exact solutions with a Melvin asymptotic behavior do, however, exist [26,27].
In Ref. [30], stationary axisymmetric electromagnetic fields surrounding a Kerr black hole were determined analytically. It was implicitly assumed that the stressenergy of the electromagnetic field is negligible not to affect the geometry of the background Kerr metric. Recall that the Kerr black hole has f (r) = Mr ( f 1 = M and f 2 = 0), so these field solutions do not extend to include the Kerr-Newman black hole nor more general black holes of the form (11) with f 2 = 0. In order to compare these field solutions with (37) we take Q = 0 in Eq. (29) of Ref. [30], which is the charge of the electromagnetic source (the background Kerr black hole is uncharged) and take Q = 0 in (37), which is the charge of the background black hole. Now, if we fix the remaining constants in Eq. (29) of Ref. [30] such that α i 1 = −B √ M 2 − a 2 /2 (≃ −BM/2 if a small), α r l = 0, α i l = 0 (l = 1), β r l = 0, and β i l = 0, this produces the leading terms proportional to r and r 2 in (37).

V. CONCLUSION
We have derived expressions for the vector potential and electromagnetic field of a rotating and nonrotating charged black hole immersed in a uniform magnetic field. The expressions are exact within the linear approximation and include all linear terms in both the rotation parameter and the magnetic field, thus introducing corrections to Wald's formulas. Since we have considered exact background black hole solutions, no special assumptions constraining the electric charge and the mass have been made. We have, however, made the implicit assumptions that the mass, electric charge, and rotation parameter are such that the solution is a black hole (not a naked solution without horizon(s)). The expressions apply to a variety of vacuum and nonvacuum solutions provided they satisfy asymptotically the linear Maxwell field equations.
As a first application we have observed the emergence of new extra force terms and evaluated some of them. Other applications will follow [29]. The generalization along with the corrections made to Wald's formulas are crucial for a consistent analysis of (un)charged-particle dynamics around black holes [4]- [6], [35]- [38], to mention but a few. From this point of view, some particledynamics analyses made in the literature have relied on Wald's formulas in cases where they do not apply.
The Kerr-Sen and Kerr-Newman-Taub-NUT black holes are charged solutions having dilaton fields resulting from dimensional reductions. Due to dilatons, the vacuum electromagnetic field equations (4) no longer describe the motion of F µν . Moreover, their metrics diverge from the form (11). It is another interesting topic to determine the vector potential of these black holes when they are immersed in a uniform magnetic field parallel to the axis of symmetry.