Showing 1–50 of 70 results for author: Noble, S

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.

The Black Hole Explorer: Motivation and Vision

Authors: Michael D. Johnson, Kazunori Akiyama, Rebecca Baturin, Bryan Bilyeu, Lindy Blackburn, Don Boroson, Alejandro Cardenas-Avendano, Andrew Chael, Chi-kwan Chan, Dominic Chang, Peter Cheimets, Cathy Chou, Sheperd S. Doeleman, Joseph Farah, Peter Galison, Ronald Gamble, Charles F. Gammie, Zachary Gelles, Jose L. Gomez, Samuel E. Gralla, Paul Grimes, Leonid I. Gurvits, Shahar Hadar, Kari Haworth, Kazuhiro Hada , et al. (43 additional authors not shown)

Abstract: We present the Black Hole Explorer (BHEX), a mission that will produce the sharpest images in the history of astronomy by extending submillimeter Very-Long-Baseline Interferometry (VLBI) to space. BHEX will discover and measure the bright and narrow "photon ring" that is predicted to exist in images of black holes, produced from light that has orbited the black hole before escaping. This discovery… ▽ More We present the Black Hole Explorer (BHEX), a mission that will produce the sharpest images in the history of astronomy by extending submillimeter Very-Long-Baseline Interferometry (VLBI) to space. BHEX will discover and measure the bright and narrow "photon ring" that is predicted to exist in images of black holes, produced from light that has orbited the black hole before escaping. This discovery will expose universal features of a black hole's spacetime that are distinct from the complex astrophysics of the emitting plasma, allowing the first direct measurements of a supermassive black hole's spin. In addition to studying the properties of the nearby supermassive black holes M87* and Sgr A*, BHEX will measure the properties of dozens of additional supermassive black holes, providing crucial insights into the processes that drive their creation and growth. BHEX will also connect these supermassive black holes to their relativistic jets, elucidating the power source for the brightest and most efficient engines in the universe. BHEX will address fundamental open questions in the physics and astrophysics of black holes that cannot be answered without submillimeter space VLBI. The mission is enabled by recent technological breakthroughs, including the development of ultra-high-speed downlink using laser communications, and it leverages billions of dollars of existing ground infrastructure. We present the motivation for BHEX, its science goals and associated requirements, and the pathway to launch within the next decade. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: Proceedings for SPIE Astronomical Telescopes and Instrumentation

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

The Dynamics of Debris Disk Creation in Neutron Star Mergers

Authors: Yossef Zenati, Julian Krolik, Leonardo Werneck, Zachariah Etienne, Scott Noble, Ariadna Murguia-Berthier, Jeremy Schnittman

Abstract: The detection of GW170817/AT2017gfo inaugurated an era of multimessenger astrophysics, in which gravitational wave and multiwavelength photon observations complement one another to provide unique insight on astrophysical systems. A broad theoretical consensus exists in which the photon phenomenology of neutron star mergers largely rests upon the evolution of the small amount of matter left on boun… ▽ More The detection of GW170817/AT2017gfo inaugurated an era of multimessenger astrophysics, in which gravitational wave and multiwavelength photon observations complement one another to provide unique insight on astrophysical systems. A broad theoretical consensus exists in which the photon phenomenology of neutron star mergers largely rests upon the evolution of the small amount of matter left on bound orbits around the black hole or massive neutron star remaining after the merger. Because this accretion disk is far from inflow equilibrium, its subsequent evolution depends very strongly on its initial state, yet very little is known about how this state is determined. Using both snapshot and tracer particle data from a numerical relativity/MHD simulation of an equal-mass neutron star merger that collapses to a black hole, we show how gravitational forces arising in a non-axisymmetric, dynamical spacetime supplement hydrodynamical effects in shaping the initial structure of the bound debris disk. The work done by hydrodynamical forces is ${\sim}10$ times greater than that due to time-dependent gravity. Although gravitational torques prior to remnant relaxation are an order of magnitude larger than hydrodynamical torques, their intrinsic sign symmetry leads to strong cancellation; as a result, hydrodynamical and gravitational torques have comparable effect. We also show that the debris disk's initial specific angular momentum distribution is sharply peaked at roughly the specific angular momentum of the merged neutron star's outer layers, a few $r_g c$, and identify the regulating mechanism. △ Less

Submitted 4 June, 2024; v1 submitted 3 April, 2024; originally announced April 2024.

Comments: 20 pages, 10 figures

Bound Debris Expulsion from Neutron Star Merger Remnants

Authors: Yossef Zenati, Julian H. Krolik, Leonardo R. Werneck, Ariadna Murguia-Berthier, Zachariah B. Etienne, Scott C. Noble, Tsvi Piran

Abstract: Many studies have found that neutron star mergers leave a fraction of the stars' mass in bound orbits surrounding the resulting massive neutron star or black hole. This mass is a site of $r-$ process nucleosynthesis and can generate a wind that contributes to a kilonova. However, comparatively little is known about the dynamics determining its mass or initial structure. Here we begin to investigat… ▽ More Many studies have found that neutron star mergers leave a fraction of the stars' mass in bound orbits surrounding the resulting massive neutron star or black hole. This mass is a site of $r-$ process nucleosynthesis and can generate a wind that contributes to a kilonova. However, comparatively little is known about the dynamics determining its mass or initial structure. Here we begin to investigate these questions, starting with the origin of the disk mass. Using tracer particle as well as discretized fluid data from numerical simulations, we identify where in the neutron stars the debris came from, the paths it takes in order to escape from the neutron stars' interiors, and the times and locations at which its orbital properties diverge from those of neighboring fluid elements that end up remaining in the merged neutron star. △ Less

Submitted 15 June, 2023; originally announced June 2023.

Comments: Comments welcome!

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

Shocks Power Tidal Disruption Events

Authors: Taeho Ryu, Julian Krolik, Tsvi Piran, Scott Noble, Mark Avara

Abstract: Accretion of debris seems to be the natural mechanism to power the radiation emitted during a tidal disruption event (TDE), in which a supermassive black hole tears apart a star. However, this requires the prompt formation of a compact accretion disk. Here, using a fully relativistic global simulation for the long-term evolution of debris in a TDE with realistic initial conditions, we show that at… ▽ More Accretion of debris seems to be the natural mechanism to power the radiation emitted during a tidal disruption event (TDE), in which a supermassive black hole tears apart a star. However, this requires the prompt formation of a compact accretion disk. Here, using a fully relativistic global simulation for the long-term evolution of debris in a TDE with realistic initial conditions, we show that at most a tiny fraction of the bound mass enters such a disk on the timescale of observed flares. To "circularize" most of the bound mass entails an increase in the binding energy of that mass by a factor $\sim 30$; we find at most an order unity change. Our simulation suggests it would take a time scale comparable to a few tens of the characteristic mass fallback time to dissipate enough energy for "circularization". Instead, the bound debris forms an extended eccentric accretion flow with eccentricity $\simeq 0.4-0.5$ by $\sim 2$ fallback times. Although the energy dissipated in shocks in this large-scale flow is much smaller than the "circularization" energy, it matches the observed radiated energy very well. Nonetheless, the impact of shocks is not strong enough to unbind initially bound debris into an outflow. △ Less

Submitted 22 August, 2023; v1 submitted 9 May, 2023; originally announced May 2023.

Comments: 20 pages, 15 figures, Submitted to ApJ, comments welcome!

Addition of tabulated equation of state and neutrino leakage support to IllinoisGRMHD

Authors: Leonardo R. Werneck, Zachariah B. Etienne, Ariadna Murguia-Berthier, Roland Haas, Federico Cipolletta, Scott C. Noble, Lorenzo Ennoggi, Federico G. Lopez Armengol, Bruno Giacomazzo, Thiago Assumpção, Joshua Faber, Tanmayee Gupte, Bernard J. Kelly, Julian H. Krolik

Abstract: We have added support for realistic, microphysical, finite-temperature equations of state (EOS) and neutrino physics via a leakage scheme to IllinoisGRMHD, an open-source GRMHD code for dynamical spacetimes in the Einstein Toolkit. These new features are provided by two new, NRPy+-based codes: NRPyEOS, which performs highly efficient EOS table lookups and interpolations, and NRPyLeakage, which imp… ▽ More We have added support for realistic, microphysical, finite-temperature equations of state (EOS) and neutrino physics via a leakage scheme to IllinoisGRMHD, an open-source GRMHD code for dynamical spacetimes in the Einstein Toolkit. These new features are provided by two new, NRPy+-based codes: NRPyEOS, which performs highly efficient EOS table lookups and interpolations, and NRPyLeakage, which implements a new, AMR-capable neutrino leakage scheme in the Einstein Toolkit. We have performed a series of strenuous validation tests that demonstrate the robustness of these new codes, particularly on the Cartesian AMR grids provided by Carpet. Furthermore, we show results from fully dynamical GRMHD simulations of single unmagnetized neutron stars, and magnetized binary neutron star mergers. This new version of IllinoisGRMHD, as well as NRPyEOS and NRPyLeakage, is pedagogically documented in Jupyter notebooks and fully open source. The codes will be proposed for inclusion in an upcoming version of the Einstein Toolkit. △ Less

Submitted 14 December, 2023; v1 submitted 30 August, 2022; originally announced August 2022.

Comments: 20 pages, 9 figures. v3 erratum: fixed typo in Eq. (21)

Journal ref: Phys. Rev. D 107, 044037 (2023)

Advancing the Landscape of Multimessenger Science in the Next Decade

Authors: Kristi Engel, Tiffany Lewis, Marco Stein Muzio, Tonia M. Venters, Markus Ahlers, Andrea Albert, Alice Allen, Hugo Alberto Ayala Solares, Samalka Anandagoda, Thomas Andersen, Sarah Antier, David Alvarez-Castillo, Olaf Bar, Dmitri Beznosko, Łukasz Bibrzyck, Adam Brazier, Chad Brisbois, Robert Brose, Duncan A. Brown, Mattia Bulla, J. Michael Burgess, Eric Burns, Cecilia Chirenti, Stefano Ciprini, Roger Clay , et al. (69 additional authors not shown)

Abstract: The last decade has brought about a profound transformation in multimessenger science. Ten years ago, facilities had been built or were under construction that would eventually discover the nature of objects in our universe could be detected through multiple messengers. Nonetheless, multimessenger science was hardly more than a dream. The rewards for our foresight were finally realized through Ice… ▽ More The last decade has brought about a profound transformation in multimessenger science. Ten years ago, facilities had been built or were under construction that would eventually discover the nature of objects in our universe could be detected through multiple messengers. Nonetheless, multimessenger science was hardly more than a dream. The rewards for our foresight were finally realized through IceCube's discovery of the diffuse astrophysical neutrino flux, the first observation of gravitational waves by LIGO, and the first joint detections in gravitational waves and photons and in neutrinos and photons. Today we live in the dawn of the multimessenger era. The successes of the multimessenger campaigns of the last decade have pushed multimessenger science to the forefront of priority science areas in both the particle physics and the astrophysics communities. Multimessenger science provides new methods of testing fundamental theories about the nature of matter and energy, particularly in conditions that are not reproducible on Earth. This white paper will present the science and facilities that will provide opportunities for the particle physics community renew its commitment and maintain its leadership in multimessenger science. △ Less

Submitted 18 March, 2022; originally announced March 2022.

Comments: 174 pages, 12 figures. Contribution to Snowmass 2021. Solicited white paper from CF07. Comments and endorsers welcome. Still accepting contributions (contact editors)

Astrophysics with the Laser Interferometer Space Antenna

Authors: Pau Amaro Seoane, Jeff Andrews, Manuel Arca Sedda, Abbas Askar, Quentin Baghi, Razvan Balasov, Imre Bartos, Simone S. Bavera, Jillian Bellovary, Christopher P. L. Berry, Emanuele Berti, Stefano Bianchi, Laura Blecha, Stephane Blondin, Tamara Bogdanović, Samuel Boissier, Matteo Bonetti, Silvia Bonoli, Elisa Bortolas, Katelyn Breivik, Pedro R. Capelo, Laurentiu Caramete, Federico Cattorini, Maria Charisi, Sylvain Chaty , et al. (134 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery… ▽ More The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe. △ Less

Submitted 25 May, 2023; v1 submitted 11 March, 2022; originally announced March 2022.

Journal ref: Living Reviews in Relativity, Volume 26, Article number: 2 (2023)

PATOKA: Simulating Electromagnetic Observables of Black Hole Accretion

Authors: George N. Wong, Ben S. Prather, Vedant Dhruv, Benjamin R. Ryan, Monika Moscibrodzka, Chi-kwan Chan, Abhishek V. Joshi, Ricardo Yarza, Angelo Ricarte, Hotaka Shiokawa, Joshua C. Dolence, Scott C. Noble, Jonathan C. McKinney, Charles F. Gammie

Abstract: The Event Horizon Telescope (EHT) has released analyses of reconstructed images of horizon-scale millimeter emission near the supermassive black hole at the center of the M87 galaxy. Parts of the analyses made use of a large library of synthetic black hole images and spectra, which were produced using numerical general relativistic magnetohydrodynamics fluid simulations and polarized ray tracing.… ▽ More The Event Horizon Telescope (EHT) has released analyses of reconstructed images of horizon-scale millimeter emission near the supermassive black hole at the center of the M87 galaxy. Parts of the analyses made use of a large library of synthetic black hole images and spectra, which were produced using numerical general relativistic magnetohydrodynamics fluid simulations and polarized ray tracing. In this article, we describe the PATOKA pipeline, which was used to generate the Illinois contribution to the EHT simulation library. We begin by describing the relevant accretion systems and radiative processes. We then describe the details of the three numerical codes we use, iharm, ipole, and igrmonty, paying particular attention to differences between the current generation of the codes and the originally published versions. Finally, we provide a brief overview of simulated data as produced by PATOKA and conclude with a discussion of limitations and future directions. △ Less

Submitted 11 March, 2022; v1 submitted 23 February, 2022; originally announced February 2022.

Comments: 42 pages, 17 figures, accepted for publication in ApJS

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)

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

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

Spin and Accretion Rate Dependence of Black Hole X-Ray Spectra

Authors: Brooks E. Kinch, Jeremy D. Schnittman, Scott C. Noble, Timothy R. Kallman, Julian H. Krolik

Abstract: We present a survey of how the spectral features of black hole X-ray binary systems depend on spin, accretion rate, viewing angle, and Fe abundance when predicted on the basis of first principles physical calculations. The power law component hardens with increasing spin. The thermal component strengthens with increasing accretion rate. The Compton bump is enhanced by higher accretion rate and low… ▽ More We present a survey of how the spectral features of black hole X-ray binary systems depend on spin, accretion rate, viewing angle, and Fe abundance when predicted on the basis of first principles physical calculations. The power law component hardens with increasing spin. The thermal component strengthens with increasing accretion rate. The Compton bump is enhanced by higher accretion rate and lower spin. The Fe K$α$ equivalent width grows sub-linearly with Fe abundance. Strikingly, the K$α$ profile is more sensitive to accretion rate than to spin because its radial surface brightness profile is relatively flat, and higher accretion rate extends the production region to smaller radii. The overall radiative efficiency is at least 30--100% greater than as predicted by the Novikov-Thorne model. △ Less

Submitted 19 May, 2021; originally announced May 2021.

Comments: 23 pages, 17 Figures, 1 Table

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

Electromagnetic Emission from a Binary Black Hole Merger Remnant in Plasma: Field Alignment and Plasma Temperature

Authors: Bernard J. Kelly, Zachariah B. Etienne, Jacob Golomb, Jeremy D. Schnittman, John G. Baker, Scott C. Noble, Geoffrey Ryan

Abstract: Comparable-mass black-hole mergers generically result in moderate to highly spinning holes, whose spacetime curvature will significantly affect nearby matter in observable ways. We investigate how the moderate spin of a post-merger Kerr black hole immersed in a plasma with initially uniform density and uniform magnetic field affects potentially observable accretion rates and energy fluxes. Varying… ▽ More Comparable-mass black-hole mergers generically result in moderate to highly spinning holes, whose spacetime curvature will significantly affect nearby matter in observable ways. We investigate how the moderate spin of a post-merger Kerr black hole immersed in a plasma with initially uniform density and uniform magnetic field affects potentially observable accretion rates and energy fluxes. Varying the initial specific internal energy of the plasma over two decades, we find very little change in steady-state mass accretion rate or Poynting luminosity, except at the lowest internal energies, where fluxes do not exhibit steady-state behavior during the simulation timescale. Fixing the internal energy and varying the initial fixed magnetic-field amplitude and orientation, we find that the steady-state Poynting luminosity depends strongly on the initial field angle with respect to the black hole spin axis, while the matter accretion rate is more stable until the field angle exceeds $\sim 45\degree$. The proto-jet formed along the black hole spin-axis conforms to a thin, elongated cylinder near the hole, while aligning with the asymptotic magnetic field at large distances. △ Less

Submitted 26 March, 2021; v1 submitted 21 October, 2020; originally announced October 2020.

Comments: 13 pages, 14 figures. Minor changes to match published version

Journal ref: Phys. Rev. D 103, 063039 (2021)

Inverse Compton Cooling in the Coronae of Simulated Black Hole Accretion Flows

Authors: Brooks E. Kinch, Scott C. Noble, Jeremy D. Schnittman, Julian H. Krolik

Abstract: We present a formulation for a local cooling function to be employed in the diffuse, hot corona region of 3D GRMHD simulations of accreting black holes. This new cooling function calculates the cooling rate due to inverse Compton scattering by considering the relevant microphysics in each cell in the corona and approximating the radiation energy density and Compton temperature there by integrating… ▽ More We present a formulation for a local cooling function to be employed in the diffuse, hot corona region of 3D GRMHD simulations of accreting black holes. This new cooling function calculates the cooling rate due to inverse Compton scattering by considering the relevant microphysics in each cell in the corona and approximating the radiation energy density and Compton temperature there by integrating over the thermal seed photon flux from the disk surface. The method either assumes ion and electron temperatures are equal (1T), or calculates them separately (2T) using an instantaneous equilibrium approach predicated on the actual relevant rate equations (Coulomb and Compton). The method is shown to be consistent with a more detailed ray-tracing calculation where the bulk of the cooling occurs, but is substantially less costly to perform. As an example, we apply these methods to a \textsc{harm3d} simulation of a $10 M_\odot$, non-spinning black hole, accreting at nominally 1\% the Eddington value. Both 1T and 2T approaches lead to increased radiative efficiency and a larger fraction of total cooling in the corona as compared to the original target-temperature cooling function used by \textsc{harm3d}, especially in the 1T case. Time-averaged post-processing reveals that the continuum spectral observations predicted from these simulations are qualitatively similar to actual X-ray binary data, especially so for the 1T approach which yields a harder power-law component ($Γ= 2.25$) compared to the 2T version ($Γ= 2.53$) △ Less

Submitted 3 September, 2020; originally announced September 2020.

Comments: 33 pages, 25 figures, submitted to ApJ

On the maximum stellar rotation to form a black hole without an accompanying luminous transient

Authors: Ariadna Murguia-Berthier, Aldo Batta, Agnieszka Janiuk, Enrico Ramirez-Ruiz, Ilya Mandel, Scott C. Noble, Rosa Wallace Everson

Abstract: The collapse of a massive star with low angular momentum content is commonly argued to result in the formation of a black hole without an accompanying bright transient. Our goal in this Letter is to understand the flow in and around a newly-formed black hole, involving accretion and rotation, via general relativistic hydrodynamics simulations aimed at studying the conditions under which infalling… ▽ More The collapse of a massive star with low angular momentum content is commonly argued to result in the formation of a black hole without an accompanying bright transient. Our goal in this Letter is to understand the flow in and around a newly-formed black hole, involving accretion and rotation, via general relativistic hydrodynamics simulations aimed at studying the conditions under which infalling material can accrete without forming a centrifugally supported structure and, as a result, generate no effective feedback. If the feedback from the black hole is, on the other hand, significant, the collapse would be halted and we suggest that the event is likely to be followed by a bright transient. We find that feedback is only efficient if the specific angular momentum of the infalling material at the innermost stable circular orbit exceeds that of geodesic circular flow at that radius by at least $\approx 20\%$. We use the results of our simulations to constrain the maximal stellar rotation rates of the disappearing massive progenitors PHL293B-LBV and N6946-BH1, and to provide an estimate of the overall rate of disappearing massive stars. We find that about a few percent of single O-type stars with measured rotational velocities are expected to spin below the critical value before collapse and are thus predicted to vanish without a trace. △ Less

Submitted 29 September, 2020; v1 submitted 20 May, 2020; originally announced May 2020.

Comments: 9 pages, 5 figures, accepted at ApJ Letters

Tidal disruptions of main sequence stars -- IV. Relativistic effects and dependence on black hole mass

Authors: Taeho Ryu, Julian Krolik, Tsvi Piran, Scott C. Noble

Abstract: Using a suite of fully relativistic hydrodynamic simulations applied to main-sequence stars with realistic internal density profiles, we examine full and partial tidal disruptions across a wide range of black hole mass ($10^{5}\leq M_{\rm BH}/\mathrm{M}_{\odot}\leq 5\times 10^{7}$) and stellar mass ($0.3 \leq M_{\star} /\mathrm{M}_{\odot}\leq 3$) as larger $M_{\rm BH}$ leads to stronger relativist… ▽ More Using a suite of fully relativistic hydrodynamic simulations applied to main-sequence stars with realistic internal density profiles, we examine full and partial tidal disruptions across a wide range of black hole mass ($10^{5}\leq M_{\rm BH}/\mathrm{M}_{\odot}\leq 5\times 10^{7}$) and stellar mass ($0.3 \leq M_{\star} /\mathrm{M}_{\odot}\leq 3$) as larger $M_{\rm BH}$ leads to stronger relativistic effects. For fixed $M_{\star}$, as $M_{\rm BH}$ increases, the ratio of the maximum pericenter distance yielding full disruptions ($\mathcal{R}_{\rm t}$) to its Newtonian prediction rises rapidly, becoming triple the Newtonian value for $M_{\rm BH} = 5\times10^{7}~{\rm M}_\odot$, while the ratio of the energy width of the stellar debris for full disruptions to the Newtonian prediction decreases steeply, resulting in a factor of two correction at $M_{\rm BH} = 5 \times 10^7~{\rm M}_\odot$. We find that for partial disruptions, the fractional remnant mass for a given ratio of the pericenter to $\mathcal{R}_{\rm t}$ is higher for larger $M_{\rm BH}$. These results have several implications. As $M_{\rm BH}$ increases above $\sim 10^7~{\rm M}_\odot$, the cross section for complete disruptions is suppressed by competition with direct capture. However, the cross section ratio for partial to complete disruptions depends only weakly on $M_{\rm BH}$. The relativistic correction to the debris energy width delays the time of peak mass-return rate and diminishes the magnitude of the peak return rate. For $M_{\rm BH} \gtrsim 10^7~{\rm M}_\odot$, the $M_{\rm BH}$-dependence of the full disruption cross section and the peak mass-return rate and time is influenced more by relativistic effects than by Newtonian dynamics. △ Less

Submitted 8 October, 2020; v1 submitted 8 January, 2020; originally announced January 2020.

Comments: 13 pages, 9 figures, 2 tables, accepted for publication in ApJ

Tidal disruptions of main sequence stars -- III. Stellar mass dependence of the character of partial disruptions

Authors: Taeho Ryu, Julian Krolik, Tsvi Piran, Scott C. Noble

Abstract: In this paper, the third in this series, we continue our study of tidal disruption events of main-sequence stars by a non-spinning $10^{6}~\rm{M}_\odot$ supermassive black hole. Here we focus on the stellar mass dependence of the outcomes of partial disruptions. As the encounter becomes weaker, the debris mass is increasingly concentrated near the outer edges of the energy distribution. As a resul… ▽ More In this paper, the third in this series, we continue our study of tidal disruption events of main-sequence stars by a non-spinning $10^{6}~\rm{M}_\odot$ supermassive black hole. Here we focus on the stellar mass dependence of the outcomes of partial disruptions. As the encounter becomes weaker, the debris mass is increasingly concentrated near the outer edges of the energy distribution. As a result, the mass fallback rate can deviate substantially from a $t^{-5/3}$ power-law, becoming more like a single peak with a tail declining as $t^{-p}$ with $p\simeq2-5$. Surviving remnants are spun-up in the prograde direction and are hotter than main sequence stars of the same mass. Their specific orbital energy is $\simeq10^{-3}\times$ that of the debris, but of either sign with respect to the black hole potential, while their specific angular momentum is close to that of the original star. Even for strong encounters, remnants have speeds at infinity relative to the black hole potential $\lesssim 300$ km s$^{-1}$, so they are unable to travel far out into the galactic bulge. The remnants most deeply bound to the black hole go through a second tidal disruption event upon their first return to pericenter; if they have not thermally relaxed, they will be completely disrupted. △ Less

Submitted 1 September, 2020; v1 submitted 8 January, 2020; originally announced January 2020.

Comments: 14 pages, 8 figures, accepted for publication in ApJ, Supersedes arXiv:1909.04041

Tidal disruptions of main sequence stars -- II. Simulation methodology and stellar mass dependence of the character of full tidal disruptions

Authors: Taeho Ryu, Julian Krolik, Tsvi Piran, Scott C. Noble

Abstract: This is the second in a series of papers presenting the results of fully general relativistic simulations of stellar tidal disruptions in which the stars' initial states are realistic main-sequence models. In the first paper (Paper I), we gave an overview of this program and discussed the principal observational implications of our work. Here we describe our calculational method and provide detail… ▽ More This is the second in a series of papers presenting the results of fully general relativistic simulations of stellar tidal disruptions in which the stars' initial states are realistic main-sequence models. In the first paper (Paper I), we gave an overview of this program and discussed the principal observational implications of our work. Here we describe our calculational method and provide details about the outcomes of full disruptions, focusing on the stellar mass dependence of the outcomes for a black hole of mass $10^{6}\rm{M}_{\odot}$. We consider eight different stellar masses, from $0.15~{\rm M}_\odot$ to $10~{\rm M}_\odot$. We find that, relative to the traditional order-of-magnitude estimate $r_{\rm t}$, the physical tidal radius of low-mass stars ($M_{\star} \lesssim 0.7~ {\rm M}_\odot$) is larger by tens of percent, while for high-mass stars ($M_{\star} \gtrsim1~ {\rm M}_\odot$) it is smaller by a factor 2--2.5. The traditional estimate of the range of energies found in the debris is $\approx 1.4\times$ too large for low-mass stars, but is a factor $\sim 2$ too small for high-mass stars; in addition, the energy distribution for high-mass stars has significant wings. For all stars undergoing tidal encounters, we find that mass-loss continues for many stellar vibration times because the black hole's tidal gravity competes with the instantaneous stellar gravity at the star's surface until the star has reached a distance from the black hole $\sim O(10)r_{\rm t}$. △ Less

Submitted 8 October, 2020; v1 submitted 8 January, 2020; originally announced January 2020.

Comments: 19 pages, 10 figures, 3 tables, accepted for publication in ApJ, Supersedes arXiv:1907.11883

Tidal Disruptions of Main Sequence Stars -- I. Observable Quantities and their Dependence on Stellar and Black Hole Mass

Authors: Taeho Ryu, Julian Krolik, Tsvi Piran, Scott C. Noble

Abstract: This paper introduces a series of papers presenting a quantitative theory for the tidal disruption of main sequence stars by supermassive black holes. Using fully general relativistic hydrodynamics simulations and MESA-model initial conditions, we explore the pericenter-dependence of tidal disruption properties for eight stellar masses ($0.15 \leq M_*/M_\odot \leq 10$) and six black hole masses (… ▽ More This paper introduces a series of papers presenting a quantitative theory for the tidal disruption of main sequence stars by supermassive black holes. Using fully general relativistic hydrodynamics simulations and MESA-model initial conditions, we explore the pericenter-dependence of tidal disruption properties for eight stellar masses ($0.15 \leq M_*/M_\odot \leq 10$) and six black hole masses ($10^5 \leq M_{BH}/M_\odot \leq 5 \times 10^7$). We present here the results most relevant to observations. The effects of internal stellar structure and relativity decouple for both the disruption cross section and the characteristic energy width of the debris. Moreover, the full disruption cross section is almost independent of $M_*$ for $M_*/M_\odot \lesssim 3$. Independent of $M_*$, relativistic effects increase the critical pericenter distance for full disruptions by up to a factor $\sim 3$ relative to the Newtonian prediction. The probability of a direct capture is also independent of $M_*$; at $M_{BH}/M_\odot \simeq 5 \times 10^6$ this probability is equal to that of a complete disruption. The width of the debris energy distribution $ΔE$ can differ from the standard estimate by factors from 0.35 to 2, depending on $M_*$ and $M_{BH}$, implying a corresponding change in the characteristic mass-return timescale. The "frozen-in approximation" is inconsistent with $ΔE$, and mass-loss continues over a long span of time. We provide analytic forms, suitable for use in both event rate estimates and parameter inference, to describe all these trends. For partial disruptions, we find a nearly-universal relation between the star's angular momentum and the fraction of $M_*$ remaining. Within the "empty loss-cone" regime, partial disruptions must precede full disruptions. These partial disruptions can drastically affect the rate and appearance of subsequent total disruptions. △ Less

Submitted 1 September, 2020; v1 submitted 8 January, 2020; originally announced January 2020.

Comments: 15 pages, 8 figures, accepted for publication in ApJ, Supersedes arXiv:1907.08205

Tidal radii of main sequence stars -- III. Partial disruptions

Authors: Taeho Ryu, Julian Krolik, Tsvi Piran, Scott C. Noble

Abstract: In this paper, the third in this series, we continue our study of tidal disruption events of main-sequence stars by a non-spinning $10^{6}~\rm{M}_\odot$ supermassive black hole. Here we focus on the outcomes of partial disruptions. As the encounter becomes weaker, the debris mass is increasingly concentrated near the outer edges of the energy distribution. As a result, the mass fallback rate can d… ▽ More In this paper, the third in this series, we continue our study of tidal disruption events of main-sequence stars by a non-spinning $10^{6}~\rm{M}_\odot$ supermassive black hole. Here we focus on the outcomes of partial disruptions. As the encounter becomes weaker, the debris mass is increasingly concentrated near the outer edges of the energy distribution. As a result, the mass fallback rate can deviate substantially from a $t^{-5/3}$ power-law, becoming more like a single peak with a tail declining as $t^{-p}$ with $p\simeq2-5$. Surviving remnants are spun-up in the prograde direction and are hotter than MS stars of the same mass. Their specific orbital energy is $\simeq10^{-3}\times$ that of the debris (but of either sign with respect to the black hole potential) while their specific angular momentum is close to that of the original star. Even for strong encounters, remnants have speeds at infinity relative to the black hole potential $\lesssim 300$ km s$^{-1}$, so they are unable to travel far out into the galactic bulge. Remnants bound to the black hole can possibly go through a second tidal disruption event. △ Less

Submitted 13 January, 2020; v1 submitted 9 September, 2019; originally announced September 2019.

Comments: We have significantly revised this paper as a response to the referee's report. So we submitted the revised version with a new title as a new paper (arXiv:2001.03503) without replacing the old version

Tidal radii of main sequence stars -- II. Simulation methodology and the character of full tidal disruptions

Authors: Taeho Ryu, Julian Krolik, Tsvi Piran, Scott Noble

Abstract: This is the second in a series of papers presenting the results of fully general relativistic simulations of stellar tidal disruptions in which the stars' initial states are realistic main-sequence models. We consider eight different stellar masses, from $0.15~{\rm M}_\odot$ to $10~{\rm M}_\odot$. In the first paper (Ryu et al. 2019a), we gave an overview of this program and discussed the principa… ▽ More This is the second in a series of papers presenting the results of fully general relativistic simulations of stellar tidal disruptions in which the stars' initial states are realistic main-sequence models. We consider eight different stellar masses, from $0.15~{\rm M}_\odot$ to $10~{\rm M}_\odot$. In the first paper (Ryu et al. 2019a), we gave an overview of this program and discussed the principal observational implications of our work. Here we describe our calculational method and provide details about the outcomes of full disruptions. We find that, relative to the traditional order-of-magnitude estimate $r_{\rm t}$, the physical tidal radius of low-mass stars is larger by tens of percent, while for high-mass stars ($M_{\star} \gtrsim1~ {\rm M}_\odot$) it is smaller by a factor $2-2.5$. The traditional estimate of the range of energies found in the debris is approximately accurate for low-mass stars, but is a factor $\sim 2$ too small for high-mass stars; in addition, the energy distribution for high-mass stars has significant wings. For all stars undergoing tidal encounters, we find that mass-loss continues for a long time because the ${\it instantaneous}$ tidal radius, the distance out to which the black hole's tidal gravity competes with the instantaneous stellar gravity at the star's surface, stays comparable to the distance to the black hole until the star has reached $O(10)~r_{\rm t}$. These findings indicate significant failings in the popular "frozen-in" approximation. △ Less

Submitted 13 January, 2020; v1 submitted 27 July, 2019; originally announced July 2019.

Comments: We significantly changed the structure of this paper (with a new title) as a response to the referee's report. We emphasize that our results have not been changed, but added more explanations for our numerical methods. We submitted the new version with the new title (arXiv:2001.03502)

GR-MHD disk winds and jets from black holes and resistive accretion disks

Authors: Christos Vourellis, Christian Fendt, Qian Qian, Scott Noble

Abstract: We perform GR-MHD simulations of outflow launching from thin accretion disks. As in the non-relativistic case, resistivity is essential for the mass loading of the disk wind. We implemented resistivity in the ideal GR-MHD code HARM3D, extending previous works (Qian et al. 2017, 2018) for larger physical grids, higher spatial resolution, and longer simulation time. We consider an initially thin, re… ▽ More We perform GR-MHD simulations of outflow launching from thin accretion disks. As in the non-relativistic case, resistivity is essential for the mass loading of the disk wind. We implemented resistivity in the ideal GR-MHD code HARM3D, extending previous works (Qian et al. 2017, 2018) for larger physical grids, higher spatial resolution, and longer simulation time. We consider an initially thin, resistive disk orbiting the black hole, threaded by a large-scale magnetic flux. As the system evolves, outflows are launched from the black hole magnetosphere and the disk surface. We mainly focus on disk outflows, investigating their MHD structure and energy output in comparison with the Poynting-dominated black hole jet. The disk wind encloses two components -- a fast component dominated by the toroidal magnetic field and a slower component dominated by the poloidal field. The disk wind transitions from sub to super-Alfvénic speed, reaching velocities $\simeq 0.1c$. We provide parameter studies varying spin parameter and resistivity level, and measure the respective mass and energy fluxes. A higher spin strengthens the $B_φ$-dominated disk wind along the inner jet. We disentangle a critical resistivity level that leads to a maximum matter and energy output for both, resulting from the interplay between re-connection and diffusion, which in combination govern the magnetic flux and the mass loading. For counter-rotating black holes the outflow structure shows a magnetic field reversal. We estimate the opacity of the inner-most accretion stream and the outflow structure around it. This stream may be critically opaque for a lensed signal, while the axial jet funnel remains optically thin. △ Less

Submitted 24 July, 2019; originally announced July 2019.

Comments: 27 pages, 31 figures, accepted for publication in ApJ

Tidal radii of main sequence stars -- I. Physical tidal radius, semi-analytic model and their implications

Authors: Taeho Ryu, Julian Krolik, Tsvi Piran, Scott Noble

Abstract: A star is tidally disrupted by a supermassive black hole when their separation is shorter than the "tidal radius". This quantity is often estimated on an order-of-magnitude basis without reference to the star's internal structure. Using MESA models for main sequence stars and fully general relativistic dynamics, we find the physical tidal radius for complete disruption $\cal{R}_t$ for a… ▽ More A star is tidally disrupted by a supermassive black hole when their separation is shorter than the "tidal radius". This quantity is often estimated on an order-of-magnitude basis without reference to the star's internal structure. Using MESA models for main sequence stars and fully general relativistic dynamics, we find the physical tidal radius for complete disruption $\cal{R}_t$ for a $10^6M_\odot$ black hole (BH). We find that across a factor $\sim20$ in stellar mass $M_*$, i.e., $0.15M_{\odot}\leq M_*\leq3M_\odot$, $\cal{R}_t\sim27\times$(BH's gravitational radius). When comparing $\cal{R}_t$ with the commonly used order-of-magnitude estimate $r_t$, we find that $\cal{R}_t\sim1.05-1.45r_t$ for $0.15M_\odot\leq M_*\leq0.5M_\odot$, but between $0.5 M_\odot$ and $1 M_\odot$, $\cal{R}_t$ drops to $\sim 0.45r_t$, and it remains at this value up to $10 M_\odot$. The near-constancy of $\cal{R}_t$ implies a weaker dependence of the full disruption rate on $M_*$ than when predicted with $r_t$. The characteristic energy width of the debris $ΔE$ ranges from $\sim1.2Δ\cal{E}$ for low-mass stars to $\sim 0.35Δ\cal{E}$ for higher-mass stars, where $Δ\cal{E}=GM_{\rm BH}R_*/\cal{R}_t^{2}$. We present analytic fits for the $M_*$ dependence of $\cal{R}_t$ and $ΔE$; these fits lead to analytic expressions for the time of peak mass fallback rate and the maximal mass fallback rate. Our results also bear on the fraction of events leading to fast or slow circularization, as well as on the character of the tidal event occurring when the remnant of a partial disruption returns to the black hole. Using a semi-analytic model, we show that $\cal{R}_t$ is primarily determined by the star's central density rather than its mean density. For high-mass stars, the full disruption rate is roughly 1/4 the partial disruption rate, while this ratio is close to unity for low-mass stars. △ Less

Submitted 13 January, 2020; v1 submitted 18 July, 2019; originally announced July 2019.

Comments: We significantly changed the structure of this paper (with a new title) enough to be considered as completely new papers as a response to the referee's report. We emphasize that our results have not been changed, but only new results (consistent with those presented in Papers I, II and III previously posted) have been added. We submitted the new version with the new title (arXiv:2001.03501)

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

The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project

Authors: Oliver Porth, Koushik Chatterjee, Ramesh Narayan, Charles F. Gammie, Yosuke Mizuno, Peter Anninos, John G. Baker, Matteo Bugli, Chi-kwan Chan, Jordy Davelaar, Luca Del Zanna, Zachariah B. Etienne, P. Chris Fragile, Bernard J. Kelly, Matthew Liska, Sera Markoff, Jonathan C. McKinney, Bhupendra Mishra, Scott C. Noble, Héctor Olivares, Ben Prather, Luciano Rezzolla, Benjamin R. Ryan, James M. Stone, Niccolò Tomei , et al. (3 additional authors not shown)

Abstract: Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic… ▽ More Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL. Agreement between the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems. △ Less

Submitted 5 August, 2019; v1 submitted 9 April, 2019; originally announced April 2019.

Comments: Accepted version for publication in ApJS (May 28th, 2019)

Multi-Messenger Astrophysics with Pulsar Timing Arrays

Authors: Luke Zoltan Kelley, Maria Charisi, Sarah Burke-Spolaor, Joseph Simon, Laura Blecha, Tamara Bogdanovic, Monica Colpi, Julie Comerford, Daniel J. D'Orazio, Massimo Dotti, Michael Eracleous, Matthew Graham, Jenny E. Greene, Zoltán Haiman, Kelly Holley-Bockelmann, Erin Kara, Bernard Kelly, S. Komossa, Shane L. Larson, Xin Liu, Chung-Pei Ma, Scott Noble, Vasileios Paschalidis, Roman R. Rafikov, Vikram Ravi , et al. (7 additional authors not shown)

Abstract: Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will… ▽ More Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts. △ Less

Submitted 18 March, 2019; originally announced March 2019.

Comments: Submitted to the Astro2020 decadal review. This is one of five core white papers written by members of the NANOGrav Collaboration, see also: J.Cordes et al., S.R.Taylor et al., X.Siemens et al., and E.Fonseca et al

X-ray follow-up of extragalactic transients

Authors: Erin Kara, Raffaella Margutti, Azadeh Keivani, Wen-fai Fong, Brad Cenko, Scott Noble, Richard Mushotzky, John Ruan, Geoffrey Ryan, Eric Burns, Daryl Haggard, Regina Caputo, Derek Fox, David Burrows

Abstract: Most violent and energetic processes in our universe, including mergers of compact objects, explosions of massive stars and extreme accretion events, produce copious amounts of X-rays. X-ray follow-up is an efficient tool for identifying transients because (1) X-rays can quickly localize transients with large error circles, and (2) X-rays reveal the nature of transients that may not have unique si… ▽ More Most violent and energetic processes in our universe, including mergers of compact objects, explosions of massive stars and extreme accretion events, produce copious amounts of X-rays. X-ray follow-up is an efficient tool for identifying transients because (1) X-rays can quickly localize transients with large error circles, and (2) X-rays reveal the nature of transients that may not have unique signatures at other wavelengths. In this white paper, we identify key science questions about several extragalactic multi-messenger and multi-wavelength transients, and demonstrate how X-ray follow-up helps answer these questions △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: White paper submitted to the Astro2020 Decadal Survey

Multimessenger science opportunities with mHz gravitational waves

Authors: John Baker, Zoltán Haiman, Elena Maria Rossi, Edo Berger, Niel Brandt, Elmé Breedt, Katelyn Breivik, Maria Charisi, Andrea Derdzinski, Daniel J. D'Orazio, Saavik Ford, Jenny E. Greene, J. Colin Hill, Kelly Holley-Bockelmann, Joey Shapiro Key, Bence Kocsis, Thomas Kupfer, Shane Larson, Piero Madau, Thomas Marsh, Barry McKernan, Sean T. McWilliams, Priyamvada Natarajan, Samaya Nissanke, Scott Noble , et al. (10 additional authors not shown)

Abstract: LISA will open the mHz band of gravitational waves (GWs) to the astronomy community. The strong gravity which powers the variety of GW sources in this band is also crucial in a number of important astrophysical processes at the current frontiers of astronomy. These range from the beginning of structure formation in the early universe, through the origin and cosmic evolution of massive black holes… ▽ More LISA will open the mHz band of gravitational waves (GWs) to the astronomy community. The strong gravity which powers the variety of GW sources in this band is also crucial in a number of important astrophysical processes at the current frontiers of astronomy. These range from the beginning of structure formation in the early universe, through the origin and cosmic evolution of massive black holes in concert with their galactic environments, to the evolution of stellar remnant binaries in the Milky Way and in nearby galaxies. These processes and their associated populations also drive current and future observations across the electromagnetic (EM) spectrum. We review opportunities for science breakthroughs, involving either direct coincident EM+GW observations, or indirect multimessenger studies. We argue that for the US community to fully capitalize on the opportunities from the LISA mission, the US efforts should be accompanied by a coordinated and sustained program of multi-disciplinary science investment, following the GW data through to its impact on broad areas of astrophysics. Support for LISA-related multimessenger observers and theorists should be sized appropriately for a flagship observatory and may be coordinated through a dedicated mHz GW research center. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: Submitted to the Astro2020 Decadal Survey call for science white papers

STROBE-X: X-ray Timing and Spectroscopy on Dynamical Timescales from Microseconds to Years

Authors: Paul S. Ray, Zaven Arzoumanian, David Ballantyne, Enrico Bozzo, Soren Brandt, Laura Brenneman, Deepto Chakrabarty, Marc Christophersen, Alessandra DeRosa, Marco Feroci, Keith Gendreau, Adam Goldstein, Dieter Hartmann, Margarita Hernanz, Peter Jenke, Erin Kara, Tom Maccarone, Michael McDonald, Michael Nowak, Bernard Phlips, Ron Remillard, Abigail Stevens, John Tomsick, Anna Watts, Colleen Wilson-Hodge , et al. (134 additional authors not shown)

Abstract: We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniqu… ▽ More We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science. △ Less

Submitted 8 March, 2019; v1 submitted 7 March, 2019; originally announced March 2019.

Comments: 50 pages, Probe class mission concept study report submitted to NASA for Astro2020 Decadal Survey

Detectability of modulated X-rays from LISA's supermassive black hole mergers

Authors: Tito Dal Canton, Alberto Mangiagli, Scott C. Noble, Jeremy Schnittman, Andrew Ptak, Antoine Klein, Alberto Sesana, Jordan Camp

Abstract: One of the central goals of LISA is the detection of gravitational waves from the merger of supermassive black holes. Contrary to stellar-mass black hole mergers, such events are expected to be rich X-ray sources due to the accretion of material from the circumbinary disks onto the black holes. The orbital dynamics before merger is also expected to modulate the resulting X-ray emission via Doppler… ▽ More One of the central goals of LISA is the detection of gravitational waves from the merger of supermassive black holes. Contrary to stellar-mass black hole mergers, such events are expected to be rich X-ray sources due to the accretion of material from the circumbinary disks onto the black holes. The orbital dynamics before merger is also expected to modulate the resulting X-ray emission via Doppler boosting in a quasi-periodic way, and in a simple phase relation with the gravitational wave from the inspiral of the black holes. Detecting the X-ray source would enable a precise and early localization of the binary, thus allowing many telescopes to observe the very moment of the merger. Although identifying the correct X-ray source in the relatively large LISA sky localization will be challenging due to the presence of many confounding point sources, the quasi-periodic modulation may aid in the identification. We explore the practical feasibility of such idea. We simulate populations of merging supermassive black holes, their detection with LISA and their X-ray lightcurves using a simple model. Taking the parameters of the X-ray Telescope on the proposed NASA Transient Astrophysics Probe, we then design and simulate an observation campaign that searches for the modulated X-ray source while LISA is still observing the inspiral of the black holes. Assuming a fiducial LISA detection rate of $10$ mergers per year at redshift closer than $3.5$, we expect a few detections of modulated X-ray counterparts over the nominal duration of the LISA mission. △ Less

Submitted 6 December, 2019; v1 submitted 4 February, 2019; originally announced February 2019.

Comments: 16 pages, 12 figures. Version accepted for publication on ApJ

Journal ref: ApJ 886 146 (2019)

Parabolic Jets from the Spinning Black Hole in M87

Authors: Masanori Nakamura, Keiichi Asada, Kazuhiro Hada, Hung-Yi Pu, Scott Noble, Chihyin Tseng, Kenji Toma, Motoki Kino, Hiroshi Nagai, Kazuya Takahashi, Juan-Carlos Algaba, Monica Orienti, Kazunori Akiyama, Akihiro Doi, Gabriele Giovannini, Marcello Giroletti, Mareki Honma, Shoko Koyama, Rocco Lico, Kotaro Niinuma, Fumie Tazaki

Abstract: The M87 jet is extensively examined by utilizing general relativistic magnetohydrodynamic (GRMHD) simulations as well as the steady axisymmetric force-free electrodynamic (FFE) solution. Quasi-steady funnel jets are obtained in GRMHD simulations up to the scale of $\sim 100$ gravitational radius ($r_{\rm g}$) for various black hole (BH) spins. As is known, the funnel edge is approximately determin… ▽ More The M87 jet is extensively examined by utilizing general relativistic magnetohydrodynamic (GRMHD) simulations as well as the steady axisymmetric force-free electrodynamic (FFE) solution. Quasi-steady funnel jets are obtained in GRMHD simulations up to the scale of $\sim 100$ gravitational radius ($r_{\rm g}$) for various black hole (BH) spins. As is known, the funnel edge is approximately determined by the following equipartitions; i) the magnetic and rest-mass energy densities and ii) the gas and magnetic pressures. Our numerical results give an additional factor that they follow the outermost parabolic streamline of the FFE solution, which is anchored to the event horizon on the equatorial plane. We also identify the matter dominated, non-relativistic corona/wind play a dynamical role in shaping the funnel jet into the parabolic geometry. We confirm a quantitative overlap between the outermost parabolic streamline of the FFE jet and the edge of jet sheath in VLBI observations at $\sim 10^{1}$-$10^{5} \, r_{\rm g}$, suggesting that the M87 jet is likely powered by the spinning BH. Our GRMHD simulations also indicate a lateral stratification of the bulk acceleration (i.e., the spine-sheath structure) as well as an emergence of knotty superluminal features. The spin characterizes the location of the jet stagnation surface inside the funnel. We suggest that the limb-brightened feature could be associated with the nature of the BH-driven jet, if the Doppler beaming is a dominant factor. Our findings can be examined with (sub-)mm VLBI observations, giving a clue for the origin of the M87 jet. △ Less

Submitted 23 October, 2018; originally announced October 2018.

Comments: 29 pages, 23 figures, accepted for publication in ApJ

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

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

PATCHWORK: A Multipatch Infrastructure for Multiphysics/Multiscale/Multiframe Fluid Simulations

Authors: Hotaka Shiokawa, Roseanne M. Cheng, Scott C. Noble, Julian H. Krolik

Abstract: We present a "multipatch" infrastructure for numerical simulation of fluid problems in which sub-regions require different gridscales, different grid geometries, different physical equations, or different reference frames. Its key element is a sophisticated client-router-server framework for efficiently linking processors supporting different regions ("patches") that must exchange boundary data. T… ▽ More We present a "multipatch" infrastructure for numerical simulation of fluid problems in which sub-regions require different gridscales, different grid geometries, different physical equations, or different reference frames. Its key element is a sophisticated client-router-server framework for efficiently linking processors supporting different regions ("patches") that must exchange boundary data. This infrastructure may be used with a wide variety of fluid dynamics codes; the only requirement is that their primary dependent variables be the same in all patches, e.g., fluid mass density, internal energy density, and velocity. Its structure can accommodate either Newtonian or relativistic dynamics. The overhead imposed by this system is both problem- and computer cluster architecture-dependent. Compared to a conventional simulation using the same number of cells and processors, the increase in runtime can be anywhere from negligible to a factor of a few; however, one of the infrastructure's advantages is that it can lead to a very large reduction in the total number of zone-updates. △ Less

Submitted 14 May, 2018; v1 submitted 19 January, 2017; originally announced January 2017.

Comments: 24 pages, 15 figures, accepted for publication in ApJ

Report number: LA-UR-18-22528

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

Scale-invariant radio jets and varying black hole spin

Authors: Monika Moscibrodzka, Heino Falcke, Scott Noble

Abstract: Compact radio cores associated with relativistic jets are often observed in both active galactic nuclei and X-ray binaries. Their radiative properties follow some general scaling laws which primarily depend on their masses and accretion rates. However, it has been suggested that the black hole spin can also strongly influence the power and radio flux of these. Here, we attempt to estimate the depe… ▽ More Compact radio cores associated with relativistic jets are often observed in both active galactic nuclei and X-ray binaries. Their radiative properties follow some general scaling laws which primarily depend on their masses and accretion rates. However, it has been suggested that the black hole spin can also strongly influence the power and radio flux of these. Here, we attempt to estimate the dependency of the radio luminosity of steady jets launched by accretion disks on black hole mass, accretion rate and spin using numerical simulations. We make use of 3D GRMHD simulations of accretion disks around low-luminosity black holes in which the jet radio emission is produced by the jet sheath. We find that the radio flux increases roughly by a factor of 6 as the back hole spin increases from a~0 to a=0.98. This is comparable to the increase in accretion power with spin, meaning that the ratio between radio jet and accretion power is hardly changing. Although our jet spine power scales as expected for the Blandford-Znajek process, the dependency of jet radio luminosity on the black hole spin is somewhat weaker. Also weakly rotating black holes can produce visible radio jets. The overall scaling of the radio emission with black hole mass and accretion rate is consistent with the scale-invariant analytical models used to explain the fundamental plane of black hole activity. Spin does not introduce a significant scatter in this model. The jet-sheath model can describe well resolved accreting systems, such as SgrA* and M87, as well as the general scaling behavior of low-luminosity black holes. Hence the model should be applicable to a wide range of radio jets in sub-Eddington black holes. The black hole spin has an effect on the production of visible radio jet, but it may not be the main driver to produce visible radio jets. An extension of our findings to powerful quasars remains speculative. △ Less

Submitted 27 October, 2016; originally announced October 2016.

Comments: 10 pages, 6 figures, A&A accepted

rHARM: Accretion and Ejection in Resistive GR-MHD

Authors: Qian Qian, Christian Fendt, Scott Noble, Matteo Bugli

Abstract: Turbulent magnetic diffusivity plays an important role for accretion disks and the launching of disk winds. We have implemented magnetic diffusivity, respective resistivity in the general relativistic MHD code HARM. This paper describes the theoretical background of our implementation, its numerical realization, our numerical tests and preliminary applications. The test simulations of the new code… ▽ More Turbulent magnetic diffusivity plays an important role for accretion disks and the launching of disk winds. We have implemented magnetic diffusivity, respective resistivity in the general relativistic MHD code HARM. This paper describes the theoretical background of our implementation, its numerical realization, our numerical tests and preliminary applications. The test simulations of the new code rHARM are compared with an analytic solution of the diffusion equation and a classical shock tube problem. We have further investigated the evolution of the magneto-rotational instability (MRI) in tori around black holes for a range of magnetic diffusivities. We find indication for a critical magnetic diffusivity (for our setup) beyond which no MRI develops in the linear regime and for which accretion of torus material to the black hole is delayed. Preliminary simulations of magnetically diffusive thin accretion disks around Schwarzschild black holes that are threaded by a large-scale poloidal magnetic field show the launching of disk winds with mass fluxes of about 50% of the accretion rate. The disk magnetic diffusivity allows for efficient disk accretion that replenishes the mass reservoir of the inner disk area and thus allows for long-term simulations of wind launching for more than 5000 time units. △ Less

Submitted 14 October, 2016; originally announced October 2016.

Comments: 21 pages, 43 figures, accepted by ApJ

Driven neutron star collapse: Type~I critical phenomena and the initial black hole mass distribution

Authors: Scott C. Noble, Matthew W. Choptuik

Abstract: We study the general relativistic collapse of neutron star (NS) models in spherical symmetry. Our initially stable models are driven to collapse by the addition of one of two things: an initially in-going velocity profile, or a shell of minimally coupled, massless scalar field that falls onto the star. Tolman-Oppenheimer-Volkoff (TOV) solutions with an initially isentropic, gamma-law equation of s… ▽ More We study the general relativistic collapse of neutron star (NS) models in spherical symmetry. Our initially stable models are driven to collapse by the addition of one of two things: an initially in-going velocity profile, or a shell of minimally coupled, massless scalar field that falls onto the star. Tolman-Oppenheimer-Volkoff (TOV) solutions with an initially isentropic, gamma-law equation of state serve as our NS models. The initial values of the velocity profile's amplitude and the star's central density span a parameter space which we have surveyed extensively and which we find provides a rich picture of the possible end states of NS collapse. This parameter space survey elucidates the boundary between Type I and Type II critical behavior in perfect fluids which coincides, on the subcritical side, with the boundary between dispersed and bound end states. For our particular model, initial velocity amplitudes greater than 0.3c are needed to probe the regime where arbitrarily small black holes can form. In addition, we investigate Type I behavior in our system by varying the initial amplitude of the initially imploding scalar field. In this case we find that the Type I critical solutions resemble TOV solutions on the 1-mode unstable branch of equilibrium solutions, and that the critical solutions' frequencies agree well with the fundamental mode frequencies of the unstable equilibria. Additionally, the critical solution's scaling exponent is shown to be well approximated by a linear function of the initial star's central density. △ Less

Submitted 9 December, 2015; originally announced December 2015.

Comments: Submitted to Phys. Rev. D., 24 pages, 25 monochrome figures. arXiv admin note: substantial text overlap with arXiv:gr-qc/0310116

Disk Emission from Magneto-hydrodynamic Simulations of Spinning Black Holes

Authors: Jeremy D. Schnittman, Julian H. Krolik, Scott C. Noble

Abstract: We present the results of a new series of global 3D relativistic magneto-hydrodynamic (MHD) simulations of thin accretion disks around spinning black holes. The disks have aspect ratios of $H/R\sim 0.05$ and spin parameters $a/M=0, 0.5, 0.9$, and $0.99$. Using the ray-tracing code Pandurata, we generate broad-band thermal spectra and polarization signatures from the MHD simulations. We find that t… ▽ More We present the results of a new series of global 3D relativistic magneto-hydrodynamic (MHD) simulations of thin accretion disks around spinning black holes. The disks have aspect ratios of $H/R\sim 0.05$ and spin parameters $a/M=0, 0.5, 0.9$, and $0.99$. Using the ray-tracing code Pandurata, we generate broad-band thermal spectra and polarization signatures from the MHD simulations. We find that the simulated spectra can be well fit with a simple, universal emissivity profile that better reproduces the behavior of the emission from the inner disk, compared to traditional analyses carried out using a Novikov-Thorne thin disk model. Lastly, we show how spectropolarization observations can be used to convincingly break the spin-inclination degeneracy well-known to the continuum fitting method of measuring black hole spin. △ Less

Submitted 2 December, 2015; originally announced December 2015.

Comments: submitted to ApJ

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)

General Relativistic Hydrodynamic Simulation of Accretion Flow from a Stellar Tidal Disruption

Authors: Hotaka Shiokawa, Julian H. Krolik, Roseanne M. Cheng, Tsvi Piran, Scott C. Noble

Abstract: We study how the matter dispersed when a supermassive black hole tidally disrupts a star joins an accretion flow. Combining a relativistic hydrodynamic simulation of the stellar disruption with a relativistic hydrodynamics simulation of the tidal debris motion, we track such a system until ~80% of the stellar mass bound to the black hole has settled into an accretion flow. Shocks near the stellar… ▽ More We study how the matter dispersed when a supermassive black hole tidally disrupts a star joins an accretion flow. Combining a relativistic hydrodynamic simulation of the stellar disruption with a relativistic hydrodynamics simulation of the tidal debris motion, we track such a system until ~80% of the stellar mass bound to the black hole has settled into an accretion flow. Shocks near the stellar pericenter and also near the apocenter of the most tightly-bound debris dissipate orbital energy, but only enough to make the characteristic radius comparable to the semi-major axis of the most-bound material, not the tidal radius as previously thought. The outer shocks are caused by post-Newtonian effects, both on the stellar orbit during its disruption and on the tidal forces. Accumulation of mass into the accretion flow is non-monotonic and slow, requiring ~3--10x the orbital period of the most tightly-bound tidal streams, while the inflow time for most of the mass may be comparable to or longer than the mass accumulation time. Deflection by shocks does, however, remove enough angular momentum and energy from some mass for it to move inward even before most of the mass is accumulated into the accretion flow. Although the accretion rate rises sharply and then decays roughly as a power-law, its maximum is ~0.1x the previous expectation, and the duration of the peak is ~5x longer than previously predicted. The geometric mean of the black hole mass and stellar mass inferred from a measured event timescale is therefore ~0.2x the value given by classical theory. △ Less

Submitted 18 January, 2015; originally announced January 2015.

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

Dynamic Fisheye Grids for Binary Black Holes Simulations

Authors: Miguel Zilhão, Scott C. Noble

Abstract: We present a new warped gridding scheme adapted to simulating gas dynamics in binary black hole spacetimes. The grid concentrates grid points in the vicinity of each black hole to resolve the smaller scale structures there, and rarefies grid points away from each black hole to keep the overall problem size at a practical level. In this respect, our system can be thought of as a "double" version of… ▽ More We present a new warped gridding scheme adapted to simulating gas dynamics in binary black hole spacetimes. The grid concentrates grid points in the vicinity of each black hole to resolve the smaller scale structures there, and rarefies grid points away from each black hole to keep the overall problem size at a practical level. In this respect, our system can be thought of as a "double" version of the fisheye coordinate system, used before in numerical relativity codes for evolving binary black holes. The gridding scheme is constructed as a mapping between a uniform coordinate system--in which the equations of motion are solved--to the distorted system representing the spatial locations of our grid points. Since we are motivated to eventually use this system for circumbinary disc calculations, we demonstrate how the distorted system can be constructed to asymptote to the typical spherical polar coordinate system, amenable to efficiently simulating orbiting gas flows about central objects with little numerical diffusion. We discuss its implementation in the Harm3d code, tailored to evolve the magnetohydrodynamics equations in curved spacetimes. We evaluate the performance of the system's implementation in Harm3d with a series of tests, such as the advected magnetic field loop test, magnetized Bondi accretion, and evolutions of hydrodynamic discs about a single black hole and about a binary black hole. Like we have done with Harm3d, this gridding scheme can be implemented in other unigrid codes as a (possibly) simpler alternative to adaptive mesh refinement. △ Less

Submitted 28 February, 2014; v1 submitted 11 September, 2013; originally announced September 2013.

Comments: 35 pages, 13 figures. Matches published version in CQG

Journal ref: 2014 Class. Quantum Grav. 31 065013

X-ray Spectra from MHD Simulations of Accreting Black Holes

Authors: Jeremy D. Schnittman, Julian H. Krolik, Scott C. Noble

Abstract: We present the results of a new global radiation transport code coupled to a general relativistic magneto-hydrodynamic simulation of an accreting, non-rotating black hole. For the first time, we are able to explain from first principles in a self-consistent way all the components seen in the X-ray spectra of stellar-mass black holes, including a thermal peak and all the features associated with st… ▽ More We present the results of a new global radiation transport code coupled to a general relativistic magneto-hydrodynamic simulation of an accreting, non-rotating black hole. For the first time, we are able to explain from first principles in a self-consistent way all the components seen in the X-ray spectra of stellar-mass black holes, including a thermal peak and all the features associated with strong hard X-ray emission: a power-law extending to high energies, a Compton reflection hump, and a broad iron line. Varying only the mass accretion rate, we are able to reproduce a wide range of X-ray states seen in most galactic black hole sources. The temperature in the corona is T_e ~ 10 keV in a boundary layer near the disk and rises smoothly to T_e >~ 100 keV in low-density regions far above the disk. Even as the disk's reflection edge varies from the horizon out to ~ 6M as the accretion rate decreases, we find that the shape of the Fe Kαline is remarkably constant. This is because photons emitted from the plunging region are strongly beamed into the horizon and never reach the observer. We have also carried out a basic timing analysis of the spectra and find that the fractional variability increases with photon energy and viewer inclination angle, consistent with the coronal hot spot model for X-ray fluctuations. △ Less

Submitted 30 April, 2013; v1 submitted 11 July, 2012; originally announced July 2012.

Comments: 22pp in emulateapj, v2: ApJ in press