Cylindrical Systems in General Relativity

Kirill A. Bronnikov1,2,3,∗ N. O. Santos 4,5,† and Anzhong Wang 6,7‡ 1 Center for Gravitation and Fundamental Metrology, VNIIMS, Ozyornaya ul. 46, Moscow 119361, Russia 2 Institute of Gravitation and Cosmology, Peoples’ Friendship University of Russia (RUDN University), ul. Miklukho-Maklaya 6, Moscow 117198, Russia 3 National Research Nuclear University “MEPhI” (Moscow Engineering Physics Institute), Kashirskoe sh. 31, Moscow 115409, Russia 4 Sorbonne Universit´e,UPMC Universit´eParis 06, LERMA, UMRS8112 CNRS, Observatoire de Paris-Meudon, 5, Place Jules Janssen, F-92195 Meudon Cedex, France 5 School of Mathematical Sciences, Queen Mary, University of London, Mile End Road, London E1 4NS, UK 6 GCAP-CASPER, Physics Department, Baylor University, Waco, TX 76798-7316, USA and 7 Institute for Theoretical Physics & Cosmology, Zhejiang University of Technology, Hangzhou, 310032, China (Dated: May 13, 2020) With the arrival of the era of gravitational wave astronomy, the strong gravitational field regime will be explored soon in various aspects. In this article, we provide a general review over cylindrical systems in Einstein’s theory of general relativity. In particular, we first review the general prop- erties, both local and global, of several important solutions of Einstein’s field equations, including the Levi-Civita and Lewis solutions and their extensions to include the cosmological constant and matter fields, and pay particular attention to properties that represent the generic features of the theory, such as the formation of the observed extragalactic jets and gravitational Faraday rota- tion. We also review studies of cylindrical wormholes, gravitational collapse and Hoop conjecture, and polarizations of gravitational waves. In addition, by rigorously defining cylindrically symmetric spacetimes, we clarify various (incorrect) claims existing in the literature, regarding to the generality of such spacetimes.

CONTENTS IV. LC solutions coupled with cosmological constant 14 I. Introduction2 A. Linet-Tian (LT) solutions 14 A. Cylindrically symmetric spacetimes3 B. Main properties of LT solutions with Λ > 0 14 B. Einstein’s field equations5 1. Global structure of the solution with 1 σ = 2 15 II. Levi-Civita (LC) vacuum spacetime6 2. Global structure of the solution with 1 A. LC solutions6 σ = − 2 15 B. Nature of the coordinates of the LC C. Main properties of LT solutions with Λ < 0 16 solutions7 1. Global structure of the solution with C. Geodesics8 σ = 0 16 D. LC spacetime as a limiting case of the γ 2. Global structure of the solution with 1 spacetime8 σ = 2 16 E. Sources producing LC spacetime9 3. Global structure of the solution with 1 arXiv:1901.06561v4 [gr-qc] 12 May 2020 σ = − 2 17 III. Static electrovacuum spacetimes 10 D. Geodesics of LT solutions 18 A. Electromagnetic fields in cylindrical 1. Planar geodesicsz ˙ = 0 18 spacetimes 10 2. Non-planar geodesicsz ˙ 6= 0 18 B. Geometry with a radial electric field 11 E. Instability of LT solutions 19 C. Geometry with a longitudinal magnetic F. Sources producing LT spacetimes 19 field 12 D. Geometry with an azimuthal magnetic V. Static cylindrical spacetimes with perfect field 13 fluids 20 A. General consideration 20 B. Solutions for p = wρ, the general case 21 1. Solutions with a regular axis 22 ∗ [email protected] C. Solutions for particular values of w = p/ρ 23 † [email protected] 1. Stiff matter, w = 1 23 ‡ anzhong [email protected]; Corresponding author 2. Disordered radiation, w = 1/3 23 2

3. A gas of cosmic strings, w = −1/3 24 later rediscovered by Einstein and Rosen in their seminal 4. A cosmological constant, w = −1 24 work on the studies of the nonlinearity of GWs in 1937 D. Some further results 24 [137], and in the same year van Stockum solved the prob- lem of a rigidly rotating infinitely long cylinder of dust, VI. Lewis vacuum spacetime 25 and found explicitly the corresponding metric [308]. In A. Stationary cylindrical vacuum spacetime 25 1957, Bonnor [43] and Weber and Wheeler [327] studied B. The parameters of the Weyl class 26 the Einstein-Rosen waves in great detials, and since then, C. Geodesics for the Weyl class 27 cylindrically symmetric spacetimes have been extensively D. The parameters of the Lewis class 27 investigated with various motivations [167, 294]. E. Sources producing Lewis spacetime 28 In particular, to understand some fundamental issues F. Rotation and translation of cylinders 29 in general relativity (GR), such as the structure of the theory [241], the nature and formation of spacetime sin- VII. Lewis vacuum spacetime coupled with a gularities [177], the cosmic censorship [257] and hoop cosmological constant 29 [302] conjectures, one often assumes certain symmetries A. Stationary cylindrical vacuum spacetime of spacetimes, because Einstien’s field equations are in with cosmological constant 29 general so complicated that it is extremely difficult to B. A source producing stationary cylindrical study them in their most general form, especially in the vacuum spacetime with cosmological times when computers had not been available. In fact, constant 30 even now it is still very difficult to study spacetimes with C. Geodesic motion, confinement and only one Killing vector (both analytically and numeri- extragalactic jet formation 31 cally). Therefore, the next step is naturally to consider spacetimes with two Killing vectors, which include cylin- VIII. Cylindrical wormholes 32 drically symmetric spacetimes [167, 294]. A. Wormhole definitions and existence The studies of cylindrically symmetric spacetimes were conditions 33 further motivated by topological defects discovered in B. Examples of static wormholes 34 1976 by Kibble [202, 203], which include cosmic strings, C. Examples of rotating wormholes 34 and the latter were once believed to provide a mechanism 1. Weak energy condition on a junction to produce the cosmic microwave background (CMB) and surface 34 large-scale structure observed in our universe [108, 309]. 2. Minkowski regions 35 However, observations of CMB have ruled out cosmic 3. A model with stiff matter, p = ρ 35 strings, formed in the context of symmetry breaking in 4. A model with an anisotropic fluid 36 grand unified theories (GUT’s), as the sources of the cosmological perturbations [26, 289, 291], and led to an IX. Gravitational collapse and emission of 2 −7 upper bound, Gµ/c . 10 [7], where µ denotes the gravitational waves 37 string’s tension, and G and c are, respectively, the New- A. The Hoop conjecture and critical collapse 37 tonian constant and the speed of light in vacuum. The B. Polarizations and Faraday rotation 38 string tension is intimately related to the energy E of the phase transition, X. Conclusions 41 ! Gµ E2 Acknowledgements 42 2 = O 2 , (1.1) c Mpl Appendix A: Static cylindrical spacetimes 42  p  where Mpl ≡ ~c/G denotes the Planck mass. Cur- References 43 2 −11 rently, the strongest constraint Gµ/c . 10 comes from the analysis of stochastic GW background, expected from cosmic strings with the latest pulsar timing arrays I. INTRODUCTION [38]. Recently, the subject has attracted further attention The investigations of cylindrically symmetric space- in the framework of string/M-Theory, the so-called cos- times can be traced back as far as to 1919 when Levi- mic superstrings [109, 132], which were formed before Civita (LC) discovered a class of solutions of Einstein’s inflation took place and stretched to macroscopic length vacuum field equations, corresponding to static cylindri- scales in the inflationary phase. During the subsequent cal spacetimes [218]. The extension of the LC space- epochs, a complicated network of various string elements times to stationary ones was obtained independently by forms [99, 100, 306]. The main phenomenological con- Lanczos in 1924 [213] and Lewis in 1932 [219]. In 1925, sequence of a string network is the emission of GWs Beck studied a class of exact solutions and interpreted [1,3, 36, 55, 115–117, 252, 267, 269, 274, 295, 307]. It them as representing the propagation of cylindrical grav- can generate bursts at cusps, kinks and junctions, as well itational waves (GWs) [24]. This class of solutions was as a stochastic GW background. Low tension strings are 3 natural in string/M-Theory, and can easily satisfy the GWs have been also studied [240]. observational bounds given above [98]. Along the above direction, another interesting phe- With the observations of the 11 GWs [4], lately the nomenon related to the nonlinear properties of GWs is studies of GW phenomena have attracted a great deal of the gravitational Faraday rotation, first discovered by Pi- attention, including their nonlinear properties, such as ran, Safier and Stark [263, 264], and later generalized memory effects [10, 35, 37, 41, 42, 53, 54, 104, 171, 176, to spacetimes with plane symmetry [313, 314, 320]. It 230, 290, 303, 325, 335, 337, 338]. Such effects might was found that the polarizations of one GW can be ro- be possibly detected by the Laser Interferometer Space tated by the other, due to their nonlinear interaction, an Antenna (LISA) [142]. analogue of the electromagnetic Faraday rotation. Since From these observations, another important result is then, detailed studies of gravitational Faraday rotation that the individual masses of such observed black hole- have been carried out in various situations, including black hole (BBHs) binaries can be much larger than what Einstein-Maxwell waves [13], rotating cylindrical GWs were expected previously both theoretically and observa- [259], and more recently, cylindrical gravitational soliton tionally [141, 149, 254], and various scenarios have been waves [305]. proposed [2, 326]. For such massive binary systems, their mHz emission of GWs from the early inspiral phase is strong enough to be visible to space-borne detectors, such A. Cylindrically symmetric spacetimes as LISA [118], TianQin [225] and Taiji [271], for several years prior to their coalescence [279]. In fact, because of It is interesting to note that, despite the extensive the poor constraints on the rates and the standard as- studies of cylindrical spacetimes in the past several sumption of 10 solar mass objects, BBHs had never been decades, in the literature there still exist confusions in considered as particularly interesting sources of space- the definition of cylindrically symmetric spacetimes and borne detectors, and multi-band GW observations had (wrong/false) claims. A reason for that is probably con- been only proposed in the context of seed or interme- nected with different kinds of evolutions that remain diate mass BBHs [9, 278]. But, now analysis of the cylindrical symmetry: motions around and along the BBH population as constrained by current LIGO/Virgo symmetry axis. In both cases the velocities of these mo- demonstrated that such space-borne detectors can see tions may depend on the radial coordinate in an arbitrary thousands of such systems, with a variety of profound manner. scientific consequences. In particular, observations of the For example, it has been claimed by various authors same signal in two different detectors provides an efficient that the Kompaneets-Jordan-Ehlers-Kundt (KJEK) independent way to cross check and validate the instru- metric [197, 206], ments, which is particularly valuable for a space-based ds2 = e2(γ−ψ) dt2 − dr2
−e2ψ (dz + ωdφ)2−W 2e−2ψdφ2, detector. (1.2) More important, multi-band detections will enhance is the most general form that describes cylindrically sym- the potential of gravity tests in the strong, dynamical metric spacetimes, where γ, ψ, ω and W are functions of field regime of merging BBHs. Massive systems will t and r only. However, it is clear that this metric does not be observed by ground based detectors with high merg- include the rotating cylindrical GW spacetimes studied ing signal/noise rates, after being tracked for years by by Mashhoon and Quevedo (MQ) [233, 234, 266], space-borne detectors in their inspiral phase. The two 2 2(γ−ψ) 2 2
2 −2ψ 2 2ψ 2 portions of signal can be combined to push the search ds = e dt − dr −W e (dφ + ωdt) −e dz . for deviations from GR. In particular, joint observa- (1.3) tions of BBHs with a total mass larger than ∼ 60 solar To clarify these issues, it is necessary first to give a masses by LIGO/Virgo and space-borne detectors can rigorous definition of cylindrically symmetric spacetimes. potentially improve current bounds on dipole emission Unfortunately, this inevitably involves the introduction from BBHs by more than six orders of magnitude [18], of some rather abstract group theory and and in certain which will impose sever constraints on various theories sense complicated geometric conceptions. For the sake of gravity [30, 111]. These include scalar-tensor theo- of reader’s convenience, in the following we shall carry ries [90, 110, 192, 195, 334, 339], screened modified grav- this out in several steps. In doing so, it will become clear ity [222, 340, 341], Einstein-dilaton Gauss-Bonnet grav- that both of the above metrics describe spacetimes with ity [333], Einstein-aether theory [163, 220, 345], theories cylindrical symmetry, and that none of them actually with parity violations [265, 343, 344]. For more details, represents the most general metric of cylindrical space- we refer readers to the review articles [31, 32]. times. In particular, to obtain metric (1.2), an additional These studies recently have gained additional mo- condition that the two-dimensional abelian Lie group G2 1 menta, due to the close relations between asymptotically acts orthogonally transitively must be imposed , in ad- symmetric theorems of soft gravitons and GW memory dition to the general definition of cylindrically symmetric [178, 296]. Such studies in general have been carried out in spacetimes with plane symmetry, but to shed some 1 new light on the nonlinear effects, recently cylindrically For the definition of a Lie group Gn, see, for example, [8, 294] 4 spacetimes, while in metric (1.3) such an assumption is there exist local coordinates (t, x, y, z), such that u = −1 α left out, but with other further constraints. N (t, x)∂t, uαdx = N(t, x)dt, ξ = ∂y, η = ∂z, and To show the above claims clearly, let us start with the 2 2 2 i j definition of the cyclical symmetry. ds = N (t, x)dt − gij(t, x)dx dx Definition 1. A spacetime (M, g) is said to have a 2k 2 2h 2 2 h 2 cyclical symmetry, if and only if the metric is invariant = e dt − e dx − W f (dy + w1dz + w2dx) under the effective smooth action SO(2) × M → M of 2i +f −1 (dz + w dx) , (1.6) the one-parameter cyclic group SO(2) [93]. 3 Definition 2. A cyclically symmetric spacetime is said to be axially symmetric, if the set of fixed points where k, h, W, f and wi are all functions of t and x (i.e., points that remain invariant under the action of the only. It can be shown that this metric includes both group) of the isometry is not empty. The set of fixed metric (1.2) and metric (1.3) as its particular cases. points is referred to as the symmetry axis [92]. In particular, when G2 acts orthogonally transitively, Note that Definition 2 implicitly assumes that there we have w2 = w3 = 0, and then (1.6) reduces to exists at least one fixed point in (M, g). It can be shown 2 2k 2 2h 2 2 h 2 −1 2i that the set of fixed points must be an autoparallel, two- ds = e dt − e dx − W f (dz + w1dφ) + f dφ . dimensional timelike surface [93, 232]. Furthermore, the (1.7) infinitesimal generator ξα of the axial symmetry is space- like in a neighborhood of the axis, and that the so-called Note that in writing the above expression, we had made elementary flatness condition holds [232], the replacement (y, z) → (z, φ). Clearly, the two-surfaces of t, x = const. is gauge-invariant under the coordinate X X gαβ ,α ,β → −1, (1.4) transformations, 4X t = f(t0, r), x = g(t0, r), (1.8) as r → 0+, which ensures the standard 2π-periodicity of the axial coordinate near the axis, where r = 0 denotes where f(t0, r) and g(t0, r) are arbitrary functions of their the location of the axis, indicated arguments. Then, using this freedom, we can α always set gt0t0 = −grr and gt0r = 0, so the resulted X ≡ |ξ ξα| + → 0, (1.5) r→0 metric takes precisely the form of (1.2). In this case, there α and ( ),α ≡ ∂/∂x . Note that in this paper we shall adopt exists a general no-go theorem regarding the existence of the conversions: Greek letters run from 0 to 3 (µ, ν, ... = black holes [319]. 0, 1, 2, 3); Latin letters run from 1 to 3 (i, j, ... = 1, 2, 3); Theorem 2. Let (M, g) be a four-dimensional Rie- and repeated indices represent sum over their individual mannian spacetime obeying Einstein’s field equations, domains. Rµν − (R/2)gµν − Λgµν = κTµν , coupled with a cosmo- To define the cylindrical symmetry, in addition to ξ logical constant Λ. Assume that the spacetime possesses we assume that the spactime (M, g) contains another two commuting spacelike orthogonally transitive Killing spacelike Killing vector η, which generates together with vectors ξ and η. Then, (M, g) contains neither outer ξ a G2 group. Then, one can prove the following: nor degenerate apparent horizons, if matter satisfies the Proposition 1. In an axially symmetric spacetime, dominant energy condition and Λ > 0. in addition to ξ if there is another Killing vector η, Note that when Λ = 0, only the existence of outer ap- then both Killing vectors commute, thus generating an parent horizons is excluded. Thus, in this case only de- Abelian G2 group [92]. generate apparent horizons are allowed. However, when With the above, we are ready to define cylindrically Λ < 0, both outer and degenerate apparent horizons can symmetric spacetimes. exist. This explains why all the topological black holes Definition 3. A spacetime (M, g) is cylindrically found so far are with Λ ≤ 0. For more details, we refer symmetric if and only if it admits a G2 on S2 group of readers to [319], and for a concrete example, see [318]. isometries containing an axial symmetry [92]. In addition, the above theorem does not require any of As mentioned above, the assumption on the existence the orbits of the two spacelike Killing vectors be closed. of two-surfaces orthogonal to the group orbits is not nec- As a result, the theorem can be equally applied to plane essary for the definition of cylindrical symmetry nor a symmetric spacetimes [166, 320]. consequence of it, but the existence of an axis is essential. When G2 acts non-orthogonally transitively, but one So, in the above definition, the orthogonal transitivity is of the two Killing vectors is hypersurface-orthogonal, the removed, and the line-element of cylindrically symmetric corresponding line-element can be obtained from (1.6) spacetimes is provided in the following theorem [311]. by setting either w1 = w2 = 0 or w1 = w3 = 0, which is Theorem 1. Suppose that a spacetime admits two referred to as the A(ii) class of solutions in [311], while spacelike commuting Killing vector fields ξ and η. Let the metric (1.7) is classified as the B(i) class of solutions. u be a hypersurface-orthogonal timelike vector field of In the current case, without loss of the generality, we only unit length, which is orthogonal to ξ and η. Then, consider the case w1 = w2 = 0, with which the metric 5

(1.6) reduces to B. Einstein’s field equations

2 2k 2 2h 2 2 h 2 −1 2i ds = e dt − e dx − W fdy + f (dz + w3dx) . The world constructed by Newton to describe his the- (1.9) ory of gravitation had much simplicity: Space is always equal to itself where particles move and act upon each Making the coordinate transformations (1.8) and using other. Furthermore, the gravitational effects propagate the gauge freedom, we can always set gt0r = gzr = 0, so with an infinitely large velocity. This concept of space the resulted metric takes the form, sounded not plausible to Einstein and he came with the brilliant idea that space is the gravitational field mean- 2 2k 2 2h 2 2  2 ds = e dt − e dr − W fdz ing the Newton space itself [270]. This needed a better i +f −1 (dφ + w dt)2 , (1.10) formulation of the field equations. Einstein, after a long 3 struggle, found that the Riemann geometry, first pro- posed by Gauss and then generalized to any dimensions after the replacement, (y, z) → (z, φ). Clearly, the above by Riemann, could properly describe the curvature of metric reduces to the MQ metric (1.3) when k = h 2. space produced by different matter distributions. By do- A class of solutions of this metric was studied recently ing this Einstein managed to express the gravitational by Chicone, Mashhoon and Santos for gravitomagnetic field equations in an arbitrary coordinate system. accelerators [101, 102, 235]. In Riemann’s geometry the parallel transport of a vec- When G acts orthogonally transitively, and both ξ 2 tor is proportional to its curvature which is described by and η are hypersurface-orthogonal and mutually orthog- a quantity called the Riemann curvature tensor R . onal, the corresponding spacetimes are described by met- αβγδ The justification of calling it curvature lies in the fact ric (1.6) with w = w = w = 0, and with the gauge 1 2 3 that it vanishes if and only if the space is flat. freedom (1.8), the metric can be further written as Einstein [136] derived the field equations in the year ds2 = e2k dt2 − dr2
− W 2 fdz2 + f −1dφ2
. (1.11) 1915 which are given by 3, This is referred to as the B(ii) class of solutions in [311] 1 Rαβ − Rgαβ = κTαβ. (1.13) and represents the Einstein-Rosen (ER) cylindrical GWs 2 [137], which are linearly polarized, that is, only the “+” They show how the space curvature, represented by the polarization exists [166]. Ricci tensor R = Rγ and its scalar R = gαβR , It should be noted that in spacetimes with cylindrical αβ αγβ αβ are brought about by the source of curvature, the matter symmetry, closed timelike curves (CTC’s) can be easily distribution given by the energy momentum tensor T . formed. To guarantee their absence, one usually also αβ In fact not only space curves but time as well and for this imposes that the condition, reason we call spacetime where Einstein field equations α dwell. The coupling constant in normalized units, the ξ ξα < 0, (1.12) κ velocity of light c = 1 and the Newton’s constant G = 1, holds in the whole spacetime. values κ ≡ 8πG/c4 = 8π. This system of equations, In addition, sometimes it is also required that the called Einstein’s field equations, constitutes a set of 10 spacetime be asymptotically flat in the radial direction, partial differential equations with respect to the metric especially for the cases in which the sources are confined gαβ. In the vacuum case where Tαβ = 0, These equations within a finite region [29]. reduce to Finally, we note that there are physically realistic situ- ations where an axis may not exist, as it may be singular Rαβ = 0. (1.14) and so not part of the manifold, or the topology of the manifold may be such that no axis exists. Barnes noted When spacetime is deprived of sources producing grav- [20] that in this case the requirement of the existence of itational fields, and thereby producing curvature, the an axis in Proposition 1 can be removed and the following spacetime is flat and given by the metric ηαβ, called holds (See also [34] and the footnote given in [20]). Minkowski metric, Proposition 2. Any two-dimensional Lie transfor- ds2 = η dxαdxβ = dt2 − dx2 − dy2 − dz2, (1.15) mation group that contains a one-dimensional subgroup αβ whose orbits are circles must be Abelian. with t being the time coordinate and x, y and z the Carte- sian coordinates. Hence, the Newtonian limit, which one

2 Note that in the 2-surfaces of z, φ = const., we can always use the gauge freedom (1.8) to write the reduced metric in 3 In this review, we shall adopt the conventions of [128]. In partic- the conformally-flat form, e2k dt2 − dr2
. But, now the sec- ular, the signature of the metric will be (+, −, −, −), and when ond part of the metric (1.10) will in general take the form, the cosmological constant Λ does not vanish, Einstein’s equations 2 h 2 −1 2i W fdz + f (dφ + w2dr + w3dt) . 1 read, Rαβ − 2 Rgαβ − Λgαβ = κTαβ . 6 expects to be obtained from Einstein field equations, is vacuum field with the following line element ds, produced when the metric is given by   2 α β 2M 2 ds = gαβdx dx = 1 − dt gαβ = ηαβ + fαβ, (1.16) r  2M −1 with f being a small deviation from the Minkowski − 1 − dr2 − r2(dθ2 + sin2 θdφ2)(2.1), αβ r spacetime. In this paper, we shall present an updated review on where the spherical coordinates are numbered x0 = t, cylindrically symmetric spacetimes. However, due to the x1 = r, x2 = θ and x3 = φ. In the Newtonian approxi- vast scope of the field carried out in the past several mation (2.1) becomes decades and the limitation of the scope of the review, it is impossible to cover all the subjects studied so far. So, g00 = 1 + 2U, (2.2) one way or another one has to make a choice on which subjects that should be included in a brief review, like the where U is the Newtonian potential and comparing to current one. Such a choice clearly contains the review- (2.1) we have ers’ bias. In addition, in this review we do not intend M to exhaust all the relevant articles even within the cho- U = − , (2.3) sen subjects, as in the information era, one can easily r find them, for example, from the list of the citations of hence the only parameter stemming from the integration relevant articles. With all these in mind, we would like of the field equations for a spherically symmetric vacuum first to offer our sincere thanks and apologies to whom spacetime is the Newtonian mass M. his/her work is not mentioned in this review. In addi- tion, there have already existed two excellent books to cover the subjects extensively [167, 294] and a brief re- A. LC solutions view on gravitational collapse in cylindrically symmetric spacetimes [237]. The second vacuum solution of Einstein field equations The rest of article is organized as follows: In the next (1.14) was obtained a few years later, in 1919 by LC four sections, we shall review the general properties, both [218], corresponding to a cylindrical vacuum spacetime, local and global, and possible sources of the Levi-Civita and can be cast in the form [Cf. Appendix A], solution (Sec. 2), coupled with an electromagnetic field (Sec. 3), with a cosmological constant (Sec. 4), or with 1 ds2 = ρ4σdt2 − ρ4σ(2σ−1)(dρ2 + dz2) − ρ2(1−2σ)dφ2, a perfect fluid (Sec. 5). Then, in the next two sections, a we shall provide a general review on the Lewis vacuum (2.4) solution (Sec. 6), and its generalization to include a cos- where the cylindrical coordinates are numbered x0 = t, mological constant (Sec. 7). Finally, In Secs. 8 and x1 = ρ, x2 = z and x3 = φ. In the Newtonian approxi- 9, we shall provide a review on wormholes, gravitational mation, correspondingly to (2.2) we have, collapse, polarizations of cylindrical GWs and gravita- 2U tional Faraday rotations. An appendix is also included, g00 = e , (2.5) in which Einstein’s field equations of static cylindrical spacetimes are presented. and compared to (2.4) for this case we have the Newto- nian potential,

U = 2σ ln ρ. (2.6) II. LEVI-CIVITA (LC) VACUUM SPACETIME Hence, from this expression we clearly see that for small Einstein’s field equations (1.13) are highly nonlinear values of σ it is the Newtonian mass per unit length as thus imposing huge difficulties in finding their solutions. produced by an infinitely long homogeneous line mass, as Furthermore, once found their solutions, even more diffi- observed by LC himself. Ever since much has been writ- culties arise in interpreting them. In spite of great efforts ten by researchers trying to grasp its physical and geo- to grasp their physical meanings, as is well shown in [293] metrical interpretations. However, this endeavor proved and [167], the majority or even almost the totality of to be difficult and uncertain. these solutions are still not well understood given their Only in 1958 did Marder [231] established that the so- high nonlinearity. A simple example like the spherical lution (2.4) contains two arbitrary independent parame- vacuum field for a point mass, presented below, which is ters, σ and a, differently from the Newtonian fields in perhaps the solution most studied in GR, still bears weird wich there is only one independent parameter σ. We also and unexpected properties as we will see in the following call the attention that for the spherical case in the rel- by studying cylindrically symmetric vacuum fields. ativistic and Newtonian theories there appears just one Schwarzschild [276] in 1916 obtained the first vacuum parameter. This fact already suggests some harder diffi- solution to (1.14), describing the spherically symmetric culties in understanding cylindrically symmetric fields. 7

In the Newtonian approximation the parameter a is as- Consider 1/2 < σ < ∞, which implies h(0) = 0 and sociated with the constant arbitrary potential that exists l(r) diverging when r → 0. Then we can interpret m as in the Newtonian solution, and can be given any value. the angular coordinate φ with topological identification However, in GR Bonnor [45] observed that a is dressed of every φ with φ + 2π, and n as an axial coordinate with a relevant global topological meaning, and cannot −∞ < z < ∞ with an = 1 by again rescaling z, and be removed by coordinate transformations. (2.10) becomes In the following we review the main properties and 2 4σ 2 2 − 2(2σ−1) 2 1 4σ(2σ−1) 2 physics that lie behind the LC spacetime, which so far ds = R Σ dt − dr − R Σ dz − R Σ dφ , have been grasped up to the present time in a large num- a ber of papers. We are aware that these results sometimes (2.13) appear contradictory since some interpretations collide where we replaced am for a. with others. In fact this is one of our main motivations In both cases (2.12) and (2.13) where σ > 0 we have to deepen into Einstein’s theory. g00 → 0 as r → 0, indicating an attractive singularity. While assuming σ < 0, we obtain g00 → ∞ as r → 0, indicating a repulsive singularity. B. Nature of the coordinates of the LC solutions The invariant quantity, under coordinate transforma- tions, built out of the Riemann curvature tensor given by R ≡ RαβγδR , called the Kretschmann scalar, is a The LC metric given by (2.4) can be written in the αβγδ good indicator of singularities. Calculating R for metrics form [186] (2.12) and (2.13) we obtain  1  1 ds2 = ρ4σdt2−ρ4σ(2σ−1) dρ2 + dm2 − ρ2(1−2σ)dn2, 64σ2(2σ − 1)2 am an R = 3 4 . (2.14) (2.7) Σ r where −∞ < t < ∞ is the time and 0 ≤ ρ < ∞ the radial From (2.14) we see that R → ∞ as r → 0, and R → 0 coordinates, and σ, am and an are constants. The nature for σ = 0, 1/2 and ∞. of the coordinates m and n, so far unspecified, depends Metric (2.12) for σ = 0 becomes upon the behaviour of the metric coefficients. Either am or a can be transformed away by a scale transformation 1 n ds2 = dt2 − dr2 − dz2 − r2dφ2, (2.15) depending upon the behaviour of the coordinates m and a n, thus leaving the metric with only two independent parameters. In order to find that behaviour we transform representing the Minkowski spacetime when a = 1 in the radius ρ into a proper length radial coordinate r by cylindrical polar coordinates r, z and φ. However, if defining (2.15) has √a > 1 there is an angle deficit 2πδ given by δ = 1 − 1/ a, producing flat spacetime everywhere ex- ρ2σ(2σ−1)dρ = dr, (2.8) cept along the axis r = 0, which is interpreted as a cos- mic string. The deficit can represent the tension of the thus obtaining string with mass per unit length µ given by µ = δ/4. If there is an angle excess a < 1 it would represent a cos- 1 2 ρ = R Σ ,R = Σr, Σ ≡ 4σ − 2σ + 1, (2.9) mic string under compression with µ < 0. Hence the constant a is directly linked to the gravitational ana- and metric (2.7) becomes log of Aharonov-Bohm effect [133]. This effect shows that gravitation depends on the topological structure of 2 2 2 2 2 ds = f(r)dt − dr − h(r)dm − l(r)dn , (2.10) spacetime giving rise to an angular deficit δ as in the elec- tromagnetic Aharonov-Bohm effect, where a (classical) with non-observable quantity (the vector potential) becomes

4σ 1 4σ(2σ−1) 1 2(1−2σ) observable (part of it) through a quantum non-local ef- f(r) = R Σ , h(r) = R Σ , l(r) = R Σ . fect. Its gravitational analog allows a (Newtonian) non- am an (2.11) observable quantity (the additional constant to the New- Consider 0 < σ < 1/2, which implies for this range tonian potential) to become observable in the relativistic h(r) diverging when r → 0 and l(0) = 0. Then we can theory through the angular deficit in strings. For a re- interpret m as the axial coordinate −∞ < z < ∞ with view in cosmic strings see [189]. am = 1 by rescaling z, and n as the angular coordinate φ In the case σ = 1/2 the two metric coefficients h and l with the topological identification of every φ with φ+2π, in (2.11) are constants, so both am and an can be set to and the metric (2.10) becomes unity. Then neither m nor n is entitled to be an angular coordinate, and the three coordinates r, m and n are 2 4σ 2 2 − 4σ(1−2σ) 2 1 2(1−2σ) 2 better visualized as Cartesian coordinates x, y and z. ds = R Σ dt − dr − R Σ dz − R Σ dφ , a Hence metric (2.7) can be written as (2.12) 2 2 2 2 2 2 where an is replaced by a. ds = z dt − dx − dy − dz , (2.16) 8 which is the static plane symmetric vacuum spacetime σ < ∞ is equivalent to the range 0 < σ < 1/2 and the obtained by Rindler [125, 169, 268]. coordinates z and φ switching their nature. In the next section we calculate the circular geodesics Hence we might have the following picture for the dif- for the above spacetime to try to get further understand- ferent values of σ. For small values of σ the metric (2.12) ing of the properties of different values that σ can attain with t and r constants describes cylindrical surfaces with for σ ≥ 0, since σ < 0 would in some way correspond to φ a periodic coordinate. As σ increases the cylindrical negative mass densities. surfaces open out and become infinite planes for σ = 1/2. For values of σ bigger that 1/2 the coordinate z becomes periodic forming new cylindrical surfaces perpendicular C. Geodesics to previous ones for 0 < σ < 1/2. Another interesting geodesic is the one that describes For the metrics (2.12) and (2.13) the circular geodesics the motion of the particle along the axis of symmetry 2 ˙2 z. These geodesics calculated with (2.12) (we restrict [121] haver ˙ =z ˙ = 0 and gtt,rt˙ + gφφ,rφ = 0, where the dot stands for differentiation with respect to the affine the calculation of this metric since (2.17) is equivalent) parameter s. The geodesic angular velocity is defined by produce ω = φ/˙ t˙, and the only nonzero component of its velocity φ √ 4σ(1 − 2σ) r˙z˙ is the tangential one, given by W = ω/ gtt with its z¨ = . (2.22) 2 α 4 Σ r modulus defined by W ≡ W Wα [184]. In the case 0 ≤ σ < 1/2, from (2.12) we obtain It means that particles increase their speed along z when distancing radially from the axis, while diminishing their 2σ ω2 = aR2(4σ−1)/Σ, (2.17) axial speed when moving radially towards the axis. This 1 − 2σ result indicates that a force parallel to the z axis appears. 2σ W 2 = . (2.18) In the flat case σ = 0 such an effect vanishes, bring- 1 − 2σ ing out its non-Newtonian nature. We further discuss this weird geodesic property in the section concerning We note that for a given σ the velocity W in (2.18) is the Lewis spacetime. the same for all circular geodesics, in agreement with the corresponding Newtonian gravitation. Furthermore, we For radial geodesics it has been shown [184] that there see that W increases monotonically with σ, that is from exist timelike particles approaching z that are reflected σ = 0 producing W = 0, to σ = 1/4 attaining W = 1, the at r = rmin and move outwards until attainingr ˙ = 0 at speed of light, and finally σ = 1/2 producing geodesics rmax repeating endlessly this trajectory. This motion is with W = ∞. For small σ and a = 1 from (2.17) and called confinement of test particles. (2.18) we obtain the Newtonian limit W = rω. In the next section we see some further limits that LC In the case 1/2 < σ < ∞, from (2.13) we obtain metric satisfies. 1 ω2 = aR8σ(1−σ)/Σ, (2.19) 2σ − 1 D. LC spacetime as a limiting case of the γ 1 spacetime W 2 = . (2.20) 2σ − 1 In cylindrical coordinates, static axisymmetric solu- With σ increasing beyond 1/2, we note from (2.20) that tions of Einstein’s vaccuum equations are given by the W diminishes, attaining W = 1 for σ = 1 and W = 0 for Weyl metric [328] σ = ∞. In other words, the circular geodesics are timelike when ds2 = e2λdt2 − e−2λ e2µ(dρ2 + dz2) + ρ2dφ2 , (2.23) either 0 < σ < 1/4 or σ > 1, are lightlike when σ = 1/4 or 1, and are spacelike when 1/4 < σ < 1. with λ(ρ, z) and µ(ρ, z) satisfying If we redefine σ by 1 1 λ,ρρ + λ,ρ + λ,zz = 0, (2.24) σ = , (2.21) ρ 4¯σ and then metric (2.13) with a rescaling of its coordinates be- 2 2 comes (2.12), hence (2.19) and (2.20) reduce to (2.17) µ,ρ = ρ(λ,ρ − λ,z), (2.25) and (2.18). This means that the parameter range 1/2 < µ,z = 2ρλ,ρλ,z, (2.26) where the comma stands for partial derivation. Synge writes [297], as the most amazing fact, that (2.24) is just 4 This W is different from the one appearing in the metric (1.2) the Laplace equation for λ in the Euclidean space. Met- or that in (1.3). ric (2.23) with the general solution of (2.24), (2.25) and 9

(2.26) is usually referred to as the γ metric, and the cor- m → ∞ the γ metric is the same as the LC metric. De- responding functions λ and µ are given by tails of these long calculations are given in [182]. Here we come to an interesting and, so far, weird result showing R + R − 2mγ e2λ = + − , (2.27) how apparently unexpected results can stem from long R+ + R− + 2m known results like the Schwarzschild and LC solutions. 2 (R + R + 2m)(R + R − 2m)γ When the density per unit length of the rod σ = γ/2 e2µ = + − + − (2.28), has the value γ = 1, or the mass density per unit length 4R+R− σ = 1/2, it becomes the Schwarzschild spherically sym- where metric spacetime, and in the limit m → ∞, it becomes the Rindler static plane symmetric vacuum spacetime. 2 2 2 R± = ρ + (z ± m) , (2.29) This is a remarkable result. For all the different limiting metrics that LC space- and γ and m are two integrations constants. These so- time can undergo see [182]. The limits for the circular lutions were first found by Bach and Weyl in 1922 [16]. geodesics of the γ spacetime to the LC spacetime are well Calculating its Newtonian potential (2.5) we obtain studied in [184]. In the next section we make a brief review of possible R + z − m U = γ ln − , (2.30) sources to the LC spacetime. R+ + z + m which corresponds to a potential due to a line segment of length 2m and mass per unit length γ/2 symmetrically E. Sources producing LC spacetime distributed along the z axis. Hence the total mass M of the line segment is M = γm. The particular case The LC spacetime, as we saw, contains two essential γ = 1 corresponds to the Schwarzschild metric. This can constants denoted by a and σ. The constant a refers be seen by taking spherical Erez-Rosen coordinates [139], to an angle deficit or excess that produces cosmological defined by strings. However, it is σ that presents the most serious obstacles to its interpretation. Indeed, for small values ρ2 = (r2 − 2mr) sin2 θ, z = (r − m) cos θ, (2.31) 0 < σ < 1/4, LC describes the spacetime generated by an infinite line mass, with σ as mass per unit coordi- and the γ metric becomes nate length. When σ = 0 the spacetime is flat. However, circular timelike geodesics only exist for 0 < σ < 1/4, be- 1 ds2 = F dt2 − Gdr2 + Hdθ2 + (r2 − 2mr) sin2 θdφ2 , coming null for σ = 1/4 and being spacelike for σ > 1/4. F (2.32) Furthermore, as the value of σ increases from 1/4 to 1/2 where the corresponding Kretschmann scalar (2.14) diminishes monotonically and vanishes when attaining σ = 1/2 im- γ  2m plying that the spacetime is flat. F = 1 − , (2.33) r Thus, the question is what does the LC spacetime rep- 2  2 γ −1 resent outside the range 0 ≤ σ ≤ 1/4 ? r − 2mr First, we observe that there are known physical sat- G = 2 , (2.34) r2 − 2mr + m2 sin θ isfactory fluid sources for the LC spacetime satisfying 2 (r2 − 2mr)γ boundary conditions for both ranges of σ (see for exam- H = 2 . (2.35) ple [188, 293]). On the other hand, the fact that the (r2 − 2mr + m2 sin2 θ)γ −1 Kretschmann scalar decreases with increasing σ may not Now we can easily check that for γ = 1 the metric (2.32) be associated with the strength of the gravitational field. reduces to the Schwarzschild metric (2.1). Instead, it could be associated with the acceleration of If we want to compare the γ metric to the LC metric, a test particle that would measure more suitably the in the limit when its length segment m → ∞, one notices strength of the gravitational field, which is the case for that by taking directly this limit in (2.27) and (2.28) the 1/4 < σ < 1/2. metric diverges. For this reason one can use the Cartan A possible interpretation of the LC solution is a space- scalars approach to obtain a finite limit. These scalars time generated by a cylinder whose radius increases with are the components of the Riemann tensor and its co- σ and tends to infinity as σ approaches 1/2. This in- variant derivatives calculated in a constant frame. Two terpretation suggests that when σ = 1/2 the cylinder metrics are equivalent if and only if there exist coordinate becomes a plane. and Lorentz transformations which transform the Cartan It might be instructive to consider a static cylinder scalars of one of the metrics into the Cartan scalars of filled with anisotropic perfect fluid and calculate its mass the other. Although the Cartan scalars provide a local per unit length by using the junction conditions to its characterization of the spacetime, global properties such exterior LC spacetime. The Whittaker formula [330] for as topological defects do not appear in them. By doing the active gravitational mass per unit length ν of a static all this procedure one can prove that locally in the limit distribution of perfect fluid with energy density µ and 10 principal stresses Pr, Pz and Pφ inside a cylinder of sur- III. STATIC ELECTROVACUUM face S is SPACETIMES Z rS √ ν = 2π (µ + Pr + Pz + Pφ) −gdr. (2.36) A. Electromagnetic fields in cylindrical 0 spacetimes Considering a static cylindrically symmetric metric In this section we discuss the gravitational fields whose ds2 = −Adt2 + dρ2 + Cdz2 + Dρ2dφ2, (2.37) only source is a free electromagnetic field described in which A, C and D are only functions of ρ, from Ein- by Maxwell’s theory. This field in curved spacetime stein’s fields equations we obtain [165] is characterized by the antisymmetric tensor Fµν = ∂ A − ∂ A (A is the vector 4-potential) satisfying   ν µ µ ν µ A,rr 1 A,r A,r 2 C,r D,r the Maxwell equations − − − − = (µ + P + P + P ) . A 2 A A ρ C D κ ρ z φ (2.38) αµ µ ∗ αµ Substituting (2.38) into (2.36) we have the simple expres- ∇αF = 4πj , ∇α F = 0, (3.1) sion µ   where ∇α is a covariant derivative, j is the electric 4- 1 A,r √ ∗ αβ αβµν √ ν = −g , (2.39) current, and F = −ε Fµν /(2 −g) is the axial 4 A αβµν S tensor dual to Fµν (ε is the LC completely antisym- where regularity conditions on the axis of symmetry [188] metric unit object). In a vacuum region of a spacetime µ have been assumed. Now taking the LC metric (2.12) for we have j = 0. the exterior spacetime of the cylindrical surface S, and Nonzero components Fµν , compatible with the geome- imposing that it satisfies Darmois’ junction conditions try (A.1), depending on the radial coordinate u only can [119, 188] with its interior spacetime described by (2.37), have three different directions: it amounts to impose the continuity of the metric func- 2 • Radial (R) fields: electric, F01(u) (so that E = tions and its derivatives on S. By doing so, from (2.39) 10 2 we obtain F01F ), and magnetic, F23(u) (so that B = F F 23). σ 23 ν = √ , (2.40) a 2 • Longitudinal (L) fields: electric, F02(u)(E = 20 2 13 where a is the constant defined in (2.12). From (2.40) F02F ), and magnetic, F13(u)(B = F13F ). we have that when a > 1, there is a topological angle • Azimuthal (A) fields: electric, F (u)(E2 = defficit, then ν < σ, and if a < 1 there is an angle excess, 03 F F 30), and magnetic, F (u)(B2 = F F 12). producing ν > σ. However, when there is no topological 03 12 12 defect, a = 1, then ν can in fact be interpreted as mass Here E and B are the absolute values of the electric per unit length of its source. Furthermore, since with field strength and magnetic induction, respectively. Self- cylindrical sources no black holes are formed, one might gravitating static, cylindrically symmetric configurations conclude that the minimum mass per unit length to form of such electromagnetic fields have been considered in a black hole satisfies the constraint ν > 1/2. This result [62, 63, 238, 272, 301] ( in which see also references to would fulfill the present knowledge of black holes that a earlier work, and more historical details can be found in lower mass limit is required for its formation. We are [294]). In our presentation we follow the lines of [62, 63]. aware that ν is model-dependent and cannot be given a The Maxwell equations follow from the least action general meaning, but nonetheless it fulfills some of the principle with the Lagrangian expected properties. A last comment on sources for the LC spacetime is µν 2 2 Fµν F E − B the fact that conformally flat static sources do not admit Lem = − ≡ . (3.2) Darmois’ matching conditions satisfied for an exterior LC 16π 8π spacetime. This result is proved in [181]. Conformal The same Lagrangian leads to the following well-known flatness of a metric means that its corresponding Weyl form of the stress-energy tensor (SET) to be used as a tensor vanishes. The interpretation of the Weyl tensor source of gravity in Einstein’s equations: is that it describes the purely gravitational characteris- tics of the source. This interesting result means that 1 T ν = −4F F να + δν F F αβ
. (3.3) a static cylindrical source deprived of its purely gravi- µ 16π µα µ αβ tational character cannot be smoothly matched to the exterior LC spacetime. It is conspicuous that for spheri- This expression leads to an important ( result: it cal symmetry this result does not stand, since there are turns out that the metric (A.1) admits only a single sort conformally flat static spherical sources matched to the of electromagnetic field (R, L, or A, defined above). Any Schwarzschild spacetime [180]. simultaneous existence of fields of different directions is 11 forbidden due to the emergence of off-diagonal compo- Equations (3.7) are trivial while (3.8) is a Liouville equa- ν nents of Tµ , whereas the left-hand side of Einstein’s equa- tion which is easily solved, giving tions (1.13) ( has only diagonal components. For exam- ple, a nonzero electric R-field (F01 6= 0) together with a µ(u) = au − γ(u), β(u) = bu − γ(u), (3.10) 3 magnetic L-field (F13 6= 0) leads to T0 6= 0.  However, simultaneous existence of electric and mag-   (1/k) sinh ku, k > 0, netic fields of the same direction is admissible, and due to  −γ(u)  the well-known electric-magnetic duality, such fields lead e = s(k, u) ≡ u, k = 0, (3.11) ν to similar contributions to Tµ (in the general case there   appears the sum E2 + B2). Therefore in what follows  (1/k) sin ku, k < 0, we consider each sort of fields (R, L, A) separately, and for definiteness, speak of electric R-fields (F01 = −F10), where a, b, k are integration constants, which can take magnetic L-fields (F13 = −F31) and magnetic A-fields any real values, and three more constants have been re- (F12 = −F21), bearing in mind that in each case mix- moved by rescaling t and z and by choosing the zero point tures of electric and magnetic fields of the same direction of the coordinate u. Furthermore, (3.9) is a first integral can be obtained by duality rotations without changing of (3.7) and (3.8), which, after substitution of (3.10) and the SET and therefore the metric. With this ansatz the (3.11) yields a simple relation between the integration Maxwell equations reduce to the first part of (3.1). In constants vacuum (jµ = 0) they lead to ab = k2 sign(k). (3.12) d √ µ1
µ1 qµ −gF = 0 ⇒ F = √ , qµ = const, du −g Thus we have a solution with three significant integra- (3.4) tion constants q, a, b, which splits into three subfamilies where µ = 0, 3, 2 for R-, L-, A-fields, respectively. according to the sign of k. To solve Einstein’s equations, we choose the harmonic radial coordinate u (A.4) in the metric (A.1), such that 1. k > 0, so that ab = k2. The metric takes the form √ α+β+γ+µ α(u) = β(u) + γ(u) + µ(u). Then −g = e = 2 2 2 2α k dt q h e , and the solutions (3.4) can be written as ds2 = − sinh2 ku e2(a+b)udu2 q2 sinh2 ku k2  01 31 21 −2α F ,F ,F = e {qR, qL, qA} , (3.5) i +e2audz2 + e2budφ2 . (3.13) where in each curly bracket one should consider a single term for each sort of field, and the constants qN have Without loss of generality we may take u ∈ R+ = (0, ∞) been renamed accordingly. In what follows, describing as the range of u, and evidently u = 0 corresponds each of these fields separately, we denote by simply q any to the axis since the cylindrical radius r ≡ eβ van- of these “charges” without a risk of confusions. ishes there. The electric field strength diverges there ν α 2 Since the tensor Tµ has zero trace, T = Tα = 0, Ein- as E ≈ |q|/u , therefore u = 0 can be interpreted as ν ν stein’s equations can be used in the form Rµ = 8πTµ , the location of a stretched electric charge. Since there ν 2γ 2 with the expressions (A.2) for Rµ. They are easily solved g00 = e ∼ 1/u → ∞, it is a repulsive singularity for by the same method in all three cases, but the resulting small test bodies. geometries are substantially different. The behavior of the metric at the other end, u → ∞, depends on the interplay of the constants a and b: while in all cases e2γ → 0, the radius eβ and the longitudinal B. Geometry with a radial electric field scale factor eµ may tend to zero, infinity or finite values, as is apparent in (3.13). In the cases where r → 0 as In this case, as one easily finds, the SET (3.3) takes u → ∞, we have to conclude that the spacetime has the form there one more axis besides u = 0, where again E → ∞, 2 but since it is the other end of the electric lines of force, it ν q −2µ−2β Tµ = e diag(1, 1, −1, −1), (3.6) is a location of a charge of opposite sign to that existing 8π at u = 0. µ where q = qR may be interpreted as a linear electric An exceptional case is a = b = −k, such that both e charge along the z axis. Due to high symmetry of the and eβ are finite at u = ∞, and this cylinder turns out SET (3.6), independent components of Einstein’s equa- to be a Killing horizon. This is confirmed by a transition ν ν −2ku tions Rµ = 8πTµ reduce to to the new radial coordinate x = 1 − e , after which the metric (3.13) takes the form β00 = µ00 = −γ00, (3.7) 2 2 2 2 2 2 00 2 2γ 4(1 − x)k q x dx q x γ = q e , (3.8) ds2 = dt2 − − dz2 + dφ2
. q2x2 16k4 1 − x 4k2 β0γ0 + β0µ0 + γ0µ0 = q2e2γ . (3.9) (3.14) 12

and consequently at spatial infinity (u → ∞) the electric field strength E may grow (if a < 0) or vanish (if a > 0). 3. k < 0, so that ab = −k2. The metric is

k2 dt2 q2 h i ds2 = − sin2 ku e2(a+b)udu2+e2audz2+e2budφ2 . q2 sin2 ku k2 (3.16) Without loss of generality the range of u can be chosen as 0 < u < π/|k|, and its both ends correspond to re- pulsive (since eγ → ∞) axial (since eβ → 0) singularities with E → ∞, to be ascribed to axial electric charges of opposite signs. The constants a and b do not affect the qualitative nature of the configuration, being only involved in exponentials which are finite and smooth in this range of u. In all cases with two axes, the 2D subspaces t = const, z = const can be imagined as warped spheres with sin- gularities at their poles.

C. Geometry with a longitudinal magnetic field

In this case, as one easily finds, the SET (3.3) takes the form q2 T ν = e−2µ−2β diag(1, −1, 1, −1), (3.17) µ 8π

FIG. 1. Carter-Penrose diagram for the inverted black hole where q = qL can be ascribed to an electric current in spacetime The lines x = 0 are repulsive singularities, the in- the azimuthal (φ) direction, like that in a solenoid. Due tersecting lines x = 1 form a Killing horizon, and x = ∞ is a to high symmetry of the SET (3.6), independent compo- temporal infinity. nents of Einstein’s equations in full similarity with (3.7)– (3.9) reduce to Here we have the axis x = 0 that coincides with u = 0; −µ00 = −γ00 = β00 = −q2e2β, the cylinder x = 1 (at which u = ∞) is a horizon, be- β0γ0 + β0µ0 + γ0µ0 = q2e2β, (3.18) yond which, at x > 1, there is a homogeneous anisotropic (Bianchi type I) cosmology where x is a temporal coor- so excluding the insignificant integration constants by dinate, and x → ∞ is a highly anisotropic temporal in- rescaling the t and z axes and choosing the zero point finity where two scale factors, eµ and eβ, tend to infinity of u, we can write the solution as while the third one, eγ , turns to zero. Since, unlike the k2 Schwarzschild and other black hole spacetimes, a static γ = −β + au, µ = −β + bu, e2β = , region is here inside the horizon, this solution was named q2 cosh2 ku “inverted black hole” [63] (Fig.1). a, b, k = const, k > 0, ab = k2. (3.19) 2. k = 0, so that by (3.12) either a or b is zero (or both). Thus the metric can be written as If we choose, for definiteness, b 6= 0, then the metric takes q2 cosh2 kuh i the form ds2 = e2audt2 − e2(a+b)udu2 − e2budz2 k2 2 2 dt 2 2 2bu 2 2 2bu 2 k2dφ2 ds = 2 2 − q u e du + dz + e dφ . (3.15) − . (3.20) q u q2 cosh2 ku

Again u ∈ R+, the axis u = 0 is a repulsive singular- The range of u is u ∈ R; moreover, without loss of gen- ity pertaining to a charged thread, while the “far end” erality we can assume a > 0 and b > 0 (making them u → ∞ is either a spatial infinity (b ≥ 0, eβ → ∞) or negative is equivalent to changing u → −u). the second axis (b < 0, eβ → 0). The latter is now an at- The circular radius eβ → 0 as u → ±∞, that is, we tracting singularity where eγ → 0 but E ∼ e−µ−β → ∞. have two axes at two ends of the u range. The end u → The case a 6= 0, b = 0 differs from the previous one ∞ is, for any k, a, b infinitely distant since the proper where a = b = 0 by a changing scale along the z axis, length integral l = R eαdu diverges. The coefficients eµ 13 and eγ also diverge there, the latter meaning that the Despite the similarity between (3.20) and (3.25), their gravitational field is repulsive from this remote axis. The interpretations are quite different due to our assumptions magnetic induction B vanishes there. of the topology of z and φ axes. As to u → −∞, the axis there is at finite proper dis- In the metric (3.20), “on the right end” u → +∞ we tance from any point with u < +∞, since R eαdu con- have eγ → ∞ (repulsion), eβ → ∞ (spatial infinity in verges. The metric behavior is more diverse there: there the sense of growing circular radius), l = R eαdu → ∞ is repulsive (if a < k) or attractive (if a > k) singularity (an infinite radial distance), and the magnetic induction according to the behavior of eγ , while the scale along z B → 0. Meanwhile, eµ → 0 (shrinking longitudinal scale) behaves oppositely to eγ : it vanishes if a < k and blows on both ends. up if a > k. In addition, B → ∞ as u → −∞ in all cases On the other end, u → −∞, l < ∞, but it is a repulsive except a = b = k. singular axis (eβ → 0, eγ → ∞) only if a < k. If a > But the most interesting geometry is observed if a = k, we have, on the contrary, eγ → 0 (attraction) but b = k: the magnetic field there is finite, and there is eβ → ∞, and the whole configuration is wormhole-like no curvature singularity at u = ∞. In this case the [81] though with a singular end at u = −∞. As in the substitution case of L-fields, B → ∞ as u → −∞ in all cases except a = b = k. |q| |q| |q| γ β eku = ρ = r , t = t¯ , z =z ¯ , (3.21) If a = b = k, then at u = −∞ both e and e are 2k 2k2 2k2 regular, and to clarify the metric properties it is helpful 2ku transforms the metric to to substitute e = y, after which the metric reads 2  2 2  2 2 2 (1 + y) 2 4k dy 2 Ky 2 2 q 2 2 2 2 2
r 2 ds = dt − − dφ − dz , ds = (1 + ρ ) dt¯ − dr − dz¯ − dφ . K y (1 + y)2 k2 (1 + ρ2)2 4k2 (3.22) K = . (3.26) This metric is regular on the axis r = 0 if q2 = 4k3 q2 and coincides with the metric of Melvin’s “geon” [236]. Now, let us further transform the x coordinate to make Otherwise there is a conical singularity at r = 0, and this the metric of (y, z) 2-surfaces explicitly conformally flat: axis behaves as a cosmic string. In the limit r → ∞ (u → ∞), as can be easily verified, Ky Z 2k(1 + y)2 all curvature invariants vanish. Thus the metric (3.22) dl2 = (dF 2 + dz2),F = dy. y,z (1 + y)2 Ky has no curvature singularities. It describes a magnetic (or (3.27) electric) field configuration supported by its own gravi- To regularize it at y = 0 we can use the notion of confor- tational field (“a magnetic universe”), with or without a mal mappings on the complex plane using the function cosmic string on its axis. W = ln Z for W (ζ) = F + iz with ζ = ξ + iη. More precisely, in terms of F and z, we substitute D. Geometry with an azimuthal magnetic field 1 η F = ln(ξ2 + η2), z = arctan ⇒ 2K ξ For an azimuthal magnetic field, the SET (3.3) takes dξ2 + dη2 dF 2 + dz2 = , (3.28) the form K(ξ2 + η2) 2 ν q −2β−2γ and since ξ2 + η2 ∼ y as y → 0, the metric is nonsingular Tµ = e diag(1, −1, −1, 1), (3.23) 8π at y = 0. However, under this transformation the range of z is π/2 < z < π/2, and the whole (ξ, η) plane 2 where q = qA may be interpreted as a steady electric R current along the z axis. The set of equations is the maps in terms of (y, z) to this interval times y ∈ R+. same as for the L-field up to the replacement φ ↔ z, To extend the solution to z ∈ R, we have to consider β ↔ µ, and the solution may be written as a Riemannian surface (ξ, η) with in infinite number of sheets and ξ = η = 0 as a branching point. k2 Thus in this case we have a self-supported magnetic γ = −µ + au, β = −µ + bu, e2µ = , q2 cosh2 ku distribution in a spacetime with a branching-point sin- 2 gularity. a, b, k = const, k > 0, ab = k , (3.24) From the above analysis it can be seen that static cylin- q2 cosh2 kuh i drically symmetric Einstein-Maxwell fields are rather di- ds2 = e2audt2 − e2(a+b)udu2 − e2budφ2 k2 verse, and the most interesting and nontrivial configura- k2dz2 tions exist in special cases where there is a kind of equi- − . (3.25) q2 cosh2 ku librium between gravity and electromagnetism: a horizon in an R-field, Melvin’s nonsingular universe with an L- As before, u ∈ R; and without loss of generality we as- field, and a system with a branching point singularity in sume a > 0 and b > 0. an A-field. 14

To understand the solutions further, sources producing we find that such spacetimes were considered in [238]. It was shown 2 2 2 2 r 2 2 that, in contrast to the vacuum case, when the electro- ds ' dt − dz − dϕ − dr , (4.5) magnetic field is present, the situation changes dramati- σ=0 C2 cally. In particular, in the case of a purely magnetic field, which represents the external spacetime of a cosmic string all the solutions can be produced by physically accept- located on the symmetry axis r = 0 with the angular de- able cylindrical thin shells, while in the case of a purely fect ∆ϕ ≡ 4µ = 1 − C−1, where µ denotes the mass per electric field, no such shells are found for any choice of unit length of the string. Note that the metric (4.4) also the free parameters involved in the solutions. becomes singular at r = rg ≡ π/(2β) for Λ > 0. As shown below, this singularity is a curvature one, which makes the interpretation of the solutions as representing IV. LC SOLUTIONS COUPLED WITH COSMOLOGICAL CONSTANT a cosmic string spacetime coupled with a positive cos- mological constant unclear. However, in the case Λ < 0, the trigonometric functions become hyperbolic ones, and A. Linet-Tian (LT) solutions this singularity does not exist any more. So, in the latter the corresponding spacetime indeed can be interpreted The Linet-Tian (LT) solutions can be cast in the form as representing a cosmic string embedded in an asymp- [124, 221, 304], totically anti-de Sitter spacetime. To study the LT solutions further, let us consider the 2 2/3 h −2(4σ2−8σ+1)/(3Σ) 2 ds = Q (r) P (r) dt cases Λ > 0 and Λ < 0, separately.

2 −P (r)2(8σ −4σ−1)/(3Σ)dz2 2 i −C−2P (r)−4(2σ +2σ−1)/(3Σ)dϕ2 − dr2, (4.1) B. Main properties of LT solutions with Λ > 0 where Σ is given by (2.9), As r → 0, we find that P (r) → r and Q(r) → r, so that (4.1) approaches to the vacuum LC solutions. Then, 1 1 P (r) ≡ tan (βr) ,Q(r) ≡ sin (2βr) , (4.2) the singularity behavior of the LT solutions at r = 0 are β 2β similar to that of the LC solutions. In particular, all √ the solutions have a curvature singularity, except for the β ≡ 3Λ/2 and C is an arbitrary real integration con- cases σ = 0 and σ = 1/2 [321]. stant. Without loss of the generality, we assumed C > 0. When σ = 1/2, similarly to the case σ = 0, the cor- It is easy to show that the above solutions reduce to the responding solution also possesses an additional Killing −1 LC solutions given by (2.10), when Λ = 0, so in the fol- vector, ξ(1/2) = C ϕ∂z − Cz∂ϕ. This corresponds to lowing we do not consider this case any more. When a rotation in the (z, ϕ)-plane. Since in the latter case Λ < 0, the trigonometric functions become hyperbolic, the curvatures of the two spacelike surfaces t, r = Con- but the resulting solutions are still well-defined and real. stant are identically zero, it is difficult to consider this All the solutions with Λ 6= 0 are Petrov type I according spacetime as having cylindrical symmetry. Instead, one to the classifications given in [294], except for the par- may extend the ϕ coordinate from the range [0, 2π] to the ticular cases, σ = −1/2, 0, 1/4, 1/2, 1, which are all range (−∞, +∞), so the resulted spacetime has a plane Petrov type D, as shown explicitly in [124]. symmetry, as in the vacuum case [124, 332]. In addition, it was also found that the solutions with However, the metric (4.1) shows that the solutions usu- σ = 1/4 + τ can be obtained from the ones with σ = ally are also singular on the hypersurface rg ≡ π/(2β). 1/4 − τ, together with the replacement [124], To see the nature of the singularities, let us first note that Q(r) → r − r and P (r) → (r − r )−1, as r → r . −1 0 0 g g g t = iC ϕ , ϕ = iCt . (4.3) Then, we find that the metric (4.1) takes the asymptotic form, Yet, the solutions also remain the same if we replace σ 2 2 with 1/(4σ) and exchange the z- and ϕ-coordinates [347], ds2 ' R4(4σ −5σ+1)/(3Σ)dt2 − R−4(2σ −σ−1)/(3Σ)dz2 similarly to the vacuum case. 2 −C−2R2(8σ +2σ−1)/(3Σ)dϕ2 − dr2, (4.6) When σ = 0, (4.1) reduces to

2 as r → rg, where R ≡ r − rg. Then, the corresponding 2 4/3 2 2
sin (βr) 2 2 ds = cos (βr) dt − dz − dϕ −drKretschmann. scalar is given by, σ=0 β2C2 cos2/3 (β r) (4.4) 64(σ − 1)2(2σ + 1)2(4σ − 1)2 R ≡ RαβγδR ' , (r ' r ). Clearly, in addition to the usual three Killing vectors αβγδ 27Σ3R4 g ξ(t) = ∂t, ξ(z) = ∂z, ξ(ϕ) = ∂ϕ, the solutions now admit (4.7) one more Killing vector ξ(0) = t∂z + z∂t, which corre- Therefore, the spacetime is always singular at r = rg, ex- sponds to a Lorentz boost in the (t, z)-plane. As r → 0, cept for the cases σ = −1/2, 1/4, 1. It can be shown that 15 all 14 independent scalars, built from the Riemann tensor in 4-dimensional Riemann spacetimes [88], have the same properties. Therefore, in the cases σ = −1/2, 1/4, 1 the singularities on the hypersurface r = rg, are coordi- nate ones, and to have geodesically complete spacetimes, the solutions need to be extended beyond this surface. Note that, in addition to the usually three Killing vec- tors, they also have a fourth Killing vector, given by −1 −1 ξ(−1/2) = C ϕ∂z − Cz∂ϕ, ξ(1/4) = C ϕ∂t + Ct∂ϕ, and ξ(1) = z∂t + t∂z, respectively. Using the same arguments as those given for the solution with σ = 1/2, the solution with σ = −1/2 can be also considered as representing plane symmetry, and the coordinate ϕ can be extended to the whole axis −∞ < ϕ < +∞. Therefore, except for the particular cases, σ = −1/2, 0, 1/4, 1/2, 1, the solutions are singular at both r = 0 and r = rg. The physics of these singularities is unclear, and FIG. 2. Carter-Penrose diagram for the LT solutions with some of them (if not all) can be considered as representing 1 σ = 2 and Λ > 0 in the region X ≤ 0, given by (4.11) and matter sources [47, 50]. (4.12). The horizontal line (X = 0) represent a spacetime When Λ > 0, in addition to (4.3), extra symmetry singularity [124]. exists [168],

4σ2−8σ+1 π  4  3Σ extended spacetime two new regions X > 1 and X < r = √ − r0, t = t0, 0 are obtained. However, the curvature singularity at 3Λ 3Λ X = 0 divides the whole spacetime into two unconnected 2 2 2σ +2σ−1 − 8σ −4σ−1  4  3Σ  4  3Σ regions X ≥ 0 and X ≤ 0. ϕ = z0, z = ϕ0,(4.8) In the region X ≤ 0, the function f(X) is always neg- 3Λ 3Λ ative and the X coordinate is timelike, while T is space- where σ ≡ (1 − 4σ0)/[4(1 − σ0)]. Because of this symme- like. Then, the spacetime is essentially time-dependent try, it was argued that z and ϕ should have the same and the singularity at X = 0 is spacelike. As X → −∞, character, and in particular should be all periodic [168]. the metric becomes asymptotically de Sitter, Due to the particular symmetry and singular behavior 2 ˜2 4βT/˜ 3 2 2 2 of the cases σ = ±1/2, in the following let us consider ds 1 ' dT −e dX +dY +dZ , (X → −∞), σ= 2 them separately. For more details, interested readers are (4.12) referred to [124, 332]. where T = e2βT/˜ 3, and X,Y,Z had been rescaled. The corresponding Carter-Penrose diagram is given by Fig.2. In the region X ≥ 0, the function f(X) is greater than 1. Global structure of the solution with σ = 1 2 zero for 0 ≤ X < 1 and less than zero for X > 1; that is, X is spacelike in the region 0 ≤ X < 1 and timelike in the When σ = 1/2, making the coordinate transformations region X > 1. On the hypersurface X = 1, it becomes [124], null, which represents a horizon. Since the spacetime sin- gularity at X = 0 now is timelike, the horizon is actually 2 2β 2β T = t, X = cos2/3(βr),Y = ϕ, Z = z, (4.9) a Cauchy horizon [177]. As X → +∞, the spacetime is 3 3C 3 also asymptotically de Sitter and approaches the same the corresponding metric takes the form, form as that given by (4.12). The corresponding Carter- Penrose diagram is given by Fig.3. " 2 # 2 9 2 dX 2 2 2 ds σ= 1 = 2 f(X)dT − − X dY + dZ , 2 4β f(X) 1 2. Global structure of the solution with σ = − 2 (4.10) where When σ = −1/2, the spacetime is singular at r = 0, 1 and is free of curvature singularity at r = rg. Thus, to f(x) ≡ − X2. (4.11) X have a geodesically complete spacetime, we need to ex- tend the solution beyond the hypersurface r = rg. To From (4.9) we can see that the region 0 ≤ r ≤ rg is make such an extension, we can introduce a new coor- 2/3 mapped to the region 0 ≤ X ≤ 1, and the point r = rg is dinate X by X = sin (βr), and then rescale the co- mapped to the point X = 0, which is singular. Extending ordinates t, z, and ϕ, we will find that the correspond- X to the range X ∈ (−∞, +∞), we find that in the ing metric takes the same form as that given by (4.12). 16

r ∈ (0, +∞). In particular, as r → ∞, we find that a b " 2 # 2 2 4βr/3 2 2 dϕ 2 ds ' C e dt − dz − − dr , (4.13) II r→∞ 0 C2

where C0 is a real constant. As first noted in [124], this is the anti-de Sitter spacetime, but written in the horo- I’ I spherical coordinates [114]. Therefore, all the LT solu- tions with a negative cosmological constant are asymp- x = 0

x = 0 totically anti-de Sitter. II’

1. Global structure of the solution with σ = 0 c d When σ = 0, the solutions are given by (4.4), i.e.,

ds2 = cosh4/3 (|β|r) dt2 − dz2
FIG. 3. Carter-Penrose diagram for the LT solutions with σ=0 σ = 1 and Λ > 0 in the region X ≥ 0, given by (4.11) 2 2 sinh (|β|r) 2 2 and (4.12). The vertical lines (X = 0) represent spacetime − dϕ − dr ,(4.14) 2 2 2/3 singularities, while the ones ad and bc (X = 1) represent |β| C cosh (|β| r) Cauchy horizons. As X → ∞ (represented by the curves ab which can be considered as representing a cosmic string and cd), the spacetime is asymptotically de Sitter [332]. located on the symmetry axis r = 0 embedded in the anti- de Sitter spacetime with its mass per unit length given by µ = 1 − C−1
/4. Clearly, when C = 1, the cosmic Thus, the solution with σ = 1/2 is actually the same as string disappears, and a regular axis is obtained. Then, that with σ = −1/2. As first noticed in [124], this is the whole spacetime is free of any kind of singularities, not expected! As a matter of fact, in the limit Λ → 0, including a conic one, and has a well-defined symmetry the solution with σ = 1/2 approaches the Rindler space axis at r = 0. In fact, the spacetime is locally flat at [268], which represents a uniform gravitational field and r = 0. So, we have an anti-de Sitter spacetime with is free of any kind of spacetime curvature singularities, cylindrical symmetry, and this solution is different from while the one with σ = −1/2 is the static Taub solution the anti-de Sitter spacetime, with plane symmetry [299], and is singular on the hyper- surface X = 0. The total mass of the Taub spacetime is  |Λ|r2   |Λ|r2 −1 ds2 = 1 + dt2 − 1 + dr2 negative, while the one of Rindler is not [125]. However, AdS 3 3 the presence of the cosmological constant makes up these −r2 dθ2 + sin φ2dφ2
, (4.15) differences and turns the two spacetimes to be identical! We would also like to call the attention to the comments as first noted by Bonnor in [49], and dubbed the on conformally flat sources given at the end of section non-uniform (NoU) AdS universe. The difference be- 2.5. tween these metrics can be seen clearly by calculating It is also interesting to note that none of these two so- their Kretschmann scalars, lutions represents black holes. In fact, there is a general αβγδ theorem for spacetimes with a positive cosmological con- K ≡ R Rαβγδ ( 2 h i stant and two spacelike Killing vectors [319], which shows 4Λ 3 − 2 tanh2(|β|r) + tanh4(|β|r) , NoU-AdS, clearly that in such cases no black holes exist, as long as = 3 8Λ2 matter fields satisfy the dominant energy condition [177]. 3 , AdS. (4.16) In fact, the AdS spacetime given by (4.15) is conformally C. Main properties of LT solutions with Λ < 0  AdS  flat Cµναβ = 0 , while the one given by (4.13) is not  Non-U-AdS  Near r = 0, the solutions with Λ < 0 have the same Cµναβ 6= 0 , where Cµναβ denotes the Weyl ten- asymptotic properties as those with Λ > 0, and all ap- sor. proach to the LC vacuum solutions. Therefore, all the LT solutions with Λ < 0 are singular at r = 0, except 1 for the two particular cases, σ = 0, 1/2. On the other 2. Global structure of the solution with σ = 2 hand, from (4.2) we can see that both P (r) and Q(r) are monotonically increasing functions of r now, so the LT When σ = 1/2, the corresponding solution also pos- −1 solutions with Λ < 0 are well-defined in the whole range sesses an additional Killing vector, ξ(1/2) = C ϕ∂z − 17 a Cz∂ϕ, which represents a rotation in the (z, ϕ)-plane, and implies that this plane is flat. Thus, one can also consider t = Const. that the solution has a plane symmetry by simply extend- ing ϕ to the range ϕ ∈ (−∞, +∞). Then, we should also need to extend r to the full range r ∈ (−∞, +∞). How- ever, the metric has a coordinate singularity at r = 0, and extension beyond it is needed, in order to obtain a geodesically complete spacetime. The extension is quite similar to the case σ = 1/2, Λ > 0. In fact, setting c 2 2|β| 2|β| T = t, X = cosh2/3(|β|r),Y = ϕ, Z = z, 3 3C 3 x = 0 (4.17) the corresponding metric takes the form, x = Const. " 2 # 2 9 2 dX 2 2 2 ds σ= 1 = 2 −f(X)dT + −X dY +dZ , 2 4β f(X) (4.18) where f(X) is given by (4.11). Comparing it with (4.10) b we can see that we can obtain one from the other by sim- ply replacing f(X) with −f(X). From the expression of FIG. 4. Carter-Penrose diagram for the LT solutions with σ = 1 and Λ < 0 in the region X ≤ 0. The vertical double line X we can see that the region 0 ≤ r < ∞ is mapped to 2 (X = 0) represents a spacetime singularity, which is timelike. the region 1 ≤ X < ∞, and region X < 1 is an extended The curve bca corresponds to X = −∞, which is also timelike region. After the extension, a spacetime curvature singu- [332]. larity appears at X = 0, which divides the whole X-axis into two parts X ≤ 0 and X ≥ 0. It can be shown that, unlike the case Λ > 0, now the spacetime is static in the a x = 0 b region X ≤ 0 and the curvature singularity at X = 0 is timelike and naked. As X → −∞, the spacetime is asymptotically anti-de Sitter, II

2 9  2 ˜ 2 2 2 ds σ= 1 ' dT − dX − dY − dZ , (4.19) 2 4|β|2X˜ 2 I’ I where X˜ = X−1. The corresponding Carter-Penrose di- agram is given by Fig.4. II’ In the region X ≥ 0, the spacetime singularity at X = 0 becomes spacelike. Except for this curvature sin- c d gularity, there is a coordinate one located at X = 1. x = 0 This coordinate singularity actually represents an event horizon. As shown in the last section, the spacetime is FIG. 5. Carter-Penrose diagram for the LT solutions with 1 asymptotically anti-de Sitter as X → +∞. The corre- σ = 2 and Λ < 0 in the region X ≥ 0. The horizontal sponding Carter-Penrose diagram is given by Fig.5. This lines (X = 0) represent spacetime singularities, while the ones is a black hole solution with plane symmetry and vanish- X = 1, represented by the straight lines ad and bc, correspond ing electromagnetic charge first found in [86], and then to event horizons. As X → ∞ (represented by the curves ac studied in detail in [332]. It should be noted that this and bd), the spacetime is asymptotically anti-de Sitter [332]. does not contradict the theorem given in [319], as now we have Λ < 0. this question, let us first introduce a new coordinate X = − sinh2/3(|β|r), and then rescale the other three, we 1 3. Global structure of the solution with σ = − 2 find that the metric takes the same form as that given by (4.18). From the expression for X we can see that Note that in this case an additional Killing vector the region 0 ≤ r < +∞ now is mapped to the region −1 ξ(−1/2) = C ϕ∂z − Cz∂ϕ also exists, and one may −∞ < X ≤ 0, while the region −∞ < r ≤ 0, is mapped consider that the spacetime also has a plane symmetry. to the region 0 ≤ X < +∞. In the region 0 ≤ r < ∞, the Hence, the range of r should be taken as r ∈ (−∞, +∞). solution represents a static spacetime with a naked singu- Then one may ask, what is the physical interpretation larity located at X = 0. The spacetime is asymptotically of the spacetime in the region r ≤ 0? To answer anti-de Sitter, and the corresponding Carter-Penrose di- 18 agram is given by Fig.4. In the region −∞ < r ≤ 0, the 1. Planar geodesics z˙ = 0 solution represents a black hole with plane symmetry, and the corresponding carter-Penrose diagram is given In this case, let us consider the null and timelike by Fig.5. Thus, similarly to the case with Λ > 0, now geodesics, separately. the solution with σ = 1/2 and the one with σ = −1/2 actually describes the same spacetime. • Null geodesics  = 0: For the null geodesics, 2 2 −2(1−4σ)/Σ when σ ≤ 1/4, if E > V∞ = (CLz) P , then outgoing particles escape to infinity for Λ < 0, while D. Geodesics of LT solutions they reach the outer singularity for Λ > 0. Incoming geodesics hit the axis for σ = 1/4, but for σ < 1/4, Geodesics describe the orbits of light rays and material they always bounce at a minimum distance r from the p p min particles, and an important physical aspect is the influ- axis, where ( 3|Λ|/2)rLCmin = ( 3|Λ|rmin/2)T and Σ/(1−4σ) ence of the cosmological constant on the stability of the rLCmin = (CLz/E) . geodesics’ orbits in an otherwise vacuum spacetime with When√ σ > 1/4, using V (r → 0) → 0 and V (r → Λ = 0. In this context spherical spacetime with Λ has π/ 3Λ) → +∞, for Λ > 0, and V (r → ∞) → V∞ = been studied in [215]. Here we summarize, in an unified p −2(1−4σ)/Σ 2 (2/ 3|Λ|) (CLz) and V (r → 0) → 0, for way, the results about the impact of Λ on the dynamics 2 Λ < 0 (and E ≤ V∞), one concludes that incoming par- of the geodesics in the LT spacetime. We recall the LT ticles approach the axis with infinite speed. In turn, out- metric (4.1) with P and Q given by going particles move with decreasing speed and negative p acceleration, attaining a maximum distance r from 2 1 3|Λ| p p max P = p T (R),Q = p S(2R),R = r, the axis, where ( 3|Λ|/2)rLCmax = ( 3|Λ|rmax/2)T 3|Λ| 3|Λ| 2 Σ/(4σ−1) with rLCmax = (E/(CLz) being the value of (4.20) rLCmax for Λ = 0. This shows, in a particularly simple where T and S are the tangent and sine functions if Λ > way, that Λ > 0 decreases the maximum distance to the 0, or the corresponding hyperbolic functions if Λ < 0. axis reached by the null particle, while Λ < 0 increases We follow closely the results obtained for Λ > 0 in 2 it. In the latter case, now with E > V∞, outgoing null [60] and Λ < 0 in [59]. We recall too√ that the range particles escape to infinity. of the radial coordinate r is r ∈ (0, π/ 3Λ) and there Proposition 1. Consider null geodesics with Λ = 0, are two singularities for Λ > 0, and the range of the Pz = 0 and Lz 6= 0. Then, bounded geodesics are: (i) radial coordinate r is r ∈ (0, ∞) and there is only one 2 unstable, if E > V∞, with the introduction of Λ < 0; singularity for Λ < 0. If Λ = 0, one simply gets P = Q = 2 (ii) stable, if E ≤ V∞, with the introduction of Λ < 0; r. (iii) stable, with the introduction of Λ > 0. The geodesics for the LT metric are given by

−2/3 2(1−8σ+4σ2)/(3Σ) • Timelike geodesics  = 1: In this case, when σ < t˙ = EQ P , (4.21) 2 1/4, for Λ = 0, since V (r) = Q2/3P −2(1−8σ+4σ )/(3Σ) + 2 −2/3 2(1−8σ+4σ2)/(3Σ) 2 2 −2(1−4σ)/Σ r˙ = Q P [E − V (r)], (4.22) (CLz) P always has a minimum, the geodesics 2 is always confined between two nonzero radii and this is z˙ = P Q−2/3P 2(1+4σ−8σ )/(3Σ), (4.23) z also the case for Λ < 0. On the contrary, if Λ > 0, ˙ 2 −2/3 −4(1−2σ−2σ2)/(3Σ) φ = C LzQ P , (4.24) then dV/dr = 0 does not always have a solution. In fact, for sufficiently large Λ, bounded geodesics become where the dot stands for differentiation with respect to √ unbounded and converge to the r = π/ 3Λ singularity. an affine parameter λ; the constants E, P and L repre- z z When σ ≥ 1/4, by similar methods, we prove that sent, respectively, the total energy of a test particle, its outgoing geodesics reach a maximum radial value before momentum along the z axis and the angular momentum bouncing back inwards to the axial singularity and in- about the z axis; and the effective potential V (r) is given coming particles hit the axis. by Proposition 2. Consider timelike geodesics with Λ = 2/3 −2(1−8σ+4σ2)/(3Σ) 2 8σ(1−σ)/Σ V (r) = Q P + Pz P 0 and Pz = 0. Then, bounded geodesics are: (i) stable, for a sufficiently small Λ > 0; (ii) unstable, for suffi- +C2L2P −2(1−4σ)/Σ, (4.25) z ciently large Λ > 0; (iii) stable, against any value of where  = 0, 1 or −1 if the geodesics are, respectively, Λ < 0. null, timelike or spacelike. We note that in (4.22) E2 = V (r) givesr ˙2 = 0 and corresponds to circular helices along cylindrical surfaces, 2. Non-planar geodesics z˙ 6= 0 or to planar circles if Pz = 0. We also note that un- bounded geodesics for Λ = 0 remain unbounded after the introduction of any Λ 6= 0 so our aim here will focus • Null geodesics  = 0: In this case, when σ < 1/4, on stability issues for bounded geodesics, where we also for Λ = 0, the geodesics are always confined between two consider Lz 6= 0. nonzero radii and the same happens for Λ > 0, since 19

V (r) has a minimum and the equation E2 = V (r) has it is possible to define unique evolutions for arbitrary 2 two real roots. For Λ < 0, if E ≥ V∞, where V∞ = perturbations that satisfy such imposed boundary con- 2 p 8σ(1−σ)/Σ 2 p −2(1−4σ)/Σ Pz (2/ 3|Λ|) + (CLz) (2/ 3|Λ|) , a ditions. Second, the spacetime is asymptotically anti-de null particle approaches the axis with decreasing negative Sitter as x → ∞, as shown by (4.13), and this boundary is acceleration and increasing speed, until it arrives at its timelike. As a result, the hypersurfaces t = Constant are minimum distance from the axis at E2 = V (r), where it not Cauchy surfaces, and there always exist null geodesics has vanishing speed. From there on, the null particle is that remain to the future of any given constant t hyper- 2 reflected escaping to infinity. If E < V∞, the equation surface, and never intersect the hypersurface [cf. Fig.4]. E2 = V (r) has two roots and the geodesics are confined Then, the future evolution of any perturbation may be between two nonzero radii. arbitrarily modified by information coming from x = ∞. When σ ≥ 1/4, for Λ = 0, the geodesics are always To overcome the second problem, following [193], Gleiser confined between two nonzero radii. If Λ < 0 and E2 ≥ found that under certain circumstances some boundary V∞, the null geodesics become unbounded. Otherwise, conditions can lead to a well-defined evolution of appro- the geodesics reach a maximum. priate initial perturbation data. With such boundary conditions, he was able to show that unstable modes al- • Timelike geodesics  = 1: In this case, there is ways exist. Then, he argued that this should be gener- always geodesic confinement in the radial direction for ically present in the evolution of arbitrary initial data, Λ ≤ 0. If Λ > 0, the geodesics are only confined for and concluded that the LT spacetimes with a negative σ <√ 1/4, in which case V (r → 0) → +∞ and V (r → cosmological constant are generically unstable under the π/ 3Λ) → +∞. If σ ≥ 1/4, then V (r → 0) is finite, linear gravitational perturbations. therefore, incoming particles hit the axis and outgoing It should be noted that the above analysis cannot be particles reach a maximum distance before turning back applied to the cases σ = ±1/2, as in these cases the so- to the axis. We summarize the results of this subsection lutions actually represent black holes, as shown in Fig.5, below. and the boundary conditions should be now at the hori- zons X = 1, instead of the symmetry axis r = 0, while Proposition 3. Consider geodesics with P 6= 0. z the ones at r = ∞ should be similar to what Gleiser con- Then, radially bounded geodesics with Λ = 0 are: sidered. So, it would be very interesting to see if such (i)stable with the introduction of any Λ > 0; (ii) stable black holes are stable or not. with the introduction of any Λ < 0 in the case of timelike Recently, the above studies were generalized to the case geodesics; (iii) unstable, if E2 ≥ V , for Λ < 0 in the ∞ Λ > 0 [160]. However, due to the complexities of the case of null geodesics. perturbation equations, a similar conclusion for the gen- eral case is absent, but in all the cases analyzed unstable E. Instability of LT solutions modes were found, which strongly suggests that all the LT spacetimes with Λ > 0 are linearly unstable under gravitational perturbations, although the problem of de- The stability of the solutions with Λ < 0 was recently termining the time evolution of arbitrary initial data still studied in [159]. Consider the kind of perturbations, remains open. (0) gµν = gµν (x) + hµν , (4.26)

(0) F. Sources producing LT spacetimes where gµν (x) is the LT solutions, and hµν = hµν (t, r, z) with   1 and x ≡ sinh(|β|r). Since the background To study the LT solutions further, it is important to depends only on x, one finds that hµν (t, r, z) takes the form, study the possible sources that could produce such space- times. In particular, the spacetime singularities are nor- Z mally considered as some kind of indications that in re- h (t, x, z) = e−(Ωt−kz)C(Ω, k)H (Ω, k, x) dΩdk, µν µν ality they should be replaced by matter sources that pro- (4.27) duce such spacetimes [47, 50]. where C(Ω, k) is determined by the initial data. It is well- Along this vein, in the case Λ > 0 the solutions are known that the spectrum of perturbations (sensitively) usually singular at both r = 0 and r = π/(2β). Then, depends on the boundary conditions, and the choice of one can see that at least two matter sources should ex- them is subtle here with the LT solutions. First, they ist in this case, and the LT solutions are only valid in are singular at the symmetry axis r = 0, except for some between these two sources. Then, it is not clear what particular cases, as shown above. These singularities are are the physical meanings of such spacetimes, if there is timelike, and there is no general rule on how to impose any. In [168], because of the symmetry given by (4.8) it boundary conditions at such singular points. Following was argued that the coordinate z is more like an angular the studies carried out for the vacuum LC solutions [158], coordinate than an infinitely long axis. So, replacing z Gleiser found that if one imposes certain physically ac- by ψ which has the range ψ ∈ [0, 2π), the LT solutions ceptable restrictions as the boundary conditions at x = 0, with Λ > 0 and σ ∈ [0, 1/4] were matched with the static 20

Einstein universe, V. STATIC CYLINDRICAL SPACETIMES WITH PERFECT FLUIDS B2 C2 ds2 = A2dt2 − 1 cos2 ∆dψ2 − 1 sin2 ∆dϕ2 − dr2, Static perfect fluid distributions with cylindrical (as ES 1 Λ Λ (4.28) well as planar) symmetry are rather widely discussed in √ the literature. This allows for approximately considering along the hypersurface r = r1, where ∆ ≡ λ (r − r0), the fields of bodies whose shape is drastically nonspher- A1,B1,C1 and r0 are real constants and are uniquely de- ical (e.g., like rods) avoiding substantial mathematical termined by the junction conditions along r = r1. How- ever, as noted in [168], the static Einstein universe can difficulties of more general geometries. In a number of replace only one of the two singularities, and to have a papers (see [210, 231, 273] and references therein) so- singularity-free spacetime, one needs to replace the re- lutions of this kind were obtained for special choices of mained singularity, too. matter equations of state (EoS) or for some restricted forms of the metric. In [140] the problem of finding such With the static Einstein universe as an example, one solutions for an unspecified EoS was reduced to a single may also consider r as an angular coordinate, so that the second-order ordinary differential equation with two un- hypersurfaces r = 0 and r = π/(2β) are identical, and known functions, and for the most frequently used EoS, the two singularities become one. In fact, setting p = wρ, w = const, to a first-order equation integrable only numerically. In [300] the problem was solved com-   pletely for disordered radiation (w = 1/3). In a more r − r0 sin = sin χ sin θ, tan(B1ψ) = tan χ cos θ, general form the problem was solved in [64, 77] both for a (4.29) fluids with an arbitrary EoS and for p = wρ. the metric (4.28) will take the form, Later there appeared a number of studies (e.g., [33, 188, 262]) discussing global properties of static, cylindri- cally symmetric perfect fluid configurations. This actu- 2 2 2h 2 2 2 2 2
i ally already began in [77], considering systems with both dsES = dT −a dχ +sin χ dθ + sin θdΦ , (4.30) cylindrical and (pseudo)planar symmetries. In the recent years, owing to new discoveries in astron- where a is a constant (the radius of the static Ein- omy, fluids with unusual EoS attracted much attention, stein universe), and T ≡ A1t and Φ ≡ C1ϕ, with including those with negative pressures and even those T ∈ (−∞, ∞), χ ∈ [0, π], θ ∈ [0, π], and Φ ∈ [0, 2π). with w = p/ρ < −1, the so-called phantom ones. Such It is clear that (4.30) represents a static universe with cylindrically symmetric distributions were considered in R × S3 topology. One can apply a similar coordinate [68]. In this section we describe the main results of these transformation to the LT solutions with the identifica- studies mostly following [68, 77]. However, some of the tions of r = 0 and r = π/(2β). Then, one can consider relationships and observations are new. the Einstein static universe as the internal spacetime of a ball with radius, say, a1, while the spacetime outside of the spherically symmetric ball is given by the LT so- A. General consideration lutions.

In addition, in [57], the authors considered the match- We write the static, cylindrically symmetric metric in ing of an elastic matter to the LT solutions and found the form (A.1) that such a matching is possible for both Λ > 0 and Λ < 0, but it is impossible for Λ = 0. On the other hand, ds2 = e2γ dt2 − e2αdx2 − e2µdz2 − e2βdϕ2, (5.1) in [58] the matching, across cylindrical surfaces, of static cylindrically symmetric conformally flat spacetimes with and choose the harmonic radial coordinate such that a cosmological constant, satisfying regularity conditions (A.4) at the axis, to an exterior LT spacetime was studied, and it was found that for Λ ≤ 0, such matching is impossible, while the one with Λ > 0 is possible. α(x) = γ(x) + β(x) + µ(x). (5.2) In the case Λ < 0, for the black holes solutions (with Einstein’s equations (1.13) for the metric (5.1) then take σ = ±1/2) it was first shown [332] that they can be pro- the following simple form: duced by gravitational collapse of type-II fluid [177]. In [347], the junctions of the LT solutions with (σ−, Λ−,C−) 00 00 0 0 0 0 0 0 0 2α β + µ − µ β − µ γ − β γ = −κT0 e , (5.3) with the ones with (σ+, Λ+,C+) were considered with a 0 0 0 0 0 0 1 2α matter shell appearing on the matching surface. µ β + µ γ + β γ = −κT1 e , (5.4) It should be noted that the matching given in [57, 347] 00 00 0 0 0 0 0 0 2 2α µ + γ − µ β − µ γ − β γ = −κT2 e , (5.5) for Λ > 0 faces the problem that a spacetime singularity 00 00 0 0 0 0 0 0 3 2α still exists at r = π/(2β) in the LT spacetime. γ + β − µ β − µ γ − β γ = −κT3 e , (5.6) 21

ν where the prime denotes d/dx, and we take Tµ as the B. Solutions for p = wρ, the general case stress-energy tensor of a perfect fluid in its comoving ref- µ −γ erence frame, so that the 4-velocity is u = ( e , 0, 0, 0): With the EoS p = wρ, w = const, the condition (5.8) ν ν ν leads to Tµ = (p + ρ)uµu − δµp = diag(ρ, −p, −p, −p). (5.7) ν −(w+1)γ(x)/w −(w+1)γ(x)/w The conservation law ∇µTµ = 0 implies ρ = ρ0 e , p = wρ0 e , (5.12) p0 + γ0(p + ρ) = 0. (5.8) for w 6= 0, where ρ0 = const > 0. At w = 0 the tensor For a cylindrically symmetric matter distribution, a (5.7) would correspond to dust which cannot be in static natural boundary condition is regularity on the symme- equilibrium under pure gravitational forces, so we assume 2 2β try axis where r ≡ e = −g33 → 0 (the coordinate w 6= 0. circles shrink to a point). In terms of the metric (5.1), From (5.5), (5.6) we have, as before, (5.10), and we this condition reads take it in the form γ = γ + O(r2), µ = µ + O(r2), e−2αr02 = 1 + O(r2), 0 0 µ(x) = β(x) + a x − ln r , (5.13) (5.9) 1 0 where the last condition is the requirement of a correct where the arbitrary constant r0 with the dimension of circumference to radius ratio (2π) for small coordinate length will determine a length scale. Another linear com- circles near the axis. As to the other end of the range bination of (5.3)–(5.6) then leads to of x, there can be a matching condition with some ex- ternal solution or some requirement to the behavior of (1 + 3w)β00 + (1 − w)γ00 = 0. (5.14) variables as r → ∞. It is necessary to note that if the external solution is the vacuum (LC) one, then match- The simpler cases 1 + 3w = 0 and 1 − w = 0 will be ing on the boundary requires p = 0, which restricts the considered separately. Assuming that these quantities choice of possible EoS. On the other hand, one could re- are nonzero, we can write quire asymptotic flatness at infinity, but this is hard to achieve since the total mass of a cylindrically symmetric (w − 1) β(x) = γ(x) + b x + b + ln r , b , b = const, configuration is, in general, infinite due to its being in- (3w + 1) 1 2 0 1 2 finitely extended in the z direction. Even the LC metric is (5.15) asymptotically flat only if it is completely flat. If cylindri- and put b2 = 0 by properly choosing the zero point of cal symmetry is approximate and valid only in a certain x. With (5.13), (5.15) and (5.2), all metric functions are neighborhood of an extended body, then the asymptotic expressed in terms of γ(x). The latter can be determined flatness requirement applies to regions far outside this using the first-order equation (5.4) after proper substitu- neighborhood. tions. Before considering particular solutions, let us make Let us denote some general observations. A If the EoS p = p(ρ) is not unspecified, there are four A = 7w2 − 6w − 1 B = , (5.16) Einstein equations for five unknowns β, γ, µ, ρ and p, or, w(3w + 1) equivalently, three independent Einstein equations and 2η(x) = Bγ(x) + 2ax. (5.17) the hydrostatic equilibrium condition (5.8), and one un- known function may be chosen arbitrarily. This makes Then (5.4) leads to the equation possible obtaining a formal general solution with an ar- bitrary function.  2  02 A 2 2η(x) 4w 2 Indeed, the difference of (5.5) and (5.6) yields 4η (x) = wρ0r e + (a1+2b1) −b1(a1+b1) , w2 κ 0 A 00 00 β = µ ⇒ µ = β + a1x + a2, a1, a2 = const. (5.10) (5.18) which is easily solved by quadratures to finally yield all Furthermore, the difference of (5.4) and (5.5) can be writ- unknown functions. Thus under the conditions w 6= ten as 0, 1, −1/3, −1/7 (A 6= 0,B is finite) the solution is com- 00 0 02 02 η − 2a1η − 2η + 2γ = 0, η ≡ β + γ, (5.11) pletely obtained analytically, though its particular form depends on the values of w and the integration constants and γ(x) is easily found by integration if we take η(x) a , b . as an arbitrary function. Thus the metric is completely 1 1 One can show that perfect fluid distributions with determined. The quantities ρ(x) and p(x) are then ob- p = wρ cannot behave as isolated systems in space, such tained directly from Einstein’s equations, and the EoS is as threads, tubes or cosmic strings. To verify that, let thus found in a parametric form. us write the corresponding conditions at spatial infinity, On a regular axis it should be finite µ while β → −∞, x → x∞, such that hence a1 6= 0, and the axis is located at x = +∞ if a1 > 0 and at x = −∞ if a < 0. Examples of solutions with a 1 β(x) regular axis will be presented below. r(x) ≡ e = ∞. (5.19) x=x∞ 22

The metric (5.1) has a flat or string spatial infinity as are x → x∞ if (w − 1) (w − 1) µ(x) = γ(x), β(x) = γ(x) + ax + ln r0, (3w + 1) (3w + 1) 2β−2α 02 |γ| < ∞, |µ| < ∞, e β → 1 − ξ, (5.20) x→x∞ (5w − 1) α(x) = γ(x) + ax + ln r0. (5.25) where ξ = const. If ξ = 0, we have asymptotic flatness, (3w + 1) at 0 < ξ < 1 there is a deficit of the angle ϕ (the surfaces The metric takes a simpler form in terms of the new z = const asymptotically behave as cones rather than radial coordinate y = eax, such that y = 0 is the axis and planes, it is a feature of cosmic string behavior), while at y = ∞ is spatial infinity. Rescaling the time coordinate 2γ ξ < 0 there is an angular excess: this is a cosmic string so that g00 = e = 1 on the axis, we obtain with negative linear density. 2 2
s0 2 (w−1)/(2w)h 2
s1 2 2 In our solutions, if w 6= −1/3, from (5.15) it follows ds = 1 + y dt − a 1 + y r0dy that γ can be finite while β → ∞ only if b1x → ∞. s i −1 + y2
2 (dz2 + r2y2dϕ2) , Assume without loss of the generality b1 > 0, so x∞ = ∞. 0 Then (5.13) shows that |µ| < ∞ requires a1 = −b1, and 4 4w(3w + 1) 4w(5w − 1) as x → ∞, η(x) ≈ b x. Then the left-hand side of (5.18) s0 = − = − , s1 = − , 1 B (7w + 1)(w − 1) (7w + 1)(w − 1) tends to a constant while the right-hand side infinitely 4w grows. This contradiction shows that our assumption s2 = − , (5.26) was wrong, and a flat or string-like asymptotic behavior 7w + 1 of our solutions is impossible. with a given in (5.24), and the density ρ and pressure p = wρ are determined by (5.12). The solution is valid for all w ∈ (0, 1) and w < −1/7. Its properties strongly depend 1. Solutions with a regular axis on w, so in what follows we discuss this dependence as well as some particular values of w. One should notice that (5.24) is meaningless if w2 + Let us find explicitly, without fixing w, perfect fluid √ configurations with a regular symmetry axis. On such 6w + 1 = 0 ⇒ w = −3 ± 8, so the constant a can- an axis we should simultaneously have β → −∞ and not be determined and remains arbitrary. It means that for these values of w (approximately −0.16 and −5.84) finite µ and γ, which is only possible if b1x → −∞ and a + b = 0. Assuming without loss of generality b = the axis with finite γ and µ cannot be regular and is 1 1 1 inevitably a conical-type singularity. −a1 = a > 0 (so that the axis corresponds to x → −∞), √ we can rewrite (5.18) in the form If w ∈ (0, 1), then both the circular radius r = −g33 and the temporal metric component g00(y) are increasing functions, so the gravitational field attracts test particles ρAr2 η02 = a2 − K e2η,K := −κ 0 . (5.21) towards the axis, while both ρ and p are positive and 4w decrease with growing r. However, with negative pressures (w < 0) things are Let us further assume that the fluid density is positive, ax hence ρ > 0, then K > 0 for the values of w much more diverse at large x (or y = e ). Let us present 0 for different w the large x behavior of the metric coeffi- 2γ 0 < w < 1, or w < −1/7. (5.22) cients g00 = e (describing an attractive or repulsive nature of gravity at large x) and r2 ≡ e2β (squared cir- cular radius), the fluid density, and the convergence or Then (5.21) is integrated to give divergence of some integral characteristics of the solution: R α a −2/B 3w+1 L = e dx (the total length along the radial direction), eη(x) = √ , eγ(x) = 1 + e2ax
a 4w , K cosh(ax) V = 2π R eα+µ+βdx (the total volume per unit length (5.23) along the symmetry axis), and where we have suppressed the emerging integration con- R α+µ+β stant by choosing the zero point of x. E = 2π ρ(x) e dx (the total fluid energy per unit The next step is to apply the third regularity condi- length along the symmetry axis). tion (5.9), taking into account that γ0(−∞) = 0 and Note that, up to a positive constant factor, β(−∞) = a. We obtain  2(1 + 7w2)  eα+µ+β ∼ exp ax , 4w2 4w2 (1 − w)(1 + 7w) 2   2  2  2 4a 2 K w +6w+1 −A ρ r w +6w+1 = aA/(4w ) ⇒ a = = κ 0 0 ,  2(1 + 8w − w2)  K 4 16w ρ eα+µ+β ∼ exp − ax . (5.27) (1 − w)(1 + 7w) (5.24) where K is defined in (5.21). This completes the solution. The large x behavior of these parameters is presented So γ(x) is given in (5.21), while the other metric functions in the Table II. 23

From the third axis regularity condition (5.9) it fol- TABLE I. The large x behavior for different w, where F lows a = 1, and the solution with a regular axis even means finite. more simplifies. It is conveniently written in terms of the coordinate r = r ex: 2γ 2 0 Parameters: e r ρ L V E 2 ds2 = eκρ0r (dt2 − dr2) − dz2 − r2dϕ2, (5.31) 0 < w < 1 ∞ ∞ 0 ∞ ∞ F 2 − ρ1r ρ = p = ρ0 e κ . (5.32) 1 1 − 3 < w < − 7 ∞ 0 ∞ F F F Consider the fluid energy distribution in stiff matter for 0 the solution (5.25), (5.26). The density is T0 = ρ = ρ0 1 at r = 0, and it quickly decays as r → ∞, being thus −1 < w < − 3 0 ∞ 0 F F F localized near the axis. Let us find the energy per unit w = −1 0 ∞ F F F F length along the z axis: ∞ Z p 2π 1 E = 2π T 0 −3g dr = = , (5.33) w < −1 0 ∞ ∞ F F F 0 0 κ 4G where 3g = − e2α+2β+2µ is the determinant of the spatial metric, and G is the Newtonian gravitational constant. Thus the total energy per unit length along z does not We see that in the range −1/3 < w < −1/7 we obtain depend on ρ and is of the order of Planck linear density. a closed type geometry with one regular axis and an- 0 other singular one. All solutions with negative pressure (w < 0) describe geometries with a finite volume of space 2. Disordered radiation, w = 1/3 per unit length along the axis, and in all cases under con- sideration the corresponding total fluid energy is finite. In the phantom range w < −1 the density grows with It is a special case of the above general solution: in growing radius, much like the phantom energy behavior particular, A = −20/9, B = −10/3. From (5.13)–(5.18) in cosmology. we find 1 1 β(x) = − γ(x) + b x + ln r , µ(x) = − γ(x) + (a + b )x, 3 1 0 3 1 1 C. Solutions for particular values of w = p/ρ 1 α(x) = γ(x) + ax + ln r , a = a + 2b (5.34) 3 0 1 1 1. Stiff matter, w = 1 where eγ(x) is found using (5.18) as:

 2 3/10 In this case, A = (7w + 1)(w − 1) = 0, and the above ρ0r √ 3  eγ(x) = κ 0 cosh2( 5hx) exp ax , scheme does not work, but (5.12)–(5.14) are valid and 3h2 5 lead to r a2 h ≡ + b1(a1 + b1). (5.35) β(x) = b1x + ln r0, µ(x) = (a1 + b1)x, 5

α(x) = γ(x) + ax + ln r0, a := a1 + 2b1, (5.28) A regular axis at x = −∞ is obtained with a1 + b1 = 0, a = b1 > 0, and the third regularity condition (5.9) where, as before, r0 = const specifies the length scale, requires a special value of a. Indeed, as x → −∞, 0 while from (5.4) we get the equation aγ (x)+(a1+b1)b1 = 2ax  2 2/5 κρ0 e for γ(x), whence (provided a 6= 0) −2α 02 2β−2α 02 12a e r = e β → 2 = 1, 5 ρ0r κρ0 2ax b1(a1 + b1) κ 0 γ(x) = e − x, α(x) = γ(x) + ax + ln r0. 2a2 a 5 ρ r2 1/7 ⇒ a = κ 0 0 . (5.36) (5.29) 12 If a = 0 (that is, a1 = −2b1), then (5.4) reduces to a In terms of y = eax, so that 0 < y < ∞, we obtain the 2 2 relation between constants, −2b1 = κρ0r0, incompatible metric (5.26) with with ρ > 0, not to be considered. β If b1 > 0, a symmetry axis ( e → 0) corresponds to s0 = 6/5, s1 = 2/5, s2 = −2/5, (5.37) x → −∞. As in the general case, to have finite γ(x) and so that µ(x) as x → −∞, we need a1 + b1 = 0, so b1 = a > 0, 2 2 6/5 2 −1 h 2 2/5 2 2 and ds = (1 + y ) dt − a (1 + y ) r0dy

µ(x) ≡ 0, β(x) = ax + ln r0, 2 −2/5 2 2 2 2 i −(1 + y ) (dz + r0y dϕ ) , κρ0 2ax κρ0 2ax γ(x) = e , α(x) = e + ax + ln r . (5.30) 2 −12/5 2a2 2a2 0 ρ = 3p = ρ0(1 + y ) . (5.38) 24

This solution coincides with the one presented in [64, 77, 4. A cosmological constant, w = −1 300]. The fluid energy density ρ is finite on the regular axis In this case, ρ = −p = Λ/κ = const. It is again a and vanishes as y → +∞. The total fluid energy per unit special case of the general solution, such that A = 12, coordinate length ∆z along the z axis is finite: B = 6, and η = 3γ + ax. We have 2πr2 Z ∞ 0 2 −1/5 β(x) = γ(x) + b x + ln r , µ(x) = γ(x) + (a + b )x, E = 3 ρ(1 + y ) ydy 1 0 1 1 a 0 2 Z ∞ 2 α(x) = 3γ(x) + ax + ln r , a = a + 2b . (5.44) 2πρ0r0 2 −13/5 5πρ0r0 0 1 1 = 3 (1 + y ) y dy = 3 .(5.39) a 0 8a The function γ(x) is found from (5.18) which now reads µ 2 −1/5 However, since e = (1+y ) vanishes at large y, this 02 2 2η 2 2 finite energy belongs to a vanishing asymptotic length η = −3Λr0 e + a1 + 3b1 + 3a1b1. (5.45) eµ∆z, and it looks meaningless to speak of its localiza- Its solution takes different forms depending on the sign tion. of Λ. This class of solutions is already discussed in detail in other coordinates in section 4, so let us here focus 3. A gas of cosmic strings, w = −1/3 on Λ > 0 and the solution with a regular axis, which corresponds to a1 + b1 = 0 and (assuming b1 = a > 0) x → −∞. We have It is a special case excluded in (5.15)–(5.18), but we √ can use (5.12)–(5.14) and write (choosing proper scales 4 a a = , eγ(x) = eµ(x) = , along the z and t axes) 3Λr2 (1 + e2ax)1/3 0 √ µ = β + a x − ln r , γ = b x, a , b , r = const, (5.40) ar eax a3/2r eax 1 0 1 1 1 0 eβ(x) = 0 , eα(x) = 0 , (5.46) (1 + e2ax)1/3 1 + e2ax where, as before, r0 is an arbitrary length scale. Substi- tuting this into (5.4), we obtain an equation for β(x): In terms of the coordinate y = eax, ranging from zero to 02 0 κρ0 4β+2ax infinity, the metric reads β + aβ + a1b1 = − 2 e , (5.41) 3r0 " 4 dt2 r2dy2 preserving the previous notation a = a + 2b . Already ds2 = − 0 1 1 3Λr2 (1 + y2)2/3 (1 + y2)2 from (5.36) we can conclude that a regular axis is im- 0 # possible in this system. Indeed, since on such an axis 1 |γ| < ∞, it should correspond to finite x. But then it is − (dz2 + r2y2dϕ2) . (5.47) (1 + y2)2/3 0 impossible to have simultaneously |µ| < ∞ and β → −∞. Equation (5.37) can be rewritten as One could also easily describe examples of solutions 4 ρ with phantom matter: w < −1, but this will hardly add η02 = a2 +4b2 − κ 0 e2η(x), η(x) := 2β(x)+ax. (5.42) 1 1 2 anything significant to the description of the more general 3r0 solution (5.23)–(5.26). It is easily solved, leading to s 3r2 h e−ax q 2β(x) 0 2 2 D. Some further results e = , h = a1 + 4b1. (5.43) 4κρ0 cosh(hx) This completes the solution. The range of x is x ∈ R. It In this subsection we mention some results on perfect is easy to verify that the nature of the geometry described fluids in static, cylindrically symmetric spacetimes, re- depends on the product a1b1: lated to the content of the present section but go beyond (a) If a1b1 > 0, then h < |a|, hence there is a (singular) it. axis (eβ → 0) on one end of the range of x and a spatial Much effort has been devoted to finding the conditions infinity (eβ → ∞) on the other. of matching cylindrical perfect fluid distributions with (b) If a1b1 < 0, then h > |a|, hence there are two external LC spacetime and an analysis of global prop- singular axes at x → ±∞. erties of the obtained configuration. This problem was (c) If a1b1 = 0, then h = |a| > 0, hence there is a already discussed in section 2, but here we find it appro- singular axis on one end and a finite radius on the other, priate to add some discussion. r = eβ → const, where, however, either µ(x) or γ(x) Thus, in particular, in [262] it is shown that matching tends to infinity, depending on which of the constants is possible for −1/2 < σ < 1/2, (σ is the LC solution a1, b1 is nonzero and which sign it has. parameter) but σ < 0 corresponds to a negative fluid Thus a gas of disordered cosmic strings can create dif- density. Global properties of cylindrical spacetimes are ferent geometries depending on the parameter values, discussed in [33], where analytical and numerical solu- which actually express different boundary conditions. tions have been obtained for an incompressible fluid. The However, none of these geometries is nonsingular. existence and uniqueness of global solutions are shown for 25 rather general fluid EoS admitting nonvanishing density VI. LEWIS VACUUM SPACETIME at zero pressure. In [188], it is shown that conformally flat internal solutions cannot be matched to an LC exte- A. Stationary cylindrical vacuum spacetime rior. An interesting result is presented in [144]: assuming The extension of the LC static cylindrically symmet- that there is a regular axis, that the fluid has a nonnega- ric vacuum spacetime to a stationary cylindrically sym- tive density and that the integral R dp/(ρ+p) from zero to metric vacuum spacetime was obtained independently by the pressure value on the axis converges, it is proven that Lanczos in 1924 [213] and Lewis in 1932 [219]. We con- at some finite radius the pressure is zero, hence the fluid sider the spacetime described by the cylindrically sym- distribution is finite in the radial direction and can be metric stationary metric matched with the LC vacuum solution. It generalizes a 2 2 µ 2 2 2 theorem from [33]. Note that for fluid distributions with ds = fdt − 2kdtdφ − e (dr + dz ) − ldφ , (6.1) w > 0 described above this integral diverges, and this where f, k, µ and l are functions only of r. The ranges of matching is impossible. On the other hand, as we saw, the coordinates are −∞ < t < ∞ for the time coordinate, fluids with w < −1/7 create spaces of finite 3-volume, 0 ≤ r < ∞ for the radial coordinate, −∞ < z < ∞ for fill the whole space, and the total fluid energy per unit the axial coordinate and 0 ≤ φ ≤ 2π for the angular length along the axis is finite. coordinate with the hypersurfaces φ = 0 and φ = 2π A number of papers discuss cylindrically symmetric being identified. The coordinates are numbered x0 = t, perfect fluid distributions, which can be charged or neu- x1 = r, x2 = z and x3 = φ. The general vacuum solution tral, in the presence of scalar or/and electromagnetic Rαβ = 0 for the metric (6.1), in the notation given by fields [64, 77, 205]. In [64, 77], it is shown how to ob- [121] and [125], is tain perfect fluid solutions with elecrtomagnetic fields of 2 n+1 c r 2 three possible directions admitted by the symmetry (see f = ar−n+1 − , k = −Af, eµ = r(n −1)/2,, section 3) as well as scalar fields using the inverse prob- n2a lem method, starting with certain combinations of metric r2 crn+1 l = − A2f, A = + b. (6.2) functions which can be chosen by hand. This freedom is f naf related to arbitrariness of the EoS and electric charge or current or scalar charge distributions. The regular The constants n, a, b and c can be either real or complex, axis conditions are also discussed there. In [205] differ- and the corresponding solutions belong to the Weyl class ent kinds of electrically and scalarly charged cylindrically or Lewis class, respectively. For the Lewis class these symmetric perfect fluid distributions are discussed both constants are given by in static and stationary spacetimes, a number of exact 1 2 a2 + ib2 solutions are found and analyzed. n = im, a = (a1 + b1) , b = , 2 a1 + ib1 It is known that dust (p = 0) without charges can- m 2 2 c = (a1 + b1), (6.3) not maintain static equilibrium in its own gravitational 2

field. However, for electrically or scalarly charged dust where m, a1, a2, b1 and b2 are real constants satisfying it is possible, and the corresponding solutions, in par- ticular, cylindrically symmetric ones, are known [72, 73, a1b2 − a2b1 = 1. (6.4) 229, 256]. According to [229, 256], such static equilib- rium configurations of arbitrary shape can only exist if The metric coefficients (6.2) with (6.3) become [125, 219] the mass density is everywhere equal (in proper units) 2 2 to the absolute value of the electric charge density: in f = r(a1 − b1) cos(m ln r) + 2ra1b1 sin(m ln r), this way, the electric repulsion exactly balances the grav- k = −r(a1a2 − b1b2) cos(m ln r) itational attraction. However, in the presence of a scalar −r(a1b2 + a2b1) sin(m ln r), field with its charge density, the situation is not so strictly 2 2 l = −r(a2 − b2) cos(m ln r) − 2ra2b2 sin(m ln r), determined, and an equilibrium is possible with different 2 relationships between the charge densities [72, 73]. Also, eµ = r−(m +1)/2. (6.5) in the framework of dilaton gravity (general relativity coupled to dust, scalar and electromagnetic fields, and A simple deduction of the Lewis metric, where one the latter two interact in a gauge-invariant way), it has does not need to consider complex parameters to obtain been shown that cylindrical black holes are possible, but the Lewis class, is given in [153]. There a physical in- cannot exist without negative energy density somewhere terpretation of the field equations is also provided, which in space [72, 73]. permits to get some understanding of the four parameters appearing in the Lewis solution (6.2). Another derivation These are some important results concerning static, of the Lewis metric is given in [226] where it is extended cylindrically symmetric perfect fluid distributions in gen- to Einstein’s spaces by including the cosmological con- eral relativity. stant. It is found that three parameters are essential, of 26 which one characterizes the local gravitational field, while Considering a cylindrical matter source for the Weyl the remaining two give information about the topological class metric consisting of a rigidly rotating anisotropic identification made to produce cylindrical symmetry. fluid, one of Einstein’s field equations can be integrated. In Newtonian physics the potential due to a cylindrical This integration produces matter source, being static or stationary, has the same dependence, that is, it depends only on one parameter, fk,r − kf,r = ξr, (6.9) the mass per unit length. For the static vacuum cylin- drical field in GR the solution is the LC metric, the one where the commas stand for differentiation with respect we studied in section 2, revealing two essential parame- to r and ξ is an integration constant. Calculating the rotation of the source as given in [121, 293] produces its ters, while for the stationary cylindrical rotating source it µ/2 has, in its usual form, four parameters reducible to three rotation magnitude given by ξ/(2fe ). Now using the essential parameters [226]. matching conditions on the surface of the source cylinder In the next subsection we study the meaning of the as given by Darmois [119] we find, parameters appearing in (6.2). ξ c ≡ − . (6.10) 2 B. The parameters of the Weyl class Note that this constant c is different from the speed of light used in other places of this article. In fact, now it The transformation [51] measures the rotation of the cylindrical source, as can be √ seen from (6.10). dτ = a(dt + bdφ), (6.6) In the Newtonian limit, the velocity term is negligible, 1 dφ¯ = − [cdt − (n − bc)dφ], (6.7) then from (6.10) c ≈ 0, and recalling (6.2), we find that n f = e2U . (6.11) casts the metric (6.1) with (6.2) into

n+1 Then, the Newtonian potential is 2 r ds2 = r1−ndτ 2 − r(n −1)/2(dr2 + dz2) − dφ¯2. (6.8) a 1 U = 2σ ln r + ln a, (6.12) This is locally the LC metric. Nevertheless, since φ = 0 2 and φ = 2π are identified, τ defined in (6.6) attains a where σ is given by periodic nature unless b = 0 [51]. On the other hand, the new coordinate φ¯ ranges from −∞ to ∞. A more detailed 1 σ = (1 − n). (6.13) account of this subject can be found in [292]. In order 4 to globally transform the Weyl class of the Lewis metric into the static LC metric we have to make b = 0 and In Newtonian theory, (6.12) is the gravitational poten- c = 0. Note that in this case, from the transformations tial of an infinite uniform line mass with mass per unit (6.6) and (6.7), τ and φ¯ become respectively true time length σ. The constant (ln a)/2 represents the constant and angular coordinates. arbitrary potential that exists in the Newtonian solution. Hence we can say that b and c are responsible for the The metric (6.2) has infinite curvature, according to its non-staticity of this family of solutions of the Lewis met- Cartan scalars, only at r = 0 for all n except n = ±1, ric. i.e., σ = 0 and 1/2. Thus the Weyl class metric has a As mentioned previously, the Cartan scalars provide singularity along the axis r = 0, then we can say that this the local characteristics of a metric. They are obtained spacetime is generated by an infinite uniform line source through the components of the Riemann tensor and its for 0 < σ < 1/2. covariant derivatives calculated in a constant frame. Two Considering the static limit for the Weyl class metric, metrics are equivalent if and only if there exist coordinate when n = 1 (σ = 0) and b = c = 0, we have from (6.2) and Lorentz transformations which transform the Cartan r2 scalars of one of the metrics into the Cartan scalars of the ds2 = dτ 2 − dr2 − dz2 − dφ2, (6.14) a other. By performing these calculations for both metrics, the LC metric and the Weyl class of the Lewis metric, which is the limit of the LC metric when σ = 0. In we obtain that both are equivalent locally and indistin- Section 2 it has been pointed that it generates strings guishable, which confirms the coordinate analysis made when a > 1 with mass per unit length µ = δ/4 and at the beginning of this subsection. Furthermore, we how a is directly linked to the gravitational analog of the showed that only the parameter n curves spacetime for Aharonov-Bohm effect [133]. both, static and stationary Weyl class, metrics. However, Considering c = 0 and n = 1 (σ = 0) in (6.2) we have we shall see that the two metrics possess very different √ r2  topological behaviour. ds2 = dτ 2 + 2b adτdφ − dr2 − dz2 − − b2a dφ2, Details of the calculations of the Cartan scalars are a given in [121]. (6.15) 27 producing a locally flat spacetime. In this case (6.15) rep- case, n = 1 or σ = 0, such an effect vanishes, exposing resents the exterior spacetime of a spinning string along its non-Newtonian nature. the axis of symmetry [196] with the same mass per unit For quasi-spherical objects it has been shown [183] that length µ = δ/4 but with angular momentum J given by positive radial acceleration can be produced along its axis √ of symmetry. b a It is also worth noticing that non-Newtonian forces, J = − , (6.16) 4 parallel to the z axis, also appear in the field of axially symmetric rotating bodies [48]. However, the force paral- for a > 1. lel to z in [48], unlike the current case, is directly related It has been shown [196] that a quantum scalar parti- to the spin of the source. For the Kerr black hole, it cle moving around a spinning cosmic string as given by is shown that particles, produced by a Penrose process, (6.15), exhibits a phase factor proportional to J, in its an- can be ejected from the ergosphere surface, covering the gular momentum. It is a reminiscence of the Aharonov- black hole through repulsive gravitational fields. In this Bohm effect. It is also worth mentioning that even if case too, unlike in the Lewis spacetime, the gravitational b = 0, an Aharonov-Bohm like effect (of a different kind) repulsion is created by the spin of the Kerr black hole. appears (as commented in the static case), since the an- These ejected particles are highly collimated and might gular momentum spectrum differs from the usual one, if be a mechanism for the observed extragalactic jets [154– only a > 1. 157, 255].

C. Geodesics for the Weyl class D. The parameters of the Lewis class

The Weyl class of solutions is given by (6.1) with (6.2) Using the transformation and all parameters being real, as shown above. Then, for circular geodesics [186]r ˙ =z ˙ = 0, we find that k dφ = dφ¯ + ωdt, ω = − , (1.20) 2 ˙ ˙2 f,rt˙ − 2k,rt˙φ − l,rφ = 0, where the dot stands for dif- l ferentiation with respect to an affine parameter s. The geodesic angular velocity is defined by ω ≡ φ/˙ t˙, and the metric (6.5) can be diagonalized. In order to have an the velocity of the test particle has only two nonzero integral coordinate transformation ω must be constant; t 3/2 √ φ therefore, from (6.5), m = 0. This implies, from (6.3), √components, W = kω/(f − fkω) and W = fω/(f − kω). For ω and W we obtain that n = 0 and c = 0. Thus the line element becomes 2 2 n r dr + dz c ± nω cr ± W
ds2 = − dt2 − √ − r(a2 − b2)dφ¯2. (1.21) 0 na 0 a2 − b2 r 2 2 ω = ,W = crn
, (6.17) 2 2 n − b(c ± nω0) 1 ± na W0 This is a particular case of the static LC metric with the where ω0 is the LC angular velocity and W0 is the LC energy density per unit length σ given by (6.13) being tangential velocity, equal to 1/4. Nevertheless, the transformation (1.20) is ¯ 1 − n a2 1 − n not global, since the new coordinate φ ranges from −∞ 2 to ∞, instead of ranging from 0 to 2π [121, 292]. ω0 = 2n ,W0 = . (6.18) 1 + n r 1 + n Considering, as in the case of the Weyl class, a rigidly rotating anisotropic fluid, one of the Einstein field equa- We note that ω and W vanish for ω0 = ∓c/n and W0 = ∓crn/na, respectively, which are equivalent to say that tions can be integrated producing (6.9). With (6.5) the the free particle in the present case is simply static. This matching conditions given by Darmois [119] yield could come about if the “centrifugal repulsion” balances ξ = −m(a2 + b2), (1.22) the gravitational attraction. 1 1 The geodetic motion of a particle along the axis of and so symmetry z for metric (6.1) produces ξ 1 − n2 r˙z˙ c = − . (1.23) z¨ = . (6.19) 2 2 r Hence, in order to have the rotation equal to zero, i.e. 2 2 It is interesting to note that for this geodesic the param- ξ = 0, we need m = 0 since a1 + b1 6= 0. eters b and c, due to the stationarity of spacetime, do Observe the difference, at this point, between the Weyl not appear in (6.19) and in fact it is indistinguishable class and the Lewis class. In the latter the vanishing of from its static limit, the LC spacetime. There is a force the rotation yields a locally LC spacetime, whereas in the that tends to damp the motion along the axis,z ¨ < 0, former the vanishing of the rotation does not necessarily whenever the particle approaches the axis,r ˙ < 0, and imply that the metric can be reduced either globally or reverses this tendency in the opposite case. In the flat locally to a static spacetime. 28

For the Cartan scalars that produce the local charac- with, teristics of a metric, we have the following results for the 1 (4w2R2 − 1)1/2 Lewis class. In the Lewis class as in the Weyl class, only N = (4w2R2 − 1)1/2, α = , the constant n appears in the Cartan scalars. Never- 2 2w3R4 2 2 2w R − 1 2 2 theless, now n must be substituted by its complex value β = , λ = e−w R . (1.28) (6.3) im. However, contrary to the Weyl class, the Car- 2w3R4 tan scalars for the Lewis class are distinguishable from For the case wR = 1/2 one obtains the relations from the LC metric, except for m = 0. Furthermore, there is the limits either from wR < 1/2 or wR > 1/2 which are no value of m for which the Cartan scalars are all zero, equal. The solution for wR < 1/2 belongs to the Weyl implying at once that the Lewis class does not include class and its real parameters n, a, b and c assume the the Minkowski spacetime as a special case. This fact im- values plies too that there must be a lower limit to the energy (α − β)2 per unit length of its source. The Cartan scalars impose n = (1 − 4w2R2)1/2, a = , a upper bound on the parameter m, given by 2α 1 α2 − β2 √ b = ± , c = N. (1.29) m ≤ 3, (1.24) α − β α since for larger values of m than this, the singularity is The solution wR > 1/2 belongs to the Lewis class and its at r = ∞ and not at r = 0. real parameters m, a1, a2, b1 and b2 assume the values

2 2 1/2 β 1 m = (4w R − 1) , a1 = , a2 = − , b1 b1 E. Sources producing Lewis spacetime 1/2 b1 = (−α) , b2 = 0. (1.30) In a fine paper, well ahead of his time, as observed by For the Weyl class we have the Newtonian mass per Bill Bonnor [46], in 1937 van Stockum [308] completely unit length given by σ = (1 − n)/4, which implies, for solved the problem of a rigidly rotating infinitely long case wR < 1/2, from (1.29), cylinder of dust, including the application of adequate 1 h i boundary conditions. The solution is a remarkable one. σ = 1 − (1 − 4w2R2)1/2 . (1.31) The metric for the interior is simple and unique, depend- 4 ing on one parameter, w in our notation. But for the Hence (1.31) establishes a lower limit for σ in the Lewis vacuum exterior, r > R, where R is the coordinate ra- class and being σ = 1/4. This value σ = 1/4 is the fron- dius of the cylinder, there are three cases depending on tier between the Weyl class metric and the Lewis class the mass per unit length of the interior. For the metric metric, at least for the rotating dust solution obtained (6.1) we have the following results. by van Stockum [308]. Case wR < 1/2: For the Lewis class metric the Cartan scalars, as it was remarked, do not admit the Minkowski spacetime. This  α2 + β2  f = −r 2β cosh(2N ln r) + sinh(2N ln r) , is in accordance with the existence of a lower limit for α σ in the van Stockum solution wR > 1/2, since with a  β  lower limit the source cannot be made a vacuum, and k = −r cosh(2N ln r) + sinh(2N ln r) , α therefore the exterior solution cannot be Minkowski. The Cartan scalars also impose a upper bound on the  r (4N 2−1)/2 r eµ = λ , l = sinh(2N ln r), (1.25) parameter m, given by R α √ m ≤ 3, (1.32) with, since for values of m larger than this, the singularity is at 1 (1 − 4w2R2)1/2 N = (1 − 4w2R2)1/2, α = , r = ∞ and not in r = 0. When we substitute this value 2 2w3R4 in (1.30), considering the equality, we have wR = 1. 2 2 1 − 2w R 2 2 β = − , λ = e−w R . (1.26) Van Stockum solution is studied at length in [46], 2w3R4 its properties concerning gravitoelectric and gravitomag- netic fields in [48], its confinement properties in [253], its Case wR > 1/2: extension to nonrigid rotation in [52]. A range of sta-  α2 − β2  tionary cylindrically symmetric perfect fluids sources are f = r 2β sin(2N ln r) + cos(2N ln r) , α presented in [293], a rigidly rotating dust coupled with a thin disk is found in [126], and an anisotropic cylindrical  β  k = r sin(2N ln r) − cos(2N ln r) , stationary source can be found in [127]. α Another source that consists of an infinitely rotating  r −(4N 2+1)/2 r thin shell was constructed in [259], in which the space- eµ = λ , l = cos(2N ln r), (1.27) R α time inside the shell is a constantly rotating Minkowski 29 spacetime so the symmetry axis is well defined, and the VII. LEWIS VACUUM SPACETIME shell satisfies all the energy conditions, the weak, strong COUPLED WITH A COSMOLOGICAL and dominant [177]. Later, this problem was studied in CONSTANT details in [123], and shown explicitly that the parame- ters a appearing in (6.2) and σ defined in (6.13) must A. Stationary cylindrical vacuum spacetime with be restricted to the ranges 0 ≤ σ ≤ 1/4 for a > 0 or cosmological constant 1/4 ≤ σ ≤ 1/2 for a < 0, in order to have no CTC’s out- side of the shell and meanwhile all the energy conditions Santos [275] in 1993 gave a set of solutions for sta- hold. tionary cylindrically symmetric spacetimes of the Lewis [219] form, published in 1932, with a cosmological con- stant. Krasinski [207] had previously given these solu- tions, but in a different form by using coordinates in F. Rotation and translation of cylinders which the static limit is hard to obtain. The solutions obtained by Santos are reexamined by MacCallum and In GR we show that the vacuum field produced by a Santos [226] in 1998 where they make some remarks on rotating mass cylinder is mathematically closely related the problem of the definition of cylindrical symmetry and to the field produced by a translating mass cylinder along they find too the essential parameters of the solutions. its axis of symmetry. Nonetheless, its physical and geo- We consider the cylindrically symmetric spacetime metrical properties differ significantly since the relativis- metric in the form tic frame dragging for rotation and translation physically differ considerably. ds2 = fdt2 − 2kdtdφ − dr2 − eµdz2 − ldφ2, (7.1) We assume the general cylindrically symmetric metric with its source translating parallel to its axis of symmetry where f, k, µ and l are functions of r only. The ranges of given by the coordinates are −∞ < t < ∞ for the time coordinate, 0 ≤ r < ∞ for the radial coordinate, −∞ < z < ∞ for ds2 = Adt2 − 2Kdtdz − Bdρ2 − Cdz2 − Bρ2dφ2, (1.33) the axial coordinate and 0 ≤ φ ≤ 2π for the angular co- ordinate with the hypersurfaces φ = 0 and φ = 2π being with the usual properties of its coordinates and A, K, identified. The general solution for Einstein’s equations B and C functions only of ρ. Making ρ = er and after for (7.1), in the notation of [226], is rescaling, the metric (1.33) can be written as 2/3n 2 f = G0 T1 exp[2(p0 − 1/3)U] ds2 = fdt2 − 2kdtdz − eµ(dr2 + dφ2) − ldz2. (1.34) 2 o −X1 exp[2(p3 − 1/3)U] ,

The general vacuum solution for (1.34), with the z and 2/3n φ interchanged is the stationary Lewis metric. So, the k = G0 X1X2 exp[2(p3 − 1/3)U] vacuum solution corresponding to (1.33) is simply the o Lewis solution with the coordinates z and φ interchanged. −T1T2 exp[2(p0 − 1/3)U] , Hence the metric coefficients in (1.34) are the same as in n o eµ = G2/3 Z2 exp[2(p − 1/3)U] , (6.2). 0 2 In spite of the mathematical similarity between the 2/3n 2 l = G X exp[2(p3 − 1/3)U] vacuum solutions for the fields produced by rotating and 0 2 o translating cylinder filled with perfect fluid they differ 2 −T2 exp[2(p0 − 1/3)U] , (7.2) substantially. Unlike the rotating case, the translating cylinder cannot be filled with pressure free dust, as there where the distinct metric solutions are given by is nothing here equivalent to a centrifugal force that would prevent the matter collapsing to the axis. The p pressure must therefore be nonzero. Furthermore, unlike cosh( 3|Λ|r) G0 = , Λ < 0, (7.3) the rotating case where matter can be rigidly rotating, p3|Λ| which means with shear free rotation, the translating exp(p3|Λ|r) matter case if translating rigidly can always be trans- G0 = p , Λ < 0, (7.4) formed to a frame where the system is static [170]. 3|Λ| The field of a cylinder of matter that is in translational sinh(p3|Λ|r) G0 = , Λ < 0, (7.5) motion along its length has not been studied in detail, p3|Λ| and any differences with the static case are unknown. It G = r, Λ = 0, (7.6) is therefore of interest to determine whether or not frames √0 are dragged by motion along the cylinder in a way similar sin( 3Λr) G0 = √ , Λ > 0, (7.7) to that in which they are dragged around it. 3Λ 30 with, respectively, The perfect fluid solution obtained by G¨odel[161] can p be expressed in comoving coordinates in the cylindrical U = arctan | sinh( 3|Λ|r)|, Λ < 0, (7.8) form U = exp(−p3|Λ|r), Λ < 0, (7.9) √ 2 2 p ds2 = dt2 + sinh2(ωr)dtdφ − dr2 − dz2 U = ln tanh[( 3|Λ|r)/2], Λ < 0, (7.10) ω U = ln r, Λ = 0, (7.11) 1 √ − [sinh2(ωr) − sinh4(ωr)]dφ2, (7.15) U = ln tan[( 3Λ)/2], Λ > 0. (7.12) ω2

The constants pj with j = 0, 2, 3 satisfy where ω is an arbitrary constant and the coordinates r, φ and z have the same range as in (7.1). The coordinate r X X 2η + 1 may be interpreted as the proper distance from a regular p = 1, p2 = , (7.13) j j 3 axis at r = 0, and the rotational symmetry of the metric j j (7.15) about the axis can be clearly seen. The fluid four with η = −1, 0, 1, 1, 1, respectively, for (7.3-7.7). Note velocity is uα = (ut, 0, 0, 0), and the constant mass den- sity ρ (the pressure p = 0) and cosmological constant Λ that for η = 0 the only real solution is pj = 1/3 for all j; other cases with η = 0 and all cases with η = −1, will are given by have complex values of p . The resulting general metric j 2 2 for a stationary cylindrically symmetris Einstein space 8πρ = 4ω , Λ = −2ω . (7.16) can be assumed with, say, p2 to be real, which leaves with the possibility of complex p0 and p3. The periodic char- This implies that the fluid density is positive and the acter of φ can be taken to be specified by the parameters cosmological constant negative. The metric has no cur- T2 and X2. Hence the essential parameters describing vature singularities and its form explicitly reduces to a the metric (7.2) are T2, X2, p2 and the cosmological con- vacuum Minkowski spacetime in the slow rotation limit, stant Λ. The inessential free parameters T1, X1, R and Z ω → 0. This solution may be interpreted as describing a may be required in order to match the metric to an inte- homogeneous spacetime with a fluid source that is rigidly rior by using the appropriate matching conditions. The rotating. In the Newtonian analogue of such a situation, case k = 0 in (7.2) is reducible to the static cylindrically relative to any point, there will exist a cylinder of finite symmetric LT solutions studied in Section 4. For further radius around the axis of rotation on which the particles analysis of the limits and properties of the metric (7.2) of the fluid would move at the speed of light. For all see [226]. points outside this cylinder, the particles of the classical Before proceeding further, we note that in [275] the fluid will be moving faster than light. Such a situation above solutions were presented in a slightly different is not possible in a relativistic theory, but this feature is form. In particularly, the function G was given by, reflected in the appearance of CTCs. In fact, it can be seen from (7.15) that curves on which      C cosh p3|Λ| r + C sinh p3|Λ| r , Λ < 0, t, r and z are constant are either spacelike or timelike ac-  1 2  cording to whether sinh2(ωr) is, respectively, less than or G = C1 + C2r, Λ = 0, √ √ √ greater than 1. They are null when ωr = 2 2. More-      1 C1 cos 3Λ r + C2 sin 3Λ r , Λ > 0, over, since φ = 2π is identified with φ = 0, these are (7.14) CTCs whenever r is larger than this value, but they are not geodesics.√ Hence we consider that for the interior where C1 and C2 are two arbitrary constants. However, cylinder ωr < 2 2. 0 using the transformation r = r − r0, where r0 is a con- Now we consider the matching of an interior cylinder stant, it can be seen that such given G can be always filled with a G¨odel fluid with its cylindrical surface de- brought to the forms given by (7.3)-(7.7) by properly scribed by its proper radius 0 < rS < r1 and for its exte- choosing r0 [226]. rior spacetime (7.1) with rS < r < ∞. Both proper radii in (7.1) and (7.15) are taken to be equal. Then, accord- ing to the Lichnerowicz matching conditions, the metric B. A source producing stationary cylindrical functions and their derivatives must be continuous across vacuum spacetime with cosmological constant the junction rS. By performing these long calculations one obtains the different parameters of (7.1) expressed Bonnor et al. [51] and Griffiths and Santos [169] have in terms of the radius rS and the cosmological constant identified a family of compound spacetimes for which the Λ. It can be shown that asymptotically the spacetime interior of an infinite cylinder of finite radius is described becomes the anti-de Sitter one, which is what one might by the G¨odelsolution [161]. They showed that this can be expect. Furthermore, this stationary system is expected matched to a vacuum exterior obtained by Santos [275]. to be qualitatively similar to the van Stockum solution The resulting spacetime is cylindrically symmetric and studied in Section 6. For the details of the calculations stationary. concerning this model see [51] and [169]. 31

C. Geodesic motion, confinement and where extragalactic jet formation 2 α = e−r ,  A  β = ± r2 + C, Ubiquitous extragalactic jets appear in active galaxies B and various models have been suggested for their origin 2  2 A Λ 2 D A [143]. A classic model is that of Blandford and Rees γ2 = r2 − e−2r + − r2 ± C , [39, 40], in which a hot plasma is assumed to be steadily 2B2 B2 B2 B 2 generated at the center of the galaxy. The central object 4r2e−2r is surrounded by a gravitationally confined rotating gas δ2 = , (7.18) D + 1 A2r2 − Λe−2r2 cloud. The thermal pressure of this gas cloud constrains 2 the outflowing relativistic plasma within two oppositely where A, B, C and D are four integration constants, and directed channels. An equilibrium flow is possible only if Λ is the cosmological constant. The energy density and the channel shape adjusts so that a nozzle forms where four-velocity of the dust fluid are given by the external pressure has dropped to one-half its central 1  2  value. The shape of the channel is assumed to adjust 2 2r µ µ ρ = A e + 2Λ , u = δt , (7.19) itself to a de Laval nozzle. This configuration, with a κ light fluid supporting a heavier fluid in a gravitational where κ is the Einstein coupling constant. In [258], it was field, is well known to be unstable. shown that three of these five parameters are completely In the studies of the geodesic motion of the van fixed by the elementary flatness condition (1.4) and the Stockum solution [308], Opher, Santos and Wang found existence of the symmetry axis (1.5), which lead to that the test particles are always confined inside a finite B = a2 + Λ,C = 0,D = Λ, (7.20) cylindrical shell, while their motion along the symmetry axis is free [253]. This motivated them to argue that where a ≡ A/2. Setting the confinement of the test particles in the radial direc- φ tion might provide another source for the extragalactic u = 2r2, ϕ = , (7.21) jet formation. Their arguments go as follows: The grav- a2 + Λ itational field produced by jets usually is negligible com- it can be shown that the above Wright solutions with the pared with the one produced by the matter at the center conditions given by (7.20) reduce precisely to the Lanc- of the galaxy. Thus, to the first-order approximation, it zos ones presented in [213], which represents the gravi- is sufficient to model those jets as made of test particles. tational field of a rigidly rotating dust cylinder coupled Also, almost all the galaxies are rotating, as a first step, with a cosmological constant, Λ. one can model the center of a galaxy by a rotating cylin- Then, the corresponding Lagrangian for test particles der. This approximation seems reasonable as long as the reads, gravitational field in the middle of the rotating galaxy is 1   concerned, though they admitted that it is indeed highly L ≡ t˙2 + 2βt˙φ˙ − δ2r˙2 − γ2φ˙2 − α2z˙2 , (7.22) simplified. Assuming that the above simplified model 2 can capture some essence of physics, one can see that the where t˙ ≡ dt/ds, etc., with s being the affine parameter. confinement can be related to the jets. Hence, the geodesic motions along z- and r-axes are given Certainly, to have the above idea really work, more by [258], realistic models should be built. Recently, such studies x2 z˙ = Pze , (7.23) were extended to the Kerr solution [154–157, 255], and 2 2 x2 found that particles can be ejected from the ergosphere x˙ = a e [V0 − V (x)] , (7.24) √ surface through the repulsive gravitational fields, which where x ≡ 2 r, and are created by the spin of the Kerr black hole. Again, 2
2 these ejected particles are highly collimated and might 2 a + Λ V (x) ≡ E2x2 + P 2ex + P 2 provide a mechanism for the observed extragalactic jets. z a2x2 φ  2  In the following, using the five-parameter class of Λ 1 − e−x  2  Wright solutions [331] as an example, we shall show the 2 2 x − E −  − Pz e , (7.25) generic properties of the confinements of test particles in a2x2 a2 + Λ rotating spacetimes. For this purpose, let us first note V ≡ E2 −  − 2EP , (7.26) that the G¨odel solution given by (7.15) belongs to one 0 φ a of the two classes of solutions found by Wright in [331]. where  = 0, 1, −1, respectively for null, timelike and The other class of Wright solutions contains five free pa- spacelike geodesics, and E, pφ and Pz are the constants rameters, and can be cast in the form [258], of motion, defined by ∂L ∂L ∂L 2 2 2 2 2 2 2 2 E ≡ ,Pz ≡ ,Pφ ≡ . (7.27) ds = dt + 2βdtdφ − δ dr − γ dφ − α dz , (7.17) ∂t˙ ∂z˙ ∂φ˙ 32

It is interesting to note that V (x) takes the form, (WEC) and its part, the Null Energy Condition (NEC),   at least near a wormhole throat (the narrowest place). 2 2 2 Λ x2 This result is known for a long time [224, 242, 322], at V (x) ' E x + Pz 1 + e , (7.28) a2x2 least for static wormholes with compact throats of finite as x  1. Thus, it is monotonically increasing when x area and for asymptotically flat spacetimes.[147]. becomes sufficiently large. Sincex ˙ 2 ≥ 0, from (7.24) and In theories extending GR, it is possible to find worm- (7.28) we can see that there must exist a maximal value, hole solutions without exotic matter, replacing it with geometric or field quantities that are absent in GR, e.g., say, xmax, so that (7.24) has solution for x ≤ xmax. That is, the geodesics are always confined inside a finite radius in Einstein-Cartan theory [67, 191], theories involving √ the Gauss-Bonnet invariant [74], multidimensional mod- r ≤ rmax ≡ xmax/ 2. However, the motion along the symmetry axis is free, as it can be seen from (7.23). els [11, 75, 199, 227, 281], etc. However, if we adhere to The above properties were noticed when studying the GR, which is justified by its excellent experimental sta- G¨odelspace-time, independently by Kundt [211] and tus, then a question of interest is whether or not it is Chandrasekhar and Wright [97], and later by Opher, San- possible to obtain phantom-free wormholes by abandon- tos and Wang [253] in the study of the van Stockum ing some conditions of the firmly established theorems. solution [308]. Combining these studies with the ones Meanwhile, asymptotic flatness is a necessary require- carried out above, as well as the ones carried out in the ment if a wormhole is expected to be observable from Kerr spacetime [154–157, 255], one can see that in ro- distant regions with small curvature. tating spacetimes confinement occurs generically in the Some phantom-free wormhole solutions in GR can be radial direction, whereas the motion in the axial direc- recalled, beginning with Kerr’s metric with a large angu- tion is free. Therefore, it is quite plausible to argue that lar momentum and Zipoy’s static vacuum solution [346] such confinements might provide a mechanism for extra- as well as their electromagnetic and scalar extensions galactic jet formation. [71, 239], but all of them possess a naked ring singularity As a matter of fact, in [157] the kinematic behavior of that bounds a disk playing the role of a throat. There- extragalactic jets was investigated. By numerical calcula- fore, despite their asymptotic flatness, they can hardly tions, it was found that a transition from non-relativistic be welcomed. Others [15, 89], being phantom-free, lack to ultrarelativistic speeds at subparsec scales exists. This asymptotic flatness. All such shortcomings may probably transition comes sooner and more abruptly than in mod- be interpreted as signs of the topological censorship. els based on magnetic paradigm, which indicates that Cylindrical symmetry seems to be a way to avoid topo- one needs a weaker magnetic field to explain observed logical censorship because, even assuming asymptotic synchrotron radiation. This ejection phenomenon was flatness in the radial (r) direction, we have no such flat- attributed to the repulsive effect of the gravitomagnetic ness in the longitudinal one, z, since at large |z| spacetime Kerr field. at given r remains equally curved as it was at small r. VIII. CYLINDRICAL WORMHOLES And, provided the system is t-independent (stationary), the same is true for the null (z, t) direction. So, there is The possibility of unusual spacetime geometries was hope to obtain phantom-free wormholes, asymptotically noticed by researchers as soon as the gravitational flat in the radial direction. field was identified with curvature in Einstein’s theory. Among numerous known cylindrical wormhole solu- Flamm [145] already in 1916 noticed that the spatial part tions in GR (see, among others, [66, 70, 78–81] and ref- of Schwarzschild’s spacetime describes a kind of tunnel erences therein) many are phantom-free but none are between two flat asymptotic regions, a configuration now asymptotically flat, for the same reason for which the called a wormhole after the suggestion of Wheeler [329] LC vacuum solution is only asymptotically flat if it is who also introduced the concept of spacetime foam on completely flat, that is, only when the entire spacetime extremely small length scales. Nowadays, wormholes, al- is Minkowiski. To solve the problem, the following con- though not yet discovered in Nature or constructed ar- struction was suggested [78]: a central wormhole region tificially, are one of actively discussed subjects in grav- with a throat, surrounded on both sides by Minkowski itational physics both due to their properties of inter- regions, matched through thin cylindrical shells Σ− and est (e.g., as possible distance shortcuts, time machines, Σ+. These shells will possess nonzero densities and pres- agents in quantum entanglements, etc.) and in connec- sures, their SETs being determined by jumps of the ex- tion with such fundamental issues as energy conditions trinsic curvature according to the well-known Darmois- and topological censorship [147, 148, 242, 310]. Possible Israel formalism. A phantom-free model will be obtained observable effects of wormholes in astronomy are studied if in the internal region and on the two surfaces all SETs in a number of papers [69, 130, 175, 204], see also ref- respect the WEC. erences therein. More general reviews can be found in In this section we will describe the conditions on the [82, 223, 224, 310]. geometry and SETs of matter for obtaining static or sta- One of the main problems in wormhole physsics in GR tionary (rotating) cylindrically symmetric wormhole con- is the necessity of “exotic”, or phantom matter for worm- figurations and present some particular configurations. hole existence, violating the Weak Energy Condition Among them are two examples of rotating phantom-free 33 wormhole models, asymptotically flat on both sides of which does not necessarily violate the WEP. Though, if the throat. pz = pφ (as, for example, in isotriopic fluids), (8.4) reads ρ − px < 0, violating the dominant energy condition if ρ > 0. Thus in the static case we can have phantom-free A. Wormhole definitions and existence wormhole solutions in the above sense but without two- conditions side asymptotic flatness. Apart from asymptotic flatness, other restrictions also We begin with the stationary cylindrically symmetric follow from combinations of (8.8), for example, the fol- metric lowing no-go theorem [81]: ds2 = e2γ(x)[dt − E(x) e−2γ(x) dϕ]2 − e2α(x)dx2 A static, cylindrically symmetric wormhole, mirror- 2µ(x) 2 2β(x) 2 symmetric with respect to its throat, cannot exist in GR − e dz − e dϕ , (8.1) 2 3 with matter whose SET satisfies the conditions T2 = T3 0 where x, z ∈ R and ϕ ∈ [0, 2π) are the radial, longitudi- and T0 ≥ 0. nal and angular coordinates. The coordinate x is speci- This concerns, in particular, isotropic fluids with any fied up to a reparametrization x → f(x), so its range de- equations of state and scalar fields Φ(x) with a vari- 2 pends on both its choice (the “gauge”) and the geometry ety of Lagrangians, such as Ls = (∂Φ) − 2V (Φ) with 2γ itself. It differs from (6.1) in notations (f 7→ e , k 7→ E, any potentials V and more general k-essences with Ls = µ 7→ 2µ, l 7→ e2β −E2 e−2γ ) and in that we here preserve f(X, Φ) where X = (∂Φ)2 and f is an arbitrary function. gauge freedom whereas in (6.1) the gauge (in the present What changes if there is rotation? To make it clear, notations) is chosen as α = µ. The quantity r = eβ has we note that rotation in spacetime with the metric (8.1) the meaning of circular radius. The off-diagonal compo- can be characterized by the angular velocity ω(x) of a nent E corresponds to rotation, unless E e−2γ = const, congruence of timelike curves defined as [78, 208, 209] in which case one can make E = 0 by redefining (resyn- 1 −2γ 0 γ−β−α chronizing) the time coordinate. If E = 0, the metric ω = 2 (E e ) e , (8.5) (8.1) is static and coincides with (A.1). The metric (8.1) is said to describe a wormhole if either where a prime stands for d/dx. Furthermore, in our (i) the circular radius r(x) = eβ(x) has a regular mini- reference frame is comoving with matter in its rota- 3 mum (the corresponding surface x = const is called an tion, then the SET component T0 = 0, hence (due to 3 r-throat) and is large or infinite far from this minimum Einstein’s equations) the Ricci tensor component R0 ∼ or (ii) the same is true for the area function a(x) = eµ+β (ω e2γ+µ)0 = 0, and [78] (its minimum defines an a-throat) [78, 81]. For certainty, −µ−2γ we choose the definition connected with the radius r(x), ω = ω0 e , ω0 = const. (8.6) looking more evident: while moving to smaller r, we ap- proach a would-be symmetry axis r = 0, but, instead, Consequently, from (8.5) we find, reach a minimum of r and observe its subsequent growth, Z paying no attention to what happens to the z scale. 2γ(x) α+β−µ−3γ E(x) = 2ω0 e e dx. (8.7) Asymptotic flatness in some limit of x implies r → ∞, finite limits of γ, µ and E (thus no rotation) and e−αr0 → 1 to provide a Euclidean shape of large circles. The nontrivial components of Einstein’s equations in the ν ν ν 1 ν It is clear that if a wormhole is asymptotically flat on form Rµ = κTeµ = κ(Tµ − 2 δµT ) are both extremes of the x range, it contains both r- and 0 −2α 00 0 0 0 2 0 a-throats. R0 = e [γ + γ (σ − α )] + 2ω = κTe0 , In the static case, assuming the most general SET of 1 −2α  00 02 0 0 2 1 matter compatible with the symmetry, R1 = e σ + σ − 2U − α σ − 2ω = κTe1 ,

ν 2 −2α 00 0 0 0 2 Tµ = diag(ρ, −px, −pz, −pφ) (8.2) R2 = e [µ + µ (σ − α )] = κTe2 , and using the expressions (A.2) for the Ricci tensor com- 3 −2α 00 0 0 0 2 3 R3 = e [β + β (σ − α )] − 2ω = κTe3 , ponents, it is easy to show [81] that at an a-throat near β+µ 0 0 −2γ 3 0
0 a minimum of e one should inevitably have R3 = G3 = E e R3 − R0 = κTe3 , (8.8) ρ < px ≤ 0, (8.3) where we use the notations which means the necessity of a negative density, violating 0 0 0 0 0 0 the WEC. Thus static, twice asymptotically flat worm- σ = β + γ + µ, U = β γ + β µ + γ µ . (8.9) holes are impossible. On the other hand, an r-throat only From (8.8) it follows [78] that the diagonal compo- implies ν ν nents of the Ricci (Rµ) and Einstein (Gµ) tensors split ρ − px − pz + pφ < 0, (8.4) into those for the static metric (8.1) with E = 0 plus 34 contributions involving ω: x > 0, we see that as u → 0 we have β → ∞ for any values of the constants. Let us look when β → ∞ at Rν = Rν + Rν , Rν = ω2 diag(2, −2, 0, −2), µ s µ ω µ ω µ large x ν ν ν ν 2 With a normal scalar field, ε = 1, we have h > 0,C 6= Gµ = sGµ + ωGµ, ωGµ = ω diag(3, −1, 1, −1), 0, x ∈ R+, and α ≈ −hx ⇒ 3β ≈ (3A + 3B − h)x, (8.10) and we obtain a wormhole if A + B > h/3. The metric ν ν where sRµ and sGµ are the static parts (see (A.2)). The can be written as ν ν  tensors sGµ and ωGµ (each separately) obey the con- 2 −2Bx 2 2 2 −2Ax 2 α ds = Y e dt − Y dx − e dz servation law ∇αGµ = 0 in terms of the static metric. Therefore, Gν / behaves as one more SET with exotic  2/3 ω µ κ 2Ax+2Bx 2 h 2 − e dφ ,Y := . (8.13) properties (the effective energy density is −3ω /κ < 0) λ sinh(hx) that favors the existence of wormholes, as confirmed by a number of examples [78, 79, 209]. Other wormhole solutions are obtained with h ≤ 0. By (8.8), if the diagonal components of Einstein’s Thus, if h < 0, then x ∈ (0, π/|h|), and all metric func- equations are solved, their only off-diagonal component tions α, β, γ, µ, being governed by sin(|h|x), tend to in- is automatically fulfilled. finity on both ends of this range, where repulsive singu- larities occur. Among these solutions, there is a symmetric subfamily: B. Examples of static wormholes it corresponds to ε = −1, C 6= 0, A = B = 0, h < 0. In accord with the above said, it has negative energy density, T 0 = −Φ02 e−2α +Λ/ . In terms of the Gaussian A number of such examples are presented in [81]. Let 0 κ coordinate l = (1/λ) log tan(hx/2) (l ∈ is a length us here mention two of them. The first is with an az- R in the radial direction), the metric in this case can be imuthal magnetic field, see (3.24), (3.25) in the case written as b > k. The solution is neither asymptotically flat nor regular, with the magnetic induction B and curvature |h| l 2/3 ds2 = −dl2 + cosh (dt2 − dz2 − dφ2). (8.14) invariants blowing up as u → −∞. λ λ Wormhole solutions have also been obtained with non- linear electromagnetic fields with Lagrangians Le = µν L(F ),F = Fµν F [83]. Their properties can be more C. Examples of rotating wormholes diverse, but there are restrictions following from the field equations. For example, a flat asymptotic region can ap- Assuming E 6= 0 and ω 6= 0, we return to (8.6)–(8.8) pear only on one side of the throat and, moreover, only at and see, in particular, that (8.6) is incompatible with the expense of a non-Maxwell behavior of L(F ) at small asymptotic flatness since finite values of γ(x) and µ(x) F [83]. imply finite ω. Accordingly, among numerous known ro- Another example concerns a massless scalar field Φ tating cylindrical wormhole solutions, none are asymp- µν with the Lagrangian Ls = εg Φ,µΦ,ν in the presence totically flat, and to obtain configurations potentially ob- of a cosmological constant Λ = −3λ2 < 0 (λ > 0), where servable by astronomers, one has to use matching with ε = 1 corresponds to a canonical scalar and ε = −1 to external regions on some surfaces Σ± as described above. a phantom one. Then the field equations for the metric Furthermore, in a phantom-free model, all involved (8.1) with E = 0, in the harmonic gauge α = β + γ + µ, SETs should respect the WEC, and it is a separate chal- have the family of solutions [81] lenge to find such solutions. For example, it has been 1 shown [66, 78] that a construction with two external Φ = Cx, γ = α − Ax, 3 Minkowski regions around a wormhole region does not 1 1 lead to phantom-free models: the SET on at least one β = α + Ax + Bx, µ = α − Bx, of the shells Σ violates the WEC, if the SET of mat- 3 3 ±  ter filling the internal wormhole region has the property  T 0 = T 3. This no-go theorem [66] was proved by ana-  (λ/h) sinh(hx), h > 0, 0 3  lyzing the matching conditions using the solution of the   equation R0 = R3. Therefore neither of wormhole so- e−α = (8.11) 0 3 λx, h = 0, lutions known for scalar fields (including those with ex-   ponential and other potentials [79]) are suitable for our   (λ/h) sin(hx), h < 0, purpose. where A, B, C, h are integration constants related by 2 2 2 2 1. Weak energy condition on a junction surface h sign h = 3(A + AB + B ) + 3εκC . (8.12) Among them, wormhole solutions are those where In the following, we will present two examples of mod- β(x) → ∞ on both ends of the range of x. Assuming els that solve the problem designated above. But before 35 that we formulate a convenient criterion of WEC fulfill- 2. Minkowski regions ment on a junction cylinder Σ (x = xs) separating two regions (x ≤ x ) and (x ≥ x ) with two different D− s D+ s If we place Minkowski regions M± around “internal” metrics of the form (8.1). The metric on Σ should be the rotating wormhole regions, we must take the Minkowski same when looking from D− or D+, therefore metric in a rotating reference frame. To do so, in the 2 2 2 2 2 2 inertial-frame metric dsM = dt − dX − dz − X dϕ we [β] = 0, [µ] = 0, [γ] = 0, [E] = 0, (8.15) substitute ϕ → ϕ+Ωt (Ω = const is the angular velocity) to obtain where, for any f(x), [f] = f(xs + 0) − f(xs − 0). The 2 2 2 2 2 2 coordinates t, z, φ are identified in the whole space, while dsM = dt − dX − dz − X (dϕ + Ωdt) . (8.22) the radial coordinates in D+ and D− may be different since all quantities involved in the matching conditions The relevant quantities in the notations of (8.1) are are insensitive to this choice. X2 The material content of Σ is found using the Darmois- eα = eµ = 1, e2γ = 1 − Ω2X2, e2β = , Israel formalism [28, 194] in terms of discontinuities of 1 − Ω2X2 the extrinsic curvature of Σ. Ω E = ΩX2, ω = . 1 − Ω2X2 b −1 b b b b c S = (8πG) [Ke ], Ke := K − δ K , (8.16) X 2 a a a a a c 0 Ω 0 1 Ω X γ = − 2 2 , β = + 2 2 . (8.23) 1 −α 0 1 − Ω X X 1 − Ω X where a, b, c = 0, 2, 3, and Kab = 2 e gab. The question b is whether Sa satisfies the WEC requirements This metric is stationary and suitable for matching if |X| < 1/|Ω|, so that the linear rotational velocity is a b S00/g00 = σs ≥ 0,Sabξ ξ ≥ 0, (8.17) smaller than c.

a where σs is the surface energy density and ξ an arbi- trary null vector on Σ. The second inequality makes the 3. A model with stiff matter, p = ρ content of the NEC as a part of the WEC, and taken together, they provide σs ≥ 0 in any reference frame on The perfect fluid with p = ρ is a maximally stiff matter Σ. compatible with causality, where the velocities of sound The condition σs ≥ 0 is easily found and reads and light coincide. This circumstance makes easier find- ing exact solutions, including wave and inhomogeneous [ e−α(β0 + µ0)] ≤ 0. (8.18) ones [65]; stiff matter also finds various applications in cosmology, see, e.g., [56, 251] and references therein. The radial coordinates x can be chosen in different ways A rotating wormhole solution with such a fluid reads in D+ and D−, and the prime denotes a coordinate deriva- [70] tive defined for each region separately. 2 a   2 The NEC should hold for any null vector ξ on Σ. 2 y dy 2 ds = P dt − dϕ − 2 2 2 − dz However, instead of directly verifying this, we can note ω0 2ω0(k + y ) that the WEC holds if in a comoving reference frame the 2 2 2 dϕ 2 density σ and the pressures p in all directions satisfy − (k + y ) , ρ0 = ω = const. (8.24) s i 2ω2 κ 0 the inequalities 0 The solution contains an arbitrary time scale factor P > σs ≥ 0, σs + pi ≥ 0. (8.19) 0 and integration constants k and ω0. It is regular in the 2 2 whole range y ∈ R, but at y > k it has g33 > 0, so that A difficulty is that Σ is in general not described in a ϕ-circles are CTCs. comoving frame, but the needed quantities can be found Consider matching with the Minkowski metric (8.22) as eigenvalues of Sab written in an orthonormal frame. In on Σ+ (y = y0 > 0, X = X0 > 0). The matching this way we obtain the following WEC criterion [70, 80]: conditions yield

a ≥ 0, a + c ≥ 0, P y k2 + y2 X2 P 2 = 1 − Ω2X2, ΩX2 = , = , ω 2ω2 P 2 p 2 2 0 0 a + c + (a − c) + 4d + 2b ≥ 0, (8.20) (8.25) where the index “zero” near X and y is omitted with- where out risk of confusion. We have six system parameters: four (P, k, ω , y) from the internal metric (8.24) and two −α 0 0 −α 0 0 0 a = −[ e (β + µ )] b = [ e (β + γ )], (Ω,X) from the external one. Choosing three of them −α 0 0 c = [ e (γ + µ )], d = −[ω] (8.21) as independent ones: X = X0 (having the dimension of 36 length), y = y0 and P , we obtain from (8.25) convenience. The coordinate x ranges from −x0 to x0. √ Since r → ∞ as x → ±x0, it is a wormhole solution, it 1 − P 2 P y 1 + P 2 √ 2 2 is symmetric with respect to the throat x = 0. However, Ω = , ω0 = , k = y 2 . X 1 − P 2X 1 − P the extremes x = ±x0 are curvature singularities, with −4 (8.26) the Kretschmann invariant ∼ |x0 − x| . As a result, the quantities a, b, c, d from (8.21) are ex- To construct an asymptotically flat configuration, we pressed as join it at some x = ±xs < x0 to Minkowski regions at X = ±X and obtain P 3y − 1 1 − P y s a = 2 , b = , 2 2 2 2 2 2 P X X X = r0,Q (x0 − x ) = 1 − Ω X =: P , p x0 + x P 2 − 1 P y − 1 + P 2 2x2 1 − P 2 = 2xx + (x2 − x2) ln , (8.31) c = , d = √ . (8.27) 0 0 0 x − x P 2X P 2X 1 − P 2 0 where the index s near x is omitted. To satisfy [µ] = 0, The factor 1/X is common and irrelevant for the inequal- we change the z scale in , obtaining −g = M 2 := ities of interest. An inspection shows that under the con- M± zz e2µ(xs) taken from the internal metric. dition It is convenient to use the ratio y = xs/x0 and to 2 − P 2 introduce the notation L = L(y) = ln[(1 + y)/(1 − y)]. y ≥ (8.28) P 3 So we have −(1−y)/2 −(1+y)/2 the inequalities (8.20) are satisfied. M = M(y) = 1 − y
1 + y
, What changes on Σ− specified by X = −X0 < 0 and  1  y = −y0 < 0? The common factor 1/X becomes neg- P 2(y) = (1 − y2) 1 − yL(y) − (1 − y2)L2(y) . (8.32) ative. But the signs of all discontinuities also change: 4 while on Σ we had [f] = f − f for any f, on Σ + out in − Calculation of the quantities a, b, c, d used in the WEC we have the opposite. Therefore, the WEC requirements requirements (8.20) on Σ gives do not change and hold with (8.28) applied to |y| instead ± of y. Thus we have obtained a completely phantom-free 1 M(y) y 1  a = − + + L(y) , wormhole model. P 2(y) x2 1 − y2 2 We also notice that (8.26) implies y2 < k2, hence y2 < 0 0 b = 1, k2 in the whole internal region, and there are no CTCs. 1 M(y) y 1  c = − 2 + 1 + 2 2 − L(y) , P (y) x0 1 − y 2 4. A model with an anisotropic fluid p1 − P 2(y) M(y) ± d = − + , (8.33) P 2(y) x2(1 − y2) One more model was obtained in [80] with the SET 0 where the insignificant factor 1/r at each of them has ν 0 −2γ 0 Tµ = ρ diag(1, −1, 1, −1) ⊕ T3 = −2ρE e , (8.29) been ignored. The expressions for a, b, c and d2 depend on two pa- describing an anisotropic fluid chosen by analogy with rameters, x and y, and coincide on Σ and Σ , hence that of a longitudinal magnetic field in the static case 0 + − the conditions (8.20) for the two surfaces also coincide. (see Section 3), which cannot be directly extended to A numerical analysis shows that [80]: involve rotation. (A full description of the anisotropic fluid formalism in the metric (8.1) may be found, e.g., in • The condition 0 < P 2(y) < 1, required by construction, [127].) holds for 0 < y . 0.564; The solution, in the gauge α = β + γ + µ, is written as • The conditions a > 0 and a > c hold in a signifi- −2γ−2µ ρ = ρ0 e , ρ0 = const > 0, cant range of x0 and y, for example, for x0 = 0.5, y ∈ (0.15, 0.47) and for x0 = 0.3, y ∈ (0.05, 0.53). 2 2 2β r0 2γ 2 2 2 r ≡ e = 2 2 2 , e = Q (x0 − x ), • The condition a + c > 0 also holds in a ccertain range Q (x0 − x ) of the parameters, e.g, for x0 = 0.5, y ∈ (0.15, 0.38) and 2µ 1−x/x0 1+x/x0 for x = 0.3, y ∈ (0.05, 0.51). e = (x0 − x) (x0 + x) , 0 Thus in the parameter space (x0, y) there is a signif- r  x + x icant range in which this asymptotically flat wormhole E = 0 2x x + (x2 − x2) ln 0 , 2 0 0 model completely satisfies the WEC. One can also verify 2x0 x0 − x |ω | that g33 < 0 in the internal region between Σ− and Σ+, x := 0 ,Q2 := ρ r2, (8.30) hence the model does not contain CTCs. 0 ρ r κ 0 0 κ 0 0 This section demonstrates the possible existence of The solution contains integration constants ω0, ρ0 and cylindrical wormholes asymptotically flat on both sides r0, and dimensionless constants x0 and Q introduced for of the throat without WEC violation. It should be noted 37 that the study of such systems is at an early stage, and will then oscillate until it settles down at some final, finite many questions are yet to be answered, including, above radius, whereby a spacetime singularity is prevented from all, possible observational signatures of such objects if forming on the symmetry axis [12]. It should be noted they can really exist. that AT considered only the case where the shell has zero In conclusion, we would like to mention that cylindrical total angular momentum and is momentarily static and thin shells in many studies were used not as auxiliary radiation-free. In more realistic situations, the space- objects needed to reach asymptotic regions but as the time has neither cylindrical symmetry nor zero angular main sources of gravity in wormhole solutions, see [138, momentum, and gravitational and particle radiations are 146, 277] and references therein. always expected to occur. In particular, when radiation is taken into account, modeled by an outgoing radiation fluid [216], Pereira and Wang found [261] that similar re- IX. GRAVITATIONAL COLLAPSE AND sults to a rotating spheroid obtained in [283] also hold EMISSION OF GRAVITATIONAL WAVES for a collapsing cylindrical shell made of counterrotating dust particles: in some cases the angular momentum of the dust particles is strong enough to halt the collapse, so Gravitational collapse of a realistic body has been one that a spacetime singularity is prevented from forming, of the most important problems in Einstein’s theory, and while in other cases it is not, and a line-like spacetime due to the complexity of Einstein’s field equations, the singularity is finally formed on the symmetry axis. problem even in simple cases, such as, spacetimes with Recently, East revisited the same problem studied by spherical symmetry, is still not well understood, and new Shapiro and Teukolsky in [282, 283] but using different phenomena keep emerging [198]. Particularly, in 1991 coordinates, and found that the final state in all cases Shapiro and Teukolsky studied numerically the problem studied is a black hole plus gravitational radiation [134]. of a dust spheroid, and found that only if the spheroid is Though initially distorted, the proper circumferences of compact enough, a black hole can be formed [282, 283]. the apparent horizons that are found do not significantly Otherwise, the collapse most likely ends with a naked sin- exceed the hoop conjecture bound [302]. gularity. Soon, Barrab´es,Israel and Letelier constructed an analytical model of a collapsing convex thin shell and found that in certain cases no apparent horizons are A. The Hoop conjecture and critical collapse formed, too [21]. Similar results were also found in more general cases [22]. However, since in all the cases con- sidered by them, the external gravitational field of the The studies of gravitational collapse with cylindrical collapsing spheroid/shell is not known, one cannot ex- symmetry have been generalized to various cases, includ- clude the formation of an event horizon outside of the ing the collapse of cylindrical shells with or without finite spheroid/shell. thickness [129, 135, 172, 187, 200, 201, 217, 244, 260, 287], In addition, Wald and Iyer proved that even in the and matter fields filled the whole collapsing spacetimes Schwarzschild spacetime there exists a family of Cauchy [61, 94–96, 150, 179, 212, 228, 243, 245–247, 250, 280, surfaces, which come arbitrarily close to the black-hole 284–286, 288, 332]. In particular, the general match- singularity at the center, but are such that there do not ing conditions of two arbitrary ER spacetimes with a exist any trapped surfaces lying within the past of any of thin shell were given in [260], and without a thin shell in these Cauchy surfaces [312]. Then, they argued that, in [129]. It is remarkable that in the latter the junction con- any spherically symmetric spacetime describing gravita- ditions can be explicitly solved. In addition, in this case tional collapse to a Schwarzschild black hole, it should be for a collapsing shear free isotropic cylindrical fluid, only possible to choose a spaceiike slicing by Cauchy surfaces a Robertson-Walker dust interior is possible, as shown in which terminate when a singularity (outside the matter) [129]. is reached, such that the region of the spacetime covered All the results obtained so far are consistent with the by the slicing contains no outer trapped surfaces. As a hoop conjecture [302]: black holes with horizons form result, the evidence that no outer trapped surfaces exist when and only when a mass M gets compacted into a region whose circumference is C . 4πM in every direction does not always mean that a naked singularity must be 5 resulted. In order to determine whether a black hole or . naked singularity occurs in the examples of Shapiro and In most of these studies, the spacetimes are described Teukolsky, it will be necessary to continue the evolution by the KJEK metric (1.2). Then, the no-go theorem re- of their spacetimes considerably further into the future. garding the existence of black holes in cylindrical space- times is applicable [319]. However, the naked singulari- To answer the above question, a simplified model was considered analytically by Apostolatos and Thorne (AT), that is, a cylindrical shell that is made of counter-rotating particles, and AT showed that the centrifugal forces as- 5 It should be noted that the circumference of a given region is sociated with an arbitrarily small amount of rotation, by not gauge-invariant, and a well-defined physical quantity of it is themselves, without the aid of any pressure, can halt the still absent in the general case, although it is not a problem in collapse at some non-zero, minimum radius, and the shell cylindrical spacetimes. 38 ties formed in cylindrical collapse may not be considered as counterexamples of the cosmic censorship conjecture (CCC) [257], as one can always argue that such space- times are not realistic, while the CCC is often referred to collapse of realistic matter [198]. Nevertheless, cylindri- cal spacetimes can be used to understand other impor- tant properties of realistic collapse, such as the role that rotation can play, and the nonlinear interaction between gravitational and matter fields. Critical phenomena in gravitational collapse of cylin- drically symmetric spacetimes are an example of such kind of studies. In 1993, Choptuik studied the non- linearity of the Einstein field equations near the threshold of black hole formation and revealed very rich phenom- ena [103], which are quite similar to critical phenom- ena in Statistical Mechanics and Quantum Field The- ory [19, 162]. In particular, by numerically studying the FIG. 6. The phase space of the dynamic system of the gravitational collapse of a massless scalar field in 3 + 1- Einstein-scalar field equations. The hypersurface S is the crit- dimensional spherically symmetric spacetimes, Choptuik ical surface of co-dimension one, which separates the basin found that the mass of such formed black holes takes the of black holes from the basin of flat spacetimes. A generic form, smooth family of initial data S[p] always passes the two basins at the critical point p = p∗ on the critical hypersurface. All ∗ γ MBH = C(p)(p − p ) , (9.1) the initial data on the hypersurface will collapse to the crit- ical solution A∗, which is a fixed point on the hypersurface where C(p) is a constant and depends on the initial data, when it has continuous self-similarity, and a fixed cycle when and p parameterizes a family of initial data in such a it has discrete self-similarity. All details of initial data are way that when p > p∗ black holes are formed, and when soon washed out during the collapsing process, and the col- p < p∗ no black holes are formed. It was shown that, lapse with initial data near the critical point will be very in contrast to C(p), the exponent γ is universal to all similar to the critical collapse. This similarity can last almost to the fixed point A∗, whereby the one unstable mode sud- the families of initial data studied, and was numerically ∗ denly draws the collapse either to form black hole or to a flat determined as γ ∼ 0.37. The solution with p = p , usu- spacetime, depending on whether p > p∗ or p < p∗ [316]. The ally called the critical solution, is found also universal. bottom bold horizontal line containing the point A∗ denotes Moreover, for the massless scalar field it is periodic, too. the final states of the initial data S[p]. Universality of the critical solution and the exponent γ, as well as the power-law scaling of the black hole mass all have given rise to the name Critical Phenomena in In addition, gravitational radiation from cosmic strings Gravitational Collapse. Choptuik’s studies were soon and their instabilities were also studied in cylindrically generalized to other matter fields, see, for example, the symmetric spacetimes [17, 23, 25, 27, 87, 164, 249, 323, review articles [173, 316] (For more recent studies of crit- 324]. ical phenomena, see, for example, [85, 174, 342], and ref- erences therein.). To understand the phenomena deeply, analytical stud- ies have been put forwards [105, 151, 152, 190]. How- ever, due to the nonlinearity of Einstein’s field equations, B. Polarizations and Faraday rotation even in spherically symmetric spacetimnes, the problem is still too complicated to be studied analytically. In As is well known, due to Birhkoff’s theorem [177], the 2003, cylindrically symmetric spacetimes with homoth- spherically symmetric vacuum spacetimes are uniquely etic self-similarity were studied and a class of exact solu- described by the Schwarzschild black hole solution, and tions to the Einstein-massless scalar field equations was it is not possible to have gravitational radiation in such found [317]. Among other things, from the analysis of spacetimes. So, one of the next simplest cases is space- their linear perturbations, it was found that there exists times with cylindrical symmetry. Studying cylindrical one solution that has one and only unstable mode. By GWs, Piran, Safier and Stark (PSS) first discovered that, definition [173, 316], this is a critical solution that may due to the nonlinear interaction between ingoing and out- sit on a boundary that separates two different basins of going cylindrical GWs, a phenomenon similar to the elec- attraction in the phase space [see Fig.6]. tromagnetic Faraday rotation exists, but now one GW Recently, a self-similar cylindrical scalar field with serves as the medium for the other, so the polarizations non-minimal coupling was studied [106, 107], and a 2- of these waves get rotated [263, 264]. The spacetimes parameter family of solutions with a regular axis was considered by PSS are the ones described by the KJEK found, and their global structure was studied. metric (1.2) with W = r, which is possible when space- 39

314]. Let√ us first introduce√ two null coordinates u ≡ (t + r)/ 2, v ≡ (t − r)/ 2, so the metric (1.2) takes the form, ds2 = 2e−M dudv − e−U eV cosh Wdφ2 − 2 sinh Wdφdz +e−V cosh Wdz2 , (9.4) where 1 ψ = − (U + V − ln cosh W) , 2 1 γ = − (M + U + V − ln cosh W) , 2 cosh1/2 W ω = −eV tanh W,W = e−U . (9.5) sinh W Introducing a null tetrad (lµ, nµ, mµ, m¯ µ)[248], via the FIG. 7. (a) A null geodesic congruence meets the SO plane relations, in the circle S. Because of the force generated by the Ψ0 µ µ µ µ GW component, the image of the circle S on the SP plane l = Bδv , n = Aδu , is a deformed into an ellipse [314]. (b) When a matter wave mµ = ζ2δµ + ζ3δµ, m¯ µ = ζ2δµ + ζ3δµ, passes through the circle S located on the SO plane, its image 2 3 2 3 (U−V )/2   on the SP plane is a uniformly-contracted circle, due to the e W W ζ2 = √ cosh + i sinh , attractive forces of matter. (c) When a Ψ0 GW with only 2 2 2 the “+” polarization passes through the circle S, its image e(U+V )/2  W W  is distorted into an ellipse with its major axis along the e2- ζ3 = √ sinh + i cosh , (9.6) direction. (d) When a Ψ0 GW with only the “×” polarization 2 2 2 passes through the circle S, its image is distorted into an ellipse but now with its major axis along the direction with M = ln(AB), one can show that each of the two null 0 µ µ 45 with respect to the e2-direction. vectors l and n defines a null geodesic congruence,

ν µ µ ν µ µ l ∇ν l = −B(ln A),vl , n ∇ν n = A(ln B),un . times are vacuum, (9.7) Thus, when A = 1 (B = 1) the null geodesics defined by 2 2(γ−ψ) 2 2
2ψ 2 2 −2ψ 2 ds = e dt − dr −e (dz + ωdφ) −r e dφ . lν (nµ) are affinely parameterized. (9.2) Without loss of generality, let us first consider the null Introducing the quantities, geodesics defined by lν . Let ηµ be the deviation vector ν between two neighbor geodesics and η lν = 0. Then, it I+ ≡ 2 (ψ,t + ψ,r) ,O+ ≡ 2 (ψ,t − ψ,r) , can be shown that the geodesic deviation is given by e2ψ e2ψ I× ≡ (ω,t + ω,r) ,O× ≡ (ω,t − ω,r) , (9.3) D2ηµ r r = −Rµ lnulρηλ = −eM Φ eµν − 2Re (Ψ ) eµν Dλ2 νλρ 00 o 0 + PSS interpreted I ,I and O ,O as representing in- µν + × + × −2Im (Ψ ) e  η , (9.8) going and outgoing cylindrical GWs, respectively. The 0 × ν indices + and × denote the two independent polariza- where Re (Ψ0) (Im (Ψ0)) denotes the real (imaginary) µ ν λ ρ µ ν tions. With the boundary conditions at the symmetry part of Ψ0 (≡ −Cµνλρl m l m ), 2Φ00 ≡ −Rµν l l axis (r = 0), I+ = O+,I× = −O×, PSS found nu- [248], and merically that the propagation of these waves displays a µν µ ν µ ν µν µ ν µ ν reflection of ingoing to outgoing waves (and vice versa), eo ≡ e2 e2 + e3 e3 , e+ ≡ e2 e2 − e3 e3 , in combination of a rotation of the polarization vector µν µ ν µ ν e× ≡ e2 e3 + e3 e2 , (9.9) between the + and × modes. √ √ Note that the definition of the + and × modes given in with e2 ≡ (m +m ¯ )/ 2, e3 ≡ −i(m − m¯ )/ 2. The Weyl (9.3) is coordinate-dependent, and the polarizations are scalar Ψ0 is interpreted as representing the transverse defined with respect to the chosen coordinates. In [313, GW component propagating alone the nν direction [298, 314], a coordinate-independent definition was given for 314]. plane GWs. However, the difference between metric (9.2) For physically realistic matter, we have Φ00 ≥ 0 [177]. and plane symmetric spacetimes only lies in the existence Equation (9.8) allows the following physical interpreta- of a symmetry axis and the closed orbit of the Killing tion. Let S0 and SP be infinitesimal 2-dimensional ele- vector ξ = ∂φ. Therefore, it can be easily generalized ments spanned by the two spacelike vectors e2 and e3, to the current case. In the following, we shall give a and orthogonal to a null geodesic C at neighbor points brief outline, and for details we refer readers to [313, O and P of C, and let S be an infinitesimal circle with 40 center O, lying in S0 [see Fig.7]. Suppose that a beam of need to find a parallel-transported basis, and then define light rays meets S0 in the circle S, then let us observe the the polarization angle with respective to this parallel- image of these light rays on SP . The first term on the transported basis. In [314], it was shown that if (e2, e3) right-hand side of (9.8) shows that matter fields always (0) is rotated by an angle ϕ0 , make the circle S contracted [see Fig.7 (b)]. The second term, corresponding to the contribution of the real part (0) 2ϕ0,u − sinh WV,u = 0, (9.14) of Ψ0, makes the circle elliptic with the main major axis along e2 [see Fig.7 (c)]. The last term on the right-hand the resulted basis is parallel-transported along the path side of (9.8), corresponding to the contribution of the of the Ψ0 wave. Thus, the angle imaginary part of Ψ0, makes the circle also elliptic but (0) with the main major axis tilted at 45 to e2 and e3 [see θ0 ≡ ϕ0 − ϕ0 , (9.15) Fig.7 (d)]. Thus, the image of these light rays on SP is an ellipse. Making the rotation in the (e2, e3)-plane with determines the polarization direction of the Ψ0 wave with an angle ϕ0 defined by respect to this parallel-transported basis. Similarly, the angle Im (Ψ0) tan 2ϕ0 = − , (9.10) (0) Re (Ψ0) θ4 ≡ ϕ4 − ϕ4 , (9.16) 0 0 and then in terms of (e2, e3), (9.8) takes the form, determines the polarization direction of the Ψ4 wave with 2 µ respect to the parallel-transported basis along the path D η µ nu ρ λ = −R l l η of the Ψ4 wave, where Dλ2 νλρ M h 0µν
1/2 0µν i (0) = −e Φ00e o − 2 Ψ0Ψ0 e + ην .(9.11) 2ϕ4,v − sinh WV,v = 0. (9.17) It follows that the main major axis of the ellipse is along Note that in the derivations of the polarizations of the 0 0 e2-direction. We call e2 the direction of the polarization cylindrical GWs, we used the null geodesic deviations. of the Ψ0 wave. The angle ϕ0 is the polarization an- One can equally use timelike geodesic deviations as done gle of this wave component with respect to the (e2, e3)- in [313] (See also [298] for the general case), and the same basis. The relative accelerations of neighboring geodesics results will be obtained, as it is expected.
1/2 Moreover, for the rotating cylindrical GWs described are proportional to Ψ0Ψ0 , which does not relate to
1/2 by metric (1.3), choosing the null tetrad any observer. Thus, Ψ0Ψ0 represents the absolute strength of the relative accelerations of neighboring rays. eψ−γ Aeψ−γ lµ = √ (1, 1, 0, −ω) , nµ = √ (1, −1, 0, −ω) , Similarly, we can consider the geodesic deviation of the 2 A 2 null congruence defined by nµ, and the geodesic deviation 1  eψ  will take the form, mµ = √ 0, 0, e−ψ, i , 2 W 2 µ D η µ ν ρ λ M  µν µν  ψ  = −R n n η = −e Φ e − 2Re (Ψ ) e µ 1 −ψ e 2 νλρ 22 o 4 + m¯ = √ 0, 0, e , −i , (9.18) Dλ W µν  2 +2Im (Ψ4) e× ην , (9.12) Pereira and Wang studied the polarizations of the rotat- µ ν λ ρ µ ν where Ψ4 ≡ −Cµνλρn m¯ n m¯ , and 2Φ22 ≡ −Rµν n n ing GWs, and found that the geodesic deviation equation [248]. The Weyl scalar Ψ4 is interpreted as representing (9.8) is replaced by [259], the transverse GW component propagating alone the lν D2ηµ direction [298, 314]. Rotating the basis (e2, e3) now with = −Rµ lν lβηα = Φ eµν + Ψ eµν an angle, Dλ2 ναβ 00 o 0 + µν µν +i (Ψ1 + Φ01) e03 + i (Ψ1 + Φ01) e13 ] ην , Im (Ψ4) tan 2ϕ4 = , (9.13) (9.19) Re (Ψ4) µ ν λ ρ where Ψ1 ≡ −Cµνλρl n l m , and 2Φ01 ≡ the geodesic deviation equation (9.12) takes a similar µ ν γ−ψ (Rµν − Rgµν /4) l m [259], A = e , form of (9.11). Thus, the angle ϕ4 denotes the polar- µν µ ν µ ν µν µ ν µ ν ization angle of the ψ4 wave component with respect to e03 ≡ e0 e3 + e3 e0 , e13 ≡ e1 e3 + e3 e1 , µ µ µ µ the (e2, e3)-basis. µ l + n µ l − n It is important to note that the basis (e , e ) defined e0 ≡ √ , e1 ≡ √ , (9.20) 2 3 2 2 above is not parallel-transported (freely falling) along µν the null geodesic congruence. As a result, the angles and e2, e3 and e+ are as defined above. It is remarkable ϕ0 and ϕ4 defined above have physical meaning only lo- to note that now the Ψ0 wave has only one polarization, cally. Thus, in order to compare the polarizations of a the “+” mode, along e2 direction, and the “×” mode now given GW at different locations along the wave path, we is absent. 41

the Ψ1 term deflects the sphere into an ellipsoid [298]. Moreover, the second term in (9.19) is absent in the time- like geodesic case. The effect of this term will make a clock “flying” with the photons slow down, in addition µ to the effect of deflecting the photons in the e3 -direction. It is interesting to note that there is a fundamental dif- ference between the time delay caused by this term and the one caused by a Lorentz boost. The latter, in partic- ular, has no contribution to geodesic deviations, timelike or null [259]. From the above analysis we can also see that for a pure Petrov type N GW propagating along the null geodesic congruence, in which only the component Ψ0 is differ- ent from zero, the GW has only one polarization state, similar to the KJEK metric (9.4) with W = 0. The differ- ence between these two cases is that in the case W = 0, the polarization angle remains the same even at different points along the wave path for the KJEK metric (9.4), while in the rotating case (1.3) with ω 6= 0, in general this is no longer true. In fact, it is easy to show that ω, eµ lν = 0, eµ lν = −W eψ−χ √r (eµ + eµ) . (9.22) 2;ν 3;ν 2 2 0 1

Thus, for the spacetimes described by metric (1.3), al- µ though e2 is parallel-transported along the null geodesic µ congruence, e3 in general is not, and is rotating with re- FIG. 8. A spherical ball consisting of photons cuts SO in the spect to a parallel-transported basis. Since the polariza- circle S with the point O as its center. The image of the ball tion angle of the Ψ0 wave remains the same with respect at the point P is turned into a spheroid with the main major µ axis along a line at 450 with respect to e in the plane spanned to e3 , the polarization direction is also rotating with re- 1 spect to the parallel-transported basis. by e1 and e3 because of the interaction of Ψ1 and Φ01, while the rays are left undeflected in the e2-direction [259]. Note that the effects of polarizations on the formation of spacetime singularities were studied for plane GWs [315], and one can generalize such studies to cylindrical To see the physical interpretation of the last two terms, GWs, and similar results should be expected to hold here, let us consider a tube along the null geodesic C. Con- too. sider a sphere consisting of photons, which will cut SO in the circle S with the point O as its center, as shown in Fig.8. Then, the last term in (9.19) will make the image X. CONCLUSIONS of the sphere at the point P as a spheroid with the main 0 µ major axis along a line at 45 with respect to e1 in the µ µ With the beginning of the era of gravitational wave plane spanned by e1 and e3 , while the rays are left un- µ astronomy [4], the strong gravitational field regime will deflected in the e2 -direction. This can be seen clearly by µ µ be soon explored observationally in various aspects. Such performing a rotation in the plane spanned by e1 and e3 , µ 0µ 0µ µ 0µ 0µ theoretical studies can be carried out analytically and/or e1 = cos αe 1 − sin αe 3 , e3 = sin αe 1 + cos αe 3 , which leads to numerically. In the former, due to the complexity of Ein- stein’s field equations, symmetries of spacetimes are often µν 0µ 0ν 0µ 0ν 0µ 0ν 0µ 0ν e13 = sin(2α)(e 1 e 1 − e 3 e 3 ) + cos(2α)(e 1 e 3 + e 3 e 1 ). imposed, such as spherical, plane and cylindrical [294]. (9.21) Although they are all ideal models, and in realistic situa- µν 0µ 0ν Thus, choosing α = π/4, we have e13 = e 1 e 1 − tions any of these symmetries may not exist, they do pro- 0µ 0ν e 3 e 3 , (α = π/4). Combining it with (9.19), we can vide attackable examples, from which some fundamental 0 0 see that the last term will make a circle in the (e 1, e 3)- issues of physics can be addressed. One of the typical plane into an ellipse with its main major axis along the cases is the Schwarzschild solution, which plays the cen- 0 0 e 1-axis, which is at 45 with respect to the e1-axis [cf. tral role in the studies of black hole physics. Another Fig.8]. case is the LC solution, discovered by LC in 1919 [218], Similar physical interpretation can be given to the sec- marked as the second exact solution to the Einstein’s ond term, but note that e0 defined by (9.20) now is time- field equations discovered historically, and has been stud- like. ied extensively in the past decades, in the efforts to un- It should be noted that in the case of timelike geodesics derstand the nonlinearity of Einstein’s theory. 42

In this article, we have first defined cylindrically sym- cently, it was found that LISA will be able to detect such 2 −17
metric spacetimes [92], and shown that the most general a background with Gµ/c & O 10 [14], while cur- form of the metric can be cast in the form (1.6)[311], rently the most stringent bounds come from pulsar tim- 2 with the conditions (1.4) and (1.5) as the existence of a ing arrays [38], which yield the upper bound Gµ/c . symmetry axis. A spacetime without a symmetry axis is O 10−11
. On the other hand, depending on the string referred to as a cyclically symmetric spacetime [20], as network model, LIGO-Virgo observations could constrain 2 −14
the axis may be singular and so not part of the manifold, it to Gµ/c . O 10 [5,6]. or the topology of the manifold may be such that no axis So far, cylindrical spacetimes have been studied mainly exists. Depending on additional assumptions, the metric in the form of the KJEK metric (1.2). As shown in In- can reduce to the KJEK metric (1.2), the rotating MQ troduction, this is only a particular case of the general metric (1.3), or the ER metric (1.11). metric (1.6). So, a natural question is to which extent Then, we have provided a general review on the main the properties obtained so far hold in more general cases? properties and possible sources of the LC solution (Sec. What are the generic features of Einstein’s theory of grav- 2), when it is coupled with an electromagnetic field (Sec. ity? which are independent of the symmetry assumed 3), or with the cosmological constant (the LT solution, but obtained from the studies of cylindrical spacetimes. Sec. 4), or with a perfect fluid (Sec. 5). Which kind of guidelines can one draw for the future The Lewis vacuum solution is the generalization of the observations of gravitational waves from such studies? static LC solution to the case with rotation, and has been Therefore, despite of the fact that cylindrical spacetimes reviewed in Sec. 6. When it is coupled with a cosmolog- have been extensively studied in the past several decades, ical constant, the solution was first studied by Lanczos and various remarkable results were obtained, there are [213], then by Wright [331], Krasimski [207], Santos [275], still many important questions that are still waiting for and Pereira, Santos and Wang [258], among many others. answers. One of the remarkable features is that the confinement of particles near the symmetry axis. Such particles are highly collimated and could provide an alternative mech- ACKNOWLEDGEMENTS anism for the formation of the observed extragalactic jets [154, 155, 157, 253, 255, 258]. A.W. would like to thank B. Mashhoon for valu- In Sec. 8, the studies of cylindrically symmetric worm- able discussions and suggestions. He is supported in holes have been reviewed, and various models have been part by the National Natural Science Foundation of presented, while in Sec. 9, the gravitational collapse China (NNSCF) with the Grants Nos. 11675145 and of cylindrically symmetric sources has been summarized, 11975203. The work of K.B. was partly supported by and particular attention has been paid to the understand- the RUDN University program 5-100 and was also per- ing of the Hoop conjecture [302] and critical phenomena formed within the framework of the Center FRPP sup- at the threshold of black hole formation [103]. In this sec- ported by MEPhI Academic Excellence Project (contract tion, a gauge-invariant definition of the polarizations of No. 02.a03.21.0005, 27.08.2013). cylindrically symmetric gravitational waves has been also given [313, 314], and the interesting features of gravita- tional Faraday rotation [263, 264] of such waves have been APPENDIX A: STATIC CYLINDRICAL reviewed. Thus, now it would be extremely interesting SPACETIMES to find some observational effects from polarizations of gravitational waves and gravitational Faraday rotation. In this Appendix, we shall present some general state- With the promise of increasing duration, observational ments and relations about the static cylindrical space- sensitivity and the number of detectors, many more grav- times. itational wave events are expected to be detected. In par- The most general static metric admitting translational ticular, LISA [118] is expected to observe tens of thou- symmetry along two mutually orthogonal axes has the sands of compact galactic binaries during its nominal four form [Cf. (1.11)], year mission lifetime [113]. Then, the properties of the strong field regime of gravitational fields will be explored ds2 = e2γ(u)dt2 − e2α(u)du2 − e2µ(u)dz2 − e2β(u)dφ2, in various directions, including the (nonlinear) interac- (A.1) tion of gravitational waves and matter fields. where u is an arbitrary admissible “radial” coordinate, In addition to GWs produced by these putative with a freedom of choosing it for one’s convenience. It is sources, another important target of GW experiments important that the coordinates z and φ appear in (A.1) is the stochastic gravitational wave background (SGWB) (and hence in the field equations) on equal grounds and of cosmological origin [112]. One of such cosmological admit different interpretations. Thus, assuming z ∈ R sources is the cosmic superstring network [91], which and φ ∈ R, we obtain a symmetry called pseudoplanar, emits GWs throughout the history of the Universe, and being an extension of planar symmetry (the latter re- generates a SGWB from the superposition of many un- quires additional isometry under rotations in the (z, φ) correlated sources, as mentioned in Introduction. Re- plane). The assumption of both z ∈ S1 and φ ∈ S1 leads 43 to toroidal symmetry. We here choose to deal with cylin- Let us here, for illustration, show how the LC vacuum ν drical symmetry, assuming that z ∈ R is the longitudinal solution is obtained using (A.4). The equations Rµ = 0 coordinate, and φ ∈ S1 is the angular one. then immediately lead to β00 = γ00 = µ00 = 0, whence Let us present for further calculations the nonzero com- ponents of the Ricci tensor for the metric (A.1) without β(u) = bu + b , γ(u) = cu + c , µ(u) = hu + h , (A.5) specifying the choice of the coordinate u, 0 0 0

0 −2α 00 0 0 0 0 0 R0 = e [γ + γ (γ − α + β + µ )] , with six integration constants, among which b0, c0, h0 1 −2α  00 00 00 02 02 can be turned to zero by changing scales along the t and R1 = e γ + µ + β + γ + µ 02 0 0 0 0  z axes and by choosing the zero point of the u coordi- +β − α (γ + µ + β ) , 1 nate. Moreover, the first-order equation G1 = 0 leads to 2 −2α 00 0 0 0 0 0 R2 = e [µ + µ (γ − α + β + µ )] , a relation between the remaining constants b, c, h: 3 −2α 00 0 0 0 0 0 R3 = e [β + β (γ − α + β + µ )] , (A.2) bc + bh + ch = 0, (A.6) where the prime denotes d/du, and the coordinates are numbered according to the scheme (0, 1, 2, 3) = 1 hence there are precisely two significant constants, in ac- (t, u, z, φ). It is also useful to write the component G1 of ν ν 1 ν cord with what was said previously. The metric takes the the Einstein tensor Gµ = Rµ − 2 δµR that contains only first-order derivatives: form

1 −2α 0 0 0 0 0 0 G1 = − e (γ µ + β γ + β µ ) . (A.3) ds2 = e2cudt2 − e2(b+c+h)udu2 − e2hudz2 − e2budφ2. (A.7) The coordinate u can be chosen in different ways, de- It is quite suitable for a further study, and its form (2.8), pending on the purpose of a study, convenience and taste, preferred for historical reasons, is obtained from (A.7) and is usually fixed by some relation between the func- using (A.6) and the following relations and notations: tions α, β, γ, µ. For example, the LC solution in the pre- vious section was written under the conditions α = µ [Cf. (2.8)] and α = 0 [Cf. (2.14)]. e(b+c)u = kρ, k = (b + c)(b+c)/(b+c+h), ν However, the components Rµ take the simplest form if 2σ = c/(b + c), a2 = (b + c)−2b/(b+c+h). (A.8) we choose the harmonic radial coordinate defined by the condition, This transformation is valid for b+c 6= 0 and also requires α = β + γ + µ. (A.4) some finite rescaling along the t and z axes.

[1] Aasi J. et al. (LIGO Scientific and VIRGO Collabora- v. D97, 102002, 2018. tions). Constraints on cosmic strings from the LIGO- [7] Ade P.A.R. et al. Planck 2015 results. XIII. Cosmolog- Virgo gravitational-wave detectors. Phys. Rev. Lett. v. ical parameters. Astron. Astrophys. v. 594, A13, 2016. 112, 131101, 2014. [8] Adams J.F. (1996). Lectures on exceptional Lie [2] B.P. Abbott, et al. (LIGO and Virgo Collaborations). groups, Chicago Lectures in Mathematics (University ASTROPHYSICAL IMPLICATIONS OF THE BI- of Chicago Press). NARY BLACK HOLE MERGER GW150914. Astro- [9] Amaro-Seoane P., Santamaria L. Detection of IMBHs phys. J. Lett. v. 818, L22, 2016. with ground-based gravitational wave observatories: A [3] Abbott B.P. et al. (LIGO & Virgo Scientific Collabo- biography of a binary of black holes, from birth to death. ration). Search for transient gravitational waves in co- Astrophys. J. v. 722, 1197, 2010. incidence with short duration radio transients during [10] Andrzejewski K., Prencel S. Memory effect, conformal 2007-2013. Phys. Rev. v. D93, 122008, 2016. symmetry and gravitational plane waves. Phys. Lett. v. [4] Abbott B.P. et al. (The LIGO Scientific and Virgo Col- B782, 421-426, 2018. laborations) GWTC-1: A gravitational-wave transient [11] Aneesh S., Bose S., Kar S. Gravitational waves from catalog of compact binary mergers observed by LIGO quasinormal modes of a class of Lorentzian wormholes. and Virgo during the first and second observing runs. Phys. Rev. v. D97, 124004, 2018. arXiv:1811.12907, 2018. [12] Apostolatos T.A., Thorne K.A. Rotation halts cylin- [5] Abbott B.P. et al. (The LIGO Scientific and Virgo Col- drical, relativistic gravitational collapse. Phys. Rev. v. laborations). Search for the isotropic stochastic back- D46, 2435-2444, 1992. ground using data from Advanced LIGO’s second ob- [13] Arafah M.R., S. Fakioglu S., Halilsoy M. Exact Fara- serving run. arXiv:1903.02886, 2019. day rotation in the cylindrical Einstein-Maxwell waves. [6] Abbott B.P. et al. (The LIGO Scientific and Virgo Col- Phys. Rev. v. D42, 437-440, 1990. laborations). Constraints on cosmic strings using data [14] Auclair P. et al. Probing the gravitational wave from the first Advanced LIGO observing run. Phys. Rev. background from cosmic strings with LISA. 44

arXiv:1909.00819, 2019. and their mathematics. AMS Notices, v. 64, 07, 2017; [15] Ay´on-BeatoE., Canfora F., Zanelli J. Analytic self- Gravitational wave memory in ΛCDM cosmology. Class. gravitating Skyrmions, cosmological bounces and AdS Quantum Grav. v. 34, 215002, 2017. wormholes. Phys. Lett. B 752, 201-205, 2016. [36] Binetruy P., Bohe A., Hertog T., Steer D.A. Gravita- [16] Bach R., Weyl H. Neue L¨osungender Einsteinschen tional wave signatures from kink proliferation on cosmic gravitationsgleichungen. B. Explizite aufstellung statis- (super-) strings. Phys. Rev. v. D82, 126007, 2010. cher axialsymmetrischer felder. Miteinem Zusattz ¨uber [37] Blanchet L., and T. Damour T. Hereditary effects in das statische Zweik¨orperproblem von H. Weyl. Math. gravitational radiation. Phys. Rev. v. D46, 4304-4319, Zeits. v. 13, 134-145, 1922. 1992. [17] Banerjee A., Banerjee N., Sen A. A. Static and nonstatic [38] Blanco-Pillado J. J., Olum K. D., Siemens X. New limits global string. Phys. Rev. v. D53, 5508-5512, 1996. on cosmic strings from gravitational wave observation. [18] Barause E., Yunes N., Chamberlain K. Theory-Agnostic Phys. Lett. v. B778, 392-396, 2017. Constraints on Black-Hole Dipole Radiation with Multi- [39] Blandford R.D., and Rees M.J. A ‘twin-exhaust’ model band Gravitational-Wave Astrophysics. Phys. Rev. for double radio sources. Mon. Not. R. Astron. Soc. v. Lett. v. 116, 241104, 2016. 169, 395-415, 1974. [19] Barenblatt G.I. Similarity, self-similarity, and inter- [40] Blandford R.D., and Rees M.J. Extragalactic double ra- mediate asymptotics. Consultants Bureau, New York, dio sources-the current observational and theoretical po- 1979. sition. Contemp. Phys. v. 16, 1, 1975. [20] Barnes A. A comment on a paper by Carot et al. Class. [41] Bondi H., Pirani F.A.E. Energy conversion by gravita- Quantum Grav. v. 17, 2605-2609, 2000. tional waves. Nature, v. 332, 212, 1988. [21] Barrab´esC., Israel W., Letelier P.S. Analytic models of [42] Bondi H., Pirani F.A.E. Gravitational waves in general nonspherical collapse, cosmic censorship and the hoop relativity. 13: Caustic property of plane waves. Proc. R. conjecture. Phys. Lett. v. A160, 41-44, 1991. Soc. Lond. v. A421, 395-410, 1989. [22] Barrab´esC., Gramain A., Lesigne E., Letelier P.S. Ge- [43] Bonnor W.B. Non-singular fields in general relativity. J. ometric inequalities and the hoop conjecture. Class. Math. Mech. v. 6, 203-214, 1957. Quantum Grav. v. 9, L105-110, 1992. [44] Bonnor W. B. A rotating dust cloud in general Relativ- [23] Battye R. A., Shellard E. P. S. Radiative back reaction ity. J. Phys. A: Math. Gen. v. 10, 1673-1677, 1977. on global strings. Phys. Rev. v. D53, 1811-1826, 1996. [45] Bonnor W. B. Solution of Einstein’s equations for a line [24] Beck G. Zur theorie binarer gravitationsfelder. Z. Phys. mass of perfect fluid. J. Phys. A: Math. Gen. v. 12, 847- v. 33, 713-728, 1925. 851, 1979. [25] Bedran M.L., de Sa E.M., Wang A., Wu Y. On the back [46] Bonnor W. B. The rigidly rotating dust cylinder. J. reaction of gravitational and particle emission and ab- Phys. A: Math. Gen. v. 13, 2121-2132, 1980. sorption from straight thick cosmic strings: a toy mode. [47] Bonnor W. B. Physical interpretation of vacuum solu- Gen. Relativ. Grav. v. 31, 1769-1776, 1999. tions of Einstein’s equations. Part I. Time-independent [26] Bennett C.L., Banday A., Gorski K.M., Hinshaw G., solutions. Gen. Relativ. Gravit. v. 24, 551-574, 1992. Jackson P., Keegstra P., Kogut A., Smoot G. F., Wilkin- [48] Bonnor W. B. The magnetic Weyl tensor and the van son D. T., Wright E. L. 4-Year COBE DMR cosmic mi- Stockum solution. Class. Quantum Grav. v. 12, 1483- crowave background observations: maps and basic re- 1489, 1995. sults. Astrophys. J. v. 464, L1-L4, 1996. [49] Bonnor W.B. A second anti-de Sitter universe. Class. [27] Bennett D. P., Bouchet F. R. Constraints on the gravity- Quantum Grav. v. 25, 225005, 2008. wave background generated by cosmic strings. Phys. [50] Bonnor W. B., Griffiths J. B., MacCallum M. A. H. Rev. v. D43, 2733-2735, 1991. Physical interpretation of vacuum solutions of Einstein’s [28] Berezin V.A., Kuzmin V.A., Tkachev I.I. Dynamics of equations. Part II. Time-dependent solutions. Gen. Rel- bubbles in general relativity. Phys. Rev. v. D36, 2919- ativ. Gravit. v. 26, 687-729, 1994. 2944, 1987. [51] Bonnor W. B., Santos N. O., MacCallum M. A. H. An [29] Berger B.K., Chrusciel P.T., Moncrief V. On “asymp- exterior for the G¨odelspacetime. Class. Quantum Grav. totically flat” spacetimes with G2-invariant Cauchy sur- v. 15, 357-366, 1998. faces. Ann. Phys. v. 237, 322-354, 1995. [52] Bonnor W. B., Steadman B. R. A vacuum exterior to [30] Berti C., et al. Testing general relativity with present Maitra’s cylindrical dust solution. Gen. Rel. Grav. v. 41, and future astrophysical observations. Class. Quantum 1381-1387, 2009. Grav. v. 32, 243001, 2015. [53] Braginsky V.P., Grishchuk L.P. Kinematic resonance [31] Berti E., Yagi K., Yunes N. Extreme gravity tests with and memory effect in free-mass gravitational antennas. gravitational waves from compact binary coalescences: Soviet Phys. JETP v. 62, 427-430, 1985 [Zh. Eksp. Teor. (I) inspiral-merger. Gen. Relativ. v. 50, 46, 2018. Fiz. v. 89, 744, 1985]. [32] Berti E., Yagi K., Yunes N. Extreme gravity tests with [54] Braginsky V.B., Thorne K.S. Gravitational-wave bursts gravitational waves from compact binary coalescences: with memory and experimental prospects. Nature, v. (II) ringdown. Gen. Relativ. v. 50, 49, 2018. 327, 123-125, 1987. [33] Bicak J., Ledvinka T., Schmidt B. G., Zofka M. Static [55] Brandenberger R., Firouzjahi H., Karouby J., Khosravi fluid cylinders and their fields: Global solutions. Class. K. Gravitational radiation by cosmic strings in a junc- Quantum Grav. v. 21, 1583-1608, 2004. tion. JCAP, v. 01, 008, 2009. [34] Bicak J., Schmidt B.G. Isometries compatible with [56] Brando G., Fabris J.C., Falciano F.T., Galkina O. Stiff gravitational radiation. J. Math. Phys. v. 25, 600-606, matter solution in Brans-Dicke theory and the general 1984. relativity limit. arXiv: 1810.07860, 2018. [35] Bieri L., Garfinkle D., Yunes N. Gravitational waves [57] Brito I., Carot J., Mena F.C., Vaz E.G.L.R. Cylindri- 45

cally symmetric static solutions of the Einstein field wormholes and energy conditions. Int. J. Mod. Phys. equations for elastic matter. J. Math. Phys. v. 53, v. A31, 1641022, 2016. 122504, 2012. [80] Bronnikov K.A., Krechet V.G. Potentially observable [58] Brito I., da Silva M.F.A., Mena F.C., Santos N.O. Con- cylindrical wormholes without exotic matter in GR. formally flat sources for the Linet-Tian spacetime. Gen. Phys. Rev. v, D99, 084051, 2019. Relativ. Grav. v. 45, 519-530, 2013. [81] Bronnikov K.A., Lemos J.P.S. Cylindrical wormholes. [59] Brito I., da Silva M. F. A., Mena F. C., Santos N. O. Phys. Rev. v. D79, 104019, 2009. Geodesics dynamics in the Linet-Tian spacetime with [82] Bronnikov K.A., Rubin S.G. Black holes, cosmology and Lambda < 0. Gen. Rel. Grav. v. 46, 1681, 2014. extra dimensions (World Scientific, Singapore, 2012). [60] Brito I., da Silva M. F. A., Mena F. C., Santos N. O. [83] Bronnikov K.A., Shikin G.N., Sibileva E.N. Self- Geodesics dynamics in the Linet-Tian spacetime with gravitating string-like configurations from nonlinear Λ > 0. Class. Quantum Grav. v. 32, 185015, 2015. electrodynamics. Grav. Cosmol. v. 9, 169-175, 2003. [61] Brito I., da Silva M.F.A., Mena F.C., Santos N. O. [84] Bronnikov K.A., Skvortsova M.V. Wormholes leading to Cylindrically symmetric inhomogeneous dust collapse extra dimensions. Grav. Cosmol. v. 22, 316-322, 2016. with a zero expansion component. Class. Quantum [85] Cai R.-G., Yang R.-Q. Multiple critical gravita- Grav. v. 34, 205005, 2017. tional collapse of charged scalar with reflecting wall. [62] Bronnikov K.A. Static, cylindrically symmetric arXiv:1602.00112, 2016. Einstein-Maxwell fields, in “Problems in gravita- [86] Cai R.-G., Zhang Y.-Z. Black plane solutions in four- tion theory and particle theory” (PGTPT), Ed. dimensional spacetimes. Phys. Rev. v. D54, 4891-4898, Staniukovich K.P., 10th issue, p. 37-50, Atomizdat, 1996. Moscow, 1979 (in Russian). [87] Caldwell R. R., Allen B. Cosmological constraints on [63] Bronnikov K.A. Inverted black holes and anisotropic cosmic-string gravitational radiation. Phys. Rev. v. D45, collapse. Russ. Phys. J. v. 22, 6, 594-600, 1979. 3447-3468, 1992. [64] Bronnikov K. A. Static fluid cylinders and plane layers [88] Campbell S. J., Wainwright J. Algebraic computing and in general relativity. J. Phys. A: Math. Gen. v. 12, 201- the Newman-Penrose formalism in general relativity. 207, 1979. Gen. Relativ. Grav. v. 8, 987-1001, 1977. [65] Bronnikov K.A. Gravitational and sound waves in stiff [89] Canfora F., Dimakis N., Paliathanasis A. Topologically matter. J. Phys. A: Math. Gen. v. 13, 3455-3463, 1980. nontrivial configurations in the 4d Einstein’s nonlinear [66] Bronnikov K.A. Rotating cylindrical wormholes: A no- σ-model system. Phys. Rev. v. D96, 025021, 2017. go theorem. J. Phys. Conf. Series v. 675, 012028, 2016. [90] Capozziello S., Corda, Ch., De Laurentis M. Stochas- [67] Bronnikov K.A. Scalar fields as sources for wormholes tic background of relic scalar gravitational waves from and regular black holes. Particles v. 2018, 1, 56-81, 2018. scalar-tensor gravity. Mod. Phys. Lett. v. A22, 2647- [68] Bronnikov K. A., Abdel-Sattar V., Chudaeva E. N., 2655, 2007. Shikin G. N. Static, cylindrically symmetric perfect fluid [91] Caprini C., Figueroa D.G. Cosmological Backgrounds of configurations. Vestnik RUDN, No. 1, 85-95, 2009. Gravitational Waves. Class. Quant. Grav. v. 35, 163001, [69] Bronnikov K.A., Baleevskikh K.A. On gravitational 2018. lensing by symmetric and asymmetric wormholes. Grav. [92] Carot J., Senovilla J.M.M. Vera R., On the definition Cosmol. v. 25, 2019 (to appear) [arXiv: 1812.05704]. of cylindrical symmetry. Class. Quantum Grav. v. 16, [70] Bronnikov K.A., Bolokhov S.V., Skvortsova M.V. Cylin- 3025-3034, 1999. drical wormholes: a search for viable phantom-free mod- [93] Carter B. The commutation property of a stationary, els in GR. In preparation, 2019. axisymmetric system. Commun. Math. Phys. v. 17, 233- [71] Bronnikov K.A., Fabris J.C. Weyl spacetimes and 238, 1970. wormholes in D-dimensional Einstein and dilaton grav- [94] Celestino J., de Oliveira H.P., Rodrigues E.L. Nonlin- ity. Class. Quantum Grav. v. 14, 831-842, 1997. ear evolution of cylindrical gravitational waves: Numer- [72] Bronnikov K. A., Fabris J. C., Silveira R., Zaslavskii ical method and physical aspects. Phys. Rev. v. D93, O. B. Dilaton gravity, (quasi-) black holes, and scalar 104018, 2016. charge. Gen. Relativ. Grav. v. 46, 1775, 2014. [95] Chakraborty S., Chakraborty S. Gravitational collapse [73] Bronnikov K. A., Fabris J. C., Silveira R., Zaslavskii of cylindrical anisotropic fluid: a source of gravitational O. B. Dilaton gravity, charged dust, and (quasi-) black waves. Gen. Relativ. Grav. v. 46, 1784, 2014. holes. Phys. Rev. v. D89, 107501, 2014. [96] Chakraborty S., Chakraborty S. A study of charged [74] Bronnikov K.A., Galiakhmetov A.M. Wormholes with- cylindrical gravitational collapse with dissipative fluid. out exotic matter in Einstein-Cartan theory. Grav. Cos- Adv. High Energy Phys. v. 2017, 8786791, 2017. mol. v. 21, 283-288, 2015. [97] Chandrasekhar S., Wright J.P. THE GEODESICS IN [75] Bronnikov K.A., Galiakhmetov A.M. Wormholes and GODEL’S UNIVERSE. Proc. Natl. Acad. Sci. v. 47, black universes without phantom fields in Einstein- 341, 1961. Cartan theory. Phys. Rev. v. D94, 124006, 2016. [98] Chernoff D.F., Flanagan E.E., Wardell B. Gravita- [76] Bronnikov K.A., Kim S.-W. Possible wormholes in a tional back-reaction on a cosmic string: Formalism. brane world. Phys. Rev. v. D67, 064027, 2003. arXiv:1808.08631, 2018. [77] Bronnikov K. A., Kovalchuk M. A. Properties of static [99] Chernoff D. F., Tye S. H. H. No Access Inflation, string fluid cylinders and plane layers in general relativity. theory and cosmic strings. Int. J. Mod. Phys. v. D24, Gen. Rel. Grav. v. 11, 343-355, 1979. 1530010, 2015. [78] Bronnikov K.A., Krechet V.G., Lemos J.P.S. Rotating [100] Chernoff D. F., Tye S. H. H. Detection of low tension cylindrical wormholes. Phys. Rev. v. D87, 084060, 2013. cosmic superstrings. JCAP, v. 05, 002, 2018. [79] Bronnikov K.A., Krechet V.G. Rotating cylindrical [101] Chicone C., Mashhoon B. Gravitomagnetic accelerators. 46

Phys. Lett. v. A375, 957-960, 2011. class. Class. Quantum Grav. v. 12, 111-118, 1995. [102] Chicone C., Mashhoon B. Gravitomagnetic jets. Phys. [123] da Silva M. F. A., Herrera L., Santos N. O., Wang A. Rev. v. D83, 064013, 2011. Rotating cylindrical shell source for Lewis spacetime. [103] Choptuik M. W. Universality and scaling in gravita- Class. Quantum Grav. v. 19, 3809-3819, 2002. tional collapse of a massless scalar field. Phys. Rev. Lett. [124] da Silva M.F.A., Wang A., Paiva F.M., Santos N.O. On v. 70, 9-12, 1993. the Levi-Civita solutions with cosmological constant. [104] Christodoulou D. Nonlinear nature of gravitation and Phys. Rev. v. D61, 044003, 2000. gravitational wave experiments. Phys. Rev. Lett. v. 67, [125] da Silva M. F. A., Wang A., Santos N. O. On the 1486-1489, 1991. sources of static plane symmetric vacuum spacetimes. [105] Cl´ement G., Fabbri A. Analytical treatment of crit- Phys. Lett. v. A244, 462-466, 1998. ical collapse in 2+1 dimensional AdS spacetime: a [126] de Araujo J.C.N., Wang A. Rigidly rotating dust in toy model. Class. Quantum Grav. v. 18, 3665-3680, general relativity. Gen. Relativ. Grav. v. 32, 1971-1980, 2001; Critical collapse in 2+1 dimensional AdS space- 2000. time: quasi-CSS solutions and linear perturbations. [127] Debbasch F., Herrera L., Pereira P. R. C. T., Santos Nucl. Phys. v. B630, 269-292, 2002. N. O. Stationary cylindrical anisotropic fluid. Gen. Rel. [106] Condron E., Nolan B.C. Collapse of a self-similar cylin- Grav. v. 38, 1825-1838, 2006. drical scalar field with non-minimal coupling. Class. [128] D’Inverno R. Introducing einstein’s relativity. Claren- Quantum Grav. v. 31, 015015, 2014. don Press. Oxford, 2003. [107] Condron E. J., Nolan B.C. Collapse of a self-similar [129] Di Prisco A., Herrera L., MacCallumM. A. H., Santos cylindrical scalar field with non-minimal coupling II: N. O. Shearfree cylindrical gravitational collapse. Phys. strong cosmic censorship. Class. Quantum Grav. v. 31, Rev. v. D80, 064031, 2009. 165018, 2014. [130] Doroshkevich A., Hansen J., Novikov I.D., Shatskiy A. [108] Copeland E. J., Kibble T. W. B. Cosmic strings and Passage of radiation through wormholes. Int. J. Mod. superstrings. Proc. Roy. Soc. Lond. v. A466, 623-657, Phys. v. D18, 1665-1691, 2009. 2010. [131] Dowker J. S. A gravitational Aharonov-Bohm effect. Il [109] Copeland E.J., Myers R.C., Polchinski J. Cosmic F and Nuovo Cimento, v. B52, 129-135, 1967. D strings. JHEP, v. 06, 013, 2004. [132] Dvali G., Vilenkin A. Formation and evolution of cosmic [110] Corda, C. The Virgo - MiniGRAIL cross correlation for D strings. JCAP, v. 03, 010, 2004. the detection of scalar gravitational waves. Mod. Phys. [133] Dowker J. S. A gravitational Aharonov-Bohm effect. Il Lett. v. A22, 1727-1735, 2007. Nuovo Cimento, v. B52, 129-135, 1967. [111] Corda, C. Interferometric Detection of Gravitational [134] East W.E. Cosmic Censorship Triumphant in waves: The Definitive Test for General Relativity. In- Spheroidal Collapse of Collisionless Matter. Phys. ter. J. Mod. Phys. v. D18, 2275-2282, 2009. Rev. Lett. v. 122, 231103, 2019. [112] Corda, Ch., Capozziello S., De Laurentis M. Tuning the [135] Echeverria F. Gravitational collapse of an infinite, cylin- stochastic background of gravitational waves with the- drical dust shell. Phys. Rev. v. D47, 2271-2282, 1993. ory and observations. AIP Conf. Proc. v. 966, 257-263, [136] Einstein A. Die feldgleichungen der gravitation. Sitz. 2008. Preuss. Akad. Wiss. 844-847, 1915. [113] Cornish N., Robson T. Galactic binary science with the [137] Einstein A., Rosen R. On gravitational waves. J. new LISA design. J. Phys. Conf. Ser. 840, 012024, 2017. Franklin Inst. v. 223, 43-54, 1937. [114] Cvetic M., Griffies S., Soleng H.H. Local and global [138] Eiroa E.F., de Celis E.R., Simeone C. Thin shells joining gravitational aspects of domain wall spacetimes. Phys. local cosmic string geometries. Eur. Phys. J. C v. 76, Rev. v. D48, 2613-2634, 1993. 546, 2016. [115] Damour T., Vilenkin A. Gravitational wave bursts from [139] Erez G., Rosen N. The gravitational field of a particle cosmic strings. Phys. Rev. Lett. v. 85, 3761-3764, 2000. possessing a multipole moment. Bull. Res. Counc. of [116] Damour T., Vilenkin A. Phys. Gravitational wave Israel, v. 8F, 47-50, 1959. bursts from cusps and kinks on cosmic strings. Phys. [140] Evans A. B. Static fluid cylinders in general relativity. Rev. v. D64, 064008, 2001. J. Phys. A: Math. Gen. v. 10, 1303-1311, 1977. [117] Damour T., Vilenkin A. Gravitational radiation from [141] Farr W.M., et al. The Mass Distribution of Stellar-Mass cosmic (super)strings: Bursts, stochastic background, Black Holes. Astrophys. J. v. 741, 103, 2011. and observational windows. Phys. Rev. v. D71, 063510, [142] Favata M. The gravitational-wave memory effect. Class. 2005. Quantum Grav. v. 27, 084036, 2010. [118] Danzmann K., et al.. Laser interferometer space an- [143] Ferrari A. MODELING EXTRAGALACTIC JETS. tenna. arXiv:1702.00786, 2017. Annu. Rev. Astron. Astrophys. v. 36, 539-98, 1998. [119] Darmois E. Les ´equationsde la gravitation einsteini- [144] Fj¨allborg M. Static cylindrically symmetric spacetimes. enne. M´emorialdes sciences math´ematiquesFacicule 25. Class. Quantum Grav. v. 24, 2253-2270, 2007. Gauttier-Villars. Paris 1927. [145] Flamm L. Beitr¨agezur Einsteinschen gravitationstheo- [120] da Silva M. F. A., Herrera L., Paiva F. M., Santos N. O. rie. Phys. Z. v. 17, 448-454, 1916. The Levi-Civita spacetime. J. Math. Phys. v. 36, 3625- [146] Forghani S. D., Mazharimousavi S. H., Halilsoy M. The 3631, 1995. stability of asymmetric cylindrical thin-shell wormholes. [121] da Silva M. F. A., Herrera L., Paiva F. M., Santos N. O. arXiv: 1807.05080, 2018. The parameters of the Lewis metric for the Weyl class. [147] Friedman J.L., Schleich K., Witt D.M. Topological cen- Gen. Rel. Grav. v. 27, 859-872, 1995. sorship. Phys. Rev. Lett. v. 71, 1486-1489 (1993); Erra- [122] da Silva M. F. A., Herrera L., Paiva F. M., Santos N. tum: ibid. v. 75, 1872, 1995. O. On the parameters of the Lewis metric for the Lewis [148] Friedman J.L, Higuchi A. Topological censorship and 47

chronology protection. Ann. Phys. v. 15, 109-128, 2006. [170] Griffiths J. B., Santos N. O. The rotation and trans- [149] Fryer C.L., et al. Compact Remnant Mass Function: lation of cylinders in General Relativity. Int. J. Mod. Dependence on the Explosion Mechanism and Metallic- Phys v. D19, 79-84, 2010. ity. Astrophys. J. v. 749, 91, 2012. [171] Grishchuk L.P., Polnarev A.G. Gravitational wave [150] Gararcia-Parrado A., Mena F.C. Gravitational radia- pulses with velocity memory. Sov. Phys. JETP, v. 69, tion and the evolution of gravitational collapse in cylin- 653-657 1989 [Zh. Eksp. Teor. Fiz. v. 96, 1153-1160, drical symmetry. arXiv:1811.10303, 2018. 1989]. [151] Garfinkle D. Exact solution for (2+1)-dimensional crit- [172] Guha S., Banerji R. Dissipative cylindrical collapse of ical collapse. Phys. Rev. v. D63, 044007, 2001. charged anisotropic fluid. Int. J. Theor. Phys. v. 53, [152] Garfinkle D., Gundlach C. Perturbations of an exact 2332-2348, 2014. solution for (2+1)-dimensional critical collapse. Phys. [173] Gundlach C., Martin-Garcia J. M. Critical phenomena Rev. v. D66, 044015, 2002. in gravitational collapse. Living Rev. Rel. v. 10, 5, 2007. [153] Gariel J., Marcilhacy G., Santos N. O. On the Lewis [174] Gundlach C., Baumgarte T.W. Critical gravitational metric. J. Math. Phys. v. 41, 7535-7543, 2000. collapse with angular momentum II: soft equations of [154] Gariel J., Marcilhacy G., Santos N. O. Unbound state. Phys. Rev. v. D97, 064006, 2018. geodesics from the ergosphere and the Messier 87 jet [175] T. Harko, Z. Kovacs, F.S.N. Lobo. Thin accretion disks profile. Astrophys. J. v. 774, 109, 2013. in stationary axisymmetric wormhole spacetimes. Phys. [155] Gariel J., Santos N. O., Silk J. Unbound geodesics from Rev. v. D79, 064001, 2009. the ergosphere and potential observability of debris from [176] Harte A.I. Strong lensing, plane gravitational waves and ultrahigh energy particle collisions. Phys. Rev. v. D90, transient flashes. Class. Quantum Grav. v. 30, 075011, 063505, 2014. 2013. [156] Gariel J., Santos N. O., Wang A. Kerr geodesics fol- [177] Hawking S. W., Ellis G. F. R. The large scale structure lowing the axis of symmetry. Gen. Rel. Grav. v. 48, 66, of spacetime. Cambridge University Press, Cambridge, 2016. England, 1973. [157] Gariel J., Santos N. O., Wang A. Observable accelera- [178] Hawking S.W., Perry M.J., Strominger A. Soft hair on tion of jets by a Kerr black hole. Gen. Rel. Grav. v. 49, black holes. Phys. Rev. Lett. v. 116, 231301, 2016; Su- 43, 2017. perrotation charge and supertranslation hair on Black [158] Gleiser R.J. Axial gravitational perturbations of an in- Holes, JHEP, v. 05, 161, 2017. finite static line source. Class. Quantum Grav. v. 32, [179] Hennig J., Neugebauer G. An analytical cylindri- 065003, 2015. cally symmetric solution for collapsing dust. arXiv:gr- [159] Gleiser R.J. Linear stability of the Linet-Tian solution qc/0409096, 2004. with negative cosmological constant. Class. Quantum [180] Herrera L., Di Prisco A., Ospino J., Fuenmayor E. Con- Grav. v. 34, 065010, 2017. formally flat anisotropic spheres in general relativity. J. [160] Gleiser R.J. Linear stability of the Linet - Tian solution Math. Phys. v. 42, 2129-2143, 2001. with positive cosmological constant. arXiv:1810.07296, [181] Herrera L., Le Denmat G., Marcilhacy G., Santos N.O. 2018. Static cylindrical symmetry and conformal flatness. Int. [161] G¨odelK. An example of a new cosmological solution J. Mod. Phys. v. D14, 657-666, 2005. of Einstein’s field equations of gravitation. Rev. Mod. [182] Herrera L., Paiva F. M., Santos N. O. The Levi-Civita Phys. v. 21, 447-450, 1949 [Gen. Rel. Grav. v. 32, 1409- spacetime as a limiting case of the γ spacetime. J. Math. 1417, 2000]. Phys. v. 40, 4064-4071, 1999. [162] Goldenfeld N. Lectures on phase transitions and the [183] Herrera L., Santos N. O. A general relativistic effect in renormalization group. Addison Wesley Publishing quasi-spherical objects as the possible origin of relativis- Company, New York, 1992. tic jets. Astrophys. Space Sci. v. 310, 251-254, 2007. [163] Gong Y.-G., Hou S.-Q., Liang D.-C., Papantonopoulos [184] Herrera L., Paiva F. M., Santos N.O. Geodesics in the E. Gravitational waves in Einstein-aether and general- γ spacetime. Int. J. Mod. Phys. v. D9, 649-659, 2000. ized TeVeS theory after GW170817. Phys. Rev. v. D97, [185] Herrera L., Santos N. O. Geodesics in Lewis spacetime. 084040, 2018. J. Math. Phys. v. 39, 3817-3827, 1998. [164] Gregory R. Nonsingular global strings. Phys. Rev. v. [186] Herrera L., Santos N.O., Teixeira A. F. F., Wang A. D54, 4955-4962, 1996. On the interpretation of cylindrically symmetric Levi- [165] Griffiths J. B. The rotation and translation of cylinders Civita spacetime for 0 ≤ σ < ∞. Class. Quantum Grav. in General Relativity. Int. J. Mod. Phys. v. D19, 79-84, v. 18, 3847-3855, 2001. 2010. [187] Herrera L., Santos N.O. Cylindrical collapse and gravi- [166] Griffiths J.B. Colliding plane waves in general relativity. tational waves. Class. Quantum Grav. v. 22, 2407-2413, Dover Publications, Inc. Mineola, New York, 2016. 2005. [167] Griffiths J. B., Podolsk´yJ. Exact spacetimes in Ein- [188] Herrera L., Ruifern´andezJ., Santos N. O. The static stein’s general relativity. Cambridge Monographs on cylinder, gyroscopes and the C-metric. Gen. Rel. Grav. Mathematical Physics. Cambridge University Press, v. 33, 515-529, 2001. 2009. [189] Hindmarsh M. B., Kibble T. W. B. Cosmic strings. Rep. [168] Griffiths J.B., Podolsky J. The Linet-Tian solution with Prog. Phys. v. 58, 477-562, 1995. a positive cosmological constant in four and higher di- [190] Hirschmann E.W., Wang A., Wu Y. Collapse of a scalar mensions. Phys. Rev. v. D81, 064015, 2010. field in 2 + 1 gravity. Class. Quant. Grav. v. 21, 1791- [169] Griffiths J. B., Santos N. O. A rotating cylinder in an 1824, 2004. asymptotically locally anti-de-Sitter background. Class. [191] Hochberg D, Visser M. Geometric structure of the Quantum Grav. v. 27, 125004, 2010. generic static traversable wormhole throat. Phys. Rev. 48

v. D56, 4745-4755, 1997. cosmological constant on trajectories of light and asso- [192] Horbatscha M.W., Burgess C.P. Cosmic black-hole hair ciated measurements in Schwarzschild de Sitter space. growth and quasar OJ287. JCAP, v. 05, 010, 2012. arXiv: 1308.4931 [gr-qc] [193] Ishibashi A., Wald R.M. Dynamics in non-globally- [216] Letelier P.S., Wang A. Singularities formed by the fo- hyperbolic static spacetimes III: Anti-de Sitter space- cusing of cylindrical null fluids. Phys. Rev. v. D49, 5105- time. Class. Quantum Grav. v. 21, 2981-3014, 2004. 5110, 1994; ibid., v. D51, 5968(E), 1995. [194] Israel W. Singular hypersurfaces and thin shells in gen- [217] Letelier P.S., Wang A. Spacc-time defects. J. Math. eral relativity. Nuovo Cimento. v. B44, 1-14, 1966; ibid. Phys. v. 36, 3023-3042, 1995. v. B48, 463(E), 1967. [218] Levi-Civita T. ds2 einsteiniani in campi newtoniani. IX: [195] Jacobson T. Primordial Black Hole Evolution in Tensor- L’analogo del potenziale logaritmico. Rend. Accad. Lin- Scalar Cosmology. Phys. Rev. Lett. v. 83, 2699, 1999. cei, v. 28, 101-109, 1919. [196] Jensen B., Ku˘cera J. On a gravitational Aharonov- [219] Lewis T. Some special solutions of the equations of ax- Bohm effect. J. Math. Phys. v. 34, 4975-4985, 1993. ially symmetric gravitational fields. Proc. R. Soc. v. [197] Jordan P., Ehlers J., Kundt W. Exact solutions of the A136, 176-192, 1932. field equations of general relativity. Akad. Wiss. Lit. [220] Lin K., Zhao X., Zhang C., Liu T., Wang B., Zhang S.- (Mainz) Abhandl. Math. Nat. Kl. Jahrg, v. 2, 21, 1960. J., Zhang X., Zhao W., Zhu T., Wang A. Gravitational [198] Joshi P.S. Global aspects in gravitation and cosmology. wave forms, polarizations, response functions and en- Clarendon, Oxford, 1993. ergy losses of triple systems in Einstein-Aether theory. [199] Kar S., Lahiri S., SenGupta S. Can extra dimensional Phys. Rev. v. D99, 023010, 2019. effects allow wormholes without exotic matter? Phys. [221] Linet B. The static metrics with cylindrical symme- Lett. v. B750, 319, 2015. try describing a model of cosmic strings. Gen. Relativ. [200] Khakshournia S. Gravitational collapse of a rotating Gravit. v. 17, 1109-1115, 1985. cylindrical null shell in the cosmic strong spacetime. In- [222] Liu T., Zhang X., Zhao W., Lin K., Zhang C., Zhang ter. J. Mod. Phys. Conf. Ser. 03, 428, 2011. S.-J., Zhao X., Zhu T., Wang, Waveforms of compact [201] Khorrami M., Mansouri R. Cylindrically symmetric thin binary inspiral gravitational radiation in screened mod- walls in general relativity. J. Math. Phys. v. 35, 951-958, ified gravity. Phys. Rev. v. D98, 083023, 2018. 1994. [223] Lobo, F.S.N. Exotic solutions in general relativity: [202] Kibble T.W.B. Topology of cosmic domains and strings. Traversable wormholes and “warp drive” spacetimes. In J. Phys. v. A9, 1387-1398, 1976. Classical and Quantum Gravity Research. Nova Science [203] Kibble T.W.B. Some implications of a cosmological Publishers: New York, NY, USA, 2008. p.1-78. phase transition. Phys. Report, v. 67, 183-199, 1980. [224] Lobo, F.S.N. (ed.) Wormholes, warp drives and energy [204] Kirillov A. A., Savelova E. P. On the value of the cos- conditions. Springer: Berlin, Germany, 2017. mological constant in a gas of virtual wormholes. Grav. [225] Luo J., et al. TianQin: a space-borne gravitational wave Cosmol. v. 19, 92-100, 2013. detector. Class. Quantum Grav. v. 33, 035010, 2016. [205] Klep´acP., J Horsk´yJ. Charged perfect fluid and scalar [226] MacCallum M. A. H., Santos N. O. Stationary and field coupled to gravity. Czech J Phys v. 51, 1177-1187, static cylindrically symmetric Einstein spaces of the 2001. Lewis form. Class. Quantum Grav. v. 15, 1627-1636, [206] Kompaneets A.S. Strong gravitational waves in vacuum, 1998. Sov. Phys. JETP, v. 7, 659-668, 1958 [Zh. Eks. Teor. Fiz. [227] Maeda H., Nozawa M. Static and symmetric wormholes v. 34, 953, 1958]. respecting energy conditions in Einstein-Gauss-Bonnet [207] Krasinski A. Solutions of the Einstein field equations gravity. Phys. Rev. v. D78, 024005, 2008. for a rotating perfect fluid II: properties of the flow- [228] Mahmood T., Shah S.M., Abbas G. Gravitational col- stationary and and vortex-homogeneous solutions. Acta lapse and expansion of charged anisotropic cylindrical Phys. Polon. v. 6, 223, 1975. source. Astrophys. Space Sci. v. 357, 56, 2015. [208] Krechet V.G. Izvestiya Vuzov. Fiz. v. No 10, 57-62, [229] Majumdar S. D. A class of exact solutions of Einstein’s 2007. field equations. Phys. Rev. v. 72, 390-398, 1947. [209] Krechet V.G., Sadovnikov D.V. Spin-spin interaction in [230] Maluf J.W., da Rocha-Neto J.F., Ulhoa S.C., Carneiro general relativity and induced geometries with nontriv- F.L. Plane gravitational waves, the kinetic energy of free ial topology. Grav. Cosmol. v. 15, 337-340, 2009. particles and the memory effect, arXiv1707.06874. [210] Krori K. D., Barua, J. Indian J. Pure Appl. Phys. v. 12, [231] Marder L. Gravitational waves in general relativity I. 818-822, 1974. Cylindrical waves. Proc. Roy. Soc. v. A244, 524-537, [211] Kundt W. Tr¨agheitsbahnen in einem von G¨odel 1958. angegebenen kosmologischen Modell. Z. Phys. v. 145, [232] Mars M., Senovilla J.M.M. Axial symmetry and confor- 611, 1956. mal Killing vectors. Class. Quantum Grav. v. 10, 1633- [212] Kurita Y., Nakao K.-i. Naked singularity resolution in 1647, 1993. cylindrical collapse. Phys. Rev. v. D73, 064022, 2006. [233] Mashhoon B., McClune J.C., Quevedo H. On a class of [213] Lanczos C. Uber¨ eine station¨arekosmologie in sinne der rotating gravitational waves. Class. Quantum Grav. v. Einsteinischen gravitationstheories. Z. Physik, v. 21, 73- 17, 533-549, 2000. 110, 1924. [234] Mashhoon B., Quevedo H. Rotating gravitational [214] Lasky P.D., Thrane E., Levin Y., Blackman J., Chen Y. waves. Phys. Lett. v. A151, 464-468, 1990. Detecting gravitational wave memory with LIGO: impli- [235] Mashhoon B., Santos N. O. Rotating cylindrical systems cations of GW150914, Phys. Rev. Lett. v. 117, 061102, and gravitomagnetism. Ann. Phys. (Berlin) v. 9, 49-63, 2016. 2000. [215] Lebedev D and Lake K 2013 On the influence of the [236] Melvin M.A. Pure magnetic and electric geons. Phys. 49

Lett. v. 8, 65-68, 1964. Geodesic motion and confinement in Lanczos spacetime. [237] Mena F.C. Cylindrical symmetric models of gravita- Class. Quantum Grav. v. 13, 1641-1654, 1996. tional collapse to black holes: a short review. Inter. J. [259] Pereira P.R.C.T., Wang A. Dynamics of rotating cylin- Mod. Phys. v. D24, 1542021, 2015. drical shells in general relativity. Gen. Relativ. Grav. v. [238] Miguelote A.Y., da Silva M.F.A., Wang A., Santos 32, 2189-2218, 2000. N.O. Levi-Civita solutions coupled with electromagnetic [260] Pereira P.R.C.T., Wang A. Collapsing and expanding fields. Class. Quantum Grav. v. 18, 4569-4588, 2001. cylindrically symmetric fields with light-like wave-fronts [239] Miranda G., Matos T. Exact rotating magnetic in general relativity. Inter. J. Mod. Phys. v. D11, 561- traversable wormholes satisfying the energy conditions. 579, 2002. arXiv: 1507.02348, 2015. [261] Pereira P.R.C.T., Wang A. Gravitational collapse of [240] Mishima T., Tomizawa S. Construction and applica- cylindrical shells made of counterrotating dust parti- tion of variations on the cylindrical gravitational waves cles. Phys. Rev. v. D62, 124001, 2000; ibid. v. D67, of Weber, Wheeler, and Bonnor. Phys. Rev. v. D96, 129902(E), 2003. 024023, 2017. [262] Philbin Thomas G. Perfect-fluid cylinders and walls - [241] Misner C.W., Thorne K.S., Wheeler J.A. Gravitation. sources for the Levi-Civita spacetime. Class. Quantum W.H. Greeman and Company. New York, 1973. Grav. v. 13, 1217-1232, 1996. [242] Morris M., Thorne K.S, Yurtsever U. Wormholes, time [263] Piran T., Safier P.N. A gravitational analogue of Fara- machines, and the Weak Energy Condition. Phys, Rev. day rotation. Nature, v. 318, 271-273, 1985. Lett. v. 61, 1446-1449, 1988. [264] Piran T., Safier P.N., Stark R.F. General numerical so- [243] Nakao K.-i., Harada T., Kurita Y., Morisawa Y. Rel- lution of cylindrical gravitational waves. Phys. Rev. v. ativistic gravitational collapse of a cylindrical shell of D32, 3101-3107, 1985. dust. II. Prog. Theor. Phys. v. 122, 521-541, 2009. [265] Qiao, J., Zhu T., Zhao W., Wang A. Waveform of grav- [244] Nakao K.-i., Ida D., Kurita Y. New instability in rel- itational waves in the ghost-free parity-violating gravi- ativistic cylindrically symmetric system. Phys. Rev. v. ties. arXiv:1909.03815, 2019. D77, 044021, 2008. [266] Quevedo H., Mashhoon B. Relativity and gravitation: [245] Nakao K.-i., Kurita Y., Morisawa Y., Harada T. Rel- Classical and quantum, Proc. 7th Conf. SILARG ed ativistic gravitational collapse of a cylindrical shell of D’Olivo J. C. et al, Singapore: World Scientific, 1991, dust. Prog. Theor. Phys. v. 117, 75-102, 2007. pp.287-93. [246] Nakao K.-i., Morisawa Y. High speed dynamics of col- [267] Regimbau T., Giampanis S., Siemens X., Mandic lapsing cylindrical dust fluid. Class. Quantum Grav. v. V. Stochastic background from cosmic (super)strings: 21, 2101-2113, 2004. Popcorn-like and (Gaussian) continuous regimes. Phys. [247] Nakao K.-i., Morisawa Y. High-speed cylindrical col- Rev. v. D85, 066001, 2012. lapse of perfect fluid. Prog. Theor. Phys. v. 113, 73-85, [268] Rindler W. Essential relativity, 2nd ed, Springer Verlag, 2005. Berlin, 1977, pp.156-164. [248] Newman E., Penrose R. An approach to gravitational [269] Ringeval C., Suyama T. Stochastic gravitational waves radiation by a method of spin coefficients. J. Math. from cosmic string loops in scaling. JCAP, v. 12, 027, Phys. v. 3, 566-578, 1962; ibid., v. 4, 998(E), 1963. 2017. [249] Nogales J.A.C., Wang A. Instability of cosmological [270] Rovelli C. Reality is not what it seems. Penguin Random event horizons of nonstatic global cosmic strings. II. House UK, 2017. Perturbations of gravitational waves and massless scalar [271] Ruan W.-H., Guo Z.-K., Cai R.-G., Zhang Y.- fields. Phys. Rev. v. D57, 6089-6093, 1998. Z. Taiji Program: Gravitational-Wave Sources. [250] Nolan B.C. Naked singularities in cylindrical collapse of arXiv:1807.09495, 2018. counterrotating dust shells. Phys. Rev. v. D65, 104006, [272] Safko J., Witten L. Some properties of cylindrically 2002. symmetric Einstein-Maxwell fileds. J. Math. Phys. v. [251] Odintsov S.D., Oikonomou V.K. The early-time cos- 12, 257-293, 1971. mology with stiff era from modified gravity. arXiv: [273] Safko J., Witten L. Some properties of cylindrically 1711.04571, 2017. symmetric gravitational field. Phys. Rev. v. D5, 293- [252] Olmez S., Mandic V., Siemens X. Gravitational-wave 300, 1972. stochastic background from kinks and cusps on cosmic [274] Sakellariadou M. Gravitational waves emitted from infi- strings. Phys. Rev. v. D81, 104028, 2010. nite strings. Phys. Rev. v. D42, 354-360, 1990 [Erratum- [253] Opher R., Santos N. O., Wang A. Geodesic motion and ibid. v. D43, 4150, 1991]. confinement in van Stockum spacetime. J. Math. Phys. [275] Santos N. O. Solution of the vacuum Einstein equa- v. 37, 1982-1990, 1996. tions with non-zero cosmological constant for a station- [254] Ozel F., et al. The Black Hole Mass Distribution in the ary cylindrically symmetric spacetime. Class. Quantum Galaxy. Astrophys. J. v. 725, 1918, 2010. Grav. v. 10, 2401-2406, 1993. [255] Pacheco J. A., Gariel J., Marcilhacy G., Santos N. O. [276] Schwarzschild K. Uber¨ das Gravitationsfeld eines Gravitational model of high-energy particles in a colli- Massenpunktes nach der Einsteinschen Theorie. Sitz. mated jet. Astrophys. J. v. 759, 125, 2012. Preuss. Akad. Wiss. 189-196, 1916. [256] Papapetrou A. A static solution of the gravitational field [277] Setare M. R., Sepehri A. Stability of cylindrical thin equations for an arbitrary charge distribution. Proc Roy shell wormhole during evolution of universe from infla- Irish Acad. v. A51, 191-204, 1947. tion to late time acceleration. arXiv:1412.8666, 2014. [257] Penrose R. Gravitational collapse: The role of general [278] Sesana A., et al. Observing gravitational waves from relativity. Riv. Nuovo Cimento, v. 1, 252-276, 1969. the first generation of black holes. Astrophys. J. Lett. [258] Pereira P. R. C. T., N O Santos, N.O., Wang A. v. 698, L129, 2009. 50

[279] Sesana A. Prospects for Multi-band Gravitational-Wave tromagnetic fields in relativity, J. Math. Phys. v. 15, Astronomy after GW150914. Phys. Rev. Lett. v. 116, 1756-1759, 1974. 231102, 2016. [302] Thorne K.S. in Magic without magic: John Archibald [280] Shaha S. M., Abbas G. Dynamics of charged bulk vis- Wheeler. ed. J. Klauder, Friemann, San Francisco, p. cous collapsing cylindrical source with heat flux. Eur. 231-258, 1972. Phys. J. v. C77, 251, 2017. [303] Thorne K.S. Gravitational-wave bursts with memory: [281] Shaikh R., and S. Kar S. Wormholes, the weak energy the Christodoulou effect. Phys. Rev. v. D45, 520-524, condition, and scalar-tensor gravity, Phys. Rev. v. D94, 1992. 024011, 2016. [304] Tian Q.J. Cosmic strings with cosmological constant. [282] Shapiro S.L., Teukolsky S.A. Formation of naked singu- Phys. Rev. v. D33, 3549-3555, 1986. larities: The violation of cosmic censorship. Phys. Rev. [305] Tomizawa S., Mishima T. New cylindrical gravitational Lett. v. 66, 994-997, 1991. soliton waves and gravitational Faraday rotation, Phys. [283] Shapiro S.L., Teukolsky S.A. Gravitational collapse of Rev. v. D90, 044036, 2014. rotating spheroids and the formation of naked singular- [306] Tye S. H. H. Brane inflation : String theory viewed from ities. Phys. Rev. v. D45, 2006-2012, 1992. the cosmos. Lect. Notes Phys. v. 737, 949-974, 2008. [284] Sharif M., Abbas G. Dynamics of non-adiabatic charged [307] Vachaspati T., Vilenkin A. Gravitational radiation from cylindrical gravitational collapse. Astrophys. Space Sci, cosmic strings. Phys. Rev. v. D31, 3052-3058, 1985. v. 335, 515-521, 2011. [308] van Stockum W. J. The gravitational field of a distribu- [285] Sharif M., Abbas G. Charged perfect fluid cylindri- tion of particles rotating about about an axis of symme- cal gravitational collapse. J. Phys. Soc. Japan, v. 80, try. Proc. Roy. Soc. Edinburgh, v. A57, 135-154, 1937. 104002, 2011. [309] Vilenkin A., Shellard E. P. S. Cosmic strings and other [286] Sharif M., Ahmad Z. High-speed cylindrical collapse of topological defects. Cambridge University Press, Cam- two perfect fluids. Gen. Relativ. Gravit. v. 39, 1331- bridge, 2000. 1344, 2007. [310] Visser M. Lorentzian wormholes: from Einstein to [287] Sharif M., Azam M. Effects of electromagnetic field on Hawking. AIP, Woodbury, 1995. the dynamical instability of cylindrical collapse. JCAP, [311] Wainwright J. Exact spatially inhomogeneous cosmolo- v. 02, 043, 2012. gies. J. Phys. A: Math. Gen. v. 14, 1131-1147, 1981. [288] Sharif M., Fatima S. Charged cylindrical collapse of [312] Wald R.W., Iyer V. Trapped surfaces in the anisotropic fluid. Gen. Relativ. Gravit. v. 43, 127-142, Schwarzschild geometry and cosmic censorship, Phys. 2011. Rev. v. D44, R3719-R3722, 1991. [289] Smoot G.F. et al. (COBE). Structure in the COBE dif- [313] Wang A. Interacting, gravitational, electromagnetic, ferential microwave radiometer first-year maps. Astro- neutrino and other waves in the context of Einstein’s phys. J. v. 396, L1-5, 1992. general theory of relativity, A dissertation submitted [290] Souriau J.-M. Ondes et radiations gravitationnelles in: to Physics Department in partial fulfillment of the re- Colloques Internationaux du CNRS, Paris, v. 220, 243- quirements for the degree of Doctor of Philosophy, the 256, 1973. University of Ioannina, Greece, 1991. [291] Spergel D.N. et al. (WMAP). Wilkinson Microwave [314] Wang A. Gravitational Faraday rotation induced from Anisotropy Probe (WMAP) Three Year Results: Im- interacting gravitational plane waves. Phys. Rev. v. plications for Cosmology. Astrophys. J. Suppl. v. 170, D44, 1120-1131, 1991. 377-408, 2007. [315] Wang A. The effect of polarization of colliding plane [292] Stachel J. Globally stationary but locally static space- gravitational waves on focusing singularities. Inter. J. times: A gravitational analog of the Aharonov-Bohm Mod. Phys. v. A6, 2273-2288, 1991. effect. Phys. Rev. v. D26, 1281-1290, 1982. [316] Wang A. Critical phenomena in gravitational collapse: [293] Stephani H. General relativity. Cambridge University The studies so far. Braz. J. Phys. v. 31, 188-197, 2001 Press, 1990, p.175. [arXiv:gr-qc/0104073]. [294] Stephani H., Kramer D., MacCallum M. A. H., Hoense- [317] Wang A. Critical collapse of cylindrically symmetric laers C., Herlt E. Exact solutions of Einstein’s field scalar field in four-dimensional Einstein’s theory of grav- equations. Cambridge Monographs on Mathematical ity. Phys. Rev. v. D68, 064006, 2003. Physics. Cambridge University Press, 2009. [318] Wang A. Comment on “absence of trapped surfaces and [295] Stott M.J., Elghozi T., Sakellariadou M. Gravitational singularities in cylindrical collapse”. Phys. Rev. v. D72, wave bursts from cosmic string cusps and pseudocusps. 108501, 2005. Phys. Rev. v. D96, 023533, 2017. [319] Wang A. No-go theorem in spacetimes with two com- [296] Strominger A. Lectures on the infrared structure of muting spacelike Killing vectors. Gen. Relativ. Gravit. gravity and gauge theory. arXiv:1703.05448, 2017. v. 37, 1919-1926, 2005. [297] Synge J. L. Relativity: the general theory. North- [320] Wang A. Interacting, gravitational, electromagnetic, Holland Publishing Company, Amsterdam, 1960, p.313. neutrino and other waves in the context of Einstein’s [298] Szekeres P. The gravitational compass. J. Math. Phys. general theory of relativity, World Scientific, 2019 v. 6, 1387-1391, 1965; ibid.,v. 7, 751(E), 1966. (ISBN: 978-981-121-148-5). [299] Taub A.H. Empty spacetimes admitting a three param- [321] Wang A., da Silva M.F.A., Santos N.O. On parameters eter group of motions. Ann. Math. v. 53, 472-490, 1951. of the Levi-Civita solution. Class. Quantum Grav. v. 14, [300] Teixeira A. F., Wolk. I., Som M. M. Exact relativistic 2417-2423, 1997. cylindrical solution of disordered radiation. Nuovo Ci- [322] Wang A., Letelier P.S. Dynamical wormholes and en- mento, v. B41, 387-396, 1977. ergy conditions. Prog. Theor. Phys. 94, L137-142, 1995. [301] Teixeira A.F., Wolk I., Som M. Generalized static elec- [323] Wang A., Nogales J.A.C. Instability of cosmological 51

event horizons of nonstatic global cosmic strings. Phys. lapse of massless scalar field in f(R) gravity. Phys. Rev. Rev. v. D56, 6217-6224, 1997. v. D94, 104013, 2016. [324] Wang A., Santos N. O. Gravitational and particle radi- [337] Zhang P.-M., Duval C., Gibbons G.W., Horvathy P.A. ation from long cosmic strings. Class. Quantum Grav. The memory effect for plane gravitational waves. Phys. v. 13, 715-722, 1996. Lett. v. B772, 743-746, 2017. [325] Wang T.-Z., Fier J., Li B.-W., L¨uG.-L., Wang Z.-J., Wu [338] Zhang P.-M., Duval C., Gibbons G.W., Horvathy Y., Wang A. Singularities of plane gravitational waves P.A. Velocity memory effect for polarized gravitational and their memory effects. arXiv:1807.09397, 2018. waves. JCAP, v. 05, 030, 2018. [326] Wang, T.-Z., Li, L., Zhu, C.-H., Wang Z.-J., Wang A., [339] Zhang X., Yu J.-M., Liu T., Zhao W., Wang A. Testing Wu Q., Liu H.-S., Lu G.-L. An Alternative Channel For Brans-Dicke gravity using the Einstein telescope. Phys. High-Mass Binary Black Holes - Dark Matter Accretion Rev. v. D95, 124008, 2017. Onto Black Hole. Astrophys. J. v. 863, 17, 2018. [340] Zhang X., Zhao W., Liu T., Lin K., Zhang C., Zhao [327] Weber J., Wheeler J.A. Reality of the cylindrical gravi- X., Zhang S.-J. Zhu T., Wang A. Angular momentum tational waves of Einstein and Rosen. Rev. Mod. Phys. loss for eccentric compact binary in screened modified v. 29, 509-515, 1957. gravity. JCAP, v. 01, 019, 2019. [328] Weyl H. Bemerkung ¨uber die axisymmetrischen [341] Zhang X., Zhao W., Liu T., Lin K., Zhang C., Zhao L¨osungen der Einsteinschen gravitationsgleichungen. X., Zhang S.-J. Zhu T., Wang A. Constraints of general Ann. Physik, v. 59 (364), 185-188, 1919. screened modified gravities from comprehensive analysis [329] Wheeler J.A. On the nature of quantum geometrody- of binary pulsars. Astrophys. J. v. 872, 121, 2019. namics. Ann. Phys. v. 2, 604-614, 1957. [342] Zhang X.-F., L¨uH. Critical behavior in a massless scalar [330] Whittaker E. T. On Gauss’ theorem and the concept of field collapse with self-interaction potential. Phys. Rev. mass in general relativity. Proc. Roy. Soc. London, v. v. D91, 044046, 2015. A149, 384-395, 1935. [343] Zhao, W., Liu T., Wen L.-Q., Zhu T., Wang A., Hu Q., [331] Wright J. P. Solutions of Einstein’s field equations for Zhou C. Model-independent test of the parity symmetry a rotating, stationary and dust-filled universe. J. Math. of gravity with gravitational waves. arXiv:1909.13007, Phys. v. 6, 103-105, 1965. 2019. [332] Wu Y., da Silva M. F. A., Santos N. O., Wang A. Topo- [344] Zhao W., Zhu T., Qiao J., Wang A. Waveform of grav- logical charged black holes in high dimensional space- itational waves in the general parity-violating gravities. times and their formation from gravitational collapse of arXiv:1909.10887, 2019. a type II fluid. Phys. Rev. v. D68, 084012, 2003. [345] Zhao X., Zhang C., Lin K., Liu T., Niu R., Wang B., [333] Yagi K., Stein L.C., Yunes N., Tanaka T. Post- Zhang S.-J., Zhang X., Zhao W., Zhu T., Wang A. Grav- Newtonian, quasicircular binary inspirals in quadratic itational waveforms and radiation powers of the triple modified gravity. Phys. Rev. v. D85, 064022, 2012. system PSR J0337+1715 in modified theories of gravity. [334] Yunes N., Yagi K., Pretorius F. Theoretical physics im- Phys. Rev. v. D100, 083012, 2019. plications of the binary black-hole mergers GW150914 [346] Zipoy D. Topology of some spheroidal metrics. J. Math. and GW151226. Phys. Rev. v. D94, 084002, 2016. Phys. v. 7, 1137-1143, 1966. [335] Zeldovitch Ya. B., Polnarev A.G. Radiation of gravita- [347] Zofka M., Bicak J. Cylindrical spacetimes with Λ 6= 0 tional waves by a cluster of superdense stars. Sov. As- and their sources. Class. Quantum Grav. v. 25, 015011, tron. v. 18, 17-23, 1974 [Astron. Zh. v. 51, 30, 1974]. 2008. [336] Zhang C.-Y., Tang Z.-Y., Wang B. Gravitational col-