Three-Dimensional General Relativistic Radiation Magnetohydrodynamical Simulation of Super-Eddington Accretion, using a new code HARMRAD with M1 Closure

Black hole (BH) accretion flows and jets are dynamic hot relativistic magnetized plasma flows whose radiative opacity can significantly affect flow structure and behavior. We describe a numerical scheme, tests, and an astrophysically relevant application using the M1 radiation closure within a new three-dimensional (3D) general relativistic (GR) radiation (R) magnetohydrodynamics (MHD) massively parallel code called HARMRAD. Our 3D GRRMHD simulation of super-Eddington accretion (about $20$ times Eddington) onto a rapidly rotating BH (dimensionless spin $j=0.9375$) shows sustained non-axisymmemtric disk turbulence, a persistent electromagnetic jet driven by the Blandford-Znajek effect, and a total radiative output consistently near the Eddington rate. The total accretion efficiency is of order $20%$, the large-scale electromagnetic jet efficiency is of order $10%$, and the total radiative efficiency that reaches large distances remains low at only order $1%$. However, the radiation jet and the electromagnetic jet both emerge from a geometrically beamed polar region, with super-Eddington isotropic equivalent luminosities. Such simulations with HARMRAD can enlighten the role of BH spin vs. disks in launching jets, help determine the origin of spectral and temporal states in x-ray binaries, help understand how tidal disruption events (TDEs) work, provide an accurate horizon-scale flow structure for M87 and other active galactic nuclei (AGN), and isolate whether AGN feedback is driven by radiation or by an electromagnetic, thermal, or kinetic wind/jet. For example, the low radiative efficiency and weak BH spin-down rate from our simulation suggest that BH growth over cosmological times to billions of solar masses by redshifts of $zsim 6-8$ is achievable even with rapidly rotating BHs and ten solar mass BH seeds.

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.