General Relativistic Magnetohydrodynamic Simulations of Magnetically Choked Accretion Flows around Black Holes

Black hole (BH) accretion flows and jets are qualitatively affected by the presence of ordered magnetic fields. We study fully three-dimensional global general relativistic magnetohydrodynamic (MHD) simulations of radially extended and thick (height $H$ to cylindrical radius $R$ ratio of $|H/R|sim 0.2--1$) accretion flows around BHs with various dimensionless spins ($a/M$, with BH mass $M$) and with initially toroidally-dominated ($phi$-directed) and poloidally-dominated ($R-z$ directed) magnetic fields. Firstly, for toroidal field models and BHs with high enough $|a/M|$, coherent large-scale (i.e. $gg H$) dipolar poloidal magnetic flux patches emerge, thread the BH, and generate transient relativistic jets. Secondly, for poloidal field models, poloidal magnetic flux readily accretes through the disk from large radii and builds-up to a natural saturation point near the BH. For sufficiently high $|a/M|$ or low $|H/R|$ the polar magnetic field compresses the inflow into a geometrically thin highly non-axisymmetric magnetically choked accretion flow (MCAF) within which the standard linear magneto-rotational instability is suppressed. The condition of a highly-magnetized state over most of the horizon is optimal for the Blandford-Znajek mechanism that generates persistent relativistic jets with $gtrsim 100$% efficiency for $|a/M|gtrsim 0.9$. A magnetic Rayleigh-Taylor and Kelvin-Helmholtz unstable magnetospheric interface forms between the compressed inflow and bulging jet magnetosphere, which drives a new jet-disk quasi-periodic oscillation (JD-QPO) mechanism. The high-frequency QPO has spherical harmonic $|m|=1$ mode period of $tausim 70GM/c^3$ for $a/Msim 0.9$ with coherence quality factors $Qgtrsim 10$. [abridged]

A Reconnection Switch to Trigger Gamma-Ray Burst Jet Dissipation

Prompt gamma-ray burst (GRB) emission requires some mechanism to dissipate an ultrarelativistic jet. Internal shocks or some form of electromagnetic dissipation are candidate mechanisms. Any mechanism needs to answer basic questions, such as what is the origin of variability, what radius does dissipation occur at, and how does efficient prompt emission occur. These mechanisms also need to be consistent with how ultrarelativistic jets form and stay baryon pure despite turbulence and electromagnetic reconnection near the compact object and despite stellar entrainment within the collapsar model. We use the latest magnetohydrodynamical models of ultrarelativistic jets to explore some of these questions in the context of electromagnetic dissipation due to the slow collisional and fast collisionless reconnection mechanisms, as often associated with Sweet-Parker and Petschek reconnection, respectively. For a highly magnetized ultrarelativistic jet and typical collapsar parameters, we find that significant electromagnetic dissipation may be avoided until it proceeds catastrophically near the jet photosphere at large radii ($rsim 10^{13}--10^{14}{rm cm}$), by which the jet obtains a high Lorentz factor ($gammasim 100--1000), has a luminosity of $L_j sim 10^{50}--10^{51}ergs$, has observer variability timescales of order 1s (ranging from 0.001-10s), achieves $gammatheta_jsim 10--20 (for opening half-angle $theta_j$) and so is able to produce jet breaks, and has comparable energy available for both prompt and afterglow emission. This reconnection switch mechanism allows for highly efficient conversion of electromagnetic energy into prompt emission and associates the observed prompt GRB pulse temporal structure with dissipation timescales of some number of reconnecting current sheets embedded in the jet.[abridged]

Stability of Relativistic Jets from Rotating, Accreting Black Holes via Fully Three-Dimensional Magnetohydrodynamic Simulations

Rotating magnetized compact objects and their accretion discs can generate strong toroidal magnetic fields driving highly magnetized plasmas into relativistic jets. Of significant concern, however, has been that a strong toroidal field in the jet should be highly unstable to the non-axisymmetric helical kink (screw) $m=1$ mode leading to rapid disruption. In addition, a recent concern has been that the jet formation process itself may be unstable due to the accretion of non-dipolar magnetic fields. We describe large-scale fully three-dimensional global general relativistic magnetohydrodynamic simulations of rapidly rotating, accreting black holes producing jets. We study both the stability of the jet as it propagates and the stability of the jet formation process during accretion of dipolar and quadrupolar fields. For our dipolar model, despite strong non-axisymmetric disc turbulence, the jet reaches Lorentz factors of $Gammasim 10$ with opening half-angle $theta_jsim 5^circ$ at $10^3$ gravitational radii without significant disruption or dissipation with only mild substructure dominated by the $m=1$ mode. On the contrary, our quadrupolar model does not produce a steady relativistic ($Gammagtrsim 3$) jet due to mass-loading of the polar regions caused by unstable polar fields. Thus, if produced, relativistic jets are roughly stable structures and may reach up to an external shock with strong magnetic fields. We discuss the astrophysical implications of the accreted magnetic geometry playing such a significant role in relativistic jet formation, and we outline avenues for future work.

Three-Dimensional Simulations of Magnetized Thin Accretion Disks around Black Holes: Stress in the Plunging Region

We describe three-dimensional general relativistic magnetohydrodynamic simulations of a geometrically thin accretion disk around a non-spinning black hole. The disk has a thickness $h/rsim0.05-0.1$ over the radial range $(2-20)GM/c^2$. In steady state, the specific angular momentum profile of the inflowing magnetized gas deviates by less than 2% from that of the standard thin disk model of Novikov & Thorne (1973). Also, the magnetic torque at the radius of the innermost stable circular orbit (ISCO) is only $sim2%$ of the inward flux of angular momentum at this radius. Both results indicate that magnetic coupling across the ISCO is relatively unimportant for geometrically thin disks.