Researcher profile

Alexander Tchekhovskoy

Alexander Tchekhovskoy contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
40works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

40 published item(s)

preprint2026arXiv

Evidence for Low Universal Equilibrium Black Hole Spin in Luminous Magnetically Arrested Disks

Relativistic collimated outflows, or jets, provide a crucial mode of active galactic nucleus feedback. Although jets extract their energy from the black hole (BH) rotation, their effect on the BH spin is poorly understood. Because the spin controls radiative and mechanical BH feedback, lack of first-principles models for spin evolution limits our ability to interpret observations, including the recent LIGO-Virgo-KAGRA spin constraints. Particularly important are luminous disks, which rapidly grow and strongly torque their BHs. Jetless and weakly magnetized standard luminous disks spin up their BHs to near-maximum spin, $a_{eq,NT}=0.998$. However, sufficient large-scale vertical magnetic flux can cause the inner disk to enter a magnetically arrested disk (MAD) state, whose jets can efficiently extract BH rotational energy and significantly spin down the BH. Lowell et al. (2024) found that nonradiative, thick MADs spin down their BHs to very low $a_{eq,MAD}^{thick}=0.07$. Their analytic model predicted that luminous, thin MADs also spin down their BHs to low $a_{eq,MAD}^{thin}\sim0.3\text{-}0.5$. To test this prediction, we perform 3D general relativistic (radiation) magnetohydrodynamic (GR(R)MHD) simulations of MADs across a wide range of BH spin ($-0.9\le{}a\le0.99$) and disk thickness ($0.03\le{}h/r\le0.3$, which corresponds to Eddington ratio, $0.35\le{}\dot{m}/\dot{m}_{Edd}\le\infty$). We find that luminous, thin MADs ($0.03\le{}h/r\le0.1$) efficiently spin down their BHs to a low universal equilibrium spin value, $a_{eq,MAD}^{thin}\approx0.3$: a maximally spinning BH ($a=1$) spins down to $a=0.5$ after accreting just $25\%$ of its initial mass. Our results follow quadratic convergence, $a_{eq,MAD}^{fit}\simeq0.3-2.7(h/r)^2\to0.3$ as $h/r\to0$, which we attribute to the aggressive cooling that renders disk thermodynamics irrelevant and magnetic forces insensitive to thermal $h/r$.

preprint2022arXiv

Black hole flares: ejection of accreted magnetic flux through 3D plasmoid-mediated reconnection

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.

preprint2022arXiv

Black hole to photosphere: 3D GRMHD simulations of collapsars reveal wobbling and hybrid composition jets

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.

preprint2022arXiv

Bridging Bondi and Event Horizon Scales: 3D GRMHD Simulations Reveal X-Shaped Radio Galaxy Morphology

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

preprint2022arXiv

On the jet-ejecta interaction in 3D GRMHD simulations of binary neutron star merger aftermath

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$).

preprint2021arXiv

Black hole to breakout: 3D GRMHD simulations of collapsar jets reveal a wide range of transients

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.

preprint2021arXiv

Magnetohydrodynamic Simulations of Active Galactic Nucleus Disks and Jets

There is a broad consensus that accretion onto supermassive black holes and consequent jet formation power the observed emission from active galactic nuclei (AGNs). However, there has been less agreement about how jets form in accretion flows, their possible relationship to black hole spin, and how they interact with the surrounding medium. There have also been theoretical concerns about instabilities in standard accretion disk models and lingering discrepancies with observational constraints. Despite seemingly successful applications to X-ray binaries, the standard accretion disk model faces a growing list of observational constraints that challenge its application to AGNs. Theoretical exploration of these questions has become increasingly reliant on numerical simulations owing to the dynamic nature of these flows and the complex interplay between hydrodynamics, magnetic fields, radiation transfer, and curved spacetime. We conclude the following: The advent of general relativistic magnetohydrodynamics (MHD) simulations has greatly improved our understanding of jet production and its dependence on black hole spin. Simulation results show both disks and jets are sensitive to the magnetic flux threading the accretion flow as well as possible misalignment between the angular momentum of the accretion flow and the black hole spin. Radiation MHD simulations are providing new insights into the stability of luminous accretion flows and highlighting the potential importance of radiation viscosity, UV opacity from atoms, and spiral density waves in AGNs.

preprint2016arXiv

Three-dimensional Relativistic MHD Simulations of Active Galactic Nuclei Jets: Magnetic Kink Instability and Fanaroff-Riley Dichotomy

Energy deposition by active galactic nuclei jets into the ambient medium can affect galaxy formation and evolution, the cooling of gas flows at the centres of galaxy clusters, and the growth of the supermassive black holes. However, the processes that couple jet power to the ambient medium and determine jet morphology are poorly understood. For instance, there is no agreement on the cause of the well-known Fanaroff-Riley (FR) morphological dichotomy of jets, with FRI jets being shorter and less stable than FRII jets. We carry out global 3D magnetohydrodynamic simulations of relativistic jets propagating through the ambient medium. We show that the flat density profiles of galactic cores slow down and collimate the jets, making them susceptible to the 3D magnetic kink instability. We obtain a critical power, which depends on the galaxy core mass and radius, below which jets become kink-unstable within the core, stall, and inflate cavities filled with relativistically-hot plasma. Jets above the critical power stably escape the core and form powerful backflows. Thus, the kink instability controls the jet morphology and can lead to the FR dichotomy. The model-predicted dependence of the critical power on the galaxy optical luminosity agrees well with observations.

preprint2015arXiv

Ab-initio pulsar magnetosphere: the role of general relativity

It has recently been demonstrated that self-consistent particle-in-cell simulations of low-obliquity pulsar magnetospheres in flat spacetime show weak particle acceleration and no pair production near the poles. We investigate the validity of this conclusion in a more realistic spacetime geometry via general-relativistic particle-in-cell simulations of the aligned pulsar magnetospheres with pair formation. We find that the addition of frame-dragging effect makes local current density along the magnetic field larger than the Goldreich-Julian value, which leads to unscreened parallel electric fields and the ignition of a pair cascade. When pair production is active, we observe field oscillations in the open field bundle which could be related to pulsar radio emission. We conclude that general relativistic effects are essential for the existence of pulsar mechanism in low obliquity rotators.

preprint2015arXiv

Core shifts, magnetic fields and magnetization of extragalactic jets

We study the effect of radio-jet core shift, which is a dependence of the position of the jet radio core on the observational frequency. We derive a new method of measuring the jet magnetic field based on both the value of the shift and the observed radio flux, which complements the standard method that assumes equipartition. Using both methods, we re-analyse the blazar sample of Zamaninasab et al. We find that equipartition is satisfied only if the jet opening angle in the radio core region is close to the values found observationally, $\simeq$0.1--0.2 divided by the bulk Lorentz factor, $Γ_{\rm j}$. Larger values, e.g., $1/Γ_{\rm j}$, would imply magnetic fields much above equipartition. A small jet opening angle implies in turn the magnetization parameter of $\ll 1$. We determine the jet magnetic flux taking into account this effect. We find that the transverse-averaged jet magnetic flux is fully compatible with the model of jet formation due to BH spin energy extraction and the accretion being a magnetically arrested disc (MAD). We calculate the jet average mass-flow rate corresponding to this model and find it consists of a substantial fraction of the mass accretion rate. This suggests the jet composition with a large fraction of baryons. We also calculate the average jet power, and find it moderately exceeds the accretion power, $\dot M c^2$, reflecting BH spin energy extraction. We find our results for radio galaxies at low Eddington ratios are compatible with MADs but require a low radiative efficiency, as predicted by standard accretion models.

preprint2015arXiv

Electron Thermodynamics in GRMHD Simulations of Low-Luminosity Black Hole Accretion

Simple assumptions made regarding electron thermodynamics often limit the extent to which general relativistic magnetohydrodynamic (GRMHD) simulations can be applied to observations of low-luminosity accreting black holes. We present, implement, and test a model that self-consistently evolves an electron entropy equation and takes into account the effects of spatially varying electron heating and relativistic anisotropic thermal conduction along magnetic field lines. We neglect the back-reaction of electron pressure on the dynamics of the accretion flow. Our model is appropriate for systems accreting at $\ll 10^{-5}$ of the Eddington rate, so radiative cooling by electrons can be neglected. It can be extended to higher accretion rates in the future by including electron cooling and proton-electron Coulomb collisions. We present a suite of tests showing that our method recovers the correct solution for electron heating under a range of circumstances, including strong shocks and driven turbulence. Our initial applications to axisymmetric simulations of accreting black holes show that (1)~physically-motivated electron heating rates yield electron temperature distributions significantly different from the constant electron to proton temperature ratios assumed in previous work, with higher electron temperatures concentrated in the coronal region between the disc and the jet; (2)~electron thermal conduction significantly modifies the electron temperature in the inner regions of black hole accretion flows if the effective electron mean free path is larger than the local scale-height of the disc (at least for the initial conditions and magnetic field configurations we study). The methods developed in this work are important for producing more realistic predictions for the emission from accreting black holes such as Sagittarius A* and M87; these applications will be explored in future work.

preprint2015arXiv

Evolution of non-spherical pulsars with plasma-filled magnetospheres

Pulsars are famous for their rotational stability. Most of them steadily spin down and display a highly repetitive pulse shape. But some pulsars experience timing irregularities such as nulling, intermittency, mode changing and timing noise. As changes in the pulse shape are often correlated with timing irregularities, precession is a possible cause of these phenomena. Whereas pulsar magnetospheres are filled with plasma, most pulsar precession studies were carried out within the vacuum approximation and neglected the effects of magnetospheric currents and charges. Recent numerical simulations of plasma-filled pulsar magnetospheres provide us with a detailed quantitative description of magnetospheric torques exerted on the pulsar surface. In this paper, we present the study of neutron star evolution using these new torque expressions. We show that they lead to (1) much slower long-term evolution of pulsar parameters and (2) much less extreme solutions for these parameters than the vacuum magnetosphere models. To facilitate the interpretation of observed pulsar timing residuals, we derive an analytic model that (1) describes the time evolution of non-spherical pulsars and (2) translates the observed pulsar timing residuals into the geometrical parameters of the pulsar. We apply this model to two pulsars with very different temporal behaviours. For the pulsar B1828-11, we demonstrate that the timing residual curves allow two pulsar geometries: one with stellar deformation pointing along the magnetic axis and one along the rotational axis. For the Crab pulsar, we use the model show that the recent observation of its magnetic and rotational axes moving away from each other can be explained by precession.

preprint2015arXiv

Horizon-Scale Lepton Acceleration in Jets: Explaining the Compact Radio Emission in M87

It has now become clear that the radio jet in the giant elliptical galaxy M87 must turn on very close to the black hole. This implies the efficient acceleration of leptons within the jet at scales much smaller than feasible by the typical dissipative events usually invoked to explain jet synchrotron emission. Here we show that the stagnation surface, the separatrix between material that falls back into the black hole and material that is accelerated outward forming the jet, is a natural site of pair formation and particle acceleration. This occurs via an inverse-Compton pair catastrophe driven by unscreened electric fields within the charge-starved region about the stagnation surface and substantially amplified by a post-gap cascade. For typical estimates of the jet properties in M87, we find excellent quantitive agreement between the predicted relativistic lepton densities and those required by recent high-frequency radio observations of M87. This mechanism fails to adequately fill a putative jet from Sagittarius A* with relativistic leptons, which may explain the lack of an obvious radio jet in the Galactic center. Finally, this process implies a relationship between the kinetic jet power and the gamma-ray luminosity of blazars, produced during the post-gap cascade.

preprint2015arXiv

Jet and disk luminosities in tidal disruption events

Tidal disruption events (TDE) in which a star is devoured by a massive black hole at a galac- tic center pose a challenge to our understanding of accretion processes. Within a month the accretion rate reaches super-Eddington levels. It then drops gradually over a time scale of a year to sub-Eddington regimes. The initially geometrically thick disk becomes a thin one and eventually an ADAF at very low accretion rates. As such, TDEs explore the whole range of accretion rates and configurations. A challenging question is what the corresponding light curves of these events are. We explore numerically the disk luminosity and the conditions within the inner region of the disk using a fully general relativistic slim disk model. Those conditions determine the magnitude of the magnetic field that engulfs the black hole and this, in turn, determines the Blandford-Znajek jet power. We estimate this power in two different ways and show that they are self-consistent. We find, as expected earlier from analytic argu- ments (Krolik & Piran 2012), that neither the disk luminosity nor the jet power follows the accretion rate throughout the disruption event. The disk luminosity varies only logarithmi- cally with the accretion rate at super-Eddington luminosities. The jet power follows initially the accretion rate but remains a constant after the transition from super- to sub- Eddington. At lower accretion rates at the end of the MAD phase the disk becomes thin and the jet may stop altogether. These new estimates of the jet power and disk luminosity that do not simply follow the mass fallback rate should be taken into account when searching for TDEs and analysing light curves of TDE candidates. Identification of some of the above mentioned transitions may enable us to estimate better TDE parameters.

preprint2015arXiv

On The Efficiency of Jet Production in Radio Galaxies

The mechanisms that produce and power relativistic jets are fundamental open questions in black hole (BH) astrophysics. In order to constrain these mechanisms, we analyze the energy efficiency of jet production &#39;eta&#39; based on archival Chandra observations of 27 nearby, low-luminosity active galactic nuclei. We obtain &#39;eta&#39; as the ratio of the jet power, inferred from the energetics of jet powered X-ray emitting cavities, to the BH mass accretion rate Mdot_BH. The standard assumption in estimating Mdot_BH is that all the gas from the Bondi radius &#39;r_B&#39; makes it down to the BH. It is now clear, however, that only a small fraction of the gas reaches the hole. To account for this effect, we use the standard disk mass-loss scaling, Mdot(r) \propto (r/r_B)^s Mdot_Bondi. This leads to much lower values of Mdot_BH and higher values of &#39;eta&#39; than in previous studies. If hot accretion flows are characterized by 0.5<s<0.6 -- on the lower end of recent theoretical and observational studies -- then dynamically-important magnetic fields near rapidly spinning BHs are necessary to account for the high eta~100-300 per cent in the sample. Moreover, values of s>0.6 are essentially ruled out, or there would be insufficient energy to power the jets. We discuss the implications of our results for the distribution of massive BH spins and the possible impact of a significant extra cold gas supply on our estimates.

preprint2015arXiv

Relativistic MHD simulations of core-collapse GRB jets: 3D instabilities and magnetic dissipation

Relativistic jets naturally occur in astrophysical systems that involve accretion onto compact objects, such as core collapse of massive stars in gamma-ray bursts (GRBs) and accretion onto supermassive black holes in active galactic nuclei (AGN). It is generally accepted that these jets are powered electromagnetically, by the magnetised rotation of a central compact object. However, how they produce the observed emission and survive the propagation for many orders of magnitude in distance without being disrupted by current-driven non-axisymmetric instabilities is the subject of active debate. We carry out time-dependent 3D relativistic magnetohydrodynamic simulations of relativistic, Poynting flux dominated jets. The jets are launched self-consistently by the rotation of a strongly magnetised central compact object. This determines the natural degree of azimuthal magnetic field winding, a crucial factor that controls jet stability. We find that the jets are susceptible to two types of instability: (i) a global, external kink mode that grows on long time scales and causes the jets to bodily bend sideways. Whereas this mode does not cause jet disruption over the simulated distances, it substantially reduces jet propagation speed. We show, via an analytic model, that the growth of the external kink mode depends on the slope of the ambient medium density profile. In flat density distributions characteristic of galactic cores, an AGN jet may stall, whereas in stellar envelopes the external kink weakens as the jet propagates outward; (ii) a local, internal kink mode that grows over short time scales and causes small-angle magnetic reconnection and conversion of about half of jet electromagnetic energy flux into heat. Based on the robustness and energetics of the internal kink mode, we suggest that this instability is the main dissipation mechanism responsible for powering GRB prompt emission.

preprint2015arXiv

Three-dimensional Analytical Description of Magnetised Winds from Oblique Pulsars

Rotating neutron stars, or pulsars, are plausibly the source of power behind many astrophysical systems, such as gamma-ray bursts, supernovae, pulsar wind nebulae and supernova remnants. In the past several years, 3D numerical simulations made it possible to compute pulsar spindown luminosity from first principles and revealed that oblique pulsar winds are more powerful than aligned ones. However, what causes this enhanced power output of oblique pulsars is not understood. In this work, using time-dependent 3D magnetohydrodynamic (MHD) and force-free simulations, we show that, contrary to the standard paradigm, the open magnetic flux, which carries the energy away from the pulsar, is laterally non-uniform. We argue that this non-uniformity is the primary reason for the increased luminosity of oblique pulsars. To demonstrate this, we construct simple analytic descriptions of aligned and orthogonal pulsar winds and combine them to obtain an accurate 3D description of the pulsar wind for any obliquity. Our approach describes both the warped magnetospheric current sheet and the smooth variation of pulsar wind properties outside of it. We find that generically the magnetospheric current sheet separates plasmas that move at mildly relativistic velocities relative to each other. This suggests that the magnetospheric reconnection is a type of driven, rather than free, reconnection. The jump in magnetic field components across the current sheet decreases with increasing obliquity, which could be a mechanism that reduces dissipation in near-orthogonal pulsars. Our analytical description of the pulsar wind can be used for constructing models of pulsar gamma-ray emission, pulsar wind nebulae, and magnetar-powered core-collapse gamma-ray bursts and supernovae.

preprint2014arXiv

Global simulations of axisymmetric radiative black hole accretion disks in general relativity with a sub-grid magnetic dynamo

We present a sub-grid model that emulates the magnetic dynamo operating in magnetized accretion disks. We have implemented this model in the general relativisic radiation magnetohydrodynamic (GRRMHD) code \koral, using results from local shearing sheet simulations of the magnetorotational instability to fix the parameters of the dynamo. With the inclusion of this dynamo, we are able to run 2D axisymmetric GRRMHD simulations of accretion disks for arbitrarily long times. The simulated disks exhibit sustained turbulence, with the poloidal and toroidal magnetic field components driven towards a state similar to that seen in 3D studies. Using this dynamo code, we present a set of long-duration global simulations of super-Eddington, optically-thick disks around non-spinning and spinning black holes. Super-Eddington disks around non-rotating black holes exhibit a surprisingly large efficiency, $η\approx0.04$, independent of the accretion rate, where we measure efficiency in terms of the total energy output, both radiation and mechanical, flowing out to infinity. Super-Eddington disks around spinning black holes are even more efficient, and appear to extract black hole rotational energy through a process similar to the Blandford-Znajek mechanism. All the simulated models are characterized by highly super-Eddington radiative fluxes collimated along the rotation axis. We also present a set of simulations that were designed to have Eddington or slightly sub-Eddington accretion rates ($\dot{M} \lesssim 2\dot M_{\rm Edd}$). None of these models reached a steady state. Instead, the disks collapsed as a result of runaway cooling, presumably because of a thermal instability.

preprint2014arXiv

Magnetic Flux of Progenitor Stars Sets Gamma-ray Burst Luminosity and Variability

Long-duration gamma-ray bursts (GRBs) are thought to come from the core-collapse of Wolf-Rayet stars. Whereas their stellar masses $M_*$ have a rather narrow distribution, the population of GRBs is very diverse, with gamma-ray luminosities $L_γ$ spanning several orders of magnitude. This suggests the existence of a &#34;hidden&#34; stellar variable whose burst-to-burst variation leads to a spread in $L_γ$. Whatever this hidden variable is, its variation should not noticeably affect the shape of GRB lightcurves, which display a constant luminosity (in a time-average sense) followed by a sharp drop at the end of the burst seen with Swift/XRT. We argue that such a hidden variable is progenitor star&#39;s large-scale magnetic flux. Shortly after the core collapse, most of stellar magnetic flux accumulates near the black hole (BH) and remains there. The flux extracts BH rotational energy and powers jets of roughly a constant luminosity, $L_j$. However, once BH mass accretion rate $\dot M$ falls below $\sim L_j/c^2$, the flux becomes dynamically important and diffuses outwards, with the jet luminosity set by the rapidly declining mass accretion rate, $L_j\sim \dot M c^2$. This provides a potential explanation for the sharp end of GRBs and the universal shape of their lightcurves. During the GRB, gas infall translates spatial variation of stellar magnetic flux into temporal variation of $L_j$. We make use of the deviations from constancy in $L_j$ to perform stellar magnetic flux &#34;tomography&#34;. Using this method, we infer the presence of magnetised tori in the outer layers of progenitor stars for GRB 920513 and GRB 940210.

preprint2014arXiv

Tidal Disruption and Magnetic Flux Capture: Powering a Jet from a Quiescent Black Hole

The transient Swift J1644+57 is believed to have been produced by an unlucky star wandering too close to a supermassive black hole (BH) leading to a tidal disruption event. This unusual flare displayed highly super-Eddington X-ray emission which likely originated in a relativistic, collimated jet. This presents challenges to modern accretion and jet theory as upper limits of prior BH activity, which we obtain from the radio afterglow of this event, imply that both the pre-disruption BH and stellar magnetic fluxes fall many orders of magnitude short of what is required to power the observed X-ray luminosity. We argue that a pre-existing, &#34;fossil&#34; accretion disc can contain a sufficient reservoir of magnetic flux and that the stellar debris stream is capable of dragging this flux into the BH. To demonstrate this, we perform local, 3D magnetohydrodynamic simulations of the disc--stream interaction and demonstrate that the interface between the two is unstable to mixing. This mixing entrains a sufficient amount of fossil disc magnetic flux into the infalling stellar debris to power the jet. We argue that the interaction with the fossil disc can have a pronounced effect on the structure and dynamics of mass fallback and likely the resulting transient. Finally, we describe possible ramifications of these interactions on unresolved problems in tidal disruption dynamics, in particular, the efficiency of debris circularization, and effects of the disruption on the preexisting black hole system. Animations online: http://goo.gl/T84tLs

preprint2014arXiv

Time evolution of pulsar obliquity angle from 3D simulations of magnetospheres

The rotational period of isolated pulsars increases over time due to the extraction of angular momentum by electromagnetic torques. These torques also change the obliquity angle $α$ between the magnetic and rotational axes. Although actual pulsar magnetospheres are plasma-filled, the time evolution of $α$ has mostly been studied for vacuum pulsar magnetospheres. In this work, we self-consistently account for the plasma effects for the first time by analysing the results of time-dependent 3D force-free and magnetohydrodynamic simulations of pulsar magnetospheres. We show that if a neutron star is spherically symmetric and is embedded with a dipolar magnetic moment, the pulsar evolves so as to minimise its spin-down luminosity: both vacuum and plasma-filled pulsars evolve toward the aligned configuration ($α=0$). However, they approach the alignment in qualitatively different ways. Vacuum pulsars come into alignment exponentially fast, with $α\propto \exp(-t/τ)$ and $τ\sim$ spindown timescale. In contrast, we find that plasma-filled pulsars align much more slowly, with $α\propto (t/τ)^{-1/2}$. We argue that the slow time evolution of obliquity of plasma-filled pulsars can potentially resolve several observational puzzles, including the origin of normal pulsars with periods of $\sim1$ second, the evidence that oblique pulsars come into alignment over a timescale of $\sim 10^7$ years, and the observed deficit, relative to an isotropic obliquity distribution, of pulsars showing interpulse emission.

preprint2013arXiv

Energy Extraction from Spinning Black Holes via Relativistic Jets

It has for long been an article of faith among astrophysicists that black hole spin energy is responsible for powering the relativistic jets seen in accreting black holes. Two recent advances have strengthened the case. First, numerical general relativistic magnetohydrodynamic simulations of accreting spinning black holes show that relativistic jets form spontaneously. In at least some cases, there is unambiguous evidence that much of the jet energy comes from the black hole, not the disk. Second, spin parameters of a number of accreting stellar-mass black holes have been measured. For ballistic jets from these systems, it is found that the radio luminosity of the jet correlates with the spin of the black hole. This suggests a causal relationship between black hole spin and jet power, presumably due to a generalized Penrose process.

preprint2013arXiv

Semi-implicit scheme for treating radiation under M1 closure in general relativistic conservative fluid dynamics codes

A numerical scheme is described for including radiation in multi-dimensional general-relativistic conservative fluid dynamics codes. In this method, a covariant form of the M1 closure scheme is used to close the radiation moments, and the radiative source terms are treated semi-implicitly in order to handle both optically thin and optically thick regimes. The scheme has been implemented in a conservative general relativistic radiation hydrodynamics code KORAL. The robustness of the code is demonstrated on a number of test problems, including radiative relativistic shock tubes, static radiation pressure supported atmosphere, shadows, beams of light in curved spacetime, and radiative Bondi accretion. The advantages of M1 closure relative to other approaches such as Eddington closure and flux-limited diffusion are discussed, and its limitations are also highlighted.

preprint2013arXiv

Swift J1644+57 gone MAD: the case for dynamically-important magnetic flux threading the black hole in a jetted tidal disruption event

The unusual transient Swift J1644+57 likely resulted from a collimated relativistic jet powered by accretion onto a massive black hole (BH) following the tidal disruption (TD) of a star. Several mysteries cloud the interpretation of this event: (1) extreme flaring and `plateau&#39; shape of the X-ray/gamma-ray light curve during the first 10 days after the gamma-ray trigger; (2) unexpected rebrightening of the forward shock radio emission months after trigger; (3) no obvious evidence for jet precession, despite misalignment typically expected between the angular momentum of the accretion disk and BH; (4) recent abrupt shut-off in jet X-ray emission after 1.5 years. Here we show that all of these seemingly disparate mysteries are naturally resolved by one assumption: the presence of strong magnetic flux Phi threading the BH. Initially, Phi is weak relative to high fall-back mass accretion rate, Mdot, and the disk and jets precess about the BH axis = our line of sight. As Mdot drops, Phi becomes dynamically important and leads to a magnetically-arrested disk (MAD). MAD naturally aligns disk and jet axis along the BH spin axis, but only after a violent rearrangement phase (jet wobbling). This explains the erratic light curve at early times and the lack of precession at later times. We use our model for Swift J1644+57 to constrain BH and disrupted star properties, finding that a solar-mass main sequence star disrupted by a relatively low mass, M~10^5-10^6 Msun, BH is consistent with the data, while a WD disruption (though still possible) is disfavored. The magnetic flux required to power Swift J1644+57 is too large to be supplied by the star itself, but it could be collected from a quiescent `fossil&#39; accretion disk present in the galactic nucleus prior to the TD. The presence (lack of) of such a fossil disk could be a deciding factor in what TD events are accompanied by powerful jets.[abridged]

preprint2013arXiv

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 $z\sim 6-8$ is achievable even with rapidly rotating BHs and ten solar mass BH seeds.

preprint2013arXiv

Time-Dependent 3D Magnetohydrodynamic Pulsar Magnetospheres: Oblique Rotators

The current state of the art in pulsar magnetosphere modeling assumes the force-free limit of magnetospheric plasma. This limit retains only partial information about plasma velocity and neglects plasma inertia and temperature. We carried out time-dependent 3D relativistic magnetohydrodynamic (MHD) simulations of oblique pulsar magnetospheres that improve upon force-free by retaining the full plasma velocity information and capturing plasma heating in strong current layers. We find rather low levels of magnetospheric dissipation, with less than 10% of pulsar spindown energy dissipated within a few light cylinder radii, and the MHD spindown that is consistent with that in force-free. While oblique magnetospheres are qualitatively similar to the rotating split-monopole force-free solution at large radii, we find substantial quantitative differences with the split-monopole, e.g., the luminosity of the pulsar wind is more equatorially concentrated than the split-monopole at high obliquities, and the flow velocity is modified by the emergence of reconnection flow directed into the current sheet.

preprint2013arXiv

Transient jet formation and state transitions from large-scale magnetic reconnection in black hole accretion discs

Magnetically arrested accretion discs (MADs), where the magnetic pressure in the inner disc is dynamically important, provide an alternative mechanism for regulating accretion to what is commonly assumed in black hole systems. We show that a global magnetic field inversion in the MAD state can destroy the jet, significantly increase the accretion rate, and move the effective inner disc edge in to the marginally stable orbit. Reconnection of the MAD field in the inner radii launches a new type of transient outflow containing hot plasma generated by magnetic dissipation. This transient outflow can be as powerful as the steady magnetically-dominated Blandford-Znajek jet in the MAD state. The field inversion qualitatively describes many of the observational features associated with the high luminosity hard to soft state transition in black hole X-ray binaries: the jet line, the transient ballistic jet, and the drop in rms variability. These results demonstrate that the magnetic field configuration can influence the accretion state directly, and hence the magnetic field structure is an important second parameter in explaining observations of accreting black holes across the mass and luminosity scales.

preprint2012arXiv

Alignment of Magnetized Accretion Disks and Relativistic Jets with Spinning Black Holes

Accreting black holes (BHs) produce intense radiation and powerful relativistic jets, which are affected by the BH&#39;s spin magnitude and direction. While thin disks might align with the BH spin axis via the Bardeen-Petterson effect, this does not apply to jet systems with thick disks. We used fully three-dimensional general relativistic magnetohydrodynamical simulations to study accreting BHs with various BH spin vectors and disk thicknesses with magnetic flux reaching saturation. Our simulations reveal a &#34;magneto-spin alignment&#34; mechanism that causes magnetized disks and jets to align with the BH spin near BHs and further away to reorient with the outer disk. This mechanism has implications for the evolution of BH mass and spin, BH feedback on host galaxies, and resolved BH images for SgrA* and M87.

preprint2012arXiv

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 ($ϕ$-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 &#34;magnetically choked accretion flow&#34; (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 $τ\sim 70GM/c^3$ for $a/M\sim 0.9$ with coherence quality factors $Q\gtrsim 10$. [abridged]

preprint2012arXiv

General Relativistic Modeling of Magnetized Jets from Accreting Black Holes

Recent advances in general relativistic magnetohydrodynamic modeling of jets offer unprecedented insights into the inner workings of accreting black holes that power the jets in active galactic nuclei (AGN) and other accretion systems. I will present the results of recent studies that determine spin-dependence of jet power and discuss the implications for the AGN radio loud/quiet dichotomy and recent observations of high jet power in a number of AGN.

preprint2012arXiv

On the Spin-Down of Intermittent Pulsars

Magnetospheres of pulsars are thought to be filled with plasma, and variations in plasma supply can affect both pulsar emission properties and spin-down rates. A number of recently discovered &#34;intermittent&#34; pulsars switch between two distinct states: an &#34;on&#34;, radio-loud state, and an &#34;off&#34;, radio-quiet state. Spin-down rates in the two states differ by a large factor, $\sim 1.5-2.5$, which is not easily understood in the context of current models. In this Letter we present self-consistent numerical solutions of &#34;on&#34; and &#34;off&#34; states of intermittent pulsar magnetospheres. We model the &#34;on&#34; state as a nearly ideal force-free magnetosphere with abundant magnetospheric plasma supply. The lack of radio emission in the &#34;off&#34; state is associated with plasma supply disruption that results in lower plasma density on the open field lines. We model the &#34;off&#34; state using nearly vacuum conditions on the open field lines and nearly ideal force-free conditions on the closed field lines, where plasma can remain trapped even in the absence of pair production. The toroidal advection of plasma in the closed zone in the &#34;off&#34; state causes spin-downs that are a factor of $\sim 2$ higher than vacuum values, and we naturally obtain a range of spin-down ratios between the &#34;on&#34; and &#34;off&#34; states, $\sim 1.2-2.9$, which corresponds to a likely range of pulsar inclination angles of $30{-}90^\circ$. We consider the implications of our model to a number of poorly understood but possibly related pulsar phenomena, including nulling, timing noise, and rotating radio transients.

preprint2012arXiv

Prograde and Retrograde Black Holes: Whose Jet is More Powerful?

The outflow efficiency (eta) from black hole (BH) accretion disc systems is known to depend upon both the BH spin (a) and the amount of large-scale magnetic flux threading the BH and disc. Semi-analytical flux-trapping models suggest retrograde BHs should trap much more large-scale magnetic flux near the BH leading to much higher eta than for prograde BHs. We self-consistently determine the amount of large-scale magnetic flux trapped by rapidly spinning (a = -0.9 and 0.9) BHs using global 3D time-dependent non-radiative general relativistic magnetohydrodynamic simulations of thick (h/r ~ 0.3-0.6) discs. We find that BH-trapped flux builds up until it is strong enough to disrupt the inner accretion disc. Contrary to prior flux-trapping models, which do not include the back-reaction of magnetic flux on the disc, our simulations show prograde BHs trap more magnetic flux, leading to about 3 times higher eta than retrograde BHs for |a| = 0.9. Both spin orientations can produce highly efficient jets, eta ~ 100%, with increasing eta for increasing disc thickness. The similarity of eta for prograde and retrograde BHs makes it challenging to infer the sign of BH spin based on jet energetics alone.

preprint2011arXiv

Constraints on Cold Magnetized Shocks in Gamma-Ray Bursts

We consider a model in which the ultra-relativistic jet in a gamma-ray burst (GRB) is cold and magnetically accelerated. We assume that the energy flux in the outflowing material is partially thermalized via internal shocks or a reverse shock, and we estimate the maximum amount of radiation that could be produced in such magnetized shocks. We compare this estimate with the available observational data on prompt gamma-ray emission in GRBs. We find that, even with highly optimistic assumptions, the magnetized jet model is radiatively too inefficient to be consistent with observations. One way out is to assume that much of the magnetic energy in the post-shock, or even pre-shock, jet material is converted to particle thermal energy by some unspecified process, and then radiated. This can increase the radiative efficiency sufficiently to fit observations. Alternatively, jet acceleration may be driven by thermal pressure rather than magnetic fields. In this case, which corresponds to the traditional fireball model, sufficient prompt GRB emission could be produced either from shocks at a large radius or from the jet photosphere closer to the center.

preprint2011arXiv

Efficient Generation of Jets from Magnetically Arrested Accretion on a Rapidly Spinning Black Hole

We describe global, 3D, time-dependent, non-radiative, general-relativistic, magnetohydrodynamic simulations of accreting black holes (BHs). The simulations are designed to transport a large amount of magnetic flux to the center, more than the accreting gas can force into the BH. The excess magnetic flux remains outside the BH, impedes accretion, and leads to a magnetically arrested disc. We find powerful outflows. For a BH with spin parameter a = 0.5, the efficiency with which the accretion system generates outflowing energy in jets and winds is eta ~ 30%. For a = 0.99, we find eta ~ 140%, which means that more energy flows out of the BH than flows in. The only way this can happen is by extracting spin energy from the BH. Thus the a = 0.99 simulation represents an unambiguous demonstration, within an astrophysically plausible scenario, of the extraction of net energy from a spinning BH via the Penrose-Blandford-Znajek mechanism. We suggest that magnetically arrested accretion might explain observations of active galactic nuclei with apparent eta ~ few x 100%.

preprint2011arXiv

Resistive Solutions for Pulsar Magnetospheres

The current state of the art in the modeling of pulsar magnetospheres invokes either the vacuum or force-free limits for the magnetospheric plasma. Neither of these limits can simultaneously account for both the plasma currents and the accelerating electric fields that are needed to explain the morphology and spectra of high-energy emission from pulsars. To better understand the structure of such magnetospheres, we combine accelerating fields and force-free solutions by considering models of magnetospheres filled with resistive plasma. We formulate Ohm&#39;s Law in the minimal velocity fluid frame and construct a family of resistive solutions that smoothly bridges the gap between the vacuum and the force-free magnetosphere solutions. The spin-down luminosity, open field line potential drop, and the fraction of open field lines all transition between the vacuum and force-free values as the plasma conductivity varies from zero to infinity. For fixed inclination angle, we find that the spin-down luminosity depends linearly on the open field line potential drop. We consider the implications of our resistive solutions for the spin down of intermittent pulsars and sub-pulse drift phenomena in radio pulsars.

preprint2010arXiv

Black Hole Spin and the Radio Loud/Quiet Dichotomy of Active Galactic Nuclei

Radio loud active galactic nuclei (AGN) are on average 1000 times brighter in the radio band compared to radio quiet AGN. We investigate whether this radio loud/quiet dichotomy can be due to differences in the spin of the central black holes that power the radio-emitting jets. Using general relativistic magnetohydrodynamic simulations, we construct steady state axisymmetric numerical models for a wide range of black hole spins (dimensionless spin parameter 0.1 <= a <= 0.9999 and a variety of jet geometries. We assume that the total magnetic flux through the black hole horizon at radius r_H(a) is held constant. If the black hole is surrounded by a thin accretion disk, we find that the total black hole power output depends approximately quadratically on the angular frequency of the hole, P \propto Ω_H^2 \propto (a/r_H)^2. We conclude that, in this scenario, differences in the black hole spin can produce power variations of only a few tens at most. However, if the disk is thick such that the jet subtends a narrow solid angle around the polar axis, then the power dependence becomes much steeper, P \propto Ω_H^4 or even \propto Ω_H^6. Power variations of 1000 are then possible for realistic black hole spin distributions. We derive an analytic solution that accurately reproduces the steeper scaling of jet power with Ω_H, and we provide a numerical fitting formula that reproduces all our simulation results. We discuss other physical effects that might contribute to the observed radio loud/quiet dichotomy of AGN.

preprint2010arXiv

Magnetohydrodynamic Simulations of Gamma-Ray Burst Jets: Beyond the Progenitor Star

Achromatic breaks in afterglow light curves of gamma-ray bursts (GRBs) arise naturally if the product of the jet&#39;s Lorentz factor γand opening angle Θ_j satisfies (γΘ_j) >> 1 at the onset of the afterglow phase, i.e., soon after the conclusion of the prompt emission. Magnetohydrodynamic (MHD) simulations of collimated GRB jets generally give (γΘ_j) <~ 1, suggesting that MHD models may be inconsistent with jet breaks. We work within the collapsar paradigm and use axisymmetric relativistic MHD simulations to explore the effect of a finite stellar envelope on the structure of the jet. Our idealized models treat the jet-envelope interface as a collimating rigid wall, which opens up outside the star to mimic loss of collimation. We find that the onset of deconfinement causes a burst of acceleration accompanied by a slight increase in the opening angle. In our fiducial model with a stellar radius equal to 10^4.5 times that of the central compact object, the jet achieves an asymptotic Lorentz factor γ~ 500 far outside the star and an asymptotic opening angle Θ_j ~ 0.04 rad ~ 2 deg, giving (γΘ_j) ~ 20. These values are consistent with observations of typical long-duration GRBs, and explain the occurrence of jet breaks. We provide approximate analytic solutions that describe the numerical results well.

preprint2010arXiv

Simulations of Magnetized Disks Around Black Holes: Effects of Black Hole Spin, Disk Thickness, and Magnetic Field Geometry

The standard general relativistic model of a razor-thin accretion disk around a black hole, developed by Novikov & Thorne (NT) in 1973, assumes the shear stress vanishes at the radius of the innermost stable circular orbit (ISCO) and that, outside the ISCO, the shear stress is produced by an effective turbulent viscosity. However, astrophysical accretion disks are not razor-thin, it is uncertain whether the shear stress necessarily vanishes at the ISCO, and the magnetic field, which is thought to drive turbulence in disks, may contain large-scale structures that do not behave like a simple local scalar viscosity. We describe three-dimensional general relativistic magnetohydrodynamic simulations of accretion disks around black holes with a range of spin parameters, and we use the simulations to assess the validity of the NT model. Our fiducial initial magnetic field consists of multiple (alternating polarity) poloidal field loops whose shape is roughly isotropic in the disk in order to match the isotropic turbulence expected in the poloidal plane. For a thin disk with an aspect ratio |h/r| ~ 0.07 around a non-spinning black hole, we find a decrease in the accreted specific angular momentum of 2.9% relative to the NT model and an excess luminosity from inside the ISCO of 3.5%. The deviations in the case of spinning black holes are of the same order. In addition, the deviations decrease with decreasing |h/r|. We therefore conclude that magnetized thin accretion disks in x-ray binaries in the thermal/high-soft spectral state ought to be well-described by the NT model, especially at luminosities below 30% of Eddington where we expect a very small disk thickness |h/r| <~ 0.05. We also discuss how the stress and the luminosity inside the ISCO depend on the assumed initial magnetic field geometry. (abridged)