Tanja Bode

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.

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.

Relativistic Mergers of Supermassive Black Holes and Their Electromagnetic Signatures

Coincident detections of electromagnetic (EM) and gravitational wave (GW) signatures from coalescence events of supermassive black holes (SMBHs) are the next observational grand challenge. Such detections will provide the means to study cosmological evolution and accretion processes associated with these gargantuan compact objects. More generally, the observations will enable testing general relativity in the strong, nonlinear regime and will provide independent cosmological measurements to high precision. Understanding the conditions under which coincidences of EM and GW signatures arise during SMBH mergers is therefore of paramount importance. As an essential step toward this goal, we present results from the first fully general relativistic, hydrodynamical study of the late inspiral and merger of equal-mass, spinning SMBH binaries in a gas cloud. We find that variable EM signatures correlated with GWs can arise in merging systems as a consequence of shocks and accretion combined with the effect of relativistic beaming. The most striking EM variability is observed for systems where spins are aligned with the orbital axis and where orbiting black holes form a stable set of density wakes, but all systems exhibit some characteristic signatures that can be utilized in searches for EM counterparts. In the case of the most massive binaries observable by the Laser Interferometer Space Antenna, calculated luminosities imply that they may be identified by EM searches to z 1, while lower mass systems and binaries immersed in low density ambient gas can only be detected in the local universe.

TIDAL DISRUPTIONS OF WHITE DWARFS FROM ULTRA-CLOSE ENCOUNTERS WITH INTERMEDIATE-MASS SPINNING BLACK HOLES

We present numerical relativity results of tidal disruptions of white dwarfs from ultra-close encounters with a spinning, intermediate-mass black hole (IMBH). These encounters require a full general relativistic treatment of gravity. We show that the disruption process and prompt accretion of the debris strongly depend on the magnitude and orientation of the black hole (BH) spin. However, the late-time accretion onto the BH follows the same decay, Ṁ ∝ t^(-5/3), estimated from Newtonian gravity disruption studies. We compute the spectrum of the disk formed from the fallback material using a slim disk model. The disk spectrum peaks in the soft X-rays and sustains Eddington luminosity for 1-3 yr after the disruption. For arbitrary BH spin orientations, the disrupted material is scattered away from the orbital plane by relativistic frame dragging, which often leads to obscuration of the inner fallback disk by the outflowing debris. The disruption events also yield bursts of gravitational radiation with characteristic frequencies of ~3.2 Hz and strain amplitudes of ~10^(–18) for galactic IMBHs. The optimistic rate of considered ultra-close disruptions is consistent with no sources found in the ROSAT all-sky survey. Future missions like Wide-Field X-ray Telescope could observe dozens of events.

MERGERS OF SUPERMASSIVE BLACK HOLES IN ASTROPHYSICAL ENVIRONMENTS

Modeling the late inspiral and merger of supermassive black holes is central to understanding accretion processes and the conditions under which electromagnetic emission accompanies gravitational waves. We use fully general relativistic, hydrodynamics simulations to investigate how electromagnetic signatures correlate with black hole spins, mass ratios, and the gaseous environment in this final phase of binary evolution. In all scenarios, we find some form of characteristic electromagnetic variability whose pattern depends on the spins and binary mass ratios. Binaries in hot accretion flows exhibit a flare followed by a sudden drop in luminosity associated with the plunge and merger, as well as quasi-periodic oscillations correlated with the gravitational waves during the inspiral. Conversely, circumbinary disk systems are characterized by a low luminosity of variable emission, suggesting challenging prospects for their detection.

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

GRB060218 AS A TIDAL DISRUPTION OF A WHITE DWARF BY AN INTERMEDIATE MASS BLACK HOLE

A highly unusual pair of a gamma-ray burst (GRB) GRB060218 and an associated supernova SN2006aj has puzzled theorists for years. A supernova shock breakout and a jet from a newborn stellar mass compact object were put forward to explain its multiwavelength signature. We propose that the source is naturally explained by another channel, a tidal disruption of a white dwarf (WD) by an intermediate mass black hole (IMBH). The tidal disruption is accompanied by a tidal pinching, which leads to the ignition of a WD and a supernova. Some debris falls back onto the IMBH, forms a disk, which quickly amplifies the magnetic field, and launches a jet. We successfully fit soft X-ray spectrum with the Comptonized blackbody emission from a jet photosphere. The optical/UV emission is consistent with self-absorbed synchrotron from the expanding jet front. The accretion rate temporal dependence _ M(t) in a tidal disruption provides a good fit to soft X-ray lightcurve. The IMBH mass is found to be about 10 4 M in three independent estimates: (1) fitting tidal disruption _ M(t) to soft X-ray lightcurve; (2) computing the jet base radius in a jet photospheric emission model; (3) inferring the central BH mass based on a host dwarf galaxy stellar mass. The supernova position is consistent with the center of the host galaxy, while low supernova ejecta mass is consistent with a WD mass. High expected rate of tidal disruptions in dwarf galaxies is consistent with one source observed by Swift satellite over several years at GRB060218 distance of 150 Mpc. The encounters with the WDs provide a lot of fuel for IMBH growth. Subject headings: accretion ‐ black hole physics ‐ gamma rays: bursts ‐ radiation mechanisms: general ‐ supernovae: general ‐ X-rays: individual (GRB060218)

Late inspiral and merger of binary black holes in scalar–tensor theories of gravity

Gravitational wave observations will probe nonlinear gravitational interactions and thus enable strong tests of Einsteins theory of general relativity. We present a numerical relativity study of the late inspiral and merger of binary black holes in scalar–tensor theories of gravity. We consider binaries inside a scalar field bubble, including in some cases a potential. We demonstrate how an evolving scalar field is able to trigger detectable differences between gravitational waves in scalar–tensor gravity and the corresponding waves in general relativity.

Properties of accretion flows around coalescing supermassive black holes

What are the properties of accretion flows in the vicinity of coalescing supermassive black holes (SBHs)? The answer to this question has direct implications on the feasibility of coincident detections of electromagnetic (EM) and gravitational wave (GW) signals from coalescences. Such detections are considered to be the next observational grand challenge that will enable testing general relativity in the strong, nonlinear regime and improve our understanding of evolution and growth of these massive compact objects. In this paper, we review the properties of the environment of coalescing binaries in the context of the circumbinary disk and hot, radiatively inefficient accretion flow models and use them to mark the extent of the parameter space spanned by this problem. We report the results from an ongoing, general relativistic, hydrodynamical study of the inspiral and merger of black holes, motivated by the latter scenario. We find that correlated EM+GW oscillations can arise during the inspiral phase followed by the gradual rise and subsequent drop-off in the light curve at the time of coalescence. While there are indications that the latter EM signature is a more robust one, a detection of either signal coincidentally with GWs would be a convincing evidence for an impending SBH binary coalescence. The observability of an EM counterpart in the hot accretion flow scenario depends on the details of a model. In the case of the most massive binaries observable by the Laser Interferometer Space Antenna, upper limits on luminosity imply that they may be identified by EM searches out to z ≈ 0.1–1. However, given the radiatively inefficient nature of the gas flow, we speculate that a majority of massive binaries may appear as low luminosity AGN in the local universe.

Binary black hole evolutions of approximate puncture initial data

Approximate solutions to the Einstein field equations are valuable tools to investigate gravitational phenomena. An important aspect of any approximation is to investigate and quantify its regime of validity. We present a study that evaluates the effects that approximate puncture initial data, based on skeleton solutions to the Einstein constraints as proposed by [G. Faye, P. Jaranowski, and G. Schaefer, Phys. Rev. D 69, 124029 (2004).], have on numerical evolutions. Using data analysis tools, we assess the effectiveness of these constraint-violating initial data for both initial and advanced LIGO and show that the matches of waveforms from skeleton data with the corresponding waveforms from constraint-satisfying initial data are > or approx. 0.97 when the total mass of the binary is > or approx. 40M{sub {center_dot}}. In addition, we demonstrate that the differences between the skeleton and the constraint-satisfying initial data evolutions, and thus waveforms, are due to negative Hamiltonian constraint violations present in the skeleton initial data located in the vicinity of the punctures. During the evolution, the skeleton data develops both Hamiltonian and momentum constraint violations that decay with time, with the binary system relaxing to a constraint-satisfying solution with black holes of smaller mass and thus different dynamics.