The Hawking mass for ellipsoidal 2-surfaces in Minkowski and Schwarzschild spacetimes Daniel Hansevi

Size: px
Start display at page:

Download "The Hawking mass for ellipsoidal 2-surfaces in Minkowski and Schwarzschild spacetimes Daniel Hansevi"

Transcription

1 Examensarbete The Hawking mass for ellipsoidal -surfaces in Minkowski and Schwarzschild spacetimes Daniel Hansevi LiTH - MAT - EX / SE

2

3 The Hawking mass for ellipsoidal -surfaces in Minkowski and Schwarzschild spacetimes Applied Mathematics, Linköpings Universitet Daniel Hansevi LiTH - MAT - EX / SE Examensarbete: 30 hp Level: D Supervisor: Göran Bergqvist, Applied Mathematics, Linköpings Universitet Examiner: Göran Bergqvist, Applied Mathematics, Linköpings Universitet Linköping: June 008

4

5 Avdelning, Institution Division, Department Matematiska Institutionen LINKÖPING SWEDEN Datum Date June 008 Språk Language Svenska/Swedish Rapporttyp Report category Licentiatavhandling ISBN ISRN x x LiTH - MAT - EX / SE Engelska/English Examensarbete C-uppsats Serietitel och serienummer ISSN D-uppsats Övrig rapport Title of series, numbering URL för elektronisk version Titel Title The Hawking mass for ellipsoidal -surfaces in Minkowski and Schwarzschild spacetimes Författare Author Daniel Hansevi Sammanfattning Abstract In general relativity, the nature of mass is non-local. However, an appropriate definition of mass at a quasi-local level could give a more detailed characterization of the gravitational field around massive bodies. Several attempts have been made to find such a definition. One of the candidates is the Hawking mass. This thesis presents a method for calculating the spin coefficients used in the expression for the Hawking mass, and gives a closed-form expression for the Hawking mass of ellipsoidal -surfaces in Minkowski spacetime. Furthermore, the Hawking mass is shown to have the correct limits, both in Minkowski and Schwarzschild, along particular foliations of leaves approaching a metric -sphere. Numerical results for Schwarzschild are also presented. Nyckelord Keyword Hawking mass, Quasi-local mass, General relativity, Ellipsoidal surface.

6 vi

7 Abstract In general relativity, the nature of mass is non-local. However, an appropriate definition of mass at a quasi-local level could give a more detailed characterization of the gravitational field around massive bodies. Several attempts have been made to find such a definition. One of the candidates is the Hawking mass. This thesis presents a method for calculating the spin coefficients used in the expression for the Hawking mass, and gives a closed-form expression for the Hawking mass of ellipsoidal -surfaces in Minkowski spacetime. Furthermore, the Hawking mass is shown to have the correct limits, both in Minkowski and Schwarzschild, along particular foliations of leaves approaching a metric -sphere. Numerical results for Schwarzschild are also presented. Keywords: Hawking mass, Quasi-local mass, General relativity, Ellipsoidal surface. Hansevi, 008. vii

8 viii

9 Contents 1 Introduction Background Purpose Chapter outline Mathematical preliminaries 3.1 Manifolds Foliations Tangent vectors forms Tensors Abstract notation Component notation Tensor algebra Tensor fields Abstract index notation Metric Curvature Covariant derivative Metric connection Parallel transportation Curvature Geodesics Tetrad formalism Newman-Penrose formalism Spin coefficients General relativity Solutions to the Einstein field equation Minkowski spacetime Schwarzschild spacetime Mass in general relativity Gravitational energy/mass Non-locality of mass Total mass of an isolated system Asymptotically flat spacetimes ADM and Bondi-Sachs mass Hansevi, 008. ix

10 x Contents 4.3 Quasi-local mass The Hawking mass Definition Interpretation Method for calculating spin coefficients Hawking mass for a -sphere In Minkowski spacetime In Schwarzschild spacetime Results Hawking mass in Minkowski spacetime Closed-form expression of the Hawking mass Limit when approaching a metric sphere Limit along a foliation Hawking mass in Schwarzschild spacetime Limit along a foliation Numerical evaluations Discussion Conclusions Future work A Maple Worksheets 37 A.1 Null tetrad A. Minkowski A.3 Schwarszchild

11 List of Figures.1 A manifold M and two overlapping coordinate patches A foliation of a manifold M Intuitive picture of a tangent space Picture of a 1-form Null cone Parallel transportation in a plane and on the surface of a sphere Deviation vector y a between two nearby geodesics λ s and λ s Spin coefficients Oblate spheroid m H (S 1 ) plotted against parameter ξ A curve given by r = 1 + ε sin θ m H ( S r ) plotted against 4 r 60 for some values of ε m H ( S r ) plotted against.3 r 16 for some values of ε m H ( S r ) plotted against.3 r 16 for ε = ω/r Hansevi, 008. xi

12 xii List of Figures

13 Chapter 1 Introduction 1.1 Background In general relativity, the nature of the gravitational field is non-local, and therefore the gravitational field energy/mass cannot be given as a pointwise density. However, there might be possible to find a satisfying definition of mass at a quasi-local level, that is, for the mass within a compact spacelike -surface. Several attempts have been made, but the task has proven difficult, and there is still no generally accepted definition. One of the candidates for a description of quasi-local mass originates from a paper about gravitational radiation written by Stephen Hawking [5] in The Hawking mass can be viewed as a measure of the bending of outgoing and ingoing light rays orthogonal to the surface of a spacelike -sphere, and it has been shown to have various desirable properties [13]. To reach an appropriate definition for quasi-local mass would certainly be of great value. It could give a more detailed characterization of the gravitational field around massive bodies, and it should be helpful for controlling errors in numerical calculations [13]. 1. Purpose Even though Hawking s expression was given for the mass contained in a spacelike -sphere, it can be calculated for a general spacelike -surface. In this thesis we will calculate the Hawking mass for spacelike ellipsoidal -surfaces, both in flat Minkowski spacetime and curved Schwarzschild spacetime. Hansevi,

14 Chapter 1. Introduction 1.3 Chapter outline Chapter This chapter provides a short introduction of the mathematics needed. We introduce the notion of a manifold that is used for the model of curved spacetime in general relativity. Structure is imposed on the manifold in the form of a covariant derivative operator and a metric tensor. The concept of geodesics as curves that are as straight as possible is introduced along with the definition of the Riemann curvature tensor which is a measure of curvature. We end the chapter with a look at the Newman-Penrose formalism and the spin coefficients which are used in the definition of the Hawking mass. Chapter 3 A very brief presentation of general relativity is given followed by the two solutions (of the field equation) of particular interest in this thesis Minkowski spacetime and Schwarzschild spacetime. Chapter 4 Discussion of mass in general relativity, and why it cannot be localized. Chapter 5 A closer look at the Hawking mass definition and interpretation. A method for calculating the spin coefficients used in the expression for the Hawking mass is presented. Chapter 6 In this chapter, we calculate the Hawking mass for ellipsoidal - surfaces in both Minkowski spacetime and Schwarzschild spacetime. The main result is the closed-form expression for the Hawking mass of an -ellipsoid in Minkowski spacetime. Furthermore, some limits of the Hawking mass are proved. Chapter 7 Conclusions and future work. Appendix Maple worksheets used for some of the calculations performed in chapter 7.

15 Chapter Mathematical preliminaries This chapter provides a short introduction of the mathematics needed for a mathematical formulation of the general theory of relativity..1 Manifolds In general relativity, spacetime is curved in the presence of mass, and gravity is a manifestation of curvature. Thus, the model of spacetime must be sufficiently general to allow curvature. An appropriate model is based on the notion of a manifold. A manifold is essentially a space that is locally similar to Euclidean space in that it can be covered by coordinate patches. Globally, however, it may have a different structure, for example, the two-dimensional surface of a sphere is a manifold. Since it is curved, compact and has finite area its global properties are different from those of the Euclidean plane, which is flat, non-compact and has infinite area. Locally, however, they share the property of being able to be covered by coordinate patches. As a mathematical structure a manifold stands on its own, but since it can be covered by coordinate patches, it can be thought of as being constructed by gluing together a number of such patches. A C k n-dimensional manifold M is a set M together with a maximal C k atlas {U α φ α }, that is, the collection of all charts (U α, φ α ) where φ α are bijective maps from subsets U α M to open subsets φ α (U α ) R n such that: (i) {U α } cover M, that is, each element in M lies in at least one U α ; (ii) if U α U β is non-empty, then the transition map (see figure.1, page 4) φ β φ 1 α : φ α (U α U β ) φ β (U α U β ), (.1) is a C k map of an open subset of R n to an open subset of R n. Each (U α, φ α ) is a local coordinate patch with coordinates x α (α = 1,..., n) defined by φ α. In the overlap U α U β of two coordinate patches (U α, φ α ) and (U β, φ β ), the coordinates x α are C k functions of the coordinates x β, and vice versa. M is said to be Hausdorff 1 if for every distinct p, q M (p q) there exist two subsets U α M and U β M such that p U α, q U β and U α U β =. 1 Felix Hausdorff ( ), German mathematician Hansevi,

16 4 Chapter. Mathematical preliminaries M p Uα φ α x α φ β φ 1 α U β φ β x β R n Figure.1: A manifold M and two overlapping coordinate patches. Furthermore, it is natural to introduce the notion of a function on M and the notion of a curve in M as follows. A (real-valued) function f on a C k manifold M is a map f: M R. It is said to be of class C r (r k) at p M if the map f φ 1 α : R n R in any coordinate patch (U α, φ α ) holding p is a C r function of the coordinates at p. A C k curve in a manifold M is a map λ of an open interval I R M such that for any coordinate patch (U α, φ α ), the map φ α λ: I R n is C k. Something that is C is usually called smooth. Accordingly, we call a C manifold a smooth manifold, a C function a smooth function and a C curve a smooth curve.. Foliations A foliation of a manifold is a decomposition of the manifold into submanifolds. These submanifolds are required to be of the same dimension, and fit together in a nice way. More precisely [7], a foliation of codimension m of an n-dimensional manifold M is a decomposition of M into a union of disjoint connected subsets {L a } a A, called the leaves of the foliation, with the property that for every point p M there is a coordinate patch (U α, φ α ) holding p such that for each leaf L a, φ α : (U α L a ) (x 1,..., x m, x m+1,..., x n ) R n, (.) where x m+1,..., x n are constants. See figure.. M φ α R n m U α p R m Figure.: A foliation of a manifold M.

17 .3. Tangent vectors 5.3 Tangent vectors A manifold can be curved and therefore has no global vector space structure. There is no natural way to add two points on a sphere and end up with a third point also on the sphere. However, a local vector space structure can be attained. A definition of a vector that only refers to the intrinsic structure of the manifold would be of great value, because such a vector would be independent of an embedding of the manifold in a space of higher dimension. There is a oneto-one correspondence between vectors and directional derivatives in Euclidean space, and since a manifold is locally similar to Euclidean space, a natural definition is provided by the notion of a vector as a differential operator. Let M be an n-dimensional manifold and let F be the collection of all smooth, real-valued functions on M. A tangent vector, or vector for short, X at a point p M, is a map X: F R such that for all f, g F and all α, β R: (i) X(αf + βg) = αx(f) + βx(g) (ii) X(fg) = f(p)x(g) + g(p)x(f) (linear); (Leibniz rule). Let (U, φ) be a coordinate patch, with coordinates x µ, holding p. For µ = 1,..., n define the map X µ : F R by X µ (f) := x µ (f φ 1 ). (.3) φ(p) It is shown in [14], that X 1,..., X n are linearly independent tangent vectors that span an n-dimensional vector space at p. We call this vector space the tangent space at p and denote it by T p (M), or just T p if the manifold is given by the context. The basis {X 1,..., X n } is called a coordinate basis and is usually denoted by { / x 1,..., / x n }. Thus, an arbitrary tangent vector X can be expressed as X = n X µ X µ =: µ=1 n µ=1 X µ x µ, (.4) where (X 1,..., X n ) R n are the components of X with respect to the coordinate basis. A tangent space at a point p in a manifold may be intuitively understood as the limiting space when smaller and smaller neighbourhoods of p are viewed at greater and greater magnification, see figure.3. T p (M) p X M Figure.3: Intuitive picture of a tangent space.

18 6 Chapter. Mathematical preliminaries.4 1-forms Let M be an n-dimensional manifold. Given a point p M, let T p be the tangent space at p. Let Tp be the space of all linear maps ω: T p R. (.5) Tp is called the dual space of T p, or the cotangent space at p, and is a vector space of dimension n. Elements of Tp are called dual vectors or covectors. They are also called 1-forms. The number which ω maps a vector X T p into, is often written as ω, X. If {e 1,..., e n } is a basis in T p, then there exists an associated dual basis {e 1,..., e n } of Tp consisting of 1-forms e 1,..., e n defined by the property e µ, e ν = δ µ ν := { 1, µ = ν 0, µ ν. (.6) To see this, for every ω T p, we define ω µ := ω, e µ for µ = 1,..., n. Let X be an arbitrary vector in T p. Then ω, X = ω, µ X µ e µ = α X µ ω, e µ = µ ω µ X µ = µ,ν ω µ X ν δ µ ν (.6) = ω µ X ν e µ, e ν µ,ν = µ ω µ e µ, ν X ν e ν = µ ω µ e µ, X. (.7) Since X was arbitrary, it follows that ω = µ ω µ e µ. (.8) Thus every 1-form can be written as a linear combination of e 1,..., e n. Given a coordinate basis { / x 1,..., / x n } for T p, the associated dual basis for T p is the basis {dx 1,..., dx n } of the so called coordinate differentials. For a given 1-form ω, there is a subspace of T p defined by all vectors X for which ω, X is constant. Therefore, a 1-form can be pictured as planes, where ω, X is the number of planes that X is piercing, see figure.4. X ω, X = 3.5 p Y ω, Y = 0 Figure.4: Picture of a 1-form.

19 .5. Tensors 7.5 Tensors General relativity is formulated in the language of tensors. Tensors summarize sets of equations succinctly and reveal structure. There are two distinct ways of introducing tensors: the abstract approach and the component approach..5.1 Abstract notation Let M be an n-dimensional manifold and let p be a point in M. The multilinear map S: Tp Tp T p T p R (.9) }{{}}{{} r factors s factors is called a tensor, of type or valence (r, s), at p, or just an (r, s)-tensor at p for short. We can see that a 1-form is a tensor of type (0, 1). Since T p is a finite dimensional vector space, it is (algebraicly) reflexive, and therefore the second (algebraic) dual space Tp is isomorphic to T p. Thus, we can identify every element in Tp with a unique element in T p and we consider a tangent vector as a tensor of type (1, 0)..5. Component notation Let (U, φ) and (U, φ ) be two overlapping coordinate patches holding a point p M, with coordinates related by x µ = x µ (x 1,..., x n ). (.10) An object with components S µ1...µr ν 1...ν s in (U, φ) and S µ 1...µ r ν 1...ν s in (U, φ ) is called an (r, s)-tensor at p under the transformation x µ x µ, if S µ 1...µ r ν 1...ν s = n µ 1,...,µ r,ν 1,...ν s=1 S µ1...µr ν 1...ν s x µ 1 x µ1... xµ r x ν1... xνs. (.11) x µr x ν 1 x ν s A (0, 1)-tensor (1-form) is often called a covariant vector, and a (1, 0)-tensor (tangent vector) is often called a contravariant vector. Since an (r, s)-tensor S depends linearly on its arguments, it is determined by its components S µ1 µr ν 1 ν s with respect to a basis. Suppose that S is a (1, )-tensor and that {e µ }, {e ν } are dual bases. Define the basis component S µ νρ := S(e µ, e ν, e ρ ) R. (.1) Let ω = µ ω µe µ Tp and X = ν Xν e ν, Y = ρ Y ρ e ρ T p. Then it follows that S(ω, X, Y) = S( µ ω µ e µ, ν X ν e ν, ρ = S(e µ, e ν, e ρ ) ω µ X ν Y ρ µ,ν,ρ Y ρ e ρ ) = µ,ν,ρ S µ νρ ω µ X ν Y ρ. (.13) The components of S satisfy the tensor transformation law (.11).

20 8 Chapter. Mathematical preliminaries.5.3 Tensor algebra Let M be an n-dimensional manifold and let p be a point in M. Assume that S and T are (r, s)-tensors at p and that S is an (r, s )-tensor at p. Furthermore, assume that ω i T p and X j T p for i = 1,..., r + r and j = 1,..., s + s. Addition of tensors (of the same type at the same point) and multiplication of a tensor by a scalar α R, are defined in the obvious way: (S + T)(ω 1,..., ω r, X 1,..., X s ) = S(ω 1,..., ω r, X 1,..., X s ) + T(ω 1,..., ω r, X 1,..., X s ); (.14) (αs)(ω 1,..., ω r, X 1,..., X s ) = α S(ω 1,..., ω r, X 1,..., X s ). (.15) With addition and scalar multiplication defined as above, the space of all (r, s)-tensors at p M, constitutes a vector space of dimension n r+s. The outer product of S and S, denoted by S S, is the (r +r, s+s )-tensor defined by, (S T)(ω 1,..., ω r+r, X 1,..., X s+s ) = S(ω 1,..., ω r, X 1,..., X s ) S (ω r+1,..., ω r+r, X s+1,..., X s+s ). (.16) The contraction with respect to the ith (1-form) and j th (tangent vector) slots is a map from an (r, s)-tensor to an (r 1, s 1)-tensor defined by, (CjS)(ω i 1,..., ω i 1, ω i+1,..., ω r ; X 1,..., X j 1, X j+1,..., X s ) = n S(..., }{{} e k,... ;..., e k,...), (.17) k=1 ith slot }{{} jth slot where {e k } and {e k } are dual bases of T p and T p, respectively..5.4 Tensor fields It is natural to define a C k tensor field of type (r, s) on a manifold M as an assignment of an (r, s)-tensor at each p M such that the components with respect to any coordinate basis are C k functions. We call a C tensor field a smooth tensor field. A vector field is a tensor field of type (1, 0), and a 1-form field is a tensor field of type (0, 1)..5.5 Abstract index notation Equations for tensor components with respect to a particular basis may only be valid in that basis. On the other hand, if we do not specify a basis, the equations we write will be true tensor equations, that is, basis-independent equations that will hold between tensors. It is convenient to introduce a notation called abstract index notation [10]. In this notation, an (r, s)-tensor S is written as S a1 ar b 1 b s, where the indices are abstract markers telling us what type of tensor it is. Assume that S is a (1, )-tensor and that T is a (3, )-tensor. In abstract index notation, we write

21 .6. Metric 9 the outer product of S and T as S a bc T def gh, and the contraction of S with respect to the first slots as S a ab. In order to distinguish between tensors written in the abstract index notation and tensors components, we write the indices of the former with lowercase latin letters and indices of the latter with lowercase greek letters, for example, S µ νρ denotes a basis component of the (1, )-tensor S a bc. Given a tensor equation written in the abstract index notation, the corresponding equation (with greek indices) holds for basis components in any basis if a summation over indices that occurs twice in a term, once as a subscript and once as a superscript, is performed..6 Metric A metric g ab on a manifold M is a non-singular symmetric tensor field. Thus, for every tangent space T p of M: (i) g ab u a v b = g ba u a v b for every u a, v b T p (symmetric); (ii) g ab u a v b = 0 for every v b T p implies that u a = 0 (non-singular). The metric has the structure of a (not necessarily positive definite) inner product on every tangent space of the manifold. If g ab u a v b = 0, then the vectors u a and v b are said to be orthogonal. For any vector v a, the metric can be viewed as a linear map g ab v b : T p R, that is, a 1-form. Since g ab is non-singular, there is a one-to-one correspondence between elements of T p and T p. Given a vector v a, we can apply the metric and get the corresponding 1-form g ab v b, usually denoted by v a in order to make the correspondence with v a notationally explicit. Thus, we can raise and lower indices on tensors by the use of the metric. Particularly, we can write the inner product of two vectors u a and v a as g ab u a v b = u a v a. (.18) Assume that we have two coordinate patches overlapping a neighbourhood of a point p M and that their coordinates at p are related by x µ = x µ (x µ ). Then the basis components of g ab are related by (.19), that is, g µ ν = µ,ν x µ x ν g µν. (.19) x µ x ν It is always possible to find an orthonormal basis {v 1 a,..., v n a } such that v ia v j a = ±δ i j. (.0) The number of basis vectors for which (.0) equals 1 and the number of basis vectors for which (.0) equals 1 are basis independent and called the signature of the metric. A metric of signature (+ ) or ( + +) is called Lorentzian. We will follow the Landau-Lifshitz 3 timelike convention [8] and use a metric of signature (+ ) for spacetime. Hendrik Antoon Lorentz ( ), Dutch physicist. 3 Lev Davidovich Landau ( ) and Evgeny Mikhailovich Lifshitz ( ), Russian physicists.

22 10 Chapter. Mathematical preliminaries With a Lorentzian metric on M, all non-zero vectors in T p can be divided into three classes. With our particular choice of signature, a vector v a is said to be timelike if v a v a > 0, null if v a v a = 0, (.1) spacelike if v a v a < 0. Thus, a Lorentzian metric defines a certain structure on each T p, called a null cone; the set of null vectors form what looks like a double cone if we suppress one spatial dimension, see figure.5. future cone null vector timelike vector spacelike vector p past cone Figure.5: Null cone..7 Curvature An intrinsic notion of curvature, that can be applied to any manifold without reference to a higher dimensional space in which it might be embedded, can be defined in terms of parallel transport. If one parallel-transports a vector around any closed path in the plane, the final vector always coincides with the initial vector. However, for a sphere, the final vector does not coincide with the initial vector when carried along the curve shown in figure.6. Based on this, we characterize the plane as flat and the sphere as curved. Once we know how to parallel transport a vector along a curve, we can use this idea to obtain an intrinsic notion of curvature of any manifold. p u q v Figure.6: Parallel transportation in a plane and on the surface of a sphere.

23 .7. Curvature 11 Since the tangent spaces at two distinct points are different vector spaces it is not meaningful to say that a vector in the first tangent space equals a vector in the latter. Thus, before we can define parallel transport, we must impose more structure on the manifold. Given a notion of a derivative operator, it is natural to define a vector to be parallel-transported if its derivative along the given curve is zero. The notion of curvature can be defined in terms of the failure of the final vector to coincide with the initial vector when parallel transported around an infinitesimal closed curve, which in turn corresponds to the lack of commutativety of derivatives..7.1 Covariant derivative A connection, or covariant derivative operator, a on a smooth manifold M assigns to every vector field x a on M, a differential operator x a a that maps an arbitrary vector field y a on M into a vector field x a a y b such that for all vector fields x a, y a, z a and functions f on M: (i) x a a (y b + z b ) = x a a y b + x a a z b (ii) x a a (fy b ) = (x a a f)y b + f(x a a y b ) (iii) x a a f = x(f) (linear); (Leibniz rule); (consistency with the notion of tangent vectors). The vector field x a a y b is called the covariant derivative of y a with respect to x a, and the (1, 1)-tensor a y b mapping x a to x a a y b the covariant derivative of y b. The definition of a can be extended to apply to any tensor field on M by the additional requirement that a when acting on contracted products should satisfy Leibniz rule..7. Metric connection A connection a is not uniquely defined by the above conditions. In [14], however, it is shown that if M is endowed with a metric g ab, then there exists a unique connection a with the properties that for all smooth functions f on M: (i) a g bc = 0 (ii) a b f b a f = 0 (compatible with metric); (torsion-free). This particular a is called the metric connection, or the Levi-Civita 4 connection, on M, and has the properties that for all smooth vector fields y a on M: a y b = a y b + Γ b acy c and a y b = a y b Γ c aby c, (.) where a is an ordinary derivative, and the Christoffel 5 symbol Γ c ab is given by Γ c ab = 1 gcd( a g bd + b g ad d g ab ). (.3) Thus, in a given coordinate patch (U, φ) with coordinates x µ, the coordinate basis components of a y b are given by y ρ x µ + 1 ( g ρσ gνσ x µ + g µσ x ν g ) µν x σ y ν. (.4) σ,ν 4 Tullio Levi-Civita ( ), Italian mathematician. 5 Elwin Bruno Christoffel ( ), German mathematician.

24 1 Chapter. Mathematical preliminaries.7.3 Parallel transportation Given a derivative operator a we can define the notion of parallel transport. A smooth vector field y a is said to be parallelly transported around a curve with tangent vector x a if the equation x a a y b = 0 (.5) is satisfied along the curve. The metric connection has the property that the inner product f = y a z a of any two smooth vector fields y a and z a remains unchanged when parallelly transported along any curve with tangent vector x a, since x a (f) = x a a (g bc y b z c ).7.4 Curvature = (x a a g }{{ bc )y b z c + (x a } a y b )g }{{} bc z c + (x a a z c )g }{{} bc y b. (.6) =0 =0 =0 Let ω a be any smooth 1-form field on M. The Riemann 6 curvature tensor R abc d is the tensor field on M defined by, ( a b b a ) ω c =: R abc d ω d, (.7) that is directly related to the failure of a vector to return to its initial value when parallel transported around a small closed curve. 7 For any smooth vector field t a on M, the corresponding expression is ( a b b a ) t c = R abd c t d. (.8) The Ricci 8 tensor is defined by contraction with respect to the second and fourth slot, R ac := R abc b, and the Ricci scalar curvature is defined by further contraction, R := R a a. The Ricci tensor and scalar curvature are used to define the Einstein tensor, which is a fundamental tensor in general relativity. G ab := R ab 1 R g ab, (.9).7.5 Geodesics Let x a be a smooth vector field on M. From the theory of ordinary differential equations, we know that given a point p M, there exists a unique curve that passes through p and has the property that for each point on the curve, the tangent vector of the curve coincides with the corresponding vector of x a. Such a curve is called an integral curve. Let x a be a smooth vector field such that x a a x b = 0. Then the integral curves of x a are called geodesics. 6 Georg Friedrich Bernhard Riemann ( ), German mathematician. 7 Some authors reverse the sign of the left-hand side in the definition. 8 Gregorio Ricci-Curbastro ( ), Italian mathematician.

25 .7. Curvature 13 Geodesics are curves that are as straight as possible, and it can be shown [14] that there is precisely one geodesic through a given point p M in a given direction x a T p. Let λ s (t) be a one-parameter family of geodesics and consider the twodimensional surface, with coordinates (t, s), spanned by λ s (t). The vector field x a = ( / t) a is tangent to the family of geodesics, thus x a a x b = 0, (.30) and the vector field y a = ( / s) a represents the deviation vector, which is the displacement from the geodesic λ s to an infinitesimally nearby geodesic λ s, see figure.7. λ s λ s x a x a y a Figure.7: Deviation vector y a between two nearby geodesics λ s and λ s. Let f be any smooth function on M. Since it follows that (x a a y b y a a x b ) b f = x a a (y b b f) y a a (x b b f) = x(y(f)) y(x(f)) = f t s f s t = 0, (.31) x a a y b = y a a x b. (.3) The relative acceleration z a, in the direction of y a, of a nearby geodesic when moving along the direction of x a, is given by z a = x c c (x b b y a ) (.3) = x c c (y b b x a ) = (x c c y b ) b x a + y b (x c c b x a ) (.8) = (y c c x b ) b x a + y b (x c b c x a ) R cbd a y b x c x d = y c c (x b b x a ) R cbd a y b x c x d (.5) = R cbd a y b x c x d. (.33) Thus, the geodesic deviation, that is, the acceleration of geodesics toward or away from each other, which is a characterization of the curvature of M, is determined by the Riemann curvature tensor. M is flat if and only if R abc d = 0.

26 14 Chapter. Mathematical preliminaries.8 Tetrad formalism The tetrad formalism (or frame formalism) is a useful technique for deriving useful and compact equations in many applications of general relativity. The idea is to use a so called tetrad basis of four linearly independent vector fields, project the relevant quantities onto the basis, and consider the equations satisfied by them. Let e i a, i = 1,, 3, 4, be smooth vector fields that are linearly independent at each point in spacetime (M, g ab ). Then the Ricci rotation coefficients are defined by γ kij := e k c e j a a e ic. (.34) It is shown in [] that γ ijk = 1 (λ ijk + λ kij λ jki ), where λ ijk = ( b e ja a e jb )e i a e k b, (.35) which is an efficient way of calculating the Ricci rotation coefficients, since there is no need to calculate any covariant derivatives..9 Newman-Penrose formalism The tetrad formalism with the choice of a particular type of null basis, introduced by Ezra Newman and Roger Penrose [9] in 196, is usually called the Newman-Penrose formalism. The basis {l a, n a, m a, m a } consists of null vectors, where l a and n a are real, m a and m a are complex conjugates, satisfying l a l a = n a n a = m a m a = m a m a = 0 (null); l a m a = l a m a = n a m a = n a m a = 0 (orthogonal); l a n a = m a m a = 1 (normalized). (.36).10 Spin coefficients In the Newman-Penrose formalism, the Ricci rotation coefficients are called spin coefficients, and are given in figure.8. κ = γ 311 κ = ν = γ 4 α = γ 14 + γ 344 ρ = γ 314 ρ = µ = γ 43 β = γ 13 + γ 343 σ = γ 313 σ = λ = γ 44 γ = γ 1 + γ 34 τ = γ 31 τ = π = γ 41 ε = γ 11 + γ 341 Figure.8: Spin coefficients. The spin coefficients ρ and ρ are of particular interest to us, since they are used in the definition of the Hawking mass.

27 Chapter 3 General relativity General relativity is the modern geometric theory of space, time and gravitation published by Albert Einstein [4] in In the theory, space and time are unified into spacetime (M, g ab ), which is represented by a smooth four-dimensional Hausdorff manifold M endowed with a Lorentzian metric g ab and a metric connection a. The presence of matter warps spacetime according to the Einstein field equation, G ab := R ab 1 R g ab = 8π T ab, (3.1) where G ab is the Einstein tensor that describes the curvature of M and T ab is the stress-energy-momentum tensor describing the distribution of matter. 1 Free particles travel along timelike geodesics and light rays travel along null geodesics. Thus gravity is a manifestation of the curvature of spacetime. Space acts on matter, telling it how to move. In turn, matter reacts back on space, telling it how to curve. 3.1 Solutions to the Einstein field equation The Einstein field equation might look simple when written with tensors. However, it constitutes a system of coupled, non-linear partial differential equations Minkowski spacetime The spacetime of special relativity Minkowski 3 spacetime is a solution to the vacuum field equation G ab = 0. It is the flat 4 spacetime (M, η ab ) given by the constant metric (η ab ) = (3.) 1 Some authors use a different sign in the definition of the Ricci tensor resulting in a minus sign in front of the right-hand side. Furthermore, we use geometrized units where the gravitational constant G, and the speed of light c, are set equal to one. Brief explanation of general relativity by the student portrayed in chapter one of [8]. 3 Hermann Minkowski ( ), Russian-born German mathematician. 4 R abc d = 0. Hansevi,

28 16 Chapter 3. General relativity In the coordinate system implied by η ab, all geodesics appear straight, that is, they take the form, x a (t) = y a + z a t Schwarzschild spacetime One of the solutions of the vacuum field equation was discovered by Karl Schwarzschild 5 in 1916, just a couple of months after Einstein published his field equation. The Schwarzschild solution is the unique solution that describes the curved spacetime exterior to a static spherically symmetric mass, such as a (non-rotating) star, planet, or black hole, and it remains one of the most important exact solutions. Schwarzschild spacetime is the curved spacetime (M, g ab ) given, in Schwarzschild coordinates, by the metric 1 M r ( ) 1 (g ab ) = 0 1 M r r 0. (3.3) r sin θ In these coordinates, the metric becomes singular at the surface r = M, which is called the event horizon. Events inside (or on) this surface cannot affect an outside observer; nothing can escape to the outside, not even light. 5 Karl Schwarzschild ( ), German physicist.

29 Chapter 4 Mass in general relativity The stress-energy-momentum tensor, T ab, is the tensor that describes the density and flux of energy and momentum in spacetime. It represents the energy due to matter and electromagnetic fields. Mass is the source of gravity, and energy is associated with mass, 1 therefore T ab is in the right-hand side of the Einstein field equation telling spacetime how to curve. 4.1 Gravitational energy/mass Imagine a system of two massive bodies at rest relative to each other. If they are far apart, then there will be a gravitational potential energy contribution that makes the total energy of the system greater than if they are close to each other. There is a difference in total energy, despite that integrating the energy densities, T ab, yields the same result in both scenarios. That energy difference is the energy attributed to the gravitational field. Since the gravitational field has energy, and therefore mass, it is a source of gravity, hence it is coupled to itself. Mathematically, this is possible because the field equation is non-linear Non-locality of mass The contribution of gravitational mass should be included in a description for total mass in a spacelike volume of spacetime. The stress-energy-momentum tensor T ab is given as a pointwise density and can be integrated over the volume. Can the gravitational mass be given as a point density that we can integrate? The answer is no. At any point in spacetime, one can always find a local coordinate system (Riemann-normal coordinates) in which all Christoffel symbol components vanish, which mean that there is no local gravitational field, hence no local gravitational mass [8, 14]. In such a coordinate system, an observer in free fall moves along a straight line, and does not feel any gravitational forces. 1 According to E = mc ; Einstein s famous equation [3]. Hansevi,

30 18 Chapter 4. Mass in general relativity 4. Total mass of an isolated system Despite the fact that the gravitational mass cannot be given as a pointwise density, there exist meaningful notions of the total mass of an isolated system Asymptotically flat spacetimes Say that we want to study the system of the two massive bodies (of section 4.1). Even though no physical system can be truly isolated from the rest of the universe, we can simplify our model by pretending that the system is isolated, ignoring the influence of distant matter. The spacetime of our simplified model will have vanishing curvature at large distances from the two bodies, and we say that it is asymptotically flat. A precise and useful, but rather technical, definition of an asymptotically flat spacetime is given in [14]. 4.. ADM and Bondi-Sachs mass In asymptotically flat spacetimes, the total mass can be determined by the asymptotic form of the metric. The ADM mass is the total mass measured at spacelike infinity, whereas the Bondi-Sachs 3 mass is measured at future null infinity [14]. 4.3 Quasi-local mass A meaningful definition of mass at a quasi-local level, that is, for the mass within a compact spacelike -surface, should have certain properties. For example, the quasi-local mass should be uniquely defined for all domains. Furthermore, it should be strictly positive (except in the flat case, where it should be equal to zero). Its limits at spacelike infinity and future null infinity should be the ADM mass and Bondi-Sachs mass, respectively. It should be monotone, that is, the mass for a domain should be greater or equal to the mass for a domain that is contained in the first. One can ask if it is possible to find a satisfying definition of mass at a quasilocal level. Several attempts have been made, for example, the Dougan-Mason mass, the Komar mass, the Penrose mass, and the Hawking mass. However, they fail to agree on the mass for a spacelike cross section of the event horizon of a Kerr 4 black hole [1]. There is still no generally accepted definition. To reach an appropriate definition for quasi-local mass would certainly be of great value. It could give a more detailed characterization of the gravitational field around massive bodies, and it should be helpful for controlling errors in numerical calculations [13]. R. Arnowitt, S. Deser and C. Misner. 3 H. Bondi and R. K. Sachs. 4 A solution to the Einstein field equation found in 1963 by Roy Kerr. It describes spacetime outside a rotating black hole.

31 Chapter 5 The Hawking mass One of many suggestions that have been made for a definition of quasi-local mass originates from a paper about gravitational radiation written by Stephen Hawking [5] in The Hawking mass has been shown to have various desirable properties, for example, the limits at spacelike infinity and future null infinity are the ADM mass and the Bondi-Sachs mass, respectively [13]. Its advantage is its simplicity, calculability and monotonicity for special families of -surfaces. In Minkowski spacetime, the Hawking mass vanish for -spheres. However, it can give negative results for general -surfaces, for example for non-convex -surfaces. 5.1 Definition Let l a and n a be, respectively, the outgoing and the ingoing null vectors orthogonal to a spacelike -sphere S and let m a and m a be tangent vectors to S. Then the Hawking mass is defined by m H (S) := Area(S) 16π ( ) ρρ ds. (5.1) π S 5. Interpretation Consider a one-parameter family of null geodesics, that is light rays, that intersects a circle in a -plane. Following the geodesics in the future direction, the optical scalars, θ and σ, introduced by Rainer Sachs [1] in 1961, can be defined as θ = Re ρ and σ = Im ρ, (5.) and interpreted as the expansion and rotation, respectively, of the circle []. Both l a and n a are orthogonal to the surface, thus both ρ and ρ are real [10]. Hence ρ and ρ measure the expansion of outgoing and ingoing null geodesics, respectively, and the Hawking mass can be viewed as a measure of the bending of outgoing and ingoing light rays orthogonal to the surface of a spacelike -sphere. Hansevi,

32 0 Chapter 5. The Hawking mass 5.3 Method for calculating spin coefficients In this section, a method for calculating the spin coefficients used in the expression for the Hawking mass is provided. In order to calculate the required spin coefficients (section.10), introduce a coordinate system such that for a given r, and t held fixed, 0 θ < π and 0 ϕ < π trace out the surface S. In the new coordinates, every vector with the first two components vanished lies in the tangent plane of S, as m a is required to do. The other two vectors, l a and n a, are parallel to the outgoing and ingoing null directions, respectively, thus it is convenient to set up a null tetrad (l a, n a, m a, m a ), by making the ansatz l a = ( A B C 0 ), n a = ( A B C 0 ), m a = ( 0 0 X iy ). (5.3) By imposing the conditions (.36) to (5.3), A, B, C, X and Y can be determined by solving the obtained equations in three steps. 1. Solve the system { ma m a = 0 (null vector) m a m a = 1 (normalization), (5.4) and pick one solution {X, Y }.. Solve the system { la l a = 0 (null vector) l a n a = 1 (normalization), (5.5) and pick the solution {A, B}, where A > 0 (and B > 0 if possible). The solution will possibly be dependent of C. 3. Determine C (and A and B if they depend on C) by solving l a m a = 0 (orthogonality). (5.6) The result of the procedure above is a null tetrad (l a, n a, m a, m a ) satisfying the conditions (.36). Calculate the spin coefficients with the use of the null tetrad by first calculating the required λ ijk given by (.35), λ 314 = ( b l a a l b ) m a m b λ 431 = ( b m a a m b ) m a l b λ 143 = ( b m a a m b ) l a m b λ 43 = ( b n a a n b ) m a m b λ 34 = ( b m a a m b ) m a n b λ 34 = ( b m a a m b ) n a m b then the spin coefficients follow easily as,, (5.7) ρ = γ 314 = (λ λ 431 λ 143 )/ ρ = γ 43 = (λ 43 + λ 34 λ 34 )/. (5.8)

33 5.4. Hawking mass for a -sphere Hawking mass for a -sphere In the following subsections, we demonstrate the method of section 5.3 by calculating the Hawking mass for a -sphere In Minkowski spacetime Using coordinates (t, r, θ, ϕ), where the spatial part (r, θ, ϕ) is written in spherical polar coordinates, the Minkowski metric (3.) can be written as, (g ab ) = r 0. (5.9) r sin θ By solving the equations (5.4), (5.5) and (5.6) in three steps, we obtain the null tetrad (l a, n a, m a, m a ), given by ( ) l a = , ( ) n a = , (5.10) ( ) m a 1 i = 0 0, r r sin θ and then it follows easily from (5.7) and (5.8) that The area of the -sphere is the familiar ρρ = 1 r. (5.11) Area(S r ) = π π 0 0 r sin θ dθ dϕ = 4πr, (5.1) and the integral of the spin coefficients S r ρρ ds = π π sin θ dθ dϕ = π. (5.13) Finally, we can confirm the well-known result that the Hawking mass vanish for all -spheres in Minkowski spacetime, since ( 4πr m H (S r ) = ) 16π π ( π) = ( r ) 1 1 = 0. (5.14) 5.4. In Schwarzschild spacetime We repeat the calculations of the previous section, but this time for a centered -sphere in Schwarzschild spacetime. The coordinates (t, r, θ, ϕ) of the

34 Chapter 5. The Hawking mass Schwarzschild metric (3.3) have the property that 0 θ < π and 0 ϕ < π trace out the surface S of a -sphere, thus we do not have to change to another coordinate system. We obtain a null tetrad (l a, n a, m a, m a ) given by l a = n a = m a = from which it follows that ( r r M ( r r M ( r ) r M 0 0, r ) r M 0 0, (5.15) r ) i, r sin θ ρρ = r M r 3. (5.16) The area of the -sphere is 4πr, and the integral of the spin coefficients S r ρρ ds = π π 0 0 r M r sin θ dθ dϕ = π(r M). (5.17) r Finally, we get the expected result for the Hawking mass of a centered -sphere in Schwarzschild spacetime, ( 4πr m H (S r ) = ) π(r M) = r ( M ) = M. (5.18) 16π π r r

35 Chapter 6 Results In this chapter, we calculate the Hawking mass, with the aid of Maple 1, for ellipsoidal -surfaces in both Minkowski and Schwarzschild spacetimes. The calculations for Minkowski are performed symbolically, and the results are presented as a theorem and two corollaries. In Schwarzschild spacetime, the Hawking mass are calculated numerically for approximately ellipsoidal -surfaces, and the results are therefore presented with diagrams. First, we prove a lemma that will be useful in section 6.1. Lemma 6.1. Assume that ξ > 1. Then arccosh ξ ξ 1 = 1 ξ 1 (ξ 1) + + O ( (ξ 1) 5/) Proof. Consider the equivalence, ξ = x + 1, x > 0 x = ξ 1, ξ > 1, (6.1) and the following Maclaurin expansions [11], x + 1 = 1 + x / x 4 /8 + O(x 6 ), (6.) ln(y + 1) = y y / + y 3 /3 y 4 /4 + y 5 /5 + O(y 6 ). (6.3) Perform a change of variable, arccosh ξ ξ 1 (def.) = ln(ξ + ξ 1) ξ 1 (6.1) = ln ( x x ) x, (6.4) and apply the Maclaurin expansions to the numerator, ln ( x x ) (6.) = ln ( 1 + x + x x4 8 + O(x6 ) ) (6.3) = x + x x4 8 + O(x6 ) 1 Mathematics software package from Waterloo Maple Inc. ( x + x x4 8 + O(x6 ) ) Hansevi,

36 4 Chapter 6. Results Divide through by x, and change back to ξ, ( x + x + + O(x4 ) ) 3 ( x + x + O(x4 ) ) ( x + O(x ) ) O(x 6 ) 5 = x x x O(x6 ). (6.5) arccosh ξ ξ 1 = 1 ξ 1 + 3(ξ 1) + O ( (ξ 1) 5/). (6.6) 6 40 Since ξ > 1, the ordo-term is actually Thus, let O ( (ξ 1) 5/) = O ( (ξ + 1) 5/ (ξ 1) 5/) = O ( (ξ 1) 5/). (6.7) f(ξ) = 1 ξ 1 + 3(ξ 1), (6.8) 6 40 and calculate the second degree Taylor expansion of f about ξ = 1. From f(1) = 1, f (1) = 1/3, and f (1) = 4/15 it follows that arccosh ξ ξ 1 = 1 ξ 1 (ξ 1) + + O ( (ξ 1) 5/). (6.9) Hawking mass in Minkowski spacetime If we let an ellipse rotate around its minor axis we get the surface of a rotationally symmetric ellipsoid called an oblate spheriod, see fig 6.1 (page 5). In Cartesian coordinates, the surface is given by the equation x + y A + z = 1. (6.10) B Let A and B depend on a variable, say r > 0, in the way given by A = ξ r and B = r. Then (6.10) is equivalent to x + y ξ + z = r, (6.11) which is the equation for an oblate spheroid, where r is the length of its minor axis. In this section, a closed-form expression for the Hawking mass within such an ellipsoid is given as a theorem with a proof. The result is also displayed as a graph in figure 6. (page 7). Furthermore, two corollaries regarding limits of the Hawking mass are proved.

37 6.1. Hawking mass in Minkowski spacetime 5 Figure 6.1: Oblate spheroid Closed-form expression of the Hawking mass Theorem 6.1. Let (M, η ab ) be Minkowski spacetime. For ξ > 1, let S r be the spacelike oblate -spheroids, in M, given by S r = { (x + y )/ξ + z = r : r > 0 }. Then the Hawking mass within S r is given by ( ) ( ) r m H Sr = 16 ξ + arccosh ξ ξ 5 ξ + arccosh ξ. ξ ξ 1 3 ξ 1 Proof. Use the method provided in section 5.3. Start by introducing a coordinate system such that for a given r, and t held fixed, θ and ϕ trace out the ellipsoidal -surface S r. This is accomplished by using the parametrized form of S r, x = ξ r sin θ cos ϕ y = ξ r sin θ sin ϕ z = r cos θ, 0 < r 0 θ < π 0 ϕ < π, (6.1) as new variables (t, r, θ, φ). Using the tensor transformation law (.19) yields the Minkowski metric (3.) in the new coordinates, (g ab ) = Ansatz (5.3) (cos θ + ξ sin θ) (ξ 1)r cos θ sin θ 0 0 (ξ 1)r cos θ sin θ r (ξ cos θ + sin θ) ξ r sin θ Solve equations (5.4), and pick a solution, say { X = 1 / ( r ξ cos θ + sin θ) Y = 1 / ( ξ r sin θ). (6.13). (6.14) Solve equations (5.5), and pick the solution given by A = 1/ and B = cos θ + ξ sin θ ξ r C (ξ 1) r cos θ sin θ C (cos θ + ξ sin θ). (6.15)

38 6 Chapter 6. Results Solve equation (5.6). The solution is given by (ξ 1) cos θ sin θ C = ξ r ξ cos θ + sin θ. (6.16) The result of this procedure is the null tetrad (l a, n a, m a, m a ), given by ( l a ξ cos = 1 θ + sin ) θ (ξ 1) cos θ sin θ ξ ξ r 0, ξ cos θ + sin θ ( n a ξ cos = 1 θ + sin θ (ξ ) 1) cos θ sin θ ξ ξ r 0, ξ cos θ + sin θ ( ) m a 1 = 0 0 r i, ξ cos θ + sin θ ξ r sin θ ( ) m a 1 = 0 0 r i. (6.17) ξ cos θ + sin θ ξ r sin θ Calculate the spin coefficients with the use of the null tetrad by first calculating the required λ ijk given by (5.7), λ 314 = ( b l a a l b )m a m b = 0, λ 431 = ( b m a a m b ) m a l b = λ 143 = ( b m a a m b )l a m b = λ 43 = ( b n a a n b ) m a m b = 0, λ 34 = ( b m a a m b )m a n b = ξ (1 + cos θ) + sin θ ξ r (ξ cos θ + sin θ) 3/, ξ (1 + cos θ) + sin θ ξ r (ξ cos θ + sin θ) 3/, ξ (1 + cos θ) + sin θ ξ r (ξ cos θ + sin θ) 3/, λ 34 = ( b m a a m b )n a m b ξ (1 + cos θ) + sin θ = ξ r (ξ cos θ + sin θ), 3/ (6.18) then the spin coefficients (5.8) follow easily as, ρ = γ 314 = ρ = γ 43 = ξ (1 + cos θ) + sin θ ξ r ( ξ cos θ + sin θ ) 3/, ξ (1 + cos θ) + sin θ ξ r ( ξ cos θ + sin θ ) 3/. (6.19) The surface area of S r is part of the expression for the Hawking mass. From the surface element ds = g θθ g ϕϕ (g θϕ ) dθ dϕ it follows that Area(S r ) = π π 0 0 ξ r sin θ ξ cos θ + sin θ dθ dϕ = π ξ r ( ξ ξ 1 + ln( ξ 1 + ξ) ) ξ 1

39 6.1. Hawking mass in Minkowski spacetime 7 ( = π ξ r ξ + arccosh ξ ), (6.0) ξ 1 Integration over S r yields π π 0 0 ρρ ds = π π 0 0 = π 6 = π 6 sin θ ( ξ (1 + cos θ) + sin θ ) ξ ( ξ cos θ + sin θ ) dθ dϕ 5/ ( ξ 3 + 7ξ) ξ arcsinh ( ξ 1) ( ξ ξ ξ ξ 1 arccosh ξ ξ 1 ). (6.1) Finally, from (6.0) and (6.1), the Hawking mass in S r follows as m H (S r ) = Area(Sr ) = r 16 ξ 16 π ( π ξ + arccosh ξ ξ 1 π π 0 0 ) ρρ ds ( ξ 5 ξ + arccosh ξ 3 ξ 1 ). (6.) As can be seen in figure 6., the Hawking mass becomes negative in Minkowski, even for a convex -surface Figure 6.: m H (S 1 ) plotted against parameter ξ Limit when approaching a metric sphere Corollary 6.1. Let (M, η ab ) be Minkowski spacetime. Given an r > 0, the Hawking mass vanishes in the limit when ξ 1 +, that is, when S r tends to a -sphere. Proof. It follows from lemma 6.1 that arccosh ξ ξ 1 = 1 ξ (ξ 1) 15 + O ( (ξ 1) 5/) = 1 + O(ξ 1), (6.3)

40 8 Chapter 6. Results and furthermore that ξ 5 3 ξ + arccosh ξ ξ 1 = ξ3 3 5 ξ O(ξ 1) 3 = ξ(ξ + 1) 3( ) ξ 1 + O(ξ 1) 3 = O(ξ 1). (6.4) Thus, it follows from (6.3), (6.4) and theorem 6.1 that m H (S r ) = ( ) r 16 ξ + arccosh ξ ξ 5 ξ + arccosh ξ ξ ξ 1 3 ξ 1 = r 1 + ξ + O(ξ 1) O(ξ 1) 0 as ξ 1 +. (6.5) ξ Limit along a foliation Corollary 6.. Let (M, η ab ) be Minkowski spacetime. For ω > 0, let {L r } r>0 be foliations of a spacelike 3-surface Ω M. Suppose that the leaves are given by L r = { (x + y )/(1 + ω/r) + z = r : r > 0 }. Then the Hawking mass vanishes in the limit along all foliations {L r } r>0. Proof. Observe that L r = S r ξ=1+ω/r, (6.6) and let ξ = 1 + ω/r in lemma 6.1. Then ξ 1 + as r, and it follows that 1 + ω/r + arccosh (1 + ω/r) (1 + ω/r) 1 = 1 + ω r + 1 ω 3r + ω 15r + O ( 1 r 5/ ) = + ω 3r + ( ω 1 ) 15r + O r ( 5/ 1 ) = + O, (6.7) r and furthermore that (1 + ω/r) 5 arccosh (1 + ω/r) (1 + ω/r) + 3 (1 + ω/r) 1 ( = 3ω 15r + ω3 1 ) 3r 3 + O r 5/ ( 1 ) = O r. (6.8) From (6.7) and theorem 6.1, it follows that r ( 1 ) ( 1 ) m H (L r ) = ω + O O r r r ( 1 ) = O 0 as r. (6.9) r

41 6.. Hawking mass in Schwarzschild spacetime 9 6. Hawking mass in Schwarzschild spacetime In this section we will study the Hawking mass in the curved Schwarzschild spacetime. A consequence of the spacetime being curved is that the required calculations tend to be more complicated. Considering this, we will calculate the Hawking mass for an approximately ellipsoidal -surface, exterior to the event horizon (r > M), that has a simple expression in spherical polar coordinates. For small ε > 0, let the curve given by r = 1 + ε sin θ, (6.30) rotate around the axis θ = 0. See figure 6.3. The surface of revolution, S, is approximately ellipsoidal. For ε = 0, it is a perfect sphere. θ = 0 θ r Figure 6.3: A curve given by r = 1 + ε sin θ. Let us introduce a coordinate system such that for a given r, and t held fixed, θ and ϕ trace out the spacelike -surface S r. This is accomplished by the following change of variables, r r 1 + ε sin θ, r > M. (6.31) In the new coordinates, the Schwarzschild metric (3.3) is given by β rα rα3 β εr α cos θ sin θ β 0 (g ab ) = 0 εr α cos θ sin θ β r (α 3 β+ε r cos θ sin θ) αβ 0, r α sin θ (6.3) where α = 1 + ε sin θ > 1 and β = r 1 + ε sin θ M > 0. (6.33) By making the ansatz (5.3), solving equations (5.4), (5.5) and (5.6), we obtain a null tetrad (l a, n a, m a, m a ) given by, l a = ( rα β α3 β + ε r cos θ sin θ rα 3 ) C 0,

42 30 Chapter 6. Results n a = m a = ( rα β ( 0 0 α3 β + ε r cos θ sin θ rα 3 αβ r α 3 β + ε r cos θ sin θ ) C 0, (6.34) ) i, rα sin θ where ε cos θ sin θ C = rα α3 β + ε r cos θ sin θ. (6.35) We calculate the spin coefficients ρ and ρ. Unfortunately, their expressions in Schwarzschild spacetime are very long, so we omit writing them out. The surface element of S r is given by ds = and the surface area by g θθ g ϕϕ (g θϕ ) dθ dϕ = r α sin θ ε r cos θ sin θ + 1 dθ dϕ, (6.36) αβ Area( S r ) = π π 0 0 Since the expression of the Hawking mass for S r, m H ( S r ) = Area( S ( r ) π π ds. (6.37) π π 0 0 ) ρρ ds, (6.38) is rather complicated, we will rely on methods like Taylor expansion and numerical integration for further investigations Limit along a foliation Analogously to corollary 6., let { L r } r>m be foliations of a spacelike 3-surface Ω M (now in Schwarzschild spacetime). Suppose that the leaves, Lr, are given by substituting ε = ω, r > M, (6.39) r into S r. It is easy to see that L r becomes more spherical the larger r gets. In this section, we will show that the limit of the Hawking mass of L r is M along all foliations { L r } r>m, that is, lim m H( L r ) = M. (6.40) r Since both ds and ρρ ds are independent of ϕ, we make the substitution (6.39), and let I 1 (ε) = ε Area( L ω/ε ) =: π π 0 Φ(ε, θ) dθ (6.41)

43 6.. Hawking mass in Schwarzschild spacetime 31 and I (ε) = π π 0 ρρ dθ =: π π 0 Ψ(ε, θ) dθ. (6.4) Under the reasonable assumption that Φ(ε, θ), Φ ε(ε, θ), Ψ(ε, θ), Ψ ε(ε, θ) and Ψ εε(ε, θ) are continuous in some neighbourhood ε < 1, it follows by Maclaurin expansion that and I 1 (ε) = π π 0 I (ε) = π Φ(0, θ) dθ + O(ε) = 4πω + O(ε), ε < 1, (6.43) π 0 Ψ(0, θ) dθ + π π 0 Ψ ε(0, θ) dθ ε + O(ε ) = π + 4πM ω ε + O(ε ), ε < 1. (6.44) Thus Area( L r ) = I 1 (ε)/ε = 4πω 1 ) ε + O( ε and π π 0 ρρ ds = π + 4πM r = 4πr + O(r), r > ω, (6.45) ( 1 ) + O r, r > ω. (6.46) By using the results (6.45) and (6.46), it follows that ( m H ( L 4πr + O(r) r ) = [ π + 4πM ( 1 )] ) + O 16π π r r = r ( 1 ) ( M ( 1 ) ) 1 + O + O r r r ( 1 ) ( ( 1 ) ) = 1 + O M + O M as r. r r (6.47) 6.. Numerical evaluations We evaluate the Hawking mass numerically with an adaptive Gaussian quadrature method. In numerical quadrature, an integral I(f) = b a f(x)dx, (6.48) is approximated by an n-point quadrature rule that has the form n Q n (f) = w i f(x i ), (6.49) i=1 where a x 1 < x < < x n b. The points x i are called nodes, and the multipliers w i are called weights. Carl Friedrich Gauss ( ), German mathematician known as princeps mathematicorum ( prince of mathematicians ).

44 3 Chapter 6. Results In Gaussian quadrature, both the nodes and the weights are optimally chosen, hence Gaussian quadrature has the highest possible accuracy for the number of nodes used. Furthermore, it is stable and Q n (f) I(f) as n [6]. In adaptive quadrature, the interval of integration is selectively refined to reflect the behavior of the particular integrand. We set M = 1, and let Maple calculate the Hawking mass numerically for some values of ε. The results are presented as graphs. See figure 6.4, 6.5 and e = e = 0.10 e = 0.15 K0.5 e = 0.0 K e = 0.5 Figure 6.4: m H ( S r ) plotted against 4 r 60 for some values of ε e = 0. e = 0. e = 0.4 e = 0.6 e = e = Figure 6.5: m H ( S r ) plotted against.3 r 16 for some values of ε Figure 6.6: m H ( S r ) plotted against.3 r 16 for ε = ω/r.

45 Chapter 7 Discussion 7.1 Conclusions In this thesis, we have derived a closed-form expression for the Hawking mass of a spacelike oblate -spheroid S r, that is, a rotationally symmetric -ellipsoid, in Minkowski spacetime. If S r is given by, S r = { (x + y )/ξ + z = r : r > 0 }. Then the Hawking mass within S r is given by ( ) ( ) r m H Sr = 16 ξ + arccosh ξ ξ 5 ξ + arccosh ξ. ξ ξ 1 3 ξ 1 From this result, we can see that the Hawking mass can be negative even for convex -surfaces in Minkowski spacetime. However, the limits along particular foliations, were shown to vanish. Furthermore, we studied the Hawking mass in Schwarzschild spacetime. Numerical calculations for approximately ellipsoidal -surfaces were done, and the results were presented with diagrams that show that the Hawking mass can be negative in Schwarzschild. It this case, the limits along particular foliations, were shown to be M, that is, equal to the Hawking mass for a centered -sphere. 7. Future work The calculations performed in this thesis can serve as a basis for similar studies of the Hawking mass in other spacetimes, for example in Reissner-Nordström which is the generalization of Schwarzschild that includes electric charge. Hansevi,

46 34 Chapter 7. Discussion

47 Bibliography [1] Bergqvist, G., Quasilocal mass for event horisons, Class. Quantum Grav. 9, , (199) [] Chandrasekhar, S., The Mathematical Theory of Black Holes, Oxford University Press, New York, 1983 [3] Einstein, A., Zur Elektrodynamik bewegter Körper, Annalen der Physik 18, , (1905) [4] Einstein, A., Die Grundlage der allgemeinen Relativitätstheorie, Annalen der Physik 49, 769-8, (1916) [5] Hawking, S. W., Gravitational Radiation in an Expanding Universe, J. Math. Phys. 9, , (1968) [6] Heath, M. T., Scientific Computing: An Introductory Survey, McGraw-Hill, 00 [7] Lawson, H. B., Jr., Foliations, Bull. Amer. Math. Soc. 80, , (1974) [8] Misner, C. W, Thorne, K. S., Wheeler, J. A, Gravitation, W. H. Freeman and Company, 1973 [9] Newman, E. T. and Penrose, R., An approach to Gravitational Radiation by a Method of Spin Coefficients, J. Math. Phys. 3, , (196) [10] Penrose, R. and Rindler, W., Spinors and space-time volume 1, Cambridge University Press, 1984 [11] Råde, L., Westergren, B., Mathematics Handbook for Science and Engineering, Studentlitteratur, Lund, 1995 [1] Sachs, R. K., Gravitational Waves in General Relativity. VI. The Outgoing Radiation Condition, Proc. Roy. Soc. (London) A 64, , (1961) [13] Szabados, L. B., Quasi-Local Energy-Momentum and Angular Momentum in GR: A Review Article, Living Rev. Relativity 7, (004), 4. [Online Article] (cited on ): [14] Wald, R. M., General Relativity, The University of Chicago Press, Chicago and London, 1984 Hansevi,

48 36 Bibliography

49 Appendix A Maple Worksheets A.1 Null tetrad # Determine a null-tetrad # Let m1=m and m=conjugate(m), # m1*m1 = m*m = 0 (null vectors) and m1*m = -1 (normalization) m1:=create([1],array(1..4,[0,0,x, I*Y])): m:=create([1],array(1..4,[0,0,x,-i*y])): sy:=solve((get_compts(prod(m1,lower(g,m1,1),[1,1])))=0,y)[1]: m1:=create([1],array(1..4,[0,0,x, I*sY])): m:=create([1],array(1..4,[0,0,x,-i*sy])): sx:=solve((get_compts(prod(m1,lower(g,m,1),[1,1])))=-1,x)[1]: m1:=create([1],array(1..4,[0,0,sx, I*eval(sY,X=sX)])): m:=create([1],array(1..4,[0,0,sx,-i*eval(sy,x=sx)])): # Let l and m, such that l*l = n*n = 0 (null vectors) # and l*n = 1 (normalization) l:=create([1],array(1..4,[a, B, C,0])): n:=create([1],array(1..4,[a,-b,-c,0])): sb:=solve((get_compts(prod(l,lower(g,l,1),[1,1])))=0,b)[1]: l:=create([1],array(1..4,[a, sb, C,0])): n:=create([1],array(1..4,[a,-sb,-c,0])): sa:=solve((get_compts(prod(l,lower(g,n,1),[1,1])))=1,a)[1]: sb:=eval(sb,a=sa): l:=create([1],array(1..4,[sa, sb, C,0])): n:=create([1],array(1..4,[sa,-sb,-c,0])): # l*m1 = l*m = n*m1 = n*m = 0 (orthogonality) sc:=solve((get_compts(prod(l,lower(g,m1,1),[1,1])))=0,c)[1]: sa:=simplify(eval(sa,c=sc)): sb:=simplify(eval(sb,c=sc)): l:=create([1],array(1..4,[sa, sb, sc,0])): n:=create([1],array(1..4,[sa,-sb,-sc,0])): # Control (result should be: 0,0,0,0,0,0,0,0,1,-1) simplify(get_compts(prod(l, lower(g,l, 1),[1,1]))), simplify(get_compts(prod(n, lower(g,n, 1),[1,1]))), simplify(get_compts(prod(m1,lower(g,m1,1),[1,1]))), simplify(get_compts(prod(m,lower(g,m,1),[1,1]))), simplify(get_compts(prod(l, lower(g,m1,1),[1,1]))), simplify(get_compts(prod(n, lower(g,m1,1),[1,1]))), simplify(get_compts(prod(l, lower(g,m,1),[1,1]))), simplify(get_compts(prod(n, lower(g,m,1),[1,1]))), simplify(get_compts(prod(l, lower(g,n, 1),[1,1]))), simplify(get_compts(prod(m1,lower(g,m,1),[1,1]))); Hansevi,

50 38 Appendix A. Maple Worksheets A. Minkowski restart:with(tensor): assume(xi::real,t::real,r::real,theta::real,phi::real): additionally(xi>=1,r>0,theta>=0,theta<=pi,phi>=0,phi<=*pi): # Metric of Minkowski spacetime ("cartesian" coordinates) coords:=[t,r,theta,phi]: h_c:=array(sparse,1..4,1..4): h_c[1,1]:=1: h_c[,]:=-1: h_c[3,3]:=-1: h_c[4,4]:=-1: # Change of variables (theta and phi parametrize ellipsoids) x[1]:=t: x[]:=xi*r*sin(theta)*cos(phi): x[3]:=xi*r*sin(theta)*sin(phi): x[4]:=r*cos(theta): g_c:=array(sparse,1..4,1..4): for i_ from 1 to 4 do for j_ from 1 to 4 do for i from 1 to 4 do for j from 1 to 4 do g_c[i_,j_]:=g_c[i_,j_] +h_c[i,j]*diff(x[i],coords[i_])*diff(x[j],coords[j_]): end do; end do; end do; end do; # Create (covariant) metric tensor g g:=create([-1,-1],eval(g_c)): # Null tetrad l:=create([1],array(1..4,[sqrt()/, sqrt()//xi*sqrt(xi***cos(theta)**+sin(theta)**), -sqrt()//xi/r*(xi**-1)*cos(theta)*sin(theta) /sqrt(xi***cos(theta)**+sin(theta)**), 0])): n:=create([1],array(1..4,[sqrt()/, -sqrt()//xi*sqrt(xi***cos(theta)**+sin(theta)**), sqrt()//xi/r*(xi**-1)*cos(theta)*sin(theta) /sqrt(xi***cos(theta)**+sin(theta)**), 0])): m1:=create([1],array(1..4,[0,0, sqrt()//r/sqrt(xi***cos(theta)**+sin(theta)**), sqrt()/*i/xi/r/sin(theta)])): m:=create([1],array(1..4,[0,0, sqrt()//r/sqrt(xi***cos(theta)**+sin(theta)**), -sqrt()/*i/xi/r/sin(theta)])): # Calculate spin coefficients e[1]:=get_compts(l): e[]:=get_compts(n): e[3]:=get_compts(m1): e[4]:=get_compts(m): f[1]:=get_compts(lower(g,l, 1)): f[]:=get_compts(lower(g,n, 1)): f[3]:=get_compts(lower(g,m1,1)): f[4]:=get_compts(lower(g,m,1)): L314:=0: L431:=0: L143:=0: L43:=0: L34:=0: L34:=0: for i from 1 to 4 do for j from 1 to 4 do L314:=L314 +(diff(f[1][i],coords[j])-diff(f[1][j],coords[i]))*e[3][i]*e[4][j]: L431:=L431 +(diff(f[3][i],coords[j])-diff(f[3][j],coords[i]))*e[4][i]*e[1][j]: L143:=L143 +(diff(f[4][i],coords[j])-diff(f[4][j],coords[i]))*e[1][i]*e[3][j]:

51 A.3. Schwarszchild 39 L43:=L43 +(diff(f[][i],coords[j])-diff(f[][j],coords[i]))*e[4][i]*e[3][j]: L34:=L34 +(diff(f[4][i],coords[j])-diff(f[4][j],coords[i]))*e[3][i]*e[][j]: L34:=L34 +(diff(f[3][i],coords[j])-diff(f[3][j],coords[i]))*e[][i]*e[4][j]: end do; end do; # rho = gamma[314], rho = gamma[43] rho_ :=simplify(1/*(l314+l431-l143)); rho_:=simplify(1/*(l43+l34-l34)); # Calculate surface area of ellipsoid area_element:=simplify(sqrt(g_c[3,3]*g_c[4,4]-g_c[3,4]**)); area:=simplify(int(int(area_element,theta=0..pi),phi=0..*pi)); # Calculate integral of spin coefficients spin_integrand:=(rho_*rho_*area_element); spin_integral:=simplify( int(eval(int(spin_integrand, theta ),theta=pi)- eval(int(spin_integrand, theta ),theta=0),phi=0..*pi)); # Calculate Hawking mass mh_:=simplify(sqrt(area/16/pi)*(1+1//pi*spin_integral)); # Simplification by hand mh:=-sqrt()/(16*sqrt(xi))*r*sqrt(xi+arccosh(xi)/sqrt(xi**-1)) *(xi*(*xi**-5)/3+arccosh(xi)/sqrt(xi**-1)); # Limit along foliation limit(eval(mh,xi=1+omega/r),r=infinity); plot(eval(mh,r=1),xi=1..4,-3..1, thickness=,axes=frame,gridlines=true,font=[times,1,1]); A.3 Schwarszchild restart:with(tensor):with(plots): assume(t::real,r::real,theta::real,phi::real,epsilon::real,omega::real,k::real,p::real): additionally(m>=0,r>*m,epsilon>=0,theta>=0,theta<=pi/,alpha>=1,beta>0,omega>0): # Metric of Schwarzschild spacetime coords:=[t,r,theta,phi]: h_c:=array(sparse,1..4,1..4): h_c[1,1]:=1-*m/r: h_c[,]:=-1/(1-*m/r): h_c[3,3]:=-r**: h_c[4,4]:=-r***sin(theta)**: # Change of variables (theta and phi parametrize ellipsoids) F:=sqrt(1+epsilon*sin(theta)**): x[1]:=t: x[]:=f*r: x[3]:=theta: x[4]:=phi: gt_c:=array(sparse,1..4,1..4): for i_ from 1 to 4 do for j_ from 1 to 4 do for i from 1 to 4 do for j from 1 to 4 do gt_c[i_,j_]:=gt_c[i_,j_] +subs({t=x[1],r=x[],theta=x[3],phi=x[4]},h_c[i,j]) *diff(x[i],coords[i_])*diff(x[j],coords[j_]): end do; end do; end do; end do; # Create (covariant) metric tensor g gt:=create([-1,-1],eval(gt_c)): AB:={alpha=sqrt(1+epsilon*sin(theta)**), beta=r*sqrt(1+epsilon*sin(theta)**)-*m}: # Simplification by hand g_c:=array(sparse,1..4,1..4): g_c[1,1]:=beta/alpha/r:

52 40 Appendix A. Maple Worksheets g_c[,]:=-alpha**3*r/beta: g_c[3,3]:=-r***(alpha**3*beta+epsilon***r*sin(theta)***cos(theta)**) /(alpha*beta): g_c[4,4]:=-alpha***r***sin(theta)**: g_c[,3]:=-epsilon*alpha*r***sin(theta)*cos(theta)/beta: g_c[3,]:=-epsilon*alpha*r***sin(theta)*cos(theta)/beta: g:=create([-1,-1],eval(g_c)): # Null tetrad l:=create([1],array(1..4,[ sqrt()/*sqrt(r*alpha/beta), sqrt()/*sqrt(alpha**3*beta+r*epsilon***cos(theta)***sin(theta)**) /(sqrt(r)*alpha**3), -sqrt()/*epsilon*sin(theta)*cos(theta) /(alpha*sqrt(r)*sqrt(alpha**3*beta +r*epsilon***cos(theta)***sin(theta)**)), 0])): n:=create([1],array(1..4,[ sqrt()/*sqrt(r*alpha/beta), -sqrt()/*sqrt(alpha**3*beta+r*epsilon***cos(theta)***sin(theta)**) /(sqrt(r)*alpha**3), sqrt()/*epsilon*sin(theta)*cos(theta) /(alpha*sqrt(r)*sqrt(alpha**3*beta +r*epsilon***cos(theta)***sin(theta)**)), 0])): m1:=create([1],array(1..4,[ 0,0, sqrt()/*sqrt(beta)*sqrt(alpha) /(r*sqrt(alpha**3*beta+r*epsilon***cos(theta)***sin(theta)**)), -sqrt()/*i/(alpha*r*sin(theta))])): m:=create([1],array(1..4,[ 0,0, sqrt()/*sqrt(beta)*sqrt(alpha) /(r*sqrt(alpha**3*beta+r*epsilon***cos(theta)***sin(theta)**)), sqrt()/*i/(alpha*r*sin(theta))])): # Calculate spin coefficients e[1]:=simplify(eval(get_compts(l), AB)): e[]:=simplify(eval(get_compts(n), AB)): e[3]:=simplify(eval(get_compts(m1),ab)): e[4]:=simplify(eval(get_compts(m),ab)): f[1]:=simplify(eval(get_compts(lower(g,l, 1)),AB)): f[]:=simplify(eval(get_compts(lower(g,n, 1)),AB)): f[3]:=simplify(eval(get_compts(lower(g,m1,1)),ab)): f[4]:=simplify(eval(get_compts(lower(g,m,1)),ab)): L314:=0: L431:=0: L143:=0: L43:=0: L34:=0: L34:=0: for i from 1 to 4 do for j from 1 to 4 do L314:=L314 +(diff(f[1][i],coords[j])-diff(f[1][j],coords[i]))*e[3][i]*e[4][j]: L431:=L431 +(diff(f[3][i],coords[j])-diff(f[3][j],coords[i]))*e[4][i]*e[1][j]: L143:=L143 +(diff(f[4][i],coords[j])-diff(f[4][j],coords[i]))*e[1][i]*e[3][j]: L43:=L43 +(diff(f[][i],coords[j])-diff(f[][j],coords[i]))*e[4][i]*e[3][j]: L34:=L34 +(diff(f[4][i],coords[j])-diff(f[4][j],coords[i]))*e[3][i]*e[][j]: L34:=L34 +(diff(f[3][i],coords[j])-diff(f[3][j],coords[i]))*e[][i]*e[4][j]: end do; end do; # rho = gamma[314], rho = gamma[43] rho_ :=simplify(1/*(l314+l431-l143));

53 A.3. Schwarszchild 41 rho_:=simplify(1/*(l43+l34-l34)); # Hawking energy mh:=sqrt(area/16/pi)*(1+1//pi*spin_integral): # Surface area of ellipsoid area_element:=simplify(eval(sqrt(g_c[3,3]*g_c[4,4]-g_c[3,4]**),ab)): area:=simplify(*pi*int(area_element,theta=0..pi, method=_gquad)): # Integral of spin coefficients spin_integral:=*pi*int(simplify(rho_*rho_)*area_element,theta=0..pi, method=_gquad): # Limit along foliation Phi:=simplify(epsilon***eval(area_element,r=omega/epsilon)): *Pi*int(simplify(eval(Phi,epsilon=0)),theta=0..Pi)/epsilon**: Psi_:=simplify(eval(simplify(rho_*rho_*area_element),r=omega/epsilon)): *Pi*int(simplify(eval(Psi_,epsilon=0)),theta=0..Pi) +*Pi*int(simplify(eval(diff(Psi_,epsilon),epsilon=0)),theta=0..Pi)*epsilon: eval(series(eval(mh,r=1/epsilon),epsilon=0,),epsilon=1/r): # Plot eps1:=[0,1/5,/5,3/5,4/5,1]: start1:=.3: base1:=1.1: for i from 1 to 6 do R:=start1*base1**0: points1[i]:=[r,evalf(eval(mh,{m=1,r=r,epsilon=eps1[i]}))]: for a from 1 to 0 do R:=start1*base1**a: points1[i]:=points1[i],[r,evalf(eval(mh,{m=1,r=r,epsilon=eps1[i]}))]: end do: end do: plots[display]({ setoptions(symbol=diamond,symbolsize=1,color=black), plot([points1[1]]),pointplot([points1[1]]), plot([points1[]]),pointplot([points1[]]), plot([points1[3]]),pointplot([points1[3]]), plot([points1[4]]),pointplot([points1[4]]), plot([points1[5]]),pointplot([points1[5]]), plot([points1[6]]),pointplot([points1[6]]) }, textplot({[16,1.0,epsilon_=0.0],[16,0.95,epsilon_=0.], [16,0.84,epsilon_=0.4],[16,0.69,epsilon_=0.6], [16,0.5,epsilon_=0.8],[16,0.9,epsilon_=1.0]}, align={right},font=[times,roman,1]), axes=boxed, axis[1]=[mode=log,gridlines=[8,thickness=1,subticks=false,color=grey]], axis[]=[gridlines=[color=grey]], view=[..16, ] ); eps:=[0,0.05,0.1,0.15,0.,0.5]: start:=4: base:=1.4: for i from 1 to 6 do R:=start*base**0: points[i]:=[r,evalf(eval(mh,{m=1,r=r,epsilon=eps[i]}))]: for a from 1 to 15 do R:=start*base**a: points[i]:=points[i],[r,evalf(eval(mh,{m=1,r=r,epsilon=eps[i]}))]: end do: end do: plots[display]({ setoptions(symbol=diamond,symbolsize=1,color=black), plot([points[1]]),pointplot([points[1]]), plot([points[]]),pointplot([points[]]),

54 4 Appendix A. Maple Worksheets plot([points[3]]),pointplot([points[3]]), plot([points[4]]),pointplot([points[4]]), plot([points[5]]),pointplot([points[5]]), plot([points[6]]),pointplot([points[6]]) }, textplot({[700,1.0,epsilon_=0.00],[700,0.9,epsilon_=0.05], [700,0.6,epsilon_=0.10],[700,0.17,epsilon_=0.15], [700,-0.45,epsilon_=0.0],[700,-1.18,epsilon_=0.5]}, align={right},font=[times,roman,1]), axes=boxed, axis[1]=[mode=log,gridlines=[1,thickness=1,subticks=false,color=grey]], axis[]=[gridlines=[color=grey]], view=[4..700, ] ); om:=[0.0,0.,0.4,0.6,0.8,1]: start3:=.3: base3:=1.1: for i from 1 to 6 do R:=start3*base3**0: points3[i]:=[r,evalf(eval(mh,{m=1,r=r,epsilon=om[i]/r}))]: for a from 1 to 0 do R:=start3*base3**a: points3[i]:=points3[i],[r,evalf(eval(mh,{m=1,r=r,epsilon=om[i]/r}))]: end do: end do: plots[display]({ setoptions(symbol=diamond,symbolsize=1,color=black), plot([points3[1]]),pointplot([points3[1]]), plot([points3[]]),pointplot([points3[]]), plot([points3[3]]),pointplot([points3[3]]), plot([points3[4]]),pointplot([points3[4]]), plot([points3[5]]),pointplot([points3[5]]), plot([points3[6]]),pointplot([points3[6]]) }, axes=boxed, axis[1]=[mode=log,gridlines=[8,thickness=1,subticks=false,color=grey]], axis[]=[gridlines=[color=grey]], view=[...16, ] );

55 Copyright The publishers will keep this document online on the Internet - or its possible replacement - for a period of 5 years from the date of publication barring exceptional circumstances. The online availability of the document implies a permanent permission for anyone to read, to download, to print out single copies for your own use and to use it unchanged for any non-commercial research and educational purpose. Subsequent transfers of copyright cannot revoke this permission. All other uses of the document are conditional on the consent of the copyright owner. The publisher has taken technical and administrative measures to assure authenticity, security and accessibility. According to intellectual property law the author has the right to be mentioned when his/her work is accessed as described above and to be protected against infringement. For additional information about the Linköping University Electronic Press and its procedures for publication and for assurance of document integrity, please refer to its WWW home page: Upphovsrätt Detta dokument hålls tillgängligt på Internet - eller dess framtida ersättare - under 5 år från publiceringsdatum under förutsättning att inga extraordinära omständigheter uppstår. Tillgång till dokumentet innebär tillstånd för var och en att läsa, ladda ner, skriva ut enstaka kopior för enskilt bruk och att anvnda det oförändrat för ickekommersiell forskning och för undervisning. Överföring av upphovsrätten vid en senare tidpunkt kan inte upphäva detta tillstånd. All annan användning av dokumentet kräver upphovsmannens medgivande. För att garantera äktheten, säkerheten och tillgängligheten finns det lösningar av teknisk och administrativ art. Upphovsmannens ideella rätt innefattar rätt att bli nämnd som upphovsman i den omfattning som god sed kräver vid användning av dokumentet på ovan beskrivna sätt samt skydd mot att dokumentet ändras eller presenteras i sådan form eller i sådant sammanhang som är kränkande för upphovsmannens litterära eller konstnärliga anseende eller egenart. För ytterligare information om Linköping University Electronic Press se förlagets hemsida c 008, Daniel Hansevi Hansevi,

An introduction to General Relativity and the positive mass theorem

An introduction to General Relativity and the positive mass theorem An introduction to General Relativity and the positive mass theorem National Center for Theoretical Sciences, Mathematics Division March 2 nd, 2007 Wen-ling Huang Department of Mathematics University of

More information

Gravitation: Tensor Calculus

Gravitation: Tensor Calculus An Introduction to General Relativity Center for Relativistic Astrophysics School of Physics Georgia Institute of Technology Notes based on textbook: Spacetime and Geometry by S.M. Carroll Spring 2013

More information

Initial-Value Problems in General Relativity

Initial-Value Problems in General Relativity Initial-Value Problems in General Relativity Michael Horbatsch March 30, 2006 1 Introduction In this paper the initial-value formulation of general relativity is reviewed. In section (2) domains of dependence,

More information

has a lot of good notes on GR and links to other pages. General Relativity Philosophy of general relativity.

has a lot of good notes on GR and links to other pages. General Relativity Philosophy of general relativity. http://preposterousuniverse.com/grnotes/ has a lot of good notes on GR and links to other pages. General Relativity Philosophy of general relativity. As with any major theory in physics, GR has been framed

More information

carroll/notes/ has a lot of good notes on GR and links to other pages. General Relativity Philosophy of general

carroll/notes/ has a lot of good notes on GR and links to other pages. General Relativity Philosophy of general http://pancake.uchicago.edu/ carroll/notes/ has a lot of good notes on GR and links to other pages. General Relativity Philosophy of general relativity. As with any major theory in physics, GR has been

More information

Exact Solutions of the Einstein Equations

Exact Solutions of the Einstein Equations Notes from phz 6607, Special and General Relativity University of Florida, Fall 2004, Detweiler Exact Solutions of the Einstein Equations These notes are not a substitute in any manner for class lectures.

More information

Curved spacetime and general covariance

Curved spacetime and general covariance Chapter 7 Curved spacetime and general covariance In this chapter we generalize the discussion of preceding chapters to extend covariance to more general curved spacetimes. 219 220 CHAPTER 7. CURVED SPACETIME

More information

Geometric inequalities for black holes

Geometric inequalities for black holes Geometric inequalities for black holes Sergio Dain FaMAF-Universidad Nacional de Córdoba, CONICET, Argentina. 3 August, 2012 Einstein equations (vacuum) The spacetime is a four dimensional manifold M with

More information

Übungen zu RT2 SS (4) Show that (any) contraction of a (p, q) - tensor results in a (p 1, q 1) - tensor.

Übungen zu RT2 SS (4) Show that (any) contraction of a (p, q) - tensor results in a (p 1, q 1) - tensor. Übungen zu RT2 SS 2010 (1) Show that the tensor field g µν (x) = η µν is invariant under Poincaré transformations, i.e. x µ x µ = L µ νx ν + c µ, where L µ ν is a constant matrix subject to L µ ρl ν ση

More information

UNIVERSITY OF DUBLIN

UNIVERSITY OF DUBLIN UNIVERSITY OF DUBLIN TRINITY COLLEGE JS & SS Mathematics SS Theoretical Physics SS TSM Mathematics Faculty of Engineering, Mathematics and Science school of mathematics Trinity Term 2015 Module MA3429

More information

Quasi-local Mass in General Relativity

Quasi-local Mass in General Relativity Quasi-local Mass in General Relativity Shing-Tung Yau Harvard University For the 60th birthday of Gary Horowtiz U. C. Santa Barbara, May. 1, 2015 This talk is based on joint work with Po-Ning Chen and

More information

Null Cones to Infinity, Curvature Flux, and Bondi Mass

Null Cones to Infinity, Curvature Flux, and Bondi Mass Null Cones to Infinity, Curvature Flux, and Bondi Mass Arick Shao (joint work with Spyros Alexakis) University of Toronto May 22, 2013 Arick Shao (University of Toronto) Null Cones to Infinity May 22,

More information

Syllabus. May 3, Special relativity 1. 2 Differential geometry 3

Syllabus. May 3, Special relativity 1. 2 Differential geometry 3 Syllabus May 3, 2017 Contents 1 Special relativity 1 2 Differential geometry 3 3 General Relativity 13 3.1 Physical Principles.......................................... 13 3.2 Einstein s Equation..........................................

More information

RELG - General Relativity

RELG - General Relativity Coordinating unit: Teaching unit: Academic year: Degree: ECTS credits: 2017 230 - ETSETB - Barcelona School of Telecommunications Engineering 749 - MAT - Department of Mathematics 748 - FIS - Department

More information

Non-existence of time-periodic dynamics in general relativity

Non-existence of time-periodic dynamics in general relativity Non-existence of time-periodic dynamics in general relativity Volker Schlue University of Toronto University of Miami, February 2, 2015 Outline 1 General relativity Newtonian mechanics Self-gravitating

More information

Quasi-local mass and isometric embedding

Quasi-local mass and isometric embedding Quasi-local mass and isometric embedding Mu-Tao Wang, Columbia University September 23, 2015, IHP Recent Advances in Mathematical General Relativity Joint work with Po-Ning Chen and Shing-Tung Yau. The

More information

General Relativity and Cosmology Mock exam

General Relativity and Cosmology Mock exam Physikalisches Institut Mock Exam Universität Bonn 29. June 2011 Theoretische Physik SS 2011 General Relativity and Cosmology Mock exam Priv. Doz. Dr. S. Förste Exercise 1: Overview Give short answers

More information

A brief introduction to Semi-Riemannian geometry and general relativity. Hans Ringström

A brief introduction to Semi-Riemannian geometry and general relativity. Hans Ringström A brief introduction to Semi-Riemannian geometry and general relativity Hans Ringström May 5, 2015 2 Contents 1 Scalar product spaces 1 1.1 Scalar products...................................... 1 1.2 Orthonormal

More information

Ask class: what is the Minkowski spacetime in spherical coordinates? ds 2 = dt 2 +dr 2 +r 2 (dθ 2 +sin 2 θdφ 2 ). (1)

Ask class: what is the Minkowski spacetime in spherical coordinates? ds 2 = dt 2 +dr 2 +r 2 (dθ 2 +sin 2 θdφ 2 ). (1) 1 Tensor manipulations One final thing to learn about tensor manipulation is that the metric tensor is what allows you to raise and lower indices. That is, for example, v α = g αβ v β, where again we use

More information

Chapter 2 General Relativity and Black Holes

Chapter 2 General Relativity and Black Holes Chapter 2 General Relativity and Black Holes In this book, black holes frequently appear, so we will describe the simplest black hole, the Schwarzschild black hole and its physics. Roughly speaking, a

More information

A rotating charged black hole solution in f (R) gravity

A rotating charged black hole solution in f (R) gravity PRAMANA c Indian Academy of Sciences Vol. 78, No. 5 journal of May 01 physics pp. 697 703 A rotating charged black hole solution in f R) gravity ALEXIS LARRAÑAGA National Astronomical Observatory, National

More information

Asymptotic Behavior of Marginally Trapped Tubes

Asymptotic Behavior of Marginally Trapped Tubes Asymptotic Behavior of Marginally Trapped Tubes Catherine Williams January 29, 2009 Preliminaries general relativity General relativity says that spacetime is described by a Lorentzian 4-manifold (M, g)

More information

Chapter 7 Curved Spacetime and General Covariance

Chapter 7 Curved Spacetime and General Covariance Chapter 7 Curved Spacetime and General Covariance In this chapter we generalize the discussion of preceding chapters to extend covariance to more general curved spacetimes. 145 146 CHAPTER 7. CURVED SPACETIME

More information

2 General Relativity. 2.1 Curved 2D and 3D space

2 General Relativity. 2.1 Curved 2D and 3D space 22 2 General Relativity The general theory of relativity (Einstein 1915) is the theory of gravity. General relativity ( Einstein s theory ) replaced the previous theory of gravity, Newton s theory. The

More information

Gravitation: Special Relativity

Gravitation: Special Relativity An Introduction to General Relativity Center for Relativistic Astrophysics School of Physics Georgia Institute of Technology Notes based on textbook: Spacetime and Geometry by S.M. Carroll Spring 2013

More information

Lecture Notes on General Relativity

Lecture Notes on General Relativity Lecture Notes on General Relativity Matthias Blau Albert Einstein Center for Fundamental Physics Institut für Theoretische Physik Universität Bern CH-3012 Bern, Switzerland The latest version of these

More information

Black Holes and Thermodynamics I: Classical Black Holes

Black Holes and Thermodynamics I: Classical Black Holes Black Holes and Thermodynamics I: Classical Black Holes Robert M. Wald General references: R.M. Wald General Relativity University of Chicago Press (Chicago, 1984); R.M. Wald Living Rev. Rel. 4, 6 (2001).

More information

BLACK HOLES (ADVANCED GENERAL RELATIV- ITY)

BLACK HOLES (ADVANCED GENERAL RELATIV- ITY) Imperial College London MSc EXAMINATION May 2015 BLACK HOLES (ADVANCED GENERAL RELATIV- ITY) For MSc students, including QFFF students Wednesday, 13th May 2015: 14:00 17:00 Answer Question 1 (40%) and

More information

A relaxation of the strangeness index

A relaxation of the strangeness index echnical report from Automatic Control at Linköpings universitet A relaxation of the strangeness index Henrik idefelt, orkel Glad Division of Automatic Control E-mail: tidefelt@isy.liu.se, torkel@isy.liu.se

More information

An Overview of Mathematical General Relativity

An Overview of Mathematical General Relativity An Overview of Mathematical General Relativity José Natário (Instituto Superior Técnico) Geometria em Lisboa, 8 March 2005 Outline Lorentzian manifolds Einstein s equation The Schwarzschild solution Initial

More information

PAPER 311 BLACK HOLES

PAPER 311 BLACK HOLES MATHEMATICAL TRIPOS Part III Friday, 8 June, 018 9:00 am to 1:00 pm PAPER 311 BLACK HOLES Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal weight. STATIONERY

More information

Rigidity of Black Holes

Rigidity of Black Holes Rigidity of Black Holes Sergiu Klainerman Princeton University February 24, 2011 Rigidity of Black Holes PREAMBLES I, II PREAMBLE I General setting Assume S B two different connected, open, domains and

More information

Mathematical Relativity, Spring 2017/18 Instituto Superior Técnico

Mathematical Relativity, Spring 2017/18 Instituto Superior Técnico Mathematical Relativity, Spring 2017/18 Instituto Superior Técnico 1. Starting from R αβµν Z ν = 2 [α β] Z µ, deduce the components of the Riemann curvature tensor in terms of the Christoffel symbols.

More information

THE INITIAL VALUE FORMULATION OF GENERAL RELATIVITY

THE INITIAL VALUE FORMULATION OF GENERAL RELATIVITY THE INITIAL VALUE FORMULATION OF GENERAL RELATIVITY SAM KAUFMAN Abstract. The (Cauchy) initial value formulation of General Relativity is developed, and the maximal vacuum Cauchy development theorem is

More information

TO GET SCHWARZSCHILD BLACKHOLE SOLUTION USING MATHEMATICA FOR COMPULSORY COURSE WORK PAPER PHY 601

TO GET SCHWARZSCHILD BLACKHOLE SOLUTION USING MATHEMATICA FOR COMPULSORY COURSE WORK PAPER PHY 601 TO GET SCHWARZSCHILD BLACKHOLE SOLUTION USING MATHEMATICA FOR COMPULSORY COURSE WORK PAPER PHY 601 PRESENTED BY: DEOBRAT SINGH RESEARCH SCHOLAR DEPARTMENT OF PHYSICS AND ASTROPHYSICS UNIVERSITY OF DELHI

More information

Introduction to General Relativity

Introduction to General Relativity Introduction to General Relativity Lectures by Igor Pesando Slides by Pietro Fré Virgo Site May 22nd 2006 The issue of reference frames Since oldest and observers antiquity the Who is at motion? The Sun

More information

Derivatives in General Relativity

Derivatives in General Relativity Derivatives in General Relativity One of the problems with curved space is in dealing with vectors how do you add a vector at one point in the surface of a sphere to a vector at a different point, and

More information

Singularity formation in black hole interiors

Singularity formation in black hole interiors Singularity formation in black hole interiors Grigorios Fournodavlos DPMMS, University of Cambridge Heraklion, Crete, 16 May 2018 Outline The Einstein equations Examples Initial value problem Large time

More information

Stationarity of non-radiating spacetimes

Stationarity of non-radiating spacetimes University of Warwick April 4th, 2016 Motivation Theorem Motivation Newtonian gravity: Periodic solutions for two-body system. Einstein gravity: Periodic solutions? At first Post-Newtonian order, Yes!

More information

General Relativity and Differential

General Relativity and Differential Lecture Series on... General Relativity and Differential Geometry CHAD A. MIDDLETON Department of Physics Rhodes College November 1, 2005 OUTLINE Distance in 3D Euclidean Space Distance in 4D Minkowski

More information

Tensor Calculus, Part 2

Tensor Calculus, Part 2 Massachusetts Institute of Technology Department of Physics Physics 8.962 Spring 2002 Tensor Calculus, Part 2 c 2000, 2002 Edmund Bertschinger. 1 Introduction The first set of 8.962 notes, Introduction

More information

Charge, geometry, and effective mass in the Kerr- Newman solution to the Einstein field equations

Charge, geometry, and effective mass in the Kerr- Newman solution to the Einstein field equations Charge, geometry, and effective mass in the Kerr- Newman solution to the Einstein field equations Gerald E. Marsh Argonne National Laboratory (Ret) 5433 East View Park Chicago, IL 60615 E-mail: gemarsh@uchicago.edu

More information

A Summary of the Black Hole Perturbation Theory. Steven Hochman

A Summary of the Black Hole Perturbation Theory. Steven Hochman A Summary of the Black Hole Perturbation Theory Steven Hochman Introduction Many frameworks for doing perturbation theory The two most popular ones Direct examination of the Einstein equations -> Zerilli-Regge-Wheeler

More information

GRAVITATION F10. Lecture Maxwell s Equations in Curved Space-Time 1.1. Recall that Maxwell equations in Lorentz covariant form are.

GRAVITATION F10. Lecture Maxwell s Equations in Curved Space-Time 1.1. Recall that Maxwell equations in Lorentz covariant form are. GRAVITATION F0 S. G. RAJEEV Lecture. Maxwell s Equations in Curved Space-Time.. Recall that Maxwell equations in Lorentz covariant form are. µ F µν = j ν, F µν = µ A ν ν A µ... They follow from the variational

More information

A5682: Introduction to Cosmology Course Notes. 2. General Relativity

A5682: Introduction to Cosmology Course Notes. 2. General Relativity 2. General Relativity Reading: Chapter 3 (sections 3.1 and 3.2) Special Relativity Postulates of theory: 1. There is no state of absolute rest. 2. The speed of light in vacuum is constant, independent

More information

carroll/notes/ has a lot of good notes on GR and links to other pages. General Relativity Philosophy of general

carroll/notes/ has a lot of good notes on GR and links to other pages. General Relativity Philosophy of general http://pancake.uchicago.edu/ carroll/notes/ has a lot of good notes on GR and links to other pages. General Relativity Philosophy of general relativity. As with any major theory in physics, GR has been

More information

Quasi-local Mass and Momentum in General Relativity

Quasi-local Mass and Momentum in General Relativity Quasi-local Mass and Momentum in General Relativity Shing-Tung Yau Harvard University Stephen Hawking s 70th Birthday University of Cambridge, Jan. 7, 2012 I met Stephen Hawking first time in 1978 when

More information

Vectors. January 13, 2013

Vectors. January 13, 2013 Vectors January 13, 2013 The simplest tensors are scalars, which are the measurable quantities of a theory, left invariant by symmetry transformations. By far the most common non-scalars are the vectors,

More information

1 The role of gravity The view of physics that is most generally accepted at the moment is that one can divide the discussion of the universe into

1 The role of gravity The view of physics that is most generally accepted at the moment is that one can divide the discussion of the universe into 1 The role of gravity The view of physics that is most generally accepted at the moment is that one can divide the discussion of the universe into two parts. First, there is the question of the local laws

More information

Solving the Geodesic Equation

Solving the Geodesic Equation Solving the Geodesic Equation Jeremy Atkins December 12, 2018 Abstract We find the general form of the geodesic equation and discuss the closed form relation to find Christoffel symbols. We then show how

More information

CALCULUS ON MANIFOLDS. 1. Riemannian manifolds Recall that for any smooth manifold M, dim M = n, the union T M =

CALCULUS ON MANIFOLDS. 1. Riemannian manifolds Recall that for any smooth manifold M, dim M = n, the union T M = CALCULUS ON MANIFOLDS 1. Riemannian manifolds Recall that for any smooth manifold M, dim M = n, the union T M = a M T am, called the tangent bundle, is itself a smooth manifold, dim T M = 2n. Example 1.

More information

GTR is founded on a Conceptual Mistake And hence Null and Void

GTR is founded on a Conceptual Mistake And hence Null and Void GTR is founded on a Conceptual Mistake And hence Null and Void A critical review of the fundamental basis of General Theory of Relativity shows that a conceptual mistake has been made in the basic postulate

More information

Astronomy 421. Lecture 24: Black Holes

Astronomy 421. Lecture 24: Black Holes Astronomy 421 Lecture 24: Black Holes 1 Outline General Relativity Equivalence Principle and its Consequences The Schwarzschild Metric The Kerr Metric for rotating black holes Black holes Black hole candidates

More information

Einstein Toolkit Workshop. Joshua Faber Apr

Einstein Toolkit Workshop. Joshua Faber Apr Einstein Toolkit Workshop Joshua Faber Apr 05 2012 Outline Space, time, and special relativity The metric tensor and geometry Curvature Geodesics Einstein s equations The Stress-energy tensor 3+1 formalisms

More information

Classical differential geometry of two-dimensional surfaces

Classical differential geometry of two-dimensional surfaces Classical differential geometry of two-dimensional surfaces 1 Basic definitions This section gives an overview of the basic notions of differential geometry for twodimensional surfaces. It follows mainly

More information

Lecture 8: Curved Spaces

Lecture 8: Curved Spaces EPGY Summer Institute Special and General Relativity 2012 Lecture 8: Curved Spaces With the necessity of curved geodesics within regions with significant energy or mass concentrations we need to understand

More information

Singularities and Causal Structure in General Relativity

Singularities and Causal Structure in General Relativity Singularities and Causal Structure in General Relativity Alexander Chen February 16, 2011 1 What are Singularities? Among the many profound predictions of Einstein s general relativity, the prediction

More information

Angular momentum and Killing potentials

Angular momentum and Killing potentials Angular momentum and Killing potentials E. N. Glass a) Physics Department, University of Michigan, Ann Arbor, Michigan 4809 Received 6 April 995; accepted for publication September 995 When the Penrose

More information

OLIVIA MILOJ March 27, 2006 ON THE PENROSE INEQUALITY

OLIVIA MILOJ March 27, 2006 ON THE PENROSE INEQUALITY OLIVIA MILOJ March 27, 2006 ON THE PENROSE INEQUALITY Abstract Penrose presented back in 1973 an argument that any part of the spacetime which contains black holes with event horizons of area A has total

More information

Accelerated Observers

Accelerated Observers Accelerated Observers In the last few lectures, we ve been discussing the implications that the postulates of special relativity have on the physics of our universe. We ve seen how to compute proper times

More information

Fisica Matematica. Stefano Ansoldi. Dipartimento di Matematica e Informatica. Università degli Studi di Udine. Corso di Laurea in Matematica

Fisica Matematica. Stefano Ansoldi. Dipartimento di Matematica e Informatica. Università degli Studi di Udine. Corso di Laurea in Matematica Fisica Matematica Stefano Ansoldi Dipartimento di Matematica e Informatica Università degli Studi di Udine Corso di Laurea in Matematica Anno Accademico 2003/2004 c 2004 Copyright by Stefano Ansoldi and

More information

Level sets of the lapse function in static GR

Level sets of the lapse function in static GR Level sets of the lapse function in static GR Carla Cederbaum Mathematisches Institut Universität Tübingen Auf der Morgenstelle 10 72076 Tübingen, Germany September 4, 2014 Abstract We present a novel

More information

Classification theorem for the static and asymptotically flat Einstein-Maxwell-dilaton spacetimes possessing a photon sphere

Classification theorem for the static and asymptotically flat Einstein-Maxwell-dilaton spacetimes possessing a photon sphere Classification theorem for the static and asymptotically flat Einstein-Maxwell-dilaton spacetimes possessing a photon sphere Boian Lazov and Stoytcho Yazadjiev Varna, 2017 Outline 1 Motivation 2 Preliminaries

More information

General Relativity and Important Physical Quantities

General Relativity and Important Physical Quantities General Relativity and Important Physical Quantities Shing-Tung Yau Harvard University 2nd LeCosPA Symposium December 14, 2015 This talk is based on joint work with Po-Ning Chen and Mu-Tao Wang. Exactly

More information

4.7 The Levi-Civita connection and parallel transport

4.7 The Levi-Civita connection and parallel transport Classnotes: Geometry & Control of Dynamical Systems, M. Kawski. April 21, 2009 138 4.7 The Levi-Civita connection and parallel transport In the earlier investigation, characterizing the shortest curves

More information

Are spacetime horizons higher dimensional sources of energy fields? (The black hole case).

Are spacetime horizons higher dimensional sources of energy fields? (The black hole case). Are spacetime horizons higher dimensional sources of energy fields? (The black hole case). Manasse R. Mbonye Michigan Center for Theoretical Physics Physics Department, University of Michigan, Ann Arbor,

More information

Notes for an Introduction to General Relativity. Fall 2011.

Notes for an Introduction to General Relativity. Fall 2011. Notes for an Introduction to General Relativity. Fall 2011. Manuel Tiglio Center for Scientific Computation and Mathematical Modeling, Department of Physics, Joint Space Sciences Institute. Maryland Center

More information

The Einstein Field Equations

The Einstein Field Equations The Einstein Field Equations on semi-riemannian manifolds, and the Schwarzschild solution Rasmus Leijon VT 2012 Examensarbete, 15hp Kandidatexamen i tillämpad matematik, 180hp Institutionen för matematik

More information

Localizing solutions of the Einstein equations

Localizing solutions of the Einstein equations Localizing solutions of the Einstein equations Richard Schoen UC, Irvine and Stanford University - General Relativity: A Celebration of the 100th Anniversary, IHP - November 20, 2015 Plan of Lecture The

More information

Properties of Traversable Wormholes in Spacetime

Properties of Traversable Wormholes in Spacetime Properties of Traversable Wormholes in Spacetime Vincent Hui Department of Physics, The College of Wooster, Wooster, Ohio 44691, USA. (Dated: May 16, 2018) In this project, the Morris-Thorne metric of

More information

arxiv:gr-qc/ v2 12 Aug 2005

arxiv:gr-qc/ v2 12 Aug 2005 Tetrads in geometrodynamics Alcides Garat 1 1. Instituto de Física, Facultad de Ciencias, Iguá 4225, esq. Mataojo, Montevideo, Uruguay. (December 7th, 2004) A new tetrad is introduced within the framework

More information

The Divergence Myth in Gauss-Bonnet Gravity. William O. Straub Pasadena, California November 11, 2016

The Divergence Myth in Gauss-Bonnet Gravity. William O. Straub Pasadena, California November 11, 2016 The Divergence Myth in Gauss-Bonnet Gravity William O. Straub Pasadena, California 91104 November 11, 2016 Abstract In Riemannian geometry there is a unique combination of the Riemann-Christoffel curvature

More information

General Birkhoff s Theorem

General Birkhoff s Theorem General Birkhoff s Theorem Amir H. Abbassi Department of Physics, School of Sciences, Tarbiat Modarres University, P.O.Box 14155-4838, Tehran, I.R.Iran E-mail: ahabbasi@net1cs.modares.ac.ir Abstract Space-time

More information

Bondi mass of Einstein-Maxwell-Klein-Gordon spacetimes

Bondi mass of Einstein-Maxwell-Klein-Gordon spacetimes of of Institute of Theoretical Physics, Charles University in Prague April 28th, 2014 scholtz@troja.mff.cuni.cz 1 / 45 Outline of 1 2 3 4 5 2 / 45 Energy-momentum in special Lie algebra of the Killing

More information

The Role of Black Holes in the AdS/CFT Correspondence

The Role of Black Holes in the AdS/CFT Correspondence The Role of Black Holes in the AdS/CFT Correspondence Mario Flory 23.07.2013 Mario Flory BHs in AdS/CFT 1 / 30 GR and BHs Part I: General Relativity and Black Holes Einstein Field Equations Lightcones

More information

1 Differentiable manifolds and smooth maps

1 Differentiable manifolds and smooth maps 1 Differentiable manifolds and smooth maps Last updated: April 14, 2011. 1.1 Examples and definitions Roughly, manifolds are sets where one can introduce coordinates. An n-dimensional manifold is a set

More information

Chapter 21. The Kerr solution The Kerr metric in Boyer-Lindquist coordinates

Chapter 21. The Kerr solution The Kerr metric in Boyer-Lindquist coordinates Chapter 21 The Kerr solution As shown in Chapter 10, the solution of Einstein s equations describing the exterior of an isolated, spherically symmetric, static object is quite simple. Indeed, the Schwarzschild

More information

Lecture 1. Introduction

Lecture 1. Introduction Lecture 1. Introduction In Differential Geometry we study spaces which are smooth enough to do calculus. These include such familiar objects as curves and surfaces in space. In its earliest manifestation,

More information

8 Symmetries and the Hamiltonian

8 Symmetries and the Hamiltonian 8 Symmetries and the Hamiltonian Throughout the discussion of black hole thermodynamics, we have always assumed energy = M. Now we will introduce the Hamiltonian formulation of GR and show how to define

More information

Aspects of Black Hole Physics

Aspects of Black Hole Physics Contents Aspects of Black Hole Physics Andreas Vigand Pedersen The Niels Bohr Institute Academic Advisor: Niels Obers e-mail: vigand@nbi.dk Abstract: This project examines some of the exact solutions to

More information

An Introduction to Kaluza-Klein Theory

An Introduction to Kaluza-Klein Theory An Introduction to Kaluza-Klein Theory A. Garrett Lisi nd March Department of Physics, University of California San Diego, La Jolla, CA 993-39 gar@lisi.com Introduction It is the aim of Kaluza-Klein theory

More information

A Brief Introduction to Mathematical Relativity

A Brief Introduction to Mathematical Relativity A Brief Introduction to Mathematical Relativity Arick Shao Imperial College London Arick Shao (Imperial College London) Mathematical Relativity 1 / 31 Special Relativity Postulates and Definitions Einstein

More information

What happens at the horizon of an extreme black hole?

What happens at the horizon of an extreme black hole? What happens at the horizon of an extreme black hole? Harvey Reall DAMTP, Cambridge University Lucietti and HSR arxiv:1208.1437 Lucietti, Murata, HSR and Tanahashi arxiv:1212.2557 Murata, HSR and Tanahashi,

More information

Tensor Analysis in Euclidean Space

Tensor Analysis in Euclidean Space Tensor Analysis in Euclidean Space James Emery Edited: 8/5/2016 Contents 1 Classical Tensor Notation 2 2 Multilinear Functionals 4 3 Operations With Tensors 5 4 The Directional Derivative 5 5 Curvilinear

More information

MATH 4030 Differential Geometry Lecture Notes Part 4 last revised on December 4, Elementary tensor calculus

MATH 4030 Differential Geometry Lecture Notes Part 4 last revised on December 4, Elementary tensor calculus MATH 4030 Differential Geometry Lecture Notes Part 4 last revised on December 4, 205 Elementary tensor calculus We will study in this section some basic multilinear algebra and operations on tensors. Let

More information

Block diagonalization of matrix-valued sum-of-squares programs

Block diagonalization of matrix-valued sum-of-squares programs Block diagonalization of matrix-valued sum-of-squares programs Johan Löfberg Division of Automatic Control Department of Electrical Engineering Linköpings universitet, SE-581 83 Linköping, Sweden WWW:

More information

Modern Physics notes Paul Fendley Lecture 35. Born, chapter III (most of which should be review for you), chapter VII

Modern Physics notes Paul Fendley Lecture 35. Born, chapter III (most of which should be review for you), chapter VII Modern Physics notes Paul Fendley fendley@virginia.edu Lecture 35 Curved spacetime black holes Born, chapter III (most of which should be review for you), chapter VII Fowler, Remarks on General Relativity

More information

Umbilic cylinders in General Relativity or the very weird path of trapped photons

Umbilic cylinders in General Relativity or the very weird path of trapped photons Umbilic cylinders in General Relativity or the very weird path of trapped photons Carla Cederbaum Universität Tübingen European Women in Mathematics @ Schloss Rauischholzhausen 2015 Carla Cederbaum (Tübingen)

More information

INVESTIGATING THE KERR BLACK HOLE USING MAPLE IDAN REGEV. Department of Mathematics, University of Toronto. March 22, 2002.

INVESTIGATING THE KERR BLACK HOLE USING MAPLE IDAN REGEV. Department of Mathematics, University of Toronto. March 22, 2002. INVESTIGATING THE KERR BLACK HOLE USING MAPLE 1 Introduction IDAN REGEV Department of Mathematics, University of Toronto March 22, 2002. 1.1 Why Study the Kerr Black Hole 1.1.1 Overview of Black Holes

More information

Chapter 3. Riemannian Manifolds - I. The subject of this thesis is to extend the combinatorial curve reconstruction approach to curves

Chapter 3. Riemannian Manifolds - I. The subject of this thesis is to extend the combinatorial curve reconstruction approach to curves Chapter 3 Riemannian Manifolds - I The subject of this thesis is to extend the combinatorial curve reconstruction approach to curves embedded in Riemannian manifolds. A Riemannian manifold is an abstraction

More information

arxiv: v2 [gr-qc] 16 Sep 2013

arxiv: v2 [gr-qc] 16 Sep 2013 An algorithm for computing geometric relative velocities through Fermi and observational coordinates Vicente J. Bolós Dpto. Matemáticas para la Economía y la Empresa, Facultad de Economía, Universidad

More information

Kadath: a spectral solver and its application to black hole spacetimes

Kadath: a spectral solver and its application to black hole spacetimes Kadath: a spectral solver and its application to black hole spacetimes Philippe Grandclément Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris F-92195 Meudon, France philippe.grandclement@obspm.fr

More information

Mass, quasi-local mass, and the flow of static metrics

Mass, quasi-local mass, and the flow of static metrics Mass, quasi-local mass, and the flow of static metrics Eric Woolgar Dept of Mathematical and Statistical Sciences University of Alberta ewoolgar@math.ualberta.ca http://www.math.ualberta.ca/~ewoolgar Nov

More information

Schwarschild Metric From Kepler s Law

Schwarschild Metric From Kepler s Law Schwarschild Metric From Kepler s Law Amit kumar Jha Department of Physics, Jamia Millia Islamia Abstract The simplest non-trivial configuration of spacetime in which gravity plays a role is for the region

More information

1 Vector fields, flows and Lie derivatives

1 Vector fields, flows and Lie derivatives Working title: Notes on Lie derivatives and Killing vector fields Author: T. Harmark Killingvectors.tex 5//2008, 22:55 Vector fields, flows and Lie derivatives Coordinates Consider an n-dimensional manifold

More information

Curved Spacetime III Einstein's field equations

Curved Spacetime III Einstein's field equations Curved Spacetime III Einstein's field equations Dr. Naylor Note that in this lecture we will work in SI units: namely c 1 Last Week s class: Curved spacetime II Riemann curvature tensor: This is a tensor

More information

Analytic Kerr Solution for Puncture Evolution

Analytic Kerr Solution for Puncture Evolution Analytic Kerr Solution for Puncture Evolution Jon Allen Maximal slicing of a spacetime with a single Kerr black hole is analyzed. It is shown that for all spin parameters, a limiting hypersurface forms

More information

Relativity Discussion

Relativity Discussion Relativity Discussion 4/19/2007 Jim Emery Einstein and his assistants, Peter Bergmann, and Valentin Bargmann, on there daily walk to the Institute for advanced Study at Princeton. Special Relativity The

More information

A GENERAL RELATIVITY WORKBOOK. Thomas A. Moore. Pomona College. University Science Books. California. Mill Valley,

A GENERAL RELATIVITY WORKBOOK. Thomas A. Moore. Pomona College. University Science Books. California. Mill Valley, A GENERAL RELATIVITY WORKBOOK Thomas A. Moore Pomona College University Science Books Mill Valley, California CONTENTS Preface xv 1. INTRODUCTION 1 Concept Summary 2 Homework Problems 9 General Relativity

More information

Newman-Penrose formalism in higher dimensions

Newman-Penrose formalism in higher dimensions Newman-Penrose formalism in higher dimensions V. Pravda various parts in collaboration with: A. Coley, R. Milson, M. Ortaggio and A. Pravdová Introduction - algebraic classification in four dimensions

More information