Rotating Melvin-like universes and wormholes in general relativity K. A. Bronnikov,a,b,c,1 V. G. Krechet,d and V. B. Oshurkod a Center for Gravitation and Fundamental Metrology, VNIIMS, Ozyornaya ul. 46, Moscow 119361, Russia b Peoples’ Friendship University of Russia (RUDN University), ul. Miklukho-Maklaya 6, Moscow 117198, Russia c National Research Nuclear University “MEPhI” (Moscow Engineering Physics Institute), Moscow, Russia d Moscow State Technological University “Stankin”, Vadkovsky per. 3A, Moscow 127055, Russia We find a family of exact solutions to the Einstein-Maxwell equations for rotating cylindrically symmetric distributions of a perfect fluid with the equation of state p = wρ ( w < 1), carrying a circular electric current in the angular direction. This current creates a magnetic field along| | the z axis. Some of the solutions describe geometries resembling that of Melvin’s static magnetic universe and contain a regular symmetry axis, while some others (in the case w> 0) describe traversable wormhole geometries which do not contain a symmetry axis. Unlike Melvin’s solution, those with rotation and a magnetic filed cannot be vacuum and require a current. The wormhole solutions admit matching with flat-space regions on both sides of the throat, thus forming a cylindrical wormhole configuration potentially visible for distant observers residing in flat or weakly curved parts of space. The thin shells, located at junctions between the inner (wormhole) and outer (flat) regions, consist of matter satisfying the Weak Energy Condition under a proper choice of the free parameters of the model, which thus forms new examples of phantom-free wormhole models in general relativity. In the limit w 1, the magnetic field tends to zero, and the wormhole model tends to the one obtained previously, where→ the source of gravity is stiff matter with the equation of state p = ρ. 1 Introduction the possible inclusion of rotation, avoiding compli- cations inherent to the more realistic axial sym- Cylindrical symmetry is the second (after the metry, not to mention the general nonsymmet- spherical one) simplest space-time symmetry mak- ric space-times. Accordingly, a great number of ing it possible to obtain numerous exact solutions in exact stationary (assuming rotation) solutions to general relativity and its extensions, characterizing the Einstein equations are known, with various local strong gravitational field configurations. One sources of gravity: the cosmological constant [5–9]; of the motivations of studying cylindrically sym- scalar fields with different self-interaction poten- metric configurations is the possible existence of tials [10–12]; rigidly or differentially rotating dust such linearly extended structures as cosmic strings [13–15], dust with electric charge [16] or a scalar arXiv:2007.01145v1 [gr-qc] 2 Jul 2020 as well as the observed cosmic jets. A large num- field [17], fluids with different equations of state, ber of static cylindrically symmetric solutions have above all, perfect fluids with p = wρ, w = const (in been obtained and studied since the advent of gen- usual notations) [18–22], some kinds of anisotropic eral relativity, including vacuum, electrovacuum, fluids [23–26] etc., see also references therein and perfect fluid and others, see reviews in [1, 2] and the reviews [1,2]. references therein. In this paper we obtain rotating counterparts Melvin’s famous solution to the Einstein- of the static cylindrically symmetric solutions to Maxwell equations, an “electric or magnetic geon” the Einstein-Maxwell equations with a longitudi- [3], is a completely regular static, cylindrically sym- nal magnetic field. It turns out that such a field metric solution with a longitudinal electric or mag- cannot exist without a source in the form of an netic field as the only source of gravity. It is a electric current, and we find solutions where such special case from a large set of static cylindrically a source is a perfect fluid with p = wρ. Many symmetric Einstein-Maxwell fields, see more details features of these solutions are quite different from in [2,4]. those of the static ones, in particular, their common An important distinguishing feature of cylindri- feature is the emergence of closed timelike curves cal symmetry as compared to the spherical one is at large radii. Also, there is a family of wormhole 1e-mail: [email protected] solution that do not have a symmetry axis but con- 2 tain a throat as a minimum of the circular radius. has the meaning of the angular velocity of a con- As in our previous studies [10,11,27,28], we try to gruence of timelike curves (vorticity) [10,29,30], make such wormholes potentially observable from 1 − ′ − − spatial infinity by joining outer flat-space regions ω = (E e 2γ) eγ β α. (5) at some junction surfaces and verify the validity of 2 the Weak Energy Condition for matter residing on Furthermore, in the reference frame comoving to these surfaces. matter as it rotates in the azimuthal (ϕ) direc- The structure of the paper is as follows. Section 3 tion, the stress-energy tensor (SET) component T0 2 briefly describes the general formalism, in Section is zero, therefore due to the Einstein equations 3 we find solutions of the field equations, in Section the Ricci tensor component R3 (ω e2γ+µ)′ = 0, 0 ∼ 4 we discuss the properties of Melvin-like solutions, which leads to [10] and in Section 5 the wormhole family. Section 6 −µ−2γ contains some concluding remarks. ω = ω0 e , ω0 = const. (6) Then, according to (5), 2 Basic relations 2γ(x) α+β−µ−3γ E(x) = 2ω0 e e dx. (7) We consider stationary cylindrically symmetric Z space-times with the metric − Note that Eqs. (4)–(7) are valid for an arbitrary ds2 = e2γ(x)[dt E(x) e 2γ(x) dϕ]2 e2α(x)dx2 − − choice of the radial coordinate x. Preserving this e2µ(x)dz2 e2β(x)dϕ2, (1) arbitrariness, we can write the nonzero components ν − − of the Ricci (Rµ) tensor as where x0 = t R, x1 = x, x2 = z R and ∈ ∈ 0 −2α ′′ ′ ′ ′ 2 x3 = ϕ [0, 2π) are the temporal, radial, longitu- R = e [γ + γ (σ α )] 2ω , ∈ 0 − − − dinal and angular (azimuthal) coordinates, respec- − ′′ ′ ′ ′ R1 = e 2α[σ + σ 2 2U α σ ] + 2ω2, tively. The variable x is here specified up to a sub- 1 − − − − ′′ ′ ′ ′ stitution x f(x), therefore its range depends on R2 = e 2α[µ + µ (σ α )], → 2 − − both the geometry itself and the “gauge” (the coor- − ′′ ′ ′ ′ R3 = e 2α[β + β (σ α )] + 2ω2, dinate condition). The off-diagonal component E 3 − − describes rotation, or the vortex component of the 0 0 −2γ 3 0 R3 = G3 = E e (R3 R0), (8) gravitational field. In the general case, this vortex − gravitational field is determined by the 4-curl of the where we are using the notations µ orthonormal tetrad field em (Greek and Latin let- ′ ′ ′ ′ ′ ′ ters are here assigned to world and tetrad indices, σ = β + γ + µ, U = β γ + β µ + γ µ . (9) respectively) [29,30]: The Einstein equations may be written in two 1 equivalent forms ωµ = εµνρσe ∂ em. (2) 2 mµ ρ σ Gν Rν 1 δν R = κT ν, or (10) Kinematically, the axial vector ωµ is the angu- µ ≡ µ − 2 µ − µ ν κ ν κ ν 1 ν lar velocity of tetrad rotation, it determines the Rµ = Tµ (Tµ δµT ). (11) − ≡− − 2 proper angular momentum density of the gravita- R being thee Ricci scalar and T the SET trace. We tional field, will mostly use the form (11), of the equations, but Sµ(g)= ωµ/κ, κ = 8πG, (3) it is also necessary to write the constraint equation from (10), which contains only first-order deriva- where G is the Newtonian gravitational constant. tives of the metric and represents a first integral of In space-times with the metric (1) we have the other equations: 1 − ′ − − − − ωµ = δµ2(E e 2γ ) eγ α β µ (4) G1 = e 2αU + ω2 = κT 1. (12) 2 1 − 1 (a prime stands for d/dx), and it appears sufficient Owing to the last line of (8) and its analogue for µ ν to consider its absolute value ω(x) = √ω ωµ that Tµ , in the Einstein equations it is sufficient to solve 3 the diagonal components, and then their only off- No-go theorem. It can be shown that a free lon- diagonal component holds automatically [10]. gitudinal magnetic field is incompatible with a non- As is evident from (8), the diagonal components static (ω = 0) metric (1). This follows from solving ν ν 6 µν of both the Ricci (Rµ ) and the Einstein (Gµ ) ten- the Maxwell equations, which, for F of the form sors split into the corresponding tensors for the (16) read static metric (the metric (1) with E = 0) plus a ′ ′ (√ gF 13) = 0, (√ gF 01) = 0, (17) contribution containing ω [10]: − − ν ν ν where Rµ = sRµ + ωRµ, ν 2 α+β+γ+µ ωRµ = ω diag( 2, 2, 0, 2), (13) g = det(g ), √ g = e . − µν − ν ν ν Gµ = sGµ + ωGµ, Equations (17) are integrated to give, respectively, ν 2 ωGµ = ω diag( 3, 1, 1, 1), (14) ′ −α−β+γ+µ ′ −α−β−γ+µ − − Φ e = h1, EΦ e = h2, (18) ν ν where sRµ and sGµ correspond to the static metric. It turns out that the tensors Gν and where h1, h2 = const = 0. From (18) it follows s µ −2γ 6 −2γ ′ ν E e = h2/h1 , whence (E e ) = 0, and ac- ωGµ (each separately) obey the conservation law α cording to (6), ω = 0. We have shown that a free αGµ = 0 in terms of this static metric.