Source author record

Roland Haas

Roland Haas appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2020arXiv

Multipolar Effective-One-Body Waveforms for Precessing Binary Black Holes: Construction and Validation

As gravitational-wave detectors become more sensitive, we will access a greater variety of signals emitted by compact binary systems, shedding light on their astrophysical origin and environment. A key physical effect that can distinguish among formation scenarios is the misalignment of the spins with the orbital angular momentum, causing the spins and the binary's orbital plane to precess. To accurately model such systems, it is crucial to include multipoles beyond the dominant quadrupole. Here, we develop the first multipolar precessing waveform model in the effective-one-body (EOB) formalism for the inspiral, merger and ringdown (IMR) of binary black holes: SEOBNRv4PHM. In the nonprecessing limit, the model reduces to SEOBNRv4HM, which was calibrated to numerical-relativity (NR) simulations, and waveforms from perturbation theory. We validate SEOBNRv4PHM by comparing it to the public catalog of 1405 precessing NR waveforms of the Simulating eXtreme Spacetimes (SXS) collaboration, and also to new 118 precessing NR waveforms, which span mass ratios 1-4 and spins up to 0.9. We stress that SEOBNRv4PHM is not calibrated to NR simulations in the precessing sector. We compute the unfaithfulness against the 1523 SXS precessing NR waveforms, and find that, for $94\%$ ($57\%$) of the cases, the maximum value, in the total mass range $20-200 M_\odot$, is below $3\%$ ($1\%$). Those numbers become $83\%$ ($20\%$) when using the IMR, multipolar, precessing phenomenological model IMRPhenomPv3HM. We investigate the impact of such unfaithfulness values with two parameter-estimation studies on synthetic signals. We also compute the unfaithfulness between those waveform models and identify in which part of the parameter space they differ the most. We validate them also against the multipolar, precessing NR surrogate model NRSur7dq4, and find that the SEOBNRv4PHM model outperforms IMRPhenomPv3HM.

preprint2019arXiv

Unequal Mass Binary Neutron Star Simulations with Neutrino Transport: Ejecta and Neutrino Emission

We present twelve new simulations of unequal mass neutron star mergers. The simulations were preformed with the SpEC code, and utilize nuclear-theory based equations of state and a two-moment gray neutrino transport scheme with an improved energy estimate based on evolving the number density. We model the neutron stars with the SFHo, LS220 and DD2 equations of state (EOS) and we study the neutrino and matter emission of all twelve models to search for robust trends between binary parameters and emission characteristics. We find that the total mass of the dynamical ejecta exceeds $0.01M_\odot$ only for SFHo with weak dependence on the mass-ratio across all models. We find that the ejecta have a broad electron fraction ($Y_e$) distribution ($\approx 0.06-0.48$), with mean $0.2$. $Y_e$ increases with neutrino irradiation over time, but decreases with increasing binary asymmetry. We also find that the models have ejecta with a broad asymptotic velocity distribution ($\approx 0.05-0.7c$). The average velocity lies in the range $0.2c - 0.3c$ and decreases with binary asymmetry. Furthermore, we find that disk mass increases with binary asymmetry and stiffness of the EOS. The $Y_e$ of the disk increases with softness of the EOS. The strongest neutrino emission occurs for the models with soft EOS. For (anti) electron neutrinos we find no significant dependence of the magnitude or angular distribution or neutrino luminosity with mass-ratio. The heavier neutrino species have a luminosity dependence on mass-ratio but an angular distribution which does not change with mass-ratio.

preprint2016arXiv

Gravitational waveforms for neutron star binaries from binary black hole simulations

Gravitational waves from binary neutron star (BNS) and black hole/neutron star (BHNS) inspirals are primary sources for detection by the Advanced Laser Interferometer Gravitational-Wave Observatory. The tidal forces acting on the neutron stars induce changes in the phase evolution of the gravitational waveform, and these changes can be used to constrain the nuclear equation of state. Current methods of generating BNS and BHNS waveforms rely on either computationally challenging full 3D hydrodynamical simulations or approximate analytic solutions. We introduce a new method for computing inspiral waveforms for BNS/BHNS systems by adding the post-Newtonian (PN) tidal effects to full numerical simulations of binary black holes (BBHs), effectively replacing the nontidal terms in the PN expansion with BBH results. Comparing a waveform generated with this method against a full hydrodynamical simulation of a BNS inspiral yields a phase difference of $<1$ radian over $\sim 15$ orbits. The numerical phase accuracy required of BNS simulations to measure the accuracy of the method we present here is estimated as a function of the tidal deformability parameter $λ$.

preprint2016arXiv

Low mass binary neutron star mergers : gravitational waves and neutrino emission

Neutron star mergers are among the most promising sources of gravitational waves for advanced ground-based detectors. These mergers are also expected to power bright electromagnetic signals, in the form of short gamma-ray bursts, infrared/optical transients, and radio emission. Simulations of these mergers with fully general relativistic codes are critical to understand the merger and post-merger gravitational wave signals and their neutrinos and electromagnetic counterparts. In this paper, we employ the SpEC code to simulate the merger of low-mass neutron star binaries (two $1.2M_\odot$ neutron stars) for a set of three nuclear-theory based, finite temperature equations of state. We show that the frequency peaks of the post-merger gravitational wave signal are in good agreement with predictions obtained from simulations using a simpler treatment of gravity. We find, however, that only the fundamental mode of the remnant is excited for long periods of time: emission at the secondary peaks is damped on a millisecond timescale in the simulated binaries. For such low-mass systems, the remnant is a massive neutron star which, depending on the equation of state, is either permanently stable or long-lived. We observe strong excitations of l=2, m=2 modes, both in the massive neutron star and in the form of hot, shocked tidal arms in the surrounding accretion torus. We estimate the neutrino emission of the remnant using a neutrino leakage scheme and, in one case, compare these results with a gray two-moment neutrino transport scheme. We confirm the complex geometry of the neutrino emission, also observed in previous simulations with neutrino leakage, and show explicitly the presence of important differences in the neutrino luminosity, disk composition, and outflow properties between the neutrino leakage and transport schemes.

preprint2016arXiv

Neutrino-Driven Convection in Core-Collapse Supernovae: High-Resolution Simulations

We present results from high-resolution semi-global simulations of neutrino-driven convection in core-collapse supernovae. We employ an idealized setup with parametrized neutrino heating/cooling and nuclear dissociation at the shock front. We study the internal dynamics of neutrino-driven convection and its role in re-distributing energy and momentum through the gain region. We find that even if buoyant plumes are able to locally transfer heat up to the shock, convection is not able to create a net positive energy flux and overcome the downwards transport of energy from the accretion flow. Turbulent convection does, however, provide a significant effective pressure support to the accretion flow as it favors the accumulation of energy, mass and momentum in the gain region. We derive an approximate equation that is able to explain and predict the shock evolution in terms of integrals of quantities such as the turbulent pressure in the gain region or the effects of non-radial motion of the fluid. We use this relation as a way to quantify the role of turbulence in the dynamics of the accretion shock. Finally, we investigate the effects of grid resolution, which we change by a factor 20 between the lowest and highest resolution. Our results show that the shallow slopes of the turbulent kinetic energy spectra reported in previous studies are a numerical artefact. Kolmogorov scaling is progressively recovered as the resolution is increased.

preprint2016arXiv

Simulations of inspiraling and merging double neutron stars using the Spectral Einstein Code

We present results on the inspiral, merger, and post-merger evolution of a neutron star - neutron star (NSNS) system. Our results are obtained using the hybrid pseudospectral-finite volume Spectral Einstein Code (SpEC). To test our numerical methods, we evolve an equal-mass system for $\approx 22$ orbits before merger. This waveform is the longest waveform obtained from fully general-relativistic simulations for NSNSs to date. Such long (and accurate) numerical waveforms are required to further improve semi-analytical models used in gravitational wave data analysis, for example the effective one body models. We discuss in detail the improvements to SpEC's ability to simulate NSNS mergers, in particular mesh refined grids to better resolve the merger and post-merger phases. We provide a set of consistency checks and compare our results to NSNS merger simulations with the independent BAM code. We find agreement between them, which increases confidence in results obtained with either code. This work paves the way for future studies using long waveforms and more complex microphysical descriptions of neutron star matter in SpEC.

preprint2015arXiv

Binary Neutron Stars with Arbitrary Spins in Numerical Relativity

We present a code to construct initial data for binary neutron star systems in which the stars are rotating. Our code, based on a formalism developed by Tichy, allows for arbitrary rotation axes of the neutron stars and is able to achieve rotation rates near rotational breakup. We compute the neutron star angular momentum through quasi-local angular momentum integrals. When constructing irrotational binary neutron stars, we find a very small residual dimensionless spin of $\sim 2\times 10^{-4}$. Evolutions of rotating neutron star binaries show that the magnitude of the stars' angular momentum is conserved, and that the spin- and orbit-precession of the stars is well described by post-Newtonian approximation. We demonstrate that orbital eccentricity of the binary neutron stars can be controlled to $\sim 0.1\%$. The neutron stars show quasi-normal mode oscillations at an amplitude which increases with the rotation rate of the stars.

preprint2015arXiv

IllinoisGRMHD: An Open-Source, User-Friendly GRMHD Code for Dynamical Spacetimes

In the extreme violence of merger and mass accretion, compact objects like black holes and neutron stars are thought to launch some of the most luminous outbursts of electromagnetic and gravitational wave energy in the Universe. Modeling these systems realistically is a central problem in theoretical astrophysics, but has proven extremely challenging, requiring the development of numerical relativity codes that solve Einstein's equations for the spacetime, coupled to the equations of general relativistic (ideal) magnetohydrodynamics (GRMHD) for the magnetized fluids. Over the past decade, the Illinois Numerical Relativity (ILNR) Group's dynamical spacetime GRMHD code has proven itself as a robust and reliable tool for theoretical modeling of such GRMHD phenomena. However, the code was written "by experts and for experts" of the code, with a steep learning curve that would severely hinder community adoption if it were open-sourced. Here we present IllinoisGRMHD, which is an open-source, highly-extensible rewrite of the original closed-source GRMHD code of the ILNR Group. Reducing the learning curve was the primary focus of this rewrite, with the goal of facilitating community involvement in the code's use and development, as well as the minimization of human effort in generating new science. IllinoisGRMHD also saves computer time, generating roundoff-precision identical output to the original code on adaptive-mesh grids, but nearly twice as fast at scales of hundreds to thousands of cores.

preprint2015arXiv

Post-merger evolution of a neutron star-black hole binary with neutrino transport

We present a first simulation of the post-merger evolution of a black hole-neutron star binary in full general relativity using an energy-integrated general relativistic truncated moment formalism for neutrino transport. We describe our implementation of the moment formalism and important tests of our code, before studying the formation phase of a disk after a black hole-neutron star merger. We use as initial data an existing general relativistic simulation of the merger of a neutron star of 1.4 solar mass with a black hole of 7 solar mass and dimensionless spin a/M=0.8. Comparing with a simpler leakage scheme for the treatment of the neutrinos, we find noticeable differences in the neutron to proton ratio in and around the disk, and in the neutrino luminosity. We find that the electron neutrino luminosity is much lower in the transport simulations, and that the remnant is less neutron-rich. The spatial distribution of the neutrinos is significantly affected by relativistic effects. Over the short timescale evolved, we do not observe purely neutrino-driven outflows. However, a small amount of material (3e-4Msun) is ejected in the polar region during the circularization of the disk. Most of that material is ejected early in the formation of the disk, and is fairly neutron rich. Through r-process nucleosynthesis, that material should produce high-opacity lanthanides in the polar region, and could thus affect the lightcurve of radioactively powered electromagnetic transients. We also show that by the end of the simulation, while the bulk of the disk is neutron-rich, its outer layers have a higher electron fraction. As that material would be the first to be unbound by disk outflows on longer timescales, the changes in Ye experienced during the formation of the disk could have an impact on the nucleosynthesis outputs from neutrino-driven and viscously-driven outflows. [Abridged]

preprint2014arXiv

Neutron star-black hole mergers with a nuclear equation of state and neutrino cooling: Dependence in the binary parameters

We present a first exploration of the results of neutron star-black hole mergers using black hole masses in the most likely range of $7M_\odot-10M_\odot$, a neutrino leakage scheme, and a modeling of the neutron star material through a finite-temperature nuclear-theory based equation of state. In the range of black hole spins in which the neutron star is tidally disrupted ($χ_{\rm BH}\gtrsim 0.7$), we show that the merger consistently produces large amounts of cool ($T\lesssim 1\,{\rm MeV}$), unbound, neutron-rich material ($M_{\rm ej}\sim 0.05M_\odot-0.20M_\odot$). A comparable amount of bound matter is initially divided between a hot disk ($T_{\rm max}\sim 15\,{\rm MeV}$) with typical neutrino luminosity $L_ν\sim 10^{53}\,{\rm erg/s}$, and a cooler tidal tail. After a short period of rapid protonization of the disk lasting $\sim 10\,{\rm ms}$, the accretion disk cools down under the combined effects of the fall-back of cool material from the tail, continued accretion of the hottest material onto the black hole, and neutrino emission. As the temperature decreases, the disk progressively becomes more neutron-rich, with dimmer neutrino emission. This cooling process should stop once the viscous heating in the disk (not included in our simulations) balances the cooling. These mergers of neutron star-black hole binaries with black hole masses $M_{\rm BH}\sim 7M_\odot-10M_\odot$ and black hole spins high enough for the neutron star to disrupt provide promising candidates for the production of short gamma-ray bursts, of bright infrared post-merger signals due to the radioactive decay of unbound material, and of large amounts of r-process nuclei.

preprint2013arXiv

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 binary's total mass is ~100-200 solar masses, 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.

preprint2013arXiv

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 Mdot(t) in a tidal disruption provides a good fit to soft X-ray lightcurve. The IMBH mass is found to be about 10^4Msun in three independent estimates: (1) fitting tidal disruption Mdot(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 150Mpc. The encounters with the WDs provide a lot of fuel for IMBH growth.

preprint2013arXiv

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 magnetised 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 hyperbolic divergence cleaning and constrained transport 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én 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 and in fully dynamical spacetimes by evolutions of a magnetised polytropic neutron star and of the collapse of a magnetised stellar core. Our results agree well with exact solutions where these are available and we demonstrate convergence. All code and input files used to generate the results are available on http://einsteintoolkit.org. This makes our work fully reproducible and provides new users with an introduction to applications of the code.

preprint2013arXiv

Scalar self-force for eccentric orbits around a Schwarzschild black hole

We revisit the problem of computing the self-force on a scalar charge moving along an eccentric geodesic orbit around a Schwarzschild black hole. This work extends previous scalar self-force calculations for circular orbits, which were based on a regular "effective" point-particle source and a full 3D evolution code. We find good agreement between our results and previous calculations based on a (1+1) time-domain code. Finally, our data visualization is unconventional: we plot the self-force through full radial cycles to create "self-force loops", which reveal many interesting features that are less apparent in standard presentations of eccentric-orbit self-force data.

preprint2013arXiv

Self-force as a cosmic censor

We examine Hubeny's scenario according to which a near-extremal Reissner-Nordström black hole can absorb a particle and be driven toward an over-extremal state in which the charge exceeds the mass, signaling the destruction of the black hole. Our analysis incorporates the particle's electromagnetic self-force and the energy radiated to infinity in the form of electromagnetic waves. With these essential ingredients, our sampling of the parameter space reveals no instances of an overcharged final state, and we conjecture that the self-force acts as a cosmic censor, preventing the destruction of a near-extremal black hole by the absorption of a charged particle. We argue, on the basis of the third law of black-hole mechanics, that this conclusion is robust and should apply to attempts to overspin a Kerr black hole.

preprint2012arXiv

Prompt Emission from Tidal Disruptions of White Dwarfs by Intermediate Mass Black Holes

We present a qualitative picture of prompt emission from tidal disruptions of white dwarfs (WD) by intermediate mass black holes (IMBH). The smaller size of an IMBH compared to a supermassive black hole and a smaller tidal radius of a WD disruption lead to a very fast event with high peak luminosity. Magnetic field is generated in situ following the tidal disruption, which leads to effective accretion. Since large-scale magnetic field is also produced, geometrically thick super-Eddington inflow leads to a relativistic jet. The dense jet possesses a photosphere, which emits quasi-thermal radiation in soft X-rays. The source can be classified as a long low-luminosity gamma-ray burst (ll-GRB). Tidal compression of a WD causes nuclear ignition, which is observable as an accompanying supernova. We suggest that GRB060218 and SN2006aj is such a pair of ll-GRB and supernova. We argue that in a flux-limited sample the disruptions of WDs by IMBHs are more frequent then the disruptions of other stars by IMBHs.

preprint2012arXiv

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. 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 spin. However, the late-time accretion onto the black hole follows the same decay, $\dot{M}$ ~ 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 yrs after the disruption. For arbitrary black hole 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 intermediate mass black holes. The optimistic rate of considered ultra-close disruptions is consistent with no sources found in ROSAT all-sky survey. The future missions like Wide-Field X-ray Telescope (WFXT) could observe dozens of events.

preprint2011arXiv

Late Inspiral and Merger of Binary Black Holes in Scalar-Tensor Theories of Gravity

Gravitational wave observations will probe non-linear gravitational interactions and thus enable strong tests of Einstein's 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 black hole binaries in an inhomogeneous scalar field, specifically binaries inside a scalar field bubble, in some cases with a potential. We calculate the emission of dipole radiation. We also show how these configurations trigger detectable differences between gravitational waves in scalar-tensor gravity and the corresponding waves in general relativity. We conclude that, barring an external mechanism to induce dynamics in the scalar field, scalar-tensor gravity binary black holes alone are not capable of awaking a dormant scalar field, and are thus observationally indistinguishable from their general relativistic counterparts.

preprint2011arXiv

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.

preprint2011arXiv

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 for 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.

preprint2011arXiv

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 article. 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.

preprint2011arXiv

Time domain calculation of the electromagnetic self-force on eccentric geodesics in Schwarzschild spacetime

I calculate the self-force acting on a particle with electric charge q moving on a generic geodesic around a Schwarzschild black hole. Using methods similar to those developed for the scalar field case discussed in a previous paper, I investigate the relative sizes of the conservative (half-advanced plus half-retarded) and dissipative (half-advanced minus half-retarded) pieces of the self-force. I also display the regularization parameters used in the mode-sum regularization scheme.

preprint2010arXiv

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 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 supermassive black hole mergers is therefore of paramount importance. As an essential step towards this goal, we present results from the first fully general relativistic, hydrodynamical study of the late inspiral and merger of equal-mass, spinning supermassive black hole 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.