(CourseWikia - Com) Interacting Gravitational
(CourseWikia - Com) Interacting Gravitational
For photocopying of material in this volume, please pay a copying fee through the Copyright Clearance
Center, Inc., 222 Rosewood Drive, Danvers, MA 01923, USA. In this case permission to photocopy
is not required from the publisher.
Printed in Singapore
v
b2530 International Strategic Relations and China’s National Security: World at the Crossroads
Preface
vii
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-fm page viii
Preface ix
that is, if such a defined angle changes in one set of coordinates, it will
change in any of others, or vice versa. Once this is made clear, we turn
to the gravitational analog of the electromagnetic Faraday rotation (1846a,
1846b) for colliding gravitational plane waves. Using the Bianchi identities,
we show explicitly that the changes of the polarization of a plane gravita-
tional wave is exactly due to the nonlinear interaction with the oppositely
moving gravitational and/or matter waves (Wang, 1991b). Without such
nonlinear interactions, the polarization of each of the gravitational plane
waves remains constant. This has been further shown by studying exactly
solutions for the spacetimes of colliding plane waves.
Among such solutions are those that describe both colliding purely grav-
itational plane waves and colliding gravitational plane waves coupled with
matter fields. The collision and interaction of such waves are interesting
in several aspects. One of them is that the nonlinearity of general relativ-
ity shows up explicitly in the failure of the principle of superposition. As
a result, colliding plane gravitational waves generically develop spacetime
singularities in the future of the collision. Killing–Cauchy horizons can be
also developed in the interacting region of the two colliding plane waves in
some particular cases (Chandrasekhar and Xanthopoulos, 1986a), but they
are not stable (Yutsever, 1987; Griffiths, 2005), and with small but generic
perturbations they ought to be turned into spacetime singularities (Clarke
and Hayward, 1989).
One may argue that such formations of spacetime singularities are due
to the high symmetry of the spacetime. In particular, with plane symmetry
the incoming waves always have infinitely large amount of energy. There-
fore, in more realistic cases the singularities might be replaced by high
curved regions. However, recently numerical studies of two colliding plane-
fronted massless particle waves in asymptotically flat spacetimes showed
that spacetime singularities are still formed, even the total energy of the
incoming waves is finite, but now they are hidden inside horizons, that is,
now black holes are formed (Pretorius and East, 2018).
The structure of this book is as follows. The introduction for the fun-
damental concepts and material are arranged so that they are contained
in two chapters, Chapters 1 and 2. Specifically, Chapter 1 includes the
introduction of some basic physical quantities, the Newman–Penrose for-
malism, the physical interpretation of the optical scalars, and the Einstein
field equations in terms of distribution-valued tensors. On the other hand,
in Chapter 2, the definition of a single plane gravitational wave is given,
and its polarization angle and amplitude are given explicitly in terms of the
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-fm page x
Weyl scalars. The spacetime singularities are also studied, and found that
spacetimes of gravitational plane waves are generally always singular when
written in the Baldwin, Jeffery and Rosen coordinates (Baldwin and Jef-
fery, 1926; Rosen, 1937), and are free of singularities only in some particular
cases (Wang et al., 2018). The relevance of such singularities to memory
effects (see, for example, Favata, 2010; Bieri, Garfinkle and Yunes, 2017)
and soft theorems of gravitons (Hawking, Perry and Strominger, 2016, 2017;
Strominger, 2017) is also mentioned.
In Chapter 3, after the Weyl and Ricci scalars, as well as the Bianchi
identities, are written in terms of distributions for the spacetimes of collid-
ing gravitational plane waves, the polarization of a gravitational plane wave
in the interaction region is defined with respect to a parallelly-transported
frame along its wave path, whereby a gravitational analog of the electro-
magnetic Faraday rotation is investigated for various types of collisions. The
nature of the singularities formed in the interaction region of the two collid-
ing plane gravitational wave due to their mutual focus, and some methods
for generating exact solutions are also discussed.
Chapters 4 is devoted to the studies of the collision and subsequent inter-
action of two pure gravitational plane waves, while Chapter 5 is devoted
to the studies of a gravitational plane wave colliding with a matter wave.
The latter can be an impulsive and/or a shock dust shell, an electromag-
netic wave, or a neutrino wave. In these chapters, we present three classes
of exact solutions of the Einstein field equations, which include most of
the known diagonal and non-diagonal solutions found so far for the colli-
sion of gravitational plane waves without or with the presence of matter
fields. The main properties of these solutions and the effects of polarization
of the colliding gravitational plane waves on the formation and nature of
singularities are investigated in detail.
In Chapter 6, we study isometries between the internal spacetimes of
black holes and the interaction regions of two colliding plane waves, and
show that the interiors of all the known black holes have a one-to-one corre-
spondence to the interaction regions of colliding plane waves. These include
the Schwarzschild, Reissner–Nordström, Kerr–Neman–NUT solutions, and
the ones with a cosmological constant. Finally, in Chapter 7, we present
our concluding remarks, and point several directions one can pursue in the
future, including some observations.
This book is based in a large part on the author’s Ph.D. dissertation
(Wang, 1991f), which was written about three decades ago. But, the mate-
rials and topics are still relevant to current studies of gravitational waves,
especially the phenomena due to the nonlinear effects of the Einstein field
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-fm page xi
Preface xi
Contents
Preface vii
xiii
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-fm page xiv
Contents xv
Bibliography 175
Index 197
b2530 International Strategic Relations and China’s National Security: World at the Crossroads
List of Figures
3.1 The projection of a spacetime for colliding plane waves onto the
(u, v)-plane. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
4.1 The soliton structure for the solutions given by Eq. (4.25). . . . 96
5.1 The projection of the colliding null dust spacetime onto the
(u, v)-plane. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 150
xvii
b2530 International Strategic Relations and China’s National Security: World at the Crossroads
Chapter 1
1
March 7, 2020 9:53 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch01 page 2
etc., with N denoting the dimension of the spacetime.1 Then, the Einstein
field equations take the form,
1
Rab − Rgab − Λgab = κTab ,
2
where Λ is the cosmological constant, Tab the energy–momentum tensor,
and κ ≡ 8πG/c4 , with G denoting the Newtonian constant. In addition,
we also define
1 1
X(ab) ≡ (Xab + Xba ), X[ab] ≡ (Xab − Xba ).
2 2
In a local coordinate basis {∂λ }, Eq. (1.6) can be written in the following
form:
T μν ;ν = 0, (1.21)
which are the equations for the conservation of energy and momentum
(stress) of the source.
March 7, 2020 9:53 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch01 page 6
d2 xμ dxν dxδ
2
+ Γμνδ = 0. (1.23)
dλ dλ dλ
When the equation for a geodesic is reduced to the form of Eq. (1.23), we
say that it is affinely parameterized. It should be noted that such defined
affine parameter λ is not unique, and still subjected to the rescaling and
shift of origin: λ̃ = αλ + λ0 , where α and λ0 are constant. Equation (1.23)
also describes the motion of a free particle.
∂ 2 xμ μ ∂x ∂x
ν δ
= −Γ , xμ = xμ (λ, w). (1.25)
∂λ2 νδ
∂λ ∂λ
We might, in general, identify λ, with the arc length on each of the geodesics.
We prefer, however, to leave λ to be defined just by Eq. (1.25) so that our
following discussions remain also valid for null geodesics.
The family of geodesics gives rise to the vector fields
∂xμ (λ, w)
tμ (λ, w) = , (1.26a)
∂λ
∂xμ (λ, w)
η μ (λ, w) = , (1.26b)
∂w
where tμ (λ, w) is the tangent vector along each geodesic, and η μ (λ, w) is
the vector that describes the deviation of two points on two infinitesimally
near geodesics, which have the same parameter value λ. The vector η μ is
usually called the geodesic deviation vector.
From Eq. (1.26b) we find that the covariant differentiation of η μ along
each geodesic is given by
(a)
where λν ≡ η (a)(σ) λμ(σ) gμν , η (μ)(ν) ≡ η μν . The components of η a represent
the spatial coordinates of a particle that moves nearby the observer, who
moves along the geodesic C [sec Fig. 1.1].
(a)
Contracting Eq. (1.29) with λμ and using Eq. (1.32), we find that the
acceleration of the particle relative to the observer is given by
d2 η (a)
= −K (a)(b) η(b) , (1.34)
dτ 2
where
K aa = 0. (1.38)
Rμν = 0, (1.39)
where
1
Eμνλρ ≡ (gμλ Sνρ + gνρ Sμλ − gνλ Sμρ − gμρ Sνλ ),
2
1
Gμνλρ ≡ (gνρ gμλ − gνλ gμρ ) R, (1.41)
12
1
Sμν ≡ Rμν − gμν R.
4
March 7, 2020 9:53 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch01 page 11
In Eq. (1.41), the tensor Sμν denotes the traceless part of the Ricci tensor.
The Weyl tensor Cμνλρ is thought of as representing the free gravitational
field (Szekeres, 1965), and has all the symmetries of the Riemann tensor
[see Eqs. (1.14a)–(1.14c)]. In addition, it is also traceless
Cλ μλν = 0. (1.42)
Combining the fact that the Weyl tensor has all the symmetries of the Rie-
mann tensor and Eq. (1.42), we can see that the Weyl tensor has only 10
independent components. At any given point of the spacetime, these com-
ponents are completely independent of the Ricci tensor. Globally, however,
the Weyl and Ricci tensors are not independent, as they are connected by
the Bianchi identities [see Eq. (1.15)]. These identities can be now written
in the form (Kundt and Trumper, 1962; Szekeres, 1966)
1
Cμνσρ; ρ = Rσ[μ;ν] − gσ[μ R;ν] , (1.43)
6
where the square bracket denotes the antisymmetrization,
1
A[μν] ≡ (Aμν − Aνμ ). (1.44)
2
The remarkable analogy between the Bianchi identities of Eq. (1.43) and
the Maxwell equations,
F μν;ν = j μ , (1.45)
suggests that the Bianchi identities represent the interaction between the
free gravitational field and matter. If we define the tensor Jμνσ as
1
Jμνσ ≡ Rσ[μ;ν] − gσ[μ R,ν] , (1.46)
6
we have
Jμνλ; λ = 0, (1.47)
J λ ;λ = 0. (1.48)
the parts of the source that do not contribute to Jμνλ are called gravita-
tionally inert. The propagation of the free gravitational field is independent
of the inert parts of the source.
An equivalent form for the decompositions of Eqs. (1.40) and (1.41) is
given by
1
Rμνλρ = Cμνλρ + [gμλ Rνρ + gνρ Rμλ − gνλ Rμρ − gμρ Rνλ ]
2
1
+ [gμρ gνλ − gμλ gνρ ] R. (1.49)
6
When the Weyl tensor Cμνλδ vanishes, the spacetime is said to be confor-
mally flat.
where lμ and nμ are real, and mμ and mμ are complex conjugates. The
orthogonality properties of these vectors are
lμ lμ = nμ nμ = mμ mμ = mμ mμ = 0,
lμ mμ = lμ mμ = nμ mμ = nμ mμ = 0, (1.51)
l nμ = −m mμ = 1.
μ μ
The tetrad indices can be raised and lowered by η (α)(β) and η(α)(β) , respec-
tively, which are given by
⎡ ⎤
0 +1 0 0
⎢+1 0 0 0⎥
η(α)(β) = η (α)(β) = ⎢⎣0
⎥. (1.52)
0 0 −1⎦
0 0 −1 0
or equivalently
Since [eμ(β) eμ(γ) ] = 0, it can be shown that the spin coefficients satisfy the
;ν
relations
Note that the calculations of the spin coefficients does not require the
calculations of the covariant derivatives. In fact, introducing the λ-symbols
via the relations
1
γ(α)(β)(γ) = [λ(α)(β)(γ) + λ(β)(γ)(α) − λ(γ)(α)(β) ]. (1.58)
2
As is evident from Eq. (1.57), the calculations of λ(α)(β)(γ) (consequently,
the calculations of γ(α)(β)(γ) ) require only the calculations of partial deriva-
tives. On the other hand, from Eq. (1.57) we find that the λ-symbols have
the properties,
where
(β)
γ(α) (γ) ≡ η (β)(δ) γ(α)(δ)(γ) , γ(α)(β)(γ),(δ) ≡ γ(α)(β)(γ);μ eμ(δ) . (1.61)
(ε)
A(α)|(β) ≡ eμ (α)Aμ;ν eν(β) = A(α),(β) − A(ε) γ(α) (β) . (1.63)
(α)
Similarly, we define A |(β) as
On the other hand, by projecting Eq. (1.49) onto the tetrad frame, we
find that the relationship among the Riemann, Weyl and Ricci tensors goes
over in the tetrad form without change, and thus is given by
1
R(α)(β)(γ)(δ) = C(α)(β)(γ)(δ) + η(α)(γ) R(β)(δ)
2
+ η(β)(δ) R(α)(γ) − η(β)(γ) R(α)(δ) − η(α)(δ) R(β)(γ)
1
+ η(α)(δ) η(β)(γ) − η(α)(γ) η(β)(δ) R, (1.67)
6
where
R ≡ η (α)(β) R(α)(β) = 2 R(0)(1) − R(2)(3) . (1.68)
Having written all of the formulas that we need in terms of the tetrad
components, we are now ready to write down the NP equations. However,
before we do so, following NP (Newman and Penrose, 1962), we introduce
the following special notations that considerably simplify the expressions
of the NP equations.
First of all, the spin coefficients are designated by
1 1
ε≡ [γ(0)(1)(0) − γ(2)(3)(0) ] = [lμ;ν nμ lν − mμ;ν mμ lν ],
2 2
1 1
α≡ [γ(0)(1)(3) − γ(2)(3)(3) ] = [lμ;ν nμ mν − mμ;ν mμ mν ],
2 2
1 1
β≡ [γ(0)(1)(2) − γ(2)(3)(2) ] = [lμ;ν nμ mν − mμ;ν mμ mν ],
2 2
1 1
γ≡ [γ(0)(1)(1) − γ(2)(3)(1) ] = [lμ;ν nμ nν − mμ;ν mμ nν ], (1.69)
2 2
and all other spin coefficients can be obtained from them by using the
symmetry given by Eq. (1.56) and the fact that the complex conjugate of
any quantity can be obtained by replacing the index 2, wherever it appears,
by the index 3, and vice versa.
As mentioned previously, the Weyl tensor has 10 independent compo-
nents at each point of spacetime. In the NP formalism, these components
are specified by five complex “scalars” as follows:
Cμνλδ ≡ −4 Ψ2 + Ψ2 [l[μ nν] l[λ nδ] + m[μ mν] m[λ mδ] ]
+ 4 Ψ2 − Ψ2 [1[μ nν] m[λ mδ] + m[μ mν] l[λ nδ] ]
− 4{Ψ0n[μ mν] n[λ mδ] + Ψ1 [l[μ nν] n[λ mδ] + n[μ mν] l[λ nδ]
+ n[μ mν] m[λ mδ] + m[μ mν] n[λ mδ] ] − Ψ2 [l[μ mν] n[λ mδ]
+ n[μ mν] l[λ mδ] ] − Ψ3 [l[μ nν] l[λ mδ] + l[μ mν] l[λ nδ]
− l[μ mν] m[λ mδ] − m[μ mν] l[λ mδ] ] + Ψ4 l[μ mν] l[λ mδ]
+ Complex conjugates}, (1.70)
where Ψ0 , . . . , Ψ4 are called the Weyl scalars, and a bar over a letter denotes
the complex conjugate, as mentioned previously.
The various terms in Eq. (1.70) have the following physical interpreta-
tions (Szekeres, 1965). The Ψ0 and Ψ1 terms represent, respectively, the
transverse and longitudinal wave components in the nμ direction, the Ψ2
term a “Coulomb” component, and the Ψ3 and Ψ4 terms the longitudi-
nal and transverse wave components in the lμ direction. By contracting
Eq. (1.70) with appropriate combinations of the null vectors, {l, n, m, m},
we find
Ψ0 ≡ −C(0)(2)(0)(2) = −Cμνλδ lμ mν lλ mδ ,
Ψ1 ≡ −C(0)(1)(0)(2) = −Cμνλδ 1μ nν lλ mδ ,
March 7, 2020 9:53 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch01 page 17
1
Ψ2 ≡ − [C(0)(1)(0)(1) − C(0)(1)(2)(3) ]
2
1
= − Cμνλδ [lμ nν lλ nδ − lμ nν mλ mδ ],
2
Ψ3 ≡ C(0)(1)(1)(3) = Cμνλδ lμ nν nλ mδ
= −Cμνλδ nμ lν nλ mδ ,
Ψ4 ≡ −C(1)(3)(1)(3) = −Cμνλδ nμ mν nλ mδ . (1.71)
Similarly, the tetrad components of the traceless Ricci tensor, S(α)(β) , can
be written by using the following notation:
1 1 1
Φ00 ≡ S(0)(0) = Sμν lμ lν = Φ00 = R(0)(0) ,
2 2 2
1 1 1
Φ01 ≡ S(0)(2) = Sμν lμ mν = Φ10 = R(0)(2) ,
2 2 2
1 1 1
Φ02 ≡ S(2)(2) = Sμν mμ mν = Φ20 = R(2)(2) ,
2 2 2
1 1
Φ11 ≡ [S + S(2)(3) ] = Sμν (lμ nν + mμ mν ) (1.72)
4 (0)(1) 4
1
= Φ11 = [R(0)(1) + R(2)(3) ],
4
1 1 1
Φ12 ≡ S(1)(2) = Sμν nμ mν = Φ21 = R(1)(2) ,
2 2 2
1 1 1
Φ22 ≡ S(1)(1) = Sμν nμ nν = Φ22 = R(1)(1) ,
2 2 2
and the trace of the Ricci tensor R(α)(β) is defined by
1 1
Λ≡− R = − [R(0)(1) − R(2)(3) ]. (1.73)
24 12
We hope that there will be no confusion between the cosmological constant
Λ used in Eq. (1.20) and the one used in Eq. (1.73), as in this book we shall
mainly consider the cases in which the cosmological constant vanishes.3
by the collision of two plane gravitational and/or electromagnetic waves (Barrabés and
Hogan, 2014, 2015; Halisoy, Mazharimousavi and Gurtug, 2014).
March 7, 2020 9:53 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch01 page 18
On the other hand, the combination of Eq. (1.67) and Eqs. (1.70)–(1.73)
gives
1
R(0)(1)(0)(1) = C(0)(1)(0)(1) − R(0)(1) + R
6
= − Ψ2 + Ψ2 − 2Φ11 + 2Λ,
1
R(0)(1)(0)(2) = C(0)(1)(0)(2) − R(0)(2) = −Ψ1 − Φ01 ,
2
1
R(0)(1)(1)(2) = C(0)(1)(1)(2) + R(1)(2) = Ψ3 + Φ12 ,
2
R(0)(1)(2)(3) = C(0)(1)(2)(3) = Ψ2 − Ψ2 ,
R(0)(2)(0)(2) = C(0)(2)(0)(2) = −Ψ0 ,
1
R(0)(2)(0)(3) = − R(0)(0) = −Φ00 ,
2
1
R(0)(2)(1)(2) = R(2)(2) = Φ02 ,
2 (1.75)
1
R(0)(2)(1)(3) = C(0)(2)(1)(3) − R = Ψ2 + 2Λ,
12
1
R(0)(2)(2)(3) = C(0)(2)(2)(3) − R(0)(2) = Ψ1 − Φ01 ,
2
R(1)(2)(1)(2) = C(1)(2)(1)(2) = −Ψ4 ,
1
R(1)(2)(1)(3) = − R(1)(1) = −Φ22 ,
2
1
R(1)(2)(2)(3) = C(1)(2)(2)(3) − R(1)(2) = Ψ3 − Φ12 ,
2
1
R(2)(3)(2)(3) = C(2)(3)(2)(3) + R(2)(3) + R
6
= − Ψ2 + Ψ2 + 2Φ11 + 2Λ.
March 7, 2020 9:53 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch01 page 19
I ≡ Rαβγδ Rαβγδ
2
= 8{3(Ψ22 + Ψ2 ) − 4 Ψ1 Ψ3 + Ψ1 Ψ3 + Ψ0 Ψ4 + Ψ0 Ψ4
− 4 Φ01 Φ12 + Φ01 Φ12 + 4Φ11 2 + 2Φ02 Φ02
+ 2Φ00 Φ22 + 12Λ2 }. (1.79)
Fig. 1.2. A null congruence meets So in the circle S. The image of the circle S
on SP is an ellipse.
By using m and m, we construct two space-like unit vectors via the relations
mμ + mμ mμ − mμ
E(2)μ = √ , E(3)μ = √ . (1.85)
2 i 2
Then, it is easy to show that
where ϕ is the angle formed by η μ and E(2)μ [Fig. 1.2(b)], and T μ is the
unit vector tangent to the circle S at O . Thus we have
E(2)μ + iE(3)μ eiϕ
mμ = √ = √ (η μ + iT μ ),
2 2
(1.87)
E(2)μ − iE(3)μ e−iϕ
mμ = √ = √ (η μ − iT μ ).
2 2
March 7, 2020 9:53 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch01 page 24
where
σ ≡ |σ|ei2ϕ0 . (1.89)
Equation (1.88) shows that if ρ and σ vanish then the change of η μ along
C is proportional to lμ . That is, the image of this null congruence on SP is
also a circle, which results from the parallel transport of the circle S from
O to P along C. However, in general, ρ and σ do not vanish. Consequently,
the first term on the right-hand side of Eq. (1.88) describes an expansion
of the circle S along C. The rate of expansion is
ρ+ρ
θ≡− . (1.90)
2
The second term describes a rotation of the circle S with a rate given by
i(ρ − ρ)
ω≡ . (1.91)
2
And the last term on the right-hand side of Eq. (1.88) depends explicitly
on the angle ϕ. It is easy to see that because of this dependence the circle
S goes over into an ellipse, and the minor axis of which forms an angle ϕ0
with respect to E(2)μ [Fig. 1.2(b)].
For the null geodesic congruence formed by nμ , the coefficients −ν, −γ,
−μ and −λ correspond to κ, ε, ρ and σ, respectively. The quantities θ, ω
and σ, as defined above, were first introduced by Sachs (1961, 1962, 1964),
and are called the optical scalars.
Finally, we note that the change of the shape of a null congruence is
due not only to the spacetime curvature, but also to the inertial field. The
latter follows from the fact that even in a flat spacetime the optical scalars
may not vanish because of the choice of the coordinate system. In order to
consider the effects only due to the spacetime curvature, we must consider
the variation of the optical scalars along geodesics (see Section 1.4). That is,
we must consider the acceleration between geodesics, instead of considering
their relative velocity.
March 7, 2020 9:53 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch01 page 25
(δ) A massless neutrino field: In general, the case for a massless neutrino
field is much more complicated than the previous ones. This is mainly due
to the fact that a neutrino field is described by a two component spinor φA ,
which satisfies the neutrino Weyl equations
μ
σA Ḃ
φA ;μ = 0, (1.100)
μ
where σA Ḃ
are the complex Pauli spin matrices (Griffiths, 1980, 1991;
Tsoubelis and Wang, 1991), and the spin indices A, B take the values 1, 2.
Then, the energy–stress tensor for a massless neutrino field takes the form
Tμν = i[σμAḂ (φA φḂ ;ν − φḂ φA ;ν ) + σνAḂ (φA φḂ ;μ − φḂ φA ;μ )]. (1.101)
OA LA = −LA OA = 1. (1.103)
In terms of Φ and Ψ and the spin coefficients, Eq. (1.98) takes the form
(Griffiths, 1976a, 1991; Tsoubelis and Wang, 1991)
Equations (1.102) and (1.103) are the basic equations for a neutrino field.
(ε) An isentropic perfect fluid: The energy–stress tensor for a perfect fluid
takes the form
where uμ is the four-velocity of the fluid, p the pressure, and μ the energy
density (where μ must not be confused with the one used for spin coeffi-
cients). Substituting Eq. (1.106) into Eq. (1.21), we obtain
μ;ν uν + (μ + p)uν;ν = 0,
(1.107)
(μ + p)uμ ;ν uν + (uμ uν − g μν ) p;ν = 0,
p = (γ − 1)μ, (1.108)
where γ is a constant, and must not be confused with the one used for spin
coefficients.
In general, matter fields need to satisfy some energy conditions. The
energy conditions for a neutrino field are discussed by Griffiths (1980).
These include the weak, dominant and strong energy conditions (Hawking
and Eillis, 1973).
March 7, 2020 9:53 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch01 page 28
(a) The weak energy condition: This condition says that the energy den-
sity measured by any observer must be non-negative. Mathematically, it is
equivalent to saying that for any time-like vector uμ we must have
T μν uμ uν ≥ 0. (1.109)
q μ qμ ≥ 0. (1.110)
(c) The strong energy condition: This is basically stated that the expan-
sion of a time-like geodesic congruence with zero vorticity will monotoni-
cally decrease along a geodesic. Mathematically, this is equivalent to require
Rμν uμ uν ≥ 0. Then, by the Einstein field equations, this implies that
1 Λ
Tμν u u ≥
μ ν
T− uν uν , (1.111)
2 κ
where uμ denotes the tangential vector of the given time-like geodesic, and
T ≡ T ν ν . It is normally said that the energy-momentum tensor Tμν satisfies
the strong energy condition, if it obeys Eq. (1.111) for Λ = 0 (Hawking and
Eillis, 1973).
equivalent to Lichnerowicz’s, and implied by, but not equivalent to, O’Brien
and Synge’s.
The non-null surface case has attracted a lot of interest with various
motivations since 1980s. In particular, the phase transitions that might
have occurred in the early epoch of the universe resulted in the formation
of topological defects, domain walls, cosmic strings, monopoles and tex-
tures, according to various field theories, including a wide variety of grand
unified theories (GUT’s) (Zel’dovish, Kobzarev and Okun, 1976; Kibble,
1976; Vilenkin, 1981). The time-varying gravitational fields associated with
the formation of cosmic strings will create particles that have contributions
to the average energy density of the Universe. These contributions could
be especially important in the early Universe (Zel’dovish, 1980; Vilenkin,
1985; Xanthopoulos, 1986a, 1986b, 1987; Economou and Tsoubelis, 1988a,
1988b; Tsoubelis, 1989a; Letelier and Wang, 1995; Wang and Santos, 1996;
Wang and Nogales, 1997; Nogales and Wang, 1998; Bronnikov, Santos and
Wang, 2019). In particular, they were once believed to provide a mecha-
nism to produce the cosmic microwave background (CMB) and large-scale
structure observed in our Universe (Vilenkin and Shellard, 2000). However,
later observations of CMB ruled out cosmic strings formed in the context
of symmetry breaking in GUT’s as the sources of the cosmological pertur-
bations (Bennett et al., 1996; Spergel et al., 2007), and led to an upper
bound (Ade et al., 2016),
Gμ
10−7 , (1.112)
c2
where μ denotes the string’s tension.
Nevertheless, the subject has attracted recently lots of attention again
in the framework of string/M-Theory, the so-called cosmic superstrings
(Dvali and Vilenkin, 2004; Copeland, Myers and Polchinski, 2004), which
were formed before inflation took place and stretched to macroscopic length
scales in the inflationary phase. During the subsequent epochs, a compli-
cated network of various string elements forms (Chernoff and Tye, 2018; and
references therein). The main phenomenological consequence of a string
network is the emission of GWs (Ringeval and Suyama, 2017; and refer-
ences therein), generating bursts at cusps, kinks and junctions, as well as a
stochastic gravitational wave background. Low tension strings are natural
in string/M-theory, and can easily satisfy the observational bounds given
above (Chernoff, Flanagan and Wardell, 2018).
On the other hand, domain walls played an important role in the the-
ory of inflation (Linde, 1984; Ipser and Sikivie, 1984; Brandenberger, 1985;
March 7, 2020 9:53 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch01 page 30
The surface Σ divides Ω into two open regions, Ω+ = {xμ : ϕ (xμ ) > 0} and
Ω− = {xμ : ϕ (xμ ) < 0}. We assume that: (i) the restrictions gμν
±
≡ gμν |Ω±
3 0
are at least C and (ii) gμν is C across Σ.
Note that for the sake of convenience in the following we consider only
the cases in which the spacetime has only one singular surface, but the
extension to the cases with multiple singular surfaces is straightforward.
The normal vector to Σ is defined as
∂ϕ
kμ ≡ = ϕ;μ . (1.114)
∂xμ
March 7, 2020 9:53 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch01 page 32
Then, we have
∂θ(ϕ)
= δ(ϕ)kμ , (1.116)
∂xμ
where δ(ϕ) denotes the Dirac delta function with support on Σ. Thus, for
a test function f of compact support, we have
√
δ(ϕ), f ≡ −g δ(ϕ)f dV = f dS = − f dS, (1.117)
Ω ∂Ω− ∂Ω+
and the following relations are valid (Gelfand and Shilov, 1964):
4 Note that instead of defining a step function as Eq. (1.115), we can use the Heaviside
step function, H(ϕ), which is unity for the non-negative arguments and otherwise zero.
But all the following results are valid for both of them, if we just simply replace one
by another. In fact, the specific value of θ(ϕ) at ϕ = 0 is irrelevant in the sense of
distributions, and does not affect the following results at all.
March 7, 2020 9:53 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch01 page 33
Then, we can see that the product F δ(ϕ) is well defined whenever F is
C 0 , and it depends only on the restriction F |Σ of F to Σ. More generally,
the product F δ (n) (ϕ) is well defined, provided that F is C n . Clearly, such
defined F δ (n) (ϕ) depends only on the values of F and its partial derivatives
of the mth order evaluated on Σ, where m ≤ n.
By using the continuity of gμν , we can write gμν as
+ −
gμν = θ(ϕ)gμν + (1 − θ(ϕ))gμν . (1.121)
Then, for a test function f we find that the distribution derivative of gμν is
+ −
∂λ gμν , f = ∂λ θ(ϕ)gμν + (1 − θ(ϕ))gμν ,f
+ −
= {θ(ϕ)gμν,λ + (1 − θ(ϕ))gμν,λ }, f
+ −
+ {gμν ∂λ θ(ϕ) + gμν ∂λ (1 − θ(ϕ))}, f . (1.122)
Since
+ −
gμν ∂λ θ(ϕ) + gμν ∂λ (1 − θ(ϕ)) , f
+
−
= ∇λ θ(ϕ), gμν f + ∇λ (1 − θ(ϕ)), gμν f
+ −
=− ∇λ gμν f dV − ∇λ gμν f dV
Ω+ Ω−
+ −
=− gμν f dSλ − (gμν f )dSλ
∂Ω+ ∂Ω−
+ −
=− gμν − gμν f dSλ = 0, (1.123)
∂Ω+
we have
+
∂λ gμν , f = {θ(ϕ)gμν,λ + (1 − θ(ϕ))gμν,λ }, f . (1.124)
On the other hand, from the definition of the connection coefficients Γμνλ
we find that
Γμνλ = θ(ϕ)Γ+μ −μ
νλ + (1 − θ(ϕ))Γνλ , (1.126)
(T μ ) ≡ θ(ϕ)T +μ + (1 − θ(ϕ))T −μ ,
D
(1.127a)
μ D
T ;ν ≡ θ(ϕ)T +μ ;ν + (1 − θ(ϕ))T −μ ;ν . (1.127b)
It must be noted that, since for any given test function f , we have
θ(ϕ)(1 − θ(ϕ)), f = [θ(ϕ)(1 − θ(ϕ))f ]dV
Ω
= (1 − θ(ϕ))f dV = 0, (1.131)
Ω+
so we can set
g ± μν = g ± μν (ϕ (xμ ), xμ ). (1.133)
Thus, we have
± ±
∂gμν 1 ∂ 2 gμν
±
gμν = 0
gμν + ϕ+ ϕ2 + O ϕ3 . (1.134)
∂ϕ ϕ=0 2 ∂ϕ2
ϕ=0
March 7, 2020 9:53 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch01 page 35
− −
where γμν ≡ [gμν,ϕ ] and γ̂μν ≡ [gμν,ϕϕ ] . Hence, we have
− 1
[Γμνλ ] = [kν γ μ λ + kλ γ μ ν − k μ γνλ ]. (1.136)
2
From the definition of the Riemann tensor given by Eq. (1.13), on the
other hand, we find
where
− −
H σ μνλ ≡ kν Γσμλ − kλ Γσμν ,
− − (1.138)
− −
I σ μνλ ≡ [Γσλδ ] Γδμν − [Γσνδ ] Γδμλ ,
and R±σ μνλ are the Riemann tensor calculated, respectively, in Ω± . Then,
the Ricci tensor is given by
where
R± μν ≡ R±σ μσν ,
1
Hμν ≡ H σ μσν = kμ kλ γ λ ν + kν kλ γ λ μ − k λ kλ γμν − kμ kν γ λ λ ,
2
1
Iμν ≡ I σ μσν = kμ kν γ λδ γλδ + 2kλ kδ γ λ μ γ δ ν − (kμ γ λ ν + kν γ λ μ )kλ γ δ δ
4
− 2γμδ γ δ ν − γ δ δ γμν kλ k λ .
(1.140)
The Ricci scalar is given by
where
R± ≡ R±ν ν ,
H ≡ H ν ν = kμ kν γ μν − k ν kν γ μ μ ,
1 1 μ δ
I ≡ Iν ν = γμδ γ δ ν − γ δ δ γμν k μ k ν + γμ γ δ − γ δμ γμδ kν k ν .
2 4
(1.142)
Now, we are ready to write down the generalized Einstein field equations,
which take the form,
1 1
Rμν − gμν R = GD μν + δ(φ) Hμν − gμν H
2 2
1
+ θ(ϕ)(1 − θ(ϕ)) Iμν − gμν I
2
= Tμν
D
+ δ(ϕ)τμν + θ(ϕ)(1 − θ(ϕ))Jμν , (1.143)
where
GD μν ≡ θ(ϕ)G+ μν + (1 − θ(ϕ))G− μν ,
(1.144)
T D μν ≡ θ(ϕ)T + μν + (1 − θ(ϕ))T − μν ,
and τμν and Jμν denote the energy–stress tensors with supports only on
the surface Σ. Equation (1.143) can be written as
G± μν = T ± μν , (1.145a)
1
Hμν − gμν H = τμν , (1.145b)
2
1
Iμν − gμν I = Jμν . (1.145c)
2
From Eq. (1.145b), we find
1
τμν = [γ δ δ k λ kλ gμν − kμ kν + kμ γ λ ν + kν γ λ μ kλ
2
− k λ kλ γμν − gμν kδ kλ γ δλ ], (1.146)
which is applicable for any kind of thin shells. Contracting it with k ν we
obtain
τμν k ν = 0. (1.147)
In the case of a null shell, Eq. (1.146) reduces to
1
τμν = [(kμ γ λ ν + kν γ λ μ )kλ − γ δ δ kμ kν − gμν kδ kλ γ δλ ]. (1.148)
2
March 7, 2020 9:53 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch01 page 37
Chapter 2
39
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch02 page 40
minor axis along λμ(2) . Taking into account the direction of relative accel-
erations, we find that λμ(a) are characterized as eigendirections (Ehlers and
Kundt, 1962). We call λμ(2) the direction of polarization with respect to uμ
(Ehlers and Kundt, 1962). The spatial projection of k μ is orthogonal to
the plane of relative accelerations, and (uμ kμ )2 is the magnitude of rela-
tive accelerations for neighboring test particles. That is, k μ determined by
Eq. (2.3) describes the direction of the propagation of the null field, and
(uμ kμ )2 represents the strength of the null field (measured by the observer
uμ ). Thus, considering Eq. (1.88) we have the following theorem.
Theorem 2.3 (Ehlers and Kundt, 1962). Null fields are character-
ized as purely transverse vacuum fields. There exists a null vector k μ such
that relative accelerations and relative rotations of inertial directions are
orthogonal to the spatial projection of k μ .
k μ = k 1/2 k̂ μ . (2.8)
then the metric for a pp-wave can be written in the following form (Ehlers
and Kundt, 1962)
(1 − A)2 2
ds2 = dt − (1 + A)4 (dx2 + dy 2 + dz 2 ), (2.14)
(1 + A)2
with A ≡ m/(2r), where r ≡ x2 + y 2 + z 2 . Starting with this form of
metric, Aichelburg and Sexl (1971) were able to obtain the gravitational
field produced by a massless particle, after first making a Lorentz boost
and then taking the massless limit. In doing so, they were able to show
that such an obtained solution takes exactly the same form as that given
by Eq. (2.10) for a plane-fronted gravitational wave.
To show the above claim, let us first consider the Lorentz boost along
the x-direction,
where γ[≡ 1/(1−v 2 )1/2 ] is the Lorentz factor. Then, the above metric takes
the form,
where
m p(1 − v 2 )
A≡ = , (2.17)
2r 2[(x̄ − v t̄ ) + (1 − v 2 )(ȳ 2 + z̄ 2 )]1/2
2
with m ≡ p(1 − v 2 )1/2 . The above metric is not well defined along the light
cone when we take the limit v → 1. To overcome this problem, Aichelburg
and Sexl introduced the new coordinates t and x via the relations,
1
− (dt − dx )2 , (2.19)
(x − vt )2 + (1 − v 2 )
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch02 page 44
2 2
where ρ2 ≡ y + z . Using the limit,
1 1
lim −
v→1 (x − vt )2 + ρ2 (1 − v 2 ) (x − vt )2 + (1 − v 2 )
= −2ln(ρ)δ(x − t ), (2.20)
we find that
2 2 2 2
lim ds2 = dt − dx − dy − dz
v→1
2 2
+ 4pδ(t − x ) ln(y + z )(dt − dx )2 , (2.21)
which takes precisely the form of Eq. (2.10) by setting (û, Z, X, Y ) =
(t − x , (t + x )/2, y , z ).
and
M = ln(AB). (2.25)
From Eq. (1.69), it can be shown that the non-vanishing spin coefficients
are given by
1 1
μ = − AU,u , γ= A[(ln B),u − (1/2)iV,u sinh W ],
2 2
(2.26)
1
λ = A[V,u cosh W + iW,u ],
2
where U,u ≡ ∂U/∂u, and so on. Then, from Eqs. (1.77a)–(1.77r) we find
that the non-vanishing Weyl and Ricci scalars are given, respectively, by
Φ22 = 0, (2.28)
where Φ22 is given explicitly in terms of the metric coefficients by Eq. (3.8d)
in Chapter 3.
On the other hand, from Eqs. (1.82) and (2.26) we see that in the
present case the null vector nμ defines a null geodesic congruence, and if B
is chosen to be constant, then the null geodesics are affinely parameterized.
Thus, when M , U , V and W are functions of the null coordinate u only,
the Petrov type N plane gravitational wave represented by Ψ4 propagates
along the null geodesics (see Fig. 2.1).
The Weyl tensor given by Eq. (1.70) now takes the following form:
Cμνλσ = −4{Ψ4l[μ mν] l[λ mσ] + Ψ4 l[μ mν] l[λ mσ] }. (2.29)
(1) (2)
If we define mμν and mμν as (Wang, 1991f)
m(1)
μν = 2l[μ E(2)ν] , m(2)
μν = 2l[μ E(3)ν] , (2.30)
where E(2)μ and E(3)μ are given by Eq. (1.85), we find that Eq. (2.29) can
be written in the following form:
1 (1) (2) (2)
Cμνλσ = − {[m(1)
μν mλσ − mμν mλσ ](Ψ4 + Ψ4 )
2
(2) (1)
+ i[m(1) (2)
μν mλσ + mμν mλσ ](Ψ4 − Ψ4 )}. (2.31)
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch02 page 46
Fig. 2.1. The spacetime for a plane gravitational wave. For a sandwich wave,
the matter can be arranged so that the regions u < 0, and u > u0 are flat, and
the hypersurface u = 0 is the leading wavefront of the wave.
we find that
(1) (1) (2)
mμν = cos ϕ4 mμν + sin ϕ4 mμν ,
(2) (1) (2)
(2.33)
mμν = − sin ϕ4 mμν + cos ϕ4 mμν ,
and
e+μνλσ = cos 2ϕ4 e+μνλσ + sin 2ϕ4 e×μνλσ ,
(2.34)
e×μνλσ = − sin 2ϕ4 e+μνλσ + cos 2ϕ4 e×μνλσ ,
where
(1) (1) (2) (2)
e+μνλσ ≡ mμν mλσ − mμν mλσ ,
(1) (2) (2) (1)
(2.35)
e×μνλσ ≡ mμν mλσ + mμν mλσ .
Thus, in terms of E(2)μ and E(3)μ , Eq. (2.31) reads
D2 ημ μ μ
= (Ψ0 Ψ0 )1/2 (nλ uλ )2 [E(2) ν
E(2) − E(3) ν
E(3) ]ην , (2.45)
Dτ 2
but now with the angle ϕ0 being defined as
Im(Ψ0 )
tan 2ϕ0 = − . (2.46)
Re(Ψ0 )
Similar to Eq. (2.43), we have
μ μ
E(2);ν nν = 0 = E(3);ν nν . (2.47)
Thus, in the latter case the metric (2.22) describes a plane gravitational
wave represented by Ψ0 with the polarization angle ϕ0 , which moves along
the null geodesic congruence defined by lμ . The angle ϕ0 is constant along
the path, along which the Ψ0 -wave propagates.
When the angle ϕ0 (ϕ4 ) is constant everywhere, we say the correspond-
ing plane gravitational wave is constantly polarized.
Harte, 2013). In addition, the passage of the GW affects not only the
position of the test particle, but also its velocity. In fact, the change of
the velocity of the particle is also permanent (Souriau, 1973; Braginsky
and Thorne, 1987; Bondi, 1957; Bondi and Pirani, 1989; Grishchuk and
Polnarev, 1989; Zhang, Duval, Gibbons and Horvathy, 2018).
When far from the sources, the emitted GWs can be well approximated
by plane GWs. The spacetimes for plane GWs can be cast in various forms,
depending on the choice of the coordinates and gauge-fixing, as showed in
Sections 2.2 and 2.4. The form of Eq. (2.22) was originally due to Baldwin,
Jeffery and Rosen (BJR) (Baldwin and Jeffery, 1926; Rosen, 1937). Despite
its several attractive features, the system of the BJR coordinates is often
singular within a finite width of a wave, and when studying the asymp-
totic behavior of the spacetime, extensions beyond this singular surface are
needed.
In this section, we point out that there exist actually two kinds of sin-
gularities in plane gravitational wave spacetimes, one represents coordinate
singularities, which can be removed by proper coordinate transformations,
and the other represents really spacetime singularities, and physical quanti-
ties, such as distortions of test particles, become infinitely large when such
singularities are approaching (Wang et al., 2018). Therefore, in the latter
these singularities already represent the boundaries of the spacetimes and
extensions beyond them are not only impossible but also not needed. Since
gravitational memory effects and soft graviton theorems are closely related
to the asymptotical behaviors of plane GW spacetimes, in the latter the
spacetimes cannot be used to study such properties.
In GR, there are powerful Hawking–Penrose theorems (Hawking and
Ellis, 1973) from which one can see that spacetimes with “physically rea-
sonable” conditions are singular. However, the theorems did not tell the
nature of the singularities, and Ellis and Schmidt (1977) classified them
into two different kinds: spacetime curvature singularities and coordinate
singularities. The former is real and cannot be removed by any coordinate
transformations of the kind,
xμ → x = ζ μ (xν ),
μ
μ, ν = 0, 1, 2, 3, (2.48)
χ + ω 2 χ = 0, (2.50)
where
1
χ ≡ e−U/2 , ω≡ V . (2.51)
2
Then, from Eq. (2.50) we can see that, for any given initial value, χ0 > 0,
there always exists a moment, say, u = us at which χ vanishes,
where α > 0, and χ̂(us ) = 0. The function χ̂(u) in general takes the form,
∞
χ̂(u) = χn (u − us )n , (2.54)
n=0
with χ0 = 0. Then, one can consider the tetrad, eμ(a) (a = 0, 1, 2, 3), defined
by
1 μ 1 μ
eμ(0) = γ0 δuμ + δ , eμ(1) = γ0 δuμ − δ ,
2γ0 v 2γ0 v (2.55)
U −V U +V
eμ(2) = e 2 δyμ , eμ(3) = e 2 δzμ ,
that is, they are unit orthogonal vectors and parallelly transported along the
time-like geodesics, defined by eμ(0) ≡ dxμ /dλ, where λ denotes the proper
time of the time-like geodesics, so that they form a freely falling frame (Ori,
2000; Nolan, 2000; Hirschmann, Wang and Wu, 2004; Sharma, Tziolas,
Wang and Wu, 2011). Projecting the Riemann tensor onto this frame, one
obtains some non-zero components of the Riemann tensor R(a)(b)(c)(d) . If
one integrates them twice along the time-like geodesics, which gives the
distortions, one finds that such integral always diverges as u → us , except
for the cases (Wang et al., 2018),
1
(i) α = or (ii) α = 1. (2.57)
2
Therefore, all the plane GW spacetimes are singular physically at the focused
point u = us , exceptions are only the ones with α = 1/2 or 1. As a result, all
the plane GW spacetimes cannot be used to study memory effects and soft
graviton theorems, except these two particular cases, as only these space-
times can be possibly extended to null infinity, whereby memory effects and
soft graviton theorems can be studied.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch02 page 52
The above results, although very simple, may have profound implica-
tions on the studies of gravitational memory effects (Favata, 2010; Bieri,
Garfinkle and Yunes, 2017) and “soft-graviton” theorems (Hawking, Perry
and Strominger, 2016, 2017; Strominger, 2017), as both of them are con-
cerned with the asymptotic behavior of the spacetimes at infinities. But
when the spacetimes become singular at the finite focusing point, such
infinities do not exist. Therefore, only the non-singular spacetimes are
relevant to them.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch03 page 53
Chapter 3
53
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch03 page 54
null coordinates, the collision of two such plane gravitational waves in the
Minkowski background can be always considered as head-on collision. Then,
the spacetimes can be arranged as follows.
For u < 0 and v > 0 (see Fig. 3.1), we assume the metric to be that
of a plane wave described by the metric (2.22) with M, U, V and W being
functions of v only. Without loss of generality, we can always choose the
coordinate v so that the hypersurface v = 0 represents the wavefront. For
a sandwich wave, we may further assume that the spacetime is flat (cor-
respondingly, M, U, V and W are constants) in the region v > v0 . In the
region v < 0 and u > 0, on the other hand, the metric is again assumed to
be that of a plane gravitational wave, but with M, U, V and W being func-
tions of u only, and the u = 0 hypersurface is the wavefront. For a sandwich
wave, it is further assumed that the region u > u0 is flat. Assuming that
the collision occurs in a flat background, then in the region u, v < 0 the
spacetime is Minkowski, since in this region the two coming waves have not
arrived, yet. In the region u > 0 and v > 0 where the two plane waves
interact, the metric coefficients M, U, V and W become functions of both
u and v. Hence, a characteristic initial value problem has been set up with
the data posed on a pair of null hypersurfaces u = 0 and v = 0 intersecting
in a space-like two-dimensional surface u = 0 = v. In principle, these initial
conditions uniquely determine the geometry in the wave interacting region
u, v > 0 (Szekeres, 1972; Xanthopoulos, 1986a; Hauser and Ernst, 1989a,
1989b, 1990). For the sake of convenience, the above various regions are
numbered as follows (see Fig. 3.1):
Fig. 3.1. The projection of a spacetime for colliding plane waves onto the (u, v)-
plane.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch03 page 55
Region I (u, v < 0): This is the region in which the two colliding plane
waves have not arrived, yet, so the spacetime in this region is flat, and the
functions M, U, V and W are constants.
Region II (u < 0, v > 0): This is the region in which a plane grav-
itational wave propagates toward the right-hand side with the v = 0
hypersurface as its wavefront. The functions M, U, V and W depend on v
only.
Region III (u > 0, v < 0): In this region, an opposite moving plane grav-
itational wave is incident, and the u = 0 hypersurface is its wavefront. In
this region M, U, V and W are functions of u only.
Region IV (u, v > 0): In this region two incoming waves interact, and
M, U, V and W are functions of both u and v.
In addition to these four regions, there are also two null hypersurfaces,
Σu and Σv , defined as Σu ≡ {xμ : u = 0} and Σv ≡ {xμ : v = 0},
respectively. At the space-like two-dimensional surface u = 0 = v, the
two plane waves collide. Across the hypersurfaces u = 0 and v = 0, the
Riemann curvature tensor, in general, suffers a shock and or an impulsive
type of discontinuity. For the collision of two sandwich waves, Regions II
and III reduce, respectively, to u < 0 and 0 < v < v0 , and 0 < u < u0 and
v < 0.
Corresponding to the metric (2.22), the non-vanishing connection coef-
ficients Γμαβ defined by Eq. (1.8) are given by
1 M V −U
Γuuu = −M,u , Γu22 = e (e cosh W ),v ,
2
1 1 M −V −U
Γu23 = − eM (e−U sinh W ),v , Γu33 = e (e cosh W ),v ,
2 2
1 M V −U
Γvvv = −M,v , Γv22 = e (e cosh W ),u ,
2
1 1 M V −U
Γv23 = − eM (e−U sinh W ),u , Γv33 = e (e cosh W ),u ,
2 2
1 1
Γ22u = (cosh2 W V,u − U,u ), Γ22v = (cosh2 W V,v − U,v ),
2 2
1 1
Γ33u = − (cosh2 W V,u + U,u ), Γ33v = − (cosh2 W V,v + U,v ),
2 2
1
Γ23u = − e−V (sinh W cosh W V,u + W,u ),
2
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch03 page 56
1
Γ23v = − e−V (sinh W cosh W V,v + W,v ),
2
1
Γ32u = eV (sinh W cosh W V,u − W,u ),
2
1
Γ32v = eV [sinh W cosh W V,v − W,v ). (3.1)
2
When W = 0, the two incoming gravitational waves have fixed polarization
directions, and are said to be collinearly (linearly) polarized. Otherwise,
they are said to be non-collinearly (nonlinearly) polarized.
1 1 1
ρ= BU,v , μ = − AU,u , ε = − B[2(ln A),v + iV,v sinh W ],
2 2 4
1 1
γ = A[2(ln B),u − iV,u sinh W ], σ = − B(V,v cosh W − iW,v ), (3.2)
4 2
1
λ = A(V,u cosh W + iW,u ),
2
and
κ = ν = τ = π = α = β = 0. (3.3)
The vanishing of the spin coefficients κ and ν implies [see Eq. (1.82)] that
the null vectors nμ and lμ defined by Eq. (2.23) still define null geodesic
congruences even in the interacting region (Region IV), as in the single
plane gravitational wave case. Moreover, if the function A is chosen to be
constant, then the null geodesics defined by lμ are affinely parameterized,
while when the function B is chosen to be constant, the null geodesics
defined by nμ are affinely parameterized.
The plane gravitational waves represented by Ψ0 and Ψ4 propagate
along the null geodesic congruences defined, respectively, by nμ and lμ .
Then, Eqs. (3.2), (1.90) and (1.91) show that after the collision the two null
geodesic congruences are no longer shear-free, although their rotations are
still zero. Therefore, after the collision, the two plane gravitational waves
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch03 page 57
will in general pick up shear. Substituting Eq. (3.3) into Eqs. (1.77a)–
(1.77r), we find
1
Ψ0 = − B 2 {V,vv cosh W + (M,v − U,v )V,v cosh W + 2 sinh W V,v W,v
2
− i[W,vv + (M,v − U,v )W,v − sinh W cosh W V,v 2 ]}, (3.7a)
1
Ψ2 = AB M,uv − U,uv + W,u W,v + cosh2 W V,u V,v
6
3
+ i cosh W (V,u W,v − V,v W,u ) , (3.7b)
2
1
Ψ4 = − A2 {V,uu cosh W + (M,u − U,u )V,u cosh W + 2 sinh W V,u W,u
2
+ i[W,uu + (M,u − U,u )W,u − sinh W cosh W V,u2 ]}, (3.7c)
and
1 2
Φ00 = B {2U,vv − U,v 2 + 2U,v M,v − W,v 2 − cosh2 W V,v 2 }, (3.8a)
4
1
Φ11 = AB{2M,uv + U,u U,v − W,u W,v − cosh2 W V,u V,v }, (3.8b)
8
1
Φ02 = AB{2 cosh W V,uv − cosh W (U,u V,v + V,u U,v )
4
+ 2 sinh W (V,u W,v + V,v W,u ) − i[2W,uv − (U,u W,v + W,u U,v )
− 2 sinh W cosh W V,u V,v ]}, (3.8c)
1 2
Φ22 = A {2U,uu − U,u 2 + 2U,u M,u − W,u 2 − cosh2 W V,u 2 }, (3.8d)
4
1
Λ = − AB{2M,uv + 4U,uv − 3U,u U,v − W,u W,v
24
− cosh2 W V,u V,v }. (3.8e)
On the other hand, from Eqs. (1.78a)–(1.78k), (3.3) and (3.5) we find
that the Bianchi identities now read
where
1A similar expression was also applied to domain walls (Wang, 1991c, 1992e, 1993;
Schmidt and Wang, 1993; Letelier and Wang, 1994) and brane worlds in string/M-theory
(Gong, Wang and Wu, 2008; Tziolas and Wang, 2008; Wang and Santos, 2008; Wu et al.,
2008; Wu, Gong and Wang, 2009; Tziolas, Wang and Wu, 20009; Devin et al., 2009; Wang
and Santos, 2010; Wang, 2010; Sharma et al., 2011), which turns out to be very useful
in the studies of singular surfaces.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch03 page 61
− F (0, v)H(v)[f + (0, v) − f − (0, v)]dSu
∂Ω+u
− F (0, v)f (0, 0)[1 − H(v)]dSu , (3.15)
∂Ω+u
f (0, 0) = 0, (3.17)
or equivalently
f,u (uH(u), vH(v)) = f,u (u, v)H(u)H(v) + f,u (u, 0)[1 − H(v)]H(u)
= f,u (u, vH(v))H(u). (3.19)
where
then the non-vanishing spin coefficients given by Eq. (3.2) can be written
as
1
ρ= B(uH(u), v)U,v (uH(u), v)H(v),
2
1
μ = − A(u, vH(v))U,u (u, vH(v))H(u),
2
1
ε = − B(uH(u), v) [ln A(uH(u), v)],v
2
i
+ V,v (uH(u), v) sinh W (uH(u), v) H(v),
2
1
γ = A(u, vH(v)) [ln B(u, vH(v))],u
2 (3.23)
i
− V,u (u, vH(v)) sinh W (u, vH(v)) H(u),
2
1
σ = − B(uH(u), v){V,v (uH(u), v) cosh W (uH(u), v)
2
− iW,v (uH(u), v)}H(v),
1
λ= A(u, vH(v)){V,u (u, vH(v)) cosh W (u, vH(v))
2
+ iW,u (u, vH(v))}H(u).
and
where ΨIV IV
i (u, v) and Φij (u, v) are the corresponding Weyl and Ricci scalars
in Region IV, obtained from Eqs. (3.7a)–(3.8e), and ΨIm Im
0 . . . Φ22 are the
impulsive part of the Riemann tensor with support on the hypersurfaces
u = 0 and v = 0, respectively, and defined by
1 2
ΨIm
0 (uH(u)) ≡ − B (uH(u), 0)[V,v (uH(u), 0) cosh W (uH(u), 0)
2
− iW,v (uH(u), 0)],
(3.26)
1
ΨIm
4 (vH(v)) ≡ − A2 (0, vH(v))[V,u (0, vH(v)) cosh W (0, vH(v))
2
+ iW,u (0, vH(v))],
and
1 2
ΦIm
00 (uH(u)) ≡ B (uH(u), 0)U,v (uH(u), 0),
2
(3.27)
1 2
ΦIm
22 (vH(v)) ≡ A (0, vH(v))U,u (0, vH(v)).
2
On the other hand, it can be shown that for any given C 1 -test function
F (u, v) we have
[1 − H(u)]H(u), F (u, v) = F (u, v)[1 − H(u)]H(u)dV
Ω
= F (u, v)[1 − H(u)]dV = 0,
Ω+u
H 2 (u), F (u.v) = F (u, v)H(u)H(u)dV
Ω
= F (u, v)H(u)dV = F (u, v)dV,
Ω+u Ω+u
2
[1 − H(u)] , F (u, v) = F (u, v)[1 − H(u)][1 − H(u)]dV
Ω
= F (u, v)[1 − H(u)]dV
Ω−u
= F (u, v)dV, (3.28)
Ω−u
or equivalently
ΨIV IV IV
4 (u > 0, vH(v)) = Ψ4 (u > 0, v > 0)H(v) + Ψ4 (u > 0, 0)[1 − H(v)],
ΨIm Im Im
4 (vH(v)) = Ψ4 (v > 0)H(v) + Ψ4 (0)[1 − H(v)], (3.30)
etc., where
ΨIV IV
4 (u > 0, 0) = lim+ Ψ4 (u > 0, v),
v→0
(3.31)
ΨIm
4 (0) = lim+ ΨIm
4 (v).
v→0
Since ΨIV Im 1 2
4 (u > 0, v > 0) and Ψ4 (v > 0) are at least C and C , respec-
tively, we can see that Ψ4 (u > 0, vH(v)) and Ψ4 (vH(v)) are C 0 across
IV Im
the hypersurface v = 0.
Similarly, it can be shown that ΦIV Im
22 (u > 0, vH(v)) and Φ22 (vH(v))
are C 0 , too, across the hypersurface v = 0, while ΨIV0 (uH(u), v > 0),
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch03 page 65
ΨIm IV Im 0
0 (uH(u)), Φ00 (uH(u), v > 0) and Φ00 (uH(u)) are C across the hyper-
surface u = 0. Thus, we have
+ΨIm Im
4 (v)H(v)δ(u) + Ψ4 (0)[1 − H(v)]δ(u)}, F (u, v)
= H(u)H(v)B∂v [ΨIV Im
4 (u, v)] + H(v)δ(u)B∂v [Ψ4 (v)], F (u, v)
= H(u)H(v)D[ΨIV
2 (u, v)]
+ ΨIV
2 (u, v)H(u)DH(v), F (u, v), (3.32)
or equivalently
Substituting Eq. (3.23) into Eqs. (3.9a)–(3.10b), and taking Eq. (3.33)
into account, we finally obtain
ΔΨIV IV IV IV IV
0 + DΦ02 = (4γ − μ)Ψ0 + 3σΨ2 + (2ε − 2ε + ρ)Φ02
+ 2σΦIV IV
11 − λΦ00 ,
3DΨIV IV IV IV IV IV
2 − 2DΦ11 + ΔΦ00 = −3λΨ0 + 9ρΨ2 + (μ + 2γ + 2γ)Φ00
− 2ρΦIV IV IV
11 − 2σΦ20 + σΦ02 ,
(3.34)
3ΔΨIV IV IV IV IV IV
2 − 2ΔΦ11 + DΦ22 = 3σΨ4 − 9μΨ2 − (ρ + 2ε + 2ε)Φ22
+ 2μΦIV IV IV
11 + 2λΦ02 − λΦ20 ,
DΨIV IV IV IV IV
4 + ΔΦ20 = (ρ − 4ε)Ψ4 − 3λΨ2 − (μ + 2γ − 2γ)Φ20 ,
− 2λΦIV IV
11 + σΦ22 , (u, v > 0),
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch03 page 66
and
DΦIV IV
11 + ΔΦ00 + 3DΛ
IV
= 2(γ + γ − μ)ΦIV IV
00 + 4ρΦ11
+σΦIV IV
02 + σΦ20 ,
(3.35)
ΔΦIV IV
11 + DΦ22 + 3ΔΛ
IV
= 2(ρ − ε − ε)ΦIV IV
22 − 4μΦ11
−λΦIV IV
20 − λΦ02 , (u, v > 0),
ΔΨIm IV Im Im
0 + BΦ02 = (4γ − μ)Ψ0 − λΦ00 ,
ΔΦIm IV IV Im Im
00 + 3BΨ2 − 2BΦ11 = (μ + 2γ + 2γ)Φ00 − 3λΨ0 , (3.36)
ΔΦIm IV
00 + BΦ11 + 3BΛ
IV
= 2(γ + γ − μ)ΦIm
00 , (u > 0, v = 0),
DΨIm IV Im Im
4 + AΦ20 = −(4ε − ρ)Ψ4 + σΦ22 ,
DΦIm IV IV Im Im
22 + 3AΨ2 − 2AΦ11 = −(ρ + 2ε + 2ε)Φ22 + 3σΨ4 , (3.37)
DΦIm IV
22 + AΦ11 + 3AΛ
IV
= −2(ε + ε − ρ)ΦIm
22 , (u = 0, v > 0),
along the null surface u = 0 with v > 0. Inside Regions I–III, the Bianchi
identities are satisfied identically.
Equation (3.34) shows that the shock parts of the two plane gravita-
tional waves represented, respectively, by ΨIV IV
0 and Ψ4 , interact with each
other through the Coulomb field Ψ2 . The components ΦIV
IV IV
02 and Φ11 of the
matter field interact with both Ψ0 and Ψ4 . That is, Φ02 and ΦIV
IV IV IV
11 are
gravitationally active to both Ψ0 and Ψ4 (Szekeres, 1965). Similarly, ΦIV
IV IV
00
is gravitationally active only to ΨIV IV IV
0 , and Φ22 only to Ψ4 . The component
ΛIV is gravitationally inert to both ΨIV 0 and Ψ4 .
IV
On the other hand, from Eqs. (3.36) and (3.37) we can see that the
component ΦIV Im Im Im
02 is gravitationally active to both Ψ0 and Ψ4 , Φ00 only to
Im Im Im
Ψ0 , and Φ22 only to Ψ4 , while all the other components of the matter
field are gravitationally inert to both ΨIm 0 and Ψ4 .
Im
From Eqs. (3.36) and (3.37) we can also see that the interaction between
purely gravitational impulsive waves is carried out through the Coulomb
field ΨIV2 (Khan and Penrose, 1971; Nutku and Halil, 1977; Wang, 1991f,
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch03 page 67
e−U = t. (3.40)
e−M
f≡ . (3.42)
8nmu2n−1 v 2m−1
which was originally found by Ernst (1968a, 1968b), where the Ernst poten-
tial E is defined via the relations
Z −1 eV sinh W eV
Z ≡ χ + iq2 , E≡ , χ≡ , q2 ≡ , (3.44)
Z +1 cosh W cosh W
η = un X + v m Y, μ = un X − v m Y,
(3.45)
X = (1 − v 2m )1/2 , Y = (1 − u2n )1/2 .
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch03 page 69
Taking the first and second derivatives of η and μ with respect to u and v,
we obtain the following useful relations:
nun−1 mv m−1
η,u = (XY − un v m ), η,v = (XY − un v m ),
Y X
n(n − 1)un−2 n2 v m u2n−2
η,uu = (XY − un v m ) − ,
Y Y3
m(m − 1)v m−2 m2 un v 2m−2
η,vv = (XY − un v m ) − ,
X X3
nmun−1 v m·1
η,uv = − η, (3.46)
XY
nun−1 mv m−1
μ,u = (XY + un v m ), μ,v = − (XY + un v m ),
Y X
n(n − 1)un−2 n2 v m u2n−2
μ,uu = (XY + un v m ) + ,
Y Y3
m(m − 1)v m−2 m2 un v 2m−2
μ,vv =− (XY + un v m ) − ,
X X3
nmu v
n−1 m−1
μ,uv = μ.
XY
In addition, we also have
Thus, we have
tη tμ
t,η = − , t,μ = − , z.η = μ, z,μ = η, (3.50)
1 − η2 1 − μ2
and
dη 2 dμ2
dt2 − dz 2 = η 2 − μ2 2
− . (3.51)
1−η 1 − μ2
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch03 page 70
Substituting Eqs. (3.44), (3.49) and (3.51) into Eq. (3.41), we obtain
dη 2 dμ2
ds2 = (η 2 − μ2 )f (η, μ) 2
−
1−η 1 − μ2
wave path and then define the polarization angle relative to this parallelly
transported basis. In this way, we can see that the changes of polariza-
tion of the wave has an absolutely physical meaning and independent of
the choice of the coordinates. For example, if the changes (relative to the
parallelly transported basis) are zero, it means that the polarization vector
μ
E(2) defined in Section 2.3 is parallel to the parallelly transported basis.
Clearly, such defined parallelism is independent of the coordinates.
To find a parallelly transported basis along the null geodesics defined
μ μ
by lμ , we make a rotation in the (E(2) , E(3) )-plane as we did in Eq. (2.32)
μ μ(0)
but with the angle φ4 now denoted by φ4 , and E(2) and E(3) by λμ(2) and
λμ(3) , respectively. Then, we find
1 (0)
λμ(2);ν lν = (sinh W V,ν − 2ϕ4,ν )lν λμ(3) ,
2
(3.58)
1 (0)
λμ(3);ν lν = − (sinh W V,ν − 2ϕ4,ν )lν λμ(2) .
2
(0)
Therefore, if the angle φ4 is chosen so that
(0)
sinh W V,v − 2ϕ4,v = 0, (3.59)
the space-like orthogonal vectors and λμ(2) λμ(3)
are parallelly transported
along the null geodesics (or the Ψ4 -wave path) defined by lμ , and the
difference
(0)
θ4 ≡ ϕ4 − ϕ4 , (3.60)
defines the angle between the polarization direction of the Ψ4 -wave and the
λμ(2) direction.
μ μ
Similarly, if the basis (E(2) , E(2) ) is rotated so that the rotated angle
(0)
ϕ0 satisfies
(0)
sinh W V,u − 2ϕ0,u = 0, (3.61)
the vectors λμ(2) and λμ(3) are parallelly transported along the Ψ0 -wave path,
and the angle
(0)
θ0 ≡ ϕ0 − ϕ0 , (3.62)
determines the polarization direction of the Ψ4 -wave relative to the
(λμ(2) , λμ(3) )-basis. From Eqs. (3.59) and (3.61), we can see that the angle
(0) (0)
φ4 (φ0 ) is constant in Region III (Region II). Without loss of generality,
(0) (0)
we can choose φ4 (φ0 ) to be zero in Region III (Region II), such that the
μ μ μ μ
(λ(2) , λ(3) )-basis coincides with the (E(2) , E(3) )-basis in these regions.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch03 page 73
On the other hand, Eq. (3.24) shows that the Ψ0 - and Ψ4 -waves, in gen-
eral, consist of two parts. We call the δ-function part the impulsive part, and
the H-function part the shock part, although the latter still includes three
different cases: smooth wave, shock wave and the wave with unbounded
wavefront (Tsoubelis and Wang, 1989). The treatment for these three cases,
however, is the same for the present problem, so we shall not distinguish
them in the following discussions. We shall use θ0sh and θ0Im to denote the
polarization angles for the shock and impulsive parts of Ψ0 , and θ4sh and θ4Im
for the shock and impulsive parts of Ψ4 , respectively. Thus, in Region II we
have θ0sh = θ0sh (v) and θ0Im = constant. That is, in Region II the shock part
in general is non-collinearly polarized, while the impulsive part is always
collinearly polarized. Along the wave path, the two angles θ0sh and θ0Im do
not change in this region
sh Im
θ0,u = 0 = θ0,u . (3.63)
Of cause, in Regions I and III the angles θ0sh and θ0Im vanish, since Ψ0
vanishes there.
Similarly, for the Ψ4 -wave, we have
sh
θ4,u = 0 = θ4Im,u , (3.64)
in Region III, and θ4sh and θ4Im vanish in Regions I and II, since Ψ4 vanishes
in these two regions.
In the interaction region (Region IV), the situation is different. From
the Bianchi identities (3.34)–(3.37), we find that
1
AΨIV
0,u = {A[4(ln B),u + U,u − 2i sinh W V,u ]Ψ0
IV
2
−3B(cosh W V,v − iW,v )ΨIV IV
2,v − 2BΦ02,v
and
1
AΨIm
0,u = {A[4(ln B),u + U,u − 2i sinh W V,u ]ΨIm
0
2
−2BΦIV Im
02 − A(cosh W V,u − iW,u )Φ00 }, (u > 0, v = 0),
(3.66)
1
BΨIm
4,v = {B[4(ln A),v + U,v + 2i sinh W V,v ]ΨIm
4
2
−2AΦIV Im
02 − B(cosh W V,v − iW,v )Φ22 }, (u = 0, v > 0).
Since from now on only the “scale invariant” terms are used, we shall drop
all of the super-indices “0” from these quantities without causing any confu-
sion. Then, substituting Eq. (3.67) into Eqs. (3.65) and (3.66), respectively,
we find
1
ΨIV
0,u = {[U,u − 2i sinh W V,u ]ΨIV
0
2
−3(cosh W V,v − iW,v )ΨIV IV
2 − 2Φ02,v
and
1
ΨIm
0,u = {[U,u − 2i sinh W V,u ]ΨIm
0
2
−2ΦIV Im
02 − (cosh W V,u − iW,u )Φ00 }, (u > 0, v = 0),
(3.69)
1
ΨIm
4,v = {[U,v + 2i sinh W V,v ]ΨIm
4
2
−2ΦIV Im
02 − (cosh W V,v − iW,v )Φ22 }, (u > 0, v = 0).
Hence, we find that
sh 1
θ0,u =− {3[cosh W V,v Im(ΨIV IV IV IV
0 Ψ2 ) + W,v Re(Ψ0 Ψ2 )]
4(ΨIV IV
0 Ψ0 )
+ 2 Im(ΨIV IV IV IV
0 Φ20,v ) + (2M,v − U,v )Im(Ψ0 Φ20 )
+ 2ΦIV IV IV
11 [cosh W V,v Im(Ψ0 ) + W,v Re(Ψ0 )]
+ ΦIV IV IV
00 [cosh W V,u Im(Ψ0 ) + W,u Re(Ψ0 )]
+ 2 Im(ΨIV IV IV IV
4 Φ02,u ) + (2M,u − U,u )Im(Ψ4 Φ02 )
+ 2ΦIV IV IV
11 [cosh W V,v Im(Ψ4 ) − W,v Re(Ψ4 )]
+ ΦIV IV IV
22 [cosh W V,v Im(Ψ0 ) − W,v Re(Ψ4 )]
+ 2 Im(ΨIm IV
0 Φ20 )}, (u > 0, v = 0),
(3.71)
1
Im
θ4,v = {[cosh W V,v Im(ΨIm Im Im
4 ) − W,v Re(Ψ4 )]Φ22
4(ΨIm Im
4 Ψ4 )
+ 2 Im(ΨIm IV
4 Φ02 )}, (u > 0, v = 0).
Equations (3.70) and (3.71) show that, due to the interaction between
the two incoming plane gravitational waves and the interaction with the
matter fields, the polarization directions of the Ψ0 - and Ψ4 -waves got
changed relative to the parallelly transported basis along their wave paths.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch03 page 76
That is, in the collinear case the polarizations of colliding plane gravita-
tional waves do not change.
To further illustrate the properties of the polarization of colliding plane
gravitational waves, in the rest of this section we restrict ourselves to several
specific cases (with W = 0), which are interesting from the point of view
of physics. Because of the symmetry shared by the two plane gravitational
waves, it is sufficient to consider only one of them, say, the Ψ0 -wave. In
addition, since we are now working in Region IV, we do not make any
more specific statements about it in the following, and understand all the
following results valid only in this region (plus its two boundaries u = 0, v >
0 and u > 0, v = 0).
Φij = 0, Λ = 0. (3.73)
sh 3
θ0,u =− {cosh W V,v Im(ΨIV IV
0 Ψ2 )
4(ΨIV IV
0 Ψ0 )
+ W,v Re(ΨIV IV
0 Ψ2 )}. (3.74)
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch03 page 77
sh
Thus, if Ψ2 = 0, then we have θ0,u = 0. That is, if Ψ2 = 0, then the
polarization of a plane gravitational wave does not change. However, it was
shown that the “Coulomb” part Ψ2 appears necessarily in the interaction
region (Region IV), due to the nonlinearity of the Einstein field equations
(Szekeres, 1972). Thus, the change of polarizations of two colliding purely
gravitational waves is purely due the nonlinear interaction between the two
plane gravitational waves.
On the other hand, from Eq. (3.68) we find
Im
θ0,u |v=0 = 0. (3.75)
That is, the impulsive plane gravitational wave does not change its polar-
ization after collision, when no matter is present in the spacetime. From
Eq. (3.36) we can see that in the present case the ΨIm
0 -wave component
does not interact with other components.
Tμv = ε1 lμ lν + ε2 nμ nν , (3.76)
which is the superposition of two pure radiation fields given by Eq. (1.94),
where ε1 and ε2 are non-negative. Equation (3.76) represents a pair of
oppositely moving null dust clouds with the energy density ε1 and ε2 ,
respectively, each of which is separately conserved (Taub, 1988a; Tsoubelis
and Wang, 1991). The corresponding non-vanishing Ricci scalars are
given by
ε2 ε1
Φ00 = , Φ22 = . (3.77)
2B 2 2A2
On the other hand, Eq. (3.25) shows that, like Ψ0 and Ψ4 , the compo-
nents Φ00 and Φ22 , in general, consist of two parts: the H-function part and
the δ-function part. The latter is supported on the hypersurfaces u = 0 or
v = 0, and usually interpreted as impulsive shells of null dust. When atten-
tion is restricted to the inside of the interaction region, only the H-function
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch03 page 78
sh 1
θ0,u =− {3[cosh W V,v Im(ΨIV IV
0 Ψ2 )
4(ΨIV IV
0 Ψ0 )
+ W,v Re(ΨIV IV IV IV
0 Ψ2 )] + Φ00 [cosh W V,u Im(Ψ0 )
+ W,u Re(ΨIV
0 )]}. (3.78)
sh
Unlike in the vacuum case, now θ0,u can be different from zero even when
IV
Ψ2 = 0, because of the presence of the last term in the right-hand side
of Eq. (3.78), which represents the interaction between ΨIV IV
0 and Φ00 [see
Eq. (3.34)].
It was shown that, when null dust is present the collision of two plane
gravitational waves does not require the Coulomb field Ψ2 to appear neces-
sarily in the interaction region (Tsoubelis and Wang, 1990). Thus, a plane
gravitational shock wave can change its polarization due to the deflection
effect. On the other hand, Eq. (3.71) now becomes
Im 1
θ0,u =− {[cosh W V,u Im(ΨIm
0 )
4(ΨIm Im
0 Ψ0 )
+ W,u Re(ΨIm Im
0 )]Φ00 }. (3.79)
Obviously, when ΦIm00 is different from zero, the polarization of the impulsive
part of Ψ0 changes after collision because of the interaction between ΨIm 0
and ΦIm
00 [see Eq. (3.36)].
sh 1
θ0,u =− {6[cosh W V,v Im(ΨIV IV
0 Ψ2 )
8(ΨIV IV
0 Ψ0 )
+ W,v Re(ΨIV IV
0 Ψ2 )]
which shows that a plane gravitational shock wave changes its polarization
due to both the nonlinear interaction between the two plane gravitational
waves and the interaction with the massless scalar field φ.
When the spacetime is filled only with a massless scalar field, we have
ΦIm Im
00 = 0 = Φ22 , (3.83)
or equivalently
which means that an impulsive plane gravitational wave does not change
its polarization when it passes through a massless scalar field, since in this
case there is no interaction between the gravitational impulsive wave and
the massless scalar field [see Eqs. (3.36) and (3.37)].
Note that Eq. (3.84) is also the condition under which the hypersurfaces
u = 0 and v = 0 are free of matter (Tsoubelis and Wang, 1989, 1990).
quantities
(0) (0)
Φ0 = BΦ0 , Φ2 = AΦ2 , (3.86)
and drop the superscript “0”, then the non-vanishing Ricci scalars are
ΦIV
00 = Φ0 Φ0 , ΦIV IV
02 = Φ0 Φ2 = Φ20 , ΦIV
22 = Φ2 Φ2 , (3.87)
and the Maxwell field equations (1.99a)–(1.99d) reduce to
2Φ0,u = (U,u − i sinh W V,u )Φ0 − (cosh W V,v − iW,v )Φ2 ,
(3.88)
2Φ2,u = (U,v + i sinh W V,v )Φ2 − (cosh W V,u + iW,u )Φ0 .
Note that in the present case all of the Ricci scalars have only the shock
part, otherwise, the Maxwell potentials Φi (or equivalently, the electromag-
netic field tensor Fμν ) will contain the square roots of δ-function, which is
not acceptable physically. Then, from Eqs. (3.70), (3.87) and (3.88) we find
sh 1
θ0,u =− {3[cosh W V,v Im(ΨIV IV
0 Ψ2 )
IV IV
4(Ψ0 Ψ0 )
+ W,v Re(ΨIV IV IV
0 Ψ2 )] + 2 Im(Ψ0 Φ0,v Φ2 )
− sinh W V,v Re(ΨIV IV
0 Φ0 Φ2 ) + 2M,v Im(Ψ0 Φ0 Φ2 )}. (3.89)
Thus, similar to the last two cases, a plane gravitational shock wave can
change its polarization when it interacts with an electromagnetic field.
On the other hand, for the impulsive part of Ψ0 , Eq. (3.71) becomes
Im Im(ΨIm
0 Φ0 Φ2 )
θ0,u =− . (3.90)
2(ΨIm Im
0 Ψ0 )
It follows that a gravitational impulsive wave can also change its polariza-
tion due to the interaction between the impulsive wave and the electromag-
netic field component ΦIV 02 .
to the case for a pure radiation field, which we have already discussed in
Subsection 3.4.2.
2 A rigorous proof of the instability for the case in which the initial waves have con-
stant aligned polarizations had to wait until 2005, when Griffiths first showed that such
horizons are indeed unstable with respect to bounded perturbations of the initial waves
(Griffiths, 2005). However, such a proof is still absent for the non-collinear case.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch03 page 82
p1 (z) = 0. (3.93)
On the other hand, Tsoubelis and Wang studied the effect of polariza-
tions of colliding plane gravitational waves on the formation and nature
of singularities (Tsoubelis and Wang, 1992). Specifically, it was found that
some astigmatic singularities are turned into anastigmatic singularities, or
vice versa, due to the interaction between different polarization modes of
colliding plane gravitational waves. Moreover, it was also found that all
of the solutions that are free of singularities in the collinear case are so
in the non-collinear case, but inversely not. That is, due to the interac-
tion between different polarization modes some singularities are turned into
Cauchy horizons.
In addition, Chandrasekhar and Xanthopoulos (1987a) studied the
effects of sources on the Cauchy horizons, and found that the coupling
of gravitational waves with an electromagnetic field does not affect in any
essential way the development of the Cauchy horizons, if the polarizations
of the colliding gravitational waves are not parallel. However, if the polar-
izations are parallel, the space-like singularity that occurs in the vacuum is
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch03 page 83
and
1 1
(ln f ),t = − + Tr(A2 + B 2 ),
t 4t
(3.95)
1
(ln f ),z = Tr(AB),
2t
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch03 page 86
N
N
μk (0) −1 (j) (k)
gab = gab − Γjk Na Nb , (a, b = 2, 3), (3.103)
t
k=1 j,k=1
where
3
(j) (0) (k)
na gab nb
Γjk ≡ ,
μj μk − t2
a,b=2
3
(k) (0)
n g
Na(j) ≡ b ab
(without summation for the index k), (3.104)
μk
b=2
3
(k) (k)
n(j)
a ≡ mb Mab (without summation for the index k),
b=2
ma (k+N/2) = m̄(k)
a . (3.105)
It must be noted that the term in the square bracket in the first equation
of Eq. (3.103) should be replaced by unit when N = 1.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch03 page 88
From the above general description of ISM (BZ), we can see that a com-
plete solution is obtained only after the system (3.98) is explicitly solved
for a given seed g (0) . However, as noticed by various authors, this in gen-
eral is not trivial. The solutions for Ψ(0) were found for several simple
diagonal seeds (Jantzen, 1980) and for a few non-diagonal seeds (Belinsky
and Francaviglia, 1982; Letelier, 1986, 1989; Kitchingham, 1986; Cespedes
and Verdaguer, 1987). Therefore, some further simplifications are usually
adopted. A common assumption is that the matrices Ψ(0) and g (0) are diag-
onal. With the latter assumption, Economou and Tsoubelis (1989; see also
Letelier, 1986) found that the corresponding expressions of Eqs. (3.103)
and (3.104) can be written in a very simple form in terms of the determi-
nants of four N × N matrices. The Economou–Tsoubelis results come from
the following observations. Integrating the trace of Eq. (3.98) along a pole
trajectory μk , we have
det(ψ (0) ) = 2wk μk (without summation for the index k), (3.106)
t (0)
(ln ψk ),ζ = (ln g22 ),ζ , (3.108)
t − μk
t (0)
(ln ψk ),ξ = (ln g22 ),ξ , (3.109)
t + μk
with
ζ ≡ t + z, ξ = t − z. (3.110)
(0)
μk Qk g22
σk ≡ , sk ≡ σk (without summation for the index k),
t ψk2
sj sk+1
Δjk ≡ (without summation for the indexes j and k),
σj σk − 1
(3.111)
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch03 page 89
where Qk are arbitrary constants, Economou and Tsoubelis found that the
new solutions can be written in the form
N
(0)
g22 = |σk | L(−1) L−1
(0) g22 ,
k=1
N
g23 = t |σk | LL−1 (0),
k=1
N
(0)
g33 = |σk |−1 |L(+1) L−1 (0)g22 , (3.112)
k=1
⎡ ⎤
(0) 2
N
N
f = Cph f (0) L(0) t−N /2
|σk |N ⎣ (σk − σ1 )−2 ⎦
k=1 l,k=1; k>l
N
× |sk |−1 ,
k=1
where
(σ σ )d s s +1
L(d) ≡ det Δ(d) , Δ(d)jk ≡ j σkj σk −1j k
, (d = 0, ±1),
(3.113)
s
L ≡ det Δ(0) − det Δ(0)jk + σj σj k .
Later, Tsoubelis and Wang found (Wang, 1991a; Tsoubelis and Wang,
1992) that the integration of Eq. (3.108) will become simpler, if we work
in the coordinates η and μ introduced in Section 3.3. Actually, introducing
the quantity Σ(0) via the relation
ψ (0) (λ, t, z) = (2wλ)1/2 diag{Σ(0) , Σ(0)−1 }, (3.114)
it was found that the system (3.98) (not restricted only along the pole
trajectories) reduces to
[(1 − η 2 )∂η + λ∂μ − 2λη∂λ ] ln Σ(0) = 1 − η 2 V (0) ,η ,
(3.115)
[(1 − μ2 )∂μ + λ∂η − 2λμ∂λ ] ln Σ(0) = (1 − μ2 )V (0) ,μ ,
and the “initial condition” of Eq. (3.115) becomes
(0)
Σ(0) (λ = 0, t, z) = eV . (3.116)
The integrability conditions for Eq. (3.115) is the Einstein vacuum field
equation for the seed function V (0) ,
[(1 − η 2 )V (0) ,η ],η − [(1 − μ2 )V (0) ,μ ], μ = 0, (W (0) = 0). (3.117)
b2530 International Strategic Relations and China’s National Security: World at the Crossroads
Chapter 4
91
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch04 page 92
scalars zero
Φij = 0, Λ = 0. (4.1)
The conditions for which the two null hypersurfaces u = 0 and v = 0 are
free of matter are
U,u (u = 0, v) = 0 = U,v (u, v = 0). (4.2)
For the special gauge of Eq. (3.40), these conditions together with Eq. (3.38)
imply
1
n, m > . (4.3)
2
In terms of η and μ, the Einstein vacuum equations (Chandrasekhar
and Ferrari, 1984) can be written in the form,
Re(Z){[ 1 − η 2 Z,η ],η − [ 1 − μ2 Z,μ ],μ }
2 2
= (1 − η 2 ) (Z,η ) − 1 − μ2 (Z,μ ) , (4.4)
and
μ η 2
2
Γ,η + 2
Γ,μ = − Z,η Z ,μ + Z ,η Z,μ , (4.5a)
1−μ 1−η (Z + Z) 2
3 1
2ηΓη + 2μΓ,μ = 2
+
1−η 1 − μ2
4 2 2
− 2
[ 1 − η 2 |Z,η | + 1 − μ2 |Z,μ | ],
(Z + Z)
(4.5b)
where
1/2
1 − η2
f≡ eΓ . (4.6)
η 2 − μ2
Equation (4.4) is equivalent to the Ernst equation (3.43) with the Ernst
potential E defined by Eq. (3.44). In terms of the Ernst potential E,
Eqs. (4.5a)–(4.5b) can be written as
μ η 2
2
Γ,η + 2
Γ,μ = − 2 E,η E,μ +E,η E,μ , (4.7a)
1−μ 1−η 2
(1 − |E| )
3 1
2ηΓ,η + 2μΓ,μ = 2
+
1−η 1 − μ2
4 2 2
− 2 [ 1 − η 2 |E,η | + 1 − μ2 |E,μ | ].
(1 − |E|2 )
(4.7b)
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch04 page 93
It can be shown that the above equations are equivalent to the ones
given by Eqs. (3.94)–(3.97).
On the other hand, Chandrasekhar and Ferrari (1984) showed that, as
we did in the stationary axisymmetric case, if we introduce a potential Q
for the function q2 in the manner,
1 − μ2 1 − η2
Q,η = q2,μ , Q,μ = q2,η , (4.8)
χ2 χ2
then the complex function Z + defined by
Z + ≡ P + iQ (4.9)
Then, the Einstein vacuum equations of Eqs. (4.4), (4.5a) and (4.5a)
reduce to
and
μ η χ,η χ,μ
2
Γ,η + 2
Γ,μ = − , (4.13a)
1−μ 1−η χ2
3 1 1 − η 2 χ2,η + 1 − μ2 χ2,μ
2ηΓ,η + 2μΓ,μ = + − . (4.13b)
1 − η2 1 − μ2 χ2
In terms of V [see Eq. (3.44)], Eq. (4.12) takes the form
1 − η 2 V,η ,η − 1 − μ2 V,μ ,μ = 0. (4.14)
V = V1 + V2 , (4.17)
such an integration has not been completely solved yet, and only for some
particular choice of V , the corresponding solutions of Γ have been found
(Griffiths, 1987; Tsoubelis and Wang, 1989; 1992; Li, 1989; Wang, 1991a).
In the following, let us consider the solutions (Tsoubelis and Wang,
1989),
V = a ln 1 − η 2 1 − μ2 − 2δ1 Q0 (η)P0 (μ) − 2δ2 P0 (η)Q0 (μ)
1−η 1−μ
= a ln 1 − η 2 1 − μ2 + δ1 ln + δ2 ln , (4.18)
1+η 1+μ
where δ1 and δ2 are two constants related to the parameters n and m
introduced in Section 3.3.
In order to complete the integration of Eqs. (4.13a) and (4.13b), we
introduce the quantity Σ by
η 2 − μ2
Γ = ln + Σ. (4.19)
(1 − η 2 )3/4 (1 − μ2 )1/4
Substituting the above expression into Eqs. (4.13a) and (4.13b), we obtain
μ η
Σ,η + Σ,μ = −V,η V,μ , (4.20a)
1 − μ2 1 − η2
2ηΣ,η + 2μΣ,μ = − 1 − η 2 (V,η )2 − 1 − μ2 (V,μ )2 . (4.20b)
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch04 page 95
It can be shown that for the function V given by Eq. (4.18), the system
of Eqs. (4.20a) and (4.20b) has the following solution (Tsoubelis and Wang,
1989):
2 2
Σ = −δ+ ln(η + μ) − δ− ln(η − μ) + (a + δ1 )2 ln(1 − η)
+ (a − δ1 )2 ln(1 + η) + (a + δ2 )2 ln(1 − μ)
+ (a − δ2 )2 ln(1 + μ) + ln C0 , (4.21)
where C0 is an arbitrary constant, and δ± are defined by
δ± ≡ δ1 ± δ2 . (4.22)
To summarize the above results, we obtain the following solutions:
C0
f = 2 2 (1 − η)
b1
(1 − μ)b2 (1 + η)c1 (1 + η)c2 ,
(η + μ) (η − μ)δ−
δ+
1/2
e−U = 1 − η 2 1 − μ2 ,
(4.23)
a 1−η
δ1
1−μ
δ2
e = 1 − η 2 1 − μ2
V
,
1+η 1+μ
W = 0,
where
1 1
bA = (a + δA )2 − , cA = (a − δA )2 − , (A = 1, 2). (4.24)
4 4
The solutions of Eq. (4.23) include most of the known diagonal (W = 0)
solutions. For example, when a = 0, we obtain Szekeres’ family of the
colliding plane gravitational wave solutions (Szekeres, 1972), which includes
the Szekeres solution (Szekeres, 1970), in which we have n = m = 2, and the
Khan–Penrose solution (Khan and Penrose, 1971), in which n = m = 1,
where n and m are defined below in Eq. (4.28). For any given a, when
δ1 = 1 and δ2 = 0, we obtain the Ferrari–Ibañez solutions (Ferrari and
Ibañez, 1987a).
In terms of t and z, the function V given by Eq. (4.23) can be written as
1+z+ (−1 − z)2 − t2
V = 2a ln t − δ+ ln
t
1−z+ (1 − z)2 − t2
− δ− ln . (4.25)
t
The last two terms in the right-hand side of Eq. (4.25) correspond to the
soliton structure, while the first term can be thought of as producing a
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch04 page 96
Fig. 4.1. The soliton structure for the solutions given by Eq. (4.25).
and Eq. (3.49) further restricts the region of validity of the solutions to that
between the two light-cones [Region IV in Fig. 4.1].
Thus, the Khan–Penrose substitutions of Eq. (3.11) actually extend
these solutions defined in the triangle ABC beyond the lines AC and BC
to Regions I–III, as shown in Fig. 4.1.
On the other hand, combining Eqs. (3.42) and (4.23) we find
In writing Eq. (4.27) we have set 8nm2−(δ+ +δ− ) C0 = 1 by using the arbi-
trariness of C0 . Equation (4.27) shows that the application of the Khan–
Penrose substitutions Eq. (3.11) will result in singularity at u = 0 or v = 0,
unless the parameters n and m are chosen so that
1 1
n= 2 , m= 2 . (4.28)
2 − δ+ 2 − δ−
Equations (4.3) and (4.28) are the restrictions on the choice of the soliton
parameters δ1 and δ2 . The solutions which do not satisfy these conditions
are not acceptable for the collision of purely plane gravitational waves.
However, further considerations from the point of view of physics (Tsoubelis
and Wang, 1989) show that the conditions of Eq. (4.3) are too weak and
that for 1/2 < m, n < 1 the extension obtained above cannot guarantee the
integrability of the Weyl scalars corresponding to the waves under collision
(Szekeres, 1972), as it will become clear in the following. Therefore, we
impose the conditions
n, m ≥ 1, (4.29)
M = U = V = 0, (4.30)
which is flat.
In Region II, where u < 0, and 0 < v < 1, we have
b1 +c2 +(1−n)/2n−1/2
e−M = (1 + v m )b2 +c1 +(1−n)/2n−1/2 (1 − v m ) ,
e−U = 1 − v 2m ,
2a+δ−
(4.32)
eV = (1 + v m )2a−δ− (1 − v m ) ,
W = 0.
e−U = 1 − u2n ,
2a−δ+ 2a+δ+
eV = (1 + un ) (1 − un ) , (4.33)
W = 0.
Substituting Eq. (4.32) into Eqs. (3.7a)–(3.7c), we find that the only
non-vanishing Weyl scalar is Ψ0 , given by
ΨII
0 (v) = −e
2U
{2am2 1 − 4a2 v 4m−2 − 12a2 m2 δ− v 3m−2
− 6am(2m − 1)v 2m−2 − m(m − 1)δ− v m−2 }. (4.34)
Similarly, ψ4III is the only non-vanishing Weyl scalar in Region III, and is
given by
ΨIII
4 (u) = −e
2U
{2an2 (1 − 4a2 )u4n−2 − 12a2 n2 δ+ u3n−2
− 6an(2n − 1)u2n−2 − n(n − 1)δ+ un−2 }. (4.35)
In Region IV, all three of the Weyl scalars Ψ0 , Ψ2 and Ψ4 do not vanish,
but the corresponding expressions are very complicated, so we shall not
present them explicitly here. Anyhow, as far as the behavior of the Weyl
scalars across the hypersurfaces u = 0 and v = 0 is concerned, we only need
to know their limits as u → 0+ and v → 0+ , which are given as follows.
When u → 0+ , we find
ΨIV II
0 (0, v) → Ψ0 (v),
nmδ+ n−1 m−1
ΨIV
2 (0, v) → u v (2av m + δ− ),
X3
(4.36)
n n−2
ΨIV
4 (0, v) → u [(n − 1)δ+ X + 6a(2n − 1)un
X2
+ 3(2n − 1)δ− v m un ],
and when v → 0+ , we have
m m−2
ΨIV
0 (u, 0) → v [(m − 1)δ− Y + 6a(2m − 1)v m
Y2
+ 3(2m − 1)δ+ v m un ],
(4.37)
nmδ− n−1 m−1
ΨIV
2 (u, 0) → u v (2aun + δ+ ),
Y3
ΨIV II
4 (u, 0) → Ψ4 (u).
Combining the above results with Eqs. (3.24) and (3.26), it can be seen
that the Weyl scalars behave as follows when we cross the null hypersurface
separating region A from region B (A → B):
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch04 page 99
(i) IV → II:
Ψ0 is continuous,
Ψ2 = ΨIV
2 H(u), (4.38)
nδ+ n−1
Ψ4 = ΨIV
4 H(u) + u δ(u),
X
(ii) IV → III:
mδ− m−1
Ψ0 = ΨIV
0 H(v) + v δ(v),
Y
Ψ2 = ΨIV
2 H(v), (4.39)
ψ4 is continuous,
(iii) II → I:
Ψ0 = ΨII
0 H(v) + mδ− v
m−1
δ(v),
(4.40)
Ψ2 , Ψ4 are continuous,
(iv) III → I:
Ψ2 , Ψ4 are continuous,
(4.41)
Ψ4 = ΨIII
4 H(u) + nδ+ u
n−1
δ(u).
Equations (4.34) and (4.40) show clearly that the parameter m determines
the type of the wave incident in Region II. Particularly, by observing the
behavior of Ψ0 across the hypersurface v = 0 in the direction II → I, we
can distinguish the following cases.
(a) m = 1. In this case, we have
Ψ0 = 6aH(v) + δ− δ(v). (4.42)
Therefore, when a = 0, the gravitational wave incident in Region II has the
form of an impulsive + shock wave. When a = 0, only the impulsive part
remains.
(b) 1 < m < 2. Then, we have
Ψ0 = m(m − 1)δ− v m−2 H(v), (4.43)
which means that the incoming gravitational wave has an unbounded wave-
front of the form v γ , with γ ∈ (0, 1).
(c) m = 2. Then, we have
Ψ0 = 2δ− H(v), (4.44)
which corresponds to a shock wave.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch04 page 100
1 2
× − (a − δ2 ) − 4au2n (δ2 + δ1 v 2m ) ,
4
3mu1−2nv 2m−1 IV
ΨIV
0 → − Ψ2 , (4.48)
2n(a + δ1 )
3nu2n−1v 1−2m IV
ΨIV
4 → − Ψ2 ,
2m(a + δ1 )
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch04 page 101
as t → 0+ , where
2
+3/4
τ = tα , α = a + δ1 , (4.51)
and
α2 − 14 α − 12 α + 12
p1 ≡ , p2 ≡ − , p3 ≡ , (4.52)
α2 + 34 α2 + 34 α2 + 34
(0)
and εμ are functions of z only. It is easy to show that the pk defined
by Eq. (4.52) satisfy the Kasner relations Eq. (3.92), but now the pk
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch04 page 102
with
b A ≡ (a + δA ) (a + δA − 1) , c A ≡ (a − δA ) (a − δA − 1), (4.58)
where q ≡ 4δ1 −a Q1 is another arbitrary constant.
Note that when η → μ the functions A and B defined by Eq. (4.57)
become unbounded unless δ− ≥ 0. Therefore, for the extension of the above
solutions beyond the hypersurfaces η = ±μ (equivalent to the u = 0 and
v = 0 hypersurfaces), we restrict the solutions of Eqs. (4.56)–(4.58) only to
the cases in which δ− ≥ 0. For δ− < 0, we write these solutions in the form
A (0)
χ = [(1 − η 2 )(1 − η 2 )]1/2 χ ,
B
(η − μ)−2δ− +1 (0) 2
q2 = −2q [χ ] , (4.59)
B
(1) (1 − η)b1 (1 − μ)b2 (1 + η)c1 (1 + μ)c2
f = Cph 2 A ,
(η + μ)δ+ (η − μ)(δ− −1)2
where
Inserting the Khan–Penrose substitutions Eq. (3.11) into the above solu-
tions and considering Eq. (3.42), we find that the extended solutions in the
pre-collision regions are given as follows.
In Region I, the metric takes the form of Eq. (4.31), which means that
the spacetime is flat.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch04 page 105
But, after replacing δ1 and δ2 by δ1 and δ2 , respectively, these are the
seed solutions given by Eqs. (4.23) and (4.24). Thus, the asymptotic behav-
ior of the solutions given by Eq. (4.69) is described by Eqs. (4.50)–(4.52) but
with α being replaced by α ≡ a + δ1 = α + 1/2. Specifically, corresponding
to Eq. (4.53), we now have
α = 0, −1. (4.71)
In the non-collinear case, the solutions given by Eqs. (4.56)–(4.58) have
the following limit as t → 0+ .
(a) α < −1/2. In this case, it can be shown that the above solutions have the
same limit as the corresponding diagonal solutions as t → 0+ . Specifically,
the solutions with α = −1 are free of spacetime curvature singularities on
the hypersurface t = 0 in both of the collinear and non-collinear cases.
(b) − 1/2 ≤ α < 0. In this case, the solutions have the limit
ds2 = ε20 dτ 2 − ε21 τ 2p1 dz 2 − ε22 τ 2p2 (dX(z)
2 2
) − ε23 τ 2p3 (dX(z)
3 2
) , (4.72)
as t → 0+ , where τ and pk are given by Eqs. (4.51) and (4.52) with α
being replaced by α = α + 1/2. Here ε2i ’s are functions of z only, and
2 3
X(z) and X(z) are z-dependent linear combinations of x2 and x3 at each
fixed point of z. Equation (4.72) shows that, relative to the corresponding
diagonal solutions, the Kasner exponents pk are not changed. Thus, in the
both collinear and non-collinear cases, the nature of singularities on the
hypersurface t = 0 is the same, and is anastigmatic. However, the Kasner
axes
2 along
which the exponents p2 and p3 are defined are rotated from
x , x3 to (X(z)
2 3
, X(z) ) at each fixed point z.
(c) 0 ≤ α ≤ 1/2. Then, the metric takes the form
ds2 = ε30 dτ 2 − ε31 τ 2p1 dz 2 − ε32 τ 2p2 (dX(z)
2 2
) − ε3 3τ 2p3 (dX(z)
3 2
) , (4.73)
as t → 0+ , but now τ and pk are defined by
2
−α+1 α(α − 1)
τ ≡ tα , p1 ≡ ,
α2 − α + 1
(4.74)
α 1−α
p2 ≡ 2
, p3 ≡ 2 .
α −α+1 α −α+1
The comparison of the above limit with the one in the collinear case shows
that the nature of the singularities on the hypersurface t = 0 is changed
in the present case. In the non-collinear case it is anastigmatic, while in
the collinear case it is astigmatic. Moreover, the Kasner axes in the present
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch04 page 108
case are also rotated. When α = 0, Eq. (4.74) gives p1 = 0. That is, the
solutions with α = 0 are free of curvature singularities on the hypersurface
t = 0.
(d) α > 1/2. Then, we find
2 2
ds2 = ε40 dτ 2 − ε41 τ 2p1 dz 2 − ε42 τ 2p2 dx2 − ε43 τ 2p3 dx3 , (4.75)
where τ and pk are given by Eqs. (4.74) and (4.75). Equation (4.75) shows
that in the present case the Kasner axes remain the same as these in the
collinear case, but the Kasner exponents are changed. When 1/2 < α < 1,
we have p1 < 0. Thus, the nature of the singularities is anastigmatic (note
that in the collinear case, it is astigmatic). When α > 1, we have p1 > 0,
which means that the nature of the singularities is astigmatic, as it is in
the collinear case. When α = 1, we have p1 = 0. Thus, the solutions with
α = 1 are free of spacetime curvature singularities on the hypersurface
t = 0, too. The above analysis shows that the non-diagonal solutions given
by Eqs. (4.56)–(4.58) are free of spacetime curvature singularities on the
hypersurface t = 0, if any one of the following conditions holds:
(i) a + δ1 = −1, (ii) a + δ1 = 0, or (iii) a + δ1 = 1. (4.76)
Comparing Eq. (4.71) with Eq. (4.76), we find that the solutions with α = 0
or −1 are free of curvature singularities in both of the collinear and non-
collinear cases, whereas the solutions with α = 1 are free of spacetime
curvature singularities only in the non-collinear case. The above difference
should be obviously attributed to the presence of the ×-polarization mode
of the plane gravitational wave moving toward the right-hand side. In other
words, the interaction between different polarization modes can change the
nature of the singularities on the hypersurface t = 0 and turn some space-
like singularities into Cauchy horizons.
Ψ2 (u, v) = ΨIV
2 (u, v)H(u)H(v),
where
ΨII IV
0 (v) = Ψ0 (0, v) = 0, ΨIII IV
4 (u) = Ψ4 (u, 0) = 0,
(4.81)
D(v) ≡ [(1 + q 2 )2 (1 − v 2 H(v))2 + 16q 2 v 2 H(v)]1/2 .
Equations (4.80) and (4.81) show that this model represents the collision
of two impulsive plane gravitational waves. Thus, the corresponding solu-
tion is either the Khan–Penrose solution (Khan and Penrose, 1971) or the
Nutku–Halil solution (Nutku and Halil, 1977). A detained analysis of the
polarizations of these two plane gravitational impulsive waves shows that
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch04 page 110
it is the Nutku–Halil solution, since the polarization angles for these two
incoming impulsive waves are different, although they are constants.
This solution was studied in detail by Chandrasekhar and Ferrari but
with a different form (Chandrasekhar and Ferrari, 1984). It was found that
the solution is singular on the hypersurface t = 0.
Case A.2: a = 1/2. In this case, the non-vanishing Weyl scalars are
given by
1
Ψ0 (u, v) = ΨIV
0 (uH(u), v)H(v) + [(1 − q 2 ) − i2q]δ(v),
Y (1 + q 2 )
1
Ψ4 (u, v) = ΨIV
4 (u, vH(v))H(u) + {[(1 − vH(v))4
CDX
+ 2q 2 [1 + 2v 2 H(v) − 3v 4 H(v)] + q 4 (1 + vH(v))4 ] (4.82)
+ i4qv 2 H(v)[(1 − v)2 − q 2 (1 + v)2 ]}δ(u),
Ψ2 (u, v) = ΨIV
2 (u, v)H(u)H(v),
where
ΨII IV
0 (v) = Ψ0 (0, v)
3 1 − v2
=− 1/2
{[16q 4 − (1 − q 2 )4 + (9 − 14q 2 )v]
(1 + q 2 ) C 2 F
− i2q(1 + q 2 )[(1 + q 2 )(5 + v 2 )v − 2(1 − q 2 )(1 + 2v 2 )]
− [(1 + q 2 )(1 + 6v 2 + v 4 ) − 4(1 − q 2 )(1 + v 2 )v]},
(4.83)
3(1 + u)
ΨIII
4 (u) = ΨN
4 (u0) = ,
(1 − u2 )2
C ≡ (1 − vH(v))2 + q 2 (1 + vH(v))2 ,
D ≡ [C 2 + 16q 2 v 2 H(v)]1/2 ,
F ≡ [(1 + q 2 )(1 + 62 + v 4 ) − 4(1 − q 2 )(1 + v 2 )v]3/2 ,
and
(1 + η)2 (1 + μ)2
ΨIV
2 (u, v) = {(1 − η)3 − q 2 (1 − μ)[2 − (1 − η)(η + μ)
XY A2
+ 2η 2 μ − q 2 (1 − 3μ2 + 2μ3 )]
iq(1 + η)(1 + μ)
+
B[(1 − η)3 (1 + μ) + q 2 (1 − μ)3 (1 + η)]
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch04 page 111
and
ΨII IV
0 (v) = Ψ0 (0, v) = 0, ΨIII IV
4 (u) = Ψ4 (u, 0) = 4kY
−4
. (4.91)
Equations (4.90) and (4.91) show that in this case the solution represents
the collision of a constantly polarized shock wave and an impulsive wave.
Comparing it with Case A.1, we find that one of the two incoming impul-
sive gravitational waves in the Nutku–Halil solution now is replaced by a
constantly polarized shock wave.
To study the singularity behavior of the solution on the hypersurface
t = 0, we find that the Weyl scalar ΨIV 2 is given by
4ku
ΨIV
2 (u, v) = {B[(1 − η 2 )4k [2 − (1 − k)η 2 − (1 + k)μ2 ]
XY t2 A2 B
− 6kq 2 (η 2 − μ2 )t4k
− q 4 (1 − μ2 )4k [2 − (1 − k)μ2 − (1 + k)η 2 ]]
− iqt2k (1 − η 2 )4k [k −1 (η 2 − μ2 )(1 + η)(1 − μ)
+ (2 + 3η + 6η 2 + 3η 3 ) − μ(3 + 7η + η 2 − η 3 )
− μ2 (1 + 2η + 5η 2 ) + 2μ3 (1 + η)]}. (4.92)
Obviously, as t → 0+ , the Weyl scalar ΨIV2 becomes unbounded. That is,
similar to the Nutku–Halil solution (Chandrasekhar and Ferrari, 1984), a
spacetime singularity is finally developed on the hypersurface t = 0 due to
the mutual focus of the two colliding gravitational waves.
Case B.2: a + δ1 = 1. In this case, the non-vanishing Weyl scalars are
given by
1
Ψ0 (u, v) = ΨIV
0 (uH(u), v)H(v) + [(1 − q 2 ) − i2q]δ(v),
Y (1 + q 2 )
Ψ4 (u, v) = ΨIV
4 (u, vH(v))H(u), (4.93)
Ψ2 (u, v) = ΨIV
2 (u, v)H(u)H(v),
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch04 page 113
and
1
ΨII IV
0 (v) = Ψ0 (0, v) = {[O(v) − q 2 P (v) − q 4 P (−v)
2 (1 − v 2 ) A2 D
+ q 6 O(−v)] − i2q(1 − v 2 )1+2k [Q(v) − 2q 2 Z(v) − q 4 Q(−v)]},
ΨIII IV
4 (u) = Ψ4 (u, 0) = 4Y
−4
[k + 9au3 + 24a2 ku5 + 2a(4a2 − 1)u6 ],
(4.94)
where
O(x) ≡ (1 + x)12k (1 − x)12 3(4 + 7x + 2x2 )
Thus, unlike the last case, the Nutku–Halil impulsive wave incident in
Region II now is replaced by a variably polarized shock + impulsive wave.
The solution in the present case is free of spacetime singularity on the
hypersurface t = 0 only when q = 0, and otherwise it is always singular.
This can be seen from the Weyl scalar ΨIV2 , which is given by
4u
ΨIV
2 (u, v) = {B[J(η, μ) − 6q 2 K(η, μ) − q 4 J(μ, η)]
XY BA2 (1 + η)(1 + μ)
− i2q(1 + η)2k+1 (1 + μ)2k+1 [L(η, μ) + 2q 2 N (η, μ)
+ q 4 L(μ, η)]}, (4.96)
where
which means that the polarization of the impulsive plane gravitational wave
with support on the hypersurface v = 0 does not change even after the
collision.
Substituting Eqs. (4.43) and (4.78) into Eq. (3.59), and then integrating
it, for v > 0 we find,
2
Im(0) 1 −1 1 + q v − 1 − q2
ϕ4 = tan
2 2q
2 2
1 + q v + 1 − q
− tan−1 + C0 . (4.102)
2q
Im(0)
From the condition that when v → 0, ϕ4 = 0, we have
4q 1 − q 2
tan C0 = 2. (4.103)
8 − (3 − q 2 )
Thus, Eqs. (4.102) and (4.103) yield,
Im(0) 1 −1 4q 1 − q 2 v 2
ϕ4 = tan H(v). (4.104)
2 (1 + q 4 ) (1 − v 2 ) + 2q 2 (1 + 3v 2 )
where C0 is a constant,
δ1 δ2
1−η 1−μ
χ(0) = [(1 − η 2 )(1 − μ2 )]a , (4.111)
1+η 1+μ
and
A ≡ (1 − η 2 )(1 − μ)2(a+δ2 ) (1 + μ)2(a−δ2 ) [(1 − η)2(a+δ1 )
+ Q1 Q2 (η − μ)2δ − (η + μ)2δ + (1 + η)2(a−δ1 ) ]2
+ (1 − μ2 )(1 − η)2(a+δ1 ) (1 + η)2(a−δ1 ) [Q1 (η + μ)2δ + (1 + μ)2(a−δ2 )
− Q2 (η − μ)2δ − (1 − μ)2(a+δ2 ) ]2 ,
B ≡ (1 − μ)2(a+δ2 ) (1 + μ)2(a−δ2 ) [(1 + η)(1 − η)2(α+δ1 )
− Q1 Q2 (1 − η)(η − μ)2δ − (η + μ)2δ+ (1 + η)2(a−δ1 ) ]2
(4.112)
+ (1 − η)2(a+δ1 ) (1 + η)2(a−δ1 ) [Q1 (1 − μ)(η + μ)2δ + (1 + μ)2(a−δ2 )
+ Q2 (1 + μ)(η − μ)2δ− (1 − μ)2(a+δ2 ) ]2 ,
C ≡ −2{Q1 (η + μ)2δ+ +1 [(1 + μ)2(a−δ2 ) (1 − η)2(a+δ1 )
+ Q22 (η − μ)4δ− (1 − μ)2(a+δ2 ) (1 + η)2(a−δ1 ) ]
+ Q2 (η − μ)2δ− +1 [(1 − μ)2(a+δ2 ) (1 − η)2(a+δ1 )
+ Q21 (η + μ)4δ+ (1 + μ)2(α−δ2 ) (1 + η)2(a−δ1 ) ]},
where Q1 and Q2 are arbitrary constants, and
1 1
ρA ≡ (a + δA )(a + δA − 2) − , rA ≡ (a − δA )(a − δA − 2) − ,
4 4
(A = 1, 2). (4.113)
Note that when η → ±μ the functions A, B and C defined by Eq. (4.112)
become unbounded unless δ± ≥ 0. Therefore, for the extension of these solu-
tions beyond the hypersurfaces η = ±μ, we restrict the solutions given in
the form of Eqs. (4.110)–(4.111) only to the cases where δ± ≥ 0. Otherwise,
we write these solutions in the following different forms.
(i) δ+ > 0, δ− < 0: In this case, the solutions given by Eqs. (4.110)–(4.111)
are written in the form
A(1) (0) C(1) (0) 2
χ= χ , q2 = (χ ) ,
B(1) B(1)
(4.114)
(1 − η)ρ1 (1 − μ)ρ2 (1 + η)T1 (1 + μ)12
f = C0 A(1) ,
(η + μ)(δ,+1)2 (η − μ)(δ−1)2
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch04 page 118
where
A(1) ≡ (η − μ)−4δ− A, B(1) ≡ (η − μ)−4δ− B,
C(1) ≡ (η − μ)−4δ− C,
(4.115)
and the functions A, B and C are given by Eq. (4.112).
(ii) δ+ < 0, δ− > 0: Then, in this case we write these solutions in the form,
A(2) (0) C(2) (0) 2
χ= χ , q2 = (χ ) ,
B(2) B(2)
(4.116)
(1 − η)ρ1 (1 − μ)ρ2 (1 + η)r1 (1 + μ)r2
f = C0 A(2) ,
(η + μ)(δ+ −1)2 (η − μ)(δ− +1)2
where the functions A(2) , B(2) and C(2) are related to the functions A, B
and C given by Eq. (4.112) via the relations,
A(2) = (η + μ)−4δ+ A, B(2) = (η + μ)−4δ+ B, C(2) = (η + μ)−4δ+ C.
(4.117)
(iii) δ± < 0: In this case, we have
A(3) (0) C(3) (0) 2
χ= χ , q2 = (χ ) ,
B(3) B(3)
(4.118)
(1 − η)ρ1 (1 − μ)ρ2 (1 + η)r1 (1 + μ)r2
f = C0 A(3) ,
(η + μ)(δ+ −1)2 (η − μ)(δ− −1)2
where
A(3) ≡ (η + μ)−4δ+ (η − μ)−4δ− A, B(3) ≡ (η + μ)−4δ+ (η − μ)−4δ− B,
(4.119)
C(3) ≡ (η + μ)−4δ+ (η − μ)−4δ− C.
Equations (4.110)–(4.119) show that in all the above cases, the function
f takes the form
f˜(η, μ)
f= , (4.120)
(η + μ)α+ (η − μ)α−
2 2
for example, in Eq. (4.110) we have α+ = (δ+ + 1) and α− = (δ− + 1) ,
etc.
The analysis carried out in the last two sections shows that for the
function f having the form of Eq. (4.120) the parameters n and m must be
chosen so that
n = (2 − α+ )−1 , m = (2 − α− )−1 . (4.121)
Without loss of generality, in the following we consider the above solutions
only for the cases in which δ± ≥ 0. The above solutions in the other cases
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch04 page 119
will have similar physical interpretations and properties. For more details,
we refer the readers to Wang (1991a).
When δ± ≥ 0, Eq. (4.121) becomes
n = [2 − (δ+ + 1)2 ]−1 , m = [2 − (δ− + 1)2 ]−1 . (4.122)
As usual, we take the solutions given by Eqs. (4.110)–(4.111) as valid
only in Region IV, then extend them back to the pre-collision regions by
means of the Khan–Penrose substitutions (3.11). Such extended solutions
represent the collisions of purely gravitational waves only in the cases where
n, m ≥ 1, which in the present case implies that
√
0 ≤ δ± < 2 − 1. (4.123)
From the discussion given in the last two sections, we can see that the
extended two-soliton solutions represent a variety of models of the collision
of two pure gravitational plane waves. For example, when n = 1 and m = 2
the solutions represent the collision of a variably polarized shock and a
variably polarized shock + impulsive wave, and so on.
Comparing with the seed solutions, we can see that in the two-soliton
case both of the two constantly (collinearly) polarized gravitational plane
waves of the seed solutions are generalized to variably (non-collinearly)
polarized ones, whereas in the one-soliton case only one of them is gen-
eralized to a variably polarized wave, while the other is still constantly
polarized.
But, these are the solutions given by Eqs. (4.23) and (4.24), if we replace
δ1 by δ1 + 1. Thus, replacing α by α = a + δ1 + 1 in Eqs. (4.50)–(4.52), we
obtain the asymptotic behavior of the solutions given by Eqs. (4.124) and
(4.125). Specifically, the conditions of Eq. (4.53) now become
2 1
(a + δ1 + 1) = , (4.126)
4
which is equivalent to
1 3
α=− , or α=− . (4.127)
2 2
On the other hand, in the non-collinear case we have Q1 Q2 = 0, which
allows us to distinguish the following cases.
(α) α < −1. Then, the solutions given by Eqs. (4.110)–(4.111) have the
same limit as their diagonal solutions given by (4.124) as t → 0+ . The
solutions with α = −3/2 are especially free of spacetime singularities on
the hypersurface t = 0 in both collinear and non-collinear cases. When
α < −3/2, we have p1 > 0. That is, the corresponding singularities are
astigmatic. When −3/2 < α < −1, we have p1 < 0, and the corresponding
singularities are anastigmatic.
as t → 0+ , where τ and pi are the same as they are defined in the corre-
sponding diagonal case. Equation (4.128) shows that the Kasner exponents
pi and the time τ are not changed relative to the corresponding diagonal
case, but the Kasner axes along which the p2 and p3 are defined are rotated
from (x2 , x3 ) to (X(z)
2 3
, X(z) ). Since in the present case p1 < 0, the nature
of the singularities on the hypersurface t = 0 is anastigmatic, as it is in the
collinear case. When α = −1/2, we have p1 = 0. Hence, the solutions with
α = −1/2 are free of spacetime curvature singularities on t = 0.
(γ) −1/2 < α < 1/2. In this case, the solutions have the limit
(3) (3)
ds2 = ε0 dτ 2 − ε1 τ 2p1 dz 2 − ε(3) 2τ 2p2 (dX(z)
2 2
)
(3)
− ε3 τ 2p3 (dX(z)
3 2
) , (4.129)
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch04 page 121
as t → 0+ , with τ and pi given exactly by Eqs. (4.50) and (4.51). Recall that
in the collinear case these solutions have the limit of Eqs. (4.49) and (4.51).
But, now the same limit is obtained with the replacement of α by α =
α + δ1 + 1. Thus, we find that in the present case the interaction between
different polarization modes of the two colliding plane gravitational waves
changes both the Kasner exponents and axes. The nature of singularities
in the present case is anastigmatic, whereas it is astigmatic in the collinear
case.
as t → 0+ , where
2
τ ≡ t(α−1) +3/4
, (4.131)
and
(α − 1)2 − 14 α − 12 3
−α
p1 ≡ , p2 ≡ , p3 ≡ 2
. (4.132)
(α − 1)2 + 34 (α − 1)2 + 34 (α − 1)2 + 34
It is easy to show that the exponents pi defined by Eq. (4.132) also satisfy
the Kasner relations (3.92).
As in the last case, the Kasner exponents as well as the Kasner axes are
all changed. In contrast to the collinear case, the nature of singularities on
t = 0 is anastigmatic. When α = 1/2, we have p1 = 0. Thus, the solutions
with α = 1/2 are free of spacetime singularities on the hypersurface t = 0
in the non-collinear case.
3 1 1 3
(i) α = − , (ii) α = − , (iii) α = , or (iv) α = , (4.134)
2 2 2 2
a A
χ= 1 − η 2 1 − μ2 ,
B
2a C
q2 = 1 − η 2 1 − μ2 , (4.135)
B
But, these are the solutions found by Ernst et al. (1987a, 1987b, 1988) by
using a different method. When δ1 = δ2 = 0, Eq. (4.122) yields n = m = 1.
Hence, from the previous analysis we can see that this subclass of solutions
in general represents the collision of variably polarized gravitational shock +
impulsive plane waves.
p2 + q 2 = 1 + l 2 . (4.138)
Substituting Eq. (4.137) into Eqs. (4.135) and (4.136), and considering the
fact that the parameter a vanishes in this case, we find that the correspond-
ing solution takes the form,
p2 1 − η 2 + q 2 1 − μ2
χ= ,
(1 − pη)2 + (l − qμ)2
2(qμ − lpη)
q2 = , (4.139)
(1 − pη)2 + (l − qμ)2
p2 1 − η 2 − q 2 1 − μ2
f = C0 1/4
,
[(1 − η 2 ) (1 − μ2 )] (η 2 − μ2 )
where C0 = 4C0 (1 − p)−2 is another arbitrary constant. The solution
given by Eq. (4.139) is the Nutku–Halil solution with the NUT parame-
ter l (Ernst, Garcia and Hauser, 1987b, 1988). Actually, it can be shown
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch04 page 124
1
E= [(pη + lqμ) + i(qμ − lpη)]. (4.140)
1 + l2
E = pη + iqμ, (4.141)
Case A.2: a = 1/2. Then, we find that the corresponding solution can be
written in the form
(1 − pη)2 + (l − qμ)2
χ= t,
p2
(1 − η 2 ) + q 2 (1 − μ2 )
p(q − lμ) 1 − η 2 − q(p − η) 1 − μ2
q2 = −2 ,
p2 (1 − η 2 ) + q 2 (1 − μ2 )
(4.142)
(1 − pη)2 + (l − qμ)2
f = C 0 ,
(η 2 − μ2 )
where C0 , p, q and l are now defined by (Economou and Tsoubelis, 1989),
It is easy to prove that p, q and l defined by Eq. (4.143) also satisfy the rela-
tion (4.138). As first noticed by Ernst, Garcia and Hauser (1987b, 1988), the
above solution is the Chandrasekhar–Xanthopoulos solution with the NUT
parameter l. In fact, setting l = 0, we shall rediscover the Chandrasekhar–
Xanthopoulos solution (Chandrasekhar and Xanthopoulos, 1986a). In this
case, the corresponding stationary axisymmetric solution is the Kerr–NUT
solution (Kramer et al., 1980).
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch04 page 125
1 + 2p η + η 2
χ= t, q2 = 2q μ,
1 − η2
(4.145)
2
1 + 2p η + η ,
f =C (q = 0),
0
η 2 − μ2
where
0 = 4C0 ,
C p2 + q 2 = 1, (4.146)
Q21 q 2
which is the Ferrari–Ibañez solution (Ferrari and Ibañez, 1988). Thus, the
Ferrari–Ibañez solution corresponds to the Taub–NUT solution (Hawking
and Ellis, 1973).
C0
C0 = , (4.147)
Q21
and then take the limit Q1 → ∞ (but keep C 0 finite), we find that the
solutions given by Eqs. (4.110)–(4.111) become
A (0)
χ= χ ,
B
(η − μ)2δ− +1 (0) 2
q2 = −2Q2 χ , (4.148)
B
(1 − η)ρ1 (1 − μ)ρ2 (1 + η)r1 (1 + μ)r2
f = C0 2 2 A,
(η + μ)(δ+ −1) (η − μ)(δ− +1)
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch04 page 126
2 1
rA = (a − δA ) − , (A = 1, 2). (4.150)
4
Equations (4.120), (4.121) and (4.148) show that the appropriate choice
of the parameters n and m now are
Following the discussion given in Section 4.2, it can be shown that this
subclass of solutions represents the same type of collision as the one-soliton
solutions given by Eqs. (4.56)–(4.58) (Tsoubelis and Wang, 1992).
In a similar fashion, we find that in the limit Q2 → ∞, the two-soliton
solutions become
A (0)
χ= χ ,
B
(η + μ)2δ+ +1 (0) 2
q2 = −2Q1 χ , (4.152)
B
(1 − η)ρ1 (1 − μ)ρ2 (1 + η)r1 (1 + μ)r2
f = C0 2 2 A,
(η + μ)(δ+ +1) (η − μ)(δ− −1)
where A, B and C are defined by
A ≡ 1 − μ2 (1 − η)2(a+δ1 ) (1 − μ)2(a+δ2 )
+ Q21 1 − η 2 (η + μ)4δ+ (1 + η)2(a−δ1 ) (1 + μ)2(a−δ2 ) ,
(4.153)
B ≡ (1 − μ)2 (1 − η)2(a+δ1 ) (1 − μ)2(a+δ2 )
+ Q21 (1 − η)2 (η + μ)4δ + (1 + η)2(a−δ1 ) (1 + μ)2(a−δ2 ) ,
Chapter 5
129
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch05 page 130
The analysis given in Section 3.2 [especially, see Eqs. (3.25) and (3.27)]
shows that, for any vacuum solution given in the interaction region
(Region IV), the extended solution obtained by the Khan–Penrose sub-
stitutions (3.11) will satisfy the Einstein vacuum equations inside the
pre-collision regions (Regions I–III). But across the hypersurfaces u = 0
and v = 0 the Ricci tensor will suffer the following discontinuities:
where
It is clear that when n = 1/2 the first term on the right-hand side of
Eq. (5.3) represents an impulsive shell of null dust with support on the
hypersurface u = 0, while when m = 1/2 the second term represents a
similar shell but with support on the hypersurface v = 0. Therefore, when
n = 1/2 and m = 1/2 the solutions given by Eqs. (4.23) and (4.24) will
represent the collision of two impulsive shells of null dust, each of which
may be accompanied by a constantly polarized gravitational plane wave.
When one of the two parameters n and m is equal to 1/2 and the other
one still satisfies the conditions (4.29), the above solutions will represent
the collision of an impulsive shell of null dust with a gravitational plane
wave.
Without loss of generality, let us consider the solutions given by
Eqs. (4.23) and (4.24) for the case where
1 1 1
(i) (n, m) = , , or (ii) (n, m) = ,≥ 1 , (5.4)
2 2 2
which implies that in all the cases to be considered there always exists an
impulsive shell of null dust that propagates along the hypersurface u = 0
toward the left in Fig. 3.1. From Eq. (4.28) we can see that n = 1/2 implies
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch05 page 132
δ+ = 0, or
1
δ1 = −δ2 , n = . (5.5)
2
Then, the combination of Eqs. (3.24), (4.23), (4.24) and (5.4) leads to the
following expressions for the non-vanishing Weyl scalars in Region IV,
3m
2 m−2 b1 (av + δ1 Y 3 )
ΨIV
0 (u, v) = 8m v
t 2
aδ1 [(2a + δ1 )v m + (a + 2δ1 )Y ]
− ,
(v m + Y )2
2b1 mv 2m−1 2aδ1 mv m−1
ΨIV
2 (u, v) = + , (5.6)
t2 Y (v m + Y )2
2b1 (δ1 v 3m + aY 3 )
ΨIV
4 (u, v) =
Y 3 t2
2aδ1 [(2a + δ1 )Y + (a + 2δ1 )v m ]
− .
Y 3 (v m + Y )2
Combining Eqs. (3.24), (3.26) with Eqs. (5.6)–(5.8), we find that the
non-vanishing Weyl scalars have the following behavior when across the
hypersurface u = 0 or v = 0.
I → III:
Ψ0 , Ψ2 are continuous,
(5.9)
ΨI–III = ΨIII (u = 0)H(u) + aδ(u).
4 4
II → IV:
Ψ0 is continuous,
ΨII–IV
2 = ΨIV
2 (u = 0, v)H(u), (5.10)
(a + kv ) m
ΨII–IV
4 = ΨIV
4 (u = 0, v)H(u) + δ(u).
X2
I → II:
Ψ2 , Ψ4 are continuous,
(5.11)
ΨI–II = ΨII (v = 0)H(v) + 2m(kv m−1 + av 2m−1 )δ(v).
0 0
III → IV:
Ψ4 is continuous,
2m(kv m−1 Y + av 2m−1 )
Ψ0III–IV = ΨIV
0 (u, v = 0)H(v) + δ(v), (5.12)
Y2
Ψ2III–IV = ΨIV
2 (u, v = 0)H(v),
where
k ≡ δ1 = −δ2 . (5.13)
To study the above solutions further, let us consider some typical cases.
Case A: m = 1/2. In this case, from Eqs. (4.28) and (5.4) we find
k = 0, (5.14)
for which we find m = 1/2, as it can be seen form Eq. (4.28). Then, setting
n = 1/2 and m = 1/2 in Eq. (5.3), we find that the non-vanishing compo-
nents of Tμν are given by
δ(u) δ(v)
Tuu = , Tvv = , (5.15)
1 − vH(v) 1 − uH(u)
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch05 page 134
which show that the corresponding solutions represent the symmetric col-
lision of a pair of impulsive shells of null dust, each of which is in general
accompanied by a constantly polarized impulsive + shock gravitational
plane wave. As a result of the collision, a Coulomb-like gravitational field
Ψ2 appears in the interaction region (Region IV) and a spacetime singular-
ity is finally developed on the hypersurface 1 = u + v.
When a = 0, Ψ0 and Ψ4 vanish everywhere, and the corresponding solu-
tion represents the collision of two pure impulsive shells of null dust. This
solution was first found by Dray and ’t Hooft (1986) and shown by Tsoubelis
(1989b) that it belongs to the Szekeres family of solutions (Szekeres, 1972).
When a = ±1/2, Eqs. (5.16) and (5.17) show that the two incoming
gravitational waves are impulsive ones. Thus, in this case the solutions
represent the collision of two impulsive shells of null dust, each of which is
accompanied by an impulsive gravitational plane waves. One of the remark-
able features in this case is that the interaction region (Region IV) is flat
and no Coulomb-like gravitational field appears in this region. Comparing
it, for example, with the Khan–Penrose solution (Khan and Penrose, 1971),
we are led to the conclusion that, when accompanied by null dust shells,
the collision of two impulsive gravitational waves does not necessarily give
rise to a Coulomb-like gravitational field in the interaction region. Since in
the present case Region IV is flat, no spacetime singularity develops on the
hypersurface 1 = u + v.
Before turning to the next case, let us note that the solutions presented
in this case were first obtained by Stoyanov (1979), but with an incor-
rect interpretation that the above solutions represented the collision of two
purely gravitational waves. The fact that Stoyanov’s interpretation cannot
be supported was first noticed by Nutku (1981) and is made clear from the
analysis presented above.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch05 page 135
δ(u)
Tuu = , Tvv = 0, (5.19)
1 − v 2 H(v)
and
1
ΨI4–III = a(4a2 − 1)H(u) + aδ(u), (5.20)
2
ΨI0–II = 6aH(v) + 2kδ(v), (5.21)
a[3v 2 + 12akv + (4a2 − 1)] a + kv
Ψ4II–IV = H(u) + δ(u), (5.22)
2(1 − v 2 )2 1 − v2
6a 2k
Ψ0III–IV = H(v) + √ δ(v). (5.23)
1−u 1−u
Case B.1: a = 0. In this case, Eqs. (5.20) and (5.21) reduce to ΨI4–III = 0
and ΨI0–II = δ(v), respectively, where = ±1. In fact, in this case the only
non-vanishing Weyl scalars are
vH(v)
Ψ0 (u, v) = δ(v), Ψ4 (u, v) = δ(u). (5.24)
1 − uH(u) 2(1 − v 2 )
Case B.2: a = ±1/2 ≡ /2. In this case, Eqs. (5.20) and (5.21) reduce to
1
ΨI–III
4 = δ(u), (5.25)
2
ΨI–II
0 = 3 H(v) + δ(v), (5.26)
where and can be chosen independently, always we have || = | | = 1.
Then, we can see that the present model represents the collision of a shock
wave with an impulsive shell of null dust, each of which is accompanied by
an impulsive gravitational wave. In this case, the formation of a spacetime
singularity along the hypersurface u = 1 − v 2 is inevitable.
Case B.3: a = −2k. Since in this case we have
3
ΨI–III
4 = − εH(u) − εδ(u), (5.27)
2
ΨI–II
0 = −6εH(v) + εδ(v), (5.28)
we find that the corresponding solutions represent the collision of two impul-
sive + shock gravitational waves, one of which is accompanied by an impul-
sive shell of null dust. The collision is such that the spacetime is free of
singularity on the hypersurface 1 − u − v 2 = 0. In fact, the Coulomb-like
gravitational field that develops in Region IV after the collision remains
finite as one approaches the hypersurface t = 1 − u − v 2 = 0, and is given by
1
ΨIV
2 (u, v) = − √ √ . (5.29)
[ 1 − u + v]2 1 − u
Thus, we have
1
ΨIV 2
2 (u = 1 − v , v) = − . (5.30)
4v 3
√
Case C: m = 2. In this case, we have k = 3/8 and that
1
ΨI4–III = a(4a2 − 1)H(u) − aδ(u), (5.31)
2
ΨI0–II = 4kH(v). (5.32)
Thus, in the present case the solutions represent the collision of a shock
gravitational wave incident from the left with an impulsive shell of null
dust incident from the right. The latter is accompanied by a gravitational
shock + impulsive wave, provided a = 0, ±1/2. The behavior of this subclass
of models after collision is similar to the ones obtained in Cases A and B
above. Specifically, the development of a spacetime singularity along the
1 − u − v 4 = 0 is inevitable, except for the cases where a = −k ± 1/2.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch05 page 137
Case D: m = 4. The shock wave incident from the left in Case C now
is replaced by a gravitational plane wave with a smooth wavefront. This
follows from the fact that in this case we have
1
ΨI4–II = a(4a2 − 1)H(u) − aδ(u), (5.34)
2
Comparing Eqs. (5.33) and (5.34) with Eqs. (5.32) and (5.31), we find that
the current models represent the same kind of collisions as the m = 2
models, except that the shock gravitational wave incident in Region II is
now replaced by a gravitational wave with a smooth wavefront. Thus, in
the following we do not go to detail, except for mentioning that in the
present case the second gravitational wave production is illustrated most
clearly, although it arises essentially in all the spacetime models constructed
above. Consider in this direction, the a = 0 case. For this particular model
Ψ4I–III = 0. Therefore, the null dust shell incident from the right is not
accompanied by any gravitational radiation. According to Eq. (5.35), on
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch05 page 138
Thus, in this case the solution represents the collision of an impulsive gravi-
tational plane wave and an impulsive shell of null dust. Since the spacetime
is flat in Region IV, the hypersurface t = 0 is free of spacetime curvature sin-
gularity, except at the two focusing points (u, v) = (0, 1) and (u, v) = (1, 0).
The study of polarization of the two impulsive gravitational waves given
by Eq. (5.44) in the next subsection will reveal that this solution is actually
the Babala solution referred as Case B.1 in Section 5.1.
Case 2: a = 1/2. In this case, if we define the functions F (x, y) and G(x, y)
as
Taking Eq. (5.38) into account, we find that this solution represents the col-
lision of a variably polarized gravitational shock + impulsive wave with an
impulsive shell of null dust, which is accompanied by a constantly polarized
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch05 page 141
gravitational impulsive wave. In this model, Region III is flat, while Region
II is curved due to the presence of the variably polarized gravitational wave.
As t → 0+ , all the non-vanishing Weyl scalars become unbounded. Thus,
the spacetime now is singular on the hypersurface t = 0.
Case 3: a = 1. Then, we find that the non-vanishing Weyl scalars are given
by
1
Ψ0 (u, v) = ΨIV
0 (uH(u), v)H(v) + {(1 − q 2 ) − i2q}δ(v),
(1 + q 2 )Y
1
Ψ4 (u, v) = ΨIV
4 (u, vH(v))H(u) + {(2 + vH(v))(1 − vH(v))8
2tAD
+ q 2 (1 + vH(v))3 [q 2 (2 − vH(v))(1 + vH(v))5
+ 4(1 + 3v 2 H(v))(1 − vH(v))3 ] (5.50)
2 4 2
+ i2qv(1 − v )H(v)[(1 − v) (1 + 4v + v )
+ q 4 (1 + v)4 (1 − 4v + v 2 )]}δ(u),
Ψ2 (u, v) = ΨIV
2 (u, v)H(u)H(v),
where
6 3
ΨIV
0 (u, v) = F (Y, v), ΨIV
4 (u, v) = F (v, Y ),
A2 D 2Y 2 A2 D (5.51)
6
ΨIV
2 (u, v) = G(Y, v),
A2 Y
and the functions D, F and G now are defined by
D(u, v) ≡ [A2 + 16q 2 v 2 t2 Y 2 ]1/2 ,
F (x, y) ≡ (x − y)10 − q 2 (x2 − y 2 )(x − y)4 (5x4 + 30x3 y + 8x2 y 2
+ 2xy 3 + 3y 4 ) − q 4 (x2 − y 2 )(x + y)4 (5x4 − 30x3 y
+ 8x2 y 2 − 2xy 3 + 3y 4 ) + q 6 (x + y)10
+ i4qx[(x − y)7 (y 2 − 2xy − x2 ) + 2q 2 y(x − y)3 (5x2 + y 2 )
(5.52)
− q 4 (x + y)7 (y 2 + 2xy − x2 )],
G(x, y) ≡ (x − y)6 − 12q 2 xy(x2 − y 2 )2 − q 4 (x + y)6
− i2q[(x − y)4 (y 2 + 4xy + x2 )
+ q 2 (x + y)4 (y 2 − 4xy + x2 )].
From Eqs. (5.3), (4.39) and (5.50), we can see that the solution in the
present case represents the collision of a variably polarized gravitational
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch05 page 142
impulsive + shock wave with the same type wave but constantly polarized.
The latter has an impulsive shell of null dust as its leading front.
As t → 0+ , Eqs. (5.51) and (5.52) show that all the non-vanishing
Weyl scalars keep finite provided q = 0. That is, the spacetime is free of
spacetime curvature singularity on the hypersurface t = 0 when q = 0.
However, when q = 0, it is easy to see that the non-vanishing Weyl scalars
become unbounded as t → 0+ . Thus, in the latter case the spacetime is
singular on t = 0. Consequently, a close relation is established between the
collision involving impulsive shells of null dust, on the one hand, and the
collision involving only gravitational waves, on the other. For the latter
type of collision, the effect of polarization of colliding plane gravitational
waves on the formation of singularities was already studied in Sections 4.2
and 4.3.
Before proceeding to the next section, we note that solutions that
represent the collision of an impulsive shell of null dust and a variably polar-
ized gravitational plane wave were first studied by Feinstein and Senovilla
(1989). Starting with the assumption that the non-diagonal term q2 (in
their notation it is denoted by ω) depends only on a null coordinate v in
the interaction region, they found a unique solution of the Einstein vacuum
equations for a given arbitrary q2 (v). After extending the solutions back
to the pre-collision regions, they found an impulsive shell of null dust to
appear on the hypersurface u = 0. It was the latter that enabled them to
interpret their solutions as representing the collision of a plane gravitational
wave with an impulsive shell of null dust.
solution are very complicated. It is difficult to see the physics. Second, from
this solution we shall obtain the similar results about the polarizations of
colliding plane gravitational waves as we shall get from the first two.
With the above in mind, let us first consider Case 1 solution. The non-
vanishing Weyl scalars in this case are given by Eq. (5.44). Since on the
hypersurface v = 0, we have W = 0 [see Eq. (5.41)]. Hence, integrating
Eq. (3.61) we find that
ϕ00 = 0. (5.53)
Eq. (5.57) will take the exact form used by Babala (1987).
Now let us turn to the solution of Case 2 given in Subsection 5.1.2. We
first note that in this case the Weyl scalars Ψ0 and Ψ4 consist of two parts:
the shock part and the impulsive part. In the following, we consider them,
separately.
Let us first consider the impulsive part. Following the discussion given
in the last case, we find that
1 2q
θ0Im = tan−1
,
2 1 − q2
(5.59)
1 2qv
θ4Im = tan−1 H(v).
2 (1 − v) + q 2 (1 + v)
Thus, the polarization angle for the impulsive part of Ψ0 remains constant
even after the collision. The reason is that, similar to the last case, ΨIm
0
does not interact with any of the other components of both gravitational
and matter fields. However, for the Ψ4 wave the situation is different. The
interaction between the impulsive part of Ψ4 and the impulsive shell of null
dust ΦIm Im
22 is such as to make the polarization angle θ4 change along the
v-axis according to Eq. (5.59).
In a similar fashion, we find that for the shock part of Ψ0 and Ψ4 the
polarization angles θ0sh and θ4sh are given, respectively, by
1 I(v, Y )
θ0sh = − tan−1
H(v),
2 J(v, Y )
(5.60)
1 I(Y, v)
θ4sh = tan−1 H(u)H(v),
2 J(Y, v)
and Marshman (1979), Tabensky and Taub (1973) and Taub (1988b), we
find that solutions of a massless scalar field can be obtained from vacuum
ones by simply setting
(M, U, V, W ) = (Mv + Ω, Uv , Vv , Wv ), (5.62)
where Mv , Uv , Vv and Wv are solutions of the Einstein vacuum equations,
and the function Ω satisfies
φ,u
Ω,u = , (5.63a)
U,u
φ,v
Ω,v = , (5.63b)
U,v
Ω,uv = φ,u φ,v , (5.63c)
where the scalar field φ satisfies the massless scalar field equation (3.80). If
we consider, on the other hand, the metric
(0) (0) (0)
ds2 = 2e−M dudv − e−U [eV (dx2 )2 + e−V (dx3 )2 ], (5.64)
we find that the corresponding Einstein vacuum equations are given by [see
Eqs. (3.8a)–(3.8e)],
2 (0) (0)
2U,vv − U,v + 2U,v M,v = V,v2 , (5.65a)
(0) (0)
2U,uu − U,u2 + 2U,u M,u = V,u2 , (5.65b)
(0)
2M,uv = U,u U,v , (5.65c)
and
U,uv − U,u U,v = 0, (5.66)
(0)
2V,uv − U,u V,v(0) − U,v V,u(0) = 0. (5.67)
Introducing the function Ω(0) by
M (0) = N (0) + Ω(0) , (5.68)
where
3
N (0) ≡
U − ln |2U,u U,v |, (5.69)
2
we find that Eqs. (5.65a)–(5.65c) can be written as
(0) (0)
(0) V,u (0) V,v
Ω,u = , Ω,v = ,
2U,u 2U,v (5.70)
(0) 1 (0) (0)
Ω,uv = V,u V,v .
2
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch05 page 148
Comparison of Eqs. (3.80) and (5.63a)–(5.63c) with Eqs. (5.67) and (5.70)
leads to the following theorem.
Theorem (Tsoubelis and Wang, 1991). Let (M (0) , U, V (0) ) be a solution
of the Einstein vacuum equations corresponding to the metric (5.64), Then
(M, U, V, W ) = (Mv + λ2 (M (0) − N (0) ) + C0 , U, Vv , Wv ), (5.71)
(0)
with N given by Eq. (5.69) is a solution
√ of the Einstein equations coupled
(0)
with a massless scalar field φ = λV / 2, where λ and C0 are constant.
Note that when Vv = Wv = 0, the Einstein vacuum equations give the
following solution for the function Mv :
Mv = N (0) + C1 , (5.72)
where C1 is another arbitrary constant.
As an application of the above theorem, we use the vacuum solutions
given by Eqs. (4.23) and (4.24) as the solutions of (M (0) , U, V (0) ), and for
the sake of convenience we set the functions Vv and Wv to be zero
Vv = Wv = 0. (5.73)
Then, the function Mv has the solution given by√Eq. (5.72). By choosing
the parameter λ introduced in Eq. (5.71) to be 2, we finally obtain the
following solutions:
1−η
φ = V (0) ≡ a ln[(1 − η 2 )(1 − μ2 )] + δ1 ln
1+η
1−μ
+ δ2 ln , (5.74)
1+μ
and
e−U = 1 − u2n − v 2n ,
2 2
e−M = (1 − η)b1 (1 − μ)b2 (1 + η)c1 (1 + μ)c2 X −2δ+ Y −2δ− , (5.75)
V = W = 0,
but now bA and cA are defined by
1 1
bA = 2(a + δA )2 − , cA = 2(a − δA )2 − , (A = 1, 2), (5.76)
4 4
and the parameters n and m are related to δ± via the relations,
1 1
n= 2 ), m= 2 ). (5.77)
2(1 − δ+ 2(1 − δ−
As usual, taking the above solutions as valid only in the interaction region
(Region IV), we can extend them to the pre-collision regions by means
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch05 page 149
Case (ii): ε < 0. In this case, φ,μ defines a space-like unit vector, say Zμ ,
via the relation
φ,μ
Zμ = . (5.84)
(−2ε)1/2
With similar considerations, we find that the metric in the present case is
given by Eq. (5.82), while the energy–stress tensor is given by
Case (iii): ε = 0: In this case, φ,μ is a null vector and Tμν takes the form
that we have been calling a null dust cloud propagating in the direction
defined by φ,μ .
The significance of the above observations will be made clear as we are
going to reveal the fact that in the solutions under considerations all the
above possibilities arise in the interaction region (Region IV) [see Fig. 5.1].
Fig. 5.1. The projection of the colliding null dust spacetime onto the (u, v)-plane.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch05 page 151
ε = eM Ruv
4nmeM tan ζ tan ξ
= [a sin 2ζ + δ+ cos ζ cos ξ + δ− sin ζ sin ξ]
t2 (sin ζ)1/n (sin ξ)1/m
×[a sin 2ξ + δ− cos ζ cos ξ + δ+ sin ζ sin ξ], (5.86)
where
4a2 1
Ruu = H(u) + δ(u), (5.89a)
[1 − u − vH(v)]2 1 − vH(v)
4a2 1
Rvv = 2
H(v) + δ(v). (5.89b)
[1 − uH(u) − v] 1 − uH(u)
After the collision, the component Ruv appears, which makes the scalar
field energetically equivalent to a perfect fluid with the energy density and
pressure given by
4a2 4a2
p=ε= 2 = . (5.90)
(1 − u − v)(8a +3)/2 t(8a2 +3)/2
The last equation shows explicitly that this class of solutions develops
a spacetime curvature singularity along the hypersurface t = 0. In fact,
the singularity happens in all the colliding wave models considered in this
section.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch05 page 152
superposition of the incoming waves after collision is such that one cannot
regard the interaction region as consisting of a fluid of a certain type. In
particular, assuming that δa < 0 with |δ/2a| < 1, from Eqs. (5.95a)–(5.95c)
we find that Ruv = 0 along the following two curves in Region IV. The first
is defined by the relation
δ
sin 2ζ = − cos(ζ − ξ), (5.96)
a
and has its origin at the point (u0 , 0) with u0 = |δ/2a|. Similarly, the second
curve is defined by
δ
sin 2ξ = − cos(ζ − ξ), (5.97)
a
and starts at (0, v0 ) where v0 = u0 . Across these curves the sign of
Ruv changes. Therefore, Region IV is split into subregions occupied by
a fluid whose type alternates among the cases (i) and (ii) considered after
Eq. (5.80), as illustrated in Fig. 5.1.
Case D: a = 0, δA = δ = 0. In this case, we find that
(1 − η)(1 − μ)
φ = δ ln , (5.98)
(1 + η)(1 + μ)
and
16n2 δ 2 2 2(n−1)
Ruu = Y u H(u), (5.99a)
t2
4δ 2 2n 1
Rvv = 2 2
u H(u)H(v) + 2n
δ(v), (5.99b)
t Y 1 − u H(u)
8nδ 2 2n−1
Ruv = u H(u)H(v), (5.99c)
t2
which show that in the present case the solutions represent the collision of
a null dust cloud with an impulsive shell of null dust. After the collision,
a current moving along the v = constant lines toward the right-hand side
is developed, which together with the Ruv component given by Eq. (5.99c)
again makes the scalar field in the interaction region energetically equivalent
to a perfect fluid with the equation of state, p = ε.
However, further considerations (Tsoubelis and Wang, 1991) showed
that the origin of the right-hand moving current is difficult to under-
stand if one sticks to the fluid picture, since before the collision there is
no incoming current moving toward the right-hand side as follows from
Eq. (5.99c). Instead, if we consider the question in terms of the scalar field,
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch05 page 154
from the Klein–Gordon equation (3.80) we can see that this is because of
the backscattering. Actually, in terms of the latter, the present model rep-
resents the collision of a scalar plane wave incident from the right-hand
side, which collides at (u, v) = (0, 0) with an impulsive wave (shell of null
dust). The focusing that results from this encounter is equivalent to the
scalar wave pulse’s entering a spacetime region where the metric depends
on both u and v. But, then the Klein–Gordon equation (3.80) implies that
right moving waves are generated, i.e. backscattering occurs. This explains
the appearance of the right moving current in Region IV in terms of the
null current picture.
where α(u) and β(v) are arbitrary functions of their indicated arguments.
On the other hand, from Eqs. (3.8a)–(3.8e), (4.20a) and (4.20b) we find
that corresponding to the metric (5.100) the condition (5.103) now reads
For the function U given by Eq. (5.105), we find that Eq. (5.103) has the
general solution
1
M = − U − A(u) − B(v), (5.107)
2
where A, B are other two arbitrary functions of their indicated arguments.
Corresponding to the metric (5.100) with the function M and U given
by Eqs. (5.107) and (5.105), it can be shown that the Ricci tensor Rμν is
given by
α A − α ˆ β B − β
Rμv = l̂μ lν + n̂μ n̂ν , (5.108)
α+β α+β
where the prime denotes the ordinary differentiation with respect to the
indicated argument. Thus, the corresponding energy–stress tensor can be
written in the form
Tμν = Tμν
L
+ Tμν
R
, (5.109)
where
α A − α
Tμν
L
= εL l̂μ ˆlv ≡ l̂μ l̂v , (5.110a)
α+β
β B − β
Tμv
R
= εR n̂μ n̂v ≡ n̂μ n̂ν . (5.110b)
α+β
Equations (5.109) and (5.110a)–(5.110b), together with Eq. (1.21), yield
Tμν
L;ν
= 0 = Tμν
R;ν
. (5.111)
It follows that, assuming both εL and εR be positive, the source given by
Eqs. (5.109) and (5.110a)–(5.110b) consists of a pair of oppositely moving
null dust clouds, one of which has the conserved energy density εL , while
the other has the conserved energy density εR .
In order to consider the solutions representing the collision of two null
dust clouds, let us assume the functions A and B introduced in Eq. (5.107)
such that
A(0) = B(0) = 0. (5.112)
We also assume that α and β, which specify the function U via Eq. (5.105),
are chosen such that
1
α(0) = β(0) = . (5.113)
2
Then, we can see that the conditions (3.22) are now satisfied, and the Khan–
Penrose substitutions (3.11) can be used. Replacing the pair (u, v) by the
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch05 page 156
one (uH(u), vH(v)) in Eqs. (5.105) and (5.107), we find that the extended
solutions are given by
eA(u)H(u)+B(v)H(v)
e−M = ,
[1 − f (u)H(u) − g(v)H(v)]1/2 (5.114)
e−U = 1 − f (u)H(u) − g(v)H(v),
where
1
mμ = √ e−U/2 (δμ2 + iδμ3 ). (5.119)
2
Provided that it satisfies the Maxwell equation (1.96), Fμν represents a null
electromagnetic field with the energy–stress tensor Tμν given by
1
Tμν
e
= Fμλ F λ ν − gμν Fρλ F λρ = 2Φ0 Φ0 n̂μ n̂ν . (5.120)
4
Comparing Eq. (5.116) with Eq. (5.120), we conclude that the second term
in the right-hand side of Eq. (5.116) can be attributed to an electromagnetic
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch05 page 157
and a pulse of neutrino wave given, respectively, with the field strengths,
1/2
2mG(v) 2
Fμν = 2v m−1 H(v)n[μ δν] , (5.138)
1 − v 2m
1/2
un−1 2nAF (u)eU
φA = H(u)eiu OA . (5.139)
2 1 − u2n
The above class of solutions covers a large variety of physically inter-
esting cases, which are obtained by choosing the corresponding parameters
appropriately. When the tetrad (p, q, r, s) vanishes and n = m = 1, for
example, we recover Griffiths’ solution (Griffiths, 1976b), which represents
the collision of a pair of constant profile shock waves made of photons
and neutrinos, respectively. When n = 2 and m = 1, the corresponding
solutions represent the collision of an electromagnetic shock wave and a
neutrino wave with smooth wavefront, etc.
Without going to more detailed analyses for each of specific models, we
turn to the issues, regarding to the ambiguity of the outcome of collisions
of matter fields, raised in Section 5.2. Choosing the parameters p, q, r, s, n
and m as
where a and δA will be chosen as the same constants that specify the
family of scalar wave solutions analyzed in the last section, we obtain a
three-parameter family of the solutions given by Eqs. (5.131)–(5.134). It
is then easy to verify that in the pre-collision regions (Regions I–III), the
metric resulting from the choice of Eq. (5.140) is identical to the one cor-
responding to the scalar wave models studied in Section 5.2, even though
the solutions differ from each other in the interaction region (Region IV).
Thus, Taub’s results about the ambiguity from planar colliding null dust
clouds with impulsive waves now extend to solutions that are free from
such impulsive components. That is, given the metric (solutions of the
Einstein field equations) in Regions I–III, one would not be able to pred-
icate what is the outcome of the collision. The colliding process becomes
determinative, however, once the energy–stress tensor or equivalently, the
type of interaction is pre-assigned in Region IV. In the present case, the
resolution of the ambiguity was resolved by constructing the corresponding
energy–stress distributions from the well-defined matter fields, neutrino and
electromagnetic.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch05 page 160
cos(au − bv)
ds2IV = 2dudv − [1 − sin2 (au) − sin2 (bv)] (dx2 )2
cos(au + bv)
cos(au + bv)
+ (dx3 )2 , (5.141)
cos(au − bv)
where a and b are two real constants. It is interesting to note that in the cur-
rent case the focusing surface, 1−sin2 (au)−sin2 (bv) = 0 (or au + bv = π/2),
is not singular, but a Cauchy horizon (Bell and Szekeres, 1974). Clarke and
Hayward (1989) showed that extensions beyond this surface are not unique,
and similar to the vacuum case (Yurtsever, 1987), it was argued that such
a horizon is also not stable (Konkowski and Helliwell, 1991; Gürses and
Halilsoy, 1982; Gurtug and Halilsoy, 2000). Another interesting aspect of
the solution is that metric (5.141) is conformally flat,
q2 2
ds2IV = [dt − dr2 − r2 (dθ2 + sin2 θdφ2 )], (5.142)
r2
where 2ab = q −2 , and the mapping between the two sets of coordinates is
given explicitly by (Griffiths, 1991),
1 y
t + r = coth sech−1 cos(au + bv) − ,
2 2q
1 −1 y
t − r = −tanh sech cos(au + bv) + , (5.143)
2 2q
π 1
θ = au − bv + , φ= x.
2 2
Using the Khan–Penrose substitutions, one can easily find the metric in
Regions I–III,
⎧
⎨ 2dudv − cos2 (au)dΣ20 , Region III,
2
ds = 2dudv − cos2 (bv)dΣ20 , Region II, (5.144)
⎩
2dudv − dΣ20 , Region I,
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch05 page 161
where dΣ20 ≡ (dx2 )2 + (dx3 )2 . Clearly, Region I is flat, while in each of the
two regions (Regions II and III) an electromagnetic shock wave is present,
Φ0 = bH(v), Φ2 = aH(u), (5.145)
while all the Weyl scalars ΨA (A = 0, 2, 4) vanish in these regions as well as
across their boundaries u = 0 and v = 0. Thus, the Bell–Szekeres solution
represents the collision of two pure electromagnetic plane waves. However,
as in the mixed cases presented in the last section, two impulsive gravita-
tional waves are created due to the collision, one moves along the hyper-
surface u = 0 and the other along the one v = 0 (Bell and Szekeres, 1974),
Ψ0 = −b tan(au)H(u)δ(v), Ψ4 = −a tan(bv)H(v)δ(u). (5.146)
In the general case, after Chandrasekhar and Xanthopoulos (1985a)
laid down the foundation of the problem, various solutions were obtained
subsequently (Griffiths, 1985, 1990; Halilsoy, 1988a, 1988b, 1989a, 1989b;
Garcı́a-Dı́az, 1988, 1989; Papacostas and Xanthopoulos, 1989; Li and Ernst,
1989; Bretón et al., 1998; Hogan, Barrabés and Bressange, 1998; Hogan
and Walsh, 2003). In particular, Chandrasekhar and Xanthopoulos (1985a)
showed explicitly how one can generalize the Khan–Penrose and Nutku–
Halii solutions to include electromagnetic plane shock waves, in which a
spacetime curvature singularity is finally formed on the focusing surface.
They also provided the generalization of the Bell–Szekeres solution to the
case with two non-collinearly polarized electromagnetic plane shock waves
(Chandrasekhar and Xanthopoulos, 1987b), and showed that a Cauchy hori-
zon is finally formed on the focusing surface.
Chapter 6
In the past three decades or so, the collisions of gravitational plane waves
coupled with various matter fields have been extensively studied, and very
important insight about the nonlinearity of the Einstein field equations has
been obtained.
One of the remarkable features is that the internal spacetimes of all
the known black holes (Stephani et al., 2009) are isometric to the interact-
ing regions of two colliding plane gravitational waves with or without the
presence of matter fields. The fundamental reason is that in both cases the
Einstein field equations can be written in terms of the Ernst equation (3.44)
(Chandrasekhar and Ferrari, 1984).1 So, in this chapter we shall present a
brief review on these fascinating phenomena.
1 In the cylindrically symmetric case, when the Abelian G2 group acts orthogonally tran-
sitively, the metric can be written in the form of Eq. (3.41), but with the x3 -coordinate
having closed orbits (Bronnikov, Santos and Wang, 2019). Then, the Einstein vacuum
field equations reduce to the “cylindrical” Ernst equation (Chandrasekhar, 1986). There-
fore, such isometries should exist between cylindrical spacetimes (Wang, 2003) and the
interiors of black holes, too.
163
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch06 page 164
r−m
η≡ , μ ≡ cos θ, x2 ≡ t, x3 ≡ mφ, C0 = m2 , (6.3)
m
we find that the metric (6.2) takes the form,
−1
2 2m 2 2m
ds = 1 − dt − 1 − dr2 − r2 (dθ2 + sin θ2 dφ2 ), (6.4)
r r
which is precisely the form of the Schwarzschild black hole solution written
in the spherically symmetric coordinates (t, r, θ, φ). Note that the focusing
hypersurface η = 1 or t = 0 corresponds to the Killing-horizon r = 2m, in
other words, now the focusing hypersurface is not singular, and instead it
represents a horizon. This is consistent with the analysis presented between
Eqs. (4.45) and (4.47). In particular, the choice of Eq. (6.1) satisfies the
non-singular condition (4.47).
It should be noted that the above equivalence holds only locally, as one
can see from Eq. (6.3), which tells us that it is the x2 -coordinate that is
identical to the t-coordinate, while the space-like coordinate μ is identical
to the angular coordinate φ, and the time-like coordinate η to r, as inside
the black hole, r becomes time-like, and t is space-like.
On the other hand, from Eq. (4.28) we find that the choice of
Eq. (6.1) yields,
m = n = 1, (6.5)
for which the incoming gravitational waves consist of two parts, the impul-
sive and shock parts, as one can see from the analysis given between
Eqs. (4.34) and (4.42). In particular, along the wavefronts u = 0 and
v = 0, the two wave components Ψ0 and Ψ4 takes the forms,
3
ΨII
0 (v) = , (Region II),
(1 + v)2 (1 − v)
(6.7)
3
ΨIII
4 (u) = , (Region III).
(1 + u)2 (1 − u)
1
a=+ , δ1 = 1, δ2 = 0, (6.8)
2
which can be obtained from the metric (6.2) by the replacement η → −η.
However, now the coordinate transformations (6.3) to the internal region
of the Schwarzschild black hole solution (6.4) become,
m−r
η≡ , μ ≡ cos θ, x2 ≡ t, x3 ≡ mφ, C0 = m2 , (6.10)
m
that is, all of them are the same as these given by Eq. (6.3) except for
η → −η. It is this difference that now the focusing hypersurface t = 0 or
η = 1 corresponds to r = 0, instead of r = 2m, as given in the last case.
Therefore, in the present case the hypersurface t = 0 (η = 1) becomes
a spacetime curvature singularity. Note that this is also consistent with
the analysis given between Eqs. (4.45) and (4.47), and in particular the
non-singular condition (4.47) is not satisfied.
From Eq. (4.28) we find that the choice of Eq. (6.8) also yields m = 1 =
n, and along the wavefronts u = 0 and v = 0, the two gravitational wave
components Ψ0 and Ψ4 take the forms,
Z −1 eV sinh W eV
Z ≡ χ + iq2 , E≡ , χ≡ , q2 ≡ . (6.14)
Z +1 cosh W cosh W
The Nutku–Halil solution represents the collision of two pure gravitational
impulsive waves with uncorrelated (non-collinear) polarizations. It is a
direct generalization of the Khan–Penrose solution (Khan and Penrose,
1971), which represents the collision of two pure gravitational impulsive
waves with correlated (collinear) polarizations, and the corresponding Ernst
potential is given by Eq. (6.13) with q = 0.2
As first noticed by Chandrasekhar and Ferrari (1984), the Kerr solution
(Kerr, 1963) also follows from the same Ernst potential (6.13) (Chan-
drasekhar, 1983). But, in the stationary case the Ernst potential is related
to the metric coefficients through the function Φ̃ and Ψ̃ (Chandrasekhar,
2 Barrabés, Bressange and Hogan showed that the Khan–Penrose and Nutku–Halil solu-
tions can be also obtained from Einstein’s vacuum field equations as an initial value
problem (Barrabés, Bressange and Hogan, 1999).
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch06 page 167
where
Δ ≡ r2 − 2M r + J 2 , ρ≡ r2 + J 2 cos2 θ. (6.20)
The above solution is exactly the Kerr solution written in the Boyer–
Lindquist coordinates with M and J being the mass and angular momen-
tum parameters (Chandrasekhar, 1983).
Since −1 ≤ η ≤ +1, from Eq. (6.18) we find that this implies
r− ≤ r ≤ r+ , (r± = M ± M 2 − J 2 ), (6.21)
that is, the interacting region of the two pure colliding plane gravitational
waves is isometric to the ergo-sphere of the Kerr black hole, in which the
time-like Killing vector ∂t becomes space-like (in addition to the space-like
Killing vector ∂φ ). As a result, in this region there are two space-like Killing
vectors, which are consistent with the requirement of spacetimes with plane
symmetry (Stephani et al., 2009), although in the Kerr spacetime, one of
them has a closed orbit [cf. Eq. (6.18)]. Therefore, such an equivalence can
be only local.
On the other hand, setting η = 1, which corresponds to the focusing
hypersurface t = 0, from Eq. (6.18) we find that
r(t = 0) = M ± M 2 − J 2 . (6.22)
Thus, if the “−” sign in Eq. (6.18) is chosen, the focusing surface corre-
sponds to the null surface r = r+ , which is the location of the event horizon
of the Kerr black hole, and when the “+” sign in Eq. (6.18) is chosen, the
focusing surface corresponds to the null surface r = r− , which is the loca-
tion of the Cauchy horizon of the Kerr black hole. In each of the two cases,
the spacetime is not singular at t = 0, and extensions beyond these surfaces
are needed. Clearly, if we require the extension is analytical, then we shall
get a global structure quite similar to that of the Kerr black hole. However,
in the case r(t = 0) = r− the focusing hypersurface is a Cauchy horizon,
which is not stable against small perturbations (Yurtsever, 1987, 1988a),
so such extensions might not be needed in this case.
It is interesting to note that the solution (6.17) is precisely the two-
soliton solution with
1
a = , δ1 = δ2 = 0, ⇒ n = m = 1, (6.23)
2
given by Eqs. (4.142) and (4.143) by setting the NUT parameter l to zero.
When l = 0, the interacting region of the colliding plane gravitational
wave is isometric to the interior of the Kerr–NUT solution. When q = 0,
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch06 page 169
the solution reduces to that given by Eq. (4.145), in which the interacting
region is isometric to the interior of the Taub–NUT solution, first studied
by Ferrari and Ibañez (1988).
Finally, we note that in the current case we have α ≡ a + δ1 = 1/2,
which satisfies one of the conditions given in Eq. (4.134), and as a result
the corresponding solution is not singular at the focusing surface t = 0.
This is well consistent with the above analysis.
region of the two colliding waves is a locally known solution of the Einstein–
Maxwell electrovacuum equations of Petrov type D, and the collision results
in the formation of a Cauchy horizon.
In addition, Griffiths and Halburd (2007) showed that part of the
C-metric spacetime inside the black hole can be also interpreted as the
interaction region of two colliding plane waves with aligned linear (collinear)
polarizations. The focusing surface t = 0 is not singular, and instead repre-
sents a Cauchy horizon. Note that the C-metric was first found by Weyl in
1917 (Weyl, 1917), and subsequently rediscovered by many authors. But,
it was Kinnersley and Walker (1970) and Bonnor (1983) who first showed
that its analytic extension represents a pair of black holes that accelerate
away from each other due to the presence of a strut along the symmetry
axis.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch07 page 171
Chapter 7
Concluding Remarks
In this book, we study the collision and nonlinear interaction of two pure
gravitational plane waves, or a gravitational plane wave with a matter wave,
or two matter waves, within the framework of Einstein’s theory of general
relativity, by using the distribution theory, first introduced by Taub to gen-
eral relativity in 1980 (Taub, 1980). The method is mathematically equal
to Israel’s junction conditions (Israel, 1966, 1967), when the singular hyper-
surface is time-like or space-like. When the surface is null, it is equivalent to
the analysis of Barrabés (see Barrabés, 1989; Barrabés and Hogan, 2003).
Thus, they can be considered as complementary to each other.
The advantage of the studies of spacetimes for colliding gravitational
and/or matter plane waves is that they can be investigated analytically,
thanks to the development of nonlinear differential equations in 1960s
(Gardner, Green, Kruskal and Miura, 1967; Whitham, 1974), and appli-
cations to Einstein’s general relativity in 1970s and 1980s (Belinsky and
Verdaguer, 2001; Griffiths and Podolský, 2009; Stephani et al., 2009),
including the soliton technique of Blelinsky and Zakharov (1978, 1979),
its generalization to Einstein–Maxwell equations (Alekseev, 1981), and
the Bäcklund transformations (Harrison, 1978, 1980; Neugebauer, 1979).
With these techniques, most of the well-known solutions were rediscov-
ered, including the Kerr–Newman–NUT black hole solutions (Belinsky and
Verdaguer, 2001).
In fact, such techniques can be applied to all the spacetimes with two
orthogonal Killing vectors, including stationary axially symmetric space-
times (Griffiths and Podolský, 2009; Stephani et al., 2009), cylindrical
spacetimes (Verdaguer, 1993; Bronnikov, Santos and Wang, 2019), and
spacetimes with plane symmetry, studied in this book.
171
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch07 page 172
With the studies of exact solutions that represent colliding plane wave
spacetimes, now it is clear that spacetime singularities are generically devel-
oped due to the nonlinear interaction of the two colliding plane waves,
as shown explicitly in Chapters 4 and 5. Killing–Cauchy horizons can be
formed in particular cases (Chandrasekhar and Xanthopoulos, 1986a), but
they are not stable against small perturbations (Yurtsever, 1987; Clarke
and Hayward, 1989; Griffiths, 2005), and are expected to be turned into
real spacetime singularities.
One might think that this is due to the high symmetry of the spacetimes,
and in particular because of the plane symmetry, the incoming waves always
have infinitely large amount of energy. However, recently it was showed
that, due to such nonlinear interactions, spacetime singularities can be still
formed in the head-on collision of axisymmetric distributions of null parti-
cles even with finite energy (Pretorius and East, 2018).
This result is extremely encouraging, and motivating various interest-
ing questions. One is regarding to the internal structure of black holes. As
shown in Chapter 6, the interiors of the most well-known black holes are iso-
metric to the interacting regions of two colliding plane waves. Is it possible
to shed lights on the internal structure of black holes by studying the head-
on collision of two gravitational and/or matter wave? What about cosmic
censorship and hoop conjectures (Penrose, 1969; Thorne, 1972)? Can we
also be able to say something about the critical phenomena of gravitational
collapse in the threshold of black hole formation (Choptuik, 1993; see also
the review articles, Wang, 2001; Gundlach and Martin-Garcia, 2007).
To be more realistic, one might like first to remove the assumption
of uniform plane-fronted waves, so one can deal with the case in which
incoming waves have finite energy, similar to what were done by Pretorius
and East (2018), and more recently by Baumgarte, Gundlach and Hilditch
(2019). Perturbative calculations of this kind already started in the early
of 1990s by D’Eath and Payne (1992a, 1992b, 1992c), and such obtained
resulted are quite consistent with these numerical ones (Pretorius and East,
2018). Recently, such studies have been also generalized to other cases,
including that the background is an expanding universe (Centrella and
Matzner, 1982; Alekseev, 2016).
The above works have been mainly related to some theoretical aspects
of Einstein’s theory. With entering the era of gravitational wave astron-
omy, one might like also to ask more observational ones, such as the detec-
tion of gravitational memory effects (Favata, 2010; Bieri, Garfinkle and
Yunes, 2017), just briefly mentioned in Section 2.5, in addition to the
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-ch07 page 173
Bibliography
175
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-bibliography page 176
Bibliography 177
Belinsky, V.A. and Zakharov, V.E. (1978). Integration of the Einstein equations
by means of the inverse scattering problem technique and construction of
exact soliton solutions, Sov. Phys. JETP 48, 985.
Belinsky, V.A. and Zakharov, V.E. (1979). Stationary gravitational solitons with
axial symmetry, Sov. Phys. JETP 50, 1.
Bell, P. and Szekeres, P. (1974). Interacting electromagnetic shock waves in gen-
eral relativity, Gen. Relativ. Grav. 5, 275.
Bennett, C.L. Banday, A. Gorski, K.M. Hinshaw, G. Jackson, P. Keegstra,
P. Kogut, A. Smoot, G.F. Wilkinson, D.T. and Wright, E.L. (1996). Four-
year COBE DMR cosmic microwave background observations: Maps and
basic results, Astrophys. J. 464, L1.
Berezin, V.A. Kuzmin, V.A. and Tkachev, I. (1987). Dynamics of bubbles in
general relativity, Phys. Rev. D 36, 2919.
Bieri, L. Garfinkle, D. and Yunes, N. (2017). Gravitational waves and their math-
ematics, AMS Notices 64, 07.
Blanchet, L. and Damour, T. (1992). Hereditary effects in gravitational radiation,
Phys. Rev. D 46, 4304.
Bondi, H. (1957). Plane gravitational waves in general relativity, Nature (London)
179, 1072.
Bondi, H. and Pirani, F.A.E. (1989). Gravitational waves in general relativity
XIII. caustic property of plane waves, Proc. R Soc. Lond. A 421, 395.
Bondi, H. Pirani, F.A.E. and Robinson, I. (1959). Gravitational waves in general
relativity III. Exact plane waves, Proc. R Soc. Lond. A 251, 519.
Bonnor, W.B. (1983). Physical interpretation of vacuum solutions of Einstein’s
equations. Part I. Time-independent solutions, Gen. Relativ. Grav. 15, 535.
Bonnor, W.B. and Vichers, P.A. ( 1981). Junction conditions in general relativity,
Gen. Relativ. Grav. 13, 29.
Boulware, D.G. (1973). Naked singularities, thin shells, and the Reissner-
Norstrom metric, Phys. Rev. D 8, 2363.
Braginsky, V.P. and Grishchuk, L.P. (1985). Kinematic resonance and memory
effect in free-mass gravitational antennas, Zh. Eksp. Teor. Fiz. 89, 744 [Sov.
Phys. JETP 62, 427].
Braginsky, V.P. and Thorne, K.P. (1987). Gravitational-wave bursts with memory
and experimental prospects, Nature 327, 123.
Brandenberger, R.H. (1985). Quantum field theory methods and inflationary uni-
verse models, Rev. Mod. Phys. 57, 1.
Braun, V. and Ovrut, B.A. (2006). Stabilizing moduli with a positive cosmological
constant in heterotic M-theory, J. High Energy Phys., 07, 035.
Bretón, N. Garcı́a, A. Macı́as, A. and Yáñez, G. (1998). Colliding plane waves in
terms of Jacobi functions, J. Math. Phys. 39, 6051.
Bretón, N. Matos, T. and Garcı́a, A. (1996). Colliding plane waves in Einstein–
Maxwell-dilaton fields, Phys. Rev. D 53, 1868.
Brinkman. H.W. (1923). On riemann spaces conformal to einstein spaces, Proc.
Natl. Acad. Sci. USA 9, 1.
Bronnikov, K.A. (1980). Gravitational and sound waves in stiff matter, J. Phys.
A: Math. Gen. 13, 3455.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-bibliography page 178
Bronnikov, K.A. Santos, N.O. and Wang, A. (2019). Cylindrical systems in general
relativity, preprint, arXiv:1901.06561 [Class. Quantum Grav. in press (2020)].
Bruckman, W. (1986). Stationary axially symmetric exterior solutions in the five-
dimensional representation of the Brans–Dicke–Jordan theory of gravitation,
Phys. Rev. D 34, 2990.
Bruckman, W. (1987). Exact axially symmetric stationary solutions of the
Kaluza–Klein–Jordan–Thiry theory, Phys. Rev. D 36, 3674.
Campbell, S.J. and Wainwright, J. (1977). Algebraic computing and the Newman-
Penrose formalism in general relativity, Gen. Relativ. Grav. 8, 987.
Carot, J. and Verdaguer, E. (1989). Generalised soliton solutions of the Weyl
class, Class. Quantum Grav. 6, 845.
Carr, B.J. and Verdaguer, E. (1984). Soliton solutions and cosmological gravita-
tional waves, Phys. Rev. D 28, 2995.
Carroll, S. (2004). Spacetime and Geometry: An Introduction to General Relativity
(Addison-Wesley, New York).
Centrella, J. and Matzner, R.A. (1982). Colliding gravitational waves in expand-
ing cosmologies, Phys. Rev. D 25, 930.
Cespedes, J. and Verdaguer, E. (1987). Gravitational wave pulse and soliton wave
collision, Phys. Rev. D 36, 2259.
Challifour, J.L. (1972). Generalized Functions and Fourier Analysis (W. A.
Benjamin, Reading, MA).
Chan, R. da Silva, M.F.A. da Rocha, J.F. and Wang, A. (2011). Radiating gravas-
tars, J. Cosmol. Astropart. Phys. 10, 013.
Chandrasekhar, S. (1983). The Mathematical Theory of Black Holes (Clarendon
Press. Oxford).
Chandrasekhar, S. (1986). Cylindrical waves in general relativity, Proc. R. Soc.
Lond. A 408, 209.
Chandrasekhar, S. and Ferrari, V. (1984). On the Nutku–Halil solution for col-
liding impulsive gravitational waves, Proc. R. Soc. Lond. A 396, 55.
Chandrasekhar, S. and Xanthopoulos, B.C. (1985a). On colliding waves in the
Einstein-Maxwell theory, Proc. R. Soc. Lond. A 398, 223.
Chandrasekhar, S. and Xanthopoulos, B.C. (1985b). On the collision of impulsive
gravitational waves when coupled with fluid motions, Proc. R. Soc. Lond. A
402, 37.
Chandrasekhar, S. and Xanthopoulos, B.C. (1985c). Some exact solutions of grav-
itational waves when coupled with fluid motions, Proc. R. Soc. Lond. A 402,
205.
Chandrasekhar, S. and Xanthopoulos, B.C. (1986a). A new type of singularity
created by colliding gravitational waves, Proc. R. Soc. Lond. A 408, 175.
Chandrasekhar, S. and Xanthopoulos, B.C. (1986b). On the collision of impulsive
gravitational waves when coupled with null dust, Proc. R. Soc. Lond. A 403,
189.
Chandrasekhar, S. and Xanthopoulos, B.C. (1987a). The effect of sources on
horizons that may develop when plane gravitational waves collide, Proc. R.
Soc. Lond. A 414, 1.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-bibliography page 179
Bibliography 179
Bibliography 181
Feinstein, A. (2002). Penrose limits, the colliding plane-wave problem and the
classical string backgrounds, Class. Quantum Grav. 19, 5353.
Feinstein, A. and Charach, C. (1986). Gravitational solitons and spatial flatness,
Class. Quantum Grav. 3, L5.
Feinstein, A. and Ibañez, J. (1989). Curvature-singularity-free solutions for col-
liding plane gravitational waves with broken u − v symmetry, Phys. Rev. D
39, 470.
Feinstein, A. MacCallum, M.A.H. and Senovilla, J.M.M. (1989). On the ambigu-
ous evolution and the production of matter in space-times with colliding
waves, Class. Quantum Grav. 6, L217.
Feinstein, A. and Senovilla, J.M.M. (1989). Collision between variably polarized
plane gravitational wave and a shell of null matter, Phys. Lett. A 138, 102.
Ferrari, V. (1988). Focusing process in the collision of gravitational plane waves,
Phys. Rev. D 37, 3061.
Ferrari, V. and Ibañez, J. (1987a). On the collision of gravitational plane waves:
a class of soliton solutions, Gen. Relativ. Grav. 19, 405.
Ferrari, V. and Ibañez, J. (1987b). A new exact solution for colliding gravitational
plane waves, Gen. Relativ. Grav. 19, 383.
Ferrari, V. and Ibañez, J. (1988). Type-D solutions describing the collision of
plane-fronted gravitational waves, Proc. R. Soc. Lond. A 417, 417.
Ferrari, V. and Ibañez, J. (1989a). On the gravitational interaction of null fields,
Class. Quantum Grav. 6, 1805.
Ferrari, V. and Ibañez, J. (1989b). Gravitational interaction of massless particles,
Phys. Lett. A 141, 233.
Ferrari, V. Ibañez, J. and Bruni, M. (1987a). Colliding gravitational waves with
non-collinear polarization: A class of soliton solutions, Phys. Lett. A 122,
459.
Ferrari, V. Ibañez, J. and Bruni, M. (1987b). Colliding plane gravitational waves:
A class of non-diagonal soliton solutions, Phys. Rev. D 36, 1053.
Frolov, V.P. (1979). The Newman–Penrose method in the theory of general
relativity, in Problems in the General Theory of Relativity and Theory of
Group Representations, ed. Basov, N. G. and translated from Russian by
Mason A. (Proceedings (Trudy) of the P. N. Lebedev Physics Institute,
Vol. 96).
Garcá-Dı́az, A. (1988). Colliding-wave generalizations of the Nutku–Halil metric
in the Einstein–Maxwell theory, Phys. Rev. Lett. 61, 507.
Garcá-Dı́az, A. (1989). Colliding wave generalization of the Ferrari–Ibañez solu-
tion in the Einstein–Maxwell theory, Phys. Lett. A 138, 370.
Gardner, C.S. Greene, J.M. Kruskal, M.D. and Miura, R.M. (1967).
Method for solving the Korteweg-deVries equation, Phys. Rev. Lett. 19,
1095.
Gelfand, I.M. and Shilov, G.E. (1964). Generalized Functions. Vol. I: Properties
and Operations (Academic Press, New York, London).
Geroch, R. (1971). A method for generating solutions of Einstein’s equations,
J. Math. Phys. 12, 918.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-bibliography page 183
Bibliography 183
Bibliography 185
Hauser, I. and Ernst, F.J. (1980b). A homogeneous Hilbert problem for the
Kinnersely–Chitre transformations of electrovac space-times, J. Math. Phys.
21, 1418.
Hauser, I. and Ernst, F.J. (1981). Proof of a Geroch conjecture, J. Math. Phys.
22, 1051.
Hauser, I. and Ernst, F.J. (1989a). Initial value problem for colliding plane grav-
itational waves. I, J. Math. Phys. 30, 872.
Hauser, I. and Ernst, F.J. (1989b). Initial value problem for colliding plane grav-
itational waves. II, J. Math. Phys. 30, 2322.
Hauser, I. and Ernst, F.J. (1990). Initial value problem for colliding plane gravi-
tational waves. III, J. Math. Phys. 31, 871.
Hawking, S.W. (1975). Particles creation by black holes, Commun. Math. Phys.
43, 199.
Hawking, S.W. and Eillis, G.F.R. (1973). The Large Scale Structure of Space-
Time, Vol. 1 (Cambridge University Press, Cambridge).
Hawking, S. W. Perry, M.J. and Strominger, A. (2016). Soft hair on black holes,
Phys. Rev. Lett. 116, 231301.
Hawking, S. W. Perry, M. J. and Strominger, A. (2017). Superrotation
charge and supertranslation hair on Black Holes, J. High Energy Phys. 05,
161.
Helliwell, T.M. and Konkowski, D.A. (1990). Electromagnetic fields in Khan–
Penrose space-time, Phys. Rev. D 41, 2507.
Herrera, L. and Santos, N.O. (2005). Cylindrical collapse and gravitational waves,
Class. Quantum Grav. 22, 2407.
Hill, C.T. Schramm, D.N. and Fry, J.N. (1989). Cosmological structure formation
from soft topological defects, Comm. Nucl. Part. Phys. 19, 25.
Hirschmann, E.W. Wang, A. and Wu, Y. (2004). Collapse of a scalar field in 2 + 1
gravity, Class. Quantum Grav. 21, 1791.
Hoenselaers, C. and Dietz, W. (1984). Solutions of Einstein’s equations: tech-
niques and results, in Proceedings of the International Seminar on Exact
Solutions of Einstein’s Equations, in Retzbach, Germany, November 14–18,
1983 (Springer-Verlag).
Hoenselaers, C. and Ernst, F.J. (1990). Matching pp-waves to the Kerr metric,
J. Math. Phys. 31, 144.
Hoenselaers, C. Kinnersley, W. and Xanthopoulos, B.C. (1979a). Generation of
asymptotically flat space-times with any number of parameters, Phys. Rev.
Lett. 42, 481.
Hoenselaers, C. Kinnersley, W. and Xanthopoulos, B.C. (1979b). Symmetries
of the stationary Einstein-Maxwell equations VI-transformations which gen-
erate asymptotically flat space-times with arbitrary multiple moments,
J. Math. Phys. 20, 2530.
Hogan, P.A. Barrabés, C. and Bressange, G.F. (1998). Colliding plane waves in
Einstein–Maxwell theory, Lett. Math. Phys. 43, 263.
Hogan, P.A. and Walsh, D.M. (2003). Collision of high frequency plane gravita-
tional and electromagnetic waves, Int. J. Mod. Phys. 12, 1459.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-bibliography page 186
Holvorcem, P.R. Letelier, P.S. and Wang, A. (1995). The interaction of outgoing
and ingoing spherically symmetric null fluids, J. Math. Phys. 36, 3663.
Hu, N. (1983). Proceedings of the Third Marcel Grossmann Meeting on General
Relativity, Shanghai, China (Science Press, Beijing).
Horava, H. and Witten, E. (1995). Heterotic and type I string dynamics from
eleven dimensions, Nucl. Phys. B 460, 506.
Horava, H. and Witten, E. (1996). Eleven-dimensional supergravity on a manifold
with boundary, Nucl. Phys. B 475, 94.
Ibañez, J. and Verdaguer, E. (1983). Soliton collision in general relativity, Phys.
Rev. Lett. 51, 1313.
Ibañez, J. and Verdaguer, E. (1986). Finite perturbations on Friedmann–
Robertson–Walker models, Astrophys. J. 306, 401.
Ipser, J. and Sikivie, P. (1984). Gravitational repulsive domain wall, Phys. Rev.
D 30, 712.
Ishihara, H. Takahashi, M. and Tomimatsu, A. (1988). Gravitational Faraday
rotation induced by a Kerr block hole, Phys. Rev. D 38, 472.
Islam, J.N. (1985). Rotating Fields in General Relativity (Cambridge University
Press, Cambridge).
Israel, W. (1966). Singular hypersurfaces and thin shells in general relativity,
Il Nuovo Cimento B 44, 1.
Israel, W. (1967). Singular hypersurfaces and thin shells in general relativity,
Il Nuovo Cimento B 48, 463(E).
Ivanov, B.V. (1998). Colliding axisymmetric pp waves, Phys. Rev. D 57, 3378.
Jantzen, R.T. (1980). Soliton solutions of the Einstein equations generated from
cosmological solutions with additional symmetry, Nuovo Cimento B 59, 287.
Jogia, S. and Griffiths, J.B. (1980). A Newman–Penrose formalism for space-tunes
with torsion, Gen. Relativ. Grav. 12, 597.
Jordan, P. Ehlers, J. and Kundt, W. (1960). Akad. Wiss. Matnz. Abn. Math.
-Nat. Kl., Jabrg 2.
Kachru, S. Schulz, M.B. and Trivedi, S. (2003). Moduli stabilization from fluxes
in a simple IIB orientifold, J. High Energy Phys. 10, 007.
Kerr, R.P. (1963). Gravitational field of a spinning mass as an example of alge-
braically special metric, Phys. Lett. 11, 237.
Khan, K.A. and Penrose, R. (1971). Scattering of two impulsive gravitational
plane waves, Nature (London) 229, 185.
Khorrami, M. and Mansouri, R. (1994). Cylindrically symmetric thin walls in
general relativity, J. Math. Phys. 35, 951.
Kibble, T.W.B. (1976). Topology of cosmic domains and strings, J. Phys. A:
Math. Gen. 9, 1387.
Kinnersley, W. (1969). Type D vacuum metrics, J. Math. Phys. 10,
1195.
Kinnersley, W. and Walker, M. (1970). Uniformly accelerating charged mass in
general relativity, Phys. Rev. D 2, 1359.
Kitchingham, D.W. (1984). The use of generating techniques for space-times with
two non-null commuting Killing vectors in vacuum and stiff perfect fluid
cosmological models, Class. Quantum Grav. 1, 677.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-bibliography page 187
Bibliography 187
Bibliography 189
Nolan, B.C. (2000). Central singularity in spherical collapse, Phys. Rev. D 62,
044015.
Nutku, Y. (1981). Comment on “collision of plane gravitational waves without
singularities”, Phys. Rev. D 24, 1040.
Nutku, Y. and Halli, M. (1977). Colliding impulsive gravitational waves, Phys.
Rev. Lett. 39, 1379.
O’Brien, S. and Synge, J.E. (1952). Jump conditions at discontinuities in general
relativity, Comm. Dublin Inst. Adv. Stud. A9.
Oliver, G. and Verdaguer, E. (1989). A family of inhomogeneous cosmological
Einstein–Rosen metrics, J. Math. Phys. 30, 442.
O’Raifeartaigh, L. (1972). General Relativity, Papers in Honour of J.L. Synge
(Clarendon Press, Oxford).
Ori, A. (2000). Strength of curvature singularities, Phys. Rev. D 61, 064016.
Paiva, F.M. and Wang, A. (1995). Geometry of planar domain walls, Phys. Rev.
D 52, 1281.
Palenta, S. and Meinel, M. (2017). A continuous Riemann-Hilbert problem for
colliding plane gravitational waves, Class. Quantum Grav. 34, 195011.
Pantaleo, X. (1979). Relativity, Quanta and Cosmology, Einstein 1879–1979 (New
York, Johnson).
Papacostas, T. and Xanthopoulos, B.C. (1989). Collisions of gravitational and
electromagnetic waves that do not develop curvature singularities, J. Math.
Phys. 30, 97.
Papapetrou, A. and Hamoui, A. (1968). Couches simples de mattere en relativite
generale, Ann. Inst. Henri Poincare IX, 179.
Penrose, R. (1960). A spinor approach to general relativity, Ann. Phys. 10, 171.
Penrose, R. (l963). Asymptotic properties of fields and space-times, Phys. Rev.
Lett. 10, 66.
Penorse, R. (1964). Conformal treatment of infinity, in Relativity, Groups and
Topology, The 1963 Les Houches Lectures, eds. DeWitt, C. and DeWitt, B.
(Gordon and Breach, New York), pp. 565–584.
Penrose, R. (1965). A remarkable property of plane waves in general relativity,
Rev. Mod. Phys. 37, 215.
Penrose, R. (1968). Structure of space-time, in Battelle Rencontres, Lectures in
Mathematical Physics, eds. DeWitt, C. and Wheeler, J.A. (Benjamin, New
York).
Penrose, R. (1969). Gravitational collapse: the role of general relativity, Riv.
Nuovo Cimento, 1, 252.
Penrose, R. (1972). The geometry of impulsive gravitational waves, in General
Relativity, Papers in Honor of J.L Synge, ed. O’Raifeartaigh, L. (Oxford
University Press, Oxford), pp. 101–115.
Pereira, P.R.C.T. and Wang, A. (2000a). Dynamic cylindrical shells with rotation
in general relativity, Gen. Relativ. Grav. 32, 2189.
Pereira, P.R.C.T. and Wang, A. (2000b). Gravitational collapse of counter-
rotating cylindrical dust shells, Phys. Rev. D 62, 124001.
Pereira, P.R.C.T. and Wang, A. (2002). Gravitational collapse of counter-rotating
cylindrical dust shells, Int. J. Mod. Phys. D 11, 561.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-bibliography page 190
Petrov, A.Z. (1955). New Methods in General Relativity (Nauka, Moscow) (in
Russian) [English edition of Petrov’s book: Einstein Spaces, Pergamon Press,
1969].
Pirani, F. (1964). Introduction to gravitational radiation theory, in Lectures in
General Relativity, Brandeis Summer Institute in Theoretical Physics, Vol. 1.
eds. Deser, S. and Ford, K. (Prentice-Hall, Englewood Cliffs).
Piran, T. and Safier, P.N. (1985). A gravitational analogue of Faraday rotation,
Nature (London) 318, 271.
Piran, T. Safier, P.N. and Katz, J. (1986). Cylindrical gravitational waves with
two degrees of freedom: an exact solution, Phys. Rev. D 34, 331.
Piran, T. Safier, P.N. and Stark, R.F. (1985). General numerical solution of cylin-
drical gravitational waves, Phys. Rev. D 32, 3101.
Pretorius, F. and East, W.E. (2018). Black hole formation from the collision of
plane-fronted gravitational waves, Phys. Rev. D 98, 084053.
Randall, L. and Sundrum, R. (1999a). Large mass hierarchy from a small extra
dimension, Phys. Rev. Lett. 83, 3370.
Randall, L. and Sundrum, R. (1999b). An alternative to compactification, Phys.
Rev. Lett. 83, 4690.
Ringeval, C. and Suyama, T. (2017). Stochastic gravitational waves from cosmic
string loops in scaling, J. Cosmol. Astropart. Phys. 12, 027.
Rocha, P. Chan, R. da Silva, M.F.A. and Wang, A. (2008a). Stable and “bounded
excursion” gravastars, and black holes in Einstein’s theory of gravity, J. Cos-
mol. Astropart. Phys. 11, 010.
Rocha, P. Miguelote, A.Y. Chan, R. da Silva, M.F.A. Santos, N.O. and Wang,
A. (2008b). Bounded excursion stable gravastars and black holes, J. Cosmol.
Astropart. Phys. 06, 025.
Rosen, N. (1937). Plane polarised waves in the general theory of relativity, Phys.
Z 12, 366.
Rosenbaum, M. Ryan, M. Urrutia, L.F. and Matzner, R.A. (1986). Colliding plane
waves in N = 1 classical supergravity, Phys. Rev. D 34, 409.
Roy, S. (2003). Accelerating cosmologies from M /string theory compactifications,
Phys. Lett. B 567, 322.
Ruffini, R. (1976). Proceedings of the First Marcel Grossmann Meeting on General
Relativity, 7–12 July 1975 (North-Holland).
Sachs, R.K. (1960). Propagation laws for null and type III gravitational waves,
Phys. Z 157, 462.
Sachs, R.K. (1961). Gravitational waves in general relativity, IV: the outgoing
radiation condition, Proc. R. Soc. Lond. A 264, 309.
Sachs, R.K. (1962). Gravitational waves in general relativity, VIII: waves in
asymptotically flat space-time, Proc. R. Soc. Lond. A 270, 103.
Sachs, R.K. (1964). Gravitational radiation, in Relativity, Groups and Topology,
eds. DeWitt, C. and DeWitt, B.S. (Gordon and Breach, New York).
Sanz, J.L. and Goicoechea, L.J. (1984). Observational and Theoretical Aspects
of Relativistic Astrophysics and Cosmology (World Scientific Publishing Co.).
Sharma, P. Tziolas, A. Wang, A. and Wu, Z.-C. (2011). Spacetime Singularities
in string and its low dimensional effective theory, Int. J. Mod. Phys. A 25,
273.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-bibliography page 191
Bibliography 191
Schmidt, H.-J. and Wang, A. (1993). Plan domain walls when coupled with the
Brans-Dicke scalar field, Phys. Rev. D 47, 4425.
Schwarz, P. (1997). Colliding axion–dilaton plane waves from black holes, Phys.
Rev. D 56, 7833.
Siegel, H.P. (1981). Evolving dust shells, Phys. Rev. D 23, 2835.
Singh, P. and Griffiths, J.B. (1990). The application of spin coefficient techniques
in the vacuum quadratic Poincare gauge field theory, Gen. Relativ. Grav. 22,
269.
Souriau, J.-M. (1973). Ondes et radiations gravitationnelles, Colloq. Int. CNRS,
Paris 220, 243.
Spergel, D.N. et al. (WMAP) (2007). Wilkinson microwave anisotropy probe
(WMAP) three year results: implications for cosmology, Astrophys. J. Suppl.
170, 377.
Stark, R.F. and Connors, P.A. (1977), Observational test for the existence of a
rotating black hole in CygX-1, Nature (London) 266, 429.
Stephani, H. Kramer, D. MacCallum, M. Hoenselaers, C. and Herlt, E. (2009).
Exact Solutions of Einstein’s Field Equations, Cambridge Monographs on
Mathematical Physics (Cambridge University Press, Cambridge).
Stoyanov, Yu.G. (1979). Interaction of plane gravitational waves, Sov. Phys.
JETP. 44, 1051.
Strominger, A. (2017). Lectures on the Infrared structure of gravity and gauge
theory, preprint, arXiv:1703.05448.
Synge, J.L. (1965). Relativity: The General Theory (North-Holland).
Szekeres, P. (1965). The gravitational compass, J. Math. Phys. 6, 1387.
Szekeres, P. (1966). On the propagation of gravitational fields in matter, J. Math.
Phys. 7, 751.
Szekeres, P. (1970). Colliding gravitational waves, Nature (London) 228, 1183.
Szekeres, P. (1972). Colliding plane gravitational waves, J. Math. Phys. 13, 286.
Tabensky, R. and Taub, A.H. (1973). Plane symmetric self-gravitating fluids with
pressure equal to energy density, Commun. Math. Phys. 29, 61.
Talbot, C.J. (1969). Newman–Penrose approach to twisting degenerate metrics,
Commun. Math. Phys. 13, 45.
Taub, A.H. (1956). Isentropic hydrodynamics in plane symmetric space-times,
Phys. Rev. 103, 454.
Taub, A.H. (1975). Spatially homogenous universe, in Proceedings of the First
Marcel Grossmann Meeting on General Relativity, ed. Ruffini, R. (North-
Holland).
Taub, A.H. (1980). Space-times with distribution valued curvature tensors,
J. Math. Phys. 21, 1423.
Taub, A.H. (1983). General relativistic hydrodynamic wave propagation, in Pro-
ceedings of the Third Marcel Grossmann Meeting On General Relativity,
August 30–September 1, 1982, Shanghai, China, ed. Hu, N. (Science Press,
Beijing).
Taub, A.H. (1988a). On the collision of planar impulsive gravitational waves,
J. Math. Phys. 29, 690.
Taub, A.H. (1988b). Collision of impulsive gravitational waves following dust
clouds, J. Math. Phys. 29, 690.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-bibliography page 192
Taub, A.H. (1990). Interaction of null dust clouds fronted by impulsive plane
gravitational waves, J. Math. Phys. 31, 664.
Thorne, K.S. (1967). The general relativistic theory of stellar structure and
dynamics, in High-Energy Astrophysics, vol. Ill, eds. DeWitt, C. Schatzman,
E. and Veron, P. (Gordon and Breach, New York), p. 414.
Thorne, K.S. (1972). In Magic without Magic: John Archibald Wheeler, ed.
J. Klauder (Freeman, San Francisco).
Thorne, K.S. (1992). Gravitational-wave bursts with memory: the Christodoulou
effect, Phys. Rev. D 45, 520.
Tippler, F.J. (1980). Singularities from colliding plane gravitational waves, Phys.
Rev. D 22, 2929.
Tomita, K. (1998). Plane-symmetric vacuum solutions with null singularities for
inhomogeneous models and colliding gravitational waves, preprint, arXiv:gr-
qc/9807033.
Townsend, P.K. and Wohlfarth, M.N.R. (2003). Accelerating cosmologies from
compactification, Phys. Rev. Lett. 91, 061302.
Tsoubelis, D. (1976). Homogeneous relativistic cosmological models of Bianchi
type IV, Ph.D. dissertation, Graduate Faculty in Physics Department, The
City University of New York.
Tsoubelis, D. (1989a). Cylindrical shells of cosmic strings, Class. Quantum Grav.
6, 101.
Tsoubelis, D. (1989b). Plane gravitational waves colliding with shells of null dust,
Class. Quantum Grav. 6, L117.
Tsoubelis, D. and Wang, A. (1989). Asymmetric collision of gravitational plane
waves: a new class of exact solutions, Gen. Relativ. Grav. 21, 807.
Tsoubelis, D. and Wang, A. (1990). Impulsive shells of null dust colliding with
gravitational plane waves, Gen. Relativ. Grav. 22, 1091.
Tsoubelis, D. and Wang, A. (1991). On the gravitational interaction of plane
symmetric null dust clouds, J. Math. Phys. 31, 1017.
Tsoubelis, D. and Wang, A. (1992). Head-on collision of gravitational plane waves
with non-collinear polarization: a new class of analytical models, J. Math.
Phys. 33, 1054.
Tziolas, A. and Wang, A. (2008). Colliding branes and formation of spacetime
singularities, Phys. Lett. B 661, 5.
Tziolas, A. Wang, A. and Wu, Z.-C. (2009). Colliding branes and formation of
spacetime singularities in string theory, J. High Energy Phys. 04, 038.
Verdaguer, E. (1984). Solitons and the generation of new cosmological solutions, in
Observational and Theoretical Aspects of Relativistic Astrophysics and Cos-
mology, eds. Sanz, J.L. and Goicoechea, L.J. (World Scientific, Singapore).
Verdaguer, E. (1987). Plane waves and soliton solutions, Il Nuovo. Cimento B
100, 787.
Verdaguer, E. (1993). Soliton solutions in spacetimes with two spacelike killing
fields, Phys. Rep. 229, 1.
Weyl, H. (1917). Zur gravitationstheorie, Ann. Phys. 54, 117.
Vladimirov, V.S. (1979). Generalized Functions in Mathematical Physics (MIR,
Moscow).
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-bibliography page 193
Bibliography 193
Vilenkin, A. (1981). Gravitational field of vacuum domain walls and strings, Phys.
Rev. D 23, 852.
Vilenkin, A. (1985). Cosmic strings and domain walls, Phys. Rep. 121, 263.
Vilenkin, A. and Shellard, E.P.S. (2000). Cosmic Strings and Other Topological
Defects (Cambridge University Press, Cambridge).
Wainwright, J. Ince, W.C.W. and Marshman, B.J. (1979). Spatially homogeneous
and inhomogeneous cosmologies with equation of state, Gen. Relativ. Grav.
10, 259.
Wald, R.M. (1984). General Relativity (University of Chicago, New York).
Wands, D. (2006). Brane-world cosmology, preprint, arXiv:gr-qc/0601078.
Wang, A. (1991a). The effect of polarization of colliding plane gravitational waves
on focusing singularities, J. Mod. Phys. A 6, 2273.
Wang, A. (1991b). Gravitational Faraday rotation induced from interacting grav-
itational plane waves, Phys. Rev. D 44, 1120.
Wang, A. (1991c). Plane walls interacting with gravitational and matter fields,
J. Math. Phys. 32, 274.
Wang, A. (1991d). Massless particle radiation of Vilenkin’s planar domain walls,
Phys. Lett. B 264, 2863.
Wang, A. (1991e). Planar domain walls emitting and absorbing electromagnetic
radiation, Phys. Rev. D 44, 1705.
Wang, A. (1991f). Interacting Gravitational, Electromagnetic, Neutrino and Other
Waves in the Context of Einstein’s General Theory of Relativity (Iaonnina
University, Greece).
Wang, A. (1992a). Gravitational interaction of plane gravitational waves and
matter shells, J. Math. Phys. 33, 1065.
Wang, A. (1992b). Planar domain walls coupled with pure radiation fields, Int.
J. Mod. Phys. A 7, 4521.
Wang, A. (1992c). Singularities of the space-time of a plane domain wall when
coupled with a massless scalar field, Phys. Lett. B 277, 49.
Wang, A. (1992d). Transparency of domain walls to electromagnetic and neutrino
waves: an exact solution, Mod. Phys. Lett. A 7, 835.
Wang, A. (1992e). Dynamics of plane-symmetric thin walls in general relativity,
Phys. Rev. D 45, 3534.
Wang, A. (1992f). On the interaction of bubbles with gravitational and matter
fields, Mod. Phys. Lett. A 7, 1779.
Wang, A. (1993). Non-trivial interaction of plane domain walls with scalar fields,
Phys. Rev. D 48, 2591.
Wang, A. (1994). Gravitational collapse of thick domain walls: an analytic model,
Mod. Phys. Lett. A 9, 3605.
Wang, A. (2000). On “Comments on collapsing null shells of matter”, J. Math.
Phys. 41, 8354.
Wang, A. (2001). Critical phenomena in gravitational collapse: The studies so far,
Braz. J. Phys. 31, 188.
Wang A. (2003). Critical collapse of cylindrically symmetric scalar field in four-
dimensional Einstein’s theory of gravity, Phys. Rev. D 68, 064006.
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-bibliography page 194
Bibliography 195
Xanthopoulos, B.C. (1986b). Cylindrical waves and cosmic strings of Petrov type
D, Phys. Rev. D 34, 3608.
Xanthopoulos, B.C. (1986c). The initial value problem for colliding gravitational
and hydrodynamic waves, J. Math. Phys. 27, 2129.
Xanthopoulos, B.C. (1987). Cosmic strings coupled with gravitational and elec-
tromagnetic waves, Phys. Rev. D 35, 3713.
Yurtsever, U. (1987). Instability of Killing-Cauchy horizons in plane-symmetric
space-times, Phys. Rev. D 36, 1662.
Yurtsever, U. (1988a). Structure of singularities produced by colliding plane grav-
itational waves, Phys. Rev. D 38, 1706.
Yurtsever, U. (1988b). New family of exact solutions for colliding plane gravita-
tional waves, Phys. Rev. D 37, 2790.
Yurtsever, U. (1989). Singularities and horizons in the collisions of gravitational
waves, Phys. Rev. D 40, 329.
Zel’dovich, Ya.B. (1980). Cosmological fluctuations produced near a singularity,
Mon. Not. R. Astron. Soc. 192, 663.
Zel’dovich, Ya.B. Kobzarev, [Link]. and Okun, L.B. (1976). Cosmological conse-
quences of a spontaneous breakdown of a discrete symmetry, Sov. Phys.
JETP 40, 1.
Zeldovitch, Ya.B. and Polnarev, A.G. (1974). Radiation of gravitational waves by
a cluster of superdense stars, Sov. Astron. 18, 17.
Zemanian, A.H. (1987). Distribution Theory and Transform Analysis: An Intro-
duction to Generalized Functions, with Applications (Dover Publications,
New York).
Zhang, P.-M. Duval, C. Gibbons, G.W. and Horvathy, P.A. (2017). The memory
effect for plane gravitational waves, Phys. Lett. B 772, 743.
Zhang, P.-M. Duval, C. Gibbons, G.W. and Horvathy, P.A. (2018). Velocity mem-
ory effect for polarized gravitational waves, J. Cosmol. Astropart. Phys. 05,
030.
b2530 International Strategic Relations and China’s National Security: World at the Crossroads
Index
197
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-index page 198
constantly polarized, 48, 91, 111, 112, electromagnetic field, 25, 79, 80, 82,
119, 131, 134, 138–140, 142 83, 85, 156
continuity, 33 electromagnetic plane waves, 129,
coordinate singularities, 49, 50 160, 161, 169
coordinate transformation, 4, 44, 46, Ellis and Schmidt, 49
49, 165 energy conditions, 27, 146, 169
cosmic acceleration, 30 energy density, 5, 27–29, 77, 145, 151,
cosmic censorship, 172 155
cosmic strings, 29 energy–momentum tensor, 2, 18, 28,
cosmic superstrings, 29 158
cosmological constant, 2, 5, 17 equation of state, 27, 84, 145, 149,
Coulomb field, 66, 78 150, 153
critical phenomena, 172 Ernst equation, 85, 92, 93, 163
curvature singularities, 49, 50, 81, 83, Ernst potential, 68, 85, 92, 124, 127,
96, 101, 107, 108, 122, 130, 161 166
cylindrical Ernst equation, 163 event horizon, 168
cylindrical spacetimes, 163, 171 expansion, 24, 28, 41, 96, 102
Index 199
neutrino field, 25–27, 80, 157 Ricci scalars, 25, 26, 45, 48, 53, 63,
Newman and Penrose, 1, 12, 15, 25 70, 76–78, 80, 92, 146
non-collinearly, 56, 73, 91, 102, 103, Rosen, 49
105, 106, 119, 129, 161
nonlinear interaction, 76, 77, 79, 142, Sachs, 40
171–173 sandwich waves, 55
NP formalism, 1, 12, 14, 16 Schmidt and Wang, 30, 60
null field, 39–41 Schwarzschild black hole, 163–165
null shell, 36 second fundamental form, 28, 31
null tetrad, 44, 56, 157 seed solutions, 103, 107, 116, 119
NUT parameter, 123, 124, 168 Sephani et al., 169
Nutku and Halil, 66, 81, 166 shear, 40, 41, 56, 57
Nutku–Halil solution, 102, 109, 110, singular surface, 28, 30, 31, 49, 60
112, 115, 123, 124, 166 soft graviton theorems, 49, 51
soliton solutions, 103
O’Brien and Synge, 28, 29 soliton structure, 95, 96
one-soliton solutions, 125, 138 soliton technique, 67
optical scalars, 1, 24 space-like singularity, 82, 83
Ori, 50, 51 spin coefficients, 13–16, 21, 26, 27, 45,
56, 62, 80
Papacostas and Xanthopoulos, 161, spinor, 26, 30, 80
169 Stark and Connors, 70
parallel transport, 1, 24 Stephani et al., 1, 10, 12, 25, 50, 83,
parallelly transported, 6, 8, 47, 53, 161, 168, 171
71, 72, 75, 76, 115 stochastic, 29
Penrose, 12, 30, 42, 172 string network, 29
Pereira and Wang, 30, 31 Strominger, 48, 52
perturbations, 29, 81, 84, 101, 168, strong energy condition, 28
172 Szekeres, 44, 66, 77, 80, 81, 95, 130,
phase transitions, 29, 129, 154 134
Piran and Safier, 71 Szekeres family, 127, 134
Pirani, 5, 12 Szekeres solution, 95, 127
plane-fronted, 41, 43, 162, 169, 172
polarization angles, 73, 110, 115 Taub, 63, 145, 146, 171
polarization direction, 56, 72, 76, 116, Taub–NUT solution, 125, 169
145 tetrad components, 9, 15, 17
polarizations, 71, 76, 77, 81, 82, 91, tetrad frame, 14, 15
93, 109, 114, 142, 143, 166, 170 Thorne, 5, 48, 172
pressure, 27, 145, 151 tidal forces, 50
Pretorius and East, 162, 172 topological defects, 29
pure radiation field, 25, 39, 40, 77, 81 Tsoubelis, 63, 130, 134, 135
Tsoubelis and Wang, 26, 63, 67, 68,
radiation field, 79 73, 77–83, 89, 94, 95, 97, 102, 103,
Randall and Sundrum, 30 126, 130, 145, 146, 148, 153, 163
real poles, 87 two-soliton solutions, 91, 116, 126,
Ricci identity, 4, 14 127
March 7, 2020 9:54 Interacting Gravitational, Electromagnetic. . . – 9in x 6in b3813-index page 201
Index 201