Jan 18, 2018 - as in the case of Harrison's TGW solution. The motion is periodic with turning points at. Â±W0, where W0 > 0 and Â±W0 are the only real...

1 downloads 7 Views 366KB Size

per le Applicazioni del Calcolo “M. Picone”, CNR, I-00185 Rome, Italy

arXiv:1801.06003v2 [gr-qc] 10 Feb 2018

2 ICRANet,

Piazza della Repubblica 10,

I-65122 Pescara, Italy 3 Department

of Mathematics,

University of Missouri, Columbia, Missouri 65211, USA 4 Department

of Physics and Astronomy,

University of Missouri, Columbia, Missouri 65211, USA 5 School

of Astronomy,

Institute for Research in Fundamental Sciences (IPM), P. O. Box 19395-5531, Tehran, Iran

(Dated: February 13, 2018) In general relativity (GR), linearized gravitational waves propagating in empty Minkowski spacetime along a fixed spatial direction have the property that the wave front is the Euclidean plane. Beyond the linear regime, exact plane waves in GR have been studied theoretically for a long time and many exact vacuum solutions of the gravitational field equations are known that represent plane gravitational waves. These have parallel rays and uniform wave fronts. It turns out, however, that GR also admits exact solutions representing gravitational waves propagating along a fixed direction that are nonplanar. The wave front is then nonuniform and the bundle of rays is twisted. We find a class of solutions representing nonplanar unidirectional gravitational waves and study some of the properties of these twisted waves. PACS numbers: 04.20.Cv, 04.30.Nk

2 Keywords: General Relativity, Exact Gravitational Waves

I.

INTRODUCTION

It is a natural consequence of a field theory of universal gravitation that gravitational radiation should exist in some form. The recent detection of gravitational waves is consistent with Einstein’s theory of gravitation; see, for instance, Ref. [1] and the references cited therein. This circumstance provides further motivation to study exact gravitational waves within the context of general relativity (GR). Indeed, exact plane wave solutions of GR have received much attention [2–6]. Many exact solutions also exist for the case of two colliding gravitational plane waves, as summarized in Ref. [5], where the main result seems to be the formation of a spacetime singularity as a consequence of the collision itself, i.e. of the mutual interaction and focusing of the two waves. Finally, special solutions corresponding to impulsive plane waves with Dirac delta-function profiles have also been studied [7]. A characteristic feature of plane waves in general relativity is that they admit a covariantly constant null vector field k, kµ;ν = 0 ,

kµ k µ = 0 ,

(1)

where k is the Killing vector field representing the propagation vector of the plane wave [3]. Let us recall that any gravitational wave propagating along a fixed direction in space—say, the z direction—always admits a null Killing propagation vector k = ∂t + ∂z , kµ k µ = 0 ,

kµ;ν + kν;µ = 0 ,

kµ;ν k ν = 0 ,

(2)

which therefore represents a non-expanding and shear-free null geodesic congruence. Furthermore, if the congruence is twist-free as well, then we recover Eq. (1) and we have a plane gravitational wave. It is however possible that the congruence is such that Tµν = k[µ;ν] = k[µ,ν] 6= 0 ,

(3)

in which case we have a twisted gravitational wave, namely, an exact unidirectional nonplanar gravitational wave in GR. Here, Tµν is the antisymmetric twist tensor. The wave ∗ † ‡

Electronic address: [email protected] Electronic address: [email protected] Electronic address: [email protected]

3 propagation vector k is orthogonal to the wave front, which are surfaces of constant u. Thus, ω 2 :=

1 1 k[µ;ν] k [µ;ν] = Tµν Tµν = 0 , 2 2

(4)

which is consistent with the notion that the twist scalar ω vanishes if and only if the null geodesic congruence is hypersurface-orthogonal [3]. Therefore, kµ = ∂u/∂xµ for plane gravitational waves, while kµ = Ω(x) ∂u/∂xµ with ∂Ω/∂xµ 6= 0 for nonplanar gravitational waves under consideration in this paper. An interesting solution which can be associated with a nonplanar gravitational wave is due to B. K. Harrison [8]. We encountered this solution within the context of cosmic jets and studied the peculiar asymptotic behavior of its timelike geodesics in our previous work [9]. The purpose of the present paper is to present a class of exact solutions of GR corresponding to twisted gravitational waves, i.e. waves propagating along a given direction (the z axis, say) with wave fronts that are nonplanar and hence nonuniform. In section II, we find a class of twisted gravitational waves with wave fronts that have negative Gaussian curvature. This class contains Harrison’s solution [8] as a special case, as discussed in section III. In section IV, we study the geodesics of a new member of this class and show that timelike and null geodesics behave as in the Harrison solution [9]. How can one observationally distinguish twisted gravitational waves from plane waves? Linearized gravitational waves in GR are transverse, a property that is present in exact plane waves as well. On the other hand, a twisted gravitational wave exhibits a longitudinal feature in addition to its transverse property. Can this longitudinal character be an observational signature of twisted gravitational waves? We address this problem in section V by comparing the solutions of the Jacobi equation in two related situations involving plane and twisted gravitational waves. Section VI contains a discussion of our results.

II.

A CLASS OF EXACT GRAVITATIONAL WAVES

In our recent work on cosmic jets [9], we considered, among other exact solutions of GR, a Ricci-flat solution due to Harrison [8] that we wrote as ds2 = −ˇ x4/3 dt2 + λ2 uˇ6/5 dx2 + xˇ−2/3 uˇ−2/5 dy 2 + xˇ4/3 dz 2 ,

(5)

4 where xˇ :=

x , T0

uˇ :=

t−z . T0

(6)

We use units such that G = c = 1 throughout this paper; moreover, the signature of the spacetime metric is +2 and greek indices run from 0 to 3, while latin indices run from 1 to 3. Here, T0 is a constant length and λ > 0 is a dimensionless parameter. The dimensionless quantity uˇ is simply proportional to the retarded null coordinate u = t − z. Equation (5) is related by means of straightforward coordinate transformations to the second degenerate solution that Harrison classified as the “III-2; D2 ” metric [8]. This spacetime is of type D in the Petrov classification [10]. The parameter λ can be set equal to unity by rescaling the spacetime coordinates. Harrison’s spacetime has curvature singularities at w = 0, where w := xˇ1/3 uˇ1/5 . Harrison’s solution describes an exact gravitational wave (GW) propagating in the z direction. The wave front is the (x, y) surface with metric dσ 2 = dx2 + xˇ−2/3 dy 2 ,

(7)

which has negative Gaussian curvature. It thus represents an exact unidirectional nonplanar gravitational radiation field; that is, a twisted gravitational wave (TGW). To investigate such TGW solutions, we consider a metric of the form ds2 = −eA(x,u) (dt2 − dz 2 ) + eB(x,u) dx2 + eC(x,u) dy 2 ,

(8)

where we assume for the sake of simplicity that all lengths are measured in units of T0 = 1. The field equations for exact gravitational waves, namely, Rµν = Λ gµν , have been worked out for the case of metric (8) in Appendix A. The Ricci-flat field equations in this case can be obtained from Eqs. (A3)–(A7) given in Appendix A by setting the cosmological constant Λ equal to zero, namely, A2x + 2 Ax Cx = 0 ,

(9)

2 Axx + A2x − Ax (Bx + Cx ) = 0 ,

(10)

2 Cxx + Cx2 − A2x − Bx Cx = 0 ,

(11)

where Ax := ∂A/∂x, etc.,

Axu + Cxu =

1 1 Ax (Bu + Cu ) + Cx (Bu − Cu ) 2 2

(12)

5 and Au (Bu + Cu ) =

1 2 (B + Cu2 ) + Buu + Cuu . 2 u

(13)

In the analysis of these equations, it proves useful to introduce A = eA ,

B = eB ,

C = eC .

For instance, in terms of these new functions, Eq. (13) takes the form 1 Bu2 Cu2 Buu Cuu Au Bu Cu =− + + 2 + + . 2 A B C 2 B C B C

(14)

(15)

It follows from Eq. (9) that either Ax = 0 in case A or Ax + 2 Cx = 0 in case B. These cases will be investigated in turn in the rest of this section.

A.

Ax = 0

This condition implies that A is only a function of u. Since dt2 − dz 2 = dudv, where v = t+ z is the advanced null coordinate, we can absorb A in the redefinition of the retarded null coordinate. Thus one may, in effect, set A = 0. Then, Eqs. (11) and (12) can be written as Bx − Cx = 2

Cxx , Cx

Bu − Cu = 2

Cxu , Cx

(16)

which can be simply integrated and the result is B = C + 2 ln Cx + 2 ln λ ,

(17)

where λ > 0 is an integration constant. The problem thus reduces to finding solutions of the system

C2 B=λ x, C 2

1 2

Bu2 Cu2 + B2 C 2

=

Buu Cuu + . B C

(18)

It is simple to find a solution using separation of variables. That is, let C = X(x) U(u) .

(19)

Then, it follows from Eq. (18) that 2 UUuu − Uu2 = 0, so that U(u) = (au + b)2 ,

(20)

6 where a and b are integration constants and the metric in this case takes the form 2 2 2 1 2 2 2 dX + X dy + dz 2 , ds = −dt + (au + b) λ X

(21)

which represents flat spacetime. We note that λ2

1 dX 2 + X dy 2 = dr 2 + r 2 dθ2 , X

(22)

where X = r 2 /(4λ2 ) and y = 2λ θ. This case is further discussed in Appendix B. Let us now return to system (18) and introduce the change of variables P=

√

B,

Q=

√

C,

(23)

since B and C are positive by definition. Then, system (18) reduces to P 2 = 4λ2 Q2x ,

Puu Quu + = 0. P Q

(24)

The resulting gravitational waves in this case all have metrics of the form ds2 = −dt2 + dz 2 + 4 λ2 Q2x dx2 + Q2 dy 2 ,

(25)

which turn out to represent plane waves, since on the wave front u must be replaced by a constant u0 , say, such that Qx (x, u0 ) and Q(x, u0 ) are then only functions of x; hence, the resulting metric on the wave front, dσ 2 = 4 λ2 dQ2 + Q2 dy 2 ,

(26)

is flat. Moreover, with k = ∂v = ∂t + ∂z , it is simple to see that Tµν = k[µ,ν] = 0, so that the rays are twist-free, as expected. System (24) can be further reduced to P = ±2λQx ,

(Q Quu )x = 0 .

(27)

Exact solutions of Q Quu = ̟(u), i.e. solutions that can be explicitly expressed in terms of familiar functions, are unknown; however, it is possible to solve this equation when ̟(u) = κ, where κ is a constant. In this case, Quu = κ/Q can be integrated once to yield Q2u = 2 κ ln Q + E(x) ,

(28)

7 where E(x) is a certain energy function. Assuming Q is a monotonically increasing function of u, we can write Z

Q(x,u)

Q(x,0)

dZ p = u, 2 κ ln Z + E(x)

(29)

where the integral can be expressed in terms of the error function. These results contain an infinite class of solutions that represent plane gravitational waves propagating along the z direction.

B.

Ax + 2 Cx = 0

It follows from this condition that A + 2 C is only a function of u. Moreover, with Ax = −2 Cx , Eqs. (10) and (11) both reduce to Bx + 3 Cx = 2 Cxx /Cx , while Eq. (12) implies that Bu + 3 Cu = 2 Cxu /Cx . Thus, B + 3 C − 2 ln Cx is simply a constant. Therefore, the field equations reduce in this case to A=

α(u) , C2

B = λ′2

Cx2 C5

(30)

and Eq. (15). Here, α(u) is an integration function and, as before, λ′ > 0 is a constant. The resulting gravitational waves in this case all have metrics of the form ds2 =

2 α(u) 2 2 ′2 Cx (−dt + dz ) + λ dx2 + C dy 2 , C2 C5

(31)

which turn out to represent nonplanar waves with wave fronts that have negative Gaussian curvature. Indeed, following the same line of argument as in Eqs. (25) and (26), the wave front has a metric of the form dσ 2 = λ′2

dC 2 + C dy 2 , C5

(32)

which has negative Gaussian curvature KG = −C 3 (x, u0 )/λ′2 , as discussed in Appendix B. Moreover, the twist tensor in this case has nonzero components T01 = −T10 = T13 = −T31 = α(u) Cx /C 3 .

Let us now return to system (30) and note that in the corresponding general TGW metric (31), it is possible to absorb α(u) in the redefinition of the retarded null coordinate. Thus we can, in effect, set α(u) = 1 in our system. Our general TGW metric turns out to be of type II in the Petrov classification.

8 Defining P and Q as in Eq. (23), we find A = Q−4 ,

P = ±2λ′

Qx , Q4

(33)

so that Eq. (15) takes the form Puu Quu Qu + = −4 P Q Q

Pu Qu + P Q

.

(34)

Using the fact that Qxu Qu Pu = −4 , P Qx Q we can reduce the problem to the following differential equation for Q, Q2 Qxuu − 4 Q Qu Qxu − (3 Q Quu − 8 Q2u ) Qx = 0 . Equation (36) turns out to be equivalent to Q2u Quu −2 4 = 0. Q3 Q x

(35)

(36)

(37)

Hence, we have

Quu Q2u − 2 = −υ(u) , Q3 Q4 where υ(u) is an integration function. Let us note that Eq. (38) can be written as Qu = −υ(u) Q . Q2 u Let us define S such that S :=

then, Eq. (39) takes the form

1 ; Q

Suu = υ(u) S −1 .

(38)

(39)

(40)

(41)

It is remarkable that we encounter here in case (B) the same type of equation that we encountered in case (A). As before, exact solutions of Eq. (41), namely, solutions expressible in terms of familiar functions, are not known, except when υ(u) is a constant k. With υ(u) = k, Eq. (41) can be integrated once to yield Su2 = 2 k ln S + E(x) ,

(42)

where E(x) is a new energy function. Assuming S is a monotonically increasing function of u, we can write Z

S(x,u) S(x,0)

dZ p = u, 2 k ln Z + E(x)

(43)

9 where the integral can be expressed in terms of the error function. We have thus demonstrated the existence of an infinite class of TGWs, namely, exact unidirectional nonplanar gravitational wave solutions of general relativity theory. Finally, as is evident from the tenor of our work, we simply attempt to develop a physical idea in this paper. On the other hand, we should mention that our solutions belong to the Kundt class [3] and are probably all known in some other coordinate systems and within certain classification schemes; in this connection, see Ref. [11] and sections 24.4.1 and 31.5.2 of Ref. [3].

III.

SIMPLE TGW SOLUTIONS

To look for manageable TGW solutions, we solve Eq. (39) for υ = 0. In view of Eq. (41), the general solution is given by Q=

1 , u f (x) + h(x)

(44)

where f (x) and h(x) are integration functions and cannot both be zero functions simultaneously. The spacetime metric can be worked out using Eq. (33) and the result is 2 dh df 2 4 2 2 ′2 4 + dx2 + (u f + h)−2 dy 2 . (45) ds = (u f + h) (−dt + dz ) + 4 λ (u f + h) u dx dx Taking advantage of the freedom in the choice of the x coordinate in Eq. (45), we have a simple TGW solution that depends upon only one arbitrary function of x. It turns out that solution (45) is of type D in the Petrov classification. Of this class of solutions, we are interested in this section in three special cases where f (x) = 0, h(x) = 0 and, finally, f (x) = 1/q and h(x) = x/q, where q is a constant. Let us now examine these cases in turn.

A.

f (x) = 0, Q = 1/h(x)

The spacetime metric in this case takes the form ds2 = −χ4 (dt2 − dz 2 ) + 4 λ′2 χ4 dχ2 + χ−2 dy 2 ,

(46)

where we have replaced h(x) by χ. The result is a spacelike Kasner metric [3]; indeed, with a new spacelike coordinate x, x :=

2 λ′ 3 χ , 3

(47)

10 and constant rescalings of the spacetime coordinates, the metric takes the standard static form [12] ds2 = −x2p1 dt2 + dx2 + x2p2 dy 2 + x2p3 dz 2 ,

(48)

with p1 + p2 + p3 = p21 + p22 + p23 = 1. In our case, p1 = p3 =

2 , 3

1 p2 = − , 3

(49)

so that we can write our solution in the convenient form ds2 = x4/3 (−dt2 + dz 2 ) + dx2 + x−2/3 dy 2 .

(50)

This particular static Kasner spacetime is of type D in the Petrov classification and its three independent Killing vectors are ∂v = ∂t + ∂z , ∂y and ∂z ; moreover, its homothetic vector field is given by t ∂t + 3 x ∂x + 4 y ∂y + z ∂z .

B.

h(x) = 0, Q = 1/[u f (x)]

The spacetime metric in this case can be expressed as ds2 = −u4 χ4 du dv + 4 λ′2 u6 χ4 dχ2 + u−2 χ−2 dy 2 ,

(51)

where we have replaced f (x) by χ. Furthermore, designating u5 /5 as u in metric (51), using Eq. (47) and constant rescalings of the coordinates finally lead to ds2 = x4/3 (−dt2 + dz 2 ) + u6/5 dx2 + x−2/3 u−2/5 dy 2 ,

(52)

which is Harrison’s TGW solution [8, 10]. Harrison’s spacetime has independent Killing vectors ∂v , ∂y and a homothetic vector field 5 t ∂t + 6 x ∂x + 12 y ∂y + 5 z ∂z . The wave front in this case has a metric of the form dσ 2 = dx2 + x−2/3 dy 2 and a Gaussian curvature KG = −4/(9 x2 ). C.

f (x) = 1/q, h(x) = x/q and Q = q/(x + u)

The spacetime metric in this case is given by ds2 =

1 4 λ′2 4 2 2 (x + u) (−dt + dz ) + (x + u)4 dx2 + q 2 (x + u)−2 dy 2 . 2 6 q q

(53)

11 Let us first rescale the y coordinate by replacing q y by y. Next, we define positive constants λ0 and λ via

2 λ′ 1 , λ := . q2 |q|3 In this way, the metric of the new TGW solution takes the final form λ0 :=

ds2 = λ0 (x + u)4 (−dt2 + dz 2 ) + λ2 (x + u)4 dx2 + (x + u)−2 dy 2 .

(54)

(55)

This is of type D in the Petrov classification and represents a nonplanar gravitational wave propagating in the z direction. Its wave front has metric dσ 2 = λ2 (x + u0 )4 dx2 + (x + u0 )−2 dy 2 ,

(56)

where u0 is a constant. The Gaussian curvature of this surface is negative as well and is given by (cf. Appendix B) 4 . (x + u0 )6 This TGW solution is further discussed in the next section. KG = −

λ2

(57)

Of the three cases we have discussed in this section, only the Harrison solution and the new solution represent propagating gravitational waves. The special Kasner solution is static and represents a limiting form of these solutions; that is, if we replace u by a constant u0 6= 0 in the metric functions of the Harrison and new solutions, then the resulting metrics are equivalent to the special Kasner metric via constant rescalings of the spacetime coordinates. Similarly, if we replace x by a constant x0 6= 0 in the metric functions of the Harrison and new solutions, we recover in the same manner the metric of the plane gravitational wave ds2 = −dt2 + u6/5 dx2 + u−2/5 dy 2 + dz 2 ,

(58)

which is a special case of the type N metrics discussed in Section V of Ref. [9] with s2 = 3/5 and s3 = −1/5. Thus our simple TGW solutions may be considered to be nonlinear superpositions of this special linearly polarized plane gravitational wave and the special static Kasner solution (50). A Ricci-flat solution in GR has four algebraically independent scalar polynomial curvature invariants given by [3] I1 = Rµνρσ Rµνρσ − iRµνρσ R∗ µνρσ

(59)

I2 = Rµνρσ Rρσαβ Rαβ µν + iRµνρσ Rρσαβ R∗ αβ µν .

(60)

and

12 The class of spacetimes under consideration here are all algebraically special of type D; hence, I13 = 12 I22. In the case of the general solution (45), we find I1 =

12 , λ′4 (uf + h)12

I2 = −

12 , λ′6 (uf + h)18

(61)

so that the timelike hypersurface u f (x) + h(x) = 0 is the curvature singularity of this spacetime. More specifically, for the new solution (55), we find I1 =

λ4

192 , (x + u)12

I2 = −

λ6

768 , (x + u)18

(62)

so that the timelike hypersurface x + u = 0 is a curvature singularity in this case. For the Harrison solution (52), we have I1 =

64 , 27 x4 u12/5

I2 = −

256 , 243 x6 u18/5

(63)

so that the timelike hypersurface w := x1/3 u1/5 = 0 is the curvature singularity of the Harrison solution. Finally, in the limiting case of the special Kasner spacetime (50), the curvature invariants are given by I1 = 64/(27 x4) and I2 = −256/(243 x6 ); hence, the timelike hypersurface x = 0 is the curvature singularity of the spacelike Kasner metric. The twisted gravitational wave solutions have nonzero twist tensor. To see this explicitly, we note that for the new solution (55), the twist tensor can be expressed as 0 1 (Tµν ) = Φ 0 0

−1 0 0

0 0 −1 , 0 0 0 1 0 0

Φ = 2λ0 (x + u)3 ,

(64)

while for the Harrison solution the only difference is that Φ = 32 x1/3 .

IV.

TIMELIKE AND NULL GEODESICS OF THE NEW TGW SOLUTION (55)

To gain insight into the nature of a twisted gravitational wave spacetime, it is useful to study the motion of free test particles in such a gravitational field. For this purpose, we choose the new TGW solution (55) and define W via W = x+u,

(65)

13 so that the metric of the new TGW can be written as ds2 = −λ0 W 4 (dt2 − dz 2 ) + λ2 W 4 dx2 + W −2 dy 2 .

(66)

This spacetime has one null and two spacelike Killing vector fields given by ∂v = ∂t + ∂z ,

∂x + ∂z ,

∂y ,

(67)

respectively, and a homothetic vector field t ∂t + x ∂x + 4 y ∂y + z ∂z . We note that

√

(68)

−g = λ0 λ W 5 , which vanishes at the curvature singularity W = 0.

The motion of free test particles and rays of radiation in this TGW spacetime involves three constants that can be obtained from the projection of the 4-velocity vector, x˙ µ = dxµ /dη, of the test particle on the null and spacelike Killing vector fields. Here, η is either the proper time along a timelike geodesic or an affine parameter along the path of the null geodesic ray. Thus we have Cv du = t˙ − z˙ = 4 , dη W

λ2 x˙ + λ0 z˙ =

C0 , W4

y˙ = Cy W 2 ,

(69)

where Cv , C0 and Cy are constants of the motion. Let us note here that a special solution of the null geodesic equation corresponds to the Killing vector field ∂v with Cv = 0, Cy = 0 and constant W along the rays that indicate the propagation of the background TGW. For other null geodesics, we require that Cv 6= 0. Likewise, we are interested in future-directed timelike geodesics; hence, we assume that Cv > 0, so that for constant z we have dt/dη > 0. Moreover, gµν x˙ µ x˙ ν = −ǫ, where ǫ = 1 or ǫ = 0, depending upon whether the geodesic is timelike or null, respectively; that is, λ0 W 4 (−t˙2 + z˙ 2 ) + λ2 W 4 x˙ 2 + W −2 y˙ 2 = −ǫ ,

ǫ = 1, 0 .

(70)

The geodesic equations of motion can be simply obtained from a Lagrangian of the form (ds/dη)2. In this way, we find that the geodesic equation for the x coordinate takes the form 2ǫ d 2 4 [λ W x] ˙ = − 3Cy2 W , dη W

(71)

where Eq. (70) has been employed. Using W 4 du/dη = Cv > 0, we can express Eq. (71) as d2 x + 2 ζ W 3 + 3 ζ′ W 5 = 0 . du2

(72)

14 where ζ and ζ ′ are constants given by ζ=

ǫ , λ2 Cv2

ζ′ =

Cy2 . λ2 Cv2

(73)

We recall that W = x + u by definition; therefore, Eq. (72) can be integrated once and the result is

dW du

2

+ ζ W 4 + ζ′ W 6 = E ,

(74)

where E must therefore be positive and is given by E=1+

λ0 C v + 2 C 0 , λ2 C v

(75)

which follows from Eqs. (69)–(74). Equation (74) bears a remarkable resemblance to Eq. (95) of Section VI of Ref. [9], which involved the geodesic equation for the Harrison TGW solution. As in [9], one can interpret Eq. (74) via a one-dimensional motion of a classical particle with positive energy E in a simple positive symmetric effective potential well of the same structure as in the case of Harrison’s TGW solution. The motion is periodic with turning points at ±W0 , where W0 > 0 and ±W0 are the only real solutions of Eq. (74) with dW/du = 0. We recall that W = 0 is the location of the spacetime singularity; therefore, the geodesic motion can start from −W0 or W0 and end up at the curvature singularity. A.

Oblique Cosmic Jet

Plane gravitational wave spacetimes have parallel rays and are of type N in the Petrov classification. This means that the four principal null directions of the Weyl tensor coincide and are all parallel to the direction of propagation of the plane wave and hence perpendicular to the uniform wave front. The timelike geodesics of these spacetimes have the peculiar property that they all asymptotically line up parallel to the direction of motion of the wave and their Lorentz factors approach infinity. This cosmic jet property was first demonstrated in Ref. [13] and further elaborated in Ref. [9]. In the case of Harrison’s nonplanar GW, the nonuniformity of the wave front led to an oblique cosmic jet, namely, the direction of the cosmic jet deviated from the direction of propagation of the wave. On the other hand, Harrison’s TGW solution, just like the new TGW solution under consideration in this section, is of type D in the Petrov classification, which means that the four principal null

15 directions in this case indicate only two directions (each with multiplicity 2): one along the direction of wave propagation and another along some oblique direction. It is possible that oblique cosmic jets occur for all TGW spacetimes. Therefore, it would be interesting to see if the same result holds for the new TGW solution as well. To this end, let us refer the motion of timelike (and null) geodesics to fiducial observers that are at rest in this spacetime. The natural tetrad frame of these static observers is given by −1/2

eˆ0 = λ0

W −2 ∂t ,

eˆ1 = λ−1 W −2 ∂x ,

eˆ2 = W ∂y ,

−1/2

eˆ3 = λ0

W −2 ∂z .

(76)

z˙ . t˙

(77)

˙ x, Projecting x˙ µ = (t, ˙ y, ˙ z) ˙ on eµ αˆ results in x˙ αˆ = Γ(1, Vx , Vy , Vz ), where Γ=

p

λ0 W 2 t˙ ,

λ x˙ Vx = √ , λ0 t˙

Vy = √

Cy , λ0 W t˙

Vz =

It follows from Eqs. (69)–(74) that " # 2 2 ǫ + Cy2 W 2 C dx λ v t˙ = + +1 , 2λ0 Cv 2W 4 λ0 du Cv dx , W 4 du y˙ = Cy W 2

x˙ =

and

" # 2 2 ǫ + Cy2 W 2 dx Cv λ + −1 . z˙ = 4 2λ0 Cv 2W λ0 du

(78) (79) (80)

(81)

Moreover, as W → 0, (dW/du)2 → E. As free test particles approach the spacetime singularity at W = 0, we find that for W → 0, a cosmic jet develops with Γ → ∞ and (Vx , Vy , Vz ) → (sin Θ± , 0, cos Θ± ) ,

(82)

where

Θ± λ √ cot =√ ± E−1 . (83) 2 λ0 In this equation, the upper (lower) sign indicates that the singularity at W = 0 is approached

from the turning point at −W0 (W0 ). The oblique character of the cosmic jet in this case is in complete correspondence with Harrison’s TGW. Are there observable differences between plane and twisted gravitational waves? To investigate this issue, we can compare and contrast the influence of these waves on the propagation of fields and on congruences of massive test particles. The case of a massless scalar field is treated in Appendix C. Tidal effects of TGWs are studied in the next section.

16 V.

JACOBI EQUATION

Imagine the world line of an arbitrary reference observer that is static in the spacetime under consideration. Let τ be the proper time and λµ αˆ (τ ) be a Fermi-Walker transported tetrad along this world line, where at each event x¯µ (τ ) we imagine all spacelike geodesic curves that issue perpendicularly from this event and generate a local hypersurface. We assume that xµ is an event on this hypersurface that can be connected to x¯µ (τ ) via a unique ˆ

spacelike geodesic of proper length ς. We assign Fermi coordinates X µˆ = (T, X i ) to event xµ , where T := τ ,

ˆ

ˆ

X i := ς ξ µ (τ ) λµ i (τ ) .

(84)

Here, ξ µ (τ ) is the unit spacelike vector at x¯µ (τ ) that is tangent to the unique geodesic ˆ

connecting x¯µ (τ ) to xµ , so that ξ µ (τ ) λµ 0 (τ ) = 0. Thus along the reference world line, ˆ

ξ µ (τ ) λµ i (τ ), for i = 1, 2, 3, are the corresponding direction cosines at proper time τ . In the Fermi coordinate system, the reference observer is permanently fixed at the spatial origin (X = 0). We are interested in the equation of motion of a neighboring free test particle relative to the reference observer in the Fermi coordinate system. Neglecting the relative velocity, the reduced geodesic equation can be expressed as ˆ

d2 X i ˆ ˆ + Ai + (Rˆ0ˆiˆ0ˆj + Aˆi Aˆj )X j = 0 , 2 dT

(85)

ˆ

where Ai is the measured 4-acceleration of the fiducial static observer and Rˆ0ˆiˆ0ˆj are the corresponding components of the tidal matrix, namely, ˆ

Ai (T ) =

Dλµ ˆ0 ˆi λµ , dτ

µ ν ρ σ Rαˆ βˆ ˆγ δˆ(T ) := Rµνρσ λ α ˆ λ βˆ λ γ ˆ λ δˆ .

(86)

For background material on equations of motion of free test particles in Fermi coordinates, we refer to our recent paper [14]; further material is contained in [15, 16] and the references cited therein. The Fermi coordinate system is generally admissible in a cylindrical region of radius |X| ∼ R in the spacetime domain around the reference world line, where R is a certain minimal radius of curvature. Let us first consider the simple case of motion in the special linearly polarized plane gravitational wave given by metric (58). In this case, the curvature singularity of spacetime occurs at u = 0. The static observers in this spacetime follow geodesics and have a natural tetrad frame that is parallel transported along their geodesic world lines, namely, λˆ0 = ∂t ,

λˆ1 = u−3/5 ∂x ,

λˆ2 = u1/5 ∂y ,

λˆ3 = ∂z .

(87)

17 Let us choose a static observer with world line (t, x, y, z) = (τ, 0, 0, 0), along which we construct a Fermi normal coordinate system. Using the results of Appendix B of Ref. [9] for the curvature of the plane wave, the equations of motion of nearby free test particles in Fermi coordinates take the form ˆ

6 d2 X 1 ˆ + X1 = 0 , 2 2 dT 25 T

(88)

ˆ

6 d2 X 2 ˆ − X2 = 0 dT 2 25 T 2 ˆ

(89) ˆ

ˆ

ˆ

and d2 X 3 /dT 2 = 0. With the initial conditions that at T = 1, (X 1 , X 2 , X 3 ) = (X0 , Y0 , Z0) ˆ

and dX i /dT = 0 for i = 1, 2, 3, we find ˆ

ˆ

X 1 = X0 (−2 T 3/5 + 3 T 2/5 ) ,

X2 =

1 Y0 (T 6/5 + 6 T −1/5 ) , 7

ˆ

X 3 = Z0 .

(90)

The plane wave in this case has the character of the plus (⊕) polarization and for T : 1 → 0, ˆ

ˆ

the particles move in the (X 1 , X 2 ) plane such that at the curvature singularity T = 0, ˆ

ˆ

ˆ

(X 1 , X 2 , X 3 ) = (0, ∞, Z0). Let us note an important limitation of our result here and in the rest of this section: Though we seek solutions of the Jacobi equation, the results are only valid so long as |X| is sufficiently small in accordance with the admissibility of Fermi normal coordinates. Thus in Eq. (90), for instance, assuming that initially |X0 | is sufficiently small, ˆ

the motion in X 2 is valid for only a very short time interval.

A.

Jacobi Equation for Harrison’s TGW

Next, we consider Harrison’s twisted gravitational wave with metric ds2 = −x4/3 dt2 + u6/5 dx2 + x−2/3 u−2/5 dy 2 + x4/3 dz 2 ,

(91)

where we have set λ = 1 with no loss in generality. This spacetime has curvature singularities at x = 0 and u = 0. The static observers in this spacetime have a natural tetrad frame given by λˆ0 = x−2/3 ∂t ,

λˆ1 = u−3/5 ∂x ,

λˆ2 = x1/3 u1/5 ∂y ,

λˆ3 = x−2/3 ∂z .

(92)

The world line of a static observer in Harrison spacetime is accelerated, namely, A = ∇λˆ0 λˆ0 =

2 ∂x . 3 x u6/5

(93)

18 Imagine that Sµ is a vector that is Fermi–Walker transported along λµ ˆ0 ; then, dSµ + Γµαβ λα ˆ0 Sβ = (A · S) λµ ˆ0 − (λˆ0 · S) Aµ . dτ

(94)

It is straightforward to check that tetrad frame (92) is Fermi–Walker propagated along λˆ0 . Let us now consider the frame components of the 4-acceleration and of the relevant components of the electric part of the Riemann tensor. We find A=

2 λˆ 3 x u3/5 1

(95)

and the symmetric and traceless tidal matrix K, Kˆiˆj = Rˆ0ˆiˆ0ˆj ,

(96)

has nonzero components, 6 2 1 1 + , 2 6/5 4/3 9x u 25 x u2 6 1 1 2 − , = − 2 6/5 4/3 9x u 25 x u2 4 , = 2 9 x u6/5 2 = Kˆ3ˆ1 = . 5/3 5 x u8/5

Kˆ1ˆ1 = − Kˆ2ˆ2 Kˆ3ˆ3 Kˆ1ˆ3

(97)

Let us now consider a static observer located at (x0 , y0, z0 ), where x0 6= 0. Then, for this 2/3

observer, t = t0 + τ /x0

2/3

and u = u0 + τ /x0 , where u0 := t0 − z0 . We establish a Fermi

normal coordinate system in the neighborhood of this observer. The equation of motion of a nearby free test particle in the Fermi system is given by Eq. (85), where Aˆi and Kˆiˆj can be obtained from Eqs. (96) and (97) with x = x0 ,

u = u0 +

T x20 /3

.

(98)

. We assume for the sake of simplicity that x0 = 1 and u0 = 0. Equation (85) can then be written as

ˆ

d2 X 1 2 1 + + dT 2 3 T 3/5

2 1 6 1 2 1 ˆ ˆ + X1 + X3 = 0 , 6/5 2 9T 25 T 5 T 8/5 ˆ 2 1 6 1 d2 X 2 ˆ − + X2 = 0 2 6/5 2 dT 9T 25 T

(99) (100)

19

FIG. 1: Numerical integration of Eqs. (99)–(101). Initial conditions at T = 1 are chosen such that ˆ

ˆ

ˆ

X i = 1 and dX i /dT = 0. Closer inspection reveals that X 1 indeed vanishes at T = 0, cf. Appendix D.

and

ˆ

d2 X 3 2 1 4 1 ˆ ˆ + X1 + X3 = 0 , 2 8/5 6/5 dT 5T 9T

(101)

which can be integrated from T = 1 to the curvature singularity at T = 0 with the boundary ˆ

ˆ

ˆ

ˆ

conditions that at T = 1, (X 1 , X 2 , X 3 ) = (X0 , Y0 , Z0 ) and dX i /dT = 0 for i = 1, 2, 3. See, for example, Figure 1. ˆ

It is the variation of X 3 with Fermi time T that is the distinguishing feature of the TGW. ˆ

To see this, let us first note that in system (99)–(101), the equation for X 2 decouples and can be treated separately. In fact, let 5 ϑ = √ T 2/5 ; 3 2

ˆ

X 2 = T 1/2 Fν (ϑ) ,

(102)

then, Eq. (100) takes the form of the modified Bessel equation of order ν = 7/4, ϑ2

d2 Fν dFν +ϑ − (ϑ2 + ν 2 ) Fν = 0 . 2 dϑ dϑ

(103)

Here, Fν (ϑ) for ν = 7/4 is a constant linear combination of I7/4 (ϑ) and K7/4 (ϑ). Moreover, as T → 0, one can show from the properties of the modified Bessel functions of the first and ˆ

second kind that X 2 diverges as T −1/5 , in agreement with Eq. (90) for the corresponding ˆ

ˆ

plane gravitational wave. Next, for the (X 1 , X 3 ) system, it proves convenient to introduce a new variable S, S = − ln T ,

(104)

20 ˆ

ˆ

so that as T : 1 → 0, we have that S : 0 → ∞. In terms of S, the (X 1 , X 3 ) system can be expressed as ˆ

ˆ

d2 X 1 dX 1 2 −7S/2 + + e + dS 2 dS 3 ˆ

2 −4S/5 6 e + 9 25

ˆ

X1 +

2 −2S/5 ˆ3 e X = 0, 5

(105)

ˆ

d2 X 3 dX 3 2 −2S/5 ˆ1 4 −4S/5 ˆ3 + + e X + e X = 0. dS 2 dS 5 9

(106)

As proved in Appendix D, all solutions of this system and their derivatives are bounded ˆ

for 0 ≤ S < ∞ and X 1 converges to zero as S → ∞. Numerical experiments reported in ˆ

Appendix D suggest that X 3 also has a limit (which may not be zero) as S → ∞. Ignoring terms involving e−2S/5 , e−4S/5 and e−7S/2 , system (105)–(106) reduces to the corresponding

system for the plane GW (58). In fact, these terms vanish for S → ∞. Thus the main ˆ

difference between the TGW and the plane wave (58) is that in time X 3 deviates from its value Z0 at S = 0 in the case of TGW. The influence of the exact plane wave (58) on test particles that are initially at rest in the Fermi system has the same transverse character as for linearized gravitational waves moving along the z direction; however, the corresponding ˆ

TGW has in addition a longitudinal influence as well resulting in the X 3 component of the ˆ

motion. This longitudinal feature of TGW and the associated temporal variation in X 3 may provide an observational signature for TGWs.

B.

Jacobi Equation for TGW (55)

To confirm the longitudinal component of motion of free test particles induced by a TGW, it is interesting to study the Jacobi equation for the TGW solution (55). In this case, the static observers’ adapted frame is given by Eq. (76). It can be shown that this tetrad frame is Fermi-Walker transported along e0 . In fact, e0 has a nonzero 4-acceleration, 2 1 1 A = ∇e0 e0 = 5 ∂x − ∂z , W λ2 λ0

(107)

where, as before, W = x + u. The projection of the 4-acceleration of static observers in this spacetime on their tetrad frame field can be expressed as ! 1 1 2 e1 − 1/2 e3 . A= 3 W λ λ0

(108)

21 Furthermore, tidal matrix (96) in this case has nonzero components given by 2(−λ0 + 3λ2 ) −6 W , λ0 λ2 2(λ0 + 3λ2 ) −6 W , = − λ0 λ2 4 = 2 W −6 , λ 6 = Kˆ3ˆ1 = 1/2 W −6 . λ0 λ

Kˆ1ˆ1 = Kˆ2ˆ2 Kˆ3ˆ3 Kˆ1ˆ3

(109)

Let us recall the fact that W = 0 corresponds to the curvature singularity in this TGW spacetime. We wish to establish a Fermi normal coordinate system in the neighborhood of a particular static observer; to this end, let us choose the observer located at xµ = (t, x0 , y0 , z0 ), where −1/2

x0 , y0 and z0 are constants and dt/dτ = λ0

W −2 in conformity with the expression for eˆ0

in Eq. (76). Let us define W to be the magnitude of W = x + u along the world line of our reference static observer, namely, W = x0 + t − z0 ,

(110)

1 dW =√ . dτ λ0 W 2

(111)

such that

Integrating this equation we find

p λ0 W 3 = 3 (τ − τ0 ) ,

(112)

where τ0 is an integration constant. To simplify matters, we choose τ0 = 0, so that

√

λ0 W 3 =

3 τ and the curvature singularity now occurs at τ = 0 along the world line of our fiducial static observer. Next, we need to evaluate the 4-acceleration and the tidal matrix along the world line of the fiducial static observer by replacing W with W in Eqs. (108) and (109). The reduced geodesic Eq. (85) then takes the form ˆ

ˆ

ˆ

1 2 2 X3 d2 X 1 2 ℓ 2 X 3 + ℓ + + + ℓ = 0, dT 2 3T 9 T2 9 T2 ˆ

(113)

ˆ

2 d2 X 2 2 2 X 3 + ℓ − =0 dT 2 9 T2

(114)

22

FIG. 2: Numerical integration of Eqs. (113)–(115) with ℓ = 1. Initial conditions are chosen so that ˆ

ˆ

ˆ

X i = 0, dX i /dT = 0 at T = 1. With these conditions X 2 (T ) = 0 identically.

and

ˆ

ˆ

ˆ

3 d2 X 3 2 1 2 X1 4 2 X − + ℓ + (1 + ℓ ) = 0, dT 2 3T 9 T2 9 T2

where ℓ :=

(115)

√

λ0 . λ

(116)

ˆ

As before, the equation for X 2 decouples from the rest and its solution is given by ˆ

X 2 = T 1/2 (ˇb+ T aˇ + ˇb− T −ˇa ) ,

(117)

where ˇb± are integration constants and a ˇ=

1√ 33 + 8 ℓ2 . 6

(118) ˆ

ˆ

Let us introduce the independent variable S = − ln T into (X 1 , X 3 ) system as before. The resulting system can be integrated using elementary methods and we find ˆ

X 1 = [A1 cos(ω1 S + ϕ1 ) + A2 cos(ω2 S + ϕ2 )] e−S/2 −

3 ℓ −S e +2

ℓ2

(119)

and 3 1 e−S , X = − A1 cos(ω1 S + ϕ1 ) + ℓ A2 cos(ω2 S + ϕ2 ) e−S/2 + 2 ℓ ℓ +2 ˆ 3

where ω1 =

1√ 7 + 8 ℓ2 , 6

ω2 =

1√ 15 + 16 ℓ2 6

(120)

(121)

23 and (A1 , ϕ1 , A2 , ϕ2 ) are constants that can be determined from the initial conditions that ˆ

ˆ

ˆ

ˆ

we need to impose on (X 1 , dX 1/dS, X 3 , dX 3 /dS) at, say, S = 0. ˆ

ˆ

It is clear from Eqs. (119)–(120) that X 1 and X 3 are bounded for 0 ≤ S < ∞ and vanish as S → ∞. For the specific TGW under consideration here, the important quantity ˆ

is the temporal variation of X 3 given by Eq. (120). This is the longitudinal signature of the ˆ

ˆ

TGW that exists in addition to the transverse components characterized by the (X 1 , X 2) components. ˆ

ˆ

Figure 2 depicts X 1 and X 3 given in Eqs. (119) and (120) for the special case of ℓ = 1 ˆ

ˆ

and initial data that at T = 1, X i = 0 and dX i /dT = 0 for i = 1, 3.

VI.

DISCUSSION

Linearized gravitational radiation can be represented as a linear superposition of monochromatic plane waves that are transverse and exhibit plus (⊕) and cross (⊗) linear polarization states. In the nonlinear regime, exact solutions of GR that represent plane gravitational waves generally correspond to linearized plane waves propagating along a fixed direction in space. On the other hand, exact gravitational waves in GR propagating along a fixed spatial direction can be nonplanar as well. This paper is about such twisted gravitational waves (TGWs). We have shown that a class of TGWs exists with wave fronts that have negative Gaussian curvature. We have investigated in some detail the properties of two such radiation fields. In particular, our study of the Jacobi equation in these spacetimes reveals a longitudinal signature that is a distinct departure from transversality. In this first discussion of TGWs, many questions remain unanswered. For instance, how can TGWs be generated by realistic sources? Are there TGWs with wave fronts that have positive Gaussian curvature? Are there TGWs with a cosmological constant Λ? Further investigation is necessary to tackle these problems.

Appendix A: Field Equations with Λ

The Einstein field equations in vacuum but with a cosmological constant Λ can be expressed as Rµν = Λ gµν .

(A1)

24 We assume a metric of the form ds2 = −eA(x,u) (dt2 − dz 2 ) + eB(x,u) dx2 + eC(x,u) dy 2 ,

(A2)

where u = t − z is the retarded null coordinate. With metric (A2), field equations (A1) can be reduced to the following five equations: A2x + 2 Ax Cx + 4 Λ eB = 0 ,

(A3)

2 Axx + A2x − Ax (Bx + Cx ) = 0 ,

(A4)

2 Cxx + Cx2 − A2x − Bx Cx = 0 ,

(A5)

where Ax := ∂A/∂x, etc.,

1 1 Ax (Bu + Cu ) + Cx (Bu − Cu ) , 2 2 1 Au (Bu + Cu ) = (Bu2 + Cu2 ) + Buu + Cuu . 2

Axu + Cxu =

(A6) (A7)

It has not been possible to find an exact solution—that is, a solution that can be expressed using familiar functions—of these equations for Λ 6= 0. Appendix B: Gaussian Curvature of the Wave Front

Suppose that the metric of the wave front is given by dσ 2 = eP (x,y) dx2 + eQ(x,y) dy 2 .

(B1)

Then, the Gaussian curvature of this surface is given by 1 1 KG = − e−P [2 Qxx − (Px − Qx )Qx ] − e−Q [2 Pyy + (Py − Qy )Py ] . 4 4

(B2)

As an application of this formula, consider the spacetime metric ds2 = −dt2 + dz 2 + u2 (eP dx2 + eQ dy 2) ,

(B3)

which has the null Killing vector field ∂v = ∂t + ∂z and the homothetic vector field t ∂t + z ∂z . It represents flat spacetime if the following field equation is satisfied eP [2 Pyy + (Py − Qy )Py ] + eQ [2 Qxx − (Px − Qx )Qx ] = 0 .

(B4)

25 That is, Eq. (B4) implies that the Gaussian curvature of the wave front vanishes by Eq. (B2); moreover, the spacetime given by ds2 = −dt2 + dz 2 + u2 (dx2 + dy 2)

(B5)

is simply flat, cf. the end of Appendix B of Ref. [9].

Appendix C: Scalar field on Harrison’s TGW Background

To distinguish TGWs from plane gravitational waves, one can explore the behavior of various perturbing fields on these backgrounds. As an example of this approach, let us consider a massless scalar field on the background of Harrison’s solution (5). We write the corresponding metric in the form ds2 = −x4/3 dudv + λ2 u6/5 dx2 + x−2/3 u−2/5 dy 2,

(C1)

where we assume that x ≥ 0 and u 6= 0. The massless scalar field equation is ✷Ψ = g µν ∇µ (∂ν Ψ) = 0 ,

(C2)

which can be written as ∂ ∂xµ

√

With metric of the form (C1), we have

−g g

√

µν

∂Ψ ∂xν

= 0.

(C3)

−g = λ x u2/5 /2 and the nonzero components of

g µν are given by g uv = g vu = −2 x−4/3 , g xx = u−6/5 /λ2 and g yy = x2/3 u2/5 . Writing Ψ = ei Kv v+i Ky y ψ(x, u) ,

(C4)

the massless scalar wave equation reduces to (x ψx )x − 2iλ2 Kv x−1/3 u4/5 [u2/5 ψu + (u2/5 ψ)u ] − λ2 Ky2 x5/3 u8/5 ψ = 0 .

(C5)

Here, we work with a complex amplitude Ψ. The wave equation is linear; therefore, the real and imaginary parts of Ψ constitute real solutions of the massless scalar wave equation. It may not be possible to express the solution of this equation in terms of familiar functions. To illustrate this point, let us assume that Ky = 0 in Eq. (C5). We then look for a solution of the form ψ(x, u) = X(x) U(u) .

(C6)

26 It is then a simple matter to find U, U(u) = U0 u

−1/5

5iC −1/5 , u exp − 2 2λ Kv

(C7)

where Kv 6= 0 and U0 as well as C are in general complex constants. However, the corresponding differential equation for X can be written as dX 1/3 d x +2CX = 0, x dx dx

(C8)

which does not appear to have a solution expressible in terms of familiar functions except for C = 0. It is a peculiar feature of Eq. (C5) that solutions exist for Kv = 0 as well. In this case, Eq. (C5) reduces to (x ψx )x − λ2 Ky2 x5/3 u8/5 ψ = 0 .

(C9)

3 ξ := λ Ky x4/3 u4/5 . 4

(C10)

Let us now define ξ,

Then, ψ = ϕ(ξ) satisfies the modified Bessel equation of zero order, namely, 1 ϕξξ + ϕξ − ϕ = 0 . ξ

(C11)

It is interesting to compare this situation with the propagation of the massless scalar field on the background of a plane gravitational wave with metric of the form ds2 = −du dv + λ2 u6/5 dx2 + u−2/5 dy 2 ,

(C12)

where the metric coefficients can be obtained from those of Eq. (C1) with x = 1. This is a special case of the metrics discussed in Section V of Ref. [9] with s2 = 3/5 and s3 = −1/5. In this case, the solution of the massless scalar field equation can be written as Ψ = ei Kv v+i Kx x+i Ky y φ(u) ,

(C13)

where φ(u) is given by φ(u) = φ0 u

−1/5

5iKx2 −1/5 5iKy2 7/5 . u + u exp − 4λKv 28Kv

Here, Kv 6= 0 and φ0 is a constant of integration.

(C14)

27 Appendix D: Behavior of solutions of Eqs. (105)–(106)

The coefficients of the linear system (105)–(106) are not all constants; therefore, the analysis of the system is expected to rely on qualitative or numerical methods. These are employed in this appendix. To this end, we transform Eqs. (105)–(106) in the usual manner ˆ

ˆ

ˆ

ˆ

into four first-order linear differential equations involving X 1 , dX 1 /dS, X 3 and dX 3 /dS. Let us first discuss the long-time asymptotics, that is the fate of solutions starting at S = 0 as S → ∞. Based on numerical experiments, it appears that all solutions are bounded and have limits as S → ∞. Moreover, the limiting behavior may be specified as an affine transformation from initial to final data such that ! ˆ ˆ 1 3 dX dX ˆ (S), (S) = (0, 0, 0) , lim X 1 (S), S→∞ dS dS while ˆ

ˆ

ˆ

lim X 3 (S) = −k1 X 1 (0) − k2

S→∞

(D1)

ˆ

dX 1 dX 3 ˆ (0) + k3 X 3 (0) + k4 (0) + k5 , dS dS

(D2)

where the ki are positive constants with approximate values k1 = 0.572 ,

k2 = 0.441 ,

k3 = 0.724 ,

k4 = 0.842 ,

k5 = 0.170 .

(D3)

One way to verify this result is to recall the basic fact from linear systems theory that the general solution is given by a particular solution plus a linear combination of fundamental solutions. A convenient choice is the particular solution starting from zero initial data and the four fundamental solutions corresponding to initial data given by the usual basis vectors of 4-dimensional Cartesian space. All that remains is numerical integration to obtain the limiting values of the particular and fundamental solutions. We find that all solutions ˆ

approach the X 3 direction with their final positions determined by the initial data. In ˆ

contrast to the treatment of the Jacobi equation for the plane GW (58), the initial X 3 component is never fixed during the evolution. However, its limiting value can be the same as its initial value by choosing the remaining components of the initial data so that the third ˆ

component of the limit vector has the value X 3 (0). This happens on a hypersurface in the space of initial data; thus, the probability of this event is zero. The boundedness of solutions is of course implied by the boundedness of a particular solution and all the fundamental solutions. So, the numerical computations suggest that this is indeed the case. For an analytic proof, one may use the variation of parameters

28 formula with respect to the constant part of the system matrix and then employ a Gronwall estimate [17]. More specifically and abstractly, consider a vector ordinary differential equation, ˜ x + ˜b(t) , x˙ = A˜ x + B(t)

(D4)

˜ is a continuous matrix function of time and ˜b(t) is where A˜ is a constant square matrix, B a continuous vector function all defined for t ≥ 0. Using the usual Euclidean norm, suppose ˜ µ and L such that for all t ≥ 0 and all vectors v˜ that there are positive constants M, K, λ, ˜

˜

˜ v| ≤ Ke−λt |˜ |B(t)˜ v|,

|etA v˜| ≤ M|˜ v |,

|˜b(t)| ≤ Le−µt .

(D5)

By variation of parameters, ˜ tA

|x(t)| ≤ |e x(0)| +

Z

t

(t−s)A˜

|e

0

˜ B(s)x(s)| ds +

Z

t 0

˜

|e(t−s)A˜b(s)| ds .

(D6)

Using the inequalities (D5), |x(t)| ≤ M|x(0)| +

Z

0

t

˜

MKe−λs |x(s)| ds +

ML ˜ (1 − e−λt ) . µ

(D7)

The last term is bounded by ML/µ. Gronwall’s inequality implies |x(t)| ≤ (M +

˜ ˜ µ + L M K(1−e−λt µ + L M K/λ˜ ML M K R t e−λs ˜ ds )/λ 0 ≤M )e e ≤M e < ∞, µ µ µ

(D8)

as required. The hypotheses are true for system (105)–(106) recast in the usual manner as a first˜ B ˜ and ˜b. In particular, the order linear system and with the natural definitions of A, four-dimensional constant-coefficient system matrix A˜ has four distinct eigenvalues: three ˜ therefore propagates fundamental negatives and one zero. The fundamental matrix exp(tA) ˜ three of which decay in norm solutions starting in the corresponding eigenspaces of A, ˜ v˜| ≤ M|˜ and one remains constant. Using linearity, the desired estimate |exp(tA) v | follows ˜ and ˜b do not involve for an appropriate choice of constant M. Required estimates for B exponentiation and are straightforward. Similar arguments can be used to prove the conjecture—based on numerical ˆ

experiments—that X 1 (T ) is not only bounded on the interval 0 < T ≤ 1, but in fact ˆ

ˆ

X 1 (T ) → 0 as T → 0+ ; or equivalently, X 1 (S) → 0 as S → ∞. Our proof to follow ˆ

ˆ

also shows that the derivatives dX 1 /dS and dX 3 /dS converge to zero exponentially fast

29 as S → ∞. Unfortunately, our argument proves only an upper bound for their rate of ˆ

convergence to zero. A lower bound would be needed to conclude that both dX 1 /dT and ˆ

dX 3 /dT blow up in absolute value as T → 0+ , a conjecture that is also supported by numerical experiments. To see the problem, recall for example that by the change of variables ˆ

ˆ

dX 1 /dT = − exp(S) dX 1 /dS. This indeterminate form, which is expected to blow up, would ˆ

do so if |dX 1/dS| does not approach zero too rapidly.

Our proof requires several steps that are merely outlined here for the sake of brevity. The underlying idea is to split off the part of the constant-coefficient system matrix A˜ corresponding to its zero eigenvalue and thus to take advantage of the exponential decay afforded by the remaining three eigenvalues: −1, −3/5 and −2/5. The boundedness of solutions already proved is also a key ingredient. ˜ to To accomplish the desired splitting, note that A˜ is diagonalizable. In fact, taking Q be the matrix whose columns are eigenvectors corresponding to the listed eigenvalues and ˜ −1 A˜Q ˜ is diagonal with the eigenvalues the remaining zero eigenvalue in that order, D := Q ˜ y and in the specified order along the main diagonal. With the change of variables x = Q using t instead of S, the first-order linear system is recast in the form ˜ −1 B(t) ˜ Qy ˜ +Q ˜ −1˜b(t) . y˙ = Dy + Q

(D9)

The first three differential equations in this vector system may be recast in the vector form ˜1 (t)z + y4 B ˜2 (t) + ˜b1 (t) , z˙ = D1 z + B

(D10)

˜2 (t) for the where the coupling to the fourth differential equation is through the term y4 B ˜2 (t). This latter matrix is simply the first three components of the last 3 × 1 matrix B ˜ −1 B(t) ˜ Q ˜ and of course B ˜1 (t) is its upper 3 × 3 diagonal block. As in the proof column of Q of boundedness, there are positive constants L and K and exponential estimates given by v| |etD1 v˜| ≤ e−2t/5 |˜

(D11)

and ˜1 (t)˜ |B v | ≤ Le−2t/5 |˜ v| ,

˜2 (t)˜ |B v | ≤ Le−2t/5 |˜ v| ,

|˜b1 (t)| ≤ Ke−7t/5 .

(D12)

Again the variation of parameters formula is employed, the triangle law estimate is made, the boundedness of y4 is noted as a corollary of our previous result, and all exponential

30 estimates are used. After some manipulation and obvious estimates, we find that Z t 2t/5 e |z(t)| ≤ |z(0)| + Lt + 1 + e−2s/5 (e2s/5 |z(s)|) ds .

(D13)

0

Gronwall’s inequality followed by a simple integral estimate implies |z(t)| ≤ (|z(0)| + Lt + 1)e5/2 e−2t/5 .

(D14)

In particular, z(t) converges to zero as t → ∞. This is enough, after some interpretation, ˆ

ˆ

to obtain the claimed limit X 1 (S) → 0 as S → ∞, or, equivalently, X 1 (T ) → 0 as T → 0+ . Acknowledgments

We are grateful to Marcello Ortaggio and Kjell Rosquist for their insightful comments on our manuscript. D.B. thanks the Italian INFN (Naples) for partial support.

[1] B. P. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration), “GW170817: Observation of gravitational waves from a binary neutron star inspiral”, Phys. Rev. Lett. 119, 161101 (2017). [2] H. W. Brinkmann, “Einstein spaces which are mapped conformally on each other”, Math. Ann. 18, 119 (1925). [3] H. Stephani, D. Kramer, M. MacCallum, C. Hoenselaers and E. Herlt, Exact Solutions of Einstein’s Field Equations (Cambridge University Press, Cambridge, UK, 2003), 2nd ed. [4] H. Bondi, F. A. E. Pirani and I. Robinson, “Gravitational waves in general relativity III. Exact plane waves”, Proc. R. Soc. A 251, 519 (1959). [5] J. B. Griffiths, Colliding Plane Waves in General Relativity (Oxford University Press, Oxford, UK, 1991). [6] J. B. Griffiths and J. Podolsk´ y, Exact Space-Times in Einstein’s General Relativity (Cambridge University Press, Cambridge, UK, 2009). [7] K. A. Khan and R. Penrose, “Scattering of two impulsive gravitational plane waves”, Nature 229, 185 (1971). [8] B. K. Harrison, “Exact three-variable solutions of the field equations of general relativity”, Phys. Rev. 116, 1285 (1959).

31 [9] D. Bini, C. Chicone and B. Mashhoon, “Anisotropic gravitational collapse and cosmic jets”, Phys. Rev. D 96, 084034 (2017) [arXiv:1708.01040 [gr-qc]] [10] R. A. d’Inverno and R. A. Russell-Clark, “Classification of the Harrison metrics”, J. Math. Phys. 12, 1258 (1971). [11] W. Kinnersley, “Type D vacuum metrics”, J. Math. Phys. 10, 1195 (1969). [12] L. D. Landau and E. M. Lifshitz, Classical Theory of Fields (Pergamon Press, Oxford, UK, 1971), p. 316. [13] D. Bini and B. Mashhoon, “Peculiar velocities in dynamic spacetimes”, Phys. Rev. D 90, 024030 (2014) [arXiv:1405.4430 [gr-qc]]. [14] D. Bini, C. Chicone and B. Mashhoon, “Relativistic Tidal Acceleration of Astrophysical Jets”, Phys. Rev. D 95, 104029 (2017) [arXiv:1703.10843 [gr-qc]]. [15] B. Mashhoon, “Tidal Radiation”, Astrophys. J. 216, 591-609 (1977). [16] C. Chicone and B. Mashhoon, “Explicit Fermi coordinates and tidal dynamics in de Sitter and G¨ odel spacetimes”, Phys. Rev. D 74, 064019 (2006). [17] C. Chicone, Ordinary Differential Equations with Applications (Springer-Verlag, New York, 2006), 2nd edn.