Showing 1–50 of 74 results for author: Campanelli, M

A Parameter Study of the Electromagnetic Signatures of an Analytical Mini-Disk Model for Supermassive Binary Black Hole Systems

Authors: Kaitlyn Porter, Scott C. Noble, Eduardo M. Gutierrez, Joaquin Pelle, Manuela Campanelli, Jeremy Schnittman, Bernard J. Kelly

Abstract: Supermassive black holes (SMBHs) are thought to be located at the centers of most galactic nuclei. When galaxies merge they form supermassive black hole binary (SMBHB) systems and these central SMBHs will also merge at later times, producing gravitational waves (GWs). Because galaxy mergers are likely gas-rich environments, SMBHBs are also potential sources of electromagnetic (EM) radiation. The E… ▽ More Supermassive black holes (SMBHs) are thought to be located at the centers of most galactic nuclei. When galaxies merge they form supermassive black hole binary (SMBHB) systems and these central SMBHs will also merge at later times, producing gravitational waves (GWs). Because galaxy mergers are likely gas-rich environments, SMBHBs are also potential sources of electromagnetic (EM) radiation. The EM signatures depend on gas dynamics, orbital dynamics, and radiation processes. The gas dynamics are governed by general relativistic magnetohydrodynamics (MHD) in a time-dependent spacetime. Numerically solving the MHD equations for a time-dependent binary spacetime is computationally expensive. Therefore, it is challenging to conduct a full exploration of the parameter space of these systems and the resulting EM signatures. We have developed an analytical accretion disk model for the mini-disks of an SMBHB system and produced images and light curves using a general relativistic ray-tracing code and a superimposed harmonic binary black hole metric. This analytical model greatly reduces the time and computational resources needed to explore these systems, while incorporating some key information from simulations. We present a parameter space exploration of the SMBHB system in which we have studied the dependence of the EM signatures on the spins of the black holes (BHs), the mass ratio, the accretion rate, the viewing angle, and the initial binary separation. Additionally, we study how the commonly used fast-light approximation affects the EM signatures and evaluate its validity in GRMHD simulations. △ Less

Submitted 4 July, 2024; originally announced July 2024.

AsterX: a new open-source GPU-accelerated GRMHD code for dynamical spacetimes

Authors: Jay V. Kalinani, Liwei Ji, Lorenzo Ennoggi, Federico G. Lopez Armengol, Lucas Timotheo Sanches, Bing-Jyun Tsao, Steven R Brandt, Manuela Campanelli, Riccardo Ciolfi, Bruno Giacomazzo, Roland Haas, Erik Schnetter, Yosef Zlochower

Abstract: We present AsterX, a novel open-source, GPU-accelerated, fully general relativistic magnetohydrodynamic (GRMHD) code designed for dynamic spacetimes in 3D Cartesian coordinates, and tailored for exascale computing. We utilize block-structured adaptive mesh refinement (AMR) through CarpetX, the new driver for the Einstein Toolkit, which is built on AMReX, a software framework for massively parallel… ▽ More We present AsterX, a novel open-source, GPU-accelerated, fully general relativistic magnetohydrodynamic (GRMHD) code designed for dynamic spacetimes in 3D Cartesian coordinates, and tailored for exascale computing. We utilize block-structured adaptive mesh refinement (AMR) through CarpetX, the new driver for the Einstein Toolkit, which is built on AMReX, a software framework for massively parallel applications. AsterX employs the Valencia formulation for GRMHD, coupled with the 'Z4c' formalism for spacetime evolution, while incorporating high resolution shock capturing schemes to accurately handle the hydrodynamics. AsterX has undergone rigorous testing in both static and dynamic spacetime, demonstrating remarkable accuracy and agreement with other codes in literature. Benchmarking the code through scaling tests on OLCF's Frontier supercomputer, we demonstrate a weak scaling efficiency of about 67%-77% on 4096 nodes compared to an 8-node performance. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 41 pages, 12 figures, to be submitted to Classical and Quantum Gravity

General Relativistic magneto-hydrodynamical simulations of accretion flows through traversable wormholes

Authors: Luciano Combi, Huan Yang, Eduardo Gutierrez, Scott C. Noble, Gustavo E. Romero, Manuela Campanelli

Abstract: We present the first dynamical model of plasma accretion onto traversable wormholes by performing General Relativistic magneto-hydrodynamical (GRMHD) simulations of the flow on both sides of the wormhole. We evolve the ideal MHD equations on a wormhole spacetime described by the spherically symmetric Simpson--Visser metric. The disk is initialized on one side of the wormhole and accretes onto the… ▽ More We present the first dynamical model of plasma accretion onto traversable wormholes by performing General Relativistic magneto-hydrodynamical (GRMHD) simulations of the flow on both sides of the wormhole. We evolve the ideal MHD equations on a wormhole spacetime described by the spherically symmetric Simpson--Visser metric. The disk is initialized on one side of the wormhole and accretes onto the throat driven by the magneto-rotational instability (MRI). We show that the inflowing plasma quickly settles in the throat and forms a hot, rotating cloud. The wormhole cloud acts as an engine in which gas coming from one side accumulates at the center, dissipates energy, and powers a mildly relativistic thermal wind toward the other side. Our novel predictions show that accreting wormholes behave very differently from black holes (BHs) in astrophysical environments. In particular, one mouth presents outflows without accretion signatures, contradicting the jet-disk symbiotic relation that holds for black holes. △ Less

Submitted 11 May, 2024; originally announced May 2024.

Comments: Accepted in PRD

Accretion onto a Supermassive Black Hole Binary Before Merger

Authors: Mark J. Avara, Julian H. Krolik, Manuela Campanelli, Scott C. Noble, Dennis Bowen, Taeho Ryu

Abstract: While supermassive binary black holes inspiral toward merger they may also experience significant accretion of matter from a surrounding disk. We study the dynamics of this system, simultaneously describing the evolving spacetime and magnetized plasma, and present the first relativistic calculation simulating two equal-mass, non-spinning black holes as they inspiral from a $20M$ ($G=c=1$) initial… ▽ More While supermassive binary black holes inspiral toward merger they may also experience significant accretion of matter from a surrounding disk. We study the dynamics of this system, simultaneously describing the evolving spacetime and magnetized plasma, and present the first relativistic calculation simulating two equal-mass, non-spinning black holes as they inspiral from a $20M$ ($G=c=1$) initial separation almost to merger, $\simeq 9M$ ($M$=binary mass). Our dynamical results imply important observational consequences: for instance, the accretion rate $\dot M$ onto the black holes first decreases and then reaches a plateau, dropping by only a factor of $\sim 3$ despite the rapid inspiral. An estimated bolometric light curve thus suggests some merging SMBBHs may be quite luminous past the predicted decoupling from the circumbinary disk. The minidisks through which the accretion reaches the black holes are very non-standard: Reynolds, not Maxwell, stresses dominate, and they oscillate between two states. In one part of the cycle, ``sloshing" streams transfer mass from one minidisk to the other through the L1 point at a rate $\sim 0.1\times$ the accretion rate, carrying kinetic energy at a rate that can be as large as the peak minidisk bolometric luminosity. We also discover that episodic accretion drives minidisks with time-varying tilts. The unsigned poloidal magnetic flux on the black hole event horizon is roughly constant at a dimensionless level $φ\sim 2-3$, but doubles just before merger; if the black holes had significant spin, this flux indicates the potential for powerful jets with variability driven by binary dynamics, another prediction of potentially unique EM signatures. This simulation is the first to employ our multipatch infrastructure \pwmhd, decreasing computational expense to $\sim 3\%$ of conventional single-grid methods' cost. △ Less

Submitted 16 January, 2024; v1 submitted 29 May, 2023; originally announced May 2023.

Comments: Comments welcome

Ameliorating the Courant-Friedrichs-Lewy condition in spherical coordinates: A double FFT filter method for general relativistic MHD in dynamical spacetimes

Authors: Liwei Ji, Vassilios Mewes, Yosef Zlochower, Lorenzo Ennoggi, Federico G. Lopez Armengol, Manuela Campanelli, Federico Cipolletta, Zachariah B. Etienne

Abstract: Numerical simulations of merging compact objects and their remnants form the theoretical foundation for gravitational wave and multi-messenger astronomy. While Cartesian-coordinate-based adaptive mesh refinement is commonly used for simulations, spherical-like coordinates are more suitable for nearly spherical remnants and azimuthal flows due to lower numerical dissipation in the evolution of flui… ▽ More Numerical simulations of merging compact objects and their remnants form the theoretical foundation for gravitational wave and multi-messenger astronomy. While Cartesian-coordinate-based adaptive mesh refinement is commonly used for simulations, spherical-like coordinates are more suitable for nearly spherical remnants and azimuthal flows due to lower numerical dissipation in the evolution of fluid angular momentum, as well as requiring fewer numbers of computational cells. However, the use of spherical coordinates to numerically solve hyperbolic partial differential equations can result in severe Courant-Friedrichs-Lewy (CFL) stability condition timestep limitations, which can make simulations prohibitively expensive. This paper addresses this issue for the numerical solution of coupled spacetime and general relativistic magnetohydrodynamics evolutions by introducing a double FFT filter and implementing it within the fully MPI-parallelized SphericalNR framework in the Einstein Toolkit. We demonstrate the effectiveness and robustness of the filtering algorithm by applying it to a number of challenging code tests, and show that it passes these tests effectively, demonstrating convergence while also increasing the timestep significantly compared to unfiltered simulations. △ Less

Submitted 9 November, 2023; v1 submitted 2 May, 2023; originally announced May 2023.

Comments: 19 pages, 15 figures, revtex4-1

Non-thermal radiation from dual jet interactions in supermassive black hole binaries

Authors: E. M. Gutiérrez, L. Combi, G. E. Romero, M. Campanelli

Abstract: Supermassive black hole binaries (SMBHBs) are natural by-products of galaxy mergers and are expected to be powerful multi-messenger sources. They can be powered by the accretion of matter and then radiate across the electromagnetic spectrum, much like normal active galactic nuclei (AGNs). Current electromagnetic observatories have a good chance of detecting and identifying these systems in the nea… ▽ More Supermassive black hole binaries (SMBHBs) are natural by-products of galaxy mergers and are expected to be powerful multi-messenger sources. They can be powered by the accretion of matter and then radiate across the electromagnetic spectrum, much like normal active galactic nuclei (AGNs). Current electromagnetic observatories have a good chance of detecting and identifying these systems in the near future. However, precise observational indicators are needed to distinguish individual AGNs from SMBHBs. In this paper, we propose a novel electromagnetic signature from SMBHBs: non-thermal emission produced by the interaction between the jets ejected by the black holes. We study close SMBHBs, which accrete matter from a circumbinary disc and the mini-discs formed around each hole. Each black hole ejects a magnetically dominated jet in the direction of its spin through the Blandford--Znajeck mechanism. We argue that in such a situation, the interaction between the jets can trigger strong magnetic reconnection events, where particles are accelerated and emit non-thermal radiation. Depending on whether the jets are aligned or misaligned, this radiation can have different periodicities. We model the evolution of the particles accelerated during the dual jet interaction and calculate their radiative output, obtaining spectra and providing estimates for the variability timescales. We finally discuss how this emission compares with that of normal AGNs. △ Less

Submitted 25 June, 2024; v1 submitted 10 January, 2023; originally announced January 2023.

Comments: 11 pages, 5 figures. Accepted to be published in MNRAS

Handing off the outcome of binary neutron star mergers for accurate and long-term post-merger simulations

Authors: Federico G. Lopez Armengol, Zachariah B. Etienne, Scott C. Noble, Bernard J. Kelly, Leonardo R. Werneck, Brendan Drachler, Manuela Campanelli, Federico Cipolletta, Yosef Zlochower, Ariadna Murguia-Berthier, Lorenzo Ennoggi, Mark Avara, Riccardo Ciolfi, Joshua Faber, Grace Fiacco, Bruno Giacomazzo, Tanmayee Gupte, Trung Ha, Julian H. Krolik, Vassilios Mewes, Richard O'Shaughnessy, Jesús M. Rueda-Becerril, Jeremy Schnittman

Abstract: We perform binary neutron star (BNS) merger simulations in full dynamical general relativity with IllinoisGRMHD, on a Cartesian grid with adaptive-mesh refinement. After the remnant black hole has become nearly stationary, the evolution of the surrounding accretion disk on Cartesian grids over long timescales (1s) is suboptimal, as Cartesian coordinates over-resolve the angular coordinates at larg… ▽ More We perform binary neutron star (BNS) merger simulations in full dynamical general relativity with IllinoisGRMHD, on a Cartesian grid with adaptive-mesh refinement. After the remnant black hole has become nearly stationary, the evolution of the surrounding accretion disk on Cartesian grids over long timescales (1s) is suboptimal, as Cartesian coordinates over-resolve the angular coordinates at large distances, and the accreting plasma flows obliquely across coordinate lines dissipating angular momentum artificially from the disk. To address this, we present the Handoff, a set of computational tools that enables the transfer of general relativistic magnetohydrodynamic (GRMHD) and spacetime data from IllinoisGRMHD to HARM3D, a GRMHD code that specializes in modeling black hole accretion disks in static spacetimes over long timescales, making use of general coordinate systems with spherical topology. We demonstrate that the Handoff allows for a smooth and reliable transition of GRMHD fields and spacetime data, enabling us to efficiently and reliably evolve BNS dynamics well beyond merger. We also discuss future plans, which involve incorporating advanced equations of state and neutrino physics into BNS simulations using the \handoff approach. △ Less

Submitted 31 October, 2022; v1 submitted 17 December, 2021; originally announced December 2021.

Comments: 27 pages, 22 figures

Journal ref: Phys. Rev. D 106, 083015 (2022)

Electromagnetic Signatures from Supermassive Binary Black Holes Approaching Merger

Authors: Eduardo M. Gutiérrez, Luciano Combi, Scott C. Noble, Manuela Campanelli, Julian H. Krolik, Federico G. López Armengol, Federico García

Abstract: We present fully relativistic predictions for the electromagnetic emission produced by accretion disks surrounding spinning and nonspinning supermassive binary black holes on the verge of merging. We use the code Bothros to post-process data from 3D general relativistic magnetohydrodynamic simulations via ray-tracing calculations. These simulations model the dynamics of a circumbinary disk and the… ▽ More We present fully relativistic predictions for the electromagnetic emission produced by accretion disks surrounding spinning and nonspinning supermassive binary black holes on the verge of merging. We use the code Bothros to post-process data from 3D general relativistic magnetohydrodynamic simulations via ray-tracing calculations. These simulations model the dynamics of a circumbinary disk and the mini-disks that form around two equal-mass black holes orbiting each other at an initial separation of 20 gravitational radii, and evolve the system for more than 10 orbits in the inspiral regime. We model the emission as the sum of thermal blackbody radiation emitted by an optically thick accretion disk and a power-law spectrum extending to hard X-rays emitted by a hot optically thin corona. We generate time-dependent spectra, images, and light curves at various frequencies to investigate intrinsic periodic signals in the emission, as well as the effects of the black hole spin. We find that prograde black hole spin makes mini-disks brighter since the smaller innermost stable circular orbit angular momentum demands more dissipation before matter plunges to the horizon. However, compared to mini-disks in larger separation binaries with spinning black holes, our mini-disks are less luminous: unlike those systems, their mass accretion rate is lower than in the circumbinary disk, and they radiate with lower efficiency because their inflow times are shorter. Compared to a single black hole system matched in mass and accretion rate, these binaries have spectra noticeably weaker and softer in the UV. Finally, we discuss the implications of our findings for the potential observability of these systems. △ Less

Submitted 15 April, 2022; v1 submitted 17 December, 2021; originally announced December 2021.

Comments: 19 pages, 15 figures; matching published version except for a clarification in Sec. 4.3

Journal ref: The Astrophysical Journal, Volume 928, Issue 2, id.137, 15 pp. (2022)

The Forward Physics Facility: Sites, Experiments, and Physics Potential

Authors: Luis A. Anchordoqui, Akitaka Ariga, Tomoko Ariga, Weidong Bai, Kincso Balazs, Brian Batell, Jamie Boyd, Joseph Bramante, Mario Campanelli, Adrian Carmona, Francesco G. Celiberto, Grigorios Chachamis, Matthew Citron, Giovanni De Lellis, Albert De Roeck, Hans Dembinski, Peter B. Denton, Antonia Di Crecsenzo, Milind V. Diwan, Liam Dougherty, Herbi K. Dreiner, Yong Du, Rikard Enberg, Yasaman Farzan, Jonathan L. Feng , et al. (56 additional authors not shown)

Abstract: The Forward Physics Facility (FPF) is a proposal to create a cavern with the space and infrastructure to support a suite of far-forward experiments at the Large Hadron Collider during the High Luminosity era. Located along the beam collision axis and shielded from the interaction point by at least 100 m of concrete and rock, the FPF will house experiments that will detect particles outside the acc… ▽ More The Forward Physics Facility (FPF) is a proposal to create a cavern with the space and infrastructure to support a suite of far-forward experiments at the Large Hadron Collider during the High Luminosity era. Located along the beam collision axis and shielded from the interaction point by at least 100 m of concrete and rock, the FPF will house experiments that will detect particles outside the acceptance of the existing large LHC experiments and will observe rare and exotic processes in an extremely low-background environment. In this work, we summarize the current status of plans for the FPF, including recent progress in civil engineering in identifying promising sites for the FPF and the experiments currently envisioned to realize the FPF's physics potential. We then review the many Standard Model and new physics topics that will be advanced by the FPF, including searches for long-lived particles, probes of dark matter and dark sectors, high-statistics studies of TeV neutrinos of all three flavors, aspects of perturbative and non-perturbative QCD, and high-energy astroparticle physics. △ Less

Submitted 25 May, 2022; v1 submitted 22 September, 2021; originally announced September 2021.

Comments: revised version, accepted by Physics Reports

Report number: BNL-222142-2021-FORE, CERN-PBC-Notes-2021-025, DESY-21-142, FERMILAB-CONF-21-452-AE-E-ND-PPD-T, KYUSHU-RCAPP-2021-01, LU TP 21-36, PITT-PACC-2118, SMU-HEP-21-10, UCI-TR-2021-22

Journal ref: Phys. Rept. 968 (2022), 1-50

Mini-disk accretion onto spinning black hole binaries: quasi-periodicities and outflows

Authors: Luciano Combi, Federico G. Lopez Armengol, Manuela Campanelli, Scott C. Noble, Mark Avara, Julian H. Krolik, Dennis Bowen

Abstract: We perform a full 3D general relativistic magnetohydrodynamical (GRMHD) simulation of an equal-mass, spinning, binary black hole approaching merger, surrounded by a circumbinary disk and with mini-disks around each black hole. For this purpose, we evolve the ideal GRMHD equations on top of an approximated spacetime for the binary that is valid in every position of space, including the black hole h… ▽ More We perform a full 3D general relativistic magnetohydrodynamical (GRMHD) simulation of an equal-mass, spinning, binary black hole approaching merger, surrounded by a circumbinary disk and with mini-disks around each black hole. For this purpose, we evolve the ideal GRMHD equations on top of an approximated spacetime for the binary that is valid in every position of space, including the black hole horizons, during the inspiral regime. We use relaxed initial data for the circumbinary disk from a previous long-term simulation, where the accretion is dominated by an $m=1$ overdensity called the lump. We compare our new spinning simulation with a previous non-spinning run, studying how spin influences the mini-disk properties. We analyze the accretion from the inner edge of the lump to the black hole, focusing on the angular momentum budget of the fluid around the mini-disks. We find that mini-disks in the spinning case have more mass over a cycle than the non-spinning case. However, in both cases, we find most of the mass received by the black holes is delivered by the direct plunging of material from the lump. We also analyze the morphology and variability of the electromagnetic fluxes and we find they share the same periodicities of the accretion rate. In the spinning case, we find that the outflows are $8$ times stronger than the non-spinning case. Our results will be useful to understand and produce realistic synthetic light curves and spectra, which can be used in future observations. △ Less

Submitted 8 February, 2022; v1 submitted 3 September, 2021; originally announced September 2021.

Comments: Accepted for publication in ApJ

Searching for solar KDAR with DUNE

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, D. Adams, M. Adinolfi, A. Aduszkiewicz, J. Aguilar, Z. Ahmad, J. Ahmed, B. Ali-Mohammadzadeh, T. Alion, K. Allison, S. Alonso Monsalve, M. Alrashed, C. Alt, A. Alton, P. Amedo, J. Anderson, C. Andreopoulos, M. Andreotti, M. P. Andrews , et al. (1157 additional authors not shown)

Abstract: The observation of 236 MeV muon neutrinos from kaon-decay-at-rest (KDAR) originating in the core of the Sun would provide a unique signature of dark matter annihilation. Since excellent angle and energy reconstruction are necessary to detect this monoenergetic, directional neutrino flux, DUNE with its vast volume and reconstruction capabilities, is a promising candidate for a KDAR neutrino search.… ▽ More The observation of 236 MeV muon neutrinos from kaon-decay-at-rest (KDAR) originating in the core of the Sun would provide a unique signature of dark matter annihilation. Since excellent angle and energy reconstruction are necessary to detect this monoenergetic, directional neutrino flux, DUNE with its vast volume and reconstruction capabilities, is a promising candidate for a KDAR neutrino search. In this work, we evaluate the proposed KDAR neutrino search strategies by realistically modeling both neutrino-nucleus interactions and the response of DUNE. We find that, although reconstruction of the neutrino energy and direction is difficult with current techniques in the relevant energy range, the superb energy resolution, angular resolution, and particle identification offered by DUNE can still permit great signal/background discrimination. Moreover, there are non-standard scenarios in which searches at DUNE for KDAR in the Sun can probe dark matter interactions. △ Less

Submitted 26 October, 2021; v1 submitted 19 July, 2021; originally announced July 2021.

Comments: 19 pages, 13 figures

Report number: FERMILAB-PUB-21-322-LBNF-ND

Journal ref: JCAP10(2021)065

HARM3D+NUC: A new method for simulating the post-merger phase of binary neutron star mergers with GRMHD, tabulated EOS and neutrino leakage

Authors: Ariadna Murguia-Berthier, Scott C. Noble, Luke F. Roberts, Enrico Ramirez-Ruiz, Leonardo R. Werneck, Michael Kolacki, Zachariah B. Etienne, Mark Avara, Manuela Campanelli, Riccardo Ciolfi, Federico Cipolletta, Brendan Drachler, Lorenzo Ennoggi, Joshua Faber, Grace Fiacco, Bruno Giacomazzo, Tanmayee Gupte, Trung Ha, Bernard J. Kelly, Julian H. Krolik, Federico G. Lopez Armengol, Ben Margalit, Tim Moon, Richard O'Shaughnessy, Jesús M. Rueda-Becerril , et al. (3 additional authors not shown)

Abstract: The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merge… ▽ More The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind. △ Less

Submitted 31 January, 2022; v1 submitted 9 June, 2021; originally announced June 2021.

Comments: 19 pages, 14 figures, published in ApJ

Mass-ratio and Magnetic Flux-Dependence of Modulated Accretion from Circumbinary Disks

Authors: Scott C. Noble, Julian H. Krolik, Manuela Campanelli, Yosef Zlochower, Bruno C. Mundim, Hiroyuki Nakano, Miguel Zilhão

Abstract: Accreting supermassive binary black holes (SMBBHs) are potential multi-messenger sources because they emit both gravitational wave and electromagnetic (EM) radiation. Past work has shown that their EM output may be periodically modulated by an asymmetric density distribution in the circumbinary disk, often called an "overdensity" or "lump;" this modulation could possibly be used to identify a sour… ▽ More Accreting supermassive binary black holes (SMBBHs) are potential multi-messenger sources because they emit both gravitational wave and electromagnetic (EM) radiation. Past work has shown that their EM output may be periodically modulated by an asymmetric density distribution in the circumbinary disk, often called an "overdensity" or "lump;" this modulation could possibly be used to identify a source as a binary. We explore the sensitivity of the overdensity to SMBBH mass ratio and magnetic flux through the accretion disk. We find that the relative amplitude of the overdensity and its associated EM periodic signal both degrade with diminishing mass ratio, vanishing altogether somewhere between 1:2 and 1:5. Greater magnetization also weakens the lump and any modulation of the light output. We develop a model to describe how lump formation results from internal stress degrading faster in the lump region than it can be rejuvenated through accretion inflow, and predicts a threshold value in specific internal stress below which lump formation should occur and which all our lump-forming simulations satisfy. Thus, detection of such a modulation would provide a constraint on both mass-ratio and magnetic flux piercing the accretion flow. △ Less

Submitted 22 March, 2021; originally announced March 2021.

Comments: 32 pages, 30 color figures. Submitted to the Astrophysical Journal. Movies may be found at https://youtube.com/playlist?list=PLNaEA0qwDBaeApzLr2oarVKiO3AFnklTV

Circumbinary Disk Accretion into Spinning Black Hole Binaries

Authors: Federico G. Lopez Armengol, Luciano Combi, Manuela Campanelli, Scott C. Noble, Julian H. Krolik, Dennis B. Bowen, Mark J. Avara, Vassilios Mewes, Hiroyuki Nakano

Abstract: Supermassive black hole binaries are likely to accrete interstellar gas through a circumbinary disk. Shortly before merger, the inner portions of this circumbinary disk are subject to general relativistic effects. To study this regime, we approximate the spacetime metric of close orbiting black holes by superimposing two boosted Kerr-Schild terms. After demonstrating the quality of this approximat… ▽ More Supermassive black hole binaries are likely to accrete interstellar gas through a circumbinary disk. Shortly before merger, the inner portions of this circumbinary disk are subject to general relativistic effects. To study this regime, we approximate the spacetime metric of close orbiting black holes by superimposing two boosted Kerr-Schild terms. After demonstrating the quality of this approximation, we carry out very long-term general relativistic magnetohydrodynamic simulations of the circumbinary disk. We consider black holes with spin dimensionless parameters of magnitude 0.9, in one simulation parallel to the orbital angular momentum of the binary, but in another anti-parallel. These are contrasted with spinless simulations. We find that, for a fixed surface mass density in the inner circumbinary disk, aligned spins of this magnitude approximately reduce the mass accretion rate by 14% and counter-aligned spins increase it by 45%, leaving many other disk properties unchanged. △ Less

Submitted 11 May, 2021; v1 submitted 30 January, 2021; originally announced February 2021.

Comments: 23 pages, 15 figures, Submitted to ApJ

Hubble Constant Measurement with GW190521 as an Eccentric Black Hole Merger

Authors: V. Gayathri, J. Healy, J. Lange, B. O'Brien, M. Szczepanczyk, I. Bartos, M. Campanelli, S. Klimenko, C. Lousto, R. O'Shaughnessy

Abstract: Gravitational wave observations can be used to accurately measure the Hubble constant $H_0$ and could help understand the present discrepancy between constraints from Type Ia supernovae and the cosmic microwave background. Neutron star mergers are primarily used for this purpose as their electromagnetic emission can be used to greatly reduce measurement uncertainties. Here we estimate $H_0$ using… ▽ More Gravitational wave observations can be used to accurately measure the Hubble constant $H_0$ and could help understand the present discrepancy between constraints from Type Ia supernovae and the cosmic microwave background. Neutron star mergers are primarily used for this purpose as their electromagnetic emission can be used to greatly reduce measurement uncertainties. Here we estimate $H_0$ using the recently observed black hole merger GW190521 and its candidate electromagnetic counterpart found by ZTF using a highly eccentric explanation of the properties of GW190521. We find that the reconstructed distance of GW190521 and the redshift of the candidate host galaxy are more consistent with standard cosmology for our eccentric model than if we reconstruct the source parameters assuming no eccentricity. We obtain $H_0=88.6^{+17.1}_{-34.3}$\,km\,s$^{-1}$Mpc$^{-1}$ for GW190521, and $H_0=73.4^{+6.9}_{-10.7}$\,km\,s$^{-1}$Mpc$^{-1}$ in combination with the results of the neutron star merger GW170817. Our results indicate that future $H_0$ computations using black hole mergers will need to account for possible eccentricity. For extreme cases, the orbital velocity of binaries in AGN disks can represent a significant systematic uncertainty. △ Less

Submitted 29 September, 2020; originally announced September 2020.

Comments: 5 pages, 3 figures

Eccentricity Estimate for Black Hole Mergers with Numerical Relativity Simulations

Authors: V. Gayathri, J. Healy, J. Lange, B. O'Brien, M. Szczepanczyk, I. Bartos, M. Campanelli, S. Klimenko, C. Lousto, R. O'Shaughnessy

Abstract: The origin of black hole mergers discovered by the LIGO and Virgo gravitational-wave observatories is currently unknown. GW190521 is the heaviest black hole merger detected so far. Its observed high mass and possible spin-induced orbital precession could arise from the binary having formed following a close encounter. An observational signature of close encounters is eccentric binary orbit; howeve… ▽ More The origin of black hole mergers discovered by the LIGO and Virgo gravitational-wave observatories is currently unknown. GW190521 is the heaviest black hole merger detected so far. Its observed high mass and possible spin-induced orbital precession could arise from the binary having formed following a close encounter. An observational signature of close encounters is eccentric binary orbit; however, this feature is currently difficult to identify due to the lack of suitable gravitational waveforms. No eccentric merger has been previously found. Here we report 611 numerical relativity simulations covering the full eccentricity range and an estimation approach to probe the eccentricity of mergers. Our set of simulations corresponds to $\sim 10^5$ waveforms, comparable to the number used in gravitational wave searches, albeit with coarser mass-ratio and spin resolution. We applied our approach to GW190521 and found that it is the most consistent with a highly eccentric ($e=0.69^{+0.17}_{-0.22}$; 90% credible level) merger within our set of waveforms. This interpretation is supported over a non-eccentric merger with $>10$ Odds ratio if $\gtrsim10\%$ of GW190521-like mergers are highly eccentric. Detectable orbital eccentricity would be evidence against an isolated binary origin, which is otherwise difficult to rule out based on observed mass and spin. △ Less

Submitted 27 January, 2022; v1 submitted 11 September, 2020; originally announced September 2020.

Comments: Published in Nature Astronomy

Supernova Neutrino Burst Detection with the Deep Underground Neutrino Experiment

Authors: DUNE collaboration, B. Abi, R. Acciarri, M. A. Acero, G. Adamov, D. Adams, M. Adinolfi, Z. Ahmad, J. Ahmed, T. Alion, S. Alonso Monsalve, C. Alt, J. Anderson, C. Andreopoulos, M. P. Andrews, F. Andrianala, S. Andringa, A. Ankowski, M. Antonova, S. Antusch, A. Aranda-Fernandez, A. Ariga, L. O. Arnold, M. A. Arroyave, J. Asaadi , et al. (949 additional authors not shown)

Abstract: The Deep Underground Neutrino Experiment (DUNE), a 40-kton underground liquid argon time projection chamber experiment, will be sensitive to the electron-neutrino flavor component of the burst of neutrinos expected from the next Galactic core-collapse supernova. Such an observation will bring unique insight into the astrophysics of core collapse as well as into the properties of neutrinos. The gen… ▽ More The Deep Underground Neutrino Experiment (DUNE), a 40-kton underground liquid argon time projection chamber experiment, will be sensitive to the electron-neutrino flavor component of the burst of neutrinos expected from the next Galactic core-collapse supernova. Such an observation will bring unique insight into the astrophysics of core collapse as well as into the properties of neutrinos. The general capabilities of DUNE for neutrino detection in the relevant few- to few-tens-of-MeV neutrino energy range will be described. As an example, DUNE's ability to constrain the $ν_e$ spectral parameters of the neutrino burst will be considered. △ Less

Submitted 29 May, 2021; v1 submitted 15 August, 2020; originally announced August 2020.

Comments: 29 pages, 17 figures; paper based on DUNE Technical Design Report. arXiv admin note: substantial text overlap with arXiv:2002.03005

Report number: FERMILAB-PUB-20-380-LBNF

AstroDance: Engaging Deaf and Hard-of-Hearing Students in Astrophysics via Multimedia Performances

Authors: J. Nordhaus, M. Campanelli, J. Bochner, T. Warfield, H. -P. Bischof, J. Noel-Storr

Abstract: The dynamics of gravitating astrophysical systems such as black holes and neutron stars are fascinatingly complex, offer some of nature's most spectacular phenomena, and capture the public's imagination in ways that few subjects can. Here, we describe {\it AstroDance}, a multi-media project to engage deaf and hard-of-hearing (DHH) students in astronomy and gravitational physics. {\it AstroDance} i… ▽ More The dynamics of gravitating astrophysical systems such as black holes and neutron stars are fascinatingly complex, offer some of nature's most spectacular phenomena, and capture the public's imagination in ways that few subjects can. Here, we describe {\it AstroDance}, a multi-media project to engage deaf and hard-of-hearing (DHH) students in astronomy and gravitational physics. {\it AstroDance} incorporates multiple means of representation of scientific concepts and was performed primarily for secondary and post-secondary audiences at $\sim$20 venues in the northeastern US prior to the historic first detection of gravitational waves. As part of the {\it AstroDance} project, we surveyed $\sim$1000 audience members roughly split evenly between hearing and DHH audience members. While both groups reported statistically equivalent high-rates of enjoyment of the performance, the DHH group reported an increase in how much they learned about science at a statistically significant rate compared to the hearing audience. Our findings suggest that multi-sensory approaches benefit both hearing and deaf audiences and enable accessible participation for broader groups. △ Less

Submitted 18 June, 2020; originally announced June 2020.

Comments: Accepted to Journal of Science Education for Students with Disabilities

Numerical relativity in spherical coordinates: A new dynamical spacetime and general relativistic MHD evolution framework for the Einstein Toolkit

Authors: Vassilios Mewes, Yosef Zlochower, Manuela Campanelli, Thomas W. Baumgarte, Zachariah B. Etienne, Federico G. Lopez Armengol, Federico Cipolletta

Abstract: We present SphericalNR, a new framework for the publicly available Einstein Toolkit that numerically solves the Einstein field equations coupled to the equations of general relativistic magnetohydrodynamics (GRMHD) in a 3+1 split of spacetime in spherical coordinates without symmetry assumptions. The spacetime evolution is performed using reference-metric versions of either the Baumgarte-Shapiro-S… ▽ More We present SphericalNR, a new framework for the publicly available Einstein Toolkit that numerically solves the Einstein field equations coupled to the equations of general relativistic magnetohydrodynamics (GRMHD) in a 3+1 split of spacetime in spherical coordinates without symmetry assumptions. The spacetime evolution is performed using reference-metric versions of either the Baumgarte-Shapiro-Shibata-Nakamura equations or the fully covariant and conformal Z4 system with constraint damping. We have developed a reference-metric version of the Valencia formulation of GRMHD with a vector potential method, guaranteeing the absence of magnetic monopoles during the evolution. In our framework, every dynamical field (both spacetime and matter) is evolved using its components in an orthonormal basis with respect to the spherical reference-metric. Furthermore, all geometric information about the spherical coordinate system is encoded in source terms appearing in the evolution equations. This allows for the straightforward extension of Cartesian high-resolution shock-capturing finite volume codes to use spherical coordinates with our framework. To this end, we have adapted GRHydro, a Cartesian finite volume GRMHD code already available in the Einstein Toolkit, to use spherical coordinates. We present the full evolution equations of the framework, as well as details of its implementation in the Einstein Toolkit. We validate SphericalNR by demonstrating it passes a variety of challenging code tests in static and dynamical spacetimes. △ Less

Submitted 5 May, 2020; v1 submitted 14 February, 2020; originally announced February 2020.

Comments: 29 pages, 13 figures, published in Physical Review D

Journal ref: Phys. Rev. D 101, 104007 (2020)

Quasi-Periodicity of Supermassive Binary Black Hole Accretion Approaching Merger

Authors: Dennis B. Bowen, Vassilios Mewes, Scott C. Noble, Mark Avara, Manuela Campanelli, Julian H. Krolik

Abstract: In this paper we continue the first ever study of magnetized mini-disks coupled to circumbinary accretion in a supermassive binary black hole (SMBBH) approaching merger reported in Bowen et al. 2018. We extend this simulation from 3 to 12 binary orbital periods. We find that relativistic SMBBH accretion acts as a resonant cavity, where quasi-periodic oscillations tied to the the frequency at which… ▽ More In this paper we continue the first ever study of magnetized mini-disks coupled to circumbinary accretion in a supermassive binary black hole (SMBBH) approaching merger reported in Bowen et al. 2018. We extend this simulation from 3 to 12 binary orbital periods. We find that relativistic SMBBH accretion acts as a resonant cavity, where quasi-periodic oscillations tied to the the frequency at which the black hole's orbital phase matches a non-linear $m=1$ density feature, or ``lump'', in the circumbinary accretion disk permeate the system. The rate of mass accretion onto each of the mini-disks around the black holes is modulated at the beat frequency between the binary frequency and the lump's mean orbital frequency, i.e., $Ω_{\rm beat} = Ω_{\rm bin} - \barΩ_{\rm lump}$, while the total mass accretion rate of this equal-mass binary is modulated at two different frequencies, $\gtrsim \barΩ_{\rm lump}$ and $\approx 2 Ω_{\rm beat}$. The instantaneous rotation rate of the lump itself is also modulated at two frequencies close to the modulation frequencies of the total accretion rate, $\barΩ_{\rm lump}$ and $2 Ω_{\rm beat}$. Because of the compact nature of the mini-disks in SMBBHs approaching merger, the inflow times within the mini-disks are comparable to the period on which their mass-supply varies, so that their masses---and the accretion rates they supply to their black holes---are strongly modulated at the same frequency. In essence, the azimuthal symmetry of the circumbinary disk is broken by the dynamics of orbits near a binary, and this $m=1$ asymmetry then drives quasi-periodic variation throughout the system, including both accretion and disk-feeding. In SMBBHs approaching merger, such time variability could introduce distinctive, increasingly rapid, fluctuations in their electromagnetic emission. △ Less

Submitted 26 April, 2019; originally announced April 2019.

Comments: 10 pages, 7 figures

On the properties of the massive binary black hole merger GW170729

Authors: Katerina Chatziioannou, Roberto Cotesta, Sudarshan Ghonge, Jacob Lange, Ken K. -Y. Ng, Juan Calderon Bustillo, James Clark, Carl-Johan Haster, Sebastian Khan, Michael Puerrer, Vivien Raymond, Salvatore Vitale, Nousha Afshari, Stanislav Babak, Kevin Barkett, Jonathan Blackman, Alejandro Bohe, Michael Boyle, Alessandra Buonanno, Manuela Campanelli, Gregorio Carullo, Tony Chu, Eric Flynn, Heather Fong, Alyssa Garcia , et al. (31 additional authors not shown)

Abstract: We present a detailed investigation into the properties of GW170729, the gravitational wave with the most massive and distant source confirmed to date. We employ an extensive set of waveform models, including new improved models that incorporate the effect of higher-order waveform modes which are particularly important for massive systems. We find no indication of spin-precession, but the inclusio… ▽ More We present a detailed investigation into the properties of GW170729, the gravitational wave with the most massive and distant source confirmed to date. We employ an extensive set of waveform models, including new improved models that incorporate the effect of higher-order waveform modes which are particularly important for massive systems. We find no indication of spin-precession, but the inclusion of higher-order modes in the models results in an improved estimate for the mass ratio of $(0.3-0.8)$ at the 90\% credible level. Our updated measurement excludes equal masses at that level. We also find that models with higher-order modes lead to the data being more consistent with a smaller effective spin, with the probability that the effective spin is greater than zero being reduced from $99\%$ to $94\%$. The 90\% credible interval for the effective spin parameter is now $(-0.01-0.50)$. Additionally, the recovered signal-to-noise ratio increases by $\sim0.3$ units compared to analyses without higher-order modes. We study the effect of common spin priors on the derived spin and mass measurements, and observe small shifts in the spins, while the masses remain unaffected. We argue that our conclusions are robust against systematic errors in the waveform models. We also compare the above waveform-based analysis which employs compact-binary waveform models to a more flexible wavelet- and chirplet-based analysis. We find consistency between the two, with overlaps of $\sim 0.9$, typical of what is expected from simulations of signals similar to GW170729, confirming that the data are well-described by the existing waveform models. Finally, we study the possibility that the primary component of GW170729 was the remnant of a past merger of two black holes and find this scenario to be indistinguishable from the standard formation scenario. △ Less

Submitted 7 November, 2019; v1 submitted 15 March, 2019; originally announced March 2019.

Comments: 14 pages, 8 figures, final published version, samples available at https://git.ligo.org/katerina.chatziioannou/gw170729hm_datarelease

Journal ref: Phys. Rev. D 100, 104015 (2019)

The second RIT binary black hole simulations catalog and its application to gravitational waves parameter estimation

Authors: James Healy, Carlos O. Lousto, Jacob Lange, Richard O'Shaughnessy, Yosef Zlochower, Manuela Campanelli

Abstract: The RIT numerical relativity group is releasing the second public catalog of black-hole-binary waveforms \url{http://ccrg.rit.edu/~RITCatalog}. This release consists of 320 accurate simulations that include 46 precessing and 274 nonprecessing binary systems with mass ratios $q=m_1/m_2$ in the range $1/6\leq q\leq1$ and individual spins up to $s/m^2=0.95$. The new catalog contains search and orderi… ▽ More The RIT numerical relativity group is releasing the second public catalog of black-hole-binary waveforms \url{http://ccrg.rit.edu/~RITCatalog}. This release consists of 320 accurate simulations that include 46 precessing and 274 nonprecessing binary systems with mass ratios $q=m_1/m_2$ in the range $1/6\leq q\leq1$ and individual spins up to $s/m^2=0.95$. The new catalog contains search and ordering tools for the waveforms based on initial parameters of the binary, trajectory information, peak radiation, and final remnant black hole properties. The final black hole remnant properties provided here can be used to model the merger of black-hole binaries from its initial configurations. The waveforms are extrapolated to infinite observer location and can be used to independently interpret gravitational wave signals from laser interferometric detectors. As an application of this waveform catalog we reanalyze the signal of GW150914 implementing parameter estimation techniques that make use of only numerical waveforms without any reference to information from phenomenological models. △ Less

Submitted 27 March, 2019; v1 submitted 8 January, 2019; originally announced January 2019.

Comments: 23 pages, 10 figures. Major extension of the GW150914 analysis

Journal ref: Phys. Rev. D 100, 024021 (2019)

Electromagnetic Emission from Supermassive Binary Black Holes Approaching Merger

Authors: Stéphane d'Ascoli, Scott C. Noble, Dennis B. Bowen, Manuela Campanelli, Julian H. Krolik, Vassilios Mewes

Abstract: We present the first fully relativistic prediction of the electromagnetic emission from the surrounding gas of a supermassive binary black hole system approaching merger. Using a ray-tracing code to post-process data from a general relativistic 3-d MHD simulation, we generate images and spectra, and analyze the viewing angle dependence of the light emitted. When the accretion rate is relatively hi… ▽ More We present the first fully relativistic prediction of the electromagnetic emission from the surrounding gas of a supermassive binary black hole system approaching merger. Using a ray-tracing code to post-process data from a general relativistic 3-d MHD simulation, we generate images and spectra, and analyze the viewing angle dependence of the light emitted. When the accretion rate is relatively high, the circumbinary disk, accretion streams, and mini-disks combine to emit light in the UV/EUV bands. We posit a thermal Compton hard X-ray spectrum for coronal emission; at high accretion rates, it is almost entirely produced in the mini-disks, but at lower accretion rates it is the primary radiation mechanism in the mini-disks and accretion streams as well. Due to relativistic beaming and gravitational lensing, the angular distribution of the power radiated is strongly anisotropic, especially near the equatorial plane. △ Less

Submitted 14 June, 2018; originally announced June 2018.

Comments: Submitted to the Astrophysical Journal, 11 color figures

Numerical relativity in spherical coordinates with the Einstein Toolkit

Authors: Vassilios Mewes, Yosef Zlochower, Manuela Campanelli, Ian Ruchlin, Zachariah B. Etienne, Thomas W. Baumgarte

Abstract: Numerical relativity codes that do not make assumptions on spatial symmetries most commonly adopt Cartesian coordinates. While these coordinates have many attractive features, spherical coordinates are much better suited to take advantage of approximate symmetries in a number of astrophysical objects, including single stars, black holes and accretion disks. While the appearance of coordinate singu… ▽ More Numerical relativity codes that do not make assumptions on spatial symmetries most commonly adopt Cartesian coordinates. While these coordinates have many attractive features, spherical coordinates are much better suited to take advantage of approximate symmetries in a number of astrophysical objects, including single stars, black holes and accretion disks. While the appearance of coordinate singularities often spoils numerical relativity simulations in spherical coordinates, especially in the absence of any symmetry assumptions, it has recently been demonstrated that these problems can be avoided if the coordinate singularities are handled analytically. This is possible with the help of a reference-metric version of the Baumgarte-Shapiro-Shibata-Nakamura formulation together with a proper rescaling of tensorial quantities. In this paper we report on an implementation of this formalism in the Einstein Toolkit. We adapt the Einstein Toolkit infrastructure, originally designed for Cartesian coordinates, to handle spherical coordinates, by providing appropriate boundary conditions at both inner and outer boundaries. We perform numerical simulations for a disturbed Kerr black hole, extract the gravitational wave signal, and demonstrate that the noise in these signals is orders of magnitude smaller when computed on spherical grids rather than Cartesian grids. With the public release of our new Einstein Toolkit thorns, our methods for numerical relativity in spherical coordinates will become available to the entire numerical relativity community. △ Less

Submitted 2 May, 2018; v1 submitted 26 February, 2018; originally announced February 2018.

Comments: 14 pages, 8 figures, published in Phys. Rev. D

Journal ref: Phys. Rev. D 97, 084059 (2018)

Targeted numerical simulations of binary black holes for GW170104

Authors: James Healy, Jacob Lange, Richard O'Shaughnessy, Carlos Lousto, Manuela Campanelli, Andrew Williamson, Yosef Zlochower, Juan Calderon Bustillo, James Clark, Christopher Evans, D. Ferguson, Sudarshan Ghonge, Karan Jani, Bhavesh Khamesra, Pablo Laguna, Deirdre M. Shoemaker, Alyssa Garcia, Michael Boyle, Daniel Hemberger, Lawrence Kidder, Prayush Kumar, Geoffrey Lovelace, Harald Pfeiffer, Mark Scheel, Saul Teukolsky

Abstract: In response to LIGO's observation of GW170104, we performed a series of full numerical simulations of binary black holes, each designed to replicate likely realizations of its dynamics and radiation. These simulations have been performed at multiple resolutions and with two independent techniques to solve Einstein's equations. For the nonprecessing and precessing simulations, we demonstrate the tw… ▽ More In response to LIGO's observation of GW170104, we performed a series of full numerical simulations of binary black holes, each designed to replicate likely realizations of its dynamics and radiation. These simulations have been performed at multiple resolutions and with two independent techniques to solve Einstein's equations. For the nonprecessing and precessing simulations, we demonstrate the two techniques agree mode by mode, at a precision substantially in excess of statistical uncertainties in current LIGO's observations. Conversely, we demonstrate our full numerical solutions contain information which is not accurately captured with the approximate phenomenological models commonly used to infer compact binary parameters. To quantify the impact of these differences on parameter inference for GW170104 specifically, we compare the predictions of our simulations and these approximate models to LIGO's observations of GW170104. △ Less

Submitted 19 December, 2017; v1 submitted 15 December, 2017; originally announced December 2017.

Comments: 11 figures, 20 pages

Journal ref: Phys. Rev. D 97, 064027 (2018)

Quasi-Periodic Behavior of Mini-Disks in Binary Black Holes Approaching Merger

Authors: Dennis B. Bowen, Vassilios Mewes, Manuela Campanelli, Scott C. Noble, Julian H. Krolik, Miguel Zilhao

Abstract: We present the first magnetohydrodynamic simulation in which a circumbinary disk around a relativistic binary black hole feeds mass to individual accretion disks ("mini-disks") around each black hole. Mass flow through the accretion streams linking the circumbinary disk to the mini-disks is modulated quasi-periodically by the streams' interaction with a nonlinear $m=1$ density feature, or "lump",… ▽ More We present the first magnetohydrodynamic simulation in which a circumbinary disk around a relativistic binary black hole feeds mass to individual accretion disks ("mini-disks") around each black hole. Mass flow through the accretion streams linking the circumbinary disk to the mini-disks is modulated quasi-periodically by the streams' interaction with a nonlinear $m=1$ density feature, or "lump", at the inner edge of the circumbinary disk: the stream supplying each mini-disk comes into phase with the lump at a frequency $0.74$ times the binary orbital frequency. Because the binary is relativistic, the tidal truncation radii of the mini-disks are not much larger than their innermost stable circular orbits; consequently, the mini-disks' inflow times are shorter than the conventional estimate and are comparable to the stream modulation period. As a result, the mini-disks are always in inflow disequilibrium, with their masses and spiral density wave structures responding to the stream's quasi-periodic modulation. The fluctuations in each mini-disk's mass are so large that as much as $75\%$ of the total mini-disk mass can be contained within a single mini-disk. Such quasi-periodic modulation of the mini-disk structure may introduce distinctive time-dependent features in the binary's electromagnetic emission. △ Less

Submitted 24 January, 2018; v1 submitted 14 December, 2017; originally announced December 2017.

Comments: 6 pages, 4 figures, Accepted for publication in ApJL

Modeling the Black hole Merger of QSO 3C 186

Authors: Carlos O. Lousto, Yosef Zlochower, Manuela Campanelli

Abstract: Recent detailed observations of the radio-loud quasar 3C 186 indicate the possibility that a supermassive recoiling black hole is moving away from the host galaxy at a speed of nearly 2100km/s. If this is the case, we can model the mass ratio and spins of the progenitor binary black hole using the results of numerical relativity simulations. We find that the black holes in the progenitor must have… ▽ More Recent detailed observations of the radio-loud quasar 3C 186 indicate the possibility that a supermassive recoiling black hole is moving away from the host galaxy at a speed of nearly 2100km/s. If this is the case, we can model the mass ratio and spins of the progenitor binary black hole using the results of numerical relativity simulations. We find that the black holes in the progenitor must have comparable masses with a mass ratio $q=m_1/m_2>1/4$ and the spin of the primary black hole must be $α_2=S_2/m_2^2>0.4$. The final remnant of the merger is bounded by $α_f>0.45$ and at least $4\%$ of the total mass of the binary system is radiated into gravitational waves. We consider four different pre-merger scenarios that further narrow those values. Assuming, for instance, a cold accretion driven merger model, we find that the binary had comparable masses with $q=0.70^{+0.29}_{-0.21}$ and the normalized spins of the larger and smaller black holes were $α_2=0.94^{+0.06}_{-0.22}$ and $α_1=0.95^{+0.05}_{-0.09}$. We can also estimate the final recoiling black hole spin $α_f=0.93^{+0.02}_{-0.03}$ and that the system radiated $9.6^{+0.8}_{-1.4}\%$ of its total mass, making the merger of those black holes the most energetic event ever observed. △ Less

Submitted 3 April, 2017; originally announced April 2017.

Comments: 4 pages, 2 figures

The RIT binary black hole simulations catalog

Authors: James Healy, Carlos O. Lousto, Yosef Zlochower, Manuela Campanelli

Abstract: The RIT numerical relativity group is releasing a public catalog of black-hole-binary waveforms. The initial release of the catalog consists of 126 recent simulations that include precessing and non precessing systems with mass ratios $q=m_1/m_2$ in the range $1/6\leq q\leq1$. The catalog contains information about the initial data of the simulation, the waveforms extrapolated to infinity, as well… ▽ More The RIT numerical relativity group is releasing a public catalog of black-hole-binary waveforms. The initial release of the catalog consists of 126 recent simulations that include precessing and non precessing systems with mass ratios $q=m_1/m_2$ in the range $1/6\leq q\leq1$. The catalog contains information about the initial data of the simulation, the waveforms extrapolated to infinity, as well as information about the peak luminosity and final remnant black hole properties. These waveforms can be used to independently interpret gravitational wave signals from laser interferometric detectors and △ Less

Submitted 9 March, 2017; originally announced March 2017.

Comments: 5 pages, 3 figures

Relativistic Dynamics and Mass Exchange in Binary Black Hole Mini-Disks

Authors: Dennis B. Bowen, Manuela Campanelli, Julian H. Krolik, Vassilios Mewes, Scott C. Noble

Abstract: We present the first exploration of relativistic gas dynamics in the immediate vicinity of binary black holes as the system inspirals close to merger in the gravitational radiation-driven regime. We focus on 2D hydrodynamical studies of comparable-mass, non-spinning systems. Relativistic effects alter the dynamics of gas in this environment in several ways. Because the gravitational potential betw… ▽ More We present the first exploration of relativistic gas dynamics in the immediate vicinity of binary black holes as the system inspirals close to merger in the gravitational radiation-driven regime. We focus on 2D hydrodynamical studies of comparable-mass, non-spinning systems. Relativistic effects alter the dynamics of gas in this environment in several ways. Because the gravitational potential between the two black holes becomes shallower than in the Newtonian regime, the mini-disks stretch toward the L1 point and the amount of gas passing back and forth between the mini-disks increases sharply with decreasing binary separation. This "sloshing" is quasi-periodically modulated at $2$ and $2.75$ times the binary orbital frequency, corresponding to timescales of hours to days for supermassive binary black holes. In addition, relativistic effects add an $m=1$ component to the tidally-driven spiral waves in the disks that are purely $m=2$ in Newtonian gravity; this component becomes dominant when the separation is $\lesssim 100$ gravitational radii. Both the sloshing and the spiral waves have the potential to create distinctive radiation features that may uniquely mark supermassive binary black holes in the relativistic regime. △ Less

Submitted 22 March, 2017; v1 submitted 7 December, 2016; originally announced December 2016.

Comments: 21 pages, 19 figures

Spinning, Precessing, Black Hole Binary Spacetime via Asymptotic Matching

Authors: Hiroyuki Nakano, Brennan Ireland, Manuela Campanelli, Eric J. West

Abstract: We briefly discuss a method to construct a global, analytic, approximate spacetime for precessing, spinning binary black holes. The spacetime construction is broken into three parts: the inner zones are the spacetimes close to each black hole, and are approximated by perturbed Kerr solutions; the near zone is far from the two black holes, and described by the post-Newtonian metric; and finally the… ▽ More We briefly discuss a method to construct a global, analytic, approximate spacetime for precessing, spinning binary black holes. The spacetime construction is broken into three parts: the inner zones are the spacetimes close to each black hole, and are approximated by perturbed Kerr solutions; the near zone is far from the two black holes, and described by the post-Newtonian metric; and finally the wave (far) zone, where retardation effects need to be taken into account, is well modeled by the post-Minkowskian metric. These individual spacetimes are then stitched together using asymptotic matching techniques to obtain a global solution that approximately satisfies the Einstein field equations. Precession effects are introduced by rotating the black hole spin direction according to the precessing equations of motion, in a way that is consistent with the global spacetime construction. △ Less

Submitted 2 August, 2016; originally announced August 2016.

Comments: 9 pages

Modeling the source of GW150914 with targeted numerical-relativity simulations

Authors: Geoffrey Lovelace, Carlos O. Lousto, James Healy, Mark A. Scheel, Alyssa Garcia, Richard O'Shaughnessy, Michael Boyle, Manuela Campanelli, Daniel A. Hemberger, Lawrence E. Kidder, Harald P. Pfeiffer, Bela Szilagyi, Saul A. Teukolsky, Yosef Zlochower

Abstract: In fall of 2015, the two LIGO detectors measured the gravitational wave signal GW150914, which originated from a pair of merging black holes. In the final 0.2 seconds (about 8 gravitational-wave cycles) before the amplitude reached its maximum, the observed signal swept up in amplitude and frequency, from 35 Hz to 150 Hz. The theoretical gravitational-wave signal for merging black holes, as predic… ▽ More In fall of 2015, the two LIGO detectors measured the gravitational wave signal GW150914, which originated from a pair of merging black holes. In the final 0.2 seconds (about 8 gravitational-wave cycles) before the amplitude reached its maximum, the observed signal swept up in amplitude and frequency, from 35 Hz to 150 Hz. The theoretical gravitational-wave signal for merging black holes, as predicted by general relativity, can be computed only by full numerical relativity, because analytic approximations fail near the time of merger. Moreover, the nearly-equal masses, moderate spins, and small number of orbits of GW150914 are especially straightforward and efficient to simulate with modern numerical-relativity codes. In this paper, we report the modeling of GW150914 with numerical-relativity simulations, using black-hole masses and spins consistent with those inferred from LIGO's measurement. In particular, we employ two independent numerical-relativity codes that use completely different analytical and numerical methods to model the same merging black holes and to compute the emitted gravitational waveform; we find excellent agreement between the waveforms produced by the two independent codes. These results demonstrate the validity, impact, and potential of current and future studies using rapid-response, targeted numerical-relativity simulations for better understanding gravitational-wave observations. △ Less

Submitted 18 July, 2016; originally announced July 2016.

Comments: 11 pages, 3 figures, submitted to Classical and Quantum Gravity

Journal ref: Class. Quantum Grav. 33, 244002 (2016)

Directly comparing GW150914 with numerical solutions of Einstein's equations for binary black hole coalescence

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, M. R. Abernathy, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, A. Allocca, P. A. Altin, S. B. Anderson , et al. (957 additional authors not shown)

Abstract: We compare GW150914 directly to simulations of coalescing binary black holes in full general relativity, accounting for all the spin-weighted quadrupolar modes, and separately accounting for all the quadrupolar and octopolar modes. Consistent with the posterior distributions reported in LVC_PE[1] (at 90% confidence), we find the data are compatible with a wide range of nonprecessing and precessing… ▽ More We compare GW150914 directly to simulations of coalescing binary black holes in full general relativity, accounting for all the spin-weighted quadrupolar modes, and separately accounting for all the quadrupolar and octopolar modes. Consistent with the posterior distributions reported in LVC_PE[1] (at 90% confidence), we find the data are compatible with a wide range of nonprecessing and precessing simulations. Followup simulations performed using previously-estimated binary parameters most resemble the data. Comparisons including only the quadrupolar modes constrain the total redshifted mass Mz \in [64 - 82M_\odot], mass ratio q = m2/m1 \in [0.6,1], and effective aligned spin χ_eff \in [-0.3, 0.2], where χ_{eff} = (S1/m1 + S2/m2) \cdot\hat{L} /M. Including both quadrupolar and octopolar modes, we find the mass ratio is even more tightly constrained. Simulations with extreme mass ratios and effective spins are highly inconsistent with the data, at any mass. Several nonprecessing and precessing simulations with similar mass ratio and χ_{eff} are consistent with the data. Though correlated, the components' spins (both in magnitude and directions) are not significantly constrained by the data. For nonprecessing binaries, interpolating between simulations, we reconstruct a posterior distribution consistent with previous results. The final black hole's redshifted mass is consistent with Mf,z between 64.0 - 73.5M_\odot and the final black hole's dimensionless spin parameter is consistent with af = 0.62 - 0.73. As our approach invokes no intermediate approximations to general relativity and can strongly reject binaries whose radiation is inconsistent with the data, our analysis provides a valuable complement to LVC_PE[1]. △ Less

Submitted 3 June, 2016; originally announced June 2016.

Comments: 29 pages. Full version of tables available electronically with this submission

Report number: LIGO P1500263

Journal ref: Phys. Rev. D 94, 064035 (2016)

Properties of the Binary Black Hole Merger GW150914

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, M. R. Abernathy, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, A. Allocca , et al. (980 additional authors not shown)

Abstract: On September 14, 2015, the Laser Interferometer Gravitational-wave Observatory (LIGO) detected a gravitational-wave transient (GW150914); we characterize the properties of the source and its parameters. The data around the time of the event were analyzed coherently across the LIGO network using a suite of accurate waveform models that describe gravitational waves from a compact binary system in ge… ▽ More On September 14, 2015, the Laser Interferometer Gravitational-wave Observatory (LIGO) detected a gravitational-wave transient (GW150914); we characterize the properties of the source and its parameters. The data around the time of the event were analyzed coherently across the LIGO network using a suite of accurate waveform models that describe gravitational waves from a compact binary system in general relativity. GW150914 was produced by a nearly equal mass binary black hole of $36^{+5}_{-4} M_\odot$ and $29^{+4}_{-4} M_\odot$; for each parameter we report the median value and the range of the 90% credible interval. The dimensionless spin magnitude of the more massive black hole is bound to be $<0.7$ (at 90% probability). The luminosity distance to the source is $410^{+160}_{-180}$ Mpc, corresponding to a redshift $0.09^{+0.03}_{-0.04}$ assuming standard cosmology. The source location is constrained to an annulus section of $610$ deg$^2$, primarily in the southern hemisphere. The binary merges into a black hole of $62^{+4}_{-4} M_\odot$ and spin $0.67^{+0.05}_{-0.07}$. This black hole is significantly more massive than any other inferred from electromagnetic observations in the stellar-mass regime. △ Less

Submitted 20 June, 2016; v1 submitted 11 February, 2016; originally announced February 2016.

Comments: 19 pages, 6 figures

Report number: LIGO-P1500218

Journal ref: Phys. Rev. Lett. 116, 241102 (2016)

Inspiralling, Non-Precessing, Spinning Black Hole Binary Spacetime via Asymptotic Matching

Authors: Brennan Ireland, Bruno C. Mundim, Hiroyuki Nakano, Manuela Campanelli

Abstract: 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 vali… ▽ More 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. △ Less

Submitted 23 March, 2016; v1 submitted 17 December, 2015; originally announced December 2015.

Comments: 21 pages, 15 figures, revtex 4.1

Journal ref: Phys. Rev. D 93, 104057 (2016)

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

Authors: Yosef Zlochower, Hiroyuki Nakano, Bruno C. Mundim, Manuela Campanelli, Scott Noble, Miguel Zilhao

Abstract: We explore how a recently developed analytical black-hole binary spacetime can be extended using numerical simulations to go beyond the slow-inspiral phase. The analytic spacetime solves the Einstein field equations approximately, with the approximation error becoming progressively smaller the more separated the binary. To continue the spacetime beyond the slow-inspiral phase, we need to transitio… ▽ More We explore how a recently developed analytical black-hole binary spacetime can be extended using numerical simulations to go beyond the slow-inspiral phase. The analytic spacetime solves the Einstein field equations approximately, with the approximation error becoming progressively smaller the more separated the binary. To continue the spacetime beyond the slow-inspiral phase, we need to transition. Such a transition was previously explored at smaller separations. Here, we perform this transition at a separation of D=20M (large enough that the analytical metric is expected to be accurate), and evolve for six orbits. We find that small constraint violations can have large dynamical effects, but these can be removed by using a constraint-damping system like the conformal covariant formulation of the Z4 system. We find agreement between the subsequent numerical spacetime and the predictions of post-Newtonian theory for the waveform and inspiral rate that is within the post-Newtonian truncation error. △ Less

Submitted 23 June, 2016; v1 submitted 1 April, 2015; originally announced April 2015.

Comments: Accepted for publication in PRD, 14 pages, 21 figures, revtex 4-1

Journal ref: Phys. Rev. D 93, 124072 (2016)

Resolving the relative influence of strong field spacetime dynamics and MHD on circumbinary disk physics

Authors: Miguel Zilhão, Scott C. Noble, Manuela Campanelli, Yosef Zlochower

Abstract: In this paper we evolve magnetized and unmagnetized circumbinary accretion disks around supermassive black hole binaries in the relativistic regime. We use a post-Newtonian expansion to construct an analytical spacetime and determine how the order of the post-Newtonian (PN) expansion affects the dynamics of the gas. We find very small differences in the late-time bulk dynamics of non-magnetized hy… ▽ More In this paper we evolve magnetized and unmagnetized circumbinary accretion disks around supermassive black hole binaries in the relativistic regime. We use a post-Newtonian expansion to construct an analytical spacetime and determine how the order of the post-Newtonian (PN) expansion affects the dynamics of the gas. We find very small differences in the late-time bulk dynamics of non-magnetized hydrodynamic evolutions between the two spacetimes down to separations of approximately $40GM/c^2$ where $M$ is the total mass of the binary. For smaller separations, the differences due to PN-order become comparable to differences caused by using initial data further from equilibrium. For magnetized gas, MHD stresses, which drives the accretion dynamics, tends to mask all higher order PN effects even at separations of $20GM/c^2$, leading to essentially the same observed electromagnetic luminosity. This implies that our calculations of the EM signal may be robust down to small binary separations. Our investigation is the first to demonstrate how the level of PN accuracy affects a circumbinary disk's evolution and informs us of the range in separation within which to trust the PN approximation for this kind of study. We also address the influence the initial conditions and binary separation have on simulation predictions. △ Less

Submitted 31 January, 2015; v1 submitted 16 September, 2014; originally announced September 2014.

Comments: 14 pages, 12 figures. v2: matches published version in PRD

The NINJA-2 project: Detecting and characterizing gravitational waveforms modelled using numerical binary black hole simulations

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the NINJA-2 Collaboration, :, J. Aasi, B. P. Abbott, R. Abbott, T. Abbott, M. R. Abernathy, T. Accadia, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, C. Affeldt, M. Agathos, N. Aggarwal, O. D. Aguiar, A. Ain, P. Ajith, A. Alemic, B. Allen, A. Allocca , et al. (867 additional authors not shown)

Abstract: The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave astrophysics communities. The purpose of NINJA is to study the ability to detect gravitational waves emitted from merging binary black holes and recover their parameters with next-generation gravitational-wave observatories. We report here on the results of… ▽ More The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave astrophysics communities. The purpose of NINJA is to study the ability to detect gravitational waves emitted from merging binary black holes and recover their parameters with next-generation gravitational-wave observatories. We report here on the results of the second NINJA project, NINJA-2, which employs 60 complete binary black hole hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. In a "blind injection challenge" similar to that conducted in recent LIGO and Virgo science runs, we added 7 hybrid waveforms to two months of data recolored to predictions of Advanced LIGO and Advanced Virgo sensitivity curves during their first observing runs. The resulting data was analyzed by gravitational-wave detection algorithms and 6 of the waveforms were recovered with false alarm rates smaller than 1 in a thousand years. Parameter estimation algorithms were run on each of these waveforms to explore the ability to constrain the masses, component angular momenta and sky position of these waveforms. We also perform a large-scale monte-carlo study to assess the ability to recover each of the 60 hybrid waveforms with early Advanced LIGO and Advanced Virgo sensitivity curves. Our results predict that early Advanced LIGO and Advanced Virgo will have a volume-weighted average sensitive distance of 300Mpc (1Gpc) for $10M_{\odot}+10M_{\odot}$ ($50M_{\odot}+50M_{\odot}$) binary black hole coalescences. We demonstrate that neglecting the component angular momenta in the waveform models used in matched-filtering will result in a reduction in sensitivity for systems with large component angular momenta. [Abstract abridged for ArXiv, full version in PDF] △ Less

Submitted 5 January, 2014; originally announced January 2014.

Comments: 52 pages, 11 figures, 7 tables

Report number: LIGO-P1300199

Journal ref: Class. Quantum Grav. 31 115004, 2014

Circumbinary MHD Accretion into Inspiraling Binary Black Holes

Authors: Scott C. Noble, Bruno C. Mundim, Hiroyuki Nakano, Julian H. Krolik, Manuela Campanelli, Yosef Zlochower, Nicolás Yunes

Abstract: As 2 black holes bound to each other in a close binary approach merger their inspiral time becomes shorter than the characteristic inflow time of surrounding orbiting matter. Using an innovative technique in which we represent the changing spacetime in the region occupied by the orbiting matter with a 2.5PN approximation and the binary orbital evolution with 3.5PN, we have simulated the MHD evolut… ▽ More As 2 black holes bound to each other in a close binary approach merger their inspiral time becomes shorter than the characteristic inflow time of surrounding orbiting matter. Using an innovative technique in which we represent the changing spacetime in the region occupied by the orbiting matter with a 2.5PN approximation and the binary orbital evolution with 3.5PN, we have simulated the MHD evolution of a circumbinary disk surrounding an equal-mass non-spinning binary. Prior to the beginning of the inspiral, the structure of the circumbinary disk is predicted well by extrapolation from Newtonian results. The binary opens a low-density gap whose radius is roughly two binary separations, and matter piles up at the outer edge of this gap as inflow is retarded by torques exerted by the binary; nonetheless, the accretion rate is diminished relative to its value at larger radius by only about a factor of 2. During inspiral, the inner edge of the disk at first moves inward in coordination with the shrinking binary, but as the orbital evolution accelerates, the rate at which the inner edge moves toward smaller radii falls behind the rate of binary compression. In this stage, the rate of angular momentum transfer from the binary to the disk slows substantially, but the net 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 in galactic nuclei could be very luminous at this stage of their evolution. If the luminosity were optically thin, it would be modulated at a frequency that is a beat between the orbital frequency of the disk's surface density maximum and the binary orbital frequency. However, a disk with sufficient surface density to be luminous should also be optically thick; as a result, the periodic modulation may be suppressed. △ Less

Submitted 4 April, 2012; originally announced April 2012.

Comments: 54 pages, color figures, submitted to the Astrophysical Journal, a high resolution version and movies can be found at http://ccrg.rit.edu/~scn/cmhdaiibh/

Study of Conformally Flat Initial Data for Highly Spinning Black Holes and their Early Evolutions

Authors: Carlos O. Lousto, Hiroyuki Nakano, Yosef Zlochower, Bruno C. Mundim, Manuela Campanelli

Abstract: We study conformally-flat initial data for an arbitrary number of spinning black holes with exact analytic solutions to the momentum constraints constructed from a linear combination of the classical Bowen-York and conformal Kerr extrinsic curvatures. The solution leading to the largest intrinsic spin, relative to the ADM mass of the spacetime epsilon_S=S/M^2_{ADM}, is a superposition with relativ… ▽ More We study conformally-flat initial data for an arbitrary number of spinning black holes with exact analytic solutions to the momentum constraints constructed from a linear combination of the classical Bowen-York and conformal Kerr extrinsic curvatures. The solution leading to the largest intrinsic spin, relative to the ADM mass of the spacetime epsilon_S=S/M^2_{ADM}, is a superposition with relative weights of Lambda=0.783 for conformal Kerr and (1-Lambda)=0.217 for Bowen-York. In addition, we measure the spin relative to the initial horizon mass M_{H_0}, and find that the quantity chi=S/M_{H_0}^2 reaches a maximum of χ^{max}=0.9856 for Lambda=0.753. After equilibration, the final black-hole spin should lie in the interval 0.9324<chi_{final}<0.9856. We perform full numerical evolutions to compute the energy radiated and the final horizon mass and spin. We find that the black hole settles to a final spin of chi_{final}^{max}=0.935 when Lambda=0.783. We also study the evolution of the apparent horizon structure of this "maximal" black hole in detail. △ Less

Submitted 14 June, 2012; v1 submitted 14 March, 2012; originally announced March 2012.

Comments: 9 pages, 8 figures

The Einstein Toolkit: A Community Computational Infrastructure for Relativistic Astrophysics

Authors: Frank Löffler, Joshua Faber, Eloisa Bentivegna, Tanja Bode, Peter Diener, Roland Haas, Ian Hinder, Bruno C. Mundim, Christian D. Ott, Erik Schnetter, Gabrielle Allen, Manuela Campanelli, Pablo Laguna

Abstract: 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, c… ▽ More 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. △ Less

Submitted 14 November, 2011; originally announced November 2011.

Comments: 62 pages, 20 figures

Intermediate-mass-ratio black hole binaries II: Modeling Trajectories and Gravitational Waveforms

Authors: Hiroyuki Nakano, Yosef Zlochower, Carlos O. Lousto, Manuela Campanelli

Abstract: We revisit the scenario of small-mass-ratio (q) black-hole binaries; performing new, more accurate, simulations of mass ratios 10:1 and 100:1 for initially nonspinning black holes. We propose fitting functions for the trajectories of the two black holes as a function of time and mass ratio (in the range 1/100 < q < 1/10$) that combine aspects of post-Newtonian trajectories at smaller orbital frequ… ▽ More We revisit the scenario of small-mass-ratio (q) black-hole binaries; performing new, more accurate, simulations of mass ratios 10:1 and 100:1 for initially nonspinning black holes. We propose fitting functions for the trajectories of the two black holes as a function of time and mass ratio (in the range 1/100 < q < 1/10$) that combine aspects of post-Newtonian trajectories at smaller orbital frequencies and plunging geodesics at larger frequencies. We then use these trajectories to compute waveforms via black hole perturbation theory. Using the advanced LIGO noise curve, we see a match of ~99.5% for the leading (l,m)=(2,2) mode between the numerical relativity and perturbative waveforms. Nonleading modes have similarly high matches. We thus prove the feasibility of efficiently generating a bank of gravitational waveforms in the intermediate-mass-ratio regime using only a sparse set of full numerical simulations. △ Less

Submitted 22 August, 2011; originally announced August 2011.

Comments: 23 pages, 35 figures, revtex4

Journal ref: Phys. Rev. D 84, 124006 (2011)

Modeling Gravitational Recoil Using Numerical Relativity

Authors: Yosef Zlochower, Manuela Campanelli, Carlos O. Lousto

Abstract: We review the developments in modeling gravitational recoil from merging black-hole binaries and introduce a new set of 20 simulations to test our previously proposed empirical formula for the recoil. The configurations are chosen to represent generic binaries with unequal masses and precessing spins. Results of these simulations indicate that the recoil formula is accurate to within a few km/s in… ▽ More We review the developments in modeling gravitational recoil from merging black-hole binaries and introduce a new set of 20 simulations to test our previously proposed empirical formula for the recoil. The configurations are chosen to represent generic binaries with unequal masses and precessing spins. Results of these simulations indicate that the recoil formula is accurate to within a few km/s in the similar mass-ratio regime for the out-of-plane recoil. △ Less

Submitted 11 November, 2010; v1 submitted 9 November, 2010; originally announced November 2010.

Comments: corrections to text, 11 pages, 1 figure

Journal ref: Class.Quant.Grav.28:114015,2011

Intermediate-mass-ratio black hole binaries: intertwining numerical and perturbative techniques

Authors: Carlos O. Lousto, Hiroyuki Nakano, Yosef Zlochower, Manuela Campanelli

Abstract: We describe in detail full numerical and perturbative techniques to compute the gravitational radiation from intermediate-mass-ratio black-hole-binary inspirals and mergers. We perform a series of full numerical simulations of nonspinning black holes with mass ratios q=1/10 and q=1/15 from different initial separations and for different finite-difference resolutions. In order to perform those full… ▽ More We describe in detail full numerical and perturbative techniques to compute the gravitational radiation from intermediate-mass-ratio black-hole-binary inspirals and mergers. We perform a series of full numerical simulations of nonspinning black holes with mass ratios q=1/10 and q=1/15 from different initial separations and for different finite-difference resolutions. In order to perform those full numerical runs, we adapt the gauge of the moving punctures approach with a variable damping term for the shift. We also derive an extrapolation (to infinite radius) formula for the waveform extracted at finite radius. For the perturbative evolutions we use the full numerical tracks, transformed into the Schwarzschild gauge, in the source terms of the Regge-Wheller-Zerilli Schwarzschild perturbations formalism. We then extend this perturbative formalism to take into account small intrinsic spins of the large black hole, and validate it by computing the quasinormal mode frequencies, where we find good agreement for spins |a/M|<0.3. Including the final spins improves the overlap functions when comparing full numerical and perturbative waveforms, reaching 99.5% for the leading (l,m)=(2,2) and (3,3) modes, and 98.3% for the nonleading (2,1) mode in the q=1/10 case, which includes 8 orbits before merger. For the q=1/15 case, we obtain overlaps near 99.7% for all three modes. We discuss the modeling of the full inspiral and merger based on a combined matching of post-Newtonian, full numerical, and geodesic trajectories. △ Less

Submitted 2 December, 2010; v1 submitted 25 August, 2010; originally announced August 2010.

Comments: 31 pages, 33 figures revtex 4, Published version. Several fixes throughout

Journal ref: Phys.Rev.D82:104057,2010

Intermediate Mass Ratio Black Hole Binaries: Numerical Relativity meets Perturbation Theory

Authors: Carlos O. Lousto, Hiroyuki Nakano, Yosef Zlochower, Manuela Campanelli

Abstract: We study black-hole binaries in the intermediate-mass-ratio regime 0.01 < q < 0.1 with a new technique that makes use of nonlinear numerical trajectories and efficient perturbative evolutions to compute waveforms at large radii for the leading and nonleading modes. As a proof-of-concept, we compute waveforms for q=1/10. We discuss applications of these techniques for LIGO/VIRGO data analysis and… ▽ More We study black-hole binaries in the intermediate-mass-ratio regime 0.01 < q < 0.1 with a new technique that makes use of nonlinear numerical trajectories and efficient perturbative evolutions to compute waveforms at large radii for the leading and nonleading modes. As a proof-of-concept, we compute waveforms for q=1/10. We discuss applications of these techniques for LIGO/VIRGO data analysis and the possibility that our technique can be extended to produce accurate waveform templates from a modest number of fully-nonlinear numerical simulations. △ Less

Submitted 13 January, 2010; originally announced January 2010.

Comments: 4 pages, 5 figures, revtex4

Journal ref: Phys.Rev.Lett.104:211101,2010

Statistical studies of Spinning Black-Hole Binaries

Authors: Carlos O. Lousto, Hiroyuki Nakano, Yosef Zlochower, Manuela Campanelli

Abstract: We study the statistical distributions of the spins of generic black-hole binaries during the inspiral and merger, as well as the distributions of the remnant mass, spin, and recoil velocity. For the inspiral regime, we start with a random uniform distribution of spin directions S1 and S2 and magnitudes S1=S2=0.97 for different mass ratios. Starting from a fiducial initial separation of ri=50m,… ▽ More We study the statistical distributions of the spins of generic black-hole binaries during the inspiral and merger, as well as the distributions of the remnant mass, spin, and recoil velocity. For the inspiral regime, we start with a random uniform distribution of spin directions S1 and S2 and magnitudes S1=S2=0.97 for different mass ratios. Starting from a fiducial initial separation of ri=50m, we perform 3.5PN evolutions down to rf=5m. At this final fiducial separation, we compute the angular distribution of the spins with respect to the final orbital angular momentum, L. We perform 16^4 simulations for six mass ratios between q=1 and q=1/16 and compute the distribution of the angles between L and Delta and L and S, directly related to recoil velocities and total angular momentum. We find a small but statistically significant bias of the distribution towards counter-alignment of both scalar products. To study the merger of black-hole binaries, we turn to full numerical techniques. We introduce empirical formulae to describe the final remnant black hole mass, spin, and recoil velocity for merging black-hole binaries with arbitrary mass ratios and spins. We then evaluate those formulae for randomly chosen directions of the individual spins and magnitudes as well as the binary's mass ratio. We found that the magnitude of the recoil velocity distribution decays as P(v) \exp(-v/2500km/s), <v>=630km/s, and sqrt{<v^2> - <v>^2}= 534km/s, leading to a 23% probability of recoils larger than 1000km/s, and a highly peaked angular distribution along the final orbital axis. The final black-hole spin magnitude show a universal distribution highly peaked at Sf/mf^2=0.73 and a 25 degrees misalignment with respect to the final orbital angular momentum. △ Less

Submitted 21 January, 2010; v1 submitted 16 October, 2009; originally announced October 2009.

Comments: 22 pages, 28 figures, abridged abstract. Major addenda in section IV and discussion

Journal ref: Phys.Rev.D81:084023,2010; Erratum-ibid.D82:129902,2010

Remnant Masses, Spins and Recoils from the Merger of Generic Black-Hole Binaries

Authors: Carlos O. Lousto, Manuela Campanelli, Yosef Zlochower, Hiroyuki Nakano

Abstract: We obtain empirical formulae for the final remnant black hole mass, spin, and recoil velocity from merging black-hole binaries with arbitrary mass ratios and spins. Our formulae are based on the mass ratio and spin dependence of the post-Newtonian expressions for the instantaneous radiated energy, linear momentum, and angular momentum, as well as the ISCO binding energy and angular momentum. The… ▽ More We obtain empirical formulae for the final remnant black hole mass, spin, and recoil velocity from merging black-hole binaries with arbitrary mass ratios and spins. Our formulae are based on the mass ratio and spin dependence of the post-Newtonian expressions for the instantaneous radiated energy, linear momentum, and angular momentum, as well as the ISCO binding energy and angular momentum. The relative weight between the different terms is fixed by amplitude parameters chosen through a least-squares fit of recently available fully nonlinear numerical simulations. These formulae can be used for statistical studies of N-body simulations of galaxy cores and clusters, and the cosmological growth of supermassive black holes. As an example, we use these formulae to obtain a universal spin magnitude distribution of merged black holes and recoil velocity distributions for dry and hot/cold wet mergers. We also revisit the long term orbital precession and resonances and discuss how they affect spin distributions before the merging regime. △ Less

Submitted 20 April, 2010; v1 submitted 22 April, 2009; originally announced April 2009.

Comments: 15 pages, 2 figures. Rewritten and added new section.

Journal ref: Class.Quant.Grav.27:114006,2010

Algebraic Classification of Numerical Spacetimes and Black-Hole-Binary Remnants

Authors: Manuela Campanelli, Carlos O. Lousto, Yosef Zlochower

Abstract: In this paper we develop a technique for determining the algebraic classification of a numerical spacetime, possibly resulting from a generic black-hole-binary merger, using the Newman-Penrose Weyl scalars. We demonstrate these techniques for a test case involving a close binary with arbitrarily oriented spins and unequal masses. We find that, post merger, the spacetime quickly approaches Petrov… ▽ More In this paper we develop a technique for determining the algebraic classification of a numerical spacetime, possibly resulting from a generic black-hole-binary merger, using the Newman-Penrose Weyl scalars. We demonstrate these techniques for a test case involving a close binary with arbitrarily oriented spins and unequal masses. We find that, post merger, the spacetime quickly approaches Petrov type II, and only approaches type D on much longer timescales. These techniques allow us to begin to explore the validity of the "no-hair theorem" for generic merging-black-hole spacetimes. △ Less

Submitted 7 May, 2009; v1 submitted 18 November, 2008; originally announced November 2008.

Comments: published version

Journal ref: Phys.Rev.D79:084012,2009

Comparison of Numerical and Post-Newtonian Waveforms for Generic Precessing Black-Hole Binaries

Authors: Manuela Campanelli, Carlos O. Lousto, Hiroyuki Nakano, Yosef Zlochower

Abstract: We compare waveforms and orbital dynamics from the first long-term, fully non-linear, numerical simulations of a generic black-hole binary configuration with post-Newtonian predictions. The binary has mass ratio q~0.8 with arbitrarily oriented spins of magnitude S_1/m_1^2~0.6 and S_2/m_2^2~0.4 and orbits 9 times prior to merger. The numerical simulation starts with an initial separation of r~11M… ▽ More We compare waveforms and orbital dynamics from the first long-term, fully non-linear, numerical simulations of a generic black-hole binary configuration with post-Newtonian predictions. The binary has mass ratio q~0.8 with arbitrarily oriented spins of magnitude S_1/m_1^2~0.6 and S_2/m_2^2~0.4 and orbits 9 times prior to merger. The numerical simulation starts with an initial separation of r~11M, with orbital parameters determined by initial 2.5PN and 3.5PN post-Newtonian evolutions of a quasi-circular binary with an initial separation of r=50M. The resulting binaries have very little eccentricity according to the 2.5PN and 3.5PN systems, but show significant eccentricities of e~0.01-0.02 and e~0.002-0.005 in the respective numerical simulations, thus demonstrating that 3.5PN significantly reduces the eccentricity of the binary compared to 2.5PN. We perform three numerical evolutions from r~11M with maximum resolutions of h=M/48,M/53.3,M/59.3, to verify numerical convergence. We observe a reasonably good agreement between the PN and numerical waveforms, with an overlap of nearly 99% for the first six cycles of the (l=2,m=+-2) modes, 91% for the (l=2,m=+-1) modes, and nearly 91% for the (l=3,m=+-3) modes. The phase differences between numerical and post-Newtonian approximations appear to be independent of the (l,m) modes considered and relatively small for the first 3-4 orbits. An advantage of the 3.5 PN model over the 2.5 PN one seems to be observed, which indicates that still higher PN order (perhaps even 4.0PN) may yield significantly better waveforms. In addition, we identify features in the waveforms likely related to precession and precession-induced eccentricity. △ Less

Submitted 1 January, 2009; v1 submitted 5 August, 2008; originally announced August 2008.

Comments: New figures, enhanced analysis, revisions throughout the paper

Journal ref: Phys.Rev.D79:084010,2009

Close encounters of three black holes

Authors: Manuela Campanelli, Carlos O. Lousto, Yosef Zlochower

Abstract: We present the first fully relativistic longterm numerical evolutions of three equal-mass black holes in a system consisting of a third black hole in a close orbit about a black-hole binary. We find that these close-three-black-hole systems have very different merger dynamics from black-hole binaries. In particular, we see complex trajectories, a redistribution of energy that can impart substant… ▽ More We present the first fully relativistic longterm numerical evolutions of three equal-mass black holes in a system consisting of a third black hole in a close orbit about a black-hole binary. We find that these close-three-black-hole systems have very different merger dynamics from black-hole binaries. In particular, we see complex trajectories, a redistribution of energy that can impart substantial kicks to one of the holes, distinctive waveforms, and suppression of the emitted gravitational radiation. We evolve two such configurations and find very different behaviors. In one configuration the binary is quickly disrupted and the individual holes follow complicated trajectories and merge with the third hole in rapid succession, while in the other, the binary completes a half-orbit before the initial merger of one of the members with the third black hole, and the resulting two-black-hole system forms a highly elliptical, well separated binary that shows no significant inspiral for (at least) the first t~1000M of evolution. △ Less

Submitted 9 April, 2008; v1 submitted 3 October, 2007; originally announced October 2007.

Comments: 4 pages, 5 figures

Journal ref: Phys.Rev.D77:101501,2008

Maximum gravitational recoil

Authors: Manuela Campanelli, Carlos O. Lousto, Yosef Zlochower, David Merritt

Abstract: Recent calculations of gravitational radiation recoil generated during black-hole binary mergers have reopened the possibility that a merged binary can be ejected even from the nucleus of a massive host galaxy. Here we report the first systematic study of gravitational recoil of equal-mass binaries with equal, but anti-aligned, spins parallel to the orbital plane. Such an orientation of the spin… ▽ More Recent calculations of gravitational radiation recoil generated during black-hole binary mergers have reopened the possibility that a merged binary can be ejected even from the nucleus of a massive host galaxy. Here we report the first systematic study of gravitational recoil of equal-mass binaries with equal, but anti-aligned, spins parallel to the orbital plane. Such an orientation of the spins is expected to maximize the recoil. We find that recoil velocity (which is perpendicular to the orbital plane) varies sinusoidally with the angle that the initial spin directions make with the initial linear momenta of each hole and scales up to a maximum of ~4000 km/s for maximally-rotating holes. Our results show that the amplitude of the recoil velocity can depend sensitively on spin orientations of the black holes prior to merger. △ Less

Submitted 27 February, 2007; v1 submitted 25 February, 2007; originally announced February 2007.

Comments: 4 pages, 4 figs, revtex4

Journal ref: Phys.Rev.Lett.98:231102,2007