Numerical Relativity Beyond General Relativity
Total Page:16
File Type:pdf, Size:1020Kb
Load more
Recommended publications
-
Numerical Relativity for Vacuum Binary Black Hole Luisa T
Numerical Relativity for Vacuum Binary Black Hole Luisa T. Buchman UW Bothell SXS collaboration (Caltech) Max Planck Institute for Gravitational Physics Albert Einstein Institute GWA NW @ LIGO Hanford Observatory June 25, 2019 Outline • What is numerical relativity? • What is its role in the detection and interpretation of gravitational waves? What is numerical relativity for binary black holes? • 3 stages for binary black hole coalescence: • inspiral • merger • ringdown • Inspiral waveform: Post-Newtonian approximations • Ringdown waveform: perturbation theory • Merger of 2 black holes: • extremely energetic, nonlinear, dynamical and violent • the strongest GW signal • warpage of spacetime • strong gravity regime The Einstein equations: • Gμν = 8πTμν • encodes all of gravity • 10 equations which would take up ~100 pages if no abbreviations or simplifications used (and if written in terms of the spacetime metric alone) • pencil and paper solutions possible only with spherical (Schwarzschild) or axisymmetric (Kerr) symmetries. • Vacuum: G = 0 μν (empty space with no matter) Numerical Relativity is: • Directly solving the full dynamical Einstein field equations using high-performance computing. • A very difficult problem: • the first attempts at numerical simulations of binary black holes was in the 1960s (Hahn and Lindquist) and mid-1970s (Smarr and others). • the full 3D problem remained unsolved until 2005 (Frans Pretorius, followed quickly by 2 other groups- Campanelli et al., Baker et al.). Solving Einstein’s Equations on a Computer • -
Covariant Hamiltonian Formalism for Field Theory: Hamilton-Jacobi
Covariant hamiltonian formalism for field theory: Hamilton-Jacobi equation on the space G Carlo Rovelli Centre de Physique Th´eorique de Luminy, CNRS, Case 907, F-13288 Marseille, EU February 7, 2008 Abstract Hamiltonian mechanics of field theory can be formulated in a generally covariant and background independent manner over a finite dimensional extended configuration space. The physical symplectic structure of the theory can then be defined over a space G of three-dimensional surfaces without boundary, in the extended configuration space. These surfaces provide a preferred over-coordinatization of phase space. I consider the covariant form of the Hamilton-Jacobi equation on G, and a canonical function S on G which is a preferred solution of the Hamilton-Jacobi equation. The application of this formalism to general relativity is equiv- alent to the ADM formalism, but fully covariant. In the quantum domain, it yields directly the Ashtekar-Wheeler-DeWitt equation. Finally, I apply this formalism to discuss the partial observables of a covariant field theory and the role of the spin networks –basic objects in quantum gravity– in the classical theory. arXiv:gr-qc/0207043v2 26 Jul 2002 1 Introduction Hamiltonian mechanics is a clean and general formalism for describing a physi- cal system, its states and its observables, and provides a road towards quantum theory. In its traditional formulation, however, the hamiltonian formalism is badly non covariant. This is a source of problems already for finite dimensional systems. For instance, the notions of state and observable are not very clean in the hamiltonian mechanics of the systems where evolution is given in para- metric form, especially if the evolution cannot be deparametrized (as in certain cosmological models). -
Linearization Instability in Gravity Theories1
Linearization Instability in Gravity Theories1 Emel Altasa;2 a Department of Physics, Middle East Technical University, 06800, Ankara, Turkey. Abstract In a nonlinear theory, such as gravity, physically relevant solutions are usually hard to find. Therefore, starting from a background exact solution with symmetries, one uses the perturbation theory, which albeit approximately, provides a lot of information regarding a physical solution. But even this approximate information comes with a price: the basic premise of a perturbative solution is that it should be improvable. Namely, by going to higher order perturbation theory, one should be able to improve and better approximate the physical problem or the solution. While this is often the case in many theories and many background solutions, there are important cases where the linear perturbation theory simply fails for various reasons. This issue is well known in the context of general relativity through the works that started in the early 1970s, but it has only been recently studied in modified gravity theories. This thesis is devoted to the study of linearization instability in generic gravity theories where there are spurious solutions to the linearized equations which do not come from the linearization of possible exact solutions. For this purpose we discuss the Taub charges, arXiv:1808.04722v2 [hep-th] 2 Sep 2018 the ADT charges and the quadratic constraints on the linearized solutions. We give the three dimensional chiral gravity and the D dimensional critical gravity as explicit examples and give a detailed ADM analysis of the topologically massive gravity with a cosmological constant. 1This is a Ph.D. -
Wave Extraction in Numerical Relativity
Doctoral Dissertation Wave Extraction in Numerical Relativity Dissertation zur Erlangung des naturwissenschaftlichen Doktorgrades der Bayrischen Julius-Maximilians-Universitat¨ Wurzburg¨ vorgelegt von Oliver Elbracht aus Warendorf Institut fur¨ Theoretische Physik und Astrophysik Fakultat¨ fur¨ Physik und Astronomie Julius-Maximilians-Universitat¨ Wurzburg¨ Wurzburg,¨ August 2009 Eingereicht am: 27. August 2009 bei der Fakultat¨ fur¨ Physik und Astronomie 1. Gutachter:Prof.Dr.Karl Mannheim 2. Gutachter:Prof.Dr.Thomas Trefzger 3. Gutachter:- der Dissertation. 1. Prufer¨ :Prof.Dr.Karl Mannheim 2. Prufer¨ :Prof.Dr.Thomas Trefzger 3. Prufer¨ :Prof.Dr.Thorsten Ohl im Promotionskolloquium. Tag des Promotionskolloquiums: 26. November 2009 Doktorurkunde ausgehandigt¨ am: Gewidmet meinen Eltern, Gertrud und Peter, f¨urall ihre Liebe und Unterst¨utzung. To my parents Gertrud and Peter, for all their love, encouragement and support. Wave Extraction in Numerical Relativity Abstract This work focuses on a fundamental problem in modern numerical rela- tivity: Extracting gravitational waves in a coordinate and gauge independent way to nourish a unique and physically meaningful expression. We adopt a new procedure to extract the physically relevant quantities from the numerically evolved space-time. We introduce a general canonical form for the Weyl scalars in terms of fundamental space-time invariants, and demonstrate how this ap- proach supersedes the explicit definition of a particular null tetrad. As a second objective, we further characterize a particular sub-class of tetrads in the Newman-Penrose formalism: the transverse frames. We establish a new connection between the two major frames for wave extraction: namely the Gram-Schmidt frame, and the quasi-Kinnersley frame. Finally, we study how the expressions for the Weyl scalars depend on the tetrad we choose, in a space-time containing distorted black holes. -
Towards a Gauge Polyvalent Numerical Relativity Code: Numerical Methods, Boundary Conditions and Different Formulations
UNIVERSITAT DE LES ILLES BALEARS Towards a gauge polyvalent numerical relativity code: numerical methods, boundary conditions and different formulations per Carles Bona-Casas Tesi presentada per a l'obtenci´o del t´ıtolde doctor a la Facultat de Ci`encies Departament de F´ısica Dirigida per: Prof. Carles Bona Garcia i Dr. Joan Mass´oBenn`assar "The fact that we live at the bottom of a deep gravity well, on the surface of a gas covered planet going around a nuclear fireball 90 million miles away and think this to be normal is obviously some indication of how skewed our perspective tends to be." Douglas Adams. Acknowledgements I would like to acknowledge everyone who ever taught me something. Specially my supervisors who, at least in the beginning, had to suffer my eyes, puzzling at their faces, saying that I had the impression that what they were explaining to me was some unreachable knowledge. It turns out that there is not such a thing and they have finally managed to make me gain some insight in this world of relativity and computing. They have also infected me the disease of worrying about calculations and stuff that most people wouldn't care about and my hair doesn't seem very happy with it, but I still love them for that. Many thanks to everyone in the UIB relativity group and to all the PhD. students who have shared some time with me. Work is less hard in their company. Special thanks to Carlos, Denis, Sasha and Dana for all their very useful comments and for letting me play with their work and codes. -
Numerical Relativity
Paper presented at the 13th Int. Conf on General Relativity and Gravitation 373 Cordoba, Argentina, 1992: Part 2, Workshop Summaries Numerical relativity Takashi Nakamura Yulcawa Institute for Theoretical Physics, Kyoto University, Kyoto 606, JAPAN In GR13 we heard many reports on recent. progress as well as future plans of detection of gravitational waves. According to these reports (see the report of the workshop on the detection of gravitational waves by Paik in this volume), it is highly probable that the sensitivity of detectors such as laser interferometers and ultra low temperature resonant bars will reach the level of h ~ 10—21 by 1998. in this level we may expect the detection of the gravitational waves from astrophysical sources such as coalescing binary neutron stars once a year or so. Therefore the progress in numerical relativity is urgently required to predict the wave pattern and amplitude of the gravitational waves from realistic astrophysical sources. The time left for numerical relativists is only six years or so although there are so many difficulties in principle as well as in practice. Apart from detection of gravitational waves, numerical relativity itself has a final goal: Solve the Einstein equations numerically for (my initial data as accurately as possible and clarify physics in strong gravity. in GRIIS there were six oral presentations and ll poster papers on recent progress in numerical relativity. i will make a brief review of six oral presenta— tions. The Regge calculus is one of methods to investigate spacetimes numerically. Brewin from Monash University Australia presented a paper Particle Paths in a Schwarzshild Spacetime via. -
The Initial Data Problem for 3+1 Numerical Relativity Part 2
The initial data problem for 3+1 numerical relativity Part 2 Eric Gourgoulhon Laboratoire Univers et Th´eories (LUTH) CNRS / Observatoire de Paris / Universit´eParis Diderot F-92190 Meudon, France [email protected] http://www.luth.obspm.fr/∼luthier/gourgoulhon/ 2008 International Summer School on Computational Methods in Gravitation and Astrophysics Asia Pacific Center for Theoretical Physics, Pohang, Korea 28 July - 1 August 2008 Eric Gourgoulhon (LUTH) Initial data problem 2 / 2 APCTP School, 31 July 2008 1 / 41 Plan 1 Helical symmetry for binary systems 2 Initial data for orbiting binary black holes 3 Initial data for orbiting binary neutron stars 4 Initial data for orbiting black hole - neutron star systems 5 References for lectures 1-3 Eric Gourgoulhon (LUTH) Initial data problem 2 / 2 APCTP School, 31 July 2008 2 / 41 Helical symmetry for binary systems Outline 1 Helical symmetry for binary systems 2 Initial data for orbiting binary black holes 3 Initial data for orbiting binary neutron stars 4 Initial data for orbiting black hole - neutron star systems 5 References for lectures 1-3 Eric Gourgoulhon (LUTH) Initial data problem 2 / 2 APCTP School, 31 July 2008 3 / 41 Helical symmetry for binary systems Helical symmetry for binary systems Physical assumption: when the two objects are sufficiently far apart, the radiation reaction can be neglected ⇒ closed orbits Gravitational radiation reaction circularizes the orbits ⇒ circular orbits Geometrical translation: spacetime possesses some helical symmetry Helical Killing vector ξ: -
3+1 Formalism and Bases of Numerical Relativity
3+1 Formalism and Bases of Numerical Relativity Lecture notes Eric´ Gourgoulhon Laboratoire Univers et Th´eories, UMR 8102 du C.N.R.S., Observatoire de Paris, Universit´eParis 7 arXiv:gr-qc/0703035v1 6 Mar 2007 F-92195 Meudon Cedex, France [email protected] 6 March 2007 2 Contents 1 Introduction 11 2 Geometry of hypersurfaces 15 2.1 Introduction.................................... 15 2.2 Frameworkandnotations . .... 15 2.2.1 Spacetimeandtensorfields . 15 2.2.2 Scalar products and metric duality . ...... 16 2.2.3 Curvaturetensor ............................... 18 2.3 Hypersurfaceembeddedinspacetime . ........ 19 2.3.1 Definition .................................... 19 2.3.2 Normalvector ................................. 21 2.3.3 Intrinsiccurvature . 22 2.3.4 Extrinsiccurvature. 23 2.3.5 Examples: surfaces embedded in the Euclidean space R3 .......... 24 2.4 Spacelikehypersurface . ...... 28 2.4.1 Theorthogonalprojector . 29 2.4.2 Relation between K and n ......................... 31 ∇ 2.4.3 Links between the and D connections. .. .. .. .. .. 32 ∇ 2.5 Gauss-Codazzirelations . ...... 34 2.5.1 Gaussrelation ................................. 34 2.5.2 Codazzirelation ............................... 36 3 Geometry of foliations 39 3.1 Introduction.................................... 39 3.2 Globally hyperbolic spacetimes and foliations . ............. 39 3.2.1 Globally hyperbolic spacetimes . ...... 39 3.2.2 Definition of a foliation . 40 3.3 Foliationkinematics .. .. .. .. .. .. .. .. ..... 41 3.3.1 Lapsefunction ................................. 41 3.3.2 Normal evolution vector . 42 3.3.3 Eulerianobservers ............................. 42 3.3.4 Gradients of n and m ............................. 44 3.3.5 Evolution of the 3-metric . 45 4 CONTENTS 3.3.6 Evolution of the orthogonal projector . ....... 46 3.4 Last part of the 3+1 decomposition of the Riemann tensor . -
Gravitational Waves and Binary Black Holes
Ondes Gravitationnelles, S´eminairePoincar´eXXII (2016) 1 { 51 S´eminairePoincar´e Gravitational Waves and Binary Black Holes Thibault Damour Institut des Hautes Etudes Scientifiques 35 route de Chartres 91440 Bures sur Yvette, France Abstract. The theoretical aspects of the gravitational motion and radiation of binary black hole systems are reviewed. Several of these theoretical aspects (high-order post-Newtonian equations of motion, high-order post-Newtonian gravitational-wave emission formulas, Effective-One-Body formalism, Numerical-Relativity simulations) have played a significant role in allowing one to compute the 250 000 gravitational-wave templates that were used for search- ing and analyzing the first gravitational-wave events from coalescing binary black holes recently reported by the LIGO-Virgo collaboration. 1 Introduction On February 11, 2016 the LIGO-Virgo collaboration announced [1] the quasi- simultaneous observation by the two LIGO interferometers, on September 14, 2015, of the first Gravitational Wave (GW) event, called GW150914. To set the stage, we show in Figure 1 the raw interferometric data of the event GW150914, transcribed in terms of their equivalent dimensionless GW strain amplitude hobs(t) (Hanford raw data on the left, and Livingston raw data on the right). Figure 1: LIGO raw data for GW150914; taken from the talk of Eric Chassande-Mottin (member of the LIGO-Virgo collaboration) given at the April 5, 2016 public conference on \Gravitational Waves and Coalescing Black Holes" (Acad´emiedes Sciences, Paris, France). The important point conveyed by these plots is that the observed signal hobs(t) (which contains both the noise and the putative real GW signal) is randomly fluc- tuating by ±5 × 10−19 over time scales of ∼ 0:1 sec, i.e. -