BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jonathan C. McKinney

Jonathan C. McKinney appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEastro-ph.GAgr-qcastro-ph.COastro-phastro-ph.IMastro-ph.SR

Catalog footprint

What is connected

33works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

33 published item(s)

preprint2022arXiv

PATOKA: Simulating Electromagnetic Observables of Black Hole Accretion

The Event Horizon Telescope (EHT) has released analyses of reconstructed images of horizon-scale millimeter emission near the supermassive black hole at the center of the M87 galaxy. Parts of the analyses made use of a large library of synthetic black hole images and spectra, which were produced using numerical general relativistic magnetohydrodynamics fluid simulations and polarized ray tracing. In this article, we describe the PATOKA pipeline, which was used to generate the Illinois contribution to the EHT simulation library. We begin by describing the relevant accretion systems and radiative processes. We then describe the details of the three numerical codes we use, iharm, ipole, and igrmonty, paying particular attention to differences between the current generation of the codes and the originally published versions. Finally, we provide a brief overview of simulated data as produced by PATOKA and conclude with a discussion of limitations and future directions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Effects of Spin on High-Energy Radiation from Accreting Black Holes

Observations of jets in X-ray binaries show a correlation between radio power and black hole spin. This correlation, if confirmed, points towards the idea that relativistic jets may be powered by the rotational energy of black holes. In order to examine this further, we perform general-relativistic radiative transport calculations on magnetically arrested accretion flows, which are known to produce powerful jets via the Blandford-Znajek (BZ) mechanism. We find that the X-ray and gamma-ray emission strongly depend on spin and inclination angle. Surprisingly, the high-energy power does not show the same dependence on spin as the BZ jet power, but instead can be understood as a redshift effect. In particular, photons observed perpendicular to the spin axis suffer little net redshift until originating from close to the horizon. Such observers see deeper into the hot, dense, highly-magnetized inner disk region. This effect is largest for rapidly rotating black holes due to a combination of frame dragging and decreasing horizon radius. While the X-ray emission is dominated by the near horizon region, the near-infrared radiation originates at larger radii. Therefore, the ratio of X-ray to near-infrared power is an observational signature of black hole spin.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Efficiency of Thin Magnetically-Arrested Disks Around Black Holes

The radiative and jet efficiencies of thin magnetized accretion disks around black holes (BHs) are affected by BH spin and the presence of a magnetic field that, when strong, could lead to large deviations from Novikov-Thorne (NT) thin disk theory. To seek the maximum deviations, we perform general relativistic magnetohydrodynamic (GRMHD) simulations of radiatively efficient thin (half-height $H$ to radius $R$ of $H/R\approx 0.10$) disks around moderately rotating BHs with $a/M=0.5$. First, our simulations, each evolved for more than $70,000r_g/c$ (gravitational radius $r_g$ and speed of light $c$), show that large-scale magnetic field readily accretes inward even through our thin disk and builds-up to the magnetically-arrested disk (MAD) state. Second, our simulations of thin MADs show the disk achieves a radiative efficiency of $η_{\rm r}\approx 15\%$ (after estimating photon capture), which is about twice the NT value of $η_{\rm r}\sim 8\%$ for $a/M=0.5$ and gives the same luminosity as a NT disk with $a/M\approx 0.9$. Compared to prior simulations with $\lesssim 10\%$ deviations, our result of an $\approx 80\%$ deviation sets a new benchmark. Building on prior work, we are now able to complete an important scaling law which suggest that observed jet quenching in the high-soft state in BH X-ray binaries is consistent with an ever-present MAD state with a weak yet sustained jet.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Jet Signatures in the Spectra of Accreting Black Holes

Jets are observed as radio emission in active galactic nuclei and during the low/hard state in X-ray binaries (XRBs), but their contribution at higher frequencies has been uncertain. We study the dynamics of jets in XRBs using the general-relativistic magnetohydrodynamic code HARM. We calculate the high-energy spectra and variability properties using a general-relativistic radiative transport code based on grmonty. We find the following signatures of jet emission (i) a significant gamma-ray peak above ~10^22 Hz, (ii) a break in the optical/UV spectrum, with a change in luminosity from L ~ nu^0 to L ~ nu, followed by another break at higher frequencies where the spectrum roughly returns to L ~ nu^0, and (iii) a pronounced synchrotron peak near or below ~10^14 Hz indicates that a significant fraction of any observed X-ray emission originates in the jet. We investigate the variability during a large-scale magnetic field inversion in which the Blandford-Znajek (BZ) jet is quenched and a new transient hot reconnecting plasmoid is launched by the reconnecting field. The ratio of the gamma-rays to X-rays changes from L_gamma / L_X > 1 in the BZ jet to L_gamma / L_X < 1 during the launching of the transient plasmoid.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Efficiency of Super-Eddington Magnetically-Arrested Accretion

The radiative efficiency of super-Eddington accreting black holes (BHs) is explored for magnetically-arrested disks (MADs), where magnetic flux builds-up to saturation near the BH. Our three-dimensional general relativistic radiation magnetohydrodynamic (GRRMHD) simulation of a spinning BH (spin $a/M=0.8$) accreting at $\sim 50$ times Eddington shows a total efficiency $\sim 50\%$ when time-averaged and total efficiency $\gtrsim 100\%$ in moments. Magnetic compression by the magnetic flux near the rotating BH leads to a thin disk, whose radiation escapes via advection by a magnetized wind and via transport through a low-density channel created by a Blandford-Znajek (BZ) jet. The BZ efficiency is sub-optimal due to inertial loading of field lines by optically thick radiation, leading to BZ efficiency $\sim 40\%$ on the horizon and BZ efficiency $\sim 5\%$ by $r\sim 400r_g$ (gravitational radii) via absorption by the wind. Importantly, radiation escapes at $r\sim 400r_g$ with efficiency $η\approx 15\%$ (luminosity $L\sim 50L_{\rm Edd}$), similar to $η\approx 12\%$ for a Novikov-Thorne thin disk and beyond $η\lesssim 1\%$ seen in prior GRRMHD simulations or slim disk theory. Our simulations show how BH spin, magnetic field, and jet mass-loading affect the radiative and jet efficiencies of super-Eddington accretion.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Resolved Magnetic-Field Structure and Variability Near the Event Horizon of Sagittarius A*

Near a black hole, differential rotation of a magnetized accretion disk is thought to produce an instability that amplifies weak magnetic fields, driving accretion and outflow. These magnetic fields would naturally give rise to the observed synchrotron emission in galaxy cores and to the formation of relativistic jets, but no observations to date have been able to resolve the expected horizon-scale magnetic-field structure. We report interferometric observations at 1.3-millimeter wavelength that spatially resolve the linearly polarized emission from the Galactic Center supermassive black hole, Sagittarius A*. We have found evidence for partially ordered fields near the event horizon, on scales of ~6 Schwarzschild radii, and we have detected and localized the intra-hour variability associated with these fields.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Soft X-ray Temperature Tidal Disruption Events from Stars on Deep Plunging Orbits

One of the puzzles associated with tidal disruption event candidates (TDEs) is that there is a dichotomy between the color temperatures of ${\rm few}\times 10^4$~K for TDEs discovered with optical and UV telescopes, and the color temperatures of ${\rm few}\times 10^5 - 10^6$~K for TDEs discovered with X-ray satellites. Here we propose that high-temperature TDEs are produced when the tidal debris of a disrupted star self-intersects relatively close to the supermassive black hole, in contrast to the more distant self-intersection that leads to lower color temperatures. In particular, we note from simple ballistic considerations that greater apsidal precession in an orbit is the key to closer self-intersection. Thus larger values of $β$, the ratio of the tidal radius to the pericenter distance of the initial orbit, are more likely to lead to higher temperatures of more compact disks which are super-Eddington and geometrically and optically thick. For a given star and $β$, apsidal precession also increases for larger black hole masses, but larger black hole masses imply a lower temperature at the Eddington luminosity. Thus the expected dependence of the temperature on the mass of the black hole is non-monotonic. We find that in order to produce a soft X-ray temperature TDE, a deep plunging stellar orbit with $β> 3$ is needed and a black hole mass of $\lesssim 5\times 10^6 M_\odot$ is favored. Although observations of TDEs are comparatively scarce and are likely dominated by selection effects, it is encouraging that both expectations are consistent with current data.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Submillimeter Quasi-Periodic Oscillations in Magnetically Choked Accretion Flows Models of Sgr A*

High-frequency quasi-periodic oscillations (QPOs) appear in general-relativistic magnetohydrodynamic simulations of magnetically choked accretion flows around rapidly rotating black holes (BHs). We perform polarized radiative transfer calculations with ASTRORAY code to explore the manifestations of these QPOs for SgrA*. We construct a simulation-based model of a radiatively inefficient accretion flow and find model parameters by fitting the mean polarized source spectrum. The simulated QPOs have a total sub-mm flux amplitude up to 5% and a linearly polarized flux amplitude up to 2%. The oscillations reach high levels of significance 10-30σand high quality factors Q~5. The oscillation period T~100M~35min corresponds to the rotation period of the BH magnetosphere that produces a trailing spiral in resolved disk images. The total flux signal is significant over noise for all tested frequencies 87GHz, 230GHz, and 857GHz and inclination angles 10deg, 37deg, and 80deg. The non-detection in the 230GHz Sub-Millimeter Array light curve is consistent with a low signal level and a low sampling rate. The presence of sub-mm QPOs in SgrA* will be better tested with the Atacama Large Millimeter Array.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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'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 "magneto-spin alignment" 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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constraining the Accretion Flow in Sgr A* by General Relativistic Dynamical and Polarized Radiative Modeling

We briefly summarize the method of simulating Sgr A* polarized sub-mm spectra from the accretion flow and fitting the observed spectrum. The dynamical flow model is based on three-dimensional general relativistic magneto hydrodynamic simulations. Fully self-consistent radiative transfer of polarized cyclo-synchrotron emission is performed. We compile a mean sub-mm spectrum of Sgr A* and fit it with the mean simulated spectra. We estimate the ranges of inclination angle theta=42-75deg, mass accretion rate Mdot=(1.4-7.0)*10^{-8}Msun/yr, and electron temperature Te=(3-4)*10^{10}K at 6M. We discuss multiple caveats in dynamical modeling, which must be resolved to make further progress.

0 citations0 reviewsNo code linkedOpen access
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 "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 $τ\sim 70GM/c^3$ for $a/M\sim 0.9$ with coherence quality factors $Q\gtrsim 10$. [abridged]

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Radiative Models of Sagittarius A* and M87 from Relativistic MHD Simulations

Ongoing millimeter VLBI observations with the Event Horizon Telescope allow unprecedented study of the innermost portion of black hole accretion flows. Interpreting the observations requires relativistic, time-dependent physical modeling. We discuss the comparison of radiative transfer calculations from general relativistic MHD simulations of Sagittarius A* and M87 with current and future mm-VLBI observations. This comparison allows estimates of the viewing geometry and physical conditions of the Sgr A* accretion flow. The viewing geometry for M87 is already constrained from observations of its large-scale jet, but, unlike Sgr A*, there is no consensus for its millimeter emission geometry or electron population. Despite this uncertainty, as long as the emission region is compact, robust predictions for the size of its jet launching region can be made. For both sources, the black hole shadow may be detected with future observations including ALMA and/or the LMT, which would constitute the first direct evidence for a black hole event horizon.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Sagittarius A* Accretion Flow and Black Hole Parameters from General Relativistic Dynamical and Polarized Radiative Modeling

We obtain estimates of Sgr A* accretion flow and black hole parameters by fitting polarized sub-mm observations with spectra computed using three-dimensional (3D) general relativistic (GR) magnetohydrodynamical (MHD) (GRMHD) simulations. Observations are compiled from averages over many epochs from reports in 29 papers for estimating the mean fluxes Fnu, linear polarization (LP) fractions, circular polarization (CP) fractions, and electric vector position angles (EVPAs). GRMHD simulations are computed with dimensionless spins a_*=0,0.5,0.7,0.9,0.98 over a 20,000M time interval. We perform fully self-consistent GR polarized radiative transfer using our new code to explore the effects of spin a_*, inclination angle θ, position angle (PA), accretion rate Mdot, and electron temperature Te (Te is reported for radius 6M). By fitting the mean sub-mm fluxes and LP/CP fractions, we obtain estimates for these model parameters and determine the physical effects that could produce polarization signatures. Our best bet model has a_*=0.5, θ=75deg, PA=115deg, Mdot=4.6*10^{-8}M_Sun/year, and Te=3.1*10^10K at 6M. The sub-mm CP is mainly produced by Faraday conversion as modified by Faraday rotation, and the emission region size at 230GHz is consistent with the VLBI size of 37microas. Across all spins, model parameters are in the ranges θ=42deg-75deg, Mdot=(1.4-7.0)*10^{-8}M_Sun/year, and Te=(3-4)*10^10K. Polarization is found both to help differentiate models and to introduce new observational constraints on the effects of the magnetic field that might not be fit by accretion models so-far considered.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The size of the jet launching region in M87

The supermassive black hole candidate at the center of M87 drives an ultra-relativistic jet visible on kiloparsec scales, and its large mass and relative proximity allow for event horizon scale imaging with very long baseline interferometry at millimeter wavelengths (mm-VLBI). Recently, relativistic magneto-hydrodynamic (MHD) simulations of black hole accretion flows have proven capable of launching magnetically-dominated jets. We construct time-dependent disc/jet models of the innermost portion of the M87 nucleus by performing relativistic radiative transfer calculations from one such simulation. We identify two types of models, jet-dominated or disc/jet, that can explain the spectral properties of M87, and use them to make predictions for current and future mm-VLBI observations. The Gaussian source size for the favored sky orientation and inclination from observations of the large-scale jet is 33-44 microarcseconds (~4-6 Schwarzschild radii) on current mm-VLBI telescopes, very similar to existing observations of Sgr A*. The black hole shadow, direct evidence of an event horizon, should be visible in future measurements using baselines between Hawaii and Mexico. Both models exhibit variability at millimeter wavelengths with factor of ~2 amplitudes on year timescales. For the low inclination of M87, the counter-jet dominates the event horizon scale millimeter wavelength emission from the jet-forming region.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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 ($r\sim 10^{13}--10^{14}{\rm cm}$), by which the jet obtains a high Lorentz factor ($γ\sim 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 $γθ_j\sim 10--20 (for opening half-angle $θ_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]

0 citations0 reviewsNo code linkedOpen access
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%.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Measuring Black Hole Spin by the Continuum-Fitting Method: Effect of Deviations from the Novikov-Thorne Disc Model

The X-ray spectra of accretion discs of eight stellar-mass black holes have been analyzed to date using the thermal continuum fitting method, and the spectral fits have been used to estimate the spin parameters of the black holes. However, the underlying model used in this method of estimating spin is the general relativistic thin-disc model of Novikov & Thorne, which is only valid for razor-thin discs. We therefore expect errors in the measured values of spin due to inadequacies in the theoretical model. We investigate this issue by computing spectra of numerically calculated models of thin accretion discs around black holes, obtained via three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations. We apply the continuum fitting method to these computed spectra to estimate the black hole spins and check how closely the values match the actual spin used in the GRMHD simulations. We find that the error in the dimensionless spin parameter is up to about 0.2 for a non-spinning black hole, depending on the inclination. For black holes with spins of 0.7, 0.9 and 0.98, the errors are up to about 0.1, 0.03 and 0.01 respectively. These errors are comparable to or smaller than those arising from current levels of observational uncertainty. Furthermore, we estimate that the GRMHD simulated discs from which these error estimates are obtained correspond to effective disc luminosities of about 0.4-0.7 Eddington, and that the errors will be smaller for discs with luminosities of 0.3 Eddington or less, which are used in the continuum-fitting method. We thus conclude that use of the Novikov-Thorne thin-disc model does not presently limit the accuracy of the continuum-fitting method of measuring black hole spin.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Slowly balding black holes

The "no hair" theorem, a key result in General Relativity, states that an isolated black hole is defined by only three parameters: mass, angular momentum, and electric charge; this asymptotic state is reached on a light-crossing time scale. We find that the "no hair" theorem is not formally applicable for black holes formed from collapse of a rotating neutron star. Rotating neutron stars can self-produce particles via vacuum breakdown forming a highly conducting plasma magnetosphere such that magnetic field lines are effectively "frozen-in" the star both before and during collapse. In the limit of no resistivity, this introduces a topological constraint which prohibits the magnetic field from sliding off the newly-formed event horizon. As a result, during collapse of a neutron star into a black hole, the latter conserves the number of magnetic flux tubes $N_B = e Φ_\infty /(πc \hbar)$, where $Φ_\infty \approx 2 π^2 B_{NS} R_{NS}^3 /(P_{\rm NS} c)$ is the initial magnetic flux through the hemispheres of the progenitor and out to infinity. We test this theoretical result via three-dimensional general relativistic plasma simulations of rotating black holes that start with a neutron star dipole magnetic field with no currents initially present outside the event horizon. The black hole's magnetosphere subsequently relaxes to the split monopole magnetic field geometry with self-generated currents outside the event horizon. The dissipation of the resulting equatorial current sheet leads to a slow loss of the anchored flux tubes, a process that balds the black hole on long resistive time scales rather than the short light-crossing time scales expected from the vacuum "no-hair" theorem.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Thin Disk Theory with a Non-Zero Torque Boundary Condition and Comparisons with Simulations

We present an analytical solution for thin disk accretion onto a Kerr black hole that extends the standard Novikov-Thorne alpha-disk in three ways: (i) it incorporates nonzero stresses at the inner edge of the disk, (ii) it extends into the plunging region, and (iii) it uses a corrected vertical gravity formula. The free parameters of the model are unchanged. Nonzero boundary stresses are included by replacing the Novikov-Thorne no torque boundary condition with the less strict requirement that the fluid velocity at the innermost stable circular orbit is the sound speed, which numerical models show to be the correct behavior for luminosities below ~30% Eddington. We assume the disk is thin so we can ignore advection. Boundary stresses scale as alpha*h and advection terms scale as h^2 (where h is the disk opening angle (h=H/r)), so the model is self-consistent when h < alpha. We compare our solution with slim disk models and general relativistic magnetohydrodynamic disk simulations. The model may improve the accuracy of black hole spin measurements.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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'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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Parsec-scale Faraday Rotation Measures from General Relativistic MHD Simulations of Active Galactic Nuclei Jets

For the first time it has become possible to compare global 3D general relativistic magnetohydrodynamic (GRMHD) jet formation simulations directly to very-long baseline interferometric multi-frequency polarization observations of the pc-scale structure of active galactic nucleus (AGN) jets. Unlike the jet emission, which requires post hoc modeling of the non-thermal electrons, the Faraday rotation measures (RMs) depend primarily upon simulated quantities and thus provide a robust way in which to confront simulations with observations. We compute RM distributions of 3D global GRMHD jet formation simulations, with which we explore the dependence upon model and observational parameters, emphasizing the signatures of structures generic to the theory of MHD jets. With typical parameters, we find that it is possible to reproduce the observed magnitudes and many of the structures found in AGN jet RMs, including the presence of transverse RM gradients. In our simulations the RMs are generated within a smooth extension of the jet itself, containing ordered toroidally dominated magnetic fields. This results in a particular bilateral morphology that is unlikely to arise due to Faraday rotation in distant foreground clouds. However, critical to efforts to probe the Faraday screen will be resolving the transverse jet structure. Therefore, the RMs of radio cores may not be reliable indicators of the properties of the rotating medium. Finally, we are able to constrain the particle content of the jet, finding that at pc-scales AGN jets are electromagnetically dominated, with roughly 2% of the comoving energy in nonthermal leptons and much less in baryons.

0 citations0 reviewsNo code linkedOpen access
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)

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Submillimeter Bump in Sgr A* from Relativistic MHD Simulations

Recent high resolution observations of the Galactic center black hole allow for direct comparison with accretion disk simulations. We compare two-temperature synchrotron emission models from three dimensional, general relativistic magnetohydrodynamic simulations to millimeter observations of Sgr A*. Fits to very long baseline interferometry and spectral index measurements disfavor the monochromatic face-on black hole shadow models from our previous work. Inclination angles \le 20 degrees are ruled out to 3 σ. We estimate the inclination and position angles of the black hole, as well as the electron temperature of the accretion flow and the accretion rate, to be i=50+35-15 degrees, ξ=-23+97-22 degrees, T_e=(5.4 +/- 3.0)x10^10 K and Mdot=(5+15-2)x10^-9 M_sun / yr respectively, with 90% confidence. The black hole shadow is unobscured in all best fit models, and may be detected by observations on baselines between Chile and California, Arizona or Mexico at 1.3mm or .87mm either through direct sampling of the visibility amplitude or using closure phase information. Millimeter flaring behavior consistent with the observations is present in all viable models, and is caused by magnetic turbulence in the inner radii of the accretion flow. The variability at optically thin frequencies is strongly correlated with that in the accretion rate. The simulations provide a universal picture of the 1.3mm emission region as a small region near the midplane in the inner radii of the accretion flow, which is roughly isothermal and has ν/ν_c ~ 1-20, where ν_c is the critical frequency for thermal synchrotron emission.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

2-D simulations of FU Orionis disk outbursts

We have developed time-dependent models of FU Ori accretion outbursts to explore the physical properties of protostellar disks. Our two-dimensional, axisymmetric models incorporate full vertical structure with a new treatment of the radiative boundary condition for the disk photosphere. We find that FU Ori-type outbursts can be explained by a slow accumulation of matter due to gravitational instability. Eventually this triggers the magnetorotational instability, which leads to rapid accretion. The thermal instability is triggered in the inner disk but this instability is not necessary for the outburst. An accurate disk vertical structure, including convection, is important for understanding the outburst behavior. Large convective eddies develop during the high state in the inner disk. The models are in agreement with Spitzer IRS spectra and also with peak accretion rates and decay timescales of observed outbursts, though some objects show faster rise timescale. We also propose that convection may account for the observed mild-supersonic turbulence and the short-timescale variations of FU Orionis objects.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

General Relativistic Magnetohydrodynamic Simulations of Jet Formation and Large-Scale Propagation from Black Hole Accretion Systems

The formation and large-scale propagation of Poynting-dominated jets produced by accreting, rapidly rotating black hole systems are studied by numerically integrating the general relativistic magnetohydrodynamic equations of motion to follow the self-consistent interaction between accretion disks and black holes. This study extends previous similar work by studying jets till $t\approx 10^4GM/c^3$ out to $r\approx 10^4GM/c^2$, by which the jet is super- fast magnetosonic and moves at a lab-frame bulk Lorentz factor of $Γ\sim 10$ with a maximum terminal Lorentz factor of $Γ_\infty\lesssim 10^3$. The radial structure of the Poynting-dominated jet is piece-wise self-similar, and fits to flow quantities along the field line are provided. Beyond the \alf surface at $r\sim 10$--$100GM/c^2$, the jet becomes marginally unstable to (at least) current-driven instabilities. Such instabilities drive shocks in the jet that limit the efficiency of magnetic acceleration and collimation. These instabilities also induce jet substructure with $3\lesssimΓ\lesssim 15$. The jet is shown to only marginally satisfy the necessary and sufficient conditions for kink instability, so this may explain how astrophysical jets can extend to large distances without completely disrupting. At large distance, the jet angular structure is Gaussian-like (or uniform within the core with sharp exponential wings) with a half-opening angle of $\approx 5^\circ$ and there is an extended component out to $\approx 27^\circ$. Unlike in some hydrodynamic simulations, the environment is found to play a negligible role in jet structure, acceleration, and collimation as long as the ambient pressure of the surrounding medium is small compared to the magnetic pressure in the jet.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Black Hole Spin Evolution

We consider a subset of the physical processes that determine the spin j = a/M of astrophysical black holes. These include: (1) Initial conditions. Recent models suggest that the collapse of supermassive stars are likely to produce black holes with j ~ 0.7. (2) Major mergers. The outcome of a nearly equal mass black hole-black hole merger is not yet known, but we review the current best guesses and analytic bounds. (3) Minor mergers. We recover the result of Blandford & Hughes that accretion of small companions with isotropically distributed orbital angular momenta results in spindown, with j ~ M^{-7/3}. (4) Accretion. We present new results from fully relativistic magnetohydrodynamic accretion simulations. These show that, at least for one sequence of flow models, spin equilibrium (dj/dt = 0) is reached for j ~ 0.9, far less than the canonical value 0.998 of Thorne that was derived in the absence of MHD effects. This equilibrium value may not apply to all accretion flows, particularly thin disks. Nevertheless, it opens the possibility that black holes that have grown primarily through accretion are not maximally rotating.

0 citations0 reviewsNo code linkedOpen access