arXiv:2007.01145v1 [gr-qc] 2 Jul 2020 .A Bronnikov, re A. K. general in wormholes and universes Melvin-like Rotating a ymtya oprdt h peia n is one spherical the to compared as symmetry cal a is cylindrically It [2,4]. static in details of more gravity. see set fields, of large Einstein-Maxwell symmetric source a from only case the special as mag- or field electric netic longitudinal a with sym- solution cylindrically static, metric regular geon” completely a magnetic is or [3], “electric an equations, and Maxwell 2] [1, in reviews see therein. others, references and electrovacuum, fluid vacuum, perfect gen- including of advent relativity, the since eral num- studied large and have obtained A solutions been symmetric jets. cylindrically static cosmic of observed ber the as strings of cosmic well as existence as structures possible extended linearly sym- the such cylindrically is studying configurations of metric motivations One the configurations. of field gravitational strong characterizing local extensions, its in and the solutions relativity exact general numerous (after obtain to possible mak- second it symmetry ing space-time the simplest one) is spherical symmetry Cylindrical Introduction 1 c d b a ainlRsac ula nvriy“Eh”(ocwEngi (Moscow “MEPhI” University Nuclear Research National epe’Finsi nvriyo usa(UNUniversity (RUDN Russia of University Friendship Peoples’ ocwSaeTcnlgclUiest Sakn,Vadkovs “Stankin”, University Technological State Moscow 1 etrfrGaiainadFnaetlMtooy VNIIMS, Metrology, Fundamental and Gravitation for Center -al [email protected] e-mail: nipratdsigihn etr fcylindri- of feature distinguishing important An Einstein- the to solution famous Melvin’s itiuin fapretfli ihteeuto fstate of equation f the equations with Einstein-Maxwell fluid the perfect to a solutions exact of of distributions family a find We xs hl oeohr i h case unive the magnetic (in static others Melvin’s some of while axis, that resembling geometries describe urn nteaglrdrcin hscretcetsamgei fi magnetic a creates current This direction. angular the in current bandpeiul,weetesuc fgaiyi tffmte with matter stiff is gravity of source the where previously, obtained reprmtr ftemdl hc hsfrsnweape fph of examples new limit forms Energ thus the Weak which In the model, relativity. jun satisfying the at matter of located of parameters shells, free thin consist The regions, space. (flat) mag of outer parts a fla curved potentia with weakly and configuration or matching wormhole rotation flat admit cylindrical with solutions a forming wormhole those thus The solution, throat, current. Melvin’s a Unlike require and axis. symmetry a a,b,c, 1 w .G Krechet, G. V. → ,temgei edtnst eo n h omoemdltnst tends model wormhole the and zero, to tends field magnetic the 1, > w )dsrb rvral omoegoere hc ontcontain not do which geometries wormhole traversable describe 0) d n .B Oshurko B. V. and orei efc udwith fluid perfect an such a of where is form solutions source the field find a a we in such and source current, that a electric out without turns exist It cannot longitudi- field. a to magnetic with nal solutions equations symmetric Einstein-Maxwell cylindrically the static the of and therein references [1,2]. also reviews see the etc., anisotropic [23–26] of kinds fluids some [18–22], notations) usual ouinta onthv ymtyai u con- wormhole but of axis symmetry family a a have not is do there curves that Also, timelike solution closed radii. of large emergence at the common their is particular, from feature in different ones, static quite the are of those solutions these of features ed[7,fliswt ieeteutoso state, of scalar equations with fluids a different perfect all, or with above [16] fluids charge [17], electric field dust with rotating dust differentially poten- [13–15], or rigidly self-interaction [10–12]; different tials with various [5–9]; fields with constant cosmological scalar known, the gravity: are of sources equations to of solutions Einstein number rotation) the great nonsymmet- (assuming a stationary general Accordingly, exact the sym- mention space-times. axial ric realistic to more not the metry, to compli- inherent avoiding rotation, cations of inclusion possible the ypr A ocw175,Russia 127055, Moscow 3A, per. ky ,u.Mkuh-aly ,Mso 118 Russia 117198, Moscow 6, Miklukho-Maklaya ul. ), p nti ae eoti oaigcounterparts rotating obtain we paper this In ern hsc nttt) ocw Russia Moscow, Institute), Physics neering zony l 6 ocw196,Russia 119361, Moscow 46, ul. Ozyornaya = wρ l ln the along eld odto ne rprcoc fthe of choice proper a under Condition y l iil o itn bevr eiigin residing observers distant for visible lly tosbtenteinr(omoe and (wormhole) inner the between ctions no-rewrhl oesi general in models wormhole antom-free h qaino state of equation the ( d | w s n oti eua symmetry regular a contain and rse -pc ein nbt ie fthe of sides both on regions t-space rrttn yidial symmetric cylindrically rotating or | < ) arigacrua electric circular a carrying 1), ei ldcno evacuum be cannot filed netic z xs oeo h solutions the of Some axis. p = p wρ = lativity p ρ , = . w h one the o wρ os (in const = Many . 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. There- ∇ ν κ longitudinal magnetic field cannot support a vortex fore, the tensor ωGµ/ may be interpreted as the SET of the vortex gravitational field. It possesses gravitational field with the metric (1). −2γ quite exotic properties (thus, the effective energy (Let us also note that in the case E e = E1 = −2γ density is 3ω2/κ < 0), which favor the existence const, the term E e dϕ is eliminated from (1) by − introducing the new time coordinate t′ = t E ϕ, of wormholes, and indeed, a number of wormhole − 1 solutions with the metric (1) have already been ob- making the metric explicitly static.) tained [10, 11, 28, 30] with sources in the form of scalar fields, isotropic or anisotropic fluids. Fur- 3.2 The fluid ther on we will obtain one more solution of this To circumvent the above no-go theorem, that is, kind, now supported by a perfect fluid and a mag- to avoid the relation E e−2γ = const, let us intro- netic field due to an electric current. duce a source of the magnetic field in the form of µ µ an electric current density J = ρeu , where ρe is 3 Solutions with a perfect fluid the effective electric charge density,2 and uµ is the 4-velocity satisfying the usual normalization condi- and a magnetic field µ tion uµu = 1. We will assume that the effective 3.1 The electromagnetic field. A no-go charge distribution is at rest in our rotating refer- theorem ence frame, so that µ −γ µ −γ Consider a longitudinal (z -directed) magnetic field, u = ( e , 0, 0, 0), J = (ρe e , 0, 0, 0). (19) corresponding to the 4-potential As usual, the electric charge conservation equa- µ tion µJ holds automatically due to the Maxwell Aµ = (0, 0, 0, Φ(x)), (15) ∇ equations F µν = J µ . ∇ν so that the only nonzero components of the Two nontrivial Maxwell equations now read Maxwell tensor Fµν = ∂µAν ∂ν Aµ are F13 = 13 ′ ′ − (√ gF ) = 0, (20) F31 = Φ (x). The nonzero contravariant compo- − − µν 1 ′ − nents F are (√ gF 01) = ρ e γ. (21) √ g − e − − ′ − − − ′ F 13 = e 2α 2βΦ , F 01 = E e 2α 2β 2γ Φ . (16) − 2If we introduce a real nonzero charge density, it becomes necessary to consider, in addition, a Coulomb electric field, The magnetic field magnitude (magnetic induction) which will make the problem hardly tractable. We therefore 2 13 B is determined by B = F F13 , so that the elec- consider an azimuthal electric current as if in a coil, in a µν 2 tromagnetic field invariant is Fµν F = 2B . neutral medium like a conductor with free electrons and ions at rest. 4

Integrating Eq. (20), we obtain, as before, where B2 = h2 e−2γ−2µ , and the only off-diagonal

′ − − component Φ = h eα+β γ µ, h = const, (22) 0 1 ′ −2γ T3 [e]= hΦ e (30) and substituting this Φ′ to Eq. (21) with (16), we −4π arrive at the following expression for ρe : does not affect the solution process, as mentioned

−2µ−3γ in a remark after Eq. (12). ρe = hω0 e . (23) Thus far all relations and expressions were writ- ten in terms of an arbitrary radial coordinate x. On the other hand, the electric charges should However, to solve the Einstein equations it is help- have a material carrier, for which we will assume a ful, at last, to choose the “gauge”, and by analogy perfect fluid with a barotropic equation of state and with our previous studies we will use the harmonic postulate a constant ratio of the effective charge radial coordinate corresponding to density ρe to energy density ρ: α = β + γ + µ. (31) ρe/ρ = A = const; p/ρ = w = const, (24) With the above expressions for the SET and p being the fluid pressure. We do not fix the value Eqs. (8), the noncoinciding components of Eqs. (11) of w but later on we will obtain a restriction on it. and (12), taking into account the expressions (6) for The fluid must obey the conservation law T ν = ∇ν µ ω and (26) for ρ, may be written as 0, which gives in our comoving reference frame ′′ 2 2β−2γ 1 + 3w 2β−γ ′ ′ γ = 2ω0 e + K e p + (p + ρ)γ = 0, (25) − 2 + Gh2 e2β, (32) which, for w = 0, leads to the expression ′′ w 1 − 6 µ = − K e2β γ + Gh2 e2β, (33) −γ(w+1)/w 2 ρ = ρ0 e , ρ0 = const. (26) ′′ − w 1 − β = 2ω2 e2β 2γ + − K e2β γ Gh2 e2β, (34) 0 2 − Comparing (23) and (26), taking into account the ′ ′ ′ ′ ′ ′ − β γ + β µ + γ µ = ω2 e2β 2γ assumption ρ /ρ = A = const, we obtain a relation 0 e − 2β−γ 2 2β between the metric coefficients e2γ and e2µ : + wK e + Gh e , (35)

− where K = κhω0/A (recall that κ = 8πG). Aρ e2µ = hω e γ(w+1)/(2w). (27) 0 0 Now, combining Eqs. (32) and (33) with (27), we arrive at the algebraic equation for γ with con- In the case w = 0 (zero pressure), the conservation stant coefficients, law (25) simply leads to γ = const. − − 4ω2(1 2w) e 2γ + (8w2 3w 1)K e γ 0 − − − 3.3 Solution of the Einstein equations + 2(4w 1)Gh2 = 0, (36) − To address the Einstein equations, let us write the from which it follows γ = const, and we can with- expressions for the SETs of the perfect fluid and out loss of generality put γ 0 by choosing a ≡ the electromagnetic field. For the fluid we have time scale. Then Eq. (27) implies µ = const which also allows us to put µ 0 by choosing the scale T ν[f]= ρ diag(1, w, w, w). (28) ≡ µ − − − along z ; from (27) we then obtain a relation among the constants: hω = Aρ , so that, in particular, For the electromagnetic field SET we have the stan- 0 0 K = κρ . dard expression 0 With eγ = eµ = 1, from Eqs. (32) and (33) we 1 find T ν[e]= F F να + 4δν F F αβ , µ 16π − µα µ αβ
2 1 2 1 Gh = (1 w)K, ω0 = (1 + w)K, (37) which in our case leads to 2 − 2 which leads to a conclusion on the range of w: B2 T ν[e]= diag(1. 1, 1, 1) T 0[e], (29) µ 8π − − ⊕ 3 1

With (37) it is directly verified that Eqs. (35) and In all cases the off-diagonal metric function E is (36) also hold. All our constant parameters may be easily obtained as expressed in terms of two of them, for example, ω0 2β and w: E(x) = 2ω0 e dx. (47) Z 1 w 2ω2 Gh2 = − ω2, κρ = 0 , 1+ w 0 0 1+ w 4 Melvin-like universes ρ A = e = 4π G(1 w2). (39) ρ p − Melvin’s electric or magnetic geon [3] is among the We see that in our system not only µ = γ = 0, most well-known static, cylindrically symmetric so- lutions to the Einstein-Maxwell equations; it is a but also both densities ρ and ρe as well as the angular velocity ω are constant. It is also of inter- special solution from a large class of static, cylin- est that the two limiting cases of the equation of drically symmetric solutions with radial, azimuthal state, w = 1 (maximum stiffness compatible with and longitudinal electric and/or magnetic fields, causality) and w = 1 (the cosmological constant) see, e.g., [1, 2, 4]. Its metric may be written in the − are excluded in the present system. In both these form [2,4] cases, static and stationary cylindrically symmetric ds2 = (1+ q2x2)2(dt2 dx2 dz2) − − solutions without an electromagnetic field are well x2 dϕ2, (48) known [2,6,7,10,12,26,28]. − (1 + q2x2)2 The remaining differential equation (34) has the where x 0, and the magnetic (let us take it for Liouville form ≥ certainty) field magnitude is ′′ 4w 2 2β β = ω e . (40) 2 2 −2 w + 1 0 B = Bz = 2q(1 + q x ) , and has the first integral with q = const characterizing the effective current that might be its source. However, this solution de- ′2 4w 2 2β 2 scribes a purely field configuration existing without β = ω0 e + k sign k, (41) w + 1 any massive matter, electric charges or currents. with k = const. The further integration depends Both the metric and the magnetic field are regu- lar on the axis x = 0. The other “end”, x , on the signs of w and k: → ∞ is infinitely far away (the distance √ gxxdx di- k − 1. w< 0, k> 0, eβ = ; . (42) verges), the magnetic field vanishesR there, and the m cosh(kx) circular radius r = √ g also tends to zero, so − ϕϕ k that the whole configuration is closed in nature, 2. w> 0, k> 0 : eβ = ; . (43) m sinh(kx) without spatial infinity, and with finite total mag- netic field energy per unit length along the z axis. 1 3. w> 0, k =0: eβ = ; . (44) As we saw in Sec. 3.1, such a free magnetic field mx cannot support a rotating counterpart of Melvin’s k solution, but Einstein-Maxwell solutions with a lon- 4. w> 0, k< 0 : eβ = | | ; . (45) m cos( k x) gitudinal magnetic field are obtained in the pres- | | ence of perfect fluids with electric currents. Let us 4 w ω2 1/2 briefly discuss their main features. where we have denoted m = | | 0 .  w + 1  In all cases under consideration, the magnetic In the previously excluded case w = 0 (dust- field is directed along the z axis and has the con- like matter), the equality γ = const is immediately stant magnitude B = h, while the metric has the obtained from (25), µ = const then follows from form (27), and as before, without loss of generality, we ds2 = (dt Edϕ)2 e2βdx2 dz2 e2βdϕ2, (49) can put µ = γ = 0. Instead of (40), we obtain − − − − ′′ β = 0 whence we can write and E is determined by Eq. (53). Note that both ϕ and x are dimensionless while t, z and eβ have eβ = r ekx, r , k = const. (46) 0 0 the dimension of length. 6

Dustlike matter, Eq. (46) of 3-dimensional space. However, as follows from (55), g is positive (hence contains CTCs) where Let us begin with the case w = 0. For E(x) we 33 tanh kx > w , that is, at large enough x , in find | | | | | | circular regions around the two poles. 2 ω0r0 2kx ω0 2 By choosing another value of the integration E = E0 + e = E0 + r , k k constant E0 one can make one of the poles free kx r = r0 e , E0 = const, (50) from CTCs, at the expense of enlarging the CTC region around the other pole. One of the poles can where E0 is an integration constant. In terms of even be made regular by a proper choice of the pa- the coordinate r, the metric reads 2 rameters. For example, choosing E0 = 2kω0/m , 2 we obtain E = 0 at x = , and it is easy to 2 2 −2 2 2 2 2 −∞ ds = dt Edϕ k dr dz r dϕ . (51) verify that the pole x = is then regular under − − − − −∞   the condition k = 1. The symmetry axis r = 0 is regular in the case 3 E0 = 0, k = 1. Also, in this case Solution 2, Eq. (43) g = E2 r2 = r2( 1+ ω2r2) (52) For w> 0, k> 0, with (43), for E(x) we find 33 − − 0 2kω changes its sign at r = ω−1 , and at larger r the E = E 0 coth kx. (56) 0 − m2 lines of constant t, r, z (coordinate circles) are time- like, thus being closed timelike curves (CTCs) vio- The metric takes the form lating causality. ds2 = (dt Edϕ)2 dz2 − − k2 (dx2 + dϕ2), (57) Solution 1, Eq. (42) − m2 sinh2(kx) For w< 0, with (42), for E(x) we calculate It is convenient to introduce the new coordinate y by substituting 2kω0 E = E + tanh kx, (53) −2kx 2k 0 2 e = 1 , (58) m − y The metric has the form after which we obtain y 2ω ds2 = (dt Edϕ)2 dz2 r2 = e2β = (y 2k). E = E 0 y. (59) − − m2 − 0− m2 k2 (dx2 + dϕ2), (54) The range x > 0 is converted to y 2k, where 2 2 ≥ − m cosh (kx) y = 2k is the axis of symmetry. The metric now where, for convenience, we have rearranged the has the form terms with dz2 and dx2 as compared to (49). dy2 ds2 = (dt Edϕ)2 For g33 , similarly to (52), again putting E0 = 0 − − my(y 2k) y − and recalling the definition of m, we obtain dz2 (y 2k)dϕ2. (60) − − m2 − 2 k 1 2 Assuming E0 = 0, for g33 it is then easy to obtain g33 = 1 tanh kx . (55) −m2  − w  the expression | | y y In this solution x R, and at both extremes g = + 2k > 0, (61) ∈ 33 m2 w x we have r = eβ 0, i.e., these are   → ±∞ → two centers of symmetry (or poles) on the (x, ϕ) 2- which means that CTCs are present everywhere, surface, or two symmetry axes from the viewpoint and actually this space-time has an incorrect sig- nature, (+ +) instead of (+ ). 3 −− − −− The axis regularity conditions require [1,2,34] finite val- However, with nonzero values of E0 it becomes ues of the curvature invariants plus local flatness (sometimes possible to get rid of CTCs in some part of space. also called “elementary flatness”) as a correct circumfer- ence to radius ratio for small circles around the axis, which Thus, choosing E0 in such a way that E = 0 at −2α ′2 in our case leads to the condition e R 1, where some y0 > 2k, we will obtain the normal sign g33 < → R = √ g33 . 0 in some range of y around y . − 0 7

Solution 3, Eq. (44) obtained in [28] for a cylindrical wormhole with stiff matter. In the case w> 0, k = 0, with (44), it is convenient As always with rotating cylindrical wormhole to use the coordinate r = 1/(mx), and then we solutions, these wormholes do not have a flat-space obtain asymptotic behavior at large x , which makes it | | 2ω0 impossible to interpret them as objects that can E = E0 + , r be observed from regions with small curvature. To dr2 ds2 = (dt Edϕ)2 dz2 r2dϕ2, (62) overcome this problem, it has been suggested [10] − − m2r2 − − to cut out of the obtained wormhole solution a reg- ular region, containing a throat, and to place it and assuming E0 = 0, we arrive at between two flat regions, thus making the whole 2 2 2ω0 2 1 w system manifestly asymptotically flat. However, to g33 = r 1 2 = r − > 0. (63) −  − m  2w interpret such a “sandwich” as a single space-time, it is necessary to identify the internal and external We again obtain a configuration with an incorrect metrics on the junction surfaces Σ and Σ− , which signature, possessing CTCs at all r. However, as + should be common for these regions. The internal in the previous case, by choosing E so that E = 0 0 region will be described in the present case by (64), at some r = r we can provide a CTC-free region 0 (65)). And since the internal metric contains rota- in a thick layer around r = r . 0 tion, the external Minkowski metric should also be taken in a rotating reference frame. 5 Wormholes Thus we take the Minkowski metric in cylindri- cal coordinates, ds2 = dt2 dX2 dz2 X2dϕ2 , β M − − − With the solution (45) for r = e , the range of and convert it to a rotating reference frame with x is x π/(2k¯), π/(2k¯) , where k¯ = k > 0, ∈ − − angular velocity Ω = const by substituting ϕ and we see that r on
both ends, confirming → → ∞ ϕ +Ωt, so that the wormhole nature of this configuration, where x = 0 is the wormhole throat (minimum of r). ds2 = dt2 dX2 dz2 X2(dϕ +Ωdt)2. (67) M − − − Substituting y = k¯ tan kx¯ , we obtain the metric in the form In the notations of (1), the relevant quantities in (67) are 2 2 2 dy ds = (dt Edϕ) X2 − − m2(k¯2 + y2) e2γ = 1 Ω2X2, e2β = , − 1 Ω2X2 k¯2 + y2 − dz2 dϕ2, (64) Ω − − m2 E =ΩX2, ω = . (68) 1 Ω2X2 where y R, and y = 0 is the throat; furthermore, − ∈ This stationary metric admits matching with the 2ω internal metric at any X < 1/ Ω , inside the “light E = E + 0 y, (65) | | | | 0 2 cylinder” X = 1/ Ω on which the linear rota- m | | | | tional velocity coincides with the speed of light. and for g33 in the case E0 = 0 (which makes the Making use of the symmetry of (64), let us as- solution symmetric with respect to y = 0) it follows sume that the internal region is y∗ Minkowski space, respectively, so that the external The expression (66) shows that CTCs are absent flat regions are X < X∗ and X > X . Matching − + around the throat, at y2 < 2k¯2w/(1 w), while at is achieved if we identify there the two metrics, so − larger y the CTCs emerge. that | | Let us note that in the limit w 1, so that the → fluid EoS tends to that of maximally stiff matter, [β] = 0, [µ] = 0, [γ] = 0, [E] = 0, (69) the magnetic field disappears (h 0 according → where, as usual, the brackets [f] denote a disconti- to (37)), and the whole solution tends to the one nuity of any function f across the surface. Under 8 the conditions (69), we can suppose that the co- condition [γ] = 0 we will rescale the time coordi- ordinates t,z,φ are the same in the whole space. nate in the internal region according to At the same time, there is no need to adjust the t = √P t,˜ P := 1 Ω2X2 (76) choice of radial coordinates on different sides of the − junction surfaces since the quantities involved in and use the new coordinate t˜, with which, instead all matching conditions are insensitive to possible of E , we must use E√P in all formulas. The re- reparametrizations of y or X . maining two conditions (69) yield Having identified the metrics, we certainly did 2ω k2 + y2 X2 not adjust their normal derivatives, whose jumps 2 0 √ ΩX = y P , 2 = , (77) are well known to determine the properties of mat- m m P ter filling a junction surface Σ and forming there a where, without risk of confusion, we omit the as- b thin shell. The SET Sa of such a thin shell is cal- terisk at X and y . With these conditions, there culated using the Darmois-Israel formalism [31–33], are four independent parameters of the system, for b and in the present case of a timelike surface, Sa is example, we can choose as such parameters b related to the extrinsic curvature Ka of Σ as 2 2ω0 w + 1 − X, y, P, n = = . (78) Sb = (8πG) 1[K˜ b ], K˜ b := Kb δb Kc, (70) m 2w a − a a a − a c The other parameters are expressed in their terms where the indices a, b, c = 0, 2, 3. The general ex- as ˜ b pressions for nonzero components of Ka for sur- √1 P ny√P faces x = const in the metric (1) are [28] Ω= − , ω0 = , X X√1 P ˜ −α+2γ ′ ′ − K00 = e (β + µ ), y2(2n 1+ P ) − k2 = − . (79) ˜ 1 −α ′ −α ′ ′ ′ 1 P K03 = 2 e E + E e (β + γ + µ ), − − Now we can calculate the quantities (75), with ˜ −α+2µ ′ ′ K22 = e (β + γ ), [f]= f f on Σ : out − in + ˜ −α+2β ′ ′ K33 = e (γ + µ ) yP 3/2 1 1 y√P − − ′ ′ ′ ′ a = − , b = − , + e α 2γ [EE E2(β + 2γ + µ )]. (71) P X X − 1 P nyP 3/2 1+ P From (71) it is straightforward to find Sa on c = − , d = − . (80) b − P X P X√1 P the surfaces Σ± . However, our interest is not in − finding these quantities themselves but, instead, a It can be easily verified that the conditions (72), a (73), (74) are satisfied as long as verification of whether or not the resulting SET Sb satisfies the WEC. Let us use for this purpose the 2 P y − , (81) necessary and sufficient conditions obtained in a ≥ P 3/2 general form in [27], see also a detailed description in full analogy with the corresponding calculation in [28]. These conditions are in [28]. a + c + (a c)2 + 4d2 0, (72) We have shown that under the condition (81) − ≥ p the WEC holds on Σ+ . Now, what changes on the 2 2 a + c + (a c) + 4d + 2b 0, (73) surface Σ− specified by X = X∗ < 0 and y = p − ≥ − y∗ < 0, where we must take [f] = f f for a + c 0, (74) − out − in ≥ any function f ? As in [28], it can be verified that where the parameters a, b, c do not change from (80) if we − ′ ′ − ′ ′ replace X with X (we denote, as before, y = y∗ > a = [ e α(β + µ )], b = [ e α(β + γ )], | | − 0). For d = [ω] there will be another expression −α ′ ′ − c = [ e (γ + µ )], d = [ω] (75) since, according to (69), Ω(Σ−) = Ω(Σ ), while − − + in the internal solution ω(Σ−)= ω(Σ ), hence on Let us discuss, for certainty, the conditions on + Σ− Σ+ : y = y∗, X = X∗ with our metrics (64), (65) and (67). Among the matching conditions (69), n y P 3/2 + 1 P d d− = | | − , [µ] = 0 holds automatically, while to satisfy the 7→ − X P √1 P | | − 9

so that d− > d , making it even easier to sat- gravitational fields with cylindrical symmetry, and | | | | isfy the WEC requirements. As a result, the WEC their asymptotic behavior making them potentially holds under the same condition (81), providing a observable from flat or weakly curved regions of wormhole model which is completely phantom-free. space is provided by joining flat regions on both We can also notice that from (79) it follows sides of the throat. Such a complex structure is nec- 2 2 2 2 y∗ < k , therefore, y < k in the whole internal essary because asymptotic flatness at large circular region, which is thus free from CTCs. radii cannot be achieved in any cylindrical solutions There is one more point to bear in mind: since with rotation. The present family of models with there is a z -directed magnetic magnetic field in a magnetic field, parametrized by the equation-of- the internal region, we must suppose that there state parameter w < 1, tends to the one obtained are some surface currents on Σ± in the ϕ direc- in [28] in the limit w 1, in which the magnetic → tion. Their values can be easily calculated using field vanishes. µν µ the Maxwell equations νF = J . Indeed, say, Let us mention that other static or stationary ∇ µν Σ+ (x = x∗) separates the region where F = 0 wormhole models with proper asymptotic behavior from the one with nonzero F µν , therefore at their and matter sources respecting the WEC have been junction we have J µ = δ(x x∗)J µ(x∗ 0), so obtained in extensions of general relativity, such − − − that the surface current is J a = J µ(x∗ 0) . as the Einstein-Cartan theory [35, 36], Einstein- − − µ=a Similarly, on Σ− (x = x∗) we obtain J a = Gauss-Bonnet gravity [37], multidimensional grav- − J µ( x∗ + 0) . In our wormhole configurations ity including brane worlds [38, 39], theories with − µ=a nonmetricity [41], Horndeski theories [40], etc. we obtain, according to (24), (23), (27) and taking into account that γ = µ 0, One can also notice that the same trick as was ≡ used with wormhole models, that is, joining a flat a 0 0 J = (J , 0, 0), J (Σ±)= hω0. (82) region taken in a rotating reference frame, can be ∓ used as well with solutions possessing a symmetry Thus the surface currents on Σ± have only the tem- axis. It is important that in all such cases the sur- poral component, i.e., they are comoving to the face to be used as a junction should not contain matter and current in the internal region. CTCs (in other words, there should be, as usual, As is the case with the internal wormhole solu- g < 0) because g < 0 in the admissible part of tion, in the limit w 1 (hence n 1) the whole 33 33 → → flat space, while g33 taken from the external and in- twice asymptotically flat construction tends to the ternal regions should be identified at the junction. one obtained in [28] with a stiff matter source. In this way one can obtain completely CTC-free models of extended cosmic strings with rotation. 6 Concluding remarks A possible observer can certainly be located far from such an extended string or wormhole config- We have obtained a family of stationary cylindri- uration and be at rest in a nonrotating frame in cally symmetric solutions to the Einstein-Maxwell flat space, other than the one used for the object equations in the presence of perfect fluids with construction. p = wρ, w < 1. Some of them (Solutions 1–3) | | An evident further development can be a search contain a symmetry axis which can be made reg- for other rotating configurations with electromag- ular by properly choosing the solution parameters. netic fields, possibly including radiation in differ- The only geometry of closed type belongs to Solu- ent directions in the spirit of [42], where radial, az- tion 1, Eqs. (42), (53), (54), (55). Unlike Melvin’s imuthal and longitudinal radiation flows were con- solution and like all other solutions with rotation, it sidered as sources of gravity in space-times with the inevitably contains a region where g33 > 0, so that metric (1). Another set of problems concerns elec- the coordinate circles parametrized by the angle ϕ trostatics in the fields of extended strings or worm- are timelike, violating causality. holes with sources including electromagnetism. As The wormhole models discussed here are of in- follows from [43], even in simpler, partly conical terest as new examples of phantom-free wormholes cylindrical geometries with thin shells electrostat- in general relativity, respecting the WEC. As in ics turns out to be rather interesting and complex. other known examples [27, 28], such a result is achieved owing to the exotic properties of vortex 10

Acknowledgments [11] K.A. Bronnikov, V.G. Krechet, “Rotating cylindri- cal wormholes and energy conditions.” Int. J. Mod. This publication was supported by the RUDN Uni- Phys. A 31, 1641022 (2016); arXiv: 1509.04665. versity program 5-100 and by the Russian Foun- [12] C. Erices and C. Martinez, “Stationary cylindri- dation for Basic Research Project 19-02-00346. cally symmetric spacetimes with a massless scalar The work of K.B. was also performed within the field and a nonpositive cosmological constant.” framework of the Center FRPP supported by Phys. Rev. D 92, 044051 (2015). MEPhI Academic Excellence Project (contract No. [13] W. J. van Stockum, “The gravitational field of a 02.a03.21.0005, 27.08.2013). The work of V.K. and distribution of particles rotating about an axis of V.O. was supported by the Ministry of Education symmetry.” Proc. Roy. Soc. Edinburgh A 57, 135 and Science of Russia in the framework of State (1937). Contract 9.1195.2017.6/7. [14] W. B. Bonnor and B. R. Steadman, “A vacuum exterior to Maitras cylindrical dust solution.” Gen. References Rel. Grav. 41, 1381 (2009). [15] B. V. Ivanov, “The general double-dust solution,” [1] H. Stephani, D. Kramer, M. A. H. MacCallum, gr-qc/0209032 C. Hoenselaers, and E. Herlt, Exact solutions of Einsteins field equations, Cambridge Monographs [16] B.V. Ivanov, “Rigidly rotating cylinders of charged on Mathematical Physics (Cambridge University dust.” Class. Quantum Grav. 19, 5131 (2002); gr- Press, 2009). qc/0207013. [2] K.A. Bronnikov, N.O. Santos, and Anzhong Wang, [17] N. O. Santos and R. P. Mondaini, “Rigidly rotating “Cylindrical systems in general relativity,” arXiv: relativistic generalized dust cylinder.” Nuovo Cim. 1901.06561. B. 72, 13 (1982) [3] M. A. Melvin, “Pure magnetic and electric geons,” [18] C. Hoenselaers and C. V. Vishveshwara, A rela- Phys. Lett. 8, 65-68 (1964). tivistically rotating fluid cylinder. Gen. Rel. Grav. 10, 43-=51 (1979). [4] K. A. Bronnikov, “Static, cylindrically symmetric Einstein-Maxwell fields,” in Problems in Gravita- [19] W. Davidson, “Barotropic perfect fluid in steady tion Theory and Particle Theory (PGTPT), Ed. K. cylindrically symmetric rotation.” Class. Quantum P. Staniukovich, 10th issue, p. 37-50 (Atomizdat, Grav. 14, 119 (1997). Moscow, 1979, in Russian). [20] W. Davidson, “A cylindrically symmetric station- [5] C. Lanczos, “´’Ueber eine station¨are Kosmologie in ary solution of Einsteins equations describing a per- Sinne der Einsteinischen Gravitationstheories.” Z. fect fluid of finite radius.” Class. Quantum Grav. Physik 21, 73 (1924). 17, 2499 (2000). [6] T. Lewis, “Some special solutions of the equations [21] D. Sklavenites, “Stationary perfect fluid cylinders.” of axially symmetric gravitational fields.” Proc. R. Class. Quantum Grav. 16, 2753 (1999). 136 Soc.. A , 176 (1932). [22] B.V. Ivanov, “On rigidly rotating perfect fluid [7] N. O. Santos, “Solution of the vacuum Einstein cylinders.” Class. Quantum Grav. 19, 3851 (2002); equations with nonzero cosmological constant for gr-qc/0205023. a stationary cylindrically symmetric spacetime.” [23] P. S. Letelier and E. Verdaguer. “Anisotropic fluid 10 Class. Quantum Grav. , 2401 (1993). with SU(2) type structure in general relativity: A [8] A. Krasinski, “Solutions of the Einstein field equa- model of localized matter.” J. Math. Phys. 28, 2431 tions for a rotating perfect fluid II: properties of the (1987). flow-stationary and and vortex-homogeneous solu- [24] L. Herrera, G. Le Denmat, G. Marcilhacy, and N. 6 tions.” Acta Phys. Pol. , 223 (1975). O. Santos. “Static cylindrical symmetry and con- [9] M. A. H. MacCallum and N. O. Santos. “Stationary formal flatness.” Int. J. Mod. Phys. D 14 (2005) and static cylindrically symmetric Einstein spaces 657; gr-qc/0411107. 15 of the Lewis form.” Class. Quantum Grav. , 1627 [25] F. Debbasch, L. Herrera, P.R.C.T. Pereira, N.O. (1998). Santos, “Stationary cylindrical anisotropic fluid.” [10] K. A. Bronnikov, V. G. Krechet, and Jos´eP. S. Gen. Rel. Grav. 38, 1825 (2006); gr-qc/0609068. Lemos, “Rotating cylindrical wormholes.” Phys. Rev. D 87, 084060 (2013); arXiv: 1303.2993. 11

[26] S. V. Bolokhov, K. A. Bronnikov, and M. V. [41] V. G. Krechet, V. B. Oshurko, and M. N. Lodi, In- Skvortsova, “Rotating cylinders with anisotropic duced nonlinearities of the scalar field and worm- fluids in general relativity,” Grav. Cosmol. 25, 122 holes in the metric-affine theory of gravity, Grav. (2019); arXiv: 1904.06727. Cosmol. 24, 186–190 (2018). [27] K. A. Bronnikov and, V. G. Krechet, “Potentially [42] K. A. Bronnikov, String clouds and radiation flows observable cylindrical wormholes without exotic as sources of gravity in static or rotating cylinders, matter in GR.” Phys. Rev. D 99, 084051 (2019); Int. J. Mod. Phys. A 35, 2040004 (2020); arXiv: arXiv: 1807.03641. 1909.09409. [28] K. A. Bronnikov, S. V. Bolokhov, and M. V. [43] E. Rub´ın de Celis and C. Simeone, Electrostat- Skvortsova, “Cylindrical wormholes: a search for ics and self-force in asymptotically flat cylindrical viable phantom-free models in GR.” Int. J. Mod. wormholes, Eur. Phys. J. C 80, 501 (2020); arXiv: Phys. D 28, 1941008 (2019); arXiv: 1903.09862. 2004.08506. [29] V. G. Krechet, Topological and physical effects of rotation and spin in the general relativistic the- ory of gravitation, Izvestiya Vuzov, Fiz. No 10, 57 (2007); Russ. Phys. J. 50, 1021 (2007). [30] V. G. Krechet and D. V. Sadovnikov, Spin-spin in- teraction in general relativity and induced geome- tries with nontrivial topology. Grav. Cosmol. 15, 337 (2009); arXiv: 0912.2181. [31] G. Darmois, Les ´equations de la gravitation ein- steinienne, in: M´emorial des Sciences Mathema- tiques (Gauthier-Villars, Paris, 1927), vol. 25. [32] W. Israel, Singular hypersurfaces and thin shells in general relativity. Nuovo Cim. B 48, 463 (1967). [33] V. A. Berezin, V. A. Kuzmin, and I. I. Tkachev, Dynamics of bubbles in general relativity. Phys. Rev. D 36, 2919 (1987). [34] K. A. Bronnikov and S. G. Rubin, Black Holes, Cosmology and Extra Dimensions (World Scien- tific, 2013). [35] K.A. Bronnikov, A.M. Galiakhmetov, Wormholes without exotic matter in Einstein-Cartan theory, Grav. Cosmol. 21, 283 (2015); arXiv: 1508.01114. [36] K.A. Bronnikov, A.M. Galiakhmetov, Wormholes without exotic matter in Einstein-Cartan the- ory, Phys. Rev. D 94, 124006 (2016); arXiv: 1607.07791. [37] H. Maeda, M. Nozawa, Static and symmet- ric wormholes respecting energy conditions in Einstein-Gauss-Bonnet gravity. Phys. Rev. D 78, 024005 (2008). [38] K.A. Bronnikov, S.-W. Kim, Possible wormholes in a brane world. Phys. Rev. D 67, 064027 (2003); gr-qc/0212112. [39] K.A. Bronnikov, M.V. Skvortsova, Wormholes leading to extra dimensions Grav. Cosmol. 22, 316 (2016). [40] Sergey V. Sushkov, R V. Korolev, Scalar wormholes with nonminimal derivative coupling, Class. Quan- tum Grav. 29, 085008 (2012); arXiv: 1111.3415.