Accurate Simulations of Binary Black Hole Mergers in Force-free Electrodynamics

We provide additional information on our recent study of the electromagnetic emission produced during the inspiral and merger of supermassive black holes when these are immersed in a force-free plasma threaded by a uniform magnetic field. As anticipated in a recent letter, our results show that although a dual-jet structure is present, the associated luminosity is ~100 times smaller than the total one, which is predominantly quadrupolar. Here we discuss the details of our implementation of the equations in which the force-free condition is not implemented at a discrete level, but rather obtained via a damping scheme which drives the solution to satisfy the correct condition. We show that this is important for a correct and accurate description of the current sheets that can develop in the course of the simulation. We also study in greater detail the three-dimensional charge distribution produced as a consequence of the inspiral and show that during the inspiral it possesses a complex but ordered structure which traces the motion of the two black holes. Finally, we provide quantitative estimates of the scaling of the electromagnetic emission with frequency, with the diffused part having a dependence that is the same as the gravitational-wave one and that scales as L^(non-coll)_(EM) ≈ Ω^((10/3)–(8/3)), while the collimated one scales as L^(coll)_(EM) ≈ Ω^((5/3)–(6/3)), thus with a steeper dependence than previously estimated. We discuss the impact of these results on the potential detectability of dual jets from supermassive black holes and the steps necessary for more accurate estimates

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.Comment: 37 pages, 6 figures, single column. Matches published versio

Black-hole horizons as probes of black-hole dynamics I: post-merger recoil in head-on collisions

The understanding of strong-field dynamics near black-hole horizons is a long-standing and challenging prob- lem in general relativity. Recent advances in numerical relativity and in the geometric characterization of black- hole horizons open new avenues into the problem. In this first paper in a series of two, we focus on the analysis of the recoil occurring in the merger of binary black holes, extending the analysis initiated in [1] with Robinson- Trautman spacetimes. More specifically, we probe spacetime dynamics through the correlation of quantities defined at the black-hole horizon and at null infinity. The geometry of these hypersurfaces responds to bulk gravitational fields acting as test screens in a scattering perspective of spacetime dynamics. Within a 3 + 1 approach we build an effective-curvature vector from the intrinsic geometry of dynamical-horizon sections and correlate its evolution with the flux of Bondi linear momentum at large distances. We employ this setup to study numerically the head-on collision of nonspinning black holes and demonstrate its validity to track the qualita- tive aspects of recoil dynamics at infinity. We also make contact with the suggestion that the antikick can be described in terms of a "slowness parameter" and how this can be computed from the local properties of the horizon. In a companion paper [2] we will further elaborate on the geometric aspects of this approach and on its relation with other approaches to characterize dynamical properties of black-hole horizons.Comment: final version published on PR

Black-hole horizons as probes of black-hole dynamics II: geometrical insights

In a companion paper [1], we have presented a cross-correlation approach to near-horizon physics in which bulk dynamics is probed through the correlation of quantities defined at inner and outer spacetime hypersurfaces acting as test screens. More specifically, dynamical horizons provide appropriate inner screens in a 3+1 setting and, in this context, we have shown that an effective-curvature vector measured at the common horizon produced in a head-on collision merger can be correlated with the flux of linear Bondi-momentum at null infinity. In this paper we provide a more sound geometric basis to this picture. First, we show that a rigidity property of dynamical horizons, namely foliation uniqueness, leads to a preferred class of null tetrads and Weyl scalars on these hypersurfaces. Second, we identify a heuristic horizon news-like function, depending only on the geometry of spatial sections of the horizon. Fluxes constructed from this function offer refined geometric quantities to be correlated with Bondi fluxes at infinity, as well as a contact with the discussion of quasi-local 4-momentum on dynamical horizons. Third, we highlight the importance of tracking the internal horizon dual to the apparent horizon in spatial 3-slices when integrating fluxes along the horizon. Finally, we discuss the link between the dissipation of the non-stationary part of the horizon's geometry with the viscous-fluid analogy for black holes, introducing a geometric prescription for a "slowness parameter" in black-hole recoil dynamics.Comment: Final version published on PR

Density Profiles of Collapsed Rotating Massive Stars Favor Long Gamma-Ray Bursts

Long-duration gamma-ray bursts (lGRBs) originate in relativistic collimated outflows -- jets -- that drill their way out of collapsing massive stars. Accurately modeling this process requires realistic stellar profiles for the jets to propagate through and break out of. Most previous studies have used simple power laws or pre-collapse models for massive stars. However, the relevant stellar profile for lGRB models is in fact that of a star after its core has collapsed to form a compact object. To self-consistently compute such a stellar profile, we use the open-source code GR1D to simulate the core-collapse process for a suite of low-metallicity, rotating, massive stellar progenitors that have undergone chemically homogeneous evolution. Our models span a range of zero-age main sequence (ZAMS) masses: $M_\mathrm{ZAMS} = 13, 18, 21, 25, 35, 40$, and $45 M_\odot$. All of these models, at the onset of core-collapse, feature steep density profiles, $\rho \propto r^{-\alpha}$ with $\alpha\approx 2.5$, which would result in jets that are inconsistent with lGRB observables. We follow the collapse of four out of our seven models until they form BHs and the other three proto-neutron stars (PNSs). We find, across all models, that the density profile outside of the newly-formed BH or PNS is well-represented by a flatter power law with $\alpha \approx 1.35{-}1.55$. Such flat density profiles are conducive to successful formation and breakout of BH-powered jets and, in fact, required to reproduce observable properties of lGRBs. Future models of lGRBs should be initialized with shallower \textit{post-collapse} stellar profiles like those presented here instead of the much steeper pre-collapse profiles that are typically used.Comment: 9 pages, 4 figures+1 table, submitted to ApJL, comments welcom

The Einstein Toolkit

The Einstein Toolkit is a community-driven software platform of core computational tools to advance and support research in relativistic astrophysics and gravitational physics. The Einstein Toolkit has been supported by NSF 2004157/2004044/2004311/2004879/2003893/1550551/1550461/1550436/1550514, Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the National Science Foundation

The Einstein Toolkit

The Einstein Toolkit is a&nbsp;community-driven software platform of core computational tools to advance and support research in relativistic astrophysics and gravitational physics. </span