Evolution of three-dimensional gravitational waves: Harmonic slicing case

Abstract

We perform numerical simulations of a three-dimensional (3D) time evolution of pure gravitational waves. We use a conformally flat and K = 0 initial condition for the evolution of the spacetime. We adopt several slicing conditions to check whether a long time integration is possible in those conditions. For the case in which the amplitude of the gravitational waves is low, a long time integration is possible by using the harmonic slice and the maximal slice, while in the geodesic slice ( a = 1) it is not possible. As in the axisymmetric case and also in the 3D case, gravitational waves with a sufficiently high amplitude collapse by their self-gravity and their final fates seem to be as black holes. In this case, the singularity avoidance property of the harmonic slice seems weak, so that it may be inappropriate for the formation problems of the black hole. By means of the gauge-invariant wave extraction technique we compute the waveform of the gravitational waves at an outer region. We find that the nonlinearity of Einstein gravity induces the higher multipole modes even if only a quadrupole mode exists initially.PACS number(s): 04.30.Nk

General relativistic magnetohydrodynamics in axisymmetric dynamical spacetimes: the X-ECHO code

“…The 3 + 1 formalism is nowadays adopted in basically all numerical schemes for general relativity (Alcubierre 2008;Baumgarte & Shapiro 2010), where the system of Einstein equations is treated like a Cauchy problem with some initial data to be evolved in time through hyperbolic equations. However, as for the solenoidal condition for the magnetic field, non-evolutionary constraints must be preserved in the numerical evolution, and computational methods for modern codes are divided into two main classes: 1) free-evolution schemes, mainly based on hyperbolic equations alone, where this problem is alleviated by appropriate reformulations of the equations (BSSN: Shibata & Nakamura 1995;Baumgarte & Shapiro 1999), eventually with the addition of propagating modes and damping terms (Z4: Bona et al 2003;Bernuzzi & Hilditch 2010); 2) fully constrained schemes, where the constraints are enforced at each timestep through the solution of elliptic equations (Bonazzola et al 2004), a more robust but computationally demanding option, since elliptic solvers are notoriously difficult to parallelize. Most of the state-of-the-art 3D codes for GRMHD in dynamical spacetimes are based on freeevolution schemes in Cartesian coordinates (Duez et al 2005;Shibata & Sekiguchi 2005;Anderson et al 2006;Giacomazzo & Rezzolla 2007;Montero et al 2008;Farris et al 2008), and have been used for gravitational collapse in the presence of magnetized plasmas Shibata et al 2006a,b;Stephens et al 2007Stephens et al , 2008, evolution of NSs (Duez et al 2006b;Liebling et al 2010), binary NS mergers (Anderson et al 2008;Liu et al 2008;Giacomazzo et al 2009Giacomazzo et al , 2011, and accreting tori around Kerr BHs (Montero et al 2010).…”

Numerical tests of evolution systems, gauge conditions, and boundary conditions for 1D colliding gravitational plane waves

“…These codes were not much more successful than ADM. Non-hyperbolic Baumgarte-Shapiro-Shibata-Nakamura (BSSN) schemes [10,11], based on conformal decomposition of the metric, have shown considerable success in improving the stability of 3D calculations for weak and strong gravitational fields and a variety of spacetime slicings [12]. Alcubierre et al [13] report that a BSSN scheme, combined with excision and certain dynamic gauge conditions, allows accurate numerical evolutions of 3D distorted dynamic black holes up to hundreds of dynamical times.…”

Gravitational Waves from Compact Bodies

“…This enterprise is rather like computing the Lamb shift to high order in powers of the fine structure constant, for comparison with experiment. The terms of leading order in the mass ratio η -μ/M are being checked by a Japanese-American consortium (Poisson, Nakamura, Sasaki, Tagoshi, Tanaka) using the Teukolsky formalism for weak perturbations of black holes (Poisson & Sasaki 1995;Shibata et al 1995). These small-77 calculations have been carried to very high post-Newtonian order for circular orbits and no spins (Tagoshi & Nakamura 1994;Tagoshi & Sasaki 1994), and from those results Cutler & Flanagan (1995) have estimated the order to which the full, finite-7/ computations must be carried in order that systematic errors in the theoretical templates will not significantly impact the information extracted from the LIGO/VIRGO observational data.…”