No Arabic abstract
We have performed 2.5-dimensional general relativistic magnetohydrodynamic (MHD) simulations of collapsars including a rotating black hole. This paper is an extension of our previous paper (Mizuno et al. 2004). The current calculation focuses on the effect of black hole rotation using general relativistic MHD with simplified microphysics. Initially, we assume that the core collapse is failed in this star. A few solar mass rotating black hole is inserted by hand into the calculation. We consider two cases, the co-rotating case and counter-rotating case with respect to the black hole rotation. Although the counter-rotating case may be unrealistic for collapsar, we perform as a maximally dragging case of magnetic field. The simulation results show the formation of a disk-like structure and the generation of a jet-like outflow near the central black hole. The jet-like outflow propagates outwardly with the twisted magnetic field and becomes collimated. We have found that the jets are generated and accelerated mainly by the magnetic field. The total jet velocity in the rotating black hole case is comparable to that of the non-rotating black hole case (Mizuno et al. 2004), 0.3c. When the rotation of the black hole is faster, the magnetic field is twisted strongly owing to the frame-dragging effect. The magnetic energy stored by the twisting magnetic field is converted to kinetic energy of the jet directly rather than propagating as an Alfven wave. Thus, as the rotation of the black hole becomes faster, the poloidal velocity of the jet becomes faster. In the rapidly rotating black hole case the jet-like outflow can be produced by the frame dragging effect only through the twisting of magnetic field even if there is no stellar rotation.
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]
General relativistic magnetohydrodynamic (GRMHD) flows along magnetic fields threading a black hole can be divided into inflow and outflow parts, according to the result of the competition between the black hole gravity and magneto-centrifugal forces along the field line. Here we present the first self-consistent, semi-analytical solution for a cold, Poynting flux-dominated (PFD) GRMHD flow, which passes all four critical (inner and outer, Alfven and fast magnetosonic) points along a parabolic streamline. By assuming that the dominating (electromagnetic) component of the energy flux per flux tube is conserved at the surface where the inflow and outflow are separated, the outflow part of the solution can be constrained by the inflow part. The semi-analytical method can provide fiducial and complementary solutions for GRMHD simulations around the rotating black hole, given that the black hole spin, global streamline, and magnetizaion (i.e., a mass loading at the inflow/outflow separation) are prescribed. For reference, we demonstrate a self-consistent result with the work by McKinney in a quantitative level.
We have studied the dynamics of an equal-mass magnetized neutron-star binary within a resistive magnetohydrodynamic (RMHD) approach in which the highly conducting stellar interior is matched to an electrovacuum exterior. Because our analysis is aimed at assessing the modifications introduced by resistive effects on the dynamics of the binary after the merger and through to collapse, we have carried out a close comparison with an equivalent simulation performed within the traditional ideal magnetohydrodynamic approximation. We have found that there are many similarities between the two evolutions but also one important difference: the survival time of the hyper massive neutron star increases in a RMHD simulation. This difference is due to a less efficient magnetic-braking mechanism in the resistive regime, in which matter can move across magnetic-field lines, thus reducing the outward transport of angular momentum. Both the RMHD and the ideal magnetohydrodynamic simulations carried here have been performed at higher resolutions and with a different grid structure than those in previous work of ours [L. Rezzolla, B. Giacomazzo, L. Baiotti, J. Granot, C. Kouveliotou, and M. A. Aloy, Astrophys. J. Letters 732, L6 (2011)], but confirm the formation of a low-density funnel with an ordered magnetic field produced by the black hole--torus system. In both regimes the magnetic field is predominantly toroidal in the highly conducting torus and predominantly poloidal in the nearly evacuated funnel. Reconnection processes or neutrino annihilation occurring in the funnel, none of which we model, could potentially increase the internal energy in the funnel and launch a relativistic outflow, which, however, is not produced in these 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.
The classical, relativistic thin-disk theory of Novikov and Thorne (NT) predicts a maximum accretion efficiency of 40% for an optically thick, radiatively efficient accretion disk around a maximally spinning black hole (BH). However, when a strong magnetic field is introduced to numerical simulations of thin disks, large deviations in efficiency are observed, in part due to mass and energy carried by jets and winds launched by the disk or BH spin. The total efficiency of accretion can be significantly enhanced beyond that predicted by NT but it has remained unclear how the radiative component is affected. In order to study the effect of a dynamically relevant large-scale magnetic field on radiatively efficient accretion, we have performed numerical 3D general relativistic - radiative - magnetohydroynamic (GRRMHD) simulations of a disk with scale height to radius ratio of $H/R~0.1$ around a moderately spinning BH (a=0.5) using the code HARMRAD. Our simulations are fully global and allow us to measure the jet, wind, and radiative properties of a magnetically arrested disk (MAD) that is kept thin via self-consistent transport of energy by radiation using the M1 closure scheme. Our fiducial disk is MAD out to a radius of ~16R_g and the majority of the total ~13% efficiency of the accretion flow is carried by a magnetically driven wind. We find that the radiative efficiency is slightly suppressed compared to NT, contrary to prior MAD GRMHD simulations with an ad hoc cooling function, but it is unclear how much of the radiation and thermal energy trapped in the outflows could ultimately escape.