Interaction between two point-like charges in nonlinear electrostatics

We consider two point-like charges in electrostatic interaction within the framework of a nonlinear model, associated with QED, that provides finiteness of their field energy. We find the common field of the two charges in a dipole-like approximation, where the separation between them R is much smaller than the observation distance r : with the linear accuracy with respect to the ratio R / r, and in the opposite approximation, where R≫r,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R\gg r,$$\end{document} up to the term quadratic in the ratio r / R. The consideration proposes the law a+bR1/3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a+b R^{1/3}$$\end{document} for the energy, when the charges are close to one another, R→0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R\rightarrow 0$$\end{document}. This leads to the singularity of the force between them to be R-2/3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R^{-2/3}$$\end{document}, which is weaker than the Coulomb law, R-2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R^{-2}$$\end{document}.


Introduction
Recently a class of nonlinear electrodynamic models was studied [1] wherein the electrostatic field of a point charge is, as usual, infinite in the point where the charge is located, but this singularity is weaker than that of the Coulomb field, so that the space integral for the energy stored in the field converges (and also the scalar potential is finite in the position of the charge [2]). This class unites Lagrangians [2][3][4] that grow with the field invariant F = B 2 − E 2 /2 faster than (−F) w , w > 3/2 (see Ref. [2] for a more subtle estimate of the boundary of the necessary growth). Thus, the simple quadratic effective Lagrangian first considered as regards a different aspect in [5] is also included into the class under consideration.
The infiniteness of the field near the charge distinguishes the class under consideration from many other models (see Refs. [6][7][8][9][10][11][12] and the references therein) with finite selfa e-mail: breev@mail.tsu.ru b e-mail: shabad@lpi.ru energy of the point charge, allied to their famous prototype, the Born-Infeld model [13], where the finiteness of the field is achieved at the cost of square-root nonanalyticity of the Lagrangian, which leads to an infinity to the Maxwell equation. The most popular application [6][7][8][9][10][11][12][14][15][16][17][18][19][20][21][22] of these models is to combine them with General Relativity in order to study their effect on the initial singularity and on the evolution of the Universe. In contrast to the Born-Infeld model, the models from the class of Ref. [1] refer to nonsingular Lagrangians that follow for instance from the Euler-Heisenberg (E-H) effective Lagrangian [23,24] of QED truncated at any finite power of its Taylor expansion in the field. This allows us to identify the self-coupling constant of the electromagnetic field with a definite combination of the electron mass and charge and to propose that such models may be used to extend QED to extreme distances, smaller than those for which QED may be thought of as a perfectly adequate theory.
More advanced approaches based on the Euler-Heisenberg Lagrangian that do not depend upon any assumption of smallness of its field argument (the background field) and do not hence appeal to expansion of the Lagrangian in powers of the background fields, have received attention, as well, under the restriction, however, that fields not-too-fast-varying in space and time are studied as solutions of the nonlinear Maxwell equations. Among the nonlinear effects studied, there are the linear and quadratic electric and magnetic responses of the vacuum with a strong constant field in it to an applied electric field [25], with the emphasis on the magneto-electric effect [26][27][28] and magnetic monopole formation [29]. Also selfinteraction of electric and magnetic dipoles was considered with the indication that the electric and magnetic moments of elementary particles are subjected to a certain electromagnetic renormalization [30] after being calculated following a strong-interaction theory, say, QCD or lattice simulations. Interaction of two laser beams against the background of a slow electromagnetic wave was studied along these lines, too [31]. The finiteness of the field energy allows one to develop a soliton view on a moving point charge [3,32].
In the present paper we are extending the consideration to cover the electrostatic problem of a system of two point charges that interact following nonlinear Maxwell equations stemming from the Lagrangian quadratic in the field invariant F. Their common field is not, of course, just a linear combination of the individual fields of the two charges. The nonlinear problem is outlined in Sect. 2, where we present the nonlinear Maxwell equations and give them a certain form in Sect. 2.2 apt for finding the approximate solutions of Sect. 3. Once the field energy is finite it is possible in principle to define the attraction or repulsion force between the charges as the derivative of the field energy with respect to the distance R between them. Contrary to standard linear electrodynamics, this is evidently not the same as the product of one charge by the field strength produced by the other! This rule holds true only if one of the charges has a much smaller value than the other.
In Sect. 3 we develop the procedure of finding the solution to the static two-body problem in two opposite approximations determined by the ratio of the distance R, to the coordinate of the observation point r 1 . Where this ratio is small, R/r 1, we find in Sect. 3.1 the leading expression for the common field, which makes the nonlinear correction to electric dipole, and the corresponding potential. In Sect. 3.2 the opposite case r/R 1 is considered, first, also in the leading approximation (Sect. 3.2.1). The simplifying circumstance that makes these approximations easy to handle is that it so happens that one needs, as a matter of fact, to solve only the second (following the classification of Ref. [33]) Maxwell equation, the one following from the least action principle, while the first one, the Bianchi identity, [∇ × E] = 0, is trivially satisfied. The situation becomes far more complicated in the next-to-leading approximation (r/R) 2 developed in Sect. 3.2.2. In developing the above approximations no assumption was made on whether the nonlinear scale determined by the self-coupling of the electromagnetic field is large or small as compared to r or R. Their use in the expression for the energy of the two charges as a function of the separation R in the limit R → 0, i.e. when the charges are so close to one another that R is much less than the nonlinearity scale, allows one to make a preliminary estimation confirmed by another approach to be reported in a separate publication, that, for small separation, the energy of the system of two point charges can be represented as a + bR 1/3 , where a and 1 Throughout, Greek indices span Minkowski space-time, Roman indices span its three-dimensional subspace. Boldfaced letters are threedimensional vectors, the same letters without boldfacing and index designate their lengths, except the coordinate vector x = r, whose length is denoted r. The scalar product is (r · R) = x i R i , the vector product is C = [r × R], C i = i jk x j R k . b are finite constants depending only on the two charges (in QED they include the electron mass and charge). Hence the force between two point-like charges goes to infinity following the law R −2/3 . This formula replaces, in the given nonlinear model, the Coulomb law R −2 for the force between two point charges.

Nonlinear Maxwell equations as they originate from QED
It is well known that QED is a nonlinear theory due to virtual electron-positron pair creation by a photon. The nonlinear Maxwell equation of QED for the electromagnetic field ten- produced by the classical source J μ (x) may be written as, see e.g. [25], Here L (F, G) is the effective Lagrangian (a function of the two field invariants F = (1/4)F μν F μν and G = (1/4)F μν F μν ), of which the generating functional of oneparticle-irreducible vertex functions, called effective action [34], is obtained by the space-time integration as (1) is the realization of the least action principle where the full action includes the standard classical, Maxwellian, electromag- in terms of the electric and magnetic fields, E and B.
Equation (1) is reliable only as long as its solutions vary but slowly in the space-time variable x μ , because we do not include the space and time derivatives of F and G as possible arguments of the functional Γ [A] treated approximately as local. This infrared, or local approximation shows itself as a rather efficient tool [25][26][27][28][29][30][31]. The calculation of one electronpositron loop with the electron propagator taken as solution to the Dirac equation in an arbitrary combination of constant electric and magnetic fields of any magnitude supplies us with a useful example of Γ [A], known as the E-H effective action [23,24]. It is valid to the lowest order in the finestructure constant α, but with no restriction imposed on the background field, except that it has no nonzero space-time derivatives. A two-loop expression of this local functional is also available [35].
The dynamical equation (1), which makes the "second pair" of Maxwell equations, may be completed by postulating also their "first pair" whose fulfillment allows for using the 4-vector potential A ν (x) for representation of the fields: . This representation is important for formulating the least action principle and quantization of the electromagnetic field. From it, Eq. (2) follows identically, unless the potential has the angular singularity like the Dirac string peculiar to a magnetic monopole. In the present paper we stick to Eq. (2), although its local opposite is not meaningless, as discussed in Ref. [29], where a magnetic charge is produced in nonlinear electrodynamics.
We are now going to separate the electrostatic case. This may be possible if the reference frame exists where all the charges are at rest, J 0 (x) = J 0 (r). (We denote r = x).) Then in this "rest frame" the spatial component of the current disappears, J (x) = 0, and the purely electric time-independent configuration F i j (r) = 0 would not contradict Eq. (1). With the magnetic field equal to zero, the invariant G = (E · B) disappears, too. In a theory even under a space reflection, to which class QED belongs, also we have ∂L(F,G) since the Lagrangian should be an even function of the pseudoscalar G. Then we are left with the equation for a static electric field E i = F i0 (x),

Generalities of solutions to nonlinear Maxwell equations
Equation (3) is seen to be the equation of motion stemming directly from the Lagrangian with the constant external charge J 0 (r) and the zero argument set for the second field invariant G. In the rest of the paper is based on this Lagrangian with the understanding that it may originate from QED as described above or, alternatively, be given ad hoc to define a certain model. In the latter case, if treated seriously as applied to short distances near a point charge where the field cannot be considered as slowly varying, in other words, beyond the applicability of the infrared approximation of QED outlined above, the Lagrangian (4) may be referred to as defining an extension of QED to short distances once L(F,0) is the E-H Lagrangian (or else its multi-loop specification) restricted to G = 0. It was shown in [1] that the important property of finiteness of the field energy of the point charge is guaranteed if L(F,0) in (4) is a polynomial of any power, obtained, for instance, by truncating the Taylor expansion of the H-E Lagrangian at any integer power of F. On the other hand, it was indicated in [3] that a weaker condition is sufficient: if L(F,0) grows with −F as (−F) w , the field energy is finite provided that w > 3/2. The derivation of this condition is given in [2] and in [4]. As a matter of fact a more subtle condition suffices: In the present paper we confine ourselves to the simplest example of the nonlinearity generated by keeping only quadratic terms in the Taylor expansion of the E-H Lagrangian in powers of the field invariant F, where the constant and linear terms are not kept, because their inclusion would contradict the correspondence principle, which does not admit changing the Maxwell Lagrangian L Max = −F for small fields. The correspondence principle is absorbed into the calculation of the E-H Lagrangian via the renormalization procedure. Finally, we shall be dealing with the model Lagrangian quartic in the field strength, with γ being a certain self-coupling coefficient with the dimensionality of the fourth power of the length, which may be taken as where e and m are the charge and mass of the electron, if L is chosen to be the E-H one-loop Lagrangian. We do not refer to this choice henceforward. Generalization to general Lagrangians can also be done in a straightforward way. The second (3) and the first (2) Maxwell equations for the electric field E with Lagrangian (5) are ∇ × E(r) = 0.
Denoting the solution of the linear Maxwell equations as E lin (r) we write the solution of (7), in the following way [25][26][27][28][29][30]: because ∇ · [∇ × (r)] = 0. The second Maxwell equation (9) may be conveniently written in the form to be exploited later, where the function ξ(x) is defined as a real solution to the cubic equation Its explicit form is given by the Cardano formula: We substitute (10) in the first Maxwell equation (8) to get where the prime designates the derivative with respect to the argument. Taking into account the relations Since the function ξ(x)/(3 − 2ξ(x)) does not have zeros on [0; ∞), Eq. (13) is equivalent to the equation In the center-symmetric case of a single point charge considered in [1][2][3]32], one has (r) = 0 as a solution to Eq. (17). This simplification makes the exact solution possible. The equality (r) = 0 holds as well in the axialsymmetric problem of two point charges within the approximations linear with respect to the ratios R/r or r/R to be considered in Sects. 3.1 and 3.2. In these cases it will be sufficient to represent the solution of the differential part of Eq. (7) in the form (9) setting (r) = 0 in it, then the first Maxwell equation (8) is fulfilled automatically. On the contrary, within the next order of (r/R) 2 the pseudovector function (r) is nontrivial, which makes the axial-symmetric "quadrupole-like" solution found in Sect. 3.2.2 for the field of two point-like charges more sophisticated.

Two-body problem
By the two point-charge problem we mean that the current j 0 (r) in (7) is the sum of delta-functions centered in the positions r = ±R of two charges q 1 and q 2 separated by the distance 2R (with the origin of coordinates x i placed in the middle between the charges), In what follows we shall be addressing this equation as accompanied by (8) for the combined field of two charges.
In what follows we shall refer to the field-energy density that in the present model (5), when there is only an electric field, reads The integral for the full energy of two charges converges, since it might diverge only when integrating over close vicinities of the charges. But in each vicinity the field of the nearest charge dominates, and we know from the previous publication [1] (also to be explained below) that the energy of a separate charge converges in the present model. When the charges are in the same point, R = 0, they make one charge q 1 + q 2 , whose energy converges, too.
3.1 Small separation r R between charges (dipole approximation) We shall be looking for the solution of (18) in the form where E (0) and E (1) are contributions of the zeroth and first order with respect to the ratio R/r 1, respectively. This strong inequality means that the observation point is far from the location of the two charges. So the result of our consideration in the present section will be an extension of the dipole field to the case that the point charges self-interact and interact nonlinearly with each other.
The zeroth-order term is spherical-symmetric, because it corresponds to two charges q 1 , q 2 in the same point that add to one charge q 1 + q 2 , where ξ(x) is the solution (12) of Eq. (12). Equation (8) is automatically fulfilled by the center-symmetric form (21).
The field E (0) is a nonlinear extension [1] of the standard Coulomb field q 1 + q 2 4πr 2 r r of the sum charge. Let us write the first-order term E (1) i in the following general axial-symmetric form, linear in the ratio R/r : where a and g are functions only of the scalar r, and the symmetry axis is fixed as the line passing through the two charges. Let us subject (22) to Eq. (8) and ∇ × E (1) = 0. The relation provided that the vectors r, R are not parallel. We shall see that with the ansatzes (22) and (21) with Eq. (9) can be satisfied with the choice (r) = 0: namely, we shall find the coefficient functions a, g from Eq.
(24) and then ascertain that Eq. (23) is obeyed by the solution.
The inhomogeneity in (24), satisfies the linear (γ = 0) limit of Eq. (18), and also (8). The inhomogeneity (25) is expanded in R/r as This is the standard monopole+dipole approximation with the understanding that d = (q 2 − q 1 )R is the dipole moment of the two charges, while the dots stand for the disregarded quadrupole and higher multipole contributions. The zeroth-order term satisfies the equation with the first term of the expansion (33) taken for the inhomogeneity. This is an algebraic (not differential) equation, cubic in the present model (5), solved explicitly for the field E (0) as a function of r in this case, but readily solved for the inverse function r (E (0) ) in any model, this solution being sufficient for many purposes. Even without solving it we see that for small r γ 1/4 the second term in the bracket dominates over the unity, therefore the asymptotic behavior in this region follows from (27) to be This compared to the Coulomb field (q 1 + q 2 )/(4π)r −2 weakened singularity is not an obstacle to convergence of both integrals in (20) for the proper field energy of the equivalent point charge q 1 + q 2 .
With the zeroth-order equation (27) fulfilled, we write a linear algebraic equation for the first-order correction E (1) from ( 24), to which the second, dipole part in (33) serves as an inhomogeneity, This equation is linear and it does not contain derivatives. We use (22) as the ansatz. After calculating we obtain two equations, along R and r, with the solutions a = − δq From (27) we obtain d dr Hence d dr With the help of this relation the derivative of (29) can be calculated to coincide with (30) times r. This proves Eq. (23) necessarily to satisfy the first Maxwell equation (8). By substituting Eqs. (29) and (30) in the decomposition (22) we finally have for the solution of both Maxwell equations up to O(R 2 /r 2 ).
The potential corresponding to the electric field (31) has the form where V 0 (r ) is the potential of the field of one charge [2]: where F(φ, m) = φ 0 (1−m sin 2 θ) −1/2 dθ is elliptic integral of the first kind.
The first term (31) is the nonlinear electric monopole field (21) substituting for the Coulomb field in the nonlinear problem under study, while the second term in (31) may be considered as giving a nonlinear correction to the electric dipole field. The lines of force and the equipotential curves of the latter field drawn under the choice of parameters corresponding to a strong nonlinearity are shown in Fig. 1.

Large separation r R between charges
The quantities that relate to the approximation valid at r R, dealt with in this section, will be written with a tilde to distinguish it from the corresponding quantities in Sect. 3.1 relating to the opposite approximation. Let us expand the inhomogeneity (25) to the first order in the ratio r/R (without assuming the smallness of R and r as compared to γ 1/4 ): The first term here has the clear meaning of the sum of two oppositely directed Coulomb fields produced in the point r = 0 by the two charges placed far from one another. The second one looks like a dipole field in the variable R with the equivalent "dipole moment" (q 2 + q 1 )r.
We are looking for a solution to Eq. (24) in the form of the expansion in powers of r/R, with the yet unknown dimensionless coefficients f, c and b being functions of R. The first Maxwell equation [∇×E] = 0 is satisfied by (34).
In the zeroth order we have the equation for in the form which implies that f /(4π R 2 ) = E (0) (R) is the function obtained from E (0) (r ) of the previous Sect. 3.1 by the substitutions r → R and q 2 + q 1 → q 2 − q 1 .

Leading (dipole-like) approximation
In the first order, the use of (35) turns Eq. (24) into a linear algebraic equation for E (1) (R), Calculating the second term in the right-hand side (the auxiliary electric field E(r) = (γ /2)E 2 (r)E(r)) with the ansatz (34), we obtain from (36) two equations for the components of E (1) along R and along r that determine the values Finally, where E (0) is the solution of Eq. (35) as a function of R.
The field (37) obviously satisfies the first Maxwell equation [∇ × E] = 0. By comparing (37) with the linear field of two charges in a similar approximation (33) we observe that in the zeroth-order term the difference (q 2 − q 1 )/R 2 of the two Coulomb fields in the point r = 0 has been replaced by the nonlinear field E (0) of the equivalent charge q 2 − q 1 , while in the first-order term the "dipole field" (q 2 + q 1 )r has been modified by two different factors in the terms parallel to r and R.
To be more general, note that the fields (37) and (31) turn into one another under the simultaneous replacement of the observation coordinate r by the separation R between the charges, and of the sum q 2 + q 1 of the charges by their difference q 2 − q 1 . The same symmetry under the interchange r ↔ R, q 2 + q 1 ↔ q 2 − q 1 certainly holds for the linear γ = 0 limits (33), (26) of Eqs. (37) and (31). This symmetry occurs, because the second Maxwell equations, within the approximations adopted in this section, r R, (36), and in the previous section, r R, Eq. (28), turn into each other under the transformation under consideration, while the first Maxwell equation is satisfied for both. As for the exact equation (24), this transformation maps it into a strange differential equation of a nonexistent theory.

Next-to-leading (quadrupole-like) approximation
In this section we are studying the next term E (2) , quadratic in the ratio r/R, extending the expansion (34). To this end we first extend the expansion (33) of the linear field (25) E lin (r) to include the corresponding term: Once, up to the first order in r/R, Eq. (34) satisfies the first Maxwell equation [∇ × E] = 0 automatically with any coefficients f, b, c we conclude, as we did in the previous section, that the curl [∇ × (r)] involved in (9) is zero to this order, Eq. (37) being the solution to Eq. (9) without this curl. This implies that the expansion of [∇ × (r)] starts with the quadratic term (r/R) 2 . Bearing in mind that the vector product [r × R] is the only pseudovector in our problem and that the action of ∇ lowers the power of r by one we look for the pseudovector in the form where Ω φ ((r · R)/(r R)) is a scalar function of the angle θ between the observation direction and the axis, on which the charges lie, cos θ = (r · R)/(r R). A straightforward calculation yields (we refer to and to the relation (∇ · ) = 0) obeyed by (39).
Similarly to the coefficient (q 1 − q 2 ) in the third (quadrupole) term in (38), the coefficients c(R) and d(R) are odd under the permutation q 1 ↔ q 2 . Note that the seeming singularity at q 1 = q 2 cancels from these coefficients due to the equality κ(R) = 1, which holds in this case. In the linear limit γ = 0 one also has κ(R) = 1, and c(R) and d(R) turn both into q 1 − q 2 , so that (47) turns into the last (quadrupole) term in the expansion (38) of the linear field.

Concluding remarks
We were working within the simplest nonlinear electrodynamics with the self-interaction of the fourth power of the electromagnetic field (5), which, if needed, may be thought of as resulting from the first nontrivial term of expansion of the Euler-Heisenberg effective Lagrangian in powers of its background field argument F, while the other field invariant is kept vanishing, G =0. In this case the coefficient γ , whose dimensionality is [length 4 ], which determines the strength of the nonlinearity, is expressed as (6) in terms of the electron mass and charge. Otherwise it may be considered to be arbitrary. Anyway, in our calculation the smallness of γ 1/4 as compared to the two other quantities r and R carrying the dimensionality of length was nowhere assumed.
We considered the electrostatic problem of interaction between two point-like charges q 1 and q 2 placed in the points r = ±R by solving the nonlinear Maxwell equation (18), which follows from the least action principle for the Lagrangian (4), together with the standard Bianchi identity (8). For a small separation between the charges, R r , we found the electric field (31) in the approximation, linear with respect to the ratio R/r , which serves as the nonlinear extension of the usual dipole field. The result for the corresponding scalar potential is Eq. (32). The lines of force and equipotential-curve pattern is shown in Fig. 1 in the configuration space r with the parameters chosen in such a way as to make the nonlinearity effect best pronounced. For large separation between the charges, R r, we found the electric field in the approximations (37), linear with respect to the ratio r/R, and (47), quadratic.
Using the two opposite representations (31) and (37) we can get a rough estimate for the behavior of the field energy (20), (19) in the asymptotic regime R → 0, where the two point charges approach each other infinitely close. According to that estimate, in this regime the energy of the system of two point charges can be represented as where a and b are finite constants depending only on the charges q 1 and q 2 , and on the self-coupling constant γ . The R-independent term a is the self-energy of the united pointlike charge with the value q 1 + q 2 . This is finite, as established in [1]. The behavior (48) is rigorously confirmed following a quite different procedure to be published elsewhere. Although the energy is finite in the limit R = 0, the force F between the two charges defined as the derivative of the energy with respect to the distance is weakly infinite: This formula replaces, in the given nonlinear model, the Coulomb law F ∼ R · R −3 for the force between two point charges. The power 2/3 here is determined by the power 2 in the self-interaction in (5).