Rio: A new computational framework for accurate initial data of binary black holes

Abstract

We present a computational framework (Rio) in the ADM 3+1 approach for numerical relativity. This work enables us to carry out high resolution calculations for initial data of two arbitrary black holes. We use the transverse conformal treatment, the Bowen-York and the puncture methods. For the numerical solution of the Hamiltonian constraint we use the domain decomposition and the spectral decomposition of Galerkin-Collocation. The nonlinear numerical code solves the set of equations for the spectral modes using the standard Newton-Raphson method, LU decomposition and Gaussian quadratures. We show the convergence of the Rio code. This code allows for easy deployment of large calculations. We show how the spin of one of the black holes is manifest in the conformal factor.

arXiv comment15 pages, 6 figures; updated version to appear in Gen. Rel. Grav

References in corpus (14)

• Observation of Gravitational Waves from a Binary Black Hole Merger
• GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral
• Multi-messenger Observations of a Binary Neutron Star Merger
• GW151226: Observation of Gravitational Waves from a 22-Solar-Mass Binary Black Hole Coalescence
• GW170104: Observation of a 50-Solar-Mass Binary Black Hole Coalescence at Redshift 0.2
• GW170814: A Three-Detector Observation of Gravitational Waves from a Binary Black Hole Coalescence
• A single-domain spectral method for black hole puncture data
• Initial Data for Numerical Relativity
• The RIT binary black hole simulations catalog
• Binary Black Hole Mergers in 3d Numerical Relativity
• Gravitational wave spectroscopy of binary neutron star merger remnants with mode stacking
• Observing binary black hole ringdowns by advanced gravitational wave detectors
• Bowen-York Type Initial Data for Binaries with Neutron Stars
• Puncture black hole initial data: a single domain Galerkin-Collocation method for trumpet and wormhole data sets

Cited by in corpus (5)

• NRPyElliptic: A Fast Hyperbolic Relaxation Elliptic Solver for Numerical Relativity, I: Conformally Flat, Binary Puncture Initial Data
• The Bondi problem revisited: a spectral domain decomposition code
• Gravitational collapse in AdS: instabilities, turbulence, and information
• Multidomain Galerkin-Collocation method: spherical collapse of scalar fields II
• New numerical framework for the generalized Baumgarte-Shapiro-Shibata-Nakamura formulation: The vacuum case for spherical symmetry