Computing Geodesics in Numerical Space Times De Computatione Distantiarum Brevissimarum in Spatiis Quattuor Dimensionum Numeris Descriptis

Total Page:16

File Type:pdf, Size:1020Kb

Computing Geodesics in Numerical Space Times De Computatione Distantiarum Brevissimarum in Spatiis Quattuor Dimensionum Numeris Descriptis Leopold-Franzens University Innsbruck Institute of Computer Science Computing Geodesics in Numerical Space Times De computatione distantiarum brevissimarum in spatiis quattuor dimensionum numeris descriptis Master Thesis Marcel Ritter Supervisor: O. Univ. Prof. Dr. Sabine Schindler March 30, 2010 * In memoriam: O. Univ. Prof. i.R. Dr. sub ausp. Manfred Ritter Who never lost his curiosity, fascination and deep passion for science and teaching even during difficult circumstances. Although he was unable to support my work in this world I know his guiding spirit was always by my side. 2 Figure 1: Geodesics in a Kerr metric. Abstract: Traditionally, tools to visualize geodesics in curved spacetimes of general rel- ativity have been specialized solutions, either tailored to a certain spacetime or limited to certain kinds of numerical data. Utilizing a Fiber Bundle data model and the Vish visualization environment, this thesis aims to solve this problem by developing an approach that is independent of the underlying numerical data. My approach allows the combination of several visualization modules, and opens the possibility of applying the computation and visu- alization of integral lines more readily to other scientific domains. These domains include computational fluid dynamics and medical imaging. Contents 1 Introduction 6 2 Theoretical Background 10 2.1 Differential Geometry . 10 2.1.1 Manifolds and Charts . 10 2.1.2 Curve . 12 2.1.3 Tangential Vector . 12 Transformation . 13 2.1.4 Covector . 13 2.1.5 Tensor Field . 14 2.1.6 Metric . 15 2.1.7 Geodesic Equation and Christoffel Symbols . 16 Geodesic Equation . 16 Christoffel Symbols . 19 2.1.8 Geodesic Deviation and Riemann Tensor . 20 2.1.9 Ricci Tensor and Scalar . 21 2.2 General Relativity . 22 2.2.1 Einstein Field Equation . 22 2.2.2 Schwarzschild Metric . 23 2.2.3 Kerr Metric . 24 2.3 Fluid Dynamics . 26 2.4 Medical Imaging . 28 3 Implementation Concepts 29 3.1 Type Traits . 29 3.2 STL Encapsulation . 33 3.3 Reference Pointers . 33 4 Modeling of Scientific Data 36 4.1 Fiber Bundle Data Model . 37 4.2 The Hierarchy Levels . 38 3 CONTENTS 4 4.2.1 Fiber Bundle . 38 4.2.2 Fiber Slice . 39 4.2.3 Fiber Grid . 40 4.2.4 Fiber Topology . 41 4.2.5 Fiber Representation . 42 4.2.6 Fiber Field . 43 4.3 Simplified Access via Selectors . 46 4.4 Data Examples . 47 4.4.1 Uniform, procedural . 48 4.4.2 Multiblock, Curvilinear . 49 4.4.3 Lines . 50 5 Vish - The Vis(h)ualization Environment 52 5.1 Development Quick Start . 53 5.1.1 Availability and Installation . 53 5.1.2 Source Code Organization and Naming Conventions . 54 5.1.3 Make Files and Compilation . 56 5.2 Scene Network . 57 5.2.1 Modules . 58 5.2.2 Data Transport and Access . 62 Creating new Parameter Types . 62 Fiber Bundle Data Access . 63 5.2.3 Rendering Modules . 65 Geometric Algebra . 66 OpenGL . 67 Using OpenGL in a Rendering Module . 67 5.2.4 Vish Scripts . 72 5.3 Caching . 75 5.4 Data Field Interpolation and Finding Local Coordinates . 77 5.4.1 UniGrid . 79 5.4.2 Multiblock . 80 5.4.3 Curvilinear . 85 Local Coordinates in one Hexahedral Cell . 86 Finding Candidates in the Grid . 89 Summarizing the Steps . 91 5.5 Basic Visualization Modules . 93 5.5.1 Coordinate Grid . 93 5.5.2 Coordinate Grid Box . 94 5.5.3 Uniform Grid Lines . 95 5.5.4 Color Map Legend . 95 5.5.5 Multiblock Outlines . 97 CONTENTS 5 6 Computation and Visualization 99 6.1 Defining Initial Conditions for Integral Lines . 100 6.1.1 Initial Positions, Seed Points . 100 Geometric Point Distributions . 101 Random Point Distribution . 105 Grid Union, Convolution and Transformation . 107 6.1.2 Defining Initial Directions . 111 Grid Subtraction . 111 6.2 Computation . 112 6.2.1 Computing First Order Integration Lines . 121 6.2.2 Computing Second Order Integration Lines . 126 6.3 Rendering Lines . 134 7 Applications 141 7.1 Visualizing Flow of CouetteFlow . 142 7.2 Visualizing Flow and Pressure in a Stirred Fluid Tank . 149 7.3 Visualizing Geodesics in a sampled Schwarzschild Metric . 160 7.4 Visualizing Geodesics in a sampled Kerr Metric . 179 7.5 Fiber Tracking in MRI Data . 205 8 Future Work 209 9 Conclusion 211 10 Fazit 213 A Acknowledgements 221 B Definition of Fiber Bundles 222 C Related Publications 223 Chapter 1 Introduction The objective of this thesis is to provide a visualization framework suitable for exploring numerical spacetimes originating from astrophysical numerical relativity. Numerical relativity is a very active research area. Having immense com- puting power available, numerical methods allow to solve the Einstein field equations, chapter 2.2.1, as is was not possible before. One specific applica- tion is the detection of gravitational waves. Gravitational waves have not yet been directly measured, but a Nobel Prize was given to Hulse and Taylor, [44], for finding a convincing indirect evidence. They measured timings in a binary pulsar. The observed frequency increase can only be explained by energy loss due to the emission of gravitational waves. Thus, the waves them- selves have not yet been measured, but their effect on the emitting systems has been observed. Gravitational waves are ripples in spacetime curvature propagating with the speed of light. \Any mass in nonspherical, nonrecti- linear motion produces gravitational waves (...), but gravitational waves are produced most copiously in events such as the coalescence of two compact stars, the merger of massive black holes, or the big bang." [36] Detectors for gravitational waves have been built in the United States in Livingston, Louisiana, in Hanford, Washington [39], in Europe in Hannover, Germany [29], and near Pisa, Italy [24]. In order to detect a gravitational wave, relative distance changes in the order of 10−22 have to be measured. Numerical analysis of situations involving strong gravity, such as the merging of black holes is of great importance. These simulations can be used to match the actually measured data of the gravitational wave detec- tors. Interactive 3D visualization techniques based purely on numerical data helps to analyze these simulations. This thesis focuses on the visualization of geodesics, the shortest (or longest) path between points in space (or space- time). They are important indicators of the structure of spacetime. 6 CHAPTER 1. INTRODUCTION 7 Geodesics in curved spacetimes have been studied before, but most work is related to the visualization of analytic spacetimes, 1992 [2], 1997 [26], 1999 [22], 2001 [15] and 2004 [28]. Computing geodesics is required for ray-tracing black holes. In 1991 Corvin Zahn at the University of T¨ubingenimplemented four dimensional ray-tracing in an analytic Schwarzschild space time, [63]. The most similar work was done already in 1992 by Steve Bryson, who im- plemented geodesic visualizations for exploration using a \boom mounted six degree of freedom head position sensitive stereo CRT system", [17]. For his setup the curved spacetimes could be given by closed formulas or also on simple uniform grids. Geodesics were also analyzed in numerical spacetimes in the 2D (axisym- metric) era. The famous \pair of pants" picture contains the event horizon in a head-on collision, along with some geodesics. The corresponding movies were made in 1995 with great effort in TV resolution, see e.g. [43]. Werner Benger at the University of Innsbruck simulated a black hole by raytracing in 1996 [4]. Andrew Hamilton implemented a real time flight simulator for a charged black hole. He uses a projective technique to compute the paths of geodesics, [35]. Spacetimes visualized in this thesis are sampled on uniform grids. How- ever, the developed infrastructure used in this work easily extends to adaptive mesh refinement (AMR) grids, which are currently used for numerical simu- lations of merging and colliding black holes, without changing the computa- tion and visualization algorithms. Moreover, the rudimentary visualization of geodesics can be enhanced with other visualization modules, for instance to show the coordinate-acceleration, equation (2.53), on the geodesics. Inspired by Bryson's seeding methods for geodesics I developed a flexible technique for the creation of seeding geometries using basic operating blocks based on the theory of fiber bundles, that can be combined to a huge variety of seeding geometries. In addition to the visualization and computation, other aspects have been addressed in this work: the data model high code re-usability high modularity provide an introduction to utilized software environments and libraries The modern scientific world uses many computational methods in different scientific domains. The requirement of collaborations increases and thus CHAPTER 1. INTRODUCTION 8 exchanging data becomes important. During my research visit at Louisiana State University I was taught that in 2005 hurricane Katrina forced scientific groups to exchange their data sets and couple different kinds of numeric simulations to predict the path of Katrina to provide warning and rescue to the public. For example, actual satellite data was coupled with simulation of the motion of air of the hurricane, which was then coupled with a simulation for wave propagation on the ocean. Such coupling is only possible if data can be exchanged without time-consuming data conversion processes between the different scientific groups. The simulation results were later gathered in a visualization combining several layers stemming from all the different domains, [14]. The thesis starts with a short theoretical part providing the necessary mathematical and physical background of general relativity, computational fluid dynamics and magnetic resonance imaging needed for the visualization and interpretation of the resulting illustrations, chapter 2. Thereafter, some C++ programming techniques used for implementation are presented in chapter 3.
Recommended publications
  • A Mathematical Derivation of the General Relativistic Schwarzschild
    A Mathematical Derivation of the General Relativistic Schwarzschild Metric An Honors thesis presented to the faculty of the Departments of Physics and Mathematics East Tennessee State University In partial fulfillment of the requirements for the Honors Scholar and Honors-in-Discipline Programs for a Bachelor of Science in Physics and Mathematics by David Simpson April 2007 Robert Gardner, Ph.D. Mark Giroux, Ph.D. Keywords: differential geometry, general relativity, Schwarzschild metric, black holes ABSTRACT The Mathematical Derivation of the General Relativistic Schwarzschild Metric by David Simpson We briefly discuss some underlying principles of special and general relativity with the focus on a more geometric interpretation. We outline Einstein’s Equations which describes the geometry of spacetime due to the influence of mass, and from there derive the Schwarzschild metric. The metric relies on the curvature of spacetime to provide a means of measuring invariant spacetime intervals around an isolated, static, and spherically symmetric mass M, which could represent a star or a black hole. In the derivation, we suggest a concise mathematical line of reasoning to evaluate the large number of cumbersome equations involved which was not found elsewhere in our survey of the literature. 2 CONTENTS ABSTRACT ................................. 2 1 Introduction to Relativity ...................... 4 1.1 Minkowski Space ....................... 6 1.2 What is a black hole? ..................... 11 1.3 Geodesics and Christoffel Symbols ............. 14 2 Einstein’s Field Equations and Requirements for a Solution .17 2.1 Einstein’s Field Equations .................. 20 3 Derivation of the Schwarzschild Metric .............. 21 3.1 Evaluation of the Christoffel Symbols .......... 25 3.2 Ricci Tensor Components .................
    [Show full text]
  • Connections on Bundles Md
    Dhaka Univ. J. Sci. 60(2): 191-195, 2012 (July) Connections on Bundles Md. Showkat Ali, Md. Mirazul Islam, Farzana Nasrin, Md. Abu Hanif Sarkar and Tanzia Zerin Khan Department of Mathematics, University of Dhaka, Dhaka 1000, Bangladesh, Email: [email protected] Received on 25. 05. 2011.Accepted for Publication on 15. 12. 2011 Abstract This paper is a survey of the basic theory of connection on bundles. A connection on tangent bundle , is called an affine connection on an -dimensional smooth manifold . By the general discussion of affine connection on vector bundles that necessarily exists on which is compatible with tensors. I. Introduction = < , > (2) In order to differentiate sections of a vector bundle [5] or where <, > represents the pairing between and ∗. vector fields on a manifold we need to introduce a Then is a section of , called the absolute differential structure called the connection on a vector bundle. For quotient or the covariant derivative of the section along . example, an affine connection is a structure attached to a differentiable manifold so that we can differentiate its Theorem 1. A connection always exists on a vector bundle. tensor fields. We first introduce the general theorem of Proof. Choose a coordinate covering { }∈ of . Since connections on vector bundles. Then we study the tangent vector bundles are trivial locally, we may assume that there is bundle. is a -dimensional vector bundle determine local frame field for any . By the local structure of intrinsically by the differentiable structure [8] of an - connections, we need only construct a × matrix on dimensional smooth manifold . each such that the matrices satisfy II.
    [Show full text]
  • Schwarzschild Solution to Einstein's General Relativity
    Schwarzschild Solution to Einstein's General Relativity Carson Blinn May 17, 2017 Contents 1 Introduction 1 1.1 Tensor Notations . .1 1.2 Manifolds . .2 2 Lorentz Transformations 4 2.1 Characteristic equations . .4 2.2 Minkowski Space . .6 3 Derivation of the Einstein Field Equations 7 3.1 Calculus of Variations . .7 3.2 Einstein-Hilbert Action . .9 3.3 Calculation of the Variation of the Ricci Tenosr and Scalar . 10 3.4 The Einstein Equations . 11 4 Derivation of Schwarzschild Metric 12 4.1 Assumptions . 12 4.2 Deriving the Christoffel Symbols . 12 4.3 The Ricci Tensor . 15 4.4 The Ricci Scalar . 19 4.5 Substituting into the Einstein Equations . 19 4.6 Solving and substituting into the metric . 20 5 Conclusion 21 References 22 Abstract This paper is intended as a very brief review of General Relativity for those who do not want to skimp on the details of the mathemat- ics behind how the theory works. This paper mainly uses [2], [3], [4], and [6] as a basis, and in addition contains short references to more in-depth references such as [1], [5], [7], and [8] when more depth was needed. With an introduction to manifolds and notation, special rel- ativity can be constructed which are the relativistic equations of flat space-time. After flat space-time, the Lagrangian and calculus of vari- ations will be introduced to construct the Einstein-Hilbert action to derive the Einstein field equations. With the field equations at hand the Schwarzschild equation will fall out with a few assumptions.
    [Show full text]
  • 2015, Volume 8
    V O L U M E 8 2015 D E PAUL UNIVERSITY Creating Knowledge THE LAS JOURNAL OF UNDERGRADUATE SCHOLARSHIP CREATING KNOWLEDGE The LAS Journal of Undergraduate Scholarship 2015 EDITOR Warren C. Schultz ART JURORS Adam Schreiber, Coordinator Laura Kina Steve Harp COPY EDITORS Stephanie Klein Rachel Pomeroy Anastasia Sasewich TABLE OF CONTENTS 6 Foreword, by Interim Dean Lucy Rinehart, PhD STUDENT RESEARCH 8 S. Clelia Sweeney Probing the Public Wound: The Serial Killer Character in True- Crime Media (American Studies Program) 18 Claire Potter Key Progressions: An Examination of Current Student Perspectives of Music School (Department of Anthropology) 32 Jeff Gwizdalski Effect of the Affordable Care Act on Insurance Coverage for Young Adults (Department of Economics) 40 Sam Okrasinski “The Difference of Woman’s Destiny”: Female Friendships as the Element of Change in Jane Austen’s Emma (Department of English) 48 Anna Fechtor Les Musulmans LGBTQ en Europe Occidentale : une communauté non reconnue (French Program, Department of Modern Languages) 58 Marc Zaparaniuk Brazil: A Stadium All Its Own (Department of Geography) 68 Erin Clancy Authority in Stone: Forging the New Jerusalem in Ethiopia (Department of the History of Art and Architecture) 76 Kristin Masterson Emmett J. Scott’s “Official History” of the African-American Experience in World War One: Negotiating Race on the National and International Stage (Department of History) 84 Lizbeth Sanchez Heroes and Victims: The Strategic Mobilization of Mothers during the 1980s Contra War (Department
    [Show full text]
  • Review Study on “The Black Hole”
    IJIRST –International Journal for Innovative Research in Science & Technology| Volume 2 | Issue 10 | March 2016 ISSN (online): 2349-6010 Review Study on “The Black Hole” Syed G. Ibrahim Department of Engineering Physics (Nanostructured Thin Film Materials Laboratory) Prof. Ram Meghe College of Engineering and Management, Badnera 444701, Maharashtra, India Abstract As a star grows old, swells, then collapses on itself, often you will hear the word “black hole” thrown around. The black hole is a gravitationally collapsed mass, from which no light, matter, or signal of any kind can escape. These exotic objects have captured our imagination ever since they were predicted by Einstein's Theory of General Relativity in 1915. So what exactly is a black hole? A black hole is what remains when a massive star dies. Not every star will become a black hole, only a select few with extremely large masses. In order to have the ability to become a black hole, a star will have to have about 20 times the mass of our Sun. No known process currently active in the universe can form black holes of less than stellar mass. This is because all present black hole formation is through gravitational collapse, and the smallest mass which can collapse to form a black hole produces a hole approximately 1.5-3.0 times the mass of the sun .Smaller masses collapse to form white dwarf stars or neutron stars. Keywords: Escape Velocity, Horizon, Schwarzschild Radius, Black Hole _______________________________________________________________________________________________________ I. INTRODUCTION Soon after Albert Einstein formulated theory of relativity, it was realized that his equations have solutions in closed form.
    [Show full text]
  • The Schwarzschild Metric and Applications 1
    The Schwarzschild Metric and Applications 1 Analytic solutions of Einstein's equations are hard to come by. It's easier in situations that e hibit symmetries. 1916: Karl Schwarzschild sought the metric describing the static, spherically symmetric spacetime surrounding a spherically symmetric mass distribution. A static spacetime is one for which there exists a time coordinate t such that i' all the components of g are independent of t ii' the line element ds( is invariant under the transformation t -t A spacetime that satis+es (i) but not (ii' is called stationary. An example is a rotating azimuthally symmetric mass distribution. The metric for a static spacetime has the form where xi are the spatial coordinates and dl( is a time*independent spatial metric. -ross-terms dt dxi are missing because their presence would violate condition (ii'. 23ote: The Kerr metric, which describes the spacetime outside a rotating ( axisymmetric mass distribution, contains a term ∝ dt d.] To preser)e spherical symmetry& dl( can be distorted from the flat-space metric only in the radial direction. In 5at space, (1) r is the distance from the origin and (2) 6r( is the area of a sphere. Let's de+ne r such that (2) remains true but (1) can be violated. Then, A,xi' A,r) in cases of spherical symmetry. The Ricci tensor for this metric is diagonal, with components S/ 10.1 /rimes denote differentiation with respect to r. The region outside the spherically symmetric mass distribution is empty. 9 The vacuum Einstein equations are R = 0. To find A,r' and B,r'# (.
    [Show full text]
  • Oil Workers' Rights Will Not Be Undermined, Minister Says
    SUBSCRIPTION WEDNESDAY, JUNE 22, 2016 RAMADAN 17, 1437 AH www.kuwaittimes.net Alghanim to More than 700 England Iran: Bahrain Ramadan TImings expand Wendy’s doctors killed progress after ‘will pay price’ Emsak: 03:04 brand into Saudi in Syrian war Slovakia draw for crackdown Fajer: 03:14 Shrooq: 04:48 Dohr: 11:50 Asr: 15:24 Maghreb: 18:51 2 7 19 13 Eshaa: 20:23 Oil workers’ rights will not Min 33º Max 47º be undermined, Minister says High Tide 02:55 & 12:12 Low Tide MPs call for new strategy for Kuwaiti investments 06:55 & 19:57 40 PAGES NO: 16912 150 FILS KUWAIT: Finance Minister and Acting Oil Minister Anas Ramadan Kareem Al-Saleh said yesterday the government will not under- mine the rights of the workers in the oil sector who Social dimension went on a three-day strike in April to demand their rights. The minister said that the oil sector proposed a of fasting, Ramadan number of initiatives to cut spending in light of the sharp drop in oil prices and the workers went on strike By Hatem Basha even before these initiatives were applied. ost people view Ramadan as a sublime peri- He said the Ministry is currently in talks with the od, where every Muslim experience sublime workers about the initiatives but insisted that the Mfeelings as hearts soften and souls tran- rights of the workers will not be touched. The minis- scend the worldly pleasures. Another dimension that ter’s comments came in response to remarks made by is often overlooked is the social dimension of fasting.
    [Show full text]
  • Friday Prime Time, April 17 4 P.M
    April 17 - 23, 2009 SPANISH FORK CABLE GUIDE 9 Friday Prime Time, April 17 4 P.M. 4:30 5 P.M. 5:30 6 P.M. 6:30 7 P.M. 7:30 8 P.M. 8:30 9 P.M. 9:30 10 P.M. 10:30 11 P.M. 11:30 BASIC CABLE Oprah Winfrey Å 4 News (N) Å CBS Evening News (N) Å Entertainment Ghost Whisperer “Save Our Flashpoint “First in Line” ’ NUMB3RS “Jack of All Trades” News (N) Å (10:35) Late Show With David Late Late Show KUTV 2 News-Couric Tonight Souls” ’ Å 4 Å 4 ’ Å 4 Letterman (N) ’ 4 KJZZ 3The People’s Court (N) 4 The Insider 4 Frasier ’ 4 Friends ’ 4 Friends 5 Fortune Jeopardy! 3 Dr. Phil ’ Å 4 News (N) Å Scrubs ’ 5 Scrubs ’ 5 Entertain The Insider 4 The Ellen DeGeneres Show (N) News (N) World News- News (N) Two and a Half Wife Swap “Burroughs/Padovan- Supernanny “DeMello Family” 20/20 ’ Å 4 News (N) (10:35) Night- Access Holly- (11:36) Extra KTVX 4’ Å 3 Gibson Men 5 Hickman” (N) ’ 4 (N) ’ Å line (N) 3 wood (N) 4 (N) Å 4 News (N) Å News (N) Å News (N) Å NBC Nightly News (N) Å News (N) Å Howie Do It Howie Do It Dateline NBC A police of cer looks into the disappearance of a News (N) Å (10:35) The Tonight Show With Late Night- KSL 5 News (N) 3 (N) ’ Å (N) ’ Å Michigan woman. (N) ’ Å Jay Leno ’ Å 5 Jimmy Fallon TBS 6Raymond Friends ’ 5 Seinfeld ’ 4 Seinfeld ’ 4 Family Guy 5 Family Guy 5 ‘Happy Gilmore’ (PG-13, ’96) ›› Adam Sandler.
    [Show full text]
  • LECTURE 6: FIBER BUNDLES in This Section We Will Introduce The
    LECTURE 6: FIBER BUNDLES In this section we will introduce the interesting class of fibrations given by fiber bundles. Fiber bundles play an important role in many geometric contexts. For example, the Grassmaniann varieties and certain fiber bundles associated to Stiefel varieties are central in the classification of vector bundles over (nice) spaces. The fact that fiber bundles are examples of Serre fibrations follows from Theorem ?? which states that being a Serre fibration is a local property. 1. Fiber bundles and principal bundles Definition 6.1. A fiber bundle with fiber F is a map p: E ! X with the following property: every ∼ −1 point x 2 X has a neighborhood U ⊆ X for which there is a homeomorphism φU : U × F = p (U) such that the following diagram commutes in which π1 : U × F ! U is the projection on the first factor: φ U × F U / p−1(U) ∼= π1 p * U t Remark 6.2. The projection X × F ! X is an example of a fiber bundle: it is called the trivial bundle over X with fiber F . By definition, a fiber bundle is a map which is `locally' homeomorphic to a trivial bundle. The homeomorphism φU in the definition is a local trivialization of the bundle, or a trivialization over U. Let us begin with an interesting subclass. A fiber bundle whose fiber F is a discrete space is (by definition) a covering projection (with fiber F ). For example, the exponential map R ! S1 is a covering projection with fiber Z. Suppose X is a space which is path-connected and locally simply connected (in fact, the weaker condition of being semi-locally simply connected would be enough for the following construction).
    [Show full text]
  • Cracking the Einstein Code: Relativity and the Birth of Black Hole Physics, with an Afterword by Roy Kerr / Fulvio Melia
    CRA C K I N G T H E E INSTEIN CODE @SZObWdWbgO\RbVS0W`bV]T0ZOQY6]ZS>VgaWQa eWbVO\/TbS`e]`RPg@]gS`` fulvio melia The University of Chicago Press chicago and london fulvio melia is a professor in the departments of physics and astronomy at the University of Arizona. He is the author of The Galactic Supermassive Black Hole; The Black Hole at the Center of Our Galaxy; The Edge of Infinity; and Electrodynamics, and he is series editor of the book series Theoretical Astrophysics published by the University of Chicago Press. The University of Chicago Press, Chicago 60637 The University of Chicago Press, Ltd., London © 2009 by The University of Chicago All rights reserved. Published 2009 Printed in the United States of America 18 17 16 15 14 13 12 11 10 09 1 2 3 4 5 isbn-13: 978-0-226-51951-7 (cloth) isbn-10: 0-226-51951-1 (cloth) Library of Congress Cataloging-in-Publication Data Melia, Fulvio. Cracking the Einstein code: relativity and the birth of black hole physics, with an afterword by Roy Kerr / Fulvio Melia. p. cm. Includes bibliographical references and index. isbn-13: 978-0-226-51951-7 (cloth: alk. paper) isbn-10: 0-226-51951-1 (cloth: alk. paper) 1. Einstein field equations. 2. Kerr, R. P. (Roy P.). 3. Kerr black holes—Mathematical models. 4. Black holes (Astronomy)—Mathematical models. I. Title. qc173.6.m434 2009 530.11—dc22 2008044006 To natalina panaia and cesare melia, in loving memory CONTENTS preface ix 1 Einstein’s Code 1 2 Space and Time 5 3 Gravity 15 4 Four Pillars and a Prayer 24 5 An Unbreakable Code 39 6 Roy Kerr 54 7 The Kerr Solution 69 8 Black Hole 82 9 The Tower 100 10 New Zealand 105 11 Kerr in the Cosmos 111 12 Future Breakthrough 121 afterword 125 references 129 index 133 PREFACE Something quite remarkable arrived in my mail during the summer of 2004.
    [Show full text]
  • Spacetimes with Singularities Ovidiu Cristinel Stoica
    Spacetimes with Singularities Ovidiu Cristinel Stoica To cite this version: Ovidiu Cristinel Stoica. Spacetimes with Singularities. An. Stiint. Univ. Ovidius Constanta, Ser. Mat., 2012, pp.26. hal-00617027v2 HAL Id: hal-00617027 https://hal.archives-ouvertes.fr/hal-00617027v2 Submitted on 30 Nov 2014 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Spacetimes with Singularities ∗†‡ Ovidiu-Cristinel Stoica Abstract We report on some advances made in the problem of singularities in general relativity. First is introduced the singular semi-Riemannian geometry for met- rics which can change their signature (in particular be degenerate). The standard operations like covariant contraction, covariant derivative, and constructions like the Riemann curvature are usually prohibited by the fact that the metric is not invertible. The things become even worse at the points where the signature changes. We show that we can still do many of these operations, in a different framework which we propose. This allows the writing of an equivalent form of Einstein’s equation, which works for degenerate metric too. Once we make the singularities manageable from mathematical view- point, we can extend analytically the black hole solutions and then choose from the maximal extensions globally hyperbolic regions.
    [Show full text]
  • The Emergence of Gravitational Wave Science: 100 Years of Development of Mathematical Theory, Detectors, Numerical Algorithms, and Data Analysis Tools
    BULLETIN (New Series) OF THE AMERICAN MATHEMATICAL SOCIETY Volume 53, Number 4, October 2016, Pages 513–554 http://dx.doi.org/10.1090/bull/1544 Article electronically published on August 2, 2016 THE EMERGENCE OF GRAVITATIONAL WAVE SCIENCE: 100 YEARS OF DEVELOPMENT OF MATHEMATICAL THEORY, DETECTORS, NUMERICAL ALGORITHMS, AND DATA ANALYSIS TOOLS MICHAEL HOLST, OLIVIER SARBACH, MANUEL TIGLIO, AND MICHELE VALLISNERI In memory of Sergio Dain Abstract. On September 14, 2015, the newly upgraded Laser Interferometer Gravitational-wave Observatory (LIGO) recorded a loud gravitational-wave (GW) signal, emitted a billion light-years away by a coalescing binary of two stellar-mass black holes. The detection was announced in February 2016, in time for the hundredth anniversary of Einstein’s prediction of GWs within the theory of general relativity (GR). The signal represents the first direct detec- tion of GWs, the first observation of a black-hole binary, and the first test of GR in its strong-field, high-velocity, nonlinear regime. In the remainder of its first observing run, LIGO observed two more signals from black-hole bina- ries, one moderately loud, another at the boundary of statistical significance. The detections mark the end of a decades-long quest and the beginning of GW astronomy: finally, we are able to probe the unseen, electromagnetically dark Universe by listening to it. In this article, we present a short historical overview of GW science: this young discipline combines GR, arguably the crowning achievement of classical physics, with record-setting, ultra-low-noise laser interferometry, and with some of the most powerful developments in the theory of differential geometry, partial differential equations, high-performance computation, numerical analysis, signal processing, statistical inference, and data science.
    [Show full text]