Person:

Narayan, Ramesh

A supermassive black hole in the nucleus of an elliptical galaxy at the centre of a cool-core group or cluster of galaxies is immersed in hot gas. Bondi accretion should occur at a rate determined by the properties of the gas at the Bondi radius and the mass of the black hole. X-ray observations of massive nearby elliptical galaxies, including M87 in the Virgo cluster, indicate a Bondi accretion rate (M_B) which roughly matches the total kinetic power of the jets, suggesting that there is a tight coupling between the jet power and the mass accretion rate. While the Bondi model considers non-rotating gas, it is likely that the external gas has some angular momentum, which previous studies have shown could decrease the accretion rate drastically. We investigate here the possibility that viscosity acts at all radii to transport angular momentum outwards so that the accretion inflow proceeds rapidly and steadily. The situation corresponds to a giant advection-dominated accretion flow (ADAF) which extends from beyond the Bondi radius down to the black hole. We find solutions of the ADAF equations in which the gas accretes at just a factor of a few less than (M_B). These solutions assume that the atmosphere beyond the Bondi radius rotates with a sub-Keplerian velocity and that the viscosity parameter is large, α≥ 0.1, both of which are reasonable for the problem at hand. The infall time of the ADAF solutions is no more than a few times the free-fall time. Thus, the accretion rate at the black hole is closely coupled to the surrounding gas, enabling tight feedback to occur. We show that jet powers of a few per cent of (M_{B}c^{2}) are expected if either a fraction of the accretion power is channelled into the jet or the black hole spin energy is tapped by a strong magnetic field pressed against the black hole by the pressure of the accretion flow. We discuss the Bernoulli parameter of the flow, the role of convection and the possibility that these as well as magnetohydrodynamic effects may invalidate the model. If the latter comes to pass, it would imply that the rough agreement between observed jet powers and the Bondi accretion rate is a coincidence and jet power is determined by factors other than the mass accretion rate.

Long-duration gamma-ray bursts (GRBs) require an engine capable of driving a jet of plasma to ultrarelativistic bulk Lorentz factors of up to several hundred and into narrow opening angles of a few degrees. We use global axisymmetric stationary solutions of magnetically dominated (force-free) ultrarelativistic jets to test whether the popular magnetic-driving paradigm can generate the required Lorentz factors and opening angles. Our global solutions are obtained via time-dependent relativistic ideal magnetodynamical numerical simulations which follow the jet from the central engine to beyond six orders of magnitude in radius. Our model is primarily motivated by the collapsar model, in which a jet is produced by a spinning black hole or neutron star and then propagates through a massive stellar envelope. We find that the size of the pre-supernova progenitor star and the radial profile of pressure inside the star determine the terminal Lorentz factor and opening angle of the jet. At the radius where the jet breaks out of the star, our well-motivated fiducial model generates a Lorentz factor γ ∼ 400 and a half-opening angle θj ∼ 2◦, consistent with observations of many longduration GRBs. Other models with slightly different parameters give γ in the range 100–5000 and θj from 0. ◦1 to 10◦, thus reproducing the range of properties inferred for GRB jets. A potentially observable feature of some of our solutions is that the maximum Poynting flux in the jet is found at θ ∼ θj with the jet power concentrated in a hollow cone, while the maximum in the Lorentz factor occurs at an angle θ substantially smaller than θj also in a hollow cone. We derive approximate analytical formulae for the radial and angular distribution of γ and the radial dependence of θj . These formulae reproduce the simulation results and allow us to predict the outcome of models beyond those simulated. We also briefly discuss applications to active galactic nuclei, X-ray binaries and short-duration GRBs.

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 centre, 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 η≈ 30 per cent. For a = 0.99, we find η≈ 140 per cent, 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 η≈ few × 100 per cent.

Radio loud active galactic nuclei (AGNs) are on average 1000 times brighter in the radio band compared to radio quiet AGNs. We investigate whether this radio loud/quiet dichotomy can be due to differences in the spin of the central black holes (BHs) that power the radio-emitting jets. Using general relativistic magnetohydrodynamic simulations, we construct steady state axisymmetric numerical models for a wide range of BH spins (dimensionless spin parameter 0.1 a 0.9999) and a variety of jet geometries. We assume that the total magnetic flux through the BH horizon at radiusrH(a) is held constant. If the BH is surrounded by a thin accretion disk, we find that the total BH power output depends approximately quadratically on the angular frequency of the hole, P ∝ Ω2 H ∝ (a/rH) 2. We conclude that, in this scenario, differences in the BH 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 ∝ Ω4 H or even ∝ Ω6 H. Power variations of 1000 are then possible for realistic BH 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 AGNs.

We consider a two-parameter family of cylindrical force-free equilibria, modeled to match numerical simulations of relativistic force-free jets. We study the linear stability of these equilibria, assuming a rigid impenetrable wall at the outer cylindrical radius Rj. Equilibria in which the Lorentz factor γ (R) increases monotonically with increasing radius R are found to be stable. On the other hand, equilibria in which γ (R) reaches a maximum value at an intermediate radius and then declines to a smaller value γj at Rj are unstable. A feature of these unstable equilibria is that poloidal field line curvature plays a prominent role in maintaining transverse force balance. The most rapidly growing mode is an m = 1 kink instability which has a growth rate ∼ (0.4/γj )(c/Rj ). The e-folding length of the equivalent convected instability is ∼2.5γjRj . For a typical jet with an opening angle θj ∼ few/γj , the mode amplitude grows only weakly with increasing distance from the base of the jet. The growth is much slower than one might expect from a naive application of the Kruskal–Shafranov stability criterion.

A new general relativistic radiation magnetohydrodynamical code KORAL is described, which employs the M1 scheme to close the radiation moment equations. The code has been successfully verified against a number of tests. Axisymmetric simulations of super-critical magnetized accretion on a non-rotating black hole (a∗ = 0.0) and a spinning black hole (a∗ = 0.9) are presented. The accretion rates in the two models are M˙ ≈ 100 ÷ 200M˙ Edd. These first general relativistic simulations of super-critical black hole accretion are potentially relevant to tidal disruption events and hyper-accreting supermassive black holes in the early universe. Both simulated models are optically and geometrically thick, and have funnels through which energy escapes in the form of relativistic gas, Poynting flux and radiative flux. The jet is significantly more powerful in the a∗ = 0.9 run. The net energy outflow rate in the two runs correspond to efficiencies of 5% (a∗ = 0) and 33% (a∗ = 0.9), as measured with respect to the mass accretion rate at the black hole. These efficiencies agree well with those measured in previous simulations of non-radiative geometrically thick disks. Furthermore, in the a∗ = 0.9 run, the outflow power appears to originate in the spinning black hole, suggesting that the associated physics is again similar in non-radiative and super-critical accretion flows. While the two simulations are efficient in terms of total energy outflow, both runs are radiatively inefficient. Their luminosities are only ∼ 1 − 10LEdd, which corresponds to a radiative efficiency ∼ 0.1%. Interestingly, most of the radiative luminosity emerges through the funnels, which subtend a very small solid angle. Therefore, measured in terms of a local radiative flux, the emitted radiation is highly super-Eddington.

We describe a set of simulations of super-critical accretion onto a non-rotating supermassive BH. The accretion flow is radiation pressure dominated and takes the form of a geometrically thick disk with twin low-density funnels around the rotation axis. For accretion rates & 10M˙ Edd, there is sufficient gas in the funnel to make this region optically thick. Radiation from the disk first flows into the funnel, after which it accelerates the optically thick funnel gas along the axis. The resulting jet is baryon-loaded and has a terminal density-weighted velocity ≈ 0.3c. Much of the radiative luminosity is converted into kinetic energy by the time the escaping gas becomes optically thin. These jets are not powered by black hole rotation or magnetic driving, but purely by radiation. Their characteristic beaming angle is ∼ 0.2 radians. For an observer viewing down the axis, the isotropic equivalent luminosity of total energy is as much as 1048 erg s−1 for a 107M BH accreting at 103 Eddington. Therefore, energetically, the simulated jets are consistent with observations of the most powerful tidal disruption events, e.g., Swift J1644. The jet velocity is, however, too low to match the Lorentz factor γ > 2 inferred in J1644. There is no such conflict in the case of other tidal disruption events. Since favorably oriented observers see isotropic equivalent luminosities that are highly superEddington, the simulated models can explain observations of ultra-luminous X-ray sources, at least in terms of luminosity and energetics, without requiring intermediate mass black holes. The spectrum remains to be worked out. Finally, since the simulated jets are baryon-loaded and have mildly relativistic velocities, they match well the jets observed in SS433. The latter are, however, more collimated than the simulated jets. This suggests that, even if magnetic fields are not important for acceleration, they may perhaps still play a role in confining the jet.

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, η ≈ 0.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 (M˙ . 2M˙ 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.

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, “fossil” 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

Unconfined relativistic outflows from rotating, magnetized compact objects are often well modeled by assuming that the field geometry is approximately a split-monopole at large radii. Earlier work has indicated that such an unconfined flow has an inefficient conversion of magnetic energy to kinetic energy. This has led to the conclusion that ideal magnetohydrodynamical (MHD) processes fail to explain observations of, e.g., the Crab pulsar wind at large radii where energy conversion appears efficient. In addition, as a model for astrophysical jets, the monopole field geometry has been abandoned in favor of externally confined jets since the latter appeared to be generically more efficient jet accelerators. We perform time-dependent axisymmetric relativistic MHD simulations in order to find steady-state solutions for a wind from a compact object endowed with a monopole field geometry. Our simulations follow the outflow for 10 orders of magnitude in distance from the compact object, which is large enough to study both the initial “acceleration zone” of the magnetized wind as well as the asymptotic “coasting zone.” We obtain the surprising result that acceleration is actually efficient in the polar region, which develops a jet despite not being confined by an external medium. Our models contain jets that have sufficient energy to account for moderately energetic long and short gamma-ray burst (GRB) events (∼1051–1052 erg), collimate into narrow opening angles (opening half-angle θj ≈ 0.03 rad), become matter-dominated at large radii (electromagnetic energy flux per unit matter energy flux σ < 1), and move at ultrarelativistic Lorentz factors (γj ∼ 200 for our fiducial model). The simulated jets have γj θj ∼ 5–15, so they are in principle capable of generating “achromatic jet breaks” in GRB afterglow light curves. By defining a “causality surface” beyond which the jet cannot communicate with a generalized “magnetic nozzle” near the axis of rotation, we obtain approximate analytical solutions for the Lorentz factor that fit the numerical solutions well. This allows us to extend our results to monopole wind models with arbitrary magnetization. Overall, our results demonstrate that the production of ultrarelativistic jets is a more robust process than previously thought.