A Unified Model for Black Hole X-Ray Binary Jets?

We have recently put forward a ‘unified’ semi-empirical model for the coupling between accretion and jet production in galactic black hole X-ray binaries. In this paper, we summarise this model and briefly discuss relevant considerations that have arisen since its publication.     

Large Scale Jets in Microquasars

Relativistic jets are now believed to be a fairly ubiquitous property of accreting compact objects, and are intimately coupled with the accretion history. Associated with rapid changes in the accretion states of the binary systems, ejections of relativistic plasma can be observed at radio frequencies on timescale of weeks before becoming undetectable. However, recent observations point to long-term effects of these ejecta on the interstellar medium with the formation of large-scale relativistic jets around binary systems. In this paper, we review the observations of these large-scale structures in microquasars, highlighting their contributions at high energies.     

Jet simulations extending radially self-similar magnetohydrodynamics models

We perform a numerical simulation of magnetohydrodynamics (MHD) radially self-similar jets, whose prototype is the Blandford & Payne analytical example. The final steady state that is reached is valid close to the rotation axis and also at large distances above the disc where the classical analytical model fails to provide physically acceptable solutions. The outflow starts with a subslow magnetosonic speed, which subsequently crosses all relevant MHD critical points and corresponding magnetosonic separatrix surfaces. The characteristics are plotted together with the Mach cones and the superfast magnetosonic outflow satisfies MHD causality. The final solution remains close enough to the analytical one, which is thus shown to be topologically stable and robust for various boundary conditions.     

Radio-loudness of active galaxies and the black hole evolution

Active galactic nuclei (AGNs) form two distinct sequences on the radio-loudness–Eddington ratio plane. The ‘upper’ sequence contains radio selected AGNs, the ‘lower’ sequence is composed mainly of optically selected AGNs. The sequences mark the upper bounds for the radio-loudness of two distinct populations of AGNs, hosted, respectively, by elliptical and disk galaxies. Both sequences show the same dependence of the radio-loudness on the Eddington ratio (an increase with decreasing Eddington ratio), which suggests that another parameter in addition to the accretion rate must play a role in determining the efficiency of jet production in AGNs. We speculate that this additional parameter is the spin of the black hole, assuming that black holes in giant elliptical galaxies have (on average) much larger spins than black holes in disc galaxies. Possible evolutionary scenarios leading to such a spin dichotomy are discussed. The galaxy-morphology related radio dichotomy breaks down at high accretion rates where the dominant fraction of luminous quasars being hosted by giant ellipticals is radio-quiet. This indicates that the production of powerful jets at high accretion rates is in most cases suppressed and, in analogy to X-ray binary systems (XRB) during high and very high states, may be intermittent. Such intermittency can be caused by switches between two different accretion modes, assuming that only during one of them an outflow from the central engine is sufficiently collimated to form a relativistic jet.     

Ring-Shaped Jets in Gamma-Ray Bursts

When the axis of a gamma-ray burst (GRB) does not coincide with the spin axis of its source, there may result a ring-shaped jet. Using some refined jet dynamics, we calculate multi-wavelength afterglow light curves for such ring-shaped jets. In the R-band we find an obvious break in the afterglow light curve due to the beaming effect and the break is affected by many parameters, such as the electron energy fraction ξe, the magnetic energy fraction ξ2B, the width of ring A0 and the medium number density n. The overall light curve can be divided into three power-law stages, I.e., an ultra-relativistic stage, an after-break stage and a deep Newtonian stage. For each stage the power-law index is larger in the ring-shaped jet than in the corresponding conical jet.     

Numerical Simulations and Astrophysical Applications of Laboratory Jets at Omega

We have conducted experiments on the Omega laser at the University of Rochester that have produced jets of supersonic Ti impacting and being deflected by a ball of high density plastic. These mm-sized jets of dense plasma are highly complex, have large Reynolds numbers, and, given sufficient time and shear, should produce a fully turbulent flow. The experiments are diagnosed with a point-projection backlighter, resulting in a single image per shot. Simulations of the 3D hydrodynamics capture the large-scale features of the experimental data fairly well while missing some of the smaller scale turbulent-like phenomena. This is to be expected given the limited characterization of the targets as well as the finite resolution of the 3D simulations. If Euler scaling holds, these experiments should model larger astrophysical jets in objects such as HH 110 where an outflow can be seen colliding with a molecular cloud. However, Euler scaling demands that not only the isothermal internal Mach numbers of the two systems be similar but also that any dissipative mechanisms, such as radiative cooling or viscous dissipation, be of equal importance relative to each other. Similar equations of state are required as well. We discuss such issues in the context of these experiments and simulations.     

Energization of interstellar media and cosmic ray production by jets from X-ray binaries

Drawing on recent estimates of the power of jets from X-ray binary systems as a function of X-ray luminosity, combined with improved estimates of the relevant  log( N )–log( L _X)  luminosity functions, we calculate the total energy input to the interstellar medium (ISM) from these objects. The input of kinetic energy to the ISM via jets is dominated by those of the black hole systems, in contrast to the radiative input, which is dominated by accreting neutron stars. Summing the energy input from black hole jets L _J in the Milky Way, we find that it is likely to correspond to ≥1 per cent of L _SNe, the time-averaged kinetic luminosity of supernovae, and ≥5 per cent of L _CR, the cosmic ray luminosity. Given uncertainties in jet power estimates, significantly larger contributions are possible. Furthermore, in elliptical galaxies with comparable distributions of low mass X-ray binaries, but far fewer supernovae, the ratio   L _J/ L _SNe  is likely to be larger by a factor of ∼5. We conclude that jets from X-ray binaries may be an important, distributed, source of kinetic energy for the ISM in the form of relativistic shocks, and as a result are likely to be a major source of cosmic rays.     

标准薄盘与径移主导吸积流的连接

用标准的龙格-库塔法求解描述黑洞吸积流的基本方程组，没有引入任何附加的能量转移机制和外边界条件，结果表明在强粘滞情况下标准薄盘可以向径移主导吸积流转变。     

Laboratory Experiments of Stellar Jets from the Perspective of an Observer

It has been two decades since astronomers first discovered that accretion disks around young stars drive highly collimated supersonic jets. Thanks to concerted efforts to understand emission line ratios from jets, we know that velocity variations dominate the heating within these flows, and motions in stellar jets, now observed in real time, are primarily radial. The fluid dynamics of the cooling zones can be complex, with interacting shocks, clumps, and instabilities that could benefit from insights into the physics that only experiments can provide. Recent laboratory experiments have reproduced jets with velocities and Mach numbers similar to those within stellar jets, and the field seems poised to make significant advances by connecting observations and theories with experiments. This article points out several aspects of stellar jets that might be clarified by such experiments.     

Jets from black hole X-ray binaries: testing, refining and extending empirical models for the coupling to X-rays

In this paper we study the relation of radio emission to X-ray spectral and variability properties for a large sample of black hole X-ray binary systems. This is done to test, refine and extend – notably into the timing properties – the previously published 'unified model' for the coupling of accretion and ejection in such sources. In 14 outbursts from 11 different sources we find that in every case the peak radio flux, on occasion directly resolved into discrete relativistic ejections, is associated with the bright hard to soft state transition near the peak of the outburst. We also note the association of the radio flaring with periods of X-ray flaring during this transition in most, but not all, of the systems. In the soft state, radio emission is in nearly all cases either undetectable or optically thin, consistent with the suppression of the core jet in these states and 'relic' radio emission from interactions of previously ejected material and the ambient medium. However, these data cannot rule out an intermittent, optically thin, jet in the soft state. In attempting to associate X-ray timing properties with the ejection events we find a close, but not exact, correspondence between phases of very low integrated X-ray variability and such ejections. In fact the data suggest that there is not a perfect one-to-one correspondence between the radio, X-ray spectral or X-ray timing properties, suggesting that they may be linked simply as symptoms of the underlying state change and not causally to one another. We further study the sparse data on the reactivation of the jet during the transition back to the hard state in decay phase of outbursts, and find marginal evidence for this in one case only. In summary we find no strong evidence against the originally proposed model, confirming and extending some aspects of it with a much larger sample, but note that several aspects remain poorly tested.     

X-ray Flashes from Off-axis Nonuniform Jets

It has been widely believed that the outflows in gamma-ray bursts are jetted and some jets may have structures like ε(θ) ∝ θ^-k. We check the possibility that X-ray flashes come from such jets. Both qualitative and quantitative analyses have shown that this model can reproduce most of the observational features of both X-ray flashes and gamma-ray bursts. Using the usual parameters of gamma-ray bursts, we have carried out numerical calculations for both uniform and nonuniform jets, of their fluxes, spectra and peak energies. It seems that nonuniform jets are more appropriate to these observational properties than uniform jets.We have also shown that in our model the observational ratio of gamma-ray bursts to X-ray flashes is about a few units.     

A toy model for magnetic extraction of energy from black hole accretion disk

A toy model for magnetic extraction of energy from black hole (BH) accretion disk is discussed by considering the restriction of the screw instability to the magnetic field configuration. Three mechanisms of extracting energy magnetically are involved. (1) The Blandford–Znajek (BZ) process is related to the open magnetic field lines connecting the BH with the astrophysical load; (2) the magnetic coupling (MC) process is related to the closed magnetic field lines connecting the BH with its surrounding disk; and (3) a new scenario (henceforth the DL process) for extracting rotational energy from the disk is related to the open field lines connecting the disk with the astrophysical load. The expressions for the electromagnetic powers and torques are derived by using the equivalent circuits corresponding to the above energy mechanisms. It turns out that the DL power is comparable with the BZ and MC powers as the BH spin approaches unity. The radiation from a quasi-steady thin disk is discussed in detail by applying the conservation laws of mass, energy and angular momentum to the regions corresponding to the MC and DL processes. In addition, the poloidal currents and the current densities in BH magnetosphere are calculated by using the equivalent circuits.     

On the local magnetorotational instability of astrophysical jets

We present the local linear stability analysis of rotating jets confined by a toroidal magnetic field. Under the thin flux tube approximation, we derive the equation of motion for slender magnetic flux tubes. In addition to the terms responsible for the conventional instability of the toroidal magnetic field, a term related to the magnetic buoyancy and a term corresponding to the differential rotation become relevant for the stability properties. We find that the rigid rotation stabilizes while the differential rotational destabilizes the jet in a way similar to the Balbus–Hawley instability. Within the frame of our local analysis, we find that if the azimuthal velocity is of the order of or higher than the Alfvén azimuthal speed, the rigidly rotating part of the jet interior can be completely stabilized, while the strong shearing instability operates in the transition layer between the rotating jet interior and the external medium. This can explain the limb-brightening effect observed in several jets. However, it is still possible to find jet equilibria that are stable all across the jet, even in the presence of differential rotation. We discuss observational consequences of these results.     

GRMHD/RMHD simulations & stability of magnetized spine-sheath relativistic jets

A new general relativistic magnetohydrodynamics (GRMHD) code “RAISHIN” used to simulate jet generation by rotating and non-rotating black holes with a geometrically thin Keplarian accretion disk finds that the jet develops a spine-sheath structure in the rotating black hole case. Spine-sheath structure and strong magnetic fields significantly modify the Kelvin-Helmholtz (KH) velocity shear driven instability. The RAISHIN code has been used in its relativistic magnetohydrodynamic (RMHD) configuration to study the effects of strong magnetic fields and weakly relativistic sheath motion, c/2, on the KH instability associated with a relativistic, γ=2.5, jet spine-sheath interaction. In the simulations sound speeds up to and Alfvén wave speeds up to ∼0.56c are considered. Numerical simulation results are compared to theoretical predictions from a new normal mode analysis of the RMHD equations. Increased stability of a weakly magnetized system resulting from c/2 sheath speeds and stabilization of a strongly magnetized system resulting from c/2 sheath speeds is found.     

A magnetohydrodynamical model for the formation of episodic jets

Episodic ejection of plasma blobs has been observed in many black hole systems. While steady, continuous jets are believed to be associated with large-scale open magnetic fields, what causes the episodic ejection of blobs remains unclear. Here by analogy with the coronal mass ejection on the Sun, we propose a magnetohydrodynamical model for episodic ejections from black holes associated with the closed magnetic fields in an accretion flow. Shear and turbulence of the accretion flow deform the field and result in the formation of a flux rope in the disc corona. Energy and helicity are accumulated and stored until a threshold is reached. The system then loses its equilibrium and the flux rope is thrust outward by the magnetic compression force in a catastrophic way. Our calculations show that for parameters appropriate for the black hole in our Galactic centre, the plasmoid can attain relativistic speeds in about 35 min.     

Electromagnetic fields in jets

The magnetic fields and energy flows in an astronomical jet described by our earlier model are calculated in detail. Though the field distribution varies with the external pressure function   p ( z )  , it depends only weakly on the other boundary conditions. Individual field lines were plotted; the lines become nearly vertical at the bottom and are twisted at the top. An animation of a field line's motion was made, which shows the line being wound up by the accretion disc's differential rotation and rising as a result of this. The distribution of Poynting flux within the jet indicates that much of the energy flows up the jet from the inside of the accretion disc but a substantial fraction flows back down to the outside.     

Time-dependent simulations of steady C-type shocks

Using a time-dependent multifluid, magnetohydrodynamic code, we calculated the structure of steady perpendicular and oblique C-type shocks in dusty plasmas. We included relevant processes to describe mass transfer between the different fluids, radiative cooling by emission lines and grain charging, and studied the effect of single- and multiple-sized grains on the shock structure. Our models are the first of oblique fast-mode molecular shocks in which such a rigorous treatment of the dust grain dynamics has been combined with a self-consistent calculation of the thermal and ionization structures including appropriate microphysics. At low densities, the grains do not play any significant rôle in the shock dynamics. At high densities, the ionization fraction is sufficiently low that dust grains are important charge and current carriers and, thus, determine the shock structure. We find that the magnetic field in the shock front has a significant rotation out of the initial upstream plane. This is most pronounced for single-sized grains and small angles of the shock normal with the magnetic field. Our results are similar to previous studies of steady C-type shocks showing that our method is efficient, rigorous and robust. Unlike the method employed in the previous most detailed treatment of dust in steady oblique fast-mode shocks, ours allow a reliable calculation even when chemical or other conditions deviate from local statistical equilibrium. We are also able to model transient phenomena.     

Jet-driven disk accretion in low luminosity AGN?*

We explore an accretion model for low luminosity AGN (LLAGN) that attributes the low radiative output to a low mass accretion rate, , rather than a low radiative efficiency. In this model, electrons are assumed to drain energy from the ions as a result of collisionless plasma microinstabilities. Consequently, the accreting gas collapses to form a geometrically thin disk at small radii and is able to cool before reaching the black hole. The accretion disk is not a standard disk, however, because the radial disk structure is modified by a magnetic torque which drives a jet and which is primarily responsible for angular momentum transport. We also include relativistic effects. We apply this model to the well known LLAGN M87 and calculate the combined disk-jet steady-state broadband spectrum. A comparison between predicted and observed spectra indicates that M87 may be a maximally spinning black hole accreting at a rate of ∼10⁻³  M _⊙ yr⁻¹. This is about 6 orders of magnitude below the Eddington rate for the same radiative efficiency. Furthermore, the total jet power inferred by our model is in remarkably good agreement with the value independently deduced from observations of the M87 jet on kiloparsec scales. ^* This paper has previously been published in Astrophysics and Space Science, vol. 310:3–4.