Arxiv:Gr-Qc/0604015 V1 4 Apr 2006 Ste Motn Olo O Etrwy Oetatmore Extract to It Ways 7]
Total Page:16
File Type:pdf, Size:1020Kb
AEI-2006-018, LSU-REL-033006 Introduction to dynamical horizons in numerical relativity 1,2, 2,† 2,‡ Erik Schnetter, ∗ Badri Krishnan, and Florian Beyer 1Center for Computation and Technology, 302 Johnston Hall, Louisiana State University, Baton Rouge, LA 70803, USA§ 2Max-Planck-Institut f¨ur Gravitationsphysik, Albert-Einstein-Institut, Am M¨uhlenberg 1, D-14476 Golm, Germany¶ (Dated: April 4, 2006) This paper presents a quasi-local method of studying the physics of dynamical black holes in nu- merical simulations. This is done within the dynamical horizon framework, which extends the earlier work on isolated horizons to time-dependent situations. In particular: (i) We locate various kinds of marginal surfaces and study their time evolution. An important ingredient is the calculation of the signature of the horizon, which can be either spacelike, timelike, or null. (ii) We generalize the calcula- tion of the black hole mass and angular momentum, which were previously defined for axisymmetric isolated horizons to dynamical situations. (iii) We calculate the source multipole moments of the black hole which can be used to verify that the black hole settles down to a Kerr solution. (iv) We also study the fluxes of energy crossing the horizon, which describes how a black hole grows as it accretes matter and/or radiation. We describe our numerical implementation of these concepts and apply them to three specific test cases, namely, the axisymmetric head-on collision of two black holes, the axisymmetric collapse of a neutron star, and a non-axisymmetric black hole collision with non-zero initial orbital angular mo- mentum. PACS numbers: 04.25.Dm, 04.70.Bw, 95.30.Sf, 97.60.Lf, I. INTRODUCTION physical information from simulations and to compare results from two different simulations performed using different coordinate systems, gauge conditions etc. This In spite of fundamental advances in our understand- can be a non-trivial task in itself, and understanding dy- ing of black holes, relatively little is known about them namical black holes is a necessary ingredient. in the fully non-perturbative, dynamical regime of gen- eral relativity. Most of our intuition regarding black The dynamics of apparent and event horizons have holes comes from studying the stationary, axisymmet- been numerically studied in the past in detail in axisym- ric Kerr-Newman solutions, and perturbations thereof. metry (see e.g. [8, 9, 10, 11, 12, 13]). We want to extend This, along with post-Newtonian calculations which this work to non-axisymmetric and non-vacuum space- treat the black hole as a point particle, are usually ade- times, and we want to emphasise non-gauge-dependent quate for understanding many astrophysical processes analysis methods. involving black holes. However, understanding the We use the formalism of dynamical horizons [14, 15] gravitational waveforms arising due to, say, the merger to study black holes. Dynamical horizons are a sig- phase of the coalescence of two black holes or the grav- nificant extension of the isolated horizon framework itational collapse of a star, will require us to go beyond [16, 17, 18, 19, 20], which models isolated stationary perturbation theory and to confront the non-linearities black holes in an otherwise dynamical spacetime. Both arXiv:gr-qc/0604015 v1 4 Apr 2006 and dynamics of the full Einstein equations. This regime these frameworks are, in turn, very closely related to may contain qualitatively new, non-perturbative fea- and motivated by the earlier work on trapping horizons tures. In this paper, we discuss an important ingredient by Hayward [21, 22, 23]. See [24, 25, 26] for reviews. The for understanding this regime, namely, the dynamics of basic aim is to study black holes quasi-locally without the black hole horizon. Numerical simulations of black reference to event horizons. holes have greatly improved in the last few years. Sim- Since the event horizon is a global, teleological con- ulations of the entire merger process, starting from the cept, we have to know the entire history of the space- last few orbits of the inspiral right up to the ringdown time in order to locate it. This implies that event hori- havebecomepossible inthe pastyear[1,2,3, 4, 5, 6,7]. It zons can never be observed experimentally and, from is then important to look for better ways to extract more the numerical relativity perspective, there is no known quasi-local condition that can be used to locate the in- tersection of the event horizon with a given Cauchy sur- face. In numerical simulations, it is instead common to ∗Electronic address: [email protected] † use marginally trapped surfaces to locate black holes on Electronic address: [email protected] a Cauchy surface in real time. Similarly, in the theoreti- ‡Electronic address: fl[email protected] §URL: http://www.cct.lsu.edu/about/focus/numerical/ cal works on quasi-local horizons mentioned above, it ¶URL: http://numrel.aei.mpg.de/; URL: http://www.aei.mpg. is shown that marginally trapped surfaces, while not de/ a substitute for event horizons, do have many useful 2 properties and can be used fruitfully to study black hole II. BASIC NOTIONS AND DEFINITIONS physics. Thus, it is likely that ideas and results from dynamical/trapping horizons can be very useful for nu- A. Trapped surfaces and apparent horizons merical relativity. Information obtained from the quasi- local horizons complements the information obtained Let S be a closed, orientable spacelike 2-surface in from the event horizon. Once a simulation is complete a 4-dimensional spacetime ( , gab). The expansion of and ready for post-processing, event horizons are use- any such surface can be definedM invariantly without any ful for studying global properties and the causal struc- reference to a time slicing of the spacetime. Since S is ture of the spacetime, and also phenomena such as the smooth, spacelike, and 2-dimensional, the set of vec- topology change of the horizon during a black hole co- tors orthogonal to it at any point form a 2-dimensional alescence. Reliable and computationally efficient codes Minkowskian vector space. Thus, we can define two lin- are now available for locating event horizons (see e.g. early independent, future-directed, null vectors ℓa and [27]). Such information cannot be obtained at the quasi- na orthogonal to S such that local level, which is instead better for tracking the physi- cal parameters and geometry of a black hole in real time. g ℓanb = 1. (2.1) ab − In particular, we consider the following applications: Note that this convention is different from that used in (i) We study the behavior of various marginally trapped [15]. We shall assume that we know a priori what the surfaces under time evolution. This leads to greater outgoing and ingoing directions on are. By conven- insights about the trapped region of a spacetime. An tion, ℓa will denote an outgoing nullM normal and na an important ingredient here is the signature of the world ingoing one. The null normals are specified only up to a tube of marginally trapped surfaces. This world tube is boost transformation known to be null for isolated horizons, and more gener- a a a 1 a ally, it can be either spacelike or timelike; we show that ℓ f ℓ , n f − n (2.2) → → both types occur frequently in numerical simulations. (ii) We give meaningful definitions for the angular mo- where f is a, positive definite, smooth function on S. All mentum, mass, and higher multipole moments for the physical quantities must be invariant under this gauge dynamical black hole. The multipole moments capture transformation. gauge invariant geometrical information regarding the The Riemannian 2-metric q˜ab on S induced by the horizon geometry, and should be useful for understand- spacetime metric gab is ing fundamental issues such as the final state of black q˜ = g + ℓ n + n ℓ . (2.3) hole collapse. For example, we would expect that after ab ab a b a b a black hole has formed and settled down, its multipole The tensor q˜b can be viewed as a projection operator on moments should be identical to the source multipoles of a to S. The null expansions are a Kerr black hole. We show that it is, in principle, pos- sible to verify this conjecture and to calculate the rate at ab ab Θ ℓ = q˜ ℓ , Θ = q˜ n . (2.4) which a black hole approaches equilibrium. (iii) We also ( ) ∇a b (n) ∇a b describe and implement methods for calculating the en- These expansions tell us how the area element of S ergy flux falling into the horizon. This gives us detailed changes as it is deformed along ℓa and na respectively. information on how black holes grow as they swallow ℓa The shear of , σ(ℓ)ab, is the symmetric trace-free part matter and radiation. of the projection of ℓ : ∇a b This paper is organized as follows. Section II sets c d 1 σ ℓ = q˜ q˜ ℓ Θ ℓ q˜ . (2.5) up notation, and summarizes the basic definitions and ( )ab a b∇(c d) − 2 ( ) ab properties of trapped surfaces and dynamical horizons. a Section III describes the various physical quantities that Similarly, the shear of n is we calculate using dynamical horizons, and also their 1 c d Θ numerical implementation. Section IV presents three σ(n)ab = q˜aq˜b (cnd) (n)q˜ab . (2.6) concrete, well known numerical examples where these ∇ − 2 concepts are applied and finally, section V discusses Note that these definitions only involve derivatives tan- some open issues and directions for further work. Un- gential to S.