Nonlinear Electrodynamics, Regular Black Holes
And Wormholes
Abstract
We consider spherically symmetric configurations in general relativity, supported by nonlinear electromagnetic fields with gaugeinvariant Lagrangians depending on the single invariant . Static black hole (BH) and solitonic solutions are briefly described, both with only an electric or magnetic charge and with both nonzero charges (the dyonic ones). It is stressed that only pure magnetic solutions can be completely nonsingular. For dyonic systems, apart from a general scheme of obtaining solutions in quadratures for an arbitrary Lagrangian function , an analytic solution is found for the truncated BornInfeld theory (depending on the invariant only). Furthermore, considering spherically symmetric metrics with two independent functions of time, we find a natural generalization of the class of wormholes found previously by Arellano and Lobo with a timedependent conformal factor. Such wormholes are shown to be only possible for some particular choices of the function , having no Maxwell weakfield limit.
Received Day Month Year
Revised Day Month Year
Keywords: General relativity, nonlinear electrodynamics, spherical symmetry, solitons, black holes, wormholes, exact solutions
PACS numbers: 04.70.Bw, 04.20.Jb, 04.20.Gz
1 Introduction
Nonlinear electrodynamics (NED) appeared in the 1930s with Born and Infeld’s effort to remove the central singularity of a point charge and the related energy divergence by generalizing Maxwell’s theory , and the version of NED put forward by Heisenberg and Euler, motivated by particle physics . It was further extended by Plebanski in the framework of special relativity, including an arbitrary function of the electromagnetic field invariants.
The interest in NED in modern studies is to a large extent motivated by the discovery that some kinds of NED appear as limiting cases of certain models of string theory . It is also clear that the real electromagnetic field should lose its linearity at high energies due to interactions with other physical fields, and NED theories may be considered as a simplified phenomenological description of these interactions. On the other hand, NED as a possible material source of gravity is able to create various nonsingular geometries of interest, in particular, regular black holes (BHs) and starlike or solitonlike configurations in the framework of general relativity (GR) and alternative theories.
Among such models, the simplest are spherically symmetric ones, where the only possible kinds of electromagnetic fields are radial electric and radial magnetic fields. Such solutions are widey discussed in the literature, beginning probably with the paper by Pellicer and Torrence where a general static solution was obtained for configurations with an electric field only. In Ref. ?, a nogo theorem was proved showing that if NED is specified by a Lagrangian function (where , and is the Maxwell tensor), there is no such function having a Maxwell weakfield limit ( as ) that a static, spherically symmetric solution of GR with an electric field has a regular center. This theorem was further extended to static dyonic configurations, with both electric and magnetic fields , and it was further shown that in numerous electric solutions describing configurations with or without horizons (that is, BH or solitonic ones), having a regular center and a ReissnerNordström (RN) behavior at large radii , there are different Lagrangian functions at large and small : at large we have whereas at small the theory is strongly nonMaxwell ( but , in agreement with the nogo theorem).
Meanwhile, purely magnetic regular configurations, both BH and solitonic ones, are possible and are readily obtained under the condition as . Electric models with the same regular metrics can be obtained from the magnetic ones using the socalled FP duality (not to be confused with the familiar electricmagnetic duality in Maxwell’s theory) that connects solutions with the same metric corresponding to different NED theories. Unlike the magnetic solutions, the electric ones suffer serious problems connected with multivaluedness of and a singular behavior of the electromagnetic fields on the branching surfaces .
Many results of interest were obtained since then, for a brief review see, e.g., Ref. ? and references therein. Among them let us point out a description of models with a kind of phase transition allowing one to circumvent the above nogo theorem , an extension of static, spherically symmetric NED solutions to GR with a nonzero cosmological constant , thermodynamic properties of regular NED BHs , cylindrically and axially symmetric regular GR/NED configurations and evolving wormhole models, the stability properties of NED BHs, and quantum effects in their fields.. One can also mention numerous studies of special cases of electric and magnetic solutions, their potential observational properties like gravitational lensing, particle motion and matter accretion in the fields of NED BHs, their counterparts in scalartensor, and multidimensional theories of gravity, inclusion of dilatonlike interactions, nonAbelian fields, constructions with thin shells, etc., but the corresponding list of references would be too long. The subject probably deserves a comprehensive review.
In this paper we discuss some recent progress concerning two subjects in the framework of NED coupled to GR: static, spherically symmetric dyonic solutions and spherically symmetric evolving wormholes. The dyonic solutions are inevitably singular at the center, as follows from the nogo theorem , but they are of interest as the first examples of this kind of solutions.. For completeness and comparison, pure electric and magnetic solutions are also briefly described.
As to wormholes as twoway tunnels or shortcuts between different universes or different, otherwise distant regions of the same universe, their possible existence and properties are widely discussed, see, e.g., Refs. ?–? for reviews. Wormholes are of interest not only as a perspective “means of transportation” but also as possible time machines or accelerators . Spherical symmetry is a natural simple framework for wormhole geometry, and most of known exact wormhole solutions in GR and alternative theories of gravity (e.g., Refs. ?–? and many others) are static, spherically symmetric. However, NED as a source of gravity in GR cannot support static wormholes because it does not provide the necessary violation of the Null Energy Condition (NEC). Only by considering evolving conformally static spacetimes it has been possible to obtain some examples of NED/GR wormhole solutions . In this paper, the approach developed in Refs. ?–? is extended to a more general class of timedependent metrics, containing two functions of time and somewhat similar to KantowskiSachs cosmologies.
After presenting some general relations valid for both static and timedependent NED/GR configurations (Section 2), in Section 3 we briefly describe all three types of static solutions: magnetic, electric and dyonic ones, following the previous papers, Refs. ? and ?. In Section 4 we obtain and discuss a new class of nonstatic spherically symmetric wormhole solutions, containing those of Arellano and Lobo as a special case. Section 5 is a conclusion.
2 Basic equations. FP duality
We start with the action
(1) 
where is the Ricci scalar, is an arbitrary functions, and units are used with . The Einstein equations can be written, as usual, in two equivalent forms
(2) 
where is the stressenergy tensor (SET), which in the theory (2) is given by ()
(3) 
The metric is taken in the general spherically symmetric form
(4) 
The only nonzero components of compatible with spherical symmetry are (a radial electric field) and (a radial magnetic field). The Maxwelllike equations and the Bianchi identities for the dual field lead to
(5) 
where and are the electric and magnetic charges, respectively.
Accordingly, the only nonzero SET components are
(6) 
and the invariant is the difference , where
(7) 
and being the absolute values of the radial electric field strength and magnetic induction, respectively, measured by an observer at rest in the reference frame under consideration.
The SET (2) has the important properties and ; the latter means the absence of radial energy flows, which is in turn related to the absence of electromagnetic monopole radiation. Taken together, these two properties define a kind of matter sometimes called Dymnikova’s vacuum , its evident vacuumlike property is that its SET structure is insensitive to any transformations of the coordinates and ; in other words, all reference frames moving in the radial direction with any velocities relative to each other are comoving to this kind of matter. Thus a spherically symmetric electromagnetic field described by NED is a particular form of Dymnikova’s vacuum.
NED with a Largangian function is known to admit a dual representation obtained by a Legendre transformation : one defines the tensor with its invariant and considers the Hamiltonianlike quantity
(8) 
as a function of ; then can be used to specify the whole theory. One has then
(9) 
with . The SET in terms of and then reads
(10) 
In a spherically symmetric spacetime with the metric (4), Eqs. (5) are rewritten in the P framework as
(11) 
We can also introduce the quantities and similar to (2):
(12) 
so that , and then the SET (10) takes the form
(13) 
Comparing (2) and (2), one can see that they coincide up to the substitutions
(14) 
where is the Hodge dual of , so that . The coincidence of the SETs means that the sets of metrics satisfying the Einstein equations (2) also coincide. It is the FP duality described in Ref. ?, which was formulated there for static, spherically symmetric systems and is extended here to nonstatic ones. It should be stressed that this duality connects configurations with the same metric but in it different NED theories. An evident exception is the Maxwell theory, where , and the FP duality turns into the conventional electricmagnetic duality.
If, however, one speaks of different formulations of the same theory, it turns out that its and formulations are not always equivalent: it is only the case if is a monotonic function, see below.
3 Static systems
If the metric is static, so that the functions depend on only, for our system it is reasonable to choose the “Schwarzschild” radial coordinate, , so that and . Then, due to the equality , we have from the Einstein equations , leading to (with a proper choice of the time scale), so that
(15) 
and the Einstein equation then leads to
(16) 
where is the energy density, and is called the mass function. It is a general relation, but it is only a part of a possible complete solution: the latter requires a knowledge of and both electric and magnetic fields.
3.1 Pure magnetic and electric solutions
Pure magnetic solutions () are obtained from Eqs. (2) and (16) quite easily. Indeed, if is specified, then, since now , the function is known from (2), and the metric function is found by integration in (16). If, on the contrary, is known (or chosen at will), then is found from (16), and is restored since . A regular center requires at small (and this leads to as ), asymptotic flatness requires , where is the Schwarzschild mass, and, if , asymptotically (A)dS solutions are obtained by adding to in (16). It is an easy way to construct regular magnetic BH and solitonic solutions, used by many authors, probably beginning with Ref. ?.
A general feature of all such regular solutions is that as both and . Moreover, the mass term contributes negatively to as long as whereas contributes positively as long as (which is necessary for getting ). Thus in regular solutions inevitably has a minimum, and the value of at this minimum — hence the existence of horizons as zeros of — depends on a relationship between and . If the mass is fixed, then at small the minimum of is negative since the solution is close to Schwarzschild’s in the whole space except small radii, . In this case each regular solution has two horizons, one of which is Schwarzschildlike, close to , and the other exists since it is necessary to return to at small to reach at . At large the mass term becomes negligible at all radii except for the asymptotic region, and a minimum of should become positive, leading to a solitonic solution. Therefore, some value of should be critical, corresponding to a double zero of , hence a single extremal horizon. This general picture is observed in all existing examples of regular static, spherically symmetric NEDGR solutions. We here illustrate it with the behavior of in the example of Ref. ? in which
(17) 
The behavior of is depicted in Fig. 1 for three values of leading to qualitatively different geometries. Their causal structures and CarterPenrose diagrams are the same as those well known for RN spacetimes, with the important difference that now the lines correspond to a regular center.
Pure electric solutions (, ) are obtained in a similar way by using the Hamiltonian form of NED, see Eqs. (9)–(2). We have simply , and specifying , we directly find and from (16). On the contrary, specifying , from (16) we determine .
A regular center requires a finite limit of as . However, in any regular asymptotically flat (or (A)dS) solution at both and , so inevitably has at least one maximum at some , violating the monotonicity of , necessary for equivalence of the and frameworks. As shown in Ref. ?, at an extremum of the Lagrangian function suffers branching, its plot forming a cusp, and different functions correspond to and . Another kind of branching occurs at extrema of , if any, and the number of Lagrangians on the way from infinity to the center equals the number of monotonicity ranges of . An example of such behavior is shown in Fig. 2, it corresponds to the electric solution with a regular center from Ref. ?, in which
.
The Hamiltonian framework might seem to be not worse than the Lagrangian one, even though the latter directly follows from the least action principle. However, it turns out that at the electromagnetic field is singular: the effective metric in which NED photons move along null geodesics, is singular at extrema of , photons are there infinitely blueshifted and can create a curvature singularity by back reaction on the metric. Thus any regular electric solution not only fails to correspond to a fixed Lagrangian but has other important undesired features.
We conclude that each choice of gives rise to both electric and magnetic solutions related by FP duality. In these solutions the NED theories are different, and if describes a geometry with a regular center and a RN asymptotic at large , only the magnetic solution is really regular.
3.2 Dyonic configurations
Let us now assume that both and are nonzero. The difficulty in finding the solutions is that now (or alternatively ) is not known explicitly. Using, for definiteness, the Lagrangian formulation of the theory, we have
(18) 
Comparing the expressions for from (2) and (16), we can write ()
(19) 
Let us specify the theory by choosing . Then Eq. (18) can be considered as either (A) an equation (in general, transcendental) for the function or (B) an expression of as a function of .
In case A, if can be found explicitly, integration of Eq. (19) (equivalent to (16)) gives the metric function , and this completes the solution.
Scheme B leads to a solution in quadratures in terms of which can be chosen as a new radial coordinate. Indeed, assuming that and are known and monotonic, so that and , Eq. (19) can be rewritten as
(20) 
(the subscript denotes ). Since the r.h.s. of (20) is known, it is straightforward to find and to pass on to the coordinate in the metric. This gives us a general scheme of finding dyonic solutions under the above conditions.
Consider two examples using scheme A. The first example is used to verify the method: the Maxwell theory, , . Then from Eq. (19) we obtain , whence and
(21) 
that is, the dyonic RN solution, as should be the case.
In the second example we assume that Eq. (18) is linear with respect to , which unambiguously leads to the truncated BornInfeld Lagrangian,
(22) 
(the full BornInfeld Lagrangian would also contain the invariant ).
Indeed, Eq. (18) is linear in only if , . Integration gives . For a Maxwell behavior () at small we put , . Denoting , we arrive at Eq. (22), and
(23) 
where
The simplest solution corresponds to the special case of a selfdual electromagnetic field, , whence , , and as in the Maxwell theory. This leads to and the dyonic RN metric with given by (21).
In the general case , and are expressed in terms of the Appell hypergeometric function :
(24)  
where , and .
Some features of the solution should be noticed. As expected, the center is singular in accord with the above nogo theorem. At large the quantities . and the energy density decay as , and the solution is asymptotically flat and approximately RN. At small , tends to a finite limit while . In a pure electric solution we also have , while in a pure magnetic one . Thus all solutions are singular at the center , but the in the pure magnetic solution the singularity (existing since the function (22) does not tend to a finite limit as ) is milder.
4 Dynamic wormholes
Apart from BH and solitonic configurations with a regular center, there can exist regular objects having no center at all, namely, wormholes and some classes of regular BHs, including the socalled black universes (BHs in which beyond the event horizon, instead of a singularity, there is an expanding universe), see, e.g., Refs. ?–?. In the static case all of them require violation of the Null Energy Condition that states . Since for the SET (2) such a difference is zero, this condition (though marginally) is observed by NED, so static wormholes are manifestly impossible in the theory (2). However, as shown by Arellano and Lobo , wormhole solutions can be obtained if we consider evolving spacetimes. We will try to extend their finding, considering a more general class of timedependent metrics.
4.1 Equations
In Ref. ? and the subsequent discussion , the metric was chosen to be static times a timedependent conformal factor. Let us assume a more general metric,
(25) 
and NED as the matter source of gravity as given in (2). The function can be absorbed by redefinition of and by redefinition of , with the result
(26) 
that substantially simplifies the Einstein equations. For this metric the nonzero components of the Ricci tensor are
(27)  
(28)  
(29)  
(30) 
The electromagnetic field equations lead to the same relations (5) and (2) as in the static metric, where now , preserving its geometric meaning of the spherical radius. As before, the component of the SET is zero, and the corresponding Einstein equation takes the form
(31) 
where dots denote and primes . Assuming , and dividing Eq. (31) by , we separate the variables and, without loss of generality, obtain
(32) 
where is the separation constant. Next, since , we have the equation , where, excluding and according to (32), we again separate the variables, obtaining
(33) 
with the separation constant assumed to be positive since we are seeking solutions with a minimum () of the function , describing a throat in 3D spatial sections .
Thus and are determined by the equations
(34)  
(35) 
whose first integrals are easily found. Specifically, for we have
(36) 
and for
(37) 
where are integration constants.
4.2 Geometry
Solutions to Eqs. (36) and (37) completely determine the metric under the ansatz (25) or (26) for any kind of matter whose SET satisfies the conditions and . As already mentioned, these conditions define the socalled Dymnikova vacuum , for which any reference frame in radial motion is comoving.
In all cases under consideration, by construction, any 3D section of spacetime contains a throat defined as a minimum of and hence of at any given time instant, whereas the global features of these spacetimes depend on the constants involved.
Some immediate observations on the time dependence of .follow directly from Eq. (35). Thus, a regular minimum of (that is, a bounce in the time evolution of ) is impossible since Eq. (35) leads to at points where and . Next, a finite limit of as is also impossible because at such a minimum the l.h.s. would be zero while the r.h.s. is finite. So the only possible way of regular evolution of is to begin or end with at infinite proper time, which is not completely excluded.
Further integration of Eqs. (36) for leads in the general case to the hypergeometric function :
(38) 
while integration of Eq. (37) for involves the error function Erf:
(39) 
Both solutions (4.2) and (4.2) are rather hard for further investigation, let us therefore restrict ourselves to some simple special cases.
Example 1: . In this special case we obtain and thus restore the isotropic nature of spacetime evolution considered in Refs. ?–?. The metric takes the form
(40) 
Assuming in Eqs. (36) and (37), we obtain the following solutioon under a proper choice of the zero points of the coordinates and :
(41) 
where we have denoted . It can be directly verified that this solution coincides with the one considered in Refs. ?–?, but in another, more preferable parametrization: in (41) the coordinate is proper time, and is proportional to the proper distance at any fixed time .^{a}^{a}aOur coordinates and the coordinates chosen by Arellano and Lobo in Ref. ? are connected by the relations , . The coordinate of Refs. ?–? should not be confused with the quantity used here, having the geometric meaning of the radius of spheres . Thus the spatial sections are evidently of wormhole nature but not asymptotically flat (the latter would require at large whereas here ), while the time evolution is different for different : at it is semiinfinite, and at it occupies a finite period between two zeros of . The zeros of correspond to bigbang type singularities.
Example 2: . In this case there is no expansion or contraction in the radial direction, and the coordinate is the true radial length in the reference frame used. The metric has the form
(42) 
Integration of Eqs. (36) gives:
(43) 
where we have denoted , (both and should be positive for Eqs. (36) to be meaningful). The properties of this solution are to a large extent the same as those of the branch in Example 1, but a significant difference is the emergence of that grows together with at large somewhat similarly to the static antide Sitter metric.
Example 3: . In this case Eqs. (36) are also easily integrated in elementary functions, but the first equation implies , that is, no throats are possible in the spacial sections of this spacetime. Since we are seeking wormhole solutions, we do not consider this case any more.
It is clear that at the throat in all cases the gradient of the spherical radius treated as a scalar function in 2D spacetime parametrized by and is timelike since while . Thus at least a certain neighborhood of the throat is a socalled Tregion, i.e., a region where may be chosen as a temporal coordinate (as happens, e.g., inside a Schwarzschild horizon), and the geometry is thus cosmological in nature, resembling KantowskiSachs spherically symmetric cosmologies . Does it mean that a Tregion covers the whole spacetime? An answer apparently depends on the parameters involved. For example, for the solution (41) we have
(44) 
Thus in models with we have , i.e., the gradient of is timelike in the whole spacetime, it is a global Tregion. On the contrary, in the case we have on the surfaces (apparent horizons) defined by the relation , and there are Rregions at , in this sense the spacetime is of black hole type.
4.3 Electromagnetic fields
Let us return to NED as the matter source of gravity. It is easy to notice that the suitable form of depends on the values of and . As follows from the expressions (2), the magnetic field is regular at all finite , while the electric field can tend to infinity not only where but also at zeros of the derivative if any.
Pure magnetic solutions are the simplest, just as in the static case. Indeed, assuming and , we have , and using (36) and (37), it is easy to verify that is expressed in terms of rather than separately in terms of and :
(45)  
(46) 
where is an arbitrary constant introduced for dimensional considerations, and its choice can is related to the arbitrariness of the constant . Since , it is straightforward to find by substituting , where :
(47)  
(48) 
where is introduced similarly to . It is the full set of NED Lagrangians that lead to magnetic solutions of GR/NED equations with the metric (26). We see that none of these possess a Maxwell behavior at small .
The corresponding electric solutions with the same metric pertain to other versions of NED, as follows from the FP duality described in Section 2. For the same metric (26) we now have , which, as we saw, is a function of , and since now , we know the function which is the same as in (47) and (48), where should be replaced with , and with , so that
(49)  
(50) 
Then, both and can in principle be found according to (14) as well as the function , and, as in the static case, there should be as many different Lagrangians in different parts of spacetime as is the number of monotonicity ranges of . From (49) and (50) we can find in terms of (either as or from the Einstein equations giving, for electric solutions, ), so that
(51)  
(52) 
and substituting . The function is also easily found since , and is known. However, it is in general hard to find and in an explicit form since finding the inverse of requires solving a transcendental equation.
Let us look what happens in the above two simple examples, and .
Example 1: , with the metric (40), (41). For magnetic wormholes, Eq. (47) gives
(53) 
a Lagrangian of the form previously used in a number of studies, see, e.g., Refs. ?, ? and references therein. For magnetic wormholes supported by NED in 2+1 dimensions such a Lagrangian was considered in Ref. ?.
For electric wormholes with the same metric the function has the form (53) with the substitution , see above. Assuming , we obtain , hence .^{b}^{b}bAccording to (2), surprisingly, the physically meaningful electric field strength is constant in this highly inhomogeneous and nonstatic spacetime. It has been asserted (e.g., in Ref. ?) that the electric field is singular at the throat in this solution; but this assertion applies to the parametrizationdependent quantity in the particular coordinates used there rather than . This strange result still does not immediately lead to a contradiction because by Eq. (51) we have also . So this solution corresponds to a constant at which the function has a certain particular value, but at other values of the function is not defined. However, since , the derivative has different values at different spacetime points, and since is fixed, it is inconsistent with any welldefined function . This makes the Lagrangian formulation of the theory illdefined for this electric solution.
Example 2: , with the metric (42), (43). For magnetic wormholes, Eq. (47) gives
(54) 
in full similarity with (53), up to the choice of the integration constants . So this is one more kind of magnetic wormhole geometry supported by NED with such Lagrangians. As to electric wormholes with the same geometry, the function for them is again illdefined for the same reasons as in Example 1.
For dyonic solutions with the same metrics, the problem of finding suitable NED formulations is more difficult and will not be considered here.
5 Concluding remarks
1. We have recalled the problem of obtaining regular static, spherically symmetric black hole solutions in GR with NED as a source of gravity and stressed the fundamental distinction between electric and magnetic fields, connected with the absence of the familiar duality inherent to the Maxwell electrodynamics. In the latter, knowing, for example, a magnetic field, it is straightforward to obtain any mixture of electric and magnetic fields by duality rotations. Unlike that, in the NED framework, FP duality connects only pure electric and pure magnetic configurations with the same metric but sourced by different versions of NED. It also happens that only pure magnetic solutions lead to completely regular configurations flat at infinity since their electric counterparts exhibit undesired features at some intermediate radii. Next, it turns out that finding dyonic configurations in NED/GR is quite a nontrivial problem, and only some special examples of such solutions are known.
2. In our search for evolving wormhole solutions, we have obtained a new family of geometries, which are in general not conformally static, except for a special case where they reduce to the known solutions of Ref. ?. Under the metric ansatz (26), these geometries represent the general solution of GR for Dymnikova’s vacuum, defined by its transformation properties in the () 2D subspace, similar to the properties of the cosmological constant in the whole 4D spacetime.^{c}^{c}cIn a sense, it is half of the whole set of solutions because we chose as the separation constant in Eq. (33) in order to obtain a minimum of . Solutions with the other sign of the separation constant certainly exist as well. A general feature of these geometries is the existence of a throat in their spatial sections. However, their time evolution in general contains cosmologicaltype singularities. It has been shown that they can neither have a bounce of nor begin or end with its finite value as . It is not excluded that they have a “remote singularity”, that is, a zero value of as , but even this opportunity may be ruled out in a future study. The existence of singularities at finite times is confirmed by examples given by Eqs. (41) and (43).
3. All that was independent of the particular choice of Dymnikova’s vacuum. If we apply NED for its implementation, it is rather easy to obtain these solutions with pure magnetic fields, but the choice of suitable NED Lagrangians is quite narrow and is expressed in Eqs. (47) and (48). None of these Lagrangians have a Maxwell weak field limit.
4. We have extended the FP duality, previously formulated for static systems, to timedependent ones and used it for comparison between the electric and magnetic evolving wormhole solutions. The electric solutions are obtained in the “Hamiltonian” formulation of NED in the same way as magnetic ones in its Lagrangian formulation. However, it is in general hard to obtain a Lagrangian for a given electric solution due to a necessity to deal with transcendental equations. Moreover, in particular examples the Lagrangian formulation is illdefined for electric solutions.
To conclude, the NED/GR system exhibits some unusual and unexpected features and deserves a further study and discussion.
Acknowledgments
I thank Milena Skvortsova, Sergei Bolokhov and Sergei Rubin for helpful discussions. The work was partly performed within the framework of the Center FRPP supported by MEPhI Academic Excellence Project (contract No. 02.a03.21.0005, 27.08.2013). The work was also funded by the RUDN University Program 5100 and by RFBR grant 160200602.
References
 1. M. Born and L. Infeld, Foundations of the new field theory, Proc. R. Soc. Lond. 144, 425 (1934).
 2. W. Heisenberg and H. Euler, Folgerungen aus der Diracschen Theorie des Positrons, Z. Phys. 98, 714 (1936).
 3. J. Plebanski, NonLinear Electrodynamics — A Study (C.I.E.A. del I.P.N., Mexico City, 1966).
 4. N. Seiberg and E. Witten, String theory and noncommutative geometry, J. High Energy Phys. 09, 032 (1999); hepth/9908142.
 5. A. A. Tseytlin, BornInfeld action, supersymmetry and string theory, hepth/9908105.
 6. R. Pellicer and R. J. Torrence, Nonlinear electrodynamics and general relativity, J. Math. Phys. 10, 17+18 (1969).

7.
K. A. Bronnikov and G. N. Shikin,
On the ReissnerNordström problem with a nonlinear electromagnetic field,
in Classical and Quantum Theory of Gravity (Trudy IF AN BSSR, Minsk, 1976),
p. 88 (in Russian);
K. A. Bronnikov, V. N. Melnikov, G. N. Shikin, and K. P. Staniukovich, Scalar, electromagnetic, and gravitational fields interaction: particlelike solutions, Ann. Phys. (N.Y.) 118, 84 (1979).  8. K. A. Bronnikov, Regular magnetic black holes and monopoles from nonlinear electrodynamics, Phys. Rev. D 63, 044005 (2001); grqc/0006014.
 9. K. A. Bronnikov, Comment on ‘Regular black hole in general relativity coupled to nonlinear electrodynamics’ , Phys. Rev. Lett. 85, 4641 (2000).
 10. K. A. Bronnikov, Dyonic configurations in nonlinear electrodynamics coupled to general relativity, Grav. Cosmol. 23, 343 (2017); arXiv: 1708.08125.
 11. A. Burinskii and S. R. Hildebrandt, New type of regular black holes and particlelike solutions from nonlinear electrodynamics, Phys. Rev. D 65, 104017 (2002); hepth/0202066 .
 12. J. Matyjasek, D. Tryniecki, and M. Klimek, Regular black holes in an asymptotically de Sitter universe, Mod. Phys. Lett. A 23, 3377 (2009); arXiv: 0809.2275.
 13. N. Bretón, Smarr’s formula for black holes with nonlinear electrodynamics, Gen. Rel. Grav. 37, 643 (2005); grqc/0405116.
 14. S. I. Kruglov, Asymptotic ReissnerNordström solution within nonlinear electrodynamics, Phys. Rev. D 94, 044026 (2016); arXiv: 1608.04275.
 15. ZhongYing Fan and Xiaobao Wang, Construction of regular black holes in general relativity, Phys. Rev. D 94, 124027 (2016); arXiv: 1610.02636.
 16. K. A. Bronnikov, G. N. Shikin, and E. N. Sibileva, Selfgravitating stringlike configurations from nonlinear electodynamics, Grav. Cosmol. 9, 169 (2003); grqc/0308002.
 17. C. Bambi and L. Modesto, Rotating regular black holes, Phys. Lett. B 721, 329 (2013); arXiv: 1302.6075.
 18. I. Dymnikova and E. Galaktionov, Regular rotating electrically charged black holes and solitons in nonlinear electrodynamics minimally coupled to gravity, Class. Quantum Grav. 32, 165015 (2015); arXiv: 1510.01353.
 19. A. V. B. Arellano and F. S. N. Lobo, Evolving wormhole geometries within nonlinear electrodynamics, Class. Quantum Grav. 23, 5811 (2006); grqc/0608003.
 20. Ch. G. Boehmer, T. Harko, and F. S. N. Lobo, Conformally symmetric traversable wormholes, Phys. Rev. D 76, 084014 (2007); arXiv: 0708.1537.
 21. A. V. B. Arellano, N. Bretón, and R GarciaSalcedo, Some properties of evolving wormhole geometries within nonlinear electrodynamics, Gen. Rel. Grav. 41, 2561 (2009); arXiv: 0804.3944.
 22. C. Moreno and O. Sarbach, Stability properties of black holes in selfgravitating nonlinear electrodynamics, Phys. Rev. D 67, 024028 (2003); grqc/0208090.
 23. N. Bretón, Stability of nonlinear magnetic black holes, Phys. Rev. D 72, 044015 (2005); hepth/0502217.
 24. Jin Li, Kai Lin and Nan Yang, Nonlinear electromagnetic quasinormal modes and Hawking radiation of a regular black hole with magnetic charge, Eur. Phys. J. C 75, 131 (2015); arXiv: 1409.5988.
 25. W. Berej and J. Matyjasek, Vacuum polarization in the spacetime of charged nonlinear black hole, Phys. Rev. D 66, 024022 (2002); grqc/0204031.
 26. J. Matyjasek, P. Sadurski, and D. Tryniecki, Inside the degenerate horizons of regular black holes, Phys. Rev. D 87, 124025 (2013); arXiv: 1304.6347.
 27. M. Visser, Lorentzian Wormholes: From Einstein to Hawking (American Institute of Physics, New York, 1995)
 28. F. S. N. Lobo, Exotic solutions in General Relativity: Traversable wormholes and warp drive spacetimes. In: Classical and Quantum Gravity Research, p. 178 (Nova Sci. Pub., 2008); ArXiv: 0710.4474.
 29. K A. Bronnikov and S. G. Rubin, Black Holes, Cosmology, and Extra Dimensions (World Scientific, 2012).
 30. F. S. N. Lobo (Editor). Wormholes, Warp Drives and Energy Conditions (Springer, 2017).
 31. M .S. Morris and K. S. Thorne, Wormholes in spacetime and their use for interstellar travel: A tool for teaching general relativity. Am. J. Phys. 56, 395 (1988).
 32. K. A. Bronnikov, Scalartensor theory and scalar charge. Acta Phys. Pol. B 4, 251 (1973).
 33. H. Ellis, Ether flow through a drainhole — a particle model in general relativity. J. Math. Phys. 14, 104 (1973).
 34. K. A. Bronnikov and A. M. Galiakhmetov, Wormholes and black universes without phantom fields in EinsteinCartan theory. Phys. Rev. D 94, 124006 (2016); arXiv: 1607.07791.
 35. G. Dotti, J. Oliva, and R. Troncoso, Static wormhole solution for higherdimensional gravity in vacuum. Phys. Rev. D 75, 024002 (2007); hepth/0607062.
 36. T. Harko, F.S.N. Lobo, M.K. Mak, and S.V. Sushkov, Gravitationally modified wormholes without exotic matter, Phys. Rev. D 87, 067504 (2013); arXiv: 1301.6878.
 37. K. A. Bronnikov and S.W. Kim, Possible wormholes in a brane world. Phys. Rev. D 67, 064027 (2003), grqc/0212112.
 38. K. A. Bronnikov and M. V. Skvortsova, Wormholes leading to extra dimensions. Grav. Cosmol. 22, 316 (2016); arXiv: 1608,04974.
 39. I.G. Dymnikova, Vacuum nonsingular black hole. Gen. Rel. Grav. 24, 235 (1992).
 40. K. A. Bronnikov, S.W. Kim and M. V. Skvortsova, The Birkhoff theorem and string clouds, Class. Quantum Grav. 33, 195006 (2016); arXiv: 1604.04905.
 41. I. H. Salazar, A. Garcia and J. Plebanski, Duality rotations and type D solutions to Einstein equations with nonlinear electromagnetic sources. J. Math. Phys. 28, 2171 (1987).
 42. L. Balart and E. C. Vagenas, Regular black hole metrics and the weak energy condition Phys. Lett. B 730, 14 (2014); arXiv: 1401.2136.
 43. M. Novello, V. A. de Lorenci, J. M. Salim, and R. Klippert, Geometrical aspects of light propagation in nonlinear electrodynamics, Phys. Rev. D 61, 045001 (2000).