Gravitational Geons on the Brane by Daniel John Kermode B. Math, University of Waterloo, 1989 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF SCIENCE in THE COLLEGE OF GRADUATE STUDIES (Interdisciplinary Graduate Studies) THE UNIVERSITY OF BRITISH COLUMBIA (Okanagan) April 2011 c Daniel John Kermode, 2011 ⃝ Abstract Gravitational geons, if they exist, represent a candidate for cosmological dark matter observations. Gravitational Geons that are nonsingular, asymptotically flat, topologically trivial spacetime solutions are not possible in General Relativity. The purpose of this thesis is to investigate whether such geons are possible on a 3 dimensional brane embedded in a 4+1 dimensional spacetime. Our approach is to identify the characteristics of a candidate spacetime solution on the brane, choose an example spacetime that exhibits those characteristics and determine if it corresponds to a gravitational geon. ii Preface The original idea to investigate gravitational geons in brane-world cosmologies is that of my supervisor Dan Vollick. Together, we collaborated to author an article entitled Gravitational Geons on the Brane which has been submitted for publication at the journal General Relativity and Gravitation. Chapter 2 of this thesis is a version of that article. The investigative effort involved in the article required extensive use of both analytical and numerical methods. Although we both participated in each effort, Professor Vollick’s contributions were largely analytical whereas mine were largely numerical. The writing of the article including the sourcing and verification of references was done by myself with extensive advice and support by Professor Vollick. iii Table of Contents Abstract . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ii Preface . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . iii Table of Contents . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . iv List of Figures . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . vi List of Symbols . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . viii Acknowledgments . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . x 1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 1.1 Differential Geometry and General Relativity . . . . . . . . . . . 2 1.1.1 The Einstein Field Equations . . . . . . . . . . . . . . . . 2 1.1.2 The Kretschmann Scalar . . . . . . . . . . . . . . . . . . 4 1.1.3 Junction Conditions . . . . . . . . . . . . . . . . . . . . 4 1.2 The Brane-World . . . . . . . . . . . . . . . . . . . . . . . . . . 6 1.3 Gravitational Geons . . . . . . . . . . . . . . . . . . . . . . . . . 12 1.4 Static Space-Time Solutions on the Brane . . . . . . . . . . . . . 16 1.4.1 Brane-World Black Holes . . . . . . . . . . . . . . . . . 16 1.4.2 Brane-World Wormholes . . . . . . . . . . . . . . . . . . 19 Gravitational Geons in Other Theories . . . . . . . . . . . . . . . 20 Gravitational Geons in Brane-World Cosmologies . . . . . . . . . . 23 2.1 2.2 24 25 1.5 2 The Candidate Space-Times . . . . . . . . . . . . . . . . . . . . Methods Of Investigation . . . . . . . . . . . . . . . . . . . . . . iv 2.2.1 The Weak Field Approximation . . . . . . . . . . . . . . 25 2.2.2 Numerical Iteration . . . . . . . . . . . . . . . . . . . . . 26 2.2.3 Direction Field Plots . . . . . . . . . . . . . . . . . . . . 27 2.2.4 Potential Singular Points . . . . . . . . . . . . . . . . . . 28 2.3 Weak Field Behaviour . . . . . . . . . . . . . . . . . . . . . . . 29 2.4 General Behaviour . . . . . . . . . . . . . . . . . . . . . . . . . 31 2.4.1 Two Critical Points – High Range . . . . . . . . . . . . . 36 2.4.2 Two Critical Points – Medium Range . . . . . . . . . . . 37 2.4.3 Two Critical Points – Low Range . . . . . . . . . . . . . 38 2.4.4 No Critical Points . . . . . . . . . . . . . . . . . . . . . . 39 Discussion And Conclusion . . . . . . . . . . . . . . . . . . . . . . . 43 3.1 The Weak Field . . . . . . . . . . . . . . . . . . . . . . . . . . . 43 3.2 Identifying m and l . . . . . . . . . . . . . . . . . . . . . . . . . 44 3.3 Suggestions for Further Analysis . . . . . . . . . . . . . . . . . . 46 3.4 Significance of the Result . . . . . . . . . . . . . . . . . . . . . . 49 Bibliography . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51 3 v List of Figures Figure 2.1 The interface of a custom application built to investigate our candidate spacetimes using a fourth-order Runge-Kutta numerical iteration. Figure 2.2 27 An example of a direction field plot showing numerical instabilities in A(r) (here plotted for l = 1, m = 1.189 and 1.25 < r < 1.42). Figure 2.3 . . . . . . . . . . . . . . . . 28 B(r) and A(r) plotted against radius r for l = 1, m = 0.01. A fourth-order numerical iteration of (2.5) is labelled A(r) and the analytical weak field result is labelled weak A(r). These two approaches give almost identical results for m l << 1. . . . . . . . . . . . . . . . . . . . . . . . . . Figure 2.4 α1 and β1 plotted against Figure 2.5 α1 and β1 plotted against Figure 2.6 A(r) iteratively plotted against radius r for l = 1 and m = 1.25 against a m l m l for the first zero of B′ r + 4B. for the second zero of . . . . . . . . B′ r + 4B. . . . . . . . direction field plot. Here, both the convergence to A = 0 and A = −1.84 at the first zero of Figure 2.7 B′ r + 4B can be clearly seen. −α1 β1 30 33 33 ≈ . . . . . . . . . . 36 Reverse iteration of A(r) plotted against radius r for l = 1 and m = 1.192 from the second zero of B′ r + 4B to the first, against a direction field plot. The solution is incompatible with the analytical result that A′ < 0 at the first zero for λ = 0 solutions. Figure 2.8 . . . . . . . . . . . . . . . . . . . . . . 37 This direction field plot of A(r) against radius r shows that A(r) is numerically unstable at the zeros of B′ r + 4B. Here we choose l = 1 and m = 1.189 such that the zeros are found at r = 1.281, 1.392. Figure 2.9 . . . . . . . . . . . . Forward iteration of A(r) plotted from r = 0 to the first zero of l = 1 and m = 1.189. A → −α1 β1 B′ r + 4B 40 to the first . . . . . . . . . . . . . . . . . . . . . . vi 39 for ≈ 1.38 as the solution approaches the first zero. Figure 2.10 Reverse iteration of A(r) plotted from the second zero of for l = 1 and m = 1.189. B′ r + 4B 40 Figure 2.11 Forward iteration of A(r) plotted from the second zero of B′ r + 4B for l = 1 and m = 1.189. . . . . . . . . . . . . . . . . . . . . . . . . . . 41 Figure 2.12 The function A(r) plotted against radius r for l = 1 and m = 1.189. A fourthorder Runge-Kutta numerical iteration was unable to navigate the zeros of B′ r +4B (at r = 1.281, 1.392), so the result is constructed piecewise between these points. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41 Figure 2.13 Function A(r) plotted against radius r. Here we choose l = 1 and chart a variety of values for m in the specified region. A fourth-order Runge-Kutta numerical iteration is used to generate each result. Figure 3.1 . . . . . . . . . . . B(r) plotted for m = 1 and 1 ≤ l ≤ 33.3. Note that all solutions have so these solutions all correspond to gravitational geons. vii m l < 42 32 27 . . . . . . . . . 45 List of Symbols The following list contains primarily symbols whose physical and mathematical meaning remain unchanged throughout this thesis. Cited page numbers refer to the first occurrence of each symbol. Rµν , (5) RAB Ricci tensor 4D/5D components (contracted from the Riemann tensor) . . . . . . . . . 2 R, (5) R Ricci scalar 4D/5D (contracted from the Ricci tensor) . . . . . . . . . . . . . . . . . . . . . . . . 2 energy momentum tensor 4D/5D components . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 Tµν , (5) TAB gµν , (5) gAB metric tensor 4D/5D components . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 Λ, Λ5 cosmological constant 4D/5D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 κ , κ5 B(r), B A(r), A O() Sµν , (5) SAB Kµν , (5) KAB Ln Rµνσ φ (5) R ABCD ∇µ (5)C ABCD Eµν Γαµν √ −g ∂α ηµν gravitational coupling constant 4D/5D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 gtt metric tensor component of static, spherically symmetric space-time . . . . . . . . 3 grr metric tensor component of static, spherically symmetric space-time . . . . . . . .3 Order of notation - terms of specified order and higher are neglected . . . . . . . . . . 4 surface energy momentum tensor 4D/5D components at a hypersurface . . . . . . . . 5 extrinsic curvature tensor 4D/5D components . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 Lie differentiation operator . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6 The Riemann Curvature Tensor 4-dimensional components . . . . . . . . . . . . . . . . . . . 9 Riemann curvature tensor 5-dimensional components . . . . . . . . . . . . . . . . . . . . . . . . 9 covariant differentiation operator . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9 Weyl curvature tensor 5-D components . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9 electric projection of the 5-D Weyl curvature tensor . . . . . . . . . . . . . . . . . . . . . . . . . 10 Christoffel symbol for the Levi-Civita connection . . . . . . . . . . . . . . . . . . . . . . . . . . . 14 square-root of the determinant of the matrix of metric tensor components . . . . . 14 partial derivative operator . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14 Minkowski space-time metric tensor components . . . . . . . . . . . . . . . . . . . . . . . . . . . 16 viii b(r), b a(r), a Bµνσ variation in B(r), defined in the weak field by b(r) = 1 − B(r) . . . . . . . . . . . . . . . .25 variation in A(r), defined in the weak field by a(r) = 1 − A(r) . . . . . . . . . . . . . . . .25 magnetic projection of the 5-D Weyl curvature tensor . . . . . . . . . . . . . . . . . . . . . . . 47 ix Acknowledgments I would like to thank my supervisor, Professor Dan Vollick for his support and patience. It was Professor Vollick’s original idea which led to the paper Gravitational Geons on the Brane which we collaborated to write and on which this thesis is based. I would also like to thank my wife Myra as well as my family and friends who have supported and encouraged me through my degree. I particularly appreciate the editorial review provided by Stephanie Halldorson. For their excellent academic instruction and unqualified support, my gratitude goes to the professors of physics and mathematics at UBCO, in particular to Professor Erik Rosolowsky, Professor Murray Neuman and Professor Sylvie Desjardins. This research was supported by the Natural Sciences and Engineering Research Council of Canada. x Chapter 1 Introduction Einstein’s general theory of relativity was published in 1915 and remains unsurpassed as a predictive mathematical model of gravitation. That having been said, general relativity also predicts the existence of space-time singularities in gravitational collapse (i.e. in black holes) and in cosmological space-times. In particle physics, quantum field theory shares the predictive success of general relativity. In quantum field theory, one obstacle to unifying gravitation with the other forces of nature is the hierarchy problem which is related to the large discrepancy between the electro-weak scale and the Planck scale. The RandallSundrum Brane-world cosmological model [1] was proposed as a new way to solve the hierarchy problem. In this investigation, we consider Randall-Sundrum type brane-world cosmologies in which the three spatial dimensions of the observable universe are modeled by a 3 dimensional membrane or 3-brane embedded in a 4+1 dimensional bulk space-time. Long before brane-world theory was introduced, John Wheeler introduced the concept of a gravitational-electromagnetic entity or geon [2] with the hope that a configuration of gravitational and electromagnetic fields could form a persistent space-time disturbance that would ultimately provide a self-consistent picture of what a particle of matter is. A gravitational geon is a non-singular configuration of the gravitational field, without horizons that persists for a long period of time. General relativity does not allow for static, topologically trivial and spherically symmetric gravitational geon solutions of the Einstein field equations. However, 1 in the Randall-Sundrum brane-world model, the 4+1 dimensional Einstein field equations in the bulk induce field equations on the brane that differ from those of general relativity. In this thesis, we consider whether there are static, topologically trivial and spherically symmetric solutions to the Einstein field equations on the brane that correspond to gravitational geons. We construct these space-time solutions on the brane in order to prove that they could exist in Randall-Sundrum type brane-world cosmologies. Our approach does not yield a general solution for gravitational geon spacetimes on the brane. However, using approximations and numerical methods, it does show that such space-time solutions exist. Although the solutions that we find cannot correspond to conventional particles of matter, they could correspond to dark matter. 1.1 Differential Geometry and General Relativity It is assumed that the reader is familiar with the mathematics of differential geometry and Einstein’s general theory of relativity. In this thesis, we work extensively with the Einstein field equations and certain ideas specific to the study of curved manifolds. The intention of this section is to state and briefly review those results that form the foundation of this thesis. For further details, see [3–5]. We adopt the convention that the gravitational constant G and the speed of light c are dimensionally adjusted to unity. 1.1.1 The Einstein Field Equations In general relativity, the Einstein field equations are 1 Rµν − gµν R = −gµν Λ + κ Tµν , 2 (1.1) where Rµν is the Ricci tensor (a contraction of the Riemann curvature tensor), R is the Ricci scalar (a contraction of the Ricci tensor), gµν is the space-time metric, Λ is the cosmological constant, κ = 8π is the gravitational coupling constant and Tµν is the energy-momentum tensor representing the source of the gravitational field. 2 The left hand side of (1.1) describes the curvature of space whereas, ignoring the cosmological constant for the moment, the right hand side describes the distribution of matter and energy. In Misner, Thorne and Wheeler, [3], this interaction is described as ”Space tells matter how to move” and ”Matter tells space how to curve”. Einstein added the cosmological constant to this equation to account for the general belief at the time that the universe is static on a cosmological scale. The space-time metric is of the form ds2 = gµν dxµ dxν . (1.2) The metric for a static, spherically symmetric space-time solution can be written in the form (adopting the - + + + signature convention for this thesis) ds2 = −B(r)dt 2 + A(r)dr2 + r2 (d θ 2 + sin2 θ d φ 2 ), (1.3) and the Ricci scalar is given by B′ B′′ − R= AB 2AB ( A′ B′ + A B ) 2 + Ar ( B′ A′ − B A ) + 2 2 − 2. 2 Ar r (1.4) In general relativity, for space-time solutions on scales much smaller than the cosmos (i.e. stars and galaxies), it is reasonable to take Λ = 0 in 1.1. Vacuum solutions have Tµν = 0 so that the Einstein field equations of general relativity reduce to Rµν = 0 (1.5) In this thesis, we are investigating vacuum solutions on the brane. Insofar as the gravitational field equations on the brane are induced from the Einstein field equations in an embedding higher dimensional space-time, the field equations on the brane do not reduce to (1.5) in general. We will see that there is a correspondence between vacuum solutions on the brane and solutions to the Einstein field equations in 3+1 dimensional general relativity with a trace-free energy-momentum tensor (i.e. with R = 0 in the absence of a cosmological constant). 3 1.1.2 The Kretschmann Scalar In our investigations, we be test candidate gravitational geon space-times. In particular, we are interested in static, spherically symmetric space-times that are non-trivial, non-singular and without horizons. For space-time metrics of the form (1.3), the Kretschmann scalar is derived from the Riemann curvature tensor as Rµνρσ Rµνρσ = 4K12 + 8K22 + 8K32 + 4K42 (1.6) where [6] 1 K1 = A ( B′′ (B′ )2 A′ B′ − − 2B 4B2 4AB ) (1.7) K2 = B′ 2ABr (1.8) K3 = −A′ 2A2 r (1.9) K4 = A−1 Ar2 (1.10) As this is a sum of squares of all components of the Riemann curvature tensor, a divergence in the Kretschmann scalar will correspond to a singularity in the spacetime. It is clear from (1.9) and (1.10) that lim A(r) = 1 + O(rn ) with n ≥ 2 to r→0 avoid divergence of the Kretschmann scalar. Similarly, from (1.8) it follows that B′ (r) lim = O(rm ) with m ≥ 1. r→0 B(r) 1.1.3 Junction Conditions In the Randall-Sundrum model, the brane is a hypersurface embedded in a higher dimensional space-time. On the brane, we also encounter embedded hypersurfaces corresponding to potential singular points in the numerical analysis of our candidate brane-world space-times. In each situation involving an embedded hypersurface, the Israel-Darmois junction conditions must be met across the hypersurface. 4 Taking y to be a Gaussian normal coordinate orthogonal to the hypersurface (and letting y = 0 at the hypersurface without loss of generality), the first junction condition is simply that there be no discontinuity in the metric lim gµν − lim− gµν = 0 y→0+ (1.11) y→0 To get the second, we first observe that, in the embedding space-time around a hypersurface, the energy momentum tensor can be decomposed [7] as + − Tµν = θ (l)Tµν + θ (−l)Tµν + δ (l)Sµν (1.12) where θ (l) is the Heaviside theta function, δ (l) is the Dirac delta function (both ± is the energy-momentum evaluated at the proper distance l from the hypersurface), Tµν tensor evaluated on either side of the hypersurface and Sµν is the surface energymomentum tensor corresponding to the hypersurface. The second Israel-Darmois junction condition states that the surface energymomentum tensor corresponding to a hypersurface is related to the discontinuity in the extrinsic curvature Kµν across the hypersurface by [3, 8] Sµν = ) 1 ( [Kµν ] − [K]hµν 8π (1.13) µ where hµν is the induced metric on the hypersurface, K = Kµ and the square brackets, indicate the discontinuity in a physical quantity crossing the hypersurface (i.e. [K] = K + − K − ). In the Randall-Sundrum brane-world cosmology, we will look at solutions to the gravitational field equations on the brane embedded in a higher dimensional bulk space-time. The solution on the hypersurface is known, but the solution in the embedding space-time is not. The effective Einstein field equations on the brane include terms involving Sµν . (1.13) is used to show that these terms are zero for static, vacuum space-time solutions on the brane. We will also examine solutions to the 3+1 dimensional Einstein field equations that can only be constructed piece-wise around embedded hypersurfaces where we find numerical instabilities. The solution on the hypersurface is unknown but the solution in the embedding space-time is. In this situation, the extrinsic curvature 5 can be evaluated using [4] 1 Kµν = Ln gµν 2 (1.14) where Ln is the Lie Derivative with respect to the normal of the hypersurface in question. The extrinsic curvature is evaluated in the limit approaching the hypersurface from either side. Should there be a discontinuity in Kµν across the hypersufrace, Sµν ∕= 0 in general. This is generally taken to mean that there is a gravitational source (other than the gravitational field) present at the hypersurface. In this thesis, we are intested in vacuum solutions and so require that gravitational geon solutions have continuity of the extrinsic curvature across every such hypersurface. 1.2 The Brane-World Einstein’s theory of general relativity is a 4 dimensional theory, the dimensions being the 3 spatial dimensions of our universe and time. Space-time models that incorporate more than 4 dimensions date back to the theory of Kaluza who showed that 5 dimensional general relativity could incorporate both the 4 dimensions of Einstein’s gravitational theory and Maxwell’s theory of electromagnetism [9]. As we only experience a 3+1 dimensional space-time, Kaluza made the assumption that the metric would not change with respect to the extra dimension. In effect, physics would take place on a 4 dimensional hypersurface in a 5 dimensional space-time [11]. Klein showed that this condition would arise naturally if the 5th dimension had a circular topology on a small enough scale [10]. That being the case, the energies of all modes above the ground state of a Fourier expansion of the periodic dependence of any physical field on the 5th dimension, could be made so high as to be unobservable. The success of this unification of gravitation and electromagnetism led physicists to theories with higher dimensions with the hope of unifying gravitation and electromagnetism with the strong and weak forces. 10 dimensional string theories and eventually the 11 dimensional M-theory emerged from these efforts. In string theory, the 6 extra dimensions are compact and small as is the extra dimension in the Kaluza-Klien theory. M-theory extends 10-dimensional string theory with one 6 extra dimension the size of which depends on the string coupling strength. M-theory introduced the idea that there could be an extra dimension whose scale is large relative to the fundamental scale. This idea forms the basis of RandallSundrum models of 5-dimensional gravity in which the three spatial dimensions of the empirical universe are a 3 dimensional membrane or 3-brane embedded in a 4+1 dimensional space-time (4+1 dimensional referring to four spatial dimensions and time). Particles and forces of the standard model are confined to the brane whereas the gravitational field can propagate through the bulk. The Randall-Sundrum models were proposed as a new mechanism for solving the hierarchy problem. The hierarchy problem is related to the fact that the fundamental energy scale at which the quantum effects of gravity become strong (known as the Planck scale) is so much larger than the energy scale at which the quantum effects of the other forces become strong. In the Randall-Sundrum model, gravitation is free to propagate through the bulk whereas the other forces are confined to the brane. The fundamental energy scale in the bulk can therefore be much smaller than the Planck scale, which remains the effective fundamental energy scale on the brane. It is the Randall and Sundrum [1, 12] single brane model which provides the mathematical framework for our investigations. The effective 3+1 dimensional Einstein field equations on the brane were derived by Shiromizu, Maeda, and Sasaki [13]. An excellent review of this work is provided by [14] (see also [15– 20]). The Einstein field equations in the 4+1 dimensional bulk space-time are (5) RAB − 1 (5) (5) gAB R = −(5) gAB Λ5 + κ52 (5) TAB 2 (1.15) where the terms are the same as for the 3+1 dimensional case (1.1) except (5) preceding a tensor or 5 following a scalar indicates the appropriate dimensionality (the convention that upper case latin indices indicate bulk tensors with lower case greek indices indicating 3+1 dimensional tensors on the brane is also used). The metric in the bulk will be (5) ds2 = (5) gAB dxA dxB 7 (1.16) If we let y be a Gaussian normal coordinate orthogonal to the brane (letting y = 0 at the brane without loss of generality) then (5) where (4) g α µν (x , y) ds2 = (4) gµν (xα , y)dxµ dxν + dy2 (1.17) corresponds to induced metric on the hypersurfaces at y = constant. The 3+1 dimensional metric on the brane is therefore (4) g α µν (x , 0) (which we will generally write as gµν ). Our objective in this section is to determine the induced Einstein field equations on the brane. We begin by taking the string theory idea that the fields of the standard model (i.e. electromagnetic, strong, weak and all matter) are confined to the brane. Gravitation is the only field that is permitted to extend into the bulk. Thus if Tµν is the 4 dimensional energy-momentum tensor of the particles and fields confined to the brane then the total energy momentum tensor on the brane is brane = Tµν − λ gµν Tµν (1.18) where λ is the brane tension which tends to prevent the gravitational field from leaking into the bulk. Insofar as this model represents a departure from standard general relativity, the accuracy to which general relativity predicts gravitation in experiment puts a lower limit on the brane tension of λ > (1 TeV)4 . In (1.17), the use of a 4-dimensional tensor in a 5-dimensional equation is fairly explicit by nature of the coordinate system used. However, we wish to decompose the the 5-dimensional energy-momentum tensor into a part that resides on the brane and those parts that do not. This requires that we include the tensor brane , which is 4-dimensional in a 5-dimensional equation. We do this by defining Tµν the 4-dimensional tensor in terms of a 5-dimensional tensor as brane Tµν = (5) brane A TAB n = (5) brane A B TAB gµ gν , (5) brane B TAB n = 0 (1.19) brane off the brane is not where nB is normal to the brane. The behaviour of (5) TAB important insofar as we decompose the bulk energy-momentum as 8 (5) brane TAB = (5) T˜AB + TAB δ (y) (1.20) where the delta function expresses the confinement of standard model particles and fields to the brane at y = 0. We can now rewrite (1.15) incorporating this to give (5) RAB − ( ) 1 (5) (5) brane gAB R = −(5) gAB Λ5 + κ52 (5) T˜AB + TAB δ (y) 2 (1.21) In order to determine the induced Einstein field equations on the brane, in addition to the induced energy-momentum tensor, we need to know the induced Riemann curvature tensor on the brane. The bulk Riemann curvature tensor is projected onto 4 dimensional hypersurfaces y = const. (with extrinsic curvature corrections) as Rµνσ ρ = (5) RABCD gµ A gν B gσ C gρ D + 2Kµ [σ Kρ ]ν (1.22) in which the square brackets denote anti-symmetrization. The change in the extrinsic curvature Kµν along the y = const. surfaces is projected from the bulk Ricci tensor via the Codazzi equation ∇ν K ν µ − ∇µ K = (5) RAB gµ A nB , (1.23) where K = K ν ν . Observing that the 5-dimensional Riemann curvature tensor can be written as (5) RABCD = (5)CACBD + ) 1 2 ((5) gA[C (5) RD]B − (5) gB[C (5) RD]A − (5) gA[C (5) gD]B (5) R, 3 6 (1.24) where (5)CACBD is the Weyl tensor (which is trace free). Further observing that (5) T˜ = 0 (i.e. there are no stresses other than the cosmological constant and those AB on the brane) then combining (1.21) and (1.22) gives an expression for the induced Einstein field equations on the brane 9 ( ) 1 1 2 1 brane brane Rµν − gµν R = − Λ5 gµν + − T gµν + KKµν T 2 2 3 µν 4 1 − Kµ α Kαν + (K αβ Kαβ − K 2 )gµν − Eµν , (1.25) 2 where Eµν = (5)CACBD nC nD gµ A gν B (1.26) is the electric projection of the 5 dimensional Weyl tensor and, by virtue of Weyl tensor symmetries, is trace-free. Now, Randall-Sundrum brane-world cosmologies are Z2 -symmetric (have mirror symmetry) so that the bulk space-time looks the same approaching the brane from one side as it does emerging from it on the other, but with the normal reversed. Consequently, we can say that the extrinsic curvature in the limit as y → 0 approaching the brane from either side satisfies lim KAB = − lim− KAB y→0+ (1.27) y→0 The Israel Darmois junction conditions for a 4+1 dimensional bulk space-time are derived by integrating (1.21) across the brane taking the limit from either side giving lim gAB − lim− gAB = 0, y→0+ lim KAB − lim− KAB = y→0+ y→0 (1.28) y→0 −κ52 ( (5) brane TAB + ) 1 brane T gAB , 3 (1.29) which, applying (1.27) and (1.18) gives the extrinsic curvature on the brane Kµν ) ( 1 2 1 = − κ5 Tµν + (λ − T )gµν , 2 3 (1.30) So we can write the contributions from the extrinsic curvature in (1.25) in terms of the energy momentum tensor. This contribution is written 10 Sµν = ( ) 1 1 1 T Tµν − Tµα T α ν + gµν 3Tαβ T αβ − T 2 . 12 4 24 (1.31) where, recalling 1.18, Tµν is the 4-dimensional energy-momentum tensor of the particles and fields confined to the brane. If there are any stresses in the bulk other than the bulk cosmological constant (i.e. if (5) T˜AB ∕= 0) , they will also affect the induced Einstein field equations on the brane. In our case we assume there are none so the induced Einstein field equations on the brane become 1 κ2 Rµν − gµν R = −gµν Λ + κ 2 Tµν + 6 Sµν − Eµν , 2 λ (1.32) where κ 2 is an effective coupling constant inherited from the fundamental coupling constant and the cosmological constant on the brane Λ= ) 1( Λ5 + κ 2 λ , 2 (1.33) will be zero as we assume that the bulk cosmological constant will be balanced by the brane tension λ . We will be working only with vacuum solutions on the brane (i.e. Tµν = 0). As a consequence, equation (1.32) simplifies to 1 Rµν − gµν R = −Eµν . 2 (1.34) We know that Eµν is trace free (i.e. Eµ µ = 0) so contracting both sides of (1.34) yields Rµ µ = Eµ µ = 0, Rµν = −Eµν , (1.35) and that, for vacuum solutions, ∇µ Eµν = 0, 11 (1.36) where ∇µ is the covariant derivative on the brane. In standard 3+1 dimensional general relativity (with no cosmological constant), the Einstein field equations for a solution to R = 0 (so not a vacuum solution in general) will reduce to Rµν = 8π Tµν , (1.37) with ∇µ Tµν = 0. We therefore make the correspondence 8π Tµν = −Eµν . (1.38) The result is that, as stated by Dadhich, Maartens, Papadopoulos and Rezania [21] A stationary general relativity solution with trace-free energy-momentum tensor gives rise to a vacuum brane-world solution in 5-dimensional gravity. The candidate gravitational geon space-time solutions that we will be investigating are stationary vacuum solutions on the brane. In standard general relativity, such solutions would have to satisfy Rµν = 0. However, in the brane-world cosmologies considered herein, we solve R = 0 to obtain non-vacuum (in general), 3+1 dimensional space-time solutions to general relativity and these will correspond to 4+1 dimensional, vacuum solutions on the brane. 1.3 Gravitational Geons In 1955, John Archibald Wheeler coined the term geon as an abbreviation for gravitational-electromagnetic entity [2]. He wrote the following description of a geon: Associated with an electromagnetic disturbance is a mass, the gravitational attraction of which under appropriate circumstances is capable of holding the disturbance together for a long time in comparison with the characteristic periods of the system. Wheeler’s intention was to provide a self-consistent picture of what a body of matter is. More than half a century later, we still do not have a consensus on what 12 that picture should be. In quantum field theory, particles of matter are considered to be field quanta but there is no accepted understanding of how the gravitational field should be quantized. General relativity provides an effective model of gravitation in terms of the geometry of space-time, but in its original form, general relativity does not explain what matter is or how the other forces of nature come about. In the years immediately following the introduction of the geon, work focused on geons that included electromagnetic waves. Regge and Wheeler introduced the idea of gravitational geons in 1957 [22]. These were defined as geons that consisted only of gravitational waves. In general relativity, the gravitational field corresponds to curvature in a 3+1 dimensional space-time. Solutions to the field equations of general relativity are considered topologically trivial if they are topologically equivalent to a Minkowski space-time (i.e. the solution can be obtained by continuously deforming a flat 3+1 dimensional space-time). In general relativity (and Kaluza’s 5 dimensional generalization [9]), there are no non-singular, topologically trivial solutions representing a gravitational field of non-vanishing mass [23, 24] (see also [25, 26] and [27] in particular for historical context). Our investigations are concerned with static, topologically trivial, non-singular gravitational geon solutions for brane-world cosmologies. We know that such solutions are not possible in 3+1 dimensional general relativity, but in brane-world cosmologies, the induced Einstein field equations on the brane differ from those of general relativity. In particular, in brane-world cosmologies, vacuum solutions on the brane correspond to solutions to R = 0 in general relativity whereas in vacuum solutions of general relativity require Rµν = 0. It is useful to briefly consider why gravitational geon solutions cannot satisfy the condition that Rµν = 0. In his 1923 theorem, Birkhoff [28] showed that the Schwarzschild metric is the unique spherically symmetric vacuum solution in general relativity. As we are interested in non-singular, spherically symmetric vacuum solutions without horizons of the form (1.3), this result is sufficient for our purposes as the Schwarzschild metric has both a horizon and is singular at its radial origin. However, Einstein and Pauli [23, 24] did not impose spherical symmetry 13 and worked with Kaluza’s 5 dimensional generalization of general relativity (which is a foundational concept for string theory and, consequently brane-world cosmologies). Their work is therefore more relevant to the study of gravitational geons in brane-world cosmologies than Birkhoff’s theorem. For that reason, I present here a simplification to the spherically symmetric case, of Einstein and Pauli’s work . In 3+1 dimensional general relativity, Einstein and Pauli considered the variation of the Ricci tensor δ Rµν = −∇α (δ Γαµν ) + ∇ν (δ Γαµα ). (1.39) We can then define the vector density ℵα by ∂α ℵ α = √ −ggµν δ Rµν , (1.40) so that ℵα = √ −g(gµα δ Γνµν − gµν δ Γαµν ). (1.41) We are interested in vacuum solutions, so we only consider variations that preserve Rµν = 0. For any such variation, we have that δ Rµν = 0 and hence ∂α ℵα = 0. (1.42) In the following treatment we consider only space-time solutions of the form (1.3) thus, gtt and grr are both functions only of r. In our notation, we continue to use lower-case Greek indices (α , β , etc.) to indicate tensors on a 3+1 dimensional space-time; we use specific lower-case latin letters (i, j, k) to indicate tensors on 3 dimensional space and we use [t, r, θ , φ ] as the indices of 3+1 dimensional space-time in spherical polar coordinates. To avoid singularities (and for consistency with the candidate space-times that we examine in this thesis), we consider only solutions for which gtt and grr are both non-zero and finite for all r ≥ 0. The variation that we are considering can be written as a coordinate transformation x˙α = xα + ξ α . 14 (1.43) Retaining only terms of the first order in the infinitesimal variation ξ α , we have δ Γαµν = −∂µ ∂ν ξ α + ∂σ (ξ α )Γσµν − ∂µ (ξ σ )Γσα ν − ∂ν (ξ σ )Γσα µ − ξ σ ∂σ Γαµν . (1.44) We are interested in coordinate transformations that preserve the static nature of the space-time solutions. Einstein and Pauli [24] considered specific transformations that lead to an integral theorem that singles out the regular solutions of Rµν = 0 and are defined by ξ t = c1t, ξ r = 0, ξ θ = 0, ξ φ = 0. (1.45) Substituting (1.44) into (1.41), (1.42) therefore becomes √ ∂i ( −ggtt Γtti ) = 0, (1.46) where i only sums over the 3 spatial dimensions since we are only considering static space-time solutions (i.e. ∂t gµν = 0). This simplifies to ∂r (√ ) −ggtt grr ∂r (gtt ) = 0. (1.47) Multiplying (1.47) through by gtt gives ) ( √ √ ∂r (gtt ) −ggtt grr ∂r (gtt ) = ∂r gtt −ggtt grr ∂r (gtt ) , (1.48) which integrating over a 3 dimensional, singularity-free region and applying Gauss’ theorem gives ∫ V √ −ggtt grr (∂r (gtt ))2 dV + ∮ S √ −ggrr ∂r (gtt )nr dS = 0. (1.49) As the surface S approaches infinity, gtt → −1 and grr → 1 so that ∮ S √ −ggrr ∂r (gtt )nr dS → 0. Therefore, 15 (1.50) ∫ V √ −ggtt grr (∂r (gtt ))2 dV → 0, as the boundary of the volume V approaches infinity. As gtt < 0, grr (1.51) √ −g > 0 and > 0 for all r ≥ 0, it must be the case that gtt = constant. For vacuum solutions we therefore have Ri j = 0, but in three dimensions this implies that the space is Euclidean. The requirement that Rµν = 0 therefore implies that there can be no static, spherically symmetric, topologically trivial and non-singular space-time solutions (i.e. gravitational geons of the type we are investigating) in 3+1 dimensional general relativity. However, in Randall-Sundrum brane-world cosmologies, vacuum space-time solutions must satisfy R = Rα α = 0 but not Rµν = 0. Thus, gravitational geon solutions of brane-world space-times may therefore be possible. 1.4 Static Space-Time Solutions on the Brane This thesis explores the possibility that space-time solutions that correspond to gravitational geons may exist in brane-world cosmologies. Insofar as there has been considerable interest in brane-world cosmological models, there are recently published articles regarding other static space-time solutions on the brane that have some aspects in common with the investigations of this thesis. The following sections briefly summarize articles that investigate solutions corresponding to black holes and wormholes on the brane. 1.4.1 Brane-World Black Holes There are a number of different ways that brane-world black holes have been considered. In 2000, Chamblin, Hawking and Reall [29] looked at the metric for a black string in a 4+1 dimensional anti de Sitter (AdS) space-time. They wrote the AdS space-time metric ds2 = e −2y l ηµν dxµ dxν + dy2 , (1.52) where ηµν is the 4-dimensional Minkowski metric and l is the AdS radius. They 16 y introduced the coordinate z = le l so that the black string metric could be written as ds2 = ) l2 ( −U (r)dt 2 +U (r)−1 dr2 + r2 (d θ 2 + sin2 θ d φ 2 ) + dz2 , 2 z (1.53) where U (r) = 1 − 2M r . They then pointed out that this metric, with an appropriate rescaling of the coordinates r and t, is in standard Schwarzschild form on the brane (at y = 0 hence z = l). A black string in the bulk thus defined gives rise to a Schwarzschild black hole on the brane. Later in the same year, Dadhich, Maartens, Papadopoulos and Rezania [21] considered the problem a different way. Instead of starting with a solution on the bulk space-time, they started (as do we) by insisting that their solution solve the induced Einstein field equations on the brane. They showed that a static, spherically symmetric black hole solution of the form (1.3) is given by −1 B(r) = A(r) = 1− ( 2M M p2 ) 1 + r ( q M˜ p2 ) 1 , r2 (1.54) where M p is the Planck mass on the brane, M˜ p is the Planck mass in the bulk, and q is a dimensionless tidal charge parameter (so named because of its correspondence to the charge parameter in the Reissner-N¨ordstrom black hole in general relativity). In general relativity, the charge parameter Q in the Reissner-N¨ordstrom black hole solution is squared and so always acts to weaken the gravitational field of the black hole (as compared to a Schwarzschild solution which is achieved in the limit as Q → 0). In their discussion of this brane-world black hole solution, the authors pointed out that the tidal charge parameter q is negative corresponding to a negative effective energy density on the brane contributed by the free gravitational field in the bulk. Consequently, the bulk effects tend to strengthen the gravitational field of the black hole as compared to a Schwarzschild black hole (achieved in the limit as q → 0). In 2002, Vollick [19] showed that solving R = 0 for space-time metrics of the form (1.3) with A(r) = B(r)−1 has the general solution 17 B(r) = A(r)−1 = 1 + α β + , r r2 (1.55) where α and β are constants. Thus the Reissner-N¨ordstrom type solution of Dadhich, Maartens, Papadopoulos and Rezania is the most general of that form. He then considered the result if the condition that A(r) = B(r)−1 is removed and B(r) = 1 − 2m r (i.e. gtt is Schwarzschild). He obtained the solution ( 2m A(r) = 1 − r )−1 [ ] 3m − 2r . λ − 2r (1.56) It was suggested that deviations from Schwarzschild geometry (which must result from the Weyl term in the induced Einstein field equations on the brane as in (1.26)) could be responsible for observations attributed to dark matter. A third distinct approach to brane-world black holes is that of Shankaranarayanan and Dadhich in 2004 [30]. They considered the possibility of non-singular black holes on the brane. They obtained a non-singular Reissner-N¨ordstrom type black hole solution in a D-dimensional, spherically symmetric space-time (where D ≥ 5) with a metric of the form ds2 = −B(r)dt 2 + A(r)dr2 + r2 dΩ2 , (1.57) where dΩ2 is the (D − 2)-dimensional angular line-element. The solution that they obtained had B(r) = A(r)−1 = 1 − 4Rg (r) , (D − 2)rD−3 (1.58) where Rg (r) = 4π ∫ r 0 ρ (x)xD−2 dx . (1.59) Having found the form of their solution, the authors proceeded to choose mass and charge distributions in order to set ρ (x) in (1.59). This led to them arriving (with some simplifying assumptions) at the expression 18 ) 16π 1 c1 ( D−1 1 − exp(−c r ) 2 (D − 2)(D − 1) rD−3 c2 ) 16π c3 1 ( 2(D−2) + 1 − exp(−c r ) 4 (D − 2)(D − 3) r2(D−3) B(r) = A(r)−1 = 1 − (1.60) where the ci ’s are constants. The induced black hole solution on the brane is non-singular and admits one horizon. In order to achieve this result, they had to relax the condition that R = 0. In our investigation, we insist on a vanishing Ricci scalar in the induced Einstein field equations on the brane (1.34). In the brane-world black hole of Shankaranarayanan and Dadhich contributions to R cannot come from the bulk Weyl tensor via Eµν as this tensor is trace-free. The non-vanishing Ricci scalar indicates a varying cosmological constant-like term, thus the gravitational geons that we investigate in this thesis differ fundamentally from their result. 1.4.2 Brane-World Wormholes As with black hole solutions on the brane, when investigating wormhole solutions on the brane there are two fundamental approaches. One can either create a wormhole solution in the bulk space-time and examine the induced solution on the brane (see for example [31]), or one can create a wormhole solution that satisfies the induced Einstein field equations on the brane. In their 2002 paper, Bronnikov and Kim [6] take the latter approach. They begin with a generic, static, spherically symmetric space-time metric on the brane, as do we (1.3). Our method is comparable to theirs insofar as they also solve for R = 0 and test for finiteness of the Kretschmann scalar (1.1.2) as a regularity criterion for the resulting space-time solutions. The brane-world wormhole solutions considered were, like our gravitational geon solutions, static, spherically symmetric, non-singular and without horizons. However, for a wormhole, the solution need only be non-singular for r > r0 for some r0 > 0. Two space-time solutions meeting this condition can be sewn together at the spherical hypersurface at r0 (which is identified as the throat). The resulting piece-wise constructed space-time forms the wormhole provided the Israel19 Darmois junction conditions (1.1.3) are met at the throat. The authors considered both symmetric and asymmetric brane-world wormhole metrics. One example given of a symmetric wormhole solution used the metric [32] ( r0 )−1 2 ds2 = −dt 2 + 1 − dr + r2 (d θ 2 + sin2 θ d φ 2 ). r (1.61) It is clear that grr → −∞ as r+ → r0 . Bronnikov and Kim assert that the cre- ation of wormholes involves the violation of the null energy condition at least in a neighbourhood of the throat. In brane-world cosmologies, a negative energy density that would provide for this violation is possible (see [33]). The surface energy-momentum tensor (1.12) associated with the throat of the wormhole may have contributions from the Weyl term in the induced field equations on the brane (1.26). The authors do not attempt to evolve the brane-world metric into the bulk in order to discover what those contributions are for any particular solution. 1.5 Gravitational Geons in Other Theories As we have seen, general relativity does not admit static, non-singular and topologically trivial gravitational geon solutions. As a consequence, most work on gravitational geons has considered concentrations of gravitational waves (see for example [34] and, for an argument that this type of solution is inadmissible as a gravitational geon in general relativity, [35, 36]) or topological geons (particles made from non-trivial spatial topology [37]). This thesis investigates whether topologically trivial gravitational geons are possible in brane-world cosmologies based on the observation that the effective Einstein field equations on the brane are different from those of 3+1 dimensional general relativity. However, brane-world space-times are not the only theories that may potentially allow the existence of topologically trivial gravitational geons. One possibility, recently investigated by Dymnikova and Galaktionov, is to assume that vacuum solutions allow a varying Ricci scalar in 3+1 dimensional general relativity. They investigated geon-like particles that they described as vacuum non-singular 20 black holes [38] in a spherically symmetric space-time with a de-Sitter centre. They examined a metric of the form (1.3) defined by −1 B(r) = A(r) ( ) 2m −r3 = 1− 1 − exp( ) r 2mr02 (1.62) where r02 = Λ3 , Λ being a cosmological constant that appears at the origin (there being no cosmological constant in the solution asymptotically). The resulting spacetime is de-Sitter at its origin (i.e. having a positive cosmological constant) but Schwarzschild asymptotically (i.e. having no cosmological constant as r → ∞). The authors showed that the solution has horizons for m greater than a critical value mcr (mcr being dependent on Λ). Solutions for m < mcr are static, spherically symmetric and topologically trivial gravitational geons. However, owing to the radially dependent cosmological constant, R ∕= 0 in general for this solution. Consequently, these solutions are not vacuum brane-world solutions of the type we are concerned with in this thesis. Another interesting example is Vollick’s 2008 paper [39] in which he looks at the possibility of stationary, asymptotically flat, non-singular and topologically trivial gravitational geons in 1+1 dimensions. He looks at two different 1+1 dimensional theories of gravity, the first being a modified Jackiw-Teitelboim theory [40, 41] of the form R + α R2 + β ∇µ ∇µ R − Λ = 8π T (1.63) where Λ is a cosmological constant (in this case Λ = 0). The second theory is based on the action L= √ −g ( ) 1 R +V (φ ) . φ (1.64) In each case, working in 1+1 dimensions and treating r as a radial-like coordinate with r ≥ 0, exact solutions that correspond to gravitational geons are found. In a published addendum [42], Vollick pointed out that each of the solutions has a jump discontinuity in d2t dτ 2 if one imposes a reflecting boundary condition at r = 0 (which is necessary if there is to be a correspondence between the derived solutions and spherically symmetrical space-time solutions in 3+1 dimensions). 21 The author removes these discontinuities for (1.63) by allowing a non-zero cosmological constant. In the original paper, it was shown that the field equations resulting from (1.64) would be solved by potentials of the form V (φ ) = A f ′ dV (φ ) = −A2 r2 f ′′ . dφ and (1.65) The choice of f (r) allows one to solve for the potential. In the addendum, the function f (r) = 1 − 2mr2 r3 + 2ml 2 (1.66) is shown to remove the jump discontinuity and has Schwarzschild behaviour at large r. The existence of gravitational geons similar to what we investigate herein have therefore been explored, but not previously in brane-world cosmologies. 22 Chapter 2 Gravitational Geons in Brane-World Cosmologies In this thesis, we investigate whether the free gravitational field in the bulk modifies the gravitational field on the brane so that gravitational geons are possible on the brane. Our approach is to look at static, spherically symmetric solutions of the form ds2 = −B(r)dt 2 + A(r)dr2 + r2 (d θ 2 + sin2 θ d φ 2 ). (2.1) The Ricci scalar for metrics of this form is B′′ B′ R= − AB 2AB ( A′ B′ + A B ) 2 + Ar ( B′ A′ − B A ) + 2 2 − 2. 2 Ar r (2.2) As explained in the introduction, a solution to R = 0 substituted into the Einstein field equations of 3+1 dimensional general relativity, will give rise to a vacuum brane-world solution in 4+1 dimensions via a correspondence between the energymomentum tensor in the 3+1 dimensional solution and the projection of the bulk Weyl tensor in the 4+1 dimensional solution. Our approach is to solve R = 0 in (2.2) for A(r). We do this generally in the weak field, but for strong field results, we choose a specific function B(r). 23 2.1 The Candidate Space-Times We were unable to find a general solution to R = 0 for (2.2). For our strong field analysis, each potential choice of B(r) represents a candidate space-time solution that may or may not correspond to gravitational geons on the brane. We are not interested in space-time solutions with singularities so we examined the Kretschmann scalar for divergent behaviour in order to narrow down our available choices. From (1.9) and (1.10) it is clear that we must have lim A(r) = r→0 1 + O(rn ) with n ≥ 2 to avoid divergence of the Kretschmann scalar. Similarly, B′ (r) from (1.8) it follows that lim = O(rm ) with m ≥ 1. r→0 B(r) In this thesis, we take B(r) = 1 − 2mr2 , r3 + 2ml 2 (2.3) for which B′ (0) = 0, B(0) = 1, both B(r) and B′ (r) are continuous for all r ≥ 0 and 2m lim B(r) = 1 − is Schwarzschild. r→∞ r In the limit as r → ∞ the gravitational field is weak and static. If we make the further assumption that any test particles are travelling at sub-relativistic speeds, we find that gtt = −(1 + 2Φ) where Φ corresponds to the Newtonian gravitational potential. Thus knowing B(r) is sufficient for us to identify the parameter m with the gravitational mass of the space-time solution. ( )2 ( )1 For (2.3) we have that B ≈ 1 − rl for r << ml 2 3 . Thus l has the effect of controlling the strength of a positive cosmological constant like term (with Λ = l32 ) that appears in the limit as r → 0 (see 3.2). To be certain that the resulting space-times are non-singular and without hori- zons, we only consider solutions A(r) and B(r) greater than zero and finite for all r ≥ 0. Considering that our choice of B(r) is not unique, having chosen (2.3) and been unsuccessful at finding a general solution, we considered other candidate solutions such as those for which B(r) = 1 − 24 2mr2 . (r + l)3 (2.4) None of the alternate expressions for B(r) that we looked at yielded a general solution to R = 0 in (1.4). We decided to focus on (2.3) and use a variety of approximate and numerical methods to determine the nature of the resulting solution. 2.2 Methods Of Investigation Our investigation concerns solutions to 3+1 dimensional general relativity for which the Ricci Scalar vanishes everywhere. As a consequence, (2.2) can be written as A′ = (2BB′′ r2 − (B′ )2 r2 + 4BB′ r + 4B2 )A − 4B2 A2 . Br(B′ r + 4B) (2.5) We did not find a general solution to (2.5) for our choice of B(r) (2.3). In the absence of a general solution, we used a variety of analytical and numerical methods to investigate the behaviour of the candidate space-times. Using a combination of the methods described in the following sections, we were able to determine for which choices of the parameters m and l in (2.3) the solutions correspond to gravitational geons. 2.2.1 The Weak Field Approximation In the absence of a cosmological constant on the brane, for B(r) = A(r) = 1, our general, spherically symmetric metric (2.1) corresponds to flat space. We therefore take the weak field limit to be B(r) = 1 + b(r) and A(r) = 1 + a(r) with ∣b(r)∣ << 1 and ∣a(r)∣ << 1. In the weak field, the differential equation that we must solve (2.5) therefore simplifies (by removing terms of second order in a, a′ , b or b′ ) to a′ = b′′ r2 + 2b′ r − 2a . 2r (2.6) With an analytical solution to this, we can determine whether there is a class of functions B(r) that correspond to gravitational geons in the weak field. For the specific choice of B(r) that we have made (2.3), we will use our analytical solution to the weak field to verify the results that we obtain using numerical methods. 25 2.2.2 Numerical Iteration In the absence of a general solution to (2.5), numerical methods proved very helpful in the investigation of the behaviour of the candidate space-times. We chose to use a fourth-order Runge-Kutta method to explore (2.5) iteratively. This is expressed given an initial value problem, which in our case is A′ = f (r, A), A(r0 ) = A0 . (2.7) The Runge-Kutta method is then given by K1 K2 K3 K4 + + + 6 3 3 6 = rn + ∆r An+1 = An + rn+1 (2.8) where ∆r is the iteration length, An+1 ≈ A(rn+1 ) and K1 = ∆r f (rn , An ) ) ( K1 ∆r K2 = ∆r f rn + , An + 2 2 ( ) ∆r K2 K3 = ∆r f rn + , An + 2 2 K4 = ∆r f (rn + ∆r, An + K3 ) . (2.9) We found that working with 3rd party software (Maple) was limiting insofar as it was slow when choosing a very small iteration length, and it proved an inflexible environment for investigating the behaviour of several variables at once. We built a custom application to address these concerns. It allowed us to choose a starting point for iteration, an iteration length, the parameters m and l as well as the range of r and the granularity for the displayed values. This produced results for certain variables of interest (see fig. 2.1). These results were then plotted using Microsoft Excel. Knowing that there are potential singular points at B = 0 and B′ r +4B = 0 (2.5), using this custom application allowed us to iteratively examine the behaviour of A 26 Figure 2.1: The interface of a custom application built to investigate our candidate spacetimes using a fourth-order Runge-Kutta numerical iteration. alongside various expressions including B, B′ r + 4B and the Kretschmann Scalar (Rµνρσ Rµνρσ ) for a variety of choices of the parameters m and l. As we continued to analyse our results, we added additional expressions to the software to provide further insight into the behaviour of the solution. It should be noted that Maple provided additional choices for methods of numerical iteration. We compared our results to those provided Maple such as the Runge-Kutta-Fehlberg adaptive numeric procedure. We found that the Maple based methods were no more successful at navigating the potential singular points in our solution. However, this comparison was invaluable for finding coding errors in our custom developed software solution as well as validating the results that we obtained. 2.2.3 Direction Field Plots Numerical iteration alone was not sufficient to obtain a complete and reliable picture of the behaviour of A. In particular, we found regions of instability where the numerical solution diverged. We used the Maple mathematics and modelling software package to produce direction field plots (for example see fig. 2.2) Whereas the iterative approach provides results for our specific initial value 27 Figure 2.2: An example of a direction field plot showing numerical instabilities in A(r) (here plotted for l = 1, m = 1.189 and 1.25 < r < 1.42). problem, the direction field plots indicate where there are critical points and how the behaviour of the solution depends on changes in initial conditions. 2.2.4 Potential Singular Points There are potential singular points in (2.5) at r = 0, B = 0 and B′ r + 4B = 0. We found solutions near these points using series approximation. Specifically, we write (2.5) in the form A′ (r) = f (r)A(r) + g(r)A2 (r) . (2.10) Letting F ′ (r) = f (r) then (2.10) has the solution [6] A(r) = e−F(r) ∫ −1 . g(r)eF (r) dr We obtain series approximations of the form 28 (2.11) ∞ ∑ αnxn+s f (x) = n=0 ∞ ∑ βm xm+t g(x) = (2.12) m=0 where x = r − r0 for some potential singular point at r = r0 , the αn and βm are real constants and s and t are integer constants. We can then solve (2.11) for lim A(x). x→0 To ensure that the resulting space-times be non-singular and without horizons, we rule out solutions for which we find points where A = 0 or A → ∞. From the Kretschmann Scalar, (1.9) in particular, we know that solutions with A′ → ∞ must also be ruled out. To that end, we differentiate our approximate solution A(x) around x = 0 to obtain lim A′ (x). x→0 With an understanding of the behaviour of A and A′ in the vicinity of a potential singular point, we can characterize the behaviour of the space-time on the spherical hypersurface corresponding to that point. 2.3 Weak Field Behaviour In the weak field limit, (2.5) simplifies to (2.6) which has the solution A = 1+ b′ r c + , 2 r (2.13) where c is a constant of integration. Noting that c = 0 to ensure that A be finite at r = 0, any choice of B(r) = 1 + b(r) with ∣b(r)∣ << 1 for which b′ r is small for all r ≥ 0 would give suitable geon solutions in the weak field. In order to analyse candidate geons in the strong field, we have chosen B(r) = 2 1 − r32mr . It is useful to solve this for A(r) in the weak field for which we obtain +2ml 2 mr5 − 4m2 l 2 r2 . (2.14) (r3 + 2ml 2 )2 ( )2 ( )2 ( )1 For r << ml 2 3 we have B ≈ 1 − rl and A ≈ 1 − rl . A(r) therefore satisfies A = 1+ the requirements that A(0) = 1 and that lim A(r) = 1 + O(rn ) with n ≥ 2. Furtherr→0 more, A(r) is non-zero and continuous for all r ≥ 0. 29 It is worth noting that for large r, A differs from Schwarzschild with A ≈ 1 + mr . Also, the conditions for the weak field from section (2.2.1) are both met for r l ∣b(r)∣ = 2mr2 << 1 and r3 + 2ml 2 (2.15) ∣a(r)∣ = mr5 − 4m2 l 2 r2 << 1 (r3 + 2ml 2 )2 (2.16) << 1 and for m l << 1. For our specific choice of B(r) (see fig. 2.3), the solutions in the weak field limit correspond to gravitational geons. We also examined the behaviour of our candidate space-times for weak fields using numerical iteration (shown in fig. 2.3 for comparison). We found that the numerical results were almost identical to the analytical results for the weak field limit. Figure 2.3: B(r) and A(r) plotted against radius r for l = 1, m = 0.01. A fourth-order numerical iteration of (2.5) is labelled A(r) and the analytical weak field result is labelled weak A(r). These two approaches give almost identical results for ml << 1. 30 2.4 General Behaviour In the absence of a general solution to (2.5) for strong gravitational fields, we examined the behaviour of A(r) using a numerical iteration. Using this method, ( )2 we were unable to begin at exactly A(0) = 1. However, A ≈ 1 − rl for small r. Thus, if we choose an iteration length ∆r so that ∆rl = 10−5 then A(∆r) ≈ 1 to within 10−10 (independent of m). Also, having made the observation from direction field plots that the solution is insensitive to small changes in the initial value of A, we chose ∆r appropriately and used A(∆r) = 1 as the starting point of the iteration. We are interested in solutions for which both A(r) and B(r) are greater than zero and finite for all r ≥ 0. The zeros of B can be determined by noting that 1 B(0) = 1, lim B(r) = 1 and that B′ = 0 at r = (4ml 2 ) 2 where B has its minimum value of r→∞ 2 1 − 31 (4 ml ) 3 . Thus ⎧ ⎨ no zeros B has one zero ⎩ two zeros if if if m l m l m l √ 27 √4 27 √4 27 4 < = > (2.17) . We proceeded to examine the numerical iteration of A(r). We found that it works well provided we choose m and l for which B′ r + 4B ∕= 0 for all r > 0 but that it is unable to navigate instabilities in the solution otherwise. The zeros of B′ r + 4B are found by observing that B′ r + 4B d ′ dr (B r + 4B) = 2 3 ( 27m 4l ) . Thus r=0 = 4, lim B′ r + 4B = 4 and that r→∞ 1 0 at r = (2ml 2 ) 3 where B′ r + 4B has its minimum value of 4 − ⎧ ⎨ no zeros ′ B r + 4B has one zero ⎩ two zeros if if if m l m l m l < = > 32 27 32 27 32 27 (2.18) . Noting that there will be choices of m and l for which there will be zeros of B′ r + 4B but no zeros of B, we investigated the behaviour of A(r) around the zeros of B′ r + 4B. In general, from (2.10), we have 31 f (r) = 2BB′′ r2 − (B′ )2 r2 + 4BB′ r + 4B2 and Br(B′ r + 4B) (2.19) −4B . r(B′ r + 4B) (2.20) g(r) = Taking x = r − r0 where B′ r + 4B r=r0 = 0 we get series approximations for (2.19) and (2.20) f (x) = α1 + α2 + α3 x + O(x2 ) and x (2.21) β1 + β2 + β3 x + O(x2 ) . x (2.22) g(x) = / {0, −1, −2, ...}, solutions to (2.11) are of the form In this case, for α1 ∈ A(x) = β1 α1 − −1 β1 α2 −β2 α1 − α1 α1 (α1 +1) x + λ ∣x∣ + O(x2 , λ x∣x∣−α1 ) (λ - const. of integration). (2.23) For α1 ∈ {0, −1, −2, ...} we get a different solution for each choice of α1 . For example, choosing α1 = −1 gives A(x) = −1 −β1 + (β2 + β1 α2 )xln∣x∣ + β1 α2 x + λ x + O(x2 ln∣x∣) (λ - const. of integration). (2.24) However, for all choices of α1 the dominant term in the denominator as x → 0 results in the same outcome that lim A(x) = x→0 { − α1 β1 0 if α1 ≤ 0 or λ = 0 if α1 ≥ 0 and λ ∕= 0 . (2.25) We calculated values for α1 and β1 for each zero of B′ r + 4B given a variety of values of m and l (see figs. 2.4 and 2.5). Having verified our results for a range of choices 0.02 < l < 750, we found that 32 Figure 2.4: α1 and β1 plotted against Figure 2.5: α1 and β1 plotted against m l m l for the first zero of B′ r + 4B. for the second zero of B′ r + 4B. at each zero the values of α1 and β1 depend on m l but not on m and l individually. Knowing that zeros of A(r) correspond to singularities in the space-time, we sought to determine where α1 = 0 in terms of ml . With x = r − r0 , we have B′ r + 4B ≈ x(B′′ r + 5B′ ) zero of B′ r + 4B, r=r0 so α1 = 2B′′ r − (B′ )2 r 4B + 4B′ + B r B′′ r + 5B′ 33 r=r0 in the vicinity of each ′′ ′ = 2B′′ r + 7B′ B r + 5B . (2.26) r=r0 Substituting (2.26) into (2.3) for α1 = 0, we get the following quadratic in r3 r6 + 14ml 2 r3 − 24m2 l 4 = 0 (2.27) 1 with one positive real root corresponding to r ≈ (1.544ml 2 ) 3 . Substituting this into B′ r + 4B = 0 yields m l ≈ 1.202. From (2.18), we know that there will be two zeros of B′ r + 4B for Figure 2.5 shows that α1 > 0 for all m l m l ≈ 1.202. at the second zero, thus at the first zero α1 > 0 when > 1.202 and α1 < 0 when ml < 1.202. The preceeding analysis shows how the value of α1 will help us to categorize the behaviour of A around the zeros of B′ r +4B given ml . As we know that A′ → ±∞ corresponds to singular behaviour in the space-time, we perform a similar analysis of the behaviour of A′ . Differentiating (2.23) gives m l β1 α2 − β2 α1 α1 λ ∣x∣−α1 − + ... α1 (α1 + 1) x A′ (x) = ( )2 β1 β1 α2 − β2 α1 − α1 − x + λ ∣x∣ + ... α1 α1 (α1 + 1) − ⎧ α1 (β2 α1 −β1 α2 ) ⎨ β12 (α1 +1) ′ lim A (x) = −sgn(α1 λ x)∞ x→0 ⎩ 0 (2.28) if α1 < −1 or λ = 0 if −1 < α1 < 1 and λ ∕= 0 (2.29) if α1 > 1 and λ ∕= 0 . We know from (2.25) that we need not consider α1 ≥ 0 (unless λ = 0), but the behaviour of A′ does depend on whether α1 is greater or less than -1 at the zeros of B′ r + 4B. We therefore sought to determine where α1 = −1 in terms of numerically found this corresponds to m l ≈ 1.191. From (2.18), we know that there will be two zeros of B′ r + 4B for Figure 2.5 shows that α1 > −1 for all m l m l m l and ≈ 1.191. at the second zero, thus at the first zero α1 > −1 when > 1.191 and α1 < −1 when ml < 1.191. At each zero, the values of α1 and β1 depend on ml but not on m and l individually. That is not true of the values of α2 and β2 and (2.29) tells us that we may m l 34 need α2 and β2 in order to categorize the behaviour of A around the first zero of B′ r + 4B for 32 27 < m l < 1.191. However, we did determine (again verifying our results for a range of choices 0.02 < l < 750) that at each zero sgn(α2 ), sgn(β2 ) and the ratio αβ22 depend on ml but not on m and l individually. We can therefore say ( ( )) α2 ) ( − α β α β 1 2 1 1 β2 β1 α2 ) = sgn depends on ml but not on m and l that sgn α1 (ββ22(αα1 −+1) β 2 (α +1) 1 1 1 1 individually. The important thing to take away from the preceeding analysis is that the behaviour of A around the zeros of B′ r + 4B can be effectively characterized in terms of ranges of values of m l. We are not concerned with solutions for which there √ exist points where B = 0 (i.e. m l ≥ 27 4 ), so this leaves five ranges of values to consider with respect to the behaviour of A. Three of these correspond to solutions in which there are two zeros of B′ r + 4B and are distinguished by values taken by a specific coefficient α1 in the series approximation (2.21). We refer to these ranges as low, medium and high referring to the relative strength of the corresponding gravitational field. ∙ Two Critical √Points – High Range: There are two zeros of B′ r + 4B for 1.202 < m l < 27 4 . Both zeros have α1 > 0. ∙ Two Critical Points – Medium Range: There are two zeros of B′ r + 4B for 1.191 < α1 > 0. m l < 1.202. The first zero has −1 < α1 < 0 and the second zero has ∙ Two Critical Points – Low Range: There are two zeros of B′ r + 4B for 32 27 < m l < 1.191. The first zero has α1 < −1 and the second zero has α1 > 0. ∙ No Critical Points: There are no zeros of B′ r + 4B for m l < 32 27 ∙ The Weak Field: See section (2.3) With the exception of the weak field analysis (which has already been covered), we look at each of the above ranges in the following sections. 35 2.4.1 Two Critical Points – High Range For 1.202 < m l < √ 27 4 , in (2.23) α1 > 0 at each zero of B′ r + 4B. We therefore considered the following two situations. λ ∕ = 0 at either zero of B ′ r + 4 B : A = 0 at either zero of B′ r + 4B when λ ∕= 0. Our analysis of the Kretschmann scalar (1.10) indicates singular behaviour of the space-time at the spherical hypersuface corresponding to either zero. λ = 0 at the first zero of B ′ r + 4 B : Noting that α1 > 0 and β1 > 0 (see fig. 2.4), at the first zero of B′ r + 4B we have A = − α1 β1 < 0 (see fig. 2.6). Since A(0) = 1 then A(r0 ) = 0 for some r0 > 0 but inside the spherical hypersurface corresponding to the first zero of B′ r + 4B. Referring again to our analysis of the Kretschmann scalar, this indicates singular behaviour in the space-time at r0 . Therefore, there can be no solutions that correspond to gravitational geons for these values of ml . Figure 2.6: A(r) iteratively plotted against radius r for l = 1 and m = 1.25 against a direction field plot. Here, both the convergence to A = 0 and A = of B′ r + 4B can be clearly seen. 36 −α1 β1 ≈ −1.84 at the first zero 2.4.2 Two Critical Points – Medium Range For 1.191 < m l < 1.202, in (2.23) −1 < α1 < 0 at the first zero of B′ r + 4B, whereas α1 > 0 at the second zero. We therefore considered the following situations. λ ∕ = 0 at the first zero of B ′ r + 4 B : A = − α1 β1 and A′ → −sgn(α1 λ x)∞ at the zero which, from (1.9), indicates singular behaviour of the space-time. λ ∕ = 0 at the second zero of B ′ r + 4 B : A = 0 at the second zero indicating singular behaviour in the space-time. λ = 0 at both zeros of B ′ r + 4 B : Noting that α1 < 0, β1 > 0, α2 < 0 and β2 > 0 at the first zero we have that A = − α1 β1 > 0 and A′ = α1 (β2 α1 −β1 α2 ) β12 (α1 +1) second zero, we have α1 > 0, β1 < 0, α2 < 0 and β2 > 0 so that A = A′ = α1 (β2 α1 −β1 α2 ) β12 (α1 +1) < 0. At the − α1 β1 > 0 and > 0. There is therefore no singular behaviour and no horizon indicated at either zero. However, direction field plots indicate that solutions with A(0) = 1 have A′ > 0 as the first zero of B′ r + 4B is approached. Iterating backward from A = second zero results in A′ − α1 β1 at the > 0 as A approaches the first zero (see fig.2.7). Figure 2.7: Reverse iteration of A(r) plotted against radius r for l = 1 and m = 1.192 from the second zero of B′ r + 4B to the first, against a direction field plot. The solution is incompatible with the analytical result that A′ < 0 at the first zero for λ = 0 solutions. 37 It is apparent that λ = 0 solutions are inconsistent with our boundary conditions. We must therefore conclude that there are no solutions that correspond to gravitational geons for these values of ml . 2.4.3 Two Critical Points – Low Range m ′ For 32 27 < l < 1.191, in (2.23) α1 < −1 at the first zero of B r + 4B, whereas α1 > 0 at the second zero. We therefore examined the behaviour of A at each zero depending on the value of the constant of integration λ . β1 α2 ) > Any choice of λ at the first zero of B ′ r + 4 B : A = −βα1 > 0 and A′ = α1 (ββ22(αα1 −+1) 1 1 1 0 at the first zero of B′ r + 4B. There is therefore no singular behaviour and no horizon indicated. λ ∕ = 0 at the second zero of B ′ r + 4 B : A = 0 at the second zero indicating singular behaviour in the space-time. λ = 0 at the second zero of B′ r + 4 B: A = − α1 β1 > 0 and A′ = α1 (β2 α1 −β1 α2 ) β12 (α1 +1) > 0 at the second zero. No singular behaviour or horizon is indicated. Our analysis shows the possibility of solutions corresponding to gravitational geons with λ = 0 at the second zero, but instability at each zero prevents us from confirming any particular solution using iterative methods (see fig.2.8). α1 (β2 α1 − β1 α2 ) −α1 However, we have been able to confirm that lim A = and lim A′ = x→0 x→0 β1 β12 (α1 + 1) around each zero for each section constructed numerically (see figs. 2.9, 2.10 and 2.11). We piecewise constructed A around each zero of B′ r +4B using the numerically constructed sections (see fig.2.12). The resulting solution is sewn together at the spherical hypersurfaces corresponding to each zero. The sewing together of two manifolds may induce a surface energy-momentum tensor (1.13) on the surface where they are joined. The extrinsic curvature Kµν depends on the metric and its first derivatives [4]. Given that A, B and B′ are all continuous in our piecewise constructed solution, there will be a non-vanishing Sµν only if A′ is discontinuous across either hypersuface corresponding to a zero of B′ r + 4B. Our analysis shows that A′ is continuous, thus implying that Sµν = 0. It should be noted that there could be contributions to Sµν from the Weyl term 38 Figure 2.8: This direction field plot of A(r) against radius r shows that A(r) is numerically unstable at the zeros of B′ r + 4B. Here we choose l = 1 and m = 1.189 such that the zeros are found at r = 1.281, 1.392. (1.26). If that were the case, Sµν = 0 in our piecewise constructed solution would imply that there is a localized gravitational source (such as matter) at the hypersurface and the stresses contributed by each cancel one another exactly. As the space-time off the brane is not known, it is not possible to check to see if such contributions exist. If there are no contributions to Sµν from the Weyl tensor, the piecewise constructed solutions correspond to gravitational geons. 2.4.4 No Critical Points For these values of m l, all solutions with A(0) = 1 correspond to gravitational geons. Choosing l = 1 we confirmed the behaviour of A using iterative numerical analysis (see fig.2.13). Also evident is the progression from the weak field solution for m l << 1 (fig.2.3) to the piecewise constructed solution for (fig.2.12). 39 32 27 < m l < 1.191 Figure 2.9: Forward iteration of A(r) plotted from r = 0 to the first zero of B′ r + 4B for l = 1 and m = 1.189. A → −α1 β1 ≈ 1.38 as the solution approaches the first zero. Figure 2.10: Reverse iteration of A(r) plotted from the second zero of B′ r + 4B to the first for l = 1 and m = 1.189. 40 Figure 2.11: Forward iteration of A(r) plotted from the second zero of B′ r + 4B for l = 1 and m = 1.189. Figure 2.12: The function A(r) plotted against radius r for l = 1 and m = 1.189. A fourthorder Runge-Kutta numerical iteration was unable to navigate the zeros of B′ r + 4B (at r = 1.281, 1.392), so the result is constructed piecewise between these points. 41 Figure 2.13: Function A(r) plotted against radius r. Here we choose l = 1 and chart a variety of values for m in the specified region. A fourth-order Runge-Kutta numerical iteration is used to generate each result. 42 Chapter 3 Discussion And Conclusion In this paper, we examined solutions to R = 0 in 3+1 dimensional general relativity that correspond to static, spherically symmetric vacuum solutions on the brane in 5 dimensional Randall-Sundrum brane-world cosmologies. We investigated these solutions for behaviour consistent with that of gravitational geons. In the absence of a general solution, we investigated the behaviour of the space-time in the weak field approximation and used a variety of numerical and analytical methods to discover the behaviour for strong field solutions. We were able to ascertain that, for a particular range of parameter choices, gravitational geons are possible. 3.1 The Weak Field In section (2.3), we found that, for metrics of the form B(r) = 1 + b(r) and A(r) = 1 + a(r) in (2.1) with weak field approximations with ∣b(r)∣ << 1 and ∣a(r)∣ << 1 for all r ≥ 0, gravitational geons will exist as long as ∣B′ r∣ << 1 for all r ≥ 0. The main result of this thesis, that gravitational geons are possible in Randall- Sundrum brane-worlds, is achieved in the weak field, independent of the choice of B(r). However, in order to quantitatively describe the potentially measurable presence of a gravitational geon, a specific choice of B(r) must be made. We chose a specific B(r) and examined the corresponding space-time solution 43 numerically for both weak and strong fields. Our numerical analysis was largely done using a custom developed software application, the implementation of which was validated in part by the very close agreement between numerical results and analytical results for weak fields (see fig. 2.3). Further analysis of the weak field solution could potentially show that some choices of B(r) are viable whereas others are not. This analysis would involve evolving the solution into the bulk (to be discussed in section (3.3)), and having done that, determining the stability of the solution with respect to perturbations. Without making weak field approximations, this work must be done numerically for each specifically chosen B(r). For the analytical weak field solution, it is possible that the evolution into the bulk could be performed without resorting to numerical methods, but the degree of inaccuracy that would be introduced due to approximation is unknown. 3.2 Identifying m and l This thesis shows that in general for brane-world space-times, static, spherically symmetric and topologically trivial gravitational geons are possible. However, for the specific B(r) that we chose for our strong field analysis (2.3), gravitational geons are only viable for particular choices of the parameters m and l. The choice of B(r) is by no means unique, but it is worth looking at what these parameters represent physically as an example of how a gravitational geon solution would manifest itself as a physical entity. 2m ). In the r→∞ r limit as r → ∞ the gravitational field is weak and static. If we make the further B(r) is Schwarzschild in the limit as r → ∞ (i.e. lim B(r) = 1 − assumption that any test particles are travelling at sub-relativistic speeds, we find that gtt = −(1 + 2Φ) where Φ corresponds to the Newtonian gravitational poten- tial. Thus knowing B(r) is sufficient for us to identify the parameter m with the gravitational mass of the space-time solution. To get a sense of the physical significance of the parameter l it is useful to consider the metric of de Sitter space in static coordinates, 44 ( ) ( ) Λ 2 Λ 2 −1 2 2 ds = − 1 − r dt + 1 − r dr + r2 (d θ 2 + sin2 θ d φ 2 ), 3 3 2 (3.1) where Λ is a cosmological constant. For our choice of B(r) we have that B ≈ ( )1 ( )2 1 − rl for r << ml 2 3 . Thus, l has the effect of controlling the strength of a positive cosmological constant like term (with Λ = 3 ) l2 that appears in the limit as r → 0. Figure 3.1: B(r) plotted for m = 1 and 1 ≤ l ≤ 33.3. Note that all solutions have m l < 32 27 so these solutions all correspond to gravitational geons. To see how l affects the space-time away from the radial center, see figure 3.1 in which B is shown for m = 1 and varying l. For a particular choice of m there will be solutions (via appropriate choice of l) that correspond to gravitational geons that 32 are radially concentrated (as ml → 27 ) or as spread out as desired (as ml → 0). If we choose B′ (r) = 0 as a radial marker with which to assign a radial scale for our gravitational geon solutions, then observing that B′ (r) = 2mr4 − 8m2 l 2 r (r3 + 2ml 2 )2 we have B′ (r) = 0 at 45 (3.2) 1 r = (4ml 2 ) 3 . (3.3) We have shown that there are no solutions for the candidate space-times corresponding to gravitational geons for B′ (r) = 0 will occur is m l > 1.191. Therefore, the lower limit on where r = 1.78 Gm , c2 (3.4) which means that in the lower limit, the radial scale of the candidate space-time solutions is of the same order as that of the Schwarzschild solution. 3.3 Suggestions for Further Analysis The conclusion that static, spherically symmetric and topologically trivial gravitational geons are possible in brane-world cosmologies has been achieved in the ideal circumstance that the 3+1 dimensional space-time on the brane is an unperturbed vacuum. We have shown that such solutions (unlike the case for 3+1 dimensional general relativity) are possible. However, further analysis is necessary to determine if the candidate space-times are likely to actually occur should our universe have a Randall-Sundrum brane-world type cosmological structure. In order to test the dynamical behaviour of our candidate space-times, a logical next analytical step would be to test for stability against small perturbations. In this thesis, we have looked in detail at the 3+1 dimensional space-time on the brane. We are confident that solving R = 0 for our candidate space-times will give rise to a vacuum brane-world solution in the 4+1 dimensional bulk [21], but the effective field equations on the brane are not a closed system. In order to test the stability of the candidate space-times, our solution must be evolved into the bulk. The projection of the bulk Weyl tensor, which factors into the induced 3+1 dimensional field equations on the brane (1.34), is governed by a coupled system of 4+1 dimensional equations obtained from the bulk Bianchi identities [15] (round brackets around indices denote symmetrization whereas square brackets denote anti-symmetrization) 46 1 Ln EAB = ∇C BC(AB) + κ52 Λ(KAB − gAB K) + KCD(4) RCADB 6 + 3KC (A EB)C − KEAB + (KAC KDB − KABKCD )KCD , Ln BABC = −2∇[A EB]C + KC D BABD − 2BCD[A KB]D , Ln (4) RABCD = −2RABE[C KD]E − ∇ABCDB + ∇ABDCB , (3.5) in which EAB is the electric part of the bulk Weyl tensor (1.26) whose projection comes into the induced Einstein field equations on the brane 1.32, KAB is the extrinsic curvature (1.14) (K = KA A ) and Bµνσ = (5)CABCD nD gµ A gν B gσ C (3.6) is the magnetic part of the bulk Weyl tensor. In (3.6), Bµνσ is expressed as a 3+1 dimensional tensor projected onto the brane from the 4+1 dimensional Weyl tensor, whereas in (3.5) BABC is expressed as a 4+1 dimensional tensor. When evolving a solution from the brane into the bulk, the 3+1 dimensional tensors B and E are equated with their 4+1 dimensional counterparts by evaluating the bulk tensors on hypersurfaces in the limit approaching the brane. We therefore define boundary conditions for B and E at the brane in order to solve (3.5) as ∇A EAB = κ54 ∇A SAB , ) ( 1 BABC = 2∇[A KB]C = κ52 ∇[A TB]C − gB]C T . 3 (3.7) We have R = 0, so for our candidate space-times, the induced field equations on the brane (1.34) become Eµν = −Rµν . (3.8) For a space-time metric of the form (2.1), the components of the Ricci tensor are [5] 47 Rtt = Rrr = Rθ θ = Rφ φ = ( ) −B′′ B′ A′ B′ B′ + + − 2A 4A rA ( ′ A ′ B) ′′ ′ ′ B B A B A − + − 2B 4B ( A B ) rA 1 r B′ A′ − − −1 + 2A B A A Rθ θ sin2 θ , (3.9) which for the weak field limit of our candidate space-times become Rtt Rrr Rθ θ Rφ φ −24m2 l 2 (−r3 + ml 2 ) (r3 + 2ml 2 )3 3m(r6 − 18r3 ml 2 + 8m2 l 4 ) = (r3 + 2ml 2 )3 2 3mr (−r6 − 6r3 ml 2 + 16m2 l 4 ) = 2(r3 + 2ml 2 )3 2 = Rθ θ sin θ . = (3.10) It is worth noting that, for strong fields, we would obtain these components numerically. In either case, the components of the Ricci tensor (and therefore E ) are non-zero in general. When a perturbation is introduced, discovering the dynamical behaviour of E (and hence the space-time both on and off the brane) depends on solving the system of equations (3.5). The intention of performing this analysis would be to determine if there are any specific choices of the parameters m and l for which the candidate space-times are dynamically stable. Of course, in the event that none of the gravitational geon solutions that we have investigated here are stable, it is understood that the form of B(r) that we have investigated is by no means a unique choice. Although there has been much research published regarding brane-world cosmologies in recent years, there are other modifications to the 3+1 dimensional Einstein field equations that may more accurately represent the universe in which we live. One such class of theories are known as f(R) theories (see [43] for a review). These theories are obtained when one replaces the Ricci scalar R in the Lagrangian 48 formulation of general relativity with a function of the Ricci scalar f (R). The Einstein-Hilbert action SEH = √ d 4 x −gR (3.11) √ d 4 x −g f (R). (3.12) 1 16π G ∫ becomes 1 S= 16π G ∫ Depending on the definition of f (R), the modified Einstein field equations may or may not admit space-time solutions corresponding to gravitational geons similar to those we have investigated here. 3.4 Significance of the Result There is currently a great deal of interest in the Randall-Sundrum brane-world, which models the 3 spatial dimensions of our universe as a brane embedded in a 4+1 dimensional space-time. In 3+1 dimensional general relativity, there are no static, spherically symmetric, topologically trivial vacuum space-time solutions other than flat, empty space. In this thesis, we have shown that the same is not true for Randall-Sundrum brane world models. On the brane, there do exist static, spherically symmetric, topologically trivial vacuum space-time solutions that have mass. The source of mass for these gravitational geons derives from the bulk Weyl tensor via its contribution to the effective Einstein field equations on the brane. Further analysis is required to show whether the candidate space-times presented here (or other candidate space-times) are likely to occur and persist in a brane-world teeming with matter and energy. Insofar as gravitational geons only interact gravitationally, they would only be observable indirectly. The possibility that gravitational geons may represent a form of dark matter has been raised before (see for example [39, 44]). The significance of the result of this thesis is therefore that, should our universe be a Randall-Sundrum brane-world, gravitational geons represent a viable explanation for dark matter observations. This result is also significant insofar as it adds to a body of work that shows 49 how space-time solutions differ between standard general relativity and RandallSundrum brane-worlds. It remains to be seen whether this body of work will ultimately provide evidence for or against the Randall-Sundrum brane-world model as an effective model of our universe. 50 Bibliography [1] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999) [arXiv:hepph/9905221]. → pages 1, 7 [2] J. A. Wheeler, Phys. Rev. 97, 511 (1955). → pages 1, 12 [3] C. W. Misner, K. S. Thorne and J. A. Wheeler, Gravitation, (Freeman, S. Francisco, 1973), pp. 551–555 → pages 2, 3, 5 [4] S. M. Carroll, Spacetime and Geometry, (Addison Wesley, S. Francisco, 2004), pp. 443–452 → pages 6, 38 [5] S. Weinberg, Gravitation and Cosmology: Principles and Applications of the General Theory of Relativity, (John Wiley and Sons, Inc., 1972), pp. 172–179 → pages 2, 47 [6] K. A. Bronnikov and S. W. Kim, Phys. Rev. D 67, 064027 (2003) [arXiv:grqc/0212112]. → pages 4, 19, 28 [7] T. Padmanabhan, Gravitation: Foundation and Frontiers, (Cambridge University Press, New York, 2010), pp. 552–554 → pages 5 [8] W. Israel, Nuovo Cim. B 44S10, 1 (1966) [Erratum-ibid. B 48, 463 (1967)] [Nuovo Cim. B 44, 1 (1966)]. → pages 5 [9] T. Kaluza, Sitzungsber. Preuss. Akad. Wiss. Berlin (Math. Phys. ) 1921, 966 (1921). → pages 6, 13 [10] O. Klein, Z. Phys. 37, 895 (1926) [Surveys High Energ. Phys. 5, 241 (1986)]. → pages 6 [11] J. M. Overduin and P. S. Wesson, Phys. Rept. 283, 303 (1997) [arXiv:grqc/9805018]. → pages 6 51 [12] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 (1999) [arXiv:hepth/9906064]. → pages 7 [13] T. Shiromizu, K. Maeda and M. Sasaki, Phys. Rev. D 62, 024012 (2000) [arXiv:gr-qc/9910076]. → pages 7 [14] R. Maartens and K. Koyama, arXiv:1004.3962 [hep-th]. → pages 7 [15] M. Sasaki, T. Shiromizu and K. Maeda, Phys. Rev. D 62, 024008 (2000) [arXiv:hep-th/9912233]. → pages 7, 46 [16] D. Langlois, Int. J. Mod. Phys. A 17, 2701 (2002) [arXiv:gr-qc/0205004]. → pages [17] R. Maartens, Frames and Gravitomagnetism, ed. J Pascual-Sanchez et al. (World Sci., 2001) pp 93-119 arXiv:gr-qc/0101059. → pages [18] E. Papantonopoulos, Lect. Notes Phys. 592, 458 (2002) [arXiv:hepth/0202044]. → pages [19] D. N. Vollick, Gen. Rel. Grav. 34, 471 (2002) [arXiv:hep-th/0005033]. → pages 17 [20] D. Wands, Class. Quant. Grav. 19, 3403 (2002) [arXiv:hep-th/0203107]. → pages 7 [21] N. Dadhich, R. Maartens, P. Papadopoulos and V. Rezania, Phys. Rev. D 62, 024008 (2000) [arXiv:hep-th/9912233]. → pages 12, 17, 46 [22] T. Regge and J. A. Wheeler, Phys. Rev. 108, 1063 (1957). → pages 13 [23] A. Einstein, mass free of singularities.,” Revista Universidad Nacional de Tucuman A2, 11 (1941). → pages 13 [24] A. Einstein,W. Pauli, Ann. of Math 44, 2 (1945). → pages 13, 15 [25] R. Serini, Rend. Lincei 27, 235 (1918). → pages 13 [26] A. Lichnerowicz, Theories Relativiste de la Gravitation et de lElectromagnetisme, (Masson, 1955) → pages 13 [27] J. van Dongen, Stud. Hist. Philos. Mod. Phys. 33, 185 (2002) [arXiv:grqc/0009087]. → pages 13 [28] G. D. Birkhoff, Relativity and Modern Physics, (Harvard University Press, Cambridge, MA, 1923) → pages 13 52 [29] A. Chamblin, S. W. Hawking and H. S. Reall, Phys. Rev. D 61, 065007 (2000) [arXiv:hep-th/9909205]. → pages 16 [30] S. Shankaranarayanan and N. Dadhich, Int. J. Mod. Phys. D 13, 1095 (2004) [arXiv:gr-qc/0306111]. → pages 18 [31] E. Rodrigo, arXiv:0708.0045 [gr-qc]. → pages 19 [32] N. Dadhich, S. Kar, S. Mukherji and M. Visser, Phys. Rev. D 65, 064004 (2002) [arXiv:gr-qc/0109069]. → pages 20 [33] D. N. Vollick, Gen. Rel. Grav. 34, 1 (2002) [arXiv:hep-th/0004064]. → pages 20 [34] D. R. Brill and J. B. Hartle, Phys. Rev. 135, B271 (1964). → pages 20 [35] F. I. Cooperstock, V. Faraoni and G. P. Perry, Mod. Phys. Lett. A 10, 359 (1995). → pages 20 [36] F. I. Cooperstock, V. Faraoni and G. P. Perry, Int. J. Mod. Phys. D 5, 375 (1996) [arXiv:gr-qc/9512025]. → pages 20 [37] R. D. Sorkin, Topological Properties and Global Structure of Space-Time, volume B138 of NATO Advanced Study Institutes Series, 249, edited by P.G. Bergmann and V. De Sabbata, D. Reidel Publishing Company, DordrechtHolland, 1986. → pages 20 [38] I. Dymnikova and E. Galaktionov, Class. Quant. Grav. 22, 2331 (2005) [arXiv:gr-qc/0409049]. → pages 21 [39] D. N. Vollick, Class. Quant. Grav. 25, 175004 (2008) [arXiv:0807.0611 [grqc]]. → pages 21, 49 [40] C. Teitelboim, Phys. Lett. B 126 (1983) 41. → pages 21 [41] R. Jackiw, Nucl. Phys. B 252, 343 (1985). → pages 21 [42] D. N. Vollick, Class. Quant. Grav. 27, 169701 (2010) [arXiv:1004.1401 [grqc]]. → pages 21 [43] T. P. Sotiriou and V. Faraoni, Rev. Mod. Phys. 82, 451 (2010) [arXiv:0805.1726 [gr-qc]]. → pages 48 [44] R. A. Sones, arXiv:gr-qc/0506011. → pages 49 53
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Gravitational geons on the brane
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Gravitational geons on the brane Kermode, Daniel John 2011
pdf
Page Metadata
Item Metadata
Title | Gravitational geons on the brane |
Creator |
Kermode, Daniel John |
Publisher | University of British Columbia |
Date | 2011 |
Date Issued | 2011-04-18 |
Description | Gravitational geons, if they exist, represent a candidate for cosmological dark matter observations. Gravitational Geons that are nonsingular, asymptotically flat, topologically trivial spacetime solutions are not possible in General Relativity. The purpose of this thesis is to investigate whether such geons are possible on a 3 dimensional brane embedded in a 4+1 dimensional spacetime. Our approach is to identify the characteristics of a candidate spacetime solution on the brane, choose an example spacetime that exhibits those characteristics and determine if it corresponds to a gravitational geon. |
Genre |
Thesis/Dissertation |
Type |
Text |
Language | Eng |
Collection |
Electronic Theses and Dissertations (ETDs) 2008+ |
Date Available | 2011-04-18 |
Provider | Vancouver : University of British Columbia Library |
DOI | 10.14288/1.0071701 |
Degree |
Master of Science - MSc |
Program |
Interdisciplinary Studies |
Affiliation |
Irving K. Barber School of Arts and Sciences (Okanagan) |
Degree Grantor | University of British Columbia |
Graduation Date | 2011-05 |
Campus |
UBCO |
Scholarly Level | Graduate |
URI | http://hdl.handle.net/2429/33721 |
Aggregated Source Repository | DSpace |
Download
- Media
- ubc_2011_spring_kermode_danny.pdf [ 860.01kB ]
- Metadata
- JSON: 1.0071701.json
- JSON-LD: 1.0071701+ld.json
- RDF/XML (Pretty): 1.0071701.xml
- RDF/JSON: 1.0071701+rdf.json
- Turtle: 1.0071701+rdf-turtle.txt
- N-Triples: 1.0071701+rdf-ntriples.txt
- Original Record: 1.0071701 +original-record.json
- Full Text
- 1.0071701.txt
- Citation
- 1.0071701.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Country | Views | Downloads |
---|---|---|
United States | 15 | 0 |
China | 9 | 20 |
Japan | 5 | 0 |
Canada | 2 | 1 |
France | 1 | 0 |
Sweden | 1 | 0 |
Finland | 1 | 0 |
City | Views | Downloads |
---|---|---|
Ashburn | 12 | 0 |
Beijing | 9 | 0 |
Tokyo | 5 | 0 |
Sunnyvale | 2 | 0 |
Vernon | 2 | 1 |
Stockholm | 1 | 0 |
Redmond | 1 | 0 |
Klaukkala | 1 | 0 |
Sollies-Ville | 1 | 0 |
{[{ mDataHeader[type] }]} | {[{ month[type] }]} | {[{ tData[type] }]} |
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.24.1-0071701/manifest