arXiv:0706.2511v1 [gr-qc] 17 Jun 2007
Chaos in Kundt type III Spacetimes
I.Sakalli∗ and M.Halilsoy†
Physics Department, Eastern Mediterranean University, G.Magosa, N.Cyprus, Mersin 10, Turkey.
(Date textdate; Received textdate; Revised textdate; Accepted textdate; Published textdate)
Abstract
We consider geodesics motion in a particular Kundt type III spacetime in which Einstein-Yang-Mills equations admit solutions. On a particular surface as constraint we project the geodesics into the (x, y) plane and treat the problem as a 2-dimensional one. Our numerical study shows that chaotic behavior emerges under reasonable conditions.
∗Electronic address: izzet.sakalli@emu.edu.tr
Kundt’s class of solutions present spacetimes with non-expanding, shear-free and twist-free null geodesic congruences [1]. Interest in this class of solutions comes from the fact that these spacetimes admit plane waves, which exhibit geometrically different properties than the pp-waves. Due to its potential application in string theory Kundt class still maintains its popularity [2, 3, 4, 5, 6]. Various Petrov types of Kundt solutions were identified, among them especially Kundt type N has been studied in detail for long time by many authors,(see for instance [7, 8, 9]). Geodesics in a specific Kundt type N was analyzed in detail by [9] , and it was shown that particular solutions obey a power-law. Recently, interest in the Kundt type III spaces has gained momentum. Firstly, Griffiths et al [10] derived and classified a complete family of Kundt type III, which admit cosmological constant Λc and/or
pure radiation Φ22 6= 0. In the absence of the cosmological constant Λc, the Kundt type
III solutions are further generalized to be the solutions of the Einstein-Yang-Mills (EYM) system by Fuster and Holten [11].
The pioneering study of proving chaos in the spacetimes of plane waves was done for pure impulsive gravitational pp-waves [12] . Recently, it has been further demonstrated that under certain conditions the emergence of chaotic motion is possible both in the spacetimes of the superposed electrovac pp-waves and in the non-Abelian plane waves of Kundt type N [13], which are the solutions of the D = 4 EYM equations. In [13], it was also pointed out that for the non-Abelian plane waves the chaotic effect of gravity dominates over the gauge field. This is due to the fact that such local fields vanish asymptotically and the chaos inherited from gravity renders the whole system chaotic.
Analysis of the geodesic motion for the D = 4 Kundt type III is also discussed by [11] in which a possible chaotic motion is highlighted in particular cases of this class.
In this Brief Report, our aim is to investigate the possibility of chaotic geodesics in Kundt class III spacetimes. For this purpose, we consider the following algebraically special line-element [3, 11],
ds2 = 2du[dv + Hdu + W dz + W dz] − 2dzdz, (1)
where H = H(u, v, z, z) is a real function while W = W (u, v, z, z) is a complex function, in general. Here our motivation is to study the geodesic motion for the case W,v = z+z−2,
W = W0 (u, z) − z2v+ z, H = H0 (u, z, z) + W 0 (u, z) + W0(u, z) z+ z v− v2 (z + z)2. (2) where W0
(u, z) is an arbitrary complex function and H0
(u, z, z) is a real function. The simplest choice of W0
(u, z) for the type III is
W0
(u, z) = g(u)z, (3)
such that the Weyl scalar Ψ3 6= 0. On the other hand, imposing the only solution (χa =
λa(u)z) on the Yang-Mills (YM) equation in which the energy density is bounded throughout
the spacetime, the solution for H0
(u, z, z) becomes
H0(u, z, z) =f (u, z) + f(u, z) (z + z) − ggzz + σ(u)(z + z)2
[ln(z + z) − 1] . (4) where f (u, z) is an arbitrary complex function and σ(u) = 2γabλa(u)λ
b
(u) + g(u)g(u) is a real function. Here, λa(u) are bounded complex functions and γ
ab is the invariant metric of
the Lie group. We note that the condition for being Kundt type III spacetimes is g(u) 6= 0. It is trivially seen that λa(u) = 0 corresponds to the vacuum solution.
Our primary interest here is to write the geodesics equations for the metric (1). Similar to the study [11], for the beginning, we eliminate the u-dependence from W0
and g by the following particular choice
W0 = z and g = 1, (5)
Next, introducing real spatial coordinates x and y by√2z = (x + iy), we get the geodesic equations as ·· u−u·2(1 − v x2) + 2 · u · x x = 0, (6) ·· x+u· 2h H,x− 1 − xv2 x− 2v x i + 2 x · uv· + 2u·x· 1 − 2 v x2 = 0, (7) ·· y+u·2hH,y+ y 1 − xv2 i − 2u·x·y x = 0, (8)
where the dot denotes d
dτ with τ being the proper time. In addition the metric condition
implies · x2+y·2− 2u·v· − 2Hu·2− 2u·x·(x − 2v x) + 2y · uy· = ǫ. (9)
where ǫ = 1, 0, −1 stands for timelike, null and spacelike geodesics, respectively. The present form of the equation set does not allow us to obtain a 2D (x, y) Hamiltonian system analogous to the previous studies [12, 13]. However, with appropriate choices of u and v surfaces it is possible to project the geodesics into the (x, y) plane in which writing a 2D Hamiltonian becomes possible. Our first intention is to shift the independent variable from τ to u as an affine parameter. If we consider a family of geodesics, which follow the light-cone coordinate u with constant rate of change in the same proper time intervals, the following assumption can be made
·
u= constant ≡ 1, (10)
Such an assumption gives rise to a condition on the v surfaces given by
v = x2
(1 − 2x′
x), (11)
Here ”′ ” denotes d
du. By this substitution into Eqs. (7) and (8), we get a 2D dynamical
system in the (x, y) plane
3x′′− H0
,x− x = 0, (12)
y′′+ H0
,y = 0, (13)
which is described by a Super-Hamiltonian [14]
H =12 Py2− P 2 x 3 + V (x, y), (14)
with the corresponding potential
V(x, y) = H0
+x
2
Let us note that the Super-Hamiltonian defined by the momenta Px= −3x′ and Py = y′
is not positive definite.
Eq. (9) stands for an energy condition in which it should be automatically satisfied by the solutions of Eqs. (12) and (13). Without loss of generality, we can assume that f and λa are independent of u. This assumption implies that σ is a positive constant. As we
mentioned before that the chaotic effect of gravity dominates over the gauge (YM) field [13], it would be sufficient to investigate chaos in vacuum, i.e. σ = 1. In other words, once the chaotic motion appears in the vacuum spacetime, the local fields could not be strong enough to avert it into a regular motion.
In general, any f = kzn (n = 0, 1, 2...), with the multiplicative factor k being an arbitrary
real parameter implies a potential
V =√2kx Re(zn ) −y 2 2 + 2x 2h ln(√2x) − 1i. (16)
which admits an integrable system for k = 0. The logarithmic term in the potential imposes a condition on the x coordinate, namely x > 0. Beside this, the case f = kz, which describes a flat space for vacuum pp-wave spacetimes [1, 12], confesses a regular motion for the geodesics particles. Contrary to the vacuum homogeneous pp-waves [12] , here the case n= 2 admits a nonintegrable dynamical system. However, the nonintegrable systems in the vacuum pp-waves were obtained for the cases with n ≥ 3, [12] .
Here we wish to study the nonintegrable system with the simplest case (n = 2), and explore whether the motion depends on the initial conditions or not. If chaos emerges in such a simplest case, intuitively it should also appear for the f functions with higher powers of z. Now, for n = 2, it can be seen that the potential (16) has various unstable points according to the range of the multiplicative constant k.
TABLE I: Unstable points depending on the value of k.
Points Saddle Repellor
n x= −1 2√2k, y = ± 1 2√2kp7 + 8 ln(−2k) o k <−12e−78 — n x= 2√2 3k LambertW 3 4k √e , y = 0o k >−12e−78 −4 3e −3 2 < k <−1 2e −7 8
Those results in the table show us that the highest possibility of the emergence of chaos corresponds to the case k < −1
2e−
7
8 , in which admits two saddle points. Particularly,
the case −4 3e −3 2 < k < −1 2e −7
stronger chaos. Conversely, the case k > −1 2e−
7
8 upon possessing one saddle point causes a
questionable chaotic motion. In order to judge the existence of the chaotic motion, we study the numerical analysis of the evolution of the test particle in the gravitational field.
We integrate numerically the equations of motion given by Eqs. (12) and (13). The initial conditions depend on 3 parameters, (x0, y0) (at u = 0) and k. For a given k value,
we may choose (x0, y0) such that keeping x0 unchanged, and checking the effect of the y0 on
the geodesic motions, while it varies. To do this, we may set x0 = c1, a real constant c1 >0,
and y0 = −3 + 18
P
j=0 j
3. For example, if we take k = 0, the solutions are trivially analytic. This
is also graphically verified in Fig. 1. Next, by considering the cases k 6= 0 the motion can lead to a chaotic motion. It is observed that the chaos has a movable character depending on the choices of x0 and y0 while k > −
1 2e−
7
8, see Fig. 2. However, when the multiplicative
constant k < −1 2e−
7
8 chaos is obvious. It is seen that the multiplicative constant k of the
function f becomes decisive for the chaotic motion. In other words, k plays the role of critical parameter for the onset of chaos. The chaotic behavior of the geodesics in the case k <−1
2e−
7
8 is illustrated in Fig. 3. Alternatively, the dynamical system can be investigated
by using the Poincar´e section method. We use the package POINCAR´E [14] to perform the numerical experiments. Fig. (4) is a demonstration of the Poincar´e section, which verifies the chaotic behavior in our dynamical system.
In conclusion, it is shown that the Kundt type III spacetimes may reveal chaotic motion under certain conditions. To our knowledge such a study did not exist in the literature before. Chaos in the spacetimes of electrovac and the specific Kundt type N with YM field was studied before [13]. This report constitutes an extension of that study. It is needless to state that the existence of the chaos in the Kundt type III spacetimes may have further implications for the particle motions in string theory and in higher dimensions.
[1] H. Stephani, D. Kramer, M. MacCallum, C. Hoenselaers, and E. Herlt, Exact Solutions of Einstein’s Field Equations, 2nd ed. (Cambridge University Press, Cambridge, 2003).
[2] R. Guven and E. Y¨or¨uk, Phys. Rev. D 54, 6413 (1996). [3] A. A. Coley, Phys.Rev.Lett. 89, 281601 (2002).
67, 104020 (2003).
[5] A. Coley, S. Hervik, and N.Pelavas, Class. Quant. Grav. 23, 3053 (2006). [6] A. Coley, A. Fuster, S. Hervik, and N.Pelavas, hep-th/0703256.
[7] R. Guven, Phys. Rev. D 19, 471 (1979).
[8] J. Podolsky and M. Ortaggio, Class. Quant. Grav. 20, 1685 (2003). [9] J. Podolsky and M. Belan, Class. Quant. Grav. 21, 2811 (2004).
[10] J. B. Griffiths, P. Docherty, and J. Podolsky, Class. Quant. Grav. 21, 207 (2004). [11] A. Fuster and J.-W. van Holten, Phys. Rev. D 72, 024011 (2005).
[12] J. Podolsky and K. Vesely, Phys. Rev. D 58, 081501 (1998). [13] I. Sakalli and M. Halilsoy, Phys. Rev. D 74, 067501 (2006).
[14] C. W. Misner, K. S. Thorne, and J. A. Wheeler, Gravitation, (Freeman, San Fransisco, 1973). [15] E. S. Ceb-Terrab and H. P. de Oliveria, Comput. Phys. Commun. 95, 171 (1996).
Figures
FIG. 1: For k = 0, 2D (x, y) plot of the geodesics. Geodesics start from ( x0 = 0.7, y0 = −3 + 18 P j=0 j 3 )
(dashed line) and move through non-intersected trajectories. The non-intersected trajectories represent the regular motion.
FIG.2: In case k > −1 2e−
7
8, k = 1 is chosen for the 2D (x, y) plot of the geodesics.The
initial positions are ( x0 = 0.7, y0 = −3 + 18 P j=0 j 3 )
(dashed line). Intersected trajectories sig-nal the existence of chaos. Two trajectories, which have different y0 initial points initially
accelerate in +x-direction contrary to the others. FIG. 3: In case k < −1 2e− 7 8, k = −0.21 (−4 3e− 3 2 < −0.21 < −1 2e− 7 8) is chosen 2D (x, y)
plot of the geodesics.The initial positions are ( x0 = 0.7, y0 = −3 + 18 P j=0 j 3 ) (dashed line). The chaotic behavior is evident from the trajectories. The symmetry in Eq. (13) shows itself along the y-axis.
FIG. 4: Poincar´e sections of (x′, x) for k = −0.21 and H = 0.2 across y = 0 KAM surface.
Some points are distributed randomly in a finite region to form a chaotic sea, however the large island surrounded by the chaotic sea indicates the existence of quasi-periodic orbits. (Here, x → q1, and x′ → p1)
k=0 –6 –4 –2 0 2 4 6 y 0.2 0.4 0.6 0.8 1 x
k=1 –6 –4 –2 0 2 4 6 y 0.2 0.4 0.6 0.8 1 x
k=–0.21 –15 –10 –5 0 5 10 15 y 5 10 15 20 x
H = -.2000000000e–2; –80 –60 –40 –20 0 20 40 60 80 p1 0 5 10 15 20 25 30 q1