-
Resolution analysis of magnetically arrested disk simulations
Authors:
Leon Sosapanta Salas,
Gibwa Musoke,
Koushik Chatterjee,
Sera Markoff,
Oliver Porth,
Matthew Liska,
Bart Ripperda
Abstract:
Polarisation measurements by the Event Horizon Telescope from M87$^{\ast}$ and Sgr A$^\ast$ suggest that there is a dynamically strong, ordered magnetic field, typical of what is expected of a magnetically arrested accretion disk (MAD). In such disks the strong poloidal magnetic field can suppress the accretion flow and cause episodic flux eruptions. Recent work shows that General Relativistic Mag…
▽ More
Polarisation measurements by the Event Horizon Telescope from M87$^{\ast}$ and Sgr A$^\ast$ suggest that there is a dynamically strong, ordered magnetic field, typical of what is expected of a magnetically arrested accretion disk (MAD). In such disks the strong poloidal magnetic field can suppress the accretion flow and cause episodic flux eruptions. Recent work shows that General Relativistic Magnetohydrodynamic (GRMHD) MAD simulations feature dynamics of turbulence and mixing instabilities that are becoming resolved at higher resolutions. We perform a convergence study of MADs exceeding the status quo by an order of magnitude in resolution. We use existing 3D simulations performed with the H-AMR code, up to resolution of 5376 x 2304 x 2304 in a logarithmic spherical-polar grid. We find consistent time-averaged disk properties across all resolutions. However, higher resolutions reveal signs of inward angular momentum transport attributed to turbulent convection, particularly evident when mixing instabilities occur at the surfaces of flux tubes during flux eruptions. Additionally, we see wave-like features in the jet sheath, which become more prominent at higher resolutions, that may induce mixing between jet and disk. At higher resolutions, we observe the sheath to be thinner, resulting in increased temperature, reduced magnetisation, and greater variability. Those differences could affect the dissipation of energy, that would eventually result in distinct observable radiative emission from high-resolution simulations. With higher resolutions, we can delve into crucial questions about horizon-scale physics and its impact on the dynamics and emission properties of larger-scale jets.
△ Less
Submitted 1 May, 2024;
originally announced May 2024.
-
Tilted Accretion Disks
Authors:
P. Chris Fragile,
Matthew Liska
Abstract:
In this chapter, we review some of the interesting consequences that tilt between the spin axis of the black hole and angular momentum axis of the accretion disk can have on the dynamics, thermodynamics, and observational appearance of accreting systems, from precessing coronae and jets to standing nozzle shocks and quasi-periodic oscillations. We begin the chapter by examining some of the reasons…
▽ More
In this chapter, we review some of the interesting consequences that tilt between the spin axis of the black hole and angular momentum axis of the accretion disk can have on the dynamics, thermodynamics, and observational appearance of accreting systems, from precessing coronae and jets to standing nozzle shocks and quasi-periodic oscillations. We begin the chapter by examining some of the reasons tilted disks are interesting as well as present arguments for how ubiquitous they may be. We then review the existing simulation results in the literature, broadly dividing them into sections on thick disks, thin disks, and magnetically arrested disks (MADs). We finish by highlighting some of the phenomenology that is unique to tilted disk simulations and discuss how this may apply to observations.
△ Less
Submitted 15 April, 2024;
originally announced April 2024.
-
Misaligned magnetized accretion flows onto spinning black holes: magneto-spin alignment, outflow power and intermittent jets
Authors:
Koushik Chatterjee,
Matthew Liska,
Alexander Tchekhovskoy,
Sera Markoff
Abstract:
Magnetic fields regulate black hole (BH) accretion, governing both the inflow and outflow dynamics. When a BH becomes saturated with large-scale vertical magnetic flux, it enters a magnetically-arrested disk (MAD) state. The dynamically-important BH magnetic flux powers highly efficient relativistic outflows (or jets) and sporadically erupts from the BH into the disk midplane. Here we explore the…
▽ More
Magnetic fields regulate black hole (BH) accretion, governing both the inflow and outflow dynamics. When a BH becomes saturated with large-scale vertical magnetic flux, it enters a magnetically-arrested disk (MAD) state. The dynamically-important BH magnetic flux powers highly efficient relativistic outflows (or jets) and sporadically erupts from the BH into the disk midplane. Here we explore the evolution of MADs when the BH and gas angular momentum are misaligned, which is expected to be more common. Using numerical simulations, we find that jets from rapidly spinning, prograde BHs force the inner accretion flow into alignment with the BH spin via the magneto-spin alignment mechanism for disks initially misaligned at $\mathcal{T}\lesssim 60^{\circ}$. Extremely misaligned MAD disks, on the other hand, exhibit intermittent jets that blow out parts of the disk to $\approx 100$ gravitational radii before collapsing, leaving behind hot cavities and magnetized filaments. These intermittent jet mechanism forms a mini-feedback cycle and could explain some cases of X-ray and radio quasi-periodic eruptions observed in dim AGN. Further, we find that (i) for BHs with low power jets, the BH spin and initial disk tilt angle changes the amount of horizon magnetic flux, and (ii) geometrically-thick, misaligned accretion flows do not undergo sustained Lense-Thirring (LT) precession. Thereby, we suggest that low-luminosity accreting BHs ($\dot{M}\ll 10^{-3} \dot{M}_{\rm Edd}$) are not likely to exhibit quasi-periodic oscillations in lightcurves due to LT precession, in agreement with observations of BH X-ray binaries and AGN in the low-hard/quiescent state. Instead, we suggest that magnetic flux eruptions can mimic precession-like motion, such as observed in the M87 jet, by driving large-scale surface waves in the jets.
△ Less
Submitted 1 November, 2023;
originally announced November 2023.
-
Magnetorotational dynamo can generate large-scale vertical magnetic fields in 3D GRMHD simulations of accreting black holes
Authors:
Jonatan Jacquemin-Ide,
François Rincon,
Alexander Tchekhovskoy,
Matthew Liska
Abstract:
Jetted astrophysical phenomena with black hole (BH) engines, including binary mergers, jetted tidal disruption events, and X-ray binaries, require a large-scale vertical magnetic field for efficient jet formation. However, a dynamo mechanism that could generate these crucial large-scale magnetic fields has not been identified and characterized. We have employed 3D global general relativistic magne…
▽ More
Jetted astrophysical phenomena with black hole (BH) engines, including binary mergers, jetted tidal disruption events, and X-ray binaries, require a large-scale vertical magnetic field for efficient jet formation. However, a dynamo mechanism that could generate these crucial large-scale magnetic fields has not been identified and characterized. We have employed 3D global general relativistic magnetohydrodynamical (MHD) simulations of accretion disks to quantify, for the first time, a dynamo mechanism that generates large-scale magnetic fields. This dynamo mechanism primarily arises from the nonlinear evolution of the magnetorotational instability (MRI). In this mechanism, large non-axisymmetric MRI-amplified shearing wave modes, mediated by the axisymmetric azimuthal magnetic field, generate and sustain the large-scale vertical magnetic field through their nonlinear interactions. We identify the advection of magnetic loops as a crucial feature, transporting the large-scale vertical magnetic field from the outer regions to the inner regions of the accretion disk. This leads to a larger characteristic size of the, now advected, magnetic field when compared to the local disk height. We characterize the complete dynamo mechanism with two timescales: one for the local magnetic field generation, $t_{\rm g}$, and one for the large-scale scale advection, $t_{\rm adv}$. Whereas the dynamo we describe is nonlinear, we explore the potential of linear mean field models to replicate its core features. Our findings indicate that traditional $α$-dynamo models, often computed in stratified shearing box simulations, are inadequate and that the effective large-scale dynamics is better described by the shear current effects or stochastic $α$-dynamos.
△ Less
Submitted 31 October, 2023;
originally announced November 2023.
-
How to Turn Jets into Cylinders near Supermassive Black Holes in 3D GRMHD Simulations
Authors:
Valeriia Rohoza,
Aretaios Lalakos,
Max Paik,
Koushik Chatterjee,
Matthew Liska,
Alexander Tchekhovskoy,
Ore Gottlieb
Abstract:
Accreting supermassive black holes (SMBHs) produce highly magnetized relativistic jets that tend to collimate gradually as they propagate outward. However, recent radio interferometric observations of the 3C 84 galaxy reveal a stunning, cylindrical jet already at several hundred SMBH gravitational radii, $r\gtrsim350r_{\rm g}$. We explore how such extreme collimation emerges via a suite of 3D gene…
▽ More
Accreting supermassive black holes (SMBHs) produce highly magnetized relativistic jets that tend to collimate gradually as they propagate outward. However, recent radio interferometric observations of the 3C 84 galaxy reveal a stunning, cylindrical jet already at several hundred SMBH gravitational radii, $r\gtrsim350r_{\rm g}$. We explore how such extreme collimation emerges via a suite of 3D general-relativistic magnetohydrodynamic (GRMHD) simulations. We consider an SMBH surrounded by a magnetized torus immersed in a constant-density ambient medium that starts at the edge of the SMBH sphere of influence, chosen to be much larger than the SMBH gravitational radius, $r_{\text{B}}=10^3r_{\text{g}}$. We find that radiatively inefficient accretion flows (e.g., M87) produce winds that collimate the jets into parabolas near the BH. After the disk winds stop collimating the jets at $r\lesssim{}r_\text{B}$, they turn conical. Once outside $r_\text{B}$, the jets run into the ambient medium and form backflows that collimate the jets into cylinders some distance beyond $r_{\text{B}}$. Interestingly, for radiatively-efficient accretion, as in 3C 84, the radiative cooling saps the energy out of the disk winds: at early times, they cannot efficiently collimate the jets, which skip the initial parabolic collimation stage, start out conical near the SMBH, and turn into cylinders already at $r\simeq300r_{\rm g}$, as observed in 3C 84. Over time, jet power remains approximately constant, whereas the mass accretion rate increases: the winds grow in strength and start to collimate the jets, which become quasi-parabolic near the base; the transition point to a nearly cylindrical jet profile moves outward while remaining inside $r_\text{B}$.
△ Less
Submitted 31 October, 2023;
originally announced November 2023.
-
Winds and Disk Turbulence Exert Equal Torques on Thick Magnetically Arrested Disks
Authors:
Vikram Manikantan,
Nicholas Kaaz,
Jonatan Jacquemin-Ide,
Gibwa Musoke,
Koushik Chatterjee,
Matthew Liska,
Alexander Tchekhovskoy
Abstract:
The conventional accretion disk lore is that magnetized turbulence is the principal angular momentum transport process that drives accretion. However, when dynamically important large-scale magnetic fields thread an accretion disk, they can produce mass and angular momentum outflows, known as winds, that also drive accretion. Yet, the relative importance of turbulent and wind-driven angular moment…
▽ More
The conventional accretion disk lore is that magnetized turbulence is the principal angular momentum transport process that drives accretion. However, when dynamically important large-scale magnetic fields thread an accretion disk, they can produce mass and angular momentum outflows, known as winds, that also drive accretion. Yet, the relative importance of turbulent and wind-driven angular momentum transport is still poorly understood. To probe this question, we analyze a long-duration ($1.2 \times 10^5 r_{\rm g}/c$) simulation of a rapidly rotating ($a=0.9$) black hole feeding from a thick ($H/r\sim0.3$), adiabatic, magnetically arrested disk (MAD), whose dynamically important magnetic field regulates mass inflow and drives both uncollimated and collimated outflows (i.e., winds and jets, respectively). By carefully disentangling the various angular momentum transport processes within the system, we demonstrate the novel result that disk winds and disk turbulence both extract roughly equal amounts of angular momentum from the disk. We find cumulative angular momentum and mass accretion outflow rates of $\dot{L}\propto r^{0.9}$ and $\dot{M}\propto r^{0.4}$, respectively. This result suggests that understanding both turbulent and laminar stresses is key to understanding the evolution of systems where geometrically thick MADs can occur, such as the hard state of X-ray binaries, low-luminosity active galactic nuclei, some tidal disruption events, and possibly gamma-ray bursts.
△ Less
Submitted 23 April, 2024; v1 submitted 17 October, 2023;
originally announced October 2023.
-
Jets with a Twist: Emergence of FR0 Jets in 3D GRMHD Simulation of Zero Angular Momentum Black Hole Accretion
Authors:
Aretaios Lalakos,
Alexander Tchekhovskoy,
Omer Bromberg,
Ore Gottlieb,
Jonatan Jacquemin-Ide,
Matthew Liska,
Haocheng Zhang
Abstract:
Spinning supermassive black holes (BHs) in active galactic nuclei (AGN) magnetically launch relativistic collimated outflows, or jets. Without angular momentum supply, such jets are thought to perish within $3$ orders of magnitude in distance from the BH, well before reaching kpc-scales. We study the survival of such jets at the largest scale separation to date, via 3D general relativistic magneto…
▽ More
Spinning supermassive black holes (BHs) in active galactic nuclei (AGN) magnetically launch relativistic collimated outflows, or jets. Without angular momentum supply, such jets are thought to perish within $3$ orders of magnitude in distance from the BH, well before reaching kpc-scales. We study the survival of such jets at the largest scale separation to date, via 3D general relativistic magnetohydrodynamic simulations of rapidly spinning BHs immersed into uniform zero-angular-momentum gas threaded by weak vertical magnetic field. We place the gas outside the BH sphere of influence, or the Bondi radius, chosen much larger than the BH gravitational radius, $R_\text{B}=10^3R_\text{g}$. The BH develops dynamically-important large-scale magnetic fields, forms a magnetically-arrested disk (MAD), and launches relativistic jets that propagate well outside $R_\text{B}$ and suppress BH accretion to $1.5\%$ of the Bondi rate, $\dot{M}_\text{B}$. Thus, low-angular-momentum accretion in the MAD state can form large-scale jets in Fanaroff-Riley (FR) type I and II galaxies. Subsequently, the disk shrinks and exits the MAD state: barely a disk (BAD), it rapidly precesses, whips the jets around, globally destroys them, and lets $5-10\%$ of $\dot{M}_\text{B}$ reach the BH. Thereafter, the disk starts rocking back and forth by angles $90-180^\circ$: the rocking accretion disk (RAD) launches weak intermittent jets that spread their energy over a large area and suppress BH accretion to $\lesssim 2 \% ~ \dot{M}_\text{B}$. Because BAD and RAD states tangle up the jets and destroy them well inside $R_\text{B}$, they are promising candidates for the more abundant, but less luminous, class of FR0 galaxies.
△ Less
Submitted 17 October, 2023;
originally announced October 2023.
-
Magnetic flux plays an important role during a BHXRB outburst in radiative 2T GRMHD simulations
Authors:
M. T. P. Liska,
N. Kaaz,
K. Chatterjee,
Razieh Emami,
Gibwa Musoke
Abstract:
Black hole (BH) X-ray binaries cycle through different spectral states of accretion over the course of months to years. Although fluctuations in the BH mass accretion rate are generally recognized as the most important component of state transitions, it is becoming increasingly evident that magnetic fields play a similarly important role. In this article, we present the first radiative two-tempera…
▽ More
Black hole (BH) X-ray binaries cycle through different spectral states of accretion over the course of months to years. Although fluctuations in the BH mass accretion rate are generally recognized as the most important component of state transitions, it is becoming increasingly evident that magnetic fields play a similarly important role. In this article, we present the first radiative two-temperature (2T) general relativistic magnetohydrodynamics (GRMHD) simulations in which an accretion disk transitions from a quiescent state at an accretion rate of $\dot{M} \sim 10^{-10} \dot{M}_{\rm Edd}$ to a hard-intermediate state at an accretion rate of $\dot{M} \sim 10^{-2} \dot{M}_{\rm Edd}$. This huge parameter space in mass accretion rate is bridged by artificially rescaling the gas density scale of the simulations. We present two jetted BH models with varying degrees of magnetic flux saturation. We demonstrate that in `Standard and Normal Evolution' models, which are unsaturated with magnetic flux, the hot torus collapses into a thin and cold accretion disk when $\dot{M} \gtrsim 5\times 10^{-3} \dot{M}_{\rm Edd}$. On the other hand, in `Magnetically Arrested Disk' models, which are fully saturated with vertical magnetic flux, the plasma remains mostly hot with substructures that condense into cold clumps of gas when $\dot{M} \gtrsim 1 \times 10^{-2} \dot{M}_{\rm Edd}$. This suggests that the spectral signatures observed during state transitions are closely tied to the level of magnetic flux saturation.
△ Less
Submitted 27 March, 2024; v1 submitted 27 September, 2023;
originally announced September 2023.
-
Synchrotron polarization signatures of surface waves in supermassive black hole jets
Authors:
Jordy Davelaar,
Bart Ripperda,
Lorenzo Sironi,
Alexander A. Philippov,
Hector Olivares,
Oliver Porth,
Bram van den Berg,
Thomas Bronzwaer,
Koushik Chatterjee,
Matthew Liska
Abstract:
Supermassive black holes in active galactic nuclei (AGN) are known to launch relativistic jets, which are observed across the entire electromagnetic spectrum and are thought to be efficient particle accelerators. Their primary radiation mechanism for radio emission is polarized synchrotron emission produced by a population of non-thermal electrons. In this Letter, we present a global general relat…
▽ More
Supermassive black holes in active galactic nuclei (AGN) are known to launch relativistic jets, which are observed across the entire electromagnetic spectrum and are thought to be efficient particle accelerators. Their primary radiation mechanism for radio emission is polarized synchrotron emission produced by a population of non-thermal electrons. In this Letter, we present a global general relativistic magnetohydrodynamical (GRMHD) simulation of a magnetically arrested disk (MAD). After the simulation reaches the MAD state, we show that waves are continuously launched from the vicinity of the black hole and propagate along the interface between the jet and the wind. At this interface, a steep gradient in velocity is present between the mildly relativistic wind and the highly relativistic jet. The interface is, therefore, a shear layer, and due to the shear, the waves generate roll-ups that alter the magnetic field configuration and the shear layer geometry. We then perform polarized radiation transfer calculations of our GRMHD simulation and find signatures of the waves in both total intensity and linear polarization, effectively lowering the fully resolved polarization fraction. The tell-tale polarization signatures of the waves could be observable by future Very Long Baseline Interferometric observations, e.g., by the next-generation Event Horizon Telescope.
△ Less
Submitted 24 October, 2023; v1 submitted 14 September, 2023;
originally announced September 2023.
-
Large-scale Evolution of Seconds-long Relativistic Jets from Black Hole-Neutron Star Mergers
Authors:
Ore Gottlieb,
Danat Issa,
Jonatan Jacquemin-Ide,
Matthew Liska,
Francois Foucart,
Alexander Tchekhovskoy,
Brian D. Metzger,
Eliot Quataert,
Rosalba Perna,
Daniel Kasen,
Matthew D. Duez,
Lawrence E. Kidder,
Harald P. Pfeiffer,
Mark A. Scheel
Abstract:
We present the first numerical simulations that track the evolution of a black hole-neutron star (BH-NS) merger from pre-merger to $r\gtrsim10^{11}\,{\rm cm}$. The disk that forms after a merger of mass ratio $q=2$ ejects massive disk winds ($3-5\times10^{-2}\,M_{\odot}$). We introduce various post-merger magnetic configurations and find that initial poloidal fields lead to jet launching shortly a…
▽ More
We present the first numerical simulations that track the evolution of a black hole-neutron star (BH-NS) merger from pre-merger to $r\gtrsim10^{11}\,{\rm cm}$. The disk that forms after a merger of mass ratio $q=2$ ejects massive disk winds ($3-5\times10^{-2}\,M_{\odot}$). We introduce various post-merger magnetic configurations and find that initial poloidal fields lead to jet launching shortly after the merger. The jet maintains a constant power due to the constancy of the large-scale BH magnetic flux until the disk becomes magnetically arrested (MAD), where the jet power falls off as $L_j\sim t^{-2}$. All jets inevitably exhibit either excessive luminosity due to rapid MAD activation when the accretion rate is high or excessive duration due to delayed MAD activation compared to typical short gamma-ray bursts (sGRBs). This provides a natural explanation for long sGRBs such as GRB 211211A but also raises a fundamental challenge to our understanding of jet formation in binary mergers. One possible implication is the necessity of higher binary mass ratios or moderate BH spins to launch typical sGRB jets. For post-merger disks with a toroidal magnetic field, dynamo processes delay jet launching such that the jets break out of the disk winds after several seconds. We show for the first time that sGRB jets with initial magnetization $σ_0>100$ retain significant magnetization ($σ\gg1$) at $r>10^{10}\,{\rm cm}$, emphasizing the importance of magnetic processes in the prompt emission. The jet-wind interaction leads to a power-law angular energy distribution by inflating an energetic cocoon whose emission is studied in a companion paper.
△ Less
Submitted 18 August, 2023; v1 submitted 26 June, 2023;
originally announced June 2023.
-
Hours-long Near-UV/Optical Emission from Mildly Relativistic Outflows in Black Hole-Neutron Star Mergers
Authors:
Ore Gottlieb,
Danat Issa,
Jonatan Jacquemin-Ide,
Matthew Liska,
Alexander Tchekhovskoy,
Francois Foucart,
Daniel Kasen,
Rosalba Perna,
Eliot Quataert,
Brian D. Metzger
Abstract:
The ongoing LIGO-Virgo-KAGRA observing run O4 provides an opportunity to discover new multi-messenger events, including binary neutron star (BNS) mergers such as GW170817, and the highly anticipated first detection of a multi-messenger black hole-neutron star (BH-NS) merger. While BNS mergers were predicted to exhibit early optical emission from mildly relativistic outflows, it has remained uncert…
▽ More
The ongoing LIGO-Virgo-KAGRA observing run O4 provides an opportunity to discover new multi-messenger events, including binary neutron star (BNS) mergers such as GW170817, and the highly anticipated first detection of a multi-messenger black hole-neutron star (BH-NS) merger. While BNS mergers were predicted to exhibit early optical emission from mildly relativistic outflows, it has remained uncertain whether the BH-NS merger ejecta provides the conditions for similar signals to emerge. We present the first modeling of early near-ultraviolet/optical emission from mildly relativistic outflows in BH-NS mergers. Adopting optimal binary properties: a mass ratio of $q=2$ and a rapidly rotating BH, we utilize numerical relativity and general relativistic magnetohydrodynamic (GRMHD) simulations to follow the binary's evolution from pre-merger to homologous expansion. We use an M1 neutrino transport GRMHD simulation to self-consistently estimate the opacity distribution in the outflows and find a bright near-ultraviolet/optical signal that emerges due to jet-powered cocoon cooling emission, outshining the kilonova emission at early time. The signal peaks at an absolute magnitude of $\sim -15$ a few hours after the merger, longer than previous estimates, which did not consider the first principles-based jet launching. By late 2024, the Rubin Observatory will have the capability to track the entire signal evolution or detect its peak up to distances of $\gtrsim1$ Gpc. In 2026, ULTRASAT will conduct all-sky surveys within minutes, detecting some of these events within $\sim 200$ Mpc. The BH-NS mergers with higher mass ratios or lower BH spins would produce shorter and fainter signals.
△ Less
Submitted 8 August, 2023; v1 submitted 26 June, 2023;
originally announced June 2023.
-
The EB-correlation in Resolved Polarized Images: Connections to Astrophysics of Black Holes
Authors:
Razieh Emami,
Sheperd S. Doeleman,
Maciek Wielgus,
Dominic Chang,
Koushik Chatterjee,
Randall Smith,
Matthew Liska,
James F. Steiner,
Angelo Ricarte,
Ramesh Narayan,
Grant Tremblay,
Douglas Finkbeiner,
Lars Hernquist,
Chi-Kwan Chan,
Lindy Blackburn,
Ben S. Prather,
Paul Tiede,
Avery E. Broderick,
Mark Vogelsberger,
Charles Alcock,
Freek Roelofs
Abstract:
We present an in-depth analysis of a newly proposed correlation function in visibility space, between the E and B modes of the linear polarization, hereafter the EB-correlation, for a set of time-averaged GRMHD simulations compared with the phase map from different semi-analytic models as well as the Event Horizon Telescope (EHT) 2017 data for M87* source. We demonstrate that the phase map of the…
▽ More
We present an in-depth analysis of a newly proposed correlation function in visibility space, between the E and B modes of the linear polarization, hereafter the EB-correlation, for a set of time-averaged GRMHD simulations compared with the phase map from different semi-analytic models as well as the Event Horizon Telescope (EHT) 2017 data for M87* source. We demonstrate that the phase map of the time-averaged EB-correlation contains novel information that might be linked to the BH spin, accretion state and the electron temperature. A detailed comparison with a semi-analytic approach with different azimuthal expansion modes shows that to recover the morphology of the real/imaginary part of the correlation function and its phase, we require higher orders of these azimuthal modes. To extract the phase features, we propose to use the Zernike polynomial reconstruction developing an empirical metric to break degeneracies between models with different BH spins that are qualitatively similar. We use a set of different geometrical ring models with various magnetic and velocity field morphologies and show that both the image space and visibility based EB-correlation morphologies in MAD simulations can be explained with simple fluid and magnetic field geometries as used in ring models. SANEs by contrast are harder to model, demonstrating that the simple fluid and magnetic field geometries of ring models are not sufficient to describe them owing to higher Faraday Rotation depths. A qualitative comparison with the EHT data demonstrates that some of the features in the phase of EB-correlation might be well explained by the current models for BH spins as well as electron temperatures, while others may require a larger theoretical surveys.
△ Less
Submitted 3 May, 2023; v1 submitted 30 April, 2023;
originally announced May 2023.
-
Millimeter Observational Signatures of Flares in Magnetically Arrested Black Hole Accretion Models
Authors:
He Jia,
Bart Ripperda,
Eliot Quataert,
Christopher J. White,
Koushik Chatterjee,
Alexander Philippov,
Matthew Liska
Abstract:
In general relativistic magneto-hydrodynamic (GRMHD) simulations, accreted magnetic flux on the black hole horizon episodically decays, during which magnetic reconnection heats up the plasma near the horizon, potentially powering high-energy flares like those observed in M87* and Sgr A*. We study the mm observational counterparts of such flaring episodes. The change in 230 GHz flux during the expe…
▽ More
In general relativistic magneto-hydrodynamic (GRMHD) simulations, accreted magnetic flux on the black hole horizon episodically decays, during which magnetic reconnection heats up the plasma near the horizon, potentially powering high-energy flares like those observed in M87* and Sgr A*. We study the mm observational counterparts of such flaring episodes. The change in 230 GHz flux during the expected high energy flares depends primarily on the efficiency of accelerating $γ\gtrsim 100$ ($T_e \gtrsim 10^{11}$ K) electrons. For models in which the electrons are heated to $T_e \sim 10^{11}$ K during flares, the hot plasma produced by reconnection significantly enhances 230 GHz emission and increases the size of the 230 GHz image. By contrast, for models in which the electrons are heated to higher temperatures (which we argue are better motivated), the reconnection-heated plasma is too hot to produce significant 230 GHz synchrotron emission, and the 230 GHz flux decreases during high energy flares. We do not find a significant change in the mm polarization during flares as long as the emission is Faraday thin. We also present expectations for the ring-shaped image as observed by the Event Horizon Telescope during flares, as well as multi-wavelength synchrotron spectra. Our results highlight several limitations of standard post-processing prescriptions for the electron temperature in GRMHD simulations. We also discuss the implications of our results for current and future observations of flares in Sgr A*, M87*, and related systems. Appendices contain detailed convergence studies with respect to resolution and plasma magnetization.
△ Less
Submitted 23 September, 2023; v1 submitted 21 January, 2023;
originally announced January 2023.
-
Probing plasma composition with the next generation Event Horizon Telescope (ngEHT)
Authors:
Razieh Emami,
Richard Anantua,
Angelo Ricarte,
Sheperd S. Doeleman,
Avery Broderick,
George Wong,
Lindy Blackburn,
Maciek Wielgus,
Ramesh Narayan,
Grant Tremblay,
Charles Alcock,
Lars Hernquist,
Randall Smith,
Matthew Liska,
Priyamvada Natarajan,
Mark Vogelsberger,
Brandon Curd,
Joana A. Kramer
Abstract:
We explore the plasma matter content in the innermost accretion disk/jet in M87* as relevant for an enthusiastic search for the signatures of anti-matter in the next generation of the Event Horizon Telescope (ngEHT). We model the impact of non-zero positron-to-electron ratio using different emission models including a constant electron to magnetic pressure (constant $β_e$ model) with a population…
▽ More
We explore the plasma matter content in the innermost accretion disk/jet in M87* as relevant for an enthusiastic search for the signatures of anti-matter in the next generation of the Event Horizon Telescope (ngEHT). We model the impact of non-zero positron-to-electron ratio using different emission models including a constant electron to magnetic pressure (constant $β_e$ model) with a population of non-thermal electrons as well as a R-beta model populated with thermal electrons. In the former case, we pick a semi-analytic fit to the force-free region of a general relativistic magnetohydrodynamic (GRMHD) simulation, while in the latter case, we analyze the GRMHD simulations directly. In both cases, positrons are being added at the post-processing level. We generate polarized images and spectra for some of these models and find out that at the radio frequencies, both of the linear and the circular polarizations get enhanced per adding pairs. On the contrary, we show that at higher frequencies a substantial positron fraction washes out the circular polarization. We report strong degeneracies between different emission models and the positron fraction, though our non-thermal models show more sensitivities to the pair fraction than the thermal models. We conclude that a large theoretical image library is indeed required to fully understand the trends probed in this study, and to place them in the context of large set of parameters which also affect polarimetric images, such as magnetic field strength, black hole spin, and detailed aspects of the electron temperature and the distribution function.
△ Less
Submitted 14 November, 2022;
originally announced November 2022.
-
Tracing the hot spot motion using the next generation Event Horizon Telescope (ngEHT)
Authors:
Razieh Emami,
Paul Tiede,
Sheperd S. Doeleman,
Freek Roelofs,
Maciek Wielgus,
Lindy Blackburn,
Matthew Liska,
Koushik Chatterjee,
Bart Ripperda,
Antonio Fuentes,
Avery Broderick,
Lars Hernquist,
Charles Alcock,
Ramesh Narayan,
Randall Smith,
Grant Tremblay,
Angelo Ricarte,
He Sun,
Richard Anantua,
Yuri Y. Kovalev,
Priyamvada Natarajan,
Mark Vogelsberger
Abstract:
We propose to trace the dynamical motion of a shearing hot spot near the SgrA* source through a dynamical image reconstruction algorithm, StarWarps. Such a hot spot may form as the exhaust of magnetic reconnection in a current sheet near the black hole horizon. A hot spot that is ejected from the current sheet into an orbit in the accretion disk may shear and diffuse due to instabilities at its bo…
▽ More
We propose to trace the dynamical motion of a shearing hot spot near the SgrA* source through a dynamical image reconstruction algorithm, StarWarps. Such a hot spot may form as the exhaust of magnetic reconnection in a current sheet near the black hole horizon. A hot spot that is ejected from the current sheet into an orbit in the accretion disk may shear and diffuse due to instabilities at its boundary during its orbit, resulting in a distinct signature. We subdivide the motion to two distinct phases; the first phase refers to the appearance of the hot spot modelled as a bright blob, followed by a subsequent shearing phase simulated as a stretched ellipse. We employ different observational arrays, including EHT(2017,2022) and the next generation event horizon telescope (ngEHTp1, ngEHT) arrays, in which few new additional sites are added to the observational array. We make dynamical image reconstructions for each of these arrays. Subsequently, we infer the hot spot phase in the first phase followed by the axes ratio and the ellipse area in the second phase. We focus on the direct observability of the orbiting hot spot in the sub-mm wavelength. Our analysis demonstrates that newly added dishes may easily trace the first phase as well as part of the second phase, before the flux is reduced substantially. The algorithm used in this work can be extended to any other types of the dynamical motion. Consequently, we conclude that the ngEHT is a key to directly observe the dynamical motions near variable sources, such as SgrA*.
△ Less
Submitted 15 November, 2022; v1 submitted 12 November, 2022;
originally announced November 2022.
-
Radiation Transport Two-Temperature GRMHD Simulations of Warped Accretion Disks
Authors:
M. T. P. Liska,
N. Kaaz,
G. Musoke,
A. Tchekhovskoy,
O. Porth
Abstract:
In many black hole systems, the accretion disk is expected to be misaligned with respect to the black hole spin axis. If the scale height of the disk is much smaller than the misalignment angle, the spin of the black hole can tear the disk into multiple, independently precessing `sub-disks'. This is most likely to happen during outbursts in black hole X-Ray binaries (BHXRBs) and in active galactic…
▽ More
In many black hole systems, the accretion disk is expected to be misaligned with respect to the black hole spin axis. If the scale height of the disk is much smaller than the misalignment angle, the spin of the black hole can tear the disk into multiple, independently precessing `sub-disks'. This is most likely to happen during outbursts in black hole X-Ray binaries (BHXRBs) and in active galactic nuclei (AGN) accreting above a few percent of the Eddington limit, because the disk becomes razor-thin. Disk tearing has the potential to explain variability phenomena including quasi-periodic oscillations (QPOs) in BHXRBs and changing-look phenomena in AGN. Here, we present the first radiative two-temperature GRMHD simulation of a strongly tilted ($65^{\circ}$) accretion disk around a $M_{BH}=10M_{\odot}$ black hole, which tears and precesses. This leads to luminosity swings between a few percent and $50 \%$ of the Eddington limit on sub-viscous timescales. Surprisingly, even where the disk is radiation pressure dominated, the accretion disk is thermally stable over $t \gtrsim 21,000 r_g/c$. This suggests warps play an important role in stabilizing the disk against thermal collapse. The disk forms two nozzle shocks perpendicular to the line of nodes where the scale height of the disk decreases $10$-fold and the electron temperature reaches $T_e \sim 10^8-10^9 K$. In addition, optically thin gas crossing the tear between the inner and outer disk gets heated to $T_e \sim 10^8 K$. This suggests that warped disks may emit a Comptonized spectrum that deviates substantially from idealized models.
△ Less
Submitted 14 July, 2023; v1 submitted 18 October, 2022;
originally announced October 2022.
-
Nozzle Shocks, Disk Tearing and Streamers Drive Rapid Accretion in 3D GRMHD Simulations of Warped Thin Disks
Authors:
Nicholas Kaaz,
Matthew T. P. Liska,
Jonatan Jacquemin-Ide,
Zachary L. Andalman,
Gibwa Musoke,
Alexander Tchekhovskoy,
Oliver Porth
Abstract:
The angular momentum of gas feeding a black hole (BH) is typically misaligned with respect to the BH spin, resulting in a tilted accretion disk. Rotation of the BH drags the surrounding space-time, manifesting as Lense-Thirring torques that lead to disk precession and warping. We study these processes by simulating a thin ($H/r=0.02$), highly tilted ($\mathcal{T}=65^\circ$) accretion disk around a…
▽ More
The angular momentum of gas feeding a black hole (BH) is typically misaligned with respect to the BH spin, resulting in a tilted accretion disk. Rotation of the BH drags the surrounding space-time, manifesting as Lense-Thirring torques that lead to disk precession and warping. We study these processes by simulating a thin ($H/r=0.02$), highly tilted ($\mathcal{T}=65^\circ$) accretion disk around a rapidly rotating ($a=0.9375$) BH at extremely high resolutions, which we performed using the general-relativistic magnetohydrodynamic (GRMHD) code H-AMR. The disk becomes significantly warped and continuously tears into two individually precessing sub-disks. We find that mass accretion rates far exceed the standard $α$-viscosity expectations. We identify two novel dissipation mechanisms specific to warped disks that are the main drivers of accretion, distinct from the local turbulent stresses that are usually thought to drive accretion. In particular, we identify extreme scale height oscillations that occur twice an orbit throughout our disk. When the scale height compresses, `nozzle' shocks form, dissipating orbital energy and driving accretion. Separate from this phenomenon, there is also extreme dissipation at the location of the tear. This leads to the formation of low-angular momentum `streamers' that rain down onto the inner sub-disk, shocking it. The addition of low angular momentum gas to the inner sub-disk causes it to rapidly accrete, even when it is transiently aligned with the BH spin and thus unwarped. These mechanisms, if general, significantly modify the standard accretion paradigm. Additionally, they may drive structural changes on much shorter timescales than expected in $α$-disks, potentially explaining some of the extreme variability observed in active galactic nuclei.
△ Less
Submitted 6 May, 2024; v1 submitted 18 October, 2022;
originally announced October 2022.
-
Relativistic Signatures of Flux Eruption Events Near Black Holes
Authors:
Zachary Gelles,
Koushik Chatterjee,
Michael Johnson,
Bart Ripperda,
Matthew Liska
Abstract:
Images of supermassive black holes produced using very long baseline interferometry provide a pathway to directly observing effects of a highly curved spacetime, such as a bright ``photon ring'' that arises from strongly lensed emission. In addition, the emission near supermassive black holes is highly variable, with bright high-energy flares regularly observed. We demonstrate that intrinsic varia…
▽ More
Images of supermassive black holes produced using very long baseline interferometry provide a pathway to directly observing effects of a highly curved spacetime, such as a bright ``photon ring'' that arises from strongly lensed emission. In addition, the emission near supermassive black holes is highly variable, with bright high-energy flares regularly observed. We demonstrate that intrinsic variability can introduce prominent associated changes in the relative brightness of the photon ring. We analyze both semianalytic toy models and GRMHD simulations with magnetic flux eruption events, showing that they each exhibit a characteristic ``loop'' in the space of relative photon ring brightness versus total flux density. For black holes viewed at high inclination, the relative photon ring brightness can change by an order of magnitude, even with variations in total flux density that are comparatively mild. We show that gravitational lensing, Doppler boosting, and magnetic field structure all significantly affect this feature, and we discuss the prospects for observing it in observations of M87* and Sgr A* with the next-generation Event Horizon Telescope.
△ Less
Submitted 24 November, 2022; v1 submitted 13 October, 2022;
originally announced October 2022.
-
Unraveling Twisty Linear Polarization Morphologies in Black Hole Images
Authors:
Razieh Emami,
Angelo Ricarte,
George N. Wong,
Daniel Palumbo,
Dominic Chang,
Sheperd S. Doeleman,
Avery Broaderick,
Ramesh Narayan,
Maciek Wielgus,
Lindy Blackburn,
Ben S. Prather,
Andrew A. Chael,
Richard Anantua,
Koushik Chatterjee,
Ivan Marti-Vidal,
Jose L. Gomez,
Kazunori Akiyama,
Matthew Liska,
Lars Hernquist,
Grant Tremblay,
Mark Vogelsberger,
Charles Alcock,
Randall Smith,
James Steiner,
Paul Tiede
, et al. (1 additional authors not shown)
Abstract:
We investigate general relativistic magnetohydrodynamic simulations (GRMHD) to determine the physical origin of the twisty patterns of linear polarization seen in spatially resolved black hole images and explain their morphological dependence on black hole spin. By characterising the observed emission with a simple analytic ring model, we find that the twisty morphology is determined by the magnet…
▽ More
We investigate general relativistic magnetohydrodynamic simulations (GRMHD) to determine the physical origin of the twisty patterns of linear polarization seen in spatially resolved black hole images and explain their morphological dependence on black hole spin. By characterising the observed emission with a simple analytic ring model, we find that the twisty morphology is determined by the magnetic field structure in the emitting region. Moreover, the dependence of this twisty pattern on spin can be attributed to changes in the magnetic field geometry that occur due to the frame dragging. By studying an analytic ring model, we find that the roles of Doppler boosting and lensing are subdominant. Faraday rotation may cause a systematic shift in the linear polarization pattern, but we find that its impact is subdominant for models with strong magnetic fields and modest ion-to-electron temperature ratios. Models with weaker magnetic fields are much more strongly affected by Faraday rotation and have more complicated emission geometries than can be captured by a ring model. However, these models are currently disfavoured by the recent EHT observations of M87*. Our results suggest that linear polarization maps can provide a probe of the underlying magnetic field structure around a black hole, which may then be usable to indirectly infer black hole spins. The generality of these results should be tested with alternative codes, initial conditions, and plasma physics prescriptions.
△ Less
Submitted 28 March, 2023; v1 submitted 3 October, 2022;
originally announced October 2022.
-
Gas Morphology of Milky Way-like Galaxies in the TNG50 Simulation: Signals of Twisting and Stretching
Authors:
Thomas K. Waters,
Colton Peterson,
Razieh Emami,
Xuejian Shen,
Lars Hernquist,
Randall Smith,
Mark Vogelsberger,
Charles Alcock,
Grant Tremblay,
Matthew Liska,
John C. Forbes,
Jorge Moreno
Abstract:
We present an in-depth analysis of gas morphologies for a sample of 25 Milky Way-like galaxies from the IllustrisTNG TNG50 simulation. We constrain the morphology of cold, warm, hot gas, and gas particles as a whole using a Local Shell Iterative Method (LSIM) and explore its observational implications by computing the hard-to-soft X-ray ratio, which ranges between $10^{-3}$-$10^{-2}$ in the inner…
▽ More
We present an in-depth analysis of gas morphologies for a sample of 25 Milky Way-like galaxies from the IllustrisTNG TNG50 simulation. We constrain the morphology of cold, warm, hot gas, and gas particles as a whole using a Local Shell Iterative Method (LSIM) and explore its observational implications by computing the hard-to-soft X-ray ratio, which ranges between $10^{-3}$-$10^{-2}$ in the inner $\sim 50 \rm kpc$ of the distribution and $10^{-5}$-$10^{-4}$ at the outer portion of the hot gas distribution. We group galaxies into three main categories: simple, stretched, and twisted. These categories are based on the radial reorientation of the principal axes of the reduced inertia tensor. We find that a vast majority ($77\%$) of the galaxies in our sample exhibit twisting patterns in their radial profiles. Additionally, we present detailed comparisons between 1) the gaseous distributions belonging to individual temperature regimes, 2) the cold gas distributions and stellar distributions, and 3) the gaseous distributions and dark matter (DM) halos. We find a strong correlation between the morphological properties of the cold gas and stellar distributions. Furthermore, we find a correlation between gaseous distributions with DM halo that increases with gas temperature, implying that we may use the warm-hot gaseous morphology as a tracer to probe the DM morphology. Finally, we show gaseous distributions exhibit significantly more prolate morphologies than the stellar distributions and DM halos, which we hypothesize is due to stellar and AGN feedback.
△ Less
Submitted 5 January, 2024; v1 submitted 3 October, 2022;
originally announced October 2022.
-
On the jet-ejecta interaction in 3D GRMHD simulations of binary neutron star merger aftermath
Authors:
Ore Gottlieb,
Serena Moseley,
Teresita Ramirez-Aguilar,
Ariadna Murguia-Berthier,
Matthew Liska,
Alexander Tchekhovskoy
Abstract:
Short $γ$-ray burst (sGRB) jets form in the aftermath of a neutron star merger, drill through disk winds and dynamical ejecta, and extend over four to five orders of magnitude in distance before breaking out of the ejecta. We present the first 3D general-relativistic magnetohydrodynamic sGRB simulations to span this enormous scale separation. They feature three possible outcomes: jet+cocoon, cocoo…
▽ More
Short $γ$-ray burst (sGRB) jets form in the aftermath of a neutron star merger, drill through disk winds and dynamical ejecta, and extend over four to five orders of magnitude in distance before breaking out of the ejecta. We present the first 3D general-relativistic magnetohydrodynamic sGRB simulations to span this enormous scale separation. They feature three possible outcomes: jet+cocoon, cocoon, and neither. Typical sGRB jets break out of the dynamical ejecta if (i) the bound ejecta's isotropic equivalent mass along the pole at the time of the BH formation is $ \lesssim10^{-4}~{\rm M_{\odot}} $, setting a limit on the delay time between the merger and BH formation, otherwise, the jets perish inside the ejecta and leave the jet-inflated cocoon to power a low-luminosity sGRB; (ii) the post-merger remnant disk contains strong large-scale vertical magnetic field, $\gtrsim10^{15}$ G; and (iii) if the jets are weak ($\lesssim10^{50}$ erg), the ejecta's isotropic equivalent mass along the pole must be small ($\lesssim10^{-2}~{\rm M_{\odot}}$). Generally, the jet structure is shaped by the early interaction with disk winds rather than the dynamical ejecta. As long as our jets break out of the ejecta, they retain a significant magnetization ($\lesssim1$), suggesting that magnetic reconnection is a fundamental property of sGRB emission. The angular structure of the outflow isotropic equivalent energy after breakout consistently features a flat core followed by a steep power-law distribution (slope $\gtrsim3$), similar to hydrodynamic jets. In the cocoon-only outcome, the dynamical ejecta broadens the outflow angular distribution and flattens it (slope $\sim1.5$).
△ Less
Submitted 23 June, 2022; v1 submitted 3 May, 2022;
originally announced May 2022.
-
Black hole to photosphere: 3D GRMHD simulations of collapsars reveal wobbling and hybrid composition jets
Authors:
Ore Gottlieb,
Matthew Liska,
Alexander Tchekhovskoy,
Omer Bromberg,
Aretaios Lalakos,
Dimitrios Giannios,
Philipp Mösta
Abstract:
Long-duration $γ$-ray bursts (GRBs) accompany the collapse of massive stars and carry information about the central engine. However, no 3D models have been able to follow these jets from their birth by a black-hole (BH) to the photosphere. We present the first such 3D general-relativity magnetohydrodynamic simulations, which span over 6 orders of magnitude in space and time. The collapsing stellar…
▽ More
Long-duration $γ$-ray bursts (GRBs) accompany the collapse of massive stars and carry information about the central engine. However, no 3D models have been able to follow these jets from their birth by a black-hole (BH) to the photosphere. We present the first such 3D general-relativity magnetohydrodynamic simulations, which span over 6 orders of magnitude in space and time. The collapsing stellar envelope forms an accretion disk, which drags inwardly the magnetic flux that accumulates around the BH, becomes dynamically important and launches bipolar jets. The jets reach the photosphere at $\sim10^{12}$ cm with an opening angle $θ_j\sim6^\circ$ and a Lorentz factor $Γ_j\lesssim 30$, unbinding $\gtrsim90\%$ of the star. We find that (i) the disk-jet system spontaneously develops misalignment relative to the BH rotational axis. As a result, the jet wobbles with an angle $θ_t\sim12^\circ$, which can naturally explain quiescent times in GRB lightcurves. The effective opening angle for detection $θ_j+θ_t$ suggests that the intrinsic GRB rate is lower by an order of magnitude than standard estimates. This suggests that successful GRBs may be rarer than currently thought and emerge in only $\sim 0.1\%$ of supernovae Ib/c, implying that jets are either not launched or choked inside most supernova Ib/c progenitors. (ii) The magnetic energy in the jet decreases due to mixing with the star, resulting in jets with a hybrid composition of magnetic and thermal components at the photosphere, where $\sim 10\%$ of the gas maintains magnetization $σ\gtrsim 0.1$. This indicates that both a photospheric component and reconnection may play a role in the prompt emission.
△ Less
Submitted 20 June, 2022; v1 submitted 26 April, 2022;
originally announced April 2022.
-
Bridging Bondi and Event Horizon Scales: 3D GRMHD Simulations Reveal X-Shaped Radio Galaxy Morphology
Authors:
Aretaios Lalakos,
Ore Gottlieb,
Nick Kaaz,
Koushik Chatterjee,
Matthew Liska,
Ian M. Christie,
Alexander Tchekhovskoy,
Irina Zhuravleva,
Elena Nokhrina
Abstract:
X-shaped radio galaxies (XRGs) produce misaligned X-shaped jet pairs and make up $\lesssim10\%$ of radio galaxies. XRGs are thought to emerge in galaxies featuring a binary supermassive black hole ($\rm SMBH$), $\rm SMBH$ merger, or large-scale ambient medium asymmetry. We demonstrate that XRG morphology can naturally form without such special, preexisting conditions. Our 3D general-relativistic m…
▽ More
X-shaped radio galaxies (XRGs) produce misaligned X-shaped jet pairs and make up $\lesssim10\%$ of radio galaxies. XRGs are thought to emerge in galaxies featuring a binary supermassive black hole ($\rm SMBH$), $\rm SMBH$ merger, or large-scale ambient medium asymmetry. We demonstrate that XRG morphology can naturally form without such special, preexisting conditions. Our 3D general-relativistic magnetohydrodynamic (GRMHD) simulation for the first time follows magnetized rotating gas from outside the $\rm SMBH$ sphere of influence of radius $R_{\rm B}$ to the $\rm SMBH$ of gravitational radius $R_{\rm g}$, at the largest scale separation $R_{\rm B}/R_{\rm g} = 10^3$ to date. Initially, our axisymmetric system of constant-density hot gas contains weak vertical magnetic field and rotates in an equatorial plane of a rapidly spinning $\rm SMBH$. We seed the gas with small-scale $2\%$-level pressure perturbations. Infalling gas forms an accretion disk, and the $\rm SMBH$ launches relativistically-magnetized collimated jets reaching well outside $R_{\rm B}$. Under the pressure of the infalling gas, the jets intermittently turn on and off, erratically wobble, and inflate pairs of cavities in different directions, resembling an X-shaped jet morphology. Synthetic X-ray images reveal multiple pairs of jet-powered shocks and cavities. Large-scale magnetic flux accumulates on the $\rm SMBH$, becomes dynamically important, and leads to a magnetically arrested disk state. The $\rm SMBH$ accretes at $2\%$ of the Bondi rate ($\dot{M}\simeq2.4\times10^{-3}M_{\odot}\,{\rm yr}^{-1}$ for M87*), and launches twin jets at $η=150\%$ efficiency. These jets are powerful enough ($P_{\rm jets}\simeq2\times10^{44}\,{\rm erg\,s}^{-1}$) to escape along the spin axis and end the short-lived jets state whose transient nature can account for the rarity of XRGs
△ Less
Submitted 16 February, 2022;
originally announced February 2022.
-
Jet Formation in 3D GRMHD Simulations of Bondi-Hoyle-Lyttleton Accretion
Authors:
Nicholas Kaaz,
Ariadna Murguia-Berthier,
Koushik Chatterjee,
Matthew Liska,
Alexander Tchekhovskoy
Abstract:
A black hole (BH) travelling through a uniform, gaseous medium is described by Bondi-Hoyle-Lyttleton (BHL) accretion. If the medium is magnetized, then the black hole can produce relativistic outflows. We performed the first 3D, general-relativistic magnetohydrodynamics simulations of BHL accretion onto rapidly rotating black holes using the code H-AMR, where we mainly varied the strength of a bac…
▽ More
A black hole (BH) travelling through a uniform, gaseous medium is described by Bondi-Hoyle-Lyttleton (BHL) accretion. If the medium is magnetized, then the black hole can produce relativistic outflows. We performed the first 3D, general-relativistic magnetohydrodynamics simulations of BHL accretion onto rapidly rotating black holes using the code H-AMR, where we mainly varied the strength of a background magnetic field that threads the medium. We found that the ensuing accretion continuously drags to the BH the magnetic flux, which accumulates near the event horizon until it becomes dynamically important. Depending on the strength of the background magnetic field, the BHs can sometimes launch relativistic jets with high enough power to drill out of the inner accretion flow, become bent by the headwind, and escape to large distances. While for stronger background magnetic fields the jets are continuously powered, at weaker field strengths they are intermittent, turning on and off depending on the fluctuating gas and magnetic flux distributions near the event horizon. We find that our jets reach extremely high efficiencies of $\sim100-300\%$, even in the absence of an accretion disk. We also calculated the drag forces exerted by the gas onto to the BH, finding that the presence of magnetic fields causes drag forces to be much less efficient than in unmagnetized BHL accretion, and sometimes become negative, accelerating the BH rather than slowing it down. Our results extend classical BHL accretion to rotating BHs moving through magnetized media and demonstrate that accretion and drag are significantly altered in this environment.
△ Less
Submitted 10 May, 2023; v1 submitted 27 January, 2022;
originally announced January 2022.
-
Formation of Magnetically Truncated Accretion Disks in 3D Radiation-Transport Two-Temperature GRMHD Simulations
Authors:
M. T. P. Liska,
G. Musoke,
A. Tchekhovskoy,
O. Porth,
A. M. Beloborodov
Abstract:
Multi-wavelength observations suggest that the accretion disk in the hard and intermediate states of X-ray binaries (XRBs) and active galactic nuclei (AGN) transitions from a cold, thin disk at large distances into a hot, thick flow close to the black hole. However, the formation, structure and dynamics of such truncated disks are poorly constrained due to the complexity of the thermodynamic, magn…
▽ More
Multi-wavelength observations suggest that the accretion disk in the hard and intermediate states of X-ray binaries (XRBs) and active galactic nuclei (AGN) transitions from a cold, thin disk at large distances into a hot, thick flow close to the black hole. However, the formation, structure and dynamics of such truncated disks are poorly constrained due to the complexity of the thermodynamic, magnetic, and radiative processes involved. We present the first radiation-transport two-temperature general relativistic magnetohydrodynamic (GRMHD) simulations of truncated disks radiating at ~35% of the Eddington luminosity with and without large-scale poloidal magnetic flux. We demonstrate that when a geometrically-thin accretion disk is threaded by large-scale net poloidal magnetic flux, it self-consistently transitions at small radii into a two-phase medium of cold gas clumps floating through a hot, magnetically dominated corona. This transition occurs at a well-defined truncation radius determined by the distance out to which the disk is saturated with magnetic flux. The average ion and electron temperatures in the semi-opaque corona reach, respectively, T_i ~ 10^10K and T_e ~ 5 10^8K. The system produces radiation, powerful collimated jets and broader winds at the total energy efficiency exceeding ~90%, the highest ever energy extraction efficiency from a spinning black hole by a radiatively efficient flow in a GRMHD simulation. This is consistent with jetted ejections observed during XRB outbursts. The two-phase medium may naturally lead to broadened iron line emission observed in the hard state.
△ Less
Submitted 9 August, 2022; v1 submitted 10 January, 2022;
originally announced January 2022.
-
Disk Tearing Leads to Low and High Frequency Quasi Periodic Oscillations in a GRMHD Simulation of a Thin Accretion Disk
Authors:
G. Musoke,
M. Liska,
O. Porth,
M. van der Klis,
A. Ingram
Abstract:
Black hole X-ray binaries (BHXRBs) display a wide range of variability phenomena, from long duration spectral state changes to short-term broadband variability and quasi-periodic oscillations (QPOs). A particularly puzzling aspect is the production of QPOs, which -- if properly understood -- could be used as a powerful diagnostic tool of black hole accretion and evolution. In this work we analyse…
▽ More
Black hole X-ray binaries (BHXRBs) display a wide range of variability phenomena, from long duration spectral state changes to short-term broadband variability and quasi-periodic oscillations (QPOs). A particularly puzzling aspect is the production of QPOs, which -- if properly understood -- could be used as a powerful diagnostic tool of black hole accretion and evolution. In this work we analyse a high resolution three-dimensional general relativistic magnetohydrodynamic simulation of a geometrically thin accretion disk which is tilted by $65^{\circ}$ with respect to the black hole spin axis. We find that the Lense-Thirring torque from the rapidly spinning 10 $M_\odot$ black hole causes several sub-disks to tear off within $\sim 10-20$ gravitational radii. Tearing occurs in cycles on timescales of seconds. During each tearing cycle the inner sub-disk precesses for 1-5 periods before it falls into the black hole. We find a precession frequency of $\sim 3\rm Hz$, consistent with observed low-frequency QPOs. In addition, we find a high frequency QPO (HFQPO) with centroid frequency of $\sim55$Hz in the power spectra of the mass-weighted radius of the inner disk. This signal is caused by radial epicyclic oscillations of a dense ring of gas at the tearing radius, which strongly suggests a corresponding modulation of the X-ray lightcurve and may thus explain some of the observed HFQPOs.
△ Less
Submitted 9 January, 2022;
originally announced January 2022.
-
Black hole flares: ejection of accreted magnetic flux through 3D plasmoid-mediated reconnection
Authors:
Bart Ripperda,
Matthew Liska,
Koushik Chatterjee,
Gibwa Musoke,
Alexander A. Philippov,
Sera B. Markoff,
Alexander Tchekhovskoy,
Ziri Younsi
Abstract:
Magnetic reconnection can power bright, rapid flares originating from the inner magnetosphere of accreting black holes. We conduct extremely high resolution ($5376\times2304\times2304$ cells) general-relativistic magnetohydrodynamics simulations, capturing plasmoid-mediated reconnection in a 3D magnetically arrested disk for the first time. We show that an equatorial, plasmoid-unstable current she…
▽ More
Magnetic reconnection can power bright, rapid flares originating from the inner magnetosphere of accreting black holes. We conduct extremely high resolution ($5376\times2304\times2304$ cells) general-relativistic magnetohydrodynamics simulations, capturing plasmoid-mediated reconnection in a 3D magnetically arrested disk for the first time. We show that an equatorial, plasmoid-unstable current sheet forms in a transient, non-axisymmetric, low-density magnetosphere within the inner few Schwarzschild radii. Magnetic flux bundles escape from the event horizon through reconnection at the universal plasmoid-mediated rate in this current sheet. The reconnection feeds on the highly-magnetized plasma in the jets and heats the plasma that ends up trapped in flux bundles to temperatures proportional to the jet's magnetization. The escaped flux bundles can complete a full orbit as low-density hot spots, consistent with Sgr A$^{*}$ observations by the GRAVITY interferometer. Reconnection near the horizon produces sufficiently energetic plasma to explain flares from accreting black holes, such as the TeV emission observed from M87. The drop in mass accretion rate during the flare, and the resulting low-density magnetosphere make it easier for very high energy photons produced by reconnection-accelerated particles to escape. The extreme resolution results in a converged plasmoid-mediated reconnection rate that directly determines the timescales and properties of the flare.
△ Less
Submitted 6 February, 2022; v1 submitted 30 September, 2021;
originally announced September 2021.
-
Black hole to breakout: 3D GRMHD simulations of collapsar jets reveal a wide range of transients
Authors:
Ore Gottlieb,
Aretaios Lalakos,
Omer Bromberg,
Matthew Liska,
Alexander Tchekhovskoy
Abstract:
We present a suite of the first 3D GRMHD collapsar simulations, which extend from the self-consistent jet launching by an accreting Kerr black hole (BH) to the breakout from the star. We identify three types of outflows, depending on the angular momentum, $ l $, of the collapsing material and the magnetic field, $ B $, on the BH horizon: (i) Subrelativistic outflow (low $ l $ and high $ B $), (ii)…
▽ More
We present a suite of the first 3D GRMHD collapsar simulations, which extend from the self-consistent jet launching by an accreting Kerr black hole (BH) to the breakout from the star. We identify three types of outflows, depending on the angular momentum, $ l $, of the collapsing material and the magnetic field, $ B $, on the BH horizon: (i) Subrelativistic outflow (low $ l $ and high $ B $), (ii) Stationary accretion shock instability (SASI; high $ l $ and low $ B $), (iii) Relativistic jets (high $ l $ and high $ B $). In the absence of jets, free-fall of the stellar envelope provides a good estimate for the BH accretion rate. Jets can substantially suppress the accretion rate, and their duration can be limited by the magnetization profile in the star. We find that progenitors with large (steep) inner density power-law indices ($ \gtrsim 2 $), face extreme challenges as gamma-ray burst (GRB) progenitors due to excessive luminosity, global time evolution in the lightcurve throughout the burst and short breakout times, inconsistent with observations. Our results suggest that the wide variety of observed explosion appearances (supernova/supernova+GRB/low-luminosity GRBs) and the characteristics of the emitting relativistic outflows (luminosity and duration) can be naturally explained by the differences in the progenitor structure. Our simulations reveal several important jet features: (i) strong magnetic dissipation inside the star, resulting in weakly magnetized jets by breakout that may have significant photospheric emission and (ii) spontaneous emergence of tilted accretion disk-jet flows, even in the absence of any tilt in the progenitor.
△ Less
Submitted 26 December, 2021; v1 submitted 29 September, 2021;
originally announced September 2021.
-
Spectral and Imaging properties of Sgr A* from High-Resolution 3D GRMHD Simulations with Radiative Cooling
Authors:
Doosoo Yoon,
Koushik Chatterjee,
Sera Markoff,
David van Eijnatten,
Ziri Younsi,
Matthew Liska,
Alexander Tchekhovskoy
Abstract:
The candidate supermassive black hole in the Galactic Centre, Sagittarius A* (Sgr A*), is known to be fed by a radiatively inefficient accretion flow (RIAF), inferred by its low accretion rate. Consequently, radiative cooling has in general been overlooked in the study of Sgr A*. However, the radiative properties of the plasma in RIAFs are poorly understood. In this work, using full 3D general-rel…
▽ More
The candidate supermassive black hole in the Galactic Centre, Sagittarius A* (Sgr A*), is known to be fed by a radiatively inefficient accretion flow (RIAF), inferred by its low accretion rate. Consequently, radiative cooling has in general been overlooked in the study of Sgr A*. However, the radiative properties of the plasma in RIAFs are poorly understood. In this work, using full 3D general-relativistic magneto-hydrodynamical simulations, we study the impact of radiative cooling on the dynamical evolution of the accreting plasma, presenting spectral energy distributions and synthetic sub-millimeter images generated from the accretion flow around Sgr A*. These simulations solve the approximated equations for radiative cooling processes self-consistently, including synchrotron, bremsstrahlung, and inverse Compton processes. We find that radiative cooling plays an increasingly important role in the dynamics of the accretion flow as the accretion rate increases: the mid-plane density grows and the infalling gas is less turbulent as cooling becomes stronger. The changes in the dynamical evolution become important when the accretion rate is larger than $10^{-8}\,M_{\odot}~{\rm yr}^{-1}$ ($\gtrsim 10^{-7} \dot{M}_{\rm Edd}$, where $\dot{M}_{\rm Edd}$ is the Eddington accretion rate). The resulting spectra in the cooled models also differ from those in the non-cooled models: the overall flux, including the peak values at the sub-mm and the far-UV, is slightly lower as a consequence of a decrease in the electron temperature. Our results suggest that radiative cooling should be carefully taken into account in modelling Sgr A* and other low-luminosity active galactic nuclei that have a mass accretion rate of $\dot{M} > 10^{-7}\,\dot{M}_{\rm Edd}$.
△ Less
Submitted 29 September, 2020;
originally announced September 2020.
-
Tidal Disruption Disks Formed and Fed by Stream-Stream and Stream-Disk Interactions in Global GRHD Simulations
Authors:
Zachary L. Andalman,
Matthew T. P. Liska,
Alexander Tchekhovskoy,
Eric R. Coughlin,
Nicholas Stone
Abstract:
When a star passes close to a supermassive black hole (BH), the BH's tidal forces rip it apart into a thin stream, leading to a tidal disruption event (TDE). In this work, we study the post-disruption phase of TDEs in general relativistic hydrodynamics (GRHD) using our GPU-accelerated code H-AMR. We carry out the first grid-based simulation of a deep-penetration TDE ($β$=7) with realistic system p…
▽ More
When a star passes close to a supermassive black hole (BH), the BH's tidal forces rip it apart into a thin stream, leading to a tidal disruption event (TDE). In this work, we study the post-disruption phase of TDEs in general relativistic hydrodynamics (GRHD) using our GPU-accelerated code H-AMR. We carry out the first grid-based simulation of a deep-penetration TDE ($β$=7) with realistic system parameters: a black-hole-to-star mass ratio of $10^6$, a parabolic stellar trajectory, and a nonzero BH spin. We also carry out a simulation of a tilted TDE whose stellar orbit is inclined relative to the BH midplane. We show that for our aligned TDE, an accretion disk forms due to the dissipation of orbital energy with $\sim$20 percent of the infalling material reaching the BH. The dissipation is initially dominated by violent self-intersections and later by stream-disk interactions near the pericenter. The self-intersections completely disrupt the incoming stream, resulting in five distinct self-intersection events separated by approximately 12 hours and a flaring in the accretion rate. We also find that the disk is eccentric with mean eccentricity e$\approx$0.88. For our tilted TDE, we find only partial self-intersections due to nodal precession near pericenter. Although these partial intersections eject gas out of the orbital plane, an accretion disk still forms with a similar accreted fraction of the material to the aligned case. These results have important implications for disk formation in realistic tidal disruptions. For instance, the periodicity in accretion rate induced by the complete stream disruption may explain the flaring events from Swift J1644+57.
△ Less
Submitted 23 November, 2021; v1 submitted 11 August, 2020;
originally announced August 2020.
-
Observational signatures of disk and jet misalignment in images of accreting black holes
Authors:
Koushik Chatterjee,
Ziri Younsi,
Matthew Liska,
Alexander Tchekhovskoy,
Sera Markoff,
Doosoo Yoon,
David van Eijnatten,
Casper Hesp,
Adam Ingram,
Michiel van der Klis
Abstract:
Black hole accretion is one of nature's most efficient energy extraction processes. When gas falls in, a significant fraction of its gravitational binding energy is either converted into radiation or flows outwards in the form of black hole-driven jets and disk-driven winds. Recently, the Event Horizon Telescope (EHT), an Earth-size sub-millimetre radio interferometer, captured the first images of…
▽ More
Black hole accretion is one of nature's most efficient energy extraction processes. When gas falls in, a significant fraction of its gravitational binding energy is either converted into radiation or flows outwards in the form of black hole-driven jets and disk-driven winds. Recently, the Event Horizon Telescope (EHT), an Earth-size sub-millimetre radio interferometer, captured the first images of M87's black hole. These images were analysed and interpreted using general-relativistic magnetohydrodynamics (GRMHD) models of accretion disks with rotation axes aligned with the black hole spin axis. However, since infalling gas is often insensitive to the black hole spin direction, misalignment between accretion disk and black hole spin may be a common occurrence in nature. In this work, we use the general-relativistic radiative transfer (GRRT) code \texttt{BHOSS} to calculate the first synthetic radio images of (highly) tilted disk/jet models generated by our GPU-accelerated GRMHD code \texttt{HAMR}. While the tilt does not have a noticeable effect on the system dynamics beyond a few tens of gravitational radii from the black hole, the warping of the disk and jet can imprint observable signatures in EHT images on smaller scales. Comparing the images from our GRMHD models to the 43 GHz and 230 GHz EHT images of M87, we find that M87 may feature a tilted disk/jet system. Further, tilted disks and jets display significant time variability in the 230 GHz flux that can be further tested by longer-duration EHT observations of M87.
△ Less
Submitted 4 September, 2020; v1 submitted 19 February, 2020;
originally announced February 2020.
-
H-AMR: A New GPU-accelerated GRMHD Code for Exascale Computing With 3D Adaptive Mesh Refinement and Local Adaptive Time-stepping
Authors:
M. T. P. Liska,
K. Chatterjee,
D. Issa,
D. Yoon,
N. Kaaz,
A. Tchekhovskoy,
D. van Eijnatten,
G. Musoke,
C. Hesp,
V. Rohoza,
S. Markoff,
A. Ingram,
M. van der Klis
Abstract:
General-relativistic magnetohydrodynamic (GRMHD) simulations have revolutionized our understanding of black hole accretion. Here, we present a graphics processing unit (GPU) accelerated GRMHD code \hammer{} with multi-faceted optimizations that, collectively, accelerate computation by 2-5 orders of magnitude for a wide range of applications. Firstly, it introduces a spherical grid with 3D adaptive…
▽ More
General-relativistic magnetohydrodynamic (GRMHD) simulations have revolutionized our understanding of black hole accretion. Here, we present a graphics processing unit (GPU) accelerated GRMHD code \hammer{} with multi-faceted optimizations that, collectively, accelerate computation by 2-5 orders of magnitude for a wide range of applications. Firstly, it introduces a spherical grid with 3D adaptive mesh refinement that operates in each of the 3 dimensions independently. This allows us to circumvent the Courant condition near the polar singularity, which otherwise cripples high-resolution computational performance. Secondly, we demonstrate that local adaptive time-stepping (LAT) on a logarithmic spherical-polar grid accelerates computation by a factor of $\lesssim10$ compared to traditional hierarchical time-stepping approaches. Jointly, these unique features lead to an effective speed of $\sim10^9$ zone-cycles-per-second-per-node on 5,400 NVIDIA V100 GPUs (i.e., 900 nodes of the OLCF Summit supercomputer). We illustrate \hammer{}'s computational performance by presenting the first GRMHD simulation of a tilted thin accretion disk threaded by a toroidal magnetic field around a rapidly spinning black hole. With an effective resolution of $13$,$440\times4$,$608\times8$,$092$ cells, and a total of $\lesssim22$ billion cells and $\sim0.65\times10^8$ timesteps, it is among the largest astrophysical simulations ever performed. We find that frame-dragging by the black hole tears up the disk into two independently precessing sub-disks. The innermost sub-disk rotation axis intermittently aligns with the black hole spin, demonstrating for the first time that such long-sought alignment is possible in the absence of large-scale poloidal magnetic fields.
△ Less
Submitted 26 October, 2022; v1 submitted 20 December, 2019;
originally announced December 2019.
-
Disc Tearing and Bardeen-Petterson Alignment in GRMHD Simulations of Highly Tilted Thin Accretion Discs
Authors:
M. Liska,
C. Hesp,
A. Tchekhovskoy,
A. Ingram,
M. van der Klis,
S. B. Markoff,
M. Van Moer
Abstract:
Luminous active galactic nuclei (AGN) and X-Ray binaries (XRBs) often contain geometrically thin, radiatively cooled accretion discs. According to theory, these are -- in many cases -- initially highly misaligned with the black hole equator. In this work, we present the first general relativistic magnetohydrodynamic simulations of very thin (h/r~0.015-0.05) accretion discs around rapidly spinning…
▽ More
Luminous active galactic nuclei (AGN) and X-Ray binaries (XRBs) often contain geometrically thin, radiatively cooled accretion discs. According to theory, these are -- in many cases -- initially highly misaligned with the black hole equator. In this work, we present the first general relativistic magnetohydrodynamic simulations of very thin (h/r~0.015-0.05) accretion discs around rapidly spinning (a~0.9) black holes and tilted by 45-65 degrees. We show that the inner regions of the discs with h/r<0.03 align with the black hole equator, though out to smaller radii than predicted by analytic work. The inner aligned and outer misaligned disc regions are separated by a sharp break in tilt angle accompanied by a sharp drop in density. We find that frame-dragging by the spinning black hole overpowers the disc viscosity, which is self-consistently produced by magnetized turbulence, tearing the disc apart and forming a rapidly precessing inner sub-disc surrounded by a slowly precessing outer sub-disc. We find that the system produces a pair of relativistic jets for all initial tilt values. At small distances the black hole launched jets precess rapidly together with the inner sub-disc, whereas at large distances they partially align with the outer sub-disc and precess more slowly. If the tearing radius can be modeled accurately in future work, emission model independent measurements of black hole spin based on precession-driven quasi-periodic oscillations may become possible.
△ Less
Submitted 10 August, 2022; v1 submitted 17 April, 2019;
originally announced April 2019.
-
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.
-
Accelerating AGN jets to parsec scales using general relativistic MHD simulations
Authors:
Koushik Chatterjee,
Matthew Liska,
Alexander Tchekhovskoy,
Sera B. Markoff
Abstract:
Accreting black holes produce collimated outflows, or jets, that traverse many orders of magnitude in distance, accelerate to relativistic velocities, and collimate into tight opening angles. Of these, perhaps the least understood is jet collimation due to the interaction with the ambient medium. In order to investigate this interaction, we carried out axisymmetric general relativistic magnetohydr…
▽ More
Accreting black holes produce collimated outflows, or jets, that traverse many orders of magnitude in distance, accelerate to relativistic velocities, and collimate into tight opening angles. Of these, perhaps the least understood is jet collimation due to the interaction with the ambient medium. In order to investigate this interaction, we carried out axisymmetric general relativistic magnetohydrodynamic simulations of jets produced by a large accretion disc, spanning over 5 orders of magnitude in time and distance, at an unprecedented resolution. Supported by such a disc, the jet attains a parabolic shape, similar to the M87 galaxy jet, and the product of the Lorentz factor and the jet half-opening angle, $γθ\ll 1$, similar to values found from very long baseline interferometry (VLBI) observations of active galactic nuclei (AGN) jets; this suggests extended discs in AGN. We find that the interaction between the jet and the ambient medium leads to the development of pinch instabilities, which produce significant radial and lateral variability across the jet by converting magnetic and kinetic energy into heat. Thus pinched regions in the jet can be detectable as radiating hotspots and may provide an ideal site for particle acceleration. Pinching also causes gas from the ambient medium to become squeezed between magnetic field lines in the jet, leading to enhanced mass-loading of the jet and potentially contributing to the spine-sheath structure observed in AGN outflows.
△ Less
Submitted 9 May, 2019; v1 submitted 5 April, 2019;
originally announced April 2019.
-
A Phase Lag between Disk and Corona in GRMHD Simulations of Precessing Tilted Accretion Disks
Authors:
M. Liska,
C. Hesp,
A. Tchekhovskoy,
A. Ingram,
M. van der Klis,
S. B. Markoff
Abstract:
In the course of its evolution, a black hole (BH) accretes gas from a wide range of directions. Given a random accretion event, the typical angular momentum of an accretion disc would be tilted by $\sim$60$^\circ$ relative to the BH spin. Misalignment causes the disc to precess at a rate that increases with BH spin and tilt angle. We present the first general-relativistic magnetohydrodynamic (GRMH…
▽ More
In the course of its evolution, a black hole (BH) accretes gas from a wide range of directions. Given a random accretion event, the typical angular momentum of an accretion disc would be tilted by $\sim$60$^\circ$ relative to the BH spin. Misalignment causes the disc to precess at a rate that increases with BH spin and tilt angle. We present the first general-relativistic magnetohydrodynamic (GRMHD) simulations spanning a full precession period of highly tilted (60$^\circ$), moderately thin ($h/r=0.1$) accretion discs around a rapidly spinning ($a\simeq0.9$) BH. While the disc and jets precess in phase, we find that the corona, sandwiched between the two, lags behind by $\gtrsim 10^{\circ}$. For spectral models of BH accretion, the implication is that hard non-thermal (corona) emission lags behind the softer (disc) emission, thus potentially explaining some properties of the hard energy lags seen in Type-C low frequency quasi-periodic oscillations in X-Ray binaries. While strong jets are unaffected by this disc-corona lag, weak jets stall when encountering the lagging corona at distances $r \sim 100$ black hole radii. This interaction may quench large-scale jet formation.
△ Less
Submitted 29 January, 2023; v1 submitted 17 January, 2019;
originally announced January 2019.
-
Bardeen-Petterson Alignment, Jets and Magnetic Truncation in GRMHD Simulations of Tilted Thin Accretion Discs
Authors:
M. Liska,
A. Tchekhovskoy,
A. Ingram,
M. van der Klis
Abstract:
Prevalent around luminous accreting black holes, thin discs are challenging to resolve in numerical simulations. When the disc and black hole angular momentum vectors are misaligned, the challenge becomes extreme, requiring adaptive meshes to follow the disc proper as it moves through the computational grid. With our new high-performance general relativistic magnetohydrodynamic (GRMHD) code H-AMR…
▽ More
Prevalent around luminous accreting black holes, thin discs are challenging to resolve in numerical simulations. When the disc and black hole angular momentum vectors are misaligned, the challenge becomes extreme, requiring adaptive meshes to follow the disc proper as it moves through the computational grid. With our new high-performance general relativistic magnetohydrodynamic (GRMHD) code H-AMR we have simulated the thinnest accretion disc to date, of aspect ratio H/R~0.03, around a rapidly spinning (a=0.9375) black hole, using a cooling function. Initially tilted at 10 degrees, the disc warps inside r~5 r_g into alignment with the black hole, where r_g is the gravitational radius. This is the first demonstration of Bardeen-Petterson alignment in MHD with viscosity self-consistently generated by magnetized turbulence. The disc develops a low-density high-viscosity (alpha_eff ~ 1.0) magnetic-pressure--dominated inner region at r<25 r_g that rapidly empties itself into the black hole. This inner region may in reality, due to thermal decoupling of ions and electrons, evaporate into a radiatively inefficient accretion flow if, as we propose, the cooling time exceeds the accretion time set by the order unity effective viscosity. We furthermore find the unexpected result that even our very thin disc can sustain large-scale vertical magnetic flux on the black hole, which launches powerful relativistic jets that carry 20-50% of the accretion power along the angular momentum vector of the outer tilted disc, providing a potential explanation for the origin of jets in radio-loud quasars.
△ Less
Submitted 10 August, 2022; v1 submitted 1 October, 2018;
originally announced October 2018.
-
Large-Scale Poloidal Magnetic Field Dynamo Leads to Powerful Jets in GRMHD Simulations of Black Hole Accretion with Toroidal Field
Authors:
M. T. P. Liska,
A. Tchekhovskoy,
E. Quataert
Abstract:
Accreting black holes (BHs) launch relativistic collimated jets, across many decades in luminosity and mass, suggesting the jet launching mechanism is universal, robust and scale-free. Theoretical models and general relativistic magnetohydrodynamic (GRMHD) simulations indicate that the key jet-making ingredient is large-scale poloidal magnetic flux. However, its origin is uncertain, and it is unkn…
▽ More
Accreting black holes (BHs) launch relativistic collimated jets, across many decades in luminosity and mass, suggesting the jet launching mechanism is universal, robust and scale-free. Theoretical models and general relativistic magnetohydrodynamic (GRMHD) simulations indicate that the key jet-making ingredient is large-scale poloidal magnetic flux. However, its origin is uncertain, and it is unknown if it can be generated in situ or dragged inward from the ambient medium. Here, we use the GPU-accelerated GRMHD code H-AMR to study global 3D BH accretion at unusually high resolutions more typical of local shearing box simulations. We demonstrate that turbulence in a radially-extended accretion disc can generate large-scale poloidal magnetic flux in situ, even when starting from a purely toroidal magnetic field. The flux accumulates around the BH till it becomes dynamically-important, leads to a magnetically arrested disc (MAD), and launches relativistic jets that are more powerful than the accretion flow. The jet power exceeds that of previous GRMHD toroidal field simulations by a factor of 10,000. The jets do not show significant kink or pinch instabilities, accelerate to $γ\sim 10$ over 3 decades in distance, and follow a collimation profile similar to the observed M87 jet.
△ Less
Submitted 10 August, 2022; v1 submitted 12 September, 2018;
originally announced September 2018.
-
Formation of Precessing Jets by Tilted Black-hole Discs in 3D General Relativistic MHD Simulations
Authors:
Matthew Liska,
Casper Hesp,
Alexander Tchekhovskoy,
Adam Ingram,
Michiel van der Klis,
Sera Markoff
Abstract:
Gas falling into a black hole (BH) from large distances is unaware of BH spin direction, and misalignment between the accretion disc and BH spin is expected to be common. However, the physics of tilted discs (e.g., angular momentum transport and jet formation) is poorly understood. Using our new GPU-accelerated code H-AMR, we performed 3D general relativistic magnetohydrodynamic simulations of til…
▽ More
Gas falling into a black hole (BH) from large distances is unaware of BH spin direction, and misalignment between the accretion disc and BH spin is expected to be common. However, the physics of tilted discs (e.g., angular momentum transport and jet formation) is poorly understood. Using our new GPU-accelerated code H-AMR, we performed 3D general relativistic magnetohydrodynamic simulations of tilted thick accretion discs around rapidly spinning BHs, at the highest resolution to date. We explored the limit where disc thermal pressure dominates magnetic pressure, and showed for the first time that, for different magnetic field strengths on the BH, these flows launch magnetized relativistic jets propagating along the rotation axis of the tilted disc (rather than of the BH). If strong large-scale magnetic flux reaches the BH, it bends the inner few gravitational radii of the disc and jets into partial alignment with the BH spin. On longer time scales, the simulated disc-jet system as a whole undergoes Lense-Thirring precession and approaches alignment, demonstrating for the first time that jets can be used as probes of disc precession. When the disc turbulence is well-resolved, our isolated discs spread out, causing both the alignment and precession to slow down.
△ Less
Submitted 5 December, 2022; v1 submitted 20 July, 2017;
originally announced July 2017.
-
Do the Herschel cold clouds in the Galactic halo embody its dark matter?
Authors:
Theo M. Nieuwenhuizen,
Erik F. G. van Heusden,
Matthew T. P. Liska
Abstract:
Recent Herschel/SPIRE maps of the Small and Large Magellanic Clouds (SMC, LMC) exhibit in each thousands of clouds. Observed at 250 microns, they must be cold, T ~ 15 K, hence the name "Herschel cold clouds" (HCCs). From the observed rotational velocity profile and the assumption of spherical symmetry, the Galactic mass density is modeled in a form close to that of an isothermal sphere. If the HCC…
▽ More
Recent Herschel/SPIRE maps of the Small and Large Magellanic Clouds (SMC, LMC) exhibit in each thousands of clouds. Observed at 250 microns, they must be cold, T ~ 15 K, hence the name "Herschel cold clouds" (HCCs). From the observed rotational velocity profile and the assumption of spherical symmetry, the Galactic mass density is modeled in a form close to that of an isothermal sphere. If the HCCs constitute a certain fraction of it, their angular size distribution has a specified shape. A fit to the data deduced from the SMC/LMC maps supports this and yields for their radius 2.5 pc, with a small change when allowing for a spread in HCC radii. There are so many HCCs that they will make up all the missing Halo mass density if there is spherical symmetry and their average mass is of order 15,000 Mo. This compares well with the Jeans mass of circa 40,000 Mo and puts forward that the HCCs are in fact Jeans clusters, constituting all the Galactic dark matter and much of its missing baryons, a conclusion deduced before from a different field of the sky (Nieuwenhuizen, Schild and Gibson 2011). A preliminary analysis of the intensities yields that the Jeans clusters themselves may consist of some billion MACHOs of a few dozen Earth masses. With a size of dozens of solar radii, they would mostly obscure stars in the LMC, SMC and towards the Galactic center, and may thus have been overlooked in microlensing.
△ Less
Submitted 26 November, 2012; v1 submitted 1 October, 2012;
originally announced October 2012.