Activation of the Blandford–Znajek mechanism in collapsing stars

The collapse of massive stars may result in the formation of accreting black holes in their interiors. The accreting stellar matter may advect substantial magnetic flux on to the black hole and promote the release of its rotational energy via magnetic stresses (the Blandford–Znajek mechanism). In this paper we explore whether this process can explain the stellar explosions and relativistic jets associated with long gamma-ray bursts. In particular, we show that the Blandford–Znajek mechanism is activated when the rest mass–energy density of matter drops below the energy density of the magnetic field in the near vicinity of the black hole (within its ergosphere). We also discuss whether such a strong magnetic field is in conflict with the rapid rotation of the stellar core required in the collapsar model, and suggest that the conflict can be avoided if the progenitor star is a component of a close binary. In this case the stellar rotation can be sustained via spin-orbital interaction. In an alternative scenario the magnetic field is generated in the accretion disc, but in this case the magnetic flux through the black hole ergosphere is not expected to be sufficiently high to explain the energetics of hypernovae by the BZ mechanism alone. However, this energy deficit can be recovered via the additional power provided by the disc.     

Models for jet power in elliptical galaxies: a case for rapidly spinning black holes

The power of jets from black holes is expected to depend on both the spin of the black hole and the structure of the accretion disc in the region of the last stable orbit. We investigate these dependencies using two different physical models for the jet power: the classical Blandford–Znajek (BZ) model and a hybrid model developed by Meier. In the BZ case, the jets are powered by magnetic fields directly threading the spinning black hole while in the hybrid model, the jet energy is extracted from both the accretion disc as well as the black hole via magnetic fields anchored to the accretion flow inside and outside the hole's ergosphere. The hybrid model takes advantage of the strengths of both the Blandford–Payne and BZ mechanisms, while avoiding the more controversial features of the latter. We develop these models more fully to account for general relativistic effects and to focus on advection-dominated accretion flows (ADAFs) for which the jet power is expected to be a significant fraction of the accreted rest mass energy.
We apply the models to elliptical galaxies, in order to see if these models can explain the observed correlation between the Bondi accretion rates and the total jet powers. For typical values of the disc viscosity parameter  α∼ 0.04 –0.3  and mass accretion rates consistent with ADAF model expectations, we find that the observed correlation requires   j ≳ 0.9  ; that is, it implies that the black holes are rapidly spinning. Our results suggest that the central black holes in the cores of clusters of galaxies must be rapidly rotating in order to drive jets powerful enough to heat the intracluster medium and quench cooling flows.     

Comptonization in the vicinity of the black hole horizon

Using a Monte Carlo method, we derive spectra arising from Comptonization taking place close to a Kerr black hole. We consider a model consisting of a hot thermal corona Comptonizing seed photons emitted by a cold accretion disc. We find that general relativistic effects are crucial for the emerging spectra in models, which involve significant contribution of radiation produced in the black hole ergosphere. As a result of this contribution, spectra of hard X-ray emission produced in the vicinity of a rapidly rotating black hole strongly depend on the inclination of the line of sight, with larger inclinations corresponding to harder spectra. Remarkably, such anisotropy could be responsible for properties of the X-ray spectra of Seyfert galaxies, which appear to be intrinsically harder in type 2 objects than in type 1, as reported recently.     

Non-evolutionarity of a reconnecting current sheet as a cause of its splitting into MHD shocks

Numerical simulations of the magnetic reconnection process in a current sheet show that, in some cases, MHD shocks appear to be attached to edges of the sheet. The appearance of the shocks may be considered to be a result of splitting of the sheet. In the present paper we suppose that this splitting takes place in consequence of non-evolutionarity of the reconnecting current sheet as a discontinuity. The problem of time evolution of small perturbations does not have a unique solution for a non-evolutionary discontinuity, and it splits into other (evolutionary) discontinuities. Such an approach allows us to determine conditions under which the splitting of the-sheet occurs. The main difficulty of this approach is that a current sheet is not reduced to a classified 1D discontinuity, because inhomogeneity of flow velocity inside the sheet is two-dimensional. To formulate the non-evolutionarity problem, we solve the linear MHD equations inside and outside the sheet and deduce linearized 1D boundary conditions at its surface. We show that for large enough conductivity, small perturbations exist which interact with the sheet as with a discontinuity. Then we obtain a non-evolutionarity criterion, with respect to these perturbations, in the form of a restriction on the flow velocity across the surface of the sheet.     

Magnetohydrodynamic properties of extragalactic jets

The theory that magnetic fields are instrumental in the formation and propagation of jets in active galactic nuclei dates back four decades. Despite a recent growing consensus on this notion stemming from the results of numerical simulations of magnetohydrodynamic (MHD) flows near black holes, the precise dynamical role of magnetic fields in observed parsec and kiloparsec jets remains uncertain. Some of the unanswered fundamental questions about extragalactic jets include the location where the flow becomes relativistic and where acceleration and collimation terminate, as well as the specifics of how the flow interacts with the ISM. Such observed properties as superluminal motions and wiggled structures based on numerical simulations to constitute the foundation of an MHD paradigm for extragalactic jets. We particularly focus our attention to the M87 jet, which is one of the best candidates to investigate relativistic outflows in extragalactic system.     

Simulations of small-scale explosive events on the Sun

Small-scale explosive events or microflares occur throughout the chromospheric network of the Sun. They are seen as sudden bursts of highly Doppler-shifted spectral lines of ions formed at temperatures in the range 2×104–5×105 K. They tend to occur near regions of cancelling photospheric magnetic fields and are thought to be directly associated with magnetic field reconnection. Recent observations have revealed that they have a bi-directional jet structure reminiscent of Petschek reconnection. In this paper compressible MHD simulations of the evolution of a current sheet to a steady Petschek, jet-like configuration are computed using the Versatile Advection Code. We obtain velocity profiles that can be compared with recent ultraviolet line-profile observations. By choosing initial conditions representative of magnetic loops in the solar corona and chromosphere, it is possible to explain the fact that jets flowing outward into the corona are more extended and appear before jets flowing towards the chromosphere. This model can reproduce the high Doppler-shifted components of the line profiles, but the brightening at low velocities, near the center of the bi-directional jet, cannot be explained by this simple MHD model.     

Accretion of a massive magnetized torus on a rotating black hole

We present numerical simulations of the axisymmetric accretion of a massive magnetized plasma torus on a rotating black hole. We use a realistic equation of state, which takes into account neutrino cooling and energy loss due to nucleus dissociations. The calculation are performed in the ideal relativistic MHD approximation using an upwind conservative scheme that is based on a linear Riemann solver and the constrained transport method to evolve the magnetic field. The gravitational attraction of the black hole is introduced via the Kerr metric in the Kerr–Schild coordinates. We simulate various magnetic field configurations and torus models, both optically thick and thin for neutrinos.We have found an effect of alternation of the magnetic field orientation in the ultrarelativistic jet formed as a result of the collapse. The calculations show evidence for heating of the wind surrounding the collapsar by the shock waves generated at the jet–wind border. It is shown that the neutrino cooling does not significantly change either the structure of the accretion flow or the total energy release of the system. The angular momentum of the accreting matter defines the time scale of the accretion. Due to the absence of the magnetic dynamo in our calculations, the initial strength and topology of the magnetic field determines the magnetization of the black hole, jet formation properties and the total energy yield. We estimate the total energy of accretion which transformed to jets as 1.3 × 10⁵² ergs which was sufficient to explain hypernova explosions like GRB 980425 or GRB 030329.     

Role of the pressure anisotropy in the relativistic pulsar wind

We generalize the hot relativistic MHD wind analysis to include the anisotropy of the pressure created in the pulsar wind by the strong magnetic field. Even with anisotropy the relativistic MHD equations integrate. In a very intense magnetic field, the motion of relativistic particles becomes rapidly one-dimensional in the direction of the field due to the very important radiative losses. Consequently, their distribution function becomes also one-dimensional and the component of the pressure, in the direction perpendicular to the magnetic field, decrease. In the limitP _⊥?0,P _∥≠0 we obtain a solution for the fluid flow which, starting at the neutron star surface, reaches smoothly infinity.     

Dynamics of relativistic jets

We discuss the structure and relativistic kinematics that develop in three spatial dimensions when a moderately hot, supersonic jet propagates into a denser background medium and encounters resistance from an oblique magnetic field. Our simulations incorporate relativistic MHD in a four-dimensional spacetime and clearly show that (a) relatively weak, oblique fields (at 1/16 of the equipartition value) have only a negligible influence on the propagating jet and they are passively pushed away by the relativistically moving head; (b) oblique fields in equipartition with the ambient plasma provide more resistance and cause bending at the jet head, but the magnitude of this deflection and the associated backflow are small compared to those identified by previous studies. The new results are understood as follows: Relativistic simulations have consistently shown that these jets are effectively heavy and so they do not suffer substantial momentum losses and are not decelerated as efficiently as their nonrelativistic counterparts. In addition, the ambient magnetic field, however strong, can be pushed aside with relative ease by the beam, provided that the degrees of freedom associated with all three spatial dimensions are followed self-consistently during the simulations. The effect is analogous to pushing Japanese “noren” or vertical Venetian blinds out of the way while the slats are allowed to bend and twist in 3-D space. Applied to relativistic extragalactic jets from blazars, the new results are encouraging since superluminal outflows exhibit bending near their sources and their environments are profoundly magnetized – but observations do not provide support for irregular kinematics such as large-scale vortical motions and pronounced reverse flows near the points of origin.     

黑洞吸积盘内边缘半径的演化与Blandford-Znajek过程

分别在薄盘和厚盘两种情况下,详细讨论了Ｂｌａｎｄｆｏｒｄ－Ｚｎａｊｅｋ过程对黑洞吸积盘内边缘半径演化的影响．结果表明,中心黑洞的无量纲角动量演化到某一临界值时;黑洞吸积盘的内边缘半径的变化率将由负交正．尽管盘的内边缘半径的极小值出现在Ｋｅｒｒ黑洞的能层的外边界上,但内边缘半径仍不可避免地会演化到Ｋｅｒｒ黑洞的能层中去．     

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.     

The jet power extracted from a magnetized accretion disc

We consider the power of a relativistic jet accelerated by the magnetic field of an accretion disc. It is found that the power extracted from the disc is mainly determined by the field strength and configuration of the field far from the disc. Comparing it with the power extracted from a rotating black hole, we find that the jet power extracted from a disc can dominate over that from the rotating black hole. However, in some cases, the jet power extracted from a rapidly rotating hole can be more important than that from the disc, even if the poloidal field threading the hole is not significantly larger than that threading the inner edge of the disc. The results imply that the radio-loudness of quasars may be governed by its accretion rate, which might be regulated by the central black hole mass. It is proposed that the different disc field generation mechanisms might be tested against observations of radio-loud quasars if their black hole masses are available.     

Steady and Time-Dependent MHD Modelling of Jets

A brief review is given of some results of our work on the construction of (I) steady and (II) time-dependent MHD models for nonrelativistic and relativistic astrophysical outflows and jets, analytically and numerically. The only available exact solutions for MHD outflows are those in separable coordinates, i.e., with the symmetry of radial or meridional self-similarity. Physically accepted solutions pass from the fast magnetosonic separatrix surface in order to satisfy MHD causality. An energetic criterion is outlined for selecting radially expanding winds from cylindrically expanding jets. Numerical simulations of magnetic self-collimation verify the conclusions of analytical steady solutions. We also propose a two-component model consisting of a wind outflow from a central object and a faster rotating outflow launched from a surrounding accretion disk which plays the role of the flow collimator. We also discuss the problem of shock formation during the magnetic collimation of wind-type outflows into jets.     

From Small-Scale Dynamo to Isotropic MHD Turbulence

We consider the problem of incompressible, forced, nonhelical, homogeneous, isotropic MHD turbulence with no mean magnetic field. This problem is essentially different from the case with externally imposed uniform mean field. There is no scale-by-scale equipartition between magnetic and kinetic energies as would be the case for the Alfvén-wave turbulence. The isotropic MHD turbulence is the end state of the turbulent dynamo which generates folded fields with small-scale direction reversals. We propose that the statistics seen in numerical simulations of isotropic MHD turbulence could be explained as a superposition of these folded fields and Alfvén-like waves that propagate along the folds.     

The onset of General Relativity: gravitationally redshifted emission lines

We study and quantify gravitational redshift by means of relativistic ray tracing simulations of emission lines. The emitter model is based on thin, Keplerian rotating rings in the equatorial plane of a rotating black hole. Emission lines are characterised by a generalized fully relativistic Doppler factor or redshift associated with the line core. Two modes of gravitational redshift, shift and distortion, become stronger with the emitting region closer to the Kerr black hole. Shifts of the line cores reveal an effect at levels of 0.0015 to 60% at gravitational radii ranging from 10⁵ to 2. The corresponding Doppler factors range from 0.999985 to 0.4048. Line shape distortion by strong gravity, i.e. very skewed and asymmetric lines occur at radii smaller than roughly ten gravitational radii. Gravitational redshift decreases with distance to the black hole but remains finite due to the asymptotical flatness of black hole space–time. The onset of gravitational redshift can be tested observationally with sufficient spectral resolution. Assuming a resolving power of ∼100000, yielding a resolution of ≈0.1 Å for optical and near‐infrared broad emission lines, the gravitational redshift can be probed out to approximately 75000 gravitational radii. In general, gravitational redshift is an indicator of black hole mass and spin as well as for the inclination angle of the emitter, e.g. an accretion disk. We suggest to do multi‐wavelength observations because all redshifted features should point towards the same central mass. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Direct numerical simulations of the Blandford–Znajek effect

The time-dependent general relativistic equations of degenerate electrodynamics are solved numerically in order to study the mechanism of the electromagnetic extraction of the rotational energy of black holes. We performed a series of 2D runs for black holes with specific angular momentum, a , from 0.1 to 0.9 and for a monopole magnetic field assuming axisymmetry. In the inner region of the wind, the solution quickly settles to a steady state with an outgoing Poynting flux. In all cases the angular velocity of the magnetic field lines is almost half the angular velocity of the black hole. Thus, at least for the configuration considered, the Blandford–Znajek mechanism operates near its maximum power output.     

A possible energy mechanism for cosmological γ-ray bursts

We suggest that an extreme Kerr black hole with a mass ∼10⁶ M_⊙, a dimensionless angular momentum     and a marginally stable orbital radius     located in a normal galaxy may produce a γ -ray burst (GRB) by capturing and disrupting a star. During the capture period, a transient accretion disc is formed and a strong transient magnetic field ∼     lasting for     may be produced at the inner boundary of the accretion disc. A large amount of rotational energy of the black hole is extracted and released in an ultrarelativistic jet with a bulk Lorentz factor Γ larger than 10³ via the Blandford–Znajek process. The relativistic jet energy can be converted into γ -radiation via an internal shock mechanism. The GRB duration should be the same as the lifetime of the strong transient magnetic field. The maximum number of sub-bursts is estimated to be     because the disc material is likely to break into pieces with a size about the thickness of the disc h at the cusp     The shortest risetime of the burst estimated from this model is ∼     The model GRB density rate is also estimated.     

Simulations of the axisymmetric magnetospheres of neutron stars

In this paper we present the results of time-dependent simulations of the dipolar axisymmetric magnetospheres of neutron stars carried out within the frameworks of both relativistic magnetohydrodynamics (MHD) and resistive force-free electrodynamics. The results of force-free simulations reveal the inability of our numerical method to accommodate the equatorial current sheets of pulsar magnetospheres, and raise a question mark about the robustness of this approach. On the other hand, the MHD approach allows us to make significant progress. We start with a non-rotating magnetically dominated dipolar magnetosphere and follow its evolution as the stellar rotation is switched on. We find that the time-dependent solution gradually approaches a steady state that is very close to the stationary solution of the pulsar equation found in 1999 by Contopoulos, Kazanas & Fendt. This result suggests that other stationary solutions that have the Y-point located well inside the light cylinder are unstable. The role of particle inertia and pressure on the structure and dynamics of MHD magnetospheres is studied in detail, as well as the potential implications of dissipative processes in the equatorial current sheet. We argue that pulsars may have differentially rotating magnetospheres which develop noticeable structural oscillations, and that this may help to explain the nature of the subpulse phenomena.     

Linkage between accretion disks and blazars

The magnetic field in an accretion disk is estimated assuming that all of the angular momentum within prescribed accretion disk radii is removed by a jet. The magnetic field estimated at the base of the jet is extrapolated to the blazar emission region using a model for a relativistic axisymmetric jet combined with some simplifying assumptions based on the relativistic nature of the flow. The extrapolated magnetic field is compared with estimates based upon the synchrotron and inverse Compton emission from three blazars, MKN 501, MKN 421 and PKS 2155-304. The magnetic fields evaluated from pure synchrotron self-Compton models are inconsistent with the magnetic fields extrapolated in this way. However, in two cases inverse Compton models in which a substantial part of the soft photon field is generated locally agree well, mainly because these models imply magnetic field strengths consistent with an important Poynting Flux component. This comparison is based on estimating the mass accretion rate from the jet energy flux. Further comparisons along these lines will be facilitated by independent estimates of the mass accretion rate in blazars and by more detailed models for jet propagation near the black hole.