Magnetic response to applied electrostatic field in external magnetic field

We show, within QED and other possible nonlinear theories, that a static charge localized in a finite domain of space becomes a magnetic dipole, if it is placed in an external (constant and homogeneous) magnetic field in the vacuum. The magnetic moment is quadratic in the charge, depends on its size and is parallel to the external field, provided the charge distribution is at least cylindrically symmetric. This magneto-electric effect is a nonlinear response of the magnetized vacuum to an applied electrostatic field. Referring to a simple example of a spherically-symmetric applied field, the nonlinearly induced current and its magnetic field are found explicitly throughout the space, the pattern of lines of force is depicted, both inside and outside the charge, which resembles that of a standard solenoid of classical magnetostatics.

section_content: With the two recent papers [1, 2] we started a series of works aimed at studying quantum electrodynamics (QED) (as well as other nonlinear Abelian theories that may be historically traced back to [3] ) under the conditions where the intrinsic nonlinearity of the theory shows itself not only as interaction of electromagnetic fields with a strong background, but also with themselves. 

Manifestations of nonlinearity of the first type mentioned have been a focus of attention during many years since the pioneering work of [4, 5] (see [6] and more recent papers [7, 8] for some reviews of the subsequent advances in that field). The strong background was generated in the corresponding studies by the constant and homogeneous electroa e-mail: tg. adorno@gmail. com; tadorno@usp. br b e-mail: gitman@dfn. if. usp. br c e-mail: shabad@lpi. ru magnetic field (note, however, Ref. [9], where a certain inhomogeneity was introduced, and Ref. [10], where it was shown that a strong nonhomogeneous magnetic field is able to produce pairs of neutral fermions from the vacuum) and by the field of a plane wave (see the review of the laser-associated research in [11] ), because in these cases the influence of the background could be exactly taken into account for arbitrarily large value of their amplitude through the use of exact solutions of the Dirac equation available for such cases. The Dirac propagators for the virtual electrons and positrons in Feynman graphs for the vacuum polarization were the agents of interaction with the background field in the intermediate state. 

In those works the varying electromagnetic fields are treated as small perturbations of the background, and only the effects linear in their amplitudes are taken into account, such as birefringence, photon capture by a magnetic field [12] [13] [14] [15], modification of the Coulomb law [16] [17] [18] [19] [20], magnetic shift of the critical charge value [21, 22],-and even the positronium collapse [23, 24] may be placed among the effects of this class. The arena of applicability of these results is mostly the study of pulsars and magnetars, possessing sufficiently large magnetic fields. 

In contrast with the above, in Refs. [1, 2] and in the present paper we are considering the effects quadratic and cubic in the amplitude of the perturbation. As a matter of fact, two important special examples of such effects were studied before, which were the processes of photon splitting in a magnetic [5, [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] and crossed [37, 38] fields, and the light-by-light scattering [39], all taken on the photon mass shell. Our goal is to deal with excitations of the vacuum, different from photons, subject to a nonlinear version of the Maxwell equations, stemming from QED or, more generally, intrinsic to any other nonlinear electrodynamics. Whereas handling many-photon matrix elements beyond the photon mass shell, necessary for addressing a general nonlinear problem, turns out to be overcomplicated, we succeeded to indicate a simple approximation describing nonlinear effects in a universal manner and independent of an expansion in powers of the fine-structure constant α, as far as QED is concerned. All kernels in the nonlinear integro-differential Maxwell equations are given in terms of the variational field-derivatives of the effective action, defined as [40] the generating functional of irreducible vertices in QED, or as an action that fixes a theory in other versions of nonlinear electrodynamics. The approximation we are dealing with is referred to as the local or infrared approximation. It assumes that the effective action functional is local, i. e., it does not depend on the space-time derivatives of the field strength. True, this assumption restricts the range of applicability to only slowly varying fields in space and time, but it enables us to efficiently advance in describing the effects cubic and quadratic in the field strength. By acting along these lines where there is no background field [2], we reproduced the known [39, 41, 42] correction to the Coulomb field which is cubic in the charge that produces it, and found cubic equations for dipole moments of selfinteracting fields of magnetic and electric dipoles that may also be viewed upon as nonlinear renormalizations of these quantities. In Ref. [1] we studied a quadratic response of the background constant and homogeneous magnetic field to an applied electric field of a static charge at rest, and we found this response to be purely magnetic. In the present paper we continue the investigation of that magneto-electric effect, and we establish that the static charge placed into a background magnetic field is a magnetic dipole with its magnetic moment proportional to the charge squared and parallel to the background field, unless the charge distribution violates the initial cylindric symmetry of the problem. (This situation resembles our results [43] [44] [45] in noncommutative electrodynamics. ) More explicit formulas for the magnetic field in the short and long ranges, for its lines of force and for the magnetic moment, are presented referring to a simple example, where the applied electrostatic field is central-symmetric and would correspond, if used in a nondispersive vacuum, to a charge distributed homogeneously inside a finite-radius sphere. 

In the rest of the present Introduction we recall the basic equations of nonlinear electrodynamics, truncated at the third power of the varying field, against a constant homogeneous background, define the kernels in them as the second-and the third-rank polarization tensors in terms of the field derivatives of the effective action, and introduce nonlinearly induced currents. We give expressions for the varying fields in terms of the currents based on the use of an eigenvector expansion of the second-rank polarization tensor and the photon Green function. In Sect. 2 the second-and third-rank polarization tensors in a magnetic field are written in the infrared approximation as expressed in terms of the second and third derivatives of the delta functions of the coordinate differences. The previously found linear-response corrections to the Coulomb field of an arbitrarily distributed static charge at large distances in a magnetic field are given for completeness, and the general structure of the quadratic response, which is purely magnetic, is analyzed in terms of the eigenmode contribution. In Sect. 3 very explicit expressions for the induced current and the magnetic field are fully elaborated using the simplest example, where the applied electric field is that of a homogeneously charged sphere of finite radius. Differential equations for the shape of the magnetic lines of force are solved, and the resulting pattern is drawn in Figs. 1 and 2 inside and outside the charge. The magnetic moment of the charge, proportional to its square, is given in terms of the background magnetic field and the radius of the charge. In Sect. 4 it is shown that spherically nonsymmetric charge distributions also are characterized by a magnetic dipole moment, which is parallel to the background magnetic field, if the charge distribution does not specify any new direction in the space.
section_name: Introduction

section_content: The exact electromagnetic field equations of QED with an external 4-current j μ are the Euler-Lagrange equations

where S Max [A] is the free Maxwell action, [A] is the effective action, and L (x) is the effective Lagrangian. By the effective action we understand in QED the generating functional of one-particle-irreducible vertices [40]. Alternatively, it may be any action defining a nonlinear electrodynamics other than QED. Due to the gauge invariance it, as a matter of fact, depends only on the field strengths F μν, and not on the 4-vector potentials A ν. Besides, only the relativistic invariant combinations

where the dual field tensor is defined as Fρσ = 1 2 ρσ λκ F λκ, of the field strengths make up the arguments of and L. In QED the effective action contains exhaustively the final information of the theory in the photon sector, and it is sub-ject to a calculation within one or another dynamic scheme or approximation, especially the perturbation theory. 

Expanding (1) in a power series of the small electromagnetic field a μ (x) = A μ (x) -A μ (x) over the external background A(x) of a constant and homogeneous magnetic field B = (∇ × A), and restricting ourselves to the nextto-leading term, the minimum action condition becomes the nonlinear Maxwell equation:

where ρα x, x, αβρ x, x, x are the second-and third-rank polarization tensor in an external magnetic field, respectively. Note that as long as the external field is constant in space-time, the polarization tensors (3) are functions on the differences of their arguments. In obtaining (2) the zeroorder power of the field a μ (x) does not appear, since the space-time-independent external field

= 0. The power series has been truncated to the next-to-leading correction (i. e., we have neglected ∼ O a 3 ). 

Defining the nonlinear current as

we may write (2) in the following way:

While solving the nonlinear set of Maxwell equations ( 5 ), (4) we should not, strictly speaking, exceed the initial accuracy. This implies that we treat the nonlinearity iteratively. To this end we divide its solution into two parts as a ν (x) = a lin ν (x) + a nl ν (x), with a lin ν (x) a nl ν (x). Then, defining the linear field a lin ν (x) as a solution to the equation

we get, as the first iteration, the result that the nonlinear correction a nl ν (x) to it is subject to the linear inhomogeneous equation

where

is the nonlinear current (4) taken at the linear value (9) of the field a ν (x) = a ν lin (x). The solution of ( 6 ) and ( 7 ) may be written as

The photon Green function in D μν (x, x ) above is defined as the inverse operator:

Its Fourier transform D ν ρ (k) with respect to the coordinate difference should satisfy the following algebraic inhomogeneous equation:

To solve this equation it is convenient to use the diagonal representation for the second-rank polarization operator in a magnetic field [46] [47] [48] [49] [50],

in terms of the mutually orthogonal 4-vectors (2) μ = ( Fk) μ, (3) μ = (Fk) μ, (4) μ = k μ, (14) where

which are the eigenvectors of the polarization operator μτ (c)τ = c (k) (c) μ, (15) the scalar functions c (k) being for its four eigenvalues, 4 (k) = 0. The eigenvalues c (k) depend on F and on any two of the three momentum-containing Lorentz invariants

, where F is taken on the external field, 2F = B 2. The solution of ( 12 ) has an arbitrary longitudinal part:

It also has a diagonal form in the same terms as (13). This propagator has three components, corresponding to separate eigenmodes. Each of them has a pole in the 4-momentum plane, where solutions of the corresponding dispersion equations lie, i. e. on the photon mass shell, defined by the equation

Now, the solutions ( 9 ), ( 10 ) may be, respectively, written as

2 Response of magnetized vacuum to static electric field in the infrared approximation

In the rest of the paper we shall be treating the equations of Sect. 2. 1 in the low-momentum-low-frequency (infrared) approximation, k μ ∼ 0, which stems from the assumption that the effective action [A] is a local functional of the field strengths F μν in the sense that it does not contain their space and time derivatives. Examples of such an action are the Heisenberg-Euler action available in the one-loop [39] and two-loop [51] approximations in QED, the Born-Infeld [52] action, etc. Within the local limit the second-and thirdrank polarization tensors (3) were calculated in [1] to give the result 2

2 The fourth-rank tensor in the same approximation is also available [2]. 

where

which expresses them in terms of the derivatives of the effective Lagrangian taken at the constant external field value, F μν = const.,

It is taken into account that once the external field is purely magnetic in a certain Lorentz frame, the invariant G for it is zero while F is positive. 

It was established in [53] that the second-rank polarization operator (19) has indeed the structure of ( 13 ) with the eigenvalues a (k) in the infrared regime being

Henceforth, we shall be dealing only with sources that are static in a reference frame, where the external field is magnetic, i. e. time-independent charges at rest in a magnetic field,

where the tilde marks the Fourier-transformed function.
section_name: Nonlinear Maxwell equations

section_content: Employing the source ( 27 ) in ( 17 ) and taking into account that at k 0 = 0 out of all the three (nontrivial) eigenvectors ( 14 ) only

μ has its zeroth component different from zero, while its spatial components disappear

(We disregarded the longitudinal part k μ, which does not contribute to the field strength. ) So, naturally, only a static electric field is produced by a static source at the linear level (9), and no magnetic field. Equation ( 28 ) is approximation independent. In the infrared limit, ( 25 ) is to be used with k 0 = 0 and

in the special frame, when the latter is so oriented that axis 3 coincides with the direction of the external magnetic field B, B = B 3 = F 12, B 1,2 = 0. Then in the coordinate space the linear potential (28) becomes

where we have used expressions for longitudinal and transverse dielectric constants from [53] (see also [54] )

where the field invariant F is henceforth taken on the external field, 2F = B 2, unlike its previous more general definition in (1). The subscript ⊥ indicates projection onto the (1, 2)-plane. The behavior of (29

where Q = q(0) = q(x)d 3 x is the total charge with the understanding that the charge density is either compactly distributed inside a volume or decreases sufficiently fast outside it, so that this integral converges. By making the change of variables ε

where the notations 1/2 were used. This anisotropic Coulomb law, resulting from the infrared limit (19) of the second-rank polarization operator, is thereby the long-distance asymptotic behavior of the electrostatic potential produced by a charge, locally distributed in space, in a constant magnetic field. 

It may also be useful to write the scalar potential (31) in an O(3)-invariant way as a function of two rotational scalars:

Previously [16] [17] [18] [19], that potential was found in QED in the whole space, the vicinity of the charge-where the potential has the Debye form-included, starting from the off-shell calculations of the full (free of the restrictive assumption k μ ∼ 0) second-rank off-shell polarization operator in a magnetic field. This was first performed within the accuracy of one fermion loop in [46] [47] [48] [49] [50]. If the one-loop Heisenberg-Euler effective Lagrangian is taken for L in (31), it makes (a corrected 3 form of) the large-distance behavior of the potential in QED.
section_name: Linear response-modified Coulomb law at large distances

section_content: When the four-vector potential in expression ( 4 ) is chosen as a μ (x) a μ lin (x) = δ 0μ a 0 lin (x), so as to carry only electrostatic field E (x) = -∇a 0 (x), the nonlinear current (4), approximated by (8) in accord with the iteration (7) (henceforth we omit the explicit indication that it is taken on linear fields), becomes [1], up to the third and fourth powers of the applied field,

after the expression (20) for the third-rank polarization tensor in the local limit has been used. Here every derivative acts on everything to the right of it. Note that in (33), E is the applied electric field, the nonlinear response to which is under consideration. Correspondingly, ( 33) is quadratic with respect to E. To the contrary, the external magnetic field B = (∇ × A), B = (2F) 1/2 enters (33) with all powers, since the coefficients in it depend on B in a complicated way according to their definitions (22) and (23). 

Let us discuss the structure of the nonlinear correction to the electromagnetic field caused by the current (33) following (10). We appeal to the representation (16) for the photon propagator in it. First we note that the mode c = 2 does not contribute, since j 0 nl = 0, and

(2) i = ( Fk) i (14) disappears, when multiplied by the Fourier transform j nl (k) δ(k 0 ) of (33). Also the zeroth components of the other two eigenvectors (14) (1,3) 0 vanish if taken at k 0 = 0. Thus we are left with (34) which indicates that the quadratic response to a static electric field is purely magnetic. 

It can be shown that (34) are in fact exact relations in the electrostatic case, independent of the infrared approximation. This implies that in that instance only modes 1 and 3 propagate the magnetic field. However, in a spherically symmetric infrared example as to be considered below we have j jnl (k) (1) j = 0, so only the term c = 3 contributes in the nonlinear electromagnetic field (34). 

If the linear vacuum polarization (in fact, the magnetization) is neglected, 4 1,3 = 0, the nonlinear magnetic field h = k×a nl obtained from (34) satisfies the standard Maxwell equation k 2 a nl = j nl. In this case this field follows from (33) [1] :

where

3 Quadratic magnetic response to a spherically symmetric applied electric field: simple example

In this section, in order to present the magneto-electric effect in its most explicit way, we shall consider the magnetic field, which is the response of the vacuum to the applied electric field, whose vector potential is chosen to be the following smooth spheric-symmetrical Coulomb-like function:

Here θ (z) is the Heaviside unit step function, defined as

and δ (z) stands for the Dirac delta function. Equation (37) supplies us with the simplest example, where the magnetic field comprising the nonlinear vacuum response can be explicitly studied, the shape of the lines of force being fully described. If not for the linear polarization, the potential (37) would be the field of the extended spherically symmetric charge distributed with the constant density ρ (r ) inside a sphere r ≤ R with the radius R:

However, it should be kept in mind that with the account of the linear vacuum polarization, the potential distribution (37) cannot be supported by any spherically symmetric charge, strictly localized in a finite space domain. To find the genuine source of the field (37), one should apply the operator in the left-hand side of (6) to it. The result looks like

This charge density is cylindrically symmetric and extends beyond the sphere, r > R, decreasing as 1/r 3, or 1/x 3 3 far away, depending on the direction. 

The particular result to be formulated in the present section that the magnetic response to an electrostatic field implies that the charge giving rise to the latter carries a magnetic dipole will be confirmed for a general charge density distribution in the next section. Only the expression for the magnetic moment will be less explicit, nor the expressions for the induced magnetic field and the shape of its lines of force in the region closer to the charge. The aim of the present section is just to detail this, appealing to the simplified example of (37). Now we proceed with the applied potential (37). It yields the following electric field:

In writing the expression for (40) we took into account the continuity of (37) and of its derivatives at r = R. As a result, Dirac delta terms stemming from the differentiation of the step function could be omitted (see the appendix for a general discussion of the subject). This simplification will be used every time functions are continuous.
section_name: Quadratic response-magneto-electric effect

section_content: Taking into account the electric field ( 39 ) and

the contributions (33) to the nonlinear current take the form

Therefore

Again, thanks to the continuity of E (x) ( 39), the expressions (42), (43) do not contain the δ-like contributions that might come from the differentiation of the Heaviside θ -function. 

But the nonlinear current is discontinuous at the surface of the sphere r = R. We shall see below that the magnetic field produced by it is still continuous. The lines of the nonlinear current ( 43 ) are circular and lie in planes orthogonal to the external magnetic field. Its density decreases as the sixth power of the distance from the center of the charge, or of the distance from the axis, parallel to the external magnetic field-for large distances. 

Equation ( 43 ) is our final expression for the current that may be of independent interest. In calculating the magnetic field produced by it in the next subsection, however, we shall not exploit it, but refer to expressions of the previous sections. 

Observe that the proportionality of the current j nl (x) above to (x × B), and hence its circular character in the transverse plane, is valid also for any spherically symmetric field distribution like (39), irrespective of the special form (40). This property means that an expansion of (33) over the (spatial part of) the eigenvectors of the polarization operator in a magnetic field (14) does not contain, in the momentum space, a nonvanishing contribution proportional to the vector b (1), whose components are

(1)

), but it is proportional to the vector b (3), such that

which is thereby the only vector contributing in the expansion of the nonlinear current in the spherically symmetric case. The three vectors b (1), b (3), and b (4) = k are mutually orthogonal. There is no contribution proportional to b (4) due to the continuity ∇j nl = 0. As for the eigenvector (2) μ, its spatial part is zero in our static case of k 0 = 0. Correspondingly, only the value c = 3 appears in the expansion of the nonlinear magnetic field (34). We do not know if this circumstance may be a general consequence of the spherical symmetry, independent of the infrared approximation.
section_name: Nonlinearly induced current

section_content: We shall calculate here the magnetic field, nonlinearly induced by the electrostatic field (37) within the infrared approximation, basing ourselves on (35), which assumes the neglect of the linear vacuum polarization. (35) In order to find the integrals
section_name: Nonlinearly induced magnetic field

section_content: it should be noted that the three parts have the following structure:

where

and u(r ) and w(r ) are the scalar coefficients in the tensor decomposition

Eur. Phys. C (2014) 74:2838

The basic angular integrals involved in ( 45 ), ( 46 ) are (we refer to (2. 222) in [56] for their values)

and

One can immediately verify that v (r ) (45) takes the form

Further, once (46) does not depend on B, one can choose B ⊥ x to find

When getting this result we counted off the angle ϑ from the radius vector x, and ϕ from B, so that B • y = By sin ϑ cos ϕ, and we took into account the relation 2π 0 cos 2 ϕdϕ = π. Analogously, by choosing B x we find from (46) that

since now B • y =By cos ϑ. Using (50) for w (r ) we get the function u (r ) in (46),

Now ( 49 ), (50), and ( 52 ) define the contributions (44) to the field h (x) according to (35). These equations are valid, generally, for an arbitrary spherically symmetric field distribution of the form (39), provided that E (r ) there decreases sufficiently fast at large r to guarantee the convergence of the remaining integrals over y. For E 2 (y) taken as (41) the remaining y-integrations in ( 49 ) and ( 51 ) can be explicitly done with the help of ( 47 ), (48), and their calculation is illustrated in the appendix. 

We obtain in this way

for (49), and

for ( 51 ), (52), and (50). With these results the integrals ( 44) have the form

where

Using these representations it is straightforward to make sure that the functions ( 55 )-( 57 ) are continuous in the point r = R. So are also all their first and second derivatives with respect to the coordinate components. Consequently, the Dirac delta functions and their derivatives, which stem from differentiation of the step functions in the calculation of ∂ i ∂ k I k (x) /4π from (44), do not contribute (see the appendix). As a result these derivatives can be left in their simplest form,

The final form of ( 59 ) is

(60)
section_name: Second contribution in

section_content: With the explicit form (60) of ∂ i ∂ k I k (x) /4π we proceed to evaluate the total magnetic field (35). Bearing in mind ( 39 ) and ( 41 ) the field h i (x) in ( 36 ) is written as

and the total magnetic field ( 35 ) has the final form

Here h in i (x) represents the total magnetic field for points inside the sphere (r < R), while the designation h out i (x) is reserved to the field outside the sphere (r > R). The total magnetic field ( 62) is continuous at r = R. The orientation of the magnetic field (62) will be revealed in the next subsection, where we present the shape of the lines of the magnetic field. 

The long-range contribution of (65), h LR i (x), behaves like a magnetic field generated by a magnetic dipole:

with µ being the equivalent magnetic dipole moment, given by
section_name: Total nonlinear magnetic field

section_content: In this section we are going to explore the dependence of the magnetic moment (67) with respect to the external applied magnetic field B. To do that one has to consider the corresponding dependence of the coefficients L FF, L GG, and L FGG on B. Such coefficients essentially depend on the model under consideration, but, confining ourselves to QED and working within the local limit approximation, they have the specific form provided by the Heisenberg-Euler effective Lagrangian [58] [59] [60] [61]. They have been considered before in [53] and due to this fact we use here the expressions previously derived there to obtain 5

where B Sch = m 2 b = Sch. Using the latter result, the magnetic moment (67) has the form

The negativity of (68) and of μ (69) indicates that the magnetic moment is directed oppositely to the background magnetic field. 

This integral does not have an analytical solution. To show the dependence of the magnetic moment with respect to the external field we plot the numerical results of the ratio -μ/λ for each given value of b within the range 10 -2 ≤ b ≤ 50 (Fig. 1 ). Although (69) does not have an analytical solution, one can estimate its asymptotic behaviors for small and large values of the external magnetic field B. In the first case, for t sufficiently small, the exponent e -t/b is approximately zero. 

Separating the integral above in two parts, where the limit is chosen such as b, one can write

since the second integral in the first line is practically zero (t is always t b). Once the function f (t) has a maximum at t = 0, one can conclude that the greatest contribution for (70) is

by which, after two integrations by parts, we obtain

Finally the asymptotic form of (69) is linear in b,

In the large-field asymptotic regime B B Sch the integral (68) decreases as -α 3π e m 2 B providing in turn a constant value to the magnetic moment (see the horizontal dot-dashed line in Fig. 1 ) μ = -λ/3.
section_name: Magnetic moment of a spherical charge

section_name: Magnetic lines of force

section_content: Here we are going to establish the form of the lines of force of (63 ), first on the inside of the sphere. To this end, in the same way as in Sect. 2. 1, we direct the axis x 3 along the magnetic field B and represent the vector (63) in the orthogonal basis of unit vectors e (1)i, e (3)i directed along the orthogonal axes 1 and 3 as with

being functions of B, independent of the coordinates x. We set x 2 = 0, since the full pattern of the lines of force is to be obtained from the one in the plane (3,1) by rotating along axis 3 due to the cylindric symmetry of the problem, so r 2 = x 2 1 + x 2 3. Equalizing the derivative dx 3 dx 1 with the ratio

provides us with the differential equation for the shape of the line of force x 3 (x 1 ). With the new notations y = x 3 /R and z = x 1 /R this differential equation follows from (74):

where

This is the so-called second-type Abel first-order differential equation with the family of solutions [57]

parametrized by the integration constant z 0. Extreme points of the lines of force given by (78) are

The corresponding extremum value y extr of the vertical coordinate turns to zero for the curve corresponding to the largest admitted value of the integration constant, (z 2 0

. This closed curve degenerates to a point. Its position at the abscissa axis is z = z foc = ( -β γ ) 1/2. We call this point the focus of the lines-of-force pattern. Larger values of the integration constant would not give rise to any line of force, since they would make y complex for any z. Therefore the integration constant can be taken within the range z foc 0 > z 0 > 0. As we let the parameter z 0 diminish down to the zero value, we pass to lines of force that go farther and farther from the horizontal axis. In the limit z 0 = 0 we reach the ultimate curve that passes through the origin z = y = 0 and coincides with the y-axis. The focal point may lie both inside and outside the sphere, depending on whether -β γ = -A D is smaller or larger than unity. Bearing in mind the asymptotic behavior at large magnetic field, B m 2 /e,

the basic quantities forming the coefficients A, D, and g (75 ) in QED (see e. g. [53] ), we find for z foc the value (7/6) 1/2 > 1 outside the sphere in this limit. 

In the limit of a pure point-like dipole β → 0, the focal point tends to the origin, and all the lines of force are squeezed between these two points. In the large-field regime (79), the constants are. 

For B → ∞, the coefficients above become

and the magnetic curves take the final form

The family of magnetic lines labeled by positive values of the integration constant z 0 in the interval 0 < z 0 < 7 √ 66 are drawn following (81) in Fig. 2. For negative z 2 0, the corresponding curves lie completely outside the sphere z 2 + y 2 = 1, rounding from outside the family presented in this figure. We are not interested in showing them, because our starting equations in this subsection belong to the interior of that sphere. For z 2 0 > 49 66 the solutions (81) are no longer real. The values taken for parametrizing the six curves in Fig. 2 are indicated in the drawing. We must mistrust those parts of the curves in Fig. 2 that belong to the exterior of the sphere, and our next task is to obtain the continuation of the magnetic lines of force to that region.
section_name: Interior

section_content: Referring to the same basis and reference frame as in the previous subsection, the magnetic field outside the sphere (65) reads

where y = x 3 /R and z = x 1 /R. The ratio h out 3 (x) / h out 1 (x) can be expressed as

Taking into account the asymptotic behavior at large magnetic field B m 2 /e (79), the coefficients A, B, D, and

then β, γ, and E read

Equating the derivative dy/dz with the ratio h out (86)
section_name: Exterior

section_content: Fig. 3 Magnetic dipole lines of a static charge in external magnetic field exampled with B = ∞. b The pattern to be trusted outside the charge, y 2 + z 2 > 1, following the solution (86). For each choice of z 0 we extract from (81) the corresponding value for y(z) at the border of the sphere. See (87) for some boundary conditions to (86). c United pattern to be trusted throughout This equation does not have analytic (closed) solutions. We found them by using numerical methods. The integration constant is fixed by the matching requirements with the pattern in Fig. 2 : we demand that solutions of (77), for each fixed z 0, have the same numerical values as (86) at the border of the sphere y 2 + z 2 = r 2 /R 2 = 1. In this way the continuous continuation of solutions of (77) to the outer region, where (86) actually holds, is achieved. Figure 3c shows the overall pattern of magnetic lines to be trusted everywhere, wherefrom the lines beyond their domains of definition have been deleted. The values of y (z) and z at the border of the sphere are listed below for each integration constant z 0 : z 0 = 0. 1 → z 0. 186, y (0. 186) 0. 983, z 0 = 0. 2 → z 0. 349, y (0. 349) 0. 973, z 0 = 0. 3 → z 0. 486, y (0. 486) 0. 874, z 0 = 0. 4 → z 0. 604, y (0. 604) 0. 797, z 0 = 0. 5 → z 0. 707, y (0. 707) 0. 707 z 0 = 0. 6 → z 0. 8, y (0. 8) 0. 6. 

(87)

The magnetic lines in 3-c) remind one very much of the standard pattern of those of a finite-thickness solenoid in classical magnetostatics.
section_name: (c) (b)

section_content: Here we search for an extension of (66) to a spherically nonsymmetric applied electric field. Such a generalization provides a more general form of the magnetic dipole moment µ. 

Let us first see how the result (66), (67) can be directly reproduced by considering the long-range behavior of the magnetic response (35) to the spherically symmetric electric field (39), (40). According to (36) the field h i (x) in the far-off domain reads

where we have restricted ourselves to the leading contribution at large r = |x|. The leading behavior of the quantities I k (x) in ( 35 ) is

provided that the integrals here converge. Then

The field h i (x) (88) falls off faster than this, namely as 1/r 4, hence its contribution into the first line of (35) can be neglected as compared to (89). So, the large-distance behavior of the nonlinear magnetic field h i (x) (35) is just (89), i. e., that of a magnetic dipole. Its magnetic dipole moment μ LD i is

The result (90) agrees with the previous result (67) in case the spherically symmetric E is specialized to (37). To ensure this, it suffices to substitute expression (61) for h i (y) into (90) and fulfill the integration, which converges both at the lower, y = 0, and the upper, y = ∞, limit. 

However, the validity of the result (90) is much wider. For instance, let us take (31) or, equivalently, (32) for the scalar potential responsible for the remote cylindrically symmetric electric field of a static extended charge Q = Ze, whose density decreases sufficiently fast at infinity, but which is otherwise arbitrary, not subject to any symmetry. Recall that this cylindrical, instead of spherical, symmetry became in Sect. 2. 1 the effect of the linear vacuum polarization in an external magnetic field. It is easy to assure oneself that when this electric field is substituted into (36) for h i (x), the resulting expression in place of (88) also decreases as 1/r 4, in the same way. Hence, we are left again with (89) for the largedistance asymptote of the nonlinearly induced magnetic field h i (x) (35). Since the electric field ( 32 ) is invariant under rotations around the external magnetic field B, the latter remains the only special direction in the space. Consequently, the magnetic moment (90) is directed along B, in the same way as for (67).
section_name: Beyond the spherical symmetry of the applied field

section_content: In this work it was shown that a static charge, apart from being a source of the customary Coulomb-like electric field, is also a nonlinear source of a magnetic field. This field is generated due to a nonlinearly induced current caused by a constant and homogeneous external magnetic field. As a result, the long-range magnetic field behaves like a magnetic field generated by a magnetic dipole moment. In other words, the extended charge has a long-range magnetic dipole character. The magnetic field lines resemble the well-known magnetic dipole structure. 

The validity of equations found here for the nonlinear magnetic response of the magnetic background to an applied electric field is restricted to fields smooth in time and space. They can be directly applied to charged large astrophysical objects, but they lead to overestimation, where small objects as charged mesons and baryons are concerned. Therefore, to make such an application reasonable, one needs to go beyond the infrared approximation. To this end QED calculations of three-photon diagrams in an external magnetic field must be efficiently exploited beyond the photon mass shell. We hope to come back to this more complicated problem in future work. 

Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited. Funded by SCOAP 3 / License Version CC BY 4. 0.
section_name: Conclusion

section_content: Performing the derivations of the potential (37) leads to Dirac delta functions and their derivatives as well. Taking into account the smoothness conditions at r = R, one is able to simplify the explicit form of some quantities under consideration. One can see that where f (r ) represents any function well-behaved at r = R. The same idea can be generalized to any function which depends on d 2 a 0 (r ) /dr 2 or higher derivatives. 

In order to evaluate the integrals I k (x), in the expression (35), for the total magnetic field, one has to evaluate v (r ), u (r ) and w (r ) (49)- (52). All of these functions can be conveniently written as sums of two other integrals. For example, we write (49) as v (r ) = (2π/r ) [v 1 (r ) + v 2 (r )] where

Similarly

Then (49) takes the final form (53). 

Besides, it should be noted that u (r ) and w (r ) can be written in a simplified form Substituting these results in (96) and using (95), the scalar functions take their final form, (54).
section_name: Appendix

section_content: Acknowledgments T. C. Adorno acknowledges the financial support of FAPESP under the process 2013/00840-9 and 2013/16592-4. He is also thankful to the Department of Physics of the University of Florida for the kind hospitality. D. Gitman thanks CNPq and FAPESP for permanent support; in addition, his work is done partially under the project 2. 3684. 2011 of Tomsk State University. A. Shabad acknowledges the support of FAPESP, Processo 2011/51867-9, and of RFBR under the Project 11-02-00685 -a. He also thanks USP for kind hospitality extended to him during his stay in Sao Paulo, Brazil, where this work was partially fulfilled. The authors are thankful to C. Costa for discussions.