Bruno C. Mundim

The Einstein Toolkit: a community computational infrastructure for relativistic astrophysics

We describe the Einstein Toolkit, a community-driven, freely accessible computational infrastructure intended for use in numerical relativity, relativistic astrophysics, and other applications. The toolkit, developed by a collaboration involving researchers from multiple institutions around the world, combines a core set of components needed to simulate astrophysical objects such as black holes, compact objects, and collapsing stars, as well as a full suite of analysis tools. The Einstein Toolkit is currently based on the Cactus framework for high-performance computing and the Carpet adaptive mesh refinement driver. It implements spacetime evolution via the BSSN evolution system and general relativistic hydrodynamics in a finite-volume discretization. The toolkit is under continuous development and contains many new code components that have been publicly released for the first time and are described in this paper. We discuss the motivation behind the release of the toolkit, the philosophy underlying its development, and the goals of the project. A summary of the implemented numerical techniques is included, as are results of numerical test covering a variety of sample astrophysical problems.

CIRCUMBINARY MAGNETOHYDRODYNAMIC ACCRETION INTO INSPIRALING BINARY BLACK HOLES

We have simulated the magnetohydrodynamic evolution of a circumbinary disk surrounding an equal-mass binary comprising two non-spinning black holes during the period in which the disk inflow time is comparable to the binary evolution time due to gravitational radiation. Both the changing spacetime and the binary orbital evolution are described by an innovative technique utilizing high-order post-Newtonian approximations. Prior to the beginning of the inspiral, the structure of the circumbinary disk is predicted well by extrapolation from Newtonian results: a gap of roughly two binary separation radii is cleared, and matter piles up at the outer edge of this gap as inflow is retarded by torques exerted by the binary; the accretion rate is roughly half its value at large radius. During inspiral, the inner edge of the disk initially moves inward in coordination with the shrinking binary, but—as the orbital evolution accelerates—the inward motion of the disk edge falls behind the rate of binary compression. In this stage, the binary torque falls substantially, but the accretion rate decreases by only 10%-20%. When the binary separation is tens of gravitational radii, the rest-mass efficiency of disk radiation is a few percent, suggesting that supermassive binary black holes could be very luminous at this stage of their evolution. Inner disk heating is modulated at a beat frequency comparable to the binary orbital frequency. However, a disk with sufficient surface density to be luminous may be optically thick, suppressing periodic modulation of the luminosity.

Error-analysis and comparison to analytical models of numerical waveforms produced by the NRAR Collaboration

The Numerical–Relativity–Analytical–Relativity (NRAR) collaboration is a joint effort between members of the numerical relativity, analytical relativity and gravitational-wave data analysis communities. The goal of the NRAR collaboration is to produce numerical-relativity simulations of compact binaries and use them to develop accurate analytical templates for the LIGO/Virgo Collaboration to use in detecting gravitational-wave signals and extracting astrophysical information from them. We describe the results of the first stage of the NRAR project, which focused on producing an initial set of numerical waveforms from binary black holes with moderate mass ratios and spins, as well as one non-spinning binary configuration which has a mass ratio of 10. All of the numerical waveforms are analysed in a uniform and consistent manner, with numerical errors evaluated using an analysis code created by members of the NRAR collaboration. We compare previously-calibrated, non-precessing analytical waveforms, notably the effective-one-body (EOB) and phenomenological template families, to the newly-produced numerical waveforms. We find that when the binarys total mass is ~100–200M⊙, current EOB and phenomenological models of spinning, non-precessing binary waveforms have overlaps above 99% (for advanced LIGO) with all of the non-precessing-binary numerical waveforms with mass ratios ≤4, when maximizing over binary parameters. This implies that the loss of event rate due to modelling error is below 3%. Moreover, the non-spinning EOB waveforms previously calibrated to five non-spinning waveforms with mass ratio smaller than 6 have overlaps above 99.7% with the numerical waveform with a mass ratio of 10, without even maximizing on the binary parameters.

The NINJA-2 catalog of hybrid post-Newtonian/numerical-relativity waveforms for non-precessing black-hole binaries

The numerical injection analysis (NINJA) project is a collaborative effort between members of the numerical-relativity and gravitational wave data-analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search and parameter-estimation algorithms using numerically generated waveforms and to foster closer collaboration between the numerical-relativity and data-analysis communities. The first NINJA project used only a small number of injections of short numerical-relativity waveforms, which limited its ability to draw quantitative conclusions. The goal of the NINJA-2 project is to overcome these limitations with long post-Newtonian—numerical-relativity hybrid waveforms, large numbers of injections and the use of real detector data. We report on the submission requirements for the NINJA-2 project and the construction of the waveform catalog. Eight numerical-relativity groups have contributed 56 hybrid waveforms consisting of a numerical portion modeling the late inspiral, merger and ringdown stitched to a post-Newtonian portion modeling the early inspiral. We summarize the techniques used by each group in constructing their submissions. We also report on the procedures used to validate these submissions, including examination in the time and frequency domains and comparisons of waveforms from different groups against each other. These procedures have so far considered only the (l, m) = (2, 2) mode. Based on these studies, we judge that the hybrid waveforms are suitable for NINJA-2 studies. We note some of the plans for these investigations.

GRHydro: a new open-source general-relativistic magnetohydrodynamics code for the Einstein toolkit

We present the new general-relativistic magnetohydrodynamics (GRMHD) capabilities of the Einstein toolkit, an open-source community-driven numerical relativity and computational relativistic astrophysics code. The GRMHD extension of the toolkit builds upon previous releases and implements the evolution of relativistic magnetized fluids in the ideal MHD limit in fully dynamical spacetimes using the same shock-capturing techniques previously applied to hydrodynamical evolution. In order to maintain the divergence-free character of the magnetic field, the code implements both constrained transport and hyperbolic divergence cleaning schemes. We present test results for a number of MHD tests in Minkowski and curved spacetimes. Minkowski tests include aligned and oblique planar shocks, cylindrical explosions, magnetic rotors, Alfv´ en waves and advected loops, as well as a set of tests designed to study the response of the divergence cleaning scheme to numerically generated monopoles. We study the code’s performance in curved spacetimes with spherical accretion onto a black hole on a fixed background spacetime

Approximate black hole binary spacetime via asymptotic matching

We construct a fully analytic, general relativistic, nonspinning black hole binary spacetime that approximately solves the vacuum Einstein equations everywhere in space and time for black holes sufficiently well separated. The metric is constructed by asymptotically matching perturbed Schwarzschild metrics near each black hole to a two-body post-Newtonian metric far from them, and a two-body post-Minkowskian metric farther still. Asymptotic matching is done without linearizing about a particular time slice, and thus it is valid dynamically and for all times, provided the binary is sufficiently well separated. This approximate global metric can be used for long dynamical evolutions of relativistic magnetohydrodynamical, circumbinary disks around inspiraling supermassive black holes to study a variety of phenomena.

Advances in Simulations of Generic Black-Hole Binaries

We review some of the recent dramatic developments in the fully nonlinear simulation of generic, highly-precessing, black-hole binaries, and introduce a new approach for generating hybrid post-Newtonian / Numerical waveforms for these challenging systems.

Inspiralling, nonprecessing, spinning black hole binary spacetime via asymptotic matching

We construct a new global, fully analytic, approximate spacetime which accurately describes the dynamics of non-precessing, spinning black hole binaries during the inspiral phase of the relativistic merger process. This approximate solution of the vacuum Einstein’s equations can be obtained by asymptotically matching perturbed Kerr solutions near the two black holes to a post-Newtonian metric valid far from the two black holes. This metric is then matched to a post-Minkowskian metric even farther out in the wave zone. The procedure of asymptotic matching is generalized to be valid on all spatial hypersurfaces, instead of a small group of initial hypersurfaces discussed in previous works. This metric is well suited for long term dynamical simulations of spinning black hole binary spacetimes prior to merger, such as studies of circumbinary gas accretion which requires hundreds of binary orbits.

Initial-data contribution to the error budget of gravitational waves from neutron-star binaries

As numerical calculations of inspiraling neutron-star binaries reach values of accuracy that are comparable with those of black-hole binaries, a fine budgeting of the various sources of error becomes increasingly important. Among such sources, the initial data are normally not accounted for, the rationale being that the error on the initial spacelike hypersurface is always far smaller than the error gained during the evolution. We here consider critically this assumption and perform a comparative analysis of the gravitational waveforms relative to essentially the same physical binary configuration when computed with two different initial-data codes, and then evolved with the same evolution code. More specifically, we consider the evolution of irrotational neutron-star binaries computed either with the pseudospectral code lorene, or with the newly developed finite-difference code cocal; both sets of initial data are subsequently evolved with the high-order-evolution code whiskythc. In this way we find that although global quantities of the system, like the mass and angular momentum, have differences of the order of

Inspiraling black-hole binary spacetimes: Challenges in transitioning from analytical to numerical techniques

\ensuremath{\lesssim}0.02%